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.     

3D-QSAR studies and molecular docking on [5-(4-amino-1<Emphasis Type="Italic">H</Emphasis>-benzoimidazol-2-yl)-furan-2-yl]-phosphonic acid derivatives as fructose-1,6-biphophatase inhibitors

Fructose-1,6-biphophatase has been regarded as a novel therapeutic target for the treatment of type 2 diabetes mellitus (T2DM). 3D-QSAR and docking studies were performed on a series of [5-(4-amino-1H-benzoimidazol-2-yl)-furan-2-yl]-phosphonic acid derivatives as fructose-1,6-biphophatase inhibitors. The CoMFA and CoMSIA models using thirty-seven molecules in the training set gave r _cv² values of 0.614 and 0.598, r ² values of 0.950 and 0.928, respectively. The external validation indicated that our CoMFA and CoMSIA models possessed high predictive powers with r ₀² values of 0.994 and 0.994, r _m² values of 0.751 and 0.690, respectively. Molecular docking studies revealed that a phosphonic group was essential for binding to the receptor, and some key features were also identified. A set of forty new analogues were designed by utilizing the results revealed in the present study, and were predicted with significantly improved potencies in the developed models. The findings can be quite useful to aid the designing of new fructose-1,6-biphophatase inhibitors with improved biological response.     

吲哚并喹啉衍生物与G-四链体相互作用的分子模拟研究

G-四链体是富含鸟嘌呤碱基的DNA序列通过氢键相互作用形成的四链螺旋结构. 通过小分子化合物诱导与稳定端粒G-四链体从而抑制端粒酶活性是一种新的抗癌策略. 为了研究一系列吲哚并喹啉衍生物与端粒G-四链体的相互作用, 探究其相互作用模式, 从而为实现基于G-四链体结构的药物合理设计提供依据, 使用分子对接的方法构建了吲哚并喹啉衍生物与G-四链体复合物结构, 在此基础上进行分子动力学模拟, 并使用线性相互作用能(LIE)方法计算了化合物与G-四链体的结合自由能. 结果表明: 化合物与G-四链体的主要相互作用方式由氢键、静电与π-π堆积作用构成, 侧链末端基团类型和侧链的长短是影响相互作用强弱的重要因素. 通过LIE方法计算的结合自由能与实验结果基本吻合, 相关度达到r²＝0.79. 并且, 基于预测的结合模式, 总结了拥有更高活性的新型吲哚并喹啉衍生物应具有的几个结构特征.     

Identification of potential CRAC channel inhibitors: Pharmacophore mapping, 3D-QSAR modelling,and molecular docking approach

Upregulation of store-operated Ca²⁺ influx via ORAI1, an integral component of the CRAC channel, is responsible for abnormal cytokine release in active rheumatoid arthritis, and therefore ORAI1 has been proposed as an attractive molecular target. In this study, we attempted to predict the mechanical insights of ORAI1 inhibitors through pharmacophore modelling, 3D-QSAR, molecular docking and free energy analysis. Various hypotheses of pharmacophores were generated and from that, a pharmacophore hypothesis with two hydrogen bond acceptors, one hydrogen bond donor and two aromatic rings (AADRR) resulted in a statistically significant 3D-QSAR model (r² = 0.84 and q² = 0.74). We believe that the obtained statistical model is a reliable QSAR model for the diverse dataset of inhibitors against the IL-2 production assay. The visualization of contours in active and inactive compounds generated from the 3D-QSAR models and molecular docking studies revealed major interaction with GLN108, HIS113 and ASP114, and interestingly, these residues are located near the Ca²⁺ selectivity filter region. Free energy binding analysis revealed that Coulomb energy, van der Waals energy and non-polar solvation terms are more favourable for ligand binding. Thus, the present study provides the physical and chemical requirements for the development of novel ORAI1 inhibitors with improved biological activity.     

Identification of novel PfDHODH inhibitors as antimalarial agents via pharmacophore-based virtual screening followed by molecular docking and in vivo antimalarial activity

Plasmodium falciparum dihydroorotate dehydrogenase (PfDHODH) catalyses the fourth reaction of de novo pyrimidine biosynthesis in parasites, and represents an important target for the treatment of malaria. In this study, we describe pharmacophore-based virtual screening combined with docking study and biological evaluation as a rational strategy for identification of novel hits as antimalarial agents. Pharmacophore models were established from known PfDHODH inhibitors using the GALAHAD module with IC₅₀ values ranging from 0.033 μM to 142 μM. The best pharmacophore model consisted of three hydrogen bond acceptor, one hydrogen bond donor and one hydrophobic features. The pharmacophore models were validated through receiver operating characteristic and Günere–Henry scoring methods. The best pharmacophore model as a 3D search query was searched against the IBS database. Several compounds with different structures (scaffolds) were retrieved as hit molecules. Among these compounds, those with a Q_FIT value of more than 81 were docked in the PfDHODH enzyme to further explore the binding modes of these compounds. In silico pharmacokinetic and toxicities were predicted for the best docked molecules. Finally, the identified hits were evaluated in vivo for their antimalarial activity in a parasite inhibition assay. The hits reported here showed good potential to become novel antimalarial agents.     

Molecular modeling studies of human immunodeficiency virus type 1 protease inhibitors using three‐dimensional quantitative structure‐activity relationship,virtual screening,and docking simulations

In this paper, two 3‐dimensional quantitative structure‐activity relationship models for 60 human immunodeficiency virus (HIV)‐1 protease inhibitors were established using random sampling analysis on molecular surface and translocation comparative molecular field vector analysis (Topomer CoMFA). The non–cross‐validation (r²), cross‐validation (q²), correlation coefficient of external validation (Q²_ext), and F of 2 models were 0.94, 0.80, 0.79, and 198.84 and 0.94, 0.72, 0.75, and 208.53, respectively. The results indicated that 2 models were reasonable and had good prediction ability. Topomer Search was used to search R groups in the ZINC database, 20 new compounds were designed, and the Topomer CoMFA model was used to predicate the biological activity. The results showed that 18 new compounds were more active than the template molecule. So the Topomer Search is effective in screening and can guide the design of new HIV/AIDS drugs. The mechanism of action was studied by molecular docking, and it showed that the protease inhibitors and Ile50, Asp25, and Arg8 sites of HIV‐1 protease have interactions. These results have provided an insight for the design of new potent inhibitors of HIV‐1 protease.     

Docking-based CoMFA and CoMSIA study of azaindole carboxylic acid derivatives as promising HIV-1 integrase inhibitors

Three-dimensional quantitative structure–activity relationship (3D-QSAR) studies were performed based on a series of azaindole carboxylic acid derivatives that had previously been reported as promising HIV-1 integrase inhibitors. Docking studies to explore the binding mode were performed based on the highly active molecule 36. The best docked conformation of molecule 36 was used as template for alignment. The comparative molecular field analysis (CoMFA) model (including steric and electrostatic fields) yielded the cross validation q ² = 0.655, non-cross validation r ² = 0.989 and predictive r ² _pred = 0.979. The best comparative molecular similarity indices analysis (CoMSIA) model (including steric, electrostatic, hydrophobic and hydrogen-bond acceptor fields) yielded the cross validation q ² = 0.719, non-cross validation r ² = 0.992 and predictive r ² _pred = 0.953. A series of new azaindole carboxylic acid derivatives were designed and the HIV-1 integrase inhibitory activities of these designed compounds were predicted based on the CoMFA and CoMSIA models.     

A three-dimensional chemical phase pharmacophore mapping,QSAR modelling and electronic feature analysis of benzofuran salicylic acid derivatives as LYP inhibitors

Lymphoid tyrosine phosphatase (LYP), encoded by the PTPN22 gene, has a critical negative regulatory role in T-cell antigen receptor (TCR) and emerged as a promising drug target for human autoimmune diseases. A five-point pharmacophore with two hydrogen bond acceptors, one hydrogen bond donor and two aromatic ring features was generated for a series of benzofuran salicylic acid derivatives as LYP inhibitors in order to elucidate their anti-autoimmune activity. The generated pharmacophore yielded a significant 3D-QSAR model with r² of 0.9146 for a training set of 27 compounds. The model also showed excellent predictive power with Q² of 0.7068 for a test set of eight compounds. The investigation of the 3D-QSAR model has revealed the structural insights which could lead to more potent analogues. The most active and inactive compounds were further subjected to electronic structure analysis using density functional theory (DFT) at B3LYP/3?21^?G level to support the 3D-QSAR predictions. The results obtained from this study are expected to be useful in the proficient design and development of benzofuran salicylic acid derivatives as inhibitors of LYP.     

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.     

Molecular modelling studies on adamantane-based Ebola virus GP-1 inhibitors using docking,pharmacophore and 3D-QSAR

The pathogenic Ebola virus (EBOV) causes a potential health risk and global spread. To date, few drugs are available for the treatment of Ebola virus disease (EVD) that allow researchers to use computational methods for designing potential drugs. The developed PHASE-based common six-point pharmacophore hypothesis (AADHPR_1) showed the necessity of two hydrogen bond acceptor features, one hydrogen bond donor feature, one hydrophobic group feature, one positively ionizable and one aromatic ring feature for further designing. We developed best 3D-QSAR models with high regression coefficients for the training (r²>0.82) and test (Q²>0.5) sets for both atoms-based and field-based 3D-QSAR models. The molecule 1A-4 (docking score = –4.711 kcal/mol) was obtained as best docked (SP mode) on Ebola virus envelope glycoprotein (PDB ID-3CSY) as compared with the standards oseltamivir (docking score = –4.39 kcal/mol) and zanamivir (docking score = –3.392 kcal/mol). The obtained ZINC hit ZINC58935541 showed a good docking score of –4.892 kcal/mol. The ZINC58935541 molecule also showed a strong binding affinity towards the receptor cavity of Ebola virus envelope glycoprotein when simulated for 1.2 ns. The good QikProp parameters reflect the fact that this molecule, upon optimization into a lead, might become a good candidate for the treatment of EVD.     

Molecular dynamics and pharmacophore modelling studies of different subtype (ALK and EGFR (T790M)) inhibitors in NSCLC

Extensively validated 3D pharmacophore models for ALK (anaplastic lymphoma kinase) and EGFR (T790M) (epithelial growth factor receptor with acquired secondary mutation) were developed. The pharmacophore model for ALK (r² = 0.96, q² = 0.692) suggested that two hydrogen bond acceptors and three hydrophobic groups arranged in 3-D space are essential for the binding affinity of ALK inhibitors. Similarly, the pharmacophore model for EGFR (T790M) (r² = 0.92, q² = 0.72) suggested that the presence of a hydrogen bond acceptor, two hydrogen bond donors and a hydrophobic group plays vital role in binding of an inhibitor of EGFR (T790M). These pharmacophore models allowed searches for novel ALK and EGFR (T790M) dual inhibitors from multiconformer 3D databases (Asinex, Chembridge and Maybridge). Finally, the eight best hits were selected for molecular dynamics simulation, to study the stability of their complexes with both proteins and final binding orientations of these molecules. After molecular dynamics simulations, one hit has been predicted to possess good binding affinity for both ALK and EGFR (T790M), which can be further investigated for its experimental in-vitro/in-vivo activities.     

Monte Carlo method based QSAR modelling of natural lipase inhibitors using hybrid optimal descriptors

Obesity is one of the most provoking health burdens in the developed countries. One of the strategies to prevent obesity is the inhibition of pancreatic lipase enzyme. The aim of this study was to build QSAR models for natural lipase inhibitors by using the Monte Carlo method. The molecular structures were represented by the simplified molecular input line entry system (SMILES) notation and molecular graphs. Three sets – training, calibration and test set of three splits – were examined and validated. Statistical quality of all the described models was very good. The best QSAR model showed the following statistical parameters: r² = 0.864 and Q² = 0.836 for the test set and r² = 0.824 and Q² = 0.819 for the validation set. Structural attributes for increasing and decreasing the activity (expressed as pIC₅₀) were also defined. Using defined structural attributes, the design of new potential lipase inhibitors is also presented. Additionally, a molecular docking study was performed for the determination of binding modes of designed molecules.     

Use of 3D QSAR to investigate the mode of binding of pyrazinones to HIV-1 RT

Abstract  

Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) based on the docked conformation were performed for 24 pyrazinone derivatives. All compounds were docked into the wild-type HIV-1 RT binding pocket and the lowest-energy docked configurations were used to construct the 3D QSAR models. The CoMFA and CoMSIA models enable good prediction of inhibition by the pyrazinones, with r_\textcv² r_{\text{cv}}^{2}  = 0.703 and 0.735. Results obtained from CoMFA and CoMSIA based on the docking conformation of the pyrazinones are, therefore, powerful means of elucidating the mode of binding of pyrazinones and suggesting the design of new potent NNRTIs.     

Molecular Modeling and 3D-QSAR Studies in 2-Aziridinyl-and 2,3-Bis(Aziridinyl)-1,4-Naphthoquinonyl Sulfonate and Acylate Derivatives as Potential Antimalarial Agents

Abstract

Malaria is still continuing to be one of the most dreadful diseases of the tropical countries particularly due to the development of resistance to the existing antimalarials. From observed, antimalarial activity of 2-aziridinyl- and 2,3-bis(aziridinyl)-1,4-naphthoquinonyl sulfonate and acylate derivatives acting through redox cycling mechanism, molecular modeling and three dimensional-quantitative structure activity relationship (3D-QSAR) studies have been carried out on a set of 63 compounds to identify important pharmacophors. Among several 3D-QSAR models generated, three models with correlation coefficient r > 0.82, match > 0.60 and chance = 0.00 have shown two common biophoric sites: one being the oxygen atom at position 1 of the naphthoquinone ring in terms of π-population, charge and electron donating ability while the second being the center of the phenyl ring in terms of its 6π-electrons. In addition to these sites, the models also share two common secondary sites: one positively contributing H-acceptor site while the second site contributing negatively in terms of steric refractivity. All these models showed good agreement between the experimental, calculated and predicted antimalarial activities.