1
|
Richardson JD, Van Lehn RC. Free Energy Analysis of Peptide-Induced Pore Formation in Lipid Membranes by Bridging Atomistic and Coarse-Grained Simulations. J Phys Chem B 2024; 128:8737-8752. [PMID: 39207202 DOI: 10.1021/acs.jpcb.4c03276] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/04/2024]
Abstract
Antimicrobial peptides (AMPs) are attractive materials for combating the antimicrobial resistance crisis because they can kill target microbes by directly disrupting cell membranes. Although thousands of AMPs have been discovered, their molecular mechanisms of action are still poorly understood. One broad mechanism for membrane disruption is the formation of membrane-spanning hydrophilic pores which can be stabilized by AMPs. In this study, we use molecular dynamics simulations to investigate the thermodynamics of pore formation in model single-component lipid membranes in the presence of one of three AMPs: aurein 1.2, melittin and magainin 2. To overcome the general challenge of modeling long time scale membrane-related behaviors, including AMP binding, clustering, and pore formation, we develop a generalizable methodology for sampling AMP-induced pore formation. This approach involves the long equilibration of peptides around a pore created with a nucleation collective variable by performing coarse-grained simulations, then backmapping equilibrated AMP-membrane configurations to all-atom resolution. We then perform all-atom simulations to resolve free energy profiles for pore formation while accurately modeling the interplay of lipid-peptide-solvent interactions that dictate pore formation free energies. Using this approach, we quantify free energy barriers for pore formation without direct biases on peptides or whole lipids, allowing us to investigate mechanisms of pore formation for these 3 AMPs that are a consequence of unbiased peptide diffusion and clustering. Further analysis of simulation trajectories then relates variations in pore lining by AMPs, AMP-induced lipid disruptions, and salt bridges between AMPs to the observed pore formation free energies and corresponding mechanisms. This methodology and mechanistic analysis have the potential to generalize beyond the AMPs in this study to improve our understanding of pore formation by AMPs and related antimicrobial materials.
Collapse
Affiliation(s)
- Joshua D Richardson
- Department of Chemical and Biological Engineering, University of Wisconsin─Madison, Madison, Wisconsin 53706, United States
| | - Reid C Van Lehn
- Department of Chemical and Biological Engineering, University of Wisconsin─Madison, Madison, Wisconsin 53706, United States
- Department of Chemistry, University of Wisconsin─Madison, Madison, Wisconsin 53706, United States
| |
Collapse
|
2
|
Klug J, Berberián MV, López Martí JM, Mayorga LS, Del Pópolo MG. Membrane Binding Strength vs Pore Formation Cost─What Drives the Membrane Permeation of Nanoparticles Coated with Cell-Penetrating Peptides? J Phys Chem B 2024; 128:937-948. [PMID: 38232319 DOI: 10.1021/acs.jpcb.3c06461] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2024]
Abstract
Cell-penetrating peptides (CPPs) enable the transport of nanoparticles through cell membranes. Using molecular simulations, we conduct an in-depth investigation into the thermodynamic forces governing the passive translocation of CPP-coated nanoparticles across lipid bilayers, contrasting their behavior with that of bare particles to dissect the contribution of the peptides. Our analysis unveils a distinctive two-stage translocation mechanism, where the adsorption energy of the particles overcomes the cost of forming a hydrophilic transmembrane pore. Proper evaluation of the translocation mechanisms is only possible when using two reaction coordinates, in particular, one that explicitly includes the density of the lipids on the binding site of the particle. An analysis of adsorption and activation free energies in terms of a simple kinetic model provides a clearer understanding of the CPP effect. Experimental validation using nonendocytic cells confirms the superior membrane permeation of CPP-coated particles. Our findings have implications for the rational design of more efficient cell-permeating particles.
Collapse
Affiliation(s)
- Joaquín Klug
- Department of Life Sciences, Faculty of Sciences, Atlantic Technological University, ATU Sligo, Ash Lane, Sligo F91 YW50, Ireland
| | - María Victoria Berberián
- Instituto de Histología y Embriología de Mendoza (IHEM) CONICET/UNCUYO. Centro Universitario, Mendoza M5502JMA, Argentina
- Instituto Interdisciplinario de Ciencias Básicas (ICB)─CONICET/UNCUYO, Padre Jorge Contreras 1300, Mendoza CP5500, Argentina
- Facultad de Ciencias Exactas y Naturales, Universidad Nacional de Cuyo, Padre Jorge Contreras 1300, Mendoza CP5500, Argentina
| | - Jesús M López Martí
- Facultad de Ciencias Exactas y Naturales, Universidad Nacional de Cuyo, Padre Jorge Contreras 1300, Mendoza CP5500, Argentina
| | - Luis S Mayorga
- Instituto de Histología y Embriología de Mendoza (IHEM) CONICET/UNCUYO. Centro Universitario, Mendoza M5502JMA, Argentina
- Facultad de Ciencias Exactas y Naturales, Universidad Nacional de Cuyo, Padre Jorge Contreras 1300, Mendoza CP5500, Argentina
| | - Mario G Del Pópolo
- Instituto Interdisciplinario de Ciencias Básicas (ICB)─CONICET/UNCUYO, Padre Jorge Contreras 1300, Mendoza CP5500, Argentina
- Facultad de Ciencias Exactas y Naturales, Universidad Nacional de Cuyo, Padre Jorge Contreras 1300, Mendoza CP5500, Argentina
| |
Collapse
|
3
|
Deng Z, You X, Lin Z, Dong X, Yuan B, Yang K. Membrane-Active Peptides Attack Cell Membranes in a Lipid-Regulated Curvature-Generating Mode. J Phys Chem Lett 2023:6422-6430. [PMID: 37432779 DOI: 10.1021/acs.jpclett.3c01115] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/13/2023]
Abstract
Membrane-active peptides (MAPs) exhibit great potential in biomedical applications due to their unique ability to overcome the cell membrane barrier. However, the interactions between MAPs and membranes are complex, and little is known about the possibility of MAP action being specific to certain types of membranes. In this study, a combination of molecular dynamics simulations and theoretical analysis was utilized to investigate the interactions between typical MAPs and realistic cell membrane systems. Remarkably, the simulations revealed that MAPs can attack membranes by generating and sensing positive mean curvature, which is dependent on lipid composition. Furthermore, theoretical calculations demonstrated that this lipid-regulated curvature-based membrane attack mechanism is an integrated result of multiple effects, including peptide-induced membrane wedge and softening effects, the lipid shape effect, the area-difference elastic effect, and the boundary edge effect of formed peptide-lipid nanodomains. This study enhances our comprehension of MAP-membrane interactions and highlights the potential for developing membrane-specific MAP-based agents.
Collapse
Affiliation(s)
- Zhixiong Deng
- Center for Soft Condensed Matter Physics and Interdisciplinary Research & School of Physical Science and Technology, Soochow University, Suzhou 215006, Jiangsu, China
| | - Xin You
- Center for Soft Condensed Matter Physics and Interdisciplinary Research & School of Physical Science and Technology, Soochow University, Suzhou 215006, Jiangsu, China
| | - Zhao Lin
- Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
| | - Xuewei Dong
- Center for Soft Condensed Matter Physics and Interdisciplinary Research & School of Physical Science and Technology, Soochow University, Suzhou 215006, Jiangsu, China
| | - Bing Yuan
- Songshan Lake Materials Laboratory, Dongguan 523808 Guangdong, China
| | - Kai Yang
- Center for Soft Condensed Matter Physics and Interdisciplinary Research & School of Physical Science and Technology, Soochow University, Suzhou 215006, Jiangsu, China
- Songshan Lake Materials Laboratory, Dongguan 523808 Guangdong, China
| |
Collapse
|
4
|
Mou Q, Xu M, Deng J, Hu N, Yang J. Studying the roles of salt ions in the pore initiation and closure stages in the biomembrane electroporation. APL Bioeng 2023; 7:026103. [PMID: 37180734 PMCID: PMC10168715 DOI: 10.1063/5.0147104] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2023] [Accepted: 04/24/2023] [Indexed: 05/16/2023] Open
Abstract
Electroporation shows great potential in biology and biomedical applications. However, there is still a lack of reliable protocol for cell electroporation to achieve a high perforation efficiency due to the unclear influence mechanism of various factors, especially the salt ions in buffer solution. The tiny membrane structure of a cell and the electroporation scale make it difficult to monitor the electroporation process. In this study, we used both molecular dynamics (MD) simulation and experimental methods to explore the influence of salt ions on the electroporation process. Giant unilamellar vesicles (GUVs) were constructed as the model, and sodium chloride (NaCl) was selected as the representative salt ion in this study. The results show that the electroporation process follows lag-burst kinetics, where the lag period first appears after applying the electric field, followed by a rapid pore expansion. For the first time, we find that the salt ion plays opposite roles in different stages of the electroporation process. The accumulation of salt ions near the membrane surface provides an extra potential to promote the pore initiation, while the charge screening effect of the ions within the pore increases the line tension of the pore to induce the instability of the pore and lead to the closure. The GUV electroporation experiments obtain qualitatively consistent results with MD simulations. This work can provide guidance for the selection of parameters for cell electroporation process.
Collapse
Affiliation(s)
- Qiongyao Mou
- Key Laboratory of Biorheological Science and Technology, Ministry of Education and Bioengineering College, Chongqing University, Chongqing 400044, China
| | - Mengli Xu
- Key Laboratory of Biorheological Science and Technology, Ministry of Education and Bioengineering College, Chongqing University, Chongqing 400044, China
| | - Jinan Deng
- Key Laboratory of Biorheological Science and Technology, Ministry of Education and Bioengineering College, Chongqing University, Chongqing 400044, China
| | - Ning Hu
- Key Laboratory of Biorheological Science and Technology, Ministry of Education and Bioengineering College, Chongqing University, Chongqing 400044, China
| | - Jun Yang
- Key Laboratory of Biorheological Science and Technology, Ministry of Education and Bioengineering College, Chongqing University, Chongqing 400044, China
| |
Collapse
|
5
|
Kasparyan G, Hub JS. Equivalence of Charge Imbalance and External Electric Fields during Free Energy Calculations of Membrane Electroporation. J Chem Theory Comput 2023; 19:2676-2683. [PMID: 37052575 DOI: 10.1021/acs.jctc.3c00065] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/14/2023]
Abstract
Electric fields across lipid membranes play important roles in physiology, medicine, and biotechnology, rationalizing the wide interest in modeling transmembrane potentials in molecular dynamics simulations. Transmembrane potentials have been implemented with external electric fields or by imposing charge imbalance between the two water compartments of a stacked double-membrane system. We compare the two methods in the context of membrane electroporation, which involves a large change of membrane structure and capacitance. We show that, given that Ewald electrostatics are defined with tinfoil boundary conditions, the two methods lead to (i) identical potentials of mean force (PMFs) of pore formation and expansion at various potentials, demonstrating that the two methods impose equivalent driving forces for large-scale transitions at membranes, and (ii) to identical polarization of water within thin water wires or open pores, suggesting that the two methods furthermore impose equivalent local electric fields. Without tinfoil boundary conditions, effects from external fields on pore formation are spuriously suppressed or even removed. Together, our study shows that both methods, external fields and charge imbalance, are well suitable for studying large-scale transitions of lipid membranes that involve changes of membrane capacitance. However, using charge imbalance is technically more challenging for maintaining a constant transmembrane potential since it requires updating of the charge imbalance as the membrane capacitance changes.
Collapse
Affiliation(s)
- Gari Kasparyan
- Theoretical Physics and Center for Biophysics, Saarland University, 66123 Saarbrücken, Germany
| | - Jochen S Hub
- Theoretical Physics and Center for Biophysics, Saarland University, 66123 Saarbrücken, Germany
| |
Collapse
|
6
|
Ouyang J, Sheng Y, Wang W. Recent Advances of Studies on Cell-Penetrating Peptides Based on Molecular Dynamics Simulations. Cells 2022; 11:cells11244016. [PMID: 36552778 PMCID: PMC9776715 DOI: 10.3390/cells11244016] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2022] [Revised: 12/09/2022] [Accepted: 12/10/2022] [Indexed: 12/14/2022] Open
Abstract
With the ability to transport cargo molecules across cell membranes with low toxicity, cell-penetrating peptides (CPPs) have become promising candidates for next generation peptide-based drug delivery vectors. Over the past three decades since the first CPP was discovered, a great deal of work has been done on the cellular uptake mechanisms and the applications for the delivery of therapeutic molecules, and significant advances have been made. But so far, we still do not have a precise and unified understanding of the structure-activity relationship of the CPPs. Molecular dynamics (MD) simulations provide a method to reveal peptide-membrane interactions at the atomistic level and have become an effective complement to experiments. In this paper, we review the progress of the MD simulations on CPP-membrane interactions, including the computational methods and technical improvements in the MD simulations, the research achievements in the CPP internalization mechanism, CPP decoration and coupling, and the peptide-induced membrane reactions during the penetration process, as well as the comparison of simulated and experimental results.
Collapse
Affiliation(s)
- Jun Ouyang
- School of Public Courses, Bengbu Medical College, Bengbu 233030, China
- Collaborative Innovation Center of Advanced Microstructures, National Laboratory of Solid State Microstructure, Department of Physics, Nanjing University, Nanjing 210093, China
| | - Yuebiao Sheng
- Collaborative Innovation Center of Advanced Microstructures, National Laboratory of Solid State Microstructure, Department of Physics, Nanjing University, Nanjing 210093, China
- High Performance Computing Center, Nanjing University, Nanjing 210093, China
- Correspondence: (Y.S.); (W.W.)
| | - Wei Wang
- Collaborative Innovation Center of Advanced Microstructures, National Laboratory of Solid State Microstructure, Department of Physics, Nanjing University, Nanjing 210093, China
- Correspondence: (Y.S.); (W.W.)
| |
Collapse
|
7
|
Alvarez AB, Caruso B, Petersen SB, Rodríguez PEA, Fidelio GD. Melittin-solid phospholipid mixed films trigger amyloid-like nano-fibril arrangements at air-water interface. BIOCHIMICA ET BIOPHYSICA ACTA. BIOMEMBRANES 2022; 1864:184048. [PMID: 36115495 DOI: 10.1016/j.bbamem.2022.184048] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/14/2022] [Revised: 08/22/2022] [Accepted: 09/06/2022] [Indexed: 06/15/2023]
Abstract
We used the Langmuir monolayers technique to study the surface properties of melittin toxin mixed with either liquid-condensed DSPC or liquid-expanded POPC phospholipids. Pure melittin peptide forms stable insoluble monolayers at the air-water interface without interacting with Thioflavin T (Th-T), a sensitive probe to detect protein amyloid formation. When melittin peptide is mixed with DSPC lipid at 50 % of peptide area proportion at the surface, we observed the formation of fibril-like structures detected by Brewster angle microscopy (BAM), but they were not observable with POPC. The nano-structures in the melittin-DSPC mixtures became Th-T positive labeling when the arrangement was observed with fluorescence microscopy. In this condition, Th-T undergoes an unexpected shift in the typical emission wavelength of this amyloid marker when a 2D fluorescence analysis is conducted. Even when reflectivity analysis of BAM imaging evidenced that these structures would correspond to the DSPC lipid component of the mixture, the interpretation of ATR-FTIR and Th-T data suggested that both components were involved in a new lipid-peptide rearrangement. These nano-fibril arrangements were also evidenced by scanning electron and atomic force microscopy when the films were transferred to a mica support. The fibril formation was not detected when melittin was mixed with the liquid-expanded POPC lipid. We postulated that DSPC lipids can dynamically trigger the process of amyloid-like nano-arrangement formation at the interface. This process is favored by the relative peptide content, the quality of the interfacial environment, and the physical state of the lipid at the surface.
Collapse
Affiliation(s)
- Alain Bolaño Alvarez
- Departamento de Química Biológica Ranwel Caputto, Facultad de Ciencias Químicas, Universidad Nacional de Córdoba, Argentina; Centro de Investigaciones en Química Biológica de Córdoba, CIQUIBIC, CONICET, Universidad Nacional de Córdoba, Argentina
| | - Benjamín Caruso
- Cátedra de Química Biológica, Departamento de Química, Facultad de Ciencias Exactas, Físicas y Naturales, Universidad Nacional de Córdoba, Argentina; Instituto de Investigaciones Biológicas y Tecnológicas (IIBYT), CONICET, Universidad Nacional de Córdoba. Córdoba, Argentina
| | | | | | - Gerardo D Fidelio
- Departamento de Química Biológica Ranwel Caputto, Facultad de Ciencias Químicas, Universidad Nacional de Córdoba, Argentina; Centro de Investigaciones en Química Biológica de Córdoba, CIQUIBIC, CONICET, Universidad Nacional de Córdoba, Argentina.
| |
Collapse
|
8
|
Guo F, Wang J, Zhou J, Qian K, Qu H, Liu P, Zhai S. All-atom molecular dynamics simulations of the combined effects of different phospholipids and cholesterol content on electroporation. RSC Adv 2022; 12:24491-24500. [PMID: 36128384 PMCID: PMC9425445 DOI: 10.1039/d2ra03895a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2022] [Accepted: 07/27/2022] [Indexed: 11/21/2022] Open
Abstract
The electroporation mechanism could be related to the composition of the plasma membrane, and the combined effect of different phospholipid molecules and cholesterol content on electroporation has rarely been studied nor conclusions drawn. In this paper, we applied all-atom molecular dynamics (MD) simulations to study the effects of phospholipids and cholesterol content on bilayer membrane electroporation. The palmitoyloleoylphosphatidylcholine (POPC) model, palmitoyloleoylphosphatidylethanolamine (POPE) model, and a 1 : 1 mixed model of POPC and POPE called PEPC, were the three basic models used. An electric field of 0.45 V nm-1 was applied to nine models, which were the three basic models, each with three different cholesterol content values of 0%, 24%, and 40%. The interfacial water molecules moved under the electric field and, once the first water bridge formed, the rest of the water molecules would dramatically flood into the membrane. The simulation showed that a rapid rise in the Z-component of the average dipole moment of the interfacial water molecules (Z-DM) indicated the occurrence of electroporation, and the same increment of Z-DM represented a similar change in the size of the water bridge. With the same cholesterol content, the formation of the first water bridge was the most rapid in the POPC model, regarding the average electroporation time (t ep), and the average t ep of the PEPC model was close to that of the POPE model. We speculate that the differences in membrane thickness and initial number of hydrogen bonds of the interfacial water molecules affect the average t ep for different membrane compositions. Our results reveal the influence of membrane composition on the electroporation mechanism at the molecular level.
Collapse
Affiliation(s)
- Fei Guo
- Institute of Ecological Safety, Chongqing University of Posts and Telecommunications Chongqing 400065 China
| | - Ji Wang
- Institute of Ecological Safety, Chongqing University of Posts and Telecommunications Chongqing 400065 China
| | - Jiong Zhou
- Institute of Ecological Safety, Chongqing University of Posts and Telecommunications Chongqing 400065 China
| | - Kun Qian
- Institute of Ecological Safety, Chongqing University of Posts and Telecommunications Chongqing 400065 China
| | - Hongchun Qu
- Institute of Ecological Safety, Chongqing University of Posts and Telecommunications Chongqing 400065 China
| | - Ping Liu
- Institute of Ecological Safety, Chongqing University of Posts and Telecommunications Chongqing 400065 China
| | - Shidong Zhai
- Institute of Ecological Safety, Chongqing University of Posts and Telecommunications Chongqing 400065 China
| |
Collapse
|
9
|
Saha S, Ratrey P, Mishra A. Association of Lasioglossin-III Antimicrobial Peptide with Model Lipid Bilayers. Biophysics (Nagoya-shi) 2022. [DOI: 10.1134/s0006350922020178] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022] Open
|
10
|
Rems L, Tang X, Zhao F, Pérez-Conesa S, Testa I, Delemotte L. Identification of electroporation sites in the complex lipid organization of the plasma membrane. eLife 2022; 11:e74773. [PMID: 35195069 PMCID: PMC8912918 DOI: 10.7554/elife.74773] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2021] [Accepted: 02/22/2022] [Indexed: 11/13/2022] Open
Abstract
The plasma membrane of a biological cell is a complex assembly of lipids and membrane proteins, which tightly regulate transmembrane transport. When a cell is exposed to strong electric field, the membrane integrity becomes transiently disrupted by formation of transmembrane pores. This phenomenon termed electroporation is already utilized in many rapidly developing applications in medicine including gene therapy, cancer treatment, and treatment of cardiac arrhythmias. However, the molecular mechanisms of electroporation are not yet sufficiently well understood; in particular, it is unclear where exactly pores form in the complex organization of the plasma membrane. In this study, we combine coarse-grained molecular dynamics simulations, machine learning methods, and Bayesian survival analysis to identify how formation of pores depends on the local lipid organization. We show that pores do not form homogeneously across the membrane, but colocalize with domains that have specific features, the most important being high density of polyunsaturated lipids. We further show that knowing the lipid organization is sufficient to reliably predict poration sites with machine learning. Additionally, by analysing poration kinetics with Bayesian survival analysis we show that poration does not depend solely on local lipid arrangement, but also on membrane mechanical properties and the polarity of the electric field. Finally, we discuss how the combination of atomistic and coarse-grained molecular dynamics simulations, machine learning methods, and Bayesian survival analysis can guide the design of future experiments and help us to develop an accurate description of plasma membrane electroporation on the whole-cell level. Achieving this will allow us to shift the optimization of electroporation applications from blind trial-and-error approaches to mechanistic-driven design.
Collapse
Affiliation(s)
- Lea Rems
- KTH Royal Institute of Technology, Dept. Applied Physics, Science for Life LaboratorySolnaSweden
- University of Ljubljana, Faculty of Electrical EngineeringLjubljanaSlovenia
| | - Xinru Tang
- KTH Royal Institute of Technology, Dept. Applied Physics, Science for Life LaboratorySolnaSweden
- University of Chinese Academy of SciencesBeijingChina
| | - Fangwei Zhao
- KTH Royal Institute of Technology, Dept. Applied Physics, Science for Life LaboratorySolnaSweden
- University of Chinese Academy of SciencesBeijingChina
| | - Sergio Pérez-Conesa
- KTH Royal Institute of Technology, Dept. Applied Physics, Science for Life LaboratorySolnaSweden
| | - Ilaria Testa
- KTH Royal Institute of Technology, Dept. Applied Physics, Science for Life LaboratorySolnaSweden
| | - Lucie Delemotte
- KTH Royal Institute of Technology, Dept. Applied Physics, Science for Life LaboratorySolnaSweden
| |
Collapse
|
11
|
Róg T, Girych M, Bunker A. Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design. Pharmaceuticals (Basel) 2021; 14:1062. [PMID: 34681286 PMCID: PMC8537670 DOI: 10.3390/ph14101062] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2021] [Revised: 10/14/2021] [Accepted: 10/15/2021] [Indexed: 11/17/2022] Open
Abstract
We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.
Collapse
Affiliation(s)
- Tomasz Róg
- Department of Physics, University of Helsinki, 00014 Helsinki, Finland;
| | - Mykhailo Girych
- Department of Physics, University of Helsinki, 00014 Helsinki, Finland;
| | - Alex Bunker
- Drug Research Program, Division of Pharmaceutical Biosciences, Faculty of Pharmacy, University of Helsinki, 00014 Helsinki, Finland;
| |
Collapse
|
12
|
Hamal P, Subasinghege Don V, Nguyenhuu H, Ranasinghe JC, Nauman JA, McCarley RL, Kumar R, Haber LH. Influence of Temperature on Molecular Adsorption and Transport at Liposome Surfaces Studied by Molecular Dynamics Simulations and Second Harmonic Generation Spectroscopy. J Phys Chem B 2021; 125:10506-10513. [PMID: 34495664 PMCID: PMC8474114 DOI: 10.1021/acs.jpcb.1c04263] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
A fundamental understanding of the kinetics and thermodynamics of chemical interactions at the phospholipid bilayer interface is crucial for developing potential drug-delivery applications. Here we use molecular dynamics (MD) simulations and surface-sensitive second harmonic generation (SHG) spectroscopy to study the molecular adsorption and transport of a small organic cation, malachite green (MG), at the surface of 1,2-dioleoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (DOPG) liposomes in water at different temperatures. The temperature-dependent adsorption isotherms, obtained by SHG measurements, provide information on adsorbate concentration, free energy of adsorption, and associated changes in enthalpy and entropy, showing that the adsorption process is exothermic, resulting in increased overall entropy. Additionally, the molecular transport kinetics are found to be more rapid under higher temperatures. Corresponding MD simulations are used to calculate the free energy profiles of the adsorption and the molecular orientation distributions of MG at different temperatures, showing excellent agreement with the experimental results.
Collapse
Affiliation(s)
- Prakash Hamal
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803-1804, United States
| | - Visal Subasinghege Don
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803-1804, United States
| | - Huy Nguyenhuu
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803-1804, United States
| | - Jeewan C Ranasinghe
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803-1804, United States
| | - Julia A Nauman
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803-1804, United States
| | - Robin L McCarley
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803-1804, United States
| | - Revati Kumar
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803-1804, United States
| | - Louis H Haber
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803-1804, United States
| |
Collapse
|
13
|
Patel SJ, Van Lehn RC. Analysis of Charged Peptide Loop-Flipping across a Lipid Bilayer Using the String Method with Swarms of Trajectories. J Phys Chem B 2021; 125:5862-5873. [PMID: 34033491 DOI: 10.1021/acs.jpcb.1c02810] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
The hydrophobic core of the lipid bilayer is conventionally considered a barrier to the translocation of charged species such that the translocation of even single ions occurs on long time scales. In contrast, experiments have revealed that some materials, including peptides, proteins, and nanoparticles, can translocate multiple charged moieties across the bilayer on experimentally relevant time scales. Understanding the molecular mechanisms underlying this behavior is challenging because resolving corresponding free-energy landscapes with molecular simulation techniques is computationally expensive. To address this challenge, we use atomistic molecular dynamics simulations with the swarms-of-trajectories (SOT) string method to analyze charge translocation pathways across single-component lipid bilayers as a function of multiple collective variables. We first demonstrate that the SOT string method can reproduce the free-energy barrier for the translocation of a charged lysine amino acid analogue in good agreement with the literature. We then obtain minimum free-energy pathways for the translocation, or flipping, of charged peptide loops across the lipid bilayer by utilizing trajectories from prior temperature-accelerated molecular dynamics (TAMD) simulations as initial configurations. The corresponding potential of mean force calculations reveal that the protonation of a central membrane-exposed aspartate residue substantially reduces the free-energy barrier for flipping charged loops by modulating the water content of the bilayer. These results provide new insight into the thermodynamics underlying loop-flipping processes and highlight how the combination of TAMD and the SOT string method can be used to understand complex charge translocation mechanisms.
Collapse
Affiliation(s)
- Samarthaben J Patel
- Department of Chemical and Biological Engineering, University of Wisconsin-Madison, Madison, Wisconsin 53706, United States
| | - Reid C Van Lehn
- Department of Chemical and Biological Engineering, University of Wisconsin-Madison, Madison, Wisconsin 53706, United States
| |
Collapse
|
14
|
Water Pores in Planar Lipid Bilayers at Fast and Slow Rise of Transmembrane Voltage. MEMBRANES 2021; 11:membranes11040263. [PMID: 33916447 PMCID: PMC8067013 DOI: 10.3390/membranes11040263] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/07/2021] [Revised: 03/28/2021] [Accepted: 03/29/2021] [Indexed: 12/12/2022]
Abstract
Basic understanding of the barrier properties of biological membranes can be obtained by studying model systems, such as planar lipid bilayers. Here, we study water pores in planar lipid bilayers in the presence of transmembrane voltage. Planar lipid bilayers were exposed to fast and slow linearly increasing voltage and current signals. We measured the capacitance, breakdown voltage, and rupture time of planar lipid bilayers composed of 1-pamitoyl 2-oleoyl phosphatidylcholine (POPC), 1-pamitoyl 2-oleoyl phosphatidylserine (POPS), and a mixture of both lipids in a 1:1 ratio. Based on the measurements, we evaluated the change in the capacitance of the planar lipid bilayer corresponding to water pores, the radius of water pores at membrane rupture, and the fraction of the area of the planar lipid bilayer occupied by water pores.planar lipid bilayer capacitance, which corresponds to water pores, water pore radius at the membrane rupture, and a fraction of the planar lipid bilayer area occupied by water pores. The estimated pore radii determining the rupture of the planar lipid bilayer upon fast build-up of transmembrane voltage are 0.101 nm, 0.110 nm, and 0.106 nm for membranes composed of POPC, POPS, and POPC:POPS, respectively. The fraction of the surface occupied by water pores at the moment of rupture of the planar lipid bilayer The fraction of an area that is occupied by water pores at the moment of planar lipid bilayer rupture is in the range of 0.1–1.8%.
Collapse
|
15
|
Hub JS. Joint Reaction Coordinate for Computing the Free-Energy Landscape of Pore Nucleation and Pore Expansion in Lipid Membranes. J Chem Theory Comput 2021; 17:1229-1239. [PMID: 33427469 DOI: 10.1021/acs.jctc.0c01134] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Topological transitions of membranes, such as pore formation or membrane fusion, play key roles in biology, biotechnology, and in medical applications. Calculating the related free-energy landscapes has been complicated by the fact that such processes involve a sequence of transitions along highly distinct directions in conformational space, making it difficult to define good reaction coordinates (RCs) for the overall process. In this study, a new RC capable of driving both pore nucleation and pore expansion in lipid membranes is presented. The potential of mean force (PMF) along the RC computed with molecular dynamics simulations provides a comprehensive view on the free-energy landscape of pore formation, including a barrier for pore nucleation; the size, free energy, and metastability of the open pore; and the energetic cost for further pore expansion against the line tension of the pore rim. The RC is illustrated by quantifying the effects of (i) simulation system size and (ii) the addition of dimethyl sulfoxide on the free-energy landscape of pore formation. PMF calculations along the RC provide mechanistic and energetic understanding of pore formation, hence they will be useful to rationalize the effects of membrane-active peptides, electric fields, and membrane composition on transmembrane pores.
Collapse
Affiliation(s)
- Jochen S Hub
- Theoretical Physics and Center for Biophysics, Saarland University, Saarbrücken D-66123, Germany
| |
Collapse
|
16
|
Dixit M, Lazaridis T. Free energy of hydrophilic and hydrophobic pores in lipid bilayers by free energy perturbation of a restraint. J Chem Phys 2021; 153:054101. [PMID: 32770888 DOI: 10.1063/5.0016682] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The free energy of pore formation in lipid bilayers has been previously calculated using a variety of reaction coordinates. Here, we use free energy perturbation of a cylindrical lipid exclusion restraint to compute the free energy profile as a function of pore radius in dimyristoylphosphatidylcholine (DMPC) and dioleoylphosphatidylcholine (DOPC) bilayers. Additionally restraining the headgroups to lie on the membrane surface allows us to also calculate the free energy profile of hydrophobic pores, i.e., cylindrical pores lined by acyl chains. For certain pore radii, the free energy of wetting of hydrophobic pores is calculated using the density bias method. It is found that wetting of hydrophobic pores becomes thermodynamically favorable at 5.0 Å for DMPC and 6.5 Å for DOPC, although significant barriers prevent spontaneous wetting of the latter on a nanosecond time scale. The free energy of transformation of hydrophilic pores to hydrophobic ones is also calculated using free energy perturbation of headgroup restraints along the bilayer normal. This quantity, along with wetting and pore growth free energies, provides complete free energy profiles as a function of radius. Pore line tension values for the hydrophilic pores obtained from the slope of the free energy profiles are 37.6 pN for DMPC and 53.7 pN for DOPC. The free energy profiles for the hydrophobic pores are analyzed in terms of elementary interfacial tensions. It is found that a positive three-phase line tension is required to explain the results. The estimated value for this three-phase line tension (51.2 pN) lies within the expected range.
Collapse
Affiliation(s)
- Mayank Dixit
- Department of Chemistry, City College of New York, 160 Convent Ave., New York, New York 10031, USA
| | - Themis Lazaridis
- Department of Chemistry, City College of New York, 160 Convent Ave., New York, New York 10031, USA
| |
Collapse
|
17
|
Rajeshwar T R, Anishkin A, Sukharev S, Vanegas JM. Mechanical Activation of MscL Revealed by a Locally Distributed Tension Molecular Dynamics Approach. Biophys J 2020; 120:232-242. [PMID: 33333032 DOI: 10.1016/j.bpj.2020.11.2274] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2020] [Revised: 11/20/2020] [Accepted: 11/30/2020] [Indexed: 02/02/2023] Open
Abstract
Membrane tension perceived by mechanosensitive (MS) proteins mediates cellular responses to mechanical stimuli and osmotic stresses, and it also guides multiple biological functions including cardiovascular control and development. In bacteria, MS channels function as tension-activated pores limiting excessive turgor pressure, with MS channel of large conductance (MscL) acting as an emergency release valve preventing cell lysis. Previous attempts to simulate gating transitions in MscL by either directly applying steering forces to the protein or by increasing the whole-system tension were not fully successful and often disrupted the integrity of the system. We present a novel, to our knowledge, locally distributed tension molecular dynamics (LDT-MD) simulation method that allows application of forces continuously distributed among lipids surrounding the channel using a specially constructed collective variable. We report reproducible and reversible transitions of MscL to the open state with measured parameters of lateral expansion and conductivity that exactly satisfy experimental values. The LDT-MD method enables exploration of the MscL-gating process with different pulling velocities and variable tension asymmetry between the inner and outer membrane leaflets. We use LDT-MD in combination with well-tempered metadynamics to reconstruct the tension-dependent free-energy landscape for the opening transition in MscL. The flexible definition of the LDT collective variable allows general application of our method to study mechanical activation of any membrane-embedded protein.
Collapse
Affiliation(s)
| | - Andriy Anishkin
- Department of Biology, University of Maryland, College Park, Maryland
| | - Sergei Sukharev
- Department of Biology, University of Maryland, College Park, Maryland
| | - Juan M Vanegas
- Department of Physics, University of Vermont, Burlington, Vermont.
| |
Collapse
|
18
|
Caparotta M, Tomes CN, Mayorga LS, Masone D. The Synaptotagmin-1 C2B Domain Is a Key Regulator in the Stabilization of the Fusion Pore. J Chem Theory Comput 2020; 16:7840-7851. [DOI: 10.1021/acs.jctc.0c00734] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Marcelo Caparotta
- Facultad de Ciencias Exactas y Naturales, Universidad Nacional de Cuyo (UNCuyo), Mendoza 5500, Argentina
| | - Claudia N. Tomes
- Facultad de Ciencias Exactas y Naturales, Universidad Nacional de Cuyo (UNCuyo), Mendoza 5500, Argentina
- Instituto de Histología y Embriología de Mendoza (IHEM) - Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Universidad Nacional de Cuyo (UNCuyo), Mendoza 5500, Argentina
| | - Luis S. Mayorga
- Facultad de Ciencias Exactas y Naturales, Universidad Nacional de Cuyo (UNCuyo), Mendoza 5500, Argentina
- Instituto de Histología y Embriología de Mendoza (IHEM) - Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Universidad Nacional de Cuyo (UNCuyo), Mendoza 5500, Argentina
| | - Diego Masone
- Instituto de Histología y Embriología de Mendoza (IHEM) - Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Universidad Nacional de Cuyo (UNCuyo), Mendoza 5500, Argentina
- Facultad de Ingeniería, Universidad Nacional de Cuyo (UNCuyo), Mendoza 5500, Argentina
| |
Collapse
|
19
|
Timmons PB, Hewage CM. HAPPENN is a novel tool for hemolytic activity prediction for therapeutic peptides which employs neural networks. Sci Rep 2020; 10:10869. [PMID: 32616760 PMCID: PMC7331684 DOI: 10.1038/s41598-020-67701-3] [Citation(s) in RCA: 42] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2020] [Accepted: 06/09/2020] [Indexed: 12/11/2022] Open
Abstract
The growing prevalence of resistance to antibiotics motivates the search for new antibacterial agents. Antimicrobial peptides are a diverse class of well-studied membrane-active peptides which function as part of the innate host defence system, and form a promising avenue in antibiotic drug research. Some antimicrobial peptides exhibit toxicity against eukaryotic membranes, typically characterised by hemolytic activity assays, but currently, the understanding of what differentiates hemolytic and non-hemolytic peptides is limited. This study leverages advances in machine learning research to produce a novel artificial neural network classifier for the prediction of hemolytic activity from a peptide's primary sequence. The classifier achieves best-in-class performance, with cross-validated accuracy of [Formula: see text] and Matthews correlation coefficient of 0.71. This innovative classifier is available as a web server at https://research.timmons.eu/happenn , allowing the research community to utilise it for in silico screening of peptide drug candidates for high therapeutic efficacies.
Collapse
Affiliation(s)
- Patrick Brendan Timmons
- UCD School of Biomolecular and Biomedical Science, UCD Centre for Synthesis and Chemical Biology, UCD Conway Institute, University College Dublin, Dublin 4, Ireland
| | - Chandralal M Hewage
- UCD School of Biomolecular and Biomedical Science, UCD Centre for Synthesis and Chemical Biology, UCD Conway Institute, University College Dublin, Dublin 4, Ireland.
| |
Collapse
|
20
|
An D, Thiyagarajan S, Antipov E, Alcott B, O’shaughnessy B. Membrane pore energetics and the pathways to membrane rupture.. [DOI: 10.1101/2020.06.29.178988] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
AbstractBiological membranes owe their strength and low permeability to the phospholipid bilayers at their core. Membrane strength is determined by the energetics and dynamics of membrane pores, whose tension-dependent nucleation and growth leads to rupture. Creation of nanoscale membrane pores is central to exocytosis, trafficking and other processes fundamental to life that require breaching of secure plasma or organelle membranes, and is the basis for biotechnologies using drug delivery, delivery of genetic material for gene editing and antimicrobial peptides. A prevailing view from seminal electroporation and membrane rupture studies is that pore growth and bilayer rupture are controlled by macroscopically long-lived metastable defect states that precede fully developed pores. It was argued that defect nucleation becomes rate-limiting at high tensions, explaining the exponential tension-dependence of rupture times [E. Evans et al., Biophys. J. 85, 2342-2350 (2003)]. Here we measured membrane pore free energies and bilayer rupture using highly coarse-grained simulations that probe very long time scales. We find no evidence of metastable pore states. At lower tensions, small hydrophobic pores mature into large hydrophilic pores on the pathway to rupture, with classical tension dependence of rupture times. Above a critical tension membranes rupture directly from a small hydrophobic pore, and rupture times depend exponentially on tension. Thus, we recover the experimentally reported regimes, but the origin of the high tension exponential regime is unrelated to macroscopically long-lived pre-pore defects. It arises because hydrophilic pores cannot exist above a critical tension, leading to radically altered pore dynamics and rupture kinetics.
Collapse
|
21
|
Ramírez PG, Del Pópolo MG, Vila JA, Longo GS. Thermodynamics of cell penetrating peptides on lipid membranes: sequence and membrane acidity regulate surface binding. Phys Chem Chem Phys 2020; 22:23399-23410. [DOI: 10.1039/d0cp02770g] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Acidic lipids respond to pH in ways that fully promote or deplete the surface accumulation of cell penetrating peptides.
Collapse
Affiliation(s)
- Pedro G. Ramírez
- Instituto de Matemática Aplicada San Luis (IMASL)
- UNSL-CONICET
- San Luis
- Argentina
| | - Mario G. Del Pópolo
- Instituto Interdisciplinario de Ciencias Básicas (ICB-CONICET) & Facultad de Ciencias Exactas y Naturales (FCEN-UNCuyo)
- Mendoza
- Argentina
| | - Jorge A. Vila
- Instituto de Matemática Aplicada San Luis (IMASL)
- UNSL-CONICET
- San Luis
- Argentina
| | - Gabriel S. Longo
- Instituto de Investigaciones Fisicoquímicas
- Teóricas y Aplicadas (INIFTA)
- UNLP-CONICET
- La Plata
- Argentina
| |
Collapse
|
22
|
Yao C, Kang Z, Yu B, Chen Q, Liu Y, Wang Q. All-Factor Analysis and Correlations on the Transmembrane Process for Arginine-Rich Cell-Penetrating Peptides. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2019; 35:9286-9296. [PMID: 31265309 DOI: 10.1021/acs.langmuir.9b01169] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Currently, arginine-rich cell-penetrating peptides (CPPs), due to their little cytotoxicity and high transmembrane efficiency, are considered as one of the important intracellular carriers. Although the mechanism of the transmembrane process for arginine-rich CPPs was proposed, the quantitative correlations and the key factors involved in this process still deserve further investigation. In this study, all-atom molecular dynamics and the umbrella sampling technique were employed to study the arginine-rich CPPs transmembrane process. In the adsorption process of CPPs from solution to the surface of the lipid bilayer, the adsorption free energy (ΔGA) is found to be linearly related to the interaction energy change (ΔEA): ΔGA = 0.0426ΔEA + 36.7, R2 = 0.92. In the CPPs transmembrane process, the transmembrane free energy barrier (ΔGB) is roughly correlated with the corresponding interaction energy change (ΔEB): ΔGB = 0.108ΔEB +135, R2 = 0.73. The multiple salt bridges of guanidinium-PO4 account for 65% of the overall interaction energy, so the increased negative charges of the lipid bilayer or more salt bridges would facilitate CPPs adsorption and transmembrane processes. Also, the increased negative charges of the lipid bilayer would reduce the amount of water to be carried into the pore and further reduce the ΔGB. The peptide backbone would not have a direct impact on transmembrane efficiency. The ΔGB is also found to be related to the length of the pore (L): ΔGB = 46.2L - 31.3, R2 = 0.92, which makes the transmembrane efficiency estimable. This work is expected to deliver an in-depth understanding and help the optimization of CPPs.
Collapse
Affiliation(s)
- Cai Yao
- Department of Chemistry and Soft Matter Research Center , Zhejiang University , Hangzhou 310027 , China
| | - Zhengzhong Kang
- Department of Chemistry and Soft Matter Research Center , Zhejiang University , Hangzhou 310027 , China
| | - Bin Yu
- Department of Chemistry and Soft Matter Research Center , Zhejiang University , Hangzhou 310027 , China
| | - Qu Chen
- School of Biological and Chemical Engineering , Zhejiang University of Science and Technology , Hangzhou 310023 , China
| | - Yingchun Liu
- Department of Chemistry and Soft Matter Research Center , Zhejiang University , Hangzhou 310027 , China
| | - Qi Wang
- Department of Chemistry and Soft Matter Research Center , Zhejiang University , Hangzhou 310027 , China
| |
Collapse
|
23
|
Kirsch SA, Böckmann RA. Coupling of Membrane Nanodomain Formation and Enhanced Electroporation near Phase Transition. Biophys J 2019; 116:2131-2148. [PMID: 31103234 PMCID: PMC6554532 DOI: 10.1016/j.bpj.2019.04.024] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2018] [Revised: 04/17/2019] [Accepted: 04/18/2019] [Indexed: 12/29/2022] Open
Abstract
Biological cells are enveloped by a heterogeneous lipid bilayer that prevents the uncontrolled exchange of substances between the cell interior and its environment. In particular, membranes act as a continuous barrier for salt and macromolecules to ensure proper physiological functions within the cell. However, it has been shown that membrane permeability strongly depends on temperature and, for phospholipid bilayers, displays a maximum at the transition between the gel and fluid phase. Here, extensive molecular dynamics simulations of dipalmitoylphosphatidylcholine bilayers were employed to characterize the membrane structure and dynamics close to phase transition, as well as its stability with respect to an external electric field. Atomistic simulations revealed the dynamic appearance and disappearance of spatially related nanometer-sized thick ordered and thin interdigitating domains in a fluid-like bilayer close to the phase transition temperature (Tm). These structures likely represent metastable precursors of the ripple phase that vanished at increased temperatures. Similarly, a two-phase bilayer with coexisting gel and fluid domains featured a thickness minimum at the interface because of splaying and interdigitating lipids. For all systems, application of an external electric field revealed a reduced bilayer stability with respect to pore formation for temperatures close to Tm. Pore formation occurred exclusively in thin interdigitating membrane nanodomains. These findings provide a link between the increased membrane permeability and the structural heterogeneity close to phase transition.
Collapse
Affiliation(s)
- Sonja A Kirsch
- Computational Biology, Friedrich-Alexander University of Erlangen-Nürnberg, Erlangen, Germany
| | - Rainer A Böckmann
- Computational Biology, Friedrich-Alexander University of Erlangen-Nürnberg, Erlangen, Germany.
| |
Collapse
|
24
|
Hong J, Lu X, Deng Z, Xiao S, Yuan B, Yang K. How Melittin Inserts into Cell Membrane: Conformational Changes, Inter-Peptide Cooperation, and Disturbance on the Membrane. Molecules 2019; 24:molecules24091775. [PMID: 31067828 PMCID: PMC6539814 DOI: 10.3390/molecules24091775] [Citation(s) in RCA: 55] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2019] [Revised: 05/02/2019] [Accepted: 05/03/2019] [Indexed: 01/27/2023] Open
Abstract
Antimicrobial peptides (AMPs), as a key component of the immune defense systems of organisms, are a promising solution to the serious threat of drug-resistant bacteria to public health. As one of the most representative and extensively studied AMPs, melittin has exceptional broad-spectrum activities against microorganisms, including both Gram-positive and Gram-negative bacteria. Unfortunately, the action mechanism of melittin with bacterial membranes, especially the underlying physics of peptide-induced membrane poration behaviors, is still poorly understood, which hampers efforts to develop melittin-based drugs or agents for clinical applications. In this mini-review, we focus on recent advances with respect to the membrane insertion behavior of melittin mostly from a computational aspect. Membrane insertion is a prerequisite and key step for forming transmembrane pores and bacterial killing by melittin, whose occurrence is based on overcoming a high free-energy barrier during the transition of melittin molecules from a membrane surface-binding state to a transmembrane-inserting state. Here, intriguing simulation results on such transition are highlighted from both kinetic and thermodynamic aspects. The conformational changes and inter-peptide cooperation of melittin molecules, as well as melittin-induced disturbances to membrane structure, such as deformation and lipid extraction, are regarded as key factors influencing the insertion of peptides into membranes. The associated intermediate states in peptide conformations, lipid arrangements, membrane structure, and mechanical properties during this process are specifically discussed. Finally, potential strategies for enhancing the poration ability and improving the antimicrobial performance of AMPs are included as well.
Collapse
Affiliation(s)
- Jiajia Hong
- Center for Soft Condensed Matter Physics and Interdisciplinary Research & School of Physical Science and Technology, Soochow University, Suzhou 215006, China.
| | - Xuemei Lu
- Center for Soft Condensed Matter Physics and Interdisciplinary Research & School of Physical Science and Technology, Soochow University, Suzhou 215006, China.
| | - Zhixiong Deng
- Center for Soft Condensed Matter Physics and Interdisciplinary Research & School of Physical Science and Technology, Soochow University, Suzhou 215006, China.
| | - Shufeng Xiao
- Center for Soft Condensed Matter Physics and Interdisciplinary Research & School of Physical Science and Technology, Soochow University, Suzhou 215006, China.
| | - Bing Yuan
- Center for Soft Condensed Matter Physics and Interdisciplinary Research & School of Physical Science and Technology, Soochow University, Suzhou 215006, China.
| | - Kai Yang
- Center for Soft Condensed Matter Physics and Interdisciplinary Research & School of Physical Science and Technology, Soochow University, Suzhou 215006, China.
| |
Collapse
|
25
|
Miyazaki Y, Okazaki S, Shinoda W. Free energy analysis of membrane pore formation process in the presence of multiple melittin peptides. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2019; 1861:1409-1419. [PMID: 30885804 DOI: 10.1016/j.bbamem.2019.03.002] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/31/2018] [Revised: 03/01/2019] [Accepted: 03/04/2019] [Indexed: 11/30/2022]
Abstract
Understanding the molecular mechanism underlying pore formation in lipid membranes by antimicrobial peptides is of great importance in biological sciences as well as in drug design applications. Melittin has been widely studied as a pore forming peptide, though the molecular mechanism for pore formation is still illusive. We examined the free energy barrier for the creation of a pore in lipid membranes with and without multiple melittin peptides. It was found that six melittin peptides significantly stabilized a pore, though a small barrier (a few kBT) for the formation still existed. With five melittin peptides or fewer, the pore formation barrier was much higher, though the established pore was in a local energy minimum. Although seven melittins effectively reduced the free energy barrier, a single melittin peptide left the pore after a long time MD simulation probably because of the overcrowded environment around the bilayer pore. Thus, it is highly selective for the number of melittin peptides to stabilize the membrane pore, as was also suggested by the line tension evaluations. The free energy cost required to insert a single melittin into the membrane is too high to explain the one-by-one insertion mechanism for pore formation, which also supports the collective melittin mechanism for pore formation.
Collapse
Affiliation(s)
- Yusuke Miyazaki
- Department of Materials Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, 464-8603, Japan
| | - Susumu Okazaki
- Department of Materials Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, 464-8603, Japan
| | - Wataru Shinoda
- Department of Materials Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, 464-8603, Japan.
| |
Collapse
|
26
|
Allolio C, Magarkar A, Jurkiewicz P, Baxová K, Javanainen M, Mason PE, Šachl R, Cebecauer M, Hof M, Horinek D, Heinz V, Rachel R, Ziegler CM, Schröfel A, Jungwirth P. Arginine-rich cell-penetrating peptides induce membrane multilamellarity and subsequently enter via formation of a fusion pore. Proc Natl Acad Sci U S A 2018; 115:11923-11928. [PMID: 30397112 PMCID: PMC6255155 DOI: 10.1073/pnas.1811520115] [Citation(s) in RCA: 141] [Impact Index Per Article: 23.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022] Open
Abstract
Arginine-rich cell-penetrating peptides do not enter cells by directly passing through a lipid membrane; they instead passively enter vesicles and live cells by inducing membrane multilamellarity and fusion. The molecular picture of this penetration mode, which differs qualitatively from the previously proposed direct mechanism, is provided by molecular dynamics simulations. The kinetics of vesicle agglomeration and fusion by an iconic cell-penetrating peptide-nonaarginine-are documented via real-time fluorescence techniques, while the induction of multilamellar phases in vesicles and live cells is demonstrated by a combination of electron and fluorescence microscopies. This concert of experiments and simulations reveals that the identified passive cell penetration mechanism bears analogy to vesicle fusion induced by calcium ions, indicating that the two processes may share a common mechanistic origin.
Collapse
Affiliation(s)
- Christoph Allolio
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, CZ-166 10 Prague 6, Czech Republic
- Institute of Physical and Theoretical Chemistry, University of Regensburg, D-93040 Regensburg, Germany
- Fritz Haber Research Center, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
- Department of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
| | - Aniket Magarkar
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, CZ-166 10 Prague 6, Czech Republic
- Faculty of Pharmacy, University of Helsinki, Helsinki 00014, Finland
| | - Piotr Jurkiewicz
- J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, 182 23 Prague 8, Czech Republic
| | - Katarína Baxová
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, CZ-166 10 Prague 6, Czech Republic
| | - Matti Javanainen
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, CZ-166 10 Prague 6, Czech Republic
| | - Philip E Mason
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, CZ-166 10 Prague 6, Czech Republic
| | - Radek Šachl
- J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, 182 23 Prague 8, Czech Republic
| | - Marek Cebecauer
- J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, 182 23 Prague 8, Czech Republic
| | - Martin Hof
- J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, 182 23 Prague 8, Czech Republic
| | - Dominik Horinek
- Institute of Physical and Theoretical Chemistry, University of Regensburg, D-93040 Regensburg, Germany
| | - Veronika Heinz
- Institute of Biophysics and Biophysical Chemistry, University of Regensburg, D-93040 Regensburg, Germany
| | - Reinhard Rachel
- Microbiology and Archaea Centre, University of Regensburg, D-93040 Regensburg, Germany
| | - Christine M Ziegler
- Institute of Biophysics and Biophysical Chemistry, University of Regensburg, D-93040 Regensburg, Germany
- Institute of Biophysics and Biophysical Chemistry, University of Regensburg, D-93040 Regensburg, Germany
| | - Adam Schröfel
- Imaging Methods Core Facility at Biocev, Faculty of Sciences, Charles University, 242 50 Vestec, Czech Republic
| | - Pavel Jungwirth
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, CZ-166 10 Prague 6, Czech Republic;
| |
Collapse
|
27
|
Allolio C, Magarkar A, Jurkiewicz P, Baxová K, Javanainen M, Mason PE, Šachl R, Cebecauer M, Hof M, Horinek D, Heinz V, Rachel R, Ziegler CM, Schröfel A, Jungwirth P. Arginine-rich cell-penetrating peptides induce membrane multilamellarity and subsequently enter via formation of a fusion pore. Proc Natl Acad Sci U S A 2018. [PMID: 30397112 DOI: 10.1073/pnas.1811520115/suppl_file/pnas.1811520115.sm01.mp4] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/19/2023] Open
Abstract
Arginine-rich cell-penetrating peptides do not enter cells by directly passing through a lipid membrane; they instead passively enter vesicles and live cells by inducing membrane multilamellarity and fusion. The molecular picture of this penetration mode, which differs qualitatively from the previously proposed direct mechanism, is provided by molecular dynamics simulations. The kinetics of vesicle agglomeration and fusion by an iconic cell-penetrating peptide-nonaarginine-are documented via real-time fluorescence techniques, while the induction of multilamellar phases in vesicles and live cells is demonstrated by a combination of electron and fluorescence microscopies. This concert of experiments and simulations reveals that the identified passive cell penetration mechanism bears analogy to vesicle fusion induced by calcium ions, indicating that the two processes may share a common mechanistic origin.
Collapse
Affiliation(s)
- Christoph Allolio
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, CZ-166 10 Prague 6, Czech Republic
- Institute of Physical and Theoretical Chemistry, University of Regensburg, D-93040 Regensburg, Germany
- Fritz Haber Research Center, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
- Department of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
| | - Aniket Magarkar
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, CZ-166 10 Prague 6, Czech Republic
- Faculty of Pharmacy, University of Helsinki, Helsinki 00014, Finland
| | - Piotr Jurkiewicz
- J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, 182 23 Prague 8, Czech Republic
| | - Katarína Baxová
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, CZ-166 10 Prague 6, Czech Republic
| | - Matti Javanainen
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, CZ-166 10 Prague 6, Czech Republic
| | - Philip E Mason
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, CZ-166 10 Prague 6, Czech Republic
| | - Radek Šachl
- J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, 182 23 Prague 8, Czech Republic
| | - Marek Cebecauer
- J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, 182 23 Prague 8, Czech Republic
| | - Martin Hof
- J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, 182 23 Prague 8, Czech Republic
| | - Dominik Horinek
- Institute of Physical and Theoretical Chemistry, University of Regensburg, D-93040 Regensburg, Germany
| | - Veronika Heinz
- Institute of Biophysics and Biophysical Chemistry, University of Regensburg, D-93040 Regensburg, Germany
| | - Reinhard Rachel
- Microbiology and Archaea Centre, University of Regensburg, D-93040 Regensburg, Germany
| | - Christine M Ziegler
- Institute of Biophysics and Biophysical Chemistry, University of Regensburg, D-93040 Regensburg, Germany
- Institute of Biophysics and Biophysical Chemistry, University of Regensburg, D-93040 Regensburg, Germany
| | - Adam Schröfel
- Imaging Methods Core Facility at Biocev, Faculty of Sciences, Charles University, 242 50 Vestec, Czech Republic
| | - Pavel Jungwirth
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, CZ-166 10 Prague 6, Czech Republic;
| |
Collapse
|
28
|
Holme M, Rashid MH, Thomas MR, Barriga HMG, Herpoldt K, Heenan RK, Dreiss CA, Bañuelos JL, Xie HN, Yarovsky I, Stevens MM. Fate of Liposomes in the Presence of Phospholipase C and D: From Atomic to Supramolecular Lipid Arrangement. ACS CENTRAL SCIENCE 2018; 4:1023-1030. [PMID: 30159399 PMCID: PMC6107861 DOI: 10.1021/acscentsci.8b00286] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/06/2018] [Indexed: 05/04/2023]
Abstract
Understanding the origins of lipid membrane bilayer rearrangement in response to external stimuli is an essential component of cell biology and the bottom-up design of liposomes for biomedical applications. The enzymes phospholipase C and D (PLC and PLD) both cleave the phosphorus-oxygen bonds of phosphate esters in phosphatidylcholine (PC) lipids. The atomic position of this hydrolysis reaction has huge implications for the stability of PC-containing self-assembled structures, such as the cell wall and lipid-based vesicle drug delivery vectors. While PLC converts PC to diacylglycerol (DAG), the interaction of PC with PLD produces phosphatidic acid (PA). Here we present a combination of small-angle scattering data and all-atom molecular dynamics simulations, providing insights into the effects of atomic-scale reorganization on the supramolecular assembly of PC membrane bilayers upon enzyme-mediated incorporation of DAG or PA. We observed that PC liposomes completely disintegrate in the presence of PLC, as conversion of PC to DAG progresses. At lower concentrations, DAG molecules within fluid PC bilayers form hydrogen bonds with backbone carbonyl oxygens in neighboring PC molecules and burrow into the hydrophobic region. This leads initially to membrane thinning followed by a swelling of the lamellar phase with increased DAG. At higher DAG concentrations, localized membrane tension causes a change in lipid phase from lamellar to the hexagonal and micellar cubic phases. Molecular dynamics simulations show that this destabilization is also caused in part by the decreased ability of DAG-containing PC membranes to coordinate sodium ions. Conversely, PLD-treated PC liposomes remain stable up to extremely high conversions to PA. Here, the negatively charged PA headgroup attracts significant amounts of sodium ions from the bulk solution to the membrane surface, leading to a swelling of the coordinated water layer. These findings are a vital step toward a fundamental understanding of the degradation behavior of PC lipid membranes in the presence of these clinically relevant enzymes, and toward the rational design of diagnostic and drug delivery technologies for phospholipase-dysregulation-based diseases.
Collapse
Affiliation(s)
- Margaret
N. Holme
- Department
of Materials, Imperial College London, London SW7 2AZ, United Kingdom
- Department of Bioengineering and Institute of Biomedical
Engineering, Imperial College London, London SW7 2AZ, United Kingdom
- Department
of Medical Biochemistry and Biophysics, Karolinska Institutet, SE-171 77 Stockholm, Sweden
| | - M. Harunur Rashid
- School
of Engineering, RMIT University, Melbourne, Victoria 3001, Australia
| | - Michael R. Thomas
- Department
of Materials, Imperial College London, London SW7 2AZ, United Kingdom
- Department of Bioengineering and Institute of Biomedical
Engineering, Imperial College London, London SW7 2AZ, United Kingdom
| | - Hanna M. G. Barriga
- Department
of Medical Biochemistry and Biophysics, Karolinska Institutet, SE-171 77 Stockholm, Sweden
| | - Karla−Luise Herpoldt
- Department
of Materials, Imperial College London, London SW7 2AZ, United Kingdom
- Department of Bioengineering and Institute of Biomedical
Engineering, Imperial College London, London SW7 2AZ, United Kingdom
| | - Richard K. Heenan
- STFC ISIS
Facility, Rutherford Appleton Laboratory, Chilton, Didcot OX11 0QX, United Kingdom
| | - Cécile A. Dreiss
- School
of Cancer and Pharmaceutical Sciences, King’s
College London, London SE1 9NH, United Kingdom
| | - José Leobardo Bañuelos
- STFC ISIS
Facility, Rutherford Appleton Laboratory, Chilton, Didcot OX11 0QX, United Kingdom
- Department
of Physics, The University of Texas at El
Paso, El Paso, Texas 79968, United States
| | - Hai-nan Xie
- Department
of Materials, Imperial College London, London SW7 2AZ, United Kingdom
- Department of Bioengineering and Institute of Biomedical
Engineering, Imperial College London, London SW7 2AZ, United Kingdom
| | - Irene Yarovsky
- School
of Engineering, RMIT University, Melbourne, Victoria 3001, Australia
- E-mail:
| | - Molly M. Stevens
- Department
of Materials, Imperial College London, London SW7 2AZ, United Kingdom
- Department of Bioengineering and Institute of Biomedical
Engineering, Imperial College London, London SW7 2AZ, United Kingdom
- Department
of Medical Biochemistry and Biophysics, Karolinska Institutet, SE-171 77 Stockholm, Sweden
- E-mail:
| |
Collapse
|
29
|
Grasso G, Muscat S, Rebella M, Morbiducci U, Audenino A, Danani A, Deriu MA. Cell penetrating peptide modulation of membrane biomechanics by Molecular dynamics. J Biomech 2018; 73:137-144. [PMID: 29631749 DOI: 10.1016/j.jbiomech.2018.03.036] [Citation(s) in RCA: 33] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2017] [Revised: 03/14/2018] [Accepted: 03/21/2018] [Indexed: 11/17/2022]
Abstract
The efficacy of a pharmaceutical treatment is often countered by the inadequate membrane permeability, that prevents drugs from reaching their specific intracellular targets. Cell penetrating peptides (CPPs) are able to route across cells' membrane various types of cargo, including drugs and nanoparticles. However, CPPs internalization mechanisms are not yet fully understood and depend on a wide variety of aspects. In this contest, the entry of a CPP into the lipid bilayer might induce molecular conformational changes, including marked variations on membrane's mechanical properties. Understanding how the CPP does influence the mechanical properties of cells membrane is crucial to design, engineer and improve new and existing penetrating peptides. Here, all atom Molecular Dynamics (MD) simulations were used to investigate the interaction between different types of CPPs embedded in a lipid bilayer of dioleoyl phosphatidylcholine (DOPC). In a greater detail, we systematically highlighted how CPP properties are responsible for modulating the membrane bending modulus. Our findings highlighted the CPP hydropathy strongly correlated with penetration of water molecules in the lipid bilayer, thus supporting the hypothesis that the amount of water each CPP can route inside the membrane is modulated by the hydrophobic and hydrophilic character of the peptide. Water penetration promoted by CPPs leads to a local decrease of the lipid order, which emerges macroscopically as a reduction of the membrane bending modulus.
Collapse
Affiliation(s)
- Gianvito Grasso
- Istituto Dalle Molle di Studi sull'Intelligenza Artificiale (IDSIA), Scuola universitaria professionale della Svizzera italiana (SUPSI), Università della Svizzera Italiana (USI), Centro Galleria 2, Manno CH-6928, Switzerland
| | - Stefano Muscat
- Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, IT-10128 Torino, Italy
| | - Martina Rebella
- Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, IT-10128 Torino, Italy
| | - Umberto Morbiducci
- Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, IT-10128 Torino, Italy
| | - Alberto Audenino
- Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, IT-10128 Torino, Italy
| | - Andrea Danani
- Istituto Dalle Molle di Studi sull'Intelligenza Artificiale (IDSIA), Scuola universitaria professionale della Svizzera italiana (SUPSI), Università della Svizzera Italiana (USI), Centro Galleria 2, Manno CH-6928, Switzerland
| | - Marco A Deriu
- Istituto Dalle Molle di Studi sull'Intelligenza Artificiale (IDSIA), Scuola universitaria professionale della Svizzera italiana (SUPSI), Università della Svizzera Italiana (USI), Centro Galleria 2, Manno CH-6928, Switzerland.
| |
Collapse
|
30
|
Perrier DL, Rems L, Boukany PE. Lipid vesicles in pulsed electric fields: Fundamental principles of the membrane response and its biomedical applications. Adv Colloid Interface Sci 2017; 249:248-271. [PMID: 28499600 DOI: 10.1016/j.cis.2017.04.016] [Citation(s) in RCA: 43] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2017] [Revised: 04/24/2017] [Accepted: 04/25/2017] [Indexed: 01/04/2023]
Abstract
The present review focuses on the effects of pulsed electric fields on lipid vesicles ranging from giant unilamellar vesicles (GUVs) to small unilamellar vesicles (SUVs), from both fundamental and applicative perspectives. Lipid vesicles are the most popular model membrane systems for studying biophysical and biological processes in living cells. Furthermore, as vesicles are made from biocompatible and biodegradable materials, they provide a strategy to create safe and functionalized drug delivery systems in health-care applications. Exposure of lipid vesicles to pulsed electric fields is a common physical method to transiently increase the permeability of the lipid membrane. This method, termed electroporation, has shown many advantages for delivering exogenous molecules including drugs and genetic material into vesicles and living cells. In addition, electroporation can be applied to induce fusion between vesicles and/or cells. First, we discuss in detail how research on cell-size GUVs as model cell systems has provided novel insight into the basic mechanisms of cell electroporation and associated phenomena. Afterwards, we continue with a thorough overview how electroporation and electrofusion have been used as versatile methods to manipulate vesicles of all sizes in different biomedical applications. We conclude by summarizing the open questions in the field of electroporation and possible future directions for vesicles in the biomedical field.
Collapse
|
31
|
Hub JS, Awasthi N. Probing a Continuous Polar Defect: A Reaction Coordinate for Pore Formation in Lipid Membranes. J Chem Theory Comput 2017; 13:2352-2366. [PMID: 28376619 DOI: 10.1021/acs.jctc.7b00106] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Various biophysical processes involve the formation of aqueous pores over lipid membranes, including processes of membrane fusion, antimicrobial peptide activity, lipid flip-flop, and membrane permeation. Reliable and efficient free-energy calculations of pore formation using molecular dynamics simulations remained challenging due to the lack of good reaction coordinates (RCs) for pore formation. We present a new RC for pore formation that probes the formation and rupture of a continuous polar defect over the membrane. Potential of mean force (PMF) calculations along the new RC rapidly converge and exhibit no hysteresis between pore-opening and pore-closing pathways, in contrast to calculations based on previous RCs. We show that restraints along the new RC may restrain the system tightly to the transition state of pore formation, rationalizing the absence of hysteresis. We observe that the PMF of pore formation in a tension-free membrane of dimyristoylphosphatidylcholine (DMPC) reveals a free-energy barrier for pore nucleation, confirming a long-hypothesized metastable prepore state. We test the influence of the lipid force field, the cutoff distance used for Lennard-Jones interactions, and the lateral membrane size on the free energies of pore formation. In contrast to PMF calculations based on previous RCs, we find that such parameters have a relatively small influence on the free energies of pore nucleation. However, the metastability of the open pore in DMPC may depend on such parameters. The RC has been implemented into an extension of the GROMACS simulation software. The new RC allows for reliable and computationally efficient free-energy calculations of pore formation in lipid membranes.
Collapse
Affiliation(s)
- Jochen S Hub
- Institute for Microbiology and Genetics, University of Göttingen , Justus-von-Liebig Weg 11, 37077 Goettingen, Germany
| | - Neha Awasthi
- Institute for Microbiology and Genetics, University of Göttingen , Justus-von-Liebig Weg 11, 37077 Goettingen, Germany
| |
Collapse
|
32
|
Woo SY, Lee H. Aggregation and insertion of melittin and its analogue MelP5 into lipid bilayers at different concentrations: effects on pore size, bilayer thickness and dynamics. Phys Chem Chem Phys 2017; 19:7195-7203. [DOI: 10.1039/c6cp06834k] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Melittin and its analogue MelP5 (five mutations T10A, R22A, K23A, R24Q, and Q26L of melittin) were simulated with lipid bilayers at different peptide/lipid molar ratios using all-atom and coarse-grained (CG) force fields.
Collapse
Affiliation(s)
- Sun Young Woo
- Department of Chemical Engineering
- Dankook University
- Yongin
- South Korea
| | - Hwankyu Lee
- Department of Chemical Engineering
- Dankook University
- Yongin
- South Korea
| |
Collapse
|
33
|
Rems L. Applicative Use of Electroporation Models. ADVANCES IN BIOMEMBRANES AND LIPID SELF-ASSEMBLY 2017. [DOI: 10.1016/bs.abl.2017.06.001] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
|
34
|
Lopes D, Jakobtorweihen S, Nunes C, Sarmento B, Reis S. Shedding light on the puzzle of drug-membrane interactions: Experimental techniques and molecular dynamics simulations. Prog Lipid Res 2017; 65:24-44. [DOI: 10.1016/j.plipres.2016.12.001] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2016] [Revised: 11/30/2016] [Accepted: 12/03/2016] [Indexed: 12/20/2022]
|
35
|
Woo SY, Lee H. Effect of lipid shape on toroidal pore formation and peptide orientation in lipid bilayers. Phys Chem Chem Phys 2017; 19:21340-21349. [DOI: 10.1039/c7cp02708g] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Disordered and thinner bilayer w/lyso-lipids; tilted orientation of peptides in bilayer w/lyso-lipids; toroidal pores stabilized by peptides and lyso-lipids.
Collapse
Affiliation(s)
- Sun Young Woo
- Department of Chemical Engineering
- Dankook University
- Yongin
- South Korea
| | - Hwankyu Lee
- Department of Chemical Engineering
- Dankook University
- Yongin
- South Korea
| |
Collapse
|
36
|
Rems L, Tarek M, Casciola M, Miklavčič D. Properties of lipid electropores II: Comparison of continuum-level modeling of pore conductance to molecular dynamics simulations. Bioelectrochemistry 2016; 112:112-24. [DOI: 10.1016/j.bioelechem.2016.03.005] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2015] [Revised: 03/24/2016] [Accepted: 03/27/2016] [Indexed: 12/18/2022]
|
37
|
Hoiles W, Gupta R, Cornell B, Cranfield C, Krishnamurthy V. The Effect of Tethers on Artificial Cell Membranes: A Coarse-Grained Molecular Dynamics Study. PLoS One 2016; 11:e0162790. [PMID: 27736860 PMCID: PMC5063460 DOI: 10.1371/journal.pone.0162790] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2015] [Accepted: 08/29/2016] [Indexed: 11/18/2022] Open
Abstract
Tethered bilayer lipid membranes (tBLMs) provide a stable platform for modeling the dynamics and order of biological membranes where the tethers mimic the cytoskeletal supports present in biological cell membranes. In this paper coarse-grained molecular dynamics (CGMD) is applied to study the effects of tethers on lipid membrane properties. Using results from the CGMD model and the overdamped Fokker-Planck equation, we show that the diffusion tensor and particle density of water in the tBLM is spatially dependent. Further, it is shown that the membrane thickness, lipid diffusion, defect density, free energy of lipid flip-flop, and membrane dielectric permittivity are all dependent on the tether density. The numerically computed results from the CGMD model are in agreement with the experimentally measured results from tBLMs containing different tether densities and lipids derived from Archaebacteria. Additionally, using experimental measurements from Escherichia coli bacteria and Saccharomyces Cerevisiae yeast tethered membranes, we illustrate how previous molecular dynamics results can be combined with the proposed model to estimate the dielectric permittivity and defect density of these membranes as a function of tether density.
Collapse
Affiliation(s)
- William Hoiles
- Department of Electrical and Computer Engineering, The University of British Columbia, Vancouver, British Columbia, Canada
| | - Rini Gupta
- Department of Electrical and Computer Engineering, The University of British Columbia, Vancouver, British Columbia, Canada
| | - Bruce Cornell
- Director of Science and Technology, Surgical Diagnostics Pty Ltd., Unit 6 30-32 Barcoo Street, Roseville, New South Wales, 2069, Australia
| | - Charles Cranfield
- School of Life Sciences, University of Technology Sydney, Broadway, New South Wales, Australia
| | - Vikram Krishnamurthy
- Electrical and Computer Engineering, Cornell University, New York, New York, United States of America
| |
Collapse
|
38
|
King MJ, Bennett AL, Almeida PF, Lee HS. Coarse-grained simulations of hemolytic peptide δ-lysin interacting with a POPC bilayer. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2016; 1858:3182-3194. [PMID: 27720634 DOI: 10.1016/j.bbamem.2016.10.004] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/10/2016] [Revised: 09/17/2016] [Accepted: 10/04/2016] [Indexed: 01/08/2023]
Abstract
δ-lysin, secreted by a Gram-positive bacterium Staphylococcus aureus, is a 26-residue membrane active peptide that shares many common features with antimicrobial peptides (AMPs). However, it possesses a few unique features that differentiate itself from typical AMPs. In particular, δ-lysin has zero net charge, even though it has many charged residues, and it preferentially lyses eukaryotic cells over bacterial cells. Here, we present the results of coarse-grained molecular dynamics simulations of δ-lysin interacting with a zwitterionic membrane over a wide range of peptide concentrations. When the peptides concentration is low, spontaneous dimerization of peptides is observed on the membrane surface, but deep insertion of peptides or pore formation was not observed. However, the calculated free energy of peptide insertion suggests that a small fraction of peptides is likely to be present inside the membrane at the peptide concentrations typically seen in dye efflux experiments. When the simulations with multiple peptides are carried out with a single pre-inserted transmembrane peptide, spontaneous pore formation occurs with a peptide-to-lipid ratio (P/L) as low as P/L=1:42. Inter-peptide salt bridges among the transmembrane peptides seem to play a role in creating compact pores with very low level of hydration. More importantly, the transmembrane peptides making up the pore are constantly pushed to the opposite side of the membrane when the mass imbalance between the two sides of membrane is significant. Thus, the pore is very dynamic, allowing multiple peptides to translocate across the membrane simultaneously.
Collapse
Affiliation(s)
- Mariah J King
- Department of Chemistry and Biochemistry, University of North Carolina Wilmington, Wilmington, NC 28403, United States
| | - Ashley L Bennett
- Department of Chemistry and Biochemistry, University of North Carolina Wilmington, Wilmington, NC 28403, United States
| | - Paulo F Almeida
- Department of Chemistry and Biochemistry, University of North Carolina Wilmington, Wilmington, NC 28403, United States
| | - Hee-Seung Lee
- Department of Chemistry and Biochemistry, University of North Carolina Wilmington, Wilmington, NC 28403, United States.
| |
Collapse
|
39
|
Hu Y, Patel S. Thermodynamics of cell-penetrating HIV1 TAT peptide insertion into PC/PS/CHOL model bilayers through transmembrane pores: the roles of cholesterol and anionic lipids. SOFT MATTER 2016; 12:6716-6727. [PMID: 27435187 DOI: 10.1039/c5sm01696g] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Efficient delivery of pharmaceutically active molecules across cellular membranes using cell penetrating peptides (CPPs), such as the cationic human immunodeficiency virus-1 trans-acting activator of transcription peptide (HIV-1 TAT), continues to attract scientific attention in drug design and disease treatment. Experimental results show that the TAT peptide is not only capable of directly penetrating the biological membrane in a passive manner, but also forming physical, membrane-spanning pores that may facilitate transport. Experiments further show that anionic lipids accelerate peptide permeation within a range of mole percentage composition. In this work, we explored the structures and translocation thermodynamics of the cationic TAT peptide across a series of DPPC/DPPS model membranes with the presence of 0-30 mol% cholesterol. We computed the potentials of the mean force by using umbrella sampling molecular dynamics simulations coupled to the Martini coarse-grained force field. We systematically investigated the roles of cholesterol and anionic lipids (membrane surface charge) in TAT peptide translocation. In qualitative agreement with experimental findings, the barrier heights were significantly reduced in the presence of anionic lipids. A toroidal hydrophilic pore was strongly suggested by membrane structure analysis. Cholesterol stabilizes the liquid-ordered (Lo) phase of membranes and increases the elastic stiffness of bilayers. Consequently, it hinders transmembrane pore formation and thus modulates solute permeability, since the liquid-ordered phase suppresses reorientation of the lipid molecules on simulation time scales. Though cholesterol contributes marginally to the total free energy associated with peptide permeation, the coordination of cholesterol to the peptide weakens more favorable peptide-lipid interactions. The addition of the anionic lipid DPPS to the neutral DPPC bilayer leads to the emergence and further enhancement of an interfacially stable state of the peptide due to the favorable peptide-anionic lipid interactions. Translocation free energy barriers decrease in lockstep with increasing DPPS composition in the model bilayers simulated. Finally, we investigated the size of hydrophilic pores emerging in our simulations, as well as the qualitative mobility of the peptide on the membrane surface.
Collapse
Affiliation(s)
- Yuan Hu
- Merck Research Laboratories, 2015 Galloping Hill Road, Kenilworth, New Jersey 07033, USA
| | - Sandeep Patel
- Department of Chemistry and Biochemistry, University of Delaware, Newark, Delaware 19716, USA.
| |
Collapse
|
40
|
Velikonja A, Kramar P, Miklavčič D, Maček Lebar A. Specific electrical capacitance and voltage breakdown as a function of temperature for different planar lipid bilayers. Bioelectrochemistry 2016; 112:132-7. [PMID: 26948707 DOI: 10.1016/j.bioelechem.2016.02.009] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2015] [Revised: 02/15/2016] [Accepted: 02/23/2016] [Indexed: 12/16/2022]
Abstract
The breakdown voltage and specific electrical capacitance of planar lipid bilayers formed from lipids isolated from the membrane of archaeon Aeropyrum pernix K1 as a function of temperature were studied and compared with data obtained previously in MD simulation studies. Temperature dependence of breakdown voltage and specific electrical capacitance was measured also for dipalmitoylphosphatidylcholine (DPPC) bilayers and bilayers formed from mixture of diphytanoylphosphocholine (DPhPC) and DPPC in ratio 80:20. The breakdown voltage of archaeal lipids planar lipid bilayers is more or less constant until 50°C, while at higher temperatures a considerable drop is observed, which is in line with the results from MD simulations. The breakdown voltage of DPPC planar lipid bilayer at melting temperature is considerably higher than in the gel phase. Specific electrical capacitance of planar lipid bilayers formed from archaeal lipids is approximately constant for temperatures up to 40°C and then gradually decreases. The difference with MD simulation predictions is discussed. Specific electrical capacitance of DPPC planar lipid bilayers in fluid phase is 1.75 times larger than that of the gel phase and it follows intermediated phases before phase transition. Increase in specific electrical capacitance while approaching melting point of DPPC is visible also for DPhPC:DPPC mixture.
Collapse
Affiliation(s)
- Aljaž Velikonja
- University of Ljubljana, Faculty of Electrical Engineering, Slovenia
| | - Peter Kramar
- University of Ljubljana, Faculty of Electrical Engineering, Slovenia
| | - Damijan Miklavčič
- University of Ljubljana, Faculty of Electrical Engineering, Slovenia
| | | |
Collapse
|
41
|
Kirsch SA, Böckmann RA. Membrane pore formation in atomistic and coarse-grained simulations. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2015; 1858:2266-2277. [PMID: 26748016 DOI: 10.1016/j.bbamem.2015.12.031] [Citation(s) in RCA: 47] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/15/2015] [Revised: 12/23/2015] [Accepted: 12/24/2015] [Indexed: 12/26/2022]
Abstract
Biological cells and their organelles are protected by ultra thin membranes. These membranes accomplish a broad variety of important tasks like separating the cell content from the outer environment, they are the site for cell-cell interactions and many enzymatic reactions, and control the in- and efflux of metabolites. For certain physiological functions e.g. in the fusion of membranes and also in a number of biotechnological applications like gene transfection the membrane integrity needs to be compromised to allow for instance for the exchange of polar molecules across the membrane barrier. Mechanisms enabling the transport of molecules across the membrane involve membrane proteins that form specific pores or act as transporters, but also so-called lipid pores induced by external fields, stress, or peptides. Recent progress in the simulation field enabled to closely mimic pore formation as supposed to occur in vivo or in vitro. Here, we review different simulation-based approaches in the study of membrane pores with a focus on lipid pore properties such as their size and energetics, poration mechanisms based on the application of external fields, charge imbalances, or surface tension, and on pores that are induced by small molecules, peptides, and lipids. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.
Collapse
Affiliation(s)
- Sonja A Kirsch
- Computational Biology, Department of Biology, Friedrich-Alexander Universität Erlangen-Nürnberg, Erlangen, Germany
| | - Rainer A Böckmann
- Computational Biology, Department of Biology, Friedrich-Alexander Universität Erlangen-Nürnberg, Erlangen, Germany.
| |
Collapse
|
42
|
Allolio C, Baxova K, Vazdar M, Jungwirth P. Guanidinium Pairing Facilitates Membrane Translocation. J Phys Chem B 2015; 120:143-53. [DOI: 10.1021/acs.jpcb.5b10404] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Christoph Allolio
- Institute
of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nám. 2, CZ-166 10 Prague 6, Czech Republic
- Institut
für Physikalische and Theoretische Chemie, Universität Regensburg, 93040 Regensburg, Germany
| | - Katarina Baxova
- Institute
of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nám. 2, CZ-166 10 Prague 6, Czech Republic
| | - Mario Vazdar
- Institut
Rudjer
Bošković, Bijenička
cesta 54, 10000 Zagreb, Croatia
| | - Pavel Jungwirth
- Institute
of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nám. 2, CZ-166 10 Prague 6, Czech Republic
- Department
of Physics, Tampere University of Technology, P.O. Box 692, FI-33101 Tampere, Finland
| |
Collapse
|