Catalytic enhancement of human carbonic anhydrase III by replacement of phenylalanine-198 with leucine

Glioblastomas with copy number gains in EGFR and RNF139 show increased expressions of carbonic anhydrase genes transformed by ENO1

Background

Prominence of glycolysis in glioblastomas may be non-specific or a feature of oncogene-related subgroups (i.e. amplified EGFR, etc.). Relationships between amplified oncogenes and expressions of metabolic genes associated with glycolysis, directly or indirectly via pH, were therefore investigated.

Methods

Using multiplex ligation-dependent probe amplification, copy numbers (CN) of 78 oncogenes were quantified in 24 glioblastomas. Related expressions of metabolic genes encoding lactate dehydrogenases (LDHA, LDHC), carbonic anhydrases (CA3, CA12), monocarboxylate transporters (SLC16A3 or MCT4, SLC16A4 or MCT5), ATP citrate lyase (ACLY), glycogen synthase1 (GYS1), hypoxia inducible factor-1A (HIF1A), and enolase1 (ENO1) were determined in 22 by RT-qPCR. To obtain supra-glycolytic levels and adjust for heterogeneity, concurrent ENO1 expression was used to mathematically transform the expression levels of metabolic genes already normalized with delta-delta crossing threshold methodology.

Results

Positive correlations with EGFR occurred for all metabolic genes. Significant differences (Wilcoxon Rank Sum) for oncogene CN gains in tumors of at least 2.00-fold versus less than 2.00-fold occurred for EGFR with CA3's expression (p < 0.03) and for RNF139 with CA12 (p < 0.004). Increased CN of XIAP associated negatively. Tumors with less than 2.00-fold CN gains differed from those with gains for XIAP with CA12 (p < 0.05). Male gender associated with CA12 (p < 0.05).

Conclusions

Glioblastomas with CN increases in EGFR had elevated CA3 expression. Similarly, tumors with RNF149 CN gains had elevated CA12 expression.

General significance

In larger studies, subgroups of glioblastomas may emerge according to oncogene-related effects on glycolysis, such as control of pH via effects on carbonic anhydrases, with prognostic and treatment implications.

•

PCR of glioblastomas show oncogene copy numbers relate to metabolic gene expressions.

•

ENO1(ENOLASE1) transformations yielded “supra-glycolytic” metabolic gene expressions.

•

EGFR, RNF139, and XIAP associated with expressions of two carbonic anhydrase genes.

•

Male gender associated (+) with the transformed expression of carbonic anhydrase 12.

•

Oncogene amplifications may aid control of pH to protect glycolysis in glioblastomas.

Enzyme-accelerated and structure-guided crystallization of calcium carbonate: role of the carbonic anhydrase in the homologous system

The calcareous spicules from sponges, e.g. from Sycon raphanus, are composed of almost pure calcium carbonate. In order to elucidate the formation of those structural skeletal elements, the function of the enzyme carbonic anhydrase (CA), isolated from this species, during the in vitro calcium carbonate-based spicule formation, was investigated. It is shown that the recombinant sponge CA substantially accelerates calcium carbonate formation in the in vitro diffusion assay. A stoichiometric calculation revealed that the turnover rate of the sponge CA during the calcification process amounts to 25 CO2s(-1) × molecule CA(-1). During this enzymatically driven process, initially pat-like particles are formed that are subsequently transformed to rhomboid/rhombohedroid crystals with a dimension of ~50 μm. The CA-catalyzed particles are smaller than those which are formed in the absence of the enzyme. The Martens hardness of the particles formed is ~4 GPa, a value which had been determined for other biogenic calcites. This conclusion is corroborated by energy-dispersive X-ray spectroscopy, which revealed that the particles synthesized are composed predominantly of the elements calcium, oxygen and carbon. Surprising was the finding, obtained by light and scanning electron microscopy, that the newly formed calcitic crystals associate with the calcareous spicules from S. raphanus in a highly ordered manner; the calcitic crystals almost perfectly arrange in an array orientation along the two opposing planes of the spicules, leaving the other two plane arrays uncovered. It is concluded that the CA is a key enzyme controlling the calcium carbonate biomineralization process, which directs the newly formed particles to existing calcareous spicular structures. It is expected that with the given tools new bioinspired materials can be fabricated.

Discovering isozyme-selective inhibitor scaffolds of human carbonic anhydrases using structural alignment and de novo drug design approaches

The development of isozyme-selective carbonic anhydrase inhibitors is currently still a great challenge. In the present study, protein-ligand complex structures were obtained by AutoDock Vina with SBR ((R)-N-(3-indol-1-yl-2-methyl-propyl)-4-sulfamoyl-benzamide) as the only inhibitor docked into the binding pockets of human isozymes CA I, II, IV, VI, IX, XII, and XIII. To make the spatial structures of complexes comparable, the co-ordinates for CA domains were reassigned based on structural alignments. With preferred docking poses of SBR been reduced to seed structures, the LigBuilder was used to build up inhibitor molecules. The results suggested that sulfonamide derivatives with naphthalene, fluorene, and acridan as the scaffold structures can be the potential isozyme-selective CAIs, especially for isozymes CA II, IV, and IX.

Transport activity of the sodium bicarbonate cotransporter NBCe1 is enhanced by different isoforms of carbonic anhydrase

Transport metabolons have been discussed between carbonic anhydrase II (CAII) and several membrane transporters. We have now studied different CA isoforms, expressed in Xenopus oocytes alone and together with the electrogenic sodium bicarbonate cotransporter 1 (NBCe1), to determine their catalytic activity and their ability to enhance NBCe1 transport activity. pH measurements in intact oocytes indicated similar activity of CAI, CAII and CAIII, while in vitro CAIII had no measurable activity and CAI only 30% of the activity of CAII. All three CA isoforms increased transport activity of NBCe1, as measured by the transport current and the rate of intracellular sodium rise in oocytes. Two CAII mutants, altered in their intramolecular proton pathway, CAII-H64A and CAII-Y7F, showed significant catalytic activity and also enhanced NBCe1 transport activity. The effect of CAI, CAII, and CAII mutants on NBCe1 activity could be reversed by blocking CA activity with ethoxyzolamide (EZA, 10 µM), while the effect of the less EZA-sensitive CAIII was not reversed. Our results indicate that different CA isoforms and mutants, even if they show little enzymatic activity in vitro, may display significant catalytic activity in intact cells, and that the ability of CA to enhance NBCe1 transport appears to depend primarily on its catalytic activity.

Atypical Protonation States in the Active Site of HIV-1 Protease: A Computational Study

The HIV protease (HIVP) is a prominent example for successful structure-based drug design. Besides its pharmaceutical impact, it is a well-studied system for which, as experimentally evidenced, protonation changes in the active site occur upon ligand binding. Therefore, it serves as an ideal candidate for a case study using our newly developed partial charge model, which was optimized toward the application of Poisson-Boltzmann based pK(a) calculations. The charge model suggests reliably experimentally determined protonation states in the active site of HIVP. Furthermore, we perform pKa calculations for two HIVP complexes with novel types of inhibitors developed and synthesized in our group. For these complexes, no experimental knowledge about the protonation states is given. For one of the compounds, containing a central pyrrolidine ring, the calculations predict that both catalytic aspartates should be deprotonated upon ligand binding.

Preorganization of the Hydroxyethylene Dipeptide Isostere: The Preferred Conformation in Solution Resembles the Conformation Bound to BACE

Conformational analysis in solution of beta-secretase inhibitors 1 and 2 by NMR spectroscopy reveals that the hydroxyethylene isostere, an apparently flexible fragment widely used as a scissile bond replacement in aspartic protease inhibitors, exists in one predominant conformation in solution. This preferred conformation is similar to that adopted by the hydroxyethylene core of 1 in complex with beta-secretase and that adopted by hydroxyethylene cores of related compounds when bound to aspartic proteases, indicating that this structural unit is preorganized in solution.

A combined QM/MM approach to protein--ligand interactions: polarization effects of the HIV-1 protease on selected high affinity inhibitors

HIV-1 protease inhibitors are one of the two widely used therapeutic agents for the treatment of HIV-infected patients. The investigation of HIV-1 protease-inhibitor interactions can provide further insight for developing new compounds that are still required due to the growing problem of drug resistance. To this end, a combined QM/MM approach was used to determine electrostatic and polarization interactions on three high affinity inhibitors, nelfinavir, mozenavir, and tipranavir. The present computational results show that explicit treatment of the polarization effect is particularly important since it can contribute as much as one-third of the total electrostatic interaction energy. Further, an amino acid decomposition analysis was applied to determine contributions of individual residues to the enzyme--inhibitor interactions. It was found that the 4-hydroxy-dihydropyrone substructure of tipranavir is especially suited for extended charge delocalization by interacting with the catalytic aspartates and isoleucines of the HIV-1 protease. The calculated electron density difference maps reaffirm and provide a means of visualizing these results.

Modeling the protonation states of the catalytic aspartates in beta-secretase

Beta-secretase (BACE) is a critical enzyme in the production of beta-amyloid, a protein that has been implicated as a potential cause of Alzheimer's disease (AD). There are two aspartic acid residues (Asp 32 and Asp 228) present in the catalytic region of BACE that can adopt multiple protonation states. The protonation state and precise location of the protons for these two residues, particularly in the presence of an inhibitor, are subjects of great interest since they have a direct bearing on the mechanism of aspartyl proteases and efforts to model beta-secretase. We have carried out full liner-scaling quantum mechanical (QM) calculations that include Poisson-Boltzmann solvation in order to identify the preferred protonation state and proton location in the presence and absence of an inhibitor. These calculations favor the monoprotonated state in the presence of ligand, and di-deprotonated state in the absence of ligand. Further the proton in the monoprotonated state is located on the inner oxygen of Asp 228. These results have implications for the catalytic mechanism of BACE and related aspartyl proteases. They also provide a reference state for the protein in structure-based modeling studies of this therapeutically important target.

Calculation of the binding affinity of beta-secretase inhibitors using the linear interaction energy method

It has been shown that the rate-limiting step in the production of beta-amyloid peptide (Abeta) is the proteolytric cleavage of the membrane-bound beta-amyloid precursor protein (APP) by beta-secretase (BACE). Since the accumulation of Abeta has been implicated as one of the key events in the progression of Alzheimer's disease, BACE has become an important therapeutic target. Recently, two crystal structures of BACE cocrystallized with the inhibitors OM99-2 and OM00-3 were published by Tang and co-workers. In addition, the Ghosh group has published binding data on a series of inhibitors based on their initial lead, OM99-2. Using this set as a basis, we have developed a model for the binding affinity of these ligands to BACE using the linear interaction energy method. The best binding affinity model for the full set of ligands had a RMSD of 1.10 kcal/mol. The best model excluding the two charged ligands had a RMSD of 0.87 kcal/mol.

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.

Carbonic anhydrase III in liver and muscle of male rats purification and properties

Cytosolic carbonic anhydrases CAI, CAII, and CAIII from liver, and CAII, and CAIII from muscle of adult male Sprague-Dawley rats were purified to homogeneity. CAIII from liver and muscle had the same amino acid composition and were immunochemically similar. Their kinetic properties at 0 degrees C were also similar. Km(CO2) was 4 mM and kcat 3x105 s(-1). Ki was 0.4 and 0.2 M for acetazolamide and NaCl, respectively. Both CAIIIs ran as single bands on SDS-electrophoresis and high-speed centrifugation, with a mol wt of 29.3 kDa. Their hydrodynamic properties suggest that CAIII is a compact, nearly spherical molecule. It contained 0.9 M zinc per M protein. In both tissues isoelectric focusing identified neutral and acidic isoforms with pIs near 7.0 and 6.3, respectively. These forms were immunologically identical and had the same amino acid composition and mol wts. The acidic forms probably represent subspecies of CAIII in different states of oxidation. CAIII is the major soluble protein in rat liver and muscle. Its function is probably to protect proteins of these tissues from oxidation catalyzed by iron-containing degradation products of haemoglobin and myoglobin. Liver CAI and CAII and muscle CAII were identical to CAI and CAII of rat erythrocytes.

Follow the protons: a low-barrier hydrogen bond unifies the mechanisms of the aspartic proteases

Seven proton transfers in five steps participate in a catalytic turnover of an aspartic protease. The Rosetta Stone for elucidating their role is a low-barrier hydrogen bond that holds the two aspartic carboxyls in a coplanar conformation. The proton of this bond shuttles between oxygens during chemical steps via hydrogen tunneling, unlike in previous proposals where it was transferred to substrate. After the release of products, both carboxyls are protonated and the bond is missing. Re-forming the bond is a significant step within a kinetic isomechanism. The bond also explains-at long last-the extremely low pK in pH profiles.

Ab initio molecular dynamics-based assignment of the protonation state of pepstatin A/HIV-1 protease cleavage site

A recent 13C NMR experiment (Smith et al. Nature Struct. Biol. 1996, 3, 946-950) on the Asp 25-Asp25' dyad in pepstatin A/HIV-1 protease measured two separate resonance lines, which were interpreted as being a singly protonated dyad. We address this issue by performing ab initio molecular dynamics calculations on models for this site accompanied by calculations of 13C NMR chemical shifts and isotopic shifts. We find that already on the picosecond time-scale the model proposed by Smith et al. is not stable and evolves toward a different monoprotonated form whose NMR pattern differs from the experimental one. We suggest, instead, a different protonation state in which both aspartic groups are protonated. Despite the symmetric protonation state, the calculated 13C NMR properties are in good agreement with the experiment. We rationalize this result using a simple valence bond model, which explains the chemical inequality of the two C sites. The model calculations, together with our calculations on the complex, allow also the rationalization of 13C NMR properties on other HIV-1 PR/inhibitor complexes. Both putative binding of the substrate to the free enzyme, which has the dyad singly protonated (Piana, S.; Carloni, P. Proteins: Struct., Funct., Genet. 2000, 39, 26-36), and pepstatin A binding to the diprotonated form are consistent with the inverse solvent isotope effect on the onset of inhibition of pepsin by pepstatin and the kinetic iso-mechanism proposed for aspartic proteases (Cho, T.-K.; Rebholz, K.; Northrop, D.B. Biochemistry 1994, 33, 9637-9642).

Structure based prediction of binding affinity of human immunodeficiency virus-1 protease inhibitors

A series of computations were performed to derive a strategy for the prediction of binding affinities of non-peptidic human immunodeficiency virus-1 (HIV-1) protease inhibitors. This paper describes the development of a 3D quantitative structure-activity relationship (3D-QSAR) methodology by using receptor information of HIV-1 protease. The docking and molecular dynamics simulations were performed on a model ligand/enzyme complex to optimize the variables involved in the generation of ligand/enzyme models. The protonation scheme of the active site aspartic acid residues of HIV-1 protease was derived from a computational study. The active site aspartate is monoprotonated with a proton placed on the OD1 atom of the ASP B25. This protocol of docking and molecular dynamics (MD) simulation was then used to derive the ligand-enzyme complexes of the molecules used in the present study. The molecular mechanics interaction descriptors were calculated from these ligand/enzyme models. A partial least squares (PLS) method was used to derive a linear correlation between the interaction descriptors and the biological activity. A good correlation was observed when the change in the energy of the ligand upon complex formation and the electrostatic contributions to the solvation energy of the ligand were included in the QSAR analysis. A highest cross-validated q2 value of 0.649 was observed. This model had a conventional r2 of 0.725, and when this model was used to predict the activity of the external test set, it produced a r2pred of 0.761. The total interaction energy was partitioned into interactions in different subsites and interactions with each of the amino acid residues of the enzyme. The PLS analysis using these descriptors helped to identify the important interactions which can be exploited for the design of HIV-1 protease inhibitors.

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.

Structure and mechanism of carbonic anhydrase

Carbonic anhydrase (CA; carbonate hydro-lyase, EC 4.2.1.1) is a zinc-containing enzyme that catalyzes the reversible hydration of carbon dioxide: CO2+ H2O<-->HCO3(-)+H+. The enzyme is the target for drugs, such as acetazolamide, methazolamide, and dichlorphenamide, for the treatment of glaucoma. There are three evolutionarily unrelated CA families, designated alpha, beta, and gamma. All known CAs from the animal kingdom are of the alpha type. There are seven mammalian CA isozymes with different tissue distributions and intracellular locations, CA I-VII. Crystal structures of human CA I and II, bovine CA III, and murine CA V have been determined. All of them have the same tertiary fold, with a central 10-stranded beta-sheet as the dominating secondary structure element. The zinc ion is located in a cone-shaped cavity and coordinated to three histidyl residues and a solvent molecule. Inhibitors bind at or near the metal center guided by a hydrogen-bonded system comprising Glu-106 and Thr-199. The catalytic mechanism of CA II has been studied in particular detail. It involves an attack of zinc-bound OH- on a CO2 molecule loosely bound in a hydrophobic pocket. The resulting zinc-coordinated HCO3- ion is displaced from the metal ion by H2O. The rate-limiting step is an intramolecular proton transfer from the zinc-bound water molecule to His-64, which serves as a proton shuttle between the metal center and buffer molecules in the reaction medium.

Constructing protein models for ligand-receptor binding thermodynamic simulations: an application to a set of peptidometic renin inhibitors

Structure-based design is the application of ligand-receptor modeling to predict the activity of a series of molecules that bind to a common receptor for which the molecular geometry is available. Successful structure-based design requires an accurate receptor model which can be economically employed in the design calculations. One goal of the work reported here has been to reduce the size of a model structure of a macromolecular receptor to allow multiple ligand-receptor molecular dynamic (MD) simulations to be computationally economical yet still provide meaningful binding thermodynamic data. A scaled-down 10 A receptor model of the enzyme renin, when subjected to an alternate atomic mass constraint, maintains the structural integrity of the composite parent crystal structure. A second goal of the work has been to develop schemes to explore and characterize the protonation states of receptors and ligand-receptor systems. Application of the charge state characterization schemes to the hydroxyethylene and statine transition state inhibitors of renin in the training set suggests a monoprotonation state of the two active-site aspartate residues, where the lone proton resides on the outer carboxylate oxygen of Asp226 is most likely. For the reduced amide transition state inhibitors an active site consisting of both aspartates in the totally ionized state, and the ligand carrying a net +1.0 charge, is most stable and consistent with experimental data.

Prediction of ligand-receptor binding thermodynamics by free energy force field (FEFF) 3D-QSAR analysis: application to a set of peptidometic renin inhibitors

A methodology is presented and applied in which the accurate estimation of ligand-receptor binding thermodynamics is achieved by formulating the calculation as a QSAR problem. When the receptor geometry is known, the free energy force field (FEFF) ligand-receptor binding energy terms can be calculated and used as independent variables in constructing FEFF 3D-QSARs. The FEFF 3D-QSAR analysis of a series of transition state inhibitors of renin was carried out. From a statistical analysis of the free energy contributions to the binding process, FEFF 3D-QSARs were constructed that reveal the change in solvation free energy upon binding and the intramolecular vacuum internal energy of the ligand in the unbound state are the most significant FEFF terms in determining the binding free energy, delta G. Other terms, such as ligand stretching, bending, and torsion energy changes, the intermolecular van der Waals interaction energy, and change in ligand conformational entropy upon binding, are also found to make significant contributions in some FEFF 3D-QSAR delta G models and in delta H and delta S binding models. Overall, a relatively small number of the thermodynamic contributions to the ligand-receptor binding process dominates the thermodynamics of binding in a given model.

Discovery of novel, non-peptide HIV-1 protease inhibitors by pharmacophore searching

Fifteen novel non-peptide HIV-1 protease inhibitors were identified by flexible 3D database pharmacophore searching of the NCI DIS 3D database. The pharmacophore query used in the search was derived directly from the X-ray determined structures of protease/inhibitor complexes. These 15 inhibitors, belonging to nine different chemical classes, are promising leads for further development. The two best inhibitors found, NSC 32180, a "dimer" of 4-hydroxycoumarin, and NSC 117027, a "tetramer" of 2-hydroxy quinone, had ID50 values of 0.32 and 0.75 microM for HIV-1 protease inhibition, respectively, and two other inhibitors had ID50 values close to 1 microM. Among the potent inhibitors, NSC 158393 not only demonstrated activity against HIV-1 protease (ID50 1.7 microM) but also exhibited promising antiviral activity in HIV-1-infected CEM-SS cells (EC50 = 11.5 microM). Validation of the pharmacophore used in the search was accomplished by conformational analysis. The binding modes of the most potent inhibitor found in our studies, NSC 32180, were predicted employing docking and molecular dynamics techniques.

Structural basis of inhibitor affinity to variants of human carbonic anhydrase II

The activities and structures of certain L198 variants of human carbonic anhydrase II (CAII) have been reported recently [Krebs, J. F., Rana, F., Dluhy, R. A., & Fierke, C. A. (1993) Biochemistry 32, 4496-4505; Nair, S. K., & Christianson, D. W. (1993) Biochemistry 32, 4506-4514]. In order to understand the structural basis of enzyme-inhibitor affinity, we now report the dissociation rate and equilibrium constants for acetazolamide and dansylamide binding to 13 variants of CAII containing substituted amino acids at position 198. These data indicate that inhibitor affinity is modulated by the hydrophobicity and charge of the 198 side chain. Furthermore, we have determined crystal structures of L198R, L198E, and L198F CAIIs complexed with the transition state analog acetazolamide. The substituted benzyl side chain of L198F CAII does not occlude the substrate association pocket, and it is therefore not surprising that this substitution has minimal effects on catalytic properties and inhibitor binding. Nevertheless, the F198 side chain undergoes a significant conformation change in order to accommodate the binding of acetazolamide; the same behavior is observed for the engineered side chain of L198R CAII. In contrast, the engineered side chain of L198E CAII does not alter its conformation upon inhibitor binding. We conclude that the mobility and hydrophobicity or residue 198 side chains affect enzyme-inhibitor (and enzyme-substrate) affinity, and these structure-function relationships are important for understanding the behavior of carbonic anhydrase isozyme III, which bears a wild-type F198 side chain.(ABSTRACT TRUNCATED AT 250 WORDS)

Kinetic and spectroscopic studies of hydrophilic amino acid substitutions in the hydrophobic pocket of human carbonic anhydrase II

The functional importance and structural determinants of a conserved hydrophobic pocket in human carbonic anhydrase II (CA II) were probed by preparing and characterizing 13 amino acid substitutions at Leu-198, situated at the mouth of the pocket. The pH dependence of the esterase activity reveals that activity decreases (up to 120-fold) as the amino acid size and charge at position 198 are varied while the pKa of the zinc-bound water molecule increases (up to 1 pH unit). Intriguingly, the pH dependence of the Leu-198-->Glu substitution is parabolic (pKas approximately 6 and 9), consistent with introduction of a general base-catalyzed mechanism. Kinetic characterization of CO2/HCO3- interconversion catalyzed by four variants (Leu-198-->Ala, His, Arg, and Glu) reveals that increasing the size of the hydrophobic pocket (Ala) does not compromise catalysis (approximately 3-fold decrease); however, substitution of charged (Arg and Glu) and larger (His) amino acids decreases kcat/KM for CO2 hydration substantially (17-fold, 19-fold, and 10-fold, respectively) but not completely. log kcat/KM for CO2 hydration, HCO3- dehydration, and p-nitrophenyl acetate hydrolysis correlates with the hydrophobicity of the residue at 198, likely reflecting desolvation or electrostatic destabilization of the ground state. The X-ray crystal structures of the Leu-198-->His, Glu, and Arg variants (Nair & Christianson, 1993) indicate that the His and Glu side chains are accommodated by minor structural reorganization leading to a wider mouth for the hydrophobic pocket while the Arg side chain blocks the pocket. Infrared spectroscopy of CO2 bound to either wild-type CA II or the Leu-198-->Arg variant indicates that the Arg substitution both decreases the affinity and alters the position of CO2 binding, suggesting that the hydrophobic pocket forms the CO2 binding site in CA II. Finally, a 1.5-fold increase (Leu-198-->Ala) and 12-fold decrease (Leu-198-->Arg) in kcat for CO2 hydration, indicative of the rate constant for intramolecular proton transfer from zinc-bound water to His-64, are likely mediated by changes in the active site solvent structure.

Refined structure of bovine carbonic anhydrase III at 2.0 A resolution

The three-dimensional structure of bovine carbonic anhydrase III (BCA III) from red skeletal muscle cells has been determined by molecular replacement methods. The structure has been refined at 2.0 A resolution by both constrained and restrained structure-factor least squares refinement. The current crystallographic R-value is 19.2% and 121 solvent molecules have so far been found associated with the protein. The structure is highly similar to the refined structure of human carbonic anhydrase II. Some differences in amino acid sequence and structure between the two isoenzymes are discussed. In BCA III, Lys 64 and Arg 91 (His 64 and Ile 91 in HCA II) are both pointing out from the active site cavity forming salt bridges with Glu 4 and Asp 72 (His 4 and Asp 72 in HCA II), respectively. However, Arg 67 and Phe 198 (Asn 67 and Leu 198 in HCA II) are oriented towards the zinc ion and significantly reduce the volume of the active site cavity. Phe 198 particularly reduces the size of the substrate binding region at the "deep water" position at the bottom of the cavity and we suggest that this is one of the major reasons for the differences in catalytic properties of isoenzyme III as compared to isozyme II.

Refined structure of the aminobenzolamide complex of human carbonic anhydrase II at 1.9 A and sulphonamide modelling of bovine carbonic anhydrase III

The binding of aminobenzolamide to human carbonic anhydrase (HCA II) has been investigated by X-ray crystallography. The three dimensional atomic structure of the enzyme inhibitor complex has been refined at 1.9 A resolution. The crystallographic R-factor is 17.8%. All inhibitor atoms are clearly possible to identify from the difference electron density map in the active site of the enzyme. The nitrogen of the sulphonamide group of the inhibitor is bound as a fourth ligand to the zinc ion, the other three are all histidyl residues. The binding conformation of the sulphonamide groups is similar to the previously described sulphonamide inhibitors. One of the oxygens of the outer sulphonamido group of the inhibitor forms a hydrogen bond to the amino group of Gln 92. The higher affinity of the benzolamide inhibitor compared with acetazolamide can be accounted for by the strong aromatic and hydrophobic interactions between the amino benzene ring of the inhibitor and the residues Phe 131 and Leu 198. In modelling studies of bovine carbonic anhydrase III (BCA III) it was evident that Phe 198 prevents an optimal interaction with sulphonamides.

Catalysis by mutants of human carbonic anhydrase II: effects of replacing hydrophobic residues 198 and 204

Previous studies shows that the replacement of Phe-198 in carbonic anhydrase III to the corresponding Leu residue found in carbonic anhydrase II caused the appearance of isozyme II-like activity (LoGrasso et al. (1991) Biochemistry 30, 8463-8470). Carbonic anhydrase II is more efficient in the catalysis of CO2 hydration by 500-fold and has an apparent pKa for this catalysis about two pKa units above that of carbonic anhydrase III. Moreover, isozyme II catalyzes the hydrolysis of 4-nitrophenyl acetate, whereas isozyme III shows no appreciable catalysis. The purpose of this work was to test the hypothesis that making the converse replacement Leu-198-->Phe as well as Leu-204-->Glu and the double replacement in carbonic anhydrase II would give the resulting mutants of isozyme II properties of isozyme III. The catalytic activities of these mutants in CO2 hydration and 4-nitrophenyl acetate hydrolysis were smaller by at most 5-fold and the pKa values for these catalyses were identical compared with wild-type isozyme II. The different effects of converse mutants of HCA II and III indicate complexity in structure not evident from their similar backbone conformations.

A T7 expression vector optimized for site-directed mutagenesis using oligodeoxyribonucleotide cassettes

Site-directed mutagenesis is widely used to examine structure/function relationships in proteins. We have designed a bacterial expression vector series which is optimized for efficient site-directed mutagenesis and subsequent protein synthesis without intervening subcloning steps. The vectors, derived from the T7 expression vectors of Studier and his collaborators [Studier et al., Methods Enzymol. 185 (1990) 60-89], are small and have a bacteriophage f1 origin of replication for production of single-stranded (ss) DNA. Both single-site mutants [using ssDNA and mutating oligodeoxyribonucleotides (oligos)] and cassette mutants (mutagenesis of a short region by inserting double-stranded oligos into unique restriction sites) are rapidly synthesized and expressed with these vectors. Vector construction and use are detailed with examples showing the expression of the sequences encoding human carbonic anhydrases II and III. Production levels of greater than 60 mg of protein per liter of culture have been obtained.