An integrated single-cell RNA-seq map of human neuroblastoma tumors and preclinical models uncovers divergent mesenchymal-like gene expression programs.

Neuroblastoma is a common pediatric cancer, where preclinical studies suggest that a mesenchymal-like gene expression program contributes to chemotherapy resistance. However, clinical outcomes remain poor, implying we need a better understanding of the relationship between patient tumor heterogeneity and preclinical models. Here, we generated single-cell RNA-seq maps of neuroblastoma cell lines, patient-derived xenograft models (PDX), and a genetically engineered mouse model (GEMM). We developed an unsupervised machine learning approach ('automatic consensus nonnegative matrix factorization' (acNMF)) to compare the gene expression programs found in preclinical models to a large cohort of patient tumors. We confirmed a weakly expressed, mesenchymal-like program in otherwise adrenergic cancer cells in some pre-treated high-risk patient tumors, but this appears distinct from the presumptive drug-resistance mesenchymal programs evident in cell lines. Surprisingly however, this weak-mesenchymal-like program was maintained in PDX and could be chemotherapy-induced in our GEMM after only 24 hours, suggesting an uncharacterized therapy-escape mechanism. Collectively, our findings improve the understanding of how neuroblastoma patient tumor heterogeneity is reflected in preclinical models, provides a comprehensive integrated resource, and a generalizable set of computational methodologies for the joint analysis of clinical and pre-clinical single-cell RNA-seq datasets.

Bone morphogenetic protein (BMP) signaling determines neuroblastoma cell fate and sensitivity to retinoic acid.

Retinoic acid (RA) is a standard-of-care neuroblastoma drug thought to be effective by inducing differentiation. Curiously, RA has little effect on primary human tumors during upfront treatment but can eliminate neuroblastoma cells from the bone marrow during post-chemo consolidation therapy-a discrepancy that has never been explained. To investigate this, we treated a large cohort of neuroblastoma cell lines with RA and observed that the most RA-sensitive cells predominantly undergo apoptosis or senescence, rather than differentiation. We conducted genome-wide CRISPR knockout screens under RA treatment, which identified BMP signaling as controlling the apoptosis/senescence vs differentiation cell fate decision and determining RA's overall potency. We then discovered that BMP signaling activity is markedly higher in neuroblastoma patient samples at bone marrow metastatic sites, providing a plausible explanation for RA's ability to clear neuroblastoma cells specifically from the bone marrow, seemingly mimicking interactions between BMP and RA during normal development.

CYP3A5 unexpectedly regulates glucose metabolism through the AKT-TXNIP-GLUT1 axis in pancreatic cancer.

CYP3A5 is a cytochrome P450 (CYP) enzyme that metabolizes drugs and contributes to drug resistance in cancer. However, it remains unclear whether CYP3A5 directly influences cancer progression. In this report, we demonstrate that CYP3A5 regulates glucose metabolism in pancreatic ductal adenocarcinoma. Multi-omics analysis showed that CYP3A5 knockdown results in a decrease in various glucose-related metabolites through its effect on glucose transport. A mechanistic study revealed that CYP3A5 enriches the glucose transporter GLUT1 at the plasma membrane by restricting the translation of TXNIP, a negative regulator of GLUT1. Notably, CYP3A5-generated reactive oxygen species were proved to be responsible for attenuating the AKT-4EBP1-TXNIP signaling pathway. CYP3A5 contributes to cell migration by maintaining high glucose uptake in pancreatic cancer. Taken together, our results, for the first time, reveal a role of CYP3A5 in glucose metabolism in pancreatic ductal adenocarcinoma and identify a novel mechanism that is a potential therapeutic target.

A CRISPR-drug perturbational map for identifying compounds to combine with commonly used chemotherapeutics.

Combination chemotherapy is crucial for successfully treating cancer. However, the enormous number of possible drug combinations means discovering safe and effective combinations remains a significant challenge. To improve this process, we conduct large-scale targeted CRISPR knockout screens in drug-treated cells, creating a genetic map of druggable genes that sensitize cells to commonly used chemotherapeutics. We prioritize neuroblastoma, the most common extracranial pediatric solid tumor, where ~50% of high-risk patients do not survive. Our screen examines all druggable gene knockouts in 18 cell lines (10 neuroblastoma, 8 others) treated with 8 widely used drugs, resulting in 94,320 unique combination-cell line perturbations, which is comparable to the largest existing drug combination screens. Using dense drug-drug rescreening, we find that the top CRISPR-nominated drug combinations are more synergistic than standard-of-care combinations, suggesting existing combinations could be improved. As proof of principle, we discover that inhibition of PRKDC, a component of the non-homologous end-joining pathway, sensitizes high-risk neuroblastoma cells to the standard-of-care drug doxorubicin in vitro and in vivo using patient-derived xenograft (PDX) models. Our findings provide a valuable resource and demonstrate the feasibility of using targeted CRISPR knockout to discover combinations with common chemotherapeutics, a methodology with application across all cancers.

Proteasome inhibition targets the KMT2A transcriptional complex in acute lymphoblastic leukemia.

Rearrangments in Histone-lysine-N-methyltransferase 2A (KMT2Ar) are associated with pediatric, adult and therapy-induced acute leukemias. Infants with KMT2Ar acute lymphoblastic leukemia (ALL) have a poor prognosis with an event-free-survival of 38%. Herein we evaluate 1116 FDA approved compounds in primary KMT2Ar infant ALL specimens and identify a sensitivity to proteasome inhibition. Upon exposure to this class of agents, cells demonstrate a depletion of histone H2B monoubiquitination (H2Bub1) and histone H3 lysine 79 dimethylation (H3K79me2) at KMT2A target genes in addition to a downregulation of the KMT2A gene expression signature, providing evidence that it targets the KMT2A transcriptional complex and alters the epigenome. A cohort of relapsed/refractory KMT2Ar patients treated with this approach on a compassionate basis had an overall response rate of 90%. In conclusion, we report on a high throughput drug screen in primary pediatric leukemia specimens whose results translate into clinically meaningful responses. This innovative treatment approach is now being evaluated in a multi-institutional upfront trial for infants with newly diagnosed ALL.

Cell type identification in spatial transcriptomics data can be improved by leveraging cell-type-informative paired tissue images using a Bayesian probabilistic model.

Spatial transcriptomics technologies have recently emerged as a powerful tool for measuring spatially resolved gene expression directly in tissues sections, revealing cell types and their dysfunction in unprecedented detail. However, spatial transcriptomics technologies are limited in their ability to separate transcriptionally similar cell types and can suffer further difficulties identifying cell types in slide regions where transcript capture is low. Here, we describe a conceptually novel methodology that can computationally integrate spatial transcriptomics data with cell-type-informative paired tissue images, obtained from, for example, the reverse side of the same tissue section, to improve inferences of tissue cell type composition in spatial transcriptomics data. The underlying statistical approach is generalizable to any spatial transcriptomics protocol where informative paired tissue images can be obtained. We demonstrate a use case leveraging cell-type-specific immunofluorescence markers obtained on mouse brain tissue sections and a use case for leveraging the output of AI annotated H&E tissue images, which we used to markedly improve the identification of clinically relevant immune cell infiltration in breast cancer tissue. Thus, combining spatial transcriptomics data with paired tissue images has the potential to improve the identification of cell types and hence to improve the applications of spatial transcriptomics that rely on accurate cell type identification.

Molecular basis of crosstalk in nuclear receptors: heterodimerization between PXR and CAR and the implication in gene regulation.

The 48 human nuclear receptors (NRs) form a superfamily of transcription factors that regulate major physiological and pathological processes. Emerging evidence suggests that NR crosstalk can fundamentally change our understanding of NR biology, but detailed molecular mechanisms of crosstalk are lacking. Here, we report the molecular basis of crosstalk between the pregnane X receptor (PXR) and constitutive androstane receptor (CAR), where they form a novel heterodimer, resulting in their mutual inhibition. PXR and CAR regulate drug metabolism and energy metabolism. Although they have been broadly perceived as functionally redundant, a growing number of reports suggests a mutual inhibitory relation, but their precise mode of coordinated action remains unknown. Using methods including RNA sequencing, small-angle X-ray scattering and crosslinking mass spectrometry we demonstrate that the mutual inhibition altered gene expression globally and is attributed to the novel PXR-CAR heterodimerization via the same interface used by each receptor to heterodimerize with its functional partner, retinoid X receptor (RXR). These findings establish an unexpected functional relation between PXR, CAR and RXR, change the perceived functional relation between PXR and CAR, open new perspectives on elucidating their role and designing approaches to regulate them, and highlight the importance to comprehensively investigate nuclear receptor crosstalk.

The chemotherapeutic CX-5461 primarily targets TOP2B and exhibits selective activity in high-risk neuroblastoma.

Survival in high-risk pediatric neuroblastoma has remained around 50% for the last 20 years, with immunotherapies and targeted therapies having had minimal impact. Here, we identify the small molecule CX-5461 as selectively cytotoxic to high-risk neuroblastoma and synergistic with low picomolar concentrations of topoisomerase I inhibitors in improving survival in vivo in orthotopic patient-derived xenograft neuroblastoma mouse models. CX-5461 recently progressed through phase I clinical trial as a first-in-human inhibitor of RNA-POL I. However, we also use a comprehensive panel of in vitro and in vivo assays to demonstrate that CX-5461 has been mischaracterized and that its primary target at pharmacologically relevant concentrations, is in fact topoisomerase II beta (TOP2B), not RNA-POL I. This is important because existing clinically approved chemotherapeutics have well-documented off-target interactions with TOP2B, which have previously been shown to cause both therapy-induced leukemia and cardiotoxicity-often-fatal adverse events, which can emerge several years after treatment. Thus, while we show that combination therapies involving CX-5461 have promising anti-tumor activity in vivo in neuroblastoma, our identification of TOP2B as the primary target of CX-5461 indicates unexpected safety concerns that should be examined in ongoing phase II clinical trials in adult patients before pursuing clinical studies in children.

Unraveling the Structural Basis of Selective Inhibition of Human Cytochrome P450 3A5.

The human cytochrome P450 (CYP) CYP3A4 and CYP3A5 enzymes metabolize more than one-half of marketed drugs. They share high structural and substrate similarity and are often studied together as CYP3A4/5. However, CYP3A5 preferentially metabolizes several clinically prescribed drugs, such as tacrolimus. Genetic polymorphism in CYP3A5 makes race-based dosing adjustment of tacrolimus necessary to minimize acute rejection after organ transplantation. Moreover, the differential tissue distribution and expression levels of CYP3A4 and CYP3A5 can aggravate toxicity during treatment. Therefore, selective inhibitors of CYP3A5 are needed to distinguish the role of CYP3A5 from that of CYP3A4 and serve as starting points for potential therapeutic development. To this end, we report the crystal structure of CYP3A5 in complex with a previously reported selective inhibitor, clobetasol propionate (CBZ). This is the first CYP3A5 structure with a type I inhibitor, which along with the previously reported substrate-free and type II inhibitor-bound structures, constitute the main CYP3A5 structural modalities. Supported by structure-guided mutagenesis analyses, the CYP3A5-CBZ structure showed that a unique conformation of the F-F' loop in CYP3A5 enables selective binding of CBZ to CYP3A5. Several polar interactions, including hydrogen bonds, stabilize the position of CBZ to interact with this unique F-F' loop conformation. In addition, functional and biophysical assays using CBZ analogs highlight the importance of heme-adjacent moieties for selective CYP3A5 inhibition. Our findings can be used to guide further development of more potent and selective CYP3A5 inhibitors.

Building a Chemical Toolbox for Human Pregnane X Receptor Research: Discovery of Agonists, Inverse Agonists, and Antagonists Among Analogs Based on the Unique Chemical Scaffold of SPA70.

Pregnane X receptor (PXR) plays roles in detoxification and other physiological processes. PXR activation may enhance drug metabolism (leading to adverse drug reactions) or inhibit inflammation. Therefore, PXR agonists, antagonists, and inverse agonists may serve as research tools and drug candidates. However, a specific PXR modulator with an associated structure-activity relationship is lacking. Based on the scaffold of specific human PXR (hPXR) antagonist SPA70 (10), we developed 81 SPA70 analogs and evaluated their receptor-binding and cellular activities. Interestingly, analogs with subtle structural differences displayed divergent cellular activities, including agonistic, dual inverse agonistic and antagonistic, antagonistic, and partial agonistic/partial antagonistic activities (as in compounds 111, 10, 97, and 42, respectively). We generated a pharmacophore model that represents 81 SPA70 analogs, and docking models that correlate strong interactions between the compounds and residues in the AF-2 helix with agonistic activity. These compounds are novel chemical tools for studying hPXR.

Mutation of a single amino acid of pregnane X receptor switches an antagonist to agonist by altering AF-2 helix positioning.

Pregnane X receptor (PXR) is activated by chemicals to transcriptionally regulate drug disposition and possibly decrease drug efficacy and increase resistance, suggesting therapeutic value for PXR antagonists. We previously reported the antagonist SPA70 and its analog SJB7, which unexpectedly is an agonist. Here, we describe another unexpected observation: mutating a single residue (W299A) within the PXR ligand-binding domain converts SPA70 to an agonist. After characterizing wild-type and W299A PXR activity profiles, we used molecular dynamics simulations to reveal that in wild-type PXR, agonists stabilize the activation function 2 (AF-2) helix in an "inward" position, but SPA70 displaces the AF-2. In W299A, however, SPA70 stabilizes the AF-2 "inward", like agonists. We validated our model by predicting the antagonist SJC2 to be a W299A agonist, which was confirmed experimentally. Our work correlates previously unobserved ligand-induced conformational changes to PXR cellular activity and, for the first time, reveals how PXR antagonists work.

KANSL2 and MBNL3 are regulators of pancreatic ductal adenocarcinoma invasion.

Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal forms of cancer. One major reason for this is that PDAC quickly metastasizes to other organs, thereby making its treatment difficult. The molecular machinery driving PDAC metastasis is still poorly understood. In this study, we applied an unbiased approach using CRISPR screening to identify genes that strongly regulate invasion (based on an in vitro assessment of their metastatic potential) in PANC-1, a PDAC cell line. Through CRISPR screening, we identified MBNL3 and KANSL2 as strong regulators of invasion in PANC-1 cells. We further validated MBNL3 and KANSL2 as regulators of PANC-1 cell invasion by using the doxycycline-inducible shRNA system. We also showed that MBNL3 and KANSL2 do not affect cell proliferation. Through our efforts, we have established a process to identify genes that regulate cell invasion and can be further investigated as potential targets for therapeutic intervention.

Clobetasol Propionate Is a Heme-Mediated Selective Inhibitor of Human Cytochrome P450 3A5.

The human cytochrome P450 (CYP) enzymes CYP3A4 and CYP3A5 metabolize most drugs and have high similarities in their structure and substrate preference. Whereas CYP3A4 is predominantly expressed in the liver, CYP3A5 is upregulated in cancer, contributing to drug resistance. Selective inhibitors of CYP3A5 are, therefore, critical to validating it as a therapeutic target. Here we report clobetasol propionate (clobetasol) as a potent and selective CYP3A5 inhibitor identified by high-throughput screening using enzymatic and cell-based assays. Molecular dynamics simulations suggest a close proximity of clobetasol to the heme in CYP3A5 but not in CYP3A4. UV-visible spectroscopy and electron paramagnetic resonance analyses confirmed the formation of an inhibitory type I heme-clobetasol complex in CYP3A5 but not in CYP3A4, thus explaining the CYP3A5 selectivity of clobetasol. Our results provide a structural basis for selective CYP3A5 inhibition, along with mechanistic insights, and highlight clobetasol as an important chemical tool for target validation.

Strategies for developing pregnane X receptor antagonists: Implications from metabolism to cancer.

Pregnane X receptor (PXR) is a ligand-activated nuclear receptor (NR) that was originally identified as a master regulator of xenobiotic detoxification. It regulates the expression of drug-metabolizing enzymes and transporters to control the degradation and excretion of endobiotics and xenobiotics, including therapeutic agents. The metabolism and disposition of drugs might compromise their efficacy and possibly cause drug toxicity and/or drug resistance. Because many drugs can promiscuously bind and activate PXR, PXR antagonists might have therapeutic value in preventing and overcoming drug-induced PXR-mediated drug toxicity and drug resistance. Furthermore, PXR is now known to have broader cellular functions, including the regulation of cell proliferation, and glucose and lipid metabolism. Thus, PXR might be involved in human diseases such as cancer and metabolic diseases. The importance of PXR antagonists is discussed in the context of the role of PXR in xenobiotic sensing and other disease-related pathways. This review focuses on the development of PXR antagonists, which has been hampered by the promiscuity of PXR ligand binding. However, substantial progress has been made in recent years, suggesting that it is feasible to develop selective PXR antagonists. We discuss the current status, challenges, and strategies in developing selective PXR antagonists. The strategies are based on the molecular mechanisms of antagonism in related NRs that can be applied to the design of PXR antagonists, primarily driven by structural information.

Discovery and Characterization of the Antimetabolite Action of Thioacetamide-Linked 1,2,3-Triazoles as Disruptors of Cysteine Biosynthesis in Gram-Negative Bacteria.

Increasing rates of drug-resistant Gram-negative (GN) infections, combined with a lack of new GN-effective antibiotic classes, are driving the need for the discovery of new agents. Bacterial metabolism represents an underutilized mechanism of action in current antimicrobial therapies. Therefore, we sought to identify novel antimetabolites that disrupt key metabolic pathways and explore the specific impacts of these agents on bacterial metabolism. This study describes the successful application of this approach to discover a new series of chemical probes, N-(phenyl)thioacetamide-linked 1,2,3-triazoles (TAT), that target cysteine synthase A (CysK), an enzyme unique to bacteria that is positioned at a key juncture between several fundamental pathways. The TAT class was identified using a high-throughput screen against Escherichia coli designed to identify modulators of pathways related to folate biosynthesis. TAT analog synthesis demonstrated a clear structure-activity relationship, and activity was confirmed against GN antifolate-resistant clinical isolates. Spontaneous TAT resistance mutations were tracked to CysK, and mode of action studies led to the identification of a false product formation mechanism between the CysK substrate O-acetyl-l-serine and the TATs. Global transcriptional responses to TAT treatment revealed that these antimetabolites impose substantial disruption of key metabolic networks beyond cysteine biosynthesis. This study highlights the potential of antimetabolite drug discovery as a promising approach to the discovery of novel GN antibiotics and the pharmacological promise of TAT CysK probes.

Genome-wide CRISPR screen reveals PSMA6 to be an essential gene in pancreatic cancer cells.

BACKGROUND: Despite its relatively low incidence, pancreatic ductal adenocarcinoma (PDAC) is a leading cause of cancer deaths because of the aggressive growth/metastasis of the tumor, the lack of early symptoms, and the poor treatment options. Basic research to identify potential therapeutic targets for PDAC is greatly needed. METHODS: We used a negative-selection genome-wide CRISPR screen to identify essential genes in the PANC-1 human pancreatic carcinoma cell line. We validated the top hits with follow-up siRNA screens, using the HPNE, HPAF-II, AsPC-1, and Mia PaCa-2 cell lines. RESULTS: The PSMA6 gene was an identified candidate hit after the CRISPR screen, siRNA validation screen, and siRNA deconvolution screen. Spheroid formation assays and flow cytometry analysis showed that PSMA6 is critical for survival in many pancreatic ductal carcinoma cell models. Lastly, as PSMA6 protein is a proteosomal subunit of the 20S core complex, we showed that bortezomib, a proteasome inhibitor, was especially toxic in PANC-1 cells. CONCLUSIONS: Further study of PSMA6 and the proteasome subunit that it encodes, along with other hits identified in our CRISPR screens, may provide valuable insights into potential therapeutic targets for PDAC.

Structural Perspectives of the CYP3A Family and Their Small Molecule Modulators in Drug Metabolism.

Cytochrome P450 enzymes function to catalyze a wide range of reactions, many of which are critically important for drug response. Members of the human cytochrome P450 3A (CYP3A) family are particularly important in drug clearance, and they collectively metabolize more than half of all currently prescribed medications. The ability of these enzymes to bind a large and structurally diverse set of compounds increases the chances of their modulating or facilitating drug metabolism in unfavorable ways. Emerging evidence suggests that individual enzymes in the CYP3A family play discrete and important roles in catalysis and disease progression. Here we review the similarities and differences among CYP3A enzymes with regard to substrate recognition, metabolism, modulation by small molecules, and biological consequence, highlighting some of those with clinical significance. We also present structural perspectives to further characterize the basis of these comparisons.

RNA interference screen identifies NAA10 as a regulator of PXR transcription.

The pregnane X receptor (PXR) is a principal xenobiotic receptor crucial in the detection, detoxification, and clearance of toxic substances from the body. PXR plays a vital role in the metabolism and disposition of drugs, and elevated PXR levels contribute to cancer drug resistance. Therefore, to modulate PXR activity and mitigate drug resistance, it is imperative to fully understand its regulation. To this end, we screened a transcription factor siRNA library in pancreatic cancer cells that express high levels of PXR. Through a comprehensive deconvolution process, we identified N-alpha-acetyltransferase 10 (NAA10) as a factor in the transcriptional machinery regulating PXR transcription. Because no one single factor has 100% operational control of PXR transcriptional regulation, our results together with other previous findings suggest that the transcriptional regulation of PXR is complex and that multiple factors contribute to the process including NAA10.

CINPA1 binds directly to constitutive androstane receptor and inhibits its activity.

The constitutive androstane receptor (CAR) and pregnane X receptor (PXR) are xenobiotic sensors that regulate the expression of drug-metabolizing enzymes and efflux transporters. CAR activation promotes drug elimination, thereby reducing therapeutic effectiveness, or causes adverse drug effects via toxic metabolites. CAR inhibitors could be used to attenuate these adverse drug effects. CAR and PXR share ligands and target genes, confounding the understanding of the regulation of receptor-specific activity. We previously identified a small-molecule inhibitor, CINPA1, that inhibits CAR (without activating PXR at lower concentrations) by altering CAR-coregulator interactions and reducing CAR recruitment to DNA response elements of regulated genes. However, solid evidence was not presented for the direct binding of CINPA1 to CAR. In this study, we demonstrate direct interaction of CINPA1 with the CAR ligand-binding domain (CAR-LBD) and identify key residues involved in such interactions through a combination of biophysical and computational methods. We found that CINPA1 resides in the ligand-binding pocket to stabilize the CAR-LBD in a more rigid, less fluid state. Molecular dynamics simulations, together with our previously reported docking model, enabled us to predict which CAR residues were critical for interactions with CINPA1. The importance of these residues for CINPA1 binding were then validated by directed mutations and testing the mutant CAR proteins in transcription reporter and coregulatory interaction assays. We demonstrated strong hydrogen bonding of CINPA1 with N165 and H203 and identified other residues involved in hydrophobic contacts with CINPA1. Overall, our data confirm that CINPA1 directly binds to CAR.