Bayesian optimization for conformer generation

Abstract Generating low-energy molecular conformers is a key task for many areas of computational chemistry, molecular modeling and cheminformatics. Most current conformer generation methods primarily focus on generating geometrically diverse conformers rather than finding the most probable or energ...

Full description

Bibliographic Details
Main Authors: Lucian Chan, Geoffrey R. Hutchison, Garrett M. Morris
Format: Article
Language:English
Published: BMC 2019-05-01
Series:Journal of Cheminformatics
Subjects:
Online Access:http://link.springer.com/article/10.1186/s13321-019-0354-7
_version_ 1811331852024152064
author Lucian Chan
Geoffrey R. Hutchison
Garrett M. Morris
author_facet Lucian Chan
Geoffrey R. Hutchison
Garrett M. Morris
author_sort Lucian Chan
collection DOAJ
description Abstract Generating low-energy molecular conformers is a key task for many areas of computational chemistry, molecular modeling and cheminformatics. Most current conformer generation methods primarily focus on generating geometrically diverse conformers rather than finding the most probable or energetically lowest minima. Here, we present a new stochastic search method called the Bayesian optimization algorithm (BOA) for finding the lowest energy conformation of a given molecule. We compare BOA with uniform random search, and systematic search as implemented in Confab, to determine which method finds the lowest energy. Energetic difference, root-mean-square deviation, and torsion fingerprint deviation are used to quantify the performance of the conformer search algorithms. In general, we find BOA requires far fewer evaluations than systematic or uniform random search to find low-energy minima. For molecules with four or more rotatable bonds, Confab typically evaluates $$10^{4}$$ 104  (median) conformers in its search, while BOA only requires $$10^{2}$$ 102 energy evaluations to find top candidates. Despite using evaluating fewer conformers, 20–40% of the time BOA finds lower-energy conformations than a systematic Confab search for molecules with four or more rotatable bonds.
first_indexed 2024-04-13T16:27:04Z
format Article
id doaj.art-bdac2ad431034711afcb83fea6fe8d7d
institution Directory Open Access Journal
issn 1758-2946
language English
last_indexed 2024-04-13T16:27:04Z
publishDate 2019-05-01
publisher BMC
record_format Article
series Journal of Cheminformatics
spelling doaj.art-bdac2ad431034711afcb83fea6fe8d7d2022-12-22T02:39:43ZengBMCJournal of Cheminformatics1758-29462019-05-0111111110.1186/s13321-019-0354-7Bayesian optimization for conformer generationLucian Chan0Geoffrey R. Hutchison1Garrett M. Morris2Department of Statistics, University of OxfordDepartment of Chemistry and Chemical Engineering, University of PittsburghDepartment of Statistics, University of OxfordAbstract Generating low-energy molecular conformers is a key task for many areas of computational chemistry, molecular modeling and cheminformatics. Most current conformer generation methods primarily focus on generating geometrically diverse conformers rather than finding the most probable or energetically lowest minima. Here, we present a new stochastic search method called the Bayesian optimization algorithm (BOA) for finding the lowest energy conformation of a given molecule. We compare BOA with uniform random search, and systematic search as implemented in Confab, to determine which method finds the lowest energy. Energetic difference, root-mean-square deviation, and torsion fingerprint deviation are used to quantify the performance of the conformer search algorithms. In general, we find BOA requires far fewer evaluations than systematic or uniform random search to find low-energy minima. For molecules with four or more rotatable bonds, Confab typically evaluates $$10^{4}$$ 104  (median) conformers in its search, while BOA only requires $$10^{2}$$ 102 energy evaluations to find top candidates. Despite using evaluating fewer conformers, 20–40% of the time BOA finds lower-energy conformations than a systematic Confab search for molecules with four or more rotatable bonds.http://link.springer.com/article/10.1186/s13321-019-0354-7Bayesian optimizationGaussian processesConformer generationRotatable bondTorsion angleConformational space
spellingShingle Lucian Chan
Geoffrey R. Hutchison
Garrett M. Morris
Bayesian optimization for conformer generation
Journal of Cheminformatics
Bayesian optimization
Gaussian processes
Conformer generation
Rotatable bond
Torsion angle
Conformational space
title Bayesian optimization for conformer generation
title_full Bayesian optimization for conformer generation
title_fullStr Bayesian optimization for conformer generation
title_full_unstemmed Bayesian optimization for conformer generation
title_short Bayesian optimization for conformer generation
title_sort bayesian optimization for conformer generation
topic Bayesian optimization
Gaussian processes
Conformer generation
Rotatable bond
Torsion angle
Conformational space
url http://link.springer.com/article/10.1186/s13321-019-0354-7
work_keys_str_mv AT lucianchan bayesianoptimizationforconformergeneration
AT geoffreyrhutchison bayesianoptimizationforconformergeneration
AT garrettmmorris bayesianoptimizationforconformergeneration