Bayesian optimization for conformer generation

Journal of Cheminformatics

Geoffrey R. Hutchison¹  Garrett M. Morris²  Lucian Chan² 

[1] Department of Chemistry and Chemical Engineering, University of Pittsburgh;Department of Statistics, University of Oxford;

关键词: Bayesian optimization;    Gaussian processes;    Conformer generation;    Rotatable bond;    Torsion angle;    Conformational space;   

DOI  :  10.1186/s13321-019-0354-7

来源: DOAJ

【 摘 要 】

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.

【 授权许可】

Unknown