Combinatorial In Silico Strategy towards Identifying Potential Hotspots during Inhibition of Structurally Identical HDAC1 and HDAC2 Enzymes for Effective Chemotherapy against Neurological Disorders

Histone deacetylases (HDACs) regulate epigenetic gene expression programs by modulating chromatin architecture and are required for neuronal development. Dysregulation of HDACs and aberrant chromatin acetylation homeostasis have been implicated in various diseases ranging from cancer to neurodegenerative disorders. Histone deacetylase inhibitors (HDACi), the small molecules interfering HDACs have shown enhanced acetylation of the genome and are gaining great attention as potent drugs for treating cancer and neurodegeneration. HDAC2 overexpression has implications in decreasing dendrite spine density, synaptic plasticity and in triggering neurodegenerative signalling. Pharmacological intervention against HDAC2 though promising also targets neuroprotective HDAC1 due to high sequence identity (94%) with former in catalytic domain, culminating in debilitating off-target effects and creating hindrance in the defined intervention. This emphasizes the need of designing HDAC2-selective inhibitors to overcome these vicious effects and for escalating the therapeutic efficacy. Here we report a top-down combinatorial in silico approach for identifying the structural variants that are substantial for interactions against HDAC1 and HDAC2 enzymes. We used extra-precision (XP)-molecular docking, Molecular Mechanics Generalized Born Surface Area (MMGBSA) for predicting affinity of inhibitors against the HDAC1 and HDAC2 enzymes. Importantly, we employed a novel in silico strategy of coupling the state-of-the-art molecular dynamics simulation (MDS) to energetically-optimized structure based pharmacophores (e-Pharmacophores) method via MDS trajectory clustering for hypothesizing the e-Pharmacophore models. Further, we performed e-Pharmacophores based virtual screening against phase database containing millions of compounds. We validated the data by performing the molecular docking and MM-GBSA studies for the selected hits among the retrieved ones. Our studies attributed inhibitor potency to the ability of forming multiple interactions and infirm potency to least interactions. Moreover, our studies delineated that a single HDAC inhibitor portrays differential features against HDAC1 and HDAC2 enzymes. The high affinity and selective HDAC2 inhibitors retrieved through e-Pharmacophores based virtual screening will play a critical role in ameliorating neurodegenerative signalling without hampering the neuroprotective isoform (HDAC1).