Understanding Missing Entropy in Coarse-Grained Systems: Addressing Issues of Representability and Transferability

Understanding dynamics in coarse-grained models. V. Extension of coarse-grained dynamics theory to non-hard sphere systems

Coarse-grained (CG) modeling has gained significant attention in recent years due to its wide applicability in enhancing the spatiotemporal scales of molecular simulations. While CG simulations, often performed with Hamiltonian mechanics, faithfully recapitulate structural correlations at equilibrium, they lead to ambiguously accelerated dynamics. In Paper I [J. Jin, K. S. Schweizer, and G. A. Voth, J. Chem. Phys. 158(3), 034103 (2023)], we proposed the excess entropy scaling relationship to understand the CG dynamics. Then, in Paper II [J. Jin, K. S. Schweizer, and G. A. Voth, J. Chem. Phys. 158(3), 034104 (2023)], we developed a theory to map the CG system into a dynamically consistent hard sphere system to analytically derive an expression for fast CG dynamics. However, many chemical and physical systems do not exhibit hard sphere-like behavior, limiting the extensibility of the developed theory. In this paper, we aim to generalize the theory to the non-hard sphere system based on the Weeks-Chandler-Andersen perturbation theory. Since non-hard sphere-like CG interactions affect the excess entropy term as it deviates from the hard sphere description, we explicitly account for the extra entropy to correct the non-hard sphere nature of the system. This approach is demonstrated for two different types of interactions seen in liquids, and we further provide a generalized description for any CG models using the generalized Gaussian CG models using Gaussian basis sets. Altogether, this work allows for extending the range and applicability of the hard sphere CG dynamics theory to a myriad of CG liquids.

Understanding the coarse-grained free energy landscape of phospholipids and their phase separation

The cell membrane exhibits lateral heterogeneity due to the preferential association among the large number of lipid species that constitute the membrane. In particular, the preferential association of cholesterol (CHOL) with saturated lipids into ordered domains has been an area of intense investigation. The large spatiotemporal scales that comprise spontaneous domain formation largely precludes computational investigation via conventional all-atom molecular dynamics. We demonstrate here that molecular coarse-grained (CG) models, obtained from the bottom-up, i.e., via statistical mechanical renormalization of atomistic models, are capable of spontaneous assembly and phase separation for two model raft-like systems, DLiPC/DPPC/CHOL and DOPC/DPPC/CHOL. The resulting bottom-up CG models exhibit spontaneous self-assembly and phase separation and recapitulate the structural correlations of the underlying atomistic models. The accuracy and fast dynamics of these CG models constitute an effective means of bypassing the limited spatiotemporal scales of atomistic simulations. As the first bottom-up CG models of lipid phase separation, the CG models in this work provide an informative analysis for further construction of bottom-up CG models transferable across a range of lipid compositions.

Thermodynamic Transferability in Coarse-Grained Force Fields Using Graph Neural Networks

Coarse-graining is a molecular modeling technique in which an atomistic system is represented in a simplified fashion that retains the most significant system features that contribute to a target output while removing the degrees of freedom that are less relevant. This reduction in model complexity allows coarse-grained molecular simulations to reach increased spatial and temporal scales compared with corresponding all-atom models. A core challenge in coarse-graining is to construct a force field that represents the interactions in the new representation in a way that preserves the atomistic-level properties. Many approaches to building coarse-grained force fields have limited transferability between different thermodynamic conditions as a result of averaging over internal fluctuations at a specific thermodynamic state point. Here, we use a graph-convolutional neural network architecture, the Hierarchically Interacting Particle Neural Network with Tensor Sensitivity (HIP-NN-TS), to develop a highly automated training pipeline for coarse-grained force fields, which allows for studying the transferability of coarse-grained models based on the force-matching approach. We show that this approach yields not only highly accurate force fields but also that these force fields are more transferable through a variety of thermodynamic conditions. These results illustrate the potential of machine learning techniques, such as graph neural networks, to improve the construction of transferable coarse-grained force fields.

Changing Your Martini Can Still Give You a Hangover

The Martini 3.0 coarse-grained force field, which was parametrized to better capture transferability in top-down coarse-grained models, is analyzed to assess its accuracy in representing thermodynamic and structural properties with respect to the underlying atomistic representation of the system. These results are compared to those obtained following the principles of statistical mechanics that start from the same underlying atomistic system. To this end, the potentials of mean force for lateral association in Martini 3.0 binary lipid bilayers are decomposed into their entropic and enthalpic components and compared to those of corresponding atomistic bilayers that have been projected onto equivalent coarse-grained mappings but evolved under the fully atomistic forces. This is accomplished by applying the reversible work theorem to lateral pair correlation functions between coarse-grained lipid beads taken at a range of different temperatures. The entropy-enthalpy decompositions provide a metric by which the underlying statistical mechanical properties of Martini can be investigated. Overall, Martini 3.0 is found to fail to properly partition entropy and enthalpy for the PMFs compared to the mapped all-atom results, despite changes made to the force field from the Martini 2.0 version. This outcome points to the fact that the development of more accurate top-down coarse-grained models such as Martini will likely necessitate temperature-dependent terms in the corresponding CG force-field; although necessary, this may not be sufficient to improve Martini. In addition to the entropy-enthalpy decompositions, Martini 3.0 produces an incorrect undulation spectrum, in particular at intermediate length scales of biophysical pertinence.

Understanding dynamics in coarse-grained models. IV. Connection of fine-grained and coarse-grained dynamics with the Stokes-Einstein and Stokes-Einstein-Debye relations

Applying an excess entropy scaling formalism to the coarse-grained (CG) dynamics of liquids, we discovered that missing rotational motions during the CG process are responsible for artificially accelerated CG dynamics. In the context of the dynamic representability between the fine-grained (FG) and CG dynamics, this work introduces the well-known Stokes-Einstein and Stokes-Einstein-Debye relations to unravel the rotational dynamics underlying FG trajectories, thereby allowing for an indirect evaluation of the effective rotations based only on the translational information at the reduced CG resolution. Since the representability issue in CG modeling limits a direct evaluation of the shear stress appearing in the Stokes-Einstein and Stokes-Einstein-Debye relations, we introduce a translational relaxation time as a proxy to employ these relations, and we demonstrate that these relations hold for the ambient conditions studied in our series of work. Additional theoretical links to our previous work are also established. First, we demonstrate that the effective hard sphere radius determined by the classical perturbation theory can approximate the complex hydrodynamic radius value reasonably well. Furthermore, we present a simple derivation of an excess entropy scaling relationship for viscosity by estimating the elliptical integral of molecules. In turn, since the translational and rotational motions at the FG level are correlated to each other, we conclude that the "entropy-free" CG diffusion only depends on the shape of the reference molecule. Our results and analyses impart an alternative way of recovering the FG diffusion from the CG description by coupling the translational and rotational motions at the hydrodynamic level.

Rigorous Progress in Coarse-Graining

Low-resolution coarse-grained (CG) models provide remarkable computational and conceptual advantages for simulating soft materials. In principle, bottom-up CG models can reproduce all structural and thermodynamic properties of atomically detailed models that can be observed at the resolution of the CG model. This review discusses recent progress in developing theory and computational methods for achieving this promise. We first briefly review variational approaches for parameterizing interaction potentials and their relationship to machine learning methods. We then discuss recent approaches for simultaneously improving both the transferability and thermodynamic properties of bottom-up models by rigorously addressing the density and temperature dependence of these potentials. We also briefly discuss exciting progress in modeling high-resolution observables with low-resolution CG models. More generally, we highlight the essential role of the bottom-up framework not only for fundamentally understanding the limitations of prior CG models but also for developing robust computational methods that resolve these limitations in practice.

Hierarchical Framework for Predicting Entropies in Bottom-Up Coarse-Grained Models

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.

Chemically Specific Systematic Coarse-Grained Polymer Model with Both Consistently Structural and Dynamical Properties

The coarse-grained (CG) model serves as a powerful tool for the simulation of polymer systems; its reliability depends on the accurate representation of both structural and dynamical properties. However, strong correlations between structural and dynamical properties on different scales and also a strong memory effect, enforced by chain connectivity between monomers in polymer systems, render developing a chemically specific systematic CG model a formidable task. In this study, we report a systematic CG approach that combines the iterative Boltzmann inversion (IBI) method and the generalized Langevin equation (GLE) dynamics. Structural properties are ensured by using conservative CG potentials derived from the IBI method. To retrieve the correct dynamical properties in the system, we demonstrate that using a combination of a Rouse-type delta function and a time-dependent short-time kernel in the GLE simulation is practically efficient. The former can be used to adjust the long-time diffusion dynamics, and the latter can be reconstructed from an iterative procedure according to the velocity autocorrelation function (ACF) from all-atomistic (AA) simulations. Taking the polystyrene as an example, we show that not only structural properties of radial distribution function, intramolecular bond, and angle distributions can be reproduced but also dynamical properties of mean-square displacement, velocity ACF, and force ACF resulted from our CG model have quantitative agreement with the reference AA model. In addition, reasonable agreements are observed in other collective properties between our GLE-CG model and the AA simulations as well.

Maximum Entropy Theory of Multiscale Coarse-Graining via Matching Thermodynamic Forces: Application to a Molecular Crystal (TATB)

The MSCG/FM (multiscale coarse-graining via force-matching) approach is an efficient supervised machine learning method to develop microscopically informed coarse-grained (CG) models. We present a theory based on the principle of maximum entropy (PME) enveloping the existing MSCG/FM approaches. This theory views the MSCG/FM method as a special case of matching the thermodynamic forces from the extended ensemble described by the set of thermodynamic (relevant) system coordinates. This set may include CG coordinates, the stress tensor, applied external fields, and so forth, and may be characterized by nonequilibrium conditions. Following the presentation of the theory, we discuss the consistent matching of both bonded and nonbonded interactions. The proposed PME formulation is used as a starting point to extend the MSCG/FM method to the constant strain ensemble, which together with the explicit matching of the bonded forces is better suited for coarse-graining anisotropic media at a submolecular resolution. The theory is demonstrated by performing the fine coarse-graining of crystalline 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), a well-known insensitive molecular energetic material, which exhibits highly anisotropic mechanical properties.

Improved Protein Dynamics and Hydration in the Martini3 Coarse-Grain Model

The Martini coarse-grain force-field has emerged as an important framework to probe cellular processes at experimentally relevant time- and length-scales. However, the recently developed version, the Martini3 force-field with the implemented Go̅ model (Martini3Go̅), as well as previous variants of the Martini model have not been benchmarked and rigorously tested for globular proteins. In this study, we consider three globular proteins, ubiquitin, lysozyme, and cofilin, and compare protein dynamics and hydration with observables from experiments and all-atom simulations. We show that the Martini3Go̅ model is able to accurately model the structural and dynamic features of small globular proteins. Overall, the structural integrity of the proteins is maintained, as validated by contact maps, radii of gyration (Rg), and SAXS profiles. The chemical shifts predicted from the ensemble sampled in the simulations are consistent with the experimental data. Further, a good match is observed in the protein-water interaction energetics, and the hydration levels of the residues are similar to atomistic simulations. However, the protein-water interaction dynamics is not accurately represented and appears to depend on the protein structural complexity, residue specificity, and water dynamics. Our work is a step toward testing and assessing the Martini3Go̅ model and provides insights into future efforts to refine Martini models with improved solvation effects and better correspondence to the underlying all-atom systems.

Insight into the Density-Dependence of Pair Potentials for Predictive Coarse-Grained Models

We investigate the temperature- and density-dependence of effective pair potentials for 1-site coarse-grained (CG) models of two industrial solvents, 1,4-dioxane and tetrahydrofuran. We observe that the calculated pair potentials are much more sensitive to density than to temperature. The generalized-Yvon-Born-Green framework reveals that this striking density-dependence reflects corresponding variations in the many-body correlations that determine the environment-mediated indirect contribution to the pair mean force. Moreover, we demonstrate, perhaps surprisingly, that this density-dependence is not important for accurately modeling the intermolecular structure. Accordingly, we adopt a density-independent interaction potential and transfer the density-dependence of the calculated pair potentials into a configuration-independent volume potential. Furthermore, we develop a single global potential that accurately models the intermolecular structure and pressure-volume equation of state across a very wide range of liquid state points. Consequently, this work provides fundamental insight into the density-dependence of effective pair potentials and also provides a significant step toward developing predictive CG models for efficiently modeling industrial solvents.

Perturbative Expansion in Reciprocal Space: Bridging Microscopic and Mesoscopic Descriptions of Molecular Interactions

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.

Gaussian representation of coarse-grained interactions of liquids: Theory, parametrization, and transferability

Coarse-grained (CG) interactions determined via bottom-up methodologies can faithfully reproduce the structural correlations observed in fine-grained (atomistic resolution) systems, yet they can suffer from limited extensibility due to complex many-body correlations. As part of an ongoing effort to understand and improve the applicability of bottom-up CG models, we propose an alternative approach to address both accuracy and transferability. Our main idea draws from classical perturbation theory to partition the hard sphere repulsive term from effective CG interactions. We then introduce Gaussian basis functions corresponding to the system's characteristic length by linking these Gaussian sub-interactions to the local particle densities at each coordination shell. The remaining perturbative long-range interaction can be treated as a collective solvation interaction, which we show exhibits a Gaussian form derived from integral equation theories. By applying this numerical parametrization protocol to CG liquid systems, our microscopic theory elucidates the emergence of Gaussian interactions in common phenomenological CG models. To facilitate transferability for these reduced descriptions, we further infer equations of state to determine the sub-interaction parameter as a function of the system variables. The reduced models exhibit excellent transferability across the thermodynamic state points. Furthermore, we propose a new strategy to design the cross-interactions between distinct CG sites in liquid mixtures. This involves combining each Gaussian in the proper radial domain, yielding accurate CG potentials of mean force and structural correlations for multi-component systems. Overall, our findings establish a solid foundation for constructing transferable bottom-up CG models of liquids with enhanced extensibility.

Understanding dynamics in coarse-grained models. III. Roles of rotational motion and translation-rotation coupling in coarse-grained dynamics

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.

Pragmatic Coarse-Graining of Proteins: Models and Applications

The molecular details involved in the folding, dynamics, organization, and interaction of proteins with other molecules are often difficult to assess by experimental techniques. Consequently, computational models play an ever-increasing role in the field. However, biological processes involving large-scale protein assemblies or long time scale dynamics are still computationally expensive to study in atomistic detail. For these applications, employing coarse-grained (CG) modeling approaches has become a key strategy. In this Review, we provide an overview of what we call pragmatic CG protein models, which are strategies combining, at least in part, a physics-based implementation and a top-down experimental approach to their parametrization. In particular, we focus on CG models in which most protein residues are represented by at least two beads, allowing these models to retain some degree of chemical specificity. A description of the main modern pragmatic protein CG models is provided, including a review of the most recent applications and an outlook on future perspectives in the field.

OpenMSCG: A Software Tool for Bottom-Up Coarse-Graining

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.

Machine learned coarse-grained protein force-fields: Are we there yet?

The successful recent application of machine learning methods to scientific problems includes the learning of flexible and accurate atomic-level force-fields for materials and biomolecules from quantum chemical data. In parallel, the machine learning of force-fields at coarser resolutions is rapidly gaining relevance as an efficient way to represent the higher-body interactions needed in coarse-grained force-fields to compensate for the omitted degrees of freedom. Coarse-grained models are important for the study of systems at time and length scales exceeding those of atomistic simulations. However, the development of transferable coarse-grained models via machine learning still presents significant challenges. Here, we discuss recent developments in this field and current efforts to address the remaining challenges.

Slicing and Dicing: Optimal Coarse-Grained Representation to Preserve Molecular Kinetics

The aim of molecular coarse-graining approaches is to recover relevant physical properties of the molecular system via a lower-resolution model that can be more efficiently simulated. Ideally, the lower resolution still accounts for the degrees of freedom necessary to recover the correct physical behavior. The selection of these degrees of freedom has often relied on the scientist's chemical and physical intuition. In this article, we make the argument that in soft matter contexts desirable coarse-grained models accurately reproduce the long-time dynamics of a system by correctly capturing the rare-event transitions. We propose a bottom-up coarse-graining scheme that correctly preserves the relevant slow degrees of freedom, and we test this idea for three systems of increasing complexity. We show that in contrast to this method existing coarse-graining schemes such as those from information theory or structure-based approaches are not able to recapitulate the slow time scales of the system.

Statistical Mechanical Design Principles for Coarse-Grained Interactions across Different Conformational Free Energy Surfaces

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.

Flow-Matching: Efficient Coarse-Graining of Molecular Dynamics without Forces

Coarse-grained (CG) molecular simulations have become a standard tool to study molecular processes on time and length scales inaccessible to all-atom simulations. Parametrizing CG force fields to match all-atom simulations has mainly relied on force-matching or relative entropy minimization, which require many samples from costly simulations with all-atom or CG resolutions, respectively. Here we present flow-matching, a new training method for CG force fields that combines the advantages of both methods by leveraging normalizing flows, a generative deep learning method. Flow-matching first trains a normalizing flow to represent the CG probability density, which is equivalent to minimizing the relative entropy without requiring iterative CG simulations. Subsequently, the flow generates samples and forces according to the learned distribution in order to train the desired CG free energy model via force-matching. Even without requiring forces from the all-atom simulations, flow-matching outperforms classical force-matching by an order of magnitude in terms of data efficiency and produces CG models that can capture the folding and unfolding transitions of small proteins.

Understanding dynamics in coarse-grained models. II. Coarse-grained diffusion modeled using hard sphere theory

The first paper of this series [J. Chem. Phys. 158, 034103 (2023)] demonstrated that excess entropy scaling holds for both fine-grained and corresponding coarse-grained (CG) systems. Despite its universality, a more exact determination of the scaling relationship was not possible due to the semi-empirical nature. In this second paper, an analytical excess entropy scaling relation is derived for bottom-up CG systems. At the single-site CG resolution, effective hard sphere systems are constructed that yield near-identical dynamical properties as the target CG systems by taking advantage of how hard sphere dynamics and excess entropy can be analytically expressed in terms of the liquid packing fraction. Inspired by classical equilibrium perturbation theories and recent advances in constructing hard sphere models for predicting activated dynamics of supercooled liquids, we propose a new approach for understanding the diffusion of molecular liquids in the normal regime using hard sphere reference fluids. The proposed "fluctuation matching" is designed to have the same amplitude of long wavelength density fluctuations (dimensionless compressibility) as the CG system. Utilizing the Enskog theory to derive an expression for hard sphere diffusion coefficients, a bridge between the CG dynamics and excess entropy is then established. The CG diffusion coefficient can be roughly estimated using various equations of the state, and an accurate prediction of accelerated CG dynamics at different temperatures is also possible in advance of running any CG simulation. By introducing another layer of coarsening, these findings provide a more rigorous method to assess excess entropy scaling and understand the accelerated CG dynamics of molecular fluids.

Understanding dynamics in coarse-grained models. I. Universal excess entropy scaling relationship

Coarse-grained (CG) models facilitate an efficient exploration of complex systems by reducing the unnecessary degrees of freedom of the fine-grained (FG) system while recapitulating major structural correlations. Unlike structural properties, assessing dynamic properties in CG modeling is often unfeasible due to the accelerated dynamics of the CG models, which allows for more efficient structural sampling. Therefore, the ultimate goal of the present series of articles is to establish a better correspondence between the FG and CG dynamics. To assess and compare dynamical properties in the FG and the corresponding CG models, we utilize the excess entropy scaling relationship. For Paper I of this series, we provide evidence that the FG and the corresponding CG counterpart follow the same universal scaling relationship. By carefully reviewing and examining the literature, we develop a new theory to calculate excess entropies for the FG and CG systems while accounting for entropy representability. We demonstrate that the excess entropy scaling idea can be readily applied to liquid water and methanol systems at both the FG and CG resolutions. For both liquids, we reveal that the scaling exponents remain unchanged from the coarse-graining process, indicating that the scaling behavior is universal for the same underlying molecular systems. Combining this finding with the concept of mapping entropy in CG models, we show that the missing entropy plays an important role in accelerating the CG dynamics.

Bottom-up Coarse-Graining: Principles and Perspectives

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.

Consistent and Transferable Force Fields for Statistical Copolymer Systems at the Mesoscale

The statistical trajectory matching (STM) method was applied successfully to derive coarse grain (CG) models for bulk properties of homopolymers. The extension of the methodology for building CG models for statistical copolymer systems is much more challenging. We present here the strategy for developing CG models for styrene-butadiene-rubber, and we compare the quality of the resulting CG force fields on the structure and thermodynamics at different chemical compositions. The CG models are used through the use of a genuine mesoscopic method called the dissipative particle dynamics method and compared to high-resolution molecular dynamics simulations. We conclude that the STM method is able to produce coarse-grained potentials that are transferable in composition by using only a few reference systems. Additionally, this methodology can be applied on any copolymer system.

Coarse-Grained Water Model Development for Accurate Dynamics and Structure Prediction

Several coarse-graining (CG) methods have been combined to develop a CG model of water capable of the accurate prediction of structure and dynamics properties. The multiscale coarse-graining (MS-CG) method based on force matching and the PDF-based coarse-graining method were used for accurate dynamics prediction. The iterative Boltzmann inversion (IBI) method was added for accurate structure representation. The approach is applied to bulk water, and the results show close reproduction of the CG structure when compared with the reference atomistic data. The combination of MS-CG and IBI methods facilitates the development of CG force fields at different temperatures based on a single MS-CG coarse-graining procedure. The dynamic properties of the CG water model closely match those obtained from the reference atomistic system. The general application of this approach to any existing coarse-graining methods is discussed.

Insights into the Phase Diagram of Pluronic L64 Using Coarse-Grained Molecular Dynamics Simulations

We investigate the concentration-dependent phase diagram of pluronic L64 in aqueous media at 300 and 320 K using coarse-grained (CG) molecular dynamics (MD) simulations. The CG model is derived by adapting the Martini model for nonbonded interactions along with the Boltzmann inversion (BI) of bonded interactions from all-atom (AA) simulations. Our derived CG model successfully captures the experimentally observed micellar-, hexagonal-, lamellar-, and polymer-rich solution phase. The end-to-end distance reveals the conformational change from an open-chain structure in the micellar phase to a folded-chain structure in the lamellar phase, increasing the orientational order. An increase in temperature leads to expulsion of water molecules from the L64 moiety, suggesting an increase in L64 hydrophobicity. Thermodynamic analysis using the two-phase thermodynamics (2PT) method suggests the entropy of the system decreases with increasing L64 concentration and the decrease in free energy (F) with temperature is mainly driven by the entropic factor (-TS). Further, the increase in aggregation number at lower concentrations and self-assembly at very high concentrations is energetically driven, whereas the change from the micellar phase to the lamellar phase with increasing L64 concentration is entropically driven. Our model provides molecular insights into L64 phases which can be further explored to design functionality-based suprastructures for drug delivery and tissue engineering applications.

Inositol Hexakisphosphate (IP6) Accelerates Immature HIV-1 Gag Protein Assembly toward Kinetically Trapped Morphologies

During the late stages of the HIV-1 lifecycle, immature virions are produced by the concerted activity of Gag polyproteins, primarily mediated by the capsid (CA) and spacer peptide 1 (SP1) domains, which assemble into a spherical lattice, package viral genomic RNA, and deform the plasma membrane. Recently, inositol hexakisphosphate (IP6) has been identified as an essential assembly cofactor that efficiently produces both immature virions in vivo and immature virus-like particles in vitro. To date, however, several distinct mechanistic roles for IP6 have been proposed on the basis of independent functional, structural, and kinetic studies. In this work, we investigate the molecular influence of IP6 on the structural outcomes and dynamics of CA/SP1 assembly using coarse-grained (CG) molecular dynamics (MD) simulations and free energy calculations. Here, we derive a bottom-up, low-resolution, and implicit-solvent CG model of CA/SP1 and IP6, and simulate their assembly under conditions that emulate both in vitro and in vivo systems. Our analysis identifies IP6 as an assembly accelerant that promotes curvature generation and fissure-like defects throughout the lattice. Our findings suggest that IP6 induces kinetically trapped immature morphologies, which may be physiologically important for later stages of viral morphogenesis and potentially useful for virus-like particle technologies.

Accurate model of liquid-liquid phase behavior of intrinsically disordered proteins from optimization of single-chain properties

Many intrinsically disordered proteins (IDPs) may undergo liquid-liquid phase separation (LLPS) and participate in the formation of membraneless organelles in the cell, thereby contributing to the regulation and compartmentalization of intracellular biochemical reactions. The phase behavior of IDPs is sequence dependent, and its investigation through molecular simulations requires protein models that combine computational efficiency with an accurate description of intramolecular and intermolecular interactions. We developed a general coarse-grained model of IDPs, with residue-level detail, based on an extensive set of experimental data on single-chain properties. Ensemble-averaged experimental observables are predicted from molecular simulations, and a data-driven parameter-learning procedure is used to identify the residue-specific model parameters that minimize the discrepancy between predictions and experiments. The model accurately reproduces the experimentally observed conformational propensities of a set of IDPs. Through two-body as well as large-scale molecular simulations, we show that the optimization of the intramolecular interactions results in improved predictions of protein self-association and LLPS.

Accurate model of liquid-liquid phase behavior of intrinsically disordered proteins from optimization of single-chain properties

Many intrinsically disordered proteins (IDPs) may undergo liquid-liquid phase separation (LLPS) and participate in the formation of membraneless organelles in the cell, thereby contributing to the regulation and compartmentalization of intracellular biochemical reactions. The phase behavior of IDPs is sequence dependent, and its investigation through molecular simulations requires protein models that combine computational efficiency with an accurate description of intramolecular and intermolecular interactions. We developed a general coarse-grained model of IDPs, with residue-level detail, based on an extensive set of experimental data on single-chain properties. Ensemble-averaged experimental observables are predicted from molecular simulations, and a data-driven parameter-learning procedure is used to identify the residue-specific model parameters that minimize the discrepancy between predictions and experiments. The model accurately reproduces the experimentally observed conformational propensities of a set of IDPs. Through two-body as well as large-scale molecular simulations, we show that the optimization of the intramolecular interactions results in improved predictions of protein self-association and LLPS.

A microcanonical approach to temperature-transferable coarse-grained models using the relative entropy

Bottom-up coarse-graining methods provide systematic tools for creating simplified models of molecular systems. However, coarse-grained (CG) models produced with such methods frequently fail to accurately reproduce all thermodynamic properties of the reference atomistic systems they seek to model and, moreover, can fail in even more significant ways when used at thermodynamic state points different from the reference conditions. These related problems of representability and transferability limit the usefulness of CG models, especially those of strongly state-dependent systems. In this work, we present a new strategy for creating temperature-transferable CG models using a single reference system and temperature. The approach is based on two complementary concepts. First, we switch to a microcanonical basis for formulating CG models, focusing on effective entropy functions rather than energy functions. This allows CG models to naturally represent information about underlying atomistic energy fluctuations, which would otherwise be lost. Such information not only reproduces energy distributions of the reference model but also successfully predicts the correct temperature dependence of the CG interactions, enabling temperature transferability. Second, we show that relative entropy minimization provides a direct and systematic approach to parameterize such classes of temperature-transferable CG models. We calibrate the approach initially using idealized model systems and then demonstrate its ability to create temperature-transferable CG models for several complex molecular liquids.

Dynamical properties across different coarse-grained models for ionic liquids

Room-temperature ionic liquids (RTILs) stand out among molecular liquids for their rich physicochemical characteristics, including structural and dynamic heterogeneity. The significance of electrostatic interactions in RTILs results in long characteristic length- and timescales, and has motivated the development of a number of coarse-grained (CG) simulation models. In this study, we aim to better understand the connection between certain CG parameterization strategies and the dynamical properties and transferability of the resulting models. We systematically compare five CG models: a model largely parameterized from experimental thermodynamic observables; a refinement of this model to increase its structural accuracy; and three models that reproduce a given set of structural distribution functions by construction, with varying intramolecular parameterizations and reference temperatures. All five CG models display limited structural transferability over temperature, and also result in various effective dynamical speedup factors, relative to a reference atomistic model. On the other hand, the structure-based CG models tend to result in more consistent cation-anion relative diffusion than the thermodynamic-based models, for a single thermodynamic state point. By linking short- and long-timescale dynamical behaviors, we demonstrate that the varying dynamical properties of the different CG models can be largely collapsed onto a single curve, which provides evidence for a route to constructing dynamically-consistent CG models of RTILs.

A Deep Graph Network-Enhanced Sampling Approach to Efficiently Explore the Space of Reduced Representations of Proteins

The limits of molecular dynamics (MD) simulations of macromolecules are steadily pushed forward by the relentless development of computer architectures and algorithms. The consequent explosion in the number and extent of MD trajectories induces the need for automated methods to rationalize the raw data and make quantitative sense of them. Recently, an algorithmic approach was introduced by some of us to identify the subset of a protein's atoms, or mapping, that enables the most informative description of the system. This method relies on the computation, for a given reduced representation, of the associated mapping entropy, that is, a measure of the information loss due to such simplification; albeit relatively straightforward, this calculation can be time-consuming. Here, we describe the implementation of a deep learning approach aimed at accelerating the calculation of the mapping entropy. We rely on Deep Graph Networks, which provide extreme flexibility in handling structured input data and whose predictions prove to be accurate and-remarkably efficient. The trained network produces a speedup factor as large as 105 with respect to the algorithmic computation of the mapping entropy, enabling the reconstruction of its landscape by means of the Wang-Landau sampling scheme. Applications of this method reach much further than this, as the proposed pipeline is easily transferable to the computation of arbitrary properties of a molecular structure.

Constructing many-body dissipative particle dynamics models of fluids from bottom-up coarse-graining

Since their emergence in the 1990s, mesoscopic models of fluids have been widely used to study complex organization and transport phenomena beyond the molecular scale. Even though these models are designed based on results from physics at the meso- and macroscale, such as fluid mechanics and statistical field theory, the underlying microscopic foundation of these models is not as well defined. This paper aims to build such a systematic connection using bottom-up coarse-graining methods. From the recently developed dynamic coarse-graining scheme, we introduce a statistical inference framework of explicit many-body conservative interaction that quantitatively recapitulates the mesoscopic structure of the underlying fluid. To further consider the dissipative and fluctuation forces, we design a novel algorithm that parameterizes these forces. By utilizing this algorithm, we derive pairwise decomposable friction kernels under both non-Markovian and Markovian limits where both short- and long-time features of the coarse-grained dynamics are reproduced. Finally, through these new developments, the many-body dissipative particle dynamics type of equations of motion are successfully derived. The methodologies developed in this work thus open a new avenue for the construction of direct bottom-up mesoscopic models that naturally bridge the meso- and macroscopic physics.

Investigating the energetic and entropic components of effective potentials across a glass transition

By eliminating unnecessary details, coarse-grained (CG) models provide the necessary efficiency for simulating scales that are inaccessible to higher resolution models. However, because they average over atomic details, the effective potentials governing CG degrees of freedom necessarily incorporate significant entropic contributions, which limit their transferability and complicate the treatment of thermodynamic properties. This work employs a dual-potential approach to consider the energetic and entropic contributions to effective interaction potentials for CG models. Specifically, we consider one- and three-site CG models for ortho-terphenyl (OTP) both above and below its glass transition. We employ the multiscale coarse-graining (MS-CG) variational principle to determine interaction potentials that accurately reproduce the structural properties of an all-atom (AA) model for OTP at each state point. We employ an energy-matching variational principle to determine an energy operator that accurately reproduces the intra- and inter-molecular energy of the AA model. While the MS-CG pair potentials are almost purely repulsive, the corresponding pair energy functions feature a pronounced minima that corresponds to contacting benzene rings. These energetic functions then determine an estimate for the entropic component of the MS-CG interaction potentials. These entropic functions accurately predict the MS-CG pair potentials across a wide range of liquid state points at constant density. Moreover, the entropic functions also predict pair potentials that quite accurately model the AA pair structure below the glass transition. Thus, the dual-potential approach appears a promising approach for modeling AA energetics, as well as for predicting the temperature-dependence of CG effective potentials.

Implicit Solvation Using the Superposition Approximation (IS-SPA): Extension to Peptides in a Polar Solvent

Efficient, accurate, and adaptable implicit solvent models remain a significant challenge in the field of molecular simulation. A recent implicit solvent model, IS-SPA, based on approximating the mean solvent force using the superposition approximation, provides a platform to achieve these goals. IS-SPA was originally developed to handle nonpolar solutes in a polar solvent and did not accurately capture polar solvation. Here, we demonstrate that IS-SPA can accurately capture polar solvation by incorporating solvent orientation and accounting for the contributions from long ranged electrostatics. Solvent orientation is approximated as that of an ideal dipole aligned in a mean electrostatic field and an analytic form of the long ranged electrostatics is derived. Parameters for the model are calculated from explicit solvent simulations of an isolated atom or molecule and include atom-based solvent densities, mean electric field functions, radially symmetric averaged Lennard-Jones forces, and multipoles of the explicit solvent model. Using these parameters, IS-SPA accounts for asymmetry of charge solvation and reproduces the explicit solvent potential of mean force of dimerization of two oppositely charged Lennard-Jones spheres in chloroform with high fidelity. Additionally, the model more accurately captures the effect of explicit solvent on the monomer and dimer configurations of alanine dipeptide in chloroform than a generalized Born or constant density dielectric model. The current version of the algorithm is expected to outperform explicit solvent simulations for aggregation of small peptides at concentrations below 150 mM, well above the typical experimental concentrations for these materials.

A new one-site coarse-grained model for water: Bottom-up many-body projected water (BUMPer). II. Temperature transferability and structural properties at low temperature

A number of studies have constructed coarse-grained (CG) models of water to understand its anomalous properties. Most of these properties emerge at low temperatures, and an accurate CG model needs to be applicable to these low-temperature ranges. However, direct use of CG models parameterized from other temperatures, e.g., room temperature, encounters a problem known as transferability, as the CG potential essentially follows the form of the many-body CG free energy function. Therefore, temperature-dependent changes to CG interactions must be accounted for. The collective behavior of water at low temperature is generally a many-body process, which often motivates the use of expensive many-body terms in the CG interactions. To surmount the aforementioned problems, we apply the Bottom-Up Many-Body Projected Water (BUMPer) CG model constructed from Paper I to study the low-temperature behavior of water. We report for the first time that the embedded three-body interaction enables BUMPer, despite its pairwise form, to capture the growth of ice at the ice/water interface with corroborating many-body correlations during the crystal growth. Furthermore, we propose temperature transferable BUMPer models that are indirectly constructed from the free energy decomposition scheme. Changes in CG interactions and corresponding structures are faithfully recapitulated by this framework. We further extend BUMPer to examine its ability to predict the structure, density, and diffusion anomalies by employing an alternative analysis based on structural correlations and pairwise potential forms to predict such anomalies. The presented analysis highlights the existence of these anomalies in the low-temperature regime and overcomes potential transferability problems.

Coarse-Grained Force Fields from the Perspective of Statistical Mechanics: Better Understanding of the Origins of a MARTINI Hangover

The popular MARTINI coarse-grained model is used as a test case to analyze the adherence of top-down coarse-grained molecular dynamics models (i.e., models primarily parametrized to match experimental results) to the known features of statistical mechanics for the underlying all-atom representations. Specifically, the temperature dependence of various pair distribution functions, and hence their underlying potentials of mean force via the reversible work theorem, are compared between MARTINI 2.0, Dry MARTINI, and all-atom simulations mapped onto equivalent coarse-grained sites for certain lipid bilayers. It is found that the MARTINI models do not completely capture the lipid structure seen in atomistic simulations as projected onto the coarse-grained mappings and that issues of accuracy and temperature transferability arise due to an incorrect enthalpy-entropy decomposition of these potentials of mean force. The potential of mean force for the association of two amphipathic helices in a lipid bilayer is also calculated, and especially at shorter ranges, the MARTINI and all-atom projection results differ substantially. The former is much less repulsive and hence will lead to a higher probability of MARTINI helix association in the MARTINI bilayer than occurs in the actual all-atom case. Additionally, the bilayer height fluctuation spectra are calculated for the MARTINI model, and compared to the all-atom results, it is found that the magnitude of thermally averaged amplitudes at intermediate length scales are quite different, pointing to a number of possible consequences for realistic modeling of membrane processes. Taken as a whole, the results presented here show disagreement in the enthalpic and entropic driving forces driving lateral structure in lipid bilayers as well as quantitative differences in association of embedded amphipathic helices, which can help direct future efforts to parametrize CG models with better agreement to the all-atom systems they aspire to represent.

Swarm-CG: Automatic Parametrization of Bonded Terms in MARTINI-Based Coarse-Grained Models of Simple to Complex Molecules via Fuzzy Self-Tuning Particle Swarm Optimization

We present Swarm-CG, a versatile software for the automatic iterative parametrization of bonded parameters in coarse-grained (CG) models, ideal in combination with popular CG force fields such as MARTINI. By coupling fuzzy self-tuning particle swarm optimization to Boltzmann inversion, Swarm-CG performs accurate bottom-up parametrization of bonded terms in CG models composed of up to 200 pseudo atoms within 4-24 h on standard desktop machines, using default settings. The software benefits from a user-friendly interface and two different usage modes (default and advanced). We particularly expect Swarm-CG to support and facilitate the development of new CG models for the study of complex molecular systems interesting for bio- and nanotechnology. Excellent performances are demonstrated using a benchmark of 9 molecules of diverse nature, structural complexity, and size. Swarm-CG is available with all its dependencies via the Python Package Index (PIP package: swarm-cg). Demonstration data are available at: www.github.com/GMPavanLab/SwarmCG.

An Information-Theory-Based Approach for Optimal Model Reduction of Biomolecules

In theoretical modeling of a physical system, a crucial step consists of the identification of those degrees of freedom that enable a synthetic yet informative representation of it. While in some cases this selection can be carried out on the basis of intuition and experience, straightforward discrimination of the important features from the negligible ones is difficult for many complex systems, most notably heteropolymers and large biomolecules. We here present a thermodynamics-based theoretical framework to gauge the effectiveness of a given simplified representation by measuring its information content. We employ this method to identify those reduced descriptions of proteins, in terms of a subset of their atoms, that retain the largest amount of information from the original model; we show that these highly informative representations share common features that are intrinsically related to the biological properties of the proteins under examination, thereby establishing a bridge between protein structure, energetics, and function.

Temperature and Phase Transferable Bottom-up Coarse-Grained Models

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.

Evidence of information limitations in coarse-grained models

Developing accurate coarse-grained (CG) models is critical for addressing long time and length scale phenomena with molecular simulations. Here, we distinguish and quantify two sources of error that are relevant to CG models in order to guide further methods development: "representability" errors, which result from the finite basis associated with the chosen functional form of the CG model and mapping operator, and "information" errors, which result from the limited kind and quantity of data supplied to the CG parameterization algorithm. We have performed a systematic investigation of these errors by generating all possible CG models of three liquids (butane, 1-butanol, and 1,3-propanediol) that conserve a set of chemically motivated locality and topology relationships. In turn, standard algorithms (iterative Boltzmann inversion, IBI, and multiscale coarse-graining, MSCG) were used to parameterize the models and the CG predictions were compared with atomistic results. For off-target properties, we observe a strong correlation between the accuracy and the resolution of the CG model, which suggests that the approximations represented by MSCG and IBI deteriorate with decreasing resolution. Conversely, on-target properties exhibit an extremely weak resolution dependence that suggests a limited role of representability errors in model accuracy. Taken together, these results suggest that simple CG models are capable of utilizing more information than is provided by standard parameterization algorithms, and that model accuracy can be improved by algorithm development rather than resorting to more complicated CG models.

Dual-potential approach for coarse-grained implicit solvent models with accurate, internally consistent energetics and predictive transferability

The dual-potential approach promises coarse-grained (CG) models that accurately reproduce both structural and energetic properties, while simultaneously providing predictive estimates for the temperature-dependence of the effective CG potentials. In this work, we examine the dual-potential approach for implicit solvent CG models that reflect large entropic effects from the eliminated solvent. Specifically, we construct implicit solvent models at various resolutions, R, by retaining a fraction 0.10 ≤ R ≤ 0.95 of the molecules from a simple fluid of Lennard-Jones spheres. We consider the dual-potential approach in both the constant volume and constant pressure ensembles across a relatively wide range of temperatures. We approximate the many-body potential of mean force for the remaining solutes with pair and volume potentials, which we determine via multiscale coarse-graining and self-consistent pressure-matching, respectively. Interestingly, with increasing temperature, the pair potentials appear increasingly attractive, while the volume potentials become increasingly repulsive. The dual-potential approach not only reproduces the atomic energetics but also quite accurately predicts this temperature-dependence. We also derive an exact relationship between the thermodynamic specific heat of an atomic model and the energetic fluctuations that are observable at the CG resolution. With this generalized fluctuation relationship, the approximate CG models quite accurately reproduce the thermodynamic specific heat of the underlying atomic model.