Many-body polarization effects and the membrane dipole potential

Free energetics and the role of water in the permeation of methyl guanidinium across the bilayer-water interface: insights from molecular dynamics simulations using charge equilibration potentials

Combining umbrella sampling molecular dynamics (MD) simulations, the weighted histogram analysis method (WHAM) for unbiasing probabilities, and polarizable charge equilibration force fields, we compute the potential of mean force for the reversible transfer of methyl guanidinium from bulk solution to the center of a model DPPC bilayer. A 5 kcal/mol minimum in the potential of mean force profile for membrane permeation suggests that the analogue will preferentially reside in the headgroup region of the lipid, qualitatively in agreement with previously published results. We find the potential of mean force for permeation to be approximately 28 kcal/mol (relative to the minimum in the headgroups), within the range of values reported for similar types of simulations using fixed-charge force fields. From analysis of the lipid structure, we find that the lipid deformation leads to a substantial destabilizing contribution to the free energy of the methyl guanidinium as it resides in the bilayer center, though this deformation allows more efficient stabilization by water defects and transient pores. Water in the bilayer core stabilizes the charged residue. The role of water in stabilizing or destabilizing the solute as it crosses the bilayer depends on bulk electrolyte concentration. In 1 M KCl solution, the water contribution to the potential of mean force is stabilizing over the entire range of the permeation coordinate, with the sole destabilizing force originating from the anionic species in solution. Conversely, methyl guanidinium experiences net destabilization from water in the absence of electrolyte. The difference in solvent contributions to permeation free energy is traced to a local effect arising from differences in water density in the bilayer-water solution interface, thus leading to starkly opposite net forces on the permeant. The origin of the local water density differential rests with the penetration of hydrated chloride anions into the solution-bilayer interface. Finally, water permeation into the bilayer is required for the deformation of individual lipid molecules and permeation of ions into the membrane. From simulations where water is first excluded from the bilayer center where methyl guanidinium is restrained and then, after equilibration, allowed to enter the bilayer, we find that in the absence of any water defects/permeation into the bilayer, the lipid headgroups do not follow the methyl guanidinium. Only when water enters the bilayer do we see deformation of individual lipid molecules to associate with the amino acid analogue at bilayer center.

Electrostatics interactions in classical simulations

Electrostatic interactions are crucial for both the accuracy and performance of atomistic biomolecular simulations. In this chapter we review well-established methods and current developments aiming at efficiency and accuracy. Specifically, we review the classical Ewald summations, particle-particle particle-method particle-method Ewald algorithms, multigrid, fast multipole, and local methods. We also highlight some recent developments targeting more accurate, yet classical, representation of the molecular charge distribution.

Polarizable force fields

This chapter provides an overview of the most common methods for including an explicit description of electronic polarization in molecular mechanics force fields: the induced point dipole, shell, and fluctuating charge models. The importance of including polarization effects in biomolecular simulations is discussed, and some of the most important achievements in the development of polarizable biomolecular force fields to date are highlighted.

Derivation and systematic validation of a refined all-atom force field for phosphatidylcholine lipids

An all-atomistic force field (FF) has been developed for fully saturated phospholipids. The parametrization has been largely based on high-level ab initio calculations in order to keep the empirical input to a minimum. Parameters for the lipid chains have been developed based on knowledge about bulk alkane liquids, for which thermodynamic and dynamic data are excellently reproduced. The FFs ability to simulate lipid bilayers in the liquid crystalline phase in a tensionless ensemble was tested in simulations of three lipids: 1,2-diauroyl-sn-glycero-3-phospocholine (DLPC), 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), and 1,2-dipalmitoyl-sn-glycero-3-phospcholine (DPPC). Computed areas and volumes per lipid, and three different kinds of bilayer thicknesses, have been investigated. Most importantly NMR order parameters and scattering form factors agree in an excellent manner with experimental data under a range of temperatures. Further, the compatibility with the AMBER FF for biomolecules as well as the ability to simulate bilayers in gel phase was demonstrated. Overall, the FF presented here provides the important balance between the hydrophilic and hydrophobic forces present in lipid bilayers and therefore can be used for more complicated studies of realistic biological membranes with protein insertions.

Measurements and implications of the membrane dipole potential

There are three kinds of membrane potentials: the surface potentials, resulting from the accumulation of charges at the membrane surfaces; the transmembrane potential, determined by imbalance of charge in the aqueous solutions; and the dipole potential, a membrane-internal potential from the dipolar components of the phospholipids and interface water. The absolute value of the dipole potential has been very difficult to measure, although its value has been estimated to be in the range of 200-1,000 mV from ion translocation rates (determined by the planar lipid bilayer method), the surface potential of lipid monolayers (determined by the lipid monolayer method), molecular-dynamics calculations, and electron scattering using cryoelectron microscopy (cryo-EM). Spectroscopy methods have also been used to monitor the dipole potential changes on the basis of the observed fluorescence changes of voltage-sensitive probes. The dipole potential accounts for the much larger permeability of a bare phospholipid membrane to anions than cations and affects the conformation and function of membrane proteins.

Determination of membrane-insertion free energies by molecular dynamics simulations

The accurate prediction of membrane-insertion probability for arbitrary protein sequences is a critical challenge to identifying membrane proteins and determining their folded structures. Although algorithms based on sequence statistics have had moderate success, a complete understanding of the energetic factors that drive the insertion of membrane proteins is essential to thoroughly meeting this challenge. In the last few years, numerous attempts to define a free-energy scale for amino-acid insertion have been made, yet disagreement between most experimental and theoretical scales persists. However, for a recently resolved water-to-bilayer scale, it is found that molecular dynamics simulations that carefully mimic the conditions of the experiment can reproduce experimental free energies, even when using the same force field as previous computational studies that were cited as evidence of this disagreement. Therefore, it is suggested that experimental and simulation-based scales can both be accurate and that discrepancies stem from disparities in the microscopic processes being considered rather than methodological errors. Furthermore, these disparities make the development of a single universally applicable membrane-insertion free energy scale difficult.

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.

The Role of Atomic Polarization in the Thermodynamics of Chloroform Partitioning to Lipid Bilayers

In spite of extensive research and use in medical practice, the precise molecular mechanism of volatile anesthetic action remains unknown. The distribution of anesthetics within lipid bilayers and potential targeting to membrane proteins is thought to be central to therapeutic function. Therefore, obtaining a molecular level understanding of volatile anesthetic partitioning into lipid bilayers is of vital importance to modern pharmacology. In this study we investigate the partitioning of the prototypical anesthetic, chloroform, into lipid bilayers and different organic solvents using molecular dynamics simulations with potential models ranging from simplified coarse-grained MARTINI to additive and polarizable CHARMM all-atom force fields. Many volatile anesthetics display significant inducible dipole moments, which correlate with their potency, yet the exact role of molecular polarizability in their stabilization within lipid bilayers remains unknown. We observe that explicit treatment of atomic polarizability makes it possible to accurately reproduce solvation free energies in solvents with different polarities, allowing for quantitative studies in heterogeneous molecular distributions, such as lipid bilayers. We calculate the free energy profiles for chloroform crossing lipid bilayers to reveal a role of polarizability in modulating chloroform partitioning thermodynamics via the chloroform-induced dipole moment and highlight competitive binding to the membrane core and toward the glycerol backbone that may have significant implications for understanding anesthetic action.

Molecular dynamics simulation of hydrated DPPC monolayers using charge equilibration force fields

We present results of molecular dynamics simulations of a model DPPC-water monolayer using charge equilibration (CHEQ) force fields, which explicitly account for electronic polarization in a classical treatment of intermolecular interactions. The surface pressure, determined as the difference between the monolayer and pure water surface tensions at 323 K, is predicted to be 22.92 ±1.29 dyne/cm, just slightly below the broad range of experimental values reported for this system. The surface tension for the DPPC-water monolayer is predicted to be 42.35 ±1.16 dyne/cm, in close agreement with the experimentally determined value of 40.9 dyne/cm. This surface tension is also consistent with the value obtained from DPPC monolayer simulations using state-of-the-art nonpolarizable force fields. The current results of simulations predict a monolayer-water potential difference relative to the pure water-air interface of 0.64 ±0.02 Volts, an improved prediction compared to the fixed-charge CHARMM27 force field, yet still overestimating the experimental range of 0.3 to 0.45 Volts. As the charge equilibration model is a purely charge-based model for polarization, the current results suggest that explicitly modeled polarization effects can offer improvements in describing interfacial electrostatics in such systems.

Recent Developments and Applications of the CHARMM force fields

Empirical force fields commonly used to describe the condensed phase properties of complex systems such as biological macromolecules are continuously being updated. Improvements in quantum mechanical (QM) methods used to generate target data, availability of new experimental target data, incorporation of new classes of compounds and new theoretical developments (eg. polarizable methods) make force-field development a dynamic domain of research. Accordingly, a number of improvements and extensions of the CHARMM force fields have occurred over the years. The objective of the present review is to provide an up-to-date overview of the CHARMM force fields. A limited presentation on the historical aspects of force fields will be given, including underlying methodologies and principles, along with a brief description of the strategies used for parameter development. This is followed by information on the CHARMM additive and polarizable force fields, including examples of recent applications of those force fields.

Influence of the membrane dipole potential on peptide binding to lipid bilayers

The implicit membrane model IMM1 is extended to include the membrane dipole potential and applied to molecular dynamics simulations of the helical peptides alamethicin, WALP23, influenza hemagglutinin fusion peptide, HIV fusion peptide, magainin, and the pre-sequence of cytochrome c oxidase subunit IV (p25). The results show that the orientation of the peptides in the membrane can be influenced by the dipole potential. The binding affinity of all peptides except for the hemagglutinin fusion peptide decreases upon increase of the dipole potential. The changes in both orientation and binding affinity are explained by the interaction of the dipole potential with the helix backbone dipole and ionic side-chains. In general, peptides that tend to insert the N-terminus in the membrane and/or have positively charged side chains will lose binding affinity upon increase of the dipole potential.

The polarizable point dipoles method with electrostatic damping: implementation on a model system

Recently, the use of polarizable force fields in Molecular Dynamics simulations has been gaining importance, since they allow a better description of heterogeneous systems compared to simple point charges force fields. Among the various techniques developed in the last years the one based on polarizable point dipoles represents one of the most used. In this paper, we review the basic technical issues of the method, illustrating the way to implement intramolecular and intermolecular damping of the electrostatic interactions, either with and without the Ewald summation method. We also show how to reduce the computational overhead for evaluating the dipoles, introducing to the state-of-the-art methods: the extended Lagrangian method and the always stable predictor corrector method. Finally we discuss the importance of screening the electrostatic interactions at short range, defending this technique against simpler approximations usually made. We compare results of density functional theory and classical force field-based Molecular Dynamics simulations of chloride in water.

Development of the CHARMM Force Field for Lipids

The development of the CHARMM additive all-atom lipid force field (FF) is traced from the early 1990's to the most recent version (C36) published in 2010. Though simulations with early versions yielded useful results, they failed to reproduce two important quantities: a zero surface tension at the experimental bilayer surface area, and the signature splitting of the deuterium order parameters in the glycerol and upper chain carbons. Systematic optimization of parameters based on high level quantum mechanical data and free energy simulations have resolved these issues, and bilayers with a wide range of lipids can be simulated in tensionless ensembles using C36. Issues associated with other all-atom lipid FFs, success and limitations in the C36 FF and ongoing developments are also discussed.

Disruption and formation of surface salt bridges are coupled to DNA binding by the integration host factor: a computational analysis

Revealing the thermodynamic driving force of protein-DNA interactions is crucial to the understanding of factors that dictate the properties and function of protein-DNA complexes. For the binding of DNA to DNA-wrapping proteins, such as the integration host factor (IHF), Record and co-workers proposed that the disruption of a large number of preexisting salt bridges is coupled with the binding process [Holbrook, J. A., et al. (2001) J. Mol. Biol. 310, 379]. To test this proposal, we have conducted explicit solvent MD simulations (multiple ∼25-50 ns trajectories for each salt concentration) to examine the behavior of charged residues in IHF, especially concerning their ability to form salt bridges at different salt concentrations. Of the 17 cationic residues noted by Record and co-workers, most are engaged in salt bridge interactions for a significant portion of the trajectories, especially in the absence of salt. This observation suggests that, from a structural point of view, their proposal is plausible. However, the complex behaviors of charged residues observed in the MD simulations also suggest that the unusual thermodynamic characteristics of IHF-DNA binding likely arise from the interplay between complex dynamics of charged residues both in and beyond the DNA binding site. Moreover, a comparison of MD simulations at different salt concentrations suggests that the strong dependence of the IHF-DNA binding enthalpy on salt concentration may not be due to a significant decrease in the number of stable salt bridges in apo IHF at high salt concentrations. In addition to the Hofmeister effects quantified in more recent studies of IHF-DNA binding, we recommend consideration of the variation of the enthalpy change of salt bridge disruption at different salt concentrations. Finally, the simulation study presented here explicitly highlights the fact that the electrostatic properties of DNA-binding proteins can be rather different in the apo and DNA-bound states, which has important implications for the design of robust methods for predicting DNA binding sites in proteins.

Electrostatics of deformable lipid membranes

It was recently demonstrated that significant local deformations of biological membranes take place due to the fields of charged peptides and ions, challenging the standard model of membrane electrostatics. The ability of ions to retain their immediate hydration environment, combined with the lack of sensitivity of permeability to ion type or even ion pairs, led us to question the extent to which hydration energetics and electrostatics control membrane ion permeation. Using the arginine analog methyl-guanidinium as a test case, we find that although hydrocarbon electronic polarizability causes dramatic changes in ion solvation free energy, as well as a significant change (approximately 0.4 V) in the membrane dipole potential, little change in membrane permeation energetics occurs. We attribute this to compensation of solvation terms from polar and polarizable nonpolar components within the membrane, and explain why the dipole potential is not fully sensed in terms of the locally deformed bilayer interface. Our descriptions provide a deeper understanding of the translocation process and allow predictions for poly-ions, ion pairs, charged lipids, and lipid flip-flop. We also report simulations of large hydrophobic-ion-like membrane defects and the ionophore valinomycin, which exhibit little membrane deformation, as well as hydrophilic defects and the ion channel gramicidin A, to provide parallels to membranes deformed by unassisted ion permeation.

Halothane solvation in water and organic solvents from molecular simulations with new polarizable potential function

The partitioning of a substrate from one phase into another is a complex process with widespread applications: from chemical technology to the pharmaceutical industry. One particularly well-known and well-studied example is 2-bromo-2-chloro-1,1,1-trifluoroethane (halothane) trafficking through the lipid bilayer. Halothane is a model volatile anesthetic known to impact functions of model lipid bilayers, altering the structure and thickness upon its partitioning from the bulk phase. A number of theoretical and experimental investigations suggest the importance of electronic polarizability, determining a preference for halothane to partition in the interfacial systems as in lipid bilayers or binary solvents. The recently published protocol for the development of polarizable force fields based on the classical Drude model has provided fresh impetus to efforts directed at understanding the molecular principles governing complex thermodynamics of the hydrophobic hydration. Here, molecular simulations were combined with free energy simulations to study solvation of halothane in polarizable water and methanol. The absolute free energy of halothane solvation in different solvents (water, methanol, and n-hexane) has been evaluated for additive and polarizable models. It was found that both additive and polarizable models provide an adequate description of the halothane solvation in high-dielectric (polar) solvents such as water, but explicit accounting for electronic polarization is imperative for a correct description of the solvation thermodynamics in nonpolar systems. To study halothane dynamics in binary mixtures, all-atom molecular dynamics (MD) simulations for halothane-methanol mixtures in a wide range of concentrations were performed alongside an analysis of structural organization, dynamics, and thermodynamic properties to dissect the molecular determinants of the halothane solvation in polar and amphiphilic liquids such as methanol. Additionally, a theoretical test of the hypothesis on the weak hydrogen bonding of halothane and methanol in the condensed phase is provided, which was presented on the basis of spectroscopic analysis of the C-H vibrations in different gas-phase complexes. The simulations performed in the condensed phase suggest that hydrophobic interactions between halothane and methanol play a dominant role in preferential solvation.

Simulating Monovalent and Divalent Ions in Aqueous Solution Using a Drude Polarizable Force Field

An accurate representation of ion solvation in aqueous solution is critical for meaningful computer simulations of a broad range of physical and biological processes. Polarizable models based on classical Drude oscillators are introduced and parametrized for a large set of monoatomic ions including cations of the alkali metals (Li(+), Na(+), K(+), Rb(+) and Cs(+)) and alkaline earth elements (Mg(2+), Ca(2+), Sr(2+) and Ba(2+)) along with Zn(2+) and halide anions (F(-), Cl(-), Br(-) and I(-)). The models are parameterized, in conjunction with the polarizable SWM4-NDP water model [Lamoureux et al., Chem. Phys. Lett. 418, 245 (2006)], to be consistent with a wide assortment of experimentally measured aqueous bulk thermodynamic properties and the energetics of small ion-water clusters. Structural and dynamic properties of the resulting ion models in aqueous solutions at infinite dilution are presented.

Polarizable intermolecular potentials for water and benzene interacting with halide and metal ions

A complete derivation of polarizable intermolecular potentials based on high-level, gas-phase quantum-mechanical calculations is proposed. The importance of appreciable accuracy together with inherent simplicity represents a significant endeavor when enhancement of existing force fields for biological systems is sought. Toward this end, symmetry-adapted perturbation theory (SAPT) can provide an expansion of the total interaction energy into physically meaningful e.g. electrostatic, induction and van der Waals terms. Each contribution can be readily compared with its counterpart in classical force fields. Since the complexity of the different intermolecular terms cannot be fully embraced using a minimalist description, it is necessary to resort to polyvalent expressions capable of encapsulating overlooked contributions from the quantum-mechanical expansion. This choice results in consistent force field components that reflect the underlying physical principles of the phenomena. This simplified potential energy function is detailed and definitive guidelines are drawn. As a proof of concept, the methodology is illustrated through a series of test cases that include the interaction of water and benzene with halide and metal ions. In each case considered, the total energy is reproduced accurately over a range of biologically relevant distances.

Polarization effects in molecular mechanical force fields

The focus here is on incorporating electronic polarization into classical molecular mechanical force fields used for macromolecular simulations. First, we briefly examine currently used molecular mechanical force fields and the current status of intermolecular forces as viewed by quantum mechanical approaches. Next, we demonstrate how some components of quantum mechanical energy are effectively incorporated into classical molecular mechanical force fields. Finally, we assess the modeling methods of one such energy component-polarization energy-and present an overview of polarizable force fields and their current applications. Incorporating polarization effects into current force fields paves the way to developing potentially more accurate, though more complex, parameterizations that can be used for more realistic molecular simulations.

Charge equilibration force fields for lipid environments: applications to fully hydrated DPPC bilayers and DMPC-embedded gramicidin A

Polarizable force fields for lipid and solvent environments are used for molecular dynamics simulations of a fully hydrated dipalmitoylphosphatidylcholine (DPPC) bilayer and gramicidin A (gA) dimer embedded in a dimyristoylphosphatidylcholine (DMPC) bilayer. The lipid bilayer is modeled using the CHARMM charge equilibration (CHEQ) polarizable force field for lipids and the TIP4P-FQ force field to represent solvent. For the DPPC bilayer system, results are compared to the same system simulated using the nonpolarizable CHARMM27r (C27r) force field and TIP3P water. Calculated atomic and electron density profiles, head group orientations as measured by the phosphorus-nitrogen vector orientation, and deuterium order parameters are found to be consistent with previous simulations and with experiment. The CHEQ model exhibits greater water penetration into the bilayer interior, as demonstrated by the potential of mean force calculated from the water density profile. This is a result of the variation of the water molecular dipole from 2.55 D in the bulk to 1.88 D in the interior. We discuss this finding in the context of previous studies (both simulation and experiment) that have investigated the extent of penetration of water into DPPC bilayers. We also discuss the effects of including explicit polarization on the water dipole moment variation as a function of distance from the bilayer. We show distributions of atomic charges over the course of the simulation since the CHEQ model allows the charges to fluctuate. We have calculated the interfacial dipole potential, which the CHEQ model predicts to be 0.95 V compared to 0.86 V as predicted by the C27r model. We also discuss dielectric permittivity profiles and the differences arising between the two models. We obtain bulk values of 72.77 for the CHEQ model (TIP4P-FQ water) and 91.22 for C27r (TIP3P), and values approaching unity in the membrane interior. Finally, we present results of simulations of gA embedded in a DMPC bilayer using the CHEQ model and discuss structural properties.

Nanoscopic description of biomembrane electrostatics: results of molecular dynamics simulations and fluorescence probing

Electrostatic fields generated on and inside biological membranes are recognized to play a fundamental role in key processes of cell functioning. Their understanding requires an adequate description on the level of elementary charges and the reconstruction of electrostatic potentials by integration over all elementary interactions. Out of all the available research tools, only molecular dynamics simulations are capable of this, extending from the atomic to the mesoscopic level of description on the required time and space scale. A complementary approach is that offered by molecular probe methods, with the application of electrochromic dyes. Highly sensitive to intermolecular interactions, they generate integrated signals arising from electric fields produced by elementary charges at the sites of their location. This review is an attempt to provide a critical analysis of these two approaches and their present and potential applications. The results obtained by both methods are consistent in that they both show an extremely complex profile of the electric field in the membrane. The nanoscopic view, with two-dimensional averaging over the bilayer plane and formal separation of the electrostatic potential into surface (Psi(s)), dipole (Psi(d)) and transmembrane (Psi(t)) potentials, is constructive in the analysis of different functional properties of membranes.

Molecular dynamics simulations of membrane channels and transporters

Membrane transport constitutes one of the most fundamental processes in all living cells with proteins as major players. Proteins as channels provide highly selective diffusive pathways gated by environmental factors, and as transporters furnish directed, energetically uphill transport consuming energy. X-ray crystallography of channels and transporters furnishes a rapidly growing number of atomic resolution structures, permitting molecular dynamics (MD) simulations to reveal the physical mechanisms underlying channel and transporter function. Ever increasing computational power today permits simulations stretching up to 1 micros, that is, to physiologically relevant time scales. Membrane protein simulations presently focus on ion channels, on aquaporins, on protein-conducting channels, as well as on various transporters. In this review we summarize recent developments in this rapidly evolving field.