GPU-Enhanced DFTB Metadynamics for Efficiently Predicting Free Energies of Biochemical Systems

24 January 2023, Version 1
This content is a preprint and has not undergone peer review at the time of posting.


Carrying out metadynamics calculations of large chemical systems with ab initio methods is computationally prohibitive due to the extensive sampling required to simulate the large degrees of freedom in these systems. To address this computational bottleneck, we utilize a GPU-enhanced density functional tight binding (DFTB) approach on a massively-parallelized cloud computing platform to efficiently calculate the thermodynamics and metadynamics of biochemical systems. To first validate our approach, we calculate free energy surfaces of alanine dipeptide and show that our GPU-enhanced DFTB calculations qualitatively agree with computationally-intensive hybrid DFT benchmarks, whereas classical force fields give significant errors. Most importantly, we show that our GPU-accelerated DFTB calculations are significantly faster than previous approaches by up to 2 orders of magnitude. To further extend our GPU-enhanced DFTB approach, we also carried out a 10-ns metadynamics simulation of remdesivir, which is prohibitively out of reach for routine DFT-based metadynamics calculations. We find that the free energy surfaces of remdesivir obtained from DFTB and classical force fields differ significantly, where the latter overestimates the internal energy contribution of high free energy states. Taken together, our benchmark tests, analyses, and extensions to large biochemical systems highlight the use of GPU-enhanced DFTB simulations for efficiently predicting the free energy surfaces/thermodynamics of large biochemical systems.


density functional tight binding
density functional theory
free energies
biochemical systems
ab initio calculations

Supplementary materials

Supporting Information
Additional comparisons of the entropic term for alanine dipeptide and remdesivir obtained from Amber-ff19SB and SCC-DFTB3 metadynamics calculations, convergence analyses of the free energy surface of remdesivir, and hybrid DFT benchmark calculations on selected remdesivir geometries.

Supplementary weblinks


Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.