@article{babin_baucom_darden_sagui_2006, title={Molecular dynamics simulations of DNA with polarizable force fields: Convergence of an ideal B-DNA structure to the crystallographic structure}, volume={110}, ISSN={["1520-5207"]}, DOI={10.1021/jp061421r}, abstractNote={We have investigated to what extent molecular dynamics (MD) simulations can reproduce DNA sequence-specific features, given different electrostatic descriptions and different cell environments. For this purpose, we have carried out multiple unrestrained MD simulations of the DNA duplex d(CCAACGTTGG)2. With respect to the electrostatic descriptions, two different force fields are studied: a traditional description based on atomic point charges and a polarizable force field. With respect to the cell environment, the difference between crystal and solution environments is emphasized, as well as the structural importance of divalent ions. By imposing the correct experimental unit cell environment, an initial configuration with two ideal B-DNA duplexes in the unit cell is shown to converge to the crystallographic structure. This convergence is measured by the appearance of sequence-dependent features that very closely resemble the crystallographic ones as well as by the decay of the all-atom root-mean-squared coordinates deviations (RMSD) with respect to the crystallographic structure. Given the appropriate crystallographic constraints, this is the first example of multiple nanosecond molecular dynamics trajectory that shows an ideal B-DNA model converging to an experimental structure, with a significant decay of RMSD.}, number={23}, journal={JOURNAL OF PHYSICAL CHEMISTRY B}, author={Babin, Volodymyr and Baucom, Jason and Darden, Thomas A. and Sagui, Celeste}, year={2006}, month={Jun}, pages={11571–11581} } @article{babin_baucom_darden_sagui_2006, title={Molecular dynamics simulations of polarizable DNA in crystal environment}, volume={106}, ISSN={0020-7608 1097-461X}, url={http://dx.doi.org/10.1002/qua.21152}, DOI={10.1002/qua.21152}, abstractNote={Abstract We have investigated the role of the electrostatic description and cell environment in molecular dynamics (MD) simulations of DNA. Multiple unrestrained MD simulations of the DNA duplex d(CCAACGTTGG) 2 have been carried out using two different force fields: a traditional description based on atomic point charges and a polarizable force field. For the time scales probed, and given the “right” distribution of divalent ions, the latter performs better than the nonpolarizable force field. In particular, by imposing the experimental unit cell environment, an initial configuration with ideal B‐DNA duplexes in the unit cell acquires sequence‐dependent features that very closely resemble the crystallographic ones. Simultaneously, the all‐atom root‐mean‐square coordinates deviation (RMSD) with respect to the crystallographic structure is seen to decay. At later times, the polarizable force field is able to maintain this lower RMSD, while the nonpolarizable force field starts to drift away. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006}, number={15}, journal={International Journal of Quantum Chemistry}, publisher={Wiley}, author={Babin, Volodymyr and Baucom, Jason and Darden, Thomas A. and Sagui, Celeste}, year={2006}, pages={3260–3269} } @article{baucom_transue_fuentes-cabrera_krahn_darden_sagui_2004, title={Molecular dynamics simulations of the d(CCAACGTTGG)(2) decamer in crystal environment: Comparison of atomic point-charge, extra-point, and polarizable force fields}, volume={121}, ISSN={["1089-7690"]}, DOI={10.1063/1.1788631}, abstractNote={Molecular dynamics simulations of the DNA duplex d(CCAACGTTGG)(2) were used to study the relationship between DNA sequence and structure in a crystal environment. Three different force fields were used: a traditional description based on atomic point charges, a polarizable force field, and an "extra-point" force field (with additional charges on extranuclear sites). It is found that all the force fields reproduce fairly well the sequence-dependent features of the experimental structure. The polarizable force field, however, provides the most accurate representation of the crystal structure and the sequence-dependent effects observed in the experiment. These results point out to the need of the inclusion of polarization for accurate descriptions of DNA.}, number={14}, journal={JOURNAL OF CHEMICAL PHYSICS}, author={Baucom, J and Transue, T and Fuentes-Cabrera, M and Krahn, JM and Darden, TA and Sagui, C}, year={2004}, month={Oct}, pages={6998–7008} }