SynRoute: A Retrosynthetic Planning Software.

Computer-assisted synthetic planning has seen major advancements that stem from the availability of large reaction databases and artificial intelligence methodologies. SynRoute is a new retrosynthetic planning software tool that uses a relatively small number of general reaction templates, currently 263, along with a literature-based reaction database to find short, practical synthetic routes for target compounds. For each reaction template, a machine learning classifier is trained using data from the Pistachio reaction database to predict whether new computer-generated reactions based on the template are likely to work experimentally in the laboratory. This reaction generation methodology is used together with a vectorized Dijkstra-like search of top-scoring routes organized by synthetic strategies for easy browsing by a synthetic chemist. SynRoute was able to find routes for an average of 83% of compounds based on selection of random subsets of drug-like compounds from the ChEMBL database. Laboratory evaluation of 12 routes produced by SynRoute, to synthesize compounds not from the previous random subsets, demonstrated the ability to produce feasible overall synthetic strategies for all compounds evaluated.

Machine Learning Models for Predicting Molecular UV-Vis Spectra with Quantum Mechanical Properties.

Accurate understanding of ultraviolet-visible (UV-vis) spectra is critical for the high-throughput synthesis of compounds for drug discovery. Experimentally determining UV-vis spectra can become expensive when dealing with a large quantity of novel compounds. This provides us an opportunity to drive computational advances in molecular property predictions using quantum mechanics and machine learning methods. In this work, we use both quantum mechanically (QM) predicted and experimentally measured UV-vis spectra as input to devise four different machine learning architectures, UVvis-SchNet, UVvis-DTNN, UVvis-Transformer, and UVvis-MPNN, and assess the performance of each method. We find that the UVvis-MPNN model outperforms the other models when using optimized 3D coordinates and QM predicted spectra as input features. This model has the highest performance for predicting UV-vis spectra with a training RMSE of 0.06 and validation RMSE of 0.08. Most importantly, our model can be used for the challenging task of predicting differences in the UV-vis spectral signatures of regioisomers.

Mega-High-Throughput Screening Platform for the Discovery of Biologically Relevant Sequence-Defined Non-Natural Polymers.

Combinatorial methods enable the synthesis of chemical libraries on scales of millions to billions of compounds, but the ability to efficiently screen and sequence such large libraries has remained a major bottleneck for molecular discovery. We developed a novel technology for screening and sequencing libraries of synthetic molecules of up to a billion compounds in size. This platform utilizes the fiber-optic array scanning technology (FAST) to screen bead-based libraries of synthetic compounds at a rate of 5 million compounds per minute (∼83 000 Hz). This ultra-high-throughput screening platform has been used to screen libraries of synthetic "self-readable" non-natural polymers that can be sequenced at the femtomole scale by chemical fragmentation and high-resolution mass spectrometry. The versatility and throughput of the platform were demonstrated by screening two libraries of non-natural polyamide polymers with sizes of 1.77M and 1B compounds against the protein targets K-Ras, asialoglycoprotein receptor 1 (ASGPR), IL-6, IL-6 receptor (IL-6R), and TNFα. Hits with low nanomolar binding affinities were found against all targets, including competitive inhibitors of K-Ras binding to Raf and functionally active uptake ligands for ASGPR facilitating intracellular delivery of a nonglycan ligand.

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.

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.

Repurposing the Ebola and Marburg Virus Inhibitors Tilorone, Quinacrine, and Pyronaridine: In Vitro Activity against SARS-CoV-2 and Potential Mechanisms.

Severe acute respiratory coronavirus 2 (SARS-CoV-2) is a newly identified virus that has resulted in over 2.5 million deaths globally and over 116 million cases globally in March, 2021. Small-molecule inhibitors that reverse disease severity have proven difficult to discover. One of the key approaches that has been widely applied in an effort to speed up the translation of drugs is drug repurposing. A few drugs have shown in vitro activity against Ebola viruses and demonstrated activity against SARS-CoV-2 in vivo. Most notably, the RNA polymerase targeting remdesivir demonstrated activity in vitro and efficacy in the early stage of the disease in humans. Testing other small-molecule drugs that are active against Ebola viruses (EBOVs) would appear a reasonable strategy to evaluate their potential for SARS-CoV-2. We have previously repurposed pyronaridine, tilorone, and quinacrine (from malaria, influenza, and antiprotozoal uses, respectively) as inhibitors of Ebola and Marburg viruses in vitro in HeLa cells and mouse-adapted EBOV in mice in vivo. We have now tested these three drugs in various cell lines (VeroE6, Vero76, Caco-2, Calu-3, A549-ACE2, HUH-7, and monocytes) infected with SARS-CoV-2 as well as other viruses (including MHV and HCoV 229E). The compilation of these results indicated considerable variability in antiviral activity observed across cell lines. We found that tilorone and pyronaridine inhibited the virus replication in A549-ACE2 cells with IC50 values of 180 nM and IC50 198 nM, respectively. We used microscale thermophoresis to test the binding of these molecules to the spike protein, and tilorone and pyronaridine bind to the spike receptor binding domain protein with K d values of 339 and 647 nM, respectively. Human Cmax for pyronaridine and quinacrine is greater than the IC50 observed in A549-ACE2 cells. We also provide novel insights into the mechanism of these compounds which is likely lysosomotropic.

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.

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.

Repurposing Pyramax®, quinacrine and tilorone as treatments for Ebola virus disease.

We have recently identified three molecules (tilorone, quinacrine and pyronaridine tetraphosphate) which all demonstrated efficacy in the mouse model of infection with mouse-adapted Ebola virus (EBOV) model of disease and had similar in vitro inhibition of an Ebola pseudovirus (VSV-EBOV-GP), suggesting they interfere with viral entry. Using a machine learning model to predict lysosomotropism these compounds were evaluated for their ability to possess a lysosomotropic mechanism in vitro. We now demonstrate in vitro that pyronaridine tetraphosphate is an inhibitor of Lysotracker accumulation in lysosomes (IC50 = 0.56 µM). Further, we evaluated antiviral synergy between pyronaridine and artesunate (Pyramax®), which are used in combination to treat malaria. Artesunate was not found to have lysosomotropic activity in vitro and the combination effect on EBOV inhibition was shown to be additive. Pyramax® may represent a unique example of the repurposing of a combination product for another disease.

Pyronaridine tetraphosphate efficacy against Ebola virus infection in guinea pig.

The recent outbreaks of the Ebola virus (EBOV) in Africa have brought global visibility to the shortage of available therapeutic options to treat patients infected with this or closely related viruses. We have recently computationally identified three molecules which have all demonstrated statistically significant efficacy in the mouse model of infection with mouse adapted Ebola virus (ma-EBOV). One of these molecules is the antimalarial pyronaridine tetraphosphate (IC50 range of 0.82-1.30 µM against three strains of EBOV and IC50 range of 1.01-2.72 µM against two strains of Marburg virus (MARV)) which is an approved drug in the European Union and used in combination with artesunate. To date, no small molecule drugs have shown statistically significant efficacy in the guinea pig model of EBOV infection. Pharmacokinetics and range-finding studies in guinea pigs directed us to a single 300 mg/kg or 600 mg/kg oral dose of pyronaridine 1hr after infection. Pyronaridine resulted in statistically significant survival of 40% at 300 mg/kg and protected from a lethal challenge with EBOV. In comparison, oral favipiravir (300 mg/kg dosed once a day) had 43.5% survival. All animals in the vehicle treatment group succumbed to disease by study day 12 (100% mortality). The in vitro metabolism and metabolite identification of pyronaridine and another of our EBOV active molecules, tilorone, suggested significant species differences which may account for the efficacy or lack thereof, respectively in guinea pig. In summary, our studies with pyronaridine demonstrates its utility for repurposing as an antiviral against EBOV and MARV.

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.

Tilorone, a Broad-Spectrum Antiviral for Emerging Viruses.

Tilorone is a 50-year-old synthetic small-molecule compound with antiviral activity that is proposed to induce interferon after oral administration. This drug is used as a broad-spectrum antiviral in several countries of the Russian Federation. We have recently described activity in vitro and in vivo against the Ebola virus. After a broad screening of additional viruses, we now describe in vitro activity against Chikungunya virus (CHIK) and Middle Eastern respiratory syndrome coronavirus (MERS-CoV).

Repurposing the antimalarial pyronaridine tetraphosphate to protect against Ebola virus infection.

Recent outbreaks of the Ebola virus (EBOV) have focused attention on the dire need for antivirals to treat these patients. We identified pyronaridine tetraphosphate as a potential candidate as it is an approved drug in the European Union which is currently used in combination with artesunate as a treatment for malaria (EC50 between 420 nM-1.14 µM against EBOV in HeLa cells). Range-finding studies in mice directed us to a single 75 mg/kg i.p. dose 1 hr after infection which resulted in 100% survival and statistically significantly reduced viremia at study day 3 from a lethal challenge with mouse-adapted EBOV (maEBOV). Further, an EBOV window study suggested we could dose pyronaridine 2 or 24 hrs post-exposure to result in similar efficacy. Analysis of cytokine and chemokine panels suggests that pyronaridine may act as an immunomodulator during an EBOV infection. Our studies with pyronaridine clearly demonstrate potential utility for its repurposing as an antiviral against EBOV and merits further study in larger animal models with the added benefit of already being used as a treatment against malaria.

Repurposing Quinacrine against Ebola Virus Infection In Vivo.

Quinacrine hydrochloride is a small-molecule, orally bioavailable drug that has been used clinically as an antimalarial and for many other applications. A machine learning model trained on Ebola virus (EBOV) screening data identified quinacrine as a potent (nanomolar) in vitro inhibitor. In the current study, quinacrine 25 mg/kg was shown to protect 70% of mice (statistically significant) from a lethal challenge with mouse-adapted EBOV with once-daily intraperitoneal dosing for 8 days.

The Natural Product Eugenol Is an Inhibitor of the Ebola Virus In Vitro.

PURPOSE: Since the 2014 Ebola virus (EBOV) outbreak in West Africa there has been considerable effort towards developing drugs to treat Ebola virus disease and yet to date there is no FDA approved treatment. This is important as at the time of writing this manuscript there is an ongoing outbreak in the Democratic Republic of the Congo which has killed over 1000. METHODS: We have evaluated a small number of natural products, some of which had shown antiviral activity against other pathogens. This is exemplified with eugenol, which is found in high concentrations in multiple essential oils, and has shown antiviral activity against feline calicivirus, tomato yellow leaf curl virus, Influenza A virus, Herpes Simplex virus type 1 and 2, and four airborne phages. RESULTS: Four compounds possessed EC50 values less than or equal to 11 µM. Of these, eugenol, had an EC50 of 1.3 µM against EBOV and is present in several plants including clove, cinnamon, basil and bay. Eugenol is much smaller and structurally unlike any compound that has been previously identified as an inhibitor of EBOV, therefore it may provide new mechanistic insights. CONCLUSION: This compound is readily accessible in bulk quantities, is inexpensive, and has a long history of human consumption, which endorses the idea for further assessment as an antiviral therapeutic. This work also suggests that a more exhaustive assessment of natural product libraries against EBOV and other viruses is warranted to improve our ability to identify compounds that are so distinct from FDA approved drugs.

Ebola Virus Bayesian Machine Learning Models Enable New in Vitro Leads.

We have previously described the first Bayesian machine learning models from FDA-approved drug screens, for identifying compounds active against the Ebola virus (EBOV). These models led to the identification of three active molecules in vitro: tilorone, pyronaridine, and quinacrine. A follow-up study demonstrated that one of these compounds, tilorone, has 100% in vivo efficacy in mice infected with mouse-adapted EBOV at 30 mg/kg/day intraperitoneal. This suggested that we can learn from the published data on EBOV inhibition and use it to select new compounds for testing that are active in vivo. We used these previously built Bayesian machine learning EBOV models alongside our chemical insights for the selection of 12 molecules, absent from the training set, to test for in vitro EBOV inhibition. Nine molecules were directly selected using the model, and eight of these molecules possessed a promising in vitro activity (EC50 < 15 µM). Three further compounds were selected for an in vitro evaluation because they were antimalarials, and compounds of this class like pyronaridine and quinacrine have previously been shown to inhibit EBOV. We identified the antimalarial drug arterolane (IC50 = 4.53 µM) and the anticancer clinical candidate lucanthone (IC50 = 3.27 µM) as novel compounds that have EBOV inhibitory activity in HeLa cells and generally lack cytotoxicity. This work provides further validation for using machine learning and medicinal chemistry expertize to prioritize compounds for testing in vitro prior to more costly in vivo tests. These studies provide further corroboration of this strategy and suggest that it can likely be applied to other pathogens in the future.

Molecular Imaging of Mesothelin-Expressing Ovarian Cancer with a Human and Mouse Cross-Reactive Nanobody.

Mesothelin is an epithelial marker highly expressed at the cell surface of cancer cells from diverse origins, including ovarian and pancreatic adenocarcinomas and mesotheliomas. Previously, we identified and characterized an antimesothelin nanobody (NbG3a) for in vitro diagnostic applications. The main goal of this research was to establish the potential of NbG3a as a molecular imaging agent. Site-specific biotinylated NbG3a (bNbG3a) was bound to streptavidin-conjugated reagents for in vitro and in vivo assays. Initially, we performed microscale thermophoresis to determine the binding affinity between bNbG3a and human ( Kd = 46 ± 8 nM) or mouse ( Kd = 4.8 ± 0.4 nM) mesothelin protein. The human and mouse cross-reactivity was confirmed by in vivo optical imaging using bNbG3a bound to fluorescent streptavidin. We also localized the binding site of nNbG3a on human mesothelin using overlapping peptide scan. NbG3a recognized an epitope within residues 21-65 of the mature membrane bound form of human mesothelin, which is part of the N-terminal region of mesothelin that is important for interactions between mesothelin on peritoneal cells and CA125 on tumor cells. Next, the bNbG3a in vivo half-life after intravenous injection in healthy mice was estimated by ELISA assay to be 5.3 ± 1.3 min. In tumor-bearing animals, fluorescent bNbG3a accumulated in a subcutaneous ovarian xenograft (A1847) and in two syngeneic, orthotopic ovarian tumors (intraovary and intraperitoneal ID8) within an hour of intravenous injection that peaked by 4 h and persisted up to 48 h. MRI analysis of bNbG3a-targeted streptavidin-labeled iron oxides showed that the MRI signal intensity decreased 1 h after injection for a subcutaneous xenograft model of ovarian cancer for bNbG3a-labeled iron oxides compared to unlabeled iron oxides. The signal intensity differences continued up to the final time point at 24 h post injection. Finally, in vivo immunofluorescence 24 or 48 h after bNbG3a intravenous injection showed bNbG3a diffuse distribution of both xenograft and syngeneic ovarian tumors, with local areas of high concentration throughout A1847 human tumor. The data support the use of NbG3a for continued preclinical development and translation to human applications for cancers that overexpress mesothelin.

Efficacy of Tilorone Dihydrochloride against Ebola Virus Infection.

Tilorone dihydrochloride (tilorone) is a small-molecule, orally bioavailable drug that is used clinically as an antiviral outside the United States. A machine-learning model trained on anti-Ebola virus (EBOV) screening data previously identified tilorone as a potent in vitro EBOV inhibitor, making it a candidate for the treatment of Ebola virus disease (EVD). In the present study, a series of in vitro ADMET (absorption, distribution, metabolism, excretion, toxicity) assays demonstrated the drug has excellent solubility, high Caco-2 permeability, was not a P-glycoprotein substrate, and had no inhibitory activity against five human CYP450 enzymes (3A4, 2D6, 2C19, 2C9, and 1A2). Tilorone was shown to have 52% human plasma protein binding with excellent plasma stability and a mouse liver microsome half-life of 48 min. Dose range-finding studies in mice demonstrated a maximum tolerated single dose of 100 mg/kg of body weight. A pharmacokinetics study in mice at 2- and 10-mg/kg dose levels showed that the drug is rapidly absorbed, has dose-dependent increases in maximum concentration of unbound drug in plasma and areas under the concentration-time curve, and has a half-life of approximately 18 h in both males and females, although the exposure was ∼2.5-fold higher in male mice. Tilorone doses of 25 and 50 mg/kg proved efficacious in protecting 90% of mice from a lethal challenge with mouse-adapted with once-daily intraperitoneal (i.p.) dosing for 8 days. A subsequent study showed that 30 mg/kg/day of tilorone given i.p. starting 2 or 24 h postchallenge and continuing through day 7 postinfection was fully protective, indicating promising activity for the treatment of EVD.

Predictive modeling targets thymidylate synthase ThyX in Mycobacterium tuberculosis.

There is an urgent need to identify new treatments for tuberculosis (TB), a major infectious disease caused by Mycobacterium tuberculosis (Mtb), which results in 1.5 million deaths each year. We have targeted two essential enzymes in this organism that are promising for antibacterial therapy and reported to be inhibited by naphthoquinones. ThyX is an essential thymidylate synthase that is mechanistically and structurally unrelated to the human enzyme. DNA gyrase is a DNA topoisomerase present in bacteria and plants but not animals. The current study set out to understand the structure-activity relationships of these targets in Mtb using a combination of cheminformatics and in vitro screening. Here, we report the identification of new Mtb ThyX inhibitors, 2-chloro-3-(4-methanesulfonylpiperazin-1-yl)-1,4-dihydronaphthalene-1,4-dione) and idebenone, which show modest whole-cell activity and appear to act, at least in part, by targeting ThyX in Mtb.

Combining Metabolite-Based Pharmacophores with Bayesian Machine Learning Models for Mycobacterium tuberculosis Drug Discovery.

Integrated computational approaches for Mycobacterium tuberculosis (Mtb) are useful to identify new molecules that could lead to future tuberculosis (TB) drugs. Our approach uses information derived from the TBCyc pathway and genome database, the Collaborative Drug Discovery TB database combined with 3D pharmacophores and dual event Bayesian models of whole-cell activity and lack of cytotoxicity. We have prioritized a large number of molecules that may act as mimics of substrates and metabolites in the TB metabolome. We computationally searched over 200,000 commercial molecules using 66 pharmacophores based on substrates and metabolites from Mtb and further filtering with Bayesian models. We ultimately tested 110 compounds in vitro that resulted in two compounds of interest, BAS 04912643 and BAS 00623753 (MIC of 2.5 and 5 µg/mL, respectively). These molecules were used as a starting point for hit-to-lead optimization. The most promising class proved to be the quinoxaline di-N-oxides, evidenced by transcriptional profiling to induce mRNA level perturbations most closely resembling known protonophores. One of these, SRI58 exhibited an MIC = 1.25 µg/mL versus Mtb and a CC50 in Vero cells of >40 µg/mL, while featuring fair Caco-2 A-B permeability (2.3 x 10-6 cm/s), kinetic solubility (125 µM at pH 7.4 in PBS) and mouse metabolic stability (63.6% remaining after 1 h incubation with mouse liver microsomes). Despite demonstration of how a combined bioinformatics/cheminformatics approach afforded a small molecule with promising in vitro profiles, we found that SRI58 did not exhibit quantifiable blood levels in mice.