High-throughput all-atom molecular dynamics simulations using distributed computing

J Chem Inf Model. 2010 Mar 22;50(3):397-403. doi: 10.1021/ci900455r.

Abstract

Although molecular dynamics simulation methods are useful in the modeling of macromolecular systems, they remain computationally expensive, with production work requiring costly high-performance computing (HPC) resources. We review recent innovations in accelerating molecular dynamics on graphics processing units (GPUs), and we describe GPUGRID, a volunteer computing project that uses the GPU resources of nondedicated desktop and workstation computers. In particular, we demonstrate the capability of simulating thousands of all-atom molecular trajectories generated at an average of 20 ns/day each (for systems of approximately 30 000-80 000 atoms). In conjunction with a potential of mean force (PMF) protocol for computing binding free energies, we demonstrate the use of GPUGRID in the computation of accurate binding affinities of the Src SH2 domain/pYEEI ligand complex by reconstructing the PMF over 373 umbrella sampling windows of 55 ns each (20.5 mus of total data). We obtain a standard free energy of binding of -8.7 +/- 0.4 kcal/mol within 0.7 kcal/mol from experimental results. This infrastructure will provide the basis for a robust system for high-throughput accurate binding affinity prediction.

Publication types

  • Research Support, Non-U.S. Gov't
  • Research Support, U.S. Gov't, Non-P.H.S.

MeSH terms

  • Humans
  • Molecular Dynamics Simulation* / economics
  • Molecular Dynamics Simulation* / trends
  • Oligopeptides / chemistry
  • Oligopeptides / metabolism*
  • Protein Binding
  • Thermodynamics
  • src Homology Domains*

Substances

  • Oligopeptides
  • phosphotyrosyl-glutamyl-glutamyl-isoleucine