Ultra-large library docking for discovering new chemotypes

Molecular determinants of antagonist interactions with chemokine receptors CCR2 and CCR5

By driving monocyte chemotaxis, the chemokine receptor CCR2 shapes inflammatory responses and the formation of tumor microenvironments. This makes it a promising target in inflammation and immuno-oncology; however, despite extensive efforts, there are no FDA-approved CCR2-targeting therapeutics. Cited challenges include the redundancy of the chemokine system, suboptimal properties of compound candidates, and species differences that confound the translation of results from animals to humans. Structure-based drug design can rationalize and accelerate the discovery and optimization of CCR2 antagonists to address these challenges. The prerequisites for such efforts include an atomic-level understanding of the molecular determinants of action of existing antagonists. In this study, using molecular docking and artificial-intelligence-powered compound library screening, we uncover the structural principles of small molecule antagonism and selectivity towards CCR2 and its sister receptor CCR5. CCR2 orthosteric inhibitors are shown to universally occupy an inactive-state-specific tunnel between receptor helices 1 and 7; we also discover an unexpected role for an extra-helical groove accessible through this tunnel, suggesting its potential as a new targetable interface for CCR2 and CCR5 modulation. By contrast, only shape complementarity and limited helix 8 hydrogen bonding govern the binding of various chemotypes of allosteric antagonists. CCR2 residues S1012.63 and V2446.36 are implicated as determinants of CCR2/CCR5 and human/mouse orthosteric and allosteric antagonist selectivity, respectively, and the role of S1012.63 is corroborated through experimental gain-of-function mutagenesis. We establish a critical role of induced fit in antagonist recognition, reveal strong chemotype selectivity of existing structures, and demonstrate the high predictive potential of a new deep-learning-based compound scoring function. Finally, this study expands the available CCR2 structural landscape with computationally generated chemotype-specific models well-suited for structure-based antagonist design.

DockOpt: A Tool for Automatic Optimization of Docking Models

Molecular docking is a widely used technique for leveraging protein structure for ligand discovery, but it remains difficult to utilize due to limitations that have not been adequately addressed. Despite some progress toward automation, docking still requires expert guidance, hindering its adoption by a broader range of investigators. To make docking more accessible, we developed a new utility called DockOpt, which automates the creation, evaluation, and optimization of docking models prior to their deployment in large-scale prospective screens. DockOpt outperforms our previous automated pipeline across all 43 targets in the DUDE-Z benchmark data set, and the generated models for 84% of targets demonstrate sufficient enrichment to warrant their use in prospective screens, with normalized LogAUC values of at least 15%. DockOpt is available as part of the Python package Pydock3 included in the UCSF DOCK 3.8 distribution, which is available for free to academic researchers at https://dock.compbio.ucsf.edu and free for everyone upon registration at https://tldr.docking.org.

Structure-Based Identification of Organoruthenium Compounds as Nanomolar Antagonists of Cannabinoid Receptors

Metal complexes exhibit a diverse range of coordination geometries, representing novel privileged scaffolds with convenient click types of preparation inaccessible for typical carbon-centered organic compounds. Herein, we explored the opportunity to identify biologically active organometallic complexes by reverse docking of a rigid, minimum-size octahedral organoruthenium scaffold against thousands of protein-binding pockets. Interestingly, cannabinoid receptor type 1 (CB1) was identified based on the docking scores and the degree of overlap between the docked organoruthenium scaffold and the hydrophobic scaffold of the cocrystallized ligand. Further structure-based optimization led to the discovery of organoruthenium complexes with nanomolar binding affinities and high selectivity toward CB2. Our work indicates that octahedral organoruthenium scaffolds may be advantageous for targeting the large and hydrophobic binding pockets and that the reverse docking approach may facilitate the discovery of novel privileged scaffolds, such as organometallic complexes, for exploring chemical space in lead discovery.

TopoFormer: Multiscale Topology-enabled Structure-to-Sequence Transformer for Protein-Ligand Interaction Predictions

Pre-trained deep Transformers have had tremendous success in a wide variety of disciplines. However, in computational biology, essentially all Transformers are built upon the biological sequences, which ignores vital stereochemical information and may result in crucial errors in downstream predictions. On the other hand, three-dimensional (3D) molecular structures are incompatible with the sequential architecture of Transformer and natural language processing (NLP) models in general. This work addresses this foundational challenge by a topological Transformer (TopoFormer). TopoFormer is built by integrating NLP and a multiscale topology techniques, the persistent topological hyperdigraph Laplacian (PTHL), which systematically converts intricate 3D protein-ligand complexes at various spatial scales into a NLP-admissible sequence of topological invariants and homotopic shapes. Element-specific PTHLs are further developed to embed crucial physical, chemical, and biological interactions into topological sequences. TopoFormer surges ahead of conventional algorithms and recent deep learning variants and gives rise to exemplary scoring accuracy and superior performance in ranking, docking, and screening tasks in a number of benchmark datasets. The proposed topological sequences can be extracted from all kinds of structural data in data science to facilitate various NLP models, heralding a new era in AI-driven discovery.

Generative artificial intelligence in drug discovery: basic framework, recent advances, challenges, and opportunities

There are two main ways to discover or design small drug molecules. The first involves fine-tuning existing molecules or commercially successful drugs through quantitative structure-activity relationships and virtual screening. The second approach involves generating new molecules through de novo drug design or inverse quantitative structure-activity relationship. Both methods aim to get a drug molecule with the best pharmacokinetic and pharmacodynamic profiles. However, bringing a new drug to market is an expensive and time-consuming endeavor, with the average cost being estimated at around $2.5 billion. One of the biggest challenges is screening the vast number of potential drug candidates to find one that is both safe and effective. The development of artificial intelligence in recent years has been phenomenal, ushering in a revolution in many fields. The field of pharmaceutical sciences has also significantly benefited from multiple applications of artificial intelligence, especially drug discovery projects. Artificial intelligence models are finding use in molecular property prediction, molecule generation, virtual screening, synthesis planning, repurposing, among others. Lately, generative artificial intelligence has gained popularity across domains for its ability to generate entirely new data, such as images, sentences, audios, videos, novel chemical molecules, etc. Generative artificial intelligence has also delivered promising results in drug discovery and development. This review article delves into the fundamentals and framework of various generative artificial intelligence models in the context of drug discovery via de novo drug design approach. Various basic and advanced models have been discussed, along with their recent applications. The review also explores recent examples and advances in the generative artificial intelligence approach, as well as the challenges and ongoing efforts to fully harness the potential of generative artificial intelligence in generating novel drug molecules in a faster and more affordable manner. Some clinical-level assets generated form generative artificial intelligence have also been discussed in this review to show the ever-increasing application of artificial intelligence in drug discovery through commercial partnerships.

Exploring the combinatorial explosion of amine-acid reaction space via graph editing

Amines and carboxylic acids are abundant chemical feedstocks that are nearly exclusively united via the amide coupling reaction. The disproportionate use of the amide coupling leaves a large section of unexplored reaction space between amines and acids: two of the most common chemical building blocks. Herein we conduct a thorough exploration of amine-acid reaction space via systematic enumeration of reactions involving a simple amine-carboxylic acid pair. This approach to chemical space exploration investigates the coarse and fine modulation of physicochemical properties and molecular shapes. With the invention of reaction methods becoming increasingly automated and bringing conceptual reactions into reality, our map provides an entirely new axis of chemical space exploration for rational property design.

Integrating QSAR modelling and deep learning in drug discovery: the emergence of deep QSAR

Quantitative structure-activity relationship (QSAR) modelling, an approach that was introduced 60 years ago, is widely used in computer-aided drug design. In recent years, progress in artificial intelligence techniques, such as deep learning, the rapid growth of databases of molecules for virtual screening and dramatic improvements in computational power have supported the emergence of a new field of QSAR applications that we term 'deep QSAR'. Marking a decade from the pioneering applications of deep QSAR to tasks involved in small-molecule drug discovery, we herein describe key advances in the field, including deep generative and reinforcement learning approaches in molecular design, deep learning models for synthetic planning and the application of deep QSAR models in structure-based virtual screening. We also reflect on the emergence of quantum computing, which promises to further accelerate deep QSAR applications and the need for open-source and democratized resources to support computer-aided drug design.

AI for targeted polypharmacology: The next frontier in drug discovery

In drug discovery, targeted polypharmacology, i.e., targeting multiple molecular targets with a single drug, is redefining therapeutic design to address complex diseases. Pre-selected pharmacological profiles, as exemplified in kinase drugs, promise enhanced efficacy and reduced toxicity. Historically, many of such drugs were discovered serendipitously, limiting predictability and efficacy, but currently artificial intelligence (AI) offers a transformative solution. Machine learning and deep learning techniques enable modeling protein structures, generating novel compounds, and decoding their polypharmacological effects, opening an avenue for more systematic and predictive multi-target drug design. This review explores the use of AI in identifying synergistic co-targets and delineating them from anti-targets that lead to adverse effects, and then discusses advances in AI-enabled docking, generative chemistry, and proteochemometric modeling of proteome-wide compound interactions, in the context of polypharmacology. We also provide insights into challenges ahead.

Curation and cheminformatics analysis of a Ugi-reaction derived library (URDL) of synthetically tractable small molecules for virtual screening application

Virtual screening (VS) is an important approach in drug discovery and relies on the availability of a virtual library of synthetically tractable molecules. Ugi reaction (UR) represents an important multi-component reaction (MCR) that reliably produces a peptidomimetic scaffold. Recent literature shows that a tactically assembled Ugi adduct can be subjected to further chemical modifications to yield a variety of rings and scaffolds, thus, renewing the interest in this old reaction. Given the reliability and efficiency of UR, we collated an UR derived library (URDL) of small molecules (total = 5773) for VS. The synthesis of the majority of URDL molecules may be carried out in 1-2 pots in a time and cost-effective manner. The detailed analysis of the average property and chemical space of URDL was also carried out using the open-source Datawarrior program. The comparison with FDA-approved oral drugs and inhibitors of protein-protein interactions (iPPIs) suggests URDL molecules are 'clean', drug-like, and conform to a structurally distinct space from the other two categories. The average physicochemical properties of compounds in the URDL library lie closer to iPPI molecules than oral drugs thus suggesting that the URDL resource can be applied to discover novel iPPI molecules. The URDL molecules consist of diverse ring systems, many of which have not been exploited yet for drug design. Thus, URDL represents a small virtual library of drug-like molecules with unexplored chemical space designed for VS. The structures of all molecules of URDL, oral drugs, and iPPI compounds are being made freely accessible as supplementary information for broader application.

Computational drug development for membrane protein targets

The application of computational biology in drug development for membrane protein targets has experienced a boost from recent developments in deep learning-driven structure prediction, increased speed and resolution of structure elucidation, machine learning structure-based design and the evaluation of big data. Recent protein structure predictions based on machine learning tools have delivered surprisingly reliable results for water-soluble and membrane proteins but have limitations for development of drugs that target membrane proteins. Structural transitions of membrane proteins have a central role during transmembrane signaling and are often influenced by therapeutic compounds. Resolving the structural and functional basis of dynamic transmembrane signaling networks, especially within the native membrane or cellular environment, remains a central challenge for drug development. Tackling this challenge will require an interplay between experimental and computational tools, such as super-resolution optical microscopy for quantification of the molecular interactions of cellular signaling networks and their modulation by potential drugs, cryo-electron microscopy for determination of the structural transitions of proteins in native cell membranes and entire cells, and computational tools for data analysis and prediction of the structure and function of cellular signaling networks, as well as generation of promising drug candidates.

Exploring Amino Acid Transporters as Therapeutic Targets for Cancer: An Examination of Inhibitor Structures, Selectivity Issues, and Discovery Approaches

Amino acid transporters are abundant amongst the solute carrier family and have an important role in facilitating the transfer of amino acids across cell membranes. Because of their impact on cell nutrient distribution, they also appear to have an important role in the growth and development of cancer. Naturally, this has made amino acid transporters a novel target of interest for the development of new anticancer drugs. Many attempts have been made to develop inhibitors of amino acid transporters to slow down cancer cell growth, and some have even reached clinical trials. The purpose of this review is to help organize the available information on the efforts to discover amino acid transporter inhibitors by focusing on the amino acid transporters ASCT2 (SLC1A5), LAT1 (SLC7A5), xCT (SLC7A11), SNAT1 (SLC38A1), SNAT2 (SLC38A2), and PAT1 (SLC36A1). We discuss the function of the transporters, their implication in cancer, their known inhibitors, issues regarding selective inhibitors, and the efforts and strategies of discovering inhibitors. The goal is to encourage researchers to continue the search and development within the field of cancer treatment research targeting amino acid transporters.

Potentiation of Antibiotic Activity of Aztreonam against Metallo-β-Lactamase-Producing Multidrug-Resistant Pseudomonas aeruginosa by 3-O-Substituted Difluoroquercetin Derivatives

The combination of aztreonam (ATM) and ceftazidime-avibactam (CAZ-AVI; CZA) has shown therapeutic potential against serine-β-lactamase (SBL)- and metallo-β-lactamase (MBL)-producing Enterobacterales. However, the ability of CZA to restore the antibiotic activity of ATM is severely limited in MBL-producing multidrug-resistant (MDR) Pseudomonas aeruginosa strains because of the myriad of intrinsic and acquired resistance mechanisms associated with this pathogen. We reasoned that the simultaneous inhibition of multiple targets associated with multidrug resistance mechanisms may potentiate the antibiotic activity of ATM against MBL-producing P. aeruginosa. During a search for the multitarget inhibitors through a molecular docking study, we discovered that di-F-Q, the previously reported efflux pump inhibitor of MDR P. aeruginosa, binds to the active sites of the efflux pump (MexB), as well as various β-lactamases, and these sites are open to the 3-O-position of di-F-Q. The 3-O-substituted di-F-Q derivatives were thus synthesized and showed hereto unknown multitarget MDR inhibitory activity against various ATM-hydrolyzing β-lactamases (AmpC, KPC, and New Delhi metallo-β-lactamase (NDM)) and the efflux pump of P. aeruginosa, presumably by forming additional hydrophobic contacts with the targets. The multitarget MDR inhibitor 27 effectively potentiated the antimicrobial activity of ATM and reduced the MIC of ATM more than four-fold in 19 out of 21 MBL-producing P. aeruginosa clinical strains, including the NDM-producing strains which were highly resistant to various combinations of ATM with β-lactamase inhibitors and/or efflux pump inhibitors. Our findings suggest that the simultaneous inhibition of multiple MDR targets might provide new avenues for the discovery of safe and efficient MDR reversal agents which can be used in combination with ATM against MBL-producing MDR P. aeruginosa.

An intriguing vision for transatlantic collaborative health data use and artificial intelligence development

Our traditional approach to diagnosis, prognosis, and treatment, can no longer process and transform the enormous volume of information into therapeutic success, innovative discovery, and health economic performance. Precision health, i.e., the right treatment, for the right person, at the right time in the right place, is enabled through a learning health system, in which medicine and multidisciplinary science, economic viability, diverse culture, and empowered patient's preferences are digitally integrated and conceptually aligned for continuous improvement and maintenance of health, wellbeing, and equity. Artificial intelligence (AI) has been successfully evaluated in risk stratification, accurate diagnosis, and treatment allocation, and to prevent health disparities. There is one caveat though: dependable AI models need to be trained on population-representative, large and deep data sets by multidisciplinary and multinational teams to avoid developer, statistical and social bias. Such applications and models can neither be created nor validated with data at the country, let alone institutional level and require a new dimension of collaboration, a cultural change with the establishment of trust in a precompetitive space. The Data for Health (#DFH23) conference in Berlin and the Follow-Up Workshop at Harvard University in Boston hosted a representative group of stakeholders in society, academia, industry, and government. With the momentum #DFH23 created, the European Health Data Space (EHDS) as a solid and safe foundation for consented collaborative health data use and the G7 Hiroshima AI process in place, we call on citizens and their governments to fully support digital transformation of medicine, research and innovation including AI.

Discovery of Nanomolar Inhibitors for Human Dihydroorotate Dehydrogenase Using Structure-Based Drug Discovery Methods

We used a structure-based drug discovery approach to identify novel inhibitors of human dihydroorotate dehydrogenase (DHODH), which is a therapeutic target for treating cancer and autoimmune and inflammatory diseases. In the case of acute myeloid leukemia, no previously discovered DHODH inhibitors have yet succeeded in this clinical application. Thus, there remains a strong need for new inhibitors that could be used as alternatives to the current standard-of-care. Our goal was to identify novel inhibitors of DHODH. We implemented prefiltering steps to omit PAINS and Lipinski violators at the earliest stages of this project. This enriched compounds in the data set that had a higher potential of favorable oral druggability. Guided by Glide SP docking scores, we found 20 structurally unique compounds from the ChemBridge EXPRESS-pick library that inhibited DHODH with IC50, DHODH values between 91 nM and 2.7 μM. Ten of these compounds reduced MOLM-13 cell viability with IC50, MOLM-13 values between 2.3 and 50.6 μM. Compound 16 (IC50, DHODH = 91 nM) inhibited DHODH more potently than the known DHODH inhibitor, teriflunomide (IC50, DHODH = 130 nM), during biochemical characterizations and presented a promising scaffold for future hit-to-lead optimization efforts. Compound 17 (IC50, MOLM-13 = 2.3 μM) was most successful at reducing survival in MOLM-13 cell lines compared with our other hits. The discovered compounds represent excellent starting points for the development and optimization of novel DHODH inhibitors.

Docking for Molecules That Bind in a Symmetric Stack with SymDOCK

Discovering ligands for amyloid fibrils, such as those formed by the tau protein, is an area of great current interest. In recent structures, ligands bind in stacks in the tau fibrils to reflect the rotational and translational symmetry of the fibril itself; in these structures, the ligands make few interactions with the protein but interact extensively with each other. To exploit this symmetry and stacking, we developed SymDOCK, a method to dock molecules that follow the protein's symmetry. For each prospective ligand pose, we apply the symmetry operation of the fibril to generate a self-interacting and fibril-interacting stack, checking that doing so will not cause a clash between the original molecule and its image. Absent a clash, we retain that pose and add the ligand-ligand van der Waals energy to the ligand's docking score (here using DOCK3.8). We can check these geometries and energies using an implementation of ANI, a neural-network-based quantum-mechanical evaluation of the ligand stacking energies. In retrospective calculations, symmetry docking can reproduce the poses of three tau PET tracers whose structures have been determined. More convincingly, in a prospective study, SymDOCK predicted the structure of the PET tracer MK-6240 bound in a symmetrical stack to AD PHF tau before that structure was determined; the docked pose was used to determine how MK-6240 fit the cryo-EM density. In proof-of-concept studies, SymDOCK enriched known ligands over property-matched decoys in retrospective screens without sacrificing docking speed and can address large library screens that seek new symmetrical stackers. Future applications of this approach will be considered.

Drug Discovery in the Field of β-Lactams: An Academic Perspective

β-Lactams are the most widely prescribed class of antibiotics that inhibit penicillin-binding proteins (PBPs), particularly transpeptidases that function in peptidoglycan synthesis. A major mechanism of antibiotic resistance is the production of β-lactamase enzymes, which are capable of hydrolyzing β-lactam antibiotics. There have been many efforts to counter increasing bacterial resistance against β-lactams. These studies have mainly focused on three areas: discovering novel inhibitors against β-lactamases, developing new β-lactams less susceptible to existing resistance mechanisms, and identifying non-β-lactam inhibitors against cell wall transpeptidases. Drug discovery in the β-lactam field has afforded a range of research opportunities for academia. In this review, we summarize the recent new findings on both β-lactamases and cell wall transpeptidases because these two groups of enzymes are evolutionarily and functionally connected. Many efforts to develop new β-lactams have aimed to inhibit both transpeptidases and β-lactamases, while several promising novel β-lactamase inhibitors have shown the potential to be further developed into transpeptidase inhibitors. In addition, the drug discovery progress against each group of enzymes is presented in three aspects: understanding the targets, screening methodology, and new inhibitor chemotypes. This is to offer insights into not only the advancement in this field but also the challenges, opportunities, and resources for future research. In particular, cyclic boronate compounds are now capable of inhibiting all classes of β-lactamases, while the diazabicyclooctane (DBO) series of small molecules has led to not only new β-lactamase inhibitors but potentially a new class of antibiotics by directly targeting PBPs. With the cautiously optimistic successes of a number of new β-lactamase inhibitor chemotypes and many questions remaining to be answered about the structure and function of cell wall transpeptidases, non-β-lactam transpeptidase inhibitors may usher in the next exciting phase of drug discovery in this field.

The Pivotal Distinction between Antagonists' and Agonists' Binding into Dopamine D4 Receptor-MD and FMO/PIEDA Studies

The dopamine D4 receptor (D4R) is a promising therapeutic target in widespread diseases, and the search for novel agonists and antagonists appears to be clinically relevant. The mechanism of binding to the receptor (R) for antagonists and agonists varies. In the present study, we conducted an in-depth computational study, teasing out key similarities and differences in binding modes, complex dynamics, and binding energies for D4R agonists and antagonists. The dynamic network method was applied to investigate the communication paths between the ligand (L) and G-protein binding site (GBS) of human D4R. Finally, the fragment molecular orbitals with pair interaction energy decomposition analysis (FMO/PIEDA) scheme was used to estimate the binding energies of L-R complexes. We found that a strong salt bridge with D3.32 initiates the inhibition of the dopamine D4 receptor. This interaction also occurs in the binding of agonists, but the change in the receptor conformation to the active state starts with interaction with cysteine C3.36. Such a mechanism may arise in the case of agonists unable to form a hydrogen bond with the serine S5.46, considered, so far, to be crucial in the activation of GPCRs. The energy calculations using the FMO/PIEDA method indicate that antagonists show higher residue occupancy of the receptor binding site than agonists, suggesting they could form relatively more stable complexes. Additionally, antagonists were characterized by repulsive interactions with S5.46 distinguishing them from agonists.

The ChEMBL Database in 2023: a drug discovery platform spanning multiple bioactivity data types and time periods

ChEMBL (https://www.ebi.ac.uk/chembl/) is a manually curated, high-quality, large-scale, open, FAIR and Global Core Biodata Resource of bioactive molecules with drug-like properties, previously described in the 2012, 2014, 2017 and 2019 Nucleic Acids Research Database Issues. Since its introduction in 2009, ChEMBL's content has changed dramatically in size and diversity of data types. Through incorporation of multiple new datasets from depositors since the 2019 update, ChEMBL now contains slightly more bioactivity data from deposited data vs data extracted from literature. In collaboration with the EUbOPEN consortium, chemical probe data is now regularly deposited into ChEMBL. Release 27 made curated data available for compounds screened for potential anti-SARS-CoV-2 activity from several large-scale drug repurposing screens. In addition, new patent bioactivity data have been added to the latest ChEMBL releases, and various new features have been incorporated, including a Natural Product likeness score, updated flags for Natural Products, a new flag for Chemical Probes, and the initial annotation of the action type for ∼270 000 bioactivity measurements.

DOCK 6: Incorporating hierarchical traversal through precomputed ligand conformations to enable large-scale docking

To allow DOCK 6 access to unprecedented chemical space for screening billions of small molecules, we have implemented features from DOCK 3.7 into DOCK 6, including a search routine that traverses precomputed ligand conformations stored in a hierarchical database. We tested them on the DUDE-Z and SB2012 test sets. The hierarchical database search routine is 16 times faster than anchor-and-grow. However, the ability of hierarchical database search to reproduce the experimental pose is 16% worse than that of anchor-and-grow. The enrichment performance is on average similar, but DOCK 3.7 has better enrichment than DOCK 6, and DOCK 6 is on average 1.7 times slower. However, with post-docking torsion minimization, DOCK 6 surpasses DOCK 3.7. A large-scale virtual screen is performed with DOCK 6 on 23 million fragment molecules. We use current features in DOCK 6 to complement hierarchical database calculations, including torsion minimization, which is not available in DOCK 3.7.

Targeting ion channels with ultra-large library screening for hit discovery

Ion channels play a crucial role in a variety of physiological and pathological processes, making them attractive targets for drug development in diseases such as diabetes, epilepsy, hypertension, cancer, and chronic pain. Despite the importance of ion channels in drug discovery, the vastness of chemical space and the complexity of ion channels pose significant challenges for identifying drug candidates. The use of in silico methods in drug discovery has dramatically reduced the time and cost of drug development and has the potential to revolutionize the field of medicine. Recent advances in computer hardware and software have enabled the screening of ultra-large compound libraries. Integration of different methods at various scales and dimensions is becoming an inevitable trend in drug development. In this review, we provide an overview of current state-of-the-art computational chemistry methodologies for ultra-large compound library screening and their application to ion channel drug discovery research. We discuss the advantages and limitations of various in silico techniques, including virtual screening, molecular mechanics/dynamics simulations, and machine learning-based approaches. We also highlight several successful applications of computational chemistry methodologies in ion channel drug discovery and provide insights into future directions and challenges in this field.

Exploring Superacid-Promoted Skeletal Reorganization of Aliphatic Nitrogen-Containing Compounds

Here we report a method to reorganize the core structure of aliphatic unsaturated nitrogen-containing substrates exploiting polyprotonation in superacid solutions. The superelectrophilic activation of N-isopropyl systems allows for the selective formal Csp3 -H activation/cyclization or homologation / functionalization of nitrogen-containing substrates. This study also reveals that this skeletal reorganization can be controlled through protonation interplay. The mechanism of this process involves an original sequence of C-N bond cleavage, isopropyl cation generation and subsequent C-N bond and C-C bond formation. This was demonstrated through in situ NMR analysis and labelling experiments, also confirmed by DFT calculations.

The rise of epitranscriptomics: recent developments and future directions

The epitranscriptomics field has undergone tremendous growth since the discovery that the RNA N6-methyladenosine (m6A) modification is reversible and is distributed throughout the transcriptome. Efforts to map RNA modifications transcriptome-wide and reshape the epitranscriptome in disease settings have facilitated mechanistic understanding and drug discovery in the field. In this review we discuss recent advancements in RNA modification detection methods and consider how these developments can be applied to gain novel insights into the epitranscriptome. We also highlight drug discovery efforts aimed at developing epitranscriptomic therapeutics for cancer and other diseases. Finally, we consider engineering of the epitranscriptome as an emerging direction to investigate RNA modifications and their causal effects on RNA processing at high specificity.

Discovery of 3-oxo-1,2,3,4-tetrahydropyrido[1,2-a]pyrazin derivatives as SARS-CoV-2 main protease inhibitors through virtual screening and biological evaluation

The COVID-19 caused by SARS-CoV-2 has led to a global pandemic that continues to impact societies and economies worldwide. The main protease (Mpro) plays a crucial role in SARS-CoV-2 replication and is an attractive target for anti-SARS-CoV-2 drug discovery. Herein, we report a series of 3-oxo-1,2,3,4-tetrahydropyrido[1,2-a]pyrazin derivatives as non-peptidomimetic inhibitors targeting SARS-CoV-2 Mpro through structure-based virtual screening and biological evaluation. Further similarity search and structure-activity relationship study led to the identification of compound M56-S2 with the enzymatic IC50 value of 4.0 μM. Moreover, the molecular simulation and predicted ADMET properties, indicated that non-peptidomimetic inhibitor M56-S2 might serve as a useful starting point for the further discovery of highly potent inhibitors targeting SARS-CoV-2 Mpro.

Design and Diversity Analysis of Chemical Libraries in Drug Discovery

Chemical libraries and compound data sets are among the main inputs to start the drug discovery process at universities, research institutes, and the pharmaceutical industry. The approach used in the design of compound libraries, the chemical information they possess, and the representation of structures, play a fundamental role in the development of studies: chemoinformatics, food informatics, in silico pharmacokinetics, computational toxicology, bioinformatics, and molecular modeling to generate computational hits that will continue the optimization process of drug candidates. The prospects for growth in drug discovery and development processes in chemical, biotechnological, and pharmaceutical companies began a few years ago by integrating computational tools with artificial intelligence methodologies. It is anticipated that it will increase the number of drugs approved by regulatory agencies shortly.

Phytochemicals Showing Antiangiogenic Effect in Pre-clinical Models and their Potential as an Alternative to Existing Therapeutics

Angiogenesis, the formation of new blood vessels from a pre-existing vascular network, is an important hallmark of several pathological conditions, such as tumor growth and metastasis, proliferative retinopathies, including proliferative diabetic retinopathy and retinopathy of prematurity, age-related macular degeneration, rheumatoid arthritis, psoriasis, and endometriosis. Putting a halt to pathology-driven angiogenesis is considered an important therapeutic strategy to slow down or reduce the severity of pathological disorders. Considering the attrition rate of synthetic antiangiogenic compounds from the lab to reaching the market due to severe side effects, several compounds of natural origin are being explored for their antiangiogenic properties. Employing pre-clinical models for the evaluation of novel antiangiogenic compounds is a promising strategy for rapid screening of antiangiogenic compounds. These studies use a spectrum of angiogenic model systems that include HUVEC two-dimensional culture, nude mice, chick chorioallantoic membrane, transgenic zebrafish, and dorsal aorta from rats and chicks, depending upon available resources. The present article emphasizes the antiangiogenic activity of the phytochemicals shown to exhibit antiangiogenic behavior in these well-defined existing angiogenic models and highlights key molecular targets. Different models help to get a quick understanding of the efficacy and therapeutics mechanism of emerging lead molecules. The inherent variability in assays and corresponding different phytochemicals tested in each study prevent their immediate utilization in clinical studies. This review will discuss phytochemicals discovered using suitable preclinical antiangiogenic models, along with a special mention of leads that have entered clinical evaluation.

HIt Discovery using docking ENriched by GEnerative Modeling (HIDDEN GEM): A novel computational workflow for accelerated virtual screening of ultra-large chemical libraries

Recent rapid expansion of make-on-demand, purchasable, chemical libraries comprising dozens of billions or even trillions of molecules has challenged the efficient application of traditional structure-based virtual screening methods that rely on molecular docking. We present a novel computational methodology termed HIDDEN GEM (HIt Discovery using Docking ENriched by GEnerative Modeling) that greatly accelerates virtual screening. This workflow uniquely integrates machine learning, generative chemistry, massive chemical similarity searching and molecular docking of small, selected libraries in the beginning and the end of the workflow. For each target, HIDDEN GEM nominates a small number of top-scoring virtual hits prioritized from ultra-large chemical libraries. We have benchmarked HIDDEN GEM by conducting virtual screening campaigns for 16 diverse protein targets using Enamine REAL Space library comprising 37 billion molecules. We show that HIDDEN GEM yields the highest enrichment factors as compared to state of the art accelerated virtual screening methods, while requiring the least computational resources. HIDDEN GEM can be executed with any docking software and employed by users with limited computational resources.

Applying atomistic neural networks to bias conformer ensembles towards bioactive-like conformations

Identifying bioactive conformations of small molecules is an essential process for virtual screening applications relying on three-dimensional structure such as molecular docking. For most small molecules, conformer generators retrieve at least one bioactive-like conformation, with an atomic root-mean-square deviation (ARMSD) lower than 1 Å, among the set of low-energy conformers generated. However, there is currently no general method to prioritise these likely target-bound conformations in the ensemble. In this work, we trained atomistic neural networks (AtNNs) on 3D information of generated conformers of a curated subset of PDBbind ligands to predict the ARMSD to their closest bioactive conformation, and evaluated the early enrichment of bioactive-like conformations when ranking conformers by AtNN prediction. AtNN ranking was compared with bioactivity-unaware baselines such as ascending Sage force field energy ranking, and a slower bioactivity-based baseline ranking by ascending Torsion Fingerprint Deviation to the Maximum Common Substructure to the most similar molecule in the training set (TFD2SimRefMCS). On test sets from random ligand splits of PDBbind, ranking conformers using ComENet, the AtNN encoding the most 3D information, leads to early enrichment of bioactive-like conformations with a median BEDROC of 0.29 ± 0.02, outperforming the best bioactivity-unaware Sage energy ranking baseline (median BEDROC of 0.18 ± 0.02), and performing on a par with the bioactivity-based TFD2SimRefMCS baseline (median BEDROC of 0.31 ± 0.02). The improved performance of the AtNN and TFD2SimRefMCS baseline is mostly observed on test set ligands that bind proteins similar to proteins observed in the training set. On a more challenging subset of flexible molecules, the bioactivity-unaware baselines showed median BEDROCs up to 0.02, while AtNNs and TFD2SimRefMCS showed median BEDROCs between 0.09 and 0.13. When performing rigid ligand re-docking of PDBbind ligands with GOLD using the 1% top-ranked conformers, ComENet ranked conformers showed a higher successful docking rate than bioactivity-unaware baselines, with a rate of 0.48 ± 0.02 compared to CSD probability baseline with a rate of 0.39 ± 0.02. Similarly, on a pharmacophore searching experiment, selecting the 20% top-ranked conformers ranked by ComENet showed higher hit rate compared to baselines. Hence, the approach presented here uses AtNNs successfully to focus conformer ensembles towards bioactive-like conformations, representing an opportunity to reduce computational expense in virtual screening applications on known targets that require input conformations.

Advances in ligand-specific biosensing for structurally similar molecules

The specificity of biological systems makes it possible to develop biosensors targeting specific metabolites, toxins, and pollutants in complex medical or environmental samples without interference from structurally similar compounds. For the last two decades, great efforts have been devoted to creating proteins or nucleic acids with novel properties through synthetic biology strategies. Beyond augmenting biocatalytic activity, expanding target substrate scopes, and enhancing enzymes' enantioselectivity and stability, an increasing research area is the enhancement of molecular specificity for genetically encoded biosensors. Here, we summarize recent advances in the development of highly specific biosensor systems and their essential applications. First, we describe the rational design principles required to create libraries containing potential mutants with less promiscuity or better specificity. Next, we review the emerging high-throughput screening techniques to engineer biosensing specificity for the desired target. Finally, we examine the computer-aided evaluation and prediction methods to facilitate the construction of ligand-specific biosensors.

Shape-Based Virtual Screening of a Billion-Compound Library Identifies Mycobacterial Lipoamide Dehydrogenase Inhibitors

Lpd (lipoamide dehydrogenase) in Mycobacterium tuberculosis (Mtb) is required for virulence and is a genetically validated tuberculosis (TB) target. Numerous screens have been performed over the last decade, yet only two inhibitor series have been identified. Recent advances in large-scale virtual screening methods combined with make-on-demand compound libraries have shown the potential for finding novel hits. In this study, the Enamine REAL library consisting of ∼1.12 billion compounds was efficiently screened using the GPU Shape screen method against Mtb Lpd to find additional chemical matter that would expand on the known sulfonamide inhibitor series. We identified six new inhibitors with IC50 in the range of 5-100 μM. While these compounds remained chemically close to the already known sulfonamide series inhibitors, some diversity was found in the cores of the hits. The two most potent hits were further validated by one-step potency optimization to submicromolar levels. The co-crystal structure of optimized analogue TDI-13537 provided new insights into the potency determinants of the series.

Docking for EP4R antagonists active against inflammatory pain

The lipid prostaglandin E2 (PGE2) mediates inflammatory pain by activating G protein-coupled receptors, including the prostaglandin E2 receptor 4 (EP4R). Nonsteroidal anti-inflammatory drugs (NSAIDs) reduce nociception by inhibiting prostaglandin synthesis, however, the disruption of upstream prostanoid biosynthesis can lead to pleiotropic effects including gastrointestinal bleeding and cardiac complications. In contrast, by acting downstream, EP4R antagonists may act specifically as anti-inflammatory agents and, to date, no selective EP4R antagonists have been approved for human use. In this work, seeking to diversify EP4R antagonist scaffolds, we computationally dock over 400 million compounds against an EP4R crystal structure and experimentally validate 71 highly ranked, de novo synthesized molecules. Further, we show how structure-based optimization of initial docking hits identifies a potent and selective antagonist with 16 nanomolar potency. Finally, we demonstrate favorable pharmacokinetics for the discovered compound as well as anti-allodynic and anti-inflammatory activity in several preclinical pain models in mice.

Structure-based discovery of novel P-glycoprotein inhibitors targeting the nucleotide binding domains

P-glycoprotein (P-gp), a membrane transport protein overexpressed in certain drug-resistant cancer cells, has been the target of numerous drug discovery projects aimed at overcoming drug resistance in cancer. Most characterized P-gp inhibitors bind at the large hydrophobic drug binding domain (DBD), but none have yet attained regulatory approval. In this study, we explored the potential of designing inhibitors that target the nucleotide binding domains (NBDs), by computationally screening a large library of 2.6 billion synthesizable molecules, using a combination of machine learning-guided molecular docking and molecular dynamics (MD). 14 of the computationally best-scoring molecules were subsequently tested for their ability to inhibit P-gp mediated calcein-AM efflux. In total, five diverse compounds exhibited inhibitory effects in the calcein-AM assay without displaying toxicity. The activity of these compounds was confirmed by their ability to decrease the verapamil-stimulated ATPase activity of P-gp in a subsequent assay. The discovery of these five novel P-gp inhibitors demonstrates the potential of in-silico screening in drug discovery and provides a new stepping point towards future potent P-gp inhibitors.

Recent advances in anti-inflammatory active components and action mechanisms of natural medicines

Inflammation is a series of reactions caused by the body's resistance to external biological stimuli. Inflammation affects the occurrence and development of many diseases. Anti-inflammatory drugs have been used widely to treat inflammatory diseases, but long-term use can cause toxic side-effects and affect human functions. As immunomodulators with long-term conditioning effects and no drug residues, natural products are being investigated increasingly for the treatment of inflammatory diseases. In this review, we focus on the inflammatory process and cellular mechanisms in the development of diseases such as inflammatory bowel disease, atherosclerosis, and coronavirus disease-2019. Also, we focus on three signaling pathways (Nuclear factor-kappa B, p38 mitogen-activated protein kinase, Janus kinase/signal transducer and activator of transcription-3) to explain the anti-inflammatory effect of natural products. In addition, we also classified common natural products based on secondary metabolites and explained the association between current bidirectional prediction progress of natural product targets and inflammatory diseases.

Advancing algorithmic drug product development: Recommendations for machine learning approaches in drug formulation

Artificial intelligence is a rapidly expanding area of research, with the disruptive potential to transform traditional approaches in the pharmaceutical industry, from drug discovery and development to clinical practice. Machine learning, a subfield of artificial intelligence, has fundamentally transformed in silico modelling and has the capacity to streamline clinical translation. This paper reviews data-driven modelling methodologies with a focus on drug formulation development. Despite recent advances, there is limited modelling guidance specific to drug product development and a trend towards suboptimal modelling practices, resulting in models that may not give reliable predictions in practice. There is an overwhelming focus on benchtop experimental outcomes obtained for a specific modelling aim, leaving the capabilities of data scraping or the use of combined modelling approaches yet to be fully explored. Moreover, the preference for high accuracy can lead to a reliance on black box methods over interpretable models. This further limits the widespread adoption of machine learning as black boxes yield models that cannot be easily understood for the purposes of enhancing product performance. In this review, recommendations for conducting machine learning research for drug product development to ensure trustworthiness, transparency, and reliability of the models produced are presented. Finally, possible future directions on how research in this area might develop are discussed to aim for models that provide useful and robust guidance to formulators.

Unrealized targets in the discovery of antibiotics for Gram-negative bacterial infections

Advances in areas that include genomics, systems biology, protein structure determination and artificial intelligence provide new opportunities for target-based antibacterial drug discovery. The selection of a 'good' new target for direct-acting antibacterial compounds is the first decision, for which multiple criteria must be explored, integrated and re-evaluated as drug discovery programmes progress. Criteria include essentiality of the target for bacterial survival, its conservation across different strains of the same species, bacterial species and growth conditions (which determines the spectrum of activity of a potential antibiotic) and the level of homology with human genes (which influences the potential for selective inhibition). Additionally, a bacterial target should have the potential to bind to drug-like molecules, and its subcellular location will govern the need for inhibitors to penetrate one or two bacterial membranes, which is a key challenge in targeting Gram-negative bacteria. The risk of the emergence of target-based drug resistance for drugs with single targets also requires consideration. This Review describes promising but as-yet-unrealized targets for antibacterial drugs against Gram-negative bacteria and examples of cognate inhibitors, and highlights lessons learned from past drug discovery programmes.

Streamlining pipeline efficiency: a novel model-agnostic technique for accelerating conditional generative and virtual screening pipelines

The discovery of potential therapeutic agents for life-threatening diseases has become a significant problem. There is a requirement for fast and accurate methods to identify drug-like molecules that can be used as potential candidates for novel targets. Existing techniques like high-throughput screening and virtual screening are time-consuming and inefficient. Traditional molecule generation pipelines are more efficient than virtual screening but use time-consuming docking software. Such docking functions can be emulated using Machine Learning models with comparable accuracy and faster execution times. However, we find that when pre-trained machine learning models are employed in generative pipelines as oracles, they suffer from model degradation in areas where data is scarce. In this study, we propose an active learning-based model that can be added as a supplement to enhanced molecule generation architectures. The proposed method uses uncertainty sampling on the molecules created by the generator model and dynamically learns as the generator samples molecules from different regions of the chemical space. The proposed framework can generate molecules with high binding affinity with [Formula: see text]a 70% improvement in runtime compared to the baseline model by labeling only [Formula: see text]30% of molecules compared to the baseline oracle.

FASTDock: A Pipeline for Allosteric Drug Discovery

Allostery is involved in innumerable biological processes and plays a fundamental role in human disease. Thus, the exploration of allosteric modulation is crucial for research on biological mechanisms and in the development of novel therapeutics. The development of small-molecule allosteric effectors can be used as tools to probe biological mechanisms of interest. One of the main limitations in targeting allosteric sites is the difficulty in uncovering them for specific receptors. Furthermore, upon discovery of novel allosteric modulation, early lead generation is made more difficult as compared to that at orthosteric sites because there is likely no information about the types of molecules that can bind at the site. In the work described here, we present a novel drug discovery pipeline, FASTDock, which allows one to uncover ligandable sites as well as small molecules that target the given site without requiring pre-existing knowledge of ligands that can bind in the targeted site. By using a hierarchical screening strategy, this method has the potential to enable high-throughput screens of an exceptionally large database of targeted ligand space.

Molecular basis of opioid receptor signaling

Opioids are used for pain management despite the side effects that contribute to the opioid crisis. The pursuit of non-addictive opioid analgesics remains unattained due to the unresolved intricacies of opioid actions, receptor signaling cascades, and neuronal plasticity. Advancements in structural, molecular, and computational tools illuminate the dynamic interplay between opioids and opioid receptors, as well as the molecular determinants of signaling pathways, which are potentially interlinked with pharmacological responses. Here, we review the molecular basis of opioid receptor signaling with a focus on the structures of opioid receptors bound to endogenous peptides or pharmacological agents. These insights unveil specific interactions that dictate ligand selectivity and likely their distinctive pharmacological profiles. Biochemical analysis further unveils molecular features governing opioid receptor signaling. Simultaneously, the synergy between computational biology and medicinal chemistry continues to expedite the discovery of novel chemotypes with the promise of yielding more efficacious and safer opioid compounds.

Integrated Molecular Modeling and Machine Learning for Drug Design

Modern therapeutic development often involves several stages that are interconnected, and multiple iterations are usually required to bring a new drug to the market. Computational approaches have increasingly become an indispensable part of helping reduce the time and cost of the research and development of new drugs. In this Perspective, we summarize our recent efforts on integrating molecular modeling and machine learning to develop computational tools for modulator design, including a pocket-guided rational design approach based on AlphaSpace to target protein-protein interactions, delta machine learning scoring functions for protein-ligand docking as well as virtual screening, and state-of-the-art deep learning models to predict calculated and experimental molecular properties based on molecular mechanics optimized geometries. Meanwhile, we discuss remaining challenges and promising directions for further development and use a retrospective example of FDA approved kinase inhibitor Erlotinib to demonstrate the use of these newly developed computational tools.

Deep Learning with Geometry-Enhanced Molecular Representation for Augmentation of Large-Scale Docking-Based Virtual Screening

Structure-based virtual screening has been a crucial tool in drug discovery for decades. However, as the chemical space expands, the existing structure-based virtual screening techniques based on molecular docking and scoring struggle to handle billion-entry ultralarge libraries due to the high computational cost. To address this challenge, people have resorted to machine learning techniques to enhance structure-based virtual screening for efficiently exploring the vast chemical space. In those cases, compounds are usually treated as sequential strings or two-dimensional topology graphs, limiting their ability to incorporate three-dimensional structural information for downstream tasks. We herein propose a novel deep learning protocol, GEM-Screen, which utilizes the geometry-enhanced molecular representation of the compounds docking to a specific target and is trained on docking scores of a small fraction of a library through an active learning strategy to approximate the docking outcome for yet nontraining entries. This protocol is applied to virtual screening campaigns against the AmpC and D4 targets, demonstrating that GEM-Screen enriches more than 90% of the hit scaffolds for AmpC in the top 4% of model predictions and more than 80% of the hit scaffolds for D4 in the same top-ranking size of library. GEM-Screen can be used in conjunction with traditional docking programs for docking of only the top-ranked compounds to avoid the exhaustive docking of the whole library, thus allowing for discovering top-scoring compounds from billion-entry libraries in a rapid yet accurate fashion.

A flexible data-free framework for structure-based de novo drug design with reinforcement learning

Contemporary structure-based molecular generative methods have demonstrated their potential to model the geometric and energetic complementarity between ligands and receptors, thereby facilitating the design of molecules with favorable binding affinity and target specificity. Despite the introduction of deep generative models for molecular generation, the atom-wise generation paradigm that partially contradicts chemical intuition limits the validity and synthetic accessibility of the generated molecules. Additionally, the dependence of deep learning models on large-scale structural data has hindered their adaptability across different targets. To overcome these challenges, we present a novel search-based framework, 3D-MCTS, for structure-based de novo drug design. Distinct from prevailing atom-centric methods, 3D-MCTS employs a fragment-based molecular editing strategy. The fragments decomposed from small-molecule drugs are recombined under predefined retrosynthetic rules, offering improved drug-likeness and synthesizability, overcoming the inherent limitations of atom-based approaches. Leveraging multi-threaded parallel simulations combined with a real-time energy constraint-based pruning strategy, 3D-MCTS achieves remarkable efficiency. At a fixed computational cost, it outperforms other state-of-the-art (SOTA) methods by producing molecules with enhanced binding affinity. Furthermore, its fragment-based approach ensures the generation of more dependable binding conformations, exhibiting a success rate 43.6% higher than that of other SOTAs. This advantage becomes even more pronounced when handling targets that significantly deviate from the training dataset. 3D-MCTS is capable of achieving thirty times more hits with high binding affinity than traditional virtual screening methods, which demonstrates the superior ability of 3D-MCTS to explore chemical space. Moreover, the flexibility of our framework makes it easy to incorporate domain knowledge during the process, thereby enabling the generation of molecules with desirable pharmacophores and enhanced binding affinity. The adaptability of 3D-MCTS is further showcased in metalloprotein applications, highlighting its potential across various drug design scenarios.

Drug discovery by targeting the protein-protein interactions involved in autophagy

Autophagy is a cellular process in which proteins and organelles are engulfed in autophagosomal vesicles and transported to the lysosome/vacuole for degradation. Protein-protein interactions (PPIs) play a crucial role at many stages of autophagy, which present formidable but attainable targets for autophagy regulation. Moreover, selective regulation of PPIs tends to have a lower risk in causing undesired off-target effects in the context of a complicated biological network. Thus, small-molecule regulators, including peptides and peptidomimetics, targeting the critical PPIs involved in autophagy provide a new opportunity for innovative drug discovery. This article provides general background knowledge of the critical PPIs involved in autophagy and reviews a range of successful attempts on discovering regulators targeting those PPIs. Successful strategies and existing limitations in this field are also discussed.

Improving drug discovery with a hybrid deep generative model using reinforcement learning trained on a Bayesian docking approximation

Generative approaches to molecular design are an area of intense study in recent years as a method to generate new pharmaceuticals with desired properties. Often though, these types of efforts are constrained by limited experimental activity data, resulting in either models that generate molecules with poor performance or models that are overfit and produce close analogs of known molecules. In this paper, we reduce this data dependency for the generation of new chemotypes by incorporating docking scores of known and de novo molecules to expand the applicability domain of the reward function and diversify the compounds generated during reinforcement learning. Our approach employs a deep generative model initially trained using a combination of limited known drug activity and an approximate docking score provided by a second machine learned Bayes regression model, with final evaluation of high scoring compounds by a full docking simulation. This strategy results in molecules with docking scores improved by 10-20% compared to molecules of similar size, while being 130 × faster than a docking only approach on a typical GPU workstation. We also show that the increased docking scores correlate with (1) docking poses with interactions similar to known inhibitors and (2) result in higher MM-GBSA binding energies comparable to the energies of known DDR1 inhibitors, demonstrating that the Bayesian model contains sufficient information for the network to learn to efficiently interact with the binding pocket during reinforcement learning. This outcome shows that the combination of the learned latent molecular representation along with the feature-based docking regression is sufficient for reinforcement learning to infer the relationship between the molecules and the receptor binding site, which suggest that our method can be a powerful tool for the discovery of new chemotypes with potential therapeutic applications.

High-throughput screening of ancient forest plant extracts shows cytotoxicity towards triple-negative breast cancer

According to the World Health Organization, women's breast cancer is among the most common cancers with 7.8 million diagnosed cases during 2016-2020 and encompasses 15 % of all female cancer-related mortalities. These mortality events from triple-negative breast cancer are a significant health issue worldwide calling for a continuous search of bioactive compounds for better cancer treatments. Historically, plants are important sources for identifying such new bioactive chemicals for treatments. Here we use high-throughput screening and mass spectrometry analyses of extracts from 100 plant species collected in Chinese ancient forests to detect novel bioactive breast cancer phytochemicals. First, to study the effects on viability of the plant extracts, we used a MTT and CCK-8 cytotoxicity assay employing triple-negative breast cancer (TNBC) MDA-MB-231 and normal epithelial MCF-10A cell lines and cell cycle arrest to estimate apoptosis using flow cytometry for the most potent three speices. Based on these analyses, the final most potent extracts were from the Amur honeysuckle (Lonicera maackii) wood/root bark and Nigaki (Picrasma quassioides) wood/root bark. Then, 5 × 106 MDA-MB-231 cells were injected subcutaneously into the right hind leg of nude mice and a tumour was allowed to grow before treatment for seven days. Subsequently, the four exposed groups received gavage extracts from Amur honeysuckle and Nigaki (Amur honeysuckle wood distilled water, Amur honeysuckle root bark ethanol, Nigaki wood ethanol or Nigaki root bark distilled water/ethanol (1:1) extracts) in phosphate-buffered saline (PBS), while the control group received only PBS. The tumour weight of treated nude mice was reduced significantly by 60.5 % within 2 weeks, while on average killing 70 % of the MDA-MB-231 breast cancer cells after 48 h treatment (MTT test). In addition, screening of target genes using the Swiss Target Prediction, STITCH, STRING and NCBI-gene database showed that the four plant extracts possess desirable activity towards several known breast cancer genes. This reflects that the extracts may kill MBD-MB-231 breast cancer cells. This is the first screening of plant extracts with high efficiency in 2 decades, showing promising results for future development of novel cancer treatments.

EasyDock: customizable and scalable docking tool

Docking of large compound collections becomes an important procedure to discover new chemical entities. Screening of large sets of compounds may also occur in de novo design projects guided by molecular docking. To facilitate these processes, there is a need for automated tools capable of efficiently docking a large number of molecules using multiple computational nodes within a reasonable timeframe. These tools should also allow for easy integration of new docking programs and provide a user-friendly program interface to support the development of further approaches utilizing docking as a foundation. Currently available tools have certain limitations, such as lacking a convenient program interface or lacking support for distributed computations. In response to these limitations, we have developed a module called EasyDock. It can be deployed over a network of computational nodes using the Dask library, without requiring a specific cluster scheduler. Furthermore, we have proposed and implemented a simple model that predicts the runtime of docking experiments and applied it to minimize overall docking time. The current version of EasyDock supports popular docking programs, namely Autodock Vina, gnina, and smina. Additionally, we implemented a supplementary feature to enable docking of boron-containing compounds, which are not inherently supported by Vina and smina, and demonstrated its applicability on a set of 55 PDB protein-ligand complexes.

A practical guide to machine-learning scoring for structure-based virtual screening

Structure-based virtual screening (SBVS) via docking has been used to discover active molecules for a range of therapeutic targets. Chemical and protein data sets that contain integrated bioactivity information have increased both in number and in size. Artificial intelligence and, more concretely, its machine-learning (ML) branch, including deep learning, have effectively exploited these data sets to build scoring functions (SFs) for SBVS against targets with an atomic-resolution 3D model (e.g., generated by X-ray crystallography or predicted by AlphaFold2). Often outperforming their generic and non-ML counterparts, target-specific ML-based SFs represent the state of the art for SBVS. Here, we present a comprehensive and user-friendly protocol to build and rigorously evaluate these new SFs for SBVS. This protocol is organized into four sections: (i) using a public benchmark of a given target to evaluate an existing generic SF; (ii) preparing experimental data for a target from public repositories; (iii) partitioning data into a training set and a test set for subsequent target-specific ML modeling; and (iv) generating and evaluating target-specific ML SFs by using the prepared training-test partitions. All necessary code and input/output data related to three example targets (acetylcholinesterase, HMG-CoA reductase, and peroxisome proliferator-activated receptor-α) are available at https://github.com/vktrannguyen/MLSF-protocol , can be run by using a single computer within 1 week and make use of easily accessible software/programs (e.g., Smina, CNN-Score, RF-Score-VS and DeepCoy) and web resources. Our aim is to provide practical guidance on how to augment training data to enhance SBVS performance, how to identify the most suitable supervised learning algorithm for a data set, and how to build an SF with the highest likelihood of discovering target-active molecules within a given compound library.

Docking for molecules that bind in a symmetric stack with SymDOCK

Discovering ligands for amyloid fibrils, such as those formed by the tau protein, is an area of much current interest. In recent structures, ligands bind in stacks in the tau fibrils to reflect the rotational and translational symmetry of the fibril itself; in these structures the ligands make few interactions with the protein but interact extensively with each other. To exploit this symmetry and stacking, we developed SymDOCK, a method to dock molecules that follow the protein's symmetry. For each prospective ligand pose, we apply the symmetry operation of the fibril to generate a self-interacting and fibril-interacting stack, checking that doing so will not cause a clash between the original molecule and its image. Absent a clash, we retain that pose and add the ligand-ligand van der Waals energy to the ligand's docking score (here using DOCK3.8). We can check these geometries and energies using an implementation of ANI, a neural network-based quantum-mechanical evaluation of the ligand stacking energies. In retrospective calculations, symmetry docking can reproduce the poses of three tau PET tracers whose structures have been determined. More convincingly, in a prospective study SymDOCK predicted the structure of the PET tracer MK-6240 bound in a symmetrical stack to AD PHF tau before that structure was determined; the docked pose was used to determine how MK-6240 fit the cryo-EM density. In proof-of-concept studies, SymDOCK enriched known ligands over property-matched decoys in retrospective screens without sacrificing docking speed, and can address large library screens that seek new symmetrical stackers. Future applications of this approach will be considered.

Precision Oncology Comes of Age: Designing Best-in-Class Small Molecules by Integrating Two Decades of Advances in Chemistry, Target Biology, and Data Science

Small-molecule drugs have enabled the practice of precision oncology for genetically defined patient populations since the first approval of imatinib in 2001. Scientific and technology advances over this 20-year period have driven the evolution of cancer biology, medicinal chemistry, and data science. Collectively, these advances provide tools to more consistently design best-in-class small-molecule drugs against known, previously undruggable, and novel cancer targets. The integration of these tools and their customization in the hands of skilled drug hunters will be necessary to enable the discovery of transformational therapies for patients across a wider spectrum of cancers.

SIGNIFICANCE

Target-centric small-molecule drug discovery necessitates the consideration of multiple approaches to identify chemical matter that can be optimized into drug candidates. To do this successfully and consistently, drug hunters require a comprehensive toolbox to avoid following the "law of instrument" or Maslow's hammer concept where only one tool is applied regardless of the requirements of the task. Combining our ever-increasing understanding of cancer and cancer targets with the technological advances in drug discovery described below will accelerate the next generation of small-molecule drugs in oncology.

Functional screening and rational design of compounds targeting GPR132 to treat diabetes

Chronic inflammation due to islet-residing macrophages plays key roles in the development of type 2 diabetes mellitus. By systematically profiling intra-islet lipid-transmembrane receptor signalling in islet-resident macrophages, we identified endogenous 9(S)-hydroxy-10,12-octadecadienoic acid-G-protein-coupled receptor 132 (GPR132)-Gi signalling as a significant contributor to islet macrophage reprogramming and found that GPR132 deficiency in macrophages reversed metabolic disorders in mice fed a high-fat diet. The cryo-electron microscopy structures of GPR132 bound with two endogenous agonists, N-palmitoylglycine and 9(S)-hydroxy-10,12-octadecadienoic acid, enabled us to rationally design both GPR132 agonists and antagonists with high potency and selectivity through stepwise translational approaches. We ultimately identified a selective GPR132 antagonist, NOX-6-18, that modulates macrophage reprogramming within pancreatic islets, decreases weight gain and enhances glucose metabolism in mice fed a high-fat diet. Our study not only illustrates that intra-islet lipid signalling contributes to islet macrophage reprogramming but also provides a broadly applicable strategy for the identification of important G-protein-coupled receptor targets in pathophysiological processes, followed by the rational design of therapeutic leads for refractory diseases such as diabetes.

Fake it until you make it? Generative de novo design and virtual screening of synthesizable molecules

Computational techniques, including virtual screening, de novo design, and generative models, play an increasing role in expediting DMTA cycles for modern molecular discovery. However, computationally proposed molecules must be synthetically feasible for laboratory testing. In this perspective, we offer a succinct introduction to the subject, and showcase typical workflows to integrate synthesis planning, synthesizability scoring, and molecule generation. Finally, we address limitations and opportunities for future research.

ViTScore: A Novel Three-Dimensional Vision Transformer Method for Accurate Prediction of Protein-Ligand Docking Poses

Protein-ligand interactions (PLIs) are essential for cellular activities and drug discovery, and due to the complexity and high cost of experimental methods, there is a great demand for computational approaches, such as protein-ligand docking, to decipher PLI patterns. One of the most challenging aspects of protein-ligand docking is to identify near-native conformations from a set of poses, but traditional scoring functions still have limited accuracy. Therefore, new scoring methods are urgently needed for methodological and/or practical implications. We present a novel deep learning-based scoring function for ranking protein-ligand docking poses based on Vision Transformer (ViT), named ViTScore. To recognize near-native poses from a set of poses, ViTScore voxelizes the protein-ligand interactional pocket into a 3D grid labeled by the occupancy contribution of atoms in different physicochemical classes. This allows ViTScore to capture the subtle differences between spatially and energetically favorable near-native poses and unfavorable non-native poses without needing extra information. After that, ViTScore will output the prediction of the root mean square deviation (rmsd) of a docking pose with reference to the native binding pose. ViTScore is extensively evaluated on diverse test sets including PDBbind2019 and CASF2016, and obtains significant improvements over existing methods in terms of RMSE, R and docking power. Moreover, the results demonstrate that ViTScore is a promising scoring function for protein-ligand docking, and it can be used to accurately identify near-native poses from a set of poses. Furthermore, the results suggest that ViTScore is a powerful tool for protein-ligand docking, and it can be used to accurately identify near-native poses from a set of poses. Additionally, ViTScore can be used to identify potential drug targets and to design new drugs with improved efficacy and safety.