Molecular Docking and 3D QSAR Research of Indolocarbazole Series as Cyclin-Dependent Kinase Inhibitors

Fifty indolocarbazole series as cyclin-dependent kinase inhibitors (CDKs) are used to establish a threedimensional quantitative structure-activity relationship (3D QSAR) model based on docking conformations resulting from the Topomer comparative molecular field analysis (Topomer CoMFA). The statistic parameters show that the cross-validation (q²), the multiple correlation coefficient of fitting (r²), and external validation statistic (Q_ext²) are 0.953, 0.968, and 0.954, respectively. It is demonstrated that this Topomer CoMFA model has good stability and prediction ability. The methodology of the fragment-based drug design (FBDD) was also used to virtually screen new CDKs by the Topomer Search technology. Four similar substitutional groups selected from the ZINC database were added to the basic scaffold. As a result, 18 new CDKs with high activities were obtained. The template molecule and new designed compounds are used to study the binding relationship between the ligands and the receptor protein with Surflex-Dock. The docking results suggest good binding interactions of the designed compounds with protein. There are several hydrogen bondings between CDKs with amino acid residues of LYS33, LYS89, ASP86, LEU83, GLU81.     

A 3D-QSAR Study of HIV-1 Integrase Inhibitors Using RASMS and Topomer CoMFA

Acquired Immunodeficiency Syndrome(AIDS) is a significant human health threat around the world. Therefore, the study of anti-human immunodeficiency virus(HIV) drug design has become an important task for today's society. In this paper, a three-dimensional quantitative structure-activity relationships study(3 D-QSAR) was conducted on 53 HIV-1 integrase inhibitors(IN) using random sampling analysis on molecular surface(RASMS) and Topomer comparative molecular field analysis(Topomer CoMFA). The multiple correlation coefficients of fitting, cross-validation, and external validation of two models were 0.926, 0.815 and 0.908 and 0.930, 0.726 and 0.855, respectively. The results indicated that two models obtained had both favorable estimation stability and good prediction capability. Topomer Search was used to search appropriate R groups from ZINC database, and 28 new compounds were designed thereby. The Topomer CoMFA model was subsequently used to predict the biological activity of these compounds, showing that 24 of the new compounds were more active than the template molecule. Ligands of the template molecule and new designed compounds were used for molecular docking to study the interaction of these compounds with the protein receptor. The results show that the ligands would form hydrogen-bonding interactions with the residues LEU58, THR83, GLN62, MET155, LYS119 and ALA154 of the protein receptor generally, thereby providing additional insights for the design of even more effective HIV/AIDS drugs.     

Combined 2D and 3D-QSAR,molecular modelling and docking studies of pyrazolodiazepinones as novel phosphodiesterase 2 inhibitors

Selective inhibition of phosphodiesterase 2 (PDE2) in cells where it is located elevates cyclic guanosine monophosphate (cGMP) and acts as novel analgesic with antinociceptive activity. Three-dimensional quantitative structure–activity relationship (QSAR) studies for pyrazolodiazepinone inhibitors exhibiting PDE2 inhibition were performed using comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA) and Topomer CoMFA, and two-dimensional QSAR study was performed using a Hologram QSAR (HQSAR) method. QSAR models were generated using training set of 23 compounds and were validated using test set of nine compounds. The optimum partial least squares (PLS) for CoMFA-Focusing, CoMSIA-SDH, Topomer CoMFA and HQSAR models exhibited good ‘leave-one-out’ cross validated correlation coefficient (q²) of 0.790, 0.769, 0.840 and 0.787, coefficient of determination (r²) of 0.999, 0.964, 0.979 and 0.980, and high predictive power (r²_pred) of 0.796, 0.833, 0.820 and 0.803 respectively. Docking studies revealed that those inhibitors able to bind to amino acid Gln859 by cGMP binding orientation called ‘glutamine-switch’, and also bind to the hydrophobic clamp of PDE2 isoform, could possess high selectivity for PDE2. From the results of all the studies, structure–activity relationships and structural requirements for binding to active site of PDE2 were established which provide useful guidance for the design and future synthesis of potent PDE2 inhibitors.     

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.

     

Three-dimensional Quantitative Structure-Activity Relationship Analyses of a Series of Butenolide ETA Antagonists

Two kinds of Three-dimensional Quantitative Structure-activity Relationship(3D-QSAR) methods,comparative molecular filed analysis(CoMFA) and comparative molecular similarity indices analysis (CoMSIA) ,were applied to analyze the structure-activity relationship of a series of 63 butenolide ETA selective antagonists with respect to their inhibition against human ETA receptor,The CoMFA and CoMSIA models were developed for the conceivable alignment of the molecules based on a template structure from the crystallized data.The statistical results from the initial orientation of the aligned molecules show that the 3D-QSAR model from CoMFA(q^2=0.543) is obviously superior to that from the conventional CoMSIA(q^2=0.407).In order to refine the model,all-space search (ASS) was applied to minimize the field sampling process.By rotating and translating the molecular aggregate within the grid systematically,all the possible samplings of the molecular fields were tested and subsequently the one with the highest q^2 was picked out .The comparison of the sensitivity of CoMFA and CoMSIA to different space orientation shows that the CoMFA q^2 values are more sensitive to the translations and rotations of the aligned molecules with respect to the lattice than those of CoMSIA.The best CoMFA model from ASS was further refined by the region focused technique.The high quality of the best model is indicated by the high corss-validated correlation and the prediction on the external test set.The CoMFA coefficient contour plots identify several key features that explain the wide range of activities,which may help us to design new effective ETA selective antagonists.     

3D－QSAR Study on Diindolylmethane and Its Analogues with Comparative Molecular Field Analysis(CoMFA)

Comparative molecular field analysis (CoMFA),a three dimensional quantitative structure-activity relationship (3D-QSAR) method was applied to a series of diindolylmethane(DIM) analogs to study the relationship between their structure and their induction of CYP 1A1-associated ethoxyresorufin-O-deethylase(EROD) activity.A DISCO model of pharmacophore was derved to guide the superposition of the compounds.The coefficient of cross-validation (q^2) and non cross-validation(r^2) for the model established by the study are 0.827 and 0.988 respectively,the value of variance ratio (F) is 103.53 and standard error estimate (SEE)is 0.044.These values indicate that the CoMFA model derived is significant and might have a good prediction for the catalytic activity of DIM compounds.As a consequence,the predicted activity values of new designed compounds were all higher than that of the reported value.     

R基团搜索技术用于硼酸三肽蛋白酶体抑制剂的分子设计

本文采用Topomer CoMFA对44个Tyropeptin硼酸三肽类蛋白酶体抑制剂进行三维定量构效关系(3D-QSAR)分析。所得最优模型的拟合、交互验证、及外部验证的复相关系数分别为0.983、0.651、0.963。采用Topomer search对ZINC数据库进行R基团的虚拟筛选,得到具有特定活性贡献的R基团,以活性最高的分子为模板过滤,得到7个R1和5个R2基团。并以此设计得到活性优于模板分子的20个新化合物。结果表明,所建立的Topomer CoMFA模型具有良好的稳定性和预测能力,基于R集团的Topomer search技术可以有效筛选,并为设计出新的蛋白酶体抑制剂提供理论依据。     

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.     

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.     

Insight into the structural requirements of benzimidazole derivatives as interleukin‐2 inducible t‐cell kinase inhibitors by computational explorations

In the present work, a set of ligand‐ and receptor‐based 3D‐QSAR models were developed to explore the structure–activity relationship of 109 benzimidazole‐based interleukin‐2‐inducible T‐cell kinase (ITK) inhibitors. In order to reveal the requisite 3D structural features impacting the biological activities, a variety of in silico modeling approaches including the comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), docking, and molecular dynamics were applied. The results showed that the ligand‐based CoMFA model (Q² = 0.552, R²_ncv = 0.908, R²_pred = 0.787, SEE = 0.252, SEP = 0.558) and CoMSIA model (Q² = 0.579, R²_ncv = 0.914, R²_pred = 0.893, SEE = 0.240, SEP = 0.538) were superior to other models with greater predictive power. In addition, a combined analysis between the 3D contour maps and docking results showed that: (1) Compounds with bulky or hydrophobic substituents near ring D and electropositive or hydrogen acceptor groups around rings C and D could increase the activity. (2) The key amino acids impacting the receptor–ligand interactions in the binding pocket are Met438, Asp500, Lys391, and Glu439. The results obtained from this work may provide helpful guidelines in design of novel benzimidazole analogs as inhibitors of ITK. © 2013 Wiley Periodicals, Inc.     

Binding conformations,QSAR, and molecular design of Alkene‐3‐quinolinecarbonitriles as Src inhibitors

A theoretical study on binding orientations and quantitative structure–activity relationship (QSAR) of a novel series of alkene‐3‐quinolinecarbonitriles acting as Src inhibitors has been carried out by using the docking study and three‐dimensional QSAR (3D‐QSAR) analyses. The appropriate binding orientations and conformations of these compounds interacting with Src kinase were revealed by the docking studies, and the established 3D‐QSAR models show significant statistical quality and satisfactory predictive ability, with high R² values and q² values: comparative molecular field analysis (CoMFA) model (q² = 0.748, R² = 0.972), comparative molecular similarity indices analysis (CoMSIA) model (q² = 0.731, R² = 0.987). The systemic external validation indicated that both CoMFA and CoMSIA models possessed high predictive powers with $ R{^2}_{\!\!\!\rm pred} $ values of 0.818 and 0.892, $ {r^2}_{\!\!\!\rm m} $ values of 0.879 and 0.886, $ {r^2}_{\!\!\!\rm m(LOO)} $ values of 0.874 and 0.874, $ r^2_{\rm m(overall)} $ values of 0.879 and 0.885, respectively. Several key structural features of the compounds responsible for inhibitory activity were discussed in detail. Based on these structural factors, eight new compounds with quite higher predicted Src‐inhibitory activities have been designed and presented. We hope these theoretical results can offer some valuable references for the pharmaceutical molecular design as well as the action mechanism analysis. © 2012 Wiley Periodicals, Inc.     

类黄酮抑制P糖蛋白的三维定量构效关系与作用模式

采用基于R基团搜索技术的Topomer CoMFA建立了30个类黄酮类P糖蛋白抑制剂的三维定量构效关系(3D-QSAR)模型, 并用包括9个样本的测试集验证模型的外部预测能力. 所得模型的拟合、 交互验证以及外部验证的复相关系数分别为r²=0.971, q²=0.728和\begin{array}{c}{r}_{\mathit{pred}}^2\end{array}=0.816. 在此基础上, 运用Surflex-dock分子对接法研究了白杨素及其异戊烯化衍生物与P糖蛋白的作用模式. 结果表明, 异戊烯化修饰可显著提高类黄酮的亲脂性, 修饰产物能更好地与P糖蛋白的疏水性口袋契合, 二者结合程度高.     

Docking Study and Three‐Dimensional Quantitative Structure‐Activity Relationship (3D‐QSAR) Analyses and Novel Molecular Design of a Series of 4‐Aminoquinazolines as Inhibitors of Aurora B Kinase

The Aurora proteins are critical regulators of major mitotic events and attractive targets for anticancer therapy. 3D‐QSAR studies based on molecular docking were performed on a dataset of 40 4‐aminoquinazolines compounds. The CoMSIA model produced significantly better results than CoMFA model, with q²=0.652 and r²=0.991. The contours analysis provides useful information about the structural requirements for 4‐aminoquinazolines for inhibiting Aurora B. Scaffold hopping method was used to generate new structures based on the maximum common substructure of the training and test set compounds. The ADMET property, binding affinity and inhibitory activity of the new designed compounds were predicted, respectively. Finally 16 compounds were identified as the novel inhibitors for Aurora B kinase.     

A selectivity study of benzenesulfonamide derivatives on human carbonic anhydrase II/IX by 3D-QSAR,Molecular Docking and Molecular Dynamics Simulation

Nowadays, different approaches have been pursued with the intent to develop sulfonamide-like carbonic anhydrase inhibitors that possess better selectivity profiles toward the different human isoforms of the enzyme. Here, we used conventional 3D-QSAR methods, including comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), and Topomer CoMFA, to construct three-dimensional quantitative structure-activity relationship (3D-QSAR) models for benzenesulfonamide derivatives as human carbonic anhydrase (hCA) II/IX inhibitors. The theoretical models had good reliability (R²>0.75) and predictability (Q²>0.55), and the contour maps could graphically present the contributions of the force fields for activity and identify the structural divergence between human carbonic anhydrase II inhibitors and human carbonic anhydrase IX inhibitors. Consequently, we explored the selectivity of inhibitor for human carbonic anhydrase II and IX through molecular docking, and the difference of activity coincides with the potential binding mode well. According to the results of the predicted values and the molecule docking, we found that the inhibitors published in the literature had stronger inhibition on the hCA IX; based on the theoretical models, we designed seven new compounds with good potential activity and reasonably good ADMET profile, which could selectively inhibit hCA IX. Molecular Dynamics Simulation showed that newly-designed compound D7 had good selectivity on hCA IX. The findings from 3D-QSAR and docking studies maybe helpful in the rational drug design of isoform-selective inhibitors.     

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.

     

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     

Theoretical studies on QSAR and mechanism of 2‐indolinone derivatives as tubulin inhibitors

The theoretical studies on three‐dimensional quantitative structure activity relationship (3D‐QSAR) and action mechanism of a series of 2‐indolinone derivatives as tubulin inhibitors against human breast cancer cell line MDA‐MB‐231 have been carried out. The established 3D‐QSAR model from the comparative molecular field analysis (CoMFA) shows not only significant statistical quality but also predictive ability, with high correlation coefficient (R² = 0.986) and cross‐validation coefficient (q² = 0.683). In particular, the appropriate binding orientations and conformations of these 2‐indolinone derivatives interacting with tubulin are located by docking study, and it is very interesting to find that the plot of the energy scores of these compounds in DOCK versus the corresponding experimental pIC₅₀ values exhibits a considerable linear correlation. Therefore, the inhibition mechanism that 2‐indolinone derivatives were regarded as tubulin inhibitors can be theoretically confirmed. Based on such an inhibition mechanism along with 3D‐QSAR results, some important factors improving the activities of these compounds were discussed in detail. These factors can be summarized as follows: the H atom adopted as substituent R₁, the substituent R₂ with higher electropositivity and smaller bulk, the substituents R₄–R₆ (on the phenyl ring) with higher electropositivity and larger bulk, and so on. These results can offer useful theoretical references for understanding the action mechanism, designing more potent inhibitors, and predicting their activities prior to synthesis. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

磺酰脲类乙酰羟基酸合成酶抑制剂Topomer CoMFA研究

磺酰脲类除草剂是一类高选择性、广谱、低毒的化合物,在世界范围内得到了广泛的应用。本文采用易位体-比较分子力场法（Topomer CoMFA）对75个磺酰脲类化合物与植物源野生型拟南芥AHAS酶的离体相互作用进行了三维定量构效关系研究,快速准确地构建了Topomer CoMFA模型,该模型具有较强的预测能力（交叉验证相关系数q2为0.890,非交叉验证相关系数r2为0.967）。此模型对测试集的10个化合物的pKi值进行预测,其预测值与实际值一致。