Contribution of bistability and noise to cell fate transitions determined by feedback opening

AbstractAlternative cell fates represent a form of non-genetic diversity, which can promote adaptation and functional specialization. It is difficult to predict the rate of the transition between two cell fates due to the strong effect of noise on feedback loops and missing parameters. We opened synthetic positive feedback loops experimentally to obtain open-loop functions. These functions allowed us to identify a deterministic model of bistability by bypassing noise and the requirement to resolve individual processes in the loop. Combining the open-loop function with kinetic measurements and reintroducing the measured noise, we were able to predict the transition rates for the feedback systems without parameter tuning. Noise in gene expression was the key determinant of the transition rates inside the bistable range. Transitions between two cell fates were also observed outside of the bistable range, evidenced by bimodality and hysteresis. In this case, a slow transient process was the rate-limiting step in the transitions. Thus, feedback opening is an effective approach to identify the determinants of cell fate transitions and to predict their rates

Steady-state fluctuations of a genetic feedback loop with fluctuating rate parameters using the unified colored noise approximation

A common model of stochastic auto-regulatory gene expression describes promoter switching via cooperative protein binding, effective protein production in the active state and dilution of proteins. Here we consider an extension of this model whereby colored noise with a short correlation time is added to the reaction rate parameters -- we show that when the size and timescale of the noise is appropriately chosen it accounts for fast reactions that are not explicitly modelled, e.g., in models with no mRNA description, fluctuations in the protein production rate can account for rapid multiple stages of nuclear mRNA processing which precede translation in eukaryotes. We show how the unified colored noise approximation can be used to derive expressions for the protein number distribution that is in good agreement with stochastic simulations. We find that even when the noise in the rate parameters is small, the protein distributions predicted by our model can be significantly different than models assuming constant reaction rates.Comment: 33 pages, 10 figure

Synergistic dual positive feedback loops established by molecular sequestration generate robust bimodal response

Feedback loops are ubiquitous features of biological networks and can produce significant phenotypic heterogeneity, including a bimodal distribution of gene expression across an isogenic cell population. In this work, a combination of experiments and computational modeling was used to explore the roles of multiple feedback loops in the bimodal, switch-like response of the Saccharomyces cerevisiae galactose regulatory network. Here, we show that bistability underlies the observed bimodality, as opposed to stochastic effects, and that two unique positive feedback loops established by Gal1p and Gal3p, which both regulate network activity by molecular sequestration of Gal80p, induce this bimodality. Indeed, systematically scanning through different single and multiple feedback loop knockouts, we demonstrate that there is always a concentration regime that preserves the system’s bimodality, except for the double deletion of GAL1 and the GAL3 feedback loop, which exhibits a graded response for all conditions tested. The constitutive production rates of Gal1p and Gal3p operate as bifurcation parameters because variations in these rates can also abolish the system’s bimodal response. Our model indicates that this second loss of bistability ensues from the inactivation of the remaining feedback loop by the overexpressed regulatory component. More broadly, we show that the sequestration binding affinity is a critical parameter that can tune the range of conditions for bistability in a circuit with positive feedback established by molecular sequestration. In this system, two positive feedback loops can significantly enhance the region of bistability and the dynamic response time

Emergent Bistability : Effects of Additive and Multiplicative Noise

Positive feedback and cooperativity in the regulation of gene expression are generally considered to be necessary for obtaining bistable expression states. Recently, a novel mechanism of bistability termed emergent bistability has been proposed which involves only positive feedback and no cooperativity in the regulation. An additional positive feedback loop is effectively generated due to the inhibition of cellular growth by the synthesized proteins. The mechanism, demonstrated for a synthetic circuit, may be prevalent in natural systems also as some recent experimental results appear to suggest. In this paper, we study the effects of additive and multiplicative noise on the dynamics governing emergent bistability. The calculational scheme employed is based on the Langevin and Fokker-Planck formalisms. The steady state probability distributions of protein levels and the mean first passage times are computed for different noise strengths and system parameters. In the region of bistability, the bimodal probability distribution is shown to be a linear combination of a lognormal and a Gaussian distribution. The variances of the individual distributions and the relative weights of the distributions are further calculated for varying noise strengths and system parameters. The experimental relevance of the model results is also pointed out.Comment: 16 pages, 11 figures, version accepted for publication in Eur. Phys. J.

Characterization of bistability and transition rates in transcriptional positive feedback loops

Positive feedback commonly displays bistability, the ability to maintain overtime in the same conditions two alternative states of activity. The presence and the range of bistability depend on ultrasensitive reactions within the loop. To investigate bistability in genetic network, we constructed synthetic feedback loops in yeast where a transcription factor activates its own expression. By measuring the presence of hysteresis behavior, which is a sign of bistability, in those loops we identified the ultrasensitive reactions supporting bistability: homodimerization and cooperative binding of transcription factor. In the absence of those reactions the feedback loop was strictly monostable and when combined an even wider range of bistability arises than when there was only a single reaction. The detection of those reactions was made possible because we introduced RNA stem-loop upstream of the coding sequence of the transcription factor to reduce its translation rate. Indeed, the initial constructs had strong growth defect due to the overexpression of the transcription factor. Next, we aimed to predict transition rates between the two states of activity. Indeed, Even though the activity converges to either of the two states in the bistable range, due to the noise arising from the low number of some chemical species, transitions between the two states occur. The prediction of those transitions is difficult as the noise is amplified by feedback loop. First, we obtained a deterministic description of the loops by the open-loop approach. By breaking the loops at the mRNA of the transcription factor, we were able to fit the main parameter values of the system and map the steady states and the bistable range. Then, we determined the transient kinetics which is the activation delay which is not inherent to feedback loop, in our case it was the slow diffusion or binding of a ligand of the transcription factor. We determined also the noise of the system by measuring the distribution of mRNA at the steady states of the feedback loops. By building a stochastic model with the information from open-loop approach and expending it and fitting its parameter values to match the transient kinetics and noise observed, we were able to predict the transition rates observed in the feedback loops. With this better understanding, we discovered that the transitions are led by either noise or slow transient kinetics depending whether the system is inside or outside in the vicinity of the bistable range, respectively. Finally, we showed that the transition rates were abruptly changing around the boundaries of the bistable region. Therefore, the bistable region can be estimated in similar feedback loops by simply measuring transition rates

Origins of Binary Gene Expression in Post-transcriptional Regulation by MicroRNAs

MicroRNA-mediated regulation of gene expression is characterised by some distinctive features that set it apart from unregulated and transcription factor-regulated gene expression. Recently, a mathematical model has been proposed to describe the dynamics of post-transcriptional regulation by microRNAs. The model explains the observations made in single cell experiments quite well. In this paper, we introduce some additional features into the model and consider two specific cases. In the first case, a non-cooperative positive feedback loop is included in the transcriptional regulation of the target gene expression. In the second case, a stochastic version of the original model is considered in which there are random transitions between the inactive and active expression states of the gene. In the first case we show that bistability is possible in a parameter regime, due to the presence of a non-linear protein decay term in the gene expression dynamics. In the second case, we derive the conditions for obtaining stochastic binary gene expression. We find that this type of gene expression is more favourable in the case of regulation by microRNAs as compared to the case of unregulated gene expression. The theoretical predictions relating to binary gene expression are experimentally testable.Comment: 10 Pages, 5 Figure

Interplay of Extrinsic and Intrinsic Cues in Cell-Fate Decisions

A cell’s decision making process is coordinated by dynamic interplay between its extracellular environment and its intracellular milieu. For example, during stem cell differentiation, fate decisions are believed to be ultimately controlled by differential expression of lineage-specific transcription factors, but cytokine receptor signals also play a crucial instructive role in addition to providing permissive proliferation and survival cues. Here, we present a minimal computational framework that integrates the intrinsic and extrinsic regulatory elements implicated in the commitment of hematopoietic progenitor cells to mature red blood cells (Chapter 2). Our model highlights the importance of bidirectional interactions between cytokine receptors and transcription factors in conferring properties such as ultrasensitivity and bistability to differentiating cells. These system-level properties can induce a switch-like characteristic during differentiation and provide robustness to the mature state. We then experimentally test predictions from this lineage commitment model in a model system for studying erythropoiesis (Chapter 3). Our experiments show that hemoglobin synthesis is highly switch-like in response to cytokine and cells undergoing lineage commitment possess memory of earlier cytokine signals. We show that erythrocyte-specific receptor and transcription factor are indeed synchronously co-upregulated and the heterogeneity in their expression is positively correlated during differentiation, confirming the presence of autofeedback and receptor-mediated positive feedback loops. To evaluate the possibility of employing this minimal topology as a synthetic “memory module” for cell engineering applications, we constructed this topology synthetically in Saccharomyces cerevisiae by integrating Arabidopsis thaliana signaling components with an endogenous yeast pathway (Chapter 4). Our experiments show that any graded and unimodal signaling pathway can be rationally rewired to achieve our desired topology and the resulting network immediately attains high ultrasensitivity and bimodality without tweaking. We further show that this topology can be tuned to regulate system dynamics such as activation/deactivation kinetics, signal amplitude, switching threshold and sensitivity. We conclude with a computational study to explore the generality of this interplay between extrinsic and intrinsic cues in hematopoiesis. We extend our minimal model analysis in Chapter 2 to examine the more complex fate decisions in bipotent and multipotent progenitors, particularly how these cells can make robust decisions in the presence of multiple extrinsic cues and intrinsic noise (Chapter 5). Our model provides support to both the instructive and stochastic theories of commitment: cell fates are ultimately driven by lineage-specific transcription factors, but cytokine signaling can strongly bias lineage commitment by regulating these inherently noisy cell-fate decisions with complex, pertinent behaviors such as ligand-mediated ultrasensitivity and robust multistability. The simulations further suggest that the kinetics of differentiation to a mature cell state can depend on the starting progenitor state as well as on the route of commitment that is chosen. Lastly, our model shows good agreement with lineage-specific receptor expression kinetics from microarray experiments and provides a computational framework that can integrate both classical and alternative commitment paths in hematopoiesis that have been observed experimentally

Bistability, Probability Transition Rate and First-Passage Time in an Autoactivating Positive-Feedback Loop

A hallmark of positive-feedback regulation is bistability, which gives rise to distinct cellular states with high and low expression levels, and that stochasticity in gene expression can cause random transitions between two states, yielding bimodal population distribution (Kaern et al., 2005, Nat Rev Genet 6: 451-464). In this paper, the probability transition rate and first-passage time in an autoactivating positive-feedback loop with bistability are investigated, where the gene expression is assumed to be disturbed by both additive and multiplicative external noises, the bimodality in the stochastic gene expression is due to the bistability, and the bistability determines that the potential of the Fokker-Planck equation has two potential wells. Our main goal is to illustrate how the probability transition rate and first-passage time are affected by the maximum transcriptional rate, the intensities of additive and multiplicative noises, and the correlation of additive and multiplicative noises. Our main results show that (i) the increase of the maximum transcription rate will be useful for maintaining a high gene expression level; (ii) the probability transition rate from one potential well to the other one will increase with the increase of the intensity of additive noise; (iii) the increase of multiplicative noise strength will increase the amount of probability in the left potential well; and (iv) positive (or negative) cross-correlation between additive and multiplicative noises will increase the amount of probability in the left (or right) potential well

Stochastic modeling of auto-regulatory genetic feedback loops

Auto-regulatory feedback loops are one of the most common network motifs. A wide variety of stochastic models have been constructed to understand how the fluctuations in protein numbers in these loops are influenced by the kinetic parameters of the main biochemical steps. These models differ according to (i) which sub-cellular processes are explicitly modelled; (ii) the modelling methodology employed (discrete, continuous or hybrid); (iii) whether they can be analytically solved for the steady-state distribution of protein numbers. We discuss the assumptions and properties of the main models in the literature, summarize our current understanding of the relationship between them and highlight some of the insights gained through modelling.Comment: 12 pages, 3 figures. Submitted to Biophysical Journa