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

## TL;DR

This paper proposes a new vaccine design using computer modeling to target multiple betacoronaviruses, including SARS-CoV-2, with a focus on long-lasting and broad immune protection.

## Contribution

A novel multi-epitope vaccine construct targeting conserved regions of betacoronaviruses, predicted to induce durable and broad-spectrum immune responses.

## Key findings

- The vaccine construct achieved 93.28% global HLA population coverage and showed favorable structural stability.
- Molecular docking and dynamics simulations confirmed strong BCR binding and complex stability.
- Immune simulations predicted elevated antibody and cytokine responses with sustained memory cell populations.

## Abstract

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.

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.

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%.

This study identifies a promising in silico multi-epitope RBD-based pan-betacoronavirus vaccine candidate with predicted broad HLA population coverage and favorable structural stability. These findings provide a computational basis for subsequent experimental validation of the construct’s immunogenicity, safety, and potential cross-protective capacity in relevant in vitro and in vivo models.