key: cord-0898453-bd1414f7
authors: Rantam, Fedik Abdul; Kharisma, Viol Dhea; Sumartono, Christrijogo; Nugraha, Jusak; Wijaya, Andi Yasmin; Susilowati, Helen; Kuncorojakti, Suryo; Nugraha, Alexander Patera
title: Molecular docking and dynamic simulation of conserved B cell epitope of SARS-CoV-2 glycoprotein Indonesian isolates: an immunoinformatic approach
date: 2021-08-16
journal: F1000Res
DOI: 10.12688/f1000research.54258.1
sha: aac2d5eaebb8913f2574bce22fb282e5ae113398
doc_id: 898453
cord_uid: bd1414f7

Background: An immunoinformatic approach may be useful to investigate the conserved region in the spike glycoprotein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Indonesia isolates. The aim of this study was to investigate Indonesian SARS-CoV-2 isolates based on B cell epitopes by targeting the conserved regions in the spike glycoprotein to trigger increased multi-variant virus neutralization and memory response for the development of vaccine seed candidates. Methods: SARS-CoV-2 spike glycoprotein gene sequences originating from Indonesia were compared with Wuhan (China), the United Kingdom, South Africa, India, the United States, and Brazil isolates obtained from the NCBI and GISAID databases. The recognition of antigens was carried out directly using B cells through the B cell receptor (BCR). An indirect B cell activation by Cluster of Differentiation (CD)4+ T cells and major histocompatibility complex (MHC)-II was predicted through the binding with human leukocyte antigen (HLA) based on IC (50 )value. In addition, vaccine allergenicity and toxicity were investigated. During the molecular complex examination, the 3D peptide structure was investigated and the lowest amount of energy formed when the vaccine candidate peptide bound to BCR and MHC-II was calculated. Results: As a result, the spike glycoprotein sequences of Indonesian SARS-CoV-2 isolates had conserved regions which were very similar to reference countries such as China, the United Kingdom, South Africa, India, the United States, and Brazil. Conclusion: It was predicted that the conserved regions could be identified as the epitope of B and T CD4+ cells that produced the peptides for vaccine candidate with antigenic, non-allergen, and non-toxic properties.

Coronavirus disease 2019 is caused by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). First identified in December 2019 in Wuhan, China, it then spread around the world and became a global health burden. COVID-19 patients may experience headaches, dizziness, coughing, and a loss of their sense of smell. 1 The virus is spread between humans when an infected individual expels droplets and airborne particles which come into contact with a healthy person. SARS-CoV-2 is a virus with RNA as its genetic material, an envelope, and in about 96% of COVID patients it is identified as BatCoV RaTG13. 2 4 The emergence of SARS-CoV-2 variants causes an increase in the rates of infection, reinfection, and avoidance of neutralization from the vaccine antibodies. 5 The COVID-19 pandemic has extended to all parts of the world including countries in the Southeast Asia region, especially Indonesia. The number of SARS-CoV-2 cases has increased rapidly in Indonesia, at around 1000 new cases per day since 2020. As of now, the distribution of SARS-CoV-2 vaccination in Indonesia still totally depends on foreign production. Thus, the development of domestic vaccines is necessary. 6 In addition, the domestic vaccine development uses various Indonesian SARS-CoV-2 isolates from the NCBI and GISAID databases. Therefore, the vaccines must include multi-variant protection, be immunogenic, and form immune memory. 1, 7 Moreover, the development of vaccines involves the use of SARS-CoV-2 glycoprotein, which plays an important role in the infection process inside the host. 7 Under normal conditions, spike glycoprotein will bind to the angiotensinconverting enzyme 2 (ACE-2) receptor as a viral attachment mechanism. This mechanism has inspired many studies that have developed vaccines focusing on the neutralization of antibodies to viral glycoproteins. 8 In addition, SARS-CoV-2 is able to mutate to form a new variant that is capable of changing the residues which can create the spike glycoprotein. This leads to the possibility that the spike glycoprotein could avoid being neutralized by antibodies produced by SARS-CoV-2 vaccine. 9 However, a previous study suggested that the development of the SARS-CoV-2 vaccine must target the conserved regions. Conserved regions can be found in every variant of SARS-CoV-2, and are the main target of vaccine design because they are able to trigger an increase in the antibody neutralization's coverage. 10 Additionally, conserved region-based vaccine development has been carried out on ZIKA, DENV and influenza viruses through an immunoinformatic approach, but is rarely found in SARS-CoV-2 vaccine development. 11 This approach was used in this study to determine the conserved region in the spike glycoprotein of Indonesian SARS-CoV-2 isolates, which can be used as vaccine seed candidates. BepiPred methods were employed in determining the B cell epitope to investigate the probability of the antigen region inducing the B cell recognition. 12 The candidate peptides that change the epitope of B cells would be expected to have antigenicity, similarity, and toxicity. Then, molecular docking and dynamic simulations would be carried out to determine the stability level of molecular complexes formed by peptides with B cell receptors (BCR). 13, 14 Furthermore, this study aimed to investigate Indonesian SARS-CoV-2 isolates based on their B cell epitopes by targeting the conserved regions in spike glycoprotein to trigger the increased multi-variant virus neutralization and memory response to develop seed candidates.

The SARS-CoV-2 spike glycoprotein gene sequences originating from Indonesian isolates and their references were obtained from NCBI and GISAID 44 databases with keyword 'SARS-COV-2 Spike Glycoprotein' and with filter 'no partial sequences'.

Then, protein alignment was performed using MEGA X (version 10.2.6, build 10210527, default settings GAP, (Open/ Extend: À2,90 & 0.00), Min Dag Length: 24, dan Max literation 16) to identify the position of the conserved region. 1, 15 The SWISS-MODEL (default settings) was applied in this study for the construction of the 3D spike glycoprotein SARS-CoV-2 structure. It utilized the construction method of homology modelling by considering sequences only where >20% of residues were homologous, and model validation through the Ramachandran plot which only considered sequences considered if >90% of residues were in the favored region. 16, 17 Additionally, comparisons of the simulation, staining selection, and structure of the SARS-CoV-2 query spike glycoprotein models with the reference sequences were completed using PyMol software version 2.5.0, 18 differentiating domain receptor, ligand and secondary protein structures by color.

The recognition of antigens is performed directly by B cells through the BCR. This consists of immunoglobulin and Cluster of Differentiation (CD) 79, which majorly contribute to signal transduction and activation of B cells in the immune response. 19, 20 The interaction between the variable region in the antigen-binding fragment (Fab) of the BCR and the epitope allow it to trigger the B cell recognition response, which can be predicted by epitope mapping analysis. 21, 22 IEDB Server Analysis Tools (default settings) were employed in this study to predict B cell epitopes based on the direct recognition mechanism through the BepiPred version 2.0 method. 8, 12, 23 CD4+ T cell epitope prediction Indirect B cell activation occurs via the peptide recognition pathway by CD4+ T cells through the major histocompatibility complex (MHC)-II. Investigating binding of a peptide with human leukocyte antigens (HLA) can predict their ability to trigger the formation of memory cells. 24, 25 MHC-II epitopes were predicted through IEDB's binding tool (default settings, locus allel MHC type II: Human, HLA-DR, with epitope length 15-mer, sorted peptide by: adjusted rank, output format: XHTML table), by binding with three HLAs (reference set) consisting of HLA-DRB1*01:01, HLA-DRB4*01:01, and HLA-DRB3*02:02. 26 Meanwhile, the peptide binding affinity for MHC-II encoded with specific HLA alleles was predicted based on the IC 50 value, which was determined through previous research by Zahroh et al. in 2011 45 ; high affinity was represented by IC 50 <50 nM, medium affinity by IC 50 <500 nM, and low affinity by IC 50 <5000 nM. 11 Vaccine profiling and peptide modeling Peptides with positive predictions (indicated by the IC 50 value of the vaccine candidate peptide of <50 nm and the percentage rank value of <0.5) for B and T CD4+ cell epitopes were tested through the VaxiJen v2.0 server (Setting as = "target organism": virus, threshold: default Output format: sequence Output) for the selection of antigenic properties. 27 Then, the level of allergenicity from the vaccine candidate peptides was predicted via the AlerCatPro 1.7 server (default settings). 28 The toxicity profile of the vaccine candidate peptides in this study was predicted using the ToxinPred server (threshold SVM: À0,1 with positive prediction (peptide but not toxin) indicated as peptide score below -0,1). 29 Lastly, the 3D peptide structure of the vaccine candidate was constructed via the PEP-FOLD 3.5 server (default settings) with de novo modeling methods, and a sample of the protein modeling result was saved in databank format. 30 Molecular docking and dynamics simulation The Cluspro 2.0 server (default settings) was applied in this study for molecular docking as the analysis aimed to determine the lowest amount of energy formed when the vaccine candidate peptide bound with BCR and MHC-II to create the molecular complex. The docking results were formed through superimposition. 13 Furthermore, the CABS-flex 2.0 server (default settings: protein rigidity: 1.0; protein restraints: ss2 3 3.8 8.0; global C-alpha restraints weight: 1.0; global side-chain restraints weight: 1.0; number of cycles: 50; cycles between trajectory: 50; temperature range: 1.40-1.40; RNG seed: 2465) was utilized in this study to validate the docking results through molecular dynamics simulation; when a molecular complex occurs, the dynamic molecular simulation can describe the fluctuations in the interaction of the amino acid residues that make up the peptide with the target protein domain (receptor). 14 

The SARS-CoV-2 spike glycoprotein sequence of Indonesian isolates that were successfully obtained from the NCBI database had the following IDs: MZ026853.1, MZ026854.1, and MZ026855.1. Additional samples were also generated from the GISAID database, which can be found in the underlying data section. 31 These gene samples were then translated via MEGA X (version 10.2.6) software to identify the amino acid residues that formed the SARS-CoV-2 spike glycoprotein. The conserved region was obtained through protein alignment with reference sequences MN908947. The 3D structures of spike glycoprotein isolates were obtained from SWISS-MODEL to determine the similarity to the structural templates and quality. The protein modeling results were displayed through the PyMol software with a cartoon display and underwent staining selection. When the results of the modeling of the spike glycoprotein SARS-CoV-2 query were completed, a comparison of the structure with the template was carried out through 3D alignment analysis. The resulting spike glycoprotein SARS-CoV-2 models of Indonesian isolates in this study had similarity values of 90-100% with the references in the database through the 3D alignment process (Figure 1 ). 32 Mapping of conserved B cell and T CD4+ epitope Indonesian isolate MZ026853.1 was chosen for SARS-COV-2 seed vaccine candidate as it is the dominant variant found in Indonesia. 33 Prediction of B cell epitope recognition was carried out on the template of Indonesian isolate MZ026853.1 using IEDB Analysis Tools with the BepiPred method. The peptides that showed a positive epitope prediction were selected based on the epitope's carried conserved region. The epitope/peptide of the B cell-based vaccine candidate refers to a score above the threshold of 0.35 (calculated based on the server's prediction results, which refers to the Hidden Markov Model and Trend Scale methods) while the T-cell-based vaccine candidate epitope/peptide refers to the IC 50 score >50 nM. The results showed that the position of the B cell epitope was identified on the SARS-CoV-2 spike glycoprotein of the Indonesian isolate and had a length of 10 to 11-mer ( Figure 2 ). 34 Furthermore, the results also revealed that all identified epitopes had a value above the threshold of 0.35, and based on the identification results, it exhibited all peptides that make up the epitope. In addition, inside the obtained B cells epitope, there were eight peptides with conserved regions from the total of 15 epitopes in this study ( Table 1) .

The prediction of CD4+ T cell epitope introduction was also carried out on the template of Indonesian isolate MZ026853.1 using IEDB Analysis Tools with the default method. In addition, the binding affinity of the vaccine candidate peptide with MHC-II encoded by specific HLA alleles was assessed based on a reference set consisting of HLA-DRB1*01:01, HLA-DRB4*01:01, and HLA-DRB3*02:02. The results showed there were 15 peptides categorized as CD4+ T cell epitopes because they had a high binding affinity for MHC-II encoded by specific HLA alleles. Moreover, there were 11 peptides that produced epitopes identified as having conserved regions and these were utilized for further analysis ( Table 2) .

There were 19 peptides which have B and T CD4+ cells with conserved regions. These were then tested for antigenicity considering the threshold value of 0.4 (determined by IEDB's default settings), which led to 13 peptides categorized as antigens because they had values higher than this. Furthermore, peptides with antigenic properties were tested for allergenicity by taking into account the threshold value of 20% (determined by IEDB's default settings), which resulted in all 13 obtained peptides being considered non-allergenic. Peptides categorized as non-allergenic then underwent toxicity prediction using the reference sequences with a threshold value of À0.8 (determined by IEDB's default settings), which incurred 10 non-toxic peptides (Table 3) . Thus, based on the analysis of vaccine property prediction, it was concluded that peptides B9, B12, L1, L2, L4, L5, L9, L11, L12 , and L13 could act as antigens and were not allergenic, as well as having non-toxic properties. These 10 peptides were then used for further analysis. 

Homology modeling is a method of determining the 3D structure of a query protein based on the availability of templates in the database. Its accuracy is around 90%, and the model is categorized as similar or homologous to the target protein if it has a modeling similarity value of >20%. 16, 17 The similarity levels are high if the superimposed query protein on the template is identified during 3D alignment. 18, 38 In this study, the spike glycoprotein structure of 20 Indonesian isolates of SARS-CoV-2 had high similarities to the isolate references such as those from China/Wuhan, the United Kingdom, South Africa, India, the United States, and Brazil with an average similarity level of 90% to 100% (See Figure 1 ). In addition, the similarities were also demonstrated between the secondary protein structure of the SARS-CoV-2 spike glycoprotein of the Indonesian isolates and the superimposed reference isolates.

The structure of the SARS-CoV-2 spike glycoprotein consists of S1 and S2 regions. S1 covers residue positions 13-685 and S2 covers positions 686-1273. 39 In the S1 region, there is a N-terminal domain as well as a receptor-binding domain (RBD) covering the sequence positions 319-541, including an integrin-binding motif (403-405), and receptor-binding motif (437-508). S2 consists of fusion peptides 1 and 2 (816-1202). 40 The conserved region is specific to the virus, is maintained in each generation of variants and has the potential to be a target for antibodies to neutralize of in vaccine design, since it is predicted to be able to protect against infections with a wide variety of variants. 8, 13, 41 Furthermore, the conserved region obtained from the identification results in this study was detected very high activity in almost all parts of the SARS-CoV-2 spike glycoprotein, especially in S1, RBD, and S2.

The recognition of B cell epitopes in the conserved region allows activation of B cells through direct and indirect pathways 19, 20 The The direct B cell activation pathway is triggered by peptide binding to the Fab region on the BCR, which allows signal transduction through CD79, and then triggers transcription factors involved in the activation, maturation, and proliferation of plasma cells, resulting in the production of IgM antibodies. 19, 46 On the other hand, the indirect pathway of activation occurs due to the representation of MHC-II by B cells, and its recognition by TCR on CD4+ T cells. 20, 47 Peptides with high binding affinity to MHC-II and encoded by specific HLA alleles are predicted to act as T cell epitopes. 11 The binding of the vaccine candidate peptide to the B cell receptor is identified by the lowest energy produced by the molecular docking complex. The lowest energy with a more negative value, is estimated to indicate the most stable and active biological response to protein initiation. 13 In this study, the activation response of B cells is predicted to be directly initiated by the molecular complex BCR_B12, and indirectly through HLA-DRB1*01:01_L5, HLA-DRB3*02:02_L9, and HLA-DRB4*01:01_L13 because they have the lowest negative energy values.

Finally, the binding stability of molecular complexes can be analyzed by molecular dynamics simulations by examining the level of fluctuation based on the RMSF value of each complex formed. 14 RMSF exhibits the movement of interacting atoms on the residues of peptides and proteins with a certain distance, while complex stability is achieved when the resulting distance has a value of ≤4Å. 48,49

The spike glycoprotein sequences of Indonesian SARS-CoV-2 isolates have a conserved region and are very similar to the reference isolates from China, the United Kingdom, South Africa, India, the United States, and Brazil. It is suggested that the conserved region of Indonesian SARS-CoV-2 isolates can be identified as an epitope of B and T CD4+ cells that produce the vaccine candidate peptides with antigenic, non-allergen, and non-toxic properties. This study recommends that peptides B12, L5, L9, and L13 can be used in the follow-up tests because they may trigger B cell activation responses via direct and indirect pathways. Lastly, the vaccine candidate is predicted to be able to initiate the production of specific antibody isotypes, which play a role in increasing the neutralization of multi-variant viruses and the formation of memory cells to prevent SARS-CoV-2 infection.

Underlying data The accession numbers for GISAID and NCBI can be found in the acknowledgments 

Faculty of Dentistry, Ibb University, Ibb, Yemen

The study entitled "Molecular docking and dynamic simulation of conserved B cell epitope of SARS-CoV-2 glycoprotein Indonesian isolates: an immunoinformatic approach" was an original article that is well organized and valuable. It checked the possibility of investigating SARS-CoV-2 in Indonesian isolates based on B cell epitopes by marking the conserved regions in the spike glycoprotein to trigger the elevated multi-variant virus neutralization and memory response for the development of vaccine seed candidates.

Methods included a comparison between SARS-CoV-2 spike glycoprotein gene sequences originating from Indonesia with Wuhan (China), the United Kingdom, South Africa, India, the United States, and Brazil isolates collected from the NCBI and GISAID databases. B cell receptor (BCR) was used to recognize antigens on the B cells. Indirect B cell activation by Cluster of Differentiation (CD)4+ T cells and major histocompatibility complex (MHC)-II was predicted through the binding with human leukocyte antigen (HLA) based on IC50 value. Vaccine allergenicity and toxicity were also examined.

Results showed similarities between Indonesian conserved regions and the compared mentioned countries, thus the predicted expectations that the conserved regions could be identified as the epitope of B and T CD4+ cells that produced the peptides for vaccine candidates with antigenic, non-allergen, and non-toxic properties.

I feel satisfied with the article's frame, methodology and results.

Just a few English corrections are to be rechecked.

Is the study design appropriate and is the work technically sound?

Distribution of COVID-19 and phylogenetic tree construction of SARS-CoV-2 in Indonesia

Retrieval and Investigation of Data on SARS-CoV-2 and COVID-19 Using Bioinformatics Approach

PubMed Abstract|Publisher Full Text

Analysis of the Genomic Distance Between Bat Coronavirus RaTG13 and SARS-CoV-2 Reveals Multiple Origins of COVID-19

PubMed Abstract|Publisher Full Text|Free Full Text

Emerging SARS-CoV-2 Variants and Impact in Global Vaccination Programs against SARS-CoV-2/COVID-19

The great escape? SARS-CoV-2 variants evading neutralizing responses

Acceptance of a COVID-19

PubMed Abstract|Publisher Full Text|Free Full Text

COVID-19 pandemic in Indonesia: Situation and challenges of rehabilitation medicine in Indonesia

Immunobioinformatics analysis and phylogenetic tree construction of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in Indonesia: spike glycoprotein gene

The British variant of the new coronavirus-19 (Sars-Cov-2) should not create a vaccine problem

PubMed Abstract|Publisher Full Text

Coronavirus epitope prediction from highly conserved region of spike protein

Immunoinformatics Approach in Designing Epitopebased Vaccine Against Meningitis-inducing Bacteria (Streptococcus pneumoniae, Neisseria meningitidis, and Haemophilus influenzae Type b). Drug target insights

Analysis of the Genome Sequence and Prediction of B-Cell Epitopes of the Envelope Protein of Middle East Respiratory Syndrome-Coronavirus

PubMed Abstract|Publisher Full Text|Free Full Text

The ClusPro web server for protein-protein docking

PubMed Abstract|Publisher Full Text|Free Full Text

CABS-flex 2.0: a web server for fast simulations of flexibility of protein structures

PubMed Abstract|Publisher Full Text|Free Full Text

Construction of epitope-based peptide vaccine against SARS-CoV-2: Immunoinformatics study

SWISS-MODEL: homology modelling of protein structures and complexes

PubMed Abstract|Publisher Full Text|Free Full Text

The exploration of medicinal plants' phytochemical compounds as potential inhibitor against human α-3 nicotinic acetylcholine receptors: The insight from computational study

PyMod: sequence similarity searches, multiple sequence-structure alignments, and homology modeling within PyMOL

PubMed Abstract|Publisher Full Text|Free Full Text

B-cell receptor: from resting state to activate

PubMed Abstract|Publisher Full Text|Free Full Text

Expression of the B-cell receptor component CD79a on immature myeloid cells contributes to their tumor promoting effects

An Introduction to B-Cell Epitope Mapping and In Silico Epitope Prediction

PubMed Abstract|Publisher Full Text|Free Full Text

Functional mapping of B-cell linear epitopes of SARS-CoV-2 in COVID-19 convalescent population

PubMed Abstract|Publisher Full Text|Free Full Text

Immunobioinformatics of Rabies Virus in Various Countries of Asia: Glycoprotein Gene

Vaccine design for H5N1 based on B-and T-cell epitope predictions

Immunoinformatic Analysis of T-and B-Cell Epitopes for SARS-CoV-2 Vaccine Design

PubMed Abstract|Publisher Full Text|Free Full Text

Human Leukocyte Antigen Susceptibility Map for Severe Acute Respiratory Syndrome Coronavirus 2

PubMed Abstract|Publisher Full Text|Free Full Text

VaxiJen Dataset of Bacterial Immunogens: An Update. Current Computer-Aided Drug Design

PubMed Abstract|Publisher Full Text

AllerCatPro-prediction of protein allergenicity potential from the protein sequence

PubMed Abstract|Publisher Full Text|Free Full Text

Peptide toxicity prediction

PubMed Abstract|Publisher Full Text

PEP-FOLD3: faster de novo structure prediction for linear peptides in solution and in complex

PubMed Abstract|Publisher Full Text|Free Full Text

Publisher Full Text

BioSuper: a web tool for the superimposition of biomolecules and assemblies with rotational symmetry

Structural basis for the recognition of SARS-CoV-2 by full-length human ACE2

PubMed Abstract|Publisher Full Text|Free Full Text

Characterization of the receptor-binding domain (RBD) of 2019 novel coronavirus: implication for development of RBD protein as a viral attachment inhibitor and vaccine

PubMed Abstract|Publisher Full Text|Free Full Text

Severe Acute Respiratory Syndrome Coronavirus 2 Spike Protein Based Novel Epitopes Induce Potent Immune Responses in vivo and Inhibit Viral Replication in vitro

IgD B cell receptors differentially respond to endogenous antigens and control B cell fate. eLife

PubMed Abstract|Publisher Full Text|Free Full Text

Role of the IgM Fc Receptor in Immunity and Tolerance

IgH chain class switch recombination: mechanism and regulation

PubMed Abstract|Publisher Full Text|Free Full Text

Distinct T helper cell dependence of memory B-cell proliferation versus plasma cell differentiation

PubMed Abstract|Publisher Full Text|Free Full Text

Structural and functional studies of Igalphabeta and its assembly with the B cell antigen receptor

PubMed Abstract|Publisher Full Text|Free Full Text

The ins and outs of MHC class II-mediated antigen processing and presentation

PubMed Abstract|Publisher Full Text|Free Full Text

Computational studies of protein-protein dissociation by statistical potential and coarse-grained simulations: a case study on interactions between colicin E9 endonuclease and immunity proteins

PubMed Abstract|Publisher Full Text|Free Full Text

Identification of a novel inhibitor of SARS-CoV-2 3CL-PRO through virtual screening and molecular dynamics simulation

Data, disease and diplomacy: GISAID's innovative contribution to global health

PubMed Abstract|Publisher Full Text|Free Full Text

Immunoinformatics Approach in Designing Epitope-based Vaccine Against Meningitis-inducing Bacteria (Streptococcus pneumoniae, Neisseria meningitidis, and Haemophilus influenzaeType b). Drug Target Insights

PubMed Abstract|Publisher Full Text|Free Full Text

We gratefully acknowledge the following Authors from the Originating laboratories responsible for obtaining the specimens and the Submitting laboratories where genetic sequence data were generated and shared via the GISAID Initiative, on which this research is based.

Are all the source data underlying the results available to ensure full reproducibility? Yes

Competing Interests: No competing interests were disclosed.

Clinical and lab studies in medicine and dentistry. Systematic reviews and meta-analyses as well.I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard. The current study, "Molecular docking and dynamic simulation of conserved B cell epitope of SARS-CoV-2 glycoprotein Indonesian isolates: an immunoinformatic approach," which employs an Immunoformatics approach, appears to be a high-quality piece of work with promising results. The use of bioinformatics tools reduces the cost and time required to predict potential epitopes. Furthermore, the current study is appealing as a first step in developing a vaccine because it focuses on the local Indonesian sequences of SARS-CoV-2.I recommend this study strongly for indexing due to the following reasons:The current study is well-presented, and the most recent literature is cited.

The study design was found sound and used immunoinformatic approaches for B-cell, T-cell ○ prediction, 3D structure development, and molecular docking are technically sound approaches.The methodology and the discussion section were both clear, explanatory, and technically sound.

The study's abstract and conclusion were found to be brief and concise.

Is the study design appropriate and is the work technically sound? Yes

If applicable, is the statistical analysis and its interpretation appropriate? I cannot comment. A qualified statistician is required.

Competing Interests: No competing interests were disclosed.

We confirm that we have read this submission and believe that we have an appropriate level of expertise to confirm that it is of an acceptable scientific standard. However, as preliminary result the author supposed to propose future direction of the next research for vaccine production.○ Molecular docking and dynamic simulation results should be confirm again using in vitro or in vivo models. It could be proven or not according to this result. It is far for realistic the candidate of molecule and the rapid test really occur in vivo.

Is the study design appropriate and is the work technically sound? Yes

Are all the source data underlying the results available to ensure full reproducibility? Yes

Competing Interests: No competing interests were disclosed.Reviewer Expertise: Cardiovascular Medicine, Molecular Cardiology, Interventional Cardiology, Pharmacogenetics I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard.The benefits of publishing with F1000Research:Your article is published within days, with no editorial bias • You can publish traditional articles, null/negative results, case reports, data notes and more • The peer review process is transparent and collaborative • Your article is indexed in PubMed after passing peer review • Dedicated customer support at every stage • For pre-submission enquiries, contact research@f1000.com