Metapredict: a fast, accurate, and easy-to-use predictor of consensus disorder and structure

Exploratory comparative transcriptomic analysis reveals potential gene targets associated with Cry1A.105 and Cry2Ab2 resistance in fall armyworm (Spodoptera frugiperda)

Genetically modified (GM) crops, expressing Bacillus thuringiensis (Bt) insecticidal toxins, have substantially transformed agriculture. Despite rapid adoption, their environmental and economic benefits face scrutiny due to unsustainable agricultural practices and the emergence of resistant pests like Spodoptera frugiperda, known as the fall armyworm (FAW). FAW's adaptation to Bt technology in corn and cotton compromises the long-term efficacy of Bt crops. To advance the understanding of the genetic foundations of resistance mechanisms, we conducted an exploratory comparative transcriptomic analysis of two divergent FAW populations. One population exhibited practical resistance to the Bt insecticidal proteins Cry1A.105 and Cry2Ab2, expressed in the genetically engineered MON-89Ø34 - 3 maize, while the other population remained susceptible to these proteins. Differential expression analysis supported that Cry1A.105 and Cry2Ab2 significantly affect the FAW physiology. A total of 247 and 254 differentially expressed genes were identified in the Cry-resistant and susceptible populations, respectively. By integrating our findings with established literature and databases, we underscored 53 gene targets potentially involved in FAW's resistance to Cry1A.105 and Cry2Ab2. In particular, we considered and discussed the potential roles of the differentially expressed genes encoding ABC transporters, G protein-coupled receptors, the P450 enzymatic system, and other Bt-related detoxification genes. Based on these findings, we emphasize the importance of exploratory transcriptomic analyses to uncover potential gene targets involved with Bt insecticidal proteins resistance, and to support the advantages of GM crops in the face of emerging challenges.

Systematic identification of transcriptional activation domains from non-transcription factor proteins in plants and yeast

Transcription factors can promote gene expression through activation domains. Whole-genome screens have systematically mapped activation domains in transcription factors but not in non-transcription factor proteins (e.g., chromatin regulators and coactivators). To fill this knowledge gap, we employed the activation domain predictor PADDLE to analyze the proteomes of Arabidopsis thaliana and Saccharomyces cerevisiae. We screened 18,000 predicted activation domains from >800 non-transcription factor genes in both species, confirming that 89% of candidate proteins contain active fragments. Our work enables the annotation of hundreds of nuclear proteins as putative coactivators, many of which have never been ascribed any function in plants. Analysis of peptide sequence compositions reveals how the distribution of key amino acids dictates activity. Finally, we validated short, "universal" activation domains with comparable performance to state-of-the-art activation domains used for genome engineering. Our approach enables the genome-wide discovery and annotation of activation domains that can function across diverse eukaryotes.

Identification of plant transcriptional activation domains

Gene expression in Arabidopsis is regulated by more than 1,900 transcription factors (TFs), which have been identified genome-wide by the presence of well-conserved DNA-binding domains. Activator TFs contain activation domains (ADs) that recruit coactivator complexes; however, for nearly all Arabidopsis TFs, we lack knowledge about the presence, location and transcriptional strength of their ADs1. To address this gap, here we use a yeast library approach to experimentally identify Arabidopsis ADs on a proteome-wide scale, and find that more than half of the Arabidopsis TFs contain an AD. We annotate 1,553 ADs, the vast majority of which are, to our knowledge, previously unknown. Using the dataset generated, we develop a neural network to accurately predict ADs and to identify sequence features that are necessary to recruit coactivator complexes. We uncover six distinct combinations of sequence features that result in activation activity, providing a framework to interrogate the subfunctionalization of ADs. Furthermore, we identify ADs in the ancient AUXIN RESPONSE FACTOR family of TFs, revealing that AD positioning is conserved in distinct clades. Our findings provide a deep resource for understanding transcriptional activation, a framework for examining function in intrinsically disordered regions and a predictive model of ADs.

Rational Tuning of the Concentration-independent Enrichment of Prion-like Domains in Stress Granules

Stress granules (SGs) are large ribonucleoprotein assemblies that form in response to acute stress in eukaryotes. SG formation is thought to be initiated by liquid-liquid phase separation (LLPS) of key proteins and RNA. These molecules serve as a scaffold for recruitment of client molecules. LLPS of scaffold proteins in vitro is highly concentration-dependent, yet biomolecular condensates in vivo contain hundreds of unique proteins, most of which are thought to be clients rather than scaffolds. Many proteins that localize to SGs contain low-complexity, prion-like domains (PrLDs) that have been implicated in LLPS and SG recruitment. The degree of enrichment of proteins in biomolecular condensates such as SGs can vary widely, but the underlying basis for these differences is not fully understood. Here, we develop a toolkit of model PrLDs to examine the factors that govern efficiency of PrLD recruitment to stress granules. Recruitment was highly sensitive to amino acid composition: enrichment in SGs could be tuned through subtle changes in hydrophobicity. By contrast, SG recruitment was largely insensitive to PrLD concentration at both a population level and single-cell level. These observations point to a model wherein PrLDs are enriched in SGs through either simple solvation effects or interactions that are effectively non-saturable even at high expression levels.

FLiPPR: A Processor for Limited Proteolysis (LiP) Mass Spectrometry Data Sets Built on FragPipe

Here, we present FLiPPR, or FragPipe LiP (limited proteolysis) Processor, a tool that facilitates the analysis of data from limited proteolysis mass spectrometry (LiP-MS) experiments following primary search and quantification in FragPipe. LiP-MS has emerged as a method that can provide proteome-wide information on protein structure and has been applied to a range of biological and biophysical questions. Although LiP-MS can be carried out with standard laboratory reagents and mass spectrometers, analyzing the data can be slow and poses unique challenges compared to typical quantitative proteomics workflows. To address this, we leverage FragPipe and then process its output in FLiPPR. FLiPPR formalizes a specific data imputation heuristic that carefully uses missing data in LiP-MS experiments to report on the most significant structural changes. Moreover, FLiPPR introduces a data merging scheme and a protein-centric multiple hypothesis correction scheme, enabling processed LiP-MS data sets to be more robust and less redundant. These improvements strengthen statistical trends when previously published data are reanalyzed with the FragPipe/FLiPPR workflow. We hope that FLiPPR will lower the barrier for more users to adopt LiP-MS, standardize statistical procedures for LiP-MS data analysis, and systematize output to facilitate eventual larger-scale integration of LiP-MS data.

An activity-specificity trade-off encoded in human transcription factors

Transcription factors (TFs) control specificity and activity of gene transcription, but whether a relationship between these two features exists is unclear. Here we provide evidence for an evolutionary trade-off between the activity and specificity in human TFs encoded as submaximal dispersion of aromatic residues in their intrinsically disordered protein regions. We identified approximately 500 human TFs that encode short periodic blocks of aromatic residues in their intrinsically disordered regions, resembling imperfect prion-like sequences. Mutation of periodic aromatic residues reduced transcriptional activity, whereas increasing the aromatic dispersion of multiple human TFs enhanced transcriptional activity and reprogramming efficiency, promoted liquid-liquid phase separation in vitro and more promiscuous DNA binding in cells. Together with recent work on enhancer elements, these results suggest an important evolutionary role of suboptimal features in transcriptional control. We propose that rational engineering of amino acid features that alter phase separation may be a strategy to optimize TF-dependent processes, including cellular reprogramming.

Posttranslational regulation of photosynthetic activity via the TOR kinase in plants

Chloroplasts are the powerhouse of the plant cell, and their activity must be matched to plant growth to avoid photooxidative damage. We have identified a posttranslational mechanism linking the eukaryotic target of rapamycin (TOR) kinase that promotes growth and the guanosine tetraphosphate (ppGpp) signaling pathway of prokaryotic origins that regulates chloroplast activity and photosynthesis in particular. We find that RelA SpoT homolog 3 (RSH3), a nuclear-encoded enzyme responsible for ppGpp biosynthesis, interacts directly with the TOR complex via a plant-specific amino-terminal region which is phosphorylated in a TOR-dependent manner. Down-regulating TOR activity causes a rapid increase in ppGpp synthesis in RSH3 overexpressors and reduces photosynthetic capacity in an RSH-dependent manner in wild-type plants. The TOR-RSH3 signaling axis therefore regulates the equilibrium between chloroplast activity and plant growth, setting a precedent for the regulation of organellar function by TOR.

MATR3's Role beyond the Nuclear Matrix: From Gene Regulation to Its Implications in Amyotrophic Lateral Sclerosis and Other Diseases

Matrin-3 (MATR3) was initially discovered as a component of the nuclear matrix about thirty years ago. Since then, accumulating studies have provided evidence that MATR3 not only plays a structural role in the nucleus, but that it is also an active protein involved in regulating gene expression at multiple levels, including chromatin organization, DNA transcription, RNA metabolism, and protein translation in the nucleus and cytoplasm. Furthermore, MATR3 may play a critical role in various cellular processes, including DNA damage response, cell proliferation, differentiation, and survival. In addition to the revelation of its biological role, recent studies have reported MATR3's involvement in the context of various diseases, including neurodegenerative and neurodevelopmental diseases, as well as cancer. Moreover, sequencing studies of patients revealed a handful of disease-associated mutations in MATR3 linked to amyotrophic lateral sclerosis (ALS), which further elevated the gene's importance as a topic of study. In this review, we synthesize the current knowledge regarding the diverse functions of MATR3 in DNA- and RNA-related processes, as well as its involvement in various diseases, with a particular emphasis on ALS.

The Rapid Evolution of De Novo Proteins in Structure and Complex

Recent studies in the rice genome-wide have established that de novo genes, evolving from noncoding sequences, enhance protein diversity through a stepwise process. However, the pattern and rate of their evolution in protein structure over time remain unclear. Here, we addressed these issues within a surprisingly short evolutionary timescale (<1 million years for 97% of Oryza de novo genes) with comparative approaches to gene duplicates. We found that de novo genes evolve faster than gene duplicates in the intrinsically disordered regions (such as random coils), secondary structure elements (such as α helix and β strand), hydrophobicity, and molecular recognition features. In de novo proteins, specifically, we observed an 8% to 14% decay in random coils and intrinsically disordered region lengths and a 2.3% to 6.5% increase in structured elements, hydrophobicity, and molecular recognition features, per million years on average. These patterns of structural evolution align with changes in amino acid composition over time as well. We also revealed higher positive charges but smaller molecular weights for de novo proteins than duplicates. Tertiary structure predictions showed that most de novo proteins, though not typically well folded on their own, readily form low-energy and compact complexes with other proteins facilitated by extensive residue contacts and conformational flexibility, suggesting a faster-binding scenario in de novo proteins to promote interaction. These analyses illuminate a rapid evolution of protein structure in de novo genes in rice genomes, originating from noncoding sequences, highlighting their quick transformation into active, protein complex-forming components within a remarkably short evolutionary timeframe.

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.

Conserved molecular recognition by an intrinsically disordered region in the absence of sequence conservation

Intrinsically disordered regions (IDRs) are critical for cellular function yet often appear to lack sequence conservation when assessed by multiple sequence alignments. This raises the question of if and how function can be encoded and preserved in these regions despite massive sequence variation. To address this question, we have applied coarse-grained molecular dynamics simulations to investigate non-specific RNA binding of coronavirus nucleocapsid proteins. Coronavirus nucleocapsid proteins consist of multiple interspersed disordered and folded domains that bind RNA. Here, we focus on the first two domains of coronavirus nucleocapsid proteins: the disordered N-terminal domain (NTD) and the folded RNA binding domain (RBD). While the NTD is highly variable across evolution, the RBD is structurally conserved. This combination makes the NTD-RBD a convenient model system for exploring the interplay between an IDR adjacent to a folded domain and how changes in IDR sequence can influence molecular recognition of a partner. Our results reveal a surprising degree of sequence-specificity encoded by both the composition and the precise order of the amino acids in the NTD. The presence of an NTD can - depending on the sequence - either suppress or enhance RNA binding. Despite this sensitivity, large-scale variation in NTD sequences is possible while certain sequence features are retained. Consequently, a conformationally-conserved dynamic and disordered RNA:protein complex is found across nucleocapsid protein orthologs despite large-scale changes in both NTD sequence and RBD surface chemistry. Taken together, these insights shed light on the ability of disordered regions to preserve functional characteristics despite their sequence variability.

Overview of AlphaFold2 and breakthroughs in overcoming its limitations

Predicting three-dimensional (3D) protein structures has been challenging for decades. The emergence of AlphaFold2 (AF2), a deep learning-based machine learning method developed by DeepMind, became a game changer in the protein folding community. AF2 can predict a protein's three-dimensional structure with high confidence based on its amino acid sequence. Accurate prediction of protein structures can dramatically accelerate our understanding of biological mechanisms and provide a solid foundation for reliable drug design. Although AF2 breaks through the barriers in predicting protein structures, many rooms remain to be further studied. This review provides a brief historical overview of the development of protein structure prediction, covering template-based, template-free, and machine learning-based methods. In addition to reviewing the potential benefits (Pros) and considerations (Cons) of using AF2, this review summarizes the diverse applications, including protein structure predictions, dynamic changes, point mutation, integration of language model and experimental data, protein complex, and protein-peptide interaction. It underscores recent advancements in efficiency, reliability, and broad application of AF2. This comprehensive review offers valuable insights into the applications of AF2 and AF2-inspired AI methods in structural biology and its potential for clinically significant drug target discovery.

Structural dynamics of the intrinsically disordered linker region of cardiac troponin T

The cardiac troponin complex, composed of troponins I, T, and C, plays a central role in regulating the calcium-dependent interactions between myosin and the thin filament. Mutations in troponin can cause cardiomyopathies; however, it is still a major challenge for the field to connect how changes in sequence affect troponin's function. Recent high-resolution structures of the thin filament revealed critical insights into the structure-function relationship of the troponin complex, but there remain large, unresolved segments of troponin, including the troponin-T linker region that is a hotspot for several cardiomyopathy mutations. This unresolved yet functionally-significant linker region has been proposed to be intrinsically disordered, with behaviors that are not well described by traditional structural approaches; however, this proposal has not been experimentally verified. Here, we used a combination of single-molecule Förster resonance energy transfer (FRET), molecular dynamics simulations, and functional reconstitution assays to investigate the troponin-T linker region. We experimentally and computationally show that in the context of both isolated troponin and the fully regulated troponin complex, the linker behaves as a dynamic, intrinsically disordered region. This region undergoes polyampholyte expansion in the presence of high salt and distinct conformational changes during the assembly of the troponin complex. We also examine the ΔE160 hypertrophic cardiomyopathy mutation in the linker, and we demonstrate that this mutation does not affect the conformational dynamics of the linker, rather it allosterically affects interactions with other subunits of the troponin complex, leading to increased molecular contractility. Taken together, our data clearly demonstrate the importance of disorder within the troponin-T linker and provide new insights into the molecular mechanisms controlling the pathogenesis of cardiomyopathies.

Control of Inflammatory Response by Tissue Microenvironment

Inflammation is an essential defense response but operates at the cost of normal functions. Whether and how the negative impact of inflammation is monitored remains largely unknown. Acidification of the tissue microenvironment is associated with inflammation. Here we investigated whether macrophages sense tissue acidification to adjust inflammatory responses. We found that acidic pH restructured the inflammatory response of macrophages in a gene-specific manner. We identified mammalian BRD4 as a novel intracellular pH sensor. Acidic pH disrupts the transcription condensates containing BRD4 and MED1, via histidine-enriched intrinsically disordered regions. Crucially, decrease in macrophage intracellular pH is necessary and sufficient to regulate transcriptional condensates in vitro and in vivo, acting as negative feedback to regulate the inflammatory response. Collectively, these findings uncovered a pH-dependent switch in transcriptional condensates that enables environmental sensing to directly control inflammation, with a broader implication for calibrating the magnitude and quality of inflammation by the inflammatory cost.

Dynamic modes of Notch transcription hubs conferring memory and stochastic activation revealed by live imaging the co-activator Mastermind

Developmental programming involves the accurate conversion of signalling levels and dynamics to transcriptional outputs. The transcriptional relay in the Notch pathway relies on nuclear complexes containing the co-activator Mastermind (Mam). By tracking these complexes in real time, we reveal that they promote the formation of a dynamic transcription hub in Notch ON nuclei which concentrates key factors including the Mediator CDK module. The composition of the hub is labile and persists after Notch withdrawal conferring a memory that enables rapid reformation. Surprisingly, only a third of Notch ON hubs progress to a state with nascent transcription, which correlates with polymerase II and core Mediator recruitment. This probability is increased by a second signal. The discovery that target-gene transcription is probabilistic has far-reaching implications because it implies that stochastic differences in Notch pathway output can arise downstream of receptor activation.

The origin recognition complex requires chromatin tethering by a hypervariable intrinsically disordered region that is functionally conserved from sponge to man

The first step toward eukaryotic genome duplication is loading of the replicative helicase onto chromatin. This 'licensing' step initiates with the recruitment of the origin recognition complex (ORC) to chromatin, which is thought to occur via ORC's ATP-dependent DNA binding and encirclement activity. However, we have previously shown that ATP binding is dispensable for the chromatin recruitment of fly ORC, raising the question of how metazoan ORC binds chromosomes. We show here that the intrinsically disordered region (IDR) of fly Orc1 is both necessary and sufficient for recruitment of ORC to chromosomes in vivo and demonstrate that this is regulated by IDR phosphorylation. Consistently, we find that the IDR confers the ORC holocomplex with ATP-independent DNA binding activity in vitro. Using phylogenetic analysis, we make the surprising observation that metazoan Orc1 IDRs have diverged so markedly that they are unrecognizable as orthologs and yet we find that these compositionally homologous sequences are functionally conserved. Altogether, these data suggest that chromatin is recalcitrant to ORC's ATP-dependent DNA binding activity, necessitating IDR-dependent chromatin tethering, which we propose poises ORC to opportunistically encircle nucleosome-free regions as they become available.

The temperature sensor TWA1 is required for thermotolerance in Arabidopsis

Plants exposed to incidences of excessive temperatures activate heat-stress responses to cope with the physiological challenge and stimulate long-term acclimation1,2. The mechanism that senses cellular temperature for inducing thermotolerance is still unclear3. Here we show that TWA1 is a temperature-sensing transcriptional co-regulator that is needed for basal and acquired thermotolerance in Arabidopsis thaliana. At elevated temperatures, TWA1 changes its conformation and allows physical interaction with JASMONATE-ASSOCIATED MYC-LIKE (JAM) transcription factors and TOPLESS (TPL) and TOPLESS-RELATED (TPR) proteins for repressor complex assembly. TWA1 is a predicted intrinsically disordered protein that has a key thermosensory role functioning through an amino-terminal highly variable region. At elevated temperatures, TWA1 accumulates in nuclear subdomains, and physical interactions with JAM2 and TPL appear to be restricted to these nuclear subdomains. The transcriptional upregulation of the heat shock transcription factor A2 (HSFA2) and heat shock proteins depended on TWA1, and TWA1 orthologues provided different temperature thresholds, consistent with the sensor function in early signalling of heat stress. The identification of the plant thermosensors offers a molecular tool for adjusting thermal acclimation responses of crops by breeding and biotechnology, and a sensitive temperature switch for thermogenetics.

Disordered clock protein interactions and charge blocks turn an hourglass into a persistent circadian oscillator

Organismal physiology is widely regulated by the molecular circadian clock, a feedback loop composed of protein complexes whose members are enriched in intrinsically disordered regions. These regions can mediate protein-protein interactions via SLiMs, but the contribution of these disordered regions to clock protein interactions had not been elucidated. To determine the functionality of these disordered regions, we applied a synthetic peptide microarray approach to the disordered clock protein FRQ in Neurospora crassa. We identified residues required for FRQ's interaction with its partner protein FRH, the mutation of which demonstrated FRH is necessary for persistent clock oscillations but not repression of transcriptional activity. Additionally, the microarray demonstrated an enrichment of FRH binding to FRQ peptides with a net positive charge. We found that positively charged residues occurred in significant "blocks" within the amino acid sequence of FRQ and that ablation of one of these blocks affected both core clock timing and physiological clock output. Finally, we found positive charge clusters were a commonly shared molecular feature in repressive circadian clock proteins. Overall, our study suggests a mechanistic purpose for positive charge blocks and yielded insights into repressive arm protein roles in clock function.

A fine kinetic balance of interactions directs transcription factor hubs to genes

Eukaryotic gene regulation relies on the binding of sequence-specific transcription factors (TFs). TFs bind chromatin transiently yet occupy their target sites by forming high-local concentration microenvironments (hubs and condensates) that increase the frequency of binding events. Despite their ubiquity, such microenvironments have been difficult to study in endogenous contexts due to technical limitations. Here, we overcome these limitations and investigate how hubs drive TF occupancy at their targets. Using a DNA binding perturbation to a hub-forming TF, Zelda, in Drosophila embryos, we find that hub properties, including the stability and frequencies of associations to targets, are key determinants of TF occupancy. Our data suggest that the targeting of these hubs is driven not just by specific DNA motif recognition, but also by a fine-tuned kinetic balance of interactions between TFs and their co-binding partners.

Labile assembly of a tardigrade protein induces biostasis

Tardigrades are microscopic animals that survive desiccation by inducing biostasis. To survive drying tardigrades rely on intrinsically disordered CAHS proteins, which also function to prevent perturbations induced by drying in vitro and in heterologous systems. CAHS proteins have been shown to form gels both in vitro and in vivo, which has been speculated to be linked to their protective capacity. However, the sequence features and mechanisms underlying gel formation and the necessity of gelation for protection have not been demonstrated. Here we report a mechanism of fibrillization and gelation for CAHS D similar to that of intermediate filament assembly. We show that in vitro, gelation restricts molecular motion, immobilizing and protecting labile material from the harmful effects of drying. In vivo, we observe that CAHS D forms fibrillar networks during osmotic stress. Fibrillar networking of CAHS D improves survival of osmotically shocked cells. We observe two emergent properties associated with fibrillization; (i) prevention of cell volume change and (ii) reduction of metabolic activity during osmotic shock. We find that there is no significant correlation between maintenance of cell volume and survival, while there is a significant correlation between reduced metabolism and survival. Importantly, CAHS D's fibrillar network formation is reversible and metabolic rates return to control levels after CAHS fibers are resolved. This work provides insights into how tardigrades induce reversible biostasis through the self-assembly of labile CAHS gels.

Mutational scanning of CRX classifies clinical variants and reveals biochemical properties of the transcriptional effector domain

Cone-Rod Homeobox, encoded by CRX, is a transcription factor (TF) essential for the terminal differentiation and maintenance of mammalian photoreceptors. Structurally, CRX comprises an ordered DNA-binding homeodomain and an intrinsically disordered transcriptional effector domain. Although a handful of human variants in CRX have been shown to cause several different degenerative retinopathies with varying cone and rod predominance, as with most human disease genes the vast majority of observed CRX genetic variants are uncharacterized variants of uncertain significance (VUS). We performed a deep mutational scan (DMS) of nearly all possible single amino acid substitution variants in CRX, using an engineered cell-based transcriptional reporter assay. We measured the ability of each CRX missense variant to transactivate a synthetic fluorescent reporter construct in a pooled fluorescence-activated cell sorting assay and compared the activation strength of each variant to that of wild-type CRX to compute an activity score, identifying thousands of variants with altered transcriptional activity. We calculated a statistical confidence for each activity score derived from multiple independent measurements of each variant marked by unique sequence barcodes, curating a high-confidence list of nearly 2,000 variants with significantly altered transcriptional activity compared to wild-type CRX. We evaluated the performance of the DMS assay as a clinical variant classification tool using gold-standard classified human variants from ClinVar, and determined that activity scores could be used to identify pathogenic variants with high specificity. That this performance could be achieved using a synthetic reporter assay in a foreign cell type, even for a highly cell type-specific TF like CRX, suggests that this approach shows promise for DMS of other TFs that function in cell types that are not easily accessible. Per-position average activity scores closely aligned to a predicted structure of the ordered homeodomain and demonstrated position-specific residue requirements. The intrinsically disordered transcriptional effector domain, by contrast, displayed a qualitatively different pattern of substitution effects, following compositional constraints without specific residue position requirements in the peptide chain. The observed compositional constraints of the effector domain were consistent with the acidic exposure model of transcriptional activation. Together, the results of the CRX DMS identify molecular features of the CRX effector domain and demonstrate clinical utility for variant classification.

Tissue Usage Preference and Intrinsically Disordered Region Remodeling of Alternative Splicing Derived Proteoforms in the Heart

A computational analysis of mass spectrometry data was performed to uncover alternative splicing derived protein variants across chambers of the human heart. Evidence for 216 non-canonical isoforms was apparent in the atrium and the ventricle, including 52 isoforms not documented on SwissProt and recovered using an RNA sequencing derived database. Among non-canonical isoforms, 29 show signs of regulation based on statistically significant preferences in tissue usage, including a ventricular enriched protein isoform of tensin-1 (TNS1) and an atrium-enriched PDZ and LIM Domain 3 (PDLIM3) isoform 2 (PDLIM3-2/ALP-H). Examined variant regions that differ between alternative and canonical isoforms are highly enriched with intrinsically disordered regions. Moreover, over two-thirds of such regions are predicted to function in protein binding and RNA binding. The analysis here lends further credence to the notion that alternative splicing diversifies the proteome by rewiring intrinsically disordered regions, which are increasingly recognized to play important roles in the generation of biological function from protein sequences.

Direct prediction of intrinsically disordered protein conformational properties from sequence

Intrinsically disordered regions (IDRs) are ubiquitous across all domains of life and play a range of functional roles. While folded domains are generally well described by a stable three-dimensional structure, IDRs exist in a collection of interconverting states known as an ensemble. This structural heterogeneity means that IDRs are largely absent from the Protein Data Bank, contributing to a lack of computational approaches to predict ensemble conformational properties from sequence. Here we combine rational sequence design, large-scale molecular simulations and deep learning to develop ALBATROSS, a deep-learning model for predicting ensemble dimensions of IDRs, including the radius of gyration, end-to-end distance, polymer-scaling exponent and ensemble asphericity, directly from sequences at a proteome-wide scale. ALBATROSS is lightweight, easy to use and accessible as both a locally installable software package and a point-and-click-style interface via Google Colab notebooks. We first demonstrate the applicability of our predictors by examining the generalizability of sequence-ensemble relationships in IDRs. Then, we leverage the high-throughput nature of ALBATROSS to characterize the sequence-specific biophysical behavior of IDRs within and between proteomes.

Disordered proteins interact with the chemical environment to tune their protective function during drying

The conformational ensemble and function of intrinsically disordered proteins (IDPs) are sensitive to their solution environment. The inherent malleability of disordered proteins combined with the exposure of their residues accounts for this sensitivity. One context in which IDPs play important roles that is concomitant with massive changes to the intracellular environment is during desiccation (extreme drying). The ability of organisms to survive desiccation has long been linked to the accumulation of high levels of cosolutes such as trehalose or sucrose as well as the enrichment of IDPs, such as late embryogenesis abundant (LEA) proteins or cytoplasmic abundant heat soluble (CAHS) proteins. Despite knowing that IDPs play important roles and are co-enriched alongside endogenous, species-specific cosolutes during desiccation, little is known mechanistically about how IDP-cosolute interactions influence desiccation tolerance. Here, we test the notion that the protective function of desiccation-related IDPs is enhanced through conformational changes induced by endogenous cosolutes. We find that desiccation-related IDPs derived from four different organisms spanning two LEA protein families and the CAHS protein family, synergize best with endogenous cosolutes during drying to promote desiccation protection. Yet the structural parameters of protective IDPs do not correlate with synergy for either CAHS or LEA proteins. We further demonstrate that for CAHS, but not LEA proteins, synergy is related to self-assembly and the formation of a gel. Our results demonstrate that functional synergy between IDPs and endogenous cosolutes is a convergent desiccation protection strategy seen among different IDP families and organisms, yet, the mechanisms underlying this synergy differ between IDP families.

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.

A multi-faceted approach to unravel coding and non-coding gene fusions and target chimeric proteins in ataxia

Ataxia represents a heterogeneous group of neurodegenerative disorders characterized by a loss of balance and coordination, often resulting from mutations in genes vital for cerebellar function and maintenance. Recent advances in genomics have identified gene fusion events as critical contributors to various cancers and neurodegenerative diseases. However, their role in ataxia pathogenesis remains largely unexplored. Our study Hdelved into this possibility by analyzing RNA sequencing data from 1443 diverse samples, including cell and mouse models, patient samples, and healthy controls. We identified 7067 novel gene fusions, potentially pivotal in disease onset. These fusions, notably in-frame, could produce chimeric proteins, disrupt gene regulation, or introduce new functions. We observed conservation of specific amino acids at fusion breakpoints and identified potential aggregate formations in fusion proteins, known to contribute to ataxia. Through AI-based protein structure prediction, we identified topological changes in three high-confidence fusion proteins-TEN1-ACOX1, PEX14-NMNAT1, and ITPR1-GRID2-which could potentially alter their functions. Subsequent virtual drug screening identified several molecules and peptides with high-affinity binding to fusion sites. Molecular dynamics simulations confirmed the stability of these protein-ligand complexes at fusion breakpoints. Additionally, we explored the role of non-coding RNA fusions as miRNA sponges. One such fusion, RP11-547P4-FLJ33910, showed strong interaction with hsa-miR-504-5p, potentially acting as its sponge. This interaction correlated with the upregulation of hsa-miR-504-5p target genes, some previously linked to ataxia. In conclusion, our study unveils new aspects of gene fusions in ataxia, suggesting their significant role in pathogenesis and opening avenues for targeted therapeutic interventions.Communicated by Ramaswamy H. Sarma.

Leucine Motifs Stabilize Residual Helical Structure in Disordered Proteins

Many examples are known of regions of intrinsically disordered proteins that fold into α-helices upon binding to their targets. These helical binding motifs (HBMs) can be partially helical also in the unbound state, and this so-called residual structure can affect binding affinity and kinetics. To investigate the underlying mechanisms governing the formation of residual helical structure, we assembled a dataset of experimental helix contents of 65 peptides containing HBM that fold-upon-binding. The average residual helicity is 17% and increases to 60% upon target binding. The helix contents of residual and target-bound structures do not correlate, however the relative location of helix elements in both states shows a strong overlap. Compared to the general disordered regions, HBMs are enriched in amino acids with high helix preference and these residues are typically involved in target binding, explaining the overlap in helix positions. In particular, we find that leucine residues and leucine motifs in HBMs are the major contributors to helix stabilization and target-binding. For the two model peptides, we show that substitution of leucine motifs to other hydrophobic residues (valine or isoleucine) leads to reduction of residual helicity, supporting the role of leucine as helix stabilizer. From the three hydrophobic residues only leucine can efficiently stabilize residual helical structure. We suggest that the high occurrence of leucine motifs and a general preference for leucine at binding interfaces in HBMs can be explained by its unique ability to stabilize helical elements.

Alternative conformations of a group 4 Late Embryogenesis Abundant protein associated to its in vitro protective activity

Late Embryogenesis Abundant (LEA) proteins are a group of intrinsically disordered proteins implicated in plant responses to water deficit. In vitro studies revealed that LEA proteins protect reporter enzymes from inactivation during low water availability. Group 4 LEA proteins constitute a conserved protein family, displaying in vitro protective capabilities. Under water deficiency or macromolecular crowding, the N-terminal of these proteins adopts an alpha-helix conformation. This region has been identified as responsible for the protein in vitro protective activity. This study investigates whether the attainment of alpha-helix conformation and/or particular amino acid residues are required for the in vitro protective activity. The LEA4-5 protein from Arabidopsis thaliana was used to generate mutant proteins. The mutations altered conserved residues, deleted specific conserved regions, or introduced prolines to hinder alpha-helix formation. The results indicate that conserved residues are not essential for LEA4-5 protective function. Interestingly, the C-terminal region was found to contribute to this function. Moreover, alpha-helix conformation is necessary for the protective activity only when the C-terminal region is deleted. Overall, LEA4-5 shows the ability to adopt alternative functional conformations under the tested conditions. These findings shed light on the in vitro mechanisms by which LEA proteins protect against water deficit stress.

Helicity of a tardigrade disordered protein contributes to its protective function during desiccation

To survive extreme drying (anhydrobiosis), many organisms, spanning every kingdom of life, accumulate intrinsically disordered proteins (IDPs). For decades, the ability of anhydrobiosis-related IDPs to form transient amphipathic helices has been suggested to be important for promoting desiccation tolerance. However, evidence empirically supporting the necessity and/or sufficiency of helicity in mediating anhydrobiosis is lacking. Here, we demonstrate that the linker region of CAHS D, a desiccation-related IDP from the tardigrade Hypsibius exemplaris, that contains significant helical structure, is the protective portion of this protein. Perturbing the sequence composition and grammar of the linker region of CAHS D, through the insertion of helix-breaking prolines, modulating the identity of charged residues, or replacement of hydrophobic amino acids with serine or glycine residues results in variants with different degrees of helical structure. Importantly, correlation of protective capacity and helical content in variants generated through different helix perturbing modalities does not show as strong a trend, suggesting that while helicity is important, it is not the only property that makes a protein protective during desiccation. These results provide direct evidence for the decades-old theory that helicity of desiccation-related IDPs is linked to their anhydrobiotic capacity.

The male pachynema-specific protein MAPS drives phase separation in vitro and regulates sex body formation and chromatin behaviors in vivo

Dynamic chromosome remodeling and nuclear compartmentalization take place during mammalian meiotic prophase I. We report here that the crucial roles of male pachynema-specific protein (MAPS) in pachynema progression might be mediated by its liquid-liquid phase separation in vitro and in cellulo. MAPS forms distinguishable liquid phases, and deletion or mutations of its N-terminal amino acids (aa) 2-9 disrupt its secondary structure and charge properties, impeding phase separation. Maps-/- pachytene spermatocytes exhibit defects in nucleus compartmentalization, including defects in forming sex bodies, altered nucleosome composition, and disordered chromatin accessibility. MapsΔ2-9/Δ2-9 male mice expressing MAPS protein lacking aa 2-9 phenocopy Maps-/- mice. Moreover, a frameshift mutation in C3orf62, the human counterpart of Maps, is correlated with nonobstructive azoospermia in a patient exhibiting pachynema arrest in spermatocyte development. Hence, the phase separation property of MAPS seems essential for pachynema progression in mouse and human spermatocytes.

Rapid and reversible dissolution of biomolecular condensates using light-controlled recruitment of a solubility tag

Biomolecular condensates are broadly implicated in both normal cellular regulation and disease. Consequently, several chemical biology and optogenetic approaches have been developed to induce phase separation of a protein of interest. However, few tools are available to perform the converse function-dissolving a condensate of interest on demand. Such a tool would aid in testing whether the condensate plays specific functional roles, a major question in cell biology and drug development. Here we report an optogenetic approach to selectively dissolve a condensate of interest in a reversible and spatially controlled manner. We show that light-gated recruitment of maltose-binding protein (MBP), a commonly used solubilizing domain in protein purification, results in rapid and controlled dissolution of condensates formed from proteins of interest. Our optogenetic MBP-based dissolution strategy (OptoMBP) is rapid, reversible, and can be spatially controlled with subcellular precision. We also provide a proof-of-principle application of OptoMBP, showing that disrupting condensation of the oncogenic fusion protein FUS-CHOP results in reversion of FUS-CHOP driven transcriptional changes. We envision that the OptoMBP system could be broadly useful for disrupting constitutive protein condensates to probe their biological functions.

Molecular dynamics simulations revealed topological frustration in the binding-wrapping process of eIF4G with eIF4E

Interaction between the cap-binding protein eIF4E and the scaffolding protein eIF4G is essential for the cap-dependent translation initiation in eukaryotes. In the Saccharomyces cerevisiae eIF4G/eIF4E complex, the intrinsically disordered eIF4E-binding domain of eIF4G folds into a bracelet-like structure upon binding to eIF4E. Aiming to unveil the molecular mechanism underlying the binding-wrapping process of eIF4G with eIF4E, we performed extensive coarse-grained molecular dynamics simulations and transition path analysis in this work. The major transition pathway revealed from our simulations showed that docking of the eIF4E-binding motif of eIF4G to the folded core of eIF4E initiates the binding process and then the disordered eIF4G wraps around the N-terminal tail of eIF4E. Additionally, we identified a minor transition pathway which indicates the involvement of topological frustration in the binding process. By manipulating the interaction strength of the wrapping contacts and the latching contacts, we further dissected factors affecting the formation of topological frustration and the binding transition kinetics. Our findings provide new clues for experimental studies on the binding mechanism of eIF4G to eIF4E in the future and exemplify the involvement of topological frustration in the binding process of intrinsically disordered proteins.

Pathogenic Variants in USH1G/SANS Alter Protein Interaction with Pre-RNA Processing Factors PRPF6 and PRPF31 of the Spliceosome

Pre-mRNA splicing is an essential process orchestrated by the spliceosome, a dynamic complex assembled stepwise on pre-mRNA. We have previously identified that USH1G protein SANS regulates pre-mRNA splicing by mediating the intranuclear transfer of the spliceosomal U4/U6.U5 tri-snRNP complex. During this process, SANS interacts with the U4/U6 and U5 snRNP-specific proteins PRPF31 and PRPF6 and regulates splicing, which is disturbed by variants of USH1G/SANS causative for human Usher syndrome (USH), the most common form of hereditary deaf-blindness. Here, we aim to gain further insights into the molecular interaction of the splicing molecules PRPF31 and PRPF6 to the CENTn domain of SANS using fluorescence resonance energy transfer assays in cells and in silico deep learning-based protein structure predictions. This demonstrates that SANS directly binds via two distinct conserved regions of its CENTn to the two PRPFs. In addition, we provide evidence that these interactions occur sequentially and a conformational change of an intrinsically disordered region to a short α-helix of SANS CENTn2 is triggered by the binding of PRPF6. Furthermore, we find that pathogenic variants of USH1G/SANS perturb the binding of SANS to both PRPFs, implying a significance for the USH1G pathophysiology.

FLiPPR: A Processor for Limited Proteolysis (LiP) Mass Spectrometry Datasets Built on FragPipe

Here, we present FLiPPR, or FragPipe LiP (limited proteolysis) Processor, a tool that facilitates the analysis of data from limited proteolysis mass spectrometry (LiP-MS) experiments following primary search and quantification in FragPipe. LiP-MS has emerged as a method that can provide proteome-wide information on protein structure and has been applied to a range of biological and biophysical questions. Although LiP-MS can be carried out with standard laboratory reagents and mass spectrometers, analyzing the data can be slow and poses unique challenges compared to typical quantitative proteomics workflows. To address this, we leverage the fast, sensitive, and accurate search and label-free quantification algorithms in FragPipe and then process its output in FLiPPR. FLiPPR formalizes a specific data imputation heuristic that carefully uses missing data in LiP-MS experiments to report on the most significant structural changes. Moreover, FLiPPR introduces a new data merging scheme (from ions to cut-sites) and a protein-centric multiple hypothesis correction scheme, collectively enabling processed LiP-MS datasets to be more robust and less redundant. These improvements substantially strengthen statistical trends when previously published data are reanalyzed with the FragPipe/FLiPPR workflow. As a final feature, FLiPPR facilitates the collection of structural metadata to identify correlations between experiments and structural features. We hope that FLiPPR will lower the barrier for more users to adopt LiP-MS, standardize statistical procedures for LiP-MS data analysis, and systematize output to facilitate eventual larger-scale integration of LiP-MS data.

Protein intrinsically disordered regions have a non-random, modular architecture

MOTIVATION

Protein sequences can be broadly categorized into two classes: those which adopt stable secondary structure and fold into a domain (i.e. globular proteins), and those that do not. The sequences belonging to this latter class are conformationally heterogeneous and are described as being intrinsically disordered. Decades of investigation into the structure and function of globular proteins has resulted in a suite of computational tools that enable their sub-classification by domain type, an approach that has revolutionized how we understand and predict protein functionality. Conversely, it is unknown if sequences of disordered protein regions are subject to broadly generalizable organizational principles that would enable their sub-classification.

RESULTS

Here, we report the development of a statistical approach that quantifies linear variance in amino acid composition across a sequence. With multiple examples, we provide evidence that intrinsically disordered regions are organized into statistically non-random modules of unique compositional bias. Modularity is observed for both low and high-complexity sequences and, in some cases, we find that modules are organized in repetitive patterns. These data demonstrate that disordered sequences are non-randomly organized into modular architectures and motivate future experiments to comprehensively classify module types and to determine the degree to which modules constitute functionally separable units analogous to the domains of globular proteins.

AVAILABILITY AND IMPLEMENTATION

The source code, documentation, and data to reproduce all figures are freely available at https://github.com/MWPlabUTSW/Chi-Score-Analysis.git. The analysis is also available as a Google Colab Notebook (https://colab.research.google.com/github/MWPlabUTSW/Chi-Score-Analysis/blob/main/ChiScore_Analysis.ipynb).

Structure of the N-RNA/P interface indicates mode of L/P recruitment to the nucleocapsid of human metapneumovirus

Human metapneumovirus (HMPV) is a major cause of respiratory illness in young children. The HMPV polymerase (L) binds an obligate cofactor, the phosphoprotein (P). During replication and transcription, the L/P complex traverses the viral RNA genome, which is encapsidated within nucleoproteins (N). An essential interaction between N and a C-terminal region of P tethers the L/P polymerase to the template. This N-P interaction is also involved in the formation of cytoplasmic viral factories in infected cells, called inclusion bodies. To define how the polymerase component P recognizes N-encapsidated RNA (N-RNA) we employed cryogenic electron microscopy (cryo-EM) and molecular dynamics simulations, coupled to activity assays and imaging of inclusion bodies in cells. We report a 2.9 Å resolution structure of a triple-complex between multimeric N, bound to both RNA and the C-terminal region of P. Furthermore, we also present cryo-EM structures of assembled N in different oligomeric states, highlighting the plasticity of N. Combined with our functional assays, these structural data delineate in molecular detail how P attaches to N-RNA whilst retaining substantial conformational dynamics. Moreover, the N-RNA-P triple complex structure provides a molecular blueprint for the design of therapeutics to potentially disrupt the attachment of L/P to its template.

JL. Additional principles that govern the release of pre-ribosomes from the nucleolus into the nucleoplasm in yeast

During eukaryotic ribosome biogenesis, pre-ribosomes travel from the nucleolus, where assembly is initiated, to the nucleoplasm and then are exported to the cytoplasm, where assembly concludes. Although nuclear export of pre-ribosomes has been extensively investigated, the release of pre-ribosomes from the nucleolus is an understudied phenomenon. Initial data indicate that unfolded rRNA interacts in trans with nucleolar components and that, when rRNA folds due to ribosomal protein (RP) binding, the number of trans interactions drops below the threshold necessary for nucleolar retention. To validate and expand on this idea, we performed a bioinformatic analysis of the protein components of the Saccharomyces cerevisiae ribosome assembly pathway. We found that ribosome biogenesis factors (RiBi factors) contain significantly more predicted trans interacting regions than RPs. We also analyzed cryo-EM structures of ribosome assembly intermediates to determine how nucleolar pre-ribosomes differ from post-nucleolar pre-ribosomes, specifically the capacity of RPs, RiBi factors, and rRNA components to interact in trans. We observed a significant decrease in the theoretical trans-interacting capability of pre-ribosomes between nucleolar and post-nucleolar stages of assembly due to the release of RiBi factors from particles and the folding of rRNA. Here, we provide a mechanism for the release of pre-ribosomes from the nucleolus.

Tutorial: a guide for the selection of fast and accurate computational tools for the prediction of intrinsic disorder in proteins

Intrinsic disorder is instrumental for a wide range of protein functions, and its analysis, using computational predictions from primary structures, complements secondary and tertiary structure-based approaches. In this Tutorial, we provide an overview and comparison of 23 publicly available computational tools with complementary parameters useful for intrinsic disorder prediction, partly relying on results from the Critical Assessment of protein Intrinsic Disorder prediction experiment. We consider factors such as accuracy, runtime, availability and the need for functional insights. The selected tools are available as web servers and downloadable programs, offer state-of-the-art predictions and can be used in a high-throughput manner. We provide examples and instructions for the selected tools to illustrate practical aspects related to the submission, collection and interpretation of predictions, as well as the timing and their limitations. We highlight two predictors for intrinsically disordered proteins, flDPnn as accurate and fast and IUPred as very fast and moderately accurate, while suggesting ANCHOR2 and MoRFchibi as two of the best-performing predictors for intrinsically disordered region binding. We link these tools to additional resources, including databases of predictions and web servers that integrate multiple predictive methods. Altogether, this Tutorial provides a hands-on guide to comparatively evaluating multiple predictors, submitting and collecting their own predictions, and reading and interpreting results. It is suitable for experimentalists and computational biologists interested in accurately and conveniently identifying intrinsic disorder, facilitating the functional characterization of the rapidly growing collections of protein sequences.

Systematic identification of conditionally folded intrinsically disordered regions by AlphaFold2

The AlphaFold Protein Structure Database contains predicted structures for millions of proteins. For the majority of human proteins that contain intrinsically disordered regions (IDRs), which do not adopt a stable structure, it is generally assumed that these regions have low AlphaFold2 confidence scores that reflect low-confidence structural predictions. Here, we show that AlphaFold2 assigns confident structures to nearly 15% of human IDRs. By comparison to experimental NMR data for a subset of IDRs that are known to conditionally fold (i.e., upon binding or under other specific conditions), we find that AlphaFold2 often predicts the structure of the conditionally folded state. Based on databases of IDRs that are known to conditionally fold, we estimate that AlphaFold2 can identify conditionally folding IDRs at a precision as high as 88% at a 10% false positive rate, which is remarkable considering that conditionally folded IDR structures were minimally represented in its training data. We find that human disease mutations are nearly fivefold enriched in conditionally folded IDRs over IDRs in general and that up to 80% of IDRs in prokaryotes are predicted to conditionally fold, compared to less than 20% of eukaryotic IDRs. These results indicate that a large majority of IDRs in the proteomes of human and other eukaryotes function in the absence of conditional folding, but the regions that do acquire folds are more sensitive to mutations. We emphasize that the AlphaFold2 predictions do not reveal functionally relevant structural plasticity within IDRs and cannot offer realistic ensemble representations of conditionally folded IDRs.

The fitness cost of spurious phosphorylation

The fidelity of signal transduction requires the binding of regulatory molecules to their cognate targets. However, the crowded cell interior risks off-target interactions between proteins that are functionally unrelated. How such off-target interactions impact fitness is not generally known, but quantifying this is required to understand the constraints faced by cell systems as they evolve. Here, we use the model organism S. cerevisiae to inducibly express tyrosine kinases. Because yeast lacks bona fide tyrosine kinases, most of the resulting tyrosine phosphorylation is spurious. This provides a suitable system to measure the impact of artificial protein interactions on fitness. We engineered 44 yeast strains each expressing a tyrosine kinase, and quantitatively analysed their phosphoproteomes. This analysis resulted in ~30,000 phosphosites mapping to ~3,500 proteins. Examination of the fitness costs in each strain revealed a strong correlation between the number of spurious pY sites and decreased growth. Moreover, the analysis of pY effects on protein structure and on protein function revealed over 1000 pY events that we predict to be deleterious. However, we also find that a large number of the spurious pY sites have a negligible effect on fitness, possibly because of their low stoichiometry. This result is consistent with our evolutionary analyses demonstrating a lack of phosphotyrosine counter-selection in species with bona fide tyrosine kinases. Taken together, our results suggest that, alongside the risk for toxicity, the cell can tolerate a large degree of non-functional crosstalk as interaction networks evolve.

Tissue Usage Preference and Intrinsically Disordered Region Remodeling of Alternative Splicing Derived Proteoforms in the Heart

A computational analysis of mass spectrometry data was performed to uncover alternative splicing derived protein variants across chambers of the human heart. Evidence for 216 non-canonical isoforms was apparent in the atrium and the ventricle, including 52 isoforms not documented on SwissProt and recovered using an RNA sequencing derived database. Among non-canonical isoforms, 29 show signs of regulation based on statistically significant preferences in tissue usage, including a ventricular enriched protein isoform of tensin-1 (TNS1) and an atrium-enriched PDZ and LIM Domain 3 (PDLIM3) isoform 2 (PDLIM3-2/ALP-H). Examined variant regions that differ between alternative and canonical isoforms are highly enriched in intrinsically disordered regions, and over two-thirds of such regions are predicted to function in protein binding and/or RNA binding. The analysis here lends further credence to the notion that alternative splicing diversifies the proteome by rewiring intrinsically disordered regions, which are increasingly recognized to play important roles in the generation of biological function from protein sequences.

Direct ionic stress sensing and mitigation by the transcription factor NFAT5

Homeostatic control of intracellular ionic strength is essential for protein, organelle and genome function, yet mechanisms that sense and enable adaptation to ionic stress remain poorly understood in animals. We find that the transcription factor NFAT5 directly senses solution ionic strength using a C-terminal intrinsically disordered region. Both in intact cells and in a purified system, NFAT5 forms dynamic, reversible biomolecular condensates in response to increasing ionic strength. This self-associative property, conserved from insects to mammals, allows NFAT5 to accumulate in the nucleus and activate genes that restore cellular ion content. Mutations that reduce condensation or those that promote aggregation both reduce NFAT5 activity, highlighting the importance of optimally tuned associative interactions. Remarkably, human NFAT5 alone is sufficient to reconstitute a mammalian transcriptional response to ionic or hypertonic stress in yeast. Thus NFAT5 is both the sensor and effector of a cell-autonomous ionic stress response pathway in animal cells.

Evaluation of the determinants for improved pluripotency induction and maintenance by engineered SOX17

An engineered SOX17 variant with point mutations within its DNA binding domain termed SOX17FNV is a more potent pluripotency inducer than SOX2, yet the underlying mechanism remains unclear. Although wild-type SOX17 was incapable of inducing pluripotency, SOX17FNV outperformed SOX2 in mouse and human pluripotency reprogramming. In embryonic stem cells, SOX17FNV could replace SOX2 to maintain pluripotency despite considerable sequence differences and upregulated genes expressed in cleavage-stage embryos. Mechanistically, SOX17FNV co-bound OCT4 more cooperatively than SOX2 in the context of the canonical SoxOct DNA element. SOX2, SOX17, and SOX17FNV were all able to bind nucleosome core particles in vitro, which is a prerequisite for pioneer transcription factors. Experiments using purified proteins and in cellular contexts showed that SOX17 variants phase-separated more efficiently than SOX2, suggesting an enhanced ability to self-organise. Systematic deletion analyses showed that the N-terminus of SOX17FNV was dispensable for its reprogramming activity. However, the C-terminus encodes essential domains indicating multivalent interactions that drive transactivation and reprogramming. We defined a minimal SOX17FNV (miniSOX) that can support reprogramming with high activity, reducing the payload of reprogramming cassettes. This study uncovers the mechanisms behind SOX17FNV-induced pluripotency and establishes engineered SOX factors as powerful cell engineering tools.

Systematic identification of transcriptional activator domains from non-transcription factor proteins in plants and yeast

Transcription factors promote gene expression via trans-regulatory activation domains. Although whole genome scale screens in model organisms (e.g. human, yeast, fly) have helped identify activation domains from transcription factors, such screens have been less extensively used to explore the occurrence of activation domains in non-transcription factor proteins, such as transcriptional coactivators, chromatin regulators and some cytosolic proteins, leaving a blind spot on what role activation domains in these proteins could play in regulating transcription. We utilized the activation domain predictor PADDLE to mine the entire proteomes of two model eukaryotes, Arabidopsis thaliana and Saccharomyces cerevisiae ( 1 ). We characterized 18,000 fragments covering predicted activation domains from >800 non-transcription factor genes in both species, and experimentally validated that 89% of proteins contained fragments capable of activating transcription in yeast. Peptides with similar sequence composition show a broad range of activities, which is explained by the arrangement of key amino acids. We also annotated hundreds of nuclear proteins with activation domains as putative coactivators; many of which have never been ascribed any function in plants. Furthermore, our library contains >250 non-nuclear proteins containing peptides with activation domain function across both eukaryotic lineages, suggesting that there are unknown biological roles of these peptides beyond transcription. Finally, we identify and validate short, 'universal' eukaryotic activation domains that activate transcription in both yeast and plants with comparable or stronger performance to state-of-the-art activation domains. Overall, our dual host screen provides a blueprint on how to systematically discover novel genetic parts for synthetic biology that function across a wide diversity of eukaryotes.

Significance Statement

Activation domains promote transcription and play a critical role in regulating gene expression. Although the mapping of activation domains from transcription factors has been carried out in previous genome-wide screens, their occurrence in non-transcription factors has been less explored. We utilize an activation domain predictor to mine the entire proteomes of Arabidopsis thaliana and Saccharomyces cerevisiae for new activation domains on non-transcription factor proteins. We validate peptides derived from >750 non-transcription factor proteins capable of activating transcription, discovering many potentially new coactivators in plants. Importantly, we identify novel genetic parts that can function across both species, representing unique synthetic biology tools.

The self-association equilibrium of DNAJA2 regulates its interaction with unfolded substrate proteins and with Hsc70

J-domain proteins tune the specificity of Hsp70s, engaging them in precise functions. Despite their essential role, the structure and function of many J-domain proteins remain largely unknown. We explore human DNAJA2, finding that it reversibly forms highly-ordered, tubular structures that can be dissociated by Hsc70, the constitutively expressed Hsp70 isoform. Cryoelectron microscopy and mutational studies reveal that different domains are involved in self-association. Oligomer dissociation into dimers potentiates its interaction with unfolded client proteins. The J-domains are accessible to Hsc70 within the tubular structure. They allow binding of closely spaced Hsc70 molecules that could be transferred to the unfolded substrate for its cooperative remodelling, explaining the efficient recovery of DNAJA2-bound clients. The disordered C-terminal domain, comprising the last 52 residues, regulates its holding activity and productive interaction with Hsc70. These in vitro findings suggest that the association equilibrium of DNAJA2 could regulate its interaction with client proteins and Hsc70.

ParSe 2.0: A web tool to identify drivers of protein phase separation at the proteome level

We have developed an algorithm, ParSe, which accurately identifies from the primary sequence those protein regions likely to exhibit physiological phase separation behavior. Originally, ParSe was designed to test the hypothesis that, for flexible proteins, phase separation potential is correlated to hydrodynamic size. While our results were consistent with that idea, we also found that many different descriptors could successfully differentiate between three classes of protein regions: folded, intrinsically disordered, and phase-separating intrinsically disordered. Consequently, numerous combinations of amino acid property scales can be used to make robust predictions of protein phase separation. Built from that finding, ParSe 2.0 uses an optimal set of property scales to predict domain-level organization and compute a sequence-based prediction of phase separation potential. The algorithm is fast enough to scan the whole of the human proteome in minutes on a single computer and is equally or more accurate than other published predictors in identifying proteins and regions within proteins that drive phase separation. Here, we describe a web application for ParSe 2.0 that may be accessed through a browser by visiting https://stevewhitten.github.io/Parse_v2_FASTA to quickly identify phase-separating proteins within large sequence sets, or by visiting https://stevewhitten.github.io/Parse_v2_web to evaluate individual protein sequences.

CEBPA phase separation links transcriptional activity and 3D chromatin hubs

Cell identity is orchestrated through an interplay between transcription factor (TF) action and genome architecture. The mechanisms used by TFs to shape three-dimensional (3D) genome organization remain incompletely understood. Here we present evidence that the lineage-instructive TF CEBPA drives extensive chromatin compartment switching and promotes the formation of long-range chromatin hubs during induced B cell-to-macrophage transdifferentiation. Mechanistically, we find that the intrinsically disordered region (IDR) of CEBPA undergoes in vitro phase separation (PS) dependent on aromatic residues. Both overexpressing B cells and native CEBPA-expressing cell types such as primary granulocyte-macrophage progenitors, liver cells, and trophectoderm cells reveal nuclear CEBPA foci and long-range 3D chromatin hubs at CEBPA-bound regions. In short, we show that CEBPA can undergo PS through its IDR, which may underlie in vivo foci formation and suggest a potential role of PS in regulating CEBPA function.

Intrinsic disorder in PRAME and its role in uveal melanoma

INTRODUCTION

The PReferentially expressed Antigen in MElanoma (PRAME) protein has been shown to be an independent biomarker for increased risk of metastasis in Class 1 uveal melanomas (UM). Intrinsically disordered proteins and regions of proteins (IDPs/IDPRs) are proteins that do not have a well-defined three-dimensional structure and have been linked to neoplastic development. Our study aimed to evaluate the presence of intrinsic disorder in PRAME and the role these structureless regions have in PRAME( +) Class 1 UM.

METHODS

A bioinformatics study to characterize PRAME's propensity for the intrinsic disorder. We first used the AlphaFold tool to qualitatively assess the protein structure of PRAME. Then we used the Compositional Profiler and a set of per-residue intrinsic disorder predictors to quantify the intrinsic disorder. The Database of Disordered Protein Prediction (D2P2) platform, IUPred, FuzDrop, fIDPnn, AUCpred, SPOT-Disorder2, and metapredict V2 allowed us to evaluate the potential functional disorder of PRAME. Additionally, we used the Search Tool for the Retrieval of Interacting Genes (STRING) to analyze PRAME's potential interactions with other proteins.

RESULTS

Our structural analysis showed that PRAME contains intrinsically disordered protein regions (IDPRs), which are structureless and flexible. We found that PRAME is significantly enriched with serine (p-value < 0.05), a disorder-promoting amino acid. PRAME was found to have an average disorder score of 16.49% (i.e., moderately disordered) across six per-residue intrinsic disorder predictors. Our IUPred analysis revealed the presence of disorder-to-order transition (DOT) regions in PRAME near the C-terminus of the protein (residues 475-509). The D2P2 platform predicted a region from approximately 140 and 175 to be highly concentrated with post-translational modifications (PTMs). FuzDrop predicted the PTM hot spot of PRAME to be a droplet-promoting region and an aggregation hotspot. Finally, our analysis using the STRING tool revealed that PRAME has significantly more interactions with other proteins than expected for randomly selected proteins of the same size, with the ability to interact with 84 different partners (STRING analysis result: p-value < 1.0 × 10-16; model confidence: 0.400).

CONCLUSION

Our study revealed that PRAME has IDPRs that are possibly linked to its functionality in the context of Class 1 UM. The regions of functionality (i.e., DOT regions, PTM sites, droplet-promoting regions, and aggregation hotspots) are localized to regions of high levels of disorder. PRAME has a complex protein-protein interaction (PPI) network that may be secondary to the structureless features of the polypeptide. Our findings contribute to our understanding of UM and suggest that IDPRs and DOT regions in PRAME may be targeted in developing new therapies for this aggressive cancer. Video Abstract.

Hierarchical TAF1-dependent co-translational assembly of the basal transcription factor TFIID

Large heteromeric multiprotein complexes play pivotal roles at every step of gene expression in eukaryotic cells. Among them, the 20-subunit basal transcription factor TFIID nucleates the RNA polymerase II preinitiation complex at gene promoters. Here, by combining systematic RNA-immunoprecipitation (RIP) experiments, single-molecule imaging, proteomics and structure-function analyses, we show that human TFIID biogenesis occurs co-translationally. We discovered that all protein heterodimerization steps happen during protein synthesis. We identify TAF1-the largest protein in the complex-as a critical factor for TFIID assembly. TAF1 acts as a flexible scaffold that drives the co-translational recruitment of TFIID submodules preassembled in the cytoplasm. Altogether, our data suggest a multistep hierarchical model for TFIID biogenesis that culminates with the co-translational assembly of the complex onto the nascent TAF1 polypeptide. We envision that this assembly strategy could be shared with other large heteromeric protein complexes.

Transcription factors interact with RNA to regulate genes

Transcription factors (TFs) orchestrate the gene expression programs that define each cell's identity. The canonical TF accomplishes this with two domains, one that binds specific DNA sequences and the other that binds protein coactivators or corepressors. We find that at least half of TFs also bind RNA, doing so through a previously unrecognized domain with sequence and functional features analogous to the arginine-rich motif of the HIV transcriptional activator Tat. RNA binding contributes to TF function by promoting the dynamic association between DNA, RNA, and TF on chromatin. TF-RNA interactions are a conserved feature important for vertebrate development and disrupted in disease. We propose that the ability to bind DNA, RNA, and protein is a general property of many TFs and is fundamental to their gene regulatory function.