Accommodation of aminoacyl-tRNA into the ribosome involves reversible excursions along multiple pathways

The Purine Bias of Coding Sequences is Determined by Physicochemical Constraints on Proteins

For this report, we analyzed protein secondary structures in relation to the statistics of three nucleotide codon positions. The purpose of this investigation was to find which properties of the ribosome, tRNA or protein level, could explain the purine bias (Rrr) as it is observed in coding DNA. We found that the Rrr pattern is the consequence of a regularity (the codon structure) resulting from physicochemical constraints on proteins and thermodynamic constraints on ribosomal machinery. The physicochemical constraints on proteins mainly come from the hydropathy and molecular weight (MW) of secondary structures as well as the energy cost of amino acid synthesis. These constraints appear through a network of statistical correlations, such as (i) the cost of amino acid synthesis, which is in favor of a higher level of guanine in the first codon position, (ii) the constructive contribution of hydropathy alternation in proteins, (iii) the spatial organization of secondary structure in proteins according to solvent accessibility, (iv) the spatial organization of secondary structure according to amino acid hydropathy, (v) the statistical correlation of MW with protein secondary structures and their overall hydropathy, (vi) the statistical correlation of thymine in the second codon position with hydropathy and the energy cost of amino acid synthesis, and (vii) the statistical correlation of adenine in the second codon position with amino acid complexity and the MW of secondary protein structures. Amino acid physicochemical properties and functional constraints on proteins constitute a code that is translated into a purine bias within the coding DNA via tRNAs. In that sense, the Rrr pattern within coding DNA is the effect of information transfer on nucleotide composition from protein to DNA by selection according to the codon positions. Thus, coding DNA structure and ribosomal machinery co-evolved to minimize the energy cost of protein coding given the functional constraints on proteins.

Sequence-dependent elongation dynamics on macrolide-bound ribosomes

The traditional view of macrolide antibiotics as plugs inside the ribosomal nascent peptide exit tunnel (NPET) has lately been challenged in favor of a more complex, heterogeneous mechanism, where drug-peptide interactions determine the fate of a translating ribosome. To investigate these highly dynamic processes, we applied single-molecule tracking of elongating ribosomes during inhibition of elongation by erythromycin of several nascent chains, including ErmCL and H-NS, which were shown to be, respectively, sensitive and resistant to erythromycin. Peptide sequence-specific changes were observed in translation elongation dynamics in the presence of a macrolide-obstructed NPET. Elongation rates were not severely inhibited in general by the presence of the drug; instead, stalls or pauses were observed as abrupt events. The dynamic pathways of nascent-chain-dependent elongation pausing in the presence of macrolides determine the fate of the translating ribosome stalling or readthrough.

Simulating movement of tRNA through the ribosome during hybrid-state formation

Biomolecular simulations provide a means for exploring the relationship between flexibility, energetics, structure, and function. With the availability of atomic models from X-ray crystallography and cryoelectron microscopy (cryo-EM), and rapid increases in computing capacity, it is now possible to apply molecular dynamics (MD) simulations to large biomolecular machines, and systematically partition the factors that contribute to function. A large biomolecular complex for which atomic models are available is the ribosome. In the cell, the ribosome reads messenger RNA (mRNA) in order to synthesize proteins. During this essential process, the ribosome undergoes a wide range of conformational rearrangements. One of the most poorly understood transitions is translocation: the process by which transfer RNA (tRNA) molecules move between binding sites inside of the ribosome. The first step of translocation is the adoption of a "hybrid" configuration by the tRNAs, which is accompanied by large-scale rotations in the ribosomal subunits. To illuminate the relationship between these rearrangements, we apply MD simulations using a multi-basin structure-based (SMOG) model, together with targeted molecular dynamics protocols. From 120 simulated transitions, we demonstrate the viability of a particular route during P/E hybrid-state formation, where there is asynchronous movement along rotation and tRNA coordinates. These simulations not only suggest an ordering of events, but they highlight atomic interactions that may influence the kinetics of hybrid-state formation. From these simulations, we also identify steric features (H74 and surrounding residues) encountered during the hybrid transition, and observe that flexibility of the single-stranded 3'-CCA tail is essential for it to reach the endpoint. Together, these simulations provide a set of structural and energetic signatures that suggest strategies for modulating the physical-chemical properties of protein synthesis by the ribosome.

Distinct functional classes of ram mutations in 16S rRNA

During decoding, the ribosome selects the correct (cognate) aminoacyl-tRNA (aa-tRNA) from a large pool of incorrect aa-tRNAs through a two-stage mechanism. In the initial selection stage, aa-tRNA is delivered to the ribosome as part of a ternary complex with elongation factor EF-Tu and GTP. Interactions between codon and anticodon lead to activation of the GTPase domain of EF-Tu and GTP hydrolysis. Then, in the proofreading stage, aa-tRNA is released from EF-Tu and either moves fully into the A/A site (a step termed "accommodation") or dissociates from the ribosome. Cognate codon-anticodon pairing not only stabilizes aa-tRNA at both stages of decoding but also stimulates GTP hydrolysis and accommodation, allowing the process to be both accurate and fast. In previous work, we isolated a number of ribosomal ambiguity (ram) mutations in 16S rRNA, implicating particular regions of the ribosome in the mechanism of decoding. Here, we analyze a representative subset of these mutations with respect to initial selection, proofreading, RF2-dependent termination, and overall miscoding in various contexts. We find that mutations that disrupt inter-subunit bridge B8 increase miscoding in a general way, causing defects in both initial selection and proofreading. Mutations in or near the A site behave differently, increasing miscoding in a codon-anticodon-dependent manner. These latter mutations may create spurious favorable interactions in the A site for certain near-cognate aa-tRNAs, providing an explanation for their context-dependent phenotypes in the cell.

Reduced model captures Mg(2+)-RNA interaction free energy of riboswitches

The stability of RNA tertiary structures depends heavily on Mg(2+). The Mg(2+)-RNA interaction free energy that stabilizes an RNA structure can be computed experimentally through fluorescence-based assays that measure Γ2+, the number of excess Mg(2+) associated with an RNA molecule. Previous explicit-solvent simulations predict that the majority of excess Mg(2+) ions interact closely and strongly with the RNA, unlike monovalent ions such as K(+), suggesting that an explicit treatment of Mg(2+) is important for capturing RNA dynamics. Here we present a reduced model that accurately reproduces the thermodynamics of Mg(2+)-RNA interactions. This model is able to characterize long-timescale RNA dynamics coupled to Mg(2+) through the explicit representation of Mg(2+) ions. KCl is described by Debye-Hückel screening and a Manning condensation parameter, which represents condensed K(+) and models its competition with condensed Mg(2+). The model contains one fitted parameter, the number of condensed K(+) ions in the absence of Mg(2+). Values of Γ2+ computed from molecular dynamics simulations using the model show excellent agreement with both experimental data on the adenine riboswitch and previous explicit-solvent simulations of the SAM-I riboswitch. This agreement confirms the thermodynamic accuracy of the model via the direct relation of Γ2+ to the Mg(2+)-RNA interaction free energy, and provides further support for the predictions from explicit-solvent calculations. This reduced model will be useful for future studies of the interplay between Mg(2+) and RNA dynamics.

Hierarchy of RNA functional dynamics

RNA dynamics play a fundamental role in many cellular functions. However, there is no general framework to describe these complex processes, which typically consist of many structural maneuvers that occur over timescales ranging from picoseconds to seconds. Here, we classify RNA dynamics into distinct modes representing transitions between basins on a hierarchical free-energy landscape. These transitions include large-scale secondary-structural transitions at >0.1-s timescales, base-pair/tertiary dynamics at microsecond-to-millisecond timescales, stacking dynamics at timescales ranging from nanoseconds to microseconds, and other "jittering" motions at timescales ranging from picoseconds to nanoseconds. We review various modes within these three different tiers, the different mechanisms by which they are used to regulate function, and how they can be coupled together to achieve greater functional complexity.

The design of a peptide sequence to inhibit HIV replication: a search algorithm combining Monte Carlo and self-consistent mean field techniques

We developed a search algorithm combining Monte Carlo (MC) and self-consistent mean field techniques to evolve a peptide sequence that has good binding capability to the anticodon stem and loop (ASL) of human lysine tRNA species, tRNA(Lys3), with the ultimate purpose of breaking the replication cycle of human immunodeficiency virus-1. The starting point is the 15-amino-acid sequence, RVTHHAFLGAHRTVG, found experimentally by Agris and co-workers to bind selectively to hypermodified tRNA(Lys3). The peptide backbone conformation is determined via atomistic simulation of the peptide-ASL(Lys3) complex and then held fixed throughout the search. The proportion of amino acids of various types (hydrophobic, polar, charged, etc.) is varied to mimic different peptide hydration properties. Three different sets of hydration properties were examined in the search algorithm to see how this affects evolution to the best-binding peptide sequences. Certain amino acids are commonly found at fixed sites for all three hydration states, some necessary for binding affinity and some necessary for binding specificity. Analysis of the binding structure and the various contributions to the binding energy shows that: 1) two hydrophilic residues (asparagine at site 11 and the cysteine at site 12) "recognize" the ASL(Lys3) due to the VDW energy, and thereby contribute to its binding specificity and 2) the positively charged arginines at sites 4 and 13 preferentially attract the negatively charged sugar rings and the phosphate linkages, and thereby contribute to the binding affinity.

Protein folding on the ribosome studied using NMR spectroscopy

NMR spectroscopy is a powerful tool for the investigation of protein folding and misfolding, providing a characterization of molecular structure, dynamics and exchange processes, across a very wide range of timescales and with near atomic resolution. In recent years NMR methods have also been developed to study protein folding as it might occur within the cell, in a de novo manner, by observing the folding of nascent polypeptides in the process of emerging from the ribosome during synthesis. Despite the 2.3 MDa molecular weight of the bacterial 70S ribosome, many nascent polypeptides, and some ribosomal proteins, have sufficient local flexibility that sharp resonances may be observed in solution-state NMR spectra. In providing information on dynamic regions of the structure, NMR spectroscopy is therefore highly complementary to alternative methods such as X-ray crystallography and cryo-electron microscopy, which have successfully characterized the rigid core of the ribosome particle. However, the low working concentrations and limited sample stability associated with ribosome-nascent chain complexes means that such studies still present significant technical challenges to the NMR spectroscopist. This review will discuss the progress that has been made in this area, surveying all NMR studies that have been published to date, and with a particular focus on strategies for improving experimental sensitivity.

Structural intermediates and folding events in the early assembly of the ribosomal small subunit

Using all-atom explicit solvent molecular dynamics (MD) simulations, we investigated the early structural intermediates of the 5' domain of the 16S rRNA in Escherichia coli upon the removal of the primary binding r-proteins S4, S17, and S20 and the secondary binding r-protein S16. Removal of each r-protein corresponded to the disappearance of subdomains with correlated dynamics. Correlation-based network analysis of the MD trajectories of the naked rRNA showed that the different subdomains are connected via multiple pathways with high betweenness. These pathways cross at the internal loop of helix 17 (h17) in the five-way junction (5WJ). The structure of the internal loop is disrupted by the binding of S17 and rescued by the addition of S16, suggesting an important function of the secondary binding protein in biasing the rRNA folding landscape toward the native basin. Using structure-based Gō simulations, we investigated the folding barriers of the lower four-way junction (4WJ) with h6, which is the primary binding site of S20 and the first to be transcribed. The folding of the 4WJ is consistent with the protection patterns observed in hydroxyl radical footprinting. Results from the all-atom simulations show that the fluctuations in the 5WJ are independent of the fluctuations in the 4WJ, suggesting that the subdomains fold independently and are stabilized by primary r-proteins.

iMODFIT: efficient and robust flexible fitting based on vibrational analysis in internal coordinates

Here, we employed the collective motions extracted from Normal Mode Analysis (NMA) in internal coordinates (torsional space) for the flexible fitting of atomic-resolution structures into electron microscopy (EM) density maps. The proposed methodology was validated using a benchmark of simulated cases, highlighting its robustness over the full range of EM resolutions and even over coarse-grained representations. A systematic comparison with other methods further showcased the advantages of this proposed methodology, especially at medium to lower resolutions. Using this method, computational costs and potential overfitting problems are naturally reduced by constraining the search in low-frequency NMA space, where covalent geometry is implicitly maintained. This method also effectively captures the macromolecular changes of a representative set of experimental test cases. We believe that this novel approach will extend the currently available EM hybrid methods to the atomic-level interpretation of large conformational changes and their functional implications.

The ribosome triggers the stringent response by RelA via a highly distorted tRNA

The bacterial stringent response links nutrient starvation with the transcriptional control of genes. This process is initiated by the stringent factor RelA, which senses the presence of deacylated tRNA in the ribosome as a symptom of amino-acid starvation to synthesize the alarmone (p)ppGpp. Here we report a cryo-EM study of RelA bound to ribosomes bearing cognate, deacylated tRNA in the A-site. The data show that RelA on the ribosome stabilizes an unusual distorted form of the tRNA, with the acceptor arm making contact with RelA and far from its normal location in the peptidyl transferase centre.

Mechanistic insights into antibiotic action on the ribosome through single-molecule fluorescence imaging

Single-molecule fluorescence imaging has provided unprecedented access to the dynamics of ribosome function, revealing transient intermediate states that are critical to ribosome activity. Imaging platforms have now been developed that are capable of probing many hundreds of molecules simultaneously at temporal and spatial resolutions approaching the sub-millisecond time and the sub-nanometer scales. These advances enable both steady- and pre-steady state measurements of individual steps in the translation process as well as processive reactions. The data generated using these methods have yielded new, quantitative structural and kinetic insights into ribosomal activity. They have also shed light on the mechanisms of antibiotic targeting the translation apparatus, revealing features of the structure-function relationship that would be difficult to obtain by other means. This review provides an overview of the types of information that can be obtained using such imaging platforms and a blueprint for using the technique to assess how small-molecule antibiotics alter macromolecular functions.

Enhanced sampling and overfitting analyses in structural refinement of nucleic acids into electron microscopy maps

Flexible fitting computational algorithms are often useful to interpret low-resolution maps of many macromolecular complexes generated by electron microscopy (EM) imaging. One such atomistic simulation technique is molecular dynamics flexible fitting (MDFF), which has been widely applied to generate structural models of large ribonucleoprotein assemblies such as the ribosome. We have previously shown that MDFF simulations of globular proteins are sensitive to the resolution of the target EM map and the strength of restraints used to preserve the secondary structure elements during fitting (Vashisth, H.; et al. Structure 2012, 20, 1453-1462). In this work, we aim to systematically examine the quality of structural models of various nucleic acids obtained via MDFF by varying the map resolution and the strength of structural restraints. We also demonstrate how an enhanced conformational sampling technique for proteins, temperature-accelerated molecular dynamics (TAMD), can be combined with MDFF for the structural refinement of nucleic acids in EM maps. Finally, we also demonstrate application of TAMD-assisted MDFF (TAMDFF) on a RNA/protein complex and suggest that TAMDFF is a viable strategy for enhanced conformational fitting in target maps of ribonucleoprotein complexes.

Connecting the kinetics and energy landscape of tRNA translocation on the ribosome

Functional rearrangements in biomolecular assemblies result from diffusion across an underlying energy landscape. While bulk kinetic measurements rely on discrete state-like approximations to the energy landscape, single-molecule methods can project the free energy onto specific coordinates. With measures of the diffusion, one may establish a quantitative bridge between state-like kinetic measurements and the continuous energy landscape. We used an all-atom molecular dynamics simulation of the 70S ribosome (2.1 million atoms; 1.3 microseconds) to provide this bridge for specific conformational events associated with the process of tRNA translocation. Starting from a pre-translocation configuration, we identified sets of residues that collectively undergo rotary rearrangements implicated in ribosome function. Estimates of the diffusion coefficients along these collective coordinates for translocation were then used to interconvert between experimental rates and measures of the energy landscape. This analysis, in conjunction with previously reported experimental rates of translocation, provides an upper-bound estimate of the free-energy barriers associated with translocation. While this analysis was performed for a particular kinetic scheme of translocation, the quantitative framework is general and may be applied to energetic and kinetic descriptions that include any number of intermediates and transition states.

Consensus among multiple approaches as a reliability measure for flexible fitting into cryo-EM data

Cryo-electron microscopy (cryo-EM) can provide low-resolution density maps of large macromolecular assemblies. As the number of structures deposited in the Protein Data Bank by fitting a high-resolution structure into a low-resolution cryo-EM map is increasing, there is a need to revise the protocols and improve the measures for fitting. A recent study suggested using a combination of multiple automated flexible fitting approaches to improve the interpretation of cryo-EM data. The current work further explores the use of multiple approaches by validating this "consensus" fitting approach and deriving a local reliability measure. Here four different flexible fitting approaches are applied for fitting an initial structure into a simulated density map of known target structure from a dataset of proteins. It is found that the models produced from different approaches often have a consensus in conformation and are also near to the target structure, whereas cases not showing consensus are away from the target. A high correlation is also observed between the RMSF profiles calculated with respect to the average and the target structures, which indicates that the relation between consensus and accuracy can also be extended to a per-residue level. Therefore, the RMSF among the fitted models is proposed as a local reliability measure, which can be used to assess the reliability of the fit at specific regions. Hence, we encourage the community to use consensus flexible fitting with different methods to report on local reliability of the resulting models and improve the interpretation of cryo-EM data.

Using enhanced sampling and structural restraints to refine atomic structures into low-resolution electron microscopy maps

For a variety of problems in structural biology, low-resolution maps generated by electron microscopy imaging are often interpreted with the help of various flexible-fitting computational algorithms. In this work, we systematically analyze the quality of final models of various proteins obtained via molecular dynamics flexible fitting (MDFF) by varying the map-resolution, strength of structural restraints, and the steering forces. We find that MDFF can be extended to understand conformational changes in lower-resolution maps if larger structural restraints and lower steering forces are used to prevent overfitting. We further show that the capabilities of MDFF can be extended by combining it with an enhanced conformational sampling method, temperature-accelerated molecular dynamics (TAMD). Specifically, either TAMD can be used to generate better starting configurations for MDFF fitting or TAMD-assisted MDFF (TAMDFF) can be performed to accelerate conformational search in atomistic simulations.

Single-molecule nanometry for biological physics

Precision measurement is a hallmark of physics but the small length scale (∼nanometer) of elementary biological components and thermal fluctuations surrounding them challenge our ability to visualize their action. Here, we highlight the recent developments in single-molecule nanometry where the position of a single fluorescent molecule can be determined with nanometer precision, reaching the limit imposed by the shot noise, and the relative motion between two molecules can be determined with ∼0.3 nm precision at ∼1 ms time resolution, as well as how these new tools are providing fundamental insights into how motor proteins move on cellular highways. We will also discuss how interactions between three and four fluorescent molecules can be used to measure three and six coordinates, respectively, allowing us to correlate the movements of multiple components. Finally, we will discuss recent progress in combining angstrom-precision optical tweezers with single-molecule fluorescent detection, opening new windows for multi-dimensional single-molecule nanometry for biological physics.

Allosteric control in a metalloprotein dramatically alters function

Metalloproteins (MPs) comprise one-third of all known protein structures. This diverse set of proteins contain a plethora of unique inorganic moieties capable of performing chemistry that would otherwise be impossible using only the amino acids found in nature. Most of the well-studied MPs are generally viewed as being very rigid in structure, and it is widely thought that the properties of the metal centers are primarily determined by the small fraction of amino acids that make up the local environment. Here we examine both theoretically and experimentally whether distal regions can influence the metal center in the diabetes drug target mitoNEET. We demonstrate that a loop (L2) 20 Å away from the metal center exerts allosteric control over the cluster binding domain and regulates multiple properties of the metal center. Mutagenesis of L2 results in significant shifts in the redox potential of the [2Fe-2S] cluster and orders of magnitude effects on the rate of [2Fe-2S] cluster transfer to an apo-acceptor protein. These surprising effects occur in the absence of any structural changes. An examination of the native basin dynamics of the protein using all-atom simulations shows that twisting in L2 controls scissoring in the cluster binding domain and results in perturbations to one of the cluster-coordinating histidines. These allosteric effects are in agreement with previous folding simulations that predicted L2 could communicate with residues surrounding the metal center. Our findings suggest that long-range dynamical changes in the protein backbone can have a significant effect on the functional properties of MPs.

Protein dynamics via computational microscope

The purpose of this overview is to provide a concise introduction to the methodology and current advances in molecular dynamics (MD) simulations. MD simulations emerged as a powerful and popular tool to study dynamic behavior of proteins and macromolecule complexes at the atomic resolution. This approach can extend static structural data, such as X-ray crystallography, into dynamic domains with realistic timescales (up to millisecond) and high precision, therefore becoming a veritable computational microscope. This perspective covers current advances and methodology in the simulation of protein folding and drug design as illustrated by several important published examples. Overall, recent progress in the simulation field points to the direction that MD will have significant impact on molecular biology and pharmaceutical science.

The origin of minus-end directionality and mechanochemistry of Ncd motors

Adaptation of molecular structure to the ligand chemistry and interaction with the cytoskeletal filament are key to understanding the mechanochemistry of molecular motors. Despite the striking structural similarity with kinesin-1, which moves towards plus-end, Ncd motors exhibit minus-end directionality on microtubules (MTs). Here, by employing a structure-based model of protein folding, we show that a simple repositioning of the neck-helix makes the dynamics of Ncd non-processive and minus-end directed as opposed to kinesin-1. Our computational model shows that Ncd in solution can have both symmetric and asymmetric conformations with disparate ADP binding affinity, also revealing that there is a strong correlation between distortion of motor head and decrease in ADP binding affinity in the asymmetric state. The nucleotide (NT) free-ADP (φ-ADP) state bound to MTs favors the symmetric conformation whose coiled-coil stalk points to the plus-end. Upon ATP binding, an enhanced flexibility near the head-neck junction region, which we have identified as the important structural element for directional motility, leads to reorienting the coiled-coil stalk towards the minus-end by stabilizing the asymmetric conformation. The minus-end directionality of the Ncd motor is a remarkable example that demonstrates how motor proteins in the kinesin superfamily diversify their functions by simply rearranging the structural elements peripheral to the catalytic motor head domain.

Proteins belonging to the kinesin superfamily are responsible for vesicle or organelle transport, spindle morphogenesis, and chromosome sorting during cell division. Interestingly, while most proteins in kinesin superfamily that share the common catalytic motor head domain have plus-end directionality along microtubules, kinesin-14 (Ncd) exhibits minus-end directionality. Despite the several circumstantial evidences on the determining factors for the motor directionality in the last decade, a comprehensive understanding of the mechanism governing the Ncd minus-end directionality is still missing. Our studies provide a clear explanation for this minus-end directionality and the associated mechanochemical cycle. Here, we modeled an Ncd motor by employing structural details available in the literature to simulate its conformational dynamics. Simulations using our structure-based model of Ncd assert that the dynamics due to a simple rearrangement of structural elements, peripheral to the catalytic motor domain, is the key player in determining both the directionality and mechanochemistry unique to Ncd motors. Although Ncd has a different directionality, it uses a similar strategy to kinesin-1 of structural adaptation of the catalytic motor domain. Therefore using the same physical principle of protein folding and very similar structural elements, motors in the kinesin superfamily are able to achieve a variety of biological function.

The dominant folding route minimizes backbone distortion in SH3

Energetic frustration in protein folding is minimized by evolution to create a smooth and robust energy landscape. As a result the geometry of the native structure provides key constraints that shape protein folding mechanisms. Chain connectivity in particular has been identified as an essential component for realistic behavior of protein folding models. We study the quantitative balance of energetic and geometrical influences on the folding of SH3 in a structure-based model with minimal energetic frustration. A decomposition of the two-dimensional free energy landscape for the folding reaction into relevant energy and entropy contributions reveals that the entropy of the chain is not responsible for the folding mechanism. Instead the preferred folding route through the transition state arises from a cooperative energetic effect. Off-pathway structures are penalized by excess distortion in local backbone configurations and contact pair distances. This energy cost is a new ingredient in the malleable balance of interactions that controls the choice of routes during protein folding.

Unraveling the dynamics of ribosome translocation

Translocation is one of the key events in translation, requiring large-scale conformational changes in the ribosome, movements of two transfer RNAs (tRNAs) across a distance of more than 20Å, and the coupled movement of the messenger RNA (mRNA) by one codon, completing one cycle of peptide-chain elongation. Translocation is catalyzed by elongation factor G (EF-G in bacteria), which hydrolyzes GTP in the process. However, how the conformational rearrangements of the ribosome actually drive the movements of the tRNAs and how EF-G GTP hydrolysis plays a role in this process are still unclear. Fluorescence methods, both single-molecule and bulk, have provided a dynamic view of translocation, allowing us to follow the different conformational changes of the ribosome in real-time. The application of electron microscopy has revealed new conformational intermediates during translocation and important structural rearrangements in the ribosome that drive tRNA movement, while computational approaches have added quantitative views of the translational pathway. These recent advances shed light on the process of translocation, providing insight on how to resolve the different descriptions of translocation in the current literature.

Substrate-specific reorganization of the conformational ensemble of CSK implicates novel modes of kinase function

Protein kinases use ATP as a phosphoryl donor for the posttranslational modification of signaling targets. It is generally thought that the binding of this nucleotide induces conformational changes leading to closed, more compact forms of the kinase domain that ideally orient active-site residues for efficient catalysis. The kinase domain is oftentimes flanked by additional ligand binding domains that up- or down-regulate catalytic function. C-terminal Src kinase (Csk) is a multidomain tyrosine kinase that is up-regulated by N-terminal SH2 and SH3 domains. Although the X-ray structure of Csk suggests the enzyme is compact, X-ray scattering studies indicate that the enzyme possesses both compact and open conformational forms in solution. Here, we investigated whether interactions with the ATP analog AMP-PNP and ADP can shift the conformational ensemble of Csk in solution using a combination of small angle x-ray scattering and molecular dynamics simulations. We find that binding of AMP-PNP shifts the ensemble towards more extended rather than more compact conformations. Binding of ADP further shifts the ensemble towards extended conformations, including highly extended conformations not adopted by the apo protein, nor by the AMP-PNP bound protein. These ensembles indicate that any compaction of the kinase domain induced by nucleotide binding does not extend to the overall multi-domain architecture. Instead, assembly of an ATP-bound kinase domain generates further extended forms of Csk that may have relevance for kinase scaffolding and Src regulation in the cell.

The Src protein kinases are integral members of numerous signaling pathways involved in cellular growth and differentiation. The master regulator of the Src family is the protein kinase Csk, which adds a phosphate to the C-terminal tail, inhibiting Src Kinase function. Proper regulation of these signaling pathways by Csk is essential as unregulated activity in these pathways is correlated with the development of various cancers and autoimmune diseases. Understanding the nature of the mechanism and structure of Csk may lead to therapeutics and a better understanding of Src signaling pathways. Conformational changes associated with nucleotide binding and release have been shown to regulate the efficiency of Src down-regulation by Csk. To obtain insights into the nature of these nucleotide-induced structural changes, we examined the conformation of Csk in solution while bound to the ATP analog AMP-PNP and product ADP using a combination of small angle x-ray scattering and molecular dynamics. Surprisingly, both nucleotides induce extended conformations of Csk compared to the apo-enzyme, suggesting a novel mode of function. Further understanding of this mode of function may aid in the design of cancer therapeutics that act by regulating Src signaling pathways by modulating the function of Csk.

Bridging the gap between ribosome structure and biochemistry by mechanistic computations

The wealth of structural and biochemical data now available for protein synthesis on the ribosome presents major new challenges for computational biochemistry. Apart from technical difficulties in modeling ribosome systems, the complexity of the overall translation cycle with a multitude of different kinetic steps presents a formidable problem for computational efforts where we have only seen the beginning. However, a range of methodologies including molecular dynamics simulations, free energy calculations, molecular docking and quantum chemical approaches have already been put to work with promising results. In particular, the combined efforts of structural biology, biochemistry, kinetics and computational modeling can lead towards a quantitative structure-based description of translation.

Mechanisms of SecM-mediated stalling in the ribosome

Nascent-peptide modulation of translation is a common regulatory mechanism of gene expression. In this mechanism, while the nascent peptide is still in the exit tunnel of the ribosome, it induces translational pausing, thereby controlling the expression of downstream genes. One example is SecM, which inhibits peptide-bond formation in the ribosome's peptidyl transferase center (PTC) during its own translation, upregulating the expression of the protein translocase SecA. Although biochemical experiments and cryo-electron microscopy data have led to the identification of some residues involved in SecM recognition, the full pathway of interacting residues that connect SecM to the PTC through the ribosome has not yet been conclusively established. Here, using the cryo-electron microscopy data, we derived the first (to our knowledge) atomic model of the SecM-stalled ribosome via molecular-dynamics flexible fitting, complete with P- and A-site tRNAs. Subsequently, we carried out simulations of native and mutated SecM-stalled ribosomes to investigate possible interaction pathways between a critical SecM residue, R163, and the PTC. In particular, the simulations reveal the role of SecM in altering the position of the tRNAs in the ribosome, and thus demonstrate how the presence of SecM in the exit tunnel induces stalling. Finally, steered molecular-dynamics simulations in which SecM was pulled toward the tunnel exit suggest how SecA interacting with SecM from outside the ribosome relieves stalling.

Structure and mechanical properties of the ribosomal L1 stalk three-way junction

The L1 stalk is a key mobile element of the large ribosomal subunit which interacts with tRNA during translocation. Here, we investigate the structure and mechanical properties of the rRNA H76/H75/H79 three-way junction at the base of the L1 stalk from four different prokaryotic organisms. We propose a coarse-grained elastic model and parameterize it using large-scale atomistic molecular dynamics simulations. Global properties of the junction are well described by a model in which the H76 helix is represented by a straight, isotropically flexible elastic rod, while the junction core is represented by an isotropically flexible spherical hinge. Both the core and the helix contribute substantially to the overall H76 bending fluctuations. The presence of wobble pairs in H76 does not induce any increased flexibility or anisotropy to the helix. The half-closed conformation of the L1 stalk seems to be accessible by thermal fluctuations of the junction itself, without any long-range allosteric effects. Bending fluctuations of H76 with a bulge introduced in it suggest a rationale for the precise position of the bulge in eukaryotes. Our elastic model can be generalized to other RNA junctions found in biological systems or in nanotechnology.

Biomolecular dynamics: order-disorder transitions and energy landscapes

While the energy landscape theory of protein folding is now a widely accepted view for understanding how relatively weak molecular interactions lead to rapid and cooperative protein folding, such a framework must be extended to describe the large-scale functional motions observed in molecular machines. In this review, we discuss (1) the development of the energy landscape theory of biomolecular folding, (2) recent advances toward establishing a consistent understanding of folding and function and (3) emerging themes in the functional motions of enzymes, biomolecular motors and other biomolecular machines. Recent theoretical, computational and experimental lines of investigation have provided a very dynamic picture of biomolecular motion. In contrast to earlier ideas, where molecular machines were thought to function similarly to macroscopic machines, with rigid components that move along a few degrees of freedom in a deterministic fashion, biomolecular complexes are only marginally stable. Since the stabilizing contribution of each atomic interaction is on the order of the thermal fluctuations in solution, the rigid body description of molecular function must be revisited. An emerging theme is that functional motions encompass order-disorder transitions and structural flexibility provides significant contributions to the free energy. In this review, we describe the biological importance of order-disorder transitions and discuss the statistical-mechanical foundation of theoretical approaches that can characterize such transitions.

Assembly of the five-way junction in the ribosomal small subunit using hybrid MD-Gō simulations

Assembly of the bacterial ribosomal small subunit (SSU) begins with the folding of the five-way junction upon interaction with the primary binding protein S4. This complex contains the largest contiguous molecular signature, which is a conserved feature in all bacterial 16S rRNAs. In a previous study, we used all-atom molecular dynamics simulations to demonstrate that the co-evolving signature in the N-terminus of S4 is intrinsically disordered and capable of accelerating the binding process through a fly casting mechanism. In this paper, comparisons between the all-atom MD simulations and FRET experiments identify multiple metastable conformations of the naked five-way junction without the presence of S4. Furthermore, we capture the simultaneous folding and binding of the five-way junction and r-protein S4 by use of a structure-based Gō potential implemented within the framework of the all-atom molecular dynamics CHARMM force field. Different folding pathways are observed for the refolding of the five-way junction upon partial binding of S4. Our simulations illustrate the complex nature of RNA folding in the presence of a protein binding partner and provide insight into the role of population shift and the induced fit mechanisms in the protein:RNA folding and binding process.

A computational investigation on the connection between dynamics properties of ribosomal proteins and ribosome assembly

Assembly of the ribosome from its protein and RNA constituents has been studied extensively over the past 50 years, and experimental evidence suggests that prokaryotic ribosomal proteins undergo conformational changes during assembly. However, to date, no studies have attempted to elucidate these conformational changes. The present work utilizes computational methods to analyze protein dynamics and to investigate the linkage between dynamics and binding of these proteins during the assembly of the ribosome. Ribosomal proteins are known to be positively charged and we find the percentage of positive residues in r-proteins to be about twice that of the average protein: Lys+Arg is 18.7% for E. coli and 21.2% for T. thermophilus. Also, positive residues constitute a large proportion of RNA contacting residues: 39% for E. coli and 46% for T. thermophilus. This affirms the known importance of charge-charge interactions in the assembly of the ribosome. We studied the dynamics of three primary proteins from E. coli and T. thermophilus 30S subunits that bind early in the assembly (S15, S17, and S20) with atomic molecular dynamic simulations, followed by a study of all r-proteins using elastic network models. Molecular dynamics simulations show that solvent-exposed proteins (S15 and S17) tend to adopt more stable solution conformations than an RNA-embedded protein (S20). We also find protein residues that contact the 16S rRNA are generally more mobile in comparison with the other residues. This is because there is a larger proportion of contacting residues located in flexible loop regions. By the use of elastic network models, which are computationally more efficient, we show that this trend holds for most of the 30S r-proteins.

Ribosomes are complex cellular machines that synthesize new proteins in the cell. The accurate and efficient assembly of ribosomal proteins (r-proteins) and ribosomal RNA (rRNA) to form a functional ribosome is important for cell growth, metabolic reactions, and other cellular processes. Additionally, some antibacterial drugs are believed to target the bacterial ribosome during its construction. Hence, ribosomal assembly has been an active research topic for many years because understanding the assembly mechanisms can provide insight into protein/RNA recognitions important in many other cellular processes, as well as optimize the development of antibacterial therapeutics. Experimental studies thus far have provided still limited understanding about the assembly process. To further understand the assembly process, we have computationally studied the dynamic properties that r-proteins exhibit during assembly and the relationship between dynamics, physical properties, and binding propensity. We observe significant charged interactions between r-proteins and rRNA. We also detect a strong correlation between contact residues and their dynamic mobilities. Protein residues contacting with rRNA are observed to be more mobile in comparison with other residues. We also relate the location of the r-protein in the fully assembled ribosome to its susceptibility for large conformational changes prior to binding.

The shadow map: a general contact definition for capturing the dynamics of biomolecular folding and function

Structure-based models (SBMs) are simplified models of the biomolecular dynamics that arise from funneled energy landscapes. We recently introduced an all-atom SBM that explicitly represents the atomic geometry of a biomolecule. While this initial study showed the robustness of the all-atom SBM Hamiltonian to changes in many of the energetic parameters, an important aspect, which has not been explored previously, is the definition of native interactions. In this study, we propose a general definition for generating atomically grained contact maps called "Shadow". The Shadow algorithm initially considers all atoms within a cutoff distance and then, controlled by a screening parameter, discards the occluded contacts. We show that this choice of contact map is not only well behaved for protein folding, since it produces consistently cooperative folding behavior in SBMs but also desirable for exploring the dynamics of macromolecular assemblies since, it distributes energy similarly between RNAs and proteins despite their disparate internal packing. All-atom structure-based models employing Shadow contact maps provide a general framework for exploring the geometrical features of biomolecules, especially the connections between folding and function.

A new understanding of the decoding principle on the ribosome

During protein synthesis, the ribosome accurately selects transfer RNAs (tRNAs) in accordance with the messenger RNA (mRNA) triplet in the decoding centre. tRNA selection is initiated by elongation factor Tu, which delivers tRNA to the aminoacyl tRNA-binding site (A site) and hydrolyses GTP upon establishing codon-anticodon interactions in the decoding centre. At the following proofreading step the ribosome re-examines the tRNA and rejects it if it does not match the A codon. It was suggested that universally conserved G530, A1492 and A1493 of 16S ribosomal RNA, critical for tRNA binding in the A site, actively monitor cognate tRNA, and that recognition of the correct codon-anticodon duplex induces an overall ribosome conformational change (domain closure). Here we propose an integrated mechanism for decoding based on six X-ray structures of the 70S ribosome determined at 3.1-3.4 Å resolution, modelling cognate or near-cognate states of the decoding centre at the proofreading step. We show that the 30S subunit undergoes an identical domain closure upon binding of either cognate or near-cognate tRNA. This conformational change of the 30S subunit forms a decoding centre that constrains the mRNA in such a way that the first two nucleotides of the A codon are limited to form Watson-Crick base pairs. When U·G and G·U mismatches, generally considered to form wobble base pairs, are at the first or second codon-anticodon position, the decoding centre forces this pair to adopt the geometry close to that of a canonical C·G pair. This by itself, or with distortions in the codon-anticodon mini-helix and the anticodon loop, causes the near-cognate tRNA to dissociate from the ribosome.

Computational studies of molecular machines: the ribosome

The past decade has produced an avalanche of experimental data on the structure and dynamics of the ribosome. Groundbreaking studies in structural biology and kinetics have placed important constraints on ribosome structural dynamics. However, a gulf remains between static structures and time dependent data. In particular, X-ray crystallography and cryo-EM studies produce static models of the ribosome in various states, but lack dynamic information. Single molecule studies produce information on the rates of transitions between these states but do not have high-resolution spatial information. Computational studies have aided in bridging this gap by providing atomic resolution simulations of structural fluctuations and transitions between configurations.

Selenocysteine insertion sequence (SECIS)-binding protein 2 alters conformational dynamics of residues involved in tRNA accommodation in 80 S ribosomes

Sec-tRNA(Sec) is site-specifically delivered at defined UGA codons in selenoprotein mRNAs. This recoding event is specified by the selenocysteine insertion sequence (SECIS) element and requires the selenocysteine (Sec)-specific elongation factor, eEFSec, and the SECIS binding protein, SBP2. Sec-tRNA(Sec) is delivered to the ribosome by eEFSec-GTP, but this ternary complex is not sufficient for Sec incorporation, indicating that its access to the ribosomal A-site is regulated. SBP2 stably associates with ribosomes, and mutagenic analysis indicates that this interaction is essential for Sec incorporation. However, the ribosomal function of SBP2 has not been elucidated. To shed light on the functional relevance of the SBP2-ribosome interaction, we screened the functional centers of the 28 S rRNA in translationally competent 80 S ribosomes using selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE). We demonstrate that SBP2 specifically alters the reactivity of specific residues in Helix 89 (H89) and expansion segment 31 (ES31). These results are indicative of a conformational change in response to SBP2 binding. Based on the known functions of H89 during translation, we propose that SBP2 allows Sec incorporation by either promoting Sec-tRNA(Sec) accommodation into the peptidyltransferase center and/or by stimulating the ribosome-dependent GTPase activity of eEFSec.

Evolutionary optimization of speed and accuracy of decoding on the ribosome

Speed and accuracy of protein synthesis are fundamental parameters for the fitness of living cells, the quality control of translation, and the evolution of ribosomes. The ribosome developed complex mechanisms that allow for a uniform recognition and selection of any cognate aminoacyl-tRNA (aa-tRNA) and discrimination against any near-cognate aa-tRNA, regardless of the nature or position of the mismatch. This review describes the principles of the selection-kinetic partitioning and induced fit-and discusses the relationship between speed and accuracy of decoding, with a focus on bacterial translation. The translational machinery apparently has evolved towards high speed of translation at the cost of fidelity.

A-site residues move independently from P-site residues in all-atom molecular dynamics simulations of the 70S bacterial ribosome

The ribosome is a large macromolecular machine, and correlated motion between residues is necessary for coordinating function across multiple protein and RNA chains. We ran two all-atom, explicit solvent molecular dynamics simulations of the bacterial ribosome and calculated correlated motion between residue pairs by using mutual information. Because of the short timescales of our simulation (ns), we expect that dynamics are largely local fluctuations around the crystal structure. We hypothesize that residues that show coupled dynamics are functionally related, even on longer timescales. We validate our model by showing that crystallographic B-factors correlate well with the entropy calculated as part of our mutual information calculations. We reveal that A-site residues move relatively independently from P-site residues, effectively insulating A-site functions from P-site functions during translation.

Quality control of mRNA decoding on the bacterial ribosome

The ribosome is a major player in providing accurate gene expression in the cell. The fidelity of substrate selection is tightly controlled throughout the translation process, including the initiation, elongation, and termination phases. Although each phase of translation involves different players, that is, translation factors and tRNAs, the general principles of selection appear surprisingly similar for very different substrates. At essentially every step of translation, differences in complex stabilities as well as induced fit are sources of selectivity. A view starts to emerge of how the ribosome uses local and global conformational switches to govern induced-fit mechanisms that ensure fidelity. This review describes the mechanisms of tRNA and mRNA selection at all phases of protein synthesis in bacteria.

The future of molecular dynamics simulations in drug discovery

Molecular dynamics simulations can now track rapid processes—those occurring in less than about a millisecond—at atomic resolution for many biologically relevant systems. These simulations appear poised to exert a significant impact on how new drugs are found, perhaps even transforming the very process of drug discovery. We predict here future results we can expect from, and enhancements we need to make in, molecular dynamics simulations over the coming 25 years, and in so doing set out several Grand Challenges for the field. In the context of the problems now facing the pharmaceutical industry, we ask how we can best address drug discovery needs of the next quarter century using molecular dynamics simulations, and we suggest some possible approaches.

In vitro and in vivo single-molecule fluorescence imaging of ribosome-catalyzed protein synthesis

Combined with the availability of highly purified, fluorescently labeled in vitro translation systems, the advent of single-molecule fluorescence imaging has ushered in a new era in high-resolution mechanistic studies of ribosome-catalyzed protein synthesis, or translation. Together with ensemble biochemical investigations of translation and structural studies of functional ribosomal complexes, in vitro single-molecule fluorescence imaging of protein synthesis is providing unique mechanistic insight into this fundamental biological process. More recently, rapidly evolving breakthroughs in fluorescence-based molecular imaging in live cells with sub-diffraction-limit spatial resolution and ever-increasing temporal resolution provide great promise for conducting mechanistic studies of translation and its regulation in living cells. Here we review the remarkable recent progress that has been made in these fields, highlight important mechanistic insights that have been gleaned from these studies thus far, and discuss what we envision lies ahead as these approaches continue to evolve and expand to address increasingly complex mechanistic and regulatory aspects of translation.

Excited states of ribosome translocation revealed through integrative molecular modeling

The dynamic nature of biomolecules leads to significant challenges when characterizing the structural properties associated with function. While X-ray crystallography and imaging techniques (such as cryo-electron microscopy) can reveal the structural details of stable molecular complexes, strategies must be developed to characterize configurations that exhibit only marginal stability (such as intermediates) or configurations that do not correspond to minima on the energy landscape (such as transition-state ensembles). Here, we present a methodology (MDfit) that utilizes molecular dynamics simulations to generate configurations of excited states that are consistent with available biophysical and biochemical measurements. To demonstrate the approach, we present a sequence of configurations that are suggested to be associated with transfer RNA (tRNA) movement through the ribosome (translocation). The models were constructed by combining information from X-ray crystallography, cryo-electron microscopy, and biochemical data. These models provide a structural framework for translocation that may be further investigated experimentally and theoretically to determine the precise energetic character of each configuration and the transition dynamics between them.

A-minor tertiary interactions in RNA kink-turns. Molecular dynamics and quantum chemical analysis

The RNA kink-turn is an important recurrent RNA motif, an internal loop with characteristic consensus sequence forming highly conserved three-dimensional structure. Functional arrangement of RNA kink-turns shows a sharp bend in the phosphodiester backbone. Among other signature interactions, kink-turns form A-minor interaction between their two stems. Most kink-turns possess extended A-minor I (A-I) interaction where adenine of the second A•G base pair of the NC-stem interacts with the first canonical pair of the C-stem (i.e., the receptor pair) via trans-sugar-edge/sugar-edge (tSS) and cis-sugar-edge/sugar-edge (cSS) interactions. The remaining kink-turns have less compact A-minor 0 (A-0) interaction with just one tSS contact. We show that kink-turns with A-I in ribosomal X-ray structures keep G═C receptor base pair during evolution while the inverted pair (C═G) is not realized. In contrast, kink-turns with A-0 in the observed structures alternate G═C and C═G base pairs in sequences. We carried out an extended set (~5 μs) of explicit-solvent molecular dynamics simulations of kink-turns to rationalize this structural/evolutionary pattern. The simulations were done using a net-neutral Na(+) cation atmosphere (with ~0.25 M cation concentration) supplemented by simulations with either excess salt KCl atmosphere or inclusion of Mg(2+). The results do not seem to depend on the treatment of ions. The simulations started with X-ray structures of several kink-turns while we tested the response of the simulated system to base substitutions, modest structural perturbations and constraints. The trends seen in the simulations reveal that the A-I/G═C arrangement is preferred over all three other structures. The A-I/C═G triple appears structurally entirely unstable, consistent with the covariation patterns seen during the evolution. The A-0 arrangements tend to shift toward the A-I pattern in simulations, which suggests that formation of the A-0 interaction is likely supported by the surrounding protein and RNA molecules. A-0 may also be stabilized by additional kink-turn nucleotides not belonging to the kink-turn consensus, as shown for the kink-turn from ribosomal Helix 15. Quantum-chemical calculations on all four A-minor triples suggest that there is a different balance of electrostatic and dispersion stabilization in the A-I/G═C and A-I/C═G triples, which may explain different behavior of these otherwise isosteric triples in the context of kink-turns.

Consensus among flexible fitting approaches improves the interpretation of cryo-EM data

Cryo-elecron microscopy (cryo-EM) can provide important structural information of large macromolecular assemblies in different conformational states. Recent years have seen an increase in structures deposited in the Protein Data Bank (PDB) by fitting a high-resolution structure into its low-resolution cryo-EM map. A commonly used protocol for accommodating the conformational changes between the X-ray structure and the cryo-EM map is rigid body fitting of individual domains. With the emergence of different flexible fitting approaches, there is a need to compare and revise these different protocols for the fitting. We have applied three diverse automated flexible fitting approaches on a protein dataset for which rigid domain fitting (RDF) models have been deposited in the PDB. In general, a consensus is observed in the conformations, which indicates a convergence from these theoretically different approaches to the most probable solution corresponding to the cryo-EM map. However, the result shows that the convergence might not be observed for proteins with complex conformational changes or with missing densities in cryo-EM map. In contrast, RDF structures deposited in the PDB can represent conformations that not only differ from the consensus obtained by flexible fitting but also from X-ray crystallography. Thus, this study emphasizes that a "consensus" achieved by the use of several automated flexible fitting approaches can provide a higher level of confidence in the modeled configurations. Following this protocol not only increases the confidence level of fitting, but also highlights protein regions with uncertain fitting. Hence, this protocol can lead to better interpretation of cryo-EM data.

Key intermolecular interactions in the E. coli 70S ribosome revealed by coarse-grained analysis

The ribosome is a very large complex that consists of many RNA and protein molecules and plays a central role in protein biosynthesis in all organisms. Extensive interactions between different molecules are critical to ribosomal functional dynamics. In this work, intermolecular interactions in the Escherichia coli 70S ribosome are investigated by coarse-grained (CG) analysis. CG models are defined to preserve dynamic domains in RNAs and proteins and to capture functional motions in the ribosome, and then the CG sites are connected by harmonic springs, and spring constants are obtained by matching the computed fluctuations to those of an all-atom molecular dynamics (MD) simulation. Those spring constants indicate how strong the interactions are between the ribosomal components, and they are in good agreement with various experimental data. Nearly all the bridges between the small and large ribosomal subunits are indicated by CG interactions with large spring constants. The head of the small subunit is very mobile because it has minimal CG interactions with the rest of the subunit; however, a large number of small subunit proteins bind to maintain the internal structure of the head. The results show a clear connection between the intermolecular interactions and the structural and functional properties of the ribosome because of the reduced complexity in domain-based CG models. The present approach also provides a useful strategy to map interactions between molecules within large biomolecular complexes since it is not straightforward to investigate these by either atomistic MD simulations or residue-based elastic network models.

Mutations of highly conserved bases in the peptidyltransferase center induce compensatory rearrangements in yeast ribosomes

Molecular dynamics simulation identified three highly conserved rRNA bases in the large subunit of the ribosome that form a three-dimensional (3D) "gate" that induces pausing of the aa-tRNA acceptor stem during accommodation into the A-site. A nearby fourth base contacting the "tryptophan finger" of yeast protein L3, which is involved in the coordinating elongation factor recruitment to the ribosome with peptidyltransfer, is also implicated in this process. To better understand the functional importance of these bases, single base substitutions as well as deletions at all four positions were constructed and expressed as the sole forms of ribosomes in yeast Saccharomyces cerevisiae. None of the mutants had strong effects on cell growth, translational fidelity, or on the interactions between ribosomes and tRNAs. However, the mutants did promote strong effects on cell growth in the presence of translational inhibitors, and differences in viability between yeast and Escherichia coli mutants at homologous positions suggest new targets for antibacterial therapeutics. Mutant ribosomes also promoted changes in 25S rRNA structure, all localized to the core of peptidyltransferase center (i.e., the proto-ribosome area). We suggest that a certain degree of structural plasticity is built into the ribosome, enabling it to ensure accurate translation of the genetic code while providing it with the flexibility to adapt and evolve.

Nascent peptide in the ribosome exit tunnel affects functional properties of the A-site of the peptidyl transferase center

The ability to monitor the nascent peptide structure and to respond functionally to specific nascent peptide sequences is a fundamental property of the ribosome. An extreme manifestation of such response is nascent peptide-dependent ribosome stalling, involved in the regulation of gene expression. The molecular mechanisms of programmed translation arrest are unclear. By analyzing ribosome stalling at the regulatory cistron of the antibiotic resistance gene ermA, we uncovered a carefully orchestrated cooperation between the ribosomal exit tunnel and the A-site of the peptidyl transferase center (PTC) in halting translation. The presence of an inducing antibiotic and a specific nascent peptide in the exit tunnel abrogate the ability of the PTC to catalyze peptide bond formation with a particular subset of amino acids. The extent of the conferred A-site selectivity is modulated by the C-terminal segment of the nascent peptide, where the third-from-last residue plays a critical role.