1
|
Sadeghi M, Rosenberger D. Dynamic framework for large-scale modeling of membranes and peripheral proteins. Methods Enzymol 2024; 701:457-514. [PMID: 39025579 DOI: 10.1016/bs.mie.2024.03.018] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/20/2024]
Abstract
In this chapter, we present a novel computational framework to study the dynamic behavior of extensive membrane systems, potentially in interaction with peripheral proteins, as an alternative to conventional simulation methods. The framework effectively describes the complex dynamics in protein-membrane systems in a mesoscopic particle-based setup. Furthermore, leveraging the hydrodynamic coupling between the membrane and its surrounding solvent, the coarse-grained model grounds its dynamics in macroscopic kinetic properties such as viscosity and diffusion coefficients, marrying the advantages of continuum- and particle-based approaches. We introduce the theoretical background and the parameter-space optimization method in a step-by-step fashion, present the hydrodynamic coupling method in detail, and demonstrate the application of the model at each stage through illuminating examples. We believe this modeling framework to hold great potential for simulating membrane and protein systems at biological spatiotemporal scales, and offer substantial flexibility for further development and parametrization.
Collapse
Affiliation(s)
- Mohsen Sadeghi
- Department of Mathematics and Computer Science, Freie Universität Berlin, Berlin, Germany.
| | | |
Collapse
|
2
|
Maeki M, Kimura N, Okada Y, Shimizu K, Shibata K, Miyazaki Y, Ishida A, Yonezawa K, Shimizu N, Shinoda W, Tokeshi M. Understanding the effects of ethanol on the liposome bilayer structure using microfluidic-based time-resolved small-angle X-ray scattering and molecular dynamics simulations. NANOSCALE ADVANCES 2024; 6:2166-2176. [PMID: 38633055 PMCID: PMC11019499 DOI: 10.1039/d3na01073b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/03/2023] [Accepted: 02/29/2024] [Indexed: 04/19/2024]
Abstract
Lipid nanoparticles (LNPs) are essential carrier particles in drug delivery systems, particularly in ribonucleic acid delivery. In preparing lipid-based nanoparticles, microfluidic-based ethanol injection may produce precisely size-controlled nanoparticles. Ethanol is critical in LNP formation and post-treatment processes and affects liposome size, structure, lamellarity, and drug-loading efficiency. However, the effects of time-dependent changes in the ethanol concentration on the structural dynamics of liposomes are not clearly understood. Herein, we investigated ethanol-induced lipid bilayer changes in liposomes on a time scale from microseconds to tens of seconds using a microfluidic-based small-angle X-ray scattering (SAXS) measurement system coupled with molecular dynamics (MD) simulations. The time-resolved SAXS measurement system revealed that single unilamellar liposomes were converted to multilamellar liposomes within 0.8 s of contact with ethanol, and the d-spacing was decreased from 6.1 (w/o ethanol) to 4.4 nm (80% ethanol) with increasing ethanol concentration. We conducted 1 μs MD simulations to understand the molecular-level structural changes in the liposomes. The MD simulations revealed that the changes in the lamellar structure caused by ethanol at the molecular level could explain the structural changes in the liposomes observed via time-resolved SAXS. Therefore, the post-treatment process to remove residual ethanol is critical in liposome formation.
Collapse
Affiliation(s)
- Masatoshi Maeki
- Division of Applied Chemistry, Faculty of Engineering, Hokkaido University Kita 13 Nishi 8, Kita-ku Sapporo 060-8628 Japan +81-11-706-6745 +81-11-706-6773 +81-11-706-6744
- JST PRESTO 4-1-8 Honcho, Kawaguchi Saitama 332-0012 Japan
- Institute of Materials Structure Science, High Energy Accelerator Research Organization (KEK) Tsukuba Ibaraki 305-0801 Japan
| | - Niko Kimura
- Graduate School of Chemical Sciences and Engineering, Hokkaido University Kita 13 Nishi 8, Kita-ku Sapporo 060-8628 Japan
| | - Yuto Okada
- Graduate School of Chemical Sciences and Engineering, Hokkaido University Kita 13 Nishi 8, Kita-ku Sapporo 060-8628 Japan
| | - Kazuki Shimizu
- Graduate School of Chemical Sciences and Engineering, Hokkaido University Kita 13 Nishi 8, Kita-ku Sapporo 060-8628 Japan
| | - Kana Shibata
- Department of Materials Chemistry, Nagoya University Chikusa-ku Nagoya 464-8603 Japan
| | - Yusuke Miyazaki
- Research Institute for Interdisciplinary Science, Okayama University Okayama 700-8530 Japan
| | - Akihiko Ishida
- Division of Applied Chemistry, Faculty of Engineering, Hokkaido University Kita 13 Nishi 8, Kita-ku Sapporo 060-8628 Japan +81-11-706-6745 +81-11-706-6773 +81-11-706-6744
| | - Kento Yonezawa
- Institute of Materials Structure Science, High Energy Accelerator Research Organization (KEK) Tsukuba Ibaraki 305-0801 Japan
| | - Nobutaka Shimizu
- Institute of Materials Structure Science, High Energy Accelerator Research Organization (KEK) Tsukuba Ibaraki 305-0801 Japan
| | - Wataru Shinoda
- Research Institute for Interdisciplinary Science, Okayama University Okayama 700-8530 Japan
| | - Manabu Tokeshi
- Division of Applied Chemistry, Faculty of Engineering, Hokkaido University Kita 13 Nishi 8, Kita-ku Sapporo 060-8628 Japan +81-11-706-6745 +81-11-706-6773 +81-11-706-6744
| |
Collapse
|
3
|
Shobhna, Dutta A, Kumari P, Kashyap HK. Stability of Cytoplasmic Membrane of Escherichia coli Bacteria in Aqueous and Ethanolic Environment. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2024; 40:2893-2906. [PMID: 38311936 DOI: 10.1021/acs.langmuir.3c02780] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/06/2024]
Abstract
The mechanism of action of any antibacterial agent or disinfectant depends largely on their interaction with the bacterial membrane. Herein, we use the SPICA (surface property fitting coarse graining) force-field and develop a coarse-grained (CG) model for the structure of the cytoplasmic membrane of Escherichia coli (E. coli) and its interaction with water and ethanol. We elucidate the impact of different concentrations of ethanol on the cytoplasmic membrane bilayers and vesicles of E. coli using the CG molecular dynamics (CG MD) simulations. Our modeling approach first focuses on the parametrization of the required force-field for POPG lipid and its interaction with water, ethanol, and POPE lipid. Subsequently, the structural stability of the E. coli bacterial membrane in the presence of high and low concentrations of ethanol is delineated. Both flat bilayers as well as vesicles of E. coli membrane were considered for the CG MD. Our results reveal that, at low ethanol concentrations (<30 mol %), the size of the E. coli vesicles increases with discernible deformations in their shapes. Because of ethanol-induced interdigitation, thinning of the E. coli vesicular membrane is also observed. However, at higher ethanol concentrations (>30 mol %), the integrity of the vesicles is lost because of deteriorating invasion of ethanol molecules into the vesicle bilayer and significant weakening of lipid-lipid interactions. At higher ethanol concentrations (40 and 70 mol %), both the multivesicle and single-vesicle bacterial membranes exhibit a similar rupturing pattern wherein the extraction of lipids from the membrane and formation of aggregates of the component lipids are observed. These aggregates consist of polar head groups of 3-5 POPE/POPG lipids with intertwined nonpolar tails.
Collapse
Affiliation(s)
- Shobhna
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| | - Ayishwarya Dutta
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| | - Pratibha Kumari
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| | - Hemant K Kashyap
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| |
Collapse
|
4
|
Zhang D, Wang Y, Safaripour M, Bellido-Aguilar DA, Van Donselaar KR, Webster DC, Croll AB, Xia W. Energy renormalization for temperature transferable coarse-graining of silicone polymer. Phys Chem Chem Phys 2024; 26:4541-4554. [PMID: 38241021 DOI: 10.1039/d3cp05969c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2024]
Abstract
The bottom-up prediction of thermodynamic and mechanical behaviors of polymeric materials based on molecular dynamics (MD) simulation is of critical importance in polymer physics. Although the atomistically informed coarse-grained (CG) model can access greater spatiotemporal scales and retain essential chemical specificity, the temperature-transferable CG model is still a big challenge and hinders widespread application of this technique. Herein, we use a silicone polymer, i.e., polydimethylsiloxane (PDMS), having an incredibly low chain rigidity as a model system, combined with an energy-renormalization (ER) approach, to systematically develop a temperature-transferable CG model. Specifically, by introducing temperature-dependent ER factors to renormalize the effective distance and cohesive energy parameters, the developed CG model faithfully preserved the dynamics, mechanical and conformational behaviors compared with the target all-atomistic (AA) model from glassy to melt regimes, which was further validated by experimental data. With the developed CG model featuring tremendously improved computational efficiency, we systematically explored the influences of cohesive interaction strength and temperature on the dynamical heterogeneity and mechanical response of polymers, where we observed consistent trends with other linear polymers with varying chain rigidity and monomeric structures. This study serves as an extension of our proposed ER approach of developing temperature transferable CG models with diverse segmental structures, highlighting the critical role of cohesive interaction strength on CG modeling of polymer dynamics and thermomechanical behaviors.
Collapse
Affiliation(s)
- Dawei Zhang
- Department of Civil, Construction, and Environmental Engineering, North Dakota State University, Fargo, ND 58108, USA
| | - Yang Wang
- Zernike Institute for Advanced Materials, University of Groningen, 9747 AG, Groningen, The Netherlands
| | - Maryam Safaripour
- Department of Coatings and Polymeric Materials, North Dakota State University, Fargo, ND 58108, USA
| | - Daniel A Bellido-Aguilar
- Department of Coatings and Polymeric Materials, North Dakota State University, Fargo, ND 58108, USA
| | | | - Dean C Webster
- Department of Coatings and Polymeric Materials, North Dakota State University, Fargo, ND 58108, USA
| | - Andrew B Croll
- Department of Physics, North Dakota State University, Fargo, ND 58108, USA
| | - Wenjie Xia
- Department of Aerospace Engineering, Iowa State University, Ames, IA 50011, USA.
| |
Collapse
|
5
|
Miyazaki Y, Shinoda W. pSPICA Force Field Extended for Proteins and Peptides. J Chem Inf Model 2024; 64:532-542. [PMID: 38156656 DOI: 10.1021/acs.jcim.3c01611] [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/03/2024]
Abstract
Many coarse-grained (CG) molecular dynamics (MD) studies have been performed to investigate biological processes involving proteins and lipids. CG force fields (FFs) in these MD studies often use implicit or nonpolar water models to reduce computational costs. CG-MD using water models cannot properly describe electrostatic screening effects owing to the hydration of ionic segments and thus cannot appropriately describe molecular events involving water channels and pores through lipid membranes. To overcome this issue, we developed a protein model in the pSPICA FF, in which a polar CG water model showing the proper dielectric response was adopted. The developed CG model greatly improved the transfer free energy profiles of charged side chain analogues across the lipid membrane. Application studies on melittin-induced membrane pores and mechanosensitive channels in lipid membranes demonstrated that CG-MDs using the pSPICA FF correctly reproduced the structure and stability of the pores and channels. Furthermore, the adsorption behavior of the highly charged nona-arginine peptides on lipid membranes changed with salt concentration, indicating the pSPICA FF is also useful for simulating protein adsorption on membrane surfaces.
Collapse
Affiliation(s)
- Yusuke Miyazaki
- Research Institute for Interdisciplinary Science, Okayama University, 3-1-1 Tsushima-naka, Kita-ku, Okayama 700-8530, Japan
| | - Wataru Shinoda
- Research Institute for Interdisciplinary Science, Okayama University, 3-1-1 Tsushima-naka, Kita-ku, Okayama 700-8530, Japan
| |
Collapse
|
6
|
Yamada T, Miyazaki Y, Harada S, Kumar A, Vanni S, Shinoda W. Improved Protein Model in SPICA Force Field. J Chem Theory Comput 2023; 19:8967-8977. [PMID: 37989551 DOI: 10.1021/acs.jctc.3c01016] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2023]
Abstract
The previous version of the SPICA coarse-grained (CG) force field (FF) protein model focused primarily on membrane proteins and successfully reproduced the dimerization free energies of several transmembrane helices and the stable structures of various membrane protein assemblies. However, that model had limited accuracy when applied to other proteins, such as intrinsically disordered proteins (IDPs) and peripheral proteins, because the dimensions of the IDPs in an aqueous solution were too compact, and protein binding on the lipid membrane surface was overstabilized. To improve the accuracy of the SPICA FF model for the simulation of such systems, in this study, we introduce protein secondary structure-dependent nonbonded interaction parameters to the backbone segments and reoptimize almost all nonbonded parameters for amino acids. The improved FF proposed here successfully reproduces the radii of gyration of various IDPs, the binding sensitivity of several peripheral membrane proteins, and the dimerization free energies of several transmembrane helices. The new model also shows improved agreement with experiments on the free energy of peptide association in water. In addition, an extensive library of nonbonded interactions between proteins and lipids, including various glycerophospholipids, sphingolipids, and cholesterol, allows the study of specific interactions between lipids and peripheral and transmembrane proteins. Hence, the new SPICA FF (version 2) proposed herein is applicable with high accuracy for simulating a wide range of protein systems.
Collapse
Affiliation(s)
- Teppei Yamada
- Graduate School of Natural Science and Technology, Okayama University, 3-1-1 Tsushima-naka, Kita-ku, Okayama 700-8530, Japan
| | - Yusuke Miyazaki
- Research Institute for Interdisciplinary Science, Okayama University, 3-1-1 Tsushima-naka, Kita-ku, Okayama 700-8530, Japan
| | - Shogo Harada
- Department of Materials Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan
| | - Ashutosh Kumar
- Department of Biology and National Center of Competence in Research Bio-inspired Materials, University of Fribourg, Chemin du Musée 10, 1700 Fribourg, Switzerland
| | - Stefano Vanni
- Department of Biology and National Center of Competence in Research Bio-inspired Materials, University of Fribourg, Chemin du Musée 10, 1700 Fribourg, Switzerland
| | - Wataru Shinoda
- Research Institute for Interdisciplinary Science, Okayama University, 3-1-1 Tsushima-naka, Kita-ku, Okayama 700-8530, Japan
| |
Collapse
|
7
|
Borges-Araújo L, Patmanidis I, Singh AP, Santos LHS, Sieradzan AK, Vanni S, Czaplewski C, Pantano S, Shinoda W, Monticelli L, Liwo A, Marrink SJ, Souza PCT. Pragmatic Coarse-Graining of Proteins: Models and Applications. J Chem Theory Comput 2023; 19:7112-7135. [PMID: 37788237 DOI: 10.1021/acs.jctc.3c00733] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
The molecular details involved in the folding, dynamics, organization, and interaction of proteins with other molecules are often difficult to assess by experimental techniques. Consequently, computational models play an ever-increasing role in the field. However, biological processes involving large-scale protein assemblies or long time scale dynamics are still computationally expensive to study in atomistic detail. For these applications, employing coarse-grained (CG) modeling approaches has become a key strategy. In this Review, we provide an overview of what we call pragmatic CG protein models, which are strategies combining, at least in part, a physics-based implementation and a top-down experimental approach to their parametrization. In particular, we focus on CG models in which most protein residues are represented by at least two beads, allowing these models to retain some degree of chemical specificity. A description of the main modern pragmatic protein CG models is provided, including a review of the most recent applications and an outlook on future perspectives in the field.
Collapse
Affiliation(s)
- Luís Borges-Araújo
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS, University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| | - Ilias Patmanidis
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Akhil P Singh
- Department of Biology, University of Fribourg, Chemin du Musée 10, Fribourg CH-1700, Switzerland
| | - Lucianna H S Santos
- Biomolecular Simulations Group, Institut Pasteur de Montevideo, Montevideo 11400, Uruguay
| | - Adam K Sieradzan
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Stefano Vanni
- Department of Biology, University of Fribourg, Chemin du Musée 10, Fribourg CH-1700, Switzerland
- Institut de Pharmacologie Moléculaire et Cellulaire, Université Côte d'Azur, Inserm, CNRS, 06560 Valbonne, France
| | - Cezary Czaplewski
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Sergio Pantano
- Biomolecular Simulations Group, Institut Pasteur de Montevideo, Montevideo 11400, Uruguay
| | - Wataru Shinoda
- Research Institute for Interdisciplinary Science, Okayama University, 3-1-1 Tsushima-naka, Kita, Okayama 700-8530, Japan
| | - Luca Monticelli
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS, University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| | - Adam Liwo
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Paulo C T Souza
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS, University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| |
Collapse
|
8
|
Koshiyama K, Nakata K. Effects of lipid saturation on bicelle to vesicle transition of a binary phospholipid mixture: a molecular dynamics simulation study. SOFT MATTER 2023; 19:7655-7662. [PMID: 37782209 DOI: 10.1039/d3sm00904a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/03/2023]
Abstract
Controlling the transition from lipid bicelles to vesicles is essential for producing engineered vesicles. We perform coarse-grained molecular dynamics (CGMD) simulations of unsaturated/saturated lipid mixtures to clarify the effects of lipid unsaturation on vesiculation at the molecular scale. The results demonstrate that vesiculation depends on the concentration of unsaturated lipids and the degree of unsaturation. The probability of vesiculation increases linearly with the apparent unsaturated lipid concentration at a low degree of unsaturation. Higher degrees of unsaturation lead to phase segregation within the binary bicelles, reducing the probability of vesiculation. A comparison between CGMD simulations and the conventional theory of vesiculation shows that the theoretical predictions of binary lipid systems must explicitly include phase segregation effects. Furthermore, simulations with biased lipid distributions reveal that vesiculation is facilitated by the preconcentration of unsaturated lipids in the core region of the bicelle but is then temporally limited as the unsaturated lipids move to the bicelle edges. These findings advance theoretical and experimental studies on binary lipid systems and promote the development of tailor-made vesicles.
Collapse
Affiliation(s)
- Kenichiro Koshiyama
- Graduate School of Technology, Industrial and Social Sciences, Tokushima University, Tokushima, 770-8506, Japan.
| | - Kazuki Nakata
- Graduate School of Sciences and Technology for Innovation, Tokushima University, Tokushima 770-8506, Japan.
| |
Collapse
|
9
|
Ma M, Song J, Dong Y, Fang W, Gao L. Structural and thermodynamic properties of bulk triglycerides and triglyceride/water mixtures reproduced using a polarizable coarse-grained model. Phys Chem Chem Phys 2023; 25:22232-22243. [PMID: 37577752 DOI: 10.1039/d3cp01839c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/15/2023]
Abstract
Triglycerides (TGs) play important roles in renewable energies, food production, medicine, and metabolism in organisms. Here, we developed a novel coarse-grained (CG) force field (FF) for triglycerides to reproduce both the structural and thermodynamic properties of bulk TGs, TG/air interfaces, and TG/water mixtures using molecular dynamics (MD) simulations. We rigorously optimized the bonded and nonbonded force parameters between the CG beads of TGs and nonbonded force parameters between TG beads and polarizable CG water beads by employing an efficient meta-multilinear interpolation parameterization algorithm recently developed by us. This CG FF performs very well in reproducing the percolating network of the TG bulk phase self-assembled in water and a variety of molecular conformations predicted by all-atom MD simulations. More importantly, it also correctly reproduces multiple experimentally measurable macroscopic thermodynamic properties, including the density and surface tensions of both the TG/air and TG/water interfaces. This paves the way for studying more complicated systems involving TGs on a large scale.
Collapse
Affiliation(s)
- Ming Ma
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, 19 Xin-Jie-Kou-Wai Street, Beijing 100875, China.
| | - Junjie Song
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, 19 Xin-Jie-Kou-Wai Street, Beijing 100875, China.
| | - Yi Dong
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, 19 Xin-Jie-Kou-Wai Street, Beijing 100875, China.
| | - Weihai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, 19 Xin-Jie-Kou-Wai Street, Beijing 100875, China.
| | - Lianghui Gao
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, 19 Xin-Jie-Kou-Wai Street, Beijing 100875, China.
| |
Collapse
|
10
|
Prabhu J, Singh AP, Vanni S. An in silico osmotic pressure approach allows characterization of pressure-area isotherms of lipid monolayers at low molecular areas. SOFT MATTER 2023; 19:3377-3385. [PMID: 37102755 PMCID: PMC10170484 DOI: 10.1039/d2sm01419j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
Surface pressure-area isotherms of lipid monolayers at the air-water interface provide essential information about the structure and mechanical behaviour of lipid membranes. These curves can be readily obtained through Langmuir trough measurements and, as such, have been collected for decades in the field of membrane biochemistry. However, it is still challenging to directly observe and understand nanoscopic features of monolayers through such experiments, and molecular dynamics (MD) simulations are generally used to provide a molecular view of such interfaces. In MD simulations, the surface pressure-area (Π-A) isotherms are generally computed using the Kirkwood-Irving formula, that relies on the evaluation of the pressure tensor. This approach, however, has intrinsic limitations when the molecular area in the monolayer is low (typically < 60 Å2 per lipid). Recently, an alternative method to compute Π-A isotherms of surfactants, based on the calculation of the three-dimensional osmotic pressure via the implementation of semipermeable barriers was proposed. In this work, we investigate the feasibility of this approach for long-chain surfactants such as phospholipids. We identify some discrepancies between the computed values and experimental results, and we propose a semi-empirical correction based on the molecular structure of the surfactants at the monolayer interface. To validate the potential of this new approach, we simulate several phosphatidylcholine and phosphatidylethanolamine lipids at various temperatures using all-atom and coarse-grained force fields, and we compute the corresponding Π-A isotherms. Our results show that the Π-A isotherms obtained using the new method are in very good agreement with experiments and far superior to the canonical pressure tensor-based method at low molecular areas. This corrected osmotic pressure method allows for accurate characterization of the molecular packing in monolayers in various physical phases.
Collapse
Affiliation(s)
- Janak Prabhu
- Department of Biology, University of Fribourg, Chemin du Musée 10, 1700 Fribourg, Switzerland.
| | - Akhil Pratap Singh
- Department of Biology, University of Fribourg, Chemin du Musée 10, 1700 Fribourg, Switzerland.
| | - Stefano Vanni
- Department of Biology, University of Fribourg, Chemin du Musée 10, 1700 Fribourg, Switzerland.
| |
Collapse
|
11
|
Xu X, Liu A, Liu S, Ma Y, Zhang X, Zhang M, Zhao J, Sun S, Sun X. Application of molecular dynamics simulation in self-assembled cancer nanomedicine. Biomater Res 2023; 27:39. [PMID: 37143168 PMCID: PMC10161522 DOI: 10.1186/s40824-023-00386-7] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2022] [Accepted: 04/21/2023] [Indexed: 05/06/2023] Open
Abstract
Self-assembled nanomedicine holds great potential in cancer theragnostic. The structures and dynamics of nanomedicine can be affected by a variety of non-covalent interactions, so it is essential to ensure the self-assembly process at atomic level. Molecular dynamics (MD) simulation is a key technology to link microcosm and macroscale. Along with the rapid development of computational power and simulation methods, scientists could simulate the specific process of intermolecular interactions. Thus, some experimental observations could be explained at microscopic level and the nanomedicine synthesis process would have traces to follow. This review not only outlines the concept, basic principle, and the parameter setting of MD simulation, but also highlights the recent progress in MD simulation for self-assembled cancer nanomedicine. In addition, the physicochemical parameters of self-assembly structure and interaction between various assembled molecules under MD simulation are also discussed. Therefore, this review will help advanced and novice researchers to quickly zoom in on fundamental information and gather some thought-provoking ideas to advance this subfield of self-assembled cancer nanomedicine.
Collapse
Affiliation(s)
- Xueli Xu
- School of Science, Shandong Jianzhu University, Jinan, 250101, China
| | - Ao Liu
- School of Science, Shandong Jianzhu University, Jinan, 250101, China
| | - Shuangqing Liu
- School of Chemistry and Pharmaceutical Engineering, Medical Science and Technology Innovation Center, Shandong First Medical University & Shandong Academy of Medical Sciences, Jinan, 250000, China
| | - Yanling Ma
- School of Chemistry and Pharmaceutical Engineering, Medical Science and Technology Innovation Center, Shandong First Medical University & Shandong Academy of Medical Sciences, Jinan, 250000, China
| | - Xinyu Zhang
- School of Chemistry and Pharmaceutical Engineering, Medical Science and Technology Innovation Center, Shandong First Medical University & Shandong Academy of Medical Sciences, Jinan, 250000, China
| | - Meng Zhang
- School of Chemistry and Pharmaceutical Engineering, Medical Science and Technology Innovation Center, Shandong First Medical University & Shandong Academy of Medical Sciences, Jinan, 250000, China
| | - Jinhua Zhao
- School of Science, Shandong Jianzhu University, Jinan, 250101, China
| | - Shuo Sun
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, 02115, USA
| | - Xiao Sun
- School of Chemistry and Pharmaceutical Engineering, Medical Science and Technology Innovation Center, Shandong First Medical University & Shandong Academy of Medical Sciences, Jinan, 250000, China.
| |
Collapse
|
12
|
Ledum M, Sen S, Li X, Carrer M, Feng Y, Cascella M, Bore SL. HylleraasMD: A Domain Decomposition-Based Hybrid Particle-Field Software for Multiscale Simulations of Soft Matter. J Chem Theory Comput 2023; 19:2939-2952. [PMID: 37130290 DOI: 10.1021/acs.jctc.3c00134] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
We present HylleraasMD (HyMD), a comprehensive implementation of the recently proposed Hamiltonian formulation of hybrid particle-field molecular dynamics. The methodology is based on a tunable, grid-independent length-scale of coarse graining, obtained by filtering particle densities in reciprocal space. This enables systematic convergence of energies and forces by grid refinement, also eliminating nonphysical force aliasing. Separating the time integration of fast modes associated with internal molecular motion from slow modes associated with their density fields, we enable the first time-reversible, energy-conserving hybrid particle-field simulations. HyMD comprises the optional use of explicit electrostatics, which, in this formalism, corresponds to the long-range potential in particle-mesh Ewald. We demonstrate the ability of HyMD to perform simulations in the microcanonical and canonical ensembles with a series of test cases, comprising lipid bilayers and vesicles, surfactant micelles, and polypeptide chains, comparing our results to established literature. An on-the-fly increase of the characteristic coarse-grain length significantly speeds up dynamics, accelerating self-diffusion and leading to expedited aggregation. Exploiting this acceleration, we find that the time scales involved in the self-assembly of polymeric structures can lie in the tens to hundreds of picoseconds instead of the multimicrosecond regime observed with comparable coarse-grained models.
Collapse
Affiliation(s)
- Morten Ledum
- Department of Chemistry and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, PO Box 1033 Blindern, 0315 Oslo, Norway
| | - Samiran Sen
- Department of Chemistry and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, PO Box 1033 Blindern, 0315 Oslo, Norway
| | - Xinmeng Li
- Department of Chemistry and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, PO Box 1033 Blindern, 0315 Oslo, Norway
| | - Manuel Carrer
- Department of Chemistry and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, PO Box 1033 Blindern, 0315 Oslo, Norway
| | - Yu Feng
- Berkeley Center for Cosmological Physics and Department of Physics, University of California, Berkeley, California 94720, United States
| | - Michele Cascella
- Department of Chemistry and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, PO Box 1033 Blindern, 0315 Oslo, Norway
| | - Sigbjørn Løland Bore
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
13
|
Petretto E, Campomanes P, Vanni S. Development of a coarse-grained model for surface-functionalized gold nanoparticles: towards an accurate description of their aggregation behavior. SOFT MATTER 2023; 19:3290-3300. [PMID: 37092690 PMCID: PMC10170483 DOI: 10.1039/d3sm00094j] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
Understanding the dispersion stability and aggregation propensity of self-assembled monolayer gold NPs at a molecular level is crucial to guide their rational design and to inform about the optimal surface functionalization for specific applications. To reach this goal, in silico modeling via coarse-grained (CG) molecular dynamics (MD) simulations is a fundamental tool to complement the information acquired from experimental studies since CG modeling allows to get a deep knowledge of the molecular interactions that take place at the nanoscale in this kind of systems. Unfortunately, current CG models of monolayer-protected AuNPs present several drawbacks that limit their accuracy in certain scenarios. We here develop a CG model that is fully compatible and extends the SPICA/SDK (Shinoda-DeVane-Klein) force field. Our model allows reproducing the behavior of AuNPs functionalized with hydrophobic as well as charged and more hydrophilic ligands. This model improves upon results obtained with previously derived CG force fields and successfully describes NPs aggregation and self-assembly in aqueous solution.
Collapse
Affiliation(s)
- Emanuele Petretto
- Department of Biology, University of Fribourg, Chemin du Musée 10, 1700 Fribourg, Switzerland.
| | - Pablo Campomanes
- Department of Biology, University of Fribourg, Chemin du Musée 10, 1700 Fribourg, Switzerland.
| | - Stefano Vanni
- Department of Biology, University of Fribourg, Chemin du Musée 10, 1700 Fribourg, Switzerland.
| |
Collapse
|
14
|
Pan F, Sun L, Li S. Dynamic Processes and Mechanical Properties of Lipid-Nanoparticle Mixtures. Polymers (Basel) 2023; 15:polym15081828. [PMID: 37111975 PMCID: PMC10144953 DOI: 10.3390/polym15081828] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2022] [Revised: 02/23/2023] [Accepted: 04/03/2023] [Indexed: 04/29/2023] Open
Abstract
In this study, we investigate the dynamic processes and mechanical properties of lipid nanoparticle mixtures in a melt via dissipation particle dynamic simulation. By investigating the distribution of nanoparticles in lamellar and hexagonal lipid matrices in equilibrium state and dynamic processes, we observe that the morphology of such composites depends not only on the geometric features of the lipid matrix but also on the concentration of nanoparticles. The dynamic processes are also demonstrated by calculating the average radius of gyration, which indicates the isotropic conformation of lipid molecules in the x-y plane and that the lipid chains are stretched in the z direction with the addition of nanoparticles. Meanwhile, we predict the mechanical properties of lipid-nanoparticle mixtures in lamellar structures by analyzing the interfacial tensions. Results show that the interfacial tension decreased with the increase in nanoparticle concentration. These results provide molecular-level information for the rational and a priori design of new lipid nanocomposites with ad hoc tailored properties.
Collapse
Affiliation(s)
- Fan Pan
- School of Data Science and Artificial Intelligence, Wenzhou University of Technology, Wenzhou 325035, China
| | - Lingling Sun
- Department of Physics, Wenzhou University, Wenzhou 325035, China
| | - Shiben Li
- Department of Physics, Wenzhou University, Wenzhou 325035, China
| |
Collapse
|
15
|
Wan M, Song J, Yang Y, Gao L, Fang W. A top-down and bottom-up combined strategy for parameterization of coarse-grained force fields for phospholipids. Phys Chem Chem Phys 2023; 25:6757-6767. [PMID: 36789502 DOI: 10.1039/d2cp05384e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2023]
Abstract
Coarse-graining (CG) molecular dynamics (MD) simulations are widely used in interpreting experimental observations and predicting assembly morphology as well as collective behaviour but also face the problem of poor accuracy. A main issue is that cross-termed interactions between different CG beads are inadequately parameterized. This work proposes a novel top-down and bottom-up combined strategy to parameterize both self- and cross-termed interactions of zwitterionic phospholipids in water solution based on a piecewise Morse potential describing nonbonded van der Waals interactions. The self-interacting force parameters were optimized by matching experimental density, heat vapourization, and surface tension in a top-down manner, while the cross-termed interactions were optimized by fitting pseudo properties obtained from atomistic simulations in a bottom-up way, including mixing density, intermolecular energy, and radial mixing coefficient. The transferability of the CG force field (FF) was confirmed by reproducing a variety of structural and thermodynamic properties of lipid membranes in both liquid and gel phases. This FF can well depict vesicle self-assembly and vesicle fusion processes. Matching pseudo properties opens a new way to develop CG FF with increased accuracy and transferability.
Collapse
Affiliation(s)
- Mingwei Wan
- Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China. .,Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, 19 Xin-Jie-Kou-Wai Street, Beijing 100875, China.
| | - Junjie Song
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, 19 Xin-Jie-Kou-Wai Street, Beijing 100875, China.
| | - Ying Yang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, 19 Xin-Jie-Kou-Wai Street, Beijing 100875, China.
| | - Lianghui Gao
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, 19 Xin-Jie-Kou-Wai Street, Beijing 100875, China.
| | - Weihai Fang
- Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China. .,Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, 19 Xin-Jie-Kou-Wai Street, Beijing 100875, China.
| |
Collapse
|
16
|
Zhang D, Xiao Q, Rahimzadeh M, Liu M, Rodriguez-Emmenegger C, Miyazaki Y, Shinoda W, Percec V. Self-Assembly of Glycerol-Amphiphilic Janus Dendrimers Amplifies and Indicates Principles for the Selection of Stereochemistry by Biological Membranes. J Am Chem Soc 2023; 145:4311-4323. [PMID: 36749951 DOI: 10.1021/jacs.3c00389] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023]
Abstract
The principles for the selection of the stereochemistry of phospholipids of biological membranes remain unclear and continue to be debated. Therefore, any new experiments on this topic may help progress in this field. To address this question, three libraries of constitutional isomeric glycerol-amphiphilic Janus dendrimers (JDs) with nonsymmetric homochiral, racemic, and symmetric achiral branching points were synthesized by an orthogonal-modular-convergent methodology. These JDs amplify self-assembly, and therefore, monodisperse vesicles known as dendrimersomes (DSs) with predictable dimensions programmed by JD concentration were assembled by rapid injection of their ethanol solution into water. DSs of homochiral JD enantiomers, racemic, including mixtures of different enantiomers, and achiral exhibited similar DS size-concentration dependence. However, the number of bilayers of DSs assembled from homochiral, achiral, and racemic JDs determined by cryo-TEM were different. Statistical analysis of the number of bilayers and coarse-grained molecular dynamics simulations demonstrated that homochiral JDs formed predominantly unilamellar DSs. Symmetric achiral JDs assembled only unilamellar DSs while racemic JDs favored multilamellar DSs. Since cell membranes are unilamellar, these results indicate a new rationale for nonsymmetric homochiral vs racemic selection. Simultaneously, these experiments imply that the symmetric achiral lipids forming more stable membrane, probably had been the preferable assemblies of prebiotic cell membranes.
Collapse
Affiliation(s)
- Dapeng Zhang
- Roy & Diana Vagelos Laboratories, Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104-6323, United States
| | - Qi Xiao
- Roy & Diana Vagelos Laboratories, Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104-6323, United States
- Institute of Computational Molecular Science, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Mehrnoush Rahimzadeh
- DWI─Leibniz Institute for Interactive Materials, Aachen 52074, Germany
- Institute of Technical and Macromolecular Chemistry, RWTH Aachen University, Aachen 52074, Germany
| | - Matthew Liu
- Roy & Diana Vagelos Laboratories, Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104-6323, United States
| | - Cesar Rodriguez-Emmenegger
- DWI─Leibniz Institute for Interactive Materials, Aachen 52074, Germany
- Institute of Technical and Macromolecular Chemistry, RWTH Aachen University, Aachen 52074, Germany
| | - Yusuke Miyazaki
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| | - Wataru Shinoda
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
| | - Virgil Percec
- Roy & Diana Vagelos Laboratories, Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104-6323, United States
| |
Collapse
|
17
|
Rogers S, Gui L, Kovalenko A, Zoni V, Carpentier M, Ramji K, Ben Mbarek K, Bacle A, Fuchs P, Campomanes P, Reetz E, Speer NO, Reynolds E, Thiam AR, Vanni S, Nicastro D, Henne WM. Triglyceride lipolysis triggers liquid crystalline phases in lipid droplets and alters the LD proteome. J Cell Biol 2022; 221:213472. [PMID: 36112368 PMCID: PMC9485706 DOI: 10.1083/jcb.202205053] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2022] [Revised: 08/10/2022] [Accepted: 08/22/2022] [Indexed: 01/08/2023] Open
Abstract
Lipid droplets (LDs) are reservoirs for triglycerides (TGs) and sterol-esters (SEs), but how these lipids are organized within LDs and influence their proteome remain unclear. Using in situ cryo-electron tomography, we show that glucose restriction triggers lipid phase transitions within LDs generating liquid crystalline lattices inside them. Mechanistically this requires TG lipolysis, which decreases the LD's TG:SE ratio, promoting SE transition to a liquid crystalline phase. Molecular dynamics simulations reveal TG depletion promotes spontaneous TG and SE demixing in LDs, additionally altering the lipid packing of the PL monolayer surface. Fluorescence imaging and proteomics further reveal that liquid crystalline phases are associated with selective remodeling of the LD proteome. Some canonical LD proteins, including Erg6, relocalize to the ER network, whereas others remain LD-associated. Model peptide LiveDrop also redistributes from LDs to the ER, suggesting liquid crystalline phases influence ER-LD interorganelle transport. Our data suggests glucose restriction drives TG mobilization, which alters the phase properties of LD lipids and selectively remodels the LD proteome.
Collapse
Affiliation(s)
- Sean Rogers
- Department of Cell Biology, University of Texas Southwestern Medical Center, Dallas, TX
| | - Long Gui
- Department of Cell Biology, University of Texas Southwestern Medical Center, Dallas, TX
| | - Anastasiia Kovalenko
- Department of Cell Biology, University of Texas Southwestern Medical Center, Dallas, TX
| | - Valeria Zoni
- Department of Biology, University of Fribourg, Fribourg, Switzerland
| | - Maxime Carpentier
- Laboratoire de Physique de l'École Normale Supérieure, École normale supérieure, Université Paris Sciences et Lettres, Centre national de la recherche scientifique, Sorbonne Université, Université de Paris, Paris, France
| | - Kamran Ramji
- Laboratoire de Physique de l'École Normale Supérieure, École normale supérieure, Université Paris Sciences et Lettres, Centre national de la recherche scientifique, Sorbonne Université, Université de Paris, Paris, France
| | - Kalthoum Ben Mbarek
- Laboratoire de Physique de l'École Normale Supérieure, École normale supérieure, Université Paris Sciences et Lettres, Centre national de la recherche scientifique, Sorbonne Université, Université de Paris, Paris, France
| | - Amelie Bacle
- Institute Jacques Monod, Centre national de la recherche scientifique, University of Paris, Paris, France
| | - Patrick Fuchs
- Laboratoire des Biomolécules, Paris, France.,Université de Paris, UFR Sciences du Vivant, Paris, France
| | - Pablo Campomanes
- Department of Biology, University of Fribourg, Fribourg, Switzerland
| | - Evan Reetz
- Department of Cell Biology, University of Texas Southwestern Medical Center, Dallas, TX
| | - Natalie Ortiz Speer
- Department of Cell Biology, University of Texas Southwestern Medical Center, Dallas, TX
| | - Emma Reynolds
- Department of Cell Biology, University of Texas Southwestern Medical Center, Dallas, TX
| | - Abdou Rachid Thiam
- Laboratoire de Physique de l'École Normale Supérieure, École normale supérieure, Université Paris Sciences et Lettres, Centre national de la recherche scientifique, Sorbonne Université, Université de Paris, Paris, France
| | - Stefano Vanni
- Department of Biology, University of Fribourg, Fribourg, Switzerland
| | - Daniela Nicastro
- Department of Cell Biology, University of Texas Southwestern Medical Center, Dallas, TX
| | - W Mike Henne
- Department of Cell Biology, University of Texas Southwestern Medical Center, Dallas, TX
| |
Collapse
|
18
|
Miyazaki Y, Shinoda W. Cooperative antimicrobial action of melittin on lipid membranes: A coarse-grained molecular dynamics study. BIOCHIMICA ET BIOPHYSICA ACTA. BIOMEMBRANES 2022; 1864:183955. [PMID: 35526599 DOI: 10.1016/j.bbamem.2022.183955] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/13/2021] [Revised: 04/09/2022] [Accepted: 04/29/2022] [Indexed: 12/29/2022]
Abstract
We conducted a series of coarse-grained molecular dynamics (CG-MD) simulations to investigate the complicated actions of melittin, which is an antimicrobial peptide (AMP) derived from honey bee venom, on a lipid membrane. To accurately simulate the AMP action, we developed and used a protein CG model as an extension of the pSPICA force field (FF), which was designed to reproduce several thermodynamic quantities and structural properties. At a low peptide-to-lipid (P/L) ratio (1/102), no defect was detected. At P/L = 1/51, toroidal pore formation was observed due to collective insertion of multiple melittin peptides from the N-termini. The pore formation was initiated by a local increase in membrane curvature in the vicinity of the peptide aggregate. At a higher P/L ratio (1/26), two more modes were detected, seemingly not controlled by the P/L ratio but by a local arrangement of melittin peptides: 1. Pore formation accompanied by lipid extraction by melittin peptides:a detergent-like mechanism. 2. A rapidly formed large pore in a significantly curved membrane: bursting. Thus, we observed three pore formation modes (toroidal pore formation, lipid extraction, and bursting) depending on the peptide concentration and local arrangement. These observations were consistent with experimental observations and hypothesized melittin modes. Through this study, we found that the local arrangements and population of melittin peptides and the area expansion rate by membrane deformation were key to the initiation of and competition among the multiple pore formation mechanisms.
Collapse
Affiliation(s)
- Yusuke Miyazaki
- Department of Materials Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan; Research Institute for Interdisciplinary Science, Okayama University, 3-1-1 Tsushima-naka, Kita-ku, Okayama 700-8530, Japan
| | - Wataru Shinoda
- Department of Materials Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan; Research Institute for Interdisciplinary Science, Okayama University, 3-1-1 Tsushima-naka, Kita-ku, Okayama 700-8530, Japan; Department of Chemistry, Faculty of Science, Okayama University, 3-1-1 Tsushima-naka, Kita-ku, Okayama 700-8530, Japan.
| |
Collapse
|
19
|
Shamaprasad P, Frame CO, Moore TC, Yang A, Iacovella CR, Bouwstra JA, Bunge AL, McCabe C. Using molecular simulation to understand the skin barrier. Prog Lipid Res 2022; 88:101184. [PMID: 35988796 PMCID: PMC10116345 DOI: 10.1016/j.plipres.2022.101184] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2022] [Revised: 08/11/2022] [Accepted: 08/12/2022] [Indexed: 11/15/2022]
Abstract
Skin's effectiveness as a barrier to permeation of water and other chemicals rests almost entirely in the outermost layer of the epidermis, the stratum corneum (SC), which consists of layers of corneocytes surrounded by highly organized lipid lamellae. As the only continuous path through the SC, transdermal permeation necessarily involves diffusion through these lipid layers. The role of the SC as a protective barrier is supported by its exceptional lipid composition consisting of ceramides (CERs), cholesterol (CHOL), and free fatty acids (FFAs) and the complete absence of phospholipids, which are present in most biological membranes. Molecular simulation, which provides molecular level detail of lipid configurations that can be connected with barrier function, has become a popular tool for studying SC lipid systems. We review this ever-increasing body of literature with the goals of (1) enabling the experimental skin community to understand, interpret and use the information generated from the simulations, (2) providing simulation experts with a solid background in the chemistry of SC lipids including the composition, structure and organization, and barrier function, and (3) presenting a state of the art picture of the field of SC lipid simulations, highlighting the difficulties and best practices for studying these systems, to encourage the generation of robust reproducible studies in the future. This review describes molecular simulation methodology and then critically examines results derived from simulations using atomistic and then coarse-grained models.
Collapse
Affiliation(s)
- Parashara Shamaprasad
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, TN 37235-1604, United States of America; Multiscale Modeling and Simulation (MuMS) Center, Vanderbilt University, Nashville, TN 37235-1604, United States of America
| | - Chloe O Frame
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, TN 37235-1604, United States of America; Multiscale Modeling and Simulation (MuMS) Center, Vanderbilt University, Nashville, TN 37235-1604, United States of America
| | - Timothy C Moore
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, TN 37235-1604, United States of America; Multiscale Modeling and Simulation (MuMS) Center, Vanderbilt University, Nashville, TN 37235-1604, United States of America
| | - Alexander Yang
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, TN 37235-1604, United States of America; Multiscale Modeling and Simulation (MuMS) Center, Vanderbilt University, Nashville, TN 37235-1604, United States of America
| | - Christopher R Iacovella
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, TN 37235-1604, United States of America; Multiscale Modeling and Simulation (MuMS) Center, Vanderbilt University, Nashville, TN 37235-1604, United States of America
| | - Joke A Bouwstra
- Division of BioTherapeutics, LACDR, Leiden University, 2333 CC Leiden, the Netherlands
| | - Annette L Bunge
- Department of Chemical and Biological Engineering, Colorado School of Mines, Golden, CO 80401, United States of America
| | - Clare McCabe
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, TN 37235-1604, United States of America; Multiscale Modeling and Simulation (MuMS) Center, Vanderbilt University, Nashville, TN 37235-1604, United States of America; School of Engineering and Physical Science, Heriot-Watt University, Edinburgh, United Kingdom.
| |
Collapse
|
20
|
Javan Nikkhah S, Vandichel M. Modeling Polyzwitterion-Based Drug Delivery Platforms: A Perspective of the Current State-of-the-Art and Beyond. ACS ENGINEERING AU 2022; 2:274-294. [PMID: 35996394 PMCID: PMC9389590 DOI: 10.1021/acsengineeringau.2c00008] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
![]()
Drug delivery platforms
are anticipated to have biocompatible and
bioinert surfaces. PEGylation of drug carriers is the most approved
method since it improves water solubility and colloid stability and
decreases the drug vehicles’ interactions with blood components.
Although this approach extends their biocompatibility, biorecognition
mechanisms prevent them from biodistribution and thus efficient drug
transfer. Recent studies have shown (poly)zwitterions to be alternatives
for PEG with superior biocompatibility. (Poly)zwitterions are super
hydrophilic, mainly stimuli-responsive, easy to functionalize and
they display an extremely low protein adsorption and long biodistribution
time. These unique characteristics make them already promising candidates
as drug delivery carriers. Furthermore, since they have highly dense
charged groups with opposite signs, (poly)zwitterions are intensely
hydrated under physiological conditions. This exceptional hydration
potential makes them ideal for the design of therapeutic vehicles
with antifouling capability, i.e., preventing undesired
sorption of biologics from the human body in the drug delivery vehicle.
Therefore, (poly)zwitterionic materials have been broadly applied
in stimuli-responsive “intelligent” drug delivery systems
as well as tumor-targeting carriers because of their excellent biocompatibility,
low cytotoxicity, insignificant immunogenicity, high stability, and
long circulation time. To tailor (poly)zwitterionic drug vehicles,
an interpretation of the structural and stimuli-responsive behavior
of this type of polymer is essential. To this end, a direct study
of molecular-level interactions, orientations, configurations, and
physicochemical properties of (poly)zwitterions is required, which
can be achieved via molecular modeling, which has become an influential
tool for discovering new materials and understanding diverse material
phenomena. As the essential bridge between science and engineering,
molecular simulations enable the fundamental understanding of the
encapsulation and release behavior of intelligent drug-loaded (poly)zwitterion
nanoparticles and can help us to systematically design their next
generations. When combined with experiments, modeling can make quantitative
predictions. This perspective article aims to illustrate key recent
developments in (poly)zwitterion-based drug delivery systems. We summarize
how to use predictive multiscale molecular modeling techniques to
successfully boost the development of intelligent multifunctional
(poly)zwitterions-based systems.
Collapse
Affiliation(s)
- Sousa Javan Nikkhah
- Department of Chemical Sciences, Bernal Institute, University of Limerick, Limerick V94 T9PX, Republic of Ireland
| | - Matthias Vandichel
- Department of Chemical Sciences, Bernal Institute, University of Limerick, Limerick V94 T9PX, Republic of Ireland
| |
Collapse
|
21
|
Wang ZS, Zhou ZH, Han L, Chen X, He HJ, Zhang Q, Xu ZC, Gong QT, Zhang L, Ma GY, Zhang L. The mechanism for lowering interfacial tension by extended surfactant containing ethylene oxide and propylene oxide groups. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.119364] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
22
|
Di Bartolo AL, Tomes CN, Mayorga LS, Masone D. Enhanced Expansion and Reduced Kiss-and-Run Events in Fusion Pores Steered by Synaptotagmin-1 C2B Domains. J Chem Theory Comput 2022; 18:4544-4554. [PMID: 35759758 DOI: 10.1021/acs.jctc.2c00424] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The fusion pore controls the release of exocytotic vesicle contents through a precise orchestration of lipids from the fusing membranes and proteins. There is a major lipid reorganization during the different stages in life of the fusion pore (membrane fusion, nucleation, and expansion) that can be scrutinized thermodynamically. In this work, using umbrella sampling simulations we describe the expansion of the fusion pore. We have calculated free energy profiles to drive a nascent, just nucleated, fusion pore to its expanded configuration. We have quantified the effects on the free energy of one and two Synaptotagmin-1 C2B domains in the cytosolic space. We show that C2B domains cumulatively reduce the cost for expansion, favoring the system to evolve toward full fusion. Finally, by conducting thousands of unbiased molecular dynamics simulations, we show that C2B domains significantly decrease the probability of kiss-and-run events.
Collapse
Affiliation(s)
- Ary Lautaro Di Bartolo
- 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), 5500 Mendoza, Argentina.,Facultad de Ciencias Exactas y Naturales, Universidad Nacional de Cuyo (UNCuyo), 5500 Mendoza, Argentina
| | - Claudia N Tomes
- 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), 5500 Mendoza, Argentina.,Facultad de Ciencias Exactas y Naturales, Universidad Nacional de Cuyo (UNCuyo), 5500 Mendoza, Argentina
| | - Luis S Mayorga
- 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), 5500 Mendoza, Argentina.,Facultad de Ciencias Exactas y Naturales, Universidad Nacional de Cuyo (UNCuyo), 5500 Mendoza, 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), 5500 Mendoza, Argentina.,Facultad de Ingeniería, Universidad Nacional de Cuyo (UNCuyo), 5500 Mendoza, Argentina
| |
Collapse
|
23
|
The lung surfactant activity probed with molecular dynamics simulations. Adv Colloid Interface Sci 2022; 304:102659. [PMID: 35421637 DOI: 10.1016/j.cis.2022.102659] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2021] [Revised: 03/18/2022] [Accepted: 03/31/2022] [Indexed: 01/17/2023]
Abstract
The surface of pulmonary alveolar subphase is covered with a mixture of lipids and proteins. This lung surfactant plays a crucial role in lung functioning. It shows a complex phase behavior which can be altered by the interaction with third molecules such as drugs or pollutants. For studying multicomponent biological systems, it is of interest to couple experimental approach with computational modelling yielding atomic-scale information. Simple two, three, or four-component model systems showed to be useful for getting more insight in the interaction between lipids, lipids and proteins or lipids and proteins with drugs and impurities. These systems were studied theoretically using molecular dynamic simulations and experimentally by means of the Langmuir technique. A better understanding of the structure and behavior of lung surfactants obtained from this research is relevant for developing new synthetic surfactants for efficient therapies, and may contribute to public health protection.
Collapse
|
24
|
Minkara MS, Josephson TR, Venteicher CL, Greenvall BR, Lindsey RK, Koenig PH, Siepmann JI. Nonane and Hexanol Adsorption in the Lamellar Phase of a Nonionic Surfactant: Molecular Simulations and Comparison to Ideal Adsorbed Solution Theory. J Phys Chem B 2022; 126:3940-3949. [PMID: 35594369 DOI: 10.1021/acs.jpcb.2c02871] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Adsorption of n-nonane/1-hexanol (C9/C6OH) mixtures into the lamellar phase formed by a 50/50 w/w triethylene glycol mono-n-decyl ether (C10E3)/water system was studied using configurational-bias Monte Carlo simulations in the osmotic Gibbs ensemble. The interactions were described by the Shinoda-Devane-Klein coarse-grained force field. Prior simulations probing single-component adsorption indicated that C9 molecules preferentially load near the center of the bilayer, increasing the bilayer thickness, whereas C6OH molecules are more likely to be found near the interface of the polar and nonpolar moieties, swelling the bilayer in the lateral dimension. Here, we extend this work to binary C9/C6OH adsorption to probe whether the difference in the spatial preferences may lead to a synergistic effect and enhanced loadings for the mixture. Comparing loading trends and the thermodynamics of binary adsorption to unary adsorption reveals that C9-C9 interactions lead to the largest enhancement, whereas C9-C6OH and C6OH-C6OH interactions are less favorable for this bilayer system. Ideal adsorbed solution theory yields satisfactory predictions of the binary loading.
Collapse
Affiliation(s)
- Mona S Minkara
- Department of Bioengineering, Northeastern University, 360 Huntington Avenue, Boston, Massachusetts 02115, United States.,Department of Chemistry and Chemical Theory Center, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| | - Tyler R Josephson
- Department of Chemistry and Chemical Theory Center, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States.,Department of Chemical, Biochemical, and Environmental Engineering, University of Maryland, Baltimore County, 1000 Hilltop Circle, Baltimore, Maryland 21250, United States
| | - Connor L Venteicher
- Department of Chemistry and Chemical Theory Center, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| | - Benjamin R Greenvall
- Department of Bioengineering, Northeastern University, 360 Huntington Avenue, Boston, Massachusetts 02115, United States
| | - Rebecca K Lindsey
- Department of Chemistry and Chemical Theory Center, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| | - Peter H Koenig
- Computational Chemistry, Modeling and Simulation, The Procter and Gamble Company, 8256 Union Centre Blvd, West Chester, Ohio 45069, United States
| | - J Ilja Siepmann
- Department of Chemistry and Chemical Theory Center, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States.,Department of Chemical Engineering and Materials Science, University of Minnesota, 421 Washington Avenue SE, Minneapolis, Minnesota 55455-0132, United States
| |
Collapse
|
25
|
Kawamoto S, Liu H, Miyazaki Y, Seo S, Dixit M, DeVane R, MacDermaid C, Fiorin G, Klein ML, Shinoda W. SPICA Force Field for Proteins and Peptides. J Chem Theory Comput 2022; 18:3204-3217. [PMID: 35413197 DOI: 10.1021/acs.jctc.1c01207] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A coarse-grained (CG) model for peptides and proteins was developed as an extension of the Surface Property fItting Coarse grAined (SPICA) force field (FF). The model was designed to examine membrane proteins that are fully compatible with the lipid membranes of the SPICA FF. A preliminary version of this protein model was created using thermodynamic properties, including the surface tension and density in the SPICA (formerly called SDK) FF. In this study, we improved the CG protein model to facilitate molecular dynamics (MD) simulations with a reproduction of multiple properties from both experiments and all-atom (AA) simulations. An elastic network model was adopted to maintain the secondary structure within a single chain. The side-chain analogues reproduced the transfer free energy profiles across the lipid membrane and demonstrated reasonable association free energy (potential of mean force) in water compared to those from AA MD. A series of peptides/proteins adsorbed onto or penetrated into the membrane simulated by the CG MD correctly predicted the penetration depths and tilt angles of peripheral and transmembrane peptides/proteins as comparable to those in the orientations of proteins in membranes (OPM) database. In addition, the dimerization free energies of several transmembrane helices within a lipid bilayer were comparable to those from experimental estimation. Application studies on a series of membrane protein assemblies, scramblases, and poliovirus capsids demonstrated the good performance of the SPICA FF.
Collapse
Affiliation(s)
- Shuhei Kawamoto
- Department of Materials Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan
| | - Huihui Liu
- Department of Materials Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan
| | - Yusuke Miyazaki
- Department of Materials Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan.,Research Institute for Interdisciplinary Science, Okayama University, 3-1-1 Tsushima-naka, Kita-ku, Okayama 700-8530, Japan
| | - Sangjae Seo
- Department of Materials Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan.,Korea Institute of Science and Technology Information, 245 Daehak-ro, Daejeon 34141, Republic of Korea
| | - Mayank Dixit
- Department of Materials Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan
| | - Russell DeVane
- Modeling & Simulation, Corporate Research & Development, The Procter and Gamble Company, West Chester, Ohio 45069, United States
| | - Christopher MacDermaid
- Institute for Computational Molecular Science, Temple University, 1925 North 12th Street, Philadelphia, Pennsylvania 19122, United States
| | - Giacomo Fiorin
- Institute for Computational Molecular Science, Temple University, 1925 North 12th Street, Philadelphia, Pennsylvania 19122, United States
| | - Michael L Klein
- Institute for Computational Molecular Science, Temple University, 1925 North 12th Street, Philadelphia, Pennsylvania 19122, United States
| | - Wataru Shinoda
- Department of Materials Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan.,Research Institute for Interdisciplinary Science, Okayama University, 3-1-1 Tsushima-naka, Kita-ku, Okayama 700-8530, Japan.,Department of Chemistry, Faculty of Science, Okayama University, 3-1-1 Tsushima-naka, Kita-ku, Okayama 700-8530, Japan
| |
Collapse
|
26
|
Lu Y, Allegri G, Huskens J. Vesicle-based artificial cells: materials, construction methods and applications. MATERIALS HORIZONS 2022; 9:892-907. [PMID: 34908080 PMCID: PMC8900604 DOI: 10.1039/d1mh01431e] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/03/2021] [Accepted: 11/19/2021] [Indexed: 05/27/2023]
Abstract
The construction of artificial cells with specific cell-mimicking functions helps to explore complex biological processes and cell functions in natural cell systems and provides an insight into the origins of life. Bottom-up methods are widely used for engineering artificial cells based on vesicles by the in vitro assembly of biomimetic materials. In this review, the design of artificial cells with a specific function is discussed, by considering the selection of synthetic materials and construction technologies. First, a range of biomimetic materials for artificial cells is reviewed, including lipid, polymeric and hybrid lipid/copolymer materials. Biomaterials extracted from natural cells are also covered in this part. Then, the formation of microscale, giant unilamellar vesicles (GUVs) is reviewed based on different technologies, including gentle hydration, electro-formation, phase transfer and microfluidic methods. Subsequently, applications of artificial cells based on single vesicles or vesicle networks are addressed for mimicking cell behaviors and signaling processes. Microreactors for synthetic biology and cell-cell communication are highlighted here as well. Finally, current challenges and future trends for the development and applications of artificial cells are described.
Collapse
Affiliation(s)
- Yao Lu
- Molecular NanoFabrication Group, Department of Molecules & Materials, MESA+ Institute, Faculty of Science and Technology, University of Twente, P.O. Box 217, 7500 AE, Enschede, The Netherlands.
| | - Giulia Allegri
- Molecular NanoFabrication Group, Department of Molecules & Materials, MESA+ Institute, Faculty of Science and Technology, University of Twente, P.O. Box 217, 7500 AE, Enschede, The Netherlands.
| | - Jurriaan Huskens
- Molecular NanoFabrication Group, Department of Molecules & Materials, MESA+ Institute, Faculty of Science and Technology, University of Twente, P.O. Box 217, 7500 AE, Enschede, The Netherlands.
| |
Collapse
|
27
|
Kashyap HK. Deciphering Ethanol-Driven Swelling, Rupturing, Aggregation, and Fusion of Lipid Vesicles Using Coarse-Grained Molecular Dynamics Simulations. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2022; 38:2445-2459. [PMID: 35167280 DOI: 10.1021/acs.langmuir.1c02763] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Traditionally, liquid ethanol is known to enhance the permeability of lipid membranes and causes vesicle aggregation and fusion. However, how the amphiphilic ethanol molecules perturb the lipid vesicles to facilitate their aggregation or fusion has not been addressed at any level of molecular simulations. Herein, not only have we developed a coarse-grained (CG) model for liquid ethanol, its aqueous mixture, and hydrated lipid membranes for molecular dynamics (MD) simulations, but also utilized it to delineate the aggregation and fusion of lipid vesicles using CG-MD simulations with multimillion particles. We have systematically parametrized the force-field for pure ethanol and its interactions with hydrated POPC and POPE model lipid membranes. In this process, we have successfully reproduced the bulk ethanol structure and concentration-dependent density of aqueous ethanol. To quantify the interaction of ethanol with lipid membranes, we have reproduced the transfer free energy of the ethanol molecule across the hydrated bilayers, and the concentration-dependent distribution of ethanol molecules across the lipid bilayers. After having acceptable force-field parameters for ethanol-membrane interactions, we have checked the effect of ethanol toward the vesicles comprising POPC lipids. We observe a rapid increase in the size of the POPC lipid vesicles with increasing amounts of ethanol up to 30 mol %. We unambiguously observe swelling and decrease in the thickness of the POPC vesicles with increasing amounts of ethanol up to 30 mol %, beyond which the vesicles begin to lose their integrity and rupture at higher mol % of ethanol. The fusion study of two vesicles demonstrates that fused vesicles can be obtained from 20 to 30 mol % of ethanol provided that they are brought closer than a critical distance at a particular mol %. The multivesicle simulations show that along with the increase in the sizes of vesicles the propensity of vesicle aggregation increases as the mol % of ethanol increases.
Collapse
Affiliation(s)
- Hemant K Kashyap
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| |
Collapse
|
28
|
Nandi S, Pyne A, Layek S, Arora C, Sarkar N. The Dietary Nutrient Trimethylamine N-Oxide Affects the Phospholipid Vesicle Membrane: Probable Route to Adverse Intake. J Phys Chem Lett 2021; 12:12411-12418. [PMID: 34939822 DOI: 10.1021/acs.jpclett.1c03201] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Trimethylamine N-oxide (TMAO), a choline-containing dietary supplement obtained from red meat, egg, and other animal resources, on excess accumulation is known to cause cardiovascular diseases (CVDs) like atherosclerosis. To understand the molecular mechanism of the pathogenesis of TMAO-induced CVDs, we have set up 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) membrane in water that mimicked the endothelial cell membrane-blood interface of the artery wall and investigated the effect of an elevated concentration of TMAO on the membrane. We found that TMAO exerts an "action at a distance" mechanism through electrostatic force of attraction that significantly alters various properties of the membrane, like hydrophobicity, lateral organization, and interfacial water dynamics, which elevates the rigidity of the membrane. Such an effect was found to be further amplified in the presence of known causes of CVDs, i.e., high content of cholesterol (Chol). Therefore, TMAO-induced membrane rigidity may restrict the intrinsic elasticity of an artery membrane, expected to be introducing "hardening of the arteries", which makes the membrane atherosclerotic.
Collapse
Affiliation(s)
- Sourav Nandi
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur, WB 721302, India
| | - Arghajit Pyne
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur, WB 721302, India
| | - Souvik Layek
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur, WB 721302, India
| | - Chirag Arora
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur, WB 721302, India
| | - Nilmoni Sarkar
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur, WB 721302, India
| |
Collapse
|
29
|
Benazieb O, Loison C, Thalmann F. Rheology of sliding leaflets in coarse-grained DSPC lipid bilayers. Phys Rev E 2021; 104:054802. [PMID: 34942802 DOI: 10.1103/physreve.104.054802] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2021] [Accepted: 10/06/2021] [Indexed: 11/07/2022]
Abstract
Amphiphilic lipid bilayers modify the friction properties of the surfaces on top of which they are deposited. In particular, the measured sliding friction coefficient can be significantly reduced compared with the native surface. We investigate in this work the friction properties of a numerical coarse-grained model of DSPC (1,2-distearoyl-sn-glycero-3-phosphocholine) lipid bilayer subject to longitudinal shear. The interleaflet friction coefficient is obtained from out-of-equilibrium pulling or from relaxation simulations. In particular, we gain access to the transient viscoelastic response of a sheared bilayer. The bilayer mechanical response is found to depend significantly on the membrane physical state, with evidence in favor of a linear response regime in the fluid but not in the gel region. The linear response validity domain is established, and the timescales appearing in the membrane response discussed.
Collapse
Affiliation(s)
- Othmene Benazieb
- Institut Charles Sadron, CNRS and University of Strasbourg, 23 rue du Loess, F-67034 Strasbourg cedex 2, France
| | - Claire Loison
- Institut Lumière Matière, UMR5306 Université Lyon 1-CNRS, Université de Lyon, 69622 Villeurbanne cedex, France
| | - Fabrice Thalmann
- Institut Charles Sadron, CNRS and University of Strasbourg, 23 rue du Loess, F-67034 Strasbourg cedex 2, France
| |
Collapse
|
30
|
Zhu Y, Bai X, Hu G. Interfacial behavior of phospholipid monolayers revealed by mesoscopic simulation. Biophys J 2021; 120:4751-4762. [PMID: 34562445 DOI: 10.1016/j.bpj.2021.09.031] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2021] [Revised: 08/28/2021] [Accepted: 09/20/2021] [Indexed: 12/01/2022] Open
Abstract
A mesoscopic model with molecular resolution is presented for dipalmitoyl phosphatidylcholine (DPPC) and palmitoyl oleoyl phosphatidylcholine (POPC) monolayer simulations at the air-water interface using many-body dissipative particle dynamics (MDPD). The parameterization scheme is rigorously based on reproducing the physical properties of water and alkane and the interfacial property of the phospholipid monolayer by comparison with experimental results. Using much less computing cost, these MDPD simulations yield a similar surface pressure-area isotherm as well as similar pressure-related morphologies as all-atom simulations and experiments. Moreover, the compressibility modulus, order parameter of lipid tails, and thickness of the phospholipid monolayer are quantitatively in line with the all-atom simulations and experiments. This model also captures the sensitive changes in the pressure-area isotherms of mixed DPPC/POPC monolayers with altered mixing ratios, indicating that the model is promising for applications with complex natural phospholipid monolayers. These results demonstrate a significant improvement of quantitative phospholipid monolayer simulations over previous coarse-grained models.
Collapse
Affiliation(s)
- Yongzheng Zhu
- Department of Engineering Mechanics, State Key Laboratory of Fluid Power and Mechatronic Systems, Zhejiang University, Hangzhou, China; The State Key Laboratory of Nonlinear Mechanics, Institute of Mechanics, Chinese Academy of Sciences, Beijing, China; School of Engineering Science, University of Chinese Academy of Sciences, Beijing, China
| | - Xuan Bai
- Department of Engineering Mechanics, State Key Laboratory of Fluid Power and Mechatronic Systems, Zhejiang University, Hangzhou, China
| | - Guoqing Hu
- Department of Engineering Mechanics, State Key Laboratory of Fluid Power and Mechatronic Systems, Zhejiang University, Hangzhou, China.
| |
Collapse
|
31
|
Ugarte La Torre D, Takada S. Modeling lipid-protein interactions for coarse-grained lipid and Cα protein models. J Chem Phys 2021; 155:155101. [PMID: 34686048 DOI: 10.1063/5.0057278] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Biological membranes that play major roles in diverse functions are composed of numerous lipids and proteins, making them an important target for coarse-grained (CG) molecular dynamics (MD) simulations. Recently, we have developed the CG implicit solvent lipid force field (iSoLF) that has a resolution compatible with the widely used Cα protein representation [D. Ugarte La Torre and S. Takada, J. Chem. Phys. 153, 205101 (2020)]. In this study, we extended it and developed a lipid-protein interaction model that allows the combination of the iSoLF and the Cα protein force field, AICG2+. The hydrophobic-hydrophilic interaction is modeled as a modified Lennard-Jones potential in which parameters were tuned partly to reproduce the experimental transfer free energy and partly based on the free energy profile normal to the membrane surface from previous all-atom MD simulations. Then, the obtained lipid-protein interaction is tested for the configuration and placement of transmembrane proteins, water-soluble proteins, and peripheral proteins, showing good agreement with prior knowledge. The interaction is generally applicable and is implemented in the publicly available software, CafeMol.
Collapse
Affiliation(s)
- Diego Ugarte La Torre
- Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto, Japan
| | - Shoji Takada
- Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto, Japan
| |
Collapse
|
32
|
Yadav HOS, Harada S, Kuo AT, Urata S, Shinoda W. Hemimicelle formation of semi-fluorocarbon chains at air–water interface: coarse-grained molecular dynamics study with an extension of the SPICA force field. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1910355] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Affiliation(s)
- Hari O. S. Yadav
- Department of Materials Chemistry, Nagoya University, Nagoya, Japan
| | - Shogo Harada
- Department of Materials Chemistry, Nagoya University, Nagoya, Japan
| | - An-Tsung Kuo
- Innovative Technology Laboratories, AGC Inc., Yokohama, Japan
| | - Shingo Urata
- Innovative Technology Laboratories, AGC Inc., Yokohama, Japan
| | - Wataru Shinoda
- Department of Materials Chemistry, Nagoya University, Nagoya, Japan
| |
Collapse
|
33
|
Zhang W, Zhang MY, Wang K, Sun R, Zhao S, Zhang Z, He YP, Yu F. Geometry transformation of ionic surfactants and adsorption behavior on water/ n-decane-interface: calculation by molecular dynamics simulation and DFT study. RSC Adv 2021; 11:28286-28294. [PMID: 35480765 PMCID: PMC9038023 DOI: 10.1039/d1ra04669a] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2021] [Accepted: 08/15/2021] [Indexed: 01/20/2023] Open
Abstract
Understanding the effect of surfactant structure on their ability to modify interfacial properties is of great scientific and industrial interest. In this work, we have synthesized four amide based ionic surfactants under acidic or basic conditions, including CTHA·HCl, CTEA·HCl, CTHA−Na+ and CTEA−Na+. Experiments have proved that the anionic surfactant with polyethylene oxide groups (CTEA−Na+) had the lowest surface tension on the water/n-decane interface. Molecular dynamics simulations have been applied to investigate the structural effect on the adsorption behavior of four different surfactants. The surface tension, interface thickness, interface formation energy, density profiles, order parameters, radial distribution function on the water/n-decane interfaces were calculated and compared. During the equilibrium states, we found that the interface configuration of two cationic surfactants are almost linear while the two anionic surfactants are changed to bending shapes due to the different positions of the hydrophilic head groups. Further DFT study and wavefunction analysis of surfactants have shown that CTEA−Na+ can form stronger vdW interactions with n-decane molecules due to a more neutral electrostatic potential distribution. Meanwhile, the introduction of polyethylene oxide groups has offered more H-bonding sites and resulted in more concentrated H-bonding interactions with water molecules. The difference of weak interactions may contribute to the conformational change and finally affect the interface properties of these ionic surfactants. The difference of weak interactions may contribute to the conformational change and finally affect the interface properties of these ionic surfactants.![]()
Collapse
Affiliation(s)
- Wannian Zhang
- Key Laboratory for Functional Material, Educational Department of Liaoning Province, School of Chemical Engineering, University of Science and Technology Liaoning Anshan 114051 P. R. China .,State Key Laboratory of Fine Chemicals, Ningbo Institute of Dalian University of Technology No. 26 Yucai Road, Jiangbei District Ningbo 315016 P. R. China .,Key Laboratory of Petrochemical Catalytic Science and Technology, Liaoning Petrochemical University Dandong Lu West 1 Fushun 113001 Liaoning P. R. China +86-2456860548
| | - Ming-Yuan Zhang
- Key Laboratory of Petrochemical Catalytic Science and Technology, Liaoning Petrochemical University Dandong Lu West 1 Fushun 113001 Liaoning P. R. China +86-2456860548
| | - Kai Wang
- State Key Laboratory of Fine Chemicals, Ningbo Institute of Dalian University of Technology No. 26 Yucai Road, Jiangbei District Ningbo 315016 P. R. China .,Key Laboratory of Petrochemical Catalytic Science and Technology, Liaoning Petrochemical University Dandong Lu West 1 Fushun 113001 Liaoning P. R. China +86-2456860548
| | - Ruixia Sun
- Key Laboratory of Petrochemical Catalytic Science and Technology, Liaoning Petrochemical University Dandong Lu West 1 Fushun 113001 Liaoning P. R. China +86-2456860548
| | - Shanlin Zhao
- Key Laboratory for Functional Material, Educational Department of Liaoning Province, School of Chemical Engineering, University of Science and Technology Liaoning Anshan 114051 P. R. China .,Key Laboratory of Petrochemical Catalytic Science and Technology, Liaoning Petrochemical University Dandong Lu West 1 Fushun 113001 Liaoning P. R. China +86-2456860548
| | - Zhiqiang Zhang
- Key Laboratory for Functional Material, Educational Department of Liaoning Province, School of Chemical Engineering, University of Science and Technology Liaoning Anshan 114051 P. R. China
| | - Yu-Peng He
- Key Laboratory for Functional Material, Educational Department of Liaoning Province, School of Chemical Engineering, University of Science and Technology Liaoning Anshan 114051 P. R. China .,State Key Laboratory of Fine Chemicals, Ningbo Institute of Dalian University of Technology No. 26 Yucai Road, Jiangbei District Ningbo 315016 P. R. China .,Key Laboratory of Petrochemical Catalytic Science and Technology, Liaoning Petrochemical University Dandong Lu West 1 Fushun 113001 Liaoning P. R. China +86-2456860548
| | - Fang Yu
- State Key Laboratory of Fine Chemicals, Ningbo Institute of Dalian University of Technology No. 26 Yucai Road, Jiangbei District Ningbo 315016 P. R. China .,Key Laboratory of Petrochemical Catalytic Science and Technology, Liaoning Petrochemical University Dandong Lu West 1 Fushun 113001 Liaoning P. R. China +86-2456860548
| |
Collapse
|
34
|
Griffiths MZ, Shinoda W. Analyzing the Role of Surfactants in the Colloidal Stability of Nanoparticles in Oil through Coarse-Grained Molecular Dynamics Simulations. J Phys Chem B 2021; 125:6315-6321. [DOI: 10.1021/acs.jpcb.1c01148] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
- Mark Z. Griffiths
- 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
|
35
|
Zoni V, Khaddaj R, Lukmantara I, Shinoda W, Yang H, Schneiter R, Vanni S. Seipin accumulates and traps diacylglycerols and triglycerides in its ring-like structure. Proc Natl Acad Sci U S A 2021; 118:e2017205118. [PMID: 33674387 PMCID: PMC7958289 DOI: 10.1073/pnas.2017205118] [Citation(s) in RCA: 51] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
Lipid droplets (LDs) are intracellular organelles responsible for lipid storage, and they emerge from the endoplasmic reticulum (ER) upon the accumulation of neutral lipids, mostly triglycerides (TG), between the two leaflets of the ER membrane. LD biogenesis takes place at ER sites that are marked by the protein seipin, which subsequently recruits additional proteins to catalyze LD formation. Deletion of seipin, however, does not abolish LD biogenesis, and its precise role in controlling LD assembly remains unclear. Here, we use molecular dynamics simulations to investigate the molecular mechanism through which seipin promotes LD formation. We find that seipin clusters TG, as well as its precursor diacylglycerol, inside its unconventional ring-like oligomeric structure and that both its luminal and transmembrane regions contribute to this process. This mechanism is abolished upon mutations of polar residues involved in protein-TG interactions into hydrophobic residues. Our results suggest that seipin remodels the membrane of specific ER sites to prime them for LD biogenesis.
Collapse
Affiliation(s)
- Valeria Zoni
- Department of Biology, University of Fribourg, 1700 Fribourg, Switzerland
| | - Rasha Khaddaj
- Department of Biology, University of Fribourg, 1700 Fribourg, Switzerland
| | - Ivan Lukmantara
- School of Biotechnology and Biomolecular Sciences, University of New South Wales, Sydney, NSW 2052, Australia
| | - Wataru Shinoda
- Department of Materials Chemistry, Nagoya University, Chikusa-ku, 464-8603 Nagoya, Japan
| | - Hongyuan Yang
- School of Biotechnology and Biomolecular Sciences, University of New South Wales, Sydney, NSW 2052, Australia
| | - Roger Schneiter
- Department of Biology, University of Fribourg, 1700 Fribourg, Switzerland
| | - Stefano Vanni
- Department of Biology, University of Fribourg, 1700 Fribourg, Switzerland;
| |
Collapse
|
36
|
Hall KW, Percec S, Shinoda W, Klein ML. Chain-End Modification: A Starting Point for Controlling Polymer Crystal Nucleation. Macromolecules 2021. [DOI: 10.1021/acs.macromol.0c02398] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Kyle Wm. Hall
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
- Institute for Computational Molecular Science, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Simona Percec
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Wataru Shinoda
- Department of Materials Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan
| | - Michael L. Klein
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
- Institute for Computational Molecular Science, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
37
|
Wan M, Song J, Yang Y, Gao L, Fang W. Development of coarse-grained force field for alcohols: an efficient meta-multilinear interpolation parameterization algorithm. Phys Chem Chem Phys 2021; 23:1956-1966. [PMID: 33464253 DOI: 10.1039/d0cp05503d] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Coarse-grained (CG) molecular dynamics are powerful tools to access a mesoscopic phenomenon and simultaneously record microscopic details, but currently the CG force fields (FFs) are still limited by low parameterization efficiency and poor accuracy especially for polar molecules. In this work, we developed a Meta-Multilinear Interpolation Parameterization (Meta-MIP) algorithm to optimize the CG FFs for alcohols. This algorithm significantly boosts parameterization efficiency by constructing on-the-fly local databases to cover the global optimal parameterization path. In specific, an alcohol molecule is mapped to a heterologous model composed of an OH bead and a hydrocarbon portion which consists of alkane beads representing two to four carbon atoms. Non-bonded potentials are described by soft Morse functions that have no tail-corrections but can still retain good continuities at truncation distance. Nearly all of the properties in terms of density, heat of vaporization, surface tension, and solvation free energy for alcohols predicted by the current FFs deviate from experimental values by less than 7%. This Meta-MIP algorithm can be readily applied to force field development for a wide variety of molecules or functional groups, in many situations including but not limited to CG FFs.
Collapse
Affiliation(s)
- Mingwei Wan
- Institution of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | | | | | | | | |
Collapse
|
38
|
Marchetto A, Si Chaib Z, Rossi CA, Ribeiro R, Pantano S, Rossetti G, Giorgetti A. CGMD Platform: Integrated Web Servers for the Preparation, Running, and Analysis of Coarse-Grained Molecular Dynamics Simulations. Molecules 2020; 25:E5934. [PMID: 33333836 PMCID: PMC7765266 DOI: 10.3390/molecules25245934] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2020] [Revised: 11/21/2020] [Accepted: 11/28/2020] [Indexed: 12/17/2022] Open
Abstract
Advances in coarse-grained molecular dynamics (CGMD) simulations have extended the use of computational studies on biological macromolecules and their complexes, as well as the interactions of membrane protein and lipid complexes at a reduced level of representation, allowing longer and larger molecular dynamics simulations. Here, we present a computational platform dedicated to the preparation, running, and analysis of CGMD simulations. The platform is built on a completely revisited version of our Martini coarsE gRained MembrAne proteIn Dynamics (MERMAID) web server, and it integrates this with other three dedicated services. In its current version, the platform expands the existing implementation of the Martini force field for membrane proteins to also allow the simulation of soluble proteins using the Martini and the SIRAH force fields. Moreover, it offers an automated protocol for carrying out the backmapping of the coarse-grained description of the system into an atomistic one.
Collapse
Affiliation(s)
- Alessandro Marchetto
- Department of Biotechnology, University of Verona, 37134 Verona, Italy; (A.M.); (C.A.R.); (R.R.)
- Computational Biomedicine, Institute for Neuroscience and Medicine (INM-9) and Institute for Advanced Simulations (IAS-5), Forschungszentrum Jülich, 52425 Jülich, Germany;
- Faculty of Mathematics, Computer Science and Natural Sciences, RWTH Aachen University, 52062 Aachen, Germany
| | - Zeineb Si Chaib
- Computational Biomedicine, Institute for Neuroscience and Medicine (INM-9) and Institute for Advanced Simulations (IAS-5), Forschungszentrum Jülich, 52425 Jülich, Germany;
- Faculty of Mathematics, Computer Science and Natural Sciences, RWTH Aachen University, 52062 Aachen, Germany
| | - Carlo Alberto Rossi
- Department of Biotechnology, University of Verona, 37134 Verona, Italy; (A.M.); (C.A.R.); (R.R.)
| | - Rui Ribeiro
- Department of Biotechnology, University of Verona, 37134 Verona, Italy; (A.M.); (C.A.R.); (R.R.)
| | - Sergio Pantano
- Institut Pasteur de Montevideo, Mataojo 2020, 11400 Montevideo, Uruguay;
| | - Giulia Rossetti
- Computational Biomedicine, Institute for Neuroscience and Medicine (INM-9) and Institute for Advanced Simulations (IAS-5), Forschungszentrum Jülich, 52425 Jülich, Germany;
- Jülich Supercomputing Center (JSC), Forschungszentrum Jülich, 52428 Jülich, Germany
- Department of Hematology, Oncology, Hemostaseology, and Stem Cell Transplantation, RWTH Aachen University, 52062 Aachen, Germany
| | - Alejandro Giorgetti
- Department of Biotechnology, University of Verona, 37134 Verona, Italy; (A.M.); (C.A.R.); (R.R.)
- Computational Biomedicine, Institute for Neuroscience and Medicine (INM-9) and Institute for Advanced Simulations (IAS-5), Forschungszentrum Jülich, 52425 Jülich, Germany;
| |
Collapse
|
39
|
Wang X, Santo KP, Neimark AV. Modeling Gas-Liquid Interfaces by Dissipative Particle Dynamics: Adsorption and Surface Tension of Cetyl Trimethyl Ammonium Bromide at the Air-Water Interface. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2020; 36:14686-14698. [PMID: 33216560 DOI: 10.1021/acs.langmuir.0c02572] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
Adsorption of surfactants at gas-liquid interfaces that causes reduction in the surface tension is a classical problem in colloid and interface science with multiple practical applications in oil and gas recovery, separations, cosmetics, personal care, and biomedicine. Here, we develop an original coarse-grained model of the liquid-gas interface within the conventional dissipative particle dynamics (DPD) framework with the goal of quantitatively predicting the surface tension in the presence of surfactants. As a practical case-study example, we explore the adsorption of the cationic surfactant cetyl trimethyl ammonium bromide (CTAB) on the air-water interface. The gas phase is modeled as a DPD fluid composed of fictitious hard-core "gas" beads with exponentially decaying repulsive potentials to prevent penetration of the liquid phase components. A rigorous parametrization scheme is proposed based on matching the bulk and interfacial properties of water and octane taken as the reference compounds. Quantitative agreement between the simulated and experimental surface tension of CTAB solutions is found for a wide range of bulk surfactant concentrations (∼10-3 to ∼1 mmol/L) with the reduction of the surface tension from ∼72 mN/m (pure water) to the limiting value of ∼37.5 mN/m at the critical micelle concentration. The gas phase DPD model with the proposed parametrization scheme can be extended and applied to modeling various gas-liquid interfaces with surfactant and lipid monolayers, such as bubble suspensions, foams, froths, etc.
Collapse
Affiliation(s)
- Xinyang Wang
- Chemical and Biochemical Engineering Department, Rutgers, The State University of New Jersey, Piscataway, New Jersey 08854, United States
- School of Resources and Civil Engineering, Northeastern University, Shenyang, Liaoning 110819, China
| | - Kolattukudy P Santo
- Chemical and Biochemical Engineering Department, Rutgers, The State University of New Jersey, Piscataway, New Jersey 08854, United States
| | - Alexander V Neimark
- Chemical and Biochemical Engineering Department, Rutgers, The State University of New Jersey, Piscataway, New Jersey 08854, United States
| |
Collapse
|
40
|
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
|
41
|
Gao P, Yang X, Tartakovsky AM. Learning Coarse-Grained Potentials for Binary Fluids. J Chem Inf Model 2020; 60:3731-3745. [PMID: 32668158 DOI: 10.1021/acs.jcim.0c00337] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
For a multiple-fluid system, CG models capable of accurately predicting the interfacial properties as a function of curvature are still lacking. In this work, we propose a new probabilistic machine learning (ML) model for learning CG potentials for binary fluids. The water-hexane mixture is selected as a typical immiscible binary liquid-liquid system. We develop a new CG force field (FF) using the Shinoda-DeVane-Klein (SDK) FF framework and compute parameters in this CG FF using the proposed probabilistic ML method. It is shown that a standard response-surface approach does not provide a unique set of parameters, as it results in a loss function with multiple shallow minima. To address this challenge, we develop a probabilistic ML approach where we compute the probability density function (PDF) of parameters that minimize the loss function. The PDF has a well-defined peak corresponding to a unique set of parameters in the CG FF that reproduces the desired properties of a liquid-liquid interface. We compare the performance of the new CG FF with several existing FFs for the water-hexane mixture, including two atomistic and three CG FFs with respect to modeling the interface structure and thermodynamic properties. It is demonstrated that the new FF significantly improves the CG model prediction of both the interfacial tension and structure for the water-hexane mixture.
Collapse
Affiliation(s)
- Peiyuan Gao
- Advanced Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Xiu Yang
- Department of Industrial and Systems Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Alexandre M Tartakovsky
- Advanced Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| |
Collapse
|
42
|
Shobhna, Kumari M, Kashyap HK. A coarse-grained model of dimethyl sulfoxide for molecular dynamics simulations with lipid membranes. J Chem Phys 2020; 153:035104. [DOI: 10.1063/5.0014614] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Affiliation(s)
- Shobhna
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| | - Monika Kumari
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| | - Hemant K. Kashyap
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| |
Collapse
|
43
|
Seo S, Murata M, Shinoda W. Pivotal Role of Interdigitation in Interleaflet Interactions: Implications from Molecular Dynamics Simulations. J Phys Chem Lett 2020; 11:5171-5176. [PMID: 32515980 DOI: 10.1021/acs.jpclett.0c01317] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
The asymmetric lipid composition in plasma membranes within the inner leaflet is not typically suitable for domain formation. Thus elucidation of the likelihood of the formation or stability of a raft-like domain in the inner leaflet is necessary. Herein we investigated the phase behavior of asymmetric membranes using coarse-grained molecular dynamics simulations. The lipid leaflet comprising dioleoylphosphatidylcholine (DOPC) and cholesterol (Chol) does not typically show well-developed domains in symmetric bilayer membranes; however, it does separate into liquid ordered (Lo) and liquid disordered (Ld) phases when the opposing leaflet containing sphingomyelin (SM), DOPC, and Chol demonstrates domain formation. We determine that interdigitated acyl chains modulated the partitioning of Chol in the opposing leaflet, resulting in phase separation. Similarly, the acyl chain length of SM within the opposing leaflet affected the phase behavior of the leaflet. Our results reveal the crucial role of interdigitation in determining the phase status in asymmetric membranes.
Collapse
Affiliation(s)
- Sangjae Seo
- Department of Materials Chemistry, Nagoya University, Nagoya 464-8603, Japan
- Korean Institute of Science and Technology Information, Daejeon 34141, Republic of Korea
| | - Michio Murata
- Graduate School of Science, Osaka University, Toyonaka, Osaka 560-0043, Japan
| | - Wataru Shinoda
- Department of Materials Chemistry, Nagoya University, Nagoya 464-8603, Japan
| |
Collapse
|
44
|
Cellular absorption of small molecules: free energy landscapes of melatonin binding at phospholipid membranes. Sci Rep 2020; 10:9235. [PMID: 32513935 PMCID: PMC7280225 DOI: 10.1038/s41598-020-65753-z] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2019] [Accepted: 05/05/2020] [Indexed: 12/28/2022] Open
Abstract
Free energy calculations are essential to unveil mechanisms at the atomic scale such as binding of small solutes and their translocation across cell membranes, eventually producing cellular absorption. Melatonin regulates biological rhythms and is directly related to carcinogenesis and neurodegenerative disorders. Free energy landscapes obtained from well-tempered metadynamics simulations precisely describe the characteristics of melatonin binding to specific sites in the membrane and reveal the role of cholesterol in free energy barrier crossing. A specific molecular torsional angle and the distance between melatonin and the center of the membrane along the normal to the membrane Z-axis have been considered as suitable reaction coordinates. Free energy barriers between two particular orientations of the molecular structure (folded and extended) have been found to be of about 18 kJ/mol for z-distances of about 1–2 nm. The ability of cholesterol to expel melatonin out of the internal regions of the membrane towards the interface and the external solvent is explained from a free energy perspective. The calculations reported here offer detailed free energy landscapes of melatonin embedded in model cell membranes and reveal microscopic information on its transition between free energy minima, including the location of relevant transition states, and provide clues on the role of cholesterol in the cellular absorption of small molecules.
Collapse
|
45
|
Damre M, Marchetto A, Giorgetti A. MERMAID: dedicated web server to prepare and run coarse-grained membrane protein dynamics. Nucleic Acids Res 2020; 47:W456-W461. [PMID: 31106328 PMCID: PMC6602572 DOI: 10.1093/nar/gkz416] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2019] [Revised: 04/26/2019] [Accepted: 05/05/2019] [Indexed: 01/06/2023] Open
Abstract
Atomistic molecular dynamics simulations of membrane proteins have been shown to be extremely useful for characterizing the molecular features underlying their function, but require high computational power, limiting the understanding of complex events in membrane proteins, e.g. ion channels gating, GPCRs activation. To overcome this issue, it has been shown that coarse-grained approaches, although requiring less computational power, are still capable of correctly describing molecular events underlying big conformational changes in biological systems. Here, we present the Martini coarse-grained membrane protein dynamics (MERMAID), a publicly available web interface that allows the user to prepare and run coarse-grained molecular dynamics (CGMD) simulations and to analyse the trajectories.
Collapse
Affiliation(s)
- Mangesh Damre
- Neurobiology, International School for Advanced Studies (SISSA), via Bonomea 265, 34136 Trieste, Italy.,Department of Biotechnology, University of Verona, Ca Vignal 1, strada Le Grazie, 15, 37134 Verona, Italy
| | - Alessandro Marchetto
- Department of Biotechnology, University of Verona, Ca Vignal 1, strada Le Grazie, 15, 37134 Verona, Italy
| | - Alejandro Giorgetti
- Department of Biotechnology, University of Verona, Ca Vignal 1, strada Le Grazie, 15, 37134 Verona, Italy.,Institute for Advanced Simulations (IAS)-5/Institute for Neuroscience and Medicine (INM)-9, Forschungszentrum Jülich, 52428 Jülich, Germany
| |
Collapse
|
46
|
MacDermaid CM, Hall KW, DeVane RH, Klein ML, Fiorin G. Coexistence of Lipid Phases Stabilizes Interstitial Water in the Outer Layer of Mammalian Skin. Biophys J 2020; 118:1588-1601. [PMID: 32101711 DOI: 10.1016/j.bpj.2020.01.044] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2019] [Revised: 01/20/2020] [Accepted: 01/27/2020] [Indexed: 12/22/2022] Open
Abstract
The lipid matrix in the outer layer of mammalian skin, the stratum corneum, has been previously investigated by multiple biophysical techniques aimed at identifying hydrophilic and lipophilic pathways of permeation. Although consensus is developing over the microscopic structure of the lipid matrix, no molecular-resolution model describes the permeability of all chemical species simultaneously. Using molecular dynamics simulations of a model mixture of skin lipids, the self-assembly of the lipid matrix lamellae has been studied. At higher humidity, the resulting lamellar phase is maintained by partitioning excess water into isolated droplets of controlled size and spatial distribution. The droplets may fuse together to form intralamellar water channels, thereby providing a pathway for the permeation of hydrophilic species. These results reconcile competing data on the outer skin's structure and broaden the scope of molecular-based methods to improve the safety of topical products and to advance transdermal drug delivery.
Collapse
Affiliation(s)
- Christopher M MacDermaid
- Institute for Computational Molecular Science and Temple Materials Institute, Philadelphia, Pennsylvania
| | - Kyle Wm Hall
- Institute for Computational Molecular Science and Temple Materials Institute, Philadelphia, Pennsylvania
| | | | - Michael L Klein
- Institute for Computational Molecular Science and Temple Materials Institute, Philadelphia, Pennsylvania
| | - Giacomo Fiorin
- Institute for Computational Molecular Science and Temple Materials Institute, Philadelphia, Pennsylvania.
| |
Collapse
|
47
|
Nandi S, Pyne A, Ghosh M, Banerjee P, Ghosh B, Sarkar N. Antagonist Effects of l-Phenylalanine and the Enantiomeric Mixture Containing d-Phenylalanine on Phospholipid Vesicle Membrane. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2020; 36:2459-2473. [PMID: 32073868 DOI: 10.1021/acs.langmuir.9b03543] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
One of the congenital flaws of metabolism, phenylketonuria (PKU), is known to be related to the self-assembly of toxic fibrillar aggregates of phenylalanine (Phe) in blood at elevated concentrations. Our experimental findings using l-phenylalanine (l-Phe) at millimolar concentration suggest the formation of fibrillar morphologies in the dry phase, which in the solution phase interact strongly with the model membrane composed of 1,2-diacyl-sn-glycero-phosphocholine (LAPC) lipid, thereby decreasing the rigidity (or increasing the fluidity) of the membrane. The hydrophobic interaction, in addition to the electrostatic attraction of Phe with the model membrane, is found to be responsible for such phenomena. On the contrary, various microscopic observations reveal that such fibrillar morphologies of l-Phe are severely ruptured in the presence of its enantiomer d-phenylalanine (d-Phe), thereby converting the fibrillar morphologies into crushed flakes. Various biophysical studies, including the solvation dynamics experiment, suggest that this l-Phe in the presence of d-Phe, when interacting with the same model membrane, now reverts the rigidity of the membrane, i.e., increases the rigidity of the membrane, which was lost due to interaction with l-Phe exclusively. Fluorescence anisotropy measurements also support this reverse rigid character of the membrane in the presence of an enantiomeric mixture of amino acids. A comprehensive understanding of the interaction of Phe with the model membrane is further pursued at the single-molecular fluorescence detection level using fluorescence correlation spectroscopy (FCS) experiments. Therefore, our experimental conclusion interprets a linear correlation between increased permeability and enhanced fluidity of the membrane in the presence of l-Phe and certifies d-Phe as a therapeutic modulator of l-Phe fibrillar morphologies. Further, the study proposes that the rigidity of the membrane lost due to interaction with l-Phe was reinstated-in fact, increased-in the presence of the enantiomeric mixture containing both d- and l-Phe.
Collapse
|
48
|
Loschwitz J, Olubiyi OO, Hub JS, Strodel B, Poojari CS. Computer simulations of protein-membrane systems. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2020; 170:273-403. [PMID: 32145948 PMCID: PMC7109768 DOI: 10.1016/bs.pmbts.2020.01.001] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
The interactions between proteins and membranes play critical roles in signal transduction, cell motility, and transport, and they are involved in many types of diseases. Molecular dynamics (MD) simulations have greatly contributed to our understanding of protein-membrane interactions, promoted by a dramatic development of MD-related software, increasingly accurate force fields, and available computer power. In this chapter, we present available methods for studying protein-membrane systems with MD simulations, including an overview about the various all-atom and coarse-grained force fields for lipids, and useful software for membrane simulation setup and analysis. A large set of case studies is discussed.
Collapse
Affiliation(s)
- Jennifer Loschwitz
- Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf, Düsseldorf, Germany; Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, Jülich, Germany
| | - Olujide O Olubiyi
- Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, Jülich, Germany; Department of Pharmaceutical Chemistry, Faculty of Pharmacy, Obafemi Awolowo University, Ile-Ife, Nigeria
| | - Jochen S Hub
- Theoretical Physics and Center for Biophysics, Saarland University, Saarbrücken, Germany
| | - Birgit Strodel
- Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf, Düsseldorf, Germany; Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, Jülich, Germany
| | - Chetan S Poojari
- Theoretical Physics and Center for Biophysics, Saarland University, Saarbrücken, Germany.
| |
Collapse
|
49
|
Gu RX, Baoukina S, Tieleman DP. Phase Separation in Atomistic Simulations of Model Membranes. J Am Chem Soc 2020; 142:2844-2856. [DOI: 10.1021/jacs.9b11057] [Citation(s) in RCA: 35] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Ruo-Xu Gu
- Centre for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive, N.W., Calgary, Alberta T2N 1N4, Canada
| | - Svetlana Baoukina
- Centre for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive, N.W., Calgary, Alberta T2N 1N4, Canada
| | - D. Peter Tieleman
- Centre for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive, N.W., Calgary, Alberta T2N 1N4, Canada
| |
Collapse
|
50
|
Susceptibility of biomembrane structure towards amphiphiles, ionic liquids, and deep eutectic solvents. ADVANCES IN BIOMEMBRANES AND LIPID SELF-ASSEMBLY 2020. [DOI: 10.1016/bs.abl.2020.02.001] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
|