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.     

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

Structure of dipalmitoylphosphatidylcholine/cholesterol bilayer at low and high cholesterol concentrations: molecular dynamics simulation.

By using molecular dynamics simulation technique we studied the changes occurring in membranes constructed of dipalmitoylphosphatidylcholine (DPPC) and cholesterol at 8:1 and 1:1 ratios. We tested two different initial arrangements of cholesterol molecules for a 1:1 ratio. The main difference between two initial structures is the average number of nearest-neighbor DPPC molecules around the cholesterol molecule. Our simulations were performed at constant temperature (T = 50 degrees C) and pressure (P = 0 atm). Durations of the runs were 2 ns. The structure of the DPPC/cholesterol membrane was characterized by calculating the order parameter profiles for the hydrocarbon chains, atom distributions, average number of gauche defects, and membrane dipole potentials. We found that adding cholesterol to membranes results in a condensing effect: the average area of membrane becomes smaller, hydrocarbon chains of DPPC have higher order, and the probability of gauche defects in DPPC tails is lower. Our results are in agreement with the data available from experiments.     

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     

Lysosomal membrane cholesterol dynamics

Although the majority of exogenous cholesterol and cholesterol ester enters the cell by LDL-receptor-mediated endocytosis and the lysosomal pathway, the assumption that cholesterol transfers out of the lysosome by rapid (minutes), spontaneous diffusion has heretofore not been tested. As shown herein, lysosomal membranes were unique among known organellar membranes in terms of cholesterol content, cholesterol dynamics, and response to cholesterol-mobilizing proteins. First, the lysosomal membrane cholesterol:phospholipid molar ratio, 0.38, was intermediate between those of the plasma membrane and other organellar membranes. Second, a fluorescence sterol exchange assay showed that the initial rate of spontaneous sterol transfer out of lysosomes and purified lysosomal membranes was extremely slow, t(1/2) >4 days. This was >100-fold longer than that reported in intact cells (2 min) and 40-60-fold longer than from any other known intracellular membrane. Third, when probed with several cholesterol-binding proteins, the initial rate of sterol transfer was maximally increased nearly 80-fold and the organization of cholesterol in the lysosomal membrane was rapidly altered. Nearly half of the essentially nonexchangeable sterol in the lysosomal membrane was converted to rapidly (t(1/2) = 6 min; fraction = 0.06) and slowly (t(1/2) = 154 min; fraction = 0.36) exchangeable sterol domains/pools. In summary, the data revealed that spontaneous cholesterol transfer out of the lysosome and lysosomal membrane was extremely slow, inconsistent with rapid spontaneous diffusion across the lysosomal membrane. In contrast, the very slow spontaneous transfer of sterol out of the lysosome and lysosomal membrane was consistent with cholesterol leaving the lysosome earlier in the endocytic process and/or with cholesterol transfer out of the lysosome being mediated by additional process(es) extrinsic to the lysosome and lysosomal membrane.     

Constant-pressure molecular dynamics investigation of cholesterol effects in a dipalmitoylphosphatidylcholine bilayer.

We report a 1.4-ns constant-pressure molecular dynamics simulation of cholesterol at 12.5 mol% in a dipalmitoylphosphatidylcholine (DPPC) bilayer at 50 degrees C and compare the results to our previous simulation of a pure DPPC bilayer. The interlamellar spacing was increased by 2.5 A in the cholesterol-containing bilayer, consistent with x-ray diffraction results, whereas the bilayer thickness was increased by only 1 A. The bilayer/water interface was more abrupt because the lipid headgroups lie flatter to fill spaces left by the cholesterol molecules. This leads to less compensation by the lipid headgroups of the oriented water contribution to the membrane dipole potential and could explain the experimentally observed increase in the magnitude of the dipole potential by cholesterol. Our calculations suggested that 12.5 mol% cholesterol does not significantly affect the conformations and packing of the hydrocarbon chains and produces only a slight reduction in the empty free volume. However, cholesterol has a significant influence on the subnanosecond time scale lipid dynamics: the diffusion constant for the center-of-mass "rattling" motion was reduced by a factor of 3, and the reorientational motion of the methylene groups was slowed along the entire length of the hydrocarbon chains.     

Molecular dynamics simulations of phospholipid bilayers with cholesterol

To investigate the microscopic interactions between cholesterol and lipids in biological membranes, we have performed a series of molecular dynamics simulations of large membranes with different levels of cholesterol content. The simulations extend to 10 ns, and were performed with hydrated dipalmitoylphosphatidylcholine (DPPC) bilayers. The bilayers contain 1024 lipids of which 0-40% were cholesterol and the rest DPPC. The effects of cholesterol on the structure and mesoscopic dynamics of the bilayer were monitored as a function of cholesterol concentration. The main effects observed are a significant ordering of the DPPC chains (as monitored by NMR type order parameters), a reduced fraction of gauche bonds, a reduced surface area per lipid, less undulations--corresponding to an increased bending modulus for the membrane, smaller area fluctuations, and a reduced lateral diffusion of DPPC-lipids as well as cholesterols.     

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 dynamics simulation of the structure of dimyristoylphosphatidylcholine bilayers with cholesterol, ergosterol, and lanosterol

Five molecular dynamics computer simulations were performed on different phospholipid:sterol membrane systems in order to study the influence of sterol structure on membrane properties. Three of these simulated bilayer systems were composed of a 1:8 sterol:phospholipid ratio, each of which employed one of the sterol molecules: cholesterol, ergosterol, and lanosterol. The two other simulations were of a bilayer with a 1:1 sterol:phospholipid ratio. These simulations employed cholesterol and lanosterol, respectively, as their sterol components. The observed differences in simulations with cholesterol and lanosterol may have their implication on the form of the phospholipid/sterol phase diagram.     

Molecular dynamics simulation of spontaneous membrane fusion during a cubic-hexagonal phase transition

We report a molecular dynamics simulation of the phase transition of monoolein from an inverted cubic phase to an inverted hexagonal phase. The transition proceeds via an intermediate structure consisting of water channels in a cubic geometry, in agreement with the predictions of the modified stalk theory (Siegel, 1999). Two mechanisms are identified by which the topology changes during the transition. Bilayer fusion proceeds via the formation of trans-monolayer contacts, whereas bilayer rupture is observed as a gradual thinning of each monolayer.     

Molecular dynamics simulation of a dye molecule in the interior of a bilayer: 1,6-diphenyl-1,3,5-hexatriene in dipalmitoylphosphatidylcholine

A molecular dynamics simulation was carried out for a dipalmitoylphosphatidylcholine (DPPC) membrane in its liquid crystalline state containing different concentrations of the dye molecule 1,6-diphenyl-1,3,5-hexatriene (DPH). From a numerical analysis of the trajectories, we obtained information concerning structural changes of the membrane due to the presence of the probe and some hydrodynamic information concerning the probe itself. The hydrodynamic properties regarding dye molecules that have been reported in this article are: rotational and translational diffusion coefficient and relaxation times. From this analysis, we estimated a range of values of 0.6-0.9 cP for the micro-viscosity in the mid-membrane. These simulations also afforded us some information regarding structural changes in the membrane as a consequence of the presence of the fluorescent dyes at different concentrations. Thus, the disorder inside the membrane, the surface area per lipid and thickness of the membrane were also investigated.     

Molecular dynamics simulation of the unfolding of individual bacteriorhodopsin helices in sodium dodecyl sulfate micelles

We report molecular dynamics simulations of the trends in the changes in secondary structure of the seven individual helices of bacteriorhodopsin when inserted into sodium dodecyl sulfate (SDS) micelles, and their dependence on the amino acid sequence. The results indicate that the partitioning of the helices in the micelles and their stability are dependent on the hydrophobicity of the transmembrane segments. Helices A, B, and E are stable and retain their initial secondary structure throughout the 100 ns simulation time. In contrast, helices C, D, F, and G show structural perturbations within the first 10 ns. The instabilities are localized near charged residues within the transmembrane segments. The overall structural instability of the helix is correlated with its partitioning to the surface of the micelle and its interaction with polar groups there. The in silico experiments were performed to complement the in vitro experiments that examined the partial denaturation of bacteriorhodopsin in SDS described in the preceding article (DOI 10.1021/bi201769z ). The simulations are consistent with the trends revealed by the experimental results but strongly underestimate the extent of helix to extended coil transformation. The reason may be either that the sampling time was not sufficiently long or, more interestingly, that interhelix residue interactions play a role in the unfolding of the helices.     

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 interaction between phospholipid membrane and pyrazine and its derivatives

We have performed a comparative molecular dynamics simulation of the diffusion process of the heterocyclic compound pyrazine and its methylated derivatives into the model membrane phospholipid bilayer. Several structural and dynamical bilayer parameters were measured, and qualitative interrelations between parameter changes and the substituted pyrazine structure were studied. The simulation results support the hypothesis that molecular mechanisms of biological effects of substituted pyrazines involve dissolution of the effector molecule in the membrane bilayer and subsequent changes in bilayer properties. This stage can provide the means for pyrazine molecules to interact with integral membrane proteins, directly or indirectly through the changed lipid environment of the protein.     

Molecular dynamics simulations on the interaction of the transmembrane NavAb channel with cholesterol and lipids in the membrane

Increased cholesterol levels are associated with multiple pathological conditions. In this work, molecular dynamics simulations were applied to observe the influence of membrane cholesterol levels on a voltage-gated sodium channel. Different lipid compositions are modeled around the channel to obtain information about the possible effects by which cholesterol influences Na_vAb channels. Cholesterol was normally not directly interacting with either the closed or inactivated conformation. Cholesterol increased lipid packing implying that it plays a crucial role in restricting lipid movement in the region around 1 nm of the channel in a 1-palmitoyl-2-oeleoyl phosphatidylcholine matrix. Our results provide the first computational indication of an indirect modulation of Na_vAb channels by membrane cholesterol.     

Molecular packing in steroid-lecithin monolayers,part II: Mixed films of cholesterol with dipalmitoylphosphatidylcholine and tetradecanoic acid

Mixed film studies of the systems cholesterol/tetradecanoic acid and cholesterol/dipalmitoylphosphatidylcholine have been carried out over the entire compositional range at 21°C. When compared on an acyl chain basis the condensing effects were found to be essentially independent of which host-lipid was utilized. The phase change of the host lipid was shifted to higher pressures, then broadened and eliminated. Maximal condensation occurred at just above 42 mol% for the cholesterol/DPPC system. In both systems the two components were initially found to be miscible at all proportions.The results are interpreted in terms of the molecular packing of cholesterol with acyl boundary layers, one significantly, one weakly affected. Maximum condensation is a result of packing that provides maximum cholesterol/acyl chain contact. Consideration is given to both long term stability of such mixed monolayers and the behaviour of the corresponding bilayer states.     

Molecular dynamics simulation of the tubulin dimer with cytostatics

Colchicin, podophylotoxin and indibulin are natural cytostatics that are used in the treatment of neoplasms. But applications of those compounds are rather restricted due to the high toxicity and low specificity. It seems very promising to investigate possibility to design new analogs of the above mentioned drugs that will possess higher cytostatic activity and less toxicity. For this purpose we see computer modeling experiments of tubulin and above mentioned compounds interaction as a powerful approach to design new artificial cytostatics with desired properties. In the current study the CHARMM software of protein-ligand interaction molecular dynamics method has been used. Particularly the following strategy has been applied. Molecules of the cytostatits have been positioned at several random positions around binding sites of tubulin and after energy minimization several binding sites have been identified on the tubulin macromolecule. In these binding sites structural changes that may be responsible for tubulin polymerization have been detected.