Fast, efficient fragment-based coordinate generation for Open Babel

Abstract

Rapidly predicting an accurate three dimensional geometry of a molecule is a crucial task for cheminformatics and across a wide range of molecular modeling. Consequently, developing a fast, accurate, and open implementation of structure prediction is necessary for reproducible cheminformatics research. We introduce a fragment-based coordinate generation implementation for Open Babel, a widely-used open source toolkit for cheminformatics. The new implementation improves speed and stereochemical accuracy, while retaining or improving accuracy of bond lengths, bond angles, and dihedral torsions. Input molecules are broken into fragments by cutting at rotatable bonds. The coordinates of fragments are set according to a fragment library, prepared from open crystallographic databases. Since the coordinates of multiple atoms are decided at once, coordinate prediction is accelerated over the previous rules-based implementation in Open Babel, as well as the widely-used distance geometry methods in RDKit. This new implementation will be beneficial for a wide range of applications, including computational property prediction in polymers, molecular materials and drug design. Electronic supplementary material The online version of this article (10.1186/s13321-019-0372-5) contains supplementary material, which is available to authorized users.

EvoMol: a flexible and interpretable evolutionary algorithm for unbiased de novo molecular generation

“…The B3LYP hybrid functional and the small basis set 3-21G are chosen to reduce the computational cost [44]. We only perform DFT optimisation on molecules conserving the same SMILES after molecular mechanics 3-D coordinates generation using Open Babel [45]. Likewise, we consider a DFT result as valid only if the SMILES remained identical after the geometric optimisation.…”

Machine learning based energy-free structure predictions of molecules, transition states, and solids

“…State of the art approaches for generating 3D molecular structures e.g., ETKDG 6 and Gen3D 7 are very efficient yet carry significant bias since they are based on mathematically rigid functional forms, empirical parameters, knowledge-based heuristic rules, and do not directly improve upon the increase of training data set sizes. While applicable to known and well-behaved regions of chemical compound space, these methods lack generality and are inherently limited when it comes to more challenging systems, such as carbene molecules or transition states (TS).…”

Machine Learning to Accelerate Screening for Marcus Reorganization Energies

“…Input files for each oligomer were created by combining the corresponding SMILES strings of its monomers (Table S1) and using OpenBabel version 3.1.0 26 to generate a 3D geometry. 27 All GFN2 calculations were performed using xTB version 6.0 15 . All DFT calculations were done using the B3LYP functional 17,18 with the 6-31G* basis set, 28 calculated with Gaussian 09, 29 for comparison with previously published internal reorganiztion energies.…”