Comparison of (13)C(alpha)H and (15)NH backbone dynamics in protein GB1

This study presents a site-resolved experimental view of backbone C(alpha)H and NH internal motions in the 56-residue immunoglobulin-binding domain of streptococcal protein G, GB1. Using (13)C(alpha)H and (15)NH NMR relaxation data [T(1), T(2), and NOE] acquired at three resonance frequencies ((1)H frequencies of 500, 600, and 800 MHz), spectral density functions were calculated as F(omega) = 2omegaJ(omega) to provide a model-independent way to visualize and analyze internal motional correlation time distributions for backbone groups in GB1. Line broadening in F(omega) curves indicates the presence of nanosecond time scale internal motions (0.8 to 5 nsec) for all C(alpha)H and NH groups. Deconvolution of F(omega) curves effectively separates overall tumbling and internal motional correlation time distributions to yield more accurate order parameters than determined by using standard model free approaches. Compared to NH groups, C(alpha)H internal motions are more broadly distributed on the nanosecond time scale, and larger C(alpha)H order parameters are related to correlated bond rotations for C(alpha)H fluctuations. Motional parameters for NH groups are more structurally correlated, with NH order parameters, for example, being larger for residues in more structured regions of beta-sheet and helix and generally smaller for residues in the loop and turns. This is most likely related to the observation that NH order parameters are correlated to hydrogen bonding. This study contributes to the general understanding of protein dynamics and exemplifies an alternative and easier way to analyze NMR relaxation data.     

Backbone dynamics of a winged helix protein and its DNA complex at different temperatures: changes of internal motions in genesis upon binding to DNA.

The dynamic properties of a winged helix protein, Genesis, and its DNA complex at different temperatures were studied. Due to the complexity of motions, the commonly used model-free formalism could not be used to reflect the dynamic properties. The reduced spectral density function mapping approach was proven to be a useful tool to describe the overall and internal motion of molecules on the picosecond to nanosecond time-scale, and conformational exchanges on the microsecond to millisecond time-scale. The local motions in DNA-free Genesis showed strong temperature dependence and the backbone dynamics of each secondary structural element responds to the temperature change differently, while the Genesis-DNA complex showed more stability with changing the temperatures. Furthermore, each DNA contact sequence of Genesis showed distinct dynamic perturbation after Genesis binds to DNA.     

A conserved π-helix plays a key role in thermoadaptation of catalysis in the glycoside hydrolase family 4

Here, we characterize the role of a π-helix in the molecular mechanisms underlying thermoadaptation in the glycoside hydrolase family 4 (GH4). The interspersed π-helix present in a subgroup is evolutionarily related to a conserved α-helix in other orthologs by a single residue insertion/deletion event. The insertional residue, Phe407, in a hyperthermophilic α-glucuronidase, makes specific interactions across the inter-subunit interface. In order to establish the sequence-structure-stability implications of the π-helix, the wild-type and the deletion variant (Δ407) were characterized. The variant showed a significant lowering of melting temperature and optimum temperature for the highest activity. Crystal structures of the proteins show a transformation of the π-helix to a continuous α-helix in the variant, identical to that in orthologs lacking this insertion. Thermodynamic parameters were determined from stability curves representing the temperature dependence of unfolding free energy. Though the proteins display maximum stabilities at similar temperatures, a higher melting temperature in the wild-type is achieved by a combination of higher enthalpy and lower heat capacity of unfolding. Comparisons of the structural changes, and the activity and thermodynamic profiles allow us to infer that specific non-covalent interactions, and the existence of residual structure in the unfolded state, are crucial determinants of its thermostability. These features permit the enzyme to balance the preservation of structure at a higher temperature with the thermodynamic stability required for optimum catalysis.     

Dynamical characterization of residual and non-native structures in a partially folded protein by (15)N NMR relaxation using a model based on a distribution of correlation times

A spectral density model based on a truncated lorentzian distribution of correlation times is used to analyze the nanosecond time-scale dynamics of the partially unfolded domain 2 of annexin I from its (15)N NMR relaxation parameters measured at three magnetic field strengths. The use of a distribution of correlation times enables the characterization of the dynamical features of the NH bonds of the protein in terms of heterogeneity of dynamical states in the nanosecond range. The variation along the sequence of the two dynamical parameters introduced, namely the center and the width of the distribution, points out the different types of residual secondary structures present in the D2 domain. Moreover, it allows a physically sensible interpretation of the dynamical behavior of the different residual helices and of the non-native structures. Also, a striking correspondence is found between the parameters obtained using an extended Lipari and Szabo model and the parameters obtained using the distribution of correlation times. This result led us to propose a specific interpretation of the model-free order parameter for internal motions in the nanosecond range in the case of unfolded states.     

13C NMR relaxation studies of RNA base and ribose nuclei reveal a complex pattern of motions in the RNA binding site for human U1A protein

The widespread importance of induced fit and order-disorder transition in RNA recognition by proteins and small molecules makes it imperative that RNA motional properties are characterized quantitatively. Until now, however, very few studies have been dedicated to the systematic characterization of RNA motion and to their changes upon protein or small-molecule binding. The U1A protein-RNA complexes provide some of the best-studied examples of the role of RNA motional changes upon protein binding. Here, we report (13)C NMR relaxation studies of base and ribose dynamics for the RNA internal loop target of human U1A protein located within the 3'-untranslated region (3'-UTR) of the mRNA coding for U1A itself. We also report the semi-quantitative analysis of both fast (nano- to picosecond) and intermediate (micro- to millisecond) motions for this paradigmatic RNA system. We measure (13)C T(1), T(1rho) and heteronuclear nuclear Overhauser effects (NOEs) for sugar and base nuclei, as well as the power dependence of T(1rho) at 500 MHz and 750 MHz, and analyze these results using the model-free formalism. The results provide a much clearer picture of the type of motions experienced by this RNA in the absence of the protein than was provided by the analysis of the structure based solely on NOEs and scalar couplings. They define a model where the RNA internal loop region "breathes" on a micro- to millisecond timescale with respect to the double-helical regions. Superimposed on this slower motion, the residues at the very tip of the loop undergo faster (nano- to picosecond) motions. We hypothesize that these motions allow the RNA to sample multiple conformations so that the protein can select a structure within the ensemble that optimizes intermolecular contacts.     

A simple model of backbone flexibility improves modeling of side-chain conformational variability

The considerable flexibility of side-chains in folded proteins is important for protein stability and function, and may have a role in mediating allosteric interactions. While sampling side-chain degrees of freedom has been an integral part of several successful computational protein design methods, the predictions of these approaches have not been directly compared to experimental measurements of side-chain motional amplitudes. In addition, protein design methods frequently keep the backbone fixed, an approximation that may substantially limit the ability to accurately model side-chain flexibility. Here, we describe a Monte Carlo approach to modeling side-chain conformational variability and validate our method against a large dataset of methyl relaxation order parameters derived from nuclear magnetic resonance (NMR) experiments (17 proteins and a total of 530 data points). We also evaluate a model of backbone flexibility based on Backrub motions, a type of conformational change frequently observed in ultra-high-resolution X-ray structures that accounts for correlated side-chain backbone movements. The fixed-backbone model performs reasonably well with an overall rmsd between computed and predicted side-chain order parameters of 0.26. Notably, including backbone flexibility leads to significant improvements in modeling side-chain order parameters for ten of the 17 proteins in the set. Greater accuracy of the flexible backbone model results from both increases and decreases in side-chain flexibility relative to the fixed-backbone model. This simple flexible-backbone model should be useful for a variety of protein design applications, including improved modeling of protein-protein interactions, design of proteins with desired flexibility or rigidity, and prediction of correlated motions within proteins.     

Factors affecting the thermodynamic stability of small asymmetric internal loops in RNA

Optical melting experiments were used to determine the thermodynamic parameters for oligoribonucleotides containing small asymmetric internal loops. The results show a broad range of thermodynamic stabilities, which depend on loop size, asymmetry, sequence, closing base pairs, and length of helix stems. Imino proton NMR experiments provide evidence for possible hydrogen bonding in GA and UU mismatches in some asymmetric loops. The stabilizing effects of GA, GG, and UU mismatches on the thermodynamic stability of internal loops vary depending on the size and asymmetry of the loop. The dependence of loop stability on Watson-Crick closing base pairs may be explained by an account of hydrogen bonds. Models are presented for approximating the free energy increments of 2 x 3 and 1 x 3 internal loops.     

Mapping of the spectral densities of N-H bond motions in eglin c using heteronuclear relaxation experiments.

A new strategy is used for studying the internal motions of proteins based on measurements of NMR relaxation parameters. The strategy yields values of the so-called spectral density functions J(omega) for N-H bond vectors. The spectral density functions are related to the distribution of frequencies contained in the rotational (overall and internal) motions of these NH bond vectors. No a priori model assumptions about the dynamics are required in this approach. The method involves measurements of six relaxation parameters consisting of 15N longitudinal relaxation rates, transverse relaxation rates of in-phase and antiphase coherence, the relaxation rates of heteronuclear 1H-15N two-spin order, the heteronuclear 1H-15N nuclear Overhauser effects, and longitudinal relaxation rates of the amide protons. The values of the spectral density functions at the five frequencies 0, omega N, omega H + omega N, omega H, and omega H - omega N are determined from the relaxation parameters using analytical relations derived previously [Peng & Wagner (1992) J. Magn. Reson. 98, 308-332]. Here, the method is applied to characterize the backbone dynamics of the 15N-enriched proteinase inhibitor eglin c, a protein of 70 residues. The values for J(0) and J(omega N = 50 MHz) vary significantly with the amino acid sequence, whereas the spectral densities at higher frequencies, J(450 MHz), J(500 MHz), and J(550 MHz), are typically much smaller and show no significant variation with the sequence. The collective behavior of the J(omega) values indicate greater internal motion for the proteinase binding loop residues and the first eight N-terminal residues. The additional internal motion in these regions is in the rate range below 450 MHz. The values of J(omega) are also compared with root mean square deviations (rmsds) of backbone atoms as obtained in NMR structure determinations. Low values of J(0) and J(omega N) are correlated with high rmsds. Spectral densities at higher frequencies, J(450 MHz), J(500 MHz), and J(550 MHz), are small and show no correlation with rmsds. A comparison with the spectral density functions obtained by fitting the experimental data to the functional dependence of the Lipari and Szabo formalism [Lipari & Szabo (1982a) J. Am. Chem. Soc. 104, 4546-4559] is made.     

Hydrogen exchange and proteolytic degradation of ribonuclease A. Similarities and distinctions of the kinetic mechanisms

The studies by IR spectroscopy of the temperature dependence of the H-D exchange rate of the RNase A peptide NH atoms permit one to characterize two types of conformation fluctuations, local and global. A comparison with the temperature dependence of the proteolytic degradation rate of RNase A shows that similar in nature fluctuations allow for the H-D exchange of NH atoms and the splitting of peptide bonds of the native protein. In the low temperature region, both processes occur through local fluctuations, by way of the EX2 mechanism, and in the high temperature region, they occur through global fluctuations with the overall denaturation desorganization of the native structure, by way of the EX1 mechanism. The biphasic dependence of the rate of H-D exchange and proteolytic degradation of RNase A on urea concentration is also explained by the combination of local and global fluctuations.     

Structural heterogeneity of type I collagen triple helix and its role in osteogenesis imperfecta

We investigated regions of different helical stability within human type I collagen and discussed their role in intermolecular interactions and osteogenesis imperfecta (OI). By differential scanning calorimetry and circular dichroism, we measured and mapped changes in the collagen melting temperature (DeltaTm) for 41 different Gly substitutions from 47 OI patients. In contrast to peptides, we found no correlations of DeltaTm with the identity of the substituting residue. Instead, we observed regular variations in DeltaTm with the substitution location in different triple helix regions. To relate the DeltaTm map to peptide-based stability predictions, we extracted the activation energy of local helix unfolding (DeltaG) from the reported peptide data. We constructed the DeltaG map and tested it by measuring the H-D exchange rate for glycine NH residues involved in interchain hydrogen bonds. Based on the DeltaTm and DeltaG maps, we delineated regional variations in the collagen triple helix stability. Two large, flexible regions deduced from the DeltaTm map aligned with the regions important for collagen fibril assembly and ligand binding. One of these regions also aligned with a lethal region for Gly substitutions in the alpha1(I) chain.     

A Link between Hinge-Bending Domain Motions and the Temperature Dependence of Catalysis in 3-Isopropylmalate Dehydrogenase

Enzyme function depends on specific conformational motions. We show that the temperature dependence of enzyme kinetic parameters can provide insight into these functionally relevant motions. While investigating the catalytic properties of IPMDH from Escherichia coli, we found that its catalytic efficiency (k_cat/K_M,IPM) for the substrate IPM has an unusual temperature dependence, showing a local minimum at ∼35°C. In search of an explanation, we measured the individual constants k_cat and K_M,IPM as a function of temperature, and found that the van 't Hoff plot of K_M,IPM shows sigmoid-like transition in the 20-40°C temperature range. By means of various measurements including hydrogen-deuterium exchange and fluorescence resonance energy transfer, we showed that the conformational fluctuations, including hinge-bending domain motions increase more steeply with temperatures >30°C. The thermodynamic parameters of ligand binding determined by isothermal titration calorimetry as a function of temperature were found to be strongly correlated to the conformational fluctuations of the enzyme. Because the binding of IPM is associated with a hinge-bending domain closure, the more intense hinge-bending fluctuations at higher temperatures increasingly interfere with IPM binding, thereby abruptly increasing its dissociation constant and leading to the observed unusual temperature dependence of the catalytic efficiency.     

Internal motional amplitudes and correlated bond rotations in an alpha-helical peptide derived from 13C and 15N NMR relaxation

Peptide GFSKAELAKARAAKRGGY folds in an alpha-helical conformation that is stabilized by formation of a hydrophobic staple motif and an N-terminal capping box (Munoz V. Blanco FJ, Serrano L, 1995, Struct Biol 2:380-385). To investigate backbone and side-chain internal motions within the helix and hydrophobic staple, residues F2, A5, L7, A8, and A10 were selectively 13C- and 15N-enriched and NMR relaxation experiments were performed in water and in water/trifluoroethanol (TFE) solution at four Larmor frequencies (62.5, 125, 150, and 200 MHz for 13C). Relaxation data were analyzed using the model free approach and an anisotropic diffusion model. In water, angular variances of motional vectors range from 10 to 20 degrees and backbone phi,psi bond rotations for helix residues A5, L7, A8, and A10 are correlated indicating the presence of Calpha-H, Calpha-Cbeta, and N-H rocking-type motions along the helix dipole axis. L7 side-chain CbetaH2 and CgammaH motions are also correlated and as motionally restricted as backbone CalphaH, suggesting considerable steric hindrance with neighboring groups. In TFE which stabilizes the fold, internal motional amplitudes are attenuated and rotational correlations are increased. For the side chain of hydrophobic staple residue F2, wobbling-in-a-cone type motions dominate in water, whereas in TFE, the Cbeta-Cgamma bond and phenyl ring fluctuate more simply about the Calpha-Cbeta bond. These data support the Daragan-Mayo model of correlated bond rotations (Daragan VA, Mayo KH, 1996, J Phys Chem 100:8378-8388) and contribute to a general understanding of internal motions in peptides and proteins.     

Comparing the thermodynamic stabilities of a related thermophilic and mesophilic enzyme

Several models have been proposed to explain the high temperatures required to denature enzymes from thermophilic organisms; some involve greater maximum thermodynamic stability for the thermophile, and others do not. To test these models, we reversibly melted two analogous protein domains in a two-state manner. E2cd is the isolated catalytic domain of cellulase E2 from the thermophile Thermomonospora fusca. CenAP30 is the analogous domain of the cellulase CenA from the mesophile Cellulomonas fimi. When reversibly denatured in a common buffer, the thermophilic enzyme E2cd had a temperature of melting (Tm) of 72.2 degrees C, a van't Hoff enthalpy of unfolding (DeltaHVH) of 190 kcal/mol, and an entropy of unfolding (DeltaSu) of 0.55 kcal/(mol*K); the mesophilic enzyme CenAP30 had a Tm of 56.4 degrees C, a DeltaHVH of 107 kcal/mol, and a DeltaSu of 0. 32 kcal/(mol*K). The higher DeltaHVH and DeltaSu values for E2cd suggest that its free energy of unfolding (DeltaGu) has a steeper dependence on temperature at the Tm than CenAP30. This result supports models that predict a greater maximum thermodynamic stability for thermophilic enzymes than for their mesophilic counterparts. This was further explored by urea denaturation. Under reducing conditions at 30 degrees C, E2cd had a concentration of melting (Cm) of 5.2 M and a DeltaGu of 11.2 kcal/mol; CenAP30 had a Cm of 2.6 M and a DeltaGu of 4.3 kcal/mol. Under nonreducing conditions, the Cm and DeltaGu of CenAP30 were increased to 4.5 M and 10.8 kcal/mol at 30 degrees C; the Cm for E2cd was increased to at least 7.4 M at 32 degrees C. We were unable to determine a DeltaGu value for E2cd under nonreducing conditions due to problems with reversibility. These data suggest that E2cd attains its greater thermal stability (DeltaTm = 15.8 degrees C) through a greater thermodynamic stability (DeltaDeltaGu = 6.9 kcal/mol) compared to its mesophilic analogue CenAP30.     

The thermodynamic stability of the proteins of the ccd plasmid addiction system

The two opponents, toxin (CcdB, LetB or LetD, protein G, LynB) and antidote (CcdA, LetA, protein H, LynA), in the plasmid addiction system ccd of the F plasmid were studied by different biophysical methods. The thermodynamic stability was measured at different temperatures combining denaturant and thermally induced unfolding. It was found that both proteins denature in a two-state equilibrium (native dimer versus unfolded monomer) and that CcdA has a significantly lower thermodynamic stability. Using a numerical model, which was developed earlier by us, and on the basis of the determined thermodynamic parameters the concentration dependence of the denaturation transition temperature was obtained for both proteins. This concentration dependence may be of physiological significance, as the concentration of both ccd addiction proteins cannot exceed a certain limit because their expression is controlled by autoregulation.The influence of DNA on the thermal stability of the two proteins was probed. It was found that cognate DNA increases the melting temperature of CcdA. In the presence of non-specific DNA the thermal stability was not changed. The melting temperature of CcdB was not influenced by the applied double-stranded oligonucleotides, neither cognate nor unspecific.     

Theory and applications of the generalized Born solvation model in macromolecular simulations.

Generalized Born (GB) models provide an attractive way to include some thermodynamic aspects of aqueous solvation into simulations that do not explicitly model the solvent molecules. Here we discuss our recent experience with this model, presenting in detail the way it is implemented and parallelized in the AMBER molecular modeling code. We compare results using the GB model (or GB plus a surface-area based "hydrophobic" term) to explicit solvent simulations for a 10 base-pair DNA oligomer, and for the 108-residue protein thioredoxin. A slight modification of our earlier suggested parameters makes the GB results more like those found in explicit solvent, primarily by slightly increasing the strength of NH [bond] O and NH [bond] N internal hydrogen bonds. Timing and energy stability results are reported, with an eye toward using these model for simulations of larger macromolecular systems and longer time scales.     

Backbone dynamics of barstar: a (15)N NMR relaxation study

Backbone dynamics of uniformly (15)N-labeled barstar have been studied at 32 degrees C, pH 6.7, by using (15)N relaxation data obtained from proton-detected 2D (1)H-(15)N NMR spectroscopy. (15)N spin-lattice relaxation rate constants (R(1)), spin-spin relaxation rate constants (R(2)), and steady-state heteronuclear (1)H-(15)N NOEs have been determined for 69 of the 86 (excluding two prolines and the N-terminal residue) backbone amide (15)N at a magnetic field strength of 14.1 Tesla. The primary relaxation data have been analyzed by using the model-free formalism of molecular dynamics, using both isotropic and axially symmetric diffusion of the molecule, to determine the overall rotational correlation time (tau(m)), the generalized order parameter (S(2)), the effective correlation time for internal motions (tau(e)), and NH exchange broadening contributions (R(ex)) for each residue. As per the axially symmetric diffusion, the ratio of diffusion rates about the unique and perpendicular axes (D( parallel)/D( perpendicular)) is 0.82 +/- 0.03. The two results have only marginal differences. The relaxation data have also been used to map reduced spectral densities for the NH vectors of these residues at three frequencies: 0, omega(H), and omega(N), where omega(H),(N) are proton and nitrogen Larmor frequencies. The value of tau(m) obtained from model-free analysis of the relaxation data is 5.2 ns. The reduced spectral density analysis, however, yields a value of 5.7 ns. The tau(m) determined here is different from that calculated previously from time-resolved fluorescence data (4.1 ns). The order parameter ranges from 0.68 to 0.98, with an average value of 0.85 +/- 0.02. A comparison of the order parameters with the X-ray B-factors for the backbone nitrogens of wild-type barstar does not show any considerable correlation. Model-free analysis of the relaxation data for seven residues required the inclusion of an exchange broadening term, the magnitude of which ranges from 2 to 9.1 s(-1), indicating the presence of conformational averaging motions only for a small subset of residues.     

Thermodynamics of the coil <==> beta-sheet transition in a membrane environment

Biologically important peptides such as the Alzheimer peptide Abeta(1-40) display a reversible random coil <==>beta-structure transition at anionic membrane surfaces. In contrast to the well-studied random coil left arrow over right arrow alpha-helix transition of amphipathic peptides, there is a dearth on information on the thermodynamic and kinetic parameters of the random coil left arrow over right arrow beta-structure transition. Here, we present a new method to quantitatively analyze the thermodynamic parameters of the membrane-induced beta-structure formation. We have used the model peptide (KIGAKI)(3) and eight analogues in which two adjacent amino acids were substituted by their d-enantiomers. The positions of the d,d pairs were shifted systematically along the three identical segments of the peptide chain. The beta-structure content of the peptides was measured in solution and when bound to anionic lipid membranes with circular dichroism spectroscopy. The thermodynamic binding parameters were determined with isothermal titration calorimetry and the binding isotherms were analysed by combining a surface partition equilibrium with the Gouy-Chapman theory. The thermodynamic parameters were found to be linearly correlated with the extent of beta-structure formation. beta-Structure formation at the membrane surface is characterized by an enthalpy change of DeltaH(beta)=-0.23 kcal/mol per residue, an entropy change of DeltaS(beta)=-0.24 cal/mol K residue and a free energy change of DeltaG(beta)=-0.15 kcal/mol residue. An increase in temperature induces an unfolding of beta-structure. The residual free energy of membrane-induced beta-structure formation is close to that of membrane-induced alpha-helix formation.     

Thermodynamic and kinetic exploration of the energy landscape of Borrelia burgdorferi OspA by native-state hydrogen exchange

We report a native-state hydrogen-exchange (HX) method to simultaneously obtain both thermodynamic and kinetic information on the formation of multiple excited states in a folding energy landscape. Our method exploits the inherent dispersion and pH dependence of the intrinsic HX rates to cover both the EX2 (thermodynamic) and EX1 (kinetic) regimes. At each concentration of denaturant, HX measurements are performed over a range of pH values. Using this strategy, we dissected Borrelia burgdorferi OspA, a predominantly beta-sheet protein containing a unique single-layer beta-sheet, into five cooperative units and postulated excited states predominantly responsible for HX. More importantly, we determined the interconversion rates between these excited states and the native state. The use of both thermodynamic and kinetic information from native-state HX enabled us to construct a folding landscape of this 28kDa protein, including local minima and maxima, and to discriminate on-pathway and off-pathway intermediates. This method, which we term EX2/EX1 HX, should be a powerful tool for characterizing the complex folding mechanisms exhibited by the majority of proteins.     

High thermal stability of 3-isopropylmalate dehydrogenase from Thermus thermophilus resulting from low DeltaC(p) of unfolding.

To characterize the thermal stability of 3-isopropylmalate dehydrogenase (IPMDH) from an extreme thermophile, Thermus thermophilus, urea-induced unfolding of the enzyme and of its mesophilic counterpart from Escherichia coli was investigated at various temperatures. The unfolding curves were analyzed with a three-state model for E.coli IPMDH and with a two-state model for T.thermophilus IPMDH, to obtain the free energy change DeltaG degrees of each unfolding process. Other thermodynamic parameters, enthalpy change DeltaH, entropy change DeltaS and heat capacity change DeltaC(p), were derived from the temperature dependence of DeltaG degrees. The main feature of the thermophilic enzyme was its lower dependence of DeltaG degrees on temperature resulting from a low DeltaC(p). The thermophilic IPMDH had a significantly lower DeltaC(p), 1.73 kcal/mol.K, than that of E.coli IPMDH (20.7 kcal/mol.K). The low DeltaC(p) of T.thermophilus IPMDH could not be predicted from its change in solvent-accessible surface area DeltaASA. The results suggested that there is a large structural difference between the unfolded state of T.thermophilus and that of E.coli IPMDH. Another responsible factor for the higher thermal stability of T.thermophilus IPMDH was the increase in the most stable temperature T(s). The DeltaG degrees maximum of T.thermophilus IPMDH was much smaller than that of E.coli IPMDH. The present results clearly demonstrated that a higher melting temperature T(m) is not necessarily accompanied by a higher DeltaG degrees maximum.     

Activity-stability relationships in extremophilic enzymes

Psychrophilic, mesophilic, and thermophilic alpha-amylases have been studied as regards their conformational stability, heat inactivation, irreversible unfolding, activation parameters of the reaction, properties of the enzyme in complex with a transition state analog, and structural permeability. These data allowed us to propose an energy landscape for a family of extremophilic enzymes based on the folding funnel model, integrating the main differences in conformational energy, cooperativity of protein unfolding, and temperature dependence of the activity. In particular, the shape of the funnel bottom, which depicts the stability of the native state ensemble, also accounts for the thermodynamic parameters of activation that characterize these extremophilic enzymes, therefore providing a rational basis for stability-activity relationships in protein adaptation to extreme temperatures.