Abstract:Molecules of the Major Histocompatibility Complex (MHC) present short protein fragments on the cell surface, an important step in T cell immune recognition. MHC-I molecules process peptides from intracellular proteins; MHC-II molecules act in antigen-presenting cells and present peptides derived from extracellular proteins. Here we show that the sequence-dependent energy landscapes of MHC-peptide binding encode class-specific nonlinearities (epistasis). MHC-I has a smooth landscape with global epistasis; the binding energy is a simple deformation of an underlying linear trait. This form of epistasis enhances the discrimination between strong-binding peptides. In contrast, MHC-II has a rugged landscape with idiosyncratic epistasis: binding depends on detailed amino acid combinations at multiple positions of the peptide sequence. The form of epistasis affects the learning of energy landscapes from training data. For MHC-I, a low-complexity problem, we derive a simple matrix model of binding energies that outperforms current models trained by machine learning. For MHC-II, higher complexity prevents learning by simple regression methods. Epistasis also affects the energy and fitness effects of mutations in antigen-derived peptides (epitopes). In MHC-I, large-effect mutations occur predominantly in anchor positions of strong-binding epitopes. In MHC-II, large effects depend on the background epitope sequence but are broadly distributed over the epitope, generating a bigger target for escape mutations due to loss of presentation. Together, our analysis shows how an energy landscape of protein-protein binding constrains the target of escape mutations from T cell immunity, linking the complexity of the molecular interactions to the dynamics of adaptive immune response. T cell immunity involves the binding of short peptides to the intracellular MHC recognition machinery. Understanding how the peptide-MHC binding energy depends on the peptide sequence is key to computationally predict immune recognition and immune escape evolution, for example, of pathogens and cancer cells. By analyzing experimental binding affinity data, we find that peptide-MHC interactions are described by nonlinear energy landscapes that depend on the recognition pathway: smooth and easy to learn for MHC class I, rugged and difficult to learn for MHC class II. The different forms of nonlinearity—or epistasis—suggest the most suitable strategy to infer these landscapes. In addition, the nonlinearity has direct implications on the energy and fitness of single mutations, generating distinct distributions of possible targets for immune escape. Together, this work establishes links between biophysical origin, nonlinear structure, learnability from data, and biological implications of peptide-MHC interactions, which can be applicable also to other protein interaction landscapes.

Energy landscapes of peptide-MHC binding

Energy landscapes of peptide-MHC binding

Mechano-regulation of Peptide-MHC Class I Conformations Determines TCR Antigen Recognition.

Quantitative analysis of MHC class II peptide exchange reveals pivotal role of peptide association rate

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

The Three-Dimensional Structure of Peptide-MHC Complexes

Trans-Allelic Model for Prediction of Peptide:MHC-II Interactions

Physical epistatic landscape of antibody binding affinity

Plasticity in the contribution of T cell receptor variable region residues to binding of peptide-HLA-A2 complexes

An Efficient Approach to the Accurate Prediction of Mutational Effects in Antigen Binding to the MHC1

MHCII-peptide presentation: an assessment of the state-of-the-art prediction methods

The electrostatic landscape of MHC-peptide binding revealed using inception networks

Ranking of binding and nonbinding peptides to MHC class I molecules using inverse folding approach: Implications for vaccine design

Experimental Structures of Antibody/MHC-I Complexes Reveal Details of Epitopes Overlooked by Computational Prediction

Predicting MHC-I ligands across alleles and species: How far can we go?

Scrutinizing human MHC polymorphism: Supertype analysis using Poisson-Boltzmann electrostatics and clustering

RLpMIEC: High-Affinity Peptide Generation Targeting Major Histocompatibility Complex-I Guided and Interpreted by Interaction Spectrum-Navigated Reinforcement Learning

Machine learning application to predict binding affinity between peptide containing non-canonical amino acids and HLA0201

NetMHCpan-4.0: Improved Peptide–MHC Class I Interaction Predictions Integrating Eluted Ligand and Peptide Binding Affinity Data

Improving generalizability for MHC-I binding peptide predictions through geometric deep learning

APE-Gen2.0: Expanding Rapid Class I Peptide–Major Histocompatibility Complex Modeling to Post-Translational Modifications and Noncanonical Peptide Geometries