Prediction and measurement of an autoregulatory genetic module

A high-throughput synthetic biology approach for studying combinatorial chromatin-based transcriptional regulation

The construction of synthetic gene circuits requires the rational combination of multiple regulatory components, but predicting their behavior can be challenging due to poorly understood component interactions and unexpected emergent behaviors. In eukaryotes, chromatin regulators (CRs) are essential regulatory components that orchestrate gene expression. Here, we develop a screening platform to investigate the impact of CR pairs on transcriptional activity in yeast. We construct a combinatorial library consisting of over 1,900 CR pairs and use a high-throughput workflow to characterize the impact of CR co-recruitment on gene expression. We recapitulate known interactions and discover several instances of CR pairs with emergent behaviors. We also demonstrate that supervised machine learning models trained with low-dimensional amino acid embeddings accurately predict the impact of CR co-recruitment on transcriptional activity. This work introduces a scalable platform and machine learning approach that can be used to study how networks of regulatory components impact gene expression.

Quantitative Dynamic Analysis of de novo Sphingolipid Biosynthesis in Arabidopsis thaliana

Sphingolipids are pivotal for plant development and stress responses. Growing interest has been directed towards fully comprehending the regulatory mechanisms of the sphingolipid pathway. We explore its de novo biosynthesis and homeostasis in Arabidopsis thaliana cell cultures, shedding light on fundamental metabolic mechanisms. Employing 15N isotope labeling and quantitative dynamic modeling approach, we developed a regularized and constraint-based Dynamic Metabolic Flux Analysis (r-DMFA) framework to predict metabolic shifts due to enzymatic changes. Our analysis revealed key enzymes such as sphingoid-base hydroxylase (SBH) and long-chain-base kinase (LCBK) to be critical for maintaining sphingolipid homeostasis. Disruptions in these enzymes were found to affect cellular viability and increase the potential for programmed cell death (PCD). Thus, this work enhances our understanding of sphingolipid metabolism and demonstrates the utility of dynamic modeling in analyzing complex metabolic pathways.

Analysis and design of single-cell experiments to harvest fluctuation information while rejecting measurement noise

Introduction: Despite continued technological improvements, measurement errors always reduce or distort the information that any real experiment can provide to quantify cellular dynamics. This problem is particularly serious for cell signaling studies to quantify heterogeneity in single-cell gene regulation, where important RNA and protein copy numbers are themselves subject to the inherently random fluctuations of biochemical reactions. Until now, it has not been clear how measurement noise should be managed in addition to other experiment design variables (e.g., sampling size, measurement times, or perturbation levels) to ensure that collected data will provide useful insights on signaling or gene expression mechanisms of interest. Methods: We propose a computational framework that takes explicit consideration of measurement errors to analyze single-cell observations, and we derive Fisher Information Matrix (FIM)-based criteria to quantify the information value of distorted experiments. Results and Discussion: We apply this framework to analyze multiple models in the context of simulated and experimental single-cell data for a reporter gene controlled by an HIV promoter. We show that the proposed approach quantitatively predicts how different types of measurement distortions affect the accuracy and precision of model identification, and we demonstrate that the effects of these distortions can be mitigated through explicit consideration during model inference. We conclude that this reformulation of the FIM could be used effectively to design single-cell experiments to optimally harvest fluctuation information while mitigating the effects of image distortion.

Enhanced Transcriptional Strength of HIV-1 Subtype C Minimizes Gene Expression Noise and Confers Stability to the Viral Latent State

Stochastic fluctuations in gene expression emanating from the HIV-1 long terminal repeat (LTR), amplified by the Tat positive feedback circuit, determine the choice between viral infection fates: active transcription (ON) or transcriptional silence (OFF). The emergence of several transcription factor binding site (TFBS) variant strains in HIV-1 subtype C (HIV-1C), especially those containing the duplication of the NF-κB motif, mandates the evaluation of the effect of enhanced transcriptional strength on gene expression noise and its influence on viral fate selection switch. Using a panel of subgenomic LTR-variant strains containing different copy numbers of the NF-κB motif (ranging from 0 to 4), we used flow cytometry, mRNA quantification, and pharmacological perturbations to demonstrate an inverse correlation between promoter strength and gene expression noise in Jurkat T cells and primary CD4⁺ T cells. The inverse correlation is consistent in clonal cell populations at constant intracellular concentrations of Tat and when NF-κB levels were regulated pharmacologically. Further, we show that strong LTRs containing at least two copies of the NF-κB motif in the enhancer establish a more stable latent state and demonstrate more rapid latency reversal than weak LTRs containing fewer motifs. We also demonstrate a cooperative binding of NF-κB to the motif cluster in HIV-1C LTRs containing two, three, or four NF-κB motifs (Hill coefficient [H] = 2.61, 3.56, and 3.75, respectively). The present work alludes to a possible evolution of the HIV-1C LTR toward gaining transcriptional strength associated with attenuated gene expression noise with implications for viral latency. IMPORTANCE Over the past two consecutive decades, HIV-1 subtype C (HIV-1C) has been undergoing directional evolution toward augmenting the transcriptional strength of the long terminal repeat (LTR) by adding more copies of the existing transcription factor binding site (TFBS) by sequence duplication. Additionally, the duplicated elements are genetically diverse, suggesting broader-range signal receptivity by variant LTRs. The HIV-1 promoter is inherently noisy, and the stochastic fluctuations in gene expression of variant LTRs may influence the active transcription (ON)/transcriptional silence (OFF) latency decisions. The evolving NF-κB motif variations of HIV-1C offer a powerful opportunity to examine how the transcriptional strength of the LTR might influence gene expression noise. Our work here shows that the augmented transcriptional strength of the HIV-1C LTR leads to concomitantly reduced gene expression noise, consequently leading to stabler latency maintenance and rapid latency reversal. The present work offers a novel lead toward appreciating the molecular mechanisms governing HIV-1 latency.

Negative autoregulation controls size scaling in confined gene expression reactions

Gene expression via transcription-translation is the most fundamental reaction to sustain biological systems, and complex reactions occur in a small compartment of living cells. There is increasing evidence that physical effects, such as molecular crowding or excluded volume effects of transcriptional-translational machinery, affect the yield of reaction products. On the other hand, transcriptional feedback that controls gene expression during mRNA synthesis is also a vital mechanism that regulates protein synthesis in cells. However, the excluded volume effect of spatial constraints on feedback regulation is not well understood. Here, we study the confinement effect on transcriptional autoregulatory feedbacks of gene expression reactions using a theoretical model. The excluded volume effects between molecules and the membrane interface suppress the gene expression in a small cell-sized compartment. We find that negative feedback regulation at the transcription step mitigates this size-induced gene repression and alters the scaling relation of gene expression level on compartment volume, approaching the regular scaling relation without the steric effect. This recovery of regular size-scaling of gene expression does not appear in positive feedback regulation, suggesting that negative autoregulatory feedback is crucial for maintaining reaction products constant regardless of compartment size in heterogeneous cell populations.

Chemical Communication in Artificial Cells: Basic Concepts, Design and Challenges

In the past decade, the focus of bottom-up synthetic biology has shifted from the design of complex artificial cell architectures to the design of interactions between artificial cells mediated by physical and chemical cues. Engineering communication between artificial cells is crucial for the realization of coordinated dynamic behaviours in artificial cell populations, which would have implications for biotechnology, advanced colloidal materials and regenerative medicine. In this review, we focus our discussion on molecular communication between artificial cells. We cover basic concepts such as the importance of compartmentalization, the metabolic machinery driving signaling across cell boundaries and the different modes of communication used. The various studies in artificial cell signaling have been classified based on the distance between sender and receiver cells, just like in biology into autocrine, juxtacrine, paracrine and endocrine signaling. Emerging tools available for the design of dynamic and adaptive signaling are highlighted and some recent advances of signaling-enabled collective behaviours, such as quorum sensing, travelling pulses and predator-prey behaviour, are also discussed.

Universally valid reduction of multiscale stochastic biochemical systems using simple non-elementary propensities

Biochemical systems consist of numerous elementary reactions governed by the law of mass action. However, experimentally characterizing all the elementary reactions is nearly impossible. Thus, over a century, their deterministic models that typically contain rapid reversible bindings have been simplified with non-elementary reaction functions (e.g., Michaelis-Menten and Morrison equations). Although the non-elementary reaction functions are derived by applying the quasi-steady-state approximation (QSSA) to deterministic systems, they have also been widely used to derive propensities for stochastic simulations due to computational efficiency and simplicity. However, the validity condition for this heuristic approach has not been identified even for the reversible binding between molecules, such as protein-DNA, enzyme-substrate, and receptor-ligand, which is the basis for living cells. Here, we find that the non-elementary propensities based on the deterministic total QSSA can accurately capture the stochastic dynamics of the reversible binding in general. However, serious errors occur when reactant molecules with similar levels tightly bind, unlike deterministic systems. In that case, the non-elementary propensities distort the stochastic dynamics of a bistable switch in the cell cycle and an oscillator in the circadian clock. Accordingly, we derive alternative non-elementary propensities with the stochastic low-state QSSA, developed in this study. This provides a universally valid framework for simplifying multiscale stochastic biochemical systems with rapid reversible bindings, critical for efficient stochastic simulations of cell signaling and gene regulation. To facilitate the framework, we provide a user-friendly open-source computational package, ASSISTER, that automatically performs the present framework.

Quantitative control of noise in mammalian gene expression by dynamic histone regulation

Fluctuation ('noise') in gene expression is critical for mammalian cellular processes. Numerous mechanisms contribute to its origins, yet the mechanisms behind large fluctuations that are induced by single transcriptional activators remain elusive. Here, we probed putative mechanisms by studying the dynamic regulation of transcriptional activator binding, histone regulator inhibitors, chromatin accessibility, and levels of mRNAs and proteins in single cells. Using a light-induced expression system, we showed that the transcriptional activator could form an interplay with dual functional co-activator/histone acetyltransferases CBP/p300. This interplay resulted in substantial heterogeneity in H3K27ac, chromatin accessibility, and transcription. Simultaneous attenuation of CBP/p300 and HDAC4/5 reduced heterogeneity in the expression of endogenous genes, suggesting that this mechanism is universal. We further found that the noise was reduced by pulse-wide modulation of transcriptional activator binding possibly as a result of alternating the epigenetic states. Our findings suggest a mechanism for the modulation of noise in synthetic and endogenous gene expression systems.

Topologies of synthetic gene circuit for optimal fold change activation

Computations widely exist in biological systems for functional regulations. Recently, incoherent feedforward loop and integral feedback controller have been implemented into Escherichia coli to achieve a robust adaptation. Here, we demonstrate that an indirect coherent feedforward loop and mutual inhibition designs can experimentally improve the fold change of promoters, by reducing the basal level while keeping the maximum activity high. We applied both designs to six different promoters in E. coli, starting with synthetic inducible promoters as a proof-of-principle. Then, we examined native promoters that are either functionally specific or systemically involved in complex pathways such as oxidative stress and SOS response. Both designs include a cascade having a repressor and a construct of either transcriptional interference or antisense transcription. In all six promoters, an improvement of up to ten times in the fold change activation was observed. Theoretically, our unitless models show that when regulation strength matches promoter basal level, an optimal fold change can be achieved. We expect that this methodology can be applied in various biological systems for biotechnology and therapeutic applications.

Designing Biological Circuits: Synthetic Biology Within the Operon Model and Beyond

In 1961, Jacob and Monod proposed the operon model of gene regulation. At the model's core was the modular assembly of regulators, operators, and structural genes. To illustrate the composability of these elements, Jacob and Monod linked phenotypic diversity to the architectures of regulatory circuits. In this review, we examine how the circuit blueprints imagined by Jacob and Monod laid the foundation for the first synthetic gene networks that launched the field of synthetic biology in 2000. We discuss the influences of the operon model and its broader theoretical framework on the first generation of synthetic biological circuits, which were predominantly transcriptional and posttranscriptional circuits. We also describe how recent advances in molecular biology beyond the operon model-namely, programmable DNA- and RNA-binding molecules as well as models of epigenetic and posttranslational regulation-are expanding the synthetic biology toolkit and enabling the design of more complex biological circuits.

From deterministic to fuzzy decision-making in artificial cells

Building autonomous artificial cells capable of homeostasis requires regulatory networks to gather information and make decisions that take time and cost energy. Decisions based on few molecules may be inaccurate but are cheap and fast. Realizing decision-making with a few molecules in artificial cells has remained a challenge. Here, we show decision-making by a bistable gene network in artificial cells with constant protein turnover. Reducing the number of gene copies from 10⁵ to about 10 per cell revealed a transition from deterministic and slow decision-making to a fuzzy and rapid regime dominated by small-number fluctuations. Gene regulation was observed at lower DNA and protein concentrations than necessary in equilibrium, suggesting rate enhancement by co-expressional localization. The high-copy regime was characterized by a sharp transition and hysteresis, whereas the low-copy limit showed strong fluctuations, state switching, and cellular individuality across the decision-making point. Our results demonstrate information processing with low-power consumption inside artificial cells.

How Do Cells Adapt? Stories Told in Landscapes

Cells adapt to changing environments. Perturb a cell and it returns to a point of homeostasis. Perturb a population and it evolves toward a fitness peak. We review quantitative models of the forces of adaptation and their visualizations on landscapes. While some adaptations result from single mutations or few-gene effects, others are more cooperative, more delocalized in the genome, and more universal and physical. For example, homeostasis and evolution depend on protein folding and aggregation, energy and protein production, protein diffusion, molecular motor speeds and efficiencies, and protein expression levels. Models provide a way to learn about the fitness of cells and cell populations by making and testing hypotheses.

Topology-dependent interference of synthetic gene circuit function by growth feedback

Growth-mediated feedback between synthetic gene circuits and host organisms leads to diverse emerged behaviors, including growth bistability and enhanced ultrasensitivity. However, the range of possible impacts of growth feedback on gene circuits remains underexplored. Here, we mathematically and experimentally demonstrated that growth feedback affects the functions of memory circuits in a network topology-dependent way. Specifically, the memory of the self-activation switch is quickly lost due to the growth-mediated dilution of the circuit products. Decoupling of growth feedback reveals its memory, manifested by its hysteresis property across a broad range of inducer concentration. On the contrary, the toggle switch is more refractory to growth-mediated dilution and can retrieve its memory after the fast-growth phase. The underlying principle lies in the different dependence of active and repressive regulations in these circuits on the growth-mediated dilution. Our results unveil the topology-dependent mechanism on how growth-mediated feedback influences the behaviors of gene circuits.

Hysteretic Genetic Circuit for Detection of Proteasomal Degradation in Mammalian Cells

Synthetic hysteretic mammalian gene circuits generating sustained cellular responses to transient perturbations provide important tools to investigate complex cellular behaviors and reprogram cells for a variety of applications, ranging from protein production to cell fate decisions. The design rules of synthetic gene circuits with controlled hysteretic behaviors, however, remain uncharacterized. To identify the criteria for achieving predictable control of hysteresis, we built a genetic circuit for detection of proteasomal degradation (Hys-Deg). The Hys-Deg circuit is based on a tetracycline-controlled transactivator (tTA) variant engineered to interface with the ubiquitin proteasome system (UPS). The tTA variant activates its own expression, generating a positive feedback loop that is triggered by expression of another tTA gene that is constitutively regulated. Guided by predictive modeling, we characterized the hysteretic response of the Hys-Deg circuit. We demonstrated that control of the hysteretic response is achieved by modulating the ratio of expression of constitutive to inducible tTA. We also showed that the system can be finely tuned through dosage of the inducer tetracycline to calibrate the circuit for detection of the desired levels of UPS activation. This study establishes the design rules for building a hysteretic genetic circuit with an autoregulatory feedback loop and provides a synthetic memory module that could be easily integrated into regulatory gene networks to study and engineer complex cellular behaviors.

Stem cell bioengineering: building from stem cell biology

New fundamental discoveries in stem cell biology have yielded potentially transformative regenerative therapeutics. However, widespread implementation of stem-cell-derived therapeutics remains sporadic. Barriers that impede the development of these therapeutics can be linked to our incomplete understanding of how the regulatory networks that encode stem cell fate govern the development of the complex tissues and organs that are ultimately required for restorative function. Bioengineering tools, strategies and design principles represent core components of the stem cell bioengineering toolbox. Applied to the different layers of complexity present in stem-cell-derived systems - from gene regulatory networks in single stem cells to the systemic interactions of stem-cell-derived organs and tissues - stem cell bioengineering can address existing challenges and advance regenerative medicine and cellular therapies.

Cell population heterogeneity driven by stochastic partition and growth optimality

A fundamental question in biology is how cell populations evolve into different subtypes based on homogeneous processes at the single cell level. Here we show that population bimodality can emerge even when biological processes are homogenous at the cell level and the environment is kept constant. Our model is based on the stochastic partitioning of a cell component with an optimal copy number. We show that the existence of unimodal or bimodal distributions depends on the variance of partition errors and the growth rate tolerance around the optimal copy number. In particular, our theory provides a consistent explanation for the maintenance of aneuploid states in a population. The proposed model can also be relevant for other cell components such as mitochondria and plasmids, whose abundances affect the growth rate and are subject to stochastic partition at cell division.

Circuit-Host Coupling Induces Multifaceted Behavioral Modulations of a Gene Switch

Quantitative modeling of gene circuits is fundamentally important to synthetic biology, as it offers the potential to transform circuit engineering from trial-and-error construction to rational design and, hence, facilitates the advance of the field. Currently, typical models regard gene circuits as isolated entities and focus only on the biochemical processes within the circuits. However, such a standard paradigm is getting challenged by increasing experimental evidence suggesting that circuits and their host are intimately connected, and their interactions can potentially impact circuit behaviors. Here we systematically examined the roles of circuit-host coupling in shaping circuit dynamics by using a self-activating gene switch as a model circuit. Through a combination of deterministic modeling, stochastic simulation, and Fokker-Planck equation formalism, we found that circuit-host coupling alters switch behaviors across multiple scales. At the single-cell level, it slows the switch dynamics in the high protein production regime and enlarges the difference between stable steady-state values. At the population level, it favors cells with low protein production through differential growth amplification. Together, the two-level coupling effects induce both quantitative and qualitative modulations of the switch, with the primary component of the effects determined by the circuit's architectural parameters. This study illustrates the complexity and importance of circuit-host coupling in modulating circuit behaviors, demonstrating the need for a new paradigm-integrated modeling of the circuit-host system-for quantitative understanding of engineered gene networks.

Cytoplasmic Amplification of Transcriptional Noise Generates Substantial Cell-to-Cell Variability

Transcription is an episodic process characterized by probabilistic bursts, but how the transcriptional noise from these bursts is modulated by cellular physiology remains unclear. Using simulations and single-molecule RNA counting, we examined how cellular processes influence cell-to-cell variability (noise). The results show that RNA noise is higher in the cytoplasm than the nucleus in ∼85% of genes across diverse promoters, genomic loci, and cell types (human and mouse). Measurements show further amplification of RNA noise in the cytoplasm, fitting a model of biphasic mRNA conversion between translation- and degradation-competent states. This multi-state translation-degradation of mRNA also causes substantial noise amplification in protein levels, ultimately accounting for ∼74% of intrinsic protein variability in cell populations. Overall, the results demonstrate how noise from transcriptional bursts is intrinsically amplified by mRNA processing, leading to a large super-Poissonian variability in protein levels.

Multiscale effects of heating and cooling on genes and gene networks

Fluctuating environments such as changes in ambient temperature represent a fundamental challenge to life. Cells must protect gene networks that protect them from such stresses, making it difficult to understand how temperature affects gene network function in general. Here, we focus on single genes and small synthetic network modules to reveal four key effects of nonoptimal temperatures at different biological scales: (i) a cell fate choice between arrest and resistance, (ii) slower growth rates, (iii) Arrhenius reaction rates, and (iv) protein structure changes. We develop a multiscale computational modeling framework that captures and predicts all of these effects. These findings promote our understanding of how temperature affects living systems and enables more robust cellular engineering for real-world applications.

Most organisms must cope with temperature changes. This involves genes and gene networks both as subjects and agents of cellular protection, creating difficulties in understanding. Here, we study how heating and cooling affect expression of single genes and synthetic gene circuits in Saccharomyces cerevisiae. We discovered that nonoptimal temperatures induce a cell fate choice between stress resistance and growth arrest. This creates dramatic gene expression bimodality in isogenic cell populations, as arrest abolishes gene expression. Multiscale models incorporating population dynamics, temperature-dependent growth rates, and Arrhenius scaling of reaction rates captured the effects of cooling, but not those of heating in resistant cells. Molecular-dynamics simulations revealed how heating alters the conformational dynamics of the TetR repressor, fully explaining the experimental observations. Overall, nonoptimal temperatures induce a cell fate decision and corrupt gene and gene network function in computationally predictable ways, which may aid future applications of engineered microbes in nonstandard temperatures.

Moment-Based Parameter Estimation for Stochastic Reaction Networks in Equilibrium

Calibrating parameters is a crucial problem within quantitative modeling approaches to reaction networks. Existing methods for stochastic models rely either on statistical sampling or can only be applied to small systems. Here, we present an inference procedure for stochastic models in equilibrium that is based on a moment matching scheme with optimal weighting and that can be used with high-throughput data like the one collected by flow cytometry. Our method does not require an approximation of the underlying equilibrium probability distribution and, if reaction rate constants have to be learned, the optimal values can be computed by solving a linear system of equations. We discuss important practical issues such as the selection of the moments and evaluate the effectiveness of the proposed approach on three case studies.

Computing with biological switches and clocks

The complex dynamics of biological systems is primarily driven by molecular interactions that underpin the regulatory networks of cells. These networks typically contain positive and negative feedback loops, which are responsible for switch-like and oscillatory dynamics, respectively. Many computing systems rely on switches and clocks as computational modules. While the combination of such modules in biological systems leads to a variety of dynamical behaviours, it is also driving development of new computing algorithms. Here we present a historical perspective on computation by biological systems, with a focus on switches and clocks, and discuss parallels between biology and computing. We also outline our vision for the future of biological computing.

Understanding Biological Regulation Through Synthetic Biology

Engineering synthetic gene regulatory circuits proceeds through iterative cycles of design, building, and testing. Initial circuit designs must rely on often-incomplete models of regulation established by fields of reductive inquiry—biochemistry and molecular and systems biology. As differences in designed and experimentally observed circuit behavior are inevitably encountered, investigated, and resolved, each turn of the engineering cycle can force a resynthesis in understanding of natural network function. Here, we outline research that uses the process of gene circuit engineering to advance biological discovery. Synthetic gene circuit engineering research has not only refined our understanding of cellular regulation but furnished biologists with a toolkit that can be directed at natural systems to exact precision manipulation of network structure. As we discuss, using circuit engineering to predictively reorganize, rewire, and reconstruct cellular regulation serves as the ultimate means of testing and understanding how cellular phenotype emerges from systems-level network function.

Underlying Principles of Natural Selection in Network Evolution: Systems Biology Approach

Systems biology is a rapidly expanding field that integrates diverse areas of science such as physics, engineering, computer science, mathematics, and biology toward the goal of elucidating the underlying principles of hierarchical metabolic and regulatory systems in the cell, and ultimately leading to predictive understanding of cellular response to perturbations. Because post-genomics research is taking place throughout the tree of life, comparative approaches offer a way for combining data from many organisms to shed light on the evolution and function of biological networks from the gene to the organismal level. Therefore, systems biology can build on decades of theoretical work in evolutionary biology, and at the same time evolutionary biology can use the systems biology approach to go in new uncharted directions. In this study, we present a review of how the post-genomics era is adopting comparative approaches and dynamic system methods to understand the underlying design principles of network evolution and to shape the nascent field of evolutionary systems biology. Finally, the application of evolutionary systems biology to robust biological network designs is also discussed from the synthetic biology perspective.

On the Adaptive Design Rules of Biochemical Networks in Evolution

Biochemical networks are the backbones of physiological systems of organisms. Therefore, a biochemical network should be sufficiently robust (not sensitive) to tolerate genetic mutations and environmental changes in the evolutionary process. In this study, based on the robustness and sensitivity criteria of biochemical networks, the adaptive design rules are developed for natural selection in the evolutionary process. This will provide insights into the robust adaptive mechanism of biochemical networks in the evolutionary process.

We find that if a mutated biochemical network satisfies the robustness and sensitivity criteria of natural selection, there is a high probability for the biochemical network to prevail during natural selection in the evolutionary process. Since there are various mutated biochemical networks that can satisfy these criteria but have some differences in phenotype, the biochemical networks increase their diversities in the evolutionary process. The robustness of a biochemical network enables co-option so that new phenotypes can be generated in evolution. The proposed robust adaptive design rules of natural selection gain much insight into the evolutionary mechanism and provide a systematic robust biochemical circuit design method of biochemical networks for biotechnological and therapeutic purposes in the future.

Nonlatching positive feedback enables robust bimodality by decoupling expression noise from the mean

Fundamental to biological decision-making is the ability to generate bimodal expression patterns where 2 alternate expression states simultaneously exist. Here, we use a combination of single-cell analysis and mathematical modeling to examine the sources of bimodality in the transcriptional program controlling HIV's fate decision between active replication and viral latency. We find that the HIV transactivator of transcription (Tat) protein manipulates the intrinsic toggling of HIV's promoter, the long terminal repeat (LTR), to generate bimodal ON-OFF expression and that transcriptional positive feedback from Tat shifts and expands the regime of LTR bimodality. This result holds for both minimal synthetic viral circuits and full-length virus. Strikingly, computational analysis indicates that the Tat circuit's noncooperative "nonlatching" feedback architecture is optimized to slow the promoter's toggling and generate bimodality by stochastic extinction of Tat. In contrast to the standard Poisson model, theory and experiment show that nonlatching positive feedback substantially dampens the inverse noise-mean relationship to maintain stochastic bimodality despite increasing mean expression levels. Given the rapid evolution of HIV, the presence of a circuit optimized to robustly generate bimodal expression appears consistent with the hypothesis that HIV's decision between active replication and latency provides a viral fitness advantage. More broadly, the results suggest that positive-feedback circuits may have evolved not only for signal amplification but also for robustly generating bimodality by decoupling expression fluctuations (noise) from mean expression levels.

Control of multidimensional systems on complex network

Multidimensional systems coupled via complex networks are widespread in nature and thus frequently invoked for a large plethora of interesting applications. From ecology to physics, individual entities in mutual interactions are grouped in families, homogeneous in kind. These latter interact selectively, through a sequence of self-consistently regulated steps, whose deeply rooted architecture is stored in the assigned matrix of connections. The asymptotic equilibrium eventually attained by the system, and its associated stability, can be assessed by employing standard nonlinear dynamics tools. For many practical applications, it is however important to externally drive the system towards a desired equilibrium, which is resilient, hence stable, to external perturbations. To this end we here consider a system made up of N interacting populations which evolve according to general rate equations, bearing attributes of universality. One species is added to the pool of interacting families and used as a dynamical controller to induce novel stable equilibria. Use can be made of the root locus method to shape the needed control, in terms of intrinsic reactivity and adopted protocol of injection. The proposed method is tested on both synthetic and real data, thus enabling to demonstrate its robustness and versatility.

Regime shifts driven by dynamic correlations in gene expression noise

Gene expression is a noisy process that leads to regime shifts between alternative steady states among individual living cells, inducing phenotypic variability. The effects of white noise on the regime shift in bistable systems have been well characterized, however little is known about such effects of colored noise (noise with nonzero correlation time). Here, we show that noise correlation time, by considering a genetic circuit of autoactivation, can have a significant effect on the regime shift between distinct phenotypic states in gene expression. We demonstrate this theoretically, using stochastic potential, stationary probability density function, and first-passage time based on the Fokker-Planck description, where the Ornstein-Uhlenbeck process is used to model colored noise. We find that an increase in noise correlation time in the degradation rate can induce a regime shift from a low to a high protein concentration state and enhance the bistable regime, while an increase in noise correlation time in the basal rate retains the bimodal distribution. We then show how cross-correlated colored noises in basal and degradation rates can induce regime shifts from a low to a high protein concentration state, but reduce the bistable regime. We also validate these results through direct numerical simulations of the stochastic differential equation. In gene expression understanding the causes of regime shift to a harmful phenotype could improve early therapeutic intervention in complex human diseases.

Intertangled stochastic motifs in networks of excitatory-inhibitory units

A stochastic model of excitatory and inhibitory interactions which bears universality traits is introduced and studied. The endogenous component of noise, stemming from finite size corrections, drives robust internode correlations that persist at large distances. Antiphase synchrony at small frequencies is resolved on adjacent nodes and found to promote the spontaneous generation of long-ranged stochastic patterns that invade the network as a whole. These patterns are lacking under the idealized deterministic scenario, and could provide hints on how living systems implement and handle a large gallery of delicate computational tasks.

Engineered microbial biosensors based on bacterial two-component systems as synthetic biotechnology platforms in bioremediation and biorefinery

Two-component regulatory systems (TCRSs) mediate cellular response by coupling sensing and regulatory mechanisms. TCRSs are comprised of a histidine kinase (HK), which serves as a sensor, and a response regulator, which regulates expression of the effector gene after being phosphorylated by HK. Using these attributes, bacterial TCRSs can be engineered to design microbial systems for different applications. This review focuses on the current advances in TCRS-based biosensors and on the design of microbial systems for bioremediation and their potential application in biorefinery.

A computational method for the investigation of multistable systems and its application to genetic switches

Background

Genetic switches exhibit multistability, form the basis of epigenetic memory, and are found in natural decision making systems, such as cell fate determination in developmental pathways. Synthetic genetic switches can be used for recording the presence of different environmental signals, for changing phenotype using synthetic inputs and as building blocks for higher-level sequential logic circuits. Understanding how multistable switches can be constructed and how they function within larger biological systems is therefore key to synthetic biology.

Results

Here we present a new computational tool, called StabilityFinder, that takes advantage of sequential Monte Carlo methods to identify regions of parameter space capable of producing multistable behaviour, while handling uncertainty in biochemical rate constants and initial conditions. The algorithm works by clustering trajectories in phase space, and iteratively minimizing a distance metric. Here we examine a collection of models of genetic switches, ranging from the deterministic Gardner toggle switch to stochastic models containing different positive feedback connections. We uncover the design principles behind making bistable, tristable and quadristable switches, and find that rate of gene expression is a key parameter. We demonstrate the ability of the framework to examine more complex systems and examine the design principles of a three gene switch. Our framework allows us to relax the assumptions that are often used in genetic switch models and we show that more complex abstractions are still capable of multistable behaviour.

Conclusions

Our results suggest many ways in which genetic switches can be enhanced and offer designs for the construction of novel switches. Our analysis also highlights subtle changes in correlation of experimentally tunable parameters that can lead to bifurcations in deterministic and stochastic systems. Overall we demonstrate that StabilityFinder will be a valuable tool in the future design and construction of novel gene networks.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-016-0375-z) contains supplementary material, which is available to authorized users.

Stochasticity in the Genotype-Phenotype Map: Implications for the Robustness and Persistence of Bet-Hedging

Nongenetic variation in phenotypes, or bet-hedging, has been observed as a driver of drug resistance in both bacterial infections and cancers. Here, we study how bet-hedging emerges in genotype-phenotype (GP) mapping through a simple interaction model: a molecular switch. We use simple chemical reaction networks to implement stochastic switches that map gene products to phenotypes, and investigate the impact of structurally distinct mappings on the evolution of phenotypic heterogeneity. Bet-hedging naturally emerges within this model, and is robust to evolutionary loss through mutations to both the expression of individual genes, and to the network itself. This robustness explains an apparent paradox of bet-hedging-why does it persist in environments where natural selection necessarily acts to remove it? The structure of the underlying molecular mechanism, itself subject to selection, can slow the evolutionary loss of bet-hedging to ensure a survival mechanism against environmental catastrophes even when they are rare. Critically, these properties, taken together, have profound implications for the use of treatment-holidays to combat bet-hedging-driven resistant disease, as the efficacy of breaks from treatment will ultimately be determined by the structure of the GP mapping.

An Orthogonal Permease-Inducer-Repressor Feedback Loop Shows Bistability

Feedback loops in biological networks, among others, enable differentiation and cell cycle progression, and increase robustness in signal transduction. In natural networks, feedback loops are often complex and intertwined, making it challenging to identify which loops are mainly responsible for an observed behavior. However, minimal synthetic replicas could allow for such identification. Here, we engineered a synthetic permease-inducer-repressor system in Saccharomyces cerevisiae to analyze if a transport-mediated positive feedback loop could be a core mechanism for the switch-like behavior in the regulation of metabolic gene networks such as the S. cerevisiae GAL system or the Escherichia coli lac operon. We characterized the synthetic circuit using deterministic and stochastic mathematical models. Similar to its natural counterparts, our synthetic system shows bistable and hysteretic behavior, and the inducer concentration range for bistability as well as the switching rates between the two stable states depend on the repressor concentration. Our results indicate that a generic permease-inducer-repressor circuit with a single feedback loop is sufficient to explain the experimentally observed bistable behavior of the natural systems. We anticipate that the approach of reimplementing natural systems with orthogonal parts to identify crucial network components is applicable to other natural systems such as signaling pathways.

Challenges in structural approaches to cell modeling

Computational modeling is essential for structural characterization of biomolecular mechanisms across the broad spectrum of scales. Adequate understanding of biomolecular mechanisms inherently involves our ability to model them. Structural modeling of individual biomolecules and their interactions has been rapidly progressing. However, in terms of the broader picture, the focus is shifting toward larger systems, up to the level of a cell. Such modeling involves a more dynamic and realistic representation of the interactomes in vivo, in a crowded cellular environment, as well as membranes and membrane proteins, and other cellular components. Structural modeling of a cell complements computational approaches to cellular mechanisms based on differential equations, graph models, and other techniques to model biological networks, imaging data, etc. Structural modeling along with other computational and experimental approaches will provide a fundamental understanding of life at the molecular level and lead to important applications to biology and medicine. A cross section of diverse approaches presented in this review illustrates the developing shift from the structural modeling of individual molecules to that of cell biology. Studies in several related areas are covered: biological networks; automated construction of three-dimensional cell models using experimental data; modeling of protein complexes; prediction of non-specific and transient protein interactions; thermodynamic and kinetic effects of crowding; cellular membrane modeling; and modeling of chromosomes. The review presents an expert opinion on the current state-of-the-art in these various aspects of structural modeling in cellular biology, and the prospects of future developments in this emerging field.

Feedback control indirect response models

A general framework is introduced for modeling pharmacodynamic processes that are subject to autoregulation, which combines the indirect response (IDR) model approach with methods from classical feedback control of engineered systems. The canonical IDR models are modified to incorporate linear combinations of feedback control terms related to the time course of the difference (the error signal) between the pharmacodynamic response and its basal value. Following the well-established approach of traditional engineering control theory, the proposed feedback control indirect response models incorporate terms proportional to the error signal itself, the integral of the error signal, the derivative of the error signal or combinations thereof. Simulations are presented to illustrate the types of responses produced by the proposed feedback control indirect response model framework, and to illustrate comparisons with other PK/PD modeling approaches incorporating feedback. In addition, four examples from literature are used to illustrate the implementation and applicability of the proposed feedback control framework. The examples reflect each of the four mechanisms of drug action as modeled by each of the four canonical IDR models and include: selective serotonin reuptake inhibitors and extracellular serotonin; histamine H2-receptor antagonists and gastric acid; growth hormone secretagogues and circulating growth hormone; β2-selective adrenergic agonists and potassium. The proposed feedback control indirect response approach may serve as an exploratory modeling tool and may provide a bridge for development of more mechanistic systems pharmacology models.

Synchronizing Genetic Oscillators by Signaling Molecules

The authors examine collective rhythms in a general multicell system with both linearly diffusive and nondiffusive couplings. The effect of coupling on synchronization through intercellular signaling in a population of Escherichia coli cells is studied. In particular, a synchronization solution is given through the auxiliary individual system for 2 types of couplings. The sufficient conditions for the global synchronization of such a coupled system are derived based on the Lyapunov function method. The authors show that an appropriate design of the coupling and the inner-linking matrix can ensure global synchronization of the coupled synthetic biological system. Moreover, they demonstrate that the dynamics of an individual cell with coupling and without coupling may be qualitatively different; one is oscillatory, and the other is steady state. The change from a nonoscillatory state to an oscillatory one is induced by appropriate coupling, which also entrains all cells to synchronization. These results establish not only a theoretical foundation but also a quantitative basis for understanding the essential cooperative dynamics, such as collective rhythms or synchronization, in a population of cells.

Nonlinear Dynamics in Gene Regulation Promote Robustness and Evolvability of Gene Expression Levels

Cellular phenotypes underpinned by regulatory networks need to respond to evolutionary pressures to allow adaptation, but at the same time be robust to perturbations. This creates a conflict in which mutations affecting regulatory networks must both generate variance but also be tolerated at the phenotype level. Here, we perform mathematical analyses and simulations of regulatory networks to better understand the potential trade-off between robustness and evolvability. Examining the phenotypic effects of mutations, we find an inverse correlation between robustness and evolvability that breaks only with nonlinearity in the network dynamics, through the creation of regions presenting sudden changes in phenotype with small changes in genotype. For genotypes embedding low levels of nonlinearity, robustness and evolvability correlate negatively and almost perfectly. By contrast, genotypes embedding nonlinear dynamics allow expression levels to be robust to small perturbations, while generating high diversity (evolvability) under larger perturbations. Thus, nonlinearity breaks the robustness-evolvability trade-off in gene expression levels by allowing disparate responses to different mutations. Using analytical derivations of robustness and system sensitivity, we show that these findings extend to a large class of gene regulatory network architectures and also hold for experimentally observed parameter regimes. Further, the effect of nonlinearity on the robustness-evolvability trade-off is ensured as long as key parameters of the system display specific relations irrespective of their absolute values. We find that within this parameter regime genotypes display low and noisy expression levels. Examining the phenotypic effects of mutations, we find an inverse correlation between robustness and evolvability that breaks only with nonlinearity in the network dynamics. Our results provide a possible solution to the robustness-evolvability trade-off, suggest an explanation for the ubiquity of nonlinear dynamics in gene expression networks, and generate useful guidelines for the design of synthetic gene circuits.

Phenotypic Variability in Synthetic Biology Applications: Dealing with Noise in Microbial Gene Expression

The stochasticity due to the infrequent collisions among low copy-number molecules within the crowded cellular compartment is a feature of living systems. Single cell variability in gene expression within an isogenic population (i.e., biological noise) is usually described as the sum of two independent components: intrinsic and extrinsic stochasticity. Intrinsic stochasticity arises from the random occurrence of events inherent to the gene expression process (e.g., the burst-like synthesis of mRNA and protein molecules). Extrinsic fluctuations reflect the state of the biological system and its interaction with the intra and extracellular environments (e.g., concentration of available polymerases, ribosomes, metabolites, and micro-environmental conditions). A better understanding of cellular noise would help synthetic biologists design gene circuits with well-defined functional properties. In silico modeling has already revealed several aspects of the network topology’s impact on noise properties; this information could drive the selection of biological parts and the design of reliably engineered pathways. Importantly, while optimizing artificial gene circuitry for industrial applications, synthetic biology could also elucidate the natural mechanisms underlying natural phenotypic variability. In this review, we briefly summarize the functional roles of noise in unicellular organisms and address their relevance to synthetic network design. We will also consider how noise might influence the selection of network topologies supporting reliable functions, and how the variability of cellular events might be exploited when designing innovative biotechnology applications.

Synthetic biology: insights into biological computation

Synthetic biology attempts to rationally engineer biological systems in order to perform desired functions. Our increasing understanding of biological systems guides this rational design, while the huge background in electronics for building circuits defines the methodology.

Robust, tunable genetic memory from protein sequestration combined with positive feedback

Natural regulatory networks contain many interacting components that allow for fine-tuning of switching and memory properties. Building simple bistable switches, synthetic biologists have learned the design principles of complex natural regulatory networks. However, most switches constructed so far are so simple (e.g. comprising two regulators) that they are functional only within a limited parameter range. Here, we report the construction of robust, tunable bistable switches in Escherichia coli using three heterologous protein regulators (ExsADC) that are sequestered into an inactive complex through a partner swapping mechanism. On the basis of mathematical modeling, we accurately predict and experimentally verify that the hysteretic region can be fine-tuned by controlling the interactions of the ExsADC regulatory cascade using the third member ExsC as a tuning knob. Additionally, we confirm that a dual-positive feedback switch can markedly increase the hysteretic region, compared to its single-positive feedback counterpart. The dual-positive feedback switch displays bistability over a 10(6)-fold range of inducer concentrations, to our knowledge, the largest range reported so far. This work demonstrates the successful interlocking of sequestration-based ultrasensitivity and positive feedback, a design principle that can be applied to the construction of robust, tunable, and predictable genetic programs to achieve increasingly sophisticated biological behaviors.

Biological Sources of Intrinsic and Extrinsic Noise in cI Expression of Lysogenic Phage Lambda

Genetically identical cells exposed to homogeneous environment can show remarkable phenotypic difference. To predict how phenotype is shaped, understanding of how each factor contributes is required. During gene expression processes, noise could arise either intrinsically in biochemical processes of gene expression or extrinsically from other cellular processes such as cell growth. In this work, important noise sources in gene expression of phage λ lysogen are quantified using models described by stochastic differential equations (SDEs). Results show that DNA looping has sophisticated impacts on gene expression noise: When DNA looping provides autorepression, like in wild type, it reduces noise in the system; When the autorepression is defected as it is in certain mutants, DNA looping increases expression noise. We also study how each gene operator affects the expression noise by changing the binding affinity between the gene and the transcription factor systematically. We find that the system shows extraordinarily large noise when the binding affinity is in certain range, which changes the system from monostable to bistable. In addition, we find that cell growth causes non-negligible noise, which increases with gene expression level. Quantification of noise and identification of new noise sources will provide deeper understanding on how stochasticity impacts phenotype.

Are Quasi-Steady-State Approximated Models Suitable for Quantifying Intrinsic Noise Accurately?

Large gene regulatory networks (GRN) are often modeled with quasi-steady-state approximation (QSSA) to reduce the huge computational time required for intrinsic noise quantification using Gillespie stochastic simulation algorithm (SSA). However, the question still remains whether the stochastic QSSA model measures the intrinsic noise as accurately as the SSA performed for a detailed mechanistic model or not? To address this issue, we have constructed mechanistic and QSSA models for few frequently observed GRNs exhibiting switching behavior and performed stochastic simulations with them. Our results strongly suggest that the performance of a stochastic QSSA model in comparison to SSA performed for a mechanistic model critically relies on the absolute values of the mRNA and protein half-lives involved in the corresponding GRN. The extent of accuracy level achieved by the stochastic QSSA model calculations will depend on the level of bursting frequency generated due to the absolute value of the half-life of either mRNA or protein or for both the species. For the GRNs considered, the stochastic QSSA quantifies the intrinsic noise at the protein level with greater accuracy and for larger combinations of half-life values of mRNA and protein, whereas in case of mRNA the satisfactory accuracy level can only be reached for limited combinations of absolute values of half-lives. Further, we have clearly demonstrated that the abundance levels of mRNA and protein hardly matter for such comparison between QSSA and mechanistic models. Based on our findings, we conclude that QSSA model can be a good choice for evaluating intrinsic noise for other GRNs as well, provided we make a rational choice based on experimental half-life values available in literature.

Phenotypic states become increasingly sensitive to perturbations near a bifurcation in a synthetic gene network

Microorganisms often exhibit a history-dependent phenotypic response after exposure to a stimulus which can be imperative for proper function. However, cells frequently experience unexpected environmental perturbations that might induce phenotypic switching. How cells maintain phenotypic states in the face of environmental fluctuations remains an open question. Here, we use environmental perturbations to characterize the resilience of phenotypic states in a synthetic gene network near a critical transition. We find that far from the critical transition an environmental perturbation may induce little to no phenotypic switching, whereas close to the critical transition the same perturbation can cause many cells to switch phenotypic states. This loss of resilience was observed for perturbations that interact directly with the gene circuit as well as for a variety of generic perturbations-such as salt, ethanol, or temperature shocks-that alter the state of the cell more broadly. We obtain qualitatively similar findings in natural gene circuits, such as the yeast GAL network. Our findings illustrate how phenotypic memory can become destabilized by environmental variability near a critical transition.

DOI:http://dx.doi.org/10.7554/eLife.07935.001

All organisms need to be able to react to the challenges thrown at them by their changing environment. Yeast, bacteria and other microbes have networks of genes that can give rise to many different traits and characteristics, which can also be referred to as phenotypes. A change in the environment can alter the activities' of the genes so that the microbes display a different phenotype. The point at which a small change in the environment can lead to a sudden switch in the phenotype is called a ‘critical transition’.

An individual microbe's history can influence the phenotype that it presents. However, it is not clear how microbes ‘remember’ their history, or how fluctuations in the environment might cause the microbe to lose the ability to store this memory and present a different phenotype instead. Here, Axelrod et al. studied phenotype memory in yeast cells grown in the laboratory. The experiments used cells that had been genetically modified to glow red in the presence of a molecule called anhydrotetracycline (or ATc) and to glow green in the absence of the molecule.

Axelrod et al. examined what effect altering the levels of this molecule would have on the phenotype produced by the cells. First, the cells were grown with no ATc present for several generations so that the cells glowed green. Next, Axelrod et al. added different amounts of ATc were added. For moderate levels of ATc the cells continued to glow green, illustrating that they ‘remembered’ their prior growth condition. However, cells exposed to higher levels of ATc lost this memory and changed color.

Next, Axelrod et al. carried out further experiments on cells exposed to ATc levels that were close to, or further away from the critical transition. At high levels of ATc (that is, close to the critical transition), many cells switched from green to red when exposed to high temperatures, salt and other changes in the environment. On the other hand, very few of the cells grown in low levels of ATc—and therefore further away from the critical transition—changed color in response to the same fluctuations in their environment.

These finding reveal that phenotype memory is less stable when yeast experience fluctuations in their environment close to a critical transition. Future work will seek to find out how salt or high temperatures can abolish phenotype memory.

DOI:http://dx.doi.org/10.7554/eLife.07935.002

Using Nonlinear Stochastic Evolutionary Game Strategy to Model an Evolutionary Biological Network of Organ Carcinogenesis Under a Natural Selection Scheme

Molecular biologists have long recognized carcinogenesis as an evolutionary process that involves natural selection. Cancer is driven by the somatic evolution of cell lineages. In this study, the evolution of somatic cancer cell lineages during carcinogenesis was modeled as an equilibrium point (ie, phenotype of attractor) shifting, the process of a nonlinear stochastic evolutionary biological network. This process is subject to intrinsic random fluctuations because of somatic genetic and epigenetic variations, as well as extrinsic disturbances because of carcinogens and stressors. In order to maintain the normal function (ie, phenotype) of an evolutionary biological network subjected to random intrinsic fluctuations and extrinsic disturbances, a network robustness scheme that incorporates natural selection needs to be developed. This can be accomplished by selecting certain genetic and epigenetic variations to modify the network structure to attenuate intrinsic fluctuations efficiently and to resist extrinsic disturbances in order to maintain the phenotype of the evolutionary biological network at an equilibrium point (attractor). However, during carcinogenesis, the remaining (or neutral) genetic and epigenetic variations accumulate, and the extrinsic disturbances become too large to maintain the normal phenotype at the desired equilibrium point for the nonlinear evolutionary biological network. Thus, the network is shifted to a cancer phenotype at a new equilibrium point that begins a new evolutionary process. In this study, the natural selection scheme of an evolutionary biological network of carcinogenesis was derived from a robust negative feedback scheme based on the nonlinear stochastic Nash game strategy. The evolvability and phenotypic robustness criteria of the evolutionary cancer network were also estimated by solving a Hamilton-Jacobi inequality - constrained optimization problem. The simulation revealed that the phenotypic shift of the lung cancer-associated cell network takes 54.5 years from a normal state to stage I cancer, 1.5 years from stage I to stage II cancer, and 2.5 years from stage II to stage III cancer, with a reasonable match for the statistical result of the average age of lung cancer. These results suggest that a robust negative feedback scheme, based on a stochastic evolutionary game strategy, plays a critical role in an evolutionary biological network of carcinogenesis under a natural selection scheme.

Roles of factorial noise in inducing bimodal gene expression

Some gene regulatory systems can exhibit bimodal distributions of mRNA or protein although the deterministic counterparts are monostable. This noise-induced bimodality is an interesting phenomenon and has important biological implications, but it is unclear how different sources of expression noise (each source creates so-called factorial noise that is defined as a component of the total noise) contribute separately to this stochastic bimodality. Here we consider a minimal model of gene regulation, which is monostable in the deterministic case. Although simple, this system contains factorial noise of two main kinds: promoter noise due to switching between gene states and transcriptional (or translational) noise due to synthesis and degradation of mRNA (or protein). To better trace the roles of factorial noise in inducing bimodality, we also analyze two limit models, continuous and adiabatic approximations, apart from the exact model. We show that in the case of slow gene switching, the continuous model where only promoter noise is considered can exhibit bimodality; in the case of fast switching, the adiabatic model where only transcriptional or translational noise is considered can also exhibit bimodality but the exact model cannot; and in other cases, both promoter noise and transcriptional or translational noise can cooperatively induce bimodality. Since slow gene switching and large protein copy numbers are characteristics of eukaryotic cells, whereas fast gene switching and small protein copy numbers are characteristics of prokaryotic cells, we infer that eukaryotic stochastic bimodality is induced mainly by promoter noise, whereas prokaryotic stochastic bimodality is induced primarily by transcriptional or translational noise.

Autoregulation and heterogeneity in expression of human Cripto-1

Cripto-1 (CR-1) is involved in various processes in embryonic development and cancer. Multiple pathways regulate CR-1 expression. Our present work demonstrates a possible positive feedback circuit where CR-1 induces its own expression. Using U-87 MG cells treated with exogenous CR-1, we show that such induction involves ALK4/SMAD2/3 pathway. Stochasticity in gene expression gives rise to heterogeneity in expression in genetically identical cells. Positive feedback increases such heterogeneity and often gives rise to two subpopulations of cells, having higher and lower expression of a gene. Using flow cytometry, we show that U-87 MG cells have a minuscule subpopulation with detectable expression of CR-1. Induction of CR-1 expression, by exogenous CR-1, increases the size of this CR-1 positive subpopulation. However, even at very high dose, most of the cells remain CR-1 negative. We show that population behavior of CR-1 induction has a signature similar to bimodal expression expected in a transcriptional circuit with positive feedback. We further show that treatment of U-87 MG cells with CR-1 leads to higher expression of drug efflux protein MDR-1 in the CR-1 positive subpopulation, indicating correlated induction of these two proteins. Positive feedback driven heterogeneity in expression of CR-1 may play crucial role in phenotypic diversification of cancer cells.