AI is a viable alternative to high throughput screening: a 318-target study

High throughput screening (HTS) is routinely used to identify bioactive small molecules. This requires physical compounds, which limits coverage of accessible chemical space. Computational approaches combined with vast on-demand chemical libraries can access far greater chemical space, provided that the predictive accuracy is sufficient to identify useful molecules. Through the largest and most diverse virtual HTS campaign reported to date, comprising 318 individual projects, we demonstrate that our AtomNet® convolutional neural network successfully finds novel hits across every major therapeutic area and protein class. We address historical limitations of computational screening by demonstrating success for target proteins without known binders, high-quality X-ray crystal structures, or manual cherry-picking of compounds. We show that the molecules selected by the AtomNet® model are novel drug-like scaffolds rather than minor modifications to known bioactive compounds. Our empirical results suggest that computational methods can substantially replace HTS as the first step of small-molecule drug discovery.

Discovery of Novel Trace Amine-Associated Receptor 5 (TAAR5) Antagonists Using a Deep Convolutional Neural Network

Trace amine-associated receptor 5 (TAAR5) is a G protein-coupled receptor that belongs to the TAARs family (TAAR1-TAAR9). TAAR5 is expressed in the olfactory epithelium and is responsible for sensing 3-methylamine (TMA). However, recent studies showed that TAAR5 is also expressed in the limbic brain regions and is involved in the regulation of emotional behaviour and adult neurogenesis, suggesting that TAAR5 antagonism may represent a novel therapeutic strategy for anxiety and depression. We used the AtomNet^® model, the first deep learning neural network for structure-based drug discovery, to identify putative TAAR5 ligands and tested them in an in vitro BRET assay. We found two mTAAR5 antagonists with low to submicromolar activity that are able to inhibit the cAMP production induced by TMA. Moreover, these two compounds also inhibited the mTAAR5 downstream signalling, such as the phosphorylation of CREB and ERK. These two hits exhibit drug-like properties and could be used to further develop more potent TAAR5 ligands with putative anxiolytic and antidepressant activity.

Predicting Productive Binding Modes for Substrates and Carbocation Intermediates in Terpene Synthases-Bornyl Diphosphate Synthase as a Representative Case

Terpene synthases comprise a family of enzymes that convert acyclic oligo-isoprenyl diphosphates to terpene natural products with complex, polycyclic carbon backbones via the generation and protection of carbocation intermediates. To accommodate this chemistry, terpene synthase active sites generally are lined with alkyl and aromatic, i.e., nonpolar, sidechains. Predicting the correct, mechanistically relevant binding modes for entire terpene synthase reaction pathways remains an unsolved challenge. Here we describe a method for identifying such modes: TerDockin, a series of protocols to predict the orientation of carbon skeletons of substrates and derived carbocations relative to the bound diphosphate group in terpene synthase active sites. Using this recipe for bornyl diphosphate synthase, we have predicted binding modes that are consistent with all current experimental observations, including the results of isotope labeling experiments and known stereoselectivity. In addition, the predicted binding modes recapitulate key findings of a seminal study involving more computationally demanding QM/MM molecular dynamics methods on part of this pathway. This work illustrates the value of the TerDockin approach as a starting point for more involved calculations and sets the stage for the rational engineering of this family of enzymes.

Changing Face: A Key Residue for the Addition of Water by Sclareol Synthase

Sclareol synthase from Salvia sclarea (SsSS) naturally acts on 8α-hydroxy-copalyl diphosphate (1), stereoselectively adding water to produce (13R)-sclareol (2a), and similarly yields hydroxylated products with manifold other such bicyclic diterpene precursors. Here a key residue for this addition of water was identified. Strikingly, substitution with glutamine switches stereochemical outcome with 1, leading to selective production of (13S)-sclareol (2b). Moreover, changes to the stereospecificity of water addition with the structurally closely-related substrate copalyl diphosphate (4) could be accomplished with alternative substitutions. Thus, this approach is expected to provide biosynthetic access to both epimers of 13-hydroxylated derivatives of manifold labdane-related diterpenes.

Elucidating Substrate Promiscuity within the FabI Enzyme Family

The rapidly growing appreciation of enzymes' catalytic and substrate promiscuity may lead to their expanded use in the fields of chemical synthesis and industrial biotechnology. Here, we explore the substrate promiscuity of enoyl-acyl carrier protein reductases (commonly known as FabI) and how that promiscuity is a function of inherent reactivity and the geometric demands of the enzyme's active site. We demonstrate that these enzymes catalyze the reduction of a wide range of substrates, particularly α,β-unsaturated aldehydes. In addition, we demonstrate that a combination of quantum mechanical hydride affinity calculations and molecular docking can be used to rapidly categorize compounds that FabI can use as substrates. The results here provide new insight into the determinants of catalysis for FabI and set the stage for the development of a new assay for drug discovery, organic synthesis, and novel biocatalysts.

Systematic Functional Analysis of Active-Site Residues in l-Threonine Dehydrogenase from Thermoplasma volcanium

Enzymes have been through millions of years of evolution during which their active-site microenvironments are fine-tuned. Active-site residues are commonly conserved within protein families, indicating their importance for substrate recognition and catalysis. In this work, we systematically mutated active-site residues of l-threonine dehydrogenase from Thermoplasma volcanium and characterized the mutants against a panel of substrate analogs. Our results demonstrate that only a subset of these residues plays an essential role in substrate recognition and catalysis and that the native enzyme activity can be further enhanced roughly 4.6-fold by a single point mutation. Kinetic characterization of mutants on substrate analogs shows that l-threonine dehydrogenase possesses promiscuous activities toward other chemically similar compounds not previously observed. Quantum chemical calculations on the hydride-donating ability of these substrates also reveal that this enzyme did not evolve to harness the intrinsic substrate reactivity for enzyme catalysis. Our analysis provides insights into connections between the details of enzyme active-site structure and specific function. These results are directly applicable to rational enzyme design and engineering.

Synthesis and Evaluation of Quinazolines as Inhibitors of the Bacterial Cell Division Protein FtsZ

The bacterial cell division protein FtsZ is one of many potential targets for the development of novel antibiotics. Recently, zantrin Z3 was shown to be a cross-species inhibitor of FtsZ; however, its specific interactions with the protein are still unknown. Herein we report the synthesis of analogues that contain a more tractable core structure and an analogue with single-digit micromolar inhibition of FtsZ's GTPase activity, which represents the most potent inhibitor of Escherichia coli FtsZ reported to date. In addition, the zantrin Z3 core has been converted to two potential photo-cross-linking reagents for proteomic studies that could shed light on the molecular interactions between FtsZ and molecules related to zantrin Z3.

Comparison of small molecule inhibitors of the bacterial cell division protein FtsZ and identification of a reliable cross-species inhibitor

FtsZ is a guanosine triphosphatase (GTPase) that mediates cytokinesis in bacteria. FtsZ is homologous in structure to eukaryotic tubulin and polymerizes in a similar head-to-tail fashion. The study of tubulin's function in eukaryotic cells has benefited greatly from specific and potent small molecule inhibitors, including colchicine and taxol. Although many small molecule inhibitors of FtsZ have been reported, none has emerged as a generally useful probe for modulating bacterial cell division. With the goal of establishing a useful and reliable small molecule inhibitor of FtsZ, a broad biochemical cross-comparison of reported FtsZ inhibitors was undertaken. Several of these molecules, including phenolic natural products, are unselective inhibitors that seem to derive their activity from the formation of microscopic colloids or aggregates. Other compounds, including the natural product viriditoxin and the drug development candidate PC190723, exhibit no inhibition of GTPase activity using protocols in this work or under published conditions. Of the compounds studied, only zantrin Z3 exhibits good levels of inhibition, maintains activity under conditions that disrupt small molecule aggregates, and provides a platform for exploration of structure-activity relationships (SAR). Preliminary SAR studies have identified slight modifications to the two side chains of this structure that modulate the inhibitory activity of zantrin Z3. Collectively, these studies will help focus future investigations toward the establishment of probes for FtsZ that fill the roles of colchicine and taxol in studies of tubulin.