High-throughput screening and Bayesian machine learning for copper-dependent inhibitors of Staphylococcus aureus

Identification of a copper-responsive small molecule inhibitor of uropathogenic Escherichia coli

Urinary tract infections (UTIs) are a major global health problem and are caused predominantly by uropathogenic Escherichia coli (UPEC). UTIs are a leading cause of prescription antimicrobial use. Incessant increase in antimicrobial resistance in UPEC and other uropathogens poses a serious threat to the current treatment practices. Copper is an effector of nutritional immunity that impedes the growth of pathogens during infection. We hypothesized that copper would augment the toxicity of select small molecules against bacterial pathogens. We conducted a small molecule screening campaign with a library of 51,098 molecules to detect hits that inhibit a UPEC ΔtolC mutant in a copper-dependent manner. A molecule, denoted as E. coli inhibitor or ECIN, was identified as a copper-responsive inhibitor of wild-type UPEC strains. Our gene expression and metal content analysis results demonstrate that ECIN works in concert with copper to exacerbate Cu toxicity in UPEC. ECIN has a broad spectrum of activity against pathogens of medical and veterinary significance including Acinetobacter baumannii, Pseudomonas aeruginosa, and methicillin-resistant Staphylococcus aureus. Subinhibitory levels of ECIN eliminate UPEC biofilm formation. Transcriptome analysis of UPEC treated with ECIN reveals induction of multiple stress response systems. Furthermore, we demonstrate that L-cysteine rescues the growth of UPEC exposed to ECIN. In summary, we report the identification and characterization of a novel copper-responsive small molecule inhibitor of UPEC.IMPORTANCEUrinary tract infection (UTI) is a ubiquitous infectious condition affecting millions of people annually. Uropathogenic Escherichia coli (UPEC) is the predominant etiological agent of UTI. However, UTIs are becoming increasingly difficult to resolve with antimicrobials due to increased antimicrobial resistance in UPEC and other uropathogens. Here, we report the identification and characterization of a novel copper-responsive small molecule inhibitor of UPEC. In addition to E. coli, this small molecule also inhibits pathogens of medical and veterinary significance including Acinetobacter baumannii, Pseudomonas aeruginosa, and methicillin-resistant Staphylococcus aureus.

A high-throughput response to the SARS-CoV-2 pandemic

Four years after the beginning of the COVID-19 pandemic, it is important to reflect on the events that have occurred during that time and the knowledge that has been gained. The response to the pandemic was rapid and highly resourced; it was also built upon a foundation of decades of federally funded basic and applied research. Laboratories in government, pharmaceutical, academic, and non-profit institutions all played roles in advancing pre-2020 discoveries to produce clinical treatments. This perspective provides a summary of how the development of high-throughput screening methods in a biosafety level 3 (BSL-3) environment at Southern Research Institute (SR) contributed to pandemic response efforts. The challenges encountered are described, including those of a technical nature as well as those of working under the pressures of an unpredictable virus and pandemic.

High-Throughput Phenotypic Screening and Machine Learning Methods Enabled the Selection of Broad-Spectrum Low-Toxicity Antitrypanosomatidic Agents

Broad-spectrum anti-infective chemotherapy agents with activity against Trypanosomes, Leishmania, and Mycobacterium tuberculosis species were identified from a high-throughput phenotypic screening program of the 456 compounds belonging to the Ty-Box, an in-house industry database. Compound characterization using machine learning approaches enabled the identification and synthesis of 44 compounds with broad-spectrum antiparasitic activity and minimal toxicity against Trypanosoma brucei, Leishmania Infantum, and Trypanosoma cruzi. In vitro studies confirmed the predictive models identified in compound 40 which emerged as a new lead, featured by an innovative N-(5-pyrimidinyl)benzenesulfonamide scaffold and promising low micromolar activity against two parasites and low toxicity. Given the volume and complexity of data generated by the diverse high-throughput screening assays performed on the compounds of the Ty-Box library, the chemoinformatic and machine learning tools enabled the selection of compounds eligible for further evaluation of their biological and toxicological activities and aided in the decision-making process toward the design and optimization of the identified lead.

Repurposing sunscreen as an antibiotic: zinc-activated avobenzone inhibits methicillin-resistant Staphylococcus aureus

Methicillin-resistant Staphylococcus aureus (MRSA) is a major healthcare concern with associated healthcare costs reaching over ${\$}$1 billion in a single year in the USA. Antibiotic resistance in S. aureus is now observed against last line of defense antibiotics, such as vancomycin, linezolid, and daptomycin. Unfortunately, high throughput drug discovery approaches to identify new antibiotics effective against MRSA have not resulted in much tangible success over the last decades. Previously, we demonstrated the feasibility of an alternative drug discovery approach, the identification of metallo-antibiotics, compounds that gain antibacterial activity only after binding to a transition metal ion and as such are unlikely to be detected in standard drug screens. We now report that avobenzone, the primary active ingredient of most sunscreens, can be activated by zinc to become a potent antibacterial compound against MRSA. Zinc-activated avobenzone (AVB-Zn) potently inhibited a series of clinical MRSA isolates [minimal inhibitory concentration (MIC): 0.62-2.5 µM], without pre-existing resistance and activity without zinc (MIC: >10 µM). AVB-Zn was also active against clinical MRSA isolates that were resistant against the commonly used zinc-salt antibiotic bacitracin. We found AVB-Zn exerted no cytotoxicity on human cell lines and primary cells. Last, we demonstrate AVB-Zn can be deployed therapeutically as lotion preparations, which showed efficacy in a mouse wound model of MRSA infection. AVB-Zn thus demonstrates Zn-activated metallo-antibiotics are a promising avenue for future drug discovery.

High-Throughput Screening of Natural Product and Synthetic Molecule Libraries for Antibacterial Drug Discovery

Due to the continued emergence of resistance and a lack of new and promising antibiotics, bacterial infection has become a major public threat. High-throughput screening (HTS) allows rapid screening of a large collection of molecules for bioactivity testing and holds promise in antibacterial drug discovery. More than 50% of the antibiotics that are currently available on the market are derived from natural products. However, with the easily discoverable antibiotics being found, finding new antibiotics from natural sources has seen limited success. Finding new natural sources for antibacterial activity testing has also proven to be challenging. In addition to exploring new sources of natural products and synthetic biology, omics technology helped to study the biosynthetic machinery of existing natural sources enabling the construction of unnatural synthesizers of bioactive molecules and the identification of molecular targets of antibacterial agents. On the other hand, newer and smarter strategies have been continuously pursued to screen synthetic molecule libraries for new antibiotics and new druggable targets. Biomimetic conditions are explored to mimic the real infection model to better study the ligand-target interaction to enable the designing of more effective antibacterial drugs. This narrative review describes various traditional and contemporaneous approaches of high-throughput screening of natural products and synthetic molecule libraries for antibacterial drug discovery. It further discusses critical factors for HTS assay design, makes a general recommendation, and discusses possible alternatives to traditional HTS of natural products and synthetic molecule libraries for antibacterial drug discovery.

Identification of Novel Arachidonic Acid 15-Lipoxygenase Inhibitors Based on the Bayesian Classifier Model and Computer-Aided High-Throughput Virtual Screening

Ferroptosis is an iron-dependent lipid peroxidative form of cell death that is distinct from apoptosis and necrosis. ALOX15, also known as arachidonic acid 15-lipoxygenase, promotes ferroptosis by converting intracellular unsaturated lipids into oxidized lipid intermediates and is an important ferroptosis target. In this study, a naive Bayesian machine learning classifier with a structure-based, high-throughput screening approach and a molecular docking program were combined to screen for three compounds with excellent target-binding potential. In the absorption, distribution, metabolism, excretion, and toxicity characterization, three candidate molecules were predicted to exhibit drug-like properties. The subsequent molecular dynamics simulations confirmed their stable binding to the targets. The findings indicated that the compounds exhibited excellent potential ALOX15 inhibitor capacity, thereby providing novel candidates for the treatment of inflammatory ischemia-related diseases caused by ferroptosis.

Machine Learning for Discovery of New ADORA Modulators

Adenosine (ADO) is an extracellular signaling molecule generated locally under conditions that produce ischemia, hypoxia, or inflammation. It is involved in modulating a range of physiological functions throughout the brain and periphery through the membrane-bound G protein-coupled receptors, called adenosine receptors (ARs) A₁AR, A_2AAR, A_2BAR, and A₃AR. These are therefore important targets for neurological, cardiovascular, inflammatory, and autoimmune diseases and are the subject of drug development directed toward the cyclic adenosine monophosphate and other signaling pathways. Initially using public data for A₁AR agonists we generated and validated a Bayesian machine learning model (Receiver Operator Characteristic of 0.87) that we used to identify molecules for testing. Three selected molecules, crisaborole, febuxostat and paroxetine, showed initial activity in vitro using the HEK293 A₁AR Nomad cell line. However, radioligand binding, β-arrestin assay and calcium influx assay did not confirm this A₁AR activity. Nevertheless, several other AR activities were identified. Febuxostat and paroxetine both inhibited orthosteric radioligand binding in the µM range for A_2AAR and A₃AR. In HEK293 cells expressing the human A_2AAR, stimulation of cAMP was observed for crisaborole (EC₅₀ 2.8 µM) and paroxetine (EC₅₀ 14 µM), but not for febuxostat. Crisaborole also increased cAMP accumulation in A_2BAR-expressing HEK293 cells, but it was weaker than at the A_2AAR. At the human A₃AR, paroxetine did not show any agonist activity at 100 µM, although it displayed binding with a K_i value of 14.5 µM, suggesting antagonist activity. We have now identified novel modulators of A_2AAR, A_2BAR and A₃AR subtypes that are clinically used for other therapeutic indications, and which are structurally distinct from previously reported tool compounds or drugs.

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 IC₅₀s 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 IC₅₀ 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 EC₅₀ 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.

Mycobacterium abscessus drug discovery using machine learning

The prevalence of infections by nontuberculous mycobacteria is increasing, having surpassed tuberculosis in the United States and much of the developed world. Nontuberculous mycobacteria occur naturally in the environment and are a significant problem for patients with underlying lung diseases such as bronchiectasis, chronic obstructive pulmonary disease, and cystic fibrosis. Current treatment regimens are lengthy, complicated, toxic and they are often unsuccessful as seen by disease recurrence. Mycobacterium abscessus is one of the most commonly encountered organisms in nontuberculous mycobacteria disease and it is the most difficult to eradicate. There is currently no systematically proven regimen that is effective for treating M. abscessus infections. Our approach to drug discovery integrates machine learning, medicinal chemistry and in vitro testing and has been previously applied to Mycobacterium tuberculosis. We have now identified several novel 1-(phenylsulfonyl)-1H-benzimidazol-2-amines that have weak activity on M. abscessus in vitro but may represent a starting point for future further medicinal chemistry optimization. We also address limitations still to be overcome with the machine learning approach for M. abscessus.

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 (K_d 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.

Molecular and Cellular Function of Transcription Factor 4 in Pitt-Hopkins Syndrome

Transcription factor 4 (TCF4, also known as ITF2 or E2-2) is a type I basic helix-loop-helix transcription factor. Autosomal dominant mutations in TCF4 cause Pitt-Hopkins syndrome (PTHS), a rare syndromic form of autism spectrum disorder. In this review, we provide an update on the progress regarding our understanding of TCF4 function at the molecular, cellular, physiological, and behavioral levels with a focus on phenotypes and therapeutic interventions. We examine upstream and downstream regulatory networks associated with TCF4 and discuss a range of in vitro and in vivo data with the aim of understanding emerging TCF4-specific mechanisms relevant for disease pathophysiology. In conclusion, we provide comments about exciting future avenues of research that may provide insights into potential new therapeutic targets for PTHS.

Comparing the Pfizer Central Nervous System Multiparameter Optimization Calculator and a BBB Machine Learning Model

The ability to calculate whether small molecules will cross the blood-brain barrier (BBB) is an important task for companies working in neuroscience drug discovery. For a decade, scientists have relied on relatively simplistic rules such as Pfizer's central nervous system multiparameter optimization models (CNS-MPO) for guidance during the drug selection process. In parallel, there has been a continued development of more sophisticated machine learning models that utilize different molecular descriptors and algorithms; however, these models represent a "black box" and are generally less interpretable. In both cases, these methods predict the ability of small molecules to cross the BBB using the molecular structure information on its own without in vitro or in vivo data. We describe here the implementation of two versions of Pfizer's algorithm (Pf-MPO.v1 and Pf-MPO.v2) and compare it with a Bayesian machine learning model of BBB penetration trained on a data set of 2296 active and inactive compounds using extended connectivity fingerprint descriptors. The predictive ability of these approaches was compared with 40 known CNS active drugs initially used by Pfizer as their positive set for validation of the Pf-MPO.v1 score. 37/40 (92.5%) compounds were predicted as active by the Bayesian model, while only 30/40 (75%) received a desirable Pf-MPO.v1 score ≥4 and 33/40 (82.5%) received a desirable Pf-MPO.v2 score ≥4, suggesting the Bayesian model is more accurate than MPO algorithms. This also indicates machine learning models are more flexible and have better predictive power for BBB penetration than simple rule sets that require multiple, accurate descriptor calculations. Our machine learning model statistics are comparable to recent published studies. We describe the implications of these findings and how machine learning may have a role alongside more interpretable methods.

Discovery of 5-Nitro-6-thiocyanatopyrimidines as Inhibitors of Cryptococcus neoformans and Cryptococcus gattii

Opportunistic infections from pathogenic fungi present a major challenge to healthcare because of a very limited arsenal of antifungal drugs, an increasing population of immunosuppressed patients, and increased prevalence of resistant clinical strains due to overuse of the few available antifungals. Cryptococcal meningitis is a life-threatening opportunistic fungal infection caused by one of two species in the Cryptococcus genus, Cryptococcus neoformans and Cryptococcus gattii. Eighty percent of cryptococcosis diseases are caused by C. neoformans that is endemic in the environment. The standard of care is limited to old antifungals, and under a high standard of care, mortality remains between 10 and 30%. We have identified a series of 5-nitro-6-thiocyanatopyrimidine antifungal drug candidates using in vitro and computational machine learning approaches. These compounds can inhibit C. neoformans growth at submicromolar levels, are effective against fluconazole-resistant C. neoformans and a clinical strain of C. gattii, and are not antagonistic with currently approved antifungals.

Flavonoids from Pterogyne nitens as Zika virus NS2B-NS3 protease inhibitors

Although the widespread epidemic of Zika virus (ZIKV) and its neurological complications are well-known there are still no approved drugs available to treat this arboviral disease or vaccine to prevent the infection. Flavonoids from Pterogyne nitens have already demonstrated anti-flavivirus activity, although their target is unknown. In this study, we virtually screened an in-house database of 150 natural and semi-synthetic compounds against ZIKV NS2B-NS3 protease (NS2B-NS3p) using docking-based virtual screening, as part of the OpenZika project. As a result, we prioritized three flavonoids from P. nitens, quercetin, rutin and pedalitin, for experimental evaluation. We also used machine learning models, built with Assay Central® software, for predicting the activity and toxicity of these flavonoids. Biophysical and enzymatic assays generally agreed with the in silico predictions, confirming that the flavonoids inhibited ZIKV protease. The most promising hit, pedalitin, inhibited ZIKV NS2B-NS3p with an IC50 of 5 μM. In cell-based assays, pedalitin displayed significant activity at 250 and 500 µM, with slight toxicity in Vero cells. The results presented here demonstrate the potential of pedalitin as a candidate for hit-to-lead (H2L) optimization studies towards the discovery of antiviral drug candidates to treat ZIKV infections.

Pyrazolyl Thioureas and Carbothioamides with an NNSN Motif against MSSA and MRSA

A novel series of copper-activatable drugs intended for use against methicillin-sensitive Staphylococcus aureus (MSSA) and methicillin-resistant S. aureus (MRSA) were synthesized, characterized, and tested against the MSSA strain Newman and the MRSA Lac strain (a USA300 strain), respectively. These drugs feature an NNSN structural motif, which enables the binding of copper. In the absence of copper, no activity against MSSA and MRSA at realistic drug concentrations was observed. Although none of the novel drug candidates exhibits a stereocenter, sub-micromolar activities against SA Newman and micromolar activities against SA Lac were observed in the presence, but not in the absence, of bioavailable copper. Copper influx is a component of cellular response to bacterial infections, which is often described as nutritional immunity.

A Machine Learning Strategy for Drug Discovery Identifies Anti-Schistosomal Small Molecules

Schistosomiasis is a chronic and painful disease of poverty caused by the flatworm parasite Schistosoma. Drug discovery for antischistosomal compounds predominantly employs in vitro whole organism (phenotypic) screens against two developmental stages of Schistosoma mansoni, post-infective larvae (somules) and adults. We generated two rule books and associated scoring systems to normalize 3898 phenotypic data points to enable machine learning. The data were used to generate eight Bayesian machine learning models with the Assay Central software according to parasite’s developmental stage and experimental time point (≤24, 48, 72, and >72 h). The models helped predict 56 active and nonactive compounds from commercial compound libraries for testing. When these were screened against S. mansoni in vitro, the prediction accuracy for active and inactives was 61% and 56% for somules and adults, respectively; also, hit rates were 48% and 34%, respectively, far exceeding the typical 1–2% hit rate for traditional high throughput screens.

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 IC₅₀'s of 14.4 nM and 64.4 nM, respectively, but does not inhibit the closely related BuChE IC₅₀ > 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.

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.

Machine Learning for Discovery of GSK3β Inhibitors

Alzheimer's disease (AD) is the most common cause of dementia, affecting approximately 35 million people worldwide. The current treatment options for people with AD consist of drugs designed to slow the rate of decline in memory and cognition, but these treatments are not curative, and patients eventually suffer complete cognitive injury. With the substantial amounts of published data on targets for this disease, we proposed that machine learning software could be used to find novel small-molecule treatments that can supplement the AD drugs currently on the market. In order to do this, we used publicly available data in ChEMBL to build and validate Bayesian machine learning models for AD target proteins. The first AD target that we have addressed with this method is the serine-threonine kinase glycogen synthase kinase 3 beta (GSK3β), which is a proline-directed serine-threonine kinase that phosphorylates the microtubule-stabilizing protein tau. This phosphorylation prompts tau to dissociate from the microtubule and form insoluble oligomers called paired helical filaments, which are one of the components of the neurofibrillary tangles found in AD brains. Using our Bayesian machine learning model for GSK3β consisting of 2368 molecules, this model produced a five-fold cross validation ROC of 0.905. This model was also used for virtual screening of large libraries of FDA-approved drugs and clinical candidates. Subsequent testing of selected compounds revealed a selective small-molecule inhibitor, ruboxistaurin, with activity against GSK3β (avg IC50 = 97.3 nM) and GSK3α (IC50 = 695.9 nM). Several other structurally diverse inhibitors were also identified. We are now applying this machine learning approach to additional AD targets to identify approved drugs or clinical trial candidates that can be repurposed as AD therapeutics. This represents a viable approach to accelerate drug discovery and do so at a fraction of the cost of traditional high throughput screening.

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.

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.

Synergistic drug combinations and machine learning for drug repurposing in chordoma

Chordoma is a devastating rare cancer that affects one in a million people. With a mean-survival of just 6 years and no approved medicines, the primary treatments are surgery and radiation. In order to speed new medicines to chordoma patients, a drug repurposing strategy represents an attractive approach. Drugs that have already advanced through human clinical safety trials have the potential to be approved more quickly than de novo discovered medicines on new targets. We have taken two strategies to enable this: (1) generated and validated machine learning models of chordoma inhibition and screened compounds of interest in vitro. (2) Tested combinations of approved kinase inhibitors already being individually evaluated for chordoma. Several published studies of compounds screened against chordoma cell lines were used to generate Bayesian Machine learning models which were then used to score compounds selected from the NIH NCATS industry-provided assets. Out of these compounds, the mTOR inhibitor AZD2014, was the most potent against chordoma cell lines (IC50 0.35 µM U-CH1 and 0.61 µM U-CH2). Several studies have shown the importance of the mTOR signaling pathway in chordoma and suggest it as a promising avenue for targeted therapy. Additionally, two currently FDA approved drugs, afatinib and palbociclib (EGFR and CDK4/6 inhibitors, respectively) demonstrated synergy in vitro (CI50 = 0.43) while AZD2014 and afatanib also showed synergy (CI50 = 0.41) against a chordoma cell in vitro. These findings may be of interest clinically, and this in vitro- and in silico approach could also be applied to other rare cancers.

Repurposing the Dihydropyridine Calcium Channel Inhibitor Nicardipine as a Nav1.8 Inhibitor In Vivo for Pitt Hopkins Syndrome

PURPOSE

Individuals with the rare genetic disorder Pitt Hopkins Syndrome (PTHS) do not have sufficient expression of the transcription factor 4 (TCF4) which is located on chromosome 18. TCF4 is a basic helix-loop-helix E protein that is critical for the normal development of the nervous system and the brain in humans. PTHS patients lacking sufficient TCF4 frequently display gastrointestinal issues, intellectual disability and breathing problems. PTHS patients also commonly do not speak and display distinctive facial features and seizures. Recent research has proposed that decreased TCF4 expression can lead to the increased translation of the sodium channel Na_v1.8. This in turn results in increased after-hyperpolarization as well as altered firing properties. We have recently identified through a drug repurposing screen an FDA approved dihydropyridine calcium antagonist nicardipine used to treat angina, which inhibited Na_v1.8.

METHODS

We have now performed behavioral testing in groups of 10 male Tcf4(± ) PTHS mice dosing by oral gavage at 3 mg/kg once a day for 3 weeks using standard methods to assess sociability, nesting, fear conditioning, self-grooming, open field and test of force.

RESULTS

Nicardipine returned this spectrum of behavioral deficits in the Tcf4(± ) PTHS mouse model to WT levels and resulted in statistically significant results.

CONCLUSIONS

These in vivo results in the well characterized Tcf4(± ) PTHS mice may suggest the potential to test this already approved drug further in a clinical study with PTHS patients or suggest the potential for use off label under compassionate use with their physician.

A copper-dependent compound restores ampicillin sensitivity in multidrug-resistant Staphylococcus aureus

Multi-drug resistant Staphylococcus aureus, including methicillin-resistant S. aureus (MRSA), has become a worldwide, major health care problem. While initially restricted to clinical settings, drug resistant S. aureus is now one of the key causative agents of community-acquired infections. We have previously demonstrated that copper dependent inhibitors (CDIs), a class of antibiotics that are only active in the presence of copper ions, are effective bactericidal agents against MRSA. A second-generation CDI, APT-6K, exerted bactericidal activity at nanomolar concentrations. At sub-bactericidal concentrations, it effectively synergized with ampicillin to reverse drug resistance in multiple MRSA strains. APT-6K had a favorable therapeutic index when tested on eukaryotic cells (TI: > 30) and, unlike some previously reported CDIs, did not affect mitochondrial activity. These results further establish inhibitors that are activated by the binding of transition metal ions as a promising class of antibiotics, and for the first time, describe their ability to reverse existing drug resistance against clinically relevant antibiotics.

Repurposing Approved Drugs as Inhibitors of Kv7.1 and Nav1.8 to Treat Pitt Hopkins Syndrome

PURPOSE

Pitt Hopkins Syndrome (PTHS) is a rare genetic disorder caused by mutations of a specific gene, transcription factor 4 (TCF4), located on chromosome 18. PTHS results in individuals that have moderate to severe intellectual disability, with most exhibiting psychomotor delay. PTHS also exhibits features of autistic spectrum disorders, which are characterized by the impaired ability to communicate and socialize. PTHS is comorbid with a higher prevalence of epileptic seizures which can be present from birth or which commonly develop in childhood. Attenuated or absent TCF4 expression results in increased translation of peripheral ion channels K_v7.1 and Na_v1.8 which triggers an increase in after-hyperpolarization and altered firing properties.

METHODS

We now describe a high throughput screen (HTS) of 1280 approved drugs and machine learning models developed from this data. The ion channels were expressed in either CHO (K_V7.1) or HEK293 (Na_v1.8) cells and the HTS used either ⁸⁶Rb⁺ efflux (K_V7.1) or a FLIPR assay (Na_v1.8).

RESULTS

The HTS delivered 55 inhibitors of K_v7.1 (4.2% hit rate) and 93 inhibitors of Na_v1.8 (7.2% hit rate) at a screening concentration of 10 μM. These datasets also enabled us to generate and validate Bayesian machine learning models for these ion channels. We also describe a structure activity relationship for several dihydropyridine compounds as inhibitors of Na_v1.8.

CONCLUSIONS

This work could lead to the potential repurposing of nicardipine or other dihydropyridine calcium channel antagonists as potential treatments for PTHS acting via Na_v1.8, as there are currently no approved treatments for this rare disorder.

Emerging Opportunities To Manipulate Metal Trafficking for Therapeutic Benefit

The indispensable requirement for metals in life processes has led to the evolution of sophisticated mechanisms that allow organisms to maintain dynamic equilibria of these ions. This dynamic control of the level, speciation, and availability of a variety of metal ions allows organisms to sustain biological processes while avoiding toxicity. When functioning properly, these mechanisms allow cells to return to their metal homeostatic set points following shifts in the metal availability or other stressors. These periods of transition, when cells are in a state of flux in which they work to regain homeostasis, present windows of opportunity to pharmacologically manipulate targets associated with metal-trafficking pathways in ways that could either facilitate a return to homeostasis and the recovery of cellular function or further push cells outside of homeostasis and into cellular distress. The purpose of this Viewpoint is to highlight emerging opportunities for chemists and chemical biologists to develop compounds to manipulate metal-trafficking processes for therapeutic benefit.