1
|
Mayorga LS, Masone D. The Secret Ballet Inside Multivesicular Bodies. ACS NANO 2024; 18:15651-15660. [PMID: 38830824 DOI: 10.1021/acsnano.4c01590] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2024]
Abstract
Lipid bilayers possess the capacity for self-assembly due to the amphipathic nature of lipid molecules, which have both hydrophobic and hydrophilic regions. When confined, lipid bilayers exhibit astonishing versatility in their forms, adopting diverse shapes that are challenging to observe through experimental means. Exploiting this adaptability, lipid structures motivate the development of bio-inspired mechanomaterials and integrated nanobio-interfaces that could seamlessly merge with biological entities, ultimately bridging the gap between synthetic and biological systems. In this work, we demonstrate how, in numerical simulations of multivesicular bodies, a fascinating evolution unfolds from an initial semblance of order toward states of higher entropy over time. We observe dynamic rearrangements in confined vesicles that reveal unexpected limit shapes of distinct geometric patterns. We identify five structures as the basic building blocks that systematically repeat under various conditions of size and composition. Moreover, we observe more complex and less frequent shapes that emerge in confined spaces. Our results provide insights into the dynamics of multivesicular systems, offering a richer understanding of how confined lipid bodies spontaneously self-organize.
Collapse
Affiliation(s)
- Luis S Mayorga
- Instituto de Histología y Embriología de Mendoza (IHEM), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Universidad Nacional de Cuyo (UNCuyo), 5500 Mendoza, Argentina
- Facultad de Ciencias Exactas y Naturales, Universidad Nacional de Cuyo (UNCuyo), 5500, Mendoza, Argentina
| | - Diego Masone
- Instituto de Histología y Embriología de Mendoza (IHEM), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Universidad Nacional de Cuyo (UNCuyo), 5500 Mendoza, Argentina
- Facultad de Ingeniería, Universidad Nacional de Cuyo (UNCuyo), 5500 Mendoza, Argentina
| |
Collapse
|
2
|
van der Pol RI, Brinkmann BW, Sevink GJA. Analyzing Lipid Membrane Defects via a Coarse-Grained to Triangulated Surface Map: The Role of Lipid Order and Local Curvature in Molecular Binding. J Chem Theory Comput 2024; 20:2888-2900. [PMID: 38537131 PMCID: PMC11008102 DOI: 10.1021/acs.jctc.4c00082] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2024] [Revised: 03/14/2024] [Accepted: 03/14/2024] [Indexed: 04/10/2024]
Abstract
Lipid packing defects are known to serve as quantitative indicators for protein binding to lipid membranes. This paper presents a protocol for mapping molecular lipid detail onto a triangulated continuum leaflet representation. Besides establishing the desired forward counterpart to the existing inverse TS2CG map, this coarse-grained to triangulated surface (CG2TS) map enables straightforward extraction of the defect characteristics for any membrane geometry found in nature. We have applied our protocol to investigate the role of local curvature and varying lipid packing on the defect constant π. We find that the defect size is greatly influenced by both factors, arguing strongly against the usual assignment of a single defect constant in the case of more realistic membrane conditions. An important discovery is that lipids in the gel phase produce larger defects, or a higher π, in domains of high (local) curvature than the same lipid in a liquid phase of any curvature. This finding suggests that membranes featuring very ordered lipid packing can bind proteins via large defects in curved regions. Finally, we propose a route for estimating defect constants directly from the standard membrane properties. Identifying the precise role of composition, lipid (tail) order, and (local) curvature in defects for the irregular lipid structures that are (temporally) present in many biological processes is instrumental for obtaining fundamental insight as well as for a rational design of membrane binding targets.
Collapse
Affiliation(s)
- Rianne
W. I. van der Pol
- Leiden
Institute of Chemistry, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Bregje W. Brinkmann
- Institute
of Environmental Sciences, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - G. J. Agur Sevink
- Leiden
Institute of Chemistry, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
3
|
Fábián B, Javanainen M. Diffusion Analyses along Mean and Gaussian-Curved Membranes with CurD. J Phys Chem Lett 2024; 15:3214-3220. [PMID: 38483514 DOI: 10.1021/acs.jpclett.4c00338] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/22/2024]
Abstract
Curved cellular membranes are both abundant and functionally relevant. While novel tomography approaches reveal the structural details of curved membranes, their dynamics pose an experimental challenge. Curvature especially affects the diffusion of lipids and macromolecules, yet neither experiments nor continuum models distinguish geometric effects from those caused by curvature-induced changes in membrane properties. Molecular simulations could excel here, yet despite community interest toward curved membranes, tools for their analysis are still lacking. Here, we satisfy this demand by introducing CurD, our novel and openly available implementation of the Vertex-oriented Triangle Propagation algorithm to the study of lipid diffusion along membranes with mean and/or Gaussian curvature. This approach, aided by our highly optimized implementation, computes geodetic distances significantly faster than conventional implementations of path-finding algorithms. Our tool, applied to coarse-grained simulations, allows for the first time the analysis of curvature effects on diffusion at size scales relevant to physiological processes such as endocytosis. Our analyses with different membrane geometries reveal that Gaussian curvature plays a surprisingly small role on lipid motion, whereas mean curvature; i.e., the packing of lipid headgroups largely dictates their mobility.
Collapse
Affiliation(s)
- Balázs Fábián
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo nám. 542/2, CZ-16000 Prague 6, Czech Republic
| | - Matti Javanainen
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo nám. 542/2, CZ-16000 Prague 6, Czech Republic
- Institute of Biotechnology, University of Helsinki, FI-00790 Helsinki, Finland
| |
Collapse
|
4
|
Kulshrestha A, Punnathanam SN, Roy R, Ayappa KG. Cholesterol catalyzes unfolding in membrane-inserted motifs of the pore forming protein cytolysin A. Biophys J 2023; 122:4068-4081. [PMID: 37740492 PMCID: PMC10598289 DOI: 10.1016/j.bpj.2023.09.005] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2023] [Revised: 08/20/2023] [Accepted: 09/13/2023] [Indexed: 09/24/2023] Open
Abstract
Plasma membrane-induced protein folding and conformational transitions play a central role in cellular homeostasis. Several transmembrane proteins are folded in the complex lipid milieu to acquire a specific structure and function. Bacterial pore forming toxins (PFTs) are proteins expressed by a large class of pathogenic bacteria that exploit the plasma membrane environment to efficiently undergo secondary structure changes, oligomerize, and form transmembrane pores. Unregulated pore formation causes ion imbalance, leading to cell death and infection. Determining the free energy landscape of these membrane-driven-driven transitions remains a challenging problem. Although cholesterol recognition is required for lytic activity of several proteins in the PFT family of toxins, the regulatory role of cholesterol for the α-PFT, cytolysin A expressed by Escherichia coli remains unexplained. In a recent free energy computation, we showed that the β tongue, a critical membrane-inserted motif of the ClyA toxin, has an on-pathway partially unfolded intermediate that refolds into the helix-turn-helix motif of the pore state. To understand the molecular role played by cholesterol, we carry out string-method-based computations in membranes devoid of cholesterol, which reveals an increase of ∼30 times in the free energy barrier for the loss of β sheet secondary structure when compared with membranes containing cholesterol. Specifically, the tyrosine-cholesterol interaction was found to be critical to creating the unfolded intermediate. Cholesterol also increases the packing and hydrophobicity of the bilayer, resulting in enhanced interactions of the bound protein before complete membrane insertion. Our study illustrates that cholesterol is critical to catalyzing and stabilizing the membrane-inserted unfolded state of the β tongue motif of ClyA, opening up fresh insights into cholesterol-assisted unfolding of membrane proteins.
Collapse
Affiliation(s)
- Avijeet Kulshrestha
- Department of Chemical Engineering, Indian Institute of Science, Bangalore, Karnataka, India
| | - Sudeep N Punnathanam
- Department of Chemical Engineering, Indian Institute of Science, Bangalore, Karnataka, India
| | - Rahul Roy
- Department of Chemical Engineering, Indian Institute of Science, Bangalore, Karnataka, India; Center for BioSystems Science and Engineering, Indian Institute of Science, Bangalore, Karnataka, India
| | - K Ganapathy Ayappa
- Department of Chemical Engineering, Indian Institute of Science, Bangalore, Karnataka, India.
| |
Collapse
|
5
|
Yesylevskyy S, Martinez-Seara H, Jungwirth P. Curvature Matters: Modeling Calcium Binding to Neutral and Anionic Phospholipid Bilayers. J Phys Chem B 2023; 127:4523-4531. [PMID: 37191140 DOI: 10.1021/acs.jpcb.3c01962] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/17/2023]
Abstract
In this work, the influence of membrane curvature on the Ca2+ binding to phospholipid bilayers is investigated by means of molecular dynamics simulations. In particular, we compared Ca2+ binding to flat, elastically buckled, or uniformly bent zwitterionic and anionic phospholipid bilayers. We demonstrate that Ca2+ ions bind preferably to the concave membrane surfaces in both types of bilayers. We also show that the membrane curvature leads to pronounced changes in Ca2+ binding including differences in free ion concentrations, lipid coordination distributions, and the patterns of ion binding to different chemical groups of lipids. Moreover, these effects differ substantially for the concave and convex membrane monolayers. Comparison between force fields with either full or scaled charges indicates that charge scaling results in reduction of the Ca2+ binding to curved phosphatidylserine bilayers, while for phosphatidylcholine membranes, calcium binds only weakly for both force fields.
Collapse
Affiliation(s)
- Semen Yesylevskyy
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo náměstí 542/2, 160 00 Prague 6, Czech Republic
- Department of Physics of Biological Systems, Institute of Physics of the National Academy of Sciences of Ukraine, Nauky Avenue 46, 03038 Kyiv, Ukraine
- Receptor.AI Incorporated, 20-22 Wenlock Road, N1 7GU London, U.K
| | - Hector Martinez-Seara
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo náměstí 542/2, 160 00 Prague 6, Czech Republic
| | - Pavel Jungwirth
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo náměstí 542/2, 160 00 Prague 6, Czech Republic
| |
Collapse
|
6
|
Blanco-Fernández G, Blanco-Fernandez B, Fernández-Ferreiro A, Otero-Espinar FJ. Lipidic lyotropic liquid crystals: Insights on biomedical applications. Adv Colloid Interface Sci 2023; 313:102867. [PMID: 36889183 DOI: 10.1016/j.cis.2023.102867] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2022] [Revised: 02/26/2023] [Accepted: 02/26/2023] [Indexed: 03/04/2023]
Abstract
Liquid crystals (LCs) possess unique physicochemical properties, translatable into a wide range of applications. To date, lipidic lyotropic LCs (LLCs) have been extensively explored in drug delivery and imaging owing to the capability to encapsulate and release payloads with different characteristics. The current landscape of lipidic LLCs in biomedical applications is provided in this review. Initially, the main properties, types, methods of fabrication and applications of LCs are showcased. Then, a comprehensive discussion of the main biomedical applications of lipidic LLCs accordingly to the application (drug and biomacromolecule delivery, tissue engineering and molecular imaging) and route of administration is examined. Further discussion of the main limitations and perspectives of lipidic LLCs in biomedical applications are also provided. STATEMENT OF SIGNIFICANCE: Liquid crystals (LCs) are those systems between a solid and liquid state that possess unique morphological and physicochemical properties, translatable into a wide range of biomedical applications. A short description of the properties of LCs, their types and manufacturing procedures is given to serve as a background to the topic. Then, the latest and most innovative research in the field of biomedicine is examined, specifically the areas of drug and biomacromolecule delivery, tissue engineering and molecular imaging. Finally, prospects of LCs in biomedicine are discussed to show future trends and perspectives that might be utilized. This article is an ampliation, improvement and actualization of our previous short forum article "Bringing lipidic lyotropic liquid crystal technology into biomedicine" published in TIPS.
Collapse
Affiliation(s)
- Guillermo Blanco-Fernández
- Pharmacology, Pharmacy and Pharmaceutical Technology Department, Faculty of Pharmacy, University of Santiago de Compostela (USC), Santiago de Compostela, Spain; Paraquasil Group, Health Research Institute of Santiago de Compostela (FIDIS), Santiago de Compostela, Spain; Institute of Materials (iMATUS), University of Santiago de Compostela (USC), Santiago de Compostela, Spain
| | - Bárbara Blanco-Fernandez
- CIBER in Bioengineering, Biomaterials and Nanomedicine, CIBER-BBN, Madrid, Spain; Institute for Bioengineering of Catalonia (IBEC), The Barcelona Institute of Science and Technology, Baldiri Reixac 10-12, Barcelona 08028, Spain.
| | - Anxo Fernández-Ferreiro
- Pharmacology Group, Health Research Institute of Santiago de Compostela (FIDIS), Santiago de Compostela, Spain; Pharmacy Department, University Clinical Hospital of Santiago de Compostela (SERGAS), Santiago de Compostela, Spain.
| | - Francisco J Otero-Espinar
- Pharmacology, Pharmacy and Pharmaceutical Technology Department, Faculty of Pharmacy, University of Santiago de Compostela (USC), Santiago de Compostela, Spain; Paraquasil Group, Health Research Institute of Santiago de Compostela (FIDIS), Santiago de Compostela, Spain; Institute of Materials (iMATUS), University of Santiago de Compostela (USC), Santiago de Compostela, Spain.
| |
Collapse
|
7
|
Krisanova N, Pozdnyakova N, Pastukhov A, Dudarenko M, Shatursky O, Gnatyuk O, Afonina U, Pyrshev K, Dovbeshko G, Yesylevskyy S, Borisova T. Amphiphilic anti-SARS-CoV-2 drug remdesivir incorporates into the lipid bilayer and nerve terminal membranes influencing excitatory and inhibitory neurotransmission. BIOCHIMICA ET BIOPHYSICA ACTA (BBA) - BIOMEMBRANES 2022; 1864:183945. [PMID: 35461828 PMCID: PMC9023372 DOI: 10.1016/j.bbamem.2022.183945] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/02/2022] [Revised: 04/03/2022] [Accepted: 04/17/2022] [Indexed: 12/05/2022]
Abstract
Remdesivir is a novel antiviral drug, which is active against the SARS-CoV-2 virus. Remdesivir is known to accumulate in the brain but it is not clear whether it influences the neurotransmission. Here we report diverse and pronounced effects of remdesivir on transportation and release of excitatory and inhibitory neurotransmitters in rat cortex nerve terminals (synaptosomes) in vitro. Direct incorporation of remdesivir molecules into the cellular membranes was shown by FTIR spectroscopy, planar phospholipid bilayer membranes and computational techniques. Remdesivir decreases depolarization-induced exocytotic release of L-[14C] glutamate and [3H] GABA, and also [3H] GABA uptake and extracellular level in synaptosomes in a dose-dependent manner. Fluorimetric studies confirmed remdesivir-induced impairment of exocytosis in nerve terminals and revealed a decrease in synaptic vesicle acidification. Our data suggest that remdesivir dosing during antiviral therapy should be precisely controlled to prevent possible neuromodulatory action at the presynaptic level. Further studies of neurotropic and membranotropic effects of remdesivir are necessary.
Collapse
|
8
|
Larsen AH. Molecular Dynamics Simulations of Curved Lipid Membranes. Int J Mol Sci 2022; 23:ijms23158098. [PMID: 35897670 PMCID: PMC9331392 DOI: 10.3390/ijms23158098] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2022] [Revised: 07/19/2022] [Accepted: 07/20/2022] [Indexed: 02/04/2023] Open
Abstract
Eukaryotic cells contain membranes with various curvatures, from the near-plane plasma membrane to the highly curved membranes of organelles, vesicles, and membrane protrusions. These curvatures are generated and sustained by curvature-inducing proteins, peptides, and lipids, and describing these mechanisms is an important scientific challenge. In addition to that, some molecules can sense membrane curvature and thereby be trafficked to specific locations. The description of curvature sensing is another fundamental challenge. Curved lipid membranes and their interplay with membrane-associated proteins can be investigated with molecular dynamics (MD) simulations. Various methods for simulating curved membranes with MD are discussed here, including tools for setting up simulation of vesicles and methods for sustaining membrane curvature. The latter are divided into methods that exploit scaffolding virtual beads, methods that use curvature-inducing molecules, and methods applying virtual forces. The variety of simulation tools allow researcher to closely match the conditions of experimental studies of membrane curvatures.
Collapse
|
9
|
Recent developments in membrane curvature sensing and induction by proteins. Biochim Biophys Acta Gen Subj 2021; 1865:129971. [PMID: 34333084 DOI: 10.1016/j.bbagen.2021.129971] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2021] [Revised: 07/11/2021] [Accepted: 07/25/2021] [Indexed: 12/22/2022]
Abstract
BACKGROUND Membrane-bound intracellular organelles have characteristic shapes attributed to different local membrane curvatures, and these attributes are conserved across species. Over the past decade, it has been confirmed that specific proteins control the large curvatures of the membrane, whereas many others due to their specific structural features can sense the curvatures and bind to the specific geometrical cues. Elucidating the interplay between sensing and induction is indispensable to understand the mechanisms behind various biological processes such as vesicular trafficking and budding. SCOPE OF REVIEW We provide an overview of major classes of membrane proteins and the mechanisms of curvature sensing and induction. We then discuss the importance of membrane elastic characteristics to induce the membrane shapes similar to intracellular organelles. Finally, we survey recently available assays developed for studying the curvature sensing and induction by many proteins. MAJOR CONCLUSIONS Recent theoretical/computational modeling along with experimental studies have uncovered fascinating connections between lipid membrane and protein interactions. However, the phenomena of protein localization and synchronization to generate spatiotemporal dynamics in membrane morphology are yet to be fully understood. GENERAL SIGNIFICANCE The understanding of protein-membrane interactions is essential to shed light on various biological processes. This further enables the technological applications of many natural proteins/peptides in therapeutic treatments. The studies of membrane dynamic shapes help to understand the fundamental functions of membranes, while the medicinal roles of various macromolecules (such as proteins, peptides, etc.) are being increasingly investigated.
Collapse
|
10
|
Yesylevskyy S, Khandelia H. EnCurv: Simple Technique of Maintaining Global Membrane Curvature in Molecular Dynamics Simulations. J Chem Theory Comput 2021; 17:1181-1193. [DOI: 10.1021/acs.jctc.0c00800] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Semen Yesylevskyy
- Department of Physics of Biological Systems, Institute of Physics of the National Academy of Sciences of Ukraine, Prospect Nauky 46, 03028 Kyiv, Ukraine
- Laboratoire Chrono Environnement UMR CNRS 6249, Université de Bourgogne Franche-Comté, 16 route de Gray, 25030 Besançon Cedex, France
| | - Himanshu Khandelia
- PHYLIFE: Physical Life Science, Department of Physics Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, Odense 5230 M, Denmark
- MEMPHYS: Center for Biomembrane Physics,
| |
Collapse
|
11
|
Aranha MP, Mukherjee D, Petridis L, Khomami B. An Atomistic Molecular Dynamics Study of Titanium Dioxide Adhesion to Lipid Bilayers. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2020; 36:1043-1052. [PMID: 31944772 DOI: 10.1021/acs.langmuir.9b03075] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
Titanium dioxide (TiO2) nanoparticles are found in an array of consumer and industrial products, and human exposure to these nanoparticles involves interaction with biological membranes. To understand the effect of the membrane lipid composition on bilayer perturbation by TiO2, we performed all-atom molecular dynamics simulations of nanosized TiO2 interacting with three single component bilayers differing only in their headgroup composition: the zwitterionic DOPC, which is overall neutral containing negatively charged phosphate and positively charged choline in its head, DOPG, which is overall anionic containing negatively charged phosphate and neutral glycerol, and the anionic DOPS, containing negatively charged phosphate attached to the hydroxyl side-chain of the amino acid, serine containing negatively charged carboxyl and positively charged ammonium. The nanoparticle adheres to all three bilayers causing a negative curvature on their top leaflet. However, the local deformation of DOPG was more pronounced than DOPC and DOPS. The anionic DOPG, which is the thinnest of the three bilayers, interacted most strongly with the TiO2. DOPS has the next strongest interaction; however, its high bending modulus enables it to resist deformation by the nanoparticle. DOPC has the weakest interaction with the nanoparticle of the three as it has the highest bending modulus and its zwitterionic head groups have strong cohesive interactions. We also observed a nonuniform response of the bilayers: the orientational order of the lipids near the nanoparticle decreases, while that of the lipids away from the nanoparticle increases. The overall thickness and bending modulus of DOPG increased upon contact with the nanoparticle owing to overall stiffening of the bilayer despite local softening, while the average structural and mechanical properties of DOPC and DOPS remain unchanged, which can be explained in part by the greater bilayer bending elasticicty of DOPC and DOPS. The above findings suggest that regions of biological membranes populated by anionic lipids with weaker bending elasticity will be more susceptible to perturbation by TiO2 nanoparticles than zwitterionic-rich regions.
Collapse
Affiliation(s)
- Michelle P Aranha
- Department of Biochemistry and Cellular and Molecular Biology , University of Tennessee 1311 Cumberland Ave , Knoxville , Tennessee 37916 , United States
- UT/ORNL Center for Molecular Biophysics , Oak Ridge National Laboratory , 1 Bethel Valley Road , Oak Ridge , Tennessee 37830 , United States
| | - Dibyendu Mukherjee
- Department of Chemical and Biomolecular Engineering , University of Tennessee , 1512 Middle Dr , Knoxville , Tennessee 37996 , United States
| | - Loukas Petridis
- Department of Biochemistry and Cellular and Molecular Biology , University of Tennessee 1311 Cumberland Ave , Knoxville , Tennessee 37916 , United States
- UT/ORNL Center for Molecular Biophysics , Oak Ridge National Laboratory , 1 Bethel Valley Road , Oak Ridge , Tennessee 37830 , United States
| | - Bamin Khomami
- Department of Chemical and Biomolecular Engineering , University of Tennessee , 1512 Middle Dr , Knoxville , Tennessee 37996 , United States
| |
Collapse
|
12
|
Yesylevskyy S, Rivel T, Ramseyer C. Curvature increases permeability of the plasma membrane for ions, water and the anti-cancer drugs cisplatin and gemcitabine. Sci Rep 2019; 9:17214. [PMID: 31748538 PMCID: PMC6868207 DOI: 10.1038/s41598-019-53952-2] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2019] [Accepted: 09/26/2019] [Indexed: 01/15/2023] Open
Abstract
In this work the permeability of a model asymmetric plasma membrane, for ions, water and the anti-cancer drugs cisplatin and gemcitabine is studied by means of all-atom molecular dynamics simulations. It is shown for the first time that permeability of the highly curved membrane increases from one to three orders of magnitude upon membrane bending depending on the compound and the sign of curvature. Our results suggest that the membrane curvature could be an important factor of drug translocation through the membrane.
Collapse
Affiliation(s)
- Semen Yesylevskyy
- Laboratoire Chrono Environnement UMR CNRS 6249, Université de Bourgogne Franche-Comté, 16 route de Gray, 25030, Besançon, Cedex, France.
- Department of Physics of Biological Systems, Institute of Physics of the National Academy of Sciences of Ukraine, Prospect Nauky 46, 03028, Kyiv, Ukraine.
| | - Timothée Rivel
- Laboratoire Chrono Environnement UMR CNRS 6249, Université de Bourgogne Franche-Comté, 16 route de Gray, 25030, Besançon, Cedex, France
| | - Christophe Ramseyer
- Laboratoire Chrono Environnement UMR CNRS 6249, Université de Bourgogne Franche-Comté, 16 route de Gray, 25030, Besançon, Cedex, France
| |
Collapse
|
13
|
Bouvier B. Curvature as a Collective Coordinate in Enhanced Sampling Membrane Simulations. J Chem Theory Comput 2019; 15:6551-6561. [DOI: 10.1021/acs.jctc.9b00716] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Benjamin Bouvier
- Laboratoire de Glycochimie, des Antimicrobiens et des Agroressources, CNRS UMR7378/Université de Picardie Jules Verne, 10, rue Baudelocque, 80039 Amiens Cedex, France
| |
Collapse
|
14
|
Santos DES, Pontes FJS, Lins RD, Coutinho K, Soares TA. SuAVE: A Tool for Analyzing Curvature-Dependent Properties in Chemical Interfaces. J Chem Inf Model 2019; 60:473-484. [DOI: 10.1021/acs.jcim.9b00569] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Denys E. S. Santos
- Departmento Química Fundamental, Universidade Federal de Pernambuco, 50740-560 Cidade Universitária, Recife, Brazil
| | - Frederico J. S. Pontes
- Departmento Química Fundamental, Universidade Federal de Pernambuco, 50740-560 Cidade Universitária, Recife, Brazil
| | - Roberto D. Lins
- Instituto Aggeu Magalhães, Fundação Oswaldo Cruz, 50740-465 Cidade Universitária, Recife, Brazil
| | - Kaline Coutinho
- Instituto de Física, Universidade de São Paulo, 05508-090 Cidade Universitária, São Paulo, Brazil
| | - Thereza A. Soares
- Departmento Química Fundamental, Universidade Federal de Pernambuco, 50740-560 Cidade Universitária, Recife, Brazil
| |
Collapse
|
15
|
Rózsa ZB, Németh LJ, Jójárt B, Nehéz K, Viskolcz B, Szőri M. Molecular Dynamics and Metadynamics Insights of 1,4-Dioxane-Induced Structural Changes of Biomembrane Models. J Phys Chem B 2019; 123:7869-7884. [DOI: 10.1021/acs.jpcb.9b04313] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Zsófia Borbála Rózsa
- Institute of Chemistry, University of Miskolc, Egyetemváros A/2, H-3515 Miskolc, Hungary
| | - Lukács József Németh
- Department of Food Engineering, Faculty of Engineering, University of Szeged, Mars tér 7, 6724 Szeged, Hungary
| | - Balázs Jójárt
- Department of Food Engineering, Faculty of Engineering, University of Szeged, Mars tér 7, 6724 Szeged, Hungary
| | - Károly Nehéz
- Department of Information Engineering, University of Miskolc, Miskolc-Egyetemváros Informatics Building, H-3515 Miskolc, Hungary
| | - Béla Viskolcz
- Institute of Chemistry, University of Miskolc, Egyetemváros A/2, H-3515 Miskolc, Hungary
| | - Milán Szőri
- Institute of Chemistry, University of Miskolc, Egyetemváros A/2, H-3515 Miskolc, Hungary
| |
Collapse
|
16
|
Rivel T, Ramseyer C, Yesylevskyy S. The asymmetry of plasma membranes and their cholesterol content influence the uptake of cisplatin. Sci Rep 2019; 9:5627. [PMID: 30948733 PMCID: PMC6449338 DOI: 10.1038/s41598-019-41903-w] [Citation(s) in RCA: 81] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2018] [Accepted: 03/18/2019] [Indexed: 12/20/2022] Open
Abstract
The composition of the plasma membrane of malignant cells is thought to influence the cellular uptake of cisplatin and to take part in developing resistance to this widespread anti-cancer drug. In this work we study the permeation of cisplatin through the model membranes of normal and cancer cells using molecular dynamics simulations. A special attention is paid to lipid asymmetry and cholesterol content of the membranes. The loss of lipid asymmetry, which is common for cancer cells, leads to a decrease in their permeability to cisplatin by one order of magnitude in comparison to the membranes of normal cells. The change in the cholesterol molar ratio from 0% to 33% also decreases the permeability of the membrane by approximately one order of magnitude. The permeability of pure DOPC membrane is 5-6 orders of magnitude higher than one of the membranes with realistic lipid composition, which makes it as an inadequate model for the studies of drug permeability.
Collapse
Affiliation(s)
- Timothée Rivel
- Laboratoire Chrono Environnement UMR CNRS 6249, Université de Bourgogne Franche-Comté, 16 route de Gray, 25030, Besançon, Cedex, France.
| | - Christophe Ramseyer
- Laboratoire Chrono Environnement UMR CNRS 6249, Université de Bourgogne Franche-Comté, 16 route de Gray, 25030, Besançon, Cedex, France
| | - Semen Yesylevskyy
- Laboratoire Chrono Environnement UMR CNRS 6249, Université de Bourgogne Franche-Comté, 16 route de Gray, 25030, Besançon, Cedex, France.,Department of Physics of Biological Systems, Institute of Physics of the National Academy of Sciences of Ukraine, Prospect Nauky 46, 03028, Kyiv, Ukraine
| |
Collapse
|
17
|
Sega M, Hantal G, Fábián B, Jedlovszky P. Pytim: A python package for the interfacial analysis of molecular simulations. J Comput Chem 2018; 39:2118-2125. [PMID: 30306571 PMCID: PMC6221047 DOI: 10.1002/jcc.25384] [Citation(s) in RCA: 49] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2018] [Revised: 05/31/2018] [Accepted: 06/02/2018] [Indexed: 11/08/2022]
Abstract
Pytim is a versatile python framework for the analysis of interfacial properties in molecular simulations. The code implements several algorithms for the identification of instantaneous interfaces of arbitrary shape, and analysis tools written specifically for the study of interfacial properties, such as intrinsic profiles. The code is written in the python language, and makes use of the numpy and scipy packages to deliver high computational performances. Pytim relies on the MDAnalysis library to analyze the trajectory file formats of popular simulation packages such as gromacs, charmm, namd, lammps or Amber, and can be used to steer OpenMM simulations. Pytim can write information about surfaces and surface atomic layers to vtk, cube, and pdb files for easy visualization. The classes of Pytim can be easily customized and extended to include new interfacial algorithms or analysis tools. The code is available as open source and is free of charge. © 2018 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Marcello Sega
- Faculty of Physics, University of Vienna, Boltzmangasse 5, Vienna A-1090, Austria
| | - György Hantal
- Faculty of Physics, University of Vienna, Boltzmangasse 5, Vienna A-1090, Austria
| | - Balázs Fábián
- Institut UTINAM (CNRS UMR 6213), Université Bourgogne Franche-Comté, 16 Route de Gray, Besançon Cedex F-25030, France.,Department of Inorganic and Analytical Chemistry, Budapest University of Technology and Economics, Szt. Gellért tér 4, Budapest H1111, Hungary
| | - Pál Jedlovszky
- Department of Chemistry, Eszterházy Károly University, Leányka utca 6, H-3300, Eger, Hungary
| |
Collapse
|
18
|
Grossman D, Sharon E, Katzav E. Shape and fluctuations of positively curved ribbons. Phys Rev E 2018; 98:022502. [PMID: 30253468 DOI: 10.1103/physreve.98.022502] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2018] [Indexed: 11/07/2022]
Abstract
We study the shape and shape fluctuations of incompatible, positively curved ribbons, with a flat reference metric and a spherelike reference curvature. Such incompatible geometry is likely to occur in many self-assembled materials and other experimental systems. Ribbons of this geometry exhibit a sharp transition between a rigid ring and an anomalously soft spring as a function of their width. As a result, the temperature dependence of these ribbons' shape is unique, exhibiting a nonmonotonic dependence of the persistence and Kuhn lengths on the temperature and width. We map the possible configuration phase space and show the existence of three phases: At high temperatures it is the ideal chain phase, where the ribbon is well described by classical models (e.g., wormlike chain model). The second phase, for cold and narrow ribbons, is the plane ergodic phase; a ribbon in this phase might be thought of as made out of segments that gyrate within an oblate spheroid with extreme aspect ratio. The third phase, for cold, wide ribbons, is a direct result of the residual stress caused by the incompatibility, called the random structured phase. A ribbon in this phase behaves on large scales as an ideal chain. However, the segments of this chain are not straight; rather they may have different shapes, mainly helices (both left and right handed) of various pitches.
Collapse
Affiliation(s)
- Doron Grossman
- Racah Institute of Physics, Hebrew University, Jerusalem 91904, Israel
| | - Eran Sharon
- Racah Institute of Physics, Hebrew University, Jerusalem 91904, Israel
| | - Eytan Katzav
- Racah Institute of Physics, Hebrew University, Jerusalem 91904, Israel
| |
Collapse
|
19
|
Bi H, Wang X, Han X, Voïtchovsky K. Impact of Electric Fields on the Nanoscale Behavior of Lipid Monolayers at the Surface of Graphite in Solution. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2018; 34:9561-9571. [PMID: 30028144 DOI: 10.1021/acs.langmuir.8b01631] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
The nanoscale organization and dynamics of lipid molecules in self-assembled membranes is central to the biological function of cells and in the technological development of synthetic lipid structures as well as in devices such as biosensors. Here, we explore the nanoscale molecular arrangement and dynamics of lipids assembled in monolayers at the surface of highly ordered pyrolytic graphite (HOPG), in different ionic solutions, and under electrical potentials. Using a combination of atomic force microscopy and fluorescence recovery after photobleaching, we show that HOPG is able to support fully formed and fluid lipid membranes, but mesoscale order and corrugations can be observed depending on the type of the lipid considered (1,2-dioleoyl- sn-glycero-3-phosphocholine, 1,2-dioleoyl- sn-glycero-3-phospho-l-serine (DOPS), and 1,2-dioleoyl-3-trimethylammoniumpropane) and the ion present (Na+, Ca2+, Cl-). Interfacial solvation forces and ion-specific effects dominate over the electrostatic changes induced by moderate electric fields (±1.0 V vs Ag/AgCl reference electrode) with particularly marked effects in the presence of calcium, and for DOPS. Our results provide insights into the interplay between the molecular, ionic, and electrostatic interactions and the formation of dynamical ordered structures in fluid lipid membranes.
Collapse
Affiliation(s)
- Hongmei Bi
- College of Science , Heilongjiang Bayi Agricultural University , Daqing 163319 , China
| | - Xuejing Wang
- State Key Laboratory of Urban Water Resource and Environment, School of Chemistry and Chemical Engineering , Harbin Institute of Technology , Harbin 150001 , China
| | - Xiaojun Han
- State Key Laboratory of Urban Water Resource and Environment, School of Chemistry and Chemical Engineering , Harbin Institute of Technology , Harbin 150001 , China
| | | |
Collapse
|
20
|
Urakami N, Jimbo T, Sakuma Y, Imai M. Molecular mechanism of vesicle division induced by coupling between lipid geometry and membrane curvatures. SOFT MATTER 2018. [PMID: 29517793 DOI: 10.1039/c7sm02188g] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/15/2023]
Abstract
We investigated the effects of lipid geometry on vesicle division using coarse grained molecular dynamics simulations. When the vesicle is composed of zero and negative spontaneous curvature lipids (ZSLs and NSLs), the difference in their molecular spontaneous curvatures destabilizes the neck of the limiting shape vesicle. In the vesicle division pathway, the neck developed into the stalk intermediates. The stalk was broken when the NSLs were expelled from the stalk. Free energy analysis shows that the coupling between the lipid geometry and the Gaussian rigidity is responsible for the observed vesicle division.
Collapse
Affiliation(s)
- Naohito Urakami
- Department of Physics and Informatics, Yamaguchi University, 1677-1 Yoshida, Yamaguchi, 753-8512, Japan.
| | | | | | | |
Collapse
|
21
|
Yesylevskyy SO, Rivel T, Ramseyer C. The influence of curvature on the properties of the plasma membrane. Insights from atomistic molecular dynamics simulations. Sci Rep 2017; 7:16078. [PMID: 29167583 PMCID: PMC5700167 DOI: 10.1038/s41598-017-16450-x] [Citation(s) in RCA: 57] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2017] [Accepted: 11/13/2017] [Indexed: 11/24/2022] Open
Abstract
In this work we conduct a systematic analysis of the influence of curvature on various properties of a realistic model of mammalian plasma membrane with asymmetric lipid content of monolayers and a realistic concentration of cholesterol. In order to do this we developed new technique for simulating membranes with the global membrane curvature restricted to any desirable value while keeping free lateral diffusion of the lipids and without introducing artifacts or perturbing the membrane structure. We show that the thickness of the hydrophobic core of the concave monolayer decreases by approximately 1.3 Å in comparison to that of the flat membrane, while the thickness of the convex monolayer does not change. The order parameter of the lipid tails decreases significantly in the certain layers of the curved membrane. The area per lipid increases in the convex monolayer and decreases in the concave monolayer. The cholesterol inclination angle decreases when the curvature of a particular monolayer changes from concave to convex. The amount of cholesterol in the minor fraction located between the membrane leaflets is zero in the membrane with positive curvature and increases to 1.7% in the flat membrane and to 2.5% in the membrane with negative curvature.
Collapse
Affiliation(s)
- Semen O Yesylevskyy
- Department of Physics of Biological Systems, Institute of Physics of the National Academy of Sciences of Ukraine, Prospect Nauky 46, Kyiv, 03680, Ukraine.
| | - Timothée Rivel
- Laboratoire Chrono Environnement UMR CNRS 6249, Université de Bourgogne Franche-Comté, 16 route de Gray, 25030, Besançon, Cedex, France
| | - Christophe Ramseyer
- Laboratoire Chrono Environnement UMR CNRS 6249, Université de Bourgogne Franche-Comté, 16 route de Gray, 25030, Besançon, Cedex, France
| |
Collapse
|
22
|
Cherniavskyi YK, Ramseyer C, Yesylevskyy SO. Interaction of C60 fullerenes with asymmetric and curved lipid membranes: a molecular dynamics study. Phys Chem Chem Phys 2015; 18:278-84. [PMID: 26608905 DOI: 10.1039/c5cp05838d] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Interaction of fullerenes with asymmetric and curved DOPC/DOPS bicelles is studied by means of coarse-grained molecular dynamics simulations. The effects caused by asymmetric lipid composition of the membrane leaflets and the curvature of the membrane are analyzed. It is shown that the aggregates of fullerenes prefer to penetrate into the membrane in the regions of the moderately positive mean curvature. Upon penetration into the hydrophobic core of the membrane fullerenes avoid the regions of the extreme positive or the negative curvature. Fullerenes increase the ordering of lipid tails, which are in direct contact with them, but do not influence other lipids significantly. Our data suggest that the effects of the membrane curvature should be taken into account in the studies concerning permeability of the membranes to fullerenes and fullerene-based drug delivery systems.
Collapse
Affiliation(s)
- Yevhen K Cherniavskyi
- Department of Physics of Biological Systems, Institute of Physics of the National Academy of Sciences of Ukraine, Prospect Nauky 46, Kiev-28, 03680, Ukraine.
| | | | | |
Collapse
|
23
|
Yesylevskyy SO. Pteros 2.0: Evolution of the fast parallel molecular analysis library for C++ and python. J Comput Chem 2015; 36:1480-8. [DOI: 10.1002/jcc.23943] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2015] [Revised: 04/26/2015] [Accepted: 04/27/2015] [Indexed: 11/12/2022]
Affiliation(s)
- Semen O. Yesylevskyy
- Department of Physics of Biological Systems; Institute of Physics, National Academy of Sciences of Ukraine; Prospect Nauki, 46, Kyiv-03039 Ukraine
| |
Collapse
|
24
|
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
|