BIreactive: A Machine-Learning Model to Estimate Covalent Warhead Reactivity

In the past decade, the pharmaceutical industry has paid closer attention to covalent drugs. Differently from standard noncovalent drugs, these compounds can exhibit peculiar properties, such as higher potency or longer duration of target inhibition with a potentially lower dosage. These properties are mainly driven by the reactive functional group present in the compound, the so-called warhead that forms a covalent bond with a specific nucleophilic amino-acid on the target. In this work, we report the possibility to combine ab initio activation energies with machine-learning to estimate covalent compound intrinsic reactivity. The idea behind this approach is to have a precise estimation of the transition state barriers, and thus of the compound reactivity, but with the speed of a machine-learning algorithm. We call this method "BIreactive". Here, we demonstrate this approach on acrylamides and 2-chloroacetamides, two warhead classes that possess different reaction mechanisms. In combination with our recently implemented truncation algorithm, we also demonstrate the possibility to use BIreactive not only for fragments but also for lead-like molecules. The generic nature of this approach allows also the extension to several other warheads. The combination of these factors makes BIreactive a valuable tool for the covalent drug discovery process in a pharmaceutical context.