Affiliations 

  • 1 Department of Pharmaceutical Chemistry, Kulliyyah of Pharmacy, International Islamic University Malaysia, Indera Mahkota, 25200 Kuantan, Pahang, Malaysia
  • 2 Department of Pharmaceutical Chemistry, Kulliyyah of Pharmacy, International Islamic University Malaysia, Indera Mahkota, 25200 Kuantan, Pahang, Malaysia; Drug Discovery and Synthetic Chemistry Research Group, Kulliyyah of Pharmacy, International Islamic University Malaysia, Indera Mahkota, 25200 Kuantan, Pahang, Malaysia. Electronic address: zalikha@iium.edu.my
  • 3 Inborn Error of Metabolism and Genetic Unit, Nutrition, Metabolism and Cardiovascular Research Centre, Institute for Medical Research, National Institutes of Health, Ministry of Health Malaysia, Section U13 Setia Alam, 40170 Shah Alam, Selangor, Malaysia
Comput Biol Chem, 2023 Dec;107:107962.
PMID: 37847978 DOI: 10.1016/j.compbiolchem.2023.107962

Abstract

Protein arginine deiminase IV (PAD4) is a potential target for diseases including rheumatoid arthritis and cancers. Currently, GSK199 is a potent, selective yet reversible PAD4 inhibitor. Its derivative, GSK106, on the other hand, was reported as an inactive compound when tested against PAD4 assay. Although they had similar skeleton, their impact towards PAD4 structural and flexibility is unknown. In order to fill the research gap, the impact of GSK199 and GSK106 binding towards PAD4 stability and flexibility is investigated via a combination of computational methods. Molecular docking indicates that GSK199 and GSK106 are capable to bind at PAD4 pocket by using its back door with -10.6 kcal/mol and -9.6 kcal/mol, respectively. The simulations of both complexes were stable throughout 100 ns. The structure of PAD4 exhibited a tighter packing in the presence of GSK106 compared to GSK199. The RMSF analysis demonstrates significant changes between the PAD4-GSK199 and PAD4-GSK106 simulations in the regions containing residues 136, 160, 220, 438, and 606. The Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) analysis shows a marked difference in binding free energies, with -11.339 kcal/mol for the PAD4-GSK199 complex and 1.063 kcal/mol for the PAD4-GSK106 complex. The hydrogen bond analysis revealed that the GSK199 and GSK106 binding to PAD4 are assisted by six hydrogen bonds and three hydrogen bonds, respectively. The visualisation of the MD simulations revealed that GSK199 remained in the PAD4 pocket, whereas GSK106 shifted away from the catalytic site. Meanwhile, molecular dockings of benzoyl arginine amide (BAEE) substrate have shown that BAEE is able to bind to PAD4 catalytic site when GSK106 was present but not when GSK199 occupied the site. Overall, combination of computational approaches successfully described the behaviour of binding pocket of PAD4 structure in the presence of the active and inactive compounds.

* Title and MeSH Headings from MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine.