1
|
Ledum M, Sen S, Bore SL, Cascella M. On the equivalence of the hybrid particle-field and Gaussian core models. J Chem Phys 2023; 158:2890484. [PMID: 37184022 DOI: 10.1063/5.0145142] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2023] [Accepted: 04/28/2023] [Indexed: 05/16/2023] Open
Abstract
Hybrid particle-field molecular dynamics is a molecular simulation strategy, wherein particles couple to a density field instead of through ordinary pair potentials. Traditionally considered a mean-field theory, a momentum and energy-conserving hybrid particle-field formalism has recently been introduced, which was demonstrated to approach the Gaussian Core model potential in the grid-converged limit. Here, we expand on and generalize the correspondence between the Hamiltonian hybrid particle-field method and particle-particle pair potentials. Using the spectral procedure suggested by Bore and Cascella, we establish compatibility to any local soft pair potential in the limit of infinitesimal grid spacing. Furthermore, we document how the mean-field regime often observed in hybrid particle-field simulations is due to the systems under consideration, and not an inherent property of the model. Considering the Gaussian filter form, in particular, we demonstrate the ability of the Hamiltonian hybrid particle-field model to recover all structural and dynamical properties of the Gaussian Core model, including solid phases, a first-order phase transition, and anomalous transport properties. We quantify the impact of the grid spacing on the correspondence, as well as the effect of the particle-field filtering length scale on the emergent particle-particle correlations.
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
| | - Sigbjørn Løland Bore
- Department of Chemistry and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, PO Box 1033 Blindern, 0315 Oslo, Norway
| | - Michele Cascella
- Department of Chemistry and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, PO Box 1033 Blindern, 0315 Oslo, Norway
| |
Collapse
|
2
|
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
|
3
|
Sun L, Wang S, Tian F, Zhu H, Dai L. Organizations of melittin peptides after spontaneous penetration into cell membranes. Biophys J 2022; 121:4368-4381. [PMID: 36199252 PMCID: PMC9703044 DOI: 10.1016/j.bpj.2022.10.002] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2022] [Revised: 08/16/2022] [Accepted: 09/30/2022] [Indexed: 12/14/2022] Open
Abstract
The antimicrobial peptide, melittin, is a potential next-generation antibiotic because melittin can spontaneously form pores in bacterial cell membranes and cause cytoplasm leakage. However, the organizations of melittin peptides in cell membranes remain elusive, which impedes the understanding of the poration mechanism. In this work, we use coarse-grained and all-atom molecular dynamics (MD) simulations to investigate the organizations of melittin peptides during and after spontaneous penetration into DPPC/POPG lipid bilayers. We find that the peptides in lipid bilayers adopt either a transmembrane conformation or a U-shaped conformation, which are referred to as T- and U-peptides, respectively. Several U-peptides and/or T-peptides aggregate to form stable pores. We analyze a T-pore consisting of four T-peptides and a U-pore consisting of three U-peptides and one T-peptide. In both pores, peptides are organized in a manner such that polar residues face inward and hydrophobic residues face outward, which stabilizes the pores and produces water channels. Compared with the U-pore, the T-pore has lower energy, larger pore diameter, and higher permeability. However, the T-pore occurs less frequently than the U-pore in our simulations, probably because the formation of the T-pore is kinetically slower than the U-pore. The stability and permeability of both pores are confirmed by 300 ns all-atom MD simulations. The peptide organizations obtained in this work should deepen the understanding of the stability, poration mechanism, and permeability of melittin, and facilitate the optimization of melittin to enhance the antibacterial ability.
Collapse
Affiliation(s)
- Liang Sun
- Department of Physics, City University of Hong Kong, Kowloon, Hong Kong, China
| | - Simin Wang
- Department of Physics, City University of Hong Kong, Kowloon, Hong Kong, China
| | - Fujia Tian
- Department of Physics, City University of Hong Kong, Kowloon, Hong Kong, China
| | - Haoqi Zhu
- Department of Physics, City University of Hong Kong, Kowloon, Hong Kong, China
| | - Liang Dai
- Department of Physics, City University of Hong Kong, Kowloon, Hong Kong, China.
| |
Collapse
|
4
|
Wu Z, Kalogirou A, De Nicola A, Milano G, Müller‐Plathe F. Atomistic hybrid
particle‐field
molecular dynamics combined with
slip‐springs
: Restoring entangled dynamics to simulations of polymer melts. J Comput Chem 2020; 42:6-18. [DOI: 10.1002/jcc.26428] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2020] [Revised: 09/07/2020] [Accepted: 09/11/2020] [Indexed: 02/03/2023]
Affiliation(s)
- Zhenghao Wu
- Eduard‐Zintl‐Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt Darmstadt Germany
| | - Andreas Kalogirou
- Eduard‐Zintl‐Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt Darmstadt Germany
| | - Antonio De Nicola
- Department of Organic Materials Science Yamagata University Yamagata‐ken Japan
| | - Giuseppe Milano
- Department of Organic Materials Science Yamagata University Yamagata‐ken Japan
| | - Florian Müller‐Plathe
- Eduard‐Zintl‐Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt Darmstadt Germany
| |
Collapse
|
5
|
Bore SL, Cascella M. Hamiltonian and alias-free hybrid particle–field molecular dynamics. J Chem Phys 2020; 153:094106. [DOI: 10.1063/5.0020733] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sigbjørn Løland Bore
- Department of Chemistry, and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, P.O. Box 1033, Blindern 0315, Oslo, Norway
| | - Michele Cascella
- Department of Chemistry, and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, P.O. Box 1033, Blindern 0315, Oslo, Norway
| |
Collapse
|
6
|
Carrer M, Škrbić T, Bore SL, Milano G, Cascella M, Giacometti A. Can Polarity-Inverted Surfactants Self-Assemble in Nonpolar Solvents? J Phys Chem B 2020; 124:6448-6458. [PMID: 32618191 PMCID: PMC8009519 DOI: 10.1021/acs.jpcb.0c04842] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
![]()
We investigate the
self-assembly process of a surfactant with inverted
polarity in water and cyclohexane using both all-atom and coarse-grained
hybrid particle-field molecular dynamics simulations. Unlike conventional
surfactants, the molecule under study, proposed in a recent experiment,
is formed by a rigid and compact hydrophobic adamantane moiety, and
a long and floppy triethylene glycol tail. In water, we report the
formation of stable inverted micelles with the adamantane heads grouping
together into a hydrophobic core and the tails forming hydrogen bonds
with water. By contrast, microsecond simulations do not provide evidence
of stable micelle formation in cyclohexane. Validating the computational
results by comparison with experimental diffusion constant and small-angle
X-ray scattering intensity, we show that at laboratory thermodynamic
conditions the mixture resides in the supercritical region of the
phase diagram, where aggregated and free surfactant states coexist
in solution. Our simulations also provide indications as to how to
escape this region to produce thermodynamically stable micellar aggregates.
Collapse
Affiliation(s)
- Manuel Carrer
- Department of Chemistry and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, P.O. Box 1033, Blindern, 0315 Oslo, Norway
| | - Tatjana Škrbić
- Department of Physics and Institute for Fundamental Science, University of Oregon, Eugene, Oregon 97403, United States.,Dipartimento di Scienze Molecolari e Nanosistemi, Università Ca' Foscari di Venezia,Campus Scientifico, Edificio Alfa, via Torino 155, 30170 Venezia Mestre, Italy
| | - Sigbjørn Løland Bore
- Department of Chemistry and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, P.O. Box 1033, Blindern, 0315 Oslo, Norway
| | - Giuseppe Milano
- Department of Organic Materials Science, Yamagata University, 4-3-16 Jonan, Yonezawa, 992-8510 Yamagata-ken, Japan.,Dipartimento di Chimica e Biologia, Università di Salerno, Via Giovanni Paolo II 132, 84084 Fisciano, Italy
| | - Michele Cascella
- Department of Chemistry and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, P.O. Box 1033, Blindern, 0315 Oslo, Norway
| | - Achille Giacometti
- Dipartimento di Scienze Molecolari e Nanosistemi, Università Ca' Foscari di Venezia,Campus Scientifico, Edificio Alfa, via Torino 155, 30170 Venezia Mestre, Italy.,European Centre for Living Technology (ECLT) Ca' Bottacin, 3911 Dorsoduro, Calle Crosera, 30123 Venice, Italy
| |
Collapse
|
7
|
Ledum M, Løland Bore S, Cascella M. Automated determination of hybrid particle-field parameters by machine learning. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1785571] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- Morten Ledum
- Department of Chemistry and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, Oslo, Norway
| | - Sigbjørn Løland Bore
- Department of Chemistry and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, Oslo, Norway
| | - Michele Cascella
- Department of Chemistry and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, Oslo, Norway
| |
Collapse
|
8
|
Bore SL, Kolli HB, De Nicola A, Byshkin M, Kawakatsu T, Milano G, Cascella M. Hybrid particle-field molecular dynamics under constant pressure. J Chem Phys 2020; 152:184908. [PMID: 32414244 DOI: 10.1063/5.0007445] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Abstract
Hybrid particle-field methods are computationally efficient approaches for modeling soft matter systems. So far, applications of these methodologies have been limited to constant volume conditions. Here, we reformulate particle-field interactions to represent systems coupled to constant external pressure. First, we show that the commonly used particle-field energy functional can be modified to model and parameterize the isotropic contributions to the pressure tensor without interfering with the microscopic forces on the particles. Second, we employ a square gradient particle-field interaction term to model non-isotropic contributions to the pressure tensor, such as in surface tension phenomena. This formulation is implemented within the hybrid particle-field molecular dynamics approach and is tested on a series of model systems. Simulations of a homogeneous water box demonstrate that it is possible to parameterize the equation of state to reproduce any target density for a given external pressure. Moreover, the same parameterization is transferable to systems of similar coarse-grained mapping resolution. Finally, we evaluate the feasibility of the proposed approach on coarse-grained models of phospholipids, finding that the term between water and the lipid hydrocarbon tails is alone sufficient to reproduce the experimental area per lipid in constant-pressure simulations and to produce a qualitatively correct lateral pressure profile.
Collapse
Affiliation(s)
- Sigbjørn Løland Bore
- Department of Chemistry, and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, P.O. Box 1033, Blindern, 0315 Oslo, Norway
| | - Hima Bindu Kolli
- Department of Chemistry, and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, P.O. Box 1033, Blindern, 0315 Oslo, Norway
| | - Antonio De Nicola
- Department of Organic Materials Science, Yamagata University, 4-3-16 Jonan, Yonezawa, Yamagata-ken 992-8510, Japan
| | - Maksym Byshkin
- Institute of Computational Science, Università Della Svizzera Italiana, 6900 Lugano, Switzerland
| | - Toshihiro Kawakatsu
- Department of Physics, Tohoku University, Aoba, Aramaki, Aoba-ku, Sendai 980-8578, Japan
| | - Giuseppe Milano
- Department of Organic Materials Science, Yamagata University, 4-3-16 Jonan, Yonezawa, Yamagata-ken 992-8510, Japan
| | - Michele Cascella
- Department of Chemistry, and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, P.O. Box 1033, Blindern, 0315 Oslo, Norway
| |
Collapse
|
9
|
Marrink SJ, Corradi V, Souza PC, Ingólfsson HI, Tieleman DP, Sansom MS. Computational Modeling of Realistic Cell Membranes. Chem Rev 2019; 119:6184-6226. [PMID: 30623647 PMCID: PMC6509646 DOI: 10.1021/acs.chemrev.8b00460] [Citation(s) in RCA: 408] [Impact Index Per Article: 81.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2018] [Indexed: 12/15/2022]
Abstract
Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.
Collapse
Affiliation(s)
- Siewert J. Marrink
- Groningen
Biomolecular Sciences and Biotechnology Institute & Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Valentina Corradi
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Paulo C.T. Souza
- Groningen
Biomolecular Sciences and Biotechnology Institute & Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Helgi I. Ingólfsson
- Biosciences
and Biotechnology Division, Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94550, United States
| | - D. Peter Tieleman
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Mark S.P. Sansom
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, U.K.
| |
Collapse
|
10
|
Sgouros AP, Lakkas AT, Megariotis G, Theodorou DN. Mesoscopic Simulations of Free Surfaces of Molten Polyethylene: Brownian Dynamics/Kinetic Monte Carlo Coupled with Square Gradient Theory and Compared to Atomistic Calculations and Experiment. Macromolecules 2018. [DOI: 10.1021/acs.macromol.8b01873] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Affiliation(s)
- A. P. Sgouros
- School of Chemical Engineering, National Technical University of Athens (NTUA), GR-15780 Athens, Greece
| | - A. T. Lakkas
- School of Chemical Engineering, National Technical University of Athens (NTUA), GR-15780 Athens, Greece
| | - G. Megariotis
- School of Chemical Engineering, National Technical University of Athens (NTUA), GR-15780 Athens, Greece
| | - D. N. Theodorou
- School of Chemical Engineering, National Technical University of Athens (NTUA), GR-15780 Athens, Greece
| |
Collapse
|
11
|
Jaschonek S, Cascella M, Gauss J, Diezemann G, Milano G. Intramolecular structural parameters are key modulators of the gel-liquid transition in coarse grained simulations of DPPC and DOPC lipid bilayers. Biochem Biophys Res Commun 2018; 498:327-333. [PMID: 29101041 DOI: 10.1016/j.bbrc.2017.10.132] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2017] [Revised: 09/29/2017] [Accepted: 10/23/2017] [Indexed: 12/11/2022]
Abstract
The capability of coarse-grained models based on the MARTINI mapping to reproduce the gel-liquid phase transition in saturated and unsaturated model lipids was investigated. We found that the model is able to reproduce a lower critical temperature for 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) with respect to 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC). Nonetheless, the appearance of a gel phase for DOPC is strictly dependent on the intramolecular parameters chosen to model its molecular structure. In particular, we show that the bending angle at the coarse-grained bead corresponding to the unsaturated carbon-carbon bond acts as an order parameter determining the temperature of the phase transition. Structural analysis of the molecular dynamics simulations runs evidences that in the gel phase, the packing of the lipophilic tails of DOPC assume a different conformation than in the liquid phase. In the latter phase, the DOPC geometry resembles that of the relaxed free molecule. DPPC:DOPC mixtures show a single phase transition temperature, indicating that the observation of a phase separation between the two lipids requires the simulation of systems with sizes much larger than the ones used here.
Collapse
Affiliation(s)
- Stefan Jaschonek
- Institut für Physikalische Chemie, Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany.
| | - Michele Cascella
- Department of Chemistry and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, Postboks 1033 Blindern, N-0315 Oslo, Norway.
| | - Jürgen Gauss
- Institut für Physikalische Chemie, Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany.
| | - Gregor Diezemann
- Institut für Physikalische Chemie, Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany.
| | - Giuseppe Milano
- Dipartimento di Chimica e Biologia, Università di Salerno, Via Giovanni Paolo II 132, I-84084 Fisciano, Italy.
| |
Collapse
|
12
|
Bore SL, Milano G, Cascella M. Hybrid Particle-Field Model for Conformational Dynamics of Peptide Chains. J Chem Theory Comput 2018; 14:1120-1130. [DOI: 10.1021/acs.jctc.7b01160] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Affiliation(s)
- Sigbjørn Løland Bore
- Department
of Chemistry and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, P.O.
Box 1033 Blindern, 0315 Oslo, Norway
| | - Giuseppe Milano
- Department
of Organic Materials Science, University of Yamagata, 4-3-16 Jonan
Yonezawa, Yamagata-ken 992-8510, Japan
| | - Michele Cascella
- Department
of Chemistry and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, P.O.
Box 1033 Blindern, 0315 Oslo, Norway
| |
Collapse
|
13
|
Qi S, Schmid F. Hybrid particle-continuum simulations coupling Brownian dynamics and local dynamic density functional theory. SOFT MATTER 2017; 13:7938-7947. [PMID: 29034937 DOI: 10.1039/c7sm01749a] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
We present a multiscale hybrid particle-field scheme for the simulation of relaxation and diffusion behavior of soft condensed matter systems. It combines particle-based Brownian dynamics and field-based local dynamics in an adaptive sense such that particles can switch their level of resolution on the fly. The switching of resolution is controlled by a tuning function which can be chosen at will according to the geometry of the system. As an application, the hybrid scheme is used to study the kinetics of interfacial broadening of a polymer blend, and is validated by comparing the results to the predictions from pure Brownian dynamics and pure local dynamics calculations.
Collapse
Affiliation(s)
- Shuanhu Qi
- Institut für Physik, Johannes Gutenberg-Universität Mainz, Staudingerweg 7, D-55099 Mainz, Germany.
| | | |
Collapse
|
14
|
Soares TA, Vanni S, Milano G, Cascella M. Toward Chemically Resolved Computer Simulations of Dynamics and Remodeling of Biological Membranes. J Phys Chem Lett 2017; 8:3586-3594. [PMID: 28707901 DOI: 10.1021/acs.jpclett.7b00493] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Cellular membranes are fundamental constituents of living organisms. Apart from defining the boundaries of the cells, they are involved in a wide range of biological functions, associated with both their structural and the dynamical properties. Biomembranes can undergo large-scale transformations when subject to specific environmental changes, including gel-liquid phase transitions, change of aggregation structure, formation of microtubules, or rupture into vesicles. All of these processes are dependent on a delicate interplay between intermolecular forces, molecular crowding, and entropy, and their understanding requires approaches that are able to capture and rationalize the details of all of the involved interactions. Molecular dynamics-based computational models at atom-level resolution are, in principle, the best way to perform such investigations. Unfortunately, the relevant spatial and time dimensionalities involved in membrane remodeling phenomena would require computational costs that are today unaffordable on a routinely basis. Such hurdles can be removed by coarse-graining the representations of the individual molecular components of the systems. This procedure anyway reduces the possibility of describing the chemical variations in the lipid mixtures composing biological membranes. New hybrid particle field multiscale approaches offer today a promising alternative to the more traditional particle-based simulations methods. By combining chemically distinguishable molecular representations with mesoscale-based computationally affordable potentials, they appear as one of the most promising ways to keep an accurate description of the chemical complexity of biological membranes and, at the same time, cover the required scales to describe remodeling events.
Collapse
Affiliation(s)
- Thereza A Soares
- Department of Fundamental Chemistry, Federal University of Pernambuco, Cidade Universitária , Recife PE 50740-560, Brazil
| | - Stefano Vanni
- Department of Biology, University of Fribourg , 1700 Fribourg, Switzerland
| | - Giuseppe Milano
- Dipartimento di Chimica e Biologia, Università di Salerno , Via Giovanni Paolo II, 132, I-84084 Fisciano, Italy
| | - Michele Cascella
- Department of Chemistry and Centre for Theoretical and Computational Chemistry (CTCC) , Sem Saelands vei 26, 0371 Oslo, Norway
| |
Collapse
|
15
|
Pizzirusso A, De Nicola A, Sevink GJA, Correa A, Cascella M, Kawakatsu T, Rocco M, Zhao Y, Celino M, Milano G. Biomembrane solubilization mechanism by Triton X-100: a computational study of the three stage model. Phys Chem Chem Phys 2017; 19:29780-29794. [DOI: 10.1039/c7cp03871b] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
The solubilization mechanism of lipid membranes in the presence of Triton X-100 (TX-100) is investigated at molecular resolution using hybrid particle field–self consistence field simulations.
Collapse
Affiliation(s)
| | - Antonio De Nicola
- Dipartimento di Chimica e Biologia
- Università di Salerno
- Fisciano
- Italy
| | - G. J. Agur Sevink
- Leiden Institute of Chemistry
- Leiden University
- 2300 RA Leiden
- The Netherlands
| | - Andrea Correa
- Department of Chemical Science
- Federico II University of Naples
- 80126 Napoli
- Italy
| | - Michele Cascella
- Department of Chemistry and Hylleraas Centre for Quantum Molecular Sciences
- University of Oslo
- 0371 Oslo
- Norway
| | | | - Mattia Rocco
- Biopolimeri e Proteomica
- Ospedale Policlinico San Martino
- Genova
- Italy
| | - Ying Zhao
- Institute of Nano-Photonics
- School of Physics and Materials Engineering
- Dalian Minzu University
- Dalian 116600
- China
| | | | - Giuseppe Milano
- Dipartimento di Chimica e Biologia
- Università di Salerno
- Fisciano
- Italy
| |
Collapse
|
16
|
Pizzirusso A, De Nicola A, Milano G. MARTINI Coarse-Grained Model of Triton TX-100 in Pure DPPC Monolayer and Bilayer Interfaces. J Phys Chem B 2016; 120:3821-32. [PMID: 27042862 DOI: 10.1021/acs.jpcb.6b00646] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The coarse-grained MARTINI model of Triton TX-100 has been validated by direct comparison of the experimental and calculated area increase in pure DPPC lipid bilayers and monolayers at water/air interfaces in the presence of surfactant and by comparison of electron density profiles calculated with more detailed atomistic models based on the CHARMM force field. Bilayer simulations have been performed and compared with monolayers and with atomistic models. The validated CG model has been employed to study the phase separation of TX-100 molecules in lipid bilayers and the effect of the lipid bilayer curvature.
Collapse
Affiliation(s)
- Antonio Pizzirusso
- Dipartimento di Chimica e Biologia, Università di Salerno , via Giovanni Paolo II, Fisciano, Salerno I-84084, Italy
| | - Antonio De Nicola
- Dipartimento di Chimica e Biologia, Università di Salerno , via Giovanni Paolo II, Fisciano, Salerno I-84084, Italy
| | - Giuseppe Milano
- Dipartimento di Chimica e Biologia, Università di Salerno , via Giovanni Paolo II, Fisciano, Salerno I-84084, Italy
| |
Collapse
|
17
|
Zhu YL, Lu ZY, Milano G, Shi AC, Sun ZY. Hybrid particle–field molecular dynamics simulation for polyelectrolyte systems. Phys Chem Chem Phys 2016; 18:9799-808. [DOI: 10.1039/c5cp06856h] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
An effective hybrid computer simulation method combining molecular dynamics and self-consistent field theory is developed by including electrostatic interactions.
Collapse
Affiliation(s)
- 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
| | - Giuseppe Milano
- Dipartimento di Chimica e Biologia
- Università degli Studi di Salerno
- Fisciano
- Italy
| | - An-Chang Shi
- Department of Physics and Astronomy
- McMaster University
- Hamilton
- Canada
| | - Zhao-Yan Sun
- State Key Laboratory of Polymer Physics and Chemistry
- Changchun Institute of Applied Chemistry
- Chinese Academy of Sciences
- Changchun 130022
- China
| |
Collapse
|
18
|
De Nicola A, Kawakatsu T, Rosano C, Celino M, Rocco M, Milano G. Self-Assembly of Triton X-100 in Water Solutions: A Multiscale Simulation Study Linking Mesoscale to Atomistic Models. J Chem Theory Comput 2015; 11:4959-71. [DOI: 10.1021/acs.jctc.5b00485] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Antonio De Nicola
- Dipartimento
di Chimica e Biologia, Università di Salerno, I-84084 via
Ponte don Melillo, Salerno, Italy
| | - Toshihiro Kawakatsu
- Department
of Physics, Tohoku University, Aoba, Aramaki, Aoba-ku, Sendai 980-8578, Japan
| | - Camillo Rosano
- Biopolimeri
e Proteomica, IRCCS AOU San Martino-IST, Istituto Nazionale per la Ricerca sul Cancro, Largo R. Benzi 10 I-16132 Genova, Italy
| | - Massimo Celino
- ENEA, C.R. Casaccia, Via Anguillarese 301, I-00123 Roma, Italy
| | - Mattia Rocco
- Biopolimeri
e Proteomica, IRCCS AOU San Martino-IST, Istituto Nazionale per la Ricerca sul Cancro, Largo R. Benzi 10 I-16132 Genova, Italy
| | - Giuseppe Milano
- Dipartimento
di Chimica e Biologia, Università di Salerno, I-84084 via
Ponte don Melillo, Salerno, Italy
- IMAST S.c.a.r.l.
Technological District in Polymer and Composite Engineering, P. leBovio 22, I-80133, Napoli, Italy
| |
Collapse
|
19
|
De Nicola A, Kawakatsu T, Milano G. Generation of Well-Relaxed All-Atom Models of Large Molecular Weight Polymer Melts: A Hybrid Particle-Continuum Approach Based on Particle-Field Molecular Dynamics Simulations. J Chem Theory Comput 2014; 10:5651-67. [PMID: 26583248 DOI: 10.1021/ct500492h] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
A procedure based on Molecular Dynamics (MD) simulations employing soft potentials derived from self-consistent field (SCF) theory (named MD-SCF) able to generate well-relaxed all-atom structures of polymer melts is proposed. All-atom structures having structural correlations indistinguishable from ones obtained by long MD relaxations have been obtained for poly(methyl methacrylate) (PMMA) and poly(ethylene oxide) (PEO) melts. The proposed procedure leads to computational costs mainly related on system size rather than to the chain length. Several advantages of the proposed procedure over current coarse-graining/reverse mapping strategies are apparent. No parametrization is needed to generate relaxed structures of different polymers at different scales or resolutions. There is no need for special algorithms or back-mapping schemes to change the resolution of the models. This characteristic makes the procedure general and its extension to other polymer architectures straightforward. A similar procedure can be easily extended to the generation of all-atom structures of block copolymer melts and polymer nanocomposites.
Collapse
Affiliation(s)
- Antonio De Nicola
- Dipartimento di Chimica e Biologia, Università degli Studi di Salerno , via Ponte don Melillo, Fisciano, Salerno I-84085, Italy.,IMAST Scarl-Technological District in Polymer and Composite Engineering, Piazza Bovio 22, Napoli, Napoli I-80133, Italy
| | - Toshihiro Kawakatsu
- Department of Physics, Tohoku University , Aoba, Aramaki, Aoba-ku, Sendai, Miyagi 980-8578, Japan
| | - Giuseppe Milano
- Dipartimento di Chimica e Biologia, Università degli Studi di Salerno , via Ponte don Melillo, Fisciano, Salerno I-84085, Italy.,IMAST Scarl-Technological District in Polymer and Composite Engineering, Piazza Bovio 22, Napoli, Napoli I-80133, Italy
| |
Collapse
|
20
|
Sarukhanyan E, De Nicola A, Roccatano D, Kawakatsu T, Milano G. Spontaneous insertion of carbon nanotube bundles inside biomembranes: A hybrid particle-field coarse-grained molecular dynamics study. Chem Phys Lett 2014. [DOI: 10.1016/j.cplett.2014.01.057] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
21
|
De Nicola A, Hezaveh S, Zhao Y, Kawakatsu T, Roccatano D, Milano G. Micellar drug nanocarriers and biomembranes: how do they interact? Phys Chem Chem Phys 2014; 16:5093-105. [DOI: 10.1039/c3cp54242d] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|