Force field validation for nucleic acid simulations: comparing energies and dynamics of a DNA dodecamer

Important questions exist regarding the quality of force fields used in molecular dynamics (MD) simulations and their interoperable use with other available MD implementations. NAMD is one of the most efficient and scalable parallel molecular dynamics codes for large-scale biomolecular simulations in the open source domain. It is the aim of this article to analyze and compare the dynamics of a benchmark DNA dodecamer d(CTTTTGCAAAAG)2 system, including its binding to a specific drug molecule arising from the use of various simulation protocols in NAMD using Amber98, with the dynamics arising from simulations of the same dodecamer using Amber98 in the AMBER package, one of the most well-established simulation codes for nucleic acids. Based upon a set of validation benchmarks, the details of which are discussed, we find that nucleic acid simulations using NAMD give meaningful results and that the essential features of the resulting dynamics are similar to those arising from the AMBER package. This sets the stage for reliable large-scale simulations of nucleic acids using NAMD.     

Spectroscopic detection of DNA quadruplexes by vibrational circular dichroism

The four-stranded G-quadruplex motif is a conformation frequently adopted by guanine-rich nucleic acids that plays an important role in biology, medicine, and nanotechnology. Although vibrational spectroscopy has been widely used to investigate nucleic acid structure, association of particular spectral features with the quadruplex structure has to date been ambiguous. In this work, experimental IR absorption and vibrational circular dichroism (VCD) spectra of the model quadruplex systems d(G)(8) and deoxyguanosine-5'-monophosphate (5'-dGMP) were analyzed using molecular dynamics (MD) and quantum-chemical modeling. The experimental spectra were unambiguously assigned to the quadruplex DNA arrangement, and several IR and VCD bands related to this structural motif were determined. Involvement of MD in the modeling was essential for realistic simulation of the spectra. The VCD signal was found to be more sensitive to dynamical structural variations than the IR signal. The combination of the spectroscopic techniques with multiscale simulations provides extended information about nucleic acid conformations and their dynamics.     

Detecting DNA mismatches with metallo-insertors: a molecular simulation study

Molecules that selectively recognize DNA mismatches (MMs) play a key role as nucleic acids probes and as chemotherapeutic agents. Metallo-insertors bind to the minor groove (mG) of double strand (ds) DNA, expelling the mismatched base pairs and acting as their π-stacking replacement. In contrast, metallo-intercalators bind to the major groove (MG) of ds DNA and π-stack to adjacent base pairs. In this study we focused on structural and energetic properties of Δ-[Rh(bpy)(2)(chrysi)](3+) (1), Δ-[Ru(bpy)(2)(ddpz)](2+) (2), and Δ-[Ru(bpy)(2)(eilatin)](2+) (3) as prototypical examples of metallo-insertors and intercalators. For all molecules we characterized both insertion and intercalation into a DNA dodecamer via force field based molecular dynamics (MD) and hybrid quantum-classical (QM/MM) MD simulations. A structural analysis of the 1-3/DNA noncovalent adducts reveals that the insertion provokes an untwist of the DNA, an opening of the mG and of the phosphate backbone in proximity of the mismatch, while the intercalation induces smaller changes of these structural parameters. This behavior appears to be correlated with the size of the inserting/intercalating ligand in proximity of the metal coordination site. Moreover, our simulations show that the different selectivity of 1 toward distinct MM types may be correlated with the thermodynamic stability of the MMs in the free DNA and with that of the corresponding insertion adduct. Understanding the factors which tune a specific insertion is of crucial importance for designing specific luminescent probes that selectively recognize MMs, as well as for developing more effective anticancer drugs active in MM repair of deficient cells lines.     

Framework Nucleic Acids for Cell Imaging and Therapy

Over the past decade,structural DNA nanotechnology has been well developed to be a promising and powerful technique to generate various nanostructures with programmability,spatial organization and biocompatibi-lity.With the advent of computer-aided tools,framework nucleic acids have been employed in a series of biomedical applications,ranging from biosensing,bioimaging,diagnosis,to therapeutics.In this review,we summarized recent advances in the construction of precisely assembled DNA nanostructures,and DNA-engineered biomimetics.We also outlined the challenges and opportunities for the translational applications of framework nucleic acids.     

Solvating,manipulating, damaging,and repairing DNA in a computer

This work highlights four different topics in modeling of DNA: (i) the importance of water and ions together with the structure and function of DNA; the hydration structure around the ions appears to be the determining factor in the ion coordination to DNA, as demonstrated in the results of our MD simulations; (ii) how MD simulations can be used to simulate single molecule manipulation experiments as a complement to reveal the structural dynamics of the studied biomolecules; (iii) how damaged DNA can be studied in computer simulations; and (iv) how repair of damaged DNA can be studied theoretically. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

An N log N approximation based on the natural organization of biomolecules for speeding up the computation of long range interactions

Presented here is a method, the hierarchical charge partitioning (HCP) approximation, for speeding up computation of pairwise electrostatic interactions in biomolecular systems. The approximation is based on multiple levels of natural partitioning of biomolecular structures into a hierarchical set of its constituent structural components. The charge distribution in each component is systematically approximated by a small number of point charges, which, for the highest level component, are much fewer than the number of atoms in the component. For short distances from the point of interest, the HCP uses the full set of atomic charges available. For long‐distance interactions, the approximate charge distributions with smaller sets of charges are used instead. For a structure consisting of N charges, the computational cost of computing the pairwise interactions via the HCP scales as O(N log N), under assumptions about the structural organization of biomolecular structures generally consistent with reality. A proof‐of‐concept implementation of the HCP shows that for large structures it can lead to speed‐up factors of up to several orders of magnitude relative to the exact pairwise O(N²) all‐atom computation used as a reference. For structures with more than 2000–3000 atoms the relative accuracy of the HCP (relative root‐mean‐square force error per atom), approaches the accuracy of the particle mesh Ewald (PME) method with parameter settings typical for biomolecular simulations. When averaged over a set of 600 representative biomolecular structures, the relative accuracies of the two methods are roughly equal. The HCP is also significantly more accurate than the spherical cutoff method. The HCP has been implemented in the freely available nucleic acids builder (NAB) molecular dynamics (MD) package in Amber tools. A 10 ns simulation of a small protein indicates that the HCP based MD simulation is stable, and that it can be faster than the spherical cutoff method. A critical benefit of the HCP approximation is that it is algorithmically very simple, and unlike the PME, the HCP is straightforward to use with implicit solvent models. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

RedMD—Reduced molecular dynamics package

We developed a software package (RedMD) to perform molecular dynamics simulations and normal mode analysis of reduced models of proteins, nucleic acids, and their complexes. With RedMD one can perform molecular dynamics simulations in a microcanonical ensemble, with Berendsen and Langevin thermostats, and with Brownian dynamics. We provide force field and topology generators which are based on the one‐bead per residue/nucleotide elastic network model and its extensions. The user can change the force field parameters with the command line options that are passed to generators. Also, the generators can be modified, for example, to add new potential energy functions. Normal mode analysis tool is available for elastic or anisotropic network models. The program is written in C and C++ languages and the structure/topology of a molecule is based on an XML format. OpenMP technology for shared‐memory architectures was used for code parallelization. The code is distributed under GNU public licence and available at http://bionano.icm.edu.pl/software/ . © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Gro2mat: A package to efficiently read gromacs output in MATLAB

Molecular dynamics (MD) simulations are a state‐of‐the‐art computational method used to investigate molecular interactions at atomic scale. Interaction processes out of experimental reach can be monitored using MD software, such as Gromacs. Here, we present the gro2mat package that allows fast and easy access to Gromacs output files from Matlab. Gro2mat enables direct parsing of the most common Gromacs output formats including the binary xtc‐format. No openly available Matlab parser currently exists for this format. The xtc reader is orders of magnitudes faster than other available pdb/ascii workarounds. Gro2mat is especially useful for scientists with an interest in quick prototyping of new mathematical and statistical approaches for Gromacs trajectory analyses. © 2014 Wiley Periodicals, Inc.     

Dendrimer building toolkit: Model building and characterization of various dendrimer architectures

We have developed a graphical user interface based dendrimer builder toolkit (DBT) which can be used to generate the dendrimer configuration of desired generation for various dendrimer architectures. The validation of structures generated by this tool was carried out by studying the structural properties of two well known classes of dendrimers: ethylenediamine cored poly(amidoamine) (PAMAM) dendrimer, diaminobutyl cored poly(propylene imine) (PPI) dendrimer. Using full atomistic molecular dynamics (MD) simulation we have calculated the radius of gyration, shape tensor and monomer density distribution for PAMAM and PPI dendrimer at neutral and high pH. A good agreement between the available simulation and experimental (small angle X‐ray and neutron scattering; SAXS, SANS) results and calculated radius of gyration was observed. With this validation we have used DBT to build another new class of nitrogen cored poly(propyl ether imine) dendrimer and study it's structural features using all atomistic MD simulation. DBT is a versatile tool and can be easily used to generate other dendrimer structures with different chemistry and topology. The use of general amber force field to describe the intra‐molecular interactions allows us to integrate this tool easily with the widely used molecular dynamics software AMBER. This makes our tool a very useful utility which can help to facilitate the study of dendrimer interaction with nucleic acids, protein and lipid bilayer for various biological applications. © 2012 Wiley Periodicals, Inc.     

Solid‐Phase Incorporation of an ATRP Initiator for Polymer–DNA Biohybrids

The combination of polymers with nucleic acids leads to materials with significantly advanced properties. To obviate the necessity and complexity of conjugating two macromolecules, a polymer initiator is described that can be directly covalently linked to DNA during solid‐phase synthesis. Polymer can then be grown from the DNA bound initiator, both in solution after the DNA‐initiator is released from the solid support as well as directly on the solid support, simplifying purification. The resulting polymer‐DNA hybrids were examined by chromatography and fluorescence methods that attested to the integrity of hybrids and the DNA. The ability to use DNA‐based supports expands the range of readily available molecules that can be used with the initiator, as exemplified by direct synthesis of a biotin polymer hybrid on solid‐support. This method expands the accessibility and range of advanced polymer biohybrid materials.     

Correlation of protein structure and dynamics to scalar couplings across hydrogen bonds

NMR-observable scalar couplings across hydrogen bonds in nucleic acids and proteins present a quantitative measure for the geometry and--by the implicit experimental time averaging--dynamics of hydrogen bonds. We have carried out in-depth molecular dynamics (MD) simulations with various force fields on three proteins: ubiquitin, the GB1 domain of protein G, and the SMN Tudor domain, for which experimental h3JNC' scalar couplings of backbone hydrogen bonds and various high-resolution X-ray structures are available. Theoretical average values for h3JNC' were calculated from the snapshots of these MD simulations either by density functional theory or by a geometric parametrization (Barfield, M. J. Am. Chem. Soc. 2002, 124, 4158-4168). No significant difference was found between the two methods. The results indicate that time-averaging using explicit water solvation in the MD simulations improves significantly the agreement between experimental and theoretical values for the lower resolution structures ubiquitin (1.8 A), Tudor domain (1.8 A), and protein G (2.1 A). Only marginal improvement is found for the high-resolution structure (1.1 A) of protein G. Hence, experimental h3JNC' values are compatible with a static, high-resolution structural model. The MD averaging of the low-resolution structures moves the averages of the rHO distance and the H...O=C angle theta closer to their respective values in the high-resolution structures, thereby improving the agreement using experimental h3JNC' data. In contrast, MD averaging with implicit water models deteriorates the agreement with experiment for all proteins. The differing behavior can be explained by an artifactual lengthening of H-bonds caused by the implicit water models.     

ALMOST: An all atom molecular simulation toolkit for protein structure determination

Almost ( al l atom mo lecular s imulation t oolkit) is an open source computational package for structure determination and analysis of complex molecular systems including proteins, and nucleic acids. Almost has been designed with two primary goals: to provide tools for molecular structure determination using various types of experimental measurements as conformational restraints, and to provide methods for the analysis and assessment of structural and dynamical properties of complex molecular systems. The methods incorporated in Almost include the determination of structural and dynamical features of proteins using distance restraints derived from nuclear Overhauser effect measurements, orientational restraints obtained from residual dipolar couplings and the structural restraints from chemical shifts. Here, we present the first public release of Almost , highlight the key aspects of its computational design and discuss the main features currently implemented. Almost is available for the most common Unix‐based operating systems, including Linux and Mac OS X. Almost is distributed free of charge under the GNU Public License, and is available both as a source code and as a binary executable from the project web site at http://www.open‐almost.org . Interested users can follow and contribute to the further development of Almost on http://sourceforge.net/projects/almost . © 2014 Wiley Periodicals, Inc.     

Evaluation of DNA Force Fields in Implicit Solvation

DNA structural deformations and dynamics are crucial to its interactions in the cell. Theoretical simulations are essential tools to explore the structure, dynamics, and thermodynamics of biomolecules in a systematic way. Molecular mechanics force fields for DNA have benefited from constant improvements during the last decades. Several studies have evaluated and compared available force fields when the solvent is modeled by explicit molecules. On the other hand, few systematic studies have assessed the quality of duplex DNA models when implicit solvation is employed. The interest of an implicit modeling of the solvent consists in the important gain in the simulation performance and conformational sampling speed. In this study, respective influences of the force field and the implicit solvation model choice on DNA simulation quality are evaluated. To this end, extensive implicit solvent duplex DNA simulations are performed, attempting to reach both conformational and sequence diversity convergence. Structural parameters are extracted from simulations and statistically compared to available experimental and explicit solvation simulation data. Our results quantitatively expose the respective strengths and weaknesses of the different DNA force fields and implicit solvation models studied. This work can lead to the suggestion of improvements to current DNA theoretical models.     

Molecular Dynamics Scoring of Protein–Peptide Models Derived from Coarse-Grained Docking

One of the major challenges in the computational prediction of protein–peptide complexes is the scoring of predicted models. Usually, it is very difficult to find the most accurate solutions out of the vast number of sometimes very different and potentially plausible predictions. In this work, we tested the protocol for Molecular Dynamics (MD)-based scoring of protein–peptide complex models obtained from coarse-grained (CG) docking simulations. In the first step of the scoring procedure, all models generated by CABS-dock were reconstructed starting from their original C-alpha trace representations to all-atom (AA) structures. The second step included geometry optimization of the reconstructed complexes followed by model scoring based on receptor–ligand interaction energy estimated from short MD simulations in explicit water. We used two well-known AA MD force fields, CHARMM and AMBER, and a CG MARTINI force field. Scoring results for 66 different protein–peptide complexes show that the proposed MD-based scoring approach can be used to identify protein–peptide models of high accuracy. The results also indicate that the scoring accuracy may be significantly affected by the quality of the reconstructed protein receptor structures.     

Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 1. Generalized Born

We present an implementation of generalized Born implicit solvent all-atom classical molecular dynamics (MD) within the AMBER program package that runs entirely on CUDA enabled NVIDIA graphics processing units (GPUs). We discuss the algorithms that are used to exploit the processing power of the GPUs and show the performance that can be achieved in comparison to simulations on conventional CPU clusters. The implementation supports three different precision models in which the contributions to the forces are calculated in single precision floating point arithmetic but accumulated in double precision (SPDP), or everything is computed in single precision (SPSP) or double precision (DPDP). In addition to performance, we have focused on understanding the implications of the different precision models on the outcome of implicit solvent MD simulations. We show results for a range of tests including the accuracy of single point force evaluations and energy conservation as well as structural properties pertainining to protein dynamics. The numerical noise due to rounding errors within the SPSP precision model is sufficiently large to lead to an accumulation of errors which can result in unphysical trajectories for long time scale simulations. We recommend the use of the mixed-precision SPDP model since the numerical results obtained are comparable with those of the full double precision DPDP model and the reference double precision CPU implementation but at significantly reduced computational cost. Our implementation provides performance for GB simulations on a single desktop that is on par with, and in some cases exceeds, that of traditional supercomputers.     

Implicit solvent simulations of DNA and DNA-protein complexes: agreement with explicit solvent vs experiment

Molecular dynamics simulations of biomolecules with implicit solvent reduce the computational cost and complexity of such simulations so that longer time scales and larger system sizes can be reached. While implicit solvent simulations of proteins have become well established, the success of implicit solvent in the simulation of nucleic acids has not been fully established to date. Results obtained in this study demonstrate that stable and efficient simulations of DNA and a protein-DNA complex can be achieved with an implicit solvent model based on continuum dielectric electrostatics. Differences in conformational sampling of DNA with two sets of atomic radii that are used to define the dielectric interface between the solute and the continuum dielectric model of the solvent are investigated. Results suggest that depending on the choice of atomic radii agreement is either closer to experimental data or to explicit solvent simulations. Furthermore, partial conformational transitions toward A-DNA conformations when salt is added within the implicit solvent framework are observed.     

Design of a colicin E7 based chimeric zinc-finger nuclease

Colicin E7 is a natural bacterial toxin. Its nuclease domain (NColE7) enters the target cell and kills it by digesting the nucleic acids. The HNH-motif as the catalytic centre of NColE7 at the C-terminus requires the positively charged N-terminal loop for the nuclease activity—offering opportunities for allosteric control in a NColE7-based artificial nuclease. Accordingly, four novel zinc finger nucleases were designed by computational methods exploiting the special structural features of NColE7. The constructed models were subjected to MD simulations. The comparison of structural stability and functional aspects showed that these models may function as safely controlled artificial nucleases. This study was complemented by random mutagenesis experiments identifying potentially important residues for NColE7 function outside the catalytic region.     

Comparison of DNA hydration patterns obtained using two distinct computational methods, molecular dynamics simulation and three-dimensional reference interaction site model theory

Because proteins and DNA interact with each other and with various small molecules in the presence of water molecules, we cannot ignore their hydration when discussing their structural and energetic properties. Although high-resolution crystal structure analyses have given us a view of tightly bound water molecules on their surface, the structural data are still insufficient to capture the detailed configurations of water molecules around the surface of these biomolecules. Thanks to the invention of various computational algorithms, computer simulations can now provide an atomic view of hydration. Here, we describe the apparent patterns of DNA hydration calculated by using two different computational methods: Molecular dynamics (MD) simulation and three-dimensional reference interaction site model (3D-RISM) theory. Both methods are promising for obtaining hydration properties, but until now there have been no thorough comparisons of the calculated three-dimensional distributions of hydrating water. This rigorous comparison showed that MD and 3D-RISM provide essentially similar hydration patterns when there is sufficient sampling time for MD and a sufficient number of conformations to describe molecular flexibility for 3D-RISM. This suggests that these two computational methods can be used to complement one another when evaluating the reliability of the calculated hydration patterns.