1. Academic Validation
  2. Asynchronous Reciprocal Coupling of Martini 2.2 Coarse-Grained and CHARMM36 All-Atom Simulations in an Automated Multiscale Framework

Asynchronous Reciprocal Coupling of Martini 2.2 Coarse-Grained and CHARMM36 All-Atom Simulations in an Automated Multiscale Framework

  • J Chem Theory Comput. 2022 Aug 9;18(8):5025-5045. doi: 10.1021/acs.jctc.2c00168.
Cesar A López 1 Xiaohua Zhang 2 Fikret Aydin 2 Rebika Shrestha 3 Que N Van 3 Christopher B Stanley 4 Timothy S Carpenter 2 Kien Nguyen 1 Lara A Patel 1 5 De Chen 3 Violetta Burns 1 Nicolas W Hengartner 1 Tyler J E Reddy 1 Harsh Bhatia 6 Francesco Di Natale 6 Timothy H Tran 3 Albert H Chan 3 Dhirendra K Simanshu 3 Dwight V Nissley 3 Frederick H Streitz 2 Andrew G Stephen 3 Thomas J Turbyville 3 Felice C Lightstone 2 Sandrasegaram Gnanakaran 1 Helgi I Ingólfsson 2 Chris Neale 1
Affiliations

Affiliations

  • 1 Theoretical Biology and Biophysics, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States.
  • 2 Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States.
  • 3 NCI RAS Initiative, Cancer Research Technology Program, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21701, United States.
  • 4 Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37830, United States.
  • 5 Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States.
  • 6 Computing Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States.
Abstract

The appeal of multiscale modeling approaches is predicated on the promise of combinatorial synergy. However, this promise can only be realized when distinct scales are combined with reciprocal consistency. Here, we consider multiscale molecular dynamics (MD) simulations that combine the accuracy and macromolecular flexibility accessible to fixed-charge all-atom (AA) representations with the sampling speed accessible to reductive, coarse-grained (CG) representations. AA-to-CG conversions are relatively straightforward because deterministic routines with unique outcomes are achievable. Conversely, CG-to-AA conversions have many solutions due to a surge in the number of degrees of freedom. While automated tools for biomolecular CG-to-AA transformation exist, we find that one popular option, called Backward, is prone to stochastic failure and the AA models that it does generate frequently have compromised protein structure and incorrect stereochemistry. Although these shortcomings can likely be circumvented by human intervention in isolated instances, automated multiscale coupling requires reliable and robust scale conversion. Here, we detail an extension to Multiscale Machine-learned Modeling Infrastructure (MuMMI), including an improved CG-to-AA conversion tool called sinceCG. This tool is reliable (∼98% weakly correlated repeat success rate), automatable (no unrecoverable hangs), and yields AA models that generally preserve protein secondary structure and maintain correct stereochemistry. We describe how the MuMMI framework identifies CG system configurations of interest, converts them to AA representations, and simulates them at the AA scale while on-the-fly analyses provide feedback to update CG parameters. Application to systems containing the peripheral membrane protein Ras and proximal components of Raf kinase on complex eight-component lipid bilayers with ∼1.5 million atoms is discussed in the context of MuMMI.

Figures
Products