Alchemical Space Exploration in Drug Design via Stochastic Expanding Boundary Optimization Search

Abstract

In this work, we propose a stochastic algorithm for the emergence of organic molecules with specific properties by steering alchemical changes, allowing the sampling of chemical space for potential drug candidates. The algorithm employs a tree-like local search capable of overcoming entrapment in local "basins" and iteratively sampling adjacent basins, creating a molecular ensemble with specific properties. The method consists of iterative stochastic "moves", which may alter the chemical constitution of a potential drug leading towards molecules with desired properties. Implementation of such a procedure necessitates the Simplified Molecular Input Line Entry System (SMILES) molecular representation. The method has been designed especially for searches in target discontinuous spaces, where local optimal solutions are expected well separated and hardly traverse. We demonstrate our method in three applications: (i) generating ensembles of organic molecules with preset octanol-water partition coefficients, (ii) identifying molecules with specific light absorption properties via quantum mechanical calculations, and (iii) optimizing drug binding free energies to protein targets using molecular docking.