Abstract

The activation of T cell immune responses, which relies on peptide antigens transported by TAP and bound to major histocompatibility complex (MHC) molecules, is recognized by T cell receptors (TCR). The quantitative prediction of MHC-epitope binding affinity can facilitate epitope screening and reduce cost and experimental efforts greatly. In this study, a comprehensive quantitative prediction method of binding affinity was established using quantitative structureactivity relationship (QSAR) modeling derived from amino acid physicochemical information. Firstly, the epitope was characterized by a set of amino acid physicochemical parameters. Secondly, the structural variables were optimized by the stepwise regression (STR). Finally, the robust quantitative models with were built by multiple linear regressions (MLR) for 31 MHC Class I subtypes. The normalized regression coefficients (NRCs) of QSAR model could demonstrate the mechanism of interaction of MHC, epitope, and TCR very well. The contribution of amino acid at each position of epitope, which was calculated by NRC, could determine which one was favorable for binding affinity or not. Therefore, the quantitative models established by STR-MLR could be used to guide virtual combinational design and high throughout screening of CTL epitope. Besides, they have many advantages, such as definite physiochemical indication, easier calculation and explanation, and good performances.

Keywords: Amino acid structural information, anchor position, epitope, Major Histocompatibility Complex (MHC), QSAR modeling, quantitative prediction.