Design and Structure of an Equilibrium Protein Folding Intermediate: A Hint into Dynamical Regions of Proteins

Partly unfolded protein conformations close to the native state may play important roles in protein function and in protein misfolding. Structural analyses of such conformations which are essential for their fully physicochemical understanding are complicated by their characteristic low populations at equilibrium. We stabilize here with a single mutation the equilibrium intermediate of apoflavodoxin thermal unfolding and determine its solution structure by NMR. It consists of a large native region identical with that observed in the X-ray structure of the wild-type protein plus an unfolded region. Small-angle X-ray scattering analysis indicates that the calculated ensemble of structures is consistent with the actual degree of expansion of the intermediate. The unfolded region encompasses discontinuous sequence segments that cluster in the 3D structure of the native protein forming the FMN cofactor binding loops and the binding site of a variety of partner proteins. Analysis of the apoflavodoxin inner interfaces reveals that those becoming destabilized in the intermediate are more polar than other inner interfaces of the protein. Natively folded proteins contain hydrophobic cores formed by the packing of hydrophobic surfaces, while natively unfolded proteins are rich in polar residues. The structure of the apoflavodoxin thermal intermediate suggests that the regions of natively folded proteins that are easily responsive to thermal activation may contain cores of intermediate hydrophobicity.     

Fesselin is a natively unfolded protein

Fesselin is a heat stable proline-rich actin binding protein. The stability, amino acid composition, and ability to bind to several proteins suggested that fesselin may be unfolded under native conditions. While the complete sequence of fesselin is unknown an analysis of a closely related protein, synaptopodin 2 from Gallus gallus, indicates that fesselin consists of a series of unstructured regions interspersed between short folded regions. To determine if fesselin is natively unfolded, we compared fesselin to a known globular protein (myosin S1) and a known unfolded protein Cad22 (the COOH terminal 22 kDa fragment of caldesmon). Fesselin, and Cad22, had larger Stokes radii than globular proteins of equivalent mass. The environments of tryptophan residues of fesselin and Cad22 were the same in the presence and absence of 6 M guanidine hydrochloride. Fesselin had a circular dichroism spectrum that was primarily random coil. Changes in pH over the range of 1.5-11.5 did not alter that spectrum. Increasing the temperature to 85 degrees C caused an increase in the degree of secondary structure. Calmodulin binding to fesselin altered the environment of the tryptophan residues so that they became less sensitive to the quencher acrylamide. These results show that fesselin is a natively unfolded protein.     

Predictors of natively unfolded proteins: unanimous consensus score to detect a twilight zone between order and disorder in generic datasets

Background  

Natively unfolded proteins lack a well defined three dimensional structure but have important biological functions, suggesting a re-assignment of the structure-function paradigm. To assess that a given protein is natively unfolded requires laborious experimental investigations, then reliable sequence-only methods for predicting whether a sequence corresponds to a folded or to an unfolded protein are of interest in fundamental and applicative studies. Many proteins have amino acidic compositions compatible both with the folded and unfolded status, and belong to a twilight zone between order and disorder. This makes difficult a dichotomic classification of protein sequences into folded and natively unfolded ones. In this work we propose an operational method to identify proteins belonging to the twilight zone by combining into a consensus score good performing single predictors of folding.     

The N-terminal domain of p53 is natively unfolded

p53 is one of the key molecules regulating cell proliferation, apoptosis and tumor suppression by integrating a wide variety of signals. The structural basis for this function is still poorly understood. p53 appears to exercise its function as a modular protein in which different functions are associated with distinct domains. Presumably, p53 contains both folded and partially structured parts. Here, we have investigated the structure of the isolated N-terminal part of p53 (amino acid residues 1-93) using biophysical techniques. We demonstrate that this domain is devoid of tertiary structure and largely missing secondary structure elements. It exhibits a large hydrodynamic radius, typical for unfolded proteins. These findings suggest strongly that the entire N-terminal part of p53 is natively unfolded under physiological conditions. Furthermore, the binding affinity to its functional antagonist Mdm2 was investigated. A comparison of the binding of human Mdm2 to the N-terminal part of p53 and full-length p53 suggests that unfolded and folded parts of p53 function synergistically.     

Prediction of natively unfolded regions in protein chain

We have shown that the ability of a protein to be in globular or in natively unfolded state (under native conditions) may be determined (besides low overall hydrophobicity and a large net charge) by such a property as the average environment density, the average number of residues enclosed at the given distance. A statistical scale of the average number of residues enclosed at the given distance for 20 types of amino acid residues in globular state has been created on the basis of 6626 protein structures. Using this scale for separation of 80 globular and 90 natively unfolded proteins we fail only in 11% of proteins (compared with 17% of errors which are observed if to use hydrophobicity scale). The present scale may be used both for prediction of form (folded or unfolded) of the native state of protein and for prediction of natively unfolded regions in protein chains. The results of comparison of our method of predicting natively unfolded regions with the other known methods show that our method has the highest fraction of correctly predicted natively unfolded regions (that is 87% and 77% if to make averaging over residues and over proteins correspondingly).     

To be folded or to be unfolded?

The lack of ordered structure in “natively unfolded” proteins raises a general question: Are there intrinsic properties of amino acid residues that are responsible for the absence of fixed structure at physiological conditions? In this article, we demonstrate that the competence of a protein to be folded or to be unfolded may be determined by the property of amino acid residues to form a sufficient number of contacts in a globular state. The expected average number of contacts per residue calculated from the amino acid sequence alone (using the average number of contacts for 20 amino acid residues in globular proteins) can be used as one of the simple indicators of natively unfolded proteins. The prediction accuracy for the sets of 80 folded and 90 natively unfolded proteins reaches 89% if the expected average number of contacts is used as a parameter and 83% in the case of hydrophobicity. An optimal set of artificial parameters for 20 amino acid residues obtained by Monte Carlo algorithm to maximally separate the sets of 90 natively unfolded and 80 folded proteins demonstrates the upper limit for prediction accuracy, which is 95%.     

Residual structure in disordered peptides and unfolded proteins from multivariate analysis and ab initio simulation of Raman optical activity data

Vibrational Raman optical activity (ROA), measured as a small difference in the intensity of Raman scattering from chiral molecules in right- and left-circularly polarized incident light, or as the intensity of a small circularly polarized component in the scattered light, is a powerful probe of the aqueous solution structure of proteins. The large number of structure-sensitive bands in protein ROA spectra makes multivariate analysis techniques such as nonlinear mapping (NLM) especially favorable for determining structural relationships between different proteins. We have previously used NLM to map a large dataset of peptide, protein, and virus ROA spectra into a readily visualizable two-dimensional space in which points close to or distant from each other, respectively, represent similar or dissimilar structures. As well as folded proteins, our dataset contains ROA spectra from many natively unfolded proteins, proteins containing both folded and unfolded domains, denatured partially structured molten globule and reduced protein states, together with folded proteins containing little or no alpha-helix or beta-sheet. In this article, the relative positions of these systems in the NLM plot are used to obtain information about any residual structure that they may contain. The striking differences between the structural propensities of proteins that are unfolded in their native states and those that are unfolded due to denaturation may be responsible for their often very different behavior, especially with regard to aggregation. An ab initio simulation of the Raman and ROA spectra of an alanine oligopeptide in the poly(L-proline) II-helical conformation confirms previous suggestions that this conformation is a significant structural element in disordered peptides and natively unfolded proteins. The use of ROA to identify and characterize proteins containing significant amounts of unfolded structure will, inter alia, be valuable in structural genomics/proteomics since unfolded sequences often inhibit crystallization.     

Why are "natively unfolded" proteins unstructured under physiologic conditions?

"Natively unfolded" proteins occupy a unique niche within the protein kingdom in that they lack ordered structure under conditions of neutral pH in vitro. Analysis of amino acid sequences, based on the normalized net charge and mean hydrophobicity, has been applied to two sets of proteins: small globular folded proteins and "natively unfolded" ones. The results show that "natively unfolded" proteins are specifically localized within a unique region of charge-hydrophobicity phase space and indicate that a combination of low overall hydrophobicity and large net charge represent a unique structural feature of "natively unfolded" proteins.     

Defining the specificity space of the human SRC homology 2 domain

Src homology 2 (SH2) domains are the largest family of interaction modules encoded by the human genome to recognize tyrosine-phosphorylated sequences and thereby play pivotal roles in transducing and controlling cellular signals emanating from protein-tyrosine kinases. Different SH2 domains select for distinct phosphopeptides, and the function of a given SH2 domain is often dictated by the specific motifs that it recognizes. Therefore, deciphering the phosphotyrosyl peptide motif recognized by an SH2 domain is the key to understanding its cellular function. Here we cloned all 120 SH2 domains identified in the human genome and determined the phosphotyrosyl peptide binding properties of 76 SH2 domains by screening an oriented peptide array library. Of these 76, we defined the selectivity for 43 SH2 domains and refined the binding motifs for another 33 SH2 domains. We identified a number of novel binding motifs, which are exemplified by the BRDG1 SH2 domain that selects specifically for a bulky, hydrophobic residue at P + 4 relative to the Tyr(P) residue. Based on the oriented peptide array library data, we developed scoring matrix-assisted ligand identification (or SMALI), a Web-based program for predicting binding partners for SH2-containing proteins. When applied to SH2D1A/SAP (SLAM-associated protein), a protein whose mutation or deletion underlies the X-linked lymphoproliferative syndrome, SMALI not only recapitulated known interactions but also identified a number of novel interacting proteins for this disease-associated protein. SMALI also identified a number of potential interactors for BRDG1, a protein whose function is largely unknown. Peptide in-solution binding analysis demonstrated that a SMALI score correlates well with the binding energy of a peptide to a given SH2 domain. The definition of the specificity space of the human SH2 domain provides both the necessary molecular basis and a platform for future exploration of the functions for SH2-containing proteins in cells.     

Intrinsic unstructuredness and abundance of PEST motifs in eukaryotic proteomes

The study of unfolded protein regions has gained importance because of their prevalence and important roles in various cellular functions. These regions have characteristically high net charge and low hydrophobicity. The amino acid sequence determines the intrinsic unstructuredness of a region and, therefore, efforts are ongoing to delineate the sequence motifs, which might contribute to protein disorder. We find that PEST motifs are enriched in the characterized disordered regions as compared with globular ones. Analysis of representative PDB chains revealed very few structures containing PEST sequences and the majority of them lacked regular secondary structure. A proteome-wide study in completely sequenced eukaryotes with predicted unfolded and folded proteins shows that PEST proteins make up a large fraction of unfolded dataset as compared with the folded proteins. Our data also reveal the prevalence of PEST proteins in eukaryotic proteomes (approximately 25%). Functional classification of the PEST-containing proteins shows an over- and under-representation in proteins involved in regulation and metabolism, respectively. Furthermore, our analysis shows that predicted PEST regions do not exhibit any preference to be localized in the C terminals of proteins, as reported earlier.     

Limited proteolysis of natively unfolded protein 4E‐BP1 in the presence of trifluoroethanol

Natively unfolded (intrinsically disordered (ID) proteins) have been attracting an increasing attention due to their involvement in many regulatory processes. Natively unfolded proteins can fold upon binding to their metabolic partners. Coupled folding and binding events usually involve only relatively short motifs (binding motifs). These binding motifs which are able to fold should have an increased propensity to form a secondary structure. The aim of the present work was to probe the conformation of the intrinsically disordered protein 4E‐BP1 in the native and partly folded states by limited proteolysis and to reveal regions with a high propensity to form an ordered structure. Trifuoroethanol (TFE) in low concentrations (up to 15 vol%) was applied to increase the helical population of protein regions with a high intrinsic propensity to fold. When forming helical structures, these regions lose mobility and become more protected from proteases than random/unfolded protein regions. Limited proteolysis followed by mass spectrometry analysis allows identification of the regions with decreased mobility in TFE solutions. Trypsin and V8 proteases were used to perform limited proteolysis of the 4E‐BP1 protein in buffer and in solutions with low TFE concentrations at 37°C and at elevated temperatures (42 and 50°C). Comparison of the results obtained with the previously established 4E‐BP1 structure and the binding motif illustrates the ability of limited proteolysis in the presence of a folding assistant (TFE) to map the regions with high and low propensities to form a secondary structure revealing potential binding motifs inside the intrinsically disordered protein. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 591–602, 2014.     

Zn(2+)-mediated structure formation and compaction of the "natively unfolded" human prothymosin alpha

Human recombinant prothymosin alpha (ProTalpha) is known to have coil-like conformation at neutral pH; i.e., it belongs to the class of "natively unfolded" proteins. By means of circular dichroism, SAXS, and ANS fluorescence, we have investigated the effect of several divalent cations on the structure of this protein. Results of these studies are consistent with the conclusion that ProTalpha conformation is unaffected by large excess of Ca(2+), Mg(2+), Mn(2+), Cu(2+), and Ni(2+). However, Zn(2+) induces compaction and considerable rearrangement of the protein structure. This means that ProTalpha can specifically interact with Zn(2+) (K(D) approximately 10(-3) M), and such interactions induce folding of the natively unfolded protein into a compact partially folded (premolten globule-like) conformation. It is possible that these structural changes may be important for the function of this protein.     

Prediction of natively unfolded regions in protein chains

Analysis showed that the globular or natively unfolded state of a protein can be inferred not only from a lower hydrophobicity or a higher charge, but also from the average environment density (average number of close residues located within a certain distance of a given one) of its residues. A database of 6626 protein structures was used to construct a statistical scale of the average number of close residues in globular structures for the 20 amino acids. The portion of false predictions in distinguishing between 80 globular and 90 natively unfolded proteins was 11% with the new scale and 17% with a hydrophobicity scale. The new scale proved suitable for predicting the folded or unfolded state for native proteins or the natively unfolded regions for protein chains. In comparisons with the available algorithms, the new method yielded the highest portion of true predictions (87 and 77% with averaging over residues and over proteins, respectively).     

SAP couples Fyn to SLAM immune receptors

SAP (SLAM-associated protein) is a small lymphocyte-specific signalling molecule that is defective or absent in patients with X-linked lymphoproliferative syndrome (XLP). Consistent with its single src homology 2 (SH2) domain architecture and unusually high affinity for SLAM (also called CD150), SAP has been suggested to function by blocking binding of SHP-2 or other SH2-containing signalling proteins to SLAM receptors. Additionally, SAP has recently been shown to be required for recruitment and activation of the Src-family kinase FynT after SLAM ligation. This signalling 'adaptor' function has been difficult to conceptualize, because unlike typical SH2-adaptor proteins, SAP contains only a single SH2 domain and lacks other recognized protein interaction domains or motifs. Here, we show that the SAP SH2 domain binds to the SH3 domain of FynT and directly couples FynT to SLAM. The crystal structure of a ternary SLAM-SAP-Fyn-SH3 complex reveals that SAP binds the FynT SH3 domain through a surface-surface interaction that does not involve canonical SH3 or SH2 binding interactions. The observed mode of binding to the Fyn-SH3 domain is expected to preclude the auto-inhibited conformation of Fyn, thereby promoting activation of the kinase after recruitment. These findings broaden our understanding of the functional repertoire of SH3 and SH2 domains.     

Prediction of "hot spots" of aggregation in disease-linked polypeptides

Background  

The polypeptides involved in amyloidogenesis may be globular proteins with a defined 3D-structure or natively unfolded proteins. The first class includes polypeptides such as β2-microglobulin, lysozyme, transthyretin or the prion protein, whereas β-amyloid peptide, amylin or α-synuclein all belong to the second class. Recent studies suggest that specific regions in the proteins act as "hot spots" driving aggregation. This should be especially relevant for natively unfolded proteins or unfolded states of globular proteins as they lack significant secondary and tertiary structure and specific intra-chain interactions that can mask these aggregation-prone regions. Prediction of such sequence stretches is important since they are potential therapeutic targets.     

Quantitative relation between intermolecular and intramolecular binding of pro-rich peptides to SH3 domains

Flexible linkers are often found to tether binding sequence motifs or connect protein domains. Here we analyze three usages of flexible linkers: 1), intramolecular binding of proline-rich peptides (PRPs) to SH3 domains for kinase regulation; 2), intramolecular binding of PRP for increasing the folding stability of SH3 domains; and 3), covalent linking of PRPs and other ligands for high-affinity bivalent binding. The basis of these analyses is a quantitative relation between intermolecular and intramolecular binding constants. This relation has the form K(i) = K(e0)p for intramolecular binding and K(e) = K(e01)K(e02)p for bivalent binding. The effective concentration p depends on the length of the linker and the distance between the linker attachment points in the bound state. Several applications illustrate the usefulness of the quantitative relation. These include intramolecular binding to the Itk SH3 domain by an internal PRP and to a circular permutant of the alpha-spectrin SH3 domain by a designed PRP, and bivalent binding to the two SH3 domains of Grb2 by two linked PRPs. These and other examples suggest that flexible linkers and sequence motifs tethered to them, like folded protein domains, are also subject to tight control during evolution.     

A novel pro-Arg motif recognized by WW domains

WW domains mediate protein-protein interactions through binding to short proline-rich sequences. Two distinct sequence motifs, PPXY and PPLP, are recognized by different classes of WW domains, and another class binds to phospho-Ser-Pro sequences. We now describe a novel Pro-Arg sequence motif recognized by a different class of WW domains using data from oriented peptide library screening, expression cloning, and in vitro binding experiments. The prototype member of this group is the WW domain of formin-binding protein 30 (FBP30), a p53-regulated molecule whose WW domains bind to Pro-Arg-rich cellular proteins. This new Pro-Arg sequence motif re-classifies the organization of WW domains based on ligand specificity, and the Pro-Arg class now includes the WW domains of FBP21 and FE65. A structural model is presented which rationalizes the distinct motifs selected by the WW domains of YAP, Pin1, and FBP30. The Pro-Arg motif identified for WW domains often overlaps with SH3 domain motifs within protein sequences, suggesting that the same extended proline-rich sequence could form discrete SH3 or WW domain complexes to transduce distinct cellular signals.     

p25alpha is flexible but natively folded and binds tubulin with oligomeric stoichiometry

p25alpha is a 219-residue protein which stimulates aberrant tubulin polymerization and is implicated in a variety of other functions. The protein has unusual secondary structure involving significant amounts of random coil, and binding to microtubules is accompanied by a large structural change, suggesting a high degree of plasticity. p25alpha has been proposed to be natively unfolded, so that folding is coupled to interaction with its physiological partners. Here we show that recombinant human p25alpha is folded under physiological conditions, since it has a well structured and solvent-sequestered aromatic environment and considerable chemical shift dispersion of amide and aliphatic protons. With increasing urea concentrations, p25alpha undergoes clear spectral changes suggesting significant loss of structure. p25alpha unfolds cooperatively in urea according to a simple two-state transition with a stability in water of approximately 5 kcal/mol. The protein behaves as a monomer and refolds with a transient on-pathway folding intermediate. However, high sensitivity to proteolytic attack and abnormal gel filtration migration behavior suggests a relatively extended structure, possibly organized in distinct domains. A deletion mutant of p25alpha lacking residues 3-43 also unfolds cooperatively and with similar stability, suggesting that the N-terminal region is largely unstructured. Both proteins undergo significant loss of structure when bound to monomeric tubulin. The stoichiometry of binding is estimated to be 3-4 molecules of tubulin per p25alpha and is not significantly affected by the deletion of residues 3-43. In conclusion, we dismiss the proposal that p25alpha is natively unfolded, although the protein is relatively flexible. This flexibility may be linked to its tubulin-binding properties.     

Lipid-binding hSH3 domains in immune cell adapter proteins

SH3 domains represent versatile scaffolds within eukaryotic cells by targeting proline-rich sequences within intracellular proteins. More recently, binding of SH3 domains to unusual peptide motifs, folded proteins or lipids has been reported. Here we show that the newly defined hSH3 domains of immune cell adapter proteins bind lipid membranes with distinct affinities. The interaction of the hSH3 domains of adhesion and degranulation promoting adapter protein (ADAP) and PRAM-1 (Promyelocytic-Retinoic acid receptor alpha target gene encoding an Adaptor Molecule-1), with phosphatidylcholine-containing liposomes is observed upon incorporation of phosphatidylserine (PS) or phosphoinositides (PIs) into the membrane bilayer. Mechanistically we show that stable association of the N-terminal, amphipathic helix with the beta-sheet scaffold favours lipid binding and that the interaction with PI(4,5)P(2)-containing liposomes is consistent with a single-site, non-cooperative binding mechanism. Functional investigations indicate that deletion of both amphipathic helices of the hSH3 domains reduces the ability of ADAP to enhance adhesion and migration in stimulated T cells.     

The structure and interactions of the proline-rich domain of ASPP2

ASPP2 is a pro-apoptotic protein that stimulates the p53-mediated apoptotic response. The C terminus of ASPP2 contains ankyrin (Ank) repeats and a SH3 domain, which mediate its interactions with numerous partner proteins such as p53, NFkappaB, and Bcl-2. It also contains a proline-rich domain (ASPP2 Pro), whose structure and function are unclear. Here we used biophysical and biochemical methods to study the structure and the interactions of ASPP2 Pro, to gain insight into its biological role. We show, using biophysical and computational methods, that the ASPP2 Pro domain is natively unfolded. We found that the ASPP2 Pro domain interacts with the ASPP2 Ank-SH3 domains, and mapped the interaction sites in both domains. Using a combination of peptide array screening, biophysical and biochemical techniques, we found that ASPP2 Ank-SH3, but not ASPP2 Pro, mediates interactions of ASPP2 with peptides derived from its partner proteins. ASPP2 Pro-Ank-SH3 bound a peptide derived from its partner protein NFkappaB weaker than ASPP2 Ank-SH3 bound this peptide. This suggested that the presence of the proline-rich domain inhibited the interactions mediated by the Ank-SH3 domains. Furthermore, a peptide from ASPP2 Pro competed with a peptide derived from NFkappaB on binding to ASPP2 Ank-SH3. Based on our results, we propose a model in which the interaction between the ASPP2 domains regulates the intermolecular interactions of ASPP2 with its partner proteins.