A Combination of 3D-QSAR, Molecular Docking and Molecular Dynamics Simulation Studies of Benzimidazole-Quinolinone Derivatives as iNOS Inhibitors

Inducible Nitric Oxide Synthase (iNOS) has been involved in a variety of diseases, and thus it is interesting to discover and optimize new iNOS inhibitors. In previous studies, a series of benzimidazole-quinolinone derivatives with high inhibitory activity against human iNOS were discovered. In this work, three-dimensional quantitative structure-activity relationships (3D-QSAR), molecular docking and molecular dynamics (MD) simulation approaches were applied to investigate the functionalities of active molecular interaction between these active ligands and iNOS. A QSAR model with R² of 0.9356, Q² of 0.8373 and Pearson-R value of 0.9406 was constructed, which presents a good predictive ability in both internal and external validation. Furthermore, a combined analysis incorporating the obtained model and the MD results indicates: (1) compounds with the proper-size hydrophobic substituents at position 3 in ring-C (R₃ substituent), hydrophilic substituents near the X₆ of ring-D and hydrophilic or H-bond acceptor groups at position 2 in ring-B show enhanced biological activities; (2) Met368, Trp366, Gly365, Tyr367, Phe363, Pro344, Gln257, Val346, Asn364, Met349, Thr370, Glu371 and Tyr485 are key amino acids in the active pocket, and activities of iNOS inhibitors are consistent with their capability to alter the position of these important residues, especially Glu371 and Thr370. The results provide a set of useful guidelines for the rational design of novel iNOS inhibitors.     

国内有关黄酮类化合物及其衍生物的2D-QSAR和3D-QSAR研究进展

天然黄酮类化合物在自然界中分布广泛,因其具有抗癌、抗氧化、抗病毒等等众多生物活性,一直以来都是人们关注和研究的热点.但由于其具有溶解性差、生物利用度低等缺点,对天然黄酮类化合物的结构修饰的研究十分必要.利用计算机药物辅助软件对黄酮类化合物进行2D-QSAR和3D-QSAR理论研究,为黄酮类化合物的分子结构修饰提供重要的理论指导.近来,2D-QSAR和3D-QSAR的研究方法已成为实验合成研究的重要依据和常规方法.     

靛玉红类衍生物的三维定量构效关系研究

采用比较分子场方法(CoMFA)对16个靛玉红类衍生物进行三维定量构效关系(3D-QSAR)研究,探讨靛玉红类衍生物的抗癌活性与其结构的关系。建立了3D-QSAR的CoMFA模型,得到了非交叉验证相关系数r2=0.987,标准偏差SD=2.141,统计方差比F=117.136,交叉验证相关系数q2=0.496。该模型合理、可信,具有一定的预测能力,表明改变取代基R1基团的体积才是提高靛玉红类衍生物抗癌活性的主要途径,为开发活性更好的抗癌药物提供一定的理论参考。     

Docking-Based 3D-QSAR Studies for 1,3,4-oxadiazol-2-one Derivatives as FAAH Inhibitors

This work aimed to construct 3D-QSAR CoMFA and CoMSIA models for a series of 31 FAAH inhibitors, containing the 1,3,4-oxadiazol-2-one moiety. The obtained models were characterized by good statistical parameters: CoMFA Q² = 0.61, R² = 0.98; CoMSIA Q² = 0.64, R² = 0.93. The CoMFA model field contributions were 54.1% and 45.9% for steric and electrostatic fields, respectively. In the CoMSIA model, electrostatic, steric, hydrogen bond donor, and hydrogen acceptor properties were equal to 34.6%, 23.9%, 23.4%, and 18.0%, respectively. These models were validated by applying the leave-one-out technique, the seven-element test set (CoMFA r²_test-set = 0.91; CoMSIA r²_test-set = 0.91), a progressive scrambling test, and external validation criteria developed by Golbraikh and Tropsha (CoMFA r²₀ = 0.98, k = 0.95; CoMSIA r²₀ = 0.98, k = 0.89). As the statistical significance of the obtained model was confirmed, the results of the CoMFA and CoMSIA field calculation were mapped onto the enzyme binding site. It gave us the opportunity to discuss the structure–activity relationship based on the ligand–enzyme interactions. In particular, examination of the electrostatic properties of the established CoMFA model revealed fields that correspond to the regions where electropositive substituents are not desired, e.g., in the neighborhood of the 1,3,4-oxadiazol-2-one moiety. This highlights the importance of heterocycle, a highly electronegative moiety in this area of each ligand. Examination of hydrogen bond donor and acceptor properties contour maps revealed several spots where the implementation of another hydrogen-bond-donating moiety will positively impact molecules’ binding affinity, e.g., in the neighborhood of the 1,3,4-oxadiazol-2-one ring. On the other hand, there is a large isopleth that refers to the favorable H-bond properties close to the terminal phenoxy group of a ligand, which means that, generally speaking, H-bond acceptors are desired in this area.     

3D-QSAR and Molecular Docking Studies on Derivatives of MK-0457, GSK1070916 and SNS-314 as Inhibitors against Aurora B Kinase

Development of anticancer drugs targeting Aurora B, an important member of the serine/threonine kinases family, has been extensively focused on in recent years. In this work, by applying an integrated computational method, including comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), homology modeling and molecular docking, we investigated the structural determinants of Aurora B inhibitors based on three different series of derivatives of 108 molecules. The resultant optimum 3D-QSAR models exhibited (q(2) = 0.605, r(2) (pred) = 0.826), (q(2) = 0.52, r(2) (pred) = 0.798) and (q(2) = 0.582, r(2) (pred) = 0.971) for MK-0457, GSK1070916 and SNS-314 classes, respectively, and the 3D contour maps generated from these models were analyzed individually. The contour map analysis for the MK-0457 model revealed the relative importance of steric and electrostatic effects for Aurora B inhibition, whereas, the electronegative groups with hydrogen bond donating capacity showed a great impact on the inhibitory activity for the derivatives of GSK1070916. Additionally, the predictive model of the SNS-314 class revealed the great importance of hydrophobic favorable contour, since hydrophobic favorable substituents added to this region bind to a deep and narrow hydrophobic pocket composed of residues that are hydrophobic in nature and thus enhanced the inhibitory activity. Moreover, based on the docking study, a further comparison of the binding modes was accomplished to identify a set of critical residues that play a key role in stabilizing the drug-target interactions. Overall, the high level of consistency between the 3D contour maps and the topographical features of binding sites led to our identification of several key structural requirements for more potency inhibitors. Taken together, the results will serve as a basis for future drug development of inhibitors against Aurora B kinase for various tumors.     

氨基吡啶衍生物的合成研究进展

综述了氨基吡啶的合成研究进展,介绍了氨基吡啶的四种制备方法:金属还原法、催化加氢还原法、电化学还原法、微生物法及氨基吡啶的代表物质之一2-氯-3-氨基吡啶的相关合成方法及新进展。     

Design,Synthesis, Biological Evaluation, 2D-QSAR Modeling,and Molecular Docking Studies of Novel 1H-3-Indolyl Derivatives as Significant Antioxidants

Novel candidates of 3-(4-(thiophen-2-yl)-pyridin/pyran/pyrimidin/pyrazol-2-yl)-1H-indole derivatives (2–12) were designed by pairing the pyridine/pyrane/pyrimidine/pyrazole heterocycles with indole and thiophene to investigate their potential activities as (2,2′-azinobis (3-ethylbenzothiazoline-6-sulfonic acid) inhibitors. The purpose of these derivatives’ modification is to create high-efficiency antioxidants, especially against ABTS, as a result of the efficiency of this set of key heterocycles in the inhibition of ROS. Herein, 2D QSAR modeling was performed to recommend the most promising members for further in vitro investigations. Furthermore, the pharmacological assay for antioxidant activity evaluation of the yielded indole-based heterocycles was tested against ABTS (2,2′-azinobis (3-ethylbenzothiazoline-6-sulfonic acid); by utilizing ascorbic acid as the standard. Candidate 10 showed higher antioxidant activity (IC₅₀ = 28.23 μg/mL) than ascorbic acid itself which achieved (IC₅₀ = 30.03 μg/mL). Moreover, molecular docking studies were performed for the newly designed and synthesized drug candidates to propose their mechanism of action as promising cytochrome c peroxidase inhibitors compared to ascorbic acid as a reference standard. Our findings could be promising in the medicinal chemistry scope for further optimization of the newly designed and synthesized compounds regarding the introduced structure-activity relationship study (SAR) in order to get a superior antioxidant lead compound in the near future.     

Studies of New Fused Benzazepine as Selective Dopamine D3 Receptor Antagonists Using 3D-QSAR, Molecular Docking and Molecular Dynamics

In recent years, great interest has been paid to the development of compounds with high selectivity for central dopamine (DA) D3 receptors, an interesting therapeutic target in the treatment of different neurological disorders. In the present work, based on a dataset of 110 collected benzazepine (BAZ) DA D3 antagonists with diverse kinds of structures, a variety of in silico modeling approaches, including comparative molecular field analysis (CoMFA), comparative similarity indices analysis (CoMSIA), homology modeling, molecular docking and molecular dynamics (MD) were carried out to reveal the requisite 3D structural features for activity. Our results show that both the receptor-based (Q(2) = 0.603, R(2) (ncv) = 0.829, R(2) (pre) = 0.690, SEE = 0.316, SEP = 0.406) and ligand-based 3D-QSAR models (Q(2) = 0.506, R(2) (ncv) =0.838, R(2) (pre) = 0.794, SEE = 0.316, SEP = 0.296) are reliable with proper predictive capacity. In addition, a combined analysis between the CoMFA, CoMSIA contour maps and MD results with a homology DA receptor model shows that: (1) ring-A, position-2 and R(3) substituent in ring-D are crucial in the design of antagonists with higher activity; (2) more bulky R(1) substituents (at position-2 of ring-A) of antagonists may well fit in the binding pocket; (3) hydrophobicity represented by MlogP is important for building satisfactory QSAR models; (4) key amino acids of the binding pocket are CYS101, ILE105, LEU106, VAL151, PHE175, PHE184, PRO254 and ALA251. To our best knowledge, this work is the first report on 3D-QSAR modeling of the new fused BAZs as DA D3 antagonists. These results might provide information for a better understanding of the mechanism of antagonism and thus be helpful in designing new potent DA D3 antagonists.     

Generation of Non-Nucleotide CD73 Inhibitors Using a Molecular Docking and 3D-QSAR Approach

Radiotherapy and chemotherapy are conventional cancer treatments. Around 60% of all patients who are diagnosed with cancer receive radio- or chemotherapy in combination with surgery during their disease. Only a few patients respond to the blockage of immune checkpoints alone, or in combination therapy, because their tumours might not be immunogenic. Under these circumstances, an increasing level of extracellular adenosine via the activation of ecto-5’-nucleotidase (CD73) and consequent adenosine receptor signalling is a typical mechanism that tumours use to evade immune surveillance. CD73 is responsible for the conversion of adenosine monophosphate to adenosine. CD73 is overexpressed in various tumour types. Hence, targetting CD73’s signalling is important for the reversal of adenosine-facilitated immune suppression. In this study, we selected a potent series of the non-nucleotide small molecule inhibitors of CD73. Molecular docking studies were performed in order to examine the binding mode of the inhibitors inside the active site of CD73 and 3D-QSAR was used to study the structure–activity relationship. The obtained CoMFA (q² = 0.844, ONC = 5, r² = 0.947) and CoMSIA (q² = 0.804, ONC = 4, r² = 0.954) models showed reasonable statistical values. The 3D-QSAR contour map analysis revealed useful structural characteristics that were needed to modify non-nucleotide small molecule inhibitors. We used the structural information from the overall docking and 3D-QSAR results to design new, potent CD73 non-nucleotide inhibitors. The newly designed CD73 inhibitors exhibited higher activity (predicted pIC₅₀) than the most active compound of all of the derivatives that were selected for this study. Further experimental studies are needed in order to validate the new CD73 inhibitors.     

3-氨基-1,2-苯并-α-吡喃酮类衍生物的合成及抗氧化研究

以2,3,4-三甲氧基苯甲醛为原料,经4步反应合成了两个未见文献报道的1,2-苯并-α-吡喃酮类衍生物,其结构经NMR和MS确证。体外抗氧化性能测试结果显示,3-苯乙酰氨基-7,8-二甲氧基-1,2-苯并-α-吡喃酮在高浓度时显示出比对照品和3-氨基-7,8-二甲氧基-1,2-苯并-α-吡喃酮更好的DPPH自由基清除活性,在低浓度时则对羟自由基表现出比对照品及3-氨基-7,8-二甲氧基-1,2-苯并-α-吡喃酮更好的抑制活性,在ABTS自由基清除实验中,低浓度的3-氨基-7,8-二甲氧基-1,2-苯并-α-吡喃酮展现出较对照品更佳的清除活性,两者具有进一步研究的价值。     

Understanding the Molecular Basis of 5-HT4 Receptor Partial Agonists through 3D-QSAR Studies

Alzheimer’s disease (AD) is a neurodegenerative disorder whose prevalence has an incidence in senior citizens. Unfortunately, current pharmacotherapy only offers symptom relief for patients with side effects such as bradycardia, nausea, and vomiting. Therefore, there is a present need to provide other therapeutic alternatives for treatments for these disorders. The 5-HT₄ receptor is an attractive therapeutic target since it has a potential role in central and peripheral nervous system disorders such as AD, irritable bowel syndrome, and gastroparesis. Quantitative structure-activity relationship analysis of a series of 62 active compounds in the 5-HT₄ receptor was carried out in the present work. The structure-activity relationship was estimated using three-dimensional quantitative structure-activity relationship (3D-QSAR) techniques based on these structures’ field molecular (force and Gaussian field). The best force-field QSAR models achieve a value for the coefficient of determination of the training set of R²_training = 0.821, and for the test set R²_test = 0.667, while for Gaussian-field QSAR the training and the test were R²_training = 0.898 and R²_test = 0.695, respectively. The obtained results were validated using a coefficient of correlation of the leave-one-out cross-validation of Q²_LOO = 0.804 and Q²_LOO = 0.886 for force- and Gaussian-field QSAR, respectively. Based on these results, novel 5-HT₄ partial agonists with potential biological activity (pEC₅₀ 8.209–9.417 for force-field QSAR and 9.111–9.856 for Gaussian-field QSAR) were designed. In addition, for the new analogues, their absorption, distribution, metabolism, excretion, and toxicity properties were also analyzed. The results show that these new derivatives also have reasonable pharmacokinetics and drug-like properties. Our findings suggest novel routes for the design and development of new 5-HT₄ partial agonists.     

Molecular Modeling Studies of 4,5-Dihydro-1H-pyrazolo[4,3-h] quinazoline Derivatives as Potent CDK2/Cyclin A Inhibitors Using 3D-QSAR and Docking

CDK2/cyclin A has appeared as an attractive drug targets over the years with diverse therapeutic potentials. A computational strategy based on comparative molecular fields analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) followed by molecular docking studies were performed on a series of 4,5-dihydro-1H-pyrazolo[4,3-h]quinazoline derivatives as potent CDK2/cyclin A inhibitors. The CoMFA and CoMSIA models, using 38 molecules in the training set, gave r(2) (cv) values of 0.747 and 0.518 and r(2) values of 0.970 and 0.934, respectively. 3D contour maps generated by the CoMFA and CoMSIA models were used to identify the key structural requirements responsible for the biological activity. Molecular docking was applied to explore the binding mode between the ligands and the receptor. The information obtained from molecular modeling studies may be helpful to design novel inhibitors of CDK2/cyclin A with desired activity.     

大豆多酚衍生物的抗氧化活性研究

染料木素与苯磺酸发生酯化反应得到染料木素-7-苯磺酸酯.利用DPPH自由基、超氧阴离子自由基、羟自由基、总还原力以及总抗氧化活性测定大豆多酚染料木素(GE)及其衍生物染料木素-7-苯磺酸酯(GB)的抗氧化活性,结果表明GE和GB均具有一定抗氧化活性,且呈量效关系,清除羟基自由基的能力为:GB>GE>Vc.     

Pharmacophore and Molecular Docking Guided 3D-QSAR Study of Bacterial Enoyl-ACP Reductase (FabI) Inhibitors

Enoyl acyl carrier protein (ACP) reductase (FabI) is a potential target for the development of antibacterial agents. Three-dimensional quantitative structure-activity relationships (3D-QSAR) for substituted formamides series of FabI inhibitors were investigated using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) techniques. Pharmacophore and molecular docking methods were used for construction of the molecular alignments. A training set of 36 compounds was performed to create the 3D-QSAR models and their external predictivity was proven using a test set of 11 compounds. Graphical interpretation of the results revealed important structural features of the formamides related to the active site of FabI. The results may be exploited for further optimization of the design of new potent FabI inhibitors.     

Combinatorial Pharmacophore-Based 3D-QSAR Analysis and Virtual Screening of FGFR1 Inhibitors

The fibroblast growth factor/fibroblast growth factor receptor (FGF/FGFR) signaling pathway plays crucial roles in cell proliferation, angiogenesis, migration, and survival. Aberration in FGFRs correlates with several malignancies and disorders. FGFRs have proved to be attractive targets for therapeutic intervention in cancer, and it is of high interest to find FGFR inhibitors with novel scaffolds. In this study, a combinatorial three-dimensional quantitative structure-activity relationship (3D-QSAR) model was developed based on previously reported FGFR1 inhibitors with diverse structural skeletons. This model was evaluated for its prediction performance on a diverse test set containing 232 FGFR inhibitors, and it yielded a SD value of 0.75 pIC₅₀ units from measured inhibition affinities and a Pearson’s correlation coefficient R² of 0.53. This result suggests that the combinatorial 3D-QSAR model could be used to search for new FGFR1 hit structures and predict their potential activity. To further evaluate the performance of the model, a decoy set validation was used to measure the efficiency of the model by calculating EF (enrichment factor). Based on the combinatorial pharmacophore model, a virtual screening against SPECS database was performed. Nineteen novel active compounds were successfully identified, which provide new chemical starting points for further structural optimization of FGFR1 inhibitors.     

油藏微生物群落结构的分子分析

利用PCR-变性梯度凝胶电泳(PCR-DGGE)技术考察本源微生物驱油前后,即营养剂注入前后微生物群落结构的变化.实验室模拟地层培养结果表明,在以N-1(主要成分为淀粉纤维素)为唯一营养剂、培养温度为48℃时,在有氧或厌氧务件下,24 h后原始水样中对驱油有利的本源茵便可以被激活,并随着培养时间的延长最终达到稳定.主要激活的对驱油有利的3种菌株为假单胞菌(Pseudomonas sp.)、梭状芽孢杆菌(Clostridium sp.)和烃类降解菌(Hydro-carbon seep bacterium).经分析,原始底层注入水中虽然含有对驱油不利的脱硫肠状菌属(Desul fotomaftculum salinum),但是在经N-1激活后的水样中均未发现该菌的存在,证明N-1不会激活原始水样中的脱硫肠状菌属.为营养体系的筛选、注入方案的优化设计和驱油效果评价提供了理论支撑.     

对苯二胺类防老剂的色谱分析和分子构型

从胺类分子的立体结构中对映异构体之间的相互转化,使分子极性发生改变,来解释对苯二胺类防老剂在色谱分析时产生双头峰、前伸峰、拖尾峰等异常色谱峰的现象,以便于保证对该类产品纯度的正确分析,以协助生产和市场的正确发展。     

高区域选择性合成3-（芳胺）-N,N-二甲基丙酰胺衍生物

以N,N-二甲基丙烯酰胺和芳胺类化合物为原料,盐酸为促进体系,经过分子间的高区域选择性加成,成功地构建出了3-（芳胺）-N,N-二甲基丙酰胺衍生物。采用1H-NMR、13C-NMR和熔点对目标产物进行了表征。研究表明该方法具有原料便宜易得、原子经济、后处理操作简单和产品纯度高等优点,为3-（芳胺）-N,N-二甲基丙酰胺衍生物的高效、快速合成提供了一条可选的策略。     

3D-QSAR Studies on Barbituric Acid Derivatives as Urease Inhibitors and the Effect of Charges on the Quality of a Model

Urease enzyme (EC 3.5.1.5) has been determined as a virulence factor in pathogenic microorganisms that are accountable for the development of different diseases in humans and animals. In continuance of our earlier study on the helicobacter pylori urease inhibition by barbituric acid derivatives, 3D-QSAR (three dimensional quantitative structural activity relationship) advance studies were performed by Comparative Molecular Field Analysis (CoMFA) and Comparative Molecular Similarity Indices Analysis (CoMSIA) methods. Different partial charges were calculated to examine their consequences on the predictive ability of the developed models. The finest developed model for CoMFA and CoMSIA were achieved by using MMFF94 charges. The developed CoMFA model gives significant results with cross-validation (q²) value of 0.597 and correlation coefficients (r²) of 0.897. Moreover, five different fields i.e., steric, electrostatic, and hydrophobic, H-bond acceptor and H-bond donors were used to produce a CoMSIA model, with q² and r² of 0.602 and 0.98, respectively. The generated models were further validated by using an external test set. Both models display good predictive power with r²_pred ≥ 0.8. The analysis of obtained CoMFA and CoMSIA contour maps provided detailed insight for the promising modification of the barbituric acid derivatives with an enhanced biological activity.     

Structural Determination of Three Different Series of Compounds as Hsp90 Inhibitors Using 3D-QSAR Modeling, Molecular Docking and Molecular Dynamics Methods

Hsp90 is involved in correcting, folding, maturation and activation of a diverse array of client proteins; it has also been implicated in the treatment of cancer in recent years. In this work, comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), molecular docking and molecular dynamics were performed on three different series of Hsp90 inhibitors to build 3D-QSAR models, which were based on the ligand-based or receptor-based methods. The optimum 3D-QSAR models exhibited reasonable statistical characteristics with averaging internal q(2) > 0.60 and external r(2) (pred) > 0.66 for Benzamide tetrahydro-4H-carbazol-4-one analogs (BT), AT13387 derivatives (AT) and Dihydroxylphenyl amides (DA). The results revealed that steric effects contributed the most to the BT model, whereas H-bonding was more important to AT, and electrostatic, hydrophobic, H-bond donor almost contributed equally to the DA model. The docking analysis showed that Asp93, Tyr139 and Thr184 in Hsp90 are important for the three series of inhibitors. Molecular dynamics simulation (MD) further indicated that the conformation derived from docking is basically consistent with the average structure extracted from MD simulation. These results not only lead to a better understanding of interactions between these inhibitors and Hsp90 receptor but also provide useful information for the design of new inhibitors with a specific activity.