1
|
Lowry GV, Giraldo JP, Steinmetz NF, Avellan A, Demirer GS, Ristroph KD, Wang GJ, Hendren CO, Alabi CA, Caparco A, da Silva W, González-Gamboa I, Grieger KD, Jeon SJ, Khodakovskaya MV, Kohay H, Kumar V, Muthuramalingam R, Poffenbarger H, Santra S, Tilton RD, White JC. Towards realizing nano-enabled precision delivery in plants. NATURE NANOTECHNOLOGY 2024; 19:1255-1269. [PMID: 38844663 DOI: 10.1038/s41565-024-01667-5] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/22/2023] [Accepted: 03/27/2024] [Indexed: 09/18/2024]
Abstract
Nanocarriers (NCs) that can precisely deliver active agents, nutrients and genetic materials into plants will make crop agriculture more resilient to climate change and sustainable. As a research field, nano-agriculture is still developing, with significant scientific and societal barriers to overcome. In this Review, we argue that lessons can be learned from mammalian nanomedicine. In particular, it may be possible to enhance efficiency and efficacy by improving our understanding of how NC properties affect their interactions with plant surfaces and biomolecules, and their ability to carry and deliver cargo to specific locations. New tools are required to rapidly assess NC-plant interactions and to explore and verify the range of viable targeting approaches in plants. Elucidating these interactions can lead to the creation of computer-generated in silico models (digital twins) to predict the impact of different NC and plant properties, biological responses, and environmental conditions on the efficiency and efficacy of nanotechnology approaches. Finally, we highlight the need for nano-agriculture researchers and social scientists to converge in order to develop sustainable, safe and socially acceptable NCs.
Collapse
Affiliation(s)
- Gregory V Lowry
- Civil and Environmental Engineering, Carnegie Mellon University, Pittsburgh, PA, USA.
| | - Juan Pablo Giraldo
- Botany and Plant Sciences, University of California, Riverside, Riverside, CA, USA.
| | - Nicole F Steinmetz
- Department of NanoEngineering, University of California San Diego, San Diego, CA, USA
- Department of Bioengineering, University of California San Diego, San Diego, CA, USA
- Department of Radiology, University of California San Diego, San Diego, CA, USA
- Center for Nano-ImmunoEngineering, University of California San Diego, San Diego, CA, USA
- Shu and K.C. Chien and Peter Farrell Collaboratory, University of California San Diego, San Diego, CA, USA
- Center for Engineering in Cancer, Institute of Engineering in Medicine, University of California San Diego, San Diego, CA, USA
- Moores Cancer Center, University of California, University of California San Diego, San Diego, CA, USA
- Institute for Materials Discovery and Design, University of California San Diego, San Diego, CA, USA
| | | | - Gozde S Demirer
- Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA, USA
| | - Kurt D Ristroph
- Agricultural and Biological Engineering, Purdue University, West Lafayette, IN, USA
| | - Gerald J Wang
- Civil and Environmental Engineering, Carnegie Mellon University, Pittsburgh, PA, USA
| | - Christine O Hendren
- Geological and Environmental Sciences, Appalachian State University, Boone, NC, USA
| | | | - Adam Caparco
- Department of NanoEngineering, University of California San Diego, San Diego, CA, USA
| | | | | | - Khara D Grieger
- Applied Ecology, North Carolina State University, Raleigh, NC, USA
| | - Su-Ji Jeon
- Botany and Plant Sciences, University of California, Riverside, Riverside, CA, USA
| | | | - Hagay Kohay
- Civil and Environmental Engineering, Carnegie Mellon University, Pittsburgh, PA, USA
| | - Vivek Kumar
- Civil and Environmental Engineering, Carnegie Mellon University, Pittsburgh, PA, USA
| | | | | | - Swadeshmukul Santra
- Department of Chemistry and Burnett School of Biomedical Sciences, University of Central Florida, Orlando, FL, USA
| | - Robert D Tilton
- Chemical Engineering and Biomedical Engineering, Carnegie Mellon University, Pittsburgh, PA, USA
| | - Jason C White
- The Connecticut Agricultural Research Station, New Haven, CT, USA
| |
Collapse
|
2
|
Jin J, Voth GA. Understanding dynamics in coarse-grained models. IV. Connection of fine-grained and coarse-grained dynamics with the Stokes-Einstein and Stokes-Einstein-Debye relations. J Chem Phys 2024; 161:034114. [PMID: 39012809 DOI: 10.1063/5.0212973] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2024] [Accepted: 06/28/2024] [Indexed: 07/18/2024] Open
Abstract
Applying an excess entropy scaling formalism to the coarse-grained (CG) dynamics of liquids, we discovered that missing rotational motions during the CG process are responsible for artificially accelerated CG dynamics. In the context of the dynamic representability between the fine-grained (FG) and CG dynamics, this work introduces the well-known Stokes-Einstein and Stokes-Einstein-Debye relations to unravel the rotational dynamics underlying FG trajectories, thereby allowing for an indirect evaluation of the effective rotations based only on the translational information at the reduced CG resolution. Since the representability issue in CG modeling limits a direct evaluation of the shear stress appearing in the Stokes-Einstein and Stokes-Einstein-Debye relations, we introduce a translational relaxation time as a proxy to employ these relations, and we demonstrate that these relations hold for the ambient conditions studied in our series of work. Additional theoretical links to our previous work are also established. First, we demonstrate that the effective hard sphere radius determined by the classical perturbation theory can approximate the complex hydrodynamic radius value reasonably well. Furthermore, we present a simple derivation of an excess entropy scaling relationship for viscosity by estimating the elliptical integral of molecules. In turn, since the translational and rotational motions at the FG level are correlated to each other, we conclude that the "entropy-free" CG diffusion only depends on the shape of the reference molecule. Our results and analyses impart an alternative way of recovering the FG diffusion from the CG description by coupling the translational and rotational motions at the hydrodynamic level.
Collapse
Affiliation(s)
- Jaehyeok Jin
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | - Gregory A Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
3
|
Malajczuk CJ, Mancera RL. An atomistic characterization of high-density lipoproteins and the conserved "LN" region of apoA-I. Biophys J 2024; 123:1116-1128. [PMID: 38555508 PMCID: PMC11079945 DOI: 10.1016/j.bpj.2024.03.039] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2023] [Revised: 03/19/2024] [Accepted: 03/27/2024] [Indexed: 04/02/2024] Open
Abstract
The physicochemical characteristics of the various subpopulations of high-density lipoproteins (HDLs) and, in particular, their surface properties determine their ability to scavenge lipids and interact with specific receptors and peptides. Five representative spheroidal HDL subpopulation models were mapped from a previously reported equilibrated coarse-grained (CG) description to an atomistic representation for subsequent molecular dynamics simulation. For each HDL model a range of finer-level analyses was undertaken, including the component-wise characterization of HDL surfaces, the average size and composition of hydrophobic surface patches, dynamic protein secondary structure monitoring, and the proclivity for solvent exposure of the proposed β-amyloid (Aβ) binding region of apolipoprotein A-I (apoA-I), "LN." This study reveals that previously characterized ellipsoidal HDL3a and HDL2a models revert to a more spherical geometry in an atomistic representation due to the enhanced conformational flexibility afforded to the apoA-I protein secondary structure, allowing for enhanced surface lipid packing and lower overall surface hydrophobicity. Indeed, the proportional surface hydrophobicity and apoA-I exposure reduced with increasing HDL size, consistent with previous characterizations. Furthermore, solvent exposure of the "LN" region of apoA-I was exclusively limited to the smallest HDL3c model within the timescale of the simulations, and typically corresponded to a distinct loss in secondary structure across the "LN" region to form part of a significant contiguous hydrophobic patch on the HDL surface. Taken together, these findings provide preliminary evidence for a subpopulation-specific interaction between HDL3c particles and circulating hydrophobic species such as Aβ via the exposed "LN" region of apoA-I.
Collapse
Affiliation(s)
- Chris J Malajczuk
- Curtin Medical School, Curtin Health Innovation Research Institute and Curtin Institute for Data Science, Curtin University, Perth, WA, Australia
| | - Ricardo L Mancera
- Curtin Medical School, Curtin Health Innovation Research Institute and Curtin Institute for Data Science, Curtin University, Perth, WA, Australia.
| |
Collapse
|
4
|
Jin J, Reichman DR. Hierarchical Framework for Predicting Entropies in Bottom-Up Coarse-Grained Models. J Phys Chem B 2024; 128:3182-3199. [PMID: 38507575 DOI: 10.1021/acs.jpcb.3c07624] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/22/2024]
Abstract
The thermodynamic entropy of coarse-grained (CG) models stands as one of the most important properties for quantifying the missing information during the CG process and for establishing transferable (or extendible) CG interactions. However, performing additional CG simulations on top of model construction often leads to significant additional computational overhead. In this work, we propose a simple hierarchical framework for predicting the thermodynamic entropies of various molecular CG systems. Our approach employs a decomposition of the CG interactions, enabling the estimation of the CG partition function and thermodynamic properties a priori. Starting from the ideal gas description, we leverage classical perturbation theory to systematically incorporate simple yet essential interactions, ranging from the hard sphere model to the generalized van der Waals model. Additionally, we propose an alternative approach based on multiparticle correlation functions, allowing for systematic improvements through higher-order correlations. Numerical applications to molecular liquids validate the high fidelity of our approach, and our computational protocols demonstrate that a reduced model with simple energetics can reasonably estimate the thermodynamic entropy of CG models without performing any CG simulations. Overall, our findings present a systematic framework for estimating not only the entropy but also other thermodynamic properties of CG models, relying solely on information from the reference system.
Collapse
Affiliation(s)
- Jaehyeok Jin
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - David R Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
5
|
Lesniewski MC, Noid WG. Insight into the Density-Dependence of Pair Potentials for Predictive Coarse-Grained Models. J Phys Chem B 2024; 128:1298-1316. [PMID: 38271676 DOI: 10.1021/acs.jpcb.3c06890] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2024]
Abstract
We investigate the temperature- and density-dependence of effective pair potentials for 1-site coarse-grained (CG) models of two industrial solvents, 1,4-dioxane and tetrahydrofuran. We observe that the calculated pair potentials are much more sensitive to density than to temperature. The generalized-Yvon-Born-Green framework reveals that this striking density-dependence reflects corresponding variations in the many-body correlations that determine the environment-mediated indirect contribution to the pair mean force. Moreover, we demonstrate, perhaps surprisingly, that this density-dependence is not important for accurately modeling the intermolecular structure. Accordingly, we adopt a density-independent interaction potential and transfer the density-dependence of the calculated pair potentials into a configuration-independent volume potential. Furthermore, we develop a single global potential that accurately models the intermolecular structure and pressure-volume equation of state across a very wide range of liquid state points. Consequently, this work provides fundamental insight into the density-dependence of effective pair potentials and also provides a significant step toward developing predictive CG models for efficiently modeling industrial solvents.
Collapse
Affiliation(s)
- Maria C Lesniewski
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - W G Noid
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| |
Collapse
|
6
|
Chauhan A, Gogoi D, Puri S, Singh A. Effect of amphiphilic polymers on phase separating binary mixtures: A DPD simulation study. J Chem Phys 2023; 159:204901. [PMID: 37991159 DOI: 10.1063/5.0173817] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2023] [Accepted: 11/02/2023] [Indexed: 11/23/2023] Open
Abstract
We present the phase separation dynamics of a binary (AB), simple fluid (SF), and amphiphilic polymer (AP) mixture using dissipative particle dynamics simulation at d = 3. We study the effect of different AP topologies, including block copolymers, ring block copolymers (RCP), and miktoarm star polymers, on the evolution morphologies, dynamic scaling functions, and length scale of the AB mixture. Our results demonstrate that the presence of APs leads to significantly different evolution morphologies in SF. However, the deviation from dynamical scaling is prominent, mainly for RCP. Typically, the characteristic length scale for SF follows the power law R(t) ∼ tϕ, where ϕ is the growth exponent. In the presence of high AP, we observe diffusive growth (ϕ → 1/3) at early times, followed by saturation in length scale (ϕ → 0) at late times. The extent of saturation varies with constraints imposed on the APs, such as topology, composition ratio, chain length, and stiffness. At lower composition ratios, the system exhibits inertial hydrodynamic growth (ϕ → 2/3) at asymptotic times without clearly exhibiting the viscous hydrodynamic regime (ϕ → 1) at earlier times in our simulations. Our results firmly establish the existence of hydrodynamic growth regimes in low surfactant-influenced phase separation kinetics of binary fluids and settle the related ambiguity in d = 3 systems.
Collapse
Affiliation(s)
- Avinash Chauhan
- Department of Physics, Indian Institute of Technology (BHU), Varanasi, Uttar Pradesh 221005, India
| | - Dorothy Gogoi
- School of Physical Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| | - Sanjay Puri
- School of Physical Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| | - Awaneesh Singh
- Department of Physics, Indian Institute of Technology (BHU), Varanasi, Uttar Pradesh 221005, India
| |
Collapse
|
7
|
Jin J, Hwang J, Voth GA. Gaussian representation of coarse-grained interactions of liquids: Theory, parametrization, and transferability. J Chem Phys 2023; 159:184105. [PMID: 37942867 DOI: 10.1063/5.0160567] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2023] [Accepted: 10/06/2023] [Indexed: 11/10/2023] Open
Abstract
Coarse-grained (CG) interactions determined via bottom-up methodologies can faithfully reproduce the structural correlations observed in fine-grained (atomistic resolution) systems, yet they can suffer from limited extensibility due to complex many-body correlations. As part of an ongoing effort to understand and improve the applicability of bottom-up CG models, we propose an alternative approach to address both accuracy and transferability. Our main idea draws from classical perturbation theory to partition the hard sphere repulsive term from effective CG interactions. We then introduce Gaussian basis functions corresponding to the system's characteristic length by linking these Gaussian sub-interactions to the local particle densities at each coordination shell. The remaining perturbative long-range interaction can be treated as a collective solvation interaction, which we show exhibits a Gaussian form derived from integral equation theories. By applying this numerical parametrization protocol to CG liquid systems, our microscopic theory elucidates the emergence of Gaussian interactions in common phenomenological CG models. To facilitate transferability for these reduced descriptions, we further infer equations of state to determine the sub-interaction parameter as a function of the system variables. The reduced models exhibit excellent transferability across the thermodynamic state points. Furthermore, we propose a new strategy to design the cross-interactions between distinct CG sites in liquid mixtures. This involves combining each Gaussian in the proper radial domain, yielding accurate CG potentials of mean force and structural correlations for multi-component systems. Overall, our findings establish a solid foundation for constructing transferable bottom-up CG models of liquids with enhanced extensibility.
Collapse
Affiliation(s)
- Jaehyeok Jin
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, 5735 S. Ellis Ave., Chicago, Illinois 60637, USA
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| | - Jisung Hwang
- Department of Statistics, The University of Chicago, 5747 S. Ellis Ave., Chicago, Illinois 60637, USA
| | - Gregory A Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, 5735 S. Ellis Ave., Chicago, Illinois 60637, USA
| |
Collapse
|
8
|
Kehrein J, Gürsöz E, Davies M, Luxenhofer R, Bunker A. Unravel the Tangle: Atomistic Insight into Ultrahigh Curcumin-Loaded Polymer Micelles. SMALL (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2023; 19:e2303066. [PMID: 37403298 DOI: 10.1002/smll.202303066] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/11/2023] [Revised: 06/22/2023] [Indexed: 07/06/2023]
Abstract
Amphiphilic ABA-triblock copolymers, comprised of poly(2-oxazoline) and poly(2-oxazine), can solubilize poorly water-soluble molecules in a structure-dependent manner forming micelles with exceptionally high drug loading. All-atom molecular dynamics simulations are conducted on previously experimentally characterized, curcumin-loaded micelles to dissect the structure-property relationships. Polymer-drug interactions for different levels of drug loading and variation in polymer structures of both the inner hydrophobic core and outer hydrophilic shell are investigated. In silico, the system with the highest experimental loading capacity shows the highest number of drug molecules encapsulated by the core. Furthermore, in systems with lower loading capacity outer A blocks show a greater extent of entanglement with the inner B blocks. Hydrogen bond analyses corroborate previous hypotheses: poly(2-butyl-2-oxazoline) B blocks, found experimentally to have reduced loading capacity for curcumin compared to poly(2-propyl-2-oxazine), establish fewer but longer-lasting hydrogen bonds. This possibly results from different sidechain conformations around the hydrophobic cargo, which is investigated by unsupervised machine learning to cluster monomers in smaller model systems mimicking different micelle compartments. Exchanging poly(2-methyl-2-oxazoline) with poly(2-ethyl-2-oxazoline) leads to increased drug interactions and reduced corona hydration; this suggests an impairment of micelle solubility or colloidal stability. These observations can help driving forward a more rational a priori nanoformulation design.
Collapse
Affiliation(s)
- Josef Kehrein
- Soft Matter Chemistry, Department of Chemistry, Faculty of Science, University of Helsinki, Helsinki, 00014, Finland
- Division of Pharmaceutical Biosciences, Drug Research Program, Faculty of Pharmacy, University of Helsinki, Helsinki, 00014, Finland
| | - Ekinsu Gürsöz
- Soft Matter Chemistry, Department of Chemistry, Faculty of Science, University of Helsinki, Helsinki, 00014, Finland
- Division of Pharmaceutical Biosciences, Drug Research Program, Faculty of Pharmacy, University of Helsinki, Helsinki, 00014, Finland
| | - Matthew Davies
- Department of Physics and Astronomy, The University of Western Ontario, 1151 Richmond Street, London, Ontario, N6A 5B7, Canada
| | - Robert Luxenhofer
- Soft Matter Chemistry, Department of Chemistry, Faculty of Science, University of Helsinki, Helsinki, 00014, Finland
| | - Alex Bunker
- Division of Pharmaceutical Biosciences, Drug Research Program, Faculty of Pharmacy, University of Helsinki, Helsinki, 00014, Finland
| |
Collapse
|
9
|
Wang Y, Wang K, Zhao X, Xu X, Sun T. Influence of pH on the self-assembly of diphenylalanine peptides: molecular insights from coarse-grained simulations. SOFT MATTER 2023; 19:5749-5757. [PMID: 37462931 DOI: 10.1039/d3sm00739a] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/03/2023]
Abstract
Nanostructures fabricated from peptide self-assemblies are attracting increasing attention owing to their possible applications in biology and nanotechnology. A known example is an aromatic dipeptide (diphenylalanine, FF) which is extracted from Alzheimer's β-amyloid polypeptide as the core recognition motif for molecular self-assembly. Many studies have been carried out to organize FF peptides into various functional ordered nanostructures. For potential applications of self-assembled FF-based nanomaterials, it becomes important to consider some influencing factors (e.g., solvents, peptide concentrations, pH, temperature, etc.) on the self-assembly process. Among these factors, the effect of pH on the self-assembly process of FF peptides into assembled nanostructures through simulation studies is the main focus of the present work. In the current study, we have investigated the assembly pathway of 1000 FF peptides and qualitatively evaluated the morphological changes of FF-based nanostructures at different pH values by performing extensive coarse-grained molecular dynamics (CG-MD) simulations. Structural analyses suggest that FF peptides can spontaneously assemble into nanotubes with different shapes under acidic, neutral and basic conditions. Based on the analysis of FF nanostructure formation pathways in different pH solutions, the self-assembly of the nanotube involves the aggregation of molecules to form a bilayer, the curling of a bilayer to form a vesicle and the transformation of a vesicle into a tubular structure. It is noted that a flat hollow columnar structure is observed as a special intermediate state during the transformation process of a vesicle-like to a tube-like structure. Energetic analysis suggests that the aggregation of FF peptides is driven by the vdW interactions but the aggregation shape is mainly affected by the electrostatic interactions. Overall, this study provides further understanding of the self-assembly behavior of aromatic short peptide derivatives in different pH solutions.
Collapse
Affiliation(s)
- Yan Wang
- Department of Physics, Zhejiang University of Science and Technology, Hangzhou, Zhejiang 310008, China.
| | - Kang Wang
- Department of Physics, Zhejiang University of Science and Technology, Hangzhou, Zhejiang 310008, China.
| | - Xinyi Zhao
- Department of Physics, Zhejiang University of Science and Technology, Hangzhou, Zhejiang 310008, China.
| | - Xiaojun Xu
- Institute of Bioinformatics and Medical Engineering, Jiangsu University of Technology, Changzhou, Jiangsu 213001, China
| | - Tingting Sun
- Department of Physics, Zhejiang University of Science and Technology, Hangzhou, Zhejiang 310008, China.
| |
Collapse
|
10
|
Jin J, Voth GA. Statistical Mechanical Design Principles for Coarse-Grained Interactions across Different Conformational Free Energy Surfaces. J Phys Chem Lett 2023; 14:1354-1362. [PMID: 36728761 PMCID: PMC9940719 DOI: 10.1021/acs.jpclett.2c03844] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2022] [Accepted: 01/27/2023] [Indexed: 06/18/2023]
Abstract
Systematic bottom-up coarse-graining (CG) of molecular systems provides a means to explore different coupled length and time scales while treating the molecular-scale physics at a reduced level. However, the configuration dependence of CG interactions often results in CG models with limited applicability for exploring the parametrized configurations. We propose a statistical mechanical theory to design CG interactions across different configurations and conditions. In order to span wide ranges of conformational space, distinct classical CG free energy surfaces for characteristic configurations are identified using molecular collective variables. The coupling interaction between different CG free energy surfaces can then be systematically determined by analogy to quantum mechanical approaches describing coupled states. The present theory can accurately capture the underlying many-body potentials of mean force in the CG variables for various order parameters applied to liquids, interfaces, and in principle proteins, uncovering the complex nature underlying the coupling interaction and imparting a new protocol for the design of predictive multiscale models.
Collapse
Affiliation(s)
| | - Gregory A. Voth
- Department of Chemistry,
Chicago Center for Theoretical Chemistry, Institute for Biophysical
Dynamics, and James Franck Institute, The
University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
11
|
Jin J, Schweizer KS, Voth GA. Understanding dynamics in coarse-grained models. II. Coarse-grained diffusion modeled using hard sphere theory. J Chem Phys 2023; 158:034104. [PMID: 36681632 DOI: 10.1063/5.0116300] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
The first paper of this series [J. Chem. Phys. 158, 034103 (2023)] demonstrated that excess entropy scaling holds for both fine-grained and corresponding coarse-grained (CG) systems. Despite its universality, a more exact determination of the scaling relationship was not possible due to the semi-empirical nature. In this second paper, an analytical excess entropy scaling relation is derived for bottom-up CG systems. At the single-site CG resolution, effective hard sphere systems are constructed that yield near-identical dynamical properties as the target CG systems by taking advantage of how hard sphere dynamics and excess entropy can be analytically expressed in terms of the liquid packing fraction. Inspired by classical equilibrium perturbation theories and recent advances in constructing hard sphere models for predicting activated dynamics of supercooled liquids, we propose a new approach for understanding the diffusion of molecular liquids in the normal regime using hard sphere reference fluids. The proposed "fluctuation matching" is designed to have the same amplitude of long wavelength density fluctuations (dimensionless compressibility) as the CG system. Utilizing the Enskog theory to derive an expression for hard sphere diffusion coefficients, a bridge between the CG dynamics and excess entropy is then established. The CG diffusion coefficient can be roughly estimated using various equations of the state, and an accurate prediction of accelerated CG dynamics at different temperatures is also possible in advance of running any CG simulation. By introducing another layer of coarsening, these findings provide a more rigorous method to assess excess entropy scaling and understand the accelerated CG dynamics of molecular fluids.
Collapse
Affiliation(s)
- Jaehyeok Jin
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - Kenneth S Schweizer
- Department of Material Science, Department of Chemistry, Department of Chemical and Biomolecular Engineering, and Materials Research Laboratory, University of Illinois, Urbana, Illinois 61801, USA
| | - Gregory A Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
12
|
Jin J, Schweizer KS, Voth GA. Understanding dynamics in coarse-grained models. I. Universal excess entropy scaling relationship. J Chem Phys 2023; 158:034103. [PMID: 36681649 DOI: 10.1063/5.0116299] [Citation(s) in RCA: 15] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Coarse-grained (CG) models facilitate an efficient exploration of complex systems by reducing the unnecessary degrees of freedom of the fine-grained (FG) system while recapitulating major structural correlations. Unlike structural properties, assessing dynamic properties in CG modeling is often unfeasible due to the accelerated dynamics of the CG models, which allows for more efficient structural sampling. Therefore, the ultimate goal of the present series of articles is to establish a better correspondence between the FG and CG dynamics. To assess and compare dynamical properties in the FG and the corresponding CG models, we utilize the excess entropy scaling relationship. For Paper I of this series, we provide evidence that the FG and the corresponding CG counterpart follow the same universal scaling relationship. By carefully reviewing and examining the literature, we develop a new theory to calculate excess entropies for the FG and CG systems while accounting for entropy representability. We demonstrate that the excess entropy scaling idea can be readily applied to liquid water and methanol systems at both the FG and CG resolutions. For both liquids, we reveal that the scaling exponents remain unchanged from the coarse-graining process, indicating that the scaling behavior is universal for the same underlying molecular systems. Combining this finding with the concept of mapping entropy in CG models, we show that the missing entropy plays an important role in accelerating the CG dynamics.
Collapse
Affiliation(s)
- Jaehyeok Jin
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - Kenneth S Schweizer
- Department of Material Science, Department of Chemistry, Department of Chemical and Biomolecular Engineering, and Materials Research Laboratory, University of Illinois, Urbana, Illinois 61801, USA
| | - Gregory A Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
13
|
Computing dissipative particle dynamics interactions to render molecular structure and temperature-dependent properties of simple liquids. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120539] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
14
|
Christofi E, Chazirakis A, Chrysostomou C, Nicolaou MA, Li W, Doxastakis M, Harmandaris VA. Deep convolutional neural networks for generating atomistic configurations of multi-component macromolecules from coarse-grained models. J Chem Phys 2022; 157:184903. [DOI: 10.1063/5.0110322] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Despite the modern advances in the available computational resources, the length and time scales of the physical systems that can be studied in full atomic detail, via molecular simulations, are still limited. To overcome such limitations, coarse-grained (CG) models have been developed to reduce the dimensionality of the physical system under study. However, to study such systems at the atomic level, it is necessary to re-introduce the atomistic details into the CG description. Such an ill-posed mathematical problem is typically treated via numerical algorithms, which need to balance accuracy, efficiency, and general applicability. Here, we introduce an efficient and versatile method for backmapping multi-component CG macromolecules of arbitrary microstructures. By utilizing deep learning algorithms, we train a convolutional neural network to learn structural correlations between polymer configurations at the atomistic and their corresponding CG descriptions, obtained from atomistic simulations. The trained model is then utilized to get predictions of atomistic structures from input CG configurations. As an illustrative example, we apply the convolutional neural network to polybutadiene copolymers of various microstructures, in which each monomer microstructure (i.e., cis-1,4, trans-1,4, and vinyl-1,2) is represented as a different CG particle type. The proposed methodology is transferable over molecular weight and various microstructures. Moreover, starting from a specific single CG configuration with a given microstructure, we show that by modifying its chemistry (i.e., CG particle types), we are able to obtain a set of well equilibrated polymer configurations of different microstructures (chemistry) than the one of the original CG configuration.
Collapse
Affiliation(s)
- Eleftherios Christofi
- Computation-based Science and Technology Research Center, The Cyprus Institute, Nicosia 2121, Cyprus
| | - Antonis Chazirakis
- Department of Mathematics and Applied Mathematics, University of Crete, Heraklion GR-71110, Greece
- Institute of Applied and Computational Mathematics, Foundation for Research and Technology–Hellas, GR-71110 Heraklion, Crete, Greece
| | - Charalambos Chrysostomou
- Computation-based Science and Technology Research Center, The Cyprus Institute, Nicosia 2121, Cyprus
| | - Mihalis A. Nicolaou
- Computation-based Science and Technology Research Center, The Cyprus Institute, Nicosia 2121, Cyprus
| | - Wei Li
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo 001-0021, Japan
| | - Manolis Doxastakis
- Department of Chemical and Biomolecular Engineering, University of Tennessee, Knoxville, Tennessee 37996, USA
| | - Vagelis A. Harmandaris
- Computation-based Science and Technology Research Center, The Cyprus Institute, Nicosia 2121, Cyprus
- Department of Mathematics and Applied Mathematics, University of Crete, Heraklion GR-71110, Greece
- Institute of Applied and Computational Mathematics, Foundation for Research and Technology–Hellas, GR-71110 Heraklion, Crete, Greece
| |
Collapse
|
15
|
Stieffenhofer M, Scherer C, May F, Bereau T, Andrienko D. Benchmarking coarse-grained models of organic semiconductors via deep backmapping. Front Chem 2022; 10:982757. [PMID: 36157043 PMCID: PMC9500322 DOI: 10.3389/fchem.2022.982757] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2022] [Accepted: 07/22/2022] [Indexed: 11/13/2022] Open
Abstract
The potential of mean force is an effective coarse-grained potential, which is often approximated by pairwise potentials. While the approximated potential reproduces certain distributions of the reference all-atom model with remarkable accuracy, important cross-correlations are typically not captured. In general, the quality of coarse-grained models is evaluated at the coarse-grained resolution, hindering the detection of important discrepancies between the all-atom and coarse-grained ensembles. In this work, the quality of different coarse-grained models is assessed at the atomistic resolution deploying reverse-mapping strategies. In particular, coarse-grained structures for Tris-Meta-Biphenyl-Triazine are reverse-mapped from two different sources: 1) All-atom configurations projected onto the coarse-grained resolution and 2) snapshots obtained by molecular dynamics simulations based on the coarse-grained force fields. To assess the quality of the coarse-grained models, reverse-mapped structures of both sources are compared revealing significant discrepancies between the all-atom and the coarse-grained ensembles. Specifically, the reintroduced details enable force computations based on the all-atom force field that yield a clear ranking for the quality of the different coarse-grained models.
Collapse
Affiliation(s)
| | | | | | - Tristan Bereau
- Van ‘t Hoff Institute for Molecular Sciences and Informatics Institute, University of Amsterdam, Amsterdam, Netherlands
| | - Denis Andrienko
- Max Planck Institute for Polymer Research, Mainz, Germany
- *Correspondence: Denis Andrienko,
| |
Collapse
|
16
|
Karibayev M, Kalybekkyzy S, Wang Y, Mentbayeva A. Molecular Modeling in Anion Exchange Membrane Research: A Brief Review of Recent Applications. Molecules 2022; 27:3574. [PMID: 35684512 PMCID: PMC9182285 DOI: 10.3390/molecules27113574] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2022] [Revised: 05/24/2022] [Accepted: 05/30/2022] [Indexed: 12/04/2022] Open
Abstract
Anion Exchange Membrane (AEM) fuel cells have attracted growing interest, due to their encouraging advantages, including high power density and relatively low cost. AEM is a polymer matrix, which conducts hydroxide (OH-) ions, prevents physical contact of electrodes, and has positively charged head groups (mainly quaternary ammonium (QA) groups), covalently bound to the polymer backbone. The chemical instability of the quaternary ammonium (QA)-based head groups, at alkaline pH and elevated temperature, is a significant threshold in AEMFC technology. This review work aims to introduce recent studies on the chemical stability of various QA-based head groups and transportation of OH- ions in AEMFC, via modeling and simulation techniques, at different scales. It starts by introducing the fundamental theories behind AEM-based fuel-cell technology. In the main body of this review, we present selected computational studies that deal with the effects of various parameters on AEMs, via a variety of multi-length and multi-time-scale modeling and simulation methods. Such methods include electronic structure calculations via the quantum Density Functional Theory (DFT), ab initio, classical all-atom Molecular Dynamics (MD) simulations, and coarse-grained MD simulations. The explored processing and structural parameters include temperature, hydration levels, several QA-based head groups, various types of QA-based head groups and backbones, etc. Nowadays, many methods and software packages for molecular and materials modeling are available. Applications of such methods may help to understand the transportation mechanisms of OH- ions, the chemical stability of functional head groups, and many other relevant properties, leading to a performance-based molecular and structure design as well as, ultimately, improved AEM-based fuel cell performances. This contribution aims to introduce those molecular modeling methods and their recent applications to the AEM-based fuel cells research community.
Collapse
Affiliation(s)
- Mirat Karibayev
- Department of Chemical & Materials Engineering, School of Engineering and Digital Sciences, Nazarbayev University, Nur-Sultan 010000, Kazakhstan;
| | - Sandugash Kalybekkyzy
- Laboratory of Advanced Materials and Systems for Energy Storage, Center for Energy and Advanced Materials Science, National Laboratory Astana, Nazarbayev University, Nur-Sultan 010000, Kazakhstan;
| | - Yanwei Wang
- Department of Chemical & Materials Engineering, School of Engineering and Digital Sciences, Nazarbayev University, Nur-Sultan 010000, Kazakhstan;
- Laboratory of Computational Materials Science for Energy Applications, Center for Energy and Advanced Materials Science, National Laboratory Astana, Nur-Sultan 010000, Kazakhstan
| | - Almagul Mentbayeva
- Department of Chemical & Materials Engineering, School of Engineering and Digital Sciences, Nazarbayev University, Nur-Sultan 010000, Kazakhstan;
- Laboratory of Advanced Materials and Systems for Energy Storage, Center for Energy and Advanced Materials Science, National Laboratory Astana, Nazarbayev University, Nur-Sultan 010000, Kazakhstan;
| |
Collapse
|
17
|
Hernández-Del-Valle M, Valencia-Expósito A, López-Izquierdo A, Casanova-Ferrer P, Tarazona P, Martín-Bermudo MD, Míguez DG. A coarse-grained approach to model the dynamics of the actomyosin cortex. BMC Biol 2022; 20:90. [PMID: 35459165 PMCID: PMC9034637 DOI: 10.1186/s12915-022-01279-2] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2021] [Accepted: 03/11/2022] [Indexed: 01/21/2023] Open
Abstract
Background The dynamics of the actomyosin machinery is at the core of many important biological processes. Several relevant cellular responses such as the rhythmic compression of the cell cortex are governed, at a mesoscopic level, by the nonlinear interaction between actin monomers, actin crosslinkers, and myosin motors. Coarse-grained models are an optimal tool to study actomyosin systems, since they can include processes that occur at long time and space scales, while maintaining the most relevant features of the molecular interactions. Results Here, we present a coarse-grained model of a two-dimensional actomyosin cortex, adjacent to a three-dimensional cytoplasm. Our simplified model incorporates only well-characterized interactions between actin monomers, actin crosslinkers and myosin, and it is able to reproduce many of the most important aspects of actin filament and actomyosin network formation, such as dynamics of polymerization and depolymerization, treadmilling, network formation, and the autonomous oscillatory dynamics of actomyosin. Conclusions We believe that the present model can be used to study the in vivo response of actomyosin networks to changes in key parameters of the system, such as alterations in the attachment of actin filaments to the cell cortex. Supplementary Information The online version contains supplementary material available at (10.1186/s12915-022-01279-2).
Collapse
Affiliation(s)
- Miguel Hernández-Del-Valle
- Centro de Biología Molecular Severo Ochoa, Universidad Autónoma de Madrid, Madrid, 28049, Spain.,IFIMAC, Fac. de Ciencias, Universidad Autónoma de Madrid, Madrid, 28049, Spain.,Instituto Nicolás Cabrera, Fac. de Ciencias, Universidad Autónoma de Madrid, Madrid, 28049, Spain.,Fisica de la Materia Condensada, Fac. de Ciencias, Universidad Autónoma de Madrid, Madrid, 28049, Spain
| | - Andrea Valencia-Expósito
- Centro Andaluz de Biología del Desarrollo, Universidad Pablo de Olavide/CSIC/JA, Carretera de Utrera km 1, Seville, 41013, Spain
| | - Antonio López-Izquierdo
- Centro de Biología Molecular Severo Ochoa, Universidad Autónoma de Madrid, Madrid, 28049, Spain.,IFIMAC, Fac. de Ciencias, Universidad Autónoma de Madrid, Madrid, 28049, Spain.,Instituto Nicolás Cabrera, Fac. de Ciencias, Universidad Autónoma de Madrid, Madrid, 28049, Spain.,Fisica de la Materia Condensada, Fac. de Ciencias, Universidad Autónoma de Madrid, Madrid, 28049, Spain
| | - Pau Casanova-Ferrer
- Centro de Biología Molecular Severo Ochoa, Universidad Autónoma de Madrid, Madrid, 28049, Spain.,IFIMAC, Fac. de Ciencias, Universidad Autónoma de Madrid, Madrid, 28049, Spain.,Instituto Nicolás Cabrera, Fac. de Ciencias, Universidad Autónoma de Madrid, Madrid, 28049, Spain.,Fisica de la Materia Condensada, Fac. de Ciencias, Universidad Autónoma de Madrid, Madrid, 28049, Spain
| | - Pedro Tarazona
- IFIMAC, Fac. de Ciencias, Universidad Autónoma de Madrid, Madrid, 28049, Spain.,Instituto Nicolás Cabrera, Fac. de Ciencias, Universidad Autónoma de Madrid, Madrid, 28049, Spain.,Fisica Teórica de la Materia Condensada, Fac. de Ciencias, Universidad Autónoma de Madrid, Madrid, 28049, Spain
| | - Maria D Martín-Bermudo
- Centro Andaluz de Biología del Desarrollo, Universidad Pablo de Olavide/CSIC/JA, Carretera de Utrera km 1, Seville, 41013, Spain
| | - David G Míguez
- Centro de Biología Molecular Severo Ochoa, Universidad Autónoma de Madrid, Madrid, 28049, Spain. .,IFIMAC, Fac. de Ciencias, Universidad Autónoma de Madrid, Madrid, 28049, Spain. .,Instituto Nicolás Cabrera, Fac. de Ciencias, Universidad Autónoma de Madrid, Madrid, 28049, Spain. .,Fisica de la Materia Condensada, Fac. de Ciencias, Universidad Autónoma de Madrid, Madrid, 28049, Spain.
| |
Collapse
|
18
|
Tolmachev D, Lukasheva N, Ramazanov R, Nazarychev V, Borzdun N, Volgin I, Andreeva M, Glova A, Melnikova S, Dobrovskiy A, Silber SA, Larin S, de Souza RM, Ribeiro MCC, Lyulin S, Karttunen M. Computer Simulations of Deep Eutectic Solvents: Challenges, Solutions, and Perspectives. Int J Mol Sci 2022; 23:645. [PMID: 35054840 PMCID: PMC8775846 DOI: 10.3390/ijms23020645] [Citation(s) in RCA: 30] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2021] [Revised: 01/02/2022] [Accepted: 01/04/2022] [Indexed: 12/13/2022] Open
Abstract
Deep eutectic solvents (DESs) are one of the most rapidly evolving types of solvents, appearing in a broad range of applications, such as nanotechnology, electrochemistry, biomass transformation, pharmaceuticals, membrane technology, biocomposite development, modern 3D-printing, and many others. The range of their applicability continues to expand, which demands the development of new DESs with improved properties. To do so requires an understanding of the fundamental relationship between the structure and properties of DESs. Computer simulation and machine learning techniques provide a fruitful approach as they can predict and reveal physical mechanisms and readily be linked to experiments. This review is devoted to the computational research of DESs and describes technical features of DES simulations and the corresponding perspectives on various DES applications. The aim is to demonstrate the current frontiers of computational research of DESs and discuss future perspectives.
Collapse
Affiliation(s)
- Dmitry Tolmachev
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
| | - Natalia Lukasheva
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
| | - Ruslan Ramazanov
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
| | - Victor Nazarychev
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
| | - Natalia Borzdun
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
| | - Igor Volgin
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
| | - Maria Andreeva
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
| | - Artyom Glova
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
| | - Sofia Melnikova
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
| | - Alexey Dobrovskiy
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
| | - Steven A. Silber
- Department of Physics and Astronomy, The University of Western Ontario, 1151 Richmond Street, London, ON N6A 5B7, Canada;
- The Centre of Advanced Materials and Biomaterials Research, The University of Western Ontario, 1151 Richmond Street, London, ON N6A 5B7, Canada
| | - Sergey Larin
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
| | - Rafael Maglia de Souza
- Departamento de Química Fundamental, Instituto de Química, Universidade de São Paulo, Avenida Professor Lineu Prestes 748, São Paulo 05508-070, Brazil; (R.M.d.S.); (M.C.C.R.)
| | - Mauro Carlos Costa Ribeiro
- Departamento de Química Fundamental, Instituto de Química, Universidade de São Paulo, Avenida Professor Lineu Prestes 748, São Paulo 05508-070, Brazil; (R.M.d.S.); (M.C.C.R.)
| | - Sergey Lyulin
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
| | - Mikko Karttunen
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
- Department of Physics and Astronomy, The University of Western Ontario, 1151 Richmond Street, London, ON N6A 5B7, Canada;
- The Centre of Advanced Materials and Biomaterials Research, The University of Western Ontario, 1151 Richmond Street, London, ON N6A 5B7, Canada
- Department of Chemistry, The University of Western Ontario, 1151 Richmond Street, London, ON N6A 5B7, Canada
| |
Collapse
|
19
|
Zhang F, Peng S, Xu J, Cai C, Zhang LJ. Morphological transitions of micelles induced by the block arrangements of copolymer blocks: Dissipative particle dynamics simulation. Phys Chem Chem Phys 2022; 24:10757-10764. [DOI: 10.1039/d2cp00617k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Polymer micelles with distinct morphologies and unique microphase separation microstructures can exhibit different properties and functions, holding the great promises for a range of biomedical applications. In current work, the...
Collapse
|
20
|
Speck T. Modeling of biomolecular machines in non-equilibrium steady states. J Chem Phys 2021; 155:230901. [PMID: 34937348 DOI: 10.1063/5.0070922] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
Numerical computations have become a pillar of all modern quantitative sciences. Any computation involves modeling-even if often this step is not made explicit-and any model has to neglect details while still being physically accurate. Equilibrium statistical mechanics guides both the development of models and numerical methods for dynamics obeying detailed balance. For systems driven away from thermal equilibrium, such a universal theoretical framework is missing. For a restricted class of driven systems governed by Markov dynamics and local detailed balance, stochastic thermodynamics has evolved to fill this gap and to provide fundamental constraints and guiding principles. The next step is to advance stochastic thermodynamics from simple model systems to complex systems with tens of thousands or even millions of degrees of freedom. Biomolecules operating in the presence of chemical gradients and mechanical forces are a prime example for this challenge. In this Perspective, we give an introduction to isothermal stochastic thermodynamics geared toward the systematic multiscale modeling of the conformational dynamics of biomolecular and synthetic machines, and we outline some of the open challenges.
Collapse
Affiliation(s)
- Thomas Speck
- Institut für Physik, Johannes Gutenberg-Universität Mainz, Staudingerweg 7-9, 55128 Mainz, Germany
| |
Collapse
|
21
|
Nasedkin A, Ermilova I, Swenson J. Atomistic molecular dynamics simulations of tubulin heterodimers explain the motion of a microtubule. EUROPEAN BIOPHYSICS JOURNAL : EBJ 2021; 50:927-940. [PMID: 34215900 PMCID: PMC8448678 DOI: 10.1007/s00249-021-01553-1] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/18/2021] [Revised: 05/24/2021] [Accepted: 06/07/2021] [Indexed: 06/13/2023]
Abstract
Microtubules are essential parts of the cytoskeleton that are built by polymerization of tubulin heterodimers into a hollow tube. Regardless that their structures and functions have been comprehensively investigated in a modern soft matter, it is unclear how properties of tubulin heterodimer influence and promote the self-assembly. A detailed knowledge of such structural mechanisms would be helpful in drug design against neurodegenerative diseases, cancer, diabetes etc. In this work atomistic molecular dynamics simulations were used to investigate the fundamental dynamics of tubulin heterodimers in a sheet and a short microtubule utilizing well-equilibrated structures. The breathing motions of the tubulin heterodimers during assembly show that the movement at the lateral interface between heterodimers (wobbling) dominates in the lattice. The simulations of the protofilament curvature agrees well with recently published experimental data, showing curved protofilaments at polymerization of the microtubule plus end. The tubulin heterodimers exposed at the microtubule minus end were less curved and displayed altered interactions at the site of sheet closure around the outmost heterodimers, which may slow heterodimer binding and polymerization, providing a potential explanation for the limited dynamics observed at the minus end.
Collapse
Affiliation(s)
- Alexandr Nasedkin
- Department of Physics, Chalmers University of Technology, SE 41296 Göteborg, Sweden
| | - Inna Ermilova
- Department of Physics, Chalmers University of Technology, SE 41296 Göteborg, Sweden
| | - Jan Swenson
- Department of Physics, Chalmers University of Technology, SE 41296 Göteborg, Sweden
| |
Collapse
|
22
|
Sengar A, Ouldridge TE, Henrich O, Rovigatti L, Šulc P. A Primer on the oxDNA Model of DNA: When to Use it, How to Simulate it and How to Interpret the Results. Front Mol Biosci 2021; 8:693710. [PMID: 34235181 PMCID: PMC8256390 DOI: 10.3389/fmolb.2021.693710] [Citation(s) in RCA: 27] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2021] [Accepted: 05/27/2021] [Indexed: 11/13/2022] Open
Abstract
The oxDNA model of Deoxyribonucleic acid has been applied widely to systems in biology, biophysics and nanotechnology. It is currently available via two independent open source packages. Here we present a set of clearly documented exemplar simulations that simultaneously provide both an introduction to simulating the model, and a review of the model's fundamental properties. We outline how simulation results can be interpreted in terms of-and feed into our understanding of-less detailed models that operate at larger length scales, and provide guidance on whether simulating a system with oxDNA is worthwhile.
Collapse
Affiliation(s)
- A. Sengar
- Centre for Synthetic Biology, Department of Bioengineering, Imperial College London, London, United Kingdom
| | - T. E. Ouldridge
- Centre for Synthetic Biology, Department of Bioengineering, Imperial College London, London, United Kingdom
| | - O. Henrich
- Department of Physics, SUPA, University of Strathclyde, Glasgow, United Kingdom
| | - L. Rovigatti
- Department of Physics, Sapienza University of Rome, Rome, Italy
- CNR Institute of Complex Systems, Sapienza University of Rome, Rome, Italy
| | - P. Šulc
- Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, Tempe, AZ, United States
- School of Molecular Sciences, Arizona State University, Tempe, AZ, United States
| |
Collapse
|
23
|
Berressem F, Scherer C, Andrienko D, Nikoubashman A. Ultra-coarse-graining of homopolymers in inhomogeneous systems. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2021; 33:254002. [PMID: 33845463 DOI: 10.1088/1361-648x/abf6e2] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/01/2020] [Accepted: 04/12/2021] [Indexed: 06/12/2023]
Abstract
We develop coarse-grained (CG) models for simulating homopolymers in inhomogeneous systems, focusing on polymer films and droplets. If the CG polymers interact solely through two-body potentials, then the films and droplets either dissolve or collapse into small aggregates, depending on whether the effective polymer-polymer interactions have been determined from reference simulations in the bulk or at infinite dilution. To address this shortcoming, we include higher order interactions either through an additional three-body potential or a local density-dependent potential (LDP). We parameterize the two- and three-body potentials via force matching, and the LDP through relative entropy minimization. While the CG models with three-body interactions fail at reproducing stable polymer films and droplets, CG simulations with an LDP are able to do so. Minor quantitative differences between the reference and the CG simulations, namely a slight broadening of interfaces accompanied by a smaller surface tension in the CG simulations, can be attributed to the deformation of polymers near the interfaces, which cannot be resolved in the CG representation, where the polymers are mapped to spherical beads.
Collapse
Affiliation(s)
- Fabian Berressem
- Institute of Physics, Johannes Gutenberg University Mainz, Staudingerweg 7, 55128 Mainz, Germany
| | - Christoph Scherer
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Denis Andrienko
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Arash Nikoubashman
- Institute of Physics, Johannes Gutenberg University Mainz, Staudingerweg 7, 55128 Mainz, Germany
| |
Collapse
|
24
|
Mayoral E, Goicochea AG. Modeling of Branched Thickening Polymers under Poiseuille Flow Gives Clues as to How to Increase a Solvent's Viscosity. J Phys Chem B 2021; 125:1692-1704. [PMID: 33544598 DOI: 10.1021/acs.jpcb.0c11087] [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
The viscosity enhancement of a solvent produced by the addition of thickening branched polymers is predicted as a function of polymer concentration, branch length and persistence length, and strength of the covalent bonding interactions. Nonequilibrium, stationary-state Poiseuille numerical simulations are performed using the dissipative particle dynamics model to obtain the viscosity of the fluid. It is found that the clustering of the polymers into aggregates increases the viscosity and that it is more strongly affected by the strength of the bonding interactions. General scaling relationships are found for the viscosity as a function of the variables studied, which are expected to be useful for the design and synthesis of new viscosifying polymers. It is argued that our results can be applied to aqueous thickeners, of importance for colloidal fluids such as paints and coatings and also for nonpolar fluids such as supercritical CO2, which is a promising nonhydraulic fracking fluid also useful in enhanced oil recovery.
Collapse
Affiliation(s)
- E Mayoral
- Instituto Nacional de Investigaciones Nucleares, Carretera México Toluca s/n, La Marquesa Ocoyoacac 52750, Estado de México, Mexico
| | - A Gama Goicochea
- Departamento de Ingeniería Química y Bioquímica, Tecnológico de Estudios Superiores de Ecatepec, Ecatepec de Morelos 55210, Estado de México, Mexico
| |
Collapse
|
25
|
Jin J, Han Y, Pak AJ, Voth GA. A new one-site coarse-grained model for water: Bottom-up many-body projected water (BUMPer). I. General theory and model. J Chem Phys 2021; 154:044104. [PMID: 33514116 PMCID: PMC7826168 DOI: 10.1063/5.0026651] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2020] [Accepted: 12/14/2020] [Indexed: 12/26/2022] Open
Abstract
Water is undoubtedly one of the most important molecules for a variety of chemical and physical systems, and constructing precise yet effective coarse-grained (CG) water models has been a high priority for computer simulations. To recapitulate important local correlations in the CG water model, explicit higher-order interactions are often included. However, the advantages of coarse-graining may then be offset by the larger computational cost in the model parameterization and simulation execution. To leverage both the computational efficiency of the CG simulation and the inclusion of higher-order interactions, we propose a new statistical mechanical theory that effectively projects many-body interactions onto pairwise basis sets. The many-body projection theory presented in this work shares similar physics from liquid state theory, providing an efficient approach to account for higher-order interactions within the reduced model. We apply this theory to project the widely used Stillinger-Weber three-body interaction onto a pairwise (two-body) interaction for water. Based on the projected interaction with the correct long-range behavior, we denote the new CG water model as the Bottom-Up Many-Body Projected Water (BUMPer) model, where the resultant CG interaction corresponds to a prior model, the iteratively force-matched model. Unlike other pairwise CG models, BUMPer provides high-fidelity recapitulation of pair correlation functions and three-body distributions, as well as N-body correlation functions. BUMPer extensively improves upon the existing bottom-up CG water models by extending the accuracy and applicability of such models while maintaining a reduced computational cost.
Collapse
Affiliation(s)
- Jaehyeok Jin
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - Yining Han
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - Alexander J. Pak
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - Gregory A. Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
26
|
Jin J, Pak AJ, Han Y, Voth GA. A new one-site coarse-grained model for water: Bottom-up many-body projected water (BUMPer). II. Temperature transferability and structural properties at low temperature. J Chem Phys 2021; 154:044105. [PMID: 33514078 PMCID: PMC7826166 DOI: 10.1063/5.0026652] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2020] [Accepted: 12/14/2020] [Indexed: 11/14/2022] Open
Abstract
A number of studies have constructed coarse-grained (CG) models of water to understand its anomalous properties. Most of these properties emerge at low temperatures, and an accurate CG model needs to be applicable to these low-temperature ranges. However, direct use of CG models parameterized from other temperatures, e.g., room temperature, encounters a problem known as transferability, as the CG potential essentially follows the form of the many-body CG free energy function. Therefore, temperature-dependent changes to CG interactions must be accounted for. The collective behavior of water at low temperature is generally a many-body process, which often motivates the use of expensive many-body terms in the CG interactions. To surmount the aforementioned problems, we apply the Bottom-Up Many-Body Projected Water (BUMPer) CG model constructed from Paper I to study the low-temperature behavior of water. We report for the first time that the embedded three-body interaction enables BUMPer, despite its pairwise form, to capture the growth of ice at the ice/water interface with corroborating many-body correlations during the crystal growth. Furthermore, we propose temperature transferable BUMPer models that are indirectly constructed from the free energy decomposition scheme. Changes in CG interactions and corresponding structures are faithfully recapitulated by this framework. We further extend BUMPer to examine its ability to predict the structure, density, and diffusion anomalies by employing an alternative analysis based on structural correlations and pairwise potential forms to predict such anomalies. The presented analysis highlights the existence of these anomalies in the low-temperature regime and overcomes potential transferability problems.
Collapse
Affiliation(s)
- Jaehyeok Jin
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - Alexander J. Pak
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - Yining Han
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - Gregory A. Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
27
|
DeLyser M, Noid WG. Bottom-up coarse-grained models for external fields and interfaces. J Chem Phys 2020; 153:224103. [PMID: 33317310 DOI: 10.1063/5.0030103] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023] Open
Abstract
Bottom-up coarse-grained (CG) models accurately describe the structure of homogeneous systems but sometimes provide limited transferability and a poor description of thermodynamic properties. Consequently, inhomogeneous systems present a severe challenge for bottom-up models. In this work, we examine bottom-up CG models for interfaces and inhomogeneous systems. We first analyze the effect of external fields upon the many-body potential of mean force. We also demonstrate that the multiscale CG (MS-CG) variational principle for modeling the external field corresponds to a generalization of the first Yvon-Born-Green equation. This provides an important connection with liquid state theory, as well as physical insight into the structure of interfaces and the resulting MS-CG models. We then develop and assess MS-CG models for a film of liquid methanol that is adsorbed on an attractive wall and in coexistence with its vapor phase. While pair-additive potentials provide unsatisfactory accuracy and transferability, the inclusion of local-density (LD) potentials dramatically improves the accuracy and transferability of the MS-CG model. The MS-CG model with LD potentials quite accurately describes the wall-liquid interface, the bulk liquid density, and the liquid-vapor interface while simultaneously providing a much improved description of the vapor phase. This model also provides an excellent description of the pair structure and pressure-density equation of state for the bulk liquid. Thus, LD potentials hold considerable promise for transferable bottom-up models that accurately describe the structure and thermodynamic properties of both bulk and interfacial systems.
Collapse
Affiliation(s)
- Michael DeLyser
- Department of Chemistry, Penn State University, University Park, Pennsylvania 16802, USA
| | - W G Noid
- Department of Chemistry, Penn State University, University Park, Pennsylvania 16802, USA
| |
Collapse
|
28
|
Rudzinski JF, Bereau T. Coarse-grained conformational surface hopping: Methodology and transferability. J Chem Phys 2020; 153:214110. [DOI: 10.1063/5.0031249] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
| | - Tristan Bereau
- Max Planck Institute for Polymer Research, 55128 Mainz, Germany
- Van ’t Hoff Institute for Molecular Sciences and Informatics Institute, University of Amsterdam, Amsterdam 1098 XH, The Netherlands
| |
Collapse
|
29
|
Bunker A, Róg T. Mechanistic Understanding From Molecular Dynamics Simulation in Pharmaceutical Research 1: Drug Delivery. Front Mol Biosci 2020; 7:604770. [PMID: 33330633 PMCID: PMC7732618 DOI: 10.3389/fmolb.2020.604770] [Citation(s) in RCA: 41] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2020] [Accepted: 11/02/2020] [Indexed: 12/12/2022] Open
Abstract
In this review, we outline the growing role that molecular dynamics simulation is able to play as a design tool in drug delivery. We cover both the pharmaceutical and computational backgrounds, in a pedagogical fashion, as this review is designed to be equally accessible to pharmaceutical researchers interested in what this new computational tool is capable of and experts in molecular modeling who wish to pursue pharmaceutical applications as a context for their research. The field has become too broad for us to concisely describe all work that has been carried out; many comprehensive reviews on subtopics of this area are cited. We discuss the insight molecular dynamics modeling has provided in dissolution and solubility, however, the majority of the discussion is focused on nanomedicine: the development of nanoscale drug delivery vehicles. Here we focus on three areas where molecular dynamics modeling has had a particularly strong impact: (1) behavior in the bloodstream and protective polymer corona, (2) Drug loading and controlled release, and (3) Nanoparticle interaction with both model and biological membranes. We conclude with some thoughts on the role that molecular dynamics simulation can grow to play in the development of new drug delivery systems.
Collapse
Affiliation(s)
- Alex Bunker
- Division of Pharmaceutical Biosciences, Drug Research Program, Faculty of Pharmacy, University of Helsinki, Helsinki, Finland
| | - Tomasz Róg
- Department of Physics, University of Helsinki, Helsinki, Finland
| |
Collapse
|
30
|
Li W, Burkhart C, Polińska P, Harmandaris V, Doxastakis M. Backmapping coarse-grained macromolecules: An efficient and versatile machine learning approach. J Chem Phys 2020; 153:041101. [DOI: 10.1063/5.0012320] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Wei Li
- Department of Chemical and Biomolecular Engineering, University of Tennessee, Knoxville, Tennessee 37996, USA
| | - Craig Burkhart
- The Goodyear Tire and Rubber Company, Akron, Ohio 44305, USA
| | - Patrycja Polińska
- Goodyear Innovation Center Luxembourg, Avenue Gordon Smith, L-7750 Colmar-Berg, Luxembourg
| | - Vagelis Harmandaris
- Department of Applied Mathematics, University of Crete, and IACM FORTH, GR-71110 Heraklion, Greece
- Computation-Based Science and Technology Research Center, The Cyprus Institute, Nicosia 2121, Cyprus
| | - Manolis Doxastakis
- Department of Chemical and Biomolecular Engineering, University of Tennessee, Knoxville, Tennessee 37996, USA
| |
Collapse
|
31
|
Haley NEC, Ouldridge TE, Mullor Ruiz I, Geraldini A, Louis AA, Bath J, Turberfield AJ. Design of hidden thermodynamic driving for non-equilibrium systems via mismatch elimination during DNA strand displacement. Nat Commun 2020; 11:2562. [PMID: 32444600 PMCID: PMC7244503 DOI: 10.1038/s41467-020-16353-y] [Citation(s) in RCA: 57] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2019] [Accepted: 04/15/2020] [Indexed: 11/09/2022] Open
Abstract
Recent years have seen great advances in the development of synthetic self-assembling molecular systems. Designing out-of-equilibrium architectures, however, requires a more subtle control over the thermodynamics and kinetics of reactions. We propose a mechanism for enhancing the thermodynamic drive of DNA strand-displacement reactions whilst barely perturbing forward reaction rates: the introduction of mismatches within the initial duplex. Through a combination of experiment and simulation, we demonstrate that displacement rates are strongly sensitive to mismatch location and can be tuned by rational design. By placing mismatches away from duplex ends, the thermodynamic drive for a strand-displacement reaction can be varied without significantly affecting the forward reaction rate. This hidden thermodynamic driving motif is ideal for the engineering of non-equilibrium systems that rely on catalytic control and must be robust to leak reactions.
Collapse
Affiliation(s)
- Natalie E C Haley
- Clarendon Laboratory, Department of Physics, University of Oxford, Parks Road, Oxford, OX1 3PU, UK
| | - Thomas E Ouldridge
- Imperial College Centre for Synthetic Biology and Department of Bioengineering, Prince Consort Road, Imperial College London, London, SW7 2AZ, UK.
| | - Ismael Mullor Ruiz
- Imperial College Centre for Synthetic Biology and Department of Bioengineering, Prince Consort Road, Imperial College London, London, SW7 2AZ, UK
| | - Alessandro Geraldini
- Rudolf Peierls Centre for Theoretical Physics, Department of Physics, University of Oxford, Keble Road, Oxford, OX1 3NP, UK
| | - Ard A Louis
- Rudolf Peierls Centre for Theoretical Physics, Department of Physics, University of Oxford, Keble Road, Oxford, OX1 3NP, UK
| | - Jonathan Bath
- Clarendon Laboratory, Department of Physics, University of Oxford, Parks Road, Oxford, OX1 3PU, UK
| | - Andrew J Turberfield
- Clarendon Laboratory, Department of Physics, University of Oxford, Parks Road, Oxford, OX1 3PU, UK.
| |
Collapse
|
32
|
Colla T, Bakhshandeh A, Levin Y. Osmotic stress and pore nucleation in charged biological nanoshells and capsids. SOFT MATTER 2020; 16:2390-2405. [PMID: 32067009 DOI: 10.1039/c9sm02532d] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
A model system is proposed to investigate the chemical equilibrium and mechanical stability of biological spherical-like nanoshells in contact with an aqueous solution with added dissociated electrolyte of a given concentration. The ionic chemical equilibrium across the permeable shell is investigated in the framework of an accurate Density Functional Theory (DFT) that incorporates electrostatic and hardcore correlations beyond the traditional mean-field (e.g., Poisson-Boltzmann) limit. The accuracy of the theory is tested by a direct comparison with Monte Carlo (MC) simulations. A simple analytical expression is then deduced which clearly highlights the entropic, electrostatic, and self-energy contributions to the osmotic stress over the shell in terms of the calculated ionic profiles. By invoking a continuum mean-field elastic approach to account for the shell surface stress upon osmotic stretching, the mechanical equilibrium properties of the shell under a wide variety of ionic strengths and surface charges are investigated. The model is further coupled to a continuum mechanical approach similar in structure to a Classical Nucleation Theory (CNT) to address the question of mechanical stability of the shells against a pore nucleation. This allows us to construct a phase diagram which delimits the mechanical stability of capsids for different ionic strengths and shell surface charges.
Collapse
Affiliation(s)
- Thiago Colla
- Instituto de Física, Universidade Federal de Ouro Preto, CEP 35400-000, Ouro Preto, MG, Brazil.
| | - Amin Bakhshandeh
- Programa de Pós-Graduação em Física, Instituto de Física e Matemática, Universidade Federal de Pelotas, Caixa Postal 354, CEP 96010-900 Pelotas, RS, Brazil.
| | - Yan Levin
- Instituto de Física, Universidade Federal do Rio Grande do Sul, Caixa Postal 15051, CEP 91501-970, Porto Alegre, RS, Brazil.
| |
Collapse
|
33
|
de Souza RM, Ratochinski RH, Karttunen M, Dias LG. Self-Assembly of Phosphocholine Derivatives Using the ELBA Coarse-Grained Model: Micelles, Bicelles, and Reverse Micelles. J Chem Inf Model 2020; 60:522-536. [PMID: 31714768 DOI: 10.1021/acs.jcim.9b00790] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
The ELBA coarse-grained force field was originally developed for lipids, and its water model is described as a single-site Lennard-Jones particle with electrostatics modeled by an embedded point-dipole, while other molecules in this force field have a three (or four)-to-one mapping scheme. Here, ELBA was applied to investigate the self-assembly processes of dodecyl-phosphocholine (DPC) micelle, 1,2-dipalmitoyl-sn-glycero-3-phosphocholine/1,2-dihexaoyl-sn-glycero-3-phosphocholine (DPPC/DHPC) bicelles, and DPPC/cyclohexane/water reverse micelles through coarse-grained molecular dynamics (MD) simulations. New parameters were obtained using a simplex algorithm-based calibration procedure to determine the Lennard-Jones parameters for cyclohexane, dodecane, and cyclohexane-dodecane cross-interactions. Density, self-diffusion coefficient, surface tension, and mixture excess volume were found to be in fair agreement with experimental data. These new parameters were used in the simulations, and the obtained structures were analyzed for shape, size, volume, and surface area. Except for the shape of DPC micelles, all other properties match well with available experimental data and all-atom simulations. Remarkably, in agreement with experiments the rodlike shape of the DPPC reverse micelle is well described by ELBA, while all-atom data in the literature predicts a disclike shape. To further check the consistency of the force field in reproducing the correct shapes of reverse micelles, additional simulations were performed doubling the system size. Two distinct reverse micelles were obtained both presenting the rodlike shape and correct aggregation number.
Collapse
Affiliation(s)
- R M de Souza
- Department of Chemistry , The University of Western Ontario , London , Ontario , Canada N6A 3K7.,Departamento de Química, FFCLRP , Universidade de São Paulo , Avenida Bandeirantes 3900 , 14040-901 Ribeirão Preto , SP , Brazil.,The Center for Advanced Materials and Biomaterials Research , The University of Western Ontario , London , Ontario , Canada N6K 3K7
| | - R H Ratochinski
- Departamento de Química, FFCLRP , Universidade de São Paulo , Avenida Bandeirantes 3900 , 14040-901 Ribeirão Preto , SP , Brazil
| | - Mikko Karttunen
- Department of Chemistry , The University of Western Ontario , London , Ontario , Canada N6A 3K7.,The Center for Advanced Materials and Biomaterials Research , The University of Western Ontario , London , Ontario , Canada N6K 3K7.,Department of Applied Mathematics , The University of Western Ontario , London , Ontario , Canada N6A 5B7
| | - L G Dias
- Departamento de Química, FFCLRP , Universidade de São Paulo , Avenida Bandeirantes 3900 , 14040-901 Ribeirão Preto , SP , Brazil
| |
Collapse
|
34
|
Honda K. A coarse-grained potential for liquid water by spherically projected hydrogen-bonding interactions. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2019.112246] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
35
|
Sicard F, Destainville N, Rousseau P, Tardin C, Manghi M. Dynamical control of denaturation bubble nucleation in supercoiled DNA minicircles. Phys Rev E 2020; 101:012403. [PMID: 32069623 DOI: 10.1103/physreve.101.012403] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2019] [Indexed: 06/10/2023]
Abstract
We examine the behavior of supercoiled DNA minicircles containing between 200 and 400 base-pairs, also named microDNA, in which supercoiling favors thermally assisted DNA denaturation bubbles of nanometer size and controls their lifetime. Mesoscopic modeling and accelerated dynamics simulations allow us to overcome the limitations of atomistic simulations encountered in such systems, and offer detailed insight into the thermodynamic and dynamical properties associated with the nucleation and closure mechanisms of long-lived thermally assisted denaturation bubbles which do not stem from bending- or torque-driven stress. Suitable tuning of the degree of supercoiling and size of specifically designed microDNA is observed to lead to the control of opening characteristic times in the millisecond range, and closure characteristic times ranging over well distinct timescales, from microseconds to several minutes. We discuss how our results can be seen as a dynamical bandwidth which might enhance selectivity for specific DNA binding proteins.
Collapse
Affiliation(s)
- François Sicard
- Department of Chemistry, King's College London, SE1 1DB London, United Kingdom
| | - Nicolas Destainville
- Laboratoire de Physique Théorique, IRSAMC, Université de Toulouse, CNRS, UPS, France
| | - Philippe Rousseau
- Laboratoire de Microbiologie et Génetique Moléculaires, Centre de Biologie Intégrative (CBI), Université de Toulouse, CNRS, UPS, France
| | - Catherine Tardin
- Institut de Pharmacologie et Biologie Structurale, Université de Toulouse, CNRS, UPS, France
| | - Manoel Manghi
- Laboratoire de Physique Théorique, IRSAMC, Université de Toulouse, CNRS, UPS, France
| |
Collapse
|
36
|
Chen T, Manz TA. A collection of forcefield precursors for metal-organic frameworks. RSC Adv 2019; 9:36492-36507. [PMID: 35539031 PMCID: PMC9075174 DOI: 10.1039/c9ra07327b] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2019] [Accepted: 10/25/2019] [Indexed: 12/14/2022] Open
Abstract
A host of important performance properties for metal-organic frameworks (MOFs) and other complex materials can be calculated by modeling statistical ensembles. The principle challenge is to develop accurate and computationally efficient interaction models for these simulations. Two major approaches are (i) ab initio molecular dynamics in which the interaction model is provided by an exchange-correlation theory (e.g., DFT + dispersion functional) and (ii) molecular mechanics in which the interaction model is a parameterized classical force field. The first approach requires further development to improve computational speed. The second approach requires further development to automate accurate forcefield parameterization. Because of the extreme chemical diversity across thousands of MOF structures, this problem is still mostly unsolved today. For example, here we show structures in the 2014 CoRE MOF database contain more than 8 thousand different atom types based on first and second neighbors. Our results showed that atom types based on both first and second neighbors adequately capture the chemical environment, but atom types based on only first neighbors do not. For 3056 MOFs, we used density functional theory (DFT) followed by DDEC6 atomic population analysis to extract a host of important forcefield precursors: partial atomic charges; atom-in-material (AIM) C6, C8, and C10 dispersion coefficients; AIM dipole and quadrupole moments; various AIM polarizabilities; quantum Drude oscillator parameters; AIM electron cloud parameters; etc. Electrostatic parameters were validated through comparisons to the DFT-computed electrostatic potential. These forcefield precursors should find widespread applications to developing MOF force fields.
Collapse
Affiliation(s)
- Taoyi Chen
- Department of Chemical & Materials Engineering, New Mexico State University Las Cruces New Mexico 88003-8001 USA
| | - Thomas A Manz
- Department of Chemical & Materials Engineering, New Mexico State University Las Cruces New Mexico 88003-8001 USA
| |
Collapse
|
37
|
Marolf DM, Jones MR. Measurement Challenges in Dynamic and Nonequilibrium Nanoscale Systems. Anal Chem 2019; 91:13324-13336. [DOI: 10.1021/acs.analchem.9b02702] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Affiliation(s)
- David M. Marolf
- Department of Chemistry, Rice University, Houston, Texas 77005, United States
| | - Matthew R. Jones
- Department of Chemistry, Rice University, Houston, Texas 77005, United States
- Department of Materials Science and Nanoengineering, Rice University, Houston, Texas 77005, United States
| |
Collapse
|
38
|
Jin J, Pak AJ, Voth GA. Understanding Missing Entropy in Coarse-Grained Systems: Addressing Issues of Representability and Transferability. J Phys Chem Lett 2019; 10:4549-4557. [PMID: 31319036 PMCID: PMC6782054 DOI: 10.1021/acs.jpclett.9b01228] [Citation(s) in RCA: 44] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/10/2023]
Abstract
Coarse-grained (CG) models facilitate efficient simulation of complex systems by integrating out the atomic, or fine-grained (FG), degrees of freedom. Systematically derived CG models from FG simulations often attempt to approximate the CG potential of mean force (PMF), an inherently multidimensional and many-body quantity, using additive pairwise contributions. However, they currently lack fundamental principles that enable their extensible use across different thermodynamic state points, i.e., transferability. In this work, we investigate the explicit energy-entropy decomposition of the CG PMF as a means to construct transferable CG models. In particular, despite its high-dimensional nature, we find for liquid systems that the entropic component to the CG PMF can similarly be represented using additive pairwise contributions, which we show is highly coupled to the CG configurational entropy. This approach formally connects the missing entropy that is lost due to the CG representation, i.e., translational, rotational, and vibrational modes associated with the missing degrees of freedom, to the CG entropy. By design, the present framework imparts transferable CG interactions across different temperatures due to the explicit definition of an additive entropic contribution. Furthermore, we demonstrate that transferability across composition state points, such as between bulk liquids and their mixtures, is also achieved by designing combining rules to approximate cross-interactions from bulk CG PMFs. Using the predicted CG model for liquid mixtures, structural correlations of the fitted CG model were found to corroborate a high-fidelity combining rule. Our findings elucidate the physical nature and compact representation of CG entropy and suggest a new approach for overcoming the transferability problem. We expect that this approach will further extend the current view of CG modeling into predictive multiscale modeling.
Collapse
|
39
|
Meneses-Juárez E, Márquez-Beltrán C, González-Melchor M. Influence of pH on the formation of a polyelectrolyte complex by dissipative particle dynamics simulation: From an extended to a compact shape. Phys Rev E 2019; 100:012505. [PMID: 31499774 DOI: 10.1103/physreve.100.012505] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2019] [Indexed: 06/10/2023]
Abstract
This work aims to investigate the influence of pH on the mechanism of assembly of macromolecules. We studied the effect of the pH on the interaction of two polyelectrolytes of opposed charge, having the same size, by means of dissipative particle dynamics method. The system consisted of a strong cationic and a weak anionic polyelectrolyte in an aqueous solution containing monovalent counterions. The analysis was made by varying the pH of the solution, which modifies the charge fraction of the weak anionic polyelectrolyte with a dissociation acid constant pKa of 5.5, while the polycation is fully charged in all the pH range used, characteristic of a strong polyelectrolyte. In order to describe the influence of pH on the complexation process, we have analyzed the pair radial distribution functions polyanion-counterion, polycation-counterion, and polyanion-polycation. The complex conformation was studied by means of the radius of gyration and the end-to-end distance of both chains as the pH varied from 1 to 14. A relevant finding obtained here was the relationship between the radial distribution functions and the counterion release from the polyelectrolytes, which leads to a reduction in the size of the complex when pH increased. Surprisingly, a transition from an extended to a compact polyelectrolyte complex was obtained when the pH reached the dissociation acid constant pKa of the weak polyelectrolyte. This systematic study can help to understand a large number of more realistic problems in biological systems such as protein complex, chromatin phase transition, or in complex systems applied in biomedical science.
Collapse
Affiliation(s)
- Efrain Meneses-Juárez
- Instituto de Física "Luis Rivera Terrazas," Benemérita Universidad Autónoma de Puebla, Apartado Postal J-48, Puebla 72570, Mexico
- Facultad de Ciencias Básicas, Ingeniería y Tecnología, Universidad Autónoma de Tlaxcala, Calzada Apizaquito S/N, Apizaco, Tlaxcala 90300, Mexico
| | - César Márquez-Beltrán
- Instituto de Física "Luis Rivera Terrazas," Benemérita Universidad Autónoma de Puebla, Apartado Postal J-48, Puebla 72570, Mexico
| | - Minerva González-Melchor
- Instituto de Física "Luis Rivera Terrazas," Benemérita Universidad Autónoma de Puebla, Apartado Postal J-48, Puebla 72570, Mexico
| |
Collapse
|
40
|
Rosenberger D, van der Vegt NFA. Relative entropy indicates an ideal concentration for structure-based coarse graining of binary mixtures. Phys Rev E 2019; 99:053308. [PMID: 31212527 DOI: 10.1103/physreve.99.053308] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2019] [Indexed: 06/09/2023]
Abstract
Many methodological approaches have been proposed to improve systematic or bottom-up coarse-graining techniques to enhance the representability and transferability of the derived interaction potentials. Transferability describes the ability of a coarse-grained (CG) model to be predictive, i.e., to describe a system at state points different from those chosen for parametrization. Whereas the representability characterizes the accuracy of a CG model to reproduce target properties of the underlying reference or fine-grained model at a given state point. In this article, we shift the focus away from methodological aspects and rather raise the question whether we can overcome the disadvantages of a given method in terms of representability and transferability by systematically selecting the state point at which the CG model gets parametrized. We answer this question by applying the inverse Monte Carlo (IMC) approach-a structure-based coarse-graining method-to derive effective interactions for binary mixtures of simple Lennard-Jones (LJ) particles, which are different in size. For such simple systems we indeed can identify a concentration where the derived potentials show the best performance in terms of structural representability and transferability. This specific concentration is identified by computing the relative entropy which quantifies the information loss between different IMC models and the reference LJ model at varying mixture compositions. Further, we show that an IMC model for mixtures of n-hexane and n-perfluorohexane shows the same trend in transferability as the IMC models for the LJ system. All derived models are more transferable in the direction of increasing concentration of the larger-sized compound.
Collapse
Affiliation(s)
- David Rosenberger
- Eduard Zintl Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, Darmstadt, 64287, Germany
| | - Nico F A van der Vegt
- Eduard Zintl Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, Darmstadt, 64287, Germany
| |
Collapse
|
41
|
Leyva Y, Martin O, Perez N, Suarez-Lezcano J, Fundora-Pozo M. The optimal size of protocells from simple entropic considerations. EUROPEAN BIOPHYSICS JOURNAL: EBJ 2019; 48:277-283. [PMID: 30899982 DOI: 10.1007/s00249-019-01359-2] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/16/2018] [Revised: 12/11/2018] [Accepted: 03/12/2019] [Indexed: 01/29/2023]
Abstract
Potential constraints on protocell size are developed from simple entropic considerations. To do that, two new different indexes as measures of their structural and dynamic order were developed and applied to an elemental model of the heterotrophic protocell. According to our results, cell size should be a key factor determining the potential of these primitive systems to evolve and consequently to support life. Our analyses also suggest that the size of the optimal vesicles could be constrained to have radii in the interval [Formula: see text], where the two extreme limits [Formula: see text] and [Formula: see text] represent the states of maximum structural order (largest accumulation of substrate inside the vesicle) and the maximum flux of entropy production, respectively. According to the above criteria, the size of the optimum vesicles falls, approximately, in the same spatial range estimated for biological living cells assuming plausible values for the second-order rate constant involved in the catabolic process. Furthermore, the existence of very small vesicles could be seriously affected by the limited efficiency, far from the theoretical limits, with which these catabolic processes may proceed in a prebiotic system.
Collapse
Affiliation(s)
- Yoelsy Leyva
- Departamento de Física, Facultad de Ciencias, Universidad de Tarapacá, Casilla 7-D, Arica, Chile.
| | - Osmel Martin
- Laboratorio de Ciencia Planetaria, Universidad Central "Marta Abreu" de las Villas, Santa Clara, Cuba
| | - Noel Perez
- Laboratorio de Ciencia Planetaria, Universidad Central "Marta Abreu" de las Villas, Santa Clara, Cuba
| | - José Suarez-Lezcano
- Escuela de Enfermería, Pontificia Universidad Católica del Ecuador Sede Esmeraldas (PUCESE), Esmeraldas, Ecuador
| | - Manuel Fundora-Pozo
- Laboratorio de Ciencia Planetaria, Universidad Central "Marta Abreu" de las Villas, Santa Clara, Cuba
| |
Collapse
|
42
|
Goodarzi F, Zendehboudi S. Effects of Salt and Surfactant on Interfacial Characteristics of Water/Oil Systems: Molecular Dynamic Simulations and Dissipative Particle Dynamics. Ind Eng Chem Res 2019. [DOI: 10.1021/acs.iecr.9b00504] [Citation(s) in RCA: 45] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Affiliation(s)
- Fatemeh Goodarzi
- Faculty of Engineering and Applied Science, Memorial University, St. John’s, NL A1B 3X5, Canada
| | - Sohrab Zendehboudi
- Faculty of Engineering and Applied Science, Memorial University, St. John’s, NL A1B 3X5, Canada
| |
Collapse
|
43
|
Lynn H, Thachuk M. Equations of motion for position-dependent coarse-grain mappings obtained with Mori-Zwanzig theory. J Chem Phys 2019; 150:024108. [DOI: 10.1063/1.5058061] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Hudson Lynn
- Department of Chemistry, University of British Columbia, Vancouver, British Columbia V6T 1Z1, Canada
| | - Mark Thachuk
- Department of Chemistry, University of British Columbia, Vancouver, British Columbia V6T 1Z1, Canada
| |
Collapse
|
44
|
Lebold KM, Noid WG. Systematic study of temperature and density variations in effective potentials for coarse-grained models of molecular liquids. J Chem Phys 2019; 150:014104. [DOI: 10.1063/1.5050509] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023] Open
Affiliation(s)
- Kathryn M. Lebold
- Department of Chemistry, Penn State University, University Park, Pennsylvania 16802, USA
| | - W. G. Noid
- Department of Chemistry, Penn State University, University Park, Pennsylvania 16802, USA
| |
Collapse
|
45
|
Bereau T, Rudzinski JF. Accurate Structure-Based Coarse Graining Leads to Consistent Barrier-Crossing Dynamics. PHYSICAL REVIEW LETTERS 2018; 121:256002. [PMID: 30608819 DOI: 10.1103/physrevlett.121.256002] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/17/2018] [Revised: 11/09/2018] [Indexed: 06/09/2023]
Abstract
Structure-based coarse graining of molecular systems offers a systematic route to reproduce the many-body potential of mean force. Unfortunately, common strategies are inherently limited by the molecular mechanics force field employed. Here, we extend the concept of multisurface dynamics, initially developed to describe electronic transitions in chemical reactions, to accurately sample the conformational ensemble of a classical system in equilibrium. In analogy to describing different electronic configurations, a surface-hopping scheme couples distinct conformational basins beyond the additivity of the Hamiltonian. The incorporation of more surfaces leads systematically toward improved cross-correlations. The resulting models naturally achieve consistent long-time dynamics for systems governed by barrier-crossing events.
Collapse
Affiliation(s)
- Tristan Bereau
- Max Planck Institute for Polymer Research, 55128 Mainz, Germany
| | | |
Collapse
|
46
|
Pannuzzo M, McDargh ZA, Deserno M. The role of scaffold reshaping and disassembly in dynamin driven membrane fission. eLife 2018; 7:39441. [PMID: 30561335 PMCID: PMC6355196 DOI: 10.7554/elife.39441] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2018] [Accepted: 12/13/2018] [Indexed: 12/11/2022] Open
Abstract
The large GTPase dynamin catalyzes membrane fission in eukaryotic cells, but despite three decades of experimental work, competing and partially conflicting models persist regarding some of its most basic actions. Here we investigate the mechanical and functional consequences of dynamin scaffold shape changes and disassembly with the help of a geometrically and elastically realistic simulation model of helical dynamin-membrane complexes. Beyond changes of radius and pitch, we emphasize the crucial role of a third functional motion: an effective rotation of the filament around its longitudinal axis, which reflects alternate tilting of dynamin’s PH binding domains and creates a membrane torque. We also show that helix elongation impedes fission, hemifission is reached via a small transient pore, and coat disassembly assists fission. Our results have several testable structural consequences and help to reconcile mutual conflicting aspects between the two main present models of dynamin fission—the two-stage and the constrictase model. When cells take up material from their surroundings, they must first transport this cargo across their outer membrane, a flexible sheet of tightly organized fat molecules that act as a barrier to the environment. Cells can achieve this by letting their membrane surround the object, pulling it inwards until it is contained in a pouch that bulges into the cell. This bag is then corded up so it splits off from the outer membrane. The ‘cord’ is a protein called dynamin, which is thought to form a tight spiral around the bag’s neck, closing it over and pinching it away. The structure of dynamin is fairly well known, and yet several theories compete to explain how it may snap the bag off the outer membrane. Here, Pannuzzo et al. have created a computer simulation that faithfully replicates the geometry and the elasticity of the membrane and of dynamin, and used it to test different ways the protein could work. The first test featured simple constriction, where the dynamin spiral contracts around the membrane to pinch it; this only separated the bag from the membrane after implausibly tight constriction. The second test added elongation, with the spiral lengthening as well as reducing its diameter, but this further reduced the ability for the protein to snap off the membrane. The final test combined constriction and rotation, whereby dynamin ‘twirls’ as it presses on the neck of the bag: this succeeded in efficiently severing the membrane once the dynamin spiral disassembled. Indeed, the simulations suggested that dynamin might start to dismantle while it constricts, without compromising its role. In fact, getting rid of excess length as the protein contracts helps to dissolve any remnants of a membrane connection. Defects in dynamin are associated with conditions such as centronuclear myopathy and Charcot‐Marie‐Tooth peripheral neuropathy. Recent research also indicates that the protein is involved in a much wider range of neurological disorders that include Alzheimer's, Parkinson's, Huntington's, and amyotrophic lateral sclerosis. The models created by Pannuzzo et al. are useful tools to understand how dynamin and similar proteins work and sometimes fail.
Collapse
Affiliation(s)
- Martina Pannuzzo
- Department of Physics, Carnegie Mellon University, Pittsburgh, United States
| | - Zachary A McDargh
- Department of Physics, Carnegie Mellon University, Pittsburgh, United States
| | - Markus Deserno
- Department of Physics, Carnegie Mellon University, Pittsburgh, United States
| |
Collapse
|
47
|
Diggins P, Liu C, Deserno M, Potestio R. Optimal Coarse-Grained Site Selection in Elastic Network Models of Biomolecules. J Chem Theory Comput 2018; 15:648-664. [PMID: 30514085 PMCID: PMC6391041 DOI: 10.1021/acs.jctc.8b00654] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Elastic network models, simple structure-based representations of biomolecules where atoms interact via short-range harmonic potentials, provide great insight into a molecule's internal dynamics and mechanical properties at extremely low computational cost. Their efficiency and effectiveness have made them a pivotal instrument in the computer-aided study of proteins and, since a few years, also of nucleic acids. In general, the coarse-grained sites, i.e. those effective force centers onto which the all-atom structure is mapped, are constructed based on intuitive rules: a typical choice for proteins is to retain only the C α atoms of each amino acid. However, a mapping strategy relying only on the atom type and not the local properties of its embedding can be suboptimal compared to a more careful selection. Here, we present a strategy in which the subset of atoms, each of which is mapped onto a unique coarse-grained site of the model, is selected in a stochastic search aimed at optimizing a cost function. The latter is taken to be a simple measure of the consistency between the harmonic approximation of an elastic network model and the harmonic model obtained through exact integration of the discarded degrees of freedom. The method is applied to two representatives of structurally very different types of biomolecules: the protein adenylate kinase and the RNA molecule adenine riboswitch. Our analysis quantifies the substantial impact that an algorithm-driven selection of coarse-grained sites can have on a model's properties.
Collapse
Affiliation(s)
- Patrick Diggins
- Department of Physics , Carnegie Mellon University , Pittsburgh , Pennsylvania 15213 , United States
| | - Changjiang Liu
- Department of Physics , Carnegie Mellon University , Pittsburgh , Pennsylvania 15213 , United States.,Department of Biophysics , University of Michigan , Ann Arbor , Michigan 48109 , United States
| | - Markus Deserno
- Department of Physics , Carnegie Mellon University , Pittsburgh , Pennsylvania 15213 , United States
| | - Raffaello Potestio
- Physics Department , University of Trento , via Sommarive, 14 I-38123 Trento , Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications , I-38123 Trento , Italy
| |
Collapse
|
48
|
Jin J, Han Y, Voth GA. Ultra-Coarse-Grained Liquid State Models with Implicit Hydrogen Bonding. J Chem Theory Comput 2018; 14:6159-6174. [PMID: 30354110 DOI: 10.1021/acs.jctc.8b00812] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Coarse-graining (CG) methodologies have been widely used to extend the time and length scales of computer simulations by averaging over the atomistic details beneath the resolution of the CG models. Despite the efficiency of CG models, important configurational information during a given process may be lost at the CG resolution. One example of this is the topology of the hydrogen bonding network in the liquid state. When the functional group that participates in hydrogen bonding (e.g., -OH in methanol) is coarse-grained into one CG site, the effective interactions of the resultant CG model are usually derived from an averaged overall trajectory and, thus, do not take into account the hydrogen bonding interactions and topologies that are present at the all-atom resolution. In order to overcome this challenge, the present study develops new ultra-coarse-grained (UCG) models that include internal states within the CG sites that participate in hydrogen bonding, where each state represents a specific configuration such as the hydrogen bonding donor or acceptor. Internal states of the UCG beads are modeled to remain in quasi-equilibrium, and the degree of mixing is controlled by utilizing the effective local density of the UCG sites. In particular, we consider two groups of UCG models with different types of hydrogen bonding motifs: chain-like and ring-like. Using five different liquid systems that contain the same fundamental functional groups as biomolecules, we demonstrate the ability of the UCG models to reproduce the structural properties that originate from the configurations beneath the resolution of the UCG model. This proposed approach can also be applied to other liquids with such specific and directional interactions, or even to complex biomolecular systems in which hydrogen bonding is critical.
Collapse
Affiliation(s)
- Jaehyeok Jin
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics , The University of Chicago , Chicago , Illinois 60637 , United States
| | - Yining Han
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics , The University of Chicago , Chicago , Illinois 60637 , United States
| | - Gregory A Voth
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics , The University of Chicago , Chicago , Illinois 60637 , United States
| |
Collapse
|
49
|
Chakraborty M, Xu C, White AD. Encoding and selecting coarse-grain mapping operators with hierarchical graphs. J Chem Phys 2018; 149:134106. [PMID: 30292213 DOI: 10.1063/1.5040114] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Coarse-grained (CG) molecular dynamics (MD) can simulate systems inaccessible to fine-grained (FG) MD simulations. A CG simulation decreases the degrees of freedom by mapping atoms from an FG representation into agglomerate CG particles. The FG to CG mapping is not unique. Research into systematic selection of these mappings is challenging due to their combinatorial growth with respect to the number of atoms in a molecule. Here we present a method of reducing the total count of mappings by imposing molecular topology and symmetry constraints. The count reduction is illustrated by considering all mappings for nearly 50 000 molecules. The resulting number of mapping operators is still large, so we introduce a novel hierarchical graphical approach which encodes multiple CG mapping operators. The encoding method is demonstrated for methanol and a 14-mer peptide. With the test cases, we show how the encoding can be used for automated selection of reasonable CG mapping operators.
Collapse
Affiliation(s)
- Maghesree Chakraborty
- Department of Chemical Engineering, University of Rochester, Rochester, New York 14627, USA
| | - Chenliang Xu
- Department of Computer Science, University of Rochester, Rochester, New York 14627, USA
| | - Andrew D White
- Department of Chemical Engineering, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
50
|
Who is left out: exploring social boundaries in entrepreneurial ecosystems. JOURNAL OF TECHNOLOGY TRANSFER 2018. [DOI: 10.1007/s10961-018-9694-0] [Citation(s) in RCA: 57] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|