Electrostatic Potential of B-DNA: Effect of Interionic Correlations

Modified Poisson-Boltzmann (MPB) equations have been numerically solved to study ionic distributions and mean electrostatic potentials around a macromolecule of arbitrarily complex shape and charge distribution. Results for DNA are compared with those obtained by classical Poisson-Boltzmann (PB) calculations. The comparisons were made for 1:1 and 2:1 electrolytes at ionic strengths up to 1 M. It is found that ion-image charge interactions and interionic correlations, which are neglected by the PB equation, have relatively weak effects on the electrostatic potential at charged groups of the DNA. The PB equation predicts errors in the long-range electrostatic part of the free energy that are only ∼1.5 kJ/mol per nucleotide even in the case of an asymmetrical electrolyte. In contrast, the spatial correlations between ions drastically affect the electrostatic potential at significant separations from the macromolecule leading to a clearly predicted effect of charge overneutralization.     

Structural dynamic aspects of protein synthesis on ribosomes

Three types of conformational changes in the translating ribosome are considered: (1) intersubunit movement (ribosome unlocking) during translocation; (2) L7/L12 stalk mobility affected by elongation factors; (3) change of tRNA residue during its transition from the A-site to the P-site. Relevant experimental data are reviewed.     

Translocation mechanism of ribosomes

The paper summarizes studies of the molecular mechanism of the dynamic function of the ribosome, i. e. translocation, performed in the author's laboratory during the past decade. The hypothesis of the locking-unlocking of the ribosomal subparticles and the kinematical model of the working ribosome, the processes of spontaneous (factor-free) and factor-dependent translocation, the sequence of events in the factor-dependent translocation, the energetics of translocation and the contribution of the elongation factors with GTP are considered. The following conclusions are made: (1) the translocation mechanism is intrinsic to the structural organization of the ribosome itself but not introduced by the protein elongation factors; (2) the transpeptidation reaction is one of the sources of energy for the work of the translocation mechanism; (3) the protein elongation factors with GTP impart additional energy to the ribosome, including that for translocation, and thus ensure excess power which is realized, in particular, in the increase of the translocation rate and its resistance against inhibitors and hindrances; (4) the promoting role of the elongation factors with GTP does not proceed by a direct conjugation of GTP hydrolysis with translocation, but through the affinity of the elongation factors to the ribosome, with a subsequent compensation of the affinity at the expense of GTP cleavage.     

Control of phosphate release from elongation factor G by ribosomal protein L7/12

Ribosomal protein L7/12 is crucial for the function of elongation factor G (EF-G) on the ribosome. Here, we report the localization of a site in the C-terminal domain (CTD) of L7/12 that is critical for the interaction with EF-G. Single conserved surface amino acids were replaced in the CTD of L7/12. Whereas mutations in helices 5 and 6 had no effect, replacements of V66, I69, K70, and R73 in helix 4 increased the Michaelis constant (KM) of EF-G.GTP for the ribosome, suggesting an involvement of these residues in EF-G binding. The mutations did not appreciably affect rapid single-round GTP hydrolysis and had no effect on tRNA translocation on the ribosome. In contrast, the release of inorganic phosphate (Pi) from ribosome-bound EF-G.GDP.Pi was strongly inhibited and became rate-limiting for the turnover of EF-G. The control of Pi release by interactions between EF-G and L7/12 appears to be important for maintaining the conformational coupling between EF-G and the ribosome for translocation and for timing the dissociation of the factor from the ribosome.     

88 Comparison of monovalent and divalent ion distributions around a DNA duplex with molecular dynamic simulation and Poisson–Boltzmann approach

Ion interactions with nucleic acids (both DNA and RNA) are an important and evolving field of investigation. Positively charged cations may interact with highly negatively charged nucleic acids via simple electrostatic interactions to help screen the electrostatic repulsion along the nucleic acids and assist their folding and/or compaction. Cations may also bind at specific sites and become integral parts of the structures, possibly playing important enzymatic roles. Two popular methods for computationally exploring a nucleic acid’s ion atmosphere are atomistic molecular dynamics (MD) simulations and the Poisson–Boltzmann (PB) equation. In general, monovalent ion results obtained from MD simulations and the PB equation agree well with experiment. However, Bai et al. (2007) observed discrepancies between experiment and the PB equation while examining the competitive binding of monovalent and divalent ions, with more significant discrepancies for divalent ions. The goal of this project was to thoroughly investigate monovalent (Na⁺) and divalent (Mg²⁺) ion distributions formed around a DNA duplex with MD simulations and the PB equation. We simulated three different cation concentrations, and matched the equilibrated bulk ion concentration for our theoretical calculations with the PB equation. Based on previous work, our Mg²⁺ ions were fully solvated, the expected state of Mg²⁺ ions when interacting with a duplex, when the production simulations began and remained throughout the simulations (Kirmizialtin, 2010; Robbins, 2012). Na⁺ ion distributions and number of Na⁺ ions within 10?Å of the DNA obtained from our two methods agreed well. However, results differed for Mg²⁺ ions, with a lower number of ions within the cut-off distance obtained from the PB equation when compared to MD simulations. The Mg²⁺ ion distributions around the DNA obtained via the two methods also differed. Based on our results, we conclude that the PB equation will systematically underestimate Mg²⁺ ions bound to DNA, and much of this deviation is attributed to dielectric saturation associated with high valency ions.     

Locking and unlocking of ribosomal motions

During the ribosomal translocation, the binding of elongation factor G (EF-G) to the pretranslocational ribosome leads to a ratchet-like rotation of the 30S subunit relative to the 50S subunit in the direction of the mRNA movement. By means of cryo-electron microscopy we observe that this rotation is accompanied by a 20 A movement of the L1 stalk of the 50S subunit, implying that this region is involved in the translocation of deacylated tRNAs from the P to the E site. These ribosomal motions can occur only when the P-site tRNA is deacylated. Prior to peptidyl-transfer to the A-site tRNA or peptide removal, the presence of the charged P-site tRNA locks the ribosome and prohibits both of these motions.     

An elongation factor G-induced ribosome rearrangement precedes tRNA-mRNA translocation

The elongation cycle of protein synthesis is completed by translocation, a rearrangement during which two tRNAs bound to the mRNA move on the ribosome. The reaction is promoted by elongation factor G (EF-G) and accelerated by GTP hydrolysis. Here we report a pre-steady-state kinetic analysis of translocation. The kinetic model suggests that GTP hydrolysis drives a conformational rearrangement of the ribosome that precedes and limits the rates of tRNA-mRNA translocation and Pi release from EF-G.GDP.Pi. The latter two steps are intrinsically rapid and take place at random. These results indicate that the energy of GTP hydrolysis is utilized to promote the ribosome rearrangement and to bias spontaneous fluctuations within the ribosome-EF-G complex toward unidirectional movement of mRNA and tRNA.     

Molecular dynamics of ribosomal elongation factors G and Tu

Translation on the ribosome is controlled by external factors. During polypeptide lengthening, elongation factors EF-Tu and EF-G consecutively interact with the bacterial ribosome. EF-Tu binds and delivers an aminoacyl-tRNA to the ribosomal A site and EF-G helps translocate the tRNAs between their binding sites after the peptide bond is formed. These processes occur at the expense of GTP. EF-Tu:tRNA and EF-G are of similar shape, share a common binding site, and undergo large conformational changes on interaction with the ribosome. To characterize the internal motion of these two elongation factors, we used 25 ns long all-atom molecular dynamics simulations. We observed enhanced mobility of EF-G domains III, IV, and V and of tRNA in the EF-Tu:tRNA complex. EF-Tu:GDP complex acquired a configuration different from that found in the crystal structure of EF-Tu with a GTP analogue, showing conformational changes in the switch I and II regions. The calculated electrostatic properties of elongation factors showed no global similarity even though matching electrostatic surface patches were found around the domain I that contacts the ribosome, and in the GDP/GTP binding region.     

Yeast translation elongation factor eEF3 promotes late stages of tRNA translocation

In addition to the conserved translation elongation factors eEF1A and eEF2, fungi require a third essential elongation factor, eEF3. While eEF3 has been implicated in tRNA binding and release at the ribosomal A and E sites, its exact mechanism of action is unclear. Here, we show that eEF3 acts at the mRNA–tRNA translocation step by promoting the dissociation of the tRNA from the E site, but independent of aminoacyl‐tRNA recruitment to the A site. Depletion of eEF3 in vivo leads to a general slowdown in translation elongation due to accumulation of ribosomes with an occupied A site. Cryo‐EM analysis of native eEF3‐ribosome complexes shows that eEF3 facilitates late steps of translocation by favoring non‐rotated ribosomal states, as well as by opening the L1 stalk to release the E‐site tRNA. Additionally, our analysis provides structural insights into novel translation elongation states, enabling presentation of a revised yeast translation elongation cycle.     

GTPases mechanisms and functions of translation factors on the ribosome

The elongation factors (EF) Tu and G and initiation factor 2 (IF2) from bacteria are multidomain GTPases with essential functions in the elongation and initiation phases of translation. They bind to the same site on the ribosome where their low intrinsic GTPase activities are strongly stimulated. The factors differ fundamentally from each other, and from the majority of GTPases, in the mechanisms of GTPase control, the timing of Pi release, and the functional role of GTP hydrolysis. EF-Tu x GTP forms a ternary complex with aminoacyl-tRNA, which binds to the ribosome. Only when a matching codon is recognized, the GTPase of EF-Tu is stimulated, rapid GTP hydrolysis and Pi release take place, EF-Tu rearranges to the GDP form, and aminoacyl-tRNA is released into the peptidyltransferase center. In contrast, EF-G hydrolyzes GTP immediately upon binding to the ribosome, stimulated by ribosomal protein L7/12. Subsequent translocation is driven by the slow dissociation of Pi, suggesting a mechano-chemical function of EF-G. Accordingly, different conformations of EF-G on the ribosome are revealed by cryo-electron microscopy. GTP hydrolysis by IF2 is triggered upon formation of the 70S initiation complex, and the dissociation of Pi and/or IF2 follows a rearrangement of the ribosome into the elongation-competent state.     

Coupling of ribosomal L1 stalk and tRNA dynamics during translation elongation

By using single-molecule fluorescence resonance energy transfer (smFRET), we observe the real-time dynamic coupling between the ribosome, labeled at the L1 stalk, and transfer RNA (tRNA). We find that an interaction between the ribosomal L1 stalk and the newly deacylated tRNA is established spontaneously upon peptide bond formation; this event involves coupled movements of the L1 stalk and tRNAs as well as ratcheting of the ribosome. In the absence of elongation factor G, the entire pretranslocation ribosome fluctuates between just two states: a nonratcheted state, with tRNAs in their classical configuration and no L1 stalk-tRNA interaction, and a ratcheted state, with tRNAs in an intermediate hybrid configuration and a direct L1 stalk-tRNA interaction. We demonstrate that binding of EF-G shifts the equilibrium toward the ratcheted state. Real-time smFRET experiments reveal that the L1 stalk-tRNA interaction persists throughout the translocation reaction, suggesting that the L1 stalk acts to direct tRNA movements during translocation.     

Hydration effects on the electrostatic potential around tuftsin

The electrostatic potential and component dielectric constants from molecular dynamics (MD) trajectories of tuftsin, a tetrapeptide with the amino acid sequence Thr–Lys–Pro–Arg in water and in saline solution are presented. The results obtained from the analysis of the MD trajectories for the total electrostatic potential at points on a grid using the Ewald technique are compared with the solution to the Poisson–Boltzmann (PB) equation. The latter was solved using several sets of dielectric constant parameters. The effects of structural averaging on the PB results were also considered. Solute conformational mobility in simulations gives rise to an electrostatic potential map around the solute dominated by the solute monopole (or lowest order multipole). The detailed spatial variation of the electrostatic potential on the molecular surface brought about by the compounded effects of the distribution of water and ions close to the peptide, solvent mobility, and solute conformational mobility are not qualitatively reproducible from a reparametrization of the input solute and solvent dielectric constants to the PB equation for a single structure or for structurally averaged PB calculations. Nevertheless, by fitting the PB to the MD electrostatic potential surfaces with the dielectric constants as fitting parameters, we found that the values that give the best fit are the values calculated from the MD trajectories. Implications of using such field calculations on the design of tuftsin peptide analogues are discussed. © 1999 John Wiley & Sons, Inc. Biopoly 50: 133–143, 1999     

Conformational changes in switch I of EF‐G drive its directional cycling on and off the ribosome

We have trapped elongation factor G (EF-G) from Escherichia coli in six, functionally defined states, representing intermediates in its unidirectional catalytic cycle, which couples GTP hydrolysis to tRNA–mRNA translocation in the ribosome. By probing EF-G with trypsin in each state, we identified a substantial conformational change involving its conserved switch I (sw1) element, which contacts the GTP substrate. By attaching FeBABE (a hydroxyl radical generating probe) to sw1, we could monitor sw1 movement (by ∼20 Å), relative to the 70S ribosome, during the EF-G cycle. In free EF-G, sw1 is disordered, particularly in GDP-bound and nucleotide-free states. On EF-G•GTP binding to the ribosome, sw1 becomes structured and tucked inside the ribosome, thereby locking GTP onto EF-G. After hydrolysis and translocation, sw1 flips out from the ribosome, greatly accelerating release of GDP and EF-G from the ribosome. Collectively, our results support a central role of sw1 in driving the EF-G cycle during protein synthesis.     

Kinetic and thermodynamic parameters for tRNA binding to the ribosome and for the translocation reaction.

Kinetic analyses of tRNA binding to the ribosome and of the translocation reaction showed the following results. 1) The activation energy for the P site binding of AcPhe-tRNA to poly(U)-programmed ribosomes is relatively high (Ea = 72 kJ mol-1; 15 mM Mg2+). If only the P site is occupied with deacylated tRNA(Phe), then the E site can be filled more easily with tRNA(Phe) (no activation energy measurable) than the A site with AcPhe-tRNA (Ea = 47 kJ mol-1; 15 mM Mg2+). 2) A ribosome with blocked P and E sites represents a standard state of the elongation cycle, in contrast to a ribosome with only a filled P site. The two states differ in that AcPhe-tRNA binding to the A site of a ribosome with prefilled P and E sites requires much higher activation energy (87 versus 47 kJ mol-1). The latter reaction simulates the allosteric transition from the post- to the pretranslocational state, whereby the tRNA(Phe) is released from the E site upon occupation of the A site (Rheinberger, H.-J., and Nierhaus, K. H. (1986) J. Biol. Chem. 261, 9133-9139). The reversed transition from the pre- to the posttranslocational state (translocation reaction) requires about the same activation energy (90 kJ mol-1). 3) Both elongation factors EF-Tu and EF-G drastically reduce the respective activation energies. 4) The rate of the A site occupation is slower than the rate of translocation in the presence of the respective elongation factors. The data suggest that the A site occupation rather than, as generally assumed, the translocation reaction is the rate-limiting step of the elongation cycle.     

Cleavage of the sarcin–ricin loop of 23S rRNA differentially affects EF-G and EF-Tu binding

Ribotoxins are potent inhibitors of protein biosynthesis and inactivate ribosomes from a variety of organisms. The ribotoxin α-sarcin cleaves the large 23S ribosomal RNA (rRNA) at the universally conserved sarcin–ricin loop (SRL) leading to complete inactivation of the ribosome and cellular death. The SRL interacts with translation factors that hydrolyze GTP, and it is important for their binding to the ribosome, but its precise role is not yet understood. We studied the effect of α-sarcin on defined steps of translation by the bacterial ribosome. α-Sarcin-treated ribosomes showed no defects in mRNA and tRNA binding, peptide-bond formation and sparsomycin-dependent translocation. Cleavage of SRL slightly affected binding of elongation factor Tu ternary complex (EF-Tu•GTP•tRNA) to the ribosome. In contrast, the activity of elongation factor G (EF-G) was strongly impaired in α-sarcin-treated ribosomes. Importantly, cleavage of SRL inhibited EF-G binding, and consequently GTP hydrolysis and mRNA–tRNA translocation. These results suggest that the SRL is more critical in EF-G than ternary complex binding to the ribosome implicating different requirements in this region of the ribosome during protein elongation.     

Allosteric interactions between the ribosomal transfer RNA-binding sites A and E

We have previously proposed a three-site model for the elongation cycle. The model is characterized by the presence of two tRNAs on the ribosome before and after translocation. We have already shown a first consequence of the model, namely that the translocation reaction is not coupled with a release of deacylated tRNA. Here we demonstrate the following conclusions. Occupation of the A site triggers the tRNA release from the E site, i.e. the A site occupation induces a drastic decrease in the affinity of the E site for deacylated tRNA. In the concentration range of deacylated tRNA in which a ribosome binds a second tRNA in addition to that one already present at the P site the deacylated tRNA does not compete for one and the same binding site with an A site ligand (AcPhe-tRNA) at 37 degrees C. It follows that the second deacylated tRNA binds to a site, the E site, which is physically distinct from the A site. When the ribosome binds a deacylated tRNA at the E site (in addition to a tRNA at the P site), the A site cannot be occupied by AcPhe-tRNA at 0 degree C and only poorly by the ternary complex elongation factor Tu . Phe-tRNA . guanyl-5'-yl imidodiphosphate. At 37 degrees C a significant A site binding is observed, with a corresponding tRNA release from the E site. In contrast, if the E site is free and only the P site occupied, the A site can bind significant amounts of charged tRNA already at 0 degree C. It follows that an occupied E site induces a low-affinity state of the A site. Thus, the ribosome always contains two high-affinity binding sites, which are A and P sites before and P and E sites after translocation. A and E sites are allosterically linked in a bidirectional manner.     

Movement of the decoding region of the 16 S ribosomal RNA accompanies tRNA translocation

The ribosome undergoes pronounced periodic conformational changes during protein synthesis. Of particular importance are those occurring around the decoding site, the region of the 16 S rRNA interacting with the mRNA-(tRNA)(2) complex. We have incorporated structural information from X-ray crystallography and nuclear magnetic resonance into cryo-electron microscopic maps of ribosomal complexes designed to capture structural changes at the translocation step of the polypeptide elongation cycle. The A-site region of the decoding site actively participates in the translocation of the tRNA from the A to the P-site upon GTP hydrolysis by elongation factor G, shifting approximately 8 A toward the P-site. This implies that elongation factor G actively pushes both the decoding site and the mRNA/tRNA complex during translocation.     

Structure of the mammalian 80S ribosome at 8.7 A resolution

In this paper, we present a structure of the mammalian ribosome determined at approximately 8.7 A resolution by electron cryomicroscopy and single-particle methods. A model of the ribosome was created by docking homology models of subunit rRNAs and conserved proteins into the density map. We then modeled expansion segments in the subunit rRNAs and found unclaimed density for approximately 20 proteins. In general, many conserved proteins and novel proteins interact with expansion segments to form an integrated framework that may stabilize the mature ribosome. Our structure provides a snapshot of the mammalian ribosome at the beginning of translation and lends support to current models in which large movements of the small subunit and L1 stalk occur during tRNA translocation. Finally, details are presented for intersubunit bridges that are specific to the eukaryotic ribosome. We suggest that these bridges may help reset the conformation of the ribosome to prepare for the next cycle of chain elongation.