Analysis of the relationships between evolvability, thermodynamics, and the functions of intrinsically disordered proteins/regions

Uncovering the structural stability of Magnaporthe oryzae effectors: a secretome-wide in silico analysis

Rice blast, caused by the ascomycete fungus Magnaporthe oryzae, is a deadly disease and a major threat to global food security. The pathogen secretes small proteinaceous effectors, virulence factors, inside the host to manipulate and perturb the host immune system, allowing the pathogen to colonize and establish a successful infection. While the molecular functions of several effectors are characterized, very little is known about the structural stability of these effectors. We analyzed a total of 554 small secretory proteins (SSPs) from the M. oryzae secretome to decipher key features of intrinsic disorder (ID) and the structural dynamics of the selected putative effectors through thorough and systematic in silico studies. Our results suggest that out of the total SSPs, 66% were predicted as effector proteins, released either into the apoplast or cytoplasm of the host cell. Of these, 68% were found to be intrinsically disordered effector proteins (IDEPs). Among the six distinct classes of disordered effectors, we observed peculiar relationships between the localization of several effectors in the apoplast or cytoplasm and the degree of disorder. We determined the degree of structural disorder and its impact on protein foldability across all the putative small secretory effector proteins from the blast pathogen, further validated by molecular dynamics simulation studies. This study provides definite clues toward unraveling the mystery behind the importance of structural distortions in effectors and their impact on plant-pathogen interactions. The study of these dynamical segments may help identify new effectors as well.Communicated by Ramaswamy H. Sarma.

The difficulty of aligning intrinsically disordered protein sequences as assessed by conservation and phylogeny

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.

Co-opting disorder into order: Intrinsically disordered proteins and the early evolution of complex multicellularity

Intrinsically disordered proteins (IDPs) are proteins that lack rigid structures yet play important roles in myriad biological phenomena. A distinguishing feature of IDPs is that they often mediate specific biological outcomes via multivalent weak cooperative interactions with multiple partners. Here, we show that several proteins specifically associated with processes that were key in the evolution of complex multicellularity in the lineage leading to the multicellular green alga Volvox carteri are IDPs. We suggest that, by rewiring cellular protein interaction networks, IDPs facilitated the co-option of ancestral pathways for specialized multicellular functions, underscoring the importance of IDPs in the early evolution of complex multicellularity.

Conservation and coevolution determine evolvability of different classes of disordered residues in human intrinsically disordered proteins

Structure, function, and evolution are interdependent properties of proteins. Diversity of protein functions arising from structural variations is a potential driving force behind protein evolvability. Intrinsically disordered proteins or regions (IDPs or IDRs) lack well-defined structure under normal physiological conditions, yet, they are highly functional. Increased occurrence of IDPs in eukaryotes compared to prokaryotes indicates strong correlation of protein evolution and disorderedness. IDPs generally have higher evolution rate compared to globular proteins. Structural pliability allows IDPs to accommodate multiple mutations without affecting their functional potential. Nevertheless, how evolutionary signals vary between different classes of disordered residues (DRs) in IDPs is poorly understood. This study addresses variation of evolutionary behavior in terms of residue conservation and intra-protein coevolution among structural and functional classes of DRs in IDPs. Analyses are performed on 579 human IDPs, which are classified based on length of IDRs, interacting partners and functional classes. We find short IDRs are less conserved than long IDRs or full IDPs. Functional classes which require flexibility and specificity to perform their activity comparatively evolve slower than others. Disorder promoting amino acids evolve faster than order promoting amino acids. Pro, Gly, Ile, and Phe have unique coevolving nature which further emphasizes on their roles in IDPs. This study sheds light on evolutionary footprints in different classes of DRs from human IDPs and enhances our understanding of the structural and functional potential of IDPs.

The biology of tardigrade disordered proteins in extreme stress tolerance

Abstract

Disordered proteins have long been known to help mediate tolerance to different abiotic stresses including freezing, osmotic stress, high temperatures, and desiccation in a diverse set of organisms. Recently, three novel families of intrinsically disordered proteins were identified in tardigrades, microscopic animals capable of surviving a battery of environmental extremes. These three families include the Cytoplasmic-, Secreted-, and Mitochondrial- Abundant Heat Soluble (CAHS, SAHS, and MAHS) proteins, which are collectively termed Tardigrade Disordered Proteins (TDPs). At the level of sequence conservation TDPs are unique to tardigrades, and beyond their high degree of disorder the CAHS, SAHS, and MAHS families do not resemble one another. All three families are either highly expressed constitutively, or significantly enriched in response to desiccation. In vivo, ex vivo, and in vitro experiments indicate functional roles for members of each TDP family in mitigating cellular perturbations induced by various abiotic stresses. What is currently lacking is a comprehensive and holistic understanding of the fundamental mechanisms by which TDPs function, and the properties of TDPs that allow them to function via those mechanisms. A quantitative and systematic approach is needed to identify precisely what cellular damage TDPs work to prevent, what sequence features are important for these functions, and how those sequence features contribute to the underlying mechanisms of protection. Such an approach will inform us not only about these fascinating proteins, but will also provide insights into how the sequence of a disordered protein can dictate its functional, structural, and dynamic properties.

Relevance of Electrostatic Charges in Compactness, Aggregation, and Phase Separation of Intrinsically Disordered Proteins

The abundance of intrinsic disorder in the protein realm and its role in a variety of physiological and pathological cellular events have strengthened the interest of the scientific community in understanding the structural and dynamical properties of intrinsically disordered proteins (IDPs) and regions (IDRs). Attempts at rationalizing the general principles underlying both conformational properties and transitions of IDPs/IDRs must consider the abundance of charged residues (Asp, Glu, Lys, and Arg) that typifies these proteins, rendering them assimilable to polyampholytes or polyelectrolytes. Their conformation strongly depends on both the charge density and distribution along the sequence (i.e., charge decoration) as highlighted by recent experimental and theoretical studies that have introduced novel descriptors. Published experimental data are revisited herein in the frame of this formalism, in a new and possibly unitary perspective. The physicochemical properties most directly affected by charge density and distribution are compaction and solubility, which can be described in a relatively simplified way by tools of polymer physics. Dissecting factors controlling such properties could contribute to better understanding complex biological phenomena, such as fibrillation and phase separation. Furthermore, this knowledge is expected to have enormous practical implications for the design, synthesis, and exploitation of bio-derived materials and the control of natural biological processes.

Mutations in disordered proteins as early indicators of nucleic acid changes triggering speciation

In this study, we analyze the role of different structural variants of proteins in the speciation processes. We separate human and mouse proteomes (taken as a reference) into three previously defined variants of disorder: ordered proteins (ORDPs), structured proteins with intrinsically disordered protein regions (IDPRs), and intrinsically disordered proteins (IDPs). Then, using the representation we call here Forsdyke plot, we study the correlation of DNA divergence with the corresponding protein (phenotypic) divergence in the three variants, comparing human and mouse coding sequences with their homologs from 26 eukaryotes. The parameters of the correlation are related to the speciation process. We find that the three variants of disordered proteins are differently related to the speciation process. Specifically, IDPs phenotypically diverge earlier than ORDPs and IDPRs. ORDPs diverge later but are phenotypically more reactive to nucleotide mutations than IDPRs and IDPs. Finally, IDPRs appear to diverge phenotypically later than IDPs, like ORDPs, but they are prone to accept mutations with rates that are similar to those of IDPs. We conclude that IDPs are involved in the early stages of the speciation process, whereas mutations in ORDPs, once speciation is initiated, accelerate phenotypic divergence.

The influence of intrinsic folding mechanism of an unfolded protein on the coupled folding-binding process during target recognition

Intrinsically disordered proteins (IDPs) are extensively involved in dynamic signaling processes which require a high association rate and a high dissociation rate for rapid binding/unbinding events and at the same time a sufficient high affinity for specific recognition. Although the coupled folding-binding processes of IDPs have been extensively studied, it is still impossible to predict whether an unfolded protein is suitable for molecular signaling via coupled folding-binding. In this work, we studied the interplay between intrinsic folding mechanisms and coupled folding-binding process for unfolded proteins through molecular dynamics simulations. We first studied the folding process of three representative IDPs with different folded structures, that is, c-Myb, AF9, and E3 rRNase. We found the folding free energy landscapes of IDPs are downhill or show low barriers. To further study the influence of intrinsic folding mechanism on the binding process, we modulated the folding mechanism of barnase via circular permutation and simulated the coupled folding-binding process between unfolded barnase permutant and folded barstar. Although folding of barnase was coupled to target binding, the binding kinetics was significantly affected by the intrinsic folding free energy barrier, where reducing the folding free energy barrier enhances binding rate up to two orders of magnitude. This accelerating effect is different from previous results which reflect the effect of structure flexibility on binding kinetics. Our results suggest that coupling the folding of an unfolded protein with no/low folding free energy barrier with its target binding may provide a way to achieve high specificity and rapid binding/unbinding kinetics simultaneously.

How order and disorder within paramyxoviral nucleoproteins and phosphoproteins orchestrate the molecular interplay of transcription and replication

In this review, we summarize computational and experimental data gathered so far showing that structural disorder is abundant within paramyxoviral nucleoproteins (N) and phosphoproteins (P). In particular, we focus on measles, Nipah, and Hendra viruses and highlight both commonalities and differences with respect to the closely related Sendai virus. The molecular mechanisms that control the disorder-to-order transition undergone by the intrinsically disordered C-terminal domain (NTAIL) of their N proteins upon binding to the C-terminal X domain (XD) of the homologous P proteins are described in detail. By having a significant residual disorder, NTAIL-XD complexes are illustrative examples of "fuzziness", whose possible functional significance is discussed. Finally, the relevance of N-P interactions as promising targets for innovative antiviral approaches is underscored, and the functional advantages of structural disorder for paramyxoviruses are pinpointed.

FOXP in Tetrapoda: Intrinsically Disordered Regions, Short Linear Motifs and their evolutionary significance

The FOXP subfamily is probably the most extensively characterized subfamily of the forkhead superfamily, playing important roles in development and homeostasis in vertebrates. Intrinsically disorder protein regions (IDRs) are protein segments that exhibit multiple physical interactions and play critical roles in various biological processes, including regulation and signaling. IDRs in proteins may play an important role in the evolvability of genetic systems. In this study, we analyzed 77 orthologous FOXP genes/proteins from Tetrapoda, regarding protein disorder content and evolutionary rate. We also predicted the number and type of short linear motifs (SLIMs) in the IDRs. Similar levels of protein disorder (approximately 70%) were found for FOXP1, FOXP2, and FOXP4. However, for FOXP3, which is shorter in length and has a more specific function, the disordered content was lower (30%). Mammals showed higher protein disorders for FOXP1 and FOXP4 than non-mammals. Specific analyses related to linear motifs in the four genes showed also a clear differentiation between FOXPs in mammals and non-mammals. We predicted for the first time the role of IDRs and SLIMs in the FOXP gene family associated with possible adaptive novelties within Tetrapoda. For instance, we found gain and loss of important phosphorylation sites in the Homo sapiens FOXP2 IDR regions, with possible implication for the evolution of human speech.

Oxytocin and arginine vasopressin receptor evolution: implications for adaptive novelties in placental mammals

Oxytocin receptor (OXTR) and arginine vasopressin receptors (AVPR1a, AVPR1b, and AVPR2) are paralogous genes that emerged through duplication events; along the evolutionary timeline, owing to speciation, numerous orthologues emerged as well. In order to elucidate the evolutionary forces that shaped these four genes in placental mammals and to reveal specific aspects of their protein structures, 35 species were selected. Specifically, we investigated their molecular evolutionary history and intrinsic protein disorder content, and identified the presence of short linear interaction motifs. OXTR seems to be under evolutionary constraint in placental mammals, whereas AVPR1a, AVPR1b, and AVPR2 exhibit higher evolutionary rates, suggesting that they have been under relaxed or experienced positive selection. In addition, we describe here, for the first time, that the OXTR, AVPR1a, AVPR1b, and AVPR2 mammalian orthologues preserve their disorder content, while this condition varies among the paralogues. Finally, our results reveal the presence of short linear interaction motifs, indicating possible functional adaptations related to physiological and/or behavioral taxa-specific traits.

An Overview of Practical Applications of Protein Disorder Prediction and Drive for Faster, More Accurate Predictions

Protein disordered regions are segments of a protein chain that do not adopt a stable structure. Thus far, a variety of protein disorder prediction methods have been developed and have been widely used, not only in traditional bioinformatics domains, including protein structure prediction, protein structure determination and function annotation, but also in many other biomedical fields. The relationship between intrinsically-disordered proteins and some human diseases has played a significant role in disorder prediction in disease identification and epidemiological investigations. Disordered proteins can also serve as potential targets for drug discovery with an emphasis on the disordered-to-ordered transition in the disordered binding regions, and this has led to substantial research in drug discovery or design based on protein disordered region prediction. Furthermore, protein disorder prediction has also been applied to healthcare by predicting the disease risk of mutations in patients and studying the mechanistic basis of diseases. As the applications of disorder prediction increase, so too does the need to make quick and accurate predictions. To fill this need, we also present a new approach to predict protein residue disorder using wide sequence windows that is applicable on the genomic scale.

Biophysics of protein evolution and evolutionary protein biophysics

The study of molecular evolution at the level of protein-coding genes often entails comparing large datasets of sequences to infer their evolutionary relationships. Despite the importance of a protein's structure and conformational dynamics to its function and thus its fitness, common phylogenetic methods embody minimal biophysical knowledge of proteins. To underscore the biophysical constraints on natural selection, we survey effects of protein mutations, highlighting the physical basis for marginal stability of natural globular proteins and how requirement for kinetic stability and avoidance of misfolding and misinteractions might have affected protein evolution. The biophysical underpinnings of these effects have been addressed by models with an explicit coarse-grained spatial representation of the polypeptide chain. Sequence-structure mappings based on such models are powerful conceptual tools that rationalize mutational robustness, evolvability, epistasis, promiscuous function performed by 'hidden' conformational states, resolution of adaptive conflicts and conformational switches in the evolution from one protein fold to another. Recently, protein biophysics has been applied to derive more accurate evolutionary accounts of sequence data. Methods have also been developed to exploit sequence-based evolutionary information to predict biophysical behaviours of proteins. The success of these approaches demonstrates a deep synergy between the fields of protein biophysics and protein evolution.

Dissecting partner recognition by an intrinsically disordered protein using descriptive random mutagenesis

In view of getting insights into the molecular determinants of the binding efficiency of intrinsically disordered proteins (IDPs), we used random mutagenesis. As a proof of concept, we chose the interaction between the intrinsically disordered C-terminal domain of the measles virus nucleoprotein (NTAIL) and the X domain (XD) of the viral phosphoprotein and assessed how amino acid substitutions introduced at random within NTAIL affect partner recognition. In contrast with directed evolution approaches, we did not apply any selection and used the gene library approach not for production purposes but for achieving a better understanding of the NTAIL/XD interaction. For that reason, and to differentiate our approach from similar approaches that make use of systematic (i.e., targeted) mutagenesis, we propose to call it "descriptive random mutagenesis" (DRM). NTAIL variants generated by error-prone PCR were picked at random in the absence of selection pressure and were characterized in terms of sequence and binding abilities toward XD. DRM not only identified determinants of NTAIL/XD interaction that were in good agreement with previous work but also provided new insights. In particular, we discovered that the primary interaction site is poorly evolvable in terms of binding abilities toward XD. We also identified a critical NTAIL residue whose role in stabilizing the NTAIL/XD complex had previously escaped detection, and we identified NTAIL regulatory sites that dampen the interaction while being located outside the primary interaction site. Results show that DRM is a valuable approach to study binding abilities of IDPs.