Atomic-scale structure and electrostatics of anionic palmitoyloleoylphosphatidylglycerol lipid bilayers with Na+ counterions

Phase transitions and spatially ordered counterion association in ionic-lipid membranes: theory versus experiment

Aqueous dispersions of phosphatidylglycerol (PG) lipids may present an anomalous chain-melting transition at low ionic strengths, as seen by different experimental techniques such as calorimetry or light scattering. The anomaly disappears at high ionic strengths or for longer acyl-chain lengths. In this article, we use a statistical model for the bilayer that distinguishes both lipid chain and headgroup states in order to compare model and experimental thermotropic and electrical properties. The effective van der Waals interactions among hydrophobic chains compete with the electrostatic repulsions between polar headgroups, which may be ionized (counterion dissociated) or electrically neutral (associated with counterions). Electric degrees of freedom introduce new thermotropic charge-ordered phases in which headgroup charges may be spatially ordered, depending on the electrolyte ionic strength, introducing a new rationale for experimental data on PGs. The thermal phases presented by the model for different chain lengths, at fixed ionic strength, compare well with an experimental phase diagram constructed on the basis of differential scanning calorimetry profiles. In the case of dispersions of DMPG (dimyristoyl phosphatidylglycerol) with added monovalent salt, the model properties reproduce the main features displayed by data from differential scanning calorimetry as well as the characteristic profile for the degree of ionization of the bilayer surface across the anomalous transition region, obtained from the theoretical interpretation of electrokinetic (conductivity and electrophoretic mobility) measurements.

The molecular basis for antimicrobial activity of pore-forming cyclic peptides

The mechanism of action of antimicrobial peptides is, to our knowledge, still poorly understood. To probe the biophysical characteristics that confer activity, we present here a molecular-dynamics and biophysical study of a cyclic antimicrobial peptide and its inactive linear analog. In the simulations, the cyclic peptide caused large perturbations in the bilayer and cooperatively opened a disordered toroidal pore, 1-2 nm in diameter. Electrophysiology measurements confirm discrete poration events of comparable size. We also show that lysine residues aligning parallel to each other in the cyclic but not linear peptide are crucial for function. By employing dual-color fluorescence burst analysis, we show that both peptides are able to fuse/aggregate liposomes but only the cyclic peptide is able to porate them. The results provide detailed insight on the molecular basis of activity of cyclic antimicrobial peptides.

Phase transitions and spatially ordered counterion association in ionic-lipid membranes: a statistical model

We propose a statistical model to account for the gel-fluid anomalous phase transitions in charged bilayer- or lamellae-forming ionic lipids. The model Hamiltonian comprises effective attractive interactions to describe neutral-lipid membranes as well as the effect of electrostatic repulsions of the discrete ionic charges on the lipid headgroups. The latter can be counterion dissociated (charged) or counterion associated (neutral), while the lipid acyl chains may be in gel (low-temperature or high-lateral-pressure) or fluid (high-temperature or low-lateral-pressure) states. The system is modeled as a lattice gas with two distinct particle types--each one associated, respectively, with the polar-headgroup and the acyl-chain states--which can be mapped onto an Ashkin-Teller model with the inclusion of cubic terms. The model displays a rich thermodynamic behavior in terms of the chemical potential of counterions (related to added salt concentration) and lateral pressure. In particular, we show the existence of semidissociated thermodynamic phases related to the onset of charge order in the system. This type of order stems from spatially ordered counterion association to the lipid headgroups, in which charged and neutral lipids alternate in a checkerboard-like order. Within the mean-field approximation, we predict that the acyl-chain order-disorder transition is discontinuous, with the first-order line ending at a critical point, as in the neutral case. Moreover, the charge order gives rise to continuous transitions, with the associated second-order lines joining the aforementioned first-order line at critical end points. We explore the thermodynamic behavior of some physical quantities, like the specific heat at constant lateral pressure and the degree of ionization, associated with the fraction of charged lipid headgroups.

Competing interactions for antimicrobial selectivity based on charge complementarity

Antimicrobial peptides (AMPs) are an evolutionary conserved component of the innate immune system and possible templates for the development of new antibiotics. An important property of antimicrobial peptides is their ability to discriminate bacterial from eucaryotic cells which is attributed to the difference in lipid composition of the outer leaflet of the plasma membrane between the two types of cells. Whereas eucaryotic cells usually expose zwitterionic lipids, procaryotic cells expose also anionic lipids which bind the cationic antimicrobial peptides electrostatically. An example is the antimicrobial peptide NK-2 which is highly cationic and favors binding to anionic membranes. In the present study, the difference in binding affinity of NK-2 for palmitoyl-oleoyl-phosphatidyl-glycerol (POPG) and palmitoyl-oleoyl-phosphatidyl-choline (POPC) is studied using molecular dynamics simulations in conjunction with a coarse grained model and thermodynamic integration, by computing the change in free energy and its components upon the transfer of NK-2 from POPC to POPG. The transfer is indeed found to be highly favorable. Interestingly, the favorable contribution from the electrostatic interaction between the peptide and the anionic lipids is overcompensated by an unfavorable contribution from the change in lipid-cation interactions due to the release of counterions from the lipids. The increase in entropy due to the release of the cations is compensated by other entropic components. The largest favorable contribution arises from the solvation of the counterions. Overall the interaction between NK-2 and POPG is not determined by a single driving force but a subtle balance of competing interactions.

Study of PEGylated lipid layers as a model for PEGylated liposome surfaces: molecular dynamics simulation and Langmuir monolayer studies

We have combined Langmuir monolayer film experiments and all-atom molecular dynamics (MD) simulation of a bilayer to study the surface structure of a PEGylated liposome and its interaction with the ionic environment present under physiological conditions. Lipids that form both gel and liquid-crystalline membranes have been used in our study. By varying the salt concentration in the Langmuir film experiment and including salt at the physiological level in the simulation, we have studied the effect of salt ions present in the blood plasma on the structure of the poly(ethylene glycol) (PEG) layer. We have also studied the interaction between the PEG layer and the lipid bilayer in both the liquid-crystalline and gel states. The MD simulation shows two clear results: (a) The Na(+) ions form close interactions with the PEG oxygens, with the PEG chains forming loops around them and (b) PEG penetrates the lipid core of the membrane for the case of a liquid-crystalline membrane but is excluded from the tighter structure of the gel membrane. The Langmuir monolayer results indicate that the salt concentration affects the PEGylated lipid system, and these results can be interpreted in a fashion that is in agreement with the results of our MD simulation. We conclude that the currently accepted picture of the PEG surface layer acting as a generic neutral hydrophilic polymer entirely outside the membrane, with its effect explained through steric interactions, is not sufficient. The phenomena we have observed may affect both the interaction between the liposome and bloodstream proteins and the liquid-crystalline-gel transition and is thus relevant to nanotechnological drug delivery device design.

Structural basis for the enhanced activity of cyclic antimicrobial peptides: the case of BPC194

We report the molecular basis for the differences in activity of cyclic and linear antimicrobial peptides. We iteratively performed atomistic molecular dynamics simulations and biophysical measurements to probe the interaction of a cyclic antimicrobial peptide and its inactive linear analogue with model membranes. We establish that, relative to the linear peptide, the cyclic one binds stronger to negatively charged membranes. We show that only the cyclic peptide folds at the membrane interface and adopts a β-sheet structure characterised by two turns. Subsequently, the cyclic peptide penetrates deeper into the bilayer while the linear peptide remains essentially at the surface. Finally, based on our comparative study, we propose a model characterising the mode of action of cyclic antimicrobial peptides. The results provide a chemical rationale for enhanced activity in certain cyclic antimicrobial peptides and can be used as a guideline for design of novel antimicrobial peptides.

Molecular dynamics simulations of the interactions of kinin peptides with an anionic POPG bilayer

We have performed molecular dynamics simulations of peptide hormone bradykinin (BK) and its fragment des-Arg9-BK in the presence of an anionic lipid bilayer, with an aim toward delineating the mechanism of action related to their bioactivity. Starting from the initial aqueous environment, both of the peptides are quickly adsorbed and stabilized on the cell surface. Whereas BK exhibits a stronger interaction with the membrane and prefers to stay on the interface, des-Arg9-BK, with the loss of C-terminal Arg, penetrates further. The heterogeneous lipid-water interface induces β-turn-like structure in the otherwise inherently flexible peptides. In the membrane-bound state, we observed C-terminal β-turn formation in BK, whereas for des-Arg9-BK, with the deletion of Arg9, turn formation occurred in the middle of the peptide. The basic Arg residues anchor the peptide to the bilayer by strong electrostatic interactions with charged lipid headgroups. Simulations with different starting orientations of the peptides with respect to the bilayer surface lead to the same observations, namely, the relative positioning of the peptides on the membrane surface, deeper penetration of the des-Arg9-BK, and the formation of turn structures. The lipid headgroups adjacent to the bound peptides become substantially tilted, causing bilayer thinning near the peptide contact region and increase the degree of disorder in nearby lipids. Again, because of hydrogen bonding with the peptide, the neighboring lipid's polar heads exhibit considerably reduced flexibility. Corroborating findings from earlier experiments, our results provide important information about how the lipid environment promotes peptide orientation/conformation and how the peptide adapts to the environment.

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.

Structure-activity relationships of an antimicrobial peptide plantaricin s from two-peptide class IIb bacteriocins

Class IIb bacteriocins are ribosomally synthesized antimicrobial peptides comprising two different peptides synergistically acting in equal amounts for optimal potency. In this study, we demonstrate for the first time potent (nanomolar) antimicrobial activity of a representative class IIb bacteriocin, plantaricin S (Pls), against four pathogenic gram-positive bacteria, including Listeria monocytogenes. The structure-activity relationships for Pls were studied using activity assays, circular dichroism (CD), and molecular dynamics (MD) simulations. The two Pls peptides and five Pls derived fragments were synthesized. The CD spectra of the Pls and selected fragments revealed helical conformations in aqueous 2,2,2-trifluoroethanol. The MD simulations showed that when the two Pls peptides are in antiparallel orientation, the helical regions interact and align, mediated by strong attraction between conserved GxxxG/AxxxA motifs. The results strongly correlate with the antimicrobial activity suggesting that helix-helix alignment of the two Pls peptides and interaction between the conserved motifs are crucial for interaction with the target cell membrane.

Headgroup mediated water insertion into the DPPC bilayer: a molecular dynamics study

Molecular dynamics simulation was performed on the 1,2-dipalmitoyl-sn-phosphocholine (DPPC) bilayer-water system using the GROMOS96 53a6 united atom force field. The transferability of force field was tested by reproducing the area per lipid within 3% accuracy from the experimental value. The simulation shows that water can penetrate much deeper inside the bilayer almost up to the starting point of the aliphatic chain. There is significant evidence from experiments that water goes deep in the DPPC bilayer, but it has not been reported from theoretical work. The mechanism of insertion of water deep inside the lipid bilayer is still not clear. In this report, for the first time, the mechanism of water insertion deep into the bilayer has been proposed. Water transport occurs by the headgroup and its first solvation shell. The trimethyl ammonium (NMe(3)) group (headgroup of DPPC) has two stable conformations at the bilayer-water interface, one outside the bilayer and another inside it. The NMe(3) group has a large clustering of water around it and takes the water molecules inside the bilayer with it during its entry into the bilayer. The water molecules penetrate into the bilayer with the help of the NMe(3) group present at the headgroup of DPPC and eventually form hydrogen bonds with carbonyl oxygen present deep inside the bilayer. Structural characteristics at the bilayer-water interface region are also reported.

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.

Potentials of mean force and permeabilities for carbon dioxide, ammonia, and water flux across a Rhesus protein channel and lipid membranes

As a member of the ubiquitous ammonium transporter/methylamine permease/Rhesus (Amt/MEP/Rh) family of membrane protein channels, the 50 kDa Rhesus channel (Rh50) has been implicated in ammonia (NH(3)) and, more recently, also in carbon dioxide (CO(2)) transport. Here we present molecular dynamics simulations of spontaneous full permeation events of ammonia and carbon dioxide across Rh50 from Nitrosomonas europaea. The simulations show that Rh50 is functional in its crystallographic conformation, without the requirement for a major conformational change or the action of a protein partner. To assess the physiological relevance of NH(3) and CO(2) permeation across Rh50, we have computed potentials of mean force (PMFs) and permeabilities for NH(3) and CO(2) flux across Rh50 and compare them to permeation through a wide range of lipid membranes, either composed of pure lipids or composed of lipids plus an increasing cholesterol content. According to the PMFs, Rh50 is expected to enhance NH(3) flux across dense membranes, such as membranes with a substantial cholesterol content. Although cholesterol reduces the intrinsic CO(2) permeability of lipid membranes, the CO(2) permeabilities of all membranes studied here are too high to allow significant Rh50-mediated CO(2) flux. The increased barrier in the PMF for water permeation across Rh50 shows that Rh50 discriminates 40-fold between water and NH(3). Thus, Rh50 channels complement aquaporins, allowing the cell to regulate water and NH(3) flux independently. The PMFs for methylamine and NH(3) are virtually identical, suggesting that methylamine provides an excellent model for NH(3) in functional experiments.

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.

Molecular order in a chromonic liquid crystal: a molecular simulation study of the anionic azo dye sunset yellow

We have carried out a detailed atomic simulation study of molecular order within a chromonic liquid crystalline material (sunset yellow) in aqueous solution. Self-assembly occurs in dilute solutions to form stacked aggregates, which show a preference for head-to-tail stacking and antiparallel dipole order. This feature is independent of solution concentration and aggregate size. Stacks are found to be dynamic entities in which rotational transitions (flips) can occur between antiparallel and parallel configurations. At a concentration matching the nematic phase of sunset yellow, the simulations show chromonic columns with a loose hexagonal packing and an intercolumn distance of 2.36 nm. Partial condensation of sodium ions occurs around a chromonic stack, with two preferred binding sites identified for sodium ions, corresponding to strong binding with the oxygens of a sulfonate group and a bridging site between a pair of molecules in a stack. A value for the free energy of binding of a molecule to a stack of 7 k(B)T was obtained for stacks of three and eight molecules, with a slightly larger value (additional 2 kJ mol(-1)) obtained for the dimer binding energy, indicating that aggregation is approximately isodesmic.

Electrically addressable, biologically relevant surface-supported bilayers

The assembly of electrically addressable, planar supported bilayers composed of biologically relevant lipids, such as those used in vesicular systems, will greatly enhance the experimental capabilities in membrane and membrane protein research. Here we assess the electrical properties of bilayers composed of a wide range of physiologically relevant lipids and lipid combinations. We demonstrate that robust, biologically relevant, planar supported bilayers with high resistance composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 25 mol % cholesterol can be constructed with high reproducibility. Furthermore, to enable studies of pore-forming peptides, which are commonly cationic, we demonstrate the construction of bilayers with biologically relevant outer leaflets incorporating up to 10 mol % negatively charged lipids. Unique features of the platform are that (1) the substrate is commercially available, atomically smooth, single-crystal silicon, (2) the polymer cushion allows for the natural incorporation of membrane proteins, and (3) the platform is highly reproducible.

Protein folding in membranes

Separation of cells and organelles by bilayer membranes is a fundamental principle of life. Cellular membranes contain a baffling variety of proteins, which fulfil vital functions as receptors and signal transducers, channels and transporters, motors and anchors. The vast majority of membrane-bound proteins contain bundles of alpha-helical transmembrane domains. Understanding how these proteins adopt their native, biologically active structures in the complex milieu of a membrane is therefore a major challenge in today's life sciences. Here, we review recent progress in the folding, unfolding and refolding of alpha-helical membrane proteins and compare the molecular interactions that stabilise proteins in lipid bilayers. We also provide a critical discussion of a detergent denaturation assay that is increasingly used to determine membrane-protein stability but is not devoid of conceptual difficulties.

Free energy profile of the interaction between a monomer or a dimer of protegrin-1 in a specific binding orientation and a model lipid bilayer

The free energies of adsorption of the monomer or dimer of the cationic beta-hairpin antimicrobial peptide protegrin-1 (PG1) in a specific binding orientation on a lipid bilayer are determined using molecular dynamics (MD) simulations and Poisson-Boltzmann calculations. The bilayer is composed of anionic palmitoyl-oleoyl-phosphatidylglycerol (POPG) and palmitoyl-oleoyl-phosphatidylethanolamine (POPE) with ratio 1:3 (POPG/POPE). PG1 is believed to kill bacteria by binding on their membranes. There, it forms pores that lyse the bacteria. Herein we focus on the thermodynamics of binding. In particular, we explore the role of counterion release from the lipid bilayer upon adsorption of either the monomeric or the dimeric form of PG1. Twenty-two 4-ns-long MD trajectories of equilibrated systems are generated to determine the free energy profiles for the monomer and dimer as a function of the distance between the peptide(s) and the membrane surface. The MD simulations are conducted at 11 different separations from the membrane for each of the two systems, one with PG1, the second with a PG1 dimer of only a specific orientation of the monomer and dimer without taking into account the change of entropy for the peptide. To calculate the potential of mean force for each peptide/membrane system, a variant of constrained MD and thermodynamic integration is used. We observed that PG1 dimer binds more favorably to the POPG/POPE membrane. A simple method for relating the free energy profile to the PG1-membrane binding constant is employed to predict a free energy of adsorption of -2.4 +/- 0.8 kcal/mol. A corresponding PG1-dimer-membrane binding constant is calculated as -3.5 +/- 1.1 kcal/mol. Free energy profiles from MD simulation were extensively analyzed and compared with results of Poisson-Boltzmann theory. We find the peptide-membrane attraction to be dominated by the entropy increase due to the release of counterions in a POPG/POPE lipid bilayer.

Structural basis of transmembrane domain interactions in integrin signaling

Cell surface receptors of the integrin family are pivotal to cell adhesion and migration. The activation state of heterodimeric alphabeta integrins is correlated to the association state of the single-pass alpha and beta transmembrane domains. The association of integrin alphaIIbbeta3 transmembrane domains, resulting in an inactive receptor, is characterized by the asymmetric arrangement of a straight (alphaIIb) and tilted (beta3) helix relative to the membrane in congruence to the dissociated structures. This allows for a continuous association interface centered on helix-helix glycine-packing and an unusual alphaIIb(GFF) structural motif that packs the conserved Phe-Phe residues against the beta3 transmembrane helix, enabling alphaIIb(D723)beta3(R995) electrostatic interactions. The transmembrane complex is further stabilized by the inactive ectodomain, thereby coupling its association state to the ectodomain conformation. In combination with recently determined structures of an inactive integrin ectodomain and an activating talin/beta complex that overlap with the alphabeta transmembrane complex, a comprehensive picture of integrin bi-directional transmembrane signaling has emerged.

Properties of hydrated excess protons near phospholipid bilayers

The behavior of the hydrated excess proton near different lipid membranes is studied with the third generation of the multistate empirical valence bond (MS-EVB3) model [Wu, Y. J.; Chen, H. N.; Wang, F.; Paesani, F.; Voth, G. A. J. Phys. Chem. B 2008, 112, 467]. Dioleoylphosphatidylcholine (DOPC), dioleoylphosphatidylethanolamine (DOPE), and dioleoylphosphatidylglycerol (DOPG) are selected as example lipids. In spite of the differences of the head groups, the molecular dynamics simulations show that all the lipid membranes have a proton-collecting antenna effect with no free energy barrier between the bulk water and interface regions. By comparison with classical hydronium model simulations, it is found that an appropriate description of proton Grotthuss shuttling and associated charge defect delocalization are necessary to obtain the correct free energy profile for the hydrated excess proton. In addition, nanosecond time scale sampling is essential to evaluate the free energy profiles, because certain slow motions are needed to stabilize the excess proton in the deep membrane interface region. It is also found that the lateral diffusion coefficients are 1 order of magnitude smaller in the interface region than in bulk water for all the lipids. These coefficients are almost the same as those of the lipid head groups. Finally, since the lipid phosphates may possibly be protonated due to the proton antenna effect of the membrane, phosphate group protonation is investigated and discussed within the MS-EVB framework.

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.

Cause and effect of melittin-induced pore formation: a computational approach

Melittin embedded in a palmitoyl oleyl phosphatidylcholine bilayer at a high peptide/lipid ratio (1:30) was simulated in the presence of explicit water and ions. The simulation results indicate the incipience of an ion-permeable water pore through collective membrane perturbation by bound peptides. The positively charged residues of melittin not only act as "anchors" but also disrupt the membrane, leading to cell lysis. A detailed analysis of the lipid tail order parameter profile depicts localized membrane perturbation. The lipids in the vicinity of the aqueous cavity adopt a tilted conformation, which allows local bilayer thinning. The prepore thus formed can be considered as the melittin-induced structural defects in the bilayer membrane. Because of the strong cationic nature, the melittin-induced prepore exhibits selectivity toward anions over cations. As Cl(-) ions entered into the prepore, they are electrostatically entrapped by positively charged residues located at its wall. The confined motion of the Cl(-) ions in the membrane interior is obvious from calculated diffusion coefficients. Moreover, reorientation of the local lipids occurs in such a way that few lipid heads along with peptide helices can line the surface of the penetrating aqueous phase. The flipping of lipids argued in favor of melittin-induced toroidal pore over a barrel-stave mechanism. Thus, our result provides atomistic level details of the mechanism of membrane disruption by antimicrobial peptide melittin.

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.

Water replacement hypothesis in atomic detail--factors determining the structure of dehydrated bilayer stacks

According to the water replacement hypothesis, trehalose stabilizes dry membranes by preventing the decrease of spacing between membrane lipids under dehydration. In this study, we use molecular-dynamics simulations to investigate the influence of trehalose on the area per lipid (APL) and related structural properties of dehydrated bilayers in atomic detail. The starting conformation of a palmitoyloleolylphosphatidylcholine lipid bilayer in excess water was been obtained by self-assembly. A series of molecular-dynamics simulations of palmitoyloleolylphosphatidylcholine with different degrees of dehydration (28.5, 11.7, and 5.4 waters per lipid) and different molar trehalose/lipid ratios (<1:1, 1:1, and >1:1) were carried out in the NPT ensemble. Water removal causes the formation of multilamellar "stacks" through periodic boundary conditions. The headgroups reorient from pointing outward to inward with dehydration. This causes changes in the electrostatic interactions between interfaces, resulting in interface interpenetration. Interpenetration creates self-spacing of the bilayers and prevents gel-phase formation. At lower concentrations, trehalose does not separate the interfaces, and acting together with self-spacing, it causes a considerable increase of APL. APL decreases at higher trehalose concentrations when the layer of sugar physically separates the interfaces. When interfaces are separated, the model confirms the water replacement hypothesis.

Binding free energy and counterion release for adsorption of the antimicrobial peptide lactoferricin B on a POPG membrane

Molecular dynamics (MD) simulations are used to study the interaction of an anionic palmitoyl-oleoyl-phosphatidylglycerol (POPG) bilayer with the cationic antimicrobial peptide bovine lactoferricin (LFCinB) in a 100 mM NaCl solution at 310 K. The interaction of LFCinB with a POPG bilayer is employed as a model system for studying the details of membrane adsorption selectivity of cationic antimicrobial peptides. Seventy eight 4 ns MD production run trajectories of the equilibrated system, with six restrained orientations of LFCinB at 13 different separations from the POPG membrane, are generated to determine the free energy profile for the peptide as a function of the distance between LFCinB and the membrane surface. To calculate the profile for this relatively large system, a variant of constrained MD and thermodynamic integration is used. A simplified method for relating the free energy profile to the LFCinB-POPG membrane binding constant is employed to predict a free energy of adsorption of -5.4+/-1.3 kcal/mol and a corresponding maximum adsorption binding force of about 58 pN. We analyze the results using Poisson-Boltzmann theory. We find the peptide-membrane attraction to be dominated by the entropy increase due to the release of counterions and polarized water from the region between the charged membrane and peptide, as the two approach each other. We contrast these results with those found earlier for adsorption of LFCinB on the mammalianlike palmitoyl-oleoyl-phosphatidylcholine membrane.

Models for phosphatidylglycerol lipids put to a structural test

Three atomistic empirical models for phosphatidylglycerol (PG) lipids are tested against structural data in the crystal and liquid crystal states. Simulations of the anhydrous crystal of dimyristoyl-phosphatidylglycerol (DMPG) show that only the CHARMM force field describes the conformation and interactions of PG head groups accurately. The other two models do not reproduce the native network of hydrogen bonds, suggesting the presence of biases in their conformational and nonbonded interaction properties. The CHARMM model is further validated in the biologically relevant liquid crystal phase by comparing experimental small-angle X-ray scattering spectra from DMPG unilamellar vesicles with data calculated from fluid bilayer simulations. The good agreement found in this model-free comparison implies that liquid crystal PG bilayers as described by CHARMM exhibit realistic bilayer thickness and lateral packing. Last, this model is used to simulate a fluid bilayer of palmitoyl-oleoyl-phosphatidylglycerol (POPG). The resulting view of the POPG bilayer structure is at variance with that proposed previously based on simulations, in particular, with respect to lateral packing of head groups and the role of counterions.

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.

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.

Membrane perturbation by the antimicrobial peptide PMAP-23: a fluorescence and molecular dynamics study

Several bioactive peptides exert their biological function by interacting with cellular membranes. Structural data on their location inside lipid bilayers are thus essential for a detailed understanding of their mechanism of action. We propose here a combined approach in which fluorescence spectroscopy and molecular dynamics (MD) simulations were applied to investigate the mechanism of membrane perturbation by the antimicrobial peptide PMAP-23. Fluorescence spectra, depth-dependent quenching experiments, and peptide-translocation assays were employed to determine the location of the peptide inside the membrane. MD simulations were performed starting from a random mixture of water, lipids and peptide, and following the spontaneous self-assembly of the bilayer. Both experimental and theoretical data indicated a peptide location just below the polar headgroups of the membrane, with an orientation essentially parallel to the bilayer plane. These findings, together with experimental results on peptide-induced leakage from large and giant vesicles, lipid flip-flop and peptide exchange between vesicles, support a mechanism of action consistent with the "carpet" model. Furthermore, the atomic detail provided by the simulations suggested the occurrence of an additional, more specific and novel mechanism of bilayer destabilization by PMAP-23, involving the unusual insertion of charged side chains into the hydrophobic core of the membrane.

Effect of lipid composition on buforin II structure and membrane entry

Buforin II is a 21-amino acid polycationic antimicrobial peptide derived from a peptide originally isolated from the stomach tissue of the Asian toad Bufo bufo gargarizans. It is hypothesized to target a wide range of bacteria by translocating into cells without membrane permeabilization and binding to nucleic acids. Previous research found that the structure and membrane interactions of buforin II are related to lipid composition. In this study, we used molecular dynamics (MD) simulations along with lipid vesicle experiments to gain insight into how buforin II interacts differently with phosphatidylcholine (PC), phosphatidylglycerol (PG), and phosphatidylethanolamine (PE) lipids. Fluorescent spectroscopic measurements agreed with the previous assertion that buforin II does not interact with pure PC vesicles. Nonetheless, the reduced entry of the peptide into anionic PG membranes versus neutral PC membranes during simulations correlates with the experimentally observed reduction in BF2 translocation through pure PG membranes. Simulations showing membrane entry into PC also provide insight into how buforin II may initially penetrate cell membranes. Our MD simulations also allowed us to consider how neutral PE lipids affect the peptide differently than PC. In particular, the peptide had a more helical secondary structure in simulations with PE lipids. A change in structure was also apparent in circular dichroism measurements. PE also reduced membrane entry in simulations, which correlates with decreased translocation in the presence of PE observed in previous studies. Together, these results provide molecular-level insight into how lipid composition can affect buforin II structure and function and will be useful in efforts to design peptides with desired antimicrobial and cell-penetrating properties.

Lipid Models for United-Atom Molecular Dynamics Simulations of Proteins

United-atom force fields for molecular dynamics (MD) simulations provide a higher computational efficiency, especially in lipid membrane simulations, with little sacrifice in accuracy, when compared to all-atom force fields. Excellent united-atom lipid models are available, but in combination with depreciated protein force fields. In this work, a united-atom model of the lipid 1,2-dipalmitoyl-sn-glycero-3-phosphocholine has been built with standard parameters of the force field GROMOS96 53a6 that reproduces the experimental area per lipid of a lipid bilayer within 3% accuracy to a value of 0.623 ± 0.011 nm(2) without the assumption of a constant surface area or the inclusion of surface pressure. In addition, the lateral self-diffusion constant and deuterium order parameters of the acyl chains are in agreement with experimental data. Furthermore, models for 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol (POPG) result in areas per lipid of 0.625 nm(2) (DMPC), 0.693 nm(2) (POPC), and 0.700 nm(2) (POPG) from 40 ns MD simulations. Experimental lateral self-diffusion coefficients are reproduced satisfactorily by the simulation. The lipid models can form the basis for molecular dynamics simulations of membrane proteins with current and future versions of united-atom protein force fields.

Calcium binding and head group dipole angle in phosphatidylserine-phosphatidylcholine bilayers

Manipulating the plasma membrane, the gateway to the cell interior, with chemical and physical agents for genetic and pharmacological therapy, and understanding the interactions of lipid membrane components with proteins and other structural and functional elements of the cell, require a detailed biomolecular membrane model. We report here progress along one path toward such a model: molecular dynamics simulations of mixed, zwitterionic-anionic, asymmetric phospholipid bilayers with monovalent and divalent cations. With phosphatidylcholine/phosphatidylserine systems, we identify temporal and concentration boundaries for equilibration of calcium with the bilayer and saturation of the calcium capacity of the membrane, we demonstrate the electrostatic- and entropic-driven associations of calcium and sodium ions with polar groups in the bilayer interface region, expressed in spatial distribution profiles and in changes in the orientation of the phospholipid head groups, and we describe for the first time simulations of dynamic, calcium-mediated adjustments in the conformation of mixed phospholipid species coresident in the same leaflet of the bilayer. The results are consistent with experimental observations and point the way to further refinement and increased realism of these molecular models.