Reverse Coarse-Graining for Equation-Free Modeling: Application to Multiscale Molecular Dynamics

Temporally Coarse-Grained All-Atom Molecular Dynamics Achieved via Stochastic Padé Approximants

A Padé approximant scheme for realizing the discrete-time evolution of the state of a many-atom system is introduced. This temporal coarse-graining scheme accounts for the underlying Newtonian physics and avoids the need for construction of spatially coarse-grained variables. Newtonian physics is incorporated through short molecular dynamics simulations at the beginning of each of the large coarse-grained timesteps. The balance between stochastic and coherent dynamics expressed by many-atom systems is captured via incorporation of the Ito formula into a Padé approximant for the time dependence of individual atom positions over large timesteps. Since the time for a many-atom system to express a characteristic ensemble of atomic velocity fluctuations is typically short relative to the characteristic time of large-scale atomic displacements, a computationally efficient and accurate temporal coarse-graining of the atom-resolved Newtonian dynamics is formulated, denoted all-atom Padé-Ito molecular dynamics (APIMD). Evolution of the system over a time step much longer than that required for standard molecular dynamics (MD) is achieved via incorporation of information from the short MD simulations into a Padé approximant extrapolation in time. The extrapolated atomic configuration is subjected to energy minimization and, when needed, thermal equilibration so as to avoid occasional unphysical close encounters deriving from the Padé approximant extrapolation and to represent configurations appropriate for the temperature of interest. APIMD is implemented and tested via comparison with traditional MD simulations of five phenomena: (1) pertussis toxin subunit deformation, (2) structural transition in a T = 1 capsid-like structure of HPV16 L1 protein, (3) coalescence of argon nanodroplets, and structural transitions in dialanine in (4) vacuum, and (5) water. Accuracy of APIMD is demonstrated using semimicroscopic descriptors (rmsd, radius of gyration, residue-residue contact maps, and densities) and the free energy. Significant computational acceleration relative to traditional molecular dynamics is illustrated.

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 Structuring and Percolation Transition in Hydrated Sulfonated Poly(ether ether ketone) Membranes

The extent of phase separation and water percolation in sulfonated membranes are the key to their performance in fuel cells. Toward this, the effect of hydration on the morphology and transport characteristics of sulfonated poly(ether ether ketone), sPEEK, membrane is investigated using atomistic molecular dynamics simulation at various hydration levels(λ: number of water molecules per sulfonate group). The evolution of local morphology is investigated using structural correlations and minimum pair distances. Transport properties are probed using mean squared displacements and diffusion coefficients. The water-sulfonate interaction in sPEEK is found to be stronger than that in Nafion, as observed in experiments. Analyses indicate the presence of narrow connected path of water and hydronium at λ = 4 and large domains, spanning half the simulation box, at λ = 15. The behavior of membrane water remains far from bulk as indicated by its diffusion coefficient. The persistence of small isolated water clusters demonstrates the extent of phase separation in sPEEK to be lesser than that in Nafion. Analyses at molecular and collective levels suggest the occurrence of a percolation transition between λ = 8 and 10, which leads to a connected network of water channels in the membrane, thereby boosting the hydronium mobility.