Engineering of protein thermodynamic, kinetic, and colloidal stability: Chemical Glycosylation with monofunctionally activated glycans

In this work we establish the relationship between chemical glycosylation and protein thermodynamic, kinetic, and colloidal stability. While there have been reports in the literature that chemical glycosylation modulates protein stability, mechanistic details still remain uncertain. To address this issue, we designed and coupled monofunctional activated glycans (lactose and dextran) to the model protein alpha-chymotrypsin (alpha-CT). This resulted in a series of glycoconjugates with variations in the glycan size and degree of glycosylation. Thermodynamic unfolding, thermal inactivation, and temperature-induced aggregation experiments revealed that chemical glycosylation increased protein thermodynamic (Delta G(25 degrees C)), kinetic (t(1/2)(45 degrees C)), and colloidal stability. These results highlight the potential of chemical glycosylation with monofunctional activated glycans as a technology for increasing the long-term stability of liquid protein formulations for industrial and biotherapeutic applications.     

Effect of PEG modification on subtilisin Carlsberg activity, enantioselectivity, and structural dynamics in 1,4-dioxane

The employment of enzymes as catalysts within organic media has traditionally been hampered by the reduced enzymatic activities when compared to catalysis in aqueous solution. Although several complementary hypotheses have provided mechanistic insights into the causes of diminished activity, further development of biocatalysts would greatly benefit from effective chemical strategies (e.g., PEGylation) to ameliorate this event. Herein we explore the effects of altering the solvent composition from aqueous buffer to 1,4-dioxane on structural, dynamical, and catalytic properties of the model enzyme subtilisin Carlsberg (SBc). Furthermore, we also investigate the effects of dissolving the enzyme in 1,4-dioxane through chemical modification with poly(ethylene)-glycol (PEG, M(W) = 20 kDa) on these enzyme properties. In 1,4-dioxane a 10(4)-fold decrease in the enzyme's catalytic activity was observed for the hydrolysis reaction of vinyl butyrate with D(2)O and a 50% decrease in enzyme structural dynamics as evidenced by reduced amide H/D exchange kinetics occurred. Attaching increasing amounts of PEG to the enzyme reversed some of the activity loss. Evaluation of the structural dynamic behavior of the PEGylated enzyme within the organic solvent revealed an increase in structural dynamics at increased PEGylation. Correlation analysis between the catalytic and structural dynamic parameters revealed that the enzyme's catalytic activity and enantioselectivity depended on the changes in protein structural dynamics within 1,4-dioxane. These results demonstrate the importance of protein structural dynamics towards regulating the catalytic behavior of enzymes within organic media.     

Influence of modulated structural dynamics on the kinetics of alpha-chymotrypsin catalysis. Insights through chemical glycosylation, molecular dynamics and domain motion analysis

Although the chemical nature of the catalytic mechanism of the serine protease alpha-chymotrypsin (alpha-CT) is largely understood, the influence of the enzyme's structural dynamics on its catalysis remains uncertain. Here we investigate whether alpha-CT's structural dynamics directly influence the kinetics of enzyme catalysis. Chemical glycosylation [Solá RJ & Griebenow K (2006) FEBS Lett 580, 1685-1690] was used to generate a series of glycosylated alpha-CT conjugates with reduced structural dynamics, as determined from amide hydrogen/deuterium exchange kinetics (k(HX)). Determination of their catalytic behavior (K(S), k(2), and k(3)) for the hydrolysis of N-succinyl-Ala-Ala-Pro-Phe p-nitroanilide (Suc-Ala-Ala-Pro-Phe-pNA) revealed decreased kinetics for the catalytic steps (k(2) and k(3)) without affecting substrate binding (K(S)) at increasing glycosylation levels. Statistical correlation analysis between the catalytic (DeltaG( not equal)k(i)) and structurally dynamic (DeltaG(HX)) parameters determined revealed that the enzyme acylation and deacylation steps are directly influenced by the changes in protein structural dynamics. Molecular modelling of the alpha-CT glycoconjugates coupled with molecular dynamics simulations and domain motion analysis employing the Gaussian network model revealed structural insights into the relation between the protein's surface glycosylation, the resulting structural dynamic changes, and the influence of these on the enzyme's collective dynamics and catalytic residues. The experimental and theoretical results presented here not only provide fundamental insights concerning the influence of glycosylation on the protein biophysical properties but also support the hypothesis that for alpha-CT the global structural dynamics directly influence the kinetics of enzyme catalysis via mechanochemical coupling between domain motions and active site chemical groups.     

Rapid cooperative two-state folding of a miniature alpha-beta protein and design of a thermostable variant

There is a great deal of interest in developing small stably folded miniature proteins. A limited number of these molecules have been described, however they typically have not been characterized in depth. In particular, almost no detailed studies of the thermodynamics and folding kinetics of these proteins have been reported. Here we describe detailed studies of the thermodynamics and kinetics of folding of a 39 residue mixed alpha-beta protein (NTL9(1-39)) derived from the N-terminal domain of the ribosomal protein L9. The protein folds cooperatively and rapidly in a two-state fashion to a native state typical of those found for normal globular proteins. At pH 5.4 in 20mM sodium acetate, 100mM NaCl the temperature of maximum stability is 6 degrees C, the t(m) is 65.3 degrees C, deltaH degrees (t(m)) is between 24.6 kcalmol(-1) and 26.3 kcalmol(-1), and deltaC(p) degrees is 0.38 kcalmol(-1)deg(-1). The thermodynamic parameters are in the range expected on the basis of per residue values determined from databases of globular proteins. H/2H exchange measurements reveal a set of amides that exchange via global unfolding, exactly as expected for a normal cooperatively folded globular protein. Kinetic measurements show that folding is two-state folding. The folding rate is 640 s(-1) and the value of deltaG degrees calculated from the folding and unfolding rates is in excellent agreement with the equilibrium value. A designed thermostable variant, generated by mutating K12 to M, was characterized and found to have a t(m) of 82 degrees C. Equilibrium and kinetic measurements demonstrate that its folding is cooperative and two-state.     

Thermal stability of proteins does not correlate with the energy of intramolecular interactions

Small monomeric proteins from mesophilic and thermophilic organisms were studied. They have close structural and physical and chemical properties but vary in thermal stability. A thermodynamic analysis of heat unfolding was made and integral enthalpy of unfolding (DeltaH(unf)), heat capacity of hydration (DeltaC(p)(hyd)) and enthalpy of hydration (DeltaH(hyd)) and of the buried surface area (DeltaASA) of nonpolar and polar groups as well as the enthalpy of disruption of intramolecular interaction (DeltaH(int) in gas phase) at 298 K were determined. The absence of correlation between protein thermostability and energetic components suggests that regulatory mechanism of protein thermal stabilization has entropic nature.     

Amino acid substitutions affecting protein dynamics in eglin C do not affect heat capacity change upon unfolding

The heat capacity change upon unfolding (deltaC(p)) is a thermodynamic parameter that defines the temperature dependence of the thermodynamic stability of proteins; however, physical basis of the heat capacity change is not completely understood. Although empirical surface area-based calculations can predict heat capacity changes reasonably well, accumulating evidence suggests that changes in hydration of those surfaces is not the only parameter contributing to the observed heat capacity changes upon unfolding. Because packing density in the protein interior is similar to that observed in organic crystals, we hypothesized that changes in protein dynamics resulting in increased rigidity of the protein structure might contribute to the observed heat capacity change upon unfolding. Using differential scanning calorimetry we characterized the thermodynamic behavior of a serine protease inhibitor eglin C and two eglin C variants with altered native state dynamics, as determined by NMR. We found no evidence of changes in deltaC(p) in either of the variants, suggesting that changes in rigidity do not contribute to the heat capacity change upon unfolding in this model system.     

Measuring the stability of partly folded proteins using TMAO

Standard methods for measuring free energy of protein unfolding by chemical denaturation require complete folding at low concentrations of denaturant so that a native baseline can be observed. Alternatively, proteins that are completely unfolded in the absence of denaturant can be folded by addition of the osmolyte trimethylamine N-oxide (TMAO), and the unfolding free energy can then be calculated through analysis of the refolding transition. However, neither chemical denaturation nor osmolyte-induced refolding alone is sufficient to yield accurate thermodynamic unfolding parameters for partly folded proteins, because neither method produces both native and denatured baselines in a single transition. Here we combine urea denaturation and TMAO stabilization as a means to bring about baseline-resolved structural transitions in partly folded proteins. For Barnase and the Notch ankyrin domain, which both show two-state equilibrium unfolding, we found that DeltaG degrees for unfolding depends linearly on TMAO concentration, and that the sensitivity of DeltaG degrees to urea (the m-value) is TMAO independent. This second observation confirms that urea and TMAO exert independent effects on stability over the range of cosolvent concentrations required to bring about baseline-resolved structural transitions. Thermodynamic parameters calculated using a global fit that assumes additive, linear dependence of DeltaG degrees on each cosolvent are similar to those obtained by standard urea-induced unfolding in the absence of TMAO. Finally, we demonstrate the applicability of this method to measurement of the free energy of unfolding of a partly folded protein, a fragment of the full-length Notch ankyrin domain.     

In vitro refolding of PEGylated lipase

Covalent modification of proteins with polyethylene glycol (PEG) has become a well established drug enhancement strategy in the biopharmaceutical industry. The general benefits of PEGylation, such as prolonged serum half-lives or reduced in vivo immunogenicity, are well known. To date, the PEGylation process has been performed with purified proteins, which often requires additional multi-step purification steps to harvest the desired PEGylate. However, it would be beneficial for bioprocessing if 'renaturation,' i.e. in vitro refolding and 'modification,' and PEGylation can be integrated, especially for inclusion body proteins. We investigated the feasibility of protein PEGylation under denaturing conditions and of protein refolding with the attached PEG molecule. Using lipase as a model protein, PEGylation occurred in 8 M urea and covalently attached PEG did not appear to hinder subsequent refolding.     

Folding and association of an extremely stable dimeric protein from Sulfolobus islandicus

ORF56 is a plasmid-encoded protein from Sulfolobus islandicus, which probably controls the copy number of the pRN1 plasmid by binding to its own promotor. The protein showed an extremely high stability in denaturant, heat, and pH-induced unfolding transitions, which can be well described by a two-state reaction between native dimers and unfolded monomers. The homodimeric character of native ORF56 was confirmed by analytical ultracentrifugation. Far-UV circular dichroism and fluorescence spectroscopy gave superimposable denaturant-induced unfolding transitions and the midpoints of both heat as well as denaturant-induced unfolding depend on the protein concentration supporting the two-state model. This model was confirmed by GdmSCN-induced unfolding monitored by heteronuclear 2D NMR spectroscopy. Chemical denaturation was accomplished by GdmCl and GdmSCN, revealing a Gibbs free energy of stabilization of -85.1 kJ/mol at 25 degrees C. Thermal unfolding was possible only above 1 M GdmCl, which shifted the melting temperature (t(m)) below the boiling point of water. Linear extrapolation of t(m) to 0 M GdmCl yielded a t(m) of 107.5 degrees C (5 microM monomer concentration). Additionally, ORF56 remains natively structured over a remarkable pH range from pH 2 to pH 12. Folding kinetics were followed by far-UV CD and fluorescence after either stopped-flow or manual mixing. All kinetic traces showed only a single phase and the two probes revealed coincident folding rates (k(f), k(u)), indicating the absence of intermediates. Apparent first-order refolding rates depend linearly on the protein concentration, whereas the unfolding rates do not. Both lnk(f) and lnk(u) depend linearly on the GdmCl concentration. Together, folding and association of homodimeric ORF56 are concurrent events. In the absence of denaturant ORF56 refolds fast (7.0 x 10(7)M(-1)s(-1)) and unfolds extremely slowly (5.7 year(-1)). Therefore, high stability is coupled to a slow unfolding rate, which is often observed for proteins of extremophilic organisms.     

Structural and thermodynamic effects of ANS binding to human interleukin-1 receptor antagonist

Although 8-anilinonaphthalene-1-sulfonic acid (ANS) is frequently used in protein folding studies, the structural and thermodynamic effects of its binding to proteins are not well understood. Using high-resolution two-dimensional NMR and human interleukin-1 receptor antagonist (IL-1ra) as a model protein, we obtained detailed information on ANS-protein interactions in the absence and presence of urea. The effects of ambient to elevated temperatures on the affinity and specificity of ANS binding were assessed from experiments performed at 25 degrees C and 37 degrees C. Overall, the affinity of ANS was lower at 37 degrees C compared to 25 degrees C, but no significant change in the site specificity of binding was observed from the chemical shift perturbation data. The same site-specific binding was evident in the presence of 5.2 M urea, well within the unfolding transition region, and resulted in selective stabilization of the folded state. Based on the two-state denaturation mechanism, ANS-dependent changes in the protein stability were estimated from relative intensities of two amide resonances specific to the folded and unfolded states of IL-1ra. No evidence was found for any ANS-induced partially denatured or aggregated forms of IL-1ra throughout the experimental conditions, consistent with a cooperative and reversible denaturation process. The NMR results support earlier observations on the tendency of ANS to interact with solvent-exposed positively charged sites on proteins. Under denaturing conditions, ANS binding appears to be selective to structured states rather than unfolded conformations. Interestingly, the binding occurs within a previously identified aggregation-critical region in IL-1ra, thus providing an insight into ligand-dependent protein aggregation.     

Comparative analysis of thermal unfolding simulations of RNA recognition motifs (RRMs) of TAR DNA-binding protein 43 (TDP-43)

TAR DNA-binding protein 43 (TDP-43) inclusions have been found in Amyotrophic lateral sclerosis (ALS) and several other neurodegenerative diseases. Many studies suggest the involvement of RNA recognition motifs (RRMs) in TDP-43 proteinopathy. To elucidate the structural stability and the unfolding dynamics of RRMs, we have carried out atomistic molecular dynamics simulations at two different temperatures (300 and 500 K). The simulations results indicate that there are distinct structural differences in the unfolding pathway between the two domains and RRM1 unfolds faster than RRM2 in accordance with the lower thermal stability found experimentally. The unfolding behaviors of secondary structures showed that the α-helix was more stable than β-sheet and structural rearrangements of β-sheets results in formation of additional α-helices. At higher temperature, RRM1 exhibit increased overall flexibility and unfolding than RRM2. The temperature-dependent free energy landscapes consist of multiple metastable states stabilized by non-native contacts and hydrogen bonds in RRM2, thus rendering the RRM2 more prone to misfolding. The structural rearrangements of RRM2 could lead to aberrant protein–protein interactions that may account for enhanced aggregation and toxicity of TDP-43. Our analysis, thus identify the structural and thermodynamic characteristics of the RRMs of TDP-43, which will serve to uncover molecular mechanisms and driving forces in TDP-43 misfolding and aggregation.     

Comparative analyses of the conformational stability of a hyperthermophilic protein and its mesophilic counterpart.

Comparison of the conformational stability of an O(6)-methylguanine-DNA methyltransferase (MGMT) from the hyperthermophilic archaeon Thermococcus kodakaraensis strain KOD1 (Tk-MGMT), and its mesophilic counterpart C-terminal Ada protein from Escherichia coli (Ec-AdaC) was performed in order to obtain information about the relationship between thermal stability and other factors, such as thermodynamic parameters, thermodynamic stability and other unfolding conditions. Tk-MGMT unfolded at Tm = 98.6 degrees C, which was 54.8 degrees C higher than the unfolding temperature of Ec-AdaC. The maximum free energy (DeltaG(max)) of the proteins were different; the value of Tk-MGMT (42.9 kJ.mol-1 at 29.5 degrees C) was 2.6 times higher than that of Ec-AdaC (16.6 kJ.mol-1 at 7.4 degrees C). The high conformational stability of Tk-MGMT was attributed to a 1.6-fold higher enthalpy value than that of Ec-AdaC. In addition, the DeltaG(max) temperature of Tk-MGMT was considerably higher (by 22.1 degrees C). The apparent heat capacity of denaturation (DeltaC(p)) of Tk-MGMT was 0.7-fold lower than that of Ec-AdaC. These three synergistic effects, increasing DeltaGmax, shifted DeltaG vs. temperature curve, and low DeltaC(p), give Tk-MGMT its thermal stability. Unfolding profiles of the two proteins, tested with four alcohols and three denaturants, showed that Tk-MGMT possessed higher stability than Ec-AdaC in all conditions studied. These results indicate that the high stability of Tk-MGMT gives resistance to chemical unfolding, in addition to thermal unfolding.     

Enthalpic and entropic contributions mediate the role of disulfide bonds on the conformational stability of interleukin-4

The role of disulfide bridges in the structure, stability, and folding pathways of proteins has been the subject of wide interest in the fields of protein design and engineering. However, the relative importance of entropic and enthalpic contributions for the stabilization of proteins provided by disulfides is not always clear. Here, we perform a detailed analysis of the role of disulfides in the conformational stability of human Interleukin-4 (IL4), a four-helix bundle protein. In order to evaluate the contribution of two out of the three disulfides to the structure and stability of IL4, two IL4 mutants, C3T-IL4 and C24T-IL4, were used. NMR and ANS binding experiments were compatible with altered dynamics and an increase of the nonpolar solvent-accessible surface area of the folded state of the mutant proteins. Chemical and thermal unfolding experiments followed by fluorescence and circular dichroism revealed that both mutant proteins have lower conformational stability than the wild-type protein. Transition temperatures of unfolding decreased 14 degrees C for C3T-IL4 and 10 degrees C for C24T-IL4, when compared to WT-IL4, and the conformational stability, at 25 degrees C, decreased 4.9 kcal/mol for C3T-IL4 and 3.2 kcal/mol for C24T-IL4. Interestingly, both the enthalpy and the entropy of unfolding, at the transition temperature, decreased in the mutant proteins. Moreover, a smaller change in heat capacity of unfolding was also observed for the mutants. Thus, disulfide bridges in IL4 play a critical role in maintaining the thermodynamic stability and core packing of the helix bundle.     

Mannosylglycerate stabilizes staphylococcal nuclease with restriction of slow β-sheet motions

Mannosylglycerate is a compatible solute typical of thermophilic marine microorganisms that has a remarkable ability to protect proteins from thermal denaturation. This ionic solute appears to be a universal stabilizing agent, but the extent of protection depends on the specific protein examined. To understand how mannosylglycerate confers protection, we have been studying its influence on the internal motions of a hyperstable staphylococcal nuclease (SNase). Previously, we found a correlation between the magnitude of protein stabilization and the restriction of fast backbone motions. We now report the effect of mannosylglycerate on the fast motions of side-chains and on the slower unfolding motions of the protein. Side-chain motions were assessed by (13)CH(3) relaxation measurements and model-free analysis while slower unfolding motions were probed by H/D exchange measurements at increasing concentrations of urea. Side-chain motions were little affected by the presence of different concentrations of mannosylglycerate or even by the presence of urea (0.25M), and show no correlation with changes in the thermodynamic stability of SNase. Native hydrogen exchange experiments showed that, contrary to reports on other stabilizing solutes, mannosylglycerate restricts local motions in addition to the global motions of the protein. The protein unfolding/folding pathway remained undisturbed in the presence of mannosylglycerate but the solute showed a specific effect on the local motions of β-sheet residues. This work reinforces the link between solute-induced stabilization and restriction of protein motions at different timescales, and shows that the solute preferentially affects specific structural elements of SNase.     

Hydration state change of proteins upon unfolding in sugar solutions

Change in hydration number of proteins upon unfolding, Deltan, was obtained from the analysis of thermal unfolding behavior of proteins in various sugar solutions with water activity, a(W), varied. By applying the reciprocal form of Wyman-Tanford equation, Deltan was determined to be 133.9, 124.1, and 139.2 per protein molecule for ribonuclease A at pH=5.5, 4.2, and 2.8, respectively, 201.4 for lysozyme at pH=5.5, and 100.1 for alpha-chymotripnogen A at pH=2.0. Among the sugars tested, reducing sugars gave the lower apparent Deltan as compared with nonreducing sugars probably because of the direct interaction of reducing terminal with amino group of proteins at a high temperature. From the knowledge of Deltan, a new thermodynamic model for protein stability was proposed with explicit consideration for hydration state change of protein upon unfolding. From this model, the contribution of a(W) was proven to be always positive for stabilization of proteins and its effect is not negligible depending on Deltan and a(W).     

Thermodynamics and mechanism of cutinase stabilization by trehalose

Trehalose has been widely used to stabilize cellular structures such as membranes and proteins. The effect of trehalose on the stability of the enzyme cutinase was studied. Thermal unfolding of cutinase reveals that trehalose delays thermal unfolding, thus increasing the temperature at the midpoint of unfolding by 7.2 degrees . Despite this stabilizing effect, trehalose also favors pathways that lead to irreversible denaturation. Stopped-flow kinetics of cutinase folding and unfolding was measured and temperature was introduced as experimental variable to assess the mechanism and thermodynamics of protein stabilization by trehalose. The main stabilizing effect of trehalose was to delay the rate constant of the unfolding of an intermediate. A full thermodynamic analysis of this step has revealed that trehalose induces the phenomenon of entropy-enthalpy compensation, but the enthalpic contribution increases more significantly leading to a net stabilizing effect that slows down unfolding of the intermediate. Regarding the molecular mechanism of stabilization, trehalose increases the compactness of the unfolded state. The conformational space accessible to the unfolded state decreases in the presence of trehalose when the unfolded state acquires residual native interactions that channel the folding of the protein. This residual structure results into less hydrophobic groups being newly exposed upon unfolding, as less water molecules are immobilized upon unfolding.     

Non-linear effects of temperature and urea on the thermodynamics and kinetics of folding and unfolding of hisactophilin

Extensive measurements and analysis of thermodynamic stability and kinetics of urea-induced unfolding and folding of hisactophilin are reported for 5-50 degrees C, at pH 6.7. Under these conditions hisactophilin has moderate thermodynamic stability, and equilibrium and kinetic data are well fit by a two-state transition between the native and the denatured states. Equilibrium and kinetic m values decrease with increasing temperature, and decrease with increasing denaturant concentration. The betaF values at different temperatures and urea concentrations are quite constant, however, at about 0.7. This suggests that the transition state for hisactophilin unfolding is native-like and changes little with changing solution conditions, consistent with a narrow free energy profile for the transition state. The activation enthalpy and entropy of unfolding are unusually low for hisactophilin, as is also the case for the corresponding equilibrium parameters. Conventional Arrhenius and Eyring plots for both folding and unfolding are markedly non-linear, but these plots become linear for constant DeltaG/T contours. The Gibbs free energy changes for structural changes in hisactophilin have a non-linear denaturant dependence that is comparable to non-linearities observed for many other proteins. These non-linearities can be fit for many proteins using a variation of the Tanford model, incorporating empirical quadratic denaturant dependencies for Gibbs free energies of transfer of amino acid constituents from water to urea, and changes in fractional solvent accessible surface area of protein constituents based on the known protein structures. Noteworthy exceptions that are not well fit include amyloidogenic proteins and large proteins, which may form intermediates. The model is easily implemented and should be widely applicable to analysis of urea-induced structural transitions in proteins.     

Structure‐Based Site‐Specific PEGylation of Fibroblast Growth Factor 2 Facilitates Rational Selection of Conjugate Sites

Polyethylene glycol modification (PEGylation) can enhance the pharmacokinetic properties of therapeutic proteins by the attachment of polyethylene glycol (PEG) to the surface of a protein to shield the protein surface from proteolytic degradation and limit aggregation. However, current PEGylation strategies often reduce biological activity, potentially as a result of steric hindrance of PEG. Overall, there are no structure‐based guidelines for selection of conjugate sites that retain optimal biological activity with improved pharmacokinetic properties. In this study, site‐specific PEGylation based on the FGF2‐FGFR1‐heparin complex structure is performed. The effects of the conjugate sites on protein function are investigated by measuring the receptor/heparin binding affinities of the modified proteins and performing assays to measure cell‐based bio‐activity and in vivo stability. Comprehensive analysis of these data demonstrates that PEGylation of FGF2 that avoids the binding sites for fibroblast growth factor receptor 1 (FGFR1) and heparin provides optimal pharmacokinetic enhancement with minimal losses to biological activity. Animal experiments demonstrate that PEGylated FGF2 exhibits greater efficacy in protecting against traumatic brain injury‐induced brain damage and neurological functions than the non‐modified FGF2. This rational structure‐based PEGylation strategy for protein modification is expected to have a major impact in the area of protein‐based therapeutics.     

Thermodynamic effects of proline introduction on protein stability

The amino acid Pro is more rigid than other naturally occurring amino acids and, in proteins, lacks an amide hydrogen. To understand the structural and thermodynamic effects of Pro substitutions, it was introduced at 13 different positions in four different proteins, leucine-isoleucine-valine binding protein, maltose binding protein, ribose binding protein, and thioredoxin. Three of the maltose binding protein mutants were characterized by X-ray crystallography to confirm that no structural changes had occurred upon mutation. In the remaining cases, fluorescence and CD spectroscopy were used to show the absence of structural change. Stabilities of wild type and mutant proteins were characterized by chemical denaturation at neutral pH and by differential scanning calorimetry as a function of pH. The mutants did not show enhanced stability with respect to chemical denaturation at room temperature. However, 6 of the 13 single mutants showed a small but significant increase in the free energy of thermal unfolding in the range of 0.3-2.4 kcal/mol, 2 mutants showed no change, and 5 were destabilized. In five of the six cases, the stabilization was because of reduced entropy of unfolding. However, the magnitude of the reduction in entropy of unfolding was typically several fold larger than the theoretical estimate of -4 cal K(-1) mol(-1) derived from the relative areas in the Ramachandran map accessible to Pro and Ala residues, respectively. Two double mutants were constructed. In both cases, the effects of the single mutations on the free energy of thermal unfolding were nonadditive.     

Contribution of the multi-turn segment in the reversible thermal stability of hyperthermophile rubredoxin: NMR thermal chemical exchange analysis of sequence hybrids

Pyrococcus furiosus (Pf) rubredoxin is the most thermostable protein characterized to date. Reflecting the complications arising from irreversible denaturation of this protein, predictions of which structural regions confer differential thermal stability have utilized kinetic stability measurements, hydrogen exchange protection factors, long range hydrogen bond NMR spin couplings, and molecular dynamics simulations, and have primarily implicated the three-stranded beta-sheet and the adjacent metal binding site. Herein, NMR chemical exchange experiments demonstrate reversible two-state unfolding at the thermal transition temperature (T(m)) for hybrids of Pf and the mesophile Clostridium pasteurianum (Cp) rubredoxins which interchange residues 14-33, the so-called multi-turn segment. This complementary pair of hybrid rubredoxins exhibits largely additive incremental thermal stabilizations vs. the parental proteins. Both stabilization free energy measurements as well as incremental T(m) values indicate that a minimum of 37% of the total differential thermal stability resides in this multi-turn segment. Such a proportionality between DeltaDeltaG and incremental T(m) values is predicted for hybrid pairs exhibiting thermodynamic additivity in which the differential stability is predominantly enthalpic.