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.

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.

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.

No abstract