Co-Evolution of Intrinsically Disordered Proteins with Folded Partners Witnessed by Evolutionary Couplings

The molecular basis for cellular function of intrinsically disordered protein regions

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

Calcium carbonate polymorph selection in fish otoliths: A key role of phosphorylation of Starmaker-like protein

Fish otoliths are calcium carbonate biominerals found in the inner ear commonly used for tracking fish biochronologies and as a model system for biomineralization. The process of fish otolith formation is biologically controlled by numerous biomacromolecules which not only affect crystal size, shape, mechanical properties, but also selection of calcium carbonate polymorph (e.g., aragonite, vaterite). The proteinaceous control over calcium carbonate polymorph selection occurs in many other species (e.g., corals, mollusks, echinoderms) but the exact mechanism of protein interactions with calcium and carbonate ions - constituents of CaCO3 - are not fully elucidated. Herein, we focus on a native Starmaker-like protein isolated from vaterite asteriscus otoliths from Cyprinus carpio. The proteomic studies show the presence of the phosphorylated protein in vaterite otoliths. In a series of in vitro mineralization experiments with Starmaker-like, we show that native phosphorylation is a crucial determinant for the selection of a crystal's polymorphic form. This is the first report showing that the switch in calcium carbonate phase depends on the phosphorylation pattern of a single isolated protein. STATEMENT OF SIGNIFICANCE: Calcium carbonate has numerous applications in industry and medicine. However, we still do not understand the mechanism of biologically driven polymorph selection which results in specific biomineral properties. Previous work on calcium carbonate biominerals showed that either several macromolecular factors or high magnesium concentration (non-physiological) are required for proper polymorph selection (e.g., in mollusk shells, corals and otoliths). In this work, we showed for the first time that protein phosphorylation is a crucial factor for controlling the calcium carbonate crystal phase. This is important because a single protein from the otolith organic matrix could switch between polymorphs depending on the phosphorylation level. It seems that protein post-translational modifications (native, not artificial) are more important for biomolecular control of crystal growth than previously considered.

A Conserved Ribosomal Protein Has Entirely Dissimilar Structures in Different Organisms

Ribosomes from different species can markedly differ in their composition by including dozens of ribosomal proteins that are unique to specific lineages but absent in others. However, it remains unknown how ribosomes acquire new proteins throughout evolution. Here, to help answer this question, we describe the evolution of the ribosomal protein msL1/msL2 that was recently found in ribosomes from the parasitic microorganism clade, microsporidia. We show that this protein has a conserved location in the ribosome but entirely dissimilar structures in different organisms: in each of the analyzed species, msL1/msL2 exhibits an altered secondary structure, an inverted orientation of the N-termini and C-termini on the ribosomal binding surface, and a completely transformed 3D fold. We then show that this fold switching is likely caused by changes in the ribosomal msL1/msL2-binding site, specifically, by variations in rRNA. These observations allow us to infer an evolutionary scenario in which a small, positively charged, de novo-born unfolded protein was first captured by rRNA to become part of the ribosome and subsequently underwent complete fold switching to optimize its binding to its evolving ribosomal binding site. Overall, our work provides a striking example of how a protein can switch its fold in the context of a complex biological assembly, while retaining its specificity for its molecular partner. This finding will help us better understand the origin and evolution of new protein components of complex molecular assemblies-thereby enhancing our ability to engineer biological molecules, identify protein homologs, and peer into the history of life on Earth.

Chasing long-range evolutionary couplings in the AlphaFold era

Coevolution between protein residues is normally interpreted as direct contact. However, the evolutionary record of a protein sequence contains rich information that may include long-range functional couplings, couplings that report on homo-oligomeric states or even conformational changes. Due to the complexity of the sequence space and the lack of structural information on various members of a protein family, it has been difficult to effectively mine the additional information encoded in a multiple sequence alignment (MSA). Here, taking advantage of the recent release of the AlphaFold (AF) database we attempt to identify coevolutionary couplings that cannot be explained simply by spatial proximity. We propose a simple computational method that performs direct coupling analysis on a MSA and searches for couplings that are not satisfied in any of the AF models of members of the identified protein family. Application of this method on 2012 protein families suggests that ~12% of the total identified coevolving residue pairs are spatially distant and more likely to be disordered than their contacting counterparts. We expect that this analysis will help improve the quality of coevolutionary distance restraints used for structure determination and will be useful in identifying potentially functional/allosteric cross-talk between distant residues.

Distinguishing features of fold-switching proteins

Though many folded proteins assume one stable structure that performs one function, a small-but-increasing number remodel their secondary and tertiary structures and change their functions in response to cellular stimuli. These fold-switching proteins regulate biological processes and are associated with autoimmune dysfunction, severe acute respiratory syndrome coronavirus-2 infection, and more. Despite their biological importance, it is difficult to computationally predict fold switching. With the aim of advancing computational prediction and experimental characterization of fold switchers, this review discusses several features that distinguish fold-switching proteins from their single-fold and intrinsically disordered counterparts. First, the isolated structures of fold switchers are less stable and more heterogeneous than single folders but more stable and less heterogeneous than intrinsically disordered proteins (IDPs). Second, the sequences of single fold, fold switching, and intrinsically disordered proteins can evolve at distinct rates. Third, proteins from these three classes are best predicted using different computational techniques. Finally, late-breaking results suggest that single folders, fold switchers, and IDPs have distinct patterns of residue-residue coevolution. The review closes by discussing high-throughput and medium-throughput experimental approaches that might be used to identify new fold-switching proteins.

Biological soft matter: intrinsically disordered proteins in liquid-liquid phase separation and biomolecular condensates

The facts that many proteins with crucial biological functions do not have unique structures and that many biological processes are compartmentalized into the liquid-like biomolecular condensates, which are formed via liquid-liquid phase separation (LLPS) and are not surrounded by the membrane, are revolutionizing the modern biology. These phenomena are interlinked, as the presence of intrinsic disorder represents an important requirement for a protein to undergo LLPS that drives biogenesis of numerous membrane-less organelles (MLOs). Therefore, one can consider these phenomena as crucial constituents of a new IDP-LLPS-MLO field. Furthermore, intrinsically disordered proteins (IDPs), LLPS, and MLOs represent a clear link between molecular and cellular biology and soft matter and condensed soft matter physics. Both IDP and LLPS/MLO fields are undergoing explosive development and generate the ever-increasing mountain of crucial data. These new data provide answers to so many long-standing questions that it is difficult to imagine that in the very recent past, protein scientists and cellular biologists operated without taking these revolutionary concepts into account. The goal of this essay is not to deliver a comprehensive review of the IDP-LLPS-MLO field but to provide a brief and rather subjective outline of some of the recent developments in these exciting fields.

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.

An enhancer sequence in the intrinsically disordered region of FtsZ promotes polymer‐guided substrate processing by ClpXP protease

The essential bacterial division protein in Escherichia coli, FtsZ, assembles into the FtsZ-ring at midcell and recruits other proteins to the division site to promote septation. A region of the FtsZ amino acid sequence that links the conserved polymerization domain to a C-terminal protein interaction site was predicted to be intrinsically disordered and has been implicated in modulating spacing and architectural arrangements of FtsZ filaments. While the majority of cell division proteins that directly bind to FtsZ engage either the polymerization domain or the C-terminal interaction site, ClpX, the recognition and unfolding component of the bacterial ClpXP proteasome, has a secondary interaction with the predicted intrinsically disordered region (IDR) of FtsZ when FtsZ is polymerized. Here, we use NMR spectroscopy and reconstituted degradation reactions in vitro to demonstrate that this linker region is indeed disordered in solution and, further, that amino acids in the IDR of FtsZ enhance the degradation in polymer-guided interactions.

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.

Mapping the transition state for a binding reaction between ancient intrinsically disordered proteins

Intrinsically disordered protein domains often have multiple binding partners. It is plausible that the strength of pairing with specific partners evolves from an initial low affinity to a higher affinity. However, little is known about the molecular changes in the binding mechanism that would facilitate such a transition. We previously showed that the interaction between two intrinsically disordered domains, NCBD and CID, likely emerged in an ancestral deuterostome organism as a low-affinity interaction that subsequently evolved into a higher-affinity interaction before the radiation of modern vertebrate groups. Here we map native contacts in the transition states of the low-affinity ancestral and high-affinity human NCBD/CID interactions. We show that the coupled binding and folding mechanism is overall similar but with a higher degree of native hydrophobic contact formation in the transition state of the ancestral complex and more heterogeneous transient interactions, including electrostatic pairings, and an increased disorder for the human complex. Adaptation to new binding partners may be facilitated by this ability to exploit multiple alternative transient interactions while retaining the overall binding and folding pathway.

Interplay of folded domains and the disordered low-complexity domain in mediating hnRNPA1 phase separation

Liquid-liquid phase separation underlies the membrane-less compartmentalization of cells. Intrinsically disordered low-complexity domains (LCDs) often mediate phase separation, but how their phase behavior is modulated by folded domains is incompletely understood. Here, we interrogate the interplay between folded and disordered domains of the RNA-binding protein hnRNPA1. The LCD of hnRNPA1 is sufficient for mediating phase separation in vitro. However, we show that the folded RRM domains and a folded solubility-tag modify the phase behavior, even in the absence of RNA. Notably, the presence of the folded domains reverses the salt dependence of the driving force for phase separation relative to the LCD alone. Small-angle X-ray scattering experiments and coarse-grained MD simulations show that the LCD interacts transiently with the RRMs and/or the solubility-tag in a salt-sensitive manner, providing a mechanistic explanation for the observed salt-dependent phase separation. These data point to two effects from the folded domains: (i) electrostatically-mediated interactions that compact hnRNPA1 and contribute to phase separation and (ii) increased solubility at higher ionic strengths mediated by the folded domains. The interplay between disordered and folded domains can modify the dependence of phase behavior on solution conditions and can obscure signatures of physicochemical interactions underlying phase separation.

Intrinsically Disordered Protein Ensembles Shape Evolutionary Rates Revealing Conformational Patterns

Intrinsically disordered proteins (IDPs) lack stable tertiary structure under physiological conditions. The unique composition and complex dynamical behaviour of IDPs make them a challenge for structural biology and molecular evolution studies. Using NMR ensembles, we found that IDPs evolve under a strong site-specific evolutionary rate heterogeneity, mainly originated by different constraints derived from their inter-residue contacts. Evolutionary rate profiles correlate with the experimentally observed conformational diversity of the protein, allowing the description of different conformational patterns possibly related to their structure-function relationships. The correlation between evolutionary rates and contact information improves when structural information is taken not from any individual conformer or the whole ensemble, but from combining a limited number of conformers. Our results suggest that residue contacts in disordered regions constrain evolutionary rates to conserve the dynamic behaviour of the ensemble and that evolutionary rates can be used as a proxy for the conformational diversity of IDPs.

Evolutionary Study of Disorder in Protein Sequences

Intrinsically disordered proteins (IDPs) contain regions lacking intrinsic globular structure (intrinsically disordered regions, IDRs). IDPs are present across the tree of life, with great variability of IDR type and frequency even between closely related taxa. To investigate the function of IDRs, we evaluated and compared the distribution of disorder content in 10,695 reference proteomes, confirming its high variability and finding certain correlation along the Euteleostomi (bony vertebrates) lineage to number of cell types. We used the comparison of orthologs to study the function of disorder related to increase in cell types, observing that multiple interacting subunits of protein complexes might gain IDRs in evolution, thus stressing the function of IDRs in modulating protein-protein interactions, particularly in the cell nucleus. Interestingly, the conservation of local compositional biases of IDPs follows residue-type specific patterns, with E- and K-rich regions being evolutionarily stable and Q- and A-rich regions being more dynamic. We provide a framework for targeted evolutionary studies of the emergence of IDRs. We believe that, given the large variability of IDR distributions in different species, studies using this evolutionary perspective are required.

Integrating disorder in globular multidomain proteins: Fuzzy sensors and the role of SH3 domains

Intrinsically disordered proteins represent about one third of eukaryotic proteins. An additional third correspond to proteins containing folded domains as well as large intrinsically disordered regions (IDR). While IDRs may represent functionally autonomous domains, in some instances it has become clear that they provide a new layer of regulation for the activity displayed by the folded domains. The sensitivity of the conformational ensembles defining the properties of IDR to small changes in the cellular environment and the capacity to modulate this response through post-translational modifications makes IDR ideal sensors enabling continuous, integrative responses to complex cellular inputs. Folded domains (FD), on the other hand, are ideal effectors, e.g. by catalyzing enzymatic reactions or participating in binary on/off switches. In this perspective review we discuss the possible role of intramolecular fuzzy complexes to integrate the very different dynamic scales of IDR and FD, inspired on the recent observations of such dynamic complexes in Src family kinases, and we explore the possible general role of the SH3 domains connecting IDRs and FD.

Recent Advances in Computational Protocols Addressing Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDP) are abundant in the human genome and have recently emerged as major therapeutic targets for various diseases. Unlike traditional proteins that adopt a definitive structure, IDPs in free solution are disordered and exist as an ensemble of conformations. This enables the IDPs to signal through multiple signaling pathways and serve as scaffolds for multi-protein complexes. The challenge in studying IDPs experimentally stems from their disordered nature. Nuclear magnetic resonance (NMR), circular dichroism, small angle X-ray scattering, and single molecule Förster resonance energy transfer (FRET) can give the local structural information and overall dimension of IDPs, but seldom provide a unified picture of the whole protein. To understand the conformational dynamics of IDPs and how their structural ensembles recognize multiple binding partners and small molecule inhibitors, knowledge-based and physics-based sampling techniques are utilized in-silico, guided by experimental structural data. However, efficient sampling of the IDP conformational ensemble requires traversing the numerous degrees of freedom in the IDP energy landscape, as well as force-fields that accurately model the protein and solvent interactions. In this review, we have provided an overview of the current state of computational methods for studying IDP structure and dynamics and discussed the major challenges faced in this field.