text:
View all Nature Research journals
Gene length as a regulator for ribosome recruitment and protein synthesis: theoretical insights
Lucas D. Fernandes ORCID: orcid.org/0000-0001-6053-22771,2,
Alessandro P. S. de Moura2 &
Luca Ciandrini ORCID: orcid.org/0000-0001-5859-27643,4
Scientific Reports volume 7, Article number: 17409 (2017) Cite this article
Protein synthesis rates are determined, at the translational level, by properties of the transcript's sequence. The efficiency of an mRNA can be tuned by varying the ribosome binding sites controlling the recruitment of the ribosomes, or the codon usage establishing the speed of protein elongation. In this work we propose transcript length as a further key determinant of translation efficiency. Based on a physical model that considers the kinetics of ribosomes advancing on the mRNA and diffusing in its surrounding, as well as mRNA circularisation and ribosome drop-off, we explain how the transcript length may play a central role in establishing ribosome recruitment and the overall translation rate of an mRNA. According to our results, the proximity of the 3′ end to the ribosomal recruitment site of the mRNA could induce a feedback in the translation process that would favour the recycling of ribosomes. We also demonstrate how this process may be involved in shaping the experimental ribosome density-gene length dependence. Finally, we argue that cells could exploit this mechanism to adjust and balance the usage of its ribosomal resources.
mRNA translation is, together with transcription, the pillar of the central dogma of molecular biology. In spite of its key role in protein synthesis, the accurate understanding of its dynamical details still remains elusive at the present time, and the sequence determinants of mRNA translation efficiency are not fully understood1,2. Initiation of translation regulates the recruitment of ribosomes and it is believed to be modulated by mRNA secondary structures3,4, while protein elongation is mainly considered to be regulated by tRNA abundances determining the pace of the ribosome5,6,7. The individual steps of translation are thought to be well understood, yet there is no reliable approach quantitatively predicting the overall protein synthesis rates for a given transcript.
A better understanding of the molecular mechanisms of mRNA translation will unravel the physiological determinants of translation efficiency. Besides, this knowledge will be extremely useful in developing applications in synthetic biology and will allow tight control on the average production of a protein and on its expression noise.
The translation efficiency of a transcript is often identified with its experimentally observed polysome state (a transcript with two or more ribosomes; monosome when there is only one ribosome residing on it), meaning that transcripts with high ribosome density are more efficiently translated8. Remarkably, many experimental observations show that the ribosome density is related to the length L of the coding sequence (CDS): the longer the mRNA, the smaller the ribosomal density. This indicates the presence of a length-dependent control of translation. As we show in Fig. 1, the observation that average ribosomal densities ρ strongly anti-correlate to CDS lengths L appears to be a conserved feature across many organisms, ranging from unicellular systems such as L. lactis 9, S. cerevisiae 10,11,12 or P. falciparum 13, to more complex organisms such as mouse and human cells14,15. The common traits in the density-length dependence suggest that this relationship is dictated by universal mechanisms underlying the translation process.
Ribosome density vs CDS length for different datasets. Blue diamonds (mice14) and fuchsia down triangles (P. falciparum 13) are individual genes, while yellow circles (HEK293T15), green squares (S. cerevisiae 10) and red triangles (S. cerevisiae 11) are length-binned data for the entire genome, with the error bars representing the standard deviation for each bin. The grey line indicates the behaviour of a power law with exponent −1.
However, this remark has been strangely overlooked in the literature (with the exception of Guo et al.16), particularly in the theoretical literature trying to provide models of mRNA translation. A few hypotheses have been proposed to justify the emergence of the length dependence of the ribosome density, which requires a regulation apparatus acting at the initiation17 or at the elongation stage12. These hypotheses have not been examined with the support of a mechanistic model and a mathematical approach. In contrast with previous studies, here we explain qualitatively and quantitatively the relationship between ribosome density and CDS length, making the point that transcript length is a critical determinant of protein synthesis rates.
In Fig. 1 we show a log-log plot of measured ribosomal densities as a function of the CDS length for different organisms. The figure suggests a power-law behaviour (L −1 is drawn as reference). However, extracting a scaling law from this kind of data is a phenomenological (and probably a too simplistic) description of this relationship that, moreover, can only be measured for a few orders of magnitude in L.
Instead, in this paper we propose a mechanistic explanation for the length dependence of translation that is found in experimental data: we describe how the proximity of the mRNA ends increases the local concentration of ribosomes close to their binding sites via a feedback mechanism, favouring their recruitment (and recycling) in short transcripts. In the last part of this work we show how this mechanism could be exploited by the cell to adjust and balance its ribosomal resources.
A stochastic model of translation
The translation of an mRNA is a three-step process, as sketched in Fig. 2A: during initiation, the ribosomal subunits are recruited and the full functional ribosome is assembled on the START codon, ready to translate the transcript; then the ribosome proceeds to elongation, assembling the protein amino acid by amino acid according to the nucleotide sequence of the mRNA; the ribosome eventually detaches when it reaches the STOP codon (termination).
Sketch of the translation process and models. The three-steps of the translation process (A): initiation (in the model approximated by a one-step process with rate α), elongation and termination (β). In the standard exclusion process (B) particles can enter the beginning of the lattice with a rate α, move from one site to the next one with rate p (provided that it is not occupied by another particle), and then exit on the last site with rate β. In this study we consider a more refined version of the model (C) in which ribosomes cover \(\ell \) sites (codons), advance one site at a time, and the unidimensional lattice is placed in a three-dimensional environment. R represents the end-to-end distance between the 5′and the 3′ region, and a is the radius of the reaction volume for initiation. The dashed grey line represents a possible diffusive trajectory of the ribosomal subunits leaving the transcript and being re-absorbed (recycled) in the reaction volume around the first site of the lattice.
We model translating ribosomes as particles moving on a unidimensional discrete track of length L representing the mRNA, as depicted in Fig. 2B. In this model particles are injected from one side of the lattice (the 5′ end of the mRNA) with a rate α, then advance one site (one codon) with rate p only if the arrival site is empty, and are removed at the last site (STOP codon) with rate β (Fig. 2B). The first step mimics the recruitment of ribosomes (initiation); elongation is given by the dynamics of ribosomes in the bulk; the exit of particles from the last site represents the termination. MacDonald and coworkers introduced this class of model at the end of the 60′s precisely in an attempt to mathematically describe the process of mRNA translation18. Since then, under the name of exclusion process, this model has been extended and thoroughly studied from a theoretical point of view; it became an emblematic framework in out-of-equilibrium physics19, for which an exact solution is known in the simplest formulation20.
In the last years, revamped extensions of the exclusion process have appeared in the literature, developed to provide more quantitative models of translation21,22,23,24, and many works are nowadays implicitly based on this framework25,26. Here we first look into a variant of the exclusion process that considers particles covering \(\ell =10\) sites of the track27, as the ribosome footprints cover around 28 nucleotides12 (see Fig. 2C). Details of the model and simulations can be found in the Materials and Methods section and in the Supplementary Material.
Translation efficiency corresponds to the ribosomal current
From the analytical solution of the exclusion process or by Monte Carlo simulations we can estimate, as a function of the initiation rate α, termination rate β and codon elongation rate p, the two quantities of interest for the translation process: the ribosomal density ρ(α, β, p), defined as the average number of ribosomes N divided by the CDS length L, and the ribosomal current J(α, β, p), defined as the average number of ribosomes advancing one site in a unit of time. Those quantities can be compared to experimental measurements of polysome profiles and protein production rates. The ribosomal current J corresponds in fact to the protein production rate per mRNA (proteins produced per unit time per mRNA), and we choose to identify it as a better descriptor for the translation efficiency.
The same analytical solution also gives the dependence of current and density on the system's parameters. The model's behaviour has been largely studied in the literature as a function of the dimensionless parameters \(\bar{\alpha }\equiv \alpha /p\) and \(\bar{\beta }\equiv \beta /p\), and depending to their values different phases can be observed (see Supplementary Material).
Since the termination rate β is not limiting translation, the translation efficiency and the density do not depend on \(\bar{\beta }\): the system is in the so-called low density phase, in which the density should always be smaller than ~0.076 consistently with Fig. 1. From the analytical solutions one can show that efficiency J and density ρ only depend on \(\bar{\alpha }\equiv \alpha /p\) in this regime (see Supplementary Material). Hence \(\rho (\alpha ,\beta ,p)=\rho (\bar{\alpha })\) and \(J(\alpha ,\beta ,p)=J(\bar{\alpha })\).
In order to make the model more realistic and compare it to experimental datasets, we need to determine the initiation rate α. The estimation of the translation initiation rate has been previously attempted for a few organisms24,28, and these studies have found a dependence of the initiation rate on the transcript length: the longer the transcript, the weaker the initiation. Here we propose a model that is able to explain this observation by coupling the translation process, in particular translation initiation, to the three-dimensional conformation of the mRNA. We will show how a feedback mechanism between the ribosomal current leaving the end of the mRNA and the initiation process, which is controlled by the polysome compaction, could induce a length-dependent initiation rate and hence a length-dependent density. Before that, we need to introduce some basic properties of the transcript's spatial conformation.
Transcript end-to-end distance depends on the polysome size
We consider the transcript as a polymer that assumes different spatial conformations and a characteristic 5′-3′ end-to-end distance R (see Fig. 2C). An undecorated mRNA (without ribosomes on it) can be considered as a polymer with a persistence length \({l}_{p}\simeq 1\) nm \(\simeq 1\) codon29, and the average end-to-end distance R can be estimated from basic principles of polymer physics. By assuming an underlying random walk one obtains that the end-to-end distance R depends on the length L of the mRNA as
$$R=\sqrt{2L{l}_{p}}\mathrm{.}$$
However, the stiffness produced by the large size \(\ell \) of the ribosomes can drastically change the persistence length of the mRNA. We assume that the persistence length of an mRNA depends on the polysome state via an average between \(\ell \) (a typical ribosome footprint) and l p (persistence length of an empty mRNA), weighted by the fraction \(f=\rho \ell \) of the transcript covered by ribosomes (at the steady state). After these considerations we write the effective persistence length of the mRNA as
$$\begin{array}{rcl}{l}_{{\rm{eff}}} & = & f\ell +\mathrm{(1}-f){l}_{p}\\ & = & {\ell }^{2}\rho +\mathrm{(1}-\rho \ell ){l}_{p},\end{array}$$
which is equal to l p when the mRNA is empty and reaches \(\ell \) when the ribosomal density attains its maximal value \(\rho =\mathrm{1/}\ell \). Substituting this value of the effective persistence length into Eq. (1) we obtain the average end-to-end distance of a polysome as a function of the density ρ and the length L:
$$R=\sqrt{2L}{[{\ell }^{2}\rho +\mathrm{(1}-\rho \ell ){l}_{p}]}^{\mathrm{1/2}}\mathrm{.}$$
This way, we have used a coarse-grained model to couple the state of the polysome to its spatial conformation. Intuitively, Eq. (3) means that a translated transcript with many ribosomes on it will be more stretched (so the distance R between 5′ and 3′ ends will be larger) compared to a situation with a small ribosomal density or an empty mRNA. We have neglected potential formation of secondary structures inside the coding region. Such structures would only slightly decrease the effective length of the sequence (a few codons) and we treat only translationally active transcripts, meaning that moving ribosomes (10–20 codons/s) continuously unfold those structures.
When initiation is limiting, the density ρ of ribosomes is fixed by \(\bar{\alpha }\), and via Eq. (3) we are hence able to determine the dependence of the end-to-end distance as a function of the initiation rate and the length of the transcript, \(R(\bar{\alpha },L)\).
Initiation can be enhanced by a length-dependent feedback mechanism
We will consider that the magnitude of the initiation rate α is determined by the concentration c of free ribosomal subunits via a first-order rate equation: α = α 0 c, with α 0 being the initiation rate constant depending, for instance, on the affinity between the ribosome and the 5′UTR binding site on the mRNA. The concentration c is the local concentration of subunits in the reaction volume of radius a around the ribosome binding site at the 5′ end of the transcript; we introduce c ∞ as the homogeneous concentration of free subunits far from the transcript. The local concentration c in the reaction volume is affected by the ribosomes terminating translation at the 3′ end of a transcript, then diffusing into the reaction volume and contributing to the abundance of free ribosomal subunits in this volume. The contribution δ R to the local concentration c due to this feedback mechanism depends on the end-to-end distance R previously calculated in Eq. (3). It can be shown for different organisms, considering typical values of transcript numbers and cytoplasm volume, that the average end-to-end distance is smaller than the average separation between transcripts. This corroborates the intrinsic assumption that the translating processes of distinct mRNAs do not interfere with each other. Thus each individual mRNA can be thought of as a sink-source system for ribosomes (the ribosome binding site representing the sink and the ribosome termination site representing the source) immersed in an environment with a constant background ribosomal concentration c ∞. Hence, the local ribosome concentration around the ribosome binding site can be written as c(R) = c ∞ + δ R .
By considering ribosomes as particles performing free diffusion when they are not bound to the mRNA (see the Supplementary Material) we can obtain a mathematical expression of the initiation rate as a function of the system's parameters. It can be shown (see Supplementary Material) that the increase δ R due to the source at the end of transcript is \({\delta }_{R}\sim J/DR\), where D is the diffusion coefficient of ribosomes (~0.04 μm2/s30). Considering that the protein production rate per transcript is J~0.1−10 proteins/s24,31, and that in E.coli \({c}_{\infty }\sim 0.5-1\cdot {10}^{3}\) free ribosomes/μm 3 32, for a typical mRNA of length ~300 codons this back-of-the envelope calculation leads to \({\delta }_{R}/{c}_{\infty }\sim 0.1-10\), meaning that the concentration increase is at least comparable with the bulk ribosomal concentration and the mechanism proposed is relevant in the biological regime.
Regarding translation as a steady state process with ribosomal density and current values given by ρ and J respectively, we find the initiation parameter \(\bar{\alpha }\) (see Supplementary Material for a complete derivation):
$$\bar{\alpha }={\bar{\alpha }}_{0}({c}_{\infty }+{\delta }_{R})={\bar{\alpha }}_{\infty }+\lambda \frac{J(\bar{\alpha })}{R(\bar{\alpha },L)},$$
where we have emphasised the dependence of the ribosomal current J on \(\bar{\alpha }\). Similarly, the end-to-end distance R also depends on \(\bar{\alpha }\) and L as shown in Eq. (3). We highlight that the parameter \({\bar{\alpha }}_{0}\) and thus \({\bar{\alpha }}_{\infty }\) and λ depend on the binding between the ribosome and the mRNA, which is supposed to be mainly regulated by mRNA secondary structures. We will consider the parameters of the model to be independent on the transcript length. This is justified by a weak correlation (Pearson r = −0.01, p-value 0.5)33 between free energies of secondary structures in the 5′UTRs and the transcript length L (see also Figure S6 in the Supplementary Material). The parameter \({\bar{\alpha }}_{\infty }\) is adimensional and \({\alpha }_{\infty }=p{\bar{\alpha }}_{\infty }\) represents the initiation rate without the feedback mechanism between the current and the initiation, or equivalently when the end-to-end distance R is large enough to make this mechanism negligible. The parameter λ characterises the strength of the feedback. It has the dimensions of a length and it corresponds to the typical separation between 5′ and 3′ below which the feedback mechanism becomes relevant. We measure this parameter in units of codon length, which roughly corresponds to 1 nm.
Consequently, the current of ribosomes leaving the end of a transcript increases the concentration of ribosomal subunits around their binding region; through modulation of the mRNA stiffness due to the ribosome load, this feedback leads to initiation rates that are strongly length-dependent. Equation (4) is an implicit equation that can be numerically solved to obtain the initiation rate, and thus the density \(\rho (\bar{\alpha })\) and the current as a function of L. To check the validity of our analytical calculations, we also developed a simulation scheme that allows us to fix the initiation rate via a self-consistent method (see Methods section).
Although Eq. (4) considers the ribosome as a single diffusing particle, we can explicitly consider the diffusion of the two ribosomal subunits. This would generate a dependence on R 2 instead of R in Eq. (4). However, the qualitative behaviour of our results does not significantly change and for the sake of simplicity we decided to present the outcomes of the theory described by Eq. (4). We include the analysis of this more refined model in the Supplementary Material (see also Figures S9 and S12).
Experimental density-length dependence emerges from initiation enhanced effects
We then compare the outcome of the model to experimental measurements of ribosome densities. The result of this analysis is shown in Fig. 3: the predicted ribosome density is quantitatively comparable to the experimental quantification, and our mechanistic model based on basic physical principles is able to capture the length dependence of the ribosome load. Our theory can therefore explain the observed length dependence of the translational properties.
Comparison between theory and experimental ribosome densities in yeast (A,B)10,11 and human embryonic kidney cells (C)15. The symbols and datasets correspond to the ones of Fig. (1). The grey lines represent the best fit of the model (the parameter values are written in each panel), while the shadow areas correspond to the regions spanned within the margin of error of the estimated \({\bar{\alpha }}_{\infty }\). Orange circles are the outcome of stochastic simulations used to test the numerical solution of the equation using the parameters obtained from the best fit.
We fit the expression of \(\rho (\bar{\alpha })\) (continuous lines in Fig. 1) and obtain the two parameters \({\bar{\alpha }}_{\infty }\) and λ for three available datasets (two yeast datasets10,11 and a human embryonic kidney cells dataset15), then check the accuracy of the solution with the stochastic simulation scheme developed as explained in the Material and Methods section. Taking the standard errors of the parameter estimation there is no significant deviation between the data and the fitted model (details can be found in the Materials and Methods). However, for the last set, there is a stronger dependence of the solution on the parameter \({\bar{\alpha }}_{\infty }\): shaded regions in Fig. 3 represent the solution considering the standard error of this parameter. This could explain the slight deviation between theory and experimental data for large mRNAs. We also emphasise that the parameter \({\bar{\alpha }}_{\infty }\), depending on the global availability of ribosomes, is supposed to be the most affected by experimental variations (for instance by growth rate dependence or cell cycle stage). The estimation of \({\bar{\alpha }}_{\infty }\) is more accurate when we take into account the diffusion of the two subunits (see Supplementary Material).
Our simulations also allow us to extract the amount of ribosomes bound to a transcript, from which we can extract the monosome:polysome ratio, which is also subject to length-effects (see Supplementary Material). By increasing the CDS length we observe a reduction of the monosome:polysome ratio following a power-law like behaviour. A recent work34 has identified, by merging polysome and ribosome profiling techniques, the amount of active monosomes, i.e. mRNAs with only one ribosome elongating the protein. Consistent with our findings, the monosome:polysome ratio also shows signs of a marked anti-correlation with the mRNA length.
Ribosome drop-off cannot alone be responsible for the density-length dependence
In this section we study the consequences of ribosome drop-off35,36, thus far neglected in our model, on the observed density-length dependence. In order to do that, we performed simulations including ribosome drop-off at a rate δ = 10−3 s−1. This is justified by the estimated drop-off rates of the order of 10−4/codon35,36, and by the codon elongation rates we considered in this paper of the order of 10 codons s−1. Figure 4 shows that the simulations of the process with drop-off (full lines) do not largely differ from the model without drop-off (dashed lines), and the deviations starts to become relevant for large sequences (order 104 codons). For such lengths the extended model is actually reproducing the experimental behaviour even better than the model without drop-off.
Model with ribosome drop-off. Green lines correspond to the Arava dataset10, while red ones correspond to the Mackay dataset11 Symbols of experimental points in yeast correspond to the ones of Fig. 1, while dashed lines represent the solutions of the model described in the previous section with the same parameters used in panels A and B of Fig. 3. The continuous lines are the outcome of simulations of our model including ribosome drop-off at a rate of 10−3 s−1, and the grey line shows the outcome of the simulations with drop-off only.
To further exclude the possibility that the length dependence rises from ribosome prematurely leaving the transcript, we simulated ribosome drop-off occurring during the translation process without the feedback mechanism that we propose (basically, we simulated a standard TASEP with particle detachment rate δ). This is represented by the grey line in Fig. 4. We were not able to obtain any length dependence for biologically relevant values of the drop-off rate (see Supplementary Figure S10), and we can therefore conclude that the behaviour observed in Fig. 1 cannot originate by ribosomal drop-off alone.
mRNA circularisation does not change the phenomenology of the model
The 5′ and 3′ end of eukaryotic mRNAs interact with each other via protein-protein interaction, for instance between the Poly(A)-binding protein PAPB and the initiation factor eIF4F bound at the 5′ cap; this coupling is believed to induce the formation of transcripts with circular structures37. However, depending on the energies at play, the transcript dynamically switches between an open, linear state, and a circularised state (see Fig. 5A). When in the circularised state, the end-to-end distance R will be only of a few nanometers (the order of magnitude of the two molecular partners supposedly involved in this interactions). We will denote as d the distance between 5′ and 3′ in the circularised state. The mRNA is found in its circularised state with a given probability P c , and in an open state with probability P o = 1 − P c , with a difference in free energies between the two states ΔG = G c − G o . The free energy difference depends on the physical parameters of the model such as \({l}_{{\rm{e}}{\rm{f}}{\rm{f}}}\), L, ε and d (as it can be seen in Eq. (S19)), and how ΔG increases as a function of these parameter is briefly discussed in the Supplementary Material.
Model with mRNA circularisation. Sketches of the two possible mRNA conformations, open and circularised, whose transitions depend on the free energy gap ΔG (A). In blue and yellow we have represented the interacting proteins (e.g. PABP and eIF4F) bound at the 3′ and 5′ ends; the black line is the end-to-end distance that is equal to d in the circularised state. We have fixed d = 5 nm in our calculation. Ribosome density computed taking into account mRNA circularisation (continuous line) is then compared to experimental data (triangles, cfr. symbols used in Fig. 1) and the previous model neglecting circularisation (dashed line) (B). The fitted parameters are \({\bar{\alpha }}_{\infty }\mathrm{=(4.7}\pm \mathrm{0.6)}\,{10}^{-3}\) s−1, λ = 7.0 ± 0.6 nm and ε = − 8.3 ± 0.4 (k B T). End-to-end distance (blue curve, right axis) and calculated probabilities P c and P o = 1 − P c as a function of the CDS length L (C).
To find how the average end-to-end distance is affected by this interaction we weight its value in the circularised and open state with their respective probabilities:
$${R}_{{\rm{circ}}}={P}_{o}R+{P}_{c}d\,\mathrm{.}$$
If P o ≈ 1 then our feedback model well approximates the mRNA translation process. While in prokaryotes we can likely assume P o = 1, this is probably an oversimplification for eukaryotes. In order to consider transcript circularisation we have now to compute how P o depends on the CDS length L and on the interaction energy ε (in k B T units) between the two ends. The intuitive explanation we used before to determine differences in local concentrations of ribosomes close to the 5′ end can be reproduced here to compute P c . For a fixed ε and a short mRNA, we expect to find circularised transcripts with a probability P c close to one; in contrast, very large transcripts should be hardly found in the circularised state. This length dependence will also contribute to the ribosome density behaviour observed in experiments (Fig. 1). We computed P c as a function of L, \({l}_{{\rm{e}}{\rm{f}}{\rm{f}}}\) and ε, which turns out to be:
$${P}_{c}=\frac{1}{1+{e}^{{\rm{\Delta }}G/{k}_{B}T}}=\frac{1}{1+[{(\frac{{l}_{{\rm{e}}{\rm{f}}{\rm{f}}}L}{{d}^{2}})}^{\frac{3}{2}}\sqrt{\frac{4\pi }{3}}-1]\,{e}^{(\frac{2{\pi }^{2}{l}_{{\rm{e}}{\rm{f}}{\rm{f}}}}{L}+\varepsilon )}}\,.$$
The details of the calculation of P c can be found in the Supplementary Material. Here ε is considered as a fitting parameter. By inserting Eq. (6) in Eq. (5) and computing the end-to-end distance to be plugged in Eq. (4) we eventually find, now as a function of \({\bar{\alpha }}_{\infty }\), λ and ε, how the initiation rate is affected by the concentration increase of ribosomes in the 5′ reaction volume with also considering transcript circularisation. We have then fitted \(\rho (\bar{\alpha })\) to the dataset we have previously used, and the outcome is shown in Fig. 5B (fitting values and other datasets can be found in the Supplementary Material). We did not find a significative difference from the best fit of the simpler model previously introduced (dashed line in Fig. 5B). In Fig. 5C we show how the probabilities of finding a circularised or open mRNA depend on the transcript length, together with the end-to-end distance R. The results agree with our intuitive explanation.
Length-dependent competition for resources
In this section we speculate on some potential applications of our model related to bacterial growth laws38,39. Specifically, we will show that cells can adjust, based on a length-discriminatory mechanism, their relative expression of genes at different ribosome concentrations. This result from our model leads us to predict and theorise a new regulation mechanism for gene expression.
We observe that changes in ribosome densities (or in translation efficiency) induced by changes in the ribosomal pool are conditional on the mRNA length (Fig. 6A). Short transcripts are less affected by the amount of available ribosomes compared to long ones, suggesting a possible mechanism to regulate the relative protein production rate at different growth rates based on transcript length only. As a proof of principle, in Fig. 6B we plot the relative expression \(\eta ({L}_{1},{L}_{2})={J}_{{L}_{1}}/{J}_{{L}_{2}}\) between the translation efficiencies of transcripts with lengths L 1 and L 2. When ribosomes are limiting, short transcripts are relatively more translated than long ones.
Effects of competition for resources (ribosomes) on the protein production rate. The ribosome density depends on the overall ribosome concentration c ∞. We show the ribosome density as a function of the transcript length for different concentrations of available ribosomes c ∞ (A). The curve denoted with c ∞ in the legend is built starting from the same parameters of Fig. 3A. We change c ∞ as described in the legend for the other curves. Short transcripts are less affected by changes in c ∞, as we also show in (B) where we plot the relative expression of transcripts η (defined in the text) as a function of c ∞. We used transcript with three different lengths (here L = 100, 500 and 2500). According to these results, ribosomal proteins that are short should be relatively more expressed under high ribosome competition regimes compared to other types of proteins.
This behaviour can be intuitively explained. The main contribution to the initiation of long mRNAs is the concentration of free ribosomes c ∞, meaning that long mRNAs are largely influenced by changes in the ribosomal pool. Short transcripts instead can more efficiently take advantage of the length dependent contribution δ R to initiation.
This constitutes a potential mechanism for regulating the relative expression between short and long genes. Ribosomal proteins are composed only of a few dozens of amino-acids and, as a consequence, our theory predicts that at least qualitatively ribosomal proteins should be proportionately more efficiently translated than longer proteins under strong ribosome competition regimes, i.e. low growth rates. This difference should decrease when ribosomal resources get less tight, and in the limit of infinite ribosomal resources, η should tend to 1 since the length dependence becomes less and less relevant.
Our model suggests a translational mechanism to relatively over express short proteins (e.g. ribosomal proteins) at the cost of longer ones.
Many aspects of translation are still puzzling researchers. Theoretical approaches propose designing principles for modulating the translation efficiency at the level of initiation4 and of elongation24, mainly based on the role of RNA secondary structures, codon bias or amino acid properties. However, when tested on synthetic constructs, the hypotheses underlying the theoretical models are often contradicted8, so that the identification of transcript-dependent determinants of translation efficiency is still debated in the literature. In this work we have identified and studied another factor modulating the translation efficiency: the length of the transcript.
According to our results, the proximity of the 3′ end to the ribosomal recruitment site of the mRNA could induce a feedback in the translation process that would favour the translation of short transcripts over long ones, as it has been shown by experiments in the last decades. We connected the emergence of the ribosome density-mRNA length dependence shown in Fig. 1 to a mechanistic model built on basic physical principles and on the properties of the translation process. Our theory then establishes a link between densities and mRNA lengths and explains experimental data with an excellent agreement without invoking an evolutionary selection of genes based on their length. As a matter of fact, a selection process towards short efficient genes to improve cellular fitness could also be conjectured40,41. Without direct experimental observation we cannot rule out this hypothesis, although this would not explain the same behaviour observed for different organisms (experimental data in Fig. 1 seems to collapse to a unique universal curve). Moreover, the poor correlation between free energy of secondary structures at the 5′ end of an mRNA and CDS length is a signature that binding sites are not significantly weaker for long genes, as it would assume an evolutionary argument. Instead our theory predicts, as an outcome, that short genes have larger initiation rates compared to long ones. A recent publication by Li et al.42. confirms that transcript length is a main determinant of translation.
Experimental results in fact suggest that translation initiation is also dependent on mRNA length10,17. Since the ribosome recruitment rate must depend on the local concentration of ribosomal subunits around the 5′UTR, we conjecture that the local concentration is modulated by the CDS length via a feedback mechanism coupling protein synthesis and initiation rates. We can roughly name this process "recycling", as subunits terminating translation will contribute to the increase of the local concentration δ R and thus be more easily re-used as sketched in Fig. 2. This coarse-grained physical model has allowed us to reproduce the scaling behaviour of experimental ribosome densities (Fig. 3) and it constitutes, to our understating, the first quantitative explanation of how translational features are affected by the transcript length. Previous works have studied the effect of particle recycling43,44,45,46 but, with the exception of Chou (2003)43, they do not explicitly compute how the recycling term is regulated by the end-to-end distance R.
The compaction of the transcript (here characterised by the end-to-end distance), also depends on its polysome state. Intuitively, an mRNA with many translating ribosomes will be more stretched than an empty mRNA. We captured this feature by introducing a ribosome density dependence on the end-to-end distance through Eq. (3). We emphasise that this is the simplest choice for coupling elongation properties and the three-dimensional conformation of the mRNA, and one could introduce more complicated relationships linking end-to-end distance, ribosome density and elongation rate; here we wanted to show that, as a proof of principle, by a feedback mechanism enhancing initiation we can reproduce experimental data very well.
We have also studied how the results change by explicitly considering the diffusion of the two ribosomal subunits, and we found no significant change (see Supplementary Material). To make the model more realistic we considered two further extensions of the model. We considered (i) ribosome drop-off and (ii) mRNA circularisation. The former brought an improvement in the comparison between data and theory for large transcripts only (see Fig. 4), while the inclusion of the latter did not lead to a particular phenomenological change of the model's outcomes (Fig. 5). This suggests that ribosome recycling, as considered in the basic model, is the fundamental element originating the length-dependence translation.
We have then speculated on how the length could be exploited to create differences in the relative expression of genes at different growth rates, here used as measure of the free ribosomes abundances. When resources are constrained, i.e. when the amount of free ribosomes c ∞ is small, competition for resources might become relevant47,48,49 and our theory predicts that the length-dependent term of the initiation rate dominates the process. In other words, when ribosomes are strongly limiting, ribosome recruitment is mainly due to recycled ribosomes, meaning that short transcripts can better capitalise the resources. This mechanism could also be a way to translationally favour the production of ribosomal proteins (which are short) in a scenario of deficiency of ribosomes. In order to formulate this hypothesis, we neglect known mechanisms responsible for ribosome biogenesis, a complex process that is beyond the scope of our work. Our conjecture has then to be interpreted in the perspective of ribosomal concentrations fixed by a certain amount of ribosome production (established, for instance, by the richness of the growth medium): for a given concentration of ribosomes we make strong predictions on how the ribosomal pool should be partitioned among the different transcripts with just a length-discrimination mechanism.
The model could be further extended to consider translation of bacterial operons: in this case, in fact, one transcript is composed of different sinks (ribosome binding sites) and sources (stop codons) of ribosomes, while here we have discussed the case of a transcript translating a single gene (with one ribosome binding site and one stop codons). Having an operon translating different genes will increase the complexity of the feedback term, and it could in principle create counter-intuitive phenomenologies.
Our findings are compared to experimental ribosomal densities, and our framework can quantitatively reproduce the measurements. We have used our model to estimate the protein production rates of synonymous genes, and the method was successful (see Figure S13). In this work we have used our model to predict the ribosome-length dependence, i.e. we have emphasised the dependence of Eq. (4) on L, but our theory predicts that there is a feedback between elongation (codon usage) and initiation, that we have exploited in Figure S13. However, to further test the model and the relevance of the regulatory mechanism we propose, it will be necessary to make fusions of a reporter gene with peptides of variable lengths, and then measure ribosome density or translation efficiency with the aim of experimentally reproducing Fig. 3 in a controlled manner. However, one should pay attention to the changes in mRNA degradation, translation elongation and initiation induced by the added nucleotides coding the fused peptides. Figure 6B also constitutes a good way to test our hypotheses. One of our predictions is the relative change of expression of short/long transcripts when changing the cellular growth rate. This could be obtained by concurrently expressing two reporter genes of different lengths, and measuring their relative expression at different growth rates obtained by changing growth medium or by different antibiotics.
The exclusion process
We base our model on the exclusion process, which is also introduced in the section Results and in Fig. 2B. More accurately, this model is known in the literature as TASEP: Totally Asymmetric Simple Exclusion Process, for which nowadays there exists a plethora of extensions applied in many different fields, from vehicular traffic to intracellular transport19. Each site of the track in Fig. 2B corresponds to a codon, and the particles can advance from site to site -provided that the next site is not occupied by another particle- mimicking the elongation process. We give a thorough description of the exclusion process and the known results in the Supplementary Material.
To simulate the dynamics of the exclusion process we used a kinetic Gillespie-like Monte Carlo as used in Ciandrini et al.24.
Fitting and numerical solutions
We substitute the expression for \(\bar{\alpha }\), Eq. (4), in the equation for the density \(\rho (\bar{\alpha })\) in the low density phase of the \(\ell \)-TASEP (see Supplementary Material). The current J is given by the J(ρ) correction in the \(\ell \)-TASEP and the end-to-end distance is R found in Eq. (3). Thus, we obtain an implicit equation \(\rho =\rho ({\bar{\alpha }}_{\infty },\lambda ,L)\) that can be numerically solved for each set of variables \(\{{\bar{\alpha }}_{\infty },\lambda ,L\}\) and used to fit the experimental data ρ exp(L) to obtain the parameters \({\bar{\alpha }}_{\infty }\) and λ for each dataset, and their standard errors (see Supplementary Material).
We have used built-in functions of Mathematica50 to obtain numerical solutions for the density and currents and to fit the three datasets used in this study.
Density and current via a self-consistent simulation scheme
Equation (4) allows us to obtain, via simulation, the values for ρ for different values of L, taking into account the feedback mechanism coupling protein synthesis and initiation rate and finite size effects (the later intrinsic to the numerical simulations). We obtain this with the following self-consistent method:
We initialise the system with an arbitrary value of α = α (0), let the system evolve until the steady state is reached and then evaluate the current J (0) and the density ρ;(0)
Compute R as in equation (3) and update α = α (1) according to equation (4) with J (0) and ρ (0) computed in (i);
Repeat the previous points for several iterations until |α (i) − α (i − 1)|/α (i) < 0.01 (in general less than 10 iterations are needed to make the algorithm converge);
The final value of α is then used to obtain the final densities and currents.
With this iteration process, for a given choice of the parameters α ∞ and λ, we can obtain the steady state density and current, which vary with the length L of the transcript, due to the joint contribution of recycling and finite-size effects.
This self-consistent method allows us to simulate the system without explicitly considering, thanks to Equation (4), the diffusion of particles when they are not bound to the lattice and the dynamics of the mRNA.
Choice of datasets
We restricted our analysis to measures made by sucrose gradient methods. We are aware that a more recent technique like ribosome profiling12 would provide ribosome densities with codon resolution, but this method does not provide an absolute ribosome density (see definition below). The length-correlation has been shown to hold in ribosome profiling experiments12. Instead of assuming arbitrary normalisation of ribosome footprints to match our theory, we analysed absolute ribosome densities that are available in the literature. For the Mackay et al. dataset11 we used the reliable subset of data.
We define the ribosome density to be the number N of translating ribosomes divided by the length L (expressed in number of codons) of the CDS. We embrace this definition instead of alternative ones (ribosomes per 100 or 1000 nucleotides) for practical reasons. Thus, the density ρ ≡ N/L is expressed in ribosomes per codons, and it can be thought of as the probability of a codon being covered by the the A-site of the ribosome (i.e., a codon being translated). For steric reasons, this density is bound by 1/\(\ell \), where \(\ell \) is the length of the ribosome footprint (in codons). For instance, \(\ell \sim 10\) in S. cerevisiae.
Gingold, H. & Pilpel, Y. Determinants of Translation Efficiency and Accuracy. Molecular Systems Biology 7, 481, https://doi.org/10.1038/msb.2011.14 (2011).
Quax, T. E. F., Claassens, N. J., Söll, D. & van der Oost, J. Codon Bias as a Means to Fine-Tune Gene Expression. Molecular Cell 59, 149–161, https://doi.org/10.1016/j.molcel.2015.05.035 (2015).
Kudla, G., Murray, A. W., Tollervey, D. & Plotkin, J. B. Coding-Sequence Determinants of Gene Expression in Escherichia coli. Science 324, 255–258, https://doi.org/10.1126/science.1170160 (2009).
ADS CAS Article PubMed PubMed Central Google Scholar
Salis, H. M., Mirsky, E. A. & Voigt, C. A. Automated design of synthetic ribosome binding sites to control protein expression. Nat Biotech 27, 946–950, https://doi.org/10.1038/nbt.1568 (2009).
Plotkin, J. B. & Kudla, G. Synonymous but not the same: the causes and consequences of codon bias. Nat Rev Genet 12, 32–42, https://doi.org/10.1038/nrg2899 (2011).
Kemp, A. J. et al. A yeast tRNA mutant that causes pseudohyphal growth exhibits reduced rates of CAG codon translation. Molecular Microbiology 87, 284–300, https://doi.org/10.1111/mmi.12096 (2013).
Gorgoni, B., Ciandrini, L., McFarland, M. R., Romano, M. C. & Stansfield, I. Identification of the mRNA targets of tRNA-specific regulation using genome-wide simulation of translation. Nucl. Acids Res. gkw630, https://doi.org/10.1093/nar/gkw630 (2016).
Li, G.-W. How do bacteria tune translation efficiency? Current Opinion in Microbiology 24, 66–71, https://doi.org/10.1016/j.mib.2015.01.001 (2015).
Picard, F. et al. Bacterial translational regulations: high diversity between all mRNAs and major role in gene expression. BMC Genomics 13, 528, https://doi.org/10.1186/1471-2164-13-528 (2012).
Arava, Y. et al. Genome-wide analysis of mRNA translation profiles in Saccharomyces cerevisiae. PNAS 100, 3889–3894, https://doi.org/10.1073/pnas.0635171100 (2003).
MacKay, V. L. et al. Gene expression analyzed by high-resolution state array analysis and quantitative proteomics: response of yeast to mating pheromone. Mol. Cell Proteomics 3, 478–489, https://doi.org/10.1074/mcp.M300129-MCP200 (2004).
Ingolia, N. T., Ghaemmaghami, S., Newman, J. R. S. & Weissman, J. S. Genome-wide analysis in vivo of translation with nucleotide resolution using ribosome profiling. Science 324, 218–223, https://doi.org/10.1126/science.1168978 (2009).
Lacsina, J. R., LaMonte, G., Nicchitta, C. V. & Chi, J.-T. Polysome profiling of the malaria parasite Plasmodium falciparum. Mol. Biochem. Parasitol. 179, 42–46, https://doi.org/10.1016/j.molbiopara.2011.05.003 (2011).
Cataldo, L., Mastrangelo, M. A. & Kleene, K. C. A quantitative sucrose gradient analysis of the translational activity of 18 mRNA species in testes from adult mice. Mol. Hum. Reprod. 5, 206–213 (1999).
Hendrickson, D. G. et al. Concordant Regulation of Translation and mRNA Abundance for Hundreds of Targets of a Human microRNA. PLOS Biol 7, e1000238, https://doi.org/10.1371/journal.pbio.1000238 (2009).
Guo, J., Lian, X., Zhong, J., Wang, T. & Zhang, G. Length-dependent translation initiation benefits the functional proteome of human cells. Mol Biosyst 11, 370–378, https://doi.org/10.1039/c4mb00462k (2015).
Arava, Y. Compaction of polyribosomal mRNA. RNA Biol 6, 399–401 (2009).
MacDonald, C. T., Gibbs, J. H. & Pipkin, A. C. Kinetics of biopolymerization on nucleic acid templates. Biopolymers 6, 1–5, https://doi.org/10.1002/bip.1968.360060102 (1968).
Chou, T., Mallick, K. & Zia, R. K. P. Non-equilibrium statistical mechanics: from a paradigmatic model to biological transport. Rep. Prog. Phys. 74, 116601, https://doi.org/10.1088/0034-4885/74/11/116601 (2011).
ADS MathSciNet Article Google Scholar
Blythe, R. A. & Evans, M. R. Nonequilibrium steady states of matrix-product form: a solver's guide. J. Phys. A: Math. Theor. 40, R333, https://doi.org/10.1088/1751-8113/40/46/R01 (2007).
ADS MathSciNet Article MATH Google Scholar
Mitarai, N., Sneppen, K. & Pedersen, S. Ribosome Collisions and Translation Efficiency: Optimization by Codon Usage and mRNA Destabilization. Journal of Molecular Biology 382, 236–245, https://doi.org/10.1016/j.jmb.2008.06.068 (2008).
Brackley, C. A., Romano, M. C. & Thiel, M. The Dynamics of Supply and Demand in mRNA Translation. PLOS Comput Biol 7, e1002203, https://doi.org/10.1371/journal.pcbi.1002203 (2011).
ADS MathSciNet CAS Article PubMed PubMed Central Google Scholar
Zia, R. K. P., Dong, J. J. & Schmittmann, B. Modeling translation in protein synthesis with TASEP: A tutorial and recent developments. J Stat Phys 144, 405, https://doi.org/10.1007/s10955-011-0183-1 (2011).
ADS MathSciNet CAS Article MATH Google Scholar
Ciandrini, L., Stansfield, I. & Romano, M. C. Ribosome traffic on mRNAs maps to gene ontology: genome-wide quantification of translation initiation rates and polysome size regulation. PLoS Comput. Biol. 9, e1002866, https://doi.org/10.1371/journal.pcbi.1002866 (2013).
Reuveni, S., Meilijson, I., Kupiec, M., Ruppin, E. & Tuller, T. Genome-Scale Analysis of Translation Elongation with a Ribosome Flow Model. PLOS Comput Biol 7, e1002127, https://doi.org/10.1371/journal.pcbi.1002127 (2011).
Tian, C. et al. Rapid Curtailing of the Stringent Response by Toxin-Antitoxin Encoded mRNases. J. Bacteriol. JB.00062–16, https://doi.org/10.1128/JB.00062-16 (2016).
Shaw, L. B., Zia, R. K. P. & Lee, K. H. Totally asymmetric exclusion process with extended objects: A model for protein synthesis. Phys. Rev. E 68, 021910, https://doi.org/10.1103/PhysRevE.68.021910 (2003).
ADS Article Google Scholar
Siwiak, M. & Zielenkiewicz, P. A Comprehensive, Quantitative, and Genome-Wide Model of Translation. PLOS Comput Biol 6, e1000865, https://doi.org/10.1371/journal.pcbi.1000865 (2010).
ADS Article PubMed PubMed Central Google Scholar
Vanzi, F., Takagi, Y., Shuman, H., Cooperman, B. S. & Goldman, Y. E. Mechanical Studies of Single Ribosome/mRNA Complexes. Biophysical Journal 89, 1909–1919, https://doi.org/10.1529/biophysj.104.056283 (2005).
Bakshi, S., Siryaporn, A., Goulian, M. & Weisshaar, J. C. Superresolution imaging of ribosomes and rna polymerase in live escherichia coli cells. Molecular Microbiology 85, 21–38, https://doi.org/10.1111/j.1365-2958.2012.08081.x (2012).
Bremer, H. & Dennis, P. P. Modulation of Chemical Composition and Other Parameters of the Cell by Growth Rate. Escherichia coli and Salmonella: cellular and molecular biology 2, 1527–1542, https://doi.org/10.1046/j.1365-2486.2000.06017.x (1987).
Arkin, A., Ross, J. & McAdams, H. H. Stochastic kinetic analysis of developmental pathway bifurcation in phage λ infected escherichia coli cells. Genetics 149, 1633–1648 (1998).
Ringnér, M. & Krogh, M. Folding free energies of 5′-utrs impact post-transcriptional regulation on a genomic scale in yeast. PLOS Computational Biology 1, 1–8, https://doi.org/10.1371/journal.pcbi.0010072 (2005).
Heyer, E. E. & Moore, M. J. Redefining the Translational Status of 80 s Monosomes. Cell 164, 757–769, https://doi.org/10.1016/j.cell.2016.01.003 (2016).
Sin, C., Chiarugi, D. & Valleriani, A. Quantitative assessment of ribosome drop-off in E. coli. Nucl. Acids Res. 44, 2528–2537, https://doi.org/10.1093/nar/gkw137 (2016).
Bonnin, P., Kern, N., Young, N. T., Stansfield, I. & Romano, M. C. Novel mrna-specific effects of ribosome drop-off on translation rate and polysome profile. PLOS Computational Biology 13, 1–38, https://doi.org/10.1371/journal.pcbi.1005555 (2017).
Wells, S. E., Hillner, P. E., Vale, R. D. & Sachs, A. B. Circularization of mrna by eukaryotic translation initiation factors. Molecular Cell 2, 135–140, https://doi.org/10.1016/S1097-2765(00)80122-7 (1998).
Scott, M., Gunderson, C. W., Mateescu, E. M., Zhang, Z. & Hwa, T. Interdependence of Cell Growth and Gene Expression: Origins and Consequences. Science 330, 1099–1102, https://doi.org/10.1126/science.1192588 (2010).
ADS CAS Article PubMed Google Scholar
Dai, X. et al. Reduction of translating ribosomes enables Escherichia coli to maintain elongation rates during slow growth. Nature Microbiology 2, 16231, https://doi.org/10.1038/nmicrobiol.2016.231 (2016).
Castillo-Davis, C. I., Mekhedov, S. L., Hartl, D. L., Koonin, E. V. & Kondrashov, F. A. Selection for short introns in highly expressed genes. Nat Genet 31, 415–418 (2002).
Eisenberg, E. & Levanon, E. Y. Human housekeeping genes are compact. Trends in Genetics 19, 362–365, https://doi.org/10.1016/S0168-9525(03)00140-9 (2003).
Li, J. J., Chew, G.-L. & Biggin, M. D. Quantitating translational control: mRNA abundance-dependent and independent contributions and the mRNA sequences that specify them. Nucleic Acids Research gkx898 (2017). https://doi.org/10.1093/nar/gkx898.
Chou, T. Ribosome Recycling, Diffusion, and mRNA Loop Formation in Translational Regulation. Biophysical Journal 85, 755–773, https://doi.org/10.1016/S0006-3495(03)74518-4 (2003).
Sharma, A. K. & Chowdhury, D. Stochastic theory of protein synthesis and polysome: Ribosome profile on a single mRNA transcript. Journal of Theoretical Biology 289, 36–46, https://doi.org/10.1016/j.jtbi.2011.08.023 (2011).
Margaliot, M. & Tuller, T. Ribosome flow model with positive feedback. Journal of the Royal Society Interface 10, 20130267, https://doi.org/10.1098/rsif.2013.0267 (2013).
Marshall, E., Stansfield, I. & Romano, M. C. Ribosome recycling induces optimal translation rate at low ribosomal availability. J R Soc Interface 11, 20140589, https://doi.org/10.1098/rsif.2014.0589 (2014).
Greulich, P., Ciandrini, L., Allen, R. J. & Romano, M. C. Mixed population of competing totally asymmetric simple exclusion processes with a shared reservoir of particles. Phys. Rev. E 85, 011142, https://doi.org/10.1103/PhysRevE.85.011142 (2012).
Shah, P., Ding, Y., Niemczyk, M., Kudla, G. & Plotkin, J. B. Rate-limiting steps in yeast protein translation. Cell 153, 1589–1601, https://doi.org/10.1016/j.cell.2013.05.049 (2013).
Raveh, A., Margaliot, M., Sontag, E. D. & Tuller, T. A model for competition for ribosomes in the cell. Journal of the Royal Society Interface 13, 20151062, https://doi.org/10.1098/rsif.2015.1062 (2016).
Mathematica. Wolfram research, inc. Version 10.2 (2015).
The authors would like to acknowledge the funding provided by the European Union Seventh Framework Programme [FP7/2007–2013] (NICHE; grant agreement 289384) (LDF). LDF also acknowledges the funding provided by the São Paulo Research Foundation (FAPESP - grant #2015/26989-4). AM was partially funded by the UK Biotechnology and Biological Research Council (BBSRC), through grant BB/N015711/1. LC would like to acknowledge Maria Carmen Romano, Jean Hausser, Marco Cosentino Lagomarsino, Jean-Charles Walter and Norbert Kern for early discussions on this work, and the CNRS for having granted him a "demi-délégation" (2017–18). We would like to dedicate this work in memory of Maxime Clusel and Vladimir Lorman.
Departamento de Entomologia e Acarologia, Escola Superior de Agricultura Luiz de Queiroz - Universidade de São Paulo (USP), 13418-900, Piracicaba/SP, Brazil
Lucas D. Fernandes
Institute for Complex Systems and Mathematical Biology, University of Aberdeen, Aberdeen, AB24 3UE, UK
Lucas D. Fernandes & Alessandro P. S. de Moura
DIMNP UMR 5235, Université de Montpellier and CNRS, F-34095, Montpellier, France
Luca Ciandrini
Laboratoire Charles Coulomb UMR5221, Université de Montpellier and CNRS, F-34095, Montpellier, France
Alessandro P. S. de Moura
L.C. and A.M. conceived and discussed the theory, L.C. and L.D.F. performed the simulations and analysed the data, L.C., L.D.F. and A.M. wrote the manuscript text and supplementary information. All authors reviewed the manuscript.
Correspondence to Luca Ciandrini.
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
LaTeX Supplementary File
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
Fernandes, L.D., Moura, A.P.S.d. & Ciandrini, L. Gene length as a regulator for ribosome recruitment and protein synthesis: theoretical insights. Sci Rep 7, 17409 (2017). https://doi.org/10.1038/s41598-017-17618-1
Monosomes actively translate synaptic mRNAs in neuronal processes
Anne Biever
, Caspar Glock
, Georgi Tushev
, Elena Ciirdaeva
, Tamas Dalmay
, Julian D. Langer
& Erin M. Schuman
Pre-emptive Quality Control of a Misfolded Membrane Protein by Ribosome-Driven Effects
Ramya Lakshminarayan
, Ben P. Phillips
, Imogen L. Binnian
, Natalia Gomez-Navarro
, Norberto Escudero-Urquijo
, Alan J. Warren
& Elizabeth A. Miller
Current Biology (2020)
So close, no matter how far: multiple paths connecting transcription to mRNA translation in eukaryotes
Boris Slobodin
& Rivka Dikstein
EMBO reports (2020)
A chemical kinetic basis for measuring translation initiation and elongation rates from ribosome profiling data
Ajeet K. Sharma
, Pietro Sormanni
, Nabeel Ahmed
, Prajwal Ciryam
, Ulrike A. Friedrich
, Günter Kramer
, Edward P. O'Brien
& Tamir Tuller
PLOS Computational Biology (2019)
Quantitative principles of cis-translational control by general mRNA sequence features in eukaryotes
Jingyi Jessica Li
, Guo-Liang Chew
& Mark Douglas Biggin
Genome Biology (2019)
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
About Scientific Reports
Guest Edited Collections
Scientific Reports Top 100 2017
Scientific Reports Top 10 2018
Editorial Board Highlights
Author Highlights
Scientific Reports ISSN 2045-2322 (online)
Close banner Close
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.
I agree my information will be processed in accordance with the Nature and Springer Nature Limited Privacy Policy.
Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing


Raymond is a Managing Partner at the Latino Business Speakers Bureau and also a Managing Partner at Chapman Farrell Group and former Chief Growth Officer for ALPFA and President of ALPFA Solutions, delivering leadership development programs for the growing Latino workforce. Raymond, a former Chief Diversity Officer at Aetna, is a Fortune 100 bilingual leader with more than 25 years of experience assessing talent needs in consumer goods, professional and financial services, association/nonprofit, and healthcare sectors. Raymond earned his B.S. and M.S. degrees from New York University and a governance certificate from Harvard University. As a volunteer, Raymond served as Chairman of the Board of Directors of the U.S. Hispanic Chamber of Commerce, Officer of the Board of Trustees of the Wheeler Clinic, and as a member of the Board of Directors of the National Hispanic Medical Association. He has served on ALPFA's Corporate Advisory Board, as President of PRIMER, as Chairman of the Corporate and Philanthropic Council of HACU, and as a member of the Board of Regents of the University of Hartford. Raymond was selected two years in a row by Hispanic Business Magazine as one of the 100 Most Influential Hispanics; by the National Society of Hispanic MBAs for the Corporate Executive Excellence Brillante Award; by Pricewaterhousecoopers as one of its exemplary alumnus and featured on its KeyWord Magazine; by Urban Professional Magazine as one of the 100 Top Executives in Corporate America, and was the featured cover story in Diversity Executive Magazine.


text:
Grey's Anatomy 15x04 Review: "Momma Knows Best" (Running Blind) [Contributor: Julia Siegel]
8:39 PM Contributor: Julia Siegel, grey's anatomy recap, s15.04 No comments
"Momma Knows Best"
Original Airdate: October 11, 2018
Calling all Grey's Anatomy fans! This is the episode you have been waiting for, as Meredith goes on her first blind date and some of the Station 19 crew visit Grey Sloan Memorial. While this isn't the best crossover, there are plenty of juicy moments that only Shonda Rhimes' shows deliver. Let's just say Teddy and Jackson are going to need a refresher when they eventually come back.
IT'S TED MOSBY!
There's only one obvious place to start with this episode — and no, the main plot isn't as important as Meredith's dating life. Patient CeCe has set up Meredith's first blind date and has only given her a name and a photo of her match. Meredith is finally happy to comply with the process and shows up at work dressed up for her date that afternoon. Naturally, she catches the eye of every character on the show, but the doctors who have been vying for her attention have the best reactions. Alex is the only person to make fun of Meredith's date, so they make a bet on how long the date will last. Both think she will hate the guy and only stay for twenty minutes; so maybe Meredith isn't as open-minded as we all hoped.
Upon arriving to the restaurant for her date, Meredith finds the guy she thinks is her date — John, played by the always wonderful Josh Radnor. While John might not be a hopeless romantic architect living in New York City, he is just as charming as Radnor's most famous character, Ted Mosby of How I Met Your Mother.
John is divorced, works in privacy technology, and is looking to find someone to fill the void in his life. It was a little surprising to see how fast he and Meredith hit it off, but their similar senses of humor and lifestyles give Meredith a glimmer of hope. They see another couple across the restaurant and start making fun of them, only to realize that the other couple is actually the two dates they were supposed to meet. As Ted Mosby would say, fate intervened and brought John and Meredith together.
They decide to text their initial dates and say that it isn't going to work out to meet up. Meredith prolongs the date and goes for a walk with John. Everything is going smoothly until John mentions that he doesn't like dating single moms because they always act desperate on dates and want to move too quickly. Of course, he had no idea that Meredith has three small children, but it's too late. When she tells him that she has kids and isn't like that, John instantly realizes he said the one thing that would stop their potential romance in its tracks. Meredith says goodbye and leaves, which is sad considering how well they were getting along.
However, I have a feeling this isn't the last we will see of John. The character was very into Meredith, and it didn't seem like he wouldn't back down that easily. Also, it would be odd for someone like Josh Radnor to only be cast for one episode. On the guest starring list at the beginning of the episode, he was listed last as "and Josh Radnor." So, like Ted Mosby, I will believe that there could be a real connection between Meredith and John, leading to their eventual reunion.
When two shows have a crossover night, I expect that more than three characters from one show will make an appearance on the other. The whole crossover thing was pretty lame this time around, since only Ben, Andy, and Dean from Station 19 make brief appearances in the Grey's Anatomy hour. Maggie and DeLuca are the only Grey's characters that extend into the Station 19 hour, so it didn't feel like a real crossover night. Instead of the characters moving between shows, it was more of both episodes sharing the same main plot. Ben, Andy, and Dean bring victims of a fire to Grey Sloan Memorial for medical care. Bailey sees Ben passing in the hallway and notices that he has a bandage on his left hand. He tells her that he burned his hand on a doorknob and promptly leaves. Next up, Andy and Dean's patient is a mother suffering from a problem with her heart. Maggie and DeLuca take over the case while Andy and Dean stay with the patient's young son.
The son wants to know everything that is happening to his mom; but when his dad arrives, he puts the kibosh on the doctors telling his kid anything. Maggie has to deliver the bad news to the father that his wife isn't going to survive and struggles with the fact that he won't tell his son that his mother is dying. Poor Maggie doesn't need any more grief on her plate, but at least she gets a message from Jackson that says he will be coming back soon. She gets some much-needed laughter when Dean hangs around the hospital just to ask her out. In typical Maggie fashion, she doesn't respond right away.
Maggie is pretty much at her limit for secrets and stress by this point. After talking to Meredith, Amelia, and Andy — who I guess is an honorary member of the sisterhood now — about the problems with the patient, Maggie decides that she can't withhold the truth anymore. She eventually blabs to Meredith that Teddy is pregnant with Owen's baby and that no one knows. I'm surprised it took Maggie this long to confess her secret, considering the immediate burden it put on her. On the dying patient front, Maggie brings the husband to his wife to say goodbye, but the man won't let his son come along. DeLuca takes the boy to the roof to talk to him and pretty much says that his mother is dying. The frightened boy flees the hospital, setting up the next hour of drama on Station 19.
Don't worry — after falling into Seattle's underground water system, the boy is eventually rescued by our favorite firefighters and learns the truth about his mom from his father. Maggie eventually turns Dean's offer down, but he leaves it as an open invitation in case she ever finds herself single.
Personally, I think Maggie should actually go out on a date with Dean. At least it would be a more stable relationship than Jackson can provide at the moment. Also, DeLuca joins the crew of the aid car to help the kid and gets his own nickname from Victoria. DeLuca mentions that he used to be a paramedic, so the stage has been set for more potential crossovers with him.
HAIL TO THE CHIEF
An overwhelmed Alex decides to take a day out of his new office to be a doctor again and treats another victim of the fire who is suffering from smoke inhalation. Young Julius doesn't want to stay in the hospital because his mother's insurance doesn't cover much. Alex and intern Qadri notice that Julius has a fixable lung condition. Julius rejects their proposed surgery because it is a previous condition uncovered by the health insurance. Alex essentially tells the kid he is being dumb by not having it fixed because it could harm him in the long run.
Richard steps in to help, but doesn't want to be the rock that Alex needs for the chief of surgery job. Richard tells Alex that he would have to be creative to treat his patient, but Alex takes it too far. He decides to cut the patient with a scalpel and pretend that he has a chest wound that needs surgical correction immediately. Alex commits insurance fraud to help his patient, much to the chagrin of Richard and Qadri. However, Alex feels that he can do whatever he wants because he is the chief.
Those famous last words will surely bite him in the butt sooner rather than later.
PARENTING WOES
The episode also continues the blossoming storyline of Owen and Amelia's young family. All is well with the on-and-off again couple, who don't want to put a label on their newly renewed status. Owen and Amelia provide the comedic relief of the episode in the form of their parenting strategies. The funny thing is that Amelia doesn't even realize that she is being an overbearing mother to teen Betty. Amelia gets upset when Betty comes home presumably high. Like the good parent she is, Amelia freaks out and yells about how a recovering addict can't do any drugs.
Naturally, Betty denies doing drugs and refuses to talk to Amelia when she is being irrational. Amelia drags Betty to work with her and Owen to force her into taking a drug test (which doesn't go over very well). Owen spends most of the episode smiling and shaking his head at Amelia's antics.
Since Amelia will stop at nothing to protect Betty, she calls in reinforcements in the form of Owen's mom. Evelyn Hunt and Amelia have a good conversation that makes Amelia realize she is a parent and needs to handle the situation appropriately so Betty will grow and learn from it. Amelia eventually gets Betty to confess that she smoked some weed, which leads to Amelia grounding her for one month. Owen laughs while Betty seethes, and the Hunt-Shepherd family continues on! With a bombshell of a secret looming, let's hope they enjoy their honeymoon phase while it lasts.

redpajama_set_name: RedPajamaCommonCrawl


This wonderful flat offers two spacious bedrooms , master with private bath, walk in closet, bright and spacious living room with picture window and high ceilings. The well-appointed kitchen with breakfast bar overlooking the living room with hardwood floors. Fabulous Colchester location just minutes to the I-89 and short drive to Burlington, UVM, shopping and airport.


Are you the talent we're looking for?
At Autarco we're looking for the best of the best to join our international team. Below are the positions we currently have available, but make sure to contact us if you're that exceptional talent we can't ignore.


Through the creation of architectural projects scaling from interior to private housing and smaller institutions, ReVærk aims to accentuate the relation between architectural space and materiality/tactility. The materials surrounding us and how these are used have great influence on our lives, physically as well as mentally, and it is this realization that has led the company to its specific dedication. It is ReVærk's intention to reintroduce natural materials and building techniques with unexplored relevance to the building industry of today. The vision is to procure the qualities of these materials and their crafts and through mindful unfolding to enrich the spaces they form for both greater aesthetic and environmental health for their occupants.


1. An Eye-Popping Inside Outline – with legitimate Airbnb Inclination. Having a remarkable outline can come in numerous styles. However, the primary concern is whether we are going for Deluxe, Current, Moderate, Bohemian, or Conventional, the outline must be better than what others are putting forth on Airbnb, or you will lose the business.
When all is said in done, the place must be professionally planned. Our organization has proficient architects. Also you need to outline , considering Airbnb, so the plan isn't the same as you would use for example as strangers do when offering a home. Be that as it may, as a rule, we are going for "eye popping".
2. Include Creation – There are "19 unique Handles" or settings in Airbnb as we call it, and if half of them are set wrong, you will really slice your deals down the middle. We can enable you to get this right. At that point there are around 3 different destinations, that you have to utilize, other than Airbnb , or you will miss out.
One tip here is to have an end of the week rate set versus an end of the week rate. You just have to call inns and imagine I am a long haul stayer at an inn, and request their valuing for the following a half year. You will see their valuing goes here and there relying upon occasions that come into town, so you will know how to value somewhat better here.
3. Speed of Promotion – the sooner the makeover, and launching of your adds – the more cash you make. Getting up and running as quick as two weeks or even one week, implies the distinction of thousands of dollars.
Most homes we do in Hollywood for instance, do $20K in addition to a month in Airbnb, so one week has a significant effect, obviously without yielding quality.
4. Deals/Appointments Group – individuals ask all day and all night, seven days seven days. Chinese individuals ask after 11 pm for instance in light of their time allotment.
5. Activities – tasks is the everyday activities and running of your Airbnb.
Everything from managing clients, managing accelerations, getting great audits (and consulting out of terrible ones), managing robberies and security stores, property support and crises, and Substantially more.
Finding a solid cleaning staff, that can likewise stock your provisions, and tally your stock, or welcome your visitors, is an absolute necessity.


Our event management team is ready to bring our world-class services to your conference.
We are innovative, we work fast, and we have the know-how, resources and experience to manage first-class events. We have extensive experience working with governments, associations, committees and the private sector. Perhaps most importantly, we only promise what we can deliver – many years' experience in major conferences and special events, excellent contacts, credentials and reputation, financial expertise and responsibility, and the ability to deliver an outstanding event every time.
As a full-service conference organizing and management provider, we can organize all aspects of a conference, including logistics management, scientific programming, financial services, registration, and on-site management and we can tailor services as required. We believe in a collaborative approach that honours and values the perspectives of multiple stakeholders. We have extensive experience working with a broad range of stakeholders in planning national knowledge exchange opportunities that address the learning needs of an equally diverse audience.
We focus on effective collaboration with the organizing committee(s) and sub-committees to ensure that the conference objectives and timelines are met.
Extra effort is always made to ensure we are communicating clearly with committees, suppliers, delegates, and everyone involved with the event.
Work with suppliers and vendors to develop innovative, creative and efficient ways to make the event a success.
We are committed to being transparent throughout the planning, budgeting, implementation and evaluation processes.
Our mutual success is built on solid relationships with everyone involved with event planning.