Proton transfer function of carbonic anhydrase: Insights from QM/MM simulations

Recent QM/MM analyses of proton transfer function of human carbonic anhydrase II (CAII) are briefly reviewed. The topics include a preliminary analysis of nuclear quadrupole coupling constant calculations for the zinc ion and more detailed analyses of microscopic pK_a of the zinc-bound water and free energy profile for the proton transfer. From a methodological perspective, our results emphasize that performing sufficient sampling is essential to the calculation of all these quantities, which reflects the well solvated nature of CAII active site. From a mechanistic perspective, our analyses highlight the importance of electrostatics in shaping the energetics and kinetics of proton transfer in CAII for its function. We argue that once the pK_a for the zinc-bound water is modulated to be in the proper range (~ 7.0), proton transfer through a relatively well solvated cavity towards/from the protein surface (His64) does not require any major acceleration. Therefore, although structural details like the length of the water wire between the donor and acceptor groups still may make a non-negligible contribution, our computational results and the framework of analysis suggest that the significance of such “fine-tuning” is likely secondary to the modulation of pK_a of the zinc-bound water. We encourage further experimental analysis with mutation of (charged) residues not in the immediate neighborhood of the zinc ion to quantitatively test this electrostatics based framework; in particular, Φ analysis based on these mutations may shed further light into the relative importance of the classical Grotthus mechanism and the “proton hole” pathway that we have proposed recently for CAII.     

Structure and catalysis by carbonic anhydrase II: role of active-site tryptophan 5

The tryptophan residue Trp5, highly conserved in the α class of carbonic anhydrases including human carbonic anhydrase II (HCA II), is positioned at the entrance of the active site cavity and forms a π-stacking interaction with the imidazole ring of the proton shuttle His64 in its outward orientation. We have observed that replacement of Trp5 in HCA II caused significant structural changes, as determined by X-ray diffraction, in the conformation of 11 residues at the N-terminus and in the orientation of the proton shuttle residue His64. Most significantly, two variants W5H and W5E HCA II had His64 predominantly outward in orientation, while W5F and wild type showed the superposition of both outward and inward orientations in crystal structures. Although Trp5 influences the orientation of the proton shuttle His64, this orientation had no significant effect on the rate constant for proton transfer near 1 μs⁻¹, determined by exchange of ¹⁸O between CO₂ and water measured by mass spectrometry. The apparent values of the pK_a of the zinc-bound water and the proton shuttle residue suggest that different active-site conformations influence the two stages of catalysis, the proton transfer stage and the interconversion of CO₂ and bicarbonate.     

Structural studies of β-carbonic anhydrase from the green alga Coccomyxa: inhibitor complexes with anions and acetazolamide

The β-class carbonic anhydrases (β-CAs) are widely distributed among lower eukaryotes, prokaryotes, archaea, and plants. Like all CAs, the β-enzymes catalyze an important physiological reaction, namely the interconversion between carbon dioxide and bicarbonate. In plants the enzyme plays an important role in carbon fixation and metabolism. To further explore the structure-function relationship of β-CA, we have determined the crystal structures of the photoautotroph unicellular green alga Coccomyxa β-CA in complex with five different inhibitors: acetazolamide, thiocyanate, azide, iodide, and phosphate ions. The tetrameric Coccomyxa β-CA structure is similar to other β-CAs but it has a 15 amino acid extension in the C-terminal end, which stabilizes the tetramer by strengthening the interface. Four of the five inhibitors bind in a manner similar to what is found in complexes with α-type CAs. Iodide ions, however, make contact to the zinc ion via a zinc-bound water molecule or hydroxide ion — a type of binding mode not previously observed in any CA. Binding of inhibitors to Coccomyxa β-CA is mediated by side-chain movements of the conserved residue Tyr-88, extending the width of the active site cavity with 1.5-1.8 Å. Structural analysis and comparisons with other α- and β-class members suggest a catalytic mechanism in which the movements of Tyr-88 are important for the CO₂-HCO₃ ^- interconversion, whereas a structurally conserved water molecule that bridges residues Tyr-88 and Gln-38, seems important for proton transfer, linking water molecules from the zinc-bound water to His-92 and buffer molecules.     

The Important Role of Active Site Water in the Catalytic Mechanism of Human Carbonic Anhydrase II – A Semiempirical MO Approach to the Hydration of CO2 †

The approach of CO₂ to a series of active site model complexes of human carbonic anhydrase II (HCAII) and its catalytic hydration to bicarbonate anion have been investigated using semiempirical MO theory (AM1). The results show that direct nucleophilic attack of zinc-bound hydroxide to the substrate carbon occurs in each model system. Further rearrangement of the bicarbonate complex thus formed via a rotation-like movement of the bicarbonate ligand can only be found in active site model systems that include at least one additional water molecule. Further refinement of the model complex by adding a methanol molecule to mimic Thr-199 makes this process almost activationless. The formation of the final bicarbonate complex by an internal (intramolecular) proton transfer is only possible in the simplest of all model systems, namely {[Im₃Zn(OH)]⁺·CO₂}. The energy of activation for this process, however, is 36.8 kcal·mol^–1 and thus too high for enzymatic catalysis. Therefore, we conclude that within the limitations of the model systems presented and the level of theory employed, the overall mechanism for the formation of the bicarbonate complex comprises an initial direct nucleophilic attack of zinc-bound hydroxide to carbon dioxide followed by a rotation-like rearrangement of the bicarbonate ligand via a penta-coordinate Zn²⁺ transition state structure, including the participation of an extra active site water molecule.Electronic Supplementary Material available.     

Proton transfer to residues of basic pK(a) during catalysis by carbonic anhydrase.

The maximal velocity in the hydration of CO(2) catalyzed by the carbonic anhydrases in well-buffered solutions is limited by an intramolecular proton transfer from zinc-bound water to acceptor groups of the enzyme and hence to buffer in solution. Stopped-flow spectrophotometry was used to accumulate evidence that this maximal velocity is affected by residues of basic pK(a), near 8 to above 9, in catalysis of the hydration of CO(2) by carbonic anhydrases III, IV, V, and VII. A mutant of carbonic anhydrase II containing the replacement His-64-->Ala, which removes the prominent histidine proton shuttle (with pK(a) near 7), allows better observation of these basic groups. We suggest this feature of catalysis is general for the human and animal carbonic anhydrases and is due to residues of basic pK(a), predominantly lysines and tyrosines more distant from the zinc than His-64, that act as proton acceptors. These groups supplement the well-studied proton transfer from zinc-bound water to His-64 in the most efficient of the carbonic anhydrases, isozymes II, IV, and VII.     

Role of hydrophilic residues in proton transfer during catalysis by human carbonic anhydrase II

Catalysis by the zinc metalloenzyme human carbonic anhydrase II (HCA II) is limited in maximal velocity by proton transfer between His64 and the zinc-bound solvent molecule. Asn62 extends into the active site cavity of HCA II adjacent to His64 and has been shown to be one of several hydrophilic residues participating in a hydrogen-bonded solvent network within the active site. We compared several site-specific mutants of HCA II with replacements at position 62 (Ala, Val, Leu, Thr, and Asp). The efficiency of catalysis in the hydration of CO 2 for the resulting mutants has been characterized by (18)O exchange, and the structures of the mutants have been determined by X-ray crystallography to 1.5-1.7 A resolution. Each of these mutants maintained the ordered water structure observed by X-ray crystallography in the active site cavity of wild-type HCA II; hence, this water structure was not a variable in comparing with wild type the activities of mutants at residue 62. Crystal structures of wild-type and N62T HCA II showed both an inward and outward orientation of the side chain of His64; however, other mutants in this study showed predominantly inward (N62A, N62V, N62L) or predominantly outward (N62D) orientations of His64. A significant role of Asn62 in HCA II is to permit two conformations of the side chain of His64, the inward and outward, that contributes to maximal efficiency of proton transfer between the active site and solution. The site-specific mutant N62D had a mainly outward orientation of His64, yet the difference in p K a between the proton donor His64 and zinc-bound hydroxide was near zero, as in wild-type HCA II. The rate of proton transfer in catalysis by N62D HCA II was 5% that of wild type, showing that His64 mainly in the outward orientation is associated with inefficient proton transfer compared with His64 in wild type which shows both inward and outward orientations. These results emphasize the roles of the residues of the hydrophilic side of the active site cavity in maintaining efficient catalysis by carbonic anhydrase.     

Human carbonic anhydrase III: structural and kinetic study of catalysis and proton transfer

The residue phenylalanine 198 (Phe 198) is a prominent cause of the lower activity of human carbonic anhydrase III (HCA III) compared with HCA II and other isozymes which have leucine at this site. We report the crystal structures of HCA III and the site-directed mutant F198L HCA III, both at 2.1 A resolution, and the enhancement of catalytic activity by exogenous proton donors containing imidazole rings. Both enzymes had a hexahistidine extension at the carboxy-terminal end, used to aid in purification, that was ordered in the crystal structures bound in the active site cavity of an adjacent symmetry-related enzyme. This observation allowed us to comment on a number of possible binding sites for imidazole and derivatives as exogenous proton donors/acceptors in catalysis by HCA III. Kinetic and structural evidence indicates that the phenyl side chain of Phe 198 in HCA III, about 5 A from the zinc, is a steric constriction in the active site, may cause altered interactions at the zinc-bound solvent, and is a binding site for the activation of catalysis by histidylhistidine. This suggests that sites of activation of the proton-transfer pathway in carbonic anhydrase are closer to the zinc than considered in previous studies.     

Proton transfer pathways in the mutant His-64-Ala of human carbonic anhydrase II

We have investigated the possible proton transfer pathways from the surface of the protein to the zinc-bound water molecule in the mutant His-64-Ala of human carbonic anhydrase II. Starting with an input of known crystallographic structures of the mutant, we model the proton pathways as hydrogen-bonded networks of proton conducting groups and bound solvent molecules. No proton path is detected in the mutant, in close agreement with the experimental observation of a 20-fold decrease in its catalytic efficiency compared to the wild-type enzyme. We also investigate in detail changes in hydration structure at the active site of the mutant and the resulting proton paths in the presence of an exogenous proton donor 4-methylimidazole (4-MI). The proton transfer pathways thus detected are correlated to the observed chemical rescue of catalytic activity by 4-MI.     

Molecular dynamics simulations of a branched tetradecasaccharide substrate in the active site of a xyloglucan endo-transglycosylase

Molecular dynamics simulations of the tetradecasaccharide XXXGXXXG in complex with the hybrid aspen xyloglucan endo-transglycosylase PttXET16-34 have been performed and analysed with respect to structure, dynamics, flexibility and ligand interactions. Notably, the charge state of the so-called ‘helper residue’ aspartate 87 (Asp87), which lies between the catalytic nucleophile [glutamate 85 (Glu85)] and general acid/base (Glu89) residues on the same beta strand, had a significant effect on PttXET16-34 active site structure. When Asp87 was deprotonated, electrostatic repulsion forced the nucleophile away from C1 of the sugar ring in subsite ? 1 and the proton–donating ability of Glu89 was also weakened due to the formation of a hydrogen bond with Asp87, whereas the protonation of Asp87 resulted in the formation of a hydrogen bond with the catalytic nucleophile and correct positioning of the catalytic machinery. The results suggest that catalysis in glycoside hydrolase family 16, and by extension clan GH-B enzymes, is optimal when the catalytic nucleophile is deprotonated for nucleophilic attack on the substrate, whereas the ‘helper residue’ and general acid/base residue are both in their conjugate acid forms to align the nucleophile and deliver a proton to the departing sugar, respectively.     

Kinetic analysis of multiple proton shuttles in the active site of human carbonic anhydrase

We have prepared a site-specific mutant of human carbonic anhydrase (HCA) II with histidine residues at positions 7 and 64 in the active site cavity. Using a different isozyme, we have placed histidine residues in HCA III at positions 64 and 67 and in another mutant at positions 64 and 7. Each of these histidine residues can act as a proton transfer group in catalysis when it is the only nonliganding histidine in the active site cavity, except His(7) in HCA III. Using an (18)O exchange method to measure rate constants for intramolecular proton transfer, we have found that inserting two histidine residues into the active site cavity of either isozyme II or III of carbonic anhydrase results in rates of proton transfer to the zinc-bound hydroxide that are antagonistic or suppressive with respect to the corresponding single mutants. The crystal structure of Y7H HCA II, which contains both His(7) and His(64) within the active site cavity, shows the conformation of the side chain of His(64) moved from its position in the wild type and hydrogen-bonded through an intervening water molecule with the side chain of His(7). This suggests a cause of decreased proton transfer in catalysis.     

Hydration of CO2 by carbonic anhydrase: intramolecular proton transfer between Zn2+-bound H2O and histidine 64 in human carbonic anhydrase II

The energy barrier for the intramolecular proton transfer between zinc-bound water and His 64 in the active site of human carbonic anhydrase II (HCA II) has been studied at the partial retention of diatomic differential overlap (PRDDO) level. The most important stabilizing factor for the intramolecular proton transfer is the zinc ion, which lowers the pKa of zinc-bound water and electrostatically repels the proton. The energy barrier of 127.5 kcal/mol for proton transfer between a water dimer is completely removed in the presence of the zinc ion. The zinc ligands, which donate electrons to the zinc ion, raise the barrier slightly to 34 kcal/mol for a 4-coordinated zinc complex including three imidazole ligands from His 94, His 96, and His 119 and to 54 kcal/mol for the 5-coordinated zinc complex including the fifth water ligand. A few model calculations indicate that these energy barriers are expected to be reduced to within experimental range (approximately 10 kcal/mol) when large basis set, correlation energies, and molecular dynamics are considered. The proton-transfer group, which functions as proton receiver in the intramolecular proton transfer, helps to attract the proton; and the partially ordered active site water molecules are important for proton relay function.     

Catalytic sites of medicinal leech enzyme destabilase-lysozyme (mlDL): Structure-function relationship

Based on the three-dimensional model of the bifunctional enzyme destabilase-lysozyme of the medicinal leech (mlDL) in complex with trimer of N-acetylglucosamine (NAG)₃ by site-directed mutagenesis method, the functional role of the group of amino acids (Glu14, Asp26, Ser29, Ser31, Lys38, His92) in manifestation of lysozyme (glycosidase, muramidase) and isopeptidase activities has been investigated by site-directed mutagenesis. The results obtained go well with hypothesis, that lysozyme active site of mlDL includes catalytic Glu14 and Asp26 residues, and isopeptidase site functions as Ser/Lys catalytic dyad presented by catalytic residues Ser29 and Lys38. Thus, among the invertebrate lysozymes, mlDL presents the first example of a bifunctional enzyme with identified position of the isopeptidase active site and localization of the corresponding catalytic residues.     

Phosphophoryn,a biomineralization template protein: pH-dependent protein folding experiments

The protein folding behavior of a polyelectrolyte protein, bovine dentine phosphophoryn (BDPP), in the pH range of 1.82–11.0 has been investigated. One- and two-dimensional nmr spectroscopy has been utilized to obtain proton spin assignments for amino acid residues in D₂O and in H₂O. One-dimensional ³¹P-nmr experiments verify the existence of three separate classes of O-phosphoserine (PSer) resonances in BDPP (α, β, χ), representing three distinct PSer residue populations at pH 6.94. By means of pH titration and ¹H-nmr, five populations of Asp residues can be identified. Three of these populations exhibit secondary inflection points on their pH titration curves that correspond to an observed pK_a of 6.17–6.95. The presence or absence of secondary inflection points for Asp populations and the ³¹P-nmr chemical shift dispersion for the three PSer residue populations indicate that BDPP may be comprised of homologous (Asp-Asp)_n. (PSer-PSer)_n, and heterologous (PSer-Asp)_n sequences arranged into polyelectrolyte cluster regions. The pH titration also revealed that certain populations of Ser, Gly, and Pro residues in BDPP exhibit pH-dependent resonance frequency shifts. The “apparent” pK, for the transition points of these frequency shifts corresponds to either the pK of Pser monophosphatc ester and/or the pK_a of Asp COOH group of BDPP polyelectrolyte regions. On the basis of these transition points, we can assign four types of Ser, Gly, or Pro-containing “intervening” regions in BDPP, based on their sensitivity to protonation and deprotonation events occurring at (Asp)_n, (PSer)_n, or (PSer-Asp)_n anionic cluster regions that flank the intervening regions. Our ¹H-ninr experiments also reveal that BDPP assumes a folded conformation at low pH. As the pH increases, this conformation undergoes several unfolding transitions as the BDPP molecule assumes more open conformations in response to increased electrostatic repulsion between polyelectrolyte anionic regions in the protein. These folding-unfolding transitions are mediated by the intervening regions, which act as “hinges” to allow the polyelectrolyte regions to fold relative to one another. © 1994 John Wiley & Sons, Inc.     

Proton transfer in catalysis and the role of proton shuttles in carbonic anhydrase

The undisputed role of His64 in proton transfer during catalysis by carbonic anhydrases in the α class has raised questions concerning the details of its mechanism. The highly conserved residues Tyr7, Asn62, and Asn67 in the active-site cavity function to fine tune the properties of proton transfer by human carbonic anhydrase II (HCA II). For example, hydrophobic residues at these positions favor an inward orientation of His64 and a low pK_a for its imidazole side chain. It appears that the predominant manner in which this fine tuning is achieved in rate constants for proton transfer is through the difference in pK_a between His64 and the zinc-bound solvent molecule. Other properties of the active-site cavity, such as inward and outward conformers of His64, appear associated with the change in ΔpK_a; however, there is no strong evidence to date that the inward and outward orientations of His64 are in themselves requirements for facile proton transfer in carbonic anhydrase.     

Marcus rate theory applied to enzymatic proton transfer

The hydration of CO(2) and the dehydration of HCO(3)(-) catalyzed by the carbonic anhydrases is accompanied by the transfer of protons between solution and the zinc-bound water molecule in the active site. This transfer is facilitated by amino acid residues of the enzyme which act as intramolecular proton shuttles; variants of carbonic anhydrase lacking such shuttle residues are enhanced or rescued in catalysis by intermolecular proton transfer from donors such as imidazole in solution. The resulting rate constants for proton transfer when compared with the values of the pK(a) of the donor and acceptor give Bronsted plots of high curvature. These data are described by Marcus theory which shows an intrinsic barrier for proton transfer from 1 to 2 kcal/mol and work terms or thermodynamic contributions to the free energy of reaction from 4 to10 kcal/mol. The interpretation of these Marcus parameters is discussed in terms of the well-studied pathway of the catalysis and structure of the enzymes.     

Shifting the optimum pH of Bacillus circulans xylanase towards acidic side by introducing arginine

Electrostatic interactions are important in protein folding, binding, flexibility, stability and function. The pH at which the enzyme is maximally active is determined by the pK_as of the active site residues, which are modulated by several factors including the change in electrostatics in its vicinity. As the acidic xylanases are important in food and animal feed industries, electrostatic interactions are introduced in Bacillus circulans xylanase to shift their pH optima towards the acidic side. Arg substitutions are made to modulate the pKas of the active site residues. Neutral residues are substituted by Arg in such a way that the substituted residue can make direct interaction with the catalytic residues. However, the mutations with other titratable residues (Asp, Arg, Lys, His, Tyr, and Ser) present in between the catalytic sites and the substituted sites are avoided. Site directed mutagenesis was conducted to confirm the strategy. The results show the shift in pH optima of the mutants towards the acidic side by 0.5–1.5 unit. Molecular dynamics simulation of the mutant V37R reveals that the decrease in activity is due to the increase in distance between the substrate oxygen atoms and catalytic glutamates.     

Activation of β- and γ-carbonic anhydrases from pathogenic bacteria with tripeptides

Six tripeptides incorporating acidic amino acid residues were prepared for investigation as activators of β- and γ-carbonic anhydrases (CAs, EC 4.2.1.1) from the pathogenic bacteria Vibrio cholerae, Mycobacterium tuberculosis, and Burkholderia pseudomallei. The primary amino acid residues that are involved in the catalytic mechanisms of these CA classes are poorly understood, although glutamic acid residues near the active site appear to be involved. The tripeptides that contain Glu or Asp residues can effectively activate VchCAβ and VchCAγ (enzymes from V. cholerae), Rv3273 CA (mtCA3, a β-CA from M. tuberculosis) and BpsCAγ (γ-CA from B. pseudomallei) at 0.21–18.1?µM levels. The position of the acidic residues in the peptide sequences can significantly affect bioactivity. For three of the enzymes, tripeptides were identified that are more effective activators than both l-Glu and l-Asp. The tripeptides are also relatively selective because they do not activate prototypical α-CAs (human carbonic anhydrases I and II). Because the role of CA activators in the pathogenicity and life cycles of these infectious bacteria are poorly understood, this study provides new molecular probes to explore such processes.     

Buffer enhancement of proton transfer in catalysis by human carbonic anhydrase III

Among the isozymes of carbonic anhydrase, isozyme III is the least efficient in the catalysis of the hydration of CO2 and was previously thought to be unaffected by proton transfer from buffers to the active site. We report that buffers of small size, especially imidazole, increase the rate of catalysis by human carbonic anhydrase III (HCA III) of (1) 18O exchange between HCO3- and water measured by membrane-inlet mass spectrometry and (2) the dehydration of HCO3- measured by stopped-flow spectrophotometry. Imidazole enhanced the rate of release of 18O-labeled water from the active site of wild-type carbonic anhydrase III and caused a much greater enhancement, up to 20-fold, for the K64H, R67H, and R67N mutants of this isozyme. Imidazole had no effect on the rate of interconversion of CO2 and HCO3- at chemical equilibrium. Steady-state measurements showed that the addition of imidazole resulted in increases in the turnover number (kcat) for the hydration of CO2 catalyzed by HCA III and for the dehydration of HCO3- catalyzed by R67N HCA III. These results are consistent with the transfer of a proton from the imidazolium cation to the zinc-bound hydroxide at the active site, a step required to regenerate the active form of enzyme in the catalytic cycle. Like isozyme II of carbonic anhydrase, isozyme III can be enhanced in catalytic rate by the presence of small molecule buffers in solution.     

Proton transfer within the active-site cavity of carbonic anhydrase III

The maximal turnover rate of CO2 hydration catalyzed by the carbonic anhydrases is limited by proton transfer steps from the zinc-bound water to solution, steps that regenerate the catalytically active zinc-bound hydroxide. Catalysis of CO2 hydration by wild-type human carbonic anhydrase III (HCA III) (k(cat) = 2 ms (-1)) is the least efficient among the carbonic anhydrases in its class, in part because it lacks an efficient proton shuttle residue. We have used site-directed mutagenesis to test positions within the active-site cavity of HCA III for their ability to carry out proton transfer by replacing various residues with histidine. Catalysis by wild-type HCA III and these six variants was determined from the initial velocity of hydration of CO2 measured by stopped-flow spectrophotometry and from the exchange of 18O between CO2 and H2O at chemical equilibrium by mass spectrometry. The results show that histidine at three positions (Lys64His, Arg67His and Phe131His) have the capacity to transfer protons during catalysis, enhancing maximal velocity of CO2 hydration and 18O exchange from 4- to 15-fold compared with wild-type HCA III. Histidine residues at the other three positions (Trp5His, Tyr7His, Phe20His) showed no firm evidence for proton transfer. These results are discussed in terms of the stereochemistry of the active-site cavity and possible proton transfer pathways.     

A proposal for the mechanism of carbonic anhydrase action

A new catalytic mechanism is proposed for the hydration of CO₂ by the zinc metalloenzyme carbonic anhydrase. This mechanism identifies the group controlling catalytic activity as an active site histidine, in which the protonated imidazole ring coordinates to zinc, losing a proton. Geometric constraints on the histidine unit make the metal-ligand bond a strained and, therefore, labile one. In the hydration reaction, the metal-bound neutral histidine moiety serves as a proton acceptor for the transient ionization of metal-bound water. Zinc-bound hydroxide attacks the carbon of the substrate to generate metal-bound bicarbonate, and the system regenerates itself by losing the elements of carbonic acid.