ConfRank: Improving GFN-FF Conformer Ranking with Pairwise Training

16 September 2024, Version 1
This content is a preprint and has not undergone peer review at the time of posting.

Abstract

Conformer ranking is a crucial task for drug discovery, with methods for generating conformers often based on molecular (meta)dynamics. These methods are constrained by the underlying force computation regarding runtime and energy ranking accuracy, limiting their effectiveness for large-scale screening applications. To address these ranking limitations, we introduce ConfRank, a machine learning-based approach that enhances conformer ranking using pairwise training. We demonstrate its performance using GFN-FF-generated conformer ensembles, leveraging the DimeNet++ architecture trained on pairs of 159 760 uncharged organic compounds from the GEOM dataset with r2SCAN-3c reference level. Instead of predicting only on single molecules, this approach captures relative energy differences between conformers, leading to a significant improvement of the overall conformational ranking, outperforming GFN-FF and GFN2-xTB. Thereby, the pairwise RMSD of the relative energy difference of two conformers can be reduced from 5.65 kcal mol−1 to 0.71 kcal mol−1 on the test dataset, allowing to correctly identify up to 81 % of all lowest lying conformers correctly (GFN-FF: 10 %, GFN2- xTB: 47 %). The ConfRank approach is cost-effective, allowing for scalable deployment on both CPU and GPU, achieving runtime accelerations by up to two orders of magnitude compared to GFN2-xTB. Out-of-sample investigations on CREST-generated conformer ensembles from the QM9 dataset and conformers taken from an extended GMTKN55 dataset show promising results for the robustness of this approach. Thereby, ranking correlation coefficient such as Spearman can be improved to 0.90 (GFN-FF: 0.39, GFN2-xTB: 0.84) reducing the probability of an incorrect sign flip in pairwise energy comparison from 32 % to 7 %. On the extended GMTKN55 subsets the pairwise MAD (RMSD) could be reduced on almost all subsets by up to 62 % (58 %) with an average improvement of 30 % (29 %). Furthermore, we motivate the usage of the pairwise training approach from a theoretical perspective, highlighting that while pairwise training can lead to a decline in single sample prediction of absolute energies for ML models, it significantly enhances conformer ranking performance. The data and models used in this study are available at https://github.com/grimme-lab/confrank.

Keywords

conformer ranking
pairwise training
machine learning
GFN-FF

Supplementary materials

Title
Description
Actions
Title
Supporting Information
Description
The Supporting Information contains additional material such as explanations on the metrics used in the manuscript, information regarding the datasets used and tabular overviews of model performances.
Actions

Supplementary weblinks

Comments

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.