Summary Nearly every enzyme undergoes a significant change in structure after binding it’s substrate. New experimental and theoretical analyses of the role of changes in HIV reverse transcriptase structure in selecting a correct substrate are presented. Atomically detailed simulations using the Milestoning method predict a rate and free energy profile of the conformational change commensurate with experimental data. A large conformational change occurring on a ms timescale locks the correct nucleotide at the active site, but promotes release of a mismatched nucleotide. The positions along the reaction coordinate that decide the yield of the reaction are not determined by the chemical step. Rather, the initial steps of weak substrate binding and protein conformational transition significantly enrich the yield of a reaction with a correct substrate, while the same steps diminish the reaction probability of an incorrect substrate.

The method of Directional Milestoning is revisited. We start from an exact and more general expression and state the conditions and validity of the memory-loss approximation. An algorithm to compute a reaction coordinate from Directional Milestoning data is presented. The reaction coordinate is calculated as a set of discrete jumps between Milestones that maximizes the flux between two stable states. As an application we consider a conformational transition in solvated Adenosine. We compare a long molecular dynamic trajectory with Directional Milestoning and discuss the differences between the maximum flux path and minimum energy coordinates.

We used single-channel electrical recordings and Langevin molecular dynamics simulations to explore the electrophoretic translocation of various beta-hairpin peptides across the staphylococcal alpha-hemolysin (alphaHL) protein pore at single-molecule resolution. The beta-hairpin peptides, which varied in their folding properties, corresponded to the C terminal residues of the B1 domain of protein G. The translocation time was strongly dependent on the electric force and was correlated with the folding features of the beta-hairpin peptides. Highly unfolded peptides entered the pore in an extended conformation, resulting in fast single-file translocation events. In contrast, the translocation of the folded beta-hairpin peptides occurred more slowly. In this case, the beta-hairpin peptides traversed the alphaHL pore in a misfolded or fully folded conformation. This study demonstrates that the interaction between a polypeptide and a beta-barrel protein pore is dependent on the folding features of the polypeptide.

RNA molecules play critical roles in many cellular processes. Traditionally viewed as genetic messengers, RNA molecules were recently discovered to have diverse functions related to gene regulation and expression. RNA also has great potential as a therapeutic and a tool for further investigation of gene regulation. Metal ions are an integral part of RNA structure and should be considered in any experimental or theoretical study of RNA. Here, we report a multidisciplinary approach that combines anomalous small-angle x-ray scattering and molecular-dynamics (MD) simulations with explicit solvent and ions around RNA. From experiment and simulation results, we find excellent agreement in the number and distribution of excess monovalent and divalent ions around a short RNA duplex. Although similar agreement can be obtained from a continuum description of the solvent and mobile ions (by solving the Poisson-Boltzmann equation and accounting for finite ion size), the use of MD is easily extended to flexible RNA systems with thermal fluctuations. Therefore, we also model a short RNA pseudoknot and find good agreement between the MD results and the experimentally derived solution structures. Surprisingly, both deviate from crystal structure predictions. These favorable comparisons of experiment and simulations encourage work on RNA in all-atom dynamic models.

Atomically detailed distributions of ions around an A-form RNA are computed. Different mixtures of monovalent and divalent ions are considered explicitly. Studies of tightly bound and of diffusive (but bound) ions around 25 base pairs RNA are conducted in explicit solvent. Replica exchange simulations provide detailed equilibrium distributions with moderate computing resources (20 nanoseconds of simulation using 64 replicas). The simulations show distinct behavior of single and doubly charged cations. Binding of Mg2+ ion includes tight binding to specific sites while Na+ binds only diffusively. The tight binding of Mg2+ is with a solvation shell while Na+ can bind directly to RNA. Negative mobile ions can be found near the RNA but must be assisted by proximate and mobile cations. At distances larger than 16Å from the RNA center, a model of RNA as charged rod in a continuum of ionic solution provides quantitative description of the ion density (the same as in atomically detailed simulation). At shorter distances, the structure of RNA (and ions) have significant impact on the pair correlation functions. Predicted binding sites of Mg2+ at the RNA surface are in accord with structures from crystallography. Electric field relaxation is investigated. The relaxation due to solution rearrangements is completed in tens of picoseconds, while the contribution of RNA tumbling continues to a few nanoseconds.