CoMFA and CoMSIA studies on C-aryl glucoside SGLT2 inhibitors as potential anti-diabetic agents

SGLT2 has become a target of therapeutic interest in diabetes research. CoMFA and CoMSIA studies were performed on C-aryl glucoside SGLT2 inhibitors (180 analogues) as potential anti-diabetic agents. Three different alignment strategies were used for the compounds. The best CoMFA and CoMSIA models were obtained by means of Distill rigid body alignment of training and test sets, and found statistically significant with cross-validated coefficients (q ²) of 0.602 and 0.618, respectively, and conventional coefficients (r ²) of 0.905 and 0.902, respectively. Both models were validated by a test set of 36 compounds giving satisfactory predicted correlation coefficients (r ² _pred) of 0.622 and 0.584 for CoMFA and CoMSIA models, respectively. A comparison was made with earlier 3D QSAR study on SGLT2 inhibitors, which shows that our 3D QSAR models are better than earlier models to predict good inhibitory activity. CoMFA and CoMSIA models generated in this work can provide useful information to design new compounds and helped in prediction of activity prior to synthesis.     

3D QSAR Modeling and Molecular Docking Studies on a Series of Triazole Analogues as Antibacterial Agents

The 3D QSAR analysis using the comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) techniques is performed on novel nalidixic acid based 1,2,4-triazole derivatives suggested earlier as antibacterial agents. The CoMFA and CoMSIA models employed for a training set of 28 compounds gives reliable values of Q² (0.53 and 0.52, respectively) and R² (0.79 and 0.85, respectively). The contour maps produced by the CoMFA and CoMSIA models are used to determine a three-dimensional quantitative structure-activity relationship. Based on the 3D QSAR contours new molecules with high predicted activities are designed. In addition, surflex-docking is performed to confirm the stability of predicted molecules in the receptor.     

3D-QSAR and docking studies on 2-arylbenzoxazole and linker-Y transthyretin amyloidogenesis inhibitors

Transthyretin (TTR), a plasma protein with a tetramer structure, could form amyloid fibril associated with several human diseases through the dissociation of tetramer and the misfolding of monomer. These amyloidogenesis can be inhibited by small molecules which bind to the central channel of TTR. A number of small molecules like 2-arylbenzoxazoles (ABZ) analogues are proposed as promising therapeutic strategy to treat amyloidosis. In this work, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) three-dimensional quantitative structure-activity relationship (3D-QSAR) and docking studies were performed on series of 2-arylbenzoxazoles (ABZ) and linker-Y analogues to investigate the inhibitory activities of TTR amyloidogenesis at atomic level. Significant correlation coefficients for ABZ series (CoMFA, r ² = 0.877, q ² = 0.431; CoMSIA, r ² = 0.836, q ² = 0.447) and those for linker-Y series (CoMFA, r ² = 0.828, q ² = 0.522; CoMSIA, r ² = 0.800, q ² = 0.493) were obtained, and the generated models were validated using test sets. In addition, docking studies on 6 compounds binding to TTR were performed to analyze the forward or reverse binding mode and interactions between molecules and TTR. These results from 3D-QSAR and docking studies have great significance for designing novel TTR amyloidogenesis inhibitors in the future.     

Computer-aided molecular design of (E)-N-Aryl-2-ethene-sulfonamide analogues as microtubule targeted agents in prostate cancer

Microtubules are tube-shaped, filamentous and cytoskeletal proteins that are essential in all eukaryotic cells. Microtubule is an attractive and promising target for anticancer agents. In this study, three-dimensional quantitative structure activity relationships (3D-QSAR) including comparative molecular field analysis, CoMFA, and comparative molecular similarity indices analysis, CoMSIA, were performed on a set of 45 (E)-N-Aryl-2-ethene-sulfonamide analogues as microtubule-targeted anti-prostate cancer agents. Automated grid potential analysis, AutoGPA module in Molecular Operating Environment 2009.10 (MOE) as a new 3D-QSAR approach with the pharmacophore-based alignment was carried out on the same dataset. AutoGPA-based 3D-QSAR model yielded better prediction parameters than CoMFA and CoMSIA. Based on the contour maps generated from the models, some key features were identified in (E)-N-Aryl-2-arylethene-sulfonamide analogues that were responsible for the anti-cancer activity. Virtual screening was performed based on pharmacophore modeling and molecular docking to identify the new inhibitors from ZINC database. Seven top ranked compounds were found based on Gold score fitness function. In silico ADMET studies were performed on compounds retrieved from virtual screening in compliance with the standard ranges.     

Design,synthesis and computational approach to study novel pyrrole scaffolds as active inhibitors of enoyl ACP reductase (InhA) and Mycobacterium tuberculosis antagonists

Novel 54 pyrrolyl acetohydrazide analogues were designed, synthesized and screened for antitubercular activity against InhA. Enoyl-ACP reductase/InhA one of the significant enzymes implicated in type II FAS (fatty acid synthase) biosynthetic pathway of bacterial outer cell membrane, in addition Mycobacterium tuberculosis H37Rv inhibition potency has proven to be one of the most promising drug target used for designing and testing against TB. In silico molecular modeling was achieved using Surflex-docking method to recognize important binding sites of the enoyl-ACP reductase. 3D-QSAR studies like CoMFA and CoMSIA approaches were studied to create 3D-QSAR depictions for InhA inhibitors. Based on docking results the synthesized molecules are oriented towards the core of the active site. The pyrolyl acetohydrazides exhibited one or two H-bonding connections with InhA enzyme. Molecule 8l (MIC 0.4 μg/mL, InhA- 70% at 50 μM) showed H-bonding connections with Tyr158 and NAD⁺ in a similar mode to that of ligand pyrrolidine carboxamide. The tested molecules also showed good antibacterial (Escherichia coli, Staphylococcus aureus) activity (MIC 0.4–25 μg/mL), while the study against A549 lung adenocarcinoma cell line confirms nontoxic nature of the reported molecules. The QSAR model from CoMFA and CoMSIA through the database configuration showed the most excellent data. The prognostic ability of CoMFA and CoMSIA representations was developed by a test set of 15 molecules that produced cross validated correlation coefficients (q²) of 0.642 and 0.701, respectively. This study gives insight into structural requirement needed for the development of more active InhA inhibitors through in silico approach.     

Homology modeling and 3D–QSAR study of benzhydrylpiperazine δ opioid receptor agonists

The binding affinity of a series of benzhydrylpiperazine δ opioid receptor agonists were pooled and evaluated by using 3D-QSAR and homology modeling/molecular docking methods. Ligand-based CoMFA and CoMSIA 3D-QSAR analyses with 46 compounds were performed on benzhydrylpiperazine analogues by taking the most active compound BW373U86 as the template. The models were generated successfully with q² value of 0.508 and r² value of 0.964 for CoMFA, and q² value of 0.530 and r² value of 0.927 for CoMSIA. The predictive capabilities of the two models were validated on the test set with R²_pred value of 0.720 and 0.814, respectively. The CoMSIA model appeared to work better in this case. A homology model of active form of δ opioid receptor was established by Swiss-Model using a reported crystal structure of active μ opioid receptor as a template, and was further optimized using nanosecond scale molecular dynamics simulation. The most active compound BW373U86 was docked to the active site of δ opioid receptor and the lowest energy binding pose was then used to identify binding residues such as s Gln105, Lys108, Leu125, Asp128, Tyr129, Leu200, Met132, Met199, Lys214, Trp274, Ile277, Ile304 and Tyr308. The docking and 3D-QSAR results showed that hydrogen bond and hydrophobic interactions played major roles in ligand-receptor interactions. Our results highlight that an approach combining structure-based homology modeling/molecular docking and ligand-based 3D-QSAR methods could be useful in designing of new opioid receptor agonists.     

Computational Analysis of CRTh2 receptor antagonist: A Ligand-based CoMFA and CoMSIA approach

CRTh2 receptor is an important mediator of inflammatory effects and has attracted much attention as a therapeutic target for the treatment of conditions such as asthma, COPD, allergic rhinitis and atopic dermatitis. In pursuit of better CRTh2 receptor antagonist agents, 3D-QSAR studies were performed on a series of 2-(2-(benzylthio)-1H-benzo[d]imidazol-1-yl) acetic acids. There is no crystal structure information available on this protein; hence in this work, ligand-based comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were performed by atom by atom matching alignment using systematic search and simulated annealing methods. The 3D-QSAR models were generated with 10 different combinations of test and training set molecules, since the robustness and predictive ability of the model is very important. We have generated 20 models for CoMFA and 100 models for CoMSIA based on two different alignments. Each model was validated with statistical cut off values such as q² > 0.4, r² > 0.5 and r²_pred > 0.5. Based on better q² and r²_pred values, the best predictions were obtained for the CoMFA (model 5 q² = 0.488, r²_pred = 0.732), and CoMSIA (model 45 q² = 0.525, r²_pred = 0.883) from systematic search conformation alignment. The high correlation between the cross-validated/predicted and experimental activities of a test set revealed that the CoMFA and CoMSIA models were robust. Statistical parameters from the generated QSAR models indicated the data is well fitted and have high predictive ability. The generated models suggest that steric, electrostatic, hydrophobic, H-bond donor and acceptor parameters are important for activity. Our study serves as a guide for further experimental investigations on the synthesis of new CRTh2 antagonist.     

2(1H)-喹啉-2,4-二酮类化合物抗小麦锈病的3D-QSAR研究

用比较分子力场分析(CoMFA)方法和比较分子相似性指数分析(CoMSIA)方法研究了21个2(1H)-喹啉-2,4-二酮类化合物抗小麦锈病的三维定量构效关系(3D-QSAR),发现用CoMFA方法可以找到最佳的3D-QSAR模型,并通过量子化学从头计算的方法研究了不同活性化合物的前线轨道及静电势分布图的差异.所得构效关系模型为发现更高活性的化合物提供理论指导.     

3D-QSAR studies of dihydropyrazole and dihydropyrrole derivatives as inhibitors of human mitotic kinesin Eg5 based on molecular docking

Human mitotic kinesin Eg5 plays an essential role in mitoses and is an interesting drug target against cancer. To find the correlation between Eg5 and its inhibitors, structure-based 3D-quantitative structure-activity relationship (QSAR) studies were performed on a series of dihydropyrazole and dihydropyrrole derivatives using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) methods. Based on the LigandFit docking results, predictive 3D-QSAR models were established, with cross-validated coefficient values (q2) up to 0.798 for CoMFA and 0.848 for CoMSIA, respectively. Furthermore, the CoMFA and CoMSIA models were mapped back to the binding sites of Eg5, which could provide a better understanding of vital interactions between the inhibitors and the kinase. Ligands binding in hydrophobic part of the inhibitor-binding pocket were found to be crucial for potent ligand binding and kinases selectivity. The analyses may be used to design more potent EG5 inhibitors and predict their activities prior to synthesis.     

CoMFA and CoMSIA studies on a new series of xanthone derivatives against the oral human epidermoid carcinoma (KB) cancer cell line

Abstract  A new series of xanthone derivatives against the oral human epidermoid carcinoma (KB) cancer cell line is examined to determine the relationship between the structural properties and the biological activity of these compounds—the 3-D quantitative structure–activity relationship (3D-QSAR)—using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). The best CoMFA and CoMSIA models were obtained using the atom-based alignment of 33 compounds, 22 training compounds and 11 tested compounds, and these give desirable statistics; those for the CoMFA standard model were: r _cv² = 0.691, r ² = 0.998, S _press = 0.178, s = 0.014 and F = 1080.765, while CoMSIA combined steric, electrostatic, hydrophobic and hydrogen-bond acceptor fields: r _cv² = 0.600, r ² = 0.988, S _press = 0.206, s = 0.034 and F = 284.433. The 3D-QSAR models calculated satisfactory test set activities. The 3D-QSAR contour plots correlated strongly with the experimental data for the binding topology. For this reason, these results would be beneficial for predicting affinities with the compounds of interest, and they are advantageous for guiding the design and synthesis of new and more effective anticancer agents. Graphical abstract   A new and more effective anticancer agent of xanthone derivatives against the oral human epidermoid carcinoma (KB) cell line, as investigated by CoMFA and CoMSIA analysis     

CoMFA and CoMSIA studies on a new series of xanthone derivatives against the oral human epidermoid carcinoma (KB) cancer cell line

Abstract  

A new series of xanthone derivatives against the oral human epidermoid carcinoma (KB) cancer cell line is examined to determine the relationship between the structural properties and the biological activity of these compounds—the 3-D quantitative structure–activity relationship (3D-QSAR)—using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). The best CoMFA and CoMSIA models were obtained using the atom-based alignment of 33 compounds, 22 training compounds and 11 tested compounds, and these give desirable statistics; those for the CoMFA standard model were: r _cv² = 0.691, r ² = 0.998, S _press = 0.178, s = 0.014 and F = 1080.765, while CoMSIA combined steric, electrostatic, hydrophobic and hydrogen-bond acceptor fields: r _cv² = 0.600, r ² = 0.988, S _press = 0.206, s = 0.034 and F = 284.433. The 3D-QSAR models calculated satisfactory test set activities. The 3D-QSAR contour plots correlated strongly with the experimental data for the binding topology. For this reason, these results would be beneficial for predicting affinities with the compounds of interest, and they are advantageous for guiding the design and synthesis of new and more effective anticancer agents.     

Improved 3D-QSAR prediction by multiple-conformational alignment: A case study on PTP1B inhibitors

Three-dimension quantitative structure activity relationship (3D-QSAR) was one of the major statistical techniques to investigate the correlation of biological activity with structural properties of candidate molecules, and the accuracy of statistic greatly depended on molecular alignment methodology. Exhaustive conformational search and successful conformational superposition could extremely improve the predictive accuracy of QSAR modeling. In this work, we proposed a solution to optimize QSAR prediction by multiple-conformational alignment methods, with a set of 40 flexible PTP1B inhibitors as case study. Three different molecular alignment methods were used for the development of 3D-QSAR models listed as following: (1) docking-based alignment (DBA); (2) pharmacophore-based alignment (PBA) and (3) co-crystallized conformer-based alignment (CCBA). Among these three alignments, it was indicated that the CCBA was the best and the fastest strategy in 3D-QSAR development, with the square correlation coefficient (r²) and cross-validated squared correlation coefficient (q²) of comparative molecular field analysis (CoMFA) were 0.992 and 0.694; the r² and q² of comparative molecular similarity indices analysis (CoMSIA) were 0.972 and 0.603, respectively. The alignment methodologies used here not only generated a robust QSAR model with useful molecular field contour maps for designing novel PTP1B inhibitors, but also provided a solution for constructing accurate 3D-QSAR model for various disease targets. Undoubtedly, such attempt in QSAR analysis would greatly help us to understand essential structural features of inhibitors required by its target, and so as to discover more promising chemical derivatives.     

In silico design novel (5-imidazol-2-yl-4-phenylpyrimidin-2-yl)[2-(2-pyridylamino)ethyl]amine derivatives as inhibitors for glycogen synthase kinase 3 based on 3D-QSAR,molecular docking and molecular dynamics simulation

Glycogen Synthase Kinase 3 (GSK-3) is a member of cellular kinase with various functions, such as glucose regulation, cellular differentiation, neuronal function and cell apoptosis. It has been proved as an important therapeutic target in type 2 diabetes mellitus and Alzheimer's disease. To better understand their structure–activity relationships and mechanism of action, an integrated computational study, including three dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking, and molecular dynamics (MD), was performed on 79 (5-Imidazol-2-yl-4-phenylpyrimidin-2-yl)[2-(2-pyridylamino)ethyl]amine GSK-3 inhibitors. In this paper, we constructed 3D-QSAR using comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) method. The results showed that the CoMFA model (q ² = 0.743, r² = 0.980) and the CoMSIA model (q² = 0.813, r² = 0.976) had stable and reliable predictive ability. The electrostatic and H-bond donor fields play important roles in the models. The contour maps of the model visually showed the relationship between the activity of compounds and their three-dimensional structure. Molecular docking was used to identify the key amino acid residues at the active site of GSK-3 and explore its binding mode with ligands. Based on 3D-QSAR models, contour maps and the binding feature between GSK-3 and inhibitor, we designed 10 novel compounds with good potential activity and ADME/T profile. Molecular dynamics simulation results validated that Ile62, Val70 and Lys85 located in the active site play a key role for GSK-3 complexed with inhibitors. These results might provide important information for designing GSK-3 inhibitors with high activity.     

Development of docking-based 3D QSAR models for the design of substituted quinoline derivatives as human dihydroorotate dehydrogenase (hDHODH) inhibitors

Abstract

This study has investigated docking-based 3D quantitative structure–activity relationships (QSARs) for a range of quinoline carboxylic acid derivatives by comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). A docking study has shown that most of the compounds formed H-bonds with Arg136 and Gln47, which have already been shown to be essential for the binding of ligands at the active site of the hydroorotate dehydrogenase adenovirus (hDHODH). Bioactive conformations of all the molecules obtained from the docking study were used for the 3D QSAR study. The best CoMFA and CoMSIA models were obtained for the training set and were found to be statistically significant, with cross-validated coefficients (q² ) of 0.672 and 0.613, r² _cv of 0.635 and 0.598 and coefficients of determination (r² ) of 0.963 and 0.896, respectively. Both models were validated by a test set of 15 compounds, giving satisfactory predicted correlation coefficients (r² _pred) of 0.824 and 0.793 for the CoMFA and CoMSIA models, respectively. From the docking-based 3D QSAR study we designed 34 novel quinoline-based compounds and performed structure-based virtual screening. Finally, in silico pharmacokinetics and toxicities were predicted for 24 of the best docked molecules. The study provides valuable information for the understanding of interactions between hDHODH and the novel compounds.     

Docking-based 3D-QSAR study of pyridyl aminothiazole derivatives as checkpoint kinase 1 inhibitors

Checkpoint kinase 1 (Chk1) is a promising target for the design of novel anticancer agents. In the present work, molecular docking simulations and three-dimensional quantitative structure–activity relationship (3D-QSAR) studies were performed on pyridyl aminothiazole derivatives as Chk1 inhibitors. AutoDock was used to determine the probable binding conformations of all the compounds inside the active site of Chk1. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) models were developed based on the docking conformations and alignments. The CoMFA model produced statistically significant results with a cross-validated correlation coefficient (q²) of 0.608 and a coefficient of determination (r²) of 0.972. The reliable CoMSIA model with q² of 0.662 and r² of 0.970 was obtained from the combination of steric, electrostatic and hydrogen bond acceptor fields. The predictive power of the models were assessed using an external test set of 14 compounds and showed reasonable external predictabilities (r²_pred) of 0.668 and 0.641 for CoMFA and CoMSIA models, respectively. The models were further evaluated by leave-ten-out cross-validation, bootstrapping and progressive scrambling analyses. The study provides valuable information about the key structural elements that are required in the rational design of potential drug candidates of this class of Chk1 inhibitors.