Protein-Protein Interactions Mediated by Intrinsically Disordered Protein Regions Are Enriched in Missense Mutations

Proteome-scale characterisation of motif-based interactome rewiring by disease mutations

Whole genome and exome sequencing are reporting on hundreds of thousands of missense mutations. Taking a pan-disease approach, we explored how mutations in intrinsically disordered regions (IDRs) break or generate protein interactions mediated by short linear motifs. We created a peptide-phage display library tiling ~57,000 peptides from the IDRs of the human proteome overlapping 12,301 single nucleotide variants associated with diverse phenotypes including cancer, metabolic diseases and neurological diseases. By screening 80 human proteins, we identified 366 mutation-modulated interactions, with half of the mutations diminishing binding, and half enhancing binding or creating novel interaction interfaces. The effects of the mutations were confirmed by affinity measurements. In cellular assays, the effects of motif-disruptive mutations were validated, including loss of a nuclear localisation signal in the cell division control protein CDC45 by a mutation associated with Meier-Gorlin syndrome. The study provides insights into how disease-associated mutations may perturb and rewire the motif-based interactome.

Interaction between Yersinia pestis Ail Outer Membrane Protein and the C-Terminal Domain of Human Vitronectin

Yersinia pestis, the causative agent of plague, is capable of evading the human immune system response by recruiting the plasma circulating vitronectin proteins, which act as a shield and avoid its lysis. Vitronectin recruitment is mediated by its interaction with the bacterial transmembrane protein Ail, protruding from the Y. pestis outer membrane. By using all-atom long-scale molecular dynamic simulations of Ail embedded in a realistic model of the bacterial membrane, we have shown that vitronectin forms a stable complex, mediated by interactions between the disordered moieties of the two proteins. The main amino acids driving the complexation have also been evidenced, thus favoring the possible rational design of specific peptides which, by inhibiting vitronectin recruitment, could act as original antibacterial agents.

Pathogenic mutations of human phosphorylation sites affect protein-protein interactions

Despite their lack of a defined 3D structure, intrinsically disordered regions (IDRs) of proteins play important biological roles. Many IDRs contain short linear motifs (SLiMs) that mediate protein-protein interactions (PPIs), which can be regulated by post-translational modifications like phosphorylation. 20% of pathogenic missense mutations are found in IDRs, and understanding how such mutations affect PPIs is essential for unraveling disease mechanisms. Here, we employ peptide-based interaction proteomics to investigate 36 disease-associated mutations affecting phosphorylation sites. Our results unveil significant differences in interactomes between phosphorylated and non-phosphorylated peptides, often due to disrupted phosphorylation-dependent SLiMs. We focused on a mutation of a serine phosphorylation site in the transcription factor GATAD1, which causes dilated cardiomyopathy. We find that this phosphorylation site mediates interaction with 14-3-3 family proteins. Follow-up experiments reveal the structural basis of this interaction and suggest that 14-3-3 binding affects GATAD1 nucleocytoplasmic transport by masking a nuclear localisation signal. Our results demonstrate that pathogenic mutations of human phosphorylation sites can significantly impact protein-protein interactions, offering insights into potential molecular mechanisms underlying pathogenesis.

Cellular targets and lysine selectivity of the HERC5 ISG15 ligase

ISG15 is a type I interferon-induced ubiquitin-like modifier that functions in innate immune responses. The major human ISG15 ligase is hHERC5, a ribosome-associated HECT E3 that broadly ISGylates proteins cotranslationally. Here, we characterized the hHERC5-dependent ISGylome and identified over 2,000 modified lysines in over 1,100 proteins in IFN-β-stimulated cells. In parallel, we compared the substrate selectivity hHERC5 to the major mouse ISG15 ligase, mHERC6, and analysis of sequences surrounding ISGylation sites revealed that hHERC5 and mHERC6 have distinct preferences for amino acid sequence context. Several features of the datasets were consistent with ISGylation of ribosome-tethered nascent chains, and mHERC6, like hHERC5, cotranslationally modified nascent polypeptides. The ISGylome datasets presented here represent the largest numbers of protein targets and modification sites attributable to a single Ub/Ubl ligase and the lysine selectivities of the hHERC5 and mHERC6 enzymes may have implications for the activities of HECT domain ligases, generally.

Dynamic nucleolar phase separation influenced by non-canonical function of LIN28A instructs pluripotent stem cell fate decisions

LIN28A is important in somatic reprogramming and pluripotency regulation. Although previous studies addressed that LIN28A can repress let-7 microRNA maturation in the cytoplasm, few focused on its role within the nucleus. Here, we show that the nucleolus-localized LIN28A protein undergoes liquid-liquid phase separation (LLPS) in mouse embryonic stem cells (mESCs) and in vitro. The RNA binding domains (RBD) and intrinsically disordered regions (IDR) of LIN28A contribute to LIN28A and the other nucleolar proteins' phase-separated condensate establishment. S120A, S200A and R192G mutations in the IDR result in subcellular mislocalization of LIN28A and abnormal nucleolar phase separation. Moreover, we find that the naive-to-primed pluripotency state conversion and the reprogramming are associated with dynamic nucleolar remodeling, which depends on LIN28A's phase separation capacity, because the LIN28A IDR point mutations abolish its role in regulating nucleolus and in these cell fate decision processes, and an exogenous IDR rescues it. These findings shed light on the nucleolar function in pluripotent stem cell states and on a non-canonical RNA-independent role of LIN28A in phase separation and cell fate decisions.

CaMKK2 and CHK1 phosphorylate human STN1 in response to replication stress to protect stalled forks from aberrant resection

Keeping replication fork stable is essential for safeguarding genome integrity; hence, its protection is highly regulated. The CTC1-STN1-TEN1 (CST) complex protects stalled forks from aberrant MRE11-mediated nascent strand DNA degradation (NSD). However, the activation mechanism for CST at forks is unknown. Here, we report that STN1 is phosphorylated in its intrinsic disordered region. Loss of STN1 phosphorylation reduces the replication stress-induced STN1 localization to stalled forks, elevates NSD, increases MRE11 access to stalled forks, and decreases RAD51 localization at forks, leading to increased genome instability under perturbed DNA replication condition. STN1 is phosphorylated by both the ATR-CHK1 and the calcium-sensing kinase CaMKK2 in response to hydroxyurea/aphidicolin treatment or elevated cytosolic calcium concentration. Cancer-associated STN1 variants impair STN1 phosphorylation, conferring inability of fork protection. Collectively, our study uncovers that CaMKK2 and ATR-CHK1 target STN1 to enable its fork protective function, and suggests an important role of STN1 phosphorylation in cancer development.

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 intrinsically disordered protein glue of the myelin major dense line: Linking AlphaFold2 predictions to experimental data

Numerous human proteins are classified as intrinsically disordered proteins (IDPs). Due to their physicochemical properties, high-resolution structural information about IDPs is generally lacking. On the other hand, IDPs are known to adopt local ordered structures upon interactions with e.g. other proteins or lipid membrane surfaces. While recent developments in protein structure prediction have been revolutionary, their impact on IDP research at high resolution remains limited. We took a specific example of two myelin-specific IDPs, the myelin basic protein (MBP) and the cytoplasmic domain of myelin protein zero (P0ct). Both of these IDPs are crucial for normal nervous system development and function, and while they are disordered in solution, upon membrane binding, they partially fold into helices, being embedded into the lipid membrane. We carried out AlphaFold2 predictions of both proteins and analysed the models in light of experimental data related to protein structure and molecular interactions. We observe that the predicted models have helical segments that closely correspond to the membrane-binding sites on both proteins. We furthermore analyse the fits of the models to synchrotron-based X-ray scattering and circular dichroism data from the same IDPs. The models are likely to represent the membrane-bound state of both MBP and P0ct, rather than the conformation in solution. Artificial intelligence-based models of IDPs appear to provide information on the ligand-bound state of these proteins, instead of the conformers dominating free in solution. We further discuss the implications of the predictions for mammalian nervous system myelination and their relevance to understanding disease aspects of these IDPs.

Where protein structure and cell diversity meet

Protein-protein interaction networks - interactomes - are charted with the hope to understand how phenotypes emerge and how they are altered in disease states. Early efforts to map interactomes have focused on the assembly of context agnostic, reference networks. However, recent studies have mapped interactomes across different cell lines and tissues, finding highly variable interactomes due to the rewiring of protein-protein interactions in different contexts. Increasing evidence points to significant links between protein structure and interactome diversity seen across cell types and tissues. We discuss how recent findings support the key role of alternative splicing and phosphorylation, two well-established regulators of protein structural and functional diversity, in defining cell type- and tissue-specific interactomes. Moreover, we show that intrinsically disordered protein regions are most favorably equipped to support interactome rewiring by acting as hubs of protein structure and function regulation.

Intrinsically Disordered Proteins: An Overview

Many proteins and protein segments cannot attain a single stable three-dimensional structure under physiological conditions; instead, they adopt multiple interconverting conformational states. Such intrinsically disordered proteins or protein segments are highly abundant across proteomes, and are involved in various effector functions. This review focuses on different aspects of disordered proteins and disordered protein regions, which form the basis of the so-called "Disorder-function paradigm" of proteins. Additionally, various experimental approaches and computational tools used for characterizing disordered regions in proteins are discussed. Finally, the role of disordered proteins in diseases and their utility as potential drug targets are explored.

Does Generic Cyclic Kinase Insert Domain of Receptor Tyrosine Kinase KIT Clone Its Native Homologue?

Receptor tyrosine kinases (RTKs) are modular membrane proteins possessing both well-folded and disordered domains acting together in ligand-induced activation and regulation of post-transduction processes that tightly couple extracellular and cytoplasmic events. They ensure the fine-turning control of signal transmission by signal transduction. Deregulation of RTK KIT, including overexpression and gain of function mutations, has been detected in several human cancers. In this paper, we analysed by in silico techniques the Kinase Insert Domain (KID), a key platform of KIT transduction processes, as a generic macrocycle (KID^GC), a cleaved isolated polypeptide (KID^C), and a natively fused TKD domain (KID^D). We assumed that these KID species have similar structural and dynamic characteristics indicating the intrinsically disordered nature of this domain. This finding means that both polypeptides, cyclic KID^GC and linear KID^C, are valid models of KID integrated into the RTK KIT and will be helpful for further computational and empirical studies of post-transduction KIT events.

Extended disordered regions of ribosome-associated NAC proteins paralogs belong only to the germline in Drosophila melanogaster

The nascent polypeptide-associated complex (NAC) consisting of α- and β-subunits is an essential ribosome-associated protein conserved in eukaryotes. NAC is a ubiquitously expressed co-translational regulator of nascent protein folding and sorting providing for homeostasis of cellular proteins. Here we report on discovering the germline-specific NACαβ paralogs (gNACs), whose β-subunits, non-distinguishable by ordinary immunodetection, are encoded by five highly homologous gene copies, while the α-subunit is encoded by a single αNAC gene. The gNAC expression is detected in the primordial embryonic and adult gonads via immunostaining. The germline-specific α and β subunits differ from the ubiquitously expressed paralogs by the extended intrinsically disordered regions (IDRs) acquired at the N- and C-termini of the coding regions, predicted to be phosphorylated. The presence of distinct phosphorylated isoforms of gNAC-β subunits is confirmed by comparing of their profiles by 2D-isoeletrofocusing resolution before and after phosphatase treatment of testis ribosomes. We revealed that the predicted S/T sites of phosphorylation in the individual orthologous IDRs of gNAC-β sequences of Drosophila species are positionally conserved despite these disordered regions are drastically different. We propose the IDR-dependent molecular crowding and specific coordination of NAC and other proteostasis regulatory factors at the ribosomes of germinal cells. Our findings imply that there may be a functional crosstalk between the germinal and ubiquitous α- and β-subunits based on assessing their depletion effects on the fly viability and gonad development.

Structural dynamics shape the fitness window of alanine:glyoxylate aminotransferase

The conformational landscape of a protein is constantly expanded by genetic variations that have a minimal impact on the function(s) while causing subtle effects on protein structure. The wider the conformational space sampled by these variants, the higher the probabilities to adapt to changes in environmental conditions. However, the probability that a single mutation may result in a pathogenic phenotype also increases. Here we present a paradigmatic example of how protein evolution balances structural stability and dynamics to maximize protein adaptability and preserve protein fitness. We took advantage of known genetic variations of human alanine:glyoxylate aminotransferase (AGT1), which is present as a common major allelic form (AGT‐Ma) and a minor polymorphic form (AGT‐Mi) expressed in 20% of Caucasian population. By integrating crystallographic studies and molecular dynamics simulations, we show that AGT‐Ma is endowed with structurally unstable (frustrated) regions, which become disordered in AGT‐Mi. An in‐depth biochemical characterization of variants from an anticonsensus library, encompassing the frustrated regions, correlates this plasticity to a fitness window defined by AGT‐Ma and AGT‐Mi. Finally, co‐immunoprecipitation analysis suggests that structural frustration in AGT1 could favor additional functions related to protein–protein interactions. These results expand our understanding of protein structural evolution by establishing that naturally occurring genetic variations tip the balance between stability and frustration to maximize the ensemble of conformations falling within a well‐defined fitness window, thus expanding the adaptability potential of the protein.

Matrin3: Disorder and ALS Pathogenesis

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disorder characterized by the degeneration of both upper and lower motor neurons in the brain and spinal cord. ALS is associated with protein misfolding and inclusion formation involving RNA-binding proteins, including TAR DNA-binding protein (TDP-43) and fused in sarcoma (FUS). The 125-kDa Matrin3 is a highly conserved nuclear DNA/RNA-binding protein that is implicated in many cellular processes, including binding and stabilizing mRNA, regulating mRNA nuclear export, modulating alternative splicing, and managing chromosomal distribution. Mutations in MATR3, the gene encoding Matrin3, have been identified as causal in familial ALS (fALS). Matrin3 lacks a prion-like domain that characterizes many other ALS-associated RNA-binding proteins, including TDP-43 and FUS, however, our bioinformatics analyses and preliminary studies document that Matrin3 contains long intrinsically disordered regions that may facilitate promiscuous interactions with many proteins and may contribute to its misfolding. In addition, these disordered regions in Matrin3 undergo numerous post-translational modifications, including phosphorylation, ubiquitination and acetylation that modulate the function and misfolding of the protein. Here we discuss the disordered nature of Matrin3 and review the factors that may promote its misfolding and aggregation, two elements that might explain its role in ALS pathogenesis.

On the potential of machine learning to examine the relationship between sequence, structure, dynamics and function of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) constitute a broad set of proteins with few uniting and many diverging properties. IDPs-and intrinsically disordered regions (IDRs) interspersed between folded domains-are generally characterized as having no persistent tertiary structure; instead they interconvert between a large number of different and often expanded structures. IDPs and IDRs are involved in an enormously wide range of biological functions and reveal novel mechanisms of interactions, and while they defy the common structure-function paradigm of folded proteins, their structural preferences and dynamics are important for their function. We here discuss open questions in the field of IDPs and IDRs, focusing on areas where machine learning and other computational methods play a role. We discuss computational methods aimed to predict transiently formed local and long-range structure, including methods for integrative structural biology. We discuss the many different ways in which IDPs and IDRs can bind to other molecules, both via short linear motifs, as well as in the formation of larger dynamic complexes such as biomolecular condensates. We discuss how experiments are providing insight into such complexes and may enable more accurate predictions. Finally, we discuss the role of IDPs in disease and how new methods are needed to interpret the mechanistic effects of genomic variants in IDPs.

Folding and Intrinsic Disorder of the Receptor Tyrosine Kinase KIT Insert Domain Seen by Conventional Molecular Dynamics Simulations

The kinase insert domain (KID) of RTK KIT is the key recruitment region for downstream signalling proteins. KID, studied by molecular dynamics simulations as a cleaved polypeptide and as a native domain fused to KIT, showed intrinsic disorder represented by a set of heterogeneous conformations. The accurate atomistic models showed that the helical fold of KID is mainly sequence dependent. However, the reduced fold of the native KID suggests that its folding is allosterically controlled by the kinase domain. The tertiary structure of KID represents a compact array of highly variable α- and 3₁₀-helices linked by flexible loops playing a principal role in the conformational diversity. The helically folded KID retains a collapsed globule-like shape due to non-covalent interactions associated in a ternary hydrophobic core. The free energy landscapes constructed from first principles-the size, the measure of the average distance between the conformations, the amount of helices and the solvent-accessible surface area-describe the KID disorder through a collection of minima (wells), providing a direct evaluation of conformational ensembles. We found that the cleaved KID simulated with restricted N- and C-ends better reproduces the native KID than the isolated polypeptide. We suggest that a cyclic, generic KID would be best suited for future studies of KID f post-transduction effects.

The Participation of the Intrinsically Disordered Regions of the bHLH-PAS Transcription Factors in Disease Development

The basic helix–loop–helix/Per-ARNT-SIM (bHLH-PAS) proteins are a family of transcription factors regulating expression of a wide range of genes involved in different functions, ranging from differentiation and development control by oxygen and toxins sensing to circadian clock setting. In addition to the well-preserved DNA-binding bHLH and PAS domains, bHLH-PAS proteins contain long intrinsically disordered C-terminal regions, responsible for regulation of their activity. Our aim was to analyze the potential connection between disordered regions of the bHLH-PAS transcription factors, post-transcriptional modifications and liquid-liquid phase separation, in the context of disease-associated missense mutations. Highly flexible disordered regions, enriched in short motives which are more ordered, are responsible for a wide spectrum of interactions with transcriptional co-regulators. Based on our in silico analysis and taking into account the fact that the functions of transcription factors can be modulated by posttranslational modifications and spontaneous phase separation, we assume that the locations of missense mutations inducing disease states are clearly related to sequences directly undergoing these processes or to sequences responsible for their regulation.

Identification of the Primary Factors Determining theSpecificity of Human VKORC1 Recognition by Thioredoxin-Fold Proteins

Redox (reduction-oxidation) reactions control many important biological processes in all organisms, both prokaryotes and eukaryotes. This reaction is usually accomplished by canonical disulphide-based pathways involving a donor enzyme that reduces the oxidised cysteine residues of a target protein, resulting in the cleavage of its disulphide bonds. Focusing on human vitamin K epoxide reductase (hVKORC1) as a target and on four redoxins (protein disulphide isomerase (PDI), endoplasmic reticulum oxidoreductase (ERp18), thioredoxin-related transmembrane protein 1 (Tmx1) and thioredoxin-related transmembrane protein 4 (Tmx4)) as the most probable reducers of VKORC1, a comparative in-silico analysis that concentrates on the similarity and divergence of redoxins in their sequence, secondary and tertiary structure, dynamics, intraprotein interactions and composition of the surface exposed to the target is provided. Similarly, hVKORC1 is analysed in its native state, where two pairs of cysteine residues are covalently linked, forming two disulphide bridges, as a target for Trx-fold proteins. Such analysis is used to derive the putative recognition/binding sites on each isolated protein, and PDI is suggested as the most probable hVKORC1 partner. By probing the alternative orientation of PDI with respect to hVKORC1, the functionally related noncovalent complex formed by hVKORC1 and PDI was found, which is proposed to be a first precursor to probe thiol-disulphide exchange reactions between PDI and hVKORC1.