Structural adaptation of the subunit interface of oligomeric thermophilic and hyperthermophilic enzymes

The activity and stability of a cold-active acylaminoacyl peptidase rely on its dimerization by domain swapping

The study of enzymes from extremophiles arouses interest in Protein Science because of the amazing solutions these proteins adopt to cope with extreme conditions. Recently solved, the structure of the psychrophilic acyl aminoacyl peptidase from Sporosarcina psychrophila (SpAAP) pinpoints a mechanism of dimerization unusual for this class of enzymes. The quaternary structure of SpAAP relies on a domain-swapping mechanism involving the N-terminal A1 helix. The A1 helix is conserved among homologous mesophilic and psychrophilic proteins and its deletion causes the formation of a monomeric enzyme, which is inactive and prone to aggregate. Here, we investigate the dimerization mechanism of SpAAP through the analysis of chimeric heterodimers where a protomer lacking the A1 helix combines with a protomer carrying the inactivated catalytic site. Our results indicate that the two active sites are independent, and that a single A1 helix is sufficient to partially recover the quaternary structure and the activity of chimeric heterodimers. Since catalytically competent protomers are unstable and inactive unless they dimerize, SpAAP reveals as an "obligomer" for both structural and functional reasons.

The type IV pilus assembly motor PilB is a robust hexameric ATPase with complex kinetics

The bacterial type IV pilus (T4P) is a versatile nanomachine that functions in pathogenesis, biofilm formation, motility, and horizontal gene transfer. T4P assembly is powered by the motor ATPase PilB which is proposed to hydrolyze ATP by a symmetrical rotary mechanism. This mechanism, which is deduced from the structure of PilB, is untested. Here, we report the first kinetic studies of the PilB ATPase, supporting co-ordination among the protomers of this hexameric enzyme. Analysis of the genome sequence of Chloracidobacterium thermophilum identified a pilB gene whose protein we then heterologously expressed. This PilB formed a hexamer in solution and exhibited highly robust ATPase activity. It displays complex steady-state kinetics with an incline followed by a decline over an ATP concentration range of physiological relevance. The incline is multiphasic and the decline signifies substrate inhibition. These observations suggest that variations in intracellular ATP concentrations may regulate T4P assembly and T4P-mediated functions in vivo in accordance with the physiological state of bacteria with unanticipated complexity. We also identified a mutant pilB gene in the genomic DNA of C. thermophilum from an enrichment culture. The mutant PilB variant, which is significantly less active, exhibited similar inhibition of its ATPase activity by high concentrations of ATP. Our findings here with the PilB ATPase from C. thermophilum provide the first line of biochemical evidence for the co-ordination among PilB protomers consistent with the symmetrical rotary model of catalysis based on structural studies.

Physical and molecular bases of protein thermal stability and cold adaptation

The molecular bases of thermal and cold stability and adaptation, which allow proteins to remain folded and functional in the temperature ranges in which their host organisms live and grow, are still only partially elucidated. Indeed, both experimental and computational studies fail to yield a fully precise and global physical picture, essentially because all effects are context-dependent and thus quite intricate to unravel. We present a snapshot of the current state of knowledge of this highly complex and challenging issue, whose resolution would enable large-scale rational protein design.

Abundance and Temperature Dependency of Protein-Protein Interaction Revealed by Interface Structure Analysis and Stability Evolution

Protein complexes are major forms of protein-protein interactions and implement essential biological functions. The subunit interface in a protein complex is related to its thermostability. Though the roles of interface properties in thermal adaptation have been investigated for protein complexes, the relationship between the interface size and the expression level of the subunits remains unknown. In the present work, we studied this relationship and found a positive correlation in thermophiles rather than mesophiles. Moreover, we found that the protein interaction strength in complexes is not only temperature-dependent but also abundance-dependent. The underlying mechanism for the observed correlation was explored by simulating the evolution of protein interface stability, which highlights the avoidance of misinteraction. Our findings make more complete the picture of the mechanisms for protein complex thermal adaptation and provide new insights into the principles of protein-protein interactions.

Low Temperature Adaptation Is Not the Opposite Process of High Temperature Adaptation in Terms of Changes in Amino Acid Composition

Previous studies focused on psychrophilic adaptation generally have demonstrated that multiple mechanisms work together to increase protein flexibility and activity, as well as to decrease the thermostability of proteins. However, the relationship between high and low temperature adaptations remains unclear. To investigate this issue, we collected the available predicted whole proteome sequences of species with different optimal growth temperatures, and analyzed amino acid variations and substitutional asymmetry in pairs of homologous proteins from related species. We found that changes in amino acid composition associated with low temperature adaptation did not exhibit a coherent opposite trend when compared with changes in amino acid composition associated with high temperature adaptation. This result indicates that during their evolutionary histories the proteome-scale evolutionary patterns associated with prokaryotes exposed to low temperature environments were distinct from the proteome-scale evolutionary patterns associated with prokaryotes exposed to high temperature environments in terms of changes in amino acid composition of the proteins.

The crystal structure of archaeal serine hydroxymethyltransferase reveals idiosyncratic features likely required to withstand high temperatures

Serine hydroxymethyltransferases (SHMTs) play an essential role in one-carbon unit metabolism and are used in biomimetic reactions. We determined the crystal structure of free (apo) and pyridoxal-5'-phosphate-bound (holo) SHMT from Methanocaldococcus jannaschii, the first from a hyperthermophile, from the archaea domain of life and that uses H₄MPT as a cofactor, at 2.83 and 3.0 Å resolution, respectively. Idiosyncratic features were observed that are likely to contribute to structure stabilization. At the dimer interface, the C-terminal region folds in a unique fashion with respect to SHMTs from eubacteria and eukarya. At the active site, the conserved tyrosine does not make a cation-π interaction with an arginine like that observed in all other SHMT structures, but establishes an amide-aromatic interaction with Asn257, at a different sequence position. This asparagine residue is conserved and occurs almost exclusively in (hyper)thermophile SHMTs. This led us to formulate the hypothesis that removal of frustrated interactions (such as the Arg-Tyr cation-π interaction occurring in mesophile SHMTs) is an additional strategy of adaptation to high temperature. Both peculiar features may be tested by designing enzyme variants potentially endowed with improved stability for applications in biomimetic processes.

Extremophilic SHMTs: from structure to biotechnology

Recent advances in molecular and structural biology have improved the availability of virtually any biocatalyst in large quantity and have also provided an insight into the detailed structure-function relationships of many of them. These results allowed the rational exploitation of biocatalysts for use in organic synthesis. In this context, extremophilic enzymes are extensively studied for their potential interest for many biotechnological and industrial applications, as they offer increased rates of reactions, higher substrate solubility, and/or longer enzyme half-lives at the conditions of industrial processes. Serine hydroxymethyltransferase (SHMT), for its ubiquitous nature, represents a suitable model for analyzing enzyme adaptation to extreme environments. In fact, many SHMT sequences from Eukarya, Eubacteria and Archaea are available in data banks as well as several crystal structures. In addition, SHMT is structurally conserved because of its critical metabolic role; consequently, very few structural changes have occurred during evolution. Our research group analyzed the molecular basis of SHMT adaptation to high and low temperatures, using experimental and comparative in silico approaches. These structural and functional studies of SHMTs purified from extremophilic organisms can help to understand the peculiarities of the enzyme activity at extreme temperatures, indicating possible strategies for rational enzyme engineering.

Conservation of functionally important global motions in an enzyme superfamily across varying quaternary structures

The α-d-phosphohexomutase superfamily comprises enzymes involved in carbohydrate metabolism that are found in all kingdoms of life. Recent biophysical studies have shown for the first time that several of these enzymes exist as dimers in solution, prompting an examination of the oligomeric state of all proteins of known structure in the superfamily (11 different proteins; 31 crystal structures) via computational and experimental analyses. We find that these proteins range in quaternary structure from monomers to tetramers, with 6 of the 11 known structures being likely oligomers. The oligomeric state of these proteins not only is associated in some cases with enzyme subgroup (i.e., substrate specificity) but also appears to depend on domain of life, with the two archaeal proteins existing as higher-order oligomers. Within the oligomers, three distinct interfaces are observed, one of which is found in both archaeal and bacterial proteins. Normal mode analysis shows that the topological arrangement of the oligomers permits domain 4 of each protomer to move independently as required for catalysis. Our analysis suggests that the advantages associated with protein flexibility in this enzyme family are of sufficient importance to be maintained during the evolution of multiple independent oligomers. This study is one of the first showing that global motions may be conserved not only within protein families but also across members of a superfamily with varying oligomeric structures.

Regulation of enzyme activity of alcohol dehydrogenase through its interactions with pyruvate-ferredoxin oxidoreductase in Thermoanaerobacter tengcongensis

Alcohol dehydrogenases (ADHs) from thermophilic microorganisms are interesting enzymes that have their potential applications in biotechnology and potentially provide insight into the mechanisms of action of thermo-tolerant proteins. The molecular mechanisms of ADHs under thermal stress in vivo have yet to be explored. Herein, we employed a proteomic strategy to survey the possible interactions of secondary-ADH (2-ADH) with other proteins in Thermoanaerobacter tengcongensis (T. tengcongensis) cultured at 75°C and found that 2-ADH, pyruvate-ferredoxin oxidoreductase (PFOR) and several glycolytic enzymes coexisted in a protein complex. Using anion exchange chromatography, the elution profile indicated that the native 2-ADH was present in two forms, PFOR-bound and PFOR-free. Immuno-precipitation and pull down analysis further validated the interactions between 2-ADH and PFOR. The kinetic behaviours of 2-ADH either in the recombinant or native form were evaluated with different substrates. The enzyme activity of 2-ADH was inhibited in a non-competitive mode by PFOR, implying the interaction of 2-ADH and PFOR negatively regulated alcohol formation. In T. tengcongensis, PFOR is an enzyme complex located at the upstream of 2-ADH in the alcohol generation pathway. These findings, therefore, offered a plausible mechanism for how alcohol metabolism is regulated by hetero-interactions between 2-ADH and PFOR, especially in anaerobic thermophiles.

Contribution of inter-subunit interactions to the thermostability of Pyrococcus furiosus citrate synthase

Using citrate synthase from the hyperthermophile Pyrococcus furiosus (PfCS) as our test molecule, we show through guanidine hydrochloride-induced unfolding that the dimer separates into folded, but inactive, monomers before individual subunit unfolding takes place. Given that forces across the dimer interface are vital for thermostability, a robust computational method was derived that uses the University of Houston Brownian Dynamics (UHBD) program to calculate both the hydrophobic and electrostatic contribution to the dimerisation energy at 100°C. The results from computational and experimental determination of the lowered stability of interface mutants were correlated, being both of the same order of magnitude and placing the mutant proteins in the same order of stability. This computational method, optimised for hyperthermophilic molecules and tested in the laboratory, after further testing on other examples, could be of widespread use in the prediction of thermostabilising mutations in other oligomeric proteins for which dissociation is the first step in unfolding.

Structure and flexibility in cold-adapted iron superoxide dismutases: the case of the enzyme isolated from Pseudoalteromonas haloplanktis

Superoxide dismutases (SODs) are metalloenzymes catalysing the dismutation of superoxide anion radicals into molecular oxygen and hydrogen peroxide. Here, we present the crystal structure of a cold-adapted Fe-SOD from the Antarctic eubacterium Pseudoalteromonas haloplanktis (PhSOD), and that of its complex with sodium azide. The structures were compared with those of the corresponding homologues having a high sequence identity with PhSOD, such as the mesophilic SOD from Escherichia coli (EcSOD) or Pseudomonas ovalis, and the psychrophilic SOD from Aliivibrio salmonicida (AsSOD). These enzymes shared a large structural similarity, such as a conserved tertiary structure and arrangement of the two monomers, an almost identical total number of inter- and intramolecular hydrogen bonds and salt bridges. However, the two cold-adapted SODs showed an increased flexibility of the active site residues with respect to their mesophilic homologues. Structural information was combined with a characterisation of the chemical and thermal stability performed by CD and fluorescence measurements. Despite of its psychrophilic origin, the denaturation temperature of PhSOD was comparable with that of the mesophilic EcSOD, whereas AsSOD showed a lower denaturation temperature. On the contrary, the values of the denaturant concentration at the transition midpoint were in line with the psychrophilic/mesophilic origin of the proteins. These data provide additional support to the hypothesis that cold-adapted enzymes achieve efficient catalysis at low temperature, by increasing the flexibility of their active site; moreover, our results underline how fine structural modifications can alter enzyme flexibility and/or stability without compromising the overall structure of typical rigid enzymes, such as SODs.

A novel scoring function for discriminating hyperthermophilic and mesophilic proteins with application to predicting relative thermostability of protein mutants

Background

The ability to design thermostable proteins is theoretically important and practically useful. Robust and accurate algorithms, however, remain elusive. One critical problem is the lack of reliable methods to estimate the relative thermostability of possible mutants.

Results

We report a novel scoring function for discriminating hyperthermophilic and mesophilic proteins with application to predicting the relative thermostability of protein mutants. The scoring function was developed based on an elaborate analysis of a set of features calculated or predicted from 540 pairs of hyperthermophilic and mesophilic protein ortholog sequences. It was constructed by a linear combination of ten important features identified by a feature ranking procedure based on the random forest classification algorithm. The weights of these features in the scoring function were fitted by a hill-climbing algorithm. This scoring function has shown an excellent ability to discriminate hyperthermophilic from mesophilic sequences. The prediction accuracies reached 98.9% and 97.3% in discriminating orthologous pairs in training and the holdout testing datasets, respectively. Moreover, the scoring function can distinguish non-homologous sequences with an accuracy of 88.4%. Additional blind tests using two datasets of experimentally investigated mutations demonstrated that the scoring function can be used to predict the relative thermostability of proteins and their mutants at very high accuracies (92.9% and 94.4%). We also developed an amino acid substitution preference matrix between mesophilic and hyperthermophilic proteins, which may be useful in designing more thermostable proteins.

Conclusions

We have presented a novel scoring function which can distinguish not only HP/MP ortholog pairs, but also non-homologous pairs at high accuracies. Most importantly, it can be used to accurately predict the relative stability of proteins and their mutants, as demonstrated in two blind tests. In addition, the residue substitution preference matrix assembled in this study may reflect the thermal adaptation induced substitution biases. A web server implementing the scoring function and the dataset used in this study are freely available at http://www.abl.ku.edu/thermorank/.