Molecular dynamics simulation of a phosphatidylglycerol membrane

Although molecular dynamics simulations are an important tool for studying membrane systems, relatively few simulations have used anionic lipids. This paper reports the first simulation of a pure phosphatidylglycerol (PG) bilayer. The properties of this equilibrated palmitoyloleoylphosphatidylglycerol membrane agree with experimental observations of PG membranes and with previous simulations of monolayers and mixed bilayers containing PG lipids. These simulations also provide interesting insights into hydrogen bonding interactions in PG membranes. This equilibrated membrane will be a useful starting point for simulations of membrane proteins interacting with PG lipids.     

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 simulation of cation-phospholipid clustering in phospholipid bilayers: possible role in stalk formation during membrane fusion

In this study, we performed all-atom long-timescale molecular dynamics simulations of phospholipid bilayers incorporating three different proportions of negatively charged lipids in the presence of K(+), Mg(2+), and Ca(2+) ions to systemically determine how membrane properties are affected by cations and lipid compositions. Our simulations revealed that the binding affinity of Ca(2+) ions with lipids is significantly stronger than that of K(+) and Mg(2+) ions, regardless of the composition of the lipid bilayer. The binding of Ca(2+) ions to the lipids resulted in bilayers having smaller lateral areas, greater thicknesses, greater order, and slower rotation of their lipid head groups, relative to those of corresponding K(+)- and Mg(2+)-containing systems. The Ca(2+) ions bind preferentially to the phosphate groups of the lipids. The complexes formed between the cations and the lipids further assembled to form various multiple-cation-centered clusters in the presence of anionic lipids and at higher ionic strength-most notably for Ca(2+). The formation of cation-lipid complexes and clusters dehydrated and neutralized the anionic lipids, creating a more-hydrophobic environment suitable for membrane aggregation. We propose that the formation of Ca(2+)-phospholipid clusters across apposed lipid bilayers can work as a "cation glue" to adhere apposed membranes together, providing an adequate configuration for stalk formation during membrane fusion.     

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.     

Molecular mechanisms of membrane fusion: Steps during phospholipid and exocytotic membrane fusion

Exocytosis is considered as four separate steps: adhesion, fusion/pore formation, pore widening, and content discharge. Experiments on both synthetic and natural membranes are presented to show each of these steps. Major differences are seen in the two fusing systems. These differences are discussed in terms of molecular mechanisms of fusion.     

Molecular dynamics simulation of dipalmitoylphosphatidylcholine membrane with cholesterol sulfate

Using the molecular dynamics simulation technique, we studied the changes occurring in a dipalmitoylphosphatidylcholine (DPPC):cholesterol (CH) membrane at 50 mol% sterol when cholesterol is replaced with cholesterol sulfate (CS). Our simulations were performed at constant pressure and temperature on a nanosecond time scale. We found that 1) the area per DPPC:CS heterodimer is greater than the area of the DPPC:CH heterodimer; 2) CS increases ordering of DPPC acyl chains, but to a lesser extent than CH; 3) the number of hydrogen bonds between DPPC and water is decreased in a CS-containing membrane, but CS forms more water hydrogen bonds than CH; and 4) the membrane dipole potential reverses its sign for a DPPC-CS membrane compared to a DPPC-CH bilayer. We also studied the changes occurring in lipid headgroup conformations and determined the location of CS molecules in the membrane. Our results are in good agreement with the data available from experiments.     

Molecular dynamics simulation of spherical-to -threadlike micelle transition in a cationic surfactant solution

The effect of sodium salicylate (NaSal) on the spherical-to-threadlike micelle shape transition in 3-hexadecyloxy-2-hydroxy-propyl trimethyl ammonium bromide (R₁₆HTAB) solution was studied using molecular dynamics simulation. The simulations were started from a preassembled infinitely long threadlike micelle of R₁₆HTAB. By analyzing the aggregation morphologies and structural details, we find that the preassembled threadlike micelle in the absence of NaSal was unstable and assembled into a spherical micelle. While in the presence of NaSal, the threadlike micelle exhibited fluctuations but remained the threadlike shape during the long simulation run. The Sal^? ions were found to penetrate inside the micelle, which promoted the junction between the surfactant and salicylate counterion. The aromatic Sal^? ions located in the surfactant headgroup region with their phenyl groups pointing toward the interior core region of the micelle. From another simulation started with two individual spherical micelles, we found that the Sal^? ions can link the two spherical micelles into a long threadlike micelle, in accordance with a mode proposed by experimental studies. Our studies showed that the H-bonds and electrostatic interactions between the Sal^? ions and the surfactants played an important role in micellar growth and stabilising the threadlike micelle.     

Effect of membrane phase transition on long-time calcium-induced fusion of phosphatidylserine vesicles

Dynamic light scattering has been used to study the temperature dependence of the extent of long-time calcium-induced fusion of sonicated vesicles composed of various natural and synthetic phosphatidylserine with different acyl chains. The vesicles of each composition are found to exhibit a peak temperature in the vicinity of which the extent of fusion shows a distinct maximum. The fusion peak temperature increases as the bilayer gel-to-liquid-crystal phase transition temperature increases. The results suggest a role played by membrane fluidity in determining fusion efficiency.     

Molecular dynamics simulation of proton transport near the surface of a phospholipid membrane

The structural and dynamical properties of a hydrated proton near the surface of DMPC membrane were studied using a molecular dynamics simulation. The proton transport between water molecules was modeled using the second generation multistate empirical valence bond model. The proton diffusion was found to be inhibited at the membrane surface. The potential of mean force for the proton adsorption to the membrane surface and its release back into the bulk water was also determined, yielding a small barrier in each direction. An efficient algorithm for Ewald summation calculations for the multistate empirical valence bond model is also introduced.     

Molecular dynamics simulation of the evolution of hydrophobic defects in one monolayer of a phosphatidylcholine bilayer: relevance for membrane fusion mechanisms

The spontaneous formation of the phospholipid bilayer underlies the permeability barrier function of the biological membrane. Tears or defects that expose water to the acyl chains are spontaneously healed by lipid lateral diffusion. However, mechanical barriers, e.g., protein aggregates held in place, could sustain hydrophobic defects. Such defects have been postulated to occur in processes such as membrane fusion. This gives rise to a new question in bilayer structure: What do the lipids do in the absence of lipid lateral diffusion to minimize the free energy of a hydrophobic defect? As a first step to understand this rather fundamental question about bilayer structure, we performed molecular dynamic simulations of up to 10 ns of a planar bilayer from which lipids have been deleted randomly from one monolayer. In one set of simulations, approximately one-half of the lipids in the defect monolayer were restrained to form a mechanical barrier. In the second set, lipids were free to diffuse around. The question was simply whether the defects caused by removing a lipid would aggregate together, forming a large hydrophobic cavity, or whether the membrane would adjust in another way. When there are no mechanical barriers, the lipids in the defect monolayer simply spread out and thin with little effect on the other intact monolayer. In the presence of a mechanical barrier, the behavior of the lipids depends on the size of the defect. When 3 of 64 lipids are removed, the remaining lipids adjust the lower one-half of their chains, but the headgroup structure changes little and the intact monolayer is unaffected. When 6 to 12 lipids are removed, the defect monolayer thins, lipid disorder increases, and lipids from the intact monolayer move toward the defect monolayer. Whereas this is a highly simplified model of a fusion site, this engagement of the intact monolayer into the fusion defect is strikingly consistent with recent results for influenza hemagglutinin mediated fusion.     

Molecular dynamics simulation of Bombolitin II in the dipalmitoylphosphatidylcholine membrane bilayer

The orientation behavior of Bombolitin II (BLT2) in the dipalmitoylphosphatidylcholine membrane bilayer was investigated by using molecular-dynamics simulation. During the 20-ns simulation, the BLT2 began to tilt and finally reached the angle of 51° from the membrane-normal. The structure of the peptide formed the amphipathic α-helical structure during the entire simulation time. The peptide tilts with its hydrophobic side faced to the hydrophobic core of the bilayer. We analyzed the mechanism of the tilting behavior of the peptide associated with the membrane in detail. The analysis showed that the hydrogen-bond interaction and the electrostatic interaction were found to exist between Lys¹² and a lipid molecule. These interactions are considered to work as an important factor in tilting the peptide to the membrane-normal.     

Spatiotemporal dynamics of membrane remodeling and fusion proteins during endocytic transport

Organelles of the endolysosomal system undergo multiple fission and fusion events to combine sorting of selected proteins to the vacuole with endosomal recycling. This sorting requires a consecutive remodeling of the organelle surface in the course of endosomal maturation. Here we dissect the remodeling and fusion machinery on endosomes during the process of endocytosis. We traced selected GFP-tagged endosomal proteins relative to exogenously added fluorescently labeled α-factor on its way from the plasma membrane to the vacuole. Our data reveal that the machinery of endosomal fusion and ESCRT proteins has similar temporal localization on endosomes, whereas they precede the retromer cargo recognition complex. Neither deletion of retromer nor the fusion machinery with the vacuole affects this maturation process, although the kinetics seems to be delayed due to ESCRT deletion. Of importance, in strains lacking the active Rab7-like Ypt7 or the vacuolar SNARE fusion machinery, α-factor still proceeds to late endosomes with the same kinetics. This indicates that endosomal maturation is mainly controlled by the early endosomal fusion and remodeling machinery but not the downstream Rab Ypt7 or the SNARE machinery. Our data thus provide important further understanding of endosomal biogenesis in the context of cargo sorting.     

Molecular dynamics simulation of dipalmitoylphosphatidylcholine modified with a MTSL nitroxide spin label in a lipid membrane

We investigate the interaction between dipalmitoylphosphatidylcholine (DPPC) and a nitroxide spin label in order to understand its influences on lipid structure and dynamics using molecular dynamics simulations. The system was modified by covalently attaching nitroxide spin labels to the headgroups of two DPPC molecules. (S-(2,2,5,5-tetramethyl-2,5-dihydro-1H-pyrrol-3-yl)methyl methanesulfonothioate) (MTSL) was used as the spin label. The label position and dynamics were analyzed as was the impact of the modified DPPC on the structure of the surrounding lipids. The modified DPPC molecules locate closer to the center of the membrane than unmodified DPPC molecules. The rotation of the spin label is unrestricted, but there are favored orientations. MTSL depresses the deuterium order parameters of the carbon atoms close to the headgroup in surrounding DPPC molecules. The spin label has no impact on order parameters of carbon atoms at the end of the lipid tails. The lateral diffusion constant of the modified DPPC is indistinguishable from unmodified DPPC molecules. These novel computational results suggest an experimental validation.     

Molecular mechanisms of flavivirus membrane fusion

Flaviviruses comprise a number of important human pathogens including yellow fever, dengue, West Nile, Japanese encephalitis and tick-borne encephalitis viruses. They are small enveloped viruses that enter cells by receptor-mediated endocytosis and release their nucleocapsid into the cytoplasm by fusing their membrane with the endosomal membrane. The fusion event is triggered by the acidic pH in the endosome and is mediated by the major envelope protein E. Based on the atomic structures of the pre- and post-fusion conformations of E, a fusion model has been proposed that includes several steps leading from the metastable assembly of E at the virion surface to membrane merger and fusion pore formation trough conversion of E into a stable trimeric post-fusion conformation. Using recombinant subviral particles of tick-borne encephalitis virus as a model, we have defined individual steps of the molecular processes underlying the flavivirus fusion mechanisms. This includes the identification of a conserved histidine as being part of the pH sensor in the fusion protein that responds to the acidic pH and thus initiates the structural transitions driving fusion.     

Molecular mechanism of mitochondrial membrane fusion

Mitochondrial fusion requires coordinated fusion of the outer and inner membranes. This process leads to exchange of contents, controls the shape of mitochondria, and is important for mitochondrial function. Two types of mitochondrial GTPases are essential for mitochondrial fusion. On the outer membrane, the fuzzy onions/mitofusin proteins form complexes in trans that mediate homotypic physical interactions between adjacent mitochondria and are likely directly involved in outer membrane fusion. Associated with the inner membrane, the OPA1 dynamin-family GTPase maintains membrane structure and is a good candidate for mediating inner membrane fusion. In yeast, Ugo1p binds to both of these GTPases to form a fusion complex, although a related protein has yet to be found in mammals. An understanding of the molecular mechanism of fusion may have implications for Charcot-Marie-Tooth subtype 2A and autosomal dominant optic atrophy, neurodegenerative diseases caused by mutations in Mfn2 and OPA1.     

Molecular mechanisms of calcium-induced membrane fusion

We have reviewed studies on calcium-induced fusion of lipid bilayer membranes and the role of synexin and other calcium-binding proteins (annexins) in membrane fusion. We have also discussed the roles of other cations, lipid phase transitions, long chain fatty acids and other fusogenic molecules. Finally, we have presented a simple molecular model for the mechanism of lipid membrane fusion, consistent with the experimental evidence and incorporating various elements proposed previously.     

A novel phase of compressed bilayers that models the prestalk transition state of membrane fusion

The force model of protein-mediated membrane fusion hypothesizes that fusion is driven by mechanical forces exerted on the membranes, but many details are unknown. Here, we investigated by x-ray diffraction the consequence of applying compressive force on a stack of membranes against the hydration barrier. We found that as the osmotic pressure increased, the lamellar phase transformed first to a new phase of tetragonal lattice (T-phase) over a narrow range of relative humidity, and then to a phase of rhombohedral lattice. The unit cell structure changed from parallel bilayers to a bent configuration with a point contact between adjacent bilayers and then to the stalk hemifusion configuration. The T-phase is discussed as a possible transition state in the membrane merging pathway of fusion. We estimate the work required to form the T-phase and the subsequent hemifusion-stalk-resembling R-phase. The work for the formation of a stalk is compatible with the energy estimated to be released by several SNARE complexes.     

Molecular dynamics simulation of a high-affinity antibody-protein complex

One nanosecond molecular dynamic (MD) simulation of anti-hen egg white lysozyme (HEL) antibody HyHEL63 (HH63) complexed with HEL reveals rigid and flexible regions of the HH63 binding site. Fifty conformations, extracted from the MD trajectory at regular time intervals were superimposed on HH63-HEL X-ray crystal structure, and the root mean squared deviations (RMSDs) and deviations in Calpha atom positions between the X-ray structure and the MD conformer were measured. Residue positions showing the large deviations in both light chain and heavy chain of the antibody were same in all the MD conformers. The residue positions showing smallest deviations were same for all the conformers in the case of light chain, whereas relatively variable in the heavy chain. Positions of large and small deviations fell in the complementarity determining regions (CDRs), for both heavy and light chains. The larger deviations were in CDR-2 of light and CDR-1 of heavy chain. Smaller deviations were in CDR-3 of light and CDR-2 and CDR-3 of heavy chains. The large and small deviating regions highlight flexible and rigid regions of HH63 binding site and suggest a mosaic binding mechanism, including both "induced fit" and preconfigured "lock-and-key" type of binding. Combined "induced fit" and "lock-and-key" binding would be a better definition for the formation of large complexes, which bury larger surface area on binding, as in the case of antibody-HEL complex. We further show that flexible regions, comprising mostly charged and polar residues, form intermolecular interactions with HEL, whereas rigid regions do not. Electrostatic complementarity between HH63 and HEL also imply optimized binding affinity. Flexible and rigid regions of a high-affinity antibody are selected during the affinity maturation of the antibody and have specific functional significance. The functional importance of local inherently flexible regions is to establish intermolecular contacts or they play a key role in molecular recognition, whereas local rigid regions provide the structural framework.