Model of ribosome translation and mRNA unwinding

Bioinformatic Assessment of Factors Affecting the Correlation between Protein Abundance and Elongation Efficiency in Prokaryotes

Protein abundance is crucial for the majority of genetically regulated cell functions to act properly in prokaryotic organisms. Therefore, developing bioinformatic methods for assessing the efficiency of different stages of gene expression is of great importance for predicting the actual protein abundance. One of these steps is the evaluation of translation elongation efficiency based on mRNA sequence features, such as codon usage bias and mRNA secondary structure properties. In this study, we have evaluated correlation coefficients between experimentally measured protein abundance and predicted elongation efficiency characteristics for 26 prokaryotes, including non-model organisms, belonging to diverse taxonomic groups The algorithm for assessing elongation efficiency takes into account not only codon bias, but also number and energy of secondary structures in mRNA if those demonstrate an impact on predicted elongation efficiency of the ribosomal protein genes. The results show that, for a number of organisms, secondary structures are a better predictor of protein abundance than codon usage bias. The bioinformatic analysis has revealed several factors associated with the value of the correlation coefficient. The first factor is the elongation efficiency optimization type-the organisms whose genomes are optimized for codon usage only have significantly higher correlation coefficients. The second factor is taxonomical identity-bacteria that belong to the class Bacilli tend to have higher correlation coefficients among the analyzed set. The third is growth rate, which is shown to be higher for the organisms with higher correlation coefficients between protein abundance and predicted translation elongation efficiency. The obtained results can be useful for further improvement of methods for protein abundance prediction.

Stochastic thermodynamics and modes of operation of a ribosome: A network theoretic perspective

The ribosome is one of the largest and most complex macromolecular machines in living cells. It polymerizes a protein in a step-by-step manner as directed by the corresponding nucleotide sequence on the template messenger RNA (mRNA) and this process is referred to as "translation" of the genetic message encoded in the sequence of mRNA transcript. In each successful chemomechanical cycle during the (protein) elongation stage, the ribosome elongates the protein by a single subunit, called amino acid, and steps forward on the template mRNA by three nucleotides called a codon. Therefore, a ribosome is also regarded as a molecular motor for which the mRNA serves as the track, its step size is that of a codon and two molecules of GTP and one molecule of ATP hydrolyzed in that cycle serve as its fuel. What adds further complexity is the existence of competing pathways leading to distinct cycles, branched pathways in each cycle, and futile consumption of fuel that leads neither to elongation of the nascent protein nor forward stepping of the ribosome on its track. We investigate a model formulated in terms of the network of discrete chemomechanical states of a ribosome during the elongation stage of translation. The model is analyzed using a combination of stochastic thermodynamic and kinetic analysis based on a graph-theoretic approach. We derive the exact solution of the corresponding master equations. We represent the steady state in terms of the cycles of the underlying network and discuss the energy transduction processes. We identify the various possible modes of operation of a ribosome in terms of its average velocity and mean rate of GTP hydrolysis. We also compute entropy production as functions of the rates of the interstate transitions and the thermodynamic cost for accuracy of the translation process.

Non-tight and tight chemomechanical couplings of biomolecular motors under hindering loads

Biomolecular motors make use of free energy released from chemical reaction (typically ATP hydrolysis) to perform mechanical motion or work. An important issue is whether a molecular motor exhibits tight or non-tight chemomechanical (CM) coupling. The tight CM coupling refers to that each ATPase activity is coupled with a mechanical step, while the non-tight CM coupling refers to that an ATPase activity is not necessarily coupled with a mechanical step. Here, we take kinesin, monomeric DNA helicase, ring-shaped hexameric DNA helicase and ribosome as examples to study this issue. Our studies indicate that some motors such as kinesin, monomeric helicase and ribosome exhibit non-tight CM coupling under hindering forces, while others such as the ring-shaped hexameric helicase exhibit tight or nearly tight CM coupling under any force. For the former, the reduction of the velocity caused by the hindering force arises mainly from the reduction of the CM coupling efficiency, while the ATPase rate is independent or nearly independent of the force. For the latter, the reduction of the velocity caused by the hindering force arises mainly from the reduction of the ATPase rate, while the CM coupling efficiency is independent or nearly independent of the force.

Mechanism of ribosome translation through mRNA secondary structures

A ribosome is a macromolecular machine that is responsible for translating the genetic codes in messenger RNA (mRNA) into polypeptide chains. It has been determined that besides translating through the single-stranded region, the ribosome can also translate through the duplex region of mRNA by unwinding the duplex. To understand the mechanism of ribosome translation through the duplex, several models have been proposed to study the dynamics of mRNA unwinding. Here, we present a comprehensive review of these models and also discuss other possible models. We evaluate each model and discuss the consistency and/or inconsistency between the theoretical results that are obtained based on each model and the available experimental data, thus determining which model is the most reasonable one to describe the mRNA unwinding mechanism and dynamics of the ribosome. Moreover, a framework for future studies in this subject is provided.

Processivity of nucleic acid unwinding and translocation by helicases

Helicases are a class of enzymes that use the chemical energy of NTP hydrolysis to drive mechanical processes such as translocation and nucleic acid (NA) strand separation. Besides the NA unwinding speed, another important factor for the helicase activity is the NA unwinding processivity. Here, we study the NA unwinding processivity with an analytical model that captures the phenomenology of the NA unwinding process. First, we study the processivity of the non-hexameric helicase that can unwind NA efficiently in the form of a monomer and the processivity of the hexameric helicase that can unwind DNA effectively, providing quantitative explanations of the available single-molecule experimental data. Then, we study the processivity of the non-hexameric helicases, in particular UvrD, in the form of a dimer and compare with that in the form of a monomer. The available single-molecule and some biochemical data showing that while UvrD monomer is a highly processive single-stranded DNA translocase it is inactive in DNA unwinding, whereas other biochemical data showing that UvrD is active in both single-stranded DNA translocation and DNA unwinding in the form of a monomer can be explained quantitatively and consistently. In addition, the recent single-molecule data are also explained quantitatively showing that constraining the 2B subdomain in closed conformation by intramolecular cross-linking can convert Rep monomer with a very poor DNA unwinding activity into a superhelicase that can unwind more than thousands of DNA base pairs processively, even against a large opposing force. Proteins 2016; 84:1590-1605. © 2016 Wiley Periodicals, Inc.

Modeling of ribosome dynamics on a ds-mRNA under an external load

Protein molecules in cells are synthesized by macromolecular machines called ribosomes. According to the recent experimental data, we reduce the complexity of the ribosome and propose a model to express its activity in six main states. Using our model, we study the translation rate in different biological relevant situations in the presence of external force and the translation through the RNA double stranded region in the absence or presence of the external force. In the present study, we give a quantitative theory for translation rate and show that the ribosome behaves more like a Brownian Ratchet motor. Our findings could shed some light on understanding behaviors of the ribosome in biological conditions.

Model of the pathway of -1 frameshifting: Long pausing

It has been characterized that the programmed ribosomal -1 frameshifting often occurs at the slippery sequence on the presence of a downstream mRNA pseudoknot. In some prokaryotic cases such as the dnaX gene of Escherichia coli, an additional stimulatory signal-an upstream, internal Shine-Dalgarno (SD) sequence-is also necessary to stimulate the efficient -1 frameshifting. However, the molecular and physical mechanism of the -1 frameshifting is poorly understood. Here, we propose a model of the pathway of the -1 translational frameshifting during ribosome translation of the dnaX -1 frameshift mRNA. With the model, the single-molecule fluorescence data (Chen et al. (2014) [29]) on the dynamics of the shunt either to long pausing or to normal translation, the tRNA transit and sampling dynamics in the long-paused rotated state, the EF-G sampling dynamics, the mean rotated-state lifetimes, etc., are explained quantitatively. Moreover, the model is also consistent with the experimental data (Yan et al. (2015) [30]) on translocation excursions and broad branching of frameshifting pathways. In addition, we present some predicted results, which can be easily tested by future optical trapping experiments.

Dwell-Time Distribution, Long Pausing and Arrest of Single-Ribosome Translation through the mRNA Duplex

Proteins in the cell are synthesized by a ribosome translating the genetic information encoded on the single-stranded messenger RNA (mRNA). It has been shown that the ribosome can also translate through the duplex region of the mRNA by unwinding the duplex. Here, based on our proposed model of the ribosome translation through the mRNA duplex we study theoretically the distribution of dwell times of the ribosome translation through the mRNA duplex under the effect of a pulling force externally applied to the ends of the mRNA to unzip the duplex. We provide quantitative explanations of the available single molecule experimental data on the distribution of dwell times with both short and long durations, on rescuing of the long paused ribosomes by raising the pulling force to unzip the duplex, on translational arrests induced by the mRNA duplex and Shine-Dalgarno(SD)-like sequence in the mRNA. The functional consequences of the pauses or arrests caused by the mRNA duplex and the SD sequence are discussed and compared with those obtained from other types of pausing, such as those induced by "hungry" codons or interactions of specific sequences in the nascent chain with the ribosomal exit tunnel.

A unified model of nucleic acid unwinding by the ribosome and the hexameric and monomeric DNA helicases

DNA helicases are enzymes that use the chemical energy to separate DNA duplex into their single-stranded forms. The ribosome, which catalyzes the translation of messenger RNAs (mRNAs) into proteins, can also unwind mRNA duplex. According to their structures, the DNA helicases can fall broadly into hexameric and monomeric forms. A puzzling issue for the monomeric helicases is that although they have similar structures, in vitro biochemical data showed convincingly that in the monomeric forms some have very weak DNA unwinding activities, some have relatively high unwinding activities while others have high unwinding activities. However, in the dimeric or oligomeric forms all of them have high unwinding activities. In addition, in the monomeric forms all of them can translocate efficiently along the single-stranded DNA (ssDNA). Here, we propose a model of the translocation along the ssDNA and DNA unwinding by the monomeric helicases, providing a consistent explanation of these in vitro experimental data. Moreover, by comparing the present model for the monomeric helicases with the model for the hexameric helicases and that for the ribosome which were proposed before, a unified model of nucleic acid unwinding by the three enzymes is proposed.

Biphasic character of ribosomal translocation and non-Michaelis-Menten kinetics of translation

We study theoretically the kinetics of mRNA translocation in the wild-type (WT) Escherichia coli ribosome, which is composed of a small 30S and large 50S subunit, and the ribosomes with mutations to some intersubunit bridges such as B1a, B4, B7a, and B8. The theoretical results reproduce well the available in vitro experimental data on the biphasic kinetics of the forward mRNA translocation catalyzed by elongation factor G (EF-G) hydrolyzing GTP, which can be best fit by the sum of two exponentials, and the monophasic kinetics of the spontaneous reverse mRNA translocation in the absence of the elongation factor, which can be best fit by a single-exponential function, in both the WT and mutant ribosomes. We show that both the mutation-induced increase in the maximal rate of the slow phase for the forward mRNA translocation and that in the rate of the spontaneous reverse mRNA translocation result from a reduction in the intrinsic energy barrier to resist the rotational movements between the two subunits, giving the same degree of increase in the two rates. The mutation-induced increase in the maximal rate of the fast phase for the forward mRNA translocation results mainly from the increase in the rate of the ribosomal unlocking, a conformational change in the ribosome that widens the mRNA channel for the mRNA translocation to take place, which could be partly due to the effect of the mutation on the intrasubunit 30S head rotation. Moreover, we study the translation rate of the WT and mutant ribosomes. It is shown that the translation rate versus the concentration of EF-G-GTP does not follow the Michaelis-Menten (MM) kinetics, which is in sharp contrast to the general property of other enzymes that the rate of the enzymatic reaction versus the concentration of a substrate follows the MM kinetics. The physical origin of this non-MM kinetics for the ribosome is revealed.

Model of EF4-induced ribosomal state transitions and mRNA translocation

EF4, a highly conserved protein present in bacteria, mitochondria and chloroplasts, can bind to both the posttranslocation and pretranslocation ribosomal complexes. When binding to the posttranslocation state, it catalyzes backward translocation to a pretranslocation state. When binding to the pretranslocation state, it catalyzes transition to another pretranslocation state that is similar and possibly identical to that resulting from the posttranslocation state bound by EF4, and competes with EF-G to regulate the elongation cycle. However, the molecular mechanism on how EF4 induces state transitions and mRNA translocation remains unclear. Here, we present both the model for state transitions induced by EF4 binding to the posttranslocation state and that by EF4 binding to the pretranslocation state, based on which we study the kinetics of EF4-induced state transitions and mRNA translocation, giving quantitative explanations of the available experimental data. Moreover, we present some predicted results on state transitions and mRNA translocation induced by EF4 binding to the pretranslocation state complexed with the mRNA containing a duplex region.

Dynamics of +1 ribosomal frameshifting

It has been well characterized that the amino acid starvation can induce +1 frameshifting. However, how the +1 frameshifting occurs has not been fully understood. Here, taking Escherichia coli RF2 programmed frameshifting as an example we present systematical analysis of the +1 frameshifting that could occur during every state-transition step in elongation phase of protein synthesis, showing that the +1 frameshifting can occur only during the period after deacylated tRNA dissociation from the posttranslocation state and before the recognition of the next "hungry" codon. The +1 frameshifting efficiency is theoretically studied, with the simple analytical solutions showing that the high efficiency is almost solely due to the occurrence of ribosome pausing which in turn results from the insufficient RF2. The analytical solutions also provide a consistent explanation of a lot of independent experimental data.

Dynamics of forward and backward translocation of mRNA in the ribosome

Translocation of the mRNA-tRNA complex in the ribosome, which is catalyzed by elongation factor EF-G, is one of critical steps in the elongation cycle of protein synthesis. Besides this conventional forward translocation, the backward translocation can also occur, which can be catalyzed by elongation factor LepA. However, the molecular mechanism of the translocation remains elusive. To understand the mechanism, here we study theoretically the dynamics of the forward translocation under various nucleotide states of EF-G and the backward translocation in the absence of and in the presence of LepA. We present a consistent explanation of spontaneous forward translocations in the absence of EF-G, the EF-G-catalyzed forward translocations in the presence of a non-hydrolysable GTP analogue and in the presence of GTP, and the spontaneous and LepA-catalyzed backward translocation. The theoretical results provide quantitative explanations of a lot of different, independent experimental data, and also provide testable predictions.

Translocation dynamics of tRNA-mRNA in the ribosome

Translocation of tRNA-mRNA complex in the ribosome is an essential step in the elongation cycle of protein synthesis. However, some important issues concerning the molecular mechanism of the tRNA-mRNA translocation catalyzed by EF-G.GTP or by EF-G.GDPNP remain controversial. For example, can EF-G.GTP selectively bind to the hybrid pretranslocation state or bind to both the non-rotated pretranslocation and the hybrid pretranslocation states? Does the greater potency of EF-G in the presence of GTP rather than GDPNP in facilitating translocation derive from the effects on transition from the classical non-rotated to hybrid state (the first step of the translocation) or on transition from the hybrid to posttranslocation state (the second step)? Here, based on our proposed model, we study theoretically the dynamics of the tRNA-mRNA translocation through the ribosome catalyzed by EF-G.GTP and by EF-G.GDPNP. By comparing our theoretical results with the available experimental data, we show that EF-G.GTP can also bind to the classical non-rotated pretranslocation state and the greater potency of GTP hydrolysis in facilitating translocation of tRNA-mRNA complex derives from its effects on the second step of the translocation process.