Intrinsically disordered regions that drive phase separation form a robustly distinct protein class.

Protein phase separation is thought to be a primary driving force for the formation of membrane-less organelles, which control a wide range of biological functions from stress response to ribosome biogenesis. Among phase-separating (PS) proteins, many have intrinsically disordered regions (IDRs) that are needed for phase separation to occur. Accurate identification of IDRs that drive phase separation is important for testing the underlying mechanisms of phase separation, identifying biological processes that rely on phase separation, and designing sequences that modulate phase separation. To identify IDRs that drive phase separation, we first curated datasets of folded, ID, and PS ID sequences. We then used these sequence sets to examine how broadly existing amino acid property scales can be used to distinguish between the three classes of protein regions. We found that there are robust property differences between the classes and, consequently, that numerous combinations of amino acid property scales can be used to make robust predictions of protein phase separation. This result indicates that multiple, redundant mechanisms contribute to the formation of phase-separated droplets from IDRs. The top-performing scales were used to further optimize our previously developed predictor of PS IDRs, ParSe. We then modified ParSe to account for interactions between amino acids and obtained reasonable predictive power for mutations that have been designed to test the role of amino acid interactions in driving protein phase separation. Collectively, our findings provide further insight into the classification of IDRs and the elements involved in protein phase separation.

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.

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.

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 (Rh ) when measured under normal conditions. To investigate sequence-modulation of IDP Rh , 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 Rh 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 (PPII ), however, modulated Rh with intricate dependence on the simulated propensities. PPII propensities at levels expected in IDPs correlated with heightened Rh sensitivity to even weak α helix propensities, indicating bias for common (φ, ψ) are important determinants of hydrodynamic size. Moreover, data show that IDP Rh 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.

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 (PPII) structure. While intrinsic PPII 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 (Rh) can be predicted from experimental PPII propensities with good agreement, even when charge-based considerations are omitted. The simulations demonstrate that Rh and chain propensity for PPII structure are linked via a simple power-law scaling relationship, which was tested using the experimental Rh of 22 IDPs covering a wide range of peptide lengths, net charge, and sequence composition. Charge effects on Rh were found to be generally weak when compared to PPII effects on Rh. Results from this study indicate that the hydrodynamic dimensions of IDPs are evidence of considerable sequence-dependent backbone propensities for PPII structure that qualitatively, if not quantitatively, match conformational propensities measured in peptides.

Structural origins of misfolding propensity in the platelet adhesive von Willebrand factor A1 domain.

The von Willebrand factor (VWF) A1 and A3 domains are structurally isomorphic yet exhibit distinct mechanisms of unfolding. The A1 domain, responsible for platelet adhesion to VWF in hemostasis, unfolds through a molten globule intermediate in an apparent three-state mechanism, while A3 unfolds by a classical two-state mechanism. Inspection of the sequences or structures alone does not elucidate the source of this thermodynamic conundrum; however, the three-state character of the A1 domain suggests that it has more than one cooperative substructure yielding two separate unfolding transitions not present in A3. We investigate the extent to which structural elements contributing to intermediate conformations can be identified using a residue-specific implementation of the structure-energy-equivalence-of-domains algorithm (SEED), which parses proteins of known structure into their constituent thermodynamically cooperative components using protein-group-specific, transfer free energies. The structural elements computed to contribute to the non-two-state character coincide with regions where Von Willebrand disease mutations induce misfolded molten globule conformations of the A1 domain. This suggests a mechanism for the regulation of rheological platelet adhesion to A1 based on cooperative flexibility of the α2 and α3 helices flanking the platelet GPIbα receptor binding interface.

Fungal metabolite ophiobolin A as a promising anti-glioma agent: In vivo evaluation, structure-activity relationship and unique pyrrolylation of primary amines.

Glioblastoma, the most common form of malignant primary brain tumor, is characterized by resistance to apoptosis, which is largely responsible for the low effectiveness of the classical chemotherapeutic approaches based on apoptosis induction in cancer cells. Previously, a fungal secondary metabolite ophiobolin A was found to have significant activity against apoptosis-resistant glioblastoma cells through the induction of a non-apoptotic cell death, thus, offering an innovative strategy to combat this type of cancer. The current work describes the results of a preliminary evaluation of ophiobolin A in an in vivo glioblastoma model and its chemical derivatization to establish first synthetically generated structure-activity relationship. The synthetic work has also led to the discovery of a unique reaction of ophiobolin A with primary amines suggesting the possibility of pyrrolylation of lysine residues on its intracellular target protein(s).

Alanine and proline content modulate global sensitivity to discrete perturbations in disordered proteins.

Molecular transduction of biological signals is understood primarily in terms of the cooperative structural transitions of protein macromolecules, providing a mechanism through which discrete local structure perturbations affect global macromolecular properties. The recognition that proteins lacking tertiary stability, commonly referred to as intrinsically disordered proteins (IDPs), mediate key signaling pathways suggests that protein structures without cooperative intramolecular interactions may also have the ability to couple local and global structure changes. Presented here are results from experiments that measured and tested the ability of disordered proteins to couple local changes in structure to global changes in structure. Using the intrinsically disordered N-terminal region of the p53 protein as an experimental model, a set of proline (PRO) and alanine (ALA) to glycine (GLY) substitution variants were designed to modulate backbone conformational propensities without introducing non-native intramolecular interactions. The hydrodynamic radius (R(h)) was used to monitor changes in global structure. Circular dichroism spectroscopy showed that the GLY substitutions decreased polyproline II (PP(II)) propensities relative to the wild type, as expected, and fluorescence methods indicated that substitution-induced changes in R(h) were not associated with folding. The experiments showed that changes in local PP(II) structure cause changes in R(h) that are variable and that depend on the intrinsic chain propensities of PRO and ALA residues, demonstrating a mechanism for coupling local and global structure changes. Molecular simulations that model our results were used to extend the analysis to other proteins and illustrate the generality of the observed PRO and alanine effects on the structures of IDPs.

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.

Convergent behavior of extended stalk regions from staphylococcal surface proteins with widely divergent sequence patterns.

Staphylococcus epidermidis and S. aureus are highly problematic bacteria in hospital settings. This stems, at least in part, from strong abilities to form biofilms on abiotic or biotic surfaces. Biofilms are well-organized multicellular aggregates of bacteria, which, when formed on indwelling medical devices, lead to infections that are difficult to treat. Cell wall-anchored (CWA) proteins are known to be important players in biofilm formation and infection. Many of these proteins have putative stalk-like regions or regions of low complexity near the cell wall-anchoring motif. Recent work demonstrated the strong propensity of the stalk region of the S. epidermidis accumulation-associated protein (Aap) to remain highly extended under solution conditions that typically induce compaction or other significant conformational changes. This behavior is consistent with the expected function of a stalk-like region that is covalently attached to the cell wall peptidoglycan and projects the adhesive domains of Aap away from the cell surface. In this study, we evaluate whether the ability to resist compaction is a common theme among stalk regions from various staphylococcal CWA proteins. Circular dichroism spectroscopy was used to examine secondary structure changes as a function of temperature and cosolvents along with sedimentation velocity analytical ultracentrifugation and SAXS to characterize structural characteristics in solution. All stalk regions tested are intrinsically disordered, lacking secondary structure beyond random coil and polyproline type II helix, and they all sample highly extended conformations. Remarkably, the Ser-Asp dipeptide repeat region of SdrC exhibited nearly identical behavior in solution when compared to the Aap Pro/Gly-rich region, despite highly divergent sequence patterns, indicating conservation of function by various distinct staphylococcal CWA protein stalk regions.

Convergent behavior of extended stalk regions from staphylococcal surface proteins with widely divergent sequence patterns.

Staphylococcus epidermidis and Staphylococcus aureus are highly problematic bacteria in hospital settings. A major challenge is their ability to form biofilms on abiotic or biotic surfaces. Biofilms are well-organized, multicellular bacterial aggregates that resist antibiotic treatment and often lead to recurrent infections. Bacterial cell wall-anchored (CWA) proteins are important players in biofilm formation and infection. Many have putative stalk-like regions or regions of low complexity near the cell wall-anchoring motif. Recent work demonstrated the strong propensity of the stalk region of S. epidermidis accumulation-associated protein (Aap) to remain highly extended under solution conditions that typically induce compaction. This behavior is consistent with the expected function of a stalk-like region that is covalently attached to the cell wall peptidoglycan and projects the adhesive domains of Aap away from the cell surface. In this study, we evaluate whether the ability to resist compaction is a common theme among stalk regions from various staphylococcal CWA proteins. Circular dichroism spectroscopy was used to examine secondary structure changes as a function of temperature and cosolvents along with sedimentation velocity analytical ultracentrifugation, size-exclusion chromatography, and SAXS to characterize structural characteristics in solution. All stalk regions tested are intrinsically disordered, lacking secondary structure beyond random coil and polyproline type II helix, and they all sample highly extended conformations. Remarkably, the Ser-Asp dipeptide repeat region of SdrC exhibited nearly identical behavior in solution when compared to the Aap Pro/Gly-rich region, despite highly divergent sequence patterns, indicating conservation of function by various distinct staphylococcal CWA protein stalk regions.

ParSe 2.0: A web tool to identify drivers of protein phase separation at the proteome level.

We have developed an algorithm, ParSe, which accurately identifies from the primary sequence those protein regions likely to exhibit physiological phase separation behavior. Originally, ParSe was designed to test the hypothesis that, for flexible proteins, phase separation potential is correlated to hydrodynamic size. While our results were consistent with that idea, we also found that many different descriptors could successfully differentiate between three classes of protein regions: folded, intrinsically disordered, and phase-separating intrinsically disordered. Consequently, numerous combinations of amino acid property scales can be used to make robust predictions of protein phase separation. Built from that finding, ParSe 2.0 uses an optimal set of property scales to predict domain-level organization and compute a sequence-based prediction of phase separation potential. The algorithm is fast enough to scan the whole of the human proteome in minutes on a single computer and is equally or more accurate than other published predictors in identifying proteins and regions within proteins that drive phase separation. Here, we describe a web application for ParSe 2.0 that may be accessed through a browser by visiting https://stevewhitten.github.io/Parse_v2_FASTA to quickly identify phase-separating proteins within large sequence sets, or by visiting https://stevewhitten.github.io/Parse_v2_web to evaluate individual protein sequences.

Mutational analysis of m-values as a strategy to identify cold-resistant substructures of the protein ensemble.

Characterizing the native ensemble of protein is an important yet difficult objective of structural biology. The structural dynamics of protein macromolecules play key roles in biological function, but the short lifetimes and low population of near-native states of the protein ensemble limit their ability to be studied directly. In part to address such issues, it was shown recently that the cooperative substructures that populate a protein ensemble could be ascertained by NMR methods performed at very cold temperatures. What is presented here is an argument that these same substructures can also be determined by denaturant-induced unfolding studies performed on protein at room temperature. Data supporting this argument are given for Staphylococcal nuclease, chymotrypsin inhibitor 2, and ubiquitin. The observation of an agreement between the thermodynamics of the protein ensemble simulated under very cold temperatures to the apparent sensitivity of the ensemble to chemical denaturants at room temperature also suggests that the overall structural-thermodynamic character of an ensemble is surprisingly robust and preserved even in the presence of strong denaturing conditions.

Evidence for the Misfolding of the A1 Domain within Multimeric von Willebrand Factor in Type 2 von Willebrand Disease.

Von Willebrand factor (VWF), an exceptionally large multimeric plasma glycoprotein, functions to initiate coagulation by agglutinating platelets in the blood stream to sites of vascular injury. This primary hemostatic function is perturbed in type 2 dysfunctional subtypes of von Willebrand disease (VWD) by mutations that alter the structure and function of the platelet GPIbα adhesive VWF A1 domains. The resulting amino acid substitutions cause local disorder and misfold the native structure of the isolated platelet GPIbα-adhesive A1 domain of VWF in both gain-of-function (type 2B) and loss-of-function (type 2M) phenotypes. These structural effects have not been explicitly observed in A1 domains of VWF multimers native to blood plasma. New mass spectrometry strategies are applied to resolve the structural effects of 2B and 2M mutations in VWF to verify the presence of A1 domain structural disorder in multimeric VWF harboring type 2 VWD mutations. Limited trypsinolysis mass spectrometry (LTMS) and hydrogen-deuterium exchange mass spectrometry (HXMS) are applied to wild-type and VWD variants of the single A1, A2, and A3 domains, an A1A2A3 tridomain fragment of VWF, plasmin-cleaved dimers of VWF, multimeric recombinant VWF, and normal VWF plasma concentrates. Comparatively, these methods show that mutations known to misfold the isolated A1 domain increase the rate of trypsinolysis and the extent of hydrogen-deuterium exchange in local secondary structures of A1 within multimeric VWF. VWD mutation effects are localized to the A1 domain without appreciably affecting the structure and dynamics of other VWF domains. The intrinsic dynamics of A1 observed in recombinant fragments of VWF are conserved in plasma-derived VWF. These studies reveal that structural disorder does occur in VWD variants of the A1 domain within multimeric VWF and provides strong support for VWF misfolding as a result of some, but not all, type 2 VWD variants.

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.

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.

The Proline/Glycine-Rich Region of the Biofilm Adhesion Protein Aap Forms an Extended Stalk that Resists Compaction.

Staphylococcus epidermidis is one of the primary bacterial species responsible for healthcare-associated infections. The most significant virulence factor for S. epidermidis is its ability to form a biofilm, which renders the bacteria highly resistant to host immune responses and antibiotic action. Intercellular adhesion within the biofilm is mediated by the accumulation-associated protein (Aap), a cell wall-anchored protein that self-assembles in a zinc-dependent manner. The C-terminal portion of Aap contains a 135-aa-long, proline/glycine-rich region (PGR) that has not yet been characterized. The region contains a set of 18 nearly identical AEPGKP repeats. Analysis of the PGR using biophysical techniques demonstrated the region is a highly extended, intrinsically disordered polypeptide with unusually high polyproline type II helix propensity. In contrast to many intrinsically disordered polypeptides, there was a minimal temperature dependence of the global conformational state of PGR in solution as measured by analytical ultracentrifugation and dynamic light scattering. Furthermore, PGR was resistant to conformational collapse or α-helix formation upon the addition of the osmolyte trimethylamine N-oxide or the cosolvent 2,2,2-trifluoroethanol. Collectively, these results suggest PGR functions as a resilient, extended stalk that projects the rest of Aap outward from the bacterial cell wall, promoting intercellular adhesion between cells in the biofilm. This work sheds light on regions of low complexity often found near the attachment point of bacterial cell wall-anchored proteins.

Ensemble-based signatures of energy propagation in proteins: a new view of an old phenomenon.

The ability of a protein to transmit the energetic effects of binding from one site to another constitutes the underlying basis for allosterism and signal transduction. Despite clear experimental evidence indicating the ability of proteins to transmit the effects of binding, the means by which this propagation is facilitated is not well understood. Using our previously developed ensemble-based description of the equilibrium, we investigated the physical basis of energy propagation and identified several fundamental and general aspects of energetic coupling between residues in a protein. First, partitioning of a conformational ensemble into four distinct sub-ensembles allows for explanation of the range of experimentally observed coupling behaviors (i.e., positive, neutral, and negative coupling between various regions of the protein structure). Second, the relative thermodynamic properties of these four sub-ensembles define the energetic coupling between residues as either positive, neutral, or negative. Third, analysis of the structural and thermodynamic features of the states within each sub-ensemble reveals significant variability. This third result suggests that a quantitative description of energy propagation in proteins requires an understanding of the structural and energetic features of more than just one or a few low-energy states, but also of many high-energy states. Such findings illuminate the difficulty in interpreting energy propagation in proteins in terms of a structural pathway that physically links coupled sites.

Molecular mechanisms of pH-driven conformational transitions of proteins: insights from continuum electrostatics calculations of acid unfolding.

The acid unfolding of staphylococcal nuclease (SNase) is very cooperative (Whitten and García-Moreno, Biochemistry 2000;39:14292-14304). As many as seven hydrogen ions (H+) are bound preferentially by the acid-unfolded state relative to the native (N) state in the pH range 3.2-3.9. To investigate the mechanism of acid unfolding, structure-based pKa calculations were performed with a variety of continuum electrostatic methods. The calculations reproduced successfully the H+ binding properties of the N state between pH 5 and 9, but they systematically overestimated the number of H+ bound upon acid unfolding. The calculated pKa values of all carboxylic residues in the N state were more depressed than they should be. The discrepancy between the observed and the calculated H+ uptake upon acid unfolding was not improved by using high protein dielectric constants, structures relaxed with molecular dynamics, or other empirical modifications implemented previously by others to maximize agreement between measured and calculated pKa values. This suggests an important role for conformational fluctuations of the backbone as important determinants of pKa values of carboxylic groups. Because no global or subglobal conformational changes have been observed previously for SNase under acidic conditions above the acid-unfolding region, these fluctuations must be local. The acid unfolding of SNase does not seem to involve the disruption of the N state by accruement of intramolecular repulsive interactions, nor the protonation of key ion paired carboxylic residues. It is more consistent with modest contributions from many H+ binding groups, with an important role for local conformational fluctuations in the coupling between H+ binding and the global structural transition.