Molecular dynamics simulation of a phospholipid membrane

We present the results of molecular dynamics (MD) simulations of a phospholipid membrane in water, including full atomic detail. The goal of the simulations was twofold: first we wanted to set up a simulation system which is able to reproduce experimental results and can serve as a model membrane in future simulations. This goal being reached it is then further possible to gain insight in to those properties that are experimentally more difficult to access. The system studied is dipalmitoylphosphatidylcholine/water, consisting of 5408 atoms. Using original force field parameters the membrane turned out to approach a gel-like state. With slight changes of the parameters, the system adopted a liquid-crystalline state. Separate 80 ps runs were performed on both the gel and liquid-crystalline systems. Comparison of MD results with reliable experimental data (bilayer repeat distance, surface area per lipid, tail order parameters, atom distributions) showed that our simulations, especially the one in the liquid-crystalline phase, can serve as a realistic model for a phospholipid membrane. Further analysis of the trajectories revealed valuable information on various properties. In the liquid-crystalline phase, the interface turns out to be quite diffuse, with water molecules penetrating into the bilayer to the position of the carbonyl groups. The 10–90% width of the interface turns out to be 1.3 nm and the width of the hydrocarbon interior 3.0 nm. The headgroup dipoles are oriented at a small angle with respect to the bilayer plane. The resulting charge distribution is almost completely cancelled by the water molecules. The electron density distribution shows a large dip in the middle of the membrane. In this part the tails are more flexible. The mean life time between dihedral transitions is 20 ps. The average number of gauche angles per tail is 3.5. The occurrence of kinks is not a significant feature.Abbreviations MD molecular dynamics - DPPC dipalmitoylphosphatidylcholine - SPC simple point charges - DPPE dipalmitoylphosphatidylethanolamine Correspondence to: H. J. C. Berendsen     

Molecular dynamics simulations of hydrophilic pores in lipid bilayers

Hydrophilic pores are formed in peptide free lipid bilayers under mechanical stress. It has been proposed that the transport of ionic species across such membranes is largely determined by the existence of such meta-stable hydrophilic pores. To study the properties of these structures and understand the mechanism by which pore expansion leads to membrane rupture, a series of molecular dynamics simulations of a dipalmitoylphosphatidylcholine (DPPC) bilayer have been conducted. The system was simulated in two different states; first, as a bilayer containing a meta-stable pore and second, as an equilibrated bilayer without a pore. Surface tension in both cases was applied to study the formation and stability of hydrophilic pores inside the bilayers. It is observed that below a critical threshold tension of approximately 38 mN/m the pores are stabilized. The minimum radius at which a pore can be stabilized is 0.7 nm. Based on the critical threshold tension the line tension of the bilayer was estimated to be approximately 3 x 10(-11) N, in good agreement with experimental measurements. The flux of water molecules through these stabilized pores was analyzed, and the structure and size of the pores characterized. When the lateral pressure exceeds the threshold tension, the pores become unstable and start to expand causing the rupture of the membrane. In the simulations the mechanical threshold tension necessary to cause rupture of the membrane on a nanosecond timescale is much higher in the case of the equilibrated bilayers, as compared with membranes containing preexisting pores.     

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.     

Molecular Dynamic Simulation of Transmembrane Pore Growth

A molecular dynamic approach was applied for simulation of dynamics of pore formation and growth in a phospholipid bilayer in the presence of an external electric field. Processing the simulation results permitted recovery of the kinetic coefficients used in the Einstein–Smoluchowski equation describing the dynamics of pore evolution. Two different models of the bilayer membrane were considered: membrane consisting of POPC and POPE lipids. The simulations permitted us to find nonempirical values of the pore energy parameters, which are compared with empirical values. It was found that the parameters are sensitive to membrane type.     

Structure and dynamics of the pore-lining helix of the nicotinic receptor: MD simulations in water, lipid bilayers, and transbilayer bundles

Multiple nanosecond duration molecular dynamics simulations on the pore-lining M2 helix of the nicotinic acetylcholine receptor reveal how its structure and dynamics change as a function of environment. In water, the M2 helix partially unfolds to form a molecular hinge in the vicinity of a central Leu residue that has been implicated in the mechanism of ion channel gating. In a phospholipid bilayer, either as a single transmembrane helix, or as part of a pentameric helix bundle, the M2 helix shows less flexibility, but still exhibits a kink in the vicinity of the central Leu. The single M2 helix tilts relative to the bilayer normal by 12 degrees, in agreement with recent solid state NMR data (Opella et al., Nat Struct Biol 6:374-379, 1999). The pentameric helix bundle, a model for the pore domain of the nicotinic receptor and for channels formed by M2 peptides in a bilayer, is remarkably stable over a 2-ns MD simulation in a bilayer, provided one adjusts the pK(A)s of ionizable residues to their calculated values (when taking their environment into account) before starting the simulation. The resultant transbilayer pore shows fluctuations at either mouth which transiently close the channel. Proteins 2000;39:47-55.     

Proton transport across transient single-file water pores in a lipid membrane studied by molecular dynamics simulations.

To test the hypothesis that water pores in a lipid membrane mediate the proton transport, molecular dynamic simulations of a phospholipid membrane, in which the formation of a water pore is induced, are reported. The probability density of such a pore in the membrane was obtained from the free energy of formation of the pore, which was computed from the average force needed to constrain the pore in the membrane. It was found that the free energy of a single file of water molecules spanning the bilayer is 108(+/-10) kJ/mol. From unconstrained molecular dynamic simulations it was further deduced that the nature of the pore is very transient, with a mean lifetime of a few picoseconds. The orientations of water molecules within the pore were also studied, and the spontaneous translocation of a turning defect was observed. The combined data allowed a permeability coefficient for proton permeation across the membrane to be computed, assuming that a suitable orientation of the water molecules in the pore allows protons to permeate the membrane relatively fast by means of a wirelike conductance mechanism. The computed value fits the experimental data only if it is assumed that the entry of the proton into the pore is not rate limiting.     

Molecular dynamics simulations of pentapeptides at interfaces: salt bridge and cation-pi interactions

Peptide-membrane interactions are important for understanding the binding, partitioning, and folding of membrane proteins; the activity of antimicrobial and fusion peptides; and a number of other processes. We describe molecular dynamics simulations (10-25 ns) of two pentapeptides Ace-WLXLL (with X = Arg or Lys side chain) (White, S. H., and Wimley, W.C. (1996) Nat. Struct. Biol. 3, 842-848) in water and three different membrane mimetic systems: (i) a water/cyclohexane interface, (ii) water-saturated octanol, and (iii) a solvated dioleoylphosphatidylcholine bilayer. A salt bridge is found between the protonated Arg or Lys side chains with the carboxyl terminus at the three interfaces. In water/cyclohexane, the salt bridge is most exposed to the water phase and least stable. In water/octanol and the lipid bilayer systems, the salt bridge once formed persists throughout the simulations. In the lipid bilayer, the salt bridge is more stable when the peptide penetrates deeper into the bilayer. In one of two peptides, a cation-pi interaction between the Arg and the Trp side chains is stable in the lipid bilayer for about 15 ns before breaking. In all cases, the conformations of the peptides are restricted by their presence at the interface and can be assigned to a few major conformational clusters. Side chains facing the water phase are most mobile. In the lipid bilayer, the peptides remain in the interface area, where they overlap with the carbonyl area of the lipid bilayer and perturb the local density profile of the bilayer. The tryptophan side chain remains in the water-lipid interface, where it interacts with the lipid choline group and forms hydrogen bonds with the ester carbonyl of the lipid and with water in the interface.     

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.     

Lipid hydration and mobility: an interplay between fluorescence solvent relaxation experiments and molecular dynamics simulations

Fluorescence solvent relaxation experiments are based on the characterization of time-dependent shifts in the fluorescence emission of a chromophore, yielding polarity and viscosity information about the chromophore’s immediate environment. A chromophore applied to a phospholipid bilayer at a well-defined location (with respect to the z-axis of the bilayer) allows monitoring of the hydration and mobility of the probed segment of the lipid molecules. Specifically, time-resolved fluorescence experiments, fluorescence quenching data and molecular dynamic (MD) simulations show that 6-lauroyl-2-dimethylaminonaphthalene (Laurdan) probes the hydration and mobility of the sn-1 acyl groups in a phosphatidylcholine bilayer. The time-dependent fluorescence shift (TDFS) of Laurdan provides information on headgroup compression and expansion induced by the addition of different amounts of cationic lipids to phosphatidylcholine bilayers. Those changes were predicted by previous MD simulations. Addition of truncated oxidized phospholipids leads to increased mobility and hydration at the sn-1 acyl level. This experimental finding can be explained by MD simulations, which indicate that the truncated chains of the oxidized lipid molecules are looping back into aqueous phase, hence creating voids below the glycerol level. Fluorescence solvent relaxation experiments are also useful in understanding salt effects on the structure and dynamics of lipid bilayers. For example, such experiments demonstrate that large anions increase hydration and mobility at the sn-1 acyl level of phosphatidylcholine bilayers, an observation which could not be explained by standard MD simulations. If polarizability is introduced into the applied force field, however, MD simulations show that big soft polarizable anions are able to interact with the hydrophilic/hydrophobic interface of the lipid bilayer, penetrating to the level probed by Laurdan, and that they expand and destabilize the bilayer making it more hydrated and mobile.     

Diffusional dynamics of an active rhodamine-labeled 1,4-dihydropyridine in sarcolemmal lipid multibilayers.

A "membrane bilayer pathway" model, involving ligand partition into the bilayer, lateral diffusion, and receptor binding has been invoked to describe the 1,4-dihydropyridine (DHP) calcium channel antagonist receptor binding mechanism. In an earlier study (Chester et al. 1987. Biophys. J. 52:1021-1030), the diffusional component of this model was examined using an active fluorescence labeled DHP calcium channel antagonist, nisoldipine-lissamine rhodamine B (Ns-R), in purified cardiac sarcolemmal (CSL) lipid multibilayers. Diffusion coefficient measurements on membrane-bound drug and phospholipid at maximum bilayer hydration yielded similar values (3.8 x 10(-8) cm2/s). However, decreases in bilayer hydration resulted in dramatically reduced diffusion coefficient values for both probes with substantially greater impact on Ns-R diffusion. These data suggested that hydration dependent diffusional differences could be a function of relative probe location along the bilayer normal. In this communication, we have addressed the relative effect of the rhodamine substituent on Ns-R diffusion complex by examining the diffusional dynamics of free rhodamine B under the same conditions used to evaluate Ns-R complex and phospholipid diffusion. X-ray diffraction studies were performed to determine the Ns-R location in the membrane and model the CSL lipid bilayer profile structure to give a rationale for the differences in probe diffusional dynamics as a function of interbilayer water space.     

Molecular dynamics of structural transitions and intercalation in DNA

The large-scale flexibility of DNA and the intercalation of actinomycin D have been studied by computer simulation using molecular dynamics. The stretching and unwinding of B and Z forms of DNA and intercalation in B-DNA were examined through molecular dynamics simulations, and the energetics of transitions were calculated by the conformational energy minimization method. The principal results of this research are as follows: (1) A dynamic conformational pathway is presented for longitudinal stretching and unwinding of the double helix to open an intercalation site. (2) Large-scale transitions are possible in both B and Z forms of DNA through a conformationally allowed kinetic pathway. (3) The stretching and untwisting of a 5′(CG)3′ step is energetically more favorable than for a GC step in B-DNA. (4) The formation of an adjacent second cavity in B-DNA requires larger energy than the formation of the first cavity, affirming the neighbor-exclusion principle of intercalation. (5) Docking an intercalated actinomycin D in the stretched structure is shown to be geometrically and energetically feasible.     

Life Cycle of an Electropore: Field-Dependent and Field-Independent Steps in Pore Creation and Annihilation

Electropermeabilization, an electric field-induced modification of the barrier functions of the cell membrane, is widely used in laboratories and increasingly in the clinic; but the mechanisms and physical structures associated with the electromanipulation of membrane permeability have not been definitively characterized. Indirect experimental observations of electrical conductance and small molecule transport as well as molecular dynamics simulations have led to models in which hydrophilic pores form in phospholipid bilayers with increased probability in the presence of an electric field. Presently available methods do not permit the direct, nanoscale examination of electroporated membranes that would confirm the existence of these structures. To facilitate the reconciliation of poration models with the observed properties of electropermeabilized lipid bilayers and cell membranes, we propose a scheme for characterizing the stages of electropore formation and resealing. This electropore life cycle, based on molecular dynamics simulations of phospholipid bilayers, defines a sequence of discrete steps in the electric field-driven restructuring of the membrane that leads to the formation of a head group-lined, aqueous pore and then, after the field is removed, to the dismantling of the pore and reassembly of the intact bilayer. Utilizing this scheme we can systematically analyze the interactions between the electric field and the bilayer components involved in pore initiation, construction and resealing. We find that the pore creation time depends strongly on the electric field gradient across the membrane interface and that the pore annihilation time is at least weakly dependent on the magnitude of the pore-initiating electric field and, in general, much longer than the pore creation time.     

Force measurements on myelin basic protein adsorbed to mica and lipid bilayer surfaces done with the atomic force microscope.

The mechanical and adhesion properties of myelin basic protein (MBP) are important for its function, namely the compaction of the myelin sheath. To get more information about these properties we used atomic force microscopy to study tip-sample interaction of mica and mixed dioleoylphosphatidylserine (DOPS) (20%)/egg phosphatidylcholine (EPC) (80%) lipid bilayer surfaces in the absence and presence of bovine MBP. On mica or DOPS/EPC bilayers a short-range repulsive force (decay length 1.0-1.3 nm) was observed during the approach. The presence of MBP always led to an attractive force between tip and sample. When retracting the tip again, force curves on mica and on lipid layers were different. While attached to the mica surface, the MBP molecules exhibited elastic stretching behavior that agreed with the worm-like chain model, yielding a persistence length of 0.5 +/- 0.25 nm and an average contour length of 53 +/- 19 nm. MBP attached to a lipid bilayer did not show elastic stretching behavior. This shows that the protein adopts a different conformation when in contact with lipids. The lipid bilayer is strongly modified by MBP attachment, indicating formation of MBP-lipid complexes and possibly disruption of the original bilayer structure.     

Studying the structure of single-component ceramide bilayers with molecular dynamics simulations using different force fields

Molecular dynamics simulations are performed on a lipid bilayer that consists of ceramide NS 24:0 in an attempt to examine several structural and physicochemical properties of the specific system. The simulations are carried out with five different force fields (OPLS, GROMOS, BERGER, CHARMM and GAFF) in order to evaluate and compare their performance in modelling lipid systems that contain ceramides. The examined properties include bilayer thickness, chain tilt, density profiles, order parameters, chain conformation, area per lipid and (intermolecular or intramolecular) hydrogen bonding between the head groups. Special focus is given to the lateral lipid arrangement. To this purpose, a method is proposed that utilises the radial distribution functions of the alkyl chains to derive quantitative information about the lateral lipid packing. In most cases, all force fields lead to similar results. For a few properties (e.g. intramolecular hydrogen bonding), there is some discrepancy between the force fields but the lack of respective experimental data does not allow an unambiguous conclusion on which force field is the most reliable.     

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.     

Interactions of liquid crystal-forming molecules with phospholipid bilayers studied by molecular dynamics simulations

Recent experiments have shown that liquid crystals can be used to image mammalian cell membranes and to amplify structural reorganization in phospholipid-laden liquid crystal-aqueous interfaces. In this work, molecular dynamics simulations were employed to explore the interactions between commonly used liquid crystal-forming molecules and phospholipid bilayers. In particular, umbrella sampling was used to obtain the potential of mean force of 4-cyano-4'-pentylbiphenyl (5CB) and 4'-(3,4-difluor-phenyl)-4-pentyl-bicylohexyl (5CF) molecules partitioning into a dipalmitoylphosphatidylcholine bilayer. In addition, results of simulations are presented for systems consisting of a fully hydrated bilayer with 5CB or 5CF molecules at the lowest (4.5 mol %) and highest (20 mol %) concentrations used in recent laboratory experiments. It is found that mesogens preferentially partition from the aqueous phase into the membrane; the potential of mean force exhibits highly favorable free energy differences for partitioning (-18 k(B)T for 5CB and -26 k(B)T for 5CF). The location and orientation of mesogens associated with the most stable free energies in umbrella sampling simulations of dilute systems were found to be consistent with those observed in liquid-crystal-rich bilayers. It is found that the presence of mesogens in the bilayer enhances the order of lipid acyl tails, and changes the spatial and orientational arrangement of lipid headgroup atoms. These effects are more pronounced at higher liquid-crystal concentrations. In comparing the behavior of 5CB and 5CF, a stronger spatial correlation (i.e., possibly leading to aggregation) is observed between 5CB molecules within a bilayer than between 5CF molecules. Also, the range of molecular orientations and positions along the bilayer normal is larger for 5CB molecules. At the same time, 5CF molecules were found to bind more strongly to lipid headgroups, thereby slowing the lateral motion of lipid molecules.     

Influence of hydrophobic mismatch on structures and dynamics of gramicidin a and lipid bilayers

Gramicidin A (gA) is a 15-amino-acid antibiotic peptide with an alternating L-D sequence, which forms (dimeric) bilayer-spanning, monovalent cation channels in biological membranes and synthetic bilayers. We performed molecular dynamics simulations of gA dimers and monomers in all-atom, explicit dilauroylphosphatidylcholine (DLPC), dimyristoylphosphatidylcholine (DMPC), dioleoylphosphatidylcholine (DOPC), and 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayers. The variation in acyl chain length among these different phospholipids provides a way to alter gA-bilayer interactions by varying the bilayer hydrophobic thickness, and to determine the influence of hydrophobic mismatch on the structure and dynamics of both gA channels (and monomeric subunits) and the host bilayers. The simulations show that the channel structure varied little with changes in hydrophobic mismatch, and that the lipid bilayer adapts to the bilayer-spanning channel to minimize the exposure of hydrophobic residues. The bilayer thickness, however, did not vary monotonically as a function of radial distance from the channel. In all simulations, there was an initial decrease in thickness within 4–5 Å from the channel, which was followed by an increase in DOPC and POPC or a further decrease in DLPC and DMPC bilayers. The bilayer thickness varied little in the monomer simulations—except one of three independent simulations for DMPC and all three DLPC simulations, where the bilayer thinned to allow a single subunit to form a bilayer-spanning water-permeable pore. The radial dependence of local lipid area and bilayer compressibility is also nonmonotonic in the first shell around gA dimers due to gA-phospholipid interactions and the hydrophobic mismatch. Order parameters, acyl chain dynamics, and diffusion constants also differ between the lipids in the first shell and the bulk. The lipid behaviors in the first shell around gA dimers are more complex than predicted from a simple mismatch model, which has implications for understanding the energetics of membrane protein-lipid interactions.     

Differential lipid packing abilities and dynamics in giant unilamellar vesicles composed of short-chain saturated glycerol-phospholipids, sphingomyelin and cholesterol

The ability of membrane components to arrange themselves heterogeneously within the bilayer induces the formation of microdomains. Much work has been devoted to mimicking domain-assembly in artificial bilayers and characterizing their physico-chemical properties. Ternary lipid mixtures composed of unsaturated phospholipids, sphingomyelin and cholesterol give rise to large, round domains. Here, we replaced the unsaturated phospholipid in the ternary mixture with sphingomyelin and cholesterol by saturated glycero-phospholipids of different chain length and characterized the critical role of cholesterol in sorting these lipids by confocal imaging and fluorescence correlation spectroscopy (FCS). More cholesterol is needed to obtain phase segregation in ternary mixtures, in which the unsaturated phospholipid is replaced by a saturated one. Finally, lipid dynamics in distinct phases is very low and astonishingly similar, thereby suggesting the poor ability of cholesterol in sorting short-chain saturated glycero-phospholipids and sphingomyelin.     

Effects of fluorescent probe NBD-PC on the structure, dynamics and phase transition of DPPC. A molecular dynamics and differential scanning calorimetry study

We present a combined theoretical (molecular dynamics, MD) and experimental (differential scanning calorimetry, DSC) study of the effect of 7-nitrobenz-2-oxa-1,3-diazol-4-yl (NBD) acyl chain-labeled fluorescent phospholipid analogs (C6-NBD-PC and C12-NBD-PC) on 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) bilayers. DSC measurements reveal that < 1 mol% of NBD-PC causes elimination of the pre-transition and a large loss of cooperativity of the main transition of DPPC. Labeling with C6-NBD-PC or C12-NBD-PC shifts the main transition temperature to lower or higher values, respectively. Following our recent report on the location and dynamics of these probes (BBA 1768 (2007) 467-478) in fluid phase DPPC, we present a detailed analysis of 100-ns MD simulations of systems containing either C6-NBD-PC or C12-NBD-PC, focused on their influence on several properties of the host bilayer. Whereas most monitored parameters are not severely affected for 1.6 mol% of probe, for the higher concentration studied (6.2 mol%) important differences are evident. In agreement with published reports, we observed that the average area per phospholipid molecule increases, whereas DPPC acyl chain order parameters decrease. Moreover, we predict that incorporation of NBD-PC should increase the electrostatic potential across the bilayer and, especially for C12-NBD-PC, slow lateral diffusion of DPPC molecules and rotational mobility of DPPC acyl chains.