Development and validation of CYP26A1 inhibition assay for high-throughput screening.

All-trans retinoic acid (atRA) is an endogenous ligand of the retinoic acid receptors, which heterodimerize with retinoid X receptors. AtRA is generated in tissues from vitamin A (retinol) metabolism to form a paracrine signal and is locally degraded by cytochrome P450 family 26 (CYP26) enzymes. The CYP26 family consists of three subtypes: A1, B1, and C1, which are differentially expressed during development. This study aims to develop and validate a high throughput screening assay to identify CYP26A1 inhibitors in a cell-free system using a luminescent P450-Glo assay technology. The assay performed well with a signal to background ratio of 25.7, a coefficient of variation of 8.9%, and a Z-factor of 0.7. To validate the assay, we tested a subset of 39 compounds that included known CYP26 inhibitors and retinoids, as well as positive and negative control compounds selected from the literature and/or the ToxCast/Tox21 portfolio. Known CYP26A1 inhibitors were confirmed, and predicted CYP26A1 inhibitors, such as chlorothalonil, prochloraz, and SSR126768, were identified, demonstrating the reliability and robustness of the assay. Given the general importance of atRA as a morphogenetic signal and the localized expression of Cyp26a1 in embryonic tissues, a validated CYP26A1 assay has important implications for evaluating the potential developmental toxicity of chemicals.

Use of Tox21 screening data to profile PFAS bioactivities on nuclear receptors, cellular stress pathways, and cytochrome p450 enzymes.

Per- and poly-fluoroalkyl substances (PFAS) are synthetic chemicals widely used in commercial products. PFAS are a global concern due to their persistence in the environment and extensive associations with adverse health outcomes. While legacy PFAS have been extensively studied, many non-legacy PFAS lack sufficient toxicity information. In this study, we first analyzed the bioactivity of PFAS using Tox21 screening data surveying more than 75 assay endpoints (e.g., nuclear receptors, stress response, and metabolism) to understand the toxicity of non-legacy PFAS and investigate potential new targets of PFAS. From the Tox21 screening data analysis, we confirmed several known PFAS targets/pathways and identified several potential novel targets/pathways of PFAS. To confirm the effect of PFAS on these novel targets/pathways, we conducted several cell- and enzyme-based assays in the follow-up studies. We found PFAS inhibited cytochromes P450s (CYPs), especially CYP2C9 with IC50 values of < 1 µM. Considering PFAS affected other targets/pathways at > 10 µM, PFAS have a higher affinity to CYP2C9. This PFAS-CYP2C9 interaction was further investigated using molecular docking analysis. The result suggested that PFAS directly bind to the active sites of CYP2C9. These findings have important implications to understand the mechanism of PFAS action and toxicity.

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.

Identification of Potent and Selective Acetylcholinesterase/Butyrylcholinesterase Inhibitors by Virtual Screening.

Acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) play important roles in human neurodegenerative disorders such as Alzheimer's disease. In this study, machine learning methods were applied to develop quantitative structure-activity relationship models for the prediction of novel AChE and BChE inhibitors based on data from quantitative high-throughput screening assays. The models were used to virtually screen an in-house collection of ∼360K compounds. The optimal models achieved good performance with area under the receiver operating characteristic curve values ranging from 0.83 ± 0.03 to 0.87 ± 0.01 for the prediction of AChE/BChE inhibition activity and selectivity. Experimental validation showed that the best-performing models increased the assay hit rate by several folds. We identified 88 novel AChE and 126 novel BChE inhibitors, 25% (AChE) and 53% (BChE) of which showed potent inhibitory effects (IC50 < 5 µM). In addition, structure-activity relationship analysis of the BChE inhibitors revealed scaffolds for chemistry design and optimization. In conclusion, machine learning models were shown to efficiently identify potent and selective inhibitors against AChE and BChE and novel structural series for further design and development of potential therapeutics against neurodegenerative disorders.

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.

Screening Method for the Identification of Compounds That Activate Pregnane X Receptor.

The pregnane X receptor (PXR) is a nuclear receptor found mainly in the liver and intestine, whose main function is to regulate the expression of drug-metabolizing enzymes and transporters. Recently, it has been noted that PXR plays critical roles in energy homeostasis, immune response, and cancer. Therefore, identifying chemicals or compounds that can modulate PXR is of great interest, as these can result in downstream toxicity or, alternatively, may have therapeutic potential. Testing one compound at a time for PXR activity would be inefficient and take thousands of hours for large compound libraries. Here, we describe a high-throughput screening method that encompasses plating and treating HepG2-CYP3A4-hPXR cells in a 1536-well plate, as well as reading and interpreting assay (e.g., luciferase reporter gene activity) endpoints. These cells are stably transfected with a human PXR expression vector and CYP3A4-promoter-driven luciferase reporter vector, allowing the identification of compounds that activate PXR through cytochrome 450 3A4. We also describe how to analyze the data from each assay and explain follow-up steps, namely pharmacological characterization and quantitative polymerase chain reaction (qPCR) assays, which can be performed to confirm results from the original screen. These methods can be used to identify and confirm hPXR activators after completion of a compound screening. Published 2022. This article is a U.S. Government work and is in the public domain in the USA. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Establishment of a high-throughput assay to identify hPXR activators Basic Protocol 2: Quantitative high-throughput screening a compound library to classify hPXR activators Basic Protocol 3: Performing pharmacological characterization and qPCR assays to confirm hPXR activators.

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.

Use of Tox21 Screening Data to Evaluate the COVID-19 Drug Candidates for Their Potential Toxic Effects and Related Pathways.

Currently, various potential therapeutic agents for coronavirus disease-2019 (COVID-19), a global pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), are being investigated worldwide mainly through the drug repurposing approach. Several anti-viral, anti-bacterial, anti-malarial, and anti-inflammatory drugs were employed in randomized trials and observational studies for developing new therapeutics for COVID-19. Although an increasing number of repurposed drugs have shown anti-SARS-CoV-2 activities in vitro, so far only remdesivir has been approved by the US FDA to treat COVID-19, and several other drugs approved for Emergency Use Authorization, including sotrovimab, tocilizumab, baricitinib, paxlovid, molnupiravir, and other potential strategies to develop safe and effective therapeutics for SARS-CoV-2 infection are still underway. Many drugs employed as anti-viral may exert unwanted side effects (i.e., toxicity) via unknown mechanisms. To quickly assess these drugs for their potential toxicological effects and mechanisms, we used the Tox21 in vitro assay datasets generated from screening ∼10,000 compounds consisting of approved drugs and environmental chemicals against multiple cellular targets and pathways. Here we summarize the toxicological profiles of small molecule drugs that are currently under clinical trials for the treatment of COVID-19 based on their in vitro activities against various targets and cellular signaling pathways.

Identification of Selective CYP3A7 and CYP3A4 Substrates and Inhibitors Using a High-Throughput Screening Platform.

Cytochrome P450 (CYP) 3A7 is one of the major xenobiotic metabolizing enzymes in human embryonic, fetal, and newborn liver. CYP3A7 expression has also been observed in a subset of the adult population, including pregnant women, as well as in various cancer patients. The characterization of CYP3A7 is not as extensive as other CYPs, and health authorities have yet to provide guidance towards DDI assessment. To identify potential CYP3A7-specific molecules, we used a P450-Glo CYP3A7 enzyme assay to screen a library of ∼5,000 compounds, including FDA-approved drugs and drug-like molecules, and compared these screening data with that from a P450-Glo CYP3A4 assay. Additionally, a subset of 1,000 randomly selected compounds were tested in a metabolic stability assay. By combining the data from the qHTS P450-Glo and metabolic stability assays, we identified several chemical features important for CYP3A7 selectivity. Halometasone was chosen for further evaluation as a potential CYP3A7-selective inhibitor using molecular docking. From the metabolic stability assay, we identified twenty-two CYP3A7-selective substrates over CYP3A4 in supersome setting. Our data shows that CYP3A7 has ligand promiscuity, much like CYP3A4. Furthermore, we have established a large, high-quality dataset that can be used in predictive modeling for future drug metabolism and interaction studies.

Improving the solubility and antileukemia activity of Wnt/ß-catenin signaling inhibitors by disrupting molecular planarity.

Leukemia cells depend on the Wnt/ß-catenin signaling pathway for their growth. Pyrvinium, a known Wnt signaling inhibitor, has demonstrated promising efficacy in the treatment of the aggressive blast phase chronic myeloid leukemia (BP-CML). We previously developed potent inhibitors 1-2 for the Wnt/ß-catenin signaling pathway. However, the further application of these compounds as anti-leukemia agents is limited by their modest anti-leukemia activity in cells and poor aqueous solubility, due to the high molecular planarity of the chemical scaffold. Here, we reported our efforts in the synthesis and in vitro evaluation of 18 new compounds (4a-r) that have been designed to disrupt the molecular planarity of the chemical scaffold. Several compounds of the series showed significantly improved anti-leukemia activity and aqueous solubility. As a highlight, compounds 4c not only maintained excellent inhibitory potency (IC50 = 1.3 nM) for Wnt signaling but also demonstrated good anti-leukemia potency (IC50 = 0.9 µM) in the CML K562 cells. Moreover, compound 4c had an aqueous solubility of 5.9 µg/mL, which is over 50-fold enhanced compared to its parents 1-2.

Targeting CAR and Nrf2 improves cyclophosphamide bioactivation while reducing doxorubicin-induced cardiotoxicity in triple-negative breast cancer treatment.

Cyclophosphamide (CPA) and doxorubicin (DOX) are key components of chemotherapy for triple-negative breast cancer (TNBC), although suboptimal outcomes are commonly associated with drug resistance and/or intolerable side effects. Through an approach combining high-throughput screening and chemical modification, we developed CN06 as a dual activator of the constitutive androstane receptor (CAR) and nuclear factor erythroid 2-related factor 2 (Nrf2). CN06 enhances CAR-induced bioactivation of CPA (a prodrug) by provoking hepatic expression of CYP2B6, while repressing DOX-induced cytotoxicity in cardiomyocytes in vitro via stimulating Nrf2-antioxidant signaling. Utilizing a multicellular coculture model incorporating human primary hepatocytes, TNBC cells, and cardiomyocytes, we show that CN06 increased CPA/DOX-mediated TNBC cell death via CAR-dependent CYP2B6 induction and subsequent conversion of CPA to its active metabolite 4-hydroxy-CPA, while protecting against DOX-induced cardiotoxicity by selectively activating Nrf2-antioxidant signaling in cardiomyocytes but not in TNBC cells. Furthermore, CN06 preserves the viability and function of human iPSC-derived cardiomyocytes by modulating antioxidant defenses, decreasing apoptosis, and enhancing the kinetics of contraction and relaxation. Collectively, our findings identify CAR and Nrf2 as potentially novel combined therapeutic targets whereby CN06 holds the potential to improve the efficacy/toxicity ratio of CPA/DOX-containing chemotherapy.

Mitochondrial Membrane Potential Assay.

Mitochondrial function, a key indicator of cell health, can be assessed through monitoring changes in mitochondrial membrane potential (MMP). Cationic fluorescent dyes are commonly used tools to assess MMP. We used a water-soluble mitochondrial membrane potential indicator (m-MPI) to detect changes in MMP in various types of cells, such as HepG2, HepaRG, and AC16 cells. A homogenous cell-based MMP assay has been optimized and performed in a 1536-well plate format, which can be used to screen several compound libraries for mitochondrial toxicity by evaluating the effects of chemical compounds on MMP.

Identification of Compounds for Butyrylcholinesterase Inhibition.

Butyrylcholinesterase (BChE) is a nonspecific cholinesterase enzyme that hydrolyzes choline-based esters. BChE plays a critical role in maintaining normal cholinergic function like acetylcholinesterase (AChE) through hydrolyzing acetylcholine (ACh). Selective BChE inhibition has been regarded as a viable therapeutic approach in Alzheimer's disease. As of now, a limited number of selective BChE inhibitors are available. To identify BChE inhibitors rapidly and efficiently, we have screened 8998 compounds from several annotated libraries against an enzyme-based BChE inhibition assay in a quantitative high-throughput screening (qHTS) format. From the primary screening, we identified a group of 125 compounds that were further confirmed to inhibit BChE activity, including previously reported BChE inhibitors (e.g., bambuterol and rivastigmine) and potential novel BChE inhibitors (e.g., pancuronium bromide and NNC 756), representing diverse structural classes. These BChE inhibitors were also tested for their selectivity by comparing their IC50 values in BChE and AChE inhibition assays. The binding modes of these compounds were further studied using molecular docking analyses to identify the differences between the interactions of these BChE inhibitors within the active sites of AChE and BChE. Our qHTS approach allowed us to establish a robust and reliable process to screen large compound collections for potential BChE inhibitors.

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.

Predictive Models to Identify Small Molecule Activators and Inhibitors of Opioid Receptors.

Opioid receptors (OPRs) are the main targets for the treatment of pain and related disorders. The opiate compounds that activate these receptors are effective analgesics but their use leads to adverse effects, and they often are highly addictive drugs of abuse. There is an urgent need for alternative chemicals that are analgesics and to reduce/avoid the unwanted effects in order to relieve the public health crisis of opioid addiction. Here, we aim to develop computational models to predict the OPR activity of small molecule compounds based on chemical structures and apply these models to identify novel OPR active compounds. We used four different machine learning algorithms to build models based on quantitative high throughput screening (qHTS) data sets of three OPRs in both agonist and antagonist modes. The best performing models were applied to virtually screen a large collection of compounds. The model predicted active compounds were experimentally validated using the same qHTS assays that generated the training data. Random forest was the best classifier with the highest performance metrics, and the mu OPR (OPRM)-agonist model achieved the best performance measured by AUC-ROC (0.88) and MCC (0.7) values. The model predicted actives resulted in hit rates ranging from 2.3% (delta OPR-agonist) to 15.8% (OPRM-agonist) after experimental confirmation. Compared to the original assay hit rate, all models enriched the hit rate by ≥2-fold. Our approach produced robust OPR prediction models that can be applied to prioritize compounds from large libraries for further experimental validation. The models identified several novel potent compounds as activators/inhibitors of OPRs that were confirmed experimentally. The potent hits were further investigated using molecular docking to find the interactions of the novel ligands in the active site of the corresponding OPR.

Profiling the Tox21 Chemical Collection for Acetylcholinesterase Inhibition.

BACKGROUND: Inhibition of acetylcholinesterase (AChE), a biomarker of organophosphorous and carbamate exposure in environmental and occupational human health, has been commonly used to identify potential safety liabilities. So far, many environmental chemicals, including drug candidates, food additives, and industrial chemicals, have not been thoroughly evaluated for their inhibitory effects on AChE activity. AChE inhibitors can have therapeutic applications (e.g., tacrine and donepezil) or neurotoxic consequences (e.g., insecticides and nerve agents). OBJECTIVES: The objective of the current study was to identify environmental chemicals that inhibit AChE activity using in vitro and in silico models. METHODS: To identify AChE inhibitors rapidly and efficiently, we have screened the Toxicology in the 21st Century (Tox21) 10K compound library in a quantitative high-throughput screening (qHTS) platform by using the homogenous cell-based AChE inhibition assay and enzyme-based AChE inhibition assays (with or without microsomes). AChE inhibitors identified from the primary screening were further tested in monolayer or spheroid formed by SH-SY5Y and neural stem cell models. The inhibition and binding modes of these identified compounds were studied with time-dependent enzyme-based AChE inhibition assay and molecular docking, respectively. RESULTS: A group of known AChE inhibitors, such as donepezil, ambenonium dichloride, and tacrine hydrochloride, as well as many previously unreported AChE inhibitors, such as chelerythrine chloride and cilostazol, were identified in this study. Many of these compounds, such as pyrazophos, phosalone, and triazophos, needed metabolic activation. This study identified both reversible (e.g., donepezil and tacrine) and irreversible inhibitors (e.g., chlorpyrifos and bromophos-ethyl). Molecular docking analyses were performed to explain the relative inhibitory potency of selected compounds. CONCLUSIONS: Our tiered qHTS approach allowed us to generate a robust and reliable data set to evaluate large sets of environmental compounds for their AChE inhibitory activity. https://doi.org/10.1289/EHP6993.

Mining of high throughput screening database reveals AP-1 and autophagy pathways as potential targets for COVID-19 therapeutics.

The recent global pandemic of the Coronavirus disease 2019 (COVID-19) caused by the new coronavirus SARS-CoV-2 presents an urgent need for the development of new therapeutic candidates. Many efforts have been devoted to screening existing drug libraries with the hope to repurpose approved drugs as potential treatments for COVID-19. However, the antiviral mechanisms of action of the drugs found active in these phenotypic screens remain largely unknown. In an effort to deconvolute the viral targets in pursuit of more effective anti-COVID-19 drug development, we mined our in-house database of approved drug screens against 994 assays and compared their activity profiles with the drug activity profile in a cytopathic effect (CPE) assay of SARS-CoV-2. We found that the autophagy and AP-1 signaling pathway activity profiles are significantly correlated with the anti-SARS-CoV-2 activity profile. In addition, a class of neurology/psychiatry drugs was found to be significantly enriched with anti-SARS-CoV-2 activity. Taken together, these results provide new insights into SARS-CoV-2 infection and potential targets for COVID-19 therapeutics, which can be further validated by in vivo animal studies and human clinical trials.

Pharmacological rescue in patient iPSC and mouse models with a rare DISC1 mutation.

We previously identified a causal link between a rare patient mutation in DISC1 (disrupted-in-schizophrenia 1) and synaptic deficits in cortical neurons differentiated from isogenic patient-derived induced pluripotent stem cells (iPSCs). Here we find that transcripts related to phosphodiesterase 4 (PDE4) signaling are significantly elevated in human cortical neurons differentiated from iPSCs with the DISC1 mutation and that inhibition of PDE4 or activation of the cAMP signaling pathway functionally rescues synaptic deficits. We further generated a knock-in mouse line harboring the same patient mutation in the Disc1 gene. Heterozygous Disc1 mutant mice exhibit elevated levels of PDE4s and synaptic abnormalities in the brain, and social and cognitive behavioral deficits. Pharmacological inhibition of the PDE4 signaling pathway rescues these synaptic, social and cognitive behavioral abnormalities. Our study shows that patient-derived isogenic iPSC and humanized mouse disease models are integral and complementary for translational studies with a better understanding of underlying molecular mechanisms.

Characterization of human pregnane X receptor activators identified from a screening of the Tox21 compound library.

The pregnane X receptor (PXR; NR1I2) is an important nuclear receptor whose main function is to regulate enzymes within drug metabolism. The main drug metabolizing enzyme regulated by PXR, cytochrome P450 (CYP) 3A4, accounts for the metabolism of nearly 50% of all marketed drugs. Recently, PXR has also been identified as playing a role in energy homeostasis, immune response, and cancer. Due to its interaction with these important roles, alongside its drug-drug interaction function, it is imperative to identify compounds which can modulate PXR. In this study, we screened the Tox21 10,000 compound collection to identify hPXR agonists using a stable hPXR-Luc HepG2 cell line. A pharmacological study in the presence of a PXR antagonist was performed to confirm the activity of the chosen potential hPXR agonists in the same cells. Finally, metabolically competent cell lines - HepaRG and HepaRG-PXR-Knockout (KO) - were used to further confirm the potential PXR activators. We identified a group of structural clusters and singleton compounds which included potentially novel hPXR agonists. Of the 21 selected compounds, 11 potential PXR activators significantly induced CYP3A4 mRNA expression in HepaRG cells. All of these compounds lost their induction when treating HepaRG-PXR-KO cells, confirming their PXR activation. Etomidoline presented as a potentially selective agonist of PXR. In conclusion, the current study has identified 11 compounds as potentially novel or not well-characterized PXR activators. These compounds should further be studied for their potential effects on drug metabolism and drug-drug interactions due to the immense implications of being a PXR agonist.