Molecular docking and 3D-QSAR studies on β-phenylalanine derivatives as dipeptidyl peptidase IV inhibitors

Three-dimensional quantitative structure-activity relationship (3D-QSAR) and molecular docking studies were carried out to explore the binding of 73 inhibitors to dipeptidyl peptidase IV (DPP-IV), and to construct highly predictive 3D-QSAR models using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). The negative logarithm of IC₅₀ (pIC₅₀) was used as the biological activity in the 3D-QSAR study. The CoMFA model was developed by steric and electrostatic field methods, and leave-one-out cross-validated partial least squares analysis yielded a cross-validated value (r_cv² {\hbox{r}}_{{\rm{cv}}}^{\rm{2}} ) of 0.759. Three CoMSIA models developed by different combinations of steric, electrostatic, hydrophobic and hydrogen-bond fields yielded significant r_cv² {\hbox{r}}_{{\rm{cv}}}^{\rm{2}} values of 0.750, 0.708 and 0.694, respectively. The CoMFA and CoMSIA models were validated by a structurally diversified test set of 18 compounds. All of the test compounds were predicted accurately using these models. The mean and standard deviation of prediction errors were within 0.33 and 0.26 for all models. Analysis of CoMFA and CoMSIA contour maps helped identify the structural requirements of inhibitors, with implications for the design of the next generation of DPP-IV inhibitors for the treatment of type 2 diabetes.     

Enaminones 8: CoMFA and CoMSIA studies on some anticonvulsant enaminones

3D-QSAR studies comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were carried out on 26 structurally diverse subcutaneous pentylenetetrazol (scPTZ) active enaminone analogues, previously synthesized in our laboratory. CoMFA and CoMSIA were employed to generate models to define the specific structural and electrostatic features essential for enhanced binding to the putative GABA receptor. The 3D-QSAR models demonstrated a reliable ability to predict the CLog P of the active anticonvulsant enaminones, resulting in a q² of 0.558 for CoMFA, and a q² of 0.698 for CoMSIA. The outcomes of the contour maps for both models provide detailed insight for the structural design of novel enaminone derivatives as potential anticonvulsant agents.     

Probing the structural requirements of A-type Aurora kinase inhibitors using 3D-QSAR and molecular docking analysis

Aurora-A, the most widely studied isoform of Aurora kinase overexpressed aberrantly in a wide variety of tumors, has been implicated in early mitotic entry, degradation of natural tumor suppressor p53 and centrosome maturation and separation; hence, potent inhibitors of Aurora-A may be therapeutically useful drugs in the treatment of various forms of cancer. Here, we report an in silico study on a group of 220 reported Aurora-A inhibitors with six different substructures. Three-dimensional quantitative structure–activity relationship (3D-QSAR) studies were carried out using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) techniques on this series of molecules. The resultant optimum 3D-QSAR models exhibited an r _cv² value of 0.404-0.582 and their predictive ability was validated using an independent test set, ending in r _pred² 0.512-0.985. In addition, docking studies were employed to explore these protein–inhibitor interactions at the molecular level. The results of 3D-QSAR and docking analyses validated each other, and the key structural requirements affecting Aurora-A inhibitory activities, and the influential amino acids involved were identified. To the best of our knowledge, this is the first report on 3D-QSAR modeling of Aurora-A inhibitors, and the results can be used to accurately predict the binding affinity of related analogues and also facilitate the rational design of novel inhibitors with more potent biological activities.     

Computational investigations of gram-negative bacteria phosphopantetheine adenylyltransferase inhibitors using 3D-QSAR,molecular docking and molecular dynamic simulations

Abstract

Phosphopantetheine adenylyltransferase (PPAT) has been recognized as a promising target to develop novel antimicrobial agents, which is a hexameric enzyme that catalyzes the penultimate step in coenzyme A biosynthesis. In this work, molecular modeling study was performed with a series of PPAT inhibitors using molecular docking, three-dimensional qualitative structure-activity relationship (3D-QSAR) and molecular dynamic (MD) simulations to reveal the structural determinants for their bioactivities. Molecular docking study was applied to understand the binding mode of PPAT with its inhibitors. Subsequently, 3D-QSAR model was constructed to find the features required for different substituents on the scaffolds. For the best comparative molecular field analysis (CoMFA) model, the Q² and R² values of which were calculated as 0.702 and 0.989, while they were calculated as 0.767 and 0.983 for the best comparative molecular similarity index analysis model. The statistical data verified the significance and accuracy of our 3D-QSAR models. Furthermore, MD simulations were carried out to evaluate the stability of the receptor–ligand contacts in physiological conditions, and the results were consistent with molecular docking studies and 3D-QSAR contour map analysis. Binding free energy was calculated with molecular mechanics generalized born surface area approach, the result of which coincided well with bioactivities and demonstrated that van der Waals accounted for the largest portion. Overall, our study provided a valuable insight for further research work on the recognition of potent PPAT inhibitors.

Communicated by Ramaswamy H. Sarma     

3D-QSAR studies of checkpoint kinase 1 inhibitors based on molecular docking and CoMFA

Three-dimensional quantitative structure–activity relationship (3D-QSAR) studies were performed on a series of substituted 1,4-dihydroindeno[1,2-c]pyrazoles inhibitors, using molecular docking and comparative molecular field analysis (CoMFA). The docking results from GOLD 3.0.1 provide a reliable conformational alignment scheme for the 3D-QSAR model. Based on the docking conformations and alignments, highly predictive CoMFA model was built with cross-validated q ² value of 0.534 and non-cross-validated partial least-squares analysis with the optimum components of six showed a conventional r ² value of 0.911. The predictive ability of this model was validated by the testing set with a conventional r ² value of 0.812. Based on the docking and CoMFA, we have identified some key features of the 1,4-dihydroindeno[1,2-c]pyrazoles derivatives that are responsible for checkpoint kinase 1 inhibitory activity. The analyses may be used to design more potent 1,4-dihydroindeno[1,2-c]pyrazoles derivatives and predict their activity prior to synthesis.     

Docking-based 3D-QSAR study for 11beta-HSD1 inhibitors

11beta-Hydroxysteroid dehydrogenase (11beta-HSD) enzymes catalyze the conversion of biologically inactive 11-ketosteroids into their active 11beta-hydroxy derivatives and vice versa. 11beta-HSD1 has been studied as a potential treatment for metabolic disease such as diabetes and obesity. To find correlation between 11beta-HSD1 and inhibitors, three-dimensional quantitative structure-activity relationship (3D-QSAR) studies were performed on 70 inhibitors, based on molecular docking conformations obtained by using FlexX-Pharm. The studies include comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). Based on the docking results, highly predictive 3D-QSAR models were developed with q(2) values of 0.543 and 0.519 for CoMFA and CoMSIA, respectively. A comparison of the 3D-QSAR field contributions with the structural features of the binding site showed good correlation between the two analyses. Therefore, these results should be useful to the prediction of the activities of new 11beta-HSD1 inhibitors.     

Molecular docking and 3D-QSAR study on 4-(1H-indazol-4-yl) phenylamino and aminopyrazolopyridine urea derivatives as kinase insert domain receptor (KDR) inhibitors

Vascular endothselial growth factor (VEGF) and its receptor tyrosine kinase VEGFR-2 or kinase insert domain receptor (KDR) have been identified as new promising targets for the design of novel anticancer agents. It is reported that 4-(1H-indazol-4-yl)phenylamino and aminopyrazolopyridine urea derivatives exhibit potent inhibitory activities toward KDR. To investigate how their chemical structures relate to the inhibitory activities and to identify the key structural elements that are required in the rational design of potential drug candidates of this class, molecular docking simulations and three-dimensional quantitative structure-activity relationship (3D-QSAR) methods were performed on 78 4-(1H-indazol-4-yl)phenylamino and aminopyrazolopyridine urea derivatives as KDR inhibitors. Surflex-dock was used to determine the probable binding conformations of all the compounds at the active site of KDR. As a result, multiple hydrophobic and hydrogen-bonding interactions were found to be two predominant factors that may be used to modulate the inhibitory activities. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) 3D-QSAR models were developed based on the docking conformations. The CoMFA model produced statistically significant results with the cross-validated correlation coefficient q² of 0.504 and the non-cross-validated correlation coefficient r² of 0.913. The best CoMSIA model was obtained from the combination of steric, electrostatic and hydrophobic fields. Its q² and r² being 0.595 and 0.947, respectively, indicated that it had higher predictive ability than the CoMFA model. The predictive abilities of the two models were further validated by 14 test compounds, giving the predicted correction coefficients r_pred² of 0.727 for CoMFA and 0.624 for CoMSIA, respectively. In addition, the CoMFA and CoMSIA models were used to guide the design of a series of new inhibitors of this class with predicted excellent activities. Thus, these models may be used as an efficient tool to predict the inhibitory activities and to guide the future rational design of 4-(1H-indazol-4-yl)phenylamino and aminopyrazolopyridine urea derivatives-based novel KDR inhibitors with potent activities.     

Receptor based 3D-QSAR to identify putative binders of Mycobacterium tuberculosis Enoyl acyl carrier protein reductase

In the current study, the applicability and scope of 3D-QSAR models (CoMFA and CoMSIA) to complement virtual screening using 3D pharmacophore and molecular docking is examined and applied to identify potential hits against Mycobacterium tuberculosis Enoyl acyl carrier protein reductase (MtENR). Initially CoMFA and CoMSIA models were developed using series of structurally related arylamides as MtENR inhibitors. Docking studies were employed to position the inhibitors into MtENR active site to derive receptor based 3D-QSAR models. Both CoMFA and CoMSIA yielded significant cross validated q² values of 0.663 and 0.639 and r² values of 0.989 and 0.963, respectively. The statistically significant models were validated by a test set of eight compounds with predictive r² value of 0.882 and 0.875 for CoMFA and CoMSIA. The contour maps from 3D-QSAR models in combination with docked binding structures help to better interpret the structure activity relationship. Integrated with CoMFA and CoMSIA predictive models structure based (3D-pharmacophore and molecular docking) virtual screening have been employed to explore potential hits against MtENR. A representative set of 20 compounds with high predicted IC₅₀ values were sorted out in the present study.     

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.     

Binding conformations and QSAR of CA-4 analogs as tubulin inhibitors

A theoretical study on the binding conformations and the quantitative structure–activity relationship (QSAR) of combretastatin A4 (CA-4) analogs as inhibitors toward tubulin has been carried out using docking analysis and comparative molecular field analysis (CoMFA). The appropriate binding orientations and conformations of these compounds interacting with tubulin were revealed by the docking study; and a 3D-QSAR model showing significant statistical quality and satisfactory predictive ability was established, in which the correlation coefficient (R²) and cross-validation coefficient (q²) were 0.955 and 0.66, respectively. The same model was further applied to predict the pIC₅₀ values for 16 congeneric compounds as external test set, and the predictive correlation coefficient R²_pred reached 0.883. Other tests on additional validations further confirmed the satisfactory predictive power of the model. In this work, it was very interesting to find that the 3D topology structure of the active site of tubulin from the docking analysis was in good agreement with the 3D-QSAR model from CoMFA for this series of compounds. Some key structural factors of the compounds responsible for cytotoxicity were reasonably presented. These theoretical results can offer useful references for understanding the action mechanism and directing the molecular design of this kind of inhibitor with improved activity.     

Design,synthesis and biological evaluation of novel oseltamivir derivatives as potent neuraminidase inhibitors

Neuraminidase (NA) is one of the particular potential targets for novel antiviral therapy. In this work, a series of neuraminidase inhibitors with the cyclohexene scaffold were studied based upon the combination of 3D-QSAR, molecular docking, and molecular dynamics techniques. The results indicate that the built 3D-QSAR models yield reliable statistical information: the correlation coefficient (r²) and cross-validation coefficient (q²) of CoMFA (comparative molecular field analysis) are 0.992 and 0.819; the r² and q² of CoMSIA (comparative molecular similarity analysis) are 0.992 and 0.863, respectively. Molecular docking and MD simulations were conducted to confirm the detailed binding mode of enzyme-inhibitor system. The new NA inhibitors had been designed, synthesized, and their inhibitory activities against group-1 neuraminidase were determined. One agent displayed excellent neuraminidase inhibition, with IC₅₀ value of 39.6?μM against NA, while IC₅₀ value for oseltamivir is 61.1?μM. This compound may be further investigated for the treatment of infection by the new type influenza virus.     

Use of molecular modeling, docking, and 3D-QSAR studies for the determination of the binding mode of benzofuran-3-yl-(indol-3-yl)maleimides as GSK-3β inhibitors

Molecular modeling and docking studies along with three-dimensional quantitative structure relationships (3D-QSAR) studies have been used to determine the correct binding mode of glycogen synthase kinase 3β (GSK-3β) inhibitors. The approaches of comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) are used for the 3D-QSAR of 51 substituted benzofuran-3-yl-(indol-3-yl)maleimides as GSK-3β inhibitors. Two binding modes of the inhibitors to the binding site of GSK-3β are investigated. The binding mode 1 yielded better 3D-QSAR correlations using both CoMFA and CoMSIA methodologies. The three-component CoMFA model from the steric and electrostatic fields for the experimentally determined pIC₅₀ values has the following statistics: R²(cv) = 0.386 nd SE(cv) = 0.854 for the cross-validation, and R² = 0.811 and SE = 0.474 for the fitted correlation. F (3,47) = 67.034, and probability of R² = 0 (3,47) = 0.000. The binding mode suggested by the results of this study is consistent with the preliminary results of X-ray crystal structures of inhibitor-bound GSK-3β. The 3D-QSAR models were used for the estimation of the inhibitory potency of two additional compounds.     

Molecular docking and 3D-QSAR studies of HIV-1 protease inhibitors

HIV-1 protease is an obligatory enzyme in the replication process of the HIV virus. The abundance of structural information on HIV-1PR has made the enzyme an attractive target for computer-aided drug design strategies. The daunting ability of the virus to rapidly generate resistant mutants suggests that there is an ongoing need for new HIV-1PR inhibitors with better efficacy profiles and reduced toxicity. In the present investigation, molecular modeling studies were performed on a series of 54 cyclic urea analogs with symmetric P2/P2′ substituents. The binding modes of these inhibitors were determined by docking. The docking results also provided a reliable conformational superimposition scheme for the 3D-QSAR studies. To gain insight into the steric, electrostatic, hydrophobic and hydrogen-bonding properties of these molecules and their influence on the inhibitory activity, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were performed. Two different alignment schemes viz. receptor-based and atom-fit alignment, were used in this study to build the QSAR models. The derived 3D-QSAR models were found to be robust with statistically significant r ² and r ² _pred values and have led to the identification of regions important for steric, hydrophobic and electronic interactions. The predictive ability of the models was assessed on a set of molecules that were not included in the training set. Superimposition of the 3D-contour maps generated from these models onto the active site of enzyme provided additional insight into the structural requirements of these inhibitors. The CoMFA and CoMSIA models were used to design some new inhibitors with improved binding affinity. Pharmacokinetic and toxicity predictions were also carried out for these molecules to gauge their ADME and safety profile. The computational results may open up new avenues for synthesis of potent HIV-1 protease inhibitors.     

3D-QSAR and molecular docking studies of azaindole derivatives as Aurora B kinase inhibitors

The Aurora kinases have been regarded as attractive targets for the development of new anticancer agents. Recently a series of azaindole derivatives with Aurora B inhibitory activities were reported. To explore the relationship between the structures of substituted azaindole derivatives and their inhibition of Aurora B, 3D-QSAR and molecular docking studies were performed on a dataset of 41 compounds. 3D-QSAR, including CoMFA and CoMSIA, were applied to identify the key structures impacting their inhibitory potencies. The CoMSIA model showed better results than CoMFA, with r ² _cv value of 0.575 and r ² value of 0.987. 3D contour maps generated from CoMFA and CoMSIA along with the docking binding structures provided enough information about the structural requirements for better activity. Based on the structure-activity relationship revealed by the present study, we have designed a set of novel Aurora B inhibitors that showed excellent potencies in the developed models. Thus, our results allowed us to design new derivatives with desired activities.     

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.     

Exploration of a binding mode of indole amide analogues as potent histone deacetylase inhibitors and 3D-QSAR analyses

Docking simulations and three-dimensional quantitative structure-activity relationship (3D-QSAR) analyses were conducted on a series of indole amide analogues as potent histone deacetylase inhibitors. The studies include comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). Selected ligands were docked into the active site of human HDAC1. Based on the docking results, a novel binding mode of indole amide analogues in the human HDAC1 catalytic core is presented, and enzyme/inhibitor interactions are discussed. The indole amide group is located in the open pocket, and anchored to the protein through a pair of hydrogen bonds with Asp99 O-atom and amide NH group on ligand. Based on the binding mode, predictive 3D-QSAR models were established, which had conventional r2 and cross-validated coefficient values (r(cv)2) up to 0.982 and 0.601 for CoMFA and 0.954 and 0.598 for CoMSIA, respectively. A comparison of the 3D-QSAR field contributions with the structural features of the binding site showed good correlation between the two analyses. The results of 3D-QSAR and docking studies validate each other and provided insight into the structural requirements for activity of this class of molecules as HDAC inhibitors. The CoMFA and CoMSIA PLS contour maps and MOLCAD-generated active site electrostatic, lipophilicity, and hydrogen-bonding potential surface maps, as well as the docking studies, provided good insights into inhibitor-HDAC interactions at the molecular level. Based on these results, novel molecules with improved activity can be designed.     

Key Features for Designing Phosphodiesterase-5 Inhibitors

Abstract

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², and non cross-validated coefficient, r², 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).     

Docking-based 3D-QSAR (CoMFA,CoMFA-RG,CoMSIA) study on hydroquinoline and thiazinan-4-one derivatives as selective COX-2 inhibitors

A series of 26 selective COX-2 inhibitors which reported previously by our laboratory was selected to generate three-dimensional quantitative structure activity relationship (3D-QSAR) model. Active conformation of each molecule was predicted by docking studies and used for molecular alignment. Activity of 20 molecules as a train set was predicted using three methods including comparative molecular field analysis (CoMFA), CoMFA region focusing (CoMFA-RG) and comparative molecular similarity index analysis (CoMSIA). The best models of CoMFA-RG and CoMSIA revealed correlation coefficients r² of 0.955 and 0.947, the leave one out cross-validation coefficients q² of 0.573 and 0.574, respectively. In addition, CoMFA-RG and CoMSIA models were validated by a test set of six molecules with predicted coefficients r²_pred of 0.644 and 0.799, respectively. Contour maps of generated models provided fruitful information about structural aspect of molecules that affected their COX-2 inhibitory activity. Based on three models results, steric and electrostatic properties are the most important factors in controlling the activity of the molecules. Results of CoMFA-RG and CoMSIA models were utilized to design new molecules. Comparison of experimental and predicted pIC₅₀ values of designed molecules indicated that CoMFA-RG had the more predictive ability.

Communicated by Ramaswamy H. Sarma     

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.