A multi-epitope pan-betacoronavirus vaccine construct predicted to induce broad-spectrum and durable immune responses: an immunoinformatics approach

Introduction: Recurrent zoonotic spillovers and continuous antigenic evolution among betacoronaviruses, including SARS-CoV, MERS-CoV, and SARS-CoV-2, highlight the urgent need for a broad-spectrum vaccine capable of eliciting cross-protective immunity. Conventional vaccines, although effective against specific strains, may be limited by antigenic mismatch and waning immunity. This study aimed to design a multi-epitope pan-betacoronavirus vaccine targeting conserved regions within the receptor-binding domain (RBD) using an integrated immunoinformatics and reverse vaccinology framework. Methods: Cytotoxic T-lymphocyte (CTL), helper T-lymphocyte (HTL), and linear B-cell epitopes were predicted and screened for antigenicity, allergenicity, toxicity, and non-homology to host proteins. Selected epitopes were assembled into a 285-amino acid multi-epitope construct using optimized linkers (AAY, GPGPG, EAAAK, and GGGGS) and human β-defensin 3 as an adjuvant. Structural modeling and refinement were performed to generate a three-dimensional vaccine model, followed by molecular docking with a B-cell receptor (BCR) Fab model using ClusPro. Molecular dynamics simulations were conducted to evaluate structural stability, and immune responses were assessed through computational immune simulation. Results: The refined vaccine construct produced a stable structural model with a C-score of -3.60. Molecular docking identified a highly ranked complex from a well-populated cluster (Cluster 1; 49 members) with a Lowest Energy score of -865.9, indicating favorable interface complementarity under the docking scoring function. Molecular dynamics simulation over 100 ns supported the structural integrity and dynamic stability of the complex, with minimal backbone deviation and sustained intermolecular interactions. Immune simulations predicted coordinated humoral and cellular responses following a simulated prime-boost regimen, including increased antibody titers, elevated IL-2 and IFN-γ levels, and sustained memory B- and T-cell populations. The selected epitope set showed an estimated global HLA population coverage of 93.28%. Conclusion: models.