A Kepler Workflow Tool for Reproducible AMBER GPU Molecular Dynamics.

With the drive toward high throughput molecular dynamics (MD) simulations involving ever-greater numbers of simulation replicates run for longer, biologically relevant timescales (microseconds), the need for improved computational methods that facilitate fully automated MD workflows gains more importance. Here we report the development of an automated workflow tool to perform AMBER GPU MD simulations. Our workflow tool capitalizes on the capabilities of the Kepler platform to deliver a flexible, intuitive, and user-friendly environment and the AMBER GPU code for a robust and high-performance simulation engine. Additionally, the workflow tool reduces user input time by automating repetitive processes and facilitates access to GPU clusters, whose high-performance processing power makes simulations of large numerical scale possible. The presented workflow tool facilitates the management and deployment of large sets of MD simulations on heterogeneous computing resources. The workflow tool also performs systematic analysis on the simulation outputs and enhances simulation reproducibility, execution scalability, and MD method development including benchmarking and validation.

[1]  Frank Noé,et al.  PyEMMA 2: A Software Package for Estimation, Validation, and Analysis of Markov Models. , 2015, Journal of chemical theory and computation.

[2]  S. Chung,et al.  Molecular dynamics study of the KcsA potassium channel. , 1999, Biophysical journal.

[3]  A. Levine,et al.  Surfing the p53 network , 2000, Nature.

[4]  D. Green,et al.  Cytoplasmic functions of the tumour suppressor p53 , 2009, Nature.

[5]  Peter M. Kasson,et al.  GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit , 2013, Bioinform..

[6]  Daniel R Roe,et al.  PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data. , 2013, Journal of chemical theory and computation.

[7]  C. Harris,et al.  The IARC TP53 database: New online mutation analysis and recommendations to users , 2002, Human mutation.

[8]  Ilkay Altintas,et al.  Provenance Collection Support in the Kepler Scientific Workflow System , 2006, IPAW.

[9]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[10]  Klaus Schulten,et al.  QwikMD — Integrative Molecular Dynamics Toolkit for Novices and Experts , 2016, Scientific Reports.

[11]  Duncan Poole,et al.  Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 2. Explicit Solvent Particle Mesh Ewald. , 2013, Journal of chemical theory and computation.

[12]  Bertram Ludäscher,et al.  Kepler: an extensible system for design and execution of scientific workflows , 2004 .

[13]  Federico D. Sacerdoti,et al.  Scalable Algorithms for Molecular Dynamics Simulations on Commodity Clusters , 2006, ACM/IEEE SC 2006 Conference (SC'06).

[14]  Mohammad M. Sultan,et al.  MSMBuilder: Statistical Models for Biomolecular Dynamics , 2016, bioRxiv.

[15]  Duncan Poole,et al.  Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 1. Generalized Born , 2012, Journal of chemical theory and computation.

[16]  Frank Noé,et al.  HTMD: High-Throughput Molecular Dynamics for Molecular Discovery. , 2016, Journal of chemical theory and computation.

[17]  Alexander D. MacKerell,et al.  CHARMM‐GUI 10 years for biomolecular modeling and simulation , 2017, J. Comput. Chem..

[18]  Ilkay Altintas,et al.  Biomedical Big Data Training Collaborative (BBDTC): An effort to bridge the talent gap in biomedical science and research , 2017, J. Comput. Sci..

[19]  Ilkay Altintas,et al.  Biomedical Big Data Training Collaborative (BBDTC): An Effort to Bridge the Talent Gap in Biomedical Science and Research , 2016, ICCS.

[20]  Rommie E. Amaro,et al.  Computational identification of a transiently open L1/S3 pocket for reactivation of mutant p53 , 2013, Nature Communications.

[21]  Rommie E. Amaro,et al.  Molecular dynamics analysis of antibody recognition and escape by human H1N1 influenza hemagglutinin. , 2015, Biophysical journal.

[22]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[23]  Nivedita Jaiswal,et al.  Molecular docking and dynamics simulation analyses unraveling the differential enzymatic catalysis by plant and fungal laccases with respect to lignin biosynthesis and degradation , 2015, Journal of biomolecular structure & dynamics.

[24]  J A McCammon,et al.  Analysis of a 10-ns molecular dynamics simulation of mouse acetylcholinesterase. , 2001, Biophysical journal.

[25]  Wei Gu,et al.  Tumor Suppression in the Absence of p53-Mediated Cell-Cycle Arrest, Apoptosis, and Senescence , 2012, Cell.

[26]  Diwakar Shukla,et al.  OpenMM 4: A Reusable, Extensible, Hardware Independent Library for High Performance Molecular Simulation. , 2013, Journal of chemical theory and computation.

[27]  Andreas W. Götz,et al.  SPFP: Speed without compromise - A mixed precision model for GPU accelerated molecular dynamics simulations , 2013, Comput. Phys. Commun..