scholarly journals Coarse-graining molecular dynamics: stochastic models with non-Gaussian force distributions

2019 ◽  
Vol 80 (1-2) ◽  
pp. 457-479 ◽  
Author(s):  
Radek Erban

Abstract Incorporating atomistic and molecular information into models of cellular behaviour is challenging because of a vast separation of spatial and temporal scales between processes happening at the atomic and cellular levels. Multiscale or multi-resolution methodologies address this difficulty by using molecular dynamics (MD) and coarse-grained models in different parts of the cell. Their applicability depends on the accuracy and properties of the coarse-grained model which approximates the detailed MD description. A family of stochastic coarse-grained (SCG) models, written as relatively low-dimensional systems of nonlinear stochastic differential equations, is presented. The nonlinear SCG model incorporates the non-Gaussian force distribution which is observed in MD simulations and which cannot be described by linear models. It is shown that the nonlinearities can be chosen in such a way that they do not complicate parametrization of the SCG description by detailed MD simulations. The solution of the SCG model is found in terms of gamma functions.


Author(s):  
Radek Erban

Molecular dynamics (MD) simulations of ions (K + , Na + , Ca 2+ and Cl − ) in aqueous solutions are investigated. Water is described using the SPC/E model. A stochastic coarse-grained description for ion behaviour is presented and parametrized using MD simulations. It is given as a system of coupled stochastic and ordinary differential equations, describing the ion position, velocity and acceleration. The stochastic coarse-grained model provides an intermediate description between all-atom MD simulations and Brownian dynamics (BD) models. It is used to develop a multiscale method which uses all-atom MD simulations in parts of the computational domain and (less detailed) BD simulations in the remainder of the domain.



2019 ◽  
Vol 5 (1) ◽  
Author(s):  
Wujie Wang ◽  
Rafael Gómez-Bombarelli

AbstractMolecular dynamics simulations provide theoretical insight into the microscopic behavior of condensed-phase materials and, as a predictive tool, enable computational design of new compounds. However, because of the large spatial and temporal scales of thermodynamic and kinetic phenomena in materials, atomistic simulations are often computationally infeasible. Coarse-graining methods allow larger systems to be simulated by reducing their dimensionality, propagating longer timesteps, and averaging out fast motions. Coarse-graining involves two coupled learning problems: defining the mapping from an all-atom representation to a reduced representation, and parameterizing a Hamiltonian over coarse-grained coordinates. We propose a generative modeling framework based on variational auto-encoders to unify the tasks of learning discrete coarse-grained variables, decoding back to atomistic detail, and parameterizing coarse-grained force fields. The framework is tested on a number of model systems including single molecules and bulk-phase periodic simulations.



2021 ◽  
Vol 12 (1) ◽  
Author(s):  
Stephan Thaler ◽  
Julija Zavadlav

AbstractIn molecular dynamics (MD), neural network (NN) potentials trained bottom-up on quantum mechanical data have seen tremendous success recently. Top-down approaches that learn NN potentials directly from experimental data have received less attention, typically facing numerical and computational challenges when backpropagating through MD simulations. We present the Differentiable Trajectory Reweighting (DiffTRe) method, which bypasses differentiation through the MD simulation for time-independent observables. Leveraging thermodynamic perturbation theory, we avoid exploding gradients and achieve around 2 orders of magnitude speed-up in gradient computation for top-down learning. We show effectiveness of DiffTRe in learning NN potentials for an atomistic model of diamond and a coarse-grained model of water based on diverse experimental observables including thermodynamic, structural and mechanical properties. Importantly, DiffTRe also generalizes bottom-up structural coarse-graining methods such as iterative Boltzmann inversion to arbitrary potentials. The presented method constitutes an important milestone towards enriching NN potentials with experimental data, particularly when accurate bottom-up data is unavailable.



2016 ◽  
Vol 110 (3) ◽  
pp. 323a
Author(s):  
Kento Inoue ◽  
Eiji Ymamoto ◽  
Daisuke Takaiwa ◽  
Kenji Yasuoka ◽  
Masuhiro Mikami


Molecules ◽  
2021 ◽  
Vol 26 (11) ◽  
pp. 3293
Author(s):  
Mateusz Zalewski ◽  
Sebastian Kmiecik ◽  
Michał Koliński

One of the major challenges in the computational prediction of protein–peptide complexes is the scoring of predicted models. Usually, it is very difficult to find the most accurate solutions out of the vast number of sometimes very different and potentially plausible predictions. In this work, we tested the protocol for Molecular Dynamics (MD)-based scoring of protein–peptide complex models obtained from coarse-grained (CG) docking simulations. In the first step of the scoring procedure, all models generated by CABS-dock were reconstructed starting from their original C-alpha trace representations to all-atom (AA) structures. The second step included geometry optimization of the reconstructed complexes followed by model scoring based on receptor–ligand interaction energy estimated from short MD simulations in explicit water. We used two well-known AA MD force fields, CHARMM and AMBER, and a CG MARTINI force field. Scoring results for 66 different protein–peptide complexes show that the proposed MD-based scoring approach can be used to identify protein–peptide models of high accuracy. The results also indicate that the scoring accuracy may be significantly affected by the quality of the reconstructed protein receptor structures.



2013 ◽  
Vol 12 (02) ◽  
pp. 1250111 ◽  
Author(s):  
HAILONG XU ◽  
QIUYU ZHANG ◽  
HEPENG ZHANG ◽  
BAOLIANG ZHANG ◽  
CHANGJIE YIN

Dissipative particle dynamics (DPD) was initially used to simulate the polystyrene/nanoparticle composite microspheres (PNCM) in this paper. The coarse graining model of PNCM was established. And the DPD parameterization of the model was represented in detail. The DPD repulsion parameters were calculated from the cohesive energy density which could be calculated by amorphous modules in Materials Studio. The equilibrium configuration of the simulated PNCM shows that the nanoparticles were actually "modified" with oleic acid and the modified nanoparticles were embedded in the bulk of polystyrene. As sodium dodecyl sulfate (SDS) was located in the interface between water and polystyrene, the hydrophilic head of SDS stretched into water while the hydrophobic tailed into polystyrene. All simulated phenomena were consistent with the experimental results in preparation of polystyrene/nanoparticles composite microspheres. The effect of surface modification of nanoparticles on its dispersion in polystyrene matrix was also studied by adjusting the interaction parameters between the OA and NP beads. The final results indicated that the nanoparticles removed from the core of composite microsphere to the surface with increase of a OA-NP . All the simulated results demonstrated that our coarse–grained model was reasonable.



Soft Matter ◽  
2018 ◽  
Vol 14 (15) ◽  
pp. 2796-2807 ◽  
Author(s):  
Andrea Catte ◽  
Mark R. Wilson ◽  
Martin Walker ◽  
Vasily S. Oganesyan

Antimicrobial action of a cationic peptide is modelled by large scale MD simulations.



Molecules ◽  
2020 ◽  
Vol 25 (19) ◽  
pp. 4413
Author(s):  
Giovanny Aguilera-Durán ◽  
Antonio Romo-Mancillas

Vitiligo is a hypopigmentary skin pathology resulting from the death of melanocytes due to the activity of CD8+ cytotoxic lymphocytes and overexpression of chemokines. These include CXCL9, CXCL10, and CXCL11 and its receptor CXCR3, both in peripheral cells of the immune system and in the skin of patients diagnosed with vitiligo. The three-dimensional structure of CXCR3 and CXCL9 has not been reported experimentally; thus, homology modeling and molecular dynamics could be useful for the study of this chemotaxis-promoter axis. In this work, a homology model of CXCR3 and CXCL9 and the structure of the CXCR3/Gαi/0βγ complex with post-translational modifications of CXCR3 are reported for the study of the interaction of chemokines with CXCR3 through all-atom (AA-MD) and coarse-grained molecular dynamics (CG-MD) simulations. AA-MD and CG-MD simulations showed the first activation step of the CXCR3 receptor with all chemokines and the second activation step in the CXCR3-CXCL10 complex through a decrease in the distance between the chemokine and the transmembrane region of CXCR3 and the separation of the βγ complex from the α subunit in the G-protein. Additionally, a general protein–ligand interaction model was calculated, based on known antagonists binding to CXCR3. These results contribute to understanding the activation mechanism of CXCR3 and the design of new molecules that inhibit chemokine binding or antagonize the receptor, provoking a decrease of chemotaxis caused by the CXCR3/chemokines axis.





Sign in / Sign up

Export Citation Format

Share Document