Zeolitic imidazolate frameworks activate endosomal Toll-like receptors and potentiate immunogenicity of SARS-CoV-2 spike protein trimer

Nanomaterials offer unique opportunities to engineer immunomodulatory activity. In this work, we report the Toll-like receptor agonist activity of a nanoscale adjuvant zeolitic imidazolate framework–8 (ZIF-8). The accumulation of ZIF-8 in endosomes and the pH-responsive release of its subunits enable selective engagement with endosomal Toll-like receptors, minimizing the risk of off-target activation. The intrinsic adjuvant properties of ZIF-8, along with the efficient delivery and biomimetic presentation of a severe acute respiratory syndrome coronavirus 2 spike protein receptor-binding domain trimer, primed rapid humoral and cell-mediated immunity in a dose-sparing manner. Our study offers insights for next-generation adjuvants that can potentially impact future vaccine development.


INTRODUCTION
Subunit protein antigens are safer and easier to manufacture than traditional whole pathogen vaccines (1-3), but they often require coadministration with an adjuvant to potentiate antigen presentation and produce effective immune outcomes (4,5).Rationally designed vaccine adjuvants aim to trigger specific receptors on antigen-presenting cells (APCs) to direct the immune system to respond in a T helper 1, T helper 2, or T helper 17 manner (6,7).Pattern recognition receptors (PRRs), particularly Toll-like receptors (TLRs), expressed mainly on APCs play a critical role in the modulation of these immune responses (8)(9)(10).PRR agonists help generate greater immune response magnitude and long-lasting memory response against vaccine targets when coadministered with subunit protein antigens (11,12).Nevertheless, only a few PRR adjuvants are licensed for use in U.S. Food and Drug Administration-approved vaccines owing to safety and tolerability concerns.However, these concerns can be potentially mitigated by novel delivery approaches including nanoparticle (NP)-based methods (13,14).
Biomimetic NP delivery strategies have been shown to enhance antibody titers even at ultralow antigen doses (4).A range of NPs have demonstrated adjuvant activity through enhancing the codelivery of antigens and adjuvants to APCs and lymph nodes (LNs) (4,15).NPs allow for passive vaccine drainage to LNs (12,14), create multivalent antigen interactions with B cells, and facilitate antigen presentation onto class I major histocompatibility complex (MHC I) molecules that enhance CD8 + T cell responses (16).Several studies have revealed that polymeric NP antigens (17,18) and liposomes (19) activate stimulatory immune pathways such as TLRs or inflammasomes in the absence of adjuvants and antigens.The intrinsic immune characteristics of NPs should be explored to avoid potential alteration in immune signal processing by APCs that result from the excessive activation of innate immune responses.Investigating NPs' immunostimulatory properties offers an opportunity to study novel carriers that not only tune the spatiotemporal distribution of vaccine antigen but also directly engage with innate signaling pathways and enable selective activation of immunologic pathways.
Zeolitic imidazolate framework-8 (ZIF-8)-a subclass of metal organic frameworks-encapsulates antigen through a biomimetic mineralization process, affording exceptional protection from biological, thermal, and chemical degradation while maintaining bioavailability (20)(21)(22).ZIF-8 is also established as an efficient vaccine delivery platform (23)(24)(25).Preclinical studies demonstrated that ZIF-8 modulates vaccine-induced immune responses and stimulates high antibody responses through the codelivery of antigen and adjuvant to APCs (23,24).However, the contribution of ZIF-8's composition to its adjuvant properties is poorly characterized.Imidazole, the basic building block of ZIF-8, is the smallest moiety in synthetic TLR agonists and can be modified to produce TLRspecific agonist or antagonist molecules, especially TLR-7 and TLR-8 (26,27).Here, we investigated how ZIF-8 NPs can modulate innate immune activation to enhance adaptive immune responses to protein antigens, using severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein receptor-binding domain (RBD) trimer as a model antigen.To better understand how ZIF-8 improves immunogenicity, we used biodistribution studies, immunological techniques, and transcriptome sequencing.
We found that ZIF-8 can passively drain to the draining LN (dLN) while also activating innate immune responses at the injection site.After uptake, ZIF-8 NPs undergo pH-responsive degradation in the endosome, allowing them to simultaneously release their payload, while the imidazole produced by their degradation activates TLRs in endosomes (TLRe).Analysis of gene expression in the dLN suggests a MyD88-dependent response and the activation of nuclear factor κB (NF-κB), leading to the production of interleukin-6 (IL-6) and type I interferon β (IFN-β) associated with antiviral responses.This coincided with increased expression of CCR-7 and costimulatory molecules such as CD80 on TLR-responsive macrophage and dendritic cell (DC) populations.When administered to mice while encapsulating RBD, ZIF-8 allows for dose sparing compared to free RBD and the ability to coencapsulate with a TLR-7 agonist adjuvant, gardiquimod (Gdq).ZIF-8 also improves the stability of RBD in storage.

ZIF-8 encapsulates Gdq and RBD trimer and exhibits pH-responsive release
ZIF-8 was prepared as described previously (22).Transmission electron microscopy (TEM) revealed cubic crystals of ZIF-8 encapsulating Gdq and RBD trimer (GR-ZIF) with an average size of 150 nm that is similar to that of ZIF-8 (Fig. 1, A and B).Likewise, the powder x-ray diffraction (PXRD) pattern of GR-ZIF showed a similar pattern to that of simulated ZIF-8 (Fig. 1C).The Fourier transform infrared (FTIR) spectrum of GR-ZIF spectrum depicted obvious stretching vibrations of C═O at 1606 cm −1 and O─H at 3451 cm −1 , which were attributed to RBD trimer (Fig. 1D).FTIR analysis demonstrated that RBD trimer is not only encapsulated within ZIF-8 but also adsorbed on the surface (Fig. 1D), which allows for the enhanced presentation of RBD trimer.The uniformity and intercalation of Gdq and RBD trimer in ZIF-8 were also validated (text S1 and fig.S1).
The acidic environment of lysosomes and endosomes in APCs is expected to induce the release of RBD trimer and Gdq from ZIF-(Fig.1, E and F).The pH-responsive RBD trimer release from ZIFwas monitored by enzyme-linked immunosorbent assay (ELISA) at physiological (pH 7.4) and lyso/endosomal acidic (pH 5. conditions.Under physiological conditions, less than 3% of RBD trimer was released over 24 hours in solution (Fig. 1E).Screening R-ZIF in a proteinaceous environment [Dulbecco's Modified Eagle Medium (DMEM) supplemented with 10% fetal bovine serum (pH 7. 2)] indicated that 6% of RBD trimer was released from ZIF-8 in 4 hours with no further release observed after 24 hours.However, a total of 43 and 60% of RBD trimer were released at acidic pH after 20 and 120 min, respectively (Fig. 1E).TEM images confirmed that low pH triggers the degradation process for GR-ZIF.Its complete decomposition after 24 hours (Fig. 1, F and G) agrees with previously reported data (22).

ZIF-8 potentiates innate immunity through the biomimetic presentation of RBD trimer
We first assessed ZIF-8 accumulation in dLN and subsequent immune activation using a series of experiments to characterize ZIF-8 biodistribution.A model cargo, XenoLight DiR, was loaded in ZIF-8 (D-ZIF) to assess the temporal localization of ZIF-8 in the dLN.Measuring Zn 2+ in dLNs using inductively coupled plasma mass spectrometry (ICP-MS) confirmed a significant increase (P < 0.01) in Zn 2+ 24 hours after injection (Fig. 2A).Flow cytometry analysis of APCs in dLN demonstrates that levels of CD45 + CD11c + associated with D-ZIF continuously increased over time (Fig. 2B).The accumulation of D-ZIF in CD45 + CD11c + DCs within 1 hour indicates intercellular drainage, while the significant increase in the proportion of CD45 + CD11c + DCs associated with D-ZIF signal at 24 hours indicates intracellular drainage.This trafficking was also observed in pathogens (28) and other adjuvants such as MF59 adjuvant (29) and in agreement with other reports for submicron particles (14,30).The activation of innate immune cells is essential for potent vaccine responses.We, therefore, investigated the capacity of ZIF-8 to activate DCs professional APCs that recognize viral antigens (31,32).Levels of CD11c + CD8 + DC and CD11c + MHC II + DC increased significantly 3 and 24 hours after immunization in dLNs (Fig. 2, D and  E), corroborating in vitro studies with DCs treated with ZIF-8 (text S2 and figs.S3 and S4).Plasmacytoid dendritic cells (pDCs) in particular participate in the first line of defense against viral infections by acting as innate effector cells, which initiate type-I IFN antiviral responses (33).They cooperate with Conventional type 1 dendritic cells (cDC.1) in LN to induce antiviral CD8 + T cells (33).Temporal trafficking of pDC in dLN revealed a significant increase in pDC levels associated with D-ZIF at 3 and 24 hours (Fig. 2C and fig.S5A).While some studies have demonstrated the ability of ZIF-8 to promote antigen and adjuvant delivery (23), the underlining immunological mechanisms are not fully understood.

pH-responsive degradation of ZIF-8 in endosomes releases imidazole subunits that directly engage with the activation of TLRe and induce proinflammatory responses
Imidazole, the subunit building block of ZIF-8, is a moiety found in many commercially available TLR agonists, such as resiquimod and Gdq (34,35).We hypothesized that APCs may be activated by ZIF-8 degradation into imidazole subunits after endocytosis, the dominant uptake mechanism by which APCs take up ZIF-8.The acidic environment of endosomes induces ZIF-8 degradation, resulting in the release of its imidazole subunits and payload.Given the role of DCs in activating immune response (36,37), ZIF-8 internalization and endosomal escape were first assessed in DCs in vitro (text S2 and fig.S3).To analyze the impact of releasing ZIF-8 subunits in endosomes upon its dissociation at acidic environment, we investigated the activation status of TLR-7, TLR-8, and TLR-9 using an in vitro quantitative TLR reporter system.While Gdq is known to selectively activate TLR-7 at any concentration and TLR-8 only at high concentration, ZIF-8 activated TLR-7 at any given concentration and TLR-8 and TLR-9 only at high concentration (text S3 and fig.S6).

ZIF-8 and GR-ZIF activate TLR expression in dLNs
We next assessed the expression of TLR-3, TLR-7, and TLR-9 in dLNs harvested 24 hours after injection from mice immunized with phosphate-buffered saline (PBS) or ZIF-8 alone using reverse transcription quantitative polymerase chain reaction (RT-qPCR).TLR-7 expression increased by four and eight times in dLN of mice immunized with low (25 μg) and high (200 μg) doses of ZIF-8, respectively (Fig. 3A).TLR-3 and TLR-9 expression increased (P < 0.05) only at a high dose of ZIF-8, indicating the importance of ZIF-8 accumulation in endosomes for the cross-activation of other TLRe (Fig. 3, B and C).To validate our results, we evaluated the expression of several target genes involved in the TLRdependent signaling pathway such as MyD88 and NF-κB.Both doses of ZIF-8 significantly (P < 0.001) increased MyD88 and NF-κB expression (P < 0.0001; Fig. 3, D and E).
Further, in studies using vaccine antigens and adjuvants, incorporating Gdq and RBD in ZIF-8 induced greater (P < 0.001) TLR-7 expression (Fig. 3I) and TLR-9 expression (fig.S7A) than RBD or Gdq alone.The activation of TLR-9 in response to R-ZIF and GR-ZIF may suggest synergy associated with combining ZIF-8 and Gdq.
The expression of MyD88 and NF-κB was also higher in GR-ZIFtreated mice compared to a high dose of ZIF-8 (fig.S7, B and C).RBD trimer induced slight increase in the expression of TLR-7 but failed to activate TLR-9, MyD88, and NF-κB (Fig. 3I and fig.S7, A to  C).The activation of multiple TLRe has the potential to drive more robust cellular and humoral responses (38).
MyD88 is an essential adaptor in most TLRs (7).Therefore, we investigated the specificity of TLRe activation in response to ZIF-8 by quantifying the expression of TLR-4.TLR-4 recognizes polysaccharides and should not be susceptible to cross-activation at high concentrations of ZIF-8 such as TLR-3, TLR-7, TLR-8, and TLR-9, which all recognize nucleic acids-further supporting the specificity of the ZIF-8 activation mechanism.The expression of TLR-4 was not affected by ZIF-8 or GR-ZIF (P > 0.05; fig.S7D), confirming the selectivity of ZIF-8 TLR activity.We also assessed gene expression of several cytokines indicative of specific TLR signaling cascades: IFN-βinduced and IL-6.ZIF-8, R-ZIF, and GR-ZIF significantly increased IFN-β (P < 0.0001) and IL-6 (P < 0.05) expression compared to PBS control, Gdq, or RBD trimer (Fig. 3, F, G, and J, and fig.S7E).
NF-κB has been identified as a key regulator of TLR-induced DC maturation.It leads to the up-regulation of CCR7, MHC II, and costimulatory molecules such as CD80 (39).We thus investigated the activation status of leukocytes in response to ZIF-8 using D-ZIF.CCR7 is the chemokine receptor required for the migration of leukocytes from peripheral tissue to LNs.A significant fraction of monocyte-derived DCs (moDCs), cDCs, subcapsular sinus macrophage (SSM), and marginal sinus macrophage (MSM) associated with D-ZIF signal in dLN was migratory (CCR7 + ) (Fig. 4A and fig.S5A).Further stratification of cDC subsets revealed that migratory DCs (m.DCs) accounted for the majority of cDCs in dLN (Fig. 4B).Migratory cDC.2 was significantly increased (P < 0.0001), suggesting that the majority of cDC.2 in dLN was of myeloid origin (Fig. 4C).Compared to PBS controls, D-ZIF induced the up-regulation of CD80 (B7-2) in B cells, SSMs, MSMs, moDCs, and migratory conventional DCs (m.cDC) (fig.S5B).A significant (P < 0.05) increase in CD80 expression was also observed in m.cDC.2 (fig.S5C).No significant expression in CD80 was associated with cDC.1 as CD8 + DCs do not respond to TLR7-mediated stimuli (40), explaining the markedly high expression of CD80 in leukocytes associated with high levels of CCR7.
Methylimidazole (HMIM) is the smallest moiety in synthetic TLR agonists and nitrogenous bases comprise both RNA and DNA (26,27).We therefore used an in vitro quantitative TLR reporter system and in vivo RT-qPCR to probe the connection between HMIM and TLR agonist activity of ZIF-8.HMIM showed a similar

of 11
activity of ZIF-8 in which it activated TLR-7 reporter cells after 24 hours at any given concentration and TLR-8 and TLR-9 reporter cells only at high concentration (fig.S6G).Zn 2+ showed no activation in TLR-7, TLR-8, or TLR-9 reporter cells (fig.S6H), indicating that ZIF-8 immunogenicity is mainly driven by HMIM degradation products.To assess the activation of TLR-7 in vivo, dLNs of mice injected with HMIM were harvested after 24 hours for RT-qPCR analysis.TLR-7 was up-regulated (P < 0.05) in response to HMIM (fig.S7F).HMIM is a small and hydrophilic molecule that is expected to diffuse rapidly to the systemic circulation.This likely contributed to the low TLR-7 expression compared to ZIF-8.Thus, ZIF-8 immunomodulatory properties are, in part, influenced by HMIM, indicating that the chemical building blocks of ZIF-8 play a major role in augmenting the immunogenicity of the carrying antigen (Fig. 3H).
To compare ZIF-8 to another vaccine delivery vehicle, we studied the possible activation of TLRe in response to mRNA-LNP.mRNA-LNP had no significant influence on the expression of TLR-7 (fig.S8C) at the used dose but significantly (P < 0.01) enhanced TLR-9 and NF-κB expression (fig.S8, D and E).Note that mRNA loaded in LNP codes for the full spike protein, which is known to contribute to the self-adjuvant properties of mRNA vaccines.Unlike ZIF-8, mRNA-LNP immunogenicity is mainly driven by moDCs and m.cDC.1 (text S4 and fig.S8).This suggests that the TLR-7-responsive properties of ZIF-8 may be somewhat unique among potential nanoscale vaccine delivery vehicles.
To assess the translatability of GR-ZIF, we examined its safety and stability.GR-ZIF did not induce reactive oxygen species production nor apoptosis at the given concentration in mice, and the body well accommodated the increase in Zn 2+ (text S5 and fig.S9).ZIF-8 improved RBD stability in both dry and solution states at room temperature and 4°C over 2 months (fig.S10 and text S6).The facile synthesis of GR-ZIF and improved stability of RBD trimer could potentially enable easier technology transfer to low-income countries and other areas with limited infrastructure.

ZIF-8 encapsulating RBD trimer promotes humoral and cellular immune responses
We hypothesized that the immunostimulatory properties of ZIF-8 NPs would make them a useful platform for subunit vaccines.First, we assessed the capacity of RBD trimer, R-ZIF, and GR-ZIF to stimulate humoral responses.Brachial dLNs were harvested 2 weeks following injection to analyze germinal center B cell (B GC ) responses.R-ZIF and GR-ZIF elicited B GC responses greater than RBD trimer (Fig. 5, B and  C, and fig.S11), suggesting the substantial impact of ZIF-8 in potentiating subunit protein immunogenicity.
We then vaccinated C57BL/6 mice with three injections of either RBD trimer (3 μg), R-ZIF, or GR-ZIF (containing 0.9 μg of RBD trimer) at week 0, 3, and 6 (Fig. 5A).To investigate RBDspecific immunoglobulin G (IgG), serum was collected 2 weeks after the first booster doses and 3 weeks after the second booster dose as illustrated in Fig. 5A.Median RBD-specific total IgG responses following both booster doses were high in R-ZIF and GR-ZIF but low in RBD trimer (Fig. 5D).Unlike RBD trimer, R-ZIF and GR-ZIF elicited balanced, high-titer, IgG1, and IgG2a serum responses (Fig. 5, E and F).While the incorporation of Gdq induced TLR-7 and IFN-β expression as indicated in Fig. 3 (I and J), R-ZIF and GR-ZIF exhibited comparable B GC responses and yielded comparable IgG1 and IgG2a titers after the first booster dose (Fig. 5, B to F). Together, these results indicate that the use of ZIF-8 in a vaccine formulation enhances antigen immunogenicity.
SARS-CoV-2-specific cellular responses were assessed by enzymelinked immunospot (ELISpot).Three weeks following the second booster dose, splenocytes were restimulated ex vivo with RBD trimer to evaluate IFN-γ-producing cells.Mice treated with GR-ZIF showed eightfold increase in the number of IFN-γ counts on average compared with PBS controls (P < 0.0001; Fig. 5F).In contrast, animals immunized with RBD trimer failed to induce IFN-γ release and, by extension, antigen-specific T cell responses (Fig. 5G), suggesting the impact of using ZIF-8 in improving memory responses.

Transcriptome analysis of dLNs in vaccinated mice reveals that GR-ZIF enhanced the activation of innate signaling pathways
To understand differences between transcriptional responses of RBD trimer and GR-ZIF and to isolate the role of ZIF-8 in tuning immune responses, RNA sequencing (RNA-seq) was performed.dLNs were collected at the end of the first week after the first booster dose for analysis.Since GR-ZIF showed rapid induction of adaptive immune responses and R-ZIF yielded B GC responses and titers similar to GR-ZIF after the first booster dose, we sought to compare transcriptional responses of GR-ZIF to RBD trimer.Principal components analysis of the transcriptomes for all vaccines revealed grouping of responses away from the PBS control (fig.S12A).Further RNA-seq global transcriptome analysis identified immune-related transcriptomic pathways affected by the immunization.The gene set enrichment analysis (GSEA) (41) showed significant up-regulation of 17 pathways in GR-ZIF-immunized mice.These pathways were related to antigen presentation on APC, T cell activation, and IFN-γ production (fig.S12B).
Unexpectedly, our differential GSEA revealed that TLR-7 was still enriched after a week of immunizing mice with ZIF-8 [false discovery rate (FDR) < 0.05; Fig. 6A and fig.S12A].Loading RBD trimer and Gdq in ZIF-8 induced better enrichment of TLR-7 (FDR < 0.05; Fig. 6A) but not in soluble RBD trimer (FDR > 0.05; Fig. 6A).To validate data obtained from the RNA-seq analysis, we performed RT-qPCR on selected genes TLR-7 and TLR-9.The expression of TLR-7 and TLR-9 for R-ZIF-and GR-ZIF-treated mice significantly increased (fig.S12, B and C), while no significant increase was observed in RBD trimer-immunized mice.Together, these data suggest an intrinsic adjuvant effect of the ZIF-8 carrier in activating innate signaling pathways, which can be further enhanced with codelivery of a TLR agonist.We next assessed the differentially expressed genes in the dLNs of mice given GR-ZIF or RBD trimer and compared to those injected with PBS.Volcano plot analysis identified significant up-regulation of 55 genes related to innate, B and T cell genes in GR-ZIF-immunized mice (Fig. 6B).Most of the genes were specifically up-regulated in GR-ZIF versus PBS but not RBD trimer versus PBS.Some of the most highly differentially expressed genes included activation-induced cytidine deaminase (AICDA) (mediates antibody class switching and somatic hypermutation in B cells) (42), germinal center associated signaling and motility (GCSAM, regulator of B cell receptor signaling) (43), secretory leukocyte peptidase inhibitor (SLPI, plays a role in regulating the activation of NF-κB and inflammatory responses) (44), and T cell receptor alpha variable 16 (TRAV16, αβ T cell receptors recognize peptide-MHCs displayed by APC) (45,46).To confirm whether R-ZIF induces a similar expression profile as GR-ZIF, AICDA and TRAV16 expression was analyzed using RT-qPCR.The RT-qPCR analysis validated the expression of AICDA and TRAV16 in the dLNs of R-ZIF-treated mice (fig.S14).Together, these data suggest an intrinsic adjuvant effect of the ZIF-8 carrier in activating innate signaling pathways, which can be further enhanced with codelivery of a TLR agonist.

A B
Spleen harvest for IFN ELISpot Serum collection Immunization overlooking their potential in activating stimulatory immune pathways that could alter their immunostimulatory properties.Therefore, it is important to understand the intrinsic immune characteristics of NPs to ensure moderate activation of inflammatory responses and enable selective activation of immunologic pathways to maximize vaccine potency.Here, we investigated the influence of ZIF-8 accumulation, distribution, retention, and efficacy on immune response, which are important determinants for successful clinical translation.We demonstrated that the chemical and physical properties of ZIF-8 modulate immune responses to SARS-CoV-2 RBD trimer at the cellular and molecular levels.Simultaneously, ZIF-8 physical properties facilitated multivalent antigen display and LN drainage in a pathogenmimetic approach.This was also observed in other U.S. Food and Drug Administration-approved NP-based vaccine platforms used to deliver SARS-CoV-2 spike protein, such as micelle-attached spike (NVX-CoV 2373) (5).
We also demonstrated that the specific chemical building blocks of ZIF-8 play a major role in augmenting the immunogenicity of the encapsulated antigen.The pH-responsive degradation of ZIF-8 into HMIM selectively engages TLR in APCs.This, in turn, activates MyD88-dependent pathways that activate NF-κB, leading to the production of IL-6 and IFN-β associated with antiviral responses and increased surface expression of CCR-7 and CD80 primarily on TLR-responsive APCs.This suggests that ZIF-8 NPs are an immunostimulatory vehicle capable of selectively activating TLRe upon cargo release.Both physical and chemical properties resulted in RBD trimer dose-sparing and efficient innate and adaptive immune responses.Future studies should address the potential impact of cellular and molecular networks activated by ZIF-8 on virus neutralizing antibody responses.Overall, our study demonstrates that probing cellular and molecular networks that drive immune responses to NPs' composition and physical properties can reveal mechanistic insights and can potentially enhance the rational development of novel vaccines.

GR-ZIF preparation and characterization
GR-ZIF was synthesized by stirring Gdq (0.1 mg ml −1 ), RBD trimer (0.1 mg ml −1 ), and HMIM (2.5 M, 0.9 ml) for 20 min.Zinc nitrate solution (0.5 M, 0.1 ml) was slowly added under mechanical agitation for 20 min.The resulting product was collected by centrifugation and washed three times with deionized water to remove any residues.ZIF-8 was synthesized by slowly adding zinc nitrate solution (0.5 M, 0.1 ml) to HMIM (2.5 M, 0.9 ml).The solution was agitated for 20 min.The supernatant of GR-ZIF was collected to calculate the loading capacity (LC) and loading efficiency (LE) of RBD trimer by ELISA and Gdq by ultraviolet-visible.D-ZIF was synthesized by stirring DiR (0.25 mg) with HMIM, followed by the slow addition of zinc nitrate solution.LE and LC of Gdq and RBD trimer in ZIF-8 were calculated as follows The zeta potential of GR-ZIF was performed using a Malvern Zetasizer Nano ZS at 25°C at pH 7.3 in aqueous solutions.PXRD measurements were performed using a Panalytical X'Pert Pro X-ray powder diffractometer at the MIT Materials Research Laboratory using the Cu Kα radiation (40 V, 40 mA, λ = 1.54056Å) in a θ-θ mode from 20° to 90° (2θ).TEM images were obtained using FEI Tecnai microscope operating at 120 kV at the MIT Materials Research Laboratory.For visualization by TEM, samples were prepared by dropping the solution on a copper grid of 300 mesh (Electron Microscopy Sciences, LC 300-Cu).To evaluate the release of RBD trimer from ZIF-8, SARS-CoV-2 spike RBD ELISA kit (Acro Biosystems) was performed.Aliquots of hydrochloric acid were added to R-ZIF (0.6 mg ml −1 ) in PBS to reach pH 5.3 at 37°C.PBS only or DMEM supplemented with 10% serum were added to the sample of pH 7.3.The supernatant of the mixture solution was obtained through centrifugation at different time points.
In vitro ZIF-8 biocompatibility CCK-8 assay was performed according to the manufacturer's protocol.Bone marrow-derived dendritic cells (BMDCs) were provided by the Koch Institute High Throughput Sciences Core.Briefly, BMDC (5 × 10 3 cells per well) were seeded onto a 96-well plate.Cells were incubated with different concentrations (100, 50, 25, 12, 6, and 3 μg ml −1 ) of ZIF-8 in 200 μl of DMEM at 37°C for 24 hours.The medium was then discarded, and the prepared culture medium containing 10% CCK-8 solution was added into each well, including a negative control of culture medium alone.After 3 hours of incubation, the absorbance was measured at 450 nm using a microplate reader (Tecan Infinite 200 microplate reader).

In vitro BMDC activation
BMDCs were seeded into a 24-well plate at 2 × 10 5 cells per well and treated with PBS, soluble RBD trimer, or GR-ZIF at ZIF-8 (10 μg ml −1 ) for 24 hours.Cells were then collected and measured for the expression of CD80 (BioLegend, 104714-BL) and CD86 (BioLegend, 105106-BL) by flow cytometry (BD FACSCanto, USA).Cell culture supernatant was also collected after 48 hours for the measurement of DC secretion of tumor necrosis factor-α by ELISA.

Uptake and intracellular processing of RBD trimer by DCs
For RBD trimer uptake study, DCs were seeded into a 24-well plate at 2 × 10 5 cells per well and incubated with soluble RBD trimer or GR-ZIF at 10 μg of ZIF-8 per well for 3 or 24 hours, followed by the measurement of fluorescence intensity among CD11c + DCs (BioLegend, 117328) by flow cytometry.For the confocal laser scanning microscopy (Olympus FV1200 CLSM) study, DCs were seeded on glass coverslips put into a 24-well plate and incubated with the soluble mixture of Alexa Fluor 647-bovine serum albumin (BSA) or BSA loaded ZIF-8) for 3 or 24 hours, followed by washing and staining with 0.1 μM LysoTracker and Hoechst 33342 (1 μg ml −1 ) at 37°C for 1 hour.Cells were then fixed using 4% paraformaldehyde solution, coated on a glass slide using the ProLong Diamond Antifade Mountant, and visualized using confocal laser scanning microscopy.

Ethics statement
All animal procedures were approved by the Institutional Animal Care and Use Committees.Six-to 8-week-old C57BL/6 female mice were purchased from Charles River Laboratories Inc. Methods were chosen to minimize pain and distress to the mice, which were observed daily by trained animal care staff.The mice were terminated using a CO 2 inhalation chamber.

In vivo time-course for LN drainage analysis
For GR-ZIF accumulation studies, C57BL/6 female mice were intradermally injected with XenoLight DiR (excitation/emission, 710/760 nm) or D-ZIF-8 in the right forelimb of hairless mice.For flow cytometry analysis, the brachial dLN was harvested and dissociated using a spleen dissociation medium according to the manufacturer's instructions.Cells were centrifuged and resuspended in a buffer containing zombie violet, CD45, CD11c (Tonbo Biosciences, 910004), CD11b (Tonbo Biosciences, 910004), CD8 (Tonbo Biosciences, 910004), and MHC II (Tonbo Biosciences, 910004).Cells were then washed and resuspended in PBS + 2% BSA.APC staining for flow cytometry Intradermal administration of D-ZIF was performed in the right forelimb of C57BL/6 mice using a tuberculin syringe.The brachial dLN was harvested and dissociated as described previously (55).The resulting single-cell suspension was stained with Live/Dead near-infrared viability stain (Thermo Fisher Scientific, Waltham, MA, USA), and fluorochrome-conjugated antibodies were adapted from Optimized Multicolor Immunofluorescence Panel 061 (OMIP 061) (55)  In vivo pharmacokinetics and histological assays C57BL/6 female mice were intradermally injected GR-ZIF-8 in the right forelimb of hairless mice.Blood, LN, kidneys, and spleen were all harvested and digested using aqua regia and incubated at 80°C for 2 hours.The obtained liquid was diluted and subjected to ICP-MS analysis at the MIT Materials Research Laboratory to measure Zn 2+ .In the histological assay, LNs were snap-frozen inside cryomolds using Tissue-Tek O.C.T. Compound mounting medium (Thermo Fisher Scientific, 23730571, Pittsburgh, PA).

Immunization
C57BL/6 female mice received primary immunization and two boosts on weeks 0, 3, and 6.Mice in each group (n = 4) were immunized intradermally at the right forelimb with PBS, RBD trimer, R-ZIF, or GR-ZIF.Mice were bled on weeks 2, 5, and 9.

Electrochemiluminescence assay for total IgG determination
Electrochemiluminescence assay plates [MesoScale Discovery (MSD) SARS-CoV-2 IgG catalog no.N05CA-1; panel 33] were designed and produced with one antigen spot in each well.Briefly, the antigens included WA1/2020 variant RBD protein.Plates were blocked with 50 μl of blocker A (1% BSA in distilled water) solution for 30 min at room temperature shaking at 700 rpm with a digital microplate shaker.Meanwhile, the serum was diluted 1:5000 in Diluent 100 (MesoScale Discovery).The plates were three times washed with 150 μl of wash buffer (0.5% Tween 20 in 1× PBS) and blotted dry.Then, 50 μl of the diluted samples were added in duplicate to the plates and set to shake at 700 rpm at room temperature for 2 hours.Secondary antibody was prepared using Jackson ImmunoResearch Rabbit Anti-Mouse IgG detection antibody (catalog no.315-005-045) conjugated to the MSD GOLD SULFO-TAG by N-hydroxysuccinimide ester chemistry per the manufacturer's guidelines (catalog no.R91AO-1).The plates were washed three times, and 50 μl of tagged secondary antibody solution diluted 1:1000 in Diluent 100 was added to each well and incubated shaking at 700 rpm at room temperature for at least 1 hour.Plates were washed again three times, 150 μl of MSD GOLD Read Buffer B was added to each well, and the plates were read immediately after on a MESO QuickPlex SQ 120 machine.Antibody titer for each sample was reported as relative light units (RLU) that were calculated as average sample RLU minus blank RLU for each sample.The limit of detection was defined as 1000 RLU for each assay.

SARS-CoV-2 anti-RBD binding titers
Plates are coated with purified recombinant SARS-CoV-2 spike S1-RBD antigen (2 μg ml −1 ) in 1× PBS and incubated at 4°C overnight.After incubation, plates were washed once with a wash buffer (0.05% Tween 20 in 1× PBS) and blocked with 350 μl of BSA block per well.The block solution was discarded after 2 to 3 hours of incubation at room temperature, and plates were blotted dry.Plates were incubated with serial dilutions of sera.Goat anti-mouse IgG1 horseradish peroxidase conjugates (ab97230) or goat anti-mouse IgG2a horseradish peroxidase conjugates (ab97245) were used as secondary antibodies, and 3,5,3′,5′-tetramethylbenzidine was used as a substrate.End-point titers were calculated as the dilution that emitted an optical density exceeding 3× background produced by serum from naïve mice.Antibody titers were log-transformed before all statistical processes and are presented end-point titer.

ELISpot assay
Splenocytes from immunized animals were isolated three months following prime dose and assessed by ELISpot (R&D Systems, EL485) for IFN-γ production according to the manufacturer's protocol.Splenocytes were added at a final concentration of 0.5 × 10 6 per well and coincubated with recall antigens (0.1 μg per well) in complete medium for 16 hours at 37°C.Following incubation, plates were washed with the provided washing buffer, and anti-mouse IFN-γ-biotin was added overnight at 4°C.Plates were washed, and streptavidin was added for a further 2 hours at room temperature.Plates were washed with washing solution, and spots were developed using bromochloroindolyl phosphate-nitro blue tetrazolium substrate.Once plates were dry, spots were counted using ELISpot reader system.

RNA isolation and total RNA-seq preparation
Mice were injected with RBD trimer (3 μg), ZIF-8, or GR-ZIF (containing 0.9 μg of RBD trimer) at week 0, 3, and 6.RNA was extracted from harvested brachial dLN (three biological replicates per condition) at week 5, after the first booster dose.Total RNA was isolated using RNA MiniPrep Kit (QIAGEN).Libraries were prepared using the ribosomal RNA depletion chemistry method, and whole-transcriptome sequencing was performed using NextSeq by Illumina.RNA-seq data were used to quantify transcripts from the mm10 mouse assembly with the Ensembl version 101 annotation using Salmon version 1.3.0(56).Gene level summaries were prepared using tximport version 1.20.0(57) running under R version 4.1 (58).Differential expression analysis was done with DESeq2 version 1.32.0 (59,60), and differentially expressed genes were defined as those having an absolute apeglm (61) log 2 fold change greater than 1 and an adjusted P value less than 0.05.Data parsing and clustering were done using Tibco Spotfire Analyst 7.6.1.Mouse genes were mapped to human orthologs using Mouse Genome Informatics (www.informatics.jax.org/)orthology report, and preranked GSEA (62) was done using javaGSEA version 2.3.0_beta_2 with msigDb version 6.2 gene sets (41).RNA-seq data are available from the Gene Expression Omnibus under accession number GSE246178.

Quantitative real-time polymerase chain reaction
Mice were injected with RBD trimer (3 μg), ZIF-8, Gdq (0.8 μg), R-ZIF or GR-ZIF (containing 0.9 μg of RBD trimer), or mRNA-LNP (containing 1 μg of mRNA) at weeks 0 and 3. RNA was extracted from harvested brachial dLN at week 5, after the first booster dose.Total RNA was isolated using RNA MiniPrep Kit (QIAGEN).cDNA was prepared with 1 μg of total using cDNA synthesis kit (QIAGEN).qPCR was performed using Roche real-time machine using cDNA template, primer mix, and SYBR green (Bio-Rad) in a final volume of 20 μl of PCR reaction.PCR primer pairs are listed in table S2.

Supplementary Materials
This PDF file includes: texts S1 to S6 Figs.S1 to S14 tables S1 and S2 References

Fig. 1 .
Fig. 1.ZIF-8 encapsulate RBD trimer and Gdq.(A and B) teM micrographs and simplified schematic of ZiF-8 (A) and GR-ZiF (B).(C) PXRd patterns of simulated ZiF-8 and GR-ZiF, demonstrating the retained crystallinity of ZiF-8 after loading RBd trimer and Gdq.(D) FtiR analyses of RBd trimer and GR-ZiF, indicating the adsorption of RBd trimer on the surface of ZiF-8.a.u., arbitrary units.(E) cumulative release profile of RBd trimer from GR-ZiF at physiological ph (7) and endosomal ph (5.5), demonstrating the ph responsive release from ZiF-8.(F and G) teM micrographs of GR-ZiF degradation over time in response to low ph (5.3).

Fig. 2 .
Fig. 2. RBD trimer presentation is promoted through the rapid GR-ZIF drainage to dLN by a cell-independent and -dependent mechanisms.(A) c57Bl/6 mice were intradermally injected with either soluble dye (Xenolight diR) or d-ZiF.icP-MS analysis of Zn 2+ in dln isolated 3 and 24 hours after intradermal injection with soluble dye or d-ZiF showing the increase in Zn 2+ in dln in d-ZiF-treated mice.ppb, parts per billion.(B to E) Flow cytometry quantification of cd45 + cd11c + (B), pdcs (c), cdc + Mhc ii + (d), and cdc + cd8 + (e) cells in dlns 1, 3, and 24 hours after injection showing the time-dependent increase in the influx of different subsets of dcs.data are shown as means ± SeM. n = 3 or 4. Statistical significance was calculated by unpaired two-tailed Student's t test or one-way analysis of variance (AnOvA): *P < 0.05, **P < 0.01, ***P < 0.001, and ****P < 0.0001.

Fig. 5 .
Fig. 5. ZIF-8 induces potent cellular and humoral immune responses at a modest dose of 0.9 μg of RBD trimer in vivo.(A) the timeline of vaccine immunization and serum collection.c57Bl/6 mice were immunized intradermally with 3 μg of soluble RBd trimer or 0.9 μg of R-ZiF or GR-ZiF and boosted on weeks 3 and 6 with the same formulations.the levels of RBd-specific igG in the sera collected on weeks 2, 5, and 9 after immunization were detected with eliSA.dlns were harvested for B Gc assessment by flow cytometry on week 2 following prime dose.Spleens were harvested on week 12 after immunization to assess memory responses by iFnγ + eliSpot.(B and C) Quantification (B) and representative flow plots (c) of B Gc cells (B220 + , Gl7 hi , and cd38 lo ) in dln.Fitc, fluorescein isothiocyanate.(D) WA1/2020 RBd-specific igG calculated after logarithmic transformation (log 10 ) of antibody levels.(E and F) RBd-specific igG1 (e) and igG2 (F) as quantified by eliSA.(G) eliSpot analysis of iFnγ spotforming cells in splenocytes after ex vivo stimulation to assess memory responses.data are shown as means ± SeM. n = 3 or 4. Statistical significance was calculated by unpaired two-tailed Student's t test or one-way AnOvA: *P < 0.05, **P < 0.01, and ***P < 0.001.

AFig. 6 .
Fig. 6.RNA-seq analyses of dLNs of C57BL/6 mice immunized intradermally with 3 μg of RBD trimer in solution, ZIF-8, or 0.9 μg of RBD trimer displayed on ZIF-8 with Gdq.One week following first booster dose, dlns were harvested for RnA-seq.(A) GSeA plots showing enrichment of tlR-7-related gene sets in RBd trimer, ZiF-8, and GR-ZiF showing the up regulation of processes related to tlR7 in ZiF-8-and GR-ZiF-treated mice.(B) volcano plot showing variance in gene expression with respect to fold change (Fc) and significance (P < 0.05) in dln of RBd trimer-and GR-ZiF-treated mice.each dot represents an individual gene: Red dots represent the up-regulated genes in GR-ZiF.
LC = [mass loaded drug ∕ mass of loaded drug + NPs] × 100 (1) LE = [mass of drug loaded ∕ mass of initial drug] × 100 (2) FACSymphony A3 (Becton Dickinson) at the MIT Koch Institute Swanson Biotechnology Center Flow Cytometry Facility.After gating on live single cells and lymphocytes, populations were identified according to the gating strategy and immunophenotypes listed in fig.S7 and table S1.