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.     

Computational explorations to gain insight into the structural features of TNF-α receptor I inhibitors

TNF-α is a crucial cytokine in the process of inflammatory diseases. The adverse effect of TNF-α is mostly mediated by interaction of TNF-α with TNF-α receptor type I (TNFR1); therefore, discovery of molecules which can bind to TNFR1 preventing TNF-α-receptor complex formation would be of great interest. In the current study, using GRID/GOLPE program, a 3D-QSAR study was conducted on a series of synthetic TNFR1 binders, which resulted in a 3D-QSAR model with appropriate power of predictivity in internal (r²?=?0.94 and q²_LOO?=?0.74) and external (r²?=?0.66 and SDEP?=?0.42) validations. The structural features of TNFR1 inhibitors essential for exerting activity were explored by analyzing the contour maps of the 3D-QSAR model showing that steric interactions and hydrogen bonds are responsible for exerting TNFR1 inhibitory activity. To propose potential chemical entities for TNFR1 inhibition, PubChem database was searched and the selected compounds were virtually tested for anti-TNFR1 activity using the generated model, resulting in two potential anti-TNFR1 compounds. Finally, the possible interactions of the compounds with TNFR1 were investigated using docking studies. The findings in the current work can pave the way for designing more potent anti-TNFR1 inhibitors.     

CoMFA,CoMSIA, Topomer CoMFA,HQSAR, molecular docking and molecular dynamics simulations study of triazine morpholino derivatives as mTOR inhibitors for the treatment of breast cancer

mTOR has become a promising target for many types of cancer like breast, lung and renal cell carcinoma. CoMFA, CoMSIA, Topomer CoMFA and HQSAR were performed on the series of 39 triazine morpholino derivatives. CoMFA analysis showed q² value of 0.735, r²_cv value of 0.722 and r²_pred value of 0.769. CoMSIA analysis (SEHD) showed q² value of 0.761, r²_cv value of 0.775 and r²_pred value of 0.651. Topomer CoMFA analysis showed q² value of 0.693, r² (conventional correlation coefficient) value of 0.940 and r²_pred value of 0.720. HQSAR analysis showed q²,r²and r²_pred values of 0.694, 0.920 and 0.750, respectively. HQSAR analysis with the combination of atomic number (A), bond type (B) and atomic connections showed q² and r² values of 0.655 and 0.891, respectively. Contour maps from all studies provided significant insights. Molecular docking studies with molecular dynamics simulations were carried out on the highly potent compound 36. Furthermore, four acridine derivatives were designed and docking results of these designed compounds showed the same interactions as that of the standard PI-103 which proved the efficiency of 3D-QSAR and MD/MS study. In future, this study might be useful prior to synthesis for the designing of novel mTOR inhibitors.     

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.     

A comparative molecular field analysis of cytochrome P450 2A5 and 2A6 inhibitors

Structure-activity relationships of 23 P450 2A5 and 2A6 inhibitors were analysed using the CoMFA [1] and GOLPE/GRID with smart region definition (SRD) [2]. The predictive power of the resulting models was validated using five compounds not belonging to the model set. All models have high internal and external predictive power and resulting 3D-QSAR models are supporting each other. Both Sybyl and GOLPE highlight properties near lactone moiety to be important for 2A5 and 2A6 inhibition. Another important feature for pIC₅₀ was the size of the substituent in the 7-positon of coumarin. The models suggest that the 2A5 binding site is larger that that of 2A6 due to larger steric regions in the CoMFA coefficient maps and corresponding GOLPE maps. In addition, the maps reveal that 2A6 disfavours negative charge near the lactone moiety of coumarin.     

Receptor-based 3D QSAR analysis of estrogen receptor ligands – merging the accuracy of receptor-based alignments with the computational efficiency of ligand-based methods

One of the major challenges in computational approaches to drug design is the accurate prediction of binding affinity of biomolecules. In the present study several prediction methods for a published set of estrogen receptor ligands are investigated and compared. The binding modes of 30 ligands were determined using the docking program AutoDock and were compared with available X-ray structures of estrogen receptor-ligand complexes. On the basis of the docking results an interaction energy-based model, which uses the information of the whole ligand-receptor complex, was generated. Several parameters were modified in order to analyze their influence onto the correlation between binding affinities and calculated ligand-receptor interaction energies. The highest correlation coefficient (r ² = 0.617, q ² _LOO = 0.570) was obtained considering protein flexibility during the interaction energy evaluation. The second prediction method uses a combination of receptor-based and 3D quantitative structure-activity relationships (3D QSAR) methods. The ligand alignment obtained from the docking simulations was taken as basis for a comparative field analysis applying the GRID/GOLPE program. Using the interaction field derived with a water probe and applying the smart region definition (SRD) variable selection, a significant and robust model was obtained (r ² = 0.991, q ² _LOO = 0.921). The predictive ability of the established model was further evaluated by using a test set of six additional compounds. The comparison with the generated interaction energy-based model and with a traditional CoMFA model obtained using a ligand-based alignment (r ² = 0.951, q ² _LOO = 0.796) indicates that the combination of receptor-based and 3D QSAR methods is able to improve the quality of the underlying model.     

3D-QSAR (CoMFA,CoMSIA, HQSAR and topomer CoMFA), MD simulations and molecular docking studies on purinylpyridine derivatives as B-Raf inhibitors for the treatment of melanoma cancer

As per the World Health Organization (WHO), cancer is the second most leading cause of death after cardiovascular diseases in worldwide with around 9.88 million total new cases and 1.08 million were observed due to skin cancer in 2018. Amongst two types of skin cancer, progression of melanoma cancer is increasing day by day due to the environmental changes than non-melanoma cancer. Most of B-Raf mutation, specifically B-RafV600E, is responsible for the progression of the melanoma cancer. Here, various 3D-QSAR techniques like comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), molecular hologram QSAR (HQSAR) and topomer CoMFA were used to design novel B-Raf inhibitors by using 28 synthetic B-Raf inhibitors. Except for topomer CoMFA model, remaining models were generated by three different alignment methods in which distil-based alignment method was found best and gave prominent statistical values. After performing N-fold statistical validation, in CoMFA, q², r² and r²_pred values were found to be 0.638, 0.969 and 0.848, respectively. Similarly, q², r² and r²_pred values were found to be 0.796, 0.978 and 0.891 in CoMSIA (SHD) and 0.761, 0.973 and 0.852 in CoMSIA (SH) by N-fold statistical validation. In HQSAR analysis, statistical values were found for q² as 0.984, r² as 0.999 and r²_pred as 0.634 with 97 as best hologram length (BHL). The results of topomer CoMFA showed the q² value of 0.663 and the r² value of 0.967. Important features of purinylpyridine were identified by contour map analysis of all 3D-QSAR techniques, which could be useful to design the novel molecules as B-Raf inhibitors for the treatment of melanoma cancer.

     

Imidazole-containing farnesyltransferase inhibitors: 3D quantitative structure–activity relationships and molecular docking

One of the most promising anticancer and recent antimalarial targets is the heterodimeric zinc-containing protein farnesyltransferase (FT). In this work, we studied a highly diverse series of 192 Abbott-initiated imidazole-containing compounds and their FT inhibitory activities using 3D-QSAR and docking, in order to gain understanding of the interaction of these inhibitors with FT to aid development of a rational strategy for further lead optimization. We report several highly significant and predictive CoMFA and CoMSIA models. The best model, composed of CoMFA steric and electrostatic fields combined with CoMSIA hydrophobic and H-bond acceptor fields, had r ² = 0.878, q ² = 0.630, and r _pred² = 0.614. Docking studies on the statistical outliers revealed that some of them had a different binding mode in the FT active site based on steric bulk and available active site space, explaining why the predicted activities differed from the experimental activities. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

靛玉红类CDK1抑制剂的同源模建、分子对接及3D-QSAR研究

细胞周期蛋白依赖性激酶1的异常表达会导致G2期的停滞及多种肿瘤的发生,故CDK1近年来已成为一个理想的治疗靶点. 本文以细胞分裂调控蛋白2的同源体为模板,同源模建了CDK1的结构,并与靛玉红类小分子抑制剂进行分子对接. 分别运用三种叠合方法进行分子叠合,并在此基础上采用Sybyl 7.1中的比较分子场分析(CoMFA)模块及Discovery Studio 3.0中的三维定量构效关系(3D-QSAR)模块(以下简称为DS)分别建立了3D-QSAR模型. 其中,将分子对接叠合与公共骨架叠合联合运用的叠合方法所得3D-QSAR模型的评价参数是最佳的(CoMFA：q²=0.681,r²=0.909,r_pred.²=0.836; DS：q²=0.579,r²=0.971,r_pred.²=0.795,其中q²为交叉验证系数,r²为非交叉验证系数). 本文的研究结果在对靛玉红类小分子进行结构修饰设计出新的CDK1抑制剂方面,可提供重要的理论基础.     

Comparative QSAR studies using HQSAR,CoMFA, and CoMSIA methods on cyclic sulfone hydroxyethylamines as BACE1 inhibitors

The inhibition of β-secretase (BACE1) is currently the main pharmacological strategy available for Alzheimer’s disease (AD). 2D QSAR and 3D QSAR analysis on some cyclic sulfone hydroxyethylamines inhibitors against β-secretase (IC₅₀: 0.002–2.75 μM) were carried out using hologram QSAR (HQSAR), comparative molecular field analysis (CoMFA), and comparative molecular similarity indices analysis (CoMSIA) methods. The best model based on the training set was generated with a HQSAR q² value of 0.693 and r² value of 0.981; a CoMFA q² value of 0.534 and r² value of 0.913; and a CoMSIA q² value of 0.512 and r² value of 0.973. In order to gain further understand of the vital interactions between cyclic sulfone hydroxyethylamines and the protease, the analysis was performed by combining the CoMFA and CoMSIA field distributions with the active sites of the BACE1. The final QSAR models could be helpful in the design and development of novel active BACE1 inhibitors.     

Molecular docking and 3D-quantitative structure activity relationship analyses of peptidyl vinyl sulfones: <Emphasis Type="Italic">Plasmodium Falciparum</Emphasis> cysteine proteases inhibitors

Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) based on three-dimensional quantitative structure–activity relationship (3D-QSAR) studies were conducted on a series (39 molecules) of peptidyl vinyl sulfone derivatives as potential Plasmodium Falciparum cysteine proteases inhibitors. Two different methods of alignment were employed: (i) a receptor-docked alignment derived from the structure-based docking algorithm GOLD and (ii) a ligand-based alignment using the structure of one of the ligands derived from a crystal structure from the PDB databank. The best predictions were obtained for the receptor-docked alignment with a CoMFA standard model (q ² = 0.696 and r ² = 0.980) and with CoMSIA combined electrostatic, and hydrophobic fields (q ² = 0.711 and r ² = 0.992). Both models were validated by a test set of nine compounds and gave satisfactory predictive r ² _pred values of 0.76 and 0.74, respectively. CoMFA and CoMSIA contour maps were used to identify critical regions where any change in the steric, electrostatic, and hydrophobic fields may affect the inhibitory activity, and to highlight the key structural features required for biological activity. Moreover, the results obtained from 3D-QSAR analyses were superimposed on the Plasmodium Falciparum cysteine proteases active site and the main interactions were studied. The present work provides extremely useful guidelines for future structural modifications of this class of compounds towards the development of superior antimalarials.     

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.     

Various atomic charge calculation schemes of CoMFA on HIF‐1 inhibitors of moracin analogs

Selection of appropriate partial charges in a molecule is crucial to derive good quantitative structure–activity relationship models. In this work, several partial atomic charges were assigned and tested in a comparative molecular field analysis (CoMFA) models. Many CoMFA models were generated for a series of hypoxia inducible factor 1 (HIF‐1) inhibitors using various partial atomic charges including charge equalization, Mülliken population analysis (MPA), natural population analysis, and electrostatic potential (ESP)‐derived charges. These atomic charges were investigated at various theoretical levels such as empirical, semiempirical, Hartree–Fock (HF), and density functional theory (DFT). Among them, Merz‐Singh‐Kollman (MK) ESP‐derived charges at the level of HF/6‐31G* gave the highest predictive q² with experimental pIC₅₀ values. With this charge scheme, a detailed analysis of CoMFA model was performed to understand the electrostatic interactions between ligand and receptor. More elaborate charge calculation schemes such as HF and DFT correlated more strongly with activity than empirical or semiempirical schemes. The choice of optimization methods was important. As geometries were fully optimized at the given levels of theory, the aligned structures were different. They differed considerably, especially for the flexible parts. This was likely the source of the substantial variation of q² values, even when the same steric factor was considered without electrostatic parameters. ESP‐derived charges were most appropriate to describe CoMFA electrostatic interactions among MPA, NBA, and ESP charges. Overall q² values vary considerably (0.8–0.5) depending on the charge schemes applied. The results demonstrate the need to consider more appropriate atomic charges rather than default CoMFA charges. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

QM/MM based 3D QSAR models for potent B-Raf inhibitors

Three dimensional (3D) quantitative structure-activity relationship studies of 37 B-Raf inhibitors, pyrazole-based derivatives, were performed. Based on the co-crystallized compound (PDB ID: 3D4Q), several alignment methods were utilized to derive reliable comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) models. Receptor-guided alignment with quantum mechanics/molecular mechanics (QM/MM) minimization led to the best CoMFA model (q ² = 0.624, r ² = 0.959). With the same alignment, a statistically reliable CoMSIA model with steric, H-bond acceptor, and hydrophobic fields was also derived (q ² = 0.590, r ² = 0.922). Both models were validated with an external test set, which gave satisfactory predictive r ² values of 0.926 and 0.878, respectively. Contour maps from CoMFA and CoMSIA models revealed important structural features responsible for increasing biological activity within the active site and explained the correlation between biological activity and receptor-ligand interactions. New fragments were identified as building blocks which can replace R1-3 groups through combinatorial screening methods. By combining these fragments a compound with a high bioactivity level prediction was found. These results can offer useful information for the design of new B-Raf inhibitors.     

Computational studies on N-phenyl pyrrole derivatives as MmpL3 inhibitors in Mycobacterium tuberculosis

The fight against tuberculosis (TB) is a time immemorial one and the emergence of new drug resistant strains of Mycobacterium tuberculosis keeps throwing new challenges to the scientific community immersed in finding mechanisms to control this dreaded disease. Computer aided drug designing (CADD) is one of the several approaches that can assist in identifying the potent actives against Mycobacterium. In this work, a series of 109 known Mycobacterial membrane proteins large 3 (MmpL3) inhibitors were pooled and atom based 3D QSAR analysis was performed to understand the structural features essential for inhibitory activity against the MmpL3, known to be a key player in transporting substances critical for cell wall integrity of Mycobacterium. The data set employed was randomly split into training set and test set molecules. The training set of 74 molecules was used to derive CoMFA and CoMSIA models that were statistically reliable (CoMFA: q²_loo = 0.53; r²_ncv = 0.93 and CoMSIA: q²_loo = 0.60; r²_ncv = 0.93). The derived models also exhibited good external predictive ability (CoMFA: r²_pred = 0.78 and CoMSIA: r²_pred = 0.79). The results are quite encouraging and information derived from these analyses was applied to design new molecules. The designed molecule showed appreciable predicted activity values and reasonably good ADMET profile. The strategy used in designing new molecules can be pursued in the hunt for new chemical entities targeting MmpL3, expanding the existing arsenal against TB.     

1,2-萘醌类化合物抑制PTP1B的三维定量构效关系研究

蛋白酪氨酸磷酸酶1B (protein tyrosine phosphatase 1B, PTP-1B)是近年来发现的治疗II型糖尿病的新靶点, 1,2-萘醌类化合物对PTP-1B有较好的抑制活性, 具有良好的药用前景. 为了设计出本类化合物抑制效果更好的分子构型, 用比较分子力场分析(CoMFA)和比较分子相似性指数分析(CoMSIA)对该类化合物进行了三维定量构效关系(3D-QSAR)的研究, 并建立了相关的预测模型. 其中, CoMFA模型的交叉验证相关系数(q²)为0.555, 非交叉验证相关系数(r²)为0.991, 标准偏差(SEE)为0.049, F值为564.910. CoMSIA模型的q²为0.558, r²为0.991, SEE为0.050, F值为542.773. 计算结果表明, 获得的CoMFA和CoMSIA模型具有良好的预测能力, 可以应用于指导该类化合物的设计.     

Three-dimensional quantitative structure-activity relationship of 4,5,6,7-tetrahydrothienopyridines analogues as glucose-6-phosphatase inhibitors

A Three-Dimensional Quantitative Structure-activity Relationship (3D-QSAR) model that correlates the biological activities with the chemical structures of a series of Glucose-6-phosphatase inhibitors, exemplified by the 4,5,6,7-tetrahydrothienopyridines derivatives, was established by means of comparative molecular field analysis (CoMFA). The resulting leave-one-out cross-validated value (q²=0.600) and non-cross-validated value (r²=0.956) indicate that the obtained pharmacophore model indeed mimics the steric and electrostatic environment, where inhibitors bind to the enzyme. Furthermore, the developed model also possesses promising predictive ability as discerned by the testing on the external test set. The analysis of the CoMFA contour map, which reveal how steric and electrostatic interactions contribute to inhibitors' bioactivities, provide us with the important information to understand the molecular nature of inhibitor-enzyme interactions and to aid in the design of more potent Glucose-6-phosphatase inhibitors.     

QSAR models of cytochrome P450 enzyme 1A2 inhibitors using CoMFA,CoMSIA and HQSAR

Quantitative structure–activity relationship (QSAR) studies were conducted on an in-house database of cytochrome P450 enzyme 1A2 inhibitors using the comparative molecular field analysis (CoMFA), comparative molecular similarity analysis (CoMSIA) and hologram QSAR (HQSAR) approaches. The database consisted of 36 active molecules featuring varied core structures. The model based on the naphthalene substructure alignment incorporating 19 molecules yielded the best model with a CoMFA cross validation value q² of 0.667 and a Pearson correlation coefficient r² of 0.976; a CoMSIA q² value of 0.616 and r² value of 0.985; and a HQSAR q² value of 0.652 and r² value of 0.917. A second model incorporating 34 molecules aligned using the benzene substructure yielded an acceptable CoMFA model with q² value of 0.5 and r² value of 0.991. Depending on the core structure of the molecule under consideration, new CYP1A2 inhibitors will be designed based on the results from these models.     

Studies of febuxostat analogues as xanthine oxidase inhibitors through 3D-QSAR,Topomer CoMFA and molecular modeling

Xanthine oxidase, a complex molybdoflavoprotein, catalyzes the hydroxylation of xanthine to uric acid, which has emerged as an important target for gout and hyperuricemia. In this work, a combination of molecular modeling methods was performed on a series of febuxostat analogues as xanthine oxidase inhibitors to establish molecular models for new drug design, including three-dimensional quantitative structure–activity relationship, topomer comparative molecular field analysis (CoMFA), molecular docking and molecular dynamic simulations. The optimal CoMFA model yielded a leave-one-out correlation coefficient (q²) of 0.841 and a non-validated correlation coefficient (r²) of 0.985. The respective q² and r² of the best comparative molecular similarity indices analysis (CoMSIA) model were 0.794 and 0.972, respectively. The Topomer CoMFA model provided a q² of 0.915 and an r² of 0.977. 3D contour maps generated from CoMFA and CoMSIA have identified several key features responsible for the inhibition activity. Molecular modeling was taken to further elucidate the proposed binding conformations of the inhibitors to the protein. The obtained results can be served as a useful guideline for designing novel febuxostat derivatives with improved activity against xanthine oxidase.