Binding and insertion of alpha-helical anti-microbial peptides in POPC bilayers studied by molecular dynamics simulations

We have performed molecular dynamics simulations of the interactions of two alpha-helical anti-microbial peptides, magainin2 and its synthetic analog of MSI-78, with palmitoyl-oleoyl-phosphatidylcholine (POPC) lipid bilayers. We used various initial positions and orientations of the peptide with respect to the lipid bilayer, including a surface-bound state parallel to the interface, a trans-membrane state, and a partially inserted state. Our 20 ns long simulations show that both magainin2 and MSI-78 are most stable in the lipid environment, with the peptide destabilized to different extents in both aqueous and lipid/water interfacial environments. We found that there are strong specific interactions between the lysine residues of the peptides and the lipid head-group regions. MSI-78, owing to its large number of lysines, shows better binding characteristics and overall stability when compared to magainin2. We also find that both peptides destabilize the bilayer environment, as observed by the increase in lipid tail disorder and the induction of local curvature on the lipid head-groups by the peptides. From all the simulations, we conclude that the hydrogen bonding interactions between the lysines of the peptides and the oxygens of the polar lipid head-groups are the strongest and determine the overall peptide binding characteristics to the lipids.     

Effective charge of melittin upon interaction with POPC vesicles

The binding of bee venom melittin to small unilamellar vesicles and large nonsonicated multilamellar bilayer membranes composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) was studied by means of circular dichroism, 31P-NMR and electrophoretic mobility. The melittin binding isotherm for small unilamellar vesicles (SUV) could be described by a partition equilibrium with Kp = (6 +/- 1).10(4) M-1. Electrostatic effects were taken into account by means of the Gouy-Chapman theory. Combining the partition equilibrium with the Gouy-Chapman analysis suggested an effective charge for melittin of Zp = 1.9, which is lower than the true electric charge of 5-6. The variation of the 31P-NMR signal of SUV showed the change in potential at the phosphodiester moiety of the lipid upon addition of melittin. This potential change was lower than that for an ion with an electrical charge of 5-6 and corresponded to a charge of 1.5. Electrophoretic mobility measurements with multilamellar vesicles confirmed the charge reduction effect. These experimental results show that the use of the simple Gouy-Chapman theory requires an effective electrical charge of the melittin which is lower than the formal charge.     

Concentration effect of cimetidine with POPC bilayer: a molecular dynamics simulation study

All-atom molecular dynamics simulations have been performed on cimetidine in the presence of a palmitoyloleoylphosphatidylcholine (POPC) bilayer. The free energy profile of a single cimetidine molecule passing across POPC bilayer displays a minimum at the interface of bilayer and water. Ten cimetidine molecules were inserted into POPC bilayer to obtain an 8 mol % drug model, and molecular dynamics results showed that cimetidine molecules reside at the polar region of POPC bilayer with sulphur atoms directing to the hydrophobic region. By comparing the one drug model with 8 mol % drug model, one can see that the central barrier to cross the membrane increases while the free energy in bulk water decreases, indicating that the ability of cimetidine passing across the POPC bilayer weakens at increased concentration. In addition, the free energy minimum shifts closer to the hydrophobic core. Our results indicate that with the increased drug concentration, it is more difficult for cimetidine to enter and pass across POPC bilayer.     

Stability of a melittin pore in a lipid bilayer: a molecular dynamics study

We have investigated the configuration and the stability of a single membrane pore bound by four melittin molecules and embedded in a fully hydrated bilayer lipid membrane. We used molecular dynamics simulations up to 5.8 ns. It is found that the initial tetrameric configuration decays with increasing time into a stable trimer and one monomer. This continuous transformation is accompanied by a lateral expansion of the aqueous pore exhibiting a final size comparable to experimental findings. The expansion-induced formation of an interface between the pore-lining acyl chains of the lipids and the pore water ("hydrophobic pore") is transformed into an energetically more favorable toroidal pore structure where some lipid heads are translocated from the rim to the central part of the interface ("hydrophilic pore"). The expansion of the pore is supported by the electrostatic repulsion among the alpha-helices. It is hypothesized that pore growth, and hence cell lysis, is induced by a melittin-mediated line tension of the pore.     

Protein-induced membrane disorder: a molecular dynamics study of melittin in a dipalmitoylphosphatidylcholine bilayer

A molecular dynamics simulation of melittin in a hydrated dipalmitoylphosphatidylcholine (DPPC) bilayer was performed. The 19, 000-atom system included a 72-DPPC phospholipid bilayer, a 26-amino acid peptide, and more than 3000 water molecules. The N-terminus of the peptide was protonated and embedded in the membrane in a transbilayer orientation perpendicular to the surface. The simulation results show that the peptide affects the lower (intracellular) layer of the bilayer more strongly than the upper (extracellular) layer. The simulation results can be interpreted as indicating an increased level of disorder and structural deformation for lower-layer phospholipids in the immediate vicinity of the peptide. This conclusion is supported by the calculated deuterium order parameters, the observed deformation at the intracellular interface, and an increase in fractional free volume. The upper layer was less affected by the embedded peptide, except for an acquired tilt relative to the bilayer normal. The effect of melittin on the surrounding membrane is localized to its immediate vicinity, and its asymmetry with respect to the two layers may result from the fact that it is not fully transmembranal. Melittin's hydrophilic C-terminus anchors it at the extracellular interface, leaving the N-terminus "loose" in the lower layer of the membrane. In general, the simulation supports a role for local deformation and water penetration in melittin-induced lysis. As for the peptide, like other membrane-embedded polypeptides, melittin adopts a significant 25 degree tilt relative to the membrane normal. This tilt is correlated with a comparable tilt of the lipids in the upper membrane layer. The peptide itself retains an overall helical structure throughout the simulation (with the exception of the three N-terminal residues), adopting a 30 degree intrahelical bend angle.     

The structure of melittin in lipid bilayer membranes

0.15 M inorganic phosphate dramatically increased the α-helix content of melittin in aqueous solution.When melittin interacted with egg yolk phosphatidylcholine liposomes in the absence of inorganic phosphate, it was converted to an α-helix rich form, as postulated by Dawson et al. (Dawson, C.R., Drake, A.F. Helliwell, J. and Hider, R.C. (1978) Biochim. Biophys. Acta 510, 75–86).     

Molecular dynamics simulation of melittin in a dimyristoylphosphatidylcholine bilayer membrane.

Molecular dynamics trajectories of melittin in an explicit dimyristoyl phosphatidylcholine (DMPC) bilayer are generated to study the details of lipid-protein interactions at the microscopic level. Melittin, a small amphipathic peptide found in bee venom, is known to have a pronounced effect on the lysis of membranes. The peptide is initially set parallel to the membrane-solution interfacial region in an alpha-helical conformation with unprotonated N-terminus. Solid-state nuclear magnetic resonance (NMR) and polarized attenuated total internal reflectance Fourier transform infrared (PATIR-FTIR) properties of melittin are calculated from the trajectory to characterize the orientation of the peptide relative to the bilayer. The residue Lys7 located in the hydrophobic moiety of the helix and residues Lys23, Arg24, Gln25, and Gln26 at the C-terminus hydrophilic form hydrogen bonds with water molecules and with the ester carbonyl groups of the lipids, suggesting their important contribution to the stability of the helix in the bilayer. Lipid acyl chains are closely packed around melittin, contributing to the stable association with the membrane. Calculated density profiles and order parameters of the lipid acyl chains averaged over the molecular dynamics trajectory indicate that melittin has effects on both layers of the membrane. The presence of melittin in the upper layer causes a local thinning of the bilayer that favors the penetration of water through the lower layer. The energetic factors involved in the association of melittin at the membrane surface are characterized using an implicit mean-field model in which the membrane and the surrounding solvent are represented as structureless continuum dielectric material. The results obtained by solving the Poisson-Bolztmann equation numerically are in qualitative agreement with the detailed dynamics. The influence of the protonation state of the N-terminus of melittin is examined. After 600 ps, the N-terminus of melittin is protonated and the trajectory is continued for 400 ps, which leads to an important penetration of water molecules into the bilayer. These observations provide insights into how melittin interacts with membranes and the mechanism by which it enhances their lysis.     

The structure of melittin in lipid bilayer membranes.

0.15 M inorganic phosphate dramatically increased the alpha-helix content of melittin in aqueous solution. When melittin interacted with egg yolk phosphatidylcholine liposomes in the absence of inorganic phosphate, it was converted to an alpha-helix rich form, as postulated by Dawson et al. (Dawson, C.R., Drake, A.F. Helliwell, J. and Hider, R.C. (1978) Biochim. Biophys. Acta 510, 75--86).     

An alamethicin channel in a lipid bilayer: molecular dynamics simulations

We present the results of 2-ns molecular dynamics (MD) simulations of a hexameric bundle of Alm helices in a 1-palmitoyl-2-oleoylphosphatidylcholine bilayer. These simulations explore the dynamic properties of a model of a helix bundle channel in a complete phospholipid bilayer in an aqueous environment. We explore the stability and conformational dynamics of the bundle in a phospholipid bilayer. We also investigate the effect on bundle stability of the ionization state of the ring of Glu18 side chains. If all of the Glu18 side chains are ionised, the bundle is unstable; if none of the Glu18 side chains are ionized, the bundle is stable. pKA calculations suggest that either zero or one ionized Glu18 is present at neutral pH, correlating with the stable form of the helix bundle. The structural and dynamic properties of water in this model channel were examined. As in earlier in vacuo simulations (Breed et al., 1996 .Biophys. J. 70:1643-1661), the dipole moments of water molecules within the pore were aligned antiparallel to the helix dipoles. This contributes to the stability of the helix bundle.     

Incorporation of highly purified melittin into phosphatidylcholine bilayer vesicles

Melittin free of phospholipase A2 was prepared. In the absence of salt this highly pure protein starts to aggregate in solution at a protein concentration of Cp greater than 10(-3) M. In high salt solution (2 M) aggregation starts at Cp greater than 10(-6) M. This was determined from the blue shift of the intrinsic fluorescence of the protein. Reinvestigation of the quenching behaviour clearly shows that self-aggregation cannot be deduced from quenching experiments using nitrate or 2,2,6,6-tetramethylpiperidine-1-oxyl as quencher. The incorporation of melittin into phosphatidylcholine bilayer vesicles was studied by fluorescence quenching and by energy-transfer experiments using 2- and 6-anthroyloxypalmitic acid as acceptor and peptide tryptophan as donor. Incorporation of melittin into small unilamellar vesicles was found to be reduced below the lipid phase transition temperature, Tt, whereas it incorporates and distributes more randomly above Tt. Cooling the temperature below Tt after incubation at T greater than Tt leads to a deeper incorporation of the peptide into the lipid bilayer due to electrostatic interaction between the lipid phosphate groups and the positively charged amino acids. This stabilizing effect is lost above Tt and melittin is extruded to the polar phase. Quenching experiments support this finding. EPR measurements clearly demonstrate that even in the presence of high amounts of melittin up to 10 mol% with respect to the lipid broadening of the phase transition curves was only observed with fatty acid spin labels, where the doxyl group is localized near the bilayer surface. The order degree of the inner part of the bilayer remains almost unchanged even in the presence of high melittin content.     

The interaction of bee melittin with lipid bilayer membranes

The influence of melittin and the related 8-26 peptide on the stability and electrical properties of bilayer lipid membranes is reported. Melittin, unlike the 8-26 peptide, has a dramatic influence on lipid membranes, causing rupture at dilute concentrations. The circular dichroism of melittin demonstrated that under physiological conditions, in water, melittin is in extended conformation, which is enhanced in aqueous ethanol. However in 'membrane-like' conditions it is essentially alpha-helical. Secondary structure predictions were used to locate possible alpha-helical nucleation centres and a model of melittin was built according to these predictions. It is postulated that melittin causes a wedge effect in membranes.     

Molecular dynamics simulation of lipid reorientation at bilayer edges

Understanding cellular membrane processes is critical for the study of events such as viral entry, neurotransmitter exocytosis, and immune activation. Supported lipid bilayers are commonly used to model these membrane processes experimentally. Despite the relative simplicity of such a system, many important structural and dynamic parameters are not experimentally observable with current techniques. Computational approaches allow the development of a high-resolution model of bilayer processes. We have performed molecular dynamics simulations of dimyristoylphosphatidylcholine (DMPC) bilayers to model the creation of bilayer gaps—a common process in bilayer patterning—and to analyze their structure and dynamics. We propose a model for gap formation in which the bilayer edges form metastable micelle-like structures on a nanosecond timescale. Molecules near edges structurally resemble lipids in ungapped bilayers but undergo small-scale motions more rapidly. These data suggest that lipids may undergo rapid local rearrangements during membrane fusion, facilitating the formation of fusion intermediates thought key to the infection cycle of viruses such as influenza, Ebola, and HIV.     

A long-lived amphiphilic fluorescent probe studied in POPC air-water monolayer and solution bilayer systems.

A novel amphiphilic fluorescent probe (Fluorazophore-L) with a strongly dipolar, nonionic azoalkane as headgroup and a palmitoyl tail has been synthesized and characterized. Pure Fluorazophore-L was found to be sufficiently amphiphilic to form stable air-water monolayers. An analysis of the surface pressure versus area suggests an area per molecule of about 34+/-2 A(2) at 29 mN m(-1). The partitioning into a lipid membrane model was quantified at the air-water interface by spreading 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) monolayers. Measurements with different molar fractions of Fluorazophore-L revealed a small but significant reduction of the mean area in the mixed monolayer. The excess free energy of mixing (-0.5+/-0.1 kT) indicated a weakly attractive interaction slightly above thermal energy, suggesting a good miscibility of the fluorescent probe within the lipid monolayer without major structural modifications. Spectroscopic measurements confirmed the incorporation of Fluorazophore-L into POPC vesicles. The fluorescence lifetime was very long (125+/-5 ns under air) with monoexponential fluorescence decays.     

Membrane protein dynamics versus environment: simulations of OmpA in a micelle and in a bilayer

The bacterial outer membrane protein OmpA is one of the few membrane proteins whose structure has been solved both by X-ray crystallography and by NMR. Crystals were obtained in the presence of detergent, and the NMR structure is of the protein in a detergent micelle. We have used 10 ns duration molecular dynamics simulations to compare the behaviour of OmpA in a detergent micelle and in a phospholipid bilayer. The dynamic fluctuations of the protein structure seem to be ca 1.5 times greater in the micelle environment than in the lipid bilayer. There are subtle differences between the nature of OmpA-detergent and OmpA-lipid interactions. As a consequence of the enhanced flexibility of the OmpA protein in the micellar environment, side-chain torsion angle changes are such as to lead to formation of a continuous pore through the centre of the OmpA molecule. This may explain the experimentally observed channel formation by OmpA.     

The transmembrane domain of Neu in a lipid bilayer: molecular dynamics simulations

The results of full-atom molecular dynamics simulations of the transmembrane domains (TMDs) of both native, and Glu⁶⁶⁴-mutant (either protonated or unprotonated) Neu in an explicit fully hydrated dimyristoylphosphatidylcholine (DMPC) lipid bilayer are presented. For the native TMD peptide, a 10.05 ns trajectory was collected, while for the mutant TMD peptides 5.05 ns trajectories were collected for each. The peptides in all three simulations display stable predominantly -helical hydrogen bonding throughout the trajectories. The only significant exception occurs near the C-terminal end of the native and unprotonated mutant TMDs just outside the level of the lipid headgroups, where -helical hydrogen bonding develops, introducing a kink in the backbone structure. However, there is no indication of the formation of a bulge within the hydrophobic region of either native or mutant peptides. Over the course of the simulation of the mutant peptide, it is found that a significant number of water molecules penetrate the hydrophobic region of the surrounding lipid molecules, effectively hydrating Glu⁶⁶⁴. If the energy cost of such water penetration is significant enough, this may be a factor in the enhanced dimerization affinity of Glu⁶⁶⁴-mutant Neu.     

Distribution of amino acids in a lipid bilayer from computer simulations

We have calculated the distribution in a lipid bilayer of small molecules mimicking 17 natural amino acids in atomistic detail by molecular dynamics simulation. We considered both charged and uncharged forms for Lys, Arg, Glu, and Asp. The results give detailed insight in the molecular basis of the preferred location and orientation of each side chain as well the preferred charge state for ionizable residues. Partitioning of charged and polar side chains is accompanied by water defects connecting the side chains to bulk water. These water defects dominate the energetic of partitioning, rather than simple partitioning between water and a hydrophobic phase. Lys, Glu, and Asp become uncharged well before reaching the center of the membrane, but Arg may be either charged or uncharged at the center of the membrane. Phe has a broad distribution in the membrane but Trp and Tyr localize strongly to the interfacial region. The distributions are useful for the development of coarse-grained and implicit membrane potentials for simulation and structure prediction. We discuss the relationship between the distribution in membranes, bulk partitioning to cyclohexane, and several amino acid hydrophobicity scales.     

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.     

Interpretation of fluctuation spectra in lipid bilayer simulations

Atomic resolution and coarse-grained simulations of dimyristoylphosphatidylcholine lipid bilayers were analyzed for fluctuations perpendicular to the bilayer using a completely Fourier-based method. We find that the fluctuation spectrum of motions perpendicular to the bilayer can be decomposed into just two parts: 1), a pure undulation spectrum proportional to q⁻⁴ that dominates in the small-q regime; and 2), a molecular density structure factor contribution that dominates in the large-q regime. There is no need for a term proportional to q⁻² that has been postulated for protrusion fluctuations and that appeared to have been necessary to fit the spectrum for intermediate q. We suggest that earlier reports of such a term were due to the artifact of binning and smoothing in real space before obtaining the Fourier spectrum. The observability of an intermediate protrusion regime from the fluctuation spectrum is discussed based on measured and calculated material constants.     

Binding and folding of melittin in the presence of ganglioside GM1 micelles

In this work we examine the binding and folding of the membrane-active peptide, melittin in the presence of ganglioside GM1 micelle. The membrane bilayer is capable of inducing folding to small proteins and peptides upon binding. Using two-dimensional NMR techniques we have shown that at low concentration, GM1 micelle is able to induce an extended helical conformation to MLT. The pulsed-field gradient diffusion NMR study indicates that the peptide partition into GM1 micelle along with about 32% binding. While looking for the binding between MLT and GM1 using saturation transfer difference NMR spectroscopy, Val5, Leu9, Thr11, Ile17, Ser18, and Trp19 have been identified as the residues that are in close proximity to GM1 micelles.