Effect of membrane tension on the electric field and dipole potential of lipid bilayer membrane

The dipole potential of lipid bilayer membrane controls the difference in permeability of the membrane to oppositely charged ions. We have combined molecular dynamics (MD) simulations and experimental studies to determine changes in electric field and electrostatic potential of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) lipid bilayer in response to applied membrane tension. MD simulations based on CHARMM36 force field showed that electrostatic potential of DOPC bilayer decreases by ~45mV in the physiologically relevant range of membrane tension values (0 to 15dyn/cm). The electrostatic field exhibits a peak (~0.8×10(9)V/m) near the water/lipid interface which shifts by 0.9? towards the bilayer center at 15dyn/cm. Maximum membrane tension of 15dyn/cm caused 6.4% increase in area per lipid, 4.7% decrease in bilayer thickness and 1.4% increase in the volume of the bilayer. Dipole-potential sensitive fluorescent probes were used to detect membrane tension induced changes in DOPC vesicles exposed to osmotic stress. Experiments confirmed that dipole potential of DOPC bilayer decreases at higher membrane tensions. These results are suggestive of a potentially new mechanosensing mechanism by which mechanically induced structural changes in the lipid bilayer membrane could modulate the function of membrane proteins by altering electrostatic interactions and energetics of protein conformational states.     

A continuum method for determining membrane protein insertion energies and the problem of charged residues

Continuum electrostatic approaches have been extremely successful at describing the charged nature of soluble proteins and how they interact with binding partners. However, it is unclear whether continuum methods can be used to quantitatively understand the energetics of membrane protein insertion and stability. Recent translation experiments suggest that the energy required to insert charged peptides into membranes is much smaller than predicted by present continuum theories. Atomistic simulations have pointed to bilayer inhomogeneity and membrane deformation around buried charged groups as two critical features that are neglected in simpler models. Here, we develop a fully continuum method that circumvents both of these shortcomings by using elasticity theory to determine the shape of the deformed membrane and then subsequently uses this shape to carry out continuum electrostatics calculations. Our method does an excellent job of quantitatively matching results from detailed molecular dynamics simulations at a tiny fraction of the computational cost. We expect that this method will be ideal for studying large membrane protein complexes.     

Membrane association and selectivity of the antimicrobial peptide NK‐2: a molecular dynamics simulation study

In an effort to better understand the initial mechanism of selectivity and membrane association of the synthetic antimicrobial peptide NK‐2, we have applied molecular dynamics simulation techniques to elucidate the interaction of the peptide with the membrane interfaces. A homogeneous dipalmitoylphosphatidylglycerol (DPPG) and a homogeneous dipalmitoylphosphatidylethanolamine (DPPE) bilayers were taken as model systems for the cytoplasmic bacterial and human erythrocyte membranes, respectively. The results of our simulations on DPPG and DPPE model membranes in the gel phase show that the binding of the peptide, which is considerably stronger for the negatively charged DPPG lipid bilayer than for the zwitterionic DPPE, is mostly governed by electrostatic interactions between negatively charged residues in the membrane and positively charged residues in the peptide. In addition, a characteristic distribution of positively charged residues along the helix facilitates a peptide orientation parallel to the membrane interface. Once the peptides reside close to the membrane surface of DPPG with the more hydrophobic side chains embedded into the membrane interface, the peptide initially disturbs the respective bilayer integrity by a decrease of the order parameter of lipid acyl chain close to the head group region, and by a slightly decrease in bilayer thickness. We found that the peptide retains a high content of helical structure on the zwitterionic membrane‐water interface, while the loss of α‐helicity is observed within a peptide adsorbed onto negatively charged lipid membranes. Copyright © 2009 European Peptide Society and John Wiley & Sons, Ltd.     

Water under the BAR

Many cellular processes require the generation of highly curved regions of cell membranes by interfacial membrane proteins. A number of such proteins are now known, and several mechanisms of curvature generation have been suggested, but so far a quantitative understanding of the importance of the various potential mechanisms remains elusive. Following previous theoretical work, we consider the electrostatic attraction that underlies the scaffold mechanism of membrane bending in the context of the N-BAR domain of amphiphysin. Analysis of atomistic molecular dynamics simulations reveals considerable water between the membrane and the positively charged concave face of the BAR, even when it is tightly bound to highly curved membranes. This results in significant screening of electrostatic interactions, suggesting that electrostatic attraction is not the main driving force behind curvature sensing, supporting recent experimental work. These results also emphasize the need for care when building coarse-grained models of protein-membrane interactions. These results are emphasized by simulations of oligomerized amphiphysin N-BARs at the atomistic and coarse-grained level. In the coarse-grained simulations, we find a strong dependence of the induced curvature on the dielectric screening.     

Membrane proteins: molecular dynamics simulations

Molecular dynamics simulations of membrane proteins are making rapid progress, because of new high-resolution structures, advances in computer hardware and atomistic simulation algorithms, and the recent introduction of coarse-grained models for membranes and proteins. In addition to several large ion channel simulations, recent studies have explored how individual amino acids interact with the bilayer or snorkel/anchor to the headgroup region, and it has been possible to calculate water/membrane partition free energies. This has resulted in a view of bilayers as being adaptive rather than purely hydrophobic solvents, with important implications, for example, for interaction between lipids and arginines in the charged S4 helix of voltage-gated ion channels. However, several studies indicate that the typical current simulations fall short of exhaustive sampling, and that even simple protein-membrane interactions require at least ca. 1mus to fully sample their dynamics. One new way this is being addressed is coarse-grained models that enable mesoscopic simulations on multi-mus scale. These have been used to model interactions, self-assembly and membrane perturbations induced by proteins. While they cannot replace all-atom simulations, they are a potentially useful technique for initial insertion, placement, and low-resolution refinement.     

Membrane bending is critical for the stability of voltage sensor segments in the membrane

The interaction between membrane proteins and the surrounding membrane is becoming increasingly appreciated for its role in regulating protein function, protein localization, and membrane morphology. In particular, recent studies have suggested that membrane deformation is needed to stably accommodate proteins harboring charged amino acids in their transmembrane (TM) region, as it is energetically prohibitive to bury charge in the hydrophobic core of the bilayer. Unfortunately, current computational methods are poorly equipped for describing such deformations, as atomistic simulations are often too short to observe large-scale membrane reorganization and most continuum approaches assume a flat membrane. Previously, we developed a method that overcomes these shortcomings by using elasticity theory to characterize equilibrium membrane distortions in the presence of a TM protein, while using traditional continuum electrostatic and nonpolar energy models to determine the energy of the protein in the membrane. Here, we linked the elastostatics, electrostatics, and nonpolar numeric solvers to permit the calculation of energies for nontrivial membrane deformations. We then coupled this procedure to a robust search algorithm that identifies optimal membrane shapes for a TM protein of arbitrary chemical composition. This advance now permits us to explore a host of biological phenomena that were beyond the scope of our original method. We show that the energy required to embed charged residues in the membrane can be highly nonadditive, and our model provides a simple mechanical explanation for this nonadditivity. Our results also predict that isolated voltage sensor segments do not insert into rigid membranes, but membrane bending dramatically stabilizes these proteins in the bilayer despite their high charge content. Additionally, we use the model to explore hydrophobic mismatch with regard to nonpolar peptides and mechanosensitive channels. Our method is in quantitative agreement with molecular dynamics simulations at a tiny fraction of the computational cost.     

Imaging the electrostatic potential of transmembrane channels: atomic probe microscopy of OmpF porin

The atomic force microscope (AFM) was used to image native OmpF porin and to detect the electrostatic potential generated by the protein. To this end the OmpF porin trimers from Escherichia coli was reproducibly imaged at a lateral resolution of approximately 0.5 nm and a vertical resolution of approximately 0.1 nm at variable electrolyte concentrations of the buffer solution. At low electrolyte concentrations the charged AFM probe not only contoured structural details of the membrane protein surface but also interacted with local electrostatic potentials. Differences measured between topographs recorded at variable ionic strength allowed mapping of the electrostatic potential of OmpF porin. The potential map acquired by AFM showed qualitative agreement with continuum electrostatic calculations based on the atomic OmpF porin embedded in a lipid bilayer at the same electrolyte concentrations. Numerical simulations of the experimental conditions showed the measurements to be reproduced quantitatively when the AFM probe was included in the calculations. This method opens a novel avenue to determine the electrostatic potential of native protein surfaces at a lateral resolution better than 1 nm and a vertical resolution of approximately 0.1 nm.     

Characterization of PROPPIN-Phosphoinositide Binding and Role of Loop 6CD in PROPPIN-Membrane Binding

PROPPINs (β-propellers that bind polyphosphoinositides) are a family of PtdIns3P- and PtdIns(3,5)P₂-binding proteins that play an important role in autophagy. We analyzed PROPPIN-membrane binding through isothermal titration calorimetry (ITC), stopped-flow measurements, mutagenesis studies, and molecular dynamics (MD) simulations. ITC measurements showed that the yeast PROPPIN family members Atg18, Atg21, and Hsv2 bind PtdIns3P and PtdIns(3,5)P₂ with high affinities in the nanomolar to low-micromolar range and have two phosphoinositide (PIP)-binding sites. Single PIP-binding site mutants have a 15- to 30-fold reduced affinity, which explains the requirement of two PIP-binding sites in PROPPINs. Hsv2 bound small unilamellar vesicles with a higher affinity than it bound large unilamellar vesicles in stopped-flow measurements. Thus, we conclude that PROPPIN membrane binding is curvature dependent. MD simulations revealed that loop 6CD is an anchor for membrane binding, as it is the region of the protein that inserts most deeply into the lipid bilayer. Mutagenesis studies showed that both hydrophobic and electrostatic interactions are required for membrane insertion of loop 6CD. We propose a model for PROPPIN-membrane binding in which PROPPINs are initially targeted to membranes through nonspecific electrostatic interactions and are then retained at the membrane through PIP binding.     

Evidence that Electrostatic Interactions between Vesicle-associated Membrane Protein 2 and Acidic Phospholipids May Modulate the Fusion of Transport Vesicles with the Plasma Membrane

The juxtamembrane domain of vesicle-associated membrane protein (VAMP) 2 (also known as synaptobrevin2) contains a conserved cluster of basic/hydrophobic residues that may play an important role in membrane fusion. Our measurements on peptides corresponding to this domain determine the electrostatic and hydrophobic energies by which this domain of VAMP2 could bind to the adjacent lipid bilayer in an insulin granule or other transport vesicle. Mutation of residues within the juxtamembrane domain that reduce the VAMP2 net positive charge, and thus its interaction with membranes, inhibits secretion of insulin granules in β cells. Increasing salt concentration in permeabilized cells, which reduces electrostatic interactions, also results in an inhibition of insulin secretion. Similarly, amphipathic weak bases (e.g., sphingosine) that reverse the negative electrostatic surface potential of a bilayer reverse membrane binding of the positively charged juxtamembrane domain of a reconstituted VAMP2 protein and inhibit membrane fusion. We propose a model in which the positively charged VAMP and syntaxin juxtamembrane regions facilitate fusion by bridging the negatively charged vesicle and plasma membrane leaflets.     

Changes in a phospholipid bilayer induced by the hydrolysis of a phospholipase A2 enzyme: a molecular dynamics simulation study

Phospholipase A2 (PLA2) enzymes are important in numerous physiological processes. Their function at lipid-water interfaces is also used as a biophysical model for protein-membrane interactions. These enzymes catalyze the hydrolysis of the sn-2 bonds of various phospholipids and the hydrolysis products are known to increase the activity of the enzymes. Here, we have applied molecular dynamics (MD) simulations to study the membrane properties in three compositionally different systems that relate to PLA2 enzyme action. One-nanosecond simulations were performed for a 1-palmitoyl-2-linoleoyl-sn-glycero-3-phosphatidylcholine (PLPC) bilayer and for two of its PLA2-hydrolyzed versions, i.e., bilayers consisting of lysophospholipids and of either free charged linoleate or free uncharged linoleic acid molecules. The results revealed loosening of the structure in the hydrolyzed bilayer due to increased mobility of the molecules in the direction normal to the bilayer. This loss of integrity due to the hydrolysis products is in accord with observations that not only the presence of hydrolysis products, but also a variety of other perturbations of the membrane may activate PLA2. Additionally, changes were observed in other structural parameters and in the electrostatic potential across the membrane-water interface. These changes are discussed in relation to the simulation methodology and the experimental observations of PLA2-hydrolyzed membranes.     

Molecular dynamics simulations of local anesthetic articaine in a lipid bilayer

In order to investigate structural and dynamical properties of local anesthetic articaine in a model lipid bilayer, a series of molecular dynamics simulations have been performed. Simulations were carried out for neutral and charged (protonated) forms of articaine inserted in fully hydrated dimyristoylphosphatidylcholine (DMPC) lipid bilayer. For comparison purpose, a fully hydrated DMPC bilayer without articaine was also simulated. The length of each simulation was 200 ns. Various properties of the lipid bilayer systems in the presence of both charged and uncharged forms of articaine taken at two different concentrations have been examined: membrane area per lipid, mass density distributions, order parameters, radial distribution functions, head group tilt, diffusion coefficients, electrostatic potential, etc, and compared with results of previous simulations of DMPC bilayer in the presence of lidocaine. It was shown that addition of both charged and neutral forms of articaine causes increase of the dipole electrostatic potential in the membrane interior.     

Molecular Dynamics Simulations of a DMPC Bilayer Using Nonadditive Interaction Models

We present a polarizable force field based on the charge-equilibration formalism for molecular dynamics simulations of phospholipid bilayers. We discuss refinement of headgroup dihedral potential parameters to reproduce ab initio conformational energies of dimethylphosphate calculated at the MP2/cc-pVTZ level of theory. We also address the refinement of electrostatic and Lennard-Jones (van der Waals) parameters to reproduce ab initio polarizabilities and water interaction energies of dimethylphosphate and tetramethylammonium. We present results of molecular dynamics simulations of a solvated dimyristoylphosphatidylcholine bilayer using this polarizable force field as well as the nonpolarizable, fixed-charge CHARMM27 and CHARMM27r force fields for comparison. Calculated atomic and electron-density profiles, deuterium order parameters, and headgroup orientations are found to be consistent with previous simulations and with experiment. Polarizable interaction models for solvent and lipid exhibit greater water penetration into the lipid interior; this is due to the variation of water molecular dipole moment from a bulk value of 2.6 Debye to a value of 1.9 Debye in the membrane interior. The reduction in the electrostatic component of the desolvation free-energy penalty allows for greater water density. The surface dipole potential predicted by the polarizable model is 0.95 V compared to the value of 0.8 V based on nonpolarizable force-field calculations. Effects of inclusion of explicit polarization are discussed in relation to water dipole moment and varying charge distributions. Dielectric permittivity profiles for polarizable and nonpolarizable interactions exhibit subtle differences arising from the nature of the individual component parameterizations; for the polarizable force field, we obtain a bulk dielectric permittivity of 79, whereas the nonpolarizable force field plateaus at 97 (the value for pure TIP3P water). In the membrane interior, both models predict unit permittivities, with the polarizable models contributing from one to two more units due to the optical dielectric (high-frequency dipole fluctuations). This contribution is a step toward the continuing development of a CHARMM (Chemistry at Harvard Molecular Mechanics) polarizable force field for simulations of biomacromolecular systems.     

How Does a Voltage Sensor Interact with a Lipid Bilayer? Simulations of a Potassium Channel Domain

The nature of voltage sensing by voltage-activated ion channels is a key problem in membrane protein structural biology. The way in which the voltage-sensor (VS) domain interacts with its membrane environment remains unclear. In particular, the known structures of Kv channels do not readily explain how a positively charged S4 helix is able to stably span a lipid bilayer. Extended (2 x 50 ns) molecular dynamics simulations of the high-resolution structure of the isolated VS domain from the archaebacterial potassium channel KvAP, embedded in zwitterionic and in anionic lipid bilayers, have been used to explore VS/lipid interactions at atomic resolution. The simulations reveal penetration of water into the center of the VS and bilayer. Furthermore, there is significant local deformation of the lipid bilayer by interactions between lipid phosphate groups and arginine side chains of S4. As a consequence of this, the electrostatic field is "focused" across the center of the bilayer.     

Reparameterization of all-atom dipalmitoylphosphatidylcholine lipid parameters enables simulation of fluid bilayers at zero tension

Molecular dynamics simulations of dipalmitoylphosphatidylcholine (DPPC) lipid bilayers using the CHARMM27 force field in the tensionless isothermal-isobaric (NPT) ensemble give highly ordered, gel-like bilayers with an area per lipid of approximately 48 A(2). To obtain fluid (L(alpha)) phase properties of DPPC bilayers represented by the CHARMM energy function in this ensemble, we reparameterized the atomic partial charges in the lipid headgroup and upper parts of the acyl chains. The new charges were determined from the electron structure using both the Mulliken method and the restricted electrostatic potential fitting method. We tested the derived charges in molecular dynamics simulations of a fully hydrated DPPC bilayer. Only the simulation with the new restricted electrostatic potential charges shows significant improvements compared with simulations using the original CHARMM27 force field resulting in an area per lipid of 60.4 +/- 0.1 A(2). Compared to the 48 A(2), the new value of 60.4 A(2) is in fair agreement with the experimental value of 64 A(2). In addition, the simulated order parameter profile and electron density profile are in satisfactory agreement with experimental data. Thus, the biologically more interesting fluid phase of DPPC bilayers can now be simulated in all-atom simulations in the NPT ensemble by employing our modified CHARMM27 force field.     

Membrane binding by MinD involves insertion of hydrophobic residues within the C-terminal amphipathic helix into the bilayer

MinD binds to phospholipid vesicles in the presence of ATP and is released by MinE, which stimulates the MinD ATPase. Membrane binding requires a short conserved C-terminal region, which has the potential to form an amphipathic helix. This finding has led to a model in which the binding of ATP regulates the formation or accessibility of this helix, which then embeds in the membrane bilayer. To test this model, we replaced each of the four hydrophobic residues within this potential helix with tryptophan or a charged residue. Introduction of a negatively charged amino acid decreased membrane binding of MinD and its ability to activate MinC. In contrast, mutants with tryptophan substitutions retained the ability to bind to the membrane and activate MinC. Fluorescence emission spectroscopy analysis of the tryptophan mutants F263W, L264W, and L267W confirmed that these tryptophan residues did insert into the hydrophobic interior of the bilayer. We conclude that membrane binding by MinD involves penetration of the hydrophobic residues within the C-terminal amphipathic helix into the hydrophobic interior of the bilayer.     

Exploring the binding dynamics of BAR proteins

We used a continuum model based on the Helfrich free energy to investigate the binding dynamics of a lipid bilayer to a BAR domain surface of a crescent-like shape of positive (e.g. I-BAR shape) or negative (e.g. F-BAR shape) intrinsic curvature. According to structural data, it has been suggested that negatively charged membrane lipids are bound to positively charged amino acids at the binding interface of BAR proteins, contributing a negative binding energy to the system free energy. In addition, the cone-like shape of negatively charged lipids on the inner side of a cell membrane might contribute a positive intrinsic curvature, facilitating the initial bending towards the crescent-like shape of the BAR domain. In the present study, we hypothesize that in the limit of a rigid BAR domain shape, the negative binding energy and the coupling between the intrinsic curvature of negatively charged lipids and the membrane curvature drive the bending of the membrane. To estimate the binding energy, the electric potential at the charged surface of a BAR domain was calculated using the Langevin-Bikerman equation. Results of numerical simulations reveal that the binding energy is important for the initial instability (i.e. bending of a membrane), while the coupling between the intrinsic shapes of lipids and membrane curvature could be crucial for the curvature-dependent aggregation of negatively charged lipids near the surface of the BAR domain. In the discussion, we suggest novel experiments using patch clamp techniques to analyze the binding dynamics of BAR proteins, as well as the possible role of BAR proteins in the fusion pore stability of exovesicles.     

NaCl interactions with phosphatidylcholine bilayers do not alter membrane structure but induce long-range ordering of ions and water

It is generally accepted that ions interact directly with lipids in biological membranes. Decades of biophysical studies on pure lipid bilayer systems have shown that only certain types of ions, most significantly large anions and multivalent cations, can fundamentally alter the structure and dynamics of lipid bilayers. It has long been accepted that at physiological concentrations NaCl ions do not alter the physical behavior or structure of bilayers composed solely of zwitterionic phosphatidylcholine (PC) lipids. Recent X-ray scattering experiments have reaffirmed this dogma, showing that below 1 M concentration, NaCl does not significantly alter bilayer structure. However, despite this history, there is an ongoing controversy within the molecular dynamics (MD) simulation community regarding NaCl/PC interactions. In particular, the CHARMM and GROMOS force fields show dramatically different behavior, including the effect on bilayer structure, surface potential, and the ability to form stable, coordinated ion–lipid complexes. Here, using long-timescale, constant-pressure simulations under the newest version of the CHARMM force field, we find that Na⁺ and Cl⁻ associate with PC head groups in a POPC bilayer with approximately equal, though weak, affinity, and that the salt has a negligible effect on bilayer structure, consistent with earlier CHARMM results and more recent X-ray data. The results suggest that interpretation of simulations where lipids interact with charged groups of any sort, including charged proteins, must be carefully scrutinized.     

Phosphatidylinositol 3-phosphate induces the membrane penetration of the FYVE domains of Vps27p and Hrs

The FYVE domain mediates the recruitment of proteins involved in membrane trafficking and cell signaling to phosphatidylinositol 3-phosphate (PtdIns(3)P)-containing membranes. To elucidate the mechanism by which the FYVE domain interacts with PtdIns(3)P-containing membranes, we measured the membrane binding of the FYVE domains of yeast Vps27p and Drosophila hepatocyte growth factor-regulated tyrosine kinase substrate and their mutants by surface plasmon resonance and monolayer penetration analyses. These measurements as well as electrostatic potential calculation show that PtdIns(3)P specifically induces the membrane penetration of the FYVE domains and increases their membrane residence time by decreasing the positive charge surrounding the hydrophobic tip of the domain and causing local conformational changes. Mutations of hydrophobic residues located close to the PtdIns(3)P-binding pocket or an Arg residue directly involved in PtdIns(3)P binding abrogated the penetration of the FYVE domains into the monolayer, the packing density of which is comparable with that of biological membranes and large unilamellar vesicles. Based on these results, we propose a mechanism of the membrane binding of the FYVE domain in which the domain first binds to the PtdIns(3)P-containing membrane by specific PtdIns(3)P binding and nonspecific electrostatic interactions, which is then followed by the PtdIns(3)P-induced partial membrane penetration of the domain.     

Biomembrane simulations of 12 lipid types using the general amber force field in a tensionless ensemble

The AMBER family of force fields is one of the most commonly used alternatives to describe proteins and drug-like molecules in molecular dynamics simulations. However, the absence of a specific set of parameters for lipids has been limiting the widespread application of this force field in biomembrane simulations, including membrane protein simulations and drug-membrane simulations. Here, we report the systematic parameterization of 12 common lipid types consistent with the General Amber Force Field (GAFF), with charge-parameters determined with RESP at the HF/6–31G(d) level of theory, to be consistent with AMBER. The accuracy of the scheme was evaluated by comparing predicted and experimental values for structural lipid properties in MD simulations in an NPT ensemble with explicit solvent in 100:100 bilayer systems. Globally, a consistent agreement with experimental reference data on membrane structures was achieved for some lipid types when using the typical MD conditions normally employed when handling membrane proteins and drug-membrane simulations (a tensionless NPT ensemble, 310?K), without the application of any of the constraints often used in other biomembrane simulations (such as the surface tension and the total simulation box area). The present set of parameters and the universal approach used in the parameterization of all the lipid types described here, as well as the consistency with the AMBER force field family, together with the tensionless NPT ensemble used, opens the door to systematic studies combining lipid components with small drug-like molecules or membrane proteins and show the potential of GAFF in dealing with biomembranes.