Virtual screening based on pharmacophore model followed by docking simulation studies in search of potential inhibitors for p38 map kinase

Discovery of novel isoform-selective histone deacetylases 5 and 9 inhibitors through combined ligand-based pharmacophore modeling, molecular mocking, and molecular dynamics simulations for cancer treatment

Class IIa histone deacetylases (HDACs) 5 and 9 play crucial roles in several human disorders such as cancer, making them important targets for drug design. Continuous research is pursed to overcome the cytotoxicity side effect that comes with the currently available broad-spectrum HDACs inhibitors. Herein, common features of active HDACs inhibitors in clinical trials and use have been calculated to generate the best pharmacophore hypothesis. Guner-Henry scoring system was used to validate the generated hypotheses. Hypo1 of HDAC5 and Hypo2 of HDAC9 exhibited the most statistically significance hypotheses. Compounds with fit value of 3 and more were examined by QuickVina 2 docking tool to calculate their binding affinity toward all class IIa HDACs. A total of 6 potential selective compounds were subjected to 100 molecular dynamics (MD) simulation to examine their binding modes. The free binding energy calculations were computed according to the MM-PBSA method. Proposed selective compounds displayed good stability with their targets and thus they may offer potent leads for the designing of HDAC5 and HDAC9 isoform selective inhibitors.

Natural compounds as potential Hsp90 inhibitors for breast cancer-Pharmacophore guided molecular modelling studies

Breast cancer is one of the major impediments affecting women globally. The ATP-dependant heat shock protein 90 (Hsp90) forms the central component of molecular chaperone machinery that predominantly governs the folding of newly synthesized peptides and their conformational maturation. It regulates the stability and function of numerous client proteins that are frequently upregulated and/or mutated in cancer cells, therefore, making Hsp90 inhibition a promising therapeutic strategy for the development of new efficacious drugs to treat breast cancer. In the present in silico investigation, a structure-based pharmacophore model was generated with hydrogen bond donor, hydrogen bond acceptor and hydrophobic features complementary to crucial residues Ala55, Lys58, Asp93, Ile96, Met98 and Thr184 directed at inhibiting the ATP-binding activity of Hsp90. Subsequently, the phytochemical dataset of 3210 natural compounds was screened to retrieve the prospective inhibitors after rigorous validation of the model pharmacophore. The retrieved 135 phytocompounds were further filtered by drug-likeness parameters including Lipinski's rule of five and ADMET properties, then investigated via molecular docking-based scoring. Molecular interactions were assessed using Genetic Optimisation for Ligand Docking program for 95 drug-like natural compounds against Hsp90 along with two clinical drugs as reference compounds - Geldanamycin and Radicicol. Docking studies revealed three phytochemicals are better than the investigated clinical drugs. The reference and hit compounds with dock scores of 48.27 (Geldanamycin), 40.90 (Radicicol), 73.04 (Hit1), 72.92 (Hit2) and 68.12 (Hit3) were further validated for their binding stability through molecular dynamics simulations. We propose that the non-macrocyclic scaffolds of three identified phytochemicals might aid in the development of novel therapeutic candidates against Hsp90-driven cancers.

Dibenzepinones, dibenzoxepines and benzosuberones based p38α MAP kinase inhibitors: Their pharmacophore modelling, 3D-QSAR and docking studies

In the present study, a series of dibenzepinones, dibenzoxepines, and benzosuberones targeting p38α MAP kinase were subjected to pharmacophore modelling, 3D-QSAR and molecular docking studies. The IC₅₀ values for these 67 compounds ranged between 0.003 and 6.80 μM. A five-point model (DDHHR.8) was generated using these compounds. This model was found to be statistically significant and was found to have high correlation (R² = 0.98), cross-validation coefficient (Q² = 0.95) and F (330) values at six component PLS factor. Tests were performed to ascertain the efficacy of the generated model. These tests included external validation, Tropsha's test for predictive ability, Y-randomisation test and domain of applicability (APD). In order to check the restrictivity of the model, enrichment studies were performed with inactive compounds by using decoy set molecules. To evaluate the effectiveness of the docking protocol, the co-crystallised ligand was extracted from the ligand-binding domain of the protein and was re-docked into the same position. Both the conformers were then superimposed, suggesting satisfactory docking parameters with an RMSD value of less than 1.0 Å (0.853 Å). A 10 ns molecular dynamics simulation confirmed the docking results of the 3UVP-ligand complex and the presumed active conformation. The outcome of the present study provides insight into the molecular features that promote bioactivity and can be exploited for the prediction of novel potent p38α MAP kinase inhibitors before carrying out their synthesis and anticancer evaluation.

Combiphore (Structure and Ligand Based Pharmacophore) - Approach for the Design of GPR40 Modulators in the Management of Diabetes

BACKGROUND

Pharmacophore mapping and molecular docking can be synergistically integrated to improve the drug design and discovery process. A rational strategy, combiphore approach, derived from the combined study of Structure and Ligand based pharmacophore has been described to identify novel GPR40 modulators.

METHODS

DISCOtech module from Discovery studio was used for the generation of the Structure and Ligand based pharmacophore models which gave hydrophobic aromatic, ring aromatic and negative ionizable as essential pharmacophoric features. The generated models were validated by screening active and inactive datasets, GH scoring and ROC curve analysis. The best model was exposed as a 3D query to screen the hits from databases like GLASS (GPCR-Ligand Association), GPCR SARfari and Mini-Maybridge. Various filters were applied to retrieve the hit molecules having good drug-like properties. A known protein structure of hGPR40 (pdb: 4PHU) having TAK-875 as ligand complex was used to perform the molecular docking studies; using SYBYL-X 1.2 software.

RESULTS AND CONCLUSION

Clustering both the models gave RMSD of 0.89. Therefore, the present approach explored the maximum features by combining both ligand and structure based pharmacophore models. A common structural motif as identified in combiphore for GPR40 modulation consists of the para-substituted phenyl propionic acid scaffold. Therefore, the combiphore approach, whereby maximum structural information (from both ligand and biological protein) is explored, gives maximum insights into the plausible protein-ligand interactions and provides potential lead candidates as exemplified in this study.

Discovery of FIXa inhibitors by combination of pharmacophore modeling, molecular docking, and 3D-QSAR modeling

Human Coagulation Factor IXa (FIXa), specifically inhibited at the initiation stage of the blood coagulation cascade, is an excellent target for developing selective and safe anticoagulants. To explore this inhibitory mechanism, 86 FIXa inhibitors were selected to generate pharmacophore models and subsequently SAR models. Both best pharmacophore model and ROC curve were built through the Receptor-Ligand Pharmacophore Generation module. CoMFA model based on molecular docking and PLS factor analysis methods were developed. Model propagations values are q² = 0.709, r² = 0.949, and r²_pred = 0.905. The satisfactory q² value of 0.609, r² value of 0.962, and r²_pred value of 0.819 for CoMSIA indicated that the CoMFA and CoMSIA models are both available to predict the inhibitory activity on FIXa. On the basis of pharmacophore modeling, molecular docking, and 3D-QSAR modeling screening, six molecules are screened as potential FIXa inhibitors.

New compounds identified through in silico approaches reduce the α-synuclein expression by inhibiting prolyl oligopeptidase in vitro

Prolyl oligopeptidase (POP) is a serine protease that is responsible for the maturation and degradation of short neuropeptides and peptide hormones. The inhibition of POP has been demonstrated in the treatment of α-synucleinopathies and several neurological conditions. Therefore, ligand-based and structure-based pharmacophore models were generated and validated in order to identify potent POP inhibitors. Pharmacophore-based and docking-based virtual screening of a drug-like database resulted in 20 compounds. The in vitro POP assays indicated that the top scoring compounds obtained from virtual screening, Hit 1 and Hit 2 inhibit POP activity at a wide range of concentrations from 0.1 to 10 µM. Moreover, treatment of the hit compounds significantly reduced the α-synuclein expression in SH-SY5Y human neuroblastoma cells, that is implicated in Parkinson’s disease. Binding modes of Hit 1 and Hit 2 compounds were explored through molecular dynamics simulations. A detailed investigation of the binding interactions revealed that the hit compounds exhibited hydrogen bond interactions with important active site residues and greater electrostatic and hydrophobic interactions compared to those of the reference inhibitors. Finally, our findings indicated the potential of the identified compounds for the treatment of synucleinopathies and CNS related disorders.