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.     

Structure and dynamics of an amphiphilic peptide in a lipid bilayer: a molecular dynamics study

A molecular dynamics simulation of a simple model membrane system composed of a single amphiphilic helical peptide (ace-K2GL16K2A-amide) in a fully hydrated 1,2-dimyristoyl-sn-glycero-3-phosphocholine bilayer was performed for a total of 1060 ps. The secondary structure of the peptide and its stability were described in terms of average dihedral angles, phi and psi, and the C alpha torsion angles formed by backbone atoms; by the average translation per residue along the helix axis; and by the intramolecular peptide hydrogen bonds. The results indicated that residues 6 through 15 remain in a stable right-handed alpha-helical conformation, whereas both termini exhibit substantial fluctuations. A change in the backbone dihedral angles for residues 16 and 17 is accompanied by the loss of two intramolecular hydrogen bonds, leading to a local but long-lived disruption of the helix. The dynamics of the peptide was characterized in terms of local and global helix motions. The local motions of the N-H bond angles were described in terms of the autocorrelation functions of P2[cos thetaNH(t, t + tau)] and reflected the different degrees of local peptide order as well as a variation in time scale for local motions. The chi1 and chi2 dihedral angles of the leucine side chains underwent frequent transitions between potential minima. No connection between the side-chain positions and their mobility was observed, however. In contrast, the lysine side chains displayed little mobility during the simulation. The global peptide motions were characterized by the tilting and bending motions of the helix. Although the peptide was initially aligned parallel to the bilayer normal, during the simulation it was observed to tilt away from the normal, reaching an angle of approximately 25 degrees by the end of the simulation. In addition, a slight bend of the helix was detected. Finally, the solvation of the peptide backbone and side-chain atoms was also investigated.     

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.     

Cooperative antimicrobial action of melittin on lipid membranes: A coarse-grained molecular dynamics study

We conducted a series of coarse-grained molecular dynamics (CG-MD) simulations to investigate the complicated actions of melittin, which is an antimicrobial peptide (AMP) derived from honey bee venom, on a lipid membrane. To accurately simulate the AMP action, we developed and used a protein CG model as an extension of the pSPICA force field (FF), which was designed to reproduce several thermodynamic quantities and structural properties. At a low peptide-to-lipid (P/L) ratio (1/102), no defect was detected. At P/L = 1/51, toroidal pore formation was observed due to collective insertion of multiple melittin peptides from the N-termini. The pore formation was initiated by a local increase in membrane curvature in the vicinity of the peptide aggregate. At a higher P/L ratio (1/26), two more modes were detected, seemingly not controlled by the P/L ratio but by a local arrangement of melittin peptides: 1. Pore formation accompanied by lipid extraction by melittin peptides:a detergent-like mechanism. 2. A rapidly formed large pore in a significantly curved membrane: bursting. Thus, we observed three pore formation modes (toroidal pore formation, lipid extraction, and bursting) depending on the peptide concentration and local arrangement. These observations were consistent with experimental observations and hypothesized melittin modes. Through this study, we found that the local arrangements and population of melittin peptides and the area expansion rate by membrane deformation were key to the initiation of and competition among the multiple pore formation mechanisms.     

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).     

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.     

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).     

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.     

Coarse-grained molecular dynamics of tetrameric transmembrane peptide bundles within a lipid bilayer

The conformations of model transmembrane peptides are studied to understand the structural and dynamical aspects of tetrameric bundles using a series of coarse grain (CG) molecular dynamics (MD) simulations since membrane proteins play a crucial role in cell function. In this work, two different amphipathic models have been constructed using similar hydrophobic/hydrophilic characteristics with two structurally distinct morphologies to evaluate the effect of roughness and hydrophilic topology on the structure of tetrameric bundles, one class that forms an ion-channel and one class that does not. Free energy calculations of typical amphipathic peptide topologies show that using a relatively smooth surface morphology allows for a stable conformation of the tetramer bundle in a diamond formation. However, the model with side chains attached to the core in order to roughen the surface has a stable square tetramer bundle which is consistent with experimental data and all-atom (AA) MD simulations. Comparisons of the CG simulations with AA MD simulations are in reasonable agreement with the formation of tetrameric homo-oligomers, partitioning within the lipid bilayer and tilt angle with respect to the bilayer normal. We concluded that a square or diamond shape tetrameric homo-oligomers could be stabilized by rational design of the peptide morphology and topology of the surface, thus allowing us to tune the permeability of the bundle or channel.     

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.     

A comparative molecular dynamics analysis of the amyloid beta-peptide in a lipid bilayer

Because the amyloid β-peptide (Aβ) functions as approximately half of the transmembrane domain of the amyloid precursor protein and interaction of Aβ with membranes is proposed to result in neurotoxicity, the association of Aβ with membranes likely is important in the etiology of Alzheimer’s disease. Atomic details of the interaction of Aβ with membranes are not accessible with most experimental techniques, but computational methods can provide this information. Here, we present the results of ten 100-ns molecular dynamics (MD) simulations of the 40-residue amyloid β-peptide (Aβ₄₀) embedded in a dipalmitoylphosphatidylcholine (DPPC) bilayer. The present study examines the effects of insertion depth, protonation state of key residues, and ionic strength on Aβ₄₀ in a DPPC bilayer. In all cases, a portion of the peptide remained embedded in the bilayer. In the case of deeper insertion depth, Aβ₄₀ adopted a near-transmembrane orientation, drawing water molecules into the bilayer to associate with its charged amino acids. In the case of shallower insertion, the most widely-accepted construct, the peptide associated strongly with the membrane-water interface and the phosphatidylcholine headgroups of the bilayer. In most cases, significant disordering of the extracellular segment of the peptide was observed, and the brief appearance of a β-strand was noted in one case. Our results compare well with a variety of experimental and computational findings. From this study, we conclude that Aβ associated with membranes is dynamic and capable of adopting a number of conformations, each of which may have significance in understanding the progression of Alzheimer’s disease.     

A molecular dynamics study of the bee venom melittin in aqueous solution, in methanol, and inserted in a phospholipid bilayer

The structural properties of melittin, a small amphipathic peptide found in the bee venom, are investigated in three different environments by molecular dynamics simulation. Long simulations have been performed for monomeric melittin solvated in water, in methanol, and shorter ones for melittin inserted in a dimyristoylphosphatidylcholine bilayer. The resulting trajectories were analysed in terms of structural properties of the peptide and compared to the available NMR data. While in water and methanol solution melittin is observed to partly unfold, the peptide retains its structure when embedded in a lipid bilayer. The latter simulation shows good agreement with the experimentally derived ³J-coupling constants. Generally, it appears that higher the stability of the helical conformation of melittin, lower is the dielectric permittivity of the environment. In addition, peptide-lipid interactions were investigated showing that the C-terminus of the peptide provides an anchor to the lipid bilayer by forming hydrogen bonds with the lipid head groups.     

Dynamical properties of a hydrated lipid bilayer from a multinanosecond molecular dynamics simulation

A fully hydrated dimiristoylphosphatidylcholine (DMPC) bilayer has been studied by a molecular dynamics simulation. The system, which consisted of 64 DMPC molecules and 1792 water molecules, was run in the NVE ensemble at a temperature of 333 K for a total of 10 ns. The resulting trajectory was used to analyze structural and dynamical quantities. The electron density, bilayer spacing, and order parameters (S(CD)), based on the AMBER forcefield and SPCE water model are in good agreement with previous calculations and experimental data. The simulation reveals evidence for two types of lateral diffusive behavior: cage hopping and that of a two-dimensional liquid. The lateral diffusion coefficient is 8 x 10(-8) cm(2)/s. We characterize the rotational motion, and find that the lipid tail rotation (D(rot_tail) = -0.04 rad(2)/ns) is slower then the head group rotation (D(rot_hg) = 2.2 rad(2)/ns), which is slower than the overall in plane (D(rot) = 3.2 rad(2)/ns) for the lipid molecule.     

Interfacing molecular dynamics and macro-scale simulations for lipid bilayer vesicles

A continuum-level model for a giant unilamellar vesicle (GUV) is bridged to a corresponding atomistic model of a dimyristoylphosphatidylcholine (DMPC) bilayer at various cholesterol concentrations via computation of the bulk modulus. The bulk modulus and other microscopically determined parameters are passed to a continuum-level model operating in time- and length-scales orders of magnitude beyond that which is accessible by atomistic-level simulation. The continuum-level simulation method used is the material point method (MPM), and the particular variation used here takes advantage of the spherical nature of many GUVs. An osmotic pressure gradient due to a solvent concentration change is incorporated into the continuum-level simulation, resulting in osmotic swelling of the vesicle. The model is then extended to treat mixtures of DMPC and cholesterol, where small domains of different composition are considered.     

Structural properties of a highly polyunsaturated lipid bilayer from molecular dynamics simulations.

The structure of a fully hydrated mixed (saturated/polyunsaturated) chain lipid bilayer in the biologically relevant liquid crystalline phase has been examined by performing a molecular dynamics study. The model membrane, a 1-stearoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine (SDPC, 18:0/22:6 PC) lipid bilayer, was investigated at constant (room) temperature and (ambient) pressure, and the results obtained in the nanosecond time scale reproduced quite well the available experimental data. Polyunsaturated fatty acids are found in high concentrations in neuronal and retinal tissues and are essential for the development of human brain function. The docosahexaenoic fatty acid, in particular, is fundamental for the proper function of the visual receptor rhodopsin. The lipid bilayer order has been investigated through the orientational order parameters. The water-lipid interface has been explored thoroughly in terms of its dimensions and the organization of the different components. Several types of interactions occurring in the system have been analyzed, specifically, the water-hydrocarbon chain, lipid-lipid and lipid-water interactions. The distribution of dihedral angles along the chains and the molecular conformations of the polyunsaturated chain of the lipids have also been studied. Special attention has been focused on the microscopic (molecular) origin of the effects of polyunsaturations on the different physical properties of membranes.     

Molecular dynamics study of aquaporin-1 water channel in a lipid bilayer

The aquaporin-1 water channel was modeled in a palmitoyl-oleoyl-phosphatidyl-choline lipid bilayer, by means of molecular dynamics simulations. Interaction of the protein with the membrane and inter-monomer interactions were analyzed. Structural features of the channel important for its biological function, including the Asn-Pro-Ala (NPA) motifs, and the diffusion of water molecules into the channels, were investigated. Simulations revealed the formation of single file water inside the channels for certain relative positions of the NPA motifs.     

Molecular dynamics study of MscL interactions with a curved lipid bilayer

Mechanosensitivity is a ubiquitous sensory mechanism found in living organisms. The simplest known mechanotransducing mechanism is found in bacteria in the form of the mechanosensitive membrane channel of large conductance, MscL. This channel has been studied extensively using a variety of methods at a functional and structural level. The channel is gated by membrane tension in the lipid bilayer alone. It serves as a safety valve protecting bacterial cells against hypoosmotic shock. MscL of Escherichia coli embedded in bilayers composed of asymmetric amounts of single-tailed and double-tailed lipids has been shown to gate spontaneously, even in the absence of membrane tension. To gain insight into the effect of the lipid membrane composition and geometry on MscL structure, a fully solvated, all-atom model of MscL in a stress-free curved bilayer composed of double- and single-tailed lipids was studied using a 9.5-ns molecular dynamics simulation. The bilayer was modeled as a domed structure accommodating the asymmetric composition of the monolayers. During the course of the simulation a spontaneous restructuring of the periplasmic loops occurred, leading to interactions between one of the loops and phospholipid headgroups. Previous experimental studies of the role of the loops agree with the observation that opening starts with a restructuring of the periplasmic loop, suggesting an effect of the curved bilayer. Because of limited resources, only one simulation of the large system was performed. However, the results obtained suggest that through the geometry and composition of the bilayer the protein structure can be affected even on short timescales.     

Cholesterol effects on a mixed-chain phosphatidylcholine bilayer: a molecular dynamics simulation study

A molecular dynamics simulation of a mono-cis-unsaturated 1-palmitoyl-2-oleoyl-phosphatidylcholine bilayer containing approximately 22 mol% of cholesterol (POPC-Chol) was carried out for 15 ns. An 8-ns trajectory was analysed to determine the effects of Chol on the membrane properties and compare it with that on the fully saturated 1,2-dimyristoyl-phosphatidylcholine bilayer containing approximately 22 mol% of Chol (DMPC-Chol). The study suggests that the experimentally observed weaker effect of Chol on the POPC than DMPC bilayer might result from a different vertical localisation of the Chol hydroxyl group (OH-Chol) in both bilayers: in the POPC-Chol bilayer, OH-Chol is placed approximately 3 A higher in the bilayer interface than in the DMPC-Chol bilayer. Because of the rigid cis double bond in the beta-chain of POPC, Chol fits worse to the POPC-Chol membrane environment and is pushed up, in effect all Chol ring atoms are, on average, located above the double bond. Both in mono-cis-unsaturated and fully saturated PC bilayers, Chol induces stronger van der Waals interactions among the chains, whereas its interactions with the chains are weak. In contrast to DMPC, the smooth alpha-face of the Chol ring lowers the order of POPC chains, whereas the rough beta-face increases the order.     

Effect of hexafluoroisopropanol alcohol on the structure of melittin: a molecular dynamics simulation study

The molecular mechanism by which HFIP stabilizes the alpha-helical structure of peptides is not well understood. In the present study, we use melittin as a model to gain insight into the details of the atomistic interactions of HFIP with the peptide. We have performed extensive comparative molecular dynamics simulations (up to 100 nsec) in the absence and in the presence of HFIP. In agreement with recent NMR experiments, the simulations show rapid loss of tertiary structure in water at pH 2 but much higher helicity in 35% HFIP. The MD simulations also indicate that melittin adopts a highly dynamic global structure in 35% HFIP solution with two alpha-helical segments sampling a wide range of angular orientations. The analysis of the HFIP distribution shows the tendency of HFIP to aggregate around the peptide, increasing the local cosolvent concentration to more than two times that in the bulk concentration. The correlation of local peptide structure with HFIP coating suggests that displacement of water at the peptide surface is the main contribution of HFIP in stabilizing the secondary structure of melittin. Finally, a stabilizing effect promoted by the presence of counter-ions was also observed in the simulations.