Structure-based virtual screening suggests inhibitors of 3-Chymotrypsin-Like Protease of SARS-CoV-2 from Vernonia amygdalina and Occinum gratissimum

Gyebi, Gideon A.; Elfiky, Abdo A.; Ogunyemi, Oludare M.; Ibrahim, Ibrahim M.; Adegunloye, Adegbenro P.; Adebayo, Joseph O.; Olaiya, Charles O.; Ocheje, Joshua; Fabusiwa, Mercy M.

Antiviral culinary plants are potential bioresources for preventive nutraceuticals and/or antiviral drugs in COVID-19. Structure-based virtual screening was undertaken to screen 173 compounds previously reported from Vernonia amygdalina and Occinum gratissimum for direct interaction with the active site of the 3-Chymotrypsin-Like Protease (3CL(pro)) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Based on docking scores and comparison with reference inhibitors, a hit-list of 10 top phytocompounds was defined, which also had strong interactions with the catalytic centre of 3CL(pro) from three related strains of coronavirus (SARS-CoV, MERS-CoV, HKU4). Among these, six compounds (neoandrographolide, vernolide, isorhamnetin, chicoric acid, luteolin, and myricetin) exhibited the highest binding tendencies to the equilibrated conformers of SARS-CoV-2 3CL(pro) in an in-depth docking analysis to 5 different representative conformations from the cluster analysis of the molecular dynamics simulation (MDS) trajectories of the protein. In silico drug-likeness analyses revealed two drug-like terpenoids viz: neoandrographolide and vernolide as promising inhibitors of SARS-CoV-2 3CL(pro). These structures were accommodated within the substrate-binding pocket; and interacted with the catalytic dyad (Cys(145) and His(41)), the oxyanion loop (residues 138–145), and the S1/S2 sub-sites of the enzyme active site through the formation of an array of hydrogen bonds and hydrophobic interactions. Molecular dynamics simulation and binding free energy calculation revealed that the terpenoid-enzyme complexes exhibit strong interactions and structural stability. Therefore, these compounds may stabilize the conformation of the flexible oxyanion loop; and thereby interfere with the tetrahedral oxyanion intermediate formation during the proteolytic activity of the enzyme.