Polio viral proteinase 2A performs several essential functions in genome replication. Its inhibition prevents viral replication, thus making it an excellent substrate for drug development. In this study, the three-dimensional structure of 2A protease was determined and optimized by homology modelling. To predict the molecular basis of the interaction of small molecular agonists, docking simulations were performed on a structurally diverse dataset of poliovirus 2A protease (PV2Apr°) inhibitors. Docking results were employed to identify high risk missense mutations that are highly damaging to the structure, as well as the function, of the protease. Intrinsic disorder regions (IDRs), drug binding sites (DBS), and protein stability changes upon mutations were also identified among them. Our results demonstrated dominant roles for Lys 15, His 20, Cys 55, Cys 57, Cys 64, Asp 108, Cys 109 and Gly 110, indicating the presence of various important drug binding sites of the protein. Upon subjecting these sites to single-nucleotide polymorphism (SNP) analysis, we observed that out of 155 high risk SNPs, 139 residues decrease the protein stability. We conclude that these missense mutations can affect the functionality of the 2A protease, and that identified protein binding sites can be directed for the attachment and inhibition of the target proteins.
Keywords: computational analysis; docking; drug binding sites; non-synonymous SNPs (nsSNPs); poliovirus 2A protease (PV2Apr°); single-nucleotide polymorphisms (SNPs); stability changes.