AutoDesigner, a De Novo Design Algorithm for Rapidly Exploring Large Chemical Space for Lead Optimization: Application to the Design and Synthesis of d-Amino Acid Oxidase Inhibitors

Comprehensive applications of the artificial intelligence technology in new drug research and development

Purpose

Target-based strategy is a prevalent means of drug research and development (R&D), since targets provide effector molecules of drug action and offer the foundation of pharmacological investigation. Recently, the artificial intelligence (AI) technology has been utilized in various stages of drug R&D, where AI-assisted experimental methods show higher efficiency than sole experimental ones. It is a critical need to give a comprehensive review of AI applications in drug R &D for biopharmaceutical field.

Methods

Relevant literatures about AI-assisted drug R&D were collected from the public databases (Including Google Scholar, Web of Science, PubMed, IEEE Xplore Digital Library, Springer, and ScienceDirect) through a keyword searching strategy with the following terms [("Artificial Intelligence" OR "Knowledge Graph" OR "Machine Learning") AND ("Drug Target Identification" OR "New Drug Development")].

Results

In this review, we first introduced common strategies and novel trends of drug R&D, followed by characteristic description of AI algorithms widely used in drug R&D. Subsequently, we depicted detailed applications of AI algorithms in target identification, lead compound identification and optimization, drug repurposing, and drug analytical platform construction. Finally, we discussed the challenges and prospects of AI-assisted methods for drug discovery.

Conclusion

Collectively, this review provides comprehensive overview of AI applications in drug R&D and presents future perspectives for biopharmaceutical field, which may promote the development of drug industry.

Discovery of novel potent drugs for influenza by inhibiting the vital function of neuraminidase via fragment-based drug design (FBDD) and molecular dynamics simulation strategies

The current work describes a fragment linking methodology to generate new neuraminidase inhibitors. A total number of 28,977 fragments from Zinc 20 have been obtained and screened for neuraminidase receptor affinity. Using Schrödinger software, the highest-scoring 270 fragment hits (with scores greater than -7.6) were subjected to fragment combining to create 100 new molecules. These 100 novel compounds were studied using XP docking to evaluate the molecular interaction modes and their binding affinity to neuraminidase receptor. The top ten molecules were selected, for ADMET, drug-likeness features. Based on these characteristics, the best four developed molecules and Zanamivir were submitted to a molecular dynamics simulation investigation to estimate their dynamics within the neuraminidase receptor using Gromacs software. All MD simulation findings show that the generated complexes are very stable when compared to the clinical inhibitor (Zanamivir). In addition, the four designed neuraminidase inhibitors formed very stable complexes with neuraminidase receptor (with total binding energies ranging from -83.50 to -107.85 Kj/mol) according to the total binding energy calculated by MM-PBSA. For the objective of developing new influenza medications, these novel molecules have the potential to be further evaluated in vitro and in vivo for influenza drug discovery.Communicated by Ramaswamy H. Sarma.

Alchemical Transformations and Beyond: Recent Advances and Real-World Applications of Free Energy Calculations in Drug Discovery

Computational methods constitute efficient strategies for screening and optimizing potential drug molecules. A critical factor in this process is the binding affinity between candidate molecules and targets, quantified as binding free energy. Among various estimation methods, alchemical transformation methods stand out for their theoretical rigor. Despite challenges in force field accuracy and sampling efficiency, advancements in algorithms, software, and hardware have increased the application of free energy perturbation (FEP) calculations in the pharmaceutical industry. Here, we review the practical applications of FEP in drug discovery projects since 2018, covering both ligand-centric and residue-centric transformations. We show that relative binding free energy calculations have steadily achieved chemical accuracy in real-world applications. In addition, we discuss alternative physics-based simulation methods and the incorporation of deep learning into free energy calculations.

AutoDesigner - Core Design, a De Novo Design Algorithm for Chemical Scaffolds: Application to the Design and Synthesis of Novel Selective Wee1 Inhibitors

The hit identification stage of a drug discovery program generally involves the design of novel chemical scaffolds with desired biological activity against the target(s) of interest. One common approach is scaffold hopping, which is the manual design of novel scaffolds based on known chemical matter. One major limitation of this approach is narrow chemical space exploration, which can lead to difficulties in maintaining or improving biological activity, selectivity, and favorable property space. Another limitation is the lack of preliminary structure-activity relationship (SAR) data around these designs, which could lead to selecting suboptimal scaffolds to advance lead optimization. To address these limitations, we propose AutoDesigner - Core Design (CoreDesign), a de novo scaffold design algorithm. Our approach is a cloud-integrated, de novo design algorithm for systematically exploring and refining chemical scaffolds against biological targets of interest. The algorithm designs, evaluates, and optimizes a vast range, from millions to billions, of molecules in silico, following defined project parameters encompassing structural novelty, physicochemical attributes, potency, and selectivity using active-learning FEP. To validate CoreDesign in a real-world drug discovery setting, we applied it to the design of novel, potent Wee1 inhibitors with improved selectivity over PLK1. Starting from a single known ligand and receptor structure, CoreDesign rapidly explored over 23 billion molecules to identify 1,342 novel chemical series with a mean of 4 compounds per scaffold. To rapidly analyze this large amount of data and prioritize chemical scaffolds for synthesis, we utilize t-Distributed Stochastic Neighbor Embedding (t-SNE) plots of in silico properties. The chemical space projections allowed us to rapidly identify a structurally novel 5-5 fused core meeting all the hit-identification requirements. Several compounds were synthesized and assayed from the scaffold, displaying good potency against Wee1 and excellent PLK1 selectivity. Our results suggest that CoreDesign can significantly speed up the hit-identification process and increase the probability of success of drug discovery campaigns by allowing teams to bring forward high-quality chemical scaffolds derisked by the availability of preliminary SAR.

Structures, dynamics, complexes, and functions: From classic computation to artificial intelligence

Computational approaches can provide highly detailed insight into the molecular recognition processes that underlie drug binding, the assembly of protein complexes, and the regulation of biological functional processes. Classical simulation methods can bridge a wide range of length- and time-scales typically involved in such processes. Lately, automated learning and artificial intelligence methods have shown the potential to expand the reach of physics-based approaches, ushering in the possibility to model and even design complex protein architectures. The synergy between atomistic simulations and AI methods is an emerging frontier with a huge potential for advances in structural biology. Herein, we explore various examples and frameworks for these approaches, providing select instances and applications that illustrate their impact on fundamental biomolecular problems.

Ubiquitin-specific protease 7 inhibitors reveal a differentiated mechanism of p53-driven anti-cancer activity

The USP7 deubiquitinase regulates proteins involved in the cell cycle, DNA repair, and epigenetics and has been implicated in cancer progression. USP7 inhibition has been pursued for the development of anti-cancer therapies. Here, we describe the discovery of potent and specific USP7 inhibitors exemplified by FX1-5303. FX1-5303 was used as a chemical probe to study the USP7-mediated regulation of p53 signaling in cells. It demonstrates mechanistic differences compared to MDM2 antagonists, a related class of anti-tumor agents that act along the same pathway. FX1-5303 synergizes with the clinically approved BCL2 inhibitor venetoclax in acute myeloid leukemia (AML) cell lines and ex vivo patient samples and leads to strong tumor growth inhibition in in vivo mouse xenograft models of multiple myeloma and AML. This work introduces new USP7 inhibitors, differentiates their mechanism of action from MDM2 inhibition, and identifies specific opportunities for their use in the treatment of AML.

Exploiting high-energy hydration sites for the discovery of potent peptide aldehyde inhibitors of the SARS-CoV-2 main protease with cellular antiviral activity

Small-molecule antivirals that prevent the replication of the SARS-CoV-2 virus by blocking the enzymatic activity of its main protease (Mpro) are and will be a tenet of pandemic preparedness. However, the peptidic nature of such compounds often precludes the design of compounds within favorable physical property ranges, limiting cellular activity. Here we describe the discovery of peptide aldehyde Mpro inhibitors with potent enzymatic and cellular antiviral activity. This structure-activity relationship (SAR) exploration was guided by the use of calculated hydration site thermodynamic maps (WaterMap) to drive potency via displacement of waters from high-energy sites. Thousands of diverse compounds were designed to target these high-energy hydration sites and then prioritized for synthesis by physics- and structure-based Free-Energy Perturbation (FEP+) simulations, which accurately predicted biochemical potencies. This approach ultimately led to the rapid discovery of lead compounds with unique SAR that exhibited potent enzymatic and cellular activity with excellent pan-coronavirus coverage.

Geometric Deep Learning for Structure-Based Ligand Design

A pervasive challenge in drug design is determining how to expand a ligand-a small molecule that binds to a target biomolecule-in order to improve various properties of the ligand. Adding single chemical groups, known as fragments, is important for lead optimization tasks, and adding multiple fragments is critical for fragment-based drug design. We have developed a comprehensive framework that uses machine learning and three-dimensional protein-ligand structures to address this challenge. Our method, FRAME, iteratively determines where on a ligand to add fragments, selects fragments to add, and predicts the geometry of the added fragments. On a comprehensive benchmark, FRAME consistently improves predicted affinity and selectivity relative to the initial ligand, while generating molecules with more drug-like chemical properties than docking-based methods currently in widespread use. FRAME learns to accurately describe molecular interactions despite being given no prior information on such interactions. The resulting framework for quality molecular hypothesis generation can be easily incorporated into the workflows of medicinal chemists for diverse tasks, including lead optimization, fragment-based drug discovery, and de novo drug design.

FEP Protocol Builder: Optimization of Free Energy Perturbation Protocols Using Active Learning

Significant improvements have been made in the past decade to methods that rapidly and accurately predict binding affinity through free energy perturbation (FEP) calculations. This has been driven by recent advances in small-molecule force fields and sampling algorithms combined with the availability of low-cost parallel computing. Predictive accuracies of ∼1 kcal mol-1 have been regularly achieved, which are sufficient to drive potency optimization in modern drug discovery campaigns. Despite the robustness of these FEP approaches across multiple target classes, there are invariably target systems that do not display expected performance with default FEP settings. Traditionally, these systems required labor-intensive manual protocol development to arrive at parameter settings that produce a predictive FEP model. Due to the (a) relatively large parameter space to be explored, (b) significant compute requirements, and (c) limited understanding of how combinations of parameters can affect FEP performance, manual FEP protocol optimization can take weeks to months to complete, and often does not involve rigorous train-test set splits, resulting in potential overfitting. These manual FEP protocol development timelines do not coincide with tight drug discovery project timelines, essentially preventing the use of FEP calculations for these target systems. Here, we describe an automated workflow termed FEP Protocol Builder (FEP-PB) to rapidly generate accurate FEP protocols for systems that do not perform well with default settings. FEP-PB uses an active-learning workflow to iteratively search the protocol parameter space to develop accurate FEP protocols. To validate this approach, we applied it to pharmaceutically relevant systems where default FEP settings could not produce predictive models. We demonstrate that FEP-PB can rapidly generate accurate FEP protocols for the previously challenging MCL1 system with limited human intervention. We also apply FEP-PB in a real-world drug discovery setting to generate an accurate FEP protocol for the p97 system. FEP-PB is able to generate a more accurate protocol than the expert user, rapidly validating p97 as amenable to free energy calculations. Additionally, through the active-learning workflow, we are able to gain insight into which parameters are most important for a given system. These results suggest that FEP-PB is a robust tool that can aid in rapidly developing accurate FEP protocols and increasing the number of targets that are amenable to the technology.

The Hitchhiker's Guide to Deep Learning Driven Generative Chemistry

This microperspective covers the most recent research outcomes of artificial intelligence (AI) generated molecular structures from the point of view of the medicinal chemist. The main focus is on studies that include synthesis and experimental in vitro validation in biochemical assays of the generated molecular structures, where we analyze the reported structures' relevance in modern medicinal chemistry and their novelty. The authors believe that this review would be appreciated by medicinal chemistry and AI-driven drug design (AIDD) communities and can be adopted as a comprehensive approach for qualifying different research outcomes in AIDD.

A Computational Physics-based Approach to Predict Unbound Brain-to-Plasma Partition Coefficient, Kp,uu

The blood-brain barrier (BBB) plays a critical role in preventing harmful endogenous and exogenous substances from penetrating the brain. Optimal brain penetration of small-molecule central nervous system (CNS) drugs is characterized by a high unbound brain/plasma ratio (K_p,uu). While various medicinal chemistry strategies and in silico models have been reported to improve BBB penetration, they have limited application in predicting K_p,uu directly. We describe a physics-based computational approach, a quantum mechanics (QM)-based energy of solvation (E-sol), to predict K_p,uu. Prospective application of this method in internal CNS drug discovery programs highlights the utility and accuracy of this new method, which showed a categorical accuracy of 79% and an R² of 0.61 from a linear regression model.

Accelerated Discovery of Macrocyclic CDK2 Inhibitor QR-6401 by Generative Models and Structure-Based Drug Design

Selective CDK2 inhibitors have the potential to provide effective therapeutics for CDK2-dependent cancers and for combating drug resistance due to high cyclin E1 (CCNE1) expression intrinsically or CCNE1 amplification induced by treatment of CDK4/6 inhibitors. Generative models that take advantage of deep learning are being increasingly integrated into early drug discovery for hit identification and lead optimization. Here we report the discovery of a highly potent and selective macrocyclic CDK2 inhibitor QR-6401 (23) accelerated by the application of generative models and structure-based drug design (SBDD). QR-6401 (23) demonstrated robust antitumor efficacy in an OVCAR3 ovarian cancer xenograft model via oral administration.

Defining Levels of Automated Chemical Design

One application area of computational methods in drug discovery is the automated design of small molecules. Despite the large number of publications describing methods and their application in both retrospective and prospective studies, there is a lack of agreement on terminology and key attributes to distinguish these various systems. We introduce Automated Chemical Design (ACD) Levels to clearly define the level of autonomy along the axes of ideation and decision making. To fully illustrate this framework, we provide literature exemplars and place some notable methods and applications into the levels. The ACD framework provides a common language for describing automated small molecule design systems and enables medicinal chemists to better understand and evaluate such systems.

Discovery of a Novel Class of d-Amino Acid Oxidase Inhibitors Using the Schrödinger Computational Platform

d-Serine is a coagonist of the N-methyl d-aspartate (NMDA) receptor, a key excitatory neurotransmitter receptor. In the brain, d-serine is synthesized from its l-isomer by serine racemase and is metabolized by the D-amino acid oxidase (DAO, DAAO). Many studies have linked decreased d-serine concentration and/or increased DAO expression and enzyme activity to NMDA dysfunction and schizophrenia. Thus, it is feasible to employ DAO inhibitors for the treatment of schizophrenia and other indications. Powered by the Schrödinger computational modeling platform, we initiated a research program to identify novel DAO inhibitors with the best-in-class properties. The program execution leveraged an hDAO FEP+ model to prospectively predict compound potency. A new class of DAO inhibitors with desirable properties has been discovered from this endeavor. Our modeling technology on this program has not only enhanced the efficiency of structure-activity relationship development but also helped to identify a previously unexplored subpocket for further optimization.