Published ahead of print on January 11, 2007, doi:10.1165/rcmb.2006-0313CB
© 2007 American Thoracic Society DOI: 10.1165/rcmb.2006-0313CB
Prediction of Cellular Immune Responses against CFTR in Patients with Cystic Fibrosis after Gene TherapyGene Therapy Program, Department of Pathology and Laboratory Medicine, Division of Transfusion Medicine, University of Pennsylvania, Philadelphia, Pennsylvania Correspondence and requests for reprints should be addressed to James M. Wilson, M.D., Ph.D., 125 S. 31st Street, TRL, Suite 2000, Philadelphia, PA 19104-3403. E-mail: wilsonjm{at}mail.med.upenn.edu Abstract
Different classes of mutations (class IVI) of the cystic fibrosis (CF) transmembrane conductance regulator (CFTR) gene are responsible for lung/pancreatic disease. The most common mutation,
Key Words: MHC ligand CFTR gene therapy
Cystic fibrosis (CF) is a debilitating human disease caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene, which encodes a protein that functions as a cAMP-regulated chloride channel. Defects of the CFTR protein result in abnormal chloride transport across the apical membranes of epithelial cells in the airways, pancreas, intestine, and vas deferens, leading to progressive lung disease, pancreatic dysfunction, elevated sweat electrolyte levels, and male infertility, respectively. More than 1,000 mutations have been identified in the CFTR gene. With reference to chloride transport dysfunction, CFTR mutations can be grouped into five classes that reflect the associated biosynthetic or functional alterations in the CFTR protein: (I) CFTR not synthesized, (II) defective processing, (III) defective regulation, (IV) defective conductance, and (V) partially defective production or processing. While for classes I and II a major fraction of the CFTR protein does not reach the apical epithelial cell surface, it is present on the cellular surface in classes III, IV, or V with some residual function.
Correlation between genotype and phenotype has been shown for pancreatic status (1) and also, more recently, for airway disease (2). Patients with genotypes that include class I or II mutations on both chromosomes have more rapid deterioration in lung function and lower survival rates compared with the other genotypes with mutations that belong to class IV or V. The most common mutation among whites,
T cell development within the thymus is heavily influenced by interaction of the
Recognition of an antigenic peptide bound to an MHC protein (peptide-MHC) by the
The type and magnitude of the immune response would depend primarily upon the underlying mutation in CFTR. Lack of expression of full length or a fragment of the endogenous protein would result in diminished central tolerance, thereby leading to the generation of an immune response toward epitopes in the missing fragments by nondeleted T cells specific for the full-length protein. A greater response toward CFTR would therefore be observed when class I mutations are involved, wherein the expression of endogenous full-length CFTR is absent. For class II mutations, the Those self-reactive T cells that were spared from thymic deletion enter into the periphery and are maintained in a quiescent state. Because the affinity threshold of TCR-peptide-MHC interaction that signals thymic deletion is lower than that for activation in the periphery, it is likely that some T cells with low avidity for self-antigens will not be activated in the periphery and, instead, remain "ignorant" of their cognate antigen. However, ignorance of the antigen cannot be relied upon to maintain peripheral tolerance because, given the proper stimulatory milieu, such antigens may no longer be ignored and could potentially initiate autoimmune responses, as observed after viral infection (9, 10). In addition to passive means of tolerance, peripheral tolerance has been shown to be maintained actively by the tolerance mechanisms exerted by regulatory T cells (11). Thus, another potential mechanism by which immune responses may be generated toward therapeutic CFTR would be due to perturbation of mechanisms that control peripheral tolerance. To induce a cellular immune response, antigens have to be processed in intracellular compartments, transported, and presented by HLA molecules before recognition by specific T cells. T cell immune responses are driven by antigenic epitopes, and hence their identification is important for understanding disease pathogenesis and etiology. There are two types of T cell epitopes, CD8 and CD4, which are recognized in the context of the MHC-I and MHC-II molecules, respectively, by the correspondent T cell types. Appropriate processing of antigen peptides must occur before their binding to the appropriate MHC molecules. Detailed understanding of antigen processing and MHC-peptide recognition and binding has led to the development of several prediction algorithms. These algorithms primarily predict MHC binding, and more recently also include proteasomal cleavage recognition, crucial for the identification of CD8T cell epitopes. Some of the databases list published HLA ligands and T cell epitopes, while others offer prediction of T cell epitopes, for both Class I and Class II alleles.
In this study, we attempt to predict the likelihood of a host specific immune response mounted against the therapeutic CFTR protein in patients with the most frequent CFTR mutation, MATERIALS AND METHODS
Epitope prediction was performed using free web-based MHC-binding prediction algorithms. The algorithms/programs that were used in this study are described in Table 1. Amino acid sequence (488528) of human CFTR (fragment containing
BIMAS ranks potential 8-mer, 9-mer, or 10-mer peptides based on a predicted half-time of dissociation to HLA class I molecules. The analysis is based on coefficient tables deduced as previously described (12).
SYFPEITHI is an updated database with previous publications on T cell epitopes and MHC ligands (13) (access via: www.syfpeithi.de) The prediction is based on published motifs (pool sequencing, natural ligands) and takes into consideration the amino acids in the anchor and auxiliary anchor positions, as well as other frequent amino acids (13). Only those MHC class I alleles for which a large body of data is available are included in the "epitope prediction" section of SYFPEITHI. A reliability of at least 80% in retrieving the most apt epitope can be expected. Thus the naturally presented epitope should be among the top-scoring 2% of all peptides predicted in 80% of all predictions. For epitope predictions using MHC class II motifs, high reliabilities usually cannot be achieved due to the more variable pocket binding behavior. A reliability of only ProPred-I is an on-line service for identifying the MHC Class-I binding regions in antigens. This is a matrix-based method that allows the prediction of MHC-binding sites in an antigenic sequence for 47 MHC class-I alleles. The matrices used in ProPred-I have been obtained from both the BIMAS server and literature (14). Propred web is an interface allow users to predict MHC Class II binding regions in antigen sequence. The server employs amino acid/position coefficient tables deduced from literature (15), in a linear prediction model (i.e., quantitative matrix-based prediction method). Rankpep predicts peptide binders to MHC-I and MHC-II molecules from protein sequence/s or sequence alignments using position-specific scoring matrices (PSSMs). In addition, it predicts those MHC-I ligands whose C-terminal end is likely to be the consequence of proteasomal cleavage (16). RESULTS
The MHC prediction algorithms yielded several potential hits, restricted to both MHC-I and MHC-II, within the CFTR fragment at positions 488528. Overall, 20 potential MHC-I and 8 potential MHC-II ligands were obtained, as shown in Tables 2 and 3, respectively. These MHC ligands are restricted by more or less frequent MHC alleles, as shown by phenotype frequency in the general U.S. white population in both Tables 2 and 3. Approximately 30% of the MHC-I ligands were hits obtained with a high score (SYFPEITHI score
While ligands predicted by BIMAS are selected on the basis of their binding affinity with the respective MHC-I allele, those predicted by SYFPEITHI are selected on motif-based matrices developed from naturally occurring ligands or T cell epitopes. ProPred-I uses quantitative matrices with the option of the proteasome cleavage filter to identify potential T cell epitopes. We used more than one algorithm to predict potential MHC ligands to take advantage of the features associated with each of these programs. Some of the peptides scored positive by more than one program. For example, peptide KENIIFGVSY, which is restricted by HLA-B*44, was identified by all four algorithms used. Rankpep employs position specific scoring matrices to predict binders and also predicts those MHC-I ligands whose C-terminal end is likely to be the result of proteasomal cleavage, a classic feature of CD8T cell epitopes. Thus, peptide(s) T/IKENIIFGV restricted for HLA-A*0201 and A*0204 represent potential CD8T cell epitopes, by virtue of their C-terminal ends.
Interestingly, the potential CD8T cell epitope "TIKENIIFGV" is restricted by HLA-A*0201, one of the most frequent class I alleles in whites (Phenotype frequency = 3850%). Since the
Prediction of MHC-II ligands is challenging because peptides binding to a single MHC-II molecule are extremely variable in length and share very limited sequence similarity (1719). The RANKPEP algorithm, therefore, uses the motif discovery program MEME to create MHC-II profiles. On average, the sensitivity of these MHC-IIspecific profiles is such that In silico T cell epitope mapping using bioinformatics, when combined with other ex silico means of evaluating MHC-peptide and T cell interaction such as tetramers and HLA transgenic mice, has proved very useful in the field of vaccinology (21). However, there are instances in which theoretical predictions cannot be supported with experimental evidence, which may be explained by sensitivity of the experimental procedures, and also low precursor frequencies of antigen-specific cells. Our results suggest that there exists a number of potential CD8- and CD4-specific T cell epitopes within the CFTR fragment encompassing the "F" residue at position 508. The possibility of the CFTR gene therapy recipient mounting an immune response to the therapeutic CFTR would depend upon activation of T cells specific for epitopes within normal CFTR containing F508 residue (1) following an escape from central tolerance due to the underlying CFTR mutation, and/or (2) due to disturbance of peripheral tolerance toward CFTR.
Antigens are processed and presented differentially depending on their entry into the cell. Endogenously produces antigens are processed via the proteasomal pathway followed by transporter associated with antigen processing (TAP)-mediated transport of the peptides to the endoplasmic reticulum, where they are loaded onto class I MHC molecules and then presented on the cell surface in context of MHC-I. Exogenously derived antigens are processed via the lysosomal endocytic pathway, after which the peptides are loaded onto MHC-II for presentation on the cell surface. Occasionally, antigens may enter the cell through an exogenous pathway and get shunted toward proteasomal cleavage followed by peptide loading on MHC-I, leading to presentation in the context of MHC-I. This alternative pathway is termed "cross presentation." In patients who are homozygous for Footnotes Grants from the NIH (PO1-NL051746, P30-DK047757) (J.M.W.) and the CF Foundation supported this work. Originally Published in Press as DOI: 10.1165/rcmb.2006-0313CB on January 11, 2006 Conflict of Interest Statement: J.F. does not have a financial relationship with a commercial entity that has an interest in the subject of this manuscript. M.L. is an inventor on a patent that is licensed to Marker Gene Technologies Inc. J.M.W. received a grant from GlaxoSmithKline (GSK) and is an inventor on a patent owned by Penn and licensed to GSK, Targeted Genetics (TGEN), and Lentigen, and also receives fees from TGEN as a result of the license. Received in original form August 23, 2006 Accepted in final form November 17, 2006
References
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||