1
|
Qutbuddin Y, Guinart A, Gavrilović S, Al Nahas K, Feringa BL, Schwille P. Light-Activated Synthetic Rotary Motors in Lipid Membranes Induce Shape Changes Through Membrane Expansion. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2024; 36:e2311176. [PMID: 38215457 DOI: 10.1002/adma.202311176] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/25/2023] [Revised: 12/16/2023] [Indexed: 01/14/2024]
Abstract
Membranes are the key structures to separate and spatially organize cellular systems. Their rich dynamics and transformations during the cell cycle are orchestrated by specific membrane-targeted molecular machineries, many of which operate through energy dissipation. Likewise, man-made light-activated molecular rotary motors have previously shown drastic effects on cellular systems, but their physical roles on and within lipid membranes remain largely unexplored. Here, the impact of rotary motors on well-defined biological membranes is systematically investigated. Notably, dramatic mechanical transformations are observed in these systems upon motor irradiation, indicative of motor-induced membrane expansion. The influence of several factors on this phenomenon is systematically explored, such as motor concentration and membrane composition., Membrane fluidity is found to play a crucial role in motor-induced deformations, while only minor contributions from local heating and singlet oxygen generation are observed. Most remarkably, the membrane area expansion under the influence of the motors continues as long as irradiation is maintained, and the system stays out-of-equilibrium. Overall, this research contributes to a comprehensive understanding of molecular motors interacting with biological membranes, elucidating the multifaceted factors that govern membrane responses and shape transitions in the presence of these remarkable molecular machines, thereby supporting their future applications in chemical biology.
Collapse
Affiliation(s)
- Yusuf Qutbuddin
- Cellular and Molecular Biophysics, Max Planck Institute of Biochemistry, 82152, Martinsried, Germany
| | - Ainoa Guinart
- Stratingh Institute for Chemistry, University of Groningen, Groningen, 9747 AG, The Netherlands
| | - Svetozar Gavrilović
- Cellular and Molecular Biophysics, Max Planck Institute of Biochemistry, 82152, Martinsried, Germany
| | - Kareem Al Nahas
- Cellular and Molecular Biophysics, Max Planck Institute of Biochemistry, 82152, Martinsried, Germany
| | - Ben L Feringa
- Stratingh Institute for Chemistry, University of Groningen, Groningen, 9747 AG, The Netherlands
| | - Petra Schwille
- Cellular and Molecular Biophysics, Max Planck Institute of Biochemistry, 82152, Martinsried, Germany
| |
Collapse
|
2
|
Del Razo MJ, Crommelin D, Bolhuis PG. Data-driven dynamical coarse-graining for condensed matter systems. J Chem Phys 2024; 160:024108. [PMID: 38193550 DOI: 10.1063/5.0177553] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2023] [Accepted: 12/05/2023] [Indexed: 01/10/2024] Open
Abstract
Simulations of condensed matter systems often focus on the dynamics of a few distinguished components but require integrating the full system. A prime example is a molecular dynamics simulation of a (macro)molecule in a solution, where the molecule(s) and the solvent dynamics need to be integrated, rendering the simulations computationally costly and often unfeasible for physically/biologically relevant time scales. Standard coarse graining approaches can reproduce equilibrium distributions and structural features but do not properly include the dynamics. In this work, we develop a general data-driven coarse-graining methodology inspired by the Mori-Zwanzig formalism, which shows that macroscopic systems with a large number of degrees of freedom can be described by a few relevant variables and additional noise and memory terms. Our coarse-graining method consists of numerical integrators for the distinguished components, where the noise and interaction terms with other system components are substituted by a random variable sampled from a data-driven model. The model is parameterized using data from multiple short-time full-system simulations, and then, it is used to run long-time simulations. Applying our methodology to three systems-a distinguished particle under a harmonic and a bistable potential and a dimer with two metastable configurations-the resulting coarse-grained models are capable of reproducing not only the equilibrium distributions but also the dynamic behavior due to temporal correlations and memory effects. Remarkably, our method even reproduces the transition dynamics between metastable states, which is challenging to capture correctly. Our approach is not constrained to specific dynamics and can be extended to systems beyond Langevin dynamics, and, in principle, even to non-equilibrium dynamics.
Collapse
Affiliation(s)
- Mauricio J Del Razo
- Department of Mathematics and Computer Science, Freie Universität Berlin, Berlin, Germany
- Van't Hoff Institute for Molecular Sciences, University of Amsterdam, PO Box 94157, 1090GD Amsterdam, The Netherlands
- Korteweg-de Vries Institute for Mathematics, University of Amsterdam, PO Box 94248, 1090GD Amsterdam, The Netherlands
- Dutch Institute for Emergent Phenomena, University of Amsterdam, Amsterdam, The Netherlands
| | - Daan Crommelin
- Korteweg-de Vries Institute for Mathematics, University of Amsterdam, PO Box 94248, 1090GD Amsterdam, The Netherlands
- Centrum Wiskunde & Informatica, 1098 XG Amsterdam, The Netherlands
| | - Peter G Bolhuis
- Van't Hoff Institute for Molecular Sciences, University of Amsterdam, PO Box 94157, 1090GD Amsterdam, The Netherlands
| |
Collapse
|
3
|
Lado-Touriño I, Cerpa-Naranjo A. Coarse-Grained Molecular Dynamics of pH-Sensitive Lipids. Int J Mol Sci 2023; 24:ijms24054632. [PMID: 36902063 PMCID: PMC10003205 DOI: 10.3390/ijms24054632] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2023] [Revised: 02/20/2023] [Accepted: 02/25/2023] [Indexed: 03/06/2023] Open
Abstract
pH-sensitive lipids represent a class of lipids that can be protonated and destabilized in acidic environments, as they become positively charged in response to low-pH conditions. They can be incorporated into lipidic nanoparticles such as liposomes, which are able to change their properties and allow specific drug delivery at the acidic conditions encountered in some pathological microenvironments. In this work, we used coarse-grained molecular-dynamic simulations to study the stability of neutral and charged lipid bilayers containing POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and various kinds of ISUCA ((F)2-(imidazol-1-yl)succinic acid)-derived lipids, which can act as pH-sensitive molecules. In order to explore such systems, we used a MARTINI-derived forcefield, previously parameterized using all-atom simulation results. We calculated the average area per lipid, the second-rank order parameter and the lipid diffusion coefficient of both lipid bilayers made of pure components and mixtures of lipids in different proportions, under neutral or acidic conditions. The results show that the use of ISUCA-derived lipids disturbs the lipid bilayer structure, with the effect being particularly marked under acidic conditions. Although more-in depth studies on these systems must be carried out, these initial results are encouraging and the lipids designed in this research could be a good basis for developing new pH-sensitive liposomes.
Collapse
|
4
|
Kanekal KH, Rudzinski JF, Bereau T. Broad chemical transferability in structure-based coarse-graining. J Chem Phys 2022; 157:104102. [DOI: 10.1063/5.0104914] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Compared to top-down coarse-grained (CG) models, bottom-up approaches are capable of offering higher structural fidelity. This fidelity results from the tight link to a higher-resolution reference, making the CG model chemically specific. Unfortunately, chemical specificity can be at odds with compound-screening strategies, which call for transferable parametrizations. Here we present an approach to reconcile bottom-up, structure-preserving CG models with chemical transferability. We consider the bottom-up CG parametrization of 3,441 C7O2 small-molecule isomers. Our approach combines atomic representations, unsupervised learning, and a large-scale extended-ensemble force-matching parametrization. We first identify a subset of 19 representative molecules, which maximally encode the local environment of all gas-phase conformers. Reference interactions between the 19 representative molecules were obtained from both homogeneous bulk liquids and various binary mixtures. An extended-ensemble parametrization over all 703 state points leads to a CG model that is both structure-based and chemically transferable. Remarkably, the resulting force field is on average more structurally accurate than single-state-point equivalents. Averaging over the extended ensemble acts as a mean-force regularizer, smoothing out both force and structural correlations that are overly specific to a single state point. Our approach aims at transferability through a set of CG bead types that can be used to easily construct new molecules, while retaining the benefits of a structure-based parametrization.
Collapse
Affiliation(s)
- Kiran H. Kanekal
- AK Kremer - Theory Group, Max Planck Institute for Polymer Research, Germany
| | | | - Tristan Bereau
- Van 't Hoff Institute for Molecular Sciences and Informatics Institute, University of Amsterdam, Netherlands
| |
Collapse
|
5
|
Singhal A, Agur Sevink GJ. The role of size and nature in nanoparticle binding to a model lung membrane: an atomistic study. NANOSCALE ADVANCES 2021; 3:6635-6648. [PMID: 36132649 PMCID: PMC9417560 DOI: 10.1039/d1na00578b] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/24/2021] [Accepted: 09/16/2021] [Indexed: 05/09/2023]
Abstract
Understanding the uptake of nanoparticles (NPs) by different types of cellular membranes plays a pivotal role in the design of NPs for medical applications and in avoiding adverse effects that result in nanotoxicity. Yet, the role of key design parameters, such as the bare NP material, NP size and surface reactivity, and the nature of NP coatings, in membrane remodelling and uptake mechanisms is still very poorly understood, particularly towards the lower range of NP dimensions that are beyond the experimental imaging resolution. The same can be said about the role of a particular membrane composition. Here, we systematically employ biased and unbiased molecular dynamics simulations to calculate the binding energy for three bare materials (Ag/SiO2/TiO2) and three NP sizes (1/3/5 nm diameter) with a representative lung surfactant membrane, and to study their binding kinetics. The calculated binding energies show that irrespective of size, Ag nanoparticles bind very strongly to the bilayer, while the NPs made of SiO2 or TiO2 experience very low to no binding. The unbiased simulations provide insight into how the NPs and membrane affect each other in terms of the solvent-accessible surface area (SASA) of the NPs and the defect types and fluidity of the membrane. Using these systematic fine-grained results in coarsening procedures will pave the way for simulations considering NP sizes that are well beyond the membrane thickness, i.e. closer to experimental dimensions, for which different binding characteristics and more significant membrane remodelling are expected.
Collapse
Affiliation(s)
- Ankush Singhal
- Leiden Institute of Chemistry, Leiden University P.O. Box 9502 2300 RA Leiden The Netherlands
| | - G J Agur Sevink
- Leiden Institute of Chemistry, Leiden University P.O. Box 9502 2300 RA Leiden The Netherlands
| |
Collapse
|
6
|
Effect of packing density of lipid vesicles on the Aβ42 fibril polymorphism. Chem Phys Lipids 2021; 236:105073. [PMID: 33675780 DOI: 10.1016/j.chemphyslip.2021.105073] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2020] [Revised: 02/23/2021] [Accepted: 03/01/2021] [Indexed: 11/22/2022]
Abstract
The aggregation of amyloid-β 1-42 (Aβ42) on lipid membranes is closely related to the pathology of Alzheimer's disease (AD). Herein, we demonstrated the effect of the packing density of lipid vesicles on the Aβ42 fibrillation kinetics and fibril morphology. We used three distinct phosphatidylcholine (PC) lipids, containing different numbers of cis-double bonds in acyl chains, and therefore, a different packing density in the lipid vesicles. Our results showed that the fibrillation of Aβ42 was greatly enhanced and the formed fibrils became shorter as the number of double bonds in lipids increased. Due to the low-density characteristics of dioleoyl phosphatidylcholine (DOPC), Aβ42 monomers were able to interact with the hydrophobic acyl chain of lipids exposed to the aqueous phase, thereby inducing rapid fibrillation and short fibril morphologies. Furthermore, the effects of the anionic lipids dioleoyl phosphatidylserine (DOPS) and dioleoyl phosphatidylglycerol (DOPG), and mixed vesicles of DOPC/DOPS and DOPC/DOPG on Aβ42 fibrillations were investigated. The tight binding of Aβ42 to the lipid head groups via electrostatic interactions was able to suppress the modulation of Aβ42 fibrillations compared to accelerated fibrillations on loosely packed membranes. Our proposed mechanism regarding the influence of lipid packing density on Aβ42 fibrillations provides an advanced understanding of lipid-associated amyloid fibrillations.
Collapse
|
7
|
Rudzinski JF, Bereau T. Coarse-grained conformational surface hopping: Methodology and transferability. J Chem Phys 2020; 153:214110. [DOI: 10.1063/5.0031249] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
| | - Tristan Bereau
- Max Planck Institute for Polymer Research, 55128 Mainz, Germany
- Van ’t Hoff Institute for Molecular Sciences and Informatics Institute, University of Amsterdam, Amsterdam 1098 XH, The Netherlands
| |
Collapse
|
8
|
Aceves-Luna H, Glossman-Mitnik D, Flores-Holguín N. Oxidation degree of a cell membrane model and its response to structural changes, a coarse-grained molecular dynamics approach. J Biomol Struct Dyn 2020; 40:1930-1941. [PMID: 33063644 DOI: 10.1080/07391102.2020.1833759] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Oxidative stress plays an essential role in the regulation of vital processes in living organisms. Reactive oxygen species can react chemically with the constituents of the cells leading to irreversible damage. The first structure of the cell in contact with the environment that surrounds it is the membrane, which protects it and allows the exchange of substances. Some signals manifest when the components of a bilayer are undergoing oxidation, like an increase in the lipid area, decrease in the thickness of the bilayer, and exchange of the oxidized groups toward the bilayer surface. In this investigation, a molecular dynamics simulation was done on a set of Dioleoylphosphatidylcholine membranes with different percentage of oxidized lipids, in order to observe the effect of the oxidation degree on the membrane structure. It was found that, as higher the concentration of oxidized lipids is, the larger the damage of the membrane. This is reflected in the increase in the lipid area and the decrease in the thickness and membrane packing. Also, it was observed that hydrophobicity inside the membrane decreases as the oxidation percentage increases.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Hugo Aceves-Luna
- Laboratorio Virtual NANOCOSMOS, Departamento de Medio Ambiente y Energía, Centro de Investigación en Materiales Avanzados, Complejo Industrial Chihuahua, Chihuahua, Chih, Mexico
| | - Daniel Glossman-Mitnik
- Laboratorio Virtual NANOCOSMOS, Departamento de Medio Ambiente y Energía, Centro de Investigación en Materiales Avanzados, Complejo Industrial Chihuahua, Chihuahua, Chih, Mexico
| | - Norma Flores-Holguín
- Laboratorio Virtual NANOCOSMOS, Departamento de Medio Ambiente y Energía, Centro de Investigación en Materiales Avanzados, Complejo Industrial Chihuahua, Chihuahua, Chih, Mexico
| |
Collapse
|
9
|
Angelescu DG. Structural behavior of amphiphilic polyion complexes interacting with saturated lipid membranes investigated by coarse-grained molecular dynamic simulations. RSC Adv 2020; 10:39204-39216. [PMID: 35518426 PMCID: PMC9057367 DOI: 10.1039/d0ra06894b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2020] [Accepted: 10/16/2020] [Indexed: 11/21/2022] Open
Abstract
Neutral polyelectrolyte complexes (PECs) made from an amphiphilic multiblock copolymer of type (AnBn)m and an oppositely charged polyion and interacting with a dipalmitoylphosphatidylcholine (DPPC) lipid membrane.
Collapse
Affiliation(s)
- Daniel G. Angelescu
- Romanian Academy
- “Ilie Murgulescu” Institute of Physical Chemistry
- 060021 Bucharest
- Romania
| |
Collapse
|
10
|
Terzi MM, Ergüder MF, Deserno M. A consistent quadratic curvature-tilt theory for fluid lipid membranes. J Chem Phys 2019; 151:164108. [PMID: 31675861 DOI: 10.1063/1.5119683] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
The tilt of a lipid molecule describes the deviation of its orientation away from the local normal of its embedding membrane. Tilt is the subleading degree of freedom after a membrane's geometry, and it becomes relevant at scales comparable to lipid bilayer thickness. Building on earlier work by Hamm and Kozlov [Eur. Phys. J. E 3, 323 (2000)], who envisioned lipid membranes as thin prestressed fluid elastic films, and Terzi and Deserno [J. Chem. Phys. 147, 084702 (2017)], who discovered a new coupling term between splay and tilt divergence, we construct a theory of membrane elasticity that is quadratic in geometry and tilt and complete at order 1/length2. We show that a general and consistent treatment of both lateral and transverse depth-dependent shear stresses creates several contributions to the elastic energy density, of which only a subset had previously been identified. Apart from the well-known penalty of lipid twist (the curl of tilt), these terms generate no qualitatively new phenomenology, but they quantitatively revise the connections between the moduli of a tilt-curvature theory and its underlying microscopic foundation. In particular, we argue that the monolayer Gaussian curvature modulus κ¯m, widely believed to be equal to the second moment of the transmonolayer stress profile, acquires a second contribution from lipid twist, which is always negative. This could resolve the long-standing conundrum that many measured values of κ¯m appeared to have a sign that violates basic stability considerations. We also show that the previously discovered novel coupling between splay and tilt divergence is not simply proportional to κ¯m but acquires its own splay-tilt coupling modulus, κst,m. We explore the predictions of our theory for various elastic moduli and their mutual interrelations and use an extensive set of existing atomistic molecular dynamics simulations for 12 different lipid types to collectively reason about such predictions. We find that bending rigidities are captured fairly well by existing theories, while reliable predictions for local moduli, especially the splay-tilt coupling modulus, remain challenging.
Collapse
Affiliation(s)
- M Mert Terzi
- Department of Physics, Carnegie Mellon University, 5000 Forbes Ave., Pittsburgh, Pennsylvania 15213, USA
| | - Muhammed F Ergüder
- Department of Physics, Carnegie Mellon University, 5000 Forbes Ave., Pittsburgh, Pennsylvania 15213, USA
| | - Markus Deserno
- Department of Physics, Carnegie Mellon University, 5000 Forbes Ave., Pittsburgh, Pennsylvania 15213, USA
| |
Collapse
|
11
|
Lebold KM, Noid WG. Dual-potential approach for coarse-grained implicit solvent models with accurate, internally consistent energetics and predictive transferability. J Chem Phys 2019; 151:164113. [PMID: 31675902 DOI: 10.1063/1.5125246] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Abstract
The dual-potential approach promises coarse-grained (CG) models that accurately reproduce both structural and energetic properties, while simultaneously providing predictive estimates for the temperature-dependence of the effective CG potentials. In this work, we examine the dual-potential approach for implicit solvent CG models that reflect large entropic effects from the eliminated solvent. Specifically, we construct implicit solvent models at various resolutions, R, by retaining a fraction 0.10 ≤ R ≤ 0.95 of the molecules from a simple fluid of Lennard-Jones spheres. We consider the dual-potential approach in both the constant volume and constant pressure ensembles across a relatively wide range of temperatures. We approximate the many-body potential of mean force for the remaining solutes with pair and volume potentials, which we determine via multiscale coarse-graining and self-consistent pressure-matching, respectively. Interestingly, with increasing temperature, the pair potentials appear increasingly attractive, while the volume potentials become increasingly repulsive. The dual-potential approach not only reproduces the atomic energetics but also quite accurately predicts this temperature-dependence. We also derive an exact relationship between the thermodynamic specific heat of an atomic model and the energetic fluctuations that are observable at the CG resolution. With this generalized fluctuation relationship, the approximate CG models quite accurately reproduce the thermodynamic specific heat of the underlying atomic model.
Collapse
Affiliation(s)
- Kathryn M Lebold
- Department of Chemistry, Penn State University, University Park, Pennsylvania 16802, USA
| | - W G Noid
- Department of Chemistry, Penn State University, University Park, Pennsylvania 16802, USA
| |
Collapse
|
12
|
Barrera EE, Machado MR, Pantano S. Fat SIRAH: Coarse-Grained Phospholipids To Explore Membrane-Protein Dynamics. J Chem Theory Comput 2019; 15:5674-5688. [PMID: 31433946 DOI: 10.1021/acs.jctc.9b00435] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
The capability to handle highly heterogeneous molecular assemblies in a consistent manner is among the greatest challenges faced when deriving simulation parameters. This is particularly the case for coarse-grained (CG) simulations in which chemical functional groups are lumped into effective interaction centers for which transferability between different chemical environments is not guaranteed. Here, we introduce the parametrization of a set of CG phospholipids compatible with the latest version of the SIRAH force field for proteins. The newly introduced lipid species include different acylic chain lengths and partial unsaturation, as well as polar and acidic head groups that show a very good reproduction of structural membrane determinants, such as areas per lipid, thickness, order parameter, etc., and their dependence with temperature. Simulation of membrane proteins showed unprecedented accuracy in the unbiased description of the thickness-dependent membrane-protein orientation in systems where this information is experimentally available (namely, the SarcoEndoplasmic Reticulum Calcium-SERCA-pump and its regulator Phospholamban). The interactions that lead to this faithful reproduction can be traced down to the single amino acid-lipid interaction level and show full agreement with biochemical data present in the literature. Finally, the present parametrization is implemented in the GROMACS and AMBER simulation packages facilitating its use by a wide portion of the biocomputing community.
Collapse
Affiliation(s)
- Exequiel E Barrera
- Biomolecular Simulations Group , Institut Pasteur de Montevideo , Mataojo 2020 , CP 11400 Montevideo , Uruguay
| | - Matías R Machado
- Biomolecular Simulations Group , Institut Pasteur de Montevideo , Mataojo 2020 , CP 11400 Montevideo , Uruguay
| | - Sergio Pantano
- Biomolecular Simulations Group , Institut Pasteur de Montevideo , Mataojo 2020 , CP 11400 Montevideo , Uruguay.,Shanghai Institute for Advanced Immunochemical Studies , ShanghaiTech University , Shanghai 201210 , China
| |
Collapse
|
13
|
Chen H, Panagiotopoulos AZ. Molecular Modeling of Surfactant Micellization Using Solvent-Accessible Surface Area. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2019; 35:2443-2450. [PMID: 30624073 DOI: 10.1021/acs.langmuir.8b03440] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
We report a new implicit solvent simulation model for studying the self-assembly of surfactants, where the hydrophobic interactions were captured by calculating the relative changes of the solvent-accessible surface area (SASA) of the hydrophobic domains. Using histogram-reweighting grand canonical Monte Carlo simulations, we demonstrate that this approach allows us to match both the experimental critical micelle concentrations (cmc) and micellar aggregation numbers simultaneously with a single phenomenological surface tension γSASA for the poly(oxyethylene) monoalkyl ether (C mE n) surfactants in aqueous solutions. Excellent transferability is observed: the same model can accurately predict the experimental cmc and aggregation numbers for the C mE n surfactants with the alkyl lengths m between 6 and 12 and the poly(oxyethylene) lengths n between 1 and 9. The SASA-based implicit solvent model put forward in this work is general and may be applied to study more complex amphiphilic systems such as surfactants with branched alkyl chains or surfactant-hydrocarbon mixtures.
Collapse
Affiliation(s)
- Hsieh Chen
- Aramco Services Company: Aramco Research Center-Boston , 400 Technology Square , Cambridge , Massachusetts 02139 , United States
| | - Athanassios Z Panagiotopoulos
- Department of Chemical and Biological Engineering , Princeton University , Princeton , New Jersey 08544 , United States
| |
Collapse
|
14
|
Angelescu DG. Coarse-grained simulation studies on the adsorption of polyelectrolyte complexes upon lipid membranes. Phys Chem Chem Phys 2019; 21:12446-12459. [DOI: 10.1039/c9cp01448a] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Conformations of a polyelectrolyte complex irreversibly bound to a zwitterionic lipid bilayer.
Collapse
Affiliation(s)
- Daniel G. Angelescu
- Romanian Academy
- “Ilie Murgulescu” Institute of Physical Chemistry
- 060021 Bucharest
- Romania
| |
Collapse
|
15
|
Kolli HB, de Nicola A, Bore SL, Schäfer K, Diezemann G, Gauss J, Kawakatsu T, Lu ZY, Zhu YL, Milano G, Cascella M. Hybrid Particle-Field Molecular Dynamics Simulations of Charged Amphiphiles in an Aqueous Environment. J Chem Theory Comput 2018; 14:4928-4937. [DOI: 10.1021/acs.jctc.8b00466] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Hima Bindu Kolli
- Department of Chemistry and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, P.O.
Box 1033, Blindern, 0315 Oslo, Norway
| | - Antonio de Nicola
- Department of Organic Materials Science, Yamagata University, 4-3-16 Jonan Yonezawa, Yamagata-ken 992-8510, Japan
| | - Sigbjørn Løland Bore
- Department of Chemistry and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, P.O.
Box 1033, Blindern, 0315 Oslo, Norway
| | - Ken Schäfer
- Institut für
Physikalische Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, 55128 Mainz, Germany
| | - Gregor Diezemann
- Institut für
Physikalische Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, 55128 Mainz, Germany
| | - Jürgen Gauss
- Institut für
Physikalische Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, 55128 Mainz, Germany
| | - Toshihiro Kawakatsu
- Department of Physics, Tohoku University, Aoba, Aramaki, Aoba-ku, Sendai 980-8578, Japan
| | - Zhong-Yuan Lu
- State Key Laboratory of Supramolecular Structure and Materials, Institute of Theoretical Chemistry, Jilin University, Changchun 130023, China
| | - You-Liang Zhu
- State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, China
| | - Giuseppe Milano
- Department of Organic Materials Science, Yamagata University, 4-3-16 Jonan Yonezawa, Yamagata-ken 992-8510, Japan
| | - Michele Cascella
- Department of Chemistry and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, P.O.
Box 1033, Blindern, 0315 Oslo, Norway
| |
Collapse
|
16
|
Abstract
On mesoscopic scales, lipid membranes are well described by continuum theories whose main ingredients are the curvature of a membrane's reference surface and the tilt of its lipid constituents. In particular, Hamm and Kozlov [Eur. Phys. J. E 3, 323 (2000)] have shown how to systematically derive such a tilt-curvature Hamiltonian based on the elementary assumption of a thin fluid elastic sheet experiencing internal lateral pre-stress. Performing a dimensional reduction, they not only derive the basic form of the effective surface Hamiltonian but also express its emergent elastic couplings as trans-membrane moments of lower-level material parameters. In the present paper, we argue, though, that their derivation unfortunately missed a coupling term between curvature and tilt. This term arises because, as one moves along the membrane, the curvature-induced change of transverse distances contributes to the area strain-an effect that was believed to be small but nevertheless ends up contributing at the same (quadratic) order as all other terms in their Hamiltonian. We illustrate the consequences of this amendment by deriving the monolayer and bilayer Euler-Lagrange equations for the tilt, as well as the power spectra of shape, tilt, and director fluctuations. A particularly curious aspect of our new term is that its associated coupling constant is the second moment of the lipid monolayer's lateral stress profile-which within this framework is equal to the monolayer Gaussian curvature modulus, κ¯m. On the one hand, this implies that many theoretical predictions now contain a parameter that is poorly known (because the Gauss-Bonnet theorem limits access to the integrated Gaussian curvature); on the other hand, the appearance of κ¯m outside of its Gaussian curvature provenance opens opportunities for measuring it by more conventional means, for instance by monitoring a membrane's undulation spectrum at short scales.
Collapse
Affiliation(s)
- M Mert Terzi
- Department of Physics, Carnegie Mellon University, 5000 Forbes Ave., Pittsburgh, Pennsylvania 15213, USA
| | - Markus Deserno
- Department of Physics, Carnegie Mellon University, 5000 Forbes Ave., Pittsburgh, Pennsylvania 15213, USA
| |
Collapse
|
17
|
Wan M, Gao L, Fang W. Implicit-solvent dissipative particle dynamics force field based on a four-to-one coarse-grained mapping scheme. PLoS One 2018; 13:e0198049. [PMID: 29795682 PMCID: PMC5967728 DOI: 10.1371/journal.pone.0198049] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2018] [Accepted: 05/12/2018] [Indexed: 11/30/2022] Open
Abstract
A new set of efficient solvent-free dissipative particle dynamics (DPD) force fields was developed for phospholipids and peptides. To enhance transferability, this model maps around four heavy atoms and their connected hydrogen atoms into a coarse-grained elementary bead based on functional group. The effective hybrid potential between any pair of beads is composed of a short-range repulsive soft-core potential that directly adopts the form of an explicit-solvent DPD model and a long-range attractive hydrophobic potential. The parameters of the attractive potentials for lipid molecules were obtained by fitting the explicit-solvent DPD simulation of one bead of any type in a water box, then finely tuning it until the bilayer membrane properties obtained in the explicit-solvent model were matched. These parameters were further extended to amino acids according to bead type. The structural and elastic properties of bilayer membranes, free energy profiles for a lipid flip-flop and amino acid analogues translocating across the membrane, and membrane pore formation induced by antimicrobial peptides obtained from this solvent-free DPD force field considerably agreed with the explicit-solvent DPD results. Importantly, the efficiency of this method is guaranteed to accelerate the assembly of vesicles composed of several thousand lipids by up to 50-fold, rendering the experimental liposome dynamics as well as membrane-peptide interactions feasible at accessible computational expense.
Collapse
Affiliation(s)
- Mingwei Wan
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing, China
| | - Lianghui Gao
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing, China
- * E-mail:
| | - Weihai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing, China
| |
Collapse
|
18
|
Langenberg M, Jackson NE, de Pablo JJ, Müller M. Role of translational entropy in spatially inhomogeneous, coarse-grained models. J Chem Phys 2018. [DOI: 10.1063/1.5018178] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Marcel Langenberg
- Institut für Theoretische Physik, Universität Göttingen, 37077 Göttingen, Germany
| | - Nicholas E. Jackson
- Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60615, USA
- Argonne National Laboratory, Lemont, Illinois 06349, USA
| | - Juan J. de Pablo
- Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60615, USA
- Argonne National Laboratory, Lemont, Illinois 06349, USA
| | - Marcus Müller
- Institut für Theoretische Physik, Universität Göttingen, 37077 Göttingen, Germany
| |
Collapse
|
19
|
Matthews JR, Shirazinejad CR, Isakson GA, Demers SME, Hafner JH. Structural Analysis by Enhanced Raman Scattering. NANO LETTERS 2017; 17:2172-2177. [PMID: 28166410 DOI: 10.1021/acs.nanolett.6b04509] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/25/2023]
Abstract
Gold nanostructures focus light to a molecular length scale at their surface, creating the possibility to visualize molecular structure. The high optical intensity leads to surface enhanced Raman scattering (SERS) from nearby molecules. SERS spectra contain information on molecular position and orientation relative to the surface but are difficult to interpret quantitatively. Here we describe a ratiometric analysis method that combines SERS and unenhanced Raman spectra with theoretical calculations of the optical field and molecular polarizability. When applied to the surfactant layer on gold nanorods, the alkane chain is found to be tilted 25° to the surface normal, which matches previous reports of the layer thickness. The analysis was also applied to fluid phase phospholipid bilayers that contain tryptophan on the surface of gold nanorods. The lipid double bond was found to be oriented normal to the bilayer and 13 Å from the nitrogen atom. Tryptophan was found to sit near the glycerol headgroup region with its indole ring 43° from the bilayer normal. This new method can determine specific interfacial structure under ambient conditions, with microscopic quantities of material, and without molecular labels.
Collapse
Affiliation(s)
- James R Matthews
- Department of Physics & Astronomy and ‡Department of Chemistry, Rice University , Houston, Texas 77251, United States
| | - Cyna R Shirazinejad
- Department of Physics & Astronomy and ‡Department of Chemistry, Rice University , Houston, Texas 77251, United States
| | - Grace A Isakson
- Department of Physics & Astronomy and ‡Department of Chemistry, Rice University , Houston, Texas 77251, United States
| | - Steven M E Demers
- Department of Physics & Astronomy and ‡Department of Chemistry, Rice University , Houston, Texas 77251, United States
| | - Jason H Hafner
- Department of Physics & Astronomy and ‡Department of Chemistry, Rice University , Houston, Texas 77251, United States
| |
Collapse
|
20
|
Sevink GJA, Schmid F, Kawakatsu T, Milano G. Combining cell-based hydrodynamics with hybrid particle-field simulations: efficient and realistic simulation of structuring dynamics. SOFT MATTER 2017; 13:1594-1623. [PMID: 28128838 DOI: 10.1039/c6sm02252a] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
We have extended an existing hybrid MD-SCF simulation technique that employs a coarsening step to enhance the computational efficiency of evaluating non-bonded particle interactions. This technique is conceptually equivalent to the single chain in mean-field (SCMF) method in polymer physics, in the sense that non-bonded interactions are derived from the non-ideal chemical potential in self-consistent field (SCF) theory, after a particle-to-field projection. In contrast to SCMF, however, MD-SCF evolves particle coordinates by the usual Newton's equation of motion. Since collisions are seriously affected by the softening of non-bonded interactions that originates from their evaluation at the coarser continuum level, we have devised a way to reinsert the effect of collisions on the structural evolution. Merging MD-SCF with multi-particle collision dynamics (MPCD), we mimic particle collisions at the level of computational cells and at the same time properly account for the momentum transfer that is important for a realistic system evolution. The resulting hybrid MD-SCF/MPCD method was validated for a particular coarse-grained model of phospholipids in aqueous solution, against reference full-particle simulations and the original MD-SCF model. We additionally implemented and tested an alternative and more isotropic finite difference gradient. Our results show that efficiency is improved by merging MD-SCF with MPCD, as properly accounting for hydrodynamic interactions considerably speeds up the phase separation dynamics, with negligible additional computational costs compared to efficient MD-SCF. This new method enables realistic simulations of large-scale systems that are needed to investigate the applications of self-assembled structures of lipids in nanotechnologies.
Collapse
Affiliation(s)
- G J A Sevink
- Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands.
| | - F Schmid
- Institut for Physik, Johannes Gutenberg Univeristat Mainz, Staudingerweg 7-9, 55128 Mainz, Germany
| | - T Kawakatsu
- Department of Physics, Tohoku University, Aoba, Aoba-ku, Sendai 980-8578, Japan
| | - G Milano
- Dipartimento di Chimica e Biologia, Universit degli Studi di Salerno, via Ponte don Melillo, Fisciano, Italy
| |
Collapse
|
21
|
Ganesan SJ, Xu H, Matysiak S. Influence of Monovalent Cation Size on Nanodomain Formation in Anionic–Zwitterionic Mixed Bilayers. J Phys Chem B 2017; 121:787-799. [DOI: 10.1021/acs.jpcb.6b10099] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Sai J. Ganesan
- Fischell
Department of Bioengineering and ‡Biophysics Program, Institute of
Physical Science and Technology, University of Maryland, College
Park, Maryland 20742, United States
| | - Hongcheng Xu
- Fischell
Department of Bioengineering and ‡Biophysics Program, Institute of
Physical Science and Technology, University of Maryland, College
Park, Maryland 20742, United States
| | - Silvina Matysiak
- Fischell
Department of Bioengineering and ‡Biophysics Program, Institute of
Physical Science and Technology, University of Maryland, College
Park, Maryland 20742, United States
| |
Collapse
|
22
|
Bereau T, Bennett WFD, Pfaendtner J, Deserno M, Karttunen M. Folding and insertion thermodynamics of the transmembrane WALP peptide. J Chem Phys 2016; 143:243127. [PMID: 26723612 DOI: 10.1063/1.4935487] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
The anchor of most integral membrane proteins consists of one or several helices spanning the lipid bilayer. The WALP peptide, GWW(LA)n (L)WWA, is a common model helix to study the fundamentals of protein insertion and folding, as well as helix-helix association in the membrane. Its structural properties have been illuminated in a large number of experimental and simulation studies. In this combined coarse-grained and atomistic simulation study, we probe the thermodynamics of a single WALP peptide, focusing on both the insertion across the water-membrane interface, as well as folding in both water and a membrane. The potential of mean force characterizing the peptide's insertion into the membrane shows qualitatively similar behavior across peptides and three force fields. However, the Martini force field exhibits a pronounced secondary minimum for an adsorbed interfacial state, which may even become the global minimum-in contrast to both atomistic simulations and the alternative PLUM force field. Even though the two coarse-grained models reproduce the free energy of insertion of individual amino acids side chains, they both underestimate its corresponding value for the full peptide (as compared with atomistic simulations), hinting at cooperative physics beyond the residue level. Folding of WALP in the two environments indicates the helix as the most stable structure, though with different relative stabilities and chain-length dependence.
Collapse
Affiliation(s)
- Tristan Bereau
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - W F Drew Bennett
- Department of Chemistry, University of Waterloo, 200 University Avenue West, Waterloo, Ontario N2L 3G1, Canada
| | - Jim Pfaendtner
- Department of Chemical Engineering, University of Washington, Seattle, Washington 98195, USA
| | - Markus Deserno
- Department of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
| | - Mikko Karttunen
- Department of Mathematics and Computer Science & Institute for Complex Molecular Systems, Eindhoven University of Technology, P.O. Box 513, MetaForum, 5600 MB Eindhoven, The Netherlands
| |
Collapse
|
23
|
Wang X, Deserno M. Determining the Lipid Tilt Modulus by Simulating Membrane Buckles. J Phys Chem B 2016; 120:6061-73. [DOI: 10.1021/acs.jpcb.6b02016] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Affiliation(s)
- Xin Wang
- Department of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Markus Deserno
- Department of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| |
Collapse
|
24
|
Ganesan SJ, Xu H, Matysiak S. Effect of lipid head group interactions on membrane properties and membrane-induced cationic β-hairpin folding. Phys Chem Chem Phys 2016; 18:17836-50. [DOI: 10.1039/c5cp07669b] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Stages in membrane induced SVS-1 folding.
Collapse
Affiliation(s)
- Sai J. Ganesan
- Fischell Department of Bioengineering
- University of Maryland
- College Park
- USA
| | - Hongcheng Xu
- Biophysics Program
- University of Maryland
- College Park
- USA
| | - Silvina Matysiak
- Fischell Department of Bioengineering
- University of Maryland
- College Park
- USA
- Biophysics Program
| |
Collapse
|
25
|
Effect of Cholesterol on the Structure of a Five-Component Mitochondria-Like Phospholipid Membrane. MEMBRANES 2015; 5:664-84. [PMID: 26529029 PMCID: PMC4704005 DOI: 10.3390/membranes5040664] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/05/2015] [Accepted: 10/16/2015] [Indexed: 12/23/2022]
Abstract
Cellular membranes have a complex phospholipid composition that varies greatly depending on the organism, cell type and function. In spite of this complexity, most structural data available for phospholipid bilayers concern model systems containing only one or two different phospholipids. Here, we examine the effect of cholesterol on the structure of a complex membrane reflecting the lipid composition of mitochondrial membranes, with five different types of headgroups (phosphatidylcholine (PC), phosphatidylethanolamine (PE), phosphatidylinositol (PI), phosphatidylserine (PS) and cardiolipin (CL)) and a variety of hydrocarbon tails. This particular system was chosen because elevated cholesterol contents in mitochondrial membranes have been linked to a breaking down of Bax-mediated membrane permeabilization and resistance to cancer treatments. High resolution electron density profiles were determined by X-ray reflectivity, while the area per phospholipid chain, Apc, and the chain order parameter, SX-ray, were determined by wide-angle X-ray scattering (WAXS). We show that chain order increases upon the addition of cholesterol, resulting in both a thickening of the lipid bilayer and a reduction in the average surface area per phospholipid chain. This effect, well known as cholesterol’s condensation effect, is similar, but not as pronounced as for single-component phospholipid membranes. We conclude by discussing the relevance of these findings for the insertion of the pro-apoptotic protein Bax in mitochondrial membranes with elevated cholesterol content.
Collapse
|
26
|
Davoudi S, Amjad‐Iranagh S, Zaeifi Yamchi M. Molecular dynamic simulation of Ca
2+
‐ATPase interacting with lipid bilayer membrane. IET Nanobiotechnol 2015; 9:85-94. [DOI: 10.1049/iet-nbt.2013.0073] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Affiliation(s)
- Samaneh Davoudi
- Chemical Engineering DepartmentAmirkabir University of TechnologyTehranIran
| | | | | |
Collapse
|
27
|
|
28
|
Deserno M. Fluid lipid membranes: From differential geometry to curvature stresses. Chem Phys Lipids 2015; 185:11-45. [DOI: 10.1016/j.chemphyslip.2014.05.001] [Citation(s) in RCA: 137] [Impact Index Per Article: 15.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2014] [Revised: 04/21/2014] [Accepted: 05/06/2014] [Indexed: 10/25/2022]
|
29
|
Enhanced Sampling of Coarse-Grained Transmembrane-Peptide Structure Formation from Hydrogen-Bond Replica Exchange. J Membr Biol 2014; 248:395-405. [DOI: 10.1007/s00232-014-9738-9] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2014] [Accepted: 09/27/2014] [Indexed: 12/14/2022]
|
30
|
Sevink GJA, Fraaije JGEM. Efficient solvent-free dissipative particle dynamics for lipid bilayers. SOFT MATTER 2014; 10:5129-5146. [PMID: 24909682 DOI: 10.1039/c4sm00297k] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
We rigorously derived effective potentials for solvent-free DPD simulation of lipid bilayers. The derivation relies on an earlier developed hybrid particle/field method and is based on the idea that the solvent is always in local equilibrium on a coarse-grained time scale, given the instantaneous templates set by the self-assembly structure. By relating the parameters in the effective implicit-solvent potentials directly to the lipid-solvent interactions and membrane properties for the explicit solvent DPD model, we constitute an efficient and general procedure for reformulating any DPD membrane model in an implicit-solvent form. Here, we determined these membrane properties for two existing DPD models, via an analysis of membrane fluctuation spectra. Equivalent single-processor implicit- and explicit-solvent calculations show the trade-mark of implicit solvent simulation: a 20-fold reduction of the total simulation time for a system containing 92% solvent. This increased efficiency enabled us to realistically simulate the spontaneous formation of a ∼20 nm diameter vesicle on a single processor overnight. We believe that this work will contribute to an enhanced computational study of large vesicles and thus a better understanding of experimental liposome dynamics.
Collapse
Affiliation(s)
- G J A Sevink
- Leiden Institute of Chemistry, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | | |
Collapse
|
31
|
Aydin F, Dutt M. Bioinspired Vesicles Encompassing Two-Tail Phospholipids: Self-Assembly and Phase Segregation via Implicit Solvent Coarse-Grained Molecular Dynamics. J Phys Chem B 2014; 118:8614-23. [DOI: 10.1021/jp503376r] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Affiliation(s)
- Fikret Aydin
- Department
of Chemical and Biochemical Engineering, Rutgers The State University of New Jersey, Piscataway, New Jersey 08854, United States
| | - Meenakshi Dutt
- Department
of Chemical and Biochemical Engineering, Rutgers The State University of New Jersey, Piscataway, New Jersey 08854, United States
| |
Collapse
|
32
|
Kar P, Gopal SM, Cheng YM, Panahi A, Feig M. Transferring the PRIMO Coarse-Grained Force Field to the Membrane Environment: Simulations of Membrane Proteins and Helix-Helix Association. J Chem Theory Comput 2014; 10:3459-3472. [PMID: 25136271 PMCID: PMC4132866 DOI: 10.1021/ct500443v] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2014] [Indexed: 12/03/2022]
Abstract
![]()
An
extension of the recently developed PRIMO coarse-grained force
field to membrane environments, PRIMO-M, is described. The membrane
environment is modeled with the heterogeneous dielectric generalized
Born (HDGB) methodology that simply replaces the standard generalized
Born model in PRIMO without further parametrization. The resulting
model was validated by comparing amino acid insertion free energy
profiles and application in molecular dynamics simulations of membrane
proteins and membrane-interacting peptides. Membrane proteins with
148–661 amino acids show stable root-mean-squared-deviations
(RMSD) between 2 and 4 Å for most systems. Transmembrane helical
peptides maintain helical shape and exhibit tilt angles in good agreement
with experimental or other simulation data. The association of two
glycophorin A (GpA) helices was simulated using replica exchange molecular
dynamics simulations yielding the correct dimer structure with a crossing
angle in agreement with previous studies. Finally, conformational
sampling of the influenza fusion peptide also generates structures
in agreement with previous studies. Overall, these findings suggest
that PRIMO-M can be used to study membrane bound peptides and proteins
and validates the transferable nature of the PRIMO coarse-grained
force field.
Collapse
Affiliation(s)
- Parimal Kar
- Department of Biochemistry and Molecular Biology and Department of Chemistry, Michigan State University , East Lansing, Michigan 48824, United States
| | - Srinivasa Murthy Gopal
- Department of Biochemistry and Molecular Biology and Department of Chemistry, Michigan State University , East Lansing, Michigan 48824, United States
| | - Yi-Ming Cheng
- Department of Biochemistry and Molecular Biology and Department of Chemistry, Michigan State University , East Lansing, Michigan 48824, United States
| | - Afra Panahi
- Departments of Chemistry and Biophysics, University of Michigan , Ann Arbor, Michigan 48109, United States
| | - Michael Feig
- Department of Biochemistry and Molecular Biology and Department of Chemistry, Michigan State University , East Lansing, Michigan 48824, United States ; Department of Biochemistry and Molecular Biology and Department of Chemistry, Michigan State University , East Lansing, Michigan 48824, United States
| |
Collapse
|
33
|
Bereau T, Wang ZJ, Deserno M. More than the sum of its parts: coarse-grained peptide-lipid interactions from a simple cross-parametrization. J Chem Phys 2014; 140:115101. [PMID: 24655203 PMCID: PMC3977883 DOI: 10.1063/1.4867465] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2013] [Accepted: 02/04/2014] [Indexed: 11/14/2022] Open
Abstract
Interfacial systems are at the core of fascinating phenomena in many disciplines, such as biochemistry, soft-matter physics, and food science. However, the parametrization of accurate, reliable, and consistent coarse-grained (CG) models for systems at interfaces remains a challenging endeavor. In the present work, we explore to what extent two independently developed solvent-free CG models of peptides and lipids--of different mapping schemes, parametrization methods, target functions, and validation criteria--can be combined by only tuning the cross-interactions. Our results show that the cross-parametrization can reproduce a number of structural properties of membrane peptides (for example, tilt and hydrophobic mismatch), in agreement with existing peptide-lipid CG force fields. We find encouraging results for two challenging biophysical problems: (i) membrane pore formation mediated by the cooperative action of several antimicrobial peptides, and (ii) the insertion and folding of the helix-forming peptide WALP23 in the membrane.
Collapse
Affiliation(s)
- Tristan Bereau
- Department of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
| | - Zun-Jing Wang
- Department of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
| | - Markus Deserno
- Department of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
| |
Collapse
|
34
|
Hu M, de Jong DH, Marrink SJ, Deserno M. Gaussian curvature elasticity determined from global shape transformations and local stress distributions: a comparative study using the MARTINI model. Faraday Discuss 2013; 161:365-82; discussion 419-59. [PMID: 23805750 DOI: 10.1039/c2fd20087b] [Citation(s) in RCA: 77] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We calculate the Gaussian curvature modulus kappa of a systematically coarse-grained (CG) one-component lipid membrane by applying the method recently proposed by Hu et al. [Biophys. J., 2012, 102, 1403] to the MARTINI representation of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC). We find the value kappa/kappa = -1.04 +/- 0.03 for the elastic ratio between the Gaussian and the mean curvature modulus and deduce kappa(m)/kappa(m) = -0.98 +/- 0.09 for the monolayer elastic ratio, where the latter is based on plausible assumptions for the distance z0 of the monolayer neutral surface from the bilayer midplane and the spontaneous lipid curvature K(0m). By also analyzing the lateral stress profile sigma0(z) of our system, two other lipid types and pertinent data from the literature, we show that determining K(0m) and kappa through the first and second moment of sigma0(z) gives rise to physically implausible values for these observables. This discrepancy, which we previously observed for a much simpler CG model, suggests that the moment conditions derived from simple continuum assumptions miss the effect of physically important correlations in the lipid bilayer.
Collapse
Affiliation(s)
- Mingyang Hu
- Carnegie Mellon University, 5000 Forbes Avenue, Pittsburgh, PA 15213, USA
| | | | | | | |
Collapse
|
35
|
Wang Y, Sigurdsson JK, Brandt E, Atzberger PJ. Dynamic implicit-solvent coarse-grained models of lipid bilayer membranes: fluctuating hydrodynamics thermostat. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 88:023301. [PMID: 24032960 DOI: 10.1103/physreve.88.023301] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/05/2012] [Indexed: 06/02/2023]
Abstract
We introduce a thermostat based on fluctuating hydrodynamics for dynamic simulations of implicit-solvent coarse-grained models of lipid bilayer membranes. We show our fluctuating hydrodynamics approach captures interesting correlations in the dynamics of lipid bilayer membranes that are missing in simulations performed using standard Langevin dynamics. Our momentum conserving thermostat accounts for solvent-mediated momentum transfer by coupling coarse-grained degrees of freedom to stochastic continuum fields that account for both the solvent hydrodynamics and thermal fluctuations. We present both a general framework and specific methods to couple the particle and continuum degrees of freedom in a manner consistent with statistical mechanics and amenable to efficient computational simulation. For self-assembled vesicles, we study the diffusivity of lipids and their spatial correlations. We find the hydrodynamic coupling yields within the bilayer interesting correlations between diffusing lipids that manifest as a vortex-like structure similar to those observed in explicit-solvent simulations. We expect the introduced fluctuating hydrodynamics methods to provide a way to extend implicit-solvent models for use in a wide variety of dynamic studies.
Collapse
Affiliation(s)
- Yaohong Wang
- Department of Mathematics, University of California Santa Barbara, Santa Barbara, California, USA
| | | | | | | |
Collapse
|
36
|
Diggins P, Deserno M. Determining the bending modulus of a lipid membrane by simulating buckling. J Chem Phys 2013; 138:214110. [DOI: 10.1063/1.4808077] [Citation(s) in RCA: 61] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
37
|
Computational Studies of Biomembrane Systems: Theoretical Considerations, Simulation Models, and Applications. FROM SINGLE MOLECULES TO NANOSCOPICALLY STRUCTURED MATERIALS 2013. [DOI: 10.1007/12_2013_258] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
38
|
Determining the Gaussian curvature modulus of lipid membranes in simulations. Biophys J 2012; 102:1403-10. [PMID: 22455923 DOI: 10.1016/j.bpj.2012.02.013] [Citation(s) in RCA: 145] [Impact Index Per Article: 12.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2011] [Revised: 01/25/2012] [Accepted: 02/08/2012] [Indexed: 11/20/2022] Open
Abstract
The Gaussian curvature modulus κ¯ of lipid bilayers likely contributes more than 100 kcal/mol to every cellular fission or fusion event. This huge impact on membrane remodeling energetics might be a factor that codetermines the complex lipid composition of biomembranes through tuning of κ¯. Yet, its value has been measured only for a handful of simple lipids, and no simulation has so far determined it better than a factor of two, rendering a systematic investigation of such enticing speculations impossible. Here we propose a highly accurate method to determine κ¯ in computer simulations. It relies on the interplay between curvature stress and edge tension of partially curved axisymmetric membrane disks and requires determining their closing probability. For a simplified lipid model we obtain κ¯ and its relation to the normal bending modulus κ for membranes differing both in stiffness and spontaneous lipid curvature. The elastic ratio κ¯/κ can be determined with a few percent statistical accuracy. Its value agrees with the scarce experimental data, and its change with spontaneous lipid curvature is compatible with theoretical expectations, thereby granting additional information on monolayer properties. We also show that an alternative determination of these elastic parameters based on moments of the lateral stress profile gives markedly different and unphysical values.
Collapse
|
39
|
Brasiello A, Crescitelli S, Milano G. Development of a coarse-grained model for simulations of tridecanoin liquid-solid phase transitions. Phys Chem Chem Phys 2011; 13:16618-28. [PMID: 21858376 DOI: 10.1039/c1cp20604d] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Novel coarse-grained models for molecular dynamics of tridecanoin melts are here proposed as results of a coarse-graining step procedure. The procedure is implemented to develop three coarse-grained models of increasing number of particle types from two to four. Force fields are computed by minimization of the deviations of appropriate distribution functions of the coarse-grained models from those of a reference atomistic one. Density, diffusivity and shear viscosity are computed by numerical simulation and compared with experimental values. The ability of each model to describe liquid-solid transitions is also analyzed. In particular, the model with four types of coarse-grained beads shows a transition from a liquid to a crystal phase.
Collapse
Affiliation(s)
- Antonio Brasiello
- University of Salerno, Department of Industrial Engineering, Fisciano (SA), Italy.
| | | | | |
Collapse
|
40
|
Kawamoto S, Takasu M, Miyakawa T, Morikawa R, Oda T, Futaki S, Nagao H. Inverted micelle formation of cell-penetrating peptide studied by coarse-grained simulation: Importance of attractive force between cell-penetrating peptides and lipid head group. J Chem Phys 2011; 134:095103. [DOI: 10.1063/1.3555531] [Citation(s) in RCA: 57] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|