Christina Bergonzo

Multidimensional Replica Exchange Molecular Dynamics Yields a Converged Ensemble of an RNA Tetranucleotide.

A necessary step to properly assess and validate the performance of force fields for biomolecules is to exhaustively sample the accessible conformational space, which is challenging for large RNA structures. Given questions regarding the reliability of modeling RNA structure and dynamics with current methods, we have begun to use RNA tetranucleotides to evaluate force fields. These systems, though small, display considerable conformational variability and complete sampling with standard simulation methods remains challenging. Here we compare and discuss the performance of known variations of replica exchange molecular dynamics (REMD) methods, specifically temperature REMD (T-REMD), Hamiltonian REMD (H-REMD), and multidimensional REMD (M-REMD) methods, which have been implemented in Amber’s accelerated GPU code. Using two independent simulations, we show that M-REMD not only makes very efficient use of emerging large-scale GPU clusters, like Blue Waters at the University of Illinois, but also is critically important in generating the converged ensemble more efficiently than either T-REMD or H-REMD. With 57.6 μs aggregate sampling of a conformational ensemble with M-REMD methods, the populations can be compared to NMR data to evaluate force field reliability and further understand how putative changes to the force field may alter populations to be in more consistent agreement with experiment.

Evaluation of enhanced sampling provided by accelerated molecular dynamics with Hamiltonian replica exchange methods.

Many problems studied via molecular dynamics require accurate estimates of various thermodynamic properties, such as the free energies of different states of a system, which in turn requires well-converged sampling of the ensemble of possible structures. Enhanced sampling techniques are often applied to provide faster convergence than is possible with traditional molecular dynamics simulations. Hamiltonian replica exchange molecular dynamics (H-REMD) is a particularly attractive method, as it allows the incorporation of a variety of enhanced sampling techniques through modifications to the various Hamiltonians. In this work, we study the enhanced sampling of the RNA tetranucleotide r(GACC) provided by H-REMD combined with accelerated molecular dynamics (aMD), where a boosting potential is applied to torsions, and compare this to the enhanced sampling provided by H-REMD in which torsion potential barrier heights are scaled down to lower force constants. We show that H-REMD and multidimensional REMD (M-REMD) combined with aMD does indeed enhance sampling for r(GACC), and that the addition of the temperature dimension in the M-REMD simulations is necessary to efficiently sample rare conformations. Interestingly, we find that the rate of convergence can be improved in a single H-REMD dimension by simply increasing the number of replicas from 8 to 24 without increasing the maximum level of bias. The results also indicate that factors beyond replica spacing, such as round trip times and time spent at each replica, must be considered in order to achieve optimal sampling efficiency.

An Improved Reaction Coordinate for Nucleic Acid Base Flipping Studies.

Base flipping is a common strategy utilized by many enzymes to gain access to the functional groups of nucleic acid bases in duplex DNA which are otherwise protected by the DNA backbone and hydrogen bonding with their partner bases. Several X-ray crystallography studies have revealed flipped conformations of nucleotides bound to enzymes. However, little is known about the base-flipping process itself, even less about the role of the enzymes. Computational studies have used umbrella sampling to elicit the free energy profile of the base-flipping process using a pseudodihedral angle to represent the reaction coordinate. In this study, we have used an unrestrained trajectory in which a flipped base spontaneously reinserted into the helix in order to evaluate and improve the previously defined pseudodihedral angle. Our modified pseudodihedral angles use a new atom selection to improve the numerical stability of the restraints and also provide better correlation to the extent of flipping observed in simulations. Furthermore, on the basis of the comparison of potential of mean force (PMF) generated using different reaction coordinates, we observed that the shape of a flipping PMF profile is strongly dependent on the definition of the reaction coordinate, even for the same data set.

Improved Force Field Parameters Lead to a Better Description of RNA Structure.

We compare the performance of two different RNA force fields in four water models in simulating the conformational ensembles r(GACC) and r(CCCC). With the increased sampling facilitated by multidimensional replica exchange molecular dynamics (M-REMD), populations are compared to NMR data to evaluate force field reliability. The combination of AMBER ff12 with vdW(bb) modifications and the OPC water model produces results in quantitative agreement with the NMR ensemble that have eluded us to date.

Active destabilization of base pairs by a DNA glycosylase wedge initiates damage recognition

Formamidopyrimidine-DNA glycosylase (Fpg) excises 8-oxoguanine (oxoG) from DNA but ignores normal guanine. We combined molecular dynamics simulation and stopped-flow kinetics with fluorescence detection to track the events in the recognition of oxoG by Fpg and its mutants with a key phenylalanine residue, which intercalates next to the damaged base, changed to either alanine (F110A) or fluorescent reporter tryptophan (F110W). Guanine was sampled by Fpg, as evident from the F110W stopped-flow traces, but less extensively than oxoG. The wedgeless F110A enzyme could bend DNA but failed to proceed further in oxoG recognition. Modeling of the base eversion with energy decomposition suggested that the wedge destabilizes the intrahelical base primarily through buckling both surrounding base pairs. Replacement of oxoG with abasic (AP) site rescued the activity, and calculations suggested that wedge insertion is not required for AP site destabilization and eversion. Our results suggest that Fpg, and possibly other DNA glycosylases, convert part of the binding energy into active destabilization of their substrates, using the energy differences between normal and damaged bases for fast substrate discrimination.

Energetic Preference of 8-oxoG Eversion Pathways in a DNA Glycosylase

Base eversion is a fundamental process in the biochemistry of nucleic acids, allowing proteins engaged in DNA repair and epigenetic modifications to access target bases in DNA. Crystal structures reveal end points of these processes, but not the pathways involved in the dynamic process of base recognition. To elucidate the pathway taken by 8-oxoguanine during base excision repair by Fpg, we calculated free energy surfaces during eversion of the damaged base through the major and minor grooves. The minor groove pathway and free energy barrier (6-7 kcal/mol) are consistent with previously reported results (Qi, Y.; Spong, M. C.; Nam, K.; Banerjee, A.; Jiralerspong, S.; Karplus, M.; Verdine, G. L. Nature 2009, 462, 762.) However, eversion of 8-oxoG through the major groove encounters a significantly lower barrier (3-4 kcal/mol) more consistent with experimentally determined rates of enzymatic sliding during lesion search (Blainey, P. C.; van Oijent, A. M.; Banerjee, A.; Verdine, G. L.; Xie, X. S. Proc. Natl. Acad. Sci. U.S.A. 2006, 103, 5752.). Major groove eversion has been suggested for other glycosylases, suggesting that in addition to function, dynamics of base eversion may also be conserved.

Divalent Ion Dependent Conformational Changes in an RNA Stem-Loop Observed by Molecular Dynamics

We compare the performance of five magnesium (Mg2+) ion models in simulations of an RNA stem loop which has an experimentally determined divalent ion dependent conformational shift. We show that despite their differences in parametrization and resulting van der Waals terms, including differences in the functional form of the nonbonded potential, when the RNA adopts its folded conformation, all models behave similarly across ten independent microsecond length simulations with each ion model. However, when the entire structure ensemble is accounted for, chelation of Mg2+ to RNA is seen in three of the five models, most egregiously and likely artifactual in simulations using a 12-6-4 model for the Lennard-Jones potential. Despite the simple nature of the fixed point-charge and van der Waals sphere models employed, and with the exception of the likely oversampled directed chelation of the 12-6-4 potential models, RNA–Mg2+ interactions via first shell water molecules are surprisingly well described by modern parameters, allowing us to observe the spontaneous conformational shift from Mg2+ free RNA to Mg2+ associated RNA structure in unrestrained molecular dynamics simulations.

Stem-Loop V of Varkud Satellite RNA Exhibits Characteristics of the Mg(2+) Bound Structure in the Presence of Monovalent Ions.

The Varkud Satellite RNA contains a self-cleaving ribozyme that has been shown to function independently of its surroundings. This 160 nucleotide ribozyme adopts a catalytically active tertiary structure that includes a kissing hairpin complex formed by stem-loop I and stem-loop V (SLV). The five-nucleotide 5′-rUGACU loop of the isolated SLV has been shown to adopt a Mg2+-dependent U-turn structure by solution NMR. This U-turn hairpin is examined here by molecular dynamics simulations in the presence of monovalent and divalent ions. Simulations confirm on an all-atom level the hypotheses for the role of the Mg2+ ions in stabilizing the loop, as well as the role of the solvent exposed U700 base. Additionally, these simulations suggest the Mg2+-free stem-loop adopts a wide range of structures, including energetically favorable structures similar to the Mg2+-bound loop structure. We propose this structure is a “gatekeeper” or precursor to Mg2+ binding when those ions are present.

Molecular modeling of nucleic acid structure.

This unit is the first in a series of four units covering the analysis of nucleic acid structure by molecular modeling. The unit provides an overview of the computer simulation of nucleic acids. Topics include the static structure model, computational graphics and energy models, the generation of an initial model, and characterization of the overall three‐dimensional structure. Curr. Protoc. Nucleic Acid Chem. 54:7.5.1‐7.5.13.

A dynamic checkpoint in oxidative lesion discrimination by formamidopyrimidine–DNA glycosylase

In contrast to proteins recognizing small-molecule ligands, DNA-dependent enzymes cannot rely solely on interactions in the substrate-binding centre to achieve their exquisite specificity. It is widely believed that substrate recognition by such enzymes involves a series of conformational changes in the enzyme–DNA complex with sequential gates favoring cognate DNA and rejecting nonsubstrates. However, direct evidence for such mechanism is limited to a few systems. We report that discrimination between the oxidative DNA lesion, 8-oxoguanine (oxoG) and its normal counterpart, guanine, by the repair enzyme, formamidopyrimidine-DNA glycosylase (Fpg), likely involves multiple gates. Fpg uses an aromatic wedge to open the Watson–Crick base pair and everts the lesion into its active site. We used molecular dynamics simulations to explore the eversion free energy landscapes of oxoG and G by Fpg, focusing on structural and energetic details of oxoG recognition. The resulting energy profiles, supported by biochemical analysis of site-directed mutants disturbing the interactions along the proposed path, show that Fpg selectively facilitates eversion of oxoG by stabilizing several intermediate states, helping the rapidly sliding enzyme avoid full extrusion of every encountered base for interrogation. Lesion recognition through multiple gating intermediates may be a common theme in DNA repair enzymes.