1
|
Zheng L, Bandara SR, Tan Z, Leal C. Lipid nanoparticle topology regulates endosomal escape and delivery of RNA to the cytoplasm. Proc Natl Acad Sci U S A 2023; 120:e2301067120. [PMID: 37364130 PMCID: PMC10318962 DOI: 10.1073/pnas.2301067120] [Citation(s) in RCA: 28] [Impact Index Per Article: 28.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2023] [Accepted: 05/04/2023] [Indexed: 06/28/2023] Open
Abstract
RNA therapeutics have the potential to resolve a myriad of genetic diseases. Lipid nanoparticles (LNPs) are among the most successful RNA delivery systems. Expanding their use for the treatment of more genetic diseases hinges on our ability to continuously evolve the design of LNPs with high potency, cellular-specific targeting, and low side effects. Overcoming the difficulty of releasing cargo from endocytosed LNPs remains a significant hurdle. Here, we investigate the fundamental properties of nonviral RNA nanoparticles pertaining to the activation of topological transformations of endosomal membranes and RNA translocation into the cytosol. We show that, beyond composition, LNP fusogenicity can be prescribed by designing LNP nanostructures that lower the energetic cost of fusion and fusion-pore formation with a target membrane. The inclusion of structurally active lipids leads to enhanced LNP endosomal fusion, fast evasion of endosomal entrapment, and efficacious RNA delivery. For example, conserving the lipid make-up, RNA-LNPs having cuboplex nanostructures are significantly more efficacious at endosomal escape than traditional lipoplex constructs.
Collapse
Affiliation(s)
- Lining Zheng
- Department of Materials Science and Engineering, University of Illinois Urbana-Champaign, Urbana, IL61801
| | - Sarith R. Bandara
- Department of Materials Science and Engineering, University of Illinois Urbana-Champaign, Urbana, IL61801
| | - Zhengzhong Tan
- Department of Materials Science and Engineering, University of Illinois Urbana-Champaign, Urbana, IL61801
| | - Cecilia Leal
- Department of Materials Science and Engineering, University of Illinois Urbana-Champaign, Urbana, IL61801
| |
Collapse
|
2
|
Siegel DP. Bicontinuous inverted cubic phase stabilization as an index of antimicrobial and membrane fusion peptide activity. BIOCHIMICA ET BIOPHYSICA ACTA. BIOMEMBRANES 2022; 1864:183815. [PMID: 34748744 DOI: 10.1016/j.bbamem.2021.183815] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/20/2021] [Revised: 10/06/2021] [Accepted: 10/29/2021] [Indexed: 02/07/2023]
Abstract
Some antimicrobial peptides (AMPs) and membrane fusion-catalyzing peptides (FPs) stabilize bicontinuous inverted cubic (QII) phases. Previous authors proposed a topological rationale: since AMP-induced pores, fusion intermediates, and QII phases all have negative Gaussian curvature (NGC), peptides which produce NGC in one structure also do it in another. This assumes that peptides change the curvature energy of the lipid membranes. Here I test this with a Helfrich curvature energy model. First, experimentally, I show that lipid systems often used to study peptide NGC have NGC without peptides at higher temperatures. To determine the net effect of an AMP on NGC, the equilibrium phase behavior of the host lipids must be determined. Second, the model shows that AMPs must make large changes in the curvature energy to stabilize AMP-induced pores. Peptide-induced changes in elastic constants affect pores and QII phase differently. Changes in spontaneous curvature affect them in opposite ways. The observed correlation between QII phase stabilization and AMP activity doesn't show that AMPs act by lowering pore curvature energy. A different rationale is proposed. In theory, AMPs could simultaneously stabilize QII phase and pores by drastically changing two particular elastic constants. This could be tested by measuring AMP effects on the individual constants. I propose experiments to do that. Unlike AMPs, FPs must make only small changes in the curvature energy to catalyze fusion. It they act in this way, their fusion activity should correlate with their ability to stabilize QII phases.
Collapse
Affiliation(s)
- D P Siegel
- Givaudan Inc., 1199 Edison Drive, Cincinnati, OH 45216, United States of America.
| |
Collapse
|
3
|
Beaven AH, Arnarez C, Lyman E, Bennett WFD, Sodt AJ. Curvature Energetics Determined by Alchemical Simulation on Four Topologically Distinct Lipid Phases. J Phys Chem B 2021; 125:1815-1824. [PMID: 33570958 PMCID: PMC9069320 DOI: 10.1021/acs.jpcb.0c09458] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The relative curvature energetics of two lipids are tested using thermodynamic integration (TI) on four topologically distinct lipid phases. Simulations use TI to switch between choline headgroup lipids (POPC; that prefers to be flat) and ethanolamine headgroup lipids (POPE; that prefer, for example, the inner monolayer of vesicles). The thermodynamical moving of the lipids between planar, inverse hexagonal (HII), cubic (QII; Pn3m space group), and vesicle topologies reveals differences in material parameters that were previously challenging to access. The methodology allows for predictions of two important lipid material properties: the difference in POPC/POPE monolayer intrinsic curvature (ΔJ0) and the difference in POPC/POPE monolayer Gaussian curvature modulus (Δκ̅m), both of which are connected to the energetics of topological variation. Analysis of the TI data indicates that, consistent with previous experiment and simulation, the J0 of POPE is more negative than POPC (ΔJ0 = -0.018 ± 0.001 Å-1). The theoretical framework extracts significant differences in κ̅m of which POPE is less negative than POPC by 2.0 to 4.0 kcal/mol. The range of these values is determined by considering subsets of the simulations, and disagreement between these subsets suggests separate mechanical parameters at very high curvature. Finally, the fit of the TI data to the model indicates that the position of the pivotal plane of curvature is not constant across topologies at high curvature. Overall, the results offer insights into lipid material properties, the limits of a single HC model, and how to test them using simulation.
Collapse
Affiliation(s)
- Andrew H Beaven
- Unit on Membrane Chemical Physics, Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, Maryland 20892, United States.,Postdoctoral Research Associate Program, National Institute of General Medical Sciences, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Clément Arnarez
- Department of Physics and Astronomy, Department of Chemistry and Biochemistry, University of Delaware, Newark, Delaware 19716, United States
| | - Edward Lyman
- Department of Physics and Astronomy, Department of Chemistry and Biochemistry, University of Delaware, Newark, Delaware 19716, United States
| | - W F Drew Bennett
- Biosciences and Biotechnology Division, Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Alexander J Sodt
- Unit on Membrane Chemical Physics, Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, Maryland 20892, United States
| |
Collapse
|
4
|
Abstract
Vesicle structures primarily embody spherical capsules composed of a single or multiple bilayers, entrapping a pool of aqueous solution in their interior. The bilayers can be synthesised by phospholipids or other amphiphiles (surfactants, block copolymers, etc.). Vesicles with broad-spectrum applications in numerous scientific disciplines, including biochemistry, biophysics, biology, and various pharmaceutical industries, have attracted widespread attention. Consequently, a multitude of protocols have been devised and proposed for their fabrication. In this review, with a motivation to derive the basic conditions for the formation of vesicles, the associated thermodynamic and kinetic aspects are comprehensively appraised. Contextually, an all-purpose overview of the underlying thermodynamics of bilayer/membrane generation and deformation, including the chemical potential of aggregates, geometric packing and the concept of elastic properties, is presented. Additionally, the current review highlights the probable, inherent mechanisms of vesicle formation under distinct modes of manufacturing. We lay focus on vesicle formation from pre-existing bilayers, as well as from bilayers, which form when lipids from an organic solvent are transferred into an aqueous medium. Furthermore, we outline the kinetic effects on vesicle formation from the lamellar phase, with and without the presence of shearing force. Wherever required, the experimental and/or theoretical outcomes, the driving forces for vesicle size selection, and various scaling laws are also reviewed, all of which facilitate an overall improved understanding of the vesicle formation mechanisms.
Collapse
Affiliation(s)
- Chandra Has
- Department of Chemical Engineering, Indian Institute of Technology Bombay, Mumbai, India
| | - Sharadwata Pan
- School of Life Sciences Weihenstephan, Technical University of Munich, Freising, Germany
| |
Collapse
|
5
|
Terzi MM, Ergüder MF, Deserno M. A consistent quadratic curvature-tilt theory for fluid lipid membranes. J Chem Phys 2019; 151:164108. [PMID: 31675861 DOI: 10.1063/1.5119683] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
The tilt of a lipid molecule describes the deviation of its orientation away from the local normal of its embedding membrane. Tilt is the subleading degree of freedom after a membrane's geometry, and it becomes relevant at scales comparable to lipid bilayer thickness. Building on earlier work by Hamm and Kozlov [Eur. Phys. J. E 3, 323 (2000)], who envisioned lipid membranes as thin prestressed fluid elastic films, and Terzi and Deserno [J. Chem. Phys. 147, 084702 (2017)], who discovered a new coupling term between splay and tilt divergence, we construct a theory of membrane elasticity that is quadratic in geometry and tilt and complete at order 1/length2. We show that a general and consistent treatment of both lateral and transverse depth-dependent shear stresses creates several contributions to the elastic energy density, of which only a subset had previously been identified. Apart from the well-known penalty of lipid twist (the curl of tilt), these terms generate no qualitatively new phenomenology, but they quantitatively revise the connections between the moduli of a tilt-curvature theory and its underlying microscopic foundation. In particular, we argue that the monolayer Gaussian curvature modulus κ¯m, widely believed to be equal to the second moment of the transmonolayer stress profile, acquires a second contribution from lipid twist, which is always negative. This could resolve the long-standing conundrum that many measured values of κ¯m appeared to have a sign that violates basic stability considerations. We also show that the previously discovered novel coupling between splay and tilt divergence is not simply proportional to κ¯m but acquires its own splay-tilt coupling modulus, κst,m. We explore the predictions of our theory for various elastic moduli and their mutual interrelations and use an extensive set of existing atomistic molecular dynamics simulations for 12 different lipid types to collectively reason about such predictions. We find that bending rigidities are captured fairly well by existing theories, while reliable predictions for local moduli, especially the splay-tilt coupling modulus, remain challenging.
Collapse
Affiliation(s)
- M Mert Terzi
- Department of Physics, Carnegie Mellon University, 5000 Forbes Ave., Pittsburgh, Pennsylvania 15213, USA
| | - Muhammed F Ergüder
- 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
|
6
|
Zelisko M, Ahmadpoor F, Gao H, Sharma P. Determining the Gaussian Modulus and Edge Properties of 2D Materials: From Graphene to Lipid Bilayers. PHYSICAL REVIEW LETTERS 2017; 119:068002. [PMID: 28949621 DOI: 10.1103/physrevlett.119.068002] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/18/2017] [Indexed: 06/07/2023]
Abstract
The dominant deformation behavior of two-dimensional materials (bending) is primarily governed by just two parameters: bending rigidity and the Gaussian modulus. These properties also set the energy scale for various important physical and biological processes such as pore formation, cell fission and generally, any event accompanied by a topological change. Unlike the bending rigidity, the Gaussian modulus is, however, notoriously difficult to evaluate via either experiments or atomistic simulations. In this Letter, recognizing that the Gaussian modulus and edge tension play a nontrivial role in the fluctuations of a 2D material edge, we derive closed-form expressions for edge fluctuations. Combined with atomistic simulations, we use the developed approach to extract the Gaussian modulus and edge tension at finite temperatures for both graphene and various types of lipid bilayers. Our results possibly provide the first reliable estimate of this elusive property at finite temperatures and appear to suggest that earlier estimates must be revised. In particular, we show that, if previously estimated properties are employed, the graphene-free edge will exhibit unstable behavior at room temperature. Remarkably, in the case of graphene, we show that the Gaussian modulus and edge tension even change sign at finite temperatures.
Collapse
Affiliation(s)
- Matthew Zelisko
- Department of Mechanical Engineering, University of Houston, Houston, Texas 77204, USA
| | - Fatemeh Ahmadpoor
- Department of Mechanical Engineering, University of Houston, Houston, Texas 77204, USA
- School of Engineering, Brown University, Providence, Rhode Island 02912, USA
| | - Huajian Gao
- School of Engineering, Brown University, Providence, Rhode Island 02912, USA
| | - Pradeep Sharma
- Department of Mechanical Engineering, University of Houston, Houston, Texas 77204, USA
- Department of Physics, University of Houston, Houston, Texas 77204, USA
| |
Collapse
|
7
|
Bubnis G, Risselada HJ, Grubmüller H. Exploiting Lipid Permutation Symmetry to Compute Membrane Remodeling Free Energies. PHYSICAL REVIEW LETTERS 2016; 117:188102. [PMID: 27834997 DOI: 10.1103/physrevlett.117.188102] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/29/2016] [Indexed: 05/26/2023]
Abstract
A complete physical description of membrane remodeling processes, such as fusion or fission, requires knowledge of the underlying free energy landscapes, particularly in barrier regions involving collective shape changes, topological transitions, and high curvature, where Canham-Helfrich (CH) continuum descriptions may fail. To calculate these free energies using atomistic simulations, one must address not only the sampling problem due to high free energy barriers, but also an orthogonal sampling problem of combinatorial complexity stemming from the permutation symmetry of identical lipids. Here, we solve the combinatorial problem with a permutation reduction scheme to map a structural ensemble into a compact, nondegenerate subregion of configuration space, thereby permitting straightforward free energy calculations via umbrella sampling. We applied this approach, using a coarse-grained lipid model, to test the CH description of bending and found sharp increases in the bending modulus for curvature radii below 10 nm. These deviations suggest that an anharmonic bending term may be required for CH models to give quantitative energetics of highly curved states.
Collapse
Affiliation(s)
- Greg Bubnis
- Department of Theoretical and Computational Biophysics, Max-Planck-Institute for Biophysical Chemistry, Göttingen 37077, Germany
| | - Herre Jelger Risselada
- Department of Theoretical and Computational Biophysics, Max-Planck-Institute for Biophysical Chemistry, Göttingen 37077, Germany
- Chemistry Department, Leibniz Institute of Surface Modification, Leipzig 04318, Germany
- Deptartment of Theoretical Physics, Georg-August University Göttingen, Göttingen 37077, Germany
| | - Helmut Grubmüller
- Department of Theoretical and Computational Biophysics, Max-Planck-Institute for Biophysical Chemistry, Göttingen 37077, Germany
| |
Collapse
|
8
|
Rangamani P, Mandadap KK, Oster G. Protein-induced membrane curvature alters local membrane tension. Biophys J 2015; 107:751-762. [PMID: 25099814 DOI: 10.1016/j.bpj.2014.06.010] [Citation(s) in RCA: 62] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2014] [Revised: 06/01/2014] [Accepted: 06/09/2014] [Indexed: 10/24/2022] Open
Abstract
Adsorption of proteins onto membranes can alter the local membrane curvature. This phenomenon has been observed in biological processes such as endocytosis, tubulation, and vesiculation. However, it is not clear how the local surface properties of the membrane, such as membrane tension, change in response to protein adsorption. In this article, we show that the partial differential equations arising from classical elastic model of lipid membranes, which account for simultaneous changes in shape and membrane tension due to protein adsorption in a local region, cannot be solved for nonaxisymmetric geometries using straightforward numerical techniques; instead, a viscous-elastic formulation is necessary to fully describe the system. Therefore, we develop a viscous-elastic model for inhomogeneous membranes of the Helfrich type. Using the newly available viscous-elastic model, we find that the lipids flow to accommodate changes in membrane curvature during protein adsorption. We show that, at the end of protein adsorption process, the system sustains a residual local tension to balance the difference between the actual mean curvature and the imposed spontaneous curvature. We also show that this change in membrane tension can have a functional impact such as altered response to pulling forces in the presence of proteins.
Collapse
Affiliation(s)
- Padmini Rangamani
- Department of Molecular and Cell Biology, University of California at Berkeley, Berkeley, California
| | - Kranthi K Mandadap
- Department of Molecular and Cell Biology, University of California at Berkeley, Berkeley, California; Department of Chemistry, University of California at Berkeley, Berkeley, California
| | - George Oster
- Department of Molecular and Cell Biology, University of California at Berkeley, Berkeley, California.
| |
Collapse
|
9
|
Deserno M. Fluid lipid membranes: From differential geometry to curvature stresses. Chem Phys Lipids 2015; 185:11-45. [DOI: 10.1016/j.chemphyslip.2014.05.001] [Citation(s) in RCA: 137] [Impact Index Per Article: 15.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2014] [Revised: 04/21/2014] [Accepted: 05/06/2014] [Indexed: 10/25/2022]
|
10
|
Hu M, de Jong DH, Marrink SJ, Deserno M. Gaussian curvature elasticity determined from global shape transformations and local stress distributions: a comparative study using the MARTINI model. Faraday Discuss 2013; 161:365-82; discussion 419-59. [PMID: 23805750 DOI: 10.1039/c2fd20087b] [Citation(s) in RCA: 77] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We calculate the Gaussian curvature modulus kappa of a systematically coarse-grained (CG) one-component lipid membrane by applying the method recently proposed by Hu et al. [Biophys. J., 2012, 102, 1403] to the MARTINI representation of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC). We find the value kappa/kappa = -1.04 +/- 0.03 for the elastic ratio between the Gaussian and the mean curvature modulus and deduce kappa(m)/kappa(m) = -0.98 +/- 0.09 for the monolayer elastic ratio, where the latter is based on plausible assumptions for the distance z0 of the monolayer neutral surface from the bilayer midplane and the spontaneous lipid curvature K(0m). By also analyzing the lateral stress profile sigma0(z) of our system, two other lipid types and pertinent data from the literature, we show that determining K(0m) and kappa through the first and second moment of sigma0(z) gives rise to physically implausible values for these observables. This discrepancy, which we previously observed for a much simpler CG model, suggests that the moment conditions derived from simple continuum assumptions miss the effect of physically important correlations in the lipid bilayer.
Collapse
Affiliation(s)
- Mingyang Hu
- Carnegie Mellon University, 5000 Forbes Avenue, Pittsburgh, PA 15213, USA
| | | | | | | |
Collapse
|
11
|
Coupling Membrane Elasticity and Structure to Protein Function. ACTA ACUST UNITED AC 2013. [DOI: 10.1016/b978-0-12-411515-6.00004-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register]
|
12
|
Computational Studies of Biomembrane Systems: Theoretical Considerations, Simulation Models, and Applications. FROM SINGLE MOLECULES TO NANOSCOPICALLY STRUCTURED MATERIALS 2013. [DOI: 10.1007/12_2013_258] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
13
|
|
14
|
Boulgaropoulos B, Rappolt M, Sartori B, Amenitsch H, Pabst G. Lipid sorting by ceramide and the consequences for membrane proteins. Biophys J 2012; 102:2031-8. [PMID: 22824266 DOI: 10.1016/j.bpj.2012.03.059] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2012] [Revised: 03/08/2012] [Accepted: 03/23/2012] [Indexed: 01/09/2023] Open
Abstract
We mimicked the effect of sphingomyelinase activity on lipid mixtures of palmitoyl-oleoyl-phosphatidylcholine, sphingomyelin, ceramide, and 10 mol % cholesterol. Using x-ray diffraction experiments in combination with osmotic stress we found, in agreement with previous studies, that ceramide induces a coexistence of L(α) and L(β) domains. A detailed structural analysis of the coexisting domains demonstrated an increase of lipid packing density and membrane thickness in the L(α) domains upon increasing overall ceramide levels. This provides evidence for a ceramide-driven accumulation of cholesterol in the L(α) domains, in support of previous reports. We further determined the bending rigidities of the coexisting domains and found that the accumulation of cholesterol in the L(α) domains stabilizes their bending rigidity, which experiences a dramatic drop in the absence of cholesterol. Deriving experimental estimates for the spontaneous curvature and Gaussian modulus of curvature, we show, using a simple geometric model for ion channels, that in this way changes in the conformational equilibrium of membrane proteins can be kept small.
Collapse
Affiliation(s)
- Beate Boulgaropoulos
- Institute of Biophysics and Nanosystems Research, Austrian Academy of Sciences, Graz, Austria
| | | | | | | | | |
Collapse
|
15
|
Determining the Gaussian curvature modulus of lipid membranes in simulations. Biophys J 2012; 102:1403-10. [PMID: 22455923 DOI: 10.1016/j.bpj.2012.02.013] [Citation(s) in RCA: 145] [Impact Index Per Article: 12.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2011] [Revised: 01/25/2012] [Accepted: 02/08/2012] [Indexed: 11/20/2022] Open
Abstract
The Gaussian curvature modulus κ¯ of lipid bilayers likely contributes more than 100 kcal/mol to every cellular fission or fusion event. This huge impact on membrane remodeling energetics might be a factor that codetermines the complex lipid composition of biomembranes through tuning of κ¯. Yet, its value has been measured only for a handful of simple lipids, and no simulation has so far determined it better than a factor of two, rendering a systematic investigation of such enticing speculations impossible. Here we propose a highly accurate method to determine κ¯ in computer simulations. It relies on the interplay between curvature stress and edge tension of partially curved axisymmetric membrane disks and requires determining their closing probability. For a simplified lipid model we obtain κ¯ and its relation to the normal bending modulus κ for membranes differing both in stiffness and spontaneous lipid curvature. The elastic ratio κ¯/κ can be determined with a few percent statistical accuracy. Its value agrees with the scarce experimental data, and its change with spontaneous lipid curvature is compatible with theoretical expectations, thereby granting additional information on monolayer properties. We also show that an alternative determination of these elastic parameters based on moments of the lateral stress profile gives markedly different and unphysical values.
Collapse
|
16
|
Jackson MB. Inferring structures of kinetic intermediates in Ca(2+)-triggered exocytosis. CURRENT TOPICS IN MEMBRANES 2012; 68:185-208. [PMID: 21771500 DOI: 10.1016/b978-0-12-385891-7.00008-8] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/24/2023]
Affiliation(s)
- Meyer B Jackson
- Department of Neuroscience, University of Wisconsin School of Medicine and Public Health, Madison, Wisconsin, USA
| |
Collapse
|
17
|
Tian F, Luo Y, Zhang X. Curvature modulates the self-assembly of amphiphilic molecules. J Chem Phys 2010; 133:144701. [DOI: 10.1063/1.3499914] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
18
|
Gofman Y, Linser S, Rzeszutek A, Shental-Bechor D, Funari SS, Ben-Tal N, Willumeit R. Interaction of an antimicrobial peptide with membranes: experiments and simulations with NKCS. J Phys Chem B 2010; 114:4230-7. [PMID: 20201501 DOI: 10.1021/jp909154y] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We used Monte Carlo simulations and biophysical measurements to study the interaction of NKCS, a derivative of the antimicrobial peptide NK-2, with a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE) membrane. The simulations showed that NKCS adsorbed on the membrane surface and the dominant conformation featured two amphipathic helices connected by a hinge region. We designed two mutants in the hinge to investigate the interplay between helicity and membrane affinity. Simulations with a Leu-to-Pro substitution showed that the helicity and membrane affinity of the mutant (NKCS-[LP]) decreased. Two Ala residues were added to NKCS to produce a sequence that is compatible with a continuous amphipathic helix structure (NKCS-[AA]), and the simulations showed that the mutant adsorbed on the membrane surface with a particularly high affinity. The circular dichroism spectra of the three peptides also showed that NKCS-[LP] is the least helical and NKCS-[AA] is the most. However, the activity of the peptides, determined in terms of their antimicrobial potency and influence on the temperature of the transition of the lipid to hexagonal phase, displayed a complex behavior: NKCS-[LP] was the least potent and had the smallest influence on the transition temperature, and NKCS was the most potent and had the largest effect on the temperature.
Collapse
Affiliation(s)
- Yana Gofman
- GKSS Research Center, 21502 Geesthacht, Germany
| | | | | | | | | | | | | |
Collapse
|
19
|
Siegel DP. Fourth-order curvature energy model for the stability of bicontinuous inverted cubic phases in amphiphile-water systems. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2010; 26:8673-8683. [PMID: 20349969 DOI: 10.1021/la904838z] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/29/2023]
Abstract
The bicontinuous inverted cubic (Q(II)) phases of amphiphiles in water have many practical applications. It is necessary to understand the stability of these phases as a function of composition and ambient conditions in order to make the best use of them. Moreover, many biomembrane lipids and some biomembrane lipid extracts form Q(II) phases. The stability of Q(II) phases in a given lipid composition is closely related to the susceptibility of that composition to membrane fusion: changes in composition that stabilize Q(II) phases usually increase the rate of membrane fusion. However, the factors determining Q(II) phase stability are not fully understood. Previously, an expression was derived for the curvature free energy of Q(II) phases with respect to that of the lamellar (L(alpha)) phase using a model for the curvature energy with terms up to fourth order in curvature as formulated by Mitov. Here this model is extended to account for the effects of water content on Q(II) phase stability. It is shown that the observed L(alpha)/Q(II) phase-transition temperature, transition enthalpy, and transition kinetics are all sensitive to water content. The same observables also become sensitive to small noncurvature energy contributions to the total free-energy difference between the Q(II) and L(alpha) phases, especially the unbinding energy in the L(alpha) phase. These predictions rationalize earlier observations of Q(II) phase formation in N-monomethylated dioleoylphosphatidylethanolamine that otherwise appear to be inconsistent. The model also provides a fundamental explanation of the hysteresis typically observed in transitions between the L(alpha) and Q(II) phases. It is an accurate model of Q(II) phase stability when the ratio of the volume fraction of the lipid in the Q(II) phase unit cell is < or = 0.5.
Collapse
Affiliation(s)
- David P Siegel
- Givaudan Inc., 1199 Edison Drive, Cincinnati, Ohio 45216, USA.
| |
Collapse
|
20
|
den Otter WK. Free energies of stable and metastable pores in lipid membranes under tension. J Chem Phys 2010; 131:205101. [PMID: 19947707 DOI: 10.1063/1.3266839] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
The free energy profile of pore formation in a lipid membrane, covering the entire range from a density fluctuation in an intact bilayer to a large tension-stabilized pore, has been calculated by molecular dynamics simulations with a coarse-grained lipid model. Several fixed elongations are used to obtain the Helmholtz free energy as a function of pore size for thermodynamically stable, metastable, and unstable pores, and the system-size dependence of these elongations is discussed. A link to the Gibbs free energy at constant tension, commonly known as the Litster model, is established by a Legendre transformation. The change of genus upon pore formation is exploited to estimate the saddle-splay modulus or Gaussian curvature modulus of the membrane leaflets. Details are provided of the simulation approach, which combines the potential of mean constraint force method with a reaction coordinate based on the local lipid density.
Collapse
Affiliation(s)
- Wouter K den Otter
- Computational Biophysics, Faculty of Science and Technology, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands.
| |
Collapse
|
21
|
Minimum membrane bending energies of fusion pores. J Membr Biol 2009; 231:101-15. [PMID: 19865786 DOI: 10.1007/s00232-009-9209-x] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2009] [Accepted: 09/24/2009] [Indexed: 10/20/2022]
Abstract
Membranes fuse by forming highly curved intermediates, culminating in structures described as fusion pores. These hourglass-like figures that join two fusing membranes have high bending energies, which can be estimated using continuum elasticity models. Fusion pore bending energies depend strongly on shape, and the present study developed a method for determining the shape that minimizes bending energy. This was first applied to a fusion pore modeled as a single surface and then extended to a more realistic model treating a bilayer as two monolayers. For the two-monolayer model, fusion pores were found to have metastable states with energy minima at particular values of the pore diameter and bilayer separation. Fusion pore energies were relatively insensitive to membrane thickness but highly sensitive to spontaneous curvature and membrane asymmetry. With symmetrical bilayers and monolayer spontaneous curvatures of -0.1 nm(-1) (a typical value) separated by 6 nm (closest distance determined by repulsive hydration forces), fusion pore formation required 43-65 kT. The pore radius of approximately 2.25 nm fell within the range estimated from conductance measurements. With bilayer separation >6 nm, fusion pore formation required less energy, suggesting that protein scaffolds can promote fusion by bending membranes toward one another. With nonzero spontaneous monolayer curvature, the shape that minimized the energy change during fusion pore formation differed from the shape that minimized its energy after it formed. Thus, a nascent fusion pore will relax spontaneously to a new shape, consistent with the experimentally observed expansion of nascent fusion pores during viral fusion.
Collapse
|
22
|
Mares T, Daniel M, Perutkova S, Perne A, Dolinar G, Iglic A, Rappolt M, Kralj-Iglic V. Role of phospholipid asymmetry in the stability of inverted hexagonal mesoscopic phases. J Phys Chem B 2009; 112:16575-84. [PMID: 19367813 DOI: 10.1021/jp805715r] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The role of phospholipid asymmetry in the transition from the lamellar (L(alpha)) to the inverted hexagonal (H(II)) phase upon the temperature increase was considered. The equilibrium configuration of the system was determined by the minimum of the free energy including the contribution of the isotropic and deviatoric bending and the interstitial energy of phospholipid monolayers. The shape and local interactions of a single lipid molecule were taken into account. The minimization with respect to the configuration of the lipid layers was performed by a numerical solution of the system of the Euler-Lagrange differential equations and by the Monte Carlo simulated annealing method. At high enough temperature, the lipid molecules attain a shape exhibiting higher intrinsic mean and deviatoric curvatures, which fits better into the H(II) phase than into the L(alpha) phase. Furthermore, the orientational ordering of lipid molecules in the curvature field expressed as the deviatoric bending provides a considerable negative contribution to the free energy, which stabilizes the nonlamellar H(II) phase. The nucleation configuration for the L(alpha)-H(II) phase transition is tuned by the isotropic and deviatoric bending energies and the interstitial energy.
Collapse
Affiliation(s)
- Tomas Mares
- Laboratory of Biomechanics, Faculty of Mechanical Engineering, Czech Technical University in Prague, Technická 4, CZ-166 07, Prague 6, Czech Republic.
| | | | | | | | | | | | | | | |
Collapse
|
23
|
Abstract
The Gaussian curvature elastic energy contribution to the energy of membrane fusion intermediates has usually been neglected because the Gaussian curvature elastic modulus, kappa, was unknown. It is now possible to measure kappa for phospholipids that form bicontinuous inverted cubic (Q(II)) phases. Here, it is shown that one can estimate kappa for lipids that do not form Q(II) phases by studying the phase behavior of lipid mixtures. The method is used to estimate kappa for several lipid compositions in excess water. The values of kappa are used to compute the curvature elastic energies of stalks and catenoidal fusion pores according to recent models. The Gaussian curvature elastic contribution is positive and similar in magnitude to the bending energy contribution: it increases the total curvature energy of all the fusion intermediates by 100 units of k(B)T or more. It is important to note that this contribution makes the predicted intermediate energies compatible with observed lipid phase behavior in excess water. An order-of-magnitude fusion rate equation is used to estimate whether the predicted stalk energies are consistent with the observed rates of stalk-mediated processes in pure lipid systems. The current theory predicts a stalk energy that is slightly too large, by approximately 30 k(B)T, to rationalize the observed rates of stalk-mediated processes in phosphatidylethanolamine or N-monomethylated dioleoylphosphatidylethanolamine systems. Despite this discrepancy, the results show that models of fusion intermediate energy are accurate enough to make semiquantitative predictions about how proteins mediate biomembrane fusion. The same rate model shows that for proteins to drive biomembrane fusion at observed rates, they have to perform mediating functions corresponding to a reduction in the energy of a purely lipidic stalk by several tens of k(B)T. By binding particular peptide sequences to the monolayer surface, proteins could lower fusion intermediate energies by altering the elastic constants of the patches of lipid monolayer that form the stalk. Here, it is shown that if peptide binding changes kappa or some other combinations of local elastic constants by only tens of percents, the stalk energy and the energy of catenoidal fusion pores would decrease by tens of k(B)T relative to the pure lipid value. This is comparable to the required mediating effect. The curvature energies of stalks and catenoidal fusion pores have almost the same dependence on monolayer elastic constants as the curvature energies of the rhombohedral and Q(II) phases; respectively. The effects of isolated fusion-relevant peptides on the energies of these intermediates can be determined by studying the effects of the peptides on the stability of rhombohedral and Q(II) phases.
Collapse
|
24
|
Rugonyi S, Biswas SC, Hall SB. The biophysical function of pulmonary surfactant. Respir Physiol Neurobiol 2008; 163:244-55. [PMID: 18632313 DOI: 10.1016/j.resp.2008.05.018] [Citation(s) in RCA: 76] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2008] [Revised: 05/23/2008] [Accepted: 05/24/2008] [Indexed: 01/12/2023]
Abstract
Pulmonary surfactant lowers surface tension in the lungs. Physiological studies indicate two key aspects of this function: that the surfactant film forms rapidly; and that when compressed by the shrinking alveolar area during exhalation, the film reduces surface tension to very low values. These observations suggest that surfactant vesicles adsorb quickly, and that during compression, the adsorbed film resists the tendency to collapse from the interface to form a 3D bulk phase. Available evidence suggests that adsorption occurs by way of a rate-limiting structure that bridges the gap between the vesicle and the interface, and that the adsorbed film avoids collapse by undergoing a process of solidification. Current models, although incomplete, suggest mechanisms that would partially explain both rapid adsorption and resistance to collapse as well as how different constituents of pulmonary surfactant might affect its behavior.
Collapse
Affiliation(s)
- Sandra Rugonyi
- Division of Biomedical Engineering, Oregon Health & Science University, Portland, OR 97239-3098, USA.
| | | | | |
Collapse
|
25
|
Hickel A, Danner-Pongratz S, Amenitsch H, Degovics G, Rappolt M, Lohner K, Pabst G. Influence of antimicrobial peptides on the formation of nonlamellar lipid mesophases. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2008; 1778:2325-33. [PMID: 18582435 DOI: 10.1016/j.bbamem.2008.05.014] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/14/2008] [Revised: 05/30/2008] [Accepted: 05/30/2008] [Indexed: 10/22/2022]
Abstract
We have studied the influence of four antimicrobial peptides of different secondary and ternary structure--melittin (Mel), protegrin-1 (PG-1), peptidyl-glycylleucine-carboxyamide (PGLa), and gramicidin S (GS)--on the lamellar-to-nonlamellar transition of palmitoyloleoyl phosphatidylethanolamine (POPE) applying differential scanning calorimetry and small-angle X-ray diffraction. None of the peptides studied led to the formation of an inverted hexagonal phase observed for pure POPE at high temperatures. Instead either cubic or lamellar phases were stabilized to different degrees. GS was most effective in inducing a cubic phase, whereas Mel fully stabilized the lamellar phase. The behavior of POPE in the presence of PG-1 and PGLa was intermediate to GS and Mel. In addition to the known role of membrane elasticity we propose two mechanisms, which cause stabilization of the lamellar phase: electrostatic repulsion and lipid/peptide pore formation. Both mechanisms prevent transmembrane contact required to form either an inverted hexagonal phase or fusion pores, as precursors of the cubic phase.
Collapse
Affiliation(s)
- Andrea Hickel
- Institute of Biophysics and Nanosystems Research, Austrian Academy of Sciences, Schmiedlstrasse 6, A 8042 Graz, Austria
| | | | | | | | | | | | | |
Collapse
|
26
|
Momot KI, Takegoshi K, Kuchel PW, Larkin TJ. Inhomogeneous NMR line shape as a probe of microscopic organization of bicontinuous cubic phases. J Phys Chem B 2008; 112:6636-45. [PMID: 18457445 DOI: 10.1021/jp8006415] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
NMR line shapes of the lipid and aqueous species in bicontinuous cubic phase (BCP) samples prepared by centrifugation are inhomogeneously broadened. The broadening of the lipid peaks is removed by magic-angle spinning (MAS). In this work, we studied the mechanism of this broadening using (1)H and (13)C NMR spectroscopy of a myverol/water BCP. It is demonstrated that the inhomogeneity possesses an intrinsic contribution that is independent of instrumental or setup factors and can be attributed to the microscopic organization of the BCP bilayer. A mechanism of the inhomogeneous broadening is proposed, which involves a spatially nonuniform diamagnetically induced magnetic field determined by the mesoscopic structure and the diamagnetic susceptibilities of the two BCP domains. The proposed mechanism does not require that molecular reorientation of the lipid be slow for the inhomogeneous broadening to survive. We discuss how this inhomogeneous broadening can be employed as a probe of compositional uniformity and microscopic organization of BCP samples.
Collapse
Affiliation(s)
- Konstantin I Momot
- School of Physical and Chemical Sciences, Queensland University of Technology, GPO Box 2434, Brisbane, Qld 4001, Australia.
| | | | | | | |
Collapse
|
27
|
Influence of the lamellar phase unbinding energy on the relative stability of lamellar and inverted cubic phases. Biophys J 2008; 94:3987-95. [PMID: 18234828 DOI: 10.1529/biophysj.107.118034] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Based on curvature energy considerations, nonbilayer phase-forming phospholipids in excess water should form stable bicontinuous inverted cubic (Q(II)) phases at temperatures between the lamellar (L(alpha)) and inverted hexagonal (H(II)) phase regions. However, the phosphatidylethanolamines (PEs), which are a common class of biomembrane phospholipids, typically display direct L(alpha)/H(II) phase transitions and may form intermediate Q(II) phases only after the temperature is cycled repeatedly across the L(alpha)/H(II) phase transition temperature, T(H), or when the H(II) phases are cooled from T > T(H). This raises the question of whether models of inverted phase stability, which are based on curvature energy alone, accurately predict the relative free energy of these phases. Here we demonstrate the important role of a noncurvature energy contribution, the unbinding energy of the L(alpha) phase bilayers, g(u), that serves to stabilize the L(alpha) phase relative to the nonlamellar phases. The planar L(alpha) phase bilayers must separate for a Q(II) phase to form and it turns out that the work of their unbinding can be larger than the curvature energy reduction on formation of Q(II) phase from L(alpha) at temperatures near the L(alpha)/Q(II) transition temperature (T(Q)). Using g(u) and elastic constant values typical of unsaturated PEs, we show that g(u) is sufficient to make T(Q) > T(H) for the latter lipids. Such systems would display direct L(alpha) --> H(II) transitions, and a Q(II) phase might only form as a metastable phase upon cooling of the H(II) phase. The g(u) values for methylated PEs and PE/phosphatidylcholine mixtures are significantly smaller than those for PEs and increase T(Q) by only a few degrees, consistent with observations of these systems. This influence of g(u) also rationalizes the effect of some aqueous solutes to increase the rate of Q(II) formation during temperature cycling of lipid dispersions. Finally, the results are relevant to protocols for determining the Gaussian curvature modulus, which substantially affects the energy of intermediates in membrane fusion and fission. Recently, two such methods were proposed based on measuring T(Q) and on measuring Q(II) phase unit cell dimensions, respectively. In view of the effect of g(u) on T(Q) that we describe here, the latter method, which does not depend on the value of g(u), is preferable.
Collapse
|
28
|
Smith KA, Jasnow D, Balazs AC. Designing synthetic vesicles that engulf nanoscopic particles. J Chem Phys 2007; 127:084703. [PMID: 17764280 DOI: 10.1063/1.2766953] [Citation(s) in RCA: 126] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We examine the interaction of a lipid bilayer membrane with a spherical particle in solution using dissipative particle dynamics, with the aim of controlling the passage of foreign objects into and out of vesicles. Parameters are chosen such that there is a favorable adhesive interaction between the membrane and the particle. Under these conditions, the membrane wraps the particle in a process resembling phagocytosis in biological cells. We find that, for a homogeneous membrane with a uniform attraction to the particle, the membrane is unable to fully wrap the particle when the adhesion strength is below a certain value. This is observed even in the limit of zero membrane tension. When the adhesion strength is increased above the threshold value, the membrane fully wraps the particle. However, the wrapped particle remains tethered to the larger membrane. We next consider an adhesive domain, or raft, in an otherwise nonadhesive membrane. We find that, when the particle is wrapped by the raft, the line tension at the raft interface promotes fission, allowing the wrapped particle to detach from the larger membrane. This mechanism could be used to allow particles to cross a vesicle membrane.
Collapse
Affiliation(s)
- Kurt A Smith
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, USA
| | | | | |
Collapse
|