- Faculty of Life Sciences and Biotechnology, Shahid Beheshti University, Tehran, Iran
Since the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) outbreak, several solutions have been proposed to manage the disease. The most viable option for controlling this virus is to produce effective vaccines. Most of the current SARS-CoV-2 vaccines have focused on the infusion spike protein. Spike exists as a trimer and plays a vital role in infecting host cells by binding to the Angiotensin-Converting Enzyme 2 (ACE2) receptor through its Receptor Binding Domain (RBD). Ferritin protein, a naturally occurring iron-storage protein, has gained attention for vaccine production due to its self-assembling property, non-toxic nature, and biocompatibility. Ferritin nanocages have recently been employed in the development of a SARS-CoV-2 vaccination eliciting not only long-term protective memory cells but also a sustained antibody response. In this study, a combination of in silico investigations including molecular docking, molecular dynamics simulations, and immune simulations were carried out to computationally model the monomeric spike protein on the ferritin nanocage as well as to evaluate its stability and interactions for the first time. The structural dynamics of the modeled complex demonstrated noticeable stability. In particular, the Receptor Binding Domain (RBD) and ferritin within the monomeric spike-ferritin complex illustrated significant stability. The lack of alterations in the secondary structure further supported the overall steadiness of the complex. The decline in the distance between ferritin and spike suggests a strong interaction over time. The cross-correlation matrices revealed that the monomeric spike and ferritin move towards each other supporting the stable interaction between spike and ferritin. Further, the orientation of monomeric spike protein within the ferritin unit facilitated the exposure of critical epitopes, specifically upward active Receptor Binding Domain (RBD), enabling effective interactions with the ACE2 receptor. The immune simulations of the model indicated high-level stimulations of both cellular and humoral immunity in the human body. It was also found that the employed model is effective regardless of the mutated spikes in different variants. These findings shed light on the current status of the SARS-CoV-2-ferritin nanoparticle vaccines and could be used as a framework for other similar vaccine designs.
Introduction
More than 773 million people have been affected, and over six million have died due to the global coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (Acuti Martellucci et al., 2020; World Health Organization, 2024; Wu et al., 2020). SARS-CoV-2 is a single-stranded RNA virus with a positive sense and large RNA genome, four structural proteins, 16 nonstructural proteins, and 9–11 accessory proteins (Alexandersen et al., 2020; Malik, 2020). The four structural proteins are spike, envelope, membrane, and nucleocapsid proteins, with the spike glycoprotein (S protein) being particularly important as it is a prominent coronavirus vaccine target (Andersen et al., 2020; Choi et al., 2021; Du et al., 2009; Hu et al., 2021; Walls et al., 2020; Wrapp et al., 2020). SARS-CoV-2 enters cells by binding to the host cellular receptor Angiotensin-Converting Enzyme 2 (ACE2) via its spike protein (Song et al., 2018; Zhao et al., 2020; Piplani et al., 2021; Jackson et al., 2022), and typically causes a lower respiratory tract infection, which can progress to severe acute respiratory syndrome and potentially multiple organ failure (Ganji et al., 2020; Tu et al., 2020; Tao et al., 2021; Yang et al., 2021). The S protein monomer, initially produced as a single polypeptide chain, is made up of a fusion peptide, two heptad repeats, an intracellular domain, an N-terminal domain, two subdomains, and a transmembrane region (Figure 1). After translation it is cleaved into S1 and S2 subunits (Li et al., 2005; Gui et al., 2017). A 25 KDa Receptor Binding Domain (RBD; residue range 334–527) within S1 interacts with the ACE2 receptor (Hoffmann et al., 2020; Walls et al., 2020). As a result, the S protein dictates SARS-COV-2’s infectivity and transmissibility (Du et al., 2009; Yang et al., 2020).
Figure 1. Structural representation of the monomeric spike SARS-CoV-2 ectodomain. (A) monomeric spike in gray cartoons and transparent surfaces. Its different parts, including Receptor Binding Domain (RBD), N-terminal domain (NTD), Fusion Peptide (FP), and Heptad Repeat 1 (HR1) are shown in red, yellow, green, and blue cartoon and transparent surfaces respectively. (B, C) subunits S1 and S2 in cyan and light pink cartoon and transparent surfaces, respectively.
Developing a safe and effective vaccine for SARS-CoV-2 and plausible emerging variants is a top priority for public health, and long-term pandemic control will necessitate one or more effective vaccines that can be widely distributed worldwide (Chen W.-H. et al., 2020; Haiou et al., 2022). Protein subunit vaccines and genetically encoded nucleic acid vaccines are the most effective for the prevention and treatment of SARS-CoV-2 disease (Amanat and Krammer, 2020). Protein nanoparticles which have been vastly used in therapeutic purposes possess a high surface area to volume ratio, enhancing their drug-holding capacity, the solubility of the drug, and their bioavailability (Bhushan et al., 2014).
Ferritin, a family of protein cages that plays a crucial function in iron storage, is highly evolutionarily ubiquitous (Andrews et al., 1992; Uchida et al., 2010; Jutz et al., 2015). Ferritin is a 24-meric cage with octahedral symmetry and a protein nanoparticle (Lawson et al., 1991; Zhang and Orner, 2011; Chakraborti and Chakrabarti, 2019). Ferritin’s surfaces are open to various types of changes, which is one of the reasons behind its utility for biological applications (Jin et al., 2019), including vaccine development. Ferritin nanocages have been utilized to display antigens from various pathogens (Bachmayer et al., 1976; Zhang and Orner, 2011; Kanekiyo et al., 2013; Stone et al., 2013; Kanekiyo et al., 2015; He et al., 2016; Chattopadhyay et al., 2017; Qi et al., 2018; Jin et al., 2019; Kamp et al., 2020; Swanson et al., 2020; Yao et al., 2020; Wuertz et al., 2021), all of which elicit a robust immune response. Nanoparticles in vaccine formulations are predicted to improve antigen stability and immunogenicity, as well as targeted distribution and gradual release (Oyewumi et al., 2010; Carmen et al., 2021; Gordon et al., 2022). Various types of nanoparticles have been utilized in vaccine development (Glück et al., 2004; Borges et al., 2008; Champion et al., 2009; Nochi et al., 2010; Zhu et al., 2010; Wang et al., 2011; Zhao et al., 2012; Bhushan et al., 2014; Zhao et al., 2014; Sliepen et al., 2015; Edwardson and Hilvert, 2019; Swanson et al., 2020). Synthetic vaccines are not only safer than attenuated or inactivated microorganisms, but they also allow for custom vaccine development (Liu et al., 2014). To produce the SARS-CoV-2 vaccine, the spike protein is displayed as an immune system stimulator on ferritin units. The ferritin-base vaccines elicited a long-lasting antibody response and a remarkable long-term memory demonstrating the efficacy of ferritin nanoparticles in preserving antibody response (Kalathiya et al., 2021; Powell et al., 2021; Wang et al., 2021). Recent research indicated that when the RBD and ferritin nanoparticles were stitched together, the resulting synthetic chimeras displayed the RBD effectively and elicited moderate to high-effect immune responses in experimental animal models (Kanekiyo et al., 2013; Qi et al., 2018; Petersen et al., 2020; Yao et al., 2020; Powell et al., 2021). In silico studies of linker and glycosylation design to improve nanoparticle vaccine constructs have been recently investigated (Kalathiya et al., 2021; Masoomi Nomandan et al., 2022). Specifically, the in silico ferritin-SARS-CoV-2 glyco-RBD nanoparticle vaccine has demonstrated significant efficacy against the virus (Masoomi Nomandan et al., 2022). These studies predominantly focused on the RBD of the spike protein (Amanat et al., 2020; Chen et al., 2020; Mulligan et al., 2020). However, experimental evidence has suggested that the potential stability of a full-length spike attached to ferritin seems to be higher than that of the RBD attached to ferritin (Powell et al., 2021). In this study, the attachment of monomeric spike to ferritin (Figure 2) was explored. A combination of molecular docking, molecular dynamics simulations, and immune simulations were carried out to assess the resulting designed vaccine in terms of structural stability and interactions.
Figure 2. Monomeric spike SARS-CoV-2 - Ferritin nanocage vaccine design. The interaction between the Receptor Binding Domain (RBD) of the monomeric spike protein and the Angiotensin-Converting Enzyme 2 (ACE2) receptor found in human cells is illustrated in the left panel. The box represents the 3D structural presentation of one monomeric spike attached to a unit of ferritin. Spike and ferritin are presented with purple and yellow cartoon and transparent surfaces, respectively.
Computational methods
3D structures of monomeric spike, ferritin and lectin
The 3D structure of the SARS-CoV-2 spike ectodomain in an open state, represented by PDB ID 6VYB (Walls et al., 2020), was chosen for the monomeric spike protein. This structure was obtained through electron microscopy, with a resolution of 3.20 Å, and was released on 2020–03-11. The selection of 6VYB was based on its representation of the spike protein in a conformation relevant to its interaction with host cells, which is crucial for eliciting an immune response against the virus antigen. The crystal structure of human L-ferritin, represented by PDB ID 2FG8 (Wang et al., 2006), was selected for the ferritin protein. This structure was determined using X-ray diffraction, with a resolution of 2.50 Å, and was released on 2006–07-04. PDB ID 2FG8 was chosen for its well-defined atomic details, which is crucial for accurate structural analysis of ferritin in the designed vaccine. The lectin protein structure was obtained from PDB ID 5B1W (Nagae et al., 2016). This structure was determined using X-ray diffraction, with a resolution of 3.05 Å, and was released on 2016–05-11.
Molecular docking of monomeric spike and ferritin
The selected structure of the monomeric spike was used to dock to one unit of ferritin with the HADDOCK server (van Zundert et al., 2016), where the HADDOCK server performed rigid body energy minimization in the initial step of the docking process. Residues at the C-terminal of monomeric spike (residues: 1,136–1,147) were chosen as a binding interface area according to previous studies (Powell et al., 2021; Wuertz et al., 2021; Gordon et al., 2022) and several docking scores (Supplementary Figure S1). One exposed loop on the ferritin structure (residues:73–91) was selected as the most plausible binding site (Masoomi Nomandan et al., 2022), which was confirmed by comparing the energetic terms and clustering of the resulting docked poses. The docked pose with the minimum HADDOCK score and lowest RMSD values was considered for model development. The poses were validated via further molecular docking and processed by the ClusPro web server (Supplementary Table S1) (Kozakov et al., 2017). Details outlining the atomistic interactions of the modeled construct are provided (Supplementary Table S2; Supplementary Figure S2). Biomolecular visualization of the complex was carried out by PyMOL (LLC, 2015).
Molecular docking of monomeric spike-ferritin complex and lectin
Lectins are a unique class of proteins or glycoproteins known for their ability to selectively recognize carbohydrate structures and form reversible linkages upon interacting with glycoconjugate complexes (Raposo et al., 2021). The molecular docking analysis examined the interaction between the monomeric spike-ferritin complex and lectin molecules, aiming to explore the potential interference of lectin with vaccine efficacy. HADDOCK server (van Zundert et al., 2016) was employed to dock full-length human lectin and monomeric spike-ferritin nanoparticle vaccine together (Supplementary Figure S3). It was observed that seven residues involved in the ferritin-lectin binding interface were also found to participate in the ferritin-monomeric spike interface. However, the most stable cluster exhibited a Root Mean Squared Deviation (RMSD) value of 4.7±0.0 Å and a HADDOCK score of 163.2±41.3. This suggests that the connection between lectin and the constructed vaccine is weak and unstable, and might not significantly impact the vaccine efficacy. These observations were further validated through ferritin and lectin docking (Supplementary Figure S4). The most stable cluster displayed an RMSD value of 20.1±0.2 Å and a HADDOCK score of 188.3±23.4. Based on these results, it is unlikely that ferritin or the monomeric spike would be bound to lectin proteins in a manner that would significantly interfere with their intended roles in the vaccine’s function.
Molecular dynamics simulations and analyses
The dynamic behavior of the monomeric spike SARS-CoV-2-ferritin complex was investigated through molecular dynamics simulations using GROMACS 2021.4 software (Hess et al., 2008). OPLS-AA/M force field (Jorgensen et al., 1996) was employed to generate parameters and topology files. The system was solvated in a cubic box of water molecules with a volume of 3,140.21 nm³, ensuring a 10 Å boundary around the proteins. TIP3P (Jorgensen et al., 1983) was used to simulate water molecules and counter ions were added to neutralize the system (13 Na + ions). A two-step approach involving 50,000 steps of steepest descent and 50,000 steps of conjugate gradient energy minimization was utilized to relieve the system from steric clashes. In the first step, the steepest descent algorithm was used with a maximum force tolerance of 500 kJ/mol/nm and an energy step size of 0.01 nm. Subsequently, the conjugate gradient algorithm was applied with a maximum force tolerance of 100 kJ/mol/nm and an energy step size of 0.01 nm, until convergence was achieved. The LINCS algorithm (Hess et al., 1997) was employed to constrain the bond lengths. Nonbonded interactions (electrostatic and VDW) were calculated with a cutoff value of 10 Å. The Particle Mesh Ewald method (Darden et al., 1993) with 0.16 nm Fourier grid spacing (Saini et al., 2019) was used for long-range electrostatic interactions. In the initial equilibration phase, the NVT ensemble was utilized to maintain a constant temperature of 310 K throughout a 100 ps simulation with a time step of 1 fs. The system was heated up from 0 K to 310 K over four steps using an annealing method. Following this, the temperature was maintained at 310 K using the V-rescale temperature coupling algorithm, a modified version of the Berendsen thermostat. To ensure structural stability, position restraints were applied to the proteins to prevent large fluctuations. Following the NVT equilibration, the system underwent further equilibration through a 250 ps simulation under an NPT ensemble at 310 K and 1 bar pressure, with a time step of 1 fs. Temperature and pressure coupling were achieved employing V-rescale and Berendsen methods, respectively. Finally, under NPT, with the V-rescale thermostat maintaining 310 K and the Parrinello-Rahman pressure coupling algorithm maintaining 1 bar, three replicates of 300 ns with a time step of 2 fs were conducted.
Bio3D package (Grant et al., 2006) in R was used for calculating the Root Mean Squared Deviation (RMSDs), Root Mean Squared Fluctuations (RMSFs), Principal Component Analysis (PCA), and Cross-Correlation of Cα atoms (DCCMs) for each replicate simulation during 300 ns of simulation time. Distances between centers of masses of spike and ferritin were calculated using VMD (Visual Molecular Dynamics) software (Humphrey et al., 1996). Secondary structure timeline analyses were calculated over 300 ns of simulation time with VMD STRIDE plugin. VMD was also utilized for calculating Solvent Accessible Surface Area (SASA). The visualization of structures was performed using PyMOL (LLC, 2015) and Figure 2 was created using BioRender.com.
Indicator for immune response simulation
The immunological simulation of the monomeric spike-ferritin vaccine was performed using the C-IMMSIM server (Rapin et al., 2010; Castiglione et al., 2021) to describe the immune response profile and immunogenicity of the antigenic peptide. The experiment was carried out by adding FASTA sequences of vaccine constructions into the simulation, which lasted approximately 12 months (a time step is about 8 h). The random seed and simulation volume was set to 12,345 and 50. Three in silico injections with no LPS were administered at time steps 1, 84, and 168, with a minimum time interval of 30 days between two injections. Additionally, to ensure the reliability of the simulation results, two additional simulations were carried out, with the random seed set to different numbers: one with the seed 3,452 and the other with 15,342 (Supplementary Figures S5, 6).
Predicting T cell epitopes on the monomeric spike protein of the designed vaccine
The peptides from monomeric spike proteins of the designed vaccine that are naturally processed by Major Histocompatibility complex (MHC) class I molecules were predicted using the Immune Epitope Database (IEDB) (Vita et al., 2019). This analysis employed the IEDB MHC prediction tool, MHC-NP (Giguère et al., 2013), which utilizes machine learning algorithms to predict peptides with potential binding affinity to MHC class I molecules. The ectodomain sequence of the monomeric spike protein of SARS-CoV-2 was utilized as the input sequence for the analysis. To ensure relevance to human immune responses, the chosen MHC source species was human. Additionally, to ensure comprehensive coverage of potential T-cell epitopes, five common Human Leukocyte Antigen (HLA)-B variants were selected as the MHC alleles. Peptide lengths ranging from 8 to 11 amino acids (all lengths) were assigned to capture a diverse range of potential epitopes. Finally, 3,620 peptides were identified for each of the examined variants, each demonstrating different affinities for MHC class I binding. Furthermore, MHC-II binding predictions were conducted to identify potential T-cell epitopes, utilizing the NetMHCIIpan EL prediction method (Reynisson et al., 2020). This analysis spanned across the 19 most common HLA-DR alleles, ensuring a broad representation of human genetic diversity. The default peptide length of 15 amino acids was employed for these predictions. This approach aimed to identify epitopes with high binding affinity to MHC-II molecules, thereby enabling investigation into vaccine efficacy across genetic variations in the host immune system.
Results
3D constructs of the monomeric spike within the ferritin nanocage
The ferritin nanoparticle vaccine was modeled with the monomeric spike antigen by molecular docking. The docking score obtained for the selected binding pose and the RMSD value were −88 and 1.1 (Å), respectively (Figure 3). The ferritin full-cage structure was assembled with 24 ferritin single units and has shown twenty monomeric spikes on ferritin structure, with four ferritin units remaining without adjuvants due to steric hindrance (Figure 4).
Figure 3. (A) Molecular docking result of the monomeric spike with ferritin unit from the HADDOCK server. (B) The best binding pose is presented in the right panel in purple surface while the ferritin unit is shown in red surface.
Figure 4. A Full nanocage model of the monomeric spike with ferritin is shown in the picture. Spikes are presented in the purple cartoon, while ferritin units are shown in the yellow cartoon.
The efficiency of the monomeric spike for inhibition of different SARS-CoV-2 variants
By examining the mutations that occurred in the spike protein of variants of concern, including Alpha, Beta, Gamma, Delta, and Omicron (Supplementary Figure S7), it was noticed that not only do none of the mutations exist in the ferritin binding interface but the center of occurrence of mutations is also far from the binding interface of ferritin in the monomeric spike vaccine model. While the mutations D1118H and H1101D, associated with the Alpha and Beta variants respectively, are in close proximity to the binding interface, it is crucial to note that they are not directly situated within it. The binding interface is specifically defined by residues within the range of 1,136–1,147, where neither D1118H nor H1101D fall within this range.
Thus, the monomeric spike-ferritin vaccine model can still be effective despite the known mutations.
Structural stability and flexibility of the monomeric spike-ferritin complex
Molecular dynamics simulations were utilized to thoroughly examine the structural stability and flexibility of the monomeric spike-ferritin complex and its constituent parts. The RMSD analysis was employed to assess the global dynamics of the monomeric spike-ferritin structure. The general trend of time-dependent RMSD plots indicated an initial increase, followed by stabilization after 200 ns (Figure 5A). The average calculated RMSD values across three replicate simulations were 16.5 Å, 27.5 Å, and 17.4 Å, respectively. Notably, the observed high RMSD values are primarily attributed to the substantial size of the monomeric spike-ferritin nanoparticle structure (1,200 a. a) (Figure 2). These observations aligned with the expectation of reduced system stability due to the size of the spike monomer, which naturally results in a higher degree of conformational changes. The distribution patterns of calculated RMSD values derived from the RMSD histogram aligned with the findings of the time-dependent RMSD plots (Figure 5B). All replicates of the overall system revealed stable RMSD values with peaks at 16 Å. That is expected due to multiple domains within its structure (Figures 1, 2).
Figure 5. (A) Time-dependent RMSD plots of the monomeric spike-ferritin complex during three replicates of 300ns MD simulations. RMSD values were calculated for Cα atoms, with the X-axis indicating the Time (ns) and the Y-axis showing RMSD values in Å. (B) The histogram depicts the distribution of RMSD values for Cα atoms during a 300ns MD simulation. Values are based on calculations from all simulation replicates.
The flexibility of the complex was further investigated through RMSF analysis, focusing on the Cα atoms per residue (Figure 6A). The average calculated RMSF values were 6.3 Å, 8.9 Å, and 8.4 Å for the three replicate simulations. These findings corroborated the insights gained from RMSD analysis and indicated that the monomeric spike-ferritin complex is highly stable. The RMSF plots revealed multiple peaks, which indicated the presence of distinct flexible loops and their dynamic motions within the structure. Specifically, the RBD region exhibited prominent peaks at residues (476–490) and (498–507) (Figure 6B), highlighting the dynamic motions within these specific loops. Similarly, fluctuations were observed in the ferritin structure, particularly at residues (155–158), suggesting the dynamic nature of the loops within its structure (Figure 6B). However, it is important to note that while these loops were flexible, they did not directly interfere in the spike-ferritin interface and thus the stability of the designed vaccine may not be affected. Loops one and 2, which were located in the interface of RBD with ACE2, were far from the monomeric spike interface with ferritin, and, therefore, would not affect the stability of the modeled vaccine. Loop three in ferritin was also located far from the interface. Thus, despite the flexibility indicated by the peaks in the RMSF plots, these loops did not influence the stability of the designed vaccine.
Figure 6. (A) The RMSF of the Cα atoms within the complex during three replicates of 300ns MD simulations. The X-axis corresponds to the residue number of the monomeric spike-ferritin, while the Y-axis shows RMSF values in Å. The dashed line separates the monomeric spike and ferritin amino acids. The first 1,058 amino acids correspond to the monomeric spike, while the subsequent 174 amino acids represent the ferritin. Dashed boxes highlight significant peaks in the RMSF plots. (B) The position of each loop associated with the peaks in the RMSF plots is presented.
The Principal Component Analysis (PCA) and its results were useful in revealing the structural changes in the monomeric spike-ferritin complex (Supplementary Figure S8). PCA, which effectively reduces the dimensionality, was performed to identify the most significant fluctuation modes of the complex. Notably, the percentage of variance captured by PC1 and PC2 of the system amounted to 73.23% and 13.88%, respectively. These findings directly demonstrated that structural fluctuations are noticeable in the systems.
To gain deeper insights into the interactions within the RBD of the monomeric spike and ferritin, separate RMSD and RMSF analyses were conducted for these components (Figure 7). Notably, three replicate simulations displayed rather consistent structural stability for RBD, with average RMSD values of 4.2 Å, 3.0 Å, and 3.6 Å (Figure 7A). Suggesting that when RBD is present within the full ectodomain spike monomer, it remains structurally stable, as evidenced by the RMSF analysis as well. The average RMSF values obtained for RBD were 1.3 Å, 1.3 Å, and 1.5 Å, respectively (Figure 7B).
Figure 7. RMSD (A) and RMSF (B) for the Cα atoms of the Receptor Binding Domain (RBD) during the simulation time. RMSD (C) and RMSF (D) for the Cα atoms of the ferritin during three replicated of 300ns MD simulations. Values are based on calculations from all simulation replicates.
Similarly, the examination of ferritin across three replicates yielded RMSD values averaging at 3.0 Å, 2.4 Å, and 2.5 Å (Figure 7C). These RMSD plots indicated the significant stability of ferritin when it creates a complex with the full ectodomain spike monomer.
Analysis of the secondary structure of the monomeric spike
To gain insights into crucial changes in protein structure during the molecular dynamics simulations, the secondary structures were calculated for the trajectories (Figure 8). The secondary structure analysis of the monomeric spike-ferritin complex did not show any significant changes in its dynamics, indicating overall stability of the complex.
Figure 8. The secondary structures of all residues in the monomeric spike-ferritin complex were analyzed during a 300 ns duration of MD simulation for all replicates (A–C).
In summary, it was demonstrated that the RBD and ferritin components within the monomeric spike-ferritin complex maintain strong stability. These findings have significant implications, suggesting that the full ectodomain spike monomer, despite its substantial size, holds promise for vaccine development against SARS-CoV-2.
In the analysis of spike-ferritin interactions, the cross-correlation matrices of the Cα atoms were calculated in all simulations (Figure 9). Positive and negative values were used to denote correlated and anti-correlated motions, respectively. Anti-correlated motions occurring between spike and ferritin supported the stable interaction between the two proteins. It was also observed that when the spike becomes compact, the anti-correlation between RBD and ferritin grows as they come closer to each other (Figure 9B, C; Figure 10E, F).
Figure 9. Residue cross-correlation matrix of monomeric spike-ferritin was calculated over three replicates (A–C) of 300 ns MD simulations for Cα atoms. The color bar on the right side displays correlation values, ranging from −1 (anti-correlated) to 1 (correlated). Pink represents anti-correlated motions and cyan indicates correlated motions. Residues of monomeric spike, RBD, and ferritin are shown with blue, yellow, and red lines respectively on the x and y-axes. Motions between RBD-ferritin (box 1) and S2-ferritin (box 2) are shown with dashed boxes.
Figure 10. Structural presentations of monomeric spike attached to one unit of ferritin. Panels (A–C) represent the initial conformations of three replicate simulations. Panels (D–F) are final conformations. spike monomer, RBD, and ferritin are shown in violet, cyan, and yellow cartoon plus transparent surfaces respectively.
Furthermore, visualizations of the dynamics indicated that the observed anti-correlated motions would not disturb the ‘‘up’’ active conformation of RBD (Figure 10). This is crucial because RBD needs to remain in an upward orientation to be accessible to ACE2.
The distance between the center of masses of ferritin and specific regions of the spike protein, including RBD, S1, and S2, was calculated over three replicates of 300ns MD simulations (Figure 11). In general, a descending trend was evident in the calculated distances during the simulation. Among calculated distances, S2-ferritin exhibited the least variation, while RBD-ferritin showed the greatest change in the overall trend. The observed decline in the distance was supported by the visualization of dynamics, indicating that the spike becomes more compact, with the S1 and S2 subunits of the spike approaching each other. However, the RBD maintains its upward conformation consistently throughout the entire simulation time (Figure 10).
Figure 11. Distance between center of masses of ferritin-spike, ferritin-RBD, ferritin-S1, and ferritin-S2 within three replicated of 300 ns MD simulations. Values are based on calculations from all simulation replicates.
Solvent Accessible Surface Area (SASA) calculations were performed on the RBD to track changes in its exposure and accessibility over time (Figure 12). The stable trend observed in the RBD SASA values suggested that despite structural changes in the spike protein, the RBD maintains its accessibility to the ACE2 receptor. This indicates that the RBD’s binding ability to the ACE2 receptor remains relatively constant, even during the compaction that happens to the spike protein.
Immune simulation of monomeric nanoparticle vaccine model
The immunogenic profile of the monomeric spike nanoparticle vaccine candidate was obtained from the C-IMMSIM server. The vaccine was injected in three doses over 2 months with 30-day intervals (Figure 13). It was shown that the modeled vaccine was able to elicit both humoral and cellular mediated immune responses. Following each vaccine dose, there was an increase in B cell population, indicating that the vaccine effectively stimulated the immune system and led to a more robust antibody response (Figure 13A). Additionally, the CD4+ T cell population peaked after the third injection (Figure 13B). Investigation of the CD8+ T cell population revealed a higher population of active cells compared to anergic cells (Figure 13C). Furthermore, the plasma cell population revealed a substantial increase (>80 cells per mm) following the third dose, suggesting robust antibody production (Figure 13D). The antibody response to the monomeric nanoparticle vaccine model (Figure 13E) exhibited a slight increase in the graphs after the first two doses. From the second dose onwards, there was a rise in antibody secretions attributed to memory T lymphocytes. An upward trend of IgM and IgG antibody titer was observed after the third injection, alongside a decline in antigen levels. Additionally, elevated levels of cytokines such as IFN-γ and IL-2 were observed (Figure 13F). Clearance of antigenic molecules was demonstrated after three doses of vaccination. There was also a significant increase in the populations of both B and T memory cells, which remained elevated for several months.
Figure 13. Immune response simulations of monomeric spike-ferritin nanoparticle vaccine with 1 month intervals. (A) B cell population, (B) CD4+ T cell population, (C) CD8+ T cell population, (D) Plasma B cell population, (E) Humoral immunity response, (F) Concentration of cytokines and interleukins.
Upon analysis of two additional simulations conducted with different random seeds (Supplementary Figures S5, S6), it became evident that the immune responses observed across all three simulations were comparable. Minor variations were noted, with indications of improved immune responses observed in some cases. Importantly, the pattern of immune system stimulation remained consistent across all three simulations.
HLA variants and immune response to monomeric nanoparticle vaccine
The adaptive immune response to SARS-CoV-2 relies on the activity of T cells, wherein MHC molecules play a crucial role in processing and presenting antigenic peptides to T cells, initiating immune responses (Lin et al., 2023). Human MHC molecules, also referred to as Human Leukocyte Antigen (HLA), are divided into two classes. MHC class I (HLA-A, HLA-B, and HLA-C) present peptides from endogenous sources to CD8+ T cells, while MHC class II molecules (HLA-DR, HLA-DP, and HLA-DQ) predominantly process exogenous peptides, presenting them to CD4+ T cells (Traherne, 2008). The diversity in HLA alleles contributes to variations in peptide processing and presentation, thereby influencing immune responses (Wolday et al., 2023).
In order to ensure the effectiveness of the designed vaccine against different variants of HLA, the prediction of T cell epitopes on the monomeric spike protein of the vaccine was performed. However, due to the high diversity in HLA alleles, assessing all variants was not feasible. Instead, the study focused on investigating the two most polymorphic MHC genes. The HLA-B gene, with 1,077 reported alleles, is recognized as the most polymorphic among MHC class I genes, while the HLA-DRB1 gene, with 669 alleles, stands out as the most polymorphic among MHC class II genes (Shiina et al., 2009; Aliseychik et al., 2018). The analysis of monomeric spike peptides processed by HLA-B variants revealed notable findings. For each of the examined variants, 3,620 peptides were predicted. These peptides were then assessed based on their affinity for binding to HLA-B molecules. The prob score distribution of monomeric spike peptides processed by five common HLA-B variants revealed distinct patterns (Supplementary Figure S9). Overall, variant HLA-B*57:01 exhibited the highest prob score distribution, with three peptides scoring as high as 1, indicating a strong probability of binding to HLA-B molecules. Following closely were variants HLA-B*53:01 and HLA-B*44:03, which also showed relatively high prob score distributions. In contrast, the low prob score distribution and median for variants HLA-B*35:01 and HLA-B*07:02 suggested a lower probability of binding processed peptides to HLA-B molecules. Remarkably, despite the diversity in HLA-B alleles, peptides with a high probability of binding to HLA-B molecules (>70%) were identified for all of five most common HLA-B variants. Further exploration involved predicting peptides that bind to MHC class II variants. A thorough examination assessed the binding affinity of predicted peptides across 13 alleles of HLA-DRB1 and its paralogues: HLA-DRB3, HLA-DRB4, and HLA-DRB5 (Supplementary Figure S10). Remarkably, a consistent pattern was observed across all alleles. While the majority of predicted peptides exhibited a binding affinity distribution ranging from 0 to 0.1, it is noteworthy that all alleles (except for HLA-DRB4*01:01) presented peptides with a high probability (>60%) of binding to MHC class II molecules. These observations underscore the robustness of the immune system’s response to the monomeric nanoparticle vaccine across both MHC class I and class II pathways, despite variations in MHC alleles.
Discussion
This study highlighted the role of in silico structural modeling in rapidly generating information on novel viruses to help predict their behavior and aid in countermeasure development. The main focus of this work was on the SARS-CoV-2 which caused the extensive pandemic of COVID-19. The structural spike protein of the virus has been revealed as the major component responsible for its pathogenicity (Du et al., 2009; Yang et al., 2020). This protein plays a key role in facilitating the virus’s entry into human cells by binding to a receptor called ACE2 (Hoffmann et al., 2020; Walls et al., 2020). Upon attachment to ACE2, the structural spike protein enables the virus to infect the cell, leading to the development of the disease. The RBD is the most immunogenic domain interacting directly with ACE2 (Fehr and Perlman, 2015; Shereen et al., 2020); thus, the focus of vaccine development in this study was on presenting this particular domain within a monomeric spike protein structure. Ferritin, primarily a cytosolic protein serving as a nanocage for iron storage, is also secreted in small amounts in the serum as an iron carrier (Andrews et al., 1992). In this study, nanoparticle ferritin was employed to display monomeric spikes as an immune system stimulant. The ferritin assembly used in this study lacked iron ions, addressing concerns about ferritin releasing iron into the bloodstream or tissues and causing potential interference. Ferritin nanoparticle vaccines, previously created to combat other infections, are now designed to fight SARS-CoV-2 and produce significant immunogenicity (Sliepen et al., 2015; Wang et al., 2017; Kalathiya et al., 2021; Powell et al., 2021; Wang et al., 2021; Wuertz et al., 2021). In this study, structural molecular modeling was employed to construct and explore the dynamics of the ferritin nanocage monomeric spike SARS-CoV-2 vaccine at the atomistic level for the first time. The vaccine was developed using the monomeric spike as a crucial viral component in its construction. Through molecular docking and consequently placing the ferritin units in an icosahedral full-cage structure, a monomeric spike model of the ferritin nanoparticle vaccine was obtained. Additionally, molecular dynamics simulations were conducted, providing a comprehensive analysis of the dynamic behavior and interactions within the monomeric spike-ferritin complex.
According to the analysis of trajectories, the monomeric spike-ferritin complex revealed overall stability and structural integrity, which are crucial for its functionality. Additionally, the native conformation of the antigen, necessary for immune recognition, remained preserved. Furthermore, the proper orientation of spike protein monomers within the ferritin unit allowed critical epitopes to be exposed to the immune system. The components of the vaccine demonstrated proper folding, maintaining their functional structures. Simulation of the immunity indicated that the monomeric spike model would stimulate cellular and humoral immunity in the human body and reduce the antigen amount. In addition, 2,458 peptides were identified with a high probability (>70%) of binding to MHC molecules of the five most common HLA-B variants, indicating their potential as epitopes for immune response stimulation. Notably, the HLA-B*57:01 variant demonstrated significantly higher peptide binding affinity compared to other HLA-B variants, presenting an opportunity for personalized vaccine design. Specifically, populations with a high frequency of the HLA-B*57:01 allele could benefit from a customized vaccine formulation designed to maximize immune response efficacy. Also, all examined HLA-DRB1 alleles exhibited peptides with a high probability (>60%) of binding to MHC class II molecules, demonstrating the potential for robust immune activation. Finally, the reliability of the modeled interfaces was also investigated by examining the mutation sites related to variants of concern, including Alpha, Beta, Gamma, and Omicron, on the spike protein. The results revealed that none of these mutations occurred at the interfaces of the model. Also, in a study that predicted possible mutation locations, three out of five mutations occurred in the previous variants, and two more mutations, including Y489 and T500, are expected to occur in the following variants. However, these threatening mutations are not only not located at the binding surface of ferritin and spike but also far from it (Mohammad et al., 2021; Alkhatib et al., 2022).
In previous studies, the efficacy of the RBD-ferritin nanoparticle vaccine’s immunogenicity has been established (Kalathiya et al., 2021; Masoomi Nomandan et al., 2022). The studies successfully developed RBD-based subunit vaccines employing ferritin protein nanocages. Their findings highlight the vaccine’s efficacy in presenting RBD, ensuring stability, and triggering a robust immune response. These efforts to develop a SARS-CoV-2 nanoparticle vaccine have predominantly focused on utilizing the RBD. However, experimental evidence has indicated that the full-length spike may be more effective (Powell et al., 2021). This current study aimed to investigate the suitability of using spike monomer rather than just RBD in the design of subunit vaccines. Utilizing the full-length ectodomain monomeric spike protein ensures that various epitopes outside the RBD region are also presented to the immune system, potentially offering protection against different variants of the virus. In this study, it was observed that the monomeric spike–ferritin complex, which is a large system with about 1,200 a. a, is highly stable compared to the RBD-ferritin vaccine construct. Immunogenicity results were comparable to those of the RBD-ferritin vaccine (Masoomi Nomandan et al., 2022), indicating that the designed vaccine is capable of producing an immune response. In general, the results obtained from this vaccine supported the potential effectiveness of utilizing the monomeric spike in nanoparticle vaccine designs.
As mentioned in several previous studies (Kalathiya et al., 2021; Powell et al., 2021), multivalent presentation of SARS-CoV-2 spike on ferritin can notably enhance the elicitation of neutralizing antibodies, constituting a viable strategy for single-dose vaccination against COVID-19. Although ferritin nanoparticle vaccines have already demonstrated promising clinical benefits against several microbial infections, including SARS-CoV-2, the molecular mechanism of these nanoparticle assemblies and the underlying atomistic interactions were poorly understood.
The SARS-CoV-2 spike protein is heavily glycosylated (Shajahan et al., 2020; Yasunori et al., 2020; Gong et al., 2021). Thus, the glycosylation of the monomeric spike in the SARS-CoV-2 nanoparticle vaccine holds significant potential to influence the immune response elicited by the designed vaccine in various ways. Glycosylation has been identified as an important factor impacting protein folding, stability, and ligand binding (Lee et al., 2015; Watanabe et al., 2019; Gong et al., 2021; Rahnama et al., 2021). Proper glycosylation of the vaccine spike protein is important for initiating a robust adaptive immune response, as it affects the folding of vaccine epitopes and antigen processing by immune cells (Reis et al., 2021). On the other hand, glycans have the potential to act as camouflage, concealing immunogenic epitopes from detection by the innate immune system (Casalino et al., 2020). This mechanism could potentially reduce the vaccine’s ability to induce effective immune responses against the virus. The simulations conducted in this research focused on structures without attached glycans. However, investigation of the N- and O-glycosylation profile of the spike protein ensured that there are no potential glycosylation sites within the ferritin-spike interface (residues: 1,147–1,136). Nevertheless, future studies could examine the presence and effects of glycans on the designed vaccine for a comprehensive understanding of vaccine immunogenicity. While in silico vaccine design offers invaluable insights, its efficacy and reliability depend on subsequent experimental confirmation through in vitro and in vivo studies. The study acknowledges the importance of integrating computational models with experimental data to advance the development of reliable and effective vaccines.
Conclusion
Ferritin nanoparticle vaccines are one of the most effective vaccines that provide optimal immunity against the virus without threatening the patient’s health. In this work, molecular docking, molecular dynamics simulation, and immune response simulations were used to construct and evaluate the self-assembled ferritin nanoparticle and monomeric spike vaccine model for the first time. The monomeric spike-ferritin complex displayed significant stability, particularly in the RBD and ferritin. This stability was further supported with minimal changes observed in the secondary structure of spike and ferritin. Despite the observed compaction of the spike protein throughout the simulations, the accessibility of the RBD remained unchanged. Strong interactions between the monomeric spike and ferritin were indicated by their observed motions toward each other. Furthermore, the sustained upward conformation of the RBD potentially facilitated interaction with ACE2. The study on mutation sites related to variants of concern found no mutations at the modeled interfaces. Additionally, while potential future mutations (Y489 and T500) are predicted in subsequent variants, they are not positioned near crucial binding sites. These findings could be crucial in improving current ferritin nanoparticle vaccines as well as future nanoparticle vaccine formulations in both in silico and in vitro investigations.
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Author contributions
FG: Data curation, Formal Analysis, Visualization, Writing–original draft. MaA: Conceptualization, Investigation, Project administration, Supervision, Writing–review and editing. ZM: Data curation, Formal Analysis, Visualization, Writing–original draft. MoA: Data curation, Formal Analysis, Writing–original draft.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2024.1403635/full#supplementary-material
References
Acuti Martellucci, C., Flacco, M. E., Cappadona, R., Bravi, F., Mantovani, L., and Manzoli, L. (2020). SARS-CoV-2 pandemic: an overview. Adv. Biol. Regul., 77, 100736. doi:10.1016/j.jbior.2020.100736
Alexandersen, S., Chamings, A., and Bhatta, T. R. (2020). SARS-CoV-2 genomic and subgenomic RNAs in diagnostic samples are not an indicator of active replication. Nat. Commun. 11 (1), 6059. doi:10.1038/s41467-020-19883-7
Aliseychik, M. P., Andreeva, T. V., and Rogaev, E. I. (2018). Immunogenetic factors of neurodegenerative diseases: the role of HLA class II. Biochem. Biokhimiia 83 (9), 1104–1116. doi:10.1134/S0006297918090122
Alkhatib, M., Salpini, R., Carioti, L., Ambrosio, F. A., D’Anna, S., Duca, L., et al. (2022). Update on SARS-CoV-2 Omicron variant of concern and its peculiar mutational profile. Microbiol. Spectr. 10 (2), e0273221. doi:10.1128/spectrum.02732-21
Amanat, F., and Krammer, F. (2020). SARS-CoV-2 vaccines: status report. Immunity, 52(4), 583–589. doi:10.1016/j.immuni.2020.03.007
Amanat, F., Stadlbauer, D., Strohmeier, S., Nguyen, T. H. O., Chromikova, V., McMahon, M., et al. (2020). A serological assay to detect SARS-CoV-2 seroconversion in humans. Nat. Med. 26 (7), 1033–1036. doi:10.1038/s41591-020-0913-5
Andersen, K. G., Rambaut, A., Lipkin, W. I., Holmes, E. C., and Garry, R. F. (2020). The proximal origin of SARS-CoV-2. Nat. Med. 26 (4), 450–452. doi:10.1038/s41591-020-0820-9
Andrews, S. C., Harrison, P. M., Yewdall, S. J., Arosio, P., Levi, S., Bottke, W., et al. (1992). Structure, function, and evolution of ferritins. J. Inorg. Biochem., 47(1), 161–174. doi:10.1016/0162-0134(92)84062-R
Bachmayer, H., Liehl, E., and Schmidt, G. (1976). Preparation and properties of a novel influenza subunit vaccine. Postgrad. Med. J. 52 (608), 360–367. doi:10.1136/pgmj.52.608.360
Bhushan, B., Kumar, S. U., Matai, I., Sachdev, A., Dubey, P., and Gopinath, P. (2014). Ferritin nanocages: a novel platform for biomedical applications. J. Biomed. Nanotechnol. 10 (10), 2950–2976. doi:10.1166/jbn.2014.1980
Borges, O., Cordeiro-da-Silva, A., Tavares, J., Santarém, N., de Sousa, A., Borchard, G., et al. (2008). Immune response by nasal delivery of hepatitis B surface antigen and codelivery of a CpG ODN in alginate coated chitosan nanoparticles. Eur. J. Pharm. Biopharm. Official J. Arbeitsgemeinschaft Fur Pharmazeutische Verfahrenstechnik e.V 69 (2), 405–416. doi:10.1016/j.ejpb.2008.01.019
Carmen, J. M., Shrivastava, S., Lu, Z., Anderson, A., Morrison, E. B., Sankhala, R. S., et al. (2021). SARS-CoV-2 ferritin nanoparticle vaccine induces robust innate immune activity driving polyfunctional spike-specific T cell responses. Npj Vaccines 6 (1), 151. doi:10.1038/s41541-021-00414-4
Casalino, L., Gaieb, Z., Goldsmith, J. A., Hjorth, C. K., Dommer, A. C., Harbison, A. M., et al. (2020). Beyond shielding: the roles of glycans in the SARS-CoV-2 spike protein. ACS central Sci. 6 (10), 1722–1734. doi:10.1021/acscentsci.0c01056
Castiglione, F., Deb, D., Srivastava, A. P., Liò, P., and Liso, A. (2021). From infection to immunity: understanding the response to SARS-CoV2 through in-silico modeling. Front. Immunol. 12, 646972. doi:10.3389/fimmu.2021.646972
Chakraborti, S., and Chakrabarti, P. (2019). Self-assembly of ferritin: structure, biological function and potential applications in nanotechnology. Adv. Exp. Med. Biol. 1174, 313–329. doi:10.1007/978-981-13-9791-2_10
Champion, C. I., Kickhoefer, V. A., Liu, G., Moniz, R. J., Freed, A. S., Bergmann, L. L., et al. (2009). A vault nanoparticle vaccine induces protective mucosal immunity. PloS One 4 (4), e5409. doi:10.1371/journal.pone.0005409
Chattopadhyay, S., Chen, J.-Y., Chen, H.-W., and Hu, C.-M. J. (2017). Nanoparticle vaccines adopting virus-like features for enhanced immune potentiation. Nanotheranostics 1, 244–260. doi:10.7150/ntno.19796
Chen, W. H., Hotez, P. J., and Bottazzi, M. E. (2020b). Potential for developing a SARS-CoV receptor-binding domain (RBD) recombinant protein as a heterologous human vaccine against coronavirus infectious disease (COVID)-19. Hum. Vaccines Immunother. 16 (6), 1239–1242. doi:10.1080/21645515.2020.1740560
Chen, W.-H., Strych, U., Hotez, P. J., and Bottazzi, M. E. (2020a). The SARS-CoV-2 vaccine pipeline: an overview. Curr. Trop. Med. Rep. 7, 61–64. doi:10.1007/s40475-020-00201-6
Choi, Y. K., Cao, Y., Frank, M., Woo, H., Park, S.-J., Yeom, M. S., et al. (2021). Structure, dynamics, receptor binding, and antibody binding of the fully glycosylated full-length SARS-CoV-2 spike protein in a viral membrane. J. Chem. Theory Comput. 17 (4), 2479–2487. doi:10.1021/acs.jctc.0c01144
Darden, T., York, D., and Pedersen, L. (1993). Particle mesh Ewald: an N log (N) method for Ewald sums in large systems. J. Chem. Phys. 98, 10089–10092. doi:10.1063/1.464397
Du, L., He, Y., Zhou, Y., Liu, S., Zheng, B.-J., and Jiang, S. (2009). The spike protein of SARS-CoV — a target for vaccine and therapeutic development. Nat. Rev. Microbiol. 7 (3), 226–236. doi:10.1038/nrmicro2090
Edwardson, T. G. W., and Hilvert, D. (2019). Virus-inspired function in engineered protein cages. J. Am. Chem. Soc. 141 (24), 9432–9443. doi:10.1021/jacs.9b03705
Fehr, A. R., and Perlman, S. (2015). Coronaviruses: an overview of their replication and pathogenesis. Methods Mol. Biol. Clift. N.J. 1282, 1–23. doi:10.1007/978-1-4939-2438-7_1
Ganji, A., Mosayebi, G., Khaki, M., and Ghazavi, A. (2020). A review of the 2019 novel coronavirus (Covid-19): immunopathogenesis, molecular biology and clinical aspects. J. Arak Univ. Med. Sci. 23 (1), 8–21. doi:10.32598/JAMS.23.1.51.5
Giguère, S., Drouin, A., Lacoste, A., Marchand, M., Corbeil, J., and Laviolette, F. (2013). MHC-NP: predicting peptides naturally processed by the MHC. J. Immunol. methods 400-401, 30–36. doi:10.1016/j.jim.2013.10.003
Glück, R., Moser, C., and Metcalfe, I. C. (2004). Influenza virosomes as an efficient system for adjuvanted vaccine delivery. Expert Opin. Biol. Ther. 4 (7), 1139–1145. doi:10.1517/14712598.4.7.1139
Gong, Y., Qin, S., Dai, L., and Tian, Z. (2021). The glycosylation in SARS-CoV-2 and its receptor ACE2. Signal Transduct. Target. Ther. 6, 396. doi:10.1038/s41392-021-00809-8
Gordon, J. M., D, K. H. A., Ines, E.-N., Aslaa, A. K. P. K., Camila, M. C., and Kayvon, M. (2022). A SARS-CoV-2 ferritin nanoparticle vaccine elicits protective immune responses in nonhuman primates. Sci. Transl. Med. 0 (0), eabi5735. doi:10.1126/scitranslmed.abi5735
Grant, B. J., Rodrigues, A. P. C., ElSawy, K. M., McCammon, J. A., and Caves, L. S. D. (2006). Bio3D: an R package for the comparative analysis of protein structures. Bioinformatics 22, 2695–2696. doi:10.1093/bioinformatics/btl461
Gui, M., Song, W., Zhou, H., Xu, J., Chen, S., Xiang, Y., et al. (2017). Cryo-electron microscopy structures of the SARS-CoV spike glycoprotein reveal a prerequisite conformational state for receptor binding. Cell Res. 27 (1), 119–129. doi:10.1038/cr.2016.152
Haiou, L., Yunjiao, Z., Meng, Z., Haizhou, W., Qiu, Z., and Jing, L. (2022). Updated approaches against SARS-CoV-2. Antimicrob. Agents Chemother. 64 (6), 004833–e520. doi:10.1128/AAC.00483-20
He, L., de Val, N., Morris, C. D., Vora, N., Thinnes, T. C., Kong, L., et al. (2016). Presenting native-like trimeric HIV-1 antigens with self-assembling nanoparticles. Nat. Commun. 7 (1), 12041. doi:10.1038/ncomms12041
Hess, B., Bekker, H., Berendsen, H. J. C., and Fraaije, J. G. E. M. (1997). LINCS: a linear constraint solver for molecular simulations. J. Comput. Chem. 18, 1463–1472. doi:10.1002/(sici)1096-987x(199709)18:12<1463::aid-jcc4>3.0.co;2-h
Hess, B., Kutzner, C., van der Spoel, D., and Lindahl, E. (2008). GROMACS 4: algorithms for highly efficient, load-balanced, and scalable molecular simulation. J. Chem. Theory Comput. 4 (3), 435–447. doi:10.1021/ct700301q
Hoffmann, M., Kleine-Weber, H., and Pöhlmann, S. (2020). A multibasic cleavage site in the spike protein of SARS-CoV-2 is essential for infection of human lung cells. Mol. Cell 78 (4), 779–784. doi:10.1016/j.molcel.2020.04.022
Hu, B., Guo, H., Zhou, P., and Shi, Z.-L. (2021). Characteristics of SARS-CoV-2 and COVID-19. Nat. Rev. Microbiol. 19 (3), 141–154. doi:10.1038/s41579-020-00459-7
Humphrey, W., Dalke, A., and Schulten, K. (1996). VMD: visual molecular dynamics. J. Mol. Graph. 14, 33–38. doi:10.1016/0263-7855(96)00018-5
Jackson, C. B., Farzan, M., Chen, B., and Choe, H. (2022). Mechanisms of SARS-CoV-2 entry into cells. Nat. Rev. Mol. Cell Biol. 23 (1), 3–20. doi:10.1038/s41580-021-00418-x
Jin, Y., He, J., Fan, K., and Yan, X. (2019). Ferritin variants: inspirations for rationally designing protein nanocarriers. Nanoscale 11 (26), 12449–12459. doi:10.1039/C9NR03823J
Jorgensen, W. L., Chandrasekhar, J., Madura, J. D., Impey, R. W., and Klein, M. L. (1983). Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 79, 926–935. doi:10.1063/1.445869
Jorgensen, W. L., Maxwell, D. S., and Tirado-Rives, J. (1996). Development and testing of the OPLS all-atom force field on conformational energetics and properties of organic liquids. J. Am. Chem. Soc. 118, 11225–11236. doi:10.1021/ja9621760
Jutz, G., van Rijn, P., Santos Miranda, B., and Böker, A. (2015). Ferritin: a versatile building block for bionanotechnology. Chem. Rev. 115 (4), 1653–1701. doi:10.1021/cr400011b
Kalathiya, U., Padariya, M., Fahraeus, R., Chakraborti, S., and Hupp, T. R. (2021). Multivalent display of SARS-CoV-2 spike (RBD domain) of COVID-19 to nanomaterial, protein ferritin nanocages. Biomolecules 11 (2), 297. doi:10.3390/biom11020297
Kamp, H. D., Swanson, K. A., Wei, R. R., Dhal, P. K., Dharanipragada, R., Kern, A., et al. (2020). Design of a broadly reactive Lyme disease vaccine. Npj Vaccines 5 (1), 33. doi:10.1038/s41541-020-0183-8
Kanekiyo, M., Bu, W., Joyce, M. G., Meng, G., Whittle, J. R. R., Baxa, U., et al. (2015). Rational design of an epstein-barr virus vaccine targeting the receptor-binding site. Cell 162 (5), 1090–1100. doi:10.1016/j.cell.2015.07.043
Kanekiyo, M., Wei, C.-J., Yassine, H. M., McTamney, P. M., Boyington, J. C., Whittle, J. R. R., et al. (2013). Self-assembling influenza nanoparticle vaccines elicit broadly neutralizing H1N1 antibodies. Nature 499 (7456), 102–106. doi:10.1038/nature12202
Kozakov, D., Hall, D. R., Xia, B., Porter, K. A., Padhorny, D., Yueh, C., et al. (2017). The ClusPro web server for protein-protein docking. Nat. Protoc. 12 (2), 255–278. doi:10.1038/nprot.2016.169
Lawson, D. M., Artymiuk, P. J., Yewdall, S. J., Smith, J. M. A., Livingstone, J. C., Treffry, A., et al. (1991). Solving the structure of human H ferritin by genetically engineering intermolecular crystal contacts. Nature 349 (6309), 541–544. doi:10.1038/349541a0
Lee, H., Qi, Y., and Im, W. (2015). Effects of N-glycosylation on protein conformation and dynamics: protein Data Bank analysis and molecular dynamics simulation study. Sci. Rep. 5, 8926. doi:10.1038/srep08926
Li, F., Li, W., Farzan, M., and Harrison, S. C. (2005). Structure of SARS coronavirus spike receptor-binding domain complexed with receptor. Sci. (New York, N.Y.) 309 (5742), 1864–1868. doi:10.1126/science.1116480
Lin, F., Lin, X., Fu, B., Xiong, Y., Zaky, M. Y., and Wu, H. (2023). Functional studies of HLA and its role in SARS-CoV-2: stimulating T cell response and vaccine development. Life Sci. 315, 121374. doi:10.1016/j.lfs.2023.121374
Liu, Z., Zhou, H., Wang, W., Tan, W., Fu, Y.-X., and Zhu, M. (2014). A novel method for synthetic vaccine construction based on protein assembly. Sci. Rep. 4 (1), 7266. doi:10.1038/srep07266
Masoomi Nomandan, S. Z., Azimzadeh Irani, M., and Hosseini, S. M. (2022). In silico design of refined ferritin-SARS-CoV-2 glyco-RBD nanoparticle vaccine. Front. Mol. Biosci. 9, 976490. doi:10.3389/fmolb.2022.976490
Mohammad, A., Valentina, S., Romina, S., Alessandra, A. F., Concetta, B. M., Luca, C., et al. (2021). SARS-CoV-2 variants and their relevant mutational profiles: update summer 2021. Microbiol. Spectr. 9 (3), 010966–e1121. doi:10.1128/Spectrum.01096-21
Mulligan, M. J., Lyke, K. E., Kitchin, N., Absalon, J., Gurtman, A., Lockhart, S., et al. (2020). Phase I/II study of COVID-19 RNA vaccine BNT162b1 in adults. Nature 586 (7830), 589–593. doi:10.1038/s41586-020-2639-4
Nagae, M., Ikeda, A., Hanashima, S., Kojima, T., Matsumoto, N., Yamamoto, K., et al. (2016). Crystal structure of human dendritic cell inhibitory receptor C-type lectin domain reveals the binding mode with N-glycan. FEBS Lett. 590 (8), 1280–1288. doi:10.1002/1873-3468.12162
Nochi, T., Yuki, Y., Takahashi, H., Sawada, S., Mejima, M., Kohda, T., et al. (2010). Nanogel antigenic protein-delivery system for adjuvant-free intranasal vaccines. Nat. Mater. 9 (7), 572–578. doi:10.1038/nmat2784
Oyewumi, M. O., Kumar, A., and Cui, Z. (2010). Nano-microparticles as immune adjuvants: correlating particle sizes and the resultant immune responses. Expert Rev. Vaccines 9 (9), 1095–1107. doi:10.1586/erv.10.89
Petersen, E., Koopmans, M., Go, U., Hamer, D. H., Petrosillo, N., Castelli, F., et al. (2020). Comparing SARS-CoV-2 with SARS-CoV and influenza pandemics. Lancet. Infect. Dis. 20 (9), e238–e244. doi:10.1016/S1473-3099(20)30484-9
Piplani, S., Singh, P. K., Winkler, D. A., and Petrovsky, N. (2021). In silico comparison of SARS-CoV-2 spike protein-ACE2 binding affinities across species and implications for virus origin. Sci. Rep. 11 (1), 13063. doi:10.1038/s41598-021-92388-5
Powell, A. E., Zhang, K., Sanyal, M., Tang, S., Weidenbacher, P. A., Li, S., et al. (2021). A single immunization with spike-functionalized ferritin vaccines elicits neutralizing antibody responses against SARS-CoV-2 in mice. ACS Central Sci. 7 (1), 183–199. doi:10.1021/acscentsci.0c01405
Qi, M., Zhang, X.-E., Sun, X., Zhang, X., Yao, Y., Liu, S., et al. (2018). Intranasal nanovaccine confers homo- and hetero-subtypic influenza protection. Small Weinheim Der Bergstrasse, Ger. 14 (13), e1703207. doi:10.1002/smll.201703207
Rahnama, S., Azimzadeh Irani, M., Amininasab, M., and Ejtehadi, M. R. (2021). S494 O-glycosylation site on the SARS-CoV-2 RBD affects the virus affinity to ACE2 and its infectivity; a molecular dynamics study. Sci. Rep. 11, 15162. doi:10.1038/s41598-021-94602-w
Rapin, N., Lund, O., Bernaschi, M., and Castiglione, F. (2010). Computational immunology meets bioinformatics: the use of prediction tools for molecular binding in the simulation of the immune system. PLOS ONE 5 (4), e9862. Accessed. doi:10.1371/journal.pone.0009862
Raposo, C. D., Canelas, A. B., and Barros, M. T. (2021). Human lectins, their carbohydrate affinities and where to find them. Biomolecules 11 (2), 188. doi:10.3390/biom11020188
Reis, C. A., Tauber, R., and Blanchard, V. (2021). Glycosylation is a key in SARS-CoV-2 infection. J. Mol. Med. Berlin, Ger. 99 (8), 1023–1031. doi:10.1007/s00109-021-02092-0
Reynisson, B., Alvarez, B., Paul, S., Peters, B., and Nielsen, M. (2020). NetMHCpan-4.1 and NetMHCIIpan-4.0: improved predictions of MHC antigen presentation by concurrent motif deconvolution and integration of MS MHC eluted ligand data. Nucleic acids Res. 48 (W1), W449–W454. doi:10.1093/nar/gkaa379
Saini, G., Dalal, V., Savita, B. K., Sharma, N., Kumar, P., and Sharma, A. K. (2019). Molecular docking and dynamic approach to virtual screen inhibitors against Esbp of Candidatus Liberibacter asiaticus. J. Mol. Graph. Model. 92, 329–340. doi:10.1016/j.jmgm.2019.08.012
Shajahan, A., Supekar, N. T., Gleinich, A. S., and Azadi, P. (2020). Deducing the N- and O-glycosylation profile of the spike protein of novel coronavirus SARSCoV-2. Glycobiology 30, 981–988. doi:10.1093/glycob/cwaa042
Shereen, M. A., Khan, S., Kazmi, A., Bashir, N., and Siddique, R. (2020). COVID-19 infection: origin, transmission, and characteristics of human coronaviruses. J. Adv. Res., 24, 91–98. doi:10.1016/j.jare.2020.03.005
Shiina, T., Hosomichi, K., Inoko, H., and Kulski, J. K. (2009). The HLA genomic loci map: expression, interaction, diversity and disease. J. Hum. Genet. 54 (1), 15–39. doi:10.1038/jhg.2008.5
Sliepen, K., Ozorowski, G., Burger, J. A., van Montfort, T., Stunnenberg, M., LaBranche, C., et al. (2015). Presenting native-like HIV-1 envelope trimers on ferritin nanoparticles improves their immunogenicity. Retrovirology 12, 82. doi:10.1186/s12977-015-0210-4
Song, W., Gui, M., Wang, X., and Xiang, Y. (2018). Cryo-EM structure of the SARS coronavirus spike glycoprotein in complex with its host cell receptor ACE2. PLoS Pathog. 14 (8), e1007236. doi:10.1371/journal.ppat.1007236
Stone, J. W., Thornburg, N. J., Blum, D. L., Kuhn, S. J., Wright, D. W., and Crowe, J. E. J. (2013). Gold nanorod vaccine for respiratory syncytial virus. Nanotechnology 24 (29), 295102. doi:10.1088/0957-4484/24/29/295102
Swanson, K. A., Rainho-Tomko, J. N., Williams, Z. P., Lanza, L., Peredelchuk, M., Kishko, M., et al. (2020). A respiratory syncytial virus (RSV) F protein nanoparticle vaccine focuses antibody responses to a conserved neutralization domain. Sci. Immunol. 5 (47), eaba6466. doi:10.1126/sciimmunol.aba6466
Tao, K., Tzou, P. L., Nouhin, J., Gupta, R. K., de Oliveira, T., Kosakovsky Pond, S. L., et al. (2021). The biological and clinical significance of emerging SARS-CoV-2 variants. Nat. Rev. Genet. 22 (12), 757–773. doi:10.1038/s41576-021-00408-x
Traherne, J. A. (2008). Human MHC architecture and evolution: implications for disease association studies. Int. J. immunogenetics 35 (3), 179–192. doi:10.1111/j.1744-313X.2008.00765.x
Tu, Y.-F., Chien, C.-S., Yarmishyn, A. A., Lin, Y.-Y., Luo, Y.-H., Lin, Y.-T., et al. (2020). A review of SARS-CoV-2 and the ongoing clinical trials. Int. J. Mol. Sci. 21, 2657. doi:10.3390/ijms21072657
Uchida, M., Kang, S., Reichhardt, C., Harlen, K., and Douglas, T. (2010). The ferritin superfamily: supramolecular templates for materials synthesis. Biochimica Biophysica Acta 1800 (8), 834–845. doi:10.1016/j.bbagen.2009.12.005
van Zundert, G. C. P., Rodrigues, J. P. G. L. M., Trellet, M., Schmitz, C., Kastritis, P. L., Karaca, E., et al. (2016). The HADDOCK2.2 web server: user-friendly integrative modeling of biomolecular complexes. J. Mol. Biol., 428(4), 720–725. doi:10.1016/j.jmb.2015.09.014
Vita, R., Mahajan, S., Overton, J. A., Dhanda, S. K., Martini, S., Cantrell, J. R., et al. (2019). The immune epitope Database (IEDB): 2018 update. Nucleic acids Res. 47 (D1), D339–D343. doi:10.1093/nar/gky1006
Walls, A. C., Park, Y.-J., Tortorici, M. A., Wall, A., McGuire, A. T., and Veesler, D. (2020). Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein. Cell 181 (2), 281–292. doi:10.1016/j.cell.2020.02.058
Wang, L. L., Xing, D., Le Van, A., Jerse, A. E., Wang, S., Kanekiyo, M., et al. (2017). Structure-based design of ferritin nanoparticle immunogens displaying antigenic loops of Neisseria gonorrhoeae. FEBS Open Bio 499 (7456), 102–106. doi:10.1038/nature12202
Wang, T., Zou, M., Jiang, H., Ji, Z., Gao, P., and Cheng, G. (2011). Synthesis of a novel kind of carbon nanoparticle with large mesopores and macropores and its application as an oral vaccine adjuvant. Eur. J. Pharm. Sci. Official J. Eur. Fed. Pharm. Sci. 44 (5), 653–659. doi:10.1016/j.ejps.2011.10.012
Wang, W., Huang, B., Zhu, Y., Tan, W., and Zhu, M. (2021). Ferritin nanoparticle-based SARS-CoV-2 RBD vaccine induces a persistent antibody response and long-term memory in mice. Cell. Mol. Immunol. 18 (3), 749–751. doi:10.1038/s41423-021-00643-6
Wang, Z., Li, C., Ellenburg, M., Soistman, E., Ruble, J., Wright, B., et al. (2006). Structure of human ferritin L chain. Acta Crystallogr. Sect. D. 62 (7), 800–806. doi:10.1107/S0907444906018294
Watanabe, Y., Bowden, T. A., Wilson, I. A., and Crispin, M. (2019). Exploitation of glycosylation in enveloped virus pathobiology. Biochimica biophysica acta. General Subj. 1863 (10), 1480–1497. doi:10.1016/j.bbagen.2019.05.012
Wolday, D., Fung, C. Y. J., Morgan, G., Casalino, S., Frangione, E., Taher, J., et al. (2023). HLA variation and SARS-CoV-2 specific antibody response. Viruses 15 (4), 906. doi:10.3390/v15040906
Wrapp, D., Wang, N., Corbett, K. S., Goldsmith, J. A., Hsieh, C.-L., Abiona, O., et al. (2020). Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation. Science 367 (6483), 1260 LP–1263. doi:10.1126/science.abb2507
Wu, D., Wu, T., Liu, Q., and Yang, Z. (2020). The SARS-CoV-2 outbreak: what we know. Int. J. Infect. Dis., 94, 44–48. doi:10.1016/j.ijid.2020.03.004
Wuertz, K. M., Barkei, E. K., Chen, W.-H., Martinez, E. J., Lakhal-Naouar, I., Jagodzinski, L. L., et al. (2021). A SARS-CoV-2 spike ferritin nanoparticle vaccine protects hamsters against Alpha and Beta virus variant challenge. Npj Vaccines 6 (1), 129. doi:10.1038/s41541-021-00392-7
Yang, B., Fan, J., Huang, J., Guo, E., Fu, Y., Liu, S., et al. (2021). Clinical and molecular characteristics of COVID-19 patients with persistent SARS-CoV-2 infection. Nat. Commun. 12 (1), 3501. doi:10.1038/s41467-021-23621-y
Yang, J., Wang, W., Chen, Z., Lu, S., Yang, F., Bi, Z., et al. (2020). A vaccine targeting the RBD of the S protein of SARS-CoV-2 induces protective immunity. Nature 586 (7830), 572–577. doi:10.1038/s41586-020-2599-8
Yao, D., Lao, F., Zhang, Z., Liu, Y., Cheng, J., Ding, F., et al. (2020). Human H-ferritin presenting RBM of spike glycoprotein as potential vaccine of SARS-CoV-2. BioRxiv 05.25, 115618. doi:10.1101/2020.05.25.115618
Yasunori, W., Daniel, W., Max, C., and Crispin, M. (2020). Site-specific glycan analysis of the SARS-CoV-2 spike. Sci. 80 369, 330–333. doi:10.1126/science.abb9983
Zhang, Y., and Orner, B. P. (2011). Self-assembly in the ferritin nano-cage protein superfamily. Int. J. Mol. Sci. 12 (8), 5406–5421. doi:10.3390/ijms12085406
Zhao, K., Chen, G., Shi, X.-M., Gao, T.-T., Li, W., Zhao, Y., et al. (2012). Preparation and efficacy of a live newcastle disease virus vaccine encapsulated in chitosan nanoparticles. PloS One 7 (12), e53314. doi:10.1371/journal.pone.0053314
Zhao, L., Seth, A., Wibowo, N., Zhao, C.-X., Mitter, N., Yu, C., et al. (2014). Nanoparticle vaccines. Vaccine 32 (3), 327–337. doi:10.1016/j.vaccine.2013.11.069
Zhao, P., Praissman, J. L., Grant, O. C., Cai, Y., Xiao, T., Rosenbalm, K. E., et al. (2020). Virus-receptor interactions of glycosylated SARS-CoV-2 spike and human ACE2 receptor. Cell Host Microbe, 28(4), 586–601. doi:10.1016/j.chom.2020.08.004
Zhu, F.-C., Zhang, J., Zhang, X.-F., Zhou, C., Wang, Z.-Z., Huang, S.-J., et al. (2010). Efficacy and safety of a recombinant hepatitis E vaccine in healthy adults: a large-scale, randomised, double-blind placebo-controlled, phase 3 trial. Lancet London, Engl. 376 (9744), 895–902. doi:10.1016/S0140-6736(10)61030-6
Keywords: nanoparticle vaccine, ferritin, SARS-CoV-2, monomeric spike, molecular docking, molecular dynamics simulation
Citation: Garmeh Motlagh F, Azimzadeh Irani M, Masoomi Nomandan SZ and Assadizadeh M (2024) Computational design and investigation of the monomeric spike SARS-CoV-2-ferritin nanocage vaccine stability and interactions. Front. Mol. Biosci. 11:1403635. doi: 10.3389/fmolb.2024.1403635
Received: 19 March 2024; Accepted: 23 May 2024;
Published: 12 June 2024.
Edited by:
Vikram Dalal, Washington University in St. Louis, United StatesReviewed by:
Padhmanand Sudhakar, Kumaraguru College of Technology, IndiaMaulikkumar P. Patel, Washington University in St. Louis, United States
Poonam Dhankhar, Cornell University, United States
Gunjan Saini, Purdue University, United States
Preeti Verma, University of Virginia, United States
Copyright © 2024 Garmeh Motlagh, Azimzadeh Irani, Masoomi Nomandan and Assadizadeh. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Maryam Azimzadeh Irani, bV9hemltemFkZWhAc2J1LmFjLmly