1
|
Zhao Y, Liu L, Hassett R, Siepel A. Model-based characterization of the equilibrium dynamics of transcription initiation and promoter-proximal pausing in human cells. Nucleic Acids Res 2023; 51:e106. [PMID: 37889042 PMCID: PMC10681744 DOI: 10.1093/nar/gkad843] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2023] [Revised: 09/13/2023] [Accepted: 09/21/2023] [Indexed: 10/28/2023] Open
Abstract
In metazoans, both transcription initiation and the escape of RNA polymerase (RNAP) from promoter-proximal pausing are key rate-limiting steps in gene expression. These processes play out at physically proximal sites on the DNA template and appear to influence one another through steric interactions. Here, we examine the dynamics of these processes using a combination of statistical modeling, simulation, and analysis of real nascent RNA sequencing data. We develop a simple probabilistic model that jointly describes the kinetics of transcription initiation, pause-escape, and elongation, and the generation of nascent RNA sequencing read counts under steady-state conditions. We then extend this initial model to allow for variability across cells in promoter-proximal pause site locations and steric hindrance of transcription initiation from paused RNAPs. In an extensive series of simulations, we show that this model enables accurate estimation of initiation and pause-escape rates. Furthermore, we show by simulation and analysis of real data that pause-escape is often strongly rate-limiting and that steric hindrance can dramatically reduce initiation rates. Our modeling framework is applicable to a variety of inference problems, and our software for estimation and simulation is freely available.
Collapse
Affiliation(s)
- Yixin Zhao
- Simons Center for Quantitative Biology, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA
| | - Lingjie Liu
- Simons Center for Quantitative Biology, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA
- Graduate Program in Genetics, Stony Brook University, Stony Brook, NY, USA
| | - Rebecca Hassett
- Simons Center for Quantitative Biology, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA
| | - Adam Siepel
- Simons Center for Quantitative Biology, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA
- Graduate Program in Genetics, Stony Brook University, Stony Brook, NY, USA
| |
Collapse
|
2
|
Keller SH, Deng H, Lim B. Regulation of the dynamic RNA Pol II elongation rate in Drosophila embryos. Cell Rep 2023; 42:113225. [PMID: 37837623 PMCID: PMC10842316 DOI: 10.1016/j.celrep.2023.113225] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2022] [Revised: 08/20/2023] [Accepted: 09/22/2023] [Indexed: 10/16/2023] Open
Abstract
An increasing number of studies have shown the key role that RNA polymerase II (RNA Pol II) elongation plays in gene regulation. We systematically examine how various enhancers, promoters, and gene body composition influence the RNA Pol II elongation rate through a single-cell-resolution live imaging assay. By using reporter constructs containing 5' MS2 and 3' PP7 repeating stem loops, we quantify the rate of RNA Pol II elongation in live Drosophila embryos. We find that promoters and exonic gene lengths have no effect on elongation rate, while enhancers and the presence of long introns may significantly change how quickly RNA Pol II moves across a gene. Furthermore, we observe in multiple constructs that the RNA Pol II elongation rate accelerates after the transcriptional onset of nuclear cycle 14 in Drosophila embryos. Our study provides a single-cell view of various mechanisms that affect the dynamic RNA Pol II elongation rate, ultimately affecting the rate of mRNA production.
Collapse
Affiliation(s)
- Samuel H Keller
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, PA 19104, USA
| | - Hao Deng
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, PA 19104, USA
| | - Bomyi Lim
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, PA 19104, USA.
| |
Collapse
|
3
|
Shi C, Yang X, Zhang J, Zhou T. Stochastic modeling of the mRNA life process: A generalized master equation. Biophys J 2023; 122:4023-4041. [PMID: 37653725 PMCID: PMC10598292 DOI: 10.1016/j.bpj.2023.08.024] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2023] [Revised: 06/29/2023] [Accepted: 08/29/2023] [Indexed: 09/02/2023] Open
Abstract
The mRNA life cycle is a complex biochemical process, involving transcription initiation, elongation, termination, splicing, and degradation. Each of these molecular events is multistep and can create a memory. The effect of this molecular memory on gene expression is not clear, although there are many related yet scattered experimental reports. To address this important issue, we develop a general theoretical framework formulated as a master equation in the sense of queue theory, which can reduce to multiple previously studied gene models in limiting cases. This framework allows us to interpret experimental observations, extract kinetic parameters from experimental data, and identify how the mRNA kinetics vary under regulatory influences. Notably, it allows us to evaluate the influences of elongation processes on mature RNA distribution; e.g., we find that the non-exponential elongation time can induce the bimodal mRNA expression and there is an optimal elongation noise intensity such that the mature RNA noise achieves the lowest level. In a word, our framework can not only provide insight into complex mRNA life processes but also bridge a dialogue between theoretical studies and experimental data.
Collapse
Affiliation(s)
- Changhong Shi
- State Key Laboratory of Respiratory Disease, School of Public Health, Guangzhou Medical University, Guangzhou, China
| | - Xiyan Yang
- School of Financial Mathematics and Statistics, Guangdong University of Finance, Guangzhou, China
| | - Jiajun Zhang
- School of Mathematics and Computational Science and Guangdong Province Key Laboratory of Computational Science, Sun Yat-Sen University, Guangzhou, China.
| | - Tianshou Zhou
- School of Mathematics and Computational Science and Guangdong Province Key Laboratory of Computational Science, Sun Yat-Sen University, Guangzhou, China.
| |
Collapse
|
4
|
Zhang J, Han X, Ma L, Xu S, Lin Y. Deciphering a global source of non-genetic heterogeneity in cancer cells. Nucleic Acids Res 2023; 51:9019-9038. [PMID: 37587722 PMCID: PMC10516630 DOI: 10.1093/nar/gkad666] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2022] [Revised: 07/09/2023] [Accepted: 08/09/2023] [Indexed: 08/18/2023] Open
Abstract
Cell-to-cell variability within a clonal population, also known as non-genetic heterogeneity, has created significant challenges for intervening with diseases such as cancer. While non-genetic heterogeneity can arise from the variability in the expression of specific genes, it remains largely unclear whether and how clonal cells could be heterogeneous in the expression of the entire transcriptome. Here, we showed that gene transcriptional activity is globally modulated in individual cancer cells, leading to non-genetic heterogeneity in the global transcription rate. Such heterogeneity contributes to cell-to-cell variability in transcriptome size and displays both dynamic and static characteristics, with the global transcription rate temporally modulated in a cell-cycle-coupled manner and the time-averaged rate being distinct between cells and heritable across generations. Additional evidence indicated the role of ATP metabolism in this heterogeneity, and suggested its implication in intrinsic cancer drug tolerance. Collectively, our work shed light on the mode, mechanism, and implication of a global but often hidden source of non-genetic heterogeneity.
Collapse
Affiliation(s)
- Jianhan Zhang
- Center for Quantitative Biology and Peking-Tsinghua Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
- The MOE Key Laboratory of Cell Proliferation and Differentiation, School of Life Sciences, Peking University, Beijing 100871, China
| | - Xu Han
- Center for Quantitative Biology and Peking-Tsinghua Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
- The MOE Key Laboratory of Cell Proliferation and Differentiation, School of Life Sciences, Peking University, Beijing 100871, China
| | - Liang Ma
- Center for Quantitative Biology and Peking-Tsinghua Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
- The MOE Key Laboratory of Cell Proliferation and Differentiation, School of Life Sciences, Peking University, Beijing 100871, China
| | - Shuhui Xu
- Center for Quantitative Biology and Peking-Tsinghua Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
- The MOE Key Laboratory of Cell Proliferation and Differentiation, School of Life Sciences, Peking University, Beijing 100871, China
| | - Yihan Lin
- Center for Quantitative Biology and Peking-Tsinghua Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
- The MOE Key Laboratory of Cell Proliferation and Differentiation, School of Life Sciences, Peking University, Beijing 100871, China
| |
Collapse
|
5
|
Szavits-Nossan J, Grima R. Uncovering the effect of RNA polymerase steric interactions on gene expression noise: Analytical distributions of nascent and mature RNA numbers. Phys Rev E 2023; 108:034405. [PMID: 37849194 DOI: 10.1103/physreve.108.034405] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2023] [Accepted: 08/24/2023] [Indexed: 10/19/2023]
Abstract
The telegraph model is the standard model of stochastic gene expression, which can be solved exactly to obtain the distribution of mature RNA numbers per cell. A modification of this model also leads to an analytical distribution of nascent RNA numbers. These solutions are routinely used for the analysis of single-cell data, including the inference of transcriptional parameters. However, these models neglect important mechanistic features of transcription elongation, such as the stochastic movement of RNA polymerases and their steric (excluded-volume) interactions. Here we construct a model of gene expression describing promoter switching between inactive and active states, binding of RNA polymerases in the active state, their stochastic movement including steric interactions along the gene, and their unbinding leading to a mature transcript that subsequently decays. We derive the steady-state distributions of the nascent and mature RNA numbers in two important limiting cases: constitutive expression and slow promoter switching. We show that RNA fluctuations are suppressed by steric interactions between RNA polymerases, and that this suppression can in some instances even lead to sub-Poissonian fluctuations; these effects are most pronounced for nascent RNA and less prominent for mature RNA, since the latter is not a direct sensor of transcription. We find a relationship between the parameters of our microscopic mechanistic model and those of the standard models that ensures excellent consistency in their prediction of the first and second RNA number moments over vast regions of parameter space, encompassing slow, intermediate, and rapid promoter switching, provided the RNA number distributions are Poissonian or super-Poissonian. Furthermore, we identify the limitations of inference from mature RNA data, specifically showing that it cannot differentiate between highly distinct RNA polymerase traffic patterns on a gene.
Collapse
Affiliation(s)
- Juraj Szavits-Nossan
- School of Biological Sciences, University of Edinburgh, Edinburgh EH9 3JH, United Kingdom
| | - Ramon Grima
- School of Biological Sciences, University of Edinburgh, Edinburgh EH9 3JH, United Kingdom
| |
Collapse
|
6
|
Weidemann DE, Holehouse J, Singh A, Grima R, Hauf S. The minimal intrinsic stochasticity of constitutively expressed eukaryotic genes is sub-Poissonian. SCIENCE ADVANCES 2023; 9:eadh5138. [PMID: 37556551 PMCID: PMC10411910 DOI: 10.1126/sciadv.adh5138] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/08/2023] [Accepted: 07/07/2023] [Indexed: 08/11/2023]
Abstract
Gene expression inherently gives rise to stochastic variation ("noise") in the production of gene products. Minimizing noise is crucial for ensuring reliable cellular functions. However, noise cannot be suppressed below a certain intrinsic limit. For constitutively expressed genes, this limit is typically assumed to be Poissonian noise, wherein the variance in mRNA numbers is equal to their mean. Here, we demonstrate that several cell division genes in fission yeast exhibit mRNA variances significantly below this limit. The reduced variance can be explained by a gene expression model incorporating multiple transcription and mRNA degradation steps. Notably, in this sub-Poissonian regime, distinct from Poissonian or super-Poissonian regimes, cytoplasmic noise is effectively suppressed through a higher mRNA export rate. Our findings redefine the lower limit of eukaryotic gene expression noise and uncover molecular requirements for achieving ultralow noise, which is expected to be important for vital cellular functions.
Collapse
Affiliation(s)
- Douglas E. Weidemann
- Department of Biological Sciences, Virginia Tech, Blacksburg, VA 24061, USA
- Fralin Life Sciences Institute, Virginia Tech, Blacksburg, VA 24061, USA
| | - James Holehouse
- The Santa Fe Institute, 1399 Hyde Park Road, Santa Fe, NM 87510, USA
| | - Abhyudai Singh
- Department of Electrical and Computer Engineering, University of Delaware, Newark, DE 19716, USA
- Department of Biomedical Engineering, University of Delaware, Newark, DE 19716, USA
| | - Ramon Grima
- School of Biological Sciences, University of Edinburgh, Edinburgh EH9 3JR, UK
| | - Silke Hauf
- Department of Biological Sciences, Virginia Tech, Blacksburg, VA 24061, USA
- Fralin Life Sciences Institute, Virginia Tech, Blacksburg, VA 24061, USA
| |
Collapse
|
7
|
Sevier SA, Hormoz S. Collective polymerase dynamics emerge from DNA supercoiling during transcription. Biophys J 2022; 121:4153-4165. [PMID: 36171726 PMCID: PMC9675029 DOI: 10.1016/j.bpj.2022.09.026] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2022] [Revised: 07/19/2022] [Accepted: 09/22/2022] [Indexed: 11/29/2022] Open
Abstract
All biological processes ultimately come from physical interactions. The mechanical properties of DNA play a critical role in transcription. RNA polymerase can over or under twist DNA (referred to as DNA supercoiling) when it moves along a gene, resulting in mechanical stresses in DNA that impact its own motion and that of other polymerases. For example, when enough supercoiling accumulates, an isolated polymerase halts, and transcription stops. DNA supercoiling can also mediate nonlocal interactions between polymerases that shape gene expression fluctuations. Here, we construct a comprehensive model of transcription that captures how RNA polymerase motion changes the degree of DNA supercoiling, which in turn feeds back into the rate at which polymerases are recruited and move along the DNA. Surprisingly, our model predicts that a group of three or more polymerases move together at a constant velocity and sustain their motion (forming what we call a polymeton), whereas one or two polymerases would have halted. We further show that accounting for the impact of DNA supercoiling on both RNA polymerase recruitment and velocity recapitulates empirical observations of gene expression fluctuations. Finally, we propose a mechanical toggle switch whereby interactions between genes are mediated by DNA twisting as opposed to proteins. Understanding the mechanical regulation of gene expression provides new insights into how endogenous genes can interact and informs the design of new forms of engineered interactions.
Collapse
Affiliation(s)
- Stuart A Sevier
- Department of Systems Biology, Harvard Medical School, Boston, Massachusetts; Department of Data Sciences, Dana-Farber Cancer Institute, Boston, Massachusetts
| | - Sahand Hormoz
- Department of Systems Biology, Harvard Medical School, Boston, Massachusetts; Department of Data Sciences, Dana-Farber Cancer Institute, Boston, Massachusetts; Broad Institute of MIT and Harvard, Cambridge, Massachusetts.
| |
Collapse
|
8
|
Szavits-Nossan J, Grima R. Mean-field theory accurately captures the variation of copy number distributions across the mRNA life cycle. Phys Rev E 2022; 105:014410. [PMID: 35193216 DOI: 10.1103/physreve.105.014410] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2021] [Accepted: 11/15/2021] [Indexed: 06/14/2023]
Abstract
We consider a stochastic model where a gene switches between two states, an mRNA transcript is released in the active state, and subsequently it undergoes an arbitrary number of sequential unimolecular steps before being degraded. The reactions effectively describe various stages of the mRNA life cycle such as initiation, elongation, termination, splicing, export, and degradation. We construct a mean-field approach that leads to closed-form steady-state distributions for the number of transcript molecules at each stage of the mRNA life cycle. By comparison with stochastic simulations, we show that the approximation is highly accurate over all the parameter space, independent of the type of expression (constitutive or bursty) and of the shape of the distribution (unimodal, bimodal, and nearly bimodal). The theory predicts that in a population of identical cells, any bimodality is gradually washed away as the mRNA progresses through its life cycle.
Collapse
Affiliation(s)
- Juraj Szavits-Nossan
- School of Biological Sciences, University of Edinburgh, Edinburgh EH9 3JH, United Kingdom
| | - Ramon Grima
- School of Biological Sciences, University of Edinburgh, Edinburgh EH9 3JH, United Kingdom
| |
Collapse
|
9
|
Koopmans L, Youk H. Predictive landscapes hidden beneath biological cellular automata. J Biol Phys 2021; 47:355-369. [PMID: 34739687 PMCID: PMC8603977 DOI: 10.1007/s10867-021-09592-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2021] [Accepted: 10/14/2021] [Indexed: 11/11/2022] Open
Abstract
To celebrate Hans Frauenfelder's achievements, we examine energy(-like) "landscapes" for complex living systems. Energy landscapes summarize all possible dynamics of some physical systems. Energy(-like) landscapes can explain some biomolecular processes, including gene expression and, as Frauenfelder showed, protein folding. But energy-like landscapes and existing frameworks like statistical mechanics seem impractical for describing many living systems. Difficulties stem from living systems being high dimensional, nonlinear, and governed by many, tightly coupled constituents that are noisy. The predominant modeling approach is devising differential equations that are tailored to each living system. This ad hoc approach faces the notorious "parameter problem": models have numerous nonlinear, mathematical functions with unknown parameter values, even for describing just a few intracellular processes. One cannot measure many intracellular parameters or can only measure them as snapshots in time. Another modeling approach uses cellular automata to represent living systems as discrete dynamical systems with binary variables. Quantitative (Hamiltonian-based) rules can dictate cellular automata (e.g., Cellular Potts Model). But numerous biological features, in current practice, are qualitatively described rather than quantitatively (e.g., gene is (highly) expressed or not (highly) expressed). Cellular automata governed by verbal rules are useful representations for living systems and can mitigate the parameter problem. However, they can yield complex dynamics that are difficult to understand because the automata-governing rules are not quantitative and much of the existing mathematical tools and theorems apply to continuous but not discrete dynamical systems. Recent studies found ways to overcome this challenge. These studies either discovered or suggest an existence of predictive "landscapes" whose shapes are described by Lyapunov functions and yield "equations of motion" for a "pseudo-particle." The pseudo-particle represents the entire cellular lattice and moves on the landscape, thereby giving a low-dimensional representation of the cellular automata dynamics. We outline this promising modeling strategy.
Collapse
Affiliation(s)
- Lars Koopmans
- Program in Applied Physics, Delft University of Technology, Delft, The Netherlands
- Department of Systems Biology, University of Massachusetts Chan Medical School, Worcester, MA, USA
| | - Hyun Youk
- Department of Systems Biology, University of Massachusetts Chan Medical School, Worcester, MA, USA.
| |
Collapse
|
10
|
Liu J, Hansen D, Eck E, Kim YJ, Turner M, Alamos S, Garcia HG. Real-time single-cell characterization of the eukaryotic transcription cycle reveals correlations between RNA initiation, elongation, and cleavage. PLoS Comput Biol 2021; 17:e1008999. [PMID: 34003867 PMCID: PMC8162642 DOI: 10.1371/journal.pcbi.1008999] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2020] [Revised: 05/28/2021] [Accepted: 04/23/2021] [Indexed: 12/23/2022] Open
Abstract
The eukaryotic transcription cycle consists of three main steps: initiation, elongation, and cleavage of the nascent RNA transcript. Although each of these steps can be regulated as well as coupled with each other, their in vivo dissection has remained challenging because available experimental readouts lack sufficient spatiotemporal resolution to separate the contributions from each of these steps. Here, we describe a novel application of Bayesian inference techniques to simultaneously infer the effective parameters of the transcription cycle in real time and at the single-cell level using a two-color MS2/PP7 reporter gene and the developing fruit fly embryo as a case study. Our method enables detailed investigations into cell-to-cell variability in transcription-cycle parameters as well as single-cell correlations between these parameters. These measurements, combined with theoretical modeling, suggest a substantial variability in the elongation rate of individual RNA polymerase molecules. We further illustrate the power of this technique by uncovering a novel mechanistic connection between RNA polymerase density and nascent RNA cleavage efficiency. Thus, our approach makes it possible to shed light on the regulatory mechanisms in play during each step of the transcription cycle in individual, living cells at high spatiotemporal resolution.
Collapse
Affiliation(s)
- Jonathan Liu
- Department of Physics, University of California at Berkeley, Berkeley, California, United States of America
| | - Donald Hansen
- Institute of Pharmacy and Molecular Biotechnology, University of Heidelberg, Heidelberg, Germany
| | - Elizabeth Eck
- Biophysics Graduate Group, University of California at Berkeley, Berkeley, California, United States of America
| | - Yang Joon Kim
- Biophysics Graduate Group, University of California at Berkeley, Berkeley, California, United States of America
| | - Meghan Turner
- Biophysics Graduate Group, University of California at Berkeley, Berkeley, California, United States of America
| | - Simon Alamos
- Department of Plant and Microbial Biology, University of California at Berkeley, Berkeley, California, United States of America
| | - Hernan G. Garcia
- Department of Physics, University of California at Berkeley, Berkeley, California, United States of America
- Biophysics Graduate Group, University of California at Berkeley, Berkeley, California, United States of America
- Department of Molecular and Cell Biology, University of California at Berkeley, Berkeley, California, United States of America
- Institute for Quantitative Biosciences-QB3, University of California at Berkeley, Berkeley, California, United States of America
| |
Collapse
|
11
|
Filatova T, Popovic N, Grima R. Statistics of Nascent and Mature RNA Fluctuations in a Stochastic Model of Transcriptional Initiation, Elongation, Pausing, and Termination. Bull Math Biol 2020; 83:3. [PMID: 33351158 PMCID: PMC7755674 DOI: 10.1007/s11538-020-00827-7] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2020] [Accepted: 10/26/2020] [Indexed: 12/20/2022]
Abstract
Recent advances in fluorescence microscopy have made it possible to measure the fluctuations of nascent (actively transcribed) RNA. These closely reflect transcription kinetics, as opposed to conventional measurements of mature (cellular) RNA, whose kinetics is affected by additional processes downstream of transcription. Here, we formulate a stochastic model which describes promoter switching, initiation, elongation, premature detachment, pausing, and termination while being analytically tractable. We derive exact closed-form expressions for the mean and variance of nascent RNA fluctuations on gene segments, as well as of total nascent RNA on a gene. We also obtain exact expressions for the first two moments of mature RNA fluctuations and approximate distributions for total numbers of nascent and mature RNA. Our results, which are verified by stochastic simulation, uncover the explicit dependence of the statistics of both types of RNA on transcriptional parameters and potentially provide a means to estimate parameter values from experimental data.
Collapse
Affiliation(s)
- Tatiana Filatova
- School of Biological Sciences, The University of Edinburgh, Edinburgh, UK.,School of Mathematics and Maxwell Institute for Mathematical Sciences, The University of Edinburgh, Edinburgh, UK
| | - Nikola Popovic
- School of Mathematics and Maxwell Institute for Mathematical Sciences, The University of Edinburgh, Edinburgh, UK
| | - Ramon Grima
- School of Biological Sciences, The University of Edinburgh, Edinburgh, UK.
| |
Collapse
|
12
|
Lammers NC, Kim YJ, Zhao J, Garcia HG. A matter of time: Using dynamics and theory to uncover mechanisms of transcriptional bursting. Curr Opin Cell Biol 2020; 67:147-157. [PMID: 33242838 PMCID: PMC8498946 DOI: 10.1016/j.ceb.2020.08.001] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2020] [Accepted: 08/03/2020] [Indexed: 12/18/2022]
Abstract
Eukaryotic transcription generally occurs in bursts of activity lasting minutes to hours; however, state-of-the-art measurements have revealed that many of the molecular processes that underlie bursting, such as transcription factor binding to DNA, unfold on timescales of seconds. This temporal disconnect lies at the heart of a broader challenge in physical biology of predicting transcriptional outcomes and cellular decision-making from the dynamics of underlying molecular processes. Here, we review how new dynamical information about the processes underlying transcriptional control can be combined with theoretical models that predict not only averaged transcriptional dynamics, but also their variability, to formulate testable hypotheses about the molecular mechanisms underlying transcriptional bursting and control.
Collapse
Affiliation(s)
- Nicholas C Lammers
- Biophysics Graduate Group, University of California at Berkeley, Berkeley, CA, USA
| | - Yang Joon Kim
- Biophysics Graduate Group, University of California at Berkeley, Berkeley, CA, USA
| | - Jiaxi Zhao
- Department of Physics, University of California at Berkeley, Berkeley, CA, USA
| | - Hernan G Garcia
- Biophysics Graduate Group, University of California at Berkeley, Berkeley, CA, USA; Department of Physics, University of California at Berkeley, Berkeley, CA, USA; Department of Molecular and Cell Biology, University of California at Berkeley, Berkeley, CA, USA; Institute for Quantitative Biosciences-QB3, University of California at Berkeley, Berkeley, CA, USA.
| |
Collapse
|