Pharmacophore models generation by catalyst and phase consensus-based virtual screening protocol against PI3Kα inhibitors

Phosphoinositide 3-kinases (PI3Ks) family has emerged as promising targets for novel therapeutic agents against neoplastic diseases. Pharmacophore and 3D-quantitative structure–activity relationship modelling were applied to study the structure–activity relationship of PI3K inhibitors. The best HypoGen pharmacophore hypothesis Hypo1 with a correlation coefficient of 0.961 consists of one hydrogen-bond acceptor, one hydrogen-bond donor and two hydrophobic features, whereas the best phase hypothesis AADRRR.378 with favourable statistics (q² = 0.7368, r² = 0.9863) has two hydrogen-bond acceptors, one hydrogen-bond donor and three ring aromatic features. Multiple methods, such as Fischer validation, molecular docking and mapping of test set molecules, were carried out to validate these pharmacophore models. Furthermore, a comparative molecular similarity indices analysis candidate hypothesis model was generated as a supplement of pharmacophore hypothesis. Detailed protein–ligand binding information obtained by Glide was utilised in compound optimisation and virtual screening. A molecular database of 133 known inhibitors and 6179 decoys was built for a screening test to quantitatively analyse various hypotheses and scoring parameters. Finally, we designed a workflow integrating HypoGen pharmacophore searching, phase pharmacophore searching and molecular docking for screening the database. With an improved criterion of enrichment factor (EF = 17.43) and ROC curve (AUC = 0.946), this workflow would provide us an original method for novel PI3K inhibitors.     

Synthesis and pharmacological evaluation of novel isoquinoline N-sulphonylhydrazones designed as ROCK inhibitors

In this study, we synthesized a new congener series of N-sulphonylhydrazones designed as candidate ROCK inhibitors using the molecular hybridization of the clinically approved drug fasudil (1) and the IKK-β inhibitor LASSBio-1524 (2). Among the synthesized compounds, the N-methylated derivative 11 (LASSBio-2065) showed the best inhibitory profile for both ROCK isoforms, with IC₅₀ values of 3.1 and 3.8?µM for ROCK1 and ROCK2, respectively. Moreover, these compounds were also active in the scratch assay performed in human breast cancer MDA-MB 231 cells and did not display toxicity in MTT and LDH assays. Molecular modelling studies provided insights into the possible binding modes of these N-sulphonylhydrazones, which present a new molecular architecture capable of being optimized and developed as therapeutically useful ROCK inhibitors.     

Pharmacophore generation,atom-based 3D-QSAR,molecular docking and molecular dynamics simulation studies on benzamide analogues as FtsZ inhibitors

FtsZ is an appealing target for the design of antimicrobial agent that can be used to defeat the multidrug-resistant bacterial pathogens. Pharmacophore modelling, molecular docking and molecular dynamics (MD) simulation studies were performed on a series of three-substituted benzamide derivatives. In the present study a five-featured pharmacophore model with one hydrogen bond acceptors, one hydrogen bond donors, one hydrophobic and two aromatic rings was developed using 97 molecules having MIC values ranging from .07 to 957 μM. A statistically significant 3D-QSAR model was obtained using this pharmacophore hypothesis with a good correlation coefficient (R² = .8319), cross validated coefficient (Q² = .6213) and a high Fisher ratio (F = 103.9) with three component PLS factor. A good correlation between experimental and predicted activity of the training (R² = .83) and test set (R² = .67) molecules were displayed by ADHRR.1682 model. The generated model was further validated by enrichment studies using the decoy test and MAE-based criteria to measure the efficiency of the model. The docking studies of all selected inhibitors in the active site of FtsZ protein showed crucial hydrogen bond interactions with Val 207, Asn 263, Leu 209, Gly 205 and Asn-299 residues. The binding free energies of these inhibitors were calculated by the molecular mechanics/generalized born surface area VSGB 2.0 method. Finally, a 15 ns MD simulation was done to confirm the stability of the 4DXD–ligand complex. On a wider scope, the prospect of present work provides insight in designing molecules with better selective FtsZ inhibitory potential.     

Pharmacophore mapping,molecular docking and QSAR studies of structurally diverse compounds as CYP2B6 inhibitors

Pharmacophore mapping, molecular docking and quantitative structure–activity relationship (QSAR) studies were carried out for a structurally diverse set of 48 compounds as CYP2B6 inhibitors. The generated best pharmacophore hypotheses from the three methods of conformer generation (FAST, BEST and conformer algorithm based on energy screening and recursive buildup) indicate the importance of two features, namely, hydrogen bond acceptor [electron-rich centre] and ring aromaticity. The distance between the two centres of the important features for ideal inhibitors varied from 5.82 to 6.03 Å. The chemometric tools used for the QSAR analysis were genetic function approximation (GFA) and genetic partial least squares. The developed QSAR models indicate the importance of an electron-rich centre, size of molecule, impact of branching and ring system and distribution of charges in the molecular surface. The docking study confirms the importance of an electron-rich centre for binding with the iron atom of the cytochrome enzyme. A GFA model with spline option was found to be the best model based on internal validation as well as the r ² _{m (overall)} criterion (Q ² = 0.772, r ² _{m (overall)} = 0.774). According to the external prediction statistics (R ² _pred = 0.876), another GFA-derived model with spline option outperforms the remaining models.     

Deciphering the crucial molecular properties of a series of Benzothiazole Hydrazone inhibitors that targets anti-apoptotic Bcl-xL protein

The Bcl-2 family proteins are the central regulators of apoptosis. Due to its predominant role in cancer progression, the Bcl-2 family proteins act as attractive therapeutic targets. Recently, molecular series of Benzothiazole Hydrazone (BH) inhibitors that exhibits drug-likeness characteristics, which selectively targets Bcl-xL have been reported. In the present study, docking was used to explore the plausible binding mode of the highly active BH inhibitor with Bcl-xL; and Molecular Dynamics (MD) simulation was applied to investigate the stability of predicted conformation over time. Furthermore, the molecular properties of the series of BH inhibitors were extensively investigated by pharmacophore based 3D-QSAR model. The docking correctly predicted the binding mode of the inhibitor inside the Bcl-xL hydrophobic groove, whereas the MD-based free energy calculation exhibited the binding strength of the complex over the time period. Furthermore, the residue decomposition analysis revealed the major energy contributing residues – F105, L108, L130, N136, and R139 – involved in complex stability. Additionally, a six-featured pharmacophore model – AAADHR.89 – was developed using the series of BH inhibitors that exhibited high survival score. The statistically significant 3D-QSAR model exhibited high correlation co-efficient (R² = .9666) and cross validation co-efficient (Q² = .9015) values obtained from PLS regression analysis. The results obtained from the current investigation might provide valuable insights for rational drug design of Bcl-xL inhibitor synthesis.     

Exploration of potential RSK2 inhibitors by pharmacophore modelling,structure-based 3D-QSAR,molecular docking study and molecular dynamics simulation

Abstract

The p90 ribosomal s6 kinase 2 (RSK2) is a promising target because of its over expression and activation in human cancer cells and tissues. Over the last few years, significant efforts have been made in order to develop RSK2 inhibitors to treat myeloma, prostatic cancer, skin cancer and etc., but with limited success so far. In this paper, pharmacophore modelling, molecular docking study and molecular dynamics (MD) simulation have been performed to explore the novel inhibitors of RSK2. Pharmacophore models were developed by 95 molecules having pIC₅₀ ranging from 4.577 to 9.000. The pharmacophore model includes one hydrogen bond acceptor (A), one hydrogen bond donor (D), one hydrophobic feature (H) and one aromatic ring (R). It is the best pharmacophore hypothesis that has the highest correlation coefficient (R² = 0.91) and cross validation coefficient (Q² = 0.71) at 5 component PLS factor. It was evaluated using enrichment analysis and the best model was used for virtual screening. The constraints used in this study were docking score, ADME properties, binding free energy estimates and IFD Score to screen the database. Ultimately, 12 hits were identified as potent and novel RSK2 inhibitors. A 15 ns molecular dynamics (MD) simulation was further employed to validate the reliability of the docking results.     

Discovery and optimization of indoles and 7-azaindoles as Rho kinase (ROCK) inhibitors (part-I)

Rho kinase (ROCK) inhibitors are potential therapeutic agents to treat disorders such as hypertension, multiple sclerosis, cancers, and glaucoma. Here, we disclose the synthesis, optimization, biological evaluation of potent indole and 7-azaindole based ROCK inhibitors that have high potency on ROCK (IC(50)=1 nM) with 740-fold selectivity over PKA (47). Moreover, 47 showed very good DMPK properties making it a good candidate for further development. Finally, docking studies with a homology model of ROCK-II were performed to rationalize the binding mode of these compounds and showed the compounds bound in both orientations to take advantage to H-bonds with Lys-121 of ROCK-II.     

A combined molecular modeling study on a series of pyrazole/isoxazole based human Hsp90α inhibitors

Inhibition of the protein chaperone Hsp90α is a promising approach for cancer therapy. In this work, a molecular modeling study combining pharmacophore model, molecular docking and three-dimensional quantitative structure-activity relationships (3D-QSAR) was performed to investigate a series of pyrazole/isoxazole scaffold inhibitors of human Hsp90α. The pharmacophore model can provide the essential features required for the biological activities of the inhibitors. The molecular docking study can give insight into the binding mode between Hsp90α and its inhibitors. 3D-QSAR based on CoMFA and CoMSIA models were performed from three different strategies for conformational selection and alignment. The receptor-based models gave the most statistically significant results with cross-validated q ² values of 0.782 and 0.829 and r ² values of 0.909 and 0.968, for CoMFA and CoMSIA respectively. Furthermore, the 3D contour maps superimposed within the binding site of Hsp90α could help to understand the pivotal interaction and the structural requirements for potent Hsp90α inhibitors. The results show 4-position of pyrazole/isoxazole ring requires bulky and hydrophobic groups, and bulky and electron repulsion substituent of 5-amides is favorable for enhancing activity. This study will be helpful for the rational design of new potent Hsp90α inhibitors.     

Hierarchical virtual screening of the dual MMP-2/HDAC-6 inhibitors from natural products based on pharmacophore models and molecular docking

The dual-target inhibitors tend to improve the response rate in treating tumors, comparing with the single-target inhibitors. Matrix metalloproteinase-2 (MMP-2) and histone deacetylase-6 (HDAC-6) are attractive targets for cancer therapy. In this study, the hierarchical virtual screening of dual MMP-2/HDAC-6 inhibitors from natural products is investigated. The pharmacophore model of MMP-2 inhibitors is built based on ligands, but the pharmacophore model of HDAC-6 inhibitors is built based on the experimental crystal structures of multiple receptor–ligand complexes. The reliability of these two pharmacophore models is validated subsequently. The hierarchical virtual screening, combining these two different pharmacophore models of MMP-2 and HDAC-6 inhibitors with molecular docking, is carried out to identify the dual MMP-2/HDAC-6 inhibitors from a database of natural products. The four potential dual MMP-2/HDAC-6 inhibitors of natural products, STOCK1 N-46177, STOCK1 N-52245, STOCK1 N-55477, and STOCK1 N-69706, are found. The studies of binding modes show that the screened four natural products can simultaneously well bind with the MMP-2 and HDAC-6 active sites by different kinds of interactions, to inhibit the MMP-2 and HDAC-6 activities. In addition, the ADMET properties of screened four natural products are assessed. These found dual MMP-2/HDAC-6 inhibitors of natural products could serve as the lead compounds for designing the new dual MMP-2/HDAC-6 inhibitors having higher biological activities by carrying out structural modifications and optimizations in the future studies.     

Multiple e-pharmacophore modelling pooled with high-throughput virtual screening,docking and molecular dynamics simulations to discover potential inhibitors of Plasmodium falciparum lactate dehydrogenase (PfLDH)

Development of new antimalarial drugs continues to be of huge importance because of the resistance of malarial parasite towards currently used drugs. Due to the reliance of parasite on glycolysis for energy generation, glycolytic enzymes have played important role as potential targets for the development of new drugs. Plasmodium falciparum lactate dehydrogenase (PfLDH) is a key enzyme for energy generation of malarial parasites and is considered to be a potential antimalarial target. Presently, there are nearly 15 crystal structures bound with inhibitors and substrate that are available in the protein data bank (PDB). In the present work, we attempted to consider multiple crystal structures with bound inhibitors showing affinity in the range of 1.4 × 10²–1.3 × 10⁶ nM efficacy and optimized the pharmacophore based on the energy involved in binding termed as e-pharmacophore mapping. A high throughput virtual screening (HTVS) combined with molecular docking, ADME predictions and molecular dynamics simulation led to the identification of 20 potential compounds which could be further developed as novel inhibitors for PfLDH.     

Insight into the structural requirements of thiophene-3-carbonitriles-based MurF inhibitors by 3D-QSAR,molecular docking and molecular dynamics study

The discovery of clinically relevant inhibitors against MurF enzyme has proven to be a challenging task. In order to get further insight into the structural features required for the MurF inhibitory activity, we performed pharmacophore and atom-based three-dimensional quantitative structure–activity relationship studies for novel thiophene-3-carbonitriles based MurF inhibitors. The five-feature pharmacophore model was generated using 48 inhibitors having IC₅₀ values ranging from 0.18 to 663?μm. The best-fitted model showed a higher coefficient of determination (R²?=?0.978), cross-validation coefficient (Q²?=?0.8835) and Pearson coefficient (0.9406) at four component partial least-squares factor. The model was validated with external data set and enrichment study. The effectiveness of the docking protocol was validated by docking the co-crystallized ligand into the catalytic pocket of MurF enzyme. Further, binding free energy calculated by the molecular mechanics generalized Born surface area approach showed that van der Waals and non-polar solvation energy terms are the main contributors to ligand binding in the active site of MurF enzyme. A 10-ns molecular dynamic simulation was performed to confirm the stability of the 3ZM6-ligand complex. Four new molecules are also designed as potent MurF inhibitors. These results provide insights regarding the development of novel MurF inhibitors with better binding affinity.     

Ligand-based and structure-based approaches in identifying ideal pharmacophore against c-Jun N-terminal kinase-3

Structure and ligand based pharmacophore modeling and docking studies carried out using diversified set of c-Jun N-terminal kinase-3 (JNK3) inhibitors are presented in this paper. Ligand based pharmacophore model (LBPM) was developed for 106 inhibitors of JNK3 using a training set of 21 compounds to reveal structural and chemical features necessary for these molecules to inhibit JNK3. Hypo1 consisted of two hydrogen bond acceptors (HBA), one hydrogen bond donor (HBD), and a hydrophobic (HY) feature with a correlation coefficient (r²) of 0.950. This pharmacophore model was validated using test set containing 85 inhibitors and had a good r² of 0.846. All the molecules were docked using Glide software and interestingly, all the docked conformations showed hydrogen bond interactions with important hinge region amino acids (Gln155 and Met149) and these interactions were compared with Hypo1 features. The results of ligand based pharmacophore model (LBPM) and docking studies are validated each other. The structure based pharmacophore model (SBPM) studies have identified additional features, two hydrogen bond donors and one hydrogen bond acceptor. The combination of these methodologies is useful in designing ideal pharmacophore which provides a powerful tool for the discovery of novel and selective JNK3 inhibitors.     

Novel ROCK inhibitors for the treatment of pulmonary arterial hypertension

A novel class of selective inhibitors of ROCK1 and ROCK2 has been identified by structural based drug design. PK/PD experiments using a set of highly selective Rho kinase inhibitors suggest that systemic Rho kinase inhibition is linked to a reversible reduction in lymphocyte counts. These results led to the consideration of topical delivery of these molecules, and to the identification of a lead molecule 7 which shows promising PK and PD in a murine model of pulmonary hypertension after intra-tracheal dosing.     

Dissecting the roles of ROCK isoforms in stress-induced cell detachment

The homologous Rho kinases, ROCK1 and ROCK2, are involved in stress fiber assembly and cell adhesion and are assumed to be functionally redundant. Using mouse embryonic fibroblasts (MEFs) derived from ROCK1^−/− and ROCK2^−/− mice, we have recently reported that they play different roles in regulating doxorubicin-induced stress fiber disassembly and cell detachment: ROCK1 is involved in destabilizing the actin cytoskeleton and cell detachment, whereas ROCK2 is required for stabilizing the actin cytoskeleton and cell adhesion. Here, we present additional insights into the roles of ROCK1 and ROCK2 in regulating stress-induced impairment of cell-matrix and cell-cell adhesion. In response to doxorubicin, ROCK1^−/− MEFs showed significant preservation of both focal adhesions and adherens junctions, while ROCK2^−/− MEFs exhibited impaired focal adhesions but preserved adherens junctions compared with the wild-type MEFs. Additionally, inhibition of focal adhesion or adherens junction formations by chemical inhibitors abolished the anti-detachment effects of ROCK1 deletion. Finally, ROCK1^−/− MEFs, but not ROCK2^−/− MEFs, also exhibited preserved central stress fibers and reduced cell detachment in response to serum starvation. These results add new insights into a novel mechanism underlying the anti-detachment effects of ROCK1 deletion mediated by reduced peripheral actomyosin contraction and increased actin stabilization to promote cell-cell and cell-matrix adhesion. Our studies further support the differential roles of ROCK isoforms in regulating stress-induced loss of central stress fibers and focal adhesions as well as cell detachment.     

In silico studies on potential MCF-7 inhibitors: a combination of pharmacophore and 3D-QSAR modeling,virtual screening,molecular docking,and pharmacokinetic analysis

Gallic acid and its derivatives exhibit a diverse range of biological applications, including anti-cancer activity. In this work, a data-set of forty-six molecules containing the galloyl moiety, and known to show anticarcinogenic activity against the MCF-7 human cancer cell line, have been chosen for pharmacophore modeling and 3D-Quantitative Structure Activity Relationship (3D-QSAR) studies. A tree-based partitioning algorithm has been used to find common pharmacophore hypotheses. The QSAR model was generated for three, four, and five featured hypotheses with increasing PLS factors and analyzed. Results for five featured hypotheses with three acceptors and two aromatic rings were the best out of all the possible combinations. On analyzing the results, the most robust (R^2?=?.8990) hypothesis with a good predictive power (Q^2?=?.7049) was found to be AAARR.35. A good external validation (R² = .6109) was also obtained. In order to design new MCF-7 inhibitors, the QSAR model was further utilized in pharmacophore-based virtual screening of a large database. The predicted IC₅₀ values of the identified potential MCF-7 inhibitors were found to lie in the micromolar range. Molecular docking into the colchicine domain of tubulin was performed in order to examine one of the probable mechanisms. This revealed various interactions between the ligand and the active site protein residues. The present study is expected to provide an effective guide for methodical development of potent MCF-7 inhibitors.     

Identification of InhA inhibitors: A combination of virtual screening,molecular dynamics simulations and quantum chemical studies

In the present work, multiple pharmacophore-based virtual screening of the SPECS natural product database was carried out to identify novel inhibitors of the validated biological target, InhA. The pharmacophore models were built from the five different groups of the co-crystallized ligands present within the active site. The generated models with the same features from each group were pooled and subjected to the test set validation, receiver–operator characteristic analysis and Güner–Henry studies. A set of five hypotheses with sensitivity > 0.5, specificity > 0.5, area under curve (AUC) > 0.7, and goodness of hit score > 0.7 were retrieved and exploited for the virtual screening. The common hits (87 molecules) obtained from these hypotheses were processed via drug-likeness filters. The filtered molecules (27 molecules) were compared for the binding modes and the top scored molecules (12 molecules) along with the reference (triclosan (TCL), docking score = ?11.65 kcal/mol) were rescored and reprioritized via molecular mechanics-generalized Born surface area approach. Eventually, the stability of reprioritized (10 molecules) docked complexes was scrutinized via molecular dynamics simulations. Moreover, the quantum chemical studies of the dynamically stable compounds (9 molecules) were performed to understand structural features essential for the activity. Overall, the protocol resulted in the recognition of nine lead compounds that can be targeted against InhA.     

Understanding nucleosomal histone and DNA interactions: a biophysical study

Histones are associated with DNA to form nucleosome essential for chromatin structure and major nuclear processes like gene regulation and expression. Histones consist of H1, H2A, H2B and H3, H4 type proteins. In the present study, combined histones from calf thymus were complexed with ct DNA and their binding affinities were measured fluorimetrically. All the five histones were resolved on SDS page and their binding with DNA was visualized. The values of biding affinities varied with pH and salt concentration. Highest affinity (4.0?×?10⁵ M^?1) was recorded at pH 6.5 in 50 mM phosphate buffer and 1.5?×?10⁴ M^?1 in 2 M NaCl at pH 7.0. The CD spectra support the highest binding affinity with maximum conformational changes at pH 7.0. The time-resolved fluorescence data recorded two life times for histone tyrosine residues at 300 nm emission in phosphate buffer pH 6.5. These life times did not show much change upon binding with DNA in buffer as well as in 2 M NaCl. The isothermal calorimetric studies yielded thermodynamic parameters ΔG, ΔH and ΔS as ?1.6?×?10⁵ cal/mol, ?1.13?×?10³ cal/mol and ?3.80 cal/mol/deg, respectively, evidencing a spontaneous exothermic reaction. The dominant binding forces in building the nucleosome are electrostatic interactions.