Condensed-phase properties of n-alkyl-amines from molecular dynamics simulations using charge equilibration force fields

Microscopic Structure of Neat Linear Alkylamine Liquids: An X-Ray Scattering and Computer Simulation Study

Linear amines, from propylamine to nonylamine, are studied under ambient conditions by X-ray scattering and molecular dynamics simulations of various force field models. The major finding is that the prepeak in alkylamines is about 1 order of magnitude weaker than that in alkanols, hence suggesting much weaker hydrogen bonding-induced clustering of the amine groups than for the hydroxyl groups. Computer simulation studies reveal that the OPLS-UA model reproduces the prepeak, but with larger amplitudes, while the GROMOS-UA and CHARMM-AA force fields show almost no prepeak. Simulations of all models show the existence of hydrogen-bonded clusters, equally confirmed by the prominent prepeak of the structure factor between the nitrogen atoms. The hydrogen bond strength, as modeled by the Coulomb association in classical force field models, is about the same order of magnitude for both systems. Then, one may ask what is the origin of the weaker prepeak in alkylamines? Simulation data reveal that the existence of the prepeak is controlled through the cancellation of the positive contributions from the charged group correlations by the negative contributions from the cross charged-uncharged correlations. The C2v symmetry of the amine headgroup hinders clustering, which favors cross correlations with the tail atoms. This is opposite to alkanols where the symmetry of the hydroxyl headgroup favors clustering and hinders cross correlations with the alkyl tail. This competition between charged and uncharged atomic groups appears as a general mechanism to explain the existence of scattering prepeaks, including their position and amplitude.

Unifying Charge-Flow Polarization Models

We show that several models where electric polarization in molecular systems is modeled by charge-flow between atoms can all be considered as different manifestations of a general underlying mathematical structure. The models can be classified according to whether they employ atomic or bond parameters and whether they employ atom/bond hardness or softness. We show that an ab initio calculated charge response kernel can be considered as the inverse screened Coulombic matrix projected onto the zero-charge subspace, and this may provide a method for deriving charge screening functions to be used in force fields. The analysis suggests that some models contain redundancies, and we argue that a parameterization of charge-flow models in terms of bond softness is preferable as it depends on local quantities and decay to zero upon bond dissociation, while bond hardness depends on global quantities and increases toward infinity upon bond dissociation.

Direct computation of parameters for accurate polarizable force fields

We present an improved electronic linear response model to incorporate polarization and charge-transfer effects in polarizable force fields. This model is a generalization of the Atom-Condensed Kohn-Sham Density Functional Theory (DFT), approximated to second order (ACKS2): it can now be defined with any underlying variational theory (next to KS-DFT) and it can include atomic multipoles and off-center basis functions. Parameters in this model are computed efficiently as expectation values of an electronic wavefunction, obviating the need for their calibration, regularization, and manual tuning. In the limit of a complete density and potential basis set in the ACKS2 model, the linear response properties of the underlying theory for a given molecular geometry are reproduced exactly. A numerical validation with a test set of 110 molecules shows that very accurate models can already be obtained with fluctuating charges and dipoles. These features greatly facilitate the development of polarizable force fields.

Charge equilibration force fields for molecular dynamics simulations of lipids, bilayers, and integral membrane protein systems

With the continuing advances in computational hardware and novel force fields constructed using quantum mechanics, the outlook for non-additive force fields is promising. Our work in the past several years has demonstrated the utility of polarizable force fields, those based on the charge equilibration formalism, for a broad range of physical and biophysical systems. We have constructed and applied polarizable force fields for lipids and lipid bilayers. In this review of our recent work, we discuss the formalism we have adopted for implementing the charge equilibration (CHEQ) method for lipid molecules. We discuss the methodology, related issues, and briefly discuss results from recent applications of such force fields. Application areas include DPPC-water monolayers, potassium ion permeation free energetics in the gramicidin A bacterial channel, and free energetics of permeation of charged amino acid analogs across the water-bilayer interface. This article is part of a Special Issue entitled: Membrane protein structure and function.

Incorporating Phase-Dependent Polarizability in Non-Additive Electrostatic Models for Molecular Dynamics Simulations of the Aqueous Liquid-Vapor Interface

We discuss a new classical water force field that explicitly accounts for differences in polarizability between liquid and vapor phases. The TIP4P-QDP (4-point transferable intermolecular potential with charge dependent-polarizability) force field is a modification of the original TIP4P-FQ fluctuating charge water force field of Rick et al.(1) that self-consistently adjusts its atomic hardness parameters via a scaling function dependent on the M-site charge. The electronegativity (χ) parameters are also scaled in order to reproduce condensed-phase dipole moments of comparable magnitude to TIP4P-FQ. TIP4P-QDP is parameterized to reproduce experimental gas-phase and select condensed-phase properties. The TIP4P-QDP water model possesses a gas phase polarizability of 1.40 Å(3) and gas-phase dipole moment of 1.85 Debye, in excellent agreement with experiment and high-level ab initio predictions. The liquid density of TIP4P-QDP is 0.9954(±0.0002) g/cm(3) at 298 K and 1 atmosphere, and the enthalpy of vaporization is 10.55(±0.12) kcal/mol. Other condensed-phase properties such as the isobaric heat capacity, isothermal compressibility, and diffusion constant are also calculated within reasonable accuracy of experiment and consistent with predictions of other current state-of-the-art water force fields. The average molecular dipole moment of TIP4P-QDP in the condensed phase is 2.641(±0.001) Debye, approximately 0.02 Debye higher than TIP4P-FQ and within the range of values currently surmised for the bulk liquid. The dielectric constant, ε = 85.8 ± 1.0, is 10% higher than experiment. This is reasoned to be due to the increase in the condensed phase dipole moment over TIP4P-FQ, which estimates ε remarkably well. Radial distribution functions for TIP4P-QDP and TIP4P-FQ show similar features, with TIP4P-QDP showing slightly reduced peak heights and subtle shifts towards larger distance interactions. Since the greatest effects of the phase-dependent polarizability are anticipated in regions with both liquid and vapor character, interfacial simulations of TIP4P-QDP were performed and compared to TIP4P-FQ, a static polarizability analog. Despite similar features in density profiles such as the position of the GDS and interfacial width, enhanced dipole moments are observed for the TIP4P-QDP interface and onset of the vapor phase. Water orientational profiles show an increased preference (over TIP4P-FQ) in the orientation of the permanent dipole vector of the molecule within the interface; an enhanced z-induced dipole moment directly results from this preference. Hydrogen bond formation is lower, on average, in the bulk for TIP4P-QDP than TIP4P-FQ. However, the average number of hydrogen bonds formed by TIP4P-QDP in the interface exceeds that of TIP4P-FQ, and observed hydrogen bond networks extend further into the gaseous region. The TIP4P-QDP interfacial potential, calculated to be -11.98(±0.08) kcal/mol, is less favorable than that for TIP4P-FQ by approximately 2% as a result of a diminished quadrupole contribution. Surface tension is calculated within a 1.3% reduction from the experimental value. Results reported demonstrate TIP4P-QDP as a model comparable to the popular TIP4P-FQ while accounting for a physical effect previously neglected by other water models. Further refinements to this model, as well as future applications are discussed.