Long time atomistic polymer trajectories from coarse grained simulations: bisphenol-A polycarbonate

Reverse Mapping of Coarse Grained Polyglutamine Conformations from PRIME20 Sampling

An inverse coarse-graining protocol is presented for generating and validating atomistic structures of large (bio-) molecules from conformations obtained via a coarse-grained sampling method. Specifically, the protocol is implemented and tested based on the (coarse-grained) PRIME20 protein model (P20/SAMC), and the resulting all-atom conformations are simulated using conventional biomolecular force fields. The phase space sampling at the coarse-grained level is performed with a stochastical approximation Monte Carlo approach. The method is applied to a series of polypeptides, specifically dimers of polyglutamine with varying chain length in aqueous solution. The majority (>70 %) of the conformations obtained from the coarse-grained peptide model can successfully be mapped back to atomistic structures that remain conformationally stable during 10 ns of molecular dynamics simulations. This work can be seen as the first step towards the overarching goal of improving our understanding of protein aggregation phenomena through simulation methods.

Construction of full-atomistic polymer amorphous structures using reverse-mapping from Kremer-Grest models

We propose a method to build full-atomistic (FA) amorphous polymer structures using reverse-mapping from coarse-grained (CG) models. In this method, three models with different resolutions are utilized, namely the CG1, CG2, and FA models. It is assumed that the CG1 model is more abstract than the CG2 model. The CG1 is utilized to equilibrate the system, and then sequential reverse-mapping procedures from the CG1 to the CG2 models and from the CG2 to the FA models are conducted. A mapping relation between the CG1 and the FA models is necessary to generate a polymer structure with a given density and radius of chains. Actually, we have used the Kremer-Grest (KG) model as the CG1 and the monomer-level CG model as the CG2 model. Utilizing the mapping relation, we have developed a scheme that constructs an FA polymer model from the KG model. In the scheme, the KG model, the monomer level CG model, and the FA model are successively constructed. The scheme is applied to polyethylene (PE), cis 1,4-polybutadiene (PB), and poly(methyl methacrylate) (PMMA). As a validation, the structures of PE and PB constructed by the scheme were carefully checked through comparison with those obtained using long-time FA molecular dynamics (MD) simulations. We found that both short- and long-range chain structures constructed by the scheme reproduced those obtained by the FA MD simulations. Then, as an interesting application, the scheme is applied to generate an entangled PMMA structure. The results showed that the scheme provides an efficient and easy way to construct amorphous structures of FA polymers.

DiAMoNDBack: Diffusion-Denoising Autoregressive Model for Non-Deterministic Backmapping of Cα Protein Traces

Coarse-grained molecular models of proteins permit access to length and time scales unattainable by all-atom models and the simulation of processes that occur on long time scales, such as aggregation and folding. The reduced resolution realizes computational accelerations, but an atomistic representation can be vital for a complete understanding of mechanistic details. Backmapping is the process of restoring all-atom resolution to coarse-grained molecular models. In this work, we report DiAMoNDBack (Diffusion-denoising Autoregressive Model for Non-Deterministic Backmapping) as an autoregressive denoising diffusion probability model to restore all-atom details to coarse-grained protein representations retaining only Cα coordinates. The autoregressive generation process proceeds from the protein N-terminus to C-terminus in a residue-by-residue fashion conditioned on the Cα trace and previously backmapped backbone and side-chain atoms within the local neighborhood. The local and autoregressive nature of our model makes it transferable between proteins. The stochastic nature of the denoising diffusion process means that the model generates a realistic ensemble of backbone and side-chain all-atom configurations consistent with the coarse-grained Cα trace. We train DiAMoNDBack over 65k+ structures from the Protein Data Bank (PDB) and validate it in applications to a hold-out PDB test set, intrinsically disordered protein structures from the Protein Ensemble Database (PED), molecular dynamics simulations of fast-folding mini-proteins from DE Shaw Research, and coarse-grained simulation data. We achieve state-of-the-art reconstruction performance in terms of correct bond formation, avoidance of side-chain clashes, and the diversity of the generated side-chain configurational states. We make the DiAMoNDBack model publicly available as a free and open-source Python package.

Multiscale equilibration of highly entangled isotropic model polymer melts

We present a computationally efficient multiscale method for preparing equilibrated, isotropic long-chain model polymer melts. As an application, we generate Kremer-Grest melts of 1000 chains with 200 entanglements and 25 000-2000 beads/chain, which cover the experimentally relevant bending rigidities up to and beyond the limit of the isotropic-nematic transition. In the first step, we employ Monte Carlo simulations of a lattice model to equilibrate the large-scale chain structure above the tube scale while ensuring a spatially homogeneous density distribution. We then use theoretical insight from a constrained mode tube model to introduce the bead degrees of freedom together with random walk conformational statistics all the way down to the Kuhn scale of the chains. This is followed by a sequence of simulations with carefully parameterized force-capped bead-spring models, which slowly introduce the local bead packing while reproducing the larger-scale chain statistics of the target Kremer-Grest system at all levels of force-capping. Finally, we can switch to the full Kremer-Grest model without perturbing the structure. The resulting chain statistics is in excellent agreement with literature results on all length scales accessible in brute-force simulations of shorter chains.

Temporally Coherent Backmapping of Molecular Trajectories From Coarse-Grained to Atomistic Resolution

Coarse-graining offers a means to extend the achievable time and length scales of molecular dynamics simulations beyond what is practically possible in the atomistic regime. Sampling molecular configurations of interest can be done efficiently using coarse-grained simulations, from which meaningful physicochemical information can be inferred if the corresponding all-atom configurations are reconstructed. However, this procedure of backmapping to reintroduce the lost atomistic detail into coarse-grain structures has proven a challenging task due to the many feasible atomistic configurations that can be associated with one coarse-grain structure. Existing backmapping methods are strictly frame-based, relying on either heuristics to replace coarse-grain particles with atomic fragments and subsequent relaxation or parametrized models to propose atomic coordinates separately and independently for each coarse-grain structure. These approaches neglect information from previous trajectory frames that is critical to ensuring temporal coherence of the backmapped trajectory, while also offering information potentially helpful to producing higher-fidelity atomic reconstructions. In this work, we present a deep learning-enabled data-driven approach for temporally coherent backmapping that explicitly incorporates information from preceding trajectory structures. Our method trains a conditional variational autoencoder to nondeterministically reconstruct atomistic detail conditioned on both the target coarse-grain configuration and the previously reconstructed atomistic configuration. We demonstrate our backmapping approach on two exemplar biomolecular systems: alanine dipeptide and the miniprotein chignolin. We show that our backmapped trajectories accurately recover the structural, thermodynamic, and kinetic properties of the atomistic trajectory data.

Coarse-Grained Modeling of Ion-Containing Polymers

Ion-containing polymers have continued to be an important research focus for several decades due to their use as an electrolyte in energy storage and conversion devices. Elucidation of connections between the mesoscopic structure and multiscale dynamics of the ions and solvent remains incompletely understood. Coarse-grained modeling provides an efficient approach for exploring the structural and dynamical properties of these soft materials. The unique physicochemical properties of such polymers are of broad interest. In this review, we summarize the current development and understanding of the structure-property relationship of ion-containing polymers and provide insights into the design of such materials determined from coarse-grained modeling and simulations accompanying significant advances in experimental strategies. We specifically concentrate on three types of ion-containing polymers: proton exchange membranes (PEMs), anion exchange membranes (AEMs), and polymerized ionic liquids (polyILs). We posit that insight into the similarities and differences in these materials will lead to guidance in the rational design of high-performance novel materials with improved properties for various power source technologies.

Accelerated Simulations of Molecular Systems through Learning of Effective Dynamics

Simulations are vital for understanding and predicting the evolution of complex molecular systems. However, despite advances in algorithms and special purpose hardware, accessing the time scales necessary to capture the structural evolution of biomolecules remains a daunting task. In this work, we present a novel framework to advance simulation time scales by up to 3 orders of magnitude by learning the effective dynamics (LED) of molecular systems. LED augments the equation-free methodology by employing a probabilistic mapping between coarse and fine scales using mixture density network (MDN) autoencoders and evolves the non-Markovian latent dynamics using long short-term memory MDNs. We demonstrate the effectiveness of LED in the Müller-Brown potential, the Trp cage protein, and the alanine dipeptide. LED identifies explainable reduced-order representations, i.e., collective variables, and can generate, at any instant, all-atom molecular trajectories consistent with the collective variables. We believe that the proposed framework provides a dramatic increase to simulation capabilities and opens new horizons for the effective modeling of complex molecular systems.

Neural Upscaling from Residue-Level Protein Structure Networks to Atomistic Structures

Coarse-graining is a powerful tool for extending the reach of dynamic models of proteins and other biological macromolecules. Topological coarse-graining, in which biomolecules or sets thereof are represented via graph structures, is a particularly useful way of obtaining highly compressed representations of molecular structures, and simulations operating via such representations can achieve substantial computational savings. A drawback of coarse-graining, however, is the loss of atomistic detail—an effect that is especially acute for topological representations such as protein structure networks (PSNs). Here, we introduce an approach based on a combination of machine learning and physically-guided refinement for inferring atomic coordinates from PSNs. This “neural upscaling” procedure exploits the constraints implied by PSNs on possible configurations, as well as differences in the likelihood of observing different configurations with the same PSN. Using a 1 μs atomistic molecular dynamics trajectory of Aβ1–40, we show that neural upscaling is able to effectively recapitulate detailed structural information for intrinsically disordered proteins, being particularly successful in recovering features such as transient secondary structure. These results suggest that scalable network-based models for protein structure and dynamics may be used in settings where atomistic detail is desired, with upscaling employed to impute atomic coordinates from PSNs.

Solvent-scaling as an alternative to coarse-graining in adaptive-resolution simulations: The adaptive solvent-scaling (AdSoS) scheme

A new approach termed Adaptive Solvent-Scaling (AdSoS) is introduced for performing simulations of a solute embedded in a fine-grained (FG) solvent region itself surrounded by a coarse-grained (CG) solvent region, with a continuous FG ↔ CG switching of the solvent resolution across a buffer layer. Instead of relying on a distinct CG solvent model, the AdSoS scheme is based on CG models defined by a dimensional scaling of the FG solvent by a factor s, accompanied by an s-dependent modulation of the atomic masses and interaction parameters. The latter changes are designed to achieve an isomorphism between the dynamics of the FG and CG models, and to preserve the dispersive and dielectric solvation properties of the solvent with respect to a solute at FG resolution. This scaling approach offers a number of advantages compared to traditional coarse-graining: (i) the CG parameters are immediately related to those of the FG model (no need to parameterize a distinct CG model); (ii) nearly ideal mixing is expected for CG variants with similar s-values (ideal mixing holding in the limit of identical s-values); (iii) the solvent relaxation timescales should be preserved (no dynamical acceleration typical for coarse-graining); (iv) the graining level N_G (number of FG molecules represented by one CG molecule) can be chosen arbitrarily (in particular, N_G = s³ is not necessarily an integer); and (v) in an adaptive-resolution scheme, this level can be varied continuously as a function of the position (without requiring a bundling mechanism), and this variation occurs at a constant number of particles per molecule (no occurrence of fractional degrees of freedom in the buffer layer). By construction, the AdSoS scheme minimizes the thermodynamic mismatch between the different regions of the adaptive-resolution system, leading to a nearly homogeneous scaled solvent density s³ρ. Residual density artifacts in and at the surface of the boundary layer can easily be corrected by means of a grid-based biasing potential constructed in a preliminary pure-solvent simulation. This article introduces the AdSoS scheme and provides an initial application to pure atomic liquids (no solute) with Lennard-Jones plus Coulomb interactions in a slab geometry.

Introducing Memory in Coarse-Grained Molecular Simulations

Preserving the correct dynamics at the coarse-grained (CG) level is a pressing problem in the development of systematic CG models in soft matter simulation. Starting from the seminal idea of simple time-scale mapping, there have been many efforts over the years toward establishing a meticulous connection between the CG and fine-grained (FG) dynamics based on fundamental statistical mechanics approaches. One of the most successful attempts in this context has been the development of CG models based on the Mori-Zwanzig (MZ) theory, where the resulting equation of motion has the form of a generalized Langevin equation (GLE) and closely preserves the underlying FG dynamics. In this Review, we describe some of the recent studies in this regard. We focus on the construction and simulation of dynamically consistent systematic CG models based on the GLE, both in the simple Markovian limit and the non-Markovian case. Some recent studies of physical effects of memory are also discussed. The Review is aimed at summarizing recent developments in the field while highlighting the major challenges and possible future directions.

Atomic-level reconstruction of biomolecules by a rigid-fragment- and local-frame-based (RF-LF) strategy

Coarse-grained (CG) model has been a powerful tool in bridging the gap between theoretical studies and experimental phenomena in biological computing field. The reconstruction from a CG model to an atomic-detail structure is especially important in CG studies of biological systems. In this work, a rigid-fragment- and local-frame-based (RF-LF) backmapping method was proposed to achieve reverse mapping from CG models to atomic-level structures. The initial atomic-level structures were further refined to yield the final backmapping ones. With the popular Martini force field, the performance of the RF-LF method was extensively examined in the CG → AA (CG to AA) backmapping of protein/DNA/RNA systems. Besides, the RF-LF method was also extended to the backmapping of the TMFF model. Numerical results illustrate that the RF-LF backmapping method is generic and parameter-free and can provide a promising way to tackle atomic-level studies in CG models.

Molecular Modeling Investigations of Sorption and Diffusion of Small Molecules in Glassy Polymers

With a wide range of applications, from energy and environmental engineering, such as in gas separations and water purification, to biomedical engineering and packaging, glassy polymeric materials remain in the core of novel membrane and state-of the art barrier technologies. This review focuses on molecular simulation methodologies implemented for the study of sorption and diffusion of small molecules in dense glassy polymeric systems. Basic concepts are introduced and systematic methods for the generation of realistic polymer configurations are briefly presented. Challenges related to the long length and time scale phenomena that govern the permeation process in the glassy polymer matrix are described and molecular simulation approaches developed to address the multiscale problem at hand are discussed.

A coarse-grain model for entangled polyethylene melts and polyethylene crystallization

The Shinoda-DeVane-Klein (SDK) model is herein demonstrated to be a viable coarse-grain model for performing molecular simulations of polyethylene (PE), affording new opportunities to advance molecular-level, scientific understanding of PE materials and processes. Both structural and dynamical properties of entangled PE melts are captured by the SDK model, which also recovers important aspects of PE crystallization phenomenology. Importantly, the SDK model can be used to represent a variety of materials beyond PE and has a simple functional form, making it unique among coarse-grain PE models. This study expands the suite of tools for studying PE in silico and paves the way for future work probing PE and PE-based composites at the molecular level.

Backmapping from Multiresolution Coarse-Grained Models to Atomic Structures of Large Biomolecules by Restrained Molecular Dynamics Simulations Using Bayesian Inference

Coarse-grained (CG) simulations have allowed access to larger length scales and longer time scales in the study of the dynamic processes of large biomolecules than all-atom (AA) molecular dynamics (MD) simulations. Backmapping from CG models to AA structures has long been studied because it enables us to gain detailed structure insights from CG simulations. Many methods first construct an AA structure from the CG model by fragments, random placement, or geometrical rules and subsequently optimize the solution via energy minimization, simulated annealing or position-restrained simulations. However, such methods may only work well on residue-level CG models and cannot consider the deviations of CG models. In this work, we describe, to the best of our knowledge, a new backmapping method based on Bayesian inference and restrained MD simulations. Restraints with log harmonic energy terms are defined according to the target CG model using the Bayesian inference in which the CG deviations can be estimated. From an initial AA structure obtained from either high-resolution experiments or homology modeling, a MD simulation with the aforementioned restraints is performed to obtain a final AA structure that is a backmapping of the target CG model. The method was validated using multiresolution CG models of the soluble extracellular region of the human epidermal growth factor receptor and was further applied to construct AA structures from CG simulations of the nucleosome core particle. The results demonstrate that our method can generate accurate AA structures of different types of biomolecules from multiple CG models with either residue-level resolution or much lower resolution than one-site-per-residue.

Energy renormalization for coarse-graining polymers having different segmental structures

Multiscale coarse-grained (CG) modeling of soft materials, such as polymers, is currently an art form because CG models normally have significantly altered dynamics and thermodynamic properties compared to their atomistic counterparts. We address this problem by exploiting concepts derived from the generalized entropy theory (GET), emphasizing the central role of configurational entropy s c in the dynamics of complex fluids. Our energy renormalization (ER) method involves varying the cohesive interaction strength in the CG models in such a way that dynamic properties related to s c are preserved. We test this ER method by applying it to coarse-graining polymer melts (i.e., polybutadiene, polystyrene, and polycarbonate), representing polymer materials having a relatively low, intermediate, and high degree of glass "fragility". We find that the ER method allows the dynamics of the atomistic polymer models to be faithfully described to a good approximation by CG models over a wide temperature range.

Polymer nucleation under high-driving force, long-chain conditions: Heat release and the separation of time scales

This study reveals important features of polymer crystal formation at high-driving forces in entangled polymer melts based on simulations of polyethylene. First and in contrast to small-molecule crystallization, the heat released during polymer crystallization does not appreciably influence structural details of early-stage, crystalline clusters (crystal nuclei). Second, early-stage polymer crystallization (crystal nucleation) can occur without substantial chain-level relaxation and conformational changes. This study's results indicate that local structures and environments guide crystal nucleation in entangled polymer melts under high-driving force conditions. Given that such conditions are often used to process polyethylene, local structures and the separation of time scales associated with crystallization and chain-level processes are anticipated to be of substantial importance to processing strategies. This study highlights new research directions for understanding polymer crystallization.

Hierarchical modelling of polystyrene melts: from soft blobs to atomistic resolution

We demonstrate that hierarchical backmapping strategies incorporating generic blob-based models can equilibrate melts of high-molecular-weight polymers, described with chemically specific, atomistic models. The central idea is first to represent polymers by chains of large soft blobs (spheres) and efficiently equilibrate the melt on large scales. Then, the degrees of freedom of more detailed models are reinserted step by step. The procedure terminates when the atomistic description is reached. Reinsertions are feasible computationally because the fine-grained melt must be re-equilibrated only locally. We consider polystyrene (PS) which is sufficiently complex to serve method development because of stereo-chemistry and bulky side groups. Our backmapping strategy bridges mesoscopic and atomistic scales by incorporating a blob-based, a moderately coarse-grained (CG), and a united-atom model of PS. We demonstrate that the generic blob-based model can be parameterised to reproduce the mesoscale properties of a specific polymer - here PS. The moderately CG model captures stereo-chemistry. To perform backmapping we improve and adjust several fine-graining techniques. We prove equilibration of backmapped PS melts by comparing their structural and conformational properties with reference data from smaller systems, equilibrated with less efficient methods.

Single molecule translocation in smectics illustrates the challenge for time-mapping in simulations on multiple scales

Understanding the connections between the characteristic dynamical time scales associated with a coarse-grained (CG) and a detailed representation is central to the applicability of the coarse-graining methods to understand molecular processes. The process of coarse graining leads to an accelerated dynamics, owing to the smoothening of the underlying free-energy landscapes. Often a single time-mapping factor is used to relate the time scales associated with the two representations. We critically examine this idea using a model system ideally suited for this purpose. Single molecular transport properties are studied via molecular dynamics simulations of the CG and atomistic representations of a liquid crystalline, azobenzene containing mesogen, simulated in the smectic and the isotropic phases. The out-of-plane dynamics in the smectic phase occurs via molecular hops from one smectic layer to the next. Hopping can occur via two mechanisms, with and without significant reorientation. The out-of-plane transport can be understood as a superposition of two (one associated with each mode of transport) independent continuous time random walks for which a single time-mapping factor would be rather inadequate. A comparison of the free-energy surfaces, relevant to the out-of-plane transport, qualitatively supports the above observations. Thus, this work underlines the need for building CG models that exhibit both structural and dynamical consistency to the underlying atomistic model.

Adaptive resolution simulations of biomolecular systems

In this review article, we discuss and analyze some recently developed hybrid atomistic-mesoscopic solvent models for multiscale biomolecular simulations. We focus on the biomolecular applications of the adaptive resolution scheme (AdResS), which allows solvent molecules to change their resolution back and forth between atomistic and coarse-grained representations according to their positions in the system. First, we discuss coupling of atomistic and coarse-grained models of salt solution using a 1-to-1 molecular mapping-i.e., one coarse-grained bead represents one water molecule-for development of a multiscale salt solution model. In order to make use of coarse-grained molecular models that are compatible with the MARTINI force field, one has to resort to a supramolecular mapping, in particular to a 4-to-1 mapping, where four water molecules are represented with one coarse-grained bead. To this end, bundled atomistic water models are employed, i.e., the relative movement of water molecules that are mapped to the same coarse-grained bead is restricted by employing harmonic springs. Supramolecular coupling has recently also been extended to polarizable coarse-grained water models with explicit charges. Since these coarse-grained models consist of several interaction sites, orientational degrees of freedom of the atomistic and coarse-grained representations are coupled via a harmonic energy penalty term. The latter aligns the dipole moments of both representations. The reviewed multiscale solvent models are ready to be used in biomolecular simulations, as illustrated in a few examples.

A Review of Multiscale Computational Methods in Polymeric Materials

Polymeric materials display distinguished characteristics which stem from the interplay of phenomena at various length and time scales. Further development of polymer systems critically relies on a comprehensive understanding of the fundamentals of their hierarchical structure and behaviors. As such, the inherent multiscale nature of polymer systems is only reflected by a multiscale analysis which accounts for all important mechanisms. Since multiscale modelling is a rapidly growing multidisciplinary field, the emerging possibilities and challenges can be of a truly diverse nature. The present review attempts to provide a rather comprehensive overview of the recent developments in the field of multiscale modelling and simulation of polymeric materials. In order to understand the characteristics of the building blocks of multiscale methods, first a brief review of some significant computational methods at individual length and time scales is provided. These methods cover quantum mechanical scale, atomistic domain (Monte Carlo and molecular dynamics), mesoscopic scale (Brownian dynamics, dissipative particle dynamics, and lattice Boltzmann method), and finally macroscopic realm (finite element and volume methods). Afterwards, different prescriptions to envelope these methods in a multiscale strategy are discussed in details. Sequential, concurrent, and adaptive resolution schemes are presented along with the latest updates and ongoing challenges in research. In sequential methods, various systematic coarse-graining and backmapping approaches are addressed. For the concurrent strategy, we aimed to introduce the fundamentals and significant methods including the handshaking concept, energy-based, and force-based coupling approaches. Although such methods are very popular in metals and carbon nanomaterials, their use in polymeric materials is still limited. We have illustrated their applications in polymer science by several examples hoping for raising attention towards the existing possibilities. The relatively new adaptive resolution schemes are then covered including their advantages and shortcomings. Finally, some novel ideas in order to extend the reaches of atomistic techniques are reviewed. We conclude the review by outlining the existing challenges and possibilities for future research.

The Perfect Glass Paradigm: Disordered Hyperuniform Glasses Down to Absolute Zero

Rapid cooling of liquids below a certain temperature range can result in a transition to glassy states. The traditional understanding of glasses includes their thermodynamic metastability with respect to crystals. However, here we present specific examples of interactions that eliminate the possibilities of crystalline and quasicrystalline phases, while creating mechanically stable amorphous glasses down to absolute zero temperature. We show that this can be accomplished by introducing a new ideal state of matter called a "perfect glass". A perfect glass represents a soft-interaction analog of the maximally random jammed (MRJ) packings of hard particles. These latter states can be regarded as the epitome of a glass since they are out of equilibrium, maximally disordered, hyperuniform, mechanically rigid with infinite bulk and shear moduli, and can never crystallize due to configuration-space trapping. Our model perfect glass utilizes two-, three-, and four-body soft interactions while simultaneously retaining the salient attributes of the MRJ state. These models constitute a theoretical proof of concept for perfect glasses and broaden our fundamental understanding of glass physics. A novel feature of equilibrium systems of identical particles interacting with the perfect-glass potential at positive temperature is that they have a non-relativistic speed of sound that is infinite.

Protein Aggregation and Molecular Crowding: Perspectives From Multiscale Simulations

Cells are extremely crowded environments, thus the use of diluted salted aqueous solutions containing a single protein is too simplistic to mimic the real situation. Macromolecular crowding might affect protein structure, folding, shape, conformational stability, binding of small molecules, enzymatic activity, interactions with cognate biomolecules, and pathological aggregation. The latter phenomenon typically leads to the formation of amyloid fibrils that are linked to several lethal neurodegenerative diseases, but that can also play a functional role in certain organisms. The majority of molecular simulations performed before the last few years were conducted in diluted solutions and were restricted both in the timescales and in the system dimensions by the available computational resources. In recent years, several computational solutions were developed to get close to physiological conditions. In this review we summarize the main computational techniques used to tackle the issue of protein aggregation both in a diluted and in a crowded environment.

Multiscale approach for the construction of equilibrated all-atom models of a poly(ethylene glycol)-based hydrogel

A multiscale modeling approach is presented for the efficient construction of an equilibrated all-atom model of a cross-linked poly(ethylene glycol) (PEG)-based hydrogel using the all-atom polymer consistent force field (PCFF). The final equilibrated all-atom model was built with a systematic simulation toolset consisting of three consecutive parts: (1) building a global cross-linked PEG-chain network at experimentally determined cross-link density using an on-lattice Monte Carlo method based on the bond fluctuation model, (2) recovering the local molecular structure of the network by transitioning from the lattice model to an off-lattice coarse-grained (CG) model parameterized from PCFF, followed by equilibration using high performance molecular dynamics methods, and (3) recovering the atomistic structure of the network by reverse mapping from the equilibrated CG structure, hydrating the structure with explicitly represented water, followed by final equilibration using PCFF parameterization. The developed three-stage modeling approach has application to a wide range of other complex macromolecular hydrogel systems, including the integration of peptide, protein, and/or drug molecules as side-chains within the hydrogel network for the incorporation of bioactivity for tissue engineering, regenerative medicine, and drug delivery applications.

Resolving Dynamic Properties of Polymers through Coarse-Grained Computational Studies

Coupled length and time scales determine the dynamic behavior of polymers and underlie their unique viscoelastic properties. To resolve the long-time dynamics it is imperative to determine which time and length scales must be correctly modeled. Here we probe the degree of coarse graining required to simultaneously retain significant atomistic details and access large length and time scales. The degree of coarse graining in turn sets the minimum length scale instrumental in defining polymer properties and dynamics. Using linear polyethylene as a model system, we probe how the coarse-graining scale affects the measured dynamics. Iterative Boltzmann inversion is used to derive coarse-grained potentials with 2-6 methylene groups per coarse-grained bead from a fully atomistic melt simulation. We show that atomistic detail is critical to capturing large-scale dynamics. Using these models we simulate polyethylene melts for times over 500  μs to study the viscoelastic properties of well-entangled polymer melts.

Simulations of inorganic-bioorganic interfaces to discover new materials: insights, comparisons to experiment, challenges, and opportunities

Natural and man-made materials often rely on functional interfaces between inorganic and organic compounds. Examples include skeletal tissues and biominerals, drug delivery systems, catalysts, sensors, separation media, energy conversion devices, and polymer nanocomposites. Current laboratory techniques are limited to monitor and manipulate assembly on the 1 to 100 nm scale, time-consuming, and costly. Computational methods have become increasingly reliable to understand materials assembly and performance. This review explores the merit of simulations in comparison to experiment at the 1 to 100 nm scale, including connections to smaller length scales of quantum mechanics and larger length scales of coarse-grain models. First, current simulation methods, advances in the understanding of chemical bonding, in the development of force fields, and in the development of chemically realistic models are described. Then, the recognition mechanisms of biomolecules on nanostructured metals, semimetals, oxides, phosphates, carbonates, sulfides, and other inorganic materials are explained, including extensive comparisons between modeling and laboratory measurements. Depending on the substrate, the role of soft epitaxial binding mechanisms, ion pairing, hydrogen bonds, hydrophobic interactions, and conformation effects is described. Applications of the knowledge from simulation to predict binding of ligands and drug molecules to the inorganic surfaces, crystal growth and shape development, catalyst performance, as well as electrical properties at interfaces are examined. The quality of estimates from molecular dynamics and Monte Carlo simulations is validated in comparison to measurements and design rules described where available. The review further describes applications of simulation methods to polymer composite materials, surface modification of nanofillers, and interfacial interactions in building materials. The complexity of functional multiphase materials creates opportunities to further develop accurate force fields, including reactive force fields, and chemically realistic surface models, to enable materials discovery at a million times lower computational cost compared to quantum mechanical methods. The impact of modeling and simulation could further be increased by the advancement of a uniform simulation platform for organic and inorganic compounds across the periodic table and new simulation methods to evaluate system performance in silico.

Transferability of a coarse-grained atactic polystyrene model: the non-bonded potential effect

In this paper, we construct an efficient and simple coarse grained (CG) model for atactic polystyrene (PS) by using a 1 : 1 mapping scheme at 463 K and 1 atm pressure and derive the corresponding bonded and non-bonded potentials in the CG force field (FF) via a direct Boltzmann inversion approach and a combined structure-based and thermodynamic quantities-based CG method, respectively.

Flexible Boundaries for Multiresolution Solvation: An Algorithm for Spatial Multiscaling in Molecular Dynamics Simulations

An algorithm is proposed for performing molecular dynamics (MD) simulations of a biomolecular solute represented at atomistic resolution surrounded by a surface layer of atomistic fine-grained (FG) solvent molecules within a bulk represented by coarse-grained (CG) solvent beads. The method, called flexible boundaries for multiresolution solvation (FBMS), is based on: (i) a three-region layering of the solvent around the solute, involving an FG layer surrounded by a mixed FG-CG buffer layer, itself surrounded by a bulk CG region; (ii) a definition of the layer boundary that relies on an effective distance to the solute surface and is thus adapted to the shape of the solute as well as adjusts to its conformational changes. The effective surface distance is defined by inverse-nth power averaging over the distances to all non-hydrogen solute atoms, and the layering is enforced by means of half-harmonic distance restraints, attractive for the FG molecules and repulsive for the CG beads. A restraint-free region at intermediate distances enables the formation of the buffer layer, where the FG and CG solvents can mix freely. The algorithm is tested and validated using the GROMOS force field and the associated FG (SPC) and CG (polarizable CGW) water models. The test systems include pure-water systems, where one FG molecule plays the role of a solute, and a deca-alanine peptide with two widely different solute shapes considered, α-helical and fully extended. In particular, as the peptide unfolds, the number of FG molecules required to fill its close-range solvation layer increases, with the additional molecules being provided by the buffer layer. Further validation involves simulations of four proteins in multiresolution FG/CG mixtures. The resulting structural, energetic, and solvation properties are found to be similar to those observed in corresponding pure FG simulations.

Biomembranes in atomistic and coarse-grained simulations

The architecture of biological membranes is tightly coupled to the localization, organization, and function of membrane proteins. The organelle-specific distribution of lipids allows for the formation of functional microdomains (also called rafts) that facilitate the segregation and aggregation of membrane proteins and thus shape their function. Molecular dynamics simulations enable to directly access the formation, structure, and dynamics of membrane microdomains at the molecular scale and the specific interactions among lipids and proteins on timescales from picoseconds to microseconds. This review focuses on the latest developments of biomembrane force fields for both atomistic and coarse-grained molecular dynamics (MD) simulations, and the different levels of coarsening of biomolecular structures. It also briefly introduces scale-bridging methods applicable to biomembrane studies, and highlights selected recent applications.

Going Backward: A Flexible Geometric Approach to Reverse Transformation from Coarse Grained to Atomistic Models

The conversion of coarse-grained to atomistic models is an important step in obtaining insight about atomistic scale processes from coarse-grained simulations. For this process, called backmapping or reverse transformation, several tools are available, but these commonly require libraries of molecule fragments or they are linked to a specific software package. In addition, the methods are usually restricted to specific molecules and to a specific force field. Here, we present an alternative method, consisting of geometric projection and subsequent force-field based relaxation. This method is designed to be simple and flexible, and offers a generic solution for resolution transformation. For simple systems, the conversion only requires a list of particle correspondences on the two levels of resolution. For special cases, such as nondefault protonation states of amino acids and virtual sites, a target particle list can be specified. The mapping uses simple building blocks, which list the particles on the different levels of resolution. For conversion to higher resolution, the initial model is relaxed with several short cycles of energy minimization and position-restrained MD. The reconstruction of an atomistic backbone from a coarse-grained model is done using a new dedicated algorithm. The method is generic and can be used to map between any two particle based representations, provided that a mapping can be written. The focus of this work is on the coarse-grained MARTINI force field, for which mapping definitions are written to allow conversion to and from the higher-resolution force fields GROMOS, CHARMM, and AMBER, and to and from a simplified three-bead lipid model. Together, these offer the possibility to simulate mesoscopic membrane structures, to be transformed to MARTINI and subsequently to an atomistic model for investigation of detailed interactions. The method was tested on a set of systems ranging from a simple, single-component bilayer to a large protein-membrane-solvent complex. The results demonstrate the efficiency and the efficacy of the new approach.