Widespread variation in molecular interactions and regulatory properties among transcription factor isoforms.

Most human Transcription factors (TFs) genes encode multiple protein isoforms differing in DNA binding domains, effector domains, or other protein regions. The global extent to which this results in functional differences between isoforms remains unknown. Here, we systematically compared 693 isoforms of 246 TF genes, assessing DNA binding, protein binding, transcriptional activation, subcellular localization, and condensate formation. Relative to reference isoforms, two-thirds of alternative TF isoforms exhibit differences in one or more molecular activities, which often could not be predicted from sequence. We observed two primary categories of alternative TF isoforms: "rewirers" and "negative regulators", both of which were associated with differentiation and cancer. Our results support a model wherein the relative expression levels of, and interactions involving, TF isoforms add an understudied layer of complexity to gene regulatory networks, demonstrating the importance of isoform-aware characterization of TF functions and providing a rich resource for further studies.

Paired yeast one-hybrid assays to detect DNA-binding cooperativity and antagonism across transcription factors.

Cooperativity and antagonism between transcription factors (TFs) can drastically modify their binding to regulatory DNA elements. While mapping these relationships between TFs is important for understanding their context-specific functions, existing approaches either rely on DNA binding motif predictions, interrogate one TF at a time, or study individual TFs in parallel. Here, we introduce paired yeast one-hybrid (pY1H) assays to detect cooperativity and antagonism across hundreds of TF-pairs at DNA regions of interest. We provide evidence that a wide variety of TFs are subject to modulation by other TFs in a DNA region-specific manner. We also demonstrate that TF-TF relationships are often affected by alternative isoform usage and identify cooperativity and antagonism between human TFs and viral proteins from human papillomaviruses, Epstein-Barr virus, and other viruses. Altogether, pY1H assays provide a broadly applicable framework to study how different functional relationships affect protein occupancy at regulatory DNA regions.

Scaling-up a fragment-based protein-protein interaction method using a human reference interaction set.

Protein-protein interactions (PPIs) are essential in understanding numerous aspects of protein function. Here, we significantly scaled and modified analyses of the recently developed all-vs-all sequencing (AVA-Seq) approach using a gold-standard human protein interaction set (hsPRS-v2) containing 98 proteins. Binary interaction analyses recovered 20 of 47 (43%) binary PPIs from this positive reference set (PRS), comparing favorably with other methods. However, the increase of 20× in the interaction search space for AVA-Seq analysis in this manuscript resulted in numerous changes to the method required for future use in genome-wide interaction studies. We show that standard sequencing analysis methods must be modified to consider the possible recovery of thousands of positives among millions of tested interactions in a single sequencing run. The PRS data were used to optimize data scaling, auto-activator removal, rank interaction features (such as orientation and unique fragment pairs), and statistical cutoffs. Using these modifications to the method, AVA-Seq recovered >500 known and novel PPIs, including interactions between wild-type fragments of tumor protein p53 and minichromosome maintenance complex proteins 2 and 5 (MCM2 and MCM5) that could be of interest in human disease.

The HTLV-1 viral oncoproteins Tax and HBZ reprogram the cellular mRNA splicing landscape.

Viral infections are known to hijack the transcription and translation of the host cell. However, the extent to which viral proteins coordinate these perturbations remains unclear. Here we used a model system, the human T-cell leukemia virus type 1 (HTLV-1), and systematically analyzed the transcriptome and interactome of key effectors oncoviral proteins Tax and HBZ. We showed that Tax and HBZ target distinct but also common transcription factors. Unexpectedly, we also uncovered a large set of interactions with RNA-binding proteins, including the U2 auxiliary factor large subunit (U2AF2), a key cellular regulator of pre-mRNA splicing. We discovered that Tax and HBZ perturb the splicing landscape by altering cassette exons in opposing manners, with Tax inducing exon inclusion while HBZ induces exon exclusion. Among Tax- and HBZ-dependent splicing changes, we identify events that are also altered in Adult T cell leukemia/lymphoma (ATLL) samples from two independent patient cohorts, and in well-known cancer census genes. Our interactome mapping approach, applicable to other viral oncogenes, has identified spliceosome perturbation as a novel mechanism coordinated by Tax and HBZ to reprogram the transcriptome.

Selenomethionine as an expressible handle for bioconjugations.

Site-selective chemical bioconjugation reactions are enabling tools for the chemical biologist. Guided by a careful study of the selenomethionine (SeM) benzylation, we have refined the reaction to meet the requirements of practical protein bioconjugation. SeM is readily introduced through auxotrophic expression and exhibits unique nucleophilic properties that allow it to be selectively modified even in the presence of cysteine. The resulting benzylselenonium adduct is stable at physiological pH, is selectively labile to glutathione, and embodies a broadly tunable cleavage profile. Specifically, a 4-bromomethylphenylacetyl (BrMePAA) linker has been applied for efficient conjugation of complex organic molecules to SeM-containing proteins. This expansion of the bioconjugation toolkit has broad potential in the development of chemically enhanced proteins.

Comprehensive characterization of protein-protein interactions perturbed by disease mutations.

Technological and computational advances in genomics and interactomics have made it possible to identify how disease mutations perturb protein-protein interaction (PPI) networks within human cells. Here, we show that disease-associated germline variants are significantly enriched in sequences encoding PPI interfaces compared to variants identified in healthy participants from the projects 1000 Genomes and ExAC. Somatic missense mutations are also significantly enriched in PPI interfaces compared to noninterfaces in 10,861 tumor exomes. We computationally identified 470 putative oncoPPIs in a pan-cancer analysis and demonstrate that oncoPPIs are highly correlated with patient survival and drug resistance/sensitivity. We experimentally validate the network effects of 13 oncoPPIs using a systematic binary interaction assay, and also demonstrate the functional consequences of two of these on tumor cell growth. In summary, this human interactome network framework provides a powerful tool for prioritization of alleles with PPI-perturbing mutations to inform pathobiological mechanism- and genotype-based therapeutic discovery.

A genome-wide positioning systems network algorithm for in silico drug repurposing.

Recent advances in DNA/RNA sequencing have made it possible to identify new targets rapidly and to repurpose approved drugs for treating heterogeneous diseases by the 'precise' targeting of individualized disease modules. In this study, we develop a Genome-wide Positioning Systems network (GPSnet) algorithm for drug repurposing by specifically targeting disease modules derived from individual patient's DNA and RNA sequencing profiles mapped to the human protein-protein interactome network. We investigate whole-exome sequencing and transcriptome profiles from ~5,000 patients across 15 cancer types from The Cancer Genome Atlas. We show that GPSnet-predicted disease modules can predict drug responses and prioritize new indications for 140 approved drugs. Importantly, we experimentally validate that an approved cardiac arrhythmia and heart failure drug, ouabain, shows potential antitumor activities in lung adenocarcinoma by uniquely targeting a HIF1α/LEO1-mediated cell metabolism pathway. In summary, GPSnet offers a network-based, in silico drug repurposing framework for more efficacious therapeutic selections.

Network Analysis of UBE3A/E6AP-Associated Proteins Provides Connections to Several Distinct Cellular Processes.

Perturbations in activity and dosage of the UBE3A ubiquitin-ligase have been linked to Angelman syndrome and autism spectrum disorders. UBE3A was initially identified as the cellular protein hijacked by the human papillomavirus E6 protein to mediate the ubiquitylation of p53, a function critical to the oncogenic potential of these viruses. Although a number of substrates have been identified, the normal cellular functions and pathways affected by UBE3A are largely unknown. Previously, we showed that UBE3A associates with HERC2, NEURL4, and MAPK6/ERK3 in a high-molecular-weight complex of unknown function that we refer to as the HUN complex (HERC2, UBE3A, and NEURL4). In this study, the combination of two complementary proteomic approaches with a rigorous network analysis revealed cellular functions and pathways in which UBE3A and the HUN complex are involved. In addition to finding new UBE3A-associated proteins, such as MCM6, SUGT1, EIF3C, and ASPP2, network analysis revealed that UBE3A-associated proteins are connected to several fundamental cellular processes including translation, DNA replication, intracellular trafficking, and centrosome regulation. Our analysis suggests that UBE3A could be involved in the control and/or integration of these cellular processes, in some cases as a component of the HUN complex, and also provides evidence for crosstalk between the HUN complex and CAMKII interaction networks. This study contributes to a deeper understanding of the cellular functions of UBE3A and its potential role in pathways that may be affected in Angelman syndrome, UBE3A-associated autism spectrum disorders, and human papillomavirus-associated cancers.

Proteome-scale Binary Interactomics in Human Cells.

Because proteins are the main mediators of most cellular processes they are also prime therapeutic targets. Identifying physical links among proteins and between drugs and their protein targets is essential in order to understand the mechanisms through which both proteins themselves and the molecules they are targeted with act. Thus, there is a strong need for sensitive methods that enable mapping out these biomolecular interactions. Here we present a robust and sensitive approach to screen proteome-scale collections of proteins for binding to proteins or small molecules using the well validated MAPPIT (Mammalian Protein-Protein Interaction Trap) and MASPIT (Mammalian Small Molecule-Protein Interaction Trap) assays. Using high-density reverse transfected cell microarrays, a close to proteome-wide collection of human ORF clones can be screened for interactors at high throughput. The versatility of the platform is demonstrated through several examples. With MAPPIT, we screened a 15k ORF library for binding partners of RNF41, an E3 ubiquitin protein ligase implicated in receptor sorting, identifying known and novel interacting proteins. The potential related to the fact that MAPPIT operates in living human cells is illustrated in a screen where the protein collection is scanned for interactions with the glucocorticoid receptor (GR) in its unliganded versus dexamethasone-induced activated state. Several proteins were identified the interaction of which is modulated upon ligand binding to the GR, including a number of previously reported GR interactors. Finally, the screening technology also enables detecting small molecule target proteins, which in many drug discovery programs represents an important hurdle. We show the efficiency of MASPIT-based target profiling through screening with tamoxifen, a first-line breast cancer drug, and reversine, an investigational drug with interesting dedifferentiation and antitumor activity. In both cases, cell microarray screens yielded known and new potential drug targets highlighting the utility of the technology beyond fundamental biology.

The transcription factor ERG recruits CCR4-NOT to control mRNA decay and mitotic progression.

Control of mRNA levels, a fundamental aspect in the regulation of gene expression, is achieved through a balance between mRNA synthesis and decay. E26-related gene (Erg) proteins are canonical transcription factors whose previously described functions are confined to the control of mRNA synthesis. Here, we report that ERG also regulates gene expression by affecting mRNA stability and identify the molecular mechanisms underlying this function in human cells. ERG is recruited to mRNAs via interaction with the RNA-binding protein RBPMS, and it promotes mRNA decay by binding CNOT2, a component of the CCR4-NOT deadenylation complex. Transcriptome-wide mRNA stability analysis revealed that ERG controls the degradation of a subset of mRNAs highly connected to Aurora signaling, whose decay during S phase is necessary for mitotic progression. Our data indicate that control of gene expression by mammalian transcription factors may follow a more complex scheme than previously anticipated, integrating mRNA synthesis and degradation.

MECP2 Is a Frequently Amplified Oncogene with a Novel Epigenetic Mechanism That Mimics the Role of Activated RAS in Malignancy.

UNLABELLED: An unbiased genome-scale screen for unmutated genes that drive cancer growth when overexpressed identified methyl cytosine-guanine dinucleotide (CpG) binding protein 2 (MECP2) as a novel oncogene. MECP2 resides in a region of the X-chromosome that is significantly amplified across 18% of cancers, and many cancer cell lines have amplified, overexpressed MECP2 and are dependent on MECP2 expression for growth. MECP2 copy-number gain and RAS family member alterations are mutually exclusive in several cancer types. The MECP2 splicing isoforms activate the major growth factor pathways targeted by activated RAS, the MAPK and PI3K pathways. MECP2 rescued the growth of a KRAS(G12C)-addicted cell line after KRAS downregulation, and activated KRAS rescues the growth of an MECP2-addicted cell line after MECP2 downregulation. MECP2 binding to the epigenetic modification 5-hydroxymethylcytosine is required for efficient transformation. These observations suggest that MECP2 is a commonly amplified oncogene with an unusual epigenetic mode of action. SIGNIFICANCE: MECP2 is a commonly amplified oncogene in human malignancies with a unique epigenetic mechanism of action. Cancer Discov; 6(1); 45-58. ©2015 AACR.This article is highlighted in the In This Issue feature, p. 1.

Protein domain-level landscape of cancer-type-specific somatic mutations.

Identifying driver mutations and their functional consequences is critical to our understanding of cancer. Towards this goal, and because domains are the functional units of a protein, we explored the protein domain-level landscape of cancer-type-specific somatic mutations. Specifically, we systematically examined tumor genomes from 21 cancer types to identify domains with high mutational density in specific tissues, the positions of mutational hotspots within these domains, and the functional and structural context where possible. While hotspots corresponding to specific gain-of-function mutations are expected for oncoproteins, we found that tumor suppressor proteins also exhibit strong biases toward being mutated in particular domains. Within domains, however, we observed the expected patterns of mutation, with recurrently mutated positions for oncogenes and evenly distributed mutations for tumor suppressors. For example, we identified both known and new endometrial cancer hotspots in the tyrosine kinase domain of the FGFR2 protein, one of which is also a hotspot in breast cancer, and found new two hotspots in the Immunoglobulin I-set domain in colon cancer. Thus, to prioritize cancer mutations for further functional studies aimed at more precise cancer treatments, we have systematically correlated mutations and cancer types at the protein domain level.

Systematic screening reveals a role for BRCA1 in the response to transcription-associated DNA damage.

BRCA1 is a breast and ovarian tumor suppressor. Given its numerous incompletely understood functions and the possibility that more exist, we performed complementary systematic screens in search of new BRCA1 protein-interacting partners. New BRCA1 functions and/or a better understanding of existing ones were sought. Among the new interacting proteins identified, genetic interactions were detected between BRCA1 and four of the interactors: TONSL, SETX, TCEANC, and TCEA2. Genetic interactions were also detected between BRCA1 and certain interactors of TONSL, including both members of the FACT complex. From these results, a new BRCA1 function in the response to transcription-associated DNA damage was detected. Specifically, new roles for BRCA1 in the restart of transcription after UV damage and in preventing or repairing damage caused by stabilized R loops were identified. These roles are likely carried out together with some of the newly identified interactors. This new function may be important in BRCA1 tumor suppression, since the expression of several interactors, including some of the above-noted transcription proteins, is repeatedly aberrant in both breast and ovarian cancers.

Systematic identification of pathological lamin A interactors.

Laminopathies are a collection of phenotypically diverse diseases that include muscular dystrophies, cardiomyopathies, lipodystrophies, and premature aging syndromes. Laminopathies are caused by >300 distinct mutations in the LMNA gene, which encodes the nuclear intermediate filament proteins lamin A and C, two major architectural elements of the mammalian cell nucleus. The genotype-phenotype relationship and the basis for the pronounced tissue specificity of laminopathies are poorly understood. Here we seek to identify on a global scale lamin A-binding partners whose interaction is affected by disease-relevant LMNA mutations. In a screen of a human genome-wide ORFeome library, we identified and validated 337 lamin A-binding proteins. Testing them against 89 known lamin A disease mutations identified 50 disease-associated interactors. Association of progerin, the lamin A isoform responsible for the premature aging disorder Hutchinson-Gilford progeria syndrome, with its partners was largely mediated by farnesylation. Mapping of the interaction sites on lamin A identified the immunoglobulin G (IgG)-like domain as an interaction hotspot and demonstrated that lamin A variants, which destabilize the Ig-like domain, affect protein-protein interactions more globally than mutations of surface residues. Analysis of a set of LMNA mutations in a single residue, which result in three phenotypically distinct diseases, identified disease-specific interactors. The results represent a systematic map of disease-relevant lamin A interactors and suggest loss of tissue-specific lamin A interactions as a mechanism for the tissue-specific appearance of laminopathic phenotypes.

An RS motif within the Epstein-Barr virus BLRF2 tegument protein is phosphorylated by SRPK2 and is important for viral replication.

Epstein-Barr virus (EBV) is a gammaherpesvirus that causes infectious mononucleosis, B cell lymphomas, and nasopharyngeal carcinoma. Many of the genes required for EBV virion morphogenesis are found in all herpesviruses, but some are specific to gammaherpesviruses. One of these gamma-specific genes, BLRF2, encodes a tegument protein that has been shown to be essential for replication in other gammaherpesviruses. In this study, we identify BLRF2 interacting proteins using binary and co-complex protein assays. Serine/Arginine-rich Protein Kinase 2 (SRPK2) was identified by both assays and was further shown to phosphorylate an RS motif in the BLRF2 C-terminus. Mutation of this RS motif (S148A+S150A) abrogated the ability of BLRF2 to support replication of a murine gammaherpesvirus 68 genome lacking the BLRF2 homolog (ORF52). We conclude that the BLRF2 RS motif is phosphorylated by SRPK2 and is important for viral replication.

Identification of FAM111A as an SV40 host range restriction and adenovirus helper factor.

The small genome of polyomaviruses encodes a limited number of proteins that are highly dependent on interactions with host cell proteins for efficient viral replication. The SV40 large T antigen (LT) contains several discrete functional domains including the LXCXE or RB-binding motif, the DNA binding and helicase domains that contribute to the viral life cycle. In addition, the LT C-terminal region contains the host range and adenovirus helper functions required for lytic infection in certain restrictive cell types. To understand how LT affects the host cell to facilitate viral replication, we expressed full-length or functional domains of LT in cells, identified interacting host proteins and carried out expression profiling. LT perturbed the expression of p53 target genes and subsets of cell-cycle dependent genes regulated by the DREAM and the B-Myb-MuvB complexes. Affinity purification of LT followed by mass spectrometry revealed a specific interaction between the LT C-terminal region and FAM111A, a previously uncharacterized protein. Depletion of FAM111A recapitulated the effects of heterologous expression of the LT C-terminal region, including increased viral gene expression and lytic infection of SV40 host range mutants and adenovirus replication in restrictive cells. FAM111A functions as a host range restriction factor that is specifically targeted by SV40 LT.

Comparative analysis of virus-host interactomes with a mammalian high-throughput protein complementation assay based on Gaussia princeps luciferase.

Comparative interactomics is a strategy for inferring potential interactions among orthologous proteins or "interologs". Herein we focus, in contrast to standard homology-based inference, on the divergence of protein interaction profiles among closely related organisms, showing that the approach can correlate specific traits to phenotypic differences. As a model, this new comparative interactomic approach was applied at a large scale to human papillomaviruses (HPVs) proteins. The oncogenic potential of HPVs is mainly determined by the E6 and E7 early proteins. We have mapped and overlapped the virus-host protein interaction networks of E6 and E7 proteins from 11 distinct HPV genotypes, selected for their different tropisms and pathologies. We generated robust and comprehensive datasets by combining two orthogonal protein interaction assays: yeast two-hybrid (Y2H), and our recently described "high-throughput Gaussia princeps protein complementation assay" (HT-GPCA). HT-GPCA detects protein interaction by measuring the interaction-mediated reconstitution of activity of a split G. princeps luciferase. Hierarchical clustering of interaction profiles recapitulated HPV phylogeny and was used to correlate specific virus-host interaction profiles with pathological traits, reflecting the distinct carcinogenic potentials of different HPVs. This comparative interactomics constitutes a reliable and powerful strategy to decipher molecular relationships in virtually any combination of microorganism-host interactions.

Viral perturbations of host networks reflect disease etiology.

Many human diseases, arising from mutations of disease susceptibility genes (genetic diseases), are also associated with viral infections (virally implicated diseases), either in a directly causal manner or by indirect associations. Here we examine whether viral perturbations of host interactome may underlie such virally implicated disease relationships. Using as models two different human viruses, Epstein-Barr virus (EBV) and human papillomavirus (HPV), we find that host targets of viral proteins reside in network proximity to products of disease susceptibility genes. Expression changes in virally implicated disease tissues and comorbidity patterns cluster significantly in the network vicinity of viral targets. The topological proximity found between cellular targets of viral proteins and disease genes was exploited to uncover a novel pathway linking HPV to Fanconi anemia.

Interpreting cancer genomes using systematic host network perturbations by tumour virus proteins.

Genotypic differences greatly influence susceptibility and resistance to disease. Understanding genotype-phenotype relationships requires that phenotypes be viewed as manifestations of network properties, rather than simply as the result of individual genomic variations. Genome sequencing efforts have identified numerous germline mutations, and large numbers of somatic genomic alterations, associated with a predisposition to cancer. However, it remains difficult to distinguish background, or 'passenger', cancer mutations from causal, or 'driver', mutations in these data sets. Human viruses intrinsically depend on their host cell during the course of infection and can elicit pathological phenotypes similar to those arising from mutations. Here we test the hypothesis that genomic variations and tumour viruses may cause cancer through related mechanisms, by systematically examining host interactome and transcriptome network perturbations caused by DNA tumour virus proteins. The resulting integrated viral perturbation data reflects rewiring of the host cell networks, and highlights pathways, such as Notch signalling and apoptosis, that go awry in cancer. We show that systematic analyses of host targets of viral proteins can identify cancer genes with a success rate on a par with their identification through functional genomics and large-scale cataloguing of tumour mutations. Together, these complementary approaches increase the specificity of cancer gene identification. Combining systems-level studies of pathogen-encoded gene products with genomic approaches will facilitate the prioritization of cancer-causing driver genes to advance the understanding of the genetic basis of human cancer.

Host-pathogen interactome mapping for HTLV-1 and -2 retroviruses.

BACKGROUND: Human T-cell leukemia virus type 1 (HTLV-1) and type 2 both target T lymphocytes, yet induce radically different phenotypic outcomes. HTLV-1 is a causative agent of Adult T-cell leukemia (ATL), whereas HTLV-2, highly similar to HTLV-1, causes no known overt disease. HTLV gene products are engaged in a dynamic struggle of activating and antagonistic interactions with host cells. Investigations focused on one or a few genes have identified several human factors interacting with HTLV viral proteins. Most of the available interaction data concern the highly investigated HTLV-1 Tax protein. Identifying shared and distinct host-pathogen protein interaction profiles for these two viruses would enlighten how they exploit distinctive or common strategies to subvert cellular pathways toward disease progression. RESULTS: We employ a scalable methodology for the systematic mapping and comparison of pathogen-host protein interactions that includes stringent yeast two-hybrid screening and systematic retest, as well as two independent validations through an additional protein interaction detection method and a functional transactivation assay. The final data set contained 166 interactions between 10 viral proteins and 122 human proteins. Among the 166 interactions identified, 87 and 79 involved HTLV-1 and HTLV-2 -encoded proteins, respectively. Targets for HTLV-1 and HTLV-2 proteins implicate a diverse set of cellular processes including the ubiquitin-proteasome system, the apoptosis, different cancer pathways and the Notch signaling pathway. CONCLUSIONS: This study constitutes a first pass, with homogeneous data, at comparative analysis of host targets for HTLV-1 and -2 retroviruses, complements currently existing data for formulation of systems biology models of retroviral induced diseases and presents new insights on biological pathways involved in retroviral infection.