Hydrogen exchange in native and alcohol forms of ubiquitin

Characterization of the interface structure of enzyme-inhibitor complex by using hydrogen-deuterium exchange and electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry

We investigated the interaction between a thiol protease inhibitor, cystatin, and its target enzyme, papain, by hydrogen-deuterium (H/D) exchange in conjunction with successive analysis by collision-induced dissociation (CID) in an rf-only hexapole ion guide with electrospray ionization-Fourier transform ion cyclotron resonance mass spectrometry (ESI-FTICR MS). The deuterium incorporation into backbone amide hydrogens of cystatin was analyzed at different time points in the presence or absence of papain, examining the mass of each fragment produced by hexapole-CID. In the absence of papain, amide hydrogens in short amino-terminal fragments, such as b10(2+) and b12(2+), were highly deuterated within 1 min. Although fewer fragments were observed for the cystatin-papain complex in the hexapole-CID spectra, significant reductions in initial deuterium content were recognized throughout the sequence of cystatin. This suggests that complex formation restricted the flexibility of the whole cystatin molecule. Detailed analyses revealed that a marked reduction in deuterium content in the region of residues 1-10 persisted for hours, suggesting that the flexible N-terminal region was tightly fixed in the binding pocket with hydrogen bonds. Our results are consistent with those of previous studies on the structure and inhibition mechanism of cystatin. We demonstrated here that enzyme-inhibitor interactions can be characterized by H/D exchange in combination with CID in a hexapole ion guide using ESI-FTICR MS rapidly and using only a small amount of sample.

Identification of side-chain clusters in protein structures by a graph spectral method

This paper presents a novel method to detect side-chain clusters in protein three-dimensional structures using a graph spectral approach. Protein side-chain interactions are represented by a labeled graph in which the nodes of the graph represent the Cbeta atoms and the edges represent the distance between the Cbeta atoms. The distance information and the non-bonded connectivity of the residues are represented in the form of a matrix called the Laplacian matrix. The constructed matrix is diagonalized and clustering information is obtained from the vector components associated with the second lowest eigenvalue and cluster centers are obtained from the vector components associated with the top eigenvalues. The method uses global information for clustering and a single numeric computation is required to detect clusters of interest. The approach has been adopted here to detect a variety of side-chain clusters and identify the residue which makes the largest number of interactions among the residues forming the cluster (cluster centers). Detecting such clusters and cluster centers are important from a protein structure and folding point of view. The crucial residues which are important in the folding pathway as determined by PhiF values (which is a measure of the effect of a mutation on the stability of the transition state of folding) as obtained from protein engineering methods, can be identified from the vector components corresponding to the top eigenvalues. Expanded clusters are detected near the active and binding site of the protein, supporting the nucleation condensation hypothesis for folding. The method is also shown to detect domains in protein structures and conserved side-chain clusters in topologically similar proteins.

Solution structure and dynamics of a designed hydrophobic core variant of ubiquitin

BACKGROUND

The recent merger of computation and protein design has resulted in a burst of success in the generation of novel proteins with native-like properties. A critical component of this coupling between theory and experiment is a detailed analysis of the structures and stabilities of designed proteins to assess and improve the accuracy of design algorithms.

RESULTS

Here we report the solution structure of a hydrophobic core variant of ubiquitin, referred to as 1D7, which was designed with the core-repacking algorithm ROC. As a measure of conformational specificity, we also present amide exchange protection factors and backbone and sidechain dynamics. The results indicate that 1D7 is similar to wild-type (WT) ubiquitin in backbone structure and degree of conformational specificity. We also observe a good correlation between experimentally determined sidechain structures and those predicted by ROC. However, evaluation of the core sidechain conformations indicates that, in general, 1D7 has more sidechains in less statistically favorable conformations than WT.

CONCLUSIONS

Our results provide an explanation for the lower stability of 1D7 compared to WT, and suggest modifications to design algorithms that may improve the accuracy with which structure and stability are predicted. The results also demonstrate that core packing can affect conformational flexibility in subtle ways that are likely to be important for the design of function and protein-ligand interactions.

The hydrogen exchange core and protein folding

A database of hydrogen-deuterium exchange results has been compiled for proteins for which there are published rates of out-exchange in the native state, protection against exchange during folding, and out-exchange in partially folded forms. The question of whether the slow exchange core is the folding core (Woodward C, 1993, Trends Biochem Sci 18:359-360) is reexamined in a detailed comparison of the specific amide protons (NHs) and the elements of secondary structure on which they are located. For each pulsed exchange or competition experiment, probe NHs are shown explicitly; the large number and broad distribution of probe NHs support the validity of comparing out-exchange with pulsed-exchange/competition experiments. There is a strong tendency for the same elements of secondary structure to carry NHs most protected in the native state, NHs first protected during folding, and NHs most protected in partially folded species. There is not a one-to-one correspondence of individual NHs. Proteins for which there are published data for native state out-exchange and theta values are also reviewed. The elements of secondary structure containing the slowest exchanging NHs in native proteins tend to contain side chains with high theta values or be connected to a turn/loop with high theta values. A definition for a protein core is proposed, and the implications for protein folding are discussed. Apparently, during folding and in the native state, nonlocal interactions between core sequences are favored more than other possible nonlocal interactions. Other studies of partially folded bovine pancreatic trypsin inhibitor (Barbar E, Barany G, Woodward C, 1995, Biochemistry 34:11423-11434; Barber E, Hare M, Daragan V, Barany G, Woodward C, 1998, Biochemistry 37:7822-7833), suggest that developing cores have site-specific energy barriers between microstates, one disordered, and the other(s) more ordered.

Predicting the protein folding nucleus from sequences [correction of a sequence]

Understanding the mechanism of protein folding would allow prediction of the three-dimensional structure from sequence data alone. It has been shown that small proteins fold in a small number of kinetic steps and that significantly populated intermediate states exist for some of them. Studies of these intermediates have demonstrated the existence of specific interactions established during the initial stages of folding. Comparison of the amino acids participating in these specific and essential interactions and constituting the folding nucleus with conserved hydrophobic positions of a given fold shows a striking correspondence. This finding opens the perspective of predicting the folding nucleus knowing only a set of divergent sequences of a protein family.

Unfolding of proteins monitored by electrospray ionization mass spectrometry: a comparison of positive and negative ion modes

Electrospray ionization (ESI) mass spectrometry (MS) in both the positive and negative ion mode has been used to study protein unfolding transitions of lysozyme, cytochrome c (cyt c), and ubiquitin in solution. As expected, ESI of unfolded lysozyme leads to the formation of substantially higher charge states than the tightly folded protein in both modes of operation. Surprisingly, the acid-induced unfolding of cyt c as well as the acid and the base-induced unfolding of ubiquitin show different behavior: In these three cases protein unfolding only leads to marginal changes in the negative ion charge state distributions, whereas in the positive ion mode pronounced shifts to higher charge states are observed. This shows that ESI MS in the negative ion mode as a method for probing conformational changes of proteins in solution should be treated with caution. The data presented in this work provide further evidence that the conformation of a protein in solution not its charge state is the predominant factor for determining the ESI charge state distribution in the positive ion mode. Furthermore, these data support the hypothesis of a recent study (Konermann and Douglas, Biochemistry 1997, 36, 12,296-12,302) which suggested that ESI in the positive ion mode is not sensitive to changes in the secondary structure of proteins but only to changes in the tertiary structure.

High-resolution, high-pressure NMR studies of proteins

Advanced high-resolution NMR spectroscopy, including two-dimensional NMR techniques, combined with high pressure capability, represents a powerful new tool in the study of proteins. This contribution is organized in the following way. First, the specialized instrumentation needed for high-pressure NMR experiments is discussed, with specific emphasis on the design features and performance characteristics of a high-sensitivity, high-resolution, variable-temperature NMR probe operating at 500 MHz and at pressures of up to 500 MPa. An overview of several recent studies using 1D and 2D high-resolution, high-pressure NMR spectroscopy to investigate the pressure-induced reversible unfolding and pressure-assisted cold denaturation of lysozyme, ribonuclease A, and ubiquitin is presented. Specifically, the relationship between the residual secondary structure of pressure-assisted, cold-denatured states and the structure of early folding intermediates is discussed.

Hydrogen exchange studies of protein structure

Hydrogen exchange techniques, with their residue-level specificity, exquisite sensitivity, and adaptability to many solution conditions, are becoming essential to the study of protein stability, folding and dynamics. Recent studies have elucidated the structures of intermediates formed transiently during protein folding and rare partially folded ensembles present at equilibrium. Analysis of hydrogen exchange mechanisms has revealed protein stability and kinetics at the level of individual residues.

A designed four helix bundle protein with native-like structure

A 108 amino acid protein was designed and constructed from a reduced alphabet of seven amino acids. The 2.9 A resolution X-ray crystal structure confirms that the protein is a four helix bundle, as it was designed to be. Hydrogen/deuterium exchange experiments reveal buried amide protons with protection factors in excess of 1 x 10(6) in the range characteristic of well protected protons in functional folded proteins (10(3)-10(8)) rather than protons in rapid exchange (0-10(2)). The protein is monomeric at 1 mM, the concentration at which the exchange experiments were undertaken, indicating that the exchange factors are due to a unique stable tertiary structure fold, and not due to any higher order quaternary structure. Thermodynamic analysis provides an estimate of the free energy of folding of -9.3 kcal mole-1 at 25 degrees C, consistent with the free energy of folding derived from the protection factors of the most protected protons, indicating that global unfolding is required for exchange of the most protected protons.

Structure of the pressure-assisted cold denatured state of ubiquitin

The pressure-assisted cold denatured state of ubiquitin in aqueous solution was investigated by high resolution NMR. Hydrogen exchange kinetics were measured for backbone amide protons in the cold denatured protein to determine its structure. In contrast to cold denatured ribonuclease A and lysozyme, cold denatured ubiquitin shows little persistent secondary structure. The behavior of ubiquitin supports the idea of a relationship between the residual structure of pressure-assisted cold-denatured states and the structure of early folding intermediates provided they exist.

Destabilisation of native tertiary structural interactions is linked to helix-induction by 2,2,2-trifluoroethanol in proteins

The effect of 2,2,2-trifluoroethanol (TFE) on the structure of an all beta-sheet protein, cardiotoxin analogue 111 (CTX III) from the Taiwan cobra (Naja naja atra) is studied. It is found that high concentrations (> 80% v/v) of TFE induced a beta-sheet to alpha-helix structural transition. It is found that in denatured and reduced CTX III (rCTX III) helical conformation is induced even upon addition of low concentrations (> 10% v/v) of TFE. Using three other proteins, namely, ribonuclease A (RNase A), lysozyme and alpha-lactalbumin, it is been observed that helix-induction by TFE is intricately linked to drastic destabilization of native tertiary structural interactions in the proteins.

Insights into protein folding from NMR

NMR has emerged as an important tool for studies of protein folding because of the unique structural insights it can provide into many aspects of the folding process. Applications include measurements of kinetic folding events and structural characterization of folding intermediates, partly folded states, and unfolded states. Kinetic information on a time scale of milliseconds or longer can be obtained by real-time NMR experiments and by quench-flow hydrogen-exchange pulse labeling. Although NMR cannot provide direct information on the very rapid processes occurring during the earliest stages of protein folding, studies of isolated peptide fragments provide insights into likely protein folding initiation events. Multidimensional NMR techniques are providing new information on the structure and dynamics of protein folding intermediates and both partly folded and unfolded states.

Mechanisms and uses of hydrogen exchange

Recent work has largely completed our understanding of the hydrogen-exchange chemistry of unstructured proteins and nucleic acids. Some of the high-energy structural fluctuations that determine the hydrogen-exchange behavior of native macromolecules have been explained; others remain elusive. A growing number of applications are exploiting hydrogen-exchange behavior to study difficult molecular systems and elicit otherwise inaccessible information on protein structure, dynamics and energetics.

Automatic recognition of hydrophobic clusters and their correlation with protein folding units

A method is described to objectively identify hydrophobic clusters in proteins of known structure. Clusters are found by examining a protein for compact groupings of side chains. Compact clusters contain seven or more residues, have an average of 65% hydrophobic residues, and usually occur in protein interiors. Although smaller clusters contain only side-chain moieties, larger clusters enclose significant portions of the peptide backbone in regular secondary structure. These clusters agree well with hydrophobic regions assigned by more intuitive methods and many larger clusters correlate with protein domains. These results are in striking contrast with the clustering algorithm of J. Heringa and P. Argos (1991, J Mol Biol 220:151-171). That method finds that clusters located on a protein's surface are not especially hydrophobic and average only 3-4 residues in size. Hydrophobic clusters can be correlated with experimental evidence on early folding intermediates. This correlation is optimized when clusters with less than nine hydrophobic residues are removed from the data set. This suggests that hydrophobic clusters are important in the folding process only if they have enough hydrophobic residues.

Properties of the scrapie prion protein: quantitative analysis of protease resistance

The disease-specific isoform of the prion protein (PrPSc) is an essential part of the infectious particle which causes spongiform degeneration in various mammalian species. PrPSc differs from PrP of normal animals (PrPc) by its relative protease resistance. The physical nature of this difference is still unknown. We analyzed the protease resistance of PrPSc quantitatively using an enzyme-linked immunofiltration assay. PrPSc was rendered completely protease-sensitive at alkaline pH or in > 1.5 M guanidinium thiocyanate (GdnSCN). Denaturation in 4 M GdnSCN completely abolished the protease resistance of PrPSc within 15 min, while denaturation in 7.2 M urea showed a slower time course. In the presence of ethanol, PrPSc was protected from denaturation by GdnSCN or alkaline pH. Denaturation curves were used to calculate the free energy (delta GD) as a function of different denaturant concentrations. Linear regression of delta GD values was used to extrapolate the free energy in the absence of denaturants (delta GH2O), yielding similar values (delta GH2O,GdnSCN = -2.3 kcal/mol; delta GH2O,urea = -3.1 kcal/mol). The linear relationship between delta GD and the denaturant concentration is suggestive of a two-state model involving the conformational change of a single protein domain. This is also reflected in the small number of side chains (11.6) additionally exposed to the solvent upon conversion of PrPSc to its protease-sensitive isoform. Our results suggest that only minor rearrangements of the structure of PrP are needed to abolish the protease resistance of PrPSc.

Thermodynamics of ubiquitin unfolding

The energetics of ubiquitin unfolding have been studied using differential scanning microcalorimetry. For the first time it has been shown directly that the enthalpy of protein unfolding is a nonlinear function of temperature. Thermodynamic parameters of ubiquitin unfolding were correlated with the structure of the protein. The enthalpy of hydrogen bonding in ubiquitin was calculated and compared to that obtained for other proteins. It appears that the energy of hydrogen bonding correlates with the average length of the hydrogen bond in a given protein structure.

Equilibrium folding studies of cellular retinoic acid binding protein, a predominantly beta-sheet protein

We have examined the conformational behavior under various unfolding conditions of a predominantly beta-sheet protein, cellular retinoic acid binding protein (CRABP). Urea unfolding-refolding of CRABP is a highly cooperative process that can be approximated by a two-state model. Acid denaturation is also cooperative and reversible and leads to a state containing nonnative residual structure: Below pH 2.6, CRABP contains a substantially larger amount of alpha-helix than under native conditions. CRABP adopts up to 75% alpha-helix in solutions containing a high percentage of 2,2,2-trifluoroethanol. The acid-denatured state of CRABP undergoes a conformational change to a state containing predominantly beta-sheet structure upon the addition of small amounts of Na2SO4. This conformational malleability may be important for the folding mechanism of CRABP. The possible implication of nonnative alpha-helical structure in the folding of CRABP is discussed.

NMR, alcohols, protein solvation and protein denaturation

The amino acid side chains of Ser, Thr and Tyr contain the hydroxyl functional group that is characteristic of alcohols, and a variety of different alcohols are used either as structure-inducing or denaturing co-solvents in physical-chemical studies of proteins. This article surveys selected aspects of the influence of protein-intrinsic and -extrinsic hydroxyl groups on studies of the solvation of peptides and proteins by nuclear magnetic resonance (NMR) spectroscopy.

Heteronuclear three-dimensional NMR spectroscopy of a partially denatured protein: the A-state of human ubiquitin

Human ubiquitin is a 76-residue protein that serves as a protein degradation signal when conjugated to another protein. Ubiquitin has been shown to exist in at least three states: native (N-state), unfolded (U-state), and, when dissolved in 60% methanol:40% water at pH 2.0, partially folded (A-state). If the A-state represents an intermediate in the folding pathway of ubiquitin, comparison of the known structure of the N-state with that of the A-state may lead to an understanding of the folding pathway. Insights into the structural basis for ubiquitin's role in protein degradation may also be obtained. To this end we determined the secondary structure of the A-state using heteronuclear three-dimensional NMR spectroscopy of uniformly 15N-enriched ubiquitin. Sequence-specific 1H and 15N resonance assignments were made for more than 90% of the residues in the A-state. The assignments were made by concerted analysis of three-dimensional 1H-15N NOESY-HMQC and TOCSY-HMQC data sets. Because of 1H chemical shift degeneracies, the increased resolution provided by the 15N dimension was critical. Analysis of short- and long-range NOEs indicated that only the first two strands of beta-sheet, comprising residues 2-17, remain in the A-state, compared to five strands in the N-state. NOEs indicative of an alpha-helix, comprising residues 25-33, were also identified. These residues were also helical in the N-state. In the N-state, residues in this helix were in contact with residues from the first two strands of beta-sheet. It is likely, therefore, that residues 1-33 comprise a folded domain in the A-state of ubiquitin. On the basis of 1H alpha chemical shifts and weak short-range NOEs, residues 34-76 do not adopt a rigid secondary structure but favor a helical conformation. This observation may be related to the helix-inducing effects of the methanol present. The secondary structure presented here differs from and is more thorough than that determined previously by two-dimensional 1H methods [Harding et al. (1991) Biochemistry, 30, 3120-3128].