Search

The Wnt/β-catenin signaling pathway plays numerous essential roles in animal development and tissue/stem cell maintenance. The activation of genes regulated by Wnt/β-catenin signaling requires the nuclear accumulation of β-catenin, a transcriptional co-activator. β-catenin is recruited to many Wnt-regulated enhancers through direct binding to T-cell factor/lymphoid enhancer factor (TCF/LEF) family transcription factors. β-catenin has previously been reported to form phase-separated biomolecular condensates (BMCs), which was implicated as a component of β-catenin's mechanism of action. This function required aromatic amino acid residues in the intrinsically disordered regions (IDRs) at the N- and C-termini of the protein. In this report, we further explore a role for β-catenin BMCs in Wnt target gene regulation. We find that β-catenin BMCs are miscible with LEF1 BMCs in vitro and in cultured cells. We characterized a panel of β-catenin mutants with different combinations of aromatic residue mutations in human cell culture and Drosophila melanogaster. Our data support a model in which aromatic residues across both IDRs contribute to BMC formation and signaling activity. Although different Wnt targets have different sensitivities to loss of β-catenin's aromatic residues, the activation of every target examined was compromised by aromatic substitution. These mutants are not defective in nuclear import or co-immunoprecipitation with several β-catenin binding partners. In addition, residues in the N-terminal IDR with no previously known role in signaling are clearly required for the activation of various Wnt readouts. Consistent with this, deletion of the N-terminal IDR results in a loss of signaling activity, which can be rescued by the addition of heterologous IDRs enriched in aromatic residues. Overall, our work supports a model in which the ability of β-catenin to form biomolecular condensates in the nucleus is tightly linked to its function as a transcriptional co-regulator., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Stewart et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Ataxin-2 (ATXN2) is a gene implicated in spinocerebellar ataxia type II (SCA2), amyotrophic lateral sclerosis (ALS) and Parkinsonism. The encoded protein is a therapeutic target for ALS and related conditions. ATXN2 (or Atx2 in insects) can function in translational activation, translational repression, mRNA stability and in the assembly of mRNP-granules, a process mediated by intrinsically disordered regions (IDRs). Previous work has shown that the LSm (Like-Sm) domain of Atx2, which can help stimulate mRNA translation, antagonizes mRNP-granule assembly. Here we advance these findings through a series of experiments on Drosophila and human Ataxin-2 proteins. Results of Targets of RNA Binding Proteins Identified by Editing (TRIBE), co-localization and immunoprecipitation experiments indicate that a polyA-binding protein (PABP) interacting, PAM2 motif of Ataxin-2 may be a major determinant of the mRNA and protein content of Ataxin-2 mRNP granules. Experiments with transgenic Drosophila indicate that while the Atx2-LSm domain may protect against neurodegeneration, structured PAM2- and unstructured IDR- interactions both support Atx2-induced cytotoxicity. Taken together, the data lead to a proposal for how Ataxin-2 interactions are remodelled during translational control and how structured and non-structured interactions contribute differently to the specificity and efficiency of RNP granule condensation as well as to neurodegeneration., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Petrauskas et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Intrinsically disordered regions (IDRs) are segments of proteins without stable three-dimensional structures. As this flexibility allows them to interact with diverse binding partners, IDRs play key roles in cell signaling and gene expression. Despite the prevalence and importance of IDRs in eukaryotic proteomes and various biological processes, associating them with specific molecular functions remains a significant challenge due to their high rates of sequence evolution. However, by comparing the observed values of various IDR-associated properties against those generated under a simulated model of evolution, a recent study found most IDRs across the entire yeast proteome contain conserved features. Furthermore, it showed clusters of IDRs with common "evolutionary signatures," i.e. patterns of conserved features, were associated with specific biological functions. To determine if similar patterns of conservation are found in the IDRs of other systems, in this work we applied a series of phylogenetic models to over 7,500 orthologous IDRs identified in the Drosophila genome to dissect the forces driving their evolution. By comparing models of constrained and unconstrained continuous trait evolution using the Brownian motion and Ornstein-Uhlenbeck models, respectively, we identified signals of widespread constraint, indicating conservation of distributed features is mechanism of IDR evolution common to multiple biological systems. In contrast to the previous study in yeast, however, we observed limited evidence of IDR clusters with specific biological functions, which suggests a more complex relationship between evolutionary constraints and function in the IDRs of multicellular organisms., Competing Interests: I have read the journal’s policy and the authors of this manuscript have the following competing interests: MBE is a founder and former member of the board of directors of PLOS., (Copyright: © 2024 Singleton, Eisen. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

One central question for cell and developmental biologists is defining how epithelial cells can change shape and move during embryonic development without tearing tissues apart. This requires robust yet dynamic connections of cells to one another, via the cell-cell adherens junction, and of junctions to the actin and myosin cytoskeleton, which generates force. The last decade revealed that these connections involve a multivalent network of proteins, rather than a simple linear pathway. We focus on Drosophila Canoe, homolog of mammalian Afadin, as a model for defining the underlying mechanisms. Canoe and Afadin are complex, multidomain proteins that share multiple domains with defined and undefined binding partners. Both also share a long carboxy-terminal intrinsically disordered region (IDR), whose function is less well defined. IDRs are found in many proteins assembled into large multiprotein complexes. We have combined bioinformatic analysis and the use of a series of canoe mutants with early stop codons to explore the evolution and function of the IDR. Our bioinformatic analysis reveals that the IDRs of Canoe and Afadin differ dramatically in sequence and sequence properties. When we looked over shorter evolutionary time scales, we identified multiple conserved motifs. Some of these are predicted by AlphaFold to be alpha-helical, and two correspond to known protein interaction sites for alpha-catenin and F-actin. We next identified the lesions in a series of eighteen canoe mutants, which have early stop codons across the entire protein coding sequence. Analysis of their phenotypes are consistent with the idea that the IDR, including the conserved motifs in the IDR, are critical for protein function. These data provide the foundation for further analysis of IDR function., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Gurley et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Intrinsically disordered proteins (IDPs) are proteins that lack a stable 3D structure but maintain a biological function. It has been frequently suggested that IDPs are difficult to align because they tend to have fewer conserved residues compared to ordered proteins, but to our knowledge this has never been directly tested. To compare the alignments of ordered proteins to IDPs, their multiple sequence alignments (MSAs) were assessed using two different methods. The first compared the similarity between MSAs produced using the same sequences but created with Clustal Omega, MAFFT, and MUSCLE. The second assessed MSAs based on how well they recapitulated the species tree. These two methods measure the "correctness" of an MSA with two different approaches; the first method measures consistency while the second measures the underlying phylogenetic signal. Proteins that contained both regions of disorder and order were analyzed along with proteins that were fully disordered and fully ordered, using nucleotide, codon and peptide sequence alignments. We observed that IDPs had less similar MSAs than ordered proteins, which is most likely linked to the lower sequence conservation in IDPs. However, comparisons of tree distances found that trees from the ordered sequence MSAs were not significantly closer to the species tree than those inferred from disordered sequence MSAs. Our results show that it is correct to say that IDPs are difficult to align on the basis of MSA consistency, but that this does not equate with alignments being of poor quality when assessed by their ability to correctly infer a species tree., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Riley et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Vertebrate kidneys contribute to homeostasis by regulating electrolyte, acid-base balance, removing toxic metabolites from blood, and preventing protein loss into the urine. Glomerular podocytes constitute the blood-urine barrier, and podocyte slit-diaphragm (SD), a modified tight junction, contributes to the glomerular permselectivity. Nephrin, KIRREL1, podocin, CD2AP, and TRPC6 are crucial members of the SD that interact with each other and contribute to the SD's structural and functional integrity. This study analyzed the distribution of these five essential SD proteins across the organisms for which the genome sequence is available. We found a diverse distribution of nephrin and KIRREL1 ranging from nematodes to higher vertebrates, whereas podocin, CD2AP, and TRPC6 are restricted to the vertebrates. Among invertebrates, nephrin and its orthologs consist of more immunoglobulin-3 domains, whereas in the vertebrates, CD80-like C2-set domains are predominant. In the case of KIRREL1 and its orthologs, more Ig domains were observed in invertebrates than vertebrates. Src Homology-3 (SH3) domain of CD2AP and SPFH domain of podocin are highly conserved among vertebrates. TRPC6 and its orthologs had conserved ankyrin repeats, TRP, and ion transport domains, except Chondrichthyes and Echinodermata, which do not possess the ankyrin repeats. Intrinsically unstructured regions (IURs) are conserved across the SD orthologs, suggesting IURs importance in the protein complexes that constitute the slit-diaphragm. For the first time, a study reports the evolutionary insights of vertebrate SD proteins and their invertebrate orthologs., Competing Interests: The authors have declared that no competing interests exist.

Protein-protein interactions are involved in a wide range of cellular processes. These interactions often involve intrinsically disordered proteins (IDPs) and protein binding domains. However, the details of IDP binding pathways are hard to characterize using experimental approaches, which can rarely capture intermediate states present at low populations. SH3 domains are common protein interaction domains that typically bind proline-rich disordered segments and are involved in cell signaling, regulation, and assembly. We hypothesized, given the flexibility of SH3 binding peptides, that their binding pathways include multiple steps important for function. Molecular dynamics simulations were used to characterize the steps of binding between the yeast Abp1p SH3 domain (AbpSH3) and a proline-rich IDP, ArkA. Before binding, the N-terminal segment 1 of ArkA is pre-structured and adopts a polyproline II helix, while segment 2 of ArkA (C-terminal) adopts a 310 helix, but is far less structured than segment 1. As segment 2 interacts with AbpSH3, it becomes more structured, but retains flexibility even in the fully engaged state. Binding simulations reveal that ArkA enters a flexible encounter complex before forming the fully engaged bound complex. In the encounter complex, transient nonspecific hydrophobic and long-range electrostatic contacts form between ArkA and the binding surface of SH3. The encounter complex ensemble includes conformations with segment 1 in both the forward and reverse orientation, suggesting that segment 2 may play a role in stabilizing the correct binding orientation. While the encounter complex forms quickly, the slow step of binding is the transition from the disordered encounter ensemble to the fully engaged state. In this transition, ArkA makes specific contacts with AbpSH3 and buries more hydrophobic surface. Simulating the binding between ApbSH3 and ArkA provides insight into the role of encounter complex intermediates and nonnative hydrophobic interactions for other SH3 domains and IDPs in general., Competing Interests: The authors have declared that no competing interests exist.

The role of electrostatic interactions and mutations that change charge states in intrinsically disordered proteins (IDPs) is well-established, but many disease-associated mutations in IDPs are charge-neutral. The Val66Met single nucleotide polymorphism (SNP) in precursor brain-derived neurotrophic factor (BDNF) is one of the earliest SNPs to be associated with neuropsychiatric disorders, and the underlying molecular mechanism is unknown. Here we report on over 250 μs of fully-atomistic, explicit solvent, temperature replica-exchange molecular dynamics (MD) simulations of the 91 residue BDNF prodomain, for both the V66 and M66 sequence. The simulations were able to correctly reproduce the location of both local and non-local secondary structure changes due to the Val66Met mutation, when compared with NMR spectroscopy. We find that the change in local structure is mediated via entropic and sequence specific effects. We developed a hierarchical sequence-based framework for analysis and conceptualization, which first identifies "blobs" of 4-15 residues representing local globular regions or linkers. We use this framework within a novel test for enrichment of higher-order (tertiary) structure in disordered proteins; the size and shape of each blob is extracted from MD simulation of the real protein (RP), and used to parameterize a self-avoiding heterogenous polymer (SAHP). The SAHP version of the BDNF prodomain suggested a protein segmented into three regions, with a central long, highly disordered polyampholyte linker separating two globular regions. This effective segmentation was also observed in full simulations of the RP, but the Val66Met substitution significantly increased interactions across the linker, as well as the number of participating residues. The Val66Met substitution replaces β-bridging between V66 and V94 (on either side of the linker) with specific side-chain interactions between M66 and M95. The protein backbone in the vicinity of M95 is then free to form β-bridges with residues 31-41 near the N-terminus, which condenses the protein. A significant role for Met/Met interactions is consistent with previously-observed non-local effects of the Val66Met SNP, as well as established interactions between the Met66 sequence and a Met-rich receptor that initiates neuronal growth cone retraction., Competing Interests: The authors have declared that no competing interests exist.

Many studies about classification and the functional annotation of intrinsically disordered proteins (IDPs) are based on either the occurrence of long disordered regions or the fraction of disordered residues in the sequence. Taking into account both criteria we separate the human proteome, taken as a case study, into three variants of proteins: i) ordered proteins (ORDPs), ii) structured proteins with intrinsically disordered regions (IDPRs), and iii) intrinsically disordered proteins (IDPs). The focus of this work is on the different functional roles of IDPs and IDPRs, which up until now have been generally considered as a whole. Previous studies assigned a large set of functional roles to the general category of IDPs. We show here that IDPs and IDPRs have non-overlapping functional spectra, play different roles in human diseases, and deserve to be treated as distinct categories of proteins. IDPs enrich only a few classes, functions, and processes: nucleic acid binding proteins, chromatin binding proteins, transcription factors, and developmental processes. In contrast, IDPRs are spread over several functional protein classes and GO annotations which they partly share with ORDPs. As regards to diseases, we observe that IDPs enrich only cancer-related proteins, at variance with previous results reporting that IDPs are widespread also in cardiovascular and neurodegenerative pathologies. Overall, the operational separation of IDPRs from IDPs is relevant towards correct estimates of the occurrence of intrinsically disordered proteins in genome-wide studies and in the understanding of the functional spectra associated to different flavors of protein disorder., Competing Interests: The authors have declared that no competing interests exist.

Intrinsic disorder is more abundant in eukaryotic than prokaryotic proteins. Methods predicting intrinsic disorder are based on the amino acid sequence of a protein. Therefore, there must exist an underlying difference in the sequences between eukaryotic and prokaryotic proteins causing the (predicted) difference in intrinsic disorder. By comparing proteins, from complete eukaryotic and prokaryotic proteomes, we show that the difference in intrinsic disorder emerges from the linker regions connecting Pfam domains. Eukaryotic proteins have more extended linker regions, and in addition, the eukaryotic linkers are significantly more disordered, 38% vs. 12-16% disordered residues. Next, we examined the underlying reason for the increase in disorder in eukaryotic linkers, and we found that the changes in abundance of only three amino acids cause the increase. Eukaryotic proteins contain 8.6% serine; while prokaryotic proteins have 6.5%, eukaryotic proteins also contain 5.4% proline and 5.3% isoleucine compared with 4.0% proline and ≈ 7.5% isoleucine in the prokaryotes. All these three differences contribute to the increased disorder in eukaryotic proteins. It is tempting to speculate that the increase in serine frequencies in eukaryotes is related to regulation by kinases, but direct evidence for this is lacking. The differences are observed in all phyla, protein families, structural regions and type of protein but are most pronounced in disordered and linker regions. The observation that differences in the abundance of three amino acids cause the difference in disorder between eukaryotic and prokaryotic proteins raises the question: Are amino acid frequencies different in eukaryotic linkers because the linkers are more disordered or do the differences cause the increased disorder?, Competing Interests: The authors have declared that no competing interests exist.

Conformational intrinsic disorder is a feature present in many virus proteins. Intrinsically disordered regions (IDRs) have weaker structural requirement than ordered regions and mutations in IDRs could have a lower impact on the virus fitness. This could favor its exploration of adaptive solutions. The potyviral protein VPg contains IDRs with determinants for adaptation to its host plant. To experimentally assess whether IDRs are more resistant to mutations than ordered regions, the biologically relevant interaction between mutant libraries of both VPg and the eukaryotic translation initiation factor 4E (eIF4E) and their respective wild type partner was examined using yeast two hybrid assay. Our data shows that VPg is significantly more robust to mutations than eIF4E and as such belongs to a particular class of intrinsically disordered proteins. This result is discussed from the standpoint of IDRs involvement in the virus adaptive processes., Competing Interests: The authors have declared that no competing interests exist.

The proteins Oskar (Osk) in Drosophila and Bucky ball (Buc) in zebrafish act as germ plasm organizers. Both proteins recapitulate germ plasm activities but seem to be unique to their animal groups. Here, we discover that Osk and Buc show similar activities during germ cell specification. Drosophila Osk induces additional PGCs in zebrafish. Surprisingly, Osk and Buc do not show homologous protein motifs that would explain their related function. Nonetheless, we detect that both proteins contain stretches of intrinsically disordered regions (IDRs), which seem to be involved in protein aggregation. IDRs are known to rapidly change their sequence during evolution, which might obscure biochemical interaction motifs. Indeed, we show that Buc binds to the known Oskar interactors Vasa protein and nanos mRNA indicating conserved biochemical activities. These data provide a molecular framework for two proteins with unrelated sequence but with equivalent function to assemble a conserved core-complex nucleating germ plasm., Competing Interests: The authors have declared that no competing interests exist.

TMEM106B was initially identified as a risk factor for FTLD, but recent studies highlighted its general role in neurodegenerative diseases. Very recently TMEM106B has also been characterized to regulate aging phenotypes. TMEM106B is a 274-residue lysosomal protein whose cytoplasmic domain functions in the endosomal/autophagy pathway by dynamically and transiently interacting with diverse categories of proteins but the underlying structural basis remains completely unknown. Here we conducted bioinformatics analysis and biophysical characterization by CD and NMR spectroscopy, and obtained results reveal that the TMEM106B cytoplasmic domain is intrinsically disordered with no well-defined three-dimensional structure. Nevertheless, detailed analysis of various multi-dimensional NMR spectra allowed defining residue-specific conformations and dynamics. Overall, the TMEM106B cytoplasmic domain is lacking of any tight tertiary packing and relatively flexible. However, several segments are populated with dynamic/nascent secondary structures and have relatively restricted backbone motions on ps-ns time scale, as indicated by their positive {1H}-15N steady-state NOE. Our study thus decodes that being intrinsically disordered may allow the TMEM106B cytoplasmic domain to dynamically and transiently interact with a variety of distinct partners., Competing Interests: The authors have declared that no competing interests exist.

In C. elegans, efficient RNA silencing requires small RNA amplification mediated by RNA-dependent RNA polymerases (RdRPs). RRF-1, an RdRP, and other Mutator complex proteins localize to Mutator foci, which are perinuclear germline foci that associate with nuclear pores and P granules to facilitate small RNA amplification. The Mutator complex protein MUT-16 is critical for Mutator foci assembly. By analyzing small deletions of MUT-16, we identify specific regions of the protein that recruit other Mutator complex components and demonstrate that it acts as a scaffolding protein. We further determine that the C-terminal region of MUT-16, a portion of which contains predicted intrinsic disorder, is necessary and sufficient to promote Mutator foci formation. Finally, we establish that MUT-16 foci have many properties consistent with a phase-separated condensate and propose that Mutator foci form through liquid-liquid phase separation of MUT-16. P granules, which contain additional RNA silencing proteins, have previously been shown to have liquid-like properties. Thus, RNA silencing in C. elegans germ cells may rely on multiple phase-separated compartments through which sorting, processing, and silencing of mRNAs occurs., Competing Interests: The authors have declared that no competing interests exist.

De novo creation of protein coding genes involves the formation of short ORFs from noncoding regions; some of these ORFs might then become fixed in the population. These orphan proteins need to, at the bare minimum, not cause serious harm to the organism, meaning that they should for instance not aggregate. Therefore, although the creation of short ORFs could be truly random, the fixation should be subjected to some selective pressure. The selective forces acting on orphan proteins have been elusive, and contradictory results have been reported. In Drosophila young proteins are more disordered than ancient ones, while the opposite trend is present in yeast. To the best of our knowledge no valid explanation for this difference has been proposed. To solve this riddle we studied structural properties and age of proteins in 187 eukaryotic organisms. We find that, with the exception of length, there are only small differences in the properties between proteins of different ages. However, when we take the GC content into account we noted that it could explain the opposite trends observed for orphans in yeast (low GC) and Drosophila (high GC). GC content is correlated with codons coding for disorder promoting amino acids. This leads us to propose that intrinsic disorder is not a strong determining factor for fixation of orphan proteins. Instead these proteins largely resemble random proteins given a particular GC level. During evolution the properties of a protein change faster than the GC level causing the relationship between disorder and GC to gradually weaken.

Yersinia pestis Caf1 is a multifunctional protein responsible for antiphagocytic activity and is a key protective antigen. It is generally conserved between globally distributed Y. pestis strains, but Y. pestis subsp. microtus biovar caucasica strains circulating within populations of common voles in Georgia and Armenia were reported to carry a single substitution of alanine to serine. We investigated polymorphism of the Caf1 sequences among other Y. pestis subsp. microtus strains, which have a limited virulence in guinea pigs and in humans. Sequencing of caf1 genes from 119 Y. pestis strains belonging to different biovars within subsp. microtus showed that the Caf1 proteins exist in three isoforms, the global type Caf1NT1 (Ala48 Phe117), type Caf1NT2 (Ser48 Phe117) found in Transcaucasian-highland and Pre-Araks natural plague foci #4-7, and a novel Caf1NT3 type (Ala48 Val117) endemic in Dagestan-highland natural plague focus #39. Both minor types are the progenies of the global isoform. In this report, Caf1 polymorphism was analyzed by comparing predicted intrinsic disorder propensities and potential protein-protein interactivities of the three Caf1 isoforms. The analysis revealed that these properties of Caf1 protein are minimally affected by its polymorphism. All protein isoforms could be equally detected by an immunochromatography test for plague at the lowest protein concentration tested (1.0 ng/mL), which is the detection limit. When compared to the classic Caf1NT1 isoform, the endemic Caf1NT2 or Caf1NT3 had lower immunoreactivity in ELISA and lower indices of self- and cross-protection. Despite a visible reduction in cross-protection between all Caf1 isoforms, our data suggest that polymorphism in the caf1 gene may not allow the carriers of Caf1NT2 or Caf1NT3 variants escaping from the Caf1NT1-mediated immunity to plague in the case of a low-dose flea-borne infection., Competing Interests: The authors have declared that no competing interests exist.

Background: Prothymosin α (ProTα) (isoform 2: iso2) is a widely distributed, small acidic protein with intracellular and extracellular-associated functions. Recently, we identified two new ProTα variants with potent anti-HIV activity from CD8+ T cells and cervicovaginal lavage. The first is a splice variant of the ProTα gene known as isoB and the second is the product of ProTα pseudogene 7 (p7). Similarly to iso2, the anti-HIV activity of both variants is mediated by type I IFN. Here we tested whether the immunomodulatory activity of isoB and p7 are also TLR4 dependent and determined their kinetic of release in response to HIV-1 infection., Methods: Type I, type III, TNF-α and IL-6 mRNA inducing activity was determined in macrophages from wild type and TLR4 knockout mice treated with recombinant ProTα variants. Supernatants from mock and HIV infected cells were analyzed by mass spectrometry in positive and negative modes for the presence of ProTα variants. In silico structural and functional analysis of ProTα variants were performed., Results: We show that both isoB and p7 upregulate IFN-β, IFN-λ1, IL-6, TNF-α and RANTES mRNAs in primary human macrophages. The potent stimulation of IFN-β by the recombinant ProTα variants in human macrophages is dependent on the TLR4 pathway, whereas the induction of TNF-α and IL-6 may also occur independently of TLR4, suggesting the interaction of ProTα variants with other signaling molecules/receptors. In silico analyses confirmed that the novel isoB and p7 variants are intrinsically disordered proteins, which lack the NLS and mass spectrometry showed release of ProTα variants within minutes post HIV-1 infection. These features are consistent with the function of ProTα variants as damage associate molecular patterns (DAMPs)., Conclusions: Our findings indicate that ProTα variants strongly inhibit viral replication mainly, but not exclusively, through TLR4 signaling and that they are released within minutes of viral infection suggesting that they may function as DAMPs.

Disordered regions and Intrinsically Disordered Proteins (IDPs) are involved in critical cellular processes and may acquire a stable three-dimensional structure only upon binding to their partners. IDPs may follow a folding-after-binding process, known as induced folding, or a folding-before-binding process, known as conformational selection. The transcription factor p53 is involved in the regulation of cellular events that arise upon stress or DNA damage. The p53 domain structure is composed of an N-terminal transactivation domain (p53TAD), a DNA Binding Domain and a tetramerization domain. The activity of TAD is tightly regulated by interactions with cofactors, inhibitors and phosphorylation. To initiate transcription, p53TAD binds to the TAZ2 domain of CBP, a co-transcription factor, and undergoes a folding and binding process, as revealed by the recent NMR structure of the complex. The activity of p53 is regulated by phosphorylation at multiple sites on the TAD domain and recent studies have shown that modifications at three residues affect the binding towards TAZ2. However, we still do not know how these phosphorylations affect the structure of the bound state and, therefore, how they regulate the p53 function. In this work, we have used computational simulations to understand how phosphorylation affects the structure of the p53TAD:TAZ2 complex and regulates the recognition mechanism. Phosphorylation has been proposed to enhance binding by direct interaction with the folded protein or by changing the unbound conformation of IDPs, for example by pre-folding the protein favoring the recognition mechanism. Here, we show an interesting turn in the p53 case: phosphorylation mainly affects the bound structure of p53TAD, highlighting the complexity of IDP protein-protein interactions. Our results are in agreement with previous experimental studies, allowing a clear picture of how p53 is regulated by phosphorylation and giving new insights into how post-translational modifications can regulate the function of IDPs.

Malaria remains a significant global health burden. The development of an effective malaria vaccine remains as a major challenge with the potential to significantly reduce morbidity and mortality. While Plasmodium spp. have been shown to contain a large number of intrinsically disordered proteins (IDPs) or disordered protein regions, the relationship of protein structure to subcellular localisation and adaptive immune responses remains unclear. In this study, we employed several computational prediction algorithms to identify IDPs at the proteome level of six Plasmodium spp. and to investigate the potential impact of protein disorder on adaptive immunity against P. falciparum parasites. IDPs were shown to be particularly enriched within nuclear proteins, apical proteins, exported proteins and proteins localised to the parasitophorous vacuole. Furthermore, several leading vaccine candidates, and proteins with known roles in host-cell invasion, have extensive regions of disorder. Presentation of peptides by MHC molecules plays an important role in adaptive immune responses, and we show that IDP regions are predicted to contain relatively few MHC class I and II binding peptides owing to inherent differences in amino acid composition compared to structured domains. In contrast, linear B-cell epitopes were predicted to be enriched in IDPs. Tandem repeat regions and non-synonymous single nucleotide polymorphisms were found to be strongly associated with regions of disorder. In summary, immune responses against IDPs appear to have characteristics distinct from those against structured protein domains, with increased antibody recognition of linear epitopes but some constraints for MHC presentation and issues of polymorphisms. These findings have major implications for vaccine design, and understanding immunity to malaria.

Background: Intrinsically disordered proteins (IDPs) lack a stable tertiary structure in isolation. Remarkably, however, a substantial portion of IDPs undergo disorder-to-order transitions upon binding to their cognate partners. Structural flexibility and binding plasticity enable IDPs to interact with a broad range of partners. However, the broader network properties that could provide additional insights into the functional role of IDPs are not known., Results: Here, we report the first comprehensive survey of network properties of IDP-induced sub-networks in multiple species from yeast to human. Our results show that IDPs exhibit greater-than-expected modularity and are connected to the rest of the protein interaction network (PIN) via proteins that exhibit the highest betweenness centrality and connect to fewer-than-expected IDP communities, suggesting that they form critical communication links from IDP modules to the rest of the PIN. Moreover, we found that IDPs are enriched at the top level of regulatory hierarchy., Conclusion: Overall, our analyses reveal coherent and remarkably conserved IDP-centric network properties, namely, modularity in IDP-induced network and a layer of critical nodes connecting IDPs with the rest of the PIN.

Bdellovibrio bacteriovorus is a δ-proteobacterium that preys upon Salmonella spp., E. coli, and other Gram-negative bacteria. Bdellovibrio can grow axenically (host-independent, HI, rare and mutation-driven) or subsist via a predatory lifecycle (host-dependent, HD, the usual case). Upon contact with prey, B. bacteriovorus enters the host periplasm from where it slowly drains the host cytosol of nutrients for its own replication. At the core of this mechanism is a retractile pilus, whose architecture is regulated by the protein Bd0108 and its interaction with the neighboring gene product Bd0109. Deletion of bd0108 results in negligible pilus formation, whereas an internal deletion (the one that instigates host-independence) causes mis-regulation of pilus length. These mutations, along with a suite of naturally occurring bd0108 mutant strains, act to control the entry to HI growth. To further study the molecular mechanism of predatory regulation, we focused on the apparent lifecycle switch protein Bd0108. Here we characterize the solution structure and dynamics of Bd0108 using nuclear magnetic resonance (NMR) spectroscopy complemented with additional biophysical methods. We then explore the interaction between Bd0108 and Bd0109 in detail utilizing isothermal titration calorimetry (ITC) and NMR spectroscopy. Together our results demonstrate that Bd0108 is an intrinsically disordered protein (IDP) and that the interaction with Bd0109 is of low affinity. Furthermore, we observe that Bd0108 retains an IDP nature while binding Bd0109. From our data we conclude that Bdellovibrio bacteriovorus utilizes an intrinsically disordered protein to regulate its pilus and control predation signaling.

Fish otoliths, biominerals composed of calcium carbonate with a small amount of organic matrix, are involved in the functioning of the inner ear. Starmaker (Stm) from zebrafish (Danio rerio) was the first protein found to be capable of controlling the formation of otoliths. Recently, a gene was identified encoding the Starmaker-like (Stm-l) protein from medaka (Oryzias latipes), a putative homologue of Stm and human dentine sialophosphoprotein. Although there is no sequence similarity between Stm-l and Stm, Stm-l was suggested to be involved in the biomineralization of otoliths, as had been observed for Stm even before. The molecular properties and functioning of Stm-l as a putative regulatory protein in otolith formation have not been characterized yet. A comprehensive biochemical and biophysical analysis of recombinant Stm-l, along with in silico examinations, indicated that Stm-l exhibits properties of a coil-like intrinsically disordered protein. Stm-l possesses an elongated and pliable structure that is able to adopt a more ordered and rigid conformation under the influence of different factors. An in vitro assay of the biomineralization activity of Stm-l indicated that Stm-l affected the size, shape and number of calcium carbonate crystals. The functional significance of intrinsically disordered properties of Stm-l and the possible role of this protein in controlling the formation of calcium carbonate crystals is discussed.

Evolution is driven by mutations, which lead to new protein functions but come at a cost to protein stability. Non-conservative substitutions are of interest in this regard because they may most profoundly affect both function and stability. Accordingly, organisms must balance the benefit of accepting advantageous substitutions with the possible cost of deleterious effects on protein folding and stability. We here examine factors that systematically promote non-conservative mutations at the proteome level. Intrinsically disordered regions in proteins play pivotal roles in protein interactions, but many questions regarding their evolution remain unanswered. Similarly, whether and how molecular chaperones, which have been shown to buffer destabilizing mutations in individual proteins, generally provide robustness during proteome evolution remains unclear. To this end, we introduce an evolutionary parameter λ that directly estimates the rate of non-conservative substitutions. Our analysis of λ in Escherichia coli, Saccharomyces cerevisiae, and Homo sapiens sequences reveals how co- and post-translationally acting chaperones differentially promote non-conservative substitutions in their substrates, likely through buffering of their destabilizing effects. We further find that λ serves well to quantify the evolution of intrinsically disordered proteins even though the unstructured, thus generally variable regions in proteins are often flanked by very conserved sequences. Crucially, we show that both intrinsically disordered proteins and highly re-wired proteins in protein interaction networks, which have evolved new interactions and functions, exhibit a higher λ at the expense of enhanced chaperone assistance. Our findings thus highlight an intricate interplay of molecular chaperones and protein disorder in the evolvability of protein networks. Our results illuminate the role of chaperones in enabling protein evolution, and underline the importance of the cellular context and integrated approaches for understanding proteome evolution. We feel that the development of λ may be a valuable addition to the toolbox applied to understand the molecular basis of evolution.

Synonymous constraint elements (SCEs) are protein-coding genomic regions with very low synonymous mutation rates believed to carry additional, overlapping functions. Thousands of such potentially multi-functional elements were recently discovered by analyzing the levels and patterns of evolutionary conservation in human coding exons. These elements provide a good opportunity to improve our understanding of how the redundant nature of the genetic code is exploited in the cell. Our premise is that the protein segments encoded by such elements might better comply with the increased functional demands if they are structurally less constrained (i.e. intrinsically disordered). To test this idea, we investigated the protein segments encoded by SCEs with computational tools to describe the underlying structural properties. In addition to SCEs, we examined the level of disorder, secondary structure, and sequence complexity of protein regions overlapping with experimentally validated splice regulatory sites. We show that multi-functional gene regions translate into protein segments that are significantly enriched in structural disorder and compositional bias, while they are depleted in secondary structure and domain annotations compared to reference segments of similar lengths. This tendency suggests that relaxed protein structural constraints provide an advantage when accommodating multiple overlapping functions in coding regions.

With the exponential increase in the number of sequenced organisms, automated annotation of proteins is becoming increasingly important. Intrinsically disordered regions are known to play a significant role in protein function. Despite their abundance, especially in eukaryotes, they are rarely used to inform function prediction systems. In this study, we extracted seven sequence features in intrinsically disordered regions and developed a scheme to use them to predict Gene Ontology Slim terms associated with proteins. We evaluated the function prediction performance of each feature. Our results indicate that the residue composition based features have the highest precision while bigram probabilities, based on sequence profiles of intrinsically disordered regions obtained from PSIBlast, have the highest recall. Amino acid bigrams and features based on secondary structure show an intermediate level of precision and recall. Almost all features showed a high prediction performance for GO Slim terms related to extracellular matrix, nucleus, RNA and DNA binding. However, feature performance varied significantly for different GO Slim terms emphasizing the need for a unique classifier optimized for the prediction of each functional term. These findings provide a first comprehensive and quantitative evaluation of sequence features in intrinsically disordered regions and will help in the development of a more informative protein function predictor.

The transcription factor Gli3 is acting mainly as a transcriptional repressor in the Sonic hedgehog signal transduction pathway. Gli3 contains a repressor domain in its N-terminus from residue G106 to E236. In this study we have characterized the intracellular structure of the Gli3 repressor domain using a combined bioinformatics and experimental approach. According to our findings the Gli3 repressor domain while being intrinsically disordered contains predicted anchor sites for partner interactions. The obvious interaction partners to test were Ski and DNA; however, with both of these the structure of Gli3 repressor domain remained disordered. To locate residues important for the repressor function we mutated several residues within the Gli3 repressor domain. Two of these, H141A and H157N, targeting predicted helical regions, significantly decreased transcriptional repression and thus identify important functional parts of the domain.