The multiscale coarse-graining method. IX. A general method for construction of three body coarse-grained force fields

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.

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.

Representability and Dynamical Consistency in Coarse-Grained Models

We address the challenge of representativity and dynamical consistency when unbonded fine-grained particles are collected together into coarse-grained particles. We implement a hybrid procedure for identifying and tracking the underlying fine-grained particles─e.g., atoms or molecules─by exchanging them between the coarse-grained particles periodically at a characteristic time. The exchange involves a back-mapping of the coarse-grained particles into fine-grained particles and a subsequent reassignment to coarse-grained particles conserving total mass and momentum. We find that an appropriate choice of the characteristic exchange time can lead to the correct effective diffusion rate of the fine-grained particles when simulated in hybrid coarse-grained dynamics. In the compressed (supercritical) fluid regime, without the exchange term, fine-grained particles remain associated with a given coarse-grained particle, leading to substantially lower diffusion rates than seen in all-atom molecular dynamics of the fine-grained particles. Thus, this work confirms the need for addressing the representativity of fine-grained particles within coarse-grained particles and offers a simple exchange mechanism so as to retain dynamical consistency between the fine- and coarse-grained scales.

Formalizing Coarse-Grained Representations of Anisotropic Interactions at Multimeric Protein Interfaces Using Virtual Sites

Molecular simulations of biomacromolecules that assemble into multimeric complexes remain a challenge due to computationally inaccessible length and time scales. Low-resolution and implicit-solvent coarse-grained modeling approaches using traditional nonbonded interactions (both pairwise and spherically isotropic) have been able to partially address this gap. However, these models may fail to capture the complex anisotropic interactions present at macromolecular interfaces unless higher-order interaction potentials are incorporated at the expense of the computational cost. In this work, we introduce an alternate and systematic approach to represent directional interactions at protein-protein interfaces by using virtual sites restricted to pairwise interactions. We show that virtual site interaction parameters can be optimized within a relative entropy minimization framework by using only information from known statistics between coarse-grained sites. We compare our virtual site models to traditional coarse-grained models using two case studies of multimeric protein assemblies and find that the virtual site models predict pairwise correlations with higher fidelity and, more importantly, assembly behavior that is morphologically consistent with experiments. Our study underscores the importance of anisotropic interaction representations and paves the way for more accurate yet computationally efficient coarse-grained simulations of macromolecular assembly in future research.

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.

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.

Hierarchical Machine Learning of Low-Resolution Coarse-Grained Free Energy Potentials

A force-matching-based method for supervised machine learning (ML) of coarse-grained (CG) free energy (FE) potentials─known as multiscale coarse-graining via force-matching (MSCG/FM)─is an efficient method to develop microscopically informed CG models that are thermodynamically and statistically equivalent to the reference microscopic models. For low-resolution models, when the coarse-graining is at supramolecular scales, objective-oriented clustering of nonbonded particles is required and the reduced description becomes a function of the clustering algorithm. In the present work, we explore the dependence of the ML of the CG Helmholtz FE potential on the clustering algorithm. We consider coarse-graining based on partitional (k-means, leading to Voronoi diagram) and hierarchical agglomerative (bottom-up) clustering algorithms common in unsupervised ML and develop theory connecting the MSCG/FM learned CG Helmholtz potential and the clustering statistics. By combining the agglomerative clustering and the MSCG/FM learning in a recursive manner, we propose an efficient ML methodology to develop the fine-to-low resolution hierarchies of the CG models. The methodology does not suffer from degrading accuracy or increased computational cost to construct larger hierarchies and as such does not impose an upper size limitation of the CG particles resulting from the extended hierarchies. The utility of the methodology is demonstrated by obtaining the bottom-up agglomerative hierarchy for liquid nitromethane from all-atom molecular dynamics (MD) simulations. For agglomerative hierarchies, we prove the existence of renormalization group transformations that indicate self-similarity and allow for learning the low-resolution MSCG/FM potentials at low computational cost by rescaling and renormalizing the certain finer-resolution members of the hierarchy. The hierarchies of the CG models can be used to carry out simulations under constant-pressure conditions.

Machine learning assisted coarse-grained molecular dynamics modeling of meso-scale interfacial fluids

A hallmark of meso-scale interfacial fluids is the multi-faceted, scale-dependent interfacial energy, which often manifests different characteristics across the molecular and continuum scale. The multi-scale nature imposes a challenge to construct reliable coarse-grained (CG) models, where the CG potential function needs to faithfully encode the many-body interactions arising from the unresolved atomistic interactions and account for the heterogeneous density distributions across the interface. We construct the CG models of both single- and two-component polymeric fluid systems based on the recently developed deep coarse-grained potential [Zhang et al., J. Chem. Phys. 149, 034101 (2018)] scheme, where each polymer molecule is modeled as a CG particle. By only using the training samples of the instantaneous force under the thermal equilibrium state, the constructed CG models can accurately reproduce both the probability density function of the void formation in bulk and the spectrum of the capillary wave across the fluid interface. More importantly, the CG models accurately predict the volume-to-area scaling transition for the apolar solvation energy, illustrating the effectiveness to probe the meso-scale collective behaviors encoded with molecular-level fidelity.

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.

A hybrid coarse-grained model for structure, solvation and assembly of lipid-like peptides

Reconstituted photosynthetic proteins which are activated upon exposure to solar energy hold enormous potential for powering future solid state devices and solar cells. The functionality and integration of these proteins into such devices has been successfully enabled by lipid-like peptides. Yet, a fundamental understanding of the organization of these peptides with respect to the photosynthetic proteins and themselves remains unknown and is critical for guiding the design of such light-activated devices. This study investigates the relative organization of one such peptide sequence V₆K₂ (V: valine and K: lysine) within assemblies. Given the expansive spatiotemporal scales associated with this study, a hybrid coarse-grained (CG) model which captures the structure, conformation and aggregation of the peptide is adopted. The CG model uses a combination of iterative Boltzmann inversion and force matching to provide insight into the relative organization of V₆K₂ in assemblies. The CG model reproduces the structure of a V₆K₂ peptide sequence along with its all atom (AA) solvation structure. The relative organization of multiple peptides in an assembly, as captured by CG simulations, is in agreement with corresponding results from AA simulations. Also, a backmapping procedure reintroduces the AA details of the peptides within the aggregates captured by the CG model to demonstrate the relative organization of the peptides. Furthermore, a large number of peptides self-assemble into an elongated micelle in the CG simulation, which is consistent with experimental findings. The coarse-graining procedure is tested for transferability to longer peptide sequences, and hence can be extended to other amphiphilic peptide sequences.

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.

Multi-body effects in a coarse-grained protein force field

The use of coarse-grained (CG) models is a popular approach to study complex biomolecular systems. By reducing the number of degrees of freedom, a CG model can explore long time- and length-scales inaccessible to computational models at higher resolution. If a CG model is designed by formally integrating out some of the system's degrees of freedom, one expects multi-body interactions to emerge in the effective CG model's energy function. In practice, it has been shown that the inclusion of multi-body terms indeed improves the accuracy of a CG model. However, no general approach has been proposed to systematically construct a CG effective energy that includes arbitrary orders of multi-body terms. In this work, we propose a neural network based approach to address this point and construct a CG model as a multi-body expansion. By applying this approach to a small protein, we evaluate the relative importance of the different multi-body terms in the definition of an accurate model. We observe a slow convergence in the multi-body expansion, where up to five-body interactions are needed to reproduce the free energy of an atomistic model.

Conservative Potentials for a Lattice-Mapped Coarse-Grained Scheme

The conservative potential, arising from a coarse-grain (CG) mapping scheme for nonbonded atomistic particles, is studied. This is a bottom-up approach from first-principles that maps atomistic particles to fluid element-like subcells whose centers lie on a regular, cubic lattice. Unlike standard CG mapping schemes, the current one uses dynamic labeling which on-the-fly changes the CG labels of the particles. The subcells can also be different sizes and shapes, in principle. Equilibrium atomistic molecular dynamics trajectories for different Lennard-Jones fluids are calculated and converted to CG ones, from which CG probability distribution functions are calculated. Correlation studies show position and mass CG variables are uncoupled in a given subcell, as are different vector components of position. Furthermore, the strongest coupling occurs with neighboring cells in specific directions, and the resulting distribution is well described by a multivariate Gaussian. This implies the CG potential has a generalized quadratic form, whose derivative can be determined analytically. A microscopic rationalization is provided for the signs and relative magnitudes of different correlation coefficients, and in some cases, a connection is made with bulk properties of the fluid. We argue the generalized quadratic form should be robust to changes in the particulars of the CG scheme, as well as the nature of the atomistic intermolecular potential. Only a few potential parameters need to be calculated from the underlying atomistic system. This is significant because it indicates the transferability of this form to other, more complex systems. This transferability will be tested in future work, where mapping schemes with fuzzy boundaries will be considered.

Ultra-coarse-graining of homopolymers in inhomogeneous systems

We develop coarse-grained (CG) models for simulating homopolymers in inhomogeneous systems, focusing on polymer films and droplets. If the CG polymers interact solely through two-body potentials, then the films and droplets either dissolve or collapse into small aggregates, depending on whether the effective polymer-polymer interactions have been determined from reference simulations in the bulk or at infinite dilution. To address this shortcoming, we include higher order interactions either through an additional three-body potential or a local density-dependent potential (LDP). We parameterize the two- and three-body potentials via force matching, and the LDP through relative entropy minimization. While the CG models with three-body interactions fail at reproducing stable polymer films and droplets, CG simulations with an LDP are able to do so. Minor quantitative differences between the reference and the CG simulations, namely a slight broadening of interfaces accompanied by a smaller surface tension in the CG simulations, can be attributed to the deformation of polymers near the interfaces, which cannot be resolved in the CG representation, where the polymers are mapped to spherical beads.

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 new one-site coarse-grained model for water: Bottom-up many-body projected water (BUMPer). I. General theory and model

Water is undoubtedly one of the most important molecules for a variety of chemical and physical systems, and constructing precise yet effective coarse-grained (CG) water models has been a high priority for computer simulations. To recapitulate important local correlations in the CG water model, explicit higher-order interactions are often included. However, the advantages of coarse-graining may then be offset by the larger computational cost in the model parameterization and simulation execution. To leverage both the computational efficiency of the CG simulation and the inclusion of higher-order interactions, we propose a new statistical mechanical theory that effectively projects many-body interactions onto pairwise basis sets. The many-body projection theory presented in this work shares similar physics from liquid state theory, providing an efficient approach to account for higher-order interactions within the reduced model. We apply this theory to project the widely used Stillinger-Weber three-body interaction onto a pairwise (two-body) interaction for water. Based on the projected interaction with the correct long-range behavior, we denote the new CG water model as the Bottom-Up Many-Body Projected Water (BUMPer) model, where the resultant CG interaction corresponds to a prior model, the iteratively force-matched model. Unlike other pairwise CG models, BUMPer provides high-fidelity recapitulation of pair correlation functions and three-body distributions, as well as N-body correlation functions. BUMPer extensively improves upon the existing bottom-up CG water models by extending the accuracy and applicability of such models while maintaining a reduced computational cost.

Atomistic detailed free-energy landscape of intrinsically disordered protein studied by multi-scale divide-and-conquer molecular dynamics simulation

Calcineurin (CaN) is a eukaryotic serine/threonine protein phosphatase activated by both Ca²⁺ and calmodulin (CaM), including intrinsically disordered region (IDR). The region undergoes folding into an α-helix form in the presence Ca²⁺ -loaded CaM. To sample the ordered structure of the IDR by conventional all atom model (AAM) molecular dynamics (MD) simulation, the IDR and Ca²⁺ -loaded CaM must be simultaneously treated. However, it is time-consuming task because the coupled folding and binding should include repeated binding and dissociation. Then, in this study, we propose novel multi-scale divide-and-conquer MD (MSDC-MD), which combines AAM-MD and coarse-grained model MD (CGM-MD). To speed up the conformation sampling, MSDC-MD simulation first treats the IDR by CGM to sample conformations from wide conformation space; then, multiple AAM-MD in a limited area is initiated using the resultant CGM conformation, which is reconstructed by homology modeling method. To investigate performance, we sampled the ordered conformation of the IDR using MSDC-MD; the root-mean-square distance (RMSD) with respect to the experimental structure was 2.23 Å.

Development of Coarse-Grained Models for Poly(4-vinylphenol) and Poly(2-vinylpyridine): Polymer Chemistries with Hydrogen Bonding

In this paper, we identify the modifications needed in a recently developed generic coarse-grained (CG) model that captured directional interactions in polymers to specifically represent two exemplary hydrogen bonding polymer chemistries-poly(4-vinylphenol) and poly(2-vinylpyridine). We use atomistically observed monomer-level structures (e.g., bond, angle and torsion distribution) and chain structures (e.g., end-to-end distance distribution and persistence length) of poly(4-vinylphenol) and poly(2-vinylpyridine) in an explicitly represented good solvent (tetrahydrofuran) to identify the appropriate modifications in the generic CG model in implicit solvent. For both chemistries, the modified CG model is developed based on atomistic simulations of a single 24-mer chain. This modified CG model is then used to simulate longer (36-mer) and shorter (18-mer and 12-mer) chain lengths and compared against the corresponding atomistic simulation results. We find that with one to two simple modifications (e.g., incorporating intra-chain attraction, torsional constraint) to the generic CG model, we are able to reproduce atomistically observed bond, angle and torsion distributions, persistence length, and end-to-end distance distribution for chain lengths ranging from 12 to 36 monomers. We also show that this modified CG model, meant to reproduce atomistic structure, does not reproduce atomistically observed chain relaxation and hydrogen bond dynamics, as expected. Simulations with the modified CG model have significantly faster chain relaxation than atomistic simulations and slower decorrelation of formed hydrogen bonds than in atomistic simulations, with no apparent dependence on chain length.

Kernel-Based Machine Learning for Efficient Simulations of Molecular Liquids

Current machine learning (ML) models aimed at learning force fields are plagued by their high computational cost at every integration time step. We describe a number of practical and computationally efficient strategies to parametrize traditional force fields for molecular liquids from ML: the particle decomposition ansatz to two- and three-body force fields, the use of kernel-based ML models that incorporate physical symmetries, the incorporation of switching functions close to the cutoff, and the use of covariant meshing to boost the training set size. Results are presented for model molecular liquids: pairwise Lennard-Jones, three-body Stillinger-Weber, and bottom-up coarse-graining of water. Here, covariant meshing proves to be an efficient strategy to learn canonically averaged instantaneous forces. We show that molecular dynamics simulations with tabulated two- and three-body ML potentials are computationally efficient and recover two- and three-body distribution functions. Many-body representations, decomposition, and kernel regression schemes are all implemented in the open-source software package VOTCA.

Direct route to reproducing pair distribution functions with coarse-grained models via transformed atomistic cross correlations

Coarse-grained (CG) models are often parameterized to reproduce one-dimensional structural correlation functions of an atomically detailed model along the degrees of freedom governing each interaction potential. While cross correlations between these degrees of freedom inform the optimal set of interaction parameters, the correlations generated from the higher-resolution simulations are often too complex to act as an accurate proxy for the CG correlations. Instead, the most popular methods determine the interaction parameters iteratively while assuming that individual interactions are uncorrelated. While these iterative methods have been validated for a wide range of systems, they also have disadvantages when parameterizing models for multicomponent systems or when refining previously established models to better reproduce particular structural features. In this work, we propose two distinct approaches for the direct (i.e., noniterative) parameterization of a CG model by adjusting the high-resolution cross correlations of an atomistic model in order to more accurately reflect correlations that will be generated by the resulting CG model. The derived models more accurately describe the low-order structural features of the underlying AA model while necessarily generating inherently distinct cross correlations compared with the atomically detailed reference model. We demonstrate the proposed methods for a one-site-per-molecule representation of liquid water, where pairwise interactions are incapable of reproducing the true tetrahedral solvation structure. We then investigate the precise role that distinct cross-correlation features play in determining the correct pair correlation functions, evaluating the importance of the placement of correlation features as well as the balance between features appearing in different solvation shells.

Model reduction of rigid-body molecular dynamics via generalized multipole potentials

Motivated by the challenges of uncertainty quantification for coarse-grained (CG) molecular dynamics, we investigate the role of perturbation theory in model reduction of classical systems. In particular, we consider the task of coarse-graining rigid bodies in the context of generalized multipole potentials that have controllable levels of accuracy relative to their atomistic counterparts. We show how the multipole framework yields a hierarchy of models that systematically connects a CG "point molecule" approximation to the exact dynamics. We use these results to understand when and how the CG models fail to describe atomistic dynamics at the trajectory level and develop asymptotic error estimates for approximate molecular potential energies. Implications for other model-reduction strategies are also discussed. Key findings of this work are that (i) omitting rotational energy introduces significant error when coarse-graining and (ii) attention to symmetry can improve accuracy of "point-molecule" approximations. Analytical derivations and numerical results support these conclusions. Relevance to nonrigid bodies is also discussed.

Compatible observable decompositions for coarse-grained representations of real molecular systems

Coarse-grained (CG) observable expressions, such as pressure or potential energy, are generally different than their fine-grained (FG, e.g., atomistic) counterparts. Recently, we analyzed this so-called "representability problem" in Wagner et al. [J. Chem. Phys. 145, 044108 (2016)]. While the issue of representability was clearly and mathematically stated in that work, it was not made clear how to actually determine CG observable expressions from the underlying FG systems that can only be simulated numerically. In this work, we propose minimization targets for the CG observables of such systems. These CG observables are compatible with each other and with structural observables. Also, these CG observables are systematically improvable since they are variationally minimized. Our methods are local and data efficient because we decompose the observable contributions. Hence, our approaches are called the multiscale compatible observable decomposition (MS-CODE) and the relative entropy compatible observable decomposition (RE-CODE), which reflect two main approaches to the "bottom-up" coarse-graining of real FG systems. The parameterization of these CG observable expressions requires the introduction of new, symmetric basis sets and one-body terms. We apply MS-CODE and RE-CODE to 1-site and 2-site CG models of methanol for the case of pressure, as well as to 1-site methanol and acetonitrile models for potential energy.

Assessing the transferability of common top-down and bottom-up coarse-grained molecular models for molecular mixtures

Assessing the performance of top-down and bottom-up coarse-graining approaches.

Coarse-graining simulation approaches for polymer melts: the effect of potential range on computational efficiency

The integral equation coarse-graining (IECG) approach is a promising high-level coarse-graining (CG) method for polymer melts, with variable resolution from soft spheres to multi CG sites, which preserves the structural and thermodynamical consistencies with the related atomistic simulations. When compared to the atomistic description, the procedure of coarse-graining results in smoother free energy surfaces, longer-ranged potentials, a decrease in the number of interaction sites for a given polymer, and more. Because these changes have competing effects on the computational efficiency of the CG model, care needs to be taken when studying the effect of coarse-graining on the computational speed-up in CG molecular dynamics simulations. For instance, treatment of long-range CG interactions requires the selection of cutoff distances that include the attractive part of the effective CG potential and force. In particular, we show how the complex nature of the range and curvature of the effective CG potential, the selection of a suitable CG timestep, the choice of the cutoff distance, the molecular dynamics algorithms, and the smoothness of the CG free energy surface affect the efficiency of IECG simulations. By direct comparison with the atomistic simulations of relatively short chain polymer melts, we find that the overall computational efficiency is highest for the highest level of CG (soft spheres), with an overall improvement of the computational efficiency being about 106-108 for various CG levels/resolutions. Therefore, the IECG method can have important applications in molecular dynamics simulations of polymeric systems. Finally, making use of the standard spatial decomposition algorithm, the parallel scalability of the IECG simulations for various levels of CG is presented. Optimal parallel scaling is observed for a reasonably large number of processors. Although this study is performed using the IECG approach, its results on the relation between the level of CG and the computational efficiency are general and apply to any properly-constructed CG model.

Understanding three-body contributions to coarse-grained force fields

Coarse-graining is a systematic reduction of the number of degrees of freedom used to describe a system of interest. Coarse-graining can be thought of as a projection on the coarse-grained degrees of freedom and is therefore dependent on the number and type of basis functions used to represent the coarse-grained force field. We show that many-body extensions of the coarse-grained force field can result in substantial changes of the two-body interactions, making them much more attractive at short distances. This interplay can be alleviated by first parametrizing the two-body potential and then fitting the additional three-body contribution to the residual forces. The approach is illustrated on liquid water where three-body interactions are essential to reproduce the structural properties, and liquid methanol where two-body interactions are sufficient to reproduce the main structural features of the atomistic system. Furthermore, we demonstrate that the structural and thermodynamic accuracy of the coarse-grained models can be controlled by varying the magnitude of the three-body interactions. Our findings motivate basis set extensions which separate the many-body contributions of different order.

DeePCG: Constructing coarse-grained models via deep neural networks

We introduce a general framework for constructing coarse-grained potential models without ad hoc approximations such as limiting the potential to two- and/or three-body contributions. The scheme, called the Deep Coarse-Grained Potential (abbreviated DeePCG), exploits a carefully crafted neural network to construct a many-body coarse-grained potential. The network is trained with full atomistic data in a way that preserves the natural symmetries of the system. The resulting model is very accurate and can be used to sample the configurations of the coarse-grained variables in a much faster way than with the original atomistic model. As an application, we consider liquid water and use the oxygen coordinates as the coarse-grained variables, starting from a full atomistic simulation of this system at the ab initio molecular dynamics level. We find that the two-body, three-body, and higher-order oxygen correlation functions produced by the coarse-grained and full atomistic models agree very well with each other, illustrating the effectiveness of the DeePCG model on a rather challenging task.

Extending the range and physical accuracy of coarse-grained models: Order parameter dependent interactions

Order parameters (i.e., collective variables) are often used to describe the behavior of systems as they capture different features of the free energy surface. Yet, most coarse-grained (CG) models only employ two- or three-body non-bonded interactions between the CG particles. In situations where these interactions are insufficient for the CG model to reproduce the structural distributions of the underlying fine-grained (FG) model, additional interactions must be included. In this paper, we introduce an approach to expand the basis sets available in the multiscale coarse-graining (MS-CG) methodology by including order parameters. Then, we investigate the ability of an additive local order parameter (e.g., density) and an additive global order parameter (i.e., distance from a hard wall) to improve the description of CG models in interfacial systems. Specifically, we study methanol liquid-vapor coexistence, acetonitrile liquid-vapor coexistence, and acetonitrile liquid confined by hard-wall plates, all using single site CG models. We find that the use of order parameters dramatically improves the reproduction of structural properties of interfacial CG systems relative to the FG reference as compared with pairwise CG interactions alone.

Coarse-Grained Simulations of Aqueous Thermoresponsive Polyethers

Thermoresponsive polymers can change structure or solubility as a function of temperature. Block co-polymers of polyethers have a response that depends on polymer molecular weight and co-polymer composition. A coarse-grained model for aqueous polyethers is developed and applied to polyethylene oxide and polyethylene oxide-polypropylene oxide-polyethylene oxide triblock co-polymers. In this model, no interaction sites on hydrogen atoms are included, no Coulombic interactions are present, and all interactions are short-ranged, treated with a combination of two- and three-body terms. Our simulations find that The triblock co-polymers tend to associate at temperatures above 350 K. The aggregation is stabilized by contact between The hydrophobic methyl groups on The propylene oxide monomers and involves a large, favorable change in entropy.

Transferable Coarse-Grained Models of Liquid-Liquid Equilibrium Using Local Density Potentials Optimized with the Relative Entropy

Bottom-up coarse-grained (CG) models are now regularly pursued to enable large length and time scale molecular simulations of complex, often macromolecular systems. However, predicting fluid phase equilibria using such models remains fundamentally challenging. A major problem stems from the typically low transferability of CG models beyond the densities and/or compositions at which they are parametrized, which is necessary if they are to describe distinct structural and thermodynamic properties unique to each phase. CG model transferability is compounded by the representation of the inherently multibody coarse interactions using pair potentials that neglect higher order effects. Here, we propose to construct transferable single site CG models of liquid mixtures by supplementing traditional CG pair interactions with local density potentials, which constitute a computationally inexpensive mean-field approach to describe many-body effects, in that site energies are modulated by the local solution environment. To illustrate the approach, we use intra- and interspecies local density potentials to develop CG models of benzene-water solutions that show impressive transferability in structural metrics (pair correlation functions, density profiles) throughout composition space, in contrast to pair-only CG representations. While further refinement may be necessary to represent more complex thermodynamic properties, like the liquid-liquid interfacial tension, the generality and improvement offered by the local density approach are highly encouraging for enabling complex phase equilibrium modeling using CG models.

Theory of wavelet-based coarse-graining hierarchies for molecular dynamics

We present a multiresolution approach to compressing the degrees of freedom and potentials associated with molecular dynamics, such as the bond potentials. The approach suggests a systematic way to accelerate large-scale molecular simulations with more than two levels of coarse graining, particularly applications of polymeric materials. In particular, we derive explicit models for (arbitrarily large) linear (homo)polymers and iterative methods to compute large-scale wavelet decompositions from fragment solutions. This approach does not require explicit preparation of atomistic-to-coarse-grained mappings, but instead uses the theory of diffusion wavelets for graph Laplacians to develop system-specific mappings. Our methodology leads to a hierarchy of system-specific coarse-grained degrees of freedom that provides a conceptually clear and mathematically rigorous framework for modeling chemical systems at relevant model scales. The approach is capable of automatically generating as many coarse-grained model scales as necessary, that is, to go beyond the two scales in conventional coarse-grained strategies; furthermore, the wavelet-based coarse-grained models explicitly link time and length scales. Furthermore, a straightforward method for the reintroduction of omitted degrees of freedom is presented, which plays a major role in maintaining model fidelity in long-time simulations and in capturing emergent behaviors.

Many-Body Coarse-Grained Interactions Using Gaussian Approximation Potentials

We introduce a computational framework that is able to describe general many-body coarse-grained (CG) interactions of molecules and use it to model the free energy surface of molecular liquids as a cluster expansion in terms of monomer, dimer, and trimer terms. The contributions to the free energy due to these terms are inferred from all-atom molecular dynamics (MD) data using Gaussian Approximation Potentials, a type of machine-learning model that employs Gaussian process regression. The resulting CG model is much more accurate than those possible using pair potentials. Though slower than the latter, our model can still be faster than all-atom simulations for solvent-free CG models commonly used in biomolecular simulations.

Coarse-grained models using local-density potentials optimized with the relative entropy: Application to implicit solvation

Bottom-up multiscale techniques are frequently used to develop coarse-grained (CG) models for simulations at extended length and time scales but are often limited by a compromise between computational efficiency and accuracy. The conventional approach to CG nonbonded interactions uses pair potentials which, while computationally efficient, can neglect the inherently multibody contributions of the local environment of a site to its energy, due to degrees of freedom that were coarse-grained out. This effect often causes the CG potential to depend strongly on the overall system density, composition, or other properties, which limits its transferability to states other than the one at which it was parameterized. Here, we propose to incorporate multibody effects into CG potentials through additional nonbonded terms, beyond pair interactions, that depend in a mean-field manner on local densities of different atomic species. This approach is analogous to embedded atom and bond-order models that seek to capture multibody electronic effects in metallic systems. We show that the relative entropy coarse-graining framework offers a systematic route to parameterizing such local density potentials. We then characterize this approach in the development of implicit solvation strategies for interactions between model hydrophobes in an aqueous environment.

On the representability problem and the physical meaning of coarse-grained models

In coarse-grained (CG) models where certain fine-grained (FG, i.e., atomistic resolution) observables are not directly represented, one can nonetheless identify indirect the CG observables that capture the FG observable's dependence on CG coordinates. Often, in these cases it appears that a CG observable can be defined by analogy to an all-atom or FG observable, but the similarity is misleading and significantly undermines the interpretation of both bottom-up and top-down CG models. Such problems emerge especially clearly in the framework of the systematic bottom-up CG modeling, where a direct and transparent correspondence between FG and CG variables establishes precise conditions for consistency between CG observables and underlying FG models. Here we present and investigate these representability challenges and illustrate them via the bottom-up conceptual framework for several simple analytically tractable polymer models. The examples provide special focus on the observables of configurational internal energy, entropy, and pressure, which have been at the root of controversy in the CG literature, as well as discuss observables that would seem to be entirely missing in the CG representation but can nonetheless be correlated with CG behavior. Though we investigate these problems in the framework of systematic coarse-graining, the lessons apply to top-down CG modeling also, with crucial implications for simulation at constant pressure and surface tension and for the interpretations of structural and thermodynamic correlations for comparison to experiment.

Bayesian calibration of coarse-grained forces: Efficiently addressing transferability

Generating and calibrating forces that are transferable across a range of state-points remains a challenging task in coarse-grained (CG) molecular dynamics. In this work, we present a coarse-graining workflow, inspired by ideas from uncertainty quantification and numerical analysis, to address this problem. The key idea behind our approach is to introduce a Bayesian correction algorithm that uses functional derivatives of CG simulations to rapidly and inexpensively recalibrate initial estimates f0 of forces anchored by standard methods such as force-matching. Taking density-temperature relationships as a running example, we demonstrate that this algorithm, in concert with various interpolation schemes, can be used to efficiently compute physically reasonable force curves on a fine grid of state-points. Importantly, we show that our workflow is robust to several choices available to the modeler, including the interpolation schemes and tools used to construct f0. In a related vein, we also demonstrate that our approach can speed up coarse-graining by reducing the number of atomistic simulations needed as inputs to standard methods for generating CG forces.

Coarse-Graining of TIP4P/2005, TIP4P-Ew, SPC/E, and TIP3P to Monatomic Anisotropic Water Models Using Relative Entropy Minimization

Coarse-grained models are becoming a competitive alternative for modeling processes that occur over time and length scales beyond the reach of fully atomistic molecular simulations. Ideally, coarse-grained models should not only achieve high computational efficiency but also provide accurate predictions and fundamental insight into the role of molecular interactions, the characteristic behavior, and properties of the system they model. In this work we derive a series of monatomic coarse-grained water models mX(REM) from the most popular atomistic water models X = TIP3P, SPC/E, TIP4P-Ew, and TIP4P/2005, using the relative entropy minimization (REM) method. Each coarse-grained water molecule is represented by a single particle that interacts through short-ranged anisotropic interactions that encourage the formation of "hydrogen-bonded" structures. We systematically investigate the features of the coarse-grained models in reproducing over 20 structural, dynamic, and thermodynamic properties of the reference atomistic water models-including the existence and locus of the characteristic density anomaly. The mX(REM) coarse-grained models reproduce quite faithfully the radial and angular distribution function of water, produce a temperature of maximum density (TMD), and stabilize the ice I crystal. Moreover, the ratio between the TMD and the melting temperature of the crystal in the mX(REM) models and liquid-ice equilibrium properties show reasonable agreement with the results of the corresponding atomistic models. The mX(REM) models, however, severely underestimate the cohesive energy of the condensed water phases. We investigate which specific limitations of the coarse-grained models arise from the REM methodology, from the monatomic nature of the models, and from the Stillinger-Weber interaction potential form. Our analysis indicates that a small compromise in the accuracy of structural properties can result in a significant increase of the overall accuracy and representability of the coarse-grained water models. We evaluate the accuracy of the atomistic and the monatomic anisotropic coarse-grained water models, including the mW water model, in reproducing experimental water properties. We find that mW and mTIP4P/2005(REM) score closer to experiment than widely used atomistic water models. We conclude that monatomic models of water with short-range, anisotropic "hydrogen-bonding" three-body interactions can be competitive in accuracy with fully atomistic models for the study of a wide range of properties and phenomena at less than 1/100th of the computational cost.

Molecular basis for the dissociation dynamics of protein A-immunoglobulin G1 complex

Staphylococcus aureus protein A (SpA) is the most popular affinity ligand for immunoglobulin G1 (IgG1). However, the molecular basis for the dissociation dynamics of SpA-IgG1 complex is unclear. Herein, coarse-grained (CG) molecular dynamics (MD) simulations with the Martini force field were used to study the dissociation dynamics of the complex. The CG-MD simulations were first verified by the agreement in the structural and interactional properties of SpA and human IgG1 (hIgG1) in the association process between the CG-MD and all-atom MD at different NaCl concentrations. Then, the CG-MD simulation studies focused on the molecular insight into the dissociation dynamics of SpA-hIgG1 complex at pH 3.0. It is found that there are four steps in the dissociation process of the complex. First, there is a slight conformational adjustment of helix II in SpA. This is followed by the phenomena that the electrostatic interactions provided by the three hot spots (Glu143, Arg146 and Lys154) of helix II of SpA break up, leading to the dissociation of helix II from the binding site of hIgG1. Subsequently, breakup of the hydrophobic interactions between helix I (Phe132, Tyr133 and His137) in SpA and hIgG1 occurs, resulting in the disengagement of helix I from its binding site of hIgG1. Finally, the non-specific interactions between SpA and hIgG1 decrease slowly till disappearance, leading to the complete dissociation of the SpA-hIgG1 complex. This work has revealed that CG-MD coupled with the Martini force field is an effective method for studying the dissociation dynamics of protein-protein complex.