151
|
Rhys NH, Duffy IB, Sowden CL, Lorenz CD, McLain SE. On the hydration of DOPE in solution. J Chem Phys 2019; 150:115104. [PMID: 30902020 DOI: 10.1063/1.5085736] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
Abstract
The atomic-scale hydration structure around the 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine (DOPE) headgroup in a chloroform/water solution has been investigated using neutron diffraction enhanced by isotopic substitution and NMR, coupled with empirical potential structure refinement and molecular dynamics simulations. The results obtained show the preferential binding sites for water molecules on the DOPE headgroups, with the most predominant interactions being with the ammonium and phosphate groups. Interestingly, the level of hydration, as well as the association of DOPE molecules, varies according to the simulation method used. The results here suggest the presence of a tight water network around these lipid headgroups that could affect the permeability of the membrane for lipid-mediated diffusion.
Collapse
Affiliation(s)
- Natasha H Rhys
- Department of Biochemistry, University of Oxford, Oxford OX1 3QU, United Kingdom
| | - Imogen B Duffy
- Department of Biochemistry, University of Oxford, Oxford OX1 3QU, United Kingdom
| | - Christopher L Sowden
- Department of Biochemistry, University of Oxford, Oxford OX1 3QU, United Kingdom
| | - Christian D Lorenz
- Department of Physics, King's College London, London WC2R 2LS, United Kingdom
| | - Sylvia E McLain
- Department of Biochemistry, University of Oxford, Oxford OX1 3QU, United Kingdom
| |
Collapse
|
152
|
Nagle JF, Venable RM, Maroclo-Kemmerling E, Tristram-Nagle S, Harper PE, Pastor RW. Revisiting Volumes of Lipid Components in Bilayers. J Phys Chem B 2019; 123:2697-2709. [PMID: 30836006 DOI: 10.1021/acs.jpcb.8b12010] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In addition to obtaining the highly precise volumes of lipids in lipid bilayers, it has been desirable to obtain the volumes of parts of each lipid, such as the methylenes and terminal methyls on the hydrocarbon chains and the head group. Obtaining such component volumes from experiment and from simulations is re-examined, first by distinguishing methods based on apparent versus partial molar volumes. Although somewhat different, both these methods give results that are counterintuitive and that differ from results obtained by a more local method that can only be applied to simulations. These comparisons reveal differences in the average methylene component volume that result in larger differences in the head group component volumes. Literature experimental volume data for unsaturated phosphocholines and for alkanes have been used and new data have been acquired for saturated phosphocholines. Data and simulations cover extended ranges of temperature to assess both the temperature and chain length dependence of the component volumes. A new method to refine the determination of component volumes is proposed that uses experimental data for different chain lengths at temperatures guided by the temperature dependence determined in simulations. These refinements enable more precise comparisons of the component volumes of different lipids and alkanes in different phases. Finally, the notion of free volume is extended to components using the Lennard-Jones radii to estimate the excluded volume of each component. This analysis reveals that head group free volumes are relatively independent of thermodynamic phase, whereas both the methylene and methyl free volumes increase dramatically when bilayers transition from gel to fluid.
Collapse
Affiliation(s)
- John F Nagle
- Department of Physics , Carnegie Mellon University , Pittsburgh , Pennsylvania 15213 , United States
| | - Richard M Venable
- Laboratory of Computational Biology , National Heart, Lung, and Blood Institute, National Institutes of Health , Bethesda , Maryland 20892 , United States
| | | | - Stephanie Tristram-Nagle
- Department of Physics , Carnegie Mellon University , Pittsburgh , Pennsylvania 15213 , United States
| | - Paul E Harper
- Department of Physics & Astronomy , Calvin College , Grand Rapids , Michigan 49546 , United States
| | - Richard W Pastor
- Laboratory of Computational Biology , National Heart, Lung, and Blood Institute, National Institutes of Health , Bethesda , Maryland 20892 , United States
| |
Collapse
|
153
|
Covino R, Hummer G, Ernst R. Integrated Functions of Membrane Property Sensors and a Hidden Side of the Unfolded Protein Response. Mol Cell 2019; 71:458-467. [PMID: 30075144 DOI: 10.1016/j.molcel.2018.07.019] [Citation(s) in RCA: 38] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2018] [Revised: 07/09/2018] [Accepted: 07/17/2018] [Indexed: 12/26/2022]
Abstract
Eukaryotic cells face the challenge of maintaining the complex composition of several coexisting organelles. The molecular mechanisms underlying the homeostasis of subcellular membranes and their adaptation during stress are only now starting to emerge. Here, we discuss three membrane property sensors of the endoplasmic reticulum (ER), namely OPI1, MGA2, and IRE1, each controlling a large cellular program impacting the lipid metabolic network. OPI1 coordinates the production of membrane and storage lipids, MGA2 regulates the production of unsaturated fatty acids required for membrane biogenesis, and IRE1 controls the unfolded protein response (UPR) to adjust ER size, protein folding, and the secretory capacity of the cell. Although these proteins use remarkably distinct sensing mechanisms, they are functionally connected via the ER membrane and cooperate to maintain membrane homeostasis. As a rationalization of the recently described mechanism of UPR activation by lipid bilayer stress, we propose that IRE1 can sense the protein-to-lipid ratio in the ER membrane to ensure a balanced production of membrane proteins and lipids.
Collapse
Affiliation(s)
- Roberto Covino
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Max-von-Laue-Strasse 3, 60438 Frankfurt am Main, Germany
| | - Gerhard Hummer
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Max-von-Laue-Strasse 3, 60438 Frankfurt am Main, Germany; Institute of Biophysics, Goethe University, 60438 Frankfurt am Main, Germany
| | - Robert Ernst
- Department of Medical Biochemistry and Molecular Biology, Saarland University, Kirrberger Str. 100, Gebäude 61.4, 66421 Homburg, Germany.
| |
Collapse
|
154
|
Pak A, Dannenhoffer-Lafage T, Madsen JJ, Voth GA. Systematic Coarse-Grained Lipid Force Fields with Semiexplicit Solvation via Virtual Sites. J Chem Theory Comput 2019; 15:2087-2100. [PMID: 30702887 PMCID: PMC6416712 DOI: 10.1021/acs.jctc.8b01033] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2018] [Indexed: 12/15/2022]
Abstract
Despite the central role of lipids in many biophysical functions, the molecular mechanisms that dictate macroscopic lipid behavior remain elusive to both experimental and computational approaches. As such, there has been much interest in the development of low-resolution, implicit-solvent coarse-grained (CG) models to dynamically simulate biologically relevant spatiotemporal scales with molecular fidelity. However, in the absence of solvent, a key challenge for CG models is to faithfully emulate solvent-mediated forces, which include both hydrophilic and hydrophobic interactions that drive lipid aggregation and self-assembly. In this work, we provide a new methodological framework to incorporate semiexplicit solvent effects through the use of virtual CG particles, which represent structural features of the solvent-lipid interface. To do so, we leverage two systematic coarse-graining approaches, multiscale coarse-graining (MS-CG) and relative entropy minimization (REM), in a hybrid fashion to construct our virtual-site CG (VCG) models. As a proof-of-concept, we focus our efforts on two lipid species, 1,2-dioleoyl- sn-glycero-3-phosphocholine (DOPC) and 1,2-dipalmitoyl- sn-glycero-3-phosphocholine (DPPC), which adopt a liquid-disordered and gel phase, respectively, at room temperature. Through our analysis, we also present, to our knowledge, the first direct comparison between the MS-CG and REM methods for a complex biomolecule and highlight each of their strengths and weaknesses. We further demonstrate that VCG models recapitulate the rich biophysics of lipids, which enable self-assembly, morphological diversity, and multiple phases. Our findings suggest that the VCG framework is a powerful approach for investigation into macromolecular biophysics.
Collapse
Affiliation(s)
- Alexander
J. Pak
- Department of Chemistry, The
University of Chicago, Chicago, Illinois 60637, United States
| | | | - Jesper J. Madsen
- Department of Chemistry, The
University of Chicago, Chicago, Illinois 60637, United States
| | - Gregory A. Voth
- Department of Chemistry, The
University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
155
|
Shahane G, Ding W, Palaiokostas M, Orsi M. Physical properties of model biological lipid bilayers: insights from all-atom molecular dynamics simulations. J Mol Model 2019; 25:76. [PMID: 30806797 DOI: 10.1007/s00894-019-3964-0] [Citation(s) in RCA: 80] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2018] [Accepted: 02/07/2019] [Indexed: 02/04/2023]
Abstract
The physical properties of lipid bilayers are sensitive to the specific type and composition of the lipids that make up the many different types of cell membranes. Studying model bilayers of representative heterogeneous compositions can provide key insights into membrane functionality. In this work, we use atomistic molecular dynamics simulations to characterize key properties in a number of bilayer membranes of varying composition. We first examine basic properties, such as lipid area, volume, and bilayer thickness, of simple, homogeneous bilayers comprising several lipid types, which are prevalent in biological membranes. Such lipids are then used in simulations of heterogeneous systems representative of bacterial, mammalian, and cancer membranes. Our analysis is especially focused on depth-dependent, transmembrane profiles; in particular, we calculate lateral pressure and dipole potential profiles, two fundamental properties which play key roles in a large number of biological functions.
Collapse
Affiliation(s)
- Ganesh Shahane
- Institute of Bioengineering, Queen Mary University of London, Mile End Road, London, E1 4NS, UK
| | - Wei Ding
- School of Engineering & Materials Science, Queen Mary University of London, Mile End Road, London, E1 4NS, UK
| | - Michail Palaiokostas
- School of Engineering & Materials Science, Queen Mary University of London, Mile End Road, London, E1 4NS, UK
| | - Mario Orsi
- Department of Applied Sciences, University of the West of England, Coldharbour Lane, Bristol, BS16 1QY, UK.
| |
Collapse
|
156
|
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]
|
157
|
Leone V, Waclawska I, Kossmann K, Koshy C, Sharma M, Prisner TF, Ziegler C, Endeward B, Forrest LR. Interpretation of spectroscopic data using molecular simulations for the secondary active transporter BetP. J Gen Physiol 2019; 151:381-394. [PMID: 30728216 PMCID: PMC6400524 DOI: 10.1085/jgp.201812111] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2018] [Revised: 11/26/2018] [Accepted: 01/11/2019] [Indexed: 11/20/2022] Open
Abstract
Mechanistic understanding of dynamic membrane proteins such as transporters, receptors, and channels requires accurate depictions of conformational ensembles, and the manner in which they interchange as a function of environmental factors including substrates, lipids, and inhibitors. Spectroscopic techniques such as electron spin resonance (ESR) pulsed electron-electron double resonance (PELDOR), also known as double electron-electron resonance (DEER), provide a complement to atomistic structures obtained from x-ray crystallography or cryo-EM, since spectroscopic data reflect an ensemble and can be measured in more native solvents, unperturbed by a crystal lattice. However, attempts to interpret DEER data are frequently stymied by discrepancies with the structural data, which may arise due to differences in conditions, the dynamics of the protein, or the flexibility of the attached paramagnetic spin labels. Recently, molecular simulation techniques such as EBMetaD have been developed that create a conformational ensemble matching an experimental distance distribution while applying the minimal possible bias. Moreover, it has been proposed that the work required during an EBMetaD simulation to match an experimentally determined distribution could be used as a metric with which to assign conformational states to a given measurement. Here, we demonstrate the application of this concept for a sodium-coupled transport protein, BetP. Because the probe, protein, and lipid bilayer are all represented in atomic detail, the different contributions to the work, such as the extent of protein backbone movements, can be separated. This work therefore illustrates how ranking simulations based on EBMetaD can help to bridge the gap between structural and biophysical data and thereby enhance our understanding of membrane protein conformational mechanisms.
Collapse
Affiliation(s)
- Vanessa Leone
- Computational Structural Biology Section, National Institutes of Neurological Disorders and Stroke, National Institutes of Health, Bethesda, MD
| | | | - Katharina Kossmann
- Institute of Biophysics and Biophysical Chemistry, University of Regensburg, Regensburg, Germany
| | - Caroline Koshy
- Max Planck Institute for Biophysics, Frankfurt am Main, Germany
| | - Monika Sharma
- Computational Structural Biology Section, National Institutes of Neurological Disorders and Stroke, National Institutes of Health, Bethesda, MD
| | - Thomas F Prisner
- Institute of Physical and Theoretical Chemistry and Center of Biomolecular Magnetic Resonance, Goethe University, Frankfurt am Main, Germany
| | - Christine Ziegler
- Institute of Biophysics and Biophysical Chemistry, University of Regensburg, Regensburg, Germany
| | - Burkhard Endeward
- Institute of Physical and Theoretical Chemistry and Center of Biomolecular Magnetic Resonance, Goethe University, Frankfurt am Main, Germany
| | - Lucy R Forrest
- Computational Structural Biology Section, National Institutes of Neurological Disorders and Stroke, National Institutes of Health, Bethesda, MD
| |
Collapse
|
158
|
Smith DJ, Klauda JB, Sodt AJ. Simulation Best Practices for Lipid Membranes [Article v1.0]. LIVING JOURNAL OF COMPUTATIONAL MOLECULAR SCIENCE 2019; 1:5966. [PMID: 36204133 PMCID: PMC9534443 DOI: 10.33011/livecoms.1.1.5966] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/06/2023]
Abstract
We establish a reliable and robust standardization of settings for practical molecular dynamics (MD) simulations of pure and mixed (single- and multi-component) lipid bilayer membranes. In lipid membranes research, particle-based molecular simulations are a powerful tool alongside continuum theory, lipidomics, and model, in vitro, and in vivo experiments. Molecular simulations can provide precise and reproducible spatiotemporal (atomic- and femtosecond-level) information about membrane structure, mechanics, thermodynamics, kinetics, and dynamics. Yet the simulation of lipid membranes can be a daunting task, given the uniqueness of lipid membranes relative to conventional liquid-liquid and solid-liquid interfaces, the immense and complex thermodynamic and statistical mechanical theory, the diversity of multiscale lipid models, limitations of modern computing power, the difficulty and ambiguity of simulation controls, finite size effects, competitive continuum simulation alternatives, and the desired application, including vesicle experiments and biological membranes. These issues can complicate an essential understanding of the field of lipid membranes, and create major bottlenecks to simulation advancement. In this article, we clarify these issues and present a consistent, thorough, and user-friendly framework for the design of state-of-the-art lipid membrane MD simulations. We hope to allow early-career researchers to quickly overcome common obstacles in the field of lipid membranes and reach maximal impact in their simulations.
Collapse
Affiliation(s)
- David J. Smith
- Department of Chemical Engineering, University of California, Santa Barbara, Santa Barbara, CA, USA
| | - Jeffery B. Klauda
- Department of Chemical and Biomolecular Engineering and Biophysics Program, University of Maryland, College Park, MD, USA
| | - Alexander J. Sodt
- Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, MD, USA
| |
Collapse
|
159
|
Doktorova M, LeVine MV, Khelashvili G, Weinstein H. A New Computational Method for Membrane Compressibility: Bilayer Mechanical Thickness Revisited. Biophys J 2019; 116:487-502. [PMID: 30665693 DOI: 10.1016/j.bpj.2018.12.016] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2018] [Revised: 11/27/2018] [Accepted: 12/20/2018] [Indexed: 01/24/2023] Open
Abstract
Because lipid bilayers can bend and stretch in ways similar to thin elastic sheets, physical models of bilayer deformation have utilized mechanical constants such as the moduli for bending rigidity (κC) and area compressibility (KA). However, the use of these models to quantify the energetics of membrane deformation associated with protein-membrane interactions, and the membrane response to stress is often hampered by the shortage of experimental data suitable for the estimation of the mechanical constants of various lipid mixtures. Although computational tools such as molecular dynamics simulations can provide alternative means to estimate KA values, current approaches suffer significant technical limitations. Here, we present a novel, to our knowledge, computational framework that allows for a direct estimation of KA values for individual bilayer leaflets. The theory is based on the concept of elasticity and derives KA from real-space analysis of local thickness fluctuations sampled in molecular dynamics simulations. We explore and validate the model on a large set of single and multicomponent bilayers of different lipid compositions and sizes, simulated at different temperatures. The calculated bilayer compressibility moduli agree with values estimated previously from experiments and those obtained from a standard computational method based on a series of constrained tension simulations. We further validate our framework in a comparison with an existing polymer brush model and confirm the polymer brush model's predicted linear relationship with proportionality coefficient of 24, using elastic parameters calculated from the simulation trajectories. The robustness of the results that emerge from the method allows us to revisit the origins of the bilayer mechanical (compressible) thickness and in particular its dependence on acyl-chain unsaturation and the presence of cholesterol.
Collapse
Affiliation(s)
- Milka Doktorova
- Tri-Institutional PhD Program in Computational Biology and Medicine, Weill Cornell Medical College, New York, New York.
| | - Michael V LeVine
- Department of Physiology and Biophysics, Weill Cornell Medicine, New York, New York; The HRH Prince Alwaleed Bin Talal Bin Abdulaziz Alsaud Institute for Computational Biomedicine, Weill Greenberg Center, New York, New York
| | - George Khelashvili
- Department of Physiology and Biophysics, Weill Cornell Medicine, New York, New York; The HRH Prince Alwaleed Bin Talal Bin Abdulaziz Alsaud Institute for Computational Biomedicine, Weill Greenberg Center, New York, New York
| | - Harel Weinstein
- Department of Physiology and Biophysics, Weill Cornell Medicine, New York, New York; The HRH Prince Alwaleed Bin Talal Bin Abdulaziz Alsaud Institute for Computational Biomedicine, Weill Greenberg Center, New York, New York
| |
Collapse
|
160
|
Mechanical properties of bilayers containing sperm sphingomyelins and ceramides with very long-chain polyunsaturated fatty acids. Chem Phys Lipids 2019; 218:178-186. [PMID: 30610838 DOI: 10.1016/j.chemphyslip.2018.12.008] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2018] [Revised: 12/17/2018] [Accepted: 12/19/2018] [Indexed: 02/02/2023]
Abstract
Sphingomyelins (SM) and ceramides (Cer) with very long chain polyunsaturated fatty acids (V) are important components of spermatozoa membranes. In this study, the mechanical properties of bilayers of SM and Cer with nonhydroxy (n-V) and 2-hydroxy (h-V) fatty acid (30:5) were studied by molecular dynamics simulation at different temperatures and in the presence and the absence of salt. From our results, it was evidenced how n-V SM and h-V SM bilayers showed similar behavior. When n-V Cer was added to a h-V SM bilayer, the Gaussian curvature modulus and Ecurve of binary bilayers decreased. This variation in the mechanical properties of the bilayer can be associated with an incipient step during the fecundation process.
Collapse
|
161
|
Mallikarjunaiah KJ, Kinnun JJ, Petrache HI, Brown MF. Flexible lipid nanomaterials studied by NMR spectroscopy. Phys Chem Chem Phys 2019; 21:18422-18457. [DOI: 10.1039/c8cp06179c] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Advances in solid-state nuclear magnetic resonance spectroscopy inform the emergence of material properties from atomistic-level interactions in membrane lipid nanostructures.
Collapse
Affiliation(s)
- K. J. Mallikarjunaiah
- Department of Chemistry and Biochemistry
- University of Arizona
- Tucson
- USA
- Department of Physics
| | - Jacob J. Kinnun
- Department of Physics
- Indiana University-Purdue University
- Indianapolis
- USA
| | - Horia I. Petrache
- Department of Physics
- Indiana University-Purdue University
- Indianapolis
- USA
| | - Michael F. Brown
- Department of Chemistry and Biochemistry
- University of Arizona
- Tucson
- USA
- Department of Physics
| |
Collapse
|
162
|
Oroskar PA, Jameson CJ, Murad S. Molecular-Level "Observations" of the Behavior of Gold Nanoparticles in Aqueous Solution and Interacting with a Lipid Bilayer Membrane. Methods Mol Biol 2019; 2000:303-359. [PMID: 31148024 DOI: 10.1007/978-1-4939-9516-5_21] [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] [Indexed: 06/09/2023]
Abstract
We use coarse-grained molecular dynamics simulations to "observe" details of interactions between ligand-covered gold nanoparticles and a lipid bilayer model membrane. In molecular dynamics simulations, one puts the individual atoms and groups of atoms of the physical system to be "observed" into a simulation box, specifies the forms of the potential energies of interactions between them (ultimately quantum based), and lets them individually move classically according to Newton's equations of motion, based on the forces arising from the assumed potential energy forms. The atoms that are chemically bonded to each other stay chemically bonded, following known potentials (force fields) that permit internal degrees of freedom (internal rotation, torsion, vibrations), and the interactions between nonbonded atoms are simplified to Lennard-Jones forms (in our case) and coulombic (where electrical charges are present) in which the parameters are previously optimized to reproduce thermodynamic properties or are based on quantum electronic calculations. The system is started out at a reasonable set of coordinates for all atoms or groups of atoms, and then permitted to develop according to the equations of motion, one small step (usually 10 fs time step) at a time, for millions of steps until the system is at a quasi-equilibrium (usually reached after hundreds of nanoseconds). We then let the system play out its motions further for many nanoseconds to observe the behavior, periodically taking snapshots (saving all positions and energies), and post-processing the snapshots to obtain various average descriptions of the system. Alkanethiols of various lengths serve as examples of hydrophobic ligands and methyl-terminated PEG with various numbers of monomer units serve as examples of hydrophilic ligands. Spherical gold particles of various diameters as well as gold nanorods form the core to which ligands are attached. The nanoparticles are characterized at the molecular level, especially the distributions of ligand configurations and their dependence on ligand length, and surface coverage. Self-assembly of the bilayer from an isotropic solution and observation of membrane properties that correspond well to experimental values validate the simulations. The mechanism of permeation of a gold NP coated with either a hydrophobic or a hydrophilic ligand, and its dependence on surface coverage, ligand length, core diameter, and core shape, is investigated. Lipid response such as lipid flip-flops, lipid extraction, and changes in order parameter of the lipid tails are examined in detail. The mechanism of permeation of a PEGylated nanorod is shown to occur by tilting, lying down, rotating, and straightening up. The nature of the information provided by molecular dynamics simulations permits understanding of the detailed behavior of gold nanoparticles interacting with lipid membranes which in turn helps to understand why some known systems work better than others and aids the design of new particles and improvement of methods for preparing existing ones.
Collapse
Affiliation(s)
- Priyanka A Oroskar
- Department of Chemical Engineering, University of Illinois at Chicago, Chicago, IL, USA
| | - Cynthia J Jameson
- Department of Chemistry, University of Illinois at Chicago, Chicago, IL, USA
| | - Sohail Murad
- Department of Chemical Engineering, University of Illinois at Chicago, Chicago, IL, USA.
- Department of Chemical Engineering, Illinois Institute of Technology, Chicago, IL, USA.
| |
Collapse
|
163
|
Yu Y, Klauda JB. Modeling Pseudomonas aeruginosa inner plasma membrane in planktonic and biofilm modes. J Chem Phys 2018; 149:215102. [DOI: 10.1063/1.5052629] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yalun Yu
- Biophysics Graduate Program, University of Maryland, College Park, Maryland 20742, USA
| | - Jeffery B. Klauda
- Biophysics Graduate Program, University of Maryland, College Park, Maryland 20742, USA
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, Maryland 20742, USA
| |
Collapse
|
164
|
Distribution of mechanical stress in the Escherichia coli cell envelope. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2018; 1860:2566-2575. [PMID: 30278180 DOI: 10.1016/j.bbamem.2018.09.020] [Citation(s) in RCA: 54] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/23/2018] [Revised: 09/14/2018] [Accepted: 09/26/2018] [Indexed: 01/05/2023]
Abstract
The cell envelope in Gram-negative bacteria comprises two distinct membranes with a cell wall between them. There has been a growing interest in understanding the mechanical adaptation of this cell envelope to the osmotic pressure (or turgor pressure), which is generated by the difference in the concentration of solutes between the cytoplasm and the external environment. However, it remains unexplored how the cell wall, the inner membrane (IM), and the outer membrane (OM) effectively protect the cell from this pressure by bearing the resulting surface tension, thus preventing the formation of inner membrane bulges, abnormal cell morphology, spheroplasts and cell lysis. In this study, we have used molecular dynamics (MD) simulations combined with experiments to resolve how and to what extent models of the IM, OM, and cell wall respond to changes in surface tension. We calculated the area compressibility modulus of all three components in simulations from tension-area isotherms. Experiments on monolayers mimicking individual leaflets of the IM and OM were also used to characterize their compressibility. While the membranes become softer as they expand, the cell wall exhibits significant strain stiffening at moderate to high tensions. We integrate these results into a model of the cell envelope in which the OM and cell wall share the tension at low turgor pressure (0.3 atm) but the tension in the cell wall dominates at high values (>1 atm).
Collapse
|
165
|
Sajadi F, Rowley CN. Simulations of lipid bilayers using the CHARMM36 force field with the TIP3P-FB and TIP4P-FB water models. PeerJ 2018; 6:e5472. [PMID: 30128211 PMCID: PMC6097494 DOI: 10.7717/peerj.5472] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2018] [Accepted: 07/28/2018] [Indexed: 12/13/2022] Open
Abstract
The CHARMM36 force field for lipids is widely used in simulations of lipid bilayers. The CHARMM family of force fields were developed for use with the mTIP3P water model. This water model has an anomalously high dielectric constant and low viscosity, which limits its accuracy in the calculation of quantities like permeability coefficients. The TIP3P-FB and TIP4P-FB water models are more accurate in terms of the dielectric constant and transport properties, which could allow more accurate simulations of systems containing water and lipids. To test whether the CHARMM36 lipid force field is compatible with the TIP3P-FB and TIP4P-FB water models, we have performed simulations of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayers. The calculated headgroup area, compressibility, order parameters, and X-ray form factors are in good agreement with the experimental values, indicating that these improved water models can be used with the CHARMM36 lipid force field without modification when calculating membrane physical properties. The water permeability predicted by these models is significantly different; the mTIP3P-model diffusion in solution and at the lipid-water interface is anomalously fast due to the spuriously low viscosity of mTIP3P-model water, but the potential of mean force of permeation is higher for the TIP3P-FB and TIP4P-FB models due to their high excess chemical potentials. As a result, the rates of water permeation calculated the FB water models are slower than the experimental value by a factor of 15-17, while simulations with the mTIP3P model only underestimate the water permeability by a factor of 3.
Collapse
Affiliation(s)
- Fatima Sajadi
- Department of Chemistry, Memorial University of Newfoundland, St. John's, NL, Canada
| | - Christopher N Rowley
- Department of Chemistry, Memorial University of Newfoundland, St. John's, NL, Canada
| |
Collapse
|
166
|
Lessen HJ, Fleming PJ, Fleming KG, Sodt AJ. Building Blocks of the Outer Membrane: Calculating a General Elastic Energy Model for β-Barrel Membrane Proteins. J Chem Theory Comput 2018; 14:4487-4497. [PMID: 29979594 PMCID: PMC6191857 DOI: 10.1021/acs.jctc.8b00377] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The outer membranes of Gram negative bacteria are the first points of contact these organisms make with their environment. Understanding how composition determines the mechanical properties of this essential barrier is of paramount importance. Therefore, we developed a new computational method to measure the elasticity of transmembrane proteins found in the outer membrane. Using all-atom molecular dynamics simulations of these proteins, we apply a set of external forces to mechanically stress the transmembrane β-barrels. Our results from four representative β-barrels show that outer membrane proteins display elastic properties that are approximately 70 to 190 times stiffer than neat lipid membranes. These findings suggest that outer membrane β-barrels are a significant source of mechanical stability in bacteria. Our all-atom approach further reveals that resistance to radial stress is encoded by a general mechanism that includes stretching of backbone hydrogen bonds and tilting of β-strands with respect to the bilayer normal. This computational framework facilitates an increased theoretical understanding of how varying lipid and protein amounts affect the mechanical properties of the bacterial outer membrane.
Collapse
Affiliation(s)
- Henry J. Lessen
- T. C. Jenkins Department of Biophysics, Johns Hopkins University, Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health
| | - Patrick J. Fleming
- T. C. Jenkins Department of Biophysics, Johns Hopkins University, Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health
| | - Karen G. Fleming
- T. C. Jenkins Department of Biophysics, Johns Hopkins University, Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health
| | - Alexander J. Sodt
- T. C. Jenkins Department of Biophysics, Johns Hopkins University, Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health
| |
Collapse
|
167
|
HTMoL: full-stack solution for remote access, visualization, and analysis of molecular dynamics trajectory data. J Comput Aided Mol Des 2018; 32:869-876. [DOI: 10.1007/s10822-018-0141-y] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2018] [Accepted: 08/01/2018] [Indexed: 10/28/2022]
|
168
|
Bassereau P, Jin R, Baumgart T, Deserno M, Dimova R, Frolov VA, Bashkirov PV, Grubmüller H, Jahn R, Risselada HJ, Johannes L, Kozlov MM, Lipowsky R, Pucadyil TJ, Zeno WF, Stachowiak JC, Stamou D, Breuer A, Lauritsen L, Simon C, Sykes C, Voth GA, Weikl TR. The 2018 biomembrane curvature and remodeling roadmap. JOURNAL OF PHYSICS D: APPLIED PHYSICS 2018; 51:343001. [PMID: 30655651 PMCID: PMC6333427 DOI: 10.1088/1361-6463/aacb98] [Citation(s) in RCA: 161] [Impact Index Per Article: 26.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/14/2023]
Abstract
The importance of curvature as a structural feature of biological membranes has been recognized for many years and has fascinated scientists from a wide range of different backgrounds. On the one hand, changes in membrane morphology are involved in a plethora of phenomena involving the plasma membrane of eukaryotic cells, including endo- and exocytosis, phagocytosis and filopodia formation. On the other hand, a multitude of intracellular processes at the level of organelles rely on generation, modulation, and maintenance of membrane curvature to maintain the organelle shape and functionality. The contribution of biophysicists and biologists is essential for shedding light on the mechanistic understanding and quantification of these processes. Given the vast complexity of phenomena and mechanisms involved in the coupling between membrane shape and function, it is not always clear in what direction to advance to eventually arrive at an exhaustive understanding of this important research area. The 2018 Biomembrane Curvature and Remodeling Roadmap of Journal of Physics D: Applied Physics addresses this need for clarity and is intended to provide guidance both for students who have just entered the field as well as established scientists who would like to improve their orientation within this fascinating area.
Collapse
Affiliation(s)
- Patricia Bassereau
- Laboratoire Physico Chimie Curie, Institut Curie, PSL Research University, CNRS UMR168, 75005 Paris, France
- Sorbonne Université, 75005 Paris, France
| | - Rui Jin
- Chemistry Department, University of Pennsylvania, Philadelphia, PA 19104-6323, United States of America
| | - Tobias Baumgart
- Chemistry Department, University of Pennsylvania, Philadelphia, PA 19104-6323, United States of America
| | - Markus Deserno
- Department of Physics, Carnegie Mellon University, Pittsburgh, PA 15213, United States of America
| | - Rumiana Dimova
- Department of Theory and Bio-Systems, Max Planck Institute of Colloids and Interfaces, Science Park Golm, 14424 Potsdam, Germany
| | - Vadim A Frolov
- Biofisika Institute (CSIC, UPV/EHU) and Department of Biochemistry and Molecular Biology, University of the Basque Country, Leioa 48940, Spain
- IKERBASQUE, Basque Foundation for Science, Bilbao 48013, Spain
| | - Pavel V Bashkirov
- Federal Research and Clinical Centre of Physical-Chemical Medicine, Moscow 119435, Russia
- A.N. Frumkin Institute of Physical Chemistry and Electrochemistry, Russian Academy of Sciences, Moscow 119071, Russia
| | - Helmut Grubmüller
- Department of Theoretical and Computational Biophysics, Max Planck Institute for Biophysical Chemistry, Göttingen, Germany
| | - Reinhard Jahn
- Department of Neurobiology, Max Planck Institute for Biophysical Chemistry, Göttingen, Germany
| | - H Jelger Risselada
- Department of Theoretical Physics, Georg-August University, Göttingen, Germany
| | - Ludger Johannes
- Cellular and Chemical Biology Unit, Institut Curie, PSL Research University, U1143 INSERM, UMR3666 CNRS, 26 rue d'Ulm, 75248 Paris Cedex 05, France
| | - Michael M Kozlov
- Sackler Faculty of Medicine, Department of Physiology and Pharmacology, Tel Aviv University
| | - Reinhard Lipowsky
- Department of Theory and Bio-Systems, Max Planck Institute of Colloids and Interfaces, Science Park Golm, 14424 Potsdam, Germany
| | | | - Wade F Zeno
- Department of Biomedical Engineering, University of Texas at Austin, Austin, TX, United States of America
| | - Jeanne C Stachowiak
- Department of Biomedical Engineering, University of Texas at Austin, Austin, TX, United States of America
- University of Texas at Austin, Institute for Cellular and Molecular Biology, Austin, TX, United States of America
| | - Dimitrios Stamou
- Bionanotechnology and Nanomedicine Laboratory, Department of Chemistry, Nano-Science Center, University of Copenhagen, Denmark
| | - Artú Breuer
- Bionanotechnology and Nanomedicine Laboratory, Department of Chemistry, Nano-Science Center, University of Copenhagen, Denmark
| | - Line Lauritsen
- Bionanotechnology and Nanomedicine Laboratory, Department of Chemistry, Nano-Science Center, University of Copenhagen, Denmark
| | - Camille Simon
- Laboratoire Physico Chimie Curie, Institut Curie, PSL Research University, CNRS UMR168, 75005 Paris, France
- Sorbonne Université, 75005 Paris, France
| | - Cécile Sykes
- Laboratoire Physico Chimie Curie, Institut Curie, PSL Research University, CNRS UMR168, 75005 Paris, France
- Sorbonne Université, 75005 Paris, France
| | - Gregory A Voth
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, IL, United States of America
| | - Thomas R Weikl
- Department of Theory and Bio-Systems, Max Planck Institute of Colloids and Interfaces, Science Park Golm, 14424 Potsdam, Germany
| |
Collapse
|
169
|
König G, Reetz MT, Thiel W. 1-Butanol as a Solvent for Efficient Extraction of Polar Compounds from Aqueous Medium: Theoretical and Practical Aspects. J Phys Chem B 2018; 122:6975-6988. [PMID: 29897756 DOI: 10.1021/acs.jpcb.8b02877] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The extraction of polar molecules from aqueous solution is a challenging task in organic synthesis. 1-Butanol has been used sporadically as an eluent for polar molecules, but it is unclear which molecular features drive its efficiency. Here, we employ free energy simulations to study the partitioning of 15 solutes between water and 1-butanol. The simulations demonstrate that the high affinity of polar molecules to the wet 1-butanol phase is associated with its nanostructure. Small inverse micelles of water are able to accommodate polar solutes and locally mimic an aqueous environment. We verify the simulations based on partition coefficients between water and 1-octanol, and include a blind prediction of the water/1-butanol partition coefficient of cyclohexane-1,2-diol. The calculations are in excellent agreement with experiment, reaching root-mean-square deviations below 0.7 kcal/mol. Actual extractions of cyclohexane-1,2-diol from buffer solutions that mimic cell lysates and suspensions in biocatalytic reactions further exemplify our findings. The yields highlight that extractions with 1-butanol can be significantly more efficient than the conventional protocol based on ethyl acetate.
Collapse
Affiliation(s)
- Gerhard König
- Max-Planck-Institut für Kohlenforschung , 45470 Mülheim an der Ruhr , Germany.,Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research, and Department of Chemistry and Chemical Biology , Rutgers University , Piscataway , New Jersey 08854 , United States
| | - Manfred T Reetz
- Max-Planck-Institut für Kohlenforschung , 45470 Mülheim an der Ruhr , Germany.,Department of Chemistry , Philipps-University Marburg , 35032 Marburg , Germany
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung , 45470 Mülheim an der Ruhr , Germany
| |
Collapse
|
170
|
Leonard AN, Pastor RW, Klauda JB. Parameterization of the CHARMM All-Atom Force Field for Ether Lipids and Model Linear Ethers. J Phys Chem B 2018; 122:6744-6754. [PMID: 29870257 DOI: 10.1021/acs.jpcb.8b02743] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
Linear ethers such as polyethylene glycol have extensive industrial and medical applications. Additionally, phospholipids containing an ether linkage between the glycerol backbone and hydrophobic tails are prevalent in human red blood cells and nerve tissue. This study uses ab initio results to revise the CHARMM additive (C36) partial-charge and dihedral parameters for linear ethers and develop parameters for the ether-linked phospholipid 1,2-di- O-hexadecyl- sn-glycero-3-phosphocholine (DHPC). The new force field, called C36e, more accurately represents the dihedral potential energy landscape and improves the densities and free energies of hydration of linear ethers. C36e allows more water to penetrate into a DHPC bilayer, increasing the surface area per lipid compared to simulations carried out with the original C36 ether parameters and improving the overall structural properties obtained from X-ray and neutron scattering. Comparison with an ester-linked DPPC bilayer (1,2-dipalmitoyl- sn-phosphatidylcholine) reveals that the ether linkage increases water organization in the headgroup region. This effect is a likely explanation for the experimentally lower water permeability of bilayers composed of ether-linked lipids.
Collapse
Affiliation(s)
| | - Richard W Pastor
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute , National Institutes of Health , Bethesda , Maryland 20892 , United States
| | | |
Collapse
|
171
|
Abstract
On mesoscopic scales, lipid membranes are well described by continuum theories whose main ingredients are the curvature of a membrane's reference surface and the tilt of its lipid constituents. In particular, Hamm and Kozlov [Eur. Phys. J. E 3, 323 (2000)] have shown how to systematically derive such a tilt-curvature Hamiltonian based on the elementary assumption of a thin fluid elastic sheet experiencing internal lateral pre-stress. Performing a dimensional reduction, they not only derive the basic form of the effective surface Hamiltonian but also express its emergent elastic couplings as trans-membrane moments of lower-level material parameters. In the present paper, we argue, though, that their derivation unfortunately missed a coupling term between curvature and tilt. This term arises because, as one moves along the membrane, the curvature-induced change of transverse distances contributes to the area strain-an effect that was believed to be small but nevertheless ends up contributing at the same (quadratic) order as all other terms in their Hamiltonian. We illustrate the consequences of this amendment by deriving the monolayer and bilayer Euler-Lagrange equations for the tilt, as well as the power spectra of shape, tilt, and director fluctuations. A particularly curious aspect of our new term is that its associated coupling constant is the second moment of the lipid monolayer's lateral stress profile-which within this framework is equal to the monolayer Gaussian curvature modulus, κ¯m. On the one hand, this implies that many theoretical predictions now contain a parameter that is poorly known (because the Gauss-Bonnet theorem limits access to the integrated Gaussian curvature); on the other hand, the appearance of κ¯m outside of its Gaussian curvature provenance opens opportunities for measuring it by more conventional means, for instance by monitoring a membrane's undulation spectrum at short scales.
Collapse
Affiliation(s)
- M Mert Terzi
- Department of Physics, Carnegie Mellon University, 5000 Forbes Ave., Pittsburgh, Pennsylvania 15213, USA
| | - Markus Deserno
- Department of Physics, Carnegie Mellon University, 5000 Forbes Ave., Pittsburgh, Pennsylvania 15213, USA
| |
Collapse
|
172
|
Guiselin B, Law JO, Chakrabarti B, Kusumaatmaja H. Dynamic Morphologies and Stability of Droplet Interface Bilayers. PHYSICAL REVIEW LETTERS 2018; 120:238001. [PMID: 29932701 DOI: 10.1103/physrevlett.120.238001] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/24/2017] [Revised: 01/24/2018] [Indexed: 06/08/2023]
Abstract
We develop a theoretical framework for understanding dynamic morphologies and stability of droplet interface bilayers (DIBs), accounting for lipid kinetics in the monolayers and bilayer, and droplet evaporation due to imbalance between osmotic and Laplace pressures. Our theory quantitatively describes distinct pathways observed in experiments when DIBs become unstable. We find that when the timescale for lipid desorption is slow compared to droplet evaporation, the lipid bilayer will grow and the droplets approach a hemispherical shape. In contrast, when lipid desorption is fast, the bilayer area will shrink and the droplets eventually detach. Our model also suggests there is a critical size below which DIBs can become unstable, which may explain experimental difficulties in miniaturizing the DIB platform.
Collapse
Affiliation(s)
- Benjamin Guiselin
- Department of Physics, Durham University, Durham DH1 3LE, United Kingdom
| | - Jack O Law
- Department of Physics, Durham University, Durham DH1 3LE, United Kingdom
| | - Buddhapriya Chakrabarti
- Department of Physics and Astronomy, University of Sheffield, Sheffield S3 7RH, United Kingdom
| | - Halim Kusumaatmaja
- Department of Physics, Durham University, Durham DH1 3LE, United Kingdom
| |
Collapse
|
173
|
Grillo DA, Albano JMR, Mocskos EE, Facelli JC, Pickholz M, Ferraro MB. Mechanical properties of drug loaded diblock copolymer bilayers: A molecular dynamics study. J Chem Phys 2018; 148:214901. [PMID: 29884038 DOI: 10.1063/1.5028377] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, we present results of coarse-grained simulations to study the encapsulation of prilocaine (PLC), both neutral and protonated, on copolymer bilayers through molecular dynamics simulations. Using a previously validated membrane model, we have simulated loaded bilayers at different drug concentrations and at low (protonated PLC) and high (neutral PLC) pH levels. We have characterized key structural parameters of the loaded bilayers in order to understand the effects of encapsulation of PLC on the bilayer structure and mechanical properties. Neutral PLC was encapsulated in the hydrophobic region leading to a thickness increase, while the protonated species partitioned between the water phase and the poly(ethylene oxide)-poly(butadiene) (PBD) interface, relaxing the PBD region and leading to a decrease in the thickness. The tangential pressures of the studied systems were calculated, and their components were decomposed in order to gain insights on their compensation. In all cases, it is observed that the loading of the membrane does not significantly decrease the stability of the bilayer, indicating that the system could be used for drug delivery.
Collapse
Affiliation(s)
- Damián A Grillo
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, Argentina
| | - Juan M R Albano
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, Argentina
| | - Esteban E Mocskos
- Departamento de Computación, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, Argentina
| | - Julio C Facelli
- Department of Biomedical Informatics, University of Utah, 421 Wakara Way, Suite 140, Salt Lake City, Utah 84108, USA
| | - Mónica Pickholz
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, Argentina
| | - Marta B Ferraro
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, Argentina
| |
Collapse
|
174
|
Doskocz J, Drabik D, Chodaczek G, Przybyło M, Langner M. Statistical Analysis of Bending Rigidity Coefficient Determined Using Fluorescence-Based Flicker-Noise Spectroscopy. J Membr Biol 2018; 251:601-608. [PMID: 29858612 DOI: 10.1007/s00232-018-0037-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2017] [Accepted: 05/29/2018] [Indexed: 10/14/2022]
Abstract
Bending rigidity coefficient describes propensity of a lipid bilayer to deform. In order to measure the parameter experimentally using flickering noise spectroscopy, the microscopic imaging is required, which necessitates the application of giant unilamellar vesicles (GUV) lipid bilayer model. The major difficulty associated with the application of the model is the statistical character of GUV population with respect to their size and the homogeneity of lipid bilayer composition, if a mixture of lipids is used. In the paper, the bending rigidity coefficient was measured using the fluorescence-enhanced flicker-noise spectroscopy. In the paper, the bending rigidity coefficient was determined for large populations of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine and 1,2-dioleoyl-sn-glycero-3-phosphocholine vesicles. The quantity of obtained experimental data allows to perform statistical analysis aiming at the identification of the distribution, which is the most appropriate for the calculation of the value of the membrane bending rigidity coefficient. It has been demonstrated that the bending rigidity coefficient is characterized by an asymmetrical distribution, which is well approximated with the gamma distribution. Since there are no biophysical reasons for that we propose to use the difference between normal and gamma fits as a measure of the homogeneity of vesicle population. In addition, the effect of a fluorescent label and types of instrumental setups on determined values has been tested. Obtained results show that the value of the bending rigidity coefficient does not depend on the type of a fluorescent label nor on the type of microscope used.
Collapse
Affiliation(s)
- Joanna Doskocz
- Department of Biomedical Engineering, Faculty of Fundamental Problems of Technology, Wrocław University of Science and Technology, Pl. Grunwaldzki 13, 50-377, Wrocław, Poland.
| | - Dominik Drabik
- Department of Biomedical Engineering, Faculty of Fundamental Problems of Technology, Wrocław University of Science and Technology, Pl. Grunwaldzki 13, 50-377, Wrocław, Poland
| | - Grzegorz Chodaczek
- Wroclaw Research Centre EIT+, ul. Stabłowicka 147, 54-066, Wrocław, Poland
| | - Magdalena Przybyło
- Department of Biomedical Engineering, Faculty of Fundamental Problems of Technology, Wrocław University of Science and Technology, Pl. Grunwaldzki 13, 50-377, Wrocław, Poland.,Lipid Systems sp. z o.o., ul. Krzemieniecka 48C, 54-613, Wrocław, Poland
| | - Marek Langner
- Department of Biomedical Engineering, Faculty of Fundamental Problems of Technology, Wrocław University of Science and Technology, Pl. Grunwaldzki 13, 50-377, Wrocław, Poland.,Lipid Systems sp. z o.o., ul. Krzemieniecka 48C, 54-613, Wrocław, Poland
| |
Collapse
|
175
|
Moths B, Witten TA. Engineering single-polymer micelle shape using nonuniform spontaneous surface curvature. Phys Rev E 2018; 97:032505. [PMID: 29776184 DOI: 10.1103/physreve.97.032505] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2017] [Indexed: 11/07/2022]
Abstract
Conventional micelles, composed of simple amphiphiles, exhibit only a few standard morphologies, each characterized by its mean surface curvature set by the amphiphiles. Here we demonstrate a rational design scheme to construct micelles of more general shape from polymeric amphiphiles. We replace the many amphiphiles of a conventional micelle by a single flexible, linear, block copolymer chain containing two incompatible species arranged in multiple alternating segments. With suitable segment lengths, the chain exhibits a condensed spherical configuration in solution, similar to conventional micelles. Our design scheme posits that further shapes are attained by altering the segment lengths. As a first study of the power of this scheme, we demonstrate the capacity to produce long-lived micelles of horseshoe form using conventional bead-spring simulations in two dimensions. Modest changes in the segment lengths produce smooth changes in the micelle's shape and stability.
Collapse
Affiliation(s)
- Brian Moths
- Department of Physics and James Franck Institute, University of Chicago, Chicago, Illinois 60637, USA
| | - T A Witten
- Department of Physics and James Franck Institute, University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
176
|
Molecular dynamics simulations of lipid nanodiscs. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2018; 1860:2094-2107. [PMID: 29729280 DOI: 10.1016/j.bbamem.2018.04.015] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/31/2017] [Revised: 04/27/2018] [Accepted: 04/28/2018] [Indexed: 01/02/2023]
Abstract
A lipid nanodisc is a discoidal lipid bilayer stabilized by proteins, peptides, or polymers on its edge. Nanodiscs have two important connections to structural biology. The first is associated with high-density lipoprotein (HDL), a particle with a variety of functionalities including lipid transport. Nascent HDL (nHDL) is a nanodisc stabilized by Apolipoprotein A-I (APOA1). Determining the structure of APOA1 and its mimetic peptides in nanodiscs is crucial to understanding pathologies related to HDL maturation and designing effective therapies. Secondly, nanodiscs offer non-detergent membrane-mimicking environments and greatly facilitate structural studies of membrane proteins. Although seemingly similar, natural and synthetic nanodiscs are different in that nHDL is heterogeneous in size, due to APOA1 elasticity, and gradually matures to become spherical. Synthetic nanodiscs, in contrast, should be homogenous, stable, and size-tunable. This report reviews previous molecular dynamics (MD) simulation studies of nanodiscs and illustrates convergence and accuracy issues using results from new multi-microsecond atomistic MD simulations. These new simulations reveal that APOA1 helices take 10-20 μs to rearrange on the nanodisc, while peptides take 2 μs to migrate from the disc surfaces to the edge. These systems can also become kinetically trapped depending on the initial conditions. For example, APOA1 was trapped in a biologically irrelevant conformation for the duration of a 10 μs trajectory; the peptides were similarly trapped for 5 μs. It therefore remains essential to validate MD simulations of these systems with experiments due to convergence and accuracy issues. This article is part of a Special Issue entitled: Emergence of Complex Behavior in Biomembranes edited by Marjorie Longo.
Collapse
|
177
|
Padariya M, Kalathiya U, Baginski M. Structural and dynamic insights on the EmrE protein with TPP + and related substrates through molecular dynamics simulations. Chem Phys Lipids 2018; 212:1-11. [DOI: 10.1016/j.chemphyslip.2017.12.004] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2017] [Revised: 11/30/2017] [Accepted: 12/22/2017] [Indexed: 12/21/2022]
|
178
|
Pezeshkian W, Khandelia H, Marsh D. Lipid Configurations from Molecular Dynamics Simulations. Biophys J 2018; 114:1895-1907. [PMID: 29694867 PMCID: PMC5937052 DOI: 10.1016/j.bpj.2018.02.016] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2017] [Revised: 02/09/2018] [Accepted: 02/13/2018] [Indexed: 01/10/2023] Open
Abstract
The extent to which current force fields faithfully reproduce conformational properties of lipids in bilayer membranes, and whether these reflect the structural principles established for phospholipids in bilayer crystals, are central to biomembrane simulations. We determine the distribution of dihedral angles in palmitoyl-oleoyl phosphatidylcholine from molecular dynamics simulations of hydrated fluid bilayer membranes. We compare results from the widely used lipid force field of Berger et al. with those from the most recent C36 release of the CHARMM force field for lipids. Only the CHARMM force field produces the chain inequivalence with sn-1 as leading chain that is characteristic of glycerolipid packing in fluid bilayers. The exposure and high partial charge of the backbone carbonyls in Berger lipids leads to artifactual binding of Na+ ions reported in the literature. Both force fields predict coupled, near-symmetrical distributions of headgroup dihedral angles, which is compatible with models of interconverting mirror-image conformations used originally to interpret NMR order parameters. The Berger force field produces rotamer populations that correspond to the headgroup conformation found in a phosphatidylcholine lipid bilayer crystal, whereas CHARMM36 rotamer populations are closer to the more relaxed crystal conformations of phosphatidylethanolamine and glycerophosphocholine. CHARMM36 alone predicts the correct relative signs of the time-average headgroup order parameters, and reasonably reproduces the full range of NMR data from the phosphate diester to the choline methyls. There is strong motivation to seek further experimental criteria for verifying predicted conformational distributions in the choline headgroup, including the 31P chemical shift anisotropy and 14N and CD3 NMR quadrupole splittings.
Collapse
Affiliation(s)
- Weria Pezeshkian
- MEMPHYS-Centre for Biomembrane Physics, University of Southern Denmark, Odense M, Denmark
| | - Himanshu Khandelia
- MEMPHYS-Centre for Biomembrane Physics, University of Southern Denmark, Odense M, Denmark
| | - Derek Marsh
- MEMPHYS-Centre for Biomembrane Physics, University of Southern Denmark, Odense M, Denmark; Max-Planck-Institut für biophysikalische Chemie, Göttingen, Germany.
| |
Collapse
|
179
|
Agmon E, Solon J, Bassereau P, Stockwell BR. Modeling the effects of lipid peroxidation during ferroptosis on membrane properties. Sci Rep 2018; 8:5155. [PMID: 29581451 PMCID: PMC5979948 DOI: 10.1038/s41598-018-23408-0] [Citation(s) in RCA: 203] [Impact Index Per Article: 33.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2017] [Accepted: 03/12/2018] [Indexed: 01/20/2023] Open
Abstract
Ferroptosis is a form of regulated cell death characterized by the accumulation of lipid hydroperoxides. There has been significant research on the pathways leading to the accumulation of oxidized lipids, but the downstream effects and how lipid peroxides cause cell death during ferroptosis remain a major puzzle. We evaluated key features of ferroptosis in newly developed molecular dynamics models of lipid membranes to investigate the biophysical consequences of lipid peroxidation, and generated hypotheses about how lipid peroxides contribute to cell death during ferroptosis.
Collapse
Affiliation(s)
- Eran Agmon
- Department of Biological Sciences, MC4846, Columbia University, 550 West 120th Street, New York, NY, 10027, USA
| | - Jérôme Solon
- Laboratoire Physico Chimie Curie, Institut Curie, PSL Research University, CNRS UMR168, 75005, Paris, France.,Sorbonne Université, 75005, Paris, France.,Center for Genomic Regulation, The Barcelona Institute of Science and Technology, Barcelona, Spain.,Universitat Pompeu Fabra, Barcelona, Spain
| | - Patricia Bassereau
- Laboratoire Physico Chimie Curie, Institut Curie, PSL Research University, CNRS UMR168, 75005, Paris, France.,Sorbonne Université, 75005, Paris, France
| | - Brent R Stockwell
- Department of Biological Sciences, MC4846, Columbia University, 550 West 120th Street, New York, NY, 10027, USA. .,Department of Chemistry, MC4846, Columbia University, 550 West 120th Street, New York, NY, 10027, USA.
| |
Collapse
|
180
|
Radanović T, Reinhard J, Ballweg S, Pesek K, Ernst R. An Emerging Group of Membrane Property Sensors Controls the Physical State of Organellar Membranes to Maintain Their Identity. Bioessays 2018; 40:e1700250. [PMID: 29574931 DOI: 10.1002/bies.201700250] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2017] [Revised: 02/06/2018] [Indexed: 11/06/2022]
Abstract
The biological membranes of eukaryotic cells harbor sensitive surveillance systems to establish, sense, and maintain characteristic physicochemical properties that ultimately define organelle identity. They are fundamentally important for membrane homeostasis and play active roles in cellular signaling, protein sorting, and the formation of vesicular carriers. Here, we compare the molecular mechanisms of Mga2 and Ire1, two sensors involved in the regulation of fatty acid desaturation and the response to unfolded proteins and lipid bilayer stress in order to identify their commonalities and specializations. We will speculate on the cellular significance of membrane property sensors in other organelles and discuss their putative mechanisms. Based on these findings, we propose membrane property sensors as an emerging class of proteins with wide implications for organelle communication and function.
Collapse
Affiliation(s)
- Toni Radanović
- Medical Faculty, Department of Medical Biochemistry and Molecular Bioloy, Saarland University, 66421 Homburg, Germany
| | - John Reinhard
- Medical Faculty, Department of Medical Biochemistry and Molecular Bioloy, Saarland University, 66421 Homburg, Germany
| | - Stephanie Ballweg
- Medical Faculty, Department of Medical Biochemistry and Molecular Bioloy, Saarland University, 66421 Homburg, Germany
| | - Kristina Pesek
- Medical Faculty, Department of Medical Biochemistry and Molecular Bioloy, Saarland University, 66421 Homburg, Germany
| | - Robert Ernst
- Medical Faculty, Department of Medical Biochemistry and Molecular Bioloy, Saarland University, 66421 Homburg, Germany
| |
Collapse
|
181
|
Masone D, Uhart M, Bustos DM. Bending Lipid Bilayers: A Closed-Form Collective Variable for Effective Free-Energy Landscapes in Quantitative Biology. J Chem Theory Comput 2018; 14:2240-2245. [PMID: 29506389 DOI: 10.1021/acs.jctc.8b00012] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
Curvature-related processes are of major importance during protein-membrane interactions. The illusive simplicity of membrane reshaping masks a complex molecular process crucial for a wide range of biological functions like fusion, endo- and exocytosis, cell division, cytokinesis, and autophagy. To date, no functional expression of a reaction coordinate capable of biasing molecular dynamics simulations to produce membrane curvature has been reported. This represents a major drawback given that the adequate identification of proper collective variables to enhance sampling is fundamental for restrained dynamics techniques. In this work, we present a closed-form equation of a collective variable that induces bending in lipid bilayers in a controlled manner, allowing for straightforward calculation of free energy landscapes of important curvature-related events, using standard methods such as umbrella sampling and metadynamics. As a direct application of the collective variable, we calculate the bending free energies of a ternary lipid bilayer in the presence and the absence of a Bin/Amphiphysin/Rvs domain with an N-terminal amphipathic helix (N-BAR), a well-known peripheral membrane protein that induces curvature.
Collapse
|
182
|
Lyu Y, Xiang N, Mondal J, Zhu X, Narsimhan G. Characterization of Interactions between Curcumin and Different Types of Lipid Bilayers by Molecular Dynamics Simulation. J Phys Chem B 2018; 122:2341-2354. [PMID: 29394060 DOI: 10.1021/acs.jpcb.7b10566] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Affiliation(s)
- Yuan Lyu
- Department of Agricultural
and Biological Engineering, Purdue University, West Lafayette, Indiana 47907, United States
| | - Ning Xiang
- Department of Agricultural
and Biological Engineering, Purdue University, West Lafayette, Indiana 47907, United States
| | - Jagannath Mondal
- Centre for Interdisciplinary Sciences, Tata Institute of Fundamental Research, 36/P, Gopanapally Village, Serilingampally Mandal, Ranga Reddy District, Hyderabad 500107, India
| | - Xiao Zhu
- Research
Computing, Rosen Center for Advanced Computing, Purdue University, West Lafayette, Indiana 47907, United States
| | - Ganesan Narsimhan
- Department of Agricultural
and Biological Engineering, Purdue University, West Lafayette, Indiana 47907, United States
| |
Collapse
|
183
|
Doktorova M, Harries D, Khelashvili G. Determination of bending rigidity and tilt modulus of lipid membranes from real-space fluctuation analysis of molecular dynamics simulations. Phys Chem Chem Phys 2018. [PMID: 28627570 DOI: 10.1039/c7cp01921a] [Citation(s) in RCA: 81] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
We have recently developed a novel computational methodology (termed RSF for Real-Space Fluctuations) to quantify the bending rigidity and tilt modulus of lipid membranes from real-space analysis of fluctuations in the tilt and splay degrees of freedom as sampled in molecular dynamics (MD) simulations. In this article, we present a comprehensive study that combines results from the application of the RSF method to a wide range of lipid bilayer systems that encompass membranes of different fluidities and sizes, including lipids with saturated and unsaturated lipid tails, single and multi-component lipid systems, as well as non-standard lipids such as the four-tailed cardiolipin. By comparing the material properties calculated with the RSF method to those obtained from experimental data and from other computational methodologies, we rigorously demonstrate the validity of our approach and show its robustness. This should allow for future applications of even more complex lipidic assemblies, whose material properties are not tractable by other computational techniques. In addition, we discuss the relationship between different definitions of the tilt modulus appearing in current literature to address some important unresolved discrepancies in the field.
Collapse
Affiliation(s)
- M Doktorova
- Tri-Institutional PhD Program in Computational Biology and Medicine, Weill Cornell Medical College, 1300 York Avenue, New York, NY 10065, USA
| | | | | |
Collapse
|
184
|
Salas-Estrada LA, Leioatts N, Romo TD, Grossfield A. Lipids Alter Rhodopsin Function via Ligand-like and Solvent-like Interactions. Biophys J 2018; 114:355-367. [PMID: 29401433 PMCID: PMC5984976 DOI: 10.1016/j.bpj.2017.11.021] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2017] [Revised: 09/20/2017] [Accepted: 11/09/2017] [Indexed: 10/18/2022] Open
Abstract
Rhodopsin, a prototypical G protein-coupled receptor, is a membrane protein that can sense dim light. This highly effective photoreceptor is known to be sensitive to the composition of its lipidic environment, but the molecular mechanisms underlying this fine-tuned modulation of the receptor's function and structural stability are not fully understood. There are two competing hypotheses to explain how this occurs: 1) lipid modulation occurs via solvent-like interactions, where lipid composition controls membrane properties like hydrophobic thickness, which in turn modulate the protein's conformational equilibrium; or 2) protein-lipid interactions are ligand-like, with specific hot spots and long-lived binding events. By analyzing an ensemble of all-atom molecular dynamics simulations of five different states of rhodopsin, we show that a local ordering effect takes place in the membrane upon receptor activation. Likewise, docosahexaenoic acid acyl tails and phosphatidylethanolamine headgroups behave like weak ligands, preferentially binding to the receptor in inactive-like conformations and inducing subtle but significant structural changes.
Collapse
Affiliation(s)
- Leslie A Salas-Estrada
- Department of Biochemistry and Biophysics, University of Rochester Medical Center, Rochester, New York
| | - Nicholas Leioatts
- Department of Theoretical and Computational Biophysics, Max Plank Institute for Biophysical Chemistry, Göttingen, Germany
| | - Tod D Romo
- Center for Integrated Research Computing, University of Rochester, Rochester, New York
| | - Alan Grossfield
- Department of Biochemistry and Biophysics, University of Rochester Medical Center, Rochester, New York.
| |
Collapse
|
185
|
Leonard AN, Simmonett AC, Pickard FC, Huang J, Venable RM, Klauda JB, Brooks BR, Pastor RW. Comparison of Additive and Polarizable Models with Explicit Treatment of Long-Range Lennard-Jones Interactions Using Alkane Simulations. J Chem Theory Comput 2018; 14:948-958. [PMID: 29268012 DOI: 10.1021/acs.jctc.7b00948] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Long-range Lennard-Jones (LJ) interactions have a significant impact on the structural and thermodynamic properties of nonpolar systems. While several methods have been introduced for the treatment of long-range LJ interactions in molecular dynamics (MD) simulations, increased accuracy and extended applicability is required for anisotropic systems such as lipid bilayers. The recently refined Lennard-Jones particle-mesh Ewald (LJ-PME) method extends the particle-mesh Ewald (PME) method to long-range LJ interactions and is suitable for use with anisotropic systems. Implementation of LJ-PME with the CHARMM36 (C36) additive and CHARMM Drude polarizable force fields improves agreement with experiment for density, isothermal compressibility, surface tension, viscosity, translational diffusion, and 13C T1 relaxation times of pure alkanes. Trends in the temperature dependence of the density and isothermal compressibility of hexadecane are also improved. While the C36 additive force field with LJ-PME remains a useful model for liquid alkanes, the Drude polarizable force field with LJ-PME is more accurate for nearly all quantities considered. LJ-PME is also preferable to the isotropic long-range correction for hexadecane because the molecular order extends to nearly 20 Å, well beyond the usual 10-12 Å cutoffs used in most simulations.
Collapse
Affiliation(s)
- Alison N Leonard
- Biophysics Program, University of Maryland , College Park, Maryland 20742, United States
| | - Andrew C Simmonett
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health , Bethesda, Maryland 20892, United States
| | - Frank C Pickard
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health , Bethesda, Maryland 20892, United States
| | - Jing Huang
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health , Bethesda, Maryland 20892, United States.,Department of Pharmaceutical Science, School of Pharmacy, University of Maryland , 20 Penn Street, Baltimore, Maryland 21201, United States
| | - Richard M Venable
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health , Bethesda, Maryland 20892, United States
| | - Jeffery B Klauda
- Biophysics Program, University of Maryland , College Park, Maryland 20742, United States.,Department of Chemical and Biomolecular Engineering, University of Maryland , College Park, Maryland 20742, United States
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health , Bethesda, Maryland 20892, United States
| | - Richard W Pastor
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health , Bethesda, Maryland 20892, United States
| |
Collapse
|
186
|
Dotson RJ, Pias SC. Reduced Oxygen Permeability upon Protein Incorporation Within Phospholipid Bilayers. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2018; 1072:405-411. [PMID: 30178379 DOI: 10.1007/978-3-319-91287-5_65] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Intracellular oxygenation is key to energy metabolism as well as tumor radiation therapy. Although integral proteins are ubiquitous in membranes, few studies have considered their effects on molecular oxygen permeability. Published experimental work with rhodopsin and bacteriorhodopsin has led to the hypothesis that integral proteins lessen membrane oxygen permeability, as well as the permeability of the lipid region. The current work uses atomistic molecular dynamics simulations to test the influence of an ungated potassium channel protein on the oxygen permeability of palmitoyloleoylphosphatidylcholine (POPC) bilayers with and without cholesterol. Consistent with experiment, whole-membrane oxygen permeability is cut in half upon adding 30 wt% potassium channel protein to POPC, and the apparent permeability of the lipid portion of the membrane decreases by 40%. Unexpectedly, oxygen is found to interact directly with the protein surface, accompanied by a 40% reduction of the apparent whole-membrane diffusion coefficient. Similar effects are seen in systems combining the potassium channel with 1:1 POPC/cholesterol, but the magnitude of permeability reduction is smaller by ~30%. Overall, the simulations indicate that integral proteins can reduce oxygen permeability by altering the diffusional path and the local diffusivity. This effect may be especially important in the protein-dense membranes of mitochondria.
Collapse
Affiliation(s)
- Rachel J Dotson
- Department of Chemistry, New Mexico Institute of Mining and Technology (New Mexico Tech), Socorro, NM, USA
| | - Sally C Pias
- Department of Chemistry, New Mexico Institute of Mining and Technology (New Mexico Tech), Socorro, NM, USA.
| |
Collapse
|
187
|
Effect of Chain Unsaturation and Temperature on Oxygen Diffusion Through Lipid Membranes from Simulations. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2018; 1072:399-404. [DOI: 10.1007/978-3-319-91287-5_64] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/03/2022]
|
188
|
Atilhan M, Costa LT, Aparicio S. On the behaviour of aqueous solutions of deep eutectic solvents at lipid biomembranes. J Mol Liq 2017. [DOI: 10.1016/j.molliq.2017.09.082] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
189
|
|
190
|
Bakarić D, Petrov D, Mouvenchery YK, Heiβler S, Oostenbrink C, Schaumann GE. Ion-induced modification of the sucrose network and its impact on melting of freeze-dried liposomes. DSC and molecular dynamics study. Chem Phys Lipids 2017; 210:38-46. [PMID: 29179944 DOI: 10.1016/j.chemphyslip.2017.11.015] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2017] [Revised: 11/14/2017] [Accepted: 11/23/2017] [Indexed: 02/07/2023]
Abstract
Disaccharides play an important role in survival of anhydrobiotic organisms during extreme environmental conditions. A key protection feature is their capability to form the hydrogen bond (HB) network in a similar fashion as the one made by water. Since various ions also affect the HB network in completely hydrated systems, it is of a great interest to understand how they impact preservation when incorporated in a disaccharide network. To address this, we employ a combination of experimental and modeling techniques to study behavior of multilamellar 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) liposomes freeze-dried with sucrose in presence of NaCl or NaH2PO4·H2O at various concentrations (0.01-1M). Differential scanning calorimetry (DSC) was employed in order to determine the cooperative unit size (CUS), the number of lipid molecules that constitute a domain of cooperative motion in the liposome, and the melting temperature (Tm). In the absence of salt CUS was estimated to be 122±12, whereas in the presence of NaCl CUS increases more (347±34 for c=1M) than for NaH2PO4·H2O (193±26 for 1M). When it comes to Tm, the situation is reversed; NaCl induces increase by about 1K, while NaH2PO4·H2O by about 10K. These findings clearly demonstrate how different interaction forces-hydrogen bonding, charge pairing, and van der Waals interactions between acyl chains-affect CUS and Tm. Their interplay and contribution of particular interaction was further analyzed with molecular dynamics (MD) simulations. This analysis demonstrated that the HB network of DMPC and sucrose is partially disrupted in the presence of NaCl ions, and even to a greater extent in the case of NaH2PO4·H2O ions. Notably, H2PO4- ions outcompete and replace the sucrose molecules at the DMPC surface, which in turn alters the nature of the DMPC-surrounding interactions, from a weaker HB-dominated to a stronger CP-dominated interaction network.
Collapse
Affiliation(s)
- Danijela Bakarić
- University of Koblenz-Landau, Institute for Environmental Sciences, Group of Environmental and Soil Chemistry, Fortstraße 7, D-76829 Landau, Germany.
| | - Dražen Petrov
- Department of Material Sciences and Process Engineering, Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences Vienna, Muthgasse 18, A-1190 Vienna, Austria
| | - Yamuna Kunhi Mouvenchery
- University of Koblenz-Landau, Institute for Environmental Sciences, Group of Environmental and Soil Chemistry, Fortstraße 7, D-76829 Landau, Germany
| | - Stefan Heiβler
- Institute for Functional Interfaces, Karlsruhe Institute for Technology, Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
| | - Chris Oostenbrink
- Department of Material Sciences and Process Engineering, Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences Vienna, Muthgasse 18, A-1190 Vienna, Austria
| | - Gabriele E Schaumann
- University of Koblenz-Landau, Institute for Environmental Sciences, Group of Environmental and Soil Chemistry, Fortstraße 7, D-76829 Landau, Germany.
| |
Collapse
|
191
|
Doktorova M, Heberle FA, Kingston RL, Khelashvili G, Cuendet MA, Wen Y, Katsaras J, Feigenson GW, Vogt VM, Dick RA. Cholesterol Promotes Protein Binding by Affecting Membrane Electrostatics and Solvation Properties. Biophys J 2017; 113:2004-2015. [PMID: 29117524 DOI: 10.1016/j.bpj.2017.08.055] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2017] [Revised: 08/07/2017] [Accepted: 08/09/2017] [Indexed: 01/28/2023] Open
Abstract
Binding of the retroviral structural protein Gag to the cellular plasma membrane is mediated by the protein's matrix (MA) domain. Prominent among MA-PM interactions is electrostatic attraction between the positively charged MA domain and the negatively charged plasma membrane inner leaflet. Previously, we reported that membrane association of HIV-1 Gag, as well as purified Rous sarcoma virus (RSV) MA and Gag, depends strongly on the presence of acidic lipids and is enhanced by cholesterol (Chol). The mechanism underlying this enhancement was unclear. Here, using a broad set of in vitro and in silico techniques we addressed molecular mechanisms of association between RSV MA and model membranes, and investigated how Chol enhances this association. In neutron scattering experiments with liposomes in the presence or absence of Chol, MA preferentially interacted with preexisting POPS-rich clusters formed by nonideal lipid mixing, binding peripherally to the lipid headgroups with minimal perturbation to the bilayer structure. Molecular dynamics simulations showed a stronger MA-bilayer interaction in the presence of Chol, and a large Chol-driven increase in lipid packing and membrane surface charge density. Although in vitro MA-liposome association is influenced by disparate variables, including ionic strength and concentrations of Chol and charged lipids, continuum electrostatic theory revealed an underlying dependence on membrane surface potential. Together, these results conclusively show that Chol affects RSV MA-membrane association by making the electrostatic potential at the membrane surface more negative, while decreasing the penalty for lipid headgroup desolvation. The presented approach can be applied to other viral and nonviral proteins.
Collapse
Affiliation(s)
- Milka Doktorova
- Tri-Institutional PhD Program in Computational Biology and Medicine, Weill Cornell Medical College, New York, New York
| | - Frederick A Heberle
- The Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, Tennessee; Joint Institute for Biological Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee; Biology and Soft Matter Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee.
| | - Richard L Kingston
- School of Biological Sciences, The University of Auckland, Auckland, New Zealand
| | - George Khelashvili
- Department of Physiology and Biophysics, Weill Cornell Medical College, New York, New York
| | - Michel A Cuendet
- Department of Physiology and Biophysics, Weill Cornell Medical College, New York, New York
| | - Yi Wen
- Department of Biochemistry and Molecular Cell Biology, Cornell University, Ithaca, New York
| | - John Katsaras
- The Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, Tennessee; Biology and Soft Matter Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee; Department of Physics and Astronomy, University of Tennessee, Knoxville, Tennessee; Shull Wollan Center, Oak Ridge National Laboratory, Oak Ridge, Tennessee
| | - Gerald W Feigenson
- Department of Biochemistry and Molecular Cell Biology, Cornell University, Ithaca, New York
| | - Volker M Vogt
- Department of Biochemistry and Molecular Cell Biology, Cornell University, Ithaca, New York
| | - Robert A Dick
- Department of Biochemistry and Molecular Cell Biology, Cornell University, Ithaca, New York.
| |
Collapse
|
192
|
Zaki AM, Carbone P. How the Incorporation of Pluronic Block Copolymers Modulates the Response of Lipid Membranes to Mechanical Stress. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2017; 33:13284-13294. [PMID: 29084428 DOI: 10.1021/acs.langmuir.7b02244] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
Abstract
We employ atomistic molecular dynamics simulations to investigate the effect that the incorporation of the nonionic amphiphilic copolymer known as Pluronic L64 has on the mechanical stability of a DPPC membrane. The simulations reveal that the incorporation of the polymer chains leads to membranes that can sustain increasing mechanical stresses. Analysis of mechanical, structural, and dynamic properties of the membrane shows that the polymer chains interact strongly with the lipids in the vicinity, restraining their mobility and imparting better mechanical stability to the membrane. The hybrid membranes under tension remain thicker, more ordered, and stiffer in comparison to their lipid analogues. Trans-bilayer lipid movements (flip-flop) are observed and appear to be triggered by the presence of the polymer chains. A careful analysis of the pore formation under high tensions reveals two distinctive mechanisms that depend on the distribution of the hydrophilic polymer blocks in the bilayer. Finally, the rate of growth of the formed membrane defects is slowed down in the presence of polymers. These findings show that Pluronic block copolymers could be exploited for the formation of optimized hybrid nanodevices with controlled elastic and dynamic properties.
Collapse
Affiliation(s)
- Afroditi Maria Zaki
- School of Chemical Engineering and Analytical Science, The University of Manchester , Oxford Road, Manchester M13 9PL, United Kingdom
| | - Paola Carbone
- School of Chemical Engineering and Analytical Science, The University of Manchester , Oxford Road, Manchester M13 9PL, United Kingdom
| |
Collapse
|
193
|
Doktorova M, Heberle FA, Kingston RL, Khelashvili G, Cuendet MA, Wen Y, Katsaras J, Feigenson GW, Vogt VM, Dick RA. Cholesterol Promotes Protein Binding by Affecting Membrane Electrostatics and Solvation Properties. Biophys J 2017. [PMID: 29117524 DOI: 10.1016/j.bpj.2017.08.055.] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022] Open
Abstract
Binding of the retroviral structural protein Gag to the cellular plasma membrane is mediated by the protein's matrix (MA) domain. Prominent among MA-PM interactions is electrostatic attraction between the positively charged MA domain and the negatively charged plasma membrane inner leaflet. Previously, we reported that membrane association of HIV-1 Gag, as well as purified Rous sarcoma virus (RSV) MA and Gag, depends strongly on the presence of acidic lipids and is enhanced by cholesterol (Chol). The mechanism underlying this enhancement was unclear. Here, using a broad set of in vitro and in silico techniques we addressed molecular mechanisms of association between RSV MA and model membranes, and investigated how Chol enhances this association. In neutron scattering experiments with liposomes in the presence or absence of Chol, MA preferentially interacted with preexisting POPS-rich clusters formed by nonideal lipid mixing, binding peripherally to the lipid headgroups with minimal perturbation to the bilayer structure. Molecular dynamics simulations showed a stronger MA-bilayer interaction in the presence of Chol, and a large Chol-driven increase in lipid packing and membrane surface charge density. Although in vitro MA-liposome association is influenced by disparate variables, including ionic strength and concentrations of Chol and charged lipids, continuum electrostatic theory revealed an underlying dependence on membrane surface potential. Together, these results conclusively show that Chol affects RSV MA-membrane association by making the electrostatic potential at the membrane surface more negative, while decreasing the penalty for lipid headgroup desolvation. The presented approach can be applied to other viral and nonviral proteins.
Collapse
Affiliation(s)
- Milka Doktorova
- Tri-Institutional PhD Program in Computational Biology and Medicine, Weill Cornell Medical College, New York, New York
| | - Frederick A Heberle
- The Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, Tennessee; Joint Institute for Biological Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee; Biology and Soft Matter Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee.
| | - Richard L Kingston
- School of Biological Sciences, The University of Auckland, Auckland, New Zealand
| | - George Khelashvili
- Department of Physiology and Biophysics, Weill Cornell Medical College, New York, New York
| | - Michel A Cuendet
- Department of Physiology and Biophysics, Weill Cornell Medical College, New York, New York
| | - Yi Wen
- Department of Biochemistry and Molecular Cell Biology, Cornell University, Ithaca, New York
| | - John Katsaras
- The Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, Tennessee; Biology and Soft Matter Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee; Department of Physics and Astronomy, University of Tennessee, Knoxville, Tennessee; Shull Wollan Center, Oak Ridge National Laboratory, Oak Ridge, Tennessee
| | - Gerald W Feigenson
- Department of Biochemistry and Molecular Cell Biology, Cornell University, Ithaca, New York
| | - Volker M Vogt
- Department of Biochemistry and Molecular Cell Biology, Cornell University, Ithaca, New York
| | - Robert A Dick
- Department of Biochemistry and Molecular Cell Biology, Cornell University, Ithaca, New York.
| |
Collapse
|
194
|
Nagao M, Kelley EG, Ashkar R, Bradbury R, Butler PD. Probing Elastic and Viscous Properties of Phospholipid Bilayers Using Neutron Spin Echo Spectroscopy. J Phys Chem Lett 2017; 8:4679-4684. [PMID: 28892394 DOI: 10.1021/acs.jpclett.7b01830] [Citation(s) in RCA: 82] [Impact Index Per Article: 11.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/08/2023]
Abstract
The elastic and viscous properties of self-assembled amphiphilic membranes dictate the intricate hierarchy of their structure and dynamics ranging from the diffusion of individual molecules to the large-scale deformation of the membrane. We previously demonstrated that neutron spin echo spectroscopy measurements of model amphiphilic membranes can access the naturally occurring submicrosecond membrane motions, such as bending and thickness fluctuations. Here we show how the experimentally measured fluctuation parameters can be used to determine the inherent membrane properties and demonstrate how membrane viscosity and compressibility modulus are influenced by lipid composition in a series of simple phosphatidylcholine bilayers with different tail lengths as a function of temperature. This approach highlights the interdependence of the bilayer elastic and viscous properties and the collective membrane dynamics and opens new avenues to investigating the mechanical properties of more complex and biologically inspired systems.
Collapse
Affiliation(s)
- Michihiro Nagao
- NIST Center for Neutron Research, National Institute of Standards and Technology , Gaithersburg, Maryland 20899, United States
- Center for Exploration of Energy and Matter, Department of Physics, Indiana University , Bloomington, Indiana 47408, United States
| | - Elizabeth G Kelley
- NIST Center for Neutron Research, National Institute of Standards and Technology , Gaithersburg, Maryland 20899, United States
| | - Rana Ashkar
- Biology and Soft Matter Division, Oak Ridge National Laboratory , Oak Ridge, Tennessee 37831, United States
| | - Robert Bradbury
- NIST Center for Neutron Research, National Institute of Standards and Technology , Gaithersburg, Maryland 20899, United States
- Center for Exploration of Energy and Matter, Department of Physics, Indiana University , Bloomington, Indiana 47408, United States
| | - Paul D Butler
- NIST Center for Neutron Research, National Institute of Standards and Technology , Gaithersburg, Maryland 20899, United States
- Department of Chemical & Biomolecular Engineering, University of Delaware , Newark, Delaware 19716, United States
| |
Collapse
|
195
|
Beaven AH, Sodt AJ, Pastor RW, Koeppe RE, Andersen OS, Im W. Characterizing Residue-Bilayer Interactions Using Gramicidin A as a Scaffold and Tryptophan Substitutions as Probes. J Chem Theory Comput 2017; 13:5054-5064. [PMID: 28870079 DOI: 10.1021/acs.jctc.7b00400] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Previous experiments have shown that the lifetime of a gramicidin A dimer channel (which forms from two nonconducting monomers) in a lipid bilayer is modulated by mutations of the tryptophan (Trp) residues at the bilayer-water interface. We explore this further using extensive molecular dynamics simulations of various gA dimer and monomer mutants at the Trp positions in phosphatidylcholine bilayers with different tail lengths. gA interactions with the surrounding bilayer are strongly modulated by mutating these Trp residues. There are three principal effects: eliminating residue hydrogen bonding ability (i.e., reducing the channel-monolayer coupling strength) reduces the extent of the bilayer deformation caused by the assembled dimeric channel; a residue's size and geometry affects its orientation, leading to different hydrogen bonding partners; and increasing a residue's hydrophobicity increases the depth of gA monomer insertion relative to the bilayer center, thereby increasing the lipid bending frustration.
Collapse
Affiliation(s)
- Andrew H Beaven
- Department of Chemistry, The University of Kansas , Lawrence, Kansas 66045, United States
| | | | | | - Roger E Koeppe
- Department of Chemistry and Biochemistry, University of Arkansas , Fayetteville, Arkansas 72701, United States
| | - Olaf S Andersen
- Department of Physiology and Biophysics, Weill Cornell Medical College , New York, New York 10065, United States
| | - Wonpil Im
- Departments of Biological Sciences and Bioengineering, Lehigh University , Bethlehem, Pennsylvania 18015, United States
| |
Collapse
|
196
|
Molugu TR, Lee S, Brown MF. Concepts and Methods of Solid-State NMR Spectroscopy Applied to Biomembranes. Chem Rev 2017; 117:12087-12132. [PMID: 28906107 DOI: 10.1021/acs.chemrev.6b00619] [Citation(s) in RCA: 70] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Concepts of solid-state NMR spectroscopy and applications to fluid membranes are reviewed in this paper. Membrane lipids with 2H-labeled acyl chains or polar head groups are studied using 2H NMR to yield knowledge of their atomistic structures in relation to equilibrium properties. This review demonstrates the principles and applications of solid-state NMR by unifying dipolar and quadrupolar interactions and highlights the unique features offered by solid-state 2H NMR with experimental illustrations. For randomly oriented multilamellar lipids or aligned membranes, solid-state 2H NMR enables direct measurement of residual quadrupolar couplings (RQCs) due to individual C-2H-labeled segments. The distribution of RQC values gives nearly complete profiles of the segmental order parameters SCD(i) as a function of acyl segment position (i). Alternatively, one can measure residual dipolar couplings (RDCs) for natural abundance lipid samples to obtain segmental SCH order parameters. A theoretical mean-torque model provides acyl-packing profiles representing the cumulative chain extension along the normal to the aqueous interface. Equilibrium structural properties of fluid bilayers and various thermodynamic quantities can then be calculated, which describe the interactions with cholesterol, detergents, peptides, and integral membrane proteins and formation of lipid rafts. One can also obtain direct information for membrane-bound peptides or proteins by measuring RDCs using magic-angle spinning (MAS) in combination with dipolar recoupling methods. Solid-state NMR methods have been extensively applied to characterize model membranes and membrane-bound peptides and proteins, giving unique information on their conformations, orientations, and interactions in the natural liquid-crystalline state.
Collapse
Affiliation(s)
- Trivikram R Molugu
- Department of Chemistry & Biochemistry and ‡Department of Physics, University of Arizona , Tucson, Arizona 85721, United States
| | - Soohyun Lee
- Department of Chemistry & Biochemistry and ‡Department of Physics, University of Arizona , Tucson, Arizona 85721, United States
| | - Michael F Brown
- Department of Chemistry & Biochemistry and ‡Department of Physics, University of Arizona , Tucson, Arizona 85721, United States
| |
Collapse
|
197
|
Nagle JF. X-ray scattering reveals molecular tilt is an order parameter for the main phase transition in a model biomembrane. Phys Rev E 2017; 96:030401. [PMID: 29346876 DOI: 10.1103/physreve.96.030401] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2017] [Indexed: 11/07/2022]
Abstract
Synchrotron diffuse x-ray scattering data reveal a dramatic softening of the molecular tilt modulus K_{θ} of the model biomembrane composed of DMPC lipids as the temperature is lowered towards the main phase transition temperature at T_{M}=24^{∘}C. Spontaneous tilt occurs below T_{M}, suggesting that tilt is a symmetry breaking order parameter. Consistent with this hypothesis, it is also found that a different lipid POPS has no spontaneous tilt below its T_{M} at 14^{∘}C and correspondingly its tilt modulus did not soften as T_{M} was approached from above. As previously known, the bending modulus K_{C} of DMPC also softens close to T_{M}, but unlike the tilt modulus, K_{C} has a maximum 3^{∘} above T_{M}, which also marks the limit of the well-known anomalous swelling regime. Tilt adds a different perspective to our previous understanding of the main phase transition in lipid bilayers.
Collapse
Affiliation(s)
- John F Nagle
- Department of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
| |
Collapse
|
198
|
Argudo D, Bethel NP, Marcoline FV, Wolgemuth CW, Grabe M. New Continuum Approaches for Determining Protein-Induced Membrane Deformations. Biophys J 2017; 112:2159-2172. [PMID: 28538153 DOI: 10.1016/j.bpj.2017.03.040] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2016] [Revised: 03/16/2017] [Accepted: 03/27/2017] [Indexed: 01/21/2023] Open
Abstract
The influence of the membrane on transmembrane proteins is central to a number of biological phenomena, notably the gating of stretch activated ion channels. Conversely, membrane proteins can influence the bilayer, leading to the stabilization of particular membrane shapes, topological changes that occur during vesicle fission and fusion, and shape-dependent protein aggregation. Continuum elastic models of the membrane have been widely used to study protein-membrane interactions. These mathematical approaches produce physically interpretable membrane shapes, energy estimates for the cost of deformation, and a snapshot of the equilibrium configuration. Moreover, elastic models are much less computationally demanding than fully atomistic and coarse-grained simulation methodologies; however, it has been argued that continuum models cannot reproduce the distortions observed in fully atomistic molecular dynamics simulations. We suggest that this failure can be overcome by using chemically and geometrically accurate representations of the protein. Here, we present a fast and reliable hybrid continuum-atomistic model that couples the protein to the membrane. We show that the model is in excellent agreement with fully atomistic simulations of the ion channel gramicidin embedded in a POPC membrane. Our continuum calculations not only reproduce the membrane distortions produced by the channel but also accurately determine the channel's orientation. Finally, we use our method to investigate the role of membrane bending around the charged voltage sensors of the transient receptor potential cation channel TRPV1. We find that membrane deformation significantly stabilizes the energy of insertion of TRPV1 by exposing charged residues on the S4 segment to solution.
Collapse
Affiliation(s)
- David Argudo
- Cardiovascular Research Institute, Department of Pharmaceutical Chemistry, University of California San Francisco, San Francisco, California
| | - Neville P Bethel
- Cardiovascular Research Institute, Department of Pharmaceutical Chemistry, University of California San Francisco, San Francisco, California
| | - Frank V Marcoline
- Cardiovascular Research Institute, Department of Pharmaceutical Chemistry, University of California San Francisco, San Francisco, California
| | - Charles W Wolgemuth
- Departments of Molecular and Cellular Biology and Physics, University of Arizona, Tucson, Arizona
| | - Michael Grabe
- Cardiovascular Research Institute, Department of Pharmaceutical Chemistry, University of California San Francisco, San Francisco, California.
| |
Collapse
|
199
|
Li H, Chowdhary J, Huang L, He X, MacKerell AD, Roux B. Drude Polarizable Force Field for Molecular Dynamics Simulations of Saturated and Unsaturated Zwitterionic Lipids. J Chem Theory Comput 2017; 13:4535-4552. [PMID: 28731702 DOI: 10.1021/acs.jctc.7b00262] [Citation(s) in RCA: 70] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Additive force fields are designed to account for induced electronic polarization in a mean-field average way, using effective empirical fixed charges. The limitation of this approximation is cause for serious concerns, particularly in the case of lipid membranes, where the molecular environment undergoes dramatic variations over microscopic length scales. A polarizable force field based on the classical Drude oscillator offers a practical and computationally efficient framework for an improved representation of electrostatic interactions in molecular simulations. Building on the first-generation Drude polarizable force field for the dipalmitoylphosphatidylcholine 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) molecule, the present effort was undertaken to improve this initial model and expand the force field to a wider range of phospholipid molecules. New lipids parametrized include dimyristoylphosphatidylcholine (DMPC), dilauroylphosphatidylcholine (DLPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), dipalmitoylphosphatidylethanolamine (DPPE), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE), and 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine (DOPE). The iterative optimization protocol employed in this effort led to lipid models that achieve a good balance between reproducing quantum mechanical data on model compound representative of phospholipids and reproducing a range of experimental condensed phase properties of bilayers. A parametrization strategy based on a restrained ensemble-maximum entropy methodology was used to help accurately match the experimental NMR order parameters in the polar headgroup region. All the parameters were developed to be compatible with the remainder of the Drude polarizable force field, which includes water, ions, proteins, DNA, and selected carbohydrates.
Collapse
Affiliation(s)
- Hui Li
- Department of Biochemistry and Molecular Biology, Gordon Center for Integrative Science, University of Chicago , Chicago, Illinois 60637, United States
| | - Janamejaya Chowdhary
- Department of Biochemistry and Molecular Biology, Gordon Center for Integrative Science, University of Chicago , Chicago, Illinois 60637, United States
| | - Lei Huang
- Department of Biochemistry and Molecular Biology, Gordon Center for Integrative Science, University of Chicago , Chicago, Illinois 60637, United States
| | - Xibing He
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore , Baltimore, Maryland 21201, United States
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore , Baltimore, Maryland 21201, United States
| | - Benoît Roux
- Department of Biochemistry and Molecular Biology, Gordon Center for Integrative Science, University of Chicago , Chicago, Illinois 60637, United States
| |
Collapse
|
200
|
Simulations of Membrane-Disrupting Peptides I: Alamethicin Pore Stability and Spontaneous Insertion. Biophys J 2017; 111:1248-1257. [PMID: 27653483 DOI: 10.1016/j.bpj.2016.08.014] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2016] [Revised: 08/09/2016] [Accepted: 08/10/2016] [Indexed: 11/23/2022] Open
Abstract
An all-atom molecular dynamics simulation of the archetype barrel-stave alamethicin (alm) pore in a 1,2-dioleoyl-sn-glycero-3-phosphocholine bilayer at 313 K indicates that ∼7 μs is required for equilibration of a preformed 6-peptide pore; the pore remains stable for the duration of the remaining 7 μs of the trajectory, and the structure factors agree well with experiment. A 5 μs simulation of 10 surface-bound alm peptides shows significant peptide unfolding and some unbinding, but no insertion. Simulations at 363 and 413 K with a -0.2 V electric field yield peptide insertion in 1 μs. Insertion is initiated by the folding of residues 3-11 into an α-helix, and mediated by membrane water or by previously inserted peptides. The stability of five alm pore peptides at 413 K with a -0.2 V electric field demonstrates a significant preference for a transmembrane orientation. Hence, and in contrast to the cationic antimicrobial peptide described in the following article, alm shows a strong preference for the inserted over the surface-bound state.
Collapse
|