Instrumentation for automated determination of protein stability

From radioactive ligands to biosensors: binding methods with olfactory proteins

In this paper, we critically review the binding protocols currently reported in the literature to measure the affinity of odorants and pheromones to soluble olfactory proteins, such as odorant-binding proteins (OBPs), chemosensory proteins (CSPs) and Niemann-Pick class C2 (NPC2) proteins. The first part contains a brief introduction on the principles of binding and a comparison of the techniques adopted or proposed so far, discussing advantages and problems of each technique, as well as their suitable application to soluble olfactory proteins. In the second part, we focus on the fluorescent binding assay, currently the most widely used approach. We analyse advantages and drawbacks, trying to identify the causes of anomalous behaviours that have been occasionally observed, and suggest how to interpret the experimental data when such events occur. In the last part, we describe the state of the art of biosensors for odorants, using soluble olfactory proteins immobilised on biochips, and discuss the possibility of using such approach as an alternative way to measure binding events and dissociation constants.

The fluorescence detected guanidine hydrochloride equilibrium denaturation of wild-type staphylococcal nuclease does not fit a three-state unfolding model

A three-state equilibrium unfolding of a protein can be difficult to detect if two of the states fail to differ in some easily measurable way. It has been unclear whether staphylococcal nuclease unfolds in a two-state fashion, with only the native and denatured states significantly populated at equilibrium, or in a three-state manner, with a well-populated intermediate. Since equilibrium unfolding experiments are commonly used to determine protein stability and the course of denaturation are followed by changes in the fluorescence which has difficulty in distinguishing various states, this is a potential problem for many proteins. Over the course of twenty years we have performed more than one hundred guanidine hydrochloride equilibrium denaturations of wild-type staphylococcal nuclease; to our knowledge, a number of denaturations unrivaled in any other protein system. A careful examination of the data from these experiments shows no sign of the behavior predicted by a three-state unfolding model. Specifically, a three-state unfolding should introduce a slight, but characteristic, non-linearity to the plot of stability versus denaturant concentration. The average residuals from this large number of repeated experiments do not show the predicted behavior, casting considerable doubt on the likelihood of a three-state unfolding for the wild-type protein. The methods used for analysis here could be applied to other protein systems to distinguish a two-state from a three-state denaturation.

Equilibrium and kinetic studies of protein cooperativity using urea-induced folding/unfolding of a Ubq-UIM fusion protein

Understanding the origins of cooperativity in proteins remains an important topic in protein folding. This study describes experimental folding/unfolding equilibrium and kinetic studies of the engineered protein Ubq-UIM, consisting of ubiquitin (Ubq) fused to the sequence of the ubiquitin interacting motif (UIM) via a short linker. Urea-induced folding/unfolding profiles of Ubq-UIM were monitored by far-UV circular dichroism and fluorescence spectroscopies and compared to those of the isolated Ubq domain. It was found that the equilibrium data for Ubq-UIM is inconsistent with a two-state model. Analysis of the kinetics of folding shows similarity in the folding transition state ensemble between Ubq and Ubq-UIM, suggesting that formation of Ubq domain is independent of UIM. The major contribution to the stabilization of Ubq-UIM, relative to Ubq, was found to be in the rates of unfolding. Moreover, it was found that the kinetic m-values for Ubq-UIM unfolding, monitored by different probes (far-UV circular dichroism and fluorescence spectroscopies), are different; thereby, further supporting deviations from a two-state behavior. A thermodynamic linkage model that involves four states was found to be applicable to the urea-induced unfolding of Ubq-UIM, which is in agreement with the previous temperature-induced unfolding study. The applicability of the model was further supported by site-directed variants of Ubq-UIM that have altered stabilities of Ubq/UIM interface and/or stabilities of individual Ubq- and UIM-domains. All variants show increased cooperativity and one variant, E43N_Ubq-UIM, appears to behave very close to an equilibrium two-state.

Protein stability by number: high-throughput and statistical approaches to one of protein science's most difficult problems

Most proteins are only barely stable, which impedes research, complicates therapeutic applications, and makes proteins susceptible to pathologically destabilizing mutations. Our ability to predict the thermodynamic consequences of even single point mutations is still surprisingly limited, and established methods of measuring stability are slow. Recent advances are bringing protein stability studies into the high-throughput realm. Some methods are based on inferential read-outs such as activity, proteolytic resistance or split-protein fragment reassembly. Other methods use miniaturization of direct measurements, such as intrinsic fluorescence, H/D exchange, cysteine reactivity, aggregation and hydrophobic dye binding (DSF). Protein engineering based on statistical analysis (consensus and correlated occurrences of amino acids) is promising, but much work remains to understand and implement these methods.

Protein stability in the presence of cosolutes

Protein scientists have long used cosolutes to study protein stability. While denaturants, such as urea, have been employed for a long time, the attention became focused more recently on protein stabilizers, including osmolytes. Here, we provide practical experimental instructions for the use of both stabilizing and denaturing osmolytes with proteins, as well as data evaluation strategies. We focus on protein stability in the presence of cosolutes and their mixtures at constant and variable temperature.

Refinement of noncalorimetric determination of the change in heat capacity, DeltaC(p), of protein unfolding and validation across a wide temperature range

The change in heat capacity, DeltaC(p), on protein unfolding has been usually determined by calorimetry. A noncalorimetric method which employs the Gibbs-Helmholtz relationship to determine DeltaC(p) has seen some use. Generally, in this method the free energy change on unfolding of the protein is determined at a variety of temperatures and the temperature at which DeltaG is zero, T(m), and change in enthalpy at T(m) are determined by thermal denaturation and DeltaC(p) is then calculated using the Gibbs-Helmholtz equation. We show here that an abbreviated method with stability determinations at just two temperatures gives values of DeltaC(p) consistent with values from free energy change on unfolding determination at a much wider range of temperatures. Further, even the free energy change on unfolding from a single solvent denaturation at the proper temperature, when coupled with the melting temperature, T(m), and the van't Hoff enthalpy, DeltaH(vH), from a thermal denaturation, gives a reasonable estimate of DeltaC(p), albeit with greater uncertainty than solvent denaturations at two temperatures. We also find that nonlinear regression of the Gibbs-Helmholtz equation as a function of stability and temperature while simultaneously fitting DeltaC(p), T(m), and DeltaH(vH) gives values for the last two parameters that are in excellent agreement with experimental values.

Effects of excluded volume upon protein stability in covalently cross-linked proteins with variable linker lengths

To explore the effects of molecular crowding and excluded volume upon protein stability, we used a series of cross-linking reagents with nine different single-cysteine mutants of staphylococcal nuclease to make covalently linked dimers. These cross-linkers ranged in length from 10.5 to 21.3 A, compelling separations which would normally be found only in the most concentrated protein solutions. The stabilities of the dimeric proteins and monomeric controls were determined by guanidine hydrochloride and thermal denaturation. Dimers with short linkers tend to exhibit pronounced three-state denaturation behavior, as opposed to the two-state behavior of the monomeric controls. Increasing linker length leads to less pronounced three-state behavior. The three-state behavior is interpreted in a three-state model where cross-linked native protein dimer, N-N, interconverts in a two-state transition with a dimer where one protein subunit is denatured, N-D. The remaining native protein in turn can denature in another two-state transition to a state, D-D, in which both tethered proteins are denatured. Three-state behavior is best explained by excluded volume effects in the denatured state. For many dimers, linkers longer than 17 A removed most three-state character. This sets a limit on the flexibility and size of the denatured state. Notably, in contradiction to theoretical predictions, these cross-linked dimers were not stabilized. The failure of these predictions is possibly due to neglect of the alteration in hydrophobic exposure that accompanies any significant reduction in the conformational space of the denatured state.

Thermal denaturations of staphylococcal nuclease wild-type and mutants monitored by fluorescence and circular dichroism are similar: lack of evidence for other than a two state thermal denaturation

It is unclear whether the thermal denaturation of staphylococcal nuclease is a two state, three state, or variable two state process. The thermal denaturation of wild-type staphylococcal nuclease was followed by tryptophan fluorescence and circular dichroism signal at 222 nm, forty-two and fourteen times, respectively. Analysis of this data using a simple two state model gave melting temperatures of 53.0+/-0.4 degrees C (fluorescence) and 52.7+/-0.6 degrees C (CD) and van't Hoff enthalpies of 82.4+/-2.6 kcal/mol and 88.6+/-4.2 kcal/mol. Ninety-seven mutants also had these parameters determined by both fluorescence and CD. The average difference between the melting temperatures was 1.05+/-0.75 degrees and the average difference between van't Hoff enthalpies was 1.6+/-4.8 kcal/mol. These very similar results for the two spectroscopic probes of structure are discussed in the context of the different models that have been proposed for nuclease denaturation. It is concluded, for most nuclease variants, that the errors introduced by a two state assumption are negligible and either virtually all helical structure is lost in any initial unfolding event or any intermediate must have low stability.

X-ray and thermodynamic studies of staphylococcal nuclease variants I92E and I92K: insights into polarity of the protein interior

We have used crystallography and thermodynamic analysis to study nuclease variants I92E and I92K, in which an ionizable side-chain is placed in the hydrophobic core of nuclease. We find that the energetic cost of burying ionizable groups is rather modest. The X-ray determinations show water molecules solvating the buried glutamic acid under cryo conditions, but not at room temperature. The lysine side-chain does not appear solvated in either case. Guanidine hydrochloride (GnHCl) denaturation of I92E and I92K, done as a function of pH and monitored by tryptophan fluorescence, showed that I92E and I92K are folded in the pH range pH 3.5-9.0 and pH 5.5-9.5, respectively. The stability of the parental protein is independent of pH over a broad range. In contrast, the stabilities of I92E and I92K exhibit a pH dependence, which is quantitatively explained by thermodynamic analysis: the PK(a) value of the buried K92 is 5.6, while that of the buried E92 is 8.65. The free energy difference between burying the uncharged and charged forms of the groups is modest, about 6 kcal/mol. We also found that epsilon(app) for I92K and I92E is in the range approximately 10-12, instead of 2-4 commonly used to represent the protein interior. Side-chains 92E and 92K were uncharged under the conditions of the X-ray experiment. Both are buried completely inside the well-defined hydrophobic core of the variant proteins without forming salt-bridges or hydrogen bonds to other functional groups of the proteins. Under cryo conditions 92E shows a chain of four water molecules, which hydrate one oxygen atom of the carboxyl group of the glutamic acid. Two other water molecules, which are present in the wild-type at all temperatures, are also connected to the water ring observed inside the hydrophobic core. The ready burial of water with an uncharged E92 raises the possibility that solvent excursions into the interior also take place in the wild-type protein, but in a random, dynamic way not detectable by crystallography. Such transient excursions could increase the average polarity, and thus epsilon(app), of the protein interior.

Mechanism of thermostabilization in a designed cold shock protein with optimized surface electrostatic interactions

Using computational and sequence analysis of bacterial cold shock proteins, we designed a protein (CspB-TB) that has the core residues of mesophilic protein from Bacillus subtilis(CspB-Bs) and altered distribution of surface charged residues. This designed protein was characterized by circular dichroism spectroscopy, and found to have secondary and tertiary structure similar to that of CspB-Bs. The activity of the CspB-TB protein as measured by the affinity to a single-stranded DNA (ssDNA) template at 25 degrees C is somewhat higher than that of CspB-Bs. Furthermore, the decrease in the apparent binding constant to ssDNA upon increase in temperature is much more pronounced for CspB-Bs than for CspB-TB. Temperature-induced unfolding (as monitored by differential scanning calorimetry and circular dichroism spectroscopy) and urea-induced unfolding experiments were used to compare the stabilities of CspB-Bs and CspB-TB. It was found that CspB-TB is approximately 20 degrees C more thermostable than CspB-Bs. The thermostabilization of CspB-TB relative to CspB-Bs is achieved by decrease in the enthalpy and entropy of unfolding without affecting their temperature dependencies, i.e. these proteins have similar heat capacity changes upon unfolding. These changes in the thermodynamic parameters result in the global stability function, i.e. Gibbs energy, deltaG(T), that is shifted to higher temperatures with only small changes in the maximum stability. Such a mechanism of thermostabilization, although predicted from the basic thermodynamic considerations, has never been identified experimentally.

Contribution of surface salt bridges to protein stability: guidelines for protein engineering

The small globular protein, ubiquitin, contains a pair of oppositely charged residues, K11 and E34, that according to the three-dimensional structure are located on the surface of this protein with a spatial orientation characteristic of a salt bridge. We investigated the strength of this salt bridge and its contribution to the global stability of the ubiquitin molecule. Using the "double mutant cycle" analysis, the strength of the pairwise interactions between K11 and E34 was estimated to be favorable by 3.6kJ/mol. Further, the salt bridge of the reverse orientation, i.e. E11/K34, can be formed and is found to have a strength (3.8kJ/mol) similar to that of the K11/E34 pair. However, the global stability of the K11/E34 variant of ubiquitin is 2.2kJ/mol higher than that of the E11/K34 variant. The difference in the contribution of the opposing salt bridge orientations to the overall stability of the ubiquitin molecule is attributed to the difference in the charge-charge interactions between residues forming the salt bridge and the rest of the ionizable groups in this protein. On the basis of these results, we concluded that surface salt bridges are stabilizing, but their contribution to the overall protein stability is strongly context-dependent, with charge-charge interactions being the largest determinant. Analysis of 16 salt bridges from six different proteins, for which detailed experimental data on energetics have been reported, support the conclusions made from the analysis of the salt bridge in ubiquitin. Implications of these findings for engineering proteins with enhanced thermostability are discussed.

Evidence for a monomeric intermediate in the reversible unfolding of F factor TraM

F factor TraM is essential for efficient bacterial conjugation, but its molecular function is not clear. Because the physical properties of TraM may provide clues to its role in conjugation, we have characterized the TraM oligomerization equilibrium. We show that the reversible unfolding transition is non-two-state, indicating the presence of at least one intermediate. Analytical ultracentrifugation experiments indicate that the first phase of unfolding involves dissociation of the tetramer into folded monomers, which are subsequently unfolded to the denatured state in the second phase. Furthermore, we show that a C-terminal domain isolated by limited proteolysis is tetrameric in solution, like the full-length protein, and that its loss of structure correlates with dissociation of the TraM tetramer. Unfolding of the individual domains indicates that the N- and C-terminal regions act cooperatively to stabilize the full-length protein. Together, these experiments suggest structural overlap of regions important for oligomerization and DNA binding. We propose that modulating the oligomerization equilibrium of TraM may regulate its essential activity in bacterial conjugation.

Changes in stability upon charge reversal and neutralization substitution in staphylococcal nuclease are dominated by favorable electrostatic effects

Single site mutations that reverse or neutralize a surface charge were made at 22 ionizable residues in staphylococcal nuclease. Unfolding free energies were obtained by guanidine hydrochloride denaturation. These data, in conjunction with previously obtained stabilities of the corresponding alanine mutants, unequivocally show that the dominant contribution to stability for virtually all of the wild-type side chains examined is the electrostatic effect associated with each residue's charged group. With only a few exceptions, these charges stabilize the native state, with an average loss of 0.5 kcal/mol of stability upon neutralization of a charge. When the charge is reversed, the average destabilization is doubled. Structure-based calculations of electrostatic free energy with the continuum method based on the finite difference solution to the linearized Poisson-Boltzmann equation reproduce the observed energetics when the polarizability in the protein interior is represented with a dielectric constant of 20. However, in some cases, large differences are found, giving insight into possible areas for improvement of the calculations. In particular, it appears that the assumptions made in the calculations about the absence of electrostatic interactions in the denatured state and the energetic consequences of dynamic fluctuations in the native state will have to be further explored.

Removal of surface charge-charge interactions from ubiquitin leaves the protein folded and very stable

The contribution of solvent-exposed charged residues to protein stability was evaluated using ubiquitin as a model protein. We combined site-directed mutagenesis and specific chemical modifications to first replace all Arg residues with Lys, followed by carbomylation of Lys-amino groups. Under the conditions in which all carboxylic groups are protonated (at pH 2), the chemically modified protein is folded and very stable (DeltaG = 18 kJ/mol). These results indicate that surface charge-charge interactions are not an essential fundamental force for protein folding and stability.

Comparing the effect on protein stability of methionine oxidation versus mutagenesis: steps toward engineering oxidative resistance in proteins

The biological activity of some proteins is known to be sensitive to oxidative damage caused by a variety of oxidants. The model protein staphylococcal nuclease was used to explore the effect on protein structural stability of oxidizing methionine to the sulfoxide form. These effects were compared with the effects of substituting methionines with isoleucine and leucine, a potential strategy for stabilizing proteins against oxidative damage. Wild-type nuclease and various mutants were oxidized with hydrogen peroxide. Stabilities of both oxidized and unoxidized proteins were determined by guanidine hydrochloride denaturation. Oxidation destabilized the wild-type protein by over 4 kcal/mol. This large loss of stability supports the idea that in some cases loss of biological activity is linked to disruption of the protein native state. Comparison of mutant protein's stability losses upon oxidation showed that methionines 65 and 98 had a much greater destabilizing effect when oxidized than methionines 26 or 32. While substitution of methionine 98 carried as great an energetic penalty as oxidation, substitution at position 65 was less disruptive than oxidation. Thus a simple substitution mutagenesis strategy to protect a protein against oxidative destabilization is practical for some methionine residues.

Increasing the thermostability of staphylococcal nuclease: implications for the origin of protein thermostability

Seven hyper-stable multiple mutants have been constructed in staphylococcal nuclease by various combinations of eight different stabilizing single mutants. The stabilities of these multiple mutants determined by guanidine hydrochloride denaturation were 3.4 to 5.6 kcal/mol higher than that of the wild-type. Their thermal denaturation midpoint temperatures were 12.6 to 22.9 deg. C higher than that of the wild-type. These are among the greatest increases in protein stability and thermal denaturation midpoint temperature relative to the wild-type yet attained. There has been great interest in understanding how proteins found in thermophilic organisms are stabilized. One frequently cited theory is that the packing of hydrophobic side-chains is improved in the cores of proteins isolated from thermophiles when compared to proteins from mesophiles. The crystal structures of four single and five multiple stabilizing mutants of staphylococcal nuclease were solved to high resolution. No large overall structural change was found, with most changes localized around the sites of mutation. Rearrangements were observed in the packing of side-chains in the major hydrophobic core, although none of the mutations was in the core. It is surprising that detailed structural analysis showed that packing had improved, with the volume of the mutant protein's hydrophobic cores decreasing as protein stability increased. Further, the number of van der Waals interactions in the entire protein showed an experimentally significant increase correlated with increasing stability. These results indicate that optimization of packing follows as a natural consequence of increased protein thermostability and that good packing is not necessarily the proximate cause of high stability. Another popular theory is that thermostable proteins have more electrostatic and hydrogen bonding interactions and these are responsible for the high stabilities. The mutants here show that increased numbers of electrostatic and hydrogen bonding interactions are not obligatory for large increases in protein stability.

High apparent dielectric constants in the interior of a protein reflect water penetration

A glutamic acid was buried in the hydrophobic core of staphylococcal nuclease by replacement of Val-66. Its pK(a) was measured with equilibrium thermodynamic methods. It was 4.3 units higher than the pK(a) of Glu in water. This increase was comparable to the DeltapK(a) of 4.9 units measured previously for a lysine buried at the same location. According to the Born formalism these DeltapK(a) are energetically equivalent to the transfer of a charged group from water to a medium of dielectric constant of 12. In contrast, the static dielectric constants of dry protein powders range from 2 to 4. In the crystallographic structure of the V66E mutant, a chain of water molecules was seen that hydrates the buried Glu-66 and links it with bulk solvent. The buried water molecules have never previously been detected in >20 structures of nuclease. The structure and the measured energetics constitute compelling and unprecedented experimental evidence that solvent penetration can contribute significantly to the high apparent polarizability inside proteins. To improve structure-based calculations of electrostatic effects with continuum methods, it will be necessary to learn to account quantitatively for the contributions by solvent penetration to dielectric effects in the protein interior.

Application of automated methods for determination of protein conformational stability

Automation of protein stability determination saves much time and often results in better data. In our laboratory we estimate that the lack of such instrumentation would increase the amount of time required to make and characterize a mutant protein by at least one-third. As with any experiment, lack of attention to important details can required repeating the experiment or, worse, give superficially good data that are flawed in some manner. Nevertheless, with a minimal investment of attention, an automated instrument can give dramatic gains in productivity.

Experimental measurement of the effective dielectric in the hydrophobic core of a protein

The dielectric inside a protein is a key physical determinant of the magnitude of electrostatic interactions in proteins. We have measured this dielectric phenomenologically, in terms of the dielectric that needs to be used with the Born equation in order to reproduce the observed pKa shifts induced by burial of an ionizable group in the hydrophobic core of a protein. Mutants of staphylococcal nuclease with a buried lysine residue at position 66 were engineered for this purpose. The pKa values of buried lysines were measured by difference potentiometry. The extent of coupling between the pKa and the global stability of the protein was evaluated by measuring pKa values in hyperstable forms of nuclease engineered to be 3.3 or 6.5 kcal mol-1 more stable than the wild type. The crystallographic structure of one mutant was determined to describe the environment of the buried lysine. The dielectrics that were measured range from 10 to 12. Published pKa values of buried ionizable residues in other proteins were analyzed in a similar fashion and the dielectrics obtained from these values are consistent with the ones measured in nuclease. These results argue strongly against the prevalent use of dielectrics of 4 or lower to describe the dielectric effect inside a protein in structure-based calculations of electrostatic energies with continuum dielectric models.

Chemically crosslinked protein dimers: stability and denaturation effects

Nine single substitution cysteine mutants of staphylococcal nuclease (nuclease) were preferentially crosslinked at the introduced cysteine residues using three different bifunctional crosslinking reagents; 1,6-bismaleimidohexane (BMH), 1,3-dibromo-2-propanol (DBP), and the chemical warfare agent, mustard gas (bis(2-chloroethyl)sulfide; mustard). BMH and mustard gas are highly specific reagents for cysteine residues, whereas DBP is not as specific. Guanidine hydrochloride (GuHCl) denaturations of the resulting dimeric proteins exhibited biphasic unfolding behavior that did not fit the two-state model of unfolding. The monofunctional reagent, epsilon-maleimidocaproic acid (MCA), was used as a control for the effects of alkylation. Proteins modified with MCA unfolded normally, showing that this unusual unfolding behavior is due to crosslinking. The data obtained from these crosslinked dimers was fitted to a three-state thermodynamic model of two successive transitions in which the individual subunits cooperatively unfold. These two unfolding transitions were very different from the unfolding of the monomeric protein. These differences in unfolding behavior can be attributed in large part to changes in the denatured state. In addition to GuHCl titrations, the crosslinked dimers were also thermally unfolded. In contrast to the GuHCl denaturations, analysis of this data fit a two-state model well, but with greatly elevated van't Hoff enthalpies in many cases. However, clear correlations between the thermal and GuHCl denaturations exist, and the differences in thermal unfolding can be rationalized by postulating interactions of the denatured crosslinked proteins.