Polarizable Empirical Force Field Research Articles

Improved Modeling of Cation-π and Anion-Ring Interactions Using the Drude Polarizable Empirical Force Field for Proteins.

Cation-π interactions are noncovalent interactions between a π-electron system and a positively charged ion that are regarded as a strong noncovalent interaction and are ubiquitous in biological systems. Similarly, though less studied, anion-ring interactions are present in proteins along with in-plane interactions of anions with aromatic rings. As these interactions are between a polarizing ion and a polarizable π system, the accuracy of the treatment of these interactions in molecular dynamics (MD) simulations using additive force fields (FFs) may be limited. In the present work, to allow for a better description of ion-π interactions in proteins in the Drude-2013 protein polarizable FF, we systematically optimized the parameters for these interactions targeting model compound quantum mechanical (QM) interaction energies with atom pair-specific Lennard-Jones parameters along with virtual particles as selected ring centroids introduced to target the QM interaction energies and geometries. Subsequently, MD simulations were performed on a series of protein structures where ion-π pairs occur to evaluate the optimized parameters in the context of the Drude-2013 FF. The resulting FF leads to a significant improvement in reproducing the ion-π pair distances observed in experimental protein structures, as well as a smaller root-mean-square differences and fluctuations of the overall protein structures from experimental structures. Accordingly, the optimized Drude-2013 protein polarizable FF is suggested for use in MD simulations of proteins where cation-π and anion-ring interactions are critical. © 2019 Wiley Periodicals, Inc.

Solvation dynamics: improved reproduction of the time-dependent Stokes shift with polarizable empirical force field chromophore models.

The inclusion of explicit polarization in molecular dynamics simulation has gained increasing interest during the last several years. An understudied area is the role of polarizability in computer simulations of solvation dynamics around chromophores, particularly for the large solutes used in experimental studies. In this work, we present fully polarizable ground and excited state force fields for the common fluorophores N-methyl-6-oxyquinolium betaine and Coumarin 153. While analyzing the solvation responses in water, methanol, and the highly viscous ionic liquid 1-ethyl-3-methylimidazolium trifluoromethanesulfonate we found that the inclusion of solute polarizability considerably increases the agreement of the obtained Stokes shift relaxation functions with experimental data. Solute polarizability slows down the inertial solvation response in the femtosecond time regime and enables the chromophore to adapt its dipole moment to the environment. Furthermore, the developed chromophore force field reproduces the solute dipole moments in both the electronic ground and excited state in environments ranging from gas phase to highly polar media correctly. Based on these studies it is anticipated that polarizable models of chromophores will lead to an improved understanding of the relationship of their environment to their spectroscopic properties.

Polarizable Empirical Force Field for Halogen-Containing Compounds Based on the Classical Drude Oscillator.

The quality of the force field is crucial to ensure the accuracy of simulations used in molecular modeling, including computer-aided drug design (CADD). To perform more accurate modeling and simulations of halogenated molecules, in this study the polarizable force field based on the classical Drude oscillator model was extended to both aliphatic and aromatic systems using halogenated ethane and benzene model compounds for the halogens F, Cl, Br, and I. The force field parameters were optimized targeting quantum mechanical dipole moments, water interactions, and molecular polarizabilities as well as experimental observables, including enthalpies of vaporization, molecular volumes, hydration free energies, and dielectric constants. The developed halogenated polarizable force field is capable of reproducing QM relative energies and geometries of both halogen bonds and halogen-hydrogen bond donor interactions at an unprecedented level due to the inclusion of a virtual particle and anisotropic atomic polarizability on the halogen and, notably, the inclusion of Lennard-Jones parameters on the halogen Drude particle. The model was validated on the basis of its ability to accurately reproduce pure solvent properties for halogenated naphthalenes and alkanes, including species analogous to those used as refrigerants. Accordingly, it is anticipated that the model will be applicable for the study of halogenated derivatives in CADD as well as in other chemical and biophysical studies.

Polarizable Empirical Force Field for Hexopyranose Monosaccharides Based on the Classical Drude Oscillator

A polarizable empirical force field based on the classical Drude oscillator is presented for the hexopyranose form of selected monosaccharides. Parameter optimization targeted quantum mechanical (QM) dipole moments, solute–water interaction energies, vibrational frequencies, and conformational energies. Validation of the model was based on experimental data on crystals, densities of aqueous-sugar solutions, diffusion constants of glucose, and rotational preferences of the exocylic hydroxymethyl of d-glucose and d-galactose in aqueous solution as well as additional QM data. Notably, the final model involves a single electrostatic model for all sixteen diastereomers of the monosaccharides, indicating the transferability of the polarizable model. The presented parameters are anticipated to lay the foundation for a comprehensive polarizable force field for saccharides that will be compatible with the polarizable Drude parameters for lipids and proteins, allowing for simulations of glycolipids and glycoproteins.

Polarizable empirical force field for acyclic polyalcohols based on the classical Drude oscillator.

A polarizable empirical force field for acyclic polyalcohols based on the classical Drude oscillator is presented. The model is optimized with an emphasis on the transferability of the developed parameters among molecules of different sizes in this series and on the condensed-phase properties validated against experimental data. The importance of the explicit treatment of electronic polarizability in empirical force fields is demonstrated in the cases of this series of molecules with vicinal hydroxyl groups that can form cooperative intra- and intermolecular hydrogen bonds. Compared to the CHARMM additive force field, improved treatment of the electrostatic interactions avoids overestimation of the gas-phase dipole moments resulting in significant improvement in the treatment of the conformational energies and leads to the correct balance of intra- and intermolecular hydrogen bonding of glycerol as evidenced by calculated heat of vaporization being in excellent agreement with experiment. Computed condensed phase data, including crystal lattice parameters and volumes and densities of aqueous solutions are in better agreement with experimental data as compared to the corresponding additive model. Such improvements are anticipated to significantly improve the treatment of polymers in general, including biological macromolecules.

Polarizable empirical force field for sulfur-containing compounds based on the classical Drude oscillator model.

Condensed-phase computational studies of molecules using molecular mechanics approaches require the use of force fields to describe the energetics of the systems as a function of structure. The advantage of polarizable force fields over nonpolarizable (or additive) models lies in their ability to vary their electronic distribution as a function of the environment. Toward development of a polarizable force field for biological molecules, parameters for a series of sulfur-containing molecules are presented. Parameter optimization was performed to reproduce quantum mechanical and experimental data for gas phase properties including geometries, conformational energies, vibrational spectra, and dipole moments as well as for condensed phase properties such as heats of vaporization, molecular volumes, and free energies of hydration. Compounds in the training set include methanethiol, ethanethiol, propanethiol, ethyl methyl sulfide, and dimethyl disulfide. The molecular volumes and heats of vaporization are in good accordance with experimental values, with the polarizable model performing better than the CHARMM22 nonpolarizable force field. Improvements with the polarizable model were also obtained for molecular dipole moments and in the treatment of intermolecular interactions as a function of orientation, in part due to the presence of lone pairs and anisotropic atomic polarizability on the sulfur atoms. Significant advantage of the polarizable model was reflected in calculation of the dielectric constants, a property that CHARMM22 systematically underestimates. The ability of this polarizable model to accurately describe a range of gas and condensed phase properties paves the way for more accurate simulation studies of sulfur-containing molecules including cysteine and methionine residues in proteins.

Structure of Large Nitrate−Water Clusters at Ambient Temperatures: Simulations with Effective Fragment Potentials and Force Fields with Implications for Atmospheric Chemistry

Structural properties of large NO(3)(-).(H(2)O)(n) (n = 15-500) clusters are studied by Monte Carlo simulations using effective fragment potentials (EFPs) and by classical molecular dynamics simulations using a polarizable empirical force field. The simulation results are analyzed with a focus on the description of hydrogen bonding and solvation in the clusters. In addition, a comparison between the electronic structure based EFP and the classical force field description of the 32 water cluster system is presented. The EFP simulations, which focused on the cases of n = 15 and 32, show an internal, fully solvated structure and a "surface adsorbed" structure for the 32 water cluster at 300 K, with the latter configuration being more probable. The internal solvated structure and the "surface adsorbed" structure differ considerably in their hydrogen bonding coordination numbers. The force field based simulations agree qualitatively with these results, and the local geometry of NO(3)(-) and solvation at the surface-adsorbed site in the force field simulations are similar to those predicted using EFPs. Differences and similarities between the description of hydrogen bonding of the anion in the two approaches are discussed. Extensive classical force field based simulations at 250 K predict that long time scale stability of "internal" NO(3)(-), which is characteristic of extended bulk aqueous interfaces, emerges only for n > 300. Ab initio Møller-Plesset perturbation theory is used to test the geometries of selected surface and interior anions for n = 32, and the results are compared to the EFP and MD simulations. Qualitatively, all approaches agree that surface structures are preferred over the interior structures for clusters of this size. The relatively large aqueous clusters of NO(3)(-) studied here are of comparable size to clusters that lead to new particle formation in air. Nitrate ions on the surface of such clusters may have significantly different photochemistry than the internal species. The possible implications of surface-adsorbed nitrate ions for atmospheric chemistry are discussed.

Polarizable empirical force field for nitrogen‐containing heteroaromatic compounds based on the classical Drude oscillator

The polarizable empirical CHARMM force field based on the classical Drude oscillator has been extended to the nitrogen-containing heteroaromatic compounds pyridine, pyrimidine, pyrrole, imidazole, indole, and purine. Initial parameters for the six-membered rings were based on benzene with nonbond parameter optimization focused on the nitrogen atoms and adjacent carbons and attached hydrogens. In the case of five-member rings, parameters were first developed for imidazole and transferred to pyrrole. Optimization of all parameters was performed against an extensive set of quantum mechanical and experimental data. Ab initio data were used for the determination of initial electrostatic parameters, the vibrational analysis, and in the optimization of the relative magnitudes of the Lennard-Jones (LJ) parameters, through computations of the interactions of dimers of model compounds, model compound-water interactions, and interactions of rare gases with model compounds. The absolute values of the LJ parameters were determined targeting experimental heats of vaporization, molecular volumes, heats of sublimation, crystal lattice parameters, and free energies of hydration. Final scaling of the polarizabilities from the gas-phase values by 0.85 was determined by reproduction of the dielectric constants of pyridine and pyrrole. The developed parameter set was extensively validated against additional experimental data such as diffusion constants, heat capacities, and isothermal compressibilities, including data as a function of temperature.

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

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.

Polarizable empirical force field for aromatic compounds based on the classical drude oscillator.

The polarizable empirical CHARMM force field based on the classical Drude oscillator has been extended to the aromatic compounds benzene and toluene. Parameters were optimized for benzene and then transferred directly to toluene, with parameters for the methyl moiety of toluene taken from the previously published work on the alkanes. Optimization of all parameters was performed against an extensive set of quantum mechanical and experimental data. Ab initio data was used for determination of the electrostatic parameters, for the vibrational analysis, and in the optimization of the relative magnitudes of the Lennard-Jones parameters. The absolute values of the Lennard-Jones parameters were determined by comparing computed and experimental heats of vaporization, molecular volumes, free energies of hydration, and dielectric constants. The newly developed parameter set was extensively tested against additional experimental data such as diffusion constants, heat capacities at constant pressure, and isothermal compressibilities including data as a function of temperature. Moreover, the structures of liquid benzene, liquid toluene, and solutions of each in water were studied. In the case of benzene, the computed and experimental total distribution function were compared, with the developed model shown to be in excellent agreement with experiment.

Polarizable Empirical Force Field for Alkanes Based on the Classical Drude Oscillator Model

Recent extensions of potential energy functions used in empirical force field calculations have involved the inclusion of electronic polarizability. To properly include this extension into a potential energy function it is necessary to systematically and rigorously optimize the associated parameters based on model compounds for which extensive experimental data are available. In the present work, optimization of parameters for alkanes in a polarizable empirical force field based on a classical Drude oscillator is presented. Emphasis is placed on the development of parameters for CH3, CH2, and CH moieties that are directly transferable to long chain alkanes, as required for lipids and other biomolecules. It is shown that a variety of quantum mechanical and experimental target data are reproduced by the polarizable model. Notable is the proper treatment of the dielectric constant of pure alkanes by the polarizable force field, a property essential for the accurate treatment of, for example, hydrophobic solvation in lipid bilayers. The present alkane force field will act as the basis for the aliphatic moieties in an extensive empirical force field for biomolecules that includes the explicit treatment of electronic polarizability.

Determination of Electrostatic Parameters for a Polarizable Force Field Based on the Classical Drude Oscillator.

A procedure to determine the electrostatic parameters has been developed for a polarizable empirical force field based on the classical Drude oscillator model. Atomic charges and polarizabilities for a given molecule of interest were derived from restrained fitting to quantum-mechanical electrostatic potentials (ESP) calculated at the B3LYP/ cc-pVDZ or B3LYP/aug-cc-pVDZ levels on grid points located on concentric Connolly surfaces. The determination of the atomic polarizabilities requires a series of perturbed ESP maps, each one representing the electronic response of the molecule in the presence of a background charge placed on Connolly surfaces primarily along chemical bonds and lone pairs. Reference values for the partial atomic charges were taken from the CHARMM27 additive all-atom force field, and those for the polarizabilities were based on adjusted Miller's ahp atomic polarizability values. The fitted values of atomic polarizabilities were scaled to reflect the reduced polarization expected for the condensed media and/or to correct for the systematic underestimation of experimental molecular polarizabilities by B3LYP calculations. Following correction of the polarizabilities, the atomic charges were adjusted to reproduce gas-phase dipole moments. The developed scheme has been tested on a set of small molecules representing functional moieties of nucleic acids. The derived electrostatic parameters have been successfully applied in a preliminary polarizable molecular dynamics simulation of a DNA octamer in a box of water with sodium counterions. Thus, this study confirms the feasibility of the use of a polarizable force field based on a classical Drude model for simulations of biomolecules in the condensed phase.