Search

Modified nucleic acids have surged as a popular therapeutic route, emphasizing the importance of nucleic acid research in drug discovery and development. Beyond well-known RNA vaccines, antisense oligonucleotides and aptamers can incorporate various modified nucleic acids to target specific biomolecules for various therapeutic activities. Molecular dynamics simulations can accelerate the design and development of these systems with noncanonical nucleic acids by observing intricate dynamic properties and relative stability on the all-atom level. However, modeling these modified systems is challenging due to the time and resources required to parametrize components outside default force field parameters. Here, we present modXNA, a tool to derive and build modified nucleotides for use with Amber force fields. Several nucleic acid systems varying in size and number of modification sites were used to evaluate the accuracy of modXNA parameters, and results indicate the dynamics and structure are preserved throughout the simulations. We detail the protocol for quantum mechanics charge derivation and describe a workflow for implementing modXNA in Amber molecular dynamics simulations, which includes updates and added features to CPPTRAJ.

We develop a multiscale coarse-grain model of the NIST Monoclonal Antibody Reference Material 8671 (NISTmAb) to enable systematic computational investigations of high-concentration physical instabilities such as phase separation, clustering, and aggregation. Our multiscale coarse-graining strategy captures atomic-resolution interactions with a computational approach that is orders of magnitude more efficient than atomistic models, assuming the biomolecule can be decomposed into one or more rigid bodies with known, fixed structures. This method reduces interactions between tens of thousands of atoms to a single anisotropic interaction site. The anisotropic interaction between unique pairs of rigid bodies is precomputed over a discrete set of relative orientations and stored, allowing interactions between arbitrarily oriented rigid bodies to be interpolated from the precomputed table during coarse-grained Monte Carlo simulations. We present this approach for lysozyme and lactoferrin as a single rigid body and for the NISTmAb as three rigid bodies bound by a flexible hinge with an implicit solvent model. This coarse-graining strategy predicts experimentally measured radius of gyration and second osmotic virial coefficient data, enabling routine Monte Carlo simulation of medically relevant concentrations of interacting proteins while retaining atomistic detail. All methodologies used in this work are available in the open-source software Free Energy and Advanced Sampling Simulation Toolkit., (© 2024 Author(s).)

The crystallizable fragment (Fc) domain of immunoglobulin subclass IgG1 antibodies is engineered for a wide variety of pharmaceutical applications. Two important structural variables in Fc constructs are the hinge region connecting the Fc to the antigen binding fragments (Fab) and the glycans present in various glycoforms. These components affect receptor binding interactions that mediate immune activation. To design new antibody drugs, a robust in silico method for linking stability to structural changes is necessary. In this work, all-atom simulations were used to compare the dynamic behavior of the four structural variants arising from presence or absence of the hinge and glycans. We expressed the simplest of these constructs, the 'minimal Fc' with no hinge and no glycans, in Escherichia coli and report its crystal structure. The 'maximal Fc' that includes full hinge and G0F/G1F glycans is based on a previously reported structure, Protein Data Bank (PDB) ID: 5VGP. These, along with two intermediate structures (with only the glycans or with only the hinge) were used to independently measure the stability effects of the two structural variables using umbrella sampling simulations. Principal component analysis (PCA) was used to determine free energy effects along the Fc's dominant mode of motion. This work provides a comprehensive picture of the effects of hinge and glycans on Fc dynamics and stability.Communicated by Ramaswamy H. Sarma.

The dosing and efficacy of chemotherapeutic drugs can be limited by toxicity caused by off-pathway reactions. One hypothesis for how such toxicity arises is via metal-catalyzed oxidative damage of cardiac myosin binding protein C (cMyBP-C) found in cardiac tissue. Previous research indicates that metal ion mediated reactive oxygen species induce high levels of protein carbonylation, changing the structure and function of this protein. In this work, we use long timescale all-atom molecular dynamics simulations to investigate the ion environment surrounding the C0 and C1 subunits of cMyBP-C responsible for actin binding. We show that divalent cations are co-localized with protein carbonylation-prone amino acid residues and that carbonylation of these residues can lead to site-specific interruption to the actin-cMyBP-C binding., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Published by Elsevier Inc.)

Numerous technologies are currently in development for use in next-generation protein sequencing platforms. A notable published approach employs fluorescently-tagged binding proteins to identity the N-terminus of immobilized peptides, in-between rounds of digestion. This approach makes use of N-terminal amino acid binder (NAAB) proteins, which would identify amino acids by chemical and shape complementarity. One source of NAABs is the ClpS protein family, which serve to recruit proteins to bacterial proteosomes based on the identity of the N-terminal amino acid. In this study, a Thermosynechococcus vestitus (also known as Thermosynechococcus elongatus) ClpS2 protein was used as the starting point for direct evolution of an NAAB with affinity and specificity for N-terminal leucine. Enriched variants were analyzed and shown to improve the interaction between the ClpS surface and the peptide chain, without increasing promiscuity. Interestingly, interactions were found that were unanticipated which favor different charged residues located at position 5 from the N-terminus of a target peptide., (© The Author(s) 2023. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Human hepatitis B virus (HBV) replication is initiated by the binding of the viral polymerase (P) to epsilon (ε), an ≈85-nucleotide (nt) cis-acting regulatory stem-loop RNA located at the 5'-end of the pre-genomic RNA (pgRNA). This interaction triggers P and pgRNA packaging and protein-primed reverse transcription and is therefore an attractive therapeutic target. Our recent nuclear magnetic resonance (NMR) structure of ε provides a useful starting point toward a detailed understanding of HBV replication, and hints at the functional importance of ε dynamics. Here, we present a detailed description of ε motions on the ps to ns and μs to ms time scales by NMR spin relaxation and relaxation dispersion, respectively. We also carried out molecular dynamics simulations to provide additional insight into ε conformational dynamics. These data outline a series of complex motions on multiple time scales within ε. Moreover, these motions occur in mostly conserved nucleotides from structural regions (i.e., priming loop, pseudo-triloop, and U43 bulge) that biochemical and mutational studies have shown to be essential for P binding, P-pgRNA packaging, protein-priming, and DNA synthesis. Taken together, our work implicates RNA dynamics as an integral feature that governs HBV replication., Competing Interests: Declaration of interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. Certain commercial equipment, instruments, and materials are identified in this article in order to specify the experimental procedure. Such identification does not imply recommendation or endorsement by the National Institute of Standards and Technology, nor does it imply that the material or equipment identified is necessarily the best available for the purpose., (Copyright © 2022. Published by Elsevier Ltd.)

We describe the conformational ensemble of the single-stranded r(UCAAUC) oligonucleotide obtained using extensive molecular dynamics (MD) simulations and Rosetta's FARFAR2 algorithm. The conformations observed in MD consist of A-form-like structures and variations thereof. These structures are not present in the pool generated using FARFAR2. By comparing with available nuclear magnetic resonance (NMR) measurements, we show that the presence of both A-form-like and other extended conformations is necessary to quantitatively explain experimental data. To further validate our results, we measure solution X-ray scattering (SAXS) data on the RNA hexamer and find that simulations result in more compact structures than observed from these experiments. The integration of simulations with NMR via a maximum entropy approach shows that small modifications to the MD ensemble lead to an improved description of the conformational ensemble. Nevertheless, we identify persisting discrepancies in matching experimental SAXS data., (© 2022 Bergonzo et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

Setting up molecular dynamics simulations from experimentally determined structures is often complicated by a variety of factors, particularly the inclusion of carbohydrates, since these have several anomer types which can be linked in a variety of ways. Here we present a stand-alone tool implemented in the widely-used software CPPTRAJ that can be used to automate building structures and generating a "ready to run" parameter and coordinate file pair. This tool automatically identifies carbohydrate anomer type, configuration, linkage, and functional groups, and performs topology modifications (e.g., renaming residue/atom names) required to build the final system using state of the art GLYCAM force field parameters. It will also generate the necessary commands for bonding carbohydrates and creating any disulfide bonds., (Published 2022. This article is a U.S. Government work and is in the public domain in the USA.)

Previous reports present different models for the stabilization of the Fc-FcγRI immune complex. Although accord exists on the importance of L235 in IgG1 and some hydrophobic contacts for complex stabilization, discord exists regarding the existence of stabilizing glycoprotein contacts between glycans of IgG1 and a conserved FG-loop ( 171 MGKHRY 176 ) of FcγRIa. Complexes formed from the FcγRIa receptor and IgG1s containing biantennary glycans with N-acetylglucosamine, galactose, and α2,6-N-acetylneuraminic terminations were measured by hydrogen-deuterium exchange mass spectrometry (HDX-MS), classified for dissimilarity with Welch's ANOVA and Games-Howell post hoc procedures, and modeled with molecular dynamics (MD) simulations. For each glycoform of the IgG1-FcγRIa complex peptic peptides of Fab, Fc and FcγRIa report distinct H/D exchange rates. MD simulations corroborate the differences in the peptide deuterium content through calculation of the percent of time that transient glycan-peptide bonds exist. These results indicate that stability of IgG1-FcγRIa complexes correlate with the presence of intermolecular glycoprotein interactions between the IgG1 glycans and the 173 KHR 175 motif within the FG-loop of FcγRIa. The results also indicate that intramolecular glycan-protein bonds stabilize the Fc region in isolated and complexed IgG1. Moreover, HDX-MS data evince that the Fab domain has glycan-protein binding contacts within the IgG1-FcγRI complex., Competing Interests: Declaration of interests The authors declare that they have no competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Published by Elsevier Ltd.)

ClpS2 is a small protein under development as a probe for selectively recognizing N-terminal amino acids of N-degron peptide fragments. To understand the structural basis of ClpS2 specificity for an N-terminal amino acid, all atom molecular dynamics (MD) simulations were conducted using the sequence of a bench-stable mutant of ClpS2, called PROSS. We predicted that a single amino acid leucine to asparagine substitution would switch the specificity of PROSS ClpS2 to an N-terminal tyrosine over the preferred phenylalanine. Experimental validation of the mutant using a fluorescent yeast-display assay showed an increase in tyrosine binding over phenylalanine, in support of the proposed hypothesis., (Published 2020. This article is a U.S. Government work and is in the public domain in the USA.)

A considerable amount of rapid-paced research is underway to combat the SARS-CoV-2 pandemic. In this work, we assess the 3D structure of the 5' untranslated region of its RNA, in the hopes that stable secondary structures can be targeted, interrupted, or otherwise measured. To this end, we have combined molecular dynamics simulations with previous Nuclear Magnetic Resonance measurements for stem loop 2 of SARS-CoV-1 to refine 3D structure predictions of that stem loop. We find that relatively short sampling times allow for loop rearrangement from predicted structures determined in absence of water or ions, to structures better aligned with experimental data. We then use molecular dynamics to predict the refined structure of the transcription regulatory leader sequence (TRS-L) region which includes stem loop 3, and show that arrangement of the loop around exchangeable monovalent potassium can interpret the conformational equilibrium determined by in-cell dimethyl sulfate (DMS) data.

Monoclonal antibodies (mAbs) represent an important platform for the development of biotherapeutic products. Most mAbs are produced in mammalian cells, but several mAbs are made in Escherichia coli , including therapeutic fragments. The NISTmAb is a well-characterized reference material made widely available to facilitate the development of both originator biologics and biosimilars. Here, when expressing NISTmAb from codon-optimized constructs in E. coli (eNISTmAb), a truncated variant of its heavy chain was observed. N-terminal protein sequencing and mutagenesis analyses indicated that the truncation resulted from an internal translation initiation from a GTG codon (encoding Val) within eNISTmAb. Using computational and biochemical approaches, we demonstrate that this translation initiates from a weak Shine-Dalgarno sequence and is facilitated by a putative ribosomal protein S1-binding site. We also observed similar internal initiation in the mAb adalimumab (the amino acid sequence of the drug Humira) when expressed in E. coli Of note, these internal initiation regions were likely an unintended result of the codon optimization for E. coli expression, and the amino acid pattern from which it is derived was identified as a Pro-Ser- X-X-X -Val motif. We discuss the implications of our findings for E. coli protein expression and codon optimization and outline possible strategies for reducing the likelihood of internal translation initiation and truncated product formation.

Molecular dynamics (MD) simulations play an important role in characterizing Ribonucleic Acid (RNA) structure, augmenting information from experimental techniques such as Nuclear Magnetic Resonance (NMR). In this work, we examine the accuracy of structural representation resulting from application of a number of explicit and implicit solvent models and refinement protocols against experimental data ranging from high density of residual dipolar coupling (RDC) restraints to completely unrestrained simulations. For a prototype A-form RNA helix, our results indicate that AMBER RNA force field with either implicit or explicit solvent can produce a realistic dynamic representation of RNA helical structure, accurately cross-validating with respect to a diverse array of NMR observables. In refinement against NMR distance restraints, modern MD force fields are found to be equally adequate, with high fidelity cross-validation to the residual dipolar couplings (RDCs) and residual chemical shift anisotropies (RCSAs), while slightly over-estimating structural order as monitored via NMR relaxation data. With restraints trimmed to encode only for base pairing information, cross-validation quality significantly deteriorates, now exhibiting a pronounced dependence on the choice of the solvent model. This deterioration is found to be partially reversible by increasing planarity restraints on the nucleobase geometry. For completely unrestrained MD simulations, the choice of water model becomes very important, with the best-performing TIP4P-Ew accurately reproducing both the RDC and RCSA data, while closely matching the NMR-derived order parameters. The information provided here will serve as a foundation for MD-based refinement of solution state NMR structures of RNA., (Published by Elsevier Inc.)

Despite their great promise as fluorescent biological probes and sensors, the structure and dynamics of Ag complexes derived from single stranded DNA (ssDNA) are less understood than their double stranded counterparts. In this work, we seek new insights into the structure of single AgNssDNA clusters using analytical ultracentrifugation (AUC), nuclear magnetic resonance spectroscopy, infrared spectroscopy and molecular dynamics simulations (MD) of a fluorescent (AgNssDNA)8+ nanocluster. The results suggest that the purified (AgNssDNA)8+ nanocluster is a mixture of predominantly Ag15 and Ag16 species that prefer two distinct long-lived conformational states: one extended, the other approaching spherical. However, the ssDNA strands within these clusters are highly mobile. Ag(i) interacts preferentially with the nucleobase rather than the phosphate backbone, causing a restructuring of the DNA strand relative to the bare DNA. Infrared spectroscopy and MD simulations of (AgNssDNA)8+ and model nucleic acid homopolymers suggest that Ag(i) has a higher affinity for cytosine over guanine bases, little interaction with adenine, and virtually none with thymine. Ag(i) shows a tendency to interact with cytosine N3 and O2 and guanine N7 and O6, opening the possibility for a Ag(i)-base bifurcated bond to act as a nanocluster nucleation and strand stabilizing site. This work provides valuable insight into nanocluster structure and dynamics which drive stability and optical properties, and additional studies using these types of characterization techniques are important for the rational design of single stranded AgDNA nanocluster sensors.

Structural information about ribonucleic acid (RNA) is lagging behind that of proteins, in part due to its high charge and conformational variability. Molecular dynamics (MD) has played an important role in describing RNA structure, complementing information from both nuclear magnetic resonance (NMR), or X-ray crystallography. We examine the impact of the choice of the empirical force field for RNA structure refinement using cross-validation against residual dipolar couplings (RDCs) as structural accuracy reporter. Four force fields, representing both the state-of-the art in RNA simulation and the most popular selections in NMR structure determination, are compared for a prototypical A-RNA helix. RNA structural accuracy is also evaluated as a function of both density and nature of input NMR data including RDCs, anisotropic chemical shifts, and distance restraints. Our results show a complex interplay between the experimental restraints and the force fields indicating two best-performing choices: high-fidelity refinement in explicit solvent, and the conformational database-derived potentials. Accuracy of RNA models closely tracks the density of 1-bond C-H RDCs, with other data types having beneficial, but smaller effects. At lower RDC density, or when refining against NOEs only, the two selected force fields are capable of accurately describing RNA helices with little or no experimental RDC data, making them available for the higher order structure assembly or better quantification of the intramolecular dynamics. Unrestrained simulations of simple RNA motifs with state-of-the art MD force fields appear to capture the flexibility inherent in nucleic acids while also maintaining a good agreement with the experimental observables.

The interactions in the tertiary structure of a ribosomal RNA fragment in the GTPase Associating Center (GAC) have been experimentally studied, but the roles of the bound and diffuse cations in its folding pathway have not yet been fully elucidated. Melting experiments have shown that the temperature of the first of the two distinguishable transitions in the unfolding pathway of the GAC RNA can be regulated by altering the magnesium concentration, yet the physical interpretation of such ion-dependent effects on folding have not been clearly understood in spite of the availability of crystal structures that depict many GAC RNA-ion interactions. Here, we use umbrella sampling and molecular dynamics (MD) simulations to provide a physical description for the first transition in this unfolding pathway, with a focus on the role of a chelated magnesium ion. Our results indicate that the presence of cations mediating the local interaction of two loops stabilizes the folded state relative to the unfolded or partially folded states. Also, our findings suggest that a bridging magnesium ion between the two loops improves the stabilizing effect. This is consistent with the multistep unfolding pathway proposed for the GAC RNA and highlights the importance of ions in the first unfolding step. The results suggest how MD simulations can provide insight into RNA unfolding pathways as a complementary approach to experiments.

Though the structure of the substrate stem loop I (SLI)-stem loop V (SLV) kissing loop junction of the Varkud Satellite ribozyme has been experimentally characterized, the dynamics of this Mg 2+ -dependent loop-loop interaction have been elusive. Specifically, each hairpin loop contains a U-turn motif, but only SLV shows a conformational shift triggered by Mg 2+ ion association. Here, we use molecular dynamics simulations to analyze the binding and dynamics of this kissing loop junction. We show that SLV acts as a scaffold, providing stability to the junction. Mg 2+ ions associate with SLV when it is part of the junction in a manner similar to when it is unbound, but there is no specificity in Mg 2+ binding for the SLI loop. This suggests that the entropic penalty of ordering the larger SLI is too high, allowing SLV to act as a scaffold for multiple substrate loop sequences., (Copyright © 2017 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

8-Oxoguanine (8-oxoG), a mutagenic DNA lesion generated under oxidative stress, differs from its precursor guanine by only two substitutions (O 8 and H 7 ). Human 8-oxoguanine glycosylase 1 (OGG1) can locate and remove 8-oxoG through extrusion and excision. To date, it remains unclear how OGG1 efficiently distinguishes 8-oxoG from a large excess of undamaged DNA bases. We recently showed that formamidopyrimidine-DNA glycosylase (Fpg), a bacterial functional analog of OGG1, can selectively facilitate eversion of oxoG by stabilizing several intermediate states, and it is intriguing whether OGG1 also employs a similar mechanism in lesion recognition. Here, we use molecular dynamics simulations to explore the mechanism by which OGG1 discriminates between 8-oxoG and guanine along the base-eversion pathway. The MD results suggest an important role for kinking of the DNA by the glycosylase, which positions DNA phosphates in a way that assists lesion recognition during base eversion. The computational predictions were validated through experimental enzyme assays on phosphorothioate substrate analogs. Our simulations suggest that OGG1 distinguishes between 8-oxoG and G using their chemical dissimilarities not only at the active site but also at earlier stages during base eversion, and this mechanism is at least partially conserved in Fpg despite a lack of structural homology. The similarity also suggests that lesion recognition through multiple gating steps may be a common theme in DNA repair. Our results provide new insight into how enzymes can exploit kinetics and DNA conformational changes to probe the chemical modifications present in DNA lesions.

We compare the performance of five magnesium (Mg(2+)) 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 Mg(2+) 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-Mg(2+) interactions via first shell water molecules are surprisingly well described by modern parameters, allowing us to observe the spontaneous conformational shift from Mg(2+) free RNA to Mg(2+) associated RNA structure in unrestrained molecular dynamics simulations.

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., (© The Author(s) 2015. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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 Mg(2+)-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 Mg(2+) ions in stabilizing the loop, as well as the role of the solvent exposed U700 base. Additionally, these simulations suggest the Mg(2+)-free stem-loop adopts a wide range of structures, including energetically favorable structures similar to the Mg(2+)-bound loop structure. We propose this structure is a "gatekeeper" or precursor to Mg(2+) binding when those ions are present.

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.

Recent modifications and improvements to standard nucleic acid force fields have attempted to fix problems and issues that have been observed as longer timescale simulations have become routine. Although previous work has shown the ability to fold the UUCG stem-loop structure, until now no group has attempted to quantify the performance of current force fields using highly converged structural populations of the tetraloop conformational ensemble. In this study, we report the use of multiple independent sets of multidimensional replica exchange molecular dynamics (M-REMD) simulations with different initial conditions to generate well-converged conformational ensembles for the tetranucleotides r(GACC) and r(CCCC), as well as the larger UUCG tetraloop motif. By generating what is to our knowledge the most complete RNA structure ensembles reported to date for these systems, we remove the coupling between force field errors and errors due to incomplete sampling, providing a comprehensive comparison between current top-performing MD force fields for RNA. Of the RNA force fields tested in this study, none demonstrate the ability to correctly identify the most thermodynamically stable structure for all three systems. We discuss the deficiencies present in each potential function and suggest areas where improvements can be made. The results imply that although "short" (nsec-μsec timescale) simulations may stay close to their respective experimental structures and may well reproduce experimental observables, inevitably the current force fields will populate alternative incorrect structures that are more stable than those observed via experiment., (© 2015 Bergonzo et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

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., (© The Author(s) 2014. Published by Oxford University Press on behalf of Nucleic Acids Research.)

This unit presents an overview of computer simulation techniques as applied to nucleic acid systems, ranging from simple in vacuo molecular modeling techniques to more complete all-atom molecular dynamics treatments that include an explicit representation of the environment. The third in a series of four units, this unit focuses on critical issues in solvation and the treatment of electrostatics. UNITS 7.5 & 7.8 introduced the modeling of nucleic acid structure at the molecular level. This included a discussion of how to generate an initial model, how to evaluate the utility or reliability of a given model, and ultimately how to manipulate this model to better understand its structure, dynamics, and interactions. Subject to an appropriate representation of the energy, such as a specifically parameterized empirical force field, the techniques of minimization and Monte Carlo simulation, as well as molecular dynamics (MD) methods, were introduced as a way of sampling conformational space for a better understanding of the relevance of a given model. This discussion highlighted the major limitations with modeling in general. When sampling conformational space effectively, difficult issues are encountered, such as multiple minima or conformational sampling problems, and accurately representing the underlying energy of interaction. In order to provide a realistic model of the underlying energetics for nucleic acids in their native environments, it is crucial to include some representation of solvation (by water) and also to properly treat the electrostatic interactions. These subjects are discussed in detail in this unit., (Copyright © 2014 John Wiley & Sons, Inc.)

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.

The last in a set of units by the same authors, this unit addresses some important remaining questions about molecular modeling of nucleic acids. The unit describes how to choose an appropriate molecular mechanics force field; how to set up and equilibrate the system for accurate simulation of a nucleic acid in an explicit solvent by molecular dynamics or Monte Carlo simulation; and how to analyze molecular dynamics trajectories., (Copyright © 2014 John Wiley & Sons, Inc.)

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.

An overview of computer simulation techniques as applied to nucleic acid systems is presented. This unit discusses methods used to treat energy and to sample representative configurations. Emphasis is placed on molecular mechanics and empirical force fields., (Copyright © 2013 John Wiley & Sons, Inc.)

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., (Copyright © 2013 John Wiley & Sons, Inc.)

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.

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.

The nudged elastic band method (NEB) can be used to find a minimum energy path between two given starting structures. This method has been available in the standard release of the Amber9 and Amber10 suite of programs. In this paper a novel implementation of this method will be discussed, in which the nudged elastic band method is applied to only a specific, user-defined subset of atoms in a particular system, returning comparable results and minimum energy pathways as the standard implementation for an alanine dipeptide test system. This allows incorporation of explicit solvent with simulated systems, which may be preferred in many cases to an implicit solvent model. From a computational standpoint, this implementation of NEB also reduces the communication overhead inside the code, resulting in better performance for larger systems.