1
|
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
|
2
|
Chaimovich M, Chaimovich A. Relative Resolution: An Analysis with the Kullback-Leibler Entropy. J Chem Theory Comput 2024; 20:2074-2087. [PMID: 38416535 DOI: 10.1021/acs.jctc.3c01052] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/29/2024]
Abstract
A novel type of a multiscale approach, called Relative Resolution (RelRes), can correctly retrieve the behavior of various nonpolar liquids while speeding up molecular simulations by almost an order of magnitude. In this approach in a single system, molecules switch their resolution in terms of their relative separation, with near neighbors interacting via fine-grained potentials, yet far neighbors interacting via coarse-grained potentials; notably, these two potentials are analytically parametrized by a multipole approximation. Our current work focuses on analyzing RelRes by relating it with the Kullback-Leibler (KL) entropy, which is a useful metric for multiscale errors. In particular, we thoroughly examine the exact and approximate versions of this informatic measure for several alkane systems. By analyzing its dependency on the system size, we devise a formula for predicting the exact KL entropy of an "infinite" system via the computation of the approximate KL entropy of an "infinitesimal" system. Demonstrating that the KL entropy can holistically capture many multiscale errors, we settle bounds for the KL entropy that ensure a sufficient representation of the structural and thermal behavior by the RelRes algorithm. This, in turn, allows the scientific community to readily determine the ideal switching distance for an arbitrary RelRes system.
Collapse
Affiliation(s)
- Mark Chaimovich
- Russian School of Mathematics, North Bethesda, Maryland 20852, United States
| | - Aviel Chaimovich
- Department of Chemical and Biological Engineering, Drexel University, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
3
|
Petix CL, Fakhraei M, Kieslich CA, Howard MP. Surrogate Modeling of the Relative Entropy for Inverse Design Using Smolyak Sparse Grids. J Chem Theory Comput 2024; 20:1538-1546. [PMID: 37703086 DOI: 10.1021/acs.jctc.3c00651] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/14/2023]
Abstract
Relative entropy minimization, a statistical-mechanics approach for finding potential energy functions that produce target structural ensembles, has proven to be a powerful strategy for the inverse design of nanoparticle self-assembly. For a given target structure, the gradient of the relative entropy with respect to the adjustable parameters of the potential energy function is computed by performing a simulation, and then these parameters are updated using iterative gradient-based optimization. Small parameter updates per iteration and many iterations can be required for numerical stability, but this incurs considerable computational expense because a new simulation must be performed to reevaluate the gradient at each iteration. Here, we investigate the use of surrogate modeling to decouple the process of minimizing the relative entropy from the computationally demanding process of determining its gradient. We approximate the relative-entropy gradient using Chebyshev polynomial interpolation on Smolyak sparse grids. Our approach potentially increases the robustness and computational efficiency of using the relative entropy for inverse design, primarily for physically informed potential energy functions that have a small number of adjustable parameters.
Collapse
Affiliation(s)
- C Levi Petix
- Department of Chemical Engineering, Auburn University, Auburn, Alabama 36849, United States
| | - Mohammadreza Fakhraei
- Department of Chemical Engineering, Auburn University, Auburn, Alabama 36849, United States
| | - Chris A Kieslich
- Department of Chemical Engineering, Auburn University, Auburn, Alabama 36849, United States
| | - Michael P Howard
- Department of Chemical Engineering, Auburn University, Auburn, Alabama 36849, United States
| |
Collapse
|
4
|
DeLuca M, Sensale S, Lin PA, Arya G. Prediction and Control in DNA Nanotechnology. ACS APPLIED BIO MATERIALS 2024; 7:626-645. [PMID: 36880799 DOI: 10.1021/acsabm.2c01045] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/08/2023]
Abstract
DNA nanotechnology is a rapidly developing field that uses DNA as a building material for nanoscale structures. Key to the field's development has been the ability to accurately describe the behavior of DNA nanostructures using simulations and other modeling techniques. In this Review, we present various aspects of prediction and control in DNA nanotechnology, including the various scales of molecular simulation, statistical mechanics, kinetic modeling, continuum mechanics, and other prediction methods. We also address the current uses of artificial intelligence and machine learning in DNA nanotechnology. We discuss how experiments and modeling are synergistically combined to provide control over device behavior, allowing scientists to design molecular structures and dynamic devices with confidence that they will function as intended. Finally, we identify processes and scenarios where DNA nanotechnology lacks sufficient prediction ability and suggest possible solutions to these weak areas.
Collapse
Affiliation(s)
- Marcello DeLuca
- Thomas Lord Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, United States
| | - Sebastian Sensale
- Department of Physics, Cleveland State University, Cleveland, Ohio 44115, United States
| | - Po-An Lin
- Thomas Lord Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, United States
| | - Gaurav Arya
- Thomas Lord Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, United States
| |
Collapse
|
5
|
Jin J, Reichman DR. Perturbative Expansion in Reciprocal Space: Bridging Microscopic and Mesoscopic Descriptions of Molecular Interactions. J Phys Chem B 2024; 128:1061-1078. [PMID: 38232134 DOI: 10.1021/acs.jpcb.3c06048] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2024]
Abstract
Determining the Fourier representation of various molecular interactions is important for constructing density-based field theories from a microscopic point of view, enabling a multiscale bridge between microscopic and mesoscopic descriptions. However, due to the strongly repulsive nature of short-ranged interactions, interparticle interactions cannot be formally defined in Fourier space, which renders coarse-grained (CG) approaches in k-space somewhat ambiguous. In this paper, we address this issue by designing a perturbative expansion of pair interactions in reciprocal space. Our perturbation theory, starting from reciprocal space, elucidates the microscopic origins underlying zeroth-order (long-range attractions) and divergent repulsive interactions from higher order contributions. We propose a systematic framework for constructing a faithful Fourier-space representation of molecular interactions, capturing key structural correlations in various systems, including simple model systems and molecular CG models of liquids. Building upon the Ornstein-Zernike equation, our approach can be combined with appropriate closure relations, and to further improve the closure approximations, we develop a bottom-up parameterization strategy for inferring the bridge function from microscopic statistics. By incorporating the bridge function into the Fourier representation, our findings suggest a systematic, bottom-up approach to performing coarse-graining in reciprocal space, leading to the systematic construction of a bottom-up classical field theory of complex aqueous systems.
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
|
6
|
Nguyen MVT, Dolph K, Delaney KT, Shen K, Sherck N, Köhler S, Gupta R, Francis MB, Shell MS, Fredrickson GH. Molecularly informed field theory for estimating critical micelle concentrations of intrinsically disordered protein surfactants. J Chem Phys 2023; 159:244904. [PMID: 38149742 PMCID: PMC10754628 DOI: 10.1063/5.0178910] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2023] [Accepted: 11/30/2023] [Indexed: 12/28/2023] Open
Abstract
The critical micelle concentration (CMC) is a crucial parameter in understanding the self-assembly behavior of surfactants. In this study, we combine simulation and experiment to demonstrate the predictive capability of molecularly informed field theories in estimating the CMC of biologically based protein surfactants. Our simulation approach combines the relative entropy coarse-graining of small-scale atomistic simulations with large-scale field-theoretic simulations, allowing us to efficiently compute the free energy of micelle formation necessary for the CMC calculation while preserving chemistry-specific information about the underlying surfactant building blocks. We apply this methodology to a unique intrinsically disordered protein platform capable of a wide variety of tailored sequences that enable tunable micelle self-assembly. The computational predictions of the CMC closely match experimental measurements, demonstrating the potential of molecularly informed field theories as a valuable tool to investigate self-assembly in bio-based macromolecules systematically.
Collapse
Affiliation(s)
- My. V. T. Nguyen
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, USA
| | - Kate Dolph
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Kris T. Delaney
- Materials Research Laboratory, University of California, Santa Barbara, California 93106, USA
| | | | | | | | - Rohini Gupta
- California Research Alliance (CARA) by BASF, Berkeley, California 94720, USA
| | | | - M. Scott Shell
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, USA
| | | |
Collapse
|
7
|
Jin J, Lee EK, Voth GA. Understanding dynamics in coarse-grained models. III. Roles of rotational motion and translation-rotation coupling in coarse-grained dynamics. J Chem Phys 2023; 159:164102. [PMID: 37870140 DOI: 10.1063/5.0167158] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2023] [Accepted: 09/29/2023] [Indexed: 10/24/2023] Open
Abstract
This paper series aims to establish a complete correspondence between fine-grained (FG) and coarse-grained (CG) dynamics by way of excess entropy scaling (introduced in Paper I). While Paper II successfully captured translational motions in CG systems using a hard sphere mapping, the absence of rotational motions in single-site CG models introduces differences between FG and CG dynamics. In this third paper, our objective is to faithfully recover atomistic diffusion coefficients from CG dynamics by incorporating rotational dynamics. By extracting FG rotational diffusion, we unravel, for the first time reported to our knowledge, a universality in excess entropy scaling between the rotational and translational diffusion. Once the missing rotational dynamics are integrated into the CG translational dynamics, an effective translation-rotation coupling becomes essential. We propose two different approaches for estimating this coupling parameter: the rough hard sphere theory with acentric factor (temperature-independent) or the rough Lennard-Jones model with CG attractions (temperature-dependent). Altogether, we demonstrate that FG diffusion coefficients can be recovered from CG diffusion coefficients by (1) incorporating "entropy-free" rotational diffusion with translation-rotation coupling and (2) recapturing the missing entropy. Our findings shed light on the fundamental relationship between FG and CG dynamics in 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
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | - Eok Kyun Lee
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, South Korea
| | - 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
|
8
|
Nadkarni I, Wu H, Aluru NR. Data-Driven Approach to Coarse-Graining Simple Liquids in Confinement. J Chem Theory Comput 2023; 19:7358-7370. [PMID: 37791529 DOI: 10.1021/acs.jctc.3c00633] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
We propose a data-driven framework for identifying coarse-grained (CG) Lennard-Jones (LJ) potential parameters in confined systems for simple liquids. Our approach involves the use of a Deep Neural Network (DNN) that is trained to approximate the solution of the Inverse Liquid State (ILST) problem for confined systems. The DNN model inherently incorporates essential physical characteristics specific to confined fluids, enabling an accurate prediction of inhomogeneity effects. By utilizing transfer learning, we predict single-site LJ potentials of simple multiatomic liquids confined in a slit-like channel, which effectively replicate both the fluid structure and molecular force of the target All-Atom (AA) system when the electrostatic interactions are not dominant. In addition, we showcase the synergy between the data-driven approach and the well-known Bottom-Up coarse-graining method utilizing Relative-Entropy (RE) Minimization. Through the sequential utilization of these two methods, the robustness of the iterative RE method is significantly augmented, leading to a remarkable enhancement in convergence.
Collapse
Affiliation(s)
- Ishan Nadkarni
- Walker Department of Mechanical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Haiyi Wu
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Narayana R Aluru
- Walker Department of Mechanical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
9
|
Peng Y, Pak AJ, Durumeric AEP, Sahrmann PG, Mani S, Jin J, Loose TD, Beiter J, Voth GA. OpenMSCG: A Software Tool for Bottom-Up Coarse-Graining. J Phys Chem B 2023; 127:8537-8550. [PMID: 37791670 PMCID: PMC10577682 DOI: 10.1021/acs.jpcb.3c04473] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2023] [Revised: 09/05/2023] [Indexed: 10/05/2023]
Abstract
The "bottom-up" approach to coarse-graining, for building accurate and efficient computational models to simulate large-scale and complex phenomena and processes, is an important approach in computational chemistry, biophysics, and materials science. As one example, the Multiscale Coarse-Graining (MS-CG) approach to developing CG models can be rigorously derived using statistical mechanics applied to fine-grained, i.e., all-atom simulation data for a given system. Under a number of circumstances, a systematic procedure, such as MS-CG modeling, is particularly valuable. Here, we present the development of the OpenMSCG software, a modularized open-source software that provides a collection of successful and widely applied bottom-up CG methods, including Boltzmann Inversion (BI), Force-Matching (FM), Ultra-Coarse-Graining (UCG), Relative Entropy Minimization (REM), Essential Dynamics Coarse-Graining (EDCG), and Heterogeneous Elastic Network Modeling (HeteroENM). OpenMSCG is a high-performance and comprehensive toolset that can be used to derive CG models from large-scale fine-grained simulation data in file formats from common molecular dynamics (MD) software packages, such as GROMACS, LAMMPS, and NAMD. OpenMSCG is modularized in the Python programming framework, which allows users to create and customize modeling "recipes" for reproducible results, thus greatly improving the reliability, reproducibility, and sharing of bottom-up CG models and their applications.
Collapse
Affiliation(s)
- Yuxing Peng
- NVIDIA
Corporation, 2788 San Tomas Expressway, Santa Clara, California 95051, United States
| | - Alexander J. Pak
- Department
of Chemical and Biological Engineering, Colorado School of Mines, Golden, Colorado 80401, United States
| | | | - Patrick G. Sahrmann
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, James Franck
Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| | - Sriramvignesh Mani
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, James Franck
Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| | - Jaehyeok Jin
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, James Franck
Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| | - Timothy D. Loose
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, James Franck
Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| | - Jeriann Beiter
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, James Franck
Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| | - Gregory A. Voth
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, James Franck
Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, United States
| |
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
|
Thaler S, Stupp M, Zavadlav J. Deep coarse-grained potentials via relative entropy minimization. J Chem Phys 2022; 157:244103. [PMID: 36586977 DOI: 10.1063/5.0124538] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022] Open
Abstract
Neural network (NN) potentials are a natural choice for coarse-grained (CG) models. Their many-body capacity allows highly accurate approximations of the potential of mean force, promising CG simulations of unprecedented accuracy. CG NN potentials trained bottom-up via force matching (FM), however, suffer from finite data effects: They rely on prior potentials for physically sound predictions outside the training data domain, and the corresponding free energy surface is sensitive to errors in the transition regions. The standard alternative to FM for classical potentials is relative entropy (RE) minimization, which has not yet been applied to NN potentials. In this work, we demonstrate, for benchmark problems of liquid water and alanine dipeptide, that RE training is more data efficient, due to accessing the CG distribution during training, resulting in improved free energy surfaces and reduced sensitivity to prior potentials. In addition, RE learns to correct time integration errors, allowing larger time steps in CG molecular dynamics simulation, while maintaining accuracy. Thus, our findings support the use of training objectives beyond FM, as a promising direction for improving CG NN potential's accuracy and reliability.
Collapse
Affiliation(s)
- Stephan Thaler
- Multiscale Modeling of Fluid Materials, Department of Engineering Physics and Computation, TUM School of Engineering and Design, Technical University of Munich, Munich, Germany
| | - Maximilian Stupp
- Multiscale Modeling of Fluid Materials, Department of Engineering Physics and Computation, TUM School of Engineering and Design, Technical University of Munich, Munich, Germany
| | - Julija Zavadlav
- Multiscale Modeling of Fluid Materials, Department of Engineering Physics and Computation, TUM School of Engineering and Design, Technical University of Munich, Munich, Germany
| |
Collapse
|
12
|
Sahu S, Schwindt NS, Coscia BJ, Shirts MR. Obtaining and Characterizing Stable Bicontinuous Cubic Morphologies and Their Nanochannels in Lyotropic Liquid Crystal Membranes. J Phys Chem B 2022; 126:10098-10110. [PMID: 36417348 DOI: 10.1021/acs.jpcb.2c06119] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
Amphiphilic monomers in polar solvents can self-assemble into lyotropic liquid crystal (LLC) bicontinuous cubic structures under the right composition and temperature conditions. After cross-linking, the resulting polymer membranes with three-dimensional (3D) continuous uniform channels are excellent candidates for filtration applications. Designing such membranes with the desired physical and chemical properties requires molecular-level understanding of the structure, which can be obtained through molecular modeling. However, building molecular models of bicontinuous cubic structures is challenging due to their narrow regime of stability and the difficulty of self-assembly of large unit cells in molecular simulations. We developed a protocol for building stable bicontinuous cubic unit cells involving both parameterization and assembly of the components. We validate the theoretical structure against experimental results for one such LLC monomer and provide insight into the structure missing in experimental data, as well as demonstrate the qualitative nature of water and solute transport through these membranes.
Collapse
Affiliation(s)
- Subin Sahu
- Department of Chemical & Biological Engineering, University of Colorado Boulder, Boulder, Colorado80309, United States
| | - Nathanael S Schwindt
- Department of Chemical & Biological Engineering, University of Colorado Boulder, Boulder, Colorado80309, United States
| | - Benjamin J Coscia
- Department of Chemical & Biological Engineering, University of Colorado Boulder, Boulder, Colorado80309, United States
| | - Michael R Shirts
- Department of Chemical & Biological Engineering, University of Colorado Boulder, Boulder, Colorado80309, United States
| |
Collapse
|
13
|
Nguyen M, Sherck N, Shen K, Edwards CER, Yoo B, Köhler S, Speros JC, Helgeson ME, Delaney KT, Shell MS, Fredrickson GH. Predicting Polyelectrolyte Coacervation from a Molecularly Informed Field-Theoretic Model. Macromolecules 2022. [DOI: 10.1021/acs.macromol.2c01759] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- My Nguyen
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, United States
| | - Nicholas Sherck
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, United States
| | - Kevin Shen
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, United States
| | - Chelsea E. R. Edwards
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, United States
- Materials Research Laboratory, University of California, Santa Barbara, California 93106, United States
| | - Brian Yoo
- BASF Corporation, Tarrytown, New York 10591, United States
| | | | - Joshua C. Speros
- California Research Alliance (CARA) by BASF, Berkeley, California 94720, United States
| | - Matthew E. Helgeson
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, United States
- Materials Research Laboratory, University of California, Santa Barbara, California 93106, United States
| | - Kris T. Delaney
- Materials Research Laboratory, University of California, Santa Barbara, California 93106, United States
| | - M. Scott Shell
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, United States
| | - Glenn H. Fredrickson
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, United States
- Materials Research Laboratory, University of California, Santa Barbara, California 93106, United States
- Department of Materials, University of California, Santa Barbara, California 93106, United States
| |
Collapse
|
14
|
Jin J, Pak AJ, Durumeric AEP, Loose TD, Voth GA. Bottom-up Coarse-Graining: Principles and Perspectives. J Chem Theory Comput 2022; 18:5759-5791. [PMID: 36070494 PMCID: PMC9558379 DOI: 10.1021/acs.jctc.2c00643] [Citation(s) in RCA: 79] [Impact Index Per Article: 39.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2022] [Indexed: 01/14/2023]
Abstract
Large-scale computational molecular models provide scientists a means to investigate the effect of microscopic details on emergent mesoscopic behavior. Elucidating the relationship between variations on the molecular scale and macroscopic observable properties facilitates an understanding of the molecular interactions driving the properties of real world materials and complex systems (e.g., those found in biology, chemistry, and materials science). As a result, discovering an explicit, systematic connection between microscopic nature and emergent mesoscopic behavior is a fundamental goal for this type of investigation. The molecular forces critical to driving the behavior of complex heterogeneous systems are often unclear. More problematically, simulations of representative model systems are often prohibitively expensive from both spatial and temporal perspectives, impeding straightforward investigations over possible hypotheses characterizing molecular behavior. While the reduction in resolution of a study, such as moving from an atomistic simulation to that of the resolution of large coarse-grained (CG) groups of atoms, can partially ameliorate the cost of individual simulations, the relationship between the proposed microscopic details and this intermediate resolution is nontrivial and presents new obstacles to study. Small portions of these complex systems can be realistically simulated. Alone, these smaller simulations likely do not provide insight into collectively emergent behavior. However, by proposing that the driving forces in both smaller and larger systems (containing many related copies of the smaller system) have an explicit connection, systematic bottom-up CG techniques can be used to transfer CG hypotheses discovered using a smaller scale system to a larger system of primary interest. The proposed connection between different CG systems is prescribed by (i) the CG representation (mapping) and (ii) the functional form and parameters used to represent the CG energetics, which approximate potentials of mean force (PMFs). As a result, the design of CG methods that facilitate a variety of physically relevant representations, approximations, and force fields is critical to moving the frontier of systematic CG forward. Crucially, the proposed connection between the system used for parametrization and the system of interest is orthogonal to the optimization used to approximate the potential of mean force present in all systematic CG methods. The empirical efficacy of machine learning techniques on a variety of tasks provides strong motivation to consider these approaches for approximating the PMF and analyzing these approximations.
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, United States
| | - 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, United States
| | - Aleksander E. P. Durumeric
- 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
| | - Timothy D. Loose
- 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
| | - 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
|
15
|
Abstract
Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.
Collapse
Affiliation(s)
- Serge Perez
- Centre de Recherche sur les Macromolecules Vegetales, University of Grenoble-Alpes, Centre National de la Recherche Scientifique, Grenoble F-38041, France
| | - Olga Makshakova
- FRC Kazan Scientific Center of Russian Academy of Sciences, Kazan Institute of Biochemistry and Biophysics, Kazan 420111, Russia
| |
Collapse
|
16
|
Wang J, in ’t Veld PJ, Robbins MO, Ge T. Effects of Coarse-Graining on Molecular Simulation of Craze Formation in Polymer Glass. Macromolecules 2022. [DOI: 10.1021/acs.macromol.1c01969] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jiuling Wang
- Department of Chemistry and Biochemistry, University of South Carolina, Columbia, South Carolina 29208, United States
| | | | - Mark O. Robbins
- Department of Physics and Astronomy, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Ting Ge
- Department of Chemistry and Biochemistry, University of South Carolina, Columbia, South Carolina 29208, United States
| |
Collapse
|
17
|
Giulini M, Rigoli M, Mattiotti G, Menichetti R, Tarenzi T, Fiorentini R, Potestio R. From System Modeling to System Analysis: The Impact of Resolution Level and Resolution Distribution in the Computer-Aided Investigation of Biomolecules. Front Mol Biosci 2021; 8:676976. [PMID: 34164432 PMCID: PMC8215203 DOI: 10.3389/fmolb.2021.676976] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2021] [Accepted: 05/06/2021] [Indexed: 12/18/2022] Open
Abstract
The ever increasing computer power, together with the improved accuracy of atomistic force fields, enables researchers to investigate biological systems at the molecular level with remarkable detail. However, the relevant length and time scales of many processes of interest are still hardly within reach even for state-of-the-art hardware, thus leaving important questions often unanswered. The computer-aided investigation of many biological physics problems thus largely benefits from the usage of coarse-grained models, that is, simplified representations of a molecule at a level of resolution that is lower than atomistic. A plethora of coarse-grained models have been developed, which differ most notably in their granularity; this latter aspect determines one of the crucial open issues in the field, i.e. the identification of an optimal degree of coarsening, which enables the greatest simplification at the expenses of the smallest information loss. In this review, we present the problem of coarse-grained modeling in biophysics from the viewpoint of system representation and information content. In particular, we discuss two distinct yet complementary aspects of protein modeling: on the one hand, the relationship between the resolution of a model and its capacity of accurately reproducing the properties of interest; on the other hand, the possibility of employing a lower resolution description of a detailed model to extract simple, useful, and intelligible information from the latter.
Collapse
Affiliation(s)
- Marco Giulini
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| | - Marta Rigoli
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| | - Giovanni Mattiotti
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| | - Roberto Menichetti
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| | - Thomas Tarenzi
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| | - Raffaele Fiorentini
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| | - Raffaello Potestio
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| |
Collapse
|
18
|
Flachmüller A, Mecking S, Peter C. Coarse grained simulation of the aggregation and structure control of polyethylene nanocrystals. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2021; 33:264001. [PMID: 33857931 DOI: 10.1088/1361-648x/abf881] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/20/2020] [Accepted: 04/15/2021] [Indexed: 06/12/2023]
Abstract
Polyethylene (PE) telechelics with carboxylate functional groups at both ends have been shown to assemble into hexagonal nanocrystal platelets with a height defined by their chain length in basic CsOH-solution. In this coarse grained (CG) simulation study we show how properties of the functional groups alter the aggregation and crystallization behavior of those telechelics. Systematic variation of the parameters of the CG model showed that important factors which control nanoparticle stability and structure are the PE chain length and the hydrophilicity and the steric demand of the head groups. To characterize the aggregation process we analyzed the number and size of the obtained aggregates as well as intramolecular order and intermolecular alignment of the polymer chains. By comparison of CG and atomistic simulation data, it could be shown that atomistic simulations representing different chemical systems can be emulated with specific, different CG parameter sets. Thus, the results from the (generic) CG simulation models can be used to explain the effect of different head groups and different counterions on the aggregation of PE telechelics and the order of the obtained nanocrystals.
Collapse
Affiliation(s)
| | - Stefan Mecking
- Department of Chemistry, Universität Konstanz, Konstanz, Germany
| | - Christine Peter
- Department of Chemistry, Universität Konstanz, Konstanz, Germany
| |
Collapse
|
19
|
Sherck N, Shen K, Nguyen M, Yoo B, Köhler S, Speros JC, Delaney KT, Shell MS, Fredrickson GH. Molecularly Informed Field Theories from Bottom-up Coarse-Graining. ACS Macro Lett 2021; 10:576-583. [PMID: 35570772 DOI: 10.1021/acsmacrolett.1c00013] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Polymer formulations possessing mesostructures or phase coexistence are challenging to simulate using atomistic particle-explicit approaches due to the disparate time and length scales, while the predictive capability of field-based simulations is hampered by the need to specify interactions at a coarser scale (e.g., χ-parameters). To overcome the weaknesses of both, we introduce a bottom-up coarse-graining methodology that leverages all-atom molecular dynamics to molecularly inform coarser field-theoretic models. Specifically, we use relative-entropy coarse-graining to parametrize particle models that are directly and analytically transformable into statistical field theories. We demonstrate the predictive capability of this approach by reproducing experimental aqueous poly(ethylene oxide) (PEO) cloud-point curves with no parameters fit to experimental data. This synergistic approach to multiscale polymer simulations opens the door to de novo exploration of phase behavior across a wide variety of polymer solutions and melt formulations.
Collapse
Affiliation(s)
- Nicholas Sherck
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, United States
| | - Kevin Shen
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, United States
- Materials Research Laboratory, University of California, Santa Barbara, California 93106, United States
| | - My Nguyen
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, United States
| | - Brian Yoo
- BASF Corporation, Tarrytown, New York 10591, United States
| | | | - Joshua C. Speros
- California Research Alliance (CARA) by BASF, Berkeley, California 94720, United States
| | - Kris T. Delaney
- Materials Research Laboratory, University of California, Santa Barbara, California 93106, United States
| | - M. Scott Shell
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, United States
| | - Glenn H. Fredrickson
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, United States
- Materials Research Laboratory, University of California, Santa Barbara, California 93106, United States
- Department of Materials, University of California, Santa Barbara, California 93106, United States
| |
Collapse
|
20
|
Sun T, Minhas V, Korolev N, Mirzoev A, Lyubartsev AP, Nordenskiöld L. Bottom-Up Coarse-Grained Modeling of DNA. Front Mol Biosci 2021; 8:645527. [PMID: 33816559 PMCID: PMC8010198 DOI: 10.3389/fmolb.2021.645527] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2020] [Accepted: 02/22/2021] [Indexed: 12/22/2022] Open
Abstract
Recent advances in methodology enable effective coarse-grained modeling of deoxyribonucleic acid (DNA) based on underlying atomistic force field simulations. The so-called bottom-up coarse-graining practice separates fast and slow dynamic processes in molecular systems by averaging out fast degrees of freedom represented by the underlying fine-grained model. The resulting effective potential of interaction includes the contribution from fast degrees of freedom effectively in the form of potential of mean force. The pair-wise additive potential is usually adopted to construct the coarse-grained Hamiltonian for its efficiency in a computer simulation. In this review, we present a few well-developed bottom-up coarse-graining methods, discussing their application in modeling DNA properties such as DNA flexibility (persistence length), conformation, "melting," and DNA condensation.
Collapse
Affiliation(s)
- Tiedong Sun
- School of Biological Sciences, Nanyang Technological University, Singapore, Singapore
| | - Vishal Minhas
- School of Biological Sciences, Nanyang Technological University, Singapore, Singapore
| | - Nikolay Korolev
- School of Biological Sciences, Nanyang Technological University, Singapore, Singapore
| | - Alexander Mirzoev
- School of Biological Sciences, Nanyang Technological University, Singapore, Singapore
| | - Alexander P. Lyubartsev
- Department of Materials and Environmental Chemistry, Stockholm University, Stockholm, Sweden
| | - Lars Nordenskiöld
- School of Biological Sciences, Nanyang Technological University, Singapore, Singapore
| |
Collapse
|
21
|
Johnson LC, Phelan FR. Dynamically consistent coarse-grain simulation model of chemically specific polymer melts via friction parameterization. J Chem Phys 2021; 154:084114. [PMID: 33639746 PMCID: PMC10075510 DOI: 10.1063/5.0034910] [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/14/2022] Open
Abstract
Coarse-grained (CG) models of polymers involve grouping many atoms in an all-atom (AA) representation into single sites to reduce computational effort yet retain the hierarchy of length and time scales inherent to macromolecules. Parameterization of such models is often via "bottom-up" methods, which preserve chemical specificity but suffer from artificially accelerated dynamics with respect to the AA model from which they were derived. Here, we study the combination of a bottom-up CG model with a dissipative potential as a means to obtain a chemically specific and dynamically correct model. We generate the conservative part of the force-field using the iterative Boltzmann inversion (IBI) method, which seeks to recover the AA structure. This is augmented with the dissipative Langevin thermostat, which introduces a single parameterizable friction factor to correct the unphysically fast dynamics of the IBI-generated force-field. We study this approach for linear polystyrene oligomer melts for three separate systems with 11, 21, and 41 monomers per chain and a mapping of one monomer per CG site. To parameterize the friction factor, target values are extracted from the AA dynamics using translational monomer diffusion, translational chain diffusion, and rotational chain motion to test the consistency of the parameterization across different modes of motion. We find that the value of the friction parameter needed to bring the CG dynamics in line with AA target values varies based on the mode of parameterization with short-time monomer translational dynamics requiring the highest values, long-time chain translational dynamics requiring the lowest values, and rotational dynamics falling in between. The friction ranges most widely for the shortest chains, and the span narrows with increasing chain length. For longer chains, a practical working value of the friction parameter may be derived from the rotational dynamics, owing to the contribution of multiple relaxation modes to chain rotation and a lack of sensitivity of the translational dynamics at these intermediate levels of friction. A study of equilibrium chain structure reveals that all chains studied are non-Gaussian. However, longer chains better approximate ideal chain dimensions than more rod-like shorter chains and thus are most closely described by a single friction parameter. We also find that the separability of the conservative and dissipative potentials is preserved.
Collapse
Affiliation(s)
- Lilian C Johnson
- Materials Science and Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| | - Frederick R Phelan
- Materials Science and Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| |
Collapse
|
22
|
Schmidt M, Schroeder I, Bauer D, Thiel G, Hamacher K. Inferring functional units in ion channel pores via relative entropy. EUROPEAN BIOPHYSICS JOURNAL : EBJ 2021; 50:37-57. [PMID: 33523249 PMCID: PMC7872957 DOI: 10.1007/s00249-020-01480-7] [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] [Subscribe] [Scholar Register] [Received: 10/21/2019] [Revised: 10/11/2020] [Accepted: 11/09/2020] [Indexed: 11/25/2022]
Abstract
Coarse-grained protein models approximate the first-principle physical potentials. Among those modeling approaches, the relative entropy framework yields promising and physically sound results, in which a mapping from the target protein structure and dynamics to a model is defined and subsequently adjusted by an entropy minimization of the model parameters. Minimization of the relative entropy is equivalent to maximization of the likelihood of reproduction of (configurational ensemble) observations by the model. In this study, we extend the relative entropy minimization procedure beyond parameter fitting by a second optimization level, which identifies the optimal mapping to a (dimension-reduced) topology. We consider anisotropic network models of a diverse set of ion channels and assess our findings by comparison to experimental results.
Collapse
Affiliation(s)
- Michael Schmidt
- Department of Physics, TU Darmstadt, Karolinenpl. 5, 64289 Darmstadt, Germany
| | - Indra Schroeder
- Department of Biology, TU Darmstadt, Schnittspahnstr. 10, 64287 Darmstadt, Germany
| | - Daniel Bauer
- Department of Biology, TU Darmstadt, Schnittspahnstr. 10, 64287 Darmstadt, Germany
| | - Gerhard Thiel
- Department of Biology, TU Darmstadt, Schnittspahnstr. 10, 64287 Darmstadt, Germany
| | - Kay Hamacher
- Department of Physics, Department of Biology, Department of Computer Science, TU Darmstadt, Schnittspahnstr. 10, 64287 Darmstadt, Germany
| |
Collapse
|
23
|
Shen K, Sherck N, Nguyen M, Yoo B, Köhler S, Speros J, Delaney KT, Fredrickson GH, Shell MS. Learning composition-transferable coarse-grained models: Designing external potential ensembles to maximize thermodynamic information. J Chem Phys 2020; 153:154116. [DOI: 10.1063/5.0022808] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Affiliation(s)
- Kevin Shen
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, USA
- Materials Research Laboratory, University of California, Santa Barbara, California 93106, USA
| | - Nicholas Sherck
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, USA
| | - My Nguyen
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, USA
| | - Brian Yoo
- BASF Corporation, Tarrytown, New York 10591, USA
| | | | - Joshua Speros
- California Research Alliance (CARA) by BASF, Berkeley, California 94720, USA
| | - Kris T. Delaney
- Materials Research Laboratory, University of California, Santa Barbara, California 93106, USA
| | - Glenn H. Fredrickson
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, USA
- Materials Research Laboratory, University of California, Santa Barbara, California 93106, USA
- Department of Materials Engineering, University of California, Santa Barbara, California 93106, USA
| | - M. Scott Shell
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, USA
| |
Collapse
|
24
|
Jin J, Yu A, Voth GA. Temperature and Phase Transferable Bottom-up Coarse-Grained Models. J Chem Theory Comput 2020; 16:6823-6842. [PMID: 32975948 DOI: 10.1021/acs.jctc.0c00832] [Citation(s) in RCA: 30] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
Despite the high fidelity of bottom-up coarse-grained (CG) approaches to recapitulate the structural correlations in atomistic simulations, the general use of bottom-up CG methods is limited because of the nontransferable nature of these CG models under different thermodynamic conditions. Because bottom-up CG potentials usually correspond to configuration-dependent free energies of the system, recent studies have focused on adjusting enthalpic or entropic contributions to account for issues with transferability. However, these approaches can require a manual adjustment of the CG interaction a priori and are usually limited to constant volume ensembles. To overcome these limitations, we construct temperature and phase transferable CG models under constant pressure by developing the ultra-coarse-graining (UCG) methodology in the mean-field limit. In the mean-field ansatz, an embedded semi-global order parameter recapitulates global changes to the system by automatically adjusting the effective CG interactions, thus bridging free energy decompositions with UCG theory. The method presented is designed to faithfully capture structural correlations under different thermodynamic conditions, using a single UCG model. Specifically, we test the applicability of the developed theory in three distinct cases: (1) different temperatures at constant pressure in liquids, (2) different temperatures across thermodynamic phases, and (3) liquid/vapor interfaces. We demonstrate that the systematic construction of both temperature and phase transferable bottom-up CG models is possible using this generalized UCG theory. Based on our findings, this approach significantly extends the transferability and applicability of the bottom-up CG theory and method.
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 Avenue, Chicago, Illinois 60637, United States
| | - Alvin Yu
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, 5735 S. Ellis Avenue, Chicago, Illinois 60637, United States
| | - 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 Avenue, Chicago, Illinois 60637, United States
| |
Collapse
|
25
|
Gkeka P, Stoltz G, Barati Farimani A, Belkacemi Z, Ceriotti M, Chodera JD, Dinner AR, Ferguson AL, Maillet JB, Minoux H, Peter C, Pietrucci F, Silveira A, Tkatchenko A, Trstanova Z, Wiewiora R, Lelièvre T. Machine Learning Force Fields and Coarse-Grained Variables in Molecular Dynamics: Application to Materials and Biological Systems. J Chem Theory Comput 2020; 16:4757-4775. [PMID: 32559068 PMCID: PMC8312194 DOI: 10.1021/acs.jctc.0c00355] [Citation(s) in RCA: 82] [Impact Index Per Article: 20.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Machine learning encompasses tools and algorithms that are now becoming popular in almost all scientific and technological fields. This is true for molecular dynamics as well, where machine learning offers promises of extracting valuable information from the enormous amounts of data generated by simulation of complex systems. We provide here a review of our current understanding of goals, benefits, and limitations of machine learning techniques for computational studies on atomistic systems, focusing on the construction of empirical force fields from ab initio databases and the determination of reaction coordinates for free energy computation and enhanced sampling.
Collapse
Affiliation(s)
- Paraskevi Gkeka
- Integrated Drug Discovery, Sanofi R&D, 91385 Chilly-Mazarin, France
| | - Gabriel Stoltz
- CERMICS, Ecole des Ponts, Marne-la-Vallée, France
- Matherials Project-Team, Inria Paris, 75012 Paris, France
| | | | - Zineb Belkacemi
- Integrated Drug Discovery, Sanofi R&D, 91385 Chilly-Mazarin, France
- CERMICS, Ecole des Ponts, Marne-la-Vallée, France
| | - Michele Ceriotti
- Laboratory of Computational Science and Modelling, Institute of Materials, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - John D Chodera
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, New York 10065, United States
| | - Aaron R Dinner
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Andrew L Ferguson
- Pritzker School of Molecular Engineering, University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637, United States
| | | | - Hervé Minoux
- Integrated Drug Discovery, Sanofi R&D, 94403 Vitry-sur-Seine, France
| | | | - Fabio Pietrucci
- UMR CNRS 7590, MNHN, Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Sorbonne Université, 75005 Paris, France
| | - Ana Silveira
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, New York 10065, United States
| | - Alexandre Tkatchenko
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Zofia Trstanova
- School of Mathematics, The University of Edinburgh, Edinburgh EH9 3FD, U.K
| | - Rafal Wiewiora
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, New York 10065, United States
| | - Tony Lelièvre
- CERMICS, Ecole des Ponts, Marne-la-Vallée, France
- Matherials Project-Team, Inria Paris, 75012 Paris, France
| |
Collapse
|
26
|
Role of Entropy in Colloidal Self-Assembly. ENTROPY 2020; 22:e22080877. [PMID: 33286648 PMCID: PMC7517482 DOI: 10.3390/e22080877] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/05/2020] [Revised: 08/03/2020] [Accepted: 08/05/2020] [Indexed: 12/13/2022]
Abstract
Entropy plays a key role in the self-assembly of colloidal particles. Specifically, in the case of hard particles, which do not interact or overlap with each other during the process of self-assembly, the free energy is minimized due to an increase in the entropy of the system. Understanding the contribution of entropy and engineering it is increasingly becoming central to modern colloidal self-assembly research, because the entropy serves as a guide to design a wide variety of self-assembled structures for many technological and biomedical applications. In this work, we highlight the importance of entropy in different theoretical and experimental self-assembly studies. We discuss the role of shape entropy and depletion interactions in colloidal self-assembly. We also highlight the effect of entropy in the formation of open and closed crystalline structures, as well as describe recent advances in engineering entropy to achieve targeted self-assembled structures.
Collapse
|
27
|
Abstract
Low resolution coarse-grained (CG) models are widely adopted for investigating phenomena that cannot be effectively simulated with all-atom (AA) models. Since the development of the many-body dissipative particle dynamics method, CG models have increasingly supplemented conventional pair potentials with one-body potentials of the local density (LD) around each site. These LD potentials appear to significantly extend the transferability of CG models, while also enabling more accurate descriptions of thermodynamic properties, interfacial phenomena, and many-body correlations. In this work, we systematically examine the properties of LD potentials. We first derive and numerically demonstrate a nontrivial transformation of pair and LD potentials that leaves the total forces and equilibrium distribution invariant. Consequently, the pair and LD potentials determined via bottom-up methods are not unique. We then investigate the sensitivity of CG models for glycerol to the weighting function employed for defining the local density. We employ the multiscale coarse-graining (MS-CG) method to simultaneously parameterize both pair and LD potentials. When employing a short-ranged Lucy function that defines the local density from the first solvation shell, the MS-CG model accurately reproduces the pair structure, pressure-density equation of state, and liquid-vapor interfacial profile of the AA model. The accuracy of the model generally decreases as the range of the Lucy function increases further. The MS-CG model provides similar accuracy when a smoothed Heaviside function is employed to define the local density from the first solvation shell. However, the model performs less well when this function acts on either longer or shorter length scales.
Collapse
Affiliation(s)
- Michael R 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
|
Pervaje AK, Walker CC, Santiso EE. Molecular simulation of polymers with a SAFT-γ Mie approach. MOLECULAR SIMULATION 2019. [DOI: 10.1080/08927022.2019.1645331] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- Amulya K. Pervaje
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC, USA
| | - Christopher C. Walker
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC, USA
| | - Erik E. Santiso
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC, USA
| |
Collapse
|
29
|
Sanyal T, Mittal J, Shell MS. A hybrid, bottom-up, structurally accurate, Go¯-like coarse-grained protein model. J Chem Phys 2019; 151:044111. [PMID: 31370551 PMCID: PMC6663515 DOI: 10.1063/1.5108761] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2019] [Accepted: 06/24/2019] [Indexed: 12/21/2022] Open
Abstract
Coarse-grained (CG) protein models in the structural biology literature have improved over the years from being simple tools to understand general folding and aggregation driving forces to capturing detailed structures achieved by actual folding sequences. Here, we ask whether such models can be developed systematically from recent advances in bottom-up coarse-graining methods without relying on bioinformatic data (e.g., protein data bank statistics). We use relative entropy coarse-graining to develop a hybrid CG but Go¯-like CG peptide model, hypothesizing that the landscape of proteinlike folds is encoded by the backbone interactions, while the sidechain interactions define which of these structures globally minimizes the free energy in a unique native fold. To construct a model capable of capturing varied secondary structures, we use a new extended ensemble relative entropy method to coarse-grain based on multiple reference atomistic simulations of short polypeptides with varied α and β character. Subsequently, we assess the CG model as a putative protein backbone forcefield by combining it with sidechain interactions based on native contacts but not incorporating native distances explicitly, unlike standard Go¯ models. We test the model's ability to fold a range of proteins and find that it achieves high accuracy (∼2 Å root mean square deviation resolution for both short sequences and large globular proteins), suggesting the strong role that backbone conformational preferences play in defining the fold landscape. This model can be systematically extended to non-natural amino acids and nonprotein polymers and sets the stage for extensions to non-Go¯ models with sequence-specific sidechain interactions.
Collapse
Affiliation(s)
- Tanmoy Sanyal
- Department of Chemical Engineering, University of California Santa Barbara, Santa Barbara, California 93106, USA
| | - Jeetain Mittal
- Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, USA
| | - M. Scott Shell
- Department of Chemical Engineering, University of California Santa Barbara, Santa Barbara, California 93106, USA
| |
Collapse
|
30
|
Lebold KM, Noid WG. Dual approach for effective potentials that accurately model structure and energetics. J Chem Phys 2019; 150:234107. [PMID: 31228924 DOI: 10.1063/1.5094330] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
Because they eliminate unnecessary degrees of freedom, coarse-grained (CG) models enable studies of phenomena that are intractable with more detailed models. For the same reason, the effective potentials that govern CG degrees of freedom incorporate entropic contributions from the eliminated degrees of freedom. Consequently, these effective potentials demonstrate limited transferability and provide a poor estimate of atomic energetics. Here, we propose a simple dual-potential approach that combines "structure-based" and "energy-based" variational principles to determine effective potentials that model free energies and potential energies, respectively, as a function of the CG configuration. We demonstrate this approach for 1-site CG models of water and methanol. We accurately sample configuration space by performing simulations with the structure-based potential. We accurately estimate average atomic energies by postprocessing the sampled configurations with the energy-based potential. Finally, the difference between the two potentials predicts a qualitatively accurate estimate for the temperature dependence of the structure-based potential.
Collapse
Affiliation(s)
- Kathryn M Lebold
- Department of Chemistry, Penn State University, University Park, Pennsylvania 16802, USA
| | - W G Noid
- Department of Chemistry, Penn State University, University Park, Pennsylvania 16802, USA
| |
Collapse
|
31
|
Jin J, Han Y, Voth GA. Coarse-graining involving virtual sites: Centers of symmetry coarse-graining. J Chem Phys 2019; 150:154103. [DOI: 10.1063/1.5067274] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Affiliation(s)
- Jaehyeok Jin
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Yining Han
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Gregory A. Voth
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
32
|
|
33
|
Walker CC, Genzer J, Santiso EE. Development of a fused-sphere SAFT-γ Mie force field for poly(vinyl alcohol) and poly(ethylene). J Chem Phys 2019; 150:034901. [DOI: 10.1063/1.5078742] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Affiliation(s)
- Christopher C. Walker
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27695, USA
| | - Jan Genzer
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27695, USA
| | - Erik E. Santiso
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27695, USA
| |
Collapse
|
34
|
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
|
35
|
Webb MA, Delannoy JY, de Pablo JJ. Graph-Based Approach to Systematic Molecular Coarse-Graining. J Chem Theory Comput 2018; 15:1199-1208. [PMID: 30557028 DOI: 10.1021/acs.jctc.8b00920] [Citation(s) in RCA: 46] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A novel methodology is introduced here to generate coarse-grained (CG) representations of molecular models for simulations. The proposed strategy relies on basic graph-theoretic principles and is referred to as graph-based coarse-graining (GBCG). It treats a given system as a molecular graph and derives a corresponding CG representation by using edge contractions to combine nodes in the graph, which correspond to atoms in the molecule, into CG sites. A key element of this methodology is that the nodes are combined according to well-defined protocols that rank-order nodes based on the underlying chemical connectivity. By iteratively performing these operations, successively coarser representations of the original atomic system can be produced to yield a systematic set of CG mappings with hierarchical resolution in an automated fashion. These capabilities are demonstrated in the context of several test systems, including toluene, pentadecane, a polysaccharide dimer, and a rhodopsin protein. In these examples, GBCG yields multiple, intuitive structures that naturally preserve the chemical topology of the system. Importantly, these representations are rendered from algorithmic implementation rather than an arbitrary ansatz, which, until now, has been the conventional approach for defining CG mapping schemes. Overall, the results presented here indicate that GBCG is efficient, robust, and unambiguous in its application, making it a valuable tool for future CG modeling.
Collapse
Affiliation(s)
- Michael A Webb
- Institute for Molecular Engineering , University of Chicago , Chicago , Illinois 60637 , United States
| | - Jean-Yves Delannoy
- Simulation, Modeling and Artificial Intelligence team , Solvay , Bristol , Pennsylvania 19007 , United States
| | - Juan J de Pablo
- Institute for Molecular Engineering , University of Chicago , Chicago , Illinois 60637 , United States.,Institute for Molecular Engineering and Materials Science Division , Argonne National Laboratory , Lemont , Illinois 60439 , United States
| |
Collapse
|
36
|
Abstract
This article gives an overview of excess-entropy scaling, the 1977 discovery by Rosenfeld that entropy determines properties of liquids like viscosity, diffusion constant, and heat conductivity. We give examples from computer simulations confirming this intriguing connection between dynamics and thermodynamics, counterexamples, and experimental validations. Recent uses in application-related contexts are reviewed, and theories proposed for the origin of excess-entropy scaling are briefly summarized. It is shown that if two thermodynamic state points of a liquid have the same microscopic dynamics, they must have the same excess entropy. In this case, the potential-energy function exhibits a symmetry termed hidden scale invariance, stating that the ordering of the potential energies of configurations is maintained if these are scaled uniformly to a different density. This property leads to the isomorph theory, which provides a general framework for excess-entropy scaling and illuminates, in particular, why this does not apply rigorously and universally. It remains an open question whether all aspects of excess-entropy scaling and related regularities reflect hidden scale invariance in one form or other.
Collapse
Affiliation(s)
- Jeppe C Dyre
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
37
|
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
|
38
|
Han Y, Jin J, Wagner JW, Voth GA. Quantum theory of multiscale coarse-graining. J Chem Phys 2018; 148:102335. [PMID: 29544317 DOI: 10.1063/1.5010270] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
Coarse-grained (CG) models serve as a powerful tool to simulate molecular systems at much longer temporal and spatial scales. Previously, CG models and methods have been built upon classical statistical mechanics. The present paper develops a theory and numerical methodology for coarse-graining in quantum statistical mechanics, by generalizing the multiscale coarse-graining (MS-CG) method to quantum Boltzmann statistics. A rigorous derivation of the sufficient thermodynamic consistency condition is first presented via imaginary time Feynman path integrals. It identifies the optimal choice of CG action functional and effective quantum CG (qCG) force field to generate a quantum MS-CG (qMS-CG) description of the equilibrium system that is consistent with the quantum fine-grained model projected onto the CG variables. A variational principle then provides a class of algorithms for optimally approximating the qMS-CG force fields. Specifically, a variational method based on force matching, which was also adopted in the classical MS-CG theory, is generalized to quantum Boltzmann statistics. The qMS-CG numerical algorithms and practical issues in implementing this variational minimization procedure are also discussed. Then, two numerical examples are presented to demonstrate the method. Finally, as an alternative strategy, a quasi-classical approximation for the thermal density matrix expressed in the CG variables is derived. This approach provides an interesting physical picture for coarse-graining in quantum Boltzmann statistical mechanics in which the consistency with the quantum particle delocalization is obviously manifest, and it opens up an avenue for using path integral centroid-based effective classical force fields in a coarse-graining methodology.
Collapse
Affiliation(s)
- Yining Han
- Department of Chemistry, James Frank Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Jaehyeok Jin
- Department of Chemistry, James Frank Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Jacob W Wagner
- Department of Chemistry, James Frank Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Gregory A Voth
- Department of Chemistry, James Frank Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
39
|
MacKay RS, Robinson JD. Aggregation of Markov flows I: theory. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2018; 376:rsta.2017.0232. [PMID: 29555805 PMCID: PMC5869541 DOI: 10.1098/rsta.2017.0232] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 01/08/2018] [Indexed: 05/29/2023]
Abstract
A Markov flow is a stationary measure, with the associated flows and mean first passage times, for a continuous-time regular jump homogeneous semi-Markov process on a discrete state space. Nodes in the state space can be eliminated to produce a smaller Markov flow which is a factor of the original one. Some improvements to the elimination methods of Wales are given. The main contribution of the paper is to present an alternative, namely a method to aggregate groups of nodes to produce a factor. The method can be iterated to make hierarchical aggregation schemes. The potential benefits are efficient computation, including recomputation to take into account local changes, and insights into the macroscopic behaviour.This article is part of the theme issue 'Hilbert's sixth problem'.
Collapse
Affiliation(s)
- R S MacKay
- Mathematics Institute and Centre for Complexity Science, University of Warwick, Coventry CV4 7AL, UK
| | - J D Robinson
- Mathematics Institute and Centre for Complexity Science, University of Warwick, Coventry CV4 7AL, UK
| |
Collapse
|
40
|
Xia W, Song J, Hansoge NK, Phelan FR, Keten S, Douglas JF. Energy Renormalization for Coarse-Graining the Dynamics of a Model Glass-Forming Liquid. J Phys Chem B 2018; 122:2040-2045. [PMID: 29400063 PMCID: PMC6217959 DOI: 10.1021/acs.jpcb.8b00321] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Coarse-grained modeling achieves the enhanced computational efficiency required to model glass-forming materials by integrating out "unessential" molecular degrees of freedom, but no effective temperature transferable coarse-graining method currently exists to capture dynamics. We address this fundamental problem through an energy-renormalization scheme, in conjunction with the localization model of relaxation relating the Debye-Waller factor ⟨u2⟩ to the structural relaxation time τ. Taking ortho-terphenyl as a model small-molecule glass-forming liquid, we show that preserving ⟨u2⟩ (at picosecond time scale) under coarse-graining by renormalizing the cohesive interaction strength allows for quantitative prediction of both short- and long-time dynamics covering the entire temperature range of glass formation. Our findings provide physical insights into the dynamics of cooled liquids and make progress for building temperature-transferable coarse-grained models that predict key properties of glass-forming materials.
Collapse
Affiliation(s)
- Wenjie Xia
- Materials Science & Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
- Center for Hierarchical Materials Design, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208-3109, United States
- Department of Civil & Environmental Engineering, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208-3109, United States
| | - Jake Song
- Department of Materials Science & Engineering, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208-3109, United States
| | - Nitin K. Hansoge
- Department of Mechanical Engineering, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208-3109, United States
| | - Frederick R. Phelan
- Materials Science & Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
| | - Sinan Keten
- Department of Civil & Environmental Engineering, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208-3109, United States
- Department of Mechanical Engineering, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208-3109, United States
| | - Jack F. Douglas
- Materials Science & Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
| |
Collapse
|
41
|
Song J, Hsu DD, Shull KR, Phelan FR, Douglas JF, Xia W, Keten S. Energy Renormalization Method for the Coarse-Graining of Polymer Viscoelasticity. Macromolecules 2018; 51:10.1021/acs.macromol.7b02560. [PMID: 30996476 PMCID: PMC6463302 DOI: 10.1021/acs.macromol.7b02560] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Developing temperature transferable coarse-grained (CG) models is essential for the computational prediction of polymeric glass-forming (GF) material behavior, but their dynamics are often greatly altered from those of all-atom (AA) models mainly because of the reduced fluid configurational entropy under coarse-graining. To address this issue, we have recently introduced an energy renormalization (ER) strategy that corrects the activation free energy of the CG polymer model by renormalizing the cohesive interaction strength ε as a function of temperature T, i.e., ε(T), thus semiempirically preserving the T-dependent dynamics of the AA model. Here we apply our ER method to consider-in addition to T-dependency-the frequency f-dependent polymer viscoelasticity. Through smallamplitude oscillatory shear molecular dynamics simulations, we show that changing the imposed oscillation f on the CG systems requires changes in ε values (i.e., ε(T, f)) to reproduce the AA viscoelasticity. By accounting for the dynamic fragility of polymers as a material parameter, we are able to predict ε(T, f) under coarse-graining in order to capture the AA viscoelasticity, and consequently the activation energy, across a wide range of T and f in the GF regime. Specifically, we showcase our achievements on two representative polymers of distinct fragilities, polybutadiene (PB) and polystyrene (PS), and show that our CG models are able to sample viscoelasticity up to the megahertz regime, which approaches state-of-the-art experimental resolutions, and capture results sampled via AA simulations and prior experiments.
Collapse
Affiliation(s)
- Jake Song
- Department of Materials Science & Engineering, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208-3109, United States
| | - David D. Hsu
- Department of Physics and Engineering, Wheaton College, 501 College Avenue, Wheaton, Illinois 60187, United States
| | - Kenneth R. Shull
- Department of Materials Science & Engineering, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208-3109, United States
| | - Frederick R. Phelan
- Materials Science & Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
| | - Jack F. Douglas
- Materials Science & Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
| | - Wenjie Xia
- Department of Civil & Environmental Engineering, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208-3109, United States
- Center for Hierarchical Materials Design, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208-3109, United States
- Materials Science & Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
| | - Sinan Keten
- Department of Civil & Environmental Engineering, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208-3109, United States
- Department of Mechanical Engineering, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208-3109, United States
| |
Collapse
|
42
|
Lemke T, Peter C. Neural Network Based Prediction of Conformational Free Energies - A New Route toward Coarse-Grained Simulation Models. J Chem Theory Comput 2017; 13:6213-6221. [DOI: 10.1021/acs.jctc.7b00864] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Tobias Lemke
- Theoretical Chemistry, University of Konstanz, 78547 Konstanz, Germany
| | - Christine Peter
- Theoretical Chemistry, University of Konstanz, 78547 Konstanz, Germany
| |
Collapse
|
43
|
Buslaev P, Gushchin I. Effects of Coarse Graining and Saturation of Hydrocarbon Chains on Structure and Dynamics of Simulated Lipid Molecules. Sci Rep 2017; 7:11476. [PMID: 28904383 PMCID: PMC5597592 DOI: 10.1038/s41598-017-11761-5] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2017] [Accepted: 08/25/2017] [Indexed: 12/14/2022] Open
Abstract
Molecular dynamics simulations are used extensively to study the processes on biological membranes. The simulations can be conducted at different levels of resolution: all atom (AA), where all atomistic details are provided; united atom (UA), where hydrogen atoms are treated inseparably of corresponding heavy atoms; and coarse grained (CG), where atoms are grouped into larger particles. Here, we study the behavior of model bilayers consisting of saturated and unsaturated lipids DOPC, SOPC, OSPC and DSPC in simulations performed using all atom CHARMM36 and coarse grained Martini force fields. Using principal components analysis, we show that the structural and dynamical properties of the lipids are similar, both in AA and CG simulations, although the unsaturated molecules are more dynamic and favor more extended conformations. We find that CG simulations capture 75 to 100% of the major collective motions, overestimate short range ordering, result in more flexible molecules and 5–7 fold faster sampling. We expect that the results reported here will be useful for comprehensive quantitative comparisons of simulations conducted at different resolution levels and for further development and improvement of CG force fields.
Collapse
Affiliation(s)
- Pavel Buslaev
- Moscow Institute of Physics and Technology, 141700, Dolgoprudny, Russia.
| | - Ivan Gushchin
- Moscow Institute of Physics and Technology, 141700, Dolgoprudny, Russia. .,Institute of Complex Systems (ICS), ICS-6: Structural Biochemistry, Research Centre Jülich, 52425, Jülich, Germany.
| |
Collapse
|
44
|
Xia W, Song J, Jeong C, Hsu DD, Phelan FR, Douglas JF, Keten S. Energy-Renormalization for Achieving Temperature Transferable Coarse-Graining of Polymer Dynamics. Macromolecules 2017; 50:10.1021/acs.macromol.7b01717. [PMID: 30996475 PMCID: PMC6463524 DOI: 10.1021/acs.macromol.7b01717] [Citation(s) in RCA: 67] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The bottom-up prediction of the properties of polymeric materials based on molecular dynamics simulation is a major challenge in soft matter physics. Coarse-grained (CG) models are often employed to access greater spatiotemporal scales required for many applications, but these models normally experience significantly altered thermodynamics and highly accelerated dynamics due to the reduced number of degrees of freedom upon coarse-graining. While CG models can be calibrated to meet certain properties at particular state points, there is unfortunately no temperature transferable and chemically specific coarse-graining method that allows for modeling of polymer dynamics over a wide temperature range. Here, we pragmatically address this problem by "correcting" for deviations in activation free energies that occur upon coarse-graining the dynamics of a model polymeric material (polystyrene). In particular, we propose a new strategy based on concepts drawn from the Adam-Gibbs (AG) theory of glass formation. Namely we renormalize the cohesive interaction strength and effective interaction length-scale parameters to modify the activation free energy. We show that this energy-renormalization method for CG modeling allows accurate prediction of atomistic dynamics over the Arrhenius regime, the non-Arrhenius regime of glass formation, and even the non-equilibrium glassy regime, thus allowing for the predictive modeling of dynamic properties of polymer over the entire range of glass formation. Our work provides a practical scheme for establishing temperature transferable coarse-grained models for predicting and designing the properties of polymeric materials.
Collapse
Affiliation(s)
- Wenjie Xia
- Materials Science & Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
- Department of Civil & Environmental Engineering, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208-3109, United States
- Center for Hierarchical Materials Design, Northwestern University, Evanston, Illinois 60208-3109, United States
| | - Jake Song
- Department of Materials Science & Engineering, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208-3109, United States
| | - Cheol Jeong
- Materials Science & Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
| | - David D. Hsu
- Department of Mechanical Engineering, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208-3109, United States
| | - Frederick R. Phelan
- Materials Science & Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
| | - Jack F. Douglas
- Materials Science & Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
| | - Sinan Keten
- Department of Civil & Environmental Engineering, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208-3109, United States
- Department of Mechanical Engineering, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208-3109, United States
| |
Collapse
|
45
|
Olsen AE, Dyre JC, Schrøder TB. Communication: Pseudoisomorphs in liquids with intramolecular degrees of freedom. J Chem Phys 2016; 145:241103. [DOI: 10.1063/1.4972860] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Andreas Elmerdahl Olsen
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - Jeppe C. Dyre
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| | - Thomas B. Schrøder
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|