Recognition pliability is coupled to structural heterogeneity: a calmodulin intrinsically disordered binding region complex

Protein interactions within regulatory networks should adapt in a spatiotemporal-dependent dynamic environment, in order to process and respond to diverse and versatile cellular signals. However, the principles governing recognition pliability in protein complexes are not well understood. We have investigated a region of the intrinsically disordered protein myelin basic protein (MBP(145-165)) that interacts with calmodulin, but that also promiscuously binds other biomolecules (membranes, modifying enzymes). To characterize this interaction, we implemented an NMR spectroscopic approach that calculates, for each conformation of the complex, the maximum occurrence based on recorded pseudocontact shifts and residual dipolar couplings. We found that the MBP(145-165)-calmodulin interaction is characterized by structural heterogeneity. Quantitative comparative analysis indicated that distinct conformational landscapes of structural heterogeneity are sampled for different calmodulin-target complexes. Such structural heterogeneity in protein complexes could potentially explain the way that transient and promiscuous protein interactions are optimized and tuned in complex regulatory networks.     

Codon usage is less optimized in eukaryotic gene segments encoding intrinsically disordered regions than in those encoding structural domains

Codon usage tends to be optimized in highly expressed genes. A plausible explanation for this phenomenon is that translational accuracy is increased in highly expressed genes with infrequent use of rare codons. Besides structural domains (SDs), eukaryotic proteins generally have intrinsically disordered regions (IDRs) that by themselves do not assume unique three-dimensional structures. As IDRs are free from structural constraint, they can probably accommodate more translational errors than SDs can. Thus, codon usage in IDRs is likely to be less optimized than that in SDs. Codon usage in all the genes of seven eukaryotes was examined in terms of both tRNA adaptation index and codon adaptation index. Different amino acid compositions in different protein regions were taken into account in calculating expected adaptation indices, to which observed indices were compared. Codon usage is less optimized in gene regions encoding IDRs than in those corresponding to SDs. The finding does not depend on whether IDRs are located at the N-terminus, in the middle, or at the C-terminus of proteins. Furthermore, the observation remains unchanged in two different algorithms used to predict IDRs in proteins. The result is consistent with the idea that IDRs tolerate more translational errors than SDs.     

SLiMPrints: conservation-based discovery of functional motif fingerprints in intrinsically disordered protein regions

Large portions of higher eukaryotic proteomes are intrinsically disordered, and abundant evidence suggests that these unstructured regions of proteins are rich in regulatory interaction interfaces. A major class of disordered interaction interfaces are the compact and degenerate modules known as short linear motifs (SLiMs). As a result of the difficulties associated with the experimental identification and validation of SLiMs, our understanding of these modules is limited, advocating the use of computational methods to focus experimental discovery. This article evaluates the use of evolutionary conservation as a discriminatory technique for motif discovery. A statistical framework is introduced to assess the significance of relatively conserved residues, quantifying the likelihood a residue will have a particular level of conservation given the conservation of the surrounding residues. The framework is expanded to assess the significance of groupings of conserved residues, a metric that forms the basis of SLiMPrints (short linear motif fingerprints), a de novo motif discovery tool. SLiMPrints identifies relatively overconstrained proximal groupings of residues within intrinsically disordered regions, indicative of putatively functional motifs. Finally, the human proteome is analysed to create a set of highly conserved putative motif instances, including a novel site on translation initiation factor eIF2A that may regulate translation through binding of eIF4E.     

Deciphering the cause of evolutionary variance within intrinsically disordered regions in human proteins

Why the intrinsically disordered regions evolve within human proteome has became an interesting question for a decade. Till date, it remains an unsolved yet an intriguing issue to investigate why some of the disordered regions evolve rapidly while the rest are highly conserved across mammalian species. Identifying the key biological factors, responsible for the variation in the conservation rate of different disordered regions within the human proteome, may revisit the above issue. We emphasized that among the other biological features (multifunctionality, gene essentiality, protein connectivity, number of unique domains, gene expression level and expression breadth) considered in our study, the number of unique protein domains acts as a strong determinant that negatively influences the conservation of disordered regions. In this context, we justified that proteins having a fewer types of domains preferably need to conserve their disordered regions to enhance their structural flexibility which in turn will facilitate their molecular interactions. In contrast, the selection pressure acting on the stretches of disordered regions is not so strong in the case of multi-domains proteins. Therefore, we reasoned that the presence of conserved disordered stretches may compensate the functions of multiple domains within a single domain protein. Interestingly, we noticed that the influence of the unique domain number and expression level acts differently on the evolution of disordered regions from that of well-structured ones.     

Proteome-wide evidence for enhanced positive Darwinian selection within intrinsically disordered regions in proteins

Background  

Understanding the adaptive changes that alter the function of proteins during evolution is an important question for biology and medicine. The increasing number of completely sequenced genomes from closely related organisms, as well as individuals within species, facilitates systematic detection of recent selection events by means of comparative genomics.     

Nonnative interactions in coupled folding and binding processes of intrinsically disordered proteins

Proteins function by interacting with other molecules, where both native and nonnative interactions play important roles. Native interactions contribute to the stability and specificity of a complex, whereas nonnative interactions mainly perturb the binding kinetics. For intrinsically disordered proteins (IDPs), which do not adopt rigid structures when being free in solution, the role of nonnative interactions may be more prominent in binding processes due to their high flexibilities. In this work, we investigated the effect of nonnative hydrophobic interactions on the coupled folding and binding processes of IDPs and its interplay with chain flexibility by conducting molecular dynamics simulations. Our results showed that the free-energy profiles became rugged, and intermediate states occurred when nonnative hydrophobic interactions were introduced. The binding rate was initially accelerated and subsequently dramatically decreased as the strength of the nonnative hydrophobic interactions increased. Both thermodynamic and kinetic analysis showed that disordered systems were more readily affected by nonnative interactions than ordered systems. Furthermore, it was demonstrated that the kinetic advantage of IDPs ("fly-casting" mechanism) was enhanced by nonnative hydrophobic interactions. The relationship between chain flexibility and protein aggregation is also discussed.     

Terminal regions of flagellin are disordered in solution

Limited proteolysis of flagellin from Salmonella typhimurium SJW1103 by subtilisin, trypsin and thermolysin results in homologous degradation patterns. The terminal regions of flagellin are very sensitive to proteolysis. These parts are degraded into small oligopeptides at the very early stage of a mild digestion that yields a relatively stable fragment with a molecular weight of 40,000. Further proteolytic degradation results in a stable 27,000 Mr fragment. The 40,000 Mr tryptic fragment has been identified as residues 67 to 446 of the flagellin sequence, while the 27,000 Mr fragment involves the 179 to 418 segment. The NH2-terminal sequence positions for the corresponding fragments produced by subtilisin are 60 and 174 for the 40,000 Mr and 27,000 Mr fragments, respectively. The fragments lost their polymerizing ability. Structural properties of flagellin and its 40,000 Mr tryptic fragment were compared by circular dichroism spectroscopy and differential scanning calorimetry. Analysis of the calorimetric melting profiles suggests that terminal parts of flagellin have no significant internal stability and they are in extensive contact with water. However, these regions contain some secondary structure, probably alpha-helices, as revealed by comparison of the circular dichroic spectra in the far-ultraviolet region. Our results indicate that, although the terminal regions of flagellin may contain some alpha-helical secondary structure of marginal stability, they have no compact ordered tertiary structure in solution. On the contrary, the central region of the molecule involves at least two compact structural units.     

Intrinsically disordered proteins display no preference for chaperone binding in vivo

Intrinsically disordered/unstructured proteins (IDPs) are extremely sensitive to proteolysis in vitro, but show no enhanced degradation rates in vivo. Their existence and functioning may be explained if IDPs are preferentially associated with chaperones in the cell, which may offer protection against degradation by proteases. To test this inference, we took pairwise interaction data from high-throughput interaction studies and analyzed to see if predicted disorder correlates with the tendency of chaperone binding by proteins. Our major finding is that disorder predicted by the IUPred algorithm actually shows negative correlation with chaperone binding in E. coli, S. cerevisiae, and metazoa species. Since predicted disorder positively correlates with the tendency of partner binding in the interactome, the difference between the disorder of chaperone-binding and non-binding proteins is even more pronounced if normalized to their overall tendency to be involved in pairwise protein–protein interactions. We argue that chaperone binding is primarily required for folding of globular proteins, as reflected in an increased preference for chaperones of proteins in which at least one Pfam domain exists. In terms of the functional consequences of chaperone binding of mostly disordered proteins, we suggest that its primary reason is not the assistance of folding, but promotion of assembly with partners. In support of this conclusion, we show that IDPs that bind chaperones also tend to bind other proteins.     

Microsecond molecular dynamics simulations of intrinsically disordered proteins involved in the oxidative stress response

Intrinsically disordered proteins (IDPs) are abundant in cells and have central roles in protein-protein interaction networks. Interactions between the IDP Prothymosin alpha (ProTα) and the Neh2 domain of Nuclear factor erythroid 2-related factor 2 (Nrf2), with a common binding partner, Kelch-like ECH-associated protein 1(Keap1), are essential for regulating cellular response to oxidative stress. Misregulation of this pathway can lead to neurodegenerative diseases, premature aging and cancer. In order to understand the mechanisms these two disordered proteins employ to bind to Keap1, we performed extensive 0.5-1.0 microsecond atomistic molecular dynamics (MD) simulations and isothermal titration calorimetry experiments to investigate the structure/dynamics of free-state ProTα and Neh2 and their thermodynamics of bindings. The results show that in their free states, both ProTα and Neh2 have propensities to form bound-state-like β-turn structures but to different extents. We also found that, for both proteins, residues outside the Keap1-binding motifs may play important roles in stabilizing the bound-state-like structures. Based on our findings, we propose that the binding of disordered ProTα and Neh2 to Keap1 occurs synergistically via preformed structural elements (PSEs) and coupled folding and binding, with a heavy bias towards PSEs, particularly for Neh2. Our results provide insights into the molecular mechanisms Neh2 and ProTα bind to Keap1, information that is useful for developing therapeutics to enhance the oxidative stress response.     

Preformed structural elements feature in partner recognition by intrinsically unstructured proteins

Intrinsically unstructured proteins (IUPs) are devoid of extensive structural order but often display signs of local and limited residual structure. To explain their effective functioning, we reasoned that such residual structure can be crucial in their interactions with their structured partner(s) in a way that preformed structural elements presage their final conformational state. To check this assumption, a database of 24 IUPs with known 3D structures in the bound state has been assembled and the distribution of secondary structure elements and backbone torsion angles have been analysed. The high proportion of residues in coil conformation and with phi, psi angles in the disallowed regions of the Ramachandran map compared to the reference set of globular proteins shows that IUPs are not fully ordered even in their bound form. To probe the effect of partner proteins on IUP folding, inherent conformational preferences of IUP sequences have been assessed by secondary structure predictions using the GOR, ALB and PROF algorithms. The accuracy of predicting secondary structure elements of IUPs is similar to that of their partner proteins and is significantly higher than the corresponding values for random sequences. We propose that strong conformational preferences mark regions in IUPs (mostly helices), which correspond to their final structural state, while regions with weak conformational preferences represent flexible linkers between them. In our interpretation, preformed elements could serve as initial contact points, the binding of which facilitates the reeling of the flexible regions onto the template. This finding implies that IUPs draw a functional advantage from preformed structural elements, as they enable their facile, kinetically and energetically less demanding, interaction with their physiological partner.     

Fatty acids covalently bound to alpha-hemolysin of Escherichia coli are involved in the molten globule conformation: implication of disordered regions in binding promiscuity

Alpha-hemolysin (HlyA) is a pore-forming toxin secreted by pathogenic strains of Escherichia coli. The toxin is synthesized as a protoxin, ProHlyA, which is matured in the cytosol to the active form by acylation at two internal lysines, K563 and K689 (HlyA). It is widely known that the presence of fatty acids is crucial for the hemolytic and cytotoxic effects of the toxin. However, no detailed physicochemical characterization of the structural changes produced by fatty acids in the soluble protein prior to membrane binding has been carried out to date. The effects of chemical denaturants, the ANS binding parameters (Kd and n) and the sensitivity to proteases were compared between the acylated and unacylated protein forms HlyA and ProHlyA. Our results are consistent with a molten globular form of the acylated protein. Moreover, because molten globule proteins are intrinsically disordered proteins, using disorder prediction analyses, we show that HlyA contains 9 regions composed of 10-30 natively disordered amino acids. We propose that this conformation induced by covalently bound fatty acids might provide HlyA with the ability to bind to a variety of molecules during its action mechanism.     

The effect of small molecules in modulating the chaperone activity of alphaB-crystallin against ordered and disordered protein aggregation

Protein aggregation can proceed via disordered or ordered mechanisms, with the latter being associated with amyloid fibril formation, which has been linked to a number of debilitating conditions including Alzheimer's, Parkinson's and Creutzfeldt-Jakob diseases. Small heat-shock proteins (sHsps), such as alphaB-crystallin, act as chaperones to prevent protein aggregation and are thought to play a key role in the prevention of protein-misfolding diseases. In this study, we have explored the potential for small molecules such as arginine and guanidine to affect the chaperone activity of alphaB-crystallin against disordered (amorphous) and ordered (amyloid fibril) forms of protein aggregation. The effect of these additives is highly dependent upon the target protein undergoing aggregation. Importantly, our results show that the chaperone action of alphaB-crystallin against aggregation of the disease-related amyloid fibril forming protein alpha-synucleinA53T is enhanced in the presence of arginine and similar positively charged compounds (such as lysine and guanidine). Thus, our results suggest that target protein identity plays a critical role in governing the effect of small molecules on the chaperone action of sHsps. Significantly, small molecules that regulate the activity of sHsps may provide a mechanism to protect cells from the toxic protein aggregation that is associated with some protein-misfolding diseases.     

The relationships among microRNA regulation, intrinsically disordered regions, and other indicators of protein evolutionary rate

Many indicators of protein evolutionary rate have been proposed, but some of them are interrelated. The purpose of this study is to disentangle their correlations. We assess the strength of each indicator by controlling for the other indicators under study. We find that the number of microRNA (miRNA) types that regulate a gene is the strongest rate indicator (a negative correlation), followed by disorder content (the percentage of disordered regions in a protein, a positive correlation); the strength of disorder content as a rate indicator is substantially increased after controlling for the number of miRNA types. By dividing proteins into lowly and highly intrinsically disordered proteins (L-IDPs and H-IDPs), we find that proteins interacting with more H-IDPs tend to evolve more slowly, which largely explains the previous observation of a negative correlation between the number of protein-protein interactions and evolutionary rate. Moreover, all of the indicators examined here, except for the number of miRNA types, have different strengths in L-IDPs and in H-IDPs. Finally, the number of phosphorylation sites is weakly correlated with the number of miRNA types, and its strength as a rate indicator is substantially reduced when other indicators are considered. Our study reveals the relative strength of each rate indicator and increases our understanding of protein evolution.     

Differences in the number of intrinsically disordered regions between yeast duplicated proteins, and their relationship with functional divergence

Background

Intrinsically disordered regions are enriched in short interaction motifs that play a critical role in many protein-protein interactions. Since new short interaction motifs may easily evolve, they have the potential to rapidly change protein interactions and cellular signaling. In this work we examined the dynamics of gain and loss of intrinsically disordered regions in duplicated proteins to inspect if changes after genome duplication can create functional divergence. For this purpose we used Saccharomyces cerevisiae and the outgroup species Lachancea kluyveri.

Principal Findings

We find that genes duplicated as part of a genome duplication (ohnologs) are significantly more intrinsically disordered than singletons (p<2.2^e-16, Wilcoxon), reflecting a preference for retaining intrinsically disordered proteins in duplicate. In addition, there have been marked changes in the extent of intrinsic disorder following duplication. A large number of duplicated genes have more intrinsic disorder than their L. kluyveri ortholog (29% for duplicates versus 25% for singletons) and an even greater number have less intrinsic disorder than the L. kluyveri ortholog (37% for duplicates versus 25% for singletons). Finally, we show that the number of physical interactions is significantly greater in the more intrinsically disordered ohnolog of a pair (p = 0.003, Wilcoxon).

Conclusion

This work shows that intrinsic disorder gain and loss in a protein is a mechanism by which a genome can also diverge and innovate. The higher number of interactors for proteins that have gained intrinsic disorder compared with their duplicates may reflect the acquisition of new interaction partners or new functional roles.     

The N-terminal domain of the Drosophila histone mRNA binding protein, SLBP, is intrinsically disordered with nascent helical structure

Stem-loop binding protein (SLBP) is a 31 kDa protein that is central to the regulation of histone mRNAs and is highly conserved in metazoans. In vertebrates, the N-terminal domain of SLBP has sequence determinants necessary for histone mRNA translation, SLBP degradation, cyclin binding, and histone mRNA import. We have used high-resolution NMR spectroscopy and circular dichroism to characterize the structural and dynamic features of this domain of SLBP from Drosophila (dSLBP). We report that the N-terminal domain of dSLBP is stably unfolded but has nascent helical structure at physiological pH and native-like solution conditions. The conformational and dynamic properties of the isolated domain are mimicked in a longer 175-residue region of the N-terminus, as well as in the full-length protein. Complete resonance assignments, secondary structure propensity, and motional properties of a 91-residue N-terminal domain (G17-K108) of dSLBP are reported here. The deviation of (1)H(alpha), (13)C(alpha), and (13)C(beta) chemical shifts from random coil reveals that there are four regions between residues I28-A45, S50-L57, S66-G75, and F91-N96 that have helical propensity. These regions also have small but positive heteronuclear NOEs, interresidue d(NN) NOEs, and small but significant protection from solvent exchange. However the lack of medium- and long-range NOEs in 3D (15)N- and (13)C-edited spectra, fast amide proton exchange rates (all greater than 1 s(-1)), and long (15)N relaxation (T(1), T(2)) times suggest that the domain from dSLBP does not adopt a well-defined tertiary fold. The backbone residual dipolar couplings (RDCs) for this domain are small and lie close to 0 Hz (+/-2 Hz) for most residues with no well-defined periodicity. The implications of this unfolded state for the function of dSLBP in regulating histone metabolism are discussed.