Methanol in water solution at 300 K

Monte Carlo simulations are presented for a methanol molecule enclosed in a cluster of 198 water molecules. We have simulated three conformations for methanol (staggered, perpendicular and eclipsed); in a fourth simulation the methanol's hydrogen atom in the hydroxyl group is free to assume any conformation. The water structure near the methyl group is similar to that reported for methane.     

Molecular dynamics simulations of bilayers containing mixtures of sphingomyelin with cholesterol and phosphatidylcholine with cholesterol

We performed six molecular dynamics simulations: three on hydrated bilayers containing pure phospholipids and three on hydrated bilayers containing mixtures of these phospholipids with cholesterol. The phospholipids in our simulations were SSM (sphingomyelin containing a saturated 18:0 acyl chain), OSM (sphingomyelin with an unsaturated 18:1 acyl chain), and POPC (palmitoyloleoylphosphatidylcholine containing one saturated and one unsaturated chain). Data from our simulations were used to study systematically the effect of cholesterol on phospholipids that differed in their headgroup and tail composition. In addition to the structural analysis, we performed an energetic analysis and observed that energies of interaction between cholesterol and neighboring SM molecules are similar to the energies of interaction between cholesterol and POPC. We also observed that the interaction energy between cholesterol and neighboring lipids cannot be used for the determination of which lipids are involved in the creation of a complex.     

Simulation of Polymer Aggregates as a Method to Investigate the Solvent Effect on Blend Complexation and the Relative Strength of H‐Bonding Groups in Blends

Simulations of polymer‐solvent and polymer‐polymer aggregates, in which the study of hydrogen bonding plays an important role, have been carried out with two blend systems. The aim was to examine the influence of the solvent on blend complexation and to compare the strength of different hydrogen bonds in a blend system. We quantified the strength of one hydrogen bond in the blend environments. For this we used the EVOCAP software, developed by our institute. It allows the building of large molecular aggregates with realistic and homogeneous densities, with an implemented positioning algorithm of the molecules under consideration and their excluded volume, and a charge equilibration method for the partial charge calculation. In the simulated aggregates the specific interaction energy of the hydrogen atom forming the hydrogen bond was a useful indicator for our studies. Through a direct correlation of this specific‐interaction energy with the strength of the hydrogen bond, we supported the experimental result that, in toluene, complex formation between poly(methyl methacrylate) (PMMA) and PSOH, a hydroxyl‐modified polystyrene, is possible, but not in tetrahydrofuran. Varying the proton‐donor polymer, also a hydroxyl‐modified polystyrene, in blends of poly(vinyl methyl ether) (PVME) with groups of different donor strength, we reconstructed the experimental row of increasing hydrogen‐bond strengths.     

Nonintercalating nanosubstrates create asymmetry between bilayer leaflets

The physical properties of lipid bilayers can be remodeled by a variety of environmental factors. Here we investigate using molecular dynamics simulations the specific effects of nanoscopic substrates or external contact points on lipid membranes. We expose palmitoyl-oleoyl phosphatidylcholine bilayers unilaterally and separately to various model nanosized substrates differing in surface hydroxyl densities. We find that a surface hydroxyl density as low as 10% is sufficient to keep the bilayer juxtaposed to the substrate. The bilayer interacts with the substrate indirectly through multiple layers of water molecules; however, despite such buffered interaction, the bilayers exhibit certain properties different from unsupported bilayers. The substrates modify transverse lipid fluctuations, charge density profiles, and lipid diffusion rates, although differently in the two leaflets, which creates an asymmetry between bilayer leaflets. Other properties that include lipid cross-sectional areas, component volumes, and order parameters are minimally affected. The extent of asymmetry that we observe between bilayer leaflets is well beyond what has been reported for bilayers adsorbed on infinite solid supports. This is perhaps because the bilayers are much closer to our nanosized finite supports than to infinite solid supports, resulting in a stronger support-bilayer electrostatic coupling. The exposure of membranes to nanoscopic contact points, therefore, cannot be considered as a simple linear interpolation between unsupported membranes and membranes supported on infinite supports. In the biological context, this suggests that the exposure of membranes to nonintercalating proteins, such as those belonging to the cytoskeleton, should not always be considered as passive nonconsequential interactions.     

Self-consistent mean-field model based on molecular dynamics: application to lipid-cholesterol bilayers

We have developed a dynamic self-consistent mean-field model, based on molecular-dynamics simulations, to study lipid-cholesterol bilayers. In this model the lipid bilayer is represented as a two-dimensional lattice field in the lipid chain order parameters, while cholesterol molecules are represented by hard rods. The motion of rods in the system is continuous and is not confined to lattice cells. The statistical mechanics of chain ordering is described by a mean field derived from an extension of a model due to Marcelja. The time evolution of the system is governed by stochastic equations. The ensemble of chain configurations required in partition sums, and the energies of interaction, are taken from atomistic level molecular-dynamics simulations of lipid bilayers. The model allows us to simulate systems 500 nm in lateral size for 20 micros time scales, or greater. We have applied the model to dipalmitoyl-phosphatidylcholine-cholesterol (Chol) bilayers at 50 degrees C for Chol concentrations between 2% and 33%. At low concentrations of Chol (2%-4%), the model predicts the formation of isolated clusters of Chol surrounded by relatively ordered lipid chains, randomly dispersed in the disordered bilayer. With increasing Chol composition, regions of Chol-induced order begin to overlap. Starting from about 11% Chol this ordering effect becomes system wide and regions unaffected by Chol are no longer detectable. From the analysis of properties of the model we conclude that the change in lipid chain order with increasing Chol concentration is continuous over the 20-mus scale of the simulations. We also conclude that at 50 degrees C no large-scale Chol-rich and Chol-depleted coexisting phase-separated regions form at any concentration. At no point in any of the simulations do we observe a higher degree of lateral organization, such as Chol-based superlattice structures.     

Design and synthesis of sphingomyelin-cholesterol conjugates and their formation of ordered membranes

A lipid raft is a cholesterol (Chol)-rich microdomain floating in a sea of lipid bilayers. Although Chol is thought to interact preferentially with sphingolipids such as sphingomyelin (SM), rather than with glycerophospholipids, the origin of the specific interaction has remained unresolved, primarily because of the high mobility of lipid molecules and weak intermolecular interactions. In this study, we synthesized SM-Chol conjugates with functionally designed linker portions to restrain Chol mobility and examined their formation of ordered membranes by a detergent insolubility assay, fluorescence anisotropy experiments, and fluorescence-quenching assay. In all of the tests, membranes prepared from the conjugates showed properties of ordered domains comparable to a SM-Chol (1:1) membrane. To gain insight into the structure of bilayers composed from the conjugates, we performed molecular dynamics simulations with 64 molecules of the conjugates, which suggested that the conjugates form a stable bilayer structure by bending at the linker portion and, mostly, reproduce the hydrogen bonds between the SM and Chol portions. These results imply that the molecular recognition between SM and Chol in an ordered domain is essentially reproduced by the conjugated molecules and, thus, demonstrates that these conjugate molecules could potentially serve as molecular probes for understanding molecular recognition in lipid rafts.     

A molecular-dynamics study of lipid bilayers: effects of the hydrocarbon chain length on permeability

In this paper, we investigate the effects of the hydrocarbon chain length of lipid molecules on the permeation process of small molecules through lipid bilayers. We perform molecular-dynamics simulations using three kinds of lipid molecules with different chain length: dilauroylphosphatidylcholine, dimyristoylphosphatidylcholine, and dipalmiltoylphosphatidylcholine. Free-energy profiles of O2, CO, NO, and water molecules are calculated by means of the cavity insertion Widom method and the probability ratio method. We show that the lipid membrane with longer chains has a larger and wider energy barrier. The local diffusion coefficients of water across the bilayers are also calculated by the force autocorrelation function method and the velocity autocorrelation function method. The local diffusion coefficients in the bilayers are not altered significantly by the chain length. We estimate the permeability coefficients of water across the three membranes according to the solubility-diffusion model; we find that the water permeability decreases modestly with increasing chain length of the lipid molecules.     

Hygroscopic growth of self-assembled layered surfactant molecules at the interface between air and organic salts

We report here the self-assembly of surfactant molecules at the interface of air and the hygroscopic quaternary ammonium salt tetrabutylammonium acetate (TBAAc). Homogeneously dissolved surfactant molecules at 100 degrees C self-assemble upon contacting air due to high moisture adsorption by the organic salt when cooling down. Highly ordered lamellar phases with different lattice spacings have been observed when surfactants with various lengths of alkyl chains were used. C(n)TMAB/TBAAc systems showed all-trans conformation of interior methylene carbons and interdigited bilayers with an average CH2 increment of 0.119 nm, while C(n)NH2/TBAAc systems showed trans/gauche mixed conformations of interior methylene carbons and bilayers with an average CH2 increment of 0.247 nm. C(n)NH2s in C(n)NH2/TBAAc formed bilayers through water-mediated intermolecular hydrogen bonds with a water layer thickness of 0.51-0.61 nm. In C(n)TAB/TBAAc, as the head group of C(n)TAB is bigger, the interdigited bilayer thickness (d-spacing) is smaller, because the bigger head groups accommodate enough space for alkyl tails to come in between them.     

Structure and dynamics of water at the interface with phospholipid bilayers

We have performed two molecular-dynamics simulations to study the structural and dynamical properties of water at the interface with phospholipid bilayers. In one of the simulations the bilayer contained neutral phospholipid molecules, dioleoylphosphatidylcholine (DOPC); in the second simulation the bilayer contained charged lipid molecules, dioleoylphosphatidylserine (DOPS). From the density profile of water we observe that water next to the DOPS bilayer is more perturbed as compared to water near the DOPC bilayer. Using an energetic criterion for the determination of hydrogen bonding we find that water molecules create strong hydrogen bonds with the headgroups of the phospholipid molecules. Due to the presence of these bonds and also due to the confinement of water, the translational and orientational dynamics of water at the interface are slowed down. The degree of slowing down of the dynamics depends upon the location of water molecules near a lipid headgroup.     

Evidence for hydroxyl radical generation during lipid (linoleate) peroxidation

The autoxidation of methyl linoleate in benzene at 37 degrees C by peroxyl radicals was found to generate hydroxyl radicals (.OH) from a secondary oxidation mechanism. The yield of hydroxyl radicals (approximately 2%) was determined by trapping these reactive radicals with benzene to give phenol. We propose that alphaC-H hydrogen abstraction from lipid hydroperoxides, the main autoxidation products, is the source of hydroxyl radicals.     

Liquid-vapor interfacial properties of water-ammonia mixtures: dependence on ammonia concentration

The equilibrium and dynamical properties of the liquid-vapor interfaces of water-ammonia mixtures are investigated by means of molecular-dynamics simulations. Altogether, we have simulated seven different systems of different concentration of ammonia. The inhomogeneous density, anisotropic orientational profiles, surface tension, and the pattern of hydrogen bonding are calculated for both water and ammonia molecules in order to characterize the location, width, thermodynamic aspects, and microscopic structure of the liquid-vapor interfaces of each of the water-ammonia systems. The dynamical aspects of the interfaces are investigated in terms of the anisotropic diffusion and dipole orientational relaxation of water and ammonia molecules. The properties of the interfaces are compared with those of the corresponding bulk phases. The present theoretical results are also compared with experimental findings wherever available.     

Neutron diffraction and computer simulation studies of D-xylose

Neutron diffraction with isotopic substitution (NDIS) experiments and molecular dynamics (MD) simulations have been used to examine the pentose D-xylose in aqueous solution. By specifically labeling D-xylose molecules with a deuterium atom at the nonexchangeable hydrogen position on C4, it was possible to extract information about the atomic structuring around just that specific position. The MD simulations were found to give satisfactory agreement with the experimental NDIS results and could be used to help interpret the scattering data in terms of the solvent structuring as well as the intramolecular hydroxyl conformations. Although the experiment is challenging and on the limit of modern instrumentation, it is possible by careful analysis, in conjunction with MD studies, to show that the conformation trans to H4 at 180 degrees is strongly disfavored, in excellent agreement with the MD results. This is the first attempt to use NDIS experiments to determine the rotameric conformation of a hydroxyl group.     

Molecular dynamics investigation of dynamical properties of phosphatidylethanolamine lipid bilayers

We describe the dynamic behavior of a 1-stearoyl-2-oleoyl-phosphatidylethanolamine (SOPE) bilayer from a 20 ns molecular dynamics simulation. The dynamics of individual molecules are characterized in terms of (2)H spin-lattice relaxation rates, nuclear overhauser enhancement spectroscopy (NOESY) cross-relaxation rates, and lateral diffusion coefficients. Additionally, we describe the dynamics of hydrogen bonding through an analysis of hydrogen bond lifetimes and the time evolution of clusters of hydrogen bonded lipids. The simulated trajectory is shown to be consistent with experimental measures of internal, intermolecular, and diffusive motion. Consistent with our analysis of SOPE structure in the companion paper, we see hydrogen bonding dominating the dynamics of the interface region. Comparison of (2)H T(1) relaxation rates for chain methylene segments in phosphatidylcholine and phosphatidylethanolamine bilayers indicates that slower motion resulting from hydrogen bonding extends at least three carbons into the hydrophobic core. NOESY cross-relaxation rates compare well with experimental values, indicating the observed hydrogen bonding dynamics are realistic. Calculated lateral diffusion rates (4 +/ -1 x 10(-8) cm(2)s) are comparable, though somewhat lower than, those determined by pulsed field gradient NMR methods.     

Interaction of hexadecylbetainate chloride with biological relevant lipids

The present work investigates the interaction of hexadecylbetainate chloride (C(16)BC), a glycine betaine-based ester with palmitoyl-oleoyl-phosphatidylcholine (POPC), sphingomyelin (SM), and cholesterol (CHOL), three biological relevant lipids present in the outer leaflet of the mammalian plasma membrane. The binding affinity and the mixing behavior between the lipids and C(16)BC are discussed based on experimental (isothermal titration calorimetry (ITC) and Langmuir film balance) and molecular modeling studies. The results show that the interaction between C(16)BC and each lipid is thermodynamically favorable and does not affect the integrity of the lipid vesicles. The primary adsorption of C(16)BC into the lipid film is mainly governed by a hydrophobic effect. Once C(16)BC is inserted in the lipid film, the polar component of the interaction energy between C(16)BC and the lipid becomes predominant. Presence of CHOL increases the affinity of C(16)BC for membrane. This result can be explained by the optimal matching between C(16)BC and CHOL within the film rather by a change of membrane fluidity due to the presence of CHOL. The interaction between C(16)BC and SM is also favorable and gives rise to highly stable monolayers probably due to hydrogen bonds between their hydrophilic groups. The interaction of C(16)BC with POPC is less favorable but does not destabilize the mixed monolayer from a thermodynamic point of view. Interestingly, for all the monolayers investigated, the exclusion surface pressures are above the presumed lateral pressure of the plasma membranes suggesting that C(16)BC would be able to penetrate into mammalian plasma membranes in vivo. These results may serve as a useful basis in understanding the interaction of C(16)BC with real membranes.     

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.     

Computational study on the properties and structure of methyl lactate

A theoretical study on the properties and molecular level structure of the very important green solvent methyl lactate is carried out in the gas phase and methanol and water solutions, with the solvent treated both explicitly and as a continuum. Torsional barriers giving rise to different conformers by rotation of the hydroxyl and methyl groups were analyzed using density functional theory (DFT) to establish the most stable conformer both in gas phase and solution. DFT computations on lactate dimers were also done to study short-range features, and the effect of the surrounding solvent on intra- and intermolecular hydrogen bonding was analyzed according to the polarizable continuum model approach. We have also studied lactate/water and lactate/methanol small clusters together with the corresponding binding energies. Moreover, classical molecular dynamics simulations (MD) were carried out to study medium- and large-range effects at lower computational cost. MD simulations at different pressure and temperature conditions on pure lactate were carried out, and mixtures with water and methanol of different compositions were also studied. Structural information, analyzed through the radial distribution functions, together with dynamic aspects of pure and mixed fluids were considered. The intramolecular hydrogen bonding ability of methyl lactate together with the possibility of homo- and hetero-intermolecular association determines the behavior of this molecule in pure fluids or in mixed.     

α-生育酚在模型生物膜中的分子动力学模拟

用分子动力学方法模拟了280, 310和350 K下α-生育酚在二豆蔻酰磷脂酰胆碱、二豆蔻酰磷脂酰乙醇胺、二硬脂酰磷脂酰胆碱和二硬脂酰磷脂酰乙醇胺双层膜中的性质, 包括了空间位置、氢键、取向和动力学性质, 取得了如下的结论. 第一, 生育酚头部的羟基一般位于脂双层亲疏水界面的下方, 升高温度将促进羟基向膜双层的中心移动, 在350 K时观察到了在上下两个单层间的翻转. 第二, 生育酚主要与磷脂的酯基形成氢键, 几乎不与磷脂酰乙醇胺的氨基形成氢键; 比较生育酚与磷脂酰胆碱和乙醇胺形成的氢键后发现, 后者更稳定. 第三, 生育酚的头部在膜中取向多变, 与膜的法线夹角不固定, 尾部的构象也很复杂. 第四, 在温度较低时, 生育酚的侧向扩散系数与磷脂的相当, 但在350 K时其扩散速度明显加快; 在垂直方向生育酚的扩散速度很慢.     

Conformational analysis of beta-glycosidic linkages in 13C-labeled glucobiosides using inter-residue scalar coupling constants

Four beta-linked glucobioses selectively (13)C labeled at C1' or C2' have been prepared. The inter-residue coupling constants, J(CH), and J(CC), have been determined and related to the solution conformations of the disaccharides using Karplus-type relationships. Relying only on the experimental coupling constants, glycosidic linkage conformation in methyl alpha-sophoroside (methyl 2-O-beta-D-glucopyranosyl-alpha-D-glucopyranoside), methyl alpha-laminarabioside (methyl 3-O-beta-D-glucopyranosyl-alpha-D-glucopyranoside), and methyl alpha-cellobioside (methyl 4-O-beta-D-glucopyranosyl-alpha-D-glucopyranoside) were found to be close to those observed in the solid state (39 degrees < phi(H) < 41 degrees , -24 degrees < psi(H) < -36 degrees ). The laminarabioside and cellobioside were found to have conformations that accommodate an intramolecular hydrogen bond to O5' that is observed in the solid state. In all compounds, the exocyclic hydroxymethyl groups retain a conformation close to that observed in unsubstituted glucose (gt/gg 1:1). Methyl alpha-gentiobioside (methyl 6-O-beta-D-glucopyranosyl-alpha-D-glucopyranoside) shows greater flexibility at the psi-torsion than the other disaccharides, but the population distribution around the C5-C6 bond is essentially unaffected by substitution. None of the O2' hydroxyl groups of the beta-D-glucopyranosyl residues in any of the disaccharides appear to be involved in inter-residue hydrogen bonding since (1)JCH, (1)JCC, and (2)JCH values sensitive to C2'-O2' rotamer distribution remain close to those observed in methyl beta-D-glucopyranoside.     

Peculiarities of lateral diffusion of lipids in three-component bilayers

The lateral diffusion of lipid molecules in macroscopically oriented bilayers of mixed dioleoyl phosphatidylcholine (DOPC), egg sphingomyeline (SM), and cholesterol (CHOL) and its dependence on cholesterol concentration and temperature was studied by NMR with pulsed field gradient. The system forms a lamellar liquid crystalline (LC) phase; in a certain range of temperatures and concentrations of cholesterol the system is separated into two subphases: a disordered LC phase (l_d) enriched with DOPC, and an ordered phase (l₀) enriched with SM. These are characterized by their own lateral diffusion coefficients (LDCs), which differ from one another by a factor of 1.5–5. The dependence of the LDCs in the phases on the cholesterol concentration was analyzed. There was no clear dependence for the disordered LC phase, but we found that LDCs tend to grow in the concentration range of 15–35 mol % of CHOL. This behavior could be due to the redistribution of lipid components as the concentration of CHOL increases, eventually leading to a rise in DOPC concentration in the lo phase. In the range of liquid-phase domains, we observed no dependence of LDCs on the diffusion time typical of the restricted diffusion regime, due to spatial restraints in the system. This could be associated with the relatively large size of the domains, and with the domain capability of lateral diffusion in a surrounding continuous phase.