QSAR and docking studies of novel β-carboline derivatives as anticancer

Quantitative structure–activity relationship (QSAR) studies were performed on β-carboline derivatives for prediction of anticancer activity. The statistically significant 2D-QSAR model having r ² = 0.726 and q ² = 0.654 with pred_r ² = 0.763 was developed by stepwise multiple linear regression method. In order to understand the structural requirement of these β-carboline derivatives, a ligand-based pharmacophore 3D-QSAR model was developed. The five-point pharmacophore hypothesis yielded a 3D-QSAR model with good partial least-square statistics results (r ² = 0.73, Q _ext ²  = 0.755, F = 67.5, SD = 0.245, RMSE = 0.241, Pearson-R = 0.883). A docking study revealed the binding orientations of these derivatives at the active site amino residues of DNA intercalate (PDB ID: 1D12). The results of 2D-QSAR, atom-based 3D-QSAR, and docking studies gave detailed structural insights as well as highlighted important binding features of β-carboline derivatives as anticancer agent which provided guidance for the rational design of novel potent anticancer agents.     

QSAR studies for some thiazolidine-2,4-dione derivatives as PIM-2 kinase inhibitors

Quantitative structure–activity relationship (QSAR) studies were performed on a series of 21 thiazolidine-2,4-dione derivatives to find the structural requirements for PIM-2 kinase inhibitory activity by two-dimensional (2D-QSAR), group-based (G-QSAR) and three-dimensional (3D-QSAR) studies. In the present study, widely used technique viz. stepwise forward–backward (SW-FB) has been applied for the development of 2D- and G-QSAR as variable selection method. The statistically significant best 2D-QSAR model was developed by partial least squares regression (PLSR) having r ² = 0.78, q ² = 0.63 with pred_r ² = 0.78. The statistically significant best G-QSAR model was developed by PLSR method having r ² = 0.89, q ² = 0.79 and pred_r ² = 0.82. The 3D-QSAR studies were performed by k-nearest neighbor molecular field analysis along with genetic algorithm method which showed q ² = 0.64 and pred_r ² = 0.94. A docking study revealed the binding orientations of these inhibitors at active site amino acid residues (PHE 43, ASP 124, ASP 182 and GLU 83) of PIM-2 enzyme (PDB ID: 3IWI). The results of this study may be useful to (medicinal) chemists to design more potent thiazolidine-2,4-dione analogs as PIM-2 kinase inhibitors.     

Molecular modeling studies of some substituted chalcone derivatives as cysteine protease inhibitors

The quantitative structure–activity relationship (QSAR) studies were performed on a series of 42 chalcone derivatives to find out the structural requirements of their antimalarial activities. The multiple linear regression (MLR) and partial least square (PLS) regression methods coupled with various feature selection methods, viz., stepwise (SW), genetic algorithm (GA) and simulated annealing (SA) were applied to derive QSAR models which were further validated for statistical significance and predictive ability by internal and external validation. The statistically significant 2D-QSAR model having r ²?=?0.8892 and q ²?=?0.7508 with pred_r ²?=?0.7403 was developed by SW-MLR and best Group-based QSAR (GQSAR) model having r ²?=?0.7884 and q ²?=?0.7038 with pred_r ²?=?0.7339 was developed by SW-PLS. Molecular field analysis was used to construct the best 3D-QSAR model using k-nearest neighbour method, showing good correlative and predictive capabilities in terms of q ²?=?0.6818 and pred_r ²?=?0.7708. A docking study revealed the binding orientations of these inhibitors at active site amino acid residues (Gln36, Cys39, Lys37, Asp35, Trp206) of falcipain-2 enzyme (PDB ID: 3BPF). Both QSAR and docking studies of such derivatives provide guidance for further lead optimization and designing of more potent antimalarial agents.     

A comprehensive structure–activity analysis 2,3,5-trisubstituted 4,5-dihydro-4-oxo-3H-imidazo [4,5-c] pyridine derivatives as angiotensin II receptor antagonists: using 2D- and 3D-QSAR approach

The QSAR studies were performed on a series of 2,3,5-trisubstituted 4,5-dihydro-4-oxo-3H-imidazo [4,5-c] pyridine derivatives as angiotensin II AT₁ receptor antagonists activity to find out the structural features requirements for the antihypertensive activity. The QSAR study was carried out on V-life Molecular Design Suite software and the derived best QSAR model by partial least square principal component regression and multiple linear regression method showed variation in biological activity. The statistically best significant model with high-correlation coefficient (r ² = 0.9425) was selected for further study and the resulted validation parameters of the model, cross-validated correlation coefficient (q ² = 0.7786 and pred_r ² = 0.8562) show the model has good predictive activities. The model showed that the parameters SdssCcount, SssNHcount, and SaaaCcount and H_Donor count are highly correlated with angiotensin II AT₁ receptor antagonists activity of 2,3,5-trisubstituted 4,5-dihydro-4-oxo-3H-imidazo [4,5-c] pyridine derivatives. Partial least square (PLS) methodology coupled with various feature selection methods viz. stepwise, simulated annealing and genetic algorithm were applied to derive 3D-QSAR models which were further validated for statistical significance and predictive ability by internal and external validation. Molecular field analysis was used to construct the best 3D-QSAR model-7 using k-nearest neighbor (kNN) method, showing good correlative and predictive capabilities in terms of q ² = 0.8316 and pred_r ² = 0.8152. Both 2D-and 3D-QSAR study of such derivatives provide guidance for further lead optimization and designing of potent anti-hypertensive agents.     

Pharmacophore modeling and 3D QSAR analysis of isothiazolidinedione derivatives as PTP1B inhibitors

The article describes the development of a robust pharmacophore model and investigation of structure activity relationship analysis of 56 isothiazolidinedione derivatives reported as PTP1B inhibitors. A six-point pharmacophore model consisting of four aromatic rings (R), one hydrogen bond donor (D) and one hydrogen bond acceptor (A) with discrete geometries as pharmacophoric features was developed and the generated pharmacophore model was used to derive a predictive 3D QSAR model for the studied dataset. The obtained 3D QSAR model has an excellent correlation coefficient value (r ² = 0.98) along with good statistical significance as shown by a high Fisher ratio (F = 428.60). The model also exhibits good predictive power confirmed by the high value of cross-validated correlation coefficient (q ² = 0.62). The QSAR model suggests that hydrophobic aromatic character is crucial for the PTP1B inhibitory activity at the R-15 site.     

Molecular docking and 3D-QSAR studies on checkpoint kinase 1 inhibitors

Checkpoint kinase 1(Chk1) is a promising target for cancer treatment. Here three-dimensional quantitative structure–activity relationship (3D-QSAR) studies were performed on 174 1,4-dihydroindeno[1,2-c]pyrazole inhibitors of Chk1 using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). Two satisfactory ligand-based QSAR models were built (CoMFA model: q ² = 0.541, r ² = 0.880, CoMSIA model: q ² = 0.590, r ² = 0.902). The docking-based studies presented a detailed understanding of interaction between the inhibitors and Chk1. The obtained QSAR models are highly predictable (CoMFA model: q ² = 0.567, r ² = 0.891, CoMSIA model: q ² = 0.596, r ² = 0.917). The models were further validated by an external testing set obtaining $ r_{\text{pred}}^{2} $ r pred 2 values 0.896 and 0.923 for CoMFA and CoMSIA, respectively. So our models might be helpful for further modification of 1,4-dihydroindeno[1,2-c]pyrazole derivatives.     

Design of novel focal adhesion kinase inhibitors using 3D-QSAR and molecular docking

Focal adhesion kinase (FAK) plays a primary role in regulating the activity of many signaling molecules. Increased FAK expression has been implicated in a series of cellular processes, including cell migration and survival. Inhibiting the activity of FAK for cancer therapy is currently under investigation. Hence, FAK and its inhibitors have been the subject of intensive research. To understand the structural factors affecting inhibitory potency, molecular docking and 3D-QSAR modeling were studied in this project. CoMFA and CoMSIA methods were used for deriving 3D-QSAR models, which were trained with 78 compounds and then were evaluated for predictive ability with additional 19 compounds. Two substructure-based 3D-QSAR models, including CoMFA model (r ² = 0.930; q ² = 0.629) and CoMSIA model (r ² = 0.952; q ² = 0.586), had a good quality to predict the biological activities of new compounds. Meanwhile, using the flexible docking strategy, two docking-based 3D-QSAR models (CoMFA with r ² = 0.914; q ² = 0.594; CoMSIA with r ² = 0.914; q ² = 0.524) were also constructed. The structure–activity relationship has been illustrated clearly by the contour maps gained from the 3D-QSAR models in combination with the docked binding structures. All the results indicated that it might be useful in the rational design of novel FAK inhibitors.     

A combined 3D QSAR and pharmacophore-based virtual screening for the identification of potent p38 MAP kinase inhibitors: an in silico approach

p38 kinase plays a vital role in inflammation mediated by tumor necrosis factor-α and interleukin-1β pathways. Inhibition of p38 kinase provides an effective way to treat inflammatory diseases. 3D-QSAR study was performed to obtain reliable comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) models for a series of p38 inhibitors with three different alignment methods (Receptor based, atom by atom matching, and pharmacophore based). Among the different alignment methods, better statistics were obtained with receptor-based alignment (CoMFA: q ² = 0.777, r ² = 0.958; CoMSIA: q ² = 0.782, r ² = 0.927). Superposing CoMFA/CoMSIA contour maps on the p38 active site gave a valuable insight to understand physical factors which are important for binding. In addition, this pharmacophore model was used as a 3D query for virtual screening against NCI database. The hit compounds were further filtered by docking and scoring, and their biological activities were predicted by CoMFA and CoMSIA models.     

Three-dimensional quantitative structure–activity relationship (3D-QSAR) analysis and molecular docking-based combined in silico rational approach to design potent and novel TRPV1 antagonists

The discovery of clinically relevant antagonists of TRPV1 for neuropathy pain therapy has proven to be a challenging task. For better understanding of the molecular interactions of antagonists with TRPV1 receptor, a series of chroman and tetrahydroquinoline ureas were analyzed by k-nearest neighbor molecular field analysis (kNN-MFA) and molecular docking. To elucidate the structural properties required for activity as TRPV1 antagonists, we report here kNN-MFA-based 3D-QSAR model for chroman and tetrahydroquinoline ureas as potent TRPV1 antagonists. Sphere exclusion method was used for dividing the compounds into training (26 compounds) and test (5 compounds) set. Overall model classification accuracy was 81.35 % (q ² = 0.8135, representing internal validation) in training set and 81.44 % (pred_r ² = 0.8144, representing external validation) in test set using stepwise forward as a method of variable selection. The stereo view of molecular rectangular grid field of 3D-QSAR using this approach showed that steric and hydrophobic effects dominantly determine binding affinities. Furthermore, the crystal structure of TRPV1 was obtained from protein data bank (PDB code 2NYJ, resolution 3.20 Å), and docking of 31 TRPV1 antagonists into putative binding sites of the TRPV1 were studies. Molecular docking was employed to explore the binding mode between these compounds and the receptor, as well as help understanding the structure–activity relationship revealed by kNN-MFA. Our QSAR model and molecular docking results corroborate with each other and propose directions for the design of new antagonists with better activity toward TRPV1.     

3D quantitative structure–activity relationship for quinoline, benzimidazole and benzofuran-based analogs as phosphodiesterases IV (PDE-IV) inhibitors

PDE-IV is one of the important targets in the treatment of asthma, COPD, and rheumatoid arthritis. In the search for novel PDE-IV inhibitors a 3D-QSAR study was performed on PDE-IV inhibitors by means of pharmacophore mapping using PHASE, Schrödinger-9. The 3D-QSAR obtained from AAHHRR-1024 hypothesis was found to be statistically significant with r ² = 0.9766 and q ² = 0.8759 with 7 PLS factors. The statistical significance of the model was confirmed by a very low value of RMSE 0.4795. The "Pearson-R" value of 0.9376 suggests a very good predictive ability of the hypothesis generated. The present study demonstrates a robust 3D-QSAR model of PDE-IV inhibitor with the help of AAHHRR-1024 hypothesis, which will help in designing novel inhibitors.