Search

Graphical abstract, Highlights • Cholesterol is an indispensable component of mammalian cell membranes. • In water, cholesterol molecules self-assemble. • Cholesterol self-assembly outside the membrane may promote pathologic processes. • Chirality of oxysterol affects its association with cholesterol in water. • Energies of sterol hydration and association in water were calculated and compared., Cholesterol (Chol) is the most prevalent sterol in the animal kingdom and an indispensable component of mammalian cell membranes. Chol content in the membrane is strictly controlled, although the oxidation of phospholipids may change the relative content of membrane Chol. An excess of it results in the formation of pure Chol microdomains in the membrane. It is likely that some Chol molecules detach from the domains and self-assemble in the aqueous environment. This may promote Chol microcrystallisation, which initiates the development of gallstones and atherosclerotic plaque. In this study, the molecular dynamics, free energy perturbation, umbrella sampling and Voronoi diagram methods are used to reveal the details of self-association of Chol and its oxidised forms (oxChol), namely 7α,β-hydroxycholesterol and 7α,β-hydroperoxycholesterol, in water. In the first part of the study the interactions between a sterol monomer and water over a short and longer timescale as well as the energy of hydration of each sterol are analysed. This helps one to understand Chol-Chol and Chol-OxChol with different chirality self-association in water better, which is analysed in the second part of the study. The Voronoi diagram approach is used to determine the relative arrangement of molecules in the dimer and, most importantly, to analyse the dehydration of the contacting surfaces of the assembling molecules. Free energy calculations indicate that Chol and 7β-hydroxycholesterol associate into the most stable dimer and that Chol-Chol is the next most stable of the five dimers studied. Employing different computational methods enables us to obtain an adequate picture of Chol-sterol self-association in water, which includes dynamic, energetic and temporal aspects of the process.

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

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

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.

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

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. Keywords: Xanthophyll, Energy barrier, Hydrogen bond, Hydrophobic effect, Molecular dynamics, Umbrella sampling

Lipid and water molecules comprising the bilayer form an integral entity owing to only weak physical interactions. At the bilayer interface, these interactions chiefly involve hydrogen bonding and charge pairing. Lipid head groups make hydrogen bonds (H-bonds) predominantly with water, whereas interlipid H-bonds and charge pairs are less numerous. Both interlipid H-bonding and charge pairing depend on the distance and relative orientation of the interacting head groups. In this computational paper, correlations are analysed between the orientation of the lipid head group and the number of interlipid interactions at the interface of a bilayer made of galactolipids, forming direct interlipid H-bonds, and of phosphatidylcholines forming interlipid charge pairs. The correlations are not strong, however, in both bilayers they show a similar trend.

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.

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.

Chloroplasts were formed by uptake of cyanobacteria into eukaryotic cells ca. 1.6 billion years ago. During evolution most of the cyanobacterial genes were transferred from the chloroplast to the nuclear genome. The rbcX gene, encoding an assembly chaperone required for Rubisco biosynthesis in cyanobacteria, was duplicated. Here we demonstrate that homologous eukaryotic chaperones (AtRbcX1 and AtRbcX2) demonstrate different affinities for the C-terminus of Rubisco large subunit and determine their crystal structures.Three-dimensional structures of AtRbcX1 and AtRbcX2 were resolved by the molecular replacement method. Equilibrium binding constants of the C-terminal RbcL peptide by AtRbcX proteins were determined by spectrofluorimetric titration. The binding mode of RbcX-RbcL was predicted using molecular dynamic simulation.We provide crystal structures of both chaperones from Arabidopsis thaliana providing the first structural insight into Rubisco assembly chaperones form higher plants. Despite the low sequence homology of eukaryotic and cyanobacterial Rubisco chaperones the eukaryotic counterparts exhibit surprisingly high similarity of the overall fold to previously determined prokaryotic structures. Modeling studies demonstrate that the overall mode of the binding of RbcL peptide is conserved among these proteins. As such, the evolution of RbcX chaperones is another example of maintaining conserved structural features despite significant drift in the primary amino acid sequence.The presented results are the approach to elucidate the role of RbcX proteins in Rubisco assembly in higher plants.

Staphylococcus aureus is a dangerous human pathogen. A number of the proteins secreted by this bacterium are implicated in its virulence, but many of the components of its secretome are poorly characterized. Strains of S. aureus can produce up to six homologous extracellular serine proteases grouped in a single spl operon. Although the SplA, SplB, and SplC proteases have been thoroughly characterized, the properties of the other three enzymes have not yet been investigated. Here, we describe the biochemical and structural characteristics of the SplD protease. The active enzyme was produced in an Escherichia coli recombinant system and purified to homogeneity. P1 substrate specificity was determined using a combinatorial library of synthetic peptide substrates showing exclusive preference for threonine, serine, leucine, isoleucine, alanine, and valine. To further determine the specificity of SplD, we used high-throughput synthetic peptide and cell surface protein display methods. The results not only confirmed SplD preference for a P1 residue, but also provided insight into the specificity of individual primed- and non-primed substrate-binding subsites. The analyses revealed a surprisingly narrow specificity of the protease, which recognized five consecutive residues (P4-P3-P2-P1-P1’) with a consensus motif of R-(Y/W)-(P/L)-(T/L/I/V)↓S. To understand the molecular basis of the strict substrate specificity, we crystallized the enzyme in two different conditions, and refined the structures at resolutions of 1.56 Å and 2.1 Å. Molecular modeling and mutagenesis studies allowed us to define a consensus model of substrate binding, and illustrated the molecular mechanism of protease specificity.