Biological and Medicinal Chemistry

Whole Genome Sequence Analysis and Homology Modelling of a 3C Like Peptidase and a Non-Structural Protein 3 of the SARS-CoV-2 Shows Protein Ligand Interaction with an Aza-Peptide and a Noncovalent Lead Inhibitor with Possible Antiviral Properties

Arun Shanker Central Research Institute for Dryland Agriculture
,

Abstract

The family of viruses belonging to Coronaviridae mainly consist of virulent pathogens that have a zoonotic property, Severe Acute Respiratory Syndrome (SARS-CoV) and Middle East Respiratory Syndrome (MERS-CoV) of this family have emerged before and now the SARS-CoV-2 has emerged in China. Characterization of spike glycoproteins, polyproteins and other viral proteins from viruses are important for vaccine development. Homology modelling of these proteins with known templates offers the opportunity to discover ligand binding sites and explore the possible antiviral properties of these protein ligand complexes. In this study we did a complete bioinformatic analysis, sequence alignment, comparison of multiple sequences and homology modelling of the SARS-CoV-2 whole genome sequences, the spike protein and the polyproteins for homology with known proteins, we also analysed receptor binding sites in these models for possible binding with ligands that exhibit antiviral properties. Our results showed that the tertiary structure of the polyprotein isolate SARS-CoV-2_HKU-SZ-001_2020 had 98.94 percent identity with SARS-Coronavirus NSP12 bound to NSP7 and NSP8 co-factors. Our results indicate that a part of the viral genome (residues 3268 -3573 in Frame 2 with 306 amino acids) of the SARS-CoV-2 virus isolate Wuhan-Hu-1 (Genbank Accession Number MN908947.3) when modelled with template 2a5i of the PDB database had 96 percent identity with a 3C like peptidase of SARS-CoV which has ability to bind with Aza-Peptide Epoxide (APE) which is known for irreversible inhibition of SARS-CoV main peptidase. Docking profile with 9 different conformations of the ligand with the protein model using Autodock Vina showed an affinity of -7.1 Kcal/mol. This region was conserved in 831 genomes of SARS-CoV-2. The part of the genome (residues 1568-1882 in Frame 2 with 315 amino acids) when modelled with template 3e9s of the PDB database had 82 percent identity with a papain-like protease/deubiquitinase which when complexed with ligand GRL0617 acts as inhibitor which can block SARS-CoV replication. Docking profile with 9 different conformation of the ligand with the protein model using Autodock Vina showed an affinity of -7.9 Kcal/mol. This region was conserved in 831 genomes of SARS-CoV-2. It is possible that these ligands can be antivirals of SARS-CoV-2.

Version notes

This is the ninth version of the MS, docking was done with Autodock Vina and the results of the docking was added in this version. Docking study results were added in the MS in Materials and Methods and Results and Discussion sections. Samriddhi Gupta who did the docking study was added in the author list in this version.

Content

Thumbnail image of Coronavirus_MS_16.pdf
download asset Coronavirus_MS_16.pdf 2 MB [opens in a new tab]

Supplementary material

Thumbnail image of Tables.docx
download asset Tables.docx 0.02 MB [opens in a new tab]
Tables
Thumbnail image of Supplementary Tables.docx
download asset Supplementary Tables.docx 0.04 MB [opens in a new tab]
Supplementary Tables
Thumbnail image of Supplementaty Fig 1.pdf
download asset Supplementaty Fig 1.pdf 0.15 MB [opens in a new tab]
Supplementaty Fig 1
Thumbnail image of Figure 11.jpg
download asset Figure 11.jpg 0.16 MB [opens in a new tab]
Figure 10
Thumbnail image of Coronavirus_MS_16.docx
download asset Coronavirus_MS_16.docx 1 MB [opens in a new tab]
Coronavirus MS 16
Thumbnail image of Figure 2.jpg
download asset Figure 2.jpg 0.04 MB [opens in a new tab]
Figure 1
Thumbnail image of Figure 3.jpg
download asset Figure 3.jpg 0.04 MB [opens in a new tab]
Figure 2
Thumbnail image of Figure 4.jpg
download asset Figure 4.jpg 0.93 MB [opens in a new tab]
Figure 3
Thumbnail image of Figure 5.jpg
download asset Figure 5.jpg 0.68 MB [opens in a new tab]
Figure 4
Thumbnail image of alignment_3C like Protein with 2a5i.jpg
download asset alignment_3C like Protein with 2a5i.jpg 0.05 MB [opens in a new tab]
Figure 5
Thumbnail image of alignment_Structural Protein_3e9s.jpg
download asset alignment_Structural Protein_3e9s.jpg 0.05 MB [opens in a new tab]
Figure 6
Thumbnail image of 2a5i_3C_protein.jpg
download asset 2a5i_3C_protein.jpg 0.08 MB [opens in a new tab]
Figure 7
Thumbnail image of Figure 9.jpg
download asset Figure 9.jpg 0.01 MB [opens in a new tab]
Figure 8
Thumbnail image of 3e9s_Structural Protein.jpg
download asset 3e9s_Structural Protein.jpg 0.08 MB [opens in a new tab]
Figure 9
Thumbnail image of Figure 12.jpg
download asset Figure 12.jpg 0.12 MB [opens in a new tab]
Figure 11
Thumbnail image of Figure 12 a.jpg
download asset Figure 12 a.jpg 0.04 MB [opens in a new tab]
Figure 12 a
Thumbnail image of Figure 12 a.jpg
download asset Figure 12 a.jpg 0.04 MB [opens in a new tab]
Figure 12 a
Thumbnail image of Figure 12 b.jpg
download asset Figure 12 b.jpg 0.04 MB [opens in a new tab]
Figure 12 b
Thumbnail image of Figure 13 a.jpg
download asset Figure 13 a.jpg 0.03 MB [opens in a new tab]
Figure 13 a
Thumbnail image of Figure 13 b.jpg
download asset Figure 13 b.jpg 0.04 MB [opens in a new tab]
Figure 13 b