Dataset Information

Polarizable Force Field for Molecular Ions Based on the Classical Drude Oscillator.

ABSTRACT: Development of accurate force field parameters for molecular ions in the context of a polarizable energy function based on the classical Drude oscillator is a crucial step toward an accurate polarizable model for modeling and simulations of biological macromolecules. Toward this goal we have undertaken a hierarchical approach in which force field parameter optimization is initially performed for small molecules for which experimental data exists that serve as building blocks of macromolecular systems. Small molecules representative of the ionic moieties of biological macromolecules include the cationic ammonium and methyl substituted ammonium derivatives, imidazolium, guanidinium and methylguanidinium, and the anionic acetate, phenolate, and alkanethiolates. In the present work, parameters for molecular ions in the context of the Drude polarizable force field are optimized and compared to results from the nonpolarizable additive CHARMM general force field (CGenFF). Electrostatic and Lennard-Jones parameters for the model compounds are developed in the context of the polarizable SWM4-NDP water model, with emphasis on assuring that the hydration free energies are consistent with previously reported parameters for atomic ions. The final parameters are shown to be in good agreement with the selected quantum mechanical (QM) and experimental target data. Analysis of the structure of water around the ions reveals substantial differences between the Drude and additive force fields indicating the important role of polarization in dictating the molecular details of aqueous solvation. The presented parameters represent the foundation for the charged functionalities in future generations of the Drude polarizable force field for biological macromolecules as well as for drug-like molecules.

SUBMITTER: Lin FY

PROVIDER: S-EPMC5975207 | biostudies-literature | 2018 May

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Polarizable Force Field for Molecular Ions Based on the Classical Drude Oscillator.

Lin Fang-Yu FY Lopes Pedro E M PEM Harder Edward E Roux Benoît B MacKerell Alexander D AD

Journal of chemical information and modeling 20180417 5

Development of accurate force field parameters for molecular ions in the context of a polarizable energy function based on the classical Drude oscillator is a crucial step toward an accurate polarizable model for modeling and simulations of biological macromolecules. Toward this goal we have undertaken a hierarchical approach in which force field parameter optimization is initially performed for small molecules for which experimental data exists that serve as building blocks of macromolecular sy ...[more]

PMID: 29624370

Similar Datasets

Project description:A polarizable force field for nucleic acid bases based on the classical Drude oscillator model is presented. Parameter optimization was performed to reproduce crystallographic geometries, crystal unit cell parameters, heats of sublimation, vibrational frequencies and assignments, dipole moments, molecular polarizabilities and quantum mechanical base-base and base-water interaction energies. The training and validation data included crystals of unsubstituted and alkyl-substituted adenine, guanine, cytosine, uracil, and thymine bases, hydrated crystals, and hydrogen bonded base pairs. Across all compounds, the RMSD in the calculated heats of sublimation is 4.1%. This equates to an improvement of more than 2.5 kcal/mol in accuracy compared to the nonpolarizable CHARMM27 force field. However, the level of agreement with experimental molecular volume decreased from 1.7% to 2.1% upon moving from the nonpolarizable to the polarizable model. The representation of dipole moments is significantly improved with the Drude polarizable force field. Unlike in additive force fields, there is no requirement for the gas-phase dipole moments to be overestimated, illustrating the ability of the Drude polarizable force field to treat accurately differently dielectric environments and indicating the improvements in the electrostatic model. Validation of the model was performed on the basis of the calculation of the gas-phase binding enthalpies of base pairs obtained via potential of mean force calculations; the additive and polarizable models both performed satisfactorily with average differences of 0.2 and 0.9 kcal/mol, respectively, and rms differences of 1.3 and 1.7 kcal/mol, respectively. Overall, considering the number of significant improvements versus the additive CHARMM force field, the incorporation of explicit polarizability into the force field for nucleic acid bases represents an additional step toward accurate computational modeling of biological systems.

Dataset Information

Polarizable Force Field for Molecular Ions Based on the Classical Drude Oscillator.

Publications

Polarizable Force Field for Molecular Ions Based on the Classical Drude Oscillator.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure