DNA Force Fields Research Articles

C-B-A Test of DNA Force Fields.

The DNA duplex may be locally strongly bent in complexes with proteins, for example, with polymerases or in a nucleosome. At such bends, the DNA helix is locally in the noncanonical forms A (with a narrow major groove and a large amount of north sugars) or C (with a narrow minor groove and a large share of BII phosphates). To model the formation of such complexes by molecular dynamics methods, the force field is required to reproduce these conformational transitions for a naked DNA. We analyzed the available experimental data on the B-C and B-A transitions under the conditions easily implemented in modeling: in an aqueous NaCl solution. We selected six DNA duplexes which conformations at different salt concentrations are known reliably enough. At low salt concentrations, poly(GC) and poly(A) are in the B-form, classical and slightly shifted to the A-form, respectively. The duplexes ATAT and GGTATACC have a strong and salt concentration dependent bias toward the A-form. The polymers poly(AC) and poly(G) take the C- and A-forms, respectively, at high salt concentrations. The reproduction of the behavior of these oligomers can serve as a test for the balance of interactions between the base stacking and the conformational flexibility of the sugar-phosphate backbone in a DNA force field. We tested the AMBER bsc1 and CHARMM36 force fields and their hybrids, and we failed to reproduce the experiment. In all the force fields, the salt concentration dependence is very weak. The known B-philicity of the AMBER force field proved to result from the B-philicity of its excessively strong base stacking. In the CHARMM force field, the B-form is a result of a fragile balance between the A-philic base stacking (especially for G:C pairs) and the C-philic backbone. Finally, we analyzed some recent simulations of the LacI-, SOX-4-, and Sac7d-DNA complex formation in the framework of the AMBER force field.

Toward Force Fields with Improved Base Stacking Descriptions.

Recent DNA force fields indicate good performance in describing flexibility and structural stability of double-stranded B-DNA. However, it is not clear how accurately base stacking interactions are represented that are critical for simulating structure formation processes and conformational changes. Based on the equilibrium nucleoside association and base pair nicking, we find that the recent Tumuc1 force field improves the description of base stacking compared to previous state-of-the-art force fields. Nevertheless, base pair stacking is still overstabilized compared to experiment. We propose a rapid method to reweight calculated free energies of stacking upon force field modifications in order to generate improved parameters. A decrease of the Lennard-Jones attraction between nucleo-bases alone appears insufficient; however, adjustments in the partial charge distribution on base atoms could help to further improve the force field description of base stacking.

An Attempt of Seeking Favorable Binding Free Energy Prediction Schemes Considering the Entropic Effect on Fis-DNA Binding.

Protein-DNA binding mechanisms in a complex manner are essential for understanding many biological processes. Over the past decades, numerous experiments and calculations have analyzed the specificity of protein-DNA binding. However, the accuracy of binding free energy prediction for multi-base DNA systems still needs to be improved. Fis is a DNA-binding protein that regulates various transcription and recombination reactions. In the present work, we tested several methods of predict binding free energy based on this system to find a favorable prediction scheme and explore the binding mechanism of Fis protein and DNA. Two solvent models (explicit and implicit solvent models) were chosen for the dynamics process, and the predicted binding free energy was more accurate under the explicit solvent model. When different Poisson-Boltzmann/Generalized Born (PB/GB) models were tested for DNA force fields (BSC1 and OL15), it was found that the binding free energy predicted by the selected OL15 force field performed better and the correlation between predicted and experimental values was improved with the increasing interior dielectric constant (Dk). Finally, using Dk = 8, the GBOBC1 model combined with interaction entropy (IE), which was calculated for entropic contribution (GBOBC1_IE_8), was screened out for the binding free energy prediction and analysis of the Fis-DNA system, and the validity of the method was further verified by testing the Cren7-DNA system. By performing conformational analysis of the minor groove, it was found that mutation of the DNA central sequence A/T to C/G and deletion of the guanine 2-amino group would change the minor groove width and thus affect the formation of the major groove, altering the interaction and atomic contact between the protein and the major groove, thus changing the binding affinity of Fis and DNA. Hopefully, the series of tests in this work can shed some light on the related studies of protein and DNA systems.

The marionette mechanism of domain–domain communication in the antagonist, agonist, and coactivator responses of the estrogen receptor

The human estrogen receptor α (hERα) is involved in the regulation of growth, development, and tissue homeostasis. Agonists that bind to the receptor's ligand-binding domain (LBD) lead to recruitment of coactivators and the enhancement of gene expression. In contrast, antagonists bind to the LBD and block the binding of coactivators thus decreasing gene expressions. In this work, we carry out simulations using the AWSEM (Associative memory, Water mediated, Structure and Energy Model)-Suite force field along with the 3SPN.2C force field for DNA to predict the structure of hERα and study its dynamics when binding to DNA and coactivators. Using simulations of antagonist-bound hERα and agonist-bound hERα by themselves and also along with bound DNA and coactivators, principal component analyses and free energy landscape analyses capture the pathway of domain-domain communication for agonist-bound hERα. This communication is mediated through the hinge domains that are ordinarily intrinsically disordered. These disordered segments manipulate the hinge domains much like the strings of a marionette as they twist in different ways when antagonists or agonists are bound to the ligand-binding domain.

Polarizable Molecular Block Model: Toward the Development of an Induced Dipole Force Field for DNA.

For flexible and highly ionized macromolecules such as DNA, it is important to correctly evaluate the intramolecular polarization in an induced dipole force field. In a proposed polarizable molecular block (PMB) model, a large molecule is divided into several molecular blocks. The atomic charges of the blocks are optimized by using the respective electrostatic potentials (ESPs) on the molecular surface. By using the capped hydrogen removal operation, the total charge of the blocks is controlled exactly to have an integer charge. The atomic polarizabilities of the blocks are optimized by using the respective polarized one-electron potentials that are the differences between ESPs with and without an external test charge. Induced dipole-charge interactions between the blocks are all included, but those interactions within the blocks are strictly excluded. All dipole-dipole interactions are included, but the damping functions are applied to the close dipole-dipole pairs. Several types of damping (simple scaling, exponential, linear, and Gaussian) are evaluated. The validity of the PMB model was verified by using trinucleotide duplexes which have A-, B-, and Z-DNA forms. The reference energies of trinucleotide duplexes including counterions (GGT3Na-ACC3Na, GAC3Na-GTC3Na, and GCG3Na-CGC3Na) are calculated using ωB97XD/aug-cc-pVDZ. All damping types reproduced well the reference interaction energies, dipole moments, and ESPs. Among them, the simple scaling with strong attenuation to 1-2 atomic pairs showed the highest stability against the polarization catastrophe. This study shows that it is possible to develop a high-quality polarizable force field by treating the intramolecular polarization on a block-by-block basis.

Tumuc1: A New Accurate DNA Force Field Consistent with High-Level Quantum Chemistry.

An accurate molecular mechanics force field forms the basis of Molecular Dynamics simulations to obtain a realistic view of the structure and dynamics of biomolecules such as DNA. Although frequently updated to improve agreement with available experimental data, DNA force fields still rely in part on parameters introduced more than 20 years ago. We have developed an entirely new DNA force field, Tumuc1, derived from quantum mechanical calculations to obtain a consistent set of bonded parameters and partial atomic charges. The performance of the force field was extensively tested on a variety of DNA molecules. It excels in accuracy of B-DNA simulations but also performs very well on other types of DNA structures and structure formation processes such as hairpin folding, duplex formation, and dynamics of DNA-protein complexes. It can complement existing force fields in order to provide an increasingly accurate description of the structure and dynamics of DNA during simulation studies.

Ethidium bromide interactions with DNA: an exploration of a classic DNA-ligand complex with unbiased molecular dynamics simulations.

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 >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.

Multipolar Description of Atom-Atom Electrostatic Interaction Energies in Single/Double-Stranded DNAs.

Molecular force field simulation is an effective method to explore the properties of DNA molecules in depth. Almost all current popular force fields calculate atom-atom electrostatic interaction energies for DNAs based on the atomic charge and dipole or quadrupole moments, without considering high-rank atomic multipole moments for more accurate electrostatics. Actually, the distribution of electrons around atomic nuclei is not spherically symmetric but is geometry dependent. In this work, a multipole expansion method that allows us to combine polarizability and anisotropy was applied. One single-stranded DNA and one double-stranded DNA were selected as pilot systems. Deoxynucleotides were cut out from pilot systems and capped by mimicking the original DNA environment. Atomic multipole moments were integrated instead of fixed-point charges to calculate atom-atom electrostatic energies to improve the accuracy of force fields for DNA simulations. Also, the applicability of modeling the behavior of both single-stranded and double-stranded DNAs was investigated. The calculation results indicated that the models can be transferred from pilot systems to test systems, which is of great significance for the development of future DNA force fields.

MARTINI-Based Protein-DNA Coarse-Grained HADDOCKing.

Modeling biomolecular assemblies is an important field in computational structural biology. The inherent complexity of their energy landscape and the computational cost associated with modeling large and complex assemblies are major drawbacks for integrative modeling approaches. The so-called coarse-graining approaches, which reduce the degrees of freedom of the system by grouping several atoms into larger “pseudo-atoms,” have been shown to alleviate some of those limitations, facilitating the identification of the global energy minima assumed to correspond to the native state of the complex, while making the calculations more efficient. Here, we describe and assess the implementation of the MARTINI force field for DNA into HADDOCK, our integrative modeling platform. We combine it with our previous implementation for protein-protein coarse-grained docking, enabling coarse-grained modeling of protein-nucleic acid complexes. The system is modeled using MARTINI topologies and interaction parameters during the rigid body docking and semi-flexible refinement stages of HADDOCK, and the resulting models are then converted back to atomistic resolution by an atom-to-bead distance restraints-guided protocol. We first demonstrate the performance of this protocol using 44 complexes from the protein-DNA docking benchmark, which shows an overall ~6-fold speed increase and maintains similar accuracy as compared to standard atomistic calculations. As a proof of concept, we then model the interaction between the PRC1 and the nucleosome (a former CAPRI target in round 31), using the same information available at the time the target was offered, and compare all-atom and coarse-grained models.

Hydrophobic unnatural base pairs show a Watson-Crick pairing in micro-second molecular dynamics simulations

Two unnatural hydrophobic nucleotides, d5SICS (2,6-dimethyl-2H-isoquiniline-1-thione) and dNaM (2-methoxy-3-methylnaphthalene), were previously replicated in vivo by a modified E. coli strand, however, a consistent structure for their pairing in terms of specific and selective directional interactions remains elusive, as data from spectroscopy experiments and simulations are inconsistent. The proposed d5SICS–dNaM pairing has been suggested to be a stacked configuration as suggested by NMR data; simulations have failed to reproduce this configuration and a Watson–Crick like pairing is observed. Previously, we focused on reproducing the d5SICS–dNaM Unnatural Base Pair (UBP) paring using an older (bsc0) AMBER force field, which was not able to correctly reproduce the experimental data. We present our efforts to reproduce the experimental pairing using the current version of the AMBER DNA force fields (OL15 and bsc1), two water models (OPC and TIP3P) and external electrostatic stabilization by Mg2+ ions. Opposite to previously reported simulations, a Watson–Crick-like pairing with no hydrogen bonds persists throughout all our results. Despite our efforts to replicate the reported stacked conformation, we cannot confirm its plausibility nor obtain a consistent structure that is independent of the neighboring nucleotides. Communicated by Ramaswamy H. Sarma

An In-Depth Look at DNA Crystals through the Prism of Molecular Dynamics Simulations

Summary X-ray crystallography is the primary tool for biomolecular structural determination. However, contacts formed through the crystal lattice are known to affect structures, especially for small and flexible molecules such as DNA oligomers, by introducing significant structural changes in comparison to solution. Furthermore, why molecules crystallize in certain symmetry groups, which role crystallization additives play, and whether they are just innocuous and unspecific crystallization catalysts remain unclear. By using one of the currently best-performing DNA force fields and applying significant computational effort, we described the nature of intermolecular forces that stabilize B-DNA crystals in various symmetry groups and solvent environments with an unprecedented level of detail. We showed a tight coupling between the lattice stability and the type of crystallization additives and that certain symmetry groups are stable only in the presence of a specific additive. Additives and crystal contacts induce small but non-negligible changes in the physical properties of DNA.

Does Cytosine Methylation Stabilize the BI Substate of DNA?

The local conformation of DNA plays a critical role in the recognition-and-binding process of transcription factors. In particular, the most common form of the DNA double helix, B-DNA, exists as a conformational equilibrium between BI and BII substates that are associated with changes in major and minor groove dimensions. It has been hypothesized that cytosine methylation, the most common epigenetic modification and one that is often found in so-called CpG islands, modulates transcription factor binding affinity by altering the BI-BII equilibrium in the neighborhood of methylation sites. Indeed, previous molecular dynamics studies have suggested that cytosine methylation leads to a uniform stabilization of the BI substate in CpG islands. In this work we first benchmark the ability of the latest Amber DNA force fields and a variety of three-point water models to reproduce BI-BII equilibrium in the Dickerson dodecamer and reduce terminal base pair fraying. Using the best-performing combinations, we then re-examine the conformational effects of cytosine methylation in a prototypical CpG island, (GC)5. We find that the BI substate is stabilized for GpC steps in such sequences, but that the BII substate is stabilized for the CpG steps. More interestingly, the Amber ff99bsc0 DNA force field used in previous studies yields opposite conclusions. These findings highlight the need for continued fixed-charge force field development for nucleic acids, as well as solution state experimental datasets involving epigenetic modification that can serve as benchmarks for simulation accuracy.

Molecular dynamics simulations on DNA flexibility: a comparative study of Amber bsc1 and bsc0 force fields

The structural flexibility of DNA plays a key role in many biological processes of DNA, such as protein-DNA interactions, DNA packaging in viruses and nucleosome positioning on genomic DNA. Some experimental techniques have been employed to investigate the structural flexibility of DNA with the combination of elastic models, but these experiments could only provide the macroscopic properties of DNA, and thus, it is still difficult to understand the corresponding microscopic mechanisms. Recently, all-atom molecular dynamics (MD) simulation has emerged as a useful tool to investigate not only the macroscopic properties of DNA, but also the microscopic description of the flexibility of DNA at an atomic level. The most important issue in all-atom MD simulations of DNA is to choose an appropriate force field for simulating DNA. Very recently, a new force field for DNA has been developed based on the last generation force field of Amber bsc0, which was named Amber bsc1. In this work, all-atom MD simulations are employed to study the flexibility of a 30-bp DNA with the force fields of Amber bsc1 and Amber bsc0 in a comparative way. Our aim of the research is to examine the improvement of the new development of force field (Amber bsc1) in the macroscopic and microscopic properties of DNA, in comparison with the corresponding experimental measurements. All the MD simulations are performed with Gromacs 4.6 and lasted with a simulation time of 600 ns. The MD trajectories are analyzed with Curves+ for the last 500 ns, since the system reaches equilibrium approximately after ~100 ns. Our results show that the new force field (Amber bsc1) can lead to the improvements in the macroscopic parameters of DNA flexibility, i.e., stretch modulus S and twist-stretch coupling D become closer to experimental measurements, while bending persistence lengths lp and torsional persistence lengths C from the two force fields (bsc1 and bsc0) are both in good agreement with experimental data. Our microscopic analyses show that the microscopic structure parameters of DNA from the MD simulation with the Amber bsc1 force field are closer to the experimental values than those with the Amber bsc0 force field, except for slide, and the obvious improvements are observed in some microscopic parameters such as twist and inclination. Our further analyses show that the improvements in macroscopic flexibility from the Amber bsc1 force field are tightly related to the microscopic parameters and their fluctuations. This study would be helpful in understanding the performances of Amber bsc1 and bsc0 force fields in the description of DNA flexibility at both macroscopic and microscopic level.

A- to B-DNA Transition in AMBER Force Fields and Its Coupling to Sugar Pucker.

The A/B transition is a basic element of DNA conformational change. Because of its involvement in the sensing of the ionic conditions by DNA and in specific protein-DNA interactions, this transition is important for biological functions of DNA. Therefore, accurate modeling of the A/B equilibrium by means of empirical force fields is of utmost interest. In this work, we examine the A/B equilibrium in three AMBER force fields, including the recent bsc1 and OL15 modifications, using much longer MD simulations than attempted before. Special attention is paid to the coupling of the A/B equilibrium with the south/north (S/N) transition of the sugar pucker. We found that none of the tested force fields provided a satisfactory description of the A/B equilibrium because the B-form was predicted to be much too stable and the A-form was predicted to be almost absent even in concentrated trifluoroethanol solutions. Based on comparison with NMR data for duplexes and single nucleosides, we hypothesize that this problem arose from the incorrect description of the S/N equilibrium of sugar pucker, where the south conformation is much too stable, thus stabilizing the B-form. Because neither the A/B equilibrium in duplexes nor the S/N equilibrium in nucleosides was described accurately, further refinements of the AMBER DNA force fields are needed.

Resolving the Conformational Dynamics of DNA with Ångstrom Resolution by Pulsed Electron-Electron Double Resonance and Molecular Dynamics.

Pulsed electron-electron double resonance (PELDOR/DEER) experiments of nucleic acids with rigid spin labels provide highly accurate distance and orientation information. Here we combine PELDOR experiments with molecular dynamics (MD) simulations to arrive at an atomistic view of the conformational dynamics of DNA. The MD simulations closely reproduce the PELDOR time traces, and demonstrate that bending, in addition to twist-stretch motions, underpin the sub-μs dynamics of DNA. PELDOR experiments correctly rank DNA force fields and resolve subtle differences in the conformational ensembles of nucleic acids, on the order of 1-2 Å. Long-range distance and angle measurements with rigid spin labels provide critical input for the refinement of computer models and the elucidation of the structure and dynamics of complex biomolecules.

Polarizable Force Field for DNA Based on the Classical Drude Oscillator: II. Microsecond Molecular Dynamics Simulations of Duplex DNA.

The structure and dynamics of DNA are governed by a sensitive balance between base stacking and pairing, hydration, and interactions with ions. Force-field models that include explicit representations of electronic polarization are capable of more accurately modeling the subtle details of these interactions versus commonly used additive force fields. In this work, we validate our recently refined polarizable force field for DNA based on the classical Drude oscillator model, in which electronic degrees of freedom are represented as negatively charged particles attached to their parent atoms via harmonic springs. The previous version of the force field, called Drude-2013, produced stable A- and B-DNA trajectories on the order of hundreds of nanoseconds, but deficiencies were identified that included weak base stacking ultimately leading to distortion of B-DNA duplexes and unstable Z-DNA. As a result of extensive refinement of base nonbonded terms and bonded parameters in the deoxyribofuranose sugar and phosphodiester backbone, we demonstrate that the new version of the Drude DNA force field is capable of simulating A- and B-forms of DNA on the microsecond time scale and the resulting conformational ensembles agree well with a broad set of experimental properties, including solution X-ray scattering profiles. In addition, simulations of Z-form duplex DNA in its crystal environment are stable on the order of 100 ns. The revised force field is to be called Drude-2017.

Polarizable Force Field for DNA Based on the Classical Drude Oscillator: I. Refinement Using Quantum Mechanical Base Stacking and Conformational Energetics.

Empirical force fields seek to relate the configuration of a set of atoms to its energy, thus yielding the forces governing its dynamics, using classical physics rather than more expensive quantum mechanical calculations that are computationally intractable for large systems. Most force fields used to simulate biomolecular systems use fixed atomic partial charges, neglecting the influence of electronic polarization, instead making use of a mean-field approximation that may not be transferable across environments. Recent hardware and software developments make polarizable simulations feasible, and to this end, polarizable force fields represent the next generation of molecular dynamics simulation technology. In this work, we describe the refinement of a polarizable force field for DNA based on the classical Drude oscillator model by targeting quantum mechanical interaction energies and conformational energy profiles of model compounds necessary to build a complete DNA force field. The parametrization strategy employed in the present work seeks to correct weak base stacking in A- and B-DNA and the unwinding of Z-DNA observed in the previous version of the force field, called Drude-2013. Refinement of base nonbonded terms and reparametrization of dihedral terms in the glycosidic linkage, deoxyribofuranose rings, and important backbone torsions resulted in improved agreement with quantum mechanical potential energy surfaces. Notably, we expand on previous efforts by explicitly including Z-DNA conformational energetics in the refinement.

Can We Execute Reliable MM-PBSA Free Energy Computations of Relative Stabilities of Different Guanine Quadruplex Folds?

The self-assembly and stability of DNA G-quadruplexes (GQs) are affected by the intrinsic stability of different GpG base steps embedded in their G-quartet stems. We have carried out MD simulations followed by MM-PBSA (molecular mechanics Poisson-Boltzmann surface area) free energy calculations on all the experimentally observed three-quartet intramolecular human telomeric GQ topologies. We also studied antiparallel GQ models with alternative syn-anti patterns of the G-quartets. We tested different ions, dihedral variants of the DNA force field, water models, and simulation lengths. In total, ∼35 μs of simulations have been carried out. The systems studied here are considerably more complete than the previously analyzed two-quartet stems. Among other effects, our computations included the stem-loop coupling and ion-ion interactions inside the stem. The calculations showed a broad agreement with the earlier predictions. However, the increase in the completeness of the system was associated with increased noise of the free energy data which could be related, for example, to the presence of long-lived loop substates and rather complex dynamics for the two bound ions inside the G-stem. As a result, the MM-PBSA data were noisy and we could not improve their quantitative convergence even by expanding the simulations to 2.5 μs long trajectories. We also suggest that the quality of MM-based free energy computations based on MD simulations of complete GQs is more sensitive to the neglect of explicit polarization effects, which, in real systems, are associated with the presence of multiple closely spaced ions inside the GQs. Thus, although the MM-PBSA procedure provides very useful insights that complement the structural-dynamics data from MD trajectories of GQs, the method is far from reaching quantitative accuracy. Our conclusions are in agreement with critical assessments of the MM-PBSA methodology available in contemporary literature for other types of problems.

CHARMM Drude Polarizable MD Simulations Reproduce Solution X-Ray Diffraction Patterns for B-DNA Sequences and Predict Differential Impact of the Li+, Na+, K+ and Rb+ Ions on DNA Conformational Properties

Recently we have presented the first generation CHARMM Drude polarizable force field for DNA, which is capable of reproducing the main conformational features of DNA in solution, such as A-to-B equilibrium and transitions between BI and BII substates. Our current efforts are directed towards further model improvement, by achieving a proper balance of the interactions among surrounding mobile ions, water and DNA. Compared to the additive (non-polarizable) models, explicit treatment of the electronic polarizability in the Drude model leads to a markedly improved description of the interplay between the ionic atmosphere and DNA conformational behavior. In particular, the Drude model is shown to more accurately reproduce counterion condensation theory predictions of DNA charge neutralization by condensed ions, as well as the experimental data on the competitive binding of Li+, Na+, K+ and Rb+ ions to DNA. The most intriguing results is that the model predicts a differential impact of these seemingly similar monovalent cations on DNA conformational properties - a phenomenon not observed in the state-of-the-art CHARMM36 and AMBER atomistic additive models. In addition, the Drude model reproduces the solution X-ray diffraction patterns for a number of B form DNA sequences at a level of accuracy similar to, or exceeding that of the above mentioned additive models. The obtained results indicate that CHARMM Drude polarizable MD simulations provide a more realistic model of the physical forces involved in the interactions of DNA with its ionic environment, offering the potential to yield new insights into salt-mediated biological processes involving DNA, such as protein-DNA recognition and chromatin folding.

An improved DNA force field for ssDNA interactions with gold nanoparticles.

The widespread applications of single-stranded DNA (ssDNA) conjugated gold nanoparticles (AuNPs) have spurred an increasing interest in the interactions between ssDNA and AuNPs. Despite extensive studies using the most sophisticated experimental techniques, the detailed molecular mechanisms still remain largely unknown. Large scale molecular dynamics (MD) simulations can thus be used to supplement experiments by providing complementary information about ssDNA-AuNP interactions. However, up to now, all modern force fields for DNA were developed based on the properties of double-stranded DNA (dsDNA) molecules, which have hydrophilic outer backbones "protecting" hydrophobic inner nucleobases from water. Without the double-helix structure of dsDNA and thus the "protection" by the outer backbone, the nucleobases of ssDNA are directly exposed to solvent, and their behavior in water is very different from that of dsDNA, especially at the interface with nanoparticles. In this work, we have improved the force field of ssDNA for use with nanoparticles, such as AuNPs, based on recent experimental results and quantum mechanics calculations. With the new improved force field, we demonstrated that a poly(A) sequence adsorbed on a AuNP surface is much more stable than a poly(T) sequence, which is consistent with recent experimental observations. On the contrary, the current standard force fields, including AMBER03, CHARMM27, and OPLSAA, all gave erroneous results as compared to experiments. The current improved force field is expected to have wide applications in the study of ssDNA with nanomaterials including AuNPs, which might help promote the development of ssDNA-based biosensors and other bionano-devices.