Effect of NaCl and KCl on phosphatidylcholine and phosphatidylethanolamine lipid membranes: insight from atomic-scale simulations for understanding salt-induced effects in the plasma membrane

Diffusive nature of xenon anesthetic changes properties of a lipid bilayer: molecular dynamics simulations

Effects of general anesthesia can be controllable by the ambient pressure. We perform molecular dynamics simulations for a 1-palmitoyl-2-oleoyl phosphatidylethanolamine lipid bilayer with or without xenon molecules by changing the pressure to elucidate the mechanism of the pressure reversal of general anesthesia. According to the diffusive nature of xenon molecules in the lipid bilayer, a decrease in the orientational order of the lipid tails, an increase in the area and volume per lipid molecule, and an increase in the diffusivity of lipid molecules are observed. We show that the properties of the lipid bilayer with xenon molecules at high pressure come close to those without xenon molecules at 0.1 MPa. Furthermore, we find that xenon molecules are concentrated in the middle of the lipid bilayer at high pressures by the pushing effect and that the diffusivity of xenon molecules is suppressed. These results suggest that the pressure reversal originates from a jamming and suppression of the diffusivity of xenon molecules in lipid bilayers.

Biophysics of lipid bilayers containing oxidatively modified phospholipids: insights from fluorescence and EPR experiments and from MD simulations

This review focuses on the influence of oxidized phosphatidylcholines (oxPCs) on the biophysical properties of model membranes and is limited to fluorescence, EPR, and MD studies. OxPCs are divided into two classes: A) hydroxy- or hydroperoxy-dieonyl phospatidylcholines, B) phospatidylcholines with oxidized and truncated chains with either aldehyde or carboxylic group. It was shown that the presence of the investigated oxPCs in phospholipid model membranes may have the following consequences: 1) decrease of the lipid order, 2) lowering of phase transition temperatures, 3) lateral expansion and thinning of the bilayer, 4) alterations of bilayer hydration profiles, 5) increased lipid mobility, 6) augmented flip-flop, 7) influence on the lateral phase organisation, and 8) promotion of water defects and, under extreme conditions (i.e. high concentrations of class B oxPCs), disintegration of the bilayer. The effects of class A oxPCs appear to be more moderate than those observed or predicted for class B. Many of the abovementioned findings are related to the ability of the oxidized chains of certain oxPCs to reorient toward the water phase. Some of the effects appear to be moderated by the presence of cholesterol. Although those biophysical alternations are found at oxPC concentrations higher than the total oxPC concentrations found under physiological conditions, certain organelles may reach such elevated oxPC concentrations locally. It is a challenge for the future to correlate the biophysics of oxidized phospholipids to metabolic studies in order to define the significance of the findings presented herein for pathophysiology. This article is part of a Special Issue entitled: Oxidized phospholipids-their properties and interactions with proteins.

Membrane Protein Simulations Using AMBER Force Field and Berger Lipid Parameters

AMBER force fields are among the most commonly used in molecular dynamics (MD) simulations of proteins. Unfortunately, they lack a specific set of lipid parameters, thus limiting its use in membrane protein simulations. In order to overcome this limitation we assessed whether the widely used united-atom lipid parameters described by Berger and co-workers could be used in conjunction with AMBER force fields in simulations of membrane proteins. Thus, free energies of solvation in water and in cyclohexane, and free energies of water to cyclohexane transfer, were computed by thermodynamic integration procedures for neutral amino acid side-chains employing AMBER99, AMBER03, and OPLS-AA amino acid force fields. In addition, MD simulations of three membrane proteins in a POPC lipid bilayer, the β2 adrenergic G protein-coupled receptor, Aquaporin-1, and the outer membrane protein Omp32, were performed with the aim of comparing the AMBER99SB/Berger combination of force fields with the OPLS-AA/Berger combination. We have shown that AMBER99SB and Berger force fields are compatible, they provide reliable free energy estimations relative to experimental values, and their combination properly describes both membrane and protein structural properties. We then suggest that the AMBER99SB/Berger combination is a reliable choice for the simulation of membrane proteins, which links the easiness of ligand parametrization and the ability to reproduce secondary structure of AMBER99SB force field with the largely validated Berger lipid parameters.

Molecular dynamics simulations of the Ca2+-pump: a structural analysis

We report large scale molecular dynamics computer simulations, ∼100 ns, of the ion pump Ca(2+)-ATPase immersed in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer. The structure simulated here, E1, one of the several conformations resolved using X-ray diffraction techniques, hosts two Ca(2+)-ions in the hydrophobic domain. Our results indicate that protonated residues lead to stronger ion-residue interactions, supporting previous conclusions regarding the sensitivity of the Ca(2+) behaviour to the protonated state of the amino acid binding sites. We also investigate how the protein perturbs the bilayer structure. We show that the POPC bilayer is ∼12% thinner than the pure bilayer, near the protein surface. This perturbation decays exponentially with the distance from the protein with a characteristic decay length of 0.8 nm. We find that the projected area per lipid also decreases near the protein. Using an analytical model we show that this change in the area is only apparent and it can be explained by considering the local curvature of the membrane. Our results indicate that the real area per lipid near the protein is not significantly modified with respect to the pure bilayer result. Further our results indicate that the local deformation of the membrane around the protein might be compatible with the enhanced protein activity observed in experiments over a narrow range of membrane thicknesses.

Purification and activity testing of the full-length YycFGHI proteins of Staphylococcus aureus

Background

The YycFG two-component regulatory system (TCS) of Staphylococcus aureus represents the only essential TCS that is almost ubiquitously distributed in Gram-positive bacteria with a low G+C-content. YycG (WalK/VicK) is a sensor histidine-kinase and YycF (WalR/VicR) is the cognate response regulator. Both proteins play an important role in the biosynthesis of the cell envelope and mutations in these proteins have been involved in development of vancomycin and daptomycin resistance.

Methodology/Principal Findings

Here we present high yield expression and purification of the full-length YycG and YycF proteins as well as of the auxiliary proteins YycH and YycI of Staphylococcus aureus. Activity tests of the YycG kinase and a mutated version, that harbours an Y306N exchange in its cytoplasmic PAS domain, in a detergent-micelle-model and a phosholipid-liposome-model showed kinase activity (autophosphorylation and phosphoryl group transfer to YycF) only in the presence of elevated concentrations of alkali salts. A direct comparison of the activity of the kinases in the liposome-model indicated a higher activity of the mutated YycG kinase. Further experiments indicated that YycG responds to fluidity changes in its microenvironment.

Conclusions/Significance

The combination of high yield expression, purification and activity testing of membrane and membrane-associated proteins provides an excellent experimental basis for further protein-protein interaction studies and for identification of all signals received by the YycFGHI system.

Aqueous solutions at the interface with phospholipid bilayers

In a sense, life is defined by membranes, because they delineate the barrier between the living cell and its surroundings. Membranes are also essential for regulating the machinery of life throughout many interfaces within the cell's interior. A large number of experimental, computational, and theoretical studies have demonstrated how the properties of water and ionic aqueous solutions change due to the vicinity of membranes and, in turn, how the properties of membranes depend on the presence of aqueous solutions. Consequently, understanding the character of aqueous solutions at their interface with biological membranes is critical to research progress on many fronts. The importance of incorporating a molecular-level description of water into the study of biomembrane surfaces was demonstrated by an examination of the interaction between phospholipid bilayers that can serve as model biological membranes. The results showed that, in addition to well-known forces, such as van der Waals and screened Coulomb, one has to consider a repulsion force due to the removal of water between surfaces. It was also known that physicochemical properties of biological membranes are strongly influenced by the specific character of the ions in the surrounding aqueous solutions because of the observation that different anions produce different effects on muscle twitch tension. In this Account, we describe the interaction of pure water, and also of aqueous ionic solutions, with model membranes. We show that a symbiosis of experimental and computational work over the past few years has resulted in substantial progress in the field. We now better understand the origin of the hydration force, the structural properties of water at the interface with phospholipid bilayers, and the influence of phospholipid headgroups on the dynamics of water. We also improved our knowledge of the ion-specific effect, which is observed at the interface of the phospholipid bilayer and aqueous solution, and its connection with the Hofmeister series. Nevertheless, despite substantial progress, many issues remain unresolved. Thus, for example, we still cannot satisfactorily explain the force of interaction between phospholipid bilayers immersed in aqueous solutions of NaI. Although we try to address many issues here, the scope of the discussion is limited and does not cover such important topics as the influence of ionic solutions on phases of bilayers, the influence of salts on the properties of Langmuir monolayers containing lipid molecules, or the influence of aqueous solutions on bilayers containing mixtures of lipids. We anticipate that the future application of more powerful experimental techniques, in combination with more advanced computational hardware, software, and theory, will produce molecular-level information about these important topics and, more broadly, will further illuminate our understanding of interfaces between aqueous solutions and biological membranes.

Force spectroscopy reveals the effect of different ions in the nanomechanical behavior of phospholipid model membranes: the case of potassium cation

How do metal cations affect the stability and structure of phospholipid bilayers? What role does ion binding play in the insertion of proteins and the overall mechanical stability of biological membranes? Investigators have used different theoretical and microscopic approaches to study the mechanical properties of lipid bilayers. Although they are crucial for such studies, molecular-dynamics simulations cannot yet span the complexity of biological membranes. In addition, there are still some experimental difficulties when it comes to testing the ion binding to lipid bilayers in an accurate way. Hence, there is a need to establish a new approach from the perspective of the nanometric scale, where most of the specific molecular phenomena take place. Atomic force microscopy has become an essential tool for examining the structure and behavior of lipid bilayers. In this work, we used force spectroscopy to quantitatively characterize nanomechanical resistance as a function of the electrolyte composition by means of a reliable molecular fingerprint that reveals itself as a repetitive jump in the approaching force curve. By systematically probing a set of bilayers of different composition immersed in electrolytes composed of a variety of monovalent and divalent metal cations, we were able to obtain a wealth of information showing that each ion makes an independent and important contribution to the gross mechanical resistance and its plastic properties. This work addresses the need to assess the effects of different ions on the structure of phospholipid membranes, and opens new avenues for characterizing the (nano)mechanical stability of membranes.

Concentration dependence of NaCl ion distributions around DPPC lipid bilayers

We study the coordination of excess NaCl to zwitterionic DPPC lipid bilayers using molecular dynamics simulations. We find that Na ions directly coordinate with the DPPC lipid carbonyl groups. As the number of excess ions increases, the number of coordinated ions increases, until it reaches a plateau at a ratio near 1 ion per every four lipids at 310 K, and 1 ion per every six lipids at 323 K. The area per lipid decreases as the number of excess ions is increased. For low number of ions per lipids (1:16 and 1:8), most Na ions are bound to the lipid carbonyls, while the Cl form an ionic cloud around the lipid choline groups. As a result of the Na binding, the lipid has an effective positive charge density. The residence time of Na ions bound to the lipid is longer than 40 ns, while Cl ions exchange faster than the nanoseconds timescale. We find that the bound Na ions replace ordered water around the carbonyls. The net linear charge density near the carbonyl groups stays positive, regardless of the presence of excess salt in the solution.

Validating affinities for ion-lipid association from simulation against experiment

Understanding biological membranes at physiological conditions requires comprehension of the interaction of lipid bilayers with sodium and potassium ions. These cations are adsorbed at palmitoyl-oleoyl-phosphatidylcholine (POPC) bilayers as indicated from previous studies. Here we compare the affinity of Na(+) and K(+) for POPC in molecular dynamics (MD) simulations with recent data from electrophoresis experiments and isothermal calorimetry (ITC) at neutral pH. NaCl and KCl were described using GROMOS or parameters matching solution activities on the basis of Kirkwood-Buff theory (KBFF), and K(+) was also described using parameters by Dang et al., all in conjunction with the Berger parameters for the lipids and the SPC water model. Apparent binding constants of GROMOS-Na(+) and KBFF-K(+) are the same within error and in good agreement with values from ITC. Although these force fields yield the same number of bound ions per number of lipids for Na(+) and K(+), they give a larger number of Na(+) ions per surface area compared to K(+), in agreement with the electrophoresis experiments, because Na(+) causes a stronger reduction in the area per lipid than K(+). The intrinsic binding constants, on the other hand, are reproduced by Dang-K(+) but overestimated by GROMOS-Na(+) and KBFF-K(+). That no ion force field reproduces the intrinsic and the apparent binding constant simultaneously arises from the fact that in MD simulations, implicitly meant to mimic neutral pH, pure PC is usually modeled with zero surface charge. In contrast, POPC at neutral conditions in experiment carries a low but significant negative surface charge and is uncharged only at acidic pH as indicated from electrophoretic mobilities. Implications for future simulation and experimental studies are discussed.

Interactions of molecular ions with model phospholipid membranes

The affinities of a series of biologically relevant ions for a hydrated phospholipid membrane were investigated using molecular dynamics simulation. Interactions of molecular ions, such as guanidinium, tetramethylammonium, and thiocyanate with the bilayer were computationally characterized for the first time. Simulations reveal strong ion specificity. On one hand, ions like guanidinium and thiocyanate adsorb relatively strongly to the headgroup region of the membrane. On the other hand, potassium or chloride interact very weakly with the phospholipids and merely act as neutralizing counterions. Calculations also show that these ions affect differently biophysical properties of the membrane, such as lipid diffusion, headgroup hydration and tilt angle.

Interaction of salicylate and a terpenoid plant extract with model membranes: reconciling experiments and simulations

We investigate the effects of two structurally similar small cyclic molecules: salicylic acid and perillic acid on a zwitterionic model lipid bilayer, and show that both molecules might have biological activity related to membrane thinning. Salicylic acid is a nonsteroidal antiinflammatory drug, some of the pharmacological properties of which arise from its interaction with the lipid bilayer component of the plasma membrane. Prior simulations show that salicylate orders zwitterionic lipid membranes. However, this is in conflict with Raman scattering and vesicle fluctuation analysis data, which suggest the opposite. We show using extensive molecular dynamics simulations, cumulatively >2.5 μs, that salicylic acid indeed disorders membranes with concomitant membrane thinning and that the conflict arose because prior simulations suffered from artifacts related to the sodium-ion induced condensation of zwitterionic lipids modeled by the Berger force field. Perillic acid is a terpenoid plant extract that has antiinfective and anticancer properties, and is extensively used in eastern medicine. We found that perillic acid causes large-scale membrane thinning and could therefore exert its antimicrobial properties via a membrane-lytic mechanism reminiscent of antimicrobial peptides. Being more amphipathic, perillic acid is more potent in disrupting lipid headgroup packing, and significantly modifies headgroup dipole orientation. Like salicylate, the membrane thinning effect of perillic acid is masked by the presence of sodium ions. As an alternative to sodium cations, we advocate the straightforward solution of using larger countercations like potassium or tetra-methyl-ammonium that will maintain electroneutrality but not interact strongly with, and thus not condense, the lipid bilayer.

Ca(2+) bridging of apposed phospholipid bilayers

In an effort to provide insight into the mechanism of Ca(2+)-induced fusion of lipid vesicles, molecular dynamics simulations in the isobaric-isothermal ensemble are used to investigate interactions of Ca(2+) with apposed lipid bilayers in close proximity. Simulations reveal the formation of a Ca(2+)-phospholipid "anhydrous complex" between apposed bilayers, whereas similar calculations performed with Na(+) display only complexation between neighboring lipids within the same bilayer. The binding of Ca(2+) to apposed phospholipids brings large regions of the bilayers into close contact (<4 Å), displacing water from phospholipid head groups in the process and creating regions of local dehydration. Dehydration of the apposed bilayers leads to ordering of the phospholipid tails, which is partially disrupted by the presence of Ca(2+)-phospholipid bridges.

Direct imaging of salt effects on lipid bilayer ordering at sub-molecular resolution

The interactions of salts with lipid bilayers are known to alter the properties of membranes and therefore influence their structure and dynamics. Sodium and calcium cations penetrate deeply into the headgroup region and bind to the lipids, whereas potassium ions only loosely associate with lipid molecules and mostly remain outside of the headgroup region. We investigated a dipalmitoylphosphatidylcholine (DPPC) bilayer in the gel phase in the presence of all three cations with a concentration of Ca²+ ions an order of magnitude smaller than the Na+ and K+ ions. Our findings indicate that the area per unit cell does not significantly change in these three salt solutions. However the lipid molecules do re-order non-isotropically under the influence of the three different cations. We attribute this reordering to a change in the highly directional intermolecular interactions caused by a variation in the dipole-dipole bonding arising from a tilt of the headgroup out of the membrane plane. Measurements in different NaCl concentrations also show a non-isotropic re-ordering of the lipid molecules.

Molecular dynamics simulations of mixed acidic/zwitterionic phospholipid bilayers

Anionic lipids are key components in the cell membranes. Many cell-regulatory and signaling mechanisms depend upon a complicated interplay between them and membrane-bound proteins. Phospholipid bilayers are commonly used as model systems in experimental or theoretical studies to gain insight into the structure and dynamics of biological membranes. We report here 200-ns-long MD simulations of pure (DMPC and DMPG) and mixed equimolar (DMPC/DMPG, DMPC/DMPS, and DMPC/DMPA) bilayers that each contain 256 lipids. The intra- and intermolecular interaction patterns in pure and mixed bilayers are analyzed and compared. The effect of monovalent ions (Na+) on the formation of salt-bridges is investigated. In particular, the number of Na(+)-mediated clusters in the presence of DMPS is higher than with DMPG and DMPA. We observe a preferential clustering of DMPS (and to some extent DMPA) lipids together rather than with DMPC molecules, which can explain the phase separation observed experimentally for DMPC/DMPS and DMPC/DMPA bilayers.

Vibrational sum frequency generation spectroscopic investigation of the interaction of thiocyanate ions with zwitterionic phospholipid monolayers at the air-water interface

Thiocyanate (SCN(-)) is a highly chaotropic anion of considerable biological significance, which interacts quite strongly with lipid interfaces. In most cases it is not exactly known if this interaction involves direct binding to lipid groups, or some type of indirect association or partitioning. Since thiocyanate is a linear ion, with a considerable dipole moment and nonspherical polarizability tensor, one should also consider its capability to adopt different or preferential orientations at lipid interfaces. In the present work, the interaction of thiocyanate anions with zwitterionic phospholipid monolayers in the liquid expanded (LE) phase is examined using surface pressure-area per molecule (pi-A(L)) isotherms and vibrational sum frequency generation (VSFG) spectroscopy. Both dipalmitoyl phosphatidylcholine (DPPC) and dimyristoyl phosphatidylethanolamine (DMPE) lipids, which form stable monolayers, have been used in this investigation, since their headgroups may be expected to interact with the electrolyte solution in different ways. The pi-A(L) isotherms of both lipids indicate a strong expansion of the monolayers when in contact with SCN(-) solutions. From the C-H stretch region of the VSFG spectra it can be deduced that the presence of the anion perturbs the conformation of the lipid chains significantly. The interfacial water structure is also perturbed in a complex way. Two distinct thiocyanate populations are detected in the CN stretch spectral region, proving that SCN(-) associates with zwitterionic phospholipids. Although this is a preliminary investigation of this complex system and more work is necessary to clarify certain points made in the discussion, a potential identification of the two SCN(-) populations and a molecular-level explanation for the observed effects of the SCN(-) on the VSFG spectra of the lipids is provided.

Mitochondrial membranes with mono- and divalent salt: changes induced by salt ions on structure and dynamics

We employ atomistic simulations to consider how mono- (NaCl) and divalent (CaCl(2)) salt affects properties of inner and outer membranes of mitochondria. We find that the influence of salt on structural properties is rather minute, only weakly affecting lipid packing, conformational ordering, and membrane electrostatic potential. The changes induced by salt are more prominent in dynamical properties related to ion binding and formation of ion-lipid complexes and lipid aggregates, as rotational diffusion of lipids is slowed down by ions, especially in the case of CaCl(2). In the same spirit, lateral diffusion of lipids is slowed down rather considerably for increasing concentration of CaCl(2). Both findings for dynamic properties can be traced to the binding of ions with lipid head groups and the related changes in interaction patterns in the headgroup region, where the binding of Na(+) and Ca(2+) ions is clearly different. The role of cardiolipins in these phenomena turns out to be important.

Ion dynamics in cationic lipid bilayer systems in saline solutions

Positively charged lipid bilayer systems are a promising class of nonviral vectors for safe and efficient gene and drug delivery. Detailed understanding of these systems is therefore not only of fundamental but also of practical biomedical interest. Here, we study bilayers comprising a binary mixture of cationic dimyristoyltrimethylammoniumpropane (DMTAP) and zwitterionic (neutral) dimyristoylphosphatidylcholine (DMPC) lipids. Using atomistic molecular dynamics simulations, we address the effects of bilayer composition (cationic to zwitterionic lipid fraction) and of NaCl electrolyte concentration on the dynamical properties of these cationic lipid bilayer systems. We find that, despite the fact that DMPCs form complexes via Na(+) ions that bind to the lipid carbonyl oxygens, NaCl concentration has a rather minute effect on lipid diffusion. We also find the dynamics of Cl(-) and Na(+) ions at the water-membrane interface to differ qualitatively. Cl(-) ions have well-defined characteristic residence times of nanosecond scale. In contrast, the binding of Na(+) ions to the carbonyl region appears to lack a characteristic time scale, as the residence time distributions displayed power-law features. As to lateral dynamics, the diffusion of Na(+) ions within the water-membrane interface consists of two qualitatively different modes of motion: very slow diffusion when ions are bound to DMPC, punctuated by fast rapid jumps when detached from the lipids. Overall, the prolonged dynamics of the Na(+) ions are concluded to be interesting for the physics of the whole membrane, especially considering its interaction dynamics with charged macromolecular surfaces.

Simulations of lipid transfer between a supported lipid bilayer and adsorbing vesicles

Recent experiments demonstrate transfer of lipid molecules between a charged, supported lipid membrane (SLB) and vesicles of opposite charge when the latter adsorb on the SLB. A simple phenomenological bead model has been developed to simulate this process. Beads were defined to be of three types, 'n', 'p', and '0', representing POPS (negatively charged), POEPC (positively charged), and POPC (neutral but zwitterionic) lipids, respectively. Phenomenological bead-bead interaction potentials and lipid transfer rate constants were used to account for the overall interaction and transfer kinetics. Using different bead mixtures in both the adsorbing vesicle and in the SLB (representing differently composed/charged vesicles and SLBs as in the reported experiments), we clarify under which circumstances a vesicle adsorbs to the SLB, and whether it, after lipid transfer and changed composition of the SLB and vesicle, desorbs back to the bulk again or not. With this model we can reproduce and provide a conceptual picture for the experimental findings.

Calculation of the electrostatic potential of lipid bilayers from molecular dynamics simulations: methodological issues

The electrostatic properties of lipid membranes are of profound importance as they are directly associated with membrane potential and, consequently, with numerous membrane-mediated biological phenomena. Here we address a number of methodological issues related to the computation of the electrostatic potential from atomic-scale molecular dynamics simulations of lipid bilayers. We discuss two slightly different forms of Poisson equation that are normally used to calculate the membrane potential: (i) a classical form when the potential and the electric field are chosen to be zero on one of the sides of a simulation box and (ii) an alternative form, when the potential is set to be the same on the opposite sides of a simulation box. Both forms differ by a position-dependent correction term, which has been shown to be proportional to the overall dipole moment of a bilayer system (for neutral systems). For symmetric bilayers we demonstrate that both approaches give essentially the same potential profiles, provided that simulations are long enough (a production run of at least 100 ns is required) and that fluctuations of the center of mass of a bilayer are properly accounted for. In contrast, for asymmetric lipid bilayers, the second approach is no longer appropriate due to a nonzero net dipole moment across a simulation box with a single asymmetric bilayer. We demonstrate that in this case the electrostatic potential can adequately be described by the classical form of Poisson equation, provided that it is employed in conjunction with tin-foil boundary conditions, which exactly balance a nonzero surface charge of a periodically replicated multibilayer system. Furthermore, we show that vacuum boundary conditions give qualitatively similar potential profiles for asymmetric lipid bilayers as compared to the conventional periodic boundaries, but accurate determination of the transmembrane potential difference is then hindered due to detachment of some water dipoles from bulk aqueous solution to vacuum.

Lipid gymnastics: evidence of complete acyl chain reversal in oxidized phospholipids from molecular simulations

In oxidative environments, biomembranes contain oxidized lipids with short, polar acyl chains. Two stable lipid oxidation products are PoxnoPC and PazePC. PoxnoPC has a carbonyl group, and PazePC has an anionic carboxyl group pendant at the end of the short, oxidized acyl chain. We have used MD simulations to explore the possibility of complete chain reversal in OXPLs in POPC-OXPL mixtures. The polar AZ chain of PazePC undergoes chain reversal without compromising the lipid bilayer integrity at concentrations up to 25% OXPL, and the carboxyl group points into the aqueous phase. Counterintuitively, the perturbation of overall membrane structural and dynamic properties is stronger for PoxnoPC than for PazePC. This is because of the overall condensing and ordering effect of sodium ions bound strongly to the lipids in the PazePC simulations. The reorientation of AZ chain is similar for two different lipid force fields. This work provides the first molecular evidence of the "extended lipid conformation" in phospholipid membranes. The chain reversal of PazePC lipids decorates the membrane interface with reactive, negatively charged functional groups. Such chain reversal is likely to exert a profound influence on the structure and dynamics of biological membranes, and on membrane-associated biological processes.

Molecular model of a cell plasma membrane with an asymmetric multicomponent composition: water permeation and ion effects

We present molecular dynamics simulations of a multicomponent, asymmetric bilayer in mixed aqueous solutions of sodium and potassium chloride. Because of the geometry of the system, there are two aqueous solution regions in our simulations: one mimics the intracellular region, and one mimics the extracellular region. Ion-specific effects are evident at the membrane/aqueous solution interface. Namely, at equal concentrations of sodium and potassium, sodium ions are more strongly adsorbed to carbonyl groups of the lipid headgroups. A significant concentration excess of potassium is needed for this ion to overwhelm the sodium abundance at the membrane. Ion-membrane interactions also lead to concentration-dependent and cation-specific behavior of the electrostatic potential in the intracellular region because of the negative charge on the inner leaflet. In addition, water permeation across the membrane was observed on a timescale of approximately 100 ns. This study represents a step toward the modeling of realistic biological membranes at physiological conditions in intracellular and extracellular environments.

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.

Statistical thermodynamics of biomembranes

An overview of the major issues involved in the statistical thermodynamic treatment of phospholipid membranes at the atomistic level is summarized: thermodynamic ensembles, initial configuration (or the physical system being modeled), force field representation as well as the representation of long-range interactions. This is followed by a description of the various ways that the simulated ensembles can be analyzed: area of the lipid, mass density profiles, radial distribution functions (RDFs), water orientation profile, deuterium order parameter, free energy profiles and void (pore) formation; with particular focus on the results obtained from our recent molecular dynamic (MD) simulations of phospholipids interacting with dimethylsulfoxide (Me(2)SO), a commonly used cryoprotective agent (CPA).

Study of the effect of Na+ and Ca2+ ion concentration on the structure of an asymmetric DPPC/DPPC + DPPS lipid bilayer by molecular dynamics simulation

A molecular dynamics simulation study of the steady and dynamic properties of an asymmetric phospholipid bilayer was carried out in the presence of sodium or calcium ions. The asymmetric lipid bilayer was seen to resemble a cellular membrane of an eukaryotic cell, which was modeled by dipalmitoylphosphatidylcholine (DPPC) and dipalmitoylphosphatidylserine (DPPS), placing the DPPS in one of the two leaflets of the lipid bilayer. From a numerical analysis of the simulated trajectories, information was obtained with atomic resolution for both membrane leaflet concerning the effect of bilayer asymmetry on different properties of the lipid/water interface, such as the translational diffusion coefficient and rotational relaxation time of the water molecules, lipid hydration, and residence time of water around different lipid atoms. In addition, information related to lipid conformation, and lipid-lipid interactions was also analyzed.

Free energy for the permeation of Na(+) and Cl(-) ions and their ion-pair through a zwitterionic dimyristoyl phosphatidylcholine lipid bilayer by umbrella integration with harmonic fourier beads

Understanding the mechanism of ion permeation across lipid bilayers is key to controlling osmotic pressure and developing new ways of delivering charged, drug-like molecules inside cells. Recent reports suggest ion-pairing as the mechanism to lower the free energy barrier for the ion permeation in disagreement with predictions from the simple electrostatic models. In this paper we quantify the effect of ion-pairing or charge quenching on the permeation of Na(+) and Cl(-) ions across DMPC lipid bilayer by computing the corresponding potentials of mean force (PMFs) using fully atomistic molecular dynamics simulations. We find that the free energy barrier to permeation reduces in the order Na(+)-Cl(-) ion-pair (27.6 kcal/mol) > Cl(-) (23.6 kcal/mol) > Na(+) (21.9 kcal/mol). Furthermore, with the help of these PMFs we derive the change in the binding free energy between the Na(+) and Cl(-) with respect to that in water as a function of the bilayer permeation depth. Despite the fact that the bilayer boosts the Na(+)-Cl(-) ion binding free energy by as high as 17.9 kcal/mol near its center, ion-pairing between such hydrophilic ions as Na(+) and Cl(-) does not assist their permeation. However, based on a simple thermodynamic cycle, we suggest that ion-pairing between ions of opposite charge and solvent philicity could enhance ion permeation. Comparison of the computed permeation barriers for Na(+) and Cl(-) ions with available experimental data supports this notion. This work establishes general computational methodology to address ion-pairing in fluid anisotropic media and details the ion permeation mechanism on atomic level.

Molecular dynamics simulations of cardiolipin bilayers

Cardiolipin is a key lipid component in the inner mitochondrial membrane, where the lipid is involved in energy production, cristae structure, and mechanisms in the apoptotic pathway. In this article we used molecular dynamics computer simulations to investigate cardiolipin and its effect on the structure of lipid bilayers. Three cardiolipin/POPC bilayers with different lipid compositions were simulated: 100, 9.2, and 0% cardiolipin. We found strong association of sodium counterions to the carbonyl groups of both lipid types, leaving in the case of 9.2% cardiolipin virtually no ions in the aqueous compartment. Although binding occurred primarily at the carbonyl position, there was a preference to bind to the carbonyl groups of cardiolipin. Ion binding and the small headgroup of cardiolipin gave a strong ordering of the hydrocarbon chains. We found significant effects in the water dipole orientation and water dipole potential which can compensate for the electrostatic repulsion that otherwise should force charged lipids apart. Several parameters relevant for the molecular structure of cardiolipin were calculated and compared with results from analyses of coarse-grained simulations and available X-ray structural data.

Determination of alkali and halide monovalent ion parameters for use in explicitly solvated biomolecular simulations

Alkali (Li(+), Na(+), K(+), Rb(+), and Cs(+)) and halide (F(-), Cl(-), Br(-), and I(-)) ions play an important role in many biological phenomena, roles that range from stabilization of biomolecular structure, to influence on biomolecular dynamics, to key physiological influence on homeostasis and signaling. To properly model ionic interaction and stability in atomistic simulations of biomolecular structure, dynamics, folding, catalysis, and function, an accurate model or representation of the monovalent ions is critically necessary. A good model needs to simultaneously reproduce many properties of ions, including their structure, dynamics, solvation, and moreover both the interactions of these ions with each other in the crystal and in solution and the interactions of ions with other molecules. At present, the best force fields for biomolecules employ a simple additive, nonpolarizable, and pairwise potential for atomic interaction. In this work, we describe our efforts to build better models of the monovalent ions within the pairwise Coulombic and 6-12 Lennard-Jones framework, where the models are tuned to balance crystal and solution properties in Ewald simulations with specific choices of well-known water models. Although it has been clearly demonstrated that truly accurate treatments of ions will require inclusion of nonadditivity and polarizability (particularly with the anions) and ultimately even a quantum mechanical treatment, our goal was to simply push the limits of the additive treatments to see if a balanced model could be created. The applied methodology is general and can be extended to other ions and to polarizable force-field models. Our starting point centered on observations from long simulations of biomolecules in salt solution with the AMBER force fields where salt crystals formed well below their solubility limit. The likely cause of the artifact in the AMBER parameters relates to the naive mixing of the Smith and Dang chloride parameters with AMBER-adapted Aqvist cation parameters. To provide a more appropriate balance, we reoptimized the parameters of the Lennard-Jones potential for the ions and specific choices of water models. To validate and optimize the parameters, we calculated hydration free energies of the solvated ions and also lattice energies (LE) and lattice constants (LC) of alkali halide salt crystals. This is the first effort that systematically scans across the Lennard-Jones space (well depth and radius) while balancing ion properties like LE and LC across all pair combinations of the alkali ions and halide ions. The optimization across the entire monovalent series avoids systematic deviations. The ion parameters developed, optimized, and characterized were targeted for use with some of the most commonly used rigid and nonpolarizable water models, specifically TIP3P, TIP4P EW, and SPC/E. In addition to well reproducing the solution and crystal properties, the new ion parameters well reproduce binding energies of the ions to water and the radii of the first hydration shells.