1
|
Delgado-Coello B, Luna-Reyes I, Méndez-Acevedo KM, Bravo-Martínez J, Montalvan-Sorrosa D, Mas-Oliva J. Analysis of cholesterol-recognition motifs of the plasma membrane Ca 2+-ATPase. J Bioenerg Biomembr 2024; 56:205-219. [PMID: 38436904 PMCID: PMC11116186 DOI: 10.1007/s10863-024-10010-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2024] [Accepted: 02/26/2024] [Indexed: 03/05/2024]
Abstract
The plasma membrane Ca2+-ATPase (PMCA) is crucial for the fine tuning of intracellular calcium levels in eukaryotic cells. In this study, we show the presence of CARC sequences in all human and rat PMCA isoforms and we performed further analysis by molecular dynamics simulations. This analysis focuses on PMCA1, containing three CARC motifs, and PMCA4, with four CARC domains. In PMCA1, two CARC motifs reside within transmembrane domains, while the third is situated at the intracellular interface. The simulations depict more stable RMSD values and lower RMSF fluctuations in the presence of cholesterol, emphasizing its potential stabilizing effect. In PMCA4, a distinct dynamic was found. Notably, the total energy differences between simulations with cholesterol and phospholipids are pronounced in PMCA4 compared to PMCA1. RMSD values for PMCA4 indicate a more energetically favorable conformation in the presence of cholesterol, suggesting a robust interaction between CARCs and this lipid in the membranes. Furthermore, RMSF analysis for CARCs in both PMCA isoforms exhibit lower values in the presence of cholesterol compared to POPC alone. The analysis of H-bond occupancy and total energy values strongly suggests the potential interaction of CARCs with cholesterol. Given the crucial role of PMCAs in physiological calcium regulation and their involvement in diverse pathological processes, this study underscores the significance of CARC motifs and their interaction with cholesterol in elucidating PMCA function. These insights into the energetic preferences associated with CARC-cholesterol interactions offer valuable implications for understanding PMCA function in maintaining calcium homeostasis and addressing potential associated pathologies.
Collapse
Affiliation(s)
- Blanca Delgado-Coello
- Departamento de Bioquímica y Biología Estructural, Instituto de Fisiología Celular, Universidad Nacional Autónoma de México, Apdo. Postal 70-243, Ciudad de México, C.P. 04510, México.
| | - Ismael Luna-Reyes
- Departamento de Bioquímica y Biología Estructural, Instituto de Fisiología Celular, Universidad Nacional Autónoma de México, Apdo. Postal 70-243, Ciudad de México, C.P. 04510, México
| | - Kevin M Méndez-Acevedo
- Departamento de Bioquímica y Biología Estructural, Instituto de Fisiología Celular, Universidad Nacional Autónoma de México, Apdo. Postal 70-243, Ciudad de México, C.P. 04510, México
- University of Cambridge Metabolic Research Laboratories and NIHR Cambridge Biomedical Research Centre, Wellcome-MRC Institute of Metabolic Science, Addenbrooke's Hospital, Cambridge, UK
| | - Jorge Bravo-Martínez
- Departamento de Fisiología, Facultad de Medicina, Universidad Nacional Autónoma de México, Ciudad de México, México
| | - Danai Montalvan-Sorrosa
- Departamento de Química de Biomacromoléculas, Instituto de Química, Universidad Nacional Autónoma de México, Ciudad de México, México
| | - Jaime Mas-Oliva
- Departamento de Bioquímica y Biología Estructural, Instituto de Fisiología Celular, Universidad Nacional Autónoma de México, Apdo. Postal 70-243, Ciudad de México, C.P. 04510, México.
| |
Collapse
|
2
|
Chiu PL, Orjuela JD, de Groot BL, Aponte-Santamaría C, Walz T. Structure and dynamics of cholesterol-mediated aquaporin-0 arrays and implications for lipid rafts. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2023.05.16.540959. [PMID: 37292626 PMCID: PMC10245776 DOI: 10.1101/2023.05.16.540959] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Aquaporin-0 (AQP0) tetramers form square arrays in lens membranes through a yet unknown mechanism, but lens membranes are enriched in sphingomyelin and cholesterol. Here, we determined electron crystallographic structures of AQP0 in sphingomyelin/cholesterol membranes and performed molecular dynamics (MD) simulations to establish that the observed cholesterol positions represent those seen around an isolated AQP0 tetramer and that the AQP0 tetramer largely defines the location and orientation of most of its associated cholesterol molecules. At a high concentration, cholesterol increases the hydrophobic thickness of the annular lipid shell around AQP0 tetramers, which may thus cluster to mitigate the resulting hydrophobic mismatch. Moreover, neighboring AQP0 tetramers sandwich a cholesterol deep in the center of the membrane. MD simulations show that the association of two AQP0 tetramers is necessary to maintain the deep cholesterol in its position and that the deep cholesterol increases the force required to laterally detach two AQP0 tetramers, not only due to protein-protein contacts but also due to increased lipid-protein complementarity. Since each tetramer interacts with four such 'glue' cholesterols, avidity effects may stabilize larger arrays. The principles proposed to drive AQP0 array formation could also underlie protein clustering in lipid rafts.
Collapse
|
3
|
Nguyen ATP, Weigle AT, Shukla D. Functional regulation of aquaporin dynamics by lipid bilayer composition. Nat Commun 2024; 15:1848. [PMID: 38418487 PMCID: PMC10901782 DOI: 10.1038/s41467-024-46027-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2023] [Accepted: 02/12/2024] [Indexed: 03/01/2024] Open
Abstract
With the diversity of lipid-protein interactions, any observed membrane protein dynamics or functions directly depend on the lipid bilayer selection. However, the implications of lipid bilayer choice are seldom considered unless characteristic lipid-protein interactions have been previously reported. Using molecular dynamics simulation, we characterize the effects of membrane embedding on plant aquaporin SoPIP2;1, which has no reported high-affinity lipid interactions. The regulatory impacts of a realistic lipid bilayer, and nine different homogeneous bilayers, on varying SoPIP2;1 dynamics are examined. We demonstrate that SoPIP2;1's structure, thermodynamics, kinetics, and water transport are altered as a function of each membrane construct's ensemble properties. Notably, the realistic bilayer provides stabilization of non-functional SoPIP2;1 metastable states. Hydrophobic mismatch and lipid order parameter calculations further explain how lipid ensemble properties manipulate SoPIP2;1 behavior. Our results illustrate the importance of careful bilayer selection when studying membrane proteins. To this end, we advise cautionary measures when performing membrane protein molecular dynamics simulations.
Collapse
Affiliation(s)
- Anh T P Nguyen
- Department of Chemical and Biomolecular Engineering, University of Illinois Urbana-Champaign, Urbana, IL, 61801, USA
| | - Austin T Weigle
- Department of Chemistry, University of Illinois Urbana-Champaign, Urbana, IL, 61801, USA
| | - Diwakar Shukla
- Department of Chemical and Biomolecular Engineering, University of Illinois Urbana-Champaign, Urbana, IL, 61801, USA.
- Center for Biophysics and Quantitative Biology, University of Illinois Urbana-Champaign, Urbana, IL, 61801, USA.
- Department of Bioengineering, University of Illinois Urbana-Champaign, Urbana, IL, 61801, USA.
- Department of Plant Biology, University of Illinois Urbana-Champaign, Urbana, IL, 61801, USA.
| |
Collapse
|
4
|
Wang N, Zhang Y, Wang W, Ye Z, Chen H, Hu G, Ouyang D. How can machine learning and multiscale modeling benefit ocular drug development? Adv Drug Deliv Rev 2023; 196:114772. [PMID: 36906232 DOI: 10.1016/j.addr.2023.114772] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2022] [Revised: 02/06/2023] [Accepted: 03/05/2023] [Indexed: 03/12/2023]
Abstract
The eyes possess sophisticated physiological structures, diverse disease targets, limited drug delivery space, distinctive barriers, and complicated biomechanical processes, requiring a more in-depth understanding of the interactions between drug delivery systems and biological systems for ocular formulation development. However, the tiny size of the eyes makes sampling difficult and invasive studies costly and ethically constrained. Developing ocular formulations following conventional trial-and-error formulation and manufacturing process screening procedures is inefficient. Along with the popularity of computational pharmaceutics, non-invasive in silico modeling & simulation offer new opportunities for the paradigm shift of ocular formulation development. The current work first systematically reviews the theoretical underpinnings, advanced applications, and unique advantages of data-driven machine learning and multiscale simulation approaches represented by molecular simulation, mathematical modeling, and pharmacokinetic (PK)/pharmacodynamic (PD) modeling for ocular drug development. Following this, a new computer-driven framework for rational pharmaceutical formulation design is proposed, inspired by the potential of in silico explorations in understanding drug delivery details and facilitating drug formulation design. Lastly, to promote the paradigm shift, integrated in silico methodologies were highlighted, and discussions on data challenges, model practicality, personalized modeling, regulatory science, interdisciplinary collaboration, and talent training were conducted in detail with a view to achieving more efficient objective-oriented pharmaceutical formulation design.
Collapse
Affiliation(s)
- Nannan Wang
- State Key Laboratory of Quality Research in Chinese Medicine, Institute of Chinese Medical Sciences (ICMS), University of Macau, Macau, China
| | - Yunsen Zhang
- State Key Laboratory of Quality Research in Chinese Medicine, Institute of Chinese Medical Sciences (ICMS), University of Macau, Macau, China
| | - Wei Wang
- State Key Laboratory of Quality Research in Chinese Medicine, Institute of Chinese Medical Sciences (ICMS), University of Macau, Macau, China
| | - Zhuyifan Ye
- State Key Laboratory of Quality Research in Chinese Medicine, Institute of Chinese Medical Sciences (ICMS), University of Macau, Macau, China
| | - Hongyu Chen
- State Key Laboratory of Quality Research in Chinese Medicine, Institute of Chinese Medical Sciences (ICMS), University of Macau, Macau, China; Faculty of Science and Technology (FST), University of Macau, Macau, China
| | - Guanghui Hu
- Faculty of Science and Technology (FST), University of Macau, Macau, China
| | - Defang Ouyang
- State Key Laboratory of Quality Research in Chinese Medicine, Institute of Chinese Medical Sciences (ICMS), University of Macau, Macau, China; Department of Public Health and Medicinal Administration, Faculty of Health Sciences (FHS), University of Macau, Macau, China.
| |
Collapse
|
5
|
Cino EA, Borbuliak M, Hu S, Tieleman DP. Lipid distributions and transleaflet cholesterol migration near heterogeneous surfaces in asymmetric bilayers. Faraday Discuss 2021; 232:103-113. [PMID: 34549760 DOI: 10.1039/d1fd00003a] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Specific and nonspecific protein-lipid interactions in cell membranes have important roles in an abundance of biological functions. We have used coarse-grained (CG) molecular dynamics (MD) simulations to assess lipid distributions and cholesterol flipping dynamics around surfaces in a model asymmetric plasma membrane containing one of six structurally distinct entities: aquaporin-1 (AQP1), the bacterial β-barrel outer membrane proteins OmpF and OmpX, the KcsA potassium channel, the WALP23 peptide and a carbon nanotube (CNT). Our findings revealed varied lipid partitioning and cholesterol flipping times around the different solutes and putative cholesterol binding sites in AQP1 and KcsA. The results suggest that protein-lipid interactions can be highly variable, and that surface-dependent lipid profiles are effectively manifested in CG simulations with the Martini force field.
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. .,Department of Biochemistry and Immunology, Federal University of Minas Gerais, Belo Horizonte, MG, Brazil
| | - Mariia Borbuliak
- Centre for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta, T2N 1N4, Canada.
| | - Shangnong Hu
- 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
|
Smith AK, Khayat E, Lockhart C, Klimov DK. Do Cholesterol and Sphingomyelin Change the Mechanism of Aβ 25-35 Peptide Binding to Zwitterionic Bilayer? J Chem Inf Model 2019; 59:5207-5217. [PMID: 31738555 DOI: 10.1021/acs.jcim.9b00763] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Using replica exchange with solute tempering all-atom molecular dynamics, we studied the equilibrium binding of Aβ25-35 peptide to the ternary bilayer composed of an equimolar mixture of dimyristoylphosphatidylcholine (DMPC), N-palmitoylsphingomyelin (PSM), and cholesterol. Binding of the same peptide to the pure DMPC bilayer served as a control. Due to significant C-terminal hydrophobic moment, binding to the ternary and DMPC bilayers promotes helical structure in the peptide. For both bilayers a polarized binding profile is observed, in which the N-terminus anchors to the bilayer surface, whereas the C-terminus alternates between unbound and inserted states. Both ternary and DMPC bilayers feature two Aβ25-35 bound states, surface bound, S, and inserted, I, separated by modest free energy barriers. Experimental data are in agreement with our results but indicate that cholesterol impact is Aβ fragment dependent. For Aβ25-35, we predict that its binding mechanism is independent of the inclusion of PSM and cholesterol into the bilayer.
Collapse
Affiliation(s)
- Amy K Smith
- School of Systems Biology , George Mason University , Manassas , Virginia 20110 , United States
| | - Elias Khayat
- School of Systems Biology , George Mason University , Manassas , Virginia 20110 , United States
| | - Christopher Lockhart
- School of Systems Biology , George Mason University , Manassas , Virginia 20110 , United States
| | - Dmitri K Klimov
- School of Systems Biology , George Mason University , Manassas , Virginia 20110 , United States
| |
Collapse
|
7
|
Corradi V, Sejdiu BI, Mesa-Galloso H, Abdizadeh H, Noskov SY, Marrink SJ, Tieleman DP. Emerging Diversity in Lipid-Protein Interactions. Chem Rev 2019; 119:5775-5848. [PMID: 30758191 PMCID: PMC6509647 DOI: 10.1021/acs.chemrev.8b00451] [Citation(s) in RCA: 245] [Impact Index Per Article: 49.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2018] [Indexed: 02/07/2023]
Abstract
Membrane lipids interact with proteins in a variety of ways, ranging from providing a stable membrane environment for proteins to being embedded in to detailed roles in complicated and well-regulated protein functions. Experimental and computational advances are converging in a rapidly expanding research area of lipid-protein interactions. Experimentally, the database of high-resolution membrane protein structures is growing, as are capabilities to identify the complex lipid composition of different membranes, to probe the challenging time and length scales of lipid-protein interactions, and to link lipid-protein interactions to protein function in a variety of proteins. Computationally, more accurate membrane models and more powerful computers now enable a detailed look at lipid-protein interactions and increasing overlap with experimental observations for validation and joint interpretation of simulation and experiment. Here we review papers that use computational approaches to study detailed lipid-protein interactions, together with brief experimental and physiological contexts, aiming at comprehensive coverage of simulation papers in the last five years. Overall, a complex picture of lipid-protein interactions emerges, through a range of mechanisms including modulation of the physical properties of the lipid environment, detailed chemical interactions between lipids and proteins, and key functional roles of very specific lipids binding to well-defined binding sites on proteins. Computationally, despite important limitations, molecular dynamics simulations with current computer power and theoretical models are now in an excellent position to answer detailed questions about lipid-protein interactions.
Collapse
Affiliation(s)
- Valentina Corradi
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Besian I. Sejdiu
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Haydee Mesa-Galloso
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Haleh Abdizadeh
- Groningen
Biomolecular Sciences and Biotechnology Institute and Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Sergei Yu. Noskov
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Siewert J. Marrink
- Groningen
Biomolecular Sciences and Biotechnology Institute and Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - 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
|
8
|
Enkavi G, Javanainen M, Kulig W, Róg T, Vattulainen I. Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance. Chem Rev 2019; 119:5607-5774. [PMID: 30859819 PMCID: PMC6727218 DOI: 10.1021/acs.chemrev.8b00538] [Citation(s) in RCA: 175] [Impact Index Per Article: 35.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2018] [Indexed: 12/23/2022]
Abstract
Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.
Collapse
Affiliation(s)
- Giray Enkavi
- Department
of Physics, University of
Helsinki, P.O. Box 64, FI-00014 Helsinki, Finland
| | - Matti Javanainen
- Department
of Physics, University of
Helsinki, P.O. Box 64, FI-00014 Helsinki, Finland
- Institute
of Organic Chemistry and Biochemistry of the Czech Academy
of Sciences, Flemingovo naḿesti 542/2, 16610 Prague, Czech Republic
- Computational
Physics Laboratory, Tampere University, P.O. Box 692, FI-33014 Tampere, Finland
| | - Waldemar Kulig
- Department
of Physics, University of
Helsinki, P.O. Box 64, FI-00014 Helsinki, Finland
| | - Tomasz Róg
- Department
of Physics, University of
Helsinki, P.O. Box 64, FI-00014 Helsinki, Finland
- Computational
Physics Laboratory, Tampere University, P.O. Box 692, FI-33014 Tampere, Finland
| | - Ilpo Vattulainen
- Department
of Physics, University of
Helsinki, P.O. Box 64, FI-00014 Helsinki, Finland
- Computational
Physics Laboratory, Tampere University, P.O. Box 692, FI-33014 Tampere, Finland
- MEMPHYS-Center
for Biomembrane Physics
| |
Collapse
|
9
|
Muller MP, Jiang T, Sun C, Lihan M, Pant S, Mahinthichaichan P, Trifan A, Tajkhorshid E. Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation. Chem Rev 2019; 119:6086-6161. [PMID: 30978005 PMCID: PMC6506392 DOI: 10.1021/acs.chemrev.8b00608] [Citation(s) in RCA: 133] [Impact Index Per Article: 26.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.
Collapse
Affiliation(s)
- Melanie P. Muller
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- College of Medicine
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Tao Jiang
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Chang Sun
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Muyun Lihan
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Shashank Pant
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Paween Mahinthichaichan
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Anda Trifan
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Emad Tajkhorshid
- NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology
- Department of Biochemistry
- Center for Biophysics and Quantitative Biology
- College of Medicine
- University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| |
Collapse
|
10
|
Hall JE, Freites JA, Tobias DJ. Experimental and Simulation Studies of Aquaporin 0 Water Permeability and Regulation. Chem Rev 2019; 119:6015-6039. [PMID: 31026155 DOI: 10.1021/acs.chemrev.9b00106] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
We begin with the history of aquaporin zero (AQP0), the most prevalent membrane protein in the eye lens, from the early days when AQP0 was a protein of unknown function known as Major Intrinsic Protein 26. We progress through its joining the aquaporin family as a water channel in its own right and discuss how regulation of its water permeability by pH and calcium came to be discovered experimentally and linked to lens homeostasis and development. We review the development of molecular dynamics (MD) simulations of lipid bilayers and membrane proteins, including aquaporins, with an emphasis on simulation studies that have elucidated the mechanisms of water conduction, selectivity, and proton exclusion by aquaporins in general. We also review experimental and theoretical progress toward understanding why mammalian AQP0 has a lower water permeability than other aquaporins and the evolution of our present understanding of how its water permeability is regulated by pH and calcium. Finally, we discuss how MD simulations have elucidated the nature of lipid interactions with AQP0.
Collapse
|
11
|
Leonard AN, Wang E, Monje-Galvan V, Klauda JB. Developing and Testing of Lipid Force Fields with Applications to Modeling Cellular Membranes. Chem Rev 2019; 119:6227-6269. [DOI: 10.1021/acs.chemrev.8b00384] [Citation(s) in RCA: 53] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
12
|
Meleleo D, Sblano C. Influence of cholesterol on human calcitonin channel formation. Possible role of sterol as molecular chaperone. AIMS BIOPHYSICS 2019. [DOI: 10.3934/biophy.2019.1.23] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
|
13
|
Klauda JB. Perspective: Computational modeling of accurate cellular membranes with molecular resolution. J Chem Phys 2018; 149:220901. [DOI: 10.1063/1.5055007] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Affiliation(s)
- Jeffery B. Klauda
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, Maryland 20742, USA
- Biophysics Graduate Program, University of Maryland, College Park, Maryland 20742, USA
| |
Collapse
|
14
|
Smith AK, Klimov DK. Molecular Dynamics Investigation of the Ternary Bilayer Formed by Saturated Phosphotidylcholine, Sphingomyelin, and Cholesterol. J Phys Chem B 2018; 122:11311-11325. [DOI: 10.1021/acs.jpcb.8b07256] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Affiliation(s)
- Amy K. Smith
- School of Systems Biology, George Mason University, Manassas, Virginia 20110, United States
| | - Dmitri K. Klimov
- School of Systems Biology, George Mason University, Manassas, Virginia 20110, United States
| |
Collapse
|
15
|
Khakbaz P, Klauda JB. Investigation of phase transitions of saturated phosphocholine lipid bilayers via molecular dynamics simulations. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2018; 1860:1489-1501. [DOI: 10.1016/j.bbamem.2018.04.014] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/21/2017] [Revised: 04/26/2018] [Accepted: 04/26/2018] [Indexed: 12/01/2022]
|
16
|
Stereospecific Interactions of Cholesterol in a Model Cell Membrane: Implications for the Membrane Dipole Potential. J Membr Biol 2018; 251:507-519. [DOI: 10.1007/s00232-018-0016-0] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2017] [Accepted: 01/25/2018] [Indexed: 12/11/2022]
|
17
|
Adams M, Wang E, Zhuang X, Klauda JB. Simulations of simple Bovine and Homo sapiens outer cortex ocular lens membrane models with a majority concentration of cholesterol. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2017; 1860:2134-2144. [PMID: 29169746 DOI: 10.1016/j.bbamem.2017.11.010] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/09/2017] [Revised: 11/06/2017] [Accepted: 11/17/2017] [Indexed: 12/31/2022]
Abstract
The lipid composition of bovine and human ocular lens membranes has been probed, and a variety of lipids have been found including phosphatidylcholine (PC), phosphatidylethanolamine (PE), sphingomyelin (SM), and cholesterol (CHOL) with cholesterol being present in particularly high concentrations. In this study, we use the all-atom CHARMM36 force field to simulate binary, ternary, and quaternary mixtures as models of the ocular lens. High concentration of cholesterol, in combination with different and varying diversity of phospholipids (PL) and sphingolipids (SL), affect the structure of the ocular lens lipid bilayer. The following analyses were done for each simulation: surface area per lipid, component surface area per lipid, deuterium order parameters (SCD), electron density profiles (EDP), membrane thickness, hydrogen bonding, radial distribution functions, clustering, and sterol tilt angle distribution. The SCD show significant bilayer alignment and packing in cholesterol-rich bilayers. The EDP show the transition from liquid crystalline to liquid ordered with the addition of cholesterol. Hydrogen bonds in our systems show the tendency for intramolecular interactions between cholesterol and fully saturated lipid tails for less complex bilayers. But with an increased number of components in the bilayer, the acyl chain of the lipids becomes a less important characteristic, and the headgroup of the lipid becomes more significant. Overall, cholesterol is the driving force of membrane structure of the ocular lens membrane where interactions between cholesterol, PL, and SL determine structure and function of the biomembrane. The goal of this work is to develop a baseline for further study of more physiologically realistic ocular lens lipid membranes. This article is part of a Special Issue entitled: Emergence of Complex Behavior in Biomembranes edited by Marjorie Longo.
Collapse
Affiliation(s)
- Mark Adams
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, MD 20742, USA
| | - Eric Wang
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, MD 20742, USA
| | - Xiaohong Zhuang
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, MD 20742, USA
| | - Jeffery B Klauda
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, MD 20742, USA; Biophysics Program, University of Maryland, College Park, MD 20742, USA.
| |
Collapse
|
18
|
Monje-Galvan V, Klauda JB. Two sterols, two bilayers: insights on membrane structure from molecular dynamics. MOLECULAR SIMULATION 2017. [DOI: 10.1080/08927022.2017.1353690] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Viviana Monje-Galvan
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, MD, USA
| | - Jeffery B. Klauda
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, MD, USA
- Biophysics Program, University of Maryland, College Park, MD, USA
| |
Collapse
|
19
|
Bandara A, Panahi A, Pantelopulos GA, Straub JE. Exploring the structure and stability of cholesterol dimer formation in multicomponent lipid bilayers. J Comput Chem 2017; 38:1479-1488. [PMID: 27761918 PMCID: PMC5398962 DOI: 10.1002/jcc.24516] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2016] [Revised: 09/27/2016] [Accepted: 09/30/2016] [Indexed: 01/13/2023]
Abstract
For 40 years, the existence and possible functional importance of cholesterol dimer formation has been discussed. Due to challenges associated with structural studies of membrane lipids, there has as yet been no direct experimental verification of the existence and relevance of the cholesterol dimer. Building on recent advances in lipid force fields for molecular simulation, in this work the structure and stability of the cholesterol dimer is characterized in POPC bilayers in absence and presence of sphingomyelin. The cholesterol dimer structural ensemble is found to consist of sub-states that reflect, but also differ from, previously proposed dimer structures. While face-to-face dimer structures predominate, no evidence is found for the existence of tail-to-tail dimers in POPC lipid bilayers. Near stoichiometric complex formation of cholesterol with sphingomyelin is found to effect cholesterol dimer structure without impacting population. Comparison with NMR-derived order parameters provide validation for the simulation model employed and conclusions drawn related to the structure and stability of cholesterol dimers in multicomponent lipid bilayers. © 2016 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Asanga Bandara
- Department of Chemistry, Boston University, 590 Commonwealth Ave., Boston, Massachusetts, 02215, United States
| | - Afra Panahi
- Department of Chemistry, Boston University, 590 Commonwealth Ave., Boston, Massachusetts, 02215, United States
| | - George A. Pantelopulos
- Department of Chemistry, Boston University, 590 Commonwealth Ave., Boston, Massachusetts, 02215, United States
| | - John E. Straub
- Department of Chemistry, Boston University, 590 Commonwealth Ave., Boston, Massachusetts, 02215, United States
| |
Collapse
|
20
|
Tong J, Wu Z, Briggs MM, Schulten K, McIntosh TJ. The Water Permeability and Pore Entrance Structure of Aquaporin-4 Depend on Lipid Bilayer Thickness. Biophys J 2017; 111:90-9. [PMID: 27410737 DOI: 10.1016/j.bpj.2016.05.039] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2016] [Accepted: 05/13/2016] [Indexed: 12/22/2022] Open
Abstract
Aquaporin-4 (AQP4), the primary water channel in glial cells of the mammalian brain, plays a critical role in water transport in the central nervous system. Previous experiments have shown that the water permeability of AQP4 depends on the cholesterol content in the lipid bilayer, but it was not clear whether changes in permeability were due to direct cholesterol-AQP4 interactions or to indirect effects caused by cholesterol-induced changes in bilayer elasticity or bilayer thickness. To determine the effects resulting only from bilayer thickness, here we use a combination of experiments and simulations to analyze AQP4 in cholesterol-free phospholipid bilayers with similar elastic properties but different hydrocarbon core thicknesses previously determined by x-ray diffraction. The channel (unit) water permeabilities of AQP4 measured by osmotic-gradient experiments were 3.5 ± 0.2 × 10(-13) cm(3)/s (mean ± SE), 3.0 ± 0.3 × 10(-13) cm(3)/s, 2.5 ± 0.2 × 10(-13) cm(3)/s, and 0.9 ± 0.1 × 10(-13) cm(3)/s in bilayers containing (C22:1)(C22:1)PC, (C20:1)(C20:1)PC, (C16:0)(C18:1)PC, and (C13:0)(C13:0)PC, respectively. Channel permeabilities obtained by molecular dynamics (MD) simulations were 3.3 ± 0.1 × 10(-13) cm(3)/s and 2.5 ± 0.1 × 10(-13) cm(3)/s in (C22:1)(C22:1)PC and (C14:0)(C14:0)PC bilayers, respectively. Both the osmotic-gradient and MD-simulation results indicated that AQP4 channel permeability decreased with decreasing bilayer hydrocarbon thickness. The MD simulations also suggested structural modifications in AQP4 in response to changes in bilayer thickness. Although the simulations showed no appreciable changes to the radius of the pore located in the hydrocarbon region of the bilayers, the simulations indicated that there were changes in both pore length and α-helix organization near the cytoplasmic vestibule of the channel. These structural changes, caused by mismatch between the hydrophobic length of AQP4 and the bilayer hydrocarbon thickness, could explain the observed differences in water permeability with changes in bilayer thickness.
Collapse
Affiliation(s)
- Jihong Tong
- Department of Cell Biology, Duke University Medical Center, Durham, North Carolina
| | - Zhe Wu
- Center for the Physics of Living Cells and Beckman Institute, University of Illinois Urbana-Champaign, Urbana, Illinois
| | - Margaret M Briggs
- Department of Cell Biology, Duke University Medical Center, Durham, North Carolina
| | - Klaus Schulten
- Center for the Physics of Living Cells and Beckman Institute, University of Illinois Urbana-Champaign, Urbana, Illinois.
| | - Thomas J McIntosh
- Department of Cell Biology, Duke University Medical Center, Durham, North Carolina.
| |
Collapse
|
21
|
Zhuang X, Ou A, Klauda JB. Simulations of simple linoleic acid-containing lipid membranes and models for the soybean plasma membranes. J Chem Phys 2017; 146:215103. [PMID: 28595398 DOI: 10.1063/1.4983655] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The all-atom CHARMM36 lipid force field (C36FF) has been tested with saturated, monounsaturated, and polyunsaturated lipids; however, it has not been validated against the 18:2 linoleoyl lipids with an unsaturated sn-1 chain. The linoleoyl lipids are common in plants and the main component of the soybean membrane. The lipid composition of soybean plasma membranes has been thoroughly characterized with experimental studies. However, there is comparatively less work done with computational modeling. Our molecular dynamics (MD) simulation results show that the pure linoleoyl lipids, 1-stearoyl-2-linoleoyl-sn-glycero-3-phosphocholine (18:0/18:2) and 1,2-dilinoleoyl-sn-glycero-3-phosphocholine (di-18:2), agree very well with the experiments, which demonstrates the accuracy of the C36FF for the computational study of soybean membranes. Based on the experimental composition, the soybean hypocotyl and root plasma membrane models are developed with each containing seven or eight types of linoleoyl phospholipids and two types of sterols (sitosterol and stigmasterol). MD simulations are performed to characterize soybean membranes, and the hydrogen bonds and clustering results demonstrate that the lipids prefer to interact with the lipids of the same/similar tail unsaturation. All the results suggest that these two soybean membrane models can be used as a basis for further research in soybean and higher plant membranes involving membrane-associated proteins.
Collapse
Affiliation(s)
- Xiaohong Zhuang
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, Maryland 20742, USA
| | - Anna Ou
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, Maryland 20742, USA
| | - Jeffery B Klauda
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, Maryland 20742, USA
| |
Collapse
|
22
|
Bera I, Klauda JB. Molecular Simulations of Mixed Lipid Bilayers with Sphingomyelin, Glycerophospholipids, and Cholesterol. J Phys Chem B 2017; 121:5197-5208. [DOI: 10.1021/acs.jpcb.7b00359] [Citation(s) in RCA: 42] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
- Indrani Bera
- Department
of Chemical and Biomolecular Engineering and ‡Biophysics Program, University of Maryland, College Park, Maryland 20742, United States
| | - Jeffery B. Klauda
- Department
of Chemical and Biomolecular Engineering and ‡Biophysics Program, University of Maryland, College Park, Maryland 20742, United States
| |
Collapse
|
23
|
Wang E, Klauda JB. Examination of Mixtures Containing Sphingomyelin and Cholesterol by Molecular Dynamics Simulations. J Phys Chem B 2017; 121:4833-4844. [DOI: 10.1021/acs.jpcb.7b01832] [Citation(s) in RCA: 38] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Affiliation(s)
- Eric Wang
- Department
of Chemical and Biomolecular Engineering and ‡Biophysics Graduate Program, University of Maryland, College Park, Maryland 20742, United States
| | - Jeffery B. Klauda
- Department
of Chemical and Biomolecular Engineering and ‡Biophysics Graduate Program, University of Maryland, College Park, Maryland 20742, United States
| |
Collapse
|
24
|
Briones R, Aponte-Santamaría C, de Groot BL. Localization and Ordering of Lipids Around Aquaporin-0: Protein and Lipid Mobility Effects. Front Physiol 2017; 8:124. [PMID: 28303107 PMCID: PMC5332469 DOI: 10.3389/fphys.2017.00124] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2016] [Accepted: 02/15/2017] [Indexed: 11/13/2022] Open
Abstract
Hydrophobic matching, lipid sorting, and protein oligomerization are key principles by which lipids and proteins organize in biological membranes. The Aquaporin-0 channel (AQP0), solved by electron crystallography (EC) at cryogenic temperatures, is one of the few protein-lipid complexes of which the structure is available in atomic detail. EC and room-temperature molecular dynamics (MD) of dimyristoylglycerophosphocholine (DMPC) annular lipids around AQP0 show similarities, however, crystal-packing and temperature might affect the protein surface or the lipids distribution. To understand the role of temperature, lipid phase, and protein mobility in the localization and ordering of AQP0-lipids, we used MD simulations of an AQP0-DMPC bilayer system. Simulations were performed at physiological and at DMPC gel-phase temperatures. To decouple the protein and lipid mobility effects, we induced gel-phase in the lipids or restrained the protein. We monitored the lipid ordering effects around the protein. Reducing the system temperature or inducing lipid gel-phase had a marginal effect on the annular lipid localization. However, restraining the protein mobility increased the annular lipid localization around the whole AQP0 surface, resembling EC. The distribution of the inter-phosphate and hydrophobic thicknesses showed that stretching of the DMPC annular layer around AQP0 surface is the mechanism that compensates the hydrophobic mismatch in this system. The distribution of the local area-per-lipid and the acyl-chain order parameters showed particular fluid- and gel-like areas that involved several lipid layers. These areas were in contact with the surfaces of higher and lower protein mobility, respectively. We conclude that the AQP0 surfaces induce specific fluid- and gel-phase prone areas. The presence of these areas might guide the AQP0 lipid sorting interactions with other membrane components, and is compatible with the squared array oligomerization of AQP0 tetramers separated by a layer of annular lipids.
Collapse
Affiliation(s)
- Rodolfo Briones
- Computational Biomolecular Dynamics Group, Max Planck Institute for Biophysical Chemistry Göttingen, Germany
| | - Camilo Aponte-Santamaría
- Molecular Biomechanics Group, Heidelberg Institute for Theoretical Studies and Interdisciplinary Center for Scientific Computing Heidelberg, Germany
| | - Bert L de Groot
- Computational Biomolecular Dynamics Group, Max Planck Institute for Biophysical Chemistry Göttingen, Germany
| |
Collapse
|
25
|
Ermilova I, Lyubartsev AP. Extension of the Slipids Force Field to Polyunsaturated Lipids. J Phys Chem B 2016; 120:12826-12842. [DOI: 10.1021/acs.jpcb.6b05422] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Affiliation(s)
- Inna Ermilova
- Department of Materials and
Environmental Chemistry, Stockholm University, SE 106 91 Stockholm, Sweden
| | - Alexander P. Lyubartsev
- Department of Materials and
Environmental Chemistry, Stockholm University, SE 106 91 Stockholm, Sweden
| |
Collapse
|
26
|
Boughter CT, Monje-Galvan V, Im W, Klauda JB. Influence of Cholesterol on Phospholipid Bilayer Structure and Dynamics. J Phys Chem B 2016; 120:11761-11772. [DOI: 10.1021/acs.jpcb.6b08574] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Christopher T. Boughter
- Department
of Chemical and Biomolecular Engineering, University of Maryland, College
Park, Maryland 20742, United States
| | - Viviana Monje-Galvan
- Department
of Chemical and Biomolecular Engineering, University of Maryland, College
Park, Maryland 20742, United States
| | - Wonpil Im
- Department
of Biological Sciences and Bioengineering Program, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Jeffery B. Klauda
- Department
of Chemical and Biomolecular Engineering, University of Maryland, College
Park, Maryland 20742, United States
- Biophysics
Program, University of Maryland, College Park, Maryland 20742, United States
| |
Collapse
|
27
|
Ghumman A, Robinson D. Modulation of the Turkey β1-Adrenergic Receptor by Membrane Rafts - Insight from Molecular Dynamics. ChemistrySelect 2016. [DOI: 10.1002/slct.201600849] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Afshan Ghumman
- School of Chemistry; University Of Nottingham; University Park Nottingham NG7 2RD UK
| | - David Robinson
- School of Chemistry; University Of Nottingham; University Park Nottingham NG7 2RD UK
| |
Collapse
|
28
|
Khakbaz P, Klauda JB. Probing the importance of lipid diversity in cell membranes via molecular simulation. Chem Phys Lipids 2015; 192:12-22. [PMID: 26260616 DOI: 10.1016/j.chemphyslip.2015.08.003] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2015] [Revised: 08/04/2015] [Accepted: 08/05/2015] [Indexed: 12/31/2022]
Abstract
Lipid membranes in prokaryotes and eukaryotes have a wide array of lipids that are necessary for proper membrane structure and function. In this paper, an introduction to lipid diversity in biology and a mini-review on how molecular simulations have been used to model biological membranes (primarily limited to one to three lipid types in most simulation-based models) is provided, which motivates the use of all-atom molecular dynamics (MD) simulations to study the effect of lipid diversity on properties of realistic membrane models of prokaryotes and eukaryotes. As an example, cytoplasmic membrane models of Escherichia coli were developed at different stages of the colony growth cycle (early-log, mid-log, stationary and overnight). The main difference between lipid compositions at each stage was the concentration of a cyclopropane-containing moiety on the sn-2 lipid acyl chain (cyC17:0). Triplicate MD simulations for each stage were run for 300 ns to study the influence of lipid diversity on the surface area per lipid, area compressibility modulus, deuterium order parameters, and electron density profiles. The overnight stage (also known as the death stage) had the highest average surface area per lipid, highest rigidity, and lowest bilayer thickness compare to other stages of E. coli cytoplasmic membrane. Although bilayer thickness did depend on the growth stage, the changes between these were small suggesting that the hydrophobic core of transmembrane proteins fit well with the membrane in all growth stages. Although it is still common practise in MD simulations of membrane proteins to use simple one- or two-component membranes, it can be important to use diverse lipid model membranes when membrane protein structure and function are influenced by changes in lipid membrane composition.
Collapse
Affiliation(s)
- Pouyan Khakbaz
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, MD 20742, USA
| | - Jeffery B Klauda
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, MD 20742, USA; Biophysics Program, University of Maryland, College Park, MD 20742, USA.
| |
Collapse
|
29
|
Grouleff J, Irudayam SJ, Skeby KK, Schiøtt B. The influence of cholesterol on membrane protein structure, function, and dynamics studied by molecular dynamics simulations. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2015; 1848:1783-95. [PMID: 25839353 DOI: 10.1016/j.bbamem.2015.03.029] [Citation(s) in RCA: 123] [Impact Index Per Article: 13.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/19/2014] [Revised: 03/24/2015] [Accepted: 03/25/2015] [Indexed: 12/12/2022]
Abstract
The plasma membrane, which encapsulates human cells, is composed of a complex mixture of lipids and embedded proteins. Emerging knowledge points towards the lipids as having a regulating role in protein function. Furthermore, insight from protein crystallography has revealed several different types of lipids intimately bound to membrane proteins and peptides, hereby possibly pointing to a site of action for the observed regulation. Cholesterol is among the lipid membrane constituents most often observed to be co-crystallized with membrane proteins, and the cholesterol levels in cell membranes have been found to play an essential role in health and disease. Remarkably little is known about the mechanism of lipid regulation of membrane protein function in health as well as in disease. Herein, we review molecular dynamics simulation studies aimed at investigating the effect of cholesterol on membrane protein and peptide properties. This article is part of a Special Issue entitled: Lipid-protein interactions.
Collapse
Affiliation(s)
- Julie Grouleff
- Department of Chemistry and Interdisciplinary Nanoscience Center, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
| | - Sheeba Jem Irudayam
- Department of Chemistry and Interdisciplinary Nanoscience Center, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
| | - Katrine K Skeby
- Department of Chemistry and Interdisciplinary Nanoscience Center, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
| | - Birgit Schiøtt
- Department of Chemistry and Interdisciplinary Nanoscience Center, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark.
| |
Collapse
|
30
|
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
|
31
|
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
|
32
|
Chaban VV, Khandelia H. Distribution of neutral lipids in the lipid droplet core. J Phys Chem B 2014; 118:11145-51. [PMID: 25188363 DOI: 10.1021/jp506693d] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
Cholesteryl esters (CEs) are a form of cholesterol (CHOL) storage in the living cells, as opposed to free CHOL. CEs are major constituents of low density lipoprotein particles. Therefore, CEs are implicated in provoking atherosclerosis. Arranged into cytoplasmic lipid droplets (LDs), CEs are stored intracellularly. They can also be transported extracellularly by means of lipoproteins. In this work, large-scale molecular dynamics (MD) simulations are used to characterize the molecular structure of LDs containing various fractions (10-50 mol %) of cholesteryl oleate (CO) with respect to triolein (TO) fraction. The simulated LDs were covered by a phospholipid monolayer formed by a mixture of 1-palmitoyl-2-oleoylphosphatidylcholine, POPC (75 mol %), and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine, POPE (25 mol %), molecules. We report that most CO molecules are located within the hydrophobic core of LDs, whereas a small fraction (0.3-1.9 mol %) penetrates the monolayer. The solubility of CO in the phospholipid monolayer is relatively small. Due to a good miscibility with TO molecules, CO forms a liquid phase inside LD at 333 K. There is long-range order in the liquid TO-CO droplet core up to 8 nm from the phospholipid interface, resulting from the structuring of hydrophilic groups. This structuring slowly decays in the direction toward the LD center of mass. No sorting of TO and CO is detected, irrespective of the molar fractions simulated. The distribution of CO within the LDs is significant in determining the rate of their hydrolysis by surface-bound enzyme lipases, and thus has a subsequent impact on the levels of CO in plasma and LDLs.
Collapse
Affiliation(s)
- Vitaly V Chaban
- MEMPHYS - Center for Biomembrane Physics, Syddansk Universitet , Odense M., 5230, Denmark
| | | |
Collapse
|
33
|
Blake HL, Robinson D. QM/MM studies of contemporary and novel membrane raft fluorescent probes. Molecules 2014; 19:10230-41. [PMID: 25029071 PMCID: PMC6271554 DOI: 10.3390/molecules190710230] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2014] [Revised: 07/09/2014] [Accepted: 07/10/2014] [Indexed: 11/22/2022] Open
Abstract
We have studied a number of contemporary and novel membrane probes, selected for their structural similarity to membrane raft components, in order to properly anchor themselves within a sphingolipid/cholesterol rich region. A QM/MM approach was adopted in order to understand the structural and electrostatic influences of fluorescence emission shifts of the probes in different lipid and solvation environments. The proposed modifications to the membrane probes have shown encouraging data relating not only to emission shifts within the membrane, but also their ability to anchor within a membrane raft domain and the stability to internalization within a membrane system.
Collapse
Affiliation(s)
- Hannah L Blake
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, UK.
| | - David Robinson
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, UK.
| |
Collapse
|
34
|
Robinson D. A Polarizable Force-Field for Cholesterol and Sphingomyelin. J Chem Theory Comput 2013; 9:2498-503. [DOI: 10.1021/ct400103e] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- David Robinson
- School of
Chemistry, University of Nottingham, University
Park, Nottingham, NG7 2RD, United Kingdom
| |
Collapse
|
35
|
Zhao Y, Vararattanavech A, Li X, Hélixnielsen C, Vissing T, Torres J, Wang R, Fane AG, Tang CY. Effects of proteoliposome composition and draw solution types on separation performance of aquaporin-based proteoliposomes: implications for seawater desalination using aquaporin-based biomimetic membranes. ENVIRONMENTAL SCIENCE & TECHNOLOGY 2013; 47:1496-1503. [PMID: 23311686 DOI: 10.1021/es304306t] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/01/2023]
Abstract
Aquaporins are a large family of water transport proteins in cell membranes. Their high water permeability and solute rejection make them potential building blocks for high-performance biomimetic membranes for desalination. In the current study, proteoliposomes were prepared using AquaporinZ from Escherichia coli cells, and their separation properties were characterized by stopped-flow measurements. The current study systematically investigated the effect of proteoliposome composition (lipid type, protein-to-lipid ratio (PLR), and the addition of cholesterol) on water permeability and NaCl retention. Among the various lipids investigated, 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC)-based proteoliposomes were found to have excellent osmotic water permeability and NaCl reflection coefficient values. Increasing the PLR of DOPC proteoliposomes up to 1:200 increased their osmotic water permeability. However, further increase in the PLR reduced the osmotic water permeability probably due to the occurrence of defects in the proteoliposomes, whereas the addition of cholesterol improved their osmotic water permeation likely due to defects sealing. The current study also investigated the effect of major dissolved ions in seawater (e.g., Mg(2+) and SO(4)(2-)) on the stability of proteoliposomes, and design criteria for aquaporin-based biomimetic membranes are proposed in the context of desalination.
Collapse
Affiliation(s)
- Yang Zhao
- Singapore Membrane Technology Centre, Nanyang Technological University, Singapore 639798
| | | | | | | | | | | | | | | | | |
Collapse
|
36
|
Lim JB, Rogaski B, Klauda JB. Update of the cholesterol force field parameters in CHARMM. J Phys Chem B 2011; 116:203-10. [PMID: 22136112 DOI: 10.1021/jp207925m] [Citation(s) in RCA: 156] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A modification of the CHARMM36 lipid force field (C36) for cholesterol, henceforth, called C36c, is reported. A fused ring compound, decalin, was used to model the steroid section of cholesterol. For decalin, C36 inaccurately predicts the heat of vaporization (~10 kJ/mol) and molar volume (~10 cc/mol), but C36c resulted in near perfect comparison with experiment. MD simulations of decalin and heptane at various compositions were run to estimate the enthalpy and volumes of mixing to compare to experiment for this simple model of cholesterol in a chain environment. Superior estimates for these thermodynamic properties were obtained with C36c versus C36. These new parameters were applied to cholesterol, and quantum mechanical calculations were performed to modify the torsional potential of an acyl chain torsion for cholesterol. This model was tested through simulations of DMPC/10% cholesterol, DMPC/30% cholesterol, and DOPC/10% cholesterol. The C36 and C36c results were similar for surface areas per lipid, deuterium order parameters, electron density profiles, and atomic form factors and generally agree well with experiment. However, C36 and C36c produced slightly different cholesterol angle distributions with C36c adopting a more perpendicular orientation with respect to the bilayer plane. The new parameters in the C36c modification should enable more accurate simulations of lipid bilayers with cholesterol, especially for those interested in the free energy of lipid flip/flop or transfer of phospholipids and/or cholesterol.
Collapse
Affiliation(s)
- Joseph B Lim
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, Maryland 20742, United States
| | | | | |
Collapse
|
37
|
Manna M, Mukhopadhyay C. Cholesterol driven alteration of the conformation and dynamics of phospholamban in model membranes. Phys Chem Chem Phys 2011; 13:20188-98. [DOI: 10.1039/c1cp21793c] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|