In Vitro Characterization and Rescue of VX Metabolism in Human Liver Microsomes.

Venomous agent X (VX) is an organophosphate acetylcholinesterase (AChE) inhibitor, and although it is one of the most toxic AChE inhibitors known, the extent of metabolism in humans is not currently well understood. The known metabolism in humans is limited to the metabolite identification from a single victim of the Osaka poisoning in 1994, which allowed for the identification of several metabolic products. VX has been reported to be metabolized in vitro by paraoxonase-1 and phosphotriesterase, although their binding constants are many orders of magnitude above the LD50, suggesting limited physiologic relevance. Using incubation with human liver microsomes (HLMs), we have now characterized the metabolism of VX and the formation of multiple metabolites as well as identified a Food and Drug Administration-approved drug [ethylenediaminetetraacetic acid (EDTA)] that enhances the metabolic rate. HLM incubation alone shows a pronounced increase in the metabolism of VX compared with buffer, suggesting that cytochrome P450-mediated metabolism of VX is occurring. We identified a biphasic decay with two distinct rates of metabolism. The enhancement of VX metabolism in multiple buffers was assessed to attempt to mitigate the effect of hydrolysis rates. The formation of VX metabolites was shown to be shifted with HLMs, suggesting a pathway enhancement over simple hydrolysis. Additionally, our investigation of hydrolysis rates in various common buffers used in biologic assays discovered dramatic differences in VX stability. The new human in vitro VX metabolic data reported points to a potential in vivo treatment strategy (EDTA) for rescue in individuals that are poisoned though enhancement of metabolism alongside existing treatments. SIGNIFICANCE STATEMENT: Venomous agent X (VX) is a potent acetylcholinesterase inhibitor and chemical weapon. To date, we do not possess a clear understanding of its metabolism in humans that would assist us in treating those exposed to it. This study now describes the human liver microsomal metabolism of VX and identifies ethylenediaminetetraacetic acid, which appears to enhance the rate of metabolism. This may provide a potential treatment option for human VX poisoning.

Sequential Contrastive and Deep Learning Models to Identify Selective Butyrylcholinesterase Inhibitors.

Butyrylcholinesterase (BChE) is a target of interest in late-stage Alzheimer's Disease (AD) where selective BChE inhibitors (BIs) may offer symptomatic treatment without the harsh side effects of acetylcholinesterase (AChE) inhibitors. In this study, we explore multiple machine learning strategies to identify BIs in silico, optimizing for precision over all other metrics. We compare state-of-the-art supervised contrastive learning (CL) with deep learning (DL) and Random Forest (RF) machine learning, across single and sequential modeling configurations, to identify the best models for BChE selectivity. We used these models to virtually screen a vendor library of 5 million compounds for BIs and tested 20 of these compounds in vitro. Seven of the 20 compounds displayed selectivity for BChE over AChE, reflecting a hit rate of 35% for our model predictions, suggesting a highly efficient strategy for modeling selective inhibition.

Human CYP2C19 Substrate and Inhibitor Characterization of Organophosphate Pesticides.

CYP2C19 is an important enzyme for organophosphate pesticide (OPP) metabolism. Because the OPPs can be both substrates and inhibitors of CYP2C19, we screened 45 OPPs for their ability to inhibit the activity of this enzyme and investigated the role of CYP2C19 in the metabolism of 22 of these molecules. We identified several nanomolar inhibitors of CYP2C19 as well as determined that thions, in general, are more potent inhibitors than oxons. We also determined that thions are readily metabolized by CYP2C19, although we saw no relationship between IC50 values and intrinsic clearance rates. This study may have implications for mitigating the risk of OPP poisoning.

Validation of Acetylcholinesterase Inhibition Machine Learning Models for Multiple Species.

Acetylcholinesterase (AChE) is an important enzyme and target for human therapeutics, environmental safety, and global food supply. Inhibitors of this enzyme are also used for pest elimination and can be misused for suicide or chemical warfare. Adverse effects of AChE pesticides on nontarget organisms, such as fish, amphibians, and humans, have also occurred as a result of biomagnifications of these toxic compounds. We have exhaustively curated the public data for AChE inhibition data and developed machine learning classification models for seven different species. Each set of models were built using up to nine different algorithms for each species and Morgan fingerprints (ECFP6) with an activity cutoff of 1 µM. The human (4075 compounds) and eel (5459 compounds) consensus models predicted AChE inhibition activity using external test sets from literature data with 81% and 82% accuracy, respectively, while the reciprocal cross (76% and 82% percent accuracy) was not species-specific. In addition, we also created machine learning regression models for human and eel AChE inhibition to return a predicted IC50 value for a queried molecule. We did observe an improved species specificity in the regression models, where a human support vector regression model of human AChE inhibition (3652 compounds) predicted the IC50s of the human test set to a better extent than the eel regression model (4930 compounds) on the same test set, based on mean absolute percentage error (MAPE = 9.73% vs 13.4%). The predictive power of these models certainly benefits from increasing the chemical diversity of the training set, as evidenced by expanding our human classification model by incorporating data from the Tox21 library of compounds. Of the 10 compounds we tested that were predicted active by this expanded model, two showed >80% inhibition at 100 µM. This machine learning approach therefore offers the ability to rapidly score massive libraries of molecules against the models for AChE inhibition that can then be selected for future in vitro testing to identify potential toxins. It also enabled us to create a public website, MegaAChE, for single-molecule predictions of AChE inhibition using these models at megaache.collaborationspharma.com.

In silico ADME/tox comes of age: twenty years later.

In the early 2000s pharmaceutical drug discovery was beginning to use computational approaches for absorption, distribution, metabolism, excretion and toxicity (ADME/Tox, also known as ADMET) prediction. This emphasis on prediction was an effort to reduce the risk of later stage failures from ADME/Tox.Much has been written in the intervening twenty plus years and significant expenditure has occurred in companies developing these in silico capabilities which can be gleaned from publications. It is therefore an appropriate time to briefly reflect on what was proposed then and what the reality is today.20 years ago, we tended to optimise bioactivity and perhaps one ADME/Tox property at a time. Previously pharmaceutical companies needed a whole infrastructure for models - in silico and in vitro experts, IT, champions on a project team, educators and management support. Now we are in the age of generative de novo design where bioactivity and many ADME/Tox properties can be optimised and large language model technologies are available.There are also some challenges such as the focus on very large molecules which may be outside of current ADME/Tox models.We provide an opportunity to look forward with the increasing public data for ADME/Tox as well as expanded types of algorithms available.

Machine Learning Models Identify New Inhibitors for Human OATP1B1.

The uptake transporter OATP1B1 (SLC01B1) is largely localized to the sinusoidal membrane of hepatocytes and is a known victim of unwanted drug-drug interactions. Computational models are useful for identifying potential substrates and/or inhibitors of clinically relevant transporters. Our goal was to generate OATP1B1 in vitro inhibition data for [3H] estrone-3-sulfate (E3S) transport in CHO cells and use it to build machine learning models to facilitate a comparison of seven different classification models (Deep learning, Adaboosted decision trees, Bernoulli naïve bayes, k-nearest neighbors (knn), random forest, support vector classifier (SVC), logistic regression (lreg), and XGBoost (xgb)] using ECFP6 fingerprints to perform 5-fold, nested cross validation. In addition, we compared models using 3D pharmacophores, simple chemical descriptors alone or plus ECFP6, as well as ECFP4 and ECFP8 fingerprints. Several machine learning algorithms (SVC, lreg, xgb, and knn) had excellent nested cross validation statistics, particularly for accuracy, AUC, and specificity. An external test set containing 207 unique compounds not in the training set demonstrated that at every threshold SVC outperformed the other algorithms based on a rank normalized score. A prospective validation test set was chosen using prediction scores from the SVC models with ECFP fingerprints and were tested in vitro with 15 of 19 compounds (84% accuracy) predicted as active (≥20% inhibition) showed inhibition. Of these compounds, six (abamectin, asiaticoside, berbamine, doramectin, mobocertinib, and umbralisib) appear to be novel inhibitors of OATP1B1 not previously reported. These validated machine learning models can now be used to make predictions for drug-drug interactions for human OATP1B1 alongside other machine learning models for important drug transporters in our MegaTrans software.

Machine Learning Models for Mycobacterium tuberculosisIn Vitro Activity: Prediction and Target Visualization.

Tuberculosis (TB) is a major global health challenge, with approximately 1.4 million deaths per year. There is still a need to develop novel treatments for patients infected with Mycobacterium tuberculosis (Mtb). There have been many large-scale phenotypic screens that have led to the identification of thousands of new compounds. Yet, there is very limited investment in TB drug discovery which points to the need for new methods to increase the efficiency of drug discovery against Mtb. We have used machine learning approaches to learn from the public Mtb data, resulting in many data sets and models with robust enrichment and hit rates leading to the discovery of new active compounds. Recently, we have curated predominantly small-molecule Mtb data and developed new machine learning classification models with 18 886 molecules at different activity cutoffs. We now describe the further validation of these Bayesian models using a library of over 1000 molecules synthesized as part of EU-funded New Medicines for TB and More Medicines for TB programs. We highlight molecular features which are enriched in these active compounds. In addition, we provide new regression and classification models that can be used for scoring compound libraries or used to design new molecules. We have also visualized these molecules in the context of known molecular targets and identified clusters in chemical property space, which may aid in future target identification efforts. Finally, we are also making these data sets publicly available, representing a significant increase to the available Mtb inhibition data in the public domain.

Multiple Computational Approaches for Predicting Drug Interactions with Human Equilibrative Nucleoside Transporter 1.

Equilibrativenucleoside transporters (ENTs) participate in the pharmacokinetics and disposition of nucleoside analog drugs. Understanding drug interactions with the ENTs may inform and facilitate the development of new drugs, including chemotherapeutics and antivirals that require access to sanctuary sites such as the male genital tract. This study created three-dimensional pharmacophores for ENT1 and ENT2 substrates and inhibitors using Kt and IC50 data curated from the literature. Substrate pharmacophores for ENT1 and ENT2 are distinct, with partial overlap of hydrogen bond donors, whereas the inhibitor pharmacophores predominantly feature hydrogen bond acceptors. Mizoribine and ribavirin mapped to the ENT1 substrate pharmacophore and proved to be substrates of the ENTs. The presence of the ENT-specific inhibitor 6-S-[(4-nitrophenyl)methyl]-6-thioinosine (NBMPR) decreased mizoribine accumulation in ENT1 and ENT2 cells (ENT1, ∼70% decrease, P = 0.0046; ENT2, ∼50% decrease, P = 0.0012). NBMPR also decreased ribavirin accumulation in ENT1 and ENT2 cells (ENT1: ∼50% decrease, P = 0.0498; ENT2: ∼30% decrease, P = 0.0125). Darunavir mapped to the ENT1 inhibitor pharmacophore and NBMPR did not significantly influence darunavir accumulation in either ENT1 or ENT2 cells (ENT1: P = 0.28; ENT2: P = 0.53), indicating that darunavir's interaction with the ENTs is limited to inhibition. These computational and in vitro models can inform compound selection in the drug discovery and development process, thereby reducing time and expense of identification and optimization of ENT-interacting compounds. SIGNIFICANCE STATEMENT: This study developed computational models of human equilibrative nucleoside transporters (ENTs) to predict drug interactions and validated these models with two compounds in vitro. Identification and prediction of ENT1 and ENT2 substrates allows for the determination of drugs that can penetrate tissues expressing these transporters.

The Antiviral Drug Tilorone Is a Potent and Selective Inhibitor of Acetylcholinesterase.

Acetylcholinesterase (AChE) is an important drug target in neurological disorders like Alzheimer's disease, Lewy body dementia, and Parkinson's disease dementia as well as for other conditions like myasthenia gravis and anticholinergic poisoning. In this study, we have used a combination of high-throughput screening, machine learning, and docking to identify new inhibitors of this enzyme. Bayesian machine learning models were generated with literature data from ChEMBL for eel and human AChE inhibitors as well as butyrylcholinesterase inhibitors (BuChE) and compared with other machine learning methods. High-throughput screens for the eel AChE inhibitor model identified several molecules including tilorone, an antiviral drug that is well-established outside of the United States, as a newly identified nanomolar AChE inhibitor. We have described how tilorone inhibits both eel and human AChE with IC50's of 14.4 nM and 64.4 nM, respectively, but does not inhibit the closely related BuChE IC50 > 50 µM. We have docked tilorone into the human AChE crystal structure and shown that this selectivity is likely due to the reliance on a specific interaction with a hydrophobic residue in the peripheral anionic site of AChE that is absent in BuChE. We also conducted a pharmacological safety profile (SafetyScreen44) and kinase selectivity screen (SelectScreen) that showed tilorone (1 µM) only inhibited AChE out of 44 toxicology target proteins evaluated and did not appreciably inhibit any of the 485 kinases tested. This study suggests there may be a potential role for repurposing tilorone or its derivatives in conditions that benefit from AChE inhibition.

Bioactivity Comparison across Multiple Machine Learning Algorithms Using over 5000 Datasets for Drug Discovery.

Machine learning methods are attracting considerable attention from the pharmaceutical industry for use in drug discovery and applications beyond. In recent studies, we and others have applied multiple machine learning algorithms and modeling metrics and, in some cases, compared molecular descriptors to build models for individual targets or properties on a relatively small scale. Several research groups have used large numbers of datasets from public databases such as ChEMBL in order to evaluate machine learning methods of interest to them. The largest of these types of studies used on the order of 1400 datasets. We have now extracted well over 5000 datasets from CHEMBL for use with the ECFP6 fingerprint and in comparison of our proprietary software Assay Central with random forest, k-nearest neighbors, support vector classification, naïve Bayesian, AdaBoosted decision trees, and deep neural networks (three layers). Model performance was assessed using an array of fivefold cross-validation metrics including area-under-the-curve, F1 score, Cohen's kappa, and Matthews correlation coefficient. Based on ranked normalized scores for the metrics or datasets, all methods appeared comparable, while the distance from the top indicated that Assay Central and support vector classification were comparable. Unlike prior studies which have placed considerable emphasis on deep neural networks (deep learning), no advantage was seen in this case. If anything, Assay Central may have been at a slight advantage as the activity cutoff for each of the over 5000 datasets representing over 570,000 unique compounds was based on Assay Central performance, although support vector classification seems to be a strong competitor. We also applied Assay Central to perform prospective predictions for the toxicity targets PXR and hERG to further validate these models. This work appears to be the largest scale comparison of these machine learning algorithms to date. Future studies will likely evaluate additional databases, descriptors, and machine learning algorithms and further refine the methods for evaluating and comparing such models.

Defending Antiviral Cationic Amphiphilic Drugs That May Cause Drug-Induced Phospholipidosis.

A recent publication in Science has proposed that cationic amphiphilic drugs repurposed for COVID-19 typically use phosholipidosis as their antiviral mechanism of action in cells but will have no in vivo efficacy. On the contrary, our viewpoint, supported by additional experimental data for similar cationic amphiphilic drugs, indicates that many of these molecules have both in vitro and in vivo efficacy with no reported phospholipidosis, and therefore, this class of compounds should not be avoided but further explored, as we continue the search for broad spectrum antivirals.

Quantum Machine Learning Algorithms for Drug Discovery Applications.

The growing quantity of public and private data sets focused on small molecules screened against biological targets or whole organisms provides a wealth of drug discovery relevant data. This is matched by the availability of machine learning algorithms such as Support Vector Machines (SVM) and Deep Neural Networks (DNN) that are computationally expensive to perform on very large data sets with thousands of molecular descriptors. Quantum computer (QC) algorithms have been proposed to offer an approach to accelerate quantum machine learning over classical computer (CC) algorithms, however with significant limitations. In the case of cheminformatics, which is widely used in drug discovery, one of the challenges to overcome is the need for compression of large numbers of molecular descriptors for use on a QC. Here, we show how to achieve compression with data sets using hundreds of molecules (SARS-CoV-2) to hundreds of thousands of molecules (whole cell screening data sets for plague and M. tuberculosis) with SVM and the data reuploading classifier (a DNN equivalent algorithm) on a QC benchmarked against CC and hybrid approaches. This study illustrates the steps needed in order to be "quantum computer ready" in order to apply quantum computing to drug discovery and to provide the foundation on which to build this field.

Development of Machine Learning Models and the Discovery of a New Antiviral Compound against Yellow Fever Virus.

Yellow fever (YF) is an acute viral hemorrhagic disease transmitted by infected mosquitoes. Large epidemics of YF occur when the virus is introduced into heavily populated areas with high mosquito density and low vaccination coverage. The lack of a specific small molecule drug treatment against YF as well as for homologous infections, such as zika and dengue, highlights the importance of these flaviviruses as a public health concern. With the advancement in computer hardware and bioactivity data availability, new tools based on machine learning methods have been introduced into drug discovery, as a means to utilize the growing high throughput screening (HTS) data generated to reduce costs and increase the speed of drug development. The use of predictive machine learning models using previously published data from HTS campaigns or data available in public databases, can enable the selection of compounds with desirable bioactivity and absorption, distribution, metabolism, and excretion profiles. In this study, we have collated cell-based assay data for yellow fever virus from the literature and public databases. The data were used to build predictive models with several machine learning methods that could prioritize compounds for in vitro testing. Five molecules were prioritized and tested in vitro from which we have identified a new pyrazolesulfonamide derivative with EC50 3.2 µM and CC50 24 µM, which represents a new scaffold suitable for hit-to-lead optimization that can expand the available drug discovery candidates for YF.

Machine Learning Models Identify Inhibitors of SARS-CoV-2.

With the rapidly evolving SARS-CoV-2 variants of concern, there is an urgent need for the discovery of further treatments for the coronavirus disease (COVID-19). Drug repurposing is one of the most rapid strategies for addressing this need, and numerous compounds have already been selected for in vitro testing by several groups. These have led to a growing database of molecules with in vitro activity against the virus. Machine learning models can assist drug discovery through prediction of the best compounds based on previously published data. Herein, we have implemented several machine learning methods to develop predictive models from recent SARS-CoV-2 in vitro inhibition data and used them to prioritize additional FDA-approved compounds for in vitro testing selected from our in-house compound library. From the compounds predicted with a Bayesian machine learning model, lumefantrine, an antimalarial was selected for testing and showed limited antiviral activity in cell-based assays while demonstrating binding (Kd 259 nM) to the spike protein using microscale thermophoresis. Several other compounds which we prioritized have since been tested by others and were also found to be active in vitro. This combined machine learning and in vitro testing approach can be expanded to virtually screen available molecules with predicted activity against SARS-CoV-2 reference WIV04 strain and circulating variants of concern. In the process of this work, we have created multiple iterations of machine learning models that can be used as a prioritization tool for SARS-CoV-2 antiviral drug discovery programs. The very latest model for SARS-CoV-2 with over 500 compounds is now freely available at www.assaycentral.org.

Exploiting machine learning for end-to-end drug discovery and development.

A variety of machine learning methods such as naive Bayesian, support vector machines and more recently deep neural networks are demonstrating their utility for drug discovery and development. These leverage the generally bigger datasets created from high-throughput screening data and allow prediction of bioactivities for targets and molecular properties with increased levels of accuracy. We have only just begun to exploit the potential of these techniques but they may already be fundamentally changing the research process for identifying new molecules and/or repurposing old drugs. The integrated application of such machine learning models for end-to-end (E2E) application is broadly relevant and has considerable implications for developing future therapies and their targeting.

Comparing Machine Learning Algorithms for Predicting Drug-Induced Liver Injury (DILI).

Drug-induced liver injury (DILI) is one the most unpredictable adverse reactions to xenobiotics in humans and the leading cause of postmarketing withdrawals of approved drugs. To date, these drugs have been collated by the FDA to form the DILIRank database, which classifies DILI severity and potential. These classifications have been used by various research groups in generating computational predictions for this type of liver injury. Recently, groups from Pfizer and AstraZeneca have collated DILI in vitro data and physicochemical properties for compounds that can be used along with data from the FDA to build machine learning models for DILI. In this study, we have used these data sets, as well as the Biopharmaceutics Drug Disposition Classification System data set, to generate Bayesian machine learning models with our in-house software, Assay Central. The performance of all machine learning models was assessed through both the internal 5-fold cross-validation metrics and prediction accuracy of an external test set of compounds with known hepatotoxicity. The best-performing Bayesian model was based on the DILI-concern category from the DILIRank database with an ROC of 0.814, a sensitivity of 0.741, a specificity of 0.755, and an accuracy of 0.746. A comparison of alternative machine learning algorithms, such as k-nearest neighbors, support vector classification, AdaBoosted decision trees, and deep learning methods, produced similar statistics to those generated with the Bayesian algorithm in Assay Central. This study demonstrates machine learning models grouped in a tool called MegaTox that can be used to predict early-stage clinical compounds, as well as recent FDA-approved drugs, to identify potential DILI.

Correction to: Tilorone: a Broad-Spectrum Antiviral Invented in the USA and Commercialized in Russia and beyond.

The mentioned Tables below had data in the bottom half (right hand side) of each table that was formatted incorrectly in the final proof version versus the submitted version.

Tilorone: a Broad-Spectrum Antiviral Invented in the USA and Commercialized in Russia and beyond.

For the last 50 years we have known of a broad-spectrum agent tilorone dihydrochloride (Tilorone). This is a small-molecule orally bioavailable drug that was originally discovered in the USA and is currently used clinically as an antiviral in Russia and the Ukraine. Over the years there have been numerous clinical and non-clinical reports of its broad spectrum of antiviral activity. More recently we have identified additional promising antiviral activities against Middle East Respiratory Syndrome, Chikungunya, Ebola and Marburg which highlights that this old drug may have other uses against new viruses. This may in turn inform the types of drugs that we need for virus outbreaks such as for the new coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Tilorone has been long neglected by the west in many respects but it deserves further reassessment in light of current and future needs for broad-spectrum antivirals.

Machine Learning Models for Estrogen Receptor Bioactivity and Endocrine Disruption Prediction.

The U.S. Environmental Protection Agency (EPA) periodically releases in vitro data across a variety of targets, including the estrogen receptor (ER). In 2015, the EPA used these data to construct mathematical models of ER agonist and antagonist pathways to prioritize chemicals for endocrine disruption testing. However, mathematical models require in vitro data prior to predicting estrogenic activity, but machine learning methods are capable of prospective prediction from the molecular structure alone. The current study describes the generation and evaluation of Bayesian machine learning models grouped by the EPA's ER agonist pathway model using multiple data types with proprietary software, Assay Central. External predictions with three test sets of in vitro and in vivo reference chemicals with agonist activity classifications were compared to previous mathematical model publications. Training data sets were subjected to additional machine learning algorithms and compared with rank normalized scores of internal five-fold cross-validation statistics. External predictions were found to be comparable or superior to previous studies published by the EPA. When assessing six additional algorithms for the training data sets, Assay Central performed similarly at a reduced computational cost. This study demonstrates that machine learning can prioritize chemicals for future in vitro and in vivo testing of ER agonism.

Comparing Machine Learning Models for Aromatase (P450 19A1).

Aromatase, or cytochrome P450 19A1, catalyzes the aromatization of androgens to estrogens within the body. Changes in the activity of this enzyme can produce hormonal imbalances that can be detrimental to sexual and skeletal development. Inhibition of this enzyme can occur with drugs and natural products as well as environmental chemicals. Therefore, predicting potential endocrine disruption via exogenous chemicals requires that aromatase inhibition be considered in addition to androgen and estrogen pathway interference. Bayesian machine learning methods can be used for prospective prediction from the molecular structure without the need for experimental data. Herein, the generation and evaluation of multiple machine learning models utilizing different sources of aromatase inhibition data are described. These models are applied to two test sets for external validation with molecules relevant to drug discovery from the public domain. In addition, the performance of multiple machine learning algorithms was evaluated by comparing internal five-fold cross-validation statistics of the training data. These methods to predict aromatase inhibition from molecular structure, when used in concert with estrogen and androgen machine learning models, allow for a more holistic assessment of endocrine-disrupting potential of chemicals with limited empirical data and enable the reduction of the use of hazardous substances.