The acidic transcription activator Gcn4 binds the mediator subunit Gal11/Med15 using a simple protein interface forming a fuzzy complex

Glycosylation Weakens Skp1 Homodimerization in Toxoplasma gondii by Interrupting a Fuzzy Interaction

Skp1/Cullin-1/F-Box protein (SCF) complexes represent a major class of E3 ubiquitin ligases responsible for proteomic control throughout eukaryotes. Target specificity is mediated by a large set of F-box proteins (FBPs) whose F-box domains interact with Skp1 in a conserved, well-organized fashion. In the social amoeba Dictyostelium, Skp1 is regulated by oxygen-dependent glycosylation which alters Skp1's FBP interactome and inhibits homodimerization that is mediated in part by an ordered interface which overlaps with that of FBPs. Based on sedimentation velocity experiments, Skp1 from the intracellular pathogen Toxoplasma gondii exhibits a homodimerization Kd comparable to that of a previously measured FBP/Skp1 interaction. Glycosylation of Skp1's disordered C-terminal region (CTR) distal to the ordered homodimer interface significantly weakens Skp1 homodimerization, an effect reproduced by CTR deletion. Replacement with a randomized CTR sequence retains high affinity excluding an extension of the ordered dimer interface. Substitution by poly serine weakens the homodimer to a degree equal to its deletion, indicating a composition dependent effect. The contribution of the CTR to Skp1 homodimerization is canceled by high salt consistent with an electrostatic mechanism. All-atom molecular dynamics simulations suggest that the CTR promotes homodimerization via charge cluster interactions. Taken together, the data indicate that glycosylation weakens homodimerization by disrupting a C-terminal fuzzy interaction that functions in tandem with the ordered dimer interface, thereby freeing Skp1 for FBP binding. Thus, the CTR contributes to Skp1/Skp1 and Skp1/FBP interactions via independent mechanisms that are each influenced by O2, indicating multiple constraints on the evolution of its sequence.

The mechanism of YAP/TAZ transactivation and dual targeting for cancer therapy

Transcriptional coactivators Yes-associated protein (YAP) and transcriptional coactivator with PDZ-binding motif (TAZ) play key roles in cancers through transcriptional outputs. However, their transactivation mechanisms remain unclear, and effective targeting strategies are lacking. Here, we show that YAP/TAZ possess a hydrophobic transactivation domain (TAD). TAD knockout prevents tumor establishment due to growth defects and enhances immune attack. Mechanistically, TADs facilitate preinitiation complex (PIC) assembly by recruiting the TATA-binding protein-associated factor 4 (TAF4)-dependent TFIID complex and enhance RNA polymerase II (Pol II) elongation through mediator complex subunit 15 (MED15)-dependent mediator recruitment for the expressions of oncogenic/immune-suppressive programs. The synthesized peptide TJ-M11 selectively disrupts TAD interactions with MED15 and TAF4, suppressing tumor growth and sensitizing tumors to immunotherapy. Our findings demonstrate that YAP/TAZ TADs exhibit dual functions in PIC assembly and Pol II elongation via hydrophobic interactions, which represent actionable targets for cancer therapy and combination immunotherapy.

Structural basis of human Mediator recruitment by the phosphorylated transcription factor Elk-1

One function of Mediator complex subunit MED23 is to mediate transcriptional activation by the phosphorylated transcription factor Elk-1, in response to the Ras-MAPK signaling pathway. Using cryogenic electron microscopy, we solve a 3.0 Å structure of human MED23 complexed with the phosphorylated activation domain of Elk-1. Elk-1 binds to MED23 via a hydrophobic sequence PSIHFWSTLSPP containing one phosphorylated residue (S383p), which forms a tight turn around the central Phenylalanine. Binding of Elk-1 induces allosteric changes in MED23 that propagate to the opposite face of the subunit, resulting in the dynamic behavior of a 19-residue segment, which alters the molecular surface of MED23. We design a specific MED23 mutation (G382F) that disrupts Elk--1 binding and consequently impairs Elk-1-dependent serum-induced activation of target genes in the Ras-Raf-MEK-ERK signaling pathway. The structure provides molecular details and insights into a Mediator subunit-transcription factor interface.

Low overlap of transcription factor DNA binding and regulatory targets

DNA sequence-specific transcription factors (TFs) modulate transcription and chromatin architecture, acting from regulatory sites in enhancers and promoters of eukaryotic genes1,2. How multiple TFs cooperate to regulate individual genes is still unclear. In yeast, most TFs are thought to regulate transcription via binding to upstream activating sequences, which are situated within a few hundred base pairs upstream of the regulated gene3. Although this model has been validated for individual TFs and specific genes, it has not been tested in a systematic way. Here we integrated information on the binding and expression targets for the near-complete set of yeast TFs and show that, contrary to expectations, there are few TFs with dedicated activator or repressor roles, and that most TFs have a dual function. Although nearly all protein-coding genes are regulated by one or more TFs, our analysis revealed limited overlap between TF binding and gene regulation. Rapid depletion of many TFs also revealed many regulatory targets that were distant from detectable TF binding sites, suggesting unexpected regulatory mechanisms. Our study provides a comprehensive survey of TF functions and offers insights into interactions between the set of TFs expressed in a single cell type and how they contribute to the complex programme of gene regulation.

Role of charges in a dynamic disordered complex between an IDP and a folded domain

Protein complexes involving intrinsically disordered proteins (IDPs) cover a continuum from IDPs that fully fold upon binding to IDPs that remain fully disordered in the complex. Here we demonstrate a case of charge-driven interactions of a folded domain with an oppositely charged IDP that remains completely disordered in the complex. Using the negatively charged and fully disordered prothymosin α and the positively charged and folded globular domain of histone H1.0, we show that they form a low-micromolar-affinity complex without fixed relative orientations or persistent contacts between specific residues. Using 25 charge variants of the globular domain, we find that the binding affinity can be modulated both by net charge and charge clustering on the folded domain, indicating some selectivity in highly charged complexes. Our results highlight that a folded protein can provide a charged surface onto which an oppositely charged IDP can bind while retaining disorder. We expect that more such complexes exist.

Boosting transcriptional activities by employing repeated activation domains in transcription factors

Enhancing the transcriptional activation activity of transcription factors (TFs) has multiple applications in organism improvement, metabolic engineering, and other aspects of plant science, but the approaches remain unclear. Here, we used gene activation assays and genetic transformation to investigate the transcriptional activities of two MYB TFs, PRODUCTION OF ANTHOCYANIN PIGMENT 1 (AtPAP1) from Arabidopsis (Arabidopsis thaliana) and EsMYBA1 from Epimedium (Epimedium sagittatum), and their synthetic variants in a range of plant species from several families. Using anthocyanin biosynthesis as a convenient readout, we discovered that homologous naturally occurring TFs showed differences in the transcriptional activation ability and that similar TFs induced large changes in the genetic program when heterologously expressed in different species. In some cases, shuffling the DNA-binding domains and transcriptional activation domains (ADs) between homologous TFs led to synthetic TFs that had stronger activation potency than the original TFs. More importantly, synthetic TFs derived from MYB, NAC, bHLH, and ethylene-insensitive3-like (EIL) family members containing tandemly repeated ADs had greatly enhanced activity compared to their natural counterparts. These findings enhance our understanding of TF activity and demonstrate that employing tandemly repeated ADs from natural TFs is a simple and widely applicable strategy to enhance the activation potency of synthetic TFs.

The Role of Med15 Sequence Features in Transcription Factor Interactions

Med15 is a general transcriptional regulator and tail module subunit within the RNA Pol II mediator complex. The Saccharomyces cerevisiae Med15 protein has a well-structured N-terminal KIX domain, three activator binding domains (ABDs) and several naturally variable polyglutamine (poly-Q) tracts (Q1, Q2, Q3) embedded in an intrinsically disordered central region, and a C-terminal mediator association domain (MAD). We investigated how the presence of ABDs and changes in length and composition of poly-Q tracts influences Med15 activity using phenotypic, gene expression, transcription factor interaction and phase separation assays of truncation, deletion, and synthetic alleles. We found that individual Med15 activities were influenced by the number of activator binding domains (ABDs) and adjacent polyglutamine tract composition. Robust Med15 activity required at least the Q1 tract and the length of that tract modulated activity in a context-dependent manner. Reduced Msn2-dependent transcriptional activation due to Med15 Q1 tract variation correlated with reduced Msn2:Med15 interaction strength, but interaction strength did not always mirror phase separation propensity. We also observed that distant glutamine tracts and Med15 phosphorylation affected the activities of the KIX domain, and interaction studies revealed that intramolecular interactions may affect some Med15-transcription factor interactions.

Mechanisms of synergistic Mediator recruitment in RNA polymerase II transcription activation revealed by single-molecule fluorescence

Transcription activators trigger transcript production by RNA Polymerase II (RNApII) via the Mediator coactivator complex. Here the dynamics of activator, Mediator, and RNApII binding at promoter DNA were analyzed using multi-wavelength single-molecule microscopy of fluorescently labeled proteins in budding yeast nuclear extract. Binding of Mediator and RNApII to the template required activator and an upstream activator sequence (UAS), but not a core promoter. While Mediator and RNApII sometimes bind as a pre-formed complex, more commonly Mediator binds first and subsequently recruits RNApII to form a preinitiation complex precursor (pre-PIC) tethered to activators on the UAS. Interestingly, Mediator occupancy has a highly non-linear response to activator concentration, and fluorescence intensity measurements show Mediator preferentially associates with templates having at least two activators bound. Statistical mechanical modeling suggests this "synergy" is not due to cooperative binding between activators, but instead occurs when multiple DNA-bound activator molecules simultaneously interact with a single Mediator.

Conservation of function without conservation of amino acid sequence in intrinsically disordered transcriptional activation domains

Protein function is canonically believed to be more conserved than amino acid sequence, but this idea is only well supported in folded domains, where highly diverged sequences can fold into equivalent 3D structures. In contrast, intrinsically disordered protein regions (IDRs) do not fold into a stable 3D structure, thus it remains unknown when and how function is conserved for IDRs that experience rapid amino acid sequence divergence. As a model system for studying the evolution of IDRs, we examined transcriptional activation domains, the regions of transcription factors that bind to coactivator complexes. We systematically identified activation domains on 502 orthologs of the transcriptional activator Gcn4 spanning 600 MY of fungal evolution. We find that the central activation domain shows strong conservation of function without conservation of sequence. This conservation of function without conservation of sequence is facilitated by evolutionary turnover (gain and loss) of key acidic and aromatic residues, the positions most important for function. This high sequence flexibility of functional orthologs mirrors the physical flexibility of the activation domain coactivator interaction interface, suggesting that physical flexibility enables evolutionary plasticity. We propose that turnover of short functional elements, sometimes individual amino acids, is a general mechanism for conservation of function without conservation of sequence during IDR evolution.

Fuzzy protein-DNA interactions and beyond: A common theme in transcription?

Gene expression regulation requires both diversity and specificity. How can these two contradictory conditions be reconciled? Dynamic DNA recognition mechanisms lead to heterogeneous bound conformations, which can be shifted by the cellular cues. Here we summarise recent experimental evidence on how fuzzy interactions contribute to chromatin remodelling, regulation of DNA replication and repair and transcription factor binding. We describe how the binding mode continuum between DNA and regulatory factors lead to variable, multisite contact patterns; polyelectrolyte competitions; on-the-fly shape readouts; autoinhibition controlled by posttranslational modifications or dynamic oligomerisation mechanisms. Increasing experimental evidence supports the rugged energy landscape of the bound protein-DNA assembly, modulation of which leads to distinct functional outcomes. Recent results suggest the evolutionary conservation of these combinatorial mechanisms with moderate sequence constraints in the malleable transcriptional machinery.

Reconciling competing models on the roles of condensates and soluble complexes in transcription factor function

Phase separation explains the exquisite spatial and temporal regulation of many biological processes, but the role of transcription factor-mediated condensates in gene regulation is contentious, requiring head-to-head comparison of competing models. Here, we focused on the prototypical yeast transcription factor Gcn4 and assessed two models for gene transcription activation, i.e., mediated via soluble complexes or transcriptional condensates. Both models rely on the ability of transcription factors and coactivators to engage in multivalent interactions. Unexpectedly, we found that propensity to form homotypic Gcn4 condensates does not correlate well with transcriptional activity. Contrary to prevailing models, binding to DNA suppresses Gcn4 phase separation. Notably, the ability of Gcn4 to form soluble complexes with coactivator subunit Med15 closely mirrored the propensity to recruit Med15 into condensates, indicating that these properties are intertwined and cautioning against interpretation of mutational data without head-to-head comparisons. However, Gcn4 variants with the highest affinity for Med15 do not function as well as expected and instead have activities that reflect their abilities to phase separate with Med15. These variants therefore indeed form cellular condensates, and those attenuate activity. Our results show that transcription factors can function as soluble complexes as well as condensates, reconciling two seemingly opposing models, and have implications for other phase-separating systems.

Revisiting the model for coactivator recruitment: Med15 can select its target sites independent of promoter-bound transcription factors

Activation domains (ADs) within transcription factors (TFs) induce gene expression by recruiting coactivators such as the Mediator complex. Coactivators lack DNA binding domains (DBDs) and are assumed to passively follow their recruiting TFs. This is supported by direct AD-coactivator interactions seen in vitro but has not yet been tested in living cells. To examine that, we targeted two Med15-recruiting ADs to a range of budding yeast promoters through fusion with different DBDs. The DBD-AD fusions localized to hundreds of genomic sites but recruited Med15 and induced transcription in only a subset of bound promoters, characterized by a fuzzy-nucleosome architecture. Direct DBD-Med15 fusions shifted DBD localization towards fuzzy-nucleosome promoters, including promoters devoid of the endogenous Mediator. We propose that Med15, and perhaps other coactivators, possess inherent promoter preference and thus actively contribute to the selection of TF-induced genes.

Context transcription factors establish cooperative environments and mediate enhancer communication

Many enhancers control gene expression by assembling regulatory factor clusters, also referred to as condensates. This process is vital for facilitating enhancer communication and establishing cellular identity. However, how DNA sequence and transcription factor (TF) binding instruct the formation of high regulatory factor environments remains poorly understood. Here we developed a new approach leveraging enhancer-centric chromatin accessibility quantitative trait loci (caQTLs) to nominate regulatory factor clusters genome-wide. By analyzing TF-binding signatures within the context of caQTLs and comparing episomal versus endogenous enhancer activities, we discovered a class of regulators, 'context-only' TFs, that amplify the activity of cell type-specific caQTL-binding TFs, that is, 'context-initiator' TFs. Similar to super-enhancers, enhancers enriched for context-only TF-binding sites display high coactivator binding and sensitivity to bromodomain-inhibiting molecules. We further show that binding sites for context-only and context-initiator TFs underlie enhancer coordination, providing a mechanistic rationale for how a loose TF syntax confers regulatory specificity.

High-throughput affinity measurements of direct interactions between activation domains and co-activators

Sequence-specific activation by transcription factors is essential for gene regulation1,2. Key to this are activation domains, which often fall within disordered regions of transcription factors3,4 and recruit co-activators to initiate transcription5. These interactions are difficult to characterize via most experimental techniques because they are typically weak and transient6,7. Consequently, we know very little about whether these interactions are promiscuous or specific, the mechanisms of binding, and how these interactions tune the strength of gene activation. To address these questions, we developed a microfluidic platform for expression and purification of hundreds of activation domains in parallel followed by direct measurement of co-activator binding affinities (STAMMPPING, for Simultaneous Trapping of Affinity Measurements via a Microfluidic Protein-Protein INteraction Generator). By applying STAMMPPING to quantify direct interactions between eight co-activators and 204 human activation domains (>1,500 K ds), we provide the first quantitative map of these interactions and reveal 334 novel binding pairs. We find that the metazoan-specific co-activator P300 directly binds >100 activation domains, potentially explaining its widespread recruitment across the genome to influence transcriptional activation. Despite sharing similar molecular properties (e.g. enrichment of negative and hydrophobic residues), activation domains utilize distinct biophysical properties to recruit certain co-activator domains. Co-activator domain affinity and occupancy are well-predicted by analytical models that account for multivalency, and in vitro affinities quantitatively predict activation in cells with an ultrasensitive response. Not only do our results demonstrate the ability to measure affinities between even weak protein-protein interactions in high throughput, but they also provide a necessary resource of over 1,500 activation domain/co-activator affinities which lays the foundation for understanding the molecular basis of transcriptional activation.

An autoinhibitory switch of the LSD1 disordered region controls enhancer silencing

Transcriptional coregulators and transcription factors (TFs) contain intrinsically disordered regions (IDRs) that are critical for their association and function in gene regulation. More recently, IDRs have been shown to promote multivalent protein-protein interactions between coregulators and TFs to drive their association into condensates. By contrast, here we demonstrate how the IDR of the corepressor LSD1 excludes TF association, acting as a dynamic conformational switch that tunes repression of active cis-regulatory elements. Hydrogen-deuterium exchange shows that the LSD1 IDR interconverts between transient open and closed conformational states, the latter of which inhibits partitioning of the protein's structured domains with TF condensates. This autoinhibitory switch controls leukemic differentiation by modulating repression of active cis-regulatory elements bound by LSD1 and master hematopoietic TFs. Together, these studies unveil alternative mechanisms by which disordered regions and their dynamic crosstalk with structured regions can shape coregulator-TF interactions to control cis-regulatory landscapes and cell fate.

Direct prediction of intermolecular interactions driven by disordered regions

Intrinsically disordered regions (IDRs) are critical for a wide variety of cellular functions, many of which involve interactions with partner proteins. Molecular recognition is typically considered through the lens of sequence-specific binding events. However, a growing body of work has shown that IDRs often interact with partners in a manner that does not depend on the precise order of the amino acid order, instead driven by complementary chemical interactions leading to disordered bound-state complexes. Despite this emerging paradigm, we lack tools to describe, quantify, predict, and interpret these types of structurally heterogeneous interactions from the underlying amino acid sequences. Here, we repurpose the chemical physics developed originally for molecular simulations to develop an approach for predicting intermolecular interactions between IDRs and partner proteins. Our approach enables the direct prediction of phase diagrams, the identification of chemically-specific interaction hotspots on IDRs, and a route to develop and test mechanistic hypotheses regarding IDR function in the context of molecular recognition. We use our approach to examine a range of systems and questions to highlight its versatility and applicability.

Fuzzy recognition by the prokaryotic transcription factor HigA2 from Vibrio cholerae

Disordered protein sequences can exhibit different binding modes, ranging from well-ordered folding-upon-binding to highly dynamic fuzzy binding. The primary function of the intrinsically disordered region of the antitoxin HigA2 from Vibrio cholerae is to neutralize HigB2 toxin through ultra-high-affinity folding-upon-binding interaction. Here, we show that the same intrinsically disordered region can also mediate fuzzy interactions with its operator DNA and, through interplay with the folded helix-turn-helix domain, regulates transcription from the higBA2 operon. NMR, SAXS, ITC and in vivo experiments converge towards a consistent picture where a specific set of residues in the intrinsically disordered region mediate electrostatic and hydrophobic interactions while "hovering" over the DNA operator. Sensitivity of the intrinsically disordered region to scrambling the sequence, position-specific contacts and absence of redundant, multivalent interactions, point towards a more specific type of fuzzy binding. Our work demonstrates how a bacterial regulator achieves dual functionality by utilizing two distinct interaction modes within the same disordered sequence.

Assessment of machine-learning predictions for the Mediator complex subunit MED25 ACID domain interactions with transactivation domains

The human Mediator complex subunit MED25 binds transactivation domains (TADs) present in various cellular and viral proteins using two binding interfaces, named H1 and H2, which are found on opposite sides of its ACID domain. Here, we use and compare deep learning methods to characterize human MED25-TAD interfaces and assess the predicted models to published experimental data. For the H1 interface, AlphaFold produces predictions with high-reliability scores that agree well with experimental data, while the H2 interface predictions appear inconsistent, preventing reliable binding modes. Despite these limitations, we experimentally assess the validity of MED25 interface predictions with the viral transcriptional activators Lana-1 and IE62. AlphaFold predictions also suggest the existence of a unique hydrophobic pocket for the Arabidopsis MED25 ACID domain.

The molecular basis for cellular function of intrinsically disordered protein regions

Intrinsically disordered protein regions exist in a collection of dynamic interconverting conformations that lack a stable 3D structure. These regions are structurally heterogeneous, ubiquitous and found across all kingdoms of life. Despite the absence of a defined 3D structure, disordered regions are essential for cellular processes ranging from transcriptional control and cell signalling to subcellular organization. Through their conformational malleability and adaptability, disordered regions extend the repertoire of macromolecular interactions and are readily tunable by their structural and chemical context, making them ideal responders to regulatory cues. Recent work has led to major advances in understanding the link between protein sequence and conformational behaviour in disordered regions, yet the link between sequence and molecular function is less well defined. Here we consider the biochemical and biophysical foundations that underlie how and why disordered regions can engage in productive cellular functions, provide examples of emerging concepts and discuss how protein disorder contributes to intracellular information processing and regulation of cellular function.

Commonly asked questions about transcriptional activation domains

Eukaryotic transcription factors activate gene expression with their DNA-binding domains and activation domains. DNA-binding domains bind the genome by recognizing structurally related DNA sequences; they are structured, conserved, and predictable from protein sequences. Activation domains recruit chromatin modifiers, coactivator complexes, or basal transcriptional machinery via structurally diverse protein-protein interactions. Activation domains and DNA-binding domains have been called independent, modular units, but there are many departures from modularity, including interactions between these regions and overlap in function. Compared to DNA-binding domains, activation domains are poorly understood because they are poorly conserved, intrinsically disordered, and difficult to predict from protein sequences. This review, organized around commonly asked questions, describes recent progress that the field has made in understanding the sequence features that control activation domains and predicting them from sequence.

Rational optimization of a transcription factor activation domain inhibitor

Transcription factors are among the most attractive therapeutic targets but are considered largely 'undruggable' in part due to the intrinsically disordered nature of their activation domains. Here we show that the aromatic character of the activation domain of the androgen receptor, a therapeutic target for castration-resistant prostate cancer, is key for its activity as transcription factor, allowing it to translocate to the nucleus and partition into transcriptional condensates upon activation by androgens. On the basis of our understanding of the interactions stabilizing such condensates and of the structure that the domain adopts upon condensation, we optimized the structure of a small-molecule inhibitor previously identified by phenotypic screening. The optimized compounds had more affinity for their target, inhibited androgen-receptor-dependent transcriptional programs, and had an antitumorigenic effect in models of castration-resistant prostate cancer in cells and in vivo. These results suggest that it is possible to rationally optimize, and potentially even to design, small molecules that target the activation domains of oncogenic transcription factors.

A Novel EYA1 Mutation Causing Alternative RNA Splicing in a Chinese Family With Branchio-Oto Syndrome: Implications for Molecular Diagnosis and Clinical Application

OBJECTIVES

Branchio-oto syndrome (BOS) primarily manifests as hearing loss, preauricular pits, and branchial defects. EYA1 is the most common pathogenic gene, and splicing mutations account for a substantial proportion of cases. However, few studies have addressed the structural changes in the protein caused by splicing mutations and potential pathogenic factors, and several studies have shown that middle-ear surgery has limited effectiveness in improving hearing in these patients. BOS has also been relatively infrequently reported in the Chinese population. This study explored the genetic etiology in the family of a proband with BOS and provided clinical treatment to improve the patient's hearing.

METHODS

We collected detailed clinical features and peripheral blood samples from the patients and unaffected individuals within the family. Pathogenic mutations were identified by whole-exome sequencing and cosegregation analysis and classified according to the American College of Medical Genetics and Genomics guidelines. Alternative splicing was verified through a minigene assay. The predicted three-dimensional protein structure and biochemical experiments were used to investigate the pathogenicity of the mutation. The proband underwent middle-ear surgery and was followed up at 1 month and 6 months postoperatively to monitor auditory improvement.

RESULTS

A novel heterozygous EYA1 splicing variant (c.1050+4 A>C) was identified and classified as pathogenic (PVS1(RNA), PM2, PP1). Skipping of exon 11 of the EYA1 pre-mRNA was confirmed using a minigene assay. This mutation may impair EYA1-SIX1 interactions, as shown by an immunoprecipitation assay. The EYA1-Mut protein exhibited cellular mislocalization and decreased protein expression in cytological experiments. Middle-ear surgery significantly improved hearing loss caused by bone-conduction abnormalities in the proband.

CONCLUSION

We reported a novel splicing variant of EYA1 in a Chinese family with BOS and revealed the potential molecular pathogenic mechanism. The significant hearing improvement observed in the proband after middle-ear surgery provides a reference for auditory rehabilitation in similar patients.

Clusters of acidic and hydrophobic residues can predict acidic transcriptional activation domains from protein sequence

Transcription factors activate gene expression in development, homeostasis, and stress with DNA binding domains and activation domains. Although there exist excellent computational models for predicting DNA binding domains from protein sequence, models for predicting activation domains from protein sequence have lagged, particularly in metazoans. We recently developed a simple and accurate predictor of acidic activation domains on human transcription factors. Here, we show how the accuracy of this human predictor arises from the clustering of aromatic, leucine, and acidic residues, which together are necessary for acidic activation domain function. When we combine our predictor with the predictions of convolutional neural network (CNN) models trained in yeast, the intersection is more accurate than individual models, emphasizing that each approach carries orthogonal information. We synthesize these findings into a new set of activation domain predictions on human transcription factors.

Pervasive, conserved secondary structure in highly charged protein regions

Understanding how protein sequences confer function remains a defining challenge in molecular biology. Two approaches have yielded enormous insight yet are often pursued separately: structure-based, where sequence-encoded structures mediate function, and disorder-based, where sequences dictate physicochemical and dynamical properties which determine function in the absence of stable structure. Here we study highly charged protein regions (>40% charged residues), which are routinely presumed to be disordered. Using recent advances in structure prediction and experimental structures, we show that roughly 40% of these regions form well-structured helices. Features often used to predict disorder-high charge density, low hydrophobicity, low sequence complexity, and evolutionarily varying length-are also compatible with solvated, variable-length helices. We show that a simple composition classifier predicts the existence of structure far better than well-established heuristics based on charge and hydropathy. We show that helical structure is more prevalent than previously appreciated in highly charged regions of diverse proteomes and characterize the conservation of highly charged regions. Our results underscore the importance of integrating, rather than choosing between, structure- and disorder-based approaches.

Genome-wide analysis of KIX gene family for organ size regulation in soybean (Glycine max L.)

The KIX domain, conserved among various nuclear and co-activator factors, acts as a binding site that interacts with other transcriptional activators and co-activators, playing a crucial role in gene expression regulation. In plants, the KIX domain is involved in plant hormone signaling, stress response regulation, cell cycle control, and differentiation, indicating its potential relevance to crop productivity. This study aims to identify and characterize KIX domains within the soybean (Glycine max L.) genome to predict their potential role in improving crop productivity. The conservation and evolutionary history of the KIX domains were explored in 59 plant species, confirming the presence of the KIX domains in diverse plants. Specifically, 13 KIX domains were identified within the soybean genome and classified into four main groups, namely GmKIX8/9, GmMED15, GmHAC, and GmRECQL, through sequence alignment, structural analysis, and phylogenetic tree construction. Association analysis was performed between KIX domain haplotypes and soybean seed-related agronomic traits using re-sequencing data from a core collection of 422 accessions. The results revealed correlations between SNP variations observed in GmKIX8-3 and GmMED15-4 and soybean seed phenotypic traits. Additionally, transcriptome analysis confirmed significant expression of the KIX domains during the early stages of soybean seed development. This study provides the first characterization of the structural, expression, genomic haplotype, and molecular features of the KIX domain in soybean, offering a foundation for functional analysis of the KIX domain in soybean and other plants.

Surprising connections between DNA binding and function for the near-complete set of yeast transcription factors

DNA sequence-specific transcription factors (TFs) modulate transcription and chromatin architecture, acting from regulatory sites in enhancers and promoters of eukaryotic genes. How TFs locate their DNA targets and how multiple TFs cooperate to regulate individual genes is still unclear. Most yeast TFs are thought to regulate transcription via binding to upstream activating sequences, situated within a few hundred base pairs upstream of the regulated gene. While this model has been validated for individual TFs and specific genes, it has not been tested in a systematic way with the large set of yeast TFs. Here, we have integrated information on the binding and expression targets for the near-complete set of yeast TFs. While we found many instances of functional TF binding sites in upstream regulatory regions, we found many more instances that do not fit this model. In many cases, rapid TF depletion affects gene expression where there is no detectable binding of that TF to the upstream region of the affected gene. In addition, for most TFs, only a small fraction of bound TFs regulates the nearby gene, showing that TF binding does not automatically correspond to regulation of the linked gene. Finally, we found that only a small percentage of TFs are exclusively strong activators or repressors with most TFs having dual function. Overall, our comprehensive mapping of TF binding and regulatory targets have both confirmed known TF relationships and revealed surprising properties of TF function.

The trans-regulatory landscape of gene networks in plants

The transcriptional effector domains of transcription factors play a key role in controlling gene expression; however, their functional nature is poorly understood, hampering our ability to explore this fundamental dimension of gene regulatory networks. To map the trans-regulatory landscape in a complex eukaryote, we systematically characterized the putative transcriptional effector domains of over 400 Arabidopsis thaliana transcription factors for their capacity to modulate transcription. We demonstrate that transcriptional effector activity can be integrated into gene regulatory networks capable of elucidating the functional dynamics underlying gene expression patterns. We further show how our characterized domains can enhance genome engineering efforts and reveal how plant transcriptional activators share regulatory features conserved across distantly related eukaryotes. Our results provide a framework to systematically characterize the regulatory role of transcription factors at a genome-scale in order to understand the transcriptional wiring of biological systems.

Towards sequence-based principles for protein phase separation predictions

The phenomenon of protein phase separation, which underlies the formation of biomolecular condensates, has been associated with numerous cellular functions. Recent studies indicate that the amino acid sequences of most proteins may harbour not only the code for folding into the native state but also for condensing into the liquid-like droplet state and the solid-like amyloid state. Here we review the current understanding of the principles for sequence-based methods for predicting the propensity of proteins for phase separation. A guiding concept is that entropic contributions are generally more important to stabilise the droplet state than they are for the native and amyloid states. Although estimating these entropic contributions has proven difficult, we describe some progress that has been recently made in this direction. To conclude, we discuss the challenges ahead to extend sequence-based prediction methods of protein phase separation to include quantitative in vivo characterisations of this process.

Alternatively spliced exon regulates context-dependent MEF2D higher-order assembly during myogenesis

During muscle cell differentiation, the alternatively spliced, acidic β-domain potentiates transcription of Myocyte-specific Enhancer Factor 2 (Mef2D). Sequence analysis by the FuzDrop method indicates that the β-domain can serve as an interaction element for Mef2D higher-order assembly. In accord, we observed Mef2D mobile nuclear condensates in C2C12 cells, similar to those formed through liquid-liquid phase separation. In addition, we found Mef2D solid-like aggregates in the cytosol, the presence of which correlated with higher transcriptional activity. In parallel, we observed a progress in the early phase of myotube development, and higher MyoD and desmin expression. In accord with our predictions, the formation of aggregates was promoted by rigid β-domain variants, as well as by a disordered β-domain variant, capable of switching between liquid-like and solid-like higher-order states. Along these lines, NMR and molecular dynamics simulations corroborated that the β-domain can sample both ordered and disordered interactions leading to compact and extended conformations. These results suggest that β-domain fine-tunes Mef2D higher-order assembly to the cellular context, which provides a platform for myogenic regulatory factors and the transcriptional apparatus during the developmental process.

Structural basis of a transcription pre-initiation complex on a divergent promoter

Most eukaryotic promoter regions are divergently transcribed. As the RNA polymerase II pre-initiation complex (PIC) is intrinsically asymmetric and responsible for transcription in a single direction, it is unknown how divergent transcription arises. Here, the Saccharomyces cerevisiae Mediator complexed with a PIC (Med-PIC) was assembled on a divergent promoter and analyzed by cryoelectron microscopy. The structure reveals two distinct Med-PICs forming a dimer through the Mediator tail module, induced by a homodimeric activator protein localized near the dimerization interface. The tail dimer is associated with ∼80-bp upstream DNA, such that two flanking core promoter regions are positioned and oriented in a suitable form for PIC assembly in opposite directions. Also, cryoelectron tomography visualized the progress of the PIC assembly on the two core promoter regions, providing direct evidence for the role of the Med-PIC dimer in divergent transcription.

Integrating deep mutational scanning and low-throughput mutagenesis data to predict the impact of amino acid variants

BACKGROUND

Evaluating the impact of amino acid variants has been a critical challenge for studying protein function and interpreting genomic data. High-throughput experimental methods like deep mutational scanning (DMS) can measure the effect of large numbers of variants in a target protein, but because DMS studies have not been performed on all proteins, researchers also model DMS data computationally to estimate variant impacts by predictors.

RESULTS

In this study, we extended a linear regression-based predictor to explore whether incorporating data from alanine scanning (AS), a widely used low-throughput mutagenesis method, would improve prediction results. To evaluate our model, we collected 146 AS datasets, mapping to 54 DMS datasets across 22 distinct proteins.

CONCLUSIONS

We show that improved model performance depends on the compatibility of the DMS and AS assays, and the scale of improvement is closely related to the correlation between DMS and AS results.

Protein interactions: anything new?

How do proteins interact in the cellular environment? Which interactions stabilize liquid-liquid phase separated condensates? Are the concepts, which have been developed for specific protein complexes also applicable to higher-order assemblies? Recent discoveries prompt for a universal framework for protein interactions, which can be applied across the scales of protein communities. Here, we discuss how our views on protein interactions have evolved from rigid structures to conformational ensembles of proteins and discuss the open problems, in particular related to biomolecular condensates. Protein interactions have evolved to follow changes in the cellular environment, which manifests in multiple modes of interactions between the same partners. Such cellular context-dependence requires multiplicity of binding modes (MBM) by sampling multiple minima of the interaction energy landscape. We demonstrate that the energy landscape framework of protein folding can be applied to explain this phenomenon, opening a perspective toward a physics-based, universal model for cellular protein behaviors.

Involvement of the SAGA and TFIID coactivator complexes in transcriptional dysregulation caused by the separation of core and tail Mediator modules

Regulation of RNA polymerase II transcription requires the concerted efforts of several multisubunit coactivator complexes, which interact with the RNA polymerase II preinitiation complex to stimulate transcription. We previously showed that separation of the Mediator core from Mediator's tail module results in modest overactivation of genes annotated as highly dependent on TFIID for expression. However, it is unclear if other coactivators are involved in this phenomenon. Here, we show that the overactivation of certain genes by Mediator core/tail separation is blunted by disruption of the Spt-Ada-Gcn5-Acetyl transferase complex through the removal of its structural Spt20 subunit, though this downregulation does not appear to completely depend on reduced Spt-Ada-Gcn5-Acetyl transferase association with the genome. Consistent with the enrichment of TFIID-dependent genes among genes overactivated by Mediator core/tail separation, depletion of the essential TFIID subunit Taf13 suppressed the overactivation of these genes when Med16 was simultaneously removed. As with Spt-Ada-Gcn5-Acetyl transferase, this effect did not appear to be fully dependent on the reduced genomic association of TFIID. Given that the observed changes in gene expression could not be clearly linked to alterations in Spt-Ada-Gcn5-Acetyl transferase or TFIID occupancy, our data may suggest that the Mediator core/tail connection is important for the modulation of Spt-Ada-Gcn5-Acetyl transferase and/or TFIID conformation and/or function at target genes.

Yeast Mediator facilitates transcription initiation at most promoters via a Tail-independent mechanism

Mediator (MED) is a conserved factor with important roles in basal and activated transcription. Here, we investigate the genome-wide roles of yeast MED by rapid depletion of its activator-binding domain (Tail) and monitoring changes in nascent transcription. Rapid Tail depletion surprisingly reduces transcription from only a small subset of genes. At most of these Tail-dependent genes, in unperturbed conditions, MED is detected at both the UASs and promoters. In contrast, at most Tail-independent genes, we find MED primarily at promoters but not at the UASs. These results suggest that MED Tail and activator-mediated MED recruitment regulates only a small subset of genes. Furthermore, we define three classes of genes that differ in PIC assembly pathways and the requirements for MED Tail, SAGA, TFIID, and BET factors Bdf1/2. Our combined results have broad implications for the roles of MED, other coactivators, and mechanisms of transcriptional regulation at different gene classes.

Transcription factors perform a 2-step search of the nucleus

Transcription factors regulate gene expression by binding to regulatory DNA and recruiting regulatory protein complexes. The DNA-binding and protein-binding functions of transcription factors are traditionally described as independent functions performed by modular protein domains. Here, I argue that genome binding can be a 2-part process with both DNA-binding and protein-binding steps, enabling transcription factors to perform a 2-step search of the nucleus to find their appropriate binding sites in a eukaryotic genome. I support this hypothesis with new and old results in the literature, discuss how this hypothesis parsimoniously resolves outstanding problems, and present testable predictions.

The Origins and Functions of De Novo Genes: Against All Odds?

"De novo" genes evolve from previously non-genic DNA. This strikes many of us as remarkable, because it seems extraordinarily unlikely that random sequence would produce a functional gene. How is this possible? In this two-part review, I first summarize what is known about the origins and molecular functions of the small number of de novo genes for which such information is available. I then speculate on what these examples may tell us about how de novo genes manage to emerge despite what seem like enormous opposing odds.

RING 3.0: fast generation of probabilistic residue interaction networks from structural ensembles

Residue interaction networks (RINs) are used to represent residue contacts in protein structures. Thanks to the advances in network theory, RINs have been proved effective as an alternative to coordinate data in the analysis of complex systems. The RING server calculates high quality and reliable non-covalent molecular interactions based on geometrical parameters. Here, we present the new RING 3.0 version extending the previous functionality in several ways. The underlying software library has been re-engineered to improve speed by an order of magnitude. RING now also supports the mmCIF format and provides typed interactions for the entire PDB chemical component dictionary, including nucleic acids. Moreover, RING now employs probabilistic graphs, where multiple conformations (e.g. NMR or molecular dynamics ensembles) are mapped as weighted edges, opening up new ways to analyze structural data. The web interface has been expanded to include a simultaneous view of the RIN alongside a structure viewer, with both synchronized and clickable. Contact evolution across models (or time) is displayed as a heatmap and can help in the discovery of correlating interaction patterns. The web server, together with an extensive help and tutorial, is available from URL: https://ring.biocomputingup.it/.

Directed mutational scanning reveals a balance between acidic and hydrophobic residues in strong human activation domains

Acidic activation domains are intrinsically disordered regions of the transcription factors that bind coactivators. The intrinsic disorder and low evolutionary conservation of activation domains have made it difficult to identify the sequence features that control activity. To address this problem, we designed thousands of variants in seven acidic activation domains and measured their activities with a high-throughput assay in human cell culture. We found that strong activation domain activity requires a balance between the number of acidic residues and aromatic and leucine residues. These findings motivated a predictor of acidic activation domains that scans the human proteome for clusters of aromatic and leucine residues embedded in regions of high acidity. This predictor identifies known activation domains and accurately predicts previously unidentified ones. Our results support a flexible acidic exposure model of activation domains in which the acidic residues solubilize hydrophobic motifs so that they can interact with coactivators. A record of this paper's transparent peer review process is included in the supplemental information.

Mediator Engineering of Saccharomyces cerevisiae To Improve Multidimensional Stress Tolerance

Saccharomyces cerevisiae is a well-performing workhorse in chemical production, which encounters complex environmental stresses during industrial processes. We constructed a multiple stress tolerance mutant, Med15^V76R/R84K, that was obtained by engineering the KIX domain of Mediator tail subunit Med15. Med15^V76R/R84K interacted with transcription factor Hap5 to improve ARV1 expression for sterol homeostasis for decreasing membrane fluidity and thereby enhancing acid tolerance. Med15^V76R/R84K interacted with transcription factor Mga2 to improve GIT1 expression for phospholipid biosynthesis for increasing membrane integrity and thereby improving oxidative tolerance. Med15^V76R/R84K interacted with transcription factor Aft1 to improve NFT1 expression for inorganic ion transport for reducing membrane permeability and thereby enhancing osmotic tolerance. Based on this Med15 mutation, Med15^V76R/R84K, the engineered S. cerevisiae strain, showed a 28.1% increase in pyruvate production in a 1.0-L bioreactor compared to that of S. cerevisiae with its native Med15. These results indicated that Mediator engineering provides a potential alternative for improving multidimensional stress tolerance in S. cerevisiae. IMPORTANCE This study identified the role of the KIX domain of Mediator tail subunit Med15 in response to acetic acid, H₂O₂, and NaCl in S. cerevisiae. Engineered KIX domain by protein engineering, the mutant strain Med15^V76R/R84K, increased multidimensional stress tolerance and pyruvate production compared with that of S. cerevisiae with its native Med15. The Med15^V76R/R84K could increase membrane related genes expression possibly by enhancing interaction with transcription factor to improve membrane physiological functions under stress conditions.

Drugging Fuzzy Complexes in Transcription

Transcription factors (TFs) are one of the most promising but underutilized classes of drug targets. The high degree of intrinsic disorder in both the structure and the interactions (i.e., “fuzziness”) of TFs is one of the most important challenges to be addressed in this context. Here, we discuss the impacts of fuzziness on transcription factor drug discovery, describing how disorder poses fundamental problems to the typical drug design, and screening approaches used for other classes of proteins such as receptors or enzymes. We then speculate on ways modern biophysical and chemical biology approaches could synergize to overcome many of these challenges by directly addressing the challenges imposed by TF disorder and fuzziness.

Molecular Determinants of Selectivity in Disordered Complexes May Shed Light on Specificity in Protein Condensates

Biomolecular condensates challenge the classical concepts of molecular recognition. The variable composition and heterogeneous conformations of liquid-like protein droplets are bottlenecks for high-resolution structural studies. To obtain atomistic insights into the organization of these assemblies, here we have characterized the conformational ensembles of specific disordered complexes, including those of droplet-driving proteins. First, we found that these specific complexes exhibit a high degree of conformational heterogeneity. Second, we found that residues forming contacts at the interface also sample many conformations. Third, we found that different patterns of contacting residues form the specific interface. In addition, we observed a wide range of sequence motifs mediating disordered interactions, including charged, hydrophobic and polar contacts. These results demonstrate that selective recognition can be realized by variable patterns of weakly defined interaction motifs in many different binding configurations. We propose that these principles also play roles in determining the selectivity of biomolecular condensates.

In Silico Modeling of the Influence of Environment on Amyloid Folding Using FOD-M Model

The role of the environment in amyloid formation based on the fuzzy oil drop model (FOD) is discussed here. This model assumes that the hydrophobicity distribution within a globular protein is consistent with a 3D Gaussian (3DG) distribution. Such a distribution is interpreted as the idealized effect of the presence of a polar solvent-water. A chain with a sequence of amino acids (which are bipolar molecules) determined by evolution recreates a micelle-like structure with varying accuracy. The membrane, which is a specific environment with opposite characteristics to the polar aquatic environment, directs the hydrophobic residues towards the surface. The modification of the FOD model to the FOD-M form takes into account the specificity of the cell membrane. It consists in "inverting" the 3DG distribution (complementing the Gaussian distribution), which expresses the exposure of hydrophobic residues on the surface. It turns out that the influence of the environment for any protein (soluble or membrane-anchored) is the result of a consensus factor expressing the participation of the polar environment and the "inverted" environment. The ratio between the proportion of the aqueous and the "reversed" environment turns out to be a characteristic property of a given protein, including amyloid protein in particular. The structure of amyloid proteins has been characterized in the context of prion, intrinsically disordered, and other non-complexing proteins to cover a wider spectrum of molecules with the given characteristics based on the FOD-M model.

Direct photoresponsive inhibition of a p53-like transcription activation domain in PIF3 by Arabidopsis phytochrome B

Photoactivated phytochrome B (PHYB) binds to antagonistically acting PHYTOCHROME-INTERACTING transcription FACTORs (PIFs) to regulate hundreds of light responsive genes in Arabidopsis by promoting PIF degradation. However, whether PHYB directly controls the transactivation activity of PIFs remains ambiguous. Here we show that the prototypic PIF, PIF3, possesses a p53-like transcription activation domain (AD) consisting of a hydrophobic activator motif flanked by acidic residues. A PIF3mAD mutant, in which the activator motif is replaced with alanines, fails to activate PIF3 target genes in Arabidopsis, validating the functions of the PIF3 AD in vivo. Intriguingly, the N-terminal photosensory module of PHYB binds immediately adjacent to the PIF3 AD to repress PIF3’s transactivation activity, demonstrating a novel PHYB signaling mechanism through direct interference of the transactivation activity of PIF3. Our findings indicate that PHYB, likely also PHYA, controls the stability and activity of PIFs via structurally separable dual signaling mechanisms.

Photoactivated phytochrome B regulates gene expression by interacting with PIF transcription factors. Here the authors show that PIF3 contains a p53-like transcription activation domain (AD) and that PHYB can directly suppress PIF3 transactivation activity by binding adjacent to the AD.

Quantification of Conformational Entropy Unravels Effect of Disordered Flanking Region in Coupled Folding and Binding

Intrinsic disorder (ID) constitutes a new dimension to the protein structure-function relationship. The ability to undergo conformational changes upon binding is a key property of intrinsically disordered proteins and remains challenging to study using conventional methods. A 1994 paper by R. S. Spolar and M. T. Record presented a thermodynamic approach for estimating changes in conformational entropy based on heat capacity changes, allowing quantification of residues folding upon binding. Here, we adapt the method for studies of intrinsically disordered proteins. We integrate additional data to provide a broader experimental foundation for the underlying relations and, based on >500 protein-protein complexes involving disordered proteins, reassess a key relation between polar and nonpolar surface area changes, previously determined using globular protein folding. We demonstrate the improved suitability of the adapted method to studies of the folded αα-hub domain RST from radical-induced cell death 1, whose interactome is characterized by ID. From extensive thermodynamic data, quantifying the conformational entropy changes upon binding, and comparison to the NMR structure, the adapted method improves accuracy for ID-based studies. Furthermore, we apply the method, in conjunction with NMR, to reveal hitherto undetected effects of interaction-motif context. Thus, inclusion of the disordered context of the DREB2A RST-binding motif induces structuring of the binding motif, resulting in major enthalpy-entropy compensation in the interaction interface. This study, also evaluating additional interactions, demonstrates the strength of the ID-adapted Spolar-Record thermodynamic approach for dissection of structural features of ID-based interactions, easily overlooked in traditional studies, and for translation of these into mechanistic knowledge.

The sequence-ensemble relationship in fuzzy protein complexes

Intrinsically disordered proteins (IDPs) interact with globular proteins through a variety of mechanisms, resulting in the structurally heterogeneous ensembles known as fuzzy complexes. While there exists a reasonable comprehension on how IDP sequence determines the unbound IDP ensemble, little is known about what shapes the structural characteristics of IDPs bound to their targets. Using a statistical thermodynamic model, we show that the target-bound ensembles are determined by a simple code that combines the IDP sequence and the distribution of IDP-target interaction hotspots. These two parameters define the conformational space of target-bound IDPs and rationalize the observed structural heterogeneity of fuzzy complexes. The presented model successfully reproduces the dynamical signatures of target-bound IDPs from the NMR relaxation experiments as well as the changes of interaction affinity and the IDP helicity induced by mutations. The model explains how the target-bound IDP ensemble adapts to mutations in order to achieve an optimal balance between conformational freedom and interaction energy. Taken together, the presented sequence-ensemble relationship of fuzzy complexes explains the different manifestations of IDP disorder in folding-upon-binding processes.

Thermodynamic Models for Assembly of Intrinsically Disordered Protein Hubs with Multiple Interaction Partners

Prevalent in diverse protein interactomes, intrinsically disordered proteins or regions (IDPs or IDRs) often drive assembly of higher-order macromolecular complexes, using multiple target-binding motifs. Such IDP hubs are suggested to process various cellular signals and coordinate relevant biological processes. However, the mechanism of assembly and functional regulation of IDP hubs remains elusive due to the challenges in dissecting their intricate protein-protein interaction networks. Here we present basic thermodynamic models for the assembly of simple IDP hubs with multiple target proteins, constructing partition functions from fundamental binding parameters. We combined these basic functions to develop advanced thermodynamic models to analyze the assembly of the Nup153 hubs interacting with multiple karyopherin β1 (Kap) molecules, critical components of nucleocytoplasmic transport. The thermodynamic analysis revealed a complex organization of the Kap binding sites within the C-terminal IDR of Nup153 including a high-affinity 1:1 interaction site and a series of low-affinity sites for binding of multiple Kaps with negative cooperativity. The negative cooperativity arises from the overlapping nature of the low-affinity sites where Kap occupies multiple dipeptide motifs. The quantitative dissection culminated in construction of the Nup153 hub ensemble, elucidating how distribution among various Kap-bound states is modulated by Kap concentration and competing nuclear proteins. In particular, the Kap occupancy of the IDR can be fine-tuned by varying the location of competition within the overlapping sites, suggesting coupling of specific nuclear processes to distinct transport activities. In general, our results demonstrate the feasibility and a potential mechanism for manifold regulation of IDP functions by diverse cellular signals.

The SAGA complex regulates early steps in transcription via its deubiquitylase module subunit USP22

The SAGA coactivator complex is essential for eukaryotic transcription and comprises four distinct modules, one of which contains the ubiquitin hydrolase USP22. In yeast, the USP22 ortholog deubiquitylates H2B, resulting in Pol II Ser2 phosphorylation and subsequent transcriptional elongation. In contrast to this H2B-associated role in transcription, we report here that human USP22 contributes to the early stages of stimulus-responsive transcription, where USP22 is required for pre-initiation complex (PIC) stability. Specifically, USP22 maintains long-range enhancer-promoter contacts and controls loading of Mediator tail and general transcription factors (GTFs) onto promoters, with Mediator core recruitment being USP22-independent. In addition, we identify Mediator tail subunits MED16 and MED24 and the Pol II subunit RBP1 as potential non-histone substrates of USP22. Overall, these findings define a role for human SAGA within the earliest steps of transcription.

Connection of core and tail Mediator modules restrains transcription from TFIID-dependent promoters

The Mediator coactivator complex is divided into four modules: head, middle, tail, and kinase. Deletion of the architectural subunit Med16 separates core Mediator (cMed), comprising the head, middle, and scaffold (Med14), from the tail. However, the direct global effects of tail/cMed disconnection are unclear. We find that rapid depletion of Med16 downregulates genes that require the SAGA complex for full expression, consistent with their reported tail dependence, but also moderately overactivates TFIID-dependent genes in a manner partly dependent on the separated tail, which remains associated with upstream activating sequences. Suppression of TBP dynamics via removal of the Mot1 ATPase partially restores normal transcriptional activity to Med16-depleted cells, suggesting that cMed/tail separation results in an imbalance in the levels of PIC formation at SAGA-requiring and TFIID-dependent genes. We propose that the preferential regulation of SAGA-requiring genes by tailed Mediator helps maintain a proper balance of transcription between these genes and those more dependent on TFIID.

Composed of over two dozen subunits, the Mediator complex plays several roles in RNA polymerase II (RNAPII) transcription in eukaryotes. In yeast, deletion of Med16, which splits Mediator into two stable subcomplexes, both increases and decreases transcript levels, suggesting that Med16 might play a repressive role. However, the direct effects of Med16 removal on RNAPII transcription have not been assessed, owing to the use of deletion mutants and measurement of steady-state RNA levels in prior studies. Here, using a combination of inducible protein depletion and analysis of nascent RNA, we find that Med16 removal 1) downregulates a small group of genes reported to be highly dependent on the SAGA complex and 2) upregulates a larger set of genes reported to be more dependent on the TFIID complex in a manner dependent on another component of Mediator. We find that artificially altering the balance of transcription pre-initiation complex (PIC) formation toward SAGA-requiring promoters and away from TFIID-dependent promoters partially restores normal transcription, indicating a contribution of altered PIC formation to the transcriptional alterations observed with Med16 loss. Taken together, our results indicate that the structural integrity of Mediator is important for maintaining balanced transcription between different gene classes.