Dataset Information

Large-Scale Analysis of 48 DNA and 48 RNA Tetranucleotides Studied by 1 ?s Explicit-Solvent Molecular Dynamics Simulations.

ABSTRACT: An understanding of how the conformational behavior of single-stranded DNAs and RNAs depend on sequence is likely to be important for attempts to rationalize the thermodynamics of nucleic acid folding. In an attempt to further our understanding of such sequence dependences, we report here the results of 192 (1 ?s) explicit-solvent molecular dynamics (MD) simulations of 48 DNA and 48 RNA tetranucleotide sequences performed using recently reported modifications to the AMBER force field. Each tetranucleotide was simulated starting from two different conformations, a fully natively stacked and a completely unstacked conformation, and populations of the various possible base stacking arrangements were analyzed. The simulations indicate that, for both DNA and RNA, the populations of fully natively stacked conformations increase linearly with increasing numbers of purines in the sequence, while the conformational entropies, computed by two complementary methods, decrease. Despite the comparatively short simulation times, the computed free energies of stacking of the 16 possible combinations of bases in the middle of the sequences are found to be in good correspondence with values reported recently from simulations of dinucleoside monophosphates using the same force field. Finally, consistent with recent reports from other groups, non-native stacking interactions, i.e., between bases that are not adjacent in sequence, are shown to be a recurring feature of the simulations; in particular, stacking interactions of bases in a i:i+2 relationship are shown to occur significantly more frequently when the intervening base is a pyrimidine. Given that the high prevalence of non-native stacking interactions is thought to be unrealistic, it appears that further parametrization work will be required before accurate conformational descriptions of single-stranded nucleic acids can be obtained with current force fields.

SUBMITTER: Schrodt MV

PROVIDER: S-EPMC4806854 | biostudies-literature | 2015 Dec

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Large-Scale Analysis of 48 DNA and 48 RNA Tetranucleotides Studied by 1 μs Explicit-Solvent Molecular Dynamics Simulations.

Schrodt Michael V MV Andrews Casey T CT Elcock Adrian H AH

Journal of chemical theory and computation 20151118 12

An understanding of how the conformational behavior of single-stranded DNAs and RNAs depend on sequence is likely to be important for attempts to rationalize the thermodynamics of nucleic acid folding. In an attempt to further our understanding of such sequence dependences, we report here the results of 192 (1 μs) explicit-solvent molecular dynamics (MD) simulations of 48 DNA and 48 RNA tetranucleotide sequences performed using recently reported modifications to the AMBER force field. Each tetra ...[more]

PMID: 26580891

Similar Datasets

Project description:Adequate sampling of conformation space remains challenging in atomistic simulations, especially if the solvent is treated explicitly. Implicit-solvent simulations can speed up conformational sampling significantly. We compare the speed of conformational sampling between two commonly used methods of each class: the explicit-solvent particle mesh Ewald (PME) with TIP3P water model and a popular generalized Born (GB) implicit-solvent model, as implemented in the AMBER package. We systematically investigate small (dihedral angle flips in a protein), large (nucleosome tail collapse and DNA unwrapping), and mixed (folding of a miniprotein) conformational changes, with nominal simulation times ranging from nanoseconds to microseconds depending on system size. The speedups in conformational sampling for GB relative to PME simulations, are highly system- and problem-dependent. Where the simulation temperatures for PME and GB are the same, the corresponding speedups are approximately onefold (small conformational changes), between ∼1- and ∼100-fold (large changes), and approximately sevenfold (mixed case). The effects of temperature on speedup and free-energy landscapes, which may differ substantially between the solvent models, are discussed in detail for the case of miniprotein folding. In addition to speeding up conformational sampling, due to algorithmic differences, the implicit solvent model can be computationally faster for small systems or slower for large systems, depending on the number of solute and solvent atoms. For the conformational changes considered here, the combined speedups are approximately twofold, ∼1- to 60-fold, and ∼50-fold, respectively, in the low solvent viscosity regime afforded by the implicit solvent. For all the systems studied, 1) conformational sampling speedup increases as Langevin collision frequency (effective viscosity) decreases; and 2) conformational sampling speedup is mainly due to reduction in solvent viscosity rather than possible differences in free-energy landscapes between the solvent models.

Project description:Although BRACO19 is a potent G-quadruplex binder, its potential for clinical usage is hindered by its low selectivity towards DNA G-quadruplex over duplex. High-resolution structures of BRACO19 in complex with neither single-stranded telomeric DNA G-quadruplexes nor B-DNA duplex are available. In this study, the binding pathway of BRACO19 was probed by 27.5 µs molecular dynamics binding simulations with a free ligand (BRACO19) to a DNA duplex and three different topological folds of the human telomeric DNA G-quadruplex (parallel, anti-parallel and hybrid). The most stable binding modes were identified as end stacking and groove binding for the DNA G-quadruplexes and duplex, respectively. Among the three G-quadruplex topologies, the MM-GBSA binding energy analysis suggested that BRACO19's binding to the parallel scaffold was most energetically favorable. The two lines of conflicting evidence plus our binding energy data suggest conformation-selection mechanism: the relative population shift of three scaffolds upon BRACO19 binding (i.e., an increase of population of parallel scaffold, a decrease of populations of antiparallel and/or hybrid scaffold). This hypothesis appears to be consistent with the fact that BRACO19 was specifically designed based on the structural requirements of the parallel scaffold and has since proven effective against a variety of cancer cell lines as well as toward a number of scaffolds. In addition, this binding mode is only slightly more favorable than BRACO19s binding to the duplex, explaining the low binding selectivity of BRACO19 to G-quadruplexes over duplex DNA. Our detailed analysis suggests that BRACO19's groove binding mode may not be stable enough to maintain a prolonged binding event and that the groove binding mode may function as an intermediate state preceding a more energetically favorable end stacking pose; base flipping played an important role in enhancing binding interactions, an integral feature of an induced fit binding mechanism.

Project description:Despite considerable advances in computing power, atomistic simulations under nonperiodic boundary conditions, with Coulombic electrostatic interactions and in systems large enough to reduce finite-size associated errors in thermodynamic quantities to within the thermal energy, are still not affordable. As a result, periodic boundary conditions, systems of microscopic size and effective electrostatic interaction functions are frequently resorted to. Ensuing artifacts in thermodynamic quantities are nowadays routinely corrected a posteriori, but the underlying configurational sampling still descends from spurious forces. The present study addresses this problem through the introduction of on-the-fly corrections to the physical forces during an atomistic molecular dynamics simulation. Two different approaches are suggested, where the force corrections are derived from special potential energy terms. In the first approach, the solvent-generated electrostatic potential sampled at a given atom site is restrained to a target value involving corrections for electrostatic artifacts. In the second approach, the long-range regime of the solvent polarization around a given atom site is restrained to the Born polarization, i.e., the solvent polarization corresponding to the ideal situation of a macroscopic system under nonperiodic boundary conditions and governed by Coulombic electrostatic interactions. The restraints are applied to the explicit-water simulation of a hydrated sodium ion, and the effect of the restraints on the structural and energetic properties of the solvent is illustrated. Furthermore, by means of the calculation of the charging free energy of a hydrated sodium ion, it is shown how the electrostatic potential restraint translates into the on-the-fly consideration of the corresponding free-energy correction terms. It is discussed how the restraints can be generalized to situations involving several solute particles. Although the present study considers a very simple system only, it is an important step toward the on-the-fly elimination of finite-size and approximate-electrostatic artifacts during atomistic molecular dynamics simulations.

Dataset Information

Large-Scale Analysis of 48 DNA and 48 RNA Tetranucleotides Studied by 1 ?s Explicit-Solvent Molecular Dynamics Simulations.

Publications

Large-Scale Analysis of 48 DNA and 48 RNA Tetranucleotides Studied by 1 μs Explicit-Solvent Molecular Dynamics Simulations.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets