Selection of Carbonic Anhydrase Variants Displayed on Phage
AROMATIC RESIDUES IN ZINC BINDING SITE ENHANCE METAL AFFINITY AND EQUILIBRATION KINETICS*

(Received for publication, May 6, 1997, and in revised form, June 4, 1997)

Jennifer A. Hunt Dagger § and Carol A. Fierke Dagger

From the Dagger  Department of Biochemistry, Duke University Medical Center, Durham, North Carolina 27710

ABSTRACT
INTRODUCTION
EXPERIMENTAL PROCEDURES
RESULTS
DISCUSSION
FOOTNOTES
REFERENCES


ABSTRACT

In all metalloenzymes, hydrophobic residues surround the metal binding site. In carbonic anhydrase II (CAII) residues Phe93, Phe95, and Trp97 flank two of the three histidines that coordinate zinc to form a hydrophobic cluster beneath the zinc binding site. A library of CAII variants differing in these hydrophobic amino acids was prepared using cassette mutagenesis, then displayed on filamentous phage, and screened for proteins retaining high zinc affinity. Wild-type CAII was enriched 20-fold by selection, and consensus residues at each position were identified from the enriched CAII variants (Ile, Phe, Leu, and Met at position 93; Ile, Leu, and Met at position 95; and Trp and Val at position 97). Highly selected variants have zinc affinity and catalytic activity nearly equal to that of wild-type CAII, indicating that the aromatic residues are not absolutely essential. However, the zinc dissociation rate constant and catalytic activity of the variants correlate with the volume of the amino acids at positions 93, 95, and 97. In summary, metalloenzyme variants displayed on phage can be selected on the basis of metal affinity; such methods will be useful for optimization of metal ion biosensors.


INTRODUCTION

Understanding the structural determinants of metal ion affinity in proteins is an important step in the design of protein metal sites. Such sites may be used to stabilize proteins, to regulate the activity of proteins, or for use in biosensors to quantify trace metal ions (see reviews in Refs. 1-3). Because of its high affinity and specificity for zinc, the His3 metal polyhedron of carbonic anhydrase (CAII)1 often has been used as a model for designing metal sites in existing proteins (4-6) and in de novo proteins such as the minibody (7) and four helical bundle protein (8, 9). While the protein metal ligands of CAII have been incorporated into these proteins with coordination geometry very similar to that observed in CAII, these designed metal sites lack the zinc avidity and catalytic activity of biological zinc sites, suggesting that further protein structural factors in CAII contribute to the metal affinity and reactivity of this enzyme.

To probe the relationship between protein structure and metal ion affinity and specificity, many of the conserved structural features near the zinc binding site of CAII have been investigated (review see Ref. 10). The refined x-ray crystal structure of CAII (11) reveals that the zinc ion is bound at the bottom of a deep active site cleft, coordinated by protein residues His94, His96, and His119. At physiological pH, a solvent hydroxide ion completes the tetrahedral coordination geometry of the bound metal. The functional and structural effects of altering the protein zinc ligands have been examined (12-14); substitution of any one of the three histidine ligands by alanine decreases the affinity for zinc by 104-fold. Additionally, the direct zinc ligands are fully saturated by hydrogen bond networks with second shell residues as follows: His94 donates a hydrogen bond to the carboxamide side chain of Gln92; His119 donates a hydrogen bond to the carboxylate side chain of Glu117; His96 donates a hydrogen bond to the backbone carbonyl oxygen of Asn244; and zinc-bound hydroxide donates a hydrogen bond to the hydroxyl side chain of Thr199 (11). This conserved network of residues that form hydrogen bonds with the His ligands enhances zinc affinity (15-17); the Q92A and E117A substitutions both decrease zinc affinity 5-10-fold and increase the zinc dissociation rate constant. Therefore, the protein structure surrounding the zinc binding site of CAII is finely tuned for high zinc affinity and slow rate constants of metal dissociation.

A further structural feature observed in all metalloproteins is that the hydrophilic direct metal ligands are embedded within a larger shell of hydrophobic groups (18). In CAII, residues Phe93, Phe95, and Trp97 flank the zinc ligands His94 and His96 in a beta -strand structure (Fig. 1), and this aromatic cluster forms packing interactions with the hydrophobic core of the protein, as visualized by x-ray crystallography (11). These bulky aromatic residues may increase metal binding and/or reactivity by sterically restricting the flexibility of the site, controlling the range of movement allowed for the residues coordinating the metal ion and affecting both KD and koff. Alternatively, the hydrophobicity of these amino acids may be important for providing an interior region with low dielectric constant to enhance electrostatic interactions (18). In this study we demonstrate that this conserved hydrophobic region is essential for increasing zinc affinity and decreasing the rate constant for zinc dissociation but has modest importance for the reactivity of the zinc-bound hydroxide. Furthermore, a correlation between the zinc dissociation rate constants and the volume of the substituted amino acid suggests that these bulky residues function mainly by decreasing the mobility of the site.


Fig. 1. Structure of the zinc binding site of wild-type CAII taken from the refined crystal structure of Hakansson et al. (11), showing the zinc tetrahedrally coordinated to His94, His96, His119 and a solvent molecule and three residues, Phe93, Phe95, and Trp97, that form an aromatic cluster underneath the zinc site. The figure was generated using MOLSCRIPT (63).
[View Larger Version of this Image (38K GIF file)]

A second goal of this work is to explore the utility of phage display technology as a means of selecting metalloenzyme variants on the basis of metal affinity. Large libraries (~109; for review, see Ref. 19) of variants can be screened using phage display; since the genetic information coding the selected protein is present in the attached phage particle, very rare variants can be identified through amplification and DNA sequencing. Phage display has proven to be a very useful tool for the in vitro development of antibodies (20). Furthermore, several enzymes, including beta -lactamase (21), trypsin (22), and alkaline phosphatase (23), have been expressed on the surface of phage and screened on the basis of their affinity for various ligands. Recently, we have demonstrated that CAII can be expressed as a fusion protein on the surface of filamentous phage in a correctly folded, active form and enriched using the affinity of the holo-enzyme for sulfonamides.2 In this study, we screen a library of CAII variants displayed on phage based on their affinity for a catalytic zinc ion and demonstrate that this method is very useful for probing the structural determinants of metal affinity and for preparing CAII variants with altered metal affinity and specificity.


EXPERIMENTAL PROCEDURES

Mutagenesis of pCANTAB-CA Phage Display Vector

In the plasmid pCANTAB-CA, the CAII gene is inserted 5' to the gene coding for the M13 minor coat protein (g3p) to create a CAII-g3p fusion protein. In this vector, the amino acid sequence Ala-Ala-Gln-Leu is inserted between the g3p signal sequence and the N terminus of CAII. At the junction between the C terminus of CAII and the N terminus of g3p, the stop codon of CAII and the lysine immediately preceding it are replaced by three Ala residues (Pharmacia Biotech Inc.).2 To facilitate cassette mutagenesis, silent mutations were added to the CAII gene producing unique KpnI and MscI sites at codons 86/87 (ggcact-ggtacc) and 102 (tggaca-tggcca), respectively, to form pCANTAB-CAII. To produce the library of CAII variants, a mutagenic cassette flanked by KpnI and MscI sites was prepared by annealing synthetic oligonucleotides composed of the coding and non-coding strands of CAII DNA from codon 87 to codon 102. This cassette contained the nucleotide sequence DBK (D = G/A/T; B = G/T/C; K = G/T) at residues 93, 95, and 97 to generate a library of 183, or 5832, DNA sequences that encode 1728 different CAII variants. Both the pCANTAB-CAII vector and the synthetic oligonucleotide were digested with the restriction enzymes KpnI and MscI, and then the cassette was ligated to the large KpnI/MscI fragment of pCANTAB-CAII using standard techniques (25, 26). To ensure that affinity enrichments were not due to contamination of the library or loss of rare variants, duplicate libraries were prepared, transformed, and screened separately.

Production of CA Phage

The library of pCANTAB-CAII plasmids was transformed into XL1-Blue cells (26), and colonies were grown for 20 h at 30 °C on SOB/agar plates (26) containing 100 µg/ml ampicillin, 12.5 µg/ml tetracycline, and 2% glucose to estimate the number of individual transformants and to prevent growth competition among the variants. In the first preparation of CAII phage library, 3,000 colonies were obtained, and 14,000 colonies were produced in a duplicate library. All of the individual transformants from both libraries were combined, suspended in 2 × yeast tryptone media (26) containing 100 µg/ml ampicillin, 12.5 µg/ml tetracycline, and 2% glucose and grown in shaker flasks at 37 °C for two doublings to reach an A600 of 0.3-0.4 before superinfection by the helper phage M13KO7 (27) at a multiplicity of infection of 100 infective colony-forming units (cfu) per cell (as estimated by A600 of 0.4 approx 108 cells/ml). After shaking at 100 rpm for 1 h, cells were transferred to glucose-free Luria Broth (LB) media containing 100 µg/ml ampicillin, 12.5 µg/ml tetracycline, and 50 µg/ml kanamycin and shaken at 300 rpm, 37 °C, for 5-6 h to produce phagemid particles. The solution was centrifuged (16,000 × g, 30 min) to remove cells. The phagemid in the supernatant was precipitated by polyethylene glycol using standard procedures (27) and resuspended in a 1/200 volume of 10 mM Tris, pH 7, 1 mM EDTA buffer (TE). The concentration of phagemid particles was estimated by infecting log-phase XL1-Blue cells with serial dilutions of the phagemid solution and plated on LB agar containing 100 µg/ml ampicillin to quantify colony-forming units. Typical phagemid yields were 2-3 × 108 cfu/ml original growth volume. DNA from individual phagemid colonies was isolated by the alkaline lysis method (26) and sequenced by the dideoxy method (28) to determine the diversity of the variant CA phage library.

Preparation of Affinity Resin and Determination of Binding Enrichments

p-Aminomethyl benzenesulfonamide was coupled to Ultralink Carboxy resin (Pierce) using carbodiimide as described previously (29).2 Affinity resin was pre-blocked by incubation with 10 mg/ml bovine serum albumin in 2 volumes of wash buffer (0.1 M Tris-SO4, pH 7.0, 0.5% Tween 20) for 1 h at 22 °C and then washed 4-5 times with 2-3 volumes of wash buffer to remove any excess bovine serum albumin. For determination of the binding enrichment, solutions of CA phagemid (109-1012 cfu) were added to 400 µl of resin in a final volume of 1 ml of wash buffer containing 1 µM ZnCl2 and incubated overnight at room temperature with rotary mixing at 7 rpm. The eluate was collected by centrifugation of the resin in a Spin-X filter tube (Costar) (1-2 min, 2,000 × g). The resin was resuspended in wash buffer and incubated for about 5 h. After 4-5 washes, CA phagemid bound to the sulfonamide resin were eluted by incubation for 16-24 h with wash buffer containing either 10 mM acetazolamide or varying concentrations of the metal ion chelator, dipicolinate (DPA). To amplify eluted phagemid for a second round of affinity selection, XL1-blue cells were infected with 103-105 cfu of phage and plated onto LB agar containing 100 µg/ml ampicillin. Fifty colonies were randomly picked; their DNA was isolated (26) and the CAII gene sequenced (28). The remainder of the transformants (about 9,000) was resuspended in liquid media and superinfected with M13KO7 phage, as described above, to generate CA phagemid.

Expression and Purification of CAII Variants

To facilitate the subcloning of CAII variants, unique KpnI and MscI sites were added to the CAII gene in the CAII expression vector pACA (30) to create the expression vector pMACA. The DNA of selected CAII variants was subcloned into pMACA by ligating the small KpnI/MscI fragment from the pCANTAB-CAII vector containing the variant sequence to the large fragment of KpnI/MscI-digested pMACA. The CAII genes of these constructed plasmids were then sequenced through both ligation sites. The plasmids encoding CAII variants were transformed into BL21(DE3) cells (31), and CAII was induced by the addition of 0.25 mM isopropyl-beta -D-thiogalactopyranoside to late log Escherichia coli BL21(DE3)pMACA cells followed by incubation at 30 °C for 5 h (32). CAII variants were purified using sequential ion exchange chromatography on DEAE-Sephacel and SP-Sepharose Fast Flow, as described previously (33).

Catalytic Activity

The kcat/KM for CAII-catalyzed p-nitrophenyl acetate (pNPA) hydrolysis was measured at 25 °C in either 50 mM MES (pH 5.5-7.0) or 50 mM Tris-SO4 (pH 7.0-9.0), 0.5 mM pNPA, 0.1- 0.3 µM CAII and the ionic strength maintained at 0.1 by the addition of Na2SO4. The initial rate of pNPA hydrolysis was monitored at 348 nm (Delta epsilon  = 5,000 M-1 cm-1) (34) and background rates of hydrolysis, measured in the presence of 50 µM acetazolamide, were subtracted to obtain the CAII-dependent rates. For some variants, the pKa was determined by fitting the data to Equation 1 using kcat/KM values measured over the pH range 5.5 to 9. For other variants, kcat/KM was measured at pH 6.5 (k6.5) and pH 8.0 (k8), and the ratio of these values was used to estimate the pKa using Equation 2.
(k<SUB><UP>cat</UP></SUB>/K<SUB>M</SUB>)<SUB><UP>obs</UP></SUB>=(k<SUB><UP>cat</UP></SUB>/K<SUB>M</SUB>)/(1+10<SUP>(<UP>p</UP>K<SUB>a</SUB><UP>−pH</UP>)</SUP>) (Eq. 1)
K<SUB>a</SUB>=[(k<SUB>6.5</SUB>/k<SUB>8</SUB>)(10<SUP><UP>−</UP>6.5</SUP>)−10<SUP><UP>−</UP>8</SUP>] /(1−k<SUB>6.5</SUB>/k<SUB>8</SUB>) (Eq. 2)

Initial rates of CO2 hydration were measured at 25 °C using a KinTek stopped-flow apparatus by the changing pH indicator (35). The reaction was monitored at 578 nm in 50 mM TAPS buffer, pH 8.5, containing 28 µM m-cresol purple, 0.1 mM EDTA, ionic strength adjusted to 0.1 with Na2SO4. The CO2 concentration was varied from 3.3 to 24 mM. For those CAII variants with rapid zinc dissociation, 0.1 mM ZnSO4 was included in the assay. Background rates measured in the absence of CAII were subtracted from the observed rates.

Zinc Dissociation Rate Constant and Dissociation Constant

All solutions were prepared in plasticware using deionized water. The rate constant for zinc dissociation was determined by diluting CAII variants (10-20 µM) into 15 mM potassium phosphate, pH 7.0, containing 35 mM EDTA at 25 °C to chelate zinc dissociating from the enzyme. Dissociation of zinc from CAII was monitored by measuring the decrease in pNPA hydrolysis activity since the apoenzyme has low catalytic activity. At various times, an aliquot of CAII was diluted 100-200-fold into assay buffer, and the esterase activity was assayed as described above. Data were fit to Equation 3 using the curve-fitting program KaleidaGraph (Synergy Software) where A0 and At are the catalytic activity at dilution and at various times after dilution, respectively. The standard errors were determined from these fits.
A<SUB>t</SUB>=A<SUB>0</SUB> e<SUP>(<UP>−</UP>kt)</SUP> (Eq. 3)

To measure the zinc dissociation constant, apo-CAII variants were prepared using Amicon diaflow filtration (36) against first 50 mM dipicolinate (DPA), pH 7.0, and then 15 mM potassium phosphate buffer, pH 7, followed by chromatography on a PD-10 column (Pharmacia) to remove excess DPA. Apo-CAII (50-80 µM) was dialyzed against varying concentrations of zinc (0-1 mM) in 0.2-3 mM DPA, 15 mM potassium phosphate buffer, pH 7.0, for 18-22 h at 30 °C (13). After equilibrium was achieved, the fraction of CAII containing a bound zinc ion ([E·Zn]) was quantified by measuring the specific pNPA hydrolysis activity, or by removing free zinc by chromatography on a PD-10 column and quantifying bound zinc using a colorimetric 4-(2-pyridylazo)resorcinol method (37). The concentration of free zinc in the dialysis buffer was calculated from the DPA-zinc stability constants at 30 °C (38). The dissociation constant and standard error were calculated using the KaleidaGraph curve-fitting program with Equation 4, varying both C and KD, where C ranged from 0.9 to 1.1.
[E · <UP>Zn</UP>]/[E]<SUB><UP>tot</UP></SUB>=C/(1+K<SUB>D</SUB>/[<UP>Zn</UP>]<SUB><UP>free</UP></SUB>) (Eq. 4)


RESULTS

Library of CAII Variants Displayed on Phage

To investigate the functional importance of a conserved hydrophobic cluster beneath the zinc binding site in CAII, we used cassette mutagenesis (25) to prepare a library of CAII variants with alterations in three residues (Phe93, Phe95, and Trp97) flanking two of the histidines (His94 and His96) that coordinate zinc (11; Fig. 1). Rather than completely randomizing the DNA sequence of each of these codons, we chose to use a semi-random mutagenesis scheme in which the codons encoding residues 93, 95, and 97 in CAII were simultaneously replaced with a mixture of codons encoding the following amino acids: Ala, Cys, Phe, Gly, Ile, Leu, Met, Arg, Ser, Thr, Val, and Trp. This scheme was chosen because it eliminates all stop codons and most charged residues (39), which are predicted to significantly destabilize the folded structure of CAII. This method creates a library of 5832 DNA sequences encoding 1728 unique CAII variants. The most rare variant is encoded by only one sequence in this coding scheme, and the most common variant is encoded by 27 sequences. These plasmids were transformed into E. coli and a library of 17,000 individual colonies were combined, as described under "Experimental Procedures." The DNA sequence of the CAII gene from 42 randomly chosen transformants was determined to assess the accuracy of the ligation reaction and the diversity of the library; 62% of the sequences contained one correctly oriented copy of the cassette. Furthermore, little nucleotide bias in residues 93, 95, and 97 was observed (first nucleotide: 42% A, 29% T, and 29% G; second nucleotide: 32% C, 32% G, and 36% T; and third nucleotide: 51% G and 49% T). One sequence encoding wild-type CAII was also observed (3.8%). This is higher than the 0.02% predicted from the randomized sequence in the starting oligonucleotide and suggests wild-type contamination from re-ligation of wild-type pCANTAB-CA when the original library was created. In the future, this problem could be alleviated by positioning a stop codon at residue 93 to prevent production of wild-type CAII. Using these data, we estimated that the library contains 10,500 (62% of 17,000) viable transformants. Based on Poisson distribution, this gives 90% confidence that all 5,832 sequences are represented in the library (40), indicating that the starting library has high diversity.

Enrichment of CAII Variants

For the purposes of selecting variants with high zinc affinity, the CAII variants were expressed as CAII-g3p fusion proteins on the external surface of phagemid, as described by Hunt and Fierke.2 These CAII phagemid were then fractionated by their ability to bind to the sulfonamide resin (described under "Experimental Procedures") where an essential interaction is the direct coordination of the sulfonamide nitrogen with the active site zinc in CAII (41). This procedure enriches CA phagemid 103-fold compared with M13 phage2 and can potentially fractionate CA variants based on their sulfonamide affinity, zinc affinity, and/or protein stability. Under the conditions of the selection, CA phagemid were added in great excess over affinity binding sites; however, the total number of binding sites (3-8 × 106)2 is much greater than the diversity of the variant library to ensure that multiple copies of desired variants are bound.

In the first round of screening, CA variants were eluted by the addition of a high affinity sulfonamide inhibitor, acetazolamide (42). This first step was designed to eliminate variants with greatly reduced zinc or sulfonamide affinity and/or severe defects in the structural integrity or stability. From this pool, 79 individual transformants were chosen at random, and the CAII gene was sequenced. In this selection step, wild-type CAII was greatly enriched (>= 10-fold) since 44% of the sequenced transformants (35 out of 79) contained the wild-type sequence compared with approx 4% wild-type sequence in the starting library. These data clearly indicate that the wild-type sequence of F93F95W97 confers significant advantages for expression of CA-g3p fusion protein and/or binding to the sulfonamide resin under the selection conditions. Additionally, the results of the first round of selection clearly indicate that positively charged residues in these positions are discriminated against in this selection. In the original variant library, arginine is the only charged amino acid possible and it comprises 5.6% of the amino acids encoded at position 93, 95, or 97. Of the variants sequenced after this first round of selection, only two contain arginine at either position 93, 95, or 97 (<1.5%), suggesting that positioning this charged residue in the hydrophobic protein core prevents correct functioning of the protein.

CAII phage from the first selection pool were amplified for a second round of screening, incubated with the sulfonamide resin, and eluted using either acetazolamide or DPA, which catalyzes the dissociation of zinc from CAII (43). Twenty-four individual transformants eluted using acetazolamide and 56 transformants eluted using dipicolinate were chosen at random, and the CAII gene of each was sequenced (Table I). The percentage of wild-type sequences was high (71%) but varied somewhat with the method of elution; 87% of the phagemid eluted with acetazolamide displayed wild-type CAII compared with 64% of the phagemid eluted with dipicolinate. These data further indicate the critical importance of the wild-type F93F95W97 sequence. In addition to wild type, several variants were identified multiple times; the sequence encoding F93M95V97 CAII was identified three times, whereas the CAII variant F93L95V97 was observed twice in the selected pool.

Table I. Amino acid sequence of CAII variants after two rounds of affinity selection


Number of variants Method of elution Amino acida
93 94 95 96 97

57 DPA, AZAb Phe His Phe His Trp
3 DPA Phe His Met His Val
2 DPA Phe His Leu His Val
1 AZA Phe His Met His Met
1 DPA Phe His Leu His Ile
1 AZA Phe His Ile His Val
1 AZA Ile His Met His Val
1 DPA Leu His Met His Val
1 DPA Met His Met His Val
1 DPA Ile His Met His Trp
1 DPA Met His Gly His Val
1 DPA Ala His Ile His Val
1 DPA Leu His Ala His Val
1 DPA Met His Leu His Trp
1 DPA Val His Ile His Trp
1 DPA Ile His Trp His Trp
1 DPA Ile His Gly His Trp
1 DPA Leu His Gly His Trp
1 DPA Met His Ala His Phe
1 DPA Leu His Gly His Phe

a Enriched amino acids are shown in boldface; wild-type amino acids are in normal type; and amino acids that are selected against are shown in italics.
b AZA, acetazolamide.

Enrichment of several amino acids compared with the starting population at each randomized position is evident (Fig. 2) and the general trend is that hydrophobic residues are enriched and hydrophilic residues are not observed. The largest enrichments (>= 3-fold) are observed for Phe, Ile, Leu, and Met at position 93, Ile, Leu, and Met at position 95, and Val and Trp at position 97. The lack of enrichment of Phe at position 95 is curious since this is the wild-type amino acid, suggesting that specific contacts between the aromatic amino acids at 93, 95, and 97 are important since we observed such a large enrichment of the wild-type sequence (20-fold). Furthermore, several amino acid residues are strongly disfavored; no Arg, Cys, Ser, or Thr residues are observed in any of the randomized positions of the selected variants, and Gly and Trp are not observed at position 93; Phe and Val are absent from position 95; and Ala, Gly, and Leu are not detected at position 97. Nonetheless, these results indicate that several similar amino acid residues can be substituted at each position without great loss of protein function, including zinc affinity and catalytic activity, as demonstrated below.


Fig. 2. Enrichment of amino acids in positions 93, 95, and 97 in CAII variants selected by affinity screening. The percent of a given amino acid at each position in the starting library (hatched lines) and in the library of variants after two rounds of sulfonamide affinity selection (solid lines) are indicated. The amino acids are listed in the order of decreasing hydrophobicity.
[View Larger Version of this Image (27K GIF file)]

Zinc Dissociation Constants of Selected CAII Variants

To evaluate the effectiveness of our selection for identifying variants with high zinc affinity and to further explore determinants of metal affinity, we used equilibrium dialysis to measure the zinc dissociation constants (KZn) for variants from the pre-selection pool and variants eluted after one and two rounds of selection (Fig. 3). After two rounds of screening, 20 CAII variants other than wild type were identified (Table I). The CAII variant observed three times in the pool (F93M95V97 CAII) has the highest zinc affinity (KZn = 1.6 pM) of all of the variants examined, comparable to the wild-type CAII zinc affinity (KZn = 0.8 pM). A similar CAII variant, F93L95V97, that is represented twice in this pool has a slightly larger KZn value of 2.2 pM. The zinc affinity of other CAII variants selected through two rounds of screening is uniformly high, with KZn values ranging from 1.6 to 12 pM (Table II). The CAII variant with the weakest measured zinc affinity, M93M95V97, contains substitutions at each of the three randomized positions. These data clearly demonstrate that aromatic residues at these positions are not absolutely required for high zinc avidity in CAII.


Fig. 3. Measurement of zinc dissociation constants for representative CAII variants from sequential rounds of screening. CAII apoenzyme (0.5 ml of approx 60 µM) was dialyzed for 20 h at 30 °C against 0.5 liters of a zinc sulfate (0-1 mM)/dipicolinate (0.2-3 mM) metal ion buffer in 15 mM potassium phosphate buffer, pH 7 (wild-type (bullet ), F93M95V97 (triangle ), F93L95V97 (black-diamond ), and A93L95W97 (open circle ) CAII). Enzyme-bound zinc was calculated from the specific catalytic activity for pNPA hydrolysis. The zinc dissociation constant was calculated from a fit of these data to Equation 4 using the KaleidaGraph curve fitting program and the values of KZn are listed in Table II.
[View Larger Version of this Image (21K GIF file)]

Table II. Zinc dissociation and dissociation rate constants of CAII variants


Selection round Variant (amino acids 93-97) KZna koffb t1/2c kond

pM h-1 × 103 days µM-1 s-1
2 (wild-type) FHFHW 0.8  ± 0.1 0.30  ± 0.06 95 0.1
2 (3 ×) FHMHV 1.6  ± 0.2 2.1  ± 0.3 14 0.35
2 FHIHV 1.6  ± 0.2 2.7  ± 0.3 11 0.47
2 (2 ×) FHLHV 2.0  ± 0.2 11  ± 2 2.5 1.5
1 FHTHW 3  ± 2 3  ± 0.2 7.2 0.35
1 FHIHS 6  ± 2 144  ± 3 0.2 6.3
2 IHMHW 8  ± 2 0.4  ± 0.1 74 0.014
2 MHLHW 10  ± 2 0.8  ± 0.1 30 0.022
2 IHMHV 11  ± 1 5.8  ± 0.4 5 0.13
2 MHMHV 12  ± 1 8.3  ± 0.6 3.6 0.18
1 MHGHW 22  ± 9 1.0  ± 0.1 39 0.018
0 SHLHM 29  ± 14 96  ± 8 0.3 0.92
1 AHLHW 30  ± 3 5.5  ± 0.3 5.0 0.05
1 MHSHW 35  ± 8 5.8  ± 1.5 5.5 0.075
0 THSHV 92  ± 41 120  ± 30 0.24 0.36
0 SHTHM 76  ± 17 280  ± 14 0.1 1.0

a Measured at pH 7.0, 30 °C, as described in the legend of Fig. 3.
b Measured at pH 7.0, 25 °C as described in the legend of Fig. 4.
c Calculated from t1/2 = ln 2/koff.
d Calculated from kon = koff/KZn.

The zinc dissociation constants of CAII variants selected after only one round of screening are weaker, ranging from 3.2 to 35 pM with an average zinc affinity of 22 pM (Fig. 3 and Table II). To evaluate the range of zinc affinities of CAII variants in the original library, three variants from the pre-selection pool containing residues that appear to be disfavored in selection were also characterized (Fig. 3 and Table II). For these CAII variants, both the zinc affinity is weaker (KZn = 30-90 pM) and the expression levels in E. coli are decreased to approx 10% that of wild-type CAII, perhaps due to lower protein stability. Both of these factors likely contribute to the negative selection of these variants using the phage display technology. Nonetheless, all of the data indicate that sulfonamide resin selection of CAII phage identifies variants with high metal affinity (Table II); however, the additional requirements that the protein be stably expressed on the surface of phage and bind sulfonamides tightly means that not all of the variants with high zinc affinity will be observed. Therefore, this technique selects proteins with both high affinity for the ligand and reasonable expression levels and protein stability.

Zinc Dissociation Rate Constants of Selected CAII Variants

In wild-type CAII, zinc not only binds tightly but the dissociation rate constant is also very slow (Table II). To probe the role of residues 93, 95, and 97 in modulating the slow zinc dissociation and association rate constants, the zinc dissociation rate constant (koff) (Fig. 5 and Table II) of a number of variants was measured. The zinc dissociation rate constants were determined from the time-dependent disappearance of CAII-bound zinc after dilution into high concentrations of EDTA to trap the dissociated zinc. Unlike dipicolinic acid which efficiently catalyzes the dissociation of zinc from wild-type CAII (43), EDTA is believed to act only by lowering the free metal ion concentration (44). Substitutions in the aromatic residues flanking the histidine ligands have much greater effects upon the zinc dissociation rate constant than upon KZn. The koff values for zinc in these variants range from the CAII variant I93M95W97 with a koff comparable to that of wild-type CAII to S93T95M97 CAII with a koff that is 1000-fold faster than that of wild-type CAII (Fig. 4 and Table II).


Fig. 5. Correlation of zinc dissociation and dissociation rate constants for CAII variants with the volume of residues in the aromatic cluster. A, the value of log KZn (Table II) for variants containing W97 (black-triangle) or other amino acids (V, S, M) at position 97 (triangle ) are plotted as a function of the additive volume (47) of the amino acids at positions 93, 95 and 97. The data are fit to a line yielding slope = -0.009 ± 0.002, r = 0.86 (- - -) and slope = -0.011 ± 0.002, r = 0.94 (---), respectively. B, the value of log koff for CAII variants in Table II (open circle ) excluding M93G95W97 CAII (bullet ) are plotted as a function of the additive volume (47) of the amino acids at positions 93, 95 and 97. The data are fit to a line yielding slope = -0.012 ± 0.001, r = 0.92.
[View Larger Version of this Image (22K GIF file)]


Fig. 4. Measurement of zinc dissociation rate constants for CAII variants. Each variant (10-20 µM final concentration; wild-type (square ), M93L95W97 (black-triangle), M93M95V97 (open circle ), S93L95M97 (black-diamond ), F93I95S97 (triangle ) and S93T95M97 (bullet ) CAII) was diluted into 15 mM potassium phosphate buffer, pH 7.0, containing 35 mM EDTA, at 25 °C. Dissociation of zinc was monitored by measuring the pNPA hydrolysis activity of aliquots withdrawn at various times. The zinc dissociation rate constant was calculated from a fit of these data to Equation 3 using the KaleidaGraph curve fitting program and the values of koff are listed in Table II.
[View Larger Version of this Image (25K GIF file)]

From these data several generalizations can be drawn. The four variants with the fastest dissociation rate constants contain a Ser residue at one of the three randomized positions. Replacing the large aromatic residues at positions 93, 95, or 97 with smaller residues may increase the conformational mobility of the beta -strand, allowing the adjacent histidines that coordinate zinc (His94 and His96) additional flexibility. Furthermore, the variants with the slowest dissociation rate constant retain Trp at position 97, and, regardless of the nature of amino acid substitutions at positions 93 and 95, the half-time for zinc dissociation of the isolated variants that retain Trp97 is greater than 5 days. In a comparison of the CAII variants I93M95W97 and I93M95V97, the values of KZn are comparable, whereas the dissociation rate constant increases 10-fold for substitution of Trp97 with Val. These data suggest that Trp97 may have an important role in anchoring the beta -strand to reduce the conformational flexibility of His94 and His96 required for zinc dissociation.

The association rate constants for zinc can be estimated from the measured KZn and koff values assuming a single association step where kon = koff/KZn (Table II). Several variants in which the wild-type residue Trp97 is retained have zinc association rate constants that are 5-10-fold slower than that of wild type CAII; substituting Ser at position 93 or 97 leads to a 10- or 60-fold increase, respectively, in the zinc association rate constant. These data further indicate that residues 93, 95, and 97 may restrict the conformational flexibility of the histidine ligands by anchoring the beta -strand on which they are located.

Catalysis and pKa of Zinc-bound Solvent

To delineate the role of the aromatic cluster beneath the zinc binding site on the reactivity and pKa of the zinc-bound solvent, we measured the pH dependence of the second-order rate constant for pNPA hydrolysis catalyzed by the CAII variants. As observed for wild-type CAII, the pH dependence of pNPA hydrolysis of all of the variants is consistent with the ionization of a single enzymic group. The pH dependence of the esterase activity of wild-type CAII, and likely these variants as well, directly reflects ionization of the zinc-water moiety to zinc-hydroxide, which acts as a nucleophile during catalysis (for review, see Ref. 45). The pKa values determined for all of the selected CAII variants examined (Table III) are the same, within experimental error, as that of wild-type CAII (6.8 ± 0.1; 46). Since the pKa of the zinc-water is indicative of the ability of the zinc to stabilize a negative charge (10), these data suggest that the substituted amino acids at positions 93, 95, and 97 do not significantly affect the electrostatic environment of the zinc binding site. However, the pH-independent kcat/KM for pNPA hydrolysis catalyzed by these CAII variants decreases 2-30-fold compared with that of wild-type CAII (Table III). Furthermore, log kcat/KM for pNPA ester hydrolysis increases roughly proportionally to the additive volume (47) or overall hydrophobicity of the substituted side chains at positions 93, 95, and 97 as measured by the transfer free energy of amino acids between octanol and water (48, 49). Similar dependence of the esterase activity on hydrophobicity has previously been observed for variations in either the substrate or amino acids Val121, Val143, and Leu198 in the active site of CAII (32, 46, 50, 51). This dependence on hydrophobicity may be due mainly to increased affinity of the ester substrate, rather than enhanced reactivity of the nucleophilic zinc-bound hydroxide.

Table III. Catalytic activity of selected carbonic anhydrase variants

The following four variants are not listed in the table because all of the catalytic parameters were within a factor of 2.4 of the wild-type CAII values: FHIHV, FHTHW, MHLHW, and AHLHW.

CAII variant, amino acids 93-97 Esterase activity,a kcat/KM CO2 hydrase activityb
kcat/KM KM kcat

M-1 s-1 µM-1 s-1 mM µs-1
FHFHW 2400  ± 20 89  ± 6 9  ± 1 0.77  ± 0.03
FHMHV 420  ± 20 73  ± 4 10  ± 1 0.69  ± 0.02
FHLHV 520  ± 10 89  ± 4 7.9  ± 0.5 0.70  ± 0.02
FHIHS 1300  ± 20 6  ± 1 64  ± 47 0.38  ± 0.22
IHMHW 700  ± 36 75  ± 2 7.6  ± 0.3 0.6  ± 0.1
IHMHV 500  ± 20 80  ± 10 8  ± 2 0.59  ± 0.05
MHMHV 280  ± 10 81  ± 3 7.8  ± 0.5 0.63  ± 0.01
MHGHW 360  ± 10 34  ± 3 9  ± 1 0.30  ± 0.02
SHLHM 640  ± 30 10  ± 1 35  ± 14 0.3  ± 0.1
MHSHW 760  ± 30 57  ± 1 10  ± 4 0.6  ± 0.1
THSHV 80  ± 6 1.2  ± 0.1 30  ± 9.1 0.04  ± 0.007
SHTHM 160  ± 14

a Measured at pH 8.0, 25 °C as described under "Experimental Procedures."
b Measured at pH 8.5, 25 °C as described under "Experimental Procedures."

There is considerable evidence that hydration of CO2 by wild-type CAII consists of two main steps (for review, see Ref. 45) as shown in Equations 5 and 6: 1) nucleophilic attack of zinc-bound hydroxide on CO2 to form enzyme-bound HCO3- followed by product dissociation resulting in the zinc-water form of the enzyme; 2) proton transfer to solvent via a proton shuttle (His64) to regenerate the zinc-hydroxide species. For wild-type CAII, kcat/KM is a combination of the association and CO2 hydration steps (kcat/KM approx  k1 as k2 > k-1). At high CO2 and buffer concentrations, intramolecular proton transfer between zinc-water and His64 becomes rate-limiting as indicated by a solvent isotope effect of 3-4 on kcat (kcat approx  k3 as k4 > k-3) (52). Therefore, KM reflects this change in rate-limiting steps, not CO2 affinity which is estimated at 0.1 M (51, 53).
E-<UP>OH<SUP>−</SUP></UP>+<UP>CO</UP><SUB>2</SUB>+<UP>H<SUB>2</SUB>O</UP> <LIM><OP><ARROW>↔</ARROW></OP><LL>k<SUB><UP>−</UP>1</SUB></LL><UL>k<SUB>1</SUB></UL></LIM> E-<UP>HCO</UP><SUP><UP>−</UP></SUP><SUB>3</SUB>+<UP>H<SUB>2</SUB>O</UP> <LIM><OP><ARROW>↔</ARROW></OP><LL>k<SUB><UP>−</UP>2</SUB></LL><UL>k<SUB>2</SUB></UL></LIM> E-<UP>OH</UP><SUB>2</SUB>+<UP>HCO</UP><SUP><UP>−</UP></SUP><SUB>3</SUB> (Eq. 5)
E-<UP>OH</UP><SUB>2</SUB>+B <LIM><OP><ARROW>↔</ARROW></OP><LL>k<SUB><UP>−</UP>3</SUB></LL><UL>k<SUB>3</SUB></UL></LIM><SUP><UP>+</UP></SUP><UP>H</UP>-E-<UP>OH</UP><SUP><UP>−</UP></SUP>+B <LIM><OP><ARROW>↔</ARROW></OP><LL>k<SUB><UP>−</UP>4</SUB></LL><UL>k<SUB>4</SUB></UL></LIM> E-<UP>OH<SUP>−</SUP></UP>+B<UP>H<SUP>+</SUP></UP> (Eq. 6)
The steady-state kinetic parameters for CO2 hydration, catalyzed by the CAII variants with substitutions in the aromatic cluster, were measured at pH 8.5, 25 °C (Table III). For 11 of the 14 variants studied, alterations in Phe93, Phe95, and Trp97 decrease kcat/KM for CO2 hydration less than 3-fold. These data clearly indicate that aromatic amino acids underneath the zinc binding pocket are not essential for high catalytic activity in CAII and that the majority of substitutions in these amino acids do not significantly alter the reactivity or position of the nucleophilic zinc-bound hydroxide. However, the selection procedure may enrich variants with high catalytic activity since the variants must retain affinity for sulfonamides which are proposed as transition state analogs (54).

Similarly, small decreases in kcat were measured for these CAII variants, indicating that the rate constant for proton transfer from zinc-water to His64 is slightly perturbed by these substitutions. However, the effect on kcat does not parallel changes in the pKa of the zinc-solvent ligand, as observed for transfer of a proton between an acid to a base in model reactions and between His64 and bases (55-57). Therefore, the decreases in kcat are likely caused by small changes in the positioning of the active site solvent/His64 proton transfer pathway, as has been observed for amino acid substitutions at positions near His64 (58, 59).

Interestingly, the variants with the largest increase in the zinc dissociation rate constant (>300-fold relative to wild-type CAII) also display the biggest decreases in CO2 hydration. For CAII variants F93I95S97, S93L95M97, and T93S95V97, kcat/KM is decreased 9-74-fold relative to wild-type CAII. If kcat/KM reflects the association and hydration of CO2 for these variants as it does in the case of wild-type CAII, then substrate association and/or reactivity is hindered in these variants, possibly due to increased conformational flexibility in the active site or re-positioning of the zinc-hydroxide. The decreases in kcat/KM for these variants are accompanied by large increases in KM so kcat is only 2-20-fold lower than that of wild-type CAII. These results indicate that CO2 association and hydration are impaired in these variants while effects on proton transfer are less severe.


DISCUSSION

Effectiveness of Phage Display

Recently, phage display technology has emerged as a powerful tool for the design and re-design of proteins (19, 20, 60). This method of rapidly screening large protein libraries for rare variants has been applied to diverse uses such as the selection of linkers for catalytic antibodies (61) and the selection of enzymes based on substrate affinity (21, 22). In this study we describe the successful adaptation of phage display methods to include the selection of metalloenzyme variants based on metal affinity.

Our results show that carbonic anhydrase variants displayed on phage can be screened on the basis of metal affinity. Using this method we have investigated the functional role of three residues, Phe93, Phe95, and Trp97, flanking amino acids His94 and His96 that coordinate zinc. This hydrophobic cluster beneath the metal binding site is highly conserved among CA isozymes (62) and is proposed to modulate the metal affinity and specificity of this protein. Conserved hydrophobic clusters observed near metal binding sites in a number of other proteins may play similar roles (18).

To investigate the role of these residues a pool of variants in which Phe93, Phe95, and Trp97 have been semi-randomly substituted was prepared and screened for the ability to bind zinc tightly. A striking result is that after two rounds of selection and elution with sulfonamides the majority of variants (87%) contained the wild-type amino acid sequence. These data indicate that the wild-type sequence is optimized to achieve a combination of high expression levels, tight zinc affinity, and catalytic activity (as indicated by high sulfonamide affinity (54)). Additional diversity in the selected pool was achieved by using a metal chelator to elute the CA phage from the column. This is understandable since this elution method should select for properties not optimized in wild-type CAII, including decreased sulfonamide affinity and increased zinc dissociation rate constants.

Additionally, non-wild-type sequences were also selected at each substituted position, pinpointing a set of consensus residues required for tight zinc binding; variants that contain the residues most highly enriched during selection were found to have zinc dissociation constants nearly identical to that of the wild-type protein. Variants isolated after only one round of selection have, on average, a lower affinity for zinc, whereas three variants that contain residues that were not enriched by selection had the lowest zinc affinity, with KZn values 70-fold higher than that of wild-type CAII. Overall, phage display proved to be a very powerful tool in selecting variants with zinc affinity in the top 2% of the range of zinc affinities through only two rounds of selection.

Function of Hydrophobic Cluster: Zinc Affinity

The enormous selection of the wild-type F93F95W97 CAII sequence observed by phage display indicates that this specific sequence confers a selective advantage over any other sequence included in this library of variants. This advantage is not due solely to high zinc affinity since variants with substitutions in residues 93, 95, and 97 are able to retain high zinc affinity; in particular, the free energy for binding zinc is decreased by <=  0.5 kcal/mol at 25 °C for variants containing Phe at position 93; Met, Ile, or Leu at position 95; and Val at position 97. The overall selection indicates a higher amino acid diversity (Fig. 2), but this includes variants with a free energy of binding that is decreased up to 1.6 kcal/mol. Of the three amino acid positions varied in this study, metal affinity is influenced most by alterations at position 93. Variants that retain the wild-type Phe93 bind zinc more tightly than any of the variants with substitutions at position 93. Additionally, substitutions at positions 95 and 97 influence zinc affinity to a lesser extent; the substitutions Phe95 right-arrow Thr (F93T95W97 versus wild-type CAII) and Trp97 right-arrow Ser (F93I95S97 versus F93I95V97) decrease zinc affinity only 4-fold.

The selection data (Fig. 2) clearly indicate that hydrophobic amino acids are preferred at all three positions. The hydrophobicity of amino acids surrounding metal sites is proposed to decrease the dielectric constant of the region to enhance electrostatic interactions (18). To investigate the characteristics of amino acids at positions 93, 95, and 97 important for high metal affinity in CAII, we have probed the dependence of the metal affinity on the combined hydrophobicity (as indicated by the free energy of transfer between octanol and water (48, 49)) and volume (47) of the amino acids at these positions. For both these parameters, a low correlation coefficient (R <=  0.74) is observed when the KZn values for all of the variants are included. However, if the proteins containing a Trp at position 97 are excluded, the log KZn shows a high correlation (R >=  0.93) with both hydrophobicity and volume (Fig. 5A). These data clearly indicate that larger and more hydrophobic amino acids in this region increase the metal affinity. On the other hand, the metal affinity of the variants containing Trp97 show a significantly decreased correlation with either hydrophobicity (R = 0.83) or volume (R = 0.86) and, furthermore, the metal affinity of these variants is lower than predicted (Fig. 5A). These data suggest that changing the size of the amino acid from Val to Trp at position 97 does not significantly enhance zinc affinity. The hydrophobicity or size of the amino acid at positions 93 and 95 may increase metal affinity mainly by decreasing the mobility of the beta -sheet to "preorganize" the protein side chains that coordinate metal and decrease the entropy loss accompanying metal binding. Similarly, calorimetric studies indicate that the hydrogen bonds between protein side chains and the histidine ligands in CAII decrease the entropy of zinc binding.3

Zinc Equilibration

In wild-type CAII, zinc not only binds tightly (1 pM), but the dissociation rate constant is very slow, and the observed rate constant for zinc association with CAII (kon; estimated from the measured KZn and koff)) is also slower than the diffusion-controlled limit (Table II). The slow zinc association rate constant suggests a two-step mechanism for the binding of zinc to wild-type CAII where the rate of formation of the metal complex is limited by an intramolecular step, such as the dissociation of inner-sphere water molecules (64-66). One possible mechanism for the association of zinc with wild-type CAII that is consistent with the properties of first coordination sphere and second shell variants (13, 15) is the rapid formation of an initial complex in which zinc is coordinated to two protein ligands, e.g. His94 and His96, followed by the slow formation of the tetrahedral complex comprised of His94, His96, His119, and a solvent molecule. The rate constant for this interconversion increases as the mobility of the site increases; for example, deletion of a hydrogen bond between His119 and Glu117 greatly increases the rate constant of zinc dissociation (15-17).

Similarly, the data in this study indicate that substitution of any one of the three aromatic residues with a small residue, such as Ser or Ala, significantly increases the rate constant of zinc dissociation (Table II). Furthermore, for all of the variants without a Trp at position 97, the association rate constant increases 1.3-10-fold more than either the dissociation rate constant or the dissociation constant, indicative of an increase in the rate-limiting intramolecular step. Finally, a high correlation between the log koff and either the additive volume (R = 0.92 excluding M93G95W97 CAII; Fig. 5B) or hydrophobicity (Refs. 48 and 49; R = 0.89) of the substituted amino acid indicates that the rate constant for zinc dissociation increases as the size or hydrophobicity of the amino acids underneath the zinc binding site decreases. In this case, 15 of the 16 CAII variants fit on a single line indicating that the Trp at position 97 is likely conserved mainly to decrease koff rather than to maintain high zinc affinity. The value of koff for one variant, M93G95W97 CAII, is slower than predicted for unknown reasons. Taken together, these data suggest that residues 93, 95, and 97 play a steric role in reducing the conformational flexibility of the beta -strand on which they are located and that flexibility of this strand containing His94 and His96 is a major factor controlling the rate constant of zinc dissociation.

Catalytic Activity

Amino acid substitutions in residues 93, 95, and 97 likely do not influence the electrostatic environment of the active site of CAII, as the pKa of the zinc-water is unaltered. However, the value of kcat/KM for the esterase activity decreases as the volume of the amino acids at these positions decreases. On the other hand, kcat/KM for CO2 hydration catalyzed by CAII variants is unaffected (<= 1.6-fold decrease) by reducing the additive volume of these three amino acids from 608 Å3 (wild-type CAII) to 466 Å3. As the volume of the amino acids at positions 93, 95, and 97 shrinks further, the value of kcat/KM decreases (by as much as 74-fold for T93S95V97 CAII, volume = 346 Å3). This break could reflect a change in the rate-limiting step of kcat/KM from CO2 association to CO2 hydration, indicating that the reactivity of zinc-hydroxide does depend on the volume or hydrophobicity of these residues. Alternatively, these data could indicate that the position of zinc and zinc-hydroxide is unaffected by substitutions in the aromatic cluster until the volume of these amino acids are significantly decreased; at this point structural changes occur that alter the ground state position of the active site nucleophile and decrease the observed reactivity.

Comparison of Phage Display Selection with Natural Selection

As deduced amino acid sequences have been determined for a number of CA proteins from different species (67), it is of interest to compare the results of the directed evolution methods of this study with data from evolutionary selection. At position 93, the wild-type residue Phe was clearly the most highly enriched residue, a result that correlates well with the finding that the six variants with highest zinc affinity retained Phe93 (Table II). Leu, Ile, and Met were also enriched at this position. Interestingly, two of these three residues (Leu and Met) occur in this position in several different CA proteins, including human CAIV, and only residues Phe, Leu, and Met are present in this position in catalytically active CA (67). In this case, phage display has selected variants that closely mimic the structures allowed through evolution. Slightly different results were observed at position 95, where residues Leu, Ile, and Met were again enriched but Phe was not observed in the non-wild-type sequences. It is not known why the wild-type residue Phe was not enriched during selection. Again, Ile and Leu occur at this position in sequences of naturally occurring CA proteins (67), although Met has not yet been observed.

Phage display enrichment at position 97 provided a very interesting result: Val was greatly enriched, the wild-type Trp was less enriched, and no other residues were selected at this position. Interestingly, in all known CA sequences Trp is found at position 97 (67). Surprisingly, of the variants characterized in this study, those containing Val at position 97 bind zinc with affinity near that of wild type, demonstrating that this residue can substitute for Trp with little loss of zinc affinity. Even variants containing Ser at position 97 bind zinc with KZn values less than 7 pM. Furthermore, the CO2 hydrase activity of variants containing Val97 is also near the wild-type value. However, as discussed above, altering position 97 greatly increases (7-1000-fold) the rate constant for zinc dissociation compared with that of wild-type CAII, which may explain why Trp97 is invariant in all known CA proteins with catalytic zinc ions. These data suggest that while Trp97 is not essential for high zinc affinity, this particular residue may have been preserved through evolution due to its role in maintaining a slow rate constant for metal dissociation from the protein. Alternatively, the stability of CAII may be significantly influenced by the amino acid at position 97 and this may be reflected in an increased zinc dissociation rate constant. These results demonstrate that phage display can be used to dissect the specific functional roles of conserved residues as well as to prepare variants with specific functions not selected by evolution.

Optimization of Metal Ion Biosensors

The ability to predict the properties of a given protein structure is necessary for rational protein design and redesign. Correlations between the zinc binding properties and residue size or hydrophobicity in CAII are of particular interest in the development of CAII-based biosensors for the detection of zinc or other specific metal ions. A CAII sensor that accurately measures levels of zinc in complex media has been designed (24, 68). However, to expand the zinc concentration range and speed of measurement of this biosensor, CAII variants with altered zinc affinities and zinc equilibration rates would be highly useful. The high correlations between the size or hydrophobicity of residues at positions 93, 95, and 97 in CAII and the zinc affinity and dissociation rate constants suggest the possibility of designing CAII variants with known zinc binding properties that can be used in a rational fashion to optimize the properties of this biosensor.


FOOTNOTES

*   This work was supported in part by the National Institutes of Health Grant GM40602 and the Office of Naval Research.The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
§   Supported in part by National Institutes of Health Postdoctoral Fellowship GM17467.
   Recipient of an American Heart Association Established Investigator Award and a David and Lucile Packard Foundation Fellowship in Science and Engineering. To whom correspondence should be addressed. Tel.: 919-684-2557; Fax: 919-684-8885.
1   The abbreviations used are: CAII, human carbonic anhydrase II; g3p, gene 3 protein of M13 phage; cfu, colony-forming units; DPA, dipicolinate; MES, 2-(N-morpholino)ethanesulfonic acid; pNPA, p-nitrophenol acetate; TAPS, 3-[tris(hydroxymethyl)methyl]aminopropanesulfonic acid.
2   J. A. Hunt and C. A. Fierke, submitted for publication.
3   M. Mahapatro, E. J. Toone, and C. A. Fierke, submitted for publication.

REFERENCES

  1. Hellinga, H. W. (1996) Curr. Opin. Biotechnol. 7, 437-441 [CrossRef][Medline] [Order article via Infotrieve]
  2. Tainer, J. A., Roberts, V. A., and Getzoff, E. D. (1992) Curr. Opin. Biotechnol. 3, 378-387 [Medline] [Order article via Infotrieve]
  3. Higaki, J. N., Fletterick, R. J., and Craik, C. S. (1992) Trends Biochem. Sci. 17, 100-104 [CrossRef][Medline] [Order article via Infotrieve]
  4. Roberts, V. A., Iverson, B. L., Iverson, S. A., Benkovic, S. J., Lerner, R. A., Getzoff, E. D., and Tainer, J. A. (1990) Proc. Natl. Acad Sci. U. S. A. 87, 6654-6658 [Abstract]
  5. Muller, H. N., and Skerra, A. (1994) Biochemistry 33, 14126-14135 [Medline] [Order article via Infotrieve]
  6. Vita, C., Roumestand, C., Toma, F., and Menez, A. (1995) Proc. Natl. Acad. Sci. U. S. A. 92, 6404-6408 [Abstract]
  7. Pessi, A., Bianchi, E., Crameri, A., Venturini, S., Tramontano, A., and Sollazo, M. (1993) Nature 362, 367-369 [CrossRef][Medline] [Order article via Infotrieve]
  8. Handel, T., and DeGrado, W. F. (1990) J. Am. Chem. Soc. 112, 6710-6711
  9. Regan, L., and Clarke, N. D. (1990) Biochemistry 29, 10878-10883 [Medline] [Order article via Infotrieve]
  10. Christianson, D. W., and Fierke, C. A. (1996) Acc. Chem. Res. 29, 331-339 [CrossRef]
  11. Hakansson, K., Carlsson, M., Svensson, L. A., and Liljas, A. (1992) J. Mol. Biol. 227, 1192-1204 [Medline] [Order article via Infotrieve]
  12. Kiefer, L. L., Ippolito, J. A., Fierke, C. A., and Christianson, D. W. (1993) J. Am. Chem. Soc. 115, 12581-12582
  13. Kiefer, L. L., and Fierke, C. A. (1994) Biochemistry 33, 15233-15240 [Medline] [Order article via Infotrieve]
  14. Ippolito, J. A., and Christianson, D. W. (1994) Biochemistry 33, 15241-15249 [Medline] [Order article via Infotrieve]
  15. Kiefer, L. L., Paterno, S. A., and Fierke, C. A. (1995) J. Am. Chem. Soc. 117, 6831-6837
  16. Lesburg, C. A., and Christianson, D. W. (1995) J. Am. Chem. Soc. 117, 6838-6844
  17. Huang, C. C., Lesburg, C. A., Kiefer, L. L., Fierke, C. A., and Christianson, D. W. (1996) Biochemistry 35, 3439-3446 [CrossRef][Medline] [Order article via Infotrieve]
  18. Yamashita, M. M., Wesson, L., Eisenman, G., and Eisenberg, D. (1990) Proc. Natl. Acad. Sci. U. S. A. 87, 5648-5652 [Abstract/Free Full Text]
  19. Marks, J. D., Hoogenboom, H. R., Griffiths, A. D., and Winter, G. (1992) J. Biol. Chem. 267, 16007-16010 [Free Full Text]
  20. Hawkins, R. E., Russell, S. J., and Winter, G. (1992) J. Mol. Biol. 226, 889-896 [Medline] [Order article via Infotrieve]
  21. Soumillion, P., Jespers, L., Bouchet, M., Marchand-Brynaert, J., Winter, G., and Fastrez, J. (1994) J. Mol. Biol. 237, 415-422 [CrossRef][Medline] [Order article via Infotrieve]
  22. Corey, D. R., Shiau, A. K., Yang, Q., Janowski, B. A., and Craik, C. S. (1993) Gene (Amst.) 128, 129-134 [CrossRef][Medline] [Order article via Infotrieve]
  23. McCafferty, J., Jackson, R. H., and Chiswell, D. (1991) Protein Eng. 4, 955-961 [Abstract]
  24. Thompson, R. B., Ge, Z., Patchan, M. W., and Fierke, C. A. (1996) J. Biomed. Optics 1, 131-137
  25. Wells, J. A., Vasser, M., and Powers, D. B. (1985) Gene (Amst.) 34, 315-323 [CrossRef][Medline] [Order article via Infotrieve]
  26. Sambrook, J., Fritsch, E. F., and Maniatis, T. (1989) Molecular Cloning: A Laboratory Manual, Vol. 2, Cold Spring Harbor Press, Cold Spring Harbor, NY
  27. Viera, J., and Messing, J. (1987) Methods Enzymol. 153, 3-11 [Medline] [Order article via Infotrieve]
  28. Sanger, F., Nicklen, S., and Coulson, A. R. (1977) Proc. Natl. Acad. Sci. U. S. A. 74, 5463-5467 [Abstract]
  29. Osborne, W. R. A., and Tashian, R. E. (1975) Anal. Biochem. 64, 297-303 [Medline] [Order article via Infotrieve]
  30. Fierke, C. A., Krebs, J. F., and Venters, R. A. (1991) in International Workshop on Carbonic Anhydrase (Botre, F., Gros, G., and Storey, B. T., eds), pp. 22-36, VCH, Weinheim, Germany
  31. Studier, F. W., and Moffat, B. A. (1986) J. Mol. Biol. 189, 113-130 [Medline] [Order article via Infotrieve]
  32. Nair, S. K., Calderone, T. L., Christianson, D. W., and Fierke, C. A. (1991) J. Biol. Chem. 266, 17320-17325 [Abstract/Free Full Text]
  33. Krebs, J. F., and Fierke, C. A. (1993) J. Biol. Chem. 268, 948-954 [Abstract/Free Full Text]
  34. Armstrong, J. McD., Myers, D. V., Verpoorte, J. A., and Edsall, J. T. (1966) J. Biol. Chem. 241, 5137-5149 [Abstract/Free Full Text]
  35. Khalifah, R. G. (1971) J. Biol. Chem. 246, 2561-2573 [Abstract/Free Full Text]
  36. Alexander, R. S., Kiefer, L. L., Fierke, C. A., and Christianson, D. W. (1993) Biochemistry 32, 1510-1518 [Medline] [Order article via Infotrieve]
  37. Hunt, J. B., Neece, S. H., Schachman, H. K., and Ginsburg, A. (1984) J. Biol. Chem. 259, 14793-14803 [Abstract/Free Full Text]
  38. Sillen, L. G., and Martell, A. E. (1964) Stability Constants of Metal Ion Complexes, Special Publication No. 17, p. 546, Chemical Society, London
  39. Arkin, A. P., and Youvan, D. C. (1992) Bio/Technology 10, 297-300 [Medline] [Order article via Infotrieve]
  40. Hutchison, C. A., III, Nordeen, S. K., Vogt, K., and Edgell, M. H. (1986) Proc. Natl. Acad. Sci. U. S. A. 83, 710-714 [Abstract]
  41. Eriksson, A. E., Kylsten, P. M., Jones, T. A., and Liljas, A. (1988) Proteins Struct. Funct. Genet. 4, 283-293 [Medline] [Order article via Infotrieve]
  42. King, R. W., and Burgen, A. S. V. (1976) Proc. R. Soc. Lond. B 193, 107-125
  43. Hunt, J. B., Rhee, M., and Storm, C. B. (1977) Anal. Biochem. 64, 297-303
  44. Lindskog, S. (1982) Adv. Inorg. Biochem. 4, 115-170
  45. Silverman, D. N., and Lindskog, S. (1988) Acc. Chem. Res. 21, 30-36
  46. Fierke, C. A., Calderone, T. L., and Krebs, J. F. (1991) Biochemistry 30, 11054-11063 [Medline] [Order article via Infotrieve]
  47. Zamyatnin, A. A. (1972) Prog. Biophys. Mol. Biol. 24, 107-123 [CrossRef][Medline] [Order article via Infotrieve]
  48. Fauchere, J. L., and Pliska, V. (1983) Eur. J. Med. Chem. 18, 369-375
  49. Eisenberg, D., and McLachlan, A. D. (1986) Nature 319, 199-203 [Medline] [Order article via Infotrieve]
  50. Pocker, Y., and Beug, M. W. (1972) Biochemistry 11, 698-707 [Medline] [Order article via Infotrieve]
  51. Krebs, J. F., Rana, F., Dluhy, R. A., and Fierke, C. A. (1993) Biochemistry 32, 4496-4505 [Medline] [Order article via Infotrieve]
  52. Steiner, H., Jonsson, B. H., and Lindskog, S. (1975) Eur. J. Biochem. 59, 253-259 [Abstract]
  53. Riepe, M. E., and Wang, J. H. (1968) J. Biol. Chem. 243, 2779-2787 [Abstract/Free Full Text]
  54. Kumar, K., King, R. W., and Carey, P. R. (1976) Biochemistry 15, 2195-2202 [Medline] [Order article via Infotrieve]
  55. Rowlett, R. S., and Silverman, D. N. (1982) J. Am. Chem. Soc. 104, 6737-6741
  56. Ren, X., Tu, C., Laipis, P. J., and Silverman, D. N. (1995) Biochemistry 34, 8492-8498 [Medline] [Order article via Infotrieve]
  57. Silverman, D. N., Tu, C., Chen, X., Tanhauser, S. M., Kresge, A. J., and Laipis, P. J. (1993) Biochemistry 32, 10757-10762 [Medline] [Order article via Infotrieve]
  58. Jackman, J. E., Merz, K. M., Jr., and Fierke, C. A. (1996) Biochemistry 51, 16421-16428 [CrossRef]
  59. Tamai, S., Waheed, A., Cody, L. B., and Sly, W. S. (1996) Proc. Natl. Acad. Sci. U. S. A. 93, 13647-13652 [Abstract/Free Full Text]
  60. Rebar, E. J., and Pabo, C. O. (1994) Science 263, 671-673 [Medline] [Order article via Infotrieve]
  61. Tang, Y., Jiang, N., Parakh, C., and Hilvert, D. (1996) J. Biol. Chem. 271, 15682-15686 [Abstract/Free Full Text]
  62. Dodgson, S. J. (1991) in The Carbonic Anhydrases: Cellular Physiology and Molecular Genetics (Dodgson, S. J., Tashian, R. E., Gros, G., and Carter, N. D., eds), Plenum Publishing Corp., New York
  63. Kraulis, P. (1991) J. Appl. Crystallogr. 24, 946-950 [CrossRef]
  64. Eigen, M., and Hammes, G. G. (1963) Adv. Enzymol. 25, 1-38 [Medline] [Order article via Infotrieve]
  65. Holyer, R. H., Hubbard, C. D., Kettle, S. F. A., and Wilkins, R. G. (1965) J. Inorg. Chem. 4, 929-935
  66. Holyer, R. H., Hubbard, C. D., Kettle, S. F. A., and Wilkins, R. G. (1966) J. Inorg. Chem. 5, 622-625
  67. Hewett-Emmett, D., and Tashian, R. E. (1996) Mol. Phylogenet. Evol. 5, 50-77 [CrossRef][Medline] [Order article via Infotrieve]
  68. Thompson, R. B., and Patchan, M. W. (1995) Anal. Biochem. 227, 123-128 [CrossRef][Medline] [Order article via Infotrieve]

©1997 by The American Society for Biochemistry and Molecular Biology, Inc.