A closer look at the active site of gamma-class carbonic anhydrases: high-resolution crystallographic studies of the carbonic anhydrase from Methanosarcina thermophila

The prototype of the gamma-class of carbonic anhydrase has been characterized from the methanogenic archaeon Methanosarcina thermophila. Previously reported kinetic studies of the gamma-class carbonic anhydrase are consistent with this enzyme having a reaction mechanism similar to that of the mammalian alpha-class carbonic anhydrase. However, the overall folds of these two enzymes are dissimilar, and apart from the zinc-coordinating histidines, the active site residues bear little resemblance to one another. The crystal structures of zinc-containing and cobalt-substituted gamma-class carbonic anhydrases from M. thermophila are reported here between 1.46 and 1.95 A resolution in the unbound form and cocrystallized with either SO(4)(2)(-) or HCO(3)(-). Relative to the tetrahedral coordination geometry seen at the active site in the alpha-class of carbonic anhydrases, the active site of the gamma-class enzyme contains additional metal-bound water ligands, so the overall coordination geometry is trigonal bipyramidal for the zinc-containing enzyme and octahedral for the cobalt-substituted enzyme. Ligands bound to the active site all make contacts with the side chain of Glu 62 in manners that suggest the side chain is likely protonated. In the uncomplexed zinc-containing enzyme, the side chains of Glu 62 and Glu 84 appear to share a proton; additionally, Glu 84 exhibits multiple conformations. This suggests that Glu 84 may act as a proton shuttle, which is an important aspect of the reaction mechanism of alpha-class carbonic anhydrases. A hydrophobic pocket on the surface of the enzyme may participate in the trapping of CO(2) at the active site. On the basis of the coordination geometry at the active site, ligand binding modes, the behavior of the side chains of Glu 62 and Glu 84, and analogies to the well-characterized alpha-class of carbonic anhydrases, a more-defined reaction mechanism is proposed for the gamma-class of carbonic anhydrases.     

Characterization of heterologously produced carbonic anhydrase from Methanosarcina thermophila.

The gene encoding carbonic anhydrase from Methanosarcina thermophila was hyperexpressed in Escherichia coli, and the heterologously produced enzyme was purified 14-fold to apparent homogeneity. The enzyme purified from E. coli has properties (specific activity, inhibitor sensitivity, and thermostability) similar to those of the authentic enzyme isolated from M. thermophila; however, a discrepancy in molecular mass suggests that the carbonic anhydrase is posttranslationally modified in either E. coli or M. thermophila. Both the authentic and heterologously produced enzymes were stable to heating at 55 degrees C for 15 min but were inactivated at higher temperatures. No esterase activity was detected with p-nitrophenylacetate as the substrate. Plasma emission spectroscopy revealed approximately 0.6 Zn per subunit. As judged from the estimated native molecular mass, the enzyme is either a trimer or a tetramer. Western blot (immunoblot) analysis of cell extract proteins from M. thermophila indicates that the levels of carbonic anhydrase are regulated in response to the growth substrate, with protein levels higher in acetate than in methanol- or trimethylamine-grown cells.     

Carbonic anhydrase inhibitors. Inhibition of the zinc and cobalt gamma-class enzyme from the archaeon Methanosarcina thermophila with anions

Anions represent the second class of inhibitors of the zinc enzyme carbonic anhydrase (CA, EC 4.2.1.1), in addition to sulfonamides, which possess clinical applications. The first inhibition study of the zinc and cobalt gamma-class enzyme from the archaeon Methanosarcina thermophila (Cam) with anions is reported here. Inhibition data of the alpha-class human isozymes hCA I and hCA II (cytosolic) as well as the membrane-bound isozyme hCA IV with a large number of anionic species such as halides, pseudohalides, bicarbonate, carbonate, nitrate, nitrite, hydrosulfide, bisulfite, and sulfate, etc., are also provided for comparison. The best Zn-Cam anion inhibitors were hydrogen sulfide and cyanate, with inhibition constants in the range of 50-90 microM, whereas thiocyanate, azide, carbonate, nitrite, and bisulfite were weaker inhibitors (K(I)s in the range of 5.8-11.7 mM). Fluoride, chloride, and sulfate do not inhibit this enzyme appreciably up to concentrations of 200 mM, whereas the substrate bicarbonate behaves as a weak inhibitor (K(I)s of 42 mM). The best Co-Cam inhibitor was carbonate, with an inhibition constant of 9 microM, followed by nitrate and bicarbonate (K(I)s in the range of 90-100 microM). The metal poisons were much more ineffective inhibitors of this enzyme, with cyanide possessing an inhibition constant of 51.5mM, whereas cyanate, thiocyanate, azide, iodide, and hydrogen sulfide showed K(I)s in the range of 2.0-6.1mM. As for Zn-Cam, fluoride, chloride, and sulfate are not inhibitors of Co-Cam. These major differences between the two gamma-CAs investigated here can be explained only in part by the different geometries of the metal ions present within their active sites.     

Evolution of structure in gamma-class carbonic anhydrase and structurally related proteins

Protein structure contains evolutionary information and it is more highly conserved than sequence. The evolution of structure in gamma-class carbonic anhydrase (gamma-CA) and its structurally related proteins (gammaCASRPs) were discussed. To obtain a reliable analysis, we defined a subset that contains all specificities and organisms as the nonredundant set using QR factorization based on the multiple structural alignment of the known crystallographic structures of gammaCASRPs with Q(H) as the structural homology measure. Then, we applied unweighted pair group method with arithmetic averages (UPGMA) to reconstruct structural phylogeny. We found that gamma-CA most likely arose through duplication events; the domain of gamma-CA underwent a process of alpha-helical content from amino-terminal end to carboxyl-terminal end of the left-handed beta-helix (LbetaH); the capacity of gamma-CA to bind Zn occurred early in evolution and only later included the ability to catalyze the reversible hydration of CO(2) efficiently for the occurrence of two loops involving Glu 62 and Glu 84, respectively, and a long helix at the carboxyl-terminal end of the LbetaH. In addition, the main conserved regions in these structures are in the structurally constrained residues of LbetaH domain, and the topology of the structural dendrogram can be rather easily understood in terms of functional diversification.     

Proposal for a hydrogen bond network in the active site of the prototypic gamma-class carbonic anhydrase

The crystal structure of Cam, the prototypic gamma-class carbonic anhydrase, reveals active site residues Gln75, Asn73, and Asn 202 previously hypothesized to participate in catalysis. These potential roles were investigated for the first time by kinetic analyses of site-specific replacement variants of the zinc and cobalt forms of Cam. Gln75 replacement variants showed large decreases in k(cat)/K(m) relative to wild-type. Further, the Gln75 variants showed a loss of the pK(a) in pH versus k(cat)/K(m) profiles previously attributed to ionization of the metal-bound water yielding the hydroxyl group attacking CO(2). These results support the previously proposed role for Gln75 in hydrogen bonding with the catalytic hydroxyl orienting it for attack on CO(2). Kinetic analyses of Asn73 variants were consistent with a role in hydrogen bonding with Gln75 to position it for optimal interaction with the catalytic hydroxyl. Kinetic analyses of Asn202 variants showed substantial decreases in k(cat)/K(m) relative to the wild-type enzyme supporting the previously hypothesized role in polarizing CO(2) and facilitating attack from the metal-bound hydroxyl. On the basis of results presented here, and previously reported structural analyses, we present a catalytic mechanism involving Gln75, Asn73, and Asn202 that also suggests a role for Glu62 not previously recognized. Finally, the results suggest that the gamma-, beta-, and alpha-class carbonic anhydrases each independently evolved variations of a fundamental hydrogen bond network essential for catalysis.     

Kinetic and spectroscopic characterization of the gamma-carbonic anhydrase from the methanoarchaeon Methanosarcina thermophila.

The zinc and cobalt forms of the prototypic gamma-carbonic anhydrase from Methanosarcina thermophila were characterized by extended X-ray absorption fine structure (EXAFS) and the kinetics were investigated using steady-state spectrophotometric and (18)O exchange equilibrium assays. EXAFS results indicate that cobalt isomorphously replaces zinc and that the metals coordinate three histidines and two or three water molecules. The efficiency of either Zn-Cam or Co-Cam for CO(2) hydration (k(cat)/K(m)) was severalfold greater than HCO(3-) dehydration at physiological pH values, a result consistent with the proposed physiological function for Cam during growth on acetate. For both Zn- and Co-Cam, the steady-state parameter k(cat) for CO(2) hydration was pH-dependent with a pK(a) of 6.5-6.8, whereas k(cat)/K(m) was dependent on two ionizations with pK(a) values of 6.7-6.9 and 8.2-8.4. The (18)O exchange assay also identified two ionizable groups in the pH profile of k(cat)/K(m) with apparent pK(a) values of 6.0 and 8.1. The steady-state parameter k(cat) (CO(2) hydration) is buffer-dependent in a saturable manner at pH 8. 2, and the kinetic analysis suggested a ping-pong mechanism in which buffer is the second substrate. The calculated rate constant for intermolecular proton transfer is 3 x 10(7) M(-1) s(-1). At saturating buffer concentrations and pH 8.5, k(cat) is 2.6-fold higher in H(2)O than in D(2)O, suggesting that an intramolecular proton transfer step is at least partially rate-determining. At high pH (pH > 8), k(cat)/K(m) is not dependent on buffer and no solvent hydrogen isotope effect was observed, consistent with a zinc hydroxide mechanism. Therefore, at high pH the catalytic mechanism of Cam appears to resemble that of human CAII, despite significant structural differences in the active sites of these two unrelated enzymes.     

Chemical rescue of proton transfer in catalysis by carbonic anhydrases in the beta- and gamma-class

Catalysis of the dehydration of HCO(3)(-) by carbonic anhydrase requires proton transfer from solution to the zinc-bound hydroxide. Carbonic anhydrases in each of the alpha, beta, and gamma classes, examples of convergent evolution, appear to have a side chain extending into the active site cavity that acts as a proton shuttle to facilitate this proton transfer, with His 64 being the most prominent example in the alpha class. We have investigated chemical rescue of mutants in two of these classes in which a proton shuttle has been replaced with a residue that does not transfer protons: H216N carbonic anhydrase from Arabidopsis thaliana (beta class) and E84A carbonic anhydrase from the archeon Methanosarcina thermophila (gamma class). A series of structurally homologous imidazole and pyridine buffers were used as proton acceptors in the activation of CO(2) hydration at steady state and as proton donors of the exchange of (18)O between CO(2) and water at chemical equilibrium. Free energy plots of the rate constants for this intermolecular proton transfer as a function of the difference in pK(a) of donor and acceptor showed extensive curvature, indicating a small intrinsic kinetic barrier for the proton transfers. Application of Marcus rate theory allowed quantitative estimates of the intrinsic kinetic barrier which were near 0.3 kcal/mol with work functions in the range of 7-11 kcal/mol for mutants in the beta and gamma class, similar to results obtained for mutants of carbonic anhydrase in the alpha class. The low values of the intrinsic kinetic barrier for all three classes of carbonic anhydrase reflect proton transfer processes that are consistent with a model of very rapid proton transfer through a flexible matrix of hydrogen-bonded solvent structures sequestered within the active sites of the carbonic anhydrases.     

Proton transport in carbonic anhydrase: Insights from molecular simulation

This article reviews the insights gained from molecular simulations of human carbonic anhydrase II (HCA II) utilizing non-reactive and reactive force fields. The simulations with a reactive force field explore protein transfer and transport via Grotthuss shuttling, while the non-reactive simulations probe the larger conformational dynamics that underpin the various contributions to the rate-limiting proton transfer event. Specific attention is given to the orientational stability of the His64 group and the characteristics of the active site water cluster, in an effort to determine both of their impact on the maximal catalytic rate. The explicit proton transfer and transport events are described by the multistate empirical valence bond (MS-EVB) method, as are alternative pathways for the excess proton charge defect to enter/leave the active site. The simulation results are interpreted in light of experimental results on the wild-type enzyme and various site-specific mutations of HCA II in order to better elucidate the key factors that contribute to its exceptional efficiency.     

Characterization of the intramolecular electron transfer pathway from 2-hydroxyphenazine to the heterodisulfide reductase from Methanosarcina thermophila

Heterodisulfide reductase (HDR) is a component of the energy-conserving electron transfer system in methanogens. HDR catalyzes the two-electron reduction of coenzyme B-S-S-coenzyme M (CoB-S-S-CoM), the heterodisulfide product of the methyl-CoM reductase reaction, to free thiols, HS-CoB and HS-CoM. HDR from Methanosarcina thermophila contains two b-hemes and two [Fe(4)S(4)] clusters. The physiological electron donor for HDR appears to be methanophenazine (MPhen), a membrane-bound cofactor, which can be replaced by a water-soluble analog, 2-hydroxyphenazine (HPhen). This report describes the electron transfer pathway from reduced HPhen (HPhenH(2)) to CoB-S-S-CoM. Steady-state kinetic studies indicate a ping-pong mechanism for heterodisulfide reduction by HPhenH(2) with the following values: k(cat) = 74 s(-1) at 25 degrees C, K(m) (HPhenH(2)) = 92 microm, K(m) (CoB-S-S-CoM) = 144 microm. Rapid freeze-quench EPR and stopped-flow kinetic studies and inhibition experiments using CO and diphenylene iodonium indicate that only the low spin heme and the high potential FeS cluster are involved in CoB-S-S-CoM reduction by HPhenH(2). Fe-S cluster disruption by mersalyl acid inhibits heme reduction by HPhenH(2), suggesting that a 4Fe cluster is the initial electron acceptor from HPhenH(2). We propose the following electron transfer pathway: HPhenH(2) to the high potential 4Fe cluster, to the low potential heme, and finally, to CoB-S-S-CoM.     

Purification and characterization of a ferredoxin from acetate-grown Methanosarcina thermophila

A ferredoxin, which functions as an electron acceptor for the CO dehydrogenase complex from Methanosarcina thermophila, was purified from acetate-grown cells. It was isolated as a trimer having a native molecular weight of approximately 16,400 and monomer molecular weight of 4,888 calculated from the amino acid composition. The ferredoxin contained 2.80 +/- 0.56 Fe atoms and 1.98 +/- 0.12 acid-labile sulfide. UV-visible absorption maxima were 395 and 295 nm with monomeric extinction coefficients of epsilon 395 = 12,800 M-1 cm-1 and epsilon 295 = 14,460 M-1 cm-1. The A395/A295 ratio ranged from 0.80 to 0.88. There were 5 cysteines per monomer but no methionine, histidine, arginine, or aromatic amino acids. The N-terminal amino acid sequence showed a 4-cysteine cluster with potential to coordinate a Fe:S center. The protein was stable for 30 min at 70 degrees C, but denatured during incubation at 85 degrees C.     

A novel carbonic anhydrase from the giant clam Tridacna gigas contains two carbonic anhydrase domains

This report describes the presence of a unique dual domain carbonic anhydrase (CA) in the giant clam, Tridacna gigas. CA plays an important role in the movement of inorganic carbon (Ci) from the surrounding seawater to the symbiotic algae that are found within the clam's tissue. One of these isoforms is a glycoprotein which is significantly larger (70 kDa) than any previously reported from animals (generally between 28 and 52 kDa). This alpha-family CA contains two complete carbonic anhydrase domains within the one protein, accounting for its large size; dual domain CAs have previously only been reported from two algal species. The protein contains a leader sequence, an N-terminal CA domain and a C-terminal CA domain. The two CA domains have relatively little identity at the amino acid level (29%). The genomic sequence spans in excess of 17 kb and contains at least 12 introns and 13 exons. A number of these introns are in positions that are only found in the membrane attached/secreted CAs. This fact, along with phylogenetic analysis, suggests that this protein represents the second example of a membrane attached invertebrate CA and it contains a dual domain structure unique amongst all animal CAs characterized to date.     

Steady-state kinetic analysis of phosphotransacetylase from Methanosarcina thermophila

Phosphotransacetylase (EC 2.3.1.8) catalyzes the reversible transfer of the acetyl group from acetyl phosphate to coenzyme A (CoA), forming acetyl-CoA and inorganic phosphate. A steady-state kinetic analysis of the phosphotransacetylase from Methanosarcina thermophila indicated that there is a ternary complex kinetic mechanism rather than a ping-pong kinetic mechanism. Additionally, inhibition patterns of products and a nonreactive substrate analog suggested that the substrates bind to the enzyme in a random order. Dynamic light scattering revealed that the enzyme is dimeric in solution.     

Ferredoxin requirement for electron transport from the carbon monoxide dehydrogenase complex to a membrane-bound hydrogenase in acetate-grown Methanosarcina thermophila

Cell extracts from acetate-grown Methanosarcina thermophila contained CO-oxidizing:H2-evolving activity 16-fold greater than extracts from methanol-grown cells. Following fractionation of cell extracts into soluble and membrane components, CO-dependent H2 evolution and CO-dependent methyl-coenzyme M methylreductase activities were only present in the soluble fraction, but addition of the membrane fraction enhanced both activities. A b-type cytochrome(s), present in the membrane fraction, was linked to a membrane-bound hydrogenase. CO-oxidizing:H2-evolving activity was reconstituted with: (i) CO dehydrogenase complex, (ii) a ferredoxin, and (iii) purified membranes with associated hydrogenase. The ferredoxin was a direct electron acceptor for the CO dehydrogenase complex. The ferredoxin also coupled CO oxidation by CO dehydrogenase complex to metronidazole reduction.     

A study of the histidine residues of human carbonic anhydrase B using 270 MHz proton magnetic resonance

Nine resonances in the 270 MHz proton magnetic resonance spectrum of human carbonic anhydrase B have been identified with imidazole C₍₂₎ protons of histidine residues, six of which are observed to titrate with pK_a values in the range 4.7 to 7.4. The behaviour of the nine resonances has been studied in the presence of the inhibitors, iodide, cyanide, acetate, hexacyanochromate, and imidazole. Measurements have also been made of the enzyme in its apo, cobalt, and mono-alkylated forms. Used in conjunction with the crystal structure, these results have enabled the tentative assignment of all nine resonances to particular histidine residues in the amino-acid sequence. Three of the active-site histidines at positions 64, 67, and 200 have low pK_a values and cannot be directly linked to the activity of the enzyme. However, the resonances assigned to the three metal-liganding histidines do exhibit changes on anion binding and with pH, which parallel changes in the esterase activity. These results are consistent with the model of an ionizable water molecule bound to the zinc ion.Linewidth measurements of the resonances of the histidine residues on the enzyme surface are used to estimate pseudo-first-order rate constants of the order of 4 × 10³ s^?1 for D⁺ exchange between imidazole N and solvent in the absence of buffer. These rates are observed to increase in the presence of small amounts of the buffers Tris and imidazole.     

Properties of intramolecular proton transfer in carbonic anhydrase III.

We investigated the efficiency of glutamic acid 64 and aspartic acid 64 as proton donors to the zinc-bound hydroxide in a series of site-specific mutants of human carbonic anhydrase III (HCA III). Rate constants for this intramolecular proton transfer, a step in the catalyzed dehydration of bicarbonate, were determined from the proton-transfer-dependent rates of release of H2 18O from the enzyme measured by mass spectrometry. The free energy plots representing these rate constants could be fit by the Marcus rate theory, resulting in an intrinsic barrier for the proton transfer of deltaG0++ = 2.2 +/- 0.5 kcal/mol, and a work function or thermodynamic contribution to the free energy of reaction wr = 10.8 +/- 0.1 kcal/mol. These values are very similar in magnitude to the Marcus parameters describing intramolecular proton transfer from His64 and His67 to the zinc-bound hydroxide in mutants of HCA III. That result and the equivalent efficiency of Glu64 and Asp64 as proton donors in the catalysis by CA III demonstrate a lack of specificity in proton transfer from these sites, which is indirect evidence of a number of proton conduction pathways through different structures of intervening water chains. The dominance of the thermodynamic contribution or work function for all of these proton transfers is consistent with the view that formation and breaking of hydrogen bonds in such water chains is a limiting factor for proton translocation.     

Crystal structure of phosphotransacetylase from the methanogenic archaeon Methanosarcina thermophila

Phosphotransacetylase (Pta) [EC 2.3.1.8] is ubiquitous in the carbon assimilation and energy-yielding pathways in anaerobic prokaryotes where it catalyzes the reversible transfer of the acetyl group from acetyl phosphate to CoA forming acetyl CoA and inorganic phosphate. The crystal structure of Pta from the methane-producing archaeon Methanosarcina thermophila, representing the first crystal structure of any Pta, was determined by multiwavelength anomalous diffraction at 2.7 A resolution. In solution and in the crystal, the enzyme forms a homodimer. Each monomer consists of two alpha/beta domains with a cleft along the domain boundary, which presumably contains the substrate binding sites. Comparison of the four monomers present in the asymmetric unit indicates substantial variations in the relative orientation of the two domains and the structure of the putative active site cleft. A search for structural homologs revealed the NADP(+)-dependent isocitrate and isopropylmalate dehydrogenases as the only homologs with a similar two-domain architecture.     

Identification of Essential Glutamates in the Acetate Kinase from Methanosarcina thermophila

Acetate kinase catalyzes the reversible phosphorylation of acetate (CH₃COO⁻ + ATPCH₃CO₂PO₃²⁻ + ADP). A mechanism which involves a covalent phosphoryl-enzyme intermediate has been proposed, and chemical modification studies of the enzyme from Escherichia coli indicate an unspecified glutamate residue is phosphorylated (J. A. Todhunter and D. L. Purich, Biochem. Biophys. Res. Commun. 60:273–280, 1974). Alignment of the amino acid sequences for the acetate kinases from E. coli (Bacteria domain), Methanosarcina thermophila (Archaea domain), and four other phylogenetically divergent microbes revealed high identity which included five glutamates. These glutamates were replaced in the M. thermophila enzyme to determine if any are essential for catalysis. The histidine-tagged altered enzymes were produced in E. coli and purified to electrophoretic homogeneity by metal affinity chromatography. Replacements of E384 resulted in either undetectable or extremely low kinase activity, suggesting E384 is essential for catalysis which supports the proposed mechanism. Replacement of E385 influenced the K_m values for acetate and ATP with only moderate decreases in k_cat, which suggests that this residue is involved in substrate binding but not catalysis. The unaltered acetate kinase was not inactivated by N-ethylmaleimide; however, replacement of E385 with cysteine conferred sensitivity to N-ethylmaleimide which was prevented by preincubation with acetate, acetyl phosphate, ATP, or ADP, suggesting that E385 is located near the active site. Replacement of E97 decreased the K_m value for acetate but not ATP, suggesting this residue is involved in binding acetate. Replacement of either E32 or E334 had no significant effects on the kinetic constants, which indicates that neither residue is essential for catalysis or significantly influences the binding of acetate or ATP.