The BioPlex Network: A Systematic Exploration of the Human Interactome.

Protein interactions form a network whose structure drives cellular function and whose organization informs biological inquiry. Using high-throughput affinity-purification mass spectrometry, we identify interacting partners for 2,594 human proteins in HEK293T cells. The resulting network (BioPlex) contains 23,744 interactions among 7,668 proteins with 86% previously undocumented. BioPlex accurately depicts known complexes, attaining 80%-100% coverage for most CORUM complexes. The network readily subdivides into communities that correspond to complexes or clusters of functionally related proteins. More generally, network architecture reflects cellular localization, biological process, and molecular function, enabling functional characterization of thousands of proteins. Network structure also reveals associations among thousands of protein domains, suggesting a basis for examining structurally related proteins. Finally, BioPlex, in combination with other approaches, can be used to reveal interactions of biological or clinical significance. For example, mutations in the membrane protein VAPB implicated in familial amyotrophic lateral sclerosis perturb a defined community of interactors.

Cdk1 Controls Global Epigenetic Landscape in Embryonic Stem Cells.

The cyclin-dependent kinase 1 (Cdk1) drives cell division. To uncover additional functions of Cdk1, we generated knockin mice expressing an analog-sensitive version of Cdk1 in place of wild-type Cdk1. In our study, we focused on embryonic stem cells (ESCs), because this cell type displays particularly high Cdk1 activity. We found that in ESCs, a large fraction of Cdk1 substrates is localized on chromatin. Cdk1 phosphorylates many proteins involved in epigenetic regulation, including writers and erasers of all major histone marks. Consistent with these findings, inhibition of Cdk1 altered histone-modification status of ESCs. High levels of Cdk1 in ESCs phosphorylate and partially inactivate Dot1l, the H3K79 methyltransferase responsible for placing activating marks on gene bodies. Decrease of Cdk1 activity during ESC differentiation de-represses Dot1l, thereby allowing coordinated expression of differentiation genes. These analyses indicate that Cdk1 functions to maintain the epigenetic identity of ESCs.

Genome-wide transcript and protein analysis highlights the role of protein homeostasis in the aging mouse heart.

Investigation of the molecular mechanisms of aging in the human heart is challenging because of confounding factors, such as diet and medications, as well as limited access to tissues from healthy aging individuals. The laboratory mouse provides an ideal model to study aging in healthy individuals in a controlled environment. However, previous mouse studies have examined only a narrow range of the genetic variation that shapes individual differences during aging. Here, we analyze transcriptome and proteome data from 185 genetically diverse male and female mice at ages 6, 12, and 18 mo to characterize molecular changes that occur in the aging heart. Transcripts and proteins reveal activation of pathways related to exocytosis and cellular transport with age, whereas processes involved in protein folding decrease with age. Additional changes are apparent only in the protein data including reduced fatty acid oxidation and increased autophagy. For proteins that form complexes, we see a decline in correlation between their component subunits with age, suggesting age-related loss of stoichiometry. The most affected complexes are themselves involved in protein homeostasis, which potentially contributes to a cycle of progressive breakdown in protein quality control with age. Our findings highlight the important role of post-transcriptional regulation in aging. In addition, we identify genetic loci that modulate age-related changes in protein homeostasis, suggesting that genetic variation can alter the molecular aging process.

Selective inhibitors of JAK1 targeting an isoform-restricted allosteric cysteine.

The Janus tyrosine kinase (JAK) family of non-receptor tyrosine kinases includes four isoforms (JAK1, JAK2, JAK3, and TYK2) and is responsible for signal transduction downstream of diverse cytokine receptors. JAK inhibitors have emerged as important therapies for immun(onc)ological disorders, but their use is limited by undesirable side effects presumed to arise from poor isoform selectivity, a common challenge for inhibitors targeting the ATP-binding pocket of kinases. Here we describe the chemical proteomic discovery of a druggable allosteric cysteine present in the non-catalytic pseudokinase domain of JAK1 (C817) and TYK2 (C838), but absent from JAK2 or JAK3. Electrophilic compounds selectively engaging this site block JAK1-dependent trans-phosphorylation and cytokine signaling, while appearing to act largely as 'silent' ligands for TYK2. Importantly, the allosteric JAK1 inhibitors do not impair JAK2-dependent cytokine signaling and are inactive in cells expressing a C817A JAK1 mutant. Our findings thus reveal an allosteric approach for inhibiting JAK1 with unprecedented isoform selectivity.

Quantitative analysis of Y-Chromosome gene expression across 36 human tissues.

Little is known about how human Y-Chromosome gene expression directly contributes to differences between XX (female) and XY (male) individuals in nonreproductive tissues. Here, we analyzed quantitative profiles of Y-Chromosome gene expression across 36 human tissues from hundreds of individuals. Although it is often said that Y-Chromosome genes are lowly expressed outside the testis, we report many instances of elevated Y-Chromosome gene expression in a nonreproductive tissue. A notable example is EIF1AY, which encodes eukaryotic translation initiation factor 1A Y-linked, together with its X-linked homolog EIF1AX Evolutionary loss of a Y-linked microRNA target site enabled up-regulation of EIF1AY, but not of EIF1AX, in the heart. Consequently, this essential translation initiation factor is nearly twice as abundant in male as in female heart tissue at the protein level. Divergence between the X and Y Chromosomes in regulatory sequence can therefore lead to tissue-specific Y-Chromosome-driven sex biases in expression of critical, dosage-sensitive regulatory genes.

The metabolic function of cyclin D3-CDK6 kinase in cancer cell survival.

D-type cyclins (D1, D2 and D3) and their associated cyclin-dependent kinases (CDK4 and CDK6) are components of the core cell cycle machinery that drives cell proliferation. Inhibitors of CDK4 and CDK6 are currently being tested in clinical trials for patients with several cancer types, with promising results. Here, using human cancer cells and patient-derived xenografts in mice, we show that the cyclin D3-CDK6 kinase phosphorylates and inhibits the catalytic activity of two key enzymes in the glycolytic pathway, 6-phosphofructokinase and pyruvate kinase M2. This re-directs the glycolytic intermediates into the pentose phosphate (PPP) and serine pathways. Inhibition of cyclin D3-CDK6 in tumour cells reduces flow through the PPP and serine pathways, thereby depleting the antioxidants NADPH and glutathione. This, in turn, increases the levels of reactive oxygen species and causes apoptosis of tumour cells. The pro-survival function of cyclin D-associated kinase operates in tumours expressing high levels of cyclin D3-CDK6 complexes. We propose that measuring the levels of cyclin D3-CDK6 in human cancers might help to identify tumour subsets that undergo cell death and tumour regression upon inhibition of CDK4 and CDK6. Cyclin D3-CDK6, through its ability to link cell cycle and cell metabolism, represents a particularly powerful oncoprotein that affects cancer cells at several levels, and this property can be exploited for anti-cancer therapy.

Fingolimod effects on the brain are mediated through biochemical modulation of bioenergetics, autophagy, and neuroinflammatory networks.

Fingolimod (FTY720) is an oral drug approved by the Food and Drug Administration (FDA) for management of multiple sclerosis (MS) symptoms, which has also shown beneficial effects against Alzheimer's (AD) and Parkinson's (PD) diseases pathologies. Although an extensive effort has been made to identify mechanisms underpinning its therapeutic effects, much remains unknown. Here, we investigated Fingolimod induced proteome changes in the cerebellum (CB) and frontal cortex (FC) regions of the brain which are known to be severely affected in MS, using a tandem mass tag (TMT) isobaric labeling-based quantitative mass-spectrometric approach to investigate the mechanism of action of Fingolimod. This study identified 6749 and 6319 proteins in CB and FC, respectively, and returned 2609 and 3086 differentially expressed proteins in mouse CB and FC, respectively, between Fingolimod treated and control groups. Subsequent bioinformatics analyses indicated a metabolic reprogramming in both brain regions of the Fingolimod treated group, where oxidative phosphorylation was upregulated while glycolysis and pentose phosphate pathway were downregulated. In addition, modulation of neuroinflammation in the Fingolimod treated group was indicated by upregulation of retrograde endocannabinoid signaling and autophagy pathways, and downregulation of neuroinflammation related pathways including neutrophil degranulation and the IL-12 mediated signaling pathway. Our findings suggest that Fingolimod may exert its protective effects on the brain by inducing metabolic reprogramming and neuroinflammation pathway modulation.

Defining the consequences of genetic variation on a proteome-wide scale.

Genetic variation modulates protein expression through both transcriptional and post-transcriptional mechanisms. To characterize the consequences of natural genetic diversity on the proteome, here we combine a multiplexed, mass spectrometry-based method for protein quantification with an emerging outbred mouse model containing extensive genetic variation from eight inbred founder strains. By measuring genome-wide transcript and protein expression in livers from 192 Diversity outbred mice, we identify 2,866 protein quantitative trait loci (pQTL) with twice as many local as distant genetic variants. These data support distinct transcriptional and post-transcriptional models underlying the observed pQTL effects. Using a sensitive approach to mediation analysis, we often identified a second protein or transcript as the causal mediator of distant pQTL. Our analysis reveals an extensive network of direct protein-protein interactions. Finally, we show that local genotype can provide accurate predictions of protein abundance in an independent cohort of collaborative cross mice.

Kinase-independent function of E-type cyclins in liver cancer.

E-type cyclins (cyclins E1 and E2) are components of the core cell cycle machinery and are overexpressed in many human tumor types. E cyclins are thought to drive tumor cell proliferation by activating the cyclin-dependent kinase 2 (CDK2). The cyclin E1 gene represents the site of recurrent integration of the hepatitis B virus in the pathogenesis of hepatocellular carcinoma, and this event is associated with strong up-regulation of cyclin E1 expression. Regardless of the underlying mechanism of tumorigenesis, the majority of liver cancers overexpress E-type cyclins. Here we used conditional cyclin E knockout mice and a liver cancer model to test the requirement for the function of E cyclins in liver tumorigenesis. We show that a ubiquitous, global shutdown of E cyclins did not visibly affect postnatal development or physiology of adult mice. However, an acute ablation of E cyclins halted liver cancer progression. We demonstrated that also human liver cancer cells critically depend on E cyclins for proliferation. In contrast, we found that the function of the cyclin E catalytic partner, CDK2, is dispensable in liver cancer cells. We observed that E cyclins drive proliferation of tumor cells in a CDK2- and kinase-independent mechanism. Our study suggests that compounds which degrade or inhibit cyclin E might represent a highly selective therapeutic strategy for patients with liver cancer, as these compounds would selectively cripple proliferation of tumor cells, while sparing normal tissues.

Retinal proteomics of experimental glaucoma model reveal intraocular pressure-induced mediators of neurodegenerative changes.

Current evidence suggests that exposure to chronically induced intraocular pressure (IOP) leads to neurodegenerative changes in the inner retina. This study aimed to determine retinal proteomic alterations in a rat model of glaucoma and compared findings with human retinal proteomics changes in glaucoma reported previously. We developed an experimental glaucoma rat model by subjecting the rats to increased IOP (9.3 ± 0.1 vs 20.8 ± 1.6 mm Hg) by weekly microbead injections into the eye (8 weeks). The retinal tissues were harvested from control and glaucomatous eyes and protein expression changes analysed using a multiplexed quantitative proteomics approach (TMT-MS3). Immunofluorescence was performed for selected protein markers for data validation. Our study identified 4304 proteins in the rat retinas. Out of these, 139 proteins were downregulated (≤0.83) while the expression of 109 proteins was upregulated (≥1.2-fold change) under glaucoma conditions (P ≤ .05). Computational analysis revealed reduced expression of proteins associated with glutathione metabolism, mitochondrial dysfunction/oxidative phosphorylation, cytoskeleton, and actin filament organisation, along with increased expression of proteins in coagulation cascade, apoptosis, oxidative stress, and RNA processing. Further functional network analysis highlighted the differential modulation of nuclear receptor signalling, cellular survival, protein synthesis, transport, and cellular assembly pathways. Alterations in crystallin family, glutathione metabolism, and mitochondrial dysfunction associated proteins shared similarities between the animal model of glaucoma and the human disease condition. In contrast, the activation of the classical complement pathway and upregulation of cholesterol transport proteins were exclusive to human glaucoma. These findings provide insights into the neurodegenerative mechanisms that are specifically affected in the retina in response to chronically elevated IOP.

An expanded mouse testis transcriptome and mass spectrometry defines novel proteins.

The testis transcriptome is exceptionally complex. Despite its complexity, previous testis transcriptome analyses relied on a reductive method for transcript identification, thus underestimating transcriptome complexity. We describe here a more complete testis transcriptome generated by combining Tuxedo, a reductive method, and spliced-RUM, a combinatorial transcript-building approach. Forty-two percent of the expanded testis transcriptome is composed of unannotated RNAs with novel isoforms of known genes and novel genes constituting 78 and 9.8% of the newly discovered transcripts, respectively. Across tissues, novel transcripts were predominantly expressed in the testis with the exception of novel isoforms which were also highly expressed in the adult ovary. Within the testis, novel isoform expression was distributed equally across all cell types while novel genes were predominantly expressed in meiotic and post-meiotic germ cells. The majority of novel isoforms retained their protein-coding potential while most novel genes had low protein-coding potential. However, a subset of novel genes had protein-coding potentials equivalent to known protein-coding genes. Shotgun mass spectrometry of round spermatid total protein identified unique peptides from four novel genes along with seven annotated non-coding RNAs. These analyses demonstrate the testis expresses a wide range of novel transcripts that give rise to novel proteins.

Pseudomonas aeruginosa Proteome under Hypoxic Stress Conditions Mimicking the Cystic Fibrosis Lung.

Pseudomonas aeruginosa is a ubiquitous Gram-negative pathogen known to inhabit hypoxic mucus plugs of cystic fibrosis (CF) patient lungs. Despite the high prevalence and related patient mortality, the protein machinery enabling the bacterium to adapt to low oxygen environment remains to be fully elucidated. We investigated this by performing both SWATH mass spectrometry and data-dependent SPS-MS3 of TMT-labeled peptides to profile the proteomes of two P. aeruginosa CF isolates, PASS2 and PASS3, and a laboratory reference strain, PAO1, grown under hypoxic stress (O2 < 1%) in media that mimic the nutrient components of the CF lung. Quantitated across all three strains were 3967 P. aeruginosa proteins, reflecting approximately 71% of predicted ORFs in PAO1 and representing the most comprehensive proteome of clinically relevant P. aeruginosa to date. Comparative analysis revealed 735, 640, and 364 proteins were altered by 2-fold or more when comparing low oxygen to aerobic growth in PAO1, PASS2, and PASS3, respectively. Strikingly, under hypoxic stress, all strains showed concurrent increased abundance of proteins required for both aerobic (cbb3-1 and cbb3-2 terminal oxidases) and anaerobic denitrification and arginine fermentation, with the two clinical isolates showing higher relative expression of proteins in these pathways. Additionally, functional annotation revealed that clinical strains portray a unique expression profile of replication, membrane biogenesis, and virulence proteins during hypoxia which may endow these bacteria with a survival advantage. These protein profiles illuminate the diversity of P. aeruginosa mechanisms to adapt to low oxygen and shows that CF isolates initiate a robust molecular response to persist under these conditions.

A multi-parametric flow cytometric assay to analyze DNA-protein interactions.

Interactions between DNA and transcription factors (TFs) guide cellular function and development, yet the complexities of gene regulation are still far from being understood. Such understanding is limited by a paucity of techniques with which to probe DNA-protein interactions. We have devised magnetic protein immobilization on enhancer DNA (MagPIE), a simple, rapid, multi-parametric assay using flow cytometric immunofluorescence to reveal interactions among TFs, chromatin structure and DNA. In MagPIE, synthesized DNA is bound to magnetic beads, which are then incubated with nuclear lysate, permitting sequence-specific binding by TFs, histones and methylation by native lysate factors that can be optionally inhibited with small molecules. Lysate protein-DNA binding is monitored by flow cytometric immunofluorescence, which allows for accurate comparative measurement of TF-DNA affinity. Combinatorial fluorescent staining allows simultaneous analysis of sequence-specific TF-DNA interaction and chromatin modification. MagPIE provides a simple and robust method to analyze complex epigenetic interactions in vitro.

Highly specific intracellular ubiquitination of a small molecule.

Ubiquitin is a small, highly conserved protein that acts as a posttranslational modification in eukaryotes. Ubiquitination of proteins frequently serves as a degradation signal, marking them for disposal by the proteasome. Here, we report a novel small molecule from a diversity-oriented synthesis library, BRD1732, that is directly ubiquitinated in cells, resulting in dramatic accumulation of inactive ubiquitin monomers and polyubiquitin chains causing broad inhibition of the ubiquitin-proteasome system. Ubiquitination of BRD1732 and its associated cytotoxicity are stereospecific and dependent upon two homologous E3 ubiquitin ligases, RNF19A and RNF19B. Our finding opens the possibility for indirect ubiquitination of a target through a ubiquitinated bifunctional small molecule, and more broadly raises the potential for posttranslational modification in trans.

The influence of signals from chilled roots on the proteome of shoot tissues in rice seedlings.

Low root temperature causes a decrease in water uptake, which leads to mineral and nutrient deficiencies with potentially decreased root and shoot growth. Differential temperature effects in plants have been studied extensively, however, the effect of root chilling on the global protein expression in shoots has not been explored. In this study, we imposed chilling temperatures on roots of rice plants while maintaining shoots at optimum atmospheric temperature. Shoot materials (growing zones and leaves) were harvested at five points over a time course of four days, including a two-day recovery period. Proteins were quantified by tandem mass tags and triple stage MS, using a method developed to overcome ratio compression in isobaric-labelled quantitation. Over 3000 proteins in each of the tissues were quantified by multiple peptides. Proteins significantly differentially expressed as compared with the control included abscisic acid-responsive and drought-associated proteins. The data also contained evidence of a possible induction of a sugar signalling pathway.

Network-based inference from complex proteomic mixtures using SNIPE.

MOTIVATION: Proteomics presents the opportunity to provide novel insights about the global biochemical state of a tissue. However, a significant problem with current methods is that shotgun proteomics has limited success at detecting many low abundance proteins, such as transcription factors from complex mixtures of cells and tissues. The ability to assay for these proteins in the context of the entire proteome would be useful in many areas of experimental biology. RESULTS: We used network-based inference in an approach named SNIPE (Software for Network Inference of Proteomics Experiments) that selectively highlights proteins that are more likely to be active but are otherwise undetectable in a shotgun proteomic sample. SNIPE integrates spectral counts from paired case-control samples over a network neighbourhood and assesses the statistical likelihood of enrichment by a permutation test. As an initial application, SNIPE was able to select several proteins required for early murine tooth development. Multiple lines of additional experimental evidence confirm that SNIPE can uncover previously unreported transcription factors in this system. We conclude that SNIPE can enhance the utility of shotgun proteomics data to facilitate the study of poorly detected proteins in complex mixtures. AVAILABILITY AND IMPLEMENTATION: An implementation for the R statistical computing environment named snipeR has been made freely available at http://genetics.bwh.harvard.edu/snipe/. CONTACT: ssunyaev@rics.bwh.harvard.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Liver membrane proteome glycosylation changes in mice bearing an extra-hepatic tumor.

Cancer is well known to be associated with alterations in membrane protein glycosylation (Bird, N. C., Mangnall, D., and Majeed, A. W. (2006) Biology of colorectal liver metastases: A review. J. Surg. Oncol. 94, 68-80; Dimitroff, C. J., Pera, P., Dall'Olio, F., Matta, K. L., Chandrasekaran, E. V., Lau, J. T., and Bernacki, R. J. (1999) Cell surface n-acetylneuraminic acid alpha2,3-galactoside-dependent intercellular adhesion of human colon cancer cells. Biochem. Biophys. Res. Commun. 256, 631-636; and Arcinas, A., Yen, T. Y., Kebebew, E., and Macher, B. A. (2009) Cell surface and secreted protein profiles of human thyroid cancer cell lines reveal distinct glycoprotein patterns. J. Proteome Res. 8, 3958-3968). Equally, it has been well established that tumor-associated inflammation through the release of pro-inflammatory cytokines is a common cause of reduced hepatic drug metabolism and increased toxicity in advanced cancer patients being treated with cytotoxic chemotherapies. However, little is known about the impact of bearing a tumor (and downstream effects like inflammation) on liver membrane protein glycosylation. In this study, proteomic and glycomic analyses were used in combination to determine whether liver membrane protein glycosylation was affected in mice bearing the Engelbreth-Holm Swarm sarcoma. Peptide IPG-IEF and label-free quantitation determined that many enzymes involved in the protein glycosylation pathway specifically; mannosidases (Man1a-I, Man1b-I and Man2a-I), mannoside N-acetylglucosaminyltransferases (Mgat-I and Mgat-II), galactosyltransferases (B3GalT-VII, B4GalT-I, B4GalT-III, C1GalT-I, C1GalT-II, and GalNT-I), and sialyltransferases (ST3Gal-I, ST6Gal-I, and ST6GalNAc-VI) were up-regulated in all livers of tumor-bearing mice (n = 3) compared with nontumor bearing controls (n = 3). In addition, many cell surface lectins: Sialoadhesin-1 (Siglec-1), C-type lectin family 4f (Kupffer cell receptor), and Galactose-binding lectin 9 (Galectin-9) were determined to be up-regulated in the liver of tumor-bearing compared with control mice. Global glycan analysis identified seven N-glycans and two O-glycans that had changed on the liver membrane proteins derived from tumor-bearing mice. Interestingly, α (2,3) sialic acid was found to be up-regulated on the liver membrane of tumor-bearing mice, which reflected the increased expression of its associated sialyltransferase and lectin receptor (siglec-1). The overall increased sialylation on the liver membrane of Engelbreth-Holm Swarm bearing mice correlates with the increased expression of their associated glycosyltransferases and suggests that glycosylation of proteins in the liver plays a role in tumor-induced liver inflammation.

Targeting protein tyrosine phosphatases for CDK6-induced immunotherapy resistance.

Elucidating the mechanisms of resistance to immunotherapy and developing strategies to improve its efficacy are challenging goals. Bioinformatics analysis demonstrates that high CDK6 expression in melanoma is associated with poor progression-free survival of patients receiving single-agent immunotherapy. Depletion of CDK6 or cyclin D3 (but not of CDK4, cyclin D1, or D2) in cells of the tumor microenvironment inhibits tumor growth. CDK6 depletion reshapes the tumor immune microenvironment, and the host anti-tumor effect depends on cyclin D3/CDK6-expressing CD8+ and CD4+ T cells. This occurs by CDK6 phosphorylating and increasing the activities of PTP1B and T cell protein tyrosine phosphatase (TCPTP), which, in turn, decreases tyrosine phosphorylation of CD3ζ, reducing the signal transduction for T cell activation. Administration of a PTP1B and TCPTP inhibitor prove more efficacious than using a CDK6 degrader in enhancing T cell-mediated immunotherapy. Targeting protein tyrosine phosphatases (PTPs) might be an effective strategy for cancer patients who resist immunotherapy treatment.