Comprehensive strategies of machine-learning-based quantitative structure-activity relationship models

An assessment of crucial structural contributors of HDAC6 inhibitors through fragment-based non-linear pattern recognition and molecular dynamics simulation approaches

Amidst the Zn2+-dependant isoforms of the HDAC family, HDAC6 has emerged as a potential target associated with an array of diseases, especially cancer and neuronal disorders like Rett's Syndrome, Alzheimer's disease, Huntington's disease, etc. Also, despite the availability of a handful of HDAC inhibitors in the market, their non-selective nature has restricted their use in different disease conditions. In this situation, the development of selective and potent HDAC6 inhibitors will provide efficacious therapeutic agents to treat different diseases. In this context, this study has been carried out to evaluate the potential structural contributors of quinazoline-cap-containing HDAC6 inhibitors via machine learning (ML), conventional classification-dependant QSAR, and MD simulation-based binding mode of interaction analysis toward HDAC6 binding. This combined conventional and modern molecular modeling study has revealed the significance of the quinazoline moiety, substitutions present at the quinazoline cap group, as well as the importance of molecular property, number of hydrogen bond donor-acceptor functions, carbon-chlorine distance that significantly affects the HDAC6 binding of these inhibitors, subsequently affecting their potency . Interestingly, the study also revealed that the substitutions such as the chloroethyl group, and bulky quinazolinyl cap group can affect the binding of the cap function with the amino acid residues present in the loops proximal to the catalytic site of HDAC6. Such contributions of cap groups can lead to both stabilization and destabilization of the cap function after occupying the hydrophobic catalytic site by the aryl hydroxamate linker-ZBG functions.

Weak links: Advancing target-based drug discovery by identifying the most vulnerable targets

Mycobacterium tuberculosis remains the most common infectious killer worldwide despite decades of antitubercular drug development. Effectively controlling the tuberculosis (TB) pandemic will require innovation in drug discovery. In this review, we provide a brief overview of the two main approaches to discovering new TB drugs-phenotypic screens and target-based drug discovery-and outline some of the limitations of each method. We then explore recent advances in genetic tools that aim to overcome some of these limitations. In particular, we highlight a novel metric to prioritize essential targets, termed vulnerability. Stratifying targets based on their vulnerability presents new opportunities for future target-based drug discovery campaigns.

Advances in machine intelligence-driven virtual screening approaches for big-data

Virtual screening (VS) is an integral and ever-evolving domain of drug discovery framework. The VS is traditionally classified into ligand-based (LB) and structure-based (SB) approaches. Machine intelligence or artificial intelligence has wide applications in the drug discovery domain to reduce time and resource consumption. In combination with machine intelligence algorithms, VS has emerged into revolutionarily progressive technology that learns within robust decision orders for data curation and hit molecule screening from large VS libraries in minutes or hours. The exponential growth of chemical and biological data has evolved as "big-data" in the public domain demands modern and advanced machine intelligence-driven VS approaches to screen hit molecules from ultra-large VS libraries. VS has evolved from an individual approach (LB and SB) to integrated LB and SB techniques to explore various ligand and target protein aspects for the enhanced rate of appropriate hit molecule prediction. Current trends demand advanced and intelligent solutions to handle enormous data in drug discovery domain for screening and optimizing hits or lead with fewer or no false positive hits. Following the big-data drift and tremendous growth in computational architecture, we presented this review. Here, the article categorized and emphasized individual VS techniques, detailed literature presented for machine learning implementation, modern machine intelligence approaches, and limitations and deliberated the future prospects.

Decoding Nanomaterial-Biosystem Interactions through Machine Learning

The interactions between biosystems and nanomaterials regulate most of their theranostic and nanomedicine applications. These nanomaterial-biosystem interactions are highly complex and influenced by a number of entangled factors, including but not limited to the physicochemical features of nanomaterials, the types and characteristics of the interacting biosystems, and the properties of the surrounding microenvironments. Over the years, different experimental approaches coupled with computational modeling have revealed important insights into these interactions, although many outstanding questions remain unanswered. The emergence of machine learning has provided a timely and unique opportunity to revisit nanomaterial-biosystem interactions and to further push the boundary of this field. This minireview highlights the development and use of machine learning to decode nanomaterial-biosystem interactions and provides our perspectives on the current challenges and potential opportunities in this field.

Transformer-Based Molecular Generative Model for Antiviral Drug Design

Since the Simplified Molecular Input Line Entry System (SMILES) is oriented to the atomic-level representation of molecules and is not friendly in terms of human readability and editable, however, IUPAC is the closest to natural language and is very friendly in terms of human-oriented readability and performing molecular editing, we can manipulate IUPAC to generate corresponding new molecules and produce programming-friendly molecular forms of SMILES. In addition, antiviral drug design, especially analogue-based drug design, is also more appropriate to edit and design directly from the functional group level of IUPAC than from the atomic level of SMILES, since designing analogues involves altering the R group only, which is closer to the knowledge-based molecular design of a chemist. Herein, we present a novel data-driven self-supervised pretraining generative model called "TransAntivirus" to make select-and-replace edits and convert organic molecules into the desired properties for design of antiviral candidate analogues. The results indicated that TransAntivirus is significantly superior to the control models in terms of novelty, validity, uniqueness, and diversity. TransAntivirus showed excellent performance in the design and optimization of nucleoside and non-nucleoside analogues by chemical space analysis and property prediction analysis. Furthermore, to validate the applicability of TransAntivirus in the design of antiviral drugs, we conducted two case studies on the design of nucleoside analogues and non-nucleoside analogues and screened four candidate lead compounds against anticoronavirus disease (COVID-19). Finally, we recommend this framework for accelerating antiviral drug discovery.

Deciphering the Lexicon of Protein Targets: A Review on Multifaceted Drug Discovery in the Era of Artificial Intelligence

Understanding protein sequence and structure is essential for understanding protein-protein interactions (PPIs), which are essential for many biological processes and diseases. Targeting protein binding hot spots, which regulate signaling and growth, with rational drug design is promising. Rational drug design uses structural data and computational tools to study protein binding sites and protein interfaces to design inhibitors that can change these interactions, thereby potentially leading to therapeutic approaches. Artificial intelligence (AI), such as machine learning (ML) and deep learning (DL), has advanced drug discovery and design by providing computational resources and methods. Quantum chemistry is essential for drug reactivity, toxicology, drug screening, and quantitative structure-activity relationship (QSAR) properties. This review discusses the methodologies and challenges of identifying and characterizing hot spots and binding sites. It also explores the strategies and applications of artificial-intelligence-based rational drug design technologies that target proteins and protein-protein interaction (PPI) binding hot spots. It provides valuable insights for drug design with therapeutic implications. We have also demonstrated the pathological conditions of heat shock protein 27 (HSP27) and matrix metallopoproteinases (MMP2 and MMP9) and designed inhibitors of these proteins using the drug discovery paradigm in a case study on the discovery of drug molecules for cancer treatment. Additionally, the implications of benzothiazole derivatives for anticancer drug design and discovery are deliberated.

Insights into the structure-activity relationship of pyrimidine-sulfonamide analogues for targeting BRAF V600E protein

B-rapidly accelerated fibrosarcoma (BRAF) V600E plays a crucial role in the progression of cutaneous melanoma. Core structures of BRAF V600E inhibitors are based on pyrimidine-sulfonamide scaffolds. Exploring the QSAR of these structures can improve our understanding of BRAF V600E inhibitor drug design. This study utilized machine learning-based QSAR to elucidate chemical substructures of pyrimidine-sulfonamide analogues that correlated to the BRAF V600E inhibitory activity. The findings indicate that the support vector regression (SVR) combined with 15 fingerprints achieved the highest statistical performances in terms of goodness-of-fit, robustness, and predictability. Nine key fingerprints from pyrimidine-sulfonamide analogues were identified to exert the BRAF V600E inhibitory activity. These key fingerprints were validated using network-based activity cliff landscape and molecular docking. Together, the developed algorithm can serve as a screening tool for designing BRAF V600E inhibitors. To further utilize this model, we deployed our developed algorithm at https://qsarlabs.com/#braf.

Unlocking the Potential of High-Quality Dopamine Transporter Pharmacological Data: Advancing Robust Machine Learning-Based QSAR Modeling

The dopamine transporter (DAT) plays a critical role in the central nervous system and has been implicated in numerous psychiatric disorders. The ligand-based approaches are instrumental to decipher the structure-activity relationship (SAR) of DAT ligands, especially the quantitative SAR (QSAR) modeling. By gathering and analyzing data from literature and databases, we systematically assemble a diverse range of ligands binding to DAT, aiming to discern the general features of DAT ligands and uncover the chemical space for potential novel DAT ligand scaffolds. The aggregation of DAT pharmacological activity data, particularly from databases like ChEMBL, provides a foundation for constructing robust QSAR models. The compilation and meticulous filtering of these data, establishing high-quality training datasets with specific divisions of pharmacological assays and data types, along with the application of QSAR modeling, prove to be a promising strategy for navigating the pertinent chemical space. Through a systematic comparison of DAT QSAR models using training datasets from various ChEMBL releases, we underscore the positive impact of enhanced data set quality and increased data set size on the predictive power of DAT QSAR models.

Computational resources and chemoinformatics for translational health research

The integration of computational resources and chemoinformatics has revolutionized translational health research. It has offered a powerful set of tools for accelerating drug discovery. This chapter overviews the computational resources and chemoinformatics methods used in translational health research. The resources and methods can be used to analyze large datasets, identify potential drug candidates, predict drug-target interactions, and optimize treatment regimens. These resources have the potential to transform the drug discovery process and foster personalized medicine research. We discuss insights into their various applications in translational health and emphasize the need for addressing challenges, promoting collaboration, and advancing the field to fully realize the potential of these tools in transforming healthcare.

Harmonization Risks and Rewards: Nano-QSAR for Agricultural Nanomaterials

This comprehensive review explores the emerging landscape of Nano-QSAR (quantitative structure-activity relationship) for assessing the risk and potency of nanomaterials in agricultural settings. The paper begins with an introduction to Nano-QSAR, providing background and rationale, and explicitly states the hypotheses guiding the review. The study navigates through various dimensions of nanomaterial applications in agriculture, encompassing their diverse properties, types, and associated challenges. Delving into the principles of QSAR in nanotoxicology, this article elucidates its application in evaluating the safety of nanomaterials, while addressing the unique limitations posed by these materials. The narrative then transitions to the progression of Nano-QSAR in the context of agricultural nanomaterials, exemplified by insightful case studies that highlight both the strengths and the limitations inherent in this methodology. Emerging prospects and hurdles tied to Nano-QSAR in agriculture are rigorously examined, casting light on important pathways forward, existing constraints, and avenues for research enhancement. Culminating in a synthesis of key insights, the review underscores the significance of Nano-QSAR in shaping the future of nanoenabled agriculture. It provides strategic guidance to steer forthcoming research endeavors in this dynamic field.

From Molecular Descriptors to Intrinsic Fish Toxicity of Chemicals: An Alternative Approach to Chemical Prioritization

The European and U.S. chemical agencies have listed approximately 800k chemicals about which knowledge of potential risks to human health and the environment is lacking. Filling these data gaps experimentally is impossible, so in silico approaches and prediction are essential. Many existing models are however limited by assumptions (e.g., linearity and continuity) and small training sets. In this study, we present a supervised direct classification model that connects molecular descriptors to toxicity. Categories can be driven by either data (using k-means clustering) or defined by regulation. This was tested via 907 experimentally defined 96 h LC50 values for acute fish toxicity. Our classification model explained ≈90% of the variance in our data for the training set and ≈80% for the test set. This strategy gave a 5-fold decrease in the frequency of incorrect categorization compared to a quantitative structure-activity relationship (QSAR) regression model. Our model was subsequently employed to predict the toxicity categories of ≈32k chemicals. A comparison between the model-based applicability domain (AD) and the training set AD was performed, suggesting that the training set-based AD is a more adequate way to avoid extrapolation when using such models. The better performance of our direct classification model compared to that of QSAR methods makes this approach a viable tool for assessing the hazards and risks of chemicals.

DELFOS-drug efficacy leveraging forked and specialized networks-benchmarking scRNA-seq data in multi-omics-based prediction of cancer sensitivity

MOTIVATION

Cancer is currently one of the most notorious diseases, with over 1 million deaths in the European Union alone in 2022. As each tumor can be composed of diverse cell types with distinct genotypes, cancer cells can acquire resistance to different compounds. Moreover, anticancer drugs can display severe side effects, compromising patient well-being. Therefore, novel strategies for identifying the optimal set of compounds to treat each tumor have become an important research topic in recent decades.

RESULTS

To address this challenge, we developed a novel drug response prediction algorithm called Drug Efficacy Leveraging Forked and Specialized networks (DELFOS). Our model learns from multi-omics data from over 65 cancer cell lines, as well as structural data from over 200 compounds, for the prediction of drug sensitivity. We also evaluated the benefits of incorporating single-cell expression data to predict drug response. DELFOS was validated using datasets with unseen cell lines or drugs and compared with other state-of-the-art algorithms, achieving a high prediction performance on several correlation and error metrics. Overall, DELFOS can effectively leverage multi-omics data for the prediction of drug responses in thousands of drug-cell line pairs.

AVAILABILITY AND IMPLEMENTATION

The DELFOS pipeline and associated data are available at github.com/MoreiraLAB/delfos.

Development and validation of machine learning models for the prediction of SH-2 containing protein tyrosine phosphatase 2 inhibitors

Discovery and development of a new drug to the market is a highly challenging and resource consuming process. Although, modern drug discovery technologies have enabled the rapid identification of lead compounds, translation of the lead compounds into successful clinical candidates remains a big challenge. In recent years, the availability of massive structural and biological data of diverse small molecules and macromolecules has helped the researchers to deep mine the multidimensional data with the help of artificial intelligence-based predictive tools to draw useful insights on the structural features of biological or therapeutic significance. The aim of this study was to utilize the available data on small molecule (SH2)-containing protein tyrosine phosphatase 2 (SHP2) inhibitors to build and develop machine learning (ML) models that can predict the SHP2 inhibitory potential of new compounds. The dataset contained 2739 unique small molecule SHP2 inhibitors obtained from the BindingDB, ChEMBL and recent literature. After curation of the data, the predictive models such as XGBoost, K nearest neighbours, neural networks were developed and validated through a tenfold cross-validation testing procedure. Out of the seven models developed, the XGBoost model showed an excellent performance with ROC AUC score of 0.96 and accuracy of 0.97 on the test data. Moreover, the Shapley Additive Explanations method was applied to assess a more in-depth understanding of the influence of variables on the model's predictions. In summary, the XGBoost model developed in this study can be useful in the identification of novel SHP2 inhibitors and therefore, can accelerate the discovery of novel therapeutics for cancer therapy.

Unlocking the Potential: Predicting Redox Behavior of Organic Molecules, from Linear Fits to Neural Networks

Redox-active organic molecules, i.e., molecules that can relatively easily accept and/or donate electrons, are ubiquitous in biology, chemical synthesis, and electronic and spintronic devices, such as solar cells and rechargeable batteries, etc. Choosing the best candidates from an essentially infinite chemical space for experimental testing in a target application requires efficient screening approaches. In this Review, we discuss modern in silico techniques for predicting reduction and oxidation potentials of organic molecules that go beyond conventional first-principles computations and thermodynamic cycles. Approaches ranging from simple linear fits based on molecular orbital energy approximation and energy difference approximation to advanced regression and neural network machine learning algorithms employing complex descriptors of molecular compositions, geometries, and electronic structures are examined in conjunction with relevant literature examples. We discuss the interplay between ab initio data and machine learning (ML), i.e., whether it is better to base predictions on low-level quantum-chemical results corrected with ML or to bypass first-principles computations entirely and instead rely on elaborate deep learning architectures. Finally, we list currently available data sets of redox-active organic molecules and their experimental and/or computed properties to facilitate the development of screening platforms and rational design of redox-active organic molecules.

Prediction of KRASG12C inhibitors using conjoint fingerprint and machine learning-based QSAR models

Kirsten rat sarcoma virus G12C (KRAS^G12C) is the major protein mutation associated with non-small cell lung cancer (NSCLC) severity. Inhibiting KRAS^G12C is therefore one of the key therapeutic strategies for NSCLC patients. In this paper, a cost-effective data driven drug design employing machine learning-based quantitative structure-activity relationship (QSAR) analysis was built for predicting ligand affinities against KRAS^G12C protein. A curated and non-redundant dataset of 1033 compounds with KRAS^G12C inhibitory activity (pIC₅₀) was used to build and test the models. The PubChem fingerprint, Substructure fingerprint, Substructure fingerprint count, and the conjoint fingerprint-a combination of PubChem fingerprint and Substructure fingerprint count-were used to train the models. Using comprehensive validation methods and various machine learning algorithms, the results clearly showed that the XGBoost regression (XGBoost) achieved the highest performance in term of goodness of fit, predictivity, generalizability and model robustness (R² = 0.81, Q²_CV = 0.60, Q²_Ext = 0.62, R² - Q²_Ext = 0.19, R²_Y-Random = 0.31 ± 0.03, Q²_Y-Random = -0.09 ± 0.04). The top 13 molecular fingerprints that correlated with the predicted pIC₅₀ values were SubFPC274 (aromatic atoms), SubFPC307 (number of chiral-centers), PubChemFP37 (≥1 Chlorine), SubFPC18 (Number of alkylarylethers), SubFPC1 (number of primary carbons), SubFPC300 (number of 1,3-tautomerizables), PubChemFP621 (N-C:C:C:N structure), PubChemFP23 (≥1 Fluorine), SubFPC2 (number of secondary carbons), SubFPC295 (number of C-ONS bonds), PubChemFP199 (≥4 6-membered rings), PubChemFP180 (≥1 nitrogen-containing 6-membered ring), and SubFPC180 (number of tertiary amine). These molecular fingerprints were virtualized and validated using molecular docking experiments. In conclusion, this conjoint fingerprint and XGBoost-QSAR model demonstrated to be useful as a high-throughput screening tool for KRAS^G12C inhibitor identification and drug design.

Design and Identification of a Novel Antiviral Affinity Peptide against Fowl Adenovirus Serotype 4 (FAdV-4) by Targeting Fiber2 Protein

Outbreaks of hydropericardium hepatitis syndrome caused by fowl adenovirus serotype 4 (FAdV-4) with a novel genotype have been reported in China since 2015, with significant economic losses to the poultry industry. Fiber2 is one of the important structural proteins on FAdV-4 virions. In this study, the C-terminal knob domain of the FAdV-4 Fiber2 protein was expressed and purified, and its trimer structure (PDB ID: 7W83) was determined for the first time. A series of affinity peptides targeting the knob domain of the Fiber2 protein were designed and synthesized on the basis of the crystal structure using computer virtual screening technology. A total of eight peptides were screened using an immunoperoxidase monolayer assay and RT-qPCR, and they exhibited strong binding affinities to the knob domain of the FAdV-4 Fiber2 protein in a surface plasmon resonance assay. Treatment with peptide number 15 (P15; WWHEKE) at different concentrations (10, 25, and 50 μM) significantly reduced the expression level of the Fiber2 protein and the viral titer during FAdV-4 infection. P15 was found to be an optimal peptide with antiviral activity against FAdV-4 in vitro with no cytotoxic effect on LMH cells up to 200 μM. This study led to the identification of a class of affinity peptides designed using computer virtual screening technology that targeted the knob domain of the FAdV-4 Fiber2 protein and may be developed as a novel potential and effective antiviral strategy in the prevention and control of FAdV-4.

Re-envisioning the design of nanomedicines: harnessing automation and artificial intelligence

INTRODUCTION

Interest in nanomedicines has surged in recent years due to the critical role they have played in the COVID-19 pandemic. Nanoformulations can turn promising therapeutic cargo into viable products through improvements in drug safety and efficacy profiles. However, the developmental pathway for such formulations is non-trivial and largely reliant on trial-and-error. Beyond the costly demands on time and resources, this traditional approach may stunt innovation. The emergence of automation, artificial intelligence (AI) and machine learning (ML) tools, which are currently underutilized in pharmaceutical formulation development, offers a promising direction for an improved path in the design of nanomedicines.

AREAS COVERED

the potential of harnessing experimental automation and AI/ML to drive innovation in nanomedicine development. The discussion centers on the current challenges in drug formulation research and development, and the major advantages afforded through the application of data-driven methods.

EXPERT OPINION

The development of integrated workflows based on automated experimentation and AI/ML may accelerate nanomedicine development. A crucial step in achieving this is the generation of high-quality, accessible datasets. Future efforts to make full use of these tools can ultimately contribute to the development of more innovative nanomedicines and improved clinical translation of formulations that rely on advanced drug delivery systems.

Guidance for good practice in the application of machine learning in development of toxicological quantitative structure-activity relationships (QSARs)

Recent years have seen a substantial growth in the adoption of machine learning approaches for the purposes of quantitative structure-activity relationship (QSAR) development. Such a trend has coincided with desire to see a shifting in the focus of methodology employed within chemical safety assessment: away from traditional reliance upon animal-intensive in vivo protocols, and towards increased application of in silico (or computational) predictive toxicology. With QSAR central amongst techniques applied in this area, the emergence of algorithms trained through machine learning with the objective of toxicity estimation has, quite naturally, arisen. On account of the pattern-recognition capabilities of the underlying methods, the statistical power of the ensuing models is potentially considerable-appropriate for the handling even of vast, heterogeneous datasets. However, such potency comes at a price: this manifesting as the general practical deficits observed with respect to the reproducibility, interpretability and generalisability of the resulting tools. Unsurprisingly, these elements have served to hinder broader uptake (most notably within a regulatory setting). Areas of uncertainty liable to accompany (and hence detract from applicability of) toxicological QSAR have previously been highlighted, accompanied by the forwarding of suggestions for "best practice" aimed at mitigation of their influence. However, the scope of such exercises has remained limited to "classical" QSAR-that conducted through use of linear regression and related techniques, with the adoption of comparatively few features or descriptors. Accordingly, the intention of this study has been to extend the remit of best practice guidance, so as to address concerns specific to employment of machine learning within the field. In doing so, the impact of strategies aimed at enhancing the transparency (feature importance, feature reduction), generalisability (cross-validation) and predictive power (hyperparameter optimisation) of algorithms, trained upon real toxicity data through six common learning approaches, is evaluated.

Application of a deep generative model produces novel and diverse functional peptides against microbial resistance

Antimicrobial resistance could threaten millions of lives in the immediate future. Antimicrobial peptides (AMPs) are an alternative to conventional antibiotics practice against infectious diseases. Despite the potential contribution of AMPs to the antibiotic's world, their development and optimization have encountered serious challenges. Cutting-edge methods with novel and improved selectivity toward resistant targets must be established to create AMPs-driven treatments. Here, we present AMPTrans-lstm, a deep generative network-based approach for the rational design of AMPs. The AMPTrans-lstm pipeline involves pre-training, transfer learning, and module identification. The AMPTrans-lstm model has two sub-models, namely, (long short-term memory) LSTM sampler and Transformer converter, which can be connected in series to make full use of the stability of LSTM and the novelty of Transformer model. These elements could generate AMPs candidates, which can then be tailored for specific applications. By analyzing the generated sequence and trained AMPs, we prove that AMPTrans-lstm can expand the design space of the trained AMPs and produce reasonable and brand-new AMPs sequences. AMPTrans-lstm can generate functional peptides for antimicrobial resistance with good novelty and diversity, so it is an efficient AMPs design tool.

MolPredictX: Online Biological Activity Predictions by Machine Learning Models

Here we report the development of MolPredictX, an innovate and freely accessible web interface for biological activity predictions of query molecules. MolPredictX utilizes in-house QSAR models to provide 27 qualitative predictions (active or inactive), and quantitative probabilities for bioactivity against parasitic (Trypanosoma and Leishmania), viral (Dengue, Sars-CoV and Hepatitis C), pathogenic yeast (Candida albicans), bacterial (Salmonella enterica and Escherichia coli), and Alzheimer disease enzymes. In this article, we introduce the methodology and usability of this webtool, highlighting its potential role in the development of new drugs against a variety of diseases. MolPredictX is undergoing continuous development and is freely available at https://www.molpredictx.ufpb.br/.

Enabling data-limited chemical bioactivity predictions through deep neural network transfer learning

The main limitation in developing deep neural network (DNN) models to predict bioactivity properties of chemicals is the lack of sufficient assay data to train the network's classification layers. Focusing on feedforward DNNs that use atom- and bond-based structural fingerprints as input, we examined whether layers of a fully trained DNN based on large amounts of data to predict one property could be used to develop DNNs to predict other related or unrelated properties based on limited amounts of data. Hence, we assessed if and under what conditions the dense layers of a pre-trained DNN could be transferred and used for the development of another DNN associated with limited training data. We carried out a quantitative study employing more than 400 pairs of assay datasets, where we used fully trained layers from a large dataset to augment the training of a small dataset. We found that the higher the correlation r between two assay datasets, the more efficient the transfer learning is in reducing prediction errors associated with the smaller dataset DNN predictions. The reduction in mean squared prediction errors ranged from 10 to 20% for every 0.1 increase in r² between the datasets, with the bulk of the error reductions associated with transfers of the first dense layer. Transfer of other dense layers did not result in additional benefits, suggesting that deeper, dense layers conveyed more specialized and assay-specific information. Importantly, depending on the dataset correlation, training sample size could be reduced by up to tenfold without any loss of prediction accuracy.

Machine learning and structure-based modeling for the prediction of UDP-glucuronosyltransferase inhibition

UDP-glucuronosyltransferases (UGTs) are responsible for 35% of the phase II drug metabolism. In this study, we focused on UGT1A1, which is a key UGT isoform. Strong inhibition of UGT1A1 may trigger adverse drug/herb-drug interactions, or result in disorders of endobiotic metabolism. Most of the current machine learning methods predicting the inhibition of drug metabolizing enzymes neglect protein structure and dynamics, both being essential for the recognition of various substrates and inhibitors. We performed molecular dynamics simulations on a homology model of the human UGT1A1 structure containing both the cofactor- (UDP-glucuronic acid) and substrate-binding domains to explore UGT conformational changes. Then, we created models for the prediction of UGT1A1 inhibitors by integrating information on UGT1A1 structure and dynamics, interactions with diverse ligands, and machine learning. These models can be helpful for further prediction of drug-drug interactions of drug candidates and safety treatments.

•

UGTs are responsible for 35% of the phase II drug metabolism reactions

•

We created machine learning models for prediction of UGT1A1 inhibitors

•

Our simulations suggested key residues of UGT1A1 involved in the substrate binding

Quantitative Structure-Toxicity Relationship in Bioactive Molecules from a Conceptual DFT Perspective

The preclinical drug discovery stage often requires a large amount of costly and time-consuming experiments using huge sets of chemical compounds. In the last few decades, this process has undergone significant improvements by the introduction of quantitative structure-activity relationship (QSAR) modelling that uses a certain percentage of experimental data to predict the biological activity/property of compounds with similar structural skeleton and/or containing a particular functional group(s). The use of machine learning tools along with it has made life even easier for pharmaceutical researchers. Here, we discuss the toxicity of certain sets of bioactive compounds towards Pimephales promelas and Tetrahymena pyriformis in terms of the global conceptual density functional theory (CDFT)-based descriptor, electrophilicity index (ω). We have compared the results with those obtained by using the commonly used hydrophobicity parameter, logP (where P is the n-octanol/water partition coefficient), considering the greater ease of computing the ω descriptor. The Human African trypanosomiasis (HAT) curing activity of 32 pyridyl benzamide derivatives is also studied against Tryphanosoma brucei. In this review article, we summarize these multiple linear regression (MLR)-based QSAR studies in terms of electrophilicity (ω, ω2) and hydrophobicity (logP, (logP)2) parameters.

Computational pharmacology: New avenues for COVID-19 therapeutics search and better preparedness for future pandemic crises

The COVID-19 pandemic created an unprecedented global healthcare emergency prompting the exploration of new therapeutic avenues, including drug repurposing. A large number of ongoing studies revealed pervasive issues in clinical research, such as the lack of accessible and organised data. Moreover, current shortcomings in clinical studies highlighted the need for a multi-faceted approach to tackle this health crisis. Thus, we set out to explore and develop new strategies for drug repositioning by employing computational pharmacology, data mining, systems biology, and computational chemistry to advance shared efforts in identifying key targets, affected networks, and potential pharmaceutical intervention options. Our study revealed that formulating pharmacological strategies should rely on both therapeutic targets and their networks. We showed how data mining can reveal regulatory patterns, capture novel targets, alert about side-effects, and help identify new therapeutic avenues. We also highlighted the importance of the miRNA regulatory layer and how this information could be used to monitor disease progression or devise treatment strategies. Importantly, our work bridged the interactome with the chemical compound space to better understand the complex landscape of COVID-19 drugs. Machine and deep learning allowed us to showcase limitations in current chemical libraries for COVID-19 suggesting that both in silico and experimental analyses should be combined to retrieve therapeutically valuable compounds. Based on the gathered data, we strongly advocate for taking this opportunity to establish robust practices for treating today's and future infectious diseases by preparing solid analytical frameworks.

Predicting bioconcentration factor and estrogen receptor bioactivity of bisphenol a and its analogues in adult zebrafish by directed message passing neural networks

The bioconcentration factor (BCF) is a key parameter for bioavailability assessment of environmental pollutants in regulatory frameworks. The comparative toxicology and mechanism of action of congeners are also of concern. However, there are limitations to acquire them by conducting field and laboratory experiments while machinelearning is emerging as a promising predictive tool to fill the gap. In this study, the Direct Message Passing Neural Network (DMPNN) was applied to predict logBCFs of bisphenol A (BPA) and its four analogues (bisphenol AF (BPAF), bisphenol B (BPB), bisphenol F (BPF) and bisphenol S (BPS)). For the test set, the Pearson correlation coefficient (PCC) and mean square error (MSE) were 0.85 and 0.52 respectively, suggesting a good predictive performance. The predicted logBCFs values by the DMPNN ranging from 0.35 (BPS) to 2.14 (BPAF) coincided well with those by the classical EPI Suite (BCFBAF model). Besides, estrogen receptor α (ERα) bioactivity of these bisphenols was also predicted well by the DMPNN, with a probability of 97.0 % (BPB) to 99.7 % (BPAF), which was validated by the extent of vitellogenin (VTG) induction in male zebrafish as a biomarker except BPS. Thus, with little need for expert knowledge, DMPNN is confirmed to be a useful tool to accurately predict logBCF and screen for estrogenic activity from molecular structures. Moreover, a gender difference was noted in the changes of three endpoints (logBCF, ER binding affinity and VTG levels), the rank order of which was BPAF > BPB > BPA > BPF > BPS consistently, and abnormal amino acid metabolism is featured as an omics signature of abnormal hormone protein expression.

hERG Blockade Prediction by Combining Site Identification by Ligand Competitive Saturation and Physicochemical Properties

Human ether-a-go-go-related gene (hERG) potassium channel is well-known contributor to drug-induced cardiotoxicity and therefore an extremely important target when performing safety assessments of drug candidates. Ligand-based approaches in connection with quantitative structure active relationships (QSAR) analyses have been developed to predict hERG toxicity. Availability of the recent published cryogenic electron microscopy (cryo-EM) structure for the hERG channel opened the prospect for using structure-based simulation and docking approaches for hERG drug liability predictions. In recent time, the idea of combining structure- and ligand-based approaches for modeling hERG drug liability has gained momentum offering improvements in predictability when compared to ligand-based QSAR practices alone. The present article demonstrates uniting the structure-based SILCS (site-identification by ligand competitive saturation) approach in conjunction with physicochemical properties to develop predictive models for hERG blockade. This combination leads to improved model predictability based on Pearson's R and percent correct (represents rank-ordering of ligands) metric for different validation sets of hERG blockers involving diverse chemical scaffold and wide range of pIC50 values. The inclusion of the SILCS structure-based approach allows determination of the hERG region to which compounds bind and the contribution of different chemical moieties in the compounds to blockade, thereby facilitating the rational ligand design to minimize hERG liability.

Development of an innovative data-driven system to generate descriptive prediction equation of dielectric constant on small sample sets

Dielectric constant (DC, ε) is a fundamental parameter in material sciences to measure polarizability of the system. In industrial processes, its value is an imperative indicator, which demonstrates the dielectric property of material and compiles information including separation information, chemical equilibrium, chemical reactivity analysis, and solubility modeling. Since, the available ε-prediction models are fairly primitive and frequently suffer from serious failures especially when deals with strong polar compounds. Therefore, we have developed a novel data-driven system to improve the efficiency and wide-range applicability of ε using in material sciences. This innovative scheme adopts the correlation distance and genetic algorithm to discriminate features’ combination and avoid overfitting. Herein, the prediction output of the single ML model as a coding to estimate the target value by simulating the layer-by-layer extraction in deep learning, and enabling instant search for the optimal combination of features is recruited. Our model established an improved correlation value of 0.956 with target as compared to the previously available best traditional ML result of 0.877. Our framework established a profound improvement, especially for material systems possessing ε value >50. In terms of interpretability, we have derived a conceptual computational equation from a minimum generating tree. Our innovative data-driven system is preferentially superior over other methods due to its application for the prediction of dielectric constants as well as for the prediction of overall micro and macro-properties of any multi-components complex.

De novo molecular design with deep molecular generative models for PPI inhibitors

We construct a protein-protein interaction (PPI) targeted drug-likeness dataset and propose a deep molecular generative framework to generate novel drug-likeness molecules from the features of the seed compounds. This framework gains inspiration from published molecular generative models, uses the key features associated with PPI inhibitors as input and develops deep molecular generative models for de novo molecular design of PPI inhibitors. For the first time, quantitative estimation index for compounds targeting PPI was applied to the evaluation of the molecular generation model for de novo design of PPI-targeted compounds. Our results estimated that the generated molecules had better PPI-targeted drug-likeness and drug-likeness. Additionally, our model also exhibits comparable performance to other several state-of-the-art molecule generation models. The generated molecules share chemical space with iPPI-DB inhibitors as demonstrated by chemical space analysis. The peptide characterization-oriented design of PPI inhibitors and the ligand-based design of PPI inhibitors are explored. Finally, we recommend that this framework will be an important step forward for the de novo design of PPI-targeted therapeutics.

Natural Products from Annonaceae as Potential Antichagasic Agents

Chagas disease, a neglected tropical disease, is endemic in 21 Latin American countries and particularly prevalent in Brazil. Chagas disease has drawn more attention in recent years due to its expansion into non-endemic areas. The aim of this work was to computationally identify and experimentally validate the natural products from an Annonaceae family as antichagasic agents. Through the ligand-based virtual screening, we identified 57 molecules with potential activity against the epimastigote form of T. cruzi. Then, 16 molecules were analyzed in the in vitro study, of which, six molecules displayed previously unknown antiepimastigote activity. We also evaluated these six molecules for trypanocidal activity. We observed that all six molecules have potential activity against the amastigote form, but no molecules were active against the trypomastigote form. 13-Epicupressic acid seems to be the most promising, as it was predicted as an active compound in the in silico study against the amastigote form of T. cruzi, in addition to having in vitro activity against the epimastigote form.

Investigation of potential descriptors of chemical compounds on prevention of nephrotoxicity via QSAR approach

Drug-induced nephrotoxicity remains a common problem after exposure to medications and diagnostic agents, which may be heightened in the kidney microenvironment and deteriorate kidney function. In this study, the toxic effects of fourteen marked drugs with the individual chemical structure were evaluated in kidney cells. The quantitative structure–activity relationship (QSAR) approach was employed to investigate the potential structural descriptors of each drug-related to their toxic effects. The most reasonable equation of the QSAR model displayed that the estimated regression coefficients such as the number of ring assemblies, three-membered rings, and six-membered rings were strongly related to toxic effects on renal cells. Meanwhile, the chemical properties of the tested compounds including carbon atoms, bridge bonds, H-bond donors, negative atoms, and rotatable bonds were favored properties and promote the toxic effects on renal cells. Particularly, more numbers of rotatable bonds were positively correlated with strong toxic effects that displayed on the most toxic compound. The useful information discovered from our regression QSAR models may help to identify potential hazardous moiety to avoid nephrotoxicity in renal preventive medicine.

Machine Learning Analysis of Essential Oils from Cuban Plants: Potential Activity against Protozoa Parasites

Essential oils (EOs) are a mixture of chemical compounds with a long history of use in food, cosmetics, perfumes, agricultural and pharmaceuticals industries. The main object of this study was to find chemical patterns between 45 EOs and antiprotozoal activity (antiplasmodial, antileishmanial and antitrypanosomal), using different machine learning algorithms. In the analyses, 45 samples of EOs were included, using unsupervised Self-Organizing Maps (SOM) and supervised Random Forest (RF) methodologies. In the generated map, the hit rate was higher than 70% and the results demonstrate that it is possible find chemical patterns using a supervised and unsupervised machine learning approach. A total of 20 compounds were identified (19 are terpenes and one sulfur-containing compound), which was compared with literature reports. These models can be used to investigate and screen for bioactivity of EOs that have antiprotozoal activity more effectively and with less time and financial cost.

A Deep Learning-Based Quantitative Structure-Activity Relationship System Construct Prediction Model of Agonist and Antagonist with High Performance

Molecular design and evaluation for drug development and chemical safety assessment have been advanced by quantitative structure–activity relationship (QSAR) using artificial intelligence techniques, such as deep learning (DL). Previously, we have reported the high performance of prediction models molecular initiation events (MIEs) on the adverse toxicological outcome using a DL-based QSAR method, called DeepSnap-DL. This method can extract feature values from images generated on a three-dimensional (3D)-chemical structure as a novel QSAR analytical system. However, there is room for improvement of this system’s time-consumption. Therefore, in this study, we constructed an improved DeepSnap-DL system by combining the processes of generating an image from a 3D-chemical structure, DL using the image as input data, and statistical calculation of prediction-performance. Consequently, we obtained that the three prediction models of agonists or antagonists of MIEs achieved high prediction-performance by optimizing the parameters of DeepSnap, such as the angle used in the depiction of the image of a 3D-chemical structure, data-split, and hyperparameters in DL. The improved DeepSnap-DL system will be a powerful tool for computer-aided molecular design as a novel QSAR system.

Why 90% of clinical drug development fails and how to improve it?

Ninety percent of clinical drug development fails despite implementation of many successful strategies, which raised the question whether certain aspects in target validation and drug optimization are overlooked? Current drug optimization overly emphasizes potency/specificity using structure‒activity-relationship (SAR) but overlooks tissue exposure/selectivity in disease/normal tissues using structure‒tissue exposure/selectivity–relationship (STR), which may mislead the drug candidate selection and impact the balance of clinical dose/efficacy/toxicity. We propose structure‒tissue exposure/selectivity–activity relationship (STAR) to improve drug optimization, which classifies drug candidates based on drug's potency/selectivity, tissue exposure/selectivity, and required dose for balancing clinical efficacy/toxicity. Class I drugs have high specificity/potency and high tissue exposure/selectivity, which needs low dose to achieve superior clinical efficacy/safety with high success rate. Class II drugs have high specificity/potency and low tissue exposure/selectivity, which requires high dose to achieve clinical efficacy with high toxicity and needs to be cautiously evaluated. Class III drugs have relatively low (adequate) specificity/potency but high tissue exposure/selectivity, which requires low dose to achieve clinical efficacy with manageable toxicity but are often overlooked. Class IV drugs have low specificity/potency and low tissue exposure/selectivity, which achieves inadequate efficacy/safety, and should be terminated early. STAR may improve drug optimization and clinical studies for the success of clinical drug development.

ER/AR Multi-Conformational Docking Server: A Tool for Discovering and Studying Estrogen and Androgen Receptor Modulators

The prediction of the estrogen receptor (ER) and androgen receptor (AR) activity of a compound is quite important to avoid the environmental exposures of endocrine-disrupting chemicals. The Estrogen and Androgen Receptor Database (EARDB, http://eardb.schanglab.org.cn/) provides a unique collection of reported ERα, ERβ, or AR protein structures and known small molecule modulators. With the user-uploaded query molecules, molecular docking based on multi-conformations of a single target will be performed. Moreover, the 2D similarity search against known modulators is also provided. Molecules predicted with a low binding energy or high similarity to known ERα, ERβ, or AR modulators may be potential endocrine-disrupting chemicals or new modulators. The server provides a tool to predict the endocrine activity for compounds of interests, benefiting for the ER and AR drug design and endocrine-disrupting chemical identification.