1
|
Kumar A, Daschakraborty S. Anomalous lateral diffusion of lipids during the fluid/gel phase transition of a lipid membrane. Phys Chem Chem Phys 2023; 25:31431-31443. [PMID: 37962400 DOI: 10.1039/d3cp04081j] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2023]
Abstract
A lipid membrane undergoes a phase transition from fluid to gel phase upon changing external thermodynamic conditions, such as decreasing temperature and increasing pressure. Extremophilic organisms face the challenge of preventing this deleterious phase transition. The main focus of their adaptive strategy is to facilitate effective temperature sensing through sensor proteins, relying on the drastic changes in packing density and membrane fluidity during the phase transition. Although the changes in packing density parameters due to the fluid/gel phase transition are studied in detail, the impact on membrane fluidity is less explored in the literature. Understanding the lateral diffusive dynamics of lipids in response to temperature, particularly during the fluid/gel phase transition, is albeit crucial. Here we have simulated the phase transition of a single component lipid membrane composed of dipalmitoylphosphatidylcholine (DPPC) lipids using a coarse-grained (CG) model and studied the changes of the structural and dynamical properties. It is observed that near the phase transition point, both fluid and gel phase domains coexist together. The dynamics remains highly non-Gaussian for a long time even when the mean square displacement reaches the Fickian regime at a much earlier time. This Fickian yet non-Gaussian diffusion (FnGD) is a characteristic of a highly heterogeneous system, previously observed for the lateral diffusion of lipids in raft mimetic membranes having liquid-ordered and liquid-disordered phases co-existing together. We have analyzed the molecular trajectories and calculated the jump-diffusion of the lipids, stemming from sudden jump translations, using a translational jump-diffusion (TJD) approach. An overwhelming contribution of the jump-diffusion of the lipids is observed suggesting anomalous diffusion of lipids during fluid/gel phase transition of the membrane. These results are important in unravelling the intricate nature of lipid diffusion during the phase transition of the membrane and open up a new possibility of investigating the most significant change of membrane properties during phase transition, which can be effectively sensed by proteins.
Collapse
Affiliation(s)
- Abhay Kumar
- Department of Chemistry, Indian Institute of Technology Patna, Bihar 801106, India.
| | | |
Collapse
|
2
|
Erimban S, Daschakraborty S. Homeoviscous Adaptation of the Lipid Membrane of a Soil Bacterium Surviving under Diurnal Temperature Variation: A Molecular Simulation Perspective. J Phys Chem B 2022; 126:7638-7650. [PMID: 36166758 DOI: 10.1021/acs.jpcb.2c01359] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A recent experiment has reported the lipidome remodeling of a soil-based plant-associated bacterium Methylobacterium extorquens due to diurnal temperature variations. The key adaptation strategy is the headgroup-specific remodeling of the acyl chain. To understand the idiosyncratic adaptation at the molecular level, we simulate the model membrane of the same bacterium using the reported lipidome compositions at four different experimental temperatures. We investigate the temperature-dependent packing density and fluidity of the membrane, the constancy of which is key to the homeoviscous adaptation. The results show that complex lipidome remodeling approximately preserves membrane properties under heat and cold stress. The headgroup-specific remodeling of the acyl chain serves to fine-tune the packing density and fluidity of the membrane at different temperatures. While lipids with strongly interacting headgroups are more abundant at higher temperatures, the lipidome is more dominated by lipids with weaker interacting headgroups at lower temperatures. This adaptation alleviates lipid membrane disruption caused by heat and cold stress. This study provides a molecular picture of the homeoviscous adaptation of the realistic lipid membrane of a soil-based bacterium.
Collapse
Affiliation(s)
- Shakkira Erimban
- Department of Chemistry, Indian Institute of Technology, Patna, Bihar 801106, India
| | | |
Collapse
|
3
|
Piskulich ZA, Cui Q. Machine Learning-Assisted Phase Transition Temperatures from Generalized Replica Exchange Simulations of Dry Martini Lipid Bilayers. J Phys Chem Lett 2022; 13:6481-6486. [PMID: 35819105 DOI: 10.1021/acs.jpclett.2c01654] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Accurate estimation of phase transition temperatures has been a longstanding challenge for molecular simulations. Recently, the generalized Replica Exchange technique for estimating phase transition temperatures has allowed for improved sampling of the phase transition; however, it requires a significant number of simultaneous replicas both inside and outside of the transition region leading to costly computational expense. In this work, the recently developed machine learning-assisted lipid phase analysis technique for learning the phase of individual lipids has been combined with generalized Replica Exchange Molecular Dynamics to reduce the overall computational expense of evaluating transition temperatures. This technique is then applied to eight different Dry Martini lipids to demonstrate its ability to describe transition temperatures as a function of chain length and tail saturation.
Collapse
Affiliation(s)
- Zeke A Piskulich
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, United States
| | - Qiang Cui
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, United States
| |
Collapse
|
4
|
Hudiyanti D, Christa SM, Mardhiyyah NH, Anugrah DSB, Widiarih T, Siahaan P. Dynamics insights into aggregation of phospholipid species with cholesterol and vitamin C. PHARMACIA 2022. [DOI: 10.3897/pharmacia.69.e81435] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
This paper provides dynamic insight into the aggregation profile of systems containing six different phospholipid species, cholesterol, and vitamin C thru Coarse-Grain Molecular Dynamics (CGMD) simulations. The simulation used 42 systems, and each system was composed of 220 molecules of each phospholipid species, a varied number of cholesterol molecules (0, 11, 22, 33, 66, 88), and 10 vitamin C molecules. The phospholipid species were DLPE, DOPE, DLiPE, DOPS, DLiPS, and DLiPC. We found curved bilayer, toroidal bilayer, concave micelle, disc-like bilayer, planar bilayer, and liposome structures in the systems during the 40 ns simulation. The systems with a ratio cholesterol:phospholipid between 15% and 40% formed liposomes regardless of the phospholipid species. Cholesterol is positioned in the liposome bilayer while vitamin C is encapsulated in the aqueous core of liposomes for all cholesterol compositions. The cholesterol influences the liposome formation of various phospholipid species and the encapsulation of vitamin C in the liposome structure.
Collapse
|
5
|
Cino EA, Tieleman DP. Curvature-based sorting of eight lipid types in asymmetric buckled plasma membrane models. Biophys J 2022; 121:2060-2068. [PMID: 35524412 DOI: 10.1016/j.bpj.2022.05.002] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2022] [Revised: 04/13/2022] [Accepted: 05/02/2022] [Indexed: 11/02/2022] Open
Abstract
Curvature is a fundamental property of biological membranes and has essential roles in cellular function. Bending of membranes can be induced by their lipid and protein compositions, as well as peripheral proteins, such as those that make up the cytoskeleton. An important aspect of membrane function is the grouping of lipid species into microdomains, or rafts, which serve as platforms for specific biochemical processes. The fluid mosaic model of membranes has evolved to recognize the importance of curvature and leaflet asymmetry, and there are efforts towards evaluating their functional roles. This work investigates the effect of curvature on the sorting of lipids in buckled asymmetric bilayers containing eight lipid types, approximating an average mammalian plasma membrane, through coarse-grained (CG) molecular dynamics (MD) simulations with the Martini force field. The simulations reveal that i) leaflet compositional asymmetry can induce curvature asymmetry, ii) lipids are sorted by curvature to different extents, and iii) curvature-based partitioning trends show moderate to strong correlations with lipid molecular volumes and head to tail bead ratios, respectively. The findings provide unique insights into the role of curvature in membrane organization, and the curvature-based sorting trends should be useful references for later investigations, and potentially interpreting the functional roles of specific lipids.
Collapse
Affiliation(s)
- Elio A Cino
- Centre for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - D Peter Tieleman
- Centre for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada.
| |
Collapse
|
6
|
Wilson KA, Wang L, O’Mara ML. Site of Cholesterol Oxidation Impacts Its Localization and Domain Formation in the Neuronal Plasma Membrane. ACS Chem Neurosci 2021; 12:3873-3884. [PMID: 34633798 DOI: 10.1021/acschemneuro.1c00395] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
Cholesterol is integral to the structure of mammalian cell membranes. Oxidation of cholesterol alters how it behaves in the membrane and influences the membrane biophysical properties. Elevated levels of oxidized cholesterol are associated with neurodegenerative diseases such as Alzheimer's disease, Parkinson's disease, multiple sclerosis, and Huntington's disease. Previous work has investigated the impact of oxidized cholesterol in the context of simple model membrane systems. However, there is a growing body of literature that shows that complex membranes possessing physiological phospholipid distributions have different properties from those of binary or trinary model membranes. In the current work, the impact of oxidized cholesterol on the biophysical properties of a complex neuronal plasma membrane is investigated using coarse-grained Martini molecular dynamics simulations. Comparison of the native neuronal membrane to neuronal membranes containing 10% tail-oxidized or 10% head-oxidized cholesterol shows that the site of oxidization changes the behavior of the oxidized cholesterol in the membrane. Furthermore, species-specific domain formation is observed between each oxidized cholesterol and minor lipid classes. Although both tail-oxidized and head-oxidized cholesterols modulate the biophysical properties of the membrane, smaller changes are observed in the complex neuronal membrane than seen in the previous work on simple binary or trinary model membranes. This work highlights the presence of compensatory effects of lipid diversity in the complex neuronal membrane. Overall, this study improves our molecular-level understanding of the effects of oxidized cholesterol on the properties of neuronal tissue and emphasizes the importance of studying membranes with realistic lipid compositions.
Collapse
Affiliation(s)
- Katie A. Wilson
- Research School of Chemistry, College of Science, The Australian National University, Canberra, Australian Capital Territory 2601, Australia
| | - Lily Wang
- Research School of Chemistry, College of Science, The Australian National University, Canberra, Australian Capital Territory 2601, Australia
| | - Megan L. O’Mara
- Research School of Chemistry, College of Science, The Australian National University, Canberra, Australian Capital Territory 2601, Australia
| |
Collapse
|
7
|
Kangarlou B, Dahanayake R, Martin IJ, Ndaya D, Wu CM, Kasi RM, Dormidontova EE, Nieh MP. Flower-like Micelles of Polyethylene Oxide End-Capped with Cholesterol. Macromolecules 2021. [DOI: 10.1021/acs.macromol.1c00896] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Behrad Kangarlou
- Materials Science Program, Institute of Materials Science, University of Connecticut, Storrs, Connecticut 06269, United States
| | - Rasika Dahanayake
- Department of Physics, University of Connecticut, Storrs, Connecticut 06269, United States
| | - Ian J. Martin
- Department of Chemistry, University of Connecticut, Storrs, Connecticut 06269, United States
| | - Dennis Ndaya
- Department of Chemistry, University of Connecticut, Storrs, Connecticut 06269, United States
| | - Chun-Ming Wu
- The Australian Nuclear Science and Technology Organisation, New Illawarra Road, Lucas Heights, New South Wales 2234, Australia
| | - Rajeswari M. Kasi
- Polymer Program, Institute of Materials Science, University of Connecticut, Storrs, Connecticut 06269, United States
- Department of Chemistry, University of Connecticut, Storrs, Connecticut 06269, United States
| | - Elena E. Dormidontova
- Polymer Program, Institute of Materials Science, University of Connecticut, Storrs, Connecticut 06269, United States
- Department of Physics, University of Connecticut, Storrs, Connecticut 06269, United States
| | - Mu-Ping Nieh
- Polymer Program, Institute of Materials Science, University of Connecticut, Storrs, Connecticut 06269, United States
- Department of Chemical and Biomolecular Engineering, University of Connecticut, Storrs, Connecticut 06269, United States
| |
Collapse
|
8
|
Franco-Ulloa S, Guarnieri D, Riccardi L, Pompa PP, De Vivo M. Association Mechanism of Peptide-Coated Metal Nanoparticles with Model Membranes: A Coarse-Grained Study. J Chem Theory Comput 2021; 17:4512-4523. [PMID: 34077229 PMCID: PMC8280734 DOI: 10.1021/acs.jctc.1c00127] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2021] [Indexed: 11/28/2022]
Abstract
Functionalized metal nanoparticles (NPs) hold great promise as innovative tools in nanomedicine. However, one of the main challenges is how to optimize their association with the cell membrane, which is critical for their effective delivery. Recent findings show high cellular uptake rates for NPs coated with the polycationic cell-penetrating peptide gH625-644 (gH), although the underlying internalization mechanism is poorly understood. Here, we use extended coarse-grained simulations and free energy calculations to study systems that simultaneously include metal NPs, peptides, lipids, and sterols. In particular, we investigate the first encounter between multicomponent model membranes and 2.5 nm metal NPs coated with gH (gHNPs), based on the evidence from scanning transmission electron microscopy. By comparing multiple membrane and (membranotropic) NP models, we found that gHNP internalization occurs by forming an intermediate state characterized by specific stabilizing interactions formed by peptide-coated nanoparticles with multicomponent model membranes. This association mechanism is mainly characterized by interactions of gH with the extracellular solvent and the polar membrane surface. At the same time, the NP core interacts with the transmembrane (cholesterol-rich) fatty phase.
Collapse
Affiliation(s)
- Sebastian Franco-Ulloa
- Molecular
Modeling and Drug Discovery Lab, Istituto
Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy
| | - Daniela Guarnieri
- Dipartimento
di Chimica e Biologia “A. Zambelli”, Università degli Studi di Salerno, Via Giovanni Paolo II, 132, Fisciano, l-84084 Salerno, Italy
| | - Laura Riccardi
- Molecular
Modeling and Drug Discovery Lab, Istituto
Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy
| | - Pier Paolo Pompa
- Nanobiointeractions
& Nanodiagnostics, Istituto Italiano
di Tecnologia, Via Morego
30, 16163 Genova, Italy
| | - Marco De Vivo
- Molecular
Modeling and Drug Discovery Lab, Istituto
Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy
| |
Collapse
|
9
|
Desikan R, Behera A, Maiti PK, Ayappa KG. Using multiscale molecular dynamics simulations to obtain insights into pore forming toxin mechanisms. Methods Enzymol 2021; 649:461-502. [PMID: 33712196 DOI: 10.1016/bs.mie.2021.01.021] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Pore forming toxins (PFTs) are virulent proteins released by several species, including many strains of bacteria, to attack and kill host cells. In this article, we focus on the utility of molecular dynamics (MD) simulations and the molecular insights gleaned from these techniques on the pore forming pathways of PFTs. In addition to all-atom simulations which are widely used, coarse-grained MARTINI models and structure-based models have also been used to study PFTs. Here, the emphasis is on methods and techniques involved while setting up, monitoring, and evaluating properties from MD simulations of PFTs in a membrane environment. We draw from several case studies to illustrate how MD simulations have provided molecular insights into protein-protein and protein-lipid interactions, lipid dynamics, conformational transitions and structures of both the oligomeric intermediates and assembled pore structures.
Collapse
Affiliation(s)
- Rajat Desikan
- Department of Chemical Engineering, Indian Institute of Science, Bengaluru, India
| | - Amit Behera
- Department of Chemical Engineering, Indian Institute of Science, Bengaluru, India
| | - Prabal K Maiti
- Centre for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bengaluru, India
| | - K Ganapathy Ayappa
- Department of Chemical Engineering, Indian Institute of Science, Bengaluru, India; Centre for Biosystems Science and Engineering, Indian Institute of Science, Bengaluru, India.
| |
Collapse
|
10
|
Souza LM, Souza FR, Reynaud F, Pimentel AS. Tuning the hydrophobicity of a coarse grained model of 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine using the experimental octanol-water partition coefficient. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2020.114132] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
|
11
|
Souza FR, Souza LMP, Pimentel AS. Recent Open Issues in Coarse Grained Force Fields. J Chem Inf Model 2020; 60:5881-5884. [PMID: 33231448 DOI: 10.1021/acs.jcim.0c01265] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
This viewpoint intends to show recent open issues of using coarse grained models in molecular dynamics simulation. It reviews the current knowledge of the comparison between experimental and simulation data of structural and physical chemical properties that depend on the hydrophilic and hydrophobic behavior of the molecule.
Collapse
Affiliation(s)
- Felipe Rodrigues Souza
- Departamento de Química, Pontifícia Universidade Católica do Rio de Janeiro, Rio de Janeiro, RJ 22453-900 Brazil
| | - Lucas Miguel Pereira Souza
- Departamento de Química, Pontifícia Universidade Católica do Rio de Janeiro, Rio de Janeiro, RJ 22453-900 Brazil
| | - Andre Silva Pimentel
- Departamento de Química, Pontifícia Universidade Católica do Rio de Janeiro, Rio de Janeiro, RJ 22453-900 Brazil
| |
Collapse
|
12
|
Siddiquee AM, Houri A, Messalea KA, Lin J, Daeneke T, Abbey B, Mechler A, Kou S. Nanoscale Probing of Cholesterol-Rich Domains in Single Bilayer Dimyristoyl-Phosphocholine Membranes Using Near-Field Spectroscopic Imaging. J Phys Chem Lett 2020; 11:9476-9484. [PMID: 33108191 DOI: 10.1021/acs.jpclett.0c02192] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Cholesterol is believed to induce the formation of membrane domains, "rafts", which are implicated in a range of natural and pathologic membrane processes. Therefore, it is important to understand the role that cholesterol plays in the formation of these structures. Here, we use label-free spectroscopic imaging to investigate cholesterol fractioning in supported bilayer membranes at nanoscale. Scattering-type scanning near-field optical microscopy (s-SNOM) was used to visualize the formation of cholesterol-induced domains in 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) membranes. Our results revealed the coexistence of phase separated domains in DMPC lipids with 10 mol % cholesterol content, whereas a mostly homogeneous bilayer was found at low (5 mol %) and high (15 mol %) cholesterol content. Near-field nano-FTIR spectroscopy was used to identify the cholesterol-rich domains based on their qualitative chemical compositions. It was determined that cholesterol binds to phosphodiester and alkyl glycerol ester moieties, likely via hydrogen bonding of the alcohol to either of the ester oxygens. The results also confirm the existence of an ideal cholesterol-lipid mixture ratio (∼15:85) with a geometrically defined packing. At lower cholesterol content there is phase separation between liquid ordered and almost neat DMPC domains. Thus, the liquid ordered phase exists at an energy minimum at a given lipid-cholesterol ratio.
Collapse
Affiliation(s)
- Arif M Siddiquee
- Department of Electronic Science, Fujian Research Center for Solid-State Lighting, Xiamen University, Xiamen 361005, China
- Department of Chemistry and Physics, La Trobe Institute for Molecular Science (LIMS), La Trobe University, Melbourne, Victoria 3086, Australia
- Australian Research Council Centre of Excellence for Advanced Molecular Imaging, Melbourne, Victoria 3086, Australia
| | - Aamd Houri
- Department of Chemistry and Physics, La Trobe Institute for Molecular Science (LIMS), La Trobe University, Melbourne, Victoria 3086, Australia
| | - Kibret A Messalea
- School of Engineering, RMIT University, Melbourne, Victoria 3001, Australia
| | - Jiao Lin
- School of Engineering, RMIT University, Melbourne, Victoria 3001, Australia
| | - Torben Daeneke
- School of Engineering, RMIT University, Melbourne, Victoria 3001, Australia
| | - Brian Abbey
- Department of Chemistry and Physics, La Trobe Institute for Molecular Science (LIMS), La Trobe University, Melbourne, Victoria 3086, Australia
- Australian Research Council Centre of Excellence for Advanced Molecular Imaging, Melbourne, Victoria 3086, Australia
| | - Adam Mechler
- Department of Chemistry and Physics, La Trobe Institute for Molecular Science (LIMS), La Trobe University, Melbourne, Victoria 3086, Australia
| | - Shanshan Kou
- Department of Chemistry and Physics, La Trobe Institute for Molecular Science (LIMS), La Trobe University, Melbourne, Victoria 3086, Australia
- Australian Research Council Centre of Excellence for Advanced Molecular Imaging, Melbourne, Victoria 3086, Australia
| |
Collapse
|
13
|
Girard M, Bereau T. Regulating Lipid Composition Rationalizes Acyl Tail Saturation Homeostasis in Ectotherms. Biophys J 2020; 119:892-899. [PMID: 32814063 DOI: 10.1016/j.bpj.2020.07.024] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2019] [Revised: 06/16/2020] [Accepted: 07/06/2020] [Indexed: 12/19/2022] Open
Abstract
Cell membranes mainly consist of lipid bilayers with an actively regulated composition. The underlying processes are still poorly understood, in particular, how the hundreds of components are controlled. Cholesterol has been found to correlate with phospholipid saturation for reasons that remain unclear. To better understand the link between cell membrane regulation and chemical composition, we establish a computational framework based on chemical reaction networks, resulting in multiple semigrand canonical ensembles. By running computer simulations, we show that regulating the chemical potential of lipid species is sufficient to reproduce the experimentally observed increase in acyl tail saturation with added cholesterol. Our model proposes a different picture of lipid regulation in which components can be regulated passively instead of actively. In this picture, phospholipid acyl tail composition naturally adapts to added molecules such as cholesterol or proteins. A comparison between regulated membranes with commonly studied ternary model membranes shows stark differences: for instance, correlation lengths and viscosities observed are independent of lipid chemical affinity.
Collapse
Affiliation(s)
- Martin Girard
- Max Planck Institute for Polymer Research, Mainz, Germany.
| | - Tristan Bereau
- Max Planck Institute for Polymer Research, Mainz, Germany; Van 't Hoff Institute for Molecular Sciences and Informatics Institute, University of Amsterdam, Amsterdam, the Netherlands
| |
Collapse
|
14
|
Mostofian B, Johnson QR, Smith JC, Cheng X. Carotenoids promote lateral packing and condensation of lipid membranes. Phys Chem Chem Phys 2020; 22:12281-12293. [PMID: 32432296 DOI: 10.1039/d0cp01031f] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Carotenoids are pigment molecules that protect biomembranes against degradation and may be involved in the formation of functional bacterial membrane microdomains. Little is known on whether different types of carotenoids have different effects on the membrane or if there is any concentration dependence of these effects. In this work, we present results from molecular dynamics simulations of phospholipid bilayers containing different amounts of either β-carotene or zeaxanthin. Both β-carotene and zeaxanthin show the ability to laterally condense the membrane lipids and reduce their inter-leaflet interactions. With increasing concentrations, both carotenoids increase the bilayer thickness and rigidity. The results reveal that carotenoids have similar effects to cholesterol on regulating the behavior of fluid-phase membranes, suggesting that they could function as sterol substitutes and confirming their potential role in the formation of functional membrane domains.
Collapse
Affiliation(s)
- Barmak Mostofian
- Center for Molecular Biophysics, Oak Ridge National Lab, Oak Ridge, TN 37830, USA.
| | | | | | | |
Collapse
|
15
|
van Hilten N, Stroh KS, Risselada HJ. Membrane Thinning Induces Sorting of Lipids and the Amphipathic Lipid Packing Sensor (ALPS) Protein Motif. Front Physiol 2020; 11:250. [PMID: 32372966 PMCID: PMC7177014 DOI: 10.3389/fphys.2020.00250] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2020] [Accepted: 03/05/2020] [Indexed: 12/16/2022] Open
Abstract
Heterogeneities (e.g., membrane proteins and lipid domains) and deformations (e.g., highly curved membrane regions) in biological lipid membranes cause lipid packing defects that may trigger functional sorting of lipids and membrane-associated proteins. To study these phenomena in a controlled and efficient way within molecular simulations, we developed an external field protocol that artificially enhances packing defects in lipid membranes by enforcing local thinning of a flat membrane region. For varying lipid compositions, we observed strong thinning-induced depletion or enrichment, depending on the lipid's intrinsic shape and its effect on a membrane's elastic modulus. In particular, polyunsaturated and lysolipids are strongly attracted to regions high in packing defects, whereas phosphatidylethanolamine (PE) lipids and cholesterol are strongly repelled from it. Our results indicate that externally imposed changes in membrane thickness, area, and curvature are underpinned by shared membrane elastic principles. The observed sorting toward the thinner membrane region is in line with the sorting expected for a positively curved membrane region. Furthermore, we have demonstrated that the amphipathic lipid packing sensor (ALPS) protein motif, a known curvature and packing defect sensor, is effectively attracted to thinner membrane regions. By extracting the force that drives amphipathic molecules toward the thinner region, our thinning protocol can directly quantify and score the lipid packing sensing of different amphipathic molecules. In this way, our protocol paves the way toward high-throughput exploration of potential defect- and curvature-sensing motifs, making it a valuable addition to the molecular simulation toolbox.
Collapse
Affiliation(s)
- Niek van Hilten
- Leiden Institute of Chemistry, Leiden University, Leiden, Netherlands
| | - Kai Steffen Stroh
- Institute for Theoretical Physics, Georg August University Göttingen, Göttingen, Germany
| | - Herre Jelger Risselada
- Leiden Institute of Chemistry, Leiden University, Leiden, Netherlands.,Institute for Theoretical Physics, Georg August University Göttingen, Göttingen, Germany
| |
Collapse
|
16
|
Fornasier F, Souza LMP, Souza FR, Reynaud F, Pimentel AS. Lipophilicity of Coarse-Grained Cholesterol Models. J Chem Inf Model 2020; 60:569-577. [DOI: 10.1021/acs.jcim.9b00830] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Franccesca Fornasier
- Departamento de Quı́mica, Pontifı́cia Universidade Católica do Rio de Janeiro Rio de Janeiro, RJ 22453-900, Brazil
| | - Lucas M. P. Souza
- Departamento de Quı́mica, Pontifı́cia Universidade Católica do Rio de Janeiro Rio de Janeiro, RJ 22453-900, Brazil
| | - Felipe R. Souza
- Departamento de Quı́mica, Pontifı́cia Universidade Católica do Rio de Janeiro Rio de Janeiro, RJ 22453-900, Brazil
| | - Franceline Reynaud
- Institut Galien Paris-Sud, CNRS 8612, Université Paris-Saclay, 92290 Châtenay, Malabry, France
- Faculdade de Farmácia, Universidade Federal do Rio de Janeiro, Rio de Janeiro, RJ 21941-902, Brazil
| | - Andre S. Pimentel
- Departamento de Quı́mica, Pontifı́cia Universidade Católica do Rio de Janeiro Rio de Janeiro, RJ 22453-900, Brazil
| |
Collapse
|
17
|
Shen H, Wu Z, Zhao K, Yang H, Deng M, Wen S. Effect of Cholesterol and 6-Ketocholestanol on Membrane Dipole Potential and Sterol Flip-Flop Motion in Bilayer Membranes. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2019; 35:11232-11241. [PMID: 31373497 DOI: 10.1021/acs.langmuir.9b01802] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
A variety of experimental and theoretical approaches have been employed to investigate the sterol flip-flop motion in lipid bilayer membranes. However, the sterol effect on the dipole potential of lipid bilayer membranes is less well studied and the influence of dipole potential on sterol flip-flop motion in lipid bilayer membranes is less well understood. In our previous works, we have demonstrated the performance of our coarse-grained (CG) model in the computation of the dipole potential. In this work, five 30 μs CG simulations of dimyristoylphosphatidylcholine (DMPC) bilayers were carried out at different sterol concentrations (in a range from 10 to 50% mole fraction). Then, a comparison was made between the effects of cholesterol (CHOL) and 6-ketocholestanol (6-KC) on the dipole potential of DMPC lipid bilayers as well as the sterol flip-flop motion. Our CG simulations show that the membrane dipole potential is impacted more significantly by 6-KC than by CHOL. This finding is consistent with recent experimental studies. Meanwhile, our work suggests that the sterol-sterol interactions (in particular, electrostatic interactions) should be critical to the formation of sterol-sterol clusters, which would hinder the sterol flip-flop motion inside the lipid bilayers. This is in support of the recent experimental study on the sterol transportation in lipid bilayer membranes.
Collapse
Affiliation(s)
- Hujun Shen
- Guizhou Provincial Key Laboratory of Computational Nano-Material Science, Guizhou Synergetic Innovation Center of Scientific Big Data for Advanced Manufacturing Technology , Guizhou Education University , No. 115, Gaoxin Road , Guiyang , Guizhou 550018 , P. R. China
- School of Information , Guizhou University of Finance and Economics , University City of Huaxi District, Guiyang , Guizhou 550025 , P. R. China
| | - Zhenhua Wu
- School of Information , Guizhou University of Finance and Economics , University City of Huaxi District, Guiyang , Guizhou 550025 , P. R. China
| | - Kun Zhao
- School of Information , Guizhou University of Finance and Economics , University City of Huaxi District, Guiyang , Guizhou 550025 , P. R. China
| | - Hengxiu Yang
- Guizhou Provincial Key Laboratory of Computational Nano-Material Science, Guizhou Synergetic Innovation Center of Scientific Big Data for Advanced Manufacturing Technology , Guizhou Education University , No. 115, Gaoxin Road , Guiyang , Guizhou 550018 , P. R. China
| | - Mingsen Deng
- Guizhou Provincial Key Laboratory of Computational Nano-Material Science, Guizhou Synergetic Innovation Center of Scientific Big Data for Advanced Manufacturing Technology , Guizhou Education University , No. 115, Gaoxin Road , Guiyang , Guizhou 550018 , P. R. China
- School of Information , Guizhou University of Finance and Economics , University City of Huaxi District, Guiyang , Guizhou 550025 , P. R. China
| | - Shuiguo Wen
- Guizhou Provincial Key Laboratory of Computational Nano-Material Science, Guizhou Synergetic Innovation Center of Scientific Big Data for Advanced Manufacturing Technology , Guizhou Education University , No. 115, Gaoxin Road , Guiyang , Guizhou 550018 , P. R. China
| |
Collapse
|
18
|
Nalakarn P, Boonnoy P, Nisoh N, Karttunen M, Wong-Ekkabut J. Dependence of fullerene aggregation on lipid saturation due to a balance between entropy and enthalpy. Sci Rep 2019; 9:1037. [PMID: 30705323 PMCID: PMC6355782 DOI: 10.1038/s41598-018-37659-4] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2018] [Accepted: 12/11/2018] [Indexed: 11/23/2022] Open
Abstract
It is well-known that fullerenes aggregate inside lipid membranes and that increasing the concentration may lead to (lethal) membrane rupture. It is not known, however, how aggregation and rupture depend on the lipid type, what physical mechanisms control this behavior and what experimental signatures detect such changes in membranes. In this paper, we attempt to answer these questions with molecular simulations, and we show that aggregation and membrane damage depend critically on the degree of saturation of the lipid acyl chains: unsaturated bonds, or "kinks", impose a subtle but crucial compartmentalization of the bilayer into core and surface regions leading to three distinct fullerene density maxima. In contrast, when the membrane has only fully saturated lipids, fullerenes prefer to be located close to the surface under the head groups until the concentration becomes too large and the fullerenes begin clustering. No clustering is observed in membranes with unsaturated lipids. The presence of "kinks" reverses the free energy balance; although the overall free energy profiles are similar, entropy is the dominant component in unsaturated bilayers whereas enthalpy controls the fully saturated ones. Fully saturated systems show two unique signatures: 1) membrane thickness behaves non-monotonously while the area per lipid increases monotonously. We propose this as a potential reason for the observations of low fullerene concentrations being effective against bacteria. 2) The fullerene-fullerene radial distribution function (RDF) shows splitting of the second peak indicating the emergence short-range order and the importance of the second-nearest neighbor interactions. Similar second peak splitting has been reported in metal glasses.
Collapse
Affiliation(s)
- Pornkamon Nalakarn
- Department of Physics, Faculty of Science, Kasetsart University, Bangkok, 10900, Thailand
- Computational Biomodelling Laboratory for Agricultural Science and Technology (CBLAST), Faculty of Science, Kasetsart University, Bangkok, 10900, Thailand
- Thailand Center of Excellence in Physics (ThEP Center), Commission on Higher Education, Bangkok, 10400, Thailand
| | - Phansiri Boonnoy
- Department of Physics, Faculty of Science, Kasetsart University, Bangkok, 10900, Thailand
- Computational Biomodelling Laboratory for Agricultural Science and Technology (CBLAST), Faculty of Science, Kasetsart University, Bangkok, 10900, Thailand
| | - Nililla Nisoh
- Department of Physics, Faculty of Science, Kasetsart University, Bangkok, 10900, Thailand
- Computational Biomodelling Laboratory for Agricultural Science and Technology (CBLAST), Faculty of Science, Kasetsart University, Bangkok, 10900, Thailand
- Thailand Center of Excellence in Physics (ThEP Center), Commission on Higher Education, Bangkok, 10400, Thailand
| | - Mikko Karttunen
- Department of Chemistry and Department of Applied Mathematics, Western University, 1151 Richmond Street, London, Ontario, N6A 5B7, Canada.
| | - Jirasak Wong-Ekkabut
- Department of Physics, Faculty of Science, Kasetsart University, Bangkok, 10900, Thailand.
- Computational Biomodelling Laboratory for Agricultural Science and Technology (CBLAST), Faculty of Science, Kasetsart University, Bangkok, 10900, Thailand.
- Thailand Center of Excellence in Physics (ThEP Center), Commission on Higher Education, Bangkok, 10400, Thailand.
- Specialized Center of Rubber and Polymer Materials for Agriculture and Industry (RPM), Faculty of Science, Kasetsart University, Bangkok, 10900, Thailand.
| |
Collapse
|
19
|
Parameterization of a coarse-grained model of cholesterol with point-dipole electrostatics. J Comput Aided Mol Des 2018; 32:1259-1271. [DOI: 10.1007/s10822-018-0164-4] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2018] [Accepted: 09/12/2018] [Indexed: 12/14/2022]
|
20
|
Chen L, Li X, Zhang Y, Chen T, Xiao S, Liang H. Morphological and mechanical determinants of cellular uptake of deformable nanoparticles. NANOSCALE 2018; 10:11969-11979. [PMID: 29904774 DOI: 10.1039/c8nr01521j] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
Understanding the interactions of nanoparticles (NPs) with cell membranes and regulating their cellular uptake processes are of fundamental importance to the design of drug delivery systems with minimum toxicity, high efficiency and long circulation time. Employing the procedure of coarse-graining, we built an elastically deformable NP model with tunable morphological and mechanical properties. We found that the cellular uptake of deformable NPs depends on their shape: an increase in the particle elasticity significantly slows the uptake rate of spherical NPs, slightly retards that of prolate NPs, and promotes the uptake of oblate NPs. The intrinsic mechanisms have been carefully investigated through analysis of the endocytic mechanisms and free energy calculations. These findings provide unique insights into how deformable NPs penetrate across cell membranes and offer novel possibilities for designing effective NP-based carriers for drug delivery.
Collapse
Affiliation(s)
- Liping Chen
- CAS Key Laboratory of Soft Matter Chemistry, iChEM (Collaborative Innovation Center of Chemistry for Energy Materials), Department of Polymer Science and Engineering, University of Science and Technology of China, Hefei, Anhui 230026, P. R. China.
| | | | | | | | | | | |
Collapse
|
21
|
Kindlein M, Elts E, Briesen H. Phospholipids in chocolate: Structural insights and mechanistic explanations of rheological behavior by coarse-grained molecular dynamics simulations. J FOOD ENG 2018. [DOI: 10.1016/j.jfoodeng.2018.02.014] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
22
|
Shen H, Deng M, Wu Z, Zhang J, Zhang Y, Gao C, Cen C. Effect of Cholesterol on Membrane Dipole Potential: Atomistic and Coarse-Grained Molecular Dynamics Simulations. J Chem Theory Comput 2018; 14:3780-3795. [DOI: 10.1021/acs.jctc.8b00092] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Hujun Shen
- Guizhou Provincial Key Laboratory of Computational Nano-Material Science, Guizhou Synergetic Innovation Center of Scientific Big Data for Advanced Manufacturing Technology, Guizhou Education University, No. 115, Gaoxin Road, Guiyang, Guizhou 550018, P. R. China
- Guizhou University of Finance and Economics, School of Information, University City of Huaxi District, Guiyang, Guizhou 550025, P. R. China
| | - Mingsen Deng
- Guizhou Provincial Key Laboratory of Computational Nano-Material Science, Guizhou Synergetic Innovation Center of Scientific Big Data for Advanced Manufacturing Technology, Guizhou Education University, No. 115, Gaoxin Road, Guiyang, Guizhou 550018, P. R. China
- Guizhou University of Finance and Economics, School of Information, University City of Huaxi District, Guiyang, Guizhou 550025, P. R. China
| | - Zhenhua Wu
- Guizhou University of Finance and Economics, School of Information, University City of Huaxi District, Guiyang, Guizhou 550025, P. R. China
| | - Jihua Zhang
- Guizhou Provincial Key Laboratory of Computational Nano-Material Science, Guizhou Synergetic Innovation Center of Scientific Big Data for Advanced Manufacturing Technology, Guizhou Education University, No. 115, Gaoxin Road, Guiyang, Guizhou 550018, P. R. China
| | - Yachao Zhang
- Guizhou Provincial Key Laboratory of Computational Nano-Material Science, Guizhou Synergetic Innovation Center of Scientific Big Data for Advanced Manufacturing Technology, Guizhou Education University, No. 115, Gaoxin Road, Guiyang, Guizhou 550018, P. R. China
| | - Chengui Gao
- Guizhou Provincial Key Laboratory of Computational Nano-Material Science, Guizhou Synergetic Innovation Center of Scientific Big Data for Advanced Manufacturing Technology, Guizhou Education University, No. 115, Gaoxin Road, Guiyang, Guizhou 550018, P. R. China
| | - Cao Cen
- Guizhou Provincial Key Laboratory of Computational Nano-Material Science, Guizhou Synergetic Innovation Center of Scientific Big Data for Advanced Manufacturing Technology, Guizhou Education University, No. 115, Gaoxin Road, Guiyang, Guizhou 550018, P. R. China
| |
Collapse
|
23
|
Jaschonek S, Cascella M, Gauss J, Diezemann G, Milano G. Intramolecular structural parameters are key modulators of the gel-liquid transition in coarse grained simulations of DPPC and DOPC lipid bilayers. Biochem Biophys Res Commun 2018; 498:327-333. [PMID: 29101041 DOI: 10.1016/j.bbrc.2017.10.132] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2017] [Revised: 09/29/2017] [Accepted: 10/23/2017] [Indexed: 12/11/2022]
Abstract
The capability of coarse-grained models based on the MARTINI mapping to reproduce the gel-liquid phase transition in saturated and unsaturated model lipids was investigated. We found that the model is able to reproduce a lower critical temperature for 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) with respect to 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC). Nonetheless, the appearance of a gel phase for DOPC is strictly dependent on the intramolecular parameters chosen to model its molecular structure. In particular, we show that the bending angle at the coarse-grained bead corresponding to the unsaturated carbon-carbon bond acts as an order parameter determining the temperature of the phase transition. Structural analysis of the molecular dynamics simulations runs evidences that in the gel phase, the packing of the lipophilic tails of DOPC assume a different conformation than in the liquid phase. In the latter phase, the DOPC geometry resembles that of the relaxed free molecule. DPPC:DOPC mixtures show a single phase transition temperature, indicating that the observation of a phase separation between the two lipids requires the simulation of systems with sizes much larger than the ones used here.
Collapse
Affiliation(s)
- Stefan Jaschonek
- Institut für Physikalische Chemie, Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany.
| | - Michele Cascella
- Department of Chemistry and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, Postboks 1033 Blindern, N-0315 Oslo, Norway.
| | - Jürgen Gauss
- Institut für Physikalische Chemie, Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany.
| | - Gregor Diezemann
- Institut für Physikalische Chemie, Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany.
| | - Giuseppe Milano
- Dipartimento di Chimica e Biologia, Università di Salerno, Via Giovanni Paolo II 132, I-84084 Fisciano, Italy.
| |
Collapse
|
24
|
Li B, Chen Q, Huang S, Liu H. Developing structure and thermodynamic properties-consistent coarse-grained model for random copolymer systems. POLYMER 2017. [DOI: 10.1016/j.polymer.2017.07.016] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
|
25
|
Fosso-Tande J, Black C, G. Aller S, Lu L, D. Hills Jr R. Simulation of lipid-protein interactions with the CgProt force field. AIMS MOLECULAR SCIENCE 2017. [DOI: 10.3934/molsci.2017.3.352] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
|
26
|
Melo MN, Ingólfsson HI, Marrink SJ. Parameters for Martini sterols and hopanoids based on a virtual-site description. J Chem Phys 2016; 143:243152. [PMID: 26723637 DOI: 10.1063/1.4937783] [Citation(s) in RCA: 104] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Abstract
Sterols play an essential role in modulating bilayer structure and dynamics. Coarse-grained molecular dynamics parameters for cholesterol and related molecules are available for the Martini force field and have been successfully used in multiple lipid bilayer studies. In this work, we focus on the use of virtual sites as a means of increasing the stability of cholesterol and cholesterol-like structures. We improve and extend the Martini parameterization of sterols in four different ways: 1-the cholesterol parameters were adapted to make use of virtual interaction sites, which markedly improves numerical stability; 2-cholesterol parameters were also modified to address reported shortcomings in reproducing correct lipid phase behavior in mixed membranes; 3-parameters for ergosterol were created and adapted from cholesterols; and 4-parameters for the hopanoid class of bacterial polycyclic molecules were created, namely, for hopane, diploptene, bacteriohopanetetrol, and for their polycyclic base structure.
Collapse
Affiliation(s)
- M N Melo
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - H I Ingólfsson
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - S J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| |
Collapse
|
27
|
Shi YZ, Jin L, Wang FH, Zhu XL, Tan ZJ. Predicting 3D Structure, Flexibility, and Stability of RNA Hairpins in Monovalent and Divalent Ion Solutions. Biophys J 2016; 109:2654-2665. [PMID: 26682822 DOI: 10.1016/j.bpj.2015.11.006] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2015] [Revised: 10/09/2015] [Accepted: 11/06/2015] [Indexed: 10/24/2022] Open
Abstract
A full understanding of RNA-mediated biology would require the knowledge of three-dimensional (3D) structures, structural flexibility, and stability of RNAs. To predict RNA 3D structures and stability, we have previously proposed a three-bead coarse-grained predictive model with implicit salt/solvent potentials. In this study, we further develop the model by improving the implicit-salt electrostatic potential and including a sequence-dependent coaxial stacking potential to enable the model to simulate RNA 3D structure folding in divalent/monovalent ion solutions. The model presented here can predict 3D structures of RNA hairpins with bulges/internal loops (<77 nucleotides) from their sequences at the corresponding experimental ion conditions with an overall improved accuracy compared to the experimental data; the model also makes reliable predictions for the flexibility of RNA hairpins with bulge loops of different lengths at several divalent/monovalent ion conditions. In addition, the model successfully predicts the stability of RNA hairpins with various loops/stems in divalent/monovalent ion solutions.
Collapse
Affiliation(s)
- Ya-Zhou Shi
- Department of Physics and Key Laboratory of Artificial Micro- and Nano-structures of the Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China
| | - Lei Jin
- Department of Physics and Key Laboratory of Artificial Micro- and Nano-structures of the Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China
| | - Feng-Hua Wang
- Engineering Training Center, Jianghan University, Wuhan, China
| | - Xiao-Long Zhu
- Department of Physics, School of Physics and Information Engineering, Jianghan University, Wuhan, China
| | - Zhi-Jie Tan
- Department of Physics and Key Laboratory of Artificial Micro- and Nano-structures of the Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China.
| |
Collapse
|
28
|
Shimizu H, Matsunaga S, Yamada T, Kobayashi T, Kawai M. Formation of Ordered Phospholipid Monolayer on a Hydrophilically Modified Au(111) Substrate. ACS NANO 2016; 10:7811-7820. [PMID: 27494363 DOI: 10.1021/acsnano.6b03421] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
The molecular arrangement of phospholipid molecules was investigated on a hydrophilically modified gold surface within an aqueous solution by scanning tunneling microscopy. By suspending phospholipid (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine, POPC) nanoparticles in the aqueous electrolyte surrounding a hydrophilically modified gold (111) substrate with 3-mercaptopropionic acid (SH-C2H4-COOH, 3-MPA), well-ordered adlattices of POPC were observed. Traces of particle fusion were visualized before formation of the adlattice. Addition of cholesterol to the suspension seems to facilitate accommodation of POPC on this surface. The observed unit cells of POPC adlattices had dimensions of 0.5 nm × 1.9-2.5 nm. By high-resolution imaging, each unit cell was discerned to be occupied by one upright POPC molecule. The POPC + cholesterol suspension also leads to formation of a flat integrated POPC layer, which may be a lipid bilayer covering the surface.
Collapse
Affiliation(s)
- Hiroaki Shimizu
- RIKEN , 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
- Department of Advanced Materials Science, The University of Tokyo , 5-1-5 Kashiwanoha, Kashiwa, Chiba 277-8561, Japan
| | - Soichiro Matsunaga
- RIKEN , 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
- Department of Advanced Materials Science, The University of Tokyo , 5-1-5 Kashiwanoha, Kashiwa, Chiba 277-8561, Japan
| | - Taro Yamada
- RIKEN , 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | | | - Maki Kawai
- RIKEN , 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
- Department of Advanced Materials Science, The University of Tokyo , 5-1-5 Kashiwanoha, Kashiwa, Chiba 277-8561, Japan
| |
Collapse
|
29
|
Wang Y, Gkeka P, Fuchs JE, Liedl KR, Cournia Z. DPPC-cholesterol phase diagram using coarse-grained Molecular Dynamics simulations. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2016; 1858:2846-2857. [PMID: 27526680 DOI: 10.1016/j.bbamem.2016.08.005] [Citation(s) in RCA: 77] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/17/2015] [Revised: 08/06/2016] [Accepted: 08/11/2016] [Indexed: 12/30/2022]
Abstract
Cholesterol-phospholipid bilayers continue to be the current state of the art in membrane models and serve as representative systems for studying the effect of cholesterol on the cell membrane. As the mixing of different lipid species requires long spatio-temporal scales, coarse-grained models have gained increasing popularity in modeling such membrane systems. In this paper, a systematic study of the MARTINI coarse-grained model for the DPPC-cholesterol binary system has been performed. We construct the phase diagram of DPPC lipid bilayers in the presence of different cholesterol concentrations and at different temperatures using coarse-grained Molecular Dynamics (MD) simulations with the MARTINI force field. The phase diagram based on the condensation effect is directly comparable to available experimental data and demonstrates qualitative agreement over all cholesterol concentrations. Self-assembled bilayers quantitatively reproduce experimental observables, such as lateral diffusion of lipids, electron density, area per lipid and lipid order parameters. The phase diagram of the DPPC-cholesterol binary system also reveals the profound effect of cholesterol on the physical properties of phospholipid bilayers such lipid order, diffusion, and fluidity. Cholesterol induces the liquid-ordered phase, which increases the fluidity of the phospholipid hydrocarbon chains above the gel to liquid-crystalline phase transition temperature and decreases it below the phase transition. The present study suggests that the MARTINI force field can be successfully used to obtain molecular level insights into cholesterol-DPPC model membranes.
Collapse
Affiliation(s)
- Yin Wang
- Department of Theoretical Chemistry and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| | - Paraskevi Gkeka
- Biomedical Research Foundation, Academy of Athens, 4 Soranou Ephessiou, 11527, Athens, Greece
| | - Julian E Fuchs
- Department of Theoretical Chemistry and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| | - Klaus R Liedl
- Department of Theoretical Chemistry and Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria.
| | - Zoe Cournia
- Biomedical Research Foundation, Academy of Athens, 4 Soranou Ephessiou, 11527, Athens, Greece.
| |
Collapse
|
30
|
Hills RD, McGlinchey N. Model parameters for simulation of physiological lipids. J Comput Chem 2016; 37:1112-8. [PMID: 26864972 PMCID: PMC5067697 DOI: 10.1002/jcc.24324] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2015] [Revised: 12/19/2015] [Accepted: 01/17/2016] [Indexed: 12/16/2022]
Abstract
Coarse grain simulation of proteins in their physiological membrane environment can offer insight across timescales, but requires a comprehensive force field. Parameters are explored for multicomponent bilayers composed of unsaturated lipids DOPC and DOPE, mixed-chain saturation POPC and POPE, and anionic lipids found in bacteria: POPG and cardiolipin. A nonbond representation obtained from multiscale force matching is adapted for these lipids and combined with an improved bonding description of cholesterol. Equilibrating the area per lipid yields robust bilayer simulations and properties for common lipid mixtures with the exception of pure DOPE, which has a known tendency to form nonlamellar phase. The models maintain consistency with an existing lipid-protein interaction model, making the force field of general utility for studying membrane proteins in physiologically representative bilayers.
Collapse
Affiliation(s)
- Ronald D Hills
- Department of Pharmaceutical Sciences, College of Pharmacy, University of New England, 716 Stevens Ave, Portland, Maine, 04103
| | - Nicholas McGlinchey
- Department of Pharmaceutical Sciences, College of Pharmacy, University of New England, 716 Stevens Ave, Portland, Maine, 04103
| |
Collapse
|
31
|
Coarse-grained bond and angle distributions from atomistic simulations: On the systematic parameterisation of lipid models. J Mol Graph Model 2016; 63:57-64. [DOI: 10.1016/j.jmgm.2015.11.009] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2015] [Revised: 11/06/2015] [Accepted: 11/09/2015] [Indexed: 11/23/2022]
|
32
|
Garg S, Castro-Roman F, Porcar L, Butler P, Bautista PJ, Krzyzanowski N, Perez-Salas U. Cholesterol solubility limit in lipid membranes probed by small angle neutron scattering and MD simulations. SOFT MATTER 2014; 10:9313-9317. [PMID: 25338228 DOI: 10.1039/c4sm01219d] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
The solubility limits of cholesterol in small unilamellar vesicles made of POPS and POPC were probed using Small Angle Neutron Scattering (SANS) and coarse grained (CG) molecular dynamics (MD) simulations. SANS, being non-invasive, allowed the direct and quantitative measurement of cholesterol in intact vesicles. Our experimental measurements reveal a 61% mole fraction solubility limit of cholesterol in POPC, consistent with previous studies. However, in POPS the solubility limit of cholesterol is found to be 73% mole fraction. Previous work reports solubility limits of cholesterol in POPS varying significantly, ranging from 36% up to 66%. The CG MD simulations are in remarkable quantitative agreement with our experimental results showing similar solubility limits. Further, neither experiments nor simulations show evidence of stable nanodomains of cholesterol in POPS membranes as suggested in some previous reports.
Collapse
Affiliation(s)
- Sumit Garg
- Materials Science Division, Argonne National Laboratory, Lemont, Illinois, USA.
| | | | | | | | | | | | | |
Collapse
|
33
|
Róg T, Vattulainen I. Cholesterol, sphingolipids, and glycolipids: what do we know about their role in raft-like membranes? Chem Phys Lipids 2014; 184:82-104. [PMID: 25444976 DOI: 10.1016/j.chemphyslip.2014.10.004] [Citation(s) in RCA: 140] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2014] [Revised: 10/24/2014] [Accepted: 10/25/2014] [Indexed: 12/14/2022]
Abstract
Lipids rafts are considered to be functional nanoscale membrane domains enriched in cholesterol and sphingolipids, characteristic in particular of the external leaflet of cell membranes. Lipids, together with membrane-associated proteins, are therefore considered to form nanoscale units with potential specific functions. Although the understanding of the structure of rafts in living cells is quite limited, the possible functions of rafts are widely discussed in the literature, highlighting their importance in cellular functions. In this review, we discuss the understanding of rafts that has emerged based on recent atomistic and coarse-grained molecular dynamics simulation studies on the key lipid raft components, which include cholesterol, sphingolipids, glycolipids, and the proteins interacting with these classes of lipids. The simulation results are compared to experiments when possible.
Collapse
Affiliation(s)
- Tomasz Róg
- Department of Physics, Tampere University of Technology, Tampere, Finland
| | - Ilpo Vattulainen
- Department of Physics, Tampere University of Technology, Tampere, Finland; MEMPHYS-Center for Biomembrane Physics, University of Southern Denmark, Odense, Denmark.
| |
Collapse
|
34
|
Koldsø H, Shorthouse D, Hélie J, Sansom MSP. Lipid clustering correlates with membrane curvature as revealed by molecular simulations of complex lipid bilayers. PLoS Comput Biol 2014; 10:e1003911. [PMID: 25340788 PMCID: PMC4207469 DOI: 10.1371/journal.pcbi.1003911] [Citation(s) in RCA: 170] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2014] [Accepted: 09/16/2014] [Indexed: 12/21/2022] Open
Abstract
Cell membranes are complex multicomponent systems, which are highly heterogeneous in the lipid distribution and composition. To date, most molecular simulations have focussed on relatively simple lipid compositions, helping to inform our understanding of in vitro experimental studies. Here we describe on simulations of complex asymmetric plasma membrane model, which contains seven different lipids species including the glycolipid GM3 in the outer leaflet and the anionic lipid, phosphatidylinositol 4,5-bisphophate (PIP2), in the inner leaflet. Plasma membrane models consisting of 1500 lipids and resembling the in vivo composition were constructed and simulations were run for 5 µs. In these simulations the most striking feature was the formation of nano-clusters of GM3 within the outer leaflet. In simulations of protein interactions within a plasma membrane model, GM3, PIP2, and cholesterol all formed favorable interactions with the model α-helical protein. A larger scale simulation of a model plasma membrane containing 6000 lipid molecules revealed correlations between curvature of the bilayer surface and clustering of lipid molecules. In particular, the concave (when viewed from the extracellular side) regions of the bilayer surface were locally enriched in GM3. In summary, these simulations explore the nanoscale dynamics of model bilayers which mimic the in vivo lipid composition of mammalian plasma membranes, revealing emergent nanoscale membrane organization which may be coupled both to fluctuations in local membrane geometry and to interactions with proteins. Cell membranes play important roles in vivo both in shielding the cell interior from the surrounding environment and in cell function through lipid components of the membrane having roles in controlling protein function, cell signaling etc. We employ molecular dynamics simulations to explore the behavior of biologically realistic membrane models. Our simulations reveal nano-domain clustering of the glycolipid GM3 and to a lesser extent of the anionic lipid phosphatidylinositol 4,5-bisphophate (PIP2). When including transmembrane proteins we are able to observe preferential interactions of known regulatory lipids (e.g. GM3, PIP2 and cholesterol) with the proteins. Membrane curvature is shown to be coupled to the local lipid composition, suggestive of a link between lipid nano-domains and membrane geometry.
Collapse
Affiliation(s)
- Heidi Koldsø
- Department of Biochemistry, University of Oxford, Oxford, United Kingdom
| | - David Shorthouse
- Department of Biochemistry, University of Oxford, Oxford, United Kingdom
| | - Jean Hélie
- Department of Biochemistry, University of Oxford, Oxford, United Kingdom
| | - Mark S. P. Sansom
- Department of Biochemistry, University of Oxford, Oxford, United Kingdom
- * E-mail:
| |
Collapse
|