Thermodynamic propensities of amino acids in the native state ensemble: implications for fold recognition.

An amino acid sequence, in the context of the solvent environment, contains all of the thermodynamic information necessary to encode a three-dimensional protein structure. To investigate the relationship between an amino acid sequence and its corresponding protein fold, a database of thermodynamic stability information was assembled that spanned 2951 residues from 44 nonhomologous proteins. This information was obtained using the COREX algorithm, which computes an ensemble-based description of the native state of a protein. It was observed that amino acid types partitioned unequally into high, medium, and low thermodynamic stability environments. Furthermore, these distributions were reproducible and were significantly different than those expected from random partitioning. To assess the structural importance of the distributions, simple fold-recognition experiments were performed based on a 3D-1D scoring matrix containing only COREX residue stability information. This procedure was able to recover amino acid sequences corresponding to correct target structures more effectively than scoring matrices derived from randomized data. High-scoring sequences were often aligned correctly with their corresponding target profiles, suggesting that calculated thermodynamic stability profiles have the potential to encode sequence information. As a control, identical fold-recognition experiments were performed on the same database of proteins using DSSP secondary structure information in the scoring matrix, instead of COREX residue stability information. The comparable performance of both approaches suggested that COREX residue stability information and secondary structure information could be of equivalent utility in more sophisticated fold-recognition techniques. The results of this work are a consequence of the idea that amino acid sequences fold not into single, rigidly stable structures but rather into thermodynamic ensembles best represented by a time-averaged structure.

Homology between O-linked GlcNAc transferases and proteins of the glycogen phosphorylase superfamily.

The O-linked GlcNAc transferases (OGTs) are a recently characterized group of largely eukaryotic enzymes that add a single beta-N-acetylglucosamine moiety to specific serine or threonine hydroxyls. In humans, this process may be part of a sugar regulation mechanism or cellular signaling pathway that is involved in many important diseases, such as diabetes, cancer, and neurodegeneration. However, no structural information about the human OGT exists, except for the identification of tetratricopeptide repeats (TPR) at the N terminus. The locations of substrate binding sites are unknown and the structural basis for this enzyme's function is not clear. Here, remote homology is reported between the OGTs and a large group of diverse sugar processing enzymes, including proteins with known structure such as glycogen phosphorylase, UDP-GlcNAc 2-epimerase, and the glycosyl transferase MurG. This relationship, in conjunction with amino acid similarity spanning the entire length of the sequence, implies that the fold of the human OGT consists of two Rossmann-like domains C-terminal to the TPR region. A conserved motif in the second Rossmann domain points to the UDP-GlcNAc donor binding site. This conclusion is supported by a combination of statistically significant PSI-BLAST hits, consensus secondary structure predictions, and a fold recognition hit to MurG. Additionally, iterative PSI-BLAST database searches reveal that proteins homologous to the OGTs form a large and diverse superfamily that is termed GPGTF (glycogen phosphorylase/glycosyl transferase). Up to one-third of the 51 functional families in the CAZY database, a glycosyl transferase classification scheme based on catalytic residue and sequence homology considerations, can be unified through this common predicted fold. GPGTF homologs constitute a substantial fraction of known proteins: 0.4% of all non-redundant sequences and about 1% of proteins in the Escherichia coli genome are found to belong to the GPGTF superfamily.

Ensemble modulation as an origin of denaturant-independent hydrogen exchange in proteins.

Native state hydrogen exchange (HX) has become a powerful tool for the analysis of conformational states that exist under native conditions. However, the interpretation of HX data in terms of conformational fluctuations is still controversial. In particular, it has been shown that many residues display exchange behavior that is independent of denaturant concentration. It has been postulated that this lack of denaturant dependence results from local fluctuations that do not expose appreciable amounts of buried surface area. Here, we use a general thermodynamic description of HX to explore the different possibilities for this behavior. We find that the denaturant dependence seen in HX experiments under native conditions is not a de facto indication of the amount of surface area exposure required for exchange. Instead, this behavior results from the relatively homogenous character of the conformational ensemble that exists under native conditions and the non-specific nature of denaturant effects. Furthermore, a comparison of the HX behavior from a stabilized mutant of Staphylococcal nuclease (SNase) with that predicted for the wild-type SNase from the COREX algorithm suggests that denaturant-independent exchange of many residues is consistent with significant (approximately 10 %) surface area exposure for this protein.

Correlation between changes in nuclear magnetic resonance order parameters and conformational entropy: molecular dynamics simulations of native and denatured staphylococcal nuclease.

Recent work has suggested that changes in NMR order parameters may quantitatively reflect changes in the conformational entropy of a protein ensemble. The extent of the mathematical relationship between local entropy changes as seen by NMR order parameters and the full protein entropy change is a complex issue. As a step towards a fuller understanding of this problem, molecular dynamics calculations of both native and denatured staphylococcal nuclease were performed. The N-H bond vector motion, in both explicit and implicit solvent, was analyzed to estimate local and global entropy changes. The calculated N-H bond vector order parameters from simulation agreed on average with experimental values for both native and denatured structures. However, the inverted-U profile of order parameters versus residue number observed experimentally for denatured nuclease was only partially reproduced by simulation of compact denatured structures. Comparisons made across the full set of simulations revealed a correlation between the N-H order parameter-based conformational entropy change and the total quasiharmonic-based conformational entropy change between the native and denatured structures. The calculations showed that about 25% of the total entropy change was reflected by changes in simulated S2 values. This result suggests that NMR-derived order parameters may be used to provide a reasonable estimate of the total conformational entropy change on protein folding.

A model of the changes in denatured state structure underlying m value effects in staphylococcal nuclease.

Hydrogen exchange kinetics were measured on the native states of wild type staphylococcal nuclease and four mutants with values of mGuHCl (defined as dDeltaG/d[guanidine hydrochloride]) ranging from 0.8 to 1.4 of the wild type value. Residues within the five-strand beta-barrel of wild type and E75A and D77A, two mutants with reduced values of m GuHCl, were significantly more protected from exchange than expected on the basis of global stability as measured by fluorescence. In contrast, mutants V23A and M26G with elevated values of mGuHCl approach a flat profile of more or less constant protection independent of position in the structure. Differences in exchange protection between the C-terminus and the beta-barrel region correlate with mGuHCl, suggesting that a residual barrel-like structure becomes more highly populated in the denatured states of m- mutants and less populated in m+ mutants. Variations in the population of such a molten globule-like structure would account for the large changes in solvent accessible surface area of the denatured state thought to underlie m value effects.

Perturbations of the denatured state ensemble: modeling their effects on protein stability and folding kinetics.

By considering the denatured state of a protein as an ensemble of conformations with varying numbers of sequence-specific interactions, the effects on stability, folding kinetics, and aggregation of perturbing these interactions can be predicted from changes in the molecular partition function. From general considerations, the following conclusions are drawn: (1) A perturbation that enhances a native interaction in denatured state conformations always increases the stability of the native state. (2) A perturbation that promotes a non-native interaction in the denatured state always decreases the stability of the native state. (3) A change in the denatured state ensemble can alter the kinetics of aggregation and folding. (4) The loss (or increase) in stability accompanying two mutations, each of which lowers (or raises) the free energy of the denatured state, will be less than the sum of the effects of the single mutations, except in cases where both mutations affect the same set of partially folded conformations. By modeling the denatured state as the ensemble of all non-native conformations of hydrophobic-polar (HP) chains configured on a square lattice, it can be shown that the stabilization obtained from enhancement of native interactions derives in large measure from the avoidance of non-native interactions in the D state. In addition, the kinetic effects of fixing single native contacts in the denatured state or imposing linear gradients in the HH contact probabilities are found, for some sequences, to significantly enhance the efficiency of folding by a simple hydrophobic zippering algorithm. Again, the dominant mechanism appears to be avoidance of non-native interactions. These results suggest stabilization of native interactions and imposition of gradients in the stability of local structure are two plausible mechanisms involving the denatured state that could play a role in the evolution of protein folding and stability.

Protein folding for realists: a timeless phenomenon.

Future research on protein folding must confront two serious dilemmas. (1) It may never be possible to observe at high resolution the very important structures that form in the first few milliseconds of the refolding reaction. (2) The energy functions used to predict structure from sequence will always be approximations of the true energy function. One strategy to resolve both dilemmas is to view protein folding from a different perspective, one that no longer emphasizes time and unique trajectories through conformation space. Instead, free energy replaces time as the reaction coordinate, and ensembles of equilibrium states of partially folded proteins are analyzed in place of trajectories of one protein chain through conformation space, either in vitro or in silico. Initial characterization of the folding of staphylococcal nuclease within this alternative conceptual framework has led to an equilibrium folding pathway with several surprising features. In addition to the finding of two bundles of four hydrophobic segments containing both native and non-native interactions, a gradient in relative stability of different substructures has been identified, with the most stable interactions located toward the amino terminus and the least stable toward the carboxy terminus. Hydrophobic bundles with up-down topology and stability gradients may be two examples of numerous tactics used by proteins to facilitate rapid folding and minimize aggregation. As NMR methods for structural analysis of partially folded proteins are refined, higher resolution descriptions of the structure and dynamics of the polypeptide chain outside the native state may provide many insights into the processes and energetics underlying the self-assembly of folded structure.

A scanning calorimetric study of unfolding equilibria in homodimeric chicken gizzard tropomyosins.

Using both circular dichroism (CD) and differential scanning calorimetry (DSC), several laboratories find that the thermal unfolding transitions of alpha alpha and beta beta homodimeric coiled coils of rabbit tropomyosin are multistate and display an overall unfolding enthalpy of near 300 kcal (mol dimer)(-1). In contrast, an extant CD study of beta beta and gamma gamma species of chicken gizzard tropomyosin concludes that their unfolding transitions are simple two-state transitions, with much smaller overall enthalpies (98 kcal mol(-1) for beta beta and 162 kcal mol(-1) for gamma gamma). However, these smaller enthalpies have been questioned, because they imply a concentration dependence of the melting temperatures that is far larger than observed by CD. We report here DSC studies of the unfolding of both beta beta and gamma gamma chicken gizzard homodimers. The results show that these transitions are very similar to those in rabbit tropomyosins in that 1) the overall unfolding enthalpy is near 300 kcal mol(-1); 2) the overall delta C(rho) values are significantly positive; 3) the various transitions are multistate, requiring at least two and as many as four domains to fit the DSC data. DSC studies are also reported on these homodimeric species of chicken gizzard tropomyosin with a single interchain disulfide cross-link. These results are also generally similar to those for the correspondingly cross-linked rabbit tropomyosins.

Thermal unfolding equilibria in homodimeric chicken gizzard tropomyosin coiled coils.

CD studies are presented on thermal unfolding of coiled-coil homodimers of two genetic variant chains of chicken gizzard tropomyosin (CG-Tm). The experiments include the effects of cross-linking both isoforms and the dependence on protein concentration of unfolding in both reduced isoforms, variables not examined in extant work. The general shapes of the unfolding curves for singly cross-linked species depend on whether the cross-link is at C190 (its site on one isoform) or at C36 (its site on the other). These curves are compared with extant ones for various cross-linked species of rabbit tropomyosin. The comparison supports the view that the unfolding behavior of cross-linked species results from a complex interaction of strain at the cross-link, local variations in structural stability, and loop entropy. The observed concentration dependence of the transition temperature for the uncross-linked (reduced) species of CG-Tm is very small (2.9 degrees C) for one variant homodimer and unobservably small (< 2 degrees C) for the other in the 100-fold concentration range (approximately 0.01-1.0 mg/mL) accessible here. These experimental values of delta Tm are much smaller than are predicted from extant values of the van't Hoff transition enthalpies, calling the latter into question.