Energetics-based protein profiling on a proteomic scale: identification of proteins resistant to proteolysis

Native states of proteins are flexible, populating more than just the unique native conformation. The energetics and dynamics resulting from this conformational ensemble are inherently linked to protein function and regulation. Proteolytic susceptibility is one feature determined by this conformational energy landscape. As an attempt to investigate energetics of proteins on a proteomic scale, we challenged the Escherichia coli proteome with extensive proteolysis and determined which proteins, if any, have optimized their energy landscape for resistance to proteolysis. To our surprise, multiple soluble proteins survived the challenge. Maltose binding protein, a survivor from thermolysin digestion, was characterized by in vitro biophysical studies to identify the physical origin of proteolytic resistance. This experimental characterization shows that kinetic stability is responsible for the unusual resistance in maltose binding protein. The biochemical functions of the identified survivors suggest that many of these proteins may have evolved extreme proteolytic resistance because of their critical roles under stressed conditions. Our results suggest that under functional selection proteins can evolve extreme proteolysis resistance by modulating their conformational energy landscapes without the need to invent new folds, and that proteins can be profiled on a proteomic scale according to their energetic properties by using proteolysis as a structural probe.     

How expression level influences the disorderness of proteins

Intrinsically disordered proteins are very common in eukaryotes and thus understanding functional roles and factors which influence protein disorderness becomes very important. In this work, we ask whether global properties not directly related to the function of the proteins, like expression level and avoidance of aggregation, influence disorderness of proteins. We found that proteins expressed at higher levels tend to be less disordered, even within the same functional class. We also found that the correlation between expression level and evolutionary rate was significantly reduced for disordered proteins indicating the role of disorderness in preventing aggregation of highly expressed proteins, which are more susceptible to misfolding due to translational errors. We reconcile these seemingly opposing results based on the observation that the correlation between expression level and disorderness was significantly less for proteins involved in binding functions, suggesting that highly expressed proteins involved in binding functions utilize disordered regions to avoid aggregation. Our results show that disorderness is not just influenced by functional properties of proteins, but also by properties not directly related to their functions like expression level and avoidance of aggregation.     

Calcium-induced folding of intrinsically disordered repeat-in-toxin (RTX) motifs via changes of protein charges and oligomerization states

Ligand-induced disorder-to-order transition plays a key role in the biological functions of many proteins that contain intrinsically disordered regions. This trait is exhibited by so-called RTX (repeat-in-toxin) motifs found in many virulence factors secreted by numerous gram-negative pathogenic bacteria: RTX proteins are natively disordered in the absence of calcium but fold upon calcium binding. The adenylate cyclase toxin (CyaA) produced by Bordetella pertussis, the causative agent of whooping cough, contains ～40 RTX motifs organized in five successive blocks separated by non-RTX flanking regions. This RTX domain mediates toxin binding to its eukaryotic cell receptor. We previously showed that the last block of the RTX domain, block V, which is critical for CyaA toxicity, exhibits the hallmarks of intrinsically disordered proteins in the absence of calcium. Moreover, the C-terminal flanking region of CyaA block V is required for its calcium-induced folding. Here, we describe a comprehensive analysis of the hydrodynamic and electrophoretic properties of several block V RTX polypeptides that differ in the presence and/or length of the flanking regions. Our results indicate that the length of the C-terminal flanking region not only controls the calcium-induced folding but also the calcium-induced multimerization of the RTX polypeptides. Moreover, we showed that calcium binding is accompanied by a strong reduction of the net charge of the RTX polypeptides. These data indicate that the disorder-to-order transition in RTX proteins is controlled by a calcium-induced change of the polypeptide charges and stabilized by multimerization.     

Prediction of Protein Binding Regions in Disordered Proteins

Many disordered proteins function via binding to a structured partner and undergo a disorder-to-order transition. The coupled folding and binding can confer several functional advantages such as the precise control of binding specificity without increased affinity. Additionally, the inherent flexibility allows the binding site to adopt various conformations and to bind to multiple partners. These features explain the prevalence of such binding elements in signaling and regulatory processes. In this work, we report ANCHOR, a method for the prediction of disordered binding regions. ANCHOR relies on the pairwise energy estimation approach that is the basis of IUPred, a previous general disorder prediction method. In order to predict disordered binding regions, we seek to identify segments that are in disordered regions, cannot form enough favorable intrachain interactions to fold on their own, and are likely to gain stabilizing energy by interacting with a globular protein partner. The performance of ANCHOR was found to be largely independent from the amino acid composition and adopted secondary structure. Longer binding sites generally were predicted to be segmented, in agreement with available experimentally characterized examples. Scanning several hundred proteomes showed that the occurrence of disordered binding sites increased with the complexity of the organisms even compared to disordered regions in general. Furthermore, the length distribution of binding sites was different from disordered protein regions in general and was dominated by shorter segments. These results underline the importance of disordered proteins and protein segments in establishing new binding regions. Due to their specific biophysical properties, disordered binding sites generally carry a robust sequence signal, and this signal is efficiently captured by our method. Through its generality, ANCHOR opens new ways to study the essential functional sites of disordered proteins.     

Interaction between Intrinsically Disordered Proteins Frequently Occurs in a Human Protein-Protein Interaction Network

Intrinsic protein disorder is a widespread phenomenon characterised by a lack of stable three-dimensional structures and is considered to play an important role in protein-protein interactions (PPIs). This study examined the genome-wide preference of disorder in PPIs by using exhaustive disorder prediction in human PPIs. We categorised the PPIs into three types (interaction between disordered proteins, interaction between structured proteins, and interaction between a disordered protein and a structured protein) with regard to the flexibility of molecular recognition and compared these three interaction types in an existing human PPI network with those in a randomised network. Although the structured regions were expected to become the identifiers for binding recognition, this comparative analysis revealed unexpected results. The occurrence of interactions between disordered proteins was significantly frequent, and that between a disordered protein and a structured protein was significantly infrequent. We found that this propensity was much stronger in interactions between nonhub proteins. We also analysed the interaction types from a functional standpoint by using GO, which revealed that the interaction between disordered proteins frequently occurred in cellular processes, regulation, and metabolic processes. The number of interactions, especially in metabolic processes between disordered proteins, was 1.8 times as large as that in the randomised network. Another analysis conducted by using KEGG pathways provided results where several signaling pathways and disease-related pathways included many interactions between disordered proteins. All of these analyses suggest that human PPIs preferably occur between disordered proteins and that the flexibility of the interacting protein pairs may play an important role in human PPI networks.     

Coupled folding-binding in a hydrophobic/polar protein model: impact of synergistic folding and disordered flanks

Coupled folding-binding is central to the function of many intrinsically disordered proteins, yet not fully understood. With a continuous three-letter protein model, we explore the free-energy landscape of pairs of interacting sequences and how it is impacted by 1), variations in the binding mechanism; and 2), the addition of disordered flanks to the binding region. In particular, we focus on two sequences, one with 16 and one with 35 amino acids, which make a stable dimeric three-helix bundle at low temperatures. Three distinct binding mechanisms are realized by altering the stabilities of the individual monomers: docking, coupled folding-binding of a single α-helix, and synergistic folding and binding. Compared to docking, the free-energy barrier for binding is reduced when the single α-helix is allowed to fold upon binding, but only marginally. A greater reduction is found for synergistic folding, which in addition results in a binding transition state characterized by very few interchain contacts. Disordered flanking chain segments attached to the α-helix sequence can, despite a negligible impact on the dimer stability, lead to a downhill free-energy surface in which the barrier for binding is eliminated.     

Molecular principles of the interactions of disordered proteins

Thorough knowledge of the molecular principles of protein-protein recognition is essential to our understanding of protein function at the cellular level. Whereas interactions of ordered proteins have been analyzed in great detail, complexes of intrinsically unstructured/disordered proteins (IUPs) have hardly been addressed so far. Here, we have collected a database of 39 complexes of experimentally verified IUPs, and compared their interfaces with those of 72 complexes of ordered, globular proteins. The characteristic differences found between the two types of complexes suggest that IUPs represent a distinct molecular implementation of the principles of protein-protein recognition. The interfaces do not differ in size, but those of IUPs cover a much larger part of the surface of the protein than for their ordered counterparts. Moreover, IUP interfaces are significantly more hydrophobic relative to their overall amino acid composition, but also in absolute terms. They rely more on hydrophobic-hydrophobic than on polar-polar interactions. Their amino acids in the interface realize more intermolecular contacts, which suggests a better fit with the partner due to induced folding upon binding that results in a better adaptation to the partner. The two modes of interaction also differ in that IUPs usually use only a single continuous segment for partner binding, whereas the binding sites of ordered proteins are more segmented. Probably, all these features contribute to the increased evolutionary conservation of IUP interface residues. These noted molecular differences are also manifested in the interaction energies of IUPs. Our approximation of these by low-resolution force-fields shows that IUPs gain much more stabilization energy from intermolecular contacts, than from folding, i.e. they use their binding energy for folding. Overall, our findings provide a structural rationale to the prior suggestions that many IUPs are specialized for functions realized by protein-protein interactions.     

The C-terminal Tails of the Bacterial Chaperonin GroEL Stimulate Protein Folding by Directly Altering the Conformation of a Substrate Protein

Many essential cellular proteins fold only with the assistance of chaperonin machines like the GroEL-GroES system of Escherichia coli. However, the mechanistic details of assisted protein folding by GroEL-GroES remain the subject of ongoing debate. We previously demonstrated that GroEL-GroES enhances the productive folding of a kinetically trapped substrate protein through unfolding, where both binding energy and the energy of ATP hydrolysis are used to disrupt the inhibitory misfolded states. Here, we show that the intrinsically disordered yet highly conserved C-terminal sequence of the GroEL subunits directly contributes to substrate protein unfolding. Interactions between the C terminus and the non-native substrate protein alter the binding position of the substrate protein on the GroEL apical surface. The C-terminal tails also impact the conformational state of the substrate protein during capture and encapsulation on the GroEL ring. Importantly, removal of the C termini results in slower overall folding, reducing the fraction of the substrate protein that commits quickly to a productive folding pathway and slowing several kinetically distinct folding transitions that occur inside the GroEL-GroES cavity. The conserved C-terminal tails of GroEL are thus important for protein folding from the beginning to the end of the chaperonin reaction cycle.     

A Large Intrinsically Disordered Region in SKIP and Its Disorder-Order Transition Induced by PPIL1 Binding Revealed by NMR

Intrinsically disordered proteins or protein regions play an important role in fundamental biological processes. During spliceosome activation, a large structural rearrangement occurs. The Prp19 complex and related factors are involved in the catalytic activation of the spliceosome. Recent mass spectrometric analyses have shown that Ski interaction protein (SKIP) and peptidylprolyl isomerase-like protein 1 (PPIL1) are Prp19-related factors that constitute the spliceosome B, B*, and C complexes. Here, we report that a highly flexible region of SKIP (SKIPN, residues 59–129) is intrinsically disordered. Upon binding to PPIL1, SKIPN undergoes a disorder-order transition. A highly conserved fragment of SKIP (residues 59–79) called the PPIL1-binding fragment (PBF) was sufficient to bind PPIL1. The structure of PBF·PPIL1 complex, solved by NMR, shows that PBF exhibits an ordered structure and interacts with PPIL1 through electrostatic and hydrophobic interactions. Three subfragments in the PBF (residues 59–67, 68–73, and 74–79) show hook-like backbone structure, and interactions between these subfragments are necessary for PBF·PPIL1 complex formation. PPIL1 is a cyclophilin family protein. It is recruited by SKIP into the spliceosome by a region other than the peptidylprolyl isomerase active site. This enables the active site of PPIL1 to remain open in the complex and still function as a peptidylprolyl cis/trans-isomerase or molecular chaperon to facilitate the folding of other proteins in the spliceosomes. The large disordered region in SKIP provides an interaction platform. Its disorder-order transition, induced by PPIL1 binding, may adapt the requirement for a large structural rearrangement occurred in the activation of spliceosome.     

Topology-based modeling of intrinsically disordered proteins: balancing intrinsic folding and intermolecular interactions

Coupled binding and folding is frequently involved in specific recognition of so-called intrinsically disordered proteins (IDPs), a newly recognized class of proteins that rely on a lack of stable tertiary fold for function. Here, we exploit topology-based Gō-like modeling as an effective tool for the mechanism of IDP recognition within the theoretical framework of minimally frustrated energy landscape. Importantly, substantial differences exist between IDPs and globular proteins in both amino acid sequence and binding interface characteristics. We demonstrate that established Gō-like models designed for folded proteins tend to over-estimate the level of residual structures in unbound IDPs, whereas under-estimating the strength of intermolecular interactions. Such systematic biases have important consequences in the predicted mechanism of interaction. A strategy is proposed to recalibrate topology-derived models to balance intrinsic folding propensities and intermolecular interactions, based on experimental knowledge of the overall residual structure level and binding affinity. Applied to pKID/KIX, the calibrated Gō-like model predicts a dominant multistep sequential pathway for binding-induced folding of pKID that is initiated by KIX binding via the C-terminus in disordered conformations, followed by binding and folding of the rest of C-terminal helix and finally the N-terminal helix. This novel mechanism is consistent with key observations derived from a recent NMR titration and relaxation dispersion study and provides a molecular-level interpretation of kinetic rates derived from dispersion curve analysis. These case studies provide important insight into the applicability and potential pitfalls of topology-based modeling for studying IDP folding and interaction in general.     

Entropic Control of an Excited Folded-Like Conformation in a Disordered Protein Ensemble

Many intrinsically disordered proteins switch between unfolded and folded-like forms in the presence of their binding partner. The possibility of a pre-equilibrium between the two macrostates is challenging to discern given the complex conformational landscape. Here, we show that CytR, a disordered DNA-binding domain, samples a folded-like excited state in its native ensemble through equilibrium multi-probe spectroscopy, kinetics and an Ising-like statistical mechanical model. The population of the excited state increases upon stabilization of the native ensemble with an osmolyte, while decreasing with increasing temperatures. A conserved proline residue, the mutation of which weakens the binding affinity to the target promoter, is found to uniquely control the population of the minor excited state. Semi-quantitative statistical mechanical modeling reveals that the conformational diffusion coefficient of disordered CytR is an order of magnitude slower than the estimates from folded domains. The osmolyte and proline mutation smoothen and roughen up the landscape, respectively, apart from modulation of populations. Our work uncovers general strategies to probe for excited structured states in disordered ensembles, and to measure and modulate the roughness of the disordered landscapes, inter-conversion rates of species and their populations.     

Effect of Src kinase phosphorylation on disordered C-terminal domain of N-methyl-D-aspartic acid (NMDA) receptor subunit GluN2B protein

NMDA receptors are ligand-gated ion channels with a regulatory intracellular C-terminal domain (CTD). In GluN2B, the CTD is the largest domain in the protein but is intrinsically disordered. The GluN2B subunit is the major tyrosine-phosphorylated protein in synapses. Src kinase phosphorylates the GluN2B CTD, but it is unknown how this affects channel activity. In disordered proteins, phosphorylation can tip the balance between order and disorder. Transitions can occur in both directions, so it is not currently possible to predict the effects of phosphorylation. We used single molecule fluorescence to characterize the effects of Src phosphorylation on GluN2B. Scanning fluorescent labeling sites throughout the domain showed no positional dependence of the energy transfer. Instead, efficiency only scaled with the separation between labeling sites suggestive of a relatively featureless conformational energy landscape. Src phosphorylation led to a general expansion of the polypeptide, which would result in greater exposure of known protein-binding sites and increase the physical separation between contiguous sites. Phosphorylation makes the CTD more like a random coil leaving open the question of how Src exerts its effects on the NMDA receptor.     

Intrinsic Disorder in Protein Interactions: Insights From a Comprehensive Structural Analysis

We perform a large-scale study of intrinsically disordered regions in proteins and protein complexes using a non-redundant set of hundreds of different protein complexes. In accordance with the conventional view that folding and binding are coupled, in many of our cases the disorder-to-order transition occurs upon complex formation and can be localized to binding interfaces. Moreover, analysis of disorder in protein complexes depicts a significant fraction of intrinsically disordered regions, with up to one third of all residues being disordered. We find that the disorder in homodimers, especially in symmetrical homodimers, is significantly higher than in heterodimers and offer an explanation for this interesting phenomenon. We argue that the mechanisms of regulation of binding specificity through disordered regions in complexes can be as common as for unbound monomeric proteins. The fascinating diversity of roles of disordered regions in various biological processes and protein oligomeric forms shown in our study may be a subject of future endeavors in this area.     

Disordered Flanks Prevent Peptide Aggregation

Natively unstructured or disordered regions appear to be abundant in eukaryotic proteins. Many such regions have been found alongside small linear binding motifs. We report a Monte Carlo study that aims to elucidate the role of disordered regions adjacent to such binding motifs. The coarse-grained simulations show that small hydrophobic peptides without disordered flanks tend to aggregate under conditions where peptides embedded in unstructured peptide sequences are stable as monomers or as part of small micelle-like clusters. Surprisingly, the binding free energy of the motif is barely decreased by the presence of disordered flanking regions, although it is sensitive to the loss of entropy of the motif itself upon binding. This latter effect allows for reversible binding of the signalling motif to the substrate. The work provides insights into a mechanism that prevents the aggregation of signalling peptides, distinct from the general mechanism of protein folding, and provides a testable hypothesis to explain the abundance of disordered regions in proteins.     

Disordered Binding Regions and Linear Motifs—Bridging the Gap between Two Models of Molecular Recognition

Intrinsically disordered proteins (IDPs) exist without the presence of a stable tertiary structure in isolation. These proteins are often involved in molecular recognition processes via their disordered binding regions that can recognize partner molecules by undergoing a coupled folding and binding process. The specific properties of disordered binding regions give way to specific, yet transient interactions that enable IDPs to play central roles in signaling pathways and act as hubs of protein interaction networks. An alternative model of protein-protein interactions with largely overlapping functional properties is offered by the concept of linear interaction motifs. This approach focuses on distilling a short consensus sequence pattern from proteins with a common interaction partner. These motifs often reside in disordered regions and are considered to mediate the interaction roughly independent from the rest of the protein. Although a connection between linear motifs and disordered binding regions has been established through common examples, the complementary nature of the two concepts has yet to be fully explored. In many cases the sequence based definition of linear motifs and the structural context based definition of disordered binding regions describe two aspects of the same phenomenon. To gain insight into the connection between the two models, prediction methods were utilized. We combined the regular expression based prediction of linear motifs with the disordered binding region prediction method ANCHOR, each specialized for either model to get the best of both worlds. The thorough analysis of the overlap of the two methods offers a bioinformatics tool for more efficient binding site prediction that can serve a wide range of practical implications. At the same time it can also shed light on the theoretical connection between the two co-existing interaction models.