Quantitation of protein-protein interactions by thermal stability shift analysis

Thermal stability shift analysis is a powerful method for examining binding interactions in proteins. We demonstrate that under certain circumstances, protein-protein interactions can be quantitated by monitoring shifts in thermal stability using thermodynamic models and data analysis methods presented in this work. This method relies on the determination of protein stabilities from thermal unfolding experiments using fluorescent dyes such as SYPRO Orange that report on protein denaturation. Data collection is rapid and straightforward using readily available real-time polymerase chain reaction instrumentation. We present an approach for the analysis of the unfolding transitions corresponding to each partner to extract the affinity of the interaction between the proteins. This method does not require the construction of a titration series that brackets the dissociation constant. In thermal shift experiments, protein stability data are obtained at different temperatures according to the affinity- and concentration-dependent shifts in unfolding transition midpoints. Treatment of the temperature dependence of affinity is, therefore, intrinsic to this method and is developed in this study. We used the interaction between maltose-binding protein (MBP) and a thermostable synthetic ankyrin repeat protein (Off7) as an experimental test case because their unfolding transitions overlap minimally. We found that MBP is significantly stabilized by Off7. High experimental throughput is enabled by sample parallelization, and the ability to extract quantitative binding information at a single partner concentration. In a single experiment, we were able to quantify the affinities of a series of alanine mutants, covering a wide range of affinities (～ 100 nM to ～ 100 μM).     

Autophosphorylation dependent destabilization of the insulin receptor kinase domain: tryptophan-1175 reports changes in the catalytic cleft

Protein kinases are regulated by conformational or chemical changes which facilitate access of substrates to the active site and promote correct orientations of catalytically essential residues and water molecules. The switch between basal and activated states of the insulin receptor's kinase domain (IRKD) results from autophosphorylation. We investigated the effects of IRKD autophosphorylation on the conformational stability by guanidine hydrochloride (GdnHCl) dependent denaturation and by iodide quenching of intrinsic fluorescence. Tryptophan residues of the recombinant soluble IRKD (residues R953-S1355) were excited at a lambdaex of 295 nm, and emission spectra were analyzed for centroid (a characteristic of average polarity of the indole rings' environments) and integrated fluorescence intensity over the lambdaem range of 310-420 nm. Denaturation profiles of both apo- and phospho-IRKD forms are complex with at least three distinct unfolding transitions. The first and last transitions were reversible and cooperative and had midpoints at 0.4 or 0.7 M GdnHCl and 2.4 or 2.7 M GdnHCl, respectively; transitions of phospho-IRKD occurred at lower GdnHCl concentrations. Calculations of free energy of unfolding suggested a loss of approximately 2.3 kcal/mol of stabilization for the first transition and approximately 1.5 kcal/mol for the third transition. Circular dichroism showed subtle changes in secondary structure over the first transition and global unfolding over the last transition. The first transition reports changes primarily in the local environment of W1175, which is near the catalytic loop and is conserved among protein tyrosine kinases. W1175 is also the dominant fluorophore of the native emission spectrum. Iodide quenching of W1175 was virtually undetectable in the apo-IRKD but significant in the phospho-IRKD, suggesting that W1175 exposure to small solutes is strongly dependent on the conformation of the activation loop. These studies indicate that autophosphorylation, while exposing the catalytic center, also produces a conformer less stable than the apoenzyme.     

Conformational transitions of adenylate kinase: switching by cracking

Conformational heterogeneity in proteins is known to often be the key to their function. We present a coarse grained model to explore the interplay between protein structure, folding and function which is applicable to allosteric or non-allosteric proteins. We employ the model to study the detailed mechanism of the reversible conformational transition of Adenylate Kinase (AKE) between the open to the closed conformation, a reaction that is crucial to the protein's catalytic function. We directly observe high strain energy which appears to be correlated with localized unfolding during the functional transition. This work also demonstrates that competing native interactions from the open and closed form can account for the large conformational transitions in AKE. We further characterize the conformational transitions with a new measure Phi(Func), and demonstrate that local unfolding may be due, in part, to competing intra-protein interactions.     

Localization of the structural change induced in tRNA fMET (Escherichia coli) by acidic pH.

We have compared the molecular mechanism of thermal unfolding for native tRNA fMet (Escherichia coli) and the denatured species produced by annealing at pH 4.3. Relaxation kinetic measurements reveal that the transitions assigned to melting of TphiC, anticodon, and acceptor stem helices at neutral pH remain essentially unaltered at pH 4.3, but the transition corresponding to coupled melting of tertiary structure and dihydrouridine helix is greatly affected. The Tm of this region is more than 20 degrees higher at pH 4.3 and it has a larger enthalpy formation than in the native state. The transition dynamics are also considerably changed. In contrast to the native structure, tRNA fMet1 and tRNA fMet3 have similar tertiary structure stabilities at pH 4.3. We conclude that the structural difference between native and acid-denatured forms is localized in the tertiary structure-dihydrouridine helix cooperative interaction region of the molecule.     

Fourier transform infrared spectroscopy suggests unfolding of loop structures precedes complete unfolding of pig citrate synthase

Pig citrate synthase (PCS) can be used as a model enzyme to gain some insight into the structural basis of protein thermostability. The thermal unfolding characteristics of the specific secondary structure elements within PCS were monitored in detail by following changes in its amide I band components. The result of our study indicates that PCS undergoes irreversible thermal denaturation. Detailed analysis reveals that the different secondary structures display a multistep transition with a major and a minor transition at different temperatures and a very small initial transition at the same temperature (30 degrees C). A plot of temperature-induced changes in (1)H-(2)H exchange, the decrease in the absorbance of the alpha-helical structures, and the increase in the absorbance of aggregated structures all have in common a multistep transition, the minor one centered at 45 degrees C and the major one around 59 degrees C. In contrast, a band that is tentatively assigned to loop structures displays these same minor and major transitions but at lower temperatures (39 and 52 degrees C, respectively). The transition, which occurs at 39-45 degrees C, is not associated with the appearance of aggregated structures. This transition may reflect a change in the tertiary structure of the protein. However, the final transition, which occurs at a higher temperature (52-59 degrees C), reflects unfolding and aggregation of the polypeptide chains. The Fourier transform infrared (FTIR) analysis suggests that PCS has a thermolabile region that unfolds first, some 7 degrees C below the main unfolding of the protein. We propose that this reflects the unfolding of the highly flexible loop segments, which in turn triggers the unfolding of the predominantly helical core structure of PCS.     

Analyzing large‐scale structural change in proteins: Comparison of principal component projection and sammon mapping

Effective analysis of large-scale conformational transitions in macromolecules requires transforming them into a lower dimensional representation that captures the dominant motions. Herein, we apply and compare two different dimensionality reduction techniques, namely, principal component analysis (PCA), a linear method, and Sammon mapping, which is nonlinear. The two methods are used to analyze four different protein transition pathways of varying complexity, obtained by using either the conjugate peak refinement method or constrained molecular dynamics. For the return-stroke in myosin, both Sammon mapping and PCA show that the conformational change is dominated by a simple rotation of a rigid body. Also, in the case of the T-->R transition in hemoglobin, both methods are able to identify the two main quaternary transition events. In contrast, in the cases of the unfolding transition of staphylococcal nuclease or the signaling switch of Ras p21, which are both more complex conformational transitions, only Sammon mapping is able to identify the distinct phases of motion.     

Non-linear qualitative signal processing for biological systems: application to the algal growth in bioreactors

We present in this paper a qualitative method to validate and monitor the structure of a non-linear model with respect to experimental data, under some hypotheses. This method is broadly independent of the analytical formulation of the model, and depends only on the qualitative structure (the signs of the Jacobian matrix). The temporal sequences of the extrema of a filtered experimental signal are compared with the transitions allowed by a graph. In particular, we show that the usual moving average of the outputs follows this transition graph. We apply this method to compare models of algal growth in a bioreactor with experimental data.     

Folding and stability of trp aporepressor from Escherichia coli

Equilibrium and kinetic studies of the urea-induced unfolding of trp aporepressor from Escherichia coli were performed to probe the folding mechanism of this intertwined, dimeric protein. The equilibrium unfolding transitions at pH 7.6 and 25 degrees C monitored by difference absorbance, fluorescence, and circular dichroism spectroscopy are coincident within experimental error. All three transitions are well described by a two-state model involving the native dimer and the unfolded monomer; the free energy of folding in the absence of denaturant and under standard-state conditions is estimated to be 23.3 +/- 0.9 kcal/mol of dimer. The midpoint of the equilibrium unfolding transition increases with increasing protein concentration in the manner expected from the law of mass action for the two-state model. We find no evidence for stable folding intermediates. Kinetic studies reveal that unfolding is governed by a single first-order reaction whose relaxation time decreases exponentially with increasing urea concentration and also decreases with increasing protein concentration in the transition zone. Refolding involves at least three phases that depend on both the protein concentration and the final urea concentration in a complex manner. The relaxation time of the slowest of these refolding phases is identical with that for the single phase in unfolding in the transition zone, consistent with the results expected for a reaction that is kinetically reversible. The two faster refolding phases are presumed to arise from slow isomerization reactions in the unfolded form and reflect parallel folding channels.     

Urea equilibrium unfolding of the major core protein of the retrovirus feline immunodeficiency virus and its tryptophan mutants

Circular dichroism and fluorescence spectroscopy have been employed to study the urea unfolding mechanism of a recombinant form of the major core protein of feline immunodeficiency virus (FIV-rp24) and its native tryptophan mutants. The equilibrium denaturation curves indicate the existence of two transitions. The first unfolding transition most likely reflects the denaturation of the carboxy-terminal region of FIV-rp24. Consequently, the second transition, where the changes in fluorescence are produced, should reflect the denaturation of the amino-terminal region. If the intermediate observed upon urea denaturation is an on-pathway species, the data described herein can reflect the sequential and independent loss of structure of the two domains that this type of proteins possesses.     

Calculation of mutational free energy changes in transition states for protein folding

Recent advances in experimental and computational methods have made it possible to determine with considerable accuracy the structures whose formation is rate limiting for the folding of some small proteins-the transition state ensemble, or TSE. We present a method to analyze and validate all-atom models of such structures. The method is based on the comparison of experimental data with the computation of the change in free energy of the TSE resulting from specific mutations. Each mutation is modeled individually in all members of an ensemble of transition state structures using a method originally developed to predict mutational changes in the stability of native proteins. We first apply this method to six proteins for which we have determined the TSEs with a technique that uses experimental mutational data (Phi-values) as restraints in the structure determination and find a highly significant correlation between the calculated free energy changes and those derived from experimental kinetic data. We then use the procedure to analyze transition state structures determined by molecular dynamics simulations of unfolding, again finding a high correlation. Finally, we use the method to estimate changes in folding rates of several hydrophobic core mutants of Fyn SH3. Taken together, these results show that the procedure developed here is a tool of general validity for analyzing, assessing, and improving the quality of the structures of transition states for protein folding.     

The stability, structural organization, and denaturation of pectate lyase C, a parallel beta-helix protein

Pectate lyase C (pelC) was the first protein in which the parallel beta-helix structure was recognized. The unique features of parallel beta-helix-containing proteins-a relatively simple topology and unusual interactions among side chains-make pelC an interesting protein to study with respect to protein folding. In this paper, we report studies of the unfolding equilibrium of pelC. PelC is unfolded reversibly by gdn-HCl at pH 7 and 5, as monitored by far- and near-UV CD and fluorescence. The coincidence of these spectroscopically detected transitions is consistent with a two-state transition at pH 7, but the three probes are not coincident at pH 5. No evidence was found for a loosely folded intermediate in the transition region at pH 5. At pH 7, the for unfolding is 12.2 kcal/mol, with the midpoint of the transition at 0.99 M gdn-HCl and m = 12.3 kcal/(mol.M). Thus, pelC is unusually stable and has an m value that is much larger than for typical globular proteins. Thermal denaturation of pelC has been studied by differential scanning calorimetry (DSC) and by CD. Although thermal denaturation is not reversible, valid thermodynamic data can be obtained for the unfolding transition. DeltaH(van't Hoff)/DeltaH(cal) is less than 1 for pHs between 5 and 8, with a maximum value of 0.91 at pH 7 decreasing to 0.85 at pH 8 and to 0.68 at pH 5. At all pHs studied, the excess heat capacity can be deconvoluted into two components corresponding to two-state transitions that are nearly coincident at pH 7, but deviate more at higher and lower pH. Thus, pelC appears to consist of two domains that interact strongly and unfold in a cooperative fashion at pH 7, but the cooperativity decreases at higher and lower pH. The crystal structure of pelC shows no obvious domain structure, however.     

Unfolding of the cold shock protein studied with biased molecular dynamics

The cold shock protein from Bacillus caldolyticus is a small beta-barrel protein that folds in a two-state mechanism. For the native protein and for several mutants, a wealth of experimental data are available on stability and folding, so that it is an optimal system to study this process. We compare data from unfolding simulations (trajectories of 5 and up to 12 ns) obtained with a bias potential at room temperature and from unbiased thermal unfolding simulations with experimental data. The unfolding patterns derived from the trajectories starting from different native-like conformations and subject to different unfolding conditions agree. The transition state found in the simulations of unfolding is close to the native structure in agreement with experiment. Moreover, a lower value of the free energy barrier of unfolding was found for the mutant R3E than for the mutant E46A and the native protein, as indicated by experimental data. The first unfolding event involves the three-stranded beta-sheet whose decomposition corresponds to the transition state. In contrast to conclusions drawn from experiments, we found that the two-stranded beta-strand forms the most stable substructure, which decomposes very late in the unfolding process. However, assuming that this structure forms very early in the folding process, our findings would not contradict the experiments but require a different interpretation of them.     

Residue-level NMR view of the urea-driven equilibrium folding transition of SUMO-1 (1-97): native preferences do not increase monotonously

SUMO-1 (1-97) is a crucial protein in the machinery of post-translational modifications. We observed by circular dichroism and fluorescence spectroscopy that urea-induced unfolding of this protein is a complex process with the possibility of occurrence of detectable intermediates along the way. The tertiary structure is completely lost around approximately 4.5 M urea with a transition mid-point at 2.53 M urea, while the secondary structure unfolding seems to show two transitions, with mid-points at 2.42 M and 5.69 M urea. We have elucidated by systematic urea titration, the equilibrium residue level structural and dynamics changes along the entire folding/unfolding transition by multidimensional NMR. With urea dilution, the protein is seen to progressively lose most of the broad beta-domain structural preferences present at 8 M urea, acquire some helical propensities at 5 M urea, and lose some of them again on further dilution of urea. Between 3 M and 2 M urea, the protein starts afresh to acquire native structural features. These observations are contrary to the conventional notion that proteins fold with monotonously increasing native-type preferences. For folding below approximately 3 M urea, the region around the alpha1 helix appears to be a potential folding initiation site. The folding seems to start with a collapse into native-like topologies, at least in parts, and is followed by formation of secondary and tertiary structure, perhaps by cooperative rearrangements. The motional characteristics of the protein show sequence-dependent variation as the concentration of urea is progressively reduced. At the sub-nanosecond level, the features are extremely unusual for denatured states, and only certain segments corresponding to the flexible regions in the native protein display these motions at the different concentrations of urea.     

Theoretical analysis of Lumry-Eyring models in differential scanning calorimetry

A theoretical analysis of several protein denaturation models (Lumry-Eyring models) that include a rate-limited step leading to an irreversibly denatured state of the protein (the final state) has been carried out. The differential scanning calorimetry transitions predicted for these models can be broadly classified into four groups: situations A, B, C, and C′. (A) The transition is calorimetrically irreversible but the rate-limited, irreversible step takes place with significant rate only at temperatures slightly above those corresponding to the transition. Equilibrium thermodynamics analysis is permissible. (B) The transition is distorted by the occurrence of the rate-limited step; nevertheless, it contains thermodynamic information about the reversible unfolding of the protein, which could be obtained upon the appropriate data treatment. (C) The heat absorption is entirely determined by the kinetics of formation of the final state and no thermodynamic information can be extracted from the calorimetric transition; the rate-determining step is the irreversible process itself. (C′) same as C, but, in this case, the rate-determining step is a previous step in the unfolding pathway. It is shown that ligand and protein concentration effects on transitions corresponding to situation C (strongly rate-limited transitions) are similar to those predicted by equilibrium thermodynamics for simple reversible unfolding models. It has been widely held in recent literature that experimentally observed ligand and protein concentration effects support the applicability of equilibrium thermodynamics to irreversible protein denaturation. The theoretical analysis reported here disfavors this claim.     

Thermal unfolding of myosin rod and light meromyosin: circular dichroism and tryptophan fluorescence studies

Rabbit skeletal myosin rod, which is the coiled-coil alpha-helical portion of myosin, contains two tryptophan residues located in the light meromyosin (LMM) portion whose fluorescence contributes 27% to the fluorescence of the entire myosin molecule. The temperature dependence of several fluorescence parameters (quantum yield, spectral position, polarization) of the rod and its LMM portion was compared to the thermal unfolding of the helix measured with circular dichroism. Rod unfolds with three major helix unfolding transitions: at 43, 47, and 53 degrees C, with the 43 and 53 degrees C transitions mainly located in the LMM region and the 47 degrees C transition mainly located in the subfragment 2 region. The fluorescence study showed that the 43 degrees C transition does not involve the tryptophan-containing region and that the 47 degrees C transition produces an intermediate with different fluorescence properties from both the completely helical and fully unfolded states. That is, although the fluorescence of the 47 degrees C intermediate is markedly quenched, the tryptophyl residues do not become appreciably exposed to solvent until the 53 degrees C transition. It is suggested that although the intermediate that is formed in the 47 degrees C transition contains an extensive region which is devoid of alpha-helix, the unfolded region is not appreciably solvated or flexible. It appears to have the properties of a collapsed nonhelical state rather than a classical random coil.     

Unassisted refolding of urea unfolded rhodanese

In vitro refolding after urea unfolding of the enzyme rhodanese (thiosulfate:cyanide sulfurtransferase, EC 2.8.1.1) normally requires the assistance of detergents or chaperonin proteins. No efficient, unassisted, reversible unfolding/folding transition has been demonstrated to date. The detergents or the chaperonin proteins have been proposed to stabilize folding intermediates that kinetically limit folding by aggregating. Based on this hypothesis, we have investigated a number of experimental conditions and have developed a protocol for refolding, without assistants, that gives evidence of a reversible unfolding transition and leads to greater than 80% recovery of native enzyme. In addition to low protein concentration (10 micrograms/ml), low temperatures are required to maximize refolding. Otherwise optimal conditions give less than 10% refolding at 37 degrees C, whereas at 10 degrees C the recovery approaches 80%. The unfolding/refolding phases of the transition curves are most similar in the region of the transition, and refolding yields are significantly reduced when unfolded rhodanese is diluted to low urea concentrations, rather than to concentrations near the transition region. This is consistent with the formation of "sticky" intermediates that can remain soluble close to the transition region. Apparently, nonnative structures, e.g. aggregates, can form rapidly at low denaturant concentrations, and their subsequent conversion to the native structure is slow.     

Local and long-range interactions in the thermal unfolding transition of bovine pancreatic ribonuclease A

This research was undertaken to distinguish between local and global unfolding in the reversible thermal denaturation of bovine pancreatic ribonclease A (RNase A). Local unfolding was monitored by steady-state and time-resolved fluorescence of nine mutants in each of which a single tryptophan was substituted for a wild-type residue. Global unfolding was monitored by far-UV circular dichroism and UV absorbance. All the mutants (except F8W and D38W) exhibited high specific enzymatic activity, and their far-UV CD spectra were very close to that of wild-type RNase A, indicating that the tryptophan substitutions did not affect the structure of any of the mutants (excluding K1W and Y92W) under folding conditions at 20 degrees C. Like wild-type RNase A, the various mutants exhibited reversible cooperative thermal unfolding transitions at pH 5, with transition temperatures 2.5-11 degrees C lower than that of the wild-type transition, as detected by far-UV CD or UV absorbance. Even at 80 degrees C, well above the cooperative transition of all the RNase A mutants, a considerable amount of secondary and tertiary structure was maintained. These studies suggest the following two-stage mechanism for the thermal unfolding transition of RNase A as the temperature is increased. First, at temperatures lower than those of the main cooperative transition, long-range interactions within the major hydrophobic core are weakened, e.g., those involving residues Phe-8 (in the N-terminal helix) and Lys-104 and Tyr-115 (in the C-terminal beta-hairpin motif). The structure of the chain-reversal loop (residues 91-95) relaxes in the same temperature range. Second, the subsequent higher-temperature cooperative unfolding transition is associated with a loss of secondary structure and additional changes in the tertiary contacts of the major hydrophobic core, e.g., those involving residues Tyr-73, Tyr-76, and Asp-38 on the other side of the molecule. The hydrophobic interactions of the C-terminal loop of the protein are enhanced by high temperature, and perhaps are responsible for the preservation of the local structural environment of Trp-124 at temperatures slightly above the major cooperative transition. The results shed new light on the thermal unfolding transitions, generally supporting the thermal unfolding hypothesis of Burgess and Scheraga, as modified by Matheson and Scheraga.     

Denaturation of replication protein A reveals an alternative conformation with intact domain structure and oligonucleotide binding activity

Replication protein A (RPA) is a heterotrimeric, multidomain, single-stranded DNA-binding protein. Using spectroscopic methods and methylene carbene-based chemical modification methods, we have identified conformational intermediates in the denaturation pathway of RPA. Intrinsic protein fluorescence studies reveal unfolding profiles composed of multiple transitions, with midpoints at 1.5, 2.7, 4.2, and 5.3 M urea. CD profiles of RPA unfolding are characterized by a single transition. RPA is stabilized with respect to the CD-monitored transition when bound to a dA15 oligonucleotide. However, oligonucleotide binding appears to exert little, if any, effect on the first fluorescence transition. Methylene carbene chemical modification, coupled with MALDI-TOF mass spectrometry analysis, was also used to monitor unfolding of several specific RPA folds of the protein. The unfolding profiles of the individual structures are characterized by single transitions similar to the CD-monitored transition. Each fold, however, unravels with different individual characteristics, suggesting significant autonomy. Based on results from chemical modification and spectroscopic analyses, we conclude the initial transition observed in fluorescence experiments represents a change in the juxtaposition of binding folds with little unraveling of the domain structures. The second transition represents the unfolding of the majority of fold structure, and the third transition observed by fluorescence correlates with the dissociation of the 70- and 32-kD subunits.     

Spectrin R16: Broad energy barrier or sequential transition states?

A number of models have been proposed to account for nonlinearity in the relation between observed rate constants for folding and/or unfolding and denaturant concentration. Where curvature is seen principally in the arm of a chevron plot, three explanations are proposed: a change in the ground state at increasing concentration of urea, movement of the transition state along a broad energy barrier, and a switch between two sequential transition states separated by an on-pathway high-energy intermediate. Here we demonstrate that the latter two models in particular can be used to describe the data for the all-alpha protein spectrin R16. Further, whatever the method of analysis, the pattern of Phi-values seen is robust; thus we would draw the same conclusions from our data set independently of the method used for analysis. While this is not a novel observation, this is the first systematic study where a comparison has been made between Phi-values calculated using the broad and sequential transition state models.