Exploring the impact of polyproline II (PII) conformational bias on the binding of peptides to the SEM-5 SH3 domain

Beta turn propensity and a model polymer scaling exponent identify intrinsically disordered phase-separating proteins

The complex cellular milieu can spontaneously demix, or phase separate, in a process controlled in part by intrinsically disordered (ID) proteins. A protein's propensity to phase separate is thought to be driven by a preference for protein-protein over protein-solvent interactions. The hydrodynamic size of monomeric proteins, as quantified by the polymer scaling exponent (v), is driven by a similar balance. We hypothesized that mean v, as predicted by protein sequence, would be smaller for proteins with a strong propensity to phase separate. To test this hypothesis, we analyzed protein databases containing subsets of proteins that are folded, disordered, or disordered and known to spontaneously phase separate. We find that the phase-separating disordered proteins, on average, had lower calculated values of v compared with their non-phase-separating counterparts. Moreover, these proteins had a higher sequence-predicted propensity for β-turns. Using a simple, surface area-based model, we propose a physical mechanism for this difference: transient β-turn structures reduce the desolvation penalty of forming a protein-rich phase and increase exposure of atoms involved in π/sp2 valence electron interactions. By this mechanism, β-turns could act as energetically favored nucleation points, which may explain the increased propensity for turns in ID regions (IDRs) utilized biologically for phase separation. Phase-separating IDRs, non-phase-separating IDRs, and folded regions could be distinguished by combining v and β-turn propensity. Finally, we propose a new algorithm, ParSe (partition sequence), for predicting phase-separating protein regions, and which is able to accurately identify folded, disordered, and phase-separating protein regions based on the primary sequence.

Conserved allosteric ensembles in disordered proteins using TROSY/anti-TROSY R2-filtered spectroscopy

Defining the role of intrinsic disorder in proteins in the myriad of biological processes with which it is involved represents a significant goal in modern biophysics. Toward this end, NMR is uniquely suited for molecular studies of dynamic and disordered regions, but studying these regions in concert with their more structured domains and binding partners presents spectroscopic challenges. Here, we investigate the interactions between the structured and disordered regions of the human glucocorticoid receptor (GR). To do this, we developed an NMR strategy that relies on a novel relaxation filter for the simultaneous study of structured and unstructured regions. Using this approach, we conducted a comparative analysis of three translational isoforms of GR containing a folded DNA-binding domain (DBD) and two disordered regions that flank the DBD, one of which varies in size in the different isoforms. Notably, we were able to assign resonances that had previously been inaccessible because of the spectral complexity of the translational isoforms, which in turn allowed us to 1) identify a region of the structured DBD that undergoes significant changes in the local chemical environment in the presence of the disordered region and 2) determine differences in the conformational ensembles of the disordered regions of the translational isoforms. Furthermore, an ensemble-based thermodynamic analysis of the isoforms reveals conserved patterns of stability within the N-terminal domain of GR that persist despite low sequence conservation. These studies provide an avenue for further investigations of the mechanistic underpinnings of the functional relevance of the translational isoforms of GR while also providing a general NMR strategy for studying systems containing both structured and disordered regions.

Structural and Energetic Characterization of the Denatured State from the Perspectives of Peptides, the Coil Library, and Intrinsically Disordered Proteins

The α and polyproline II (PPII) basins are the two most populated regions of the Ramachandran map when constructed from the protein coil library, a widely used denatured state model built from the segments of irregular structure found in the Protein Data Bank. This indicates the α and PPII conformations are dominant components of the ensembles of denatured structures that exist in solution for biological proteins, an observation supported in part by structural studies of short, and thus unfolded, peptides. Although intrinsic conformational propensities have been determined experimentally for the common amino acids in short peptides, and estimated from surveys of the protein coil library, the ability of these intrinsic conformational propensities to quantitatively reproduce structural behavior in intrinsically disordered proteins (IDPs), an increasingly important class of proteins in cell function, has thus far proven elusive to establish. Recently, we demonstrated that the sequence dependence of the mean hydrodynamic size of IDPs in water and the impact of heat on the coil dimensions, provide access to both the sequence dependence and thermodynamic energies that are associated with biases for the α and PPII backbone conformations. Here, we compare results from peptide-based studies of intrinsic conformational propensities and surveys of the protein coil library to those of the sequence-based analysis of heat effects on IDP hydrodynamic size, showing that a common structural and thermodynamic description of the protein denatured state is obtained.

Impact of Heat on Coil Hydrodynamic Size Yields the Energetics of Denatured State Conformational Bias

Conformational equilibria in the protein denatured state have key roles regulating folding, stability, and function. The extent of conformational bias in the protein denatured state under folding conditions, however, has thus far proven elusive to quantify, particularly with regard to its sequence dependence and energetic character. To better understand the structural preferences of the denatured state, we analyzed both the sequence dependence to the mean hydrodynamic size of disordered proteins in water and the impact of heat on the coil dimensions, showing that the sequence dependence and thermodynamic energies associated with intrinsic biases for the α and polyproline II (PPII) backbone conformations can be obtained. Experiments that evaluate how the hydrodynamic size changes with compositional changes in the protein reveal amino acid specific preferences for PPII that are in good quantitative agreement with calorimetry-measured values from unfolded peptides and those inferred by survey of the protein coil library. At temperatures above 25 °C, the denatured state follows the predictions of a PPII-dominant ensemble. Heat effects on coil hydrodynamic size indicate the α bias is comparable to the PPII bias at cold temperatures. Though historically thought to give poor resolution to structural details, the hydrodynamic size of the unfolded state is found to be an effective reporter on the extent of the biases for the α and PPII backbone conformations.

Sequence Reversal Prevents Chain Collapse and Yields Heat-Sensitive Intrinsic Disorder

Sequence patterns of charge, hydrophobicity, hydrogen bonding, and other amino acid physicochemical properties contribute to mechanisms of protein folding, but how sequence composition and patterns influence the conformational dynamics of the denatured state ensemble is not fully understood. To investigate structure-sequence relationships in the denatured state, we reversed the sequence of staphylococcal nuclease and characterized its structure, thermodynamic character, and hydrodynamic radius using circular dichroism spectroscopy, dynamic light scattering, analytical ultracentrifugation, and size-exclusion chromatography as a function of temperature. The macromolecular size of "Retro-nuclease" is highly expanded in solution with characteristics similar to biological intrinsically disordered proteins. In contradistinction to a disordered state, Retro-nuclease exhibits a broad sigmoid transition of its hydrodynamic dimensions as temperature is increased, indicating a thermodynamically controlled compaction. Counterintuitively, the magnitude of these temperature-induced hydrodynamic changes exceed that observed from thermal denaturation of folded unaltered staphylococcal nuclease. Undetectable by calorimetry and intrinsic tryptophan fluorescence, the lack of heat capacity or fluorescence changes throughout the thermal transition indicate canonical hydrophobic collapse did not drive the Retro-nuclease structural transitions. Temperature-dependent circular dichroism spectroscopy performed on Retro-nuclease and computer simulations correlate to temperature sensitivity in the intrinsic sampling of backbone conformations for polyproline II and α-helix. The experimental results indicate a role for sequence direction in mediating the collapse of the polypeptide chain, whereas the simulation trends illustrate the generality of the observed heat effects on disordered protein structure.

Intrinsic α helix propensities compact hydrodynamic radii in intrinsically disordered proteins

Proteins that lack tertiary stability under normal conditions, known as intrinsically disordered, exhibit a wide range of biological activities. Molecular descriptions for the biology of intrinsically disordered proteins (IDPs) consequently rely on disordered structural models, which in turn require experiments that assess the origins to structural features observed. For example, while hydrodynamic size is mostly insensitive to sequence composition in chemically denatured proteins, IDPs show strong sequence-specific effects in the hydrodynamic radius (R_h ) when measured under normal conditions. To investigate sequence-modulation of IDP R_h , disordered ensembles generated by a hard sphere collision model modified with a structure-based parameterization of the solution energetics were used to parse the contributions of net charge, main chain dihedral angle bias, and excluded volume on hydrodynamic size. Ensembles for polypeptides 10-35 residues in length were then used to establish power-law scaling relationships for comparison to experimental R_h from 26 IDPs. Results showed the expected outcomes of increased hydrodynamic size from increases in excluded volume and net charge, and compaction from chain-solvent interactions. Chain bias representing intrinsic preferences for α helix and polyproline II (PP_II ), however, modulated R_h with intricate dependence on the simulated propensities. PP_II propensities at levels expected in IDPs correlated with heightened R_h sensitivity to even weak α helix propensities, indicating bias for common (φ, ψ) are important determinants of hydrodynamic size. Moreover, data show that IDP R_h can be predicted from sequence with good accuracy from a small set of physicochemical properties, namely intrinsic conformational propensities and net charge. Proteins 2017; 85:296-311. © 2016 Wiley Periodicals, Inc.

Phosphorylation Increases Persistence Length and End-to-End Distance of a Segment of Tau Protein

Intrinsically disordered regions of proteins, which lack unique tertiary structure under physiological conditions, are enriched in phosphorylation sites and in significant local bias toward the polyproline II conformation. The overrepresented coincidence of this posttranslational regulatory signal and local conformational bias within unstructured regions raises a question: can phosphorylation serve to manipulate the conformational preferences of a disordered protein? In this study, we use time-resolved fluorescence resonance energy transfer and a, to our knowledge, novel data analysis method to directly measure the end-to-end distance distribution of a phosphorylatable peptide derived from the human microtubule associated protein tau. Our results show that phosphorylation at threonine or serine extends the end-to-end distance and increases the effective persistence length of the tested model peptides. Unexpectedly, the extension is independent of salt concentration, suggestive of a nonelectrostatic origin. The phosphorylation extension and stiffening effect provides a peptide-scale physical interpretation for the posttranslational regulation of the highly abundant protein-protein interactions found in disordered proteins, as well as a potential insight into the regulatory mechanism of the tau protein's microtubule binding activity.

Hydrodynamic Radii of Intrinsically Disordered Proteins Determined from Experimental Polyproline II Propensities

The properties of disordered proteins are thought to depend on intrinsic conformational propensities for polyproline II (PP_II) structure. While intrinsic PP_II propensities have been measured for the common biological amino acids in short peptides, the ability of these experimentally determined propensities to quantitatively reproduce structural behavior in intrinsically disordered proteins (IDPs) has not been established. Presented here are results from molecular simulations of disordered proteins showing that the hydrodynamic radius (R_h) can be predicted from experimental PP_II propensities with good agreement, even when charge-based considerations are omitted. The simulations demonstrate that R_h and chain propensity for PP_II structure are linked via a simple power-law scaling relationship, which was tested using the experimental R_h of 22 IDPs covering a wide range of peptide lengths, net charge, and sequence composition. Charge effects on R_h were found to be generally weak when compared to PP_II effects on R_h. Results from this study indicate that the hydrodynamic dimensions of IDPs are evidence of considerable sequence-dependent backbone propensities for PP_II structure that qualitatively, if not quantitatively, match conformational propensities measured in peptides.

Molecular models of disordered protein structures are needed to elucidate the functional mechanisms of intrinsically disordered proteins, a class of proteins implicated in many disease pathologies and human health issues. Several studies have measured intrinsic conformational propensities for polyproline II helix, a key structural motif of disordered proteins, in short peptides. Whether or not these experimental polyproline II propensities, which vary by amino acid type, reproduce structural behavior in intrinsically disordered proteins has yet to be demonstrated. Presented here are simulation results showing that polyproline II propensities from short peptides accurately describe sequence-dependent variability in the hydrodynamic dimensions of intrinsically disordered proteins. Good agreement was observed from a simple molecular model even when charge-based considerations were ignored, predicting that global organization of disordered protein structure is strongly dependent on intrinsic conformational propensities and, for many intrinsically disordered proteins, modulated only weakly by coulombic effects.

Neighboring residue effects in terminally blocked dipeptides: implications for residual secondary structures in intrinsically unfolded/disordered proteins

For nuclear magnetic resonance (NMR)-based protein structure determinations, the random coil chemical shifts are very important because the secondary and tertiary protein structure predictions become possible by examining deviations of measured chemical shifts from those reference chemical shift values. In addition, neighboring residue effects on chemical shifts and J-coupling constants are crucial in understanding the nature of conformational propensities exhibited by unfolded or intrinsically disordered proteins. We recently reported the 1D NMR results for a complete set of terminally blocked dipeptides (Oh KI, Jung YS, Hwang GS, Cho M. J Biomol NMR 2012;53:25-41), but the NMR resonance assignments were not possible so that the average chemical shifts and J-coupling constants were only considered. In the present work, to thoroughly investigate the neighboring residue effects and random coil chemical shifts we extend the previous studies with 2D NMR, and measured all the (3) J(HNHα) values and H(α) and H(N) chemical shifts of the same set of terminally blocked dipeptides that are free from structural effects like secondary structure, hydrogen-bond, long-range backbone, and side-chain interactions. In particular, the preceding and following residue effects on amino-acid backbone conformational propensities are revealed and directly compared with previous works on either short peptides or empirical chemical shift database.

Temperature effects on the hydrodynamic radius of the intrinsically disordered N-terminal region of the p53 protein

Intrinsically disordered proteins (IDPs) are often characterized in terms of the hydrodynamic radius, Rh . The Rh of IDPs are known to depend on fractional proline content and net charge, where increased numbers of proline residues and increased net charge cause larger Rh . Though sequence and charge effects on the Rh of IDPs have been studied, the temperature sensitivity has been noted only briefly. Reported here are Rh measurements in the temperature range of 5-75°C for the intrinsically disordered N-terminal region of the p53 protein, p53(1-93). Of note, the Rh of this protein fragment was highly sensitive to temperature, decreasing from 35 Å at 5°C to 26 Å at 75°C. Computer generated simulations of conformationally dynamic and disordered polypeptide chains were performed to provide a hypothesis for the heat-induced compaction of p53(1-93) structure, which was opposite to the heat-induced increase in Rh observed for a model folded protein. The simulations demonstrated that heat caused Rh to trend toward statistical coil values for both proteins, indicating that the effects of heat on p53(1-93) structure could be interpreted as thermal denaturation. The simulation data also predicted that proline content contributed minimally to the native Rh of p53(1-93), which was confirmed by measuring Rh for a substitution variant that had all 22 proline residues changed for glycine.

Evolutionary conservation of the polyproline II conformation surrounding intrinsically disordered phosphorylation sites

Intrinsically disordered (ID) proteins function in the absence of a unique stable structure and appear to challenge the classic structure-function paradigm. The extent to which ID proteins take advantage of subtle conformational biases to perform functions, and whether signals for such mechanism can be identified in proteome-wide studies is not well understood. Of particular interest is the polyproline II (PII) conformation, suggested to be highly populated in unfolded proteins. We experimentally determine a complete calorimetric propensity scale for the PII conformation. Projection of the scale into representative eukaryotic proteomes reveals significant PII bias in regions coding for ID proteins. Importantly, enrichment of PII in ID proteins, or protein segments, is also captured by other PII scales, indicating that this enrichment is robustly encoded and universally detectable regardless of the method of PII propensity determination. Gene ontology (GO) terms obtained using our PII scale and other scales demonstrate a consensus for molecular functions performed by high PII proteins across the proteome. Perhaps the most striking result of the GO analysis is conserved enrichment (P < 10(-8) ) of phosphorylation sites in high PII regions found by all PII scales. Subsequent conformational analysis reveals a phosphorylation-dependent modulation of PII, suggestive of a conserved "tunability" within these regions. In summary, the application of an experimentally determined polyproline II (PII) propensity scale to proteome-wide sequence analysis and gene ontology reveals an enrichment of PII bias near disordered phosphorylation sites that is conserved throughout eukaryotes.

Temperature and urea have opposing impacts on polyproline II conformational bias

The native states of globular proteins have been accessed in atomic detail by X-ray crystallography and nuclear magnetic resonance spectroscopy, yet characterization of denatured proteins beyond global metrics has proven to be elusive. Denatured proteins have been observed to exhibit global geometric properties of a random coil polymer. However, this does not preclude the existence of nonrandom, local conformational bias that may be significant for protein folding and function. Indeed, circular dichroism (CD) spectroscopy and other methods have suggested that the denatured state contains considerable local bias to the polyproline II (PII) conformation. Here, we develop predictive models to determine the extent that temperature and the chemical denaturant urea modulate PII propensity. In agreement with our predictive model, PII propensity is observed experimentally to decrease with an increase in temperature. Conversely, urea appears to promote the PII conformation as determined by CD and isothermal titration calorimetry. Importantly, the calorimetric data are in quantitative agreement with a model that predicts the stability of the PII helix relative to other denatured state conformations based upon solvent accessible surface area and experimentally measured Gibbs transfer free energies. The ability of urea to promote the PII conformation can be attributed to the favorable interaction of urea with the peptide backbone. Thus, perturbing denatured states by temperature or cosolutes has subtle, yet opposing, impacts on local PII conformational biases. These results have implications for protein folding as well as for the function of signaling proteins that bind proline-rich targets in globular or intrinsically disordered proteins.

Conformational heterogeneity within the LID domain mediates substrate binding to Escherichia coli adenylate kinase: function follows fluctuations

Proteins exist as dynamic ensembles of molecules, implying that protein amino acid sequences evolved to code for both the ground state structure as well as the entire energy landscape of excited states. Accumulating theoretical and experimental evidence suggests that enzymes use such conformational fluctuations to facilitate allosteric processes important for substrate binding and possibly catalysis. This phenomenon can be clearly demonstrated in Escherichia coli adenylate kinase, where experimentally observed local unfolding of the LID subdomain, as opposed to a more commonly postulated rigid-body opening motion, is related to substrate binding. Because "entropy promoting" glycine mutations designed to increase specifically the local unfolding of the LID domain also affect substrate binding, changes in the excited energy landscape effectively tune the function of this enzyme without changing the ground state structure or the catalytic site. Thus, additional thermodynamic information, above and beyond the single folded structure of an enzyme-substrate complex, is likely required for a full and quantitative understanding of how enzymes work.

Interfacial water molecules in SH3 interactions: a revised paradigm for polyproline recognition

In spite of its biomedical relevance, polyproline recognition is still not fully understood. The disagreement between the current description of SH3 (Src homology 3) complexes and their thermodynamic behaviour calls for a revision of the SH3-binding paradigm. Recently, Abl-SH3 was demonstrated to recognize its ligands by a dual binding mechanism involving a robust network of water-mediated hydrogen bonds that complements the canonical hydrophobic interactions. The systematic analysis of the SH3 structural database in the present study reveals that this dual binding mode is universal to SH3 domains. Tightly bound buried-interfacial water molecules were found in all SH3 complexes studied mediating the interaction between the peptide ligand and the domain. Moreover, structural waters were also identified in a high percentage of the free SH3 domains. A detailed analysis of the pattern of water-mediated interactions enabled the identification of conserved hydration sites in the polyproline-recognition region and the establishment of relationships between hydration profiles and the sequence of both ligands and SH3 domains. Water-mediated interactions were also systematically observed in WW (protein–protein interaction domain containing two conserved tryptophan residues), UEV (ubiquitin-conjugating enzyme E2 variant) and EVH-1 [Ena/VASP (vasodilator-stimulated phosphoprotein) homology 1] structures. The results of the present study clearly indicate that the current description of proline-rich sequence recognition by protein–protein interaction modules is incomplete and insufficient for a correct understanding of these systems. A new binding paradigm is required that includes interfacial water molecules as relevant elements in polyproline recognition.

Strategies for the thermodynamic characterization of linked binding/local folding reactions within the native state application to the LID domain of adenylate kinase from Escherichia coli

Conformational fluctuations in proteins have emerged as an important aspect of biological function, having been linked to processes ranging from molecular recognition and catalysis to allostery and signal transduction. In spite of the realization of their importance, however, the connections between fluctuations and function have largely been empirical, even when they have been quantitative. Part of the problem in understanding the role of fluctuations in function is the fact that the mere existence of fluctuations complicates the interpretation of classic mutagenesis approaches. Namely, mutagenesis, which is typically targeted to an internal position (to elicit an effect), will change the fluctuations as well as the structure of the native state. Decoupling these effects is essential to an unambiguous understanding of the role of fluctuations in function. Here, we use a mutation strategy that targets surface-exposed sites in flexible parts of the molecule for mutation to glycine. Such mutations leave the ground-state structure unaffected. As a result, we can assess the nature of the fluctuations, develop a quantitative model relating fluctuations to function (in this case, molecular recognition), and unambiguously resolve the probabilities of the fluctuating states. We show that when this approach is applied to Escherichia coli adenylate kinase (AK), unique thermodynamic and structural insights are obtained, even when classic mutagenesis approaches targeted to the same region yield ambiguous results.

Characterizing the role of ensemble modulation in mutation-induced changes in binding affinity

Protein conformational fluctuations are key contributors to biological function, mediating important processes such as enzyme catalysis, molecular recognition, and allosteric signaling. To better understand the role of conformational fluctuations in substrate/ligand recognition, we analyzed, experimentally and computationally, the binding reaction between an SH3 domain and the recognition peptide of its partner protein. The fluctuations in this SH3 domain were enumerated by using an algorithm based on the hard sphere collision model, and the binding energetics resulting from these fluctuations were calculated using a structure-based energy function parametrized to solvent accessible surface areas. Surprisingly, this simple model reproduced the effects of mutations on the experimentally determined SH3 binding energetics, within the uncertainties of the measurements, indicating that conformational fluctuations in SH3, and in particular the RT loop region, are structurally diverse and are well-approximated by the randomly configured states. The mutated positions in SH3 were distant to the binding site and involved Ala and Gly substitutions of solvent exposed positions in the RT loop. To characterize these fluctuations, we applied principal coordinate analysis to the computed ensembles, uncovering the principal modes of conformational variation. It is shown that the observed differences in binding affinity between each mutant, and thus the apparent coupling between the mutated sites, can be described in terms of the changes in these principal modes. These results indicate that dynamic loops in proteins can populate a broad conformational ensemble and that a quantitative understanding of molecular recognition requires consideration of the entire distribution of states.

Denatured-state energy landscapes of a protein structural database reveal the energetic determinants of a framework model for folding

Position-specific denatured-state thermodynamics were determined for a database of human proteins by use of an ensemble-based model of protein structure. The results of modeling denatured protein in this manner reveal important sequence-dependent thermodynamic properties in the denatured ensembles as well as fundamental differences between the denatured and native ensembles in overall thermodynamic character. The generality and robustness of these results were validated by performing fold-recognition experiments, whereby sequences were matched with their respective folds based on amino acid propensities for the different energetic environments in the protein, as determined through cluster analysis. Correlation analysis between structure and energetic information revealed that sequence segments destined for beta-sheet in the final native fold are energetically more predisposed to a broader repertoire of states than are sequence segments destined for alpha-helix. These results suggest that within the subensemble of mostly unstructured states, the energy landscapes are dominated by states in which parts of helices adopt structure, whereas structure formation for sequences destined for beta-strand is far less probable. These results support a framework model of folding, which suggests that, in general, the denatured state has evolutionarily evolved to avoid low-energy conformations in sequences that ultimately adopt beta-strand. Instead, the denatured state evolved so that sequence segments that ultimately adopt alpha-helix and coil will have a high intrinsic structure formation capability, thus serving as potential nucleation sites.