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.

Machine learning-aided search for ligands of P2Y6 and other P2Y receptors.

The P2Y6 receptor, activated by uridine diphosphate (UDP), is a target for antagonists in inflammatory, neurodegenerative, and metabolic disorders, yet few potent and selective antagonists are known to date. This prompted us to use machine learning as a novel approach to aid ligand discovery, with pharmacological evaluation at three P2YR subtypes: initially P2Y6 and subsequently P2Y1 and P2Y14. Relying on extensive published data for P2Y6R agonists, we generated and validated an array of classification machine learning model using the algorithms deep learning (DL), adaboost classifier (ada), Bernoulli NB (bnb), k-nearest neighbors (kNN) classifier, logistic regression (lreg), random forest classifier (rf), support vector classification (SVC), and XGBoost (XGB) classifier models, and the common consensus was applied to molecular selection of 21 diverse structures. Compounds were screened using human P2Y6R-induced functional calcium transients in transfected 1321N1 astrocytoma cells and fluorescent binding inhibition at closely related hP2Y14R expressed in CHO cells. The hit compound ABBV-744, an experimental anticancer drug with a 6-methyl-7-oxo-6,7-dihydro-1H-pyrrolo[2,3-c]pyridine scaffold, had multifaceted interactions with the P2YR family: hP2Y6R inhibition in a non-surmountable fashion, suggesting that noncompetitive antagonism, and hP2Y1R enhancement, but not hP2Y14R binding inhibition. Other machine learning-selected compounds were either weak (experimental anti-asthmatic drug AZD5423 with a phenyl-1H-indazole scaffold) or inactive in inhibiting the hP2Y6R. Experimental drugs TAK-593 and GSK1070916 (100 µM) inhibited P2Y14R fluorescent binding by 50% and 38%, respectively, and all other compounds by < 20%. Thus, machine learning has led the way toward revealing previously unknown modulators of several P2YR subtypes that have varied effects.

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.

Preventing AI From Creating Biochemical Threats.

We have previously applied our machine learning models for bioactivity and toxicity along with a generative algorithm to develop VX and tens of thousands of analogues. The publication brought attention to the ease of designing chemical warfare agents. In this Viewpoint, we discuss 10 recommendations to prevent future biochemical threats.

The protein disulfide isomerase inhibitor 3-methyltoxoflavin inhibits Chikungunya virus.

Chikungunya virus (CHIKV) is the etiological agent of chikungunya fever, a (re)emerging arbovirus infection, that causes severe and often persistent arthritis, as well as representing a serious health concern worldwide for which no antivirals are currently available. Despite efforts over the last decade to identify and optimize new inhibitors or to reposition existing drugs, no compound has progressed to clinical trials for CHIKV and current prophylaxis is based on vector control, which has shown limited success in containing the virus. Our efforts to rectify this situation were initiated by screening 36 compounds using a replicon system and ultimately identified the natural product derivative 3-methyltoxoflavin with activity against CHIKV using a cell-based assay (EC50 200 nM, SI = 17 in Huh-7 cells). We have additionally screened 3-methyltoxoflavin against a panel of 17 viruses and showed that it only additionally demonstrated inhibition of the yellow fever virus (EC50 370 nM, SI = 3.2 in Huh-7 cells). We have also showed that 3-methyltoxoflavin has excellent in vitro human and mouse microsomal metabolic stability, good solubility and high Caco-2 permeability and it is not likely to be a P-glycoprotein substrate. In summary, we demonstrate that 3-methyltoxoflavin has activity against CHIKV, good in vitro absorption, distribution, metabolism and excretion (ADME) properties as well as good calculated physicochemical properties and may represent a valuable starting point for future optimization to develop inhibitors for this and other related viruses.

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.

Lack of an Effect of Polysorbate 80 on Intestinal Drug Permeability in Humans.

PURPOSE: Despite no broad, direct evidence in humans, there is a potential concern that surfactants alter active or passive drug intestinal permeation to modulate oral drug absorption. The purpose of this study was to investigate the impact of the surfactant polysorbate 80 on active and passive intestinal drug absorption in humans. METHODS: The human (n = 12) pharmacokinetics (PK) of three probe substrates of intestinal absorption, valacyclovir, chenodeoxycholic acid (CDCA), and enalaprilat, were assessed. Endogenous bile acid levels were assessed as a secondary measure of transporter and microbiota impact. RESULTS: Polysorbate 80 did not inhibit peptide transporter 1 (PepT1)- or apical sodium bile acid transporter (ASBT)-mediated PK of valacyclovir and CDCA, respectively. Polysorbate 80 did not increase enalaprilat absorption. Modest increases in unconjugated secondary bile acid Cmax ratios suggest a potential alteration of the in vivo intestinal microbiota by polysorbate 80. CONCLUSIONS: Polysorbate 80 did not alter intestinal membrane fluidity or cause intestinal membrane disruption. This finding supports regulatory relief of excipient restrictions for Biopharmaceutics Classification System-based biowaivers.

Discovery of New Zika Protease and Polymerase Inhibitors through the Open Science Collaboration Project OpenZika.

The Zika virus (ZIKV) is a neurotropic arbovirus considered a global threat to public health. Although there have been several efforts in drug discovery projects for ZIKV in recent years, there are still no antiviral drugs approved to date. Here, we describe the results of a global collaborative crowdsourced open science project, the OpenZika project, from IBM's World Community Grid (WCG), which integrates different computational and experimental strategies for advancing a drug candidate for ZIKV. Initially, molecular docking protocols were developed to identify potential inhibitors of ZIKV NS5 RNA-dependent RNA polymerase (NS5 RdRp), NS3 protease (NS2B-NS3pro), and NS3 helicase (NS3hel). Then, a machine learning (ML) model was built to distinguish active vs inactive compounds for the cytoprotective effect against ZIKV infection. We performed three independent target-based virtual screening campaigns (NS5 RdRp, NS2B-NS3pro, and NS3hel), followed by predictions by the ML model and other filters, and prioritized a total of 61 compounds for further testing in enzymatic and phenotypic assays. This yielded five non-nucleoside compounds which showed inhibitory activity against ZIKV NS5 RdRp in enzymatic assays (IC50 range from 0.61 to 17 µM). Two compounds thermally destabilized NS3hel and showed binding affinity in the micromolar range (Kd range from 9 to 35 µM). Moreover, the compounds LabMol-301 inhibited both NS5 RdRp and NS2B-NS3pro (IC50 of 0.8 and 7.4 µM, respectively) and LabMol-212 thermally destabilized the ZIKV NS3hel (Kd of 35 µM). Both also protected cells from death induced by ZIKV infection in in vitro cell-based assays. However, while eight compounds (including LabMol-301 and LabMol-212) showed a cytoprotective effect and prevented ZIKV-induced cell death, agreeing with our ML model for prediction of this cytoprotective effect, no compound showed a direct antiviral effect against ZIKV. Thus, the new scaffolds discovered here are promising hits for future structural optimization and for advancing the discovery of further drug candidates for ZIKV. Furthermore, this work has demonstrated the importance of the integration of computational and experimental approaches, as well as the potential of large-scale collaborative networks to advance drug discovery projects for neglected diseases and emerging viruses, despite the lack of available direct antiviral activity and cytoprotective effect data, that reflects on the assertiveness of the computational predictions. The importance of these efforts rests with the need to be prepared for future viral epidemic and pandemic outbreaks.

Multiple approaches to repurposing drugs for neuroblastoma.

Neuroblastoma (NB) is the second leading extracranial solid tumor of early childhood with about two-thirds of cases presenting before the age of 5, and accounts for roughly 15 percent of all pediatric cancer fatalities in the United States. Treatments against NB are lacking, resulting in a low survival rate in high-risk patients. A repurposing approach using already approved or clinical stage compounds can be used for diseases for which the patient population is small, and the commercial market limited. We have used Bayesian machine learning, in vitro cell assays, and combination analysis to identify molecules with potential use for NB. We demonstrated that pyronaridine (SH-SY5Y IC50 1.70 µM, SK-N-AS IC50 3.45 µM), BAY 11-7082 (SH-SY5Y IC50 0.85 µM, SK-N-AS IC50 1.23 µM), niclosamide (SH-SY5Y IC50 0.87 µM, SK-N-AS IC50 2.33 µM) and fingolimod (SH-SY5Y IC50 4.71 µM, SK-N-AS IC50 6.11 µM) showed cytotoxicity against NB. As several of the molecules are approved drugs in the US or elsewhere, they may be repurposed more readily for NB treatment. Pyronaridine was also tested in combinations in SH-SY5Y cells and demonstrated an antagonistic effect with either etoposide or crizotinib. Whereas when crizotinib and etoposide were combined with each other they had a synergistic effect in these cells. We have also described several analogs of pyronaridine to explore the structure-activity relationship against cell lines. We describe multiple molecules demonstrating cytotoxicity against NB and the further evaluation of these molecules and combinations using other NB cells lines and in vivo models will be important in the future to assess translational potential.

Chalcones from Angelica keiskei (ashitaba) inhibit key Zika virus replication proteins.

Zika virus (ZIKV) is a dangerous human pathogen and no antiviral drugs have been approved to date. The chalcones are a group of small molecules that are found in a number of different plants, including Angelica keiskei Koidzumi, also known as ashitaba. To examine chalcone anti-ZIKV activity, three chalcones, 4-hydroxyderricin (4HD), xanthoangelol (XA), and xanthoangelol-E (XA-E), were purified from a methanol-ethyl acetate extract from A. keiskei. Molecular and ensemble docking predicted that these chalcones would establish multiple interactions with residues in the catalytic and allosteric sites of ZIKV NS2B-NS3 protease, and in the allosteric site of the NS5 RNA-dependent RNA-polymerase (RdRp). Machine learning models also predicted 4HD, XA and XA-E as potential anti-ZIKV inhibitors. Enzymatic and kinetic assays confirmed chalcone inhibition of the ZIKV NS2B-NS3 protease allosteric site with IC50s from 18 to 50 µM. Activity assays also revealed that XA, but not 4HD or XA-E, inhibited the allosteric site of the RdRp, with an IC50 of 6.9 µM. Finally, we tested these chalcones for their anti-viral activity in vitro with Vero cells. 4HD and XA-E displayed anti-ZIKV activity with EC50 values of 6.6 and 22.0 µM, respectively, while XA displayed relatively weak anti-ZIKV activity with whole cells. With their simple structures and relative ease of modification, the chalcones represent attractive candidates for hit-to-lead optimization in the search of new anti-ZIKV therapeutics.

A critical overview of computational approaches employed for COVID-19 drug discovery.

COVID-19 has resulted in huge numbers of infections and deaths worldwide and brought the most severe disruptions to societies and economies since the Great Depression. Massive experimental and computational research effort to understand and characterize the disease and rapidly develop diagnostics, vaccines, and drugs has emerged in response to this devastating pandemic and more than 130 000 COVID-19-related research papers have been published in peer-reviewed journals or deposited in preprint servers. Much of the research effort has focused on the discovery of novel drug candidates or repurposing of existing drugs against COVID-19, and many such projects have been either exclusively computational or computer-aided experimental studies. Herein, we provide an expert overview of the key computational methods and their applications for the discovery of COVID-19 small-molecule therapeutics that have been reported in the research literature. We further outline that, after the first year the COVID-19 pandemic, it appears that drug repurposing has not produced rapid and global solutions. However, several known drugs have been used in the clinic to cure COVID-19 patients, and a few repurposed drugs continue to be considered in clinical trials, along with several novel clinical candidates. We posit that truly impactful computational tools must deliver actionable, experimentally testable hypotheses enabling the discovery of novel drugs and drug combinations, and that open science and rapid sharing of research results are critical to accelerate the development of novel, much needed therapeutics for COVID-19.

Remdesivir and EIDD-1931 Interact with Human Equilibrative Nucleoside Transporters 1 and 2: Implications for Reaching SARS-CoV-2 Viral Sanctuary Sites.

Equilibrative nucleoside transporters (ENTs) are present at the blood-testis barrier (BTB), where they can facilitate antiviral drug disposition to eliminate a sanctuary site for viruses detectable in semen. The purpose of this study was to investigate ENT-drug interactions with three nucleoside analogs, remdesivir, molnupiravir, and molnupiravir's active metabolite, ß-d-N4-hydroxycytidine (EIDD-1931), and four non-nucleoside molecules repurposed as antivirals for coronavirus disease 2019 (COVID-19). The study used three-dimensional pharmacophores for ENT1 and ENT2 substrates and inhibitors and Bayesian machine learning models to identify potential interactions with these transporters. In vitro transport experiments demonstrated that remdesivir was the most potent inhibitor of ENT-mediated [3H]uridine uptake (ENT1 IC50: 39 µM; ENT2 IC50: 77 µM), followed by EIDD-1931 (ENT1 IC50: 259 µM; ENT2 IC50: 467 µM), whereas molnupiravir was a modest inhibitor (ENT1 IC50: 701 µM; ENT2 IC50: 851 µM). Other proposed antivirals failed to inhibit ENT-mediated [3H]uridine uptake below 1 mM. Remdesivir accumulation decreased in the presence of 6-S-[(4-nitrophenyl)methyl]-6-thioinosine (NBMPR) by 30% in ENT1 cells (P = 0.0248) and 27% in ENT2 cells (P = 0.0054). EIDD-1931 accumulation decreased in the presence of NBMPR by 77% in ENT1 cells (P = 0.0463) and by 64% in ENT2 cells (P = 0.0132), which supported computational predictions that both are ENT substrates that may be important for efficacy against COVID-19. NBMPR failed to decrease molnupiravir uptake, suggesting that ENT interaction is likely inhibitory. Our combined computational and in vitro data can be used to identify additional ENT-drug interactions to improve our understanding of drugs that can circumvent the BTB. SIGNIFICANCE STATEMENT: This study identified remdesivir and EIDD-1931 as substrates of equilibrative nucleoside transporters 1 and 2. This provides a potential mechanism for uptake of these drugs into cells and may be important for antiviral potential in the testes and other tissues expressing these transporters.

Predicting Drug Interactions with Human Equilibrative Nucleoside Transporters 1 and 2 Using Functional Knockout Cell Lines and Bayesian Modeling.

Equilibrative nucleoside transporters (ENTs) 1 and 2 facilitate nucleoside transport across the blood-testis barrier (BTB). Improving drug entry into the testes with drugs that use endogenous transport pathways may lead to more effective treatments for diseases within the reproductive tract. In this study, CRISPR/CRISPR-associated protein 9 was used to generate HeLa cell lines in which ENT expression was limited to ENT1 or ENT2. We characterized uridine transport in these cell lines and generated Bayesian models to predict interactions with the ENTs. Quantification of [3H]uridine uptake in the presence of the ENT-specific inhibitor S-(4-nitrobenzyl)-6-thioinosine (NBMPR) demonstrated functional loss of each transporter. Nine nucleoside reverse-transcriptase inhibitors and 37 nucleoside/heterocycle analogs were evaluated to identify ENT interactions. Twenty-one compounds inhibited uridine uptake and abacavir, nevirapine, ticagrelor, and uridine triacetate had different IC50 values for ENT1 and ENT2. Total accumulation of four identified inhibitors was measured with and without NBMPR to determine whether there was ENT-mediated transport. Clofarabine and cladribine were ENT1 and ENT2 substrates, whereas nevirapine and lexibulin were ENT1 and ENT2 nontransported inhibitors. Bayesian models generated using Assay Central machine learning software yielded reasonably high internal validation performance (receiver operator characteristic > 0.7). ENT1 IC50-based models were generated from ChEMBL; subvalidations using this training data set correctly predicted 58% of inhibitors when analyzing activity by percent uptake and 63% when using estimated-IC50 values. Determining drug interactions with these transporters can be useful in identifying and predicting compounds that are ENT1 and ENT2 substrates and can thereby circumvent the BTB through this transepithelial transport pathway in Sertoli cells. SIGNIFICANCE STATEMENT: This study is the first to predict drug interactions with equilibrative nucleoside transporter (ENT) 1 and ENT2 using Bayesian modeling. Novel CRISPR/CRISPR-associated protein 9 functional knockouts of ENT1 and ENT2 in HeLa S3 cells were generated and characterized. Determining drug interactions with these transporters can be useful in identifying and predicting compounds that are ENT1 and ENT2 substrates and can circumvent the blood-testis barrier through this transepithelial transport pathway in Sertoli cells.

UV-adVISor: Attention-Based Recurrent Neural Networks to Predict UV-Vis Spectra.

Ultraviolet-visible (UV-Vis) absorption spectra are routinely collected as part of high-performance liquid chromatography (HPLC) analysis systems and can be used to identify chemical reaction products by comparison to the reference spectra. Here, we present UV-adVISor as a new computational tool for predicting the UV-Vis spectra from a molecule's structure alone. UV-Vis prediction was approached as a sequence-to-sequence problem. We utilized Long-Short Term Memory and attention-based neural networks with Extended Connectivity Fingerprint Diameter 6 or molecule SMILES to generate predictive models for the UV spectra. We have produced two spectrum datasets (dataset I, N = 949, and dataset II, N = 2222) using different compound collections and spectrum acquisition methods to train, validate, and test our models. We evaluated the prediction accuracy of the complete spectra by the correspondence of wavelengths of absorbance maxima and with a series of statistical measures (the best test set median model parameters are in parentheses for model II), including RMSE (0.064), R2 (0.71), and dynamic time warping (DTW, 0.194) of the entire spectrum curve. Scrambling molecule structures with the experimental spectra during training resulted in a degraded R2, confirming the utility of the approaches for prediction. UV-adVISor is able to provide fast and accurate predictions for libraries of compounds.

Cationic Compounds with SARS-CoV-2 Antiviral Activity and Their Interaction with Organic Cation Transporter/Multidrug and Toxin Extruder Secretory Transporters.

In the wake of the COVID-19 pandemic, drug repurposing has been highlighted for rapid introduction of therapeutics. Proposed drugs with activity against SARS-CoV-2 include compounds with positive charges at physiologic pH, making them potential targets for the organic cation secretory transporters of kidney and liver, i.e., the basolateral organic cation transporters, OCT1 and OCT2; and the apical multidrug and toxin extruders, MATE1 and MATE2-K. We selected several compounds proposed to have in vitro activity against SARS-CoV-2 (chloroquine, hydroxychloroquine, quinacrine, tilorone, pyronaridine, cetylpyridinium, and miramistin) to test their interaction with OCT and MATE transporters. We used Bayesian machine learning models to generate predictions for each molecule with each transporter and also experimentally determined IC50 values for each compound against labeled substrate transport into CHO cells that stably expressed OCT2, MATE1, or MATE2-K using three structurally distinct substrates (atenolol, metformin and 1-methyl-4-phenylpyridinium) to assess the impact of substrate structure on inhibitory efficacy. For the OCTs substrate identity influenced IC50 values, although the effect was larger and more systematic for OCT2. In contrast, inhibition of MATE1-mediated transport was largely insensitive to substrate identity. Unlike MATE1, inhibition of MATE2-K was influenced, albeit modestly, by substrate identity. Maximum unbound plasma concentration/IC50 ratios were used to identify potential clinical DDI recommendations; all the compounds interacted with the OCT/MATE secretory pathway, most with sufficient avidity to represent potential DDI issues for secretion of cationic drugs. This should be considered when proposing cationic agents as repurposed antivirals. SIGNIFICANCE STATEMENT: Drugs proposed as potential COVID-19 therapeutics based on in vitro activity data against SARS-CoV-2 include compounds with positive charges at physiological pH, making them potential interactors with the OCT/MATE renal secretory pathway. We tested seven such molecules as inhibitors of OCT1/2 and MATE1/2-K. All the compounds blocked transport activity regardless of substrate used to monitor activity. Suggesting that plasma concentrations achieved by normal clinical application of the test agents could be expected to influence the pharmacokinetics of selected cationic drugs.

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.