Noise is often indispensable to key cellular activities, such as gene expression, necessitating t... more Noise is often indispensable to key cellular activities, such as gene expression, necessitating the use of stochastic models to capture its dynamics. The chemical master equation (CME) is a commonly used stochastic model of Kolmogorov forward equations that describe how the probability distribution of a chemically reacting system varies with time. Finding analytic solutions to the CME can have benefits, such as expediting simulations of multiscale biochemical reaction networks and aiding the design of distributional responses. However, analytic solutions are rarely known. A recent method of computing analytic stationary solutions relies on gluing simple state spaces together recursively at one or two states. We explore the capabilities of this method and introduce algorithms to derive analytic stationary solutions to the CME. We first formally characterize state spaces that can be constructed by performing single-state gluing of paths, cycles or both sequentially. We then study stoc...
Noise is indispensible to key cellular activities, including gene expression coordination and pro... more Noise is indispensible to key cellular activities, including gene expression coordination and probabilistic differentiation. Stochastic models, such as the chemical master equation (CME), are essential to model noise in the levels of cellular components. In the CME framework, each state is associated with the molecular counts of all component species, and specifies the probability for the system to have that set of molecular counts. Analytic solutions to the CME are rarely known but can bring exciting benefits. For instance, simulations of biochemical reaction networks that are multiscale in time can be sped up tremendously by incorporating analytic solutions of the slow time-scale dynamics. Analytic solutions also enable the design of stationary distributions with properties such as the modality of the distribution, the mean expression level, and the level of noise. One way to derive the analytic steady state response of a biochemical reaction network was recently proposed by (Mély...
This thesis advances our understanding of three important aspects of biological systems engineeri... more This thesis advances our understanding of three important aspects of biological systems engineering: analysis, design, and computational methods. First, biological circuit design is necessary to engineer biological systems that behave consistently and follow our design specifications. We contribute by formulating and solving novel problems in stochastic biological circuit design. Second, computational methods for solving biological systems are often limited by the nonlinearity and high dimensionality of the system’s dynamics. This problem is particularly extreme for the parameter identification of stochastic, nonlinear systems. Thus, we develop a method for parameter identification that relies on data-driven stochastic model reduction. Finally, biological system analysis encompasses understanding the stability, performance, and robustness of these systems, which is critical for their implementation. We analyze a sequestration feedback motif for implementing biological control. First...
Many biological and physical systems exhibit population-density dependent transitions to synchron... more Many biological and physical systems exhibit population-density dependent transitions to synchronized oscillations in a process often termed "dynamical quorum sensing". Synchronization frequently arises through chemical communication via signaling molecules distributed through an external media. We study a simple theoretical model for dynamical quorum sensing: a heterogenous population of limit-cycle oscillators diffusively coupled through a common media. We show that this model exhibits a rich phase diagram with four qualitatively distinct mechanisms fueling population-dependent transitions to global oscillations, including a new type of transition we term "dynamic death". We derive a single pair of analytic equations that allows us to calculate all phase boundaries as a function of population density and show that the model reproduces many of the qualitative features of recent experiments of BZ catalytic particles as well as synthetically engineered bacteria.
An ongoing area of study in synthetic biology has been the design and construction of synthetic c... more An ongoing area of study in synthetic biology has been the design and construction of synthetic circuits that maintain homeostasis at the population level. Here, we are interested in designing a synthetic control circuit that regulates the total cell population and the relative ratio between cell strains in a culture containing two different cell strains. We have developed a dual feedback control strategy that uses two separate control loops to achieve the two functions respectively. By combining both of these control loops, we have created a population regulation circuit where both the total population size and relative cell type ratio can be set by reference signals. The dynamics of the regulation circuit show robustness and adaptation to perturbations in cell growth rate and changes in cell numbers. The control architecture is general and could apply to any organism for which synthetic biology tools for quorum sensing, comparison between outputs, and growth control are available.
SummaryIntegral control is commonly used in mechanical and electrical systems to ensure perfect a... more SummaryIntegral control is commonly used in mechanical and electrical systems to ensure perfect adaptation. A proposed design of integral control for synthetic biological systems employs the sequestration of two biochemical controller species. The unbound amount of controller species captures the integral of the error between the current and the desired state of the system. However, implementing integral control inside bacterial cells using sequestration feedback has been challenging due to the controller molecules being degraded and diluted. Furthermore, integral control can only be achieved under stability conditions that not all sequestration feedback networks fulfill. In this work, we give guidelines for ensuring stability and good performance (small steady-state error) in sequestration feedback networks. Our guidelines provide simple tuning options to obtain a flexible and practical biological implementation of sequestration feedback control. Using tools and metrics from contro...
It has been an ongoing scientific debate whether biological parameters are conserved across exper... more It has been an ongoing scientific debate whether biological parameters are conserved across experimental setups with different media, pH values, and other experimental conditions. Our work explores this question using Bayesian probability as a rigorous framework to assess the biological context of parameters in a model of the cell growth controller in You et al. When this growth controller is uninduced, the E. coli cell population grows to carrying capacity; however, when the circuit is induced, the cell population growth is regulated to remain well below carrying capacity. This growth control controller regulates the E. coli cell population by cell–cell communication using the signaling molecule AHL and by cell death using the bacterial toxin CcdB.To evaluate the context dependence of parameters such as the cell growth rate, the carrying capacity, the AHL degradation rate, the leakiness of AHL, the leakiness of toxin CcdB, and the IPTG induction factor, we collect experimental data...
SummaryFeedback regulation is pervasive in biology at both the organismal and cellular level. In ... more SummaryFeedback regulation is pervasive in biology at both the organismal and cellular level. In this article, we explore the properties of a particular biomolecular feedback mechanism implemented using the sequestration binding of two molecules. Our work develops an analytic framework for understanding the hard limits, performance tradeoffs, and architectural properties of this simple model of biological feedback control. Using tools from control theory, we show that there are simple parametric relationships that determine both the stability and the performance of these systems in terms of speed, robustness, steady-state error, and leakiness. These findings yield a holistic understanding of the behavior of sequestration feedback and contribute to a more general theory of biological control systems.
The finite state projection (FSP) approach to solving the chemical master equation (CME) has enab... more The finite state projection (FSP) approach to solving the chemical master equation (CME) has enabled successful inference of discrete stochastic models to predict single-cell gene regulation dynamics. Unfortunately, the FSP approach is highly computationally intensive for all but the simplest models, an issue that is highly problematic when parameter inference and uncertainty quantification takes enormous numbers of parameter evaluations. To address this issue, we propose two new computational methods for the Bayesian inference of stochastic gene expression parameters given single-cell experiments. First, we present an adaptive scheme to improve parameter proposals for Metropolis-Hastings sampling using full FSP-based likelihood evaluations. We then formulate and verify an Adaptive Delayed Acceptance Metropolis-Hastings (ADAMH) algorithm to utilize with reduced Krylov-basis projections of the FSP. We test and compare both algorithms on three example models and simulated data to show...
Stochastic gene expression poses an important challenge for engineering robust behaviors in a het... more Stochastic gene expression poses an important challenge for engineering robust behaviors in a heterogeneous cell population. Cells address this challenge by operating on distributions of cellular responses generated by noisy processes. Similarly, a previously published temporal logic gate considers the distribution of responses across a cell population under chemical inducer pulsing events. The design uses a system of two integrases to engineer an E. coli strain with four DNA states that records the temporal order of two chemical signal events. The heterogeneous cell population response was used to infer the timing and duration of the two chemical signals for a small set of events. Here we use the temporal logic gate system to address the problem of extracting information about chemical signal events. We use the heterogeneous cell population response to infer whether any event has occurred or not and also to infer its properties such as timing and amplitude. Bayesian inference provi...
Noise is often indispensable to key cellular activities, such as gene expression, necessitating t... more Noise is often indispensable to key cellular activities, such as gene expression, necessitating the use of stochastic models to capture its dynamics. The chemical master equation (CME) is a commonly used stochastic model that describes how the probability distribution of a chemically reacting system varies with time. Knowing analytic solutions to the CME can have benefits, such as expediting simulations of multiscale biochemical reaction networks and aiding the design of distributional responses. However, analytic solutions are rarely known. A recent method of computing analytic stationary solutions relies on gluing simple state spaces together recursively at one or two states. We explore the capabilities of this method and introduce algorithms to derive analytic stationary solutions to the CME. We first formally characterise state spaces that can be constructed by performing single-state gluing of paths, cycles or both sequentially. We then study stochastic biochemical reaction net...
Journal of the Royal Society, Interface, May 1, 2017
Noise is often indispensable to key cellular activities, such as gene expression, necessitating t... more Noise is often indispensable to key cellular activities, such as gene expression, necessitating the use of stochastic models to capture its dynamics. The chemical master equation (CME) is a commonly used stochastic model of Kolmogorov forward equations that describe how the probability distribution of a chemically reacting system varies with time. Finding analytic solutions to the CME can have benefits, such as expediting simulations of multiscale biochemical reaction networks and aiding the design of distributional responses. However, analytic solutions are rarely known. A recent method of computing analytic stationary solutions relies on gluing simple state spaces together recursively at one or two states. We explore the capabilities of this method and introduce algorithms to derive analytic stationary solutions to the CME. We first formally characterize state spaces that can be constructed by performing single-state gluing of paths, cycles or both sequentially. We then study stoc...
An ongoing area of study in synthetic biology has been the design and construction of synthetic c... more An ongoing area of study in synthetic biology has been the design and construction of synthetic circuits that maintain homeostasis at the population level. Here, we are interested in designing a synthetic control circuit that regulates the total cell population and the relative ratio between cell strains in a culture containing two different cell strains. We have developed a dual feedback control strategy that uses two separate control loops to achieve the two functions respectively. By combining both of these control loops, we have created a population regulation circuit where both the total population size and relative cell type ratio can be set by reference signals. The dynamics of the regulation circuit show robustness and adaptation to perturbations in cell growth rate and changes in cell numbers. The control architecture is general and could apply to any organism for which synthetic biology tools for quorum sensing, comparison between outputs, and growth control are available.
Many biological and physical systems exhibit population-density dependent transitions to synchron... more Many biological and physical systems exhibit population-density dependent transitions to synchronized oscillations in a process often termed "dynamical quorum sensing". Synchronization frequently arises through chemical communication via signaling molecules distributed through an external medium. We study a simple theoretical model for dynamical quorum sensing: a heterogenous population of limit-cycle oscillators diffusively coupled through a common medium. We show that this model exhibits a rich phase diagram with four qualitatively distinct physical mechanisms that can lead to a loss of coherent population-level oscillations, including a novel mechanism arising from effective time-delays introduced by the external medium. We derive a single pair of analytic equations that allow us to calculate phase boundaries as a function of population density and show that the model reproduces many of the qualitative features of recent experiments on BZ catalytic particles as well as synthetically engineered bacteria.
Noise is often indispensable to key cellular activities, such as gene expression, necessitating t... more Noise is often indispensable to key cellular activities, such as gene expression, necessitating the use of stochastic models to capture its dynamics. The chemical master equation (CME) is a commonly used stochastic model of Kolmogorov forward equations that describe how the probability distribution of a chemically reacting system varies with time. Finding analytic solutions to the CME can have benefits, such as expediting simulations of multiscale biochemical reaction networks and aiding the design of distributional responses. However, analytic solutions are rarely known. A recent method of computing analytic stationary solutions relies on gluing simple state spaces together recursively at one or two states. We explore the capabilities of this method and introduce algorithms to derive analytic stationary solutions to the CME. We first formally characterize state spaces that can be constructed by performing single-state gluing of paths, cycles or both sequentially. We then study stoc...
Noise is indispensible to key cellular activities, including gene expression coordination and pro... more Noise is indispensible to key cellular activities, including gene expression coordination and probabilistic differentiation. Stochastic models, such as the chemical master equation (CME), are essential to model noise in the levels of cellular components. In the CME framework, each state is associated with the molecular counts of all component species, and specifies the probability for the system to have that set of molecular counts. Analytic solutions to the CME are rarely known but can bring exciting benefits. For instance, simulations of biochemical reaction networks that are multiscale in time can be sped up tremendously by incorporating analytic solutions of the slow time-scale dynamics. Analytic solutions also enable the design of stationary distributions with properties such as the modality of the distribution, the mean expression level, and the level of noise. One way to derive the analytic steady state response of a biochemical reaction network was recently proposed by (Mély...
This thesis advances our understanding of three important aspects of biological systems engineeri... more This thesis advances our understanding of three important aspects of biological systems engineering: analysis, design, and computational methods. First, biological circuit design is necessary to engineer biological systems that behave consistently and follow our design specifications. We contribute by formulating and solving novel problems in stochastic biological circuit design. Second, computational methods for solving biological systems are often limited by the nonlinearity and high dimensionality of the system’s dynamics. This problem is particularly extreme for the parameter identification of stochastic, nonlinear systems. Thus, we develop a method for parameter identification that relies on data-driven stochastic model reduction. Finally, biological system analysis encompasses understanding the stability, performance, and robustness of these systems, which is critical for their implementation. We analyze a sequestration feedback motif for implementing biological control. First...
Many biological and physical systems exhibit population-density dependent transitions to synchron... more Many biological and physical systems exhibit population-density dependent transitions to synchronized oscillations in a process often termed "dynamical quorum sensing". Synchronization frequently arises through chemical communication via signaling molecules distributed through an external media. We study a simple theoretical model for dynamical quorum sensing: a heterogenous population of limit-cycle oscillators diffusively coupled through a common media. We show that this model exhibits a rich phase diagram with four qualitatively distinct mechanisms fueling population-dependent transitions to global oscillations, including a new type of transition we term "dynamic death". We derive a single pair of analytic equations that allows us to calculate all phase boundaries as a function of population density and show that the model reproduces many of the qualitative features of recent experiments of BZ catalytic particles as well as synthetically engineered bacteria.
An ongoing area of study in synthetic biology has been the design and construction of synthetic c... more An ongoing area of study in synthetic biology has been the design and construction of synthetic circuits that maintain homeostasis at the population level. Here, we are interested in designing a synthetic control circuit that regulates the total cell population and the relative ratio between cell strains in a culture containing two different cell strains. We have developed a dual feedback control strategy that uses two separate control loops to achieve the two functions respectively. By combining both of these control loops, we have created a population regulation circuit where both the total population size and relative cell type ratio can be set by reference signals. The dynamics of the regulation circuit show robustness and adaptation to perturbations in cell growth rate and changes in cell numbers. The control architecture is general and could apply to any organism for which synthetic biology tools for quorum sensing, comparison between outputs, and growth control are available.
SummaryIntegral control is commonly used in mechanical and electrical systems to ensure perfect a... more SummaryIntegral control is commonly used in mechanical and electrical systems to ensure perfect adaptation. A proposed design of integral control for synthetic biological systems employs the sequestration of two biochemical controller species. The unbound amount of controller species captures the integral of the error between the current and the desired state of the system. However, implementing integral control inside bacterial cells using sequestration feedback has been challenging due to the controller molecules being degraded and diluted. Furthermore, integral control can only be achieved under stability conditions that not all sequestration feedback networks fulfill. In this work, we give guidelines for ensuring stability and good performance (small steady-state error) in sequestration feedback networks. Our guidelines provide simple tuning options to obtain a flexible and practical biological implementation of sequestration feedback control. Using tools and metrics from contro...
It has been an ongoing scientific debate whether biological parameters are conserved across exper... more It has been an ongoing scientific debate whether biological parameters are conserved across experimental setups with different media, pH values, and other experimental conditions. Our work explores this question using Bayesian probability as a rigorous framework to assess the biological context of parameters in a model of the cell growth controller in You et al. When this growth controller is uninduced, the E. coli cell population grows to carrying capacity; however, when the circuit is induced, the cell population growth is regulated to remain well below carrying capacity. This growth control controller regulates the E. coli cell population by cell–cell communication using the signaling molecule AHL and by cell death using the bacterial toxin CcdB.To evaluate the context dependence of parameters such as the cell growth rate, the carrying capacity, the AHL degradation rate, the leakiness of AHL, the leakiness of toxin CcdB, and the IPTG induction factor, we collect experimental data...
SummaryFeedback regulation is pervasive in biology at both the organismal and cellular level. In ... more SummaryFeedback regulation is pervasive in biology at both the organismal and cellular level. In this article, we explore the properties of a particular biomolecular feedback mechanism implemented using the sequestration binding of two molecules. Our work develops an analytic framework for understanding the hard limits, performance tradeoffs, and architectural properties of this simple model of biological feedback control. Using tools from control theory, we show that there are simple parametric relationships that determine both the stability and the performance of these systems in terms of speed, robustness, steady-state error, and leakiness. These findings yield a holistic understanding of the behavior of sequestration feedback and contribute to a more general theory of biological control systems.
The finite state projection (FSP) approach to solving the chemical master equation (CME) has enab... more The finite state projection (FSP) approach to solving the chemical master equation (CME) has enabled successful inference of discrete stochastic models to predict single-cell gene regulation dynamics. Unfortunately, the FSP approach is highly computationally intensive for all but the simplest models, an issue that is highly problematic when parameter inference and uncertainty quantification takes enormous numbers of parameter evaluations. To address this issue, we propose two new computational methods for the Bayesian inference of stochastic gene expression parameters given single-cell experiments. First, we present an adaptive scheme to improve parameter proposals for Metropolis-Hastings sampling using full FSP-based likelihood evaluations. We then formulate and verify an Adaptive Delayed Acceptance Metropolis-Hastings (ADAMH) algorithm to utilize with reduced Krylov-basis projections of the FSP. We test and compare both algorithms on three example models and simulated data to show...
Stochastic gene expression poses an important challenge for engineering robust behaviors in a het... more Stochastic gene expression poses an important challenge for engineering robust behaviors in a heterogeneous cell population. Cells address this challenge by operating on distributions of cellular responses generated by noisy processes. Similarly, a previously published temporal logic gate considers the distribution of responses across a cell population under chemical inducer pulsing events. The design uses a system of two integrases to engineer an E. coli strain with four DNA states that records the temporal order of two chemical signal events. The heterogeneous cell population response was used to infer the timing and duration of the two chemical signals for a small set of events. Here we use the temporal logic gate system to address the problem of extracting information about chemical signal events. We use the heterogeneous cell population response to infer whether any event has occurred or not and also to infer its properties such as timing and amplitude. Bayesian inference provi...
Noise is often indispensable to key cellular activities, such as gene expression, necessitating t... more Noise is often indispensable to key cellular activities, such as gene expression, necessitating the use of stochastic models to capture its dynamics. The chemical master equation (CME) is a commonly used stochastic model that describes how the probability distribution of a chemically reacting system varies with time. Knowing analytic solutions to the CME can have benefits, such as expediting simulations of multiscale biochemical reaction networks and aiding the design of distributional responses. However, analytic solutions are rarely known. A recent method of computing analytic stationary solutions relies on gluing simple state spaces together recursively at one or two states. We explore the capabilities of this method and introduce algorithms to derive analytic stationary solutions to the CME. We first formally characterise state spaces that can be constructed by performing single-state gluing of paths, cycles or both sequentially. We then study stochastic biochemical reaction net...
Journal of the Royal Society, Interface, May 1, 2017
Noise is often indispensable to key cellular activities, such as gene expression, necessitating t... more Noise is often indispensable to key cellular activities, such as gene expression, necessitating the use of stochastic models to capture its dynamics. The chemical master equation (CME) is a commonly used stochastic model of Kolmogorov forward equations that describe how the probability distribution of a chemically reacting system varies with time. Finding analytic solutions to the CME can have benefits, such as expediting simulations of multiscale biochemical reaction networks and aiding the design of distributional responses. However, analytic solutions are rarely known. A recent method of computing analytic stationary solutions relies on gluing simple state spaces together recursively at one or two states. We explore the capabilities of this method and introduce algorithms to derive analytic stationary solutions to the CME. We first formally characterize state spaces that can be constructed by performing single-state gluing of paths, cycles or both sequentially. We then study stoc...
An ongoing area of study in synthetic biology has been the design and construction of synthetic c... more An ongoing area of study in synthetic biology has been the design and construction of synthetic circuits that maintain homeostasis at the population level. Here, we are interested in designing a synthetic control circuit that regulates the total cell population and the relative ratio between cell strains in a culture containing two different cell strains. We have developed a dual feedback control strategy that uses two separate control loops to achieve the two functions respectively. By combining both of these control loops, we have created a population regulation circuit where both the total population size and relative cell type ratio can be set by reference signals. The dynamics of the regulation circuit show robustness and adaptation to perturbations in cell growth rate and changes in cell numbers. The control architecture is general and could apply to any organism for which synthetic biology tools for quorum sensing, comparison between outputs, and growth control are available.
Many biological and physical systems exhibit population-density dependent transitions to synchron... more Many biological and physical systems exhibit population-density dependent transitions to synchronized oscillations in a process often termed "dynamical quorum sensing". Synchronization frequently arises through chemical communication via signaling molecules distributed through an external medium. We study a simple theoretical model for dynamical quorum sensing: a heterogenous population of limit-cycle oscillators diffusively coupled through a common medium. We show that this model exhibits a rich phase diagram with four qualitatively distinct physical mechanisms that can lead to a loss of coherent population-level oscillations, including a novel mechanism arising from effective time-delays introduced by the external medium. We derive a single pair of analytic equations that allow us to calculate phase boundaries as a function of population density and show that the model reproduces many of the qualitative features of recent experiments on BZ catalytic particles as well as synthetically engineered bacteria.
Uploads
Papers by Ania Baetica