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.

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. This article is protected by copyright. All rights reserved.

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.

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.

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.

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.

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.

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.

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.

Using the COREX/BEST server to model the native-state ensemble

Protein structures under normal conditions exist as ensembles of interconverting, transient microstates. A computer algorithm known as COREX/BEST (Biology using Ensemble-based Structural Thermodynamics) was developed to model microstate structures and describe the native ensembles of proteins in statistical thermodynamic terms. This algorithm has been tested extensively and validated through experimental comparisons examining a range of biophysical and functional phenomena, such as structural cooperativity, pH-dependent stability, and cold denaturation. Here, we describe a Web-based implementation of the COREX/BEST algorithm, called the COREX/BEST Server, and demonstrate how to use this online resource to characterize the structural and thermodynamic properties of the native protein ensemble.

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.

Thermal unfolding of the N-terminal region of p53 monitored by circular dichroism spectroscopy

It has been estimated that 30% of eukaryotic protein and 70% of transcription factors are intrinsically disordered (ID). The biochemical significance of proteins that lack stable tertiary structure, however, is not clearly understood, largely owing to an inability to assign well-defined structures to specific biological tasks. In an attempt to investigate the structural character of ID protein, we have measured the circular dichroism spectrum of the N-terminal region of p53 over a range of temperatures and solution conditions. p53 is a well-studied transcription factor that has a proline-rich N-terminal ID region containing two activation domains. High proline content is a property commonly associated with ID, and thus p53 may be a good model system for investigating the biochemical importance of ID. The spectra presented here suggest that the N-terminal region of p53 may adopt an ordered structure under physiological conditions and that this structure can be thermally unfolded in an apparent two-state manner. The midpoint temperature for this thermal unfolding of the N-terminal region of p53 was at the near-physiological temperature of 39°C, suggesting the possibility of a physiological role for the observed structural equilibrium.

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.

Thermodynamic principles for the engineering of pH-driven conformational switches and acid insensitive proteins

The general thermodynamic principles behind pH driven conformational transitions of biological macromolecules are well understood. What is less obvious is how they can be used to engineer pH switches in proteins. The acid unfolding of staphylococcal nuclease (SNase) was used to illustrate different factors that can affect pH-driven conformational transitions. Acid unfolding is a structural transition driven by preferential H(+) binding to the acid unfolded state (U) over the native (N) state of a protein. It is the result of carboxylic groups that titrate with more normal pK(a) values in the U state than in the N state. Acid unfolding profiles of proteins reflect a balance between electrostatic and non-electrostatic contributions to stability. Several strategies were used in attempts to turn SNase into an acid insensitive protein: (1) enhancing global stability of the protein with mutagenesis or with osmolytes, (2) use of high salt concentrations to screen Coulomb interactions, (3) stabilizing the N state through specific anion effects, (4) removing Asp or Glu residues that titrate with depressed pK(a) values in the N state, and (5) removing basic residues that might have strong repulsive interactions in the N state at low pH. The only effective way to engineer acid resistance in SNase is not through modulation of pK(a) values of Asp/Glu but by enhancing the global stability of the protein. Modulation of pH-driven conformational transitions by selective manipulation of the electrostatic component of the switch is an extremely difficult undertaking.

The role of protein conformational fluctuations in allostery, function, and evolution

It is now well-known that proteins exist at equilibrium as ensembles of conformational states rather than as unique static structures. Here we review from an ensemble perspective important biological effects of such spontaneous fluctuations on protein allostery, function, and evolution. However, rather than present a thorough literature review on each subject, we focus instead on connecting these phenomena through the ensemble-based experimental, theoretical, and computational investigations from our laboratory over the past decade. Special emphasis is given to insights that run counter to some of the prevailing ideas that have emerged over the past 40 years of structural biology research. For instance, when proteins are viewed as conformational ensembles rather than as single structures, the commonly held notion of an allosteric pathway as an obligate series of individual structural distortions loses its meaning. Instead, allostery can result from energetic linkage between distal sites as one Boltzmann distribution of states transitions to another. Additionally, the emerging principles from this ensemble view of proteins have proven surprisingly useful in describing the role of intrinsic disorder in inter-domain communication, functional adaptation mediated by mutational control of fluctuations, and evolutionary conservation of the energetics of protein stability.

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.

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

The left-handed polyproline II helical structure (P(II)) is observed to be a dominant conformation in the disordered states of protein and small polypeptide chains, even when no prolines are present in the sequence. Recently, in work by Ferreon and Hilser, the energetics associated with Ala and Gly substitutions at a surface exposed proline site were determined calorimetrically by measuring the binding energetics of Sos peptide variants to the C-terminal Src Homology 3 domain of SEM-5. The results were interpreted as a significant conformational bias toward the bound conformation (i.e., P(II)), even when the ligand is unbound. That study was not able to determine, however, whether the conformational bias of the peptides could be explained in terms other than that of a P(II) preference. Here, we test, using a computer algorithm based on the hard sphere collision (HSC) model, the notion of whether a bias in the unbound states of the peptide ligands is specific for the P(II) conformation, or if a bias to any other region of (phi, psi) space can also result in the same observed binding energetics. The results of these computer simulations indicate that, of the regions of (phi, psi) modeled for bias in the small peptides, only the bias to the P(II) conformation, and at rates of bias similar to the experimentally observed rates, quantitatively reproduced the experimental binding energetics.

Ligand effects on the protein ensemble: unifying the descriptions of ligand binding, local conformational fluctuations, and protein stability

Detailed description of the structural and physical basis of allostery, cooperativity, and other manifestations of long-range communication between binding sites in proteins remains elusive. Here we describe an ensemble-based structural-thermodynamic model capable of treating explicitly the coupling between ligand binding reactions, local fluctuations in structure, and global conformational transitions. The H(+) binding reactions of staphylococcal nuclease and the effects of pH on its stability were used to illustrate the properties of proteins that can be described quantitatively with this model. Each microstate in the native ensemble was modeled to have dual structural character; some regions were treated as folded and retained the same atomic geometry as in the crystallographic structure while other regions were treated thermodynamically as if they were unfolded. Two sets of pK(a) values were used to describe the affinity of each H(+) binding site. One set, calculated with a standard continuum electrostatics method, describes H(+) binding to sites in folded parts of the protein. A second set of pK(a) values, obtained from model compounds in water, was used to describe H(+) binding to sites in unfolded regions. An empirical free energy function, parameterized to reproduce folding thermodynamics measured by differential scanning calorimetry, was used to calculate the probability of each microstate. The effects of pH on the distribution of microstates were determined by the H(+) binding properties of each microstate. The validity of the calculations was established by comparison with a number of different experimental observables.

Functional residues serve a dominant role in mediating the cooperativity of the protein ensemble

Conformational fluctuations in proteins have emerged as a potentially important aspect of biological function, although the precise relationship and the implications have yet to be fully explored. Numerous studies have reported that the binding of ligand can influence fluctuations. However, the role of the binding site in mediating these fluctuations is not known. Of particular interest is whether in addition to serving as structural scaffolds for recognition and catalysis, active-site residues may also play a role in modulating the cooperative network. To address this question, we employ an experimentally validated ensemble-based description of proteins to elucidate the extent to which perturbations at different sites can influence the cooperative network in the protein. Applying this method to a database of test proteins, it is found statistically that binding sites are located in regions most able to affect the cooperative network, even for cooperative interactions between residues distant to the binding sites. This indicates that the conformational manifold under native conditions is determined by the network of cooperative interactions within the protein and suggests that proteins have evolved to use these conformational fluctuations in carrying out their functions. Furthermore, because the energetic coupling pattern calculated for each protein is robust and relatively insensitive to sequence, these studies further suggest that binding sites evolved in regions of the protein that are inherently poised to take advantage of the fluctuations in the native structure.

Revealing the nature of the native state ensemble through cold denaturation

Recent advances in NMR methodology have enabled the structural analysis of proteins at temperatures far below the freezing point of water, thus opening a window to the cold denaturation process. Although the phenomenon of cold denaturation has been known since the mid-1970s, the freezing point of water has prevented detailed and structurally resolved studies without application of additional significant perturbations of the protein ensemble. As a result, the cold-denatured state and the process of cold denaturation have gone largely unstudied. Here, the structural and thermodynamic basis of cold denaturation is explored with emphasis placed on the insights that are uniquely ascertained from low-temperature studies. It is shown that the noncooperative cold-induced unfolding of protein results in the population of partially folded states that cannot be accessed by other techniques. The structurally resolved view of the cold denaturation process therefore can provide direct access to the cooperative substructures within the protein molecule and provide an unprecedented structurally resolved picture of the states that comprise the native state ensemble.

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.

Low pH immobilizes and kills human leukocytes and prevents transmission of cell-associated HIV in a mouse model

Background

Both cell-associated and cell-free HIV virions are present in semen and cervical secretions of HIV-infected individuals. Thus, topical microbicides may need to inactivate both cell-associated and cell-free HIV to prevent sexual transmission of HIV/AIDS. To determine if the mild acidity of the healthy vagina and acid buffering microbicides would prevent transmission by HIV-infected leukocytes, we measured the effect of pH on leukocyte motility, viability and intracellular pH and tested the ability of an acidic buffering microbicide (BufferGel^®) to prevent the transmission of cell-associated HIV in a HuPBL-SCID mouse model.

Methods

Human lymphocyte, monocyte, and macrophage motilities were measured as a function of time and pH using various acidifying agents. Lymphocyte and macrophage motilities were measured using video microscopy. Monocyte motility was measured using video microscopy and chemotactic chambers. Peripheral blood mononuclear cell (PBMC) viability and intracellular pH were determined as a function of time and pH using fluorescent dyes. HuPBL-SCID mice were pretreated with BufferGel, saline, or a control gel and challenged with HIV-1-infected human PBMCs.

Results

Progressive motility was completely abolished in all cell types between pH 5.5 and 6.0. Concomitantly, at and below pH 5.5, the intracellular pH of PBMCs dropped precipitously to match the extracellular medium and did not recover. After acidification with hydrochloric acid to pH 4.5 for 60 min, although completely immotile, 58% of PBMCs excluded ethidium homodimer-1 (dead-cell dye). In contrast, when acidified to this pH with BufferGel, a microbicide designed to maintain vaginal acidity in the presence of semen, only 4% excluded dye at 10 min and none excluded dye after 30 min. BufferGel significantly reduced transmission of HIV-1 in HuPBL-SCID mice (1 of 12 infected) compared to saline (12 of 12 infected) and a control gel (5 of 7 infected).

Conclusion

These results suggest that physiologic or microbicide-induced acid immobilization and killing of infected white blood cells may be effective in preventing sexual transmission of cell-associated HIV.

Local conformational fluctuations can modulate the coupling between proton binding and global structural transitions in proteins

Local conformational fluctuations in proteins can affect the coupling between ligand binding and global structural transitions. This finding was established by monitoring quantitatively how the population distribution in the ensemble of microstates of staphylococcal nuclease was affected by proton binding. Analysis of acid unfolding and proton-binding data with an ensemble-based model suggests that local fluctuations: (i) can be effective modulators of ligand-binding affinities, (ii) are important determinants of the cooperativity of ligand-driven global structural transitions, and (iii) are well represented thermodynamically as local unfolding processes. These studies illustrate how an ensemble-based description of proteins can be used to describe quantitatively the interdependence of local conformational fluctuations, ligand-binding processes, and global structural transitions. This level of understanding of the relationship between conformation, energy, and dynamics is required for a detailed mechanistic understanding of allostery, cooperativity, and other complex functional and regulatory properties of macromolecules.

Thermodynamic mechanism and consequences of the polyproline II (PII) structural bias in the denatured states of proteins

A quantitative characterization of the structure and energy of the denatured states of proteins represents the cornerstone to a molecular-level understanding of both protein stability and fold specificity. Recent studies have revealed a significant bias in unstructured peptides toward the polyproline II (P(II)) conformation, even when no prolines are present in the sequence. This indicates that the P(II) conformation is a dominant component of the denatured states of proteins, although a quantitative description of the component enthalpy and entropy functions associated with this conformation (i.e., the thermodynamic mechanism) has thus far proven elusive. An experimental system has been designed that, when analyzed with high-precision isothermal titration calorimetry, provides direct access to the residue-specific thermodynamics of the P(II) structure formation in disordered proteins and peptides. Here, it is shown that the P(II) bias is driven by a favorable and significant enthalpy (Deltah) of -1.7 kcal mol(-1) residue(-1), which is partially offset by an unfavorable entropy (TDeltas) of -0.7 kcal mol(-1) residue(-1), relative to the ensemble of disordered conformations of the molecule. In addition to impacting dramatically the interpretation of thermal denaturation experiments, these experimental values form the framework of a quantitative energetic description of the denatured states of proteins.

The origin of pH-dependent changes in m-values for the denaturant-induced unfolding of proteins

Denaturant-induced unfolding is one of the most prevalent means of evaluating the structural stability of proteins and of determining the energetic consequences of mutations or changes in solution conditions. In spite of the widespread use of this approach, controversies and inconsistencies still persist with regard to the interpretation of the results of such studies. For example, most proteins show either a significant increase or a decrease (as much as 100 %) in the denaturant-dependence of the free energy of unfolding (i.e. the m-value) under increasingly acidic conditions. The pH dependence of the m-value is given different interpretations depending on whether the m-values increase or decrease with decreasing pH. In cases where m-values decrease, the decrease is attributed to the presence of an intermediate that becomes transiently stabilized during the unfolding transition at low pH. Cases where m-values increase as pH is lowered are usually interpreted in terms of an increase in the amount of surface area exposed by the denatured state at low pH. We have developed a general thermodynamic model that accounts for both types of behavior in terms of an intermediate that is populated throughout the unfolding transition. The model provides a unified framework for explaining both types of observed behavior, and the validity of the model was tested through the analysis of the pH dependence of m-values of staphylococcal nuclease. According to the model, the observed increase in m-values with decreasing pH is consistent with the existence of an intermediate that is populated during urea and guanidine unfolding. The intermediate becomes less populated during the unfolding transition at lower pH values giving rise to the apparent increase in m-values. These results argue that the prevailing interpretation need not apply to all proteins.

pH dependence of stability of staphylococcal nuclease: evidence of substantial electrostatic interactions in the denatured state

The pH dependence of stability of staphylococcal nuclease was studied with two independent equilibrium thermodynamic approaches. First, by measurement of stability in the pH range 9 to 3.5 by fluorescence-monitored denaturation with urea (Delta), GdnHCl (Delta), and heat (Delta). Second, by numerical integration of H(+) titration curves (Delta) measured potentiometrically under native (100 mM KCl) and unfolding (6.0 M GdnHCl) conditions. The pH dependence of stability described by Delta, Delta, and Delta was comparable but significantly different from the one described by Delta. The decrease in Delta between pH 9 and pH 4 was 4 kcal/mol greater than the decrease in Delta, Delta, and Delta in the same pH range. In 6 M GdnHCl, all the ionizable groups titrated with the pK(a) values of model compounds. Therefore, Delta represents the free energy difference between the native state (N) and an ensemble of unstructured, or expanded, and highly screened conformations. In contrast, the shallower pH dependence of stability described by Delta and by Delta between pH 9 and 5 was consistent with the titration of histidines with depressed, nativelike pK(a) values in the denatured state (D). These depressed pK(a) values likely reflect long-range electrostatic interactions with the other 29 basic groups and are a consequence of the compact character of the D state. The steep change in Delta and Delta at pH < 5 suggests that near pH 5 the structural and thermodynamic character of the D state shifts toward a state in which acidic residues titrate with normal pK(a) values, presumably because the electrostatic interactions with basic residues are lost, maybe as a consequence of an expansion.