Dissecting subunit interfaces in homodimeric proteins

The subunit interfaces of 122 homodimers of known three-dimensional structure are analyzed and dissected into sets of surface patches by clustering atoms at the interface; 70 interfaces are single-patch, the others have up to six patches, often contributed by different structural domains. The average interface buries 1,940 A2 of the surface of each monomer, contains one or two patches burying 600-1,600 A2, is 65% nonpolar and includes 18 hydrogen bonds. However, the range of size and of hydrophobicity is wide among the 122 interfaces. Each interface has a core made of residues with atoms buried in the dimer, surrounded by a rim of residues with atoms that remain accessible to solvent. The core, which constitutes 77% of the interface on average, has an amino acid composition that resembles the protein interior except for the presence of arginine residues, whereas the rim is more like the protein surface. These properties of the interfaces in homodimers, which are permanent assemblies, are compared to those of protein-protein complexes where the components associate after they have independently folded. On average, subunit interfaces in homodimers are twice larger than in complexes, and much less polar due to the large fraction belonging to the core, although the amino acid compositions of the cores are similar in the two types of interfaces.     

The Subunit Interfaces of Weakly Associated Homodimeric Proteins

We analyzed subunit interfaces in 315 homodimers with an X-ray structure in the Protein Data Bank, validated by checking the literature for data that indicate that the proteins are dimeric in solution and that, in the case of the “weak” dimers, the homodimer is in equilibrium with the monomer. The interfaces of the 42 weak dimers, which are smaller by a factor of 2.4 on average than in the remainder of the set, are comparable in size with antibody-antigen or protease-inhibitor interfaces. Nevertheless, they are more hydrophobic than in the average transient protein-protein complex and similar in amino acid composition to the other homodimer interfaces. The mean numbers of interface hydrogen bonds and hydration water molecules per unit area are also similar in homodimers and transient complexes. Parameters related to the atomic packing suggest that many of the weak dimer interfaces are loosely packed, and we suggest that this contributes to their low stability. To evaluate the evolutionary selection pressure on interface residues, we calculated the Shannon entropy of homologous amino acid sequences at 60% sequence identity. In 93% of the homodimers, the interface residues are better conserved than the residues on the protein surface. The weak dimers display the same high degree of interface conservation as other homodimers, but their homologs may be heterodimers as well as homodimers. Their interfaces may be good models in terms of their size, composition, and evolutionary conservation for the labile subunit contacts that allow protein assemblies to share and exchange components, allosteric proteins to undergo quaternary structure transitions, and molecular machines to operate in the cell.     

Distinct roles of overlapping and non-overlapping regions of hub protein interfaces in recognition of multiple partners

Cellular functions of an organism are maintained by protein-protein interactions. Those proteins that bind multiple partners asynchronously (date hub proteins) are important to make the interaction network coordinated. It is known that many date hub proteins bind different partners at overlapping (OV) interfaces. To understand how OV interfaces of date hub proteins can recognize multiple partners, we analyzed the difference between OV and non-overlapping (Non-OV) regions of interfaces involved in the binding of different partners. By using the structures of 16 date hub proteins with various interaction partners (ranging from 5 to 33), we compared buried surface area, compositions of amino acid residues and secondary structures, and side-chain orientations. It was found that buried interface residues are important for recognizing multiple partners, while exposed interface residues are important for determining specificity to a particular ligand. In addition, our analyses reveal that residue compositions in OV and Non-OV regions are different and that residues in OV region show diverse side-chain torsion angles to accommodate binding to multiple targets.     

A dissection of the protein-protein interfaces in icosahedral virus capsids

We selected 49 icosahedral virus capsids whose crystal structures are reported in the Protein Data Bank. They belong to the T=1, T=3, pseudo T=3 and other lattice types. We identified in them 779 unique interfaces between pairs of subunits, all repeated by icosahedral symmetry. We analyzed the geometric and physical chemical properties of these interfaces and compared with interfaces in protein-protein complexes and homodimeric proteins, and with crystal packing contacts. The capsids contain one to 16 subunits implicated in three to 66 unique interfaces. Each subunit loses 40-60% of its accessible surface in contacts with an average of 8.5 neighbors. Many of the interfaces are very large with a buried surface area (BSA) that can exceed 10,000 A(2), yet 39% are small with a BSA<800 A(2) comparable to crystal packing contacts. Pairwise capsid interfaces overlap, so that one-third of the residues are part of more than one interface. Those with a BSA>800 A(2) resemble homodimer interfaces in their chemical composition. Relative to the protein surface, they are non-polar, enriched in aliphatic residues and depleted of charged residues, but not of neutral polar residues. They contain one H-bond per about 200 A(2) BSA. Small capsid interfaces (BSA<800 A(2)) are only slightly more polar. They have a similar amino acid composition, but they bury fewer atoms and contain fewer H-bonds for their size. Geometric parameters that estimate the quality of the atomic packing suggest that the small capsid interfaces are loosely packed like crystal packing contacts, whereas the larger interfaces are close-packed as in protein-protein complexes and homodimers. We discuss implications of these findings on the mechanism of capsid assembly, assuming that the larger interfaces form first to yield stable oligomeric species (capsomeres), and that medium-size interfaces allow the stepwise addition of capsomeres to build larger intermediates.     

Dissecting protein-protein recognition sites

The recognition sites in 70 pairwise protein-protein complexes of known three-dimensional structure are dissected in a set of surface patches by clustering atoms at the interface. When the interface buries <2000 A2 of protein surface, the recognition sites usually form a single patch on the surface of each component protein. In contrast, larger interfaces are generally multipatch, with at least one pair of patches that are equivalent in size to a single-patch interface. Each recognition site, or patch within a site, contains a core made of buried interface atoms, surrounded by a rim of atoms that remain accessible to solvent in the complex. A simple geometric model reproduces the number and distribution of atoms within a patch. The rim is similar in composition to the rest of the protein surface, but the core has a distinctive amino acid composition, which may help in identifying potential protein recognition sites on single proteins of known structures.     

Dissection, residue conservation, and structural classification of protein-DNA interfaces

The basic DNA-binding modules of 128 protein-DNA interfaces have been analyzed. Although these are less planar, like the protein-protein interfaces, the protein-DNA interfaces can also be dissected into core regions in which all the fully-buried atoms are located, and rim regions having atoms with residual accessibilities. The sequence entropy of the core residues is smaller than those in the rim, indicating that the former are better conserved and possibly contribute more towards the binding free energy, as has been implicated in protein-protein interactions. On the protein side, 1014 A(2) of the surface is buried of which 63% belong to the core. There are some differences in the propensities of residues to occur in the core and the rim. In the DNA strands, the nucleotide(s) containing fully-buried atoms in all three components usually occupy central positions of the binding region. A new classification scheme for the interfaces has been introduced based on the composition of secondary structural elements of residues and the results compared with the conventional classification of DNA-binding proteins, as well as the protein class of the molecule. It appears that a common framework may be developed to understand both protein-protein and protein-DNA interactions.     

Conformational equilibria and rates of localized motion within hepatitis B virus capsids

Functional analysis of hepatitis B virus (HBV) core particles has associated a number of biological roles with the C terminus of the capsid protein. One set of functions require the C terminus to be on the exterior of the capsid, while others place this domain on the interior. According to the crystal structure of the capsid, this segment is strictly internal to the capsid shell and buried at a protein-protein interface. Using kinetic hydrolysis, a form of protease digestion assayed by SDS-PAGE and mass spectrometry, the structurally and biologically important C-terminal region of HBV capsid protein assembly domain (Cp149, residues 1-149) has been shown to be dynamic in both dimer and capsid forms. HBV is an enveloped virus with a T = 4 icosahedral core that is composed of 120 copies of a homodimer capsid protein. Free dimer and assembled capsid forms of the protein are readily hydrolyzed by trypsin and thermolysin, around residues 127-128, indicating that this region is dynamic and exposed to the capsid surface. The measured conformational equilibria have an opposite temperature dependence between free dimer and assembled capsid. This work helps to explain the previously described allosteric regulation of assembly and functional properties of a buried domain. These observations make a critical connection between structure, dynamics, and function: made possible by the first quantitative measurements of conformational equilibria and rates of conversion between protein conformers for a megaDalton complex.     

YopD Self-assembly and Binding to LcrV Facilitate Type III Secretion Activity by Yersinia pseudotuberculosis

YopD-like translocator proteins encoded by several Gram-negative bacteria are important for type III secretion-dependent delivery of anti-host effectors into eukaryotic cells. This probably depends on their ability to form pores in the infected cell plasma membrane, through which effectors may gain access to the cell interior. In addition, Yersinia YopD is a negative regulator essential for the control of effector synthesis and secretion. As a prerequisite for this functional duality, YopD may need to establish molecular interactions with other key T3S components. A putative coiled-coil domain and an α-helical amphipathic domain, both situated in the YopD C terminus, may represent key protein-protein interaction domains. Therefore, residues within the YopD C terminus were systematically mutagenized. All 68 mutant bacteria were first screened in a variety of assays designed to identify individual residues essential for YopD function, possibly by providing the interaction interface for the docking of other T3S proteins. Mirroring the effect of a full-length yopD gene deletion, five mutant bacteria were defective for both yop regulatory control and effector delivery. Interestingly, all mutations clustered to hydrophobic amino acids of the amphipathic domain. Also situated within this domain, two additional mutants rendered YopD primarily defective in the control of Yop synthesis and secretion. Significantly, protein-protein interaction studies revealed that functionally compromised YopD variants were also defective in self-oligomerization and in the ability to engage another translocator protein, LcrV. Thus, the YopD amphipathic domain facilitates the formation of YopD/YopD and YopD/LcrV interactions, two critical events in the type III secretion process.     

Role of β-lactamase residues in a common interface for binding the structurally unrelated inhibitory proteins BLIP and BLIP-II

The β-lactamase inhibitory proteins (BLIPs) are a model system for examining molecular recognition in protein-protein interactions. BLIP and BLIP-II are structurally unrelated proteins that bind and inhibit TEM-1 β-lactamase. Both BLIPs share a common binding interface on TEM-1 and make contacts with many of the same TEM-1 surface residues. BLIP-II, however, binds TEM-1 over 150-fold tighter than BLIP despite the fact that it has fewer contact residues and a smaller binding interface. The role of eleven TEM-1 amino acid residues that contact both BLIP and BLIP-II was examined by alanine mutagenesis and determination of the association (k_on) and dissociation (k_off) rate constants for binding each partner. The substitutions had little impact on association rates and resulted in a wide range of dissociation rates as previously observed for substitutions on the BLIP side of the interface. The substitutions also had less effect on binding affinity for BLIP than BLIP-II. This is consistent with the high affinity and small binding interface of the TEM-1-BLIP-II complex, which predicts per residue contributions should be higher for TEM-1 binding to BLIP-II versus BLIP. Two TEM-1 residues (E104 and M129) were found to be hotspots for binding BLIP while five (L102, Y105, P107, K111, and M129) are hotspots for binding BLIP-II with only M129 as a common hotspot for both. Thus, although the same TEM-1 surface binds to both BLIP and BLIP-II, the distribution of binding energy on the surface is different for the two target proteins, that is, different binding strategies are employed.     

Statistics and Physical Origins of pK and Ionization State Changes upon Protein-Ligand Binding

This work investigates statistical prevalence and overall physical origins of changes in charge states of receptor proteins upon ligand binding. These changes are explored as a function of the ligand type (small molecule, protein, and nucleic acid), and distance from the binding region. Standard continuum solvent methodology is used to compute, on an equal footing, pK changes upon ligand binding for a total of 5899 ionizable residues in 20 protein-protein, 20 protein-small molecule, and 20 protein-nucleic acid high-resolution complexes. The size of the data set combined with an extensive error and sensitivity analysis allows us to make statistically justified and conservative conclusions: in 60% of all protein-small molecule, 90% of all protein-protein, and 85% of all protein-nucleic acid complexes there exists at least one ionizable residue that changes its charge state upon ligand binding at physiological conditions (pH = 6.5). Considering the most biologically relevant pH range of 4-8, the number of ionizable residues that experience substantial pK changes (ΔpK > 1.0) due to ligand binding is appreciable: on average, 6% of all ionizable residues in protein-small molecule complexes, 9% in protein-protein, and 12% in protein-nucleic acid complexes experience a substantial pK change upon ligand binding. These changes are safely above the statistical false-positive noise level. Most of the changes occur in the immediate binding interface region, where approximately one out of five ionizable residues experiences substantial pK change regardless of the ligand type. However, the physical origins of the change differ between the types: in protein-nucleic acid complexes, the pK values of interface residues are predominantly affected by electrostatic effects, whereas in protein-protein and protein-small molecule complexes, structural changes due to the induced-fit effect play an equally important role. In protein-protein and protein-nucleic acid complexes, there is a statistically significant number of substantial pK perturbations, mostly due to the induced-fit structural changes, in regions far from the binding interface.     

An ultra‐high affinity protein–protein interface displaying sequence‐robustness

Protein–protein interactions are crucial in biology and play roles in for example, the immune system, signaling pathways, and enzyme regulation. Ultra‐high affinity interactions (K _d <0.1 nM) occur in these systems, however, structures and energetics behind stability of ultra‐high affinity protein–protein complexes are not well understood. Regulation of the starch debranching barley limit dextrinase (LD) and its endogenous cereal type inhibitor (LDI) exemplifies an ultra‐high affinity complex (K _d of 42 pM). In this study the LD–LDI complex is investigated to unveil how robust the ultra‐high affinity is to LDI sequence variation at the protein–protein interface and whether alternative sequences can retain the ultra‐high binding affinity. The interface of LD–LDI was engineered using computational protein redesign aiming at identifying LDI variants predicted to retain ultra‐high binding affinity. These variants present a very diverse set of mutations going beyond conservative and alanine substitutions typically used to probe interfaces. Surface plasmon resonance analysis of the LDI variants revealed that high affinity of LD–LDI requires interactions of several residues at the rim of the protein interface, unlike the classical hotspot arrangement where key residues are found at the center of the interface. Notably, substitution of interface residues in LDI, including amino acids with functional groups different from the wild‐type, could occur without loss of affinity. This demonstrates that ultra‐high binding affinity can be conferred without hotspot residues, thus making complexes more robust to mutational drift in evolution. The present mutational analysis also demonstrates how energetic coupling can emerge between residues at large distances at the interface.     

Non‐interacting surface solvation and dynamics in protein–protein interactions

Protein–protein interactions control a plethora of cellular processes, including cell proliferation, differentiation, apoptosis, and signal transduction. Understanding how and why proteins interact will inevitably lead to novel structure‐based drug design methods, as well as design of de novo binders with preferred interaction properties. At a structural and molecular level, interface and rim regions are not enough to fully account for the energetics of protein–protein binding, even for simple lock‐and‐key rigid binders. As we have recently shown, properties of the global surface might also play a role in protein–protein interactions. Here, we report on molecular dynamics simulations performed to understand solvent effects on protein–protein surfaces. We compare properties of the interface, rim, and non‐interacting surface regions for five different complexes and their free components. Interface and rim residues become, as expected, less mobile upon complexation. However, non‐interacting surface appears more flexible in the complex. Fluctuations of polar residues are always lower compared with charged ones, independent of the protein state. Further, stable water molecules are often observed around polar residues, in contrast to charged ones. Our analysis reveals that (a) upon complexation, the non‐interacting surface can have a direct entropic compensation for the lower interface and rim entropy and (b) the mobility of the first hydration layer, which is linked to the stability of the protein–protein complex, is influenced by the local chemical properties of the surface. These findings corroborate previous hypotheses on the role of the hydration layer in shielding protein–protein complexes from unintended protein–protein interactions. Proteins 2015; 83:445–458. © 2014 Wiley Periodicals, Inc.     

The distribution of residues in a polypeptide sequence is a determinant of aggregation optimized by evolution

It has been shown that the propensity of a protein to form amyloid-like fibrils can be predicted with high accuracy from the knowledge of its amino acid sequence. It has also been suggested, however, that some regions of the sequences are more important than others in determining the aggregation process. Here, we have addressed this issue by constructing a set of “sequence scrambled” variants of the first 29 residues of horse heart apomyoglobin (apoMb_1-29), in which the sequence was modified while maintaining the same amino acid composition. The clustering of the most amyloidogenic residues in one region of the sequence was found to cause a marked increase of the elongation rate (k_agg) and a remarkable shortening of the lag phase (t_lag) of the fibril growth, as determined by far-UV circular dichroism and thioflavin T fluorescence. We also show that taking explicitly into consideration the presence of aggregation-promoting regions in the predictive methods results in a quantitative agreement between the theoretical and observed k_agg and t_lag values of the apoMb_1-29 variants. These results, together with a comparison between homologous segments from the family of globins, indicate the existence of a negative selection against the clustering of highly amyloidogenic residues in one or few regions of polypeptide sequences.     

KFC2: a knowledge-based hot spot prediction method based on interface solvation, atomic density, and plasticity features

Hot spots constitute a small fraction of protein-protein interface residues, yet they account for a large fraction of the binding affinity. Based on our previous method (KFC), we present two new methods (KFC2a and KFC2b) that outperform other methods at hot spot prediction. A number of improvements were made in developing these new methods. First, we created a training data set that contained a similar number of hot spot and non-hot spot residues. In addition, we generated 47 different features, and different numbers of features were used to train the models to avoid over-fitting. Finally, two feature combinations were selected: One (used in KFC2a) is composed of eight features that are mainly related to solvent accessible surface area and local plasticity; the other (KFC2b) is composed of seven features, only two of which are identical to those used in KFC2a. The two models were built using support vector machines (SVM). The two KFC2 models were then tested on a mixed independent test set, and compared with other methods such as Robetta, FOLDEF, HotPoint, MINERVA, and KFC. KFC2a showed the highest predictive accuracy for hot spot residues (True Positive Rate: TPR = 0.85); however, the false positive rate was somewhat higher than for other models. KFC2b showed the best predictive accuracy for hot spot residues (True Positive Rate: TPR = 0.62) among all methods other than KFC2a, and the False Positive Rate (FPR = 0.15) was comparable with other highly predictive methods.     

Crystal structure of glucose-6-phosphate isomerase from Thermus thermophilus HB8 showing a snapshot of active dimeric state

Glucose-6-phosphate isomerase (GPI) is a glycolytic enzyme with ill-defined oligomeric state. In order to obtain insight into the correlation between oligomerization and the catalytic function of this enzyme, the crystal structure of GPI from the extreme thermophile Thermus thermophilus HB8 (TtGPI) has been determined at 1.95 Å resolution. The crystallographic asymmetric unit contains an apparent dimer. The core fold of protomer and the interprotomer spatial arrangement of the dimer are similar to those of already reported crystal structures of other GPIs. The active site is located on the dimer interface, and putative catalytic residues are well conserved among the GPIs. These results suggest that the observed dimeric state of TtGPI in the crystal is biologically relevant and that this enzyme uses a common catalytic mechanism for the isomerase reaction. Gel-filtration chromatography, chemical cross-linking, sedimentation equilibrium by analytical ultracentrifugation, and dynamic light-scattering experiments indicate that TtGPI exists in a dynamic equilibrium between monomeric and dimeric states in solution. Several factors potentially contributing to the thermal stability of TtGPI protomer were identified: (i) a decrease in denaturation entropy by the shorter polypeptide length and by amino acid composition, including the increased number of proline residues and a higher arginine-to-lysine ratio; (ii) a larger number of ion pairs; and (iii) a reduction in cavity volume. From these results, it is suggested that transient dimer formation is sufficient for the catalytic function and that the TtGPI protomer itself has intrinsically higher thermal stability.     

Hydrophobic effect on the stability and folding of a hyperthermophilic protein

Ribonuclease HII from hyperthermophile Thermococcus kodakaraensis (Tk-RNase HII) is a kinetically robust monomeric protein. The conformational stability and folding kinetics of Tk-RNase HII were measured for nine mutant proteins in which a buried larger hydrophobic side chain is replaced by a smaller one (Leu/Ile to Ala). The mutant proteins were destabilized by 8.9 to 22.0 kJ mol^− 1 as compared with the wild-type protein. The removal of each -CH₂- group burial decreased the stability by 5.1 kJ mol^− 1 on average in the mutant proteins of Tk-RNase HII examined. This is comparable with the value of 5.3 kJ mol^− 1 obtained from experiments for proteins from organisms growing at moderate temperature. We conclude that the hydrophobic residues buried inside protein molecules contribute to the stabilization of hyperthermophilic proteins to a similar extent as proteins at normal temperature. In the folding experiments, the mutant proteins of Tk-RNase HII examined exhibited faster unfolding compared with the wild-type protein. These results indicate that the buried hydrophobic residues strongly contribute to the kinetic robustness of Tk-RNase HII. This is the first report that provides a practical cause of slow unfolding of hyperthermostable proteins.     

Protein-protein interaction site predictions with three-dimensional probability distributions of interacting atoms on protein surfaces

Protein-protein interactions are key to many biological processes. Computational methodologies devised to predict protein-protein interaction (PPI) sites on protein surfaces are important tools in providing insights into the biological functions of proteins and in developing therapeutics targeting the protein-protein interaction sites. One of the general features of PPI sites is that the core regions from the two interacting protein surfaces are complementary to each other, similar to the interior of proteins in packing density and in the physicochemical nature of the amino acid composition. In this work, we simulated the physicochemical complementarities by constructing three-dimensional probability density maps of non-covalent interacting atoms on the protein surfaces. The interacting probabilities were derived from the interior of known structures. Machine learning algorithms were applied to learn the characteristic patterns of the probability density maps specific to the PPI sites. The trained predictors for PPI sites were cross-validated with the training cases (consisting of 432 proteins) and were tested on an independent dataset (consisting of 142 proteins). The residue-based Matthews correlation coefficient for the independent test set was 0.423; the accuracy, precision, sensitivity, specificity were 0.753, 0.519, 0.677, and 0.779 respectively. The benchmark results indicate that the optimized machine learning models are among the best predictors in identifying PPI sites on protein surfaces. In particular, the PPI site prediction accuracy increases with increasing size of the PPI site and with increasing hydrophobicity in amino acid composition of the PPI interface; the core interface regions are more likely to be recognized with high prediction confidence. The results indicate that the physicochemical complementarity patterns on protein surfaces are important determinants in PPIs, and a substantial portion of the PPI sites can be predicted correctly with the physicochemical complementarity features based on the non-covalent interaction data derived from protein interiors.     

Prediction of Protein-Protein Interfaces on G-Protein β Subunits Reveals a Novel Phospholipase C β2 Binding Domain

Gβ subunits from heterotrimeric G-proteins (guanine nucleotide-binding proteins) directly bind diverse proteins, including effectors and regulators, to modulate a wide array of signaling cascades. These numerous interactions constrained the evolution of the molecular surface of Gβ. Although mammals contain five Gβ genes comprising two classes (Gβ1-like and Gβ5-like), plants and fungi have a single ortholog, and organisms such as Caenorhabditis elegans and Drosophila melanogaster contain one copy from each class. A limited number of crystal structures of complexes containing Gβ subunits and complementary biochemical data highlight specific sites within Gβs needed for protein interactions. It is difficult to determine from these interaction sites what, if any, additional regions of the Gβ molecular surface comprise interaction interfaces essential to Gβ's role as a nexus in numerous signaling cascades. We used a comparative evolutionary approach to identify five known and eight previously unknown putative interfaces on the surface of Gβ. We show that one such novel interface occurs between Gβ and phospholipase C β2 (PLC-β2), a mammalian Gβ interacting protein. Substitutions of residues within this Gβ-PLC-β2 interface reduce the activation of PLC-β2 by Gβ1, confirming that our de novo comparative evolutionary approach predicts previously unknown Gβ-protein interfaces. Similarly, we hypothesize that the seven remaining untested novel regions contribute to putative interfaces for other Gβ interacting proteins. Finally, this comparative evolutionary approach is suitable for application to any protein involved in a significant number of protein-protein interactions.     

Crystal structure and biochemical characterization of a manganese superoxide dismutase from Chaetomium thermophilum

A manganese superoxide dismutase from the thermophilic fungus Chaetomium thermophilum (CtMnSOD) was expressed in Pichia pastoris and purified to homogeneity. Its optimal temperature was 60 °C with approximately 75% of its activity retained after incubation at 70 °C for 60 min. Recombinant yeast cells carrying C. thermophilum mnsod gene exhibited higher stress resistance to salt and oxidative stress-inducing agents than control yeast cells. In an effort to provide structural insights, CtMnSOD was crystallized and its structure was determined at 2.0 Å resolution. The overall architecture of CtMnSOD was found similar to other MnSODs with highest structural similarities obtained against a MnSOD from the thermotolerant fungus Aspergillus fumigatus. In order to explain its thermostability, structural and sequence analysis of CtMnSOD with other MnSODs was carried out. An increased number of charged residues and an increase in the number of intersubunit salt bridges and the Thr:Ser ratio were identified as potential reasons for the thermostability of CtMnSOD.