Molecular simulations of the ribosome and associated translation factors

The role of ribosomal protein networks in ribosome dynamics

Accurate protein synthesis requires ribosomes to integrate signals from distant functional sites and execute complex dynamics. Despite advances in understanding ribosome structure and function, two key questions remain: how information is transmitted between these distant sites, and how ribosomal movements are synchronized? We recently highlighted the existence of ribosomal protein networks, likely evolved to participate in ribosome signaling. Here, we investigate the relationship between ribosomal protein networks and ribosome dynamics. Our findings show that major motion centers in the bacterial ribosome interact specifically with r-proteins, and that ribosomal RNA exhibits high mobility around each r-protein. This suggests that periodic electrostatic changes in the context of negatively charged residues (Glu and Asp) induce RNA-protein 'distance-approach' cycles, controlling key ribosomal movements during translocation. These charged residues play a critical role in modulating electrostatic repulsion between RNA and proteins, thus coordinating ribosomal dynamics. We propose that r-protein networks synchronize ribosomal dynamics through an 'electrostatic domino' effect, extending the concept of allostery to the regulation of movements within supramolecular assemblies.

Stabilization Mechanism of Initiator Transfer RNA in the Small Ribosomal Subunit from Coarse-Grained Molecular Simulations

Proteins play a variety of roles in biological phenomena in cells. Proteins are synthesized by the ribosome, which is a large molecular complex composed of proteins and nucleic acids. Among the many molecules involved in the process of protein synthesis, tRNA is one of the essential molecules. In this study, coarse-grained molecular dynamics simulations were performed to understand how the tRNA molecule is stabilized in the ribosome, and the free energy along the dissociation path of the tRNA was calculated. We found that some ribosomal proteins, which are components of the ribosome, are involved in the stabilization of the tRNA. The positively charged amino acid residues in the C-terminal region of the ribosomal proteins are particularly important for stabilization. These findings contribute to our understanding of the molecular evolution of protein synthesis in terms of the ribosome, which is a universal component of life.

Conformational changes of ribosomes during translation elongation resolved by molecular dynamics simulations

Molecular dynamics simulations have emerged as a powerful set of tools to unravel the intricate dynamics of ribosomes during protein synthesis. Recent advancements in this field have enabled simulations to delve deep into the conformational rearrangements of ribosomes and associated factors, providing invaluable insights into the intricacies of translation. Emphasis on simulations has recently been on translation elongation, such as tRNA selection, translocation, and ribosomal head-swivel motions. These studies have offered crucial structural interpretations of how genetic information is faithfully translated into proteins. This review outlines recent discoveries concerning ribosome conformational changes occurring during translation elongation, as elucidated through molecular dynamics simulations.

Geometric differences in the ribosome exit tunnel impact the escape of small nascent proteins

The exit tunnel is the subcompartment of the ribosome that contains the nascent polypeptide chain and, as such, is involved in various vital functions, including regulation of translation and protein folding. As the geometry of the tunnel shows important differences across species, we focus on key geometrical features of eukaryote and prokaryote tunnels. We used a simple coarse-grained molecular dynamics model to study the role of the tunnel geometry in the post-translational escape of short proteins (short open reading frames [sORFs]) with lengths ranging from 6 to 56 amino acids. We found that the probability of escape for prokaryotes is one for all but the 12-mer chains. Moreover, proteins of this length have an extremely low escape probability in eukaryotes. A detailed examination of the associated single trajectories and energy profiles showed that these variations can be explained by the interplay between the protein configurational space and the confinement effects introduced by the constriction sites of the ribosome exit tunnel. For certain lengths, either one or both of the constriction sites can lead to the trapping of the protein in the "pocket" regions preceding these sites. As the distribution of existing sORFs indicates some bias in length that is consistent with our findings, we finally suggest that the constraints imposed by the tunnel geometry have impacted the evolution of sORFs.

Binding of the peptide deformylase on the ribosome surface modulates the exit tunnel interior

Proteosynthesis on ribosomes is regulated at many levels. Conformational changes of the ribosome, possibly induced by external factors, may transfer over large distances and contribute to the regulation. The molecular principles of this long-distance allostery within the ribosome remain poorly understood. Here, we use structural analysis and atomistic molecular dynamics simulations to investigate peptide deformylase (PDF), an enzyme that binds to the ribosome surface near the ribosomal protein uL22 during translation and chemically modifies the emerging nascent peptide. Our simulations of the entire ribosome-PDF complex reveal that the PDF undergoes a swaying motion on the ribosome surface at the submicrosecond timescale. We show that the PDF affects the conformational dynamics of parts of the ribosome over distances of more than 5 nm. Using a supervised-learning algorithm, we demonstrate that the exit tunnel is influenced by the presence or absence of PDF. Our findings suggest a possible effect of the PDF on the nascent peptide translocation through the ribosome exit tunnel.

Reversible Unwrapping Algorithm for Constant-Pressure Molecular Dynamics Simulations

Molecular simulation technologies have afforded researchers a unique look into the nanoscale interactions driving physical processes. However, a limitation for molecular dynamics (MD) simulations is that they must be performed on finite-sized systems in order to map onto computational resources. To minimize artifacts arising from finite-sized simulation systems, it is common practice for MD simulations to be performed with periodic boundary conditions (PBCs). However, in order to calculate specific physical properties, such as mean square displacements to calculate diffusion coefficients, continuous particle trajectories where the atomic movements are continuous and do not jump between cell faces are required. In these cases, modifying atomic coordinates through unwrapping schemes is an essential post-processing tool to remove these jumps. Here, two established trajectory unwrapping schemes are applied to 1 μs wrapped trajectories for a small water box and lysozyme in water. The existing schemes can result in spurious diffusion coefficients, long bonds within unwrapped molecules, and inconsistent atomic coordinates when coordinates are rewrapped after unwrapping. We determine that prior unwrapping schemes do not account for changing periodic box dimensions and introduce an additional correction term to the existing displacement unwrapping scheme to correct for these artifacts. We also demonstrate that the resulting algorithm is a hybrid between the existing heuristic and displacement unwrapping schemes. After treatment using this new unwrapping scheme, molecular geometries are correct even after long simulations. In anticipation for longer MD trajectories, we develop implementations for this new scheme in multiple PBC handling tools.

Hyper-swivel head domain motions are required for complete mRNA-tRNA translocation and ribosome resetting

Translocation of messenger RNA (mRNA) and transfer RNA (tRNA) substrates through the ribosome during protein synthesis, an exemplar of directional molecular movement in biology, entails a complex interplay of conformational, compositional, and chemical changes. The molecular determinants of early translocation steps have been investigated rigorously. However, the elements enabling the ribosome to complete translocation and reset for subsequent protein synthesis reactions remain poorly understood. Here, we have combined molecular simulations with single-molecule fluorescence resonance energy transfer imaging to gain insights into the rate-limiting events of the translocation mechanism. We find that diffusive motions of the ribosomal small subunit head domain to hyper-swivelled positions, governed by universally conserved rRNA, can maneuver the mRNA and tRNAs to their fully translocated positions. Subsequent engagement of peptidyl-tRNA and disengagement of deacyl-tRNA from mRNA, within their respective small subunit binding sites, facilitate the ribosome resetting mechanism after translocation has occurred to enable protein synthesis to resume.

Probing Interplays between Human XBP1u Translational Arrest Peptide and 80S Ribosome

The ribosome stalling mechanism is a crucial biological process, yet its atomistic underpinning is still elusive. In this framework, the human XBP1u translational arrest peptide (AP) plays a central role in regulating the unfolded protein response (UPR) in eukaryotic cells. Here, we report multimicrosecond all-atom molecular dynamics simulations designed to probe the interactions between the XBP1u AP and the mammalian ribosome exit tunnel, both for the wild type AP and for four mutant variants of different arrest potencies. Enhanced sampling simulations allow investigating the AP release process of the different variants, shedding light on this complex mechanism. The present outcomes are in qualitative/quantitative agreement with available experimental data. In conclusion, we provide an unprecedented atomistic picture of this biological process and clear-cut insights into the key AP-ribosome interactions.

The ugly, bad, and good stories of large-scale biomolecular simulations

Molecular modeling of large biomolecular assemblies exemplifies a disruptive area holding both promises and contentions. Propelled by peta and exascale computing, several simulation methodologies have now matured into user-friendly tools that are successfully employed for modeling viruses, membranous nano-constructs, and key pieces of the genetic machinery. We present three unifying biophysical themes that emanate from some of the most recent multi-million atom simulation endeavors. Despite connecting molecular changes with phenotypic outcomes, the quality measures of these simulations remain questionable. We discuss the existing and upcoming strategies for constructing representative ensembles of large systems, how new computing technologies will boost this area, and make a point that integrative modeling guided by experimental data is the future of biomolecular computations.

Folding of VemP into translation-arresting secondary structure is driven by the ribosome exit tunnel

The ribosome is a fundamental biomolecular complex that synthesizes proteins in cells. Nascent proteins emerge from the ribosome through a tunnel, where they may interact with the tunnel walls or small molecules such as antibiotics. These interactions can cause translational arrest with notable physiological consequences. Here, we studied the arrest caused by the regulatory peptide VemP, which is known to form α-helices inside the ribosome tunnel near the peptidyl transferase center under specific conditions. We used all-atom molecular dynamics simulations of the entire ribosome and circular dichroism spectroscopy to study the driving forces of helix formation and how VemP causes the translational arrest. To that aim, we compared VemP dynamics in the ribosome tunnel with its dynamics in solution. We show that the VemP peptide has a low helical propensity in water and that the propensity is higher in mixtures of water and trifluorethanol. We propose that helix formation within the ribosome is driven by the interactions of VemP with the tunnel and that a part of VemP acts as an anchor. This anchor might slow down VemP progression through the tunnel enabling α-helix formation, which causes the elongation arrest.

Searching for Low Probability Opening Events in a DNA Sliding Clamp

The β subunit of E. coli DNA polymererase III is a DNA sliding clamp associated with increasing the processivity of DNA synthesis. In its free form, it is a circular homodimer structure that can accomodate double-stranded DNA in a nonspecific manner. An open state of the clamp must be accessible before loading the DNA. The opening mechanism is still a matter of debate, as is the effect of bound DNA on opening/closing kinetics. We use a combination of atomistic, coarse-grained, and enhanced sampling strategies in both explicit and implicit solvents to identify opening events in the sliding clamp. Such simulations of large nucleic acid and their complexes are becoming available and are being driven by improvements in force fields and the creation of faster computers. Different models support alternative opening mechanisms, either through an in-plane or out-of-plane opening event. We further note some of the current limitations, despite advances, in modeling these highly charged systems with implicit solvent.

Comparative Analysis of Structural and Dynamical Features of Ribosome Upon Association With mRNA Reveals Potential Role of Ribosomal Proteins

Ribosomes play a critical role in maintaining cellular proteostasis. The binding of messenger RNA (mRNA) to the ribosome regulates kinetics of protein synthesis. To generate an understanding of the structural, mechanistic, and dynamical features of mRNA recognition in the ribosome, we have analysed mRNA-protein interactions through a structural comparison of the ribosomal complex in the presence and absence of mRNA. To do so, we compared the 3-Dimensional (3D) structures of components of the two assembly structures and analysed their structural differences because of mRNA binding, using elastic network models and structural network-based analysis. We observe that the head region of 30S ribosomal subunit undergoes structural displacement and subunit rearrangement to accommodate incoming mRNA. We find that these changes are observed in proteins that lie far from the mRNA-protein interface, implying allostery. Further, through perturbation response scanning, we show that the proteins S13, S19, and S20 act as universal sensors that are sensitive to changes in the inter protein network, upon binding of 30S complex with mRNA and other initiation factors. Our study highlights the significance of mRNA binding in the ribosome complex and identifies putative allosteric sites corresponding to alterations in structure and/or dynamics, in regions away from mRNA binding sites in the complex. Overall, our work provides fresh insights into mRNA association with the ribosome, highlighting changes in the interactions and dynamics of the ribosome assembly because of the binding.

Mechanistic studies of non-canonical amino acid mutagenesis

Over the past decade, harnessing the cellular protein synthesis machinery to incorporate non-canonical amino acids (ncAAs) into tailor-made peptides has significantly advanced many aspects of molecular science. More recently, groundbreaking progress in our ability to engineer this machinery for improved ncAA incorporation has led to significant enhancements of this powerful tool for biology and chemistry. By revealing the molecular basis for the poor or improved incorporation of ncAAs, mechanistic studies of ncAA incorporation by the protein synthesis machinery have tremendous potential for informing and directing such engineering efforts. In this chapter, we describe a set of complementary biochemical and single-molecule fluorescence assays that we have adapted for mechanistic studies of ncAA incorporation. Collectively, these assays provide data that can guide engineering of the protein synthesis machinery to expand the range of ncAAs that can be incorporated into peptides and increase the efficiency with which they can be incorporated, thereby enabling the full potential of ncAA mutagenesis technology to be realized.

Nascent Folding of Proteins Across the Three Domains of Life

We study the nascent behavior of three model coarse-grained proteins in six rigid all-atom structures representing ribosomes that come from three domains of life. The synthesis of the proteins is implemented as a growth process. The geometry of the exit tunnel is quantified and shown to differ between the domains of life: both in volume and the size of constriction sites. This results in different characteristic times of capture within the tunnel and various probabilities of the escape. One of the proteins studied is the bacterial YibK which is knotted in its native state. A fraction of the trajectories results in knotting and the probability of doing so is largest for the bacterial ribosomes. Relaxing the condition of the rigidness of the ribosomes should result in a better avoidance of trapping and better proper folding.

How Does the Ribosome Fold the Proteome?

Folding of polypeptides begins during their synthesis on ribosomes. This process has evolved as a means for the cell to maintain proteostasis, by mitigating the risk of protein misfolding and aggregation. The capacity to now depict this cellular feat at increasingly higher resolution is providing insight into the mechanistic determinants that promote successful folding. Emerging from these studies is the intimate interplay between protein translation and folding, and within this the ribosome particle is the key player. Its unique structural properties provide a specialized scaffold against which nascent polypeptides can begin to form structure in a highly coordinated, co-translational manner. Here, we examine how, as a macromolecular machine, the ribosome modulates the intrinsic dynamic properties of emerging nascent polypeptide chains and guides them toward their biologically active structures.

Pathway regulation mechanism revealed by cotranslational folding of villin headpiece subdomain HP35

Cotranslational folding is one of the most important features of protein folding in vivo. Although many studies have shown that the folding pathways of cotranslational folding are different from free folding in vitro, the detailed mechanism of folding dynamics is lacking. Here we combine all-atom molecular simulations with an ideal ribosome tunnel model to investigate the cotranslational folding of villin headpiece subdomain HP35. By comparing the folding dynamics between cotranslational folding and free folding, we found that cotranslational folding tends to fold along the pathway that is easier to fold into native state in the latter. In addition, the roles of the ribosome tunnel and sequential folding are analyzed separately. Our results show that the ribosome can prevent the untimely folding of the C segment of HP35 to reduce the non-native interactions, while the translation speed can regulate the amounts of native and non-native interactions and the balance between them. Overall, these results give insights into the general mechanisms of cotranslational protein folding.

Piecewise All-Atom SMD Simulations Reveal Key Secondary Structures in Luciferase Unfolding Pathway

Although the folding of single-domain proteins is well characterized theoretically and experimentally, the folding of large multidomain proteins is less well known. Firefly luciferase, a 550 residue three-domain protein, has been commonly used as a substrate to study chaperone reactions and as a model system for the study of folding of long polypeptide chains, including related phenomena such as cotranslational folding. Despite being characterized by various experimental techniques, the atomic-level contributions of various secondary structures of luciferase to its fold's mechanical stability remain unknown. Here, we developed a piecewise approach for all-atom steered molecular dynamics simulations to examine specific secondary structures that resist mechanical unfolding while minimizing the amount of computational resources required by the large water box of standard all-atom steered molecular dynamics simulations. We validated the robustness of this approach with a small NI3C protein and used our approach to elucidate the specific secondary structures that provide the largest contributions to luciferase mechanostability. In doing so, we show that piecewise all-atom steered molecular dynamics simulations can provide novel atomic resolution details regarding mechanostability and can serve as a platform for novel mutagenesis studies as well as a point for comparison with high-resolution force spectroscopy experiments.

Advanced Methods for Studying Structure and Interactions of Macrolide Antibiotics

Macrolide antibiotics are macrocyclic compounds that are clinically used and prescribed for the treatment of upper and lower respiratory tract infections. They inhibit the synthesis of bacterial proteins by reversible binding to the 23S rRNA at or near the peptidyl transferase center. However, their excellent antibacterial profile was largely compromised by the emergence of bacterial resistance. Today, fighting resistance to antibiotics is one of the greatest challenges in medicinal chemistry. Considering various physicochemical properties of macrolides, understanding their structure and interactions with macromolecular targets is crucial for the design of new antibiotics efficient against resistant pathogens. The solid-state structures of some macrolide-ribosome complexes have recently been solved, throwing new light on the macrolide binding mechanisms. On the other hand, a combination of NMR spectroscopy and molecular modeling calculations can be applied to study free and bound conformations in solution. In this article, a description of advanced physicochemical methods for elucidating the structure and interactions of macrolide antibiotics in solid state and solution will be provided, and their principal advantages and drawbacks will be discussed.

Structural Heterogeneities of the Ribosome: New Frontiers and Opportunities for Cryo-EM

The extent of ribosomal heterogeneity has caught increasing interest over the past few years, as recent studies have highlighted the presence of structural variations of the ribosome. More precisely, the heterogeneity of the ribosome covers multiple scales, including the dynamical aspects of ribosomal motion at the single particle level, specialization at the cellular and subcellular scale, or evolutionary differences across species. Upon solving the ribosome atomic structure at medium to high resolution, cryogenic electron microscopy (cryo-EM) has enabled investigating all these forms of heterogeneity. In this review, we present some recent advances in quantifying ribosome heterogeneity, with a focus on the conformational and evolutionary variations of the ribosome and their functional implications. These efforts highlight the need for new computational methods and comparative tools, to comprehensively model the continuous conformational transition pathways of the ribosome, as well as its evolution. While developing these methods presents some important challenges, it also provides an opportunity to extend our interpretation and usage of cryo-EM data, which would more generally benefit the study of molecular dynamics and evolution of proteins and other complexes.

Recent Advances in Understanding Biological GTP Hydrolysis through Molecular Simulation

GTP hydrolysis is central to biology, being involved in regulating a wide range of cellular processes. However, the mechanisms by which GTPases hydrolyze this critical reaction remain controversial, with multiple mechanistic possibilities having been proposed based on analysis of experimental and computational data. In this mini-review, we discuss advances in our understanding of biological GTP hydrolysis based on recent computational studies and argue in favor of solvent-assisted hydrolysis as a conserved mechanism among GTPases. A concrete understanding of the fundamental mechanisms by which these enzymes facilitate GTP hydrolysis will have significant impact both for drug discovery efforts and for unraveling the role of oncogenic mutations.

tRNA Dissociation from EF-Tu after GTP Hydrolysis: Primary Steps and Antibiotic Inhibition

In each round of ribosomal translation, the translational GTPase elongation factor Tu (EF-Tu) delivers a transfer RNA (tRNA) to the ribosome. After successful decoding, EF-Tu hydrolyzes GTP, which triggers a conformational change that ultimately results in the release of the tRNA from EF-Tu. To identify the primary steps of these conformational changes and how they are prevented by the antibiotic kirromycin, we employed all-atom explicit-solvent molecular dynamics simulations of the full ribosome-EF-Tu complex. Our results suggest that after GTP hydrolysis and Pi release, the loss of interactions between the nucleotide and the switch 1 loop of EF-Tu allows domain D1 of EF-Tu to rotate relative to domains D2 and D3 and leads to an increased flexibility of the switch 1 loop. This rotation induces a closing of the D1-D3 interface and an opening of the D1-D2 interface. We propose that the opening of the D1-D2 interface, which binds the CCA tail of the tRNA, weakens the crucial EF-Tu-tRNA interactions, which lowers tRNA binding affinity, representing the first step of tRNA release. Kirromycin binds within the D1-D3 interface, sterically blocking its closure, but does not prevent hydrolysis. The resulting increased flexibility of switch 1 explains why it is not resolved in kirromycin-bound structures.

Assembly of heteropolymers via a network of reaction coordinates

In biochemistry, heteropolymers encoding biological information are assembled out of equilibrium by sequentially incorporating available monomers found in the environment. Current models of polymerization treat monomer incorporation as a sequence of discrete chemical reactions between intermediate metastable states. In this paper, we use ideas from reaction rate theory and describe nonequilibrium assembly of a heteropolymer via a continuous reaction coordinate. Our approach allows for estimating the copy error and incorporation speed from the Gibbs free energy landscape of the process. We apply our theory to several examples from a simple reaction characterized by a free energy barrier to more complex cases incorporating error correction mechanisms, such as kinetic proofreading.

Principles of tRNAAla Selection by Alanyl-tRNA Synthetase Based on the Critical G3·U70 Base Pair

Throughout evolution, the presence of a single G3·U70 mismatch in the acceptor stem of tRNA^Ala is the major determinant for aminoacylation with alanine by alanyl-tRNA synthetase (AlaRS). Recently reported crystal structures of the complexes AlaRS-tRNA^Ala/G3·U70 and AlaRS-tRNA^Ala/A3·U70 suggest two very different conformations, representing a reactive and a nonreactive state, respectively. On the basis of these structures, it has been proposed that the G3·U70 base pair guides the -CCA end of the tRNA acceptor stem into the active site of AlaRS, thereby enabling aminoacylation. The crystal structures open up the possibility of directly computing the energetics of tRNA specificity by AlaRS. We have carried out molecular dynamics free-energy simulations to quantitatively estimate tRNA discrimination by AlaRS, focusing on the mutations of the single critical base pair G3·U70 to uncover the energetics underlying the accuracy of tRNA selection. The calculations show that the reactive complex is highly selective in favor of the cognate tRNA^Ala/G3·U70 over its noncognate analogues (A3·U70/G3·C70/A3·C70). In contrast, the nonreactive complex is predicted to be unselective between tRNA^Ala/G3·U70 and tRNA^Ala/A3·U70. Utilizing our calculated relative binding free energies, we show how a simple three-step kinetic scheme for aminoacylation, involving both an initial nonspecific binding step and a subsequent transition to a selective reactive complex, accounts for the observed kinetics of the process.

Role of Computational Methods in Going beyond X-ray Crystallography to Explore Protein Structure and Dynamics

Protein structural biology came a long way since the determination of the first three-dimensional structure of myoglobin about six decades ago. Across this period, X-ray crystallography was the most important experimental method for gaining atomic-resolution insight into protein structures. However, as the role of dynamics gained importance in the function of proteins, the limitations of X-ray crystallography in not being able to capture dynamics came to the forefront. Computational methods proved to be immensely successful in understanding protein dynamics in solution, and they continue to improve in terms of both the scale and the types of systems that can be studied. In this review, we briefly discuss the limitations of X-ray crystallography in studying protein dynamics, and then provide an overview of different computational methods that are instrumental in understanding the dynamics of proteins and biomacromolecular complexes.