Computational construction of a glycoprotein multi-epitope subunit vaccine candidate for old and new South-African SARS-CoV-2 virus strains

Abstract

The discovery of a new SARS-CoV-2 virus strain in South Africa presents a major public health threat, therefore contributing to increased infections and transmission rates during the second wave of the global pandemic. This study lays the groundwork for the development of a novel subunit vaccine candidate from the circulating strains of South African SARS-CoV-2 and provides an understanding of the molecular epidemiological trend of the circulating strains. A total of 475 whole-genome nucleotide sequences from South Africa submitted between December 1, 2020 and February 15, 2021 available at the GISAID database were retrieved based on its size, coverage level and hosts. To obtain the distribution of the clades and lineages of South African SARS-CoV-2 circulating strains, the metadata of the sequence retrieved were subjected to an epidemiological analysis. There was a prediction of the cytotoxic T lymphocytes (CTL), Helper T cells (HTL) and B-cell epitopes. Furthermore, there was allergenicity, antigenicity and toxicity predictions on the epitopes. The analysis of the physicochemical properties of the vaccine construct was performed; the secondary structure, tertiary structure and B-cell 3D conformational structure of the vaccine construct were predicted. Also, molecular binding simulations and dynamics simulations were adopted in the prediction of the vaccine construct's stability and binding affinity with TLRs. Result obtained from the metadata analysis indicated lineage B.1.351 to be in higher circulation among various circulating strains of SARS-CoV-2 in South Africa and GH has the highest number of circulating clades. The construct of the novel vaccine was antigenic, non-allergenic and non-toxic. The Instability index (II) score and aliphatic index were estimated as 41.74 and 78.72 respectively. The computed half-life in mammalian reticulocytes was 4.4 h in vitro, for yeast and in E. coli was >20 h and >10 h in vivo respectively. The grand average of hydropathicity (GRAVY) score is estimated to be −0.129, signifying the hydrophilic nature of the protein. The molecular docking indicates that the vaccine construct has a high binding affinity towards the TLRs with TLR 3 having the highest binding energy (−1203.2 kcal/mol) and TLR 9 with the lowest (−1559.5 kcal/mol). These results show that the vaccine construct is promising and should be evaluated using animal model.