Molecular dynamics studies of the molecular structure and interactions of cholesterol superlattices and random domains in an unsaturated phosphatidylcholine bilayer membrane

The effect of the molecular organization of lipid components on the properties of the bilayer membrane has been a topic of increasing interest. Several experimental and theoretical studies have suggested that cholesterol is not randomly distributed in the fluid-state lipid bilayer but forms nanoscale domains. Several cholesterol-enriched nanodomain structures have been proposed, including rafts, regular or maze arrays, complexes, and superlattices. At present, the molecular mechanisms by which lipid composition influences the formation and stability of lipid nanodomains remain unclear. In this study, we have used molecular dynamics (MD) simulations to investigate the effects of the molecular organization of cholesterol--superlattice versus random--on the structure of and interactions between lipids and water in lipid bilayers of cholesterol and 1-palmitoyl-2-oleoylphosphatidylcholine (cholesterol/POPC) at a fixed cholesterol mole fraction of 0.40. On the basis of four independent replicates of 200-ns MD simulations for a superlattice or random bilayer, statistically significant differences were observed in the lipid structural parameters, area per lipid, density profile, and acyl chain order profile, as well as the hydrogen bonding between various pairs (POPC and water, cholesterol and water, and POPC and cholesterol). The time evolution of the radial distribution of the cholesterol hydroxy oxygen suggests that the lateral distribution of cholesterol in the superlattice bilayer is more stable than that in the random bilayer. Furthermore, the results indicate that a relatively long simulation time, more than 100 ns, is required for these two-component bilayers to reach equilibrium and that this time is influenced by the initial lateral distribution of lipid components.     

Computational analysis of local membrane properties

In the field of biomolecular simulations, dynamics of phospholipid membranes is of special interest. A number of proteins, including channels, transporters, receptors and short peptides are embedded in lipid bilayers and tightly interact with phospholipids. While the experimental measurements report on the spatial and/or temporal average membrane properties, simulation results are not restricted to the average properties. In the current study, we present a collection of methods for an efficient local membrane property calculation, comprising bilayer thickness, area per lipid, deuterium order parameters, Gaussian and mean curvature. The local membrane property calculation allows for a direct mapping of the membrane features, which subsequently can be used for further analysis and visualization of the processes of interest. The main features of the described methods are highlighted in a number of membrane systems, namely: a pure dimyristoyl-phosphatidyl-choline (DMPC) bilayer, a fusion peptide interacting with a membrane, voltage-dependent anion channel protein embedded in a DMPC bilayer, cholesterol enriched bilayer and a coarse grained simulation of a curved palmitoyl-oleoyl-phosphatidyl-choline lipid membrane. The local membrane property analysis proves to provide an intuitive and detailed view on the observables that are otherwise interpreted as averaged bilayer properties.     

Accelerated molecular dynamics simulations of protein folding

Folding of four fast‐folding proteins, including chignolin, Trp‐cage, villin headpiece and WW domain, was simulated via accelerated molecular dynamics (aMD). In comparison with hundred‐of‐microsecond timescale conventional molecular dynamics (cMD) simulations performed on the Anton supercomputer, aMD captured complete folding of the four proteins in significantly shorter simulation time. The folded protein conformations were found within 0.2–2.1 Å of the native NMR or X‐ray crystal structures. Free energy profiles calculated through improved reweighting of the aMD simulations using cumulant expansion to the second‐order are in good agreement with those obtained from cMD simulations. This allows us to identify distinct conformational states (e.g., unfolded and intermediate) other than the native structure and the protein folding energy barriers. Detailed analysis of protein secondary structures and local key residue interactions provided important insights into the protein folding pathways. Furthermore, the selections of force fields and aMD simulation parameters are discussed in detail. Our work shows usefulness and accuracy of aMD in studying protein folding, providing basic references in using aMD in future protein‐folding studies. © 2015 Wiley Periodicals, Inc.     

Structural, dynamic, and electrostatic properties of fully hydrated DMPC bilayers from molecular dynamics simulations accelerated with graphical processing units (GPUs)

We present results of molecular dynamics simulations of fully hydrated DMPC bilayers performed on graphics processing units (GPUs) using current state-of-the-art non-polarizable force fields and a local GPU-enabled molecular dynamics code named FEN ZI. We treat the conditionally convergent electrostatic interaction energy exactly using the particle mesh Ewald method (PME) for solution of Poisson's Equation for the electrostatic potential under periodic boundary conditions. We discuss elements of our implementation of the PME algorithm on GPUs as well as pertinent performance issues. We proceed to show results of simulations of extended lipid bilayer systems using our program, FEN ZI. We performed simulations of DMPC bilayer systems consisting of 17,004, 68,484, and 273,936 atoms in explicit solvent. We present bilayer structural properties (atomic number densities, electron density profiles), deuterium order parameters (S(CD)), electrostatic properties (dipole potential, water dipole moments), and orientational properties of water. Predicted properties demonstrate excellent agreement with experiment and previous all-atom molecular dynamics simulations. We observe no statistically significant differences in calculated structural or electrostatic properties for different system sizes, suggesting the small bilayer simulations (less than 100 lipid molecules) provide equivalent representation of structural and electrostatic properties associated with significantly larger systems (over 1000 lipid molecules). We stress that the three system size representations will have differences in other properties such as surface capillary wave dynamics or surface tension related effects that are not probed in the current study. The latter properties are inherently dependent on system size. This contribution suggests the suitability of applying emerging GPU technologies to studies of an important class of biological environments, that of lipid bilayers and their associated integral membrane proteins. We envision that this technology will push the boundaries of fully atomic-resolution modeling of these biological systems, thus enabling unprecedented exploration of meso-scale phenomena (mechanisms, kinetics, energetics) with atomic detail at commodity hardware prices.     

Effects of lipid chain length on molecular interactions between paclitaxel and phospholipid within model biomembranes

Molecular interactions between an anticancer drug, paclitaxel, and phosphatidylcholine (PC) of various chain lengths were investigated in the present work by the Langmuir film balance technique and differential scanning calorimetry (DSC). Both the lipid monolayer at the air-water interface and lipid bilayer vesicles (liposomes) were employed as model biological cell membranes. Measurement and analysis of the surface pressure versus molecular area curves of the mixed monolayers of phospholipids and paclitaxel under various molar ratio showed that phospholipids and paclitaxel formed a nonideal miscible system at the interface. Paclitaxel exerted an area-condensing effect on the lipid monolayer at small molecular surface areas and an area-expanding effect at large molecular areas, which could be explained by the intermolecular forces and geometric accommodation between the two components. Paclitaxel and phospholipids could form thermodynamically stable monolayer systems: the stability increased with the chain length in the order DMPC (C14:0)>DPPC (C16:0)>DSPC (C18:0). Investigation of paclitaxel penetration into the pure lipid monolayer showed that DMPC had a higher ability to incorporate paclitaxel and the critical surface pressure for paclitaxel penetration also increased with the chain length in the order DMPC>DPPC>DSPC. A similar trend was testified by DSC studies on vesicles of the mixed paclitaxel/phospholipids bilayer. Paclitaxel showed the greatest interaction with DMPC while little interaction could be measured in the paclitaxel/DSPC liposomes. Paclitaxel caused broadening of the main phase transition without significant change at the peak melting temperature of the phospholipid bilayers, which demonstrated that paclitaxel was localized in the outer hydrophobic cooperative zone of the bilayer. The interaction between paclitaxel and phospholipid was nonspecific and the dominant factor in this interaction was the van der Waals force or hydrophobic force. As the result of the lower net van der Waals interaction between hydrocarbon chains for the shorter acyl chains, paclitaxel interacted more readily with phospholipids of shorter chain length, which also increased the bilayer intermolecular spacing.     

Molecular characterization of gel and liquid-crystalline structures of fully hydrated POPC and POPE bilayers

Molecular dynamics simulations were used for a comprehensive study of the structural properties of monounsaturated POPC and POPE bilayers in the gel and liquid-crystalline state at a number of temperatures, ranging from 250 to 330 K. Though the chemical structures of POPC and POPE are largely similar (choline versus ethanolamine headgroup), their transformation processes from a gel to a liquid-crystalline state are contrasting. In the similarities, the lipid tails for both systems are tilted below the phase transition and become more random above the phase transition temperature. The average area per lipid and bilayer thickness were found less sensitive to phase transition changes as the unsaturated tails are able to buffer reordering of the bilayer structure, as observed from hysteresis loops in annealing simulations. For POPC, changes in the structural properties such as the lipid tail order parameter, hydrocarbon trans-gauche isomerization, lipid tail tilt-angle, and level of interdigitation identified a phase transition at about 270 K. For POPE, three temperature ranges were identified, in which the lower one (270-280 K) was associated with a pre-transition state and the higher (290-300 K) with the post-transition state. In the pre-transition state, there was a significant increase in the number of gauche arrangements formed along the lipid tails. Near the main transition (280-290 K), there was a lowering of the lipid order parameters and a disappearance of the tilted lipid arrangement. In the post-transition state, the carbon atoms along the lipid tails became less hindered as their density profiles showed uniform distributions. This study also demonstrates that atomistic simulations of current lipid force fields are capable of capturing the phase transition behavior of lipid bilayers, providing a rich set of molecular and structural information at and near the main transition state.     

Molecular dynamics simulations of ternary membrane mixture: phosphatidylcholine, phosphatidic acid, and cholesterol

A ternary mixture of 1-palmitoyl-2-oleoyl phosphatidylcholine (POPC), 1-palmitoyl-2-oleoyl phosphatidic acid (POPA), and cholesterol (CHOL) works effectively for a functional conformation of nicotinic acetylcholine receptor (nAChR) that can undergo agonist-induced conformation changes, but POPC alone can stabilize only a desensitized state of nAChR. To gain insights into the lipid mixture that has strong impact to nAChR functions, we performed more than 50 ns all atom molecular dynamic (MD) simulations at 303 K on a fully hydrated bilayer consisting of 240 POPC, 80 POPA, and 80 CHOL (3:1:1). The MD simulation revealed various interactions between different types of molecular pairs that ultimately regulated lipid organization. The heterogeneous interactions among three different constituents resulted in a broad spectrum of lipid properties, including extensive distributions of average area per lipid and varied lipid ordering as a function of lipid closeness to CHOL. Higher percentage of POPA than POPC had close association with CHOL, which coincided with relatively higher ordering of POPA molecules in their acyl chains near lipid head groups. Lower fraction of gauche dihedrals was also found in the same region of POPA. Although the CHOL molecules had the effects on the enhancement of surrounding lipid order, relatively low lipid order parameters and high fraction of gauche bonds were observed in the ternary mixture. Collectively, these results suggest that the dynamical structure of the ternary system could be determinant for a functional nAChR.     

Molecular dynamics simulations of rhodopsin in different one-component lipid bilayers

Four 20 ns molecular dynamic simulations of rhodopsin embedded in different one-component lipid bilayers have been carried out to ascertain the importance of membrane lipids on the protein structure. Specifically, dimyristoyl phosphatidylcholine (DMPC), dipalmitoyl phosphatidylcholine (DPPC), palmitoyl oleoyl phosphatidylcholine (POPC), and palmitoyl linoleyl phosphatidylcholine (PLPC) lipid bilayers have been considered for the present work. The results reported here provide information on the hydrophobic matching between the protein and the bilayer and about the differential effects of the protein on the thickness of the different membranes. Furthermore, a careful analysis of the individual protein-lipid interactions permits the identification of residues that exhibit permanent interactions with atoms of the lipid environment that may putatively act as hooks of the protein to the membrane. The analysis of the trajectories also provides information about the effect of the bilayer on the protein structure, including secondary structural elements, salt bridges, and rigid-body motions.     

Electrochemical and PM-IRRAS studies of the effect of cholesterol on the properties of the headgroup region of a DMPC bilayer supported at a Au(111) electrode

Polarization modulation infrared reflection absorption spectroscopy (PM-IRRAS) was employed to investigate the interaction of cholesterol with the headgroups of dimyristoylphosphatidycholine (DMPC) molecules under a static electric field. DMPC/cholesterol (7:3 molar ratio) mixtures form a bilayer on a Au(111) electrode surface by fusion and spreading of small unilamellar vesicles. PM-IRRAS experiments provided detailed information concerning the conformation and hydration of headgroups of DMPC bilayers in the presence and absence of 30% cholesterol. The presence of 30% cholesterol increases the space between the headgroups of DMPC molecules and hence increases the hydration of the DMPC/cholesterol mixed bilayer. The conformational state of the headgroups of DMPC molecules in the mixed bilayer is also significantly changed. The phosphate group is closer to the surface compared with the pure DMPC bilayer. The conformation of the -O-C-C-N moiety changes from gauche to trans in the presence of cholesterol.     

Combination of COSMOmic and molecular dynamics simulations for the calculation of membrane–water partition coefficients

The importance of membrane‐water partition coefficients led to the recent extension of the conductor‐like screening model for realistic solvation (COSMO‐RS) to micelles and biomembranes termed COSMOmic. Compared to COSMO‐RS, this new approach needs structural information to account for the anisotropy of colloidal systems. This information can be obtained from molecular dynamics (MD) simulations. In this work, we show that this combination of molecular methods can efficiently be used to predict partition coefficients with good agreement to experimental data and enables screening studies. However, there is a discrepancy between the amount of data generated by MD simulations and the structural information needed for COSMOmic. Therefore, a new scheme is presented to extract data from MD trajectories for COSMOmic calculations. In particular, we show how to calculate the system structure from MD, the influence of lipid conformers, the relation to the COSMOmic layer size, and the water/lipid ratio impact. For a 1,2‐dimyristoyl‐sn‐glycero‐3‐phosphocholine (DMPC) bilayer, 66 partition coefficients for various solutes were calculated. Further, 52 partition coefficients for a 1‐palmitoyl‐2‐oleoyl‐sn‐glycero‐3‐phosphocholine (POPC) bilayer system were calculated. All these calculations were compared to experimental data. © 2013 Wiley Periodicals, Inc.     

Systematic implicit solvent coarse graining of dimyristoylphosphatidylcholine lipids

We have used systematic structure‐based coarse graining to derive effective site–site potentials for a 10‐site coarse‐grained dimyristoylphosphatidylcholine (DMPC) lipid model and investigated their state point dependence. The potentials provide for the coarse‐grained model the same site–site radial distribution functions, bond and angle distributions as those computed in atomistic simulations carried out at four different lipid–water molar ratios. It was shown that there is a non‐negligible dependence of the effective potentials on the concentration at which they were generated, which is also manifested in the properties of the lipid bilayers simulated using these potentials. Thus, effective potentials computed at low lipid concentration favor to more condensed and ordered structure of the bilayer with lower average area per lipid, while potentials obtained at higher lipid concentrations provide more fluid‐like structure. The best agreement with the reference data and experiment was achieved using the set of potentials derived from atomistic simulations at 1:30 lipid:water molar ratio providing fully saturated hydration of DMPC lipids. Despite theoretical limitations of pairwise coarse‐grained potentials expressed in their state point dependence, all the resulting potentials provide a stable bilayer structure with correct partitioning of different lipid groups across the bilayer as well as acceptable values of the average lipid area, compressibility and orientational ordering. In addition to bilayer simulations, the model has proven its robustness in modeling of self‐aggregation of lipids from randomly dispersed solution to ordered bilayer structures, bicelles, and vesicles. © 2014 Wiley Periodicals, Inc.     

Perturbation of phospholipid bilayer properties by ethanol at a high concentration

Atomistic molecular dynamics (MD) simulations have been carried out at 30 degrees C on a fully hydrated liquid crystalline lamellar phase of dimyrystoylphosphatidylcholine (DMPC) lipid bilayer with embedded ethanol molecules at 1:1 composition, as well as on the pure bilayer phase. The ethanol molecules are found to exhibit a preference to occupy regions near the upper part of the lipid acyl chains and the phosphocholine headgroups. The calculations revealed that the phosphocholine headgroup dipoles (P- --> N+) of the lipids prefer to orient more toward the aqueous layer in the presence of ethanol. It is noticed that the ethanol molecules modify the dynamic properties of both lipids as well as the water molecules in the hydration layer of the lipid headgroups. Both the in-plane "rattling" and out-of-plane "protrusion" motions of the lipids have been found to increase in the presence of ethanol. Most importantly, it is observed that the water molecules within the hydration layer of the lipid headgroups exhibit faster translational and rotational motions in the presence of ethanol. This arises due to faster dynamics of hydrogen bonds between lipid headgroups and water in the presence of ethanol.     

Surface planar bilayers of phospholipids used in protein membrane reconstitution: an atomic force microscopy study

In this work, using atomic force microscopy (AFM), we have studied the influence of the temperature on the properties of the surface planar bilayers (SPBs) formed with: (i) the total lipid extract of Escherichia coli; (ii) 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 1,2-dimyristoyl-sn-glycero-3-phosphoethanolamine (DMPC) (1:1, mol/mol); and, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanol-amine (POPE) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol (POPG) (3:1, mol/mol). According to the height profile analysis we performed, the height of the SPBs of DMPC:POPC were temperature dependent. Separated domains were observed in the SPBs of the POPE:POPG mixture and the E. coli lipid extract. The implication of those domains for the correct insertion of membrane proteins into proteoliposomes is discussed.     

Comprehensive Characterization of the Self-Folding Cavitand Dynamics

The conformational equilibria and guest exchange process of a resorcin[4]arene derived self-folding cavitand receptor have been characterized in detail by molecular dynamics simulations (MD) and ¹H EXSY NMR experiments. A multi-timescale strategy for exploring the fluxional behaviour of this system has been constructed, exploiting conventional MD and accelerated MD (aMD) techniques. The use of aMD allows the reconstruction of the folding/unfolding process of the receptor by sampling high-energy barrier processes unattainable by conventional MD simulations. We obtained MD trajectories sampling events occurring at different timescales from ns to s: 1) rearrangement of the directional hydrogen bond seam stabilizing the receptor, 2) folding/unfolding of the structure transiting partially open intermediates, and 3) guest departure from different folding stages. Most remarkably, reweighing of the biased aMD simulations provided kinetic barriers that are in very good agreement with those determined experimentally by ¹H NMR. These results constitute the first comprehensive characterization of the complex dynamic features of cavitand receptors. Our approach emerges as a valuable rational design tool for synthetic host-guest systems     

Molecular dynamics simulations of the Ca2+-pump: a structural analysis

We report large scale molecular dynamics computer simulations, ～100 ns, of the ion pump Ca(2+)-ATPase immersed in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer. The structure simulated here, E1, one of the several conformations resolved using X-ray diffraction techniques, hosts two Ca(2+)-ions in the hydrophobic domain. Our results indicate that protonated residues lead to stronger ion-residue interactions, supporting previous conclusions regarding the sensitivity of the Ca(2+) behaviour to the protonated state of the amino acid binding sites. We also investigate how the protein perturbs the bilayer structure. We show that the POPC bilayer is ～12% thinner than the pure bilayer, near the protein surface. This perturbation decays exponentially with the distance from the protein with a characteristic decay length of 0.8 nm. We find that the projected area per lipid also decreases near the protein. Using an analytical model we show that this change in the area is only apparent and it can be explained by considering the local curvature of the membrane. Our results indicate that the real area per lipid near the protein is not significantly modified with respect to the pure bilayer result. Further our results indicate that the local deformation of the membrane around the protein might be compatible with the enhanced protein activity observed in experiments over a narrow range of membrane thicknesses.     

Interaction of cetylpyridinium chloride with giant lipid vesicles

The interaction of cationic surfactant cetylpyridinium chloride, CPC, with giant lipid vesicles prepared from 1-palmitoyl-2-oleoylphosphatidylcholine, POPC, was examined at various concentrations of the lipid component. The lipid concentration was determined by a spectrophotometric method. The potentiometric method based on surfactant-selective electrode was used for the determination of surfactant concentration in the external water solution. From these results, moles of surfactant incorporated in the membrane per mole of lipid (parameter beta) and two kinds of partition coefficients were calculated. Their values were found to be considerably larger than the available literature data. A three stage process of surfactant-induced solubilization of lipid vesicles was observed. First, stable mixed bilayers form, which become saturated with CPC at a value beta(sat) larger than 0.8, which then gradually disintegrate. Just prior to the breakdown of the vesicular structure, formation of ellipsoidal vesicles was observed by optical microscopy. This phenomenon was attributed to the cooperative incorporation of surfactant into the bilayer. Fluorescence measurements have shown that the second stage in the solubilization process of POPC by the C16 chain-length surfactant does not involve mixed micelles. These are formed only in the third stage, which is the complete solubilization of POPC bilayers. The corresponding critical micellization concentration decreases with increasing concentration of the lipid component.     

Modulating membrane properties: the effect of trehalose and cholesterol on a phospholipid bilayer

The protective properties of trehalose on cholesterol-containing lipid dipalmitoylphosphatidylcholine (DPPC) bilayers are studied through molecular simulations. The ability of the disaccharide to interact with the phospholipid headgroups and stabilize the membrane persists even at high cholesterol concentrations and restricts some of the changes to the structure that would otherwise be imposed by cholesterol molecules. Predictions of bilayer properties such as area per lipid, tail ordering, and chain conformation support the notion that the disaccharide decreases the main melting transition in these multicomponent model membranes, which correspond more closely to common biological systems than pure bilayers. Molecular simulations indicate that the membrane dynamics are slowed considerably by the presence of trehalose, indicating that high sugar concentrations would serve to avert possible phase separations that could arise in mixed phospholipid systems. Various time correlation functions suggest that the character of the modifications in lipid dynamics induced by trehalose and cholesterol is different in the hydrophilic and hydrophobic regions of the membrane.     

Combined Monte Carlo and molecular dynamics simulation of hydrated 18:0 sphingomyelin-cholesterol lipid bilayers

We have carried out atomic level molecular dynamics and Monte Carlo simulations of hydrated 18:0 sphingomyelin (SM)-cholesterol (CHOL) bilayers at temperatures of 20 and 50 degrees C. The simulated systems each contained 266 SM, 122 CHOL, and 11861 water molecules. Each simulation was run for 10 ns under semi-isotropic pressure boundary conditions. The particle-mesh Ewald method was used for long-range electrostatic interactions. Properties of the systems were calculated over the final 3 ns. We compare the properties of 20 and 50 degrees C bilayer systems with each other, with experimental data, and with experimental and simulated properties of pure SM bilayers and dipalmitoyl phospatidyl choline (DPPC)-CHOL bilayers. The simulations reveal an overall similarity of both systems, despite the 30 degrees C temperature difference which brackets the pure SM main phase transition. The area per molecule, lipid chain order parameter profiles, atom distributions, and electron density profiles are all very similar for the two simulated systems. Consistent with simulations from our lab and others, we find strong intramolecular hydrogen bonding in SM molecules between the phosphate ester oxygen and the hydroxyl hydrogen atoms. We also find that cholesterol hydroxyl groups tend to form hydrogen bonds primarily with SM carbonyl, methyl, and amide moieties and to a lesser extent methyl and hydroxyl oxygens.