Inverse bifurcation analysis: application to simple gene systems

BACKGROUND: Bifurcation analysis has proven to be a powerful method for understanding the qualitative behavior of gene regulatory networks. In addition to the more traditional forward problem of determining the mapping from parameter space to the space of model behavior, the inverse problem of determining model parameters to result in certain desired properties of the bifurcation diagram provides an attractive methodology for addressing important biological problems. These include understanding how the robustness of qualitative behavior arises from system design as well as providing a way to engineer biological networks with qualitative properties. RESULTS: We demonstrate that certain inverse bifurcation problems of biological interest may be cast as optimization problems involving minimal distances of reference parameter sets to bifurcation manifolds. This formulation allows for an iterative solution procedure based on performing a sequence of eigen-system computations and one-parameter continuations of solutions, the latter being a standard capability in existing numerical bifurcation software. As applications of the proposed method, we show that the problem of maximizing regions of a given qualitative behavior as well as the reverse engineering of bistable gene switches can be modelled and efficiently solved

Timing Cellular Decision Making Under Noise via Cell–Cell Communication

Many cellular processes require decision making mechanisms, which must act reliably even in the unavoidable presence of substantial amounts of noise. However, the multistable genetic switches that underlie most decision-making processes are dominated by fluctuations that can induce random jumps between alternative cellular states. Here we show, via theoretical modeling of a population of noise-driven bistable genetic switches, that reliable timing of decision-making processes can be accomplished for large enough population sizes, as long as cells are globally coupled by chemical means. In the light of these results, we conjecture that cell proliferation, in the presence of cell–cell communication, could provide a mechanism for reliable decision making in the presence of noise, by triggering cellular transitions only when the whole cell population reaches a certain size. In other words, the summation performed by the cell population would average out the noise and reduce its detrimental impact

A model for improving microbial biofuel production using a synthetic feedback loop

Cells use feedback to implement a diverse range of regulatory functions. Building synthetic feedback control systems may yield insight into the roles that feedback can play in regulation since it can be introduced independently of native regulation, and alternative control architectures can be compared. We propose a model for microbial biofuel production where a synthetic control system is used to increase cell viability and biofuel yields. Although microbes can be engineered to produce biofuels, the fuels are often toxic to cell growth, creating a negative feedback loop that limits biofuel production. These toxic effects may be mitigated by expressing efflux pumps that export biofuel from the cell. We developed a model for cell growth and biofuel production and used it to compare several genetic control strategies for their ability to improve biofuel yields. We show that controlling efflux pump expression directly with a biofuel-responsive promoter is a straightforward way of improving biofuel production. In addition, a feed forward loop controller is shown to be versatile at dealing with uncertainty in biofuel production rates

Stochastic analysis of the GAL genetic switch in Saccharomyces cerevisiae: Modeling and experiments reveal hierarchy in glucose repression

<p>Abstract</p> <p>Background</p> <p>Transcriptional regulation involves protein-DNA and protein-protein interactions. Protein-DNA interactions involve reactants that are present in low concentrations, leading to stochastic behavior. In addition, multiple regulatory mechanisms are typically involved in transcriptional regulation. In the <it>GAL </it>regulatory system of <it>Saccharomyces cerevisiae</it>, the inhibition of glucose is accomplished through two regulatory mechanisms: one through the transcriptional repressor Mig1p, and the other through regulating the amount of transcriptional activator Gal4p. However, the impact of stochasticity in gene expression and hierarchy in regulatory mechanisms on the phenotypic level is not clearly understood.</p> <p>Results</p> <p>We address the question of quantifying the effect of stochasticity inherent in these regulatory mechanisms on the performance of various genes under the regulation of Mig1p and Gal4p using a dynamic stochastic model. The stochastic analysis reveals the importance of both the mechanisms of regulation for tight expression of genes in the <it>GAL </it>network. The mechanism involving Gal4p is the dominant mechanism, yielding low variability in the expression of <it>GAL </it>genes. The mechanism involving Mig1p is necessary to maintain the switch-like response of certain <it>GAL </it>genes. The number of binding sites for Mig1p and Gal4p further influences the expression of the genes, with extra binding sites lowering the variability of expression. Our experiments involving growth on various substrates show that the trends predicted in mean expression and its variability are transmitted to the phenotypic level.</p> <p>Conclusion</p> <p>The mechanisms involved in the transcriptional regulation and their variability set up a hierarchy in the phenotypic response to growth on various substrates. Structural motifs, such as the number of binding sites and the mechanism of regulation, determine the level of stochasticity and eventually, the phenotypic response.</p

Asynchronous combinatorial action of four regulatory factors activates Bcl11b for T cell commitment

During T cell development, multipotent progenitors relinquish competence for other fates and commit to the T cell lineage by turning on Bcl11b, which encodes a transcription factor. To clarify lineage commitment mechanisms, we followed developing T cells at the single-cell level using Bcl11b knock-in fluorescent reporter mice. Notch signaling and Notch-activated transcription factors collaborate to activate Bcl11b expression irrespectively of Notch-dependent proliferation. These inputs work via three distinct, asynchronous mechanisms: an early locus 'poising' function dependent on TCF-1 and GATA-3, a stochastic-permissivity function dependent on Notch signaling, and a separate amplitude-control function dependent on Runx1, a factor already present in multipotent progenitors. Despite their necessity for Bcl11b expression, these inputs act in a stage-specific manner, providing a multitiered mechanism for developmental gene regulation

Noise Reduction by Diffusional Dissipation in a Minimal Quorum Sensing Motif

Cellular interactions are subject to random fluctuations (noise) in quantities of interacting molecules. Noise presents a major challenge for the robust function of natural and engineered cellular networks. Past studies have analyzed how noise is regulated at the intracellular level. Cell–cell communication, however, may provide a complementary strategy to achieve robust gene expression by enabling the coupling of a cell with its environment and other cells. To gain insight into this issue, we have examined noise regulation by quorum sensing (QS), a mechanism by which many bacteria communicate through production and sensing of small diffusible signals. Using a stochastic model, we analyze a minimal QS motif in Gram-negative bacteria. Our analysis shows that diffusion of the QS signal, together with fast turnover of its transcriptional regulator, attenuates low-frequency components of extrinsic noise. We term this unique mechanism “diffusional dissipation” to emphasize the importance of fast signal turnover (or dissipation) by diffusion. We further show that this noise attenuation is a property of a more generic regulatory motif, of which QS is an implementation. Our results suggest that, in a QS system, an unstable transcriptional regulator may be favored for regulating expression of costly proteins that generate public goods

Genome-wide analysis of differential transcriptional and epigenetic variability across human immune cell types

Abstract Background A healthy immune system requires immune cells that adapt rapidly to environmental challenges. This phenotypic plasticity can be mediated by transcriptional and epigenetic variability. Results We apply a novel analytical approach to measure and compare transcriptional and epigenetic variability genome-wide across CD14+CD16− monocytes, CD66b+CD16+ neutrophils, and CD4+CD45RA+ naïve T cells from the same 125 healthy individuals. We discover substantially increased variability in neutrophils compared to monocytes and T cells. In neutrophils, genes with hypervariable expression are found to be implicated in key immune pathways and are associated with cellular properties and environmental exposure. We also observe increased sex-specific gene expression differences in neutrophils. Neutrophil-specific DNA methylation hypervariable sites are enriched at dynamic chromatin regions and active enhancers. Conclusions Our data highlight the importance of transcriptional and epigenetic variability for the key role of neutrophils as the first responders to inflammatory stimuli. We provide a resource to enable further functional studies into the plasticity of immune cells, which can be accessed from: http://blueprint-dev.bioinfo.cnio.es/WP10/hypervariability

Consequences of Eukaryotic Enhancer Architecture for Gene Expression Dynamics, Development, and Fitness

The regulatory logic of time- and tissue-specific gene expression has mostly been dissected in the context of the smallest DNA fragments that, when isolated, recapitulate native expression in reporter assays. It is not known if the genomic sequences surrounding such fragments, often evolutionarily conserved, have any biological function or not. Using an enhancer of the even-skipped gene of Drosophila as a model, we investigate the functional significance of the genomic sequences surrounding empirically identified enhancers. A 480 bp long “minimal stripe element” is able to drive even-skipped expression in the second of seven stripes but is embedded in a larger region of 800 bp containing evolutionarily conserved binding sites for required transcription factors. To assess the overall fitness contribution made by these binding sites in the native genomic context, we employed a gene-replacement strategy in which whole-locus transgenes, capable of rescuing even-skipped- lethality to adulthood, were substituted for the native gene. The molecular phenotypes were characterized by tagging Even-skipped with a fluorescent protein and monitoring gene expression dynamics in living embryos. We used recombineering to excise the sequences surrounding the minimal enhancer and site-specific transgenesis to create co-isogenic strains differing only in their stripe 2 sequences. Remarkably, the flanking sequences were dispensable for viability, proving the sufficiency of the minimal element for biological function under normal conditions. These sequences are required for robustness to genetic and environmental perturbation instead. The mutant enhancers had measurable sex- and dose-dependent effects on viability. At the molecular level, the mutants showed a destabilization of stripe placement and improper activation of downstream genes. Finally, we demonstrate through live measurements that the peripheral sequences are required for temperature compensation. These results imply that seemingly redundant regulatory sequences beyond the minimal enhancer are necessary for robust gene expression and that “robustness” itself must be an evolved characteristic of the wild-type enhancer

Control of Stochastic Gene Expression by Host Factors at the HIV Promoter

The HIV promoter within the viral long terminal repeat (LTR) orchestrates many aspects of the viral life cycle, from the dynamics of viral gene expression and replication to the establishment of a latent state. In particular, after viral integration into the host genome, stochastic fluctuations in viral gene expression amplified by the Tat positive feedback loop can contribute to the formation of either a productive, transactivated state or an inactive state. In a significant fraction of cells harboring an integrated copy of the HIV-1 model provirus (LTR-GFP-IRES-Tat), this bimodal gene expression profile is dynamic, as cells spontaneously and continuously flip between active (Bright) and inactive (Off) expression modes. Furthermore, these switching dynamics may contribute to the establishment and maintenance of proviral latency, because after viral integration long delays in gene expression can occur before viral transactivation. The HIV-1 promoter contains cis-acting Sp1 and NF-κB elements that regulate gene expression via the recruitment of both activating and repressing complexes. We hypothesized that interplay in the recruitment of such positive and negative factors could modulate the stability of the Bright and Off modes and thereby alter the sensitivity of viral gene expression to stochastic fluctuations in the Tat feedback loop. Using model lentivirus variants with mutations introduced in the Sp1 and NF-κB elements, we employed flow cytometry, mRNA quantification, pharmacological perturbations, and chromatin immunoprecipitation to reveal significant functional differences in contributions of each site to viral gene regulation. Specifically, the Sp1 sites apparently stabilize both the Bright and the Off states, such that their mutation promotes noisy gene expression and reduction in the regulation of histone acetylation and deacetylation. Furthermore, the NF-κB sites exhibit distinct properties, with κB site I serving a stronger activating role than κB site II. Moreover, Sp1 site III plays a particularly important role in the recruitment of both p300 and RelA to the promoter. Finally, analysis of 362 clonal cell populations infected with the viral variants revealed that mutations in any of the Sp1 sites yield a 6-fold higher frequency of clonal bifurcation compared to that of the wild-type promoter. Thus, each Sp1 and NF-κB site differentially contributes to the regulation of viral gene expression, and Sp1 sites functionally “dampen” transcriptional noise and thereby modulate the frequency and maintenance of this model of viral latency. These results may have biomedical implications for the treatment of HIV latency