We compared four conformational search algorithms: one systematic method as implemented by Confab in Open Babel ("Confab"); and three stochastic search methods: uniform random search ("Uniform"), plus two variants of Bayesian optimization algorithm ("BOA") each with different acquisition functions: BOA with expected improvement ("EI"), and BOA with Gaussian Process lower confidence bound ("LCB"). We used the Python package, GPyOpt, for the Bayesian optimization algorithm variants and numpy to generate random numbers between 0 and 2π for the uniform search. Pybel is used to drive the torsion angles of the molecules for both uniform random search and BOA. This means that the molecule's bond lengths, bond angles, and ring systems remain unchanged throughout these searches. Moreover, it is possible to sample torsion angles that generate steric clashes in the stochastic searches, which will return high energies.