Osmolyte perturbation reveals conformational equilibria in spin‐labeled proteins

Recent evidence suggests that proteins at equilibrium can exist in a manifold of conformational substates, and that these substates play important roles in protein function. Therefore, there is great interest in identifying regions in proteins that are in conformational exchange. Electron paramagnetic resonance spectra of spin‐labeled proteins containing the nitroxide side chain (R1) often consist of two (or more) components that may arise from slow exchange between conformational substates (lifetimes > 100 ns). However, crystal structures of proteins containing R1 have shown that multicomponent spectra can also arise from equilibria between rotamers of the side chain itself. In this report, it is shown that these scenarios can be distinguished by the response of the system to solvent perturbation with stabilizing osmolytes such as sucrose. Thus, site‐directed spin labeling (SDSL) emerges as a new tool to explore slow conformational exchange in proteins of arbitrary size, including membrane proteins in a native‐like environment. Moreover, equilibrium between substates with even modest differences in conformation is revealed, and the simplicity of the method makes it suitable for facile screening of multiple proteins. Together with previously developed strategies for monitoring picosecond to millisecond backbone dynamics, the results presented here expand the timescale over which SDSL can be used to explore protein flexibility.     

Structural heterogeneity in protein crystals

Extensive conformational heterogeneity is reported in highly refined crystallographic models for the proteins crambin, erabutoxin, myohemerythrin, and lamprey hemoglobin. From 6% to 13% of the amino acid side chains of these four proteins are seen in multiple, discrete conformations. Most common are flexible side chains on the molecular surface, but structural heterogeneity occasionally extends to buried side chains or to the polypeptide backbone. A few instances of sequence heterogeneity are also very clear. Numerous solvent sites are multiplets, and at high resolution, multiple, mutually exclusive solvent networks are observed. The proteins have been studied with X-ray diffraction data extending to spacings of from 0.945 to 2.0 A. The extensive heterogeneity observed here provides detailed, accurate structures for conformational substates of these molecules and sets a lower bound on the number of substates accessible to each protein molecule in solution. Electron density is missing or very weak for only a few side chains in these protein crystals, revealing a strong preference for discrete over continuous conformational perturbations. The results at very high resolution further suggest that even rather small conformational fluctuations produce discrete substates and that unresolved conformers are accommodated in increased atomic thermal parameters.     

Osmolytes modulate conformational exchange in solvent‐exposed regions of membrane proteins

Site‐directed spin labeling (SDSL) was used to investigate local structure and conformational exchange in two bacterial outer‐membrane TonB‐dependent transporters, BtuB and FecA. Protecting osmolytes, such as polyethylene glycols (PEGs) are known to modulate a substrate‐dependent conformational equilibrium in the energy coupling motif (Ton box) of BtuB. Here, we demonstrate that a segment that is N‐terminal to the Ton box in BtuB, is in conformational exchange between ordered and disordered states with or without substrate. Protecting osmolytes shift this equilibrium to favor the more ordered, folded state. However, a segment of BtuB that is C‐terminal to the Ton box that is not solvent exposed is insensitive to PEGs. Protecting osmolytes also modulate a conformational equilibrium in the Ton box of FecA, with larger molecular weight PEGs producing the largest shifts in the conformational free energy. These data indicate that solvent‐exposed regions of these transporters undergo conformational exchange and that regions of these transporters that are involved in protein–protein interactions sample multiple conformational substates. The sensitivity to solute provides an explanation for differences seen between two high‐resolution structures of BtuB, which each likely represent one conformation from a subset of states that are normally sampled by the protein. This work also illustrates how SDSL and osmolytes may be used to characterize and quantitate conformational equilibria in membrane proteins.     

Hydrogen bond dynamics in membrane protein function

Changes in inter-helical hydrogen bonding are associated with the conformational dynamics of membrane proteins. The function of the protein depends on the surrounding lipid membrane. Here we review through specific examples how dynamical hydrogen bonds can ensure an elegant and efficient mechanism of long-distance intra-protein and protein-lipid coupling, contributing to the stability of discrete protein conformational substates and to rapid propagation of structural perturbations. This article is part of a Special Issue entitled: Protein Folding in Membranes.     

PABMB Lecture. Protein dynamics, folding and misfolding: from basic physical chemistry to human conformational diseases.

Proteins exhibit a variety of motions ranging from amino acid side-chain rotations to the motions of large domains. Recognition of their conformational flexibility has led to the view that protein molecules undergo fast dynamic interconversion between different conformational substates. This proposal has received support from a wide variety of experimental techniques and from computer simulations of protein dynamics. More recently, studies of the subunit dissociation of oligomeric proteins induced by hydrostatic pressure have shown that the characteristic times for subunit exchange between oligomers and for interconversion between different conformations may be rather slow (hours or days). In such cases, proteins cannot be treated as an ensemble of rapidly interconverting conformational substates, but rather as a persistently heterogeneous population of different long-lived conformers. This is reminiscent of the deterministic behavior exhibited by macroscopic bodies, and may have important implications for our understanding of protein folding and biological functions. Here, we propose that the deterministic behavior of proteins may be closely related to the genesis of conformational diseases, a class of pathological conditions that includes transmissible spongiform encephalopathies, Alzheimer's disease and other amyloidosis.     

Hydrogen exchange monitored by MALDI-TOF mass spectrometry for rapid characterization of the stability and conformation of proteins

Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) has been used to monitor hydrogen exchange on entire proteins. Two alternative methods have been used to carry out the hydrogen exchange studies, exchanging deuteron (H to D experiments) or proton (D to H experiments). In the former case, the use of a deuterated matrix has made possible to overcome back-exchange problems and attain reproducible results. The methods presented have been used to determine the slow exchange core of the potato carboxypeptidase inhibitor in different folding states, and to differentially compare the activation domain of human procarboxypeptidase A2 versus three site-directed mutants of different conformational stability. In this work, we show that by using MALDI-TOF MS to monitor hydrogen exchange in entire proteins, it is possible to rapidly check the folding state of a protein and characterize mutational effects on protein conformation and stability, while requiring minimal amounts of sample.     

A statistical mechanical model for hydrogen exchange in globular proteins.

We develop a statistical mechanical theory for the mechanism of hydrogen exchange in globular proteins. Using the HP lattice model, we explore how the solvent accessibilities of chain monomers vary as proteins fluctuate from their stable native conformations. The model explains why hydrogen exchange appears to involve two mechanisms under different conditions of protein stability; (1) a “global unfolding” mechanism by which all protons exchange at a similar rate, approaching that of the denatured protein, and (2) a “stable-state” mechanism by which protons exchange at rates that can differ by many orders of magnitude. There has been some controversy about the stable-state mechanism: does exchange take place inside the protein by solvent penetration, or outside the protein by the local unfolding of a subregion? The present model indicates that the stable-state mechanism of exchange occurs through an ensemble of conformations, some of which may bear very little resemblance to the native structure. Although most fluctuations are small-amplitude motions involving solvent penetration or local unfolding, other fluctuations (the conformational distant relatives) can involve much larger transient excursions to completely different chain folds.     

Backbone dynamics and energetics of a calmodulin domain mutant exchanging between closed and open conformations.

Previous studies have suggested that the Ca2+-saturated E140Q mutant of the C-terminal domain of calmodulin exhibits equilibrium exchange between "open" and "closed" conformations similar to those of the Ca2+-free and Ca2+-saturated states of wild-type calmodulin. The backbone dynamics of this mutant were studied using15N spin relaxation experiments at three different temperatures. Measurements at each temperature of the15N rate constants for longitudinal and transverse auto-relaxation, longitudinal and transverse cross-correlation relaxation, and the1H-15N cross-relaxation afforded unequivocal identification of conformational exchange processes on microsecond to millisecond time-scales, and characterization of fast fluctuations on picosecond to nanosecond time-scales using model-free approaches. The results show that essentially all residues of the protein are involved in conformational exchange. Generalized order parameters of the fast internal motions indicate that the conformational substates are well folded, and exclude the possibility that the exchange involves a significant population of unfolded or disordered species. The temperature dependence of the order parameters offers qualitative estimates of the contribution to the heat capacity from fast fluctuations of the protein backbone, revealing significant variation between the well-ordered secondary structure elements and the more flexible regions. The temperature dependence of the conformational exchange contributions to the transverse auto-relaxation rate constants directly demonstrates that the microscopic exchange rate constants are greater than 2.7x10(3)s-1at 291 K. The conformational exchange contributions correlate with the chemical shift differences between the Ca2+-free and Ca2+-saturated states of the wild-type protein, thereby substantiating that the conformational substates are similar to the open and closed states of wild-type calmodulin. Taking the wild-type chemical shifts to represent the conformational substates of the mutant and populations estimated previously, the microscopic exchange rate constants could be estimated as 2x10(4)to 3x10(4)s-1at 291 K for a subset of residues. The temperature depen dence of the exchange allows the characterization of apparent energy barriers of the conformational transition, with results suggesting a complex process that does not correspond to a single global transition between substates.     

Conformational substates of myoglobin detected by extrinsic dynamic fluorescence studies

The extent of conformational substates of two apomyoglobins, i.e., sperm whale and tuna apomyoglobin, was investigated by examining the fluorescence decay in the frequency domain of the extrinsic fluorophore TNS [6-(p-toluidino)-2-naphthalenesulfonic acid] bound to the heme binding site. Data analysis was performed in terms of a continuous, unimodal lifetime distribution having a Lorentzian shape. The results were compared with those for the free fluorophore in an isotropic nonviscous solvent. The incorporation of TNS into the protein matrix resulted in a broadening of the lifetime distribution due to the microenvironmental heterogeneity generated by structural fluctuations. The larger width of lifetime distribution observed for TNS bound to tuna apomyoglobin was related to a more extended conformational space accessible to the fluorophore in this protein compared to sperm whale myoglobin. A temperature increase from 15 to 40 degrees C produced a further broadening of the lifetime distributions of TNS bound to both proteins. This result can be explained by assuming the existence of conformational substates at high energy content or separated by high energy barriers, which are not populated at low temperature. The overall picture emerging from the reported data is that the lifetime distributions of TNS bound to apomyoglobins are determined largely by the number of conformational substates accessible to the protein matrix and, to a lesser extent, by the interconversion rates among these states.     

Hierarchical Conformational Analysis of Native Lysozyme Based on Sub-Millisecond Molecular Dynamics Simulations

Hierarchical organization of free energy landscape (FEL) for native globular proteins has been widely accepted by the biophysics community. However, FEL of native proteins is usually projected onto one or a few dimensions. Here we generated collectively 0.2 milli-second molecular dynamics simulation trajectories in explicit solvent for hen egg white lysozyme (HEWL), and carried out detailed conformational analysis based on backbone torsional degrees of freedom (DOF). Our results demonstrated that at micro-second and coarser temporal resolutions, FEL of HEWL exhibits hub-like topology with crystal structures occupying the dominant structural ensemble that serves as the hub of conformational transitions. However, at 100ns and finer temporal resolutions, conformational substates of HEWL exhibit network-like topology, crystal structures are associated with kinetic traps that are important but not dominant ensembles. Backbone torsional state transitions on time scales ranging from nanoseconds to beyond microseconds were found to be associated with various types of molecular interactions. Even at nanoseconds temporal resolution, the number of conformational substates that are of statistical significance is quite limited. These observations suggest that detailed analysis of conformational substates at multiple temporal resolutions is both important and feasible. Transition state ensembles among various conformational substates at microsecond temporal resolution were observed to be considerably disordered. Life times of these transition state ensembles are found to be nearly independent of the time scales of the participating torsional DOFs.     

Reduced temperature dependence of collective conformational opening in a hyperthermophile rubredoxin.

Spatially localized differences in the conformational dynamics of the rubredoxins from the hyperthermophile Pyrococcus furiosus (Pf) and the mesophile Clostridium pasteurianum (Cp) are monitored via amide exchange measurements. As shown previously for the hyperthermophile protein, nearly all backbone amides of the Cp rubredoxin exhibit EX(2) hydrogen exchange kinetics with conformational opening rates of >1 s(-)(1). Significantly slower amide exchange is observed for Pf rubredoxin in the region surrounding the metal site and the proximal end of the three-stranded beta-sheet, while for the rest of the structure, the exchange rates at 23 degrees C are similar for both proteins. For the multiple-turn region comprising residues 14-32 in both rubredoxins, the uniformity of both the exchange rate constants and the values of the activation energy at the slowly exchanging sites is consistent with a model of solvent exposure via a subglobal cooperative conformational opening. In contrast to the common expectation of increased rigidity in the hyperthermophile proteins, below room temperature Pf rubredoxin exhibits a larger apparent flexibility in this multiple-turn region. The smaller enthalpy for the conformational opening process of this region in Pf rubredoxin reflects the much weaker temperature dependence of the underlying conformational equilibrium in the hyperthermophile protein.     

Biological insights from hydrogen exchange mass spectrometry

Over the past two decades, hydrogen exchange mass spectrometry (HXMS) has achieved the status of a widespread and routine approach in the structural biology toolbox. The ability of hydrogen exchange to detect a range of protein dynamics coupled with the accessibility of mass spectrometry to mixtures and large complexes at low concentrations result in an unmatched tool for investigating proteins challenging to many other structural techniques. Recent advances in methodology and data analysis are helping HXMS deliver on its potential to uncover the connection between conformation, dynamics and the biological function of proteins and complexes. This review provides a brief overview of the HXMS method and focuses on four recent reports to highlight applications that monitor structure and dynamics of proteins and complexes, track protein folding, and map the thermodynamics and kinetics of protein unfolding at equilibrium. These case studies illustrate typical data, analysis and results for each application and demonstrate a range of biological systems for which the interpretation of HXMS in terms of structure and conformational parameters provides unique insights into function. This article is part of a Special Issue entitled: Mass spectrometry in structural biology.     

Microsecond exchange of internal water molecules in bacteriorhodopsin

The proton-conducting pathway of bacteriorhodopsin (BR) contains at least nine internal water molecules that are thought to be key players in the proton translocation mechanism. Here, we report the results of a multinuclear (1H, 2H, 17O) magnetic relaxation dispersion (MRD) study with the primary goal of determining the rate of exchange of these internal water molecules with bulk water. This rate is of interest in current attempts to elucidate the molecular details of the proton translocation mechanism. The relevance of water exchange kinetics is underscored by recent crystallographic findings of substantial variations in the number and locations of internal water molecules during the photocycle. Moreover, internal water exchange is believed to be governed by conformational fluctuations in the protein and can therefore provide information about the thermal accessibility of functionally important conformational substates. The present 2H and 17O MRD data show that at least seven water molecules, or more if they are orientationally disordered, in BR have residence times (inverse exchange rate constant) in the range 0.1-10 micros at 277 K. At least five of these water molecules have residence times in the more restrictive range 0.1-0.5 micros. These results show that most or all of the deeply buried water molecules in BR exchange on a time-scale that is short compared to the rate-limiting step in the photocycle. The MRD measurements were performed on BR solubilized in micelles of octyl glucoside. From the MRD data, the rotational correlation time of detergent-solubilized BR was determined to 35 ns at 300 K, consistent with a monomeric protein in complex with about 150 detergent molecules. The solubilized protein was found to be stable in the dark for at least eight months at 277 K.     

Peptide Modulation of Class I Major Histocompatibility Complex Protein Molecular Flexibility and the Implications for Immune Recognition

T cells use the αβ T cell receptor (TCR) to recognize antigenic peptides presented by class I major histocompatibility complex proteins (pMHCs) on the surfaces of antigen-presenting cells. Flexibility in both TCRs and peptides plays an important role in antigen recognition and discrimination. Less clear is the role of flexibility in the MHC protein; although recent observations have indicated that mobility in the MHC can impact TCR recognition in a peptide-dependent fashion, the extent of this behavior is unknown. Here, using hydrogen/deuterium exchange, fluorescence anisotropy, and structural analyses, we show that the flexibility of the peptide binding groove of the class I MHC protein HLA-A*0201 varies significantly with different peptides. The variations extend throughout the binding groove, impacting regions contacted by TCRs as well as other activating and inhibitory receptors of the immune system. Our results are consistent with statistical mechanical models of protein structure and dynamics, in which the binding of different peptides alters the populations and exchange kinetics of substates in the MHC conformational ensemble. Altered MHC flexibility will influence receptor engagement, impacting conformational adaptations, entropic penalties associated with receptor recognition, and the populations of binding-competent states. Our results highlight a previously unrecognized aspect of the “altered self” mechanism of immune recognition and have implications for specificity, cross-reactivity, and antigenicity in cellular immunity.     

Solvent exchange of buried water and hydrogen exchange of peptide NH groups hydrogen bonded to buried waters in bovine pancreatic trypsin inhibitor

Solvent exchange of 18O-labeled buried water in bovine pancreatic trypsin inhibitor (BPTI), trypsin, and trypsin-BPTI complex is measured by high-precision isotope ratio mass spectrometry. Buried water is labeled by equilibration of the protein in 18O-enriched water. Protein samples are then rapidly dialyzed against water of normal isotope composition by gel filtration and stored. The exchangeable 18O label eluting with the protein in 10-300 s is determined by an H2O-CO2 equilibration technique. Exchange of buried waters with solvent water is complete before 10-15 s in BPTI, trypsin, and BPTI-trypsin, as well as in lysozyme and carboxypeptidase measured as controls. When in-exchange dialysis and storage are carried out at pH greater than or equal to 2.5, trypsin-BPTI and trypsin, but not free BPTI, have the equivalent of one 18O atom that exchanges slowly (after 300 s and before several days). This oxygen is probably covalently bound to a specific site in trypsin. When in-exchange dialysis and storage are carried out at pH 1.1, the equivalent of three to seven 18O atoms per molecule is associated with the trypsin-BPTI complex, apparently due to nonspecific covalent 18O labeling of carboxyl groups at low pH. In addition to 18O exchange of buried waters, the hydrogen isotope exchange of buried NH groups H bonded to buried waters was also measured. Their base-catalyzed exchange rate constants are on the order of NH groups that in the crystal are exposed to solvent (static accessibility greater than 0) and hydrogen-bonded main chain O, and their pH min is similar to that for model compounds. The pH dependence of their exchange rate constants suggests that direct exchange with water may significantly contribute to their observed exchange rate.     

Unraveling the dynamics of protein interactions with quantitative mass spectrometry

Knowledge of structure and dynamics of proteins and protein complexes is important to unveil the molecular basis and mechanisms involved in most biological processes. Protein complex dynamics can be defined as the changes in the composition of a protein complex during a cellular process. Protein dynamics can be defined as conformational changes in a protein during enzyme activation, for example, when a protein binds to a ligand or when a protein binds to another protein. Mass spectrometry (MS) combined with affinity purification has become the analytical tool of choice for mapping protein-protein interaction networks and the recent developments in the quantitative proteomics field has made it possible to identify dynamically interacting proteins. Furthermore, hydrogen/deuterium exchange MS is emerging as a powerful technique to study structure and conformational dynamics of proteins or protein assemblies in solution. Methods have been developed and applied for the identification of transient and/or weak dynamic interaction partners and for the analysis of conformational dynamics of proteins or protein complexes. This review is an overview of existing and recent developments in studying the overall dynamics of in vivo protein interaction networks and protein complexes using MS-based methods.     

Redox-dependent conformational selection in a Cys4Fe2S2 ferredoxin

Putidaredoxin (Pdx), a Cys4Fe2S2 ferredoxin from Pseudomonas putida, exhibits redox-dependent binding to its physiological redox partner, cytochrome P450(cam) (CYP101), with the reduced form of Pdx (Pdx(r)) binding with greater affinity to oxidized camphor-bound CYP101 than the oxidized form, Pdx(o). It has been previously shown that Pdx(o) is more dynamic than Pdx(r) on all accessible time scales, and it has been proposed that Pdx(r) samples only a fraction of the conformational substates populated by Pdx(o) on a time average. It is postulated that the ensemble subset populated by Pdx(r) is the same subset that binds CYP101, providing a mechanism for coupling the Pdx oxidation state to binding affinity for CYP101. Evidence from a variety of sources, including redox-dependent shifts of 15N and 13C resonances, indicates that the metal cluster binding loop of Pdx is the primary determinant of redox-dependent conformational selection. Patterns of paramagnetic effects suggest that the metal cluster binding loop contracts around the metal cluster upon reduction, possibly due to the strengthening of hydrogen bonds between the sulfur atoms of the metal cluster and the surrounding polypeptide NH and OH groups. Effects of this perturbation are then transmitted mechanically to other affected regions of the protein. A specific mutation has been introduced into the metal binding loop of Pdx, G40N, that slows conformational exchange sufficiently that the ensemble of conformational substates in Pdx(o) are directly observable as severe broadenings or splittings in affected NMR resonances. Many of the residues most affected by the mutation also show significant exchange contributions to 15N T(2) relaxation in wild-type Pdx(o). As predicted, G40N Pdx(r) shows a collapse of many of these multiplets and broadened lines to form much sharper resonances that are essentially identical to those observed in wild-type Pdx(r), indicating that Pdx(r) occupies fewer conformational substates than does Pdx(o). This is the first direct observation of such redox-dependent ensembles at slow exchange on the chemical shift time scale. These results confirm that conformational selection within the Fe2S2 cluster binding loop is the primary source of redox-dependent changes in protein dynamics in Pdx.     

Solvent accessibility in folded proteins. Studies of hydrogen exchange in trypsin.

In a native protein, the exchange of a peptide amide proton with solvent occurs by one of two pathways, either directly from the folded protein, or via unfolding, exchange taking place from the unfolded protein. From the thermal unfolding rate constants, the contribution of unfolding to the over-all kinetics as a function of solvent and temperature has been determined. Exchange involving unfolding of the protein is characterized by a high activation energy, in the range of 50 to 60 Cal per mol. The activiation energy (Eapp) of the rates of exchange directly from the folded protein is approximately 20 to 25 Cal per mol. Because for the proton transfer step, Eapp approximately equal to 20 Cal per mol, the activation energy for any contributing protein conformational process(es) is approximately equal to 0 to 5 Cal per mol. Most, if not all, of the peptide amide protons in a folded protein can exchange directly with solvent without the protein unfolding. The number of "slowly" exchanging protons at a given condition of pH and temperature is not related to a discrete structural unit, but rather to the distribution of observed rates within the broader distribution of actual rates. The large attenuation of hydrogen exchange rates in folded proteins, resulting in a distribution of first order rates over 6 orders of magnitude, is primarily due to the effects of restricted solvent accessibility of labile protons in the three-dimensional structure. Any protein conformational process, such as protein fluctuations, invoked to explain the solvent accessibility must be of low activation energy and attenuated by ethanol and other co-solvents (Woodward, C. K., Ellis, L. M., and Rosenberg, A. (1974) J. Biol. Chem. 250, 440-444).     

Correlation between calculated local stability and hydrogen exchange rates in proteins

The attempt is made to find new correlations between local structural characteristics of proteins and the hydrogen exchange rates of their individual main-chain amides, and to relate such correlations to possible mechanisms of hydrogen exchange. It is found that in bovine pancreatic trypsin inhibitor (BPTI) the surface area buried by a particular residue and its neighbors correlates with the exchange rate of the main-chain amide of that residue. As the area buried by a particular fragment can be associated with the stabilization of the protein structure by this fragment, the correlation suggests a role for the energetics of the local unfolding in the mechanism of hydrogen exchange. Calculations based on the assumption that the exchange mechanism involves local unfolding lead to quantitative agreement between the calculated and experimentally measured exchange rates for 80% of the amides of BPTI that are buried or hydrogen bonded to the main-chain or to internal water molecules. The same degree of correlation is found between the calculated exchange rates and partial exchange data for ribonuclease S, hen lysozyme and cytochrome c. A similarly strong correlation is found between calculated exchange rates and the exchange rates of ribonuclease A determined by neutron diffraction in the crystal. The criteria of correlation are, however, less stringent in this case because of the experimental errors, which are larger than for solution data. It is suggested that the observed correlation be used for predictions of hydrogen exchange rates in proteins.     

A partially buried site in homologous HPr proteins is not optimized for stability

The energetic consequences of site-specific replacement of a residue at a partially buried site in the two homologous HPr proteins from Escherichia coli and Bacillus subtilis is described. We determined previously that the replacement of a partially buried Lys residue with Glu at position 49 in E.coli HPr increased the conformational stability of the protein substantially because the side-chain of the latter residue could act as a hydrogen-bond acceptor. Here, we extend this analysis to other side-chains with different chemical properties and abilities to form hydrogen bonds to compare the properties of this position in the backgrounds of two different homologous HPr proteins. We find that the variants with polar residues that can form a tertiary hydrogen bond with a nearby site in the protein are more stable than either hydrophobic residues or polar residues that become buried yet are incapable of forming a new hydrogen bond. Furthermore, the protein with the wild-type residue in each HPr variant is not among the most stable of the proteins studied. These results suggest a general strategy for designing variants in which the overall stability of a protein can be modulated in a defined fashion.