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.     

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.     

Molecular dynamics simulation of a dipalmitoylphosphatidylcholine bilayer with NaCl

Molecular dynamics simulations are performed on two hydrated dipalmitoylphosphatidylcholine bilayer systems: one with pure water and one with added NaCl. Due to the rugged nature of the membrane/electrolyte interface, ion binding to the membrane surface is characterized by the loss of ion hydration. Using this structural characterization, binding of Na(+) and Cl(-) ions to the membrane is observed, although the binding of Cl(-) is seen to be slightly weaker than that of Na(+). Dehydration is seen to occur to a different extent for each type of ion. In addition, the excess binding of Na(+) gives rise to a net positive surface charge density just outside the bilayer. The positive density produces a positive electrostatic potential in this region, whereas the system without salt shows an electrostatic potential of zero.     

Molecular dynamics simulation of DPPC bilayer in DMSO

We performed molecular dynamics simulations on dipalmitoylphosphatidylcholine (DPPC)/dimethylsulfoxide (DMSO) system that has the same lipid:solvent weight ratio as in our previous simulation done on DPPC/water. We did not observe a large change in the size of DPPC membrane when the solvent was changed from water to DMSO. Also, we did not observe that a large number of DMSO molecules is permeating into the membrane, as it was suggested to explain the observed change in the bilayer repeat period. We found that the surface potential reverses its sign when water is replaced by DMSO. Based on the results from our simulations, we propose that the repulsion force acting between membranes is reduced when DMSO is added to solvent water and therefore membrane surfaces approach closer to each other and the extra solvent is removed into excess solution.     

Molecular dynamics simulation of neuropeptide B and neuropeptide W in the dipalmitoylphosphatidylcholine membrane bilayer

GPR7 and GPR8 are recently deorphanized G-protein-coupled receptors that are implicated in the regulation of neuroendocrine function, feeding behavior, and energy homeostasis. Neuropeptide B (NPB) and neuropeptide W (NPW) are two membrane-bound hypothalamic peptides, which specifically antagonize GPR7 and GPR8. Despite years of research, an accurate estimation of structure and molecular recognition of these neuropeptide systems still remains elusive. Herein, we investigated the structure, orientation, and interaction of NPB and NPW in a dipalmitoylphosphatidylcholine bilayer using long-range molecular dynamics (MD) simulation. During 30-ns simulation, membrane-embedded helical axes of NPB and NPW tilted 30 and 15°, respectively, from the membrane normal in order to overcome possible hydrophobic mismatch with the lipid bilayer. The calculation of various structural parameters indicated that NPW is more rigid and compact as compared to NPB. Qualitatively, the peptides exhibited flexible N-terminal (residues 1–12) and rigid C-terminal α-helical parts (residues 13–21), confirming previous NMR data. A strong electrostatic attraction between C-termini and headgroup atoms caused translocation of the peptides towards lower leaflet of the bilayer. The stabilizing hydrogen bonds (H-bonds) between phosphate groups and Trp1, Lys3, and Arg15 of the peptides played important roles for membrane anchoring. MD simulations of Alanine (Ala) mutants revealed that WYK->Ala variant of NPB/NPW lacked crucial H-bond interactions with phospholipid headgroups and also caused severe misfolding in NPB. Altogether, the knowledge of preferred structural fold and interaction of neuropeptides within the membrane bilayer will be useful to develop synthetic agonist or antagonist peptides for GPR7 and GPR8.     

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 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 studies of lipid bilayer systems

The main structural element of biological membranes is a liquid-crystalline lipid bilayer. Other constituents, i.e. proteins, sterols and peptides, either intercalate into or loosely attach to the bilayer. We applied a molecular dynamics simulation method to study membrane systems at various levels of compositional complexity. The studies were started from simple lipid bilayers containing a single type phosphatidylcholine (PC) and water molecules (PC bilayers). As a next step, cholesterol (Chol) molecules were introduced to the PC bilayers (PC-Chol bilayers). These studies provided detailed information about the structure and dynamics of the membrane/water interface and the hydrocarbon chain region in bilayers built of various types of PCs and Chol. This enabled studies of membrane systems of higher complexity. They included the investigation of an integral membrane protein in its natural environment of a PC bilayer, and the antibacterial activity of magainin-2. The latter study required the construction of a model bacterial membrane which consisted of two types of phospholipids and counter ions. Whenever published experimental data were available, the results of the simulations were compared with them.     

Molecular dynamics simulation of GM1 in phospholipid bilayer

We report molecular dynamics simulation of fully hydrated lipid bilayer of dimyristoyl phosphatidyl choline (DMPC) at room temperature with ganglioside GM1 attached to it in the upper layer under periodic boundary conditions. The simulation results indicate that the presence of a single GM1 molecule has local effects on the bilayer. Three sugar residues (GalNAc-Gal-Glc) of the pentasaccharide head group of GM1 remain on the lipid surface where as the NeuNAc residue extends out in the aqueous layer. The radial distribution functions suggest ordering of water molecules near the glycerol and carboxyl group of the sialic acid in the upper layer. One of the ceramide chains of GM1, the sphingosine chain, folds up and is stacked under the sugar residues lying on the surface. The other ceramide chain is inserted into the lipid bilayer. The arrangement of the polar head group as well as the acyl chains of the lipids which are immediate neighbours of the GM1 are modified compared to the non-neighbour ones and others at the lower layer. The time average conformation of GM1-pentasaccharide is stabilized by a number of inter residue hydrogen bonds that were observed experimentally. The trajectory average conformation of GM1-pentasaccharide was docked on to the cholera toxin molecule and the minimized complex reveals alternative binding modes between the toxin and the GM1-pentasaccharide moiety. The results of these simulation studies might help to understand the structure and nature of the effects of GM1 on the membrane at atomic resolution.     

Molecular dynamics of the KcsA K(+) channel in a bilayer membrane

Molecular dynamics (MD) simulations of an atomic model of the KcsA K(+) channel embedded in an explicit dipalmitoylphosphatidylcholine (DPPC) phospholipid bilayer solvated by a 150 mM KCl aqueous salt solution are performed and analyzed. The model includes the KcsA K(+) channel, based on the recent crystallographic structure of, Science. 280:69-77), 112 DPPC, K(+) and Cl(-) ions, as well as over 6500 water molecules for a total of more than 40,000 atoms. Three K(+) ions are explicitly included in the pore. Two are positioned in the selectivity filter on the extracellular side and one in the large water-filled cavity. Different starting configurations of the ions and water molecules in the selectivity filter are considered, and MD trajectories are generated for more than 4 ns. The conformation of KcsA is very stable in all of the trajectories, with a global backbone root mean square (RMS) deviation of less than 1.9 A with respect to the crystallographic structure. The RMS atomic fluctuations of the residues surrounding the selectivity filter on the extracellular side of the channel are significantly lower than those on the intracellular side. The motion of the residues with aromatic side chains surrounding the selectivity filter (Trp(67), Trp(68), Tyr(78), and Tyr(82)) is anisotropic with the smallest RMS fluctuations in the direction parallel to the membrane plane. A concerted dynamic transition of the three K(+) ions in the pore is observed, during which the K(+) ion located initially in the cavity moves into the narrow part of the selectivity filter, while the other two K(+) ions move toward the extracellular side. A single water molecule is stabilized between each pair of ions during the transition, suggesting that each K(+) cation translocating through the narrow pore is accompanied by exactly one water molecule, in accord with streaming potential measurements (, Biophys. J. 55:367-371). The displacement of the ions is coupled with the structural fluctuations of Val(76) and Gly(77), in the selectivity filter, as well as the side chains of Glu(71), Asp(80), and Arg(89), near the extracellular side. Thus the mechanical response of the channel structure at distances as large as 10-20 A from the ions in the selectivity filter appears to play an important role in the concerted transition.     

Molecular dynamics of a vasopressin V2 receptor in a phospholipid bilayer membrane

Molecular dynamics simulations were carried out for a V2 receptor (V2R) model embedded in a dimyristoylphosphatidylcholine (DMPC) bilayer. Both free and ligand-bound states of V2R were modeled. Our initial V2R model was obtained using a rule-based automated method for GPCR modeling and refined using constrained simulated annealing in vacuo. The docking site of the native vasopressin ligand was selected and justified upon consideration of ligand-receptor interactions and structure-activity data. The primary purpose of this work was to investigate the usefulness of MD simulation of an integral membrane protein like a GPCR receptor, upon inclusion of a carefully parameterized surrounding lipid membrane and water. Physical properties of the system were evaluated and compared with the fully hydrated pure DMPC bilayer membrane. The solvation interactions, individual lipid-protein interaction and fluctuations of the protein, the lipid, and water were analyzed in detail. As expected, the membrane-spanning helices of the protein fluctuate less than the peripheral loops do. The protein appears to disturb the local lipid structure. Simulations were carried out using AMBER 4.1 package upon constant number-pressure-temperature (NPT) conditions on massively parallel computers Cray T3E and IBM SP2.     

Molecular dynamics simulation of lipid reorientation at bilayer edges

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

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

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

Molecular dynamics simulation of lipid packing and mobility in bilayer membranes

A model of lipid bilayer membrane in water has been developed. Parameters have been selected that allow molecular dynamics simulation of lipid bilayers in the all-atom approximation. The calculated indices of packing and mobility of lipid molecules for the liquid crystalline state of the bilayer agree well with the experimental data. Based on the model of the liquid crystalline state of the membrane, a system in the gel-like state has been constructed. The gel-state model reproduces well the packing of lipids in real bilayers, whereas the mobility of molecules proves to be overestimated.     

Molecular dynamics simulation of dipalmitoylphosphatidylserine bilayer with Na+ counterions

We performed a molecular dynamics simulation of dipalmitoylphosphatidylserine (DPPS) bilayer with Na+ counterions. We found that hydrogen bonding between the NH group and the phosphate group leads to a reduction in the area per headgroup when compared to the area in dipalmitoylphosphatidylcholine bilayer. The Na+ ions bind to the oxygen in the carboxyl group of serine, thus giving rise to a dipolar bilayer similar to dipalmitoylphosphatidylethanolamine bilayer. The results of the simulation show that counterions play a crucial role in determining the structural and electrostatic properties of DPPS bilayer.     

Glycophorin-induced cholesterol-phospholipid domains in dimyristoylphosphatidylcholine bilayer vesicles.

Glycophorin has been incorporated into unilamellar cholesterol-containing dimyristoylphosphatidylcholine vesicles that were reconstituted by the freeze and thaw technique. Evidence was obtained for a protein-induced structural reorganization of these mixed membranes. By differential scanning calorimetry, we were able to construct a phase diagram for the phospholipid/cholesterol mixture consisting of a liquid-ordered, a solid-ordered, and a liquid-disordered phase. Glycophorin at low molar fractions (XG less than 3 X 10(-3)) increases the relative amount of lipid in the liquid-ordered phase, which is interpreted as an enrichment of cholesterol in the vicinity of the protein. The formation of such steroid-enriched domains could be demonstrated directly by electron paramagnetic resonance using a spin-labeled cholesterol analogue. A drastic increase of the spin-spin interaction of the labeled steroid was observed in the presence of glycophorin.     

Transmembrane localization of cis-isomers of zeaxanthin in the host dimyristoylphosphatidylcholine bilayer membrane

The effects of the 9-cis and 13-cis isomers of zeaxanthin on the molecular organization and dynamics of dimyristoylphosphatidylcholine (DMPC) membranes were investigated using conventional and saturation recovery EPR observations of the 1-palmitoyl-2-(14-doxylstearoyl)phosphatidylcholine (14-PC) spin label. The results were compared with the effects caused by the all-trans isomer of zeaxanthin. Effects on membrane fluidity, order, hydrophobicity, and the oxygen transport parameter were monitored at the center of the fluid phase DMPC membrane. The local diffusion-solubility product of oxygen molecules (oxygen transport parameter) in the membrane center, studied by saturation-recovery EPR, decreased by 47% and 27% by including 10 mol% 13-cis and 9-cis zeaxanthin, respectively; whereas, incorporation of all-trans zeaxanthin decreased this parameter by only 11%. At a zeaxanthin-to-DMPC mole ratio of 1:9, all investigated isomers decreased the membrane fluidity and increased the alkyl chain order in the membrane center. They also increased the hydrophobicity of the membrane interior. The effects of these isomers of zeaxanthin on the membrane properties mentioned above increase as: all-trans<9-cis相似文献   

Molecular dynamics simulation of packing and mobility of lipids in bilayer membranes

A model of DSPC lipid membrane in gel and liquid-crystalline states has been developed. The parameters have been determined that enable one to calculate the molecular dynamics of lipid bilayers in the full-atromic approximation. The parameters of packing and mobility of lipid molecules for the liquid crystalline state of the bilayer have been calculated. The values agree well with experimental data. Based on the model of the liquid crystalline state of the membrane, a system in the gel-like state has been constructed. The model of the gel-like state reproduces well the packing of lipids in real bilayers, whereas the mobility of molecules in the gel-like state was found to be overestimated.     

Molecular dynamics of leucine and dopamine transporter proteins in a model cell membrane lipid bilayer

The dopamine transporter (DAT) operates via facilitated diffusion, harnessing an inward Na⁺ gradient to drive dopamine from the extracellular synaptic cleft to the neuron interior. The DAT is relevant to central nervous system disorders such as Parkinson disease and attention‐deficit hyperactivity disorder and is the primary site of action for the abused psychostimulants cocaine and amphetamines. Crystallization of a DAT homolog, the bacterial leucine transporter LeuT, provided the first reliable 3‐D DAT template. Here, the LeuT crystal structure and the DAT molecular model have been combined with their respective substrates, leucine and dopamine, in lipid bilayer molecular dynamics simulations toward tracking substrate movement along the protein's substrate/ion permeation pathway. Specifically, movement of residue pairs that comprise the “external gate” was followed as a function of substrate presence. The transmembrane (TM) 1 arginine‐TM 10 aspartate strut formed less readily in DAT compared with LeuT, with or without substrate present. For LeuT but not DAT, the addition of substrate enhanced the chances of forming the TM 1‐10 bridge. Also, movement of the fourth extracellular loop EL‐4 in the presence of substrate was more pronounced for DAT, the EL‐4 unwinding to a degree. The overall similarity between the LeuT and DAT molecular dynamics simulations indicated that LeuT was a legitimate model to guide DAT structure‐function predictions. There were, nevertheless, differences significant enough to allow for DAT‐unique insights, which may include how cocaine, methylphenidate (Ritalin, NIDA Drug Supply, Rockville, MD), and other DAT blockers are not recognized as substrates even though they can access the primary substrate binding pocket. Proteins 2010. © 2009 Wiley‐Liss, Inc.