Search

Cardiolipins (CL) are uniquely structured double lipids that are universally found in membranes which couple electron transport and phosphorylation. The proposed roles of CL in these membranes include mainly two aspects: the effects on the structure and dynamics of the membranes' lipid component and the interplay with membrane-associated proteins. In previous work, we have studied both of these aspects by utilizing atomistic molecular dynamics simulations, starting with the general effects of CL on membrane properties and the interplay with ions [1], further moving on to the interactions with the membrane-embedded respiratory complex cyt bc1 [2]. The latter model system that is of greater interest comprises the entire cyt bc1 dimer of the purple photosynthetic bacterium Rhodobacter capsulatus embedded in a lipid bilayer, whose lipid composition mimics that of the inner mitochondrial membrane. Intriguingly CLs were observed to diffuse spontaneously to the dimer interface and to the immediate vicinity of the catalytic Qi-sites [2]. This observation is in agreement with experimental data, as CLs are indeed located close to the Qi-sites in several X-ray crystal structures of the complex. Importantly, our observations support the proposed role of CL in delivering protons for the non-reduced substrate forms in the active site.In ongoing work that we discuss here we focus more specifically on the roles of individual components of the proposed proton uptake pathway (CL, water, and individual protein residues) and on the atom-level reaction mechanism in the binding pocket. To this end, further MD simulations, QM calculations, and cite-directed mutagenesis experiments were employed. We also discuss whether there is a plausible pathway for substrate movement between the active sites through the lipid-filled insides of the complex, and the role of oxidative stress in cyt bc1 behavior.References[1] Poyry et al. J.Phys.Chem.B,113, 15513(2009).[2] Poyry et al. BBA,1827, 769(2013).

In vitro toxicological studies together with atomistic molecular dynamics simulations show that occupational co-exposure with C-60 fullerene may strengthen the health effects of organic industrial chemicals. The chemicals studied are acetophenone, benzaldehyde, benzyl alcohol, m-cresol, and toluene which can be used with fullerene as reagents or solvents in industrial processes. Potential co-exposure scenarios include a fullerene dust and organic chemical vapor, or a fullerene solution aerosolized in workplace air. Unfiltered and filtered mixtures of C-60 and organic chemicals represent different co-exposure scenarios in in vitro studies where acute cytotoxicity and immunotoxicity of C-60 and organic chemicals are tested together and alone by using human THP-1-derived macrophages. Statistically significant co-effects are observed for an unfiltered mixture of benzaldehyde and C-60 that is more cytotoxic than benzaldehyde alone, and for a filtered mixture of m-cresol and C-60 that is slightly less cytotoxic than m-cresol. Hydrophobicity of chemicals correlates with co-effects when secretion of pro-inflammatory cytokines IL-1 beta and TNF-alpha is considered. Complementary atomistic molecular dynamics simulations reveal that C-60 co-aggregates with all chemicals in aqueous environment. Stable aggregates have a fullerene-rich core and a chemical-rich surface layer, and while essentially all C-60 molecules aggregate together, a portion of organic molecules remains in water.

A molecular dynamics (MD) simulation of a fully hydrated, liquid-crystalline dimyristoylphosphatidylcholine (DMPC)-Chol bilayer membrane containing approximately 22 mol% Chol was carried out for 4.3 ns. The bilayer reached thermal equilibrium after 2.3 ns of MD simulation. A 2.0-ns trajectory generated during 2.3-4.3 ns of MD simulation was used for analyses to determine the effects of Chol on the membrane/water interfacial region. In this region, 70% of Chol molecules are linked to DMPC molecules via short-distance interactions, where the Chol hydroxyl group (OH-Chol) is 1) charge paired to methyl groups of the DMPC choline moiety ( approximately 34%), via the hydroxyl oxygen atom (Och); 2) water bridged to carbonyl ( approximately 19%) and nonester phosphate ( approximately 14%) oxygen atoms, via both Och and the hydroxyl hydrogen atom (Hch); and 3) directly hydrogen (H) bonded to carbonyl ( approximately 11%) and nonester phosphate ( approximately 5%) oxygen atoms, via Hch ( approximately 17% of DMPC-Chol links are multiple). DMPC's gamma-chain carbonyl oxygen atom is involved in 44% of water bridges and 51% of direct H bonds formed between DMPC and Chol. On average, a Chol molecule forms 0.9 links with DMPC molecules, while a DMPC molecule forms 2.2 and 0.3 links with DMPC and Chol molecules, respectively. OH-Chol makes hydrogen bonds with 1.1 water molecules, preferentially via Hch. The average number of water molecules H bonded to the DMPC headgroup is increased by 7% in the presence of Chol. These results indicate that inclusion of Chol decreases interlipid links and increases hydration in the polar region of the membrane.

A 15-ns molecular dynamics simulation of the fully hydrated liquid-crystalline dimyristoylphosphatidylcholine-cholesterol (DMPC-Chol) bilayer containing approximately 22 mol% Chol was carried out. The generated trajectory was analysed to investigate the mechanism of the Chol condensing effect on DMPC hydrocarbon chains and the influence of Chol on the chain packing in the membrane. Chol was found to induce stronger van der Waals interactions among the chains, whereas its interactions with the chains were weak. In the DMPC-Chol bilayer, as in the DMPC bilayer, DMPC chains were regularly packed around a chosen chain but around a Chol molecule they were not. DMPC gamma chains made closer contacts with Chol than the beta chains.

Using extensive atomistic simulations, we show that there is a single experimentally accessible parameter--the sterol tilt--that can be used to determine a sterol's capability to induce order, and thus to promote, e.g., formation of lipid rafts. The observations also facilitate designing new effective sterols.

Poly(ethylene glycol) (PEG) is a polymer commonly used for functionalization of drug molecules to increase their bloodstream lifetime, hence efficacy. However, the interactions between the PEGylated drugs and biomembranes are not clearly understood. In this study, we employed atomic-scale molecular dynamics (MD) simulations to consider the behavior of two drug molecules functionalized with PEG (tetraphenylporphyrin used in cancer phototherapy and biochanin A belonging to the isoflavone family) in the presence of a lipid bilayer. The commonly held view is that functionalization of a drug molecule with a polymer acts as an entropic barrier, inhibiting the penetration of the drug molecule through a cell membrane. Our results indicate that in the bloodstream there is an additional source of electrostatic repulsive interactions between the PEGylated drugs and the lipid bilayer. Both the PEG chain and lipids can bind Na(+) ions, thus effectively becoming positively charged molecules. This leads to an extra repulsive effect resulting from the presence of salt in the bloodstream. Thus, our study sheds further light on the role of PEG in drug delivery.

The molecular lipid composition of exosomes is largely unknown. In this study, sophisticated shotgun and targeted molecular lipidomic assays were performed for in-depth analysis of the lipidomes of the metastatic prostate cancer cell line, PC-3, and their released exosomes. This study, based in the quantification of approximately 280 molecular lipid species, provides the most extensive lipid analysis of cells and exosomes to date. Interestingly, major differences were found in the lipid composition of exosomes compared to parent cells. Exosomes show a remarkable enrichment of distinct lipids, demonstrating an extraordinary discrimination of lipids sorted into these microvesicles. In particular, exosomes are highly enriched in glycosphingolipids, sphingomyelin, cholesterol, and phosphatidylserine (mol% of total lipids). Furthermore, lipid species, even of classes not enriched in exosomes, were selectively included in exosomes. Finally, it was found that there is an 8.4-fold enrichment of lipids per mg of protein in exosomes. The detailed lipid composition provided in this study may be useful to understand the mechanism of exosome formation, release and function. Several of the lipids enriched in exosomes could potentially be used as cancer biomarkers.

Cytochrome (cyt) bc 1 is a multi-subunit membrane protein complex that is a vital component of the respiratory and photosynthetic electron transfer chains both in bacteria and eukaryotes. Although the complex’s dimer structure has been solved using X-ray crystallography, it has not yet been studied in large-scale classical molecular dynamics (MD) simulations. In part, this is due to lack of suitable force field parameters, centered atomic point charges in particular, for the complex’s prosthetic redox centers. Accurate redox center charges are needed to depict realistically the inter-molecular interactions at different redox stages of the cyt bc 1 complex. Accordingly, here we present high-precision atomic point charges for the metal centers of the cyt bc 1 complex of Rhodobacter capsulatus derived from extensive density functional theory calculations, fitted using the restrained electrostatic potential methodology and combined with the CHARMM force field parameters. We also provide the Hartree–Fock charges for all substrate forms (quinol, quinone, and semiquinone) and the inhibitors antimycin and stigmatellin of the bacterial bc 1 complex. The accuracy of the parameterization scheme was verified by running a 200-ns MD simulation encompassing the entire cyt bc 1 complex embedded in a lipid bilayer and solvated with explicit water. The results indicate that these meticulously derived parameters are ready for running extensive MD simulations encompassing all biologically relevant stages of the cyt bc 1 complex reaction cycle.

This work investigates the problem of selective binding of the chosen β-blockers towards human β1 and β2-adrenergic receptors (βARs). The selective blockade of the cardiac β1 has important medical applications. However, the current understanding of this prominent medical problem is still very limited. The recently published crystal structures of β1 and the β2 [1,2] offer an opportunity to understand the mechanism of β1-selective binding with atomic accuracy. Numerous computational studies have explored the binding pockets or conformational properties of these receptors through homology modeling, structure-based modeling, fragment-based screening and molecular dynamics simulations [3-5]. Yet, as far as the authors know, thus far there have been no computational studies focusing exclusively on the problem of the selective blockade of the β1 subtype. The computations carried out and reported here fill this gap. Surprisingly, our simulations show that water plays a fundamental role in the binding site of both βARs, being particularly important in determining the selectivity of binding.1. T. Warne et al., (2008) Structure of a beta1-adrenergic G-protein-coupled receptor. Nature 454, 486-91.2. V. Cherezov et al., (2007) High-resolution crystal structure of an engineered human beta2-adrenergic G protein-coupled receptor. Science 318, 1258-65.3. K. Kaszuba et al., (2010) Molecular dynamics simulations reveal fundamental role of water as factor determining affinity of binding of beta-blocker nebivolol to beta(2)-adrenergic receptor. J. Phys. Chem. B. B 114, 8374-86.4. A. Ivetacet et al., (2010) Mapping the druggable allosteric space of g-protein coupled receptors: a fragment-based molecular dynamics approach. CB&DD 76, 201-17.5. S. Vanni et al., (2011) Predicting Novel Binding Modes of Agonists to β Adrenergic Receptors using All-Atom Molecular Dynamics Simulations. Plos Comput Biol 7, 1, ee1001053.

Glycosylation where glycans are enzymatically attached to proteins is a common post-translational modification among many membrane protein families. The glycans associated with glycosylation serve a variety of structural and functional roles in protein structures, such as recognition and cell-cell adhesion. Determination of glycosylation in membrane protein structures is a complicated matter, though, since due to methodological limitations and the soft nature of glycans, the structure of the crystallized membrane proteins does not include the information arising from glycosylation. Here we use extensive atomistic molecular dynamics simulations together with experimental ligand binding assays and antibodies (1) to show for the N-glycosylated human epidermal growth factor receptor (EGFR) in biomimetic lipid bilayers that N-glycosylation is a crucial determinant of EGFR conformation, and specifically the orientation of the EGFR ectodomain relative to the membrane. It is expected that similar conclusions as to the importance of glycosylation in protein conformational behavior can be drawn for a number of other membrane proteins, too.(1) K. Kaszuba et al., Proc. Natl. Acad. Sci. USA 112, 4334 (2015).

Simulations based on molecular dynamics (MD) are an important tool in structural biology. In the field of lipid membrane research, MD represents a major tool as it grant access to atomic level features of the membrane difficult to access by experimental means. In this chapter we describe step by step how to simulate a membrane in silico. The explanation comprehends from the initial set up of a membrane system to its analysis passing through the simulation itself. Most of the tasks described here are performed using the Gromacs simulation suite, a widely used software for molecular modeling and simulations. In the selected example the attention is focused on a single component lipid bilayer composed of DPPC, one of the most studied lipid species.

The interactions of the neurotransmitter dopamine, and its precursor l-dopa, with membrane lipids were investigated through a set of molecular dynamic simulations with all atom resolution. The results obtained indicate that both dopamine and l-dopa have a pronounced association with the lipid head groups, predominantly mediated through H-bonds. As a result the molecules are anchored to the interfacial region of the membrane. The strength of this interaction is dependent on lipid composition - the presence of phosphatidylserine leads to an increase in the strength of this interaction, resulting in an H-bond network with a lifetime much longer than the timescale of our simulations. Also, bilayers that include sphingomieline and cholesterol interact strongly with dopamine and l-dopa. We postulate that the high membrane association that we have observed for both dopamine and l-dopa could have the following effects: 1) when on the plasma membrane exterior, favour the availability of these compounds for cell membrane uptake processes and, 2) when on an internal membrane surface, accentuate the importance of membrane-bound metabolizing enzymes over their soluble counterparts.

Molecular dynamics simulation has been used to study the speci fi c interactions between poly(ethylene glycol) (PEG) and three drug molecules for which PEG is used to aid delivery: paclitaxel and piroxicam, where PEG is a carrier agent, and hematoporphyrin, where PEG is covalently attached to form a “ stealth shield ” . Simulating at physiological salt concentration, we found no evidence of any speci fi c interaction between paclitaxel or piroxicam with PEG, but found a strong interaction for the case of hematoporphyrin. This interaction is lipophilic in nature, between the nonpolar (CH 2 ) 2 groups of the PEG and the porphin ring of the hematoporphyrin. This interaction was found to be strong enough that the PEG aggregated to the hematoporphyrin, independent of whether or not it was covalently bound. Interestingly, when the simulation was repeated in absence of salt we found evidence of this interaction being weakened. This led us to hypothesize a previously unforeseen mechanism: interaction with salt cations cause the PEG to coil around the salt ions, each ion binding to many PEG oxygens, increasing the exposure of the nonpolar ethylene groups, thus increasing the e ff ective hydrophobicity of PEG. The Hydrophobic ethylene groups of the PEG chains adhere strongly to the hydrophobic porphin ring. Experiments involving absorption spectra measurements were conducted, and these results also indicated that presence of salt at physiological level increases the e ff ective attractive interaction between PEG and hematoporphyrin. Taken together, our results demonstrate that while PEG, due to its solubility in both polar and nonpolar solvents, may act as a dissolution aid for paclitaxel and piroxicam, of the three drug molecules studied it will only have a protective role for the case of the hematoporphyrin.

Drug nanocarriers are often derivatized with targeting moieties to achieve site specific delivery, however, the results from this approach have, as yet, not reached expectations. We have tested a new phage display based targeting moiety, the activated endothelium targeting peptide (AETP), for its vascular endothelium directed targeting efficiency, when anchored to a PEGylated liposome via maleimide chemistry. Our results have, however, not shown any evidence of improved targeting. We have hypothesized that the failure of the AETP moiety is due to its availability to target receptors being restricted, as a result of steric hindrance due to the PEG polymer, and possibly affinity for bloodstream proteins, particularly human serum albumin (HSA). In this context, molecular modeling was used to contrast the properties of the AETP moiety to those of the RGD targeting peptide, already found to be effective in previous trials. Our molecular dynamics simulation results indicate the AETP moiety is located within the PEG layer, and its hydrophobic nature causes it to be obscured by PEG to a greater extent than the more hydrophilic RGD targeting peptide. Protein–ligand docking results indicated similar affinities for HSA of both the AETP moiety and a PEG fragment, and a significantly lower affinity for the RGD peptide. We know of no means to investigate this experimentally with atomic level resolution, thus our use of computational methods to investigate this can be seen as a new tool for rational design in nanomedicine.

We use atom-scale molecular dynamics simulations to clarify the role of glycosphingolipids in the dynamics of cholesterol-rich lipid rafts. To this end, we consider lipid membranes that contain varying. amounts of galactosylceramide (GalCer), sphingomyelin, cholesterol, and phosphatidylcholine. The results indicate that increasing the portion of GalCer molecules greatly slows down the lateral diffusion, Only 5-10 mol % of GalCer causes a decrease of almost an order of magnitude compared to corresponding membranes without GalCer. The slowing down is not related to interdigitation, which becomes weaker with increasing GalCer concentration. Instead, the decrease in diffusion is found to correlate with the increasing number of hydrogen bonds formed between GalCer and the phospholipid molecules, which is also observed to have other effects, such as to increase the friction between the membrane leaflets.

Theoretical studies predict hydrophobic matching between transmembrane domains of proteins and bilayer lipids to be a physical mechanism by which membranes laterally self-organize. We now experimentally study the direct consequences of mismatching of transmembrane peptides of different length with bilayers of different thicknesses at the molecular level. In both model membranes and simulations we show that cholesterol critically constrains structural adaptations at the peptide-lipid interface under mismatch. These constraints translate into a sorting potential and lead to selective lateral segregation of peptides and lipids according to their hydrophobic length.

We used atomistic simulations to study the membrane-bound form of catechol-O-methyltransferase (MB-COMT). In particular we investigated the 26-residue transmembrane a-helical segment of MB-COMT together with the 24-residue fragment that links the transmembrane component to the main protein unit that was not included in our model. In numerous independent simulations we observed the formation of a salt bridge between ARC 27 and GLU40. The salt bridge closed the flexible loop that formed in the linker and kept it in the vicinity of the membrane-water interface. All simulations supported this conclusion that the linker has a clear affinity for the interface and preferentially arranges its residues to reside next to the membrane, without a tendency to relocate into the water phase. Furthermore, an extensive analysis of databases for sequences of membrane proteins that have a single transmembrane helical segment brought about an interesting view that the flexible loop observed in our work can be a common structural element in these types of proteins. In the same spirit we close the article by discussing the role of salt bridges in the formation of three-dimensional structures of membrane proteins that exhibit a single transmembrane helix.

Identification of early mechanisms that may lead from obesity towards complications such as metabolic syndrome is of great interest. Here we performed lipidomic analyses of adipose tissue in twin pairs discordant for obesity but still metabolically compensated. In parallel we studied more evolved states of obesity by investigating a separated set of individuals considered to be morbidly obese. Despite lower dietary polyunsaturated fatty acid intake, the obese twin individuals had increased proportions of palmitoleic and arachidonic acids in their adipose tissue, including increased levels of ethanolamine plasmalogens containing arachidonic acid. Information gathered from these experimental groups was used for molecular dynamics simulations of lipid bilayers combined with dependency network analysis of combined clinical, lipidomics, and gene expression data. The simulations suggested that the observed lipid remodeling maintains the biophysical properties of lipid membranes, at the price, however, of increasing their vulnerability to inflammation. Conversely, in morbidly obese subjects, the proportion of plasmalogens containing arachidonic acid in the adipose tissue was markedly decreased. We also show by in vitro Elovl6 knockdown that the lipid network regulating the observed remodeling may be amenable to genetic modulation. Together, our novel approach suggests a physiological mechanism by which adaptation of adipocyte membranes to adipose tissue expansion associates with positive energy balance, potentially leading to higher vulnerability to inflammation in acquired obesity. Further studies will be needed to determine the cause of this effect.

Cytochrome (cyt) bc1 complex has two active sites where energy-sustaining electron/proton transfer reactions occur. To determine the exact transfer steps involving the substrate quinone, the cyt bc1 complex of Rhodobacter capsulatus was studied using atomistic molecular dynamics simulations. The computations give clear-cut ideas how the proton transfer might take place at both of the active sites. The simulations strongly suggest that water-mediated interactions are needed to achieve quinol oxidation and to prevent short-circuit suppression at the Qo-site. In the case of the Qi-site, the simulations emphasize the role of negatively charged cardiolipin in quinone reduction. Together with site-directed mutagenesis this in silico approach provides long-awaited information regarding the structure-function relationship of the cyt bc1 complex.

Sphingomyelins (SMs) and sterols are important constituents of the plasma membrane and have also been identified as major lipid components in membrane rafts. Using SM analogs with decreasing headgroup methylation, we systemically analyzed the effect of headgroup size on membrane properties and interactions with cholesterol. An increase in headgroup size resulted in a decrease in the main phase transition. Atom-scale molecular-dynamics simulations were in agreement with the fluorescence anisotropy experiments, showing that molecular areas increased and acyl chain order decreased with increasing headgroup size. Furthermore, the transition temperatures were constantly higher for SM headgroup analogs compared to corresponding phosphatidylcholine headgroup analogs. The sterol affinity for phospholipid bilayers was assessed using a sterol-partitioning assay and an increased headgroup size increased sterol affinity for the bilayer, with a higher sterol affinity for SM analogs as compared to phosphatidylcholine analogs. Moreover, the size of the headgroup affected the formation and composition of cholesterol-containing ordered domains. Palmitoyl-SM (the largest headgroup) seemed to attract more cholesterol into ordered domains than the other SM analogs with smaller headgroups. The ordering and condensing effect of cholesterol on membrane lipids was also largest for palmitoyl-SM as compared to the smaller SM analogs. The results show that the size of the SM headgroup is crucially important for SM-SM and SM-sterol interactions. Our results further emphasize that interfacial electrostatic interactions are important for stabilizing cholesterol interactions with SMs.

International audience; Interactions between small organic molecules and lipid or cell membranes are important because of their role in the distribution of biologically active substances inside the membrane and their permeation through the cell membranes. In the current paper, we have explored the effect of the attachment of long hydrocarbon tails on the behavior of small organic molecule inside the lipid membrane. Naphthalene with two decyloxy groups attached at the opposite sites of the ring (2,6-bis(decyloxy)naphthalene, 3) was synthesized and incorporated into phosphatidylcholine (PC) vesicles. Fluorescence methods as well as molecular dynamic (MD) simulations were used to estimate the position, orientation, and migration of compound 3 in PC bilayer. It was found that the naphthalene ring of compound 3 resides in the upper acyl chain region of the bilayer and the hydrocarbon tails are directed to the center of the bilayer. As was shown with cryotransmission electron microscopy (cryo-TEM), such lipidlike conformation enables compound 3 to be incorporated into liposomes at a very high content without their disintegration. Moreover, compound 3 can migrate from one leaflet to other. The mechanism of this process is, however, different from that characteristic of the flip-flop event of lipid molecules in the membrane. Finally, the possible application of compound 3 as a rotational molecular probe for monitoring fluidity of liposomal membrane in the acyl side chain region was checked by studies of the effect of cholesterol on the fluorescence anisotropy of 3.

1 - Article; Phospholipid-based reverse micelles are composed of branched cylinders. Their branching points are known to attract themselves and to slide along branches. The rate of this sliding is governed by the lifetime of H(D)-bonded water bridges between phospholipid molecules. This lifetime is increased when the water is deuterated. On condition that the water contains at least 40 D atoms%, water/dipalmitoylphosphatidylcholine (DPPC)/deuterated pyridine reverse micelles with the composition 1.11:250 (v/v) have been shown to self-organize into a liquid crystal in the 310-316 K temperature range. The mechanism of this self-organization is unraveled by following the FTIR and (1)H NMR spectra of more concentrated micelles upon heating. During the preparation of micelles, pyridine-(D(+))H(+) ions are formed. They give rise to hydron transfers, under the influence of the DPPC electric charges, evidenced by two broad FTIR absorptions above (BB1) and below (13132) the nu(C-0) stretch. These hydron transfers occur along strong (D(+))H(+) bonds of pyridinium ions with pyridine (BB1) and DPPC C=O groups (BB2). The proton transfers at the interface of micelles, relayed in the continuous pyridine medium, create a tenuous link between separated micelles, thus facilitating their organization. Upon heating, DPPC heads shrink and DPPC chains expand to make wedge-shaped DPPC molecules. The micelles then change in shape: cylinders constrict and enclosed water drifts towards branching points, which swell. Branching points of neighboring micelles come into contact. Due to the deuteration of water these contacts are prolonged and H bonds are formed between DPPC molecules located in each branching point. Upon storage at 39 degrees C, these branching points fuse. The lateral diffusion of DPPC molecules becomes free, as evidenced by a narrowing of all (1)H NMR resonances. Upon further heating, reorganization into a liquid crystal occurs.

We used the crystal structure of prolyl oligopeptidase (POP) with bound Z-pro-prolinal (ZPP) inhibitor (Protein Data Bank (PDB) structure 1QFS) to perform an intensive molecular dynamics study of the POP-ZPP complex. We performed 100 ns of simulation with the hemiacetal bond, through which the ZPP is bound to the POP, removed in order to better investigate the binding cavity environment. From basic analysis, measuring the radius of gyration, root mean square deviation, solvent accessible surface area and definition of the secondary structure of protein, we determined that the protein structure is highly stable and maintains its structure over the entire simulation time. This demonstrates that such long time simulations can be performed without the protein structure losing stability. We found that water bridges and hydrogen bonds play a negligible role in binding the ZPP thus indicating the importance of the hemiacetal bond. The two domains of the protein are bound by a set of approximately 12 hydrogen bonds, specific to the particular POP protein. Keywords: prolyl oligopeptidase, POP, molecular dynamics, hydrogen bonds, protein dynamics, inhibitor ligands, hemiacetal bond

Despite the large number of possible glycerol-based phospholipids, biological membranes contain only a small number of them. For example, double bonds in acyl chains are preferably located in the sn-2 chain. The question that emerges is: Why? We have addressed this question through atomistic simulations by considering pure one-component bilayers comprising monounsaturated glycerophospholipids [1-stearoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (SOPC) and 1-oleoyl-2-stearoyl-sn-glycero-3-phosphatidylcholine (OSPC)] and membranes of these lipids mixed with cholesterol. By considering the cases in which an individual double bond is in either the sn-1 or the sn-2 chain, we elucidated how membrane properties depend on this intrinsic feature. We found small but systematic differences in all structural and dynamic membrane properties that we considered. It turns out that the differences are driven by two factors: the mismatch in the acyl chain lengths and the interaction of the double bond in the acyl chains with the cholesterol off-plane methyl groups. The results highlight the fact that unsaturated sn-2 chains lead to more disordered membranes than systems with unsaturated sn-1 chains. The differences between the two isomers are enhanced when cholesterol is present as a result of the interaction of the off-plane cholesterol methyl groups with the double-bond carbon segments in the lipid acyl chains.

Physicochemical properties of heavy water (D2O) differ to some extent from those of normal water. Substituting D2O for H2O has been shown to affect the structural and dynamic properties of proteins, but studies of its effects on lipid bilayers are scarce. In this paper, the atomic level molecular dynamics (MD) simulation method was used to determine the effects of this substitution on the properties of a dipalmitoylphosphatidylcholine (DPPC) bilayer and its hydrating water. MD simulations of two DPPC bilayers, one fully hydrated with H2O and the other with D2O, were carried out for over 50 ns. For H2O, the simple point charge (SPC) model was used, and for D2O, the extended SPC-HW model was employed. Analyses of the simulation trajectories indicate that several properties of the membrane core and the membrane/water interface are affected by replacing H2O by D2O. However, the time-averaged properties, such as membrane compactness, acyl chain order, and numbers of PC-water H (D)-bonds and PC-PC water bridges, are much less affected than time-resolved properties. In particular, the lifetimes of these interactions are much longer for D2O molecules than for H2O ones. These longer lifetimes results in a slightly better ordering of the D2O molecules and average self-diffusion, which is 50% slower compared with the H2O molecules. This large isotope effect has been assigned to the repercussions of the longer lived D-bonding to DPPC headgroups insofar as all water molecules sense the presence of the DPPC bilayer.

Without any exaggeration, cholesterol is one of the most important lipid species in eukaryotic cells. Its effects on cellular membranes and functions range from purely mechanistic to complex metabolic ones, besides which it is also a precursor of the sex hormones (steroids) and several vitamins. In this review, we discuss the biophysical effects of cholesterol on the lipid bilayer, in particular the ordering and condensing effects, concentrating on the molecular level or inter-atomic interactions perspective, starting from two-component systems and proceeding to many-component ones e.g., modeling lipid rafts. Particular attention is paid to the roles of the methyl groups in the cholesterol ring system, and their possible biological function. Although our main research methodology is computer modeling, in this review we make extensive comparisons between experiments and different modeling approaches. Keywords: Sterol; Lipid bilayer; Unsaturated lipid; Sphingolipid; Lipid raft; Molecular dynamics simulation

We have employed atomistic molecular dynamics simulations to investigate the effect of double-bond parametrization on lipid membrane properties. As models, we use one-component membranes composed of glycerol-based phosphatidylcholines (PCs) with monounsaturated acyl chains, and we complement these studies by additional PC/cholesterol simulations. We compare differences between double-bond parametrizations by varying the position of the double bond systematically along the lipid hydrocarbon chains. The results give rise for concern: They indicate that the double-bond description may change not only the quantitative but also the qualitative nature of membrane behavior. In particular, we find that the double-bond description which accounts for skew states in the vicinity of a double bond predicts a maximum in membrane disorder, when the double bond resides at the middle of an acyl chain, in agreement with experiments. The more commonly used description which does not accommodate skew states, however, predicts membrane disorder to decrease monotonically as the double bond is shifted from the glycerol backbone to the end of an acyl chain. The results highlight the importance of properly describing double bonds especially in many-component membranes, where the interplay of different molecule types is difficult to predict on intuitive grounds.

The structural and dynamical properties of lipid membranes rich in phospholipids and cholesterol are known to be strongly affected by the unsaturation of lipid acyl chains. We show that not only unsaturation but also the position of a double bond has a pronounced effect on membrane properties. We consider how cholesterol interacts with phosphatidylcholines comprising two 18-carbon long monounsaturated acyl chains, where the position of the double bond is varied systematically along the acyl chains. Atomistic molecular dynamics simulations indicate that when the double bond is not in contact with the cholesterol ring, and especially with the C18 group on its rough β-side, the membrane properties are closest to those of the saturated bilayer. However, any interaction between the double bond and the ring promotes membrane disorder and fluidity. Maximal disorder is found when the double bond is located in the middle of a lipid acyl chain, the case most commonly found in monounsaturated acyl chains of phospholipids. The results suggest a cholesterol-mediated lipid selection mechanism in eukaryotic cell membranes. With saturated lipids, cholesterol promotes the formation of highly ordered raft-like membrane domains, whereas domains rich in unsaturated lipids with a double bond in the middle remain highly fluid despite the presence of cholesterol.

Despite extensive studies, the remarkable structure-function relationship of cholesterol in cellular membranes has remained rather elusive. This is exemplified by the fact that the membrane properties of cholesterol are distinctly different from those of many other sterols. Here we elucidate this issue through atomic-scale simulations of desmosterol and 7-dehydrocholesterol (7DHC), which are immediate precursors of cholesterol in its two distinct biosynthetic pathways. While desmosterol and 7DHC differ from cholesterol only by one additional double bond, we find that their influence on saturated lipid bilayers is substantially different from cholesterol. The capability to form ordered regions in a saturated (dipalmitoyl-phosphatidylcholine) membrane is given by cholesterol>7DHC>desmosterol, indicating the important role of cholesterol in saturated lipid environments. For comparison, in an unsaturated (dioleoyl-phosphatidylcholine) bilayer, the membrane properties of all sterols were found to be essentially identical. Our studies indicate that the different membrane ordering properties of sterols can be characterized by a single experimentally accessible parameter, the sterol tilt. The smaller the tilt, the more ordered are the lipids around a given sterol. The molecular level mechanisms responsible for tilt modulation are found to be related to changes in local packing around the additional double bonds.

Cholesterol is an indispensable molecule in mammalian cell membranes. To truly understand its role in the functions of membranes, it is essential to unravel cholesterol's structure-function relationship determined by underlying molecular interactions. For this purpose, we elaborate on this issue by considering the previously proposed idea that cholesterol's effects on a number of physical properties of membranes have been optimized during the evolution by removal of its excess methyl groups from the alpha-face of cholesterol, thus "smoothening" the structure. Consequently, the methyl groups still attached to cholesterol are one of the most intriguing structural features of the molecule. An obvious question arises: Why do these methyl groups still exist, and could cholesterol properties be further optimized by their removal? Because of the nature of the biosynthetic pathways of cholesterol, and the evidence of decreased interactions between sterols and lipid acyl chains when methyl groups are present, it seems plausible that removal of the methyl groups might indeed lead to stronger ordering and condensing effects of the cholesterol molecule. Atomic-scale molecular dynamics simulations of numerous modified sterols embedded in saturated lipid bilayers demonstrate, however, that the issue is more subtle. The analysis reveals a complex interplay between the lipid acyl chains and the structural details of cholesterol. Changes in cholesterol structure typically do not improve its performance in terms of promoting membrane order. This view is substantiated by a detailed analysis of the simulation data. In particular, it highlights the importance of the methyl group C18 for cholesterol properties. The C18 group resides between the third and fourth ring of cholesterol on its "rough" beta-side, and the results provide compelling evidence that C18 is crucial for the proper orientation of the sterol. More generally, the data provide insight into the role of the methyl groups of cholesterol.

Phosphatidylcholines (PCs) are among the most common phospholipids in plasma membranes. Their structural and dynamic properties are known to be strongly affected by unsaturation of lipid hydrocarbon chains, yet the role of the exact positions of the double bonds is poorly understood. In this work, we shed light on this matter through atomic-scale molecular dynamics simulations of eight different one-component lipid bilayers comprised of PCs with 18 carbons in their acyl chains. By introducing a single double bond in each acyl chain and varying its position in a systematic manner, we elucidate the effects of a double bond on various membrane properties. Studies in the fluid phase show that a number of membrane properties depend on the double bond position. In particular, when the double bond in an acyl chain is located close to the membrane-water interface, the area per lipid is considerably larger than that found for a saturated lipid. Further, when the double bond is shifted from the interfacial region toward membrane center, the area per lipid is observed to increase and have a maximum when the double bond is in the middle of the chain. Beyond this point, the surface area decreases systematically as the double bond approaches membrane center. These changes in area per lipid are accompanied by corresponding changes in membrane thickness and ordering of the chains. Further changes are observed in the tilt angles of the chains, membrane hydration together with changes in the number of gauche conformations, and direct head group interactions. All of these effects can be associated with changes in acyl chain conformations and local effects of the double bond on the packing of the surrounding atoms.

Although sterols constitute one of the most important molecular species in cells, the reasons for their structure-function relationships in lipid membranes are not well understood. The main objective of this work is to elucidate the recently suggested possibility that the ordering and condensing effects of sterols on phospholipid membranes are related to the smoothness of a sterol. We focus on cholesterol, which has two methyl groups attached to its beta-face, and compare its properties to those of demethylated cholesterol (Dchol), from which the two methyl groups have been removed. Atomic-scale molecular dynamics simulations of lipid membranes comprised of saturated lipids and sterols, either cholesterol or Dchol, provide compelling evidence that despite its smoother structure, the ordering and condensing effects of Dchol are less effective than those of cholesterol. The ordering capability of both cholesterol and Dchol is highly asymmetric with respect to their ring structure, but whereas cholesterol favors the alpha-face, Dchol favors the beta-face. The origin and implications of this difference are analyzed in detail. The picture that emerges from this study supports a view that the two methyl groups at the steroid ring system of cholesterol play an important role in cholesterol-lipid interactions by reducing sterol tilt in the bilayer and hence allowing for an optimal orientation for cholesterol.

Anionic palmitoyloleoylphosphatidylglycerol (POPG) is one of the most abundant lipids in nature, yet its atomic-scale properties have not received significant attention. Here we report extensive 150-ns molecular dynamics simulations of a pure POPG lipid membrane with sodium counterions. It turns out that the average area per lipid of the POPG bilayer under physiological conditions is approximately 19% smaller than that of a bilayer built from its zwitterionic phosphatidylcholine analog, palmitoyloleoylphosphatidylcholine. This suggests that there are strong attractive interactions between anionic POPG lipids, which overcome the electrostatic repulsion between negative charges of PG headgroups. We demonstrate that interlipid counterion bridges and strong intra- and intermolecular hydrogen bonding play a key role in this seemingly counterintuitive behavior. In particular, the substantial strength and stability of ion-mediated binding between anionic lipid headgroups leads to complexation of PG molecules and ions and formation of large PG-ion clusters that act in a concerted manner. The ion-mediated binding seems to provide a possible molecular-level explanation for the low permeability of PG-containing bacterial membranes to organic solvents: highly polar interactions at the water/membrane interface are able to create a high free energy barrier for hydrophobic molecules such as benzene.

Spin-labeled stearic acid species are commonly used for electron paramagnetic resonance (EPR) studies of cell membranes to investigate phase transitions, fluidity, and other physical properties. In this paper, we use large-scale molecular dynamics simulations to investigate the position and behavior of nitroxide spin labels attached to stearic acid molecules in dipalmitoylphosphatidylcholine (DPPC) bilayers. The results of these studies are potentially very important for the interpretation of EPR spectra, which rely on assumptions about the position of the label in the membrane. Additionally, we investigate the effect of chirality and ionization of the carboxyl group of the label. For a non-ionized species, we observe that spin-label molecules are even able to make flip-flop transitions between the leaflets of the bilayer. Such transitions have been previously observed only in very rare cases in molecular simulations.

Stearoylsphingomyelin (SSM) bilayers containing 0, 22, and 50 mol % cholesterol (Chol) and a pentadecanoyl-stearoylphosphatidylcholine (15SPC) bilayer containing 22 mol % Chol were molecular dynamics simulated at two temperatures (37 degrees C and 60 degrees C). 15SPC is the best PC equivalent of SSM. The Chol effect on the SSM bilayer differs significantly from that on the 15SPC bilayer. At the same temperature and Chol content, H-bonding of Chol with SSM is more extensive than with 15SPC. SSM-Chol H-bonding anchors the OH group of Chol in the lower regions of the SSM-Chol bilayer interface. Such a location strengthens the influence of Chol on the SSM chains. In effect, the phase of the SSM-Chol bilayer containing 22 mol % Chol at 37 degrees C is shifted from the gel to the liquid-ordered phase, and the bilayer displays similar properties below and above the main phase-transition temperature for a pure SSM bilayer of approximately 45 degrees C. In contrast, due to a higher location, Chol is not able to change the phase of the 15SPC-Chol bilayer, which at 37 degrees C remains in the gel phase. Chol affects both the core and interface of the SSM bilayer. With increasing Chol content, the order of SSM chains and hydration of SSM headgroups increase, whereas polar interactions between lipids decrease.

A series of multi-nanosecond molecular dynamics (MD) simulations of wild-type cytochrome c and its spin-labeled variants with the methanethiosulfonate moiety attached at position C102 were performed (1) to elucidate the effect of the spin probe presence on the protein structure and (2) to describe the structure and dynamics of the spin-label moiety. Comparisons with the reference crystal structure of cytochrome c (PDB entry: 1YCC) indicate that the protein secondary structure is well preserved during simulations of the wild-type cytochrome c but slightly changed in simulations of the cytochrome c labeled at position C102. At the time scale covered in our simulations, the spin label exhibits highly dynamical behavior. The number of observed distinct conformations of the spin label moiety is between 3 and 13. The spin probe was found to form short-lived hydrogen bonds with the protein. Temporary hydrophobic interactions between the probe and the protein were also detected. The MD simulations directly show that the disulfide bond in the tether linking a spin probe with a protein strongly influence the behavior of the nitroxide group. The conformational flexibility and interaction with the protein are different for each of the two low energy conformations of the disulfide bond.

Atomistic molecular dynamics simulations of fully hydrated 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC), 1-palmitoyl-2-oleoyl-phosphatidylethanolamine (POPE), and 1-palmitoyl-2-oleoyl-phosphatidylglycerol (POPG) bilayers in the liquid-crystalline state were carried out to investigate the effect of different lipid headgroups on the dynamics of water at the bilayer surface in short 80 ps time scales. Results obtained in these studies show that the hydrogen bonding amine group of POPE and the glycerol group of POPG slow water motion more than the equivalent choline group of POPC. Therefore, it is surprising that the effect of a POPC bilayer surface on water dynamics is similar to that of POPE and POPG bilayers. That result is due to a much higher number of water molecules interacting with the choline group of POPC than hydrogen-bonded molecules interacting with amine or glycerol groups of POPE and POPG.

We report evidence of dense, ordered nanodomains in single-component fluid lipid bilayers. Our atomic-scale molecular dynamics simulations suggest that the area available to a lipid acyl chain exhibits large fluctuations, resulting in denser and sparser domains. The sizes of the dense domains can be up to approximately 10 nm, and their lifetimes are of the order of approximately 10 ns. In addition, our simulations suggest that domains of lipids with highly ordered acyl chains form predominantly within the dense regions, their sizes ranging from a few chains up to a few nanometers, and with lifetimes between approximately 10 ps-10 ns. These observations shed light on the origin of experimentally observed fluctuations, as well as on the mechanisms of phase transitions in lipid membranes.

Desmosterol is an immediate precursor of cholesterol in the Bloch pathway of sterol synthesis and an abundant membrane lipid in specific cell types. The significance of the difference between the two sterols, an additional double bond at position C24 in the tail of desmosterol, is not known. Here, we provide evidence that the biophysical and functional characteristics of the two sterols differ and that this is because the double bond at C24 significantly weakens the sterol ordering potential. In model membranes, desmosterol was significantly weaker than cholesterol in promoting the formation or stability of ordered domains, and in mammalian cell membranes, desmosterol associated less avidly than cholesterol with detergent-resistant membranes. Atomic scale molecular dynamics simulations showed that the double bond gives rise to additional stress in the tail, creating a rigid structure between C24 and C27 and favoring tilting of desmosterol distinct from cholesterol. Functional effects of desmosterol in cell membranes were assessed upon acutely exchanging approximately 70% of cholesterol to desmosterol. This led to impaired raft-dependent signaling via the insulin receptor, whereas non-raft-dependent protein secretion was not affected. We suggest that the choice of cholesterol synthesis route may provide a physiological mechanism to modulate raft-dependent functions in cells.

A 15-ns molecular dynamics simulation of the fully hydrated dimyristoylphosphatidylcholine-cholesterol (DMPC-Chol) bilayer containing approximately 22 mol% Chol was carried out. An 8-ns trajectory was analysed to investigate the effect of Chol on the chain packing in the bilayer core. While the packing of DMPC chains on the smooth alpha-face side of the Chol ring is similar to that in the pure DMPC bilayer, the packing on the rough beta-face side is less regular and less tight. Two methyl groups located on the Chol beta-face disturb the packing; in effect, van der Waals (vdW) interactions between Chol rings and DMPC chains are weaker than the ones between sole DMPC chains. VdW interactions between an alkyl chain of DMPC and an isooctyl tail of Chol are similarly strong as those between two DMPC chains.

Magainin-2 is a natural peptide that kills bacteria at concentrations that are harmless to animal cells. Molecular modelling methods were applied to investigate basic mechanisms of magainin-2 amide (M2a) membrane selectivity. Three computer models of a lipid matrix of either animal or bacterial membrane containing M2a were constructed and simulated. Specific interactions between membrane lipids and M2a peptides, responsible for M2a selectivity are elucidated.

Molecular dynamics simulations of two monounsaturated phosphatidylcholine (PC) bilayers made of 1-palmitoyl-2-oleoyl-PC (POPC; cis-unsaturated) and 1-palmitoyl-2-elaidoyl-PC (PEPC; trans-unsaturated) were carried out to investigate the effect of a double bond in the PC beta-chain and its conformation on the bilayer core. Four nanosecond trajectories were used for analyses. A fully saturated 1,2-dimyristoyl-PC (DMPC) bilayer was used as a reference system. In agreement with experimental data, this study shows that properties of the PEPC bilayer are more similar to those of the DMPC than to the POPC bilayer. The differences between POPC and PEPC bilayers may be attributed to the different ranges of angles covered by the torsion angles beta10 and beta12 of the single bonds next to the double bond in the oleoyl (O) and elaidoyl (E) chains. Broader distributions of beta10 and beta12 in the E chain than in the O chain make the E chain more flexible. In effect, the packing of chains in the PEPC bilayer is similar to that in the DMPC bilayer, whereas that in the POPC bilayer is looser than that in the DMPC bilayer. The effect of the cis-double bond on torsions at the beginning of the O chain (beta4 and beta5) is similar to that of cholesterol on these torsions in a myristoyl chain.

Molecular dynamics (MD) simulations complement experimental methods in studies of the structure and dynamics of lipid bilayers. The choice of algorithms employed in this computational method represents a trade-off between the accuracy and real calculation time. The largest portion of the simulation time is devoted to calculation of long-range electrostatic interactions. To speed-up evaluation of these interactions, various approximations have been used. The most common ones are the truncation of long-range interactions with the use of cut-offs, and the particle-mesh Ewald (PME) method. In this study, several multi-nanosecond cut-off and PME simulations were performed to establish the influence of the simulation protocol on the bilayer properties. Two bilayers were used. One consisted of neutral phosphatidylcholine molecules. The other was a mixed lipid bilayer consisting of neutral phosphatidylethanolamine and negatively charged phosphatidylglycerol molecules. The study shows that the cut-off simulation of a bilayer containing charge molecules generates artefacts; in particular the mobility and order of the charged molecules are vastly different from those determined experimentally. In the PME simulation, the bilayer properties are in general agreement with experimental data. The cut-off simulation of bilayers containing only uncharged molecules does not generate artefacts, nevertheless, the PME simulation gives generally better agreement with experimental data.

International audience; We present a new implementation of the program nMoldyn, which has been developed for the computation and decomposition of neutron scattering intensities from Molecular Dynamics trajectories (Comp. Phys. Commun 1995, 91, 191-214). The new implementation extends the functionality of the original version, provides a much more convenient user interface (both graphical/interactive and batch), and can be used as a tool set for implementing new analysis modules. This was made possible by the use of a high-level language, Python, and of modern object-oriented programming techniques. The quantities that can be calculated by nMoldyn are the mean-square displacement, the velocity autocorrelation function as well as its Fourier transform (the density of states) and its memory function, the angular velocity autocorrelation function and its Fourier transform, the reorientational correlation function, and several functions specific to neutron scattering: the coherent and incoherent intermediate scattering functions with their Fourier transforms, the memory function of the coherent scattering function, and the elastic incoherent structure factor. The possibility to compute memory function is a new and powerful feature that allows to relate simulation results to theoretical studies.

Epicholesterol (Echol) is an epimeric form of cholesterol (Chol). A molecular dynamics simulation of the fully hydrated dimyristoylphosphatidylcholine-Echol (DMPC-Echol) bilayer membrane containing approximately 22 mol % of Echol was carried out for 5 ns. A 3-ns trajectory generated between 2 and 5 ns of molecular dynamics simulation was used for analyses to determine the effects of Echol on the membrane properties. As reference systems, pure DMPC and mixed DMPC-Chol bilayers were used. The study shows that Echol, like Chol, changes the organization of the bilayer/water interface and increases membrane order and condensation, but to a lesser degree. Effects of both sterols are based on the same atomic level mechanisms; their different strength arises from different vertical localizations of Echol and Chol hydroxyl groups in the membrane/water interface.