These are preliminary reports that have not been peer-reviewed. They should not be regarded as conclusive, guide clinical practice/health-related behavior, or be reported in news media as established information. For more information, please see our FAQs.
2 files

TorsionNet: A Deep Neural Network to Rapidly Predict Small Molecule Torsion Energy Profiles with the Accuracy of Quantum Mechanics

submitted on 23.12.2020, 14:31 and posted on 24.12.2020, 10:54 by Brajesh Rai, Vishnu Sresht, Qingyi Yang, Rayomond J. Unwalla, Meihua Tu, Alan M. Mathiowetz, Gregory A. Bakken

TorsionNet: A Deep Neural Network to Rapidly Predict Small Molecule Torsion Energy Profiles with the Accuracy of Quantum Mechanics

Brajesh K. Rai*,1, Vishnu Sresht1, Qingyi Yang2, Ray Unwalla2, Meihua Tu2, Alan M. Mathiowetz2, and Gregory A. Bakken3

1Simulation and Modeling Sciences and 2Medicine Design, Pfizer Worldwide Research Development and Medical, 610 Main Street, Cambridge, Massachusetts 02139, United States

3Digital, Pfizer, Eastern Point Road, Groton, Connecticut 06340, United States


Fast and accurate assessment of small molecule dihedral energetics is crucial for molecular design and optimization in medicinal chemistry. Yet, accurate prediction of torsion energy profiles remains a challenging task as current molecular mechanics methods are limited by insufficient coverage of druglike chemical space and accurate quantum mechanical (QM) methods are too expensive. To address this limitation, we introduce TorsionNet, a deep neural network (DNN) model specifically developed to predict small molecule torsion energy profiles with QM-level accuracy. We applied active learning to identify nearly 50k fragments (with elements H, C, N, O, F, S, and Cl) that maximized the coverage of our corporate library and leveraged massively parallel cloud computing resources to perform DFT torsion scan of these fragments, generating a training dataset of 1.2 million DFT energies. By training TorsionNet on this dataset, we obtain a model that can rapidly predict the torsion energy profile of typical druglike fragments with DFT-level accuracy. Importantly, our method also provides a direct estimate of the uncertainty in the predicted profiles without any additional calculations. In this report, we show that TorsionNet can reliably identify the preferred dihedral geometries observed in crystal structures. We also present practical applications of TorsionNet that demonstrate how consideration of DNN-based strain energy leads to substantial improvement in existing lead discovery and design workflows. A benchmark dataset (TorsionNet500) comprising 500 chemically diverse fragments with DFT torsion profiles (12k DFT-optimized geometries and energies) has been created and is made freely available.


Email Address of Submitting Author


Pfizer R&D



ORCID For Submitting Author


Declaration of Conflict of Interest

Authors declare no conflict of interest

Version Notes

initial submission