Quantum chemistry analysis of the mechanism of action of proteolytic enzymes. III. Proton transport in serine proteases

The model system for the proton transfer on the amide atom of the substrate leaving group based on the existence of "charge relay system" in the serine type proteases was analysed by the CNDO/2 method. The unfitness of this model to explain the action mechanism of serine proteases was shown. The model system for proton transfer with the water molecule as the intermediate acceptor of the Ser-195 proton was suggested and analysed by the same method. The acylation activation barrier of this system was shown to localize on the stage of synchronous transfer of the Ser-195 alcoholic proton and the water molecule proton hydrogen bound to the His-57 N epsilon 2-atom on the water molecule oxygen atom and the N epsilon 2-atom, respectively. The protonation of substrate in the case of the model system with the water molecule as the intermediate acceptor of proton was demonstrated to begin before the completion of the tetrahedral intermediate substance and the protonated from of the tetrahedral intermediate was shown to form only. A hypothesis considering the role of this water molecule as the nucleophilic reagent on the deacylation stage is presented.     

Proton transport behavior through the influenza A M2 channel: insights from molecular simulation

The structural properties of the influenza A virus M2 transmembrane channel in dimyristoylphosphatidylcholine bilayer for each of the four protonation states of the proton-gating His-37 tetrad and their effects on proton transport for this low-pH activated, highly proton-selective channel are studied by classical molecular dynamics with the multistate empirical valence-bond (MS-EVB) methodology. The excess proton permeation free energy profile and maximum ion conductance calculated from the MS-EVB simulation data combined with the Poisson-Nernst-Planck theory indicates that the triply protonated His-37 state is the most likely open state via a significant side-chain conformational change of the His-37 tetrad. This proposed open state of M2 has a calculated proton permeation free energy barrier of 7 kcal/mol and a maximum conductance of 53 pS compared to the experimental value of 6 pS. By contrast, the maximum conductance for Na(+) is calculated to be four orders of magnitude lower, in reasonable agreement with the experimentally observed proton selectivity. The pH value to activate the channel opening is estimated to be 5.5 from dielectric continuum theory, which is also consistent with experimental results. This study further reveals that the Ala-29 residue region is the primary binding site for the antiflu drug amantadine (AMT), probably because that domain is relatively spacious and hydrophobic. The presence of AMT is calculated to reduce the proton conductance by 99.8% due to a significant dehydration penalty of the excess proton in the vicinity of the channel-bound AMT.     

Conformation of cobrotoxin in aqueous solution as studied by nuclear magnetic resonance.

The 270-MHz proton NMR spectra of cobrotoxin from Naja naja atra were observed in 2H2O solution. The pKa value (5.93) of His-32 is slightly lower than the pKa value (6.65) of the reference model of N-acetylhistidine methylamide, because of the electrostatic interaction with Arg-33 and Asp-31. The pKa value (5.3--5.4) of His-4 is appreciably low, because of the interaction with the positively charged guanidino group possibly of Arg-59. The hydrogen-deuterium exchange rates in 2H2O solution were measured of cobrotoxin and imidazole-bearing models. The second-order rate constants of N-acetylhistidine methylamide, N-acetylhistidine and imidazole acetic acid satisfy the Br?nsted relation. With reference to this Br?nsted relation, the imidazole ring of His-32 is confirmed to be exposed. The imidazole ring of His-4 is also exposed and the exchange rate is excessively promoted by the presence possibly of Arg-59 in the proximity. All the methyl proton resonances are assigned to amino-acid types, by conventional double-resonance method and more effectively by the spin-echo double-resonance method. Eight methyl proton resonances are identified as due to the gamma and/or delta-methyl groups of Val-46, Leu-1, Ile-50 and Ile-52 residues. The proximity of aromatic ring protons and methyl protons is elucidated by the analyses of nulcear Overhauser effect enhancements. The aromatic proton resonances of Trp-29 are affected by the ionizable groups of Asp-31, His-32 and Tyr-35. The methyl groups of Ile-50 are in the proximity to the aromatic ring of Trp-29 and the methyl groups of Ile-52 are in the proximity to Tyr-25. The highest-field methyl proton resonance is due to a threonine residue in the proximity to His-4. The appreciable temperature-dependent chemical shift of this methyl proton resonance suggests a temperature-dependent local conformational equilibrium around the His-4 residue of the first loop of the cobrotoxin molecule.     

QM/MM study of the active site of free papain and of the NMA-papain complex.

Hybrid quantum mechanical/molecular mechanical (QM/MM) calculations using restricted and unrestricted Hartree-Fock and B3LYP ab initio (QM) and Amber force field (MM), respectively, have been applied to study the catalytic site of papain in both free and substrate bonded forms. Ab initio geometry optimizations have been performed for the active site of papain and the N-methyl-acetamide (NMA)-papain complex within the molecular mechanical treatment of the protein environment. A covalent tetrahedral intermediate structure could be obtained only when the amide N atom of the substrate molecule was protonated through a proton transfer from the His-159 in the catalytic site. Our results support the previous assumption that a proton transfer from His-159 to the amide N atom of the substrate occurs prior to or concerted with the nucleophilic attack of the Cys-25 sulfur atom to the carbonyl group of the substrate. The electron correlation effect will reduce the proton transfer barrier. Therefore, this proton transfer can be easily observed in the B3LYP/6-31G* calculations. The HF/6-31G* method overestimates the reaction barrier against this proton transfer. The sulfur atom of Cys-25 and the imidazole ring of His-159 are found to be coplanar in the free form of the enzyme. However, the rotation of the imidazole ring of His-159 was observed during the formation of the tetrahedral intermediate. Without the papain environment, the coplanar thiolate-imidazolium ion pair RS-...ImH+ is much less stable than the neutral form of RSH....Im. Within the protein environment, however, the thiolate-imidazolium ion pair becomes more stable than its neutral form by 4.1 and 0.4 kcal/mol in HF/6-31G* and B3LYP/6-31G* calculations, respectively. The barrier of proton transfer from S-H group of Cys-25 to the imidazole ring of His-159 was reduced from 22.0 kcal/mol to 15.2 kcal/mol by the protein environment in HF/6-31G* calculations. This barrier is found to be much smaller (2.5 kcal/mol) in B3LYP/6-31G* calculations.     

Stereochemical aspects of peptide hydrolysis catalyzed by serine proteases of the chymotrypsin type

The steric course of peptide hydrolysis catalyzed by serine proteases has been studied on the basis of the available, extensive structural data and taking into account the stereoelectronic theory of Deslongchamps (Heterocycles, 7, 1271 (1977)). These studies allowed elucidation of the structure of intermediates, in particular of the tetrahedral intermediate, and of the main structural events taking place during catalysis. They reveal a difficulty inherent in the generally accepted mechanism of peptide hydrolysis: protonation of the leaving nitrogen in the configuration arising from nucleophilic attack of Ser-195 on the carbonyl carbon cannot take place internally from His-57. Two alternative mechanisms are discussed which are compatible with all implications of the stereoelectronic theory. The main features of the more probable mechanism are: (i) a conformational change allowing the imidazole ring of His-57 to occupy two distinct positions; in one position a proton is abstracted from O^γ of Ser-195, and in the other this proton is donated to the leaving nitrogen; (ii) a configurational change (inversion) of the pyramidal leaving nitrogen reorienting the lone-pair orbital developed during nucleophilic attack; in one orientation CO bond breaking, and in the other CN bond breaking, is allowed. This inversion process confers on the nitrogen the property of a switch controlling the breakdown of the tetrahedral intermediate.     

Proton transfer roles of lysine 64 and glutamic acid 64 replacing histidine 64 in the active site of human carbonic anhydrase II.

The CO2 hydration activities of cloned human carbonic anhydrase II (carbonate hydro-lyase, EC 4.2.1.1) and variants with Lys, Glu, Gln or Ala replacing His at sequence position 64 have been measured in a variety of different buffers in the pH range 6-9. The variants with Lys-64, Gln-64 and Ala-64 showed non-Michaelis-Menten behavior under some conditions, apparent substrate inhibition being prominent near pH 9. However, asymptotic Michaelis-Menten parameters could be estimated for the limit of low substrate concentrations. All variants show distinct buffer specificities, and imidazole derivatives, Ches and phosphate buffers yield higher kcat values that Bicine, Taps and Mops buffers under otherwise similar conditions. These results are interpreted in terms of different pathways for a rate-limiting proton transfer. In unmodified enzyme, the very high catalytic activity depends on His-64 functioning as an efficient proton transfer group, but this pathway is not available in the variants with Gln-64 and Ala-64. Imidazoles, Ches and phosphate are thought to participate in a metal center-to-buffer proton transfer pathway, whereas Bicine, Taps, Mops and Mes appear to lack this capacity, so that the rate-limiting proton transfer occurs in a metal center-to-bulk water pathway for these variants. The Lys-64 and Glu-64 variants give significantly higher kcat values in Taps, Mops and Mes buffers than the Ala-64 and Gln-64 variants. The pH dependencies of these kcat values are compatible with the hypothesis that Lys-64 and Glu-64 can function as proton transfer groups. Thus, at pH near 9, Lys-64 appears to be only 5-times less efficient than His-64, while Glu-64 is inefficient. At pH 6, Lys-64 is an inefficient proton transfer group, but Glu-64 is only 2-3-times less efficient than His-64. The data indicate that Lys-64 and Glu-64 have pKa values near 8 and below 6, respectively.     

Role of protein motions on proton transfer pathways in human carbonic anhydrase II

We report here a theoretical study on the formation of long-range proton transfer pathways in proteins due to side chain conformational fluctuations of amino acid residues and reorganization of interior hydration positions. The proton transfer pathways in such systems may be modeled as fluctuating hydrogen-bonded networks with both short- and long-lived connections between the networked nodes, the latter being formed by polar protein atoms and water molecules. It is known that these fluctuations may extend over several decades of time ranging from a few femtoseconds to a few milliseconds. We have shown in this article how the use of a variety of theoretical methods may be utilized to detect a generic set of pathways and assess the feasibility of forming one or more transient connections. We demonstrate the application of these methods to the enzyme human carbonic anhydrase II and its mutants. Our results reveal several alternative pathways in addition to the one mediated by His-64. We also probe at length the mechanism of key conformational fluctuations contributing to the formation of the detected pathways.     

Conformational mobility of His-64 in the Thr-200----Ser mutant of human carbonic anhydrase II

The three-dimensional structure of the Thr-200----Ser (T200S) mutant of human carbonic anhydrase II (CAII) has been determined by X-ray crystallographic methods at 2.1-A resolution. This particular mutant of CAII exhibits CO2 hydrase activity that is comparable to that of the wild-type enzyme with a 2-fold stabilization of the E.HCO3- complex and esterase activity that is 4-fold greater than that of the wild-type enzyme. The structure of the mutant enzyme reveals no significant local changes accompanying the conservative T200S substitution, but an important nonlocal structural change is evident: the side chain of catalytic residue His-64 rotates away from the active site by 105 degrees about chi 1 and apparently displaces a water molecule. The displaced water molecule is present in the wild-type enzyme; however, the electron density into which this water is built is interpretable as an alternate conformation of His-64 with 10-20% occupancy. The rate constants for proton transfer from the zinc-water ligand to His-64 and from His-64 to bulk solvent are maintained in the T200S variant; therefore, if His-64 is conformationally mobile about chi 1 and/or chi 2 during catalysis, compensatory changes in solvent configuration must sustain efficient proton transfer.     

NMR study of structure and electron transfer mechanism of Pseudomonas aeruginosa azurin

The nuclear spin-spin and spin-lattice relaxation times of the C epsilon 1-proton of His-35 and the C delta 2-proton of His-46 of reduced Pseudomonas aeruginosa azurin have been determined at 298 and 320 K and at pH 4.5 and 9.0 at various concentrations of total azurin and in the presence of varying amounts of oxidized azurin. The relaxation times appear strongly influenced by the electron self-exchange reaction between oxidized and reduced protein. The T1 data of the His-35 proton have been analyzed according to the "fast-exchange limit," while the "slow-exchange limit" appears to obtain for the T2 data of the His-46 proton. Analysis of the proton relaxation data yields values of the electron self-exchange rate constants of (9.6 +/- 0.7) X 10(5) M-1 S-1 (pH 4.5) and (7.0 +/- 1.3) X 10(5) M-1 S-1 (pH 9.0) at 298 K. The dipolar correlation time amounts to 1-2.5 ns in the temperature range of 298-320 K. A Fermi-contact interaction of about 100 mG for the C delta 2-proton of His-46 is compatible with the experimental observations. The pH-induced conformational changes lead to variations on the order of about 1 A in the distance from the copper to the His-35 protons. The data implicate the "hydrophobic patch" around His-117 as the site of electron transfer in the self-exchange reaction of the azurin.     

A proton nuclear magnetic resonance study on the solution structure of crotamine

Proton nuclear magnetic resonance (NMR) spectra of crotamine, a myotoxic protein from a Brazilian rattlesnake (Crotalus durissus terrificus), have been analyzed. All the aromatic proton resonances have been assigned to amino acid types, and those from Tyr-1, Phe-12, and Phe-25 to the individual residues. ThepH dependence of the chemical shifts of the aromatic proton resonances indicates that Tyr-1 and one of the two histidines (His-5 or His-10) are in close proximity. A conformational transition takes place at acidicpH, together with immobilization of Met-28 and His-5 or His-10. Two sets of proton resonances have been observed for He-17 and His-5 or His-10, which suggests the presence of two structural states for the crotamine molecule in solution.     

Reaction of morphinone reductase with 2-cyclohexen-1-one and 1-nitrocyclohexene: proton donation, ligand binding, and the role of residues Histidine 186 and Asparagine 189

Morphinone reductase (MR) catalyzes the NADH-dependent reduction of alpha/beta unsaturated carbonyl compounds in a reaction similar to that catalyzed by Old Yellow Enzyme (OYE1). The two enzymes are related at the sequence and structural levels, but key differences in active site architecture exist which have major implications for the reaction mechanism. We report detailed kinetic and solution NMR data for wild-type MR and two mutant forms in which residues His-186 and Asn-189 have been exchanged for alanine residues. We show that both residues are involved in the binding of the reducing nicotinamide coenzyme NADH and also the binding of the oxidizing substrates 2-cyclohexen-1-one and 1-nitrocyclohexene. Reduction of 2-cyclohexen-1-one by FMNH(2) is concerted with proton transfer from an unknown proton donor in the active site. NMR spectroscopy and flavin reoxidation studies with 2-cyclohexen-1-one are consistent with His-186 being unprotonated in oxidized, reduced, and ligand-bound MR, suggesting that His-186 is not the key proton donor required for the reduction of 2-cyclohexen-1-one. Hydride transfer is decoupled from proton transfer with 1-nitrocyclohexene as oxidizing substrate, and unlike with OYE1 the intermediate nitronate species produced after hydride transfer from FMNH(2) is not converted to 1-nitrocyclohexane. The work highlights key mechanistic differences in the reactions catalyzed by MR and OYE1 and emphasizes the need for caution in inferring mechanistic similarities in structurally related proteins.     

pH dependence of tritium exchange with the C-2 protons of the histidines in bovine trypsin.

At pH 8.9 and 37 degrees C the half-times for tritium exchange with the C-2 protons of the histidines of trypsin are 73 days for His-57, and greater than 1000 days for His-40 and His-91. These half-times are much longer than the half-life of exchange for the C-2 proton of free histidine (2.8 days at pD 8.2), and longer than any previously reported half-time of exchange at pH greater than 8. These very low rates of exchange are discussed with reference to the refined structure of trypsin. The tritium exchange of His-57 depends on an apparent pKa of 6.6. This pKa may represent the pKa of the imidazole of His-57 in an inactive conformation of the enzyme.     

Potentiometric determination of ionizations at the active site of papain.

The ionization behavior of groups at the active site of papain was determined from the pH dependence of the difference of proton content of papain and the methylthio derivative of the thiol group at the active site of papain (papain-S-SCH3). This difference in proton content was determined directly by two independent methods. One method involved potentiometric measurements of the protons released and demethylthiolation of papain-S-SCH3 with dithiothreitol, as a function of pH. The other method involved analogous measurements of the protons released on methylthiolation of papain with methyl methanethiosulfonate. The methylthio pH-difference titrations generated by these measurements indicate that ionization of the thiol group at the active site of papain is linked to the ionization of His-159. The pK of the thiol group changes from 3.3 to 7.6 on deprotonation of His-159 at 29 degrees C/20.05. Similarly, the pK of His-159 shifts from 4.3 to 8.5 when the active site thiol group is deprotonated. The microscopic ionization constants determined in this work for Cys-25 and His-159 indicate that equilibrium constant for transfer of the proton from Cys-25 to His-159 is 8--12, and that in the physiological pH range the active site thiol group exists mainly as a thiol anion.     

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.     

The pH dependence of serpin-proteinase complex dissociation reveals a mechanism of complex stabilization involving inactive and active conformational states of the proteinase which are perturbable by calcium

Serpin family protein proteinase inhibitors trap proteinases at the acyl-intermediate stage of cleavage of the serpin as a proteinase substrate by undergoing a dramatic conformational change, which is thought to distort the proteinase active site and slow deacylation. To investigate the extent to which proteinase catalytic function is defective in the serpin-proteinase complex, we compared the pH dependence of dissociation of several serpin-proteinase acyl-complexes with that of normal guanidinobenzoyl-proteinase acyl-intermediate complexes. Whereas the apparent rate constant for dissociation of guanidinobenzoyl-proteinase complexes (k(diss, app)) showed a pH dependence characteristic of His-57 catalysis of complex deacylation, the pH dependence of k(diss, app) for the serpin-proteinase complexes showed no evidence for His-57 involvement in complex deacylation and was instead characteristic of a hydroxide-mediated deacylation similar to that observed for the hydrolysis of tosylarginine methyl ester. Hydroxylamine enhanced the rate of serpin-proteinase complex dissociation but with a rate constant for nucleophilic attack on the acyl bond several orders of magnitude slower than that of hydroxide, implying limited accessibility of the acyl bond in the complex. The addition of 10-100 mm Ca(2+) ions stimulated up to 80-fold the dissociation rate constant of several serpin-trypsin complexes in a saturable manner at neutral pH and altered the pH dependence to a pattern characteristic of His-57-catalyzed complex deacylation. These results support a mechanism of kinetic stabilization of serpin-proteinase complexes wherein the complex is trapped as an acyl-intermediate by a serpin conformational change-induced inactivation of the proteinase catalytic function, but suggest that the inactive proteinase conformation in the complex is in equilibrium with an active proteinase conformation that can be stabilized by the preferential binding of an allosteric ligand such as Ca(2+).     

Two conformational states of Glu242 and pKas in bovine cytochrome c oxidase.

Cytochrome c oxidase (CcO) is the terminal enzyme in the respiratory electron transport chain of aerobic organisms. It catalyses the reduction of atmospheric oxygen to water, and couples this reaction to proton pumping across the membrane; this process generates the electrochemical gradient that subsequently drives the synthesis of ATP. The molecular details of the mechanism by which electron transfer is coupled to proton pumping in CcO is poorly understood. Recent calculations from our group indicate that His291, a ligand of the Cu(B) center of the enzyme, may play the role of the pumping element. In this paper we describe calculations in which a DFT/continuum electrostatic method is used to explore the coupling of the conformational changes of Glu242 residue, the main proton donor of both chemical and pump protons, to its pKa, and the pKa of His291, a putative proton loading site of our pumping model. The computations are done for several redox states of metal centers, different protonation states of Glu242 and His291, and two well-defined conformations of the Glu242 side chain. Thus, in addition to equilibrium redox/protonation states of the catalytic cycle, we also examine the transient and intermediate states. Different dielectric models are employed to investigate the robustness of the results, and their viability in the light of the proposed proton pumping mechanism of CcO. The main results are in agreement with the experimental measurements and support the proposed pumping mechanism. Additionally, the present calculations indicate a possibility of gating through conformational changes of Glu242; namely, in the pumping step, we find that Glu242 needs to be reprotonated before His291 can eject a proton to the P-site of membrane. As a result, the reprotonation of Glu can control proton release from the proton loading site.     

Conditionally lethal mutants of bacteriophage T4 defective in production of a transfer RNA

The two buried carboxyls (Asp-102 and Asp-194) in both chymotrypsin and chymotrypsinogen are ionized at pH values greater than 4.2 and may be ionized even as low as pH 3.This was demonstrated by coupling most of the surface carboxyis of the proteins by a carbodi-imide with glycinamide or semicarbazide to diminish the groups ionizing at low pH and then titrating the proton uptake on denaturation by sodium dodecyl sulphate between pH 3.0 and 4.6. At pH values greater than 4.2 all unblocked carboxyls are ionized. The proton uptake during the conformational change on denaturation was determined by a stopped-flow procedure and found to be about 2H⁺/mol between pH 3.0 and 3.6. The rate constant for the uptake of protons is the same as that for the exposure of tryptophan and lies in the tens of millisecond region.The buried negative charge at the active site appears to be mainly on Asp-102 rather than on His-57, the pK_a of which must be raised by the buried charge. This enhances its efficacy as a base catalyst in the “charge relay system”.The presence of an intact charge relay system in the inactive zymogen illustrates the importance of stereochemical fit between enzyme and substrate. Enzyme catalysis could hardly be mediated by a catalyst which is uniquely reactive in the absence of correct enzyme-substrate orientation as this would be inconsistent with its specificity.     

Photoinduced electron transfer between the Rieske iron-sulfur protein and cytochrome c(1) in the Rhodobacter sphaeroides cytochrome bc(1) complex. Effects of pH,temperature, and driving force

Electron transfer from the Rieske iron-sulfur protein to cytochrome c(1) (cyt c(1)) in the Rhodobacter sphaeroides cytochrome bc(1) complex was studied using a ruthenium dimer complex, Ru(2)D. Laser flash photolysis of a solution containing reduced cyt bc(1), Ru(2)D, and a sacrificial electron acceptor results in oxidation of cyt c(1) within 1 micros, followed by electron transfer from the iron-sulfur center (2Fe-2S) to cyt c(1) with a rate constant of 80,000 s(-1). Experiments were carried out to evaluate whether the reaction was rate-limited by true electron transfer, proton gating, or conformational gating. The temperature dependence of the reaction yielded an enthalpy of activation of +17.6 kJ/mol, which is consistent with either rate-limiting conformational gating or electron transfer. The rate constant was nearly independent of pH over the range pH 7 to 9.5 where the redox potential of 2Fe-2S decreases significantly due to deprotonation of His-161. The rate constant was also not greatly affected by the Rieske iron-sulfur protein mutations Y156W, S154A, or S154A/Y156F, which decrease the redox potential of 2Fe-2S by 62, 109, and 159 mV, respectively. It is concluded that the electron transfer reaction from 2Fe-2S to cyt c(1) is controlled by conformational gating.     

Proton nuclear magnetic resonance characterization of phospholipase A2 from Laticauda semifasciata

The molecular properties of phospholipases (PLases) A2 I and A2 III from a sea snake, Laticauda semifasciata, have been characterized by gel-filtration, as well as proton NMR, CD, UV absorption, and fluorescence spectroscopic methods. PLase A2 I exists as a monomer in aqueous solution in the presence or in the absence of Ca2+. The dissociation constants of the Ca2+-enzyme complexes have been determined for the two enzymes. The 270-mHz proton NMR spectra of PLases A2 I and A2 III have been measured, and the aromatic proton resonances of His-21 and His-48 in the active site have been assigned. By analyzing the pH dependence of the chemical shifts of the histidine proton resonances, pKa values have been determined for His-21 and His-48 with and without Ca2+. The conformational transitions have been found to take place at low pH or at high temperature (at approximately 65 degrees C). Fluorescence change of PLase A2 I upon addition of substrate analogs suggests that Trp-70 in PLase A2 I is involved in the binding to micellar substrates. The lack of Trp-70 in PLase A2 III is probably related to the low enzymatic activity as compared with that of PLase A2 I.     

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.