Integration of virtual screening and computational simulation identifies photodynamic therapeutics against human Protoporphyrinogen Oxidase IX (hPPO)

Abstract

Photodynamic therapy (PDT) is a rapidly evolving area of cancer management against solid tumors. PDT is either administrated by injecting photosensitizer (porphyrins) or by accumulation of intracellular protoporphyrin IX via the inhibition of human Protoporphyrinogen Oxidase IX (hPPO). In this study, novel inhibitors of hPPO have been investigated by integrating virtual screening, molecular docking, and molecular dynamics (MD) simulation. A ligand-based pharmacophore was generated from a training set of 22 inhibitors of hPPO. The selected pharmacophore had four chemical features including three hydrogen bond acceptors and one hydrophobic. The pharmacophore was characterized by highest correlation coefficient of 0.96, cost difference of 53.20, and lowest root mean square deviation of 0.73. The resultant pharmacophore was validated by Fischer’s Randomization and Test Set Validation methods. The validated pharmacophore was used as a 3D query to screen chemical databases including NCI, Asinex, Chembridge, and Maybridge. The screening of chemical databases and the subsequent application of Lipinski’s Rule of Five, and ADMET Assessment Test, retrieved 1176 drug-like compounds. The drug-like compounds were subjected to molecular docking studies in the active site of hPPO to eliminate false positive hits and to elucidate their true binding orientation. Top three candidate molecules with high docking scores and hydrogen bond interactions with catalytic active residues were selected as best candidate inhibitors against hPPO. The binding stability of selected candidate inhibitors was evaluated by MD simulation. The MD simulation of hits portrayed strong hydrogen bonds and key hydrophobic interactions with catalytic active residues of hPPO including R59, R97, G159, G332 and flavin moiety of FAD (coenzyme of hPPO). Our study predicts three hit compounds against hPPO, which could possibly accumulate high concentration of protoporphyrinogen-IX, and thereby acting as an intracellular photosensitizer against tumor cells through photodynamic therapy. Keywords: Photodynamic therapy (PDT), hPPO Inhibition, Virtual screening, Pharmacophore modeling, Molecular docking simulation, Molecular dynamics (MD) simulation