Dataset Information

Polarizable empirical force field for the primary and secondary alcohol series based on the classical Drude model.

ABSTRACT: A polarizable empirical force field based on the classical Drude oscillator has been developed for the aliphatic alcohol series. The model is optimized with emphasis on condensed-phase properties and is validated against a variety of experimental data. Transferability of the developed parameters is emphasized by the use of a single electrostatic model for the hydroxyl group throughout the alcohol series. Aliphatic moiety parameters were transferred from the polarizable alkane parameter set, with only the Lennard-Jones parameters on the carbon in methanol optimized. The developed model yields good agreement with pure solvent properties with the exception of the heats of vaporization of 1-propanol and 1-butanol, which are underestimated by approximately 6%; special LJ parameters for the oxygen in these two molecules that correct for this limitation are presented. Accurate treatment of the free energies of aqueous solvation required the use of atom-type specific O(alcohol)-O(water) LJ interaction terms, with specific terms used for the primary and secondary alcohols. With respect to gas phase properties the polarizable model overestimates experimental dipole moments and quantum mechanical interaction energies with water by approximately 10 and 8 %, respectively, a significant improvement over 44 and 46 % overestimations of the corresponding properties in the CHARMM22 fixed-charge additive model. Comparison of structural properties of the polarizable and additive models for the pure solvents and in aqueous solution shows significant differences indicating atomic details of intermolecular interactions to be sensitive to the applied force field. The polarizable model predicts pure solvent and aqueous phase dipole moment distributions for ethanol centered at 2.4 and 2.7 D, respectively, a significant increase over the gas phase value of 1.8 D, whereas in a solvent of lower polarity, benzene, a value of 1.9 is obtained. The ability of the polarizable model to yield changes in dipole moment as well as the reproduction of a range of condensed phase properties indicates its utility in the study of the properties of alcohols in a variety of condensed phase environments as well as representing an important step in the development of a comprehensive force field for biological molecules.

SUBMITTER: Anisimov VM

PROVIDER: S-EPMC2542883 | biostudies-literature | 2007

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Polarizable empirical force field for the primary and secondary alcohol series based on the classical Drude model.

Anisimov Victor M VM Vorobyov Igor V IV Roux Benoît B Mackerell Alexander D AD

Journal of chemical theory and computation 20070101 6

A polarizable empirical force field based on the classical Drude oscillator has been developed for the aliphatic alcohol series. The model is optimized with emphasis on condensed-phase properties and is validated against a variety of experimental data. Transferability of the developed parameters is emphasized by the use of a single electrostatic model for the hydroxyl group throughout the alcohol series. Aliphatic moiety parameters were transferred from the polarizable alkane parameter set, with ...[more]

PMID: 18802495

Similar Datasets

Project description:The CHARMM Drude-2013 polarizable force field (FF) was developed to include the explicit treatment of induced electronic polarizability, resulting in a more accurate description of the electrostatic interactions in molecular dynamics (MD) simulations. While the Drude-2013 protein FF has shown success in improving the folding properties of α-helical peptides and to reproduce experimental observables in simulations up to 1 μs, some limitations were noted regarding the stability of β-sheet structures in simulations longer than 100 ns as well as larger deviations from crystal structures in simulations of a number of proteins compared to the additive CHARMM36 protein FF. The origin of the instability has been identified and appears to be primarily due to overestimated atomic polarizabilities and induced dipole-dipole interactions on the Cβ, Cγ, and Cδ side chain atoms. To resolve this and other issues, a number of aspects of the model were revisited, resulting in Drude-2019 protein FF. Backbone parameters were optimized targeting the conformational properties of the (Ala)5 peptide in solution along with gas phase properties of the alanine dipeptide. Dipeptides that contain N-acetylated and N'-methylamidated termini, excluding Gly, Pro, and Ala, were used as models to optimize the atomic polarizabilities and Thole screening factors on selected Cβ, Cγ, and Cδ carbons by targeting quantum mechanical (QM) dipole moments and molecular polarizabilities. In addition, to obtain better conformational properties, side chain χ1 and χ2 dihedral parameters were optimized targeting QM data for the respective side chain dipeptide conformations as well as Protein Data Bank survey data based on the χ1, χ2 sampling from Hamiltonian replica-exchange MD simulations of (Ala)4-X-(Ala)4 in solution, where X is the amino acid of interest. Further improvements include optimizing nonbonded interactions between charged residues to reproduce QM interaction energies of the charged-protein model compounds and experimental osmotic pressures. Validation of the optimized Drude protein FF includes MD simulations of a collection of peptides and proteins including β-sheet structures, as well as transmembrane ion channels. Results showed that the updated Drude-2019 protein FF yields smaller overall root-mean-square differences of proteins as compared to the additive CHARMM36m and Drude-2013 FFs as well as similar or improved agreement with experimental NMR properties, allowing for long time scale simulation studies of proteins and more complex biomolecular systems in conjunction with the remainder of the Drude polarizable FF.

Dataset Information

Polarizable empirical force field for the primary and secondary alcohol series based on the classical Drude model.

Publications

Polarizable empirical force field for the primary and secondary alcohol series based on the classical Drude model.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets