A kinetic model of protein synthesis. Application to hemoglobin synthesis and translational control

Comprehensive translational profiling and STE AI uncover rapid control of protein biosynthesis during cell stress

Translational control is important in all life, but it remains a challenge to accurately quantify. When ribosomes translate messenger (m)RNA into proteins, they attach to the mRNA in series, forming poly(ribo)somes, and can co-localize. Here, we computationally model new types of co-localized ribosomal complexes on mRNA and identify them using enhanced translation complex profile sequencing (eTCP-seq) based on rapid in vivo crosslinking. We detect long disome footprints outside regions of non-random elongation stalls and show these are linked to translation initiation and protein biosynthesis rates. We subject footprints of disomes and other translation complexes to artificial intelligence (AI) analysis and construct a new, accurate and self-normalized measure of translation, termed stochastic translation efficiency (STE). We then apply STE to investigate rapid changes to mRNA translation in yeast undergoing glucose depletion. Importantly, we show that, well beyond tagging elongation stalls, footprints of co-localized ribosomes provide rich insight into translational mechanisms, polysome dynamics and topology. STE AI ranks cellular mRNAs by absolute translation rates under given conditions, can assist in identifying its control elements and will facilitate the development of next-generation synthetic biology designs and mRNA-based therapeutics.

High-throughput and proteome-wide discovery of endogenous biomolecular condensates

Phase separation inside mammalian cells regulates the formation of the biomolecular condensates that are related to gene expression, signalling, development and disease. However, a large population of endogenous condensates and their candidate phase-separating proteins have yet to be discovered in a quantitative and high-throughput manner. Here we demonstrate that endogenously expressed biomolecular condensates can be identified across a cell's proteome by sorting proteins across varying oligomeric states. We employ volumetric compression to modulate the concentrations of intracellular proteins and the degree of crowdedness, which are physical regulators of cellular biomolecular condensates. The changes in degree of the partition of proteins into condensates or phase separation led to varying oligomeric states of the proteins, which can be detected by coupling density gradient ultracentrifugation and quantitative mass spectrometry. In total, we identified 1,518 endogenous condensate proteins, of which 538 have not been reported before. Furthermore, we demonstrate that our strategy can identify condensate proteins that respond to specific biological processes.

First-principles model of optimal translation factors stoichiometry

Enzymatic pathways have evolved uniquely preferred protein expression stoichiometry in living cells, but our ability to predict the optimal abundances from basic properties remains underdeveloped. Here, we report a biophysical, first-principles model of growth optimization for core mRNA translation, a multi-enzyme system that involves proteins with a broadly conserved stoichiometry spanning two orders of magnitude. We show that predictions from maximization of ribosome usage in a parsimonious flux model constrained by proteome allocation agree with the conserved ratios of translation factors. The analytical solutions, without free parameters, provide an interpretable framework for the observed hierarchy of expression levels based on simple biophysical properties, such as diffusion constants and protein sizes. Our results provide an intuitive and quantitative understanding for the construction of a central process of life, as well as a path toward rational design of pathway-specific enzyme expression stoichiometry.

Spurious regulatory connections dictate the expression-fitness landscape of translation factors

During steady-state cell growth, individual enzymatic fluxes can be directly inferred from growth rate by mass conservation, but the inverse problem remains unsolved. Perturbing the flux and expression of a single enzyme could have pleiotropic effects that may or may not dominate the impact on cell fitness. Here, we quantitatively dissect the molecular and global responses to varied expression of translation termination factors (peptide release factors, RFs) in the bacterium Bacillus subtilis. While endogenous RF expression maximizes proliferation, deviations in expression lead to unexpected distal regulatory responses that dictate fitness reduction. Molecularly, RF depletion causes expression imbalance at specific operons, which activates master regulators and detrimentally overrides the transcriptome. Through these spurious connections, RF abundances are thus entrenched by focal points within the regulatory network, in one case located at a single stop codon. Such regulatory entrenchment suggests that predictive bottom-up models of expression-fitness landscapes will require near-exhaustive characterization of parts.

Ligand-dependent tRNA processing by a rationally designed RNase P riboswitch

We describe a synthetic riboswitch element that implements a regulatory principle which directly addresses an essential tRNA maturation step. Constructed using a rational in silico design approach, this riboswitch regulates RNase P-catalyzed tRNA 5′-processing by either sequestering or exposing the single-stranded 5′-leader region of the tRNA precursor in response to a ligand. A single base pair in the 5′-leader defines the regulatory potential of the riboswitch both in vitro and in vivo. Our data provide proof for prior postulates on the importance of the structure of the leader region for tRNA maturation. We demonstrate that computational predictions of ligand-dependent structural rearrangements can address individual maturation steps of stable non-coding RNAs, thus making them amenable as promising target for regulatory devices that can be used as functional building blocks in synthetic biology.

Computational resources and strategies to assess single-molecule dynamics of the translation process in S. cerevisiae

This work provides a systematic and comprehensive overview of available resources for the molecular-scale modelling of the translation process through agent-based modelling. The case study is the translation in Saccharomyces cerevisiae, one of the most studied yeasts. The data curation workflow encompassed structural information about the yeast (i.e. the simulation environment), and the proteins, ribonucleic acids and other types of molecules involved in the process (i.e. the agents). Moreover, it covers the main process events, such as diffusion (i.e. motion of molecules in the environment) and collision efficiency (i.e. interaction between molecules). Data previously determined by wet-lab techniques were preferred, resorting to computational predictions/extrapolations only when strictly necessary. The computational modelling of the translation processes is of added industrial interest, since it may bring forward knowledge on how to control such phenomena and enhance the production of proteins of interest in a faster and more efficient manner.

Predictive biophysical modeling and understanding of the dynamics of mRNA translation and its evolution

mRNA translation is the fundamental process of decoding the information encoded in mRNA molecules by the ribosome for the synthesis of proteins. The centrality of this process in various biomedical disciplines such as cell biology, evolution and biotechnology, encouraged the development of dozens of mathematical and computational models of translation in recent years. These models aimed at capturing various biophysical aspects of the process. The objective of this review is to survey these models, focusing on those based and/or validated on real large-scale genomic data. We consider aspects such as the complexity of the models, the biophysical aspects they regard and the predictions they may provide. Furthermore, we survey the central systems biology discoveries reported on their basis. This review demonstrates the fundamental advantages of employing computational biophysical translation models in general, and discusses the relative advantages of the different approaches and the challenges in the field.

eIF4F: a retrospective

The original purification of the heterotrimeric eIF4F was published over 30 years ago (Grifo, J. A., Tahara, S. M., Morgan, M. A., Shatkin, A. J., and Merrick, W. C. (1983) J. Biol. Chem. 258, 5804-5810). Since that time, numerous studies have been performed with the three proteins specifically required for the translation initiation of natural mRNAs, eIF4A, eIF4B, and eIF4F. These have involved enzymatic and structural studies of the proteins and a number of site-directed mutagenesis studies. The regulation of translation exhibited through the mammalian target of rapamycin (mTOR) pathway is predominately seen as the phosphorylation of 4E-BP, an inhibitor of protein synthesis that functions by binding to the cap binding subunit of eIF4F (eIF4E). A hypothesis that requires the disassembly of eIF4F during translation initiation to yield free subunits (eIF4A, eIF4E, and eIF4G) is presented.

Control of ribosome traffic by position-dependent choice of synonymous codons

Messenger RNA (mRNA) encodes a sequence of amino acids by using codons. For most amino acids, there are multiple synonymous codons that can encode the amino acid. The translation speed can vary from one codon to another, thus there is room for changing the ribosome speed while keeping the amino acid sequence and hence the resulting protein. Recently, it has been noticed that the choice of the synonymous codon, via the resulting distribution of slow- and fast-translated codons, affects not only on the average speed of one ribosome translating the mRNA but also might have an effect on nearby ribosomes by affecting the appearance of 'traffic jams' where multiple ribosomes collide and form queues. To test this 'context effect' further, we here investigate the effect of the sequence of synonymous codons on the ribosome traffic by using a ribosome traffic model with codon-dependent rates, estimated from experiments. We compare the ribosome traffic on wild-type (WT) sequences and sequences where the synonymous codons were swapped randomly. By simulating translation of 87 genes, we demonstrate that the WT sequences, especially those with a high bias in codon usage, tend to have the ability to reduce ribosome collisions, hence optimizing the cellular investment in the translation apparatus. The magnitude of such reduction of the translation time might have a significant impact on the cellular growth rate and thereby have importance for the survival of the species.

A generic rate equation for catalysed, template-directed polymerisation

Biosynthetic networks link to growth and reproduction processes through template-directed synthesis of macromolecules such as polynucleotides and polypeptides. No rate equation exists that captures this link in a way that it can effectively be incorporated into a single computational model of the overall process. This paper describes the derivation of such a generic steady-state rate equation for catalysed, template-directed polymerisation reactions with varying monomer stoichiometry and varying chain length. The derivation is based on a classical Michaelis-Menten mechanism with template binding and an arbitrary number of chain elongation steps that produce a polymer composed of an arbitrary number of monomer types. The rate equation only requires the identity of the first dimer in the polymer sequence; for the remainder only the monomer composition needs be known. Further simplification of a term in the denominator yielded an equation requiring no positional information at all, only the monomer composition of the polymer; this equation still gave an excellent estimate of the reaction rate provided that either the monomer concentrations are at least half-saturating, or the polymer is very long.

Kinetic signatures of microRNA modes of action

MicroRNAs (miRNAs) are key regulators of all important biological processes, including development, differentiation, and cancer. Although remarkable progress has been made in deciphering the mechanisms used by miRNAs to regulate translation, many contradictory findings have been published that stimulate active debate in this field. Here we contribute to this discussion in three ways. First, based on a comprehensive analysis of the existing literature, we hypothesize a model in which all proposed mechanisms of microRNA action coexist, and where the apparent mechanism that is detected in a given experiment is determined by the relative values of the intrinsic characteristics of the target mRNAs and associated biological processes. Among several coexisting miRNA mechanisms, the one that will effectively be measurable is that which acts on or changes the sensitive parameters of the translation process. Second, we have created a mathematical model that combines nine known mechanisms of miRNA action and estimated the model parameters from the literature. Third, based on the mathematical modeling, we have developed a computational tool for discriminating among different possible individual mechanisms of miRNA action based on translation kinetics data that can be experimentally measured (kinetic signatures). To confirm the discriminatory power of these kinetic signatures and to test our hypothesis, we have performed several computational experiments with the model in which we simulated the coexistence of several miRNA action mechanisms in the context of variable parameter values of the translation.

Eukaryotic initiation factor 4F: a vulnerability of tumor cells

Protein synthesis is a complex, tightly regulated process in eukaryotic cells and its deregulation is a hallmark of many cancers. Translational control occurs primarily at the rate-limiting initiation step, where ribosomal subunits are recruited to template mRNAs through the concerted action of several eukaryotic initiation factors (eIFs). One factor that interacts with both the mRNA and ribosomes, and appears limiting for translation is eIF4F, a complex composed of the cap-binding protein, eIF4E; the scaffold protein, eIF4G; and the ATP-dependent DEAD-box helicase, eIF4A. eIF4E appears to play an important role in tumor initiation and progression since its overexpression can cooperate with oncogenes to accelerate transformation in cell lines and animal models, and its levels are elevated in many human cancers. This, therefore, represents a vulnerability for transformed cells, and presents an opportunity for therapeutic intervention. In this review, we discuss approaches for targeting eIF4F activity.

Mathematical and Computational Modelling of Ribosomal Movement and Protein Synthesis: an overview

Translation or protein synthesis consists of a complex system of chemical reactions, which ultimately result in decoding of the mRNA and the production of a protein. The complexity of this reaction system makes it difficult to quantitatively connect its input parameters (such as translation factor or ribosome concentrations, codon composition of the mRNA, or energy availability) to output parameters (such as protein synthesis rates or ribosome densities on mRNAs). Mathematical and computational models of translation have now been used for nearly five decades to investigate translation, and to shed light on the relationship between the different reactions in the system. This review gives an overview over the principal approaches used in the modelling efforts, and summarises some of the major findings that were made.

A systems view of the protein expression process

Many biological processes are regulated by changing the concentration and activity of proteins. The presence of a protein at a given subcellular location at a given time with a certain conformation is the result of an apparently sequential process. The rate of protein formation is influenced by chromatin state, and the rates of transcription, translation, and degradation. There is an exquisite control system where each stage of the process is controlled both by seemingly unregulated proteins as well as through feedbacks mediated by RNA and protein products. Here we review the biological facts and mathematical models for each stage of the protein production process. We conclude that advances in experimental techniques leading to a detailed description of the process have not been matched by mathematical models that represent the details of the process and facilitate analysis. Such an exercise is the first step towards development of a framework for a systems biology analysis of the protein production process.

ELECTRONIC SUPPLEMENTARY MATERIAL

The online version of this article (doi:10.1007/s11693-011-9088-1) contains supplementary material, which is available to authorized users.

Adenine-containing codons enhance protein synthesis by promoting mRNA binding to ribosomal 30S subunits provided that specific tRNAs are not exhausted

Adenines downstream of the initiation codon promote protein synthesis; however, some adenine-containing codons (AGA, AGG and AUA) at early positions inhibit protein synthesis when cognate tRNA is exhausted. It has also been reported, although not convincingly, the presence of adenines enhancing mRNA binding to the ribosome. To understand these apparent inconsistencies we analyzed the effect of these codons in mRNA-ribosome binding strength, mRNA stability, the production of peptidyl-tRNA (pep-tRNA) and protein synthesis. Constructs harboring lacZ derivatives were obtained by site directed mutagenesis where tandems of GGG, AGG, AGA, ATA and AAA codons were inserted at codon positions 2-3 and 3-4. Codons containing more adenines, irrespective of being common or rare, (AAA, ATA and AGA) promoted a higher synthesis of β-galactosidase (β-gal) in comparison with those rich in guanines (GGG and AGG) in a wild type transcription-translation system. Full-length mRNAs were also detected when the adenine-rich constructs were expressed in wild type cells. Under conditions where the pool of tRNAs is readily exhausted (pep-tRNA hydrolase defective cells), the adenine-rich lacZ derivatives caused a stronger and general inhibition of protein synthesis and cell growth. With the exception of the ATA lacZ derivative, only plasmid constructs containing hungry codons generated pep-tRNA (AGA and to a lesser extent AGG) in Pth defective cells. Codons containing more adenines clearly promoted lacZ mRNA binding to 30S subunit. The GGG lacZ mRNA showed a moderate increase in binding when mRNA secondary structures were disrupted by heating mRNAs before the binding assay which agrees with the lacZ mRNA secondary structures predicted with MFOLD. Altogether, these results indicate that mRNA binding to ribosome plays a major role in the enhancement of translation by adenine-rich codons irrespective of codon usage. This effect is naturally expressed in wild type systems and depends on adenine content, in contrast to the inhibition caused after over-expressing the lacZ derivatives containing rare codons in Pth defective cells.

Analysing GCN4 translational control in yeast by stochastic chemical kinetics modelling and simulation

Background

The yeast Saccharomyces cerevisiae responds to amino acid starvation by inducing the transcription factor Gcn4. This is mainly mediated via a translational control mechanism dependent upon the translation initiation eIF2·GTP·Met-tRNA_i^Met ternary complex, and the four short upstream open reading frames (uORFs) in its 5' mRNA leader. These uORFs act to attenuate GCN4 mRNA translation under normal conditions. During amino acid starvation, levels of ternary complex are reduced. This overcomes the GCN4 translation attenuation effect via a scanning/reinitiation control mechanism dependent upon uORF spacing.

Results

Using published experimental data, we have developed and validated a probabilistic formulation of GCN4 translation using the Chemical Master Equation (Model 1). Model 1 explains GCN4 translation's nonlinear dependency upon uORF placements, and predicts that an as yet unidentified factor, which was proposed to regulate GCN4 translation under some conditions, only has pronounced effects upon GCN4 translation when intercistronic distances are unnaturally short. A simpler Model 2 that does not include this unidentified factor could well represent the regulation of a natural GCN4 mRNA. Using parameter values optimised for this algebraic Model 2, we performed stochastic simulations by Gillespie algorithm to investigate the distribution of ribosomes in different sections of GCN4 mRNA under distinct conditions. Our simulations demonstrated that ribosomal loading in the 5'-untranslated region is mainly determined by the ratio between the rates of 5'-initiation and ribosome scanning, but was not significantly affected by rate of ternary complex binding. Importantly, the translation rate for codons starved of cognate tRNAs is predicted to be the most significant contributor to the changes in ribosomal loading in the coding region under repressing and derepressing conditions.

Conclusions

Our integrated probabilistic Models 1 and 2 explained GCN4 translation and helped to elucidate the role of a yet unidentified factor. The ensuing stochastic simulations evaluated different factors that may impact on the translation of GCN4 mRNA, and integrated translation status with ribosomal density.

Coupling between codon usage, translation and protein export in Escherichia coli

Proteins destined for export via the Sec-dependent pathway are synthesized with a short N-terminal signal peptide. A requirement for export is that the proteins are in a translocationally competent state. This is a loosely folded state that allows the protein to pass through the SecYEG apparatus and pass into the periplasm. In order to maintain pre-secretory proteins in an export-competent state, there are many factors that slow the folding of the pre-secretory protein in the cytoplasm. These include cytoplasmic chaperones, such as SecB, and the signal recognition particle, which bind the pre-secretory protein and direct it to the cytoplasmic membrane for export. Recently, evidence has been published that non-optimal codons in the signal sequence are important for a time-critical early event to allow the correct folding of pre-secretory proteins. This review details the recent developments in folding of the signal peptide and the pre-secretory protein.

A quantitative model for mRNA translation in Saccharomyces cerevisiae

Messenger RNA (mRNA) translation is an essential step in eukaryotic gene expression that contributes to the regulation of this process. We describe a deterministic model based on ordinary differential equations that describe mRNA translation in Saccharomyces cerevisiae. This model, which was parameterized using published data, was developed to examine the kinetic behaviour of translation initiation factors in response to amino acid availability. The model predicts that the abundance of the eIF1-eIF3-eIF5 complex increases under amino acid starvation conditions, suggesting a possible auxiliary role for these factors in modulating translation initiation in addition to the known mechanisms involving eIF2. Our analyses of the robustness of the mRNA translation model suggest that individual cells within a randomly generated population are sensitive to external perturbations (such as changes in amino acid availability) through Gcn2 signalling. However, the model predicts that individual cells exhibit robustness against internal perturbations (such as changes in the abundance of translation initiation factors and kinetic parameters). Gcn2 appears to enhance this robustness within the system. These findings suggest a trade-off between the robustness and performance of this biological network. The model also predicts that individual cells exhibit considerable heterogeneity with respect to their absolute translation rates, due to random internal perturbations. Therefore, averaging the kinetic behaviour of cell populations probably obscures the dynamic robustness of individual cells. This highlights the importance of single-cell measurements for evaluating network properties.

Analysis and prediction of translation rate based on sequence and functional features of the mRNA

Protein concentrations depend not only on the mRNA level, but also on the translation rate and the degradation rate. Prediction of mRNA's translation rate would provide valuable information for in-depth understanding of the translation mechanism and dynamic proteome. In this study, we developed a new computational model to predict the translation rate, featured by (1) integrating various sequence-derived and functional features, (2) applying the maximum relevance & minimum redundancy method and incremental feature selection to select features to optimize the prediction model, and (3) being able to predict the translation rate of RNA into high or low translation rate category. The prediction accuracies under rich and starvation condition were 68.8% and 70.0%, respectively, evaluated by jackknife cross-validation. It was found that the following features were correlated with translation rate: codon usage frequency, some gene ontology enrichment scores, number of RNA binding proteins known to bind its mRNA product, coding sequence length, protein abundance and 5′UTR free energy. These findings might provide useful information for understanding the mechanisms of translation and dynamic proteome. Our translation rate prediction model might become a high throughput tool for annotating the translation rate of mRNAs in large-scale.

Synonymous but not the same: the causes and consequences of codon bias

Despite their name, synonymous mutations have significant consequences for cellular processes in all taxa. As a result, an understanding of codon bias is central to fields as diverse as molecular evolution and biotechnology. Although recent advances in sequencing and synthetic biology have helped to resolve longstanding questions about codon bias, they have also uncovered striking patterns that suggest new hypotheses about protein synthesis. Ongoing work to quantify the dynamics of initiation and elongation is as important for understanding natural synonymous variation as it is for designing transgenes in applied contexts.

A mathematical modelling framework for elucidating the role of feedback control in translation termination

Translation is the final stage of gene expression where messenger RNA is used as a template for protein polymerization from appropriate amino acids. Release of the completed protein requires a release factor protein acting at the termination/stop codon to liberate it. In this paper we focus on a complex feedback control mechanism involved in the translation and synthesis of release factor proteins, which has been observed in different systems. These release factor proteins are involved in the termination stage of their own translation. Further, mutations in the release factor gene can result in a premature stop codon. In this case translation can result either in early termination and the production of a truncated protein or readthrough of the premature stop codon and production of the complete release factor protein. Thus during translation of the release factor mRNA containing a premature stop codon, the full length protein negatively regulates its production by its action on a premature stop codon, while positively regulating its production by its action on the regular stop codon. This paper develops a mathematical modelling framework to investigate this complex feedback control system involved in translation. A series of models is established to carefully investigate the role of individual mechanisms and how they work together. The steady state and dynamic behaviour of the resulting models are examined both analytically and numerically.

Ribosome collisions and translation efficiency: optimization by codon usage and mRNA destabilization

Individual mRNAs are translated by multiple ribosomes that initiate translation with an interval of a few seconds. The ribosome speed is codon dependent, and ribosome queuing has been suggested to explain specific data for translation of some mRNAs in vivo. By modeling the stochastic translation process as a traffic problem, we here analyze conditions and consequences of collisions and queuing. The model allowed us to determine the on-rate (0.8 to 1.1 initiations/s) and the time (1 s) the preceding ribosome occludes initiation for Escherichia coli lacZ mRNA in vivo. We find that ribosome collisions and queues are inevitable consequences of a stochastic translation mechanism that reduce the translation efficiency substantially on natural mRNAs. The cells minimize collisions by having its mRNAs being unstable and by a highly selected codon usage in the start of the mRNA. The cost of mRNA breakdown is offset by the concomitant increase in translation efficiency.

Minigene-like inhibition of protein synthesis mediated by hungry codons near the start codon

Rare AGA or AGG codons close to the initiation codon inhibit protein synthesis by a tRNA-sequestering mechanism as toxic minigenes do. To further understand this mechanism, a parallel analysis of protein synthesis and peptidyl-tRNA accumulation was performed using both a set of lacZ constructs where AGAAGA codons were moved codon by codon from +2, +3 up to +7, +8 positions and a series of 3–8 codon minigenes containing AGAAGA codons before the stop codon. β-Galactosidase synthesis from the AGAAGA lacZ constructs (in a Pth defective in vitro system without exogenous tRNA) diminished as the AGAAGA codons were closer to AUG codon. Likewise, β-galactosidase expression from the reporter +7 AGA lacZ gene (plus tRNA, 0.25 μg/μl) waned as the AGAAGAUAA minigene shortened. Pth counteracted both the length-dependent minigene effect on the expression of β-galactosidase from the +7 AGA lacZ reporter gene and the positional effect from the AGAAGA lacZ constructs. The +2, +3 AGAAGA lacZ construct and the shortest +2, +3 AGAAGAUAA minigene accumulated the highest percentage of peptidyl-tRNA^Arg4. These observations lead us to propose that hungry codons at early positions, albeit with less strength, inhibit protein synthesis by a minigene-like mechanism involving accumulation of peptidyl-tRNA.

Nonassembled human chorionic gonadotropin subunits and alphaalpha-homodimers use fast-track processing in the secretory pathway in contrast to alphabeta-heterodimers

In multimeric glycoproteins, like glycoprotein hormones, mutual subunit interactions are required for correct folding, assembly, and transport in the secretory pathway. However, character and time course of these interactions need further elucidation. The influence of the glycoprotein hormone alpha-subunit (GPHalpha) on the folding of the human chorionic gonadotropin (hCG) beta-subunit (hCGbeta) in hCG alphabeta-heterodimers was investigated in [(35)S]Met/Cys-labeled JEG-3 cells. Completeness of disulfide bridge formation during the time course of folding was estimated by labeling with [(3)H]N-ethylmaleinimide of free thiol groups not yet consumed. Subunit association took place between immature hCGbeta (high (3)H/(35)S ratio) and almost completely folded GPHalpha. Analysis revealed a highly dynamic maturation process comprising of at least eight main hCGbeta folding intermediates (molecular masses from 107 to 28 kDa) that could be micro-preparatively isolated and characterized. These hCGbeta variants developed while being associated with GPHalpha. The 107-kDa variant was identified as a complex with calnexin. In contrast to hCG alphabeta-heterodimers, free nonassociated hCGbeta, free large GPHalpha, and GPHalphaalpha homodimers showed a fast-track-like processing in the secretory pathway. At 10 min before hCG secretion, sialylation of these variants had already been completed in the late Golgi, whereas hCG alphabeta-heterodimers had still not arrived medial Golgi. This shows that the GPHalpha in the hCG alphabeta-heterodimers decelerates the maturation of the hCGbeta portion in the heterodimer complex. This results in a postponed approval of hCG alphabeta-heterodimers by the endoplasmic reticulum quality control unlike GPHalphaalpha homodimers, free hCGbeta, and GPHalpha subunits.

Traffic of interacting ribosomes: effects of single-machine mechanochemistry on protein synthesis

Many ribosomes simultaneously move on the same messenger RNA (mRNA), each synthesizing separately a copy of the same protein. In contrast to the earlier models, here we develop a "unified" theoretical model that not only incorporates the mutual exclusions of the interacting ribosomes, but also describes explicitly the mechanochemistry of each of these macromolecular machines during protein synthesis. Using analytical and numerical techniques of nonequilibrium statistical mechanics, we analyze the rates of protein synthesis and the spatiotemporal organization of the ribosomes in this model. We also predict how these properties would change with the changes in the rates of the various chemomechanical processes in each ribosome. Finally, we illustrate the power of this model by making experimentally testable predictions on the rates of protein synthesis and the density profiles of the ribosomes on some mRNAs in E-coli.

A model of protein translation including codon bias, nonsense errors, and ribosome recycling

We present and analyse a model of protein translation at the scale of an individual messenger RNA (mRNA) transcript. The model we develop is unique in that it incorporates the phenomena of ribosome recycling and nonsense errors. The model conceptualizes translation as a probabilistic wave of ribosome occupancy traveling down a heterogeneous medium, the mRNA transcript. Our results show that the heterogeneity of the codon translation rates along the mRNA results in short-scale spikes and dips in the wave. Nonsense errors attenuate this wave on a longer scale while ribosome recycling reinforces it. We find that the combination of nonsense errors and codon usage bias can have a large effect on the probability that a ribosome will completely translate a transcript. We also elucidate how these forces interact with ribosome recycling to determine the overall translation rate of an mRNA transcript. We derive a simple cost function for nonsense errors using our model and apply this function to the yeast (Saccharomyces cervisiae) genome. Using this function we are able to detect position dependent selection on codon bias which correlates with gene expression levels as predicted a priori. These results indirectly validate our underlying model assumptions and confirm that nonsense errors can play an important role in shaping codon usage bias.

Model-based inference of gene expression dynamics from sequence information

A dynamic model of prokaryotic gene expression is developed that makes considerable use of gene sequence information. The main contribution arises from the fact that the combined gene expression model allows us to access the impact of altering a nucleotide sequence on the dynamics of gene expression rates mechanistically. The high level of detail of the mathematical model is considered as an important step towards bringing together the tremendous amount of biological in-depth knowledge that has been accumulated at the molecular level, using a systems level analysis (in the sense of a bottom-up, inductive approach). This enables to the model to provide highly detailed insights into the various steps of the protein expression process and it allows us to access possible targets for model-based design. Taken as a whole, the mathematical gene expression model presented in this study provides a comprehensive framework for a thorough analysis of sequence-related effects on the stages of mRNA synthesis, mRNA degradation and ribosomal translation, as well as their nonlinear interconnectedness. Therefore, it may be useful in the rational design of recombinant bacterial protein synthesis systems, the modulation of enzyme activities in pathway design, in vitro protein biosynthesis, and RNA-based vaccination.

Studies on the posttranscriptional site of cAMP action in the regulation of the synthesis of tyrosine aminotransferase

Synthesis of L-tyrosine:2-oxoglutarate aminotransferase (EC 2.6.1.5) can be induced by N6,O2'-dibutyryl-adenosine 3',5'-monophosphate (Bt2cAMP) in Reuber H35 cell cultures. Actinomycin D fails to block this induction which indicates a target for Bt2cAMP at a posttranscriptional level. We have determined the influence of Bt2cAMP on several translational events during the tyrosine aminotransferase synthesis with the following results. (1) The number of nascent tyrosine aminotransferase chains increased, whereas no effect was measured on the number of nascent total protein chains. (2) The rate of elongation along the tyrosine aminotransferase mRNA and total mRNA is not enhanced by Bt2cAMP. (3) The induced synthesis of tyrosine aminotransferase is more sensitive to the inhibition of elongation. We conclude from our results that Bt2cAMP induces the synthesis of tyrosine aminotransferase by an increase in the rate of initiation on the tyrosine aminotransferase mRNA.

Hsp70 and Hsp90--a relay team for protein folding

Molecular chaperones are a functionally defined set of proteins which assist the structure formation of proteins in vivo. Without certain protective mechanisms, such as binding nascent polypeptide chains by molecular chaperones, cellular protein concentrations would lead to misfolding and aggregation. In the mammalian system, the molecular chaperones Hsp70 and Hsp90 are involved in the folding and maturation of key regulatory proteins, like steroid hormone receptors, transcription factors, and kinases, some of which are involved in cancer progression. Hsp70 and Hsp90 form a multichaperone complex, in which both are connected by a third protein called Hop. The connection of and the interplay between the two chaperone machineries is of crucial importance for cell viability. This review provides a detailed view of the Hsp70 and Hsp90 machineries, their cofactors and their mode of regulation. It summarizes the current knowledge in the field, including the ATP-dependent regulation of the Hsp70/Hsp90 multichaperone cycle and elucidates the complex interplay and their synergistic interaction.

Totally asymmetric exclusion process with extended objects: a model for protein synthesis

The process of protein synthesis in biological systems resembles a one dimensional driven lattice gas in which the particles have spatial extent, covering more than one lattice site. We expand the well studied totally asymmetric exclusion process, in which particles typically cover a single lattice site, to include cases with extended objects. Exact solutions can be determined for a uniform closed system. We analyze the uniform open system through two approaches. First, a continuum limit produces a modified diffusion equation for particle density profiles. Second, an extremal principle based on domain wall theory accurately predicts the phase diagram and currents in each phase. Finally, we briefly consider approximate approaches to a nonuniform open system with quenched disorder in the particle hopping rates and compare these approaches with Monte Carlo simulations.

The pair of arginine codons AGA AGG close to the initiation codon of the lambda int gene inhibits cell growth and protein synthesis by accumulating peptidyl-tRNAArg4

To analyse the mechanism by which rare codons near the initiation codon inhibit cell growth and protein synthesis, we used the bacteriophage lambda int gene or early codon substitution derivatives. The lambda int gene has a high frequency of rare ATA, AGA and AGG codons; two of them (AGA AGG) located at positions 3 and 4 of the int open reading frame (ORF). Escherichia coli pth (rap) cells, which are defective in peptidyl-tRNA hydrolase (Pth) activity, are more susceptible to the inhibitory effects of int expression as compared with wild-type cells. Cell growth and Int protein synthesis were enhanced by overexpression of Pth and tRNAArg4 cognate to AGG and AGA but not of tRNAIle2a specific for ATA. The increase of Int protein synthesis also takes place when the rare arginine codons AGA and AGG at positions 3 and 4 are changed to common arginine CGT or lysine AAA codons but not to rare isoleucine ATA codons. In addition, overexpression of int in Pth defective cells provokes accumulation of peptidyl-tRNAArg4 in the soluble fraction. Therefore, cell growth and Int synthesis inhibition may be due to ribosome stalling and premature release of peptidyl-tRNAArg4 from the ribosome at the rare arginine codons of the first tandem, which leads to cell starvation for the specific tRNA.

Translation at higher than an optimal level interferes with coupling at an intercistronic junction

In pairs of adjacent genes co-transcribed on bacterial polycistronic mRNAs, translation of the first coding region frequently functions as a positive factor to couple translation to the distal coding region. Coupling efficiencies vary over a wide range, but synthesis of both gene products at similar levels is common. We report the results of characterizing an unusual gene pair, in which only about 1% of the translational activity from the upstream gene is transmitted to the distal gene. The inefficient coupling was unexpected because the upstream gene is highly translated, the distal initiation site has weak but intrinsic ability to bind ribosomes, and the AUG is only two nucleotides beyond the stop codon for the upstream gene. The genes are those in the filamentous phage IKe genome, which encode the abundant single-stranded DNA binding protein (gene V) and the minor coat protein that caps one tip of the phage (gene VII). Here, we have used chimeras between the related phage IKe and f1 sequences to localize the region responsible for inefficient coupling. It mapped upstream from the intercistronic region containing the gene V stop codon and the gene VII initiation site, indicating that low coupling efficiency is associated with gene V. The basis for inefficient coupling emerged when coupling efficiency was found to increase as gene V translation was decreased below the high wild-type level. This was achieved by lowering the rate of elongation and by decreasing the efficiency of suppression at an amber codon within the gene. Increasing the strength of the Shine-Dalgarno interaction with 16S rRNA at the gene VII start also increased coupling efficiency substantially. In this gene pair, upstream translation thus functions in an unprecedented way as a negative factor to limit downstream expression. We interpret the results as evidence that translation in excess of an optimal level in an upstream gene interferes with coupling in the intercistronic junction.

Improved heterologous expression of human glutathione transferase A4-4 by random silent mutagenesis of codons in the 5' region

Glutathione transferase A4-4 (GST A4-4) is involved in the detoxication of lipid peroxidation products such as alkenals. The human enzyme has been heterologously expressed in Escherichia coli, but for more extensive characterization of the enzyme the expression level had to be elevated. A clone providing up to 8-fold higher yields was created, by screening an expression library with random silent mutations in the 5' region of the cDNA encoding GST A4-4.

Metabolic control analysis of monoclonal antibody synthesis

A general route for protein synthesis in eukaryotic cells has been proposed and applied to monoclonal antibody (MAb) synthesis. It takes into account transcription of the gene, binding of ribosomes to mRNA, and polypeptide elongation including binding to SRP (signal recognition particles) and SRP-receptor, competing translocation, folding and glycosylation, assembly of the heavy and light chains in a tetrameric protein and Golgi processing and secretion. A comprehensive model was built on the basis of the proposed pathway. The model takes into account the mechanism of each step. Metabolic control analysis (MCA) principles were applied to the general pathway using the proposed model, and control coefficients were calculated. The results show a shared flux control (of both pathway flux and flux ratio at the branch) among different steps, i.e., transcription, folding, glycosylation, translocation and building blocks synthesis. The steps sharing the control depend on the concentration of building blocks, pathway flux and levels of OST (oligosacharyl transferase), BiP (heavy chain binding protein) and PDI (protein disulfide isomerase). Model predictions compare well with experimental data for MAb synthesis, explaining the control structure of the route and the heterogeneity of the product and also addressing future targets for improvement of the production rate of MAbs.

Phosphorylation of p70(S6k) correlates with increased skeletal muscle mass following resistance exercise

High-resistance exercise training results in an increase in muscle wet mass and protein content. To begin to address the acute changes following a single bout of high-resistance exercise, a new model has been developed. Training rats twice a week for 6 wk resulted in 13.9 and 14.4% hypertrophy in the extensor digitorum longus (EDL) and tibialis anterior (TA) muscles, respectively. Polysome profiles after high-resistance lengthening contractions suggest that the rate of initiation is increased. The activity of the 70-kDa S6 protein kinase (p70(S6k)), a regulator of translation initiation, is also increased following high-resistance lengthening contractions (TA, 363 +/- 29%; EDL, 353 +/- 39%). Furthermore, the increase in p70(S6k) activity 6 h after exercise correlates with the percent change in muscle mass after 6 wk of training (r = 0.998). The tight correlation between the activation of p70(S6k) and the long-term increase in muscle mass suggests that p70(S6k) phosphorylation may be a good marker for the phenotypic changes that characterize muscle hypertrophy and may play a role in load-induced skeletal muscle growth.

Posttranscriptional control of gene expression in yeast

Studies of the budding yeast Saccharomyces cerevisiae have greatly advanced our understanding of the posttranscriptional steps of eukaryotic gene expression. Given the wide range of experimental tools applicable to S. cerevisiae and the recent determination of its complete genomic sequence, many of the key challenges of the posttranscriptional control field can be tackled particularly effectively by using this organism. This article reviews the current knowledge of the cellular components and mechanisms related to translation and mRNA decay, with the emphasis on the molecular basis for rate control and gene regulation. Recent progress in characterizing translation factors and their protein-protein and RNA-protein interactions has been rapid. Against the background of a growing body of structural information, the review discusses the thermodynamic and kinetic principles that govern the translation process. As in prokaryotic systems, translational initiation is a key point of control. Modulation of the activities of translational initiation factors imposes global regulation in the cell, while structural features of particular 5' untranslated regions, such as upstream open reading frames and effector binding sites, allow for gene-specific regulation. Recent data have revealed many new details of the molecular mechanisms involved while providing insight into the functional overlaps and molecular networking that are apparently a key feature of evolving cellular systems. An overall picture of the mechanisms governing mRNA decay has only very recently begun to develop. The latest work has revealed new information about the mRNA decay pathways, the components of the mRNA degradation machinery, and the way in which these might relate to the translation apparatus. Overall, major challenges still to be addressed include the task of relating principles of posttranscriptional control to cellular compartmentalization and polysome structure and the role of molecular channelling in these highly complex expression systems.

Growing up in a dangerous environment: a network of multiple targeting and folding pathways for nascent polypeptides in the cytosol

The first events in the lives of proteins are the most hazardous. Starting at the ribosome, nascent polypeptides undergo complex folding processes endangered by aggregation reactions. Proteins with organellar destinations require correct targeting to the translocation machineries and prevention from premature folding. The high precision and speed of these processes is ensured by a cystosolic system consisting of molecular chaperones, folding catalysts and targeting factors. This review focuses on the interactions of this system with nascent polypeptides and discusses new concepts for protein folding in the cytosol. It is proposed that folding and targeting are promoted by a flexible network of multiple unassisted and assisted pathways.

Codon pair utilization biases influence translational elongation step times

Two independent assays capable of measuring the relative in vivo translational step times across a selected codon pair in a growing polypeptide in the bacterium Escherichia coli have been employed to demonstrate that codon pairs observed in protein coding sequences more frequently than predicted (over-represented codon pairs) are translated slower than pairs observed less frequently than expected (under-represented codon pairs). These results are consistent with the findings that translational step times are influenced by codon context and that these context effects are related to the compatabilities of adjacent tRNA isoacceptor molecules on the surface of a translating ribosome. These results also support our previous suggestion that the frequency of one codon next to another has co-evolved with the structure and abundance of tRNA isoacceptors in order to control the rates of translational step times without imposing additional constraints on amino acid sequences or protein structures.

Dietary fat elevates hepatic apoA-I production by increasing the fraction of apolipoprotein A-I mRNA in the translating pool

Elevated plasma high density lipoprotein cholesterol (HDL-C) levels are associated with a decreased risk for coronary heart disease. Ironically, diets enriched in saturated fat and cholesterol (HF/HC diets), which tend to accelerate atherosclerotic processes by increasing LDL cholesterol levels, also raise HDL-C. We have recently reported, using a human apoA-I (hapoA-1) transgenic mouse model, that the elevation of HDL-C by a HF/HC diet is attributable, in part, to an increase in the hepatic production of hapoA-1. To further define the hepatocellular processes associated with this induction, we have prepared primary hepatocytes from hapoA-1 transgenic mice. Rates of hapoA-1 secretion were 40% greater from cells prepared from animals fed the HF/HC relative to a low fat-low cholesterol (LF/LC) control diet. The abundance of hapoA-1 mRNA in these cells was similar between hepatocytes prepared from the HF/HC and LF/LC diet fed animals, suggesting a post-transcriptional mechanism that does not involve mRNA stability. Inhibition of secretion using brefeldin A revealed an increase in cellular hapoA-1 accumulation. Thus, the HF/HC diet apparently affects hepatic hapoA-1 production via a mechanism that is manifest prior to the exit of newly synthesized hapoA-1 from the Golgi. Pulse-chase experiments revealed a 39% greater peak hapoA-1 synthesis, with no difference in the degradation of total labeled hapoA-1 protein, as a result of the HF/HC diet feeding. Finally, resolution of liver S10 extracts via sucrose density sedimentation and metrizamide density equilibrium gradient centrifugation analyses both revealed similar increases (31 and 24%, respectively) in the relative percentage of hapoA-1 mRNA associated with the translating polysomal fractions as a result of the HF/HC feeding. Together, these data suggest that the HF/HC diet affects hepatic hapoA-1 production via a specific modulation in the relative amount of hapoA-1 mRNA in the polysomal pool. These observations provide an opportunity to explore a new mechanism regulating apoA-1 production and might lead to the development of novel therapies to elevate plasma HDL-C levels.

Decrease in heart peptide initiation during head-down tilt may be modulated by HSP-70

This study examines the mechanism of the rapid decrease in cardiac muscle protein synthesis during rodent hindlimb non-weight bearing. Polysomes isolated from rat hearts 8 h after suspension show less RNA in the polysome pool and a shift in polysome size toward fewer ribosomes per mRNA; 18 h after suspension, the size shift persists, but the amount of RNA in the polysome pool returns to control values. These data are consistent with a decrease in the rate of initiation of protein synthesis. At both 8 and 12 h of suspension, the cardiac polysomes show a 78 and 93% increase association with the nascent polypeptide chaperone protein 70-kDa heat-shock cognate/heat-shock protein (HSC/HSP-70), respectively, that persists after 7 days of non-weight bearing. Because the dissociation of HSC/HSP-70 from unfolded protein can be modulated by ATP, we measured the adenosine nucleotide pools and found a 53% decrease in ATP levels after 18 h of suspension. We propose a mechanism in which a shift of HSC/HSP-70 to the nascent polypeptide indirectly inhibits protein synthesis initiation.

Efficient translation of poly(A)-deficient mRNAs in Saccharomyces cerevisiae

The polyadenylate tail of eukaryotic mRNAs is thought to influence various metabolic phenomena including mRNA stability, translation initiation, and nucleo-cytoplasmic transport. We have analyzed the fate of mRNAs following inactivation of poly(A) polymerase in Saccharomyces cerevisiae containing a temperature-sensitive, lethal mutation (pap1-1) in the gene for poly(A) polymerase (PAP1). Inactivation of poly(A) polymerase (Pap1) by shifting cells to the nonpermissive temperature resulted in the loss of at least 80% of measurable poly(A) within 60 min. Northern blot analysis revealed the disappearance of some mRNAs (CYH2 and HIS4) consistent with a role for poly(A) tails in mRNA stability. However, other mRNAs (TCM1, PAB1, ACT1, and HTB2) accumulate as poly(A)-deficient (A < approximately 25) transcripts as defined by an inability to bind oligo(dT)-cellulose. Sucrose density gradient analysis of polyribosomes revealed a twofold reduction in the amount of each size class of polyribosomes in shifted cells and a commensurate increase in free ribosomes. However, poly(A)-deficient mRNAs in shifted cells remain associated with the same size polyribosomes as poly(A)+ mRNAs in unshifted cells, indicating normal initiation of translation. RNase mapping of transcripts from pap1-1 cells revealed PAB1 mRNA to be poly(A)- whereas TCM1 exists as equal amounts of poly(A)- and poly(A)+ mRNA 60 min after shift. Interestingly, both of these classes of TCM1 mRNA appear in similar amounts in each polyribosome fraction indicating that ribosomes may not distinguish between them. These findings suggest that under conditions of excess translational capacity, poly(A)- and poly(A)+ mRNAs may initiate translation with comparable efficiencies.

Soleus muscle nascent polypeptide chain elongation slows protein synthesis rate during non-weight-bearing activity

Protein synthesis rate of the soleus muscle decreases rapidly during non-weight-bearing activity. We isolated polysomes from 18-h non-weight-bearing soleus muscle to investigate the mechanism of this phenomenon. The distribution of polysomal alpha-actin mRNA and 18S rRNA on sucrose density gradients shows that polysomes shift to larger sizes (more ribosomes per mRNA) during non-weight-bearing activity. Furthermore, RNA is mobilized into the polysome pool of the non-weight-bearing soleus muscle; these data indicate that initiation of protein synthesis is not rate limiting. We explain these results as the slowing of nascent polypeptide chain elongation, such that there is a "traffic jam" of ribosomes on the mRNAs, increasing the number of ribosomes per mRNA while, at the same time, decreasing protein synthesis rate. In support of this hypothesis, myoblasts treated with a low dose of cycloheximide (a specific elongation inhibitor) show a similar shift in polysome size. A numerical model of protein synthesis further shows that elongation is more effective than initiation and termination in affecting protein synthesis and polysome size. We conclude that the non-weight-bearing-induced decrease in postural muscle protein synthesis rate is initially caused by slowing of nascent polypeptide chain elongation.

Importance of mRNA folding and start codon accessibility in the expression of genes in a ribosomal protein operon of Escherichia coli

The trmD operon of Escherichia coli consists of the genes for the ribosomal protein (r-protein) S16, a 21 kilodalton protein (21K) of unknown function, the tRNA(m1G37)methyltransferase (TrmD), and r-protein L19, in that order. The synthesis of the 21K and TrmD proteins is 12 and 40-fold lower, respectively, than that of the two r-proteins, although the corresponding parts of the mRNA are equally abundant. This translational control of expression of at least the 21K protein gene (21K), is mediated by a negative control element located between codons 18 and 50 of 21K. Here, we present evidence for a model in which mRNA sequences up to around 100 nucleotides downstream from the start codon of 21K fold back and base-pair to the 21K translation initiation region, thereby decreasing the translation initiation frequency. Mutations in the internal negative control element of 21K that would prevent the formation of the proposed mRNA secondary structure over both the Shine-Dalgarno (SD) sequence and the start codon increased expression up to about 20-fold, whereas mutations that would disrupt the base-pairing with the SD-sequence had only relatively small effects on expression. In addition, the expression increased 12-fold when the stop codon of the preceding gene, rpsP, was moved next to the SD-sequence of 21K allowing the ribosomes to unfold the postulated mRNA secondary structure. The expression increased up to 150-fold when that stop codon change was combined with the internal negative control element base-substitutions that derepressed translation about 20-fold. The negative control element of 21K does not seem to be responsible for the low expression of the trmD gene located downstream. However, a similar negative control element native to trmD can explain at least partly the low expression of trmD. Possibly, the two mRNA secondary structures function to decouple translation of 21K and trmD from that of the respective upstream cistron in order to achieve their independent regulation.

Growth rate dependence of global amino acid composition

The global amino acid composition of bacteria growing in different media has been studied. The data reveal significant changes in the amino acid composition in the growth rate range between 0.5 and 2.1 doublings per hour at 37 degrees C. The changes are consistent with a progressive simplification of the protein population and mRNA pools as the growth rates increase.