The Key Parameters that Govern Translation Efficiency

How total mRNA influences cell growth

Experimental observations tracing back to the 1960s imply that ribosome quantities play a prominent role in determining a cell's growth. Nevertheless, in biologically relevant scenarios, growth can also be influenced by the levels of mRNA and RNA polymerase. Here, we construct a quantitative model of biosynthesis providing testable scenarios for these situations. The model explores a theoretically motivated regime where RNA polymerases compete for genes and ribosomes for transcripts and gives general expressions relating growth rate, mRNA concentrations, ribosome, and RNA polymerase levels. On general grounds, the model predicts how the fraction of ribosomes in the proteome depends on total mRNA concentration and inspects an underexplored regime in which the trade-off between transcript levels and ribosome abundances sets the cellular growth rate. In particular, we show that the model predicts and clarifies three important experimental observations, in budding yeast and Escherichia coli bacteria: i) that the growth-rate cost of unneeded protein expression can be affected by mRNA levels, ii) that resource optimization leads to decreasing trends in mRNA levels at slow growth, and iii) that ribosome allocation may increase, stay constant, or decrease, in response to transcription-inhibiting antibiotics. Since the data indicate that a regime of joint limitation may apply in physiological conditions and not only to perturbations, we speculate that this regime is likely self-imposed.

Weak Pinning and Long-Range Anticorrelated Motion of Phase Boundaries in Driven Diffusive Systems

We show that domain walls separating coexisting extremal current phases in driven diffusive systems exhibit complex stochastic dynamics with a subdiffusive temporal growth of position fluctuations due to long-range anticorrelated current fluctuations and a weak pinning at long times. This weak pinning manifests itself in a saturated width of the domain wall position fluctuations that increases sublinearly with the system size. As a function of time t and system size L, the width w(t,L) has a scaling behavior w(t,L)=L^{3/4}f(t/L^{9/4}), with f(u) constant for u≫1 and f(u)∼u^{1/3} for u≪1. An Orstein-Uhlenbeck process with long-range anticorrelated noise is shown to capture this scaling behavior. The exponent 9/4 is a new dynamical exponent for relaxation processes in driven diffusive systems.

Codon optimality modulates protein output by tuning translation initiation

The impact of synonymous codon choice on protein output has important implications for understanding endogenous gene expression and design of synthetic mRNAs. Previously, we used a neural network model to design a series of synonymous fluorescent reporters whose protein output in yeast spanned a seven-fold range corresponding to their predicted translation speed. Here, we show that this effect is not due primarily to the established impact of slow elongation on mRNA stability, but rather, that an active mechanism further decreases the number of proteins made per mRNA. We combine simulations and careful experiments on fluorescent reporters to argue that translation initiation is limited on non-optimally encoded transcripts. Using a genome-wide CRISPRi screen to discover factors modulating the output from non-optimal transcripts, we identify a set of translation initiation factors including multiple subunits of eIF3 whose depletion restored protein output of a non-optimal reporter. Our results show that codon usage can directly limit protein production, across the full range of endogenous variability in codon usage, by limiting translation initiation.

Model-based characterization of the equilibrium dynamics of transcription initiation and promoter-proximal pausing in human cells

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.

Single-file transport of binary hard-sphere mixtures through periodic potentials

Single-file transport occurs in various scientific fields, including diffusion through nanopores, nanofluidic devices, and cellular processes. We here investigate the impact of polydispersity on particle currents for single-file Brownian motion of hard spheres when they are driven through periodic potentials by a constant drag force. Through theoretical analysis and extensive Brownian dynamics simulations, we unveil the behavior of particle currents for random binary mixtures. The particle currents show a recurring pattern in dependence of the hard-sphere diameters and mixing ratio. We explain this recurrent behavior by showing that a basic unit cell exists in the space of the two hard-sphere diameters. Once the behavior of an observable inside the unit cell is determined, it can be inferred for any diameter. The overall variation of particle currents with the mixing ratio and hard-sphere diameters is reflected by their variation in the limit where the system is fully covered by hard spheres. In this limit, the currents can be predicted analytically. Our analysis explains the occurrence of pronounced maxima and minima of the currents by changes in the effective potential barrier for the center-of-mass motion.

Bayesian inference of polymerase dynamics over the exclusion process

Transcription is a complex phenomenon that permits the conversion of genetic information into phenotype by means of an enzyme called RNA polymerase, which erratically moves along and scans the DNA template. We perform Bayesian inference over a paradigmatic mechanistic model of non-equilibrium statistical physics, i.e. the asymmetric exclusion processes in the hydrodynamic limit, assuming a Gaussian process prior for the polymerase progression rate as a latent variable. Our framework allows us to infer the speed of polymerases during transcription given their spatial distribution, while avoiding the explicit inversion of the system's dynamics. The results, which show processing rates strongly varying with genomic position and minor role of traffic-like congestion, may have strong implications for the understanding of gene expression.

choros: correction of sequence-based biases for accurate quantification of ribosome profiling data

Ribosome profiling quantifies translation genome-wide by sequencing ribosome-protected fragments, or footprints. Its single-codon resolution allows identification of translation regulation, such as ribosome stalls or pauses, on individual genes. However, enzyme preferences during library preparation lead to pervasive sequence artifacts that obscure translation dynamics. Widespread over- and under-representation of ribosome footprints can dominate local footprint densities and skew estimates of elongation rates by up to five fold. To address these biases and uncover true patterns of translation, we present choros, a computational method that models ribosome footprint distributions to provide bias-corrected footprint counts. choros uses negative binomial regression to accurately estimate two sets of parameters: (i) biological contributions from codon-specific translation elongation rates; and (ii) technical contributions from nuclease digestion and ligation efficiencies. We use these parameter estimates to generate bias correction factors that eliminate sequence artifacts. Applying choros to multiple ribosome profiling datasets, we are able to accurately quantify and attenuate ligation biases to provide more faithful measurements of ribosome distribution. We show that a pattern interpreted as pervasive ribosome pausing near the beginning of coding regions is likely to arise from technical biases. Incorporating choros into standard analysis pipelines will improve biological discovery from measurements of translation.

Probiotic-fermented Pueraria lobata (Willd.) Ohwi alleviates alcoholic liver injury by enhancing antioxidant defense and modulating gut microbiota

BACKGROUND

Pueraria lobata (Willd.) Ohwi (PL) has been used in China to detoxify alcohol and protect the liver for millennia, though its mechanism of liver protection has not been elucidated. However, fermentation is considered to be one of the effective ways to enhance the efficacy of traditional Chinese medicine. The aim of this study was to investigate the hepatoprotective mechanism of probiotic-fermented PL (FPL). Sprague Dawley rats were administered with FPL followed by gavage of alcohol for seven consecutive days; following that, liver injury levels were evaluated in rats.

RESULTS

FPL ameliorated lipid accumulation and inflammation levels in rats. Meanwhile, the levels of ethanol dehydrogenase, acetaldehyde dehydrogenase, and cytochrome P4502E1 were elevated by FPL treatment. It was observed that the levels of catalase, superoxide dismutase, and glutathione peroxidase were elevated, and the expression of nuclear transcriptional factor (erythroid-derived 2)-like 2 (Nrf2) and heme oxygenase-1 genes and proteins were increased by FPL treatment, demonstrating that the Nrf2-mediated signal pathway was activated. Furthermore, FPL restored the composition of the gut microbiota with an increase in the abundances of Firmicutes and Lactobacillus and a decrease in the abundances of Bacteroidota and Akkermansia. Additionally, a strong correlation was found between the gut microbiota and the antioxidant parameters.

CONCLUSION

The results indicate that FPL possesses an excellent protective effect in alcoholic liver injury. Our findings are beneficial to the development of hepatoprotective nutraceuticals for alcoholics. © 2022 Society of Chemical Industry.

End-to-end pipeline for differential analysis of pausing in ribosome profiling data

Ribosome profiling is a powerful technique which maps the distribution of ribosomes along mRNAs to analyze translation genome-wide. Ribosome density can be affected by multiple factors, such as changes to translation initiation or elongation rates. We describe the application of a metric for identifying genes rate-limited by these rates by analyzing the relative distribution of ribosome footprints along transcripts. This protocol also details two sample analyses comparing gene translation efficiencies and the distribution of ribosome densities on downloadable datasets.

For complete details on the use and execution of this protocol, please refer to Flanagan et al. (2022).

•

Pipeline for measuring translation efficiency using ribosome profiling data

•

Quantification of differential ribosome distributions along mRNAs

•

Determine whether changes are consistent with altered initiation or elongation

Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

FMRP-dependent production of large dosage-sensitive proteins is highly conserved

Mutations in FMR1 are the most common heritable cause of autism spectrum disorder. FMR1 encodes an RNA-binding protein, FMRP, which binds to long, autism-relevant transcripts and is essential for normal neuronal and ovarian development. In contrast to the prevailing model that FMRP acts to block translation elongation, we previously found that FMRP activates the translation initiation of large proteins in Drosophila oocytes. We now provide evidence that FMRP-dependent translation is conserved and occurs in the mammalian brain. Our comparisons of the mammalian cortex and Drosophila oocyte ribosome profiling data show that translation of FMRP-bound mRNAs decreases to a similar magnitude in FMRP-deficient tissues from both species. The steady-state levels of several FMRP targets were reduced in the Fmr1 KO mouse cortex, including a ∼50% reduction of Auts2, a gene implicated in an autosomal dominant autism spectrum disorder. To distinguish between effects on elongation and initiation, we used a novel metric to detect the rate-limiting ribosome stalling. We found no evidence that FMRP target protein production is governed by translation elongation rates. FMRP translational activation of large proteins may be critical for normal human development, as more than 20 FMRP targets including Auts2 are dosage sensitive and are associated with neurodevelopmental disorders caused by haploinsufficiency.

Slow nucleosome dynamics set the transcriptional speed limit and induce RNA polymerase II traffic jams and bursts

Nucleosomes are recognized as key regulators of transcription. However, the relationship between slow nucleosome unwrapping dynamics and bulk transcriptional properties has not been thoroughly explored. Here, an agent-based model that we call the dynamic defect Totally Asymmetric Simple Exclusion Process (ddTASEP) was constructed to investigate the effects of nucleosome-induced pausing on transcriptional dynamics. Pausing due to slow nucleosome dynamics induced RNAPII convoy formation, which would cooperatively prevent nucleosome rebinding leading to bursts of transcription. The mean first passage time (MFPT) and the variance of first passage time (VFPT) were analytically expressed in terms of the nucleosome rate constants, allowing for the direct quantification of the effects of nucleosome-induced pausing on pioneering polymerase dynamics. The mean first passage elongation rate γ(h_c, h_o) is inversely proportional to the MFPT and can be considered to be a new axis of the ddTASEP phase diagram, orthogonal to the classical αβ-plane (where α and β are the initiation and termination rates). Subsequently, we showed that, for β = 1, there is a novel jamming transition in the αγ-plane that separates the ddTASEP dynamics into initiation-limited and nucleosome pausing-limited regions. We propose analytical estimates for the RNAPII density ρ, average elongation rate v, and transcription flux J and verified them numerically. We demonstrate that the intra-burst RNAPII waiting times t_in follow the time-headway distribution of a max flux TASEP and that the average inter-burst interval tIBI¯ correlates with the index of dispersion D_e. In the limit γ→0, the average burst size reaches a maximum set by the closing rate h_c. When α≪1, the burst sizes are geometrically distributed, allowing large bursts even while the average burst size NB¯ is small. Last, preliminary results on the relative effects of static and dynamic defects are presented to show that dynamic defects can induce equal or greater pausing than static bottle necks.

To perform specific functions, cells must express specific genes by copying the information in DNA into RNA via transcription. Structural proteins called nucleosomes are spaced every 200 base pairs along the length of a strand of DNA and play a crucial function in the regulation of gene activity by tightly binding DNA strands and condensing them into heterochromatin, preventing transcription by RNA polymerase II (RNAPII). Even on active genes where nucleosomes are loosely attached to DNA strands, the wrapping and unwrapping of nucleosomes pause transcription as RNAPII passes by. Previous mathematical models of transcription have compared this biological process to traffic on a one lane highway without obstructions. In contrast, our proposed model simulates transcription like traffic in a grid system where nucleosomes can be thought of as pedestrians or other vehicles crossing the road at regularly spaced intersections. Just as side street traffic and pedestrian crossings can cause cars to form convoys and cause jams limiting the max speed in an area, nucleosomes can cause RNAPII to form convoys that lead to bursts of mRNA production and limit the average polymerase flux through the gene.

Analysis of Ribosome Profiling Data

Ribosome profiling methods are based on high-throughput sequencing of ribosome-protected mRNA footprints and allow to study in detail translational changes. Bioinformatic and statistical tools are necessary to analyze sequencing data. Here, we describe our developed methods for a fast and reliable quality control of ribosome profiling data, to efficiently visualize ribosome positions and to estimate ribosome speed in an unbiased way. The methodology described here is applicable to several genetic and environmental conditions including stress and are based on the R package RiboVIEW and calculation of quantitative estimates of local and global translation speed, based on a biophysical model of translation dynamics.

Density- and elongation speed-dependent error correction in RNA polymerization

Backtracking of RNA polymerase (RNAP) is an important pausing mechanism during DNA transcription that is part of the error correction process that enhances transcription fidelity. We model the backtracking mechanism of RNA polymerase, which usually happens when the polymerase tries to incorporate a noncognate or "mismatched" nucleotide triphosphate. Previous models have made simplifying assumptions such as neglecting the trailing polymerase behind the backtracking polymerase or assuming that the trailing polymerase is stationary. We derive exact analytic solutions of a stochastic model that includes locally interacting RNAPs by explicitly showing how a trailing RNAP influences the probability that an error is corrected or incorporated by the leading backtracking RNAP. We also provide two related methods for computing the mean times for error correction and incorporation given an initial local RNAP configuration. Using these results, we propose an effective interacting-RNAP lattice that can be readily simulated.

A deterministic model for non-monotone relationship between translation of upstream and downstream open reading frames

Totally asymmetric simple exclusion process (TASEP) modelling was shown to offer a parsimonious explanation for the experimentally confirmed ability of a single upstream open reading frames (uORFs) to upregulate downstream translation during the integrated stress response. As revealed by numerical simulations, the model predicts that reducing the density of scanning ribosomes upstream of certain uORFs increases the flow of ribosomes downstream. To gain a better insight into the mechanism which ensures the non-monotone relation between the upstream and downstream flows, in this work, we propose a phenomenological deterministic model approximating the TASEP model of the translation process. We establish the existence of a stationary solution featuring the decreasing density along the uORF for the deterministic model. Further, we find an explicit non-monotone relation between the upstream ribosome density and the downstream flow for the stationary solution in the limit of increasing uORF length and increasingly leaky initiation. The stationary distribution of the TASEP model, the stationary solution of the deterministic model and the explicit limit are compared numerically.

Quantitative elucidation of associations between nucleotide identity and physicochemical properties of amino acids and the functional insight

Studies on codon property would deepen our understanding of the origin of primitive life and enlighten biotechnical application. Here, we proposed a quantitative measurement of codon-amino acid association and found that seven out of 13 physicochemical properties have stronger associations with the nucleotide identity at the second codon position, indicating that protein structure and function may associate more closely with it than the other two sites. When extending the effect of codon-amino acid association to protein level, it was found that the correlation between the second codon position (measured by the relative frequencies of nucleobase T and A at this codon site) and hydrophobicity (by the form of GRAVY value) became stronger with 96% genomes having R > 0.90 and p < 1e-60. Furthermore, we revealed that informational genes encoding proteins have lower GRAVY values than operational proteins (p < 3e-37) in both prokaryotic and eukaryotic genomes. The above results reveal a complete link from codon identity (A₂ versus T₂) to amino acid property (hydrophilic versus hydrophobic) and then to protein functions (informational versus operational). Hence, our work may help to understand how the nucleotide sequence determines protein function.

Zebrafish as a Model to Evaluate a CRISPR/Cas9-Based Exon Excision Approach as a Future Treatment Option for EYS-Associated Retinitis Pigmentosa

Retinitis pigmentosa (RP) is an inherited retinal disease (IRD) with an overall prevalence of 1 in 4000 individuals. Mutations in EYS (Eyes shut homolog) are among the most frequent causes of non-syndromic autosomal recessively inherited RP and act via a loss-of-function mechanism. In light of the recent successes for other IRDs, we investigated the therapeutic potential of exon skipping for EYS-associated RP. CRISPR/Cas9 was employed to generate zebrafish from which the region encompassing the orthologous exons 37-41 of human EYS (eys exons 40-44) was excised from the genome. The excision of these exons was predicted to maintain the open reading frame and to result in the removal of exactly one Laminin G and two EGF domains. Although the eys^Δexon40-44 transcript was found at levels comparable to wild-type eys, and no unwanted off-target modifications were identified within the eys coding sequence after single-molecule sequencing, Eys^Δexon40-44 protein expression could not be detected. Visual motor response experiments revealed that eys^Δexon40-44 larvae were visually impaired and histological analysis revealed a progressive degeneration of the retinal outer nuclear layer in these zebrafish. Altogether, the data obtained in our zebrafish model currently provide no indications for the skipping of EYS exons 37-41 as an effective future treatment strategy for EYS-associated RP.

Resistance Exercise, Aging, Disuse, and Muscle Protein Metabolism

Skeletal muscle is the organ of locomotion, its optimal function is critical for athletic performance, and is also important for health due to its contribution to resting metabolic rate and as a site for glucose uptake and storage. Numerous endogenous and exogenous factors influence muscle mass. Much of what is currently known regarding muscle protein turnover is owed to the development and use of stable isotope tracers. Skeletal muscle mass is determined by the meal- and contraction-induced alterations of muscle protein synthesis and muscle protein breakdown. Increased loading as resistance training is the most potent nonpharmacological strategy by which skeletal muscle mass can be increased. Conversely, aging (sarcopenia) and muscle disuse lead to the development of anabolic resistance and contribute to the loss of skeletal muscle mass. Nascent omics-based technologies have significantly improved our understanding surrounding the regulation of skeletal muscle mass at the gene, transcript, and protein levels. Despite significant advances surrounding the mechanistic intricacies that underpin changes in skeletal muscle mass, these processes are complex, and more work is certainly needed. In this article, we provide an overview of the importance of skeletal muscle, describe the influence that resistance training, aging, and disuse exert on muscle protein turnover and the molecular regulatory processes that contribute to changes in muscle protein abundance. © 2021 American Physiological Society. Compr Physiol 11:2249-2278, 2021.

Polysomally protected viruses

It is conceivable that an RNA virus could use a polysome, that is, a string of ribosomes covering the RNA strand, to protect the genetic material from degradation inside a host cell. This paper discusses how such a virus might operate, and how its presence might be detected by ribosome profiling. There are two possible forms for such apolysomally protected virus, depending upon whether just the forward strand or both the forward and complementary strands can be encased by ribosomes (these will be termed type 1 and type 2, respectively). It is argued that in the type 2 case the viral RNA would evolve anambigrammaticproperty, whereby the viral genes are free of stop codons in a reverse reading frame (with forward and reverse codons aligned). Recent observations of ribosome profiles of ambigrammatic narnavirus sequences are consistent with our predictions for the type 2 case.

Full-length ribosome density prediction by a multi-input and multi-output model

Translation elongation is regulated by a series of complicated mechanisms in both prokaryotes and eukaryotes. Although recent advance in ribosome profiling techniques has enabled one to capture the genome-wide ribosome footprints along transcripts at codon resolution, the regulatory codes of elongation dynamics are still not fully understood. Most of the existing computational approaches for modeling translation elongation from ribosome profiling data mainly focus on local contextual patterns, while ignoring the continuity of the elongation process and relations between ribosome densities of remote codons. Modeling the translation elongation process in full-length coding sequence (CDS) level has not been studied to the best of our knowledge. In this paper, we developed a deep learning based approach with a multi-input and multi-output framework, named RiboMIMO, for modeling the ribosome density distributions of full-length mRNA CDS regions. Through considering the underlying correlations in translation efficiency among neighboring and remote codons and extracting hidden features from the input full-length coding sequence, RiboMIMO can greatly outperform the state-of-the-art baseline approaches and accurately predict the ribosome density distributions along the whole mRNA CDS regions. In addition, RiboMIMO explores the contributions of individual input codons to the predictions of output ribosome densities, which thus can help reveal important biological factors influencing the translation elongation process. The analyses, based on our interpretable metric named codon impact score, not only identified several patterns consistent with the previously-published literatures, but also for the first time (to the best of our knowledge) revealed that the codons located at a long distance from the ribosomal A site may also have an association on the translation elongation rate. This finding of long-range impact on translation elongation velocity may shed new light on the regulatory mechanisms of protein synthesis. Overall, these results indicated that RiboMIMO can provide a useful tool for studying the regulation of translation elongation in the range of full-length CDS.

Translation elongation is a process in which amino acids are linked into proteins by ribosomes in cells. Translation elongation rates along the mRNAs are not constant, and are regulated by a series of mechanisms, such as codon rarity and mRNA stability. In this study, we modeled the translation elongation process at a full-length coding sequence level and developed a deep learning based approach to predict the translation elongation rates from mRNA sequences, through extracting the regulatory codes of elongation rates from the contextual sequences. The analyses, based on our interpretable metric named codon impact score, for the first time (to the best of our knowledge), revealed that in addition to the neighboring codons of the ribosomal A sites, the remote codons may also have an important impact on the translation elongation rates. This new finding may stimulate additional experiments and shed light on the regulatory mechanisms of protein synthesis.

Molecular motor traffic with a slow binding site

We discuss how the presence of a slow binding site in molecular motor traffic gives rise to defect-induced "traffic jams" that have properties different from those of the well-studied boundary-induced jams that originate from an imbalance between initiation and termination. To this end we analyze in detail the stationary distribution of a lattice gas model for traffic of molecular motors with a defect. In particular, we obtain analytically the exact spatial distribution of motors, the probability distribution of the random position of the molecular traffic jam and we report unexpected spatial anticorrelations between local molecular motor densities near the defect.

EGGTART: A tool to visualize the dynamics of biophysical transport under the inhomogeneous l-TASEP

The totally asymmetric simple exclusion process (TASEP), which describes the stochastic dynamics of interacting particles on a lattice, has been actively studied over the past several decades and applied to model important biological transport processes. Here, we present a software package, called EGGTART (Extensive GUI gives TASEP-realization in Real Time), which quantifies and visualizes the dynamics associated with a generalized version of the TASEP with an extended particle size and heterogeneous jump rates. This computational tool is based on analytic formulas obtained from deriving and solving the hydrodynamic limit of the process. It allows an immediate quantification of the particle density, flux, and phase diagram, as a function of a few key parameters associated with the system, which would be difficult to achieve via conventional stochastic simulations. Our software should therefore be of interest to biophysicists studying general transport processes and can in particular be used in the context of gene expression to model and quantify mRNA translation of different coding sequences.

What determines eukaryotic translation elongation: recent molecular and quantitative analyses of protein synthesis

Protein synthesis from mRNA is an energy-intensive and tightly controlled cellular process. Translation elongation is a well-coordinated, multifactorial step in translation that undergoes dynamic regulation owing to cellular state and environmental determinants. Recent studies involving genome-wide approaches have uncovered some crucial aspects of translation elongation including the mRNA itself and the nascent polypeptide chain. Additionally, these studies have fuelled quantitative and mathematical modelling of translation elongation. In this review, we provide a comprehensive overview of the key determinants of translation elongation. We discuss consequences of ribosome stalling or collision, and how the cells regulate translation in case of such events. Next, we review theoretical approaches and widely used mathematical models that have become an essential ingredient to interpret complex molecular datasets and study translation dynamics quantitatively. Finally, we review recent advances in live-cell reporter and related analysis techniques, to monitor the translation dynamics of single cells and single-mRNA molecules in real time.

Inferring efficiency of translation initiation and elongation from ribosome profiling

One of the main goals of ribosome profiling is to quantify the rate of protein synthesis at the level of translation. Here, we develop a method for inferring translation elongation kinetics from ribosome profiling data using recent advances in mathematical modelling of mRNA translation. Our method distinguishes between the elongation rate intrinsic to the ribosome’s stepping cycle and the actual elongation rate that takes into account ribosome interference. This distinction allows us to quantify the extent of ribosomal collisions along the transcript and identify individual codons where ribosomal collisions are likely. When examining ribosome profiling in yeast, we observe that translation initiation and elongation are close to their optima and traffic is minimized at the beginning of the transcript to favour ribosome recruitment. However, we find many individual sites of congestion along the mRNAs where the probability of ribosome interference can reach \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$50\%$\end{document}. Our work provides new measures of translation initiation and elongation efficiencies, emphasizing the importance of rating these two stages of translation separately.

Current-density relation in the exclusion process with dynamic obstacles

We investigate the totally asymmetric simple exclusion process (TASEP) in the presence of obstacles that dynamically bind and unbind from the lattice. The model is motivated by biological processes such as transcription in the presence of DNA-binding proteins. Similar models have been studied before using the mean-field approximation, but the exact relation between the particle current and density remains elusive. Here, we first show using extensive Monte Carlo simulations that the current-density relation in this model assumes a quasiparabolic form similar to that of the ordinary TASEP without obstacles. We then attempt to explain this relation using exact calculations in the limit of low and high density of particles. Our results suggest that the symmetric, quasiparabolic current-density relation arises through a nontrivial cancellation of higher-order terms, similarly as in the standard TASEP.

Dynamics of ribosomes in mRNA translation under steady- and nonsteady-state conditions

Recent advances in DNA sequencing and fluorescence imaging have made it possible to monitor the dynamics of ribosomes actively engaged in messenger RNA (mRNA) translation. Here, we model these experiments within the inhomogeneous totally asymmetric simple exclusion process (TASEP) using realistic kinetic parameters. In particular, we present analytic expressions to describe the following three cases: (a) translation of a newly transcribed mRNA, (b) translation in the steady state and, specifically, the dynamics of individual (tagged) ribosomes, and (c) runoff translation after inhibition of translation initiation. In cases (b) and (c) we develop an effective medium approximation to describe many-ribosome dynamics in terms of a single tagged ribosome in an effective medium. The predictions are in good agreement with stochastic simulations.

The yeast exoribonuclease Xrn1 and associated factors modulate RNA polymerase II processivity in 5' and 3' gene regions

mRNA levels are determined by the balance between mRNA synthesis and decay. Protein factors that mediate both processes, including the 5'-3' exonuclease Xrn1, are responsible for a cross-talk between the two processes that buffers steady-state mRNA levels. However, the roles of these proteins in transcription remain elusive and controversial. Applying native elongating transcript sequencing (NET-seq) to yeast cells, we show that Xrn1 functions mainly as a transcriptional activator and that its disruption manifests as a reduction of RNA polymerase II (Pol II) occupancy downstream of transcription start sites. By combining our sequencing data and mathematical modeling of transcription, we found that Xrn1 modulates transcription initiation and elongation of its target genes. Furthermore, Pol II occupancy markedly increased near cleavage and polyadenylation sites in xrn1Δ cells, whereas its activity decreased, a characteristic feature of backtracked Pol II. We also provide indirect evidence that Xrn1 is involved in transcription termination downstream of polyadenylation sites. We noted that two additional decay factors, Dhh1 and Lsm1, seem to function similarly to Xrn1 in transcription, perhaps as a complex, and that the decay factors Ccr4 and Rpb4 also perturb transcription in other ways. Interestingly, the decay factors could differentiate between SAGA- and TFIID-dominated promoters. These two classes of genes responded differently to XRN1 deletion in mRNA synthesis and were differentially regulated by mRNA decay pathways, raising the possibility that one distinction between these two gene classes lies in the mechanisms that balance mRNA synthesis with mRNA decay.