Reduced Protein Stability of 11 Pathogenic Missense STXBP1/MUNC18-1 Variants and Improved Disease Prediction.

BACKGROUND: Pathogenic variants in STXBP1/MUNC18-1 cause severe encephalopathies that are among the most common in genetic neurodevelopmental disorders. Different molecular disease mechanisms have been proposed, and pathogenicity prediction is limited. In this study, we aimed to define a generalized disease concept for STXBP1-related disorders and improve prediction. METHODS: A cohort of 11 disease-associated and 5 neutral variants (detected in healthy individuals) were tested in 3 cell-free assays and in heterologous cells and primary neurons. Protein aggregation was tested using gel filtration and Triton X-100 insolubility. PRESR (predicting STXBP1-related disorder), a machine learning algorithm that uses both sequence- and 3-dimensional structure-based features, was developed to improve pathogenicity prediction using 231 known disease-associated variants and comparison to our experimental data. RESULTS: Disease-associated variants, but none of the neutral variants, produced reduced protein levels. Cell-free assays demonstrated directly that disease-associated variants have reduced thermostability, with most variants denaturing around body temperature. In addition, most disease-associated variants impaired SNARE-mediated membrane fusion in a reconstituted assay. Aggregation/insolubility was observed for none of the variants in vitro or in neurons. PRESR outperformed existing tools substantially: Matthews correlation coefficient = 0.71 versus <0.55. CONCLUSIONS: These data establish intrinsic protein instability as the generalizable, primary cause for STXBP1-related disorders and show that protein-specific ortholog and 3-dimensional information improve disease prediction. PRESR is a publicly available diagnostic tool.

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.

Functional analysis of peripheral and intratumoral neoantigen-specific TCRs identified in a patient with melanoma.

BACKGROUND: Neoantigens derived from somatic mutations correlate with therapeutic responses mediated by treatment with immune checkpoint inhibitors. Neoantigens are therefore highly attractive targets for the development of therapeutic approaches in personalized medicine, although many aspects of their quality and associated immune responses are not yet well understood. In a case study of metastatic malignant melanoma, we aimed to perform an in-depth characterization of neoantigens and respective T-cell responses in the context of immune checkpoint modulation. METHODS: Three neoantigens, which we identified either by immunopeptidomics or in silico prediction, were investigated using binding affinity analyses and structural simulations. We isolated seven T-cell receptors (TCRs) from the patient's immune repertoire recognizing these antigens. TCRs were compared in vitro by multiparametric analyses including functional avidity, multicytokine secretion, and cross-reactivity screenings. A xenograft mouse model served to study in vivo functionality of selected TCRs. We investigated the patient's TCR repertoire in blood and different tumor-related tissues over 3 years using TCR beta deep sequencing. RESULTS: Selected mutated peptide ligands with proven immunogenicity showed similar binding affinities to the human leukocyte antigen complex and comparable disparity to their wild-type counterparts in molecular dynamic simulations. Nevertheless, isolated TCRs recognizing these antigens demonstrated distinct patterns in functionality and frequency. TCRs with lower functional avidity showed at least equal antitumor immune responses in vivo. Moreover, they occurred at high frequencies and particularly demonstrated long-term persistence within tumor tissues, lymph nodes and various blood samples associated with a reduced activation pattern on primary in vitro stimulation. CONCLUSIONS: We performed a so far unique fine characterization of neoantigen-specific T-cell responses revealing defined reactivity patterns of neoantigen-specific TCRs. Our data highlight qualitative differences of these TCRs associated with function and longevity of respective T cells. Such features need to be considered for further optimization of neoantigen targeting including adoptive T-cell therapies using TCR-transgenic T cells.

Exchange catalysis by tapasin exploits conserved and allele-specific features of MHC-I molecules.

The repertoire of peptides presented by major histocompatibility complex class I (MHC-I) molecules on the cell surface is tailored by the ER-resident peptide loading complex (PLC), which contains the exchange catalyst tapasin. Tapasin stabilizes MHC-I molecules and promotes the formation of stable peptide-MHC-I (pMHC-I) complexes that serve as T cell antigens. Exchange of suboptimal by high-affinity ligands is catalyzed by tapasin, but the underlying mechanism is still elusive. Here we analyze the tapasin-induced changes in MHC-I dynamics, and find the catalyst to exploit two essential features of MHC-I. First, tapasin recognizes a conserved allosteric site underneath the α2-1-helix of MHC-I, 'loosening' the MHC-I F-pocket region that accomodates the C-terminus of the peptide. Second, the scoop loop11-20 of tapasin relies on residue L18 to target the MHC-I F-pocket, enabling peptide exchange. Meanwhile, tapasin residue K16 plays an accessory role in catalysis of MHC-I allotypes bearing an acidic F-pocket. Thus, our results provide an explanation for the observed allele-specificity of catalyzed peptide exchange.

Major histocompatibility complex (MHC) class I and class II proteins: impact of polymorphism on antigen presentation.

The major histocompatibility complex (MHC) loci are amongst the most polymorphic regions in the genomes of vertebrates. In the human population, thousands of MHC gene variants (alleles) exist that translate into distinct allotypes equipped with overlapping but unique peptide binding profiles. Understanding the differential structural and dynamic properties of MHC alleles and their interaction with critical regulators of peptide exchange bears the potential for more personalized strategies of immune modulation in the context of HLA-associated diseases.

Distinct editing functions of natural HLA-DM allotypes impact antigen presentation and CD4+ T cell activation.

Classical human leukocyte antigen (HLA) molecules of the major histocompatibility class II (MHCII) complex present peptides for the development, surveillance and activation of CD4+ T cells. The nonclassical MHCII-like protein HLA-DM (DM) catalyzes the exchange and loading of peptides onto MHCII molecules, thereby shaping MHCII immunopeptidomes. Natural variations of DM in both chains of the protein (DMA and DMB) have been hypothesized to impact peptide presentation, but no evidence for altered function has been reported. Here we define the presence of DM allotypes in human populations covered by the 1000 Genomes Project and probe their activity. The functional properties of several allotypes are investigated and show strong enhancement of peptide-induced T cell activation for a particular combination of DMA and DMB. Biochemical evidence suggests a broader pH activity profile for the new variant relative to that of the most commonly expressed DM allotype. Immunopeptidome analysis indicates that the compartmental activity of the new DM heterodimer extends beyond the late endosome and suggests that the natural variation of DM has profound effects on adaptive immunity when antigens bypass the canonical processing pathway.

Geometrical characterization of T cell receptor binding modes reveals class-specific binding to maximize access to antigen.

Recognition of antigenic peptides bound to major histocompatibility complex (MHC) proteins by αß T cell receptors (TCRs) is a hallmark of T cell mediated immunity. Recent data suggest that variations in TCR binding geometry may influence T cell signaling, which could help explain outliers in relationships between physical parameters such as TCR-pMHC binding affinity and T cell function. Traditionally, TCR binding geometry has been described with simple descriptors such as the crossing angle, which quantifies what has become known as the TCR's diagonal binding mode. However, these descriptors often fail to reveal distinctions in binding geometry that are apparent through visual inspection. To provide a better framework for relating TCR structure to T cell function, we developed a comprehensive system for quantifying the geometries of how TCRs bind peptide/MHC complexes. We show that our system can discern differences not clearly revealed by more common methods. As an example of its potential to impact biology, we used it to reveal differences in how TCRs bind class I and class II peptide/MHC complexes, which we show allow the TCR to maximize access to and "read out" the peptide antigen. We anticipate our system will be of use in not only exploring these and other details of TCR-peptide/MHC binding interactions, but also addressing questions about how TCR binding geometry relates to T cell function, as well as modeling structural properties of class I and class II TCR-peptide/MHC complexes from sequence information. The system is available at https://tcr3d.ibbr.umd.edu/tcr_com or for download as a script.

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 ß2-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.

FBP21's C-Terminal Domain Remains Dynamic When Wrapped around the c-Sec63 Unit of Brr2 Helicase.

Based on our recent finding that FBP21 regulates human Brr2 helicase activity involved in the activation of the spliceosomal B-complex, we investigated the structural and dynamic contribution of FBP21 to the interaction. By using NMR spectroscopy, we could show that the 50 C-terminal residues of FBP21 (FBP21326-376), which are sufficient to fully form the interaction with the C-terminal Sec63 unit of Brr2 (Brr2C-Sec63), adopt a random-coil conformation in their unbound state. Upon interaction with Brr2C-Sec63, 42 residues of FBP21326-376 cover the large binding site on Brr2C-Sec63 in an extended conformation. Short charged motifs are steering complex formation, still allowing the bound state to retain dynamics. Based on fragment docking in combination with experimental restraints, we present models of the complex structure. The FBP21326-376/Brr2C-Sec63 interaction thus presents an example of an intrinsically disordered protein/ordered-protein interaction in which a large binding site provides high specificity and, in combination with conformational disorder, displays a relatively high affinity.

Quantification of HLA-DM-Dependent Major Histocompatibility Complex of Class II Immunopeptidomes by the Peptide Landscape Antigenic Epitope Alignment Utility.

The major histocompatibility complex of class II (MHCII) immunopeptidome represents the repertoire of antigenic peptides with the potential to activate CD4+ T cells. An understanding of how the relative abundance of specific antigenic epitopes affects the outcome of T cell responses is an important aspect of adaptive immunity and offers a venue to more rationally tailor T cell activation in the context of disease. Recent advances in mass spectrometric instrumentation, computational power, labeling strategies, and software analysis have enabled an increasing number of stratified studies on HLA ligandomes, in the context of both basic and translational research. A key challenge in the case of MHCII immunopeptidomes, often determined for different samples at distinct conditions, is to derive quantitative information on consensus epitopes from antigenic peptides of variable lengths. Here, we present the design and benchmarking of a new algorithm [peptide landscape antigenic epitope alignment utility (PLAtEAU)] allowing the identification and label-free quantification (LFQ) of shared consensus epitopes arising from series of nested peptides. The algorithm simplifies the complexity of the dataset while allowing the identification of nested peptides within relatively short segments of protein sequences. Moreover, we apply this algorithm to the comparison of the ligandomes of cell lines with two different expression levels of the peptide-exchange catalyst HLA-DM. Direct comparison of LFQ intensities determined at the peptide level is inconclusive, as most of the peptides are not significantly enriched due to poor sampling. Applying the PLAtEAU algorithm for grouping of the peptides into consensus epitopes shows that more than half of the total number of epitopes is preferentially and significantly enriched for each condition. This simplification and deconvolution of the complex and ambiguous peptide-level dataset highlights the value of the PLAtEAU algorithm in facilitating robust and accessible quantitative analysis of immunopeptidomes across cellular contexts. In silico analysis of the peptides enriched for each HLA-DM expression conditions suggests a higher affinity of the pool of peptides isolated from the high DM expression samples. Interestingly, our analysis reveals that while for certain autoimmune-relevant epitopes their presentation increases upon DM expression others are clearly edited out from the peptidome.

Author Correction: Protein-peptide association kinetics beyond the seconds timescale from atomistic simulations.

In the original version of this Article, the Acknowledgement section omitted financial support from the Deutsche Forschungsgemeinschaft grant SFB 958/A4. This error has now been corrected in both the PDF and HTML versions of the Article.

Protein-peptide association kinetics beyond the seconds timescale from atomistic simulations.

Understanding and control of structures and rates involved in protein ligand binding are essential for drug design. Unfortunately, atomistic molecular dynamics (MD) simulations cannot directly sample the excessively long residence and rearrangement times of tightly binding complexes. Here we exploit the recently developed multi-ensemble Markov model framework to compute full protein-peptide kinetics of the oncoprotein fragment 25-109Mdm2 and the nano-molar inhibitor peptide PMI. Using this system, we report, for the first time, direct estimates of kinetics beyond the seconds timescale using simulations of an all-atom MD model, with high accuracy and precision. These results only require explicit simulations on the sub-milliseconds timescale and are tested against existing mutagenesis data and our own experimental measurements of the dissociation and association rates. The full kinetic model reveals an overall downhill but rugged binding funnel with multiple pathways. The overall strong binding arises from a variety of conformations with different hydrophobic contact surfaces that interconvert on the milliseconds timescale.

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.

Release from endoplasmic reticulum matrix proteins controls cell surface transport of MHC class I molecules.

The anterograde transport of secretory proteins from the endoplasmic reticulum (ER) to the plasma membrane is a multi-step process. Secretory proteins differ greatly in their transport rates to the cell surface, but the contribution of each individual step to this difference is poorly understood. Transport rates may be determined by protein folding, chaperone association in the ER, access to ER exit sites (ERES) and retrieval from the ER-Golgi intermediate compartment or the cis-Golgi to the ER. We have used a combination of folding and trafficking assays to identify the differential step in the cell surface transport of two natural allotypes of the murine major histocompatibility complex (MHC) class I peptide receptor, H-2D(b) and H-2K(b) . We find that a novel pre-ER exit process that acts on the folded lumenal part of MHC class I molecules and that drastically limits their access to ERES accounts for the transport difference of the two allotypes. Our observations support a model in which the cell surface transport of MHC class I molecules and other type I transmembrane proteins is governed by the affinity of all their folding and maturation states to the proteins of the ER matrix.

F pocket flexibility influences the tapasin dependence of two differentially disease-associated MHC Class I proteins.

The human MHC class I protein HLA-B*27:05 is statistically associated with ankylosing spondylitis, unlike HLA-B*27:09, which differs in a single amino acid in the F pocket of the peptide-binding groove. To understand how this unique amino acid difference leads to a different behavior of the proteins in the cell, we have investigated the conformational stability of both proteins using a combination of in silico and experimental approaches. Here, we show that the binding site of B*27:05 is conformationally disordered in the absence of peptide due to a charge repulsion at the bottom of the F pocket. In agreement with this, B*27:05 requires the chaperone protein tapasin to a greater extent than the conformationally stable B*27:09 in order to remain structured and to bind peptide. Taken together, our data demonstrate a method to predict tapasin dependence and physiological behavior from the sequence and crystal structure of a particular class I allotype. Also watch the Video Abstract.

Proline substitution independently enhances H-2D(b) complex stabilization and TCR recognition of melanoma-associated peptides.

The immunogenicity of H-2D(b) (D(b)) restricted epitopes can be significantly increased by substituting peptide position 3 to a proline (p3P). The p3P modification enhances MHC stability without altering the conformation of the modified epitope allowing for T-cell cross-reactivity with the native peptide. The present study reveals how specific interactions between p3P and the highly conserved MHC heavy chain residue Y159 increase the stability of D(b) in complex with an optimized version of the melanoma-associated epitope gp10025-33 . Furthermore, the p3P modification directly increased the affinity of the D(b)/gp10025-33 -specific T-cell receptor (TCR) pMel. Surprisingly, the enhanced TCR binding was independent from the observed increased stability of the optimized D(b)/gp10025-33 complex and from the interactions formed between p3P and Y159, indicating a direct effect of the p3P modification on TCR recognition.