A high-affinity human TCR-like antibody detects celiac disease gluten peptideтАУMHC complexes and inhibits T cell activation
TCR-like antibodies tackle celiac disease
Ingestion of gluten-containing food triggers the gastrointestinal symptoms of celiac disease in patients with CD4+ T cells specific for deamidated gluten peptides presented by disease-associated HLA-DQ class II MHC molecules. Frick et al. used phage display technology to screen for TCR-like antibodies specific for an immunodominant gluten peptide bound by HLA-DQ2.5. Antibody engineering optimized affinity and binding stability, yielding an improved TCR-like antibody that structurally mimicked the TCR interface with gluten peptideтАУMHC complexes. These TCR-like antibodies blocked activation and proliferation of gluten-responsive human CD4+ T cells both in vitro and in DQ2.5 transgenic mice. TCR-like antibodies that block immunodominant epitope recognition have potential as personalized medicine treatments for blunting gluten-activated T cell responses without compromising effector functions provided by other T cells.
Abstract
Antibodies specific for peptides bound to human leukocyte antigen (HLA) molecules are valuable tools for studies of antigen presentation and may have therapeutic potential. Here, we generated human T cell receptor (TCR)тАУlike antibodies toward the immunodominant signature gluten epitope DQ2.5-glia-╬▒2 in celiac disease (CeD). Phage display selection combined with secondary targeted engineering was used to obtain highly specific antibodies with picomolar affinity. The crystal structure of a Fab fragment of the lead antibody 3.C11 in complex with HLA-DQ2.5:DQ2.5-glia-╬▒2 revealed a binding geometry and interaction mode highly similar to prototypic TCRs specific for the same complex. Assessment of CeD biopsy material confirmed disease specificity and reinforced the notion that abundant plasma cells present antigen in the inflamed CeD gut. Furthermore, 3.C11 specifically inhibited activation and proliferation of gluten-specific CD4+ T cells in vitro and in HLA-DQ2.5 humanized mice, suggesting a potential for targeted intervention without compromising systemic immunity.
INTRODUCTION
Antibodies (Abs) with specificity for peptide-MHC (pMHC), so-called T cell receptor (TCR)тАУlike Abs, have been successfully used to detect and quantify peptide presentation on cells, and several studies suggest therapeutic potential by different modes of action, including inhibition of pathogenic T cells and killing mechanisms to delete antigen-presenting cells (APCs) (1тАУ3). As soluble reagents, Abs are often preferred over TCRs because of their increased stability and higher affinity (4тАУ8). We have previously and in this study generated TCR-like Abs specific for pMHC complexes implicated in celiac disease (CeD) (9).
CeD is an inflammatory autoimmune-like condition of the small intestine caused by immune reactions to dietary gluten proteins (10). The disease is driven by CD4+ T cells that recognize deamidated gluten peptides in the context of the disease-associated human leukocyte antigen (HLA)тАУDQ molecules HLA-DQ2.5 (DQA1*05/DQB1*02), HLA-DQ2.2 (DQA1*02:01/DQB1*02:02), and HLA-DQ8 (DQA1*03/DQB1*03:02) (11). Out of these HLA variants, HLA-DQ2.5 is the most strongly associated with CeD and is detected in around 90% of patients with CeD compared with about 20% in healthy controls (11). Deamidation of gluten peptides is mediated by the enzyme transglutaminase 2 (TG2) and results in conversion of Gln to negatively charged Glu at specific sites in polypeptides transforming proteolytically stable, but immunologically inert, peptides into pathogenic T cell epitopes (12, 13). A range of gluten T cell epitopes has been characterized, but four immunodominant epitopes derived from ╬▒-gliadin and ╧Й-gliadin are particularly prominent in the context of HLA-DQ2.5, namely, DQ2.5-glia-╬▒1a (PFPQPELPY), DQ2.5-glia-╬▒2 (PQPELPYPQ), DQ2.5-glia-╧Й1 (PFPQPEQPF), and DQ2.5-glia-╧Й2 (PQPEQPFPW). T cell responses against these epitopes are found in most HLA-DQ2.5+ patients and are considered to orchestrate tissue destruction in the small intestine and autoantibody production (14тАУ16). The only currently available treatment for CeD is lifelong adherence to a gluten-free diet. Development of alternative treatments is sought after because of poor patient compliance with the dietary restrictions, unavailability of strictly gluten-free food, and the feared transformation of uncomplicated CeD to refractory CeD (17, 18). Although it is still unclear how the destructive T cells emerge in CeD, recent findings have suggested a role for environmental factors including possible T cell cross-reactivity with microbial antigens (19, 20) and highlighted a need for further investigations into the roles of different MHC-presented antigens in CeD.
Recently, we reported the generation and use of Abs selected on HLA-DQ2.5:DQ2.5-glia-╬▒1a. We unexpectedly found that gut plasma cells (PCs) are the most abundant gluten peptideтАУpresenting cells in the inflamed small intestine of patients with CeD and that they express both HLA class II and T cell costimulatory molecules (9). In the current study, we describe the generation and affinity maturation of human TCR-like Abs specific for HLA-DQ2.5 in complex with DQ2.5-glia-╬▒2. We have previously shown that this complex elicits a prototypic signature CeD T cell response, which appears to be shared across all patients (21, 22). Our strategy was based on phage display selection of a fully human na├пve Ab library combined with secondary complementarity determining region (CDR)тАУtargeted optimization and gave rise to highly specific binders with picomolar monomeric affinities toward HLA-DQ2.5:DQ2.5-glia-╬▒2. The cocrystal structure of a Fab fragment of the lead Ab 3.C11 in complex with HLA-DQ2.5-glia-╬▒2 revealed an about diagonal docking geometry with respect to the peptide and a footprint highly similar to prototypic CeD patientтАУderived TCRs with the same specificity. The 3.C11 Ab readily stained a population of PCs from inflamed CeD lamina propria biopsy material, whereas the scarce CD11c+ and CD14+ dendritic cells (DCs) and macrophages (Mfs) from the corresponding material stained less prominently. These results confirm and extend our previous observation on the central role of the B cell compartment in this tissue (9). We further observed potent and strictly HLA- and epitope-dependent inhibition of T cell activation, both using CeD-derived CD4+ T cells in vitro, as well as primary CD4+ T cells in vivo upon oral administration of gluten peptide in HLA-DQ2.5 humanized mice, suggesting that the lead Ab 3.C11 has potential to be used for CeD-specific immunotherapy.
RESULTS
Primary selection of Abs specific for HLA-DQ2.5 with bound DQ2.5-glia-╬▒2
To generate human Abs specific for HLA-DQ2.5 in complex with the CeD epitope DQ2.5-glia-╬▒2, a human na├пve single-chain fragment variable (scFv) phage display library (23) was panned against soluble, recombinant pMHC. The method used to select the primary lead was based on our previous selection of a HLA-DQ2.5:DQ2.5-glia-╬▒1aтАУspecific Ab (9). After three rounds of selection, we assessed antigen reactivity of the selection output in a polyclonal enzyme-linked immunosorbent assay (ELISA) and observed increased and preferential binding to the target (fig. S1A). We then reformatted the selection output from the phagemid to a vector for soluble scFv expression and screened 70 single clones for target binding by ELISA (Fig. 1A). A total of 49 single clones bound target preferentially, and sequence analysis identified 14 unique sequences (Fig. 1A and fig. S1B). Five of the clones were enriched in the selection output. The V gene segment usage of the single clones was dominated by IGHV1-69, which paired with a diverse set of IGKV/LV segments (fig. S1B and table S1). To characterize the selected scFvs and choose a lead clone, we expressed all unique clones in Escherichia coli and directly compared periplasmic fractions in ELISA for target binding (Fig. 1B). Next, we performed pilot surface plasmon resonance (SPR) measurements using purified scFv (fig. S1C). On the basis of reactivity profiles and apparent affinities, we chose a lead clone, termed 206. When reformatted to full-length human immunoglobulin G1 (hIgG1), 206 retained binding to HLA-DQ2.5:DQ2.5-glia-╬▒2, albeit weakly (Fig. 1C and fig. S1D). No binding to the highly similar epitope HLA-DQ2.5:DQ2.5-glia-╧Й2 was observed (Fig. 1D). To accurately determine the monomeric affinity, a Fab version of 206 was constructed and SPR analysis estimated the monomeric affinity to KD 240 ┬▒ 20 nM with a high off-rate (2.4 ├Ч 10тИТ1 sтИТ1) (Fig. 1E).
(A) Selection output after three rounds of selection was batch-cloned into a vector for soluble scFv expression, and single clones were expressed in the E. coli periplasm followed by screening for binding to HLA-DQ2.5:DQ2.5-glia-╬▒2 (target) and HLA-DQ2.5:CLIP2 (background) in ELISA. ScFv anti-2-phenyloxazol-5-one (phOx) binding to BSA-phOx and empty E. coli were included as positive and negative controls, respectively. The ratios of target/background binding are represented as dots for each clone. Positive clones were sequenced, and colors represent clones found repeatedly. Gray dots denote clones with unknown sequence. (B) Fourteen unique clones were expressed, and periplasmic fractions were analyzed for binding to related gluten-pMHC complexes in ELISA. E. coli XL1-Blue were included as a negative control (ctrl), and pMHC capture levels were controlled with the HLA-DQ2 Ab 2.12.E11. Error bars illustrate means ┬▒ SD of duplicates (n = 2). A450nm, absorbance at 405 nm. (C) The lead candidate 206 was reformatted to hIgG1, and binding to a larger peptide panel of pMHC complexes was tested in ELISA. Error bars illustrate means ┬▒ SD of duplicates (n = 2). (D) Structural alignment of DQ2.5-glia-╬▒2 [PDB code: 4OZH (27)] and the related DQ2.5-glia-╧Й2 epitope (model). Differing positions are underlined. (E) The monomeric affinity of Fab 206 for HLA-DQ2.5:DQ2.5-glia-╬▒2 was determined by SPR using single-cycle kinetics and fitting a 1:1 binding model to the data (n = 2, representative sensorgram shown). RU, resonance units.
CDR-targeted secondary library design and selection of high-affinity gluten pMHCтАУspecific Abs
Given the low affinity of the 206 clone that was inadequate, for example, in detecting peptide on pulsed model APC (Fig. 6A), we sought to improve the Ab-pMHC interaction by targeted engineering. We generated focused libraries with sequence randomization and increased lengths of either the CDR H1 or H3 loop based on structural models assuming a TCR-like binding mode and the anticipated importance of the CDR H3 loop in target binding. We have previously found high valence (HV) pIX phage display to be very efficient at isolating Ab clones with high affinity and stability (23). Thus, the libraries were prepared accordingly and selected using a stringent strategy aimed at parallel identification of specific binders with high affinity and high thermostability in two separate arms (Fig. 2A). After an initial low-stringency round (R1), the libraries were split into a competition branch (CDR H1 and H3 selected separately, libraries denoted H1C and H3C) and a thermostability branch (CDR H1 and H3 pooled for selection, library denoted тАЬH1/H3TтАЭ) for a highly stringent R2 with low target concentration followed by a nonstringent R3 loosely based on the hammer-hug selection protocol (24). In R2 of the competition branch, scFvs were displayed at low valence (LV). In the thermal branch, scFvs were displayed at HV and subjected to a heat challenge at temperatures that induced unfolding of the parent clone (fig. S2A). This was done before selection to aggregate and remove unstable library members (6, 25). In R3, we aimed at recovering and amplifying selected binders and therefore increased antigen concentration. For the thermal branch, we included a second heat challenge and displayed scFvs at LV to favor clones with high monomeric affinities. The stringent competition in R2 of the competition branch resulted in close to no selection output in the CDR H3C library, and it was therefore discontinued. To determine antigen reactivity of the output, we performed a polyclonal phage ELISA (fig. S2B). The libraries showed signs of enrichment of binders, and we continued to screen single clones from the R3 output as soluble scFv (Fig. 2B) and as scFv displayed on phage (Fig. 2C) by ELISA. Clones with preferential binding to the target were present in both selection branches, and sequence analysis revealed that 66 of 73 were unique DNA sequences, and 64 were unique at the amino acid level (fig. S2C). All positive clones came from the CDR H1 library only, and most of the unique DNA sequences (56 of 66) were derived from the library with a loop length increased by two residues with the remaining clones (10 of 66) stemming from the pool with a length increase of three residues.
(A) Schematic overview of the selection strategy. LV and HV displays were achieved by packaging with the helper phage M13K07 or DeltaPhage, respectively. After R1, the libraries were split into a competition branch and a thermostability branch. (B) Selection outputs after R3 were screened as soluble scFvs to assess binding to target pMHC and HLA-DQ2.5:CLIP2 (background) in ELISA, and ratios were calculated. Each dot represents one clone. Gray dots denote unknown sequences, black dots denote unique single sequences, and colors represent enriched amino acid sequences. The H1/H3T library containing pooled CDR H1 and H3 libraries was selected in the thermostability branch, whereas the H1C library (CDR H1 only) was selected in the competition branch. (C) R3 selection outputs were screened in phage format and are represented as in (B). (D) Purified, monomeric Fab fragments were tested for binding to a panel of HLA-DQ2.5:peptide complexes in ELISA. Error bars illustrate means ┬▒ SD of duplicates (n = 2). (E) Sequence alignment of the mother clone 206 and the high-affinity offspring containing CDR H1 mutations (bold red). IGHV gene segment usage and numbering according to the international ImMunoGeneTics information system (IMGT). The IGHV1-69 V gene of clone 206 is in germline configuration.
On the basis of target binding in screening and enrichment of sequence features, we chose six clones for large-scale Fab expression in human embryonic kidney 293E cells. These were analyzed regarding their peptide specificity in ELISA, and all were found to bind their HLA-DQ2.5:DQ2.5-glia-╬▒2 target specifically (Fig. 2D and fig. S2D). The Fabs cross-reacted neither with the homologous HLA-DQ2.5:DQ2.5-glia-╧Й2 complex, which differs at three peptide positions only (p5, p7, and p9), nor with HLA-DQ2.5 with any of DQ2.5-glia-╧Й1, DQ2.5-glia-╬▒1a, or class IIтАУassociated invariant chain peptide 2 (CLIP2). Thus, the chosen candidates all bound specifically to their target despite the relatively broad selection of CDR H1 phenotypes (Fig. 2E).
Biophysical characterization of affinity-matured pMHC-specific Abs
We then performed binding analysis by SPR and ranked the six Abs on the basis of their off-rates (Figs. 2E and 3A and table S2). Strongly reduced off-rates were observed for all clones tested with clone 3.C11 exhibiting the lowest off-rate. In concordance with the improved (lower) off-rates, 3.C11 had a strong improvement in affinity with a KD 88 ┬▒ 8 pM (Fig. 3B and table S2). This is a 2700-fold improvement compared with the parent clone. 3.C11 was then expressed as full-length hIgG1 and tested for specific binding in ELISA (Fig. 3C). In agreement with previous results, 3.C11 bound exclusively to HLA-DQ2.5:DQ2.5-glia-╬▒2. Thus, the high-affinity Abs maintained the high specificity of the mother clone.
(A) Fab fragments were ranked on the basis of off-rates in SPR with the lead clone 3.C11 highlighted in red and the parent clone using dotted lines (n = 2). (B) Representative sensorgrams of Fab 3.C11 (n = 3). (C) 3.C11 was reformatted to full-length hIgG1 and analyzed in ELISA against a panel of related soluble peptide:HLA-DQ2.5 complexes. Error bars illustrate means ┬▒ SD of duplicates (n = 3). (D) Melting temperatures of mother clone Fab 206 (squared) and the affinity-matured Fab fragments with the lead candidate 3.C11 highlighted in red. Error bars illustrate means ┬▒ SD of triplicates (n = 1).
We next assessed the thermostability of all Fab fragments by determining their melting temperatures using nanoscale differential scanning fluorimetry (nanoDSF) (Fig. 3D). The high-affinity variants unexpectedly had slightly lower melting temperatures than their mother clone. The lead clone 3.C11 had the highest thermostabilities out of the engineered clones. Because of its favorable biophysical properties, we selected 3.C11 as a lead candidate for further characterization.
Structural basis for epitope specificity of 3.C11
To understand the molecular basis for the specificity and affinity of 3.C11, we crystallized the complex of the 3.C11 Fab fragment and HLA-DQ2.5:DQ2.5-glia-╬▒2 and solved its structure to 2.4-├Е resolution (table S3). The structure revealed that the 3.C11 Fab bound HLA-DQ2.5:DQ2.5-glia-╬▒2 with a тАЬTCR-likeтАЭ docking mode (Fig. 4, A and B) (26), which was unexpectedly similar to that observed for the prototypic HLA-DQ2.5:DQ2.5-glia-╬▒2тАУspecific TRAV26-1/TRBV7-2+ TCRs D2, S16, and JR5.1 (27). To characterize the TCR-like features of 3.C11, we compared the structure of the 3.C11 ternary complex with that of the S16 TCR (Fig. 4, C and D). The 3.C11 Fab engaged the peptide with all CDR loops (Fig. 4A) with an angle of 57┬░ between the VL and VH center of mass positions and the axis of the peptide (Fig. 4B). Thus, the VL and VH domains of 3.C11 (Fig. 4B) occupied the place of the canonical TCR V╬▒ and V╬▓ domain binding sites, respectively (Fig. 4D). The overall positioning of 3.C11 on HLA-DQ2.5:DQ2.5-glia-╬▒2 was similar to the TCR in the S16 ternary complex, albeit with the 3.C11 VL and VH center of mass positions offset toward the HLA-DQ2.5 ╬▒ chain helix by 3.7 and 6.1 ├Е compared with the S16 V╬▒ and V╬▓ center of mass positions, respectively (Fig. 4, B and D). The layout of the CDR loops relative to the pMHC and the footprint of 3.C11 (Fig. 4B) also broadly resembled that of the S16 TCR (Fig. 4D), such that 3.C11 and S16 shared a large number of pMHC contact residues.
(A) Overview of the Fab 3.C11-HLA-DQ2.5:DQ2.5-glia-╬▒2 complex. HLA-DQ2.5 ╬▒ and ╬▓ chains are light green and light yellow, respectively; the peptide is dark gray and shown in stick representation; 3.C11 H and L chains are brown and gray, respectively. CDRs and framework are colored as indicated in the figure. (B) CDR loop arrangement (top) and footprint (bottom) of the Fab 3.C11-HLA-DQ2.5:DQ2.5-glia-╬▒2 complex. In the footprint, HLA-DQ2.5:HLA-DQ2.5-glia-╬▒2 is shown in surface representation with the peptide in dark gray. Atoms contacted by the Fab are colored according to the nearest CDR loop or framework residue. Black dots and line represent the center of mass positions of the Ab variable domains and illustrate the approximate docking angle. (C) Overview and (D) CDR loop arrangement and footprint of the TCR S16-HLA-DQ2.5:DQ2.5-glia-╬▒2 complex illustrated for comparison (27).
The total buried surface area (BSA) of 3.C11 and HLA-DQ2.5:DQ2.5-glia-╬▒2 was 2385 ├Е2, which is moderately larger than that of the prototypic TCRs S16 (2130 ├Е2), D2 (2145 ├Е2), and JR5.1 (1765 ├Е2). The increase was mainly attributable to additional pMHC contacts made by the 3.C11 CDR H1 and CDR H2. With a BSA contribution of 76.2%, the 3.C11 H chain dominated the pMHC interface with relatively even contributions made by CDR H3 (27.5%), CDR H2 (22.7%), and CDR H1 (22.4%). The smaller interface of the L chain (23.8%) featured CDR L1 (11.2%) as the largest contributor (Fig. 4, B and D). The dominant role of the H chain in the 3.C11 pMHC interface broadly reflected the relative contributions of the ╬▓ chains observed for TCRs S16 (59.7%; Fig. 4D), D2 (62.8%), and JR5.1 (73.6%). Accordingly, the moderately larger BSA of 3.C11 compared with the TCRs contributes to the vastly higher affinity of 3.C11. The 3.C11 Fab-pMHC interface comprises 5 H-bonds and three salt bridges, which is numerically comparable to the 12 H-bonds present in the pMHC interface of S16 TCR. To further investigate the basis for the vastly higher affinity of 3.C11 compared with the S16 TCR, we calculated the geometric surface complementarity of the interfaces using shape complementarity (SC) (28), which resulted in nearly identical SC values for 3.C11 (SC = 0.691) and S16 (SC = 0.684). To investigate the charge complementarity of the complexes, we calculated the electrostatic surface potentials of the 3.C11 (fig. S3A) and S16 (fig. S3B) complexes using the Adaptive Poisson-Boltzmann Solver (29). The results revealed an overall similar polarity of the surface potential on 3.C11 and S16. 3.C11 showed an additional positively charged patch in the area of CDR H1, which was opposed by an overall weakly charged area on the pMHC (fig. S3A) and is therefore unlikely to drive the vast difference in affinities.
In accordance with the positional shift toward the HLA-DQ2.5 ╬▒ chain, compared with the S16 TCR, 3.C11 formed an extensive interface with the HLA-DQ2.5 ╬▒ chain (48.3%) and comparatively smaller interfaces with the HLA-DQ2.5 ╬▓ chain (35.3%) and the peptide (16.5%) (Fig. 4, B and D). The CDR H2 and CDR H3 loops made the largest contribution to the interface with the HLA-DQ2.5 ╬▒ chain. Here, the CDR H3 loop contacted the inner side of the peptide binding cleft, whereas the CDR H2 cupped the top and outer side of the HLA-DQ2.5 ╬▒ chain helix (Fig. 5A). The hydrophobic core of the CDR H2тАУHLA-DQ2.5 ╬▒ chain interface was composed of germline residues unique to IGHV1-69 (Ile57, Ile59, and Phe62), which suggests that the preferential IGHV1-69 selection as observed is in part driven by these interactions. Moreover, key elements of the 3.C11 Fab-pMHC interface were unexpectedly similar to features observed for the S16 TCRтАУHLA-DQ2.5:DQ2.5-glia-╬▒2 interface. Namely, the peptide residues directly contacted by 3.C11 (p2, p5, p6, p7, p8, p9, and p11) (Fig. 5B) included all peptide residues contacted by S16 (p2, p5, p6, p7, and p8) (Fig. 5C), and the interaction pattern of 3.C11 with the peptide (CDR L1 тЮФ p2; CDR H3 тЮФ p5, p6, p7; CDR H1 тЮФ p8, p9, p11) closely resembled that observed in the S16 complex (CDR ╬▒1 and CDR ╬▒3 тЮФ p2; CDR ╬▓3 тЮФ p5, p6, p7; CDR ╬▓1 тЮФ p8). Both the CDR H3 loop of 3.C11 (Fig. 4B) and the CDR ╬▓3 loop of S16 (Fig. 4D) formed brace-like structure across the central portion of the peptide and bridge across the peptide binding cleft (Fig. 5A).