Yes SIR! On the structure-inactivity relationships in drug discovery

Regiodivergent Radical-Relay Alkene Dicarbofunctionalization

Herein, we report a regiodivergent 1,n-dicarbofunctionalization of unactivated olefins enabled by a Ni-catalyzed radical relay that forges both C(sp3)-C(sp3) and C(sp2)-C(sp3) linkages, even at long-range. Initial studies support an intertwined scenario resulting from the merger of an atom-transfer radical addition (ATRA) and a chain-walking event, with site-selectivity being dictated by a judicious choice of the ligand backbone.

Toward structure-multiple activity relationships (SMARts) using computational approaches: A polypharmacological perspective

In the current era of biological big data, which are rapidly populating the biological chemical space, in silico polypharmacology drug design approaches help to decode structure-multiple activity relationships (SMARts). Current computational methods can predict or categorize multiple properties simultaneously, which aids the generation, identification, curation, prioritization, optimization, and repurposing of molecules. Computational methods have generated opportunities and challenges in medicinal chemistry, pharmacology, food chemistry, toxicology, bioinformatics, and chemoinformatics. It is anticipated that computer-guided SMARts could contribute to the full automatization of drug design and drug repurposing campaigns, facilitating the prediction of new biological targets, side and off-target effects, and drug-drug interactions.

De Novo Design of Inhibitors of DNA Methyltransferase 1: A Critical Comparison of Ligand- and Structure-Based Approaches

Designing and developing inhibitors against the epigenetic target DNA methyltransferase (DNMT) is an attractive strategy in epigenetic drug discovery. DNMT1 is one of the epigenetic enzymes with significant clinical relevance. Structure-based de novo design is a drug discovery strategy that was used in combination with similarity searching to identify a novel DNMT inhibitor with a novel chemical scaffold and warrants further exploration. This study aimed to continue exploring the potential of de novo design to build epigenetic-focused libraries targeted toward DNMT1. Herein, we report the results of an in-depth and critical comparison of ligand- and structure-based de novo design of screening libraries focused on DNMT1. The newly designed chemical libraries focused on DNMT1 are freely available on GitHub.

KLSD: a kinase database focused on ligand similarity and diversity

Due to the similarity and diversity among kinases, small molecule kinase inhibitors (SMKIs) often display multi-target effects or selectivity, which have a strong correlation with the efficacy and safety of these inhibitors. However, due to the limited number of well-known popular databases and their restricted data mining capabilities, along with the significant scarcity of databases focusing on the pharmacological similarity and diversity of SMIKIs, researchers find it challenging to quickly access relevant information. The KLIFS database is representative of specialized application databases in the field, focusing on kinase structure and co-crystallised kinase-ligand interactions, whereas the KLSD database in this paper emphasizes the analysis of SMKIs among all reported kinase targets. To solve the current problem of the lack of professional application databases in kinase research and to provide centralized, standardized, reliable and efficient data resources for kinase researchers, this paper proposes a research program based on the ChEMBL database. It focuses on kinase ligands activities comparisons. This scheme extracts kinase data and standardizes and normalizes them, then performs kinase target difference analysis to achieve kinase activity threshold judgement. It then constructs a specialized and personalized kinase database platform, adopts the front-end and back-end separation technology of SpringBoot architecture, constructs an extensible WEB application, handles the storage, retrieval and analysis of the data, ultimately realizing data visualization and interaction. This study aims to develop a kinase database platform to collect, organize, and provide standardized data related to kinases. By offering essential resources and tools, it supports kinase research and drug development, thereby advancing scientific research and innovation in kinase-related fields. It is freely accessible at: http://ai.njucm.edu.cn:8080.

Benchmarking compound activity prediction for real-world drug discovery applications

Identifying active compounds for target proteins is fundamental in early drug discovery. Recently, data-driven computational methods have demonstrated promising potential in predicting compound activities. However, there lacks a well-designed benchmark to comprehensively evaluate these methods from a practical perspective. To fill this gap, we propose a Compound Activity benchmark for Real-world Applications (CARA). Through carefully distinguishing assay types, designing train-test splitting schemes and selecting evaluation metrics, CARA can consider the biased distribution of current real-world compound activity data and avoid overestimation of model performances. We observed that although current models can make successful predictions for certain proportions of assays, their performances varied across different assays. In addition, evaluation of several few-shot training strategies demonstrated different performances related to task types. Overall, we provide a high-quality dataset for developing and evaluating compound activity prediction models, and the analyses in this work may inspire better applications of data-driven models in drug discovery.

Benford's Law and distributions for better drug design

INTRODUCTION

Modern drug discovery incorporates various tools and data, heralding the beginning of the data-driven drug design (DD) era. The distributions of chemical and physical data used for Artificial Intelligence (AI)/Machine Learning (ML) and to drive DD have thus become highly important to be understood and used effectively.

AREAS COVERED

The authors perform a comprehensive exploration of the statistical distributions driving the data-intensive era of drug discovery, including Benford's Law in AI/ML-based DD.

EXPERT OPINION

As the relevance of data-driven discovery escalates, we anticipate meticulous scrutiny of datasets utilizing principles like Benford's Law to enhance data integrity and guide efficient resource allocation and experimental planning. In this data-driven era of the pharmaceutical and medical industries, addressing critical aspects such as bias mitigation, algorithm effectiveness, data stewardship, effects, and fraud prevention are essential. Harnessing Benford's Law and other distributions and statistical tests in DD provides a potent strategy to detect data anomalies, fill data gaps, and enhance dataset quality. Benford's Law is a fast method for data integrity and quality of datasets, the backbone of AI/ML and other modeling approaches, proving very useful in the design process.

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.

Indole-Based Compounds as Potential Drug Candidates for SARS-CoV-2

The COVID-19 pandemic has posed a significant threat to society in recent times, endangering human health, life, and economic well-being. The disease quickly spreads due to the highly infectious SARS-CoV-2 virus, which has undergone numerous mutations. Despite intense research efforts by the scientific community since its emergence in 2019, no effective therapeutics have been discovered yet. While some repurposed drugs have been used to control the global outbreak and save lives, none have proven universally effective, particularly for severely infected patients. Although the spread of the disease is generally under control, anti-SARS-CoV-2 agents are still needed to combat current and future infections. This study reviews some of the most promising repurposed drugs containing indolyl heterocycle, which is an essential scaffold of many alkaloids with diverse bio-properties in various biological fields. The study also discusses natural and synthetic indole-containing compounds with anti-SARS-CoV-2 properties and computer-aided drug design (in silico studies) for optimizing anti-SARS-CoV-2 hits/leads.

Exploring chemical space for lead identification by propagating on chemical similarity network

Motivation

Lead identification is a fundamental step to prioritize candidate compounds for downstream drug discovery process. Machine learning (ML) and deep learning (DL) approaches are widely used to identify lead compounds using both chemical property and experimental information. However, ML or DL methods rarely consider compound similarity information directly since ML and DL models use abstract representation of molecules for model construction. Alternatively, data mining approaches are also used to explore chemical space with drug candidates by screening undesirable compounds. A major challenge for data mining approaches is to develop efficient data mining methods that search large chemical space for desirable lead compounds with low false positive rate.

Results

In this work, we developed a network propagation (NP) based data mining method for lead identification that performs search on an ensemble of chemical similarity networks. We compiled 14 fingerprint-based similarity networks. Given a target protein of interest, we use a deep learning-based drug target interaction model to narrow down compound candidates and then we use network propagation to prioritize drug candidates that are highly correlated with drug activity score such as IC50. In an extensive experiment with BindingDB, we showed that our approach successfully discovered intentionally unlabeled compounds for given targets. To further demonstrate the prediction power of our approach, we identified 24 candidate leads for CLK1. Two out of five synthesizable candidates were experimentally validated in binding assays. In conclusion, our framework can be very useful for lead identification from very large compound databases such as ZINC.

Diversity and Chemical Space Characterization of Inhibitors of the Epigenetic Target G9a: A Chemoinformatics Approach

G9a is a histone-lysine methyltransferase that performs the mono- and dimethylation of lysine 9 at histone 3 of the nucleosome. It belongs to the SET PKMT family, and its methylations are related to promoter repression and activation. G9a is a promising epigenetic target. Despite the fact that there are several G9a inhibitors under development, there are no compounds in clinical use due to adverse in vivo ADMET (absorption, distribution, metabolism, excretion, and toxicity) issues. The goal of this study is to discuss the exploration, characterization, and analysis of the chemical space of 409 G9a inhibitors reported in a large public database. Exploring the chemical space of the inhibitors led to the quantification of their structural diversity based on molecular scaffolds and structural fingerprints of different designs. As part of the analysis, the G9a inhibitors were compared with commercial libraries focused on epigenetic targets. The findings of this work will help in the development of, in a follow-up study, predictive models to identify G9a inhibitors. This study also points out the relevance of screening commercial libraries to expand the epigenetic relevant chemical space, in particular, G9a inhibitors.

Virtual Screening Strategy to Identify Retinoic Acid-Related Orphan Receptor γt Modulators

Molecular docking is a key method used in virtual screening (VS) campaigns to identify small-molecule ligands for drug discovery targets. While docking provides a tangible way to understand and predict the protein-ligand complex formation, the docking algorithms are often unable to separate active ligands from inactive molecules in practical VS usage. Here, a novel docking and shape-focused pharmacophore VS protocol is demonstrated for facilitating effective hit discovery using retinoic acid receptor-related orphan receptor gamma t (RORγt) as a case study. RORγt is a prospective target for treating inflammatory diseases such as psoriasis and multiple sclerosis. First, a commercial molecular database was flexibly docked. Second, the alternative docking poses were rescored against the shape/electrostatic potential of negative image-based (NIB) models that mirror the target's binding cavity. The compositions of the NIB models were optimized via iterative trimming and benchmarking using a greedy search-driven algorithm or brute force NIB optimization. Third, a pharmacophore point-based filtering was performed to focus the hit identification on the known RORγt activity hotspots. Fourth, free energy binding affinity evaluation was performed on the remaining molecules. Finally, twenty-eight compounds were selected for in vitro testing and eight compounds were determined to be low μM range RORγt inhibitors, thereby showing that the introduced VS protocol generated an effective hit rate of ~29%.

Towards Decoding Hepatotoxicity of Approved Drugs through Navigation of Multiverse and Consensus Chemical Spaces

Drug-induced liver injury (DILI) is the principal reason for failure in developing drug candidates. It is the most common reason to withdraw from the market after a drug has been approved for clinical use. In this context, data from animal models, liver function tests, and chemical properties could complement each other to understand DILI events better and prevent them. Since the chemical space concept improves decision-making drug design related to the prediction of structure-property relationships, side effects, and polypharmacology drug activity (uniquely mentioning the most recent advances), it is an attractive approach to combining different phenomena influencing DILI events (e.g., individual "chemical spaces") and exploring all events simultaneously in an integrated analysis of the DILI-relevant chemical space. However, currently, no systematic methods allow the fusion of a collection of different chemical spaces to collect different types of data on a unique chemical space representation, namely "consensus chemical space." This study is the first report that implements data fusion to consider different criteria simultaneously to facilitate the analysis of DILI-related events. In particular, the study highlights the importance of analyzing together in vitro and chemical data (e.g., topology, bond order, atom types, presence of rings, ring sizes, and aromaticity of compounds encoded on RDKit fingerprints). These properties could be aimed at improving the understanding of DILI events.

Synthesis of non-symmetric N-benzylbispidinol amides and study of their inhibitory activity against the main protease of the SARS-CoV-2 virus

Based on the data obtained by molecular modeling of the non-covalent interaction of non-symmetric N-benzylbispidin-9-ol amides with the active site of the main protease 3CLpro of the SARS-CoV-2 virus, a series of compounds was synthesized, and their inhibitory activity against 3CLpro was studied and compared with that of the known inhibitor ML188 (IC₅₀ = 1.56±0.55 µmol L^-1). It was found that only compound 1g containing the 1,4-dihydroindeno[1,2-c]pyrazole fragment showed moderate activity (IC₅₀ = 100±5.7µmol L^-1) and was characterized by the highest calculated binding energy among the studied bispidine derivatives according to molecular docking data.

Computer Science Technology in Natural Products Research: A Review of Its Applications and Implications

Computational approaches to drug development are rapidly growing in popularity and have been used to produce significant results. Recent developments in information science have expanded databases and chemical informatics knowledge relating to natural products. Natural products have long been well-studied, and a large number of unique structures and remarkable active substances have been reported. Analyzing accumulated natural product knowledge using emerging computational science techniques is expected to yield more new discoveries. In this article, we discuss the current state of natural product research using machine learning. The basic concepts and frameworks of machine learning are summarized. Natural product research that utilizes machine learning is described in terms of the exploration of active compounds, automatic compound design, and application to spectral data. In addition, efforts to develop drugs for intractable diseases will be addressed. Lastly, we discuss key considerations for applying machine learning in this field. This paper aims to promote progress in natural product research by presenting the current state of computational science and chemoinformatics approaches in terms of its applications, strengths, limitations, and implications for the field.

Toward Quantitative Models in Safety Assessment: A Case Study to Show Impact of Dose-Response Inference on hERG Inhibition Models

Due to challenges with historical data and the diversity of assay formats, in silico models for safety-related endpoints are often based on discretized data instead of the data on a natural continuous scale. Models for discretized endpoints have limitations in usage and interpretation that can impact compound design. Here, we present a consistent data inference approach, exemplified on two data sets of Ether-à-go-go-Related Gene (hERG) K+ inhibition data, for dose-response and screening experiments that are generally applicable for in vitro assays. hERG inhibition has been associated with severe cardiac effects and is one of the more prominent safety targets assessed in drug development, using a wide array of in vitro and in silico screening methods. In this study, the IC₅₀ for hERG inhibition is estimated from diverse historical proprietary data. The IC₅₀ derived from a two-point proprietary screening data set demonstrated high correlation (R = 0.98, MAE = 0.08) with IC_50s derived from six-point dose-response curves. Similar IC₅₀ estimation accuracy was obtained on a public thallium flux assay data set (R = 0.90, MAE = 0.2). The IC₅₀ data were used to develop a robust quantitative model. The model's MAE (0.47) and R² (0.46) were on par with literature statistics and approached assay reproducibility. Using a continuous model has high value for pharmaceutical projects, as it enables rank ordering of compounds and evaluation of compounds against project-specific inhibition thresholds. This data inference approach can be widely applicable to assays with quantitative readouts and has the potential to impact experimental design and improve model performance, interpretation, and acceptance across many standard safety endpoints.

Chemoinformatics and artificial intelligence colloquium: progress and challenges in developing bioactive compounds

We report the main conclusions of the first Chemoinformatics and Artificial Intelligence Colloquium, Mexico City, June 15-17, 2022. Fifteen lectures were presented during a virtual public event with speakers from industry, academia, and non-for-profit organizations. Twelve hundred and ninety students and academics from more than 60 countries. During the meeting, applications, challenges, and opportunities in drug discovery, de novo drug design, ADME-Tox (absorption, distribution, metabolism, excretion and toxicity) property predictions, organic chemistry, peptides, and antibiotic resistance were discussed. The program along with the recordings of all sessions are freely available at https://www.difacquim.com/english/events/2022-colloquium/ .

DrugTax: package for drug taxonomy identification and explainable feature extraction

DrugTax is an easy-to-use Python package for small molecule detailed characterization. It extends a previously explored chemical taxonomy making it ready-to-use in any Artificial Intelligence approach. DrugTax leverages small molecule representations as input in one of their most accessible and simple forms (SMILES) and allows the simultaneously extraction of taxonomy information and key features for big data algorithm deployment. In addition, it delivers a set of tools for bulk analysis and visualization that can also be used for chemical space representation and molecule similarity assessment. DrugTax is a valuable tool for chemoinformatic processing and can be easily integrated in drug discovery pipelines. DrugTax can be effortlessly installed via PyPI (https://pypi.org/project/DrugTax/) or GitHub (https://github.com/MoreiraLAB/DrugTax).

Predicting the Potentiometric Sensitivity of Membrane Sensors Based on Modified Diphenylphosphoryl Acetamide Ionophores with QSPR Modeling

While potentiometric, plasticized membrane sensors are known as convenient, portable and inexpensive analytical instruments, their development is time- and resource-consuming, with a poorly predictable outcome. In this study, we investigated the applicability of the QSPR (quantitative structure-property relationship) method for predicting the potentiometric sensitivity of plasticized polymeric membrane sensors, using the ionophore chemical structure as model input. The QSPR model was based on the literature data on sensitivity, from previously studied, structurally similar ionophores, and it has shown reasonably good metrics in relating ionophore structures to their sensitivities towards Cu²⁺, Cd²⁺ and Pb²⁺. The model predictions for four newly synthesized diphenylphosphoryl acetamide ionophores were compared with real potentiometric experimental data for these ionophores, and satisfactory agreement was observed, implying the validity of the proposed approach.