Sequence and structural comparison of thermophilic phosphoglycerate kinases with a mesophilic equivalent

Carbohydrate metabolism in Archaea: current insights into unusual enzymes and pathways and their regulation

The metabolism of Archaea, the third domain of life, resembles in its complexity those of Bacteria and lower Eukarya. However, this metabolic complexity in Archaea is accompanied by the absence of many "classical" pathways, particularly in central carbohydrate metabolism. Instead, Archaea are characterized by the presence of unique, modified variants of classical pathways such as the Embden-Meyerhof-Parnas (EMP) pathway and the Entner-Doudoroff (ED) pathway. The pentose phosphate pathway is only partly present (if at all), and pentose degradation also significantly differs from that known for bacterial model organisms. These modifications are accompanied by the invention of "new," unusual enzymes which cause fundamental consequences for the underlying regulatory principles, and classical allosteric regulation sites well established in Bacteria and Eukarya are lost. The aim of this review is to present the current understanding of central carbohydrate metabolic pathways and their regulation in Archaea. In order to give an overview of their complexity, pathway modifications are discussed with respect to unusual archaeal biocatalysts, their structural and mechanistic characteristics, and their regulatory properties in comparison to their classic counterparts from Bacteria and Eukarya. Furthermore, an overview focusing on hexose metabolic, i.e., glycolytic as well as gluconeogenic, pathways identified in archaeal model organisms is given. Their energy gain is discussed, and new insights into different levels of regulation that have been observed so far, including the transcript and protein levels (e.g., gene regulation, known transcription regulators, and posttranslational modification via reversible protein phosphorylation), are presented.

Molecular basis of the unusual catalytic preference for GDP/GTP in Entamoeba histolytica 3-phosphoglycerate kinase

Phosphoglycerate kinase (EC 2.7.2.3) catalyzes reversible phosphoryl transfer from 1,3-bisphosphoglycerate to ADP to synthesize 3-phosphoglycerate and ATP during glycolysis. Phosphoglycerate kinases from several sources can use GDP/GTP as alternative substrates to ADP/ATP; however, the maximal velocities (V(m)) reached with the guanine nucleotides are approximately 50% of those displayed with the adenine nucleotides. By contrast, Entamoeba histolytica phosphoglycerate kinase (EC 2.7.2.10) is the only reported phosphoglycerate kinase displaying higher activity with GDP/GTP and lower affinities for the adenine nucleotides. To elucidate the molecular basis of the Entamoeba histolytica phosphoglycerate kinase selectivity for GDP/GTP, a conformational analysis was carried out on a homology model based on crystallographic structures of yeast and pig phosphoglycerate kinases. Some amino acid residues involved in the purine ring binding site not previously described were detected. Accordingly, Y239, E309 and V311 were replaced by site-directed mutagenesis in the Entamoeba histolytica phosphoglycerate kinase gene for the corresponding amino acid residues present in the adenine nucleotide-dependent phosphoglycerate kinases and the recombinant proteins were purified. Kinetic analysis of the enzymes showed that the single mutants Y239F, E309Q, E309M and V311L increased their catalytic efficiencies (V(m)/K(m)) with ADP/ATP as a result of both, increased V(m) and decreased K(m) values. Furthermore, a higher catalytic efficiency in the double mutant Y239F/E309M was achieved, which was mainly due to an increased affinity for ADP/ATP with a concomitant diminished affinity for GDP/GTP. The main Entamoeba histolytica phosphoglycerate kinase amino acid residues involved in the selectivity for guanine nucleotides were thus identified.

Molecular characterization of the 3-phosphoglycerate kinase gene (PGK1) from the methylotrophic yeast Pichia pastoris

We report the cloning of the 3-phosphoglycerate kinase gene (PGK1) from the methylotrophic yeast Pichia pastoris by a PCR approach. The coding sequence of the PGK1 gene comprises 1251 bp with the potential to encode a polypeptide of 416 amino acid residues, which shows high identity to homologous proteins from other yeasts. The promoter region of this gene (P(PGK1)) contains regulatory cis-elements found in other PGK1 genes, such as TATA box, CT-rich block and a heat shock element. In the 3' downstream region we identified a tripartite element 5'-TAG-TAGT-TTT-3', which is supposed to be important for transcription termination. As in other yeasts, the PGK1 gene from P. pastoris is present as a single-copy gene. Northern blot analysis revealed that the gene is transcribed as a 1.5 kb mRNA; when cells are grown on glucose the levels of this mRNA are increased two-fold in comparison to cells grown on glycerol. The transcriptional regulation of this gene by the carbon source was further confirmed when the alpha-amylase gene from Bacillus subtilis was placed under the control of P(PGK1): higher levels of expression were obtained when cells were grown on glucose as compared to glycerol and methanol. Preliminary results related to the strength of P(PGK1) show that it represents a potential alternative to constitutive heterologous expression in P. pastoris.

Distribution and phylogenies of enzymes of the Embden-Meyerhof-Parnas pathway from archaea and hyperthermophilic bacteria support a gluconeogenic origin of metabolism

Enzymes of the gluconeogenic/glycolytic pathway (the Embden-Meyerhof-Parnas (EMP) pathway), the reductive tricarboxylic acid cycle, the reductive pentose phosphate cycle and the Entner-Doudoroff pathway are widely distributed and are often considered to be central to the origins of metabolism. In particular, several enzymes of the lower portion of the EMP pathway (the so-called trunk pathway), including triosephosphate isomerase (TPI; EC 5.3.1.1), glyceraldehyde-3-phosphate dehydrogenase (GAPDH; EC 1.2.1.12/13), phosphoglycerate kinase (PGK; EC 2.7.2.3) and enolase (EC 4.2.1.11), are extremely well conserved and universally distributed among the three domains of life. In this paper, the distribution of enzymes of gluconeogenesis/glycolysis in hyperthermophiles--microorganisms that many believe represent the least evolved organisms on the planet--is reviewed. In addition, the phylogenies of the trunk pathway enzymes (TPIs, GAPDHs, PGKs and enolases) are examined. The enzymes catalyzing each of the six-carbon transformations in the upper portion of the EMP pathway, with the possible exception of aldolase, are all derived from multiple gene sequence families. In contrast, single sequence families can account for the archaeal and hyperthermophilic bacterial enzyme activities of the lower portion of the EMP pathway. The universal distribution of the trunk pathway enzymes, in combination with their phylogenies, supports the notion that the EMP pathway evolved in the direction of gluconeogenesis, i.e., from the bottom up.

Evolutionary relatedness between glycolytic enzymes most frequently occurring in genomes

More than 100 sequenced genomes were searched for genes coding for the enzymes involved in glycolysis in an effort to find the most frequently occurring ones. Triosephosphate isomerase (TIM), glyceraldehyde-3-phosphate dehydrogenase (GAPD), phosphoglycerate kinase (PGK) and enolase (ENOL) were found to be present in 90 investigated genomes all together. The final set consisted of 80 prokaryotic and 10 eukaryotic genomes. Of the 80 prokaryotic genomes, 73 were from Bacteria, 7 from Archaea. Two microbial genomes were also from Eucarya (yeasts). Eight genomes of nonmicrobial origin were included for comparison. The amino acid sequences of TIMs, GAPDs, PGKs and ENOLs were collected and aligned, and their individual as well as concatenated evolutionary trees were constructed and discussed. The trees clearly demonstrate a closer relatedness between Eucarya and Archaea (especially the concatenated tree) but they do not support the hypothesis that eukaryotic glycolytic enzymes should be closely related to their alpha-proteobacterial counterparts. Phylogenetic analyses further reveal that although the taxonomic groups (e.g., alpha-proteobacteria, gamma-proteobacteria, firmicutes, actinobacteria, etc.) form their more or less compact clusters in the trees, the inter-clade relationships between the trees are not conserved at all. On the other hand, several examples of conservative relatedness separating some clades of the same taxonomic groups were observed, e.g., Buchnera along with Wigglesworthia and the rest of gamma-proteobacteria, or mycoplasmas and the rest of firmicutes. The results support the view that these glycolytic enzymes may have their own evolutionary history.

Protein conformer selection by sequence-dependent packing contacts in crystals of 3-phosphoglycerate kinase

In several crystal structures of 3-phosphoglycerate kinase (PGK), the two domains occupy different relative positions. It is intriguing that the two extreme (open and closed) conformations have never been observed for the enzyme from the same species. Furthermore, in certain cases, these different crystalline conformations represent the enzyme-ligand complex of the same composition, such as the ternary complex containing either the substrate 3-phosphoglycerate (3-PG) and beta,gamma-imido-adenosine-5'-triphosphate (AMP-PNP), an analogue of the substrate MgATP, or 3-PG and the product MgADP. Thus, the protein conformation in the crystal is apparently determined by the origin of the isolated enzyme: PGK from pig muscle has only been crystallized in open conformation, whereas PGK from either Thermotoga maritima or Trypanosoma brucei has only been reported in closed conformations. A systematic analysis of the underlying sequence differences at the crucial hinge regions of the molecule and in the protein-protein contact surfaces in the crystal, in two independent pairs of open and closed states, have revealed that 1) sequential differences around the molecular hinges do not explain the appearance of fundamentally different conformations and 2) the species-specific intermolecular contacts between the nonconserved residues are responsible for stabilizing one conformation over the other in the crystalline state. A direct relationship between the steric position of the contacts in the three-dimensional structure and the conformational state of the protein has been demonstrated.

Factors influencing the stability of alpha-helices and beta-strands in thermophilic ribonuclease H

Understanding the influence of structural parameters is crucial to enhance the thermal stability of proteins. In this work, the stability (deltaG) of residues in different secondary structures of Ribonuclease H (RNase H) has been analyzed with 48 amino acid properties. The properties reflecting hydrophobicity show a good correlation with stability. Further, the linear distribution of surrounding hydrophobicity in alpha-helices, obtained from the three dimensional structure of thermophilic RNase H, agrees well with experimental deltaG values. Moreover, the stability parameters correlate better in alpha-helices than those did in beta-strand segments. Multiple regression analysis, incorporating combinations of three properties from among all possible combinations of the 48 properties, increased the correlation coefficient to 0.77.

A 1.8 A resolution structure of pig muscle 3-phosphoglycerate kinase with bound MgADP and 3-phosphoglycerate in open conformation: new insight into the role of the nucleotide in domain closure

3-phosphoglycerate kinase (PGK) is a typical kinase with two structural domains. The domains each bind one of the two substrates, 3-phosphoglycerate (3-PG) and MgATP. For the phospho-transfer reaction to take place the substrates must be brought closer by a hinge-bending domain closure. Open and closed structures of the enzyme with different relative domain positions have been determined from different species, but a comprehensive description of this conformational transition is yet to be attained. Crystals of pig muscle PGK in complex with MgADP and 3-phosphoglycerate were grown under the conditions which have previously resulted in crystals of the closed, catalytically competent conformation of Trypanosoma brucei PGK. The X-ray structure of the pig muscle ternary complex was determined at 1.8 A and the model was refined to R=20.8% and Rfree=24.1%. Contrary to expectation, however, it represents an essentially open conformation compared to that of T. brucei PGK. In addition, the beta-phosphate group of ADP is mobile in the new structure, in contrast to its well-defined position in T. brucei PGK. An extensive comparison of the ternary complexes from these remote species has been carried out in order to establish general differences between the two conformations and is reported here. A second pair of the open and closed structures was also compared. These analyses have made it possible to define several characteristic changes which accompany the structural transition, in addition to those identified previously: (1) the operation of a hinge at beta-strand L in the inter-domain region which greatly affects the relative domain positions; (2) the rearrangement and movement of helix 8, regulated through the interactions with the nucleotide phosphate; and (3) the existence of another hinge between helix 14 and the rest of the C-terminal part of the chain, which allows fine adjustment of the N-domain position. The main hinge at beta-strand L acts in concert with the C-terminal hinge at helix 7 described previously. Simultaneous interactions of the nucleotide phosphate groups with the loop that precedes helix 8, beta-strand J and the N terminus of helix 13 are required for propagation of the nucleotide effect towards the beta-strand L molecular hinge. A detailed description of the role of nucleotide binding in the hinge operation is presented.

Distributions of structural features contributing to thermostability in mesophilic and thermophilic alpha/beta barrel glycosyl hydrolases

Analysis of the structural basis for thermostability in proteins has come mainly from pairwise comparisons of mesophilic and thermophilic structures and has often yielded conflicting results. Interpretation of these results would be enhanced by knowing the normal range of features found for mesophilic proteins. In order to provide the average and distribution values of structural features among similar mesophilic proteins, we compared the amino acid composition, solvent accessible surface area, hydrogen bonds, number of ion pairs, and thermal factors of 22 structures of alpha/beta barrel glycosyl hydrolases. These distributions are then compared to values from seven alpha/beta barrel glycosyl hydrolases from thermophilic organisms. We find that the distribution of each structural feature is broad within the mesophilic proteins and illustrates the difficulty of making pairwise comparisons of mesophiles to thermophiles where differences for individual proteins may be within the normal range for the group. In comparing mesophiles to thermophiles as a group, we find that thermophilic structures have fewer glycines in a particular region of the structure and higher thermal factors at room temperature. These results suggest the basis for thermostability may be related to protein motion rather than to static features of protein structure.

Crystal structure of the glyceraldehyde-3-phosphate dehydrogenase from the hyperthermophilic archaeon Sulfolobus solfataricus

The enzyme glyceraldehyde-3-phosphate dehydrogenase (GAPDH) from the archaea shows low sequence identity (16-20%) with its eubacterial and eukaryotic counterparts. The crystal structure of the apo GAPDH from Sulfolobus solfataricus has been determined by multiple isomorphous replacement at 2.05 A resolution. The enzyme has several differences in secondary structure when compared with eubacterial GAPDHs, with an overall increase in the number of alpha-helices. There is a relocation of the active-site residues within the catalytic domain of the enzyme. The thermostability of the S. solfataricus enzyme can be attributed to a combination of an ion pair cluster and an intrasubunit disulphide bond.