Stepwise development of structure–activity relationship of diverse PARP-1 inhibitors through comparative and validated in silico modeling techniques and molecular dynamics simulation

Abstract

Inhibitors of poly (ADP-ribose) polymerase-1 (PARP-1) enzyme are useful for the treatment of various diseases including cancer. Comparative in silico studies were performed on different ligand-based (2D-QSAR, Kernel-based partial least square (KPLS) analysis, Pharmacophore Search Engine (PHASE) pharmacophore mapping), and structure-based (molecular docking, MM-GBSA analyses, Gaussian-based 3D-QSAR analyses on docked poses) modeling techniques to explore the structure–activity relationship of a diverse set of PARP-1 inhibitors. Two-dimensional (2D)-QSAR highlighted the importance of charge topological index (JGI7), fractional polar surface area (Jurs_FPSA3), and connectivity index (CIC2) along with different molecular fragments. Favorable and unfavorable fingerprints were demonstrated in KPLS analysis, whereas important pharmacophore features (one acceptor, one donor, and two ring aromatic) along with favorable and unfavorable field effects were demonstrated in PHASE-based pharmacophore model. MM-GBSA analyses revealed significance of different polar, non-polar, and solvation energies. Docking-based alignment of ligands was used to perform Gaussian-based 3D-QSAR study that further demonstrated importance of different field effects. Overall, it was found that polar interactions (hydrogen bonding, bridged hydrogen bonding, and pi–cation) play major roles for higher activity. Steric groups increase the total contact surface area but it should have higher fractional polar surface area to adjust solvation energy. Structure-based pharmacophore mapping spotted the positive ionizable feature of ligands as the most important feature for discriminating highly active compounds from inactives. Molecular dynamics simulation, conducted on highly active ligands, described the dynamic behaviors of the protein complexes and supported the interpretations obtained from other modeling analyses. The current study may be useful for designing PARP-1 inhibitors.

Keywords:

• PARP-1 inhibitors
• KPLS
• MM-GBSA
• pharmacophore mapping
• molecular dynamics

Acknowledgement

Authors are thankful to the All India Council for Technical Education (AICTE), New Delhi, Council for Scientific and Industrial Research (CSIR), New Delhi and University Grants Commission (UGC), New Delhi for providing financial supports. One of the authors (AKH) thanks Council for Scientific and Industrial Research (CSIR), New Delhi for providing a Senior Research Fellowship (SRF). We are also thankful to the authority of Jadavpur University and University of Calcutta for providing us the facilities required for the work.