含磷嘧啶类CDK9抑制剂的分子对接、3D-QSAR和分子动力学模拟

选取64个具有潜力的含磷嘧啶类细胞周期依赖性蛋白激酶(CDK9)小分子抑制剂,采用分子对接方法研究了该类小分子与CDK9的结合作用,结果表明,分子构象、氢键形成、疏水性和氨基酸残基Cys106在此类抑制剂与CDK9的结合过程中具有重要作用.在配体叠合的基础上,运用比较分子力场分析(Co MFA)、比较分子相似性指数分析(Co MSIA)和Topomer Co MFA(T-COMFA)研究了分子结构与抑制活性的关系,发现由训练集立体场、静电场和疏水场组合的Co MSIA模型为最优模型,其内部交叉验证相关系数(Q2=0.557)、非交叉验证相关系数(R2=0.959)和外部预测相关系数(r2=0.863)具有统计学意义,该模型的三维等值线图直观显示了化合物的活性与其三维结构的关系.根据这些结果设计了10个具有新结构的含磷嘧啶类化合物,分子对接和分子动力学模拟结果表明,新化合物和CDK9的结合模式与原化合物64相同,自由能分析从理论上证明了新化合物64d的CDK9抑制活性优于化合物64,并且显示含磷基团与残基Asp109的静电场能在化合物与CDK9作用过程中有重要作用.     

Integration of genetic virtual screening patterns and latent multivariate modeling techniques for QSAR optimization based on combinations and/or interactions between peptides and proteins

Both the concept and the model of snug quantitative structure-activity relationship (QSAR) were pro-posed and developed for molecular design through constructing QSAR based on some known mode of receptor/ligand interactions. Many disadvantages of traditional models can be avoided by using the proposed method because the traditional models only determined upon molecular structural features in sample sets themselves. A genetic virtual screening of peptide/protein combinations (GVSPPC) is proposed for the first time by utilizing this idea to examine peptide/protein affinity activities. A genetic algorithm (GA) was developed for screening combinative targets with an interaction mode for virtual receptors. GVSPPC succeeds in disposing difficulties in rational QSAR,in order to search for the ligand/receptor interactions on conditions of unknown structures. Some bioactive oligo-/poly-peptide systems covering 58 angiotensin converting enzyme (ACE) inhibitors and 18 double site mutation residues in camel antibody protein cAb-Lys3 were investigated by GVSPPC with satisfactory results (R 2 cu>0.91,Q 2 cv > 0.86,ERMS=0.19-0.95),respectively,which demonstrates that GVSPPC is more inter-pretable in the ligand-receptor interaction than the traditional QSAR method.     

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.     

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.     

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.     

Design, Synthesis and Evaluation of 6-Oxo-1, 6-dihydropyrimidine-2,5-dicarboxamide Derivatives as MMP 13 Inhibitors

We designed and synthsized a series of novel 6-oxo-1,6-dihydropyrimidine-2,5-dicarboxamide derivatives and evaluated their inhibition effects on MMP 3, MMP 12 and MMP 13. The pharmacological results show that they have potent and highly selective activity of inhibiting MMP 13.     

Computational Investigations about the Effects of Hetero‐molecular Aggregation on Bioactivities: a Case of Neonicotinoids and Water

Molecular aggregation state of bioactive compounds plays a key role in bio‐interactive procedure. Diverse aggregation states of bioactive compounds contribute to different biological or chemical properties. Water‐bridge, as the simple hetero‐molecular aggregation, has been found bridging the binding between many bioactive compounds and their targets through hydrogen bonding network, e.g. in the recognition of neonicotinoids with insect nAChRs. To better understanding the roles of water‐bridge on bioactivities of compounds, an approach of hetero‐dimeric aggregation with water was proposed. Quantitative structure‐activity relationship (QSAR) and pharmacophore modeling investigations were applied on 19 neonicotinoids, as well as their aggregates with water. The aggregate‐based CoMSIA, PHASE and linear QSAR models presented better statistical significance and predictabilities than the monomer ones, which indicated that the bioactivities correlated with the aggregate properties and water bridged hydrogen bond of the active site. All results revealed the essential roles of water‐bridge in ligand recognition, which should be considered in future ligand design and optimization.     

Inhibition of Metalloproteinase Activity by Chinese Formulations Used to Treat Inflammatory Diseases

Introduction Matrixmetalloproteinases(MMPs)areaclassof zinc requiringextracellularendopeptidasesthatcanto getherdegradeallcomponentsoftheextracellularma trixandbasementmembranes[1—3].Theyplayimpor tantrolesinconnectivetissueremodeling,occurringin normalb…     

Studies on molecular mechanism between SHP2 and pyridine derivatives by 3D-QSAR,molecular docking and MD simulations

BackgroundSrc homology 2 (SH2)-containing protein tyrosine phosphatase 2 (SHP2) as a major phosphatase would affect the development of tumors by regulating several cellular processes, and is a significant potential target for cancer treatment.MethodsIn the present work, a series of pyridine derivatives possessing a wide range of inhibitory activity was employed to investigate the structural requirements by developing three dimensional quantitative structure–activity relationship (3D-QSAR) models using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) methods. The results show that CoMFA (R²_cv = 0.646, R²_pred = 0.5587) and CoMSIA (R²_cv = 0.777, R²_pred = 0.7131) have excellent stability and predictability. The relationship between the inhibitory activity and structure of the inhibitors was analyzed by the derived contour maps. Furthermore, the QSAR models were validated by molecular docking and molecular dynamics simulations, which were also applied to reveal the potential molecular mechanism of these inhibitors.FindingsIt was found that Arg110, Asn216, Thr218, Thr252 and Pro490 play a crucial role in stabilizing the inhibitors. Additionally, MM/PBSA calculations provided the binding free energy were also conducted to explain the discrepancy of binding activities. Overall, the outcomes of this work could provide useful information and theoretical guidance for the development of novel and potent SHP2 inhibitors.     

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     

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.     

(13)C NMR quantitative spectrometric data-activity relationship (QSDAR) models of steroids binding the aromatase enzyme

Five quantitative spectroscopic data-activity relationships (QSDAR) models for 50 steroidal inhibitors binding to aromatase enzyme have been developed based on simulated (13)C nuclear magnetic resonance (NMR) data. Three of the models were based on comparative spectral analysis (CoSA), and the two other models were based on comparative structurally assigned spectral analysis (CoSASA). A CoSA QSDAR model based on five principal components had an explained variance (r(2)) of 0.78 and a leave-one-out (LOO) cross-validated variance (q(2)) of 0.71. A CoSASA model that used the assigned (13)C NMR chemical shifts from a steroidal backbone at five selected positions gave an r(2) of 0.75 and a q(2) of 0.66. The (13)C NMR chemical shifts from atoms in the steroid template position 9, 6, 3, and 7 each had correlations greater than 0.6 with the relative binding activity to the aromatase enzyme. All five QSDAR models had explained and cross-validated variances that were better than the explained and cross-validated variances from a five structural parameter quantitative structure-activity relationship (QSAR) model of the same compounds. QSAR modeling suffers from errors introduced by the assumptions and approximations used in partial charges, dielectric constants, and the molecular alignment process of one structural conformation. One postulated reason that the variances of QSDAR models are better than the QSAR models is that (13)C NMR spectral data, based on quantum mechanical principles, are more reflective of binding than the QSAR model's calculated electrostatic potentials and molecular alignment process. The QSDAR models provide a rapid, simple way to model the steroid inhibitor activity in relation to the aromatase enzyme.     

酮醇酸还原异构酶(KARI)与其抑制剂间相互作用的分子模拟研究

采用MCCE, Autodock及密度泛函方法对酮醇酸还原异构酶(KARI)与其抑制剂间相互作用进行了研究. 计算结果表明, KARI活性位点中的Mg^2＋在活性位点残基的离子化状态、活性位点的静电性质以及与抑制剂结合等方面起重要的作用; 同时, 抑制剂在结合方式、前线轨道布居及静电势等方面与酶促反应中间体(HOIV)具有一定程度的相似性; 可电离的羧基是当前发现的靶向KARI抑制剂一个重要的结构特征, 进一步推广可认为潜在的抑制剂应该拥有可电离成负电荷的功能团. 在药物设计中考虑到以上结论, 将有利于发现和改造靶向KARI的抑制剂.