Structure determination of the glutamate dehydrogenase from the hyperthermophile Thermococcus litoralis and its comparison with that from Pyrococcus furiosus

Structural Evolution of Primate Glutamate Dehydrogenase 2 as Revealed by In Silico Predictions and Experimentally Determined Structures

Glutamate dehydrogenase (GDH) interconverts glutamate to a-ketoglutarate and ammonia, interconnecting amino acid and carbohydrate metabolism. In humans, two functional GDH genes, GLUD1 and GLUD2, encode for hGDH1 and hGDH2, respectively. GLUD2 evolved from retrotransposition of the GLUD1 gene in the common ancestor of modern apes. These two isoenzymes are involved in the pathophysiology of human metabolic, neoplastic, and neurodegenerative disorders. The 3D structures of hGDH1 and hGDH2 have been experimentally determined; however, no information is available about the path of GDH2 structure changes during primate evolution. Here, we compare the structures predicted by the AlphaFold Colab method for the GDH2 enzyme of modern apes and their extinct primate ancestors. Also, we analyze the individual effect of amino acid substitutions emerging during primate evolution. Our most important finding is that the predicted structure of GDH2 in the common ancestor of apes was the steppingstone for the structural evolution of primate GDH2s. Two changes with a strong functional impact occurring at the first evolutionary step, Arg443Ser and Gly456Ala, had a destabilizing and stabilizing effect, respectively, making this step the most important one. Subsequently, GDH2 underwent additional modifications that fine-tuned its enzymatic properties to adapt to the functional needs of modern-day primate tissues.

Understanding Life at High Temperatures: Relationships of Molecular Channels in Enzymes of Methanogenic Archaea and Their Growth Temperatures

Analyses of protein structures have shown the existence of molecular channels in enzymes from Prokaryotes. Those molecular channels suggest a critical role of spatial voids in proteins, above all, in those enzymes functioning under high temperature. It is expected that these spaces within the protein structure are required to access the active site and to maximize availability and thermal stability of their substrates and cofactors. Interestingly, numerous substrates and cofactors have been reported to be highly temperature-sensitive biomolecules. Methanogens represent a singular phylogenetic group of Archaea that performs anaerobic respiration producing methane during growth. Methanogens inhabit a variety of environments including the full range of temperatures for the known living forms. Herein, we carry out a dimensional analysis of molecular tunnels in key enzymes of the methanogenic pathway from methanogenic Archaea growing optimally over a broad temperature range. We aim to determine whether the dimensions of the molecular tunnels are critical for those enzymes from thermophiles. Results showed that at increasing growth temperature the dimensions of molecular tunnels in the enzymes methyl-coenzyme M reductase and heterodisulfide reductase become increasingly restrictive and present strict limits at the highest growth temperatures, i.e., for hyperthermophilic methanogens. However, growth at lower temperature allows a wide dimensional range for the molecular spaces in these enzymes. This is in agreement with previous suggestions on a potential major role of molecular tunnels to maintain biomolecule stability and activity of some enzymes in microorganisms growing at high temperatures. These results contribute to better understand archaeal growth at high temperatures. Furthermore, an optimization of the dimensions of molecular tunnels would represent an important adaptation required to maintain the activity of key enzymes of the methanogenic pathway for those methanogens growing optimally at high temperatures.

Mapping the Intramolecular Communications among Different Glutamate Dehydrogenase States Using Molecular Dynamics

Glutamate dehydrogenase (GDH) is a ubiquitous enzyme that catalyzes the reversible oxidative deamination of glutamate to α-ketoglutarate. It acts as an important branch-point enzyme between carbon and nitrogen metabolisms. Due to the multifaceted roles of GDH in cancer, hyperinsulinism/hyperammonemia, and central nervous system development and pathologies, tight control of its activity is necessitated. To date, several GDH structures have been solved in its closed form; however, intrinsic structural information in its open and apo forms are still deficient. Moreover, the allosteric communications and conformational changes taking place in the three different GDH states are not well studied. To mitigate these drawbacks, we applied unbiased molecular dynamic simulations (MD) and network analysis to three different GDH states i.e., apo, active, and inactive forms, for investigating their modulatory mechanisms. In this paper, based on MD and network analysis, crucial residues important for signal transduction, conformational changes, and maps of information flow among the different GDH states were elucidated. Moreover, with the recent findings of allosteric modulators, an allosteric wiring illustration of GDH intramolecular signal transductions would be of paramount importance to obtain the process of this enzyme regulation. The structural insights gained from this study will pave way for large-scale screening of GDH regulators and could support researchers in the design and development of new and potent GDH ligands.

Structural Studies of Glutamate Dehydrogenase (Isoform 1) From Arabidopsis thaliana, an Important Enzyme at the Branch-Point Between Carbon and Nitrogen Metabolism

Glutamate dehydrogenase (GDH) releases ammonia in a reversible NAD(P)⁺-dependent oxidative deamination of glutamate that yields 2-oxoglutarate (2OG). In current perception, GDH contributes to Glu homeostasis and plays a significant role at the junction of carbon and nitrogen assimilation pathways. GDHs are members of a superfamily of ELFV (Glu/Leu/Phe/Val) amino acid dehydrogenases and are subdivided into three subclasses, based on coenzyme specificity: NAD⁺-specific, NAD⁺/NADP⁺ dual-specific, and NADP⁺-specific. We determined in this work that the mitochondrial AtGDH1 isozyme from A. thaliana is NAD⁺-specific. Altogether, A. thaliana expresses three GDH isozymes (AtGDH1-3) targeted to mitochondria, of which AtGDH2 has an extra EF-hand motif and is stimulated by calcium. Our enzymatic assays of AtGDH1 established that its sensitivity to calcium is negligible. In vivo the AtGDH1-3 enzymes form homo- and heterohexamers of varied composition. We solved the crystal structure of recombinant AtGDH1 in the apo-form and in complex with NAD⁺ at 2.59 and 2.03 Å resolution, respectively. We demonstrate also that both in the apo form and in 1:1 complex with NAD⁺, it forms D ₃-symmetric homohexamers. A subunit of AtGDH1 consists of domain I, which is involved in hexamer formation and substrate binding, and of domain II which binds coenzyme. Most of the subunits in our crystal structures, including those in NAD⁺ complex, are in open conformation, with domain II forming a large (albeit variable) angle with domain I. One of the subunits of the AtGDH1-NAD⁺ hexamer contains a serendipitous 2OG molecule in the active site, causing a dramatic (∼25°) closure of the domains. We provide convincing evidence that the N-terminal peptide preceding domain I is a mitochondrial targeting signal, with a predicted cleavage site for mitochondrial processing peptidase (MPP) at Leu17-Leu18 that is followed by an unexpected potassium coordination site (Ser27, Ile30). We also identified several MPD [(+/-)-2-methyl-2,4-pentanediol] binding sites with conserved sequence. Although AtGDH1 is insensitive to MPD in our assays, the observation of druggable sites opens a potential for non-competitive herbicide design.

Molecular Tunnels in Enzymes and Thermophily: A Case Study on the Relationship to Growth Temperature

Developments in protein expression, analysis and computational capabilities are decisively contributing to a better understanding of the structure of proteins and their relationship to function. Proteins are known to be adapted to the growth rate of microorganisms and some microorganisms (named (hyper)thermophiles) thrive optimally at high temperatures, even above 100 °C. Nevertheless, some biomolecules show great instability at high temperatures and some of them are universal and required substrates and cofactors in multiple enzymatic reactions for all (both mesophiles and thermophiles) living cells. Only a few possibilities have been pointed out to explain the mechanisms that thermophiles use to successfully thrive under high temperatures. As one of these alternatives, the role of molecular tunnels or channels in enzymes has been suggested but remains to be elucidated. This study presents an analysis of channels in proteins (i.e., substrate tunnels), comparing two different protein types, glutamate dehydrogenase and glutamine phosphoribosylpyrophosphate amidotransferase, which are supposed to present a different strategy on the requirement for substrate tunnels with low and high needs for tunneling, respectively. The search and comparison of molecular tunnels in these proteins from microorganisms thriving optimally from 15 °C to 100 °C suggested that those tunnels in (hyper)thermophiles are required and optimized to specific dimensions at high temperatures for the enzyme glutamine phosphoribosylpyrophosphate amidotransferase. For the enzyme glutamate dehydrogenase, a reduction of empty spaces within the protein could explain the optimization at increasing temperatures. This analysis provides further evidence on molecular channeling as a feasible mechanism in hyperthermophiles with multiple relevant consequences contributing to better understand how they live under those extreme conditions.

Structural basis for the catalytic mechanism and α-ketoglutarate cooperativity of glutamate dehydrogenase

Glutamate dehydrogenase (GDH) is a key enzyme connecting carbon and nitrogen metabolism in all living organisms. Despite extensive studies on GDHs from both prokaryotic and eukaryotic organisms in the last 40 years, the structural basis of the catalytic features of this enzyme remains incomplete. This study reports the structural basis of the GDH catalytic mechanism and allosteric behavior. We determined the first high-resolution crystal structures of glutamate dehydrogenase from the fungus Aspergillus niger (AnGDH), a unique NADP⁺-dependent allosteric enzyme that is forward-inhibited by the formation of mixed disulfide. We determined the structures of the active enzyme in its apo form and in binary/ternary complexes with bound substrate (α-ketoglutarate), inhibitor (isophthalate), coenzyme (NADPH), or two reaction intermediates (α-iminoglutarate and 2-amino-2-hydroxyglutarate). The structure of the forward-inhibited enzyme (fiAnGDH) was also determined. The hexameric AnGDH had three open subunits at one side and three partially closed protomers at the other, a configuration not previously reported. The AnGDH hexamers having subunits with different conformations indicated that its α-ketoglutarate-dependent homotropic cooperativity follows the Monod-Wyman-Changeux (MWC) model. Moreover, the position of the water attached to Asp-154 and Gly-153 defined the previously unresolved ammonium ion-binding pocket, and the binding site for the 2'-phosphate group of the coenzyme was also better defined by our structural data. Additional structural and mutagenesis experiments identified the residues essential for coenzyme recognition. This study reveals the structural features responsible for positioning α-ketoglutarate, NADPH, ammonium ion, and the reaction intermediates in the GDH active site.

Purification, crystallization and preliminary X-ray diffraction analysis of NADP-dependent glutamate dehydrogenase from Aspergillus niger

Glutamate dehydrogenase (GDH) catalyzes the NAD-dependent or NADP-dependent oxidative deamination of L-glutamate to 2-oxoglutarate and ammonia. This important reversible reaction establishes the link between carbon and nitrogen metabolism. In this study, Aspergillus niger NADP-GDH (AnGDH) has been overexpressed and purified. Purified AnGDH, with a high specific activity of 631.1 units per milligram of protein, was crystallized and the crystal diffracted to 2.9 Å resolution using a home X-ray source. Preliminary analysis of the X-ray diffraction data showed that the crystal belonged to space group R32, with unit-cell parameters a=b=173.8, c=241.5 Å, α=β=90, γ=120°. The crystals exhibited an unusually high solvent content (83.0%) and had only one molecule in the asymmetric unit. Initial phases were obtained by molecular replacement, and model building and structure refinement of AnGDH are in progress.

Biochemical characterization of two glutamate dehydrogenases with different cofactor specificities from a hyperthermophilic archaeon Pyrobaculum calidifontis

Two putative glutamate dehydrogenase (GDH) genes (pcal_1031 and pcal_1606) were found in a sulfur-dependent hyperthermophilic archaeon, Pyrobaculum calidifontis. The two genes were then expressed in Escherichia coli, and both of the recombinant gene products showed GDH activity. The two enzymes were then purified to homogeneity and characterized in detail. Although both purified GDHs had a hexameric structure and neither exhibited allosteric regulation, they showed different coenzyme specificities: one was specific for NAD(+), the other for NADP(+) and different heat activation mechanisms. In addition, there was little difference in the kinetic constants, optimal temperature, thermal stability, optimal pH and pH stability between the two enzymes. The overall sequence identity between the two proteins was very high (81%), but was not high in the region recognizing the 2' position of the adenine ribose moiety, which is responsible for coenzyme specificity. This is the first report on the identification of two GDHs with different coenzyme specificities from a single hyperthermophilic archaeon and the definition of their basic in vitro properties.

L-aspartate dehydrogenase: features and applications

L-amino acid dehydrogenases are a group of enzymes that catalyze the reversible oxidative deamination of L-amino acids to their corresponding 2-oxoacids, using either nicotinamide adenine dinucleotide (NAD(+)) or nicotinamide adenine dinucleotide phosphate (NADP(+)) as cofactors. These enzymes have been studied widely because of their potential applications in the synthesis of amino acids for use in production of pharmaceutical peptides, herbicides and insecticides, in biosensors or diagnostic kits, and development of coenzyme regeneration systems for industrial processes. This article presents a review of the currently available data about the recently discovered amino acid dehydrogenase superfamily member L-aspartate dehydrogenase (L-AspDH), their relevant catalytic properties and speculated physiological roles, and potential for biotechnological applications. The proposed classification of L-AspDH on the basis of bioinformatic information and potential role in vivo into NadB (NAD biosynthesis-related) and non-NadB type is unique. In particular, the mesophilic non-NadB type L-AspDH is a novel group of amino acid dehydrogenases with great promise as potential industrial biocatalysts owing to their relatively high catalytic properties at room temperature. Considering that only a few L-AspDH homologs have been characterized so far, identification and prodigious enzymological research of the new members will be necessary to shed light on the gray areas pertaining to these enzymes.

Crystal structure of NAD+-dependent Peptoniphilus asaccharolyticus glutamate dehydrogenase reveals determinants of cofactor specificity

Glutamate dehydrogenases (EC 1.4.1.2-4) catalyse the oxidative deamination of l-glutamate to α-ketoglutarate using NAD(P) as a cofactor. The bacterial enzymes are hexamers and each polypeptide consists of an N-terminal substrate-binding (Domain I) followed by a C-terminal cofactor-binding segment (Domain II). The reaction takes place at the junction of the two domains, which move as rigid bodies and are presumed to narrow the cleft during catalysis. Distinct signature sequences in the nucleotide-binding domain have been linked to NAD(+) vs. NADP(+) specificity, but they are not unambiguous predictors of cofactor preferences. Here, we have determined the crystal structure of NAD(+)-specific Peptoniphilus asaccharolyticus glutamate dehydrogenase in the apo state. The poor quality of native crystals was resolved by derivatization with selenomethionine, and the structure was solved by single-wavelength anomalous diffraction methods. The structure reveals an open catalytic cleft in the absence of substrate and cofactor. Modeling of NAD(+) in Domain II suggests that a hydrophobic pocket and polar residues contribute to nucleotide specificity. Mutagenesis and isothermal titration calorimetry studies of a critical glutamate at the P7 position of the core fingerprint confirms its role in NAD(+) binding. Finally, the cofactor binding site is compared with bacterial and mammalian enzymes to understand how the amino acid sequences and three-dimensional structures may distinguish between NAD(+) vs. NADP(+) recognition.

Functional dissection of a trigger enzyme: mutations of the bacillus subtilis glutamate dehydrogenase RocG that affect differentially its catalytic activity and regulatory properties

Any signal transduction requires communication between a sensory component and an effector. Some enzymes engage in signal perception and transduction, as well as in catalysis, and these proteins are known as "trigger" enzymes. In this report, we detail the trigger properties of RocG, the glutamate dehydrogenase of Bacillus subtilis. RocG not only deaminates the key metabolite glutamate to form alpha-ketoglutarate but also interacts directly with GltC, a LysR-type transcription factor that regulates glutamate biosynthesis from alpha-ketoglutarate, thus linking the two metabolic pathways. We have isolated mutants of RocG that separate the two functions. Several mutations resulted in permanent inactivation of GltC as long as a source of glutamate was present. These RocG proteins have lost their ability to catabolize glutamate due to a strongly reduced affinity for glutamate. The second class of mutants is exemplified by the replacement of aspartate residue 122 by asparagine. This mutant protein has retained enzymatic activity but has lost the ability to control the activity of GltC. Crystal structures of glutamate dehydrogenases that permit a molecular explanation of the properties of the various mutants are presented. Specifically, we may propose that D122N replacement affects the surface of RocG. Our data provide evidence for a correlation between the enzymatic activity of RocG and its ability to inactivate GltC, and thus give insights into the mechanism that couples the enzymatic activity of a trigger enzyme to its regulatory function.

Temperature dependence of protein motions in a thermophilic dihydrofolate reductase and its relationship to catalytic efficiency

We report hydrogen deuterium exchange by mass spectrometry (HDX-MS) as a function of temperature in a thermophilic dihydrofolate reductase from Bacillus stearothermophilus (Bs-DHFR). Protein stability, probed with circular dichroism, established an accessible temperature range of 10 degrees C to 55 degrees C for the interrogation of HDX-MS. Although both the rate and extent of HDX are sensitive to temperature, the majority of peptides showed rapid kinetics of exchange, allowing us to focus on plateau values for the maximal extent of exchange at each temperature. Arrhenius plots of the ratio of hydrogens exchanged at 5 h normalized to the number of exchangeable hydrogens vs. 1/T provides an estimate for the apparent enthalpic change of local unfolding, DeltaH degrees (unf(avg)). Most regions in the enzyme show DeltaH degrees (unf(avg)) < or = 2.0 kcal/mol, close to the value of kT; by contrast, significantly elevated values for DeltaH degrees (unf(avg)) are observed in regions within the core of protein that contain the cofactor and substrate-binding sites. Our technique introduces a new strategy for probing the temperature dependence of local protein unfolding within native proteins. These findings are discussed in the context of the demonstrated role for nuclear tunneling in hydride transfer from NADPH to dihydrofolate, and relate the observed enthalpic changes to two classes of motion, preorganization and reorganization, that have been proposed to control the efficiency of hydrogenic wave function overlap. Our findings suggest that the enthalpic contribution to the heavy atom environmental reorganizations controlling the hydrogenic wave function overlap will be dominated by regions of the protein proximal to the bound cofactor and substrate.

Thermal stability of isocitrate dehydrogenase from Archaeoglobus fulgidus studied by crystal structure analysis and engineering of chimers

Isocitrate dehydrogenase from Archaeoglobus fulgidus (AfIDH) has an apparent melting temperature (T(m)) of 98.5 degrees C. To identify the structural features involved in thermal stabilization of AfIDH, the structure was solved to 2.5 A resolution. AfIDH was strikingly similar to mesophilic IDH from Escherichia coli (EcIDH) and displayed almost the same number of ion pairs and ionic networks. However, two unique inter-domain networks were present in AfIDH; one three-membered ionic network between the large and the small domain and one four-membered ionic network between the clasp and the small domain. The latter ionic network was presumably reduced in size when the clasp domain of AfIDH was swapped with that of EcIDH and the T (m) decreased by 18 degrees C. Contrarily, EcIDH was only stabilized by 4 degrees C by the clasp domain of AfIDH, a result probably due to the introduction of a unique inter-subunit aromatic cluster in AfIDH that may strengthen the dimeric interface in this enzyme. A unique aromatic cluster was identified close to the N-terminus of AfIDH that could provide additional stabilization of this region. Common and unique heat adaptive traits of AfIDH with those recently observed for hyperthermophilic IDH from Aeropyrum pernix (ApIDH) and Thermotoga maritima (TmIDH) are discussed herein.

Explanation of the stability of thermophilic proteins based on unique micromorphology

Two mesophilic/thermophilic variants of the G-domain of the elongation factor Tu were studied via molecular dynamics simulations. By analyzing the simulation data via the Voronoi space tessellation, we have found that the two proteins have the same macromolecular packing, while the water-exposed surface area is larger for the thermophile. A larger coordination with water is probably due to a peculiar corrugation of the exposed surface of this species. From an enthalpic point of view, the thermophile shows a larger number of intramolecular hydrogen bonds, stronger electrostatic interactions, and a flatter free-energy landscape. Overall, the data suggest that the specific hydration state enhances macromolecular fluctuations but, at the same time, increases thermal stability.

The crystal structure of Plasmodium falciparum glutamate dehydrogenase, a putative target for novel antimalarial drugs

Plasmodium falciparum is the main causative agent of tropical malaria, the most severe parasitic disease in the world. Growing resistance of Plasmodia towards available drugs is an increasing problem in countries where malaria is endemic. As Plasmodia are sensitive to oxidative stress, augmenting this in the parasite represents a promising principle for the development of novel antimalarial drugs. The NADP-dependent glutamate dehydrogenase (GDH) of P.falciparum is largely responsible for the production of NADPH in the parasite, which in turn serves as electron source for the antioxidative enzymes glutathione reductase and thioredoxin reductase. As GDH does not occur in the host erythrocyte, GDH is a particularly attractive target for drug therapy. The three-dimensional structure of P.falciparum GDH in the unligated state has been determined by X-ray crystallography to a resolution of 2.7A. Compared to the mammalian enzymes, two amino acid residues are exchanged in the putative active site of the parasite GDH. The most obvious differences between parasite and human GDH are the subunit interfaces of the hexameric proteins. In the parasite protein, several salt-bridges mediate contacts between the subunits whereas in the human enzyme these interactions are mainly of hydrophobic nature. Furthermore, P.falciparum GDH possesses a unique N-terminal extension that does not occur in any other GDH sequence so far studied. These findings might be exploited for the design of peptidomimetics capable of disrupting the oligomeric organisation of the parasite enzyme.

Asparaginyl deamidation in two glutamate dehydrogenase isoenzymes from Saccharomyces cerevisiae

The non-enzymatic deamidation of asparaginyl residues is a major source of spontaneous damage of several proteins under physiological conditions. In many cases, deamidation and isoaspartyl formation alters the biological activity or stability of the native polypeptide. Rates of deamidation of particular residues depend on many factors including protein structure and solvent exposure. Here, we investigated the spontaneous deamidation of the two NADP-glutamate dehydrogenase isoenzymes from Saccharomyces cerevisiae, which have different kinetic properties and are differentially expressed in this yeast. Our results show that Asn54, present in Gdh3p but missing in the GDH1-encoded homologue, is readily deamidated in vitro under alkaline conditions. Relative to the native enzyme, deamidated Gdh3p shows reduced protein stability. The different deamidation rates of the two isoenzymes could explain to some extent, the relative in vivo instability of the allosteric Gdh3p enzyme, compared to that of Gdh1p. It is thus possible that spontaneous asparaginyl modification could play a role in the metabolic regulation of ammonium assimilation and glutamate biosynthesis.

Design and Construction of Highly Stable, Protease-resistant Chimeric Avidins

The chicken avidin gene family consists of avidin and seven separate avidin-related genes (AVRs) 1-7. Avidin protein is a widely used biochemical tool, whereas the other family members have only recently been produced as recombinant proteins and characterized. In our previous study, AVR4 was found to be the most stable biotin binding protein thus far characterized (T(m) = 106.4 degrees C). In this study, we studied further the biotin-binding properties of AVR4. A decrease in the energy barrier between the biotin-bound and unbound state of AVR4 was observed when compared with that of avidin. The high resolution structure of AVR4 facilitated comparison of the structural details of avidin and AVR4. In the present study, we used the information obtained from these comparative studies to transfer the stability and functional properties of AVR4 to avidin. A chimeric avidin protein, ChiAVD, containing a 21-amino acid segment of AVR4 was found to be significantly more stable (T(m) = 96.5 degrees C) than native avidin (T(m) = 83.5 degrees C), and its biotin-binding properties resembled those of AVR4. Optimization of a crucial subunit interface of avidin by an AVR4-inspired point mutation, I117Y, significantly increased the thermostability of the avidin mutant (T(m) = 97.5 degrees C) without compromising its high biotin-binding properties. By combining these two modifications, a hyperthermostable ChiAVD(I117Y) was constructed (T(m) = 111.1 degrees C). This study provides an example of rational protein engineering in which another member of the protein family has been utilized as a source in the optimization of selected properties.

The first crystal structure of hyperthermostable NAD-dependent glutamate dehydrogenase from Pyrobaculum islandicum

The extremely thermostable NAD-dependent glutamate dehydrogenase (NAD-GluDH) from Pyrobaculum islandicum, a member of the Crenarchaeota, was crystallized, and its 3D structure has been determined by X-ray diffraction methods. The homohexameric structure of Pb. islandicum glutamate dehydrogenase (Pis-GluDH) was solved and refined at a resolution of 2.9A with a crystallographic R-factor of 19.9% (Rfree 26.0%). The structure indicates that each subunit consists of two domains separated by a deep cleft containing an active site. The secondary structural elements and catalytically important residues of the enzyme were highly conserved among the NAD(P)-dependent GluDHs from other sources. A structural comparison of Pis-GluDH with other NAD(P)-dependent GluDHs suggests that a significant difference in the alpha8-loop-alpha9 region of this enzyme is associated with its coenzyme specificity. From the analysis of the 3D structure, hydrophobic interactions between intersubunits were found to be important features for the enzyme oligomerization. It has been reported that Pis-GluDH is highly thermostable, like the GluDH of the hyperthermophilic archaeum Pyrococcus furiosus, and the increase in the intersubunit ion pair networks is responsible for the extreme thermostability of the Pc. furiosus enzyme. However, the number of intersubunit ion pairs in the Pis-GluDH molecules is much smaller than those of the Pc. furiosus GluDH. The number of hydrophobic interactions at the intersubunit interfaces were increased and responsible for the extremely high thermostability. This indicates that the major molecular strategy for high thermostability of the GluDHs may be different for each hyperthermophile.

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.

The upper temperature of life--where do we draw the line?

A new archaeal isolate has been reported that is capable of growing at up to 121 degrees C. The hyperthermophile, dubbed strain 121, grows chemoautotrophically using formate as an electron donor and FeIII as an electron acceptor and is closely related to members of the archaeal genera Pyrodictium and Pyrobaculum. Although the reported maximum growth temperature of strain 121 is 8 degrees C higher than the previous record holder (Pyrolobus fumarii; Tmax = 113 degrees C), the two organisms have virtually the same optimal growth temperatures.

An evolutionary route to xylanase process fitness

Directed evolution technologies were used to selectively improve the stability of an enzyme without compromising its catalytic activity. In particular, this article describes the tandem use of two evolution strategies to evolve a xylanase, rendering it tolerant to temperatures in excess of 90 degrees C. A library of all possible 19 amino acid substitutions at each residue position was generated and screened for activity after a temperature challenge. Nine single amino acid residue changes were identified that enhanced thermostability. All 512 possible combinatorial variants of the nine mutations were then generated and screened for improved thermal tolerance under stringent conditions. The screen yielded eleven variants with substantially improved thermal tolerance. Denaturation temperature transition midpoints were increased from 61 degrees C to as high as 96 degrees C. The use of two evolution strategies in combination enabled the rapid discovery of the enzyme variant with the highest degree of fitness (greater thermal tolerance and activity relative to the wild-type parent).

Using motif-based methods in multiple genome analyses: a case study comparing orthologous mesophilic and thermophilic proteins

Protein motifs represent highly conserved regions within protein families and are generally accepted to describe critical regions required for protein stability and/or function. In this comprehensive analysis, we present a robust, unique approach to identify and compare corresponding mesophilic and thermophilic sequence motifs between all orthologous proteins within 44 microbial genomes. Motif similarity is determined through global sequence alignment of mesophilic and thermophilic motif pairs, which are identified by a greedy algorithm. Our results reveal only modest correlation between motif and overall sequence similarity, highlighting the rationale of motif-based approaches in comprehensive multigenome comparisons. Conserved mutations reflect previously suggested physiochemical principles for conferring thermostability. Additionally, comparisons between corresponding mesophilic and thermophilic motif pairs provide key biochemical insights related to thermostability and can be used to test the evolutionary robustness of individual structural comparisons. We demonstrate the ability of our unique approach to provide key insights in two examples: the TATA-box binding protein and glutamate dehydrogenase families. In the latter example, conserved mutations hint at novel origins leading to structural stability differences within the hexamer structures. Additionally, we present amino acid composition data and average protein length comparisons for all 44 microbial genomes.

Heat effect on the structure and activity of the recombinant glutamate dehydrogenase from a hyperthermophilic archaeon Pyrococcus horikoshii

Glutamate dehydrogenase from Pyrococcus horikoshii (Pho-GDH) was cloned and overexpressed in Escherichia coli. The cloned enzyme with His-tag was purified to homogeneity by affinity chromatography and shown to be a hexamer enzyme of 290+/-8 kDa (subunit mass 48 kDa). Its optimal pH and temperature were 7.6 and 90 degrees C, respectively. The purified enzyme has outstanding thermostability (the half-life for thermal inactivation at 100 degrees C was 4 h). The enzyme shows strict specificity for 2-oxoglutarate and L-glutamate and requires NAD(P)H and NADP as cofactors but it does not reveal activity on NAD as cofactor. K(m) values of the recombinant enzyme are comparable for both substrates: 0.2 mM for L-glutamate and 0.53 mM for 2-oxoglutarate. The enzyme was activated by heating at 80 degrees C for 1 h, which was accompanied by the formation of its active conformation. Circular dichroism and fluorescence spectra show that the active conformation is heat-inducible and time-dependent.

Structural and thermodynamic studies on a salt-bridge triad in the NADP-binding domain of glutamate dehydrogenase from Thermotoga maritima: cooperativity and electrostatic contribution to stability

Cooperative interactions within ion-pair networks of hyperthermostable proteins are thought to be a major determinant for extreme protein stability. While the favorable thermodynamic contributions of optimized electrostatics in general as well as those of pairwise interactions have been documented, cooperativity between pairwise interactions has not yet been studied thermodynamically in proteins from hyperthermophiles. In this study we use the isolated cofactor binding domain of glutamate dehydrogenase from the hyperthermophilic bacterium Thermotoga maritima to analyze pairwise and cooperative interactions within the salt-bridge triad Arg190-Glu231-Lys193. The X-ray structure of the domain was solved at 1.43 A and reveals the salt-bridge network with surrounding solvent molecules in detail. All three participating charges in the network were mutated to alanine in all combinations. The X-ray structure of the variant lacking all three charges reveals that the removal of the side chains has no effect on the overall conformation of the protein. Using solvent denaturation and thermodynamic cycles, the interaction energies between each pair of residues in the network were determined in the presence and in the absence of the third residue. Both the Arg190-Glu231 ion pair and the Lys193-Glu231 salt bridge in the absence of the third residue, contribute favorably to the free energy for unfolding of the domain in urea. Using guanidinium chloride as denaturant reveals a strong cooperativity between the two ion-pair interactions, the presence of the second ion pair converts the first interaction from destabilizing into stabilizing by as much as 1.09 kcal/mol. The different energetics of the salt-bridge triad in urea and GdmCl are discussed with reference to the observed anion binding in the crystal structure at high ionic strength and their possible role in a highly charged, high-temperature environment such as the cytoplasm of hyperthermophiles.

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.

Structural basis for thermophilic protein stability: structures of thermophilic and mesophilic malate dehydrogenases

The three-dimensional structure of four malate dehydrogenases (MDH) from thermophilic and mesophilic phototropic bacteria have been determined by X-ray crystallography and the corresponding structures compared. In contrast to the dimeric quaternary structure of most MDHs, these MDHs are tetramers and are structurally related to tetrameric malate dehydrogenases from Archaea and to lactate dehydrogenases. The tetramers are dimers of dimers, where the structures of each subunit and the dimers are similar to the dimeric malate dehydrogenases. The difference in optimal growth temperature of the corresponding organisms is relatively small, ranging from 32 to 55 degrees C. Nevertheless, on the basis of the four crystal structures, a number of factors that are likely to contribute to the relative thermostability in the present series have been identified. It appears from the results obtained, that the difference in thermostability between MDH from the mesophilic Chlorobium vibrioforme on one hand and from the moderate thermophile Chlorobium tepidum on the other hand is mainly due to the presence of polar residues that form additional hydrogen bonds within each subunit. Furthermore, for the even more thermostable Chloroflexus aurantiacus MDH, the use of charged residues to form additional ionic interactions across the dimer-dimer interface is favored. This enzyme has a favorable intercalation of His-Trp as well as additional aromatic contacts at the monomer-monomer interface in each dimer. A structural alignment of tetrameric and dimeric prokaryotic MDHs reveal that structural elements that differ among dimeric and tetrameric MDHs are located in a few loop regions.

The crystal structure of a hyper-thermophilic carboxylesterase from the archaeon Archaeoglobus fulgidus

The crystal structure of AFEST, a novel hyper-thermophilic carboxylesterase from the archaeon Archaeoglobus fulgidus, complexed with a sulphonyl derivative, has been determined and refined to 2.2 A resolution. This enzyme, which has recently been classified as a member of the hormone- sensitive-lipase (H) group of the esterase/lipase superfamily, presents a canonical alpha/beta hydrolase core, shielded on the C-terminal side by a cap region composed of five alpha-helices. It contains the catalytic triad Ser160, His285 and Asp255, whereby the nucleophile is covalently modified and the oxyanion hole formed by Gly88, Gly89 and Ala161. A structural comparison of AFEST with its mesophilic and thermophilic homologues, Brefeldin A esterase from Bacillus subtilis (BFAE) and EST2 from Alicyclobacillus acidocaldarius, reveals an increase in the number of intramolecular ion pairs and secondary structure content, as well as a significant reduction in loop extensions and ratio of hydrophobic to charged surface area. The variety of structural differences suggests possible strategies for thermostabilization of lipases and esterases with potential industrial applications.

Increasing the thermostability of Flavobacterium meningosepticum glycerol kinase by changing Ser329 to Asp in the subunit interface region

The thermostability enhancement of Flavobacterium meningosepticum glycerol kinase (FGK) by random mutagenesis in the subunit interface region was investigated. A single Escherichia coli transformant, which produced a more thermostable glycerol kinase than the parent enzyme, was obtained. The nucleotide sequence of the gene of the mutant enzyme (FGK2615) was determined, and the four amino acid replacements were identified as Glu327 to Asp, Ser329 to Asp, Thr330 to Ala and Ser334 to Lys. Although the properties of FGK2615 were fundamentally similar to those of the parent enzyme, the thermostability and Km for ATP had changed. The thermostability of FGK2615 was apparently increased; the temperature at which the enzyme activity is inactivated by 50% for a 30-min incubation of FGK2615 was determined to be 72.1 degrees C which was 3.1 degrees C higher than that of the parent FGK. Four additional mutants each having a single amino acid replacement (Glu327 to Asp, Ser329 to Asp, Thr330 to Ala and Ser334 to Lys) were prepared and their thermostability and Km for substrates were evaluated. The effect of the substitution of Ser329 to Asp is discussed.

Enzymes from extremophiles

The industrial application of enzymes that can withstand harsh conditions has greatly increased over the past decade. This is mainly a result of the discovery of novel enzymes from extremophilic microorganisms. Recent advances in the study of extremozymes point to the acceleration of this trend. In particular, enzymes from thermophilic organisms have found the most practical commercial use to date because of their overall inherent stability. This has also led to a greater understanding of stability factors involved in adaptation of these enzymes to their unusual environments.

Large-scale domain movements and hydration structure changes in the active-site cleft of unligated glutamate dehydrogenase from Thermococcus profundus studied by cryogenic X-ray crystal structure analysis and small-angle X-ray scattering

Here we describe the large-scale domain movements and hydration structure changes in the active-site cleft of unligated glutamate dehydrogenase. Glutamate dehydrogenase from Thermococcus profundus is composed of six identical subunits of M(r) 46K, each with two distinct domains of roughly equal size separated by a large active-site cleft. The enzyme in the unligated state was crystallized so that one hexamer occupied a crystallographic asymmetric unit, and the crystal structure of the hexamer was solved and refined at a resolution of 2.25 A with a crystallographic R-factor of 0.190. In that structure, the six subunits displayed significant conformational variations with respect to the orientations of the two domains. The variation was most likely explained as a hinge-bending motion caused by small changes in the main chain torsion angle of the residue composing a loop connecting the two domains. Small-angle X-ray scattering profiles both at 293 and 338 K suggested that the apparent molecular size of the hexamer was slightly larger in solution than in the crystalline state. These results led us to the conclusion that (i) the spontaneous domain motion was the property of the enzyme in solution, (ii) the domain motion was trapped in the crystallization process through different modes of crystal contacts, and (iii) the magnitude of the motion in solution was greater than that observed in the crystal structure. The present cryogenic diffraction experiment enabled us to identify 1931 hydration water molecules around the hexamer. The hydration structures around the subunits exhibited significant changes in accord with the degree of the domain movement. In particular, the hydration water molecules in the active-site cleft were rearranged markedly through migrations between specific hydration sites in coupling strongly with the domain movement. We discussed the cooperative dynamics between the domain motion and the hydration structure changes in the active site of the enzyme. The present study provides the first example of a visualized hydration structure varying transiently with the dynamic movements of enzymes and may form a new concept of the dynamics of multidomain enzymes in solution.

A new class of glutamate dehydrogenases (GDH). Biochemical and genetic characterization of the first member, the AMP-requiring NAD-specific GDH of Streptomyces clavuligerus

A new class of glutamate dehydrogenase (GDH) is reported. The GDH of Streptomyces clavuligerus was purified to homogeneity and characterized. It has a native molecular mass of 1,100 kDa and exists as an alpha(6) oligomeric structure composed of 183-kDa subunits. GDH, which requires AMP as an essential activator, shows a maximal rate of catalysis in 100 mm phosphate buffer, pH 7.0, at 30 degrees C. Under these conditions, GDH displayed hyperbolic behavior toward ammonia (K(m), 33 mm) and sigmoidal responses to changes in alpha-ketoglutarate (S(0.5) 1.3 mm; n(H) 1.50) and NADH (S(0.5) 20 microm; n(H) 1.52) concentrations. Aspartate and asparagine were found to be allosteric activators. This enzyme is inhibited by an excess of NADH or NH(4)(+), by some tricarboxylic acid cycle intermediates and by ATP. This GDH seems to be a catabolic enzyme as indicated by the following: (i) it is NAD-specific; (ii) it shows a high value of K(m) for ammonia; and (iii) when S. clavuligerus was cultured in minimal medium containing glutamate as the sole source of carbon and nitrogen, a 5-fold increase in specific activity of GDH was detected compared with cultures provided with glycerol and ammonia. GDH has 1,651 amino acids, and it is encoded by a DNA fragment of 4,953 base pairs (gdh gene). It shows strong sequence similarity to proteins encoded by unidentified open reading frames present in the genomes of species belonging to the genera Mycobacterium, Rickettsia, Pseudomonas, Vibrio, Shewanella, and Caulobacter, suggesting that it has a broad distribution. The GDH of S. clavuligerus is the first member of a class of GDHs included in a subfamily of GDHs (large GDHs) whose catalytic requirements and evolutionary implications are described and discussed.

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.