Abstract
Bats serve as reservoirs for numerous zoonotic viruses, yet they typically remain asymptomatic owing to their unique immune system. Of particular significance is the MHC-I in bats, which plays crucial role in anti-viral response and exhibits polymorphic amino acid (AA) insertions. This study demonstrated that both 5AA and 3AA insertions enhance the thermal stability of the bat MHC-I complex and enrich the diversity of bound peptides in terms of quantity and length distribution, by stabilizing the 310 helix, a region prone to conformational changes during peptide loading. However, the mismatched insertion could diminish the stability of bat pMHC-I. We proposed that a suitable insertion may help bat MHC-I adapt to high body temperatures during flight while enhancing antiviral responses. Moreover, this site-specific insertions may represent a strategy of evolutionary adaptation of MHC-I molecules to fluctuations in body temperature, as similar insertions have been found in other lower vertebrates.
Similar content being viewed by others
Introduction
COVID-19 has resulted in immense global devastation and loss, with its origins traced back to a bat-borne coronavirus1,2. Bats harbor a multitude of viruses, including rabies3, Hendra virus (HEV)4,5, Nipah virus6,7, severe acute respiratory syndrome coronavirus8,9,10, Middle East respiratory syndrome coronavirus11,12,13 and Ebola virus14, all of which pose a threat to human and other mammals. Due to their diverse species composition15,16, widespread distribution17, and unique ability for sustained flight, bats are considered significant hosts of zoonotic viruses, contributing to the emergence of infectious diseases. Astonishingly, despite being carriers of numerous lethal viruses, bats do not display obvious pathological symptoms18,19,20,21. This remarkable resilience can be attributed to their exceptional immune system, which combines heightened host defense responses with immune tolerance, effectively combating viral infections while suppressing excessive inflammatory reactions22,23. Interestingly, bats have evolved immune tolerance not only in response to pathogen challenges, but also to adapt to the cellular and DNA damages induced by the intense metabolism and high body temperatures during flight. Bats have elevated expression of heat shock proteins (HSPs) to withstand high temperatures and oxidative stress associated with flight. Furthermore, bats have dampened the activity of the pattern recognition receptor STING, lost the PYHIN gene family, and suppressed the NLR-family pyrin domain containing 3 (NLRP3) and the inflammatory cytokine IL-1β, thus curbing the inflammation triggered by flight-induced DNA damage22. These traits should also help bats tolerate lethal viral infections, as the damage and death caused by infections are often also closely associated with high fevers and inflammation.
While research on adaptive immunity in bats remains limited, the major histocompatibility complex I (MHC-I) molecule stands out as a subject of interest24,25,26,27,28. The MHC-I molecule plays a crucial role in antiviral immunity by presenting antigenic peptides to activate specific CD8+ T cells to clear virus-infected cells, while its downregulation can activate NK cells and trigger inflammatory responses29,30,31. Bats possess distinctive characteristic in their MHC-I molecule, with prevalent 3- or 5-amino acid (3AA or 5AA) insertions at the N-terminal end of the peptide binding groove (PBG), which may generate a unique antiviral immune response. Peptidomics and crystal structure analyses of Ptal-N*01:01 have shed light on the impact of the 3AA insertion (MDL52-54) on peptide binding24,25,27. The 3AA insertion facilitates the formation of a salt-bridge network between the N-terminal end of the PBG of Ptal-N*01:01 and peptides containing aspartic acid as their first residue (P1-D), thereby enhancing their binding. Additionally, Ptal-N*01:01 is capable of binding long peptides of up to 15 amino acids. These studies indicated that the 3AA insertion can augment the binding of peptides to bat MHC-I molecules.
There are still some questions about the insertions in bat MHC-I. First, 5AA insertions are more prevalent in bat MHC-I alleles and their roles remain unknown. Second, the mechanism for how insertion affects bat MHC-I is not fully characterized. Current studies suggested that the effect of the insertion on bat MHC-I depends on the subsequent occurrence of a pair of residues with opposite charges. However, this pair is not present in all bat MHC-I molecules, suggesting that there should be other possibilities. Interestingly, the similar insertions in MHC-I have also been observed in ectothermic animals (bats and marsupials) and reptiles, all which have large fluctuations in body temperature28,32. All these insertions, including 5AA insertions, occur within the relatively unstable 310 helix of MHC-I, a region prone to conformational changes33,34. Based on these observations, it reasonable to speculate that these insertions are associated with the thermal stability of the peptide/MHC-I complex (pMHC-I).
In this study, we uncovered the peptide binding motif (PBM) and structures of the bat MHC-I molecule Mylu-B*67:01 and its mutant Mylu-B*67:01ΔMQQPW, featuring a 5AA insertion (MQQPW52-56). We demonstrated that both 5AA and 3AA insertions can enhance the formation and thermal stability of pMHC-I by stabilizing the 310 helix and may help pMHC-I to tolerate the high body temperatures during flight, but the insertion and the inserted bat MHC-I allele need to match each other. Furthermore, AlphaFold2 (AF2) predictions indicated that the effects of different insertions on bat MHC-I were similar to those observed in Mylu-B*67:01 or Ptal-N*01:01. These findings showed that the insertions could improove the thermal stability of bat pMHC-I, and significantly enhance the peptide binding capacity of bat MHC-I, including the quantity, length, and stability of bound peptides. We hypothesized that these inserts serve a dual purposes: first, to help pMHC-I tolerate the high temperatures of bat flight and avoid undue dissociation to cause unwanted immune responses; and second, to improve peptide presentation and clear virus-infected cells more efficiently.
Results
The peptide binding motif and crystal structure of Mylu-B*67:01 with MQQPW insertion
The distinguishing feature that sets bat MHC-I molecules apart from those of other mammals lies in the presence of three or five AA insertions within the α1 region of the heavy chain (HC)26,28. Among the 131 available bat MHC-I sequences in GenBank, the 5AA insertions were found to be the most prevalent, accounting for 72% (Supplementary Fig. 1a, b). These inserted sequences exhibit notable diversity, with MQQPW, MEGPW, and MEQPW being more commonly observed (Supplementary Fig. 1b). Notably, unlike the 3AA insertions, which are also observed in marsupial and reptile MHC-I molecules, the 5AA insertion appears to be unique to bat MHC-I. In addition, previous studies have shown that the 3AA insertion in the bat Ptal-N*01:01 enhances binding to peptides with aspartic acid as the first residue (P1-D), but requires the participation of aspartic acid at position 59 (D59) and arginine at position 65 (R65) in the formation of a salt bridge network24,25.
To explore the influence of the 5AA insertion on the structure and function of bat MHC-I, we focused on Mylu-B*67:01. Mylu-B*67:01 contains the most common 5AA insertion sequence, MQQPW, and also exhibits the D61/R67 pairing. The peptide-binding motif (PBM) of Mylu-B*67:01 was determined using our previously established RPLD-MS technique24. Initially, Mylu-B*67:01 and bat β2m were co-refolded in vitro with a library of random 9 peptides (R9Ps), followed by purification to obtain the pMylu-B*67:01-R9Ps complexes (Fig. 1a). The peptides bound within the complexes were eluted with a weak acid, sequenced using the De Novo LS-MS/MS technique, and ultimately, the PBM of Mylu-B*67:01 was visualized using Icelogo35 (Fig. 1b). The PBM exhibited a clear preference for Mylu-B*67:01 to bind peptides with a proline at the P2 position (P2-P), while hydrophobic amino acids with long side chains, like phenylalanine (F), leucine or isoleucine (denoted as L due to indistinguishability in mass spectra), methionine (M), or valine (V), were predominantly observed at the PΩ position (Fig. 1b). It is noteworthy that D/E occupies a relatively high proportion of the P1 position, although not as pronounced as P2-P. Three peptides matching the PBM were screened and synthesized from the spike protein of COVID-19, and their binding to Mylu-B*67:01 was confirmed through in vitro refolding (Fig. 1c), which validated the accuracy of the obtained PBM. Among these three peptides, P1-FPQSAPHGV exhibited highest refolding efficiency (Fig. 1c) and ultimately successfully formed crystals in complex with Mylu-B*67:01 (pMylu-B*67:01-P1). The crystal structures of the pMylu-B*67:01-P1 complex achieves a resolution of 2.2 Å. Detailed crystal parameters can be found in Table 1.
The crystal structure of pMylu-B*67:01-P1 unveils the intriguing consequences of the MQQPW insertion, forming an additional delicate helix nestled between the 310 Helix and α1 helix of HC (Fig. 1d). The inserted MQQPW resides at the junction of the 310 Helix and the emerging helix, creating a dainty loop connecting these two helical structures. The MQQPW insertion maintains a considerable distance from the peptide and does not engage in direct interactions (Fig. 1e). The presence of the R67 side chain induces the closure of the A pocket, thereby obstructing the exposure of the P1-F side chain. Notably, despite the close proximity of D61 and R67, no salt bridge was formed between them (Fig. 1e).
The preference for P2-P in the B pocket of Mylu-B*67:01 aligns with observed binding characteristics in other pMHC-I structures (57-61). Residues at positions 71 and 72 emerge as pivotal determinants in shaping the B pocket preference (Fig. 1f). The large aromatic side chain of F/Y72 poses a hindrance to the deep insertion of the P2 residue side chain into the pocket, while I71 establishes a favorable hydrophobic interaction environment. Consequently, P2-P aligns exceedingly well within this B-pocket configuration, emerging as a remarkable anchoring residue. Y7 is highly conserved in the vast majority of MHC-I molecules36, and the other constituent residues are blocked by F/Y72 from contacting P2 residue, so none of them are determinants of B pocket preference.
The F pocket of Mylu-B*67:01 closely resembles the monkey MHC-I molecule Mamu-A*01:0137, with a slight difference in a single residue (I/L100). Both molecules show a preference for binding PΩ-L, indicating similar anchoring residue preferences (Fig. 1g). Additionally, the F pocket of Mylu-B*67:01 also demonstrates affinity for other long-side chain hydrophobic residues (F/M/V), similar to the observed binding patterns in pHLA-A24var and pHLA-B*14:02 structures38,39. Detailed analysis reveals significant similarity in residues 86, 100, and 121 within the F pocket, along with conserved sites across most MHC-I molecules (89/123/128/129/148/151/152). These sites contribute to shaping the size and phobic landscape of the F pocket, creating an ideal environment for binding PΩ-F/L/M/V.
The inserted MQQPW enhances peptide binding but is not solely dependent on the assistance of D61/R67
To examine how the insertion influences the peptide presentation by Mylu-B*67:01, a mutant version of Mylu-B*67:01, named Mylu-B*67:01ΔMQQPW, was created by removing the MQQPW sequence. In vitro refolding experiments demonstrated that this mutant could still bind to the COVID-19th P1-P3 peptides, albeit with slightly reduced complex formation compared to the wild-type HC (Supplementary Fig. 2). This difference in complex formation became more pronounced when co-refolding with R9Ps (Fig. 2a). RPLD-MS analysis revealed that the PBM of Mylu-B*67:01ΔMQQPW was similar to that of the wild-type HC (Fig. 2b). However, the preference for P2 and PΩ anchoring residues was significantly diminished, and the difference of P1 residues was relatively obvious (Fig. 2c). Based on previous studies of Ptal-N*01:01, it was hypothesized that the MQQPW insertion might affect the formation of a salt bridge between D59 and R65, consequently impacting the anchoring of P1 residues. To explore this hypothesis, the COVID-19 P1 peptide was mutated at the P1-F position, generating P1-T and P1-E mutants that were more compatible with the P1 residue preference of Mylu-B*67:01ΔMQQPW (Fig. 2d). However, both wild-type and mutant COVID-19 P1 peptides exhibited better binding to Mylu-B*67:01 than Mylu-B*67:01ΔMQQPW (Fig. 2d). The crystal structures of wild-type and mutant P1 peptides in complex with Mylu-B*67:01 and Mylu-B*67:01ΔMQQPW were solved, except pMylu-B*67:01ΔMQQPW-P1-E which may be attributed to by its low thermal stability (Fig. 2d). The details of these crystal structures were listed in Table 1. The binding patterns of peptides P1-F and P1-T with both Mylu-B*67:01 and the variant Mylu-B*67:01ΔMQQPW remain largely unchanged, with the primary structural variations localized around the region where the MQQPW sequence has been inserted (Fig. 2e). We refined the structures after the deletion of MQQOW and generated unbiased 2Fo-Fc electron density maps to confirm that the structural changes caused by the insertion are credible (Fig. 2f). Clear electron density maps enable us to conduct further comparisons and analyses.
The structures show that Mylu-B*67:01ΔMQQPW has a classical 310 helix conformation, whereas the insertion of MQQPW introduces a new minor helix and slightly delays the formation of the α1 helix after 61D (Fig. 3a). The insertion of MQQPW brings D61 closer to R67 (Fig. 3b), but the salt bridge network is only formed in the presence of P1-E (Fig. 3a). The salt bridge network induced by P1-E makes a significant shift in the minor helix formed by MQQPW (Fig. 3c). The interactions between P1-E and 61D/67 R is similar to that of Ptal-N*01:01 with an inserted MDL (Fig. 3a, d)24,25. From the surface pattern, the change in the position of 61D allows it to participate in the formation of the A pocket and to form a salt bridge network with P1-E via 67 R (Fig. 3e). Combined with the PBM and CD results (Fig. 1b, 2d), we believed that 5AA insertion, like 3AA insertion, can also enhance the binding of bat MHC-I molecules to peptides bearing P1 acidic residue with the help of 61D/67 R.
At the same time, some evidence also implied that the impact of MQQPW insertion is not entirely dependent on the salt bridge network formed by 61D/67 R. MQQPW insertion could enhance the binding of peptide P1-F and peptide P1-T to Mylu-B*67:01 (Fig. 2d), but 61D and 67 R do not form salt bridges in their complex structures. (Fig. 3a). The P1-E mutant peptide was less thermally stable with Mylu-B*67:01 compared to the original P1-F peptide (Fig. 2d), suggesting that the acidic residues are not the only preference for the P1 position. This was also supported by the PBM comparison of Mylu-B*67:01 and Mylu-B*67:01ΔMQQPW (Fig. 2c). Additionally, it is worth noting that many bat MHC-I sequences lack the corresponding charged residue pairing after the insertion of 3AA or 5AA (Supplementary Fig. 1c). The pairing of D59 and R65 accounts for 50% of the 3AA insertion sequences. The pairing of this residue combination (D61 and R67) accounts for about 17% of the 5AA insertion sequence, and the pairing of residues with the potential to form a salt bridge accounts for 31% (Supplementary Fig. 1c). Together, these findings indicated that the insertion of 5AA or 3AA in bat MHC-I could also enhance binding to peptides in a manner not involving the 61D/67 R, especially to peptides with non-acidic residues at the P1 position.
Insertion can enhance the stability of bat pMHC-I by strengthening the 310 helix region
The MQQPW insertion is precisely located within the 310 helix region, which is the most unstable part of the α-helix structure of the MHC-I molecule33. The B-factor can be used to characterize the stability of the structure. We performed Min-Max normalization of the B factors in the peptide binding regions of Mylu-B*67:01 and Mylu-B*67:01ΔMQQPW to compare the effect of MQQPW insertion on the B factors of the relevant regions (Fig. 4a). In Mylu-B*67:01ΔMQQPW, the 310 helix region is one of the regions with the highest B factors, and in Mylu-B*67:01, the insertion is able to reduce the B factors of neighboring residues. The B-factor putty model can visualize how the B-factors change in the structure, and the comparative analysis revealed that the insertion of MQQPW reduced the B-factor of the 310 helix region, indicating an improvement in the structural stability of this region (Fig. 4b). This was further supported by the analysis of hydrogen bonding in this region (Fig. 4c). In Mylu-B*67:01, there are a total of 12 intra-chain hydrogen bonds in this region to maintain the structure, whereas in Mylu-B*67:01ΔMQQPW, there are only four hydrogen bonds. The conformation of W52 changes significantly during peptide exchange34, suggesting that improving the stability of W52 contributes to the stabilization of the complex. Within the helix formed by the MQQPW insertion, the inserted W58 forms a hydrogen bond with W52 on the 310 helix, which helps to strengthen the link between the two smaller helices.
Subsequently, we examined whether the 3AA insertion had a similar effect on Ptal-N*01:01. Their B factors in the peptide binding regions were compared in (Fig. 4d, e). The analysis revealed that the MDL insertion, similar to MQQPW, reduce the B factor of the 310 helix region, but through a different mechanism. The insertion of MDL extends the 310 helix and forms 6 hydrogen bonds to stabilize the structure, which is four more than in the ΔMDL mutant. (Fig. 4f). The inserted L55 forms a hydrogen bond with W52 on the 310 helix. So we believed that both 3AA and 5AA insertions can enhance the stability of bat pMHC-I by strengthening the 310 helix region, independent of D61/R67.
Insertion may help bat pMHC-I to tolerate high body temperature during flight and improve the peptide binding ability, but requires a match to the inserted allele
Both 3AA and 5AA insertions have been shown to enhance the thermal stability of bat pMHC-I with specific peptides25 (Fig. 2d). However, relying solely on the thermal stability measurements of a few pMHC-I individuals may not adequately reflect the overall impact of the insertion. To address this, we assessed the thermal stabilities of mixed pMylu-B*67:01-R9Ps and pMylu-B*67:01ΔMQQPW-R9Ps. The circular dichroism (CD) analysis revealed that the insertion of MQQPW increased the Tm value of the mixed pMHC-Is by approximately 6.5 °C (Fig. 5a), which closely matches the magnitude of temperature increase observed during bat flight. We further examined the changes in these mixed bat pMHC-Is after incubation at the body temperature experienced by bats during flight (41 °C) for 30 min and 60 min (Fig. 5b). As expected, 41 °C exceeded the Tm value of the mixed bat pMHC-Is, resulting in the dissociation of more than half of the pMHC-Is, consistent with the CD results. However, at both time points, a higher proportion of pMylu-B*67:01-R9Ps remained compared to pMylu-B*67:01ΔMQQPW-R9Ps, indicating that the insertion of MQQPW indeed improves the tolerance of bat pMHC-Is to high body temperature. Moreover, the thermal stability of MHC-I molecules plays a crucial role in peptide presentation40,41. The insertion of MQQPW leads to structural solidification, thereby enhancing the peptide-binding capacity of bat MHC-I. In vitro experiments demonstrated that Mylu-B*67:01 exhibited a significantly higher binding ability for both short and long peptides compared to Mylu-B*67:01ΔMQQPW (Fig. 5c).
The insertion of MDL also enhanced the thermal stability of pPtal-N*01:01-R9Ps and improved its tolerance to 41 °C, although to a lesser extent compared to MQQPW (Fig. 5d, e). The MDL insertion also significantly increased the number and length of peptides that Ptal-N*01:01 could bind (Fig. 5f).
Previous studies have indicated that the insertion of MDL does not enhance the peptide-binding capacity of HLA-A*02:0125. Additionally, insertions in bat MHC-I alleles exhibit different sequences and lengths (Supplementary Fig. 4a, b), suggesting that the effects of insertions may be specific to certain species or alleles. To investigate this hypothesis, we selected the bat MHC-I allele Mylu-B*57:01, which does not possess any amino acid insertion. We introduced MDL and MQQPW insertions into Mylu-B*57:01 and compared their impact on peptide binding (Fig. 5g). The results demonstrated that only MDL was able to enhance the binding of Mylu-B*57:01 to random nonapeptides and slightly improve the thermal stability of pMHC-I. In contrast, MQQPW did not exhibit any enhancing effect and significantly reduced the thermal stability of pMHC-I (Fig. 5h).
These findings indicated that both the 3AA and 5AA insertions may enhance the ability of bat pMHC-I to tolerate high body temperature during flight and improve the peptide binding, but the insertion and the inserted bat MHC-I allele need to match each other.
Analysis of the impact of diverse insertions on the structure of bat MHC-I with AlphaFold2
Investigating the effect of diverse insertions observed in numerous bat MHC-I alleles by crystallographic methods, as done for Mylu-B*67:01 and Ptal-N*01:01, is an impractical approach. The emergence of AlphaFold2 (AF2) has significantly enhanced the accuracy of structure prediction42, offering a potential solution to solve this problem. Prior to AF2, predictive methods like Swiss Model could predict the structure of MHC-I, but their accuracies was limited by sequence similarity between the target sequence and template, particularly for MHC-I molecules with substantial sequence variations across species, such as Mylu-B*67:01 (Supplementary Fig. 3). Moreover, previous tools struggled to accurately predict the complex structure of pMHC-I, particularly the peptide-MHC-I binding. DeepMind, utilizing AF2, has developed a novel tool called AF-Multimer43 to predict complex structures. The main question to address is whether AF2 and AF-Multimer can effectively predict the structure of pMHC-I, particularly in unknown animal pMHC-I. We utilized AF2 to predict the structures of Mylu-B*67:01 HC and Mylu-B*67:01ΔMQQPW HC, as well as the structures of pMylu-B*67:01 and pMylu-B*67:01ΔMQQPW using AF-Multimer (Fig. 6a). AF2 accurately prediced the formation of the additional small helix by MQQPW insertion with a high confidence level (pLDDT > 90) comparable to that of the rest of the structure. The predicted structures closely resemble the actual structures (Fig. 6b), exhibiting RMSD values of 0.723 Å for pMylu-B*67:01 and 0.707 Å for pMylu-B*67:01ΔMQQPW. Remarkably, AF-Multimer’s prediction of the peptide binding conformation closely matched the true conformation. The above results indicated that AF2 can be used to predict the structure of bat MHC-I. Consequently, we employed AF2 to investigate the various AA insertions in bat MHC-I alleles (Fig. 6c). The predicted structures of diversve 5AA insertions and 3AA insertions are similar to MDL and MQQPW respectively.
Sequence alignment revealed that among the 5AA insertions, W58 was highly conserved, except in four sequences (XP_0244236.1/XP_036985990.1/XP_037023937.1/XP_0194562) (Supplementary Fig. 4). However, predictions indicated that the α2 helix of these four bat MHC-Is varied significantly from the classical MHC-I structure, making the accuracy of their sequences uncertain (Supplementary Fig. 4). The remaining 5AA insertions could form a secondary helix and stabilize the 310 helix through a hydrogen bond between W58 and W52, similar to Mylu-B*67:01 (Fig. 6c). Furthermore, all 3AA insertions could extend the 310 helix, and the 55th residue could maintain the helix by forming a hydrogen bond with W52, as observed in the case of Ptal-N*01:01 (Fig. 6c). The predictions also indicated that salt bridges could form when residues with opposite charges were paired at positions 61/67 (5AA insertion number) or 59/65 (3AA insertion number). Therefore, we suggested that despite sequence differences, insertions of the same length have similar effects on bat MHC-I and that all insertions should have the ability to enhance 310 helix stability.
Discussion
The MHC-I molecule plays a key role in antiviral immunity, and the presence of prevalent and polymorphic 5AA or 3AA insertions in bat MHC-I molecules has attracted particular attention32. However, existing studies have examined 3AA insertions with Ptal-N*01:01 and have not addressed 5AA insertions, which are unique and more prevalent in bat MHC-I molecules. Here, we investigated the 5AA insertion with Mylu-B*67:01. The MQQPW insertion introduces a new small helix and pushes D61 closer to R67, allowing them to form a salt bridge network with the P1 acidic residue of the bound peptide (Fig. 3). This is similar to the effect of MDL insertion on Ptal-N*01:01 (Fig. 3d), suggesting that with the help of a specific pair of oppositely charged residues after insertion, both 5AA and 3AA insertions enhance the binding of bat MHC-I to peptides bearing P1 acidic residues24,25.
Nevertheless, the above findings do not fully explain the effects of insertion on bat MHC-I. First, the MQQPW insertion also enhances the binding of Mylu-B*67:01 to the peptides with P1 non-acidic residues, but no salt bridge is formed between 61D and 67 R (Figs. 2, 3 and Supplementary Fig. 2). Second, the oppositely charged residue pairing is not present in all bat MHC-I sequences with the insertion (Supplementary Fig. 1c). By comparing the insertion-induced structural changes of bat pMHC-I with the B factor putty mode, a presentation that indicates structural stability, we found a possible new explanation. Both 5AA and 3AA insertions can reduce the B factors of the 310 helix region (Fig. 4), which is an unstable region of the antigen-binding domain that needs to flip its conformation during peptide assembly33,34. The 3AA insertion extends the 310 helix and increases its stability, while the 5AA insertion forms an additional small helix that links to the 310 helix. Both insertions bring more intra-chain hydrogen bonds and interact with W52 on the 310 helix, which undergoes a significant deflection during peptide assembly34. We believed that these interactions reduce the conformational flexibility of the 310 helix and increase the difficulty of peptide dissociation, thus enhancing the thermal stability of bat pMHC-I. The new interpretation we proposed here is more general, as it not only involves the insertion itself and applies to all binding peptides.
The structural stability of bat pMHC-I is strengthened by the insertion, which should be reflected in its increased thermal stability. This can be determined by in vitro thermal incubation experiments, and the determination of mixed complexes formed by MHC-I with multiple peptides can avoid bias due to the use of a few specific peptides and better reflect the overall differences. Jappe, E. C. et al. previously used mixed pMHC-Is captured from cells40, whereas we used mixed complexes of MHC-I with random peptides formed by in vitro refolding. The results showed that either 3AA or 5AA insertion can improve the thermal stability of bat pMHC-I (Fig. 5a, d) and reduce degradation upon incubation at 41 °C (Fig. 5b, e), which is the body temperature that bats can reach in flight We also confirmed that the mixed pMHC-Is can be used to obtain definite Tm values by CD assay. This will help to assess the overall peptide binding capacity of MHC-I, especially when studying MHC-I molecules in animals with limited conditions.
In addition to enhancing the binding of peptides with P1 acidic residues, the insertions also enhance the ability of bat MHC-I molecules to present diverse peptides. Deleting the insertion reduces the number of random peptides that can bind to Mylu-B*67:01 or Ptal-N*01:01 (Fig. 5c, f). This insertion does not affect pockets of anchoring residues, such as the B and F pockets, but can prevent peptide dissociation by strengthening the 310 helix region, thereby increasing the overall number of bound peptides. More eluted peptides helps the identification of the PBM of Mylu-B*67:01 by RPLD-MS (Fig. 1a). The peptidome eluted from the cells showed that Ptal-N*01:01 can present peptides up to 15 AA24. Our in vitro results showed that Ptal-N*01:01 could bind more 12-16AA peptides than the deletion mutant (Fig. 5f). Interestingly, Mylu-B*67:01 might have a strong binding capacity for long peptides, as the number of bound long peptides did not drop significantly compared to 9AA peptides (Fig. 5c). We speculated that the insertion improves the peptide binding ability of bat MHC-I molecules and helps bats produce effective antiviral T-cell responses to clear viral infections.
The wide diversity of insertions in bat MHC-I molecules (Supplementary Fig. 1) and their different effects on the same bat MHC-I molecule (Fig. 5g, h) suggest that the insertions and the inserted bat MHC-I have undergone selection during evolution. Recent advances in protein structure prediction allowed us to study diverse insertions on a large scale. Using AF2 and AF2-Multimer, which can accurately predict the structure of pMylu-B*67:01, we predicted the structure of all available bat MHC-I molecules. We did not used the structures of Ptal-N*01:01 and Mylu-B*67:01 to train AF2, but the predicted structures show that all diverse insertions can stabilize the 310 helix, similar to Ptal-N*01:01 and Mylu-B*67:01. AF2 and AF2-Mulimer have shown unique advantages in the study of HLA44,45, and our study suggested that they can also be applied in the study of MHC-I structure in unknown animals to help identify epitopes and develop vaccines.
The reason why bats carry multiple viruses (e.g., coronaviruses) but do not develop the disease may be due to the ability of their immune system to enhance the host defense response while suppressing excessive inflammatory responses caused by the intense physiological metabolism and high temperatures during flight. Some immune molecules associated with inflammation, such as STING, NLRP3, and IL-1β, are dampened or suppressed so as not to react to the damage caused flight. On the other hand, high fever and a strong inflammatory response are important factors in the pathological damage caused by viral infections, so bats are also very resistant to viral infections. It is an interesting question how the bat immune system adapts to flight and affects the antiviral response. Our study suggested that the insertions of MHC-I maybe related to this issue. We proposed the following hypothesis (Fig. 7): To adapt to the extreme temperature fluctuations experienced during flight, bat MHC-I molecules have developed additional insertions at the 310 helix. These insertions serve to enhance the stability of this region, which is prone to conformational changes. By preserving the stability of pMHC-I, these insertions may help prevent unnecessary inflammation, such as the activation of NK cells triggered by pMHC-I degradation29,30,31. These insertions not only contribute to the stability of the complex during flight but also enhance the antigen presentation ability of bat MHC-I molecules, which may result in a stronger T-cell immune response and improved clearance of viral infections. The dual role of the insertion in bat MHC-I appears to be in line with the features of the bat immune system and favors bats carrying viruses without getting sick. However, further research is needed to fully understand the implications of these insertions and their role in bat immune responses.
In fact, the 3AA insertion is also present in MHC-I molecules of lower animals such as reptiles and marsupials, and the effect on peptide binding is similar to the case of Ptal-N*01:0125. Bats and marsupials are heterotherms, and reptiles are ectotherm, which is less able to regulate body temperature than birds and higher mammals, and therefore have a great variation in body temperature. The presence of insertions in their MHC-I enhances the thermal stability of pMHC-Is, probably to adapt to the large changes in body temperature and to maintain the stability of MHC-I molecular function. To date, the 5AA insertions have been only found in bat MHC-Is, which are more prevalent and highly diverse compared to the 3AA insertions (Supplementary Fig. 1). In suitable inserted bat MHC-I, 5AA insertion may have a stronger enhancement on the thermal stability of pMHC-I (Fig. 5a) compared to insertion of 3AA (Fig. 5d). Bats, due to the intense temperature rise during flight, may had experienced stronger selection pressure than other lower animals, leading to 5AA insertions and higher diversity. These implied that the insertions in the MHC-Is of lower vertebrates might have evolved in response to their drastic changes in body temperature.
Methods
Sequence retrieval analysis and synthesis
The sequences of 131 MHC-I (including predicted genes) of bats were retrieved from the NCBI database (https://www.ncbi.nlm.nih.gov/) (Supplementary Table 1), and the sequences of human MHC-I were retrieved from the Immunopolymorphism Database (IPD) (www.ebi.ac.uk/ipd/mhc), and sequence comparisons were performed using MEGA-X46 and Jalview47.
Among the MHC I sequences used in this paper, the extracellular region expression plasmids of Ptal-N*01:01 (GenBank No. KT987929), Ptal-N*01:01ΔMDL, and bat β2m (GenBank No. XP_006920478.1)24 was previously constructed in our laboratory. The amplification product of the extracellular structural domain of Mylu-B*67:01 (GenBank number: XP_006107180.2), Mylu-B*57:01(GenBank number: XP_023604095.1) was cloned into the pET-21a vector (Novagen). To investigate the functions of Met52, Gln53, Gln54, Pro55, and Trp56, overlapping PCR was used (using the following primers:5’-TGCTCCTTGGGTTGAACAGGAAGATCCGGAATA-3’,5’-TGTTCAACCCAAGGAGCACGAGGTTCGGCTTTAGG-3’.), the five amino acid deletion mutant sequence Mylu-B*67:01ΔMQQPW was constructed and cloned into the pET-21a vector (Novagen).To investigate the effect of three and five amino acid insertions on bat MHC I molecules, we constructed two mutant sequences: MYLU-B-57-3 and MYLU-B-57-5 using overlapping PCR(using the following primers: MYLU-B-57-3:5’-GGATGGATTTAATG
GAAAGAGAGGAACCGGGTTA-3’,5’-TTCCATTAAATCCATCCAAGGGGCACGTGGCT-3’,MYLU-B-57-5:5’-ATGCAGCAGCCGTGGATGGAAAGAGAGGAACCGGGTTA-3’,5’- ATCCACGGCTGCTGCATCCAAGGGGCACGTGGC −3’.) and cloned them into the pET-21a vector (Novagen).
Peptide synthesis
The peptides with known sequences used in this experiment are listed in Supplementary Table 2. Eleven synthetic random peptide libraries, Ran_7Xsplitted -Ran_17Xsplitted (where X is a random amino acid other than cysteine). Synthesize as before to ensure its randomness48. All peptides were synthesized and purified to 99% by reverse-phase high-performance liquid chromatography (HPLC) and mass spectrometry (SciLight Biotechnology). All peptides are stored as a lyophilized powder in a −80 °C refrigerator and dissolved in dimethyl sulfoxide (DMSO) before use.
The constructed plasmids were transformed into E. coli strain BL21 (DE3). After prokaryotic expression of the plasmid and purification to inclusion bodies, the purified protein was dissolved in 6 M guanidine hydrochloride at a final concentration of 30 mg/ml and stored at −20 °C as described previously49,50. MHC heavy chain, β2m, and peptide were added slowly dropwise to the refolding buffer (100 mM Tris-HCl (pH 8.0), 2 mM EDTA, 400 mM L-Arg-HCl (JS0276, HongKong JiSiEnBei International Trade Co., Limited), 0.5 mM oxidized glutathione and 5 mM reduced glutathione) in a molar ratio of 1:1:3. After proteins were magnetically stirred in a revival buffer system for 24 h at 4 °C, the complexes were concentrated and exchanged into a buffer of 20 mM Tris-HCl (pH 8.0) and 50 mM NaCl, and the complexes were further purified by gel filtration chromatography and anion-exchange chromatography (Superdex 200 16/60 column and Resource Q anion exchange column, GE Healthcare)51.
Isolation and Identification of High-Affinity Peptides
The peptide-containing fraction of the complex was acidified with 0.2 N acetic acid at 65 °C for 30 min as described previously24,51. Then, the peptide-containing fraction was concentrated through a 3 kDa filter and desalted using a desalting column (WAT094225, waters). The desalting procedure is shown below: 1. The column was activated twice using 200 µl of 0.5% TFA, and 80% ACN 2. Equilibrate the column twice with 200 µl of 0.5% TFA. 3. Replace the collection tube, load the sample onto the column, and repeat. 4. Wash the sample with 200 µl of 0.5% TFA. 5. Replace the collection tube and elute the sample from the column with 200 µl of 0.5% TFA, 80% ACN. 6. Spin dry and store for use.
The EasyNano LC 1000 system (Thermo Fisher Scientific, San Jose, CA) was used to fractionate the peptides. Firstly, the desalted peptides were loaded onto an in-house constructed trap column (5-μm pore size, 150-μm i.d. × 3-cm length, 120 Å) before being separated on a custom-made C18 column (3-μm pore size, 75-μm i.d. × 15-cm length, 100 Å) with a flow rate of 450 μl/min. A linear gradient spanning 60 min was performed, starting with 3% B (0.1% formic acid in acetonitrile [v/v]/97% A (0.1% formic acid in H2O [v/v]) and increasing to 6% B in 8 min, 22% B in 37 min, 35% B in 8 min, and 100% B in 2 min before holding at 100% B for 5 min. The Q Exactive HF (Thermo Fisher Scientific, Bremen) acquired the MS data in data-dependent acquisition mode, with the top 20 precursors by intensity from mass range m/z 300 to 1800 being sequentially fragmented using higher-energy collisional dissociation at a normalized collision energy of 27. A dynamic exclusion time of 20 s was used, with automatic gain control set to 3e6 and 1e for MS1 and MS2, respectively, and a resolution of 120 and 30 K for MS1 and MS224,51.
Peptides were isolated from each spectrum (FDR = 1%) via de novo sequencing in Peaks Studio software based on the spectrum information. The following parameters were set: enzyme was set to nonspecific, variable modifications of oxidation (M)/deamidation (N, Q) were adjusted, peptide mass tolerance was approximately ±10 ppm, and fragment mass tolerance was set to 0.02 Da. After being adjusted by the detection threshold (score ≥ 50). All the data were deposited and conformed to community standards52.All nine-amino acid peptides were selected, and their binding motifs were mapped using Icelogo software. The appropriate antigenic peptide in the spike protein of COVID-19 was selected based on the obtained peptide binding motifs.
Determination of protein thermostability using CD spectroscopy
The sample concentration of all complexes to be measured was adjusted to 0.2 mg/ml.CD spectra at 218 nm were measured on a Chirascan spectrometer (Applied Photophysics) using a thermostatically controlled cuvette at temperature intervals of 0.2 °C at an ascending rate of 20 to 80 °C. The spread fraction (%) is expressed as (θ - θ a)/(θ a - θb), where θ a and θb are the average residual ellipticity values in the fully folded and fully expanded states, respectively. The degeneracy curves were generated by nonlinear fitting with OriginPro 8.053 (OriginLab). Tm was calculated by fitting the data to the degeneracy curves and using the deformation-determining derivatives.
Determination of protein thermostability using an in vitro heating method
At the end of the refolding of the p/MHC I molecules, the complexes were purified in the first step using a gel filtration chromatography column (GE), the target complex molecules were collected and the samples were concentrated to 2 ml using a 3 kDa ultrafiltration concentration tube. Each sample was divided into equal triplicates and incubated at 41 °C for 0 min, 30 min, and 60 min, respectively, and then subjected to gel filtration chromatography, and the harvested target complex molecules were quantified.
Crystallization, data collection, and processing
Three peptides, FPQSAPHGV/TPQSAPHGV/EPQSAPHGV, were utilized to form crystals with Mylu-B*67:01/Mylu-B*67:01ΔMQQPW and β2m, respectively. The crystals were initially screened using the sitting-drop method and the Crystallization Screening Kit (Hampton Research). Purified protein complexes were diluted to 2 and 4 mg/ml with molecular sieve buffer and then mixed with the crystallization screening kit buffer in a 1:1 volume ratio. pMylu-B*67:01/P9-FPQ, pMylu-B*67:01/P9-EPQ, and pMylu-B*67:01/P9-TPQ were grown as relatively thin slices in 0.2 M ammonium acetate, 0.1 M BIS-TRIS pH 7.1, and 20% w/v polyethylene glycol 3350. pPtal-N*01:01ΔMDL/HEV-1 were grown as relatively thick slices in 0.2 M sodium citrate tribasic dihydrate, 0.1 M BIS-TRIS pH 8.3, and 20% w/v polyethylene glycol 3350. Diffraction data were collected at Beamline BL18U1/BL19U of the Shanghai Synchrotron Radiation Facility using an R-AXIS IV + + imaging plate detector. The collected diffraction data were indexed and integrated using iMosflm54, and scaling and merging were performed using the CCP4i suite55.
The collected data were processed using CCP4i to determine the structure and calculate the phase via MOLREP56 and phasers57 (employing the molecular substitution method58 and using Ptal-N*01:01 as a structural template) to obtain the model and coordinates. After model building in COOT59 and several rounds of refinement in Refmac560, further refinement was conducted using phenix refine61. Subsequently, the MolProbity62 tools were utilized in Phenix for assessing model quality63 (refer to Structure data table). Structure-correlation diagram analysis and plotting were performed using PyMOL64 (Schrödinger, LLC).
Structure prediction using AlphaFold
Using 131 sequences of bat MHC class I molecules (with a signal peptide and transmembrane region removed) and bat β2m molecules obtained by downloading from the NCBI database, protein-protein interaction prediction was performed with the help of AlphaFold242 (https://colab.research.google.com/github/sokrypton/ColabFold/blob/main/AlphaFold2.ipynb) and AlphaFold-Multimer43 (https://github.com/deepmind/alphafold) online server, AlphaFold-Multimer uses a weighted combination of pTM and ipTM as the confidence indicator of the models, and finally selecting the model with the highest pLDDT value was used as our template for structure comparison. The prediction structure confidence level is specified in the following equation:
Statistics and Reproducibility
The relevant experiments covered in this manuscript are not statistically relevant, and the data obtained regarding the in vitro refolding experiments were repeated more than three times, with consistent results obtained each time.
Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.
Data availability
The coordinate and structure factors presented in this article have been submitted to the Protein Data Bank65 (https://www.rcsb.org/) under accession numbers 8HSM,8HSO,8HSW,8HT1, and 8HT9. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE66 (https://www.ebi.ac.uk/pride/) partner repository with the dataset identifier PXD047090,PXD047017,and PXD046952. All other data are available from the corresponding author on reasonable request.
References
Zhou, P. et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature 579, 270–273 (2020).
Andersen, K. G., Rambaut, A., Lipkin, W. I., Holmes, E. C. & Garry, R. F. The proximal origin of SARS-CoV-2. Nat. Med 26, 450–452 (2020).
Metchnikoff, E., Weinberg, M., Pozerski, E., Distaso, A. & Berthelot, A. Rousettes et microbes. Ann. Inst. Pasteur (Paris) 23, 61 (1909).
Halpin, K., Young, P. L., Field, H. E. & Mackenzie, J. S. Isolation of Hendra virus from pteropid bats: a natural reservoir of Hendra virus. J. Gen. Virol. 81, 1927–1932 (2000).
Young, P. L. et al. Serologic evidence for the presence in Pteropus bats of a paramyxovirus related to equine morbillivirus. Emerg. Infect. Dis. 2, 239–240 (1996).
Rahman, S. A. et al. Characterization of Nipah Virus from Naturally Infected Pteropus vampyrus Bats, Malaysia. Emerg. Infect. Dis. 16, 1990–1993 (2010).
Chua, K. B. et al. Isolation of Nipah virus from Malaysian Island flying-foxes. Microbes Infect. 4, 145–151 (2002).
Ge., X.-Y. et al. Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor. Nature 503, 535–538, https://doi.org/10.1038/nature12711 (2013).
Li, W. et al. Bats are natural reservoirs of SARS-like coronaviruses. Science 310, 676–679 (2005).
Lau, S. K. et al. Severe acute respiratory syndrome coronavirus-like virus in Chinese horseshoe bats. Proc. Natl Acad. Sci. USA 102, 14040–14045 (2005).
Wang, Q. et al. Bat origins of MERS-CoV supported by bat coronavirus HKU4 usage of human receptor CD26. Cell Host Microbe 16, 328–337 (2014).
van Boheemen, S. et al. Genomic characterization of a newly discovered coronavirus associated with acute respiratory distress syndrome in humans. mBio 3, https://doi.org/10.1128/mBio.00473-12 (2012).
Müller, M. A. et al. Human Coronavirus EMC Does Not Require the SARS-Coronavirus Receptor and Maintains Broad Replicative Capability in Mammalian Cell Lines. mBio 3, e00515–00512 (2012).
Leroy, E. M. et al. Fruit bats as reservoirs of Ebola virus. Nature 438, 575–576 (2005).
Upham, N. e. a. ASM’s Mammal Diversity Database, Accessed on 02/27/2023).
Simmons, N. B. & Cirranello, A. L. Bat Species of the World: A Taxonomic and Geographic Database, Accessed on 02/27/2023).
Nowak, R. M. & Walker, E. P. Walker’s bats of the world. (Johns Hopkins University Press, 1995).
Munster, V. J. et al. Replication and shedding of MERS-CoV in Jamaican fruit bats (Artibeus jamaicensis). Sci. Rep. 6, 21878 (2016).
Watanabe, S. et al. Bat coronaviruses and experimental infection of bats, the Philippines. Emerg. Infect. Dis. 16, 1217–1223 (2010).
Middleton, D. J. et al. Experimental Nipah virus infection in pteropid bats (Pteropus poliocephalus). J. Comp. Pathol. 136, 266–272 (2007).
Swanepoel, R. et al. Experimental inoculation of plants and animals with Ebola virus. Emerg. Infect. Dis. 2, 321–325 (1996).
Irving, A. T., Ahn, M., Goh, G., Anderson, D. E. & Wang, L. F. Lessons from the host defences of bats, a unique viral reservoir. Nature 589, 363–370 (2021).
Banerjee, A. et al. Novel Insights Into Immune Systems of Bats. Front. Immunol. 11, https://doi.org/10.3389/fimmu.2020.00026 (2020).
Qu, Z. et al. Structure and Peptidome of the Bat MHC Class I Molecule Reveal a Novel Mechanism Leading to High-Affinity Peptide Binding. J. Immunol. 202, 3493–3506 (2019).
Lu., D. et al. Peptide presentation by bat MHC class I provides new insight into the antiviral immunity of bats. PLoS Biol 17, https://doi.org/10.1371/journal.pbio.3000436 (2019).
Abduriyim, S., Zou, D. H. & Zhao, H. B. Origin and evolution of the major histocompatibility complex class I region in eutherian mammals. Ecol. Evolution 9, 7861–7874 (2019).
Wynne, J. W. et al. Characterization of the Antigen Processing Machinery and Endogenous Peptide Presentation of a Bat MHC Class I Molecule. J. Immunol. 196, 4468–4476 (2016).
Ng, J. H. et al. Evolution and comparative analysis of the bat MHC-I region. Sci. Rep. 6, 21256 (2016).
Shang, Q. N. et al. Functional Competence of NK Cells via the KIR/MHC Class I Interaction Correlates with DNAM-1 Expression. J. Immunol. 208, 492–500 (2022).
Yokoyama, W. M. & Kim, S. Licensing of natural killer cells by self-major histocompatibility complex class I. Immunological Rev. 214, 143–154 (2006).
Anfossi, N. et al. Human NK cell education by inhibitory receptors for MHC class I. Immunity 25, 331–342 (2006).
Wang, P. Y. et al. Peptide Presentations of Marsupial MHC Class I Visualize Immune Features of Lower Mammals Paralleled with Bats. J. Immunol. 207, 2167 (2021).
Wieczorek, M. et al. Major Histocompatibility Complex (MHC) Class I and MHC Class II Proteins: Conformational Plasticity in Antigen Presentation. Front Immunol. 8, 292 (2017).
Mage, M. G. et al. The peptide-receptive transition state of MHC class I molecules: insight from structure and molecular dynamics. J. Immunol. 189, 1391–1399 (2012).
Maddelein, D. et al. The iceLogo web server and SOAP service for determining protein consensus sequences. Nucleic Acids Res 43, W543–546, (2015).
Wu, Y. et al. The Structure of a Peptide-Loaded Shark MHC Class I Molecule Reveals Features of the Binding between β(2)-Microglobulin and H Chain Conserved in Evolution. J. Immunol. 207, 308–321 (2021).
Chu, F. et al. First glimpse of the peptide presentation by rhesus macaque MHC class I: crystal structures of Mamu-A*01 complexed with two immunogenic SIV epitopes and insights into CTL escape. J. Immunol. 178, 944–952 (2007).
Shi, Y., Qi, J., Iwamoto, A. & Gao, G. F. Plasticity of human CD8 alpha alpha binding to peptide-HLA-A*2402. Mol. Immunol. 48, 2198–2202 (2011).
Kumar, P. et al. Structural basis for T cell alloreactivity among three HLA-B14 and HLA-B27 antigens. J. Biol. Chem. 284, 29784–29797 (2009).
Jappe, E. C. et al. Thermostability profiling of MHC-bound peptides: a new dimension in immunopeptidomics and aid for immunotherapy design. Nat. Commun. 11, 6305 (2020).
Garstka, M. A. et al. The first step of peptide selection in antigen presentation by MHC class I molecules. Proc. Natl Acad. Sci. USA 112, 1505–1510 (2015).
Jumper, J. et al. Highly accurate protein structure prediction with AlphaFold. Nature 596, 583–589 (2021).
Evans, R. et al. Protein complex prediction with AlphaFold-Multimer. bioRxiv https://doi.org/10.1101/2021.10.04.463034 (2021).
Shen, Y., Parks, J. M. & Smith, J. C. HLA Class I Supertype Classification Based on Structural Similarity. J. Immunol. 210, 103–114 (2023).
Motmaen, A. et al. Peptide-binding specificity prediction using fine-tuned protein structure prediction networks. Proc. Natl Acad. Sci. USA 120, e2216697120 (2023).
Kumar, S., Stecher, G., Li, M., Knyaz, C. & Tamura, K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol. Biol. Evol. 35, 1547–1549 (2018).
Waterhouse, A. M., Procter, J. B., Martin, D. M., Clamp, M. & Barton, G. J. Jalview Version 2–a multiple sequence alignment editor and analysis workbench. Bioinformatics 25, 1189–1191 (2009).
Lam, K. S. et al. A new type of synthetic peptide library for identifying ligand-binding activity. Nature 354, 82–84 (1991).
Zhang, L. et al. A Newly Recognized Pairing Mechanism of the α- and β-Chains of the Chicken Peptide-MHC Class II Complex. J. Immunol. 204, 1630–1640 (2019).
Liang, R. et al. Major Histocompatibility Complex Class I (FLA-E*01801) Molecular Structure in Domestic Cats Demonstrates Species-Specific Characteristics in Presenting Viral Antigen Peptides. J. Virol. 92, https://doi.org/10.1128/JVI.01631-17 (2018).
Wei, X., Wang, S., Wang, S., Xie, X. & Zhang, N. Structure and Peptidomes of Swine MHC Class I with Long Peptides Reveal the Cross-Species Characteristics of the Novel N-Terminal Extension Presentation Mode. J. Immunol. 208, 480–491 (2022).
Lill, J. R. et al. Minimal Information About an Immuno-Peptidomics Experiment (MIAIPE). Proteomics 18, e1800110 (2018).
Tobita, T. et al. A role for the P1 anchor residue in the thermal stability of MHC class II molecule I-Ab. Immunol. Lett. 85, 47–52 (2003).
Battye, T. G., Kontogiannis, L., Johnson, O., Powell, H. R. & Leslie, A. G. iMOSFLM: a new graphical interface for diffraction-image processing with MOSFLM. Acta Crystallogr D. Biol. Crystallogr 67, 271–281 (2011).
Collaborative Computational Project, N. The CCP4 suite: programs for protein crystallography. Acta Crystallogr D. Biol. Crystallogr 50, 760–763, (1994).
Vagin, A. & Teplyakov, A. Molecular replacement with MOLREP. Acta Crystallogr D. Biol. Crystallogr 66, 22–25 (2010).
McCoy, A. J. Solving structures of protein complexes by molecular replacement with Phaser. Acta Crystallogr D. Biol. Crystallogr 63, 32–41 (2007).
Rossmann, M. G. The molecular replacement method. Acta Crystallogr A 46, 73–82 (1990).
Emsley, P., Lohkamp, B., Scott, W. G. & Cowtan, K. Features and development of Coot. Acta Crystallogr D. Biol. Crystallogr 66, 486–501 (2010).
Murshudov, G. N. et al. REFMAC5 for the refinement of macromolecular crystal structures. Acta Crystallogr D. Biol. Crystallogr 67, 355–367 (2011).
Afonine, P. V. et al. Towards automated crystallographic structure refinement with phenix.refine. Acta Crystallogr D. Biol. Crystallogr 68, 352–367 (2012).
Williams, C. J. et al. MolProbity: More and better reference data for improved all-atom structure validation. Protein Sci. 27, 293–315 (2018).
Adams, P. D. et al. PHENIX: a comprehensive Python-based system for macromolecular structure solution. Acta Crystallogr D. Biol. Crystallogr 66, 213–221 (2010).
The PyMOL Molecular Graphics System, v. Version 2.0 (2023.03.18).
Berman, H. M. et al. The Protein Data Bank. Nucleic Acids Res. 28, 235–242 (2000).
Perez-Riverol, Y. et al. The PRIDE database resources in 2022: a hub for mass spectrometry-based proteomics evidences. Nucleic Acids Res 50, D543–d552 (2022).
Acknowledgements
We thank the staff at the Shanghai Synchrotron Radiation Facility of China for their technical assistance during data collection. This work was supported financially by Major Science and Technology Project of Liaoning Province (2020JH1/10200003). Beijing Science and Technology Project (Z221100006422009). The 2115 Talent Development Program of China Agricultural University.
Author information
Authors and Affiliations
Contributions
N.Z. designed the experiments. S.W. completed all experiments. Z.Q. and X.W. completed part of the crystal structure analysis. L.Z. and S.W predicted the structure of bat MHC-I. L.D. purified the protein. S.W. and N.Z. wrote the manuscript. Funding was provided by N.Z. S.W. and N.Z. analyzed the data and performed statistical analysis.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Inclusion and Ethics
Our research complies with the “Inclusion and Ethics” Statement.
Peer review
Peer review information
Communications Biology thanks Anthony Purcell and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Primary Handling Editors: Ingrid Span and Christina Karlsson Rosenthal. A peer review file is available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Wang, S., Zheng, L., Wei, X. et al. Amino acid insertion in Bat MHC-I enhances complex stability and augments peptide presentation. Commun Biol 7, 586 (2024). https://doi.org/10.1038/s42003-024-06292-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s42003-024-06292-5
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.