The Theory of Ultra-Coarse-Graining. 3. Coarse-Grained Sites with Rapid Local Equilibrium of Internal States

Evaluation of Essential Dynamics and Fixed-Length Coarse Graining for Multidomain Proteins

For multiscale modeling of biomolecules, reliable coarse-grained (CG) models can offer great potential to simulate larger temporal and spatial scales than traditional all-atom (AA) models. In this study, we explore the essential dynamics coarse graining (EDCG) and fixed-length coarse graining (FLCG) approaches for constructing highly coarse-grained models for multidomain proteins (MDPs), with 1 to 10 amino acid residues per CG site. In the studies of 13 MDPs, our data indicate that both EDCG and FLCG can preserve the protein dynamics of MDPs. FLCG, which restricts an equal number of residues in each CG site, represents an excellent approximation to EDCG and a straightforward approach for coarse-graining MDPs. Furthermore, FLCG is tested with a class B G-protein-coupled receptor protein, and the agreement with prior experiments suggests its general application to various MDPs in different environments or conditions. Finally, we demonstrate another application of FLCG through progressive backmapping, showcasing the ability to recover from lower-resolution CG models (6 residues/CG site) to higher-resolution ones (1 residue/CG site). These promising outcomes underscore the broad applicability of FLCG to construct highly or ultra-coarse-grained models of complex biomolecules for multiscale simulations.

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.

Utilizing Machine Learning to Greatly Expand the Range and Accuracy of Bottom-Up Coarse-Grained Models through Virtual Particles

Coarse-grained (CG) models parametrized using atomistic reference data, i.e., "bottom up" CG models, have proven useful in the study of biomolecules and other soft matter. However, the construction of highly accurate, low resolution CG models of biomolecules remains challenging. We demonstrate in this work how virtual particles, CG sites with no atomistic correspondence, can be incorporated into CG models within the context of relative entropy minimization (REM) as latent variables. The methodology presented, variational derivative relative entropy minimization (VD-REM), enables optimization of virtual particle interactions through a gradient descent algorithm aided by machine learning. We apply this methodology to the challenging case of a solvent-free CG model of a 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) lipid bilayer and demonstrate that introduction of virtual particles captures solvent-mediated behavior and higher-order correlations which REM alone cannot capture in a more standard CG model based only on the mapping of collections of atoms to the CG sites.

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.

Enhancing pressure consistency and transferability of structure-based coarse-graining

Coarse-graining, which models molecules with coarse-grained (CG) beads, allows molecular dynamics simulations to be applied to systems with large length and time scales while preserving the essential molecular structure. However, CG models generally have insufficient representability and transferability. A commonly used method to resolve this problem is multi-state iterative Boltzmann inversion (MS-IBI) with pressure correction, which matches both the structural properties and pressures at different thermodynamic states between CG and all-atom (AA) simulations. Nevertheless, this method is usually effective only in a narrow pressure range. In this paper, we propose a modified CG scheme to overcome this limitation. We find that the fundamental reason for this limitation is that CG beads at close distances are ellipsoids rather than isotropically compressed spheres, as described in conventional CG models. Hence, we propose a method to compensate for such differences by slightly modifying the radial distribution functions (RDFs) derived from AA simulations and using the modified RDFs as references for pressure-corrected MS-IBI. We also propose a method to determine the initial non-bonded potential using both the target RDF and pressure. Using n-dodecane as a case study, we demonstrate that the CG model developed using our scheme reproduces the RDFs and pressures over a wide range of pressure states, including three reference low-pressure states and two test high-pressure states. The proposed scheme allows for accurate CG simulations of systems in which pressure or density varies with time and/or position.

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 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.

Physical Characterization of Triolein and Implications for Its Role in Lipid Droplet Biogenesis

Lipid droplets (LDs) are neutral lipid-storing organelles surrounded by a phospholipid (PL) monolayer. At present, how LDs are formed in the endoplasmic reticulum (ER) bilayer is poorly understood. In this study, we present a revised all-atom (AA) triolein (TG) model, the main constituent of the LD core, and characterize its properties in a bilayer membrane to demonstrate the implications of its behavior in LD biogenesis. In bilayer simulations, TG resides at the surface, adopting PL-like conformations (denoted in this work as SURF-TG). Free energy sampling simulation results estimate the barrier for TG relocating from the bilayer surface to the bilayer center to be ∼2 kcal/mol in the absence of an oil lens. SURF-TG is able to modulate membrane properties by increasing PL ordering, decreasing bending modulus, and creating local negative curvature. The other neutral lipid, dioleoyl-glycerol (DAG), also reduces the membrane bending modulus and populates negative curvature regions. A phenomenological coarse-grained (CG) model is also developed to observe larger-scale SURF-TG-mediated membrane deformation. CG simulations confirm that TG nucleates between the bilayer leaflets at a critical concentration when SURF-TG is evenly distributed. However, when one monolayer contains more SURF-TG, the membrane bends toward the other leaflet, followed by TG nucleation if a concentration is higher than the critical threshold. The central conclusion of this study is that SURF-TG is a negative curvature inducer, as well as a membrane modulator. To this end, a model is proposed in which the accumulation of SURF-TG in the luminal leaflet bends the ER bilayer toward the cytosolic side, followed by TG nucleation.

From System Modeling to System Analysis: The Impact of Resolution Level and Resolution Distribution in the Computer-Aided Investigation of Biomolecules

The ever increasing computer power, together with the improved accuracy of atomistic force fields, enables researchers to investigate biological systems at the molecular level with remarkable detail. However, the relevant length and time scales of many processes of interest are still hardly within reach even for state-of-the-art hardware, thus leaving important questions often unanswered. The computer-aided investigation of many biological physics problems thus largely benefits from the usage of coarse-grained models, that is, simplified representations of a molecule at a level of resolution that is lower than atomistic. A plethora of coarse-grained models have been developed, which differ most notably in their granularity; this latter aspect determines one of the crucial open issues in the field, i.e. the identification of an optimal degree of coarsening, which enables the greatest simplification at the expenses of the smallest information loss. In this review, we present the problem of coarse-grained modeling in biophysics from the viewpoint of system representation and information content. In particular, we discuss two distinct yet complementary aspects of protein modeling: on the one hand, the relationship between the resolution of a model and its capacity of accurately reproducing the properties of interest; on the other hand, the possibility of employing a lower resolution description of a detailed model to extract simple, useful, and intelligible information from the latter.

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.

Structure and dynamics of responsive colloids with dynamical polydispersity

Dynamical polydispersity in single-particle properties, for example a fluctuating particle size, shape, charge density, etc, is intrinsic to responsive colloids (RCs), such as biomacromolecules or microgels, but is typically not resolved in coarse-grained mesoscale simulations. Here, we present Brownian dynamics simulations of suspensions of RCs modeling soft hydrogel colloids, for which the size of the individual particles is an explicitly resolved (Gaussian) degree of freedom and dynamically responds to the local interacting environment. We calculate the liquid structure, emergent size distributions, long-time diffusion, and property (size) relaxation kinetics for a wide range of densities and intrinsic property relaxation times in the canonical ensemble. Comparison to interesting reference cases, such as conventional polydisperse suspensions with a frozen parent distribution, or conventional monodisperse systems interacting with an effective pair potential for one fixed size, shows a significant spread in the structure and dynamics. The differences, most apparent in the high density regimes, are due to many-body correlations and the dynamical coupling between property and translation in RC systems, not explicitly accounted for in the conventional treatments. In particular, the translational diffusion in the RC systems is surprisingly close to the free (single RC) diffusion, mainly due to a cancellation of crowding and size compression effects. We show that an effective monodisperse pair potential can be constructed that describes the many-body correlations reasonably well by convoluting the RC pair potential with the density-dependent emergent size distributions and using a mean effective diffusion constant.

Compressive and Tensile Deformations Alter ATP Hydrolysis and Phosphate Release Rates in Actin Filaments

Actin filament networks in eukaryotic cells are constantly remodeled through nucleotide state controlled interactions with actin binding proteins, leading to macroscopic structures such as bundled filaments, branched filaments, and so on. The nucleotide (ATP) hydrolysis, phosphate release, and polymerization/depolymerization reactions that lead to the formation of these structures are correlated with the conformational fluctuations of the actin subunits at the molecular scale. The resulting structures generate and experience varying levels of force and impart cells with several functionalities such as their ability to move, divide, transport cargo, etc. Models that explicitly connect the structure to reactions are essential to elucidate a fundamental level of understanding of these processes. In this regard, a bottom-up Ultra-Coarse-Grained (UCG) model of actin filaments that can simulate ATP hydrolysis, inorganic phosphate release (Pi), and depolymerization reactions is presented in this work. In this model, actin subunits are represented using coarse-grained particles that evolve in time and undergo reactions depending on the conformations sampled. The reactions are represented through state transitions, with each state represented by a unique effective coarse-grained potential. Effects of compressive and tensile strains on the rates of reactions are then analyzed. Compressive strains tend to unflatten the actin subunits, reduce the rate of ATP hydrolysis, and increase the Pi release rate. On the other hand, tensile strain flattens subunits, increases the rate of ATP hydrolysis, and decrease the Pi release rate. Incorporating these predictions into a Markov State Model highlighted that strains alter the steady-state distribution of subunits with ADPPi and ADP nucleotide, thus identifying possible additional factors underlying the cooperative binding of regulatory proteins to actin filaments.

Coarse-Grained Simulation of Mechanical Properties of Single Microtubules With Micrometer Length

Microtubules are one of the most important components in the cytoskeleton and play a vital role in maintaining the shape and function of cells. Because single microtubules are some micrometers long, it is difficult to simulate such a large system using an all-atom model. In this work, we use the newly developed convolutional and K-means coarse-graining (CK-CG) method to establish an ultra-coarse-grained (UCG) model of a single microtubule, on the basis of the low electron microscopy density data of microtubules. We discuss the rationale of the micro-coarse-grained microtubule models of different resolutions and explore microtubule models up to 12-micron length. We use the devised microtubule model to quantify mechanical properties of microtubules of different lengths. Our model allows mesoscopic simulations of micrometer-level biomaterials and can be further used to study important biological processes related to microtubule function.

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.

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.

Double-Well Ultra-Coarse-Grained Model to Describe Protein Conformational Transitions

The double-well model is usually used to describe the conformational transition between two states of a protein. Since conformational changes usually occur within a relatively large time scale, coarse-grained models are often used to accelerate the dynamic process due to their inexpensive computational cost. In this work, we develop a double-well ultra-coarse-grained (DW-UCG) model to describe the conformational transitions of the adenylate kinase, glutamine-binding protein, and lactoferrin. The coarse-grained simulation results show that the DW-UCG model of adenylate kinase captures the crucial intermediate states in the LID-closing and NMP-closing pathways, reflecting the key secondary structural changes in the conformational transition. A comparison of the different DW-UCG models of adenylate kinase indicates that an appropriate choice of bead resolution could generate the free energy landscape that is comparable to that from the residue-based model. The coarse-grained simulations for the glutamine-binding protein and lactoferrin also demonstrate that the DW-UCG model is valid in reproducing the correct two-state behavior for their functional study, which indicates the potential application of the DW-UCG model in investigating the mechanism of conformational changes of large proteins.

Protein-ligand binding with the coarse-grained Martini model

The detailed understanding of the binding of small molecules to proteins is the key for the development of novel drugs or to increase the acceptance of substrates by enzymes. Nowadays, computer-aided design of protein–ligand binding is an important tool to accomplish this task. Current approaches typically rely on high-throughput docking essays or computationally expensive atomistic molecular dynamics simulations. Here, we present an approach to use the recently re-parametrized coarse-grained Martini model to perform unbiased millisecond sampling of protein–ligand interactions of small drug-like molecules. Remarkably, we achieve high accuracy without the need of any a priori knowledge of binding pockets or pathways. Our approach is applied to a range of systems from the well-characterized T4 lysozyme over members of the GPCR family and nuclear receptors to a variety of enzymes. The presented results open the way to high-throughput screening of ligand libraries or protein mutations using the coarse-grained Martini model.

Computer-aided design of protein-ligand binding is important for the development of novel drugs. Here authors present an approach to use the recently re-parametrized coarse-grained Martini model to perform unbiased millisecond sampling of protein-ligand binding interactions of small drug-like molecules.

Interfacial solvation and slow transport of hydrated excess protons in non-ionic reverse micelles

This work employs molecular dynamics simulations to investigate the solvation and transport properties of hydrated excess protons (with a hydronium-like core structure) in non-ionic Igepal CO-520 reverse micelles of various sizes in a non-polar solvent. Multiscale Reactive Molecular Dynamics (MS-RMD) simulations were used to describe vehicular and hopping diffusion during the proton transport process. As detailed herein, an excess proton shows a marked tendency to localize in the interfacial region of micellar water pools. Slow proton transport was observed which becomes faster with increasing micellar size. Further analysis reveals that the slow diffusion of an excess proton is a combined result of slow water diffusion and the low proton hopping rate. This study also confirms that a low proton hopping rate in reverse micelles stems from the interfacial solvation of hydrated excess protons and the immobilization of interfacial water. The low water density in the interfacial region makes it difficult to form a complete hydrogen bond network near the hydrated excess proton, and therefore locks in the orientation of hydrated proton cations. The immobilization of the interfacial water also slows the relaxation of the overall hydrogen bond network.

Determining Optimal Coarse-Grained Representation for Biomolecules Using Internal Cluster Validation Indexes

The development of ultracoarse-grained models for large biomolecules needs to derive the optimal number of coarse-grained (CG) sites to represent the targets. In this work, we propose to use the statistical internal cluster validation indexes to determine the optimal number of CG sites that are optimized based on the essential dynamics coarse-graining method. The calculated curves of Calinski-Harabasz and Silhouette Coefficient indexes exhibit the extrema corresponding to the similar CG numbers. The calculated ratios of the optimal CG numbers to the residue numbers of fine-grained models are in the range from 4 to 2. The comparison of the stability of index results indicates that Calinski-Harabasz index is the better choice to determine the optimal CG representation in coarse-graining. © 2019 Wiley Periodicals, Inc.

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

Coarse-grained (CG) models facilitate efficient simulation of complex systems by integrating out the atomic, or fine-grained (FG), degrees of freedom. Systematically derived CG models from FG simulations often attempt to approximate the CG potential of mean force (PMF), an inherently multidimensional and many-body quantity, using additive pairwise contributions. However, they currently lack fundamental principles that enable their extensible use across different thermodynamic state points, i.e., transferability. In this work, we investigate the explicit energy-entropy decomposition of the CG PMF as a means to construct transferable CG models. In particular, despite its high-dimensional nature, we find for liquid systems that the entropic component to the CG PMF can similarly be represented using additive pairwise contributions, which we show is highly coupled to the CG configurational entropy. This approach formally connects the missing entropy that is lost due to the CG representation, i.e., translational, rotational, and vibrational modes associated with the missing degrees of freedom, to the CG entropy. By design, the present framework imparts transferable CG interactions across different temperatures due to the explicit definition of an additive entropic contribution. Furthermore, we demonstrate that transferability across composition state points, such as between bulk liquids and their mixtures, is also achieved by designing combining rules to approximate cross-interactions from bulk CG PMFs. Using the predicted CG model for liquid mixtures, structural correlations of the fitted CG model were found to corroborate a high-fidelity combining rule. Our findings elucidate the physical nature and compact representation of CG entropy and suggest a new approach for overcoming the transferability problem. We expect that this approach will further extend the current view of CG modeling into predictive multiscale modeling.

A new method for the construction of coarse-grained models of large biomolecules from low-resolution cryo-electron microscopy data

The rapid development of cryo-electron microscopy (cryo-EM) has led to the generation of significant low-resolution electron density data of biomolecules. However, the atomistic details of huge biomolecules usually cannot be obtained because it is very difficult to construct all-atom models for MD simulations. Thus, it is still a challenge to make use of the rich low-resolution cryo-EM data for computer simulation and functional study. In this study, we proposed a new method called Convolutional and K-means Coarse-Graining (CK-CG) for the efficient coarse-graining of large biological systems. Using the CK-CG method, we could directly map the cryo-EM data into coarse-grained (CG) beads. Furthermore, the CG beads were parameterized with an empirical harmonic potential to construct a new CG model. We subjected the CK-CG models of the fibrillar protein assemblies F-actin and collagen to external forces in pulling dynamic simulations to assess their mechanical response. The agreement between the estimated tensile stiffness between CG models and experiments demonstrates the validity of the CK-CG method. Thus, our method provides a practical strategy for the direct construction of a structural model from low-resolution data for biological function studies.

Coarse-Grained Simulation of Full-Length Integrin Activation

Integrin conformational dynamics are critical to their receptor and signaling functions in many cellular processes, including spreading, adhesion, and migration. However, assessing integrin conformations is both experimentally and computationally challenging because of limitations in resolution and dynamic sampling. Thus, structural changes that underlie transitions between conformations are largely unknown. Here, focusing on integrin αvβ₃, we developed a modified form of the coarse-grained heterogeneous elastic network model (hENM), which allows sampling conformations at the onset of activation by formally separating local fluctuations from global motions. Both local fluctuations and global motions are extracted from all-atom molecular dynamics simulations of the full-length αvβ₃ bent integrin conformer, but whereas the former are incorporated in the hENM as effective harmonic interactions between groups of residues, the latter emerge by systematically identifying and treating weak interactions between long-distance domains with flexible and anharmonic connections. The new hENM model allows integrins and single-point mutant integrins to explore various conformational states, including the initiation of separation between α- and β-subunit cytoplasmic regions, headpiece extension, and legs opening.

Systematic Coarse-Grained Lipid Force Fields with Semiexplicit Solvation via Virtual Sites

Despite the central role of lipids in many biophysical functions, the molecular mechanisms that dictate macroscopic lipid behavior remain elusive to both experimental and computational approaches. As such, there has been much interest in the development of low-resolution, implicit-solvent coarse-grained (CG) models to dynamically simulate biologically relevant spatiotemporal scales with molecular fidelity. However, in the absence of solvent, a key challenge for CG models is to faithfully emulate solvent-mediated forces, which include both hydrophilic and hydrophobic interactions that drive lipid aggregation and self-assembly. In this work, we provide a new methodological framework to incorporate semiexplicit solvent effects through the use of virtual CG particles, which represent structural features of the solvent-lipid interface. To do so, we leverage two systematic coarse-graining approaches, multiscale coarse-graining (MS-CG) and relative entropy minimization (REM), in a hybrid fashion to construct our virtual-site CG (VCG) models. As a proof-of-concept, we focus our efforts on two lipid species, 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), which adopt a liquid-disordered and gel phase, respectively, at room temperature. Through our analysis, we also present, to our knowledge, the first direct comparison between the MS-CG and REM methods for a complex biomolecule and highlight each of their strengths and weaknesses. We further demonstrate that VCG models recapitulate the rich biophysics of lipids, which enable self-assembly, morphological diversity, and multiple phases. Our findings suggest that the VCG framework is a powerful approach for investigation into macromolecular biophysics.

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.

Multiscale simulation of actin filaments and actin-associated proteins

Actin is an important cytoskeletal protein that serves as a building block to form filament networks that span across the cell. These networks are orchestrated by a myriad of other cytoskeletal entities including the unbranched filament-forming protein formin and branched network-forming protein complex Arp2/3. Computational models have been able to provide insights into many important structural transitions that are involved in forming these networks, and into the nature of interactions essential for actin filament formation and for regulating the behavior of actin-associated proteins. In this review, we summarize a subset of such models that focus on the atomistic features and those that can integrate atomistic features into a larger picture in a multiscale fashion.

Advances in coarse-grained modeling of macromolecular complexes

Recent progress in coarse-grained (CG) molecular modeling and simulation has facilitated an influx of computational studies on biological macromolecules and their complexes. Given the large separation of length-scales and time-scales that dictate macromolecular biophysics, CG modeling and simulation are well-suited to bridge the microscopic and mesoscopic or macroscopic details observed from all-atom molecular simulations and experiments, respectively. In this review, we first summarize recent innovations in the development of CG models, which broadly include structure-based, knowledge-based, and dynamics-based approaches. We then discuss recent applications of different classes of CG models to explore various macromolecular complexes. Finally, we conclude with an outlook for the future in this ever-growing field of biomolecular modeling.

Ultra-Coarse-Grained Liquid State Models with Implicit Hydrogen Bonding

Coarse-graining (CG) methodologies have been widely used to extend the time and length scales of computer simulations by averaging over the atomistic details beneath the resolution of the CG models. Despite the efficiency of CG models, important configurational information during a given process may be lost at the CG resolution. One example of this is the topology of the hydrogen bonding network in the liquid state. When the functional group that participates in hydrogen bonding (e.g., -OH in methanol) is coarse-grained into one CG site, the effective interactions of the resultant CG model are usually derived from an averaged overall trajectory and, thus, do not take into account the hydrogen bonding interactions and topologies that are present at the all-atom resolution. In order to overcome this challenge, the present study develops new ultra-coarse-grained (UCG) models that include internal states within the CG sites that participate in hydrogen bonding, where each state represents a specific configuration such as the hydrogen bonding donor or acceptor. Internal states of the UCG beads are modeled to remain in quasi-equilibrium, and the degree of mixing is controlled by utilizing the effective local density of the UCG sites. In particular, we consider two groups of UCG models with different types of hydrogen bonding motifs: chain-like and ring-like. Using five different liquid systems that contain the same fundamental functional groups as biomolecules, we demonstrate the ability of the UCG models to reproduce the structural properties that originate from the configurations beneath the resolution of the UCG model. This proposed approach can also be applied to other liquids with such specific and directional interactions, or even to complex biomolecular systems in which hydrogen bonding is critical.

Insights into the Cooperative Nature of ATP Hydrolysis in Actin Filaments

Actin filaments continually assemble and disassemble within a cell. Assembled filaments "age" as a bound nucleotide ATP within each actin subunit quickly hydrolyzes followed by a slower release of the phosphate Pi, leaving behind a bound ADP. This subtle change in nucleotide state of actin subunits affects filament rigidity as well as its interactions with binding partners. We present here a systematic multiscale ultra-coarse-graining approach that provides a computationally efficient way to simulate a long actin filament undergoing ATP hydrolysis and phosphate-release reactions while systematically taking into account available atomistic details. The slower conformational changes and their dependence on the chemical reactions are simulated with the ultra-coarse-graining model by assigning internal states to the coarse-grained sites. Each state is represented by a unique potential surface of a local heterogeneous elastic network. Internal states undergo stochastic transitions that are coupled to conformations of the underlying molecular system. The model reproduces mechanical properties of the filament and allows us to study whether conformational fluctuations in actin subunits produce cooperative filament aging. We find that the nucleotide states of neighboring subunits modulate the reaction kinetics, implying cooperativity in ATP hydrolysis and Pi release. We further systematically coarse grain the system into a Markov state model that incorporates assembly and disassembly, facilitating a direct comparison with previously published models. We find that cooperativity in ATP hydrolysis and Pi release significantly affects the filament growth dynamics only near the critical G-actin concentration, whereas far from it, both cooperative and random mechanisms show similar growth dynamics. In contrast, filament composition in terms of the bound nucleotide distribution varies significantly at all monomer concentrations studied. These results provide new insights, to our knowledge, into the cooperative nature of ATP hydrolysis and Pi release and the implications it has for actin filament properties, providing novel predictions for future experimental studies.

Adaptive resolution simulations coupling atomistic water to dissipative particle dynamics

Multiscale methods are the most efficient way to address the interlinked spatiotemporal scales encountered in soft matter and molecular liquids. In the literature reported hybrid approaches span from quantum to atomistic, coarse-grained, and continuum length scales. In this article, we present the hybrid coupling of the molecular dynamics (MD) and dissipative particle dynamics (DPD) methods, bridging the micro- and mesoscopic descriptions. The interfacing is performed within the adaptive resolution scheme (AdResS), which is a linear momentum conserving coupling technique. Our methodology is hence suitable to simulate fluids on the micro/mesoscopic scale, where hydrodynamics plays an important role. The presented approach is showcased for water at ambient conditions. The supramolecular coupling is enabled by a recently developed clustering algorithm SWINGER that assembles, disassembles, and reassembles clusters as needed during the course of the simulation. This allows for a seamless coupling between standard atomistic MD and DPD models. The developed framework can be readily applied to various applications in the fields of materials and life sciences, e.g., simulations of phospholipids and polymer melts, or to study the red blood cells behavior in normal and disease states.

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.

Ultra-Coarse-Grained Models Allow for an Accurate and Transferable Treatment of Interfacial Systems

Interfacial systems are fundamentally important in many processes. However, constructing coarse-grained (CG) models for such systems is a significant challenge due to their inhomogeneous nature. This problem is made worse due to the generally nontransferable nature of the interactions in CG models across different phases. In this paper, we address these challenges by systematically constructing ultra-coarse-grained (UCG) models for interfaces, in which the CG sites are allowed to have internal states. We find that a multiscale coarse-grained (MS-CG) representation of a single CG site model fails to identify the directionality of a molecule and is unable to reproduce the correct phase coexistence for aspherical molecules. In contrast with conventional MS-CG models, the UCG methodology allows chemical and environmental changes to be captured by modulating the interactions between internal states. In this work, we design the internal states to depend on local particle density to distinguish different phases in liquid/vapor or liquid/liquid interfaces. These UCG models are able to capture phase coexistence and recapitulate structures, notably at state points in which the MS-CG method yields poor results. Interestingly, effective pairwise forces and potentials from the UCG models are almost identical to those of the bulk liquids that correspond to each phase, indicating that the UCG approach can provide transferable interactions. This approach is expected to be applicable to other systems that exhibit phase coexistence and also to complex macromolecular systems by modulating interactions based on local density or other order parameters to unravel the complex nature underlying heterogeneous system boundaries.

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.

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.