1
|
Chauhan G, Bremer A, Dar F, Mittag T, Pappu RV. Crowder titrations enable the quantification of driving forces for macromolecular phase separation. Biophys J 2024; 123:1376-1392. [PMID: 37717144 PMCID: PMC11163301 DOI: 10.1016/j.bpj.2023.09.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2023] [Revised: 09/03/2023] [Accepted: 09/13/2023] [Indexed: 09/18/2023] Open
Abstract
Macromolecular solubility is an important contributor to the driving forces for phase separation. Formally, the driving forces in a binary mixture comprising a macromolecule dissolved in a solvent can be quantified in terms of the saturation concentration, which is the threshold macromolecular concentration above which the mixture separates into coexisting dense and dilute phases. In addition, the second virial coefficient, which measures the effective strength of solvent-mediated intermolecular interactions provides direct assessments of solvent quality. The sign and magnitude of second virial coefficients will be governed by a combination of solution conditions and the nature of the macromolecule of interest. Here, we show, using a combination of theory, simulation, and in vitro experiments, that titrations of crowders, providing they are true depletants, can be used to extract the intrinsic driving forces for macromolecular phase separation. This refers to saturation concentrations in the absence of crowders and the second virial coefficients that quantify the magnitude of the incompatibility between macromolecules and the solvent. Our results show how the depletion-mediated attractions afforded by crowders can be leveraged to obtain comparative assessments of macromolecule-specific, intrinsic driving forces for phase separation.
Collapse
Affiliation(s)
- Gaurav Chauhan
- Department of Biomedical Engineering and Center for Biomolecular Condensates, James F. McKelvey School of Engineering, Washington University in St. Louis, St. Louis, Missouri
| | - Anne Bremer
- Department of Structural Biology, St. Jude Children's Research Hospital, Memphis, Tennessee
| | - Furqan Dar
- Department of Biomedical Engineering and Center for Biomolecular Condensates, James F. McKelvey School of Engineering, Washington University in St. Louis, St. Louis, Missouri
| | - Tanja Mittag
- Department of Structural Biology, St. Jude Children's Research Hospital, Memphis, Tennessee
| | - Rohit V Pappu
- Department of Biomedical Engineering and Center for Biomolecular Condensates, James F. McKelvey School of Engineering, Washington University in St. Louis, St. Louis, Missouri.
| |
Collapse
|
2
|
Mondal S, Cui Q. Sequence Sensitivity in Membrane Remodeling by Polyampholyte Condensates. J Phys Chem B 2024; 128:2087-2099. [PMID: 38407041 DOI: 10.1021/acs.jpcb.3c08149] [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: 02/27/2024]
Abstract
Intrinsically disordered peptides (IDPs) have been found to undergo liquid-liquid phase separation (LLPS) and produce complex coacervates that play numerous regulatory roles in the cell. Recent experimental studies have discovered that LLPS at or near the membrane surface helps in the biomolecular organization during signaling events and can significantly alter the membrane morphology. However, the molecular mechanism and microscopic details of such processes still remain unclear. Here we study the effect of polyampholyte and polyelectrolyte condensation on two different anionic membranes, as they represent a majority of naturally occurring IDPs. The polyampholytes are fifty-residue polymers, made of glutamate(E) and lysine(K) with different charge patterns. The polyelectrolytes are separate chains of E25 and K25. We first calibrate the MARTINI v3.0 force field and then perform long-time-scale coarse-grained molecular dynamics simulations. We find that condensates formed by all the polyampholytes get adsorbed on the membrane. However, the strong polyampholytes (i.e., blocky sequences) can remodel the membranes more prominently than the weaker ones (i.e., scrambled sequences). Condensates formed by the blocky sequences induce a significant negative curvature (∼0.1 nm-1) and local demixing of lipids, whereas those by the scrambled sequences tend to wet the membrane to a greater extent without generating significant curvature or demixing. We perform several microscopic analyses to characterize the nature of the interaction between membranes and these condensates. Our analyses of interaction energetics reveal that membrane remodeling and/or wetting are favored by enhanced interactions between polyampholytes with lipids and the counterions.
Collapse
Affiliation(s)
- Sayantan Mondal
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Qiang Cui
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
- Department of Physics, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
- Department of Biomedical Engineering, Boston University, 44 Cummington Mall, Boston, Massachusetts 02215, United States
| |
Collapse
|
3
|
Macias E, Travesset A. Hydrogen Bond Network Disruption by Hydration Layers in Water Solutions with Salt and Hydrogen-Bonding Polymers (PEO). J Phys Chem B 2023. [PMID: 37478338 DOI: 10.1021/acs.jpcb.3c02505] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/23/2023]
Abstract
A mean field theory model describing the interaction of ion hydration layers with the network of hydrogen bonds of both water and the nonionic polymer poly(ethylene oxide) (PEO) is presented. The predictions of the model for types and statistics of hydrogen bonds, the number of water molecules bound to PEO, or their dependence on temperature are successfully verified from all-atom simulations at different NaCl and PEO concentrations. Furthermore, our simulations show that the binding of cations to PEO increases monotonically with salt concentration, in agreement with recent experimental results, through a mechanism in which the sum of the number of bound water and cations is independent of salt concentration. The model introduced is general and can describe any salt or hydrogen-bond-forming polymer.
Collapse
Affiliation(s)
- Elizabeth Macias
- Department of Physics and Astronomy, Iowa State University and Ames Lab, Ames, Iowa 50011, United States
| | - Alex Travesset
- Department of Physics and Astronomy, Iowa State University and Ames Lab, Ames, Iowa 50011, United States
| |
Collapse
|
4
|
Chauhan G, Bremer A, Dar F, Mittag T, Pappu RV. Crowder titrations enable the quantification of driving forces for macromolecular phase separation. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.07.03.547544. [PMID: 37461587 PMCID: PMC10350001 DOI: 10.1101/2023.07.03.547544] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/22/2023]
Abstract
Macromolecular solubility is an important contributor to the driving forces for phase separation. Formally, the driving forces in a binary mixture comprising a macromolecule dissolved in a solvent can be quantified in terms of the saturation concentration, which is the threshold macromolecular concentration above which the mixture separates into coexisting dense and dilute phases. Additionally, the second virial coefficient, which measures the effective strength of solvent-mediated intermolecular interactions provides direct assessments of solvent quality. The sign and magnitude of second virial coefficients will be governed by a combination of solution conditions and the nature of the macromolecule of interest. Here, we show, using a combination of theory, simulation, and in vitro experiments, that titrations of crowders, providing they are true depletants, can be used to extract the intrinsic driving forces for macromolecular phase separation. This refers to saturation concentrations in the absence of crowders and the second virial coefficients that quantify the magnitude of the incompatibility between macromolecules and the solvent. Our results show how the depletion-mediated attractions afforded by crowders can be leveraged to obtain comparative assessments of macromolecule-specific, intrinsic driving forces for phase separation. SIGNIFICANCE Phase separation has emerged as a process of significant relevance to sorting macromolecules into distinct compartments, thereby enabling spatial and temporal control over cellular matter. Considerable effort is being invested into uncovering the driving forces that enable the separation of macromolecular solutions into coexisting phases. At its heart, this process is governed by the balance of macromolecule-solvent, inter-macromolecule, and solvent-solvent interactions. We show that the driving forces for phase separation, including the coefficients that measure interaction strengths between macromolecules, can be extracted by titrating the concentrations of crowders that enable macromolecules to phase separate at lower concentrations. Our work paves the way to leverage specific categories of measurements for quantitative characterizations of driving forces for phase separation.
Collapse
|
5
|
Alessandri R, Grünewald F, Marrink SJ. The Martini Model in Materials Science. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2021; 33:e2008635. [PMID: 33956373 DOI: 10.1002/adma.202008635] [Citation(s) in RCA: 48] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/22/2020] [Revised: 03/15/2021] [Indexed: 06/12/2023]
Abstract
The Martini model, a coarse-grained force field initially developed with biomolecular simulations in mind, has found an increasing number of applications in the field of soft materials science. The model's underlying building block principle does not pose restrictions on its application beyond biomolecular systems. Here, the main applications to date of the Martini model in materials science are highlighted, and a perspective for the future developments in this field is given, particularly in light of recent developments such as the new version of the model, Martini 3.
Collapse
Affiliation(s)
- Riccardo Alessandri
- Zernike Institute for Advanced Materials and Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Nijenborgh 4, Groningen, 9747AG, The Netherlands
| | - Fabian Grünewald
- Zernike Institute for Advanced Materials and Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Nijenborgh 4, Groningen, 9747AG, The Netherlands
| | - Siewert J Marrink
- Zernike Institute for Advanced Materials and Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Nijenborgh 4, Groningen, 9747AG, The Netherlands
| |
Collapse
|
6
|
Effectiveness of coarse graining degree and speedup on the dynamic properties of homopolymer. J Mol Model 2021; 27:55. [PMID: 33511476 DOI: 10.1007/s00894-020-04661-5] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2020] [Accepted: 12/20/2020] [Indexed: 10/22/2022]
Abstract
Evaluation of effective coarse graining (CG) degree and reasonable speedup relative to all-atomistic (AA) model was conducted to provide a basis for building appropriate larger-scale model. The reproducibility of atomistic conformation and temperature transferability both act as the analysis criteria to resolve the maximum acceptable CG degree. Taking short- and long time spans into account simultaneously in the estimation of computational speedup, a dynamic scaling factor is accessible in fitting mean squared displacement ratio of CG to AA as an exponential function. Computing loss in parallel running is an indispensable component in acceleration, which was also added in the evaluation. Subsequently, a quantified prediction of CG speedup arises as a multiplication of dynamic scaling factor, computing loss, time step, and the square of reduction in the number of degrees of freedom. Polyethylene oxide was adopted as a reference system to execute the direct Boltzmann inversion and iterative Boltzmann inversion. Bonded and non-bonded potentials were calculated in CG models with 1~4 monomers per bead. The effective CG degree was determined as two at the most with a speedup of four orders magnitude over AA in this study. Determination of effectiveness CG degree and the corresponding speedup prediction provide available tools in larger spatiotemporal-scale calculations.
Collapse
|
7
|
Sadeghi MS, Moghbeli MR, Goddard WA. A coarse-grain force field based on quantum mechanics (CGq FF) for molecular dynamics simulation of poly(ethylene glycol)-block-poly(ε-caprolactone) (PEG-b-PCL) micelles. Phys Chem Chem Phys 2020; 22:24028-24040. [PMID: 33078174 DOI: 10.1039/d0cp04364h] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
In order to provide the means to predict from molecular dynamics (MD) simulations the structures of copolymer-based micelles in solution, we developed coarse grain force field (CGq FF) parameters for poly(ethylene glycol) (PEG) and for poly(ε-caprolactone) (PCL). A key advance here is the use of quantum mechanics to train the parameters describing the non-bonded (NB) interactions between the CG beads. The functional forms are the same as the MARTINI CG FF so standard MD codes can be used. Our CGq FF describes well the experimentally observed properties for the polymer-air and polymer-water interfaces, indicating the accuracy of the NB interactions. The structural properties (density, radius of gyration (Rg), and end-to-end distance (h)) match both experiment and all atom (AA) simulations. We illustrate the application of this CGq FF by following the formation of a spherical micelle from 250 chains of PEG23-b-PCL9 diblock copolymer, each block with molecular weight of 1000 Daltons (10 500 beads, corresponding to 123 250 atoms), in a water box with 119 139 water beads (426 553 water molecules).
Collapse
Affiliation(s)
- Maryam S Sadeghi
- Smart Polymers and Nanocomposites Research Group, School of Chemical Engineering, Iran University of Science and Technology, Tehran 16846-13114, Iran.
| | - Mohammad Reza Moghbeli
- Smart Polymers and Nanocomposites Research Group, School of Chemical Engineering, Iran University of Science and Technology, Tehran 16846-13114, Iran.
| | - William A Goddard
- Materials and Process Simulation Center (MSC), California Institute of Technology, Pasadena, California 91125, USA.
| |
Collapse
|
8
|
Ma SM, Zhao L, Wang YL, Zhu YL, Lu ZY. The coarse-grained models of poly(ethylene oxide) and poly(propylene oxide) homopolymers and poloxamers in big multipole water (BMW) and MARTINI frameworks. Phys Chem Chem Phys 2020; 22:15976-15985. [PMID: 32632434 DOI: 10.1039/d0cp01006e] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Polyethylene oxide (PEO) and poly(propylene oxide) (PPO), especially their tri-block copolymers PEO-PPO-PEO (poloxamers), have a broad range of applications in biotechnology and medical science. Understanding their specific interactions with biomembranes is the key to unveil the unique features of poloxamers either as membrane-healing or membrane pore-forming agents. Based on the coarse-graining convention of the MARTINI force field and the big multipole water (BMW) model, which has a three charged site topology and can reproduce the correct dipole moment of four-water clusters, we generated coarse-grained (CG) models with analytical and numerical potentials for PEO and PPO homopolymers and poloxamers in dilute solution. The effective bonded interaction potentials between CG beads were determined from the probability distributions of bond lengths, angles and dihedrals that are determined from atomistic simulations. The nonbonded interaction parameters were fine-tuned to reproduce the conformational properties of atomistic PEO and PPO homopolymers and poloxamers via extensive CG simulations of PEO and PPO homopolymers and poloxamers in a BMW water environment. The reported CG models provide a promising framework for a comprehensive understanding of the microstructural, conformational, and dynamic properties of poloxamers and their delicate interactions with other species in an explicit water environment.
Collapse
Affiliation(s)
- Su-Min Ma
- State Key Laboratory of Supramolecular Structure and Materials, Institute of Theoretical Chemistry, Jilin University, Changchun 130023, China.
| | - Li Zhao
- College of Life Sciences, Jilin University, Changchun 130012, China
| | - Yong-Lei Wang
- Department of Materials and Environmental Chemistry, Arrhenius Laboratory, Stockholm University, SE-10691, Stockholm, Sweden.
| | - You-Liang Zhu
- State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, China
| | - Zhong-Yuan Lu
- State Key Laboratory of Supramolecular Structure and Materials, Institute of Theoretical Chemistry, Jilin University, Changchun 130023, China.
| |
Collapse
|
9
|
Lee H. Molecular Simulations of PEGylated Biomolecules, Liposomes, and Nanoparticles for Drug Delivery Applications. Pharmaceutics 2020; 12:E533. [PMID: 32531886 PMCID: PMC7355693 DOI: 10.3390/pharmaceutics12060533] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2020] [Revised: 06/04/2020] [Accepted: 06/08/2020] [Indexed: 12/17/2022] Open
Abstract
Since the first polyethylene glycol (PEG)ylated protein was approved by the FDA in 1990, PEGylation has been successfully applied to develop drug delivery systems through experiments, but these experimental results are not always easy to interpret at the atomic level because of the limited resolution of experimental techniques. To determine the optimal size, structure, and density of PEG for drug delivery, the structure and dynamics of PEGylated drug carriers need to be understood close to the atomic scale, as can be done using molecular dynamics simulations, assuming that these simulations can be validated by successful comparisons to experiments. Starting with the development of all-atom and coarse-grained PEG models in 1990s, PEGylated drug carriers have been widely simulated. In particular, recent advances in computer performance and simulation methodologies have allowed for molecular simulations of large complexes of PEGylated drug carriers interacting with other molecules such as anticancer drugs, plasma proteins, membranes, and receptors, which makes it possible to interpret experimental observations at a nearly atomistic resolution, as well as help in the rational design of drug delivery systems for applications in nanomedicine. Here, simulation studies on the following PEGylated drug topics will be reviewed: proteins and peptides, liposomes, and nanoparticles such as dendrimers and carbon nanotubes.
Collapse
Affiliation(s)
- Hwankyu Lee
- Department of Chemical Engineering, Dankook University, Yongin 16890, Korea
| |
Collapse
|
10
|
Das M, Dahal U, Mesele O, Liang D, Cui Q. Molecular Dynamics Simulation of Interaction between Functionalized Nanoparticles with Lipid Membranes: Analysis of Coarse-Grained Models. J Phys Chem B 2019; 123:10547-10561. [DOI: 10.1021/acs.jpcb.9b08259] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Affiliation(s)
- Mitradip Das
- School of Chemical Sciences, National Institute of Science Education and Research, Khordha, Odisha, India, 752050
- Homi Bhabha National Institute, Training School
Complex, Anushaktinagar, Mumbai, Maharashtra, India, 400094
| | - Udaya Dahal
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Oluwaseun Mesele
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Dongyue Liang
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin−Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Qiang Cui
- Departments of Chemistry, Physics and Biomedical Engineering, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| |
Collapse
|
11
|
McDonnell MT, Keffer DJ. Reactive molecular dynamics simulations of an excess proton in polyethylene glycol-water solutions. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1557328] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Affiliation(s)
| | - David J. Keffer
- Department of Materials Science and Engineering, University of Tennessee, Knoxville, TN, USA
| |
Collapse
|
12
|
Sousa SF, Peres J, Coelho M, Vieira TF. Analyzing PEGylation through Molecular Dynamics Simulations. ChemistrySelect 2018. [DOI: 10.1002/slct.201800855] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Sérgio F. Sousa
- UCIBIO@REQUIMTE; BioSIM; Departamento de Biomedicina; Faculdade de Medicina da Universidade do Porto, Alameda Professor Hernâni Monteiro; 4200-319, Porto Portugal
| | - Joana Peres
- LEPABE; Faculdade de Engenharia; Universidade do Porto, Porto; Portugal
| | - Manuel Coelho
- LEPABE; Faculdade de Engenharia; Universidade do Porto, Porto; Portugal
| | - Tatiana F. Vieira
- LEPABE; Faculdade de Engenharia; Universidade do Porto, Porto; Portugal
| |
Collapse
|
13
|
Son CY, McDaniel JG, Cui Q, Yethiraj A. Conformational and Dynamic Properties of Poly(ethylene oxide) in BMIM+BF4–: A Microsecond Computer Simulation Study Using ab Initio Force Fields. Macromolecules 2018. [DOI: 10.1021/acs.macromol.8b01002] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Chang Yun Son
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin—Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Jesse G. McDaniel
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin—Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Qiang Cui
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin—Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Arun Yethiraj
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin—Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| |
Collapse
|
14
|
De Silva CC, Leophairatana P, Ohkuma T, Koberstein JT, Kremer K, Mukherji D. Sequence transferable coarse-grained model of amphiphilic copolymers. J Chem Phys 2018; 147:064904. [PMID: 28810779 DOI: 10.1063/1.4997638] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
Polymer properties are inherently multi-scale in nature, where delicate local interaction details play a key role in describing their global conformational behavior. In this context, deriving coarse-grained (CG) multi-scale models for polymeric liquids is a non-trivial task. Further complexities arise when dealing with copolymer systems with varying microscopic sequences, especially when they are of amphiphilic nature. In this work, we derive a segment-based generic CG model for amphiphilic copolymers consisting of repeat units of hydrophobic (methylene) and hydrophilic (ethylene oxide) monomers. The system is a simulation analogue of polyacetal copolymers [S. Samanta et al., Macromolecules 49, 1858 (2016)]. The CG model is found to be transferable over a wide range of copolymer sequences and also to be consistent with existing experimental data.
Collapse
Affiliation(s)
| | | | - Takahiro Ohkuma
- Central Research Division, Bridgestone Corporation, Kodaira, Tokyo 187-8531, Japan
| | - Jeffrey T Koberstein
- Department of Chemical Engineering, Columbia University, Columbia, New York 10027, USA
| | - Kurt Kremer
- Max-Planck Institut für Polymerforschung, Ackermannweg 10, 55128 Mainz, Germany
| | - Debashish Mukherji
- Max-Planck Institut für Polymerforschung, Ackermannweg 10, 55128 Mainz, Germany
| |
Collapse
|
15
|
Barden DR, Vashisth H. Parameterization and atomistic simulations of biomimetic membranes. Faraday Discuss 2018; 209:161-178. [DOI: 10.1039/c8fd00047f] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We have developed CHARMM force-field compatible parameters and conducted all-atom explicit-solvent MD simulations of biomimetic membranes composed of block copolymers of poly(butadiene), poly(isoprene), and poly(ethylene oxide).
Collapse
Affiliation(s)
| | - Harish Vashisth
- Department of Chemical Engineering
- University of New Hampshire
- Durham
- USA
| |
Collapse
|
16
|
Wu C. Multiscale Modeling Scheme for Simulating Polymeric Melts: Application to Poly(Ethylene Oxide). MACROMOL THEOR SIMUL 2017. [DOI: 10.1002/mats.201700066] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- C. Wu
- Hunan Provincial Key Laboratory of Fine Ceramics and Powder Materials; Hunan University of Humanities Science & Technology; Dixing Road 487, Louxing District Loudi 417000 Hunan Province P. R. China
- College of Materials and Environment Engineering; Hunan University of Humanities Science & Technology; Dixing Road 487, Louxing District Loudi 417000 Hunan Province P. R. China
| |
Collapse
|
17
|
Chudoba R, Heyda J, Dzubiella J. Temperature-Dependent Implicit-Solvent Model of Polyethylene Glycol in Aqueous Solution. J Chem Theory Comput 2017; 13:6317-6327. [DOI: 10.1021/acs.jctc.7b00560] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Richard Chudoba
- Institut
für Physik, Humboldt-Universität zu Berlin, Newtonstraße
15, D-12489 Berlin, Germany
- Institut
für Weiche Materie und Funktionale Materialen, Helmholtz-Zentrum Berlin, Hahn-Meitner-Platz 1, D-14109 Berlin, Germany
| | - Jan Heyda
- Department
of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, CZ-16628 Praha 6, Czech Republic
| | - Joachim Dzubiella
- Institut
für Physik, Humboldt-Universität zu Berlin, Newtonstraße
15, D-12489 Berlin, Germany
- Institut
für Weiche Materie und Funktionale Materialen, Helmholtz-Zentrum Berlin, Hahn-Meitner-Platz 1, D-14109 Berlin, Germany
| |
Collapse
|
18
|
Taddese T, Carbone P. Effect of Chain Length on the Partition Properties of Poly(ethylene oxide): Comparison between MARTINI Coarse-Grained and Atomistic Models. J Phys Chem B 2017; 121:1601-1609. [PMID: 28151665 DOI: 10.1021/acs.jpcb.6b10858] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The MARTINI coarse-grained beads are parameterized to match the partition coefficients of several organic molecules in different solvents. Here, we test the method when modeling the partitioning properties of poly(ethylene oxide) between solvents of different polarities. We show that, among the existing models, the latest model developed by Lee and co-workers [ Lee , H. ; Pastor , R. W. J. Phys. Chem. B 2011 , 115 , 7830 - 7837 ] is the one that most successfully reproduces the hydration free energy of short oligomers, although it predicts highly negative solvation free energies in octanol and hexane. We develop a new CG model matching the solvation free energy of the monomer in different solvents and propose a simple method to select the Lennard-Jones parameters that reproduce the desired partition coefficients. The model correctly reproduces water/hexane partition properties for oligomers up to 10 monomers but still suffers from a transferability problem for larger molecular weight.
Collapse
Affiliation(s)
- Tseden Taddese
- School of Chemical Engineering and Analytical Science, The University of Manchester , Oxford Road, Manchester M13 9PL, U.K
| | - Paola Carbone
- School of Chemical Engineering and Analytical Science, The University of Manchester , Oxford Road, Manchester M13 9PL, U.K
| |
Collapse
|
19
|
Lu J, Jacobson LC, Perez Sirkin YA, Molinero V. High-Resolution Coarse-Grained Model of Hydrated Anion-Exchange Membranes that Accounts for Hydrophobic and Ionic Interactions through Short-Ranged Potentials. J Chem Theory Comput 2016; 13:245-264. [PMID: 28068769 DOI: 10.1021/acs.jctc.6b00874] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Jibao Lu
- Department
of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, United States
| | - Liam C. Jacobson
- Department
of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, United States
| | - Yamila A. Perez Sirkin
- Departamento
de Química Inorgánica, Analítica y Química
Física, and INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires C1428EHA, Argentina
| | - Valeria Molinero
- Department
of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, United States
| |
Collapse
|
20
|
Oroskar PA, Jameson CJ, Murad S. Rotational behaviour of PEGylated gold nanorods in a lipid bilayer system. Mol Phys 2016. [DOI: 10.1080/00268976.2016.1248515] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Priyanka A. Oroskar
- Department of Chemical Engineering, University of Illinois at Chicago, Chicago, IL, USA
| | - Cynthia J. Jameson
- Department of Chemical Engineering, University of Illinois at Chicago, Chicago, IL, USA
- Department of Chemistry, University of Illinois at Chicago, Chicago, IL, USA
| | - Sohail Murad
- Department of Chemical Engineering, University of Illinois at Chicago, Chicago, IL, USA
- Department of Chemical and Biological Engineering, Illinois Institute of Technology, Chicago, IL, USA
| |
Collapse
|
21
|
Bunker A, Magarkar A, Viitala T. Rational design of liposomal drug delivery systems, a review: Combined experimental and computational studies of lipid membranes, liposomes and their PEGylation. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2016; 1858:2334-2352. [DOI: 10.1016/j.bbamem.2016.02.025] [Citation(s) in RCA: 103] [Impact Index Per Article: 12.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/21/2015] [Revised: 02/09/2016] [Accepted: 02/10/2016] [Indexed: 01/22/2023]
|
22
|
Oroskar PA, Jameson CJ, Murad S. Simulated Permeation and Characterization of PEGylated Gold Nanoparticles in a Lipid Bilayer System. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2016; 32:7541-7555. [PMID: 27399834 DOI: 10.1021/acs.langmuir.6b01740] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
PEGylated gold nanoparticles are considered suitable nanocarriers for use in biomedical applications and targeted drug delivery systems. In our previous investigation with the alkanethiol-functionalized gold nanoparticle, we found that permeation across a protein-free phospholipid membrane resulted in damaging effects of lipid displacement and water and ion leakage. In the present study, we carry out a series of coarse-grained molecular simulations to explore permeation of lipid bilayer systems by a PEGylated gold nanoparticle, especially at the bulk-liquid-lipid interface as well as the interface between the two lipid leaflets. Initially, we examine molecular-level details of a PEGylated gold nanoparticle (constructed from cycled annealing) in water and find a distribution of ligand configurations (from mushroom to brush states) present in nanoparticles with medium to high surface coverage. We also find that the characteristic properties of the PEGylated gold nanoparticle do not change when it is placed in a salt solution. In our permeation studies, we investigate events of water and ion penetration as well as lipid translocation while varying the ligand length, nanoparticle surface coverage, and ion concentration gradient of our system. Results from our studies show the following: (1) The number of water molecules in the interior of the membrane during ligand-coated nanoparticle permeation increases with PEGn-SH surface coverage, ligand length, and permeation velocity but is not sensitive to the ion concentration gradient. (2) Lipid molecules do not leave the membrane; instead they complete trans-bilayer lipid flip-flop with longer ligands and higher surface coverages. (3) The lack of formation of stable water pores prevents ion translocation. (4) The PEGylated nanoparticle causes less damage to the membrane overall due to favorable interactions with the lipid headgroups which may explain why experimentalists observe endocytosis of PEGylated nanocarriers in vivo.
Collapse
Affiliation(s)
- Priyanka A Oroskar
- Department of Chemical Engineering, University of Illinois at Chicago , 810 South Clinton Street, Chicago, Illinois 60607, United States
| | - Cynthia J Jameson
- Department of Chemical Engineering, University of Illinois at Chicago , 810 South Clinton Street, Chicago, Illinois 60607, United States
- Department of Chemistry, University of Illinois at Chicago , 845 West Taylor Street, Chicago, Illinois 60607, United States
| | - Sohail Murad
- Department of Chemical Engineering, University of Illinois at Chicago , 810 South Clinton Street, Chicago, Illinois 60607, United States
- Department of Chemical and Biological Engineering, Illinois Institute of Technology , 3300 South Federal Street, Chicago, Illinois 60616, United States
| |
Collapse
|
23
|
Li X, Murthy NS, Becker ML, Latour RA. Multiscale approach for the construction of equilibrated all-atom models of a poly(ethylene glycol)-based hydrogel. Biointerphases 2016; 11:021002. [PMID: 27013229 PMCID: PMC4808066 DOI: 10.1116/1.4944829] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2016] [Revised: 03/09/2016] [Accepted: 03/15/2016] [Indexed: 11/17/2022] Open
Abstract
A multiscale modeling approach is presented for the efficient construction of an equilibrated all-atom model of a cross-linked poly(ethylene glycol) (PEG)-based hydrogel using the all-atom polymer consistent force field (PCFF). The final equilibrated all-atom model was built with a systematic simulation toolset consisting of three consecutive parts: (1) building a global cross-linked PEG-chain network at experimentally determined cross-link density using an on-lattice Monte Carlo method based on the bond fluctuation model, (2) recovering the local molecular structure of the network by transitioning from the lattice model to an off-lattice coarse-grained (CG) model parameterized from PCFF, followed by equilibration using high performance molecular dynamics methods, and (3) recovering the atomistic structure of the network by reverse mapping from the equilibrated CG structure, hydrating the structure with explicitly represented water, followed by final equilibration using PCFF parameterization. The developed three-stage modeling approach has application to a wide range of other complex macromolecular hydrogel systems, including the integration of peptide, protein, and/or drug molecules as side-chains within the hydrogel network for the incorporation of bioactivity for tissue engineering, regenerative medicine, and drug delivery applications.
Collapse
Affiliation(s)
- Xianfeng Li
- Department of Bioengineering, Clemson University, Clemson, South Carolina 29634
| | - N Sanjeeva Murthy
- New Jersey Center for Biomaterials, Rutgers University, Piscataway, New Jersey 08854
| | - Matthew L Becker
- Departments of Polymer Science and Biomedical Engineering, The University of Akron, Akron, Ohio 44325
| | - Robert A Latour
- Department of Bioengineering, Clemson University, Clemson, South Carolina 29634
| |
Collapse
|
24
|
Xie F, Turesson M, Jansson M, Skepö M, Forsman J. A simple and versatile implicit solvent model for polyethylene glycol in aqueous solution at room temperature. POLYMER 2016. [DOI: 10.1016/j.polymer.2015.12.034] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
25
|
Mantha S, Yethiraj A. Conformational Properties of Sodium Polystyrenesulfonate in Water: Insights from a Coarse-Grained Model with Explicit Solvent. J Phys Chem B 2015; 119:11010-8. [DOI: 10.1021/acs.jpcb.5b01700] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Sriteja Mantha
- Theoretical
Chemistry Institute, Department of Chemistry, University of Wisconsin, 1101 University
Avenue, Madison, Wisconsin 53706, United States
| | - Arun Yethiraj
- Theoretical
Chemistry Institute, Department of Chemistry, University of Wisconsin, 1101 University
Avenue, Madison, Wisconsin 53706, United States
| |
Collapse
|
26
|
Choi E, Yethiraj A. Conformational Properties of a Polymer in an Ionic Liquid: Computer Simulations and Integral Equation Theory of a Coarse-Grained Model. J Phys Chem B 2014; 119:9091-7. [DOI: 10.1021/jp508876q] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Eunsong Choi
- Department
of Physics, University of Wisconsin, Madison, Wisconsin 53706, United States
| | - Arun Yethiraj
- Department
of Chemistry, University of Wisconsin, Madison, Wisconsin 53706, United States
| |
Collapse
|
27
|
Chaudhari MI, Pratt LR, Paulaitis ME. Loop-Closure and Gaussian Models of Collective Structural Characteristics of Capped PEO Oligomers in Water. J Phys Chem B 2014; 119:8863-7. [DOI: 10.1021/jp504244x] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- M. I. Chaudhari
- Department of Chemical and
Biomolecular Engineering, Tulane University, New Orleans, Louisiana 70118, United States
| | - L. R. Pratt
- Department of Chemical and
Biomolecular Engineering, Tulane University, New Orleans, Louisiana 70118, United States
| | - M. E. Paulaitis
- Department
of Chemical and
Biomolecular Engineering, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|