Cheol-Min Ghim

Genetic noise control via protein oligomerization

BackgroundGene expression in a cell entails random reaction events occurring over disparate time scales. Thus, molecular noise that often results in phenotypic and population-dynamic consequences sets a fundamental limit to biochemical signaling. While there have been numerous studies correlating the architecture of cellular reaction networks with noise tolerance, only a limited effort has been made to understand the dynamic role of protein-protein interactions.ResultsWe have developed a fully stochastic model for the positive feedback control of a single gene, as well as a pair of genes (toggle switch), integrating quantitative results from previous in vivo and in vitro studies. In particular, we explicitly account for the fast binding-unbinding kinetics among proteins, RNA polymerases, and the promoter/operator sequences of DNA. We find that the overall noise-level is reduced and the frequency content of the noise is dramatically shifted to the physiologically irrelevant high-frequency regime in the presence of protein dimerization. This is independent of the choice of monomer or dimer as transcription factor and persists throughout the multiple model topologies considered. For the toggle switch, we additionally find that the presence of a protein dimer, either homodimer or heterodimer, may significantly reduce its random switching rate. Hence, the dimer promotes the robust function of bistable switches by preventing the uninduced (induced) state from randomly being induced (uninduced).ConclusionThe specific binding between regulatory proteins provides a buffer that may prevent the propagation of fluctuations in genetic activity. The capacity of the buffer is a non-monotonic function of association-dissociation rates. Since the protein oligomerization per se does not require extra protein components to be expressed, it provides a basis for the rapid control of intrinsic or extrinsic noise. The stabilization of regulatory circuits and epigenetic memory in general is of direct implications to organism fitness. Our results also suggest possible avenues for the design of synthetic gene circuits with tunable robustness for a wide range of engineering purposes.

The Fixation Probability of Rare Mutators in Finite Asexual Populations

A mutator is an allele that increases the mutation rate throughout the genome by disrupting some aspect of DNA replication or repair. Mutators that increase the mutation rate by the order of 100-fold have been observed to spontaneously emerge and achieve high frequencies in natural populations and in long-term laboratory evolution experiments with Escherichia coli. In principle, the fixation of mutator alleles is limited by (i) competition with mutations in wild-type backgrounds, (ii) additional deleterious mutational load, and (iii) random genetic drift. Using a multiple-locus model and employing both simulation and analytic methods, we investigate the effects of these three factors on the fixation probability Pfix of an initially rare mutator as a function of population size N, beneficial and deleterious mutation rates, and the strength of mutations s. Our diffusion-based approximation for Pfix successfully captures effects ii and iii when selection is fast compared to mutation (\batchmode \documentclass[fleqn,10pt,legalpaper]{article} \usepackage{amssymb} \usepackage{amsfonts} \usepackage{amsmath} \pagestyle{empty} \begin{document} \(\mathrm{{\mu}}/s{\ll}1\) \end{document}). This enables us to predict the conditions under which mutators will be evolutionarily favored. Surprisingly, our simulations show that effect i is typically small for strong-effect mutators. Our results agree semiquantitatively with existing laboratory evolution experiments and suggest future experimental directions.

Systems biology of Microbial Communities

Microbes exist naturally in a wide range of environments in communities where their interactions are significant, spanning the extremes of high acidity and high temperature environments to soil and the ocean. We present a practical discussion of three different approaches for modeling microbial communities: rate equations, individual-based modeling, and population dynamics. We illustrate the approaches with detailed examples. Each approach is best fit to different levels of system representation, and they have different needs for detailed biological input. Thus, this set of approaches is able to address the operation and function of microbial communities on a wide range of organizational levels.

Selected heterozygosity at cis-regulatory sequences increases the expression homogeneity of a cell population in humans

BackgroundExamples of heterozygote advantage in humans are scarce and limited to protein-coding sequences. Here, we attempt a genome-wide functional inference of advantageous heterozygosity at cis-regulatory regions.ResultsThe single-nucleotide polymorphisms bearing the signatures of balancing selection are enriched in active cis-regulatory regions of immune cells and epithelial cells, the latter of which provide barrier function and innate immunity. Examples associated with ancient trans-specific balancing selection are also discovered. Allelic imbalance in chromatin accessibility and divergence in transcription factor motif sequences indicate that these balanced polymorphisms cause distinct regulatory variation. However, a majority of these variants show no association with the expression level of the target gene. Instead, single-cell experimental data for gene expression and chromatin accessibility demonstrate that heterozygous sequences can lower cell-to-cell variability in proportion to selection strengths. This negative correlation is more pronounced for highly expressed genes and consistently observed when using different data and methods. Based on mathematical modeling, we hypothesize that extrinsic noise from fluctuations in transcription factor activity may be amplified in homozygotes, whereas it is buffered in heterozygotes. While high expression levels are coupled with intrinsic noise reduction, regulatory heterozygosity can contribute to the suppression of extrinsic noise.ConclusionsThis mechanism may confer a selective advantage by increasing cell population homogeneity and thereby enhancing the collective action of the cells, especially of those involved in the defense systems in humans.

Hypothetical biomolecular probe based on a genetic switch with tunable symmetry and stability

BackgroundGenetic switches are ubiquitous in nature, frequently associated with the control of cellular functions and developmental programs. In the realm of synthetic biology, it is of great interest to engineer genetic circuits that can change their mode of operation from monostable to bistable, or even to multistable, based on the experimental fine-tuning of readily accessible parameters. In order to successfully design robust, bistable synthetic circuits to be used as biomolecular probes, or understand modes of operation of such naturally occurring circuits, we must identify parameters that are key in determining their characteristics.ResultsHere, we analyze the bistability properties of a general, asymmetric genetic toggle switch based on a chemical-reaction kinetic description. By making appropriate approximations, we are able to reduce the system to two coupled differential equations. Their deterministic stability analysis and stochastic numerical simulations are in excellent agreement. Drawing upon this general framework, we develop a model of an experimentally realized asymmetric bistable genetic switch based on the LacI and TetR repressors. By varying the concentrations of two synthetic inducers, doxycycline and isopropyl β-D-1-thiogalactopyranoside, we predict that it will be possible to repeatedly fine-tune the mode of operation of this genetic switch from monostable to bistable, as well as the switching rates over many orders of magnitude, in an experimental setting. Furthermore, we find that the shape and size of the bistability region is closely connected with plasmid copy number.ConclusionsBased on our numerical calculations of the LacI-TetR asymmetric bistable switch phase diagram, we propose a generic work-flow for developing and applying biomolecular probes: Their initial state of operation should be specified by controlling inducer concentrations, and dilution due to cellular division would turn the probes into memory devices in which information could be preserved over multiple generations. Additionally, insights from our analysis of the LacI-TetR system suggest that this particular system is readily available to be employed in this kind of probe.