Advancing Biological Understanding and Therapeutics Discovery with Small-Molecule Probes.

Small-molecule probes can illuminate biological processes and aid in the assessment of emerging therapeutic targets by perturbing biological systems in a manner distinct from other experimental approaches. Despite the tremendous promise of chemical tools for investigating biology and disease, small-molecule probes were unavailable for most targets and pathways as recently as a decade ago. In 2005, the NIH launched the decade-long Molecular Libraries Program with the intent of innovating in and broadening access to small-molecule science. This Perspective describes how novel small-molecule probes identified through the program are enabling the exploration of biological pathways and therapeutic hypotheses not otherwise testable. These experiences illustrate how small-molecule probes can help bridge the chasm between biological research and the development of medicines but also highlight the need to innovate the science of therapeutic discovery.

High-Throughput Screening to Advance In Vitro Toxicology: Accomplishments, Challenges, and Future Directions.

Traditionally, chemical toxicity is determined by in vivo animal studies, which are low throughput, expensive, and sometimes fail to predict compound toxicity in humans. Due to the increasing number of chemicals in use and the high rate of drug candidate failure due to toxicity, it is imperative to develop in vitro, high-throughput screening methods to determine toxicity. The Tox21 program, a unique research consortium of federal public health agencies, was established to address and identify toxicity concerns in a high-throughput, concentration-responsive manner using a battery of in vitro assays. In this article, we review the advancements in high-throughput robotic screening methodology and informatics processes to enable the generation of toxicological data, and their impact on the field; further, we discuss the future of assessing environmental toxicity utilizing efficient and scalable methods that better represent the corresponding biological and toxicodynamic processes in humans.

Discovering selective antiferroptotic inhibitors of the 15LOX/PEBP1 complex noninterfering with biosynthesis of lipid mediators.

Programmed ferroptotic death eliminates cells in all major organs and tissues with imbalanced redox metabolism due to overwhelming iron-catalyzed lipid peroxidation under insufficient control by thiols (Glutathione (GSH)). Ferroptosis has been associated with the pathogenesis of major chronic degenerative diseases and acute injuries of the brain, cardiovascular system, liver, kidneys, and other organs, and its manipulation offers a promising new strategy for anticancer therapy. This explains the high interest in designing new small-molecule-specific inhibitors against ferroptosis. Given the role of 15-lipoxygenase (15LOX) association with phosphatidylethanolamine (PE)-binding protein 1 (PEBP1) in initiating ferroptosis-specific peroxidation of polyunsaturated PE, we propose a strategy of discovering antiferroptotic agents as inhibitors of the 15LOX/PEBP1 catalytic complex rather than 15LOX alone. Here we designed, synthesized, and tested a customized library of 26 compounds using biochemical, molecular, and cell biology models along with redox lipidomic and computational analyses. We selected two lead compounds, FerroLOXIN-1 and 2, which effectively suppressed ferroptosis in vitro and in vivo without affecting the biosynthesis of pro-/anti-inflammatory lipid mediators in vivo. The effectiveness of these lead compounds is not due to radical scavenging or iron-chelation but results from their specific mechanisms of interaction with the 15LOX-2/PEBP1 complex, which either alters the binding pose of the substrate [eicosatetraenoyl-PE (ETE-PE)] in a nonproductive way or blocks the predominant oxygen channel thus preventing the catalysis of ETE-PE peroxidation. Our successful strategy may be adapted to the design of additional chemical libraries to reveal new ferroptosis-targeting therapeutic modalities.

A high dose KRP203 induces cytoplasmic vacuoles associated with altered phosphoinositide segregation and endosome expansion.

In animal cells, vacuoles are absent, but can be induced by diseases and drugs. While phosphoinositides are critical for membrane trafficking, their role in the formation of these vacuoles remains unclear. The immunosuppressive KRP203/Mocravimod, which antagonizes sphingosine-1-phosphate receptors, has been identified as having novel multimodal activity against phosphoinositide kinases. However, the impact of this novel KRP203 activity is unknown. Here, we show that KRP203 disrupts the spatial organization of phosphoinositides and induces extensive vacuolization in tumor cells and immortalized fibroblasts. The KRP203-induced vacuoles are primarily from endosomes, and augmented by inhibition of PIKFYVE and VPS34. Conversely, overexpression of PTEN decreased KRP203-induced vacuole formation. Furthermore, V-ATPase inhibition completely blunted KRP203-induced vacuolization, pointing to a critical requirement of the endosomal maturation process. Importantly, nearly a half of KRP203-induced vacuoles are significantly decorated with PI4P, a phosphoinositide typically enriched at the plasma membrane and Golgi. These results suggest a model that noncanonical spatial reorganization of phosphoinositides by KRP203 alters the endosomal maturation process, leading to vacuolization. Taken together, this study reveals a previously unrecognized bioactivity of KRP203 as a vacuole-inducing agent and its unique mechanism of phosphoinositide modulation, providing a new insight of phosphoinositide regulation into vacuolization-associated diseases and their molecular pathologies.

A versatile polypharmacology platform promotes cytoprotection and viability of human pluripotent and differentiated cells.

Human pluripotent stem cells (hPSCs) are capable of extensive self-renewal yet remain highly sensitive to environmental perturbations in vitro, posing challenges to their therapeutic use. There is an urgent need to advance strategies that ensure safe and robust long-term growth and functional differentiation of these cells. Here, we deployed high-throughput screening strategies to identify a small-molecule cocktail that improves viability of hPSCs and their differentiated progeny. The combination of chroman 1, emricasan, polyamines, and trans-ISRIB (CEPT) enhanced cell survival of genetically stable hPSCs by simultaneously blocking several stress mechanisms that otherwise compromise cell structure and function. CEPT provided strong improvements for several key applications in stem-cell research, including routine cell passaging, cryopreservation of pluripotent and differentiated cells, embryoid body (EB) and organoid formation, single-cell cloning, and genome editing. Thus, CEPT represents a unique poly-pharmacological strategy for comprehensive cytoprotection, providing a rationale for efficient and safe utilization of hPSCs.

A quantitative high-throughput screen identifies compounds that lower expression of the SCA2-and ALS-associated gene ATXN2.

CAG repeat expansions in the ATXN2 (ataxin-2) gene can cause the autosomal dominant disorder spinocerebellar ataxia type 2 (SCA2) as well as increase the risk of ALS. Abnormal molecular, motor, and neurophysiological phenotypes in SCA2 mouse models are normalized by lowering ATXN2 transcription, and reduction of nonmutant Atxn2 expression has been shown to increase the life span of mice overexpressing the TDP-43 (transactive response DNA-binding protein 43 kDa) ALS protein, demonstrating the potential benefits of targeting ATXN2 transcription in humans. Here, we describe a quantitative high-throughput screen to identify compounds that lower ATXN2 transcription. We screened 428,759 compounds in a multiplexed assay using an ATXN2-luciferase reporter in human embryonic kidney 293 (HEK-293) cells and identified a diverse set of compounds capable of lowering ATXN2 transcription. We observed dose-dependent reductions of endogenous ATXN2 in HEK-293 cells treated with procillaridin A, 17-dimethylaminoethylamino-17-demethoxygeldanamycin (17-DMAG), and heat shock protein 990 (HSP990), known inhibitors of HSP90 and Na+/K+-ATPases. Furthermore, HEK-293 cells expressing polyglutamine-expanded ATXN2-Q58 treated with 17-DMAG had minimally detectable ATXN2, as well as normalized markers of autophagy and endoplasmic reticulum stress, including STAU1 (Staufen 1), molecular target of rapamycin, p62, LC3-II (microtubule-associated protein 1A/1B-light chain 3II), CHOP (C/EBP homologous protein), and phospho-eIF2α (eukaryotic initiation factor 2α). Finally, bacterial artificial chromosome ATXN2-Q22 mice treated with 17-DMAG or HSP990 exhibited highly reduced ATXN2 protein abundance in the cerebellum. Taken together, our study demonstrates inhibition of HSP90 or Na+/K+-ATPases as potentially effective therapeutic strategies for treating SCA2 and ALS.

Novel Transmission-Blocking Antimalarials Identified by High-Throughput Screening of Plasmodium berghei Ookluc.

Safe and effective malaria transmission-blocking chemotherapeutics would allow a community-level approach to malaria control and eradication efforts by targeting the mosquito sexual stage of the parasite life cycle. However, only a single drug, primaquine, is currently approved for use in reducing transmission, and drug toxicity limits its widespread implementation. To address this limitation in antimalarial chemotherapeutics, we used a recently developed transgenic Plasmodium berghei line, Ookluc, to perform a series of high-throughput in vitro screens for compounds that inhibit parasite fertilization, the initial step of parasite development within the mosquito. Screens of antimalarial compounds, approved drug collections, and drug-like molecule libraries identified 185 compounds that inhibit parasite maturation to the zygote form. Seven compounds were further characterized to block gametocyte activation or to be cytotoxic to formed zygotes. These were further validated in mosquito membrane-feeding assays using Plasmodium falciparum and P. vivax. This work demonstrates that high-throughput screens using the Ookluc line can identify compounds that are active against the two most relevant human Plasmodium species and provides a list of compounds that can be explored for the development of new antimalarials to block transmission.

Multimodal action of KRP203 on phosphoinositide kinases in vitro and in cells.

Increased phosphoinositide signaling is commonly associated with cancers. While "one-drug one-target" has been a major drug discovery strategy for cancer therapy, a "one-drug multi-targets" approach for phosphoinositide enzymes has the potential to offer a new therapeutic approach. In this study, we sought a new way to target phosphoinositides metabolism. Using a high-throughput phosphatidylinositol 5-phosphate 4-kinase-alpha (PI5P4Kα) assay, we have identified that the immunosuppressor KRP203/Mocravimod induces a significant perturbation in phosphoinositide metabolism in U87MG glioblastoma cells. Despite high sequence similarity of PI5P4K and PI4K isozymes, in vitro kinase assays showed that KRP203 activates some (e.g., PI5P4Kα, PI4KIIß) while inhibiting other phosphoinositide kinases (e.g., PI5P4Kß, γ, PI4KIIα, class I PI3K-p110α, δ, γ). Furthermore, KRP203 enhances PI3P5K/PIKFYVE's substrate selectivity for phosphatidylinositol (PI) while preserving its selectivity for PI(3)P. At cellular levels, 3 h of KRP203 treatment induces a prominent increase of PI(3)P and moderate increase of PI(5)P, PI(3,5)P2, and PI(3,4,5)P3 levels in U87MG cells. Collectively, the finding of multimodal activity of KRP203 towards multi-phosphoinositide kinases may open a novel basis to modulate cellular processes, potentially leading to more effective treatments for diseases associated with phosphoinositide signaling pathways.

In vitro profiling of pesticides within the Tox21 10K compound library for bioactivity and potential toxicity.

Pesticides include a diverse class of toxic chemicals, often having numerous modes of actions when used in agriculture against targeted organisms to control insect infestation, halt unwanted vegetation, and prevent the spread of disease. In this study, the in vitro assay activity of pesticides within the Tox21 10K compound library were examined. The assays in which pesticides showed significantly more activities than non-pesticide chemicals revealed potential targets and mechanisms of action for pesticides. Furthermore, pesticides that showed promiscuous activity against many targets and cytotoxicity were identified, which warrant further toxicological evaluation. Several pesticides were shown to require metabolic activation, demonstrating the importance of introducing metabolic capacity to in vitro assays. Overall, the activity profiles of pesticides highlighted in this study can contribute to the knowledge gaps surrounding pesticide mechanisms and to the better understanding of the on- and off-target organismal effects of pesticides.

Predictive Models for Human Cytochrome P450 3A7 Selective Inhibitors and Substrates.

Inappropriate use of prescription drugs is potentially more harmful in fetuses/neonates than in adults. Cytochrome P450 (CYP) 3A subfamily undergoes developmental changes in expression, such as a transition from CYP3A7 to CYP3A4 shortly after birth, which provides a potential way to distinguish medication effects on fetuses/neonates and adults. The purpose of this study was to build first-in-class predictive models for both inhibitors and substrates of CYP3A7/CYP3A4 using chemical structure analysis. Three metrics were used to evaluate model performance: area under the receiver operating characteristic curve (AUC-ROC), balanced accuracy (BA), and Matthews correlation coefficient (MCC). The performance varied for each CYP3A7/CYP3A4 inhibitor/substrate model depending on the data set type, model type, rebalancing method, and specific feature set. For the active inhibitor/substrate data set, the optimal models achieved AUC-ROC values ranging from 0.77 ± 0.01 to 0.84 ± 0.01. For the selective inhibitor/substrate data set, the optimal models achieved AUC-ROC values ranging from 0.72 ± 0.02 to 0.79 ± 0.04. The predictive power of the optimal models was validated by compounds with known potencies as CYP3A7/CYP3A4 inhibitors or substrates. In addition, we identified structural features significant for CYP3A7/CYP3A4 selective or common inhibitors and substrates. In summary, the top performing models can be further applied as a tool to rapidly evaluate the safety and efficacy of new drugs separately for fetuses/neonates and adults. The significant structural features could guide the design of new therapeutic drugs as well as aid in the optimization of existing medicine for fetuses/neonates.

TCRD and Pharos 2021: mining the human proteome for disease biology.

In 2014, the National Institutes of Health (NIH) initiated the Illuminating the Druggable Genome (IDG) program to identify and improve our understanding of poorly characterized proteins that can potentially be modulated using small molecules or biologics. Two resources produced from these efforts are: The Target Central Resource Database (TCRD) (http://juniper.health.unm.edu/tcrd/) and Pharos (https://pharos.nih.gov/), a web interface to browse the TCRD. The ultimate goal of these resources is to highlight and facilitate research into currently understudied proteins, by aggregating a multitude of data sources, and ranking targets based on the amount of data available, and presenting data in machine learning ready format. Since the 2017 release, both TCRD and Pharos have produced two major releases, which have incorporated or expanded an additional 25 data sources. Recently incorporated data types include human and viral-human protein-protein interactions, protein-disease and protein-phenotype associations, and drug-induced gene signatures, among others. These aggregated data have enabled us to generate new visualizations and content sections in Pharos, in order to empower users to find new areas of study in the druggable genome.

Therapeutic candidates for the Zika virus identified by a high-throughput screen for Zika protease inhibitors.

When Zika virus emerged as a public health emergency there were no drugs or vaccines approved for its prevention or treatment. We used a high-throughput screen for Zika virus protease inhibitors to identify several inhibitors of Zika virus infection. We expressed the NS2B-NS3 Zika virus protease and conducted a biochemical screen for small-molecule inhibitors. A quantitative structure-activity relationship model was employed to virtually screen ∼138,000 compounds, which increased the identification of active compounds, while decreasing screening time and resources. Candidate inhibitors were validated in several viral infection assays. Small molecules with favorable clinical profiles, especially the five-lipoxygenase-activating protein inhibitor, MK-591, inhibited the Zika virus protease and infection in neural stem cells. Members of the tetracycline family of antibiotics were more potent inhibitors of Zika virus infection than the protease, suggesting they may have multiple mechanisms of action. The most potent tetracycline, methacycline, reduced the amount of Zika virus present in the brain and the severity of Zika virus-induced motor deficits in an immunocompetent mouse model. As Food and Drug Administration-approved drugs, the tetracyclines could be quickly translated to the clinic. The compounds identified through our screening paradigm have the potential to be used as prophylactics for patients traveling to endemic regions or for the treatment of the neurological complications of Zika virus infection.

The AKT modulator A-443654 reduces α-synuclein expression and normalizes ER stress and autophagy.

Accumulation of α-synuclein is a main underlying pathological feature of Parkinson's disease and α-synucleinopathies, for which lowering expression of the α-synuclein gene (SNCA) is a potential therapeutic avenue. Using a cell-based luciferase reporter of SNCA expression we performed a quantitative high-throughput screen of 155,885 compounds and identified A-443654, an inhibitor of the multiple functional kinase AKT, as a potent inhibitor of SNCA. HEK-293 cells with CAG repeat expanded ATXN2 (ATXN2-Q58 cells) have increased levels of α-synuclein. We found that A-443654 normalized levels of both SNCA mRNA and α-synuclein monomers and oligomers in ATXN2-Q58 cells. A-443654 also normalized levels of α-synuclein in fibroblasts and iPSC-derived dopaminergic neurons from a patient carrying a triplication of the SNCA gene. Analysis of autophagy and endoplasmic reticulum stress markers showed that A-443654 successfully prevented α-synuclein toxicity and restored cell function in ATXN2-Q58 cells, normalizing the levels of mTOR, LC3-II, p62, STAU1, BiP, and CHOP. A-443654 also decreased the expression of DCLK1, an inhibitor of α-synuclein lysosomal degradation. Our study identifies A-443654 and AKT inhibition as a potential strategy for reducing SNCA expression and treating Parkinson's disease pathology.

Controlled Proteolysis of an Essential Virulence Determinant Dictates Infectivity of Lyme Disease Pathogens.

The Borrelia burgdorferi BB0323 protein undergoes a complex yet poorly defined proteolytic maturation event that generates N-terminal and C-terminal proteins with essential functions in cell growth and infection. Here, we report that a borrelial protease, B. burgdorferi high temperature requirement A protease (BbHtrA), cleaves BB0323 between asparagine (N) and leucine (L) at positions 236 and 237, while the replacement of these residues with alanine in the mutant protein prevents its cleavage, despite preserving its normal secondary structure. The N-terminal BB0323 protein binds BbHtrA, but its cleavage site mutant displays deficiency in such interaction. An isogenic borrelial mutant with NL-to-AA substitution in BB0323 (referred to as Bbbb0323NL) maintains normal growth yet is impaired for infection of mice or transmission from infected ticks. Notably, the BB0323 protein is still processed in Bbbb0323NL, albeit with lower levels of mature N-terminal BB0323 protein and multiple aberrantly processed polypeptides, which could result from nonspecific cleavages at other asparagine and leucine residues in the protein. The lack of infectivity of Bbbb0323NL is likely due to the impaired abundance or stoichiometry of a protein complex involving BB0238, another spirochete protein. Together, these studies highlight that a precise proteolytic event and a particular protein-protein interaction, involving multiple borrelial virulence determinants, are mutually inclusive and interconnected, playing essential roles in the infectivity of Lyme disease pathogens.

Prediction of drug-induced liver injury and cardiotoxicity using chemical structure and in vitro assay data.

Drug-induced liver injury (DILI) and cardiotoxicity (DICT) are major adverse effects triggered by many clinically important drugs. To provide an alternative to in vivo toxicity testing, the U.S. Tox21 consortium has screened a collection of ∼10K compounds, including drugs in clinical use, against >70 cell-based assays in a quantitative high-throughput screening (qHTS) format. In this study, we compiled reference compound lists for DILI and DICT and compared the potential of Tox21 assay data with chemical structure information in building prediction models for human in vivo hepatotoxicity and cardiotoxicity. Models were built with four different machine learning algorithms (e.g., Random Forest, Naïve Bayes, eXtreme Gradient Boosting, and Support Vector Machine) and model performance was evaluated by calculating the area under the receiver operating characteristic curve (AUC-ROC). Chemical structure-based models showed reasonable predictive power for DILI (best AUC-ROC = 0.75 ± 0.03) and DICT (best AUC-ROC = 0.83 ± 0.03), while Tox21 assay data alone only showed better than random performance. DILI and DICT prediction models built using a combination of assay data and chemical structure information did not have a positive impact on model performance. The suboptimal predictive performance of the assay data is likely due to insufficient coverage of an adequately predictive number of toxicity mechanisms. The Tox21 consortium is currently expanding coverage of biological response space with additional assays that probe toxicologically important targets and under-represented pathways that may improve the prediction of in vivo toxicity such as DILI and DICT.

Synergistic and Antagonistic Drug Combinations against SARS-CoV-2.

Antiviral drug development for coronavirus disease 2019 (COVID-19) is occurring at an unprecedented pace, yet there are still limited therapeutic options for treating this disease. We hypothesized that combining drugs with independent mechanisms of action could result in synergy against SARS-CoV-2, thus generating better antiviral efficacy. Using in silico approaches, we prioritized 73 combinations of 32 drugs with potential activity against SARS-CoV-2 and then tested them in vitro. Sixteen synergistic and eight antagonistic combinations were identified; among 16 synergistic cases, combinations of the US Food and Drug Administration (FDA)-approved drug nitazoxanide with remdesivir, amodiaquine, or umifenovir were most notable, all exhibiting significant synergy against SARS-CoV-2 in a cell model. However, the combination of remdesivir and lysosomotropic drugs, such as hydroxychloroquine, demonstrated strong antagonism. Overall, these results highlight the utility of drug repurposing and preclinical testing of drug combinations for discovering potential therapies to treat COVID-19.

Identification of environmental chemicals that activate p53 signaling after in vitro metabolic activation.

Currently, approximately 80,000 chemicals are used in commerce. Most have little-to-no toxicity information. The U.S. Toxicology in the 21st Century (Tox21) program has conducted a battery of in vitro assays using a quantitative high-throughput screening (qHTS) platform to gain toxicity information on environmental chemicals. Due to technical challenges, standard methods for providing xenobiotic metabolism could not be applied to qHTS assays. To address this limitation, we screened the Tox21 10,000-compound (10K) library, with concentrations ranging from 2.8 nM to 92 µM, using a p53 beta-lactamase reporter gene assay (p53-bla) alone or with rat liver microsomes (RLM) or human liver microsomes (HLM) supplemented with NADPH, to identify compounds that induce p53 signaling after biotransformation. Two hundred and seventy-eight compounds were identified as active under any of these three conditions. Of these 278 compounds, 73 gave more potent responses in the p53-bla assay with RLM, and 2 were more potent in the p53-bla assay with HLM compared with the responses they generated in the p53-bla assay without microsomes. To confirm the role of metabolism in the differential responses, we re-tested these 75 compounds in the absence of NADPH or with heat-attenuated microsomes. Forty-four compounds treated with RLM, but none with HLM, became less potent under these conditions, confirming the role of RLM in metabolic activation. Further evidence of biotransformation was obtained by measuring the half-life of the parent compounds in the presence of microsomes. Together, the data support the use of RLM in qHTS for identifying chemicals requiring biotransformation to induce biological responses.

Compound screening in cell-based models of tau inclusion formation: Comparison of primary neuron and HEK293 cell assays.

The hallmark pathological features of Alzheimer's disease (AD) brains are senile plaques, comprising ß-amyloid (Aß) peptides, and neuronal inclusions formed from tau protein. These plaques form 10-20 years before AD symptom onset, whereas robust tau pathology is more closely associated with symptoms and correlates with cognitive status. This temporal sequence of AD pathology development, coupled with repeated clinical failures of Aß-directed drugs, suggests that molecules that reduce tau inclusions have therapeutic potential. Few tau-directed drugs are presently in clinical testing, in part because of the difficulty in identifying molecules that reduce tau inclusions. We describe here two cell-based assays of tau inclusion formation that we employed to screen for compounds that inhibit tau pathology: a HEK293 cell-based tau overexpression assay, and a primary rat cortical neuron assay with physiological tau expression. Screening a collection of ∼3500 pharmaceutical compounds with the HEK293 cell tau aggregation assay, we obtained only a low number of hit compounds. Moreover, these compounds generally failed to inhibit tau inclusion formation in the cortical neuron assay. We then screened the Prestwick library of mostly approved drugs in the cortical neuron assay, leading to the identification of a greater number of tau inclusion inhibitors. These included four dopamine D2 receptor antagonists, with D2 receptors having previously been suggested to regulate tau inclusions in a Caenorhabditis elegans model. These results suggest that neurons, the cells most affected by tau pathology in AD, are very suitable for screening for tau inclusion inhibitors.

A Universal and High-Throughput Proteomics Sample Preparation Platform.

Major advances have been made to improve the sensitivity of mass analyzers, spectral quality, and speed of data processing enabling more comprehensive proteome discovery and quantitation. While focus has recently begun shifting toward robust proteomics sample preparation efforts, a high-throughput proteomics sample preparation is still lacking. We report the development of a highly automated universal 384-well plate sample preparation platform with high reproducibility and adaptability for extraction of proteins from cells within a culture plate. Digestion efficiency was excellent in comparison to a commercial digest peptide standard with minimal sample loss while improving sample preparation throughput by 20- to 40-fold (the entire process from plated cells to clean peptides is complete in ∼300 min). Analysis of six human cell types, including two primary cell samples, identified and quantified ∼4,000 proteins for each sample in a single high-performance liquid chromatography (HPLC)-tandem mass spectrometry injection with only 100-10K cells, thus demonstrating universality of the platform. The selected protein was further quantified using a developed HPLC-multiple reaction monitoring method for HeLa digests with two heavy labeled internal standard peptides spiked in. Excellent linearity was achieved across different cell numbers indicating a potential for target protein quantitation in clinical research.

Development of Robust Quantitative Structure-Activity Relationship Models for CYP2C9, CYP2D6, and CYP3A4 Catalysis and Inhibition.

Cytochrome P450 enzymes are responsible for the metabolism of >75% of marketed drugs, making it essential to identify the contributions of individual cytochromes P450 to the total clearance of a new candidate drug. Overreliance on one cytochrome P450 for clearance levies a high risk of drug-drug interactions; and considering that several human cytochrome P450 enzymes are polymorphic, it can also lead to highly variable pharmacokinetics in the clinic. Thus, it would be advantageous to understand the likelihood of new chemical entities to interact with the major cytochrome P450 enzymes at an early stage in the drug discovery process. Typical screening assays using human liver microsomes do not provide sufficient information to distinguish the specific cytochromes P450 responsible for clearance. In this regard, we experimentally assessed the metabolic stability of ∼5000 compounds for the three most prominent xenobiotic metabolizing human cytochromes P450, i.e., CYP2C9, CYP2D6, and CYP3A4, and used the data sets to develop quantitative structure-activity relationship models for the prediction of high-clearance substrates for these enzymes. Screening library included the NCATS Pharmaceutical Collection, comprising clinically approved low-molecular-weight compounds, and an annotated library consisting of drug-like compounds. To identify inhibitors, the library was screened against a luminescence-based cytochrome P450 inhibition assay; and through crossreferencing hits from the two assays, we were able to distinguish substrates and inhibitors of these enzymes. The best substrate and inhibitor models (balanced accuracies ∼0.7), as well as the data used to develop these models, have been made publicly available (https://opendata.ncats.nih.gov/adme) to advance drug discovery across all research groups. SIGNIFICANCE STATEMENT: In drug discovery and development, drug candidates with indiscriminate cytochrome P450 metabolic profiles are considered advantageous, since they provide less risk of potential issues with cytochrome P450 polymorphisms and drug-drug interactions. This study developed robust substrate and inhibitor quantitative structure-activity relationship models for the three major xenobiotic metabolizing cytochromes P450, i.e., CYP2C9, CYP2D6, and CYP3A4. The use of these models early in drug discovery will enable project teams to strategize or pivot when necessary, thereby accelerating drug discovery research.