The Thermodynamic Mechanism of Peptide-MHC Class II Complex Formation Is a Determinant of Susceptibility to HLA-DM

Predicting reverse-bound peptide conformations in MHC Class II with PANDORA

Recent discoveries have transformed our understanding of peptide binding in Major Histocompatibility Complex (MHC) molecules, showing that peptides, for some MHC class II alleles, can bind in a reverse orientation (C-terminus to N-terminus) and can still effectively activate CD4+ T cells. These finding challenges established concepts of immune recognition and suggests new pathways for therapeutic intervention, such as vaccine design. We present an updated version of PANDORA, which, to the best of our knowledge, is the first tool capable of modeling reversed-bound peptides. Modeling these peptides presents a unique challenge due to the limited structural data available for these orientations in existing databases. PANDORA has overcome this challenge through integrative modeling using algorithmically reversed peptides as templates. We have validated the new PANDORA feature through two targeted experiments, achieving an average backbone binding-core L-RMSD value of 0.63 Å. Notably, it maintained low RMSD values even when using templates from different alleles and peptide sequences. Our results suggest that PANDORA will be an invaluable resource for the immunology community, aiding in the development of targeted immunotherapies and vaccine design.

MHC-II dynamics are maintained in HLA-DR allotypes to ensure catalyzed peptide exchange

Presentation of antigenic peptides by major histocompatibility complex class II (MHC-II) proteins determines T helper cell reactivity. The MHC-II genetic locus displays a large degree of allelic polymorphism influencing the peptide repertoire presented by the resulting MHC-II protein allotypes. During antigen processing, the human leukocyte antigen (HLA) molecule HLA-DM (DM) encounters these distinct allotypes and catalyzes exchange of the placeholder peptide CLIP by exploiting dynamic features of MHC-II. Here, we investigate 12 highly abundant CLIP-bound HLA-DRB1 allotypes and correlate dynamics to catalysis by DM. Despite large differences in thermodynamic stability, peptide exchange rates fall into a target range that maintains DM responsiveness. A DM-susceptible conformation is conserved in MHC-II molecules, and allosteric coupling between polymorphic sites affects dynamic states that influence DM catalysis. As exemplified for rheumatoid arthritis, we postulate that intrinsic dynamic features of peptide-MHC-II complexes contribute to the association of individual MHC-II allotypes with autoimmune disease.

Incorporating the Molecular Mimicry of Environmental Antigens into the Causality of Autoimmune Hepatitis

Molecular mimicry between foreign and self-antigens has been implicated as a cause of autoimmune hepatitis in experimental models and cross-reacting antibodies in patients. This review describes the experimental and clinical evidence for molecular mimicry as a cause of autoimmune hepatitis, indicates the limitations and uncertainties of this premise, and encourages investigations that assess diverse environmental antigens as sources of disease-relevant molecular mimics. Pertinent articles were identified in PubMed using multiple search phrases. Several pathogens have linear or conformational epitopes that mimic the self-antigens of autoimmune hepatitis. The occurrence of an acute immune-mediated hepatitis after vaccination for severe acute respiratory syndrome (SARS)-associated coronavirus 2 (SARS-CoV-2) has suggested that vaccine-induced peptides may mimic disease-relevant tissue antigens. The intestinal microbiome is an under-evaluated source of gut-derived antigens that could also engage in molecular mimicry. Chaperone molecules may enhance the pathogenicity of molecular mimics, and they warrant investigation. Molecular mimics of immune dominant epitopes within cytochrome P450 IID6, the autoantigen most closely associated with autoimmune hepatitis, should be sought in diverse environmental antigens and assessed for pathogenicity. Avoidance strategies, dietary adjustments, vaccine improvement, and targeted manipulation of the intestinal microbiota may emerge as therapeutic possibilities. In conclusion, molecular mimicry may be a missing causality of autoimmune hepatitis. Molecular mimics of key immune dominant epitopes of disease-specific antigens must be sought in diverse environmental antigens. The ubiquity of molecular mimicry compels rigorous assessments of peptide mimics for immunogenicity and pathogenicity in experimental models. Molecular mimicry may complement epigenetic modifications as causative mechanisms of autoimmune hepatitis.

Partnering for the major histocompatibility complex class II and antigenic determinant requires flexibility and chaperons

Cytotoxic, or helper T cells recognize antigen via T cell receptors (TCRs) that can see their target antigen as short sequences of peptides bound to the groove of proteins of major histocompatibility complex (MHC) class I, and class II respectively. For MHC class II epitope selection from exogenous pathogens or self-antigens, participation of several accessory proteins, molecular chaperons, processing enzymes within multiple vesicular compartments is necessary. A major contributing factor is the MHC class II structure itself that uniquely offers a dynamic and flexible groove essential for epitope selection. In this review, I have taken a historical perspective focusing on the flexibility of the MHC II molecules as the driving force in determinant selection and interactions with the accessory molecules in antigen processing, HLA-DM and HLA-DO.

How Does B Cell Antigen Presentation Affect Memory CD4 T Cell Differentiation and Longevity?

Dendritic cells are the antigen presenting cells that process antigens effectively and prime the immune system, a characteristic that have gained them the spotlights in recent years. B cell antigen presentation, although less prominent, deserves equal attention. B cells select antigen experienced CD4 T cells to become memory and initiate an orchestrated genetic program that maintains memory CD4 T cells for life of the individual. Over years of research, we have demonstrated that low levels of antigens captured by B cells during the resolution of an infection render antigen experienced CD4 T cells into a quiescent/resting state. Our studies suggest that in the absence of antigen, the resting state associated with low-energy utilization and proliferation can help memory CD4 T cells to survive nearly throughout the lifetime of mice. In this review we would discuss the primary findings from our lab as well as others that highlight our understanding of B cell antigen presentation and the contributions of the MHC Class II accessory molecules to this outcome. We propose that the quiescence induced by the low levels of antigen presentation might be a mechanism necessary to regulate long-term survival of CD4 memory T cells and to prevent cross-reactivity to autoantigens, hence autoimmunity.

MHCBI: a pipeline for calculating peptide-MHC binding energy using semi-empirical quantum mechanical methods with explicit/implicit solvent models

Experimentally estimating peptide-major histocompatibility complex (pMHC) binding affinity has been quite challenging due to the many receptors and the many potential ligands implicated in it. We have thus proposed a straightforward computational methodology considering the different mechanisms involved in pMHC binding to facilitate studying such receptor-ligand interactions. We have developed a pipeline using semi-empirical quantum mechanical methods for calculating pMHC class I and II molecules' binding energy (BE). This pipeline can systematize the methodology for calculating pMHC system BE, enabling the rational design of T-cell epitopes to be used as pharmaceuticals and vaccines.

Impact of Structural Observables From Simulations to Predict the Effect of Single-Point Mutations in MHC Class II Peptide Binders

The prediction of peptide binders to Major Histocompatibility Complex (MHC) class II receptors is of great interest to study autoimmune diseases and for vaccine development. Most approaches predict the affinities using sequence-based models trained on experimental data and multiple alignments from known peptide substrates. However, detecting activity differences caused by single-point mutations is a challenging task. In this work, we used interactions calculated from simulations to build scoring matrices for quickly estimating binding differences by single-point mutations. We modelled a set of 837 peptides bound to an MHC class II allele, and optimized the sampling of the conformations using the Rosetta backrub method by comparing the results to molecular dynamics simulations. From the dynamic trajectories of each complex, we averaged and compared structural observables for each amino acid at each position of the 9°mer peptide core region. With this information, we generated the scoring-matrices to predict the sign of the binding differences. We then compared the performance of the best scoring-matrix to different computational methodologies that range in computational costs. Overall, the prediction of the activity differences caused by single mutated peptides was lower than 60% for all the methods. However, the developed scoring-matrix in combination with existing methods reports an increase in the performance, up to 86% with a scoring method that uses molecular dynamics.

Comparison of HLA ligand elution data and binding predictions reveals varying prediction performance for the multiple motifs recognized by HLA-DQ2.5

Binding prediction tools are commonly used to identify peptides presented on MHC class II molecules. Recently, a wealth of data in the form of naturally eluted ligands has become available and discrepancies between ligand elution data and binding predictions have been reported. Quantitative metrics for such comparisons are currently lacking. In this study, we assessed how efficiently MHC class II binding predictions can identify naturally eluted peptides, and investigated instances with discrepancies between the two methods in detail. We found that, in general, MHC class II eluted ligands are predicted to bind to their reported restriction element with high affinity. But, for several studies reporting an increased number of ligands that were not predicted to bind, we found that the reported MHC restriction was ambiguous. Additional analyses determined that most of the ligands predicted to not bind, are predicted to bind other co‐expressed MHC class II molecules. For selected alleles, we addressed discrepancies between elution data and binding predictions by experimental measurements and found that predicted and measured affinities correlate well. For DQA1*05:01/DQB1*02:01 (DQ2.5) however, binding predictions did miss several peptides that were determined experimentally to be binders. For these peptides and several known DQ2.5 binders, we determined key residues for conferring DQ2.5 binding capacity, which revealed that DQ2.5 utilizes two different binding motifs, of which only one is predicted effectively. These findings have important implications for the interpretation of ligand elution data and for the improvement of MHC class II binding predictions.

HLA-DM catalytically enhances peptide dissociation by sensing peptide-MHC class II interactions throughout the peptide-binding cleft

Human leukocyte antigen-DM (HLA-DM) is an integral component of the major histocompatibility complex class II (MHCII) antigen-processing and -presentation pathway. HLA-DM shapes the immune system by differentially catalyzing peptide exchange on MHCII molecules, thereby editing the peptide-MHCII (pMHCII) repertoire by imposing a bias on the foreign and self-derived peptide cargos that are presented on the cell surface for immune surveillance and tolerance induction by CD4⁺ T cells. To better understand DM selectivity, here we developed a real-time fluorescence anisotropy assay to delineate the pMHCII intrinsic stability, DM-binding affinity, and catalytic turnover, independent kinetic parameters of HLA-DM enzymatic activity. We analyzed prominent pMHCII contacts by differentiating the kinetic parameters in pMHCII homologs, observing that peptide interactions throughout the MHCII-binding cleft influence both the rate of peptide dissociation from the DM-pMHCII catalytic complex and the binding affinity of HLA-DM for a pMHCII. We show that the intrinsic stability of a pMHCII linearly correlates with DM catalytic turnover, but is nonlinearly correlated with its binding affinity. Surprisingly, interactions at the peptides N terminus up to and including MHCII position one (P1) anchor affected the catalytic turnover, suggesting that the active DM-pMHCII catalytic complex operates on pMHCII complexes with full peptide occupancy. Furthermore, interactions at the peptide C terminus modulated DM-binding affinity, suggesting distal communication between peptide interactions with the MHCII and the DM-pMHCII binding interface. Our results imply an intimate linkage between the DM-pMHCII interface and peptide-MHCII interactions throughout the peptide-binding cleft.

How a Proposed Hypothesis during My PhD Training Shaped My Career

In this memoir-style essay, I have narrated the evolution of my scientific career, as deeply influenced by my PhD training and the mentorship of Professor Eli Sercarz. Starting in his lab, and continuing to my own laboratory, many of the questions we have pursued link in some way to Eli's ideas. In this essay, I have summarized the path that I followed after graduating from his lab and highlight findings along the way. I apologize to my colleagues whose work was not discussed here due to the nature of this review and space limitations.

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

T cell receptor (TCR) recognition of antigenic peptides bound and presented by class I major histocompatibility complex (MHC) proteins underlies the cytotoxic immune response to diseased cells. Crystallographic structures of TCR-peptide/MHC complexes have demonstrated how TCRs simultaneously interact with both the peptide and the MHC protein. However, it is increasingly recognized that, beyond serving as a static platform for peptide presentation, the physical properties of class I MHC proteins are tuned by different peptides in ways that are not always structurally visible. These include MHC protein motions, or dynamics, which are believed to influence interactions with a variety of MHC-binding proteins, including not only TCRs, but other activating and inhibitory receptors as well as components of the peptide loading machinery. Here, we investigated the mechanisms by which peptides tune the dynamics of the common class I MHC protein HLA-A2. By examining more than 50 lengthy molecular dynamics simulations of HLA-A2 presenting different peptides, we identified regions susceptible to dynamic tuning, including regions in the peptide binding domain as well as the distal α3 domain. Further analyses of the simulations illuminated mechanisms by which the influences of different peptides are communicated throughout the protein, and involve regions of the peptide binding groove, the β₂-microglobulin subunit, and the α3 domain. Overall, our results demonstrate that the class I MHC protein is a highly tunable peptide sensor whose physical properties vary considerably with bound peptide. Our data provides insight into the underlying principles and suggest a role for dynamically driven allostery in the immunological function of MHC proteins.

Class II MHC antigen processing in immune tolerance and inflammation

Presentation of peptide antigens by MHC-II proteins is prerequisite to effective CD4 T cell tolerance to self and to recognition of foreign antigens. Antigen uptake and processing pathways as well as expression of the peptide exchange factors HLA-DM and HLA-DO differ among the various professional and non-professional antigen-presenting cells and are modulated by cell developmental state and activation. Recent studies have highlighted the importance of these cell-specific factors in controlling the source and breadth of peptides presented by MHC-II under different conditions. During inflammation, increased presentation of selected self-peptides has implications for maintenance of peripheral tolerance and autoimmunity.

Dynamics Determine Signaling in a Multicomponent System Associated with Rheumatoid Arthritis

Strategies that target multiple components are usually required for treatment of diseases originating from complex biological systems. The multicomponent system consisting of the DR4 major histocompatibility complex type II molecule, the glycopeptide CII259-273 from type II collagen, and a T-cell receptor is associated with development of rheumatoid arthritis (RA). We introduced non-native amino acids and amide bond isosteres into CII259-273 and investigated the effect on binding to DR4 and the subsequent T-cell response. Molecular dynamics simulations revealed that complexes between DR4 and derivatives of CII259-273 were highly dynamic. Signaling in the overall multicomponent system was found to depend on formation of an appropriate number of dynamic intramolecular hydrogen bonds between DR4 and CII259-273, together with the positioning of the galactose moiety of CII259-273 in the DR4 binding groove. Interestingly, the system tolerated modifications at several positions in CII259-273, indicating opportunities to use analogues to increase our understanding of how rheumatoid arthritis develops and for evaluation as vaccines to treat RA.

Exploring the nature of the H-bonds between the human class II MHC protein, HLA-DR1 (DRB*0101) and the influenza virus hemagglutinin peptide, HA306-318, using the quantum theory of atoms in molecules

The nature of the H-bonds between the human protein HLA-DR1 (DRB*0101) and the hemagglutinin peptide HA306-318 has been studied using the Quantum Theory of Atoms in Molecules for the first time. We have found four H-bond groups: one conventional CO··HN bond group and three nonconventional CO··HC, π··HC involving aromatic rings and HN··HC_aliphatic groups. The calculated electron density at the determined H-bond critical points suggests the follow protein pocket binding trend: P1 (2,311) >> P9 (1.109) > P4 (0.950) > P6 (0.553) > P7 (0.213) which agrees and reveal the nature of experimental findings, showing that P1 produces by a long way the strongest binding of the HLA-DR1 human protein molecule with the peptide backbone as consequence of the vast number of H-bonds in the P1 area and at the same time the largest specific binding of the peptide Tyr308 residue with aromatic residues located at the binding groove floor. The present results suggest the topological analysis of the electronic density as a valuable tool that allows a non-arbitrary partition of the pockets binding energy via the calculated electron density at the determined critical points.

Peptide and Peptide-Dependent Motions in MHC Proteins: Immunological Implications and Biophysical Underpinnings

Structural biology of peptides presented by class I and class II MHC proteins has transformed immunology, impacting our understanding of fundamental immune mechanisms and allowing researchers to rationalize immunogenicity and design novel vaccines. However, proteins are not static structures as often inferred from crystallographic structures. Their components move and breathe individually and collectively over a range of timescales. Peptides bound within MHC peptide-binding grooves are no exception and their motions have been shown to impact recognition by T cell and other receptors in ways that influence function. Furthermore, peptides tune the motions of MHC proteins themselves, which impacts recognition of peptide/MHC complexes by other proteins. Here, we review the motional properties of peptides in MHC binding grooves and discuss how peptide properties can influence MHC motions. We briefly review theoretical concepts about protein motion and highlight key data that illustrate immunological consequences. We focus primarily on class I systems due to greater availability of data, but segue into class II systems as the concepts and consequences overlap. We suggest that characterization of the dynamic “energy landscapes” of peptide/MHC complexes and the resulting functional consequences is one of the next frontiers in structural immunology.

Major Histocompatibility Complex (MHC) Class I and MHC Class II Proteins: Conformational Plasticity in Antigen Presentation

Antigen presentation by major histocompatibility complex (MHC) proteins is essential for adaptive immunity. Prior to presentation, peptides need to be generated from proteins that are either produced by the cell's own translational machinery or that are funneled into the endo-lysosomal vesicular system. The prolonged interaction between a T cell receptor and specific pMHC complexes, after an extensive search process in secondary lymphatic organs, eventually triggers T cells to proliferate and to mount a specific cellular immune response. Once processed, the peptide repertoire presented by MHC proteins largely depends on structural features of the binding groove of each particular MHC allelic variant. Additionally, two peptide editors-tapasin for class I and HLA-DM for class II-contribute to the shaping of the presented peptidome by favoring the binding of high-affinity antigens. Although there is a vast amount of biochemical and structural information, the mechanism of the catalyzed peptide exchange for MHC class I and class II proteins still remains controversial, and it is not well understood why certain MHC allelic variants are more susceptible to peptide editing than others. Recent studies predict a high impact of protein intermediate states on MHC allele-specific peptide presentation, which implies a profound influence of MHC dynamics on the phenomenon of immunodominance and the development of autoimmune diseases. Here, we review the recent literature that describe MHC class I and II dynamics from a theoretical and experimental point of view and we highlight the similarities between MHC class I and class II dynamics despite the distinct functions they fulfill in adaptive immunity.

Major Histocompatibility Complex (MHC) Class I and MHC Class II Proteins: Conformational Plasticity in Antigen Presentation

Antigen presentation by major histocompatibility complex (MHC) proteins is essential for adaptive immunity. Prior to presentation, peptides need to be generated from proteins that are either produced by the cell’s own translational machinery or that are funneled into the endo-lysosomal vesicular system. The prolonged interaction between a T cell receptor and specific pMHC complexes, after an extensive search process in secondary lymphatic organs, eventually triggers T cells to proliferate and to mount a specific cellular immune response. Once processed, the peptide repertoire presented by MHC proteins largely depends on structural features of the binding groove of each particular MHC allelic variant. Additionally, two peptide editors—tapasin for class I and HLA-DM for class II—contribute to the shaping of the presented peptidome by favoring the binding of high-affinity antigens. Although there is a vast amount of biochemical and structural information, the mechanism of the catalyzed peptide exchange for MHC class I and class II proteins still remains controversial, and it is not well understood why certain MHC allelic variants are more susceptible to peptide editing than others. Recent studies predict a high impact of protein intermediate states on MHC allele-specific peptide presentation, which implies a profound influence of MHC dynamics on the phenomenon of immunodominance and the development of autoimmune diseases. Here, we review the recent literature that describe MHC class I and II dynamics from a theoretical and experimental point of view and we highlight the similarities between MHC class I and class II dynamics despite the distinct functions they fulfill in adaptive immunity.

Peptidomic analysis of type 1 diabetes associated HLA-DQ molecules and the impact of HLA-DM on peptide repertoire editing

HLA-DM and class II associated invariant chain (Ii) are key cofactors in the MHC class II (MHCII) antigen processing pathway. We used tandem mass spectrometry sequencing to directly interrogate the global impact of DM and Ii on the repertoire of MHCII-bound peptides in human embryonic kidney 293T cells expressing HLA-DQ molecules in the absence or presence of these cofactors. We found that Ii and DM have a major impact on the repertoire of peptides presented by DQ1 and DQ6, with the caveat that this technology is not quantitative. The peptide repertoires of type 1 diabetes (T1D) associated DQ8, DQ2, and DQ8/2 are altered to a lesser degree by DM expression, and these molecules share overlapping features in their peptide binding motifs that are distinct from control DQ1 and DQ6 molecules. Peptides were categorized into DM-resistant, DM-dependent, or DM-sensitive groups based on the mass spectrometry data, and representative peptides were tested in competitive binding assays and peptide dissociation rate experiments with soluble DQ6. Our data support the conclusion that high intrinsic stability of DQ-peptide complexes is necessary but not sufficient to confer resistance to DM editing, and provide candidate parameters that may be useful in predicting the sensitivity of T-cell epitopes to DM editing.

MHC class II complexes sample intermediate states along the peptide exchange pathway

The presentation of peptide-MHCII complexes (pMHCIIs) for surveillance by T cells is a well-known immunological concept in vertebrates, yet the conformational dynamics of antigen exchange remain elusive. By combining NMR-detected H/D exchange with Markov modelling analysis of an aggregate of 275 microseconds molecular dynamics simulations, we reveal that a stable pMHCII spontaneously samples intermediate conformations relevant for peptide exchange. More specifically, we observe two major peptide exchange pathways: the kinetic stability of a pMHCII's ground state defines its propensity for intrinsic peptide exchange, while the population of a rare, intermediate conformation correlates with the propensity of the HLA-DM-catalysed pathway. Helix-destabilizing mutants designed based on our model shift the exchange behaviour towards the HLA-DM-catalysed pathway and further allow us to conceptualize how allelic variation can shape an individual's MHC restricted immune response.

MHCII proteins bind and present both foreign and self-antigens to potentially activate CD4+ T cells via cognate T cell receptors (TCRs) during the adaptive immune response. Here, the authors combine NMR-detected H/D exchange with Markov modelling analysis to shed light on the dynamics of MHCII peptide exchange.

Molecular dynamics at the receptor level of immunodominant myelin oligodendrocyte glycoprotein 35-55 epitope implicated in multiple sclerosis

Multiple Sclerosis (MS) is a common autoimmune disease whereby myelin is destroyed by the immune system. The disease is triggered by the stimulation of encephalitogenic T-cells via the formation of a trimolecular complex between the Human Leukocyte Antigen (HLA), an immunodominant epitope of myelin proteins and T-cell Receptor (TCR). Myelin Oligodendrocyte Glycoprotein (MOG) is located on the external surface of myelin and has been implicated in MS induction. The immunodominant 35-55 epitope of MOG is widely used for in vivo biological evaluation and immunological studies that are related with chronic Experimental Autoimmune Encephalomyelitis (EAE, animal model of MS), inflammatory diseases and MS. In this report, Molecular Dynamics (MD) simulations were used to explore the interactions of MOG35-55 at the receptor level. A detailed mapping of the developed interactions during the creation of the trimolecular complex is reported. This is the first attempt to gain an understanding of the molecular recognition of the MOG35-55 epitope by the HLA and TCR receptors. During the formation of the trimolecular complex, the residues Arg(41) and Arg(46) of MOG35-55 have been confirmed to serve as TCR anchors while Tyr(40) interacts with HLA. The present structural findings indicate that the Arg at positions 41 and 46 is a key residue for the stimulation of the encephalitogenic T-cells.

A step-by-step overview of the dynamic process of epitope selection by major histocompatibility complex class II for presentation to helper T cells

T cell antigen receptors (TCRs) expressed on cytotoxic or helper T cells can only see their specific target antigen as short sequences of peptides bound to the groove of proteins of major histocompatibility complex (MHC) class I, and class II respectively. In addition to the many steps, several participating proteins, and multiple cellular compartments involved in the processing of antigens, the MHC structure, with its dynamic and flexible groove, has perfectly evolved as the underlying instrument for epitope selection. In this review, I have taken a step-by-step, and rather historical, view to describe antigen processing and determinant selection, as we understand it today, all based on decades of intense research by hundreds of laboratories.