Donor and recipient human leukocyte antigen (HLA) incompatibility is the main immunological barrier to successful organ transplantation. HLA matching has a beneficial effect on solid organ transplant outcomes, including on graft function and on graft and patient survival, although the magnitude of this effect varies in different organs . Due to the extensive polymorphism of the HLA system and the limited size of the donor organ pool, HLA matching is difficult to achieve and most recipients receive HLA mismatched grafts. Moreover, the benefits of minimizing the number of donor–recipient HLA incompatibilities must be balanced against the, often competing, interests of ensuring equity of access to transplantation and of reducing the detrimental effects of prolonging the length of time on the transplant waiting list. Consequently, even though HLA matching is incorporated into many deceased-donor organ allocation schemes, the significance placed on improving HLA compatibility has progressively diminished in favor of increasing reliance on immunosuppression therapies to improve graft outcomes. It is evident, however, that HLA incompatible allografts necessitate the use of heavier immunosuppression, a major cause of recipient morbidity and mortality, and increase the risk of sensitization which can severely limit the opportunity for repeat organ transplantation.
The principles of HLA matching for organ allocation and of evaluating the immunological risk associated with a particular donor–recipient HLA combination have remained largely unchanged for many decades. Histocompatibility assessment, currently focuses almost entirely on HLA-A, HLA-B and HLA-DR loci, is based on counting HLA antigenic differences at the serological level and is predicated on the assumption that all mismatches within an HLA locus are of equal significance to transplant outcomes. However, it is now increasingly being recognized that the conventional approach to assessment of HLA compatibility is inadequate. Central to this notion is the realization that the capacity of a donor HLA to induce immune responses is not an intrinsic property of the mismatched alloantigen but that it is critically dependent on the HLA phenotype of the recipient. Recent advances in molecular sequence technology, increasing availability of high quality, crystallographically resolved, HLA structures and computational approaches to studying B-cell and T-cell epitopes offer exciting opportunities for a more rational approach to determining HLA compatibility and to defining HLA-related immunological risk. In this article, we review existing strategies for assessing HLA immunogenicity, including amino acid sequence, structural and physicochemical approaches, focusing on both B-cell and T-cell pathways of alloreactivity. We also discuss current challenges and opportunities in the field and consider future research directions that might transform clinical practice to maximize the benefits of transplantation.
ASSESSMENT OF HUMAN LEUKOCYTE ANTIGEN IMMUNOGENICITY AT THE AMINO ACID SEQUENCE LEVEL
The immunogenicity of nonself-antigenic determinants can be defined by their capacity to induce a host immune response. In the context of transplantation, HLA immunogenicity might be considered at the cellular (T-cell alloreactivity) and/or the humoral (B-cell alloreactivity leading to the development of HLA-specific antibodies) level, although it is evident that these immune responses are interdependent and intrinsically linked. The complexity of T-cell allorecognition and the absence of a readily available and easily interpretable assay to detect allospecific T cells in transplant patients have, so far, hampered the successful development of theoretical algorithms to predict the potential of donor HLA to induce T-cell alloreactivity . In contrast, the detection of alloantibody responses against HLA expressed on donor tissue has driven the field of HLA histocompatibility. The differential capacity of donor HLA to generate alloantibody responses and the favorable graft outcomes in the presence of ‘permissible’ HLA mismatches have long been recognized [3,4], driving attempts to create algorithms that define the potential of donor HLA to induce humoral alloresponses. HLA immunogenicity in the context of humoral alloimmunity is the focus of the present review.
The relationship between amino acid polymorphisms on donor HLA with alloantibody responses and kidney transplant outcomes was recognized over 2 decades ago and these observations have formed the basis of subsequent HLA immunogenicity algorithms [5,6]. Most significant in this regard has been the description of HLAMatchmaker by Duquesnoy et al. a computer algorithm that aims to determine the capacity of donor HLA to induce humoral alloimmunity by evaluating differences in the number and location of amino acid polymorphisms at continuous (triplets) and discontinuous (eplets) positions between donor and recipient HLA molecules [7,8]. The fundamental premise of HLAMatchmaker is that eplets, small patches of mismatched amino acids on or near the molecular surface of HLA, are potential immunogenic epitopes that define the specificity of HLA-specific antibodies. To verify this hypothesis, HLAMatchmaker has been used to analyze alloantibody binding profiles in patient sera with a view to describing HLA-specific antibody reactivity based on reactive eplets rather than HLA specificities (reviewed in [9,10]); such eplets, termed ‘antibody-verified’ eplets, are recorded in a publically accessible registry (http://www.epregistry.com.br/). It is evident, however, that this approach is not straight-forward [confounded by the presence of multiple alloantibodies in a patient's serum, variations in single-antigen bead assay analysis – including in the use of mean fluorescence intensity (MFI) cutoff values etc.] and often leads to complex theoretical interpretation, such as the hypothetical requirement for combinations of multiple eplets or of ‘self’ eplets, to explain alloantibody reactivity . It is not, therefore, clear whether B-cell epitopes are adequately defined by HLAMatchmaker eplets and the feasibility of an epitope (eplet)-based approach to donor–recipient matching, and to organ sharing, as has recently been suggested, remains doubtful .
The principal hypothesis underpinning theoretical approaches to predicting the risk of development of HLA-specific antibody is that HLA allorecognition by recipient B-cells is more likely the more ‘different’ the donor HLA is compared with recipient HLA molecules. In this regard, HLA immunogenicity algorithms should aim to quantify a measure of ‘dissimilarity’ between donor and recipient HLA. The HLAMatchmaker algorithm has been used successfully in this context as multiple studies have shown an association between the total number of eplets present on HLA class I and II mismatches (termed ‘eplet load’) and the risk of development of HLA-specific antibodies [13▪,14–16] (recently reviewed in ). Many studies in this field are often confounded by their retrospective design and the lack of multivariate modeling. Wiebe et al., however, have recently performed a relatively large, single-center, prospective study that benefited from high-resolution HLA typing, information on immunosuppression levels and patient adherence and from serial posttransplant antibody monitoring. Their analysis demonstrated that donor HLA-DR and HLA-DQ eplet load, along with tacrolimus trough levels, were independent predictors of development of donor-specific antibody (DSA) [17▪▪]. The clinical utility of this molecular definition of donor–recipient HLA compatibility has been highlighted by studies demonstrating an association between donor HLA eplet load and the development of transplant glomerulopathy, acute kidney graft rejection and graft loss, as well as chronic lung allograft dysfunction and pediatric heart transplant loss [18–21]. It is important to emphasize that the majority of the aforementioned studies on risk of humoral alloimmunity investigated the overall immunogenicity score (eplet load) of mismatches within an HLA locus (the sum of eplets for one or two HLA mismatches within the locus) and the likelihood of a locus-specific alloantibody response (to any of the mismatched alloantigens). The immunological basis of this type of analysis might be inconsistent with the underlying hypothesis in so far as, in cases in which a high eplet load reflects the sum of a low eplet and a high eplet HLA mismatch within the same locus, the alloantibody response might be directed against the low eplet donor HLA. Given that a molecular definition of HLA compatibility requires discrimination between low vs. high immunogenicity single HLA mismatches, future studies should aim to address this potential confounder.
An alternative approach to predicting HLA immunogenicity by assessment of the ‘dissimilarity’ between donor and recipient HLA, based on information derived from amino acid sequence analysis, has also been described. The Cambridge HLA immunogenicity algorithm performs interlocus (for HLA class I) and intralocus (for HLA class II) comparisons between donor and recipient HLA to enumerate all mismatched amino acid sequence polymorphisms on donor HLA which are then scored according to their physicochemical properties (hydrophobicity and electrostatic charge). Analysis of alloantibody responses in highly sensitized patients awaiting kidney transplantation, using this approach, showed that the capacity of donor HLA to induce a humoral response depends not only on the number of its polymorphic residues but also on the physicochemical nature of their side chains [22,23]. In these studies, consideration of only surface accessible amino acid polymorphisms did not improve the power of the model to predict alloantibody responses. Assessment of the electrostatic charge of amino acid polymorphisms on donor HLA-A and HLA-B alloantigens has been shown to correlate independently (accounting for the number of amino acid mismatches) with risk of DSA development in both highly sensitized patients and against HLA class I mismatches expressed on a failed kidney allograft [13▪,23]. However, due to the correlation between the number of amino acid polymorphisms expressed on a donor HLA and their overall physicochemical disparity, assessment of donor–recipient HLA physicochemical differences at the sequence level has not been shown to provide an advantage, over simple enumeration of amino acid mismatches, in assessing the risk of DSA responses against HLA class II alloantigens [13▪,22].
The aforementioned studies demonstrate that assessment of donor–recipient HLA mismatch at the amino acid sequence level enables quantification of the degree of mismatch improving the precision of immunological risk assessment with DSA development as the immune response readout. The potential benefit of using a particular computational algorithm to assess HLA immunogenicity was recently examined in a cohort of 596 renal transplant patients that were prospectively followed for DSA development against HLA-DR and HLA-DQ mismatches, accounting for immunosuppression therapy and recipient nonadherence [24▪]. This study showed that all donor HLA immunogenicity scoring methods, eplet mismatch load or number of amino acid mismatches or sequence-based electrostatic mismatch score (EMS), were significant multivariate correlates of DSA development outperforming conventional HLA mismatch grade for this purpose. No advantage was demonstrated in using one approach over another and each method provided equivalent assessment of immunological risk associated with donor HLA class II mismatches. As might be expected, there was strong correlation between molecular mismatch scores (R2 = 0.85–0.96) highlighting the fact they all reflect differences in donor–recipient amino acid sequence (Fig. 1).
ASSESSMENT OF HUMAN LEUKOCYTE ANTIGEN IMMUNOGENICITY AT THE TERTIARY STRUCTURE LEVEL
Even though amino acid sequence-based methods have been successful in predicting the risk of humoral alloimmunity after transplantation and in quantifying donor–recipient HLA incompatibility, sequence information alone provides limited insight into the structural aspects of HLA allorecognition by recipient B-cell receptors. Given the conformational nature of the majority of B-cell epitopes, a fully structural approach might be required to capture and characterize their immunogenic potential. Epitope–paratope interactions are largely governed by short-range electrostatic interactions, such as van der Waals forces, hydrogen bonds and salt-bridges , which are important determinants of the affinity and specificity of antibody–antigen binding . It has also been suggested that the process of affinity maturation involves optimization of electrostatic interactions in the B-cell receptor–antigen binding site [26–29]. Previous studies have shown that HLA B-cell epitopes can be described accurately by their unique surface electrostatic potential profiles that help explain serological patterns of HLA-specific antibody binding [30–32]. Given that recipient B-cells do not recognize self HLA molecules, it might be hypothesized that donor HLA with disparate electrostatic potential profiles compared with recipient HLA molecules might be recognized more efficiently by recipient B-cell receptors leading to improved selection and survival of differentiated B-cells during the process of affinity maturation in the germinal center. Indeed, recent insights into the mechanisms that determine the fate decision of proliferating, antigen-activated B-cells at the pregerminal center stage suggested that B-cells with higher affinity to their antigen presented more HLA-peptide to and made longer lasting contact with cognate T follicular helper cells at the B-cell–T-cell border in secondary lymphoid organs, resulting in more T-cell help and differentiation into germinal center B-cells [33,34]. These observations have recently led to the development of a novel HLA immunogenicity algorithm that enables quantitative comparison of surface electrostatic potential properties between donor and recipient HLA at the tertiary structure level (EMS-three–dimensional; a description of the methodology, applied to the study of HLA B-cell epitopes, is provided in ). This approach was applied to examine humoral alloimmune responses in healthy females subjected to a standardized injection of donor lymphocytes from their male partner. Preliminary results, published in abstract form [35,36] showed a strong association between the EMS-three-dimensional of donor HLA and donor-specific alloantibody development (Mallon et al., article in submission). The algorithm was also used to analyze donor–recipient HLA compatibility in a large cohort of kidney transplant recipients, and a preliminary report suggested a strong association between EMS-three-dimensional and risk of graft failure using a multivariate Cox regression model of graft survival . This is a developing field of research and further work is required to determine the potential of this approach to provide a fully structural description of HLA immunogenicity and whether it might enable improved assessment of donor–recipient histocompatibility compared to conventional HLA matching and to HLA amino acid sequence-based algorithms.
T-CELL APPROACHES TO PREDICTING HUMAN LEUKOCYTE ANTIGEN IMMUNOGENICITY
The above-mentioned approaches to determining the immunogenicity of HLA have predominantly focused on the B-cell component of humoral alloimmunity. Proliferation and differentiation of antigen-specific naïve B-cells into memory B-cells and long-lived plasma cells requires T-cell help through linked recognition of antigenic peptides presented in the context of B-cell HLA class II molecules . The implication of this mechanism is that the capacity of recipient HLA class II molecules to present donor-HLA-derived peptides is a critical determinant of the risk of DSA development, and this was supported by observational studies that showed an association between the recipient HLA-DR phenotype and humoral alloresponses to donor HLA class I alloantigens [40,41]. This concept has recently been applied into a computational algorithm that enables enumeration of putative CD4+ T-cell epitopes derived from donor HLA (predicted indirectly recognizable HLA epitopes or PIRCHE-II) . PIRCHE implements the NetMHCIIpan program which predicts peptide binding affinity into the groove of an HLA class II molecule based on a neural network method trained on a set of quantitative peptide binding data [43,44]. For a given peptide length, the program identifies potential core nonamer (9-mer) peptides contained within it and outputs the binding affinity of the highest scoring nonamer. In the PIRCHE implementation, predicted nonameric binding cores must be polymorphic compared with the recipient's own HLA molecules. The algorithm enables processing of donor HLA classes I and II alloantigens to derive putative candidate peptide fragments but assessment of peptide binding affinity is currently limited to recipient HLA-DR molecules. Donor-derived candidate peptides with predicted binding of IC50 (half maximal inhibitory concentration) less than 1000 nmol/l to recipient HLA-DR are summed to produce the PIRCHE score. The number of PIRCHE peptides derived from HLA class I mismatches was shown to correlate in univariate analyses with development of DSA after graft nephrectomy and after pregnancy, although there was a significant overlap in PIRCHE scores between immunogenic and nonimmunogenic donor HLA [42,45]. In a more recent study that investigated a large cohort of kidney transplant patients (n = 2787), the overall PIRCHE score for the total number of donor–recipient HLA-A, HLA-B, HLA-C, HLA-DR and HLA-DQB mismatches (sum of PIRCHE scores derived from every HLA mismatch in a given donor–recipient combination) was an independent predictor of risk of development of DSA (defined as alloantibody against any of the mismatched donor HLA; area-under-the-curve: 0.641), after adjustment for HLAMatchmaker eplet mismatching and for conventional HLA mismatching (there was no data on posttransplant immunosuppression regimens and on patient adherence) [46▪]. However, the same independent effect was not confirmed in a similar analysis of a subgroup of these patients (n = 1247) with the most complete pretransplant and posttransplant antibody monitoring in which the total number of HLA eplets was the best predictor of risk of DSA responses. It was notable that after analysis of locus-specific alloantibodies, the PIRCHE score did not correlate with risk of development of DSA against individual HLA-A, HLA-B, HLA-C, HLA-DR and HLA-DQB alloantigens when two HLA mismatches were present in the locus and, therefore, further studies are needed to determine whether PIRCHE can be used to predict the immunogenicity of individual HLA mismatches. When the effect of this approach on kidney graft outcomes was analyzed in this study, and more recently in an independent patient cohort [47▪], PIRCHE mismatching correlated with graft survival, although it is not clear whether this effect is independent of conventional HLA matching and/or eplet mismatching (only univariate analysis was performed in the Lachmann et al.[46▪] study and a forward stepwise selection multivariable model – that might not be appropriate in cases of multicollinearity – was used in the Geneugelijk et al.[47▪] study). The challenge of confirming an independent effect of a novel HLA scoring system on graft survival, over and above that observed with conventional HLA matching, is better appreciated considering the significant correlation between HLA scoring variables (e.g. there was high correlation between PIRCHE scores and HLAMatchmaker scores – Spearman's Rho of 0.75 – and between PIRCHE scores and conventional HLA mismatches in the study by Lachmann et al.[46▪]). To illustrate the latter observation, we examined the relationship between the number of possible peptide epitopes (unique, polymorphic nonameric binding cores) that can be derived from donor HLA (irrespective of predicted binding affinity to recipient HLA-DR), the number of PIRCHE epitopes (with predicted binding to recipient HLA-DR of IC50 <1000 nmol/l) and the number of amino acid polymorphisms present in HLA-A, HLA-B, HLA-C, HLA-DR, HLA-DQ and HLA-DP mismatches within a local cohort of 182 donor–recipient pairs. As shown in Fig. 2, there was high correlation between HLA scoring systems (although there was variation of PIRCHE scores for the same number of amino acid polymorphisms) suggesting that, to an extent, the current PIRCHE approach reflects donor HLA polymorphism at the amino sequence level. Overall, PIRCHE is an interesting approach for assessing potential CD4+ T-cell epitopes derived from donor HLA and how this might impact on risk of humoral alloimmunity after transplantation. Further improvements of CD4+ T-cell alloreactivity prediction algorithms might arise from consideration of peptide presentation by recipient HLA-DQ and/or HLA-DP and from experimental validation of predicted epitopes that could inform the theoretical design of computational algorithms (e.g. by investigating the predicted peptide binding affinity cutoff values that best reflect experimental observations).
CONCLUDING REMARKS AND FUTURE DIRECTIONS
There is currently great need to develop and implement clinically applicable algorithms that enable better assessment of donor–recipient HLA incompatibility and of the associated immunological risk, both at a population and at the individual patient level, to maximize the benefits of transplantation. It is evident from the aforementioned studies that significant progress has been made, but there are a number of challenges to overcome to facilitate a sea change in the field. Significant in this regard is the heterogeneity of studies investigating HLA immunogenicity that limits interpretation and synthesis of available evidence. This heterogeneity pertains to a number of factors that introduce uncertainty, including the level of HLA typing resolution (serological, one-field, two-field), Luminex single antigen bead output interpretation (e.g. MFI cutoff values used for antibody detection), availability of prospective DSA monitoring, adjustment for significant confounders (e.g. immunosuppression regimen, patient adherence) and methodological model used for determining the association between HLA immunogenicity score and DSA response (e.g. immunogenicity score of individual HLA mismatches and DSA response to the same antigen vs. total immunogenicity score at an HLA locus level and alloantibody response to any of the locus mismatches vs. total immunogenicity score of all donor–recipient HLA mismatches and alloantibody response to any of these mismatches). International collaboration and setting up of multicenter consortia will be required to address these sources of bias and facilitate the design of large, appropriately powered, datasets to investigate the predictive power of HLA immunogenicity algorithms and how they might be applied in the clinical setting. Future studies should focus on investigating how HLA immunogenicity algorithms might be applied to improve organ allocation policies (provided that a clear benefit over conventional HLA matching at the antigen level is demonstrated) and model how any changes might impact equity of access to transplantation. In addition, prospective studies will be needed to evaluate the potential of HLA immunogenicity analysis to assess immunological risk and facilitate clinical decision making at individual patient level . There is also an impetus on continuing improvement of currently available algorithms and/or development of new approaches. This relates to both B-cell and T-cell HLA immunogenicity algorithms which might be further improved by better definition of the structure and physicochemical properties of HLA B-cell epitopes (e.g. through experimental resolution of alloantibody-HLA structures and computational approaches ), more accurate prediction of T-cell epitopes (including experimental epitope validation and cross-fertilization from relevant fields [50–52]), and from development of combined approaches that incorporate all aspects of alloreactivity.
Financial support and sponsorship
The current study was supported by the Cambridge NIHR Biomedical Research Centre and the NIHR Blood and Transplant Research Unit in Organ Donation and Transplantation at the University of Cambridge in collaboration with Newcastle University and in partnership with NHS Blood and Transplant (NHSBT). The views expressed are those of the authors and not necessarily those of the NHS, the NIHR, the Department of Health or NHSBT. V.K. was supported by an Evelyn Trust Grant and an NIHR Post-Doctoral Fellowship (PDF-2016-09-065).
Conflicts of interest
There are no conflicts of interest.
REFERENCES AND RECOMMENDED READING
Papers of particular interest, published within the annual period of review, have been highlighted as:
▪ of special interest
▪▪ of outstanding interest
1. Zachary AA, Leffell MS. HLA mismatching strategies for solid organ transplantation – a balancing act. Front Immunol 2016; 7:575.
2. Claas FH, Roelen DL, Mulder A, et al. Differential immunogenicity of HLA class I alloantigens for the humoral versus the cellular immune response: ‘towards tailor-made HLA mismatching’. Hum Immunol 2006; 67:424.
3. Maruya E, Takemoto S, Terasaki PI. HLA matching: identification of permissible HLA mismatches. Clin Transpl 1993; 511–520.
4. Doxiadis II, Smits JM, Schreuder GM, et al. Association between specific HLA combinations and probability of kidney allograft loss: the taboo concept. Lancet 1996; 348:850–853.
5. Barbetti AA, Park MS, Terasaki PI, et al. HLA class II epitope detection by serology. Clin Transpl 1990; 533–565.
6. Takemoto S, Terasaki PI. HLA epitopes and graft survival. Clin Transpl 1991; 363–383.
7. Duquesnoy RJ. A structurally based approach to determine HLA compatibility at the humoral immune level. Hum Immunol 2006; 67:847–862.
8. Duquesnoy RJ. HLAMatchmaker: a molecularly based algorithm for histocompatibility determination. I. Description of the algorithm. Hum Immunol 2002; 63:339–352.
9. Duquesnoy RJ. Reflections on HLA epitope-based matching for transplantation. Front Immunol 2016; 7:469.
10. Duquesnoy RJ, Marrari M. HLAMatchmaker-based definition of structural human leukocyte antigen epitopes detected by alloantibodies. Curr Opin Organ Transplant 2009; 14:403–409.
11. Marrari M, Mostecki J, Mulder A, et al. Human monoclonal antibody reactivity with human leukocyte antigen class I epitopes defined by pairs of mismatched eplets and self-eplets. Transplantation 2010; 90:1468–1472.
12. Duquesnoy RJ. Are we ready for epitope-based HLA matching in clinical organ transplantation? Transplantation 2017; 101:1755–1765.
13▪. Kosmoliaptsis V, Mallon DH, Chen Y, et al. Alloantibody responses after renal transplant failure can be better predicted by donor-recipient HLA amino acid sequence and physicochemical disparities than conventional HLA matching. Am J Transplant 2016; 16:2139–2147.
14. Kosmoliaptsis V, Bradley JA, Sharples LD, et al. Predicting the immunogenicity of human leukocyte antigen class I alloantigens using structural epitope analysis determined by HLAMatchmaker. Transplantation 2008; 85:1817–1825.
15. Wiebe C, Nickerson P. Acceptable mismatching at the class II epitope level: the Canadian experience. Curr Opin Organ Transplant 2014; 19:442–446.
16. Wiebe C, Pochinco D, Blydt-Hansen TD, et al. Class II HLA epitope matching – a strategy to minimize de novo donor-specific antibody development and improve outcomes. Am J Transplant 2013; 13:3114–3122.
17▪▪. Wiebe C, Rush DN, Nevins TE, et al. Class II eplet mismatch modulates tacrolimus trough levels required to prevent donor-specific antibody development. J Am Soc Nephrol 2017; 28:3353–3362.
18. Sapir-Pichhadze R, Tinckam K, Quach K, et al. HLA-DR and -DQ eplet mismatches and transplant glomerulopathy: a nested case–control study. Am J Transplant 2015; 15:137–148.
19. Wiebe C, Nevins TE, Robiner WN, et al. The synergistic effect of class II HLA epitope-mismatch and nonadherence on acute rejection and graft survival. Am J Transplant 2015; 15:2197–2202.
20. Sullivan PM, Warner P, Kemna MS, et al. HLA molecular epitope mismatching and long-term graft loss in pediatric heart transplant recipients. J Heart Lung Transplant 2015; 34:950–957.
21. Walton DC, Hiho SJ, Cantwell LS, et al. HLA matching at the eplet level protects against chronic lung allograft dysfunction. Am J Transplant 2016; 16:2695–2703.
22. Kosmoliaptsis V, Sharples LD, Chaudhry AN, et al. Predicting HLA class II alloantigen immunogenicity from the number and physiochemical properties of amino acid polymorphisms. Transplantation 2011; 91:183–190.
23. Kosmoliaptsis V, Chaudhry AN, Sharples LD, et al. Predicting HLA class I alloantigen immunogenicity from the number and physiochemical properties of amino acid polymorphisms. Transplantation 2009; 88:791–798.
24▪. Wiebe C, Kosmoliaptsis V, Pochinco D, et al. A comparison of HLA molecular mismatch methods to determine HLA immunogenicity. Transplantation 2018; doi: 10.1097/TP.0000000000002117. [Epub ahead of print].
25. Persson BA, Jonsson B, Lund M. Enhanced protein steering: cooperative electrostatic and van der Waals forces in antigen–antibody complexes. J Phys Chem B 2009; 113:10459–10464.
26. Sinha N, Smith-Gill SJ. Electrostatics in protein binding and function. Curr Protein Pept Sci 2002; 3:601–614.
27. Nakamura H. Roles of electrostatic interaction in proteins. Q Rev Biophys 1996; 29:1–90.
28. Chong LT, Duan Y, Wang L, et al. Molecular dynamics and free-energy calculations applied to affinity maturation in antibody 48G7. Proc Natl Acad Sci U S A 1999; 96:14330–14335.
29. Lippow SM, Wittrup KD, Tidor B. Computational design of antibody-affinity improvement beyond in vivo maturation. Nat Biotechnol 2007; 25:1171–1176.
30. Mallon DH, Bradley JA, Winn PJ, et al. Three-dimensional structural modelling and calculation of electrostatic potentials of HLA Bw4 and Bw6 epitopes to explain the molecular basis for alloantibody binding: toward predicting HLA antigenicity and immunogenicity. Transplantation 2015; 99:385–390.
31. Mallon DH, Bradley JA, Taylor CJ, Kosmoliaptsis V. Structural and electrostatic analysis of HLA B-cell epitopes: inference on immunogenicity and prediction of humoral alloresponses. Curr Opin Organ Transplant 2014; 19:420–427.
32. Kosmoliaptsis V, Dafforn TR, Chaudhry AN, et al. High-resolution, three-dimensional modeling of human leukocyte antigen class I structure and surface electrostatic potential reveals the molecular basis for alloantibody binding epitopes. Hum Immunol 2011; 72:1049–1059.
33. Schwickert TA, Victora GD, Fooksman DR, et al. A dynamic T cell-limited checkpoint regulates affinity-dependent B cell entry into the germinal center. J Exp Med 2011; 208:1243–1252.
34. Allen CD, Okada T, Tang HL, Cyster JG. Imaging of germinal center selection events during affinity maturation. Science 2007; 315:528–531.
35. Mallon DH, Kling C, Bradley JA, et al. Computational scoring system to predict HLA immunogenicity. Lancet 2016; 387:S68.
36. Mallon DH, Kling C, Robb M, et al. Validation of a computational scoring system for predicting HLA immunogenicity based on quantification of surface electrostatic potential differences between donor and recipient HLA molecules. Transpl Int 2017; 30:31.
37. Mallon DH, Kling C, Robb M, et al. Predicting humoral alloimmunity from differences in donor-recipient HLA surface electrostatic potential. bioRxiv 2018; doi: 294066.
38. Kosmoliaptsis V, Mallon DH, Fuggle S, et al. A novel computational HLA matching algorithm for improving donor-recipient histocompatibility and graft outcomes after kidney transplantation. Am J Transplant 2017; 17:420.
39. Ali JM, Bolton EM, Bradley JA, Pettigrew GJ. Allorecognition pathways in transplant rejection and tolerance. Transplantation 2013; 96:681–688.
40. Fuller TC, Fuller A. The humoral immune response against an HLA class I allodeterminant correlates with the HLA-DR phenotype of the responder. Transplantation 1999; 68:173–182.
41. Dankers MK, Roelen DL, Nagelkerke NJ, et al. The HLA-DR phenotype of the responder is predictive of humoral response against HLA class I antigens. Hum Immunol 2004; 65:13–19.
42. Geneugelijk K, Hönger G, van Deutekom HW, et al. Predicted indirectly recognizable HLA epitopes presented by HLA-DRB1 are related to HLA antibody formation during pregnancy. Am J Transplant 2015; 15:3112–3122.
43. Andreatta M, Karosiene E, Rasmussen M, et al. Accurate pan-specific prediction of peptide-MHC class II binding affinity with improved binding core identification. Immunogenetics 2015; 67:641–650.
44. Nielsen M, Lundegaard C, Blicher T, et al. Quantitative predictions of peptide binding to any HLA-DR molecule of known sequence: NetMHCIIpan. PLoS Comput Biol 2008; 4:e1000107.
45. Otten HG, Calis JJ, Kesmir C, et al. Predicted indirectly recognizable HLA epitopes presented by HLA-DR correlate with the de novo development of donor-specific HLA IgG antibodies after kidney transplantation. Hum Immunol 2013; 74:290–296.
46▪. Lachmann N, Niemann M, Reinke P, et al. Donor-recipient matching based on predicted indirectly recognizable HLA epitopes independently predicts the incidence of de novo donor-specific HLA antibodies following renal transplantation. Am J Transplant 2017; 17:3076–3086.
47▪. Geneugelijk K, Niemann M, Drylewicz J, et al. PIRCHE-II is related to graft failure after kidney transplantation. Front Immunol 2018; 9:321.
48. Wiebe C, Nickerson P. Strategic use of epitope matching to improve outcomes. Transplantation 2016; 100:2048–2052.
49. Potocnakova L, Bhide M, Pulzova LB. An introduction to B-cell epitope mapping and in silico epitope prediction. J Immunol Res 2016; 2016:6760830.
50. Sanchez-Trincado JL, Gomez-Perosanz M, Reche PA. Fundamentals and methods for T- and B-cell epitope prediction. J Immunol Res 2017; 2017:2680160.
51. Wang P, Sidney J, Kim Y, et al. Peptide binding predictions for HLA DR, DP and DQ molecules. BMC Bioinformatics 2010; 11:568.
52. Agudelo WA, Patarroyo ME. Quantum chemical analysis of MHC-peptide interactions for vaccine design. Mini Rev Med Chem 2010; 10:746–758.