Elucidating Compound Mechanism of Action by Network Perturbation Analysis.

Genome-wide identification of the mechanism of action (MoA) of small-molecule compounds characterizing their targets, effectors, and activity modulators represents a highly relevant yet elusive goal, with critical implications for assessment of compound efficacy and toxicity. Current approaches are labor intensive and mostly limited to elucidating high-affinity binding target proteins. We introduce a regulatory network-based approach that elucidates genome-wide MoA proteins based on the assessment of the global dysregulation of their molecular interactions following compound perturbation. Analysis of cellular perturbation profiles identified established MoA proteins for 70% of the tested compounds and elucidated novel proteins that were experimentally validated. Finally, unknown-MoA compound analysis revealed altretamine, an anticancer drug, as an inhibitor of glutathione peroxidase 4 lipid repair activity, which was experimentally confirmed, thus revealing unexpected similarity to the activity of sulfasalazine. This suggests that regulatory network analysis can provide valuable mechanistic insight into the elucidation of small-molecule MoA and compound similarity.

Identification of causal genetic drivers of human disease through systems-level analysis of regulatory networks.

Identification of driver mutations in human diseases is often limited by cohort size and availability of appropriate statistical models. We propose a framework for the systematic discovery of genetic alterations that are causal determinants of disease, by prioritizing genes upstream of functional disease drivers, within regulatory networks inferred de novo from experimental data. We tested this framework by identifying the genetic determinants of the mesenchymal subtype of glioblastoma. Our analysis uncovered KLHL9 deletions as upstream activators of two previously established master regulators of the subtype, C/EBPß and C/EBPδ. Rescue of KLHL9 expression induced proteasomal degradation of C/EBP proteins, abrogated the mesenchymal signature, and reduced tumor viability in vitro and in vivo. Deletions of KLHL9 were confirmed in > 50% of mesenchymal cases in an independent cohort, thus representing the most frequent genetic determinant of the subtype. The method generalized to study other human diseases, including breast cancer and Alzheimer's disease.

TET2 guards against unchecked BATF3-induced CAR T cell expansion.

Further advances in cell engineering are needed to increase the efficacy of chimeric antigen receptor (CAR) and other T cell-based therapies1-5. As T cell differentiation and functional states are associated with distinct epigenetic profiles6,7, we hypothesized that epigenetic programming may provide a means to improve CAR T cell performance. Targeting the gene that encodes the epigenetic regulator ten-eleven translocation 2 (TET2)8 presents an interesting opportunity as its loss may enhance T cell memory9,10, albeit not cause malignancy9,11,12. Here we show that disruption of TET2 enhances T cell-mediated tumour rejection in leukaemia and prostate cancer models. However, loss of TET2 also enables antigen-independent CAR T cell clonal expansions that may eventually result in prominent systemic tissue infiltration. These clonal proliferations require biallelic TET2 disruption and sustained expression of the AP-1 factor BATF3 to drive a MYC-dependent proliferative program. This proliferative state is associated with reduced effector function that differs from both canonical T cell memory13,14 and exhaustion15,16 states, and is prone to the acquisition of secondary somatic mutations, establishing TET2 as a guardian against BATF3-induced CAR T cell proliferation and ensuing genomic instability. Our findings illustrate the potential of epigenetic programming to enhance T cell immunity but highlight the risk of unleashing unchecked proliferative responses.

An extensive microRNA-mediated network of RNA-RNA interactions regulates established oncogenic pathways in glioblastoma.

By analyzing gene expression data in glioblastoma in combination with matched microRNA profiles, we have uncovered a posttranscriptional regulation layer of surprising magnitude, comprising more than 248,000 microRNA (miR)-mediated interactions. These include ∼7,000 genes whose transcripts act as miR "sponges" and 148 genes that act through alternative, nonsponge interactions. Biochemical analyses in cell lines confirmed that this network regulates established drivers of tumor initiation and subtype implementation, including PTEN, PDGFRA, RB1, VEGFA, STAT3, and RUNX1, suggesting that these interactions mediate crosstalk between canonical oncogenic pathways. siRNA silencing of 13 miR-mediated PTEN regulators, whose locus deletions are predictive of PTEN expression variability, was sufficient to downregulate PTEN in a 3'UTR-dependent manner and to increase tumor cell growth rates. Thus, miR-mediated interactions provide a mechanistic, experimentally validated rationale for the loss of PTEN expression in a large number of glioma samples with an intact PTEN locus.

Axonal precursor miRNAs hitchhike on endosomes and locally regulate the development of neural circuits.

Various species of non-coding RNAs (ncRNAs) are enriched in specific subcellular compartments, but the mechanisms orchestrating their localization and their local functions remain largely unknown. We investigated both aspects using the elongating retinal ganglion cell axon and its tip, the growth cone, as models. We reveal that specific endogenous precursor microRNAs (pre-miRNAs) are actively trafficked to distal axons by hitchhiking primarily on late endosomes/lysosomes. Upon exposure to the axon guidance cue semaphorin 3A (Sema3A), pre-miRNAs are processed specifically within axons into newly generated miRNAs, one of which, in turn, silences the basal translation of tubulin beta 3 class III (TUBB3), but not amyloid beta precursor protein (APP). At the organismal level, these mature miRNAs are required for growth cone steering and a fully functional visual system. Overall, our results uncover a novel mode of ncRNA transport from one cytosolic compartment to another within polarized cells. They also reveal that newly generated miRNAs are critical components of a ncRNA-based signaling pathway that transduces environmental signals into the structural remodeling of subcellular compartments.

An evolutionarily conserved Lhx2-Ldb1 interaction regulates the acquisition of hippocampal cell fate and regional identity.

The protein co-factor Ldb1 regulates cell fate specification by interacting with LIM-homeodomain (LIM-HD) proteins in a tetrameric complex consisting of an LDB:LDB dimer that bridges two LIM-HD molecules, a mechanism first demonstrated in the Drosophila wing disc. Here, we demonstrate conservation of this interaction in the regulation of mammalian hippocampal development, which is profoundly defective upon loss of either Lhx2 or Ldb1 Electroporation of a chimeric construct that encodes the Lhx2-HD and Ldb1-DD (dimerization domain) in a single transcript cell-autonomously rescues a comprehensive range of hippocampal deficits in the mouse Ldb1 mutant, including the acquisition of field-specific molecular identity and the regulation of the neuron-glia cell fate switch. This demonstrates that the LHX:LDB complex is an evolutionarily conserved molecular regulatory device that controls complex aspects of regional cell identity in the developing brain.

Knowledge-based classification of fine-grained immune cell types in single-cell RNA-Seq data.

Single-cell RNA sequencing (scRNA-Seq) is an emerging strategy for characterizing immune cell populations. Compared to flow or mass cytometry, scRNA-Seq could potentially identify cell types and activation states that lack precise cell surface markers. However, scRNA-Seq is currently limited due to the need to manually classify each immune cell from its transcriptional profile. While recently developed algorithms accurately annotate coarse cell types (e.g. T cells versus macrophages), making fine distinctions (e.g. CD8+ effector memory T cells) remains a difficult challenge. To address this, we developed a machine learning classifier called ImmClassifier that leverages a hierarchical ontology of cell type. We demonstrate that its predictions are highly concordant with flow-based markers from CITE-seq and outperforms other tools (+15% recall, +14% precision) in distinguishing fine-grained cell types with comparable performance on coarse ones. Thus, ImmClassifier can be used to explore more deeply the heterogeneity of the immune system in scRNA-Seq experiments.

Discovery of a new flavin N5-adduct in a tyrosine to phenylalanine variant of d-Arginine dehydrogenase.

d-Arginine dehydrogenase from Pseudomonas aeruginosa (PaDADH) catalyzes the flavin-dependent oxidation of d-arginine and other d-amino acids. Here, we report the crystal structure at 1.29 Å resolution for PaDADH-Y249F expressed and co-crystallized with d-arginine. The overall structure of PaDADH-Y249F resembled PaDADH-WT, but the electron density for the flavin cofactor was ambiguous, suggesting the presence of modified flavins. Electron density maps and mass spectrometric analysis confirmed the presence of both N5-(4-guanidino-oxobutyl)-FAD and 6-OH-FAD in a single crystal of PaDADH-Y249F and helped with the further refinement of the X-ray crystal structure. The versatility of the reduced flavin is apparent in the PaDADH-Y249F structure and is evidenced by the multiple functions it can perform in the same active site.

A Single-Point Mutation in d-Arginine Dehydrogenase Unlocks a Transient Conformational State Resulting in Altered Cofactor Reactivity.

Proteins are inherently dynamic, and proper enzyme function relies on conformational flexibility. In this study, we demonstrated how an active site residue changes an enzyme's reactivity by modulating fluctuations between conformational states. Replacement of tyrosine 249 (Y249) with phenylalanine in the active site of the flavin-dependent d-arginine dehydrogenase yielded an enzyme with both an active yellow FAD (Y249F-y) and an inactive chemically modified green FAD, identified as 6-OH-FAD (Y249F-g) through various spectroscopic techniques. Structural investigation of Y249F-g and Y249F-y variants by comparison to the wild-type enzyme showed no differences in the overall protein structure and fold. A closer observation of the active site of the Y249F-y enzyme revealed an alternative conformation for some active site residues and the flavin cofactor. Molecular dynamics simulations probed the alternate conformations observed in the Y249F-y enzyme structure and showed that the enzyme variant with FAD samples a metastable conformational state, not available to the wild-type enzyme. Hybrid quantum/molecular mechanical calculations identified differences in flavin electronics between the wild type and the alternate conformation of the Y249F-y enzyme. The computational studies further indicated that the alternate conformation in the Y249F-y enzyme is responsible for the higher spin density at the C6 atom of flavin, which is consistent with the formation of 6-OH-FAD in the variant enzyme. The observations in this study are consistent with an alternate conformational space that results in fine-tuning the microenvironment around a versatile cofactor playing a critical role in enzyme function.

Cupid: simultaneous reconstruction of microRNA-target and ceRNA networks.

We introduce a method for simultaneous prediction of microRNA-target interactions and their mediated competitive endogenous RNA (ceRNA) interactions. Using high-throughput validation assays in breast cancer cell lines, we show that our integrative approach significantly improves on microRNA-target prediction accuracy as assessed by both mRNA and protein level measurements. Our biochemical assays support nearly 500 microRNA-target interactions with evidence for regulation in breast cancer tumors. Moreover, these assays constitute the most extensive validation platform for computationally inferred networks of microRNA-target interactions in breast cancer tumors, providing a useful benchmark to ascertain future improvements.

Urinary levels of triclosan and triclocarban in several Asian countries, Greece and the USA: Association with oxidative stress.

Triclosan (TCS) and Triclocarban (TCC) are widely used as antimicrobial preservatives in personal care products (PCPs). Because of their potential for endocrine disrupting effects, human exposure to these chemicals is a concern. Biomonitoring studies of human exposure to TCS and TCC have shown widespread exposure of populations in western European countries and the USA. However, exposure to TCC and TCS by populations in Asian countries is less well known. In this study, concentrations of TCS and TCC were determined in human urine collected from seven Asian countries (China, India, Korea, Kuwait, Japan, Saudi Arabia, and Vietnam), and Greece and the USA. A total of 430 urine samples were analyzed for TCS and TCC, of which 355 (83%) and 82 (19%), respectively, contained measurable levels of these chemicals. The overall geometric mean [GM] concentrations of TCS and TCC, were 1.36 and 0.03ng/mL, respectively. The highest mean concentration of TCS was found in urine from China (100ng/mL) and the lowest concentration was found in urine from Vietnam (2.34ng/mL). We also analyzed urinary 8-OHdG, a marker of oxidative stress, to elucidate the association with TCS and TCC levels for samples from Saudi Arabia (n=130) and a positive correlation between Ln-transformed TCC levels and 8-OHdG was found, although this was not statistically significant. This is the first study to report urinary levels of TCS and TCC in several Asian countries, especially for Vietnam, Kuwait, and Japan.

Associations of Vitamin D Levels and Vitamin D Receptor Genotypes with Patient-Reported Outcome/Disease Activity in Patients with Rheumatoid Arthritis.

BACKGROUND: The aim of this study is to evaluate the prevalence of vitamin-D insufficiency and vitamin-D receptor (VDR) polymorphisms in rheumatoid arthritis (RA) patients and its association with disease activity and patient reported outcomes (PROs). METHODS: Eighty-two individuals were included in a cross-sectional study (41 RA patients, 41 controls). Prior to assessment, each patient completed a PRO questionnaire. Serum vitamin-D levels and genotyping for VDR were assessed. Vitamin-D deficient patients received vitamin-D supplementation. Re-assessment of disease activity (DAS28) was performed after 9-months. RESULTS: Low vitamin-D levels were more frequent in RA patients (p < 0.01). A negative, but insignificant, association with DAS-28 score was identified; whereas, there was a significant negative association with the PROs (p < 0.01). Vitamin-D supplementation was associated with significant improvement in the patients' scores for pain, fatigue, global assessment, physical disability, and quality of life. In contrast to the control group, the frequency of the recessive TaqI and FokI genotypes was higher in RA patients. CONCLUSIONS: In RA patients, serum vitamin-D level was significantly and inversely associated with both PROs and disease activity. The TaqI and FokI fragment length polymorphisms of VDR significantly contributed to the risk of RA. Having a significant positive impact on patient reported outcomes, vitamin-D supplementation may have a role in RA management.

A new mitochondrial pool of cyclin E, regulated by Drp1, is linked to cell-density-dependent cell proliferation.

The regulation and function of the crucial cell cycle regulator cyclin E (CycE) remains elusive. Unlike other cyclins, CycE can be uniquely controlled by mitochondrial energetics, the exact mechanism being unclear. Using mammalian cells (in vitro) and Drosophila (in vivo) model systems in parallel, we show that CycE can be directly regulated by mitochondria through its recruitment to the organelle. Active mitochondrial bioenergetics maintains a distinct mitochondrial pool of CycE (mtCycE) lacking a key phosphorylation required for its degradation. Loss of the mitochondrial fission protein dynamin-related protein 1 (Drp1, SwissProt O00429 in humans) augments mitochondrial respiration and elevates the mtCycE pool allowing CycE deregulation, cell cycle alterations and enrichment of stem cell markers. Such CycE deregulation after Drp1 loss attenuates cell proliferation in low-cell-density environments. However, in high-cell-density environments, elevated MEK-ERK signaling in the absence of Drp1 releases mtCycE to support escape of contact inhibition and maintain aberrant cell proliferation. Such Drp1-driven regulation of CycE recruitment to mitochondria might be a mechanism to modulate CycE degradation during normal developmental processes as well as in tumorigenic events.

Amine oxidation by d-arginine dehydrogenase in Pseudomonas aeruginosa.

d-Arginine dehydrogenase from Pseudomonas aeruginosa (PaDADH) is a flavin-dependent oxidoreductase, which is part of a novel two-enzyme racemization system that functions to convert d-arginine to l-arginine. PaDADH contains a noncovalently linked FAD that shows the highest activity with d-arginine. The enzyme exhibits broad substrate specificity towards d-amino acids, particularly with cationic and hydrophobic d-amino acids. Biochemical studies have established the structure and the mechanistic properties of the enzyme. The enzyme is a true dehydrogenase because it displays no reactivity towards molecular oxygen. As established through solvent and multiple kinetic isotope studies, PaDADH catalyzes an asynchronous CH and NH bond cleavage via a hydride transfer mechanism. Steady-state kinetic studies with d-arginine and d-histidine are consistent with the enzyme following a ping-pong bi-bi mechanism. As shown by a combination of crystallography, kinetic and computational data, the shape and flexibility of loop L1 in the active site of PaDADH are important for substrate capture and broad substrate specificity.

Association of Paraoxonase 1 Polymorphism and Serum 25-Hydroxyvitamin D with the Risk of Cardiovascular Disease in Patients with Rheumatoid Arthritis.

BACKGROUND: Patients with rheumatoid arthritis (RA) have significantly increased cardiovascular (CV) morbidity and mortality that are not accounted for by traditional risk factors alone. Paraoxonase 1 (PON1) and 25-hydroxyvitamin D have been shown to be involved in the pathogenesis of CV diseases. Objective: This study aimed to investigate PON1 gene polymorphism and serum 25-hydroxyvitamin D concentrations in RA patients, and to determine their association with CV risk in RA. METHODS: Serum samples from 46 RA patients and 45 healthy controls were tested for PON1 R192Q genotypes and serum vitamin D concentrations. The cardiovascular risks were assessed by Q-risk. Lipoprotein cholesterol levels, traditional CV risk factors, medication use, and RA disease activity status were also assessed. RESULTS: PON1 polymorphism and low serum 25-hydroxyvitamin D were significantly associated with increased CV risk (p < 0.05). Compared to patients with either the PON1 QQ genotype or the QR genotype, patients with the RR genotype demonstrated decreased CV risk on multivariate analysis, controlling for traditional CV risk factors, C-reactive protein levels, prednisone use, and cholesterol-lowering medication use (p < 0.05). CONCLUSIONS: There was a relationship of the genetic determinants of paraoxonase 1 (PON1 192) and serum 25-hydroxyvitamin D to CV risk in RA patients. Paired measurement of paraoxonase 1 genotype and serum 25-hydroxyvitamin D can be used as biomarkers of CV risk in RA patients.

Urinary biomarkers of exposure to 57 xenobiotics and its association with oxidative stress in a population in Jeddah, Saudi Arabia.

Oxidative stress arises from excessive free radicals in the body and is a trigger for numerous diseases, such as cancer and atherosclerosis. Elevated exposure to environmental chemicals can contribute to oxidative stress. The association between exposure to xenobiotics and oxidative stress, however, has rarely been studied. In this study, urinary concentrations of 57 xenobiotics (antimicrobials, parabens, bisphenols, benzophenones, and phthalates metabolites) were determined in a population from Jeddah, Saudi Arabia, to delineate association with the oxidative stress biomarker, 8-hydroxy-2'-deoxyguanosine (8OHDG). We collected 130 urine samples and analyzed for 57 xenobiotics using liquid chromatography-tandem mass spectrometry (LC/MS/MS) methods. The association between unadjusted and creatinine- or specific gravity-adjusted concentrations of xenobiotics and 8OHDG was examined by Pearson correlations and multiple regression analysis. High concentrations of mCPP (a metabolite of di-n-octyl phthalate; DnOP) and mCMHP (a metabolite of diethylhexyl phthalate; DEHP) were found in urine. In addition, the concentrations of bisphenol S (BPS) were higher than those of bisphenol A (BPA). The concentrations of metabolites of DEHP, phthalic acid, BPA, BPS, and methyl-protocatechuic acid were significantly associated with 8OHDG. This is the first biomonitoring study to report exposure of the Saudi population to a wide range of environmental chemicals and provides evidence that environmental chemical exposures contribute to oxidative stress.

Direct ChIP-Seq significance analysis improves target prediction.

BACKGROUND: Chromatin immunoprecipitation followed by sequencing of protein-bound DNA fragments (ChIP-Seq) is an effective high-throughput methodology for the identification of context specific DNA fragments that are bound by specific proteins in vivo. Despite significant progress in the bioinformatics analysis of this genome-scale data, a number of challenges remain as technology-dependent biases, including variable target accessibility and mappability, sequence-dependent variability, and non-specific binding affinity must be accounted for. RESULTS AND DISCUSSION: We introduce a nonparametric method for scoring consensus regions of aligned immunoprecipitated DNA fragments when appropriate control experiments are available. Our method uses local models for null binding; these are necessary because binding prediction scores based on global models alone fail to properly account for specialized features of genomic regions and chance pull downs of specific DNA fragments, thus disproportionally rewarding some genomic regions and decreasing prediction accuracy. We make no assumptions about the structure or amplitude of bound peaks, yet we show that our method outperforms leading methods developed using either global or local null hypothesis models for random binding. We test prediction performance by comparing analyses of ChIP-seq, ChIP-chip, motif-based binding-site prediction, and shRNA assays, showing high reproducibility, binding-site enrichment in predicted target regions, and functional regulation of predicted targets. CONCLUSIONS: Given appropriate controls, a direct nonparametric method for identifying transcription-factor targets from ChIP-Seq assays may lead to both higher sensitivity and higher specificity, and should be preferred or used in conjunction with methods that use parametric models for null binding.

Noninvasive In Vivo Imaging of T-Cells during Cancer Immunotherapy Using Rare-Earth Nanoparticles.

Only a minority of patients respond positively to cancer immunotherapy, and addressing this variability is an active area of immunotherapy research. Infiltration of tumors by immune cells is one of the most significant prognostic indicators of response and disease-free survival. However, the ability to noninvasively sample the tumor microenvironment for immune cells remains limited. Imaging in the near-infrared-II region using rare-earth nanocrystals is emerging as a powerful imaging tool for high-resolution deep-tissue imaging. In this paper, we demonstrate that these nanoparticles can be used for noninvasive in vivo imaging of tumor-infiltrating T-cells in a highly aggressive melanoma tumor model. We present nanoparticle synthesis and surface modification strategies for the generation of small, ultrabright, and biocompatible rare-earth nanocrystals necessary for deep tissue imaging of rare cell types. The ability to noninvasively monitor the immune contexture of a tumor during immunotherapy could lead to early identification of nonresponding patients in real time, leading to earlier interventions and better outcomes.

Absence of Tec family kinases interleukin-2 inducible T cell kinase (Itk) and Bruton's tyrosine kinase (Btk) severely impairs Fc epsilonRI-dependent mast cell responses.

Mast cells are critical effector cells in the pathophysiology of allergic asthma and other IgE-mediated diseases. The Tec family of tyrosine kinases Itk and Btk serve as critical signal amplifiers downstream of antigen receptors. Although both kinases are expressed and activated in mast cells following FcεRI stimulation, their individual contributions are not clear. To determine whether these kinases play unique and/or complementary roles in FcεRI signaling and mast cell function, we generated Itk and Btk double knock-out mice. Analyses of these mice show decreased mast cell granularity and impaired passive systemic anaphylaxis responses. This impaired response is accompanied by a significant elevation in serum IgE in Itk/Btk double knock-out mice. In vitro analyses of bone marrow-derived mast cells (BMMCs) indicated that Itk/Btk double knock-out BMMCs are defective in degranulation and cytokine secretion responses downstream to FcεRI activation. These responses were accompanied by a significant reduction in PLCγ2 phosphorylation and severely impaired calcium responses in these cells. This defect also results in altered NFAT1 nuclear localization in double knock-out BMMCs. Network analysis suggests that although they may share substrates, Itk plays both positive and negative roles, while Btk primarily plays a positive role in mast cell FcεRI-induced cytokine secretion.