An immunoinformatics approach to study the epitopes of SARS-CoV-2 helicase, Nsp13

Vaccines are administered worldwide to control on-going coronavirus disease-19 (COVID-19) pandemic caused by SARS-CoV-2. Vaccine efficacy is largely contributed by the epitopes present on the viral proteins and their alteration might help emerging variants to escape host immune surveillance. Therefore, this study was designed to study SARS-CoV-2 Nsp13 protein, its epitopes and evolution.

Clustal Omega was used to identify mutations in Nsp13 protein. Secondary structure and disorder score was predicted by CFSSP and PONDR-VSL2 webservers. Protein stability was predicted by DynaMut webserver. B cell epitopes were predicted by IEDB DiscoTope 2.0 tools and their 3D structures were represented by discovery studio. Antigenicity and allergenicity of epitopes were predicted by Vaxijen2.0 and AllergenFPv.1.0. Physiochemical properties of epitopes were predicted by Toxinpred, HLP webserver tool.

Our data revealed 182 mutations in Nsp13 among Indian SARS-CoV-2 isolates, which were characterised by secondary structure and per-residue disorderness, stability and dynamicity predictions. To correlate the functional impact of these mutations, we characterised the most prominent B cell and T cell epitopes contributed by Nsp13. Our data revealed twenty-one epitopes, which exhibited antigenicity, stability and interactions with MHC class-I and class-II molecules. Subsequently, the physiochemical properties of these epitopes were analysed. Furthermore, eighteen mutations reside in these Nsp13 epitopes.

We report appearance of eighteen mutations in the predicted twenty-one epitopes of Nsp13. Among these, at least seven epitopes closely matches with the functionally validated epitopes. Altogether, our study shows the pattern of evolution of Nsp13 epitopes and their probable implications.