Molecular dynamics simulations of cardiolipin bilayers

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

Microscopic structure of phospholipid bilayers: comparison between molecular dynamics simulations and wide-angle X-ray spectra

We present results of molecular dynamics simulations of fully hydrated dipalmitoylphosphatidylcholine and dimyristoylphosphatidylcholine bilayers in the disordered liquid crystalline phase (Lalpha) and compare them to wide-angle X-ray scattering experiments. Though we find a generally good agreement between the simulated and experimental spectra, there are some deviations whose origin has been investigated by a reparametrization of the aliphatic chains' force field. A detailed analysis of the various contribution to the X-ray spectra shows that a non-negligible contribution to the total scattered intensity comes from the headgroups and the head-tail cross correlation.     

Molecular dynamics simulations of rhodopsin in different one-component lipid bilayers

Four 20 ns molecular dynamic simulations of rhodopsin embedded in different one-component lipid bilayers have been carried out to ascertain the importance of membrane lipids on the protein structure. Specifically, dimyristoyl phosphatidylcholine (DMPC), dipalmitoyl phosphatidylcholine (DPPC), palmitoyl oleoyl phosphatidylcholine (POPC), and palmitoyl linoleyl phosphatidylcholine (PLPC) lipid bilayers have been considered for the present work. The results reported here provide information on the hydrophobic matching between the protein and the bilayer and about the differential effects of the protein on the thickness of the different membranes. Furthermore, a careful analysis of the individual protein-lipid interactions permits the identification of residues that exhibit permanent interactions with atoms of the lipid environment that may putatively act as hooks of the protein to the membrane. The analysis of the trajectories also provides information about the effect of the bilayer on the protein structure, including secondary structural elements, salt bridges, and rigid-body motions.     

Molecular dynamics simulation of the formation, structure, and dynamics of small phospholipid vesicles

Here, we use coarse grained molecular dynamics (MD) simulations to study the spontaneous aggregation of dipalmitoylphosphatidylcholine (DPPC) lipids into small unilamellar vesicles. We show that the aggregation process occurs on a nanosecond time scale, with bicelles and cuplike vesicles formed at intermediate stages. Formation of hemifused vesicles is also observed at higher lipid concentration. With either 25% dipalmitoylphosphatidylethanolamine (DPPE) or lysoPC mixed into the system, the final stages of the aggregation process occur significantly faster. The structure of the spontaneously formed vesicles is analyzed in detail. Microsecond simulations of isolated vesicles reveal significant differences in the packing of the lipids between the inner and outer monolayers, and between PC, PE, and lysoPC. Due to the small size of the vesicles they remain almost perfectly spherical, undergoing very limited shape fluctuations or bilayer undulations. The lipid lateral diffusion rate is found to be faster in the outer than in the inner monolayer. The water permeability coefficient of the pure DPPC vesicles is of the order of 10(-)(3) cm s(-)(1), in agreement with experimental measurements.     

Molecular dynamics simulations of neutral chlorpromazine in zwitterionic phospholipid monolayers

Molecular dynamics simulations have been performed to investigate the interactions between chlorpromazine (CPZ), a neuroleptic drug used in the treatment of psychiatric disorders, and dipalmitoylphosphatidylcholine (DPPC), a zwitterionic phospholipid, in Langmuir monolayers. The results from simulations carried out at different monolayer surface densities were able to capture important features of the CPZ-lipid interaction. We find that neutral (unprotonated) CPZ is preferentially located in the lipid tail region of the phospholipids, in little contact with the aqueous phase, and that the orientation of its rigid ring structure and tail conformation vary with lipid surface density. CPZ is found to promote ordering of the lipid tails for all surface densities because of a reduction in the effective surface area per lipid upon addition of the drug. Similar effects have been observed in previous studies of cholesterol in DPPC monolayers, in which lipid tails were seen to order around the solute. This feature, however, is quite distinct from what we observe for the most dense monolayer considered here (area per lipid of 50 A(2)), for which we find that CPZ promotes a local distortion of the lipid tails in its immediate vicinity and a concomitant ordering of lipid tails located further away from the solute. This view is further supported by the results obtained for an approximated nonlinear vibrational sum frequency generation susceptibility, which showed greater tail disorder close to CPZ.     

Molecular dynamics simulations of bilayers containing mixtures of sphingomyelin with cholesterol and phosphatidylcholine with cholesterol

We performed six molecular dynamics simulations: three on hydrated bilayers containing pure phospholipids and three on hydrated bilayers containing mixtures of these phospholipids with cholesterol. The phospholipids in our simulations were SSM (sphingomyelin containing a saturated 18:0 acyl chain), OSM (sphingomyelin with an unsaturated 18:1 acyl chain), and POPC (palmitoyloleoylphosphatidylcholine containing one saturated and one unsaturated chain). Data from our simulations were used to study systematically the effect of cholesterol on phospholipids that differed in their headgroup and tail composition. In addition to the structural analysis, we performed an energetic analysis and observed that energies of interaction between cholesterol and neighboring SM molecules are similar to the energies of interaction between cholesterol and POPC. We also observed that the interaction energy between cholesterol and neighboring lipids cannot be used for the determination of which lipids are involved in the creation of a complex.     

Improved dissipative particle dynamics simulations of lipid bilayers

The authors introduce a new parameterization for the dissipative particle dynamics simulations of lipid bilayers. In this parameterization, the conservative pairwise forces between beads of the same type in two different hydrophobic chains are chosen to be less repulsive than the water-water interaction, but the intrachain bead interactions are the same as the water-water interaction. For a certain range of parameters, the new bilayer can only be stretched up to 30% before it ruptures. Membrane tension, density profiles, and the in-plane lipid diffusion coefficient of the new bilayer are discussed in detail. They find two kinds of finite size effects that influence the membrane tension: lateral finite size effects, for which larger membranes rupture at a smaller stretch, and transverse finite size effects, for which tensionless bilayers are more compact in larger systems. These finite size effects become rather small when the simulation box is sufficiently large.     

Comparative molecular dynamics study of ether- and ester-linked phospholipid bilayers

The lipid membranes found in archaea have high bilayer stability and low permeability. The molecular structure of their constituent lipids is characterized by ether-linked, branched hydrophobic chains, whereas the conventional lipids obtained from eukaryotic or eubacterial sources have ester linked straight chains. In order to elucidate the influence of the ether linkage, instead of an ester one, on the physical properties of the lipid bilayers, we have carried out comparative 10 ns molecular dynamics simulations of diphytanyl phosphatidylcholine (ether-DPhPC) and diphytanoyl phosphatidylcholine (ester-DPhPC) bilayers in water, respectively. We analyze bilayer structures, hydration of the lipids, membrane dipole potentials, and free energy profiles of water and oxygen across the bilayers. We observe that the membrane dipole potential for the ether-DPhPC bilayer, which arises mainly from the ether linkage, is about half of that of the ester-DPhPC. The calculated free energy barrier for a water molecule in the ether-DPhPC bilayer system is slightly higher than that in the ester-DPhPC counterpart, which is in accord with experimental data.     

Molecular dynamics simulations of PAMAM dendrimer-induced pore formation in DPPC bilayers with a coarse-grained model

We have performed 0.5-micros-long molecular dynamics (MD) simulations of 0%, 50%, and 100% acetylated third- (G3) and fifth-generation (G5) polyamidoamine (PAMAM) dendrimers in dipalmitoylphosphatidylcholine (DPPC) bilayers with explicit water using the coarse-grained (CG) model developed by Marrink et al. (J.Phys. Chem. B 2004, 108, 750-760), but with long-range electrostatic interactions included. Radii of gyration of the CG G5 dendrimers are 1.99-2.32 nm, close to those measured in the experiments by Prosa et al. (J. Polym. Sci. 1997, 35, 2913-2924) and atomistic simulations by Lee et al. (J. Phys. Chem. B 2006, 110, 4014-4019). Starting with the dendrimer initially positioned near the bilayer, we find that positively charged un-acetylated G3 and 50%-acetylated and un-acetylated G5 dendrimers insert themselves into the bilayer, and only un-acetylated G5 dendrimer induces hole formation at 310 K, but not at 277 K, which agrees qualitatively with experimental observations of Hong et al. (Bioconj. Chem. 2004, 15, 774-782) and Mecke et al. (Langmuir 2005, 21, 10348-10354). At higher salt concentration (approximately 500 mM NaCl), un-acetylated G5 dendrimer does not insert into the bilayer. The results suggest that with inclusion of long-range electrostatic interactions into coarse-grained models, realistic MD simulation of membrane-disrupting effects of nanoparticles at the microsecond time scale is now possible.     

Structure of docosahexaenoic acid-containing phospholipid bilayers as studied by (2)H NMR and molecular dynamics simulations.

Polyunsaturated phospholipids are known to be important with regard to the biological functions of essential fatty acids, for example, involving neural tissues such as the brain and retina. Here we have employed two complementary structural methods for the study of polyunsaturated bilayer lipids, viz. deuterium ((2)H) NMR spectroscopy and molecular dynamics (MD) computer simulations. Our research constitutes one of the first applications of all-atom MD simulations to polyunsaturated lipids containing docosahexaenoic acid (DHA; 22:6 cis-Delta(4,7,10,13,16,19)). Structural features of the highly unsaturated, mixed-chain phospholipid, 1-palmitoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine (PDPC), have been studied in the liquid-crystalline (L(alpha)) state and compared to the less unsaturated homolog, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC). The (2)H NMR spectra of polyunsaturated bilayers are dramatically different from those of less unsaturated phospholipid bilayers. We show how use of MD simulations can aid in interpreting the complex (2)H NMR spectra of polyunsaturated bilayers, in conjunction with electron density profiles determined from small-angle X-ray diffraction studies. This work clearly demonstrates preferred helical and angle-iron conformations of the polyunsaturated chains in liquid-crystalline bilayers, which favor chain extension while maintaining bilayer flexibility. The presence of relatively long, extended fatty acyl chains may be important for solvating the hydrophobic surfaces of integral membrane proteins, such as rhodopsin. In addition, the polyallylic DHA chains have a tendency to adopt back-bended (hairpin-like) structures, which increase the interfacial area per lipid. Finally, the material properties have been analyzed in terms of the response of the bilayer to mechanical stress. Simulated bilayers of phospholipids containing docosahexaenoic acid were less sensitive to the applied surface tension than were saturated phospholipids, possibly implying a decrease in membrane elasticity (area elastic modulus, bending rigidity). The above features distinguish DHA-containing lipids from saturated or monounsaturated lipids and may be important for their biological modes of action.     

Interface water dynamics and porating electric fields for phospholipid bilayers

Lipid bilayers, normally a barrier to charged species and large molecules, are permeabilized by electric fields, a phenomenon exploited by cell biologists and geneticists for porating and transfecting cells and tissues. Recent molecular simulation studies have advanced our understanding of electroporation, but the relative contributions of atomically local details (interface water and headgroup dipole and counterion configurations) and medium- and long-range electrostatic gradients and changes in membrane structural shifts to the initiating conditions and mechanisms of pore formation remain unclear. Molecular dynamics simulations of electroporation in several lipid systems presented here reveal the effects of lipid hydrocarbon tail length and composition on the magnitude of the field required for poration and on the location of the initial sites of field-driven water intrusion into the bilayer. Minimum porating external fields of 260 mV nm(-1), 280 mV nm(-1), 320 mV nm(-1), and 380 mV nm(-1) were found for 1,2-dilauroyl-sn-glycero-3-phosphatidylcholine (DLPC), 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine (DPPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC), and 1,2-dioleoyl- sn-glycero-3-phosphatidylcholine (DOPC) bilayers, respectively, and correlated most strongly with the bilayer thickness. These phospholipid systems share several common features including a wide, dynamic distribution of the headgroup dipole angle with the bilayer normal ranging from 0 to 155 degrees that is only slightly shifted in a porating electric field, and similar electric field-induced shifts in water dipole orientation, although the mean water dipole moment profile at the aqueous-membrane interface is more sensitive to the electric field for DOPC than for the other phospholipids. The location of pore initiation, at the anode- or cathode-facing leaflet, varies with the composition of the bilayer and correlates with a change in the polarity of the localized electric field at the interface.     

Structure and dynamics of water at the interface with phospholipid bilayers

We have performed two molecular-dynamics simulations to study the structural and dynamical properties of water at the interface with phospholipid bilayers. In one of the simulations the bilayer contained neutral phospholipid molecules, dioleoylphosphatidylcholine (DOPC); in the second simulation the bilayer contained charged lipid molecules, dioleoylphosphatidylserine (DOPS). From the density profile of water we observe that water next to the DOPS bilayer is more perturbed as compared to water near the DOPC bilayer. Using an energetic criterion for the determination of hydrogen bonding we find that water molecules create strong hydrogen bonds with the headgroups of the phospholipid molecules. Due to the presence of these bonds and also due to the confinement of water, the translational and orientational dynamics of water at the interface are slowed down. The degree of slowing down of the dynamics depends upon the location of water molecules near a lipid headgroup.     

Isomerization dynamics of photochromic spiropyran molecular switches in phospholipid bilayers

Transient absorption spectroscopy was used to investigate the dynamics of the photochromic indolinobenzospiropyran reaction in toluene solution and in phosphatidylcholine bilayers (1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), and 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC)). After excitation with UV light, colorless (R/S)-2-(3',3'-dimethyl-6-nitro-3'H-spiro[chromene-2,2'-indol]-1'-yl)ethanol derivatives are converted to colored merocyanine products in high yield; Phi = 0.45 in DMPC liposomes. We find that the reaction occurs in the bilayer aliphatic region in the gel (P(beta)(')) and liquid (L(alpha)) phases. The Arrhenius activation energy for the isomerization in DMPC bilayers was approximately 3.5 times larger in the liquid phase (L(alpha), E(a) = 26.0 +/- 1.0 kJ mol(-1)) than that in the gel phase (P(beta)('), E(a) = 7.3 +/- 1.6 kJ mol(-1)). Analysis of the isomerization rate constant temperature dependence allows an estimation of the bilayer viscosity and free volume properties in the L(alpha) phase.     

Molecular dynamics simulations of extensional,sheet and unidirectional flow

Non-equilibrium molecular dynamics, NEMD, has been used to determine the technologically important tensile viscosities. Simulations of elongational, sheet and unidirectional viscoelastic responses via the “differences-in-trajectories” technique have been made on Lennard-Jones liquid at ?* = 0.8442 and T* = 0.72, and ?* = 1.01304 and T* = 1.26. At finite strain rates (·γ* ≥ 0.1) shear thinning behavior is evident, in particular for unidirectional (“poker chip”) extensional flow. At γ ≥ 0.4 steady state unidirectional viscosities were not achieved before severe material deterioration. These trends are the result of extremely complex structural reorganisation well within the viscoelastic regime, which involved coupled component relaxation moduli along distortion direction.     

Area compressibility and buckling of amphiphilic bilayers in molecular dynamics simulations

The elastic modulus or area compressibility of a membrane is routinely calculated in molecular dynamics simulations as the proportionality constant relating surface tension and projected surface area. Recent studies, however, have revealed a marked system size dependence of these moduli, which we attribute to the neglect of thermal undulations in the area calculation. We discuss several methods, based on the Helfrich model and on numerical triangulation, to remedy this situation, and find a satisfying agreement between them. The Helfrich model also quantitatively describes a buckling transition observed for compressed bilayers.     

Origin of short-range repulsion between hydrated phospholipid bilayers: a computer simulation study

The grand canonical Monte Carlo technique is used to simulate the pressure-distance dependence for supported dilauroylphosphatidylethanolamine (DLPE) membranes. The intra- and intermolecular interactions in the system are described with a combination of an AMBER-based force field for DLPE and a TIP4P model for water. To improve the balance between the pair interactions of like and unlike molecules, the water-lipid interaction potentials are scaled to reproduce the hydration level and intermembrane separation at full hydration. It is found that the short-range water-mediated repulsion originates from the hydration component of the intermembrane pressure, whereas the direct interaction between the membranes remains attractive throughout the pressure range studied (0-5 kbar).     

Molecular dynamics simulations of halide glasses

Summary Halide glasses have been extensively studied in recent years because of their potential application as infrared transmitting fibre optic materials. They are believed to be more ionic than glasses based on silica and should therefore be more amenable to molecular dynamics simulation using simple two-body potentials. Here the main features of structural models derived using such techniques are described and compared with available structural data. Possible future applications of this approach are outlined.     

Molecular dynamics simulations of peptide-surface interactions

Proteins, which are bioactive molecules, adsorb on implants placed in the body through complex and poorly understood mechanisms and directly influence biocompatibility. Molecular dynamics modeling using empirical force fields provides one of the most direct methods of theoretically analyzing the behavior of complex molecular systems and is well-suited for the simulation of protein adsorption behavior. To accurately simulate protein adsorption behavior, a force field must correctly represent the thermodynamic driving forces that govern peptide residue-surface interactions. However, since existing force fields were developed without specific consideration of protein-surface interactions, they may not accurately represent this type of molecular behavior. To address this concern, we developed a host-guest peptide adsorption model in the form of a G(4)-X-G(4) peptide (G is glycine, X is a variable residue) to enable determination of the contributions to adsorption free energy of different X residues when adsorbed to functionalized Au-alkanethiol self-assembled monolayers (SAMs). We have previously reported experimental results using surface plasmon resonance (SPR) spectroscopy to measure the free energy of peptide adsorption for this peptide model with X = G and K (lysine) on OH and COOH functionalized SAMs. The objectives of the present research were the development and assessment of methods to calculate adsorption free energy using molecular dynamics simulations with the GROMACS force field for these same peptide adsorption systems, with an oligoethylene oxide (OEG) functionalized SAM surface also being considered. By comparing simulation results to the experimental results, the accuracy of the selected force field to represent the behavior of these molecular systems can be evaluated. From our simulations, the G(4)-G-G(4) and G(4)-K-G(4) peptides showed minimal to no adsorption to the OH SAM surfaces and the G(4)-K-G(4) showed strong adsorption to the COOH SAM surface, which is in agreement with our SPR experiments. Contrary to our experimental results, however, the simulations predicted a relatively strong adsorption of G(4)-G-G(4) peptide to the COOH SAM surface. In addition, both peptides were unexpectedly predicted to adsorb to the OEG surface. These findings demonstrate the need for GROMACS force field parameters to be rebalanced for the simulation of peptide adsorption behavior on SAM surfaces. The developed methods provide a direct means of assessing, modifying, and validating force field performance for the simulation of peptide and protein adsorption to surfaces, without which little confidence can be placed in the simulation results that are generated with these types of systems.     

Interleaflet interaction and asymmetry in phase separated lipid bilayers: molecular dynamics simulations

In order to investigate experimentally inaccessible, molecular-level detail regarding interleaflet interaction in membranes, we have run an extensive series of coarse-grained molecular dynamics simulations of phase separated lipid bilayers. The simulations are motivated by differences in lipid and cholesterol composition in the inner and outer leaflets of biological membranes. Over the past several years, this phenomenon has inspired a series of experiments in model membrane systems which have explored the effects of lipid compositional asymmetry in the two leaflets. The simulations are directed at understanding one potential consequence of compositional asymmetry, that being regions of bilayers where liquid-ordered (L(o)) domains in one leaflet are opposite liquid-disordered (L(d)) domains in the other leaflet (phase asymmetry). The simulated bilayers are of two sorts: 1) Compositionally symmetric leaflets where each of the two leaflets contains an identical, phase separated (L(o)/L(d)) mixture of cholesterol, saturated and unsaturated phospholipid; and 2) Compositionally asymmetric leaflets, where one leaflet contains a phase separated (L(o)/L(d)) mixture while the other contains only unsaturated lipid, which on its own would be in the L(d) phase. In addition, we have run simulations where the lengths of the saturated lipid chains as well as the mole ratios of the three lipid components are varied. Collectively, we report on three types of interleaflet coupling within a bilayer. First, we show the effects of compositional asymmetry on acyl chain tilt and order, lipid rotational dynamics, and lateral diffusion in regions of leaflets that are opposite L(o) domains. Second, we show substantial effects of compositional asymmetry on local bilayer curvature, with the conclusion that phase separated leaflets resist curvature, while inducing large degrees of curvature in an opposing L(d) leaflet. Finally, in compositionally symmetric, phase separated bilayers, we find phase asymmetry (domain antiregistration) between the two leaflets occurs as a consequence of mismatched acyl chain-lengths in the saturated and unsaturated lipids.