Variability in mRNA translation: a random matrix theory approach

The rate of mRNA translation depends on the initiation, elongation, and termination rates of ribosomes along the mRNA. These rates depend on many “local” factors like the abundance of free ribosomes and tRNA molecules in the vicinity of the mRNA molecule. All these factors are stochastic and their experimental measurements are also noisy. An important question is how protein production in the cell is affected by this considerable variability. We develop a new theoretical framework for addressing this question by modeling the rates as identically and independently distributed random variables and using tools from random matrix theory to analyze the steady-state production rate. The analysis reveals a principle of universality: the average protein production rate depends only on the of the set of possible values that the random variable may attain. This explains how total protein production can be stabilized despite the overwhelming stochasticticity underlying cellular processes.

[1]  Eytan Domany,et al.  An exact solution of a one-dimensional asymmetric exclusion model with open boundaries , 1992 .

[2]  Michael Margaliot,et al.  A model for competition for ribosomes in the cell , 2015, Journal of The Royal Society Interface.

[3]  Michael Margaliot,et al.  Ribosome flow model with positive feedback , 2013, Journal of The Royal Society Interface.

[4]  Michal Ziv-Ukelson,et al.  Composite effects of gene determinants on the translation speed and density of ribosomes , 2011, Genome Biology.

[5]  Eran Gazit,et al.  The RNA Polymerase Flow Model of Gene Transcription , 2014, IEEE Transactions on Biomedical Circuits and Systems.

[6]  L. S. Churchman,et al.  Genome-wide profiling of RNA polymerase transcription at nucleotide resolution in human cells with native elongating transcript sequencing , 2016, Nature Protocols.

[7]  Tamir Tuller,et al.  Efficient Manipulations of Synonymous Mutations for Controlling Translation Rate: An Analytical Approach , 2012, J. Comput. Biol..

[8]  Nathan Morris,et al.  Codon Optimality Is a Major Determinant of mRNA Stability , 2015, Cell.

[9]  J. Russell,et al.  Energetics of bacterial growth: balance of anabolic and catabolic reactions , 1995 .

[10]  P. Freire-Pritchett,et al.  Ribosome collisions trigger cis-acting feedback inhibition of translation initiation , 2020, eLife.

[11]  R. Hegde,et al.  The ASC-1 Complex Disassembles Collided Ribosomes , 2020, Molecular cell.

[12]  Wen-Chyuan Yueh.  EIGENVALUES OF SEVERAL TRIDIAGONAL MATRICES , 2005 .

[13]  Michael Margaliot,et al.  A deterministic mathematical model for bidirectional excluded flow with Langmuir kinetics , 2016, PloS one.

[14]  J. Magnus.  On Differentiating Eigenvalues and Eigenvectors , 1985, Econometric Theory.

[15]  C. Kimchi-Sarfaty,et al.  Understanding the contribution of synonymous mutations to human disease , 2011, Nature Reviews Genetics.

[16]  Isaac Meilijson,et al.  Genome-Scale Analysis of Translation Elongation with a Ribosome Flow Model , 2011, PLoS Comput. Biol..

[17]  T. Lambert.  Antibiotics that affect the ribosome. , 2012, Revue scientifique et technique.

[18]  A. Pipkin,et al.  Kinetics of biopolymerization on nucleic acid templates , 1968, Biopolymers.

[19]  F. Spitzer.  Interaction of Markov processes , 1970 .

[20]  Hamid Teimouri,et al.  Dynamics of translation can determine the spatial organization of membrane-bound proteins and their mRNA , 2017, Proceedings of the National Academy of Sciences.

[21]  Tamir Tuller,et al.  New computational model for miRNA-mediated repression reveals novel regulatory roles of miRNA bindings inside the coding region , 2020, Bioinform..

[22]  M. Margaliot,et al.  A deterministic model for one-dimensional excluded flow with local interactions , 2017, PloS one.

[23]  Michael Margaliot,et al.  Ribosome Flow Model with Different Site Sizes , 2019, SIAM J. Appl. Dyn. Syst..

[24]  R. Takors,et al.  Protein production in Escherichia coli is guided by the trade-off between intracellular substrate availability and energy cost , 2019, Microbial Cell Factories.

[25]  T. Tuller,et al.  Properties and determinants of codon decoding time distributions , 2014, BMC Genomics.

[26]  Mads Kærn,et al.  Noise in eukaryotic gene expression , 2003, Nature.

[27]  P. Tomançak,et al.  Global Analysis of mRNA Localization Reveals a Prominent Role in Organizing Cellular Architecture and Function , 2007, Cell.

[28]  J. H. Gibbs,et al.  Concerning the kinetics of polypeptide synthesis on polyribosomes , 1969 .

[29]  M. Kunitski,et al.  Double-slit photoelectron interference in strong-field ionization of the neon dimer , 2018, Nature Communications.

[30]  Tamir Tuller,et al.  The effect of tRNA levels on decoding times of mRNA codons , 2014, Nucleic acids research.

[31]  F. Buttgereit,et al.  A hierarchy of ATP-consuming processes in mammalian cells. , 1995, The Biochemical journal.

[32]  Daniel N. Wilson.  Ribosome-targeting antibiotics and mechanisms of bacterial resistance , 2013, Nature Reviews Microbiology.

[33]  J. Davies,et al.  Molecular Biology of the Cell , 1983, Bristol Medico-Chirurgical Journal.

[34]  B. Derrida,et al.  Exact solution of a 1d asymmetric exclusion model using a matrix formulation , 1993 .

[35]  R. Green,et al.  Ribosomopathies: There’s strength in numbers , 2017, Science.

[36]  Michael Margaliot,et al.  Analyzing Linear Communication Networks Using the Ribosome Flow Model , 2015, 2015 IEEE International Conference on Computer and Information Technology; Ubiquitous Computing and Communications; Dependable, Autonomic and Secure Computing; Pervasive Intelligence and Computing.

[37]  Michael Margaliot,et al.  Stability Analysis of the Ribosome Flow Model , 2012, IEEE/ACM Transactions on Computational Biology and Bioinformatics.

[38]  Y. Pilpel,et al.  An Evolutionarily Conserved Mechanism for Controlling the Efficiency of Protein Translation , 2010, Cell.

[39]  M. Margaliot,et al.  Modeling and Analyzing the Flow of Molecular Machines in Gene Expression , 2018 .

[40]  Fabian J Theis,et al.  Computational analysis of cell-to-cell heterogeneity in single-cell RNA-sequencing data reveals hidden subpopulations of cells , 2015, Nature Biotechnology.

[41]  Magnus Johansson,et al.  Sequence-dependent elongation dynamics on macrolide-bound ribosomes. , 2014, Cell reports.

[42]  U. Alon,et al.  Central dogma rates and the trade-off between precision and economy in gene expression , 2019, Nature Communications.

[43]  M. Margaliot,et al.  Sensitivity of mRNA Translation , 2014, Scientific Reports.

[44]  E. Nudler,et al.  RNA polymerase and the ribosome: the close relationship. , 2013, Current opinion in microbiology.

[45]  F. Dyson.  The Dynamics of a Disordered Linear Chain , 1953 .

[46]  Rati Sharma.  Extrinsic noise acts to lower protein production at higher translation initiation rates , 2020, bioRxiv.

[47]  Tamir Tuller,et al.  Estimation of ribosome profiling performance and reproducibility at various levels of resolution , 2016, Biology Direct.

[48]  P. Mahadevan,et al.  An overview , 2007, Journal of Biosciences.

[49]  Thomas E Gorochowski,et al.  A Minimal Model of Ribosome Allocation Dynamics Captures Trade-offs in Expression between Endogenous and Synthetic Genes. , 2016, ACS synthetic biology.

[50]  M. Margaliot,et al.  Optimal Down Regulation of mRNA Translation , 2016, Scientific Reports.

[51]  Michael Margaliot,et al.  On Minimizing the Maximal Characteristic Frequency of a Linear Chain , 2017, IEEE Transactions on Automatic Control.

[52]  Elfed Lewis,et al.  FPGA based Real time 'secure' body temperature monitoring suitable for WBSN 2015 IEEE International Conference on Computer and Information Technology; Ubiquitous Computing and Communications; Dependable, Autonomic and Secure Computing; Pervasive Intelligence and Computing , 2015 .

[53]  T. Tuller,et al.  Evidence of translation efficiency adaptation of the coding regions of the bacteriophage lambda , 2017, DNA research : an international journal for rapid publication of reports on genes and genomes.

[54]  Michael Margaliot,et al.  Entrainment to Periodic Initiation and Transition Rates in a Computational Model for Gene Translation , 2014, PloS one.

[55]  B. Schmittmann,et al.  Modeling Translation in Protein Synthesis with TASEP: A Tutorial and Recent Developments , 2011, 1108.3312.

[56]  B. M. Fulk.  MATH , 1992 .

[57]  V. Imbert,et al.  Structure–function insights reveal the human ribosome as a cancer target for antibiotics , 2016, Nature Communications.

[58]  S. Janicki,et al.  Seeing Is Believing: Visualizing Transcriptional Dynamics in Single Cells , 2014, Journal of cellular physiology.

[59]  J. W. Silverstein,et al.  Spectral Analysis of Large Dimensional Random Matrices , 2009 .

[60]  Michael Margaliot,et al.  Networks of ribosome flow models for modeling and analyzing intracellular traffic , 2018, Scientific Reports.

[61]  Vadim N. Gladyshev,et al.  Ribonuclease selection for ribosome profiling , 2016, Nucleic acids research.

[62]  Tamir Tuller,et al.  Predictive biophysical modeling and understanding of the dynamics of mRNA translation and its evolution , 2016, Nucleic acids research.

[63]  Michael Margaliot,et al.  Optimal Translation Along a Circular mRNA , 2017, Scientific Reports.

[64]  K. J. Mukherjee,et al.  Genome engineering for improved recombinant protein expression in Escherichia coli , 2014, Microbial Cell Factories.

[65]  Tamir Tuller,et al.  Quantifying the Effect of Ribosomal Density on mRNA Stability , 2014, PloS one.

[66]  Michael Margaliot,et al.  Ribosome Flow Model on a Ring , 2015, IEEE/ACM Transactions on Computational Biology and Bioinformatics.

[67]  A. Arkin,et al.  Stochastic mechanisms in gene expression. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[68]  Michael Margaliot,et al.  On the Ribosomal Density that Maximizes Protein Translation Rate , 2016, PloS one.

[69]  Michael Margaliot,et al.  Checkable Conditions for Contraction After Small Transients in Time and Amplitude , 2017 .

[70]  M. Barna,et al.  The Discovery of Ribosome Heterogeneity and Its Implications for Gene Regulation and Organismal Life. , 2018, Molecular cell.

[71]  T. Tuller,et al.  Novel insights into gene expression regulation during meiosis revealed by translation elongation dynamics , 2019, npj Systems Biology and Applications.

[72]  Edward P O'Brien,et al.  Determinants of translation speed are randomly distributed across transcripts resulting in a universal scaling of protein synthesis times. , 2018, Physical review. E.

[73]  M. Margaliot,et al.  Maximizing protein translation rate in the non-homogeneous ribosome flow model: a convex optimization approach , 2014, Journal of The Royal Society Interface.

[74]  Nicholas T. Ingolia,et al.  Seeing translation , 2016, Science.

[75]  J. Russell,et al.  Energetics of bacterial growth: balance of anabolic and catabolic reactions. , 1995, Microbiological reviews.

[76]  T. von der Haar,et al.  Mathematical and Computational Modelling of Ribosomal Movement and Protein Synthesis: an overview , 2012, Computational and structural biotechnology journal.

[77]  J. Paulsson.  Summing up the noise in gene networks , 2004, Nature.

[78]  Stijn Sonneveld,et al.  Heterogeneity in mRNA Translation. , 2020, Trends in cell biology.

[79]  Tamir Tuller,et al.  Computational analysis of the oscillatory behavior at the translation level induced by mRNA levels oscillations due to finite intracellular resources , 2018, PLoS Comput. Biol..

[80]  Eduardo D. Sontag,et al.  Mediating Ribosomal Competition by Splitting Pools , 2020, IEEE Control Systems Letters.

[81]  J. Derisi,et al.  Single-cell proteomic analysis of S. cerevisiae reveals the architecture of biological noise , 2006, Nature.

[82]  Nicholas T. Ingolia.  Ribosome profiling: new views of translation, from single codons to genome scale , 2014, Nature Reviews Genetics.

[83]  K. Martin,et al.  mRNA Localization: Gene Expression in the Spatial Dimension , 2009, Cell.

[84]  T. Tuller,et al.  Multiple roles of the coding sequence 5′ end in gene expression regulation , 2014, Nucleic acids research.

[85]  Tobias von der Haar,et al.  Mathematical and Computational Modelling of Ribosomal Movement and Protein Synthesis: an overview. , 2012 .

[86]  C. Fonseca,et al.  Eigenpairs of a family of tridiagonal matrices: three decades later , 2020, Acta Mathematica Hungarica.

[87]  J. Graham,et al.  Architecture and function , 1993 .

[88]  Florence Besse,et al.  Translational control of localized mRNAs: restricting protein synthesis in space and time , 2008, Nature Reviews Molecular Cell Biology.

[89]  Debashish Chowdhury,et al.  Stochastic Transport in Complex Systems: From Molecules to Vehicles , 2010 .

[90]  A. Giraldez,et al.  Ribosome Profiling Shows That miR-430 Reduces Translation Before Causing mRNA Decay in Zebrafish , 2012, Science.

[91]  R. Hendler,et al.  A Model for Protein Synthesis , 1962, Nature.

[92]  Kelvin H. Lee,et al.  Totally asymmetric exclusion process with extended objects: a model for protein synthesis. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[93]  Anna Feldman,et al.  The extent of ribosome queuing in budding yeast , 2018, PLoS Comput. Biol..

[94]  Samie R Jaffrey,et al.  Reading, writing and erasing mRNA methylation , 2019, Nature Reviews Molecular Cell Biology.

[95]  T. Tuller,et al.  Widespread non-modular overlapping codes in the coding regions , 2020, Physical biology.

[96]  M. Margaliot,et al.  Ribosome flow model with extended objects , 2017, Journal of The Royal Society Interface.

[97]  W. Martin,et al.  The energetics of genome complexity , 2010, Nature.