Molecular dynamics study on the stabilization of dehydrated lipid bilayers with glucose and trehalose

In this study, we use molecular dynamics simulations to investigate and compare the interactions of DPPC bilayers with and without saccharides (glucose or trehalose) under dehydrated conditions. Results from the simulations indicate that unilamellar bilayers lose their structural integrity under dehydrated conditions in the absence of saccharides; however, in the presence of either glucose or trehalose, the bilayers maintain their stability. Hydrogen bond analysis shows that the saccharide molecules displace a significant amount of water surrounding the lipid headgroups. At the same time, the additional hydrogen bonds formed between water and saccharide molecules help to maintain a hydration layer on the lipid bilayer interface. On the basis of the hydrogen bond distributions, trehalose forms more hydrogen bonds with the lipids than glucose, and it is less likely to interact with neighboring saccharide molecules. These results suggest that the interaction between the saccharide and lipid molecules through hydrogen bonds is an essential component of the mechanism for the stabilization of lipid bilayers.     

Molecular dynamics simulation study of interaction between model rough hydrophobic surfaces

We study some aspects of hydrophobic interaction between molecular rough and flexible model surfaces. The model we use in this work is based on a model we used previously (Eun, C.; Berkowitz, M. L. J. Phys. Chem. B 2009, 113, 13222-13228), when we studied the interaction between model patches of lipid membranes. Our original model consisted of two graphene plates with attached polar headgroups; the plates were immersed in a water bath. The interaction between such plates can be considered as an example of a hydrophilic interaction. In the present work, we modify our previous model by removing the charge from the zwitterionic headgroups. As a result of this procedure, the plate character changes: it becomes hydrophobic. By separating the total interaction (or potential of mean force, PMF) between plates into the direct and the water-mediated interactions, we observe that the latter changes from repulsive to attractive, clearly emphasizing the important role of water as a medium. We also investigate the effect of roughness and flexibility of the headgroups on the interaction between plates and observe that roughness enhances the character of the hydrophobic interaction. The presence of a dewetting transition in a confined space between charge-removed plates confirms that the interaction between plates is strongly hydrophobic. In addition, we notice that there is a shallow local minimum in the PMF in the case of the charge-removed plates. We find that this minimum is associated with the configurational changes that flexible headgroups undergo as the two plates are brought together.     

Molecular dynamics simulation study of the water-mediated interaction between zwitterionic and charged surfaces

We calculated the potential of mean force (PMF) for the interaction between a model zwitterionic bilayer and a model charged bilayer. To understand the role of water, we separated the PMF into two components: one due to direct interaction and the other due to water-mediated interaction. In our calculations, we observed that water-mediated interaction is attractive at larger distances and repulsive at shorter. The calculation of the entropic and enthalpic contributions to the solvent-mediated components of the PMF showed that attraction is entropically dominant, while repulsion is dominated by the enthalpy.     

Molecular dynamics simulation of an archaeal lipid bilayer with sodium chloride

We have performed molecular dynamics simulations of a bilayer formed by the synthetic archaeal lipid, diphytanyl phosphatidylcholine, in NaCl electrolyte solution at four different concentrations (0-4 M) to investigate how structural and dynamic properties of the model archaeal membrane are changed due to the ionic strength of the solution. The archaeal lipid bilayer shows minor changes in their physical properties, indicating the unusual high stability of the membrane against salt, though small reductions of molecular area and lateral diffusion of the lipid are detected at the highest electrolyte concentration of 4 M. Sodium ions penetrate to the ether-rich region, where the ions are likely bound to the ether oxygen in the sn-1 chain rather than to that in the sn-2 chain. The observed salt bridges among two or three neighboring lipids account for the small reduction in the molecular area. The bound ions together with the counter (chloride) ions give rise to a diffusive electric double layer; as a result, the membrane dipole potential is slightly increased with increasing NaCl concentration.     

The interaction of lipid modified pseudopeptides with lipid membranes

We have studied the structure of two lipopeptides based on the simple dipeptide building block L-Phe-D-Oxd. These peptides have been reported previously to form fiber-like materials. The lipopeptides synthesized here had the structures C(n)(2)H((2n+1))CO-L-Phe-D-Oxd-OBn or C(n)(2)H((2n+1))CO-D-Phe-L-Oxd-OBn with n = 5 or 11. Addition of the N-terminal lipid modification did not cause a major disturbance of the structures these molecules form. The lipid modifications themselves showed highly rigid structures as inferred from solid-state (2)H NMR. The peptide backbone showed (13)C NMR chemical shifts in agreement with β-sheet secondary structure. Addition of a lipid modification to the N-terminus is a common motif in biology to attach proteins to the membrane. Therefore, we also investigated the lipopeptides in the presence of synthetic POPC bilayers. Two different molecular species were detected under these circumstances: (i) lipopeptide monomers that showed chain order parameters similar to those of the host membrane, (ii) lipopeptide aggregates that exhibited very similar structures and dynamics as the crystalline aggregates. Overall, the lipopeptides showed a well defined and rigid secondary structure that is in agreement with fibrillar aggregates previously detected for those peptides without the lipid modification.     

Molecular dynamics simulation and linear interaction energy study of d-Glu-based inhibitors of the MurD ligase

The biosynthetic pathway of the bacterial peptidoglycan, where MurD is an enzyme involved at the intracellular stage of its construction, represents a collection of highly selective macromolecular targets for novel antibacterial drug design. In this study as part of our investigation of the MurD bacterial target two recently discovered classes of the MurD ligase inhibitors were investigated resulting from the lead optimization phases of the N-sulfonamide d-Glu MurD inhibitors. Molecular dynamics simulations, based on novel structural data, in conjunction with the linear interaction energy (LIE) method suggested the transferability of our previously obtained LIE coefficients to further d-Glu based classes of MurD inhibitors. Analysis of the observed dynamical behavior of these compounds in the MurD active site was supported by static drug design techniques. These results complement the current knowledge of the MurD inhibitory mechanism and provide valuable support for the d-Glu paradigm of the inhibitor design.     

Molecular dynamics simulation of sub-and supercritical water with a new interaction potential

A molecular dynamics experiment was performed for water under sub-and supercritical conditions with a new interaction potential including the nonelectrostatic H-bond component. The internal energy, self-diffusion coefficient, mean number of H-bonds per water molecule, and distributions of molecules according to the number of H-bonds were calculated over a wide temperature range at pressures of 50 and 100 MPa. The temperature dependences of these properties were analyzed.     

Molecular dynamics simulation of a lipid diamond cubic phase.

This paper presents the first atomistic simulation of a cubic membrane phase. Using the molecular dynamics simulation technique both the global and the local organization of glycerolmonoolein molecules inside the diamond cubic phase are studied. Multinanosecond simulations reveal that the center of the cubic bilayer remains close to the infinite periodic minimal surface that describes the diamond geometry. We further show that the equilibrium structure of the surfactant molecules inside the cubic phase is very similar to their structure inside a simulated lamellar bilayer. The small differences arise from the packing constraints of the surfactants within the cubic phase which has an area per surfactant that increases toward the bilayer center.     

Molecular dynamics simulation study on the interaction of KRN 7000 and three analogues with human CD1d

Many analogues of KRN 7000, a synthetic glycolipid (α-galactosylceramide) exhibiting immuno-stimulatory activity and antitumor properties, were previously synthesized and tested in order to understand the reasons for the resulting biological activity and Th1/Th2 cytokine profile. Principles have been established for the interaction of such glycolipids with the human CD1d molecule but the exact mechanism by which different ligands with the same polar head elicit distinct biological responses remains unclear. Based on these experiments and on the available crystal structures, protein-ligand interactions are explored using molecular dynamics simulations. Hydrogen bond interactions are examined with regard to the polar group orientation. The influence of modulations on the dynamic behavior of the CD1d-glycolipid complex is addressed. Overall, our data support the mechanism by which the shortening of the α-GalCer sphingosine chain causes a significant twist of the CD1d α1 helix structure from residue Phe84 that affects the position of CD1d residues involved in the TCR recognition.     

Molecular dynamics simulation study of superhydrated perdeuterated natrolite using a new interaction potential model

To test a new interaction potential, molecular dynamics simulations of zeolite natrolite were performed for the structures under ambient conditions hydrated by perdeuterated water and at high pressure (1.87 GPa) in the superhydrated phase, which were recently studied by neutron diffraction. The experimental structures were reproduced with reasonable accuracy, and the hydrogen bond features are discussed. As in ordinary natrolite, a flip motion of water molecules around the HOH bisector is found, which, together with translational oscillations, gives rise to transient hydrogen bonds between water molecules, which do not appear from experimental equilibrium coordinates. The dynamics of water molecules can explain some problems encountered in refining the experimental structure. Vibrational spectra of natrolite containing perdeuterated water, which are not yet measured, were simulated, and their qualitative trend is discussed.     

Investigation of interaction of Leu-enkephalin with lipid membranes

Enkephalins are peptides with morphine-like activity. To achieve their biological function, they must be transported from an aqueous phase to the lipid-rich environment of their membrane bound receptor proteins. In our study, zeta potential (ZP) method was used to detect the association of Leu-enkephalin and Leu-enkephalinamide with phospholipid liposomes constituted from egg-phosphatidylcholine (EPC), dioleoyl-phosphatidylethanolamine (DOPE), cholesterol (Chol), sphingomyelin (SM) as well as soybean phospholipid (SBPL). Transfer of the peptides over lipid membranes was examined by electrophysiology technique (ET) and fluorescence spectroscopy (FS), and further confirmed using 4-fluoro-7-nitrobenzofurazan (NBD-F) labeled Leu-enkephalin (NBD-F-enkephalin) with confocal laser scanning microscopy method (CLSM). Results of zeta potential showed that enkephalinamide associated with lipid membranes and gradually saturated on the membranes either hydrophobically or electrostatically or both. Data from electrophysiology technique indicated that Leu-enkephalin could cause transmembrane currents, suggesting the transfer of peptides across lipid membranes. Transfer examined by fluorescence spectroscopy implied that it could be separated into three steps, adsorption, transportation and desorption, which was afterward reaffirmed by confocal laser scanning microscopy. Transfer efficiencies of enkephalin across SBPL, EPC/DOPE, EPC/DOPE/SM, EPC/SM and EPC/Chol lipid bilayer membranes were evaluated with ET and CLSM experiments. Results showed that the addition of either sphingomyelin or cholesterol, or negatively charged lipid in lipid membrane composition could lower the transfer efficiency.     

Stabilization of bilayer lipid membranes on solid supports by trehalose

Using the electrostriction method the effect of the glucose and trehalose on the elasticity modulus perpendicular to the membrane plane, E_⊥, and the electrical capacitance, C, of supported bilayer lipid membranes (s-BLM) formed on the freshly cut tip of Teflon-coated Ag wire was studied. Addition of saccharides into the electrolyte resulted in a decrease in the elasticity modulus of the s-BLM formed from the soybean phosphatidylcholine in n-hexadecane, while the capacitance was increased. In addition, the trehalose has a considerable stabilizing effect on the above parameters of the s-BLM. Treatment of the s-BLM in an electrolyte containing 0.3 M of the trehalose allowed storage of the s-BLM under dry conditions and under refrigeration, with the subsequent recovery of membrane parameters after the wire had been dipped into the electrolyte.     

Interaction of zervamicin IIB with lipid bilayers. Molecular dynamics study

In this work we have studied the interaction of zervamicin IIB (ZrvIIB) with the model membranes of eukaryotes and prokaryotes using all-atom molecular dynamics. In all our simulations zervamicin molecule interacted only with lipid headgroups but did not penetrate the hydrophobic core of the bilayers. During the interaction with the prokaryotic membrane zervamicin placed by its N-termini towards the lipids and rotated at an angle of 40° relatively to the bilayer surface. In the case of eukaryotic membrane zervamicin stayed in the water and located parallel to the membrane surface. We compared hydrogen bonds between peptide and lipids and concluded that interactions of ZrvIIB with prokaryotic membrane are stronger than those with eukaryotic one. Also it was shown that two zervamicin molecules formed dimer and penetrated deeper in the area of lipid headgroups.     

Molecular dynamics simulation study of the dynamics of fluids at solid-liquid interfaces

The dynamics of fluids at solid-liquid interfaces is investigated. In particular, we consider a simple Lennard-Jones fluid as well as a melt of hexadecane chains. For the Lennard-Jones fluid, the numerical results are compared with analytical calculations based on the diffusion equation, which shows that the numerical results can very well by described by the solution of the diffusion equation for reflecting surfaces. The diffusion coefficient is practically independent of the position within the film, although the fluid is inhomogeneous perpendicular to the surface. In contrast, the dynamics of the centers of mass of hexadecane molecules perpendicular to repulsive surfaces is severely slowed down due to their extended and anisotropic nature and cannot be described by a single particle diffusion equation.     

The effect of cationic gemini surfactants upon lipid membranes. An experimental and molecular dynamics simulation study

Gemini surfactants possess interesting interfacial and aggregation properties that have prompted comprehensive studies and successful applications in a wide variety of fields. However, a systematic study on the effect of gemini tail and spacer length upon the organization of lipid membranes has not been presented so far. In this study, we analyze the action of dicationic alkylammonium bromide gemini surfactants on DPPC liposomes, the latter employed as a model of lipid membranes. Differential scanning calorimetry results indicate that the surfactants presenting shorter tails (12 carbons) induce a decrease in the overall order of the bilayer, while those with longer tails (16 and 18 carbons) lead to the formation of more ordered structures. The respective influence on the degree of lipid order transverse to the bilayer was additionally studied resorting to a detailed fluorescence anisotropy study. In this case, it is observed that among the shorter tail surfactants, those with longer spacers (6 and 10 carbons) are responsible for a more pronounced disrupting effect upon the membrane, especially close to the lipid polar heads. Molecular dynamics simulation supports the most important findings and provides insight into the mechanism that governs this interaction. Accordingly, the interplay between tail and spacer length accounts for the differential vertical positioning of the gemini molecules and atom-density in the core of the bilayer, that provide a rationale for the experimental observations.     

Molecular dynamics study of perovskite structures with modified interatomic interaction potentials

The structure of compounds with the perovskite structure ABX₃ (A and B are cations, X are anions O^2—, F^—, Cl^—, Br^—, and I^—), which are widely used in engineering due to unique electrical, optical, and photovoltaic properties, has been considered. Hybrid organic—inorganic halide perovskites important for photovoltaics of a new generation are worth mentioning; they contain cations of organic nitrogen bases as monovalent cations. A molecular dynamics (MD) study of the CaTiO₃ base structure (Ca²⁺, Ti⁴⁺, and O^2—) has been performed in order to develop the methodology of computer simulation and optimization of the shape and parameters of atomic potentials for perovskite systems.     

Molecular dynamics simulation on the interaction between polymer inhibitors and anhydrite surface

Investigation on the microscopic interaction between polymer inhibitors and calcium sulfate will be helpful for understanding its scale inhibition mechanism and can provide a theoretical guidance to developing new scale inhibitors. In this work, molecular dynamics simulations with COMPASS force field have been performed to simulate the interaction between hydrolyzed polymaleic anhydride (HPMA), polyaspartic acid (PASP), polyepoxysuccinic acid (PESA), polyacrylic acid (PAA) and the (001) and (020) surfaces of anhydrite (AD) crystal with and without water. The results show that the sequence of binding energies between four polymer inhibitors and AD (001) and (020) with water is PESA > PASP > HPMA > PAA. The binding energy of the same polymer inhibitor on AD (001) is smaller than that on AD (020). Water molecules weaken the deformations of HPMA and PAA but aggravate those of PASP and PESA. Natural bond orbital (NBO) charges of the repeat units of polymer inhibitors were calculated by B3LYP/6‐31G* method. The Coulomb interaction is formed between the O atoms of polymer inhibitors and the Ca atoms of AD crystal. The system of polymer–AD is mainly contributed from the non‐bonding interaction. Polymer inhibitors do not interact directly with AD crystal, but indirectly through the interactions between inhibitor–H₂O and H₂O–AD, i.e. water molecules participate in scale inhibition of polymer inhibitors to AD crystal. Water molecules cannot be ignored when the interaction models are constructed, i.e. solvent effect cannot be ignored. Copyright © 2015 John Wiley & Sons, Ltd.     

Molecular dynamics study on the interaction of a mithramycin dimer with a decanucleotide duplex

The complex of a minor groove binding drug mithramycin (MTR) and the self-complementary d(TAGCTAGCTA) 10-mer duplex was investigated by molecular dynamics (MD) simulations using the AMBER 7.0 suite of programs. There is one disaccharide and trisaccharide segment projecting from opposite ends of an aglycone chromophore of MTR. A MTR dimer complex (MTR)2Mg2+ is formed in the presence of a coordinated ion Mg2+. A NMR solution structure of two (MTR)2Mg2+ complexes bound with one DNA duplex, namely, the 2:1 duplex complex, was taken as the starting structure for the MD simulation. The partial charge on each atom was calculated using the multiple-RESP fitting procedure, and all of the missing parameters in the Parm99 force field used were adapted comparably from the literature. The length of the MD simulation was 5 ns, and the binding free energy for the formation of a 1:1 or 2:1 duplex complex was determined from the last 4 ns of the simulation. The binding free energies were decomposed to components of the contributions from different energy types, and the changes in the helical parameters of the bound DNA duplex plus the glycosidic linkages between sugar residues of the bound MTR dimer were determined. It was found that binding of the first (MTR)2Mg2+ complex with the DNA duplex to form a 1:1 duplex complex does not cause stiffening of the duplex especially in the unoccupied site of the duplex. However, the overall flexibility of the DNA duplex is reduced substantially once the second (MTR)2Mg2+ complex is bound with the unoccupied site to form the 2:1 duplex complex. The van der Waals interactions were found to be dominant in the central part of the DNA duplex where sugar residues from each bound (MTR)2Mg2+ complex were inwardly pointing and the corresponding minor groove was widened.     

Molecular dynamics simulation with stochastically constrained pressure

The authors propose a new algorithm for molecular dynamics simulation. The method includes a Monte Carlo scheme for incrementing the dilation rate in the equations of motion. The new algorithm needs no extra computation and the dynamics of the system preserves its continuity. Application of this approach is very advantageous for models where the derivation and the computation of the pressure is time consuming. The authors present results of model calculations.     

Molecular dynamics simulation study on the structural stabilities of polyglutamine peptides

It is known that Huntington's disease patients commonly have glutamine (Q) repeat sequences longer than a critical length in the coding area of Huntingtin protein in their genes. As the polyglutamine (polyQ) region becomes longer than the critical length, the disease occurs and Huntingtin protein aggregates, both in vitro and in vivo, as suggested by experimental and clinical data. The determination of polyglutamine structure is thus very important for elucidation of the aggregation and disease mechanisms. Here, we perform molecular dynamics calculations on the stability of the structure based on the β-helix structure suggested by Perutz et al. (2002) [Perutz, M.F., Finch, J.T., Berriman, J., Lesk, A., 2002. Amyloid fibers are water-filled nanotubes. Proc. Natl. Acad. Sci. USA 99, 5591]. We ensure that perfect hydrogen bonds are present between main chains of the β-helix based on the previous studies, and perform simulations of stretches with 20, 25, 30, 37 and 40 glutamine residues (20Q, 25Q, 30Q, 37Q and 40Q) for the Perutz models with 18.5 and 20 residues per turn (one coil). Our results indicate that the structure becomes more stable with the increase of repeated number of Q, and there is a critical Q number of around 30, above which the structure of the Perutz model is kept stable. In contrast to previous studies, we started molecular dynamics simulations from conformations in which the hydrogen bonds are firmly formed between stacked main chains. This has rendered the initial β-helix structures of polyQ much more stable for longer time, as compared to those proposed previously. Model calculations for the initial structures of polyQ dimer and tetramer have also been carried out to study a possible mechanism for aggregation.