The class II MHC protein HLA-DR1 in complex with an endogenous peptide: implications for the structural basis of the specificity of peptide binding

In-silico development of a novel TLR2-mediating multi-epitope vaccine against Mycobacterium tuberculosis

Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb), still remains one of the leading causes of mortality worldwide. The elusive nature of this pathogen and its ability to develop drug resistance makes it a serious threat to global health. BCG, the only preventive vaccine for TB, has a limited efficacy and provides partial protection against the disease. A new effective recombinant vaccine capable of producing a stronger and more comprehensive immune response is required to address this global threat. In the present study, we adopted an in-silico approach to develop a multi-epitope vaccine by screening 198 "regulatory proteins" of Mtb H37Rv strain. Epitopes generated from these proteins were screened on the basis of antigenicity, cytokine profile, allergenicity, toxicity, conservancy and population coverage. Selected epitopes were docked with predominant MHC alleles that were used to develop a vaccine construct using suitable linkers and adjuvant. The construct was subjected to homology modelling, tertiary structure validation and refinement and was eventually docked with Toll-like receptor 2 receptor. Molecular dynamic simulation studies revealed stable interactions between the vaccine construct and TLR-2 receptor. The construct also displayed a high probability to elicit a protective immune response involving both humoral and cell-mediated components. In conclusion, the findings suggest that the constructed vaccine has the potential to induce a robust immune response against Mtb. However, further in-vitro and in-vivo studies are required to assess the safety, efficacy, and long-term protective effects of the vaccine construct.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-025-00322-8.

Structural insights into human MHC-II association with invariant chain

The loading of processed peptides on to major histocompatibility complex II (MHC-II) molecules for recognition by T cells is vital to cell-mediated adaptive immunity. As part of this process, MHC-II associates with the invariant chain (Ii) during biosynthesis in the endoplasmic reticulum to prevent premature peptide loading and to serve as a scaffold for subsequent proteolytic processing into MHC-II-CLIP. Cryo-electron microscopy structures of full-length Human Leukocyte Antigen-DR (HLA-DR) and HLA-DQ complexes associated with Ii, resolved at 3.0 to 3.1 Å, elucidate the trimeric assembly of the HLA/Ii complex and define atomic-level interactions between HLA, Ii transmembrane domains, loop domains, and class II-associated invariant chain peptides (CLIP). Together with previous structures of MHC-II peptide loading intermediates DO and DM, our findings complete the structural path governing class II antigen presentation.

A new chromosome-scale duck genome shows a major histocompatibility complex with several expanded multigene families

BACKGROUND

The duck (Anas platyrhynchos) is one of the principal natural hosts of influenza A virus (IAV), harbors almost all subtypes of IAVs and resists to many IAVs which cause extreme virulence in chicken and human. However, the response of duck's adaptive immune system to IAV infection is poorly characterized due to lack of a detailed gene map of the major histocompatibility complex (MHC).

RESULTS

We herein reported a chromosome-scale Beijing duck assembly by integrating Nanopore, Bionano, and Hi-C data. This new reference genome SKLA1.0 covers 40 chromosomes, improves the contig N50 of the previous duck assembly with highest contiguity (ZJU1.0) of more than a 5.79-fold, surpasses the chicken and zebra finch references in sequence contiguity and contains a complete genomic map of the MHC. Our 3D MHC genomic map demonstrated that gene family arrangement in this region was primordial; however, families such as AnplMHCI, AnplMHCIIβ, AnplDMB, NKRL (NK cell receptor-like genes) and BTN underwent gene expansion events making this area complex. These gene families are distributed in two TADs and genes sharing the same TAD may work in a co-regulated model.

CONCLUSIONS

These observations supported the hypothesis that duck's adaptive immunity had been optimized with expanded and diversified key immune genes which might help duck to combat influenza virus. This work provided a high-quality Beijing duck genome for biological research and shed light on new strategies for AIV control.

PANDORA v2.0: Benchmarking peptide-MHC II models and software improvements

T-cell specificity to differentiate between self and non-self relies on T-cell receptor (TCR) recognition of peptides presented by the Major Histocompatibility Complex (MHC). Investigations into the three-dimensional (3D) structures of peptide:MHC (pMHC) complexes have provided valuable insights of MHC functions. Given the limited availability of experimental pMHC structures and considerable diversity of peptides and MHC alleles, it calls for the development of efficient and reliable computational approaches for modeling pMHC structures. Here we present an update of PANDORA and the systematic evaluation of its performance in modelling 3D structures of pMHC class II complexes (pMHC-II), which play a key role in the cancer immune response. PANDORA is a modelling software that can build low-energy models in a few minutes by restraining peptide residues inside the MHC-II binding groove. We benchmarked PANDORA on 136 experimentally determined pMHC-II structures covering 44 unique αβ chain pairs. Our pipeline achieves a median backbone Ligand-Root Mean Squared Deviation (L-RMSD) of 0.42 Å on the binding core and 0.88 Å on the whole peptide for the benchmark dataset. We incorporated software improvements to make PANDORA a pan-allele framework and improved the user interface and software quality. Its computational efficiency allows enriching the wealth of pMHC binding affinity and mass spectrometry data with 3D models. These models can be used as a starting point for molecular dynamics simulations or structure-boosted deep learning algorithms to identify MHC-binding peptides. PANDORA is available as a Python package through Conda or as a source installation at https://github.com/X-lab-3D/PANDORA.

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.

Ancient duplication, coevolution, and selection at the MHC class IIA and IIB genes of birds

Introduction

The Major Histocompatibility Complex (MHC) of vertebrates is a dynamically evolving multigene family primarily responsible for recognizing non-self peptide antigens and triggering a pathogen-specific adaptive immune response. In birds, the MHC was previously thought to evolve via concerted evolution with high degree of gene homogenization and the rapid loss of orthology. However, the discovery of two ancient avian MHC-IIB gene lineages (DAB1 and DAB2) originating before the radiation of extant birds indicated that despite the action of concerted evolution, orthology may be detectable for long evolutionary periods.

Methods

Here, we take advantage of rapidly accumulating digital genomic resources to search for the signal of an ancient duplication at the avian MHC-IIA genes, as well as to compare phylogenetic distribution and selection between MHC-IIA and IIB gene lineages.

Results

The analysis of MHC sequences from over 230 species representing ca. 70 bird families revealed the presence of two ancient MHC-IIA gene lineages (DAA1 and DAA2) and showed that their phylogenetic distribution matches exactly the distribution of DAB1 and DAB2 lineages, suggesting tight coevolution. The early post-duplication divergence of DAA1 and DAA2 was driven by positive selection fixing radical amino acid differences within the membrane-proximal domain and, most probably, being functionally related to the interactions between α2 and β2 chains of the MHC-II heterodimer. We detected no evidence for an overall (gene-wide) relaxation or intensification of selection at either DAA1/DAB1 or DAA2/DAB2, but codon-specific differences in selection signature were found at the peptide-binding sites between the two gene lineages, perhaps implying specialization to different pathogen regimes.

Discussion

Our results suggest that specific pairing of MHC-II α and β chains may have an adaptive significance, a conclusion that advances knowledge on the macroevolution of the avian MHC-II and opens exciting novel directions for future research.

Reverse vaccinology and immunoinformatics approaches to design multi-epitope based vaccine against oncogenic KRAS

Mutant KRAS-induced tumorigenesis is highly involved in the progression of pancreatic, lung, and breast cancer. Comparatively, KRAS G12D and KRAS G12C are the most frequent mutations that promote cancer progression and aggressiveness. Although KRAS mutant inhibitors exhibit significant therapeutic potential, day by day, they are becoming resistant among patients. Multi-epitope based cancer vaccines are a promising alternative strategy that induces an immune response against tumor antigens. In the present study, we have designed, constructed, and validated a novel multi-epitope vaccine construct against KRAS G12D and G12C mutants using reverse vaccinology and immunoinformatics approaches. In addition, the vaccine construct was structurally refined and showed significant physiochemical properties, and could induce an immune response. Furthermore, the optimized vaccine construct was cloned into a pET‑28a (+) expression vector through in silico cloning. Conclusively, the multi-epitope vaccine construct is structurally stable, soluble, antigenic, non‑allergic, and non‑toxic. Further, it has to be studied in in vitro and in vivo to evaluate its therapeutic efficacy against KRAS-mutated cancers in the near future.

HLA-II immunopeptidome profiling and deep learning reveal features of antigenicity to inform antigen discovery

CD4+ T cell responses are exquisitely antigen specific and directed toward peptide epitopes displayed by human leukocyte antigen class II (HLA-II) on antigen-presenting cells. Underrepresentation of diverse alleles in ligand databases and an incomplete understanding of factors affecting antigen presentation in vivo have limited progress in defining principles of peptide immunogenicity. Here, we employed monoallelic immunopeptidomics to identify 358,024 HLA-II binders, with a particular focus on HLA-DQ and HLA-DP. We uncovered peptide-binding patterns across a spectrum of binding affinities and enrichment of structural antigen features. These aspects underpinned the development of context-aware predictor of T cell antigens (CAPTAn), a deep learning model that predicts peptide antigens based on their affinity to HLA-II and full sequence of their source proteins. CAPTAn was instrumental in discovering prevalent T cell epitopes from bacteria in the human microbiome and a pan-variant epitope from SARS-CoV-2. Together CAPTAn and associated datasets present a resource for antigen discovery and the unraveling genetic associations of HLA alleles with immunopathologies.

In-silico design and evaluation of an epitope-based serotype-independent promising vaccine candidate for highly cross-reactive regions of pneumococcal surface protein A

BACKGROUND

The pathogenicity of pneumococcus with high morbidity, mortality, and multi-drug resistance patterns has been increasing. The limited coverage of the licensed polysaccharide-based vaccines and the replacement of the non-vaccine serotypes are the main reasons for producing a successful serotype-independent vaccine. Pneumococcal surface protein A (PspA) is an extremely important virulence factor and an interesting candidate for conserved protein-based pneumococcal vaccine classified into two prominent families containing five clades. PspA family-elicited immunity is clade-dependent, and the level of the PspA cross-reactivity is restricted to the same family.

METHODS

To cover and overcome the clade-dependent immunity of the PspAs in this study, we designed and tested a PspA_1-5c+p vaccine candidate composed of the highest immunodominant coverage of B- and T-cell epitope truncated domain of each clade focusing on two cross-reactive B and C regions of the PspAs. The antigenicity, toxicity, physicochemical properties, 3D structure prediction, stability and flexibility of the designed protein using molecular dynamic (MD) simulation, molecular docking of the construct withHLADRB1*(01:01) and human lactoferrin N-lop, and immune simulation were assessed using immunoinformatics tools. In the experimental section, after intraperitoneal immunization of the mice with Alum adjuvanted recombinant PspA_1-5c+p, we evaluated the immune response, cross-reactivity, and functionality of the Anti-PspA_1-5c+p antibody using ELISA, Opsonophagocytic killing activity, and serum bactericidal assay.

RESULTS

For the first time, this work suggested a novel PspA-based vaccine candidate using immunoinformatics tools. The designed PspA_1-5c+p protein is predicted to be highly antigenic, non-toxic, soluble, stable with low flexibility in MD simulation, and able to stimulate both humoral and cellular immune responses. The designed protein also could interact strongly with HLADRB1*(01:01) and human lactoferrin N-lop in the docking study. Our immunoinformatics predictions were validated using experimental data. Results showed that the anti-PspA_1-5c+p IgG not only had a high titer with strong and same cross-reactivity coverage against all pneumococcal serotypes used but also had high and effective bioactivity for pneumococcal clearance using complement system and phagocytic cells.

CONCLUSION

Our findings elucidated the potential application of the PspA_1-5c+p vaccine candidate as a serotype-independent pneumococcal vaccine with a strong cross-reactivity feature. Further in-vitro and in-vivo investigations against other PspA clades should be performed to confirm the full protection of the PspA_1-5c+p vaccine candidate.

HLA Class I Supertype Classification Based on Structural Similarity

HLA class I proteins, a critical component in adaptive immunity, bind and present intracellular Ags to CD8+ T cells. The extreme polymorphism of HLA genes and associated peptide binding specificities leads to challenges in various endeavors, including neoantigen vaccine development, disease association studies, and HLA typing. Supertype classification, defined by clustering functionally similar HLA alleles, has proven helpful in reducing the complexity of distinguishing alleles. However, determining supertypes via experiments is impractical, and current in silico classification methods exhibit limitations in stability and functional relevance. In this study, by incorporating three-dimensional structures we present a method for classifying HLA class I molecules with improved breadth, accuracy, stability, and flexibility. Critical for these advances is our finding that structural similarity highly correlates with peptide binding specificity. The new classification should be broadly useful in peptide-based vaccine development and HLA-disease association studies.

An immunoinformatic approach to assessing the immunogenic capacity of alpha-neurotoxins in elapid snake venoms

Introduction: Most elapid snakes produce venoms that contain alpha-neurotoxins (α-NTXs), which are proteins that cause post-synaptic blockade and paralysis in snakebite envenoming. However, existing elapid antivenoms are known for their low potency in neutralizing the neurotoxic activity of α-NTXs, while the immunological basis has not been elucidated. Methods: In this study, a structure-based major histocompatibility complex II (MHCII) epitope predictor of horse (Equus caballus), complemented with DM-editing determinant screening algorithm was adopted to assess the immunogenicity of α-NTXs in the venoms of major Asiatic elapids (Naja kaouthia, Ophiophagus hannah, Laticauda colubrina, Hydrophis schistosus, Hydrophis curtus). Results: The scoring metric M₂R, representing the relative immunogenic performance of respective α-NTXs, showed all α-NTXs have an overall low M₂R of <0.3, and most of the predicted binders feature non-optimal P1 anchor residues. The M₂R scores correlate strongly (R² = 0.82) with the potency scores (p-score) generated based on the relative abundances of α-NTXs and the neutralization potency of commercial antivenoms. Discussion: The immunoinformatic analysis indicates that the inferior antigenicity of α-NTXs is not only due to their small molecular size but also the subpar immunogenicity affected by their amino acid composition. Structural modification with conjugation and synthetic epitope as immunogen may potentially enhance the immunogenicity for improved antivenom potency against α-NTXs of elapid snakes.

Stereo electronic principles for selecting fully-protective, chemically-synthesised malaria vaccines

Major histocompatibility class II molecule-peptide-T-cell receptor (MHCII-p-TCR) complex-mediated antigen presentation for a minimal subunit-based, multi-epitope, multistage, chemically-synthesised antimalarial vaccine is essential for inducing an appropriate immune response. Deep understanding of this MHCII-p-TCR complex’s stereo-electronic characteristics is fundamental for vaccine development. This review encapsulates the main principles for achieving such epitopes’ perfect fit into MHC-II human (HLADRβ̞1*) or Aotus (Aona DR) molecules. The enormous relevance of several amino acids’ physico-chemical characteristics is analysed in-depth, as is data regarding a 26.5 ± 2.5Å distance between the farthest atoms fitting into HLA-DRβ1* structures’ Pockets 1 to 9, the role of polyproline II-like (PPIIL) structures having their O and N backbone atoms orientated for establishing H-bonds with specific HLA-DRβ1*-peptide binding region (PBR) residues. The importance of residues having specific charge and orientation towards the TCR for inducing appropriate immune activation, amino acids’ role and that of structures interfering with PPIIL formation and other principles are demonstrated which have to be taken into account when designing immune, protection-inducing peptide structures (IMPIPS) against diseases scourging humankind, malaria being one of them.

Immunoinformatics-Aided Analysis of RSV Fusion and Attachment Glycoproteins to Design a Potent Multi-Epitope Vaccine

Respiratory syncytial virus (RSV) usually causes respiratory tract infections of upper airways in infants and young children. Despite recent medical advances, no approved vaccine is available to control RSV infections. Therefore, we conducted an immunoinformatics study to design and evaluate a potential multi-epitope vaccine against RSV. Sequence-based analyses of the glycoproteins F and G revealed a total of eight CD8 T-cell and three CD4 T-cell epitopes after considering antigenicity, binding affinity and other parameters. Molecular docking analysis confirmed that these T-cell epitopes developed strong structural associations with HLA allele(s). By integrating these prioritized epitopes with linkers and a cholera toxin-derived adjuvant, a multi-epitope vaccine was designed. The developed vaccine was found to be stable, non-allergenic, flexible and antigenic. Molecular docking analysis revealed a striking mean HADDOCK score (-143.3) of top-ranked vaccine-TLR cluster and a Gibbs free energy change (ΔG) value of -11.3 kcal mol-1. As per computational immune simulation results, the vaccine generated a high titer of antibodies (especially IgM) and effector T-cells. Also, codon optimization and in silico cloning ensured the increased expression of vaccine in Escherichia coli. Altogether, we anticipate that the multi-epitope vaccine reported in this study will stimulate humoral and cellular responses against RSV infection, subject to follow-up experimental validation.

Discovery of an ancient MHC category with both class I and class II features

Two classes of major histocompatibility complex (MHC) molecules, MHC class I and MHC class II, constitute the basis of our elaborate, adaptive immune system as antigen-presenting molecules. They perform distinct, critical functions: especially, MHC class I in case of antivirus and antitumor defenses, and MHC class II, in case of effective antibody responses. This important class diversification has long been enigmatic, as vestiges of the evolutionary molecular changes have not been found. The revealed ancient MHC category represents a plausible intermediate group between the two classes, and the data suggest that class II preceded class I in molecular evolution. Fundamental understanding of the molecular evolution of MHC molecules should contribute to understanding the basis of our complex biological defense system.

Two classes of major histocompatibility complex (MHC) molecules, MHC class I and class II, play important roles in our immune system, presenting antigens to functionally distinct T lymphocyte populations. However, the origin of this essential MHC class divergence is poorly understood. Here, we discovered a category of MHC molecules (W-category) in the most primitive jawed vertebrates, cartilaginous fish, and also in bony fish and tetrapods. W-category, surprisingly, possesses class II–type α- and β-chain organization together with class I–specific sequence motifs for interdomain binding, and the W-category α2 domain shows unprecedented, phylogenetic similarity with β₂-microglobulin of class I. Based on the results, we propose a model in which the ancestral MHC class I molecule evolved from class II–type W-category. The discovery of the ancient MHC group, W-category, sheds a light on the long-standing critical question of the MHC class divergence and suggests that class II type came first.

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.

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.

Mining the Mycobacterium tuberculosis proteome for identification of potential T-cell epitope based vaccine candidates

Identification of protective antigens for designing a high-efficacy tuberculosis vaccine is the need of the hour. Till date only 7% of the Mycobacterium tuberculosis proteome has been explored for discovering antigens capable of activating T-cell responses. Therefore, it becomes crucial to screen the remaining Mycobacterium tuberculosis proteome for more immunodominant T-cell epitopes. An extensive knowledge of the epitopes recognized by our immune system can aid this process of finding potential T cell antigens for development of a better TB vaccine. In the present in-silico study, 237 proteins belonging to the 'virulence, detoxification, and adaptation' category of Mycobacterium tuberculosis proteome were targeted for T-cell epitope screening. 50825 MHC Class I and 49357 MHC Class II epitopes were generated using NetMHC3.4 and IEDB servers respectively and tested for their antigenicity and cytokine stimulation. The highest antigenic epitopes were analyzed for their world population coverage and epitope conservancy. Molecular docking and molecular dynamics simulation studies were performed to corroborate the binding affinities and structural stability of the peptide-MHC complexes. We predicted a total of 3 MHC Class I (ILLKMCWPA, FAVGMNVYV, and SLAGNSAKV) and 7 MHC Class II (DLTIGFFLHIPFPPV, RPDLTIGFFLHIPFP, LTIGFFLHIPFPPVE, VLVFALVVALVYLQF, LVFALVVALVYLQFR, PNLVAARFIQLTPVY, and LVLVFALVVALVYLQ) epitopes that can be promising vaccine candidates. These predicted epitopes belong to 6 distinct proteins: Rv0169 (mce1a), Rv3490 (ostA), Rv3496 (mce4D), Rv1085c, Rv0563 (HtpX), Rv3497c (mce4C). All these proteins are expressed at different stages in the life cycle of Mycobacterium tuberculosis and thus, the predicted epitopes could be employed as candidates for designing a multistage-multiepitopic vaccine.

Design of an epitope-based peptide vaccine against the SARS-CoV-2: a vaccine-informatics approach

The recurrent and recent global outbreak of the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has turned into a global concern which has infected more than 42 million people all over the globe, and this number is increasing in hours. Unfortunately, no vaccine or specific treatment is available, which makes it more deadly. A vaccine-informatics approach has shown significant breakthrough in peptide-based epitope mapping and opens the new horizon in vaccine development. In this study, we have identified a total of 15 antigenic peptides [including thymus cells (T-cells) and bone marrow or bursa-derived cells] in the surface glycoprotein (SG) of SARS-CoV-2 which is nontoxic and nonallergenic in nature, nonallergenic, highly antigenic and non-mutated in other SARS-CoV-2 virus strains. The population coverage analysis has found that cluster of differentiation 4 (CD4+) T-cell peptides showed higher cumulative population coverage over cluster of differentiation 8 (CD8+) peptides in the 16 different geographical regions of the world. We identified 12 peptides ((LTDEMIAQY, WTAGAAAYY, WMESEFRVY, IRASANLAA, FGAISSVLN, VKQLSSNFG, FAMQMAYRF, FGAGAALQI, YGFQPTNGVGYQ, LPDPSKPSKR, QTQTNSPRRARS and VITPGTNTSN) that are $80\hbox{--} 90\%$ identical with experimentally determined epitopes of SARS-CoV, and this will likely be beneficial for a quick progression of the vaccine design. Moreover, docking analysis suggested that the identified peptides are tightly bound in the groove of human leukocyte antigen molecules which can induce the T-cell response. Overall, this study allows us to determine potent peptide antigen targets in the SG on intuitive grounds, which opens up a new horizon in the coronavirus disease (COVID-19) research. However, this study needs experimental validation by in vitro and in vivo.

Conserved HLA binding peptides from five non-structural proteins of SARS-CoV-2-An in silico glance

Coronavirus Disease 2019 (COVID-19) is a dangerous global threat that has no clinically approved treatment yet. Bioinformatics represent an outstanding approach to reveal key immunogenic regions in viral proteins. Here, five severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) non-structural proteins (NSPs) (NSP7, NSP8, NSP9, NSP12, and NSP13) were screened to identify potential human leukocyte antigen (HLA) binding peptides. These peptides showed robust viral antigenicity, immunogenicity, and a marked interaction with HLA alleles. Interestingly, several peptides showed affinity by HLA class I (HLA-I) alleles that commonly activates to natural killer (NK) cells. Notably, HLA biding peptides are conserved among SARS-CoV-2, severe acute respiratory syndrome coronavirus (SARS-CoV), and Middle Eastern respiratory syndrome coronavirus (MERS-CoV). Interestingly, HLA-I and HLA class II (HLA-II) binding peptides induced humoral and cell-mediated responses after in silico vaccination. These results may open further in vitro and in vivo investigations to develop novel therapeutic strategies against coronaviral infections.

HLA loci predisposing to immune TTP in Japanese: potential role of the shared ADAMTS13 peptide bound to different HLA-DR

Immune-mediated thrombotic thrombocytopenic purpura (iTTP) is a rare autoimmune disorder caused by neutralizing anti-ADAMTS13 autoantibodies. In white individuals, HLA allele DRB1*11 is a predisposing factor for iTTP, whereas DRB1*04 is a protective factor. However, the role of HLA in Asians is unclear. In this study, we analyzed 10 HLA loci using next-generation sequencing in 52 Japanese patients with iTTP, and the allele frequency in the iTTP group was compared with that in a Japanese control group. We identified the following HLA alleles as predisposing factors for iTTP in the Japanese population: DRB1*08:03 (odds ratio [OR], 3.06; corrected P [Pc] = .005), DRB3/4/5*blank (OR, 2.3; Pc = .007), DQA1*01:03 (OR, 2.25; Pc = .006), and DQB1*06:01 (OR,: 2.41; Pc = .003). The estimated haplotype consisting of these 4 alleles was significantly more frequent in the iTTP group than in the control group (30.8% vs 6.0%; Pc &lt; .001). DRB1*15:01 and DRB5*01:01 were weak protective factors for iTTP (OR, 0.23; Pc = .076; and OR, 0.23, Pc = .034, respectively). On the other hand, DRB1*11 and DRB1*04 were not associated with iTTP in the Japanese. These findings indicated that predisposing and protective factors for iTTP differ between Japanese and white individuals. HLA-DR molecules encoded by DRB1*08:03 and DRB1*11:01 have different peptide-binding motifs, but interestingly, bound to the shared ADAMTS13 peptide in an in silico prediction model.

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.

Exploring T & B-cell epitopes and designing multi-epitope subunit vaccine targeting integration step of HIV-1 lifecycle using immunoinformatics approach

Till now, AIDS, caused by the human immunodeficiency virus (HIV) is still a severe health problem worldwide. It weakens the immune system by targeting the T-helper cells. Specifically, the severity of the pandemic HIV-1 makes the emergence of an enduring effective vaccine against HIV-1. Therefore, we have applied a series of immunoinformatics approaches to the four conserved domains of HIV-1 integrase (IN) proteins to design an effective multi-epitope based subunit vaccine which might induce a competent immunity against HIV-1. Therefore, we have selected three peptide fragments that contained all overlapping epitopes (35 CD4⁺, 8 CD8⁺ T-cell epitopes, and 3 B-cell epitopes) where the epitopes had a high conservancy score. The cumulative population coverage for combined CD8⁺ and CD4⁺ T-cell epitopes and their respective HLA-alleles were found as 98.03% in the world which is also followed by East Asia (96.24%), South Asia (96.31%), North Africa (96.14%), North America (98.99%), and Europe (98.80%). The proposed vaccine composed by an adjuvant (β-defensin) at the N-terminal site of the vaccine constructs and three peptide fragments where the adjuvant was fused by EAAAK linker and the peptide fragments were fused by GPGPG linker. The designed final vaccine construct (length: 159 amino acid) was found to be antigenic and non-allergic, which indicates its safety. The vaccine construct was found as good antigenic, stable, higher thermostable, and hydrophilic in nature. The codon adaptation and in silico cloning ensured the high expression rate of the vaccine constructs in E. coli K12 with CAI value of 1.0. Finally, the binding affinity of the vaccine constructs with the immune receptor TLR3 was confirmed by the lowest energy score of -1026.8 evaluated by molecular docking. However, the proposed in silico vaccine construct needs experimental validation for assuring the safety and immunogenicity profile which will ensure an active immunity against HIV-1.

Ligandomes obtained from different HLA-class II-molecules are homologous for N- and C-terminal residues outside the peptide-binding cleft

Human CD4+ T lymphocytes play an important role in inducing potent immune responses. T cells are activated and stimulated by peptides presented in human leucocyte antigen (HLA)-class II molecules. These HLA-class II molecules typically present peptides of between 12 and 20 amino acids in length. The region that interacts with the HLA molecule, designated as the peptide-binding core, is highly conserved in the residues which anchor the peptide to the molecule. In addition, as these peptides are the product of proteolytic cleavages, certain conserved residues may be expected at the N- and C-termini outside the binding core. To study whether similar conserved residues are present in different cell types, potentially harbouring different proteolytic enzymes, the ligandomes of HLA-DRB1*03:01/HLA-DRB > 1 derived from two different cell types (dendritic cells and EBV-transformed B cells) were identified with mass spectrometry and the binding core and N- and C-terminal residues of a total of 16,568 peptides were analysed using the frequencies of the amino acids in the human proteome. Similar binding motifs were found as well as comparable conservations in the N- and C-terminal residues. Furthermore, the terminal conservations of these ligandomes were compared to the N- and C-terminal conservations of the ligandome acquired from dendritic cells homozygous for HLA-DRB1*04:01. Again, comparable conservations were evident with only minor differences. Taken together, these data show that there are conservations in the terminal residues of peptides, presumably the result of the activity of proteases involved in antigen processing.

Molecular and Structural Characterization of MHC Class II β Genes Reveals High Diversity in the Cold-Adapted Icefish Chionodraco hamatus

This study reports the presence of two distinct MHC class II β genes in the Antarctic icefish Chionodraco hamatus, belonging to the classical (ChhaDAB) and nonclassical (ChhaDBB) evolutionary lineages. By the application of targeted sequencing approach, a remarkable molecular diversity in the exon 2 sequence of the highly expressed gene ChhaDAB has been observed, resulting in an estimate of 92 different variants translated in 87 different peptides from 54 analysed icefish individuals. A highly conservative estimate, based on a 95% sequence identity threshold clustering, translate this variability in 41 different peptide clusters belonging to four different clades and showing the signature of different kinds of selection. In stark contrast, the poorly expressed ChhaDBB gene displayed a very low level of molecular diversity within exon 2, in agreement with expectations for a nonclassical MHC class II β gene.

Interactions of HLA-DR and Topoisomerase I Epitope Modulated Genetic Risk for Systemic Sclerosis

The association of systemic sclerosis with anti-Topoisomerase 1 antibody (ATASSc) with specific alleles of human leukocyte antigen (HLA)-DR has been observed among various ethnics. The anti-Topoisomerase 1 antibody is a common autoantibody in SSc with diffuse cutaneous scleroderma, which is one of the clinical subtypes of SSc. On the other hand, an immunodominant peptide of topoisomerase 1 (Top1) self-protein (residues 349-368) was reported to have strong association with ATASSc. In this study, molecular dynamics simulation was performed on the complexes of Top1 peptide with various HLA-DR subtypes divided into ATASSc-associated alleles (HLA-DRB1*08:02, HLA-DRB1*11:01 and HLA-DRB1*11:04), suspected allele (HLA-DRB5*01:02), and non-associated allele (HLA-DRB1*01:01). The unique interaction for each system was compared to the others in terms of dynamical behaviors, binding free energies and solvation effects. Our results showed that three HLA-DR/Top1 complexes of ATASSc association mostly exhibited high protein stability and increased binding efficiency without solvent interruption, in contrast to non-association. The suspected case (HLA-DRB5*01:02) binds Top1 as strongly as the ATASSc association case, which implied a highly possible risk for ATASSc development. This finding might support ATASSc development mechanism leading to a guideline for the treatment and avoidance of pathogens like Top1 self-peptide risk for ATASSc.

Subtractive proteomics and immunoinformatics revealed novel B-cell derived T-cell epitopes against Yersinia enterocolitica: An etiological agent of Yersiniosis

Yersinia enterocolitica is the third most common cause of gastrointestinal manifestations in Europe. Statistically, every year the pathogen accounts for 640 hospitalizations, 117,000 illnesses, and 35 deaths in the United States. The associated mortality rate of the pathogen is 50% and is virtually resistant to penicillin G, ampicillin and cephalotin. The development of new and effective therapeutic procedures is urgently needed to counter the multi-drug-resistant phenotypes imposed by the said pathogen. Based on subtractive reverse vaccinology and immunoinformatics approaches, we have successfully predicted novel antigenic peptide vaccine candidates against Y. enterocolitica. The pipeline revealed two isoforms of ompC family; meoA (ompC) and ompC2 as promising vaccine targets. Protein-protein interactions elaborated the involvement of target candidates in the major biological pathways of the pathogen. The predicted 9-mer B-cell derived T-cell epitope of proteins are found to be virulent, antigenic, non-allergic, surface exposed and conserved in all nine completely sequenced strains of the pathogen. Molecular docking predicts deep and stable binding of the epitopes in the binding pocket of the most predominant allele in human population-the DRB1*0101. These epitopes of target proteins could provide the foundation for the development of an epitope-driven vaccine against Y. enterocolitica.

Proteome-wide identification of epitope-based vaccine candidates against multi-drug resistant Proteus mirabilis

Proteus mirabilis is one of the important pathogens of urinary tract and exhibits resistance to multiple drugs. Development of vaccine tends to be the most promising and cost-effective remedy against the said pathogen. Herein, we implement a combinatorial approach for screening proteins harboring potential broad-spectrum antigenic epitopes in the proteome of P. mirabilis. The targets are host non-homologous, essential and virulent, and have localization in the extracellular and outer membrane. Immuno-informatics revealed antigenic, surface exposed and broad-spectrum B-cell derived T-cell epitopes for three membrane usher family candidates: AtfC, PMI2533 and PMI1466, which could evoke a substantial immune response. Protein-protein interactions of targeted three proteins have shown their involvement in biologically significant pathways indispensable for the growth and survival of the pathogen. The antigenic epitopes are conserved among all completely annotated strains and docked deeply in the binding cavity of the most prevalent allele-DRB1*0101 in human population. Future work is necessary to characterize the shortlisted proteins and epitopes for immune protection in animal models.

The N-terminal region of photocleavable peptides that bind HLA-DR1 determines the kinetics of fragment release

Major Histocompatibility Complex class II (MHC-II) molecules bind peptides and present them to receptors on CD4+ T cells as part of the immune system's surveillance of pathogens and malignancy. In the absence of peptide, MHC-II equilibrates between peptide-receptive and peptide-averse conformations. The conversion between these forms has been postulated to be important in regulating cellular antigen presentation but has been difficult to study. In order to generate the MHC-II molecule HLA-DR1 in the peptide-receptive form, we designed and tested a series of photocleavable peptides that included the UV-sensitive 3-amino-3-(2-nitrophenyl)-propionate amino acid analog. They were intended to bind tightly to the HLA-DR1 MHC molecule, but to generate low-affinity fragments after UV exposure that would be released to yield HLA-DR1 in the peptide-receptive conformation. We were able to identify photocleavable peptides that bound tightly to HLA-DR1 and generated the peptide-receptive conformation after UV exposure. However, slow release of photocleaved peptide fragments from the binding site limited the rate of binding of an incoming labeled peptide and complicated kinetic measurements of the individual steps of the overall peptide binding reaction. Modification of the N-terminal region of the photocleavable peptide to reduce MHC-II pocket or H-bonding interactions allowed for generation of the peptide receptive form immediately after UV exposure with peptide fragments neither retained within the site nor interfering with binding of an incoming peptide. However this was achieved only at the expense of a substantial reduction in overall peptide binding affinity, and these peptides had such weak interaction with HLA-DR1 that they were easily exchanged by incoming peptide without UV exposure. These results show that photocleavable peptides can be used to generate peptide-receptive HLA-DR1 and to facilitate peptide exchange in generation of specific peptide-MHC-II complexes, but that usage of these peptides for kinetic studies can be constrained by slow fragment release.

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

Major histocompatibility complex class two (MHC-II) molecules are trans-membrane proteins and key components of the cellular immune system. Upon recognition of foreign peptides expressed on the MHC-II binding groove, CD4⁺ T cells mount an immune response against invading pathogens. Therefore, mechanistic identification and knowledge of physicochemical features that govern interactions between peptides and MHC-II molecules is useful for the design of effective epitope-based vaccines, as well as for understanding of immune responses. In this article, we present a comprehensive trans-allelic prediction model, a generalized version of our previous biophysical model, that can predict peptide interactions for all three human MHC-II loci (HLA-DR, HLA-DP, and HLA-DQ), using both peptide sequence data and structural information of MHC-II molecules. The advantage of this approach over other machine learning models is that it offers a simple and plausible physical explanation for peptide–MHC-II interactions. We train the model using a benchmark experimental dataset and measure its predictive performance using novel data. Despite its relative simplicity, we find that the model has comparable performance to the state-of-the-art method, the NetMHCIIpan method. Focusing on the physical basis of peptide–MHC binding, we find support for previous theoretical predictions about the contributions of certain binding pockets to the binding energy. In addition, we find that binding pocket P5 of HLA-DP, which was not previously considered as a primary anchor, does make strong contribution to the binding energy. Together, the results indicate that our model can serve as a useful complement to alternative approaches to predicting peptide–MHC interactions.

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.

Binding Patterns Associated Aß-HSP60 p458 Conjugate to HLA-DR-DRB Allele of Human in Alzheimer's Disease: An In Silico Approach

Alzheimer's disease (AD) is a complex, irreversible, progressive brain disorder, which diminishes memory in a slow pace and thinking skills; ranked third by experts. It is a complex disorder that involves numerous cellular and subcellular alterations. The pathogenesis of AD is still unknown, but for better understanding, we proposed an in silico analysis to find out the binding patterns associated with HSP60. Several experimental conclusions have been drawn to understand the actual mechanism behind the forming of aggregation due to misfolding. Protein misfolding disorder is experimentally identified by the accumulation of protein aggregates at the intracellular or extracellular region of brain that adversely affects the cell functioning by disrupting the connection between the cells and ultimately leading to cell death. To unravel the mystery behind the mechanism of AD through computational approach, the current proposal shows the designing of Aß-HSP60 p458 conjugate followed by secondary structure analysis, which is further targeted to HLA-DR-DRB allele of human. The antigenicity of Aß (1-42) peptide is the major concern in our study predicted through PVS server, which provides an insight into the immunogenic behavior of Aß peptide. The mechanism involved in the interaction of HSP60-Aß conjugate with HLA-DR-DRB allele considering the fact that Aß (1-42) is highly immunogenic in human and interactions evoked highly robust T-cell response through MHC class II binding predictions. It was assisted by molecular dynamics simulation of predicted HSP60 structure followed by validation through Ramachandran plot analysis and protein-protein interaction of Aß (1-42) with HSP60.

Granuloma genes in sarcoidosis: what is new?

PURPOSE OF REVIEW

Nonnecrotizing granulomas in the affected organ are the hallmark of sarcoidosis. This review summarizes most recent genetic findings in sarcoidosis with a focus on genes that might influence granuloma formation or resolution. Specific results in multiple ethnic groups and certain clinical subphenotypes, such as extra-pulmonary organ involvement, are discussed.

RECENT FINDINGS

Associations of genetic variants in antigen-presenting molecules (HLA-DRB1) were shown to confer risk to sarcoidosis and certain disease phenotypes in populations of different ethnic origins. Specific DRB1 alleles, such as *0301 and *0302, appear to confer protection against chronic disease, but in an ethnic-specific manner illustrating the extensive genetic heterogeneity and complexity at this locus. Mechanistic studies of putative sarcoid antigens lend further credence to a role of HLA-DRB1 in disease pathogenesis. With relevance to granuloma formation, genes involved in apoptotic processes and immune cell activation were further confirmed (ANXA11 and BTNL2) in multiple ethnicities; others were newly identified (XAF1). Linking mechanism to clinical application, a TNF variant was shown to correlate with anti-TNF response in sarcoidosis patients.

SUMMARY

The investigation of known and novel risk variants for sarcoidosis and specific clinical phenotypes in various ethnicities highlights the genetic complexity of the disease. Detailed subanalysis of disease phenotypes revealed the potential for prediction of extra-pulmonary organ involvement and therapy response based on the patient's genotype.

Matched Peptides: Tuning Matched Molecular Pair Analysis for Biopharmaceutical Applications

Biopharmaceuticals hold great promise for the future of drug discovery. Nevertheless, rational drug design strategies are mainly focused on the discovery of small synthetic molecules. Herein we present matched peptides, an innovative analysis technique for biological data related to peptide and protein sequences. It represents an extension of matched molecular pair analysis toward macromolecular sequence data and allows quantitative predictions of the effect of single amino acid substitutions on the basis of statistical data on known transformations. We demonstrate the application of matched peptides to a data set of major histocompatibility complex class II peptide ligands and discuss the trends captured with respect to classical quantitative structure–activity relationship approaches as well as structural aspects of the investigated protein–peptide interface. We expect our novel readily interpretable tool at the interface of cheminformatics and bioinformatics to support the rational design of biopharmaceuticals and give directions for further development of the presented methodology.

Association of HLA-DRB1 with Sarcoidosis Susceptibility and Progression in African Americans

HLA-DRB1 is a sarcoidosis risk gene, and the *03:01 allele is strongly associated with disease resolution in European sarcoidosis cases. Whereas the HLA-DRB1 variation is associated with sarcoidosis susceptibility in African Americans, DRB1 risk alleles are not as well defined, and associations with disease resolution have not been studied. Associations between genotyped and imputed HLA-DRB1 alleles and disease susceptibility/resolution were evaluated in a sample of 1,277 African-American patients with sarcoidosis and 1,467 control subjects. In silico binding assays were performed to assess the functional significance of the associated alleles. Increased disease susceptibility was associated with the HLA-DRB1 alleles *12:01 (odds ratio [OR], 2.11; 95% confidence interval [CI], 1.65-2.69; P = 3.2 × 10(-9)) and *11:01 (OR, 1.69; 95% CI, 1.42-2.01; P = 3.0 × 10(-9)). The strongest protective association was found with *03:01 (OR, 0.56; 95% CI, 0.44-0.73; P = 1.0 × 10(-5)). The African-derived allele *03:02 was associated with decreased risk of persistent radiographic disease (OR, 0.52; 95% CI, 0.37-0.72; P = 1.3 × 10(-4)), a finding consistent across the three component studies comprising the analytic sample. The DRB1*03:01 association with disease persistence was dependent upon local ancestry, with carriers of at least one European allele at DRB1 at a decreased risk of persistent disease (OR, 0.36; 95% CI, 0.14-0.94; P = 0.037). Results of in silico binding analyses showed that DRB1*03:01 consistently demonstrated the highest binding affinities for six bacterial peptides previously found in sarcoidosis granulomas, whereas *12:01 displayed the lowest binding affinities. This study has identified DRB1*03:01 and *03:02 as novel alleles associated with disease susceptibility and course in African Americans. Further investigation of DRB1*03 alleles may uncover immunologic factors that favor sarcoidosis protection and resolution among African Americans.

Crossreactivity to vinculin and microbes provides a molecular basis for HLA-based protection against rheumatoid arthritis

The HLA locus is the strongest risk factor for anti-citrullinated protein antibody (ACPA)(+) rheumatoid arthritis (RA). Despite considerable efforts in the last 35 years, this association is poorly understood. Here we identify (citrullinated) vinculin, present in the joints of ACPA(+) RA patients, as an autoantigen targeted by ACPA and CD4(+) T cells. These T cells recognize an epitope with the core sequence DERAA, which is also found in many microbes and in protective HLA-DRB1*13 molecules, presented by predisposing HLA-DQ molecules. Moreover, these T cells crossreact with vinculin-derived and microbial-derived DERAA epitopes. Intriguingly, DERAA-directed T cells are not detected in HLA-DRB1*13(+) donors, indicating that the DERAA epitope from HLA-DRB1*13 mediates (thymic) tolerance in these donors and explaining the protective effects associated with HLA-DRB1*13. Together our data indicate the involvement of pathogen-induced DERAA-directed T cells in the HLA-RA association and provide a molecular basis for the contribution of protective/predisposing HLA alleles.

Susceptibility to HLA-DM protein is determined by a dynamic conformation of major histocompatibility complex class II molecule bound with peptide

HLA-DM mediates the exchange of peptides loaded onto MHCII molecules during antigen presentation by a mechanism that remains unclear and controversial. Here, we investigated the sequence and structural determinants of HLA-DM interaction. Peptides interacting nonoptimally in the P1 pocket exhibited low MHCII binding affinity and kinetic instability and were highly susceptible to HLA-DM-mediated peptide exchange. These changes were accompanied by conformational alterations detected by surface plasmon resonance, SDS resistance assay, antibody binding assay, gel filtration, dynamic light scattering, small angle x-ray scattering, and NMR spectroscopy. Surprisingly, all of those changes could be reversed by substitution of the P9 pocket anchor residue. Moreover, MHCII mutations outside the P1 pocket and the HLA-DM interaction site increased HLA-DM susceptibility. These results indicate that a dynamic MHCII conformational determinant rather than P1 pocket occupancy is the key factor determining susceptibility to HLA-DM-mediated peptide exchange and provide a molecular mechanism for HLA-DM to efficiently target unstable MHCII-peptide complexes for editing and exchange those for more stable ones.

HLA-DRα1 constructs block CD74 expression and MIF effects in experimental autoimmune encephalomyelitis

CD74, the cell-surface form of the MHC class II invariant chain, is a key inflammatory factor that is involved in various immune-mediated diseases as part of the macrophage migration inhibitory factor (MIF) binding complex. However, little is known about the natural regulators of CD74 in this context. In order to study the role of the HLA-DR molecule in regulating CD74, we used the HLA-DRα1 domain, which was shown to bind to and downregulate CD74 on CD11b(+) monocytes. We found that DRα1 directly inhibited binding of MIF to CD74 and blocked its downstream inflammatory effects in the spinal cord of mice with experimental autoimmune encephalomyelitis (EAE). Potency of the DRα1 domain could be destroyed by trypsin digestion but enhanced by addition of a peptide extension (myelin oligodendrocyte glycoprotein [MOG]-35-55 peptide) that provided secondary structure not present in DRα1. These data suggest a conformationally sensitive determinant on DRα1-MOG that is responsible for optimal binding to CD74 and antagonism of MIF effects, resulting in reduced axonal damage and reversal of ongoing clinical and histological signs of EAE. These results demonstrate natural antagonist activity of DRα1 for MIF that was strongly potentiated by the MOG peptide extension, resulting in a novel therapeutic, DRα1-MOG-35-55, that within the limitations of the EAE model may have the potential to treat autoimmune diseases such as multiple sclerosis.

A novel method to measure HLA-DM-susceptibility of peptides bound to MHC class II molecules based on peptide binding competition assay and differential IC(50) determination

HLA-DM (DM) functions as a peptide editor that mediates the exchange of peptides loaded onto MHCII molecules by accelerating peptide dissociation and association kinetics. The relative DM-susceptibility of peptides bound to MHCII molecules correlates with antigen presentation and immunodominance hierarchy, and measurement of DM-susceptibility has been a key effort in this field. Current assays of DM-susceptibility, based on differential peptide dissociation rates measured for individually labeled peptides over a long time base, are difficult and cumbersome. Here, we present a novel method to measure DM-susceptibility based on peptide binding competition assays performed in the presence and absence of DM, reported as a delta-IC(50) (change in 50% inhibition concentration) value. We simulated binding competition reactions of peptides with various intrinsic and DM-catalyzed kinetic parameters and found that under a wide range of conditions the delta-IC(50) value is highly correlated with DM-susceptibility as measured in off-rate assay. We confirmed experimentally that DM-susceptibility measured by delta-IC(50) is comparable to that measured by traditional off-rate assay for peptides with known DM-susceptibility hierarchy. The major advantage of this method is that it allows simple, fast and high throughput measurement of DM-susceptibility for a large set of unlabeled peptides in studies of the mechanism of DM action and for identification of CD4+ T cell epitopes.

A novel family of human leukocyte antigen class II receptors may have its origin in archaic human species

HLA class II α and β chains form receptors for antigen presentation to CD4(+) T cells. Numerous pairings of class II α and β subunits from the wide range of haplotypes and isotypes may form, but most of these combinations, in particular those produced by isotype mixing, yielded mismatched dimers. It is unclear how selection of functional receptors is achieved. At the atomic level, it is not known which interactions of class II residues regulate selection of matched αβ heterodimers and the evolutionary origin of matched isotype mixed dimer formation. In this study we investigated assembly of isotype-mixed HLA class II α and β heterodimers. Assembly and carbohydrate maturation of various HLA-class II isotype-mixed α and β subunits was dependent on the groove binding section of the invariant chain (Ii). By mutation of polymorphic DPβ sequences, we identified two motifs, Lys-69 and GGPM-(84-87), that are engaged in Ii-dependent assembly of DPβ with DRα. We identified five members of a family of DPβ chains containing Lys-69 and GGPM 84-87, which assemble with DRα. The Lys/GGPM motif is present in the DPβ sequence of the Neanderthal genome, and this ancient sequence is related to the human allele DPB1*0401. By site-directed mutagenesis, we inspected Neanderthal amino acid residues that differ from the DPB1*0401 allele and aimed to determine whether matched heterodimers are formed by assembly of DPβ mutants with DRα. Because the *0401 allele is rare in the sub-Saharan population but frequent in the European population, it may have arisen in modern humans by admixture with Neanderthals in Europe.

Molecular dynamics simulations to provide insights into epitopes coupled to the soluble and membrane-bound MHC-II complexes

Epitope recognition by major histocompatibility complex II (MHC-II) is essential for the activation of immunological responses to infectious diseases. Several studies have demonstrated that this molecular event takes place in the MHC-II peptide-binding groove constituted by the α and β light chains of the heterodimer. This MHC-II peptide-binding groove has several pockets (P1-P11) involved in peptide recognition and complex stabilization that have been probed through crystallographic experiments and in silico calculations. However, most of these theoretical calculations have been performed without taking into consideration the heavy chains, which could generate misleading information about conformational mobility both in water and in the membrane environment. Therefore, in absence of structural information about the difference in the conformational changes between the peptide-free and peptide-bound states (pMHC-II) when the system is soluble in an aqueous environment or non-covalently bound to a cell membrane, as the physiological environment for MHC-II is. In this study, we explored the mechanistic basis of these MHC-II components using molecular dynamics (MD) simulations in which MHC-II was previously co-crystallized with a small epitope (P₇) or coupled by docking procedures to a large (P₂₂) epitope. These MD simulations were performed at 310 K over 100 ns for the water-soluble (MHC-II_w, MHC-II-P_7w, and MHC-II-P_22w) and 150 ns for the membrane-bound species (MHC-II_m, MHC-II-P_7m, and MHC-II-P_22m). Our results reveal that despite the different epitope sizes and MD simulation environments, both peptides are stabilized primarily by residues lining P1, P4, and P6-7, and similar noncovalent intermolecular energies were observed for the soluble and membrane-bound complexes. However, there were remarkably differences in the conformational mobility and intramolecular energies upon complex formation, causing some differences with respect to how the two peptides are stabilized in the peptide-binding groove.

Structural and dynamical insights on HLA-DR2 complexes that confer susceptibility to multiple sclerosis in Sardinia: a molecular dynamics simulation study

Sardinia is a major Island in the Mediterranean with a high incidence of multiple sclerosis, a chronic autoimmune inflammatory disease of the central nervous system. Disease susceptibility in Sardinian population has been associated with five alleles of major histocompatibility complex (MHC) class II DRB1 gene. We performed 120 ns of molecular dynamics simulation on one predisposing and one protective alleles, unbound and in complex with the two relevant peptides: Myelin Basic Protein and Epstein Barr Virus derived peptide. In particular we focused on the MHC peptide binding groove dynamics. The predisposing allele was found to form a stable complex with both the peptides, while the protective allele displayed stability only when bound with myelin peptide. The local flexibility of the MHC was probed dividing the binding groove into four compartments covering the well known peptide anchoring pockets. The predisposing allele in the first half cleft exhibits a narrower and more rigid groove conformation in the presence of myelin peptide. The protective allele shows a similar behavior, while in the second half cleft it displays a narrower and more flexible groove conformation in the presence of viral peptide. We further characterized these dynamical differences by evaluating H-bonds, hydrophobic and stacking interaction networks, finding striking similarities with super-type patterns emerging in other autoimmune diseases. The protective allele shows a defined preferential binding to myelin peptide, as confirmed by binding free energy calculations. All together, we believe the presented molecular analysis could help to design experimental assays, supports the molecular mimicry hypothesis and suggests that propensity to multiple sclerosis in Sardinia could be partly linked to distinct peptide-MHC interaction and binding characteristics of the antigen presentation mechanism.

Endoplasmic reticulum stress sensing in the unfolded protein response

Secretory and transmembrane proteins enter the endoplasmic reticulum (ER) as unfolded proteins and exit as either folded proteins in transit to their target organelles or as misfolded proteins targeted for degradation. The unfolded protein response (UPR) maintains the protein-folding homeostasis within the ER, ensuring that the protein-folding capacity of the ER meets the load of client proteins. Activation of the UPR depends on three ER stress sensor proteins, Ire1, PERK, and ATF6. Although the consequences of activation are well understood, how these sensors detect ER stress remains unclear. Recent evidence suggests that yeast Ire1 directly binds to unfolded proteins, which induces its oligomerization and activation. BiP dissociation from Ire1 regulates this oligomeric equilibrium, ultimately modulating Ire1's sensitivity and duration of activation. The mechanistic principles of ER stress sensing are the focus of this review.