AlphaFold and Implications for Intrinsically Disordered Proteins

Accurate predictions of the three-dimensional structures of proteins from their amino acid sequences have come of age. AlphaFold, a deep learning-based approach to protein structure prediction, shows remarkable success in independent assessments of prediction accuracy. A significant epoch in structural bioinformatics was the structural annotation of over 98% of protein sequences in the human proteome. Interestingly, many predictions feature regions of very low confidence, and these regions largely overlap with intrinsically disordered regions (IDRs). That over 30% of regions within the proteome are disordered is congruent with estimates that have been made over the past two decades, as intense efforts have been undertaken to generalize the structure–function paradigm to include the importance of conformational heterogeneity and dynamics. With structural annotations from AlphaFold in hand, there is the temptation to draw inferences regarding the “structures” of IDRs and their interactomes. Here, we offer a cautionary note regarding the misinterpretations that might ensue and highlight efforts that provide concrete understanding of sequence-ensemble-function relationships of IDRs. This perspective is intended to emphasize the importance of IDRs in sequence-function relationships (SERs) and to highlight how one might go about extracting quantitative SERs to make sense of how IDRs function.     

Power Law Distribution Defines Structural Disorder as a Structural Element Directly Linked with Function

Although intrinsically disordered proteins are prevalent and functionally important, it has never been asked whether structural disorder should be considered as a separate structural category on its own or merely as a lack of secondary and/or tertiary structure. We address this issue by showing that its length distribution in the human proteome follows a power law, with many short regions but also a significant incidence of very long disordered regions. This behavior is in sharp contrast with that of conventional secondary structural elements and is highly reminiscent of the distribution of tertiary structural units in proteins. We interpret this finding by the direct functional involvement of disorder, which distinguishes it from secondary structural elements and endows it with tertiary structural attributes.     

A decade and a half of protein intrinsic disorder: Biology still waits for physics

The abundant existence of proteins and regions that possess specific functions without being uniquely folded into unique 3D structures has become accepted by a significant number of protein scientists. Sequences of these intrinsically disordered proteins (IDPs) and IDP regions (IDPRs) are characterized by a number of specific features, such as low overall hydrophobicity and high net charge which makes these proteins predictable. IDPs/IDPRs possess large hydrodynamic volumes, low contents of ordered secondary structure, and are characterized by high structural heterogeneity. They are very flexible, but some may undergo disorder to order transitions in the presence of natural ligands. The degree of these structural rearrangements varies over a very wide range. IDPs/IDPRs are tightly controlled under the normal conditions and have numerous specific functions that complement functions of ordered proteins and domains. When lacking proper control, they have multiple roles in pathogenesis of various human diseases. Gaining structural and functional information about these proteins is a challenge, since they do not typically “freeze” while their “pictures are taken.” However, despite or perhaps because of the experimental challenges, these fuzzy objects with fuzzy structures and fuzzy functions are among the most interesting targets for modern protein research. This review briefly summarizes some of the recent advances in this exciting field and considers some of the basic lessons learned from the analysis of physics, chemistry, and biology of IDPs.     

Cell communication using intrinsically disordered proteins: what can syndecans say?

Because intrinsically disordered proteins are incapable of forming unique tertiary structures in isolation, their interaction with partner structures enables them to play important roles in many different biological functions. Therefore, such proteins are usually multifunctional, and their ability to perform their major function, as well as accessory functions, depends on the characteristics of a given interaction. The present paper demonstrates, using predictions from two programs, that the transmembrane proteoglycans syndecans are natively disordered because of their diverse functions and large number of interaction partners. Syndecans perform multiple functions during development, damage repair, tumor growth, angiogenesis, and neurogenesis. By mediating the binding of a large number of extracellular ligands to their receptors, these proteoglycans trigger a cascade of reactions that subsequently regulate various cell processes: cytoskeleton formation, proliferation, differentiation, adhesion, and migration. The occurrences of 20 amino acids in syndecans 1–4 from 25 animals were compared with those in 17 animal proteomes. Gly?+?Ala, Thr, Glu, and Pro were observed to predominate in the syndecans, contributing to the lack of an ordered structure. In contrast, there were many fewer amino acids in syndecans that promote an ordered structure, such as Cys, Trp, Asn, and His. In addition, a region rich in Asp has been identified between two heparan sulfate-binding sites in the ectodomains, and a region rich in Lys has been identified in the conserved C1 site of the cytoplasmic domain. These particular regions play an essential role in the various functions of syndecans due to their lack of structure.     

Intrinsically disordered regions mediate macromolecular assembly of the Slit diaphragm proteins associated with Nephrotic syndrome

The glomerular filtration barrier (GFB) of the kidney plays an instrumental role in preventing the excretion of large molecules and ensuring the formation of ultra-filtrated urine. Podocytes are essential components of GFB that provide epithelial coverage to the fine glomerular capillaries. Slit-diaphragm (SD) that forms the sole contact between adjacent foot-processes of the podocytes consists of multimeric protein assemblies. SD serves as a molecular sieve and confers size and charge-selective barrier. Nephrin, podocin, TRPC6, and CD2AP are some of the key proteins that constitute the SD. Mutations in these proteins are implicated in nephrotic syndrome and congenital nephropathies which are characterised by heavy proteinuria. The mechanism of how mutations in these proteins predispose to proteinuria is not known. Furthermore, the structural details of proteins that constitute SD are largely unknown. In this study, we built models for nephrin, CD2AP, podocin, and TRPC6 followed by docking and molecular dynamics simulations of the complex of these proteins. We speculate that the interfacial residues of SD proteins form a macromolecular complex through intrinsically disordered regions thereby conferring architectural stability to the SD, which is critical for glomerular permselectivity.     

Torches,Candles, Lamps,Lanterns, Flashlights,Spotlights, Night Vision Goggles… You Need Them All to See in Darkness

Articles assembled in the second part of this Special Issue describe some experimental and computational approaches for the structural and functional characterization of intrinsically disordered proteins. Since these tools represent specialized gear for the focused analysis of various aspects of dark proteome, they can be viewed as torches, candles, lamps, lanterns, flashlights, spotlights, night vision goggles, and other means needed to see in darkness.     

Molecular dynamics simulation of intrinsically disordered proteins

Intrinsically disordered proteins are biomolecules that do not have a definite 3D structure; therefore, their dynamical simulation cannot start from a known list of atomistic positions, such as a Protein Data Bank file. We describe a method to start a computer simulation of these proteins. The first step of the procedure is the creation of a multi-rod configuration of the molecule, derived from its primary sequence. This structure is dynamically evolved in vacuo until its gyration radius reaches the experimental average value; at this point solvent molecules, in explicit or implicit implementation, are added to the protein and a regular molecular dynamics simulation follows. We have applied this procedure to the simulation of tau, one of the largest totally disordered proteins.     

An N-terminal, 830 residues intrinsically disordered region of the cytoskeleton-regulatory protein supervillin contains Myosin II- and F-actin-binding sites

Supervillin, the largest member of the villin/gelsolin family, is a cytoskeleton regulating, peripheral membrane protein. Supervillin increases cell motility and promotes invasive activity in tumors. Major cytoskeletal interactors, including filamentous actin and myosin II, bind within the unique supervillin amino terminus, amino acids 1–830. The structural features of this key region of the supervillin polypeptide are unknown. Here, we utilize circular dichroism and bioinformatics sequence analysis to demonstrate that the N-terminal part of supervillin forms an extended intrinsically disordered region (IDR). Our combined data indicate that the N-terminus of human and bovine supervillin sequences (positions 1–830) represents an IDR, which is the largest IDR known to date in the villin/gelsolin family. Moreover, this result suggests a potentially novel mechanism of regulation of myosin II and F-actin via the intrinsically disordered N-terminal region of hub protein supervillin.     

Backbone-driven collapse in unfolded protein chains

Collapse of unfolded protein chains is an early event in folding. It affects structural properties of intrinsically disordered proteins, which take a considerable fraction of the human proteome. Collapse is generally believed to be driven by hydrophobic forces imposed by the presence of nonpolar amino acid side chains. Contributions from backbone hydrogen bonds to protein folding and stability, however, are controversial. To date, the experimental dissection of side-chain and backbone contributions has not yet been achieved because both types of interactions are integral parts of protein structure. Here, we realized this goal by applying mutagenesis and chemical modification on a set of disordered peptides and proteins. We measured the protein dimensions and kinetics of intra-chain diffusion of modified polypeptides at the level of individual molecules using fluorescence correlation spectroscopy, thereby avoiding artifacts commonly caused by aggregation of unfolded protein material in bulk. We found no contributions from side chains to collapse but, instead, identified backbone interactions as a source sufficient to form globules of native-like dimensions. The presence of backbone hydrogen bonds decreased polypeptide water solubility dramatically and accelerated the nanosecond kinetics of loop closure, in agreement with recent predictions from computer simulation. The presence of side chains, instead, slowed loop closure and modulated the dimensions of intrinsically disordered domains. It appeared that the transient formation of backbone interactions facilitates the diffusive search for productive conformations at the early stage of folding and within intrinsically disordered proteins.     

Structural plasticity and the evolution of antibody affinity and specificity

The germline precursor to the ferrochelatase antibody 7G12 was found to bind the polyether jeffamine in addition to its cognate hapten N-methylmesoporphyrin. A comparison of the X-ray crystal structures of the ligand-free germline Fab and its complex with either hapten or jeffamine reveals that the germline antibody undergoes significant conformational changes upon the binding of these two structurally distinct ligands, which lead to increased antibody-ligand complementarity. The five somatic mutations introduced during affinity maturation lead to enhanced binding affinity for hapten and a loss in affinity for jeffamine. Moreover, a comparison of the crystal structures of the germline and affinity-matured antibodies reveals that somatic mutations not only fix the optimal binding site conformation for the hapten, but also introduce interactions that interfere with the binding of non-hapten molecules. The structural plasticity of this germline antibody and the structural effects of the somatic mutations that result in enhanced affinity and specificity for hapten likely represent general mechanisms used by the immune response, and perhaps primitive proteins, to evolve high affinity, selective receptors for so many distinct chemical structures.     

Mutational analysis of the affinity maturation of antibody 48G7

The affinity maturation of antibody 48G7 from its germline predecessor 48G7g has been studied at a molecular level through a combination of structural and biochemical means. Each of the nine somatic mutations accumulated during affinity maturation has been assessed for gain or loss of function in both the germline and affinity-matured antibodies. Individual somatic mutations were found to be either positive or neutral in their effects on affinity for hapten JWJ1, with a marked context-dependence for some sites of mutation. In a number of cases significant cooperativity was found between pairs of somatically mutated residues. Interpretation of the structural changes introduced by many of the point mutations has been possible due to the availability of high-resolution crystal structures of 48G7g and 48G7, and mechanisms by which these structural changes may result in enhanced affinity for hapten have been identified. Precise dissection of structure-function relationships in this system provides additional insights into the role of cooperativity in the evolution of antibody affinity. Comparison of 48G7 with previously characterized systems provides a varied view of the structure-function mechanisms by which the humoral immune system produces large increases in affinity.     

Mapping the amino acid properties of constituent nucleoporins onto the yeast nuclear pore complex

Visualization of molecular structures aids in the understanding of structural and functional roles of biological macromolecules. Macromolecular transport between the cell nucleus and cytoplasm is facilitated by the nuclear pore complex (NPC). The ring structure of the NPC is large and contains several distinct proteins (nucleoporins) which function as a selective gate for the passage of certain molecules into and out of the nucleus. In this note we demonstrate the utility of a python code that allows direct mapping of the physiochemical properties of the constituent nucleoporins on the scaffold of the yeast NPC׳s cytoplasmic view. We expect this tool to be useful for researchers to visualize the NPC based on their physiochemical properties and how it alters when specific mutations are introduced in one or more of the nucleoporins. The code developed using Python is available freely from the authors.     

Structural disorder: a tool for housekeeping proteins performing tissue-specific interactions

An interaction between a pair of proteins unique for a particular tissue is denoted as a tissue-specific interaction (TSI). Tissue-specific (TS) proteins always perform TSIs with a limited number of interacting partners. However, it has been claimed that housekeeping (HK) proteins frequently take part in TSIs. This is actually an unusual phenomenon. How a single HK protein mediates TSIs – remains an interesting yet an unsolved question. We have hypothesized that HK proteins have attained a high degree of structural flexibility to modulate TSIs efficiently. We have observed that HK proteins are selected to be intrinsically disordered compared to TS proteins. Therefore, the purposeful adaptation of structural disorder brings out special advantages for HK proteins compared to TS proteins. We have demonstrated that TSIs may play vital roles in shaping the molecular adaptation of disordered regions within HK proteins. We also have noticed that HK proteins, mediating a huge number of TSIs, have a greater portion of their interacting interfaces overlapped with the adjacent disordered segment. Moreover, these HK proteins, mediating TSIs, preferably adapt single domain (SD). We have concluded that HK proteins adapt a high degree of structural flexibility to mediate TSIs. Besides, having a SD along with structural flexibility is more economic than maintaining multiple domains with a rigid structure. This assists them in attaining various structural conformations upon binding to their partners, thereby designing an economically optimum molecular system.     

Structural Capacitance in Protein Evolution and Human Diseases

Canonical mechanisms of protein evolution include the duplication and diversification of pre-existing folds through genetic alterations that include point mutations, insertions, deletions, and copy number amplifications, as well as post-translational modifications that modify processes such as folding efficiency and cellular localization. Following a survey of the human mutation database, we have identified an additional mechanism that we term “structural capacitance,” which results in the de novo generation of microstructure in previously disordered regions. We suggest that the potential for structural capacitance confers select proteins with the capacity to evolve over rapid timescales, facilitating saltatory evolution as opposed to gradualistic canonical Darwinian mechanisms. Our results implicate the elements of protein microstructure generated by this distinct mechanism in the pathogenesis of a wide variety of human diseases. The benefits of rapidly furnishing the potential for evolutionary change conferred by structural capacitance are consequently counterbalanced by this accompanying risk. The phenomenon of structural capacitance has implications ranging from the ancestral diversification of protein folds to the engineering of synthetic proteins with enhanced evolvability.     

Biological function in a non-native partially folded state of a protein

As structural flexibility is known to be required for enzyme catalysis and pattern recognition and a significant fraction of eukaryotic proteins appear to be unfolded or contain unstructured regions, biological activity of conformational states distinct from fully folded structures could be more common than previously thought. By applying a procedure that allows the recovery of enzymatic activity to be monitored in real time, we show that a non-native state populated transiently during folding of the acylphosphatase from Sulfolobus solfataricus is enzymatically active. The structural characterization of this partially folded state reveals that enzymatic activity is possible even if the catalytic site is structurally heterogeneous, whereas the remainder of the structure acts as a scaffold. These results extend the spectrum of biological functions carried out in the absence of a folded state to include enzyme catalysis.     

Structural divergence is more extensive than sequence divergence for a family of intrinsically disordered proteins

The p53 transactivation domain (p53TAD) is an intrinsically disordered protein (IDP) domain that undergoes coupled folding and binding when interacting with partner proteins like the E3 ligase, MDM2, and the 70 kDa subunit of replication protein A, RPA70. The secondary structure and dynamics of six closely related mammalian homologues of p53TAD were investigated using nuclear magnetic resonance (NMR) spectroscopy. Differences in both transient secondary structure and backbone dynamics were observed for the homologues. Many of these differences were localized to the binding sites for MDM2 and RPA70. The amount of transient helical secondary structure observed for the MDM2 binding site was lower for the dog and mouse homologues, compared with human, and the amount of transient helical secondary structure observed for the RPA70 binding site was higher for guinea pig and rabbit, compared with human. Differences in the amount of transient helical secondary structure observed for the MDM2 binding site were directly related to amino acid substitutions occurring on the solvent exposed side of the amphipathic helix that forms during the p53TAD/MDM2 interaction. Differences in the amount of transient helical secondary structure were not as easily explained for the RPA70 binding site because of its extensive sequence divergence. Clustering analysis shows that the divergence in the transient secondary structure of the p53TAD homologues exceeds the amino acid sequence divergence. In contrast, strong correlations were observed between the backbone dynamics of the homologues and the sequence identity matrix, suggesting that the dynamic behavior of IDPs is a conserved evolutionary feature. Proteins 2013; 81:1686–1698. © 2013 Wiley Periodicals, Inc.     

Investigation of deleterious effects of nsSNPs in the POT1 gene: a structural genomics-based approach to understand the mechanism of cancer development

Protection of telomere 1 (POT1) is one of the key components of shelterin complex, implicated in maintaining the telomere homeostasis, and thus stability of the eukaryotic genome. A large number of non-synonymous single nucleotide polymorphisms (nsSNPs) in the POT1 gene have been reported to cause varieties of human diseases, including cancer. In recent years, a number of mutations in POT1 has been markedly increased, and interpreting the effect of these large numbers of mutations to understand the mechanism of associated diseases seems impossible using experimental approaches. Herein, we employ varieties of computational methods such as PROVEAN, PolyPhen-2, SIFT, PoPMuSiC, SDM2, STRUM, and MAESTRO to identify the effects of 387 nsSNPs on the structure and function of POT1 protein. We have identified about 183 nsSNPs as deleterious and termed them as “high-confidence nsSNPs.” Distribution of these high-confidence nsSNPs demonstrates that the mutation in oligonucleotide binding domain 1 is highly deleterious (one in every three nsSNPs), and high-confidence nsSNPs show a strong correlation with residue conservation. The structure analysis provides a detailed insights into the structural changes occurred in consequence of conserved mutations which lead to the cancer progression. This study, for the first time, offers a newer prospective on the role of POT1 mutations on the structure, function, and their relation to associated diseases.     

Effect of natural polymorphism on structure and function of the Yersinia pestis outer membrane porin F (OmpF protein): a computational study

The Yersinia pestis outer membrane porin F (OmpF) is a transmembrane protein located in the outer membrane of this Gram-negative bacterium which is the causative agent of plague, where it plays a significant role in controlling the selective permeability of the membrane. The amino acid sequences of OmpF proteins from 48 Y. pestis strains representing all currently available phylogenetic groups of this Gram-negative bacterium were recently deduced. Comparison of these amino acid sequences revealed that the OmpF can be present in four isoforms, the pestis-pestis type, and the pestis-microtus types I, II, and III. OmpF of the most recent pestis-pestis type has an alanine residue at the position 148, where all the pestis-microtus types have threonine there (T148A polymorphism). The variability of different pestis-microtus types is caused by an additional polymorphism at the 193rd position, where the OmpFs of the pestis-microtus type II and type III have isoleucine-glycine (IG⁺₁₉₃) or isoleucine-glycine-isoleucine-glycine (IGIG⁺₁₉₃) insertions, respectively (IG⁺₁₉₃ and IGIG⁺₁₉₃ polymorphism). To investigate potential effects of these sequence polymorphisms on the structural properties of the OmpF protein, we conducted multi-level computational analysis of its isoforms. Analysis of the I-TASSER-generated 3D-models revealed that the Yersinia OmpF is very similar to other non-specific enterobacterial porins. The T148A polymorphism affected a loop located in the external vestibule of the OmpF channel, whereas IG⁺₁₉₃ and IGIG⁺₁₉₃ polymorphisms affected one of its β-strands. Our analysis also suggested that polymorphism has moderate effect on the predicted local intrinsic disorder predisposition of OmpF, but might have some functional implementations.     

Germline fitness-based scoring of cancer mutations

A key goal in cancer research is to find the genomic alterations that underlie malignant cells. Genomics has proved successful in identifying somatic variants at a large scale. However, it has become evident that a typical cancer exhibits a heterogenous mutation pattern across samples. Cases where the same alteration is observed repeatedly seem to be the exception rather than the norm. Thus, pinpointing the key alterations (driver mutations) from a background of variations with no direct causal link to cancer (passenger mutations) is difficult. Here we analyze somatic missense mutations from cancer samples and their healthy tissue counterparts (germline mutations) from the viewpoint of germline fitness. We calibrate a scoring system from protein domain alignments to score mutations and their target loci. We show first that this score predicts to a good degree the rate of polymorphism of the observed germline variation. The scoring is then applied to somatic mutations. We show that candidate cancer genes prone to copy number loss harbor mutations with germline fitness effects that are significantly more deleterious than expected by chance. This suggests that missense mutations play a driving role in tumor suppressor genes. Furthermore, these mutations fall preferably onto loci in sequence neighborhoods that are high scoring in terms of germline fitness. In contrast, for somatic mutations in candidate onco genes we do not observe a statistically significant effect. These results help to inform how to exploit germline fitness predictions in discovering new genes and mutations responsible for cancer.