Search

Purpose/Aim : The goal of this study is to reveal how age-related changes in phospholipid (PL) composition in the fiber cell plasma membranes of the human eye lens affect the cholesterol (Chol) content at which Chol bilayer domains (CBDs) and Chol crystals start to form. Materials and Methods : Saturation-recovery electron paramagnetic resonance with spin-labeled cholesterol analogs and differential scanning calorimetry were used to determine the Chol contents at which CBDs and cholesterol crystals, respectively, start to form in in membranes made of the major PL constituents of the plasma membrane of the human eye lens fiber cells. To preserve compositional homogeneity throughout the membrane suspension, the lipid multilamellar dispersions investigated in this work were prepared using a rapid solvent exchange method. The cholesterol content changed from 0 to 75 mol%. Results : The saturation recovery electron paramagnetic resonance results show that CBDs start to form at 33, 50, 46, and 48 mol% Chol in the phosphatidylethanolamine, phosphatidylcholine, phosphatidylserine, and sphingomyelin bilayers, respectively. The differential scanning calorimetry results show that Chol crystals start to form at 50, 66, 70, and 66 mol% Chol in the phosphatidylethanolamine, phosphatidylcholine, phosphatidylserine, and sphingomyelin bilayers, respectively. Conclusions : These results, as well those of our previous studies, indicate that the formation of CBDs precedes the formation of Chol crystals in all of the studied systems, and the appearance of each depends on the type of PL forming the bilayer. These findings contribute to a better understanding of the molecular mechanisms involved in the regulation of Chol-dependent processes in eye lens fiber cell membranes.

Background: Activated coagulation factor X (FXa) is the serine protease component of prothrombinase, the physiological activator of prothrombin. Factor X Nottingham (A404T) and Taunton (R405G) are two naturally occurring mutations, identified in families with a bleeding phenotype., Objective: To characterize these FX variants functionally., Methods: The activity and inhibition of recombinant FX variants were quantified in plasma-based and pure component assays., Results: The prothrombin times in FX-depleted plasma supplemented with FX Nottingham and Taunton were greatly increased compared to that of wild-type (WT) FX. Kinetic investigations of activated variants in the prothrombinase complex showed k cat /K m reduced ~50-fold and ~5-fold, respectively, explaining the prolonged prothrombin time (PT). The substituted residues are located in the protease domain Na + -binding loop, important for the activity of FXa, as well as its inhibition. Both FXa Nottingham and Taunton showed reduced affinity for Na + . Plasma-based thrombin generation assays triggered with 1 pmol/L tissue factor (TF) demonstrated only small differences in activities compared to WT FX, but large reductions at 10 pmol/L TF. Severely reduced inhibition of both FXa Nottingham and Taunton by tissue factor pathway inhibitor (TFPI) and antithrombin (AT), was shown in pure-component FXa inhibition assays. Factor Xa Nottingham and Taunton produced higher amounts of thrombin than WT FXa in pure-component prothrombinase assays in the presence of TFPI and AT, explaining the results from the plasma-based assay., Conclusions: Factor X Nottingham and Taunton both display decreased proteolytic activity. However, their reduced activity in plasma triggered by low TF can be rescued by decreased inhibition by the natural FXa inhibitors, TFPI and AT., (© 2019 International Society on Thrombosis and Haemostasis.)

Background: Activated protein C (APC)-mediated inactivation of factor (F)Va is greatly enhanced by protein S. For inactivation to occur, a trimolecular complex among FVa, APC, and protein S must form on the phospholipid membrane. However, direct demonstration of complex formation has proven elusive., Objectives: To elucidate the nature of the phospholipid-dependent interactions among APC, protein S, and FVa., Methods: We evaluated binding of active site blocked APC to phospholipid-coated magnetic beads in the presence and absence of protein S and/or FVa. The importance of protein S and FV residues were evaluated functionally., Results: Activated protein C alone bound weakly to phospholipids. Protein S mildly enhanced APC binding to phospholipid surfaces, whereas FVa did not. However, FVa together with protein S enhanced APC binding (>14-fold), demonstrating formation of an APC/protein S/FVa complex. C4b binding protein-bound protein S failed to enhance APC binding, agreeing with its reduced APC cofactor function. Protein S variants (E36A and D95A) with reduced APC cofactor function exhibited essentially normal augmentation of APC binding to phospholipids, but diminished APC/protein S/FVa complex formation, suggesting involvement in interactions dependent upon FVa. Similarly, FVa Nara (W1920R), an APC-resistant FV variant, also did not efficiently incorporate into the trimolecular complex as efficiently as wild-type FVa. FVa inactivation assays suggested that the mutation impairs its affinity for phospholipid membranes and with protein S within the complex., Conclusions: FVa plays a central role in the formation of its inactivation complex. Furthermore, membrane proximal interactions among FVa, APC, and protein S are essential for its cofactor function., (© 2019 The Authors. Journal of Thrombosis and Haemostasis published by Wiley Periodicals, Inc. on behalf of International Society on Thrombosis and Haemostasis.)

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.

Rationale: Increased expression of CLIC4 (chloride intracellular channel 4) is a feature of endothelial dysfunction in pulmonary arterial hypertension, but its role in disease pathology is not fully understood., Objective: To identify CLIC4 effectors and evaluate strategies targeting CLIC4 signaling in pulmonary hypertension., Methods and Results: Proteomic analysis of CLIC4-interacting proteins in human pulmonary artery endothelial cells identified regulators of endosomal trafficking, including Arf6 (ADP ribosylation factor 6) GTPase activating proteins and clathrin, while CLIC4 overexpression affected protein regulators of vesicular trafficking, lysosomal function, and inflammation. CLIC4 reduced BMPRII (bone morphogenetic protein receptor II) expression and signaling as a result of Arf6-mediated reduction in gyrating clathrin and increased lysosomal targeting of the receptor. BMPRII expression was restored by Arf6 siRNA, Arf inhibitor Sec7 inhibitor H3 (SecinH3), and inhibitors of clathrin-mediated endocytosis but was unaffected by chloride channel inhibitor, indanyloxyacetic acid 94 or Arf1 siRNA. The effects of CLIC4 on NF-κB (nuclear factor-kappa B), HIF (hypoxia-inducible factor), and angiogenic response were prevented by Arf6 siRNA and SecinH3. Sugen/hypoxia mice and monocrotaline rats showed elevated expression of CLIC4, activation of Arf6 and NF-κB, and reduced expression of BMPRII in the lung. These changes were established early during disease development. Lung endothelium-targeted delivery of CLIC4 siRNA or treatment with SecinH3 attenuated the disease, reduced CLIC4/Arf activation, and restored BMPRII expression in the lung. Endothelial colony-forming cells from idiopathic pulmonary hypertensive patients showed upregulation of CLIC4 expression and Arf6 activity, suggesting potential importance of this pathway in the human condition., Conclusions: Arf6 is a novel effector of CLIC4 and a new therapeutic target in pulmonary hypertension.

Invasive Streptococcus pyogenes infections are rare, with often-unexplained severity. Prompt diagnosis is desirable, as deaths can occur rapidly following onset and there is an increased, but preventable, risk to contacts. Here, proteomic analyses of clinical samples from invasive human S. pyogenes infections were undertaken to determine if novel diagnostic targets could be detected, and to augment our understanding of disease pathogenesis. Fluid samples from 17 patients with confirmed invasive S. pyogenes infection (empyema, septic arthritis, necrotising fasciitis) were analysed by proteomics for streptococcal and human proteins; 16/17 samples had detectable S. pyogenes DNA. Nineteen unique S. pyogenes proteins were identified in just 6/17 samples, and 15 of these were found in a single pleural fluid sample including streptococcal inhibitor of complement, trigger factor, and phosphoglycerate kinase. In contrast, 469 human proteins were detected in patient fluids, 177 (38%) of which could be identified as neutrophil proteins, including alpha enolase and lactotransferrin which, together, were found in all 17 samples. Our data suggest that streptococcal proteins are difficult to detect in infected fluid samples. A vast array of human proteins associated with leukocyte activity are, however, present in samples that deserve further evaluation as potential biomarkers of infection.

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., (Copyright © 2017 Elsevier B.V. All rights reserved.)

Lipid composition determines membrane properties, and cholesterol plays a major role in this determination as it regulates membrane fluidity and permeability, as well as induces the formation of coexisting phases and domains in the membrane. Biological membranes display a very diverse lipid composition, the lateral organization of which plays a crucial role in regulating a variety of membrane functions. We hypothesize that, during biological evolution, membranes with a particular cholesterol content were selected to perform certain functions in the cells of eukaryotic organisms. In this review, we discuss the major membrane properties induced by cholesterol, and their relationship to certain membrane functions.

Tissue factor pathway inhibitor (TFPI), the main inhibitor of initiation of coagulation, exerts an important anticoagulant role through the factor Xa (FXa)-dependent inhibition of tissue factor/factor VIIa. Protein S is a TFPI cofactor, enhancing the efficiency of FXa inhibition. TFPI can also inhibit prothrombinase assembly by directly interacting with coagulation factor V (FV), which has been activated by FXa. Because full-length TFPI associates with FV in plasma, we hypothesized that FV may influence TFPI inhibitory function. Using pure component FXa inhibition assays, we found that although FV alone did not influence TFPI-mediated FXa inhibition, it further enhanced TFPI in the presence of protein S, resulting in an ∼8-fold reduction in K i compared with TFPI alone. A FV variant (R709Q/R1018Q/R1545Q, FV ΔIIa ) that cannot be cleaved/activated by thrombin or FXa also enhanced TFPI-mediated inhibition of FXa ∼12-fold in the presence of protein S. In contrast, neither activated FV nor recombinant B-domain-deleted FV could enhance TFPI-mediated inhibition of FXa in the presence of protein S, suggesting a functional contribution of the B domain. Using TFPI and protein S variants, we show further that the enhancement of TFPI-mediated FXa inhibition by protein S and FV depends on a direct protein S/TFPI interaction and that the TFPI C-terminal tail is not essential for this enhancement. In FXa-catalyzed prothrombin activation assays, both FV and FV ΔIIa (but not activated FV) enhanced TFPI function in the presence of protein S. These results demonstrate a new anticoagulant (cofactor) function of FV that targets the early phase of coagulation before prothrombinase assembly., (© 2017 by The American Society for Biochemistry and Molecular Biology, Inc.)

Xanthones are tricyclic compounds of natural or synthetic origin exhibiting a broad spectrum of therapeutic activities. Three synthetic xanthone derivatives (KS1, KS2, and KS3) with properties typical for nonsteroidal anti-inflammatory drugs (NSAID) were objects of the presented model study. NSAIDs are in common use however; several of them exhibit gastric toxicity predominantly resulting from their direct interactions with the outermost lipid layer of the gastric mucosa that impair its hydrophobic barrier property. Among the studied xanthones, gastric toxicity of only KS2 has been determined in previous pharmacological studies, and it is low. In this study, carried out using X-ray diffraction and computer simulation, a palmitoyloleoylphosphatidylcholine-cholesterol bilayer (POPC-Chol) was used as a model of a hydrophobic layer of lipids protecting gastric mucosa as POPC and Chol are the main lipids in human mucus. X-ray diffraction data were used to validate the computer model. The aim of the study was to assess potential gastric toxicity of the xanthones by analysing their atomic level interactions with lipids, ions, and water in the lipid bilayer and their effect on the bilayer physicochemical properties. The results show that xanthones have small effect on the bilayer properties except for its rigidity whereas their interactions with water, ions, and lipids depend on their protonation state and for a given state, are similar for all the xanthones. As gastric toxicity of KS2 is low, based on MD simulations one can predict that toxicity of KS1 and KS3 is also low., (Copyright © 2016 Elsevier B.V. All rights reserved.)

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., (Copyright © 2016 Elsevier B.V. All rights reserved.)

The range of exoproteins and core exoproteome of 14 Staphylococcus aureus isolates representing major lineages associated with asymptomatic carriage and clinical disease in the UK was identified by MS proteomics using a combined database incorporating sequences derived from 39 S. aureus genomes. In all, 632 different proteins were identified and, of these, only 52 (8 %) were found in all 14 isolates whereas 144 (23 %) were found in just a single isolate. Comparison of the observed mass of each protein (based on migration by SDS-PAGE) with its predicted mass (based on amino acid sequence) suggested that 95 % of the proteins identified were not subject to any major post-translational modification. Migration of 5 % of the proteins was not as expected: 1 % of the proteins migrated at a mass greater than predicted, while 4 % appeared to have undergone proteolytic cleavage; these included SsaA2, Aur, SspP, Ebh as well as BlaR1, MecR1, FsH, OatA and LtaS. Intriguingly, a truncated SasG was produced by a single CC8 USA300-like strain. The analysis provided evidence of the marked heterogeneity in protein expression by S. aureus in broth, while yielding a core but narrow common exoproteome.

Improvements in crystallographic hardware and software have allowed automated structure-solution pipelines to approach a near-`one-click' experience for the initial determination of macromolecular structures. However, in many cases the resulting initial model requires a laborious, iterative process of refinement and validation. A new method has been developed for the automatic modeling of side-chain conformations that takes advantage of rotamer-prediction methods in a crystallographic context. The algorithm, which is based on deterministic dead-end elimination (DEE) theory, uses new dense conformer libraries and a hybrid energy function derived from experimental data and prior information about rotamer frequencies to find the optimal conformation of each side chain. In contrast to existing methods, which incorporate the electron-density term into protein-modeling frameworks, the proposed algorithm is designed to take advantage of the highly discriminatory nature of electron-density maps. This method has been implemented in the program Fitmunk, which uses extensive conformational sampling. This improves the accuracy of the modeling and makes it a versatile tool for crystallographic model building, refinement and validation. Fitmunk was extensively tested on over 115 new structures, as well as a subset of 1100 structures from the PDB. It is demonstrated that the ability of Fitmunk to model more than 95% of side chains accurately is beneficial for improving the quality of crystallographic protein models, especially at medium and low resolutions. Fitmunk can be used for model validation of existing structures and as a tool to assess whether side chains are modeled optimally or could be better fitted into electron density. Fitmunk is available as a web service at http://kniahini.med.virginia.edu/fitmunk/server/ or at http://fitmunk.bitbucket.org/.

Trans unsaturated lipids are uncommon in nature. In the human diet, they occur as natural products of ruminal bacteria or from industrial food processing like hydrogenation of vegetable oils. Consumption of trans unsaturated lipids has been shown to have a negative influence on human health; in particular, the risk of cardiovascular disease is higher when the amount of trans unsaturated lipids in the diet is elevated. In this study, we first performed quantum mechanical calculations to specifically and accurately parameterize cis and trans mono-unsaturated lipids and subsequently validated the newly derived parameter set. Then, we carried out molecular dynamics (MD) simulations of lipid bilayers composed of cis or trans unsaturated lipids with and without cholesterol. Our results show that trans mono-unsaturated chains are more flexible than cis mono-unsaturated chains due to lower barriers for rotation around the single bonds next to the trans double bond than those next to the cis double bond. In effect, interactions between cholesterol and trans unsaturated chains are stronger than cis unsaturated chains, which results in a higher ordering effect of cholesterol in trans unsaturated bilayers., (Copyright © 2015. Published by Elsevier Ireland Ltd.)

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., (Copyright © 2015 Elsevier B.V. and Société Française de Biochimie et Biologie Moléculaire (SFBBM). All rights reserved.)

Immunity to common bacteria requires the generation of antibodies that promote opsonophagocytosis and neutralise toxins. Pooled human immunoglobulin is widely advocated as an adjunctive treatment for clinical Streptococcus pyogenes infection however, the protein targets of the reagent remain ill defined. Affinity purification of the anti-streptococcal antibodies present within pooled immunoglobulin resulted in the generation of an IgG preparation that promoted opsonophagocytic killing of S. pyogenes in vitro and provided passive immunity in vivo. Isolation of the streptococcal surface proteins recognised by pooled human immunoglobulin permitted identification and ranking of 94 protein antigens, ten of which were reproducibly identified across four contemporary invasive S. pyogenes serotypes (M1, M3, M12 and M89). The data provide novel insight into the action of pooled human immunoglobulin during invasive S. pyogenes infection, and demonstrate a potential route to enhance the efficacy of antibody based therapies.

In this data article we provide topologies and force field parameters files for molecular dynamics simulations of lipids in the OPLS-aa force field using the GROMACS package. This is the first systematic parameterization of lipid molecules in this force field. Topologies are provided for four phosphatidylcholines: saturated DPPC, mono-cis unsaturated POPC and DOPC, and mono-trans unsaturated PEPC. Parameterization of the phosphatidylcholines was achieved in two steps: first, we supplemented the OPLS force field parameters for DPPC with new parameters for torsion angles and van der Waals parameters for the carbon and hydrogen atoms in the acyl chains, as well as new partial atomic charges and parameters for torsion angles in the phosphatidylcholine and glycerol moieties [1]. Next, we derived parameters for the cis and trans double bonds and the neighboring them single bonds [2]. Additionally, we provide GROMACS input files with parameters describing simulation conditions (md.mdp), which are strongly recommended to be used with these lipids models. The data are associated with the research article "Cis and trans unsaturated phosphatidylcholine bilayers: a molecular dynamics simulation study" [2] and provided as supporting materials.

Lipid A is the most chemically invariant part of lipopolysaccharides (LPS). Both lipid A and LPS constitute the external layer of the outer membrane of Gram-negative bacteria. E. coli-specific hexacyl lipid A (ECLA) forms stable bilayers in the presence of sodium or magnesium cations. To characterize biologically relevant properties of the ECLA bilayer, and in particular its water/membrane interface, 800 ns molecular dynamics (MD) simulations of fully hydrated bilayers made of ECLA at 50 °C (i.e., 6 °C above the main phase transition) were performed. The validation of the computer model for the ECLA bilayer was performed using available experimental data. The overall good agreement with the data was found. An ECLA molecule makes on average ∼1.3 ion-mediated bridges with neighboring lipid molecules. The average number of interlipid hydrogen bonds is 2.7. The abundance of such intermolecular links results in tight packing of ECLA molecules in the bilayer and explains the relatively small value of the surface area per lipid (1.515 nm(2)).

Protein S functions as a cofactor for tissue factor pathway inhibitor (TFPI) and activated protein C (APC). The sex hormone binding globulin (SHBG)-like region of protein S, consisting of two laminin G-like domains (LG1 and LG2), contains the binding site for C4b-binding protein (C4BP) and TFPI. Furthermore, the LG-domains are essential for the TFPI-cofactor function and for expression of full APC-cofactor function. The aim of the current study was to localise functionally important interaction sites in the protein S LG-domains using amino acid substitutions. Four protein S variants were created in which clusters of surface-exposed amino acid residues within the LG-domains were substituted. All variants bound normally to C4BP and were fully functional as cofactors for APC in plasma and in pure component assays. Two variants, SHBG2 (E612A, I614A, F265A, V393A, H453A), involving residues from both LG-domains, and SHBG3 (K317A, I330A, V336A, D365A) where residues in LG1 were substituted, showed 50-60 % reduction in enhancement of TFPI in FXa inhibition assays. For SHBG3 the decreased TFPI cofactor function was confirmed in plasma based thrombin generation assays. Both SHBG variants bound to TFPI with decreased affinity in surface plasmon resonance experiments. The TFPI Kunitz 3 domain is known to contain the interaction site for protein S. Using in silico analysis and protein docking exercises, preliminary models of the protein S SHBG/TFPI Kunitz domain 3 complex were created. Based on a combination of experimental and in silico data we propose a binding site for TFPI on protein S, involving both LG-domains.

Molecular dynamics simulations of 1,2-di-O-acyl-3-O-β-D-galactopyranosyl-sn-glycerol (MGDG) and 1,2-dioleoyl-sn-glycero-3-phosphatidylcholine (DOPC) bilayers were carried out to compare the effect of the lipid head group's chemical structure on the dynamics and orientational order of the water molecules hydrating the bilayer. The effect of the bilayers on the diffusion of water is strong for the neighbouring water molecules i.e., those located not further than 4 Å from any bilayer atom. This is because the neighbouring water molecules are predominantly hydrogen bonded to the lipid oxygen atoms and their mobility is limited to a confined spatial volume. The choline group of DOPC and the galactose group of MGDG affect water diffusion less than the polar groups located deeper in the bilayer interface, and similarly. The latter is an unexpected result since interactions of water with these groups have a vastly different origin. The least affected by the bilayer lipids is the lateral diffusion of unbound water in the bilayer plane (x,y-plane) - it is because the diffusion is not confined by the periodic boundary conditions, whereas that perpendicular to the plane is. Interactions of water molecules with lipid groups also enforce certain orientations of water dipole moments. The profile of an average water orientation along the bilayer normal for the MGDG bilayer differs from that for the DOPC bilayer. In the DOPC bilayer, the ordering effect of the lipid head groups extends further into the water phase than in the MGDG bilayer, whereas inside the bilayer/water interface, ordering of the water dipoles in the MGDG bilayer is higher. It is possible that differences in the profiles of an average water orientation across the bilayer in the DOPC and MGDG bilayers are responsible for differences in the lateral pressure profiles of these bilayers.

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).

OPLS All-Atom (OPLS/AA) is a generic all-atom force field which was fine-tuned to accurately reproduce condensed phase properties of organic liquids. Its application in modeling of lipid membranes is, however, limited mainly due to the inability to correctly describe phase behavior and organization of the hydrophobic core of the model lipid bilayers. Here we report new OPLS/AA parameters for n-pentadecane, methyl acetate, and dimethyl phosphate anion. For the new force field parameters, we show very good agreement between calculated and numerous reference data, including liquid density, enthalpy of vaporization, free energy of hydration, and selected transport properties. The new OPLS/AA parameters have been used in successful submicrosecond MD simulations of bilayers made of bacterial glycolipids whose results will be published elsewhere shortly.

A study into the effects of amorphous nano-SiO2 particles on A549 lung epithelial cells was undertaken using proteomics to understand the interactions that occur and the biological consequences of exposure of lung to nanoparticles. Suitable conditions for treatment, where A549 cells remained viable for the exposure period, were established by following changes in cell morphology, flow cytometry, and MTT reduction. Label-free proteomics was used to estimate the relative level of proteins from their component tryptic peptides detected by mass spectrometry. It was found that A549 cells tolerated treatment with 100 µg/ml nano-SiO2 in the presence of 1.25% serum for at least 4 h. After this time detrimental changes in cell morphology, flow cytometry, and MTT reduction were evident. Proteomics performed after 4 h indicated changes in the expression of 47 proteins. Most of the proteins affected fell into four functional groups, indicating that the most prominent cellular changes were those that affected apoptosis regulation (e.g. UCP2 and calpain-12), structural reorganisation and regulation of actin cytoskeleton (e.g. PHACTR1), the unfolded protein response (e.g. HSP 90), and proteins involved in protein synthesis (e.g. ribosomal proteins). Treatment with just 10 µg/ml nano-SiO2 particles in serum-free medium resulted in a rapid deterioration of the cells and in medium containing 10% serum the cells were resistant to up to 1000 µg/ml nano-SiO2 particles, suggesting interaction of serum components with the nanoparticles. A variety of serum proteins were found which bound to nano-SiO2 particles, the most prominent of which were albumin, apolipoprotein A-I, hemoglobin, vitronectin and fibronectin. The use of a proteomics platform, with appropriately designed experimental conditions, enabled the early biological perturbations induced by nano-SiO2 in a model target cell system to be identified. The approach facilitates the design of more focused test systems for use in tiered evaluations of nanomaterials.

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.

Group A streptococcal isolates of serotype M18 are historically associated with epidemic waves of pharyngitis and the non-suppurative immune sequela rheumatic fever. The serotype is defined by a unique, highly encapsulated phenotype, yet the molecular basis for this unusual colony morphology is unknown. Here we identify a truncation in the regulatory protein RocA, unique to and conserved within our serotype M18 GAS collection, and demonstrate that it underlies the characteristic M18 capsule phenotype. Reciprocal allelic exchange mutagenesis of rocA between M18 GAS and M89 GAS demonstrated that truncation of RocA was both necessary and sufficient for hyper-encapsulation via up-regulation of both precursors required for hyaluronic acid synthesis. Although RocA was shown to positively enhance covR transcription, quantitative proteomics revealed RocA to be a metabolic regulator with activity beyond the CovR/S regulon. M18 GAS demonstrated a uniquely protuberant chain formation following culture on agar that was dependent on excess capsule and the RocA mutation. Correction of the M18 rocA mutation reduced GAS survival in human blood, and in vivo naso-pharyngeal carriage longevity in a murine model, with an associated drop in bacterial airborne transmission during infection. In summary, a naturally occurring truncation in a regulator explains the encapsulation phenotype, carriage longevity and transmissibility of M18 GAS, highlighting the close interrelation of metabolism, capsule and virulence.

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., (© 2012 The Authors. Journal of Neurochemistry © 2012 International Society for Neurochemistry.)

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., (Copyright © 2011. Published by Elsevier B.V.)

Background: Human skin has the capacity to metabolise foreign chemicals (xenobiotics), but knowledge of the various enzymes involved is incomplete. A broad-based unbiased proteomics approach was used to describe the profile of xenobiotic metabolising enzymes present in human skin and hence indicate principal routes of metabolism of xenobiotic compounds. Several in vitro models of human skin have been developed for the purpose of safety assessment of chemicals. The suitability of these epidermal models for studies involving biotransformation was assessed by comparing their profiles of xenobiotic metabolising enzymes with those of human skin., Methodology/principal Findings: Label-free proteomic analysis of whole human skin (10 donors) was applied and analysed using custom-built PROTSIFT software. The results showed the presence of enzymes with a capacity for the metabolism of alcohols through dehydrogenation, aldehydes through dehydrogenation and oxidation, amines through oxidation, carbonyls through reduction, epoxides and carboxylesters through hydrolysis and, of many compounds, by conjugation to glutathione. Whereas protein levels of these enzymes in skin were mostly just 4-10 fold lower than those in liver and sufficient to support metabolism, the levels of cytochrome P450 enzymes were at least 300-fold lower indicating they play no significant role. Four epidermal models of human skin had profiles very similar to one another and these overlapped substantially with that of whole skin., Conclusions/significance: The proteomics profiling approach was successful in producing a comprehensive analysis of the biotransformation characteristics of whole human skin and various in vitro skin models. The results show that skin contains a range of defined enzymes capable of metabolising different classes of chemicals. The degree of similarity of the profiles of the in vitro models indicates their suitability for epidermal toxicity testing. Overall, these results provide a rational basis for explaining the fate of xenobiotics in skin and will aid chemical safety testing programmes.

Lutein is present in the human retina and lens, where it plays a protective role. As lutein is associated with the lipid matrix of biomembranes, the role depends on its membrane location. Experimental studies predicted two orientations of lutein in a phosphatidylcholine (PC) bilayer: vertical and horizontal. Using a molecular dynamics simulation, we observed, in two different PC bilayers, both orientations of lutein, and in each bilayer, a single change from vertical to horizontal orientation or vice versa. Both orientations were stabilized by hydrogen bonding of lutein OH groups with mainly carbonyl but also phosphate oxygen atoms of PC.

We used atomistic simulations to study the membrane-bound form of catechol-O-methyltransferase (MB-COMT). In particular we investigated the 26-residue transmembrane α-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 ARG27 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., (© 2011 American Chemical Society)

This communication presents a molecular dynamics simulation study of a bilayer consisting of 128 dioleoyl-sn-glycero-3-phosphocholine molecules, which focusses on the center-of-mass diffusion of the lipid molecules parallel to the membrane plane. The analysis of the simulation results is performed within the framework of the generalized Langevin equation and leads to a consistent picture of subdiffusion. The mean square displacement of the lipid molecules evolves as ∝ t(α), with α between 0.5 and 0.6, and the fractional diffusion coefficient is close to the experimental value for a similar system obtained by fluorescence correlation spectroscopy. We show that the long-time tails of the lateral velocity autocorrelation function and the associated memory function agree well with exact results which have been recently derived by asymptotic analysis [G. Kneller, J. Chem. Phys. 134, 224106 (2011)]. In this context, we define characteristic time scales for these two quantities., (© 2011 American Institute of Physics)

Objectives: Accurate differentiation between benign and malignant causes of biliary obstruction remains challenging and reliable biomarkers are urgently needed. Bile is a potential source of such biomarkers. Our aim was to apply a proteomic approach to identify a potential biomarker in bile that differentiates between malignant and benign disease, and to assess its diagnostic accuracy. Neutrophil gelatinase-associated lipocalin (NGAL) is multi-functional protein, released from activated neutrophils, with roles in inflammation, immune function, and carcinogenesis. It has not previously been described in bile., Methods: Bile, urine, and serum were collected prospectively from 38 patients undergoing endoscopic retrograde cholangiopancreatography ("discovery" cohort); 22 had benign and 16 had malignant pancreatobiliary disease. Initially, label-free proteomics and immunoblotting were performed in samples from a subset of these patients. Enzyme-linked immunosorbent assay was then performed for NGAL as a potential biomarker on all samples in this cohort. The diagnostic performance of biliary NGAL was then validated in a second, independent group ("validation" cohort) of 21 patients with pancreatobiliary disease (benign n=14, malignant n=7)., Results: NGAL levels were significantly raised in bile from the malignant disease group, compared with bile from the benign disease group in the discovery cohort (median 1,556 vs. 480 ng/ml, P=0.007). Biliary NGAL levels had a receiver operating characteristic area under curve of 0.76, sensitivity 94%, specificity 55%, positive predictive value 60%, and negative predictive value 92% for distinguishing malignant from benign causes. Biliary NGAL was independent of serum biochemistry and carbohydrate antigen 19-9 (CA 19-9) in differentiating between underlying benign and malignant disease. No significant differences in serum and urine NGAL levels were found between benign and malignant disease. Combining biliary NGAL and serum CA 19-9 improved diagnostic accuracy for malignancy (sensitivity 85%, specificity 82%, positive predictive value 79%, and negative predictive value 87%). The diagnostic accuracy of biliary NGAL was confirmed in the second independent validation cohort., Conclusions: NGAL in bile is a novel potential biomarker to help distinguish benign from malignant biliary obstruction.

We have combined Langmuir monolayer film experiments and all-atom molecular dynamics (MD) simulation of a bilayer to study the surface structure of a PEGylated liposome and its interaction with the ionic environment present under physiological conditions. Lipids that form both gel and liquid-crystalline membranes have been used in our study. By varying the salt concentration in the Langmuir film experiment and including salt at the physiological level in the simulation, we have studied the effect of salt ions present in the blood plasma on the structure of the poly(ethylene glycol) (PEG) layer. We have also studied the interaction between the PEG layer and the lipid bilayer in both the liquid-crystalline and gel states. The MD simulation shows two clear results: (a) The Na(+) ions form close interactions with the PEG oxygens, with the PEG chains forming loops around them and (b) PEG penetrates the lipid core of the membrane for the case of a liquid-crystalline membrane but is excluded from the tighter structure of the gel membrane. The Langmuir monolayer results indicate that the salt concentration affects the PEGylated lipid system, and these results can be interpreted in a fashion that is in agreement with the results of our MD simulation. We conclude that the currently accepted picture of the PEG surface layer acting as a generic neutral hydrophilic polymer entirely outside the membrane, with its effect explained through steric interactions, is not sufficient. The phenomena we have observed may affect both the interaction between the liposome and bloodstream proteins and the liquid-crystalline-gel transition and is thus relevant to nanotechnological drug delivery device design., (© 2011 American Chemical Society)

A high percentage of people treated with a long-term nonsteroidal anti-inflammatory drug (NSAID) therapy suffer NSAID-induced gastrointestinal-tract-related side effects. A current hypothesis states that the side effects are related to the topical action of NSAID molecules on gastric mucus that lowers its resistance to luminal acid. The main lipids in human mucus are palmitoyloleoylphosphatidylcholine (POPC) and cholesterol (Chol). In this study, both X-ray diffraction and molecular dynamics (MD) simulation methods were employed to investigate the effects of selected NSAIDs in protonated and deprotonated states on the structural parameters of a POPC-Chol bilayer. The drugs were three commonly used NSAIDs with apparently different gastric toxicity: ketoprofen (KET), aspirin (ASP), and piroxicam (PXM). Both methods revealed that the effects of the NSAIDs on the POPC-Chol bilayer parameters were moderate and only slightly differentiated among the drugs. Much larger differences among the drugs were noticed in their interactions with interfacial water and Na(+) as well as with the polar groups of POPC and Chol, mainly via H-bonds. Of the three NSAIDs, KET interacted with POPC and water the most extensively, whereas ASP interacted with Chol and Na(+) more than did the other two. Interactions of PXM with POPC and Chol polar groups as well as with water and Na(+) were limited.

We performed 200 ns MD simulations of phosphatidylcholine (PC) bilayers in the liquid crystalline (L(α)) and gel (L(β)') states to compare the properties of the water-membrane interfaces in these two thermotropic bilayer phases. Our simulations show that the membrane phase determines the behavior of water, ions, and PC head groups. When the membrane was in the gel phase, we found partial dehydration (fewer PC-water interactions), particularly in the carbonyl groups region, as well as an almost complete lack of ionic penetration into this region as compared with the bilayer in the liquid-crystalline phase. In the latter case, there is an exchange of Na(+) ions between the water layer and the interfacial region. This is mainly due to the fact that the most stable binding of Na(+) in the liquid-crystalline bilayer is to the carbonyl groups. The lack of Na(+) binding to the carbonyl groups in the gel phase bilayer can be explained by the more compact structure of the bilayer in this phase.

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.

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 beta-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.

A molecular dynamics simulation study of four lipid bilayers with inserted trans-membrane helical fragment of epithelial growth factor (EGF) receptor (EGF peptide) was performed. The lipid bilayers differ in their lipid composition and consist of (i) unsaturated phosphatidylcholine (palmitoyloleoylphosphatidylcholine, POPC), (ii) POPC and 20 mol% of cholesterol (Chol), (iii) sphingomyelin (SM) and 20 mol% of Chol, and (iv) SM and 50 mol% of Chol. Only 1 out of 26 residues in the EGF-peptide sequence is polar (Thr). The hydrophobic thickness of each bilayer is different but shorter than the length of the peptide and so, due to hydrophobic mismatch, the inserted peptide is tilted in each bilayer. Additionally, in the POPC bilayer, which is the thinnest, the peptide loses its helical structure in a short three-amino acid fragment. This facilitates bending of the peptide and burying all hydrophobic amino acids inside the membrane core (Figure 1(b)). Bilayer lipid composition affects interactions between the peptide and lipids in the membrane core. Chol increases packing of atoms relative to the peptide side chains, and thus increases van der Waals interactions. On average, the packing around the peptide is higher in SM-based bilayers than POPC-based bilayers but for certain amino acids, packing depends on their position relative to the bilayer center. In the bilayer center, packing is higher in POPC-based bilayers, while in regions closer to the interface packing is higher in SM-based bilayers. In general, amino acids with larger side chains interact strongly with lipids, and thus the peptide sequence is important for the pattern of interactions at different membrane depths. This pattern closely resembles the shape of recently published lateral pressure profiles [Ollila et alJ. Struct. Biol. DOI:10.1016/j.jsb.2007.01.012].

The 3alpha-hydroxyl group is a characteristic structural element of all membrane sterol molecules, while the 3-ketone group is more typically found in steroid hormones. In this work, we investigate the effect of substituting the hydroxyl group in cholesterol with the ketone group to produce ketosterone. Extensive atomistic molecular dynamics simulations of saturated lipid membranes with either cholesterol or ketosterone show that, like cholesterol, ketosterone increases membrane order and induces condensation. However, the effect of ketosterone on membrane properties is considerably weaker than that of cholesterol. This is largely due to the unstable positioning of ketosterone at the membrane-water interface, which gives rise to a small but significant number of flip-flop transitions, where ketosterone is exchanged between membrane leaflets. This is remarkable, as flip-flop motions of sterol molecules have not been previously reported in analogous lipid bilayer simulations. In the same context, ketosterone is found to be more tilted with respect to the membrane normal than cholesterol. The atomic level mechanism responsible for the increase of the steroid tilt and the promotion of flip-flops is the decrease in polar interactions at the membrane-water interface. Interactions between lipids or water and the ketone group are found to be weaker than in the case of the hydroxyl group, which allows ketosterone to penetrate through the hydrocarbon region of a membrane.

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.

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.

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 molecular dynamics (MD) simulation of the fully hydrated bilayer made of 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) and containing beta-carotene (beta-Car) molecules was carried out as a complementary approach to experimental techniques to investigate the orientation of beta-Car in the lipid membrane as well as its influence on the bilayer properties. The bilayer reached thermal equilibrium after 1200 ps of MD simulation and the productive run was carried out for 2800 ps. The results indicate that the carotene rings are located in the region occupied by the carbonyl groups of the POPC gamma-chain with no trace of penetration towards the centre of the bilayer. Carotene exhibits an ordering effect on both the beta- and the gamma-chain. While the fully saturated gamma-chain is affected evenly along, the order of the mono-unsaturated beta-chain is modified mainly below the double bond. In general, a high value of the order parameter and the chain tilt in the range from 11.4 degrees to 26.7 degrees were observed for the beta-Car molecules. However, for chain segment adjacent to methyl groups the value of the order parameter is low and the tilt angle is close to 75 degrees . Moreover, the probability of trans conformation being generally close to 1.0 along the beta-Car chain is reduced for these segments. Our MD simulation study suggests two pools of the preferential orientation of beta-Car: a slightly bent structure corresponding to a small chain tilt angle and a rather stretched structure that corresponds to a higher chain tilt. The results are discussed in the light of experimental findings.

Phosphatidylethanolamine (PE) and phosphatidylglycerol (PG) are the main lipid components of the inner bacterial membrane. A computer model for such a membrane was built of palmitoyloleoyl PE (POPE) and palmitoyloleoyl PG (POPG) in the proportion 3:1, and sodium ions (Na+) to neutralize the net negative charge on each POPG (POPE-POPG bilayer). The bilayer was simulated for 25 ns. A final 10-ns trajectory fragment was used for analyses. In the bilayer interfacial region, POPEs and POPGs interact readily with one another via intermolecular hydrogen (H) bonds and water bridges. POPE is the main H-bond donor in either PEPE or PEPG H-bonds; PGPG H-bonds are rarely formed. Almost all POPEs are H-bonded and/or water bridged to either POPE or POPG but PE-PG links are favored. In effect, the atom packing in the near-the-interface regions of the bilayer core is tight. Na+ does not bind readily to lipids, and interlipid links via Na+ are not numerous. Although POPG and POPE comprise one bilayer, their bilayer properties differ. The average surface area per POPG is larger and the average vertical location of the POPG phosphate group is lower than those of POPE. Also, the alkyl chains of POPG are more ordered and less densely packed than the POPE chains. The main conclusion of this study is that in the PE-PG bilayer PE interacts more strongly with PG than with PE. This is a likely molecular-level event behind a regulating mechanism developed by the bacteria to control its membrane permeability and stability consisting in changes of the relative PG/PE concentration in the membrane.

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.

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.

G protein-coupled receptors relay diverse extracellular signals into cells via a common mechanism, involving activation of cytosol G proteins. The mechanism underlies the actions of approximately 50% of all drugs. In this work, we focus on simulating three protein-ligand complexes of the neurohypophyseal hormone analog 4-OH-phenylacetyl- D-Y(Me)FQNRPR-NH2 (I) with the human V1a, V2 and oxytocin receptors. The peptideI is a potent selective V1a receptor antagonist. To obtain relaxed models of the complexes, the following techniques were used: docking ofI into the vasopressin V1a, V2 and oxytocin receptor models, optimization of the geometry of the resulting complexes and molecular dynamics in a fully hydrated 1-palmitoyl-2-oleoyl- sn-glycero-3-phosphatidylcholine lipid bilayer. The results of the simulations allow us to draw some conclusions about the ligand selectivity to V1aR.