DISCUSSION PAPER: MOLECULAR STRUCTURES AND CONFORMATIONS OF THE PHOSPHOLIPIDS AND SPHINGOMYELINS

Insight into the Molecular Mechanism of Surface Interactions of Phosphatidylcholines─Langmuir Monolayer Study Complemented with Molecular Dynamics Simulations

Mutual interactions between components of biological membranes are pivotal for maintaining their proper biophysical properties, such as stability, fluidity, or permeability. The main building blocks of biomembranes are lipids, among which the most important are phospholipids (mainly phosphatidylcholines (PCs)) and sterols (mainly cholesterol). Although there is a plethora of reports on interactions between PCs, as well as between PCs and cholesterol, their molecular mechanism has not yet been fully explained. Therefore, to resolve this issue, we carried out systematic investigations based on the classical Langmuir monolayer technique complemented with molecular dynamics simulations. The studies involved systems containing 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) analogues possessing in the structure one or two polar functional groups similar to those of DPPC. The interactions and rheological properties of binary mixtures of DPPC analogues with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and cholesterol were compared with reference systems (DPPC/POPC and DPPC/cholesterol). This pointed to the importance of the ternary amine group in PC/cholesterol interactions, while in PC mixtures, the phosphate group played a key role. In both cases, the esterified glycerol group had an effect on the magnitude of interactions. The obtained results are crucial for establishing structure-property relationships as well as for designing substitutes for natural lipids.

Stacks of monogalactolipid bilayers can transform into a lattice of water channels

The lipid matrix of thylakoid membranes is a lamellar bilayer, but under a certain condition it can convert locally into a nonlamellar structure. This is possible because one of the main membrane lipids, MGDG, promotes the formation of an inverse hexagonal phase. Here, the spontaneous transformation of aligned hydrated MGDG bilayers into nonlamellar structures is investigated using all-atom molecular dynamics simulation. Previous studies have demonstrated that MGDG polar head groups connect vertically across the interface. In this study, the evolution of the system's initial structure into a lattice of water channels and contacted surfaces created by numerous vertical MGDG connections depended on the width of the hydrating water layers. These widths controlled the bilayers' ability to bend, which was a prerequisite for channel formation. Locally, an intensive exchange of MGDG molecules between apposing bilayer leaflets occurred, although a stable semi-toroidal stalk did not develop.

Lipid/water interface of galactolipid bilayers in different lyotropic liquid-crystalline phases

In this study, carried out using computational methods, the organisation of the lipid/water interface of bilayers composed of galactolipids with both α-linolenoyl acyl chains is analysed and compared in three different lyotropic liquid-crystalline phases. These systems include the monogalactosyldiglyceride (MGDG) and digalactosyldiglyceride (DGDG) bilayers in the lamellar phase, the MGDG double bilayer during stalk phase formation and the inverse hexagonal MGDG phase. For each system, lipid-water and direct and water-mediated lipid-lipid interactions between the lipids of one bilayer leaflet and those of two apposing leaflets at the onset of new phase (stalk) formation, are identified. A network of interactions between DGDG molecules and its topological properties are derived and compared to those for the MGDG bilayer.

Lutein and Zeaxanthin in the Lipid Bilayer-Similarities and Differences Revealed by Computational Studies

Lutein and zeaxanthin are two similar carotenoids of the xanthophyll subgroup. Carotenoids are synthesized almost entirely by plants but are also present in significant amounts in animals. They are essential components of the lipid matrix of biomembranes, and one of their functions is to protect cells from light radiation, free radicals and oxidative stress. Carotenoids, depending on their chemical structure, can locate at various positions and in different orientations in the bilayer. Xanthophylls (XAN) are polar and in the bilayer are positionally restricted. In the case of lutein and zeaxanthin, whose both ionone rings are hydroxy-substituted and as such are anchored in the lipid bilayer interfaces, the position is generally transmembrane. However, both experimental and computer modelling studies indicate that lutein can also locate horizontally below the bilayer interface. This location has never been observed for zeaxanthin. To find a molecular-level explanation for the difference in the orientations of the XAN molecules in the bilayer, a number of phosphatidylcholine-XAN bilayers were constructed and molecular dynamics (MD) simulated for 1.1 µs each. The all-trans XAN molecules were initially placed either parallel or perpendicular to the bilayer surface. With the exception of one lutein, the horizontally placed molecules adopted the transmembrane orientation within 100-600 ns. On the basis of detailed analyses of the XAN orientations and the numbers and lifetimes of their interactions in the bilayer, a plausible explanation is offered as to why a lutein molecule may remain in the horizontal orientation while zeaxanthin does not. Contrary to common believe, lutein horizontal orientation is not related to the ε-ring rotation around the C6'-C7' bond.

CHARMM36 Lipid Force Field with Explicit Treatment of Long-Range Dispersion: Parametrization and Validation for Phosphatidylethanolamine, Phosphatidylglycerol, and Ether Lipids

Long-range Lennard-Jones (LJ) interactions have been incorporated into the CHARMM36 (C36) lipid force field (FF) using the LJ particle-mesh Ewald (LJ-PME) method in order to remove the inconsistency of bilayer and monolayer properties arising from the exclusion of long-range dispersion [Yu, Y.; Semi-automated Optimization of the CHARMM36 Lipid Force Field to Include Explicit Treatment of Long-Range Dispersion. J. Chem. Theory Comput. 2021, 10.1021/acs.jctc.0c01326. (preceding article in this issue)]. The new FF is denoted C36/LJ-PME. While the first optimization was based on three phosphatidylcholines (PCs), this work extends the validation and parametrization to more lipids including PC, phosphatidylethanolamine (PE), phosphatidylglycerol (PG), and ether lipids. The agreement with experimental structure data is excellent for PC, PE, and ether lipids. C36/LJ-PME also compares favorably with scattering data of PG bilayers but less so with NMR deuterium order parameters of 1,2-dimyristoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (DMPG) at 303.15 K, indicating a need for future optimization regarding PG-specific parameters. Frequency dependence of NMR T1 spin-lattice relaxation times is well-described by C36/LJ-PME, and the overall agreement with experiment is comparable to C36. Lipid diffusion is slower than C36 due to the added long-range dispersion causing a higher viscosity, although it is still too fast compared to experiment after correction for periodic boundary conditions. When using a 10 Å real-space cutoff, the simulation speed of C36/LJ-PME is roughly equal to C36. While more lipids will be incorporated into the FF in the future, C36/LJ-PME can be readily used for common lipids and extends the capability of the CHARMM FF by supporting monolayers and eliminating the cutoff dependence.

Asymmetric Spontaneous Intercalation of Lutein into a Phospholipid Bilayer, a Computational Study

Lutein, a hydroxylated carotenoid, is a pigment synthesised by plants and bacteria. Animals are unable to synthesise lutein, nevertheless, it is present in animal tissues, where its only source is dietary intake. Both in plants and animals, carotenoids are associated mainly with membranes where they carry out important physiological functions. Their trafficking to and insertion into membranes are not well recognised due to experimental difficulties. In this paper, a computational approach is used to elucidate details of the dynamics and energetics of lutein intercalation from the water to the phospholipid bilayer phase. The dynamics is studied using molecular dynamics simulation, and the energetics using umbrella sampling. Lutein spontaneous insertion into the bilayer and translocation across it proceed via formation of hydrogen bonds between its hydroxyl groups and water and/or phospholipid oxygen atoms as well as desolvation of its polyene chain. As lutein molecule is asymmetric, its bilayer intercalation is also asymmetric. The course of events and timescale of the intercalation are different from those of helical peptides. The time of full lutein intercalation ranges from 20 to 100 ns and its final orientation is predominately vertical. Nevertheless, some lutein molecules are in the final horizontal position and some aggregate in the water phase and remain there for the whole simulation time. The highest energy barrier for the intercalation process is ~2.2 kcal/mol and the energy gain is ~18 kcal/mol. The results obtained for lutein can be applied to other xanthophylls and molecules of a similar structure.

•

Lutein as an amphiphilic molecule readily intercalates into a lipid bilayer.

•

MD simulations with dense sampling revealed details of the intercalation process.

•

Time of full lutein membrane intercalation ranges from 20 to 100 ns.

•

Lutein intercalation from the ε ring end is less probable than from the β ring end.

•

Horizontal position of lutein in the bilayer is less probable than vertical.

Lipid Configurations from Molecular Dynamics Simulations

The extent to which current force fields faithfully reproduce conformational properties of lipids in bilayer membranes, and whether these reflect the structural principles established for phospholipids in bilayer crystals, are central to biomembrane simulations. We determine the distribution of dihedral angles in palmitoyl-oleoyl phosphatidylcholine from molecular dynamics simulations of hydrated fluid bilayer membranes. We compare results from the widely used lipid force field of Berger et al. with those from the most recent C36 release of the CHARMM force field for lipids. Only the CHARMM force field produces the chain inequivalence with sn-1 as leading chain that is characteristic of glycerolipid packing in fluid bilayers. The exposure and high partial charge of the backbone carbonyls in Berger lipids leads to artifactual binding of Na+ ions reported in the literature. Both force fields predict coupled, near-symmetrical distributions of headgroup dihedral angles, which is compatible with models of interconverting mirror-image conformations used originally to interpret NMR order parameters. The Berger force field produces rotamer populations that correspond to the headgroup conformation found in a phosphatidylcholine lipid bilayer crystal, whereas CHARMM36 rotamer populations are closer to the more relaxed crystal conformations of phosphatidylethanolamine and glycerophosphocholine. CHARMM36 alone predicts the correct relative signs of the time-average headgroup order parameters, and reasonably reproduces the full range of NMR data from the phosphate diester to the choline methyls. There is strong motivation to seek further experimental criteria for verifying predicted conformational distributions in the choline headgroup, including the 31P chemical shift anisotropy and 14N and CD3 NMR quadrupole splittings.

Is the cholesterol bilayer domain a barrier to oxygen transport into the eye lens?

In the eye lens, the oxygen partial pressure is very low and the cholesterol (Chol) content in cell membranes is very high. Disturbance of these quantities results in cataract development. In human lens membranes, both bulk phospholipid-Chol domains and the pure Chol bilayer domains (CBDs) were experimentally detected. It is hypothesized that the CBD constitutes a significant barrier to oxygen transport into the lens. Transmembrane profiles of the oxygen diffusion-concentration product, obtained with electron paramagnetic resonance spin-labeling methods, allow evaluation of the oxygen permeability (P_M) of phospholipid membranes but not the CBD. Molecular dynamics simulation can independently provide components of the product across any bilayer domain, thus allowing evaluation of the P_M across the CBD. Therefore, to test the hypothesis, MD simulation was used. Three bilayers containing palmitoyl-oleoyl-phosphorylcholine (POPC) and Chol were built. The pure Chol bilayer modeled the CBD, the 1:1 POPC-Chol bilayer modeled the bulk membrane in which the CBD is embedded, and the POPC bilayer was a reference. To each model, 200 oxygen molecules were added. After equilibration, the oxygen concentration and diffusion profiles were calculated for each model and multiplied by each other. From the respective product profiles, the P_M of each bilayer was calculated. Favorable comparison with experimental data available only for the POPC and POPC-Chol bilayers validated these bilayer models and allowed the conclusion that oxygen permeation across the CBD is ~10 smaller than across the bulk membrane, supporting the hypothesis that the CBD is a barrier to oxygen transport into the eye lens.

Drude Polarizable Force Field for Molecular Dynamics Simulations of Saturated and Unsaturated Zwitterionic Lipids

Additive force fields are designed to account for induced electronic polarization in a mean-field average way, using effective empirical fixed charges. The limitation of this approximation is cause for serious concerns, particularly in the case of lipid membranes, where the molecular environment undergoes dramatic variations over microscopic length scales. A polarizable force field based on the classical Drude oscillator offers a practical and computationally efficient framework for an improved representation of electrostatic interactions in molecular simulations. Building on the first-generation Drude polarizable force field for the dipalmitoylphosphatidylcholine 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) molecule, the present effort was undertaken to improve this initial model and expand the force field to a wider range of phospholipid molecules. New lipids parametrized include dimyristoylphosphatidylcholine (DMPC), dilauroylphosphatidylcholine (DLPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), dipalmitoylphosphatidylethanolamine (DPPE), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE), and 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine (DOPE). The iterative optimization protocol employed in this effort led to lipid models that achieve a good balance between reproducing quantum mechanical data on model compound representative of phospholipids and reproducing a range of experimental condensed phase properties of bilayers. A parametrization strategy based on a restrained ensemble-maximum entropy methodology was used to help accurately match the experimental NMR order parameters in the polar headgroup region. All the parameters were developed to be compatible with the remainder of the Drude polarizable force field, which includes water, ions, proteins, DNA, and selected carbohydrates.

Computer modelling studies of the bilayer/water interface

This review summarises high resolution studies on the interface of lamellar lipid bilayers composed of the most typical lipid molecules which constitute the lipid matrix of biomembranes. The presented results were obtained predominantly by computer modelling methods. Whenever possible, the results were compared with experimental results obtained for similar systems. The first and main section of the review is concerned with the bilayer-water interface and is divided into four subsections. The first describes the simplest case, where the interface consists only of lipid head groups and water molecules and focuses on interactions between the lipid heads and water molecules; the second describes the interface containing also mono- and divalent ions and concentrates on lipid-ion interactions; the third describes direct inter-lipid interactions. These three subsections are followed by a discussion on the network of direct and indirect inter-lipid interactions at the bilayer interface. The second section summarises recent computer simulation studies on the interactions of antibacterial membrane active compounds with various models of the bacterial outer membrane. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

A computer model of a polyunsaturated monogalactolipid bilayer

1,2-di-O-acyl-3-O-β-D-galactopyranosyl-sn-glycerol (MGDG) is the main lipid component of thylakoid membranes of higher plants and algae. This monogalactolipid is thought of as a non-bilayer lipid but actually it can form both lamellar and nonlamellar phases. In this study, molecular dynamics (MD) simulations of the fully hydrated di-18:3 MGDG bilayer in the lamellar phase were carried out at 310 and 295 K for 200 and 450 ns, respectively, using the GROMACS 4 software package and OPLS-AA force field. At both temperatures, the lamellar phase of the systems was stable. The pure di-18:3 MGDG bilayer is the first step towards creating a computer model of the lipid matrix of the thylakoid membrane and the main aim of this study was to validate the computer model of di-18:3 MGDG in the bilayer and also to assess the properties of the bilayer. However, only a few of the predicted properties could be compared with those derived experimentally and in other MD simulations because of insufficient amount of such data. Thus, direct validation of the MGDG bilayer proved difficult. Therefore, in the validation process also an indirect approach was used, in which a computer model of the 1,2-dioleoyl-sn-glycero-3-phosphatidylcholine (DOPC) bilayer simulated at the same temperatures using the same force field as the MGDG bilayer was assessed. Successful validation of the DOPC bilayer parameterized in the OPLS-AA force field and similar properties of the MGDG molecules in the pure 18:3 MGDG and in binary 18:3 MGDG-PC bilayers indicate that the computer model of the MDGD molecule is faithful and the MGDG bilayer is representative on the time scales covered in these MD simulations.

Refined OPLS all-atom force field for saturated phosphatidylcholine bilayers at full hydration

We report parametrization of dipalmitoyl-phosphatidylcholine (DPPC) in the framework of the Optimized Parameters for Liquid Simulations all-atom (OPLS-AA) force field. We chose DPPC as it is one of the most studied phospholipid species and thus has plenty of experimental data necessary for model validation, and it is also one of the highly important and abundant lipid types, e.g., in lung surfactant. Overall, PCs have not been previously parametrized in the OPLS-AA force field; thus, there is a need to derive its bonding and nonbonding parameters for both the polar and nonpolar parts of the molecule. In the present study, we determined the parameters for torsion angles in the phosphatidylcholine and glycerol moieties and in the acyl chains, as well the partial atomic charges. In these calculations, we used three methods: (1) Hartree-Fock (HF), (2) second order Møller-Plesset perturbation theory (MP2), and (3) density functional theory (DFT). We also tested the effect of the polar environment by using the polarizable continuum model (PCM), and for acyl chains the van der Waals parameters were also adjusted. In effect, six parameter sets were generated and tested on a DPPC bilayer. Out of these six sets, only one was found to be able to satisfactorily reproduce experimental data for the lipid bilayer. The successful DPPC model was obtained from MP2 calculations in an implicit polar environment (PCM).

Comparative computer simulation study of cholesterol in hydrated unary and binary lipid bilayers and in an anhydrous crystal

Models created with molecular dynamics simulations are used to compare the organization and dynamics of cholesterol (Chol) molecules in three different environments: (1) a hydrated pure Chol bilayer that models the Chol bilayer domain, which is a pure Chol domain embedded in the bulk membrane; (2) a 2-palmitoyl-3-oleoyl-d-glycerol-1-phosphorylcholine bilayer saturated with cholesterol (POPC-Chol50) that models the bulk membrane; (3) a Chol crystal. The computer model of the hydrated pure Chol bilayer is stable on the microsecond time scale. Some structural characteristics of Chol molecules in the Chol bilayer are similar to those in the POPC-Chol50 bilayer (e.g., tilt of Chol rings and chains), while others are similar to those in Chol crystals (e.g., surface area per Chol, bilayer thickness). The key result of this study is that the Chol bilayer has, unexpectedly, a dynamic structure, with Chol mobility similar to that in the POPC-Chol50 bilayer though slower. This is the major difference compared to Chol crystals, where Chol molecules are immobile. Also, water accessibility to Chol-OH groups in the Chol bilayer is not limited. On average, each Chol molecule makes 2.3 hydrogen bonds with water in the Chol bilayer, compared with 1.7 hydrogen bonds in the POPC-Col50 bilayer.

Orientation and conformation of lipids in crystals of transmembrane proteins

Orientational order parameters and individual dihedral torsion angles are evaluated for phospholipid and glycolipid molecules that are resolved in X-ray structures of integral transmembrane proteins in crystals. The order parameters of the lipid chains and glycerol backbones in protein crystals are characterised by a much wider distribution of orientational order than is found in fluid lipid bilayers and reconstituted lipid-protein membranes. This indicates that the lipids that are resolved in crystals of membrane proteins are mostly not representative of the entire lipid-protein interface. Much of the chain configurational disorder of the membrane-bound lipids in crystals arises from C-C bonds in energetically disallowed skew conformations. This suggests configurational heterogeneity of the lipids at a single binding site: eclipsed conformations occur also in the glycerol backbone torsion angles and the C-C torsion angles of the lipid head groups. Conformations of the lipid glycerol backbone in protein crystals are not restricted to the gauche C1-C2 rotamers found invariably in phospholipid bilayer crystals. Lipid head-group conformations in the protein crystals also do not conform solely to the bent-down conformation, with gauche-gauche configuration of the phosphodiester, that is characteristic of phospholipid bilayer membranes. Stereochemical violations in the protein-bound lipids are evidenced by ester carboxyl groups in non-planar configurations, and even in the cis configuration. Some lipids have the incorrect enantiomeric configuration of the glycerol backbone, and many of the branched methyl groups in the phytanyl chains associated with bacteriorhodopsin have the incorrect S configuration.

Headgroup conformations of phospholipids from molecular dynamics simulation: sampling challenges and comparison to experiment

The preferred conformations of the glycerol region of a phospholipid have been explored using replica exchange molecular dynamics (MD) simulations and compared with the results of standard MD approaches and with experiment. We found that due to isomerization rates in key torsions that are slow on the timescale of atomistic MD simulations, standard MD is not able to produce accurate equilibrium conformer distributions from reasonable trajectory lengths (e.g., on the 100 ns) timescale. Replica exchange MD, however, results in quite efficient sampling due to the rapid increase in isomerization rate with temperature. The equilibrium distributions obtained from replica exchange MD have been compared with the results of experimental nuclear magnetic resonance observations. This comparison suggests that the sampling approach demonstrated here is a valuable tool that can be used in evaluating force fields for molecular simulation of lipids.

Saturation with cholesterol increases vertical order and smoothes the surface of the phosphatidylcholine bilayer: a molecular simulation study

Molecular dynamics (MD) simulations of a mono-cis-unsaturated 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayer and a POPC bilayer containing 50mol% cholesterol (POPC-Chol50) were carried out for 200ns to compare the spatial organizations of the pure POPC bilayer and the POPC bilayer saturated with Chol. The results presented here indicate that saturation with Chol significantly narrows the distribution of vertical positions of the center-of-mass of POPC molecules and POPC atoms in the bilayer. In the POPC-Chol50 bilayer, the same moieties of the lipid molecules are better aligned at a given bilayer depth, forming the following clearly separated membrane regions: the polar headgroup, the rigid core consisting of steroid rings and upper fragments of the acyl chains, and the fluid hydrocarbon core consisting of Chol chains and the lower fragments of POPC chains. The membrane surface of the POPC-Chol50 bilayer is smooth. The results have biological significance because the POPC-Chol50 bilayer models the bulk phospholipid portion of the fiber-cell membrane in the eye lens. It is hypothesized that in the eye lens cholesterol-induced smoothing of the membrane surface decreases light-scattering and helps to maintain lens transparency.

Thermal phase behavior of DMPG bilayers in aqueous dispersions as revealed by 2H- and 31P-NMR

The synthetic lipid 1,2-dimyristoyl-sn-3-phosphoglycerol (DMPG), when dispersed in water/NaCl exhibits a complex phase behavior caused by its almost unlimited swelling in excess water. Using deuterium ((2)H)- and phosphorus ((31)P)-NMR we have studied the molecular properties of DMPG/water/NaCl dispersions as a function of lipid and NaCl concentration. We have measured the order profile of the hydrophobic part of the lipid bilayer with deuterated DMPG while the orientation of the phosphoglycerol headgroup was deduced from the (31)P NMR chemical shielding anisotropy. At temperatures >30 °C we observe well-resolved (2)H- and (31)P NMR spectra not much different from other liquid crystalline bilayers. From the order profiles it is possible to deduce the average length of the flexible fatty acyl chain. Unusual spectra are obtained in the temperature interval of 20-25 °C, indicating one or several phase transitions. The most dramatic changes are seen at low lipid concentration and low ionic strength. Under these conditions and at 25 °C, the phosphoglycerol headgroup rotates into the hydrocarbon layer and the hydrocarbon chains show larger flexing motions than at higher temperatures. The orientation of the phosphoglycerol headgroup depends on the bilayer surface charge and correlates with the degree of dissociation of DMPG-Na(+). The larger the negative surface charge, the more the headgroup rotates toward the nonpolar region.

Conformational study of methylphosphocholine: a prototype for phospholipid headgroups in membranes

Phospholipid bilayers constitute the largest structural component of cell membranes, in which choline phospholipids are abundant. In this study, through a theoretical sampling on a methylphosphocholine (MePC) potential energy surface, a set of conformers was selected as a prototype for the membrane phospholipid head. We performed a detailed conformational study of such a prototype, both as an isolated moiety and in a solvated system. We used the polarizable continuum model (PCM) to account for solvation effects. We used a quantum-mechanical methodology based on density functional theory (DFT) and the 6-31G(d,p) basis set for the calculations. Through this methodology we were able to obtain a set of conformations that presented a mirror-image pattern, in good agreement with the experimental geometric values for the different phosphocholine derivatives. Potential curves for the main parameters of the dihedral space of MePC were obtained and are provided to guide future force-field parameterizations.

Update of the CHARMM all-atom additive force field for lipids: validation on six lipid types

A significant modification to the additive all-atom CHARMM lipid force field (FF) is developed and applied to phospholipid bilayers with both choline and ethanolamine containing head groups and with both saturated and unsaturated aliphatic chains. Motivated by the current CHARMM lipid FF (C27 and C27r) systematically yielding values of the surface area per lipid that are smaller than experimental estimates and gel-like structures of bilayers well above the gel transition temperature, selected torsional, Lennard-Jones and partial atomic charge parameters were modified by targeting both quantum mechanical (QM) and experimental data. QM calculations ranging from high-level ab initio calculations on small molecules to semiempirical QM studies on a 1,2-dipalmitoyl-sn-phosphatidylcholine (DPPC) bilayer in combination with experimental thermodynamic data were used as target data for parameter optimization. These changes were tested with simulations of pure bilayers at high hydration of the following six lipids: DPPC, 1,2-dimyristoyl-sn-phosphatidylcholine (DMPC), 1,2-dilauroyl-sn-phosphatidylcholine (DLPC), 1-palmitoyl-2-oleoyl-sn-phosphatidylcholine (POPC), 1,2-dioleoyl-sn-phosphatidylcholine (DOPC), and 1-palmitoyl-2-oleoyl-sn-phosphatidylethanolamine (POPE); simulations of a low hydration DOPC bilayer were also performed. Agreement with experimental surface area is on average within 2%, and the density profiles agree well with neutron and X-ray diffraction experiments. NMR deuterium order parameters (S(CD)) are well predicted with the new FF, including proper splitting of the S(CD) for the aliphatic carbon adjacent to the carbonyl for DPPC, POPE, and POPC bilayers. The area compressibility modulus and frequency dependence of (13)C NMR relaxation rates of DPPC and the water distribution of low hydration DOPC bilayers also agree well with experiment. Accordingly, the presented lipid FF, referred to as C36, allows for molecular dynamics simulations to be run in the tensionless ensemble (NPT), and is anticipated to be of utility for simulations of pure lipid systems as well as heterogeneous systems including membrane proteins.

Differential Effects of Cholesterol, Ergosterol and Lanosterol on a Dipalmitoyl Phosphatidylcholine Membrane: A Molecular Dynamics Simulation Study

Lipid raft/domain formation may arise as a result of the effects of specific sterols on the physical properties of membranes. Here, using molecular dynamics simulation, we examine the effects of three closely-related sterols, ergosterol, cholesterol, and lanosterol, at a biologically relevant concentration (40 mol %) on the structural properties of a model dipalmitoyl phosphatidylcholine (DPPC) membrane at 309 and 323 K. All three sterols are found to order the DPPC acyl tails and condense the membrane relative to the DPPC liquid-phase membrane, but each one does this to a significantly different degree. The smooth alpha-face of ergosterol, together with the presence of tail unsaturation in this sterol, leads to closer interaction of ergosterol with the lipids and closer packing of the lipids with each other, so ergosterol has a higher condensing effect on the membrane, as reflected by the area per lipid. Moreover, ergosterol induces a higher proportion of trans lipid conformers, a thicker membrane, and higher lipid order parameters and is aligned more closely with the membrane normal. Ergosterol also positions itself closer to the bilayer/water interface. In contrast, the rough alpha-face of lanosterol leads to a less close interaction of the steroid ring system with the phospholipid acyl chains, and so lanosterol orders, straightens, and packs the lipid acyl chains less well and is less closely aligned with the membrane normal. Furthermore, lanosterol lies closer to the relatively disordered membrane center than do the other sterols. The behavior of cholesterol in all the above respects is intermediate between that of lanosterol and ergosterol. The findings here may explain why ergosterol is the most efficient of the three sterols at promoting the liquid-ordered phase and lipid domain formation and may also furnish part of the explanation as to why cholesterol is evolutionarily preferred over lanosterol in higher-vertebrate plasma membranes.

Lipid conformation in crystalline bilayers and in crystals of transmembrane proteins

Dihedral torsion angles evaluated for the phospholipid molecules resolved in the X-ray structures of transmembrane proteins in crystals are compared with those of phospholipids in bilayer crystals, and with the phospholipid conformations in fluid membranes. Conformations of the lipid glycerol backbone in protein crystals are not restricted to the gauche C1-C2 rotamers found invariably in phospholipid bilayer crystals. Lipid headgroup conformations in protein crystals also do not conform solely to the bent-down conformation, with gauche-gauche configuration of the phospho-diester, that is characteristic of phospholipid bilayer membranes. This suggests that the lipids that are resolved in crystals of membrane proteins are not representative of the entire lipid-protein interface. Much of the chain configurational disorder of the membrane-bound lipids in crystals arises from energetically disallowed skew conformations. This indicates a configurational heterogeneity in the lipids at a single binding site: eclipsed conformations occur also in some glycerol backbone torsion angles and C-C torsion angles in the lipid headgroups. Stereochemical violations in the protein-bound lipids are evidenced by one-third of the ester carboxyl groups in non-planar configurations, and certain of the carboxyls in the cis configuration. Some of the lipid structures in protein crystals have the incorrect enantiomeric configuration of the glycerol backbone, and many of the branched methyl groups in structures of the phytanyl chains associated with bacteriorhodopsin crystals are in the incorrect S-configuration.

Empirical force fields for biological macromolecules: overview and issues

Empirical force field-based studies of biological macromolecules are becoming a common tool for investigating their structure-activity relationships at an atomic level of detail. Such studies facilitate interpretation of experimental data and allow for information not readily accessible to experimental methods to be obtained. A large part of the success of empirical force field-based methods is the quality of the force fields combined with the algorithmic advances that allow for more accurate reproduction of experimental observables. Presented is an overview of the issues associated with the development and application of empirical force fields to biomolecular systems. This is followed by a summary of the force fields commonly applied to the different classes of biomolecules; proteins, nucleic acids, lipids, and carbohydrates. In addition, issues associated with computational studies on "heterogeneous" biomolecular systems and the transferability of force fields to a wide range of organic molecules of pharmacological interest are discussed.

Sphingomyelin and cholesterol: from membrane biophysics and rafts to potential medical applications

The preferential sphingomyelin-cholesterol interaction which results from the structure and the molecular properties of these two lipids seems to be the physicochemical basis for the formation and maintenance of cholesterol/sphingolipid-enriched nano- and micro-domains (referred to as membrane "rafts") in the plane of plasma and other organelle (i.e., Golgi) membranes. This claim is supported by much experimental evidence and also by theoretical considerations. However, although there is a large volume of information about these rafts regarding their lipid and protein composition, their size, and their dynamics, there is still much to be clarified on these issues, as well as on how rafts are formed and maintained. It is well accepted now that the lipid phase of the rafts is the liquid ordered (LO) phase. However, other (non-raft) parts of the membrane may also be in the LO phase. There are indications that the raft LO phase domains are more tightly packed than the non-raft LO phase, possibly due to intermolecular hydrogen bonding involving sphingolipid and cholesterol. This also explains why the former are detergent-resistant membranes (DRM), while the non-raft LO phase domains are detergent-soluble (sensitive) membranes (DSM). Recent findings suggest that protein-protein interactions such as cross-linking can be controlled by protein distribution between raft and non-raft domains, and, as well, these interactions affect raft size distribution. The cholesterol/sphingomyelin-enriched rafts seem to be involved in many biological processes, mediated by various receptors, as exemplified by various lipidated glycosylphosphatidylinositol (GPI)- and acyl chain-anchored proteins that reside in the rafts. The rafts serve as signaling platforms in the cell. Various pathogens (viruses and toxins) utilize the raft domains on the host cell membrane as a port of entry, site of assembly (viruses), and port of exit (viral budding). Existence and maintenance of cholesterol-sphingomyelin rafts are dependent on the level of membrane cholesterol and sphingomyelin. This explains why reduction of cholesterol level--either through reverse cholesterol transport, using cholesterol acceptors such as beta-cyclodextrin, or through cholesterol biosynthesis inhibition using statins--interferes with many processes which involve rafts and can be applied to treating raft-related infections and diseases. Detailed elucidation of raft structure and function will improve understanding of biological membrane composition-structure-function relationships and also may serve as a new avenue for the development of novel treatments for major diseases, including viral infections, neurodegenerative diseases (Alzheimer's), atherosclerosis, and tumors.

Orientational changes in dipalmitoyl phosphatidyl glycerol Langmuir monolayers

Dipalmitoyl phosphatidyl glycerol (DPPG) as Langmuir monolayers at the air/water interface was investigated by means of surface pressure measurements in addition to Brewster angle microscopy (BAM) during film compression/expansion. A characteristic phase transition region appeared in the course of surface pressure-area (pi-A) isotherms for monolayers spread on alkaline water or buffer subphase, while on neutral or acidic water the plateau region was absent. This phase transition region was attributed to the ionization of DPPG monolayer. It has been postulated that the ionization of the phosphatidyl glycerol group leads to its increased solvation, which probably provokes both a change in the orientation of the polar group and its deeper penetration into bulk phase. Film compression along the transition region provokes the dehydration of polar groups and subsequent change of their conformation, thus causing the DPPG molecules to emerge up to the interface. Quantitative Brewster angle microscopy (BAM) measurements revealed that along the liquid-expanded to liquid-condensed phase transition the thickness of the ionized DPPG monolayer increases by 4.2 A as a result of the conformational changes of the ionized polar groups, which tend to emerge from the bulk subphase up to the surface.

Both isoforms of mammalian phosphatidylinositol transfer protein are capable of binding and transporting sphingomyelin

The structurally related mammalian alpha and beta isoforms of phosphatidylinositol (PtdIns) transfer protein (PITP) bind reversibly a single phospholipid molecule, preferably PtdIns or phosphatidylcholine (PtdCho), and transport that lipid between membrane surfaces. PITPbeta, but not PITPalpha, is reported extensively in the scientific literature to exhibit the additional capacity to bind and transport sphingomyelin (CerPCho). We undertook a detailed investigation of the lipid binding and transfer specificity of the soluble mammalian PITP isoforms. We employed a variety of donor and acceptor membrane lipid compositions to determine the sensitivity of recombinant rat PITPalpha and PITPbeta isoforms toward PtdIns, PtdCho, CerPCho, and phosphatidate (PtdOH). Results indicated often striking differences in protein-phospholipid and protein-membrane interactions. We demonstrated unequivocally that both isoforms were capable of binding and transferring CerPCho; we confirmed that the beta isoform was the more active. The order of transfer specific activity was similar for both isoforms: PtdIns>PtdCho>CerPCho>>PtdOH. Independently, we verified the binding of CerPCho to both isoforms by showing an increase in holoprotein isoelectric point following the exchange of protein-bound phosphatidylglycerol for membrane-associated CerPCho. We conclude that PITPalpha and PITPbeta are able to bind and transport glycero- and sphingophospholipids.

Molecular dynamics generation of nonarbitrary membrane models reveals lipid orientational correlations

This report addresses the following problems associated with the generation of computer models of phospholipid bilayer membranes using molecular dynamics simulations: arbitrary initial structures and short equilibration periods, an Ewald-induced strong coupling of phospholipids, uncertainty regarding which value should be used for surface tension to alleviate the problem of the small size of the membrane, and simultaneous realization of both order parameters and the surface area. We generated a computer model of the liquid-crystalline L-alpha-dimyristoylphosphatidylcholine (DMPC) bilayer, starting from a configuration based on a crystal structure (rather than from an arbitrary structure). To break the crystalline structure, a 20-ps high-temperature pulse of 510 K (but not 450 or 480 K) was effective. The system finally obtained is an all-atom model, with Ewald summation to evaluate Coulombic interactions and a constant surface tension of 35 dynes/cm/water-membrane interface, equilibrated for 12 ns (over 50 ns total calculation time), which reproduces all of the experimentally observed parameters examined in this work. Furthermore, this model shows the presence of significant orientational correlations between neighboring alkyl chains and between shoulder vectors (which show the orientations of the lipids about their long axes) of neighboring DMPCs.

Molecular dynamics and (2)H-NMR study of the influence of an amphiphilic peptide on membrane order and dynamics

A molecular dynamics simulation of a fully hydrated model membrane consisting of 12 molecules of 1, 2-dimyristoyl-sn-glycero-3-phosphocholine, one amphiphilic peptide with the sequence acetyl-Lys-Lys-Gly-Leu(16)-Lys-Lys-Ala-amide, and 593 water molecules was performed for 1.06 ns (Belohorcova, K., J. H. Davis, T. B. Woolf, and B. Roux. 1997. Biophys. J. 73:3039-3055). The analysis presented here is primarily focused on the phospholipid component and the results are compared with experimental (2)H-NMR studies of the lipid component of mixtures of the same peptide and lipid at a molar ratio of 1:32, and with earlier studies of closely related peptide/lipid mixtures. The phospholipid chain and headgroup isomer populations and isomerization rates compare favorably with previous simulations and experimental measurements. Of particular interest is the effect of the peptide on the phospholipid headgroup and hydrocarbon chain orientational order calculated from the simulation, which also agree well with experimental measurements performed on this and closely related systems. Comparison of the experimental results with the simulations not only shows that there is significant agreement between the two methods, but also provides new insight into the effect of the peptide on the lipid dynamics. In particular, these results confirm that a membrane spanning peptide has little effect on lipid chain order, and bilayer thickness if its hydrophobic length closely matches the lipid hydrocarbon thickness. In addition, we find that the peptide can have a strong ordering effect if it is longer than the lipid hydrophobic thickness.

Short-chain phospholipids as detergents

The physico-chemical properties of short-chain phosphatidylcholine are reviewed to the extent that its biological activity as a mild detergent can be rationalized. Long-chain diacylphosphatidylcholines are typical membrane phospholipids that form preferentially smectic lamellar phases (bilayers) when dispersed in water. In contrast, the preferred phase of the short-chain analogues dispersed in excess water is the micellar phase. The preferred conformation and the dynamics of short-chain phosphatidylcholines in the monomeric and micellar state present in H(2)O are discussed. The motionally averaged conformation of short-chain phosphatidylcholines is then compared to the single-crystal structures of membrane lipids. The main conclusion emerging is that in terms of preferred conformation and motional averaging short-chain phosphatidylcholines closely resemble their long-chain analogues. The dispersing power of short-chain phospholipids is emphasized in the second part of the review. Evidence is presented to show that this class of compounds is superior to most other detergents used in the solubilization of membrane proteins and the reconstitution of the solubilized proteins to artificial membrane systems (proteoliposomes). The prominent feature of the solubilization/reconstitution of integral membrane proteins by short-chain PC is the retention of the native protein structure and hence the protein function. Due to their special detergent-like properties, short-chain PC lend themselves very well not only to membrane solubilization but also to the purification of integral membrane proteins. The retention of the native protein structure in the solubilized state, i.e. in mixed micelles consisting of the integral membrane protein, intrinsic membrane lipids and short-chain PC, is rationalized. It is hypothesized that short-chain PC interacts primarily with the lipid bilayer of a membrane and very little if at all with the membrane proteins. In this way, the membrane protein remains associated with its preferred intrinsic membrane lipids and retains its native structure and its function.

Conformational studies of sphingolipids by NMR spectroscopy. I. Dihydrosphingomyelin

The conformational features of dihydrosphingomyelin (DHSM), the major phospholipid of human lens membranes, were investigated by 1H and 31P nuclear magnetic resonance spectroscopy. Several postulates emerge from the observed trends: (a) in partially hydrated samples of DHSM in CDCl3 above 13 mM, at which lipid-lipid interactions prevail, the amide proton is mostly involved in intermolecular H-bonds that link neighboring phospholipids through bridging water molecules. In the absence of water, the NH group is involved in an intramolecular H-bond that restricts the mobility of the phosphate group. (b) In the monomeric form of the lipid molecule, the amide proton of the major conformer is bound intramolecularly with one of the anionic and/or ester oxygens of the phosphate group. A minor conformer may also be present in which the NH proton participates in an intramolecular H-bond linking to the OH group of the sphingoid base. (c) Complete hydration leads to an extension of the head group as water molecules bind to the phosphate and NH groups via H-bonds, thus disrupting the intramolecular H-bonds prevalent at low concentrations.

1H and (13)C NMR of multilamellar dispersions of polyunsaturated (22:6) phospholipids

The polyunsaturated fatty acid docosahexaenoic acid (DHA) makes up approximately 50% of the lipid chains in the retinal rod outer segment disk membranes and a large fraction of the lipid chains in the membranes of neuronal tissues. There is an extensive literature concerned with the dietary requirements for essential fatty acids and the importance of DHA to human health, but relatively little research has been done on the physical properties of this important molecule. Using (1)H and (13)C MAS NMR measurements of dispersions of 1-palmitoyl-2-docosahexaenoyl-phosphatidylcholine in excess phosphate buffer, we have unambiguously assigned most of the resonances in both the (1)H and (13)C NMR spectra. We were able to use cross-polarization spectroscopy to follow the transfer of polarization from specific (1)H nuclei not only to their directly bonded (13)C but also to those (13)C that are in close proximity, even though they are not directly bonded. Cross-peaks in two-dimensional cross-polarization spectra revealed a close association between the choline headgroup and at least part of the DHA chain but not with the palmitate chain. Finally, we examined the dynamics of the different parts of this lipid molecule, using rotating frame spin-lattice relaxation measurements, and found that methylene groups of both chains experience important motions with correlation times in the 10-micros range, with those for the palmitate chain being approximately 50% longer than those of the DHA chain. The choline headgroup and the chain terminal groups have significantly shorter correlation times, and that part of the dipolar interaction that is fluctuating at these correlation times is significantly smaller for these groups than it is for the palmitate and DHA chain methylenes.

Molecular dynamics simulations of biomembrane models

A molecular force field dedicated to molecular dynamics simulation of biomembranes was developed. It was parameterized on model compounds related to phospholipids and was able to reproduce at the same time structures, energies, and vibrational spectra. Cross terms in the potential energy function were introduced by solving the redundancy problem among internal coordinates. This force field was used in the 400-ps molecular dynamics simulation of a hydrated bilayer in the gel and liquid crystal phases. The conformational properties of the polar head groups were in particular agreement with the experimental observations using Raman scattering.

Molecular dynamics simulation of conformational flexibility of alamethicin fragments in aqueous and membranous environment

We present here results on molecular dynamics (MD) simulation on two fragments of channel forming antibiotic peptide Alamethicin, containing isoamino butyric acid (Aib). Simulations are carried out in aqueous and membranous environment in a bilayer of 39 molecules of Dimyristoyl phosphatidyl choline (DMPC). The peptides Boc-Pro-Aib-Ala-Aib-OBzl (Alam 1) and Boc-Leu-Aib-Pro-OBzl (Alam 2) were simulated from their crystallographic coordinates. The bilayers were built from two different conformations (A and B) of DMPC reported in crystal data. The P-N dipoles were arranged hexagonally with surface area per lipid molecule 66.5 A degrees 2 and P-P separation across the bilayer 34 A degrees. They were hydrated by 28.6 and 25.5 water molecules per DMPC molecule. Simulations are done using AMBER 4.0 package in constant number volume temperature (NVT) condition for 100 pico seconds (ps) in aqueous environment and 250 ps of equilibrated bilayer. Geometric parameters of lipids as: bilayer thickness, order parameter of the chains, transfraction of chain torsional angles were monitored. We also monitored geometric parameters of the peptides as backbone torsional angles, distances amongst C alpha atoms, angles between C alpha atoms, movement of center of gravity (CG) along and perpendicular to bilayer normal. We find that membrane bilayer is slightly disturbed due to the presence of peptides. In case of alam 2 in water angles phi 1 and phi 3 showed larger variation in water compared to same in the bilayer. The peptide conformation is more stable in DMPC bilayer. However the peptides showed movement along and perpendicular to bilayer normal. This we believe is due to hydrophobic nature of these peptides.

260 ps molecular dynamics simulation of substance P with hydrated dimyristoyl phosphatidyl choline bilayer

We present here results on 260 pico seconds (ps) molecular dynamics (MD) simulation of substance P (SP) in hydrated bilayer of dimyristoyl phosphatidyl choline (DMPC) (39 molecules of DMPC with 776 water molecules). 260 ps MD simulation has been carried out in 0.001 ps time interval with united atom force field, using AMBER 4.0 package. Non bonded pair list was updated every 20 cycles using 12.5 Angstrom cut off distance. Analysis of MD data is done using our package ANALMD. The obtained models are presented using graphics package RASMOL. All simulations, analysis of MD data and graphics is done on INDIGO-2, R-4400 extreme graphics work station. Our results show no systematic change in order parameter, but reduction in transfraction of the chain torsional angles, compared to our earlier results on MD simulation on hydrated DMPC bilayer without SP. C-terminal and central peptide residues adopt partial helical conformation. Helix type as classified on the basis of H-bonds is between alpha and 3(10). The peptide backbone shows flexibility during heating runs. Later, it is stabilized and there was not much change in the spatial position of the backbone. Lipid matrix serves the role of immobilization of the peptide backbone in a preferred conformation.

Phases and phase transitions of the sphingolipids

LIPIDAT is a computerized database providing access to the wealth of information scattered throughout the literature concerning synthetic and biologically derived polar lipid polymorphic and mesomorphic phase behavior. Herein, we present a review of the LIPIDAT data subset referring to sphingolipids together with an analysis of these data. It includes data collected over a 40-year period and consists of 867 records obtained from 112 articles in 25 different journals. An analysis of these data has allowed us to identify trends in hydrated sphingolipid phase behavior reflecting differences in fatty acyl chain length, saturation and hydroxylation, head group type, and sphingoid base identity. Information on the mesomorphism of biologically-derived and dry sphingolipids is also presented. This review includes 161 references.