Dynamically Driven Allostery in MHC Proteins: Peptide-Dependent Tuning of Class I MHC Global Flexibility

Enhanced T cell receptor specificity through framework engineering

Development of T cell receptors (TCRs) as immunotherapeutics is hindered by inherent TCR cross-reactivity. Engineering more specific TCRs has proven challenging, as unlike antibodies, improving TCR affinity does not usually improve specificity. Although various protein design approaches have been explored to surmount this, mutations in TCR binding interfaces risk broadening specificity or introducing new reactivities. Here we explored if TCR specificity could alternatively be tuned through framework mutations distant from the interface. Studying the 868 TCR specific for the HIV SL9 epitope presented by HLA-A2, we used deep mutational scanning to identify a framework mutation above the mobile CDR3β loop. This glycine to proline mutation had no discernable impact on binding affinity or functional avidity towards the SL9 epitope but weakened recognition of SL9 escape variants and led to fewer responses in a SL9-derived positional scanning library. In contrast, an interfacial mutation near the tip of CDR3α that also did not impact affinity or functional avidity towards SL9 weakened specificity. Simulations indicated that the specificity-enhancing mutation functions by reducing the range of loop motions, limiting the ability of the TCR to adjust to different ligands. Although our results are likely to be TCR dependent, using framework engineering to control TCR loop motions may be a viable strategy for improving the specificity of TCR-based immunotherapies.

Structural and physical features that distinguish tumor-controlling from inactive cancer neoepitopes

Neoepitopes arising from amino acid substitutions due to single nucleotide polymorphisms are targets of T cell immune responses to cancer and are of significant interest in the development of cancer vaccines. However, understanding the characteristics of rare protective neoepitopes that truly control tumor growth has been a challenge, due to their scarcity as well as the challenge of verifying true, neoepitope-dependent tumor control in humans. Taking advantage of recent work in mouse models that circumvented these challenges, here, we compared the structural and physical properties of neoepitopes that range from fully protective to immunologically inactive. As neoepitopes are derived from self-peptides that can induce immune tolerance, we studied not only how the various neoepitopes differ from each other but also from their wild-type counterparts. We identified multiple features associated with protection, including features that describe how neoepitopes differ from self as well as features associated with recognition by diverse T cell receptor repertoires. We demonstrate both the promise and limitations of neoepitope structural analysis and predictive modeling and illustrate important aspects that can be incorporated into neoepitope prediction pipelines.

Conformational plasticity of RAS Q61 family of neoepitopes results in distinct features for targeted recognition

The conformational landscapes of peptide/human leucocyte antigen (pHLA) protein complexes encompassing tumor neoantigens provide a rationale for target selection towards autologous T cell, vaccine, and antibody-based therapeutic modalities. Here, using complementary biophysical and computational methods, we characterize recurrent RAS55-64 Q61 neoepitopes presented by the common HLA-A*01:01 allotype. We integrate sparse NMR restraints with Rosetta docking to determine the solution structure of NRASQ61K/HLA-A*01:01, which enables modeling of other common RAS55-64 neoepitopes. Hydrogen/deuterium exchange mass spectrometry experiments alongside molecular dynamics simulations reveal differences in solvent accessibility and conformational plasticity across a panel of common Q61 neoepitopes that are relevant for recognition by immunoreceptors. Finally, we predict binding and provide structural models of NRASQ61K antigens spanning the entire HLA allelic landscape, together with in vitro validation for HLA-A*01:191, HLA-B*15:01, and HLA-C*08:02. Our work provides a basis to delineate the solution surface features and immunogenicity of clinically relevant neoepitope/HLA targets for cancer therapy.

Be2+ Causes Hypersensitivity but Mg2+ and Ca2+ Do Not─Favorable Metal Coordination Is the Key for Differential Allosteric Modulation and Binding Affinities

Although the ion selectivity of metalloproteins has been well established, selective metal antigen recognition by immunoproteins remains elusive. One such case is the recognition of the Be2+ ion against its heavier congeners, Mg2+ and Ca2+, by the human leukocyte antigen immunoprotein (HLA-DP2), leading to immunotoxicity. Integrating with our previous mechanistic study on Be2+ toxicity, herein, we have explored the basis of characteristic nontoxicity of Mg2+ and Ca2+ ions despite their in vivo abundance. The ion binding cleft of the HLA-DP2-peptide complex is composed of four acidic residues, p4D and p7E from the peptide and β26E and β69E from the protein. While the tetrahedral coordination site of the smaller Be2+ ion is located deep inside the cavity, hexa- to octa-coordination sites of Mg2+ and Ca2+ ions are located closer to the protein surface. The intrinsic high coordination number of Mg2+/Ca2+ ions induces allosteric modifications on the HLA-DP2_M2 surface, which are atypical for TCR recognition. Furthermore, the lower binding energy of larger Mg2+ and Ca2+ ions with the cavity residues can be correlated to the lower charge density and reduced covalent bonding nature as compared to those of the smaller Be2+ ion. In short, weak binding of Mg2+ and Ca2+ ions and the unfavorable allosteric surface modifications are probably the major determinants for the absence of Mg2+/Ca2+ ion-mediated hypersensitivity in humans.

Coarse-grained molecular dynamics-guided immunoinformatics to explain the binder and non-binder classification of Cytotoxic T-cell epitope for SARS-CoV-2 peptide-based vaccine discovery

Epitope-based peptide vaccine can elicit T-cell immunity against SARS-CoV-2 to clear the infection. However, finding the best epitope from the whole antigen is challenging. A peptide screening using immunoinformatics usually starts from MHC-binding peptide, immunogenicity, cross-reactivity with the human proteome, to toxicity analysis. This pipeline classified the peptides into three categories, i.e., strong-, weak-, and non-binder, without incorporating the structural aspect. For this reason, the molecular detail that discriminates the binders from non-binder is interesting to be investigated. In this study, five CTL epitopes against HLA-A*02:01 were identified from the coarse-grained molecular dynamics-guided immunoinformatics screening. The strong binder showed distinctive activities from the non-binder in terms of structural and energetic properties. Furthermore, the second residue from the nonameric peptide was most important in the interaction with HLA-A*02:01. By understanding the nature of MHC-peptide interaction, we hoped to improve the chance of finding the best epitope for a peptide vaccine candidate.

Elucidation of binding mechanism, affinity, and complex structure between mWT1 tumor-associated antigen peptide and HLA-A*24:02

We have applied our advanced computational and experimental methodologies to investigate the complex structure and binding mechanism of a modified Wilms' Tumor 1 (mWT1) protein epitope to the understudied Asian-dominant allele HLA-A*24:02 (HLA-A24) in aqueous solution. We have applied our developed multicanonical molecular dynamics (McMD)-based dynamic docking method to analyze the binding pathway and mechanism, which we verified by comparing the highest probability structures from simulation with our experimentally solved x-ray crystal structure. Subsequent path sampling MD simulations elucidated the atomic details of the binding process and indicated that first an encounter complex is formed between the N-terminal's positive charge of the 9-residue mWT1 fragment peptide and a cluster of negative residues on the surface of HLA-A24, with the major histocompatibility complex (MHC) molecule preferring a predominantly closed conformation. The peptide first binds to this closed MHC conformation, forming an encounter complex, after which the binding site opens due to increased entropy of the binding site, allowing the peptide to bind to form the native complex structure. Further sequence and structure analyses also suggest that although the peptide loading complex would help with stabilizing the MHC molecule, the binding depends in a large part on the intrinsic affinity between the MHC molecule and the antigen peptide. Finally, our computational tools and analyses can be of great benefit to study the binding mechanism of different MHC types to their antigens, where it could also be useful in the development of higher affinity variant peptides and for personalized medicine.

Crystal structures of MHC class I complexes reveal the elusive intermediate conformations explored during peptide editing

Studies have suggested that MHC class I (MHC I) molecules fluctuate rapidly between numerous conformational states and these motions support peptide sampling. To date, MHC I intermediates are largely uncharacterized experimentally and remain elusive. Here, we present x-ray crystal structures of HLA-B8 loaded with 20mer peptides that show pronounced distortions at the N-terminus of the groove. Long stretches of N-terminal amino acid residues are missing in the electron density maps creating an open-ended groove. Our structures also reveal highly unusual features in MHC I-peptide interaction at the N-terminus of the groove. Molecular dynamics simulations indicate that the complexes have varying degrees of conformational flexibility in a manner consistent with the structures. We suggest that our structures have captured the remarkable molecular dynamics of MHC I-peptide interaction. The visualization of peptide-dependent conformational motions in MHC I is a major step forward in our conceptual understanding of dynamics in high-affinity peptide selection.

Visualising tapasin- and TAPBPR-assisted editing of major histocompatibility complex class-I immunopeptidomes

Which peptides are selected for presentation by major histocompatibility complex class-I (MHC-I) molecules is a key determinant of successful immune responses. Peptide selection is co-ordinated by the tapasin and TAP Binding PRotein (TAPBPR) proteins, which ensure MHC-I molecules preferentially acquire high-affinity-binding peptides. New structural analyses have offered insight into how tapasin achieves this function within the peptide-loading complex (PLC) (comprising the Transporter associated with Antigen Presentation (TAP) peptide transporter, tapasin-ERp57, MHC-I and calreticulin), and how TAPBPR performs a peptide editing function independently of other molecules. The new structures reveal nuances in how tapasin and TAPBPR interact with MHC-I, and how calreticulin and ERp57 complement tapasin to exploit the plasticity of MHC-I molecules to achieve peptide editing.

Usefulness of Docking and Molecular Dynamics in Selecting Tumor Neoantigens to Design Personalized Cancer Vaccines: A Proof of Concept

Personalized cancer vaccines based on neoantigens are a new and promising treatment for cancer; however, there are still multiple unresolved challenges to using this type of immunotherapy. Among these, the effective identification of immunogenic neoantigens stands out, since the in silico tools used generate a significant portion of false positives. Inclusion of molecular simulation techniques can refine the results these tools produce. In this work, we explored docking and molecular dynamics to study the association between the stability of peptide-HLA complexes and their immunogenicity, using as a proof of concept two HLA-A2-restricted neoantigens that were already evaluated in vitro. The results obtained were in accordance with the in vitro immunogenicity, since the immunogenic neoantigen ASTN1 remained bound at both ends to the HLA-A2 molecule. Additionally, molecular dynamic simulation suggests that position 1 of the peptide has a more relevant role in stabilizing the N-terminus than previously proposed. Likewise, the mutations may have a "delocalized" effect on the peptide-HLA interaction, which means that the mutated amino acid influences the intensity of the interactions of distant amino acids of the peptide with the HLA. These findings allow us to propose the inclusion of molecular simulation techniques to improve the identification of neoantigens for cancer vaccines.

Partial peptide dissociation and binding groove plasticity in two major histocompatibility complex class I alleles - differences between alleles versus force field and sampling effects

Major histocompatibility complex class I (MHC I) reports a cell's health status by presenting antigenic peptides inside its binding groove. However, MHC I binding grooves can differ largely in their plasticity, from binding grooves that are conformationally stable by themselves to those that require a high-affinity peptide to be bound to attain conformational stability. These latter MHC I alleles are dependent on the C-terminus of the peptide that stabilizes the F-pocket region of their binding grooves. It has remained unclear to what extent a peptide-MHC I complex (pMHC I) can tolerate the (intermittent) partial dissociation of high-affinity peptides, especially of the peptide's N-terminus. Using bias exchange umbrella sampling (BEUS), a technique to achieve enhanced sampling in molecular dynamics (MD) simulations, we obtained the free-energy profiles of the N-terminal dissociation of a respective high-affinity peptide from HLA-B*35:01 and HLA-B*44:02, two alleles on opposite ends of the scale regarding binding groove plasticity. The potential of mean force (PMF) for HLA-B*35:01 was calculated for two different sets of starting structures and is compared with a PMF obtained previously with a different force field to disentangle allele differences from force field and sampling effects. For both alleles, the free-energy profiles indicate that the peptide N-terminus dissociates in a substantial fraction of the pMHC I, suggesting that their crystal structures with fully bound peptides only partially capture the dynamic conformational ensemble of pMHC I in solution, and thus in the cell.

In solution, high-affinity peptides are likely to dissociate partially from two alleles of major histocompatibility complex I. Despite very similar free-energy profiles, two molecular dynamics force fields predict different underlying mechanisms.

Phenol sensing in nature is modulated via a conformational switch governed by dynamic allostery

The NtrC family of proteins senses external stimuli and accordingly stimulates stress and virulence pathways via activation of associated σ⁵⁴-dependent RNA polymerases. However, the structural determinants that mediate this activation are not well understood. Here, we establish using computational, structural, biochemical, and biophysical studies that MopR, an NtrC protein, harbors a dynamic bidirectional electrostatic network that connects the phenol pocket to two distal regions, namely the “G-hinge” and the “allosteric linker.” While the G-hinge influences the entry of phenol into the pocket, the allosteric linker passes the signal to the downstream ATPase domain. We show that phenol binding induces a rewiring of the electrostatic connections by eliciting dynamic allostery and demonstrates that perturbation of the core relay residues results in a complete loss of ATPase stimulation. Furthermore, we found a mutation of the G-hinge, ∼20 Å from the phenol pocket, promotes altered flexibility by shifting the pattern of conformational states accessed, leading to a protein with 7-fold enhanced phenol binding ability and enhanced transcriptional activation. Finally, we conducted a global analysis that illustrates that dynamic allostery-driven conserved community networks are universal and evolutionarily conserved across species. Taken together, these results provide insights into the mechanisms of dynamic allostery-mediated conformational changes in NtrC sensor proteins.

CD8 T cell function and cross-reactivity explored by stepwise increased peptide-HLA versus TCR affinity

Recruitment and activation of CD8 T cells occur through specific triggering of T cell receptor (TCR) by peptide-bound human leucocyte antigen (HLA) ligands. Within the generated trimeric TCR-peptide:HLA complex, the molecular binding affinities between peptide and HLA, and between TCR and peptide:HLA both impact T cell functional outcomes. However, how their individual and combined effects modulate immunogenicity and overall T cell responsiveness has not been investigated systematically. Here, we established two panels of human tumor peptide variants differing in their affinity to HLA. For precise characterization, we developed the “blue peptide assay”, an upgraded cell-based approach to measure the peptide:HLA affinity. These peptide variants were then used to investigate the cross-reactivity of tumor antigen-specific CD8 T cell clonotypes derived from blood of cancer patients after vaccination with either the native or an affinity-optimized Melan-A/MART-1 epitope, or isolated from tumor infiltrated lymph nodes (TILNs). Vaccines containing the native tumor epitope generated T cells with better functionality, and superior cross-reactivity against potential low affinity escape epitopes, as compared to T cells induced by vaccines containing an HLA affinity-optimized epitope. Comparatively, Melan-A/MART-1-specific TILN cells displayed functional and cross-reactive profiles that were heterogeneous and clonotype-dependent. Finally, we took advantage of a collection of T cells expressing affinity-optimized NY-ESO-1-specific TCRs to interrogate the individual and combined impact of peptide:HLA and TCR-pHLA affinities on overall CD8 T cell responses. We found profound and distinct effects of both biophysical parameters, with additive contributions and absence of hierarchical dominance. Altogether, the biological impact of peptide:HLA and TCR-pHLA affinities on T cell responses was carefully dissected in two antigenic systems, frequently targeted in human cancer immunotherapy. Our technology and stepwise comparison open new insights into the rational design and selection of vaccine-associated tumor-specific epitopes and highlight the functional and cross-reactivity profiles that endow T cells with best tumor control capacity.

Peptide-dependent tuning of major histocompatibility complex motional properties and the consequences for cellular immunity

T cell receptors (TCRs) and other receptors of the immune system recognize peptides presented by class I or class II major histocompatibility complex (MHC) proteins. Although we generally distinguish between the MHC protein and its peptide, at an atomic level the two form a structural composite, which allows peptides to influence MHC properties and vice versa. One consequence is the peptide-dependent tuning of MHC structural dynamics, which contributes to protein structural adaptability and influences how receptors identify and bind targets. Peptide-dependent tuning of MHC protein dynamics can impact processes such as antigenicity, TCR cross-reactivity, and T cell repertoire selection. Motional tuning extends beyond the binding groove, influencing peptide selection and exchange, as well as interactions with other immune receptors. Here, we review recent findings showing how peptides can affect the dynamic and adaptable nature of MHC proteins. We highlight consequences for immunity and demonstrate how MHC proteins have evolved to be highly sensitive dynamic reporters, with broad immunological consequences.

Dynamics of peptide loading into major histocompatibility complex class I molecules chaperoned by TAPBPR

Major histocompatibility complex class I (MHC-I) molecules display antigenic peptides on the cell surface for T cell receptor scanning, thereby activating the immune response. Peptide loading into MHC-I molecules is thus a critical step during the antigen presentation process. Chaperone TAP-binding protein related (TAPBPR) plays a critical role in promoting high-affinity peptide loading into MHC-I, by discriminating against the low-affinity ones. However, the complete peptide loading dynamics into TAPBPR-bound MHC-I is still elusive. Here, we constructed kinetic network models based on hundreds of short-time MD simulations with an aggregated simulation time of ∼21.7 μs, and revealed, at atomic level, four key intermediate states of one antigenic peptide derived from melanoma-associated MART-1/Melan-A protein during its loading process into TAPBPR-bound MHC-I. We find that the TAPBPR binding at the MHC-I pocket-F can substantially reshape the distant pocket-B via allosteric regulations, which in turn promotes the following peptide N-terminal loading. Intriguingly, the partially loaded peptide could profoundly weaken the TAPBPR-MHC stability, promoting the dissociation of the TAPBPR scoop-loop (SL) region from the pocket-F to a more solvent-exposed conformation. Structural inspections further indicate that the peptide loading could remotely affect the SL binding site through both allosteric perturbations and direct contacts. In addition, another structural motif of TAPBPR, the jack hairpin region, was also found to participate in mediating the peptide editing. Our study sheds light on the detailed molecular mechanisms underlying the peptide loading process into TAPBPR-bound MHC-I and pinpoints the key structural factors responsible for dictating the peptide-loading dynamics.

Unconventional modes of peptide-HLA-I presentation change the rules of TCR engagement

The intracellular proteome of virtually every nucleated cell in the body is continuously presented at the cell surface via the human leukocyte antigen class I (HLA-I) antigen processing pathway. This pathway classically involves proteasomal degradation of intracellular proteins into short peptides that can be presented by HLA-I molecules for interrogation by T-cell receptors (TCRs) expressed on the surface of CD8+ T cells. During the initiation of a T-cell immune response, the TCR acts as the T cell's primary sensor, using flexible loops to mould around the surface of the pHLA-I molecule to identify foreign or dysregulated antigens. Recent findings demonstrate that pHLA-I molecules can also be highly flexible and dynamic, altering their shape according to minor polymorphisms between different HLA-I alleles, or interactions with different peptides. These flexible presentation modes have important biological consequences that can, for example, explain why some HLA-I alleles offer greater protection against HIV, or why some cancer vaccine approaches have been ineffective. This review explores how these recent findings redefine the rules for peptide presentation by HLA-I molecules and extend our understanding of the molecular mechanisms that govern TCR-mediated antigen discrimination.

l- to d-Amino Acid Substitution in the Immunodominant LCMV-Derived Epitope gp33 Highlights the Sensitivity of the TCR Recognition Mechanism for the MHC/Peptide Structure and Dynamics

Presentation of pathogen-derived epitopes by major histocompatibility complex I (MHC-I) can lead to the activation and expansion of specific CD8⁺ T cell clones, eventually resulting in the destruction of infected target cells. Altered peptide ligands (APLs), designed to elicit immunogenicity toward a wild-type peptide, may affect the overall stability of MHC-I/peptide (pMHC) complexes and modulate the recognition by T cell receptors (TCR). Previous works have demonstrated that proline substitution at position 3 (p3P) of different MHC-restricted epitopes, including the immunodominant LCMV-derived epitope gp33 and escape variants, may be an effective design strategy to increase epitope immunogenicity. These studies hypothesized that the p3P substitution increases peptide rigidity, facilitating TCR binding. Here, molecular dynamics simulations indicate that the p3P modification rigidifies the APLs in solution predisposing them for the MHC-I loading as well as once bound to H-2D^b, predisposing them for TCR binding. Our results also indicate that peptide position 6, key for interaction of H-2D^b/gp33 with the TCR P14, takes a suboptimal conformation before as well as after binding to the TCR. Analyses of H-2D^b in complex with APLs, in which position 6 was subjected to an l- to d-amino acid modification, revealed small conformational changes and comparable pMHC thermal stability. However, the l- to d-modification reduced significantly the binding to P14 even in the presence of the p3P modification. Our combined data highlight the sensitivity of the TCR for the conformational dynamics of pMHC and provide further tools to dissect and modulate TCR binding and immunogenicity via APLs.

N-Terminal-Driven Binding Mechanism of an Antigen Peptide to Human Leukocyte Antigen-A*2402 Elucidated by Multicanonical Molecular Dynamic-Based Dynamic Docking and Path Sampling Simulations

We have applied our advanced multicanonical molecular dynamics (McMD)-based dynamic docking methodology to investigate the binding mechanism of an HIV-1 Nef protein epitope to the Asian-dominant allele human leukocyte antigen (HLA)-A*2402. Even though pMHC complex formation [between a Major histocompatibility complex (MHC) class I molecule, which is encoded by an HLA allele, and an antigen peptide] is one of the fundamental processes of the adaptive human immune response, its binding mechanism has not yet been well studied, partially due to the high allelic variation of HLAs in the population. We have used our developed McMD-based dynamic docking method and have successfully reproduced the native complex structure, which is located near the free energy global minimum. Subsequent path sampling MD simulations elucidated the atomic details of the binding process and indicated that the peptide binding is initially driven by the highly positively charged N-terminus of the peptide that is attracted to the various negatively charged residues on the MHC molecule's surface. Upon nearing the pocket, the second tyrosine residue of the peptide anchors the peptide by strongly binding to the B-site of the MHC molecule via hydrophobic driven interactions, resulting in a very strong bound complex structure. Our methodology can be effectively used to predict the bound complex structures between MHC molecules and their antigens to study their binding mechanism in close detail, which would help with the development of new vaccines against cancers, as well as viral infections such as HIV and COVID-19.

An unbiased approach to defining bona fide cancer neoepitopes that elicit immune-mediated cancer rejection

Identification of neoepitopes that are effective in cancer therapy is a major challenge in creating cancer vaccines. Here, using an entirely unbiased approach, we queried all possible neoepitopes in a mouse cancer model and asked which of those are effective in mediating tumor rejection and, independently, in eliciting a measurable CD8 response. This analysis uncovered a large trove of effective anticancer neoepitopes that have strikingly different properties from conventional epitopes and suggested an algorithm to predict them. It also revealed that our current methods of prediction discard the overwhelming majority of true anticancer neoepitopes. These results from a single mouse model were validated in another antigenically distinct mouse cancer model and are consistent with data reported in human studies. Structural modeling showed how the MHC I-presented neoepitopes had an altered conformation, higher stability, or increased exposure to T cell receptors as compared with the unmutated counterparts. T cells elicited by the active neoepitopes identified here demonstrated a stem-like early dysfunctional phenotype associated with effective responses against viruses and tumors of transgenic mice. These abundant anticancer neoepitopes, which have not been tested in human studies thus far, can be exploited for generation of personalized human cancer vaccines.

The Structure of a Peptide-Loaded Shark MHC Class I Molecule Reveals Features of the Binding between β2-Microglobulin and H Chain Conserved in Evolution

Cartilaginous fish are the most primitive extant species with MHC molecules. Using the nurse shark, the current study is, to the best of our knowledge, the first to present a peptide-loaded MHC class I (pMHC-I) structure for this class of animals. The overall structure was found to be similar between cartilaginous fish and bony animals, showing remarkable conservation of interactions between the three pMHC-I components H chain, β₂-microglobulin (β₂-m), and peptide ligand. In most previous studies, relatively little attention was given to the details of binding between the H chain and β₂-m, and our study provides important new insights. A pronounced conserved feature involves the insertion of a large β₂-m F56+W60 hydrophobic knob into a pleat of the β-sheet floor of the H chain α1α2 domain, with the knob being surrounded by conserved residues. Another conserved feature is a hydrogen bond between β₂-m Y10 and a proline in the α3 domain of the H chain. By alanine substitution analysis, we found that the conserved β₂-m residues Y10, D53, F56, and W60-each binding the H chain-are required for stable pMHC-I complex formation. For the β₂-m residues Y10 and F56, such observations have not been reported before. The combined data indicate that for stable pMHC-I complex formation β₂-m should not only bind the α1α2 domain but also the α3 domain. Knowing the conserved structural features of pMHC-I should be helpful for future elucidations of the mechanisms of pMHC-I complex formation and peptide editing.

Structural Comparison Between MHC Classes I and II; in Evolution, a Class-II-Like Molecule Probably Came First

Structures of peptide-loaded major histocompatibility complex class I (pMHC-I) and class II (pMHC-II) complexes are similar. However, whereas pMHC-II complexes include similar-sized IIα and IIβ chains, pMHC-I complexes include a heavy chain (HC) and a single domain molecule β₂-microglobulin (β₂-m). Recently, we elucidated several pMHC-I and pMHC-II structures of primitive vertebrate species. In the present study, a comprehensive comparison of pMHC-I and pMHC-II structures helps to understand pMHC structural evolution and supports the earlier proposed—though debated—direction of MHC evolution from class II-type to class I. Extant pMHC-II structures share major functional characteristics with a deduced MHC-II-type homodimer ancestor. Evolutionary establishment of pMHC-I presumably involved important new functions such as (i) increased peptide selectivity by binding the peptides in a closed groove (ii), structural amplification of peptide ligand sequence differences by binding in a non-relaxed fashion, and (iii) increased peptide selectivity by syngeneic heterotrimer complex formation between peptide, HC, and β₂-m. These new functions were associated with structures that since their establishment in early pMHC-I have been very well conserved, including a shifted and reorganized P1 pocket (aka A pocket), and insertion of a β₂-m hydrophobic knob into the peptide binding domain β-sheet floor. A comparison between divergent species indicates better sequence conservation of peptide binding domains among MHC-I than among MHC-II, agreeing with more demanding interactions within pMHC-I complexes. In lungfishes, genes encoding fusions of all MHC-IIα and MHC-IIβ extracellular domains were identified, and although these lungfish genes presumably derived from classical MHC-II, they provide an alternative mechanistic hypothesis for how evolution from class II-type to class I may have occurred.

TAPBPR promotes antigen loading on MHC-I molecules using a peptide trap

Chaperones Tapasin and TAP-binding protein related (TAPBPR) perform the important functions of stabilizing nascent MHC-I molecules (chaperoning) and selecting high-affinity peptides in the MHC-I groove (editing). While X-ray and cryo-EM snapshots of MHC-I in complex with TAPBPR and Tapasin, respectively, have provided important insights into the peptide-deficient MHC-I groove structure, the molecular mechanism through which these chaperones influence the selection of specific amino acid sequences remains incompletely characterized. Based on structural and functional data, a loop sequence of variable lengths has been proposed to stabilize empty MHC-I molecules through direct interactions with the floor of the groove. Using deep mutagenesis on two complementary expression systems, we find that important residues for the Tapasin/TAPBPR chaperoning activity are located on a large scaffolding surface, excluding the loop. Conversely, loop mutations influence TAPBPR interactions with properly conformed MHC-I molecules, relevant for peptide editing. Detailed biophysical characterization by solution NMR, ITC and FP-based assays shows that the loop hovers above the MHC-I groove to promote the capture of incoming peptides. Our results suggest that the longer loop of TAPBPR lowers the affinity requirements for peptide selection to facilitate peptide loading under conditions and subcellular compartments of reduced ligand concentration, and to prevent disassembly of high-affinity peptide-MHC-I complexes that are transiently interrogated by TAPBPR during editing.

Venus flytrap or pas de trois? The dynamics of MHC class I molecules

The peptide binding site of major histocompatibility complex (MHC) class I molecules is natively unfolded when devoid of peptides. Peptide binding stabilizes the structure and slows the dynamics, but peptide-specific and subtype-specific motions influence, and are influenced by, interaction with assembly chaperones, the T cell receptor, and other class I-binding proteins. The molecular mechanisms of cooperation between peptide, class I heavy chain, and beta-2 microglobulin are insufficiently known but are being elucidated by nuclear magnetic resonance and other modern methods. It appears that micropolymorphic clusters of charged amino acids, often hidden in the molecule interior, determine the dynamics and thus chaperone dependence, cellular fate, and disease association of class I.

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.

Dynamic allostery controls the peptide sensitivity of the Ly49C natural killer receptor

Using a variety of activating and inhibitory receptors, natural killer (NK) cells protect against disease by eliminating cells that have downregulated class I major histocompatibility complex (MHC) proteins, such as in response to cell transformation or viral infection. The inhibitory murine NK receptor Ly49C specifically recognizes the class I MHC protein H-2K^b. Unusual among NK receptors, Ly49C exhibits a peptide-dependent sensitivity to H-2K^b recognition, which has not been explained despite detailed structural studies. To gain further insight into Ly49C peptide sensitivity, we examined Ly49C recognition biochemically and through the lens of dynamic allostery. We found that the peptide sensitivity of Ly49C arises through small differences in H-2K^b-binding affinity. Although molecular dynamics simulations supported a role for peptide-dependent protein dynamics in producing these differences in binding affinity, calorimetric measurements indicated an enthalpically as opposed to entropically driven process. A quantitative linkage analysis showed that this emerges from peptide-dependent dynamic tuning of electrostatic interactions across the Ly49C–H-2K^b interface. We propose a model whereby different peptides alter the flexibility of H-2K^b, which in turn changes the strength of electrostatic interactions across the protein–protein interface. Our results provide a quantitative assessment of how peptides alter Ly49C-binding affinity, suggest the underlying mechanism, and demonstrate peptide-driven allostery at work in class I MHC proteins. Lastly, our model provides a solution for how dynamic allostery could impact binding of some, but not all, class I MHC partners depending on the structural and chemical composition of the interfaces.

MHC I assembly and peptide editing - chaperones, clients, and molecular plasticity in immunity

Peptides presented on MHC I molecules allow the immune system to detect diseased cells. The displayed peptides typically stem from proteasomal degradation of cytoplasmic proteins and are translocated into the ER lumen where they are trimmed and loaded onto MHC I. Peptide translocation is carried out by the transporter associated with antigen processing, which forms the central building block of a dynamic assembly called the peptide-loading complex (PLC). By coordinating peptide transfer with MHC I loading and peptide optimization, the PLC is a linchpin in the adaptive immune system. Peptide loading and optimization is catalyzed by the PLC component tapasin and the PLC-independent TAPBPR, two MHC I-dedicated enzymes chaperoning empty or suboptimally loaded MHC I and selecting stable peptide-MHC I complexes in a process called peptide editing or proofreading. Recent structural and functional studies of peptide editing have dramatically improved our understanding of this pivotal event in antigen processing/presentation. This review is dedicated to Vincenzo Cerundolo (1959-2020) for his pioneering work in the field of antigen processing/presentation.

Common allotypes of ER aminopeptidase 1 have substrate-dependent and highly variable enzymatic properties

Polymorphic variation of immune system proteins can drive variability of individual immune responses. Endoplasmic reticulum aminopeptidase 1 (ERAP1) generates antigenic peptides for presentation by major histocompatibility complex class I molecules. Coding SNPs in ERAP1 have been associated with predisposition to inflammatory rheumatic disease and shown to affect functional properties of the enzyme, but the interplay between combinations of these SNPs as they exist in allotypes has not been thoroughly explored. We used phased genotype data to estimate ERAP1 allotype frequency in 2504 individuals across five major human populations, generated highly pure recombinant enzymes corresponding to the ten most common ERAP1 allotypes, and systematically characterized their in vitro enzymatic properties. We find that ERAP1 allotypes possess a wide range of enzymatic activities, up to 60-fold, whose ranking is substrate dependent. Strikingly, allotype 10, previously associated with Behçet’s disease, is consistently a low-activity outlier, suggesting that a significant percentage of individuals carry a subactive ERAP1 gene. Enzymatic analysis revealed that ERAP1 allotypes can differ in both catalytic efficiency and substrate affinity, differences that can change intermediate accumulation in multistep trimming reactions. Alterations in efficacy of an allosteric inhibitor that targets the regulatory site suggest that allotypic variation influences the communication between the regulatory and the active site. Our work defines the wide landscape of ERAP1 activity in human populations and demonstrates how common allotypes can induce substrate-dependent variability in antigen processing, thus contributing, in synergy with major histocompatibility complex haplotypes, to immune response variability and predisposition to chronic inflammatory conditions.

Structurally silent peptide anchor modifications allosterically modulate T cell recognition in a receptor-dependent manner

Presentation of peptides by class I MHC proteins underlies T cell immune responses to pathogens and cancer. The association between peptide binding affinity and immunogenicity has led to the engineering of modified peptides with improved MHC binding, with the hope that these peptides would be useful for eliciting cross-reactive immune responses directed toward their weak binding, unmodified counterparts. Increasing evidence, however, indicates that T cell receptors (TCRs) can perceive such anchor-modified peptides differently than wild-type (WT) peptides, although the scope of discrimination is unclear. We show here that even modifications at primary anchors that have no discernible structural impact can lead to substantially stronger or weaker T cell recognition depending on the TCR. Surprisingly, the effect of peptide anchor modification can be sensed by a TCR at regions distant from the site of modification, indicating a through-protein mechanism in which the anchor residue serves as an allosteric modulator for TCR binding. Our findings emphasize caution in the use and interpretation of results from anchor-modified peptides and have implications for how anchor modifications are accounted for in other circumstances, such as predicting the immunogenicity of tumor neoantigens. Our data also highlight an important need to better understand the highly tunable dynamic nature of class I MHC proteins and the impact this has on various forms of immune recognition.

Markov state modeling reveals alternative unbinding pathways for peptide-MHC complexes

Peptide binding to major histocompatibility complexes (MHCs) is a central component of the immune system, and understanding the mechanism behind stable peptide-MHC binding will aid the development of immunotherapies. While MHC binding is mostly influenced by the identity of the so-called anchor positions of the peptide, secondary interactions from nonanchor positions are known to play a role in complex stability. However, current MHC-binding prediction methods lack an analysis of the major conformational states and might underestimate the impact of secondary interactions. In this work, we present an atomically detailed analysis of peptide-MHC binding that can reveal the contributions of any interaction toward stability. We propose a simulation framework that uses both umbrella sampling and adaptive sampling to generate a Markov state model (MSM) for a coronavirus-derived peptide (QFKDNVILL), bound to one of the most prevalent MHC receptors in humans (HLA-A24:02). While our model reaffirms the importance of the anchor positions of the peptide in establishing stable interactions, our model also reveals the underestimated importance of position 4 (p4), a nonanchor position. We confirmed our results by simulating the impact of specific peptide mutations and validated these predictions through competitive binding assays. By comparing the MSM of the wild-type system with those of the D4A and D4P mutations, our modeling reveals stark differences in unbinding pathways. The analysis presented here can be applied to any peptide-MHC complex of interest with a structural model as input, representing an important step toward comprehensive modeling of the MHC class I pathway.

An Engineered T Cell Receptor Variant Realizes the Limits of Functional Binding Modes

T cell receptors (TCRs) orchestrate cellular immunity by recognizing peptides presented by a range of major histocompatibility complex (MHC) proteins. Naturally occurring TCRs bind the composite peptide/MHC surface, recognizing peptides that are structurally and chemically compatible with the TCR binding site. Here we describe a molecularly evolved TCR variant that binds the human class I MHC protein HLA-A2 independent of the bound peptide, achieved by a drastic perturbation of the TCR binding geometry that places the molecule far from the peptide binding groove. This unique geometry is unsupportive of normal T cell signaling. A substantial divergence between affinity measurements in solution and in two dimensions between proximal cell membranes leads us to attribute the lack of signaling to steric hindrance that limits binding in the confines of a cell-cell interface. Our results provide an example of how receptor binding geometry can impact T cell function and provide further support for the view that germline-encoded residues in TCR binding loops evolved to drive productive TCR recognition and signaling.

Atomistic structure and dynamics of the human MHC-I peptide-loading complex

The major histocompatibility complex class-I (MHC-I) peptide-loading complex (PLC) is a cornerstone of the human adaptive immune system, being responsible for processing antigens that allow killer T cells to distinguish between healthy and compromised cells. Based on a recent low-resolution cryo-electron microscopy (cryo-EM) structure of this large membrane-bound protein complex, we report an atomistic model of the PLC and study its conformational dynamics on the multimicrosecond time scale using all-atom molecular dynamics (MD) simulations in an explicit lipid bilayer and water environment (1.6 million atoms in total). The PLC has a layered structure, with two editing modules forming a flexible protein belt surrounding a stable, catalytically active core. Tapasin plays a central role in the PLC, stabilizing the MHC-I binding groove in a conformation reminiscent of antigen-loaded MHC-I. The MHC-I-linked glycan steers a tapasin loop involved in peptide editing toward the binding groove. Tapasin conformational dynamics are also affected by calreticulin through a conformational selection mechanism that facilitates MHC-I recruitment into the complex.

Capturing the Flexibility of a Protein-Ligand Complex: Binding Free Energies from Different Enhanced Sampling Techniques

Enhanced sampling techniques are a promising approach to obtain reliable binding free-energy profiles for flexible protein-ligand complexes from molecular dynamics (MD) simulations. To put four popular enhanced sampling techniques to a biologically relevant and challenging test, we studied the partial dissociation of an antigenic peptide from the Major Histocompatibility Complex I (MHC I) HLA-B*35:01 to systematically investigate the performance of umbrella sampling (US), replica exchange with solute tempering 2 (REST2), bias exchange umbrella sampling (BEUS, or replica-exchange umbrella sampling), and well-tempered metadynamics (MTD). With regard to the speed of sampling and convergence, the peptide-MHC I complex (pMHC I) under study showcases intrinsic strengths and weaknesses of the four enhanced sampling techniques used. We found that BEUS can best handle the sampling challenges that arise from the coexistence of an enthalpically and an entropically stabilized free-energy minimum in the pMHC I under study. These findings might also be relevant for other flexible biomolecular systems with competing enthalpically and entropically stabilized minima.

T cell receptor interactions with human leukocyte antigen govern indirect peptide selectivity for the cancer testis antigen MAGE-A4

T cell-mediated immunity is governed primarily by T cell receptor (TCR) recognition of peptide-human leukocyte antigen (pHLA) complexes and is essential for immunosurveillance and disease control. This interaction is generally stabilized by interactions between the HLA surface and TCR germline-encoded complementarity-determining region (CDR) loops 1 and 2, whereas peptide selectivity is guided by direct interactions with the TCR CDR3 loops. Here, we solved the structure of a newly identified TCR in complex with a clinically relevant peptide derived from the cancer testis antigen melanoma antigen-A4 (MAGE-A4). The TCR bound pHLA in a position shifted toward the peptide's N terminus. This enabled the TCR to achieve peptide selectivity via an indirect mechanism, whereby the TCR sensed the first residue of the peptide through HLA residue Trp-167, which acted as a tunable gateway. Amino acid substitutions at peptide position 1 predicted to alter the HLA Trp-167 side-chain conformation abrogated TCR binding, indicating that this indirect binding mechanism is essential for peptide recognition. These findings extend our understanding of the molecular rules that underpin antigen recognition by TCRs and have important implications for the development of TCR-based therapies.

Peptide cargo tunes a network of correlated motions in human leucocyte antigens

Most biomolecular interactions are typically thought to increase the (local) rigidity of a complex, for example, in drug‐target binding. However, detailed analysis of specific biomolecular complexes can reveal a more subtle interplay between binding and rigidity. Here, we focussed on the human leucocyte antigen (HLA), which plays a crucial role in the adaptive immune system by presenting peptides for recognition by the αβ T‐cell receptor (TCR). The role that the peptide plays in tuning HLA flexibility during TCR recognition is potentially crucial in determining the functional outcome of an immune response, with obvious relevance to the growing list of immunotherapies that target the T‐cell compartment. We have applied high‐pressure/temperature perturbation experiments, combined with molecular dynamics simulations, to explore the drivers that affect molecular flexibility for a series of different peptide–HLA complexes. We find that different peptide sequences affect peptide–HLA flexibility in different ways, with the peptide cargo tuning a network of correlated motions throughout the pHLA complex, including in areas remote from the peptide‐binding interface, in a manner that could influence T‐cell antigen discrimination.

Multifunctional Chaperone and Quality Control Complexes in Adaptive Immunity

The fundamental process of adaptive immunity relies on the differentiation of self from nonself. Nucleated cells are continuously monitored by effector cells of the immune system, which police the peptide status presented via cell surface molecules. Recent integrative structural approaches have provided insights toward our understanding of how sophisticated cellular machineries shape such hierarchical immune surveillance. Biophysical and structural achievements were invaluable for defining the interconnection of many key factors during antigen processing and presentation, and helped to solve several conundrums that persisted for many years. In this review, we illuminate the numerous quality control machineries involved in different steps during the maturation of major histocompatibility complex class I (MHC I) proteins, from their synthesis in the endoplasmic reticulum to folding and trafficking via the secretory pathway, optimization of antigenic cargo, final release to the cell surface, and engagement with their cognate receptors on cytotoxic T lymphocytes.

Molecular mechanism of Be2+-ion binding to HLA-DP2: tetrahedral coordination, conformational changes and multi-ion binding

Be small and positive: the smaller size and higher charge of the Be²⁺-ion results in strong binding between the M2 peptide and the β-chain of HLA-DP2, which induces conformational changes at the periphery suitable for TCR binding.

Metal-triggered conformational reorientation of a self-peptide bound to a disease-associated HLA-B*27 subtype

Conformational changes of major histocompatibility complex (MHC) antigens have the potential to be recognized by T cells and may arise from polymorphic variation of the MHC molecule, the binding of modifying ligands, or both. Here, we investigated whether metal ions could affect allele-dependent structural variation of the two minimally distinct human leukocyte antigen (HLA)-B*27:05 and HLA-B*27:09 subtypes, which exhibit differential association with the rheumatic disease ankylosing spondylitis (AS). We employed NMR spectroscopy and X-ray crystallography coupled with ensemble refinement to study the AS-associated HLA-B*27:05 subtype and the AS-nonassociated HLA-B* 27:09 in complex with the self-peptide pVIPR (RRKWRRWHL). Both techniques revealed that pVIPR exhibits a higher degree of flexibility when complexed with HLA-B*27:05 than with HLA-B*27:09. Furthermore, we found that the binding of the metal ion Cu²⁺ or Ni²⁺, but not Mn²⁺, Zn²⁺, or Hg²⁺, affects the structure of a pVIPR-bound HLA-B*27 molecule in a subtype-dependent manner. In HLA-B*27:05, the metals triggered conformational reorientations of pVIPR, but no such structural changes were observed in the HLA-B*27:09 subtype, with or without bound metal ion. These observations provide the first demonstration that not only major histocompatibility complex class II, but also class I, molecules can undergo metal ion-induced conformational alterations. Our findings suggest that metals may have a role in triggering rheumatic diseases such as AS and also have implications for the molecular basis of metal-induced hypersensitivities and allergies.