Crystal structure and amide H/D exchange of binary complexes of alcohol dehydrogenase from Bacillus stearothermophilus: insight into thermostability and cofactor binding

The crystal structure of NAD(+)-dependent alcohol dehydrogenase from Bacillus stearothermophilus strain LLD-R (htADH) was determined using X-ray diffraction data at a resolution of 2.35 A. The structure of homotetrameric htADH is highly homologous to those of bacterial and archaeal homotetrameric alcohol dehydrogenases (ADHs) and also to the mammalian dimeric ADHs. There is one catalytic zinc atom and one structural zinc atom per enzyme subunit. The enzyme was crystallized as a binary complex lacking the nicotinamide adenine dinucleotide (NAD(+)) cofactor but including a zinc-coordinated substrate analogue trifluoroethanol. The binary complex structure is in an open conformation similar to ADH structures without the bound cofactor. Features important for the thermostability of htADH are suggested by a comparison with a homologous mesophilic enzyme (55% identity), NAD(+)-dependent alcohol dehydrogenase from Escherichia coli. To gain insight into the conformational change triggered by NAD(+) binding, amide hydrogen-deuterium exchange of htADH, in the presence and absence of NAD(+), was studied by HPLC-coupled electrospray mass spectrometry. When the deuteron incorporation of the protein-derived peptides was analyzed, it was found that 9 of 21 peptides show some decrease in the level of deuteron incorporation upon NAD(+) binding, and another 4 peptides display slower exchange rates. With one exception (peptide number 8), none of the peptides that are altered by bound NAD(+) are in contact with the alcohol-substrate-binding pocket. Furthermore, peptides 5 and 8, which are located outside the NAD(+)-binding pocket, are notable by displaying changes upon NAD(+) binding. This suggests that the transition from the open to the closed conformation caused by cofactor binding has some long-range effects on the protein structure and dynamics.     

Important inter-residue contacts for enhancing the thermal stability of thermophilic proteins

Proteins from thermophilic organisms exhibit high thermal stability, but have structures that are very similar to their mesophilic homologues. In order to gain insight into the basis of thermostability, we have analyzed the medium- and long-range contacts in mesophilic and thermophilic proteins of 16 different families. We found that the thermophiles prefer to have contacts between residues with hydrogen-bond-forming capability. Apart from hydrophobic contacts, more contacts are observed between polar and non-polar residues in thermophiles than mesophiles. Residue-wise analysis showed that Tyr has good contacts with several other residues, and Cys has considerably higher long-range contacts in thermophiles compared with mesophiles. Furthermore, the residues occurring in the range of 31-34 residues apart in the sequence contribute significant long-range contacts to the stability of thermophilic proteins.     

Discrimination of mesophilic and thermophilic proteins using machine learning algorithms

Discriminating thermophilic proteins from their mesophilic counterparts is a challenging task and it would help to design stable proteins. In this work, we have systematically analyzed the amino acid compositions of 3075 mesophilic and 1609 thermophilic proteins belonging to 9 and 15 families, respectively. We found that the charged residues Lys, Arg, and Glu as well as the hydrophobic residues, Val and Ile have higher occurrence in thermophiles than mesophiles. Further, we have analyzed the performance of different methods, based on Bayes rules, logistic functions, neural networks, support vector machines, decision trees and so forth for discriminating mesophilic and thermophilic proteins. We found that most of the machine learning techniques discriminate these classes of proteins with similar accuracy. The neural network-based method could discriminate the thermophiles from mesophiles at the five-fold cross-validation accuracy of 89% in a dataset of 4684 proteins. Moreover, this method is tested with 325 mesophiles in Xylella fastidosa and 382 thermophiles in Aquifex aeolicus and it could successfully discriminate them with the accuracy of 91%. These accuracy levels are better than other methods in the literature and we suggest that this method could be effectively used to discriminate mesophilic and thermophilic proteins.     

Patterns of temperature adaptation in proteins from Methanococcus and Bacillus

It has long been known that amino acid substitutions in proteins of organisms living at moderate and high temperatures (mesophiles and thermophiles, respectively) are not all symmetrical; for example, more aligned sites have lysine in mesophiles and arginine in thermophiles than have the opposite pattern. This is generally taken to indicate that certain amino acids are favored over others by selection at different temperatures. Previous comparisons of protein sequences from mesophiles and thermophiles have used relatively small numbers of sequences from a diverse array of species, meaning that only the most common amino acid substitutions could be examined and any taxon-specific patterns would be obscured. Here, we compare a large number of proteins between mesophiles and thermophiles in the archaeal genus Methanococcus and the bacterial genus Bacillus. Each genus exhibits dramatically asymmetrical substitution patterns for many pairs of amino acids. There are several pairs of amino acids for which one amino acid is favored in thermophilic Bacillus and the other is favored in thermophilic Methanococcus; this appears to result from the higher G + C content of the DNA of thermophilic Bacillus, a complication not seen in Methanococcus.     

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.     

The stability of thermophilic proteins: a study based on comprehensive genome comparison

We address the question of the thermal stability of proteins in thermophiles through comprehensive genome comparison, focussing on the occurrence of salt bridges. We compared a set of 12 genomes (from four thermophilic archaeons, one eukaryote, six mesophilic eubacteria, and one thermophilic eubacteria). Our results showed that thermophiles have a greater content of charged residues than mesophiles, both at the overall genomic level and in alpha helices. Furthermore, we found that in thermophiles the charged residues in helices tend to be preferentially arranged with a 1–4 helical spacing and oriented so that intra-helical charge pairs agree with the helix dipole. Collectively, these results imply that intra-helical salt bridges are more prevalent in thermophiles than mesophiles and thus suggest that they are an important factor stabilizing thermophilic proteins. We also found that the proteins in thermophiles appear to be somewhat shorter than those in mesophiles. However, this later observation may have more to do with evolutionary relationships than with physically stabilizing factors. In all our statistics we were careful to controls for various biases. These could have, for instance, arisen due to repetitive or duplicated sequences. In particular, we repeated our calculation using a variety of random and directed sampling schemes. One of these involved making a "stratified sample," a representative cross-section of the genomes derived from a set of 52 orthologous proteins present roughly once in each genome. For another sample, we focused on the subset of the 52 orthologs that had a known 3D structure. This allowed us to determine the frequency of tertiary as well as main-chain salt bridges. Our statistical controls supported our overall conclusion about the prevalence of salt bridges in thermophiles in comparison to mesophiles. Electronic Publication     

Protein thermal stability: insights from atomic displacement parameters (B values)

The factors contributing to the thermal stability of proteins from thermophilic origins are matters of intense debate and investigation. Thermophilic proteins are thought to possess better packed interiors than their mesophilic counterparts, leading to lesser overall flexibility and a corresponding reduction in surface-to-volume ratio. These observations prompted an analysis of B values reported in high-resolution X-ray crystal structures of mesophilic and thermophilic proteins. In this analysis, the following aspects were addressed: (1) frequency distribution of normalized B values (B' factors) over all the proteins and for individual amino acids; (2) amino acid compositions in high B value regions of polypeptide chains; (3) variation in the B values from core to the surface of proteins in terms of their radius of gyration; and (4) degree of dispersion of normalized B values in spheres around the Calpha atoms. The analysis revealed that (1) Ser and Thr have lesser flexibility in thermophiles than in mesophiles, (2) the proportion of Glu and Lys in high B value regions of thermophiles is higher and that of Ser and Thr is lower and (3) the dispersion of B values within spheres at Calpha atoms is similar in mesophiles and thermophiles. These observations reflect plausible differences in the dynamics of thermophilic and mesophilic proteins and suggest amino acid substitutions that are likely to change thermal stability.     

Simple efficient methods for the isolation of malate dehydrogenase from thermophilic and mesophilic bacteria.

Malate dehydrogenase from a number of bacteria drawn from several genera and representing the mesophilic, moderately thermophilic and extremely thermophilic classes was isolated by procedures which involve only a small number of steps (in most cases only two), of which the key one is affinity chromatography on 5'-AMP--Sepharose and/or on NAD+--hexane--agarose. Electrophoretic analysis of the native enzymes in polyacrylamide gel and of the denaturated enzymes in sodium dodecyl sulphate/polyacrylamide gel revealed no significant protein impurity in the purified preparations. The yields ranged from about 40% to over 80%. The malate dehydrogenases from the extreme thermophiles and from some of the moderate thermophiles are appreciably less efficient catalytically than their mesophilic homologues.     

Structural differences between thermophilic and mesophilic membrane proteins

The evolutionary adaptations of thermophilic water‐soluble proteins required for maintaining stability at high temperature have been extensively investigated. Little is known about the adaptations in membrane proteins, however. Here, we compare many properties of mesophilic and thermophilic membrane protein structures, including side‐chain burial, packing, hydrogen bonding, transmembrane kinks, loop lengths, hydrophobicity, and other sequence features. Most of these properties are quite similar between mesophiles and thermophiles although we observe a slight increase in side‐chain burial and possibly a slight decrease in the frequency of transmembrane kinks in thermophilic membrane protein structures. The most striking difference is the increased hydrophobicity of thermophilic transmembrane helices, possibly reflecting more stringent hydrophobicity requirements for membrane partitioning at high temperature. In agreement with prior work examining transmembrane sequences, we find that thermophiles have an increase in small residues (Gly, Ala, Ser, and Val) and a strong suppression of Cys. We also find a relative dearth of most strongly polar residues (Asp, Asn, Glu, Gln, and Arg). These results suggest that in thermophiles, there is significant evolutionary pressure to offload destabilizing polar amino acids, to decrease the entropy cost of side chain burial, and to eliminate thermally sensitive amino acids.     

Crystal structure of the autocatalytic initiator of glycogen biosynthesis,glycogenin

The crystal structure of a medium-chain NAD(H)-dependent alcohol dehydrogenase (ADH) from an archaeon has been solved by multiwavelength anomalous diffraction, using a selenomethionine-substituted enzyme. The protein (SsADH), extracted from the hyperthermophilic organism Sulfolobus solfataricus, is a homo-tetramer with a crystallographic 222 symmetry. Despite the low level of sequence identity, the overall fold of the monomer is similar to that of the other homologous ADHs of known structure. However, a significant difference is the orientation of the catalytic domain relative to the coenzyme-binding domain that results in a larger interdomain cleft. At the bottom of this cleft, the catalytic zinc ion is coordinated tetrahedrally and lacks the zinc-bound water molecule that is usually found in ADH apoform structures. The fourth coordination position is indeed occupied by a Glu residue, as found in bacterial tetrameric ADHs. Other differences are found in the architecture of the substrate pocket whose entrance is more restricted than in other ADHs. SsADH is the first tetrameric ADH X-ray structure containing a second zinc ion playing a structural role. This latter metal ion shows a peculiar coordination, with a glutamic acid residue replacing one of the four cysteine ligands that are highly conserved throughout the structural zinc-containing dimeric ADHs.     

Role of cation-pi interactions to the stability of thermophilic proteins

Elucidating the factors responsible for exhibiting extreme thermal stability of thermophilic proteins is very important for an understanding of the mechanism of protein stability, as well as to design stable proteins. In this work, we have analyzed the influence of cation-pi interactions to enhance the stability from mesophilic to thermophilic proteins. The favorable residue pairs forming such a system of interactions have been brought out. We found that the Tyr has a greater number of such interactions with Lys in thermophilic proteins. Specifically, the same Lys would experience a greater number of cation-pi interactions with several Tyr residues in thermophiles. On the other hand, the influence of Phe in making cation-pi interactions is higher in mesophiles than in thermophiles. Further, a network of cation-pi interactions are maintained by Lys in thermophiles, whereas Arg plays a major role in mesophilic proteins. Moreover, atoms that have a substantial positive charge in both Lys and Arg make a more significant contribution for cation-pi interactions than do cationic group atoms.     

How Do Thermophilic Proteins and Proteomes Withstand High Temperature?

We attempt to understand the origin of enhanced stability in thermophilic proteins by analyzing thermodynamic data for 116 proteins, the largest data set achieved to date. We compute changes in entropy and enthalpy at the convergence temperature where different driving forces are maximally decoupled, in contrast to the majority of previous studies that were performed at the melting temperature. We find, on average, that the gain in enthalpy upon folding is lower in thermophiles than in mesophiles, whereas the loss in entropy upon folding is higher in mesophiles than in thermophiles. This implies that entropic stabilization may be responsible for the high melting temperature, and hints at residual structure or compactness of the denatured state in thermophiles. We find a similar trend by analyzing a homologous set of proteins classified based only on the optimum growth temperature of the organisms from which they were extracted. We find that the folding free energy at the temperature of maximal stability is significantly more favorable in thermophiles than in mesophiles, whereas the maximal stability temperature itself is similar between these two classes. Furthermore, we extend the thermodynamic analysis to model the entire proteome. The results explain the high optimal growth temperature in thermophilic organisms and are in excellent quantitative agreement with full thermal growth rate data obtained in a dozen thermophilic and mesophilic organisms.     

Patterns of temperature adaptation in proteins from the bacteria Deinococcus radiodurans and Thermus thermophilus

Asymmetrical patterns of amino acid substitution in proteins of organisms living at moderate and high temperatures (mesophiles and thermophiles, respectively) are generally taken to indicate selection favoring different amino acids at different temperatures due to their biochemical properties. If that were the case, comparisons of different pairs of mesophilic and thermophilic taxa would exhibit similar patterns of substitutional asymmetry. A previous comparison of mesophilic versus thermophilic Methanococcus with mesophilic versus thermophilic Bacillus revealed several pairs of amino acids for which one amino acid was favored in thermophilic Bacillus and the other was favored in thermophilic Methanococcus. Most of this could be explained by the higher G+C content of the DNA of thermophilic Bacillus, a phenomenon not seen in the Methanococcus comparison. Here, I compared the mesophilic bacterium Deinococcus radiodurans and its thermophilic relative Thermus thermophilus, which are similar in G+C content. Of the 190 pairs of amino acids, 83 exhibited significant substitutional asymmetry, consistent with the pervasive effects of selection. Most of these significantly asymmetrical pairs of amino acids were asymmetrical in the direction predicted from the Methanococcus data, consistent with thermal adaptation resulting from universal biochemical properties of the amino acids. However, 12 pairs of amino acids exhibited asymmetry significantly different from and in the opposite direction of that found in the Methanococcus comparison, and 21 pairs of amino acids exhibited asymmetry that was significantly different from that found in the Bacillus comparison and could not be explained by the greater G+C content in thermophilic Bacillus. This suggests that selection due to universal biochemical properties of the amino acids and differences in G+C content are not the only causes of substitutional asymmetry between mesophiles and thermophiles. Instead, selection on taxon-specific properties of amino acids, such as their metabolic cost, may play a role in causing asymmetrical patterns of substitution.     

Using a strategy based on the concept of convergent evolution to identify residue substitutions responsible for thermal adaptation

Factors that are related to thermostability of proteins have been extensively studied in recent years, especially by comparing thermophiles and mesophiles. However, most of them are global characters. It is still not clear how to identify specific residues or fragments which may be more relevant to protein thermostability. Moreover, some of the differences among the thermophiles and mesophiles may be due to phylogenetic differences instead of thermal adaptation. To resolve these problems, I adopted a strategy to identify residue substitutions evolved convergently in thermophiles or mesophiles. These residues may therefore be responsible for thermal adaptation. Four classes of genomes were utilized in this study, including thermophilic archaea, mesophilic archaea, thermophilic bacteria, and mesophilic bacteria. For most clusters of orthologous groups (COGs) with sequences from all of these four classes of genomes, I can identify specific residues or fragments that may potentially be responsible for thermal adaptation. Functional or structural constraints (represented as sequence conservation) were suggested to have higher impact on thermal adaptation than secondary structure or solvent accessibility does. I further compared thermophilic archaea and mesophilic bacteria, and found that the most diverged fragments may not necessarily correspond to the thermostability-determining ones. The usual approach to compare thermophiles and mesophiles without considering phylogenetic relationships may roughly identify sequence features contributing to thermostability; however, to specifically identify residue substitutions responsible for thermal adaptation, one should take sequence evolution into consideration.     

Genesis of Drosophila ADH: the shaping of the enzymatic activity from a SDR ancestor

Drosophila alcohol dehydrogenase (ADH) is an NAD(H)-dependent oxidoreductase that catalyzes the oxidation of alcohols and aldehydes. Structurally and biochemically distinct from all the reported ADHs (typically, the mammalian medium-chain dehydrogenase/reductase-ethanol-metabolizing enzyme), it stands as the only small-alcohol transforming system that has originated from a short-chain dehydrogenase/reductase (SDR) ancestor. The crystal structures of the apo, binary (E.NAD(+)) and three ternary (E.NAD(+).acetone, E.NAD(+).3-pentanone and E.NAD(+).cyclohexanone) forms of Drosophila lebanonensis ADH have allowed us to infer the structural and kinetic features accounting for the generation of the ADH activity within the SDR lineage.     

Natural alcohol exposure: is ethanol the main substrate for alcohol dehydrogenases in animals?

Alcohol dehydrogenase (ADH) activity is widely distributed in all phyla. In animals, three non-homologous NAD(P)(+)-dependent ADH protein families are reported. These arose independently throughout evolution and possess different structures and mechanisms of reaction: type I (medium-chain) ADHs are zinc-containing enzymes and comprise the most studied group in vertebrates; type II (short-chain) ADHs lack metal cofactor and have been extensively studied in Drosophila; and type III ADHs are iron-dependent/-activated enzymes that were initially identified only in microorganisms. The presence of these different ADHs in animals has been assumed to be a consequence of chronic exposure to ethanol. By far the most common natural source of ethanol is fermentation of fruit sugars by yeast, and available data support that this fruit trait evolved in concert with the characteristics of their frugivorous seed dispersers. Therefore, if the presence of ADHs in animals evolved as an adaptive response to dietary ethanol exposure, then it can be expected that the enzymogenesis of these enzymes began after the appearance of angiosperms with fleshy fruits, because substrate availability must precede enzyme selection. In this work, available evidence supporting this possibility is discussed. Phylogenetic analyses reveal that type II ADHs suffered several duplications, all of these restricted to flies (order Diptera). Induction of type II Adh by ethanol exposure, a positive correlation between ADH activity and ethanol resistance, and the fact that flies and type II Adh diversification occurred in concert with angiosperm diversification, strongly suggest that type II ADHs were recruited to allow larval flies to exploit new restricted niches with high ethanol content. In contrast, phyletic distribution of types I and III ADHs in animals showed that these appeared before angiosperms and land plants, independently of ethanol availability. Because these enzymes are not induced by ethanol exposure and possess a high affinity and/or catalytic efficiency for non-ethanol endogenous substrates, it can be concluded that the participation of types I and III ADHs in ethanol metabolism can be considered as incidental, and not adaptive.     

Regulation of coenzyme utilization by mitochondrial NAD(P)-dependent malic enzyme

1. Skeletal muscle mitochondrial NAD(P)-dependent malic enzyme [EC 1.1.1. 39, L-malate:NAD+ oxidoreductase (decarboxylating)] from herring could use both coenzymes, NAD and NADP, in a similar manner. 2. The coenzyme preference of mitochondrial NAD(P)-dependent malic enzyme was probed using dual wavelength spectroscopy and pairing the natural coenzymes, NAD or NADP with their respective thionicotinamide analogues, s-NADP or s-NAD, that have absorbance maxima in reduced forms at 400 nm. 3. s-NAD and s-NADP were found to be good alternate substrates for NAD(P)-dependent malic enzyme, the apparent Km values for the thioderivatives were similar to those of the corresponding natural coenzymes. 4. ATP produced greater inhibition of the NAD or s-NAD linked reactions than of the NADP or s-NADP-linked reactions of skeletal muscle mitochondrial NAD(P)-dependent malic enzyme. 5. At 5 mM malate concentration and in the presence of 2 mM ATP the NADP-linked reaction is favoured and the activity ratios, V(s-NADP)/V(NAD) or V(NADP)/V(s-NAD), are 6 and 26, respectively.     

Archaeal RecA homologues: different response to DNA-damaging agents in mesophilic and thermophilic Archaea

Two archaeal proteins, RadA and RadB, share similarity with the RecA/Rad51 family of recombinases, with RadA being the functional homologue. We have studied and compared the RadA and RadB proteins of mesophilic and thermophilic Archaea. In growing cells, RadA levels are similar in mesophilic Methanococcus species and the hyperthermophile Methanococcus jannaschii. Treatment of cells with mutagenic agents (methylmethane sulfonate or UV light) increased the expression of RadA (as evidenced by higher levels of both mRNA and protein) in all organisms tested, but the increase was greater in the mesophiles than in the thermophiles M. jannaschii and Sulfolobus solfataricus. Recombinantly expressed RadA proteins from the mesophile M. voltae and the thermophile M. jannaschii were similar in their ATPase- and DNA-binding activities. All the data are consistent with proposals that RadA plays the same role as eukaryotic Rad51. Surprisingly, the data also suggested that the thermophiles do not need more RadA protein or activity than the mesophiles. On the other hand, RadB is not coregulated with RadA, and its role remains unclear. Neither RadA nor RadB from a mesophile or from a thermophile rescued the UV-sensitive phenotype of an Escherichia coli recA- host.