Intramolecular Domain Movements of Free and Bound pMHC and TCR Proteins: A Molecular Dynamics Simulation Study

Identification of an HLA-A*11:01-restricted neoepitope of mutant PIK3CA and its specific T cell receptors for cancer immunotherapy targeting hotspot driver mutations

Hotspot driver mutations presented by human leukocyte antigens might be recognized by anti-tumor T cells. Based on their advantages of tumor-specificity and immunogenicity, neoantigens derived from hotspot mutations, such as PIK3CAH1047L, may serve as emerging targets for cancer immunotherapies. NetMHCpan V4.1 was utilized for predicting neoepitopes of PIK3CA hotspot mutation. Using in vitro stimulation, antigen-specific T cells targeting the HLA-A*11:01-restricted PIK3CA mutation were isolated from healthy donor-derived peripheral blood mononuclear cells. T cell receptors (TCRs) were cloned using single-cell PCR and sequencing. Their functionality was assessed through T cell activation markers, cytokine production and cytotoxic response to cancer cell lines pulsed with peptides or transduced genes of mutant PIK3CA. Immunogenic mutant antigens from PIK3CA and their corresponding CD8+ T cells were identified. These PIK3CA mutation-specific CD8+ T cells were subsequently enriched, and their TCRs were isolated. The TCR clones exhibited mutation-specific and HLA-restricted reactivity, demonstrating varying degrees of functional avidity. Identified TCR genes were transferred into CD8+ Jurkat cells and primary T cells deficient of endogenous TCRs. TCR-expressing cells demonstrated specific recognition and reactivity against the PIK3CAH1047L peptide presented by HLA-A*11:01-expressing K562 cells. Furthermore, mutation-specific TCR-T cells demonstrated an elevation in cytokine production and profound cytotoxic effects against HLA-A*11:01+ malignant cell lines harboring PIK3CAH1047L. Our data demonstrate the immunogenicity of an HLA-A*11:01-restricted PIK3CA hotspot mutation and its targeting therapeutic potential, together with promising candidates of TCR-T cell therapy.

Unraveling the Allosteric Communication Mechanisms in T-Cell Receptor-Peptide-Loaded Major Histocompatibility Complex Dynamics Using Molecular Dynamics Simulations: An Approach Based on Dynamic Cross Correlation Maps and Residue Interaction Energy Calculations

Antigen presentation by major histocompatibility complex (MHC) proteins to T-cell receptors (TCRs) plays a crucial role in triggering the adaptive immune response. Most of our knowledge on TCR-peptide-loaded major histocompatibility complex (pMHC) interaction stemmed from experiments yielding static structures, yet the dynamic aspects of this molecular interaction are equally important to understand the underlying molecular mechanisms and to develop treatment strategies against diseases such as cancer and autoimmune diseases. To this end, computational biophysics studies including all-atom molecular dynamics simulations have provided useful insights; however, we still lack a basic understanding of an overall allosteric mechanism that results in conformational changes in the TCR and subsequent T-cell activation. Previous hydrogen-deuterium exchange and nuclear magnetic resonance studies provided clues regarding these molecular mechanisms, including global rigidification and allosteric effects on the constant domain of TCRs away from the pMHC interaction site. Here, we show that molecular dynamics simulations can be used to identify how this overall rigidification may be related to the allosteric communication within TCRs upon pMHC interaction via essential dynamics and nonbonded residue-residue interaction energy analyses. The residues taking part in the rigidification effect are highlighted with an intricate analysis on residue interaction changes, which lead to a detailed outline of the complex formation event. Our results indicate that residues of the Cβ domain of TCRs show significant differences in their nonbonded interactions upon complex formation. Moreover, the dynamic cross correlations between these residues are also increased, in line with their nonbonded interaction energy changes. Altogether, our approach may be valuable for elucidating intramolecular allosteric changes in the TCR structure upon pMHC interaction in molecular dynamics simulations.

Molecular dynamics simulations provide molecular insights into the role of HLA-B51 in Behçet's disease pathogenesis

Behçet's disease is an inflammatory disorder of unknown etiology. Genetic tendency has an important role in its pathogenesis, and HLA-B51, a class I MHC antigen, has been recognized as the strongest susceptibility factor for Behçet's disease. Despite the confirmation of the association of HLA-B51 with Behçet's disease in different populations, its pathogenic mechanisms remain elusive. HLA-B51 differs in only two amino acids from HLA-B52, other split antigen of HLA-B5, which is not associated with Behçet's disease. These two amino acids are located in the B pocket of the antigen-binding groove, which occupies the second amino acids of the bound peptides. To understand the nature of the HLA-peptide interactions, differences in structure and dynamics of two HLA alleles were investigated by molecular dynamics simulations using YAYDGKDYI, LPRSTVINI, and IPYQDLPHL peptides. For HLA-B51, all bound peptides fluctuated to larger extent than HLA-B52. Free energy profiles of unbinding process for YAYDGKDYI by steered molecular dynamics simulations showed that unbinding from HLA-B52 results in greater free energy differences than HLA-B51. These results suggest the possibility of an instability of HLA-B51 associated with the repertoire of peptides, and this finding may provide significant insight to its pathogenic role in Behçet's disease.