Investigating protein unfolding kinetics by pulse proteolysis

Investigation of protein unfolding kinetics of proteins in crude samples may provide many exciting opportunities to study protein energetics under unconventional conditions. As an effort to develop a method with this capability, we employed “pulse proteolysis” to investigate protein unfolding kinetics. Pulse proteolysis has been shown to be an effective and facile method to determine global stability of proteins by exploiting the difference in proteolytic susceptibilities between folded and unfolded proteins. Electrophoretic separation after proteolysis allows monitoring protein unfolding without protein purification. We employed pulse proteolysis to determine unfolding kinetics of E. coli maltose binding protein (MBP) and E. coli ribonuclease H (RNase H). The unfolding kinetic constants determined by pulse proteolysis are in good agreement with those determined by circular dichroism. We then determined an unfolding kinetic constant of overexpressed MBP in a cell lysate. An accurate unfolding kinetic constant was successfully determined with the unpurified MBP. Also, we investigated the effect of ligand binding on unfolding kinetics of MBP using pulse proteolysis. On the basis of a kinetic model for unfolding of MBP•maltose complex, we have determined the dissociation equilibrium constant (K_d) of the complex from unfolding kinetic constants, which is also in good agreement with known K_d values of the complex. These results clearly demonstrate the feasibility and the accuracy of pulse proteolysis as a quantitative probe to investigate protein unfolding kinetics.     

Unfolding rates of globular proteins determined by kinetics of proteolysis

A convenient method for the determination of unfolding rates of small globular proteins under physiological conditions was developed using digestion with proteases. The apparent first-order rate constants for digestion of lysozyme with thermolysin and with Pronase at pH 8 and 50 degrees C were shown to be saturated with increases of concentrations of these proteases. The maximum rate constants extrapolated were identical in digestions with two different proteases, and were found to be equal to the unfolding rate constant of lysozyme. Similarly, the unfolding rate constant of RNase A at pH 8 and 50 degrees C, and those of lysozyme, RNase A and beta-lactoglobulin at pH 8 and 40 degrees C, were determined by the digestion method. Thus, it was shown that digestion by proteases proceeds mainly via the unfolded state of proteins.     

Probing the high energy states in proteins by proteolysis

Unless the native conformation has an unstructured region, proteases cannot effectively digest a protein under native conditions. Digestion must occur from a higher energy form, when at least some part of the protein is exposed to solvent and becomes accessible by proteases. Monitoring the kinetics and denaturant dependence of proteolysis under native conditions yields insight into the mechanism of proteolysis as well as these high-energy conformations. We propose here a generalized approach to exploit proteolysis as a tool to probe high-energy states in proteins. This "native state proteolysis" experiment was carried out on Escherichia coli ribonuclease HI. Mass spectrometry and N-terminal sequencing showed that thermolysin cleaves the peptide bond between Thr92 and Ala93 in an extended loop region of the protein. By comparing the proteolysis rate of the folded protein and a peptidic substrate mimicking the sequence at the cleavage site, the energy required to reach the susceptible state (Delta G(proteolysis)) was determined. From the denaturant dependence of Delta G(proteolysis), we determined that thermolysin digests this protein through a local fluctuation, i.e. localized unfolding with minimal change in solvent assessable surface area. Proteolytic susceptibilities of proteins are discussed based on the finding of this local fluctuation mechanism for proteolysis under native conditions.     

Local fluctuations vs. global unfolding of proteins investigated by limited proteolysis

The structure of a protein molecule consists of both rigid and flexible sections to satisfy the demands for stability and catalysis. Because the flexibility of a protein segment is indispensable for a proteolytic attack, limited proteolysis is a superb tool to analyse both confined local fluctuations and global unfolding events in proteins. While the identification of the primary cleavage products allows the assignment of the flexible regions to the primary structure, the kinetics of proteolytic degradation enables differentiation between local fluctuations in the native protein molecule and the global unfolding process during denaturation. Modifications of the amino acid sequence in the concerned regions can tune proteolytic susceptibility and alter protein stability. In the present paper, we summarise our results on native-state and unfolded-state proteolysis of ribonuclease A (RNase A) and the effect of mutations in the detected flexible regions on the stability and unfolding of the RNase A molecule.     

Local fluctuations vs. global unfolding of proteins investigated by limited proteolysis

The structure of a protein molecule consists of both rigid and flexible sections to satisfy the demands for stability and catalysis. Because the flexibility of a protein segment is indispensable for a proteolytic attack, limited proteolysis is a superb tool to analyse both confined local fluctuations and global unfolding events in proteins. While the identification of the primary cleavage products allows the assignment of the flexible regions to the primary structure, the kinetics of proteolytic degradation enables differentiation between local fluctuations in the native protein molecule and the global unfolding process during denaturation. Modifications of the amino acid sequence in the concerned regions can tune proteolytic susceptibility and alter protein stability. In the present paper, we summarise our results on native-state and unfolded-state proteolysis of ribonuclease A (RNase A) and the effect of mutations in the detected flexible regions on the stability and unfolding of the RNase A molecule.     

Mapping Transient Partial Unfolding by Protein Engineering and Native-State Proteolysis

Transient partial unfolding of proteins under native conditions may have significant consequences in the biochemical and biophysical properties of proteins. Native-state proteolysis offers a facile way to investigate the thermodynamic and kinetic accessibilities of partially unfolded forms (cleavable forms) under native conditions. However, determination of the structure of the cleavable form, which is populated only transiently, remains challenging. Although in some cases partially cleaved products from proteolysis provide information on the structure of this elusive form, proteolysis of many proteins does not accumulate detectable intermediates. Here, we describe a systematic approach to determining structures of cleavable forms by protein engineering and native-state proteolysis. By devising φ_c analysis, which is analogous to conventional φ analysis, we have determined the structure of the cleavable form of Escherichia coli maltose-binding protein (MBP), which does not accumulate any partially cleaved products. We mutated 10 buried residues in MBP to alanine and determined φ_c values from the effects of the mutations on global stability and proteolytic susceptibility. The result of this analysis suggests that two C-terminal helices in MBP are unfolded in their cleavable form. The effect of ligand binding on proteolytic susceptibility and C-terminal deletion mutations also confirms the proposed structure. Our approach and methodology are generally applicable not only in elucidating the mechanism of proteolysis but also in investigating other important processes involving partial unfolding under native conditions such as protein misfolding and aggregation.     

Probing Stability and Dynamics of Proteins by Protease Digestion I: Comparison of Protease Susceptibility and Thermal Stability of Cytochromes c

Abstract

Protease susceptibility of homologous proteins in their native conformations was studied. This work aims to establish a broad and quantitative basis for the utilization of protease digestion to analyze the local stability of native proteins. Using high-performance liquid chromatography (HPLC) the time course of the proteolytic degradation of intact proteins was quantitatively traced. Rapid separation of peptide fragments with HPLC made possible the elucidation of sequential digestion originating from the cleavage at a very few sites which are locally unstable in the protein structure. Using four serine proteases, chymotrypsin, trypsin, elastase and subtilisin BPN', we found some common trends in proteolysis for a group of proteins of the cytochrome c family. By comparing of the proteolysis and thermal denaturation with ten homologous cytochromes c extracted from horse, beef, Candida krusei, Saccharomyces cerevisiae, chicken, tuna, pigeon, rabbit, dog and rat, protease susceptibility was related to locally unfolding states intrinsic to the native conformation.     

Towards a consistent modeling of protein thermodynamic and kinetic cooperativity: how applicable is the transition state picture to folding and unfolding?

To what extent do general features of folding/unfolding kinetics of small globular proteins follow from their thermodynamic properties? To address this question, we investigate a new simplified protein chain model that embodies a cooperative interplay between local conformational preferences and hydrophobic burial. The present four-helix-bundle 55mer model exhibits protein-like calorimetric two-state cooperativity. It rationalizes native-state hydrogen exchange observations. Our analysis indicates that a coherent, self-consistent physical account of both the thermodynamic and kinetic properties of the model leads naturally to the concept of a native state ensemble that encompasses considerable conformational fluctuations. Such a multiple-conformation native state is seen to involve conformational states similar to those revealed by native-state hydrogen exchange. Many of these conformational states are predicted to lie below native baselines commonly used in interpreting calorimetric data. Folding and unfolding kinetics are studied under a range of intrachain interaction strengths as in experimental chevron plots. Kinetically determined transition midpoints match well with their thermodynamic counterparts. Kinetic relaxations are found to be essentially single-exponential over an extended range of model interaction strengths. This includes the entire unfolding regime and a significant part of a folding regime with a chevron rollover, as has been observed for real proteins that fold with non-two-state kinetics. The transition state picture of protein folding and unfolding is evaluated by comparing thermodynamic free energy profiles with actual kinetic rates. These analyses suggest that some chevron rollovers may arise from an internal frictional effect that increasingly impedes chain motions with more native conditions, rather than being caused by discrete deadtime folding intermediates or shifts of the transition state peak as previously posited.     

Imaging real-time proteolysis of single collagen I molecules with an atomic force microscope.

The dynamic process of synthesis and degradation of extracellular matrix molecules, including various collagens, is important in normal physiological functions and pathological conditions. Existing models of collagen enzymatic degradation reactions are derived from bulk biochemical assays. In this study, we have imaged in real-time individual collagen I molecules and their proteolysis by Clostridium histolyticum collagenases in phosphate-buffered saline (PBS) with atomic force microscopy (AFM). We have also imaged the likely binding and unbinding of collagenase molecules to single triple-helical collagen I molecules and subsequent proteolysis of subsets of the collagen molecules. The proteolysis of collagen molecules was inhibited by reduced calcium and acidification. Results from AFM study of collagen proteolysis are consistent with SDS-PAGE biochemical assays. The real-time proteolysis of single collagen I molecules followed simple Michaelis-Menton kinetics previously derived from bulk biochemical assays. This is the first report of imaging real-time proteolysis of single macromolecules and its inhibition on a molecular scale. A strong correspondence between the kinetics of proteolysis of single collagen molecules and the kinetics of proteolysis derived from bulk biochemical assays will have a wide applicability in examining real-time enzymatic reactions and their regulation at single molecule structural level. Such real-time study of single molecule proteolysis could provide a better understanding of the interactions between proteases and target proteins as well as proteases and protease inhibitors.     

Cytochrome c(553), a small heme protein that lacks misligation in its unfolded state, folds with rapid two-state kinetics

Cytochrome c(553) (cyt c(553)) from Desulfovibrio vulgaris is a small helical heme protein that displays apparent two-state equilibrium-unfolding behavior. The covalently attached heme is low-spin, ligated by Met and His residues, in the native state but becomes high-spin upon unfolding at pH 7. Here, we show that in contrast to other c-type heme proteins, where misligations in the unfolded states are prominent, cyt c(553) refolding kinetics at pH 7 proceeds rapidly without detectable intermediates. The extrapolated folding rate constant in water for oxidized cyt c(553) matches exactly that predicted from the cyt c(553) native-state topology: 5300 s(-1 )(experimental) versus 5020 s(-1) (predicted). We therefore conclude that the presence of the oxidized cofactor does not affect the intrinsic formation speed of the cyt c(553 )structural motif.     

Site-selective chemical protein glycosylation protects from autolysis and proteolytic degradation

Glycosylation is often cited as having a stabilizing effect upon proteins with respect to proteolysis, thermolysis and other forms of degradation. We present here a model study on an autolytic protease that has been chemically glycosylated to produce single glycoforms. The resulting glycosylated enzymes are more stable with respect to their own autolytic degradation and that by other proteases. Kinetic parameters for protease activity with respect to the degradation of small-molecule amide substrate reveal no significant change in inherent activity thereby suggesting that reduced autolysis and proteolysis are a consequence of stabilization, perhaps by steric blockade of cleavage points or alteration of local unfolding kinetics. Variation in glycan identity suggests that greater glycan size leads to greater stabilization.     

Effect of signal peptide on the stability and folding kinetics of maltose binding protein

While the role of the signal sequence in targeting proteins to specific subcellular compartments is well characterized, there are fewer studies that characterize its effects on the stability and folding kinetics of the protein. We report a detailed characterization of the folding kinetics and thermodynamic stabilities of maltose binding protein (MBP) and its precursor form, preMBP. Isothermal GdmCl and urea denaturation as a function of temperature and thermal denaturation studies have been carried out to compare stabilities of the two proteins. preMBP was found to be destabilized by about 2-6 kcal/mol (20-40%) with respect to MBP. Rapid cleavage of the signal peptide by various proteases shows that the signal peptide is accessible in the native form of preMBP. The observed rate constant of the major slow phase in folding was decreased 5-fold in preMBP relative to MBP. The rate constants of unfolding were similar at 25 degrees C, but preMBP also exhibited a large burst phase change in unfolding that was absent in MBP. At 10 degrees C, preMBP exhibited a higher unfolding rate than MBP as well as a large burst phase. The appreciable destabilization of MBP by signal peptide is functionally relevant, because it enhances the likelihood of finding the protein in an unfolded translocation-competent form and may influence the interactions of the protein with the translocation machinery. Destabilization is likely to result from favorable interactions between the hydrophobic signal peptide and other hydrophobic regions that are exposed in the unfolded state.     

Stability and unfolding of reduced Escherichia coli glutaredoxin 2: a monomeric structural homologue of the glutathione transferase family

Glutaredoxin 2 (Grx2) from Escherichia coli is monomeric and an atypical glutaredoxin that takes part in the monothiol deglutathionylation of proteins. Unlike its orthologs, Grx2 is a larger molecule with a canonical glutathione transferase (GST) fold that consists of two structurally distinct domains, an N-terminal glutaredoxin domain and a C-terminal alpha-helical domain. While GSTs are dimeric proteins, the conformational stability and unfolding kinetics of Grx2 were investigated to establish the contribution made by the domain interface to the stability of the tertiary structure of GST-like proteins without any influence from quaternary interactions. Equilibrium unfolding transitions for Grx2, using urea as a denaturant, are monophasic and exhibit coincidence of the fluorescence and CD data indicative of a concerted loss or formation of tertiary and secondary structure. The data fit well to a two-state N <--> U model with no evidence that an intermediate is being formed. The experimental m value [2.7 kcal mol (-1) (M urea) (-1)] is in excellent agreement with a predicted value of 2.5 kcal mol (-1) (M urea) (-1) based on the amount of surface area expected to become exposed during unfolding. These findings provide evidence that the two structurally distinct domains of Grx2 behave as a single cooperative folding unit. The unfolding kinetics are complex which, as a result of native-state heterogeneity, are characterized by two observable unfolding reactions that occur in parallel. A major population representing one distinct nativelike form unfolds on a fast track to denatured Grx2 with cis-Pro49. This is followed by a spectroscopically silent cis-trans proline isomerization reaction as determined by interrupted unfolding experiments. A minor population representing the other distinct nativelike form unfolds slowly with its rate being limited by an undetermined structural isomerization reaction. Further, there is no evidence indicating that unfolding proceeds via a high-energy intermediate that might suggest independent unfolding of the two nonidentical domains in Grx2. The kinetics data are, therefore, consistent with the existence of cooperativity between the domains, in agreement with the equilibrium data.     

Shear flow-induced detachment kinetics of Dictyostelium discoideum cells from solid substrate

Using Dictyostelium discoideum as a model organism of specific and nonspecific adhesion, we studied the kinetics of shear flow-induced cell detachment. For a given cell, detachment occurs for values of the applied hydrodynamic stress above a threshold. Cells are removed from the substrate with an apparent first-order rate constant that strongly depends on the applied stress. The threshold stress depends on cell size and physicochemical properties of the substrate, but is not affected by depolymerization of the actin and tubulin cytoskeleton. In contrast, the kinetics of cell detachment is almost independent of cell size, but is strongly affected by a modification of the substrate and the presence of an intact actin cytoskeleton. These results are interpreted in the framework of a peeling model. The threshold stress and the cell-detachment rate measure the local equilibrium energy and the dissociation rate constant of the adhesion bridges, respectively.     

Effects of sucrose on conformational equilibria and fluctuations within the native-state ensemble of proteins

Osmolytes increase the thermodynamic conformational stability of proteins, shifting the equilibrium between native and denatured states to favor the native state. However, their effects on conformational equilibria within native-state ensembles of proteins remain controversial. We investigated the effects of sucrose, a model osmolyte, on conformational equilibria and fluctuations within the native-state ensembles of bovine pancreatic ribonuclease A and S and horse heart cytochrome c. In the presence of sucrose, the far- and near-UV circular dichroism spectra of all three native proteins were slightly altered and indicated that the sugar shifted the native-state ensemble toward species with more ordered, compact conformations, without detectable changes in secondary structural contents. Thermodynamic stability of the proteins, as measured by guanidine HCl-induced unfolding, increased in proportion to sucrose concentration. Native-state hydrogen exchange (HX) studies monitored by infrared spectroscopy showed that addition of 1 M sucrose reduced average HX rate constants at all degrees of exchange of the proteins, for which comparison could be made in the presence and absence of sucrose. Sucrose also increased the exchange-resistant core regions of the proteins. A coupling factor analysis relating the free energy of HX to the free energy of unfolding showed that sucrose had greater effects on large-scale than on small-scale fluctuations. These results indicate that the presence of sucrose shifts the conformational equilibria toward the most compact protein species within native-state ensembles, which can be explained by preferential exclusion of sucrose from the protein surface.     

Increased proteolytic resistance of ribonuclease A by protein engineering.

Although highly stable toward unfolding, native ribonuclease A is known to be cleaved by unspecific proteases in the flexible loop region near Ala20. With the aim to create a protease-resistant ribonuclease A, Ala20 was substituted for Pro by site-directed mutagenesis. The resulting mutant enzyme was nearly identical to the wild-type enzyme in the near-UV and far-UV circular dichroism spectra, in its activity to 2',3'-cCMP and in its thermodynamic stability. However, the proteolytic resistance to proteinase K and subtilisin Carlsberg was extremely increased. Pseudo-first-order rate constants of proteolysis, determined by densitometric analysis of the bands of intact protein in SDS-PAGE, decreased by two orders of magnitude. In contrast, the rate constant of proteolysis with elastase was similar to that of the wild-type enzyme. These differences can be explained by the analysis of the fragments occurring in proteolysis with elastase. Ser21-Ser22 was identified as the main primary cleavage site in the degradation of the mutant enzyme by elastase. Obviously, this bond is not cleavable by proteinase K or subtilisin Carlsberg. The results demonstrate the high potential of a single mutation in protein stabilization to proteolytic degradation.     

Probing membrane protein unfolding with pulse proteolysis

Technical challenges have greatly impeded the investigation of membrane protein folding and unfolding. To develop a new tool that facilitates the study of membrane proteins, we tested pulse proteolysis as a probe for membrane protein unfolding. Pulse proteolysis is a method to monitor protein folding and unfolding, which exploits the significant difference in proteolytic susceptibility between folded and unfolded proteins. This method requires only a small amount of protein and, in many cases, may be used with unpurified proteins in cell lysates. To evaluate the effectiveness of pulse proteolysis as a probe for membrane protein unfolding, we chose Halobacterium halobium bacteriorhodopsin (bR) as a model system. The denaturation of bR in SDS has been investigated extensively by monitoring the change in the absorbance at 560 nm (A₅₆₀). In this work, we demonstrate that denaturation of bR by SDS results in a significant increase in its susceptibility to proteolysis by subtilisin. When pulse proteolysis was applied to bR incubated in varying concentrations of SDS, the remaining intact protein determined by electrophoresis shows a cooperative transition. The midpoint of the cooperative transition (C_m) shows excellent agreement with that determined by A₅₆₀. The C_m values determined by pulse proteolysis for M56A and Y57A bRs are also consistent with the measurements made by A₅₆₀. Our results suggest that pulse proteolysis is a quantitative tool to probe membrane protein unfolding. Combining pulse proteolysis with Western blotting may allow the investigation of membrane protein unfolding in situ without overexpression or purification.     

Characterization of the unfolding of ribonuclease a by a pulsed hydrogen exchange study: evidence for competing pathways for unfolding

The unfolding of ribonuclease A was studied in 5.2 M guanidine hydrochloride at pH 8 and 10 degrees C using multiple optical probes, native-state hydrogen exchange (HX), and pulse labeling by hydrogen exchange. First, native-state HX studies were used to demonstrate that the protein exists in two slowly interconverting forms under equilibrium native conditions: a predominant exchange-incompetent N form and an alternative ensemble of conformations, N(I), in which some amide hydrogens are fully exposed to exchange. Pulsed HX studies indicated that, during unfolding, the rates of exposure to exchange with solvent protons were similar for all backbone NH probe protons. It is shown that two parallel routes of unfolding are available to the predominant N conformation as soon as it encounters strong unfolding conditions. A fraction of molecules appears to rapidly form N(I) on one route. On the other route an exchange-incompetent intermediate state ensemble, I(U)(2), is formed. The kinetics of unfolding measured by far-UV circular dichroism (CD) were faster than those measured by near-UV CD and intrinsic tyrosine fluorescence of the protein. The logarithms of the rate constants of the unfolding reaction measured by all three optical probes also showed a nonlinear dependence on GdnHCl concentration. All of the data suggest that N(I) and I(U)(2) are nativelike in their secondary and tertiary structures. While N(I) unfolds directly to the fully exchange-competent unfolded state (U), I(U)(2) forms another intermediate I(U)(3) which then unfolds to U. I(U)(3) is devoid of all native alpha-helical secondary structure and has only 30% of the tertiary interactions still intact. Since the rates of global unfolding measured by near-UV CD and fluorescence agree well with the rates of exposure determined for all of the backbone NH probe protons, it appears that the rate-limiting step for the unfolding of RNase A is the dissolution of the entire native tertiary structure and penetration of water into the hydrophobic core.     

Vanadate inhibits the ATP-dependent degradation of proteins in reticulocytes without affecting ubiquitin conjugation

Reticulocytes contain a nonlysosomal, ATP-dependent system for degrading abnormal proteins and normal proteins during cell maturation. Vanadate, which inhibits several ATPases including the ATP-dependent proteases in Escherichia coli and liver mitochondria, also markedly reduced the ATP-dependent degradation of proteins in reticulocyte extracts. At low concentrations (K1 = 50 microM), vanadate inhibited the ATP-dependent hydrolysis of [3H]methylcasein and denatured 125I-labeled bovine serum albumin, but it did not reduce the low amount of proteolysis seen in the absence of ATP. This inhibition by vanadate was rapid in onset, reversed by dialysis, and was not mimicked by molybdate. Vanadate inhibits proteolysis at an ATP-stimulated step which is independent of the ATP requirement for ubiquitin conjugation to protein substrates. When the amino groups on casein and bovine serum albumin were covalently modified so as to prevent their conjugation to ubiquitin, the derivatized proteins were still degraded by an ATP-stimulated process that was inhibited by vanadate. In addition, vanadate did not reduce the ATP-dependent conjugation of 125I-ubiquitin to endogenous reticulocyte proteins, although it markedly inhibited their degradation. In intact reticulocytes vanadate also inhibited the degradation of endogenous proteins and of abnormal proteins containing amino acid analogs. This effect was rapid and reversible; however, vanadate also reduced protein synthesis and eventually lowered ATP levels in the intact cells. Vanadate (10 mM) has also been reported to decrease intralysosomal proteolysis in hepatocytes. However, in liver extracts this effect on lysosomal proteases required high concentrations of vanadate (K1 = 500 microM) and was also observed with molybdate, unlike the inhibition of ATP-dependent proteolysis in reticulocytes.     

The ligand affinity of proteins measured by isothermal denaturation kinetics

An isothermal denaturation kinetic method was developed for identifying potential ligands of proteins and measuring their affinity. The method is suitable for finding ligands specific toward proteins of unknown function and for large-scale drug screening. It consists of analyzing the kinetics of isothermal denaturation of the protein-with and without the presence of potential specific ligands-as measured by long-wavelength fluorescent dyes whose quantum yield increases when bound to hydrophobic regions exposed upon unfolding of the proteins. The experimental procedure was developed using thymidylate kinase and stromelysin as target proteins. The kinetics of thermal unfolding of both of these enzymes were consistent with a pathway of two consecutive first-order rate-limiting steps. Reflecting the stabilizing effect of protein/ligand complexes, the presence of specific ligands decreased the value of the rate constants of both steps in a dose-dependent manner. The dependence of the rate constants on ligand concentration obeyed a simple binding isotherm, the analysis of which yielded an accurate equilibrium constant for ligand binding. The method was validated by comparing its results with those obtained under the same conditions by steady-state fluorescence spectroscopy, circular dichroism, and uv spectrophotometry: The corresponding rate constants were comparable for each of the analytical detection methods.