(Received for publication, July 17, 1995; and in revised form, September 6, 1995)
From the
We constructed a bacteriophage-displayed library containing
randomized mutations at H chain residues 30-35 of the
anti-digoxin antibody 26-10 Fab to investigate sequence
constraints necessary for high affinity binding in an antibody of known
crystal structure. Phage were selected by panning against digoxin and
three C-16-substituted analogues. All antigen-positive mutants selected
using other analogues also bound digoxin. Among 73 antigen-positive
clones, 26 different nucleotide sequences were found. The majority of
Fabs had high affinity for digoxin (K
3.4
10
M
) despite wide sequence diversity. Two
mutants displayed affinities 2- and 4-fold higher than the parental
antibody.
Analysis of the statistical distribution of sequences showed that highest affinity binding occurred with a restricted set of amino acid substitutions at positions H33-35. All clones save two retained the parental Asn-H35, which contacts hapten and hydrogen bonds to other binding site residues in the parental structure. Positions H30-32 display remarkable diversity, with 10-14 different substitutions for each residue, consistent with high affinity binding. Thus complementarity can be retained and even improved despite diversity in the conformation of the N-terminal portion of the H-CDR1 loop.
Although our understanding of the mechanisms for antibody (Ab) ()diversification has increased rapidly from studies at the
DNA level using monoclonal antibodies (mAbs), we have insufficient
tools to accurately relate the primary structure of hypervariable
regions to the three-dimensional structure dictating binding
specificity. Wu and Kabat (1) proposed that the hypervariable
regions located within variable (V) regions of heavy (H) and light (L)
chains were the sites of amino acid residues conferring Ab specificity
and folded together to form Ab combining sites. Subsequently, x-ray
crystallographic analyses on hapten binding and anti-protein Fab
fragments (reviewed in (2, 3, 4, 5) ) revealed that the
antigen combining cavities are indeed composed of the H and L chain
hypervariable or complementarity-determining regions (CDRs). Certain
CDR sequences have been shown, in general, to function in the context
of different framework sequences, as the framework three-dimensional
folding is remarkably uniform(6, 7) . Moreover, there
is evidence for a small repertoire of homologous
``canonical'' CDR loop conformations (8, 9) . The correlation between antibody-combining
site structure and binding function has been studied intensively using
mutagenesis, modeling, and three-dimensional structure determination.
The model system we use for studies of antigen-Ab complementarity
through protein engineering is based on a set of murine mAbs to the
cardiac glycoside digoxin. Several features of these antibodies and the
haptens lend themselves to this objective. Cardiac glycosides consist
of a cardenolide steroid ring structure with a C-14 -OH, attached
-glycosidic moieties at C-3 and a 17
unsaturated
five-membered lactone ring (Fig. 1). Digoxin is a relatively
rigid large hydrophobic hapten without charged groups. There are
hundreds of structurally related natural and synthetic cardenolide
analogues of known crystal structure or stereochemistry (reviewed in (10) ), permitting studies of fine specificity. In addition,
anti-digoxin antibodies have relatively high affinity compared to
antibodies against smaller haptens.
Figure 1: Schematic representation of the digoxin (digoxigenin tridigitoxoside) numbering system and steroid ring nomenclature. The digoxin analogues used in panning are substituted at position 16 in the D ring (see Table 4).
The three-dimensional structure of the high affinity anti-digoxin 26-10 Fab in the uncomplexed state and complexed with digoxin has been determined(11) . Neither the antibody nor digoxin undergoes any significant conformational change upon forming the complex. Neither hydrogen bonds nor salt bridges are formed between 26-10 and digoxin. High affinity binding occurs solely by shape complementarity involving extensive nonpolar interactions between antibody and the hydrophobic hapten.
We used in vitro selection to isolate variants of
26-10 with altered affinity or hapten binding
specificity(12, 13) . Among a set of 11 independent
spontaneous mutants of 26-10, mutation at several H chain
residues affected affinity. Some of these proved to be hapten contact
residues, as shown by the crystal structure(11) . Among these
spontaneous variants, recurrent mutations at H chain residue 35 were
observed(13) . ()The complementarity between the
26-10 Asn-H35 residue and digoxin was analyzed further by
constructing mutants using site-directed mutagenesis in a hybridoma
system(13) . In order to advance our understanding of the
complementarity between Ab 26-10 and digoxin, we turned to a
phage display system.
The demonstration that both Fv (14) and Fab (15) can be produced correctly folded in bacteria by expression in the periplasmic space was a critical step, which has accelerated experimental antibody-combining site engineering. Display of diverse libraries of V regions as Fab or single-chain Fv (sFv) on the surface of filamentous bacteriophage has permitted selection of rare antibodies with altered binding affinity (16, 17, 18) . Here we use a phage-displayed 26-10 Fab library, randomized in the region of H chain CDR1 and selection with digoxin and C-16-substituted digoxin analogues, to probe the 26-10 combining site.
The 26-10 hybridoma cell line
(2a,
) was generated from the fusion of Sp2/0 cells with
spleen cells of an A/J mouse immunized with digoxin-human serum
albumin(21) . The 26-10 H and L chain cDNAs were prepared
from total RNA using standard methods(22) .
Oligonucleotide primers (1-4, Table 1and Fig. 2) were used in a polymerase chain reaction (Taq polymerase from Perkin Elmer, Norwalk, CT) with the 26-10 cDNA template to modify the ends of 26-10 DNA for insertion into the pComb3 vector (Fig. 2A). Primer 1 introduces a MunI site into the first two codons of 26-10 H chain DNA without altering the amino acid sequence. Primer 2 inserts an SpeI site 3` to nucleotides encoding Arg-228 in the hinge region of the 26-10 H chain. Primer 3 introduces an AatII site in the 5` terminus of the L chain. Primer 4 contains stop codons and an XbaI site inserted immediately 3` of Cys-L214, the final L chain codon. Mutagenic oligonucleotides were synthesized using an Applied Biosystems 380B DNA synthesizer.
Figure 2: Diagram of the construction of the vector pComb3-26-10-1 and generation of randomly mutated residues H30-35. A, the portion of the pComb3 bacterial expression vector (20) containing the lacZ promoters, pelB leader sequences, insertion sites for the Fd portion of the H chain and the L chain, and the GGGGS linker at the M13/coat protein III (cpIII) fusion site. The MunI and AatII restriction enzyme sites were introduced to provide sites for insertion of 26-10 H and L chains that would not alter the wild type amino acid sequence. The resulting first four residues of the H chain (EVQL) and the first two of the L chain (DV) and Ala(A) encoded by the pelB leader are indicated with the single letter codes. Numbered arrows refer to the oligonucleotides used in these manipulations (see Table 1). B, the 26-10 H chain V region with the oligonucleotides used for alteration of the polymerase chain reaction template and randomization of amino acids H30 to H35. Nucleotide numbering is in italics.
The expression vector pComb3 was altered using primers 5-8 (Table 1), complementary to and overlapping those used to modify 26-10, in order to retain the precise N-terminal amino acid sequences of 26-10 H chain (EVQL) and L chain (DV). Two new unique restriction enzyme sites (MunI and AatII) were inserted into pComb3 as a result of these manipulations (Fig. 2A).
The vector and the antibody encoding DNA
were each cut with MunI, SpeI, AatII, and XbaI (New England Biolabs, Beverly, MA) (Fig. 2A). The resulting fragments were purified from
agarose gels by electroelution (Stratagene, La Jolla, CA), ligated, and
transformed into XL1-Blue Escherichia coli. Plasmid clones
were screened for the presence of the desired pComb3 and 26-10
restriction enzyme sites. The sequence of the 26-10 Fab and
altered pComb3 coding portions of the resulting construct, designated
pComb3-26-10-1, was verified using dideoxy sequencing
(Sequenase, U.S. Biochemical, Corp.) and
[-
S]thioadenosine 5`-triphosphate (DuPont
NEN).
The vector designated pComb3-26-10-2 is a template with HCDR1 alterations designed to prevent parental 26-10 contamination of 26-10 libraries randomized in HCDR1 and consequent competition in biopanning. The vector pComb3-26-10-2 has an altered HCDR1 created by oligonucleotide-directed mutagenesis employing the method of Eckstein (23) (Kit RPN1523, Amersham Corp.). The mutations insert a unique AgeI site, a Pro-H34 (CCG), and stop codon (TAA) in H chain CDR1 using oligonucleotide primer 9 ( Table 1and Fig. 2B). The new restriction enzyme site allows linearization and elimination of unwanted plasmids, while the amino acid changes render the 26-10 protein product nonparental and nonfunctional.
Random mutations were introduced into six codons of the H chain V region of 26-10, including ``framework'' position 30 and CDR1 positions 31-35 using primers 10 and 1 (Table 1). Primer 10 is antisense 5`-CTT TCC ATG GCT CTG CCT GAC CCA SNN SNN SNN SNN SNN SNN GAA TAT GTA TCC AGA AGA CTT G-3`, where N represents an equimolar mixture of all four dNTPs and S represents an equimolar mixture of dCTP and dGTP. The NcoI site is underlined. Primer 1 contains the 5` terminus of the 26-10 H chain and the MunI site (Table 1, Fig. 2B). The 147-bp product of this reaction encodes the 5` V region of 26-10 H chain through CDR1 and the adjacent NcoI site at nucleotide position 504 bp.
This 147-bp fragment was digested with NcoI and MunI, and the resulting 114-bp fragment was dialyzed and concentrated using a Centricon 30 filter unit (Amicon, Beverly, MA), ligated to a NcoI-MunI-digested pComb3-26-10-2 vector, and introduced by electroporation into E. coli XL1-Blue F` cells. Infection of the XL1-Blue library of CDR1 mutants with VCSM13 helper phage generated a library of phage with surface Fab (XL1-Blue F` and VCSM13 are from Stratagene, La Jolla, CA). Phage were recovered and concentrated by polyethylene glycol/NaCl precipitation from bacterial supernatants(20) . Bacteriophage yield was quantitated by titration on lawns of XL1-Blue F` bacteria, and phagemid was quantitated by postinfection XL1-Blue F` colony formation on LB (1.0% tryptone, 0.5% yeast extract, 1.0% NaCl, pH 7.0)/agar/carbenicillin plates(22, 24) .
Phage were analyzed after four successive rounds of biopanning. Phagemid was grown without VCSM13 helper, and DNA was isolated from a 50-ml culture using a Plasmid Midi Kit (Qiagen, Chatsworth, CA). An aliquot was sequenced to estimate the complexity of the sequences in the mutated CDR1. A 1-µg sample was cut with NheI and SpeI and religated (see Fig. 2A) to remove the gene III sequence that serves as a membrane anchor(20, 24) . Fabs secreted from cells containing these phagemids terminate at Arg-228 in the hinge region of 26-10 and contain a gene III-encoded Thr and Ser just 5` to the stop codon. After religation, the DNA was transformed into XL1-Blue cells and plated on LB/carbenicillin plates. Bacterial colonies were isolated and screened for Fab production and antigen binding.
The specificity of Fab for different
cardiac glycosides was determined using a competition radioimmunoassay
based on the affinity assay described
above(12, 13, 29) , with the addition of 0.5
µg/ml goat anti-mouse Fab antibody to each sample to retain Fab on
the filter. Stock solutions of cardiac glycosides
(10M in pyridine) were diluted in PBSA to
result in 10
to 10
M in
the final mixture. Each mutant Fab was competed with digoxin as an
internal control. The values reported are ratios of molar
concentrations of inhibitor required to give 50% inhibition relative to
the molar concentration of digoxin that gave 50% inhibition.
Gln
derived from the native codon or from the suppressed TAG stop codon in
XL1-Blue E. coli was included in the same group.
is the difference between the observed frequency of an amino acid and
its expected frequency as standard deviation according to the following
formula.
On-line formulae not verified for accuracy
where Obs(X) is the number of occurrences of amino acid X in the selected sequences, P(X) is the expected
frequency of amino acid X occurring (1/32, 2/32, or 3/32), and n is the total number of sequences. Values greater than
2 are considered significant(30) .
DNA restriction fragments containing sequences encoding the murine anti-digoxin antibody 26-10 Fab were introduced into the vector pComb3 (20) , which had been modified to reproduce the 5` amino acids of 26-10 H and L chains (Fig. 2). The complete DNA sequences of the 26-10 insertion and the pComb3 boundaries were verified as correct. When transformed into E. coli XL1-Blue cells, the plasmid pComb3-26-10-1 produced 26-10 Fab on the surface of bacteriophage M13 as a gene III fusion protein(31) , detected by phage ELISA.
Bacterial 26-10 Fab was characterized by comparison with 26-10 Fab prepared enzymatically from hybridoma protein. The secreted bacterial 26-10 Fab was produced by excision of the gene III coding fragment of the phagemid by NheI/SpeI digestion (Fig. 2) (20) , followed by religation and transformation into E. coli XL1-Blue. Bacterial Fab was indistinguishable from enzymatically prepared 26-10 Fab by Western blot using isotype-specific reagents and by ELISA utilizing glycoside-BSA conjugates (data not shown). Because periplasmic space preparations and secreted preparations both yielded about 1 mg/liter of 26-10 Fab, the more easily prepared secreted samples were used thereafter.
Bacterial 26-10 Fab was purified on a ouabain-amine-Sepharose column (32) . This preparation was indistinguishable from enzymatically prepared 26-10 Fab by SDS-polyacrylamide gel electrophoresis. Amino acid sequence analysis of bacterial 26-10 Fab resulted in a mixture of H and L chain amino acids in a 1:1 ratio with residues identical to the bona fide hybridoma antibody 26-10(12) . This result indicated that cleavage of the pelB leader did not result in a product of different length and that the N-terminal sequences of 26-10 H chain and L chain in pComb3 had been correctly modified.
In
preliminary biopanning experiments the phagemid library proved
contaminated with unmutated phagemid having the same nucleotide
sequence as wild-type 26-10. In order to avoid contamination, a
nonfunctional template designated pComb3-26-10-2 was
designed. pComb3-26-10-2 was used to create a library of
mutants randomized at H chain positions 30-35 inclusive. We
randomized the H chain CDR1 amino acids 31-35 as well as the
nominal framework residue 30. Randomization was achieved using nucleic
acid mixtures containing all four dNTPs(N) or dCTP and dGTP (S), in the
order NNS at the three codon positions of the H chain 30-35 amino
acids. This results in a complexity of 6.4 10
amino
acids from 1.1
10
different codon combinations.
Phage aliquots were panned against four cardiac glycosides differing
at position 16 (Fig. 1) on the cardenolide moiety (16-H,
digoxin; 16-OH, gitoxin; 16-OCHO, formylgitoxin or gitaloxin; and
16-OCOCH, 16-acetylgitoxin) (Table 4). These
analogues were chosen because HCDR1 residues are close to position 16
in the co-crystal structure of 26-10-digoxin (Fig. 3) (11) .
Figure 3:
Statistical distribution of amino acids at
H chain positions 30-35 among a set of unique clones recovered
after panning using four different cardiac glycosides. The
corresponding amino acid sequences are listed in Table 2.
Position numbers are at the upper right for each panel. is the standard deviation of the difference
between the observed frequency and the expected
frequency(30) .
Initially, the library was panned against digoxin-BSA and 16-acetylgitoxin-BSA. Following four rounds of panning, 20 clones for each cardiac glycoside were isolated, and the secreted products were tested for binding to digoxin, 16-acetylgitoxin, and goat anti-mouse Fab (Table 2). Using an ELISA, all digoxin-selected clones secreted Fab and bound digoxin-BSA. All digoxin-selected mutants except mutants D4-3 (see Table 2, Footnote a, for nomenclature) and D4-12 bound 16-acetylgitoxin. Nineteen of 20 of the 16-acetylgitoxin-selected clones contained Fab and bound both cardiac glycoside conjugates.
DNA sequence was obtained for the entire H chain V region of each mutant Fab. There were no sequence alterations outside the mutated region, and all selected clones had G or C in the third codon position. After four rounds of panning, 15 unique sequences were present in the 20 clones selected by digoxin-BSA and only two unique sequences in the 19 clones selected by 16-acetylgitoxin. In every instance, clones selected more than once had identical nucleic acid sequences.
After four rounds of panning, the Asn(N) residue at position H35 was present in 13 of 15 unique sequences selected by digoxin-BSA and both sequences selected by 16-acetylgitoxin (Table 2). The exceptions were Val-H35 in clone D4-3 and Ser-H35 in clone D4-15. Clone D4-3 differs from clone D4-11 only at position H35. Clone D4-15 has Asn at position H34, while other digoxin selected clones have Phe, Ile, Val, or Tyr at position H34. The wild-type Met-H34 was not found in any mutant. Position H33 was predominantly Phe, Tyr, or Trp. Positions H30-32 demonstrated a less restricted set of amino acids. At position H32 we observed Phe, Gln, Arg, Ser, Thr, Gly, Met, Ala, or Leu. Positions H30 and H31 contained the most varied sets of amino acids found in the mutated segment, although there is some preference for the non-wild-type Pro at position H30 and the wild-type Asp at H31.
The parental 26-10 sequence was not panned from this library, although the sequence chosen five times by digoxin-BSA selection (RDFYYN) differs by only two residues from 26-10. The wild-type residues Thr-H30 and Met-H34 do not appear in any of the mutants selected by these two cardiac glycosides.
Another aliquot of the same library was grown and selected against gitoxin-BSA and 16-formylgitoxin-BSA. Four rounds of panning produced 6/20 clones ELISA-positive for gitoxin binding and 0/20-positive for 16-formylgitoxin binding. All were producing approximately 1 µg of Fab/ml of media. Two additional rounds of panning were carried out, and 10 clones were analyzed for each analogue from pan 6. All pan 6 gitoxin-BSA-selected clones were positive for Fab production and bound both gitoxin and digoxin. Eight of 10 16-formylgitoxin-BSA clones selected from pan 6 were positive for Fab production, 16-formylgitoxin binding, and digoxin binding.
DNA sequence analysis of the gitoxin-BSA and 16-formylgitoxin-BSA pan 6 clones revealed a limited set of sequences. Seven of 10 of the gitoxin-selected clones were RDFYYN, two contained GERFFN, and one was SKRYIN. Of the 16-formylgitoxin clones sequenced, six of eight had the sequence SHSYIN and two were TRYWFN. Using 16-formylgitoxin-BSA, we did not observe the mutant RDFYYN, which dominated the output of the final pans in the experiments using the other three cardiac glycosides, even though the phage preparation used was an aliquot of the same library.
To seek a more diverse sample of positive clones, 10 clones from pan 5 on gitoxin-BSA and 16-formylgitoxin-BSA were analyzed. All 10 gitoxin-BSA-selected clones were ELISA-positive for gitoxin and digoxin, while six of 10 16-formylgitoxin-BSA clones were positive for 16-formylgitoxin and digoxin binding (Table 2). The corresponding amino acid sequences are shown in Table 2. Pan 5 contained a more diverse group of clones than pan 6. Two of the gitoxin-selected clones from pan 6 were also selected in pan 5 (shown by asterisk), and four additional pan 5 sequences were found. All antigen-binding clones had Asn-H35. The distributions of amino acids at positions 33 and 34 were similar to those in the digoxin-BSA-selected group. Clones from the 16-formylgitoxin pan 5 group contained both clones seen in pan 6 and two additional unique clones. Shown for comparison in the table are sequences from clones that did not bind either digoxin or 16-formylgitoxin. Two of these (mutants F5-1 and F5-7) had a residue other than Asn at position H35, and two others (mutants F5-2 and F5-8) had a pair of adjacent prolines.
The most
prevalent amino acids at each position in unique clones selected by all
four analogues shown in Table 2are (H30-35)
PD(F/Y)(F/W/Y)(I/F/Y)(N). This ``consensus'' is based upon
the statistical distribution of amino acids at each substituted
position (Fig. 3), and includes any amino acid that appears with
a frequency that differs from that expected by 2
at any
position(30) .
Affinities for digoxin of hybridoma and
bacterial 26-10 Fab and certain mutants from pans 4 and 6 (Table 3) were measured. Bacterial 26-10 Fab had the same
affinity (5.3 10
M
) as
enzymatically prepared 26-10 Fab (5.4
10
M
) in this assay. Mutant Fabs had
affinities ranging from 6.0
10
M
to 2.2
10
M
. The mutants (A4-19 and
A4-20, Table 3) with the highest affinities for digoxin
(2.2
10
M
and 9.7
10
M
) bore the
sequences PSFYYN and RDFYYN, respectively. RDFYYN is the clone that was
prevalent among the sequences in the final pan of all analogues except
16-formylgitoxin (Table 2). The highest affinity clone
(A4-19) was rescued only once, from panning employing
16-acetylgitoxin. All clones with affinities in the range of that of
26-10 had Asn-H35. Clone D4-3 had an affinity reduced
5000-fold compared with parental and a Val at position H35. Clone
D4-11 had an identical sequence to D4-3, except that
D4-11 had Asn at position H35 and an affinity for digoxin of 6.6
10
M
. Clone D4-15
contained Ser at position H35 as well as four other differences from
parental, yet had an affinity reduced only 3-fold. Clone G6-2 had
reduced affinity (3.0
10
M
) and was the only clone containing Lys.
Clone D4-4 had Ser at position H33 instead of Tyr compared with
the parental and had reduced affinity for digoxin (9
10
M
). Among 20 different
mutant clones analyzed, 17 had affinities for digoxin equal to or
greater than 1.8
10
M
.
All mutants with Tyr-H34 had affinities for digoxin equal to or greater
than parental 26-10.
The specificity of mutant Fab for different cardiac glycosides was compared using a competition assay (Table 4). Binding to the other cardiac glycosides was not enhanced compared with digoxin. Binding of mutants F6-1 and F6-2 to gitoxin was reduced 4-9-fold compared with wild-type 26-10, and binding of mutant G6-2 to 16-acetylgitoxin was reduced 30-fold.
To determine whether there was another characteristic besides relative avidity for cardiac glycoside-BSA that served as a basis for mutant selection, a series of mixing experiments were performed (Table 5). In the first experiment, A4-20 phage was mixed with wild-type 26-10. Following a single round of panning, the only sequence detectable was that of A4-20 (data not shown). In a second experiment, we mixed clones A4-19, G6-1, G6-2, F6-1, F6-2, and 26-10. After three rounds of panning, the batch sequences from digoxin-BSA, 16-acetylgitoxin-BSA, and 16-formylgitoxin-BSA pannings were exclusively A4-19. The only exception was the sample panned against gitoxin-BSA, which showed a mixture of G6-2 and A4-19. The unselected control showed no predominant sequence (data not shown).
In a third experiment, G6-1, G6-2, F6-1, F6-2, and 26-10 were mixed. DNA sequences from the pooled output of the third pan against each glycoside revealed selection bias, prompting sequencing of 10 individual clones from each analogue (Table 5). Nine of 10 of the unselected clones were G6-1, suggesting that a growth advantage exists for G6-1. Seven of 10 clones selected by gitoxin-BSA and 16-formylgitoxin-BSA were G6-2 despite the low affinity for digoxin of G6-2 (Table 3). Clones selected by digoxin and 16-acetylgitoxin showed no predominant sequence. Parental 26-10 was observed only twice among 10 sequences selected by digoxin-BSA. A4-19 and A4-20 dominated in mixing experiments and had the highest affinities for digoxin. However, certain clones, such as G6-2, were able to outcompete the others in panning against specific analogues.
To
determine panning efficiency, the number of colonies bound to the
glycoside-BSA well minus the number bound to BSA alone (output) divided
by the number of colonies in the aliquot of input bacteriophage was
measured. The 26-10 phage and the negative control 36-65
phage bound gitoxin, 16-acetylgitoxin, and 16-formylgitoxin conjugates
nearly equal to BSA alone. Phage 26-10 bound digoxin-BSA with a
panning efficiency of 2.8 10
. In contrast,
A4-3, A4-19, G6-1, G6-2, F6-1, and
F6-2 all showed panning efficiencies for digoxin at least an
order of magnitude higher than 26-10 phage. Panning efficiencies
of the mutants for the other glycoside conjugates were generally
similar for digoxin and gitoxin and lower for 16-acetylgitoxin and
16-formylgitoxin but did not provide a basis for preferential selection
of one mutant over another by a particular glycoside conjugate.
We previously demonstrated that single amino acid changes in crystallographically defined contact and noncontact residues of the high affinity anti-digoxin monoclonal antibody 26-10 can cause profound changes in affinity and alterations in specificity for digoxin and its analogues. The primary structures and binding characteristics of a set of spontaneous variants and engineered mutants of 26-10 were determined(12, 13, 29, 33) . The importance of the identity of position H35 (Asn) to digoxin binding was indicated by these studies, in which all substitutions for Asn at position H35 reduced affinity(13) . In the structure of the complex of digoxin and 26-10 there is close complementarity between the side chain of Asn-H35 and digoxin(11) . Tyr-H33 is also a digoxin contact residue (Fig. 4).
Figure 4: Stereo views of H-CDR1 and digoxin from the x-ray crystal structure of the 26-10 Fab-digoxin complex(11) . The digoxin molecule is on the left with the lactone at the bottom. At the top is shown a single attached sugar (digoxigenin monodigitoxoside). Tyr-H33 and Asn-H35 of 26-10 CDR1 make van der Waals contact with the hapten.
In order to examine a wider variety of mutants that potentially alter digoxin binding, we used phage display of 26-10 Fab and randomized H chain residues 30-35 including H chain CDR1 (residues 31-35). By enriching and selecting mutants by panning on antigen-coated surfaces, we attempted to investigate several questions. To what degree do structural constraints on digoxin binding limit sequence diversity in the CDR loop? Is there more than one solution for sequences in this segment consistent with high affinity digoxin binding? Are there mutants with enhanced affinity or altered specificity relative to the wild-type 26-10?
Bacteriophage libraries were selected by binding to digoxin-protein conjugates and analogues of digoxin conjugated to proteins, to probe the structural requirements for antigen binding. We chose three digoxin analogues that bear substitutions at the cardenolide C-16 position (see Fig. 1) of varying size, with corresponding effects on binding to 26-10 as compared with digoxin(12, 13, 29) . We displayed Fab rather than sFv, as reduced binding may be observed for certain sFv constructs, which perforce contain a synthetic interdomain peptide linker(34) .
The pComb3 vector (20, 24) dictates certain amino acids at the N terminus of both H and L chains when DNA coding for Fab is inserted at the existing restriction sites. These sequences differ from those of 26-10 nucleotides encoding the first four amino acids of the H chain and first two of the L chain (Fig. 2). We previously demonstrated that H chain N-terminal sequence differences of anti-digoxin antibodies cause reduction in affinity for digoxin(35, 36) . These results were confirmed by mutagenesis in which short truncations and point mutations of the N terminus altered affinity for digoxin. A single engineered mutation in the H chain signal peptide sequence at the -2 position of 26-10 resulted in an antibody with three additional residues at the N terminus associated with a 100-fold reduction in affinity for digoxin. In order to avoid the possible confounding effect of N-terminal sequence differences on antigen binding of 26-10 and 26-10 mutants, we modified pComb3 to code for the 26-10 H and L N-terminal sequences. Because 26-10 H and L chains were synthesized with bacterial leader sequences ordinarily associated with specific N-terminal H and L chain sequence contexts and because the N-terminal amino acids encoded by the original vector were changed to correspond to those in 26-10, it was important to verify proper cleavage of the leader sequence(37, 38) . Amino acid sequence analysis of the bacterial Fab indicated that the leader sequence was cleaved correctly.
The bacterial 26-10 Fab was indistinguishable from enzymatically prepared Fab on the basis of SDS-polyacrylamide gel electrophoresis, Western blots, ELISA for digoxin binding, and affinity determination (Table 3). The 26-10 Fab affinity values reported here using a secondary antibody for Fab immobilization are slightly lower than those reported for 26-10 Ig in the absence of a secondary antibody(33) .
The H-CDR1 loop of antibody 26-10 contains residues at positions H26, H27, H29, and H34, also identified in a canonical H-CDR1 loop described by Chothia et al.(8) that are thought to be important for packing H-CDR1 against the framework region. However, 26-10 lacks Arg at position H94, which is presumed to be important in maintaining the H-CDR1 loop structure. In addition, such a canonical loop conformation would be incompatible with the structure of 26-10 (11) because of interactions of H-CDR1 with residues in H-CDR2 and H-CDR3. Pro-H52 in H-CDR2 would clash with the conserved aromatic at H29 in the ``canonical'' structure, Tyr-H53 would clash with the side chain of Thr-H30, and Ser-H96 in H-CDR3 would clash with the side chain of Tyr-H32 (Phe in 26-10).
There is no evidence from the electron density maps of a well
defined conformation for residues H27-30, although the possible
backbone conformations are somewhat restricted by the end point
connections at H26 and H31. This disorder is also evident in the
structure of the uncomplexed 26-10 Fab (11) as well as in
the structure of the 26-10 nonbinding mutant R9. ()The
crystallographic disorder indicates either that there are several
different conformations for that region or that the region is highly
mobile.
We analyzed 26 unique mutants obtained from the library by
panning against digoxin and three analogues (Table 2). All
mutants specifically bound digoxin in an ELISA, and were therefore
included in a statistical analysis of amino acids at each position (Fig. 3). As 26 unique clones constitute a relatively small
sample size, the absence of a particular amino acid at any position is
not statistically significant. One hundred fifty different clones are
the minimum needed for the absence of an amino acid designated by 1 in
32 codons to be significant (greater than 2). On the other
hand, the frequent presence of a given amino acid at a particular
position can be significant. For a sample including 26 clones, an amino
acid must be present at least three, five, or six times to differ from
the mean by greater than 2
, depending on whether it is
encoded by one, two, or three codons, respectively.
A distinct pattern of amino acid usage emerged from the statistical distribution of amino acids at positions H30-35 (Fig. 3). The most dramatic was the conservation at position H35 of Asn, which contacts digoxin in the crystal structure. This result was not entirely unexpected, as in earlier experiments employing site-directed mutagenesis at position H35, all substitutions tested had significantly lower affinity for digoxin, indicating that Asn was essential for optimal complementarity at this position(13) .
In the
co-crystal structure of 26-10 with digoxin(11) , the side
chain of Asn-H35 contacts the digoxin D ring atoms C-16 and C-17 and
the lactone (atoms C-20, C-21, and C-22) (Fig. 1). In addition,
Asn-H35 forms hydrogen bonds with Tyr-H47 and Ser-H95, each of which
also contacts hapten. Asn was found at position H35 in every high
affinity clone isolated from the phage libraries. Clone D4-3
(Val-H35) had an affinity of 6 10
M
, and clone D4-15 (Ser-H35) had an
affinity for digoxin of 1.8
10
M
. Clone D4-15 is the only mutant
with Asn at position H34. Asn-H34 cannot functionally replace Asn-H35
and maintain loop integrity, but Ser-H35 could partially fulfill the
role of Asn-H35 by making a hydrogen bond with either Ser-H95 or
Tyr-H47.
At position H34, the wild-type residue methionine was not observed in any mutant. The predominant residues were Ile and the aromatic residues Phe and Tyr. In the wild-type structure, the side chain at residue H34 points away from the binding cavity and packs into the interior. Aromatic substitution at this position (particularly Tyr) would clash with Pro-H52 and Trp-H32, resulting in some rearrangement of the local structure. However, all four mutants containing Tyr-H34 (Table 3) had affinities equal to or greater than parental 26-10, suggesting that such rearrangement can enhance complementarity. The phage-selected consensus sequence (see below) differs from the parental 26-10 (TDFYMN) only at positions 34 and 30. Met-H34 and Thr-H30 occur in more than 50% of murine H chains (39) and may reflect conserved germ line residues. Three base changes are required to convert Met-H34 (ATG) to either codon for Tyr (TAT, TAC). The affinity maturation process in vivo that results in high affinity Abs such as 26-10 is constrained by the germ line V region gene sequences, which are subjected to random somatic mutation and antigen-driven selection. Such putative constraints are not operative during in vitro selection of phage-displayed mutants (such as A4-19 and A-20), which can display even higher affinities for digoxin. Enhancement of affinity of high affinity Abs using phage display selection was previously reported(40) , as was improvement in affinity using site-directed mutagenesis based on crystal structures of secondary response mAb (41) .
The other contact residue in CDR1, Tyr-H33 (Fig. 4), is predominantly selected among the mutants (Fig. 3), but can be replaced with Phe or Trp, consistent with the importance of a planar hydrophobic (aromatic) residue at this position. The side chain of Tyr-H33 contacts the C-3, C-4, C-7, and C-15 atoms of digoxin.
Amino acid residues at positions H30-32 demonstrate much greater sequence variability than those at positions H33-35 (Table 2). The statistically significant predominant amino acids are Pro-H30, Asp-H31, and Phe- or Tyr-H32 (Fig. 3). However, 12 different amino acids are found at position H30, 14 at position H31, and 10 at position H32, all of which are permissible with retention of digoxin binding. Residues H27-32 exhibit disorder in the crystal structure of 26-10-digoxin (11) . This disorder is paralleled by the promiscuity in amino acid sequence usage at residues H30-32 among the mutants. For most Fabs, the side chains at positions H30 and H31 are solvent-exposed and have wide variations in their conformation. Since this part of H-CDR1 is not well-defined in the 26-10 crystal structure, one might not expect substitution at these positions to have much effect on affinity. Nonetheless, the observation that mutant G6-2 (SKRYIN) has an affinity for digoxin reduced by more than an order of magnitude compared with mutant D4-2 (PGRYIN), which differs only at positions H30 and H31, suggests that there are restrictions upon certain substitutions or sets of substitutions at these positions. Position H32 is partially solvent-exposed, with side chain conformations clustered in other Fab structures. In 26-10, the wild-type Phe and also Tyr occur frequently, but many other substitutions are allowed.
Although we found a consensus sequence for positions H30-35 (PD(F/Y)(F/W/Y)(F/I/Y)N) (Fig. 3), no single mutant Fab exhibits the entire consensus sequence (Table 2). Several mutants differ from the consensus by only a single amino acid within positions 30-32 (mutants D4-8, D4-17, A4-19, and A4-20). Such a consensus sequence accounts only for the prevalence of amino acids at individual positions, independent of possible effects of each upon the other.
Among digoxin-binding clones, proline was not observed at positions H32-35 inclusive, yet it occurs nine times at positions H30 and H31 (Table 2). Proline residues restrict backbone conformation, and may disturb backbone folding necessary for complementarity at positions H33-35. The permissiveness at residues H30 and H31 is illustrated by mutant D4-6, which binds digoxin despite proline at both positions. However, a double proline present in the mutants F5-2 at H33-34 and F5-8 at H31-32 (Table 2) is associated with failure to bind digoxin, despite the presence of Asn at H35.
The data indicate that H-CDR1 is more
plastic than we anticipated. Although the diversity of side chain use
is consistent with the disordered crystal structure in this region of
26-10, it is remarkable that there are many amino acid
combinations that result in Fabs with similar affinities for digoxin.
Two Fabs (A19 and A20) have affinities greater than 26-10 (4.1-
and 1.8-fold, respectively), and 14 of 20 mutants (Table 2) have
affinities 3.4
10
M
. Thus, the variability at positions
H30-32 does not prevent digoxin binding but must provide a stable
framework for the contacting residues H33-35, which are the main
arbiters of affinity and specificity in this region.
We wondered
whether the absence of the wild-type 26-10 sequence among
selected mutants could be accounted for through considerations of
library size, preferential growth of certain mutants, or difference in
panning efficiency. Not all amino acid sequences are represented
equally in an NNS-substituted library, with the ratio of various amino
acids being 1:2:3, depending on codon usage. In a library with six
randomized positions, the likelihood of finding an amino acid encoded
by 1 of 32 codons at all positions is 1 in 1.1 10
.
Wild-type 26-10 would be expected to occur once in 5.4
10
clones. Because the library contains 1
10
different nucleotide sequences, it is not surprising that
26-10 was not observed. The probability of identifying
26-10 among the panned mutants may be further reduced by lack of
selection due to inefficient expression of 26-10 on phage as a
fusion protein or uneven clone distribution due to growth advantages of
individual clones. Mixing experiments revealed that 26-10 does
not compete as well with some other mutants with similar or even lower
affinities for digoxin in competitive binding to antigen-coated
surfaces (Table 5), although 26-10 was readily enriched
from mixtures (1:1000) with an anti-tetanus clone. Fab 26-10
showed a panning efficiency at least an order of magnitude lower than
the other tested phage, sufficient to account for its absence after
panning against digoxin-BSA. These observations raise the possibility
that clones of equal or higher affinity for cardiac glycosides may not
have been recruited from libraries panned against cardiac glycoside-BSA
conjugates.
It is also possible that carrier effects during panning are operative in selecting certain mutants. However, as digoxin is coupled to BSA through the terminal digitoxose moiety, and the crystal structures of the 26-10-digoxin complex have shown that the sugars are not involved in binding, a role for carrier in affecting the panning results is unlikely.
The selection of clone G6-2 from
the bacteriophage library, despite its low affinity, is reproduced in
the mixing experiment (Table 5), and therefore does not directly
reflect the affinity of this mutant for digoxin or gitoxin. Clone
G6-2 predominates in selection from the mixture of five clones,
and clone G6-1 appears to display a distinct growth advantage
based on its predominance in the unpanned control. The specificity
experiments (Table 4) indicate that G6-1 and G6-2
Fab, both isolated following panning against gitoxin-BSA, bind gitoxin
in the same range as native 26-10, A4-19, and A4-20,
but better than the 16-formylgitoxin-selected clones. A4-19 has
higher affinity for digoxin (Table 3) and 16-acetylgitoxin (Table 4) than A4-20, yet it was not among the clones
selected following panning against digoxin-BSA and occurs in only one
of 19 clones selected in pans against 16-acetylgitoxin (while
A4-20 accounts for 18 of 19). Furthermore, A4-19
predominated in mixing experiments. Perhaps A4-19 may not have
been selected as frequently as A4-20 due to underrepresentation
in the original library. The close complementarity between antibody
26-10 and digoxin in the region of the 16 position in the D ring
of digoxin (Fig. 4) results in 26-10 distinguishing
substituents of different sizes at C-16. Digoxin C-16 contacts H-CDR1
as well as H-CDR3(11) . The 26-10 antibody binds gitoxin
(16--OH) 3-5-fold more weakly than digoxin. Larger
substitutions at C-16 would cause more disruption of the combining
site, and analogues with formyl and acetyl groups at C-16 were shown to
demonstrate further reductions in affinity for wild-type 26-10 (Table 4)(29) .
All 26-10 mutants bound digoxin best and bound the other three analogues in the same decreasing order: gitoxin, 16-formylgitoxin, and 16-acetylgitoxin (Table 4). Three of six mutants tested demonstrated specificity shifts. Clone G6-2 binds 16-acetylgitoxin 30-fold more weakly than digoxin as compared with wild-type 26-10. However, Fab G6-2 has reduced affinity for digoxin. Clones F6-1 and F6-2 selected by 16-formylgitoxin had moderate reduction in relative affinity for gitoxin. The affinity data (Table 2) indicating the frequent occurrence of high affinity clones of diverse sequences, taken together with the modest specificity shifts, indicate that, at least in experiments employing random mutagenesis limited to H-CDR1, no mutants with significant improvement in relative binding for C-16-substituted analogues can be obtained without concomitant drastic reduction of affinity for hapten. Thus, antibody 26-10 is relatively ``optimized'' in this region to bind digoxin. However, the results for a few mutants (A419 and A420) indicate that this affinity can be improved.