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

Novel Structures of Type 1 Glyceraldehyde-3-phosphate Dehydrogenase from Escherichia coli Provide New Insights into the Mechanism of Generation of 1,3-Bisphosphoglyceric Acid

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) is a highly conserved enzyme involved in the ubiquitous process of glycolysis and presents a loop (residues 208-215 of Escherichia coli GAPDH) in two alternative conformations (I and II). It is uncertain what triggers this loop rearrangement, as well as which is the precise site from which phosphate attacks the thioacyl intermediate precursor of 1,3-bisphosphoglycerate (BPG). To clarify these uncertainties, we determined the crystal structures of complexes of wild-type GAPDH (WT) with NAD and phosphate or G3P, and of essentially inactive GAPDH mutants (C150S, H177A), trapping crystal structures for the thioacyl intermediate or for ternary complexes with NAD and either phosphate, BPG, or G3P. Analysis of these structures reported here lead us to propose that phosphate is located in the "new Pi site" attacks the thioester bond of the thioacyl intermediate to generate 1,3-bisphosphoglyceric acid (BPG). In the structure of the thioacyl intermediate, the mobile loop is in conformation II in subunits O, P, and R, while both conformations coexist in subunit Q. Moreover, only the Q subunit hosts bound NADH. In the R subunit, only the pyrophosphate part of NADH is well defined, and NADH is totally absent from the O and P subunits. Thus, the change in loop conformation appears to occur after NADH is produced, before NADH is released. In addition, two new D-glyceraldehyde-3-phosphate (G3P) binding forms are observed in WT.NAD.G3P and C150A+H177A.NAD.G3P. In summary, this paper improves our understanding of the GAPDH catalytic mechanism, particularly regarding BPG formation.

Systems biology of the modified branched Entner-Doudoroff pathway in Sulfolobus solfataricus

Sulfolobus solfataricus is a thermoacidophilic Archaeon that thrives in terrestrial hot springs (solfatares) with optimal growth at 80°C and pH 2–4. It catabolizes specific carbon sources, such as D-glucose, to pyruvate via the modified Entner-Doudoroff (ED) pathway. This pathway has two parallel branches, the semi-phosphorylative and the non-phosphorylative. However, the strategy of S.solfataricus to endure in such an extreme environment in terms of robustness and adaptation is not yet completely understood. Here, we present the first dynamic mathematical model of the ED pathway parameterized with quantitative experimental data. These data consist of enzyme activities of the branched pathway at 70°C and 80°C and of metabolomics data at the same temperatures for the wild type and for a metabolic engineered knockout of the semi-phosphorylative branch. We use the validated model to address two questions: 1. Is this system more robust to perturbations at its optimal growth temperature? 2. Is the ED robust to deletion and perturbations? We employed a systems biology approach to answer these questions and to gain further knowledge on the emergent properties of this biological system. Specifically, we applied deterministic and stochastic approaches to study the sensitivity and robustness of the system, respectively. The mathematical model we present here, shows that: 1. Steady state metabolite concentrations of the ED pathway are consistently more robust to stochastic internal perturbations at 80°C than at 70°C; 2. These metabolite concentrations are highly robust when faced with the knockout of either branch. Connected with this observation, these two branches show different properties at the level of metabolite production and flux control. These new results reveal how enzyme kinetics and metabolomics synergizes with mathematical modelling to unveil new systemic properties of the ED pathway in S.solfataricus in terms of its adaptation and robustness.

Improving the 'tool box' for robust industrial enzymes

The speed of sequencing of microbial genomes and metagenomes is providing an ever increasing resource for the identification of new robust biocatalysts with industrial applications for many different aspects of industrial biotechnology. Using 'natures catalysts' provides a sustainable approach to chemical synthesis of fine chemicals, general chemicals such as surfactants and new consumer-based materials such as biodegradable plastics. This provides a sustainable and 'green chemistry' route to chemical synthesis which generates no toxic waste and is environmentally friendly. In addition, enzymes can play important roles in other applications such as carbon dioxide capture, breakdown of food and other waste streams to provide a route to the concept of a 'circular economy' where nothing is wasted. The use of improved bioinformatic approaches and the development of new rapid enzyme activity screening methodology can provide an endless resource for new robust industrial biocatalysts.This mini-review will discuss several recent case studies where industrial enzymes of 'high priority' have been identified and characterised. It will highlight specific hydrolase enzymes and recent case studies which have been carried out within our group in Exeter.

D-Glyceraldehyde-3-Phosphate Dehydrogenase Structure and Function

Aside from its well-established role in glycolysis, glyceraldehyde-3-phosphate dehydrogenase (GAPDH) has been shown to possess many key functions in cells. These functions are regulated by protein oligomerization , biomolecular interactions, post-translational modifications , and variations in subcellular localization . Several GAPDH functions and regulatory mechanisms overlap with one another and converge around its role in intermediary metabolism. Several structural determinants of the protein dictate its function and regulation. GAPDH is ubiquitously expressed and is found in all domains of life. GAPDH has been implicated in many diseases, including those of pathogenic, cardiovascular, degenerative, diabetic, and tumorigenic origins. Understanding the mechanisms by which GAPDH can switch between its functions and how these functions are regulated can provide insights into ways the protein can be modulated for therapeutic outcomes.

Archaeal Enzymes and Applications in Industrial Biocatalysts

Archaeal enzymes are playing an important role in industrial biotechnology. Many representatives of organisms living in “extreme” conditions, the so-called Extremophiles, belong to the archaeal kingdom of life. This paper will review studies carried by the Exeter group and others regarding archaeal enzymes that have important applications in commercial biocatalysis. Some of these biocatalysts are already being used in large scale industrial processes for the production of optically pure drug intermediates and amino acids and their analogues. Other enzymes have been characterised at laboratory scale regarding their substrate specificity and properties for potential industrial application. The increasing availability of DNA sequences from new archaeal species and metagenomes will provide a continuing resource to identify new enzymes of commercial interest using both bioinformatics and screening approaches.

Structural and catalytic effects of proline substitution and surface loop deletion in the extended active site of human carbonic anhydrase II

Bioengineering of a thermophilic enzyme starting from a mesophilic scaffold has proven to be a significant challenge, as several stabilizing elements have been proposed to be the foundation of thermal stability, including disulfide bridges, surface loop reduction, ionic pair networks, proline substitutions and aromatic clusters. This study emphasizes the effect of increasing the rigidity of human carbonic anhydrase II (HCA II; EC 4.2.1.1) via incorporation of proline residues at positions 170 and 234, which are located in surface loops that are able to accommodate restrictive main-chain conformations without rearrangement of the surrounding peptide backbone. Additionally, the effect of the compactness of HCA II was examined by deletion of a surface loop (residues 230-240) that had been previously identified as a possible source of thermal stability for the hyperthermophilic carbonic anhydrase isolated from the bacterium Sulfurihydrogenibium yellowstonense YO3AOP1. Differential scanning calorimetry analysis of these HCA II variants revealed that these structural modifications had a minimum effect on the thermal stability of the enzyme, while kinetic studies showed unexpected effects on the catalytic efficiency and proton transfer rates. X-ray crystallographic analysis of these HCA II variants showed that the electrostatic potential and configuration of the highly acidic loop (residues 230-240) play an important role in its high catalytic activity. Based on these observations and previous studies, a picture is emerging of the various components within the general structural architecture of HCA II that are key to stability. These elements may provide blueprints for rational thermal stability engineering of other enzymes.

DATABASE

Structural data have been submitted to the Protein Data Bank under accession numbers 4QK1 (K170P), 4QK2 (E234P) and 4QK3 (Δ230-240).

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.

Dynamic oligomeric properties

This chapter provides a foundation for further research into the relationship between dynamic oligomeric properties and functional diversity. The structural basis that underlies the conformational sub-states of the GAPDH oligomer is discussed. The issue of protein stability is given a thorough analysis, since it is well-established that the primary strategy for protein oligomerization is to stabilize conformation. Several factors that affect oligomerization are described, including chemical modification by synthetic reagents. The effects of native substrates and coenzymes are also discussed. The curious feature of chloride ions having a de-stabilizing effect on native GAPDH structure is described. Additionally, the role of adenine dinucleotides in tetramer-dimer equilibrium dynamics is suggested to be a major part of the physiological regulation of GAPDH structure and function. This chapter also contends that a vast amount of useful information can come from comparative analyses of diverse species, particularly regarding protein stability and subunit-subunit interaction. Lastly, the concept of domain exchange is introduced as a means of understanding the stabilization of dynamic oligomers, suggesting that inter-subunit contacts may also be a way of masking docking sites to other proteins.

A pursuit of lineage-specific and niche-specific proteome features in the world of archaea

Background

Archaea evoke interest among researchers for two enigmatic characteristics –a combination of bacterial and eukaryotic components in their molecular architectures and an enormous diversity in their life-style and metabolic capabilities. Despite considerable research efforts, lineage- specific/niche-specific molecular features of the whole archaeal world are yet to be fully unveiled. The study offers the first large-scale in silico proteome analysis of all archaeal species of known genome sequences with a special emphasis on methanogenic and sulphur-metabolising archaea.

Results

Overall amino acid usage in archaea is dominated by GC-bias. But the environmental factors like oxygen requirement or thermal adaptation seem to play important roles in selection of residues with no GC-bias at the codon level. All methanogens, irrespective of their thermal/salt adaptation, show higher usage of Cys and have relatively acidic proteomes, while the proteomes of sulphur-metabolisers have higher aromaticity and more positive charges. Despite of exhibiting thermophilic life-style, korarchaeota possesses an acidic proteome. Among the distinct trends prevailing in COGs (Cluster of Orthologous Groups of proteins) distribution profiles, crenarchaeal organisms display higher intra-order variations in COGs repertoire, especially in the metabolic ones, as compared to euryarchaea. All methanogens are characterised by a presence of 22 exclusive COGs.

Conclusions

Divergences in amino acid usage, aromaticity/charge profiles and COG repertoire among methanogens and sulphur-metabolisers, aerobic and anaerobic archaea or korarchaeota and nanoarchaeota, as elucidated in the present study, point towards the presence of distinct molecular strategies for niche specialization in the archaeal world.

Biochemical characterization of glyceraldehyde-3-phosphate dehydrogenase from Thermococcus kodakarensis KOD1

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) plays an essential role in glycolysis by catalyzing the conversion of D-glyceraldehyde 3-phosphate (D-G3P) to 1,3-diphosphoglycerate using NAD(+) as a cofactor. In this report, the GAPDH gene from the hyperthermophilic archaeon Thermococcus kodakarensis KOD1 (GAPDH-tk) was cloned and the protein was purified to homogeneity. GAPDH-tk exists as a homotetramer with a native molecular mass of 145 kDa; the subunit molecular mass was 37 kDa. GAPDH-tk is a thermostable protein with a half-life of 5 h at 80-90°C. The apparent K (m) values for NAD(+) and D-G3P were 77.8 ± 7.5 μM and 49.3 ± 3.0 μM, respectively, with V (max) values of 45.1 ± 0.8 U/mg and 59.6 ± 1.3 U/mg, respectively. Transmission electron microscopy (TEM) and image processing confirmed that GAPDH-tk has a tetrameric structure. Interestingly, GAPDH-tk migrates as high molecular mass forms (~232 kDa and ~669 kDa) in response to oxidative stress.

Structure of glyceraldehyde-3-phosphate dehydrogenase from the archaeal hyperthermophile Methanocaldococcus jannaschii

The X-ray crystal structure of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) from the hyperthermophilic archaeon Methanocaldococcus jannaschii (Mj-GAPDH) was determined to 1.81 A resolution. The crystal belonged to space group C222(1), with unit-cell parameters a = 83.4, b = 152.0, c = 118.6 A. The structure was solved by molecular replacement and was refined to a final R factor of 17.1% (R(free) = 19.8%). The final structure included the cofactor NADP(+) at the nucleotide-binding site and featured unoccupied inorganic and substrate phosphate-binding sites. A comparison with GAPDH structures from mesophilic sources suggested that Mj-GAPDH is stabilized by extensive electrostatic interactions between the C-terminal alpha-helices and various distal loop regions, which are likely to contribute to thermal stability. The key phosphate-binding residues in the active site of Mj-GAPDH are conserved in other archaeal GAPDH proteins. These residues undergo a conformational shift in response to occupancy of the inorganic phosphate site.

Structure of the extended diarrhea-inducing domain of rotavirus enterotoxigenic protein NSP4

Rotavirus nonstructural protein 4 (NSP4) is a multidomainal and multifunctional protein and is recognized as the first virus-encoded enterotoxin. Extensive efforts to crystallize the complete cytoplasmic tail (CT), which exhibits all the known biological functions, have been unsuccessful, and to date, the structure of only a synthetic peptide corresponding to amino acids (aa) 95-137 has been reported. Recent studies indicate that the interspecies-variable domain (ISVD) from aa 135 to 141 as well as the extreme C-terminus are critical determinants of virus virulence and the diarrhea-inducing ability of the protein. Among the five NSP4 genotypes identified, those belonging to genotypes A1, B and C possess either a proline at position 138 or a glycine at 140, while those of A2, D and E lack these residues in the ISVD, suggesting conformational differences in this region among different NSP4s. Here, we examined the crystallization properties of several deletion mutants and report the structure of a recombinant mutant, NSP4:95-146, lacking the N-terminal 94 and C-terminal 29 aa, from SA11 (A1) and I321 (A2) at 1.67 and 2.7 A, respectively. In spite of the high resolution of one of the structures, electron density for the C-terminal 9 residues could not be seen for either of the mutants, and the crystal packing resulted in the creation of a clear empty space for this region. Extension of the unstructured C-terminus beyond aa 146 hindered crystallization under the experimental conditions. The present structure revealed significant differences from that of the synthetic peptide in the conformation of amino acids at the end of the helix as well as the crystal packing owing to the additional space required to accommodate the un structured virulence-determining region. The crystal structure and secondary structure prediction of the NSP4:95-146 mutants from different genotypes suggest that the region C-terminal to aa 137 in all the NSP4 proteins is likely to be unstructured, and this might be of structural and biological functional significance.

Unusual pathways and enzymes of central carbohydrate metabolism in Archaea

Sugar-utilizing hyperthermophilic and halophilic Archaea degrade glucose and glucose polymers to acetate or to CO2 using O2, nitrate, sulfur or sulfate as electron acceptors. Comparative analyses of glycolytic pathways in these organisms indicate a variety of differences from the classical Emden-Meyerhof and Entner-Doudoroff pathways that are operative in Bacteria and Eukarya, respectively. The archaeal pathways are characterized by the presence of numerous novel enzymes and enzyme families that catalyze, for example, the phosphorylation of glucose and of fructose 6-phosphate, the isomerization of glucose 6-phosphate, the cleavage of fructose 1,6-bisphosphate, the oxidation of glyceraldehyde 3-phosphate and the conversion of acetyl-CoA to acetate. Recent major advances in deciphering the complexity of archaeal central carbohydrate metabolism were gained by combination of classical biochemical and genomic-based approaches.

A new branch in the family: structure of aspartate-beta-semialdehyde dehydrogenase from Methanococcus jannaschii

The structure of aspartate-beta-semialdehyde dehydrogenase (ASADH) from Methanococcus jannaschii has been determined to 2.3 angstroms resolution using multiwavelength anomalous diffraction (MAD) phasing of a selenomethionine-substituted derivative to define a new branch in the family of ASADHs. This new structure has a similar overall fold and domain organization despite less than 10% conserved sequence identity with the bacterial enzymes. However, the entire repertoire of functionally important active site amino acid residues is conserved, suggesting an identical catalytic mechanism but with lower catalytic efficiency. A new coenzyme-binding conformation and dual NAD/NADP coenzyme specificity further distinguish this archaeal branch from the bacterial ASADHs. Several structural differences are proposed to account for the dramatically enhanced thermostability of this archaeal enzyme. Finally, the intersubunit communication channel connecting the active sites in the bacterial enzyme dimer has been disrupted in the archaeal ASADHs by amino acid changes that likely prevent the alternating sites reactivity previously proposed for the bacterial ASADHs.

Posttranslational protein modification in Archaea

One of the first hurdles to be negotiated in the postgenomic era involves the description of the entire protein content of the cell, the proteome. Such efforts are presently complicated by the various posttranslational modifications that proteins can experience, including glycosylation, lipid attachment, phosphorylation, methylation, disulfide bond formation, and proteolytic cleavage. Whereas these and other posttranslational protein modifications have been well characterized in Eucarya and Bacteria, posttranslational modification in Archaea has received far less attention. Although archaeal proteins can undergo posttranslational modifications reminiscent of what their eucaryal and bacterial counterparts experience, examination of archaeal posttranslational modification often reveals aspects not previously observed in the other two domains of life. In some cases, posttranslational modification allows a protein to survive the extreme conditions often encountered by Archaea. The various posttranslational modifications experienced by archaeal proteins, the molecular steps leading to these modifications, and the role played by posttranslational modification in Archaea form the focus of this review.

The structure and function of Xenopus NO38-core, a histone chaperone in the nucleolus

Xenopus NO38 is an abundant nucleolar chaperone and a member of the nucleoplasmin (Np) family. Here, we report high-resolution crystal structures of the N-terminal domain of NO38, as a pentamer and a decamer. As expected, NO38 shares the Np family fold. In addition, NO38- and Np-core pentamers each use highly conserved residues and numerous waters to form their respective decamers. Further studies show that NO38 and Np each bind equal amounts of the four core histones. However, NO38 prefers the (H3-H4)(2) tetramer, while Np probably prefers H2A-H2B dimers. We also show that NO38 and Np will each bind noncognate histones when the preferred partner is absent. We suggest that these chaperones must form decamers in order to bind histones and differentiate between histone tetramers and dimers. When taken together, these data imply that NO38 may function as a histone chaperone in the nucleolus.

Isocitrate dehydrogenase from the hyperthermophile Aeropyrum pernix: X-ray structure analysis of a ternary enzyme-substrate complex and thermal stability

Isocitrate dehydrogenase from Aeropyrum pernix (ApIDH) is a homodimeric enzyme that belongs to the beta-decarboxylating dehydrogenase family and is the most thermostable IDH identified. It catalyzes the NADP+ and metal-dependent oxidative decarboxylation of isocitrate to alpha-ketoglutarate. We have solved the crystal structures of a native ApIDH at 2.2 A, a pseudo-native ApIDH at 2.1 A, and of ApIDH in complex with NADP+, Ca2+ and d-isocitrate at 2.3 A. The pseudo-native ApIDH is in complex with etheno-NADP+ which was located at the surface instead of in the active site revealing a novel adenine-nucleotide binding site in ApIDH. The native and the pseudo-native ApIDHs were found in an open conformation, whereas one of the subunits of the ternary complex was closed upon substrate binding. The closed subunit showed a domain rotation of 19 degrees compared to the open subunit. The binding of isocitrate in the closed subunit was identical with that of the binary complex of porcine mitochondrial IDH, whereas the binding of NADP+ was similar to that of the ternary complex of IDH from Escherichiacoli. The reaction mechanism is likely to be conserved in the different IDHs. A proton relay chain involving at least five solvent molecules, the 5'-phosphate group of the nicotinamide-ribose and a coupled lysine-tyrosine pair in the active site, is postulated as essential in both the initial and the final steps of the catalytic reaction of IDH. ApIDH was found to be highly homologous to the mesophilic IDHs and was subjected to a comparative analysis in order to find differences that could explain the large difference in thermostability. Mutational studies revealed that a disulfide bond at the N terminus and a seven-membered inter-domain ionic network at the surface are major determinants for the higher thermostability of ApIDH compared to EcIDH. Furthermore, the total number of ion pairs was dramatically higher in ApIDH compared to the mesophilic IDHs if a cutoff of 4.2 A was used. A calculated net charge of only +1 compared to -19 and -25 in EcIDH and BsIDH, respectively, suggested a high degree of electrostatic optimization, which is known to be an important determinant for increased thermostability.

Effects of a novel disulfide bond and engineered electrostatic interactions on the thermostability of azurin

Identification and evaluation of factors important for thermostability in proteins is a growing research field with many industrial applications. This study investigates the effects of introducing a novel disulfide bond and engineered electrostatic interactions with respect to the thermostability of holo azurin from Pseudomonas aeruginosa. Four mutants were selected on the basis of rational design and novel temperature-dependent atomic displacement factors from crystal data collected at elevated temperatures. The atomic displacement parameters describe the molecular movement at higher temperatures. The thermostability was evaluated by optical spectroscopy as well as by differential scanning calorimetry. Although azurin has a high inherent stability, the introduction of a novel disulfide bond connecting a flexible loop with small alpha-helix (D62C/K74C copper-containing mutant), increased the T(m) by 3.7 degrees C compared with the holo protein. Furthermore, three mutants were designed to introduce electrostatic interactions, K24R, D23E/K128R, and D23E/K128R/K24R. Mutant K24R stabilizes loops between two separate beta-strands and D23E/K128R was selected to stabilize the C-terminus of azurin. Furthermore, D23E/K128R/K24R was selected to reflect the combination of the electrostatic interactions in D23E/K128R and K24R. The mutants involving electrostatic interactions had a minor effect on the thermostability. The crystal structures of the copper-containing mutants D62C/K74C and K24R have been determined to 1.5 and 1.8 A resolution. In addition the crystal structure of the zinc-loaded mutant D62C/K74C has also been completed to 1.8 A resolution. These structures support the selected design and provide valuable information for evaluating effects of the modifications on the thermostability of holo azurin.

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.

Structure and Function of a Regulated Archaeal Triosephosphate Isomerase Adapted to High Temperature

Triosephophate isomerase (TIM) is a dimeric enzyme in eucarya, bacteria and mesophilic archaea. In hyperthermophilic archaea, however, TIM exists as a tetramer composed of monomers that are about 10% shorter than other eucaryal and bacterial TIM monomers. We report here the crystal structure of TIM from Thermoproteus tenax, a hyperthermophilic archaeon that has an optimum growth temperature of 86 degrees C. The structure was determined from both a hexagonal and an orthorhombic crystal form to resolutions of 2.5A and 2.3A, and refined to R-factors of 19.7% and 21.5%, respectively. In both crystal forms, T.tenax TIM exists as a tetramer of the familiar (betaalpha)(8)-barrel. In solution, however, and unlike other hyperthermophilic TIMs, the T.tenax enzyme exhibits an equilibrium between inactive dimers and active tetramers, which is shifted to the tetramer state through a specific interaction with glycerol-1-phosphate dehydrogenase of T.tenax. This observation is interpreted in physiological terms as a need to reduce the build-up of thermolabile metabolic intermediates that would be susceptible to destruction by heat. A detailed structural comparison with TIMs from organisms with growth optima ranging from 15 degrees C to 100 degrees C emphasizes the importance in hyperthermophilic proteins of the specific location of ionic interactions for thermal stability rather than their numbers, and shows a clear correlation between the reduction of heat-labile, surface-exposed Asn and Gln residues with thermoadaptation. The comparison confirms the increase in charged surface-exposed residues at the expense of polar residues.

The crystal structure of the tryptophan synthase beta subunit from the hyperthermophile Pyrococcus furiosus. Investigation of stabilization factors

The structure of the tryptophan synthase beta2 subunit (Pfbeta2) from the hyperthermophile, Pyrococcus furiosus, was determined by X-ray crystallographic analysis at 2.2 A resolution, and its stability was examined by DSC. This is the first report of the X-ray structure of the tryptophan synthase beta2 subunit alone, although the structure of the tryptophan synthase alpha2beta2 complex from Salmonella typhimurium has already been reported. The structure of Pfbeta2 was essentially similar to that of the beta2 subunit (Stbeta2) in the alpha2beta2 complex from S. typhimurium. The sequence alignment with secondary structures of Pfbeta and Stbeta in monomeric form showed that six residues in the N-terminal region and three residues in the C-terminal region were deleted in Pfbeta, and one residue at Pro366 of Stbeta and at Ile63 of Pfbeta was inserted. The denaturation temperature of Pfbeta2 was higher by 35 degrees C than the reported values from mesophiles at approximately pH 8. On the basis of structural information on both proteins, the analyses of the contributions of each stabilization factor indicate that: (a) the higher stability of Pfbeta2 is not caused by either a hydrophobic interaction or an increase in ion pairs; (b) the number of hydrogen bonds involved in the main chains of Pfbeta is greater by about 10% than that of Stbeta, indicating that the secondary structures of Pfbeta are more stabilized than those of Stbeta and (c) the sequence of Pfbeta seems to be better fitted to an ideally stable structure than that of Stbeta, as assessed from X-ray structure data.

Template-assisted filament growth by parallel stacking of tau

Tau filaments are found in >20 neurodegenerative diseases. Yet, because of their enormous molecular weights and poor tendency to form highly ordered 3D crystal lattices, they have evaded high-resolution structure determination. Here, we studied 25 derivatized tau mutants by using electron paramagnetic resonance and fluorescence spectroscopy to report structural details of tau filaments. Based on strong spin exchange and pyrene excimer formation of core residues, we find that individual tau proteins form single molecule layers along the fiber axis that perfectly stack on top of each other by in-register, parallel alignment of beta-strands. We suggest a model of filament growth wherein the existing filament serves as a template for the incoming, unfolded tau molecule, resulting in a new structured layer with maximized hydrogen-bonded contact surface and side-chain stacking.

An electrostatic basis for the stability of thermophilic proteins

Two factors provide key contributions to the stability of thermophilic proteins relative to their mesophilic homologues: electrostatic interactions of charged residues in the folded state and the dielectric response of the folded protein. The dielectric response for proteins in a "thermophilic series" globally modulates the thermal stability of its members, with the calculated dielectric constant for the protein increasing from mesophiles to hyperthermophiles. This variability results from differences in the distribution of charged residues on the surface of the protein, in agreement with structural and genetic observations. Furthermore, the contribution of electrostatic interactions to the stability of the folded state is more favorable for thermophilic proteins than for their mesophilic homologues. This leads to the conclusion that electrostatic interactions play an important role in determining the stability of proteins at high temperatures. The interplay between electrostatic interactions and dielectric response also provides further rationalization for the enhanced stability of thermophilic proteins with respect to cold-denaturation. Taken together, the distribution of charged residues and their fluctuations have been shown to be factors in modulating protein stability over the entire range of biologically relevant temperatures.

The structure of an alcohol dehydrogenase from the hyperthermophilic archaeon Aeropyrum pernix

The structure of the recombinant medium chain alcohol dehydrogenase (ADH) from the hyperthermophilic archaeon Aeropyrum pernix has been solved by the multiple anomalous dispersion technique using the signal from the naturally occurring zinc ions. The enzyme is a tetramer with 222 point group symmetry. The ADH monomer is formed from a catalytic and a cofactor-binding domain, with the overall fold similar to previously solved ADH structures. The 1.62 A resolution A.pernix ADH structure is that of the holo form, with the cofactor NADH bound into the cleft between the two domains. The electron density found in the active site has been interpreted to be octanoic acid, which has been shown to be an inhibitor of the enzyme. This inhibitor is positioned with its carbonyl oxygen atom forming the fourth ligand of the catalytic zinc ion. The structural zinc ion of each monomer is present at only partial occupancy and in its absence a disulfide bond is formed. The enhanced thermal stability of the A.pernix ADH is thought to arise primarily from increased ionic and hydrophobic interactions on the subunit interfaces.

Metal and redox modulation of cysteine protein function

In biological systems, the amino acid cysteine combines catalytic activity with an extensive redox chemistry and unique metal binding properties. The interdependency of these three aspects of the thiol group permits the redox regulation of proteins and metal binding, metal control of redox activity, and ligand control of metal-based enzyme catalysis. Cysteine proteins are therefore able to act as "redox switches," to sense concentrations of oxidative stressors and unbound zinc ions in the cytosol, to provide a "storage facility" for excess metal ions, to control the activity of metalloproteins, and to take part in important regulatory and signaling pathways. The diversity of cysteine's multiple roles in vivo is equally as fascinating as it is promising for future biochemical and pharmacological research.

Dual coenzyme specificity of photosynthetic glyceraldehyde-3-phosphate dehydrogenase interpreted by the crystal structure of A4 isoform complexed with NAD

Photosynthetic glyceraldehyde-3-phosphate dehydrogenase (GAPDH) of Spinacia oleracea belongs to a wide group of GAPDHs found in most organisms displaying oxygenic photosynthesis, including cyanobacteria, green and red algae, and higher plants. As a major catalytic difference with respect to glycolytic GAPDH, photosynthetic GAPDH exhibits dual cofactor specificity toward pyridine nucleotides with a preference for NADP(H). Here we report the crystal structure of NAD-complexed recombinant A(4)-GAPDH (NAD-A(4)-GAPDH) from Spinacia oleracea, expressed in Escherichia coli. Its superimposition onto native A(4)-GAPDH complexed with NADP (NADP-A(4)-GAPDH) pinpoints specific conformational changes resulting from cofactor replacement. In photosynthetic NAD-A(4)-GAPDH, the side chain of Asp32 is oriented toward the coenzyme to interact with the adenine ribose diol, similar to glycolytic GAPDHs (NAD-specific). On the contrary, in NADP-A(4)-GAPDH Asp32 moves away to accommodate the additional 2'-phosphate group of the coenzyme and to minimize electrostatic repulsion. Asp32 rotation is allowed by the presence of the small residue Ala40, conserved in most photosynthetic GAPDHs, replacing bulky amino acid side chains in glycolytic GAPDHs. While in NADP-A(4)-GAPDH two amino acids, Thr33 and Ser188, are involved in hydrogen bonds with the 2'-phosphate group of NADP, in the NAD-complexed enzyme these interactions are lacking. The crystallographic structure of NAD-A(4)-GAPDH highlights that four residues, Thr33, Ala40, Ser188, and Ala187 (Leu, Leu, Pro, and Leu respectively, in glycolytic Bacillus stearothermophilus GAPDH sequence) are of primary importance for the dual cofactor specificity of photosynthetic GAPDH. These modifications seem to trace the minimum evolutionary route for a primitive NAD-specific GAPDH to be converted into the NADP-preferring enzyme of oxygenic photosynthetic organisms.

The crystal structure of Drosophila NLP-core provides insight into pentamer formation and histone binding

The nucleoplasmin-like protein from Drosophila (dNLP) functions as a chaperone for core histones and may remodel chromatin in embryos. We now report the crystal structure of a dNLP-core pentamer at 1.5 A resolution. The monomer has an eight-stranded, beta barrel topology that is similar to nucleoplasmin (Np). However, a signature beta hairpin is tucked in along the lateral surface of the dNLP-core pentamer, while it extends outward in the Np-core decamer. Drosophila NLP and Np both assemble histone octamers. This process may require each chaperone to form a decamer, which would create symmetric binding sites for the histones. Conformational differences between dNLP and Np may reflect their different oligomeric states, while a conserved, nonpolar subunit interface may allow conformational plasticity during histone binding.

Stepwise adaptations of citrate synthase to survival at life's extremes. From psychrophile to hyperthermophile

The crystal structure of citrate synthase from the thermophilic Archaeon Sulfolobus solfataricus (optimum growth temperature = 85 degrees C) has been determined, extending the number of crystal structures of citrate synthase from different organisms to a total of five that span the temperature range over which life exists (from psychrophile to hyperthermophile). Detailed structural analysis has revealed possible molecular mechanisms that determine the different stabilities of the five proteins. The key to these mechanisms is the precise structural location of the additional interactions. As one ascends the temperature ladder, the subunit interface of this dimeric enzyme and loop regions are reinforced by complex electrostatic interactions, and there is a reduced exposure of hydrophobic surface. These observations reveal a progressive pattern of stabilization through multiple additional interactions at solvent exposed, loop and interfacial regions.

Comparative analysis of thermoadaptation within the archaeal glyceraldehyde-3-phosphate dehydrogenases from mesophilic Methanobacterium bryantii and thermophilic Methanothermus fervidus

To gain insight into the molecular determinants of thermoadaptation within the family of archaeal glyceraldehyde-3-phosphate dehydrogenases (GAPDH), a homology-based 3-D model of the mesophilic GAPDH from Methanobacterium bryantii was built and compared with the crystal structure of the thermophilic GAPDH from Methanothermus fervidus. The homotetrameric model of the holoenzyme was initially assembled from identical subunits completed with NADP molecules. The structure was then refined by energy minimization and simulated-annealing procedures. PROCHECK and the 3-D profile method were used to appraise the model reliability. Striking molecular features underlying the difference in stability between the enzymes were deduced from their structural comparison. First, both the increase in hydrophobic contacts and the decrease in accessibility to the protein core were shown to discriminate in favor of the thermophilic enzyme. Besides, but to a lesser degree, the number of ion pairs involved in cooperative clusters appeared to correlate with thermostability. Finally, the decreased stability of the mesophilic enzyme was also predicted to proceed from both the lack of charge-dipole interactions within alpha-helices and the enhanced entropy of unfolding due to an increase in chain flexibility. Thus, archaeal GAPDHs appear to be governed by thermoadaptation rules that differ in some aspects from those previously observed within their eubacterial counterparts.

Ionic network at the C-terminus of the beta-glycosidase from the hyperthermophilic archaeon Sulfolobus solfataricus: Functional role in the quaternary structure thermal stabilization

Biochemical, crystallographic, and computational data support the hypothesis that electrostatic interactions are among the dominant forces in stabilizing hyperthermophilic proteins. The thermostable beta-glycosidase from the hyperthermophile Sulfolobus solfataricus (Ssbeta-gly) is an interesting model system for the study of protein adaptation to high temperatures. The largest ion-pair network of Ssbeta-gly is located at the tetrameric interface of the molecule; in this paper, key residues in this region were modified by site-directed mutagenesis and the stability of the mutants was analyzed by kinetics of thermal denaturation. All mutations produced faster enzyme inactivation, suggesting that the C-terminal ionic network prevents the dissociation into monomers, which is the limiting step in the mechanism of Ssbeta-gly inactivation. Moreover, the calculated reaction order showed that the mechanism of inactivation depends on the mutation introduced, suggesting that intermediates maintaining enzymatic activity are produced during the inactivation transition of some, but not all, mutants. Molecular models of each mutant allow us to rationalize the experimental evidence and give support to the current theories on the mechanism of ion pair stabilization in proteins from hyperthermophiles.

Crystal structure of the non-regulatory A(4 )isoform of spinach chloroplast glyceraldehyde-3-phosphate dehydrogenase complexed with NADP

Here, we report the first crystal structure of a photosynthetic glyceraldehyde-3-phosphate dehydrogenase (GAPDH) complexed with NADP. The enzyme, purified from spinach chloroplasts, is constituted of a single type of subunit (A) arranged in homotetramers. It shows non-regulated NADP-dependent and NAD-dependent activities, with a preference for NADP. The structure has been solved to 3.0 A resolution by molecular replacement. The crystals belong to space group C222 with three monomers in the asymmetric unit. One of the three monomers generates a tetramer using the space group 222 point symmetry and a very similar tetramer is generated by the other two monomers, related by a non-crystallographic symmetry, using a crystallographic 2-fold axis. The protein reveals a large structural homology with known GAPDHs both in the cofactor-binding domain and in regions of the catalytic domain. Like all other GAPDHs investigated so far, the A(4)-GAPDH belongs to the Rossmann fold family of dehydrogenases. However, unlike most dehydrogenases of this family, the adenosine 2'-phosphate group of NADP does not form a salt-bridge with any positively charged residue in its surroundings, being instead set in place by hydrogen bonds with a threonine residue belonging to the Rossmann fold and a serine residue located in the S-loop of a symmetry-related monomer. While increasing our knowledge of an important photosynthetic enzyme, these results contribute to a general understanding of NADP versus NAD recognition in pyridine nucleotide-dependent enzymes. Although the overall structure of A(4)-GAPDH is similar to that of the cytosolic GAPDH from bacteria and eukaryotes, the chloroplast tetramer is peculiar, in that it can actually be considered a dimer of dimers, since monomers are bound in pairs by a disulphide bridge formed across Cys200 residues. This bridge is not found in other cytosolic or chloroplast GAPDHs from animals, bacteria, or plants other than spinach.

The architecture of parallel beta-helices and related folds

Three-dimensional structures have been determined of a large number of proteins characterized by a repetitive fold where each of the repeats (coils) supplies a strand to one or more parallel beta-sheets. Some of these proteins form superfamilies of proteins, which have probably arisen by divergent evolution from a common ancestor. The classical example is the family including four families of pectinases without obviously related primary sequences, the phage P22 tailspike endorhamnosidase, chrondroitinase B and possibly pertactin from Bordetella pertusis. These show extensive stacking of similar residues to give aliphatic, aromatic and polar stacks such as the asparagine ladder. This suggests that coils can be added or removed by duplication or deletion of the DNA corresponding to one or more coils and explains how homologous proteins can have different numbers of coils. This process can also account for the evolution of other families of proteins such as the beta-rolls, the leucine-rich repeat proteins, the hexapeptide repeat family, two separate families of beta-helical antifreeze proteins and the spiral folds. These families need not be related to each other but will share features such as relative untwisted beta-sheets, stacking of similar residues and turns between beta-strands of approximately 90 degrees often stabilized by hydrogen bonding along the direction of the parallel beta-helix.Repetitive folds present special problems in the comparison of structures but offer attractive targets for structure prediction. The stacking of similar residues on a flat parallel beta-sheet may account for the formation of amyloid with beta-strands at right-angles to the fibril axis from many unrelated peptides.

Chloroplast glyceraldehyde-3-phosphate dehydrogenase contains a single disulfide bond located in the C-terminal extension to the B subunit

Mass mapping analysis based on cyanylation and CN-induced cleavage indicates that the two cysteine residues in the C-terminal extension of the B subunit of the light-activated pea leaf chloroplast glyceraldehyde-3-phosphate dehydrogenase form a disulfide bond. No evidence was found for a disulfide bond in the A subunit, nor was there any indication of a second disulfide bond in the B subunit. The availability of the structure of the extended glyceraldehyde-3-phosphate dehydrogenase from the archaeon Sulfolobus solfataricus allows modeling of the B subunit. As modeled, the two cysteine residues in the extension are positioned to form an interdomain disulfide cross-link.

Entropic stabilization of the tryptophan synthase alpha-subunit from a hyperthermophile, Pyrococcus furiosus. X-ray analysis and calorimetry

The structure of the tryptophan synthase alpha-subunit from Pyrococcus furiosus was determined by x-ray analysis at 2.0-A resolution, and its stability was examined by differential scanning calorimetry. Although the structure of the tryptophan synthase alpha(2)beta(2) complex from Salmonella typhimurium has been already determined, this is the first report of the structure of the alpha-subunit alone. The alpha-subunit from P. furiosus (Pf-alpha-subunit) lacked 12 and 6 residues at the N and C termini, respectively, and one residue each in two loop regions as compared with that from S. typhimurium (St-alpha-subunit), resulting in the absence of an N-terminal helix and the shortening of a C-terminal helix. The structure of the Pf-alpha-subunit was essentially similar to that of the St-alpha-subunit in the alpha(2)beta(2) complex. The differences between both structures were discussed in connection with the higher stability of the Pf-alpha-subunit and the complex formation of the alpha- and beta-subunits. Calorimetric results indicated that the Pf-alpha-subunit has extremely high thermostability and that its higher stability is caused by an entropic effect. On the basis of structural information of both proteins, we analyzed the contributions of each stabilization factor and could conclude that hydrophobic interactions in the protein interior do not contribute to the higher stability of the Pf-alpha-subunit. Rather, the increase in ion pairs, decrease in cavity volume, and entropic effects due to shortening of the polypeptide chain play important roles in extremely high stability in Pf-alpha-subunit.

Tiny TIM: a small, tetrameric, hyperthermostable triosephosphate isomerase

Comparative structural studies on proteins derived from organisms with growth optima ranging from 15 to 100 degrees C are beginning to shed light on the mechanisms of protein thermoadaptation. One means of sustaining hyperthermostability is for proteins to exist in higher oligomeric forms than their mesophilic homologues. Triosephosphate isomerase (TIM) is one of the most studied enzymes, whose fold represents one of nature's most common protein architectures. Most TIMs are dimers of approximately 250 amino acid residues per monomer. Here, we report the 2.7 A resolution crystal structure of the extremely thermostable TIM from Pyrococcus woesei, a hyperthermophilic archaeon growing optimally at 100 degrees C, representing the first archaeal TIM structure. P. woesei TIM exists as a tetramer comprising monomers of only 228 amino acid residues. Structural comparisons with other less thermostable TIMs show that although the central beta-barrel is largely conserved, severe pruning of several helices and truncation of some loops give rise to a much more compact monomer in the small hyperthermophilic TIM. The classical TIM dimer formation is conserved in P. woesei TIM. The extreme thermostability of PwTIM appears to be achieved by the creation of a compact tetramer where two classical TIM dimers interact via an extensive hydrophobic interface. The tetramer is formed through largely hydrophobic interactions between some of the pruned helical regions. The equivalent helical regions in less thermostable dimeric TIMs represent regions of high average temperature factor. The PwTIM seems to have removed these regions of potential instability in the formation of the tetramer. This study of PwTIM provides further support for the role of higher oligomerisation states in extreme thermal stabilisation.

Hyperthermophilic enzymes: sources, uses, and molecular mechanisms for thermostability

Enzymes synthesized by hyperthermophiles (bacteria and archaea with optimal growth temperatures of > 80 degrees C), also called hyperthermophilic enzymes, are typically thermostable (i.e., resistant to irreversible inactivation at high temperatures) and are optimally active at high temperatures. These enzymes share the same catalytic mechanisms with their mesophilic counterparts. When cloned and expressed in mesophilic hosts, hyperthermophilic enzymes usually retain their thermal properties, indicating that these properties are genetically encoded. Sequence alignments, amino acid content comparisons, crystal structure comparisons, and mutagenesis experiments indicate that hyperthermophilic enzymes are, indeed, very similar to their mesophilic homologues. No single mechanism is responsible for the remarkable stability of hyperthermophilic enzymes. Increased thermostability must be found, instead, in a small number of highly specific alterations that often do not obey any obvious traffic rules. After briefly discussing the diversity of hyperthermophilic organisms, this review concentrates on the remarkable thermostability of their enzymes. The biochemical and molecular properties of hyperthermophilic enzymes are described. Mechanisms responsible for protein inactivation are reviewed. The molecular mechanisms involved in protein thermostabilization are discussed, including ion pairs, hydrogen bonds, hydrophobic interactions, disulfide bridges, packing, decrease of the entropy of unfolding, and intersubunit interactions. Finally, current uses and potential applications of thermophilic and hyperthermophilic enzymes as research reagents and as catalysts for industrial processes are described.

Thermostability and thermoactivity of citrate synthases from the thermophilic and hyperthermophilic archaea, Thermoplasma acidophilum and Pyrococcus furiosus

Citrate synthases from Thermoplasma acidophilum (optimal growth at 55 degrees C) and Pyrococcus furiosus (100 degrees C) are homo-dimeric enzymes that show a high degree of structural homology with each other, and thermostabilities commensurate with the environmental temperatures in which their host cells are found. A comparison of their atomic structures with citrate synthases from mesophilic and psychrophilic organisms has indicated the potential importance of inter-subunit contacts for thermostability, and here we report the construction and analysis of site-directed mutants of the two citrate synthases to investigate the contribution of these interactions. Three sets of mutants were made: (a) chimeric mutants where the large (inter-subunit contact) and small (catalytic) domains of the T. acidophilum and P. furiosus enzymes were swapped; (b) mutants of the P. furiosus citrate synthase where the inter-subunit ionic network is disrupted; and (c) P. furiosus citrate synthase mutants in which the C-terminal arms that wrap around their partner subunits have been deleted. All three sets of mutant enzymes were expressed as recombinant proteins in Escherichia coli and were found to be catalytically active. Kinetic parameters and the dependence of catalytic activity on temperature were determined, and the stability of each enzyme was analysed by irreversible thermal inactivation experiments. The chimeric mutants indicate that the thermostability of the whole enzyme is largely determined by the origin of the large, inter-subunit domain, whereas the dependence of catalytic activity on temperature is a function of the small domain. Disruption of the inter-subunit ionic network and prevention of the C-terminal interactions both generated enzymes that were substantially less thermostable. Taken together, these data demonstrate the crucial importance of the subunit contacts to the stability of these oligomeric enzymes. Additionally, they also provide a clear distinction between thermostability and thermoactivity, showing that stability is necessary for, but does not guarantee, catalytic activity at elevated temperatures.

Crystal structure of saccharopine reductase from Magnaporthe grisea, an enzyme of the alpha-aminoadipate pathway of lysine biosynthesis

BACKGROUND

The biosynthesis of the essential amino acid lysine in higher fungi and cyanobacteria occurs via the alpha-aminoadipate pathway, which is completely different from the lysine biosynthetic pathway found in plants and bacteria. The penultimate reaction in the alpha-aminoadipate pathway is catalysed by NADPH-dependent saccharopine reductase. We set out to determine the structure of this enzyme as a first step in exploring the structural biology of fungal lysine biosynthesis.

RESULTS

We have determined the three-dimensional structure of saccharopine reductase from the plant pathogen Magnaporthe grisea in its apo form to 2.0 A resolution and as a ternary complex with NADPH and saccharopine to 2.1 A resolution. Saccharopine reductase is a homodimer, and each subunit consists of three domains, which are not consecutive in amino acid sequence. Domain I contains a variant of the Rossmann fold that binds NADPH. Domain II folds into a mixed seven-stranded beta sheet flanked by alpha helices and is involved in substrate binding and dimer formation. Domain III is all-helical. The structure analysis of the ternary complex reveals a large movement of domain III upon ligand binding. The active site is positioned in a cleft between the NADPH-binding domain and the second alpha/beta domain. Saccharopine is tightly bound to the enzyme via a number of hydrogen bonds to invariant amino acid residues.

CONCLUSIONS

On the basis of the structure of the ternary complex of saccharopine reductase, an enzymatic mechanism is proposed that includes the formation of a Schiff base as a key intermediate. Despite the lack of overall sequence homology, the fold of saccharopine reductase is similar to that observed in some enzymes of the diaminopimelate pathway of lysine biosynthesis in bacteria. These structural similarities suggest an evolutionary relationship between two different major families of amino acid biosynthetic pathway, the glutamate and aspartate families.

Structural analysis of glyceraldehyde 3-phosphate dehydrogenase from Escherichia coli: direct evidence of substrate binding and cofactor-induced conformational changes

The crystal structures of gyceraldehyde 3-phosphate dehydrogenase (GAPDH) from Escherichia coli have been determined in three different enzymatic states, NAD(+)-free, NAD(+)-bound, and hemiacetal intermediate. The NAD(+)-free structure reported here has been determined from monoclinic and tetragonal crystal forms. The conformational changes in GAPDH induced by cofactor binding are limited to the residues that bind the adenine moiety of NAD(+). Glyceraldehyde 3-phosphate (GAP), the substrate of GAPDH, binds to the enzyme with its C3 phosphate in a hydrophilic pocket, called the "new P(i)" site, which is different from the originally proposed binding site for inorganic phosphate. This observed location of the C3 phosphate is consistent with the flip-flop model proposed for the enzyme mechanism [Skarzynski, T., Moody, P. C., and Wonacott, A. J. (1987) J. Mol. Biol. 193, 171-187]. Via incorporation of the new P(i) site in this model, it is now proposed that the C3 phosphate of GAP initially binds at the new P(i) site and then flips to the P(s) site before hydride transfer. A superposition of NAD(+)-bound and hemiacetal intermediate structures reveals an interaction between the hydroxyl oxygen at the hemiacetal C1 of GAP and the nicotinamide ring. This finding suggests that the cofactor NAD(+) may stabilize the transition state oxyanion of the hemiacetal intermediate in support of the flip-flop model for GAP binding.