No Arabic abstract
We assess the impact of cell cycle noise on gene circuit dynamics. For bistable genetic switches and excitable circuits, we find that transitions between metastable states most likely occur just after cell division and that this concentration effect intensifies in the presence of transcriptional delay. We explain this concentration effect with a 3-states stochastic model. For genetic oscillators, we quantify the temporal correlations between daughter cells induced by cell division. Temporal correlations must be captured properly in order to accurately quantify noise sources within gene networks.
Although accumulation of molecular damage is suggested to be an important molecular mechanism of aging, a quantitative link between the dynamics of damage accumulation and mortality of species has so far remained elusive. To address this question, we examine stability properties of a generic gene regulatory network (GRN) and demonstrate that many characteristics of aging and the associated population mortality rate emerge as inherent properties of the critical dynamics of gene regulation and metabolic levels. Based on the analysis of age-dependent changes in gene-expression and metabolic profiles in Drosophila melanogaster, we explicitly show that the underlying GRNs are nearly critical and inherently unstable. This instability manifests itself as aging in the form of distortion of gene expression and metabolic profiles with age, and causes the characteristic increase in mortality rate with age as described by a form of the Gompertz law. In addition, we explain late-life mortality deceleration observed at very late ages for large populations. We show that aging contains a stochastic component, related to accumulation of regulatory errors in transcription/translation/metabolic pathways due to imperfection of signaling cascades in the network and of responses to environmental factors. We also establish that there is a strong deterministic component, suggesting genetic control. Since mortality in humans, where it is characterized best, is strongly associated with the incidence of age-related diseases, our findings support the idea that aging is the driving force behind the development of chronic human diseases.
In contrast to engineering applications, in which the structure of control laws are designed to satisfy prescribed function requirements, in biology it is often necessary to infer gene-circuit function from incomplete data on gene-circuit structure. By using the feed-forward loop as a model system, this paper introduces a technique for classifying gene-circuit function given gene-circuit structure. In simulations performed on a comprehensive set of models that span a broad range of parameter space, some designs are robust, producing one unique type of functional response regardless of parameter selection. Other designs may exhibit a variety of functional responses, depending upon parameter values. We conclude that, although some feed-forward loop models have designs that lend themselves to unique function inference, others have designs for which the function type may be uncertain.
Recent experiments showed that multiple copies of the molecular machine RNA polymerase (RNAP) can efficiently synthesize mRNA collectively in the active state of the promoter. However, environmentally-induced promoter repression results in long-distance antagonistic interactions that drastically reduce the speed of RNAPs and cause a quick arrest of mRNA synthesis. The mechanism underlying this transition between cooperative and antagonistic dynamics remains poorly understood. In this Letter, we introduce a continuum deterministic model for the translocation of RNAPs, where the speed of an RNAP is coupled to the local DNA supercoiling as well as the density of RNAPs on the gene. We assume that torsional stress experienced by individual RNAPs is exacerbated by high RNAP density on the gene and that transcription factors act as physical barriers to the diffusion of DNA supercoils. We show that this minimal model exhibits two transcription modes mediated by the torsional stress: a fluid mode when the promoter is active and a torsionally stressed mode when the promoter is repressed, in quantitative agreement with experimentally observed dynamics of co-transcribing RNAPs. Our work provides an important step towards understanding the collective dynamics of molecular machines involved in gene expression.
Cell internalization of a blastomere, namely gastrulation, is a common and significant milestone during development of metazoans from worm to human, which generates multiple embryonic layers with distinct cell fates and spatial organizations. Although many molecular activities (e.g., cell polarization, asymmetrical intercellular adhesion, and apical actomyosin cortex contraction) have been revealed to facilitate this morphogenetic process, in this paper, we focus on gastrulation of the worm Caenorhabditis elegans and demonstrate that even a simple mechanical system, like a group of cells with isotropic repulsive and attractive interactions, can experience such internalization behavior spontaneously when dividing within a confined space. In principle, when the total cell number exceeds a threshold, a double-layer structure acquires lower potential energy and longer neighbor distance than the single-layer one. Besides, both mechanical analysis and simulation suggest that the cells with a large size or placed near a small-curvature boundary are easier to internalize. Last but not least, extra regulation on a limited part of cells to internalize autonomously can stabilize this process against motional noise. Our work successfully recaptures many key characteristics in worm gastrulation by mechanical modeling and provides a novel and rational interpretation on how this phenomenon emerges and is optimally programed.
We analyze gene co-expression network under the random matrix theory framework. The nearest neighbor spacing distribution of the adjacency matrix of this network follows Gaussian orthogonal statistics of random matrix theory (RMT). Spectral rigidity test follows random matrix prediction for a certain range, and deviates after wards. Eigenvector analysis of the network using inverse participation ratio (IPR) suggests that the statistics of bulk of the eigenvalues of network is consistent with those of the real symmetric random matrix, whereas few eigenvalues are localized. Based on these IPR calculations, we can divide eigenvalues in three sets; (A) The non-degenerate part that follows RMT. (B) The non-degenerate part, at both ends and at intermediate eigenvalues, which deviate from RMT and expected to contain information about {it important nodes} in the network. (C) The degenerate part with $zero$ eigenvalue, which fluctuates around RMT predicted value. We identify nodes corresponding to the dominant modes of the corresponding eigenvectors and analyze their structural properties.