High performance computing in biology: multimillion atom simulations of nanoscale systems

Computational methods have been used in biology for sequence analysis (bioinformatics), all-atom simulation (molecular dynamics and quantum calculations), and more recently for modeling biological networks (systems biology). Of these three techniques, all-atom simulation is currently the most computationally demanding, in terms of compute load, communication speed, and memory load. Breakthroughs in electrostatic force calculation and dynamic load balancing have enabled molecular dynamics simulations of large biomolecular complexes. Here, we report simulation results for the ribosome, using approximately 2.64 million atoms, the largest all-atom biomolecular simulation published to date. Several other nano-scale systems with different numbers of atoms were studied to measure the performance of the NAMD molecular dynamics simulation program on the Los Alamos National Laboratory Q Machine. We demonstrate that multimillion atom systems represent a 'sweet spot' for the NAMD code on large supercomputers. NAMD displays an unprecedented 85% parallel scaling efficiency for the ribosome system on 1024 CPUs. We also review recent targeted molecular dynamics simulations of the ribosome that prove useful for studying conformational changes of this large biomolecular complex in atomic detail.     

单分子荧光共振能量转移技术在核糖体移位研究中的进展

核糖体是蛋白质的"合成工厂",也是临床上多种抗菌药物的作用靶点,因此,深入理解细菌核糖体的蛋白质翻译机制意义重大.蛋白质翻译是通过多步骤相互协调、多组分精细配合来实现高保真和精确调控.核糖体在mRNA上的移位作为翻译过程中最重要的事件之一,需要核糖体大规模的构象重排以及tRNA2-mRNA沿着核糖体的精确移动.在细菌中,移位是由延伸因子EF-G催化GTP水解来驱动的.近年来,单分子荧光共振能量技术(smFRET)的发展使得人们可以探究单个tRNA分子移位的动力学过程并实时观测核糖体的构象变化.本文首先介绍了smFRET技术的原理及特点,对其在核糖体结构动态及tRNA移位研究中的应用进行了较为系统的总结,并对其应用前景进行了展望.     

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.     

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.     

The Role of L1 Stalk-tRNA Interaction in the Ribosome Elongation Cycle

The ribosomal L1 stalk is a mobile structure implicated in directing tRNA movement during translocation through the ribosome. This article investigates three aspects of L1 stalk-tRNA interaction. First, by combining data from cryo electron microscopy, X-ray crystallography, and molecular dynamics simulations through the molecular dynamics flexible fitting method, we obtained atomic models of different tRNAs occupying the hybrid P/E state interacting with the L1 stalk. These models confirm the assignment of fluorescence resonance energy transfer states from previous single-molecule investigations of L1 stalk dynamics. Second, the models reconcile how initiator tRNA^fMet interacts less strongly with the L1 stalk compared to elongator tRNA^Phe, as seen in previous single-molecule experiments. Third, results from a simulation of the entire ribosome in which the L1 stalk is moved from a half-closed conformation to its open conformation are found to support the hypothesis that L1 stalk opening is involved in tRNA release from the ribosome.     

Structural dynamics of the ribosome

Protein synthesis is inherently a dynamic process, requiring both small-scale and large-scale movements of tRNA and mRNA. It has long been suspected that these movements might be coupled to conformational changes in the ribosome, and in its RNA moieties in particular. Recently, the nature of ribosome structural dynamics has begun to emerge from a combination of approaches, most notably cryo-EM, X-ray crystallography, and FRET. Ribosome movement occurs both on a grand scale, as in the intersubunit rotational movements that are coupled to tRNA-mRNA translocation, and in intricate localized rearrangements such as those that accompany codon-anticodon recognition and peptide bond formation. In spite of much progress, our understanding of the mechanics of translation is now beset with countless new questions, reflecting the vast molecular architecture of the ribosome itself.     

Direct Measurements of Erythromycin’s Effect on Protein Synthesis Kinetics in Living Bacterial Cells

Macrolide antibiotics, such as erythromycin, bind to the nascent peptide exit tunnel (NPET) of the bacterial ribosome and modulate protein synthesis depending on the nascent peptide sequence. Whereas in vitro biochemical and structural methods have been instrumental in dissecting and explaining the molecular details of macrolide-induced peptidyl-tRNA drop-off and ribosome stalling, the dynamic effects of the drugs on ongoing protein synthesis inside live bacterial cells are far less explored. In the present study, we used single-particle tracking of dye-labeled tRNAs to study the kinetics of mRNA translation in the presence of erythromycin, directly inside live Escherichia coli cells. In erythromycin-treated cells, we find that the dwells of elongator tRNA^Phe on ribosomes extend significantly, but they occur much more seldom. In contrast, the drug barely affects the ribosome binding events of the initiator tRNA^fMet. By overexpressing specific short peptides, we further find context-specific ribosome binding dynamics of tRNA^Phe, underscoring the complexity of erythromycin’s effect on protein synthesis in bacterial cells.     

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.     

Dynamics of the translational machinery

The recent growth in single molecule studies of translation has provided an insight into the molecular mechanism of ribosomal function. Single molecule fluorescence approaches allowed direct observation of the structural rearrangements occurring during translation and revealed dynamic motions of the ribosome and its ligands. These studies demonstrated how ligand binding affects dynamics of the ribosome, and the role of the conformational sampling in large-scale rearrangements intrinsic to translation elongation. The application of time-resolved cryo-electron microscopy revealed new conformational intermediates during back-translocation providing an insight into ribosomal dynamics from an alternative perspective. Recent developments permitted examination of conformational and compositional dynamics of the ribosome in real-time through multiple cycles of elongation at the single molecule level. The zero-mode waveguide approach allowed direct observation of the compositional dynamics of tRNA occupancy on the elongating ribosome. The emergence of single molecule in vivo techniques provided insights into the mechanism and regulation of translation at the organismal level.     

Energetics of codon-anticodon recognition on the small ribosomal subunit

Recent crystal structures of the small ribosomal subunit have made it possible to examine the detailed energetics of codon recognition on the ribosome by computational methods. The binding of cognate and near-cognate anticodon stem loops to the ribosome decoding center, with mRNA containing the Phe UUU and UUC codons, are analyzed here using explicit solvent molecular dynamics simulations together with the linear interaction energy (LIE) method. The calculated binding free energies are in excellent agreement with experimental binding constants and reproduce the relative effects of mismatches in the first and second codon position versus a mismatch at the wobble position. The simulations further predict that the Leu2 anticodon stem loop is about 10 times more stable than the Ser stem loop in complex with the Phe UUU codon. It is also found that the ribosome significantly enhances the intrinsic stability differences of codon-anticodon complexes in aqueous solution. Structural analysis of the simulations confirms the previously suggested importance of the universally conserved nucleotides A1492, A1493, and G530 in the decoding process.     

Intrinsic molecular properties of the protein-protein bridge facilitate ratchet-like motion of the ribosome

The ribosomal intersubunit bridges maintain the overall architecture of the ribosome and thereby play a pivotal role in the dynamics of translation. The only protein-protein bridge, b1b, is formed by the two proteins, S13 and L5 of the small and large ribosomal subunits, respectively. B1b absorbs the largest movement during ratchet-like motion, and its two proteins reorganize in different constellations during this motion of the ribosome.Our results in this study of b1b in the Escherichia coli 70S ribosome suggest that the intrinsic molecular features of the bridging proteins allow the bridge to modulate the ratchet-like motion in a controlled manner. Additionally, another large subunit protein, L31, seems to participate with S13 and L5 in the formation, dynamics, and stabilization of this bridge.     

A noncanonical binding site of chloramphenicol revealed via molecular dynamics simulations

Chloramphenicol, an antibiotic belonging to the family of amphenicols, is an inhibitor of translation. On the basis of X–ray structural analysis of the binding of chloramphenicol to free bacterial ribosomes, the chloramphenicol action mechanism that consists in preventing the binding of aminoacyl-tRNA to the A–site of the large subunit of the ribosome was adopted. However, the known structures of chloramphenicol complexes with bacterial ribosomes poorly explain the results of the experiments on the chemical modification of 23S rRNA, the resistance to chloramphenicol caused by mutations in 23S rRNA and, which is particularly important, the selectivity of chloramphenicol in suppression of translation, depending on the amino acid sequence of the nascent peptide. In the present study the putative structure of the chloramphenicol complex with a bacterial ribosome in the A,A/P,P–state has been obtained by molecular dynamics simulations methods. The proposed structure of the complex allows us to explain the results of biochemical studies of the interaction of chloramphenicol with the bacterial ribosome.     

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.     

Flexible fitting of atomic structures into electron microscopy maps using molecular dynamics

A novel method to flexibly fit atomic structures into electron microscopy (EM) maps using molecular dynamics simulations is presented. The simulations incorporate the EM data as an external potential added to the molecular dynamics force field, allowing all internal features present in the EM map to be used in the fitting process, while the model remains fully flexible and stereochemically correct. The molecular dynamics flexible fitting (MDFF) method is validated for available crystal structures of protein and RNA in different conformations; measures to assess and monitor the fitting process are introduced. The MDFF method is then used to obtain high-resolution structures of the E. coli ribosome in different functional states imaged by cryo-EM.     

Non-bulk-like solvent behavior in the ribosome exit tunnel

As nascent proteins are synthesized by the ribosome, they depart via an exit tunnel running through the center of the large subunit. The exit tunnel likely plays an important part in various aspects of translation. Although water plays a key role in many bio-molecular processes, the nature of water confined to the exit tunnel has remained unknown. Furthermore, solvent in biological cavities has traditionally been characterized as either a continuous dielectric fluid, or a discrete tightly bound molecule. Using atomistic molecular dynamics simulations, we predict that the thermodynamic and kinetic properties of water confined within the ribosome exit tunnel are quite different from this simple two-state model. We find that the tunnel creates a complex microenvironment for the solvent resulting in perturbed rotational dynamics and heterogenous dielectric behavior. This gives rise to a very rugged solvation landscape and significantly retarded solvent diffusion. We discuss how this non-bulk-like solvent is likely to affect important biophysical processes such as sequence dependent stalling, co-translational folding, and antibiotic binding. We conclude with a discussion of the general applicability of these results to other biological cavities.     

Refinement of a low-resolution crystal structure to better understand erythromycin interactions on large ribosomal subunit

Macrolides are a group of diverse class of naturally occurring and synthetic antibiotics made of macrocyclic-lactone ring carrying one or more sugar moieties linked to various atoms of the lactone ring. These macrolides selectively bind to a single high affinity site on the prokaryotic 50S ribosomal subunit, making them highly effective towards a wide range of bacterial pathogens. The understanding of binding between macrolides and ribosome serves a good basis in elucidating how they work at the molecular level and these findings would be important in rational drug design. Here, we report refinement of reconstructed PDB structure of erythromycin-ribosome system using molecular dynamics (MD) simulation. Interesting findings were observed in this refinement stage that could improve the understanding of the binding of erythromycin A (ERYA) onto the 50S subunit. The results showed ERYA was highly hydrated and water molecules were found to be important in bridging hydrogen bond at the binding pocket during the simulation time. ERYA binding to ribosome was also strengthened by hydrogen bond network and hydrophobic interactions between the antibiotic and the ribosome. Our MD simulation also demonstrated direct interaction of ERYA with Domains II, V and with C1773 (U1782EC), a residue in Domain IV that has yet been described of its role in ERYA binding. It is hoped that this refinement will serve as a starting model for a further enhancement of our understanding towards the binding of ERYA to ribosome.     

Single‐molecule dynamics of the molecular chaperone trigger factor in living cells

In bacteria, trigger factor (TF) is the molecular chaperone that interacts with the ribosome to assist the folding of nascent polypeptides. Studies in vitro have provided insights into the function and mechanism of TF. Much is to be elucidated, however, about how TF functions in vivo. Here, we use single‐molecule tracking, in combination with genetic manipulations, to study the dynamics and function of TF in living E. coli cells. We find that TF, besides interacting with the 70S ribosome, may also bind to ribosomal subunits and form TF‐polypeptide complexes that may include DnaK/DnaJ proteins. The TF‐70S ribosome interactions are highly dynamic inside cells, with an average residence time of ～0.2 s. Our results confirm that the signal recognition particle weakens TF's interaction with the 70S ribosome, and further identify that this weakening mainly results from a change in TF's binding to the 70S ribosome, rather than its unbinding. Moreover, using photoconvertible bimolecular fluorescence complementation, we selectively probe TF₂ dimers in the cell and show that TF₂ does not bind to the 70S ribosome but is involved in the post‐translational interactions with polypeptides. These findings contribute to the fundamental understanding of molecular chaperones in assisting protein folding in living cells.     

Single-molecule study of viomycin's inhibition mechanism on ribosome translocation

Viomycin belongs to the tuberactinomycin family of antibiotics against tuberculosis. However, its inhibition mechanism remains elusive. Although it is clear that viomycin inhibits the ribosome intersubunit ratcheting, there are contradictory reports about whether the antibiotic viomycin stabilizes the tRNA hybrid or classical state. By using a single-molecule FRET method to directly observe the tRNA dynamics relative to ribosomal protein L27, we have found that viomycin trapped the hybrid state within certain ribosome subgroups but did not significantly suppress the tRNA dynamics. The persistent fluctuation of tRNA implied that tRNA motions were decoupled from the ribosome intersubunit ratcheting. Viomycin also promoted peptidyl-tRNA fluctuation in the posttranslocation complex, implying that, in addition to acylated P-site tRNA, the decoding center also played an important role of ribosome locking after translocation. Therefore, viomycin inhibits translocation by trapping the hybrid state in the pretranslocation complex and disturbing the stability of posttranslocation complex. Our results imply that ribosome translocation is possibly a synergistic process of multiple decoupled local dynamics.     

The ribosome as an allosterically regulated molecular machine

The ribosome as a complex molecular machine undergoes significant conformational rearrangements during the synthesis of polypeptide chains of proteins. In this review, information obtained using various experimental methods on the internal consistency of such rearrangements is discussed. It is demonstrated that allosteric regulation involves all the main stages of the operation of the ribosome and connects functional elements remote by tens and even hundreds of angstroms. Data obtained using Förster resonance energy transfer (FRET) show that translocation is controlled in general by internal mechanisms of the ribosome, and not by the position of the ligands. Chemical probing data revealed the relationship of such remote sites as the decoding, peptidyl transferase, and GTPase centers of the ribosome. Nevertheless, despite the large amount of experimental data accumulated to date, many details and mechanisms of these phenomena are still not under- stood. Analysis of these data demonstrates that the development of new approaches is necessary for deciphering the mech- anisms of allosteric regulation of the operation of the ribosome.     

Applications of the molecular dynamics flexible fitting method

In recent years, cryo-electron microscopy (cryo-EM) has established itself as a key method in structural biology, permitting the structural characterization of large biomolecular complexes in various functional states. The data obtained through single-particle cryo-EM has recently seen a leap in resolution thanks to landmark advances in experimental and computational techniques, resulting in sub-nanometer resolution structures being obtained routinely. The remaining gap between these data and revealing the mechanisms of molecular function can be closed through hybrid modeling tools that incorporate known atomic structures into the cryo-EM data. One such tool, molecular dynamics flexible fitting (MDFF), uses molecular dynamics simulations to combine structures from X-ray crystallography with cryo-EM density maps to derive atomic models of large biomolecular complexes. The structures furnished by MDFF can be used subsequently in computational investigations aimed at revealing the dynamics of the complexes under study. In the present work, recent applications of MDFF are presented, including the interpretation of cryo-EM data of the ribosome at different stages of translation and the structure of a membrane-curvature-inducing photosynthetic complex.