Insights into the dynamics of DMSO in phosphatidylcholine bilayers.

The solvation effects of dimethyl sulfoxide (DMSO) on the phase stability of dimyristoylphosphatidylcholine (DMPC) have been fully characterized using differential scanning calorimetry (DSC) and fluorescence spectroscopy with 1,6-diphenyl-1,3,5-hexatriene (DPH). The temperatures of the sub-, pre-, and main transitions of DMPC were found to increase linearly with increasing mole fraction of DMSO up to mole fraction X=0.13 DMSO/H(2)O. Beyond X=0.13, the pre-transition peak started to merge with the peak representing the main transition. Simultaneously, the subtransition peak began to disappear as its transition temperature also decreased. At X=0.18, with both the subtransition and pre-transition absent, the main transition between the planar gel and the liquid-crystalline phase was observed at 30.3 degrees C. Transition enthalpy values indicated that the subgel, planar gel and rippled gel phases are most stable at X=0.11, 0.16 and 0.20 DMSO/H(2)O, respectively. This demonstrates that DMSO exerts distinct effects on each respective phase and corresponding transition. Temperature-dependent fluorescence emission scans show an increase in hydration as the system proceeds from the subgel phase all the way to the liquid-crystalline phase and correlated well with the effects of DMSO on the transition temperatures of DMPC observed in our calorimetry data. Initial observations for the sub- and main transition are further confirmed by fluorescence anisotropy using DPH as a probe. The results illustrate the differences in the microviscosity of each phase and how DMSO affects the phase transitions. Ultimately, our results suggest the most likely mechanism governing the biological actions of DMSO may involve the regulation of the solvation effects of water on the phospholipid bilayer.     

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

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

Molecular packing in 1-hexanol-DMPC bilayers studied by molecular dynamics simulation

The structure and molecular packing density of a "mismatched" solute, 1-hexanol, in lipid membranes of dimyristoyl phosphatidylcholine (DMPC) was studied by molecular dynamics simulations. We found that the average location and orientation of the hexanol molecules matched earlier experimental data on comparable systems. The local density or molecular packing in DMPC-hexanol was elucidated through the average Voronoi volumes of all heavy (non-hydrogen) atoms. Analogous analysis was conducted on trajectories from simulations of pure 1-hexanol and pure (hydrated) DMPC bilayers. The results suggested a positive volume change, DeltaV(m), of 4 cm(3) mol(-1) hexanol partitioned at 310 K in good accordance with experimental values. Analysis of the apparent volumes of each component in the pure and mixed states further showed that DeltaV(m) reflects a balance between a substantial increase in the packing density of the alcohol upon partitioning and an even stronger loosening in the packing of the lipid. Furthermore, analysis of Voronoi volumes along the membrane normal identifies a distinctive depth dependence of the changes in molecular packing. The outer (interfacial) part of the lipid acyl chains (up to C8) is stretched by about 4%. Concomitantly, the average lateral area per chain decreases and these two effects compensate so that the overall packing density in the outer region, where the hexanol molecules are located, remains practically constant. The core of the bilayer (C9-C13) is slightly thinned. The average lateral area per chain in this region expands, resulting in a looser packing density. The net effect in the core is a 2-3% decrease in density corresponding to a total volume increase of approximately 14 cm(3) mol(-1) hexanol partitioned.     

Folding study of an Aib-rich peptide in DMSO by molecular dynamics simulations.

To evaluate the ability of molecular dynamics (MD) simulations using atomic force-fields to correctly predict stable folded conformations of a peptide in solution, we show results from MD simulations of the reversible folding of an octapeptide rich in alpha-aminoisobutyric acid (2-amino-2-methyl-propanoic acid, Aib) solvated in di-methyl-sulfoxide (DMSO). This solvent generally prevents the formation of secondary structure, whereas Aib-rich peptides show a high propensity to form secondary structural elements, in particular 3(10)- and alpha-helical structures. Aib is, moreover, achiral, so that Aib-rich peptides can form left- or right-handed helices depending on the overall composition of the peptide, the temperature, and the solvation conditions. This makes the system an interesting case to study the ensembles of peptide conformations as a function of temperature by MD simulation. Simulations involving the folding and unfolding of the peptide were performed starting from two initial structures, a right-handed alpha-helical structure and an extended structure, at three temperatures, 298 K, 340 K, and 380 K, and the results are compared with experimental nuclear magnetic resonance (NMR) data measured at 298 K and 340 K. The simulations generally reproduce the available experimental nuclear Overhauser effect (NOE) data, even when a wide range of conformations is sampled at each temperature. The importance of adequate statistical sampling in order to reliably interpret the experimental data is discussed.     

Disparate molecular dynamics of plasmenylcholine and phosphatidylcholine bilayers

The molecular dynamics of binary dispersions of plasmenylcholine/cholesterol and phosphatidylcholine/cholesterol were quantified by electron spin resonance (ESR) and deuterium magnetic resonance (2H NMR) spectroscopy. The order parameter of both 5-doxylstearate (5DS) and 16-doxylstearate (16DS) was larger in vesicles comprised of plasmenylcholine in comparison to phosphatidylcholine at all temperatures studied (e.g., S = 0.592 vs. 0.487 for 5DS and 0.107 vs. 0.099 for 16DS, respectively, at 38 degrees C). Similarly, the order parameter of plasmenylcholine vesicles was larger than that of phosphatidylcholine vesicles utilizing either spin-labeled phosphatidylcholine or spin-labeled plasmenylcholine as probes of molecular motion. The ratio of the low-field to the midfield peak height in ESR spectra of 16-doxylstearate containing moieties (i.e., spin-labeled plasmenylcholine and phosphatidylcholine) was lower in plasmenylcholine vesicles (0.93 +/- 0.01) in comparison to phosphatidylcholine vesicles (1.03 +/- 0.01). 2H NMR spectroscopy demonstrated that the order parameter of plasmenylcholine was greater than that of phosphatidylcholine for one of the two diastereotopic deuterons located at the C-2 carbon of the sn-2 fatty acyl chain. The spin-lattice relaxation times for deuterated plasmenylcholine and phosphatidylcholine in binary mixtures containing 0-50 mol % cholesterol varied nonmonotonically as a function of cholesterol concentration and were different for each phospholipid subclass. Taken together, the results indicate that the vinyl ether linkage in the proximal portion of the sn-1 aliphatic chain of plasmenylcholine has substantial effects on the molecular dynamics of membrane bilayers both locally and at sites spatially distant from the covalent alteration.     

简述调控线粒体形态变化的分子机制

线粒体是细胞内高度动态变化的细胞器,其在细胞内不断运动、融合、分裂并形成动态平衡的网状结构。线粒体的融合和分裂是由多种蛋白精确调控完成。Mfns／Fz01P控制线粒体外膜的融合,而Mgmlp／OPA1则参与线粒体内膜的融合;Dnm1p／Drp1、Fis1p／Fis1和Mdv1p介导线粒体分裂的调控。线粒体形态对于细胞维持正常生理代谢和机体发育起着重要的作用,一旦调控出现障碍会导致严重的疾病。     

Acetylcholine receptor pore permeability studied by molecular dynamics simulation

A dynamic model of the closed-state pore of an acetylcholine receptor (five M2 α-helices stabilized with a (CH₂)₁₀₅ ring) is used to examine the migration of uncharged and charged probe particles equivalent to a hexahydrated sodium ion (van der Waals diameter 7.27 Å) propelled by varied external force along the channel axis. Ion movement through the pore is hindered by steric constraints and electrostatic interactions. The van der Waals gate is formed by helix residues 13′ (A-Val255, B-Val261, C-Val269, D-Val255, and E-Ile264), whereas the negatively charged residues in the upper part of the channel are important for ion selectivity.     

Effect of magnesium sulfate in oxidized lipid bilayers properties by using molecular dynamics

Magnesium sulfate (MgSO₄) has been used as a protector agent for many diseases related to oxidative stress. The effect of MgSO₄ on the oxidized lipid bilayer has not yet been studied using molecular dynamics calculations. In this work, the effects of oxidation were evaluated by using a POPC membrane model at different concentrations of its aldehyde (-CHO) and hydroperoxide (-OOH) derivatives with and without MgSO₄. Several quantitative and qualitative properties were evaluated, such as membrane thickness, area per lipid, area compressibility modulus, snapshots after simulation finish, density distributions, time evolutions of oxidized group positions, and radial distributions of oxidized group concerning Mg. Results indicate that in the absence of MgSO₄ the mobility of oxidized groups, particularly –CHO, toward the surface interface is high. At a low oxidation level of the bilayer there is an increase in the compressibility modulus as compared to the unoxidized bilayer. MgSO₄, at a low oxidation level, tends to lessen the oxidation effects by lowering the dispersion in the distribution of oxidized species toward the membrane surface and the water region. However, MgSO₄ does not change the trends of decreasing membrane thickness and area compressibility modulus and increasing area per lipid upon oxidation. In this regard, MgSO₄ diminishes the electrostatic long-distance attractive interactions between the oxidized groups and the charged headgroups of the interface, owing to the Mg⁺² and SO₄^-2 screening effects and an electrostatic stabilization of the headgroups, preventing the pore formation, which is well-known to occur in oxidized membranes.     

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

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

Structural characterization of cyclic kallidin analogues in DMSO by nuclear magnetic resonance and molecular dynamics.

The conformational properties in DMSO of two head-to-tail cyclic analogues of kallidin ([Lys(0)]-bradykinin, KL) as well as those of the corresponding linear peptides were studied by NMR and molecular dynamics (MD) simulations. The modifications in the sequence were introduced at position 6, resulting in the four peptides, [Tyr(6)]-KL (YKL), [Trp(6)]-KL (WKL), cyclo-([Tyr(6)]-KL) (YCKL) and cyclo-([Trp(6)]-KL) (WCKL).The linear WKL analogue was significantly more potent than kallidin on rat duodenum preparations, whereas YKL was significantly less potent. Both cyclic peptides, YCKL and WCKL displayed similar activity, lower than that of the linear analogues and also of cyclo-KL.The two linear analogues display high conformational flexibility in DMSO. In the predominant conformer, for both peptides, all three X-Pro bonds adopt a trans configuration. Three out of four conformers present in YCKL and WCKL were completely assigned. The configurations at the X-Pro bonds are the same for the two analogues. All cyclic conformers show a cis configuration in at least one X-Pro bond and always opposite configuration for the two consecutive X-Pro bonds.The NOE-restrained MD calculations resulted in the detection of several elements of secondary structure in each of the conformers. Such elements are described and their possible relevance to biological activity is discussed.     

Binding and insertion of alpha-helical anti-microbial peptides in POPC bilayers studied by molecular dynamics simulations

We have performed molecular dynamics simulations of the interactions of two alpha-helical anti-microbial peptides, magainin2 and its synthetic analog of MSI-78, with palmitoyl-oleoyl-phosphatidylcholine (POPC) lipid bilayers. We used various initial positions and orientations of the peptide with respect to the lipid bilayer, including a surface-bound state parallel to the interface, a trans-membrane state, and a partially inserted state. Our 20 ns long simulations show that both magainin2 and MSI-78 are most stable in the lipid environment, with the peptide destabilized to different extents in both aqueous and lipid/water interfacial environments. We found that there are strong specific interactions between the lysine residues of the peptides and the lipid head-group regions. MSI-78, owing to its large number of lysines, shows better binding characteristics and overall stability when compared to magainin2. We also find that both peptides destabilize the bilayer environment, as observed by the increase in lipid tail disorder and the induction of local curvature on the lipid head-groups by the peptides. From all the simulations, we conclude that the hydrogen bonding interactions between the lysines of the peptides and the oxygens of the polar lipid head-groups are the strongest and determine the overall peptide binding characteristics to the lipids.     

Interpretation of small angle X-ray measurements guided by molecular dynamics simulations of lipid bilayers

Reconstruction and interpretation of lipid bilayer structure from X-ray scattering often rely on assumptions regarding the molecular distributions across the bilayer. It is usually assumed that changes in head-head spacings across the bilayer, as measured from electron density profiles, equal the variations in hydrocarbon thicknesses. One can then determine the structure of a bilayer by comparison to the known structure of a lipid with the same headgroup. Here we examine this procedure using simulated electron density profiles for the benchmark lipids DMPC and DPPC. We compare simulation and experiment in both real and Fourier space to address two main aspects: (i) the measurement of head-head spacings from relative electron density profiles, and (ii) the determination of the absolute scale for these profiles. We find supporting evidence for the experimental procedure, thus explaining the robustness and consistency of experimental structural results derived from electron density profiles. However, we also expose potential pitfalls in the Fourier reconstruction that are due to the limited number of scattering peaks. Volumetric analysis of simulated bilayers allows us to propose an improved, yet simple method for scale determination. In this way we are able to remove some of the restrictions imposed by limited scattering data in constructing reliable electron density profiles.     

Naratriptan aggregation in lipid bilayers: perspectives from molecular dynamics simulations

In order to understand the interaction between naratriptan and a fully hydrated bilayer of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidyl-choline (POPC), we carried out molecular dynamics simulations. The simulations were performed considering neutral and protonated ionization states, starting from different initial conditions. At physiological pH, the protonated state of naratriptan is predominant. It is expected that neutral compounds could have larger membrane partition than charged compounds. However, for the specific case of triptans, it is difficult to study neutral species in membranes experimentally, making computer simulations an interesting tool. When the naratriptan molecules were originally placed in water, they partitioned between the bilayer/water interface and water phase, as has been described for similar compounds. From this condition, the drugs displayed low access to the hydrophobic environment, with no significant effects on bilayer organization. The molecules anchored in the interface, due mainly to the barrier function of the polar and oriented lipid heads. On the other hand, when placed inside the bilayer, both neutral and protonated naratriptan showed self-aggregation in the lipid tail environment. In particular, the protonated species exhibited a pore-like structure, dragging water through this environment.

Open image in new window

Graphical Abstract Different behaviour of Naratriptan and Sumatriptan, when the drugs were originally placed in the lipid core

     

Hemolytic mechanism of dioscin proposed by molecular dynamics simulations

Saponins are a class of compounds containing a triterpenoid or steroid core with some attached carbohydrate modules. Many saponins cause hemolysis. However, the hemolytic mechanism of saponins at the molecular level is not yet fully understood. In an attempt to explore this issue, we have studied dioscin—a saponin with high hemolytic activity—through extensive molecular dynamics (MD) simulations. Firstly, all-atom MD simulations of 8 ns duration were conducted to study the stability of the dioscin–cholesterol complex and the cholesterol–cholesterol complex in water and in decane, respectively. MM-GB/SA computations indicate that the dioscin–cholesterol complex is energetically more favorable than the cholesterol–cholesterol complex in a non-polar environment. Next, several coarse-grained MD simulations of 400 ns duration were conducted to directly observe the distribution of multiple dioscin molecules on a DPPC-POPC-PSM-CHOL lipid bilayer. Our results indicate that dioscin can penetrate into the lipid bilayer, accumulate in the lipid raft micro-domain, and then bind cholesterol. This leads to the destabilization of lipid raft and consequent membrane curvature, which may eventually result in the hemolysis of red cells. This possible mechanism of hemolysis can well explain some experimental observations on hemolysis.     

Modulation of glycophorin A transmembrane helix interactions by lipid bilayers: molecular dynamics calculations

Starting from the glycophorin A dimer structure determined by NMR, we performed simulations of both dimer and monomer forms in explicit lipid bilayers with constant normal pressure, lateral area, and temperature using the CHARMM potential. Analysis of the trajectories in four different lipids reveals how lipid chain length and saturation modulate the structural and energetic properties of transmembrane helices. Helix tilt, helix-helix crossing angle, and helix accessible volume depend on lipid type in a manner consistent with hydrophobic matching concepts: the most relevant lipid property appears to be the bilayer thickness. Although the net helix-helix interaction enthalpy is strongly attractive, analysis of residue-residue interactions reveals significant unfavorable electrostatic repulsion between interfacial glycine residues previously shown to be critical for dimerization. Peptide volume is nearly conserved upon dimerization in all lipid types, indicating that the monomeric helices pack equally well with lipid as dimer helices do with one another. Enthalpy calculations indicate that the helix-environment interaction energy is lower in the dimer than in the monomer form, when solvated by unsaturated lipids. In all lipid environments there is a marked preference for lipids to interact with peptide predominantly through one rather than both acyl chains. Although our trajectories are not long enough to allow a full thermodynamic treatment, these results demonstrate that molecular dynamics simulations are a powerful method for investigating the protein-protein, protein-lipid, and lipid-lipid interactions that determine the structure, stability and dynamics of transmembrane alpha-helices in membranes.     

Detergent-resistant, ceramide-enriched domains in sphingomyelin/ceramide bilayers

When cell membranes are treated with Triton X-100 or other detergents at 4 degrees C, a nonsolubilized fraction can often be recovered, the "detergent-resistant membranes", that is not found when detergent treatment takes place at 37 degrees C. Detergent-resistant membranes may be related in some cases to membrane "rafts". However, several basic aspects of the formation of detergent-resistant membranes are poorly understood. To answer some of the relevant questions, a simple bilayer composition that would mimic detergent-resistant membranes was required. The screening of multiple lipid compositions has shown that the binary mixture egg sphingomyelin/egg ceramide (SM/Cer) exhibits the required detergent resistance. In detergent-free membranes composed of different mixtures of SM and Cer (5-30 mol % of Cer) differential scanning calorimetry, fluorescence spectroscopy, and fluorescence microscopy experiments reveal the presence of discrete, Cer-enriched gel domains in a broad temperature range. In particular, at temperatures below SM phase transition ( approximately 40 degrees C) two gel (respectively Cer-rich and SM-rich) phases are directly observed using fluorescence microscopy. Although pure SM membranes are fully solubilized by Triton X-100 at room temperature, 5 mol % Cer is also enough to induce detergent resistance, even with a large detergent excess and lengthy equilibration times. Short-chain Cers do not give rise to detergent resistance. SM/Cer mixtures containing up to 30 mol % Cer become fully soluble at approximately 50 degrees C, i.e., well above the gel-fluid transition temperature of SM. The combined results of temperature-dependent solubilization and differential scanning calorimetry reveal that SM-rich domains are preferentially solubilized over the Cer-rich ones as soon as the former melt (i.e., at approximately 40 degrees C). As a consequence, at temperatures allowing only partial solubilization, the nonsolubilized residue is enriched in Cer with respect to the original bilayer composition. Fluorescence microscopy of giant unilamellar vesicles at room temperature clearly shows that SM-rich domains are preferentially solubilized over the Cer-rich ones and that the latter become more rigid and extensive as a consequence of the detergent effects. These observations may be relevant to the phenomena of sphingomyelinase-dependent signaling, generation of "raft platforms", and detergent-resistant cell membranes.     

Cationic DMPC/DMTAP lipid bilayers: molecular dynamics study

Cationic lipid membranes are known to form compact complexes with DNA and to be effective as gene delivery agents both in vitro and in vivo. Here we employ molecular dynamics simulations for a detailed atomistic study of lipid bilayers consisting of a mixture of cationic dimyristoyltrimethylammonium propane (DMTAP) and zwitterionic dimyristoylphosphatidylcholine (DMPC). Our main objective is to examine how the composition of the DMPC/DMTAP bilayers affects their structural and electrostatic properties in the liquid-crystalline phase. By varying the mole fraction of DMTAP, we have found that the area per lipid has a pronounced nonmonotonic dependence on the DMTAP concentration, with a minimum around the point of equimolar DMPC/DMTAP mixture. We show that this behavior has an electrostatic origin and is driven by the interplay between positively charged TAP headgroups and the zwitterionic phosphatidylcholine (PC) heads. This interplay leads to considerable reorientation of PC headgroups for an increasing DMTAP concentration, and gives rise to major changes in the electrostatic properties of the lipid bilayer, including a significant increase of total dipole potential across the bilayer and prominent changes in the ordering of water in the vicinity of the membrane. Moreover, chloride counterions are bound mostly to PC nitrogens implying stronger screening of PC heads by Cl ions compared to TAP headgroups. The implications of these findings are briefly discussed.     

Constant surface tension molecular dynamics simulations of lipid bilayers with trehalose

Surface areas and fluctuations evaluated from 50 ns molecular dynamics simulations of fully hydrated dipalmitoylphosphatidylcholine (DPPC) bilayers in a 1:2 trehalose:lipid ratio carried out at surface tensions 10, 17 and 25 dyn/cm/leaflet are compared with those of pure bilayers under the same conditions. Trehalose increases the surface area, as consistent with the surface tension lowering observed in simulations at constant area. The system bulk elastic modulus K _b  = 1.5 ± 0.3 × 10¹⁰ dyn/cm². It is independent of bilayer surface area and trehalose content within statistical error. In contrast, the area elastic modulus K _a shows a strong area dependence. At 64 Å²/lipid (the experimental surface area), K _a  = 138 ± 26 dyn/cm for a pure DPPC bilayer and 82 ± 10 dyn/cm for one with trehalose; i.e. trehalose increases fluidity of the bilayer surface at this area per lipid.