Search

Abstract 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.

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., (© 2023 The Authors. Protein Science published by Wiley Periodicals LLC on behalf of The Protein Society.)

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., (© 2023 The Authors. Protein Science published by Wiley Periodicals LLC on behalf of The Protein Society.)

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., Competing Interests: Conflict of interest The authors declare that they have no conflicts of interest with the contents of this article., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

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 π/sp 2 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., Competing Interests: Conflict of interest The authors declare that they have no conflicts of interest with the contents of this article., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

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.

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., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

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 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., (Copyright © 2018 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

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., (© 2016 Wiley Periodicals, Inc.)

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., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

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.

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)., (Copyright © 2015 Elsevier Ltd. All rights reserved.)

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., (Copyright © 2015 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

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., (© 2014 Wiley Periodicals, Inc.)

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., (Copyright © 2013 Wiley Periodicals, Inc.)

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.

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., (Copyright © 2012 The Protein Society.)

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., (Copyright © 2011 Wiley Periodicals, Inc.)

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., (Copyright © 2011. Published by Elsevier B.V.)

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., (Copyright © 2011 Elsevier B.V. All rights reserved.)

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.

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.

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.

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.

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.

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.

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., (2006 Wiley-Liss, Inc.)

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., ((c) 2005 Wiley-Liss, Inc.)

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.

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.

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., (Copyright 2001 Academic Press.)

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.