Search

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein binds the receptor angiotensin converting enzyme 2 (ACE2) and drives virus-host membrane fusion through refolding of its S2 domain. Whereas the S1 domain contains high sequence variability, the S2 domain is conserved and is a promising pan-betacoronavirus vaccine target. We applied cryo-electron tomography to capture intermediates of S2 refolding and understand inhibition by antibodies to the S2 stem-helix. Subtomogram averaging revealed ACE2 dimers cross-linking spikes before transitioning into S2 intermediates, which were captured at various stages of refolding. Pan-betacoronavirus neutralizing antibodies targeting the S2 stem-helix bound to and inhibited refolding of spike prehairpin intermediates. Combined with molecular dynamics simulations, these structures elucidate the process of SARS-CoV-2 entry and reveal how pan-betacoronavirus S2-targeting antibodies neutralize infectivity by arresting prehairpin intermediates.

The mitoribosome translates mitochondrial mRNAs and regulates energy conversion that is a signature of aerobic life forms. We present a 2.2 Å resolution structure of human mitoribosome together with validated mitoribosomal RNA (rRNA) modifications, including aminoacylated CP-tRNA Val . The structure shows how mitoribosomal proteins stabilise binding of mRNA and tRNA helping to align it in the decoding center, whereas the GDP-bound mS29 stabilizes intersubunit communication. Comparison between different states, with respect to tRNA position, allowed us to characterize a non-canonical L1 stalk, and molecular dynamics simulations revealed how it facilitates tRNA transitions in a way that does not require interactions with rRNA. We also report functionally important polyamines that are depleted when cells are subjected to an antibiotic treatment. The structural, biochemical, and computational data illuminate the principal functional components of the translation mechanism in mitochondria and provide a description of the structure and function of the human mitoribosome., (© 2024. The Author(s).)

The ribosome is a prototypical assembly that can be used to establish general principles and techniques for the study of biological molecular machines. Motivated by the fact that the dynamics of every biomolecule is governed by an underlying energy landscape, there has been great interest to understand and quantify ribosome energetics. In the present review, we will focus on theoretical and computational strategies for probing the interactions that shape the energy landscape of the ribosome, with an emphasis on more recent studies of the elongation cycle. These efforts include the application of quantum mechanical methods for describing chemical kinetics, as well as classical descriptions to characterize slower (microsecond to millisecond) large-scale (10-100 Å) rearrangements, where motion is described in terms of diffusion across an energy landscape. Together, these studies provide broad insights into the factors that control a diverse range of dynamical processes in this assembly., (© 2023 The Authors. Biopolymers published by Wiley Periodicals LLC.)

Protein synthesis by the ribosome requires large-scale rearrangements of the 'small' subunit (SSU; ∼1 MDa), including inter- and intra-subunit rotational motions. However, with nearly 2000 structures of ribosomes and ribosomal subunits now publicly available, it is exceedingly difficult to design experiments based on analysis of all known rotation states. To overcome this, we developed an approach where the orientation of each SSU head and body is described in terms of three angular coordinates (rotation, tilt and tilt direction) and a single translation. By considering the entire RCSB PDB database, we describe 1208 fully-assembled ribosome complexes and 334 isolated small subunits, which span >50 species. This reveals aspects of subunit rearrangements that are universal, and others that are organism/domain-specific. For example, we show that tilt-like rearrangements of the SSU body (i.e. 'rolling') are pervasive in both prokaryotic and eukaryotic (cytosolic and mitochondrial) ribosomes. As another example, domain orientations associated with frameshifting in bacteria are similar to those found in eukaryotic ribosomes. Together, this study establishes a common foundation with which structural, simulation, single-molecule and biochemical efforts can more precisely interrogate the dynamics of this prototypical molecular machine., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Conformational changes are an essential feature for the function of some dynamic proteins. Understanding the mechanism of such motions may allow us to identify important properties, which may be directly related to the regulatory function of a protein. Also, this knowledge may be employed for a rational design of drugs that can shift the balance between active and inactive conformations, as well as affect the kinetics of the activation process. Here, the conformational changes in carboxyl-terminal Src kinase, the major catalytic repressor to the Src family of kinases, was investigated, and it was proposed as a functionally related hypothesis. A Cα Structure-Based Model (Cα-SBM) was applied to provide a description of the overall conformational landscape and further analysis complemented by detailed molecular dynamics simulations. As a first approach to Cα-SBM simulations, reversible transitions between active (closed) and inactive (open) forms were modeled as fluctuations between these two energetic basins. It was found that, in addition to the interdomain Carboxyl-terminal SRC Kinase (Csk) correlated motions, a conformational change in the αC helix is required for a complete conformational transition. The result reveals this as an important region of transition control and domain coordination. Restrictions in the αC helix region of the Csk protein were performed, and the analyses showed a direct correlation with the global conformational changes, with this location being propitious for future studies of ligands. Also, the Src Homology 3 (SH3) and SH3 plus Src Homology 2 (SH2) domains were excluded for a direct comparison with experimental results previously published. Simulations where the SH3 was deleted presented a reduction of the transitions during the simulations, while the SH3-SH2 deletion vanishes the Csk transitions, corroborating the experimental results mentioned and linking the conformational changes with the catalytic functionality of Csk. The study was complemented by the introduction of a known kinase inhibitor close to the Csk αC helix region where its consequences for the kinetic behavior and domain displacement of Csk were verified through detailed molecular dynamics. The findings describe the mechanisms involving the Csk αC helix for the transitions and also support the dynamic correlation between SH3 and SH2 domains against the Csk lobes and how local energetic restrictions or interactions in the Csk αC helix can play an important role for long-range motions. The results also allow speculation if the Csk activity is restricted to one specific conformation or a consequence of a state transition, this point being a target for future studies. However, the αC helix is revealed as a potential region for rational drug design.

The ribosome is a complex biomolecular machine that utilizes large-scale conformational rearrangements to synthesize proteins. For example, during the elongation cycle, the "head" domain of the ribosomal small subunit (SSU) is known to undergo transient rotation events that allow for movement of tRNA molecules (i.e., translocation). While the head may exhibit rigid-body-like properties, the precise relationship between experimentally accessible probes and multidimensional rotations has yet to be established. To address this gap, we perform molecular dynamics simulations of the translocation step of the elongation cycle in the ribosome, where the SSU head spontaneously undergoes rotation and tilt-like motions. With this data set (1250 simulated events), we used statistical and information-theory-based measures to identify possible single-molecule probes that can isolate SSU head rotation and head tilting. This analysis provides a molecular interpretation for previous single-molecule measurements, while establishing a framework for the design of next-generation experiments that may precisely probe the mechanistic and kinetic aspects of the ribosome.

Protein synthesis involves a complex series of large-scale conformational changes in the ribosome. While long-lived intermediate states of these processes can be characterized by experiments, computational methods can be used to identify the interactions that contribute to the rate-limiting free-energy barriers. To this end, we use a simplified energetic model to perform molecular dynamics (MD) simulations of aminoacyl-tRNA (aa-tRNA) accommodation on the ribosome. While numerous studies have probed the energetics of the early stages of accommodation, we focus on the final stage of accommodation, where the 3'-CCA tail of aa-tRNA enters the peptidyl transferase center (PTC). These simulations show how a distinct intermediate is induced by steric confinement of the tail, immediately before it completes accommodation. Multiple pathways for 3'-CCA tail accommodation can be quantitatively distinguished, where the tail enters the PTC by moving past a pocket enclosed by Helix 89, 90, and 92, or through an alternate route formed by Helix 93 and the P-site tRNA. C2573, located within Helix 90, is shown to provide the largest contribution to this late-accommodation steric barrier, such that sub-Å perturbations to this residue can alter the time scale of tail accommodation by nearly an order of magnitude. In terms of biological function, these calculations suggest how this late-stage sterically induced barrier may contribute to tRNA proofreading by the ribosome.

Proper ionic concentrations are required for the functional dynamics of RNA and ribonucleoprotein (RNP) assemblies. While experimental and computational techniques have provided many insights into the properties of chelated ions, less is known about the energetic contributions of diffuse ions to large-scale conformational rearrangements. To address this, we present a model that is designed to quantify the influence of diffuse monovalent and divalent ions on the dynamics of biomolecular assemblies. This model employs all-atom (non-H) resolution and explicit ions, where effective potentials account for hydration effects. We first show that the model accurately predicts the number of excess Mg 2+ ions for prototypical RNA systems, at a level comparable to modern coarse-grained models. We then apply the model to a complete ribosome and show how the balance between diffuse Mg 2+ and K + ions can control the dynamics of tRNA molecules during translation. The model predicts differential effects of diffuse ions on the free-energy barrier associated with tRNA entry and the energy of tRNA binding to the ribosome. Together, this analysis reveals the direct impact of diffuse ions on the dynamics of an RNP assembly.

Evaluation of immunogenic epitopes for universal vaccine development in the face of ongoing SARS-CoV-2 evolution remains a challenge. Herein, we investigate the genetic and structural conservation of an immunogenically relevant epitope (C662-C671) of spike (S) protein across SARS-CoV-2 variants to determine its potential utility as a broad-spectrum vaccine candidate against coronavirus diseases. Comparative sequence analysis, structural assessment, and molecular dynamics simulations of C662-C671 epitope were performed. Mathematical tools were employed to determine its mutational cost. We found that the amino acid sequence of C662-C671 epitope is entirely conserved across the observed major variants of SARS-CoV-2 in addition to SARS-CoV. Its conformation and accessibility are predicted to be conserved, even in the highly mutated Omicron variant. Costly mutational rate in the context of energy expenditure in genome replication and translation can explain this strict conservation. These observations may herald an approach to developing vaccine candidates for universal protection against emergent variants of coronavirus., (© The Author(s) 2022. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Estrogen receptor α (ERα) is a regulatory protein that can access a set of distinct structural configurations. ERα undergoes extensive remodeling as it interacts with different agonists and antagonists, as well as transcription activation and repression factors. Moreover, breast cancer tumors resistant to hormone therapy have been associated with the imbalance between the active and inactive ERα states. Cancer-activating mutations in ERα play a crucial role in this imbalance and can promote the progression of cancer. However, the rate of this progression can also be increased by dysregulated pH in the tumor microenvironment. Many molecular aspects of the process of activation of ERα that can be affected by these pH changes and mutations are still unclear. Thus, we applied computational and experimental techniques to explore the activation process dynamics of ER for environments with different pHs and in the presence of one of the most recurrent cancer-activating mutations, D538G. Our results indicated that the effect of the pH increase associated with the D538G mutation promoted a robust stabilization of the active state of ER. We were also able to determine the main protein regions that have the most potential to influence the activation process under different pH conditions, which may provide targets of future therapeutics for the treatment of hormone-resistant breast cancer tumors. Finally, the approach used here can be applied for proteins associated with the proliferation of other cancer types, which can also have their function affected by small pH changes.

Applying simulations with structure-based G o ¯ - like models has proven to be an effective strategy for investigating the factors that control biomolecular dynamics. The common element of these models is that some (or all) of the intra/inter-molecular interactions are explicitly defined to stabilize an experimentally determined structure. To facilitate the development and application of this broad class of models, we previously released the SMOG 2 software package. This suite allows one to easily customize and distribute structure-based (i.e., SMOG) models for any type of polymer-ligand system. The force fields generated by SMOG 2 may then be used to perform simulations in highly optimized MD packages, such as Gromacs, NAMD, LAMMPS, and OpenMM. Here, we describe extensions to the software and demonstrate the capabilities of the most recent version (SMOG v2.4.2). Changes include new tools that aid user-defined customization of force fields, as well as an interface with the OpenMM simulation libraries (OpenSMOG v1.1.0). The OpenSMOG module allows for arbitrary user-defined contact potentials and non-bonded potentials to be employed in SMOG models, without source-code modifications. To illustrate the utility of these advances, we present applications to systems with millions of atoms, long polymers and explicit ions, as well as models that include non-structure-based (e.g., AMBER-based) energetic terms. Examples include large-scale rearrangements of the SARS-CoV-2 Spike protein, the HIV-1 capsid with explicit ions, and crystallographic lattices of ribosomes and proteins. In summary, SMOG 2 and OpenSMOG provide robust support for researchers who seek to develop and apply structure-based models to large and/or intricate biomolecular systems., (© 2021 The Protein Society.)

Protein synthesis in the cell is controlled by an elaborate sequence of conformational rearrangements in the ribosome. The composition of a ribosome varies by species, though they typically contain ∼ 50-100 RNA and protein molecules. While advances in structural techniques have revolutionized our understanding of long-lived conformational states, a vast range of transiently visited configurations can not be directly observed. In these cases, computational/simulation methods can be used to understand the mechanical properties of the ribosome. Insights from these approaches can then help guide next-generation experimental measurements. In this short review, we discuss theoretical strategies that have been deployed to quantitatively describe the energetics of collective rearrangements in the ribosome. We focus on efforts to probe large-scale subunit rotation events, which involve the coordinated displacement of large numbers of atoms (tens of thousands). These investigations are revealing how the molecular structure of the ribosome encodes the mechanical properties that control large-scale dynamics., Competing Interests: Conflict of interestThe authors declare no competing interests., (© International Union for Pure and Applied Biophysics (IUPAB) and Springer-Verlag GmbH Germany, part of Springer Nature 2021.)

In this session, experts in molecular biophysics described the dynamics of biopolymers across a wide range of length and time scales. This discussion highlighted numerous techniques that span from highly detailed simulations, to coarse-grained theoretical models, as well as high-resolution structural analysis. The topics were equally diverse, where there was discussion of biological processes at small (individual atoms), intermediate (assemblies) and very large scales (phase separation)., Competing Interests: Conflict of interestThe authors declare no competing interests., (© International Union for Pure and Applied Biophysics (IUPAB) and Springer-Verlag GmbH Germany, part of Springer Nature 2021.)

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is highly contagious, and transmission involves a series of processes that may be targeted by vaccines and therapeutics. During transmission, host cell invasion is controlled by a large-scale (200-300 Å) conformational change of the Spike protein. This conformational rearrangement leads to membrane fusion, which creates transmembrane pores through which the viral genome is passed to the host. During Spike-protein-mediated fusion, the fusion peptides must be released from the core of the protein and associate with the host membrane. While infection relies on this transition between the prefusion and postfusion conformations, there has yet to be a biophysical characterization reported for this rearrangement. That is, structures are available for the endpoints, though the intermediate conformational processes have not been described. Interestingly, the Spike protein possesses many post-translational modifications, in the form of branched glycans that flank the surface of the assembly. With the current lack of data on the pre-to-post transition, the precise role of glycans during cell invasion has also remained unclear. To provide an initial mechanistic description of the pre-to-post rearrangement, an all-atom model with simplified energetics was used to perform thousands of simulations in which the protein transitions between the prefusion and postfusion conformations. These simulations indicate that the steric composition of the glycans can induce a pause during the Spike protein conformational change. We additionally show that this glycan-induced delay provides a critical opportunity for the fusion peptides to capture the host cell. In contrast, in the absence of glycans, the viral particle would likely fail to enter the host. This analysis reveals how the glycosylation state can regulate infectivity, while providing a much-needed structural framework for studying the dynamics of this pervasive pathogen., Competing Interests: ED, JO, PW No competing interests declared, (© 2021, Dodero-Rojas et al.)

Development of effective vaccines against coronavirus disease 2019 (COVID-19) is a global imperative. Rapid immunization of the entire human population against a widespread, continually evolving, and highly pathogenic virus is an unprecedented challenge, and different vaccine approaches are being pursued. Engineered filamentous bacteriophage (phage) particles have unique potential in vaccine development due to their inherent immunogenicity, genetic plasticity, stability, cost-effectiveness for large-scale production, and proven safety profile in humans. Herein we report the development and initial evaluation of two targeted phage-based vaccination approaches against SARS-CoV-2: dual ligand peptide-targeted phage and adeno-associated virus/phage (AAVP) particles. For peptide-targeted phage, we performed structure-guided antigen design to select six solvent-exposed epitopes of the SARS-CoV-2 spike (S) protein. One of these epitopes displayed on the major capsid protein pVIII of phage induced a specific and sustained humoral response when injected in mice. These phage were further engineered to simultaneously display the peptide CAKSMGDIVC on the minor capsid protein pIII to enable their transport from the lung epithelium into the systemic circulation. Aerosolization of these "dual-display" phage into the lungs of mice generated a systemic and specific antibody response. In the second approach, targeted AAVP particles were engineered to deliver the entire S protein gene under the control of a constitutive CMV promoter. This induced tissue-specific transgene expression, stimulating a systemic S protein-specific antibody response in mice. With these proof-of-concept preclinical experiments, we show that both targeted phage- and AAVP-based particles serve as robust yet versatile platforms that can promptly yield COVID-19 vaccine prototypes for translational development., Competing Interests: Competing interest statement: D.I.S., W.A., and R.P. are listed as inventors on a patent application related to this technology (International Patent Application no. PCT/US2020/053758, entitled “Targeted Pulmonary Delivery Compositions and Methods Using Same”). Provisional patent application nos. 63/048, 279, and 63/161,136, entitled “Enhancing Immune Responses Through Targeted Antigen Expression,” have also been filed on the technology and intellectual property reported here. PhageNova Bio has licensed these intellectual properties and D.I.S., F.H.F.T., C.M., V.J.Y., T.L.S., S.K.L., W.A., and R.P. may be entitled to standard royalties. R.P., S.K.L., and W.A. are founders and equity stockholders of PhageNova Bio. S.K.L. is a board member and R.P. is chief scientific officer and a paid consultant of PhageNova Bio. V.J.Y. is currently a full-time employee of PhageNova Bio. R.P. and W.A. are founders and equity shareholders of MBrace Therapeutics; R.P. serves as the chief scientific officer and W.A. is a member of the scientific advisory board at MBrace Therapeutics. F.I.S. is currently a full-time employee of MBrace Therapeutics. These arrangements are managed in accordance with the established institutional conflict-of-interest policies of Rutgers, The State University of New Jersey., (Copyright © 2021 the Author(s). Published by PNAS.)

Protein synthesis by the ribosome is coordinated by an intricate series of large-scale conformational rearrangements. Structural studies can provide information about long-lived states, however biological kinetics are controlled by the intervening free-energy barriers. While there has been progress describing the energy landscapes of bacterial ribosomes, very little is known about the energetics of large-scale rearrangements in eukaryotic systems. To address this topic, we constructed an all-atom model with simplified energetics and performed simulations of subunit rotation in the yeast ribosome. In these simulations, the small subunit (SSU; ~1MDa) undergoes spontaneous and reversible rotations (~8°). By enabling the simulation of this rearrangement under equilibrium conditions, these calculations provide initial insights into the molecular factors that control dynamics in eukaryotic ribosomes. Through this, we are able to identify specific inter-subunit interactions that have a pronounced influence on the rate-limiting free-energy barrier. We also show that, as a result of changes in molecular flexibility, the thermodynamic balance between the rotated and unrotated states is temperature-dependent. This effect may be interpreted in terms of differential molecular flexibility within the rotated and unrotated states. Together, these calculations provide a foundation, upon which the field may begin to dissect the energetics of these complex molecular machines., Competing Interests: Conflicts of Interest: The authors declare no conflict of interest.

The ribosome is a biomolecular machine that undergoes multiple large-scale structural rearrangements during protein elongation. Here, we focus on a conformational rearrangement during translocation, known as P/E hybrid-state formation. Using a model that explicitly represents all non-hydrogen atoms, we simulated more than 120 spontaneous transitions, where the tRNA molecule is displaced between the P and E sites of the large subunit. In addition to predicting a free-energy landscape that is consistent with previous experimental observations, the simulations reveal how a six-residue gate-like region can limit P/E formation, where sub-angstrom structural perturbations lead to an order-of-magnitude change in kinetics. Thus, this precisely defined set of residues represents a novel target that may be used to control functional dynamics in bacterial ribosomes. This theoretical analysis establishes a direct relationship between ribosome structure and large-scale dynamics, and it suggests how next-generation experiments may precisely dissect the energetics of hybrid formation on the ribosome.

Molecular dynamics techniques provide numerous strategies for investigating biomolecular energetics, though quantitative analysis is often only accessible for relatively small (frequently monomeric) systems. To address this limit, we use simulations in combination with a simplified energetic model to study complex rearrangements in a large assembly. We use cryo-EM reconstructions to simulate the DNA packaging-associated 3 nm expansion of the protein shell of an initially assembled phage T7 capsid (called procapsid or capsid I). This is accompanied by a disorder-order transition and expansion-associated externalization displacement of the 420 N-terminal tails of the shell proteins. For the simulations, we use an all-atom structure-based model (1.07 million atoms), which is specifically designed to probe the influence of molecular sterics on dynamics. We find that the rate at which the N-terminal tails undergo translocation depends heavily on their position within hexons and pentons. Specifically, trans-shell displacements of the hexon E subunits are the most frequent and hexon A subunits are the least frequent. The simulations also implicate numerous tail translocation intermediates during tail translocation that involve topological traps, as well as sterically induced barriers. The presented study establishes a foundation for understanding the precise relationship between molecular structure and phage maturation.

Modifications in the tRNA anticodon loop, adjacent to the three-nucleotide anticodon, influence translation fidelity by stabilizing the tRNA to allow for accurate reading of the mRNA genetic code. One example is the N1-methylguanosine modification at guanine nucleotide 37 (m 1 G37) located in the anticodon loop andimmediately adjacent to the anticodon nucleotides 34, 35, 36. The absence of m 1 G37 in tRNA Pro causes +1 frameshifting on polynucleotide, slippery codons. Here, we report structures of the bacterial ribosome containing tRNA Pro bound to either cognate or slippery codons to determine how the m 1 G37 modification prevents mRNA frameshifting. The structures reveal that certain codon-anticodon contexts and the lack of m 1 G37 destabilize interactions of tRNA Pro with the P site of the ribosome, causing large conformational changes typically only seen during EF-G-mediated translocation of the mRNA-tRNA pairs. These studies provide molecular insights into how m 1 G37 stabilizes the interactions of tRNA Pro with the ribosome in the context of a slippery mRNA codon., Competing Interests: EH, SH, SS, TM, RG, PW, CD No competing interests declared, (© 2020, Hoffer et al.)

Electrokinetic translocation of biomolecules through solid-state nanopores represents a label-free single-molecule technique that may be used to measure biomolecular structure and dynamics. Recent investigations have attempted to distinguish individual transfer RNA (tRNA) species based on the associated pore translocation times, ion-current noise, and blockage currents. By manufacturing sufficiently smaller pores, each tRNA is required to undergo a deformation to translocate. Accordingly, differences in nanopore translocation times and distributions may be used to infer the mechanical properties of individual tRNA molecules. To bridge our understanding of tRNA structural dynamics and nanopore measurements, we apply molecular dynamics simulations using a simplified "structure-based" energetic model. Calculating the free-energy landscape for distinct tRNA species implicates transient unfolding of the terminal RNA helix during nanopore translocation. This provides a structural and energetic framework for interpreting current experiments, which can aid the design of methods for identifying macromolecules using nanopores., (Copyright © 2020 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Understanding how enzymes achieve their tremendous catalytic power is a major question in biochemistry. Greater understanding is also needed for enzyme engineering applications. In many cases, enzyme efficiency and specificity depend on residues not in direct contact with the substrate, termed remote residues. This work focuses on Escherichia coli ornithine transcarbamoylase (OTC), which plays a central role in amino acid metabolism. OTC has been reported to undergo an induced-fit conformational change upon binding its first substrate, carbamoyl phosphate (CP), and several residues important for activity have been identified. Using computational methods based on the computed chemical properties from theoretical titration curves, sequence-based scores derived from evolutionary history, and protein surface topology, residues important for catalytic activity were predicted. The roles of these residues in OTC activity were tested by constructing mutations at predicted positions, followed by steady-state kinetics assays and substrate binding studies with the variants. First-layer mutations R57A and D231A, second-layer mutation H272L, and third-layer mutation E299Q, result in 57- to 450-fold reductions in kcat/KM with respect to CP and 44- to 580-fold reductions with respect to ornithine. Second-layer mutations D140N and Y160S also reduce activity with respect to ornithine. Most variants had decreased stability relative to wild-type OTC, with variants H272L, H272N, and E299Q having the greatest decreases. Variants H272L, E299Q, and R57A also show compromised CP binding. In addition to direct effects on catalytic activity, effects on overall protein stability and substrate binding were observed that reveal the intricacies of how these residues contribute to catalysis., Competing Interests: The authors have declared that no competing interests exist.

Protein folding occurs in a high dimensional phase space, and the representation of the associated energy landscape is nontrivial. A widely applied approach to studying folding landscapes is to describe the dynamics along a small number of reaction coordinates. However, other strategies involve more elaborate analysis of the complex phase space. There have been many attempts to obtain a more detailed representation of all available conformations for a given system. In this work, we address this problem using a metric based on internal distances between amino acids to describe the differences between any two conformations. Using an effective projection method, we are able to go beyond the typical one-dimensional representation and provide intuitive two dimensional visualizations of the landscape. We refer to this method as the energy landscape visualization method (ELViM). We have applied this methodology using a C α structure-based model to study the folding of two well-known proteins: SH3 domain and protein-A. Our visualization method yields a detailed description of the folding process, making possible the identification of transition state regions, and establishing the paths that lead to the native state. For SH3, we have analyzed structural differences in the distribution of folding routes. The competition between the native and mirror structures in protein A is also discussed. Finally, the method is applied to study conformational changes in the protein elongation factor thermally unstable. Distinct features of ELViM are that it does not require or assume a reaction coordinate, and it does not require analysis of kinetic aspects of the system.

The stochastic drift-diffusion (DrDiff) theory is an approach used to characterize the dynamical properties of simulation data. With new features in transition times analyses, the framework characterized the thermodynamic free-energy profile [F(Q)], the folding time (τ f ), and transition path time (τ TP ) by determining the coordinate-dependent drift-velocity [v(Q)] and diffusion [D(Q)] coefficients from trajectory time traces. In order to explore the DrDiff approach and to tune it with two other methods (Bayesian analysis and fep1D algorithm), a numerical integration of the Langevin equation with known D(Q) and F(Q) was performed and the inputted coefficients were recovered with success by the diffusion models. DrDiff was also applied to investigate the prion protein (PrP) kinetics and thermodynamics by analyzing folding/unfolding simulations. The protein structure-based model, the well-known Go¯-model, was employed in a coarse-grained C α level to generate long constant-temperature time series. PrP was chosen due to recent experimental single-molecule studies in D and τ TP that stressed the importance and the difficulty of probing these quantities and the rare transition state events related to prion misfolding and aggregation. The PrP thermodynamic double-well F(Q) profile, the "X" shape of τ f (T), and the linear shape of τ TP (T) were predicted with v(Q) and D(Q) obtained by the DrDiff algorithm. With the advance of single-molecule techniques, the DrDiff framework might be a useful ally for determining kinetic and thermodynamic properties by analyzing time observables of biomolecular systems. The code is freely available at https://github.com/ronaldolab/DrDiff.

In recent years, there has been a growing interest to quantify the energy landscape that governs ribosome dynamics. However, in order to quantitatively integrate theoretical predictions and experimental measurements, it is essential that one has a detailed understanding of the associated diffusive properties. Here, all-atom explicit-solvent simulations (50 μs of aggregate sampling) predict that the diffusion coefficient of a tRNA molecule will depend on its position within the ribosome. Specifically, during aa-tRNA accommodation (i.e., the process by which tRNA enters the ribosome), the apparent diffusion coefficient decreases by approximately an order of magnitude. By comparing these to values obtained with an energetically "smooth" model, we show that the observed nonuniform behavior likely arises from electrostatic and solvation interactions between the tRNA and ribosome. These calculations also reveal the hierarchical character of ribosomal energetics, where steric interactions induce a large-scale free-energy barrier, and short-scale roughness determines the rate of diffusive movement across the landscape.

With the broad accessibility of high-performance computing resources, the significance of a molecular dynamics simulation is now rarely limited by hardware and/or software availability. Rather, the scientific value of each calculation is determined by the principles that underlie the theoretical model. The current review addresses this topic in the context of simplified models applied to large-scale (∼20-100 Å) dynamics in the ribosome. Specifically, we focus on applications of the "SMOG" class of structure-based models, which can be used to simulate spontaneous (i.e. non-targeted) conformational rearrangements in complex molecular assemblies. Here, we aim to provide an entry-level assessment of the methods, which can help bridge conceptual and communication gaps between the experimental and computational communities. In addition, inspecting the strategies that have been deployed previously can provide guidelines for future computational investigations into the relationship between structure, energetics and dynamics in other assemblies., (Copyright © 2019 Elsevier Inc. All rights reserved.)

The accurate expression of proteins requires the ribosome to efficiently undergo elaborate conformational rearrangements. The most dramatic of these motions is subunit rotation, which is necessary for tRNA molecules to transition between ribosomal binding sites. While rigid-body descriptions provide a qualitative picture of the process, obtaining quantitative mechanistic insights requires one to account for the relationship between molecular flexibility and collective dynamics. Using simulated rotation events, we assess the quality of experimentally accessible measures for describing the collective displacement of the ∼4000-residue small subunit. For this, we ask whether each coordinate is able to identify the underlying free-energy barrier and transition state ensemble (TSE). We find that intuitive structurally motivated coordinates (e.g., rotation angle, interprotein distances) can distinguish between the endpoints, though they are poor indicators of barrier-crossing events, and they underestimate the free-energy barrier. In contrast, coordinates based on intersubunit bridges can identify the TSE. We additionally verify that the committor probability for the putative TSE configurations is 0.5, a hallmark feature of any transition state. In terms of structural properties, these calculations implicate a transition state in which flexibility allows for asynchronous rearrangements of the bridges, as the ribosome adopts a partially rotated orientation. This provides a theoretical foundation, upon which experimental techniques may precisely quantify the energy landscape of the ribosome.

Over the last 20 years, the application of structure-based (Gō-like) models has ranged from protein folding with coarse-grained models to all-atom representations of large-scale molecular assemblies. While there are many variants that may be employed, the common feature of these models is that some (or all) of the stabilizing energetic interactions are defined based on the knowledge of a particular experimentally obtained conformation. With the generality of this approach, there was a need for a versatile computational platform for designing and implementing this class of models. To this end, the SMOG 2 software package provides an easy-to-use interface, where the user has full control of the model parameters. This software allows the user to edit XML-formatted files in order to provide definitions of new structure-based models. SMOG 2 reads these "template" files and maps the interactions onto specific structures, which are provided in PDB format. The force field files produced by SMOG 2 may then be used to perform simulations with a variety of popular molecular dynamics suites. In this chapter, we describe some of the key features of the SMOG 2 package, while providing examples and strategies for applying these techniques to complex (often large-scale) molecular assemblies, such as the ribosome.

Elongation factor Tu (EF-Tu) is a three-domain protein that is responsible for delivering aminoacyl-tRNA (aa-tRNA) molecules to the ribosome. During the delivery process, EF-Tu undergoes a large-scale (~50Å) conformational transition that results in rearrangement of domain I, relative to the II/III superdomain. Despite the central role of EF-Tu during protein synthesis, little is known about the structural and energetic properties of this reordering process. To study the physical-chemical properties of domain motion, we constructed a multi-basin structure-based (i.e., Gō-like) model, with which we have simulated hundreds of spontaneous conformational rearrangements. By analyzing the statistical properties of these events, we show that EF-Tu is likely to adopt a disordered intermediate ensemble during this transition. We further show that this disordered intermediate will favor a specific sequence of conformational substeps when bound to the ribosome, and the disordered ensemble can influence the kinetics of the incoming aa-tRNA molecule. Overall, this study highlights the dynamic nature of EF-Tu by revealing a relationship between conformational disorder and biological function., (© 2018 Wiley Periodicals, Inc.)

A major challenge in the study of biomolecular assemblies is to identify reaction coordinates that precisely monitor conformational rearrangements. This is central to the interpretation of single-molecule fluorescence resonance energy transfer measurements, where the observed dynamics depends on the labeling strategy. As an example, different probes of subunit rotation in the ribosome have provided qualitatively distinct descriptions. In one study, changes in fluorescence suggested that the 30S body undergoes a single rotation/back-rotation cycle during the process of mRNA-tRNA translocation. In contrast, an alternate assay implicated the presence of reversible rotation events before completing translocation. For future single-molecule experiments to unambiguously measure the relationship between subunit rotation and translocation, it is necessary to rationalize these conflicting descriptions. To this end, we have simulated hundreds of spontaneous subunit rotation events (≈8°) using a residue-level coarse-grained model of the ribosome. We analyzed nine different reaction coordinates and found that the apparently inconsistent measurements are likely a consequence of ribosomal flexibility. Further, we propose a metric for quantifying the degree of energetic coupling between experimentally measured degrees of freedom and subunit rotation. This analysis provides a physically grounded framework that can guide the development of more precise single-molecule techniques., (Copyright © 2017 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

The ribosome is a large ribonucleoprotein complex that is responsible for the production of proteins in all organisms. Accommodation is the process by which an incoming aminoacyl-tRNA (aa-tRNA) molecule binds the ribosomal A site, and its kinetics has been implicated in the accuracy of tRNA selection. In addition to rearrangements in the aa-tRNA molecule, the L11 stalk can undergo large-scale anisotropic motions during translation. To explore the potential impact of this protruding region on the rate of aa-tRNA accommodation, we used molecular dynamics simulations with a simplified model to evaluate the free energy as a function of aa-tRNA position. Specifically, these calculations describe the transition between A/T and elbow-accommodated (EA) configurations (∼20 Å displacement). We find that the free-energy barrier associated with elbow accommodation is proportional to the degree of mobility exhibited by the L11 stalk. That is, when L11 is more rigid, the free-energy barrier height is decreased. This effect arises from the ability of L11 to confine, and thereby destabilize, the A/T ensemble. In addition, when elongation factor Tu (EF-Tu) is present, the A/T ensemble is further destabilized in an L11-dependent manner. These results provide a framework that suggests how next-generation experiments may precisely control the dynamics of the ribosome.

Proteins are structurally dynamic macromolecules, and it is challenging to quantify the conformational properties of their native state in solution. Nanopores can be efficient tools to study proteins in a solution environment. In this method, an electric field induces electrophoretic and/or electro-osmotic transport of protein molecules through a nanopore slightly larger than the protein molecule. High-bandwidth ion current measurement is used to detect the transit of each protein molecule. First, our measurements reveal a correlation between the mean current blockade amplitude and the radius of gyration for each protein. Next, we find a correlation between the shape of the current signal amplitude distributions and the protein fluctuation as obtained from molecular dynamics simulations. Further, the magnitude of the structural fluctuations, as probed by experiments and simulations, correlates with the ratio of α-helix to β-sheet content. We highlight the resolution of our measurements by resolving two states of calmodulin, a canonical protein that undergoes a conformational change in response to calcium binding.

Proteins are synthesized by the joint action of the ribosome and tRNA molecules, where the rate of synthesis can be affected by numerous factors, such as the concentration of tRNA, the binding affinity of tRNA for the ribosome, or post-transcriptional modifications. Here, we expand this range of contributors by demonstrating how differences in tRNA structure can give rise to tRNA species-specific dynamics in the ribosome. To show this, we perform simulations of A/P hybrid-state formation for two tRNA species (tRNA Phe and tRNA Leu ), which differ in the size of their variable loops (VLs). These calculations reveal that steric interactions between the VL and the ribosomal A-site finger (ASF, i.e., H38 of 23S rRNA) can directly modulate the free-energy landscape for each tRNA species. We also find that tRNA and ASF motions are highly correlated, where fluctuations of the ASF are predictive of tRNA transition events. Finally, by introducing perturbations to the model, we demonstrate that ASF flexibility is a determinant of the rate of A/P hybrid-state formation.

For decades, protein folding and functional dynamics have been described in terms of diffusive motion across an underlying energy landscape. With continued advances in structural biology and high-performance computing, the field is positioned to extend these approaches to large biomolecular assemblies. Through the application of energy landscape techniques to the ribosome, one may work towards establishing a comprehensive description of the dynamics, which will bridge theoretical concepts and experimental observations. In this perspective, we discuss a few of the challenges that will need to be addressed as we extend the application of landscape principles to the ribosome.

It has long been recognized that the thermodynamics of mRNA-tRNA base pairing is insufficient to explain the high fidelity and efficiency of aminoacyl-tRNA (aa-tRNA) selection by the ribosome. To rationalize this apparent inconsistency, Hopfield proposed that the ribosome may improve accuracy by utilizing a multi-step kinetic proofreading mechanism. While biochemical, structural and single-molecule studies have provided a detailed characterization of aa-tRNA selection, there is a limited understanding of how the physical-chemical properties of the ribosome enable proofreading. To this end, we probe the role of EF-Tu during aa-tRNA accommodation (the proofreading step) through the use of energy landscape principles, molecular dynamics simulations and kinetic models. We find that the steric composition of EF-Tu can reduce the free-energy barrier associated with the first step of accommodation: elbow accommodation. We interpret this effect within an extended kinetic model of accommodation and show how EF-Tu can contribute to efficient and accurate proofreading.

In order to quantitatively describe the energetics of biomolecular rearrangements, it is necessary to identify reaction coordinates that accurately capture the relevant transition events. Here, we perform simulations of A-site tRNA movement (∼20 Å) during hybrid-state formation in the ribosome and quantify the ability of interatomic distances to capture the transition state ensemble. Numerous coordinates are found to be accurate indicators of the transition state, allowing tRNA rearrangements to be described as diffusion across a one-dimensional free-energy surface. In addition to providing insights into the physical-chemical relationship between biomolecular structure and dynamics, these results can help enable single-molecule techniques to probe the free-energy landscape of the ribosome.

Molecular dynamics simulations with coarse-grained or simplified Hamiltonians have proven to be an effective means of capturing the functionally important long-time and large-length scale motions of proteins and RNAs. Originally developed in the context of protein folding, structure-based models (SBMs) have since been extended to probe a diverse range of biomolecular processes, spanning from protein and RNA folding to functional transitions in molecular machines. The hallmark feature of a structure-based model is that part, or all, of the potential energy function is defined by a known structure. Within this general class of models, there exist many possible variations in resolution and energetic composition. SMOG 2 is a downloadable software package that reads user-designated structural information and user-defined energy definitions, in order to produce the files necessary to use SBMs with high performance molecular dynamics packages: GROMACS and NAMD. SMOG 2 is bundled with XML-formatted template files that define commonly used SBMs, and it can process template files that are altered according to the needs of each user. This computational infrastructure also allows for experimental or bioinformatics-derived restraints or novel structural features to be included, e.g. novel ligands, prosthetic groups and post-translational/transcriptional modifications. The code and user guide can be downloaded at http://smog-server.org/smog2.

Translocation of mRNA and tRNA through the ribosome is associated with large-scale rearrangements of the head domain in the 30S ribosomal subunit. To elucidate the relationship between 30S head dynamics and mRNA-tRNA displacement, we apply molecular dynamics simulations using an all-atom structure-based model. Here we provide a statistical analysis of 250 spontaneous transitions between the A/P-P/E and P/P-E/E ensembles. Consistent with structural studies, the ribosome samples a chimeric ap/P-pe/E intermediate, where the 30S head is rotated ∼18°. It then transiently populates a previously unreported intermediate ensemble, which is characterized by a ∼10° tilt of the head. To identify the origins of head tilting, we analyse 781 additional simulations in which specific steric features are perturbed. These calculations show that head tilting may be attributed to specific steric interactions between tRNA and the 30S subunit (PE loop and protein S13). Taken together, this study demonstrates how molecular structure can give rise to large-scale collective rearrangements.

The ever-increasing capacity of computing resources has extended ribosome calculations from the study of small-scale fluctuations to large-scale barrier-crossing processes. As the field of computational/theoretical biophysics shifts focus to large-scale conformational transitions, there is a growing need for a systematic framework to interpret and analyze ribosome dynamics. To this end, energy landscape principles, largely developed for the study of biomolecular folding, have proven to be invaluable. These tools not only provide a foundation for describing simulations but can be used to reconcile experimental results, as well. In this review, I will discuss recent efforts to employ computational methods to reveal the characteristics of the ribosome's landscape and how these studies can help guide a new generation of experiments that more closely probe the underlying energetics. As a result of these investigations, general principles about ribosome function are beginning to emerge, including that: (1) small-scale fluctuations are the result of structure, rather than detailed energetics, (2) molecular flexibility leads to entropically favored rearrangements, and (3) tRNA dynamics may be accurately described as diffusive movement across an energy landscape.

RNA is highly sensitive to the ionic environment and typically requires Mg(2+) to form compact structures. There is a need for models capable of describing the ion atmosphere surrounding RNA with quantitative accuracy. We present a model of RNA electrostatics and apply it within coarse-grained molecular dynamics simulation. The model treats Mg(2+) ions explicitly to account for ion-ion correlations neglected by mean-field theories. Since mean-field theories capture KCl well, it is treated implicitly by a generalized Manning counterion condensation model. The model extends Manning condensation to deal with arbitrary RNA conformations, nonlimiting KCl concentrations, and the ion inaccessible volume of RNA. The model is tested against experimental measurements of the excess Mg(2+) associated with the RNA, Γ(2+), because Γ(2+) is directly related to the Mg(2+)-RNA interaction free energy. The excellent agreement with experiment demonstrates that the model captures the ionic dependence of the RNA free energy landscape.

As our understanding of biological dynamics continues to be refined, it is becoming clear that biomolecules can undergo transitions between ordered and disordered states as they execute functional processes. From a computational perspective, studying disorder events poses a challenge, as they typically occur on long timescales, and the associated molecules are often large (i.e., hundreds of residues). These size and time requirements make it advantageous to use computationally inexpensive models to characterize large-scale dynamics, where more highly detailed models can provide information about individual sub-steps associated with function. To reduce computational demand, one often uses a coarse-grained representation of the molecule or a simplified description of the energetics. In order to use simpler models to identify transient disorder in RNA and proteins, it is imperative that these models can accurately capture structural fluctuations about folded configurations, as well as the overall stability of each molecule. Here, we explore a class of simplified model for which all non-hydrogen atoms are explicitly represented. We find that this model can provide a consistent description of protein folding and native-basin dynamics for several representative biomolecules. We additionally show that the native-basin fluctuations of tRNA and the ribosome are robust to variations in the model. Finally, the extended variable loop in tRNAIle is predicted to be very dynamic, which may facilitate biologically-relevant rearrangements. Together, this study provides a foundation that will aid in the application of simplified models to study disorder during function in ribonucleoprotein (RNP) assemblies.

Protein folding was the first area of molecular biology for which a systematic statistical-mechanical analysis of dynamics was developed. As a result, folding is described as a process by which a disordered protein chain diffuses across a high-dimensional energy landscape and finally reaches the folded ensemble. Folding studies have produced countless theoretical concepts that are generalizable to other biomolecular processes, such as the functional dynamics of molecular assemblies. Common themes in folding and function include the dominant role of excluded volume, that a balance between energetic roughness and geometrical effects guides dynamics, and that folding/functional landscapes are relatively smooth. Here, we discuss how insights into protein folding have been applied to investigate the functional dynamics of biomolecular assemblies., (Copyright © 2014 Elsevier Ltd. All rights reserved.)

To reveal the molecular determinants of biological function, one seeks to characterize the interactions that are formed in conformational and chemical transition states. In other words, what interactions govern the molecule's energy landscape? To accomplish this, it is necessary to determine which degrees of freedom can unambiguously identify each transition state. Here, we perform simulations of large-scale aminoacyl-transfer RNA (aa-tRNA) rearrangements during accommodation on the ribosome and project the dynamics along experimentally accessible atomic distances. From this analysis, we obtain evidence for which coordinates capture the correct number of barrier-crossing events and accurately indicate when the aa-tRNA is on a transition path. Although a commonly used coordinate in single-molecule experiments performs poorly, this study implicates alternative coordinates along which rearrangements are accurately described as diffusive movements across a one-dimensional free-energy profile. From this, we provide the theoretical foundation required for single-molecule techniques to uncover the energy landscape governing aa-tRNA selection by the ribosome., (Copyright © 2014 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Influenza hemagglutinin (HA), a homotrimeric glycoprotein crucial for membrane fusion, undergoes a large-scale structural rearrangement during viral invasion. X-ray crystallography has shown that the pre- and postfusion configurations of HA2, the membrane-fusion subunit of HA, have disparate secondary, tertiary, and quaternary structures, where some regions are displaced by more than 100 Å. To explore structural dynamics during the conformational transition, we studied simulations of a minimally frustrated model based on energy landscape theory. The model combines structural information from both the pre- and postfusion crystallographic configurations of HA2. Rather than a downhill drive toward formation of the central coiled-coil, we discovered an order-disorder transition early in the conformational change as the mechanism for the release of the fusion peptides from their burial sites in the prefusion crystal structure. This disorder quickly leads to a metastable intermediate with a broken threefold symmetry. Finally, kinetic competition between the formation of the extended coiled-coil and C-terminal melting results in two routes from this intermediate to the postfusion structure. Our study reiterates the roles that cracking and disorder can play in functional molecular motions, in contrast to the downhill mechanical interpretations of the "spring-loaded" model proposed for the HA2 conformational transition.

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., (Copyright © 2014 Biophysical Society. Published by Elsevier Inc. All rights reserved.)