3D-QSAR and molecular docking studies of 2-pyrimidinecarbonitrile derivatives as inhibitors against falcipain-3

The three dimensional-quantitative structure activity relationship (3D-QSAR) studies were performed on a series of falcipain-3 inhibitors using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) techniques. A training set containing 42 molecules served to establish the QSAR models. The optimum CoMFA and CoMSIA models obtained for the training set were statistically significant with cross-validated correlation coefficients r(cv)(2) (q(2)) of 0.549 and 0.608, and conventional correlation coefficients (r(2)) of 0.976 and 0.932, respectively. An independent test set of 12 molecules validated the external predictive power of both models with predicted correlation coefficients (r(pred)(2)) for CoMFA and CoMSIA as 0.697 and 0.509, respectively. The docking of inhibitors into falcipain-3 active site using GOLD software revealed the vital interactions and binding conformation of the inhibitors. The CoMFA and CoMSIA field contour maps agree well with the structural characteristics of the binding pocket of falcipain-3 active site, which suggests that the information rendered by 3D-QSAR models and the docking interactions can provide guidelines for the development of improved falcipain-3 inhibitors.     

Pharmacophore modeling, virtual screening and 3D-QSAR studies of 5-tetrahydroquinolinylidine aminoguanidine derivatives as sodium hydrogen exchanger inhibitors

Sodium hydrogen exchanger (SHE) inhibitor is one of the most important targets in treatment of myocardial ischemia. In the course of our research into new types of non-acylguanidine, SHE inhibitory activities of 5-tetrahydroquinolinylidine aminoguanidine derivatives were used to build pharmacophore and 3D-QSAR models. Genetic Algorithm Similarity Program (GASP) was used to derive a 3D pharmacophore model which was used in effective alignment of data set. Eight molecules were selected on the basis of structure diversity to build 10 different pharmacophore models. Model 1 was considered as the best model as it has highest fitness score compared to other nine models. The obtained model contained two acceptor sites, two donor atoms and one hydrophobic region. Pharmacophore modeling was followed by substructure searching and virtual screening. The best CoMFA model, representing steric and electrostatic fields, obtained for 30 training set molecules was statistically significant with cross-validated coefficient (q(2)) of 0.673 and conventional coefficient (r(2)) of 0.988. In addition to steric and electrostatic fields observed in CoMFA, CoMSIA also represents hydrophobic, hydrogen bond donor and hydrogen bond acceptor fields. CoMSIA model was also significant with cross-validated coefficient (q(2)) and conventional coefficient (r(2)) of 0.636 and 0.986, respectively. Both models were validated by an external test set of eight compounds and gave satisfactory prediction (r(pred)(2)) of 0.772 and 0.701 for CoMFA and CoMSIA models, respectively. This pharmacophore based 3D-QSAR approach provides significant insights that can be used to design novel, potent and selective SHE inhibitors.     

Three-dimensional quantitative structure-activity relationship (3 D-QSAR) and docking studies on (benzothiazole-2-yl) acetonitrile derivatives as c-Jun N-terminal kinase-3 (JNK3) inhibitors

Three-dimensional quantitative structure-activity relationship (3D-QSAR) models were developed for 44 (benzothiazole-2-yl) acetonitrile derivatives, inhibiting c-Jun N-terminal kinase-3 (JNK3). It includes molecular field analysis (MFA) and receptor surface analysis (RSA). The QSAR model was developed using 34 compounds and its predictive ability was assessed using a test set of 10 compounds. The predictive 3D-QSAR models have conventional r2 values of 0.849 and 0.766 for MFA and RSA, respectively; while the cross-validated coefficient r(cv)2 values of 0.616 and 0.605 for MFA and RSA, respectively. The results of the QSAR model were further compared with a structure-based analysis using docking studies with crystal structure of JNK3. Ligands bind in the ATP pocket and the hydrogen bond with GLN155 was found to be crucial for selectivity among other kinases. The results of 3D-QSAR and docking studies validate each other and hence, the combination of both methodologies provides a powerful tool directed to the design of novel and selective JNK3 inhibitors.     

Three-dimensional quantitative structure-activity relationship (3D-QSAR) studies of tricyclic oxazolidinones as antibacterial agents

Oxazolidinones exemplified by eprezolid and linezolid are a new class of antibacterials that are active against Gram positive and anaerobic bacteria including methicillin-resistant Staphylococcus aureus (MRSA), methicillin-resistant Staphylococcus epidermidis (MRSE) and vancomycin resistant enterococci (VRE). In an effort to have a better antibacterial agent in the oxazolidinone class, we have performed three-dimensional quantitative structure-activity relationship (3D-QSAR) studies for a series of tricyclic oxazolidinones. 3D-QSAR studies were performed using the Comparative Molecular Field Analysis (CoMFA) and Comparative Molecular Similarity Indices Analysis (CoMSIA) procedures. These studies were performed using 42 compounds; the QSAR model was developed using a training set of 33 compounds. The predictive ability of the QSAR model was assessed using a test set of 9 compounds. The predictive 3D-QSAR models have conventional r(2) values of 0.975 and 0.940 for CoMFA and CoMSIA respectively; similarly, cross-validated coefficient q(2) values of 0.523 and 0.557 for CoMFA and CoMSIA, respectively, were obtained. The CoMFA 3D-QSAR model performed better than the CoMSIA model.     

Molecular docking and 3D-QSAR studies of Yersinia protein tyrosine phosphatase YopH inhibitors

Three-dimensional quantitative structure-activity relationship (QSAR) studies were conducted on two classes of recently explored compounds with known YopH inhibitory activities. Docking studies were employed to position the inhibitors into the YopH active site to determine the probable binding conformation. Good correlations between the predicated binding free energies and the inhibitory activities were found for two subsets of phosphate mimetics: alpha-ketocarboxylic acid and squaric acid (R2=0.70 and 0.68, respectively). The docking results also provided a reliable conformational alignment scheme for 3D-QSAR modeling. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were performed based on the docking conformations, giving q2 of 0.734 and 0.754 for CoMFA and CoMSIA models, respectively. The 3D-QSAR models were significantly improved after removal of an outlier (q2=0.829 for CoMFA and q2=0.837 for CoMSIA). The predictive ability of the models was validated using a set of compounds that were not included in the training set. Mapping the 3D-QSAR models to the active site of YopH provides new insight into the protein-inhibitor interactions for this enzyme. These results should be applicable to the prediction of the activities of new YopH inhibitors, as well as providing structural implications for designing potent and selective YopH inhibitors as antiplague agents.     

Molecular docking and 3D-QSAR studies on inhibitors of DNA damage signaling enzyme human PARP-1

Poly (ADP-ribose) polymerase-1 (PARP-1) operates in a DNA damage signaling network. Molecular docking and three dimensional-quantitative structure activity relationship (3D-QSAR) studies were performed on human PARP-1 inhibitors. Docked conformation obtained for each molecule was used as such for 3D-QSAR analysis. Molecules were divided into a training set and a test set randomly in four different ways, partial least square analysis was performed to obtain QSAR models using the comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). Derived models showed good statistical reliability that is evident from their r2, q2(loo) and r2(pred) values. To obtain a consensus for predictive ability from all the models, average regression coefficient r2(avg) was calculated. CoMFA and CoMSIA models showed a value of 0.930 and 0.936, respectively. Information obtained from the best 3D-QSAR model was applied for optimization of lead molecule and design of novel potential inhibitors.     

Understanding the structure-activity and structure-selectivity correlation of cyclic guanine derivatives as phosphodiesterase-5 inhibitors by molecular docking, CoMFA and CoMSIA analyses

Molecular docking and 3D-QSAR analyses were performed to understand how PDE5 and PDE6 interact with a series of (49) cyclic guanine derivatives. Using the conformations of the compounds revealed by molecular docking, CoMFA and CoMSIA analyses resulted in the first quantitative structure-activity relationship (QSAR) and first quantitative structure-selectivity relationship (QSSR) models (with high cross-validated correlation coefficient q(2) and conventional correlation coefficient r(2) values) for predicting the inhibitory activity against PDE5 and the selectivity against PDE6. The high q(2) and r(2) values, along with further testing, indicate that the obtained 3D-QSAR and 3D-QSSR models will be valuable in predicting both the inhibitory activity and selectivity of cyclic guanine derivatives for these protein targets. A set of 3D contour plots drawn based on the 3D-QSAR and 3D-QSSR models reveal some useful clues to improve both the activity and selectivity by modifying structures of the compounds. It has been demonstrated that both the steric and electrostatic factors should appropriately be taken into account in future rational design and development of more active and more selective PDE5 inhibitors for the therapeutic treatment of erectile dysfunction (ED).     

New potent 5-substituted benzofuroxans as inhibitors of Trypanosoma cruzi growth: quantitative structure-activity relationship studies

Benzofuroxan derivatives have been shown to inhibit the growth of Trypanosoma cruzi, the etiological agent of Chagas' disease. Therefore, 2D- and 3D-QSAR models of their in vitro antichagasic activity were developed. Six new derivatives were synthesized to complete a final set of 26 structurally diverse benzofuroxans. The 2D-QSAR model (r = 0.939, r(adj)(2) = 0.849) was generated using multiple regression analysis of tabulated substituents' physicochemical properties and indicator variables. In addition, a 3D-QSAR model (r(2) = 0.997, q(2) = 0.802) was obtained using a comparative molecular field analysis (CoMFA). Due to the well-known benzofuroxan tautomerism, in both approaches (2D- and 3D-QSAR) it was necessary to include an indicator variable to consider the N-oxide position (I(6)). This parameter was established using low-temperature NMR experiments. Both QSAR models identified the electrophilic character of the substituent alpha-atom as a requirement for activity. Further support was found using a density functional theory (DFT) approach.     

Computational studies of COX-2 inhibitors: 3D-QSAR and docking

The 3D-QSAR (three-dimensional quantitative structure-activity relationships) studies for 88 selective COX-2 (cyclooxygenase-2) inhibitors belonging to three chemical classes (triaryl rings, diaryl cycloalkanopyrazoles, and diphenyl hydrazides) were conducted using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). Partial least squares analysis produced statistically significant models with q(2) values of 0.84 and 0.79 for CoMFA and CoMSIA, respectively. The binding energies calculated from flexible docking were correlated with inhibitory activities by the least-squares fit method. The three chemical classes of inhibitors showed reasonable internal predictability (r(2)=0.51, 0.49, and 0.54), but the sulfonyl-containing inhibitors demonstrated distinctively low binding energy compared to the others. The electrostatic interaction energy between the Arg513 of the COX-2 active site and sulfonyl group of the triaryl rings seemed to have the responsibility for difference in binding energy. Comparative binding energy (COMBINE) analyses gave q(2) values of 0.64, 0.63, and 0.50 for triaryl rings, diaryl cycloalkanopyrazoles, and diphenyl hydrazides, respectively. In this COMBINE model, some protein residues were highlighted as particularly important for inhibitory activity. The combination of ligand-based and structure-based models provided an improved understanding in the interaction between the three chemical classes and the COX-2.     

3D-QSAR studies on fluroquinolones derivatives as inhibitors for tuberculosis

A quantitative structure activity relationship (QSAR) study was performed on the fluroquinolones known to have anti-tuberculosis activity. The 3D-QSAR models were generated using stepwise variable selection of the four methods - multiple regression (MR), partial least square regression (PLSR), principal component regression (PCR) and artificial neural networks (kNN-MFA). The statistical result showed a significant correlation coefficient q(2) (90%) for MR model and an external test set of (pred_r(2)) -1.7535, though the external predictivity showed to improve using kNN-MFA method with pred_r(2) of -0.4644. Contour maps showed that steric effects dominantly determine the binding affinities. The QSAR models may lead to a better understanding of the structural requirements of anti-tuberculosis compounds and also help in the design of novel molecules.     

3D-QSAR studies on tripeptide aldehyde inhibitors of proteasome using CoMFA and CoMSIA methods

The ubiquitin-proteasome pathway plays a crucial role in the regulation of many physiological processes and in the development of a number of major human diseases, such as cancer, Alzheimer's, Parkinson's, diabetes, etc. As a new target, the study on the proteasome inhibitors has received much attention recently. Three-dimensional quantitative structure-activity relationship (3D-QSAR) studies using comparative molecule field analysis (CoMFA) and comparative molecule similarity indices analysis (CoMSIA) techniques were applied to analyze the binding affinity of a set of tripeptide aldehyde inhibitors of 20S proteasome. The optimal CoMFA and CoMSIA models obtained for the training set were all statistically significant with cross-validated coefficients (q(2)) of 0.615, 0.591 and conventional coefficients (r(2)) of 0.901, 0.894, respectively. These models were validated by a test set of eight molecules that were not included in the training set. The predicted correlation coefficients (r(2)) of CoMFA and CoMSIA are 0.944 and 0.861, respectively. The CoMFA and CoMSIA field contour maps agree well with the structural characteristics of the binding pocket of beta5 subunit of 20S proteasome, which suggests that the 3D-QSAR models built in this paper can be used to guide the development of novel inhibitors of 20S proteasome.     

Pharmacophore generation and atom-based 3D-QSAR of novel quinoline-3-carbonitrile derivatives as Tpl2 kinase inhibitors

Tumour progression locus-2 (Tpl2) is a serine/threonine kinase, which regulates the expression of tumour necrosis factor α. The article describes the development of a robust pharmacophore model and the investigation of structure-activity relationship analysis of quinoline-3-carbonitrile derivatives reported for Tpl2 kinase inhibition. A five point pharmacophore model (ADRRR) was developed and used to derive a predictive atom-based 3-dimensional quantitative structure activity relationship (3D-QSAR) model. The obtained 3D-QSAR model has an excellent correlation coefficient value (r(2)= 0.96), Fisher ratio (F = 131.9) and exhibited good predictive power (q(2) = 0.79). The QSAR model suggests that the inclusion of hydrophobic substituents will enhance the Tpl2 kinase inhibition. In addition, H-bond donating groups, negative ionic groups and electron withdrawing groups positively contribute to the Tpl2 kinase inhibition. Further, pharmacophoric model was validated by the receiver operating characteristic curve analysis and was employed for virtual screening to identify six potential Tpl2 kinase inhibitors. The findings of this study provide a set of guidelines for designing compounds with better Tpl2 kinase inhibitory potency.     

Design of EGFR kinase inhibitors: a ligand-based approach and its confirmation with structure-based studies

Three-dimensional quantitative structure-activity relationship (3D-QSAR) models were developed for 100 anilinoquinazolines, inhibiting epidermal growth factor receptor (EGFR) kinase. The studies included molecular field analysis (MFA) and receptor surface analysis (RSA). The cross-validated r2 (r2cv) values are 0.81 and 0.79 for MFA and RSA, respectively. The predictive ability of these models was validated by 28 test set molecules. The results of the best QSAR model were further compared with structure-based investigations using docking studies with the crystal structure of EGFR kinase domain. The results helped to understand the nature of substituents at the 6- and 7-positions, thereby providing new guidelines for the design of novel inhibitors.     

Key features for designing phosphodiesterase-5 inhibitors

Phosphodiesterase superfamily is the key regulator of 3',5'-cyclic guanosine monophosphate (cGMP) decomposition in human body. Phosphodiesterase-5 (PDE-5) inhibitors, sildenafil, vardenafil and tadalafil, are well known oral treatment for males with erectile dysfunction. To investigate the inhibitory effects of traditional Chinese medicine (TCM) compounds to PDE-5, we performed both ligand-based and structure-based studies on this topic. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) studies were conducted to construct three dimensional quantitative structure-activity relationship (3D-QSAR) models of series of known PDE-5 inhibitors. The predictive models had cross-validated, q(2), and non cross-validated coefficient, r(2), values of 0.791 and 0.948 for CoMFA and 0.724 and 0.908 for CoMSIA. These two 3D-QSAR models were used to predict activity of TCM compounds. Docking simulations were performed to further analyze the binding mode of training set and TCM compounds. A putative binding model was proposed based on CoMFA and CoMSIA contour maps and docking simulations; formation of pi-stacking, water bridge and specific hydrogen bonding were deemed important interactions between ligands and PDE-5. Of our TCM compounds, engeletin, satisfied our binding model, and hence, emerged as PDE-5 inhibitor candidate. Using this study as an example, we demonstrated that docking should be conducted for qualitative purposes, such as identifying protein characteristics, rather than for quantitative analyses that rank compound efficacy based on results of scoring functions. Prediction of compound activity should be reserved for QSAR analyses, and scoring functions and docking scores should be used for preliminary screening of TCM database (http://tcm.cmu.edu.tw/index.php).     

Quantitative structure-activity analysis of 5-arylidene-2,4-thiazolidinediones as aldose reductase inhibitors

Design of aldose reductase (ALR2) inhibitors has received considerable attention. Aldose reductase inhibitors, when administered from the onset of hyperglycemia, prevent the progression of polyol accumulation-linked complications. The feasibility that inhibition of aldose reductase provides a pharmacologically direct treatment for diabetic complications that is independent of the control of blood sugar levels has spurred the development of structurally diverse aldose reductase inhibitors. In this work, we report quantitative structure-activity relationship (QSAR) analysis performed by 3D-QSAR analysis, Hansch analysis, and Fujita-Ban analysis on a series of 5-arylidene-2,4-thiazolidinediones as aldose reductase inhibitors. The 2D & 3D-QSAR models were generated using 18 compounds and Fujita-Ban analysis models were obtained using 23 compounds. The predictive ability of the resulting 2D and 3D models was evaluated against a test set of 5 compounds. Analyses of results from the present QSAR study inferred that 3rd position of the phenyl ring and acetic acid substitution at N-position of thiazolidinediones play a key role in the aldose reductase inhibitory activity.