Z-DNA duplexes are a particularly complicated test case for current force fields. We performed a set of explicit solvent molecular dynamics (MD) simulations with various AMBER force field parametrizations including our recent refinements of the ε/ζ and glycosidic torsions. None of these force fields described the ZI/ZII and other backbone substates correctly, and all of them underpredicted the population of the important ZI substate. We show that this underprediction can be attributed to an inaccurate potential for the sugar-phosphate backbone torsion angle β. We suggest a refinement of this potential, β(OL1), which was derived using our recently introduced methodology that includes conformation-dependent solvation effects. The new potential significantly increases the stability of the dominant ZI backbone substate and improves the overall description of the Z-DNA backbone. It also has a positive (albeit small) impact on another important DNA form, the antiparallel guanine quadruplex (G-DNA), and improves the description of the canonical B-DNA backbone by increasing the population of BII backbone substates, providing a better agreement with experiment. We recommend using β(OL1) in combination with our previously introduced corrections, εζ(OL1) and χ(OL4), (the combination being named OL15) as a possible alternative to the current β torsion potential for more accurate modeling of nucleic acids.

The utility of molecular dynamics (MD) simulations to model biomolecular structure, dynamics, and interactions has witnessed enormous advances in recent years due to the availability of optimized MD software and access to significant computational power, including GPU multicore computing engines and other specialized hardware. This has led researchers to routinely extend conformational sampling times to the microsecond level and beyond. The extended sampling time has allowed the community not only to converge conformational ensembles through complete sampling but also to discover deficiencies and overcome problems with the force fields. Accuracy of the force fields is a key component, along with sampling, toward being able to generate accurate and stable structures of biopolymers. The Amber force field for nucleic acids has been used extensively since the 1990s, and multiple artifacts have been discovered, corrected, and reassessed by different research groups. We present a direct comparison of two of the most recent and state-of-the-art Amber force field modifications, bsc1 and OL15, that focus on accurate modeling of double-stranded DNA. After extensive MD simulations with five test cases and two different water models, we conclude that both modifications are a remarkable improvement over the previous bsc0 force field. Both force field modifications show better agreement when compared to experimental structures. To ensure convergence, the Drew–Dickerson dodecamer (DDD) system was simulated using 100 independent MD simulations, each extended to at least 10 μs, and the independent MD simulations were concatenated into a single 1 ms long trajectory for each combination of force field and water model. This is significantly beyond the time scale needed to converge the conformational ensemble of the internal portions of a DNA helix absent internal base pair opening. Considering all of the simulations discussed in the current work, the MD simulations performed to assess and validate the current force fields and water models aggregate over 14 ms of simulation time. The results suggest that both the bsc1 and OL15 force fields render average structures that deviate significantly less than 1 Å from the average experimental structures. This can be compared to similar but less exhaustive simulations with the CHARMM 36 force field that aggregate to the ∼90 μs time scale and also perform well but do not produce structures as close to the DDD NMR average structures (with root-mean-square deviations of 1.3 Å) as the newer Amber force fields. On the basis of these analyses, any future research involving double-stranded DNA simulations using the Amber force fields should employ the bsc1 or OL15 modification.

Background The structure and dynamics of DNA are critically related to its function. Molecular dynamics (MD) simulations augment experiment by providing detailed information about the atomic motions. However, to date the simulations have not been long enough for convergence of the dynamics and structural properties of DNA. Methods MD simulations performed with AMBER using the ff99SB force field with the parmbsc0 modifications, including ensembles of independent simulations, were compared to long timescale MD performed with the specialized Anton MD engine on the B-DNA structure d(GCACGAACGAACGAACGC). To assess convergence, the decay of the average RMSD values over longer and longer time intervals was evaluated in addition to assessing convergence of the dynamics via the Kullback-Leibler divergence of principal component projection histograms. Results These MD simulations —including one of the longest simulations of DNA published to date at ~44 μs—surprisingly suggest that the structure and dynamics of the DNA helix, neglecting the terminal base pairs, are essentially fully converged on the ~1–5 μs timescale. Conclusions We can now reproducibly converge the structure and dynamics of B-DNA helices, omitting the terminal base pairs, on the μs time scale with both the AMBER and CHARMM C36 nucleic acid force fields. Results from independent ensembles of simulations starting from different initial conditions, when aggregated, match the results from long timescale simulations on the specialized Anton MD engine. General Significance With access to large-scale GPU resources or the specialized MD engine “Anton” it is possibly for a variety of molecular systems to reproducibly and reliably converge the conformational ensemble of sampled structures.

DNA helices display a rich tapestry of motion on both short (< 100 ns) and long (> 1 ms) timescales. However, with the exception of mismatched or damaged DNA, experimental measures indicate that motions in the 1 µs to 1 ms range are effectively absent, which is often attributed to difficulties in measuring motions in this time range. We hypothesized that these motions have not been measured because there is effectively no motion on this timescale, as this provides a means to distinguish faithful Watson-Crick base paired DNA from damaged DNA. The absence of motion on this timescale would present a “static” DNA sequence-specific structure that matches the encounter timescales of proteins, thereby facilitating recognition. Here we report long timescale (~10-44 µs) molecular dynamics simulations of a B-DNA duplex structure that addresses this hypothesis using both an “Anton” machine and large ensembles of AMBER GPU simulations.

The family of anticancer complexes that include the transition metal copper known as Casiopeínas® shows promising results. Two of these complexes are currently in clinical trials. The interaction of these compounds with DNA has been observed experimentally and several hypotheses regarding the mechanism of action have been developed, and these include the generation of reactive oxygen species, phosphate hydrolysis and/or base-pair intercalation. To advance in the understanding on how these ligands interact with DNA, we present a molecular dynamics study of 21 Casiopeínas with a DNA dodecamer using 10 μs of simulation time for each compound. All the complexes were manually inserted into the minor groove as the starting point of the simulations. The binding energy of each complex and the observed representative type of interaction between the ligand and the DNA is reported. With this extended sampling time, we found that four of the compounds spontaneously flipped open a base pair and moved inside the resulting cavity and four compounds formed stacking interactions with the terminal base pairs. The complexes that formed the intercalation pocket led to more stable interactions.

The nuclear receptor ligand-binding domain (LBD) is a highly dynamic entity. Crystal structures have defined multiple low-energy LBD structural conformations of the activation function-2 (AF-2) co-regulator-binding surface, yet it remains unclear how ligand binding influences the number and population of conformations within the AF-2 structural ensemble. Here, we present a nuclear receptor co-regulator-binding surface structural ensemble in solution, viewed through the lens of fluorine-19 (19F) nuclear magnetic resonance (NMR) and molecular simulations, and the response of this ensemble to ligands, co-regulator peptides and heterodimerization. We correlate the composition of this ensemble with function in peroxisome proliferator-activated receptor-γ (PPARγ) utilizing ligands of diverse efficacy in co-regulator recruitment. While the co-regulator surface of apo PPARγ and partial-agonist-bound PPARγ is characterized by multiple thermodynamically accessible conformations, the full and inverse-agonist-bound PPARγ co-regulator surface is restricted to a few conformations which favor coactivator or corepressor binding, respectively.

Visualization of double stranded DNA in gels with the binding of the fluorescent dye ethidium bromide has been a basic experimental technique in any molecular biology laboratory for &gt;40 years. The interaction between ethidium and double stranded DNA has been observed to be an intercalation between base pairs with strong experimental evidence. This presents a unique opportunity for computational chemistry and biomolecular simulation techniques to benchmark and assess their models in order to see if the theory can reproduce experiments and ultimately provide new insights. We present molecular dynamics simulations of the interaction of ethidium with two different double stranded DNA models. The first model system is the classic sequence d(CGCGAATTCGCG)2 also known as the Drew–Dickerson dodecamer. We found that the ethidium ligand binds mainly stacked on, or intercalated between, the terminal base pairs of the DNA with little to no interaction with the inner base pairs. As the intercalation at the terminal CpG steps is relatively rapid, the resultant DNA unwinding, rigidification, and increased stability of the internal base pair steps inhibits further intercalation. In order to reduce these interactions and to provide a larger groove space, a second 18-mer DNA duplex system with the sequence d(GCATGAACGAACGAACGC) was tested. We computed molecular dynamics simulations for 20 independent replicas with this sequence, each with ∼27 μs of sampling time. Results show several spontaneous intercalation and base-pair eversion events that are consistent with experimental observations. The present work suggests that extended MD simulations with modern DNA force fields and optimized simulation codes are allowing the ability to reproduce unbiased intercalation events that we were not able to previously reach due to limits in computing power and the lack of extensively tested force fields and analysis tools.

The aim of this work is to describe the specific recognition site between DNA and an anticancer copper complex by means of computational methods. Molecular dynamics were used to find the preferred site of binding between selected DNA chains and [Cu(2,2'-bipyridine)(acetylacetonate)(H(2)O)](+) (Cas). Full DFT optimizations of selected geometries extracted from simulations, followed by a topological analysis of electron density allowed us to define the specific interactions inside the recognition site. Cas links deoxyribose-phosphate by a coordination bond between metal and the phosphate group. There are several C-H···π, O···π(C) and O···π(N) interactions between the sugar group and the aromatic ligand of Cas. The results indicate that the adduct Cas-deoxyribose-phosphate may be an initial step toward the hydrolysis of DNA chains. Overall, this study provides insights into the initial step of the action mechanism of copper complexes as apoptosis-inducing agents and provides guidelines for the design of this kind of drugs.