γ-分泌酶抑制剂的药效团模型构建

利用Catalyst软件系统, 选择具有较高体外抑制活性的苯并二氮(艹卓)类化合物作为训练集, 经计算机建模, 构象优化, 由Catalyst系统构建出药效团模型. 并结合γ-分泌酶的作用机制等因素, 筛选出一个含有一个芳环中心, 一个疏水中心和两个氢键受体的具有较好预测能力(RMS=0.366343, Correl=0.95535, Weight=1.17389, Config=18.8671)的药效团模型. 该模型的建立有助于设计及合成新型结构的γ-分泌酶抑制剂.     

p53-MDM2结合抑制剂药效团模型的构建

采用Catalyst软件, 选择5类共24个p53-MDM2结合抑制剂作为训练集, 经计算机建模、构象优化, 由Catalyst系统构建出药效团模型, 并对药效团进行有效性分析, 结合已知的p53-MDM2结合抑制剂的结构信息, 筛选得到含有一个芳环中心、三个疏水中心和一个氢键受体的具有较好预测能力(Correl=0.941, Config=17.530, 吟cost=150.830)的药效团模型.     

Pharmacophore modelling and electronic feature analysis of hydroxamic acid derivatives,the HIV integrase inhibitors

Hydroxamic acid derivatives with metal ion binding properties were collected from the literature to generate a pharmacophore and 3D-QSAR model for HIV strand transfer inhibition. The derived pharmacophore model (AAAHRR) recognizes both metal ion binding site and hydrophobic group. The QSAR model generated using this hypothesis expressed statistical significance (r ² = 0.971 for the training set and q ² = 0.913 for the test set). The ability of this pharmacophore model to retrieve other metal ion binding inhibitors was examined by screening the ChemBank database (ligandinfo) incorporated with 10 known strand transfer inhibitors. The studied favourable and unfavourable contours of chemical features (H-bond donor, acceptor and hydrophobic sites) revealed the role of hydrophobic substitution at the fluorobenzene ring and cyclization of the metal ion binding hydroxamic acid in effective integrase inhibition. Analysis of the frontier orbitals, HOMO and LUMO revealed that the nucleophilic / electrophilic interactions depend on the significant overlapping observed at the azaindole and hydroxamic acid groups. In essence, the generated pharmacophore model is competent enough to disclose the essential site-specific interactions involved in the inhibition of HIV integrase, and hence can be used in virtual screening to identify novel scaffolds as leads with increased anti-viral potency.     

EGFR酪氨酸激酶抑制剂的药效团模型构建

以80个作用方式相同, 分子结构特征不同的表皮生长因子受体酪氨酸激酶(EGFR TK)竞争性抑制剂为训练集, 利用计算机药物辅助软件Catalyst, 构建不同的药效团模型, 并结合酪氨酸激酶的作用位点等因素, 筛选出一个含有两个芳环中心, 一个疏水中心和一个阳离子基团的具有较好预测能力(RMS=0.438, Correl=0.908, Weight=1.52, Config=17.36)的药效团模型, 为设计和合成新型结构的EGFR TK抑制剂提供参考.     

A pharmacophore model for NK2 antagonist comprising compounds from several structurally diverse classes

A neurokinin 2 (NK2) antagonist pharmacophore model has been developed on the basis of five non-peptide antagonists from several structurally diverse classes. To evaluate the pharmacophore model, another 20 antagonists were fitted to the model. By use of exhaustive conformational analysis (MMFFs force field and the GB/SA hydration model) and least-squares molecular superimposition studies, 23 of the 25 antagonists were fitted to the model in a low energy conformation with a low RMS value. The pharmacophore model is described by four pharmacophore elements: Three hydrophobic groups and a hydrogen bond donor represented as a vector. The hydrophobic groups are generally aromatic rings, but this is not a requirement. The antagonists bind in an extended conformation with two aromatic rings in a parallel displaced and tilted conformation. The model was able to explain the enantioselectivity of SR48968 and GR159897.     

Identification of potential influenza virus endonuclease inhibitors through virtual screening based on the 3D-QSAR model

Influenza endonucleases have appeared as an attractive target of antiviral therapy for influenza infection. With the purpose of designing a novel antiviral agent with enhanced biological activities against influenza endonuclease, a three-dimensional quantitative structure-activity relationships (3D-QSAR) model was generated based on 34 influenza endonuclease inhibitors. The comparative molecular similarity index analysis (CoMSIA) with a steric, electrostatic and hydrophobic (SEH) model showed the best correlative and predictive capability (q ² = 0.763, r ² = 0.969 and F = 174.785), which provided a pharmacophore composed of the electronegative moiety as well as the bulky hydrophobic group. The CoMSIA model was used as a pharmacophore query in the UNITY search of the ChemDiv compound library to give virtual active compounds. The 3D-QSAR model was then used to predict the activity of the selected compounds, which identified three compounds as the most likely inhibitor candidates.     

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.     

新型酪氨酸激酶小分子抑制的三维药效团研究

通过CATALYST软件包得到了两类HER2抑制的三维药效团模型。尽管亚苄基丙二腈化合物和3-取代吲哚啉-2-酮系列化合物具有完全不同的骨架结构,但得到的药效团却具有共同的特性,这表明当这两类抑制剂和受体发生相互作用时,采用了相似的结合模式。共同的药效团模型包括一个氢键受体,一个氢键给体,一个脂肪类疏水团以及一个芳香类疏水团。根据药效团模型,我们还进行了三维构效关系的研究,结果表明得到的药效团模型具有很好的预测能力(线性回归系数R≈0.96)。药效团模型对于研究酪氨酸激酶小分子抑制剂的结构与活性关系,以及评估和预测此类未知化合物活性具人重要的意义。     

Toward the identification of a reliable 3D QSAR pharmacophore model for the CCK2 receptor antagonism

The present study describes application of computational approaches to identify a validated and reliable 3D QSAR pharmacophore model for the CCK-2R antagonism through integrated ligand and structure based studies using anthranilic sulfonamide and 1,3,4-benzotriazepine based CCK-2R antagonists. The best hypothesis consisted five features viz. two aliphatic hydrophobic, one aromatic hydrophobic, one H-bond acceptor, and one ring aromatic feature with an excellent correlation for 34 training set (r2(training) = 0.83) and 58 test set compounds (r2(test) = 0.74). This model was validated through F-test and docking studies at the active site of the plausible CCK-2R where the 99% significance and well corroboration with the pharmacophore model respectively describes the model's reliability. The model also predicts well to other known clinically effective CCK-2R antagonists. Therefore, the developed model may useful in finding new scaffolds that may aid in design and develop new chemical entities (NCEs) as potent CCK-2R antagonists before their synthesis.     

人类免疫缺陷病毒整合酶二酮酸类抑制剂的三维药效团构建

应用遗传算法相似性程序(GASP), 以作用于I型人类免疫缺陷病毒(human immun-odeficiency virus type 1, HIV-1)整合酶(IN)的二酮酸类(diketoacids, DKAs)抑制剂构建药效团模型. 所选训练集分子均具有可靠的类药性特征及DKAs药效团特征. 尝试将抑制剂与药效团叠合后的构象和抑制剂与IN的对接构象进行叠合, 得到药效团模型与分子对接构象中IN残基的相对位置, 并基于抑制剂的药效团模型特征与周围IN氨基酸残基位置的匹配情况进行药效团特征的修改. 所得最优药效团由1个疏水特征、3对氢键特征和1个氢键供体特征组成. 该药效团的命中物质量(goodness of hit, GH)为0.56, 产出率(Y)达63.6%, 假阳性率(FP)为0.41%. 该药效团具有较好的置信度, 产出率较高而假阳性率较低, 可用于数据库搜索发现新的具有DKAs药效团特征的活性化合物, 也可为先导化合物的改造提供帮助.     

药效团检索设计新的HIV-1蛋白酶抑制剂

通过对自建的未开发化合物三维结构库进行药效团检索,得到了4个对HIV-1蛋白酶抑制活化的化合物,通过构象分析发现包含药效团的构象处于优势构象,而且4个结构都含有带两个邻位羟基的苯环和一个间位羰基的药效团以及公共子结构。通过计算发现它们的疏水参数都很小。在考虑满足包含药效团的结构特征和有适中的疏水参数两个因素的前提下,设计出了新的具有潜在抑制HIV-1蛋白酶活性的化合物。它们的结构都比检索得到的四个化合物更为简单,因此易于合成。     

Pharmacophore-based structure optimization of angiotensin converting enzyme inhibitory peptide

Chemical feature based pharmacophore models were generated for an angiotensin converting enzyme(ACE) inhibitory peptide using the Discovery Studio 2.0 pharmacophore modeling approach. The pharmacophore hypothesis selected has five features(one negative ionizable region,one hydrogen bond donor,one hydrogen bond acceptor and two hydrophobic functional groups). Additionally,ACE inhibitory hexapeptide previously obtained from silkworm pupae protein was optimized to target the ACE based on the selected pharmacophore. The results suggest that tri-peptide(thr-val-phe) may be structural determinant of ACE activity. Docking studies further provided confidence for the validity of the selected pharmacophore model to perform structure optimization of the ACE inhibitory peptide.     

Structure-based virtual screening of dipeptidyl peptidase 4 inhibitors and their in vitro analysis

Type 2 diabetes mellitus (T2DM) is one of the most widely prevalent metabolic disorders with no cure to date thus remains the most challenging task in the current drug discovery. Therefore, the only strategy to control diabetes prevalence is to develop novel efficacious therapeutics. Dipeptidyl Peptidase 4 (DPP-4) inhibitors are currently used as anti-diabetic drugs for the inhibition of incretins. This study aims to construct the chemical feature based on pharmacophore models for dipeptidyl peptidase IV. The structure-based pharmacophore modeling has been employed to evaluate new inhibitors of DPP-4. A four-featured pharmacophore model was developed from crystal structure of DPP-4 enzyme with 4-(2-aminoethyl) benzenesulfonyl fluoride in its active site via pharmacophore constructing tool of Molecular Operating Environment (MOE) consisting F1 Hyd (hydrophobic region), F2 Hyd|Cat|Don (hydrophobic cationic and donor region), F3 Acc (acceptor region) and F4 Hyd (hydrophobic region). The generated pharmacophore model was used for virtual screening of in-house compound library (the available compounds which were used for initial screening to get the few compounds for the current studies). The resultant selected compounds, after virtual screening were further validated using in vitro assay. Furthermore, structure-activity relationship was carried out for the compounds possessing significant inhibition potential after docking studies. The binding free energy of analogs was evaluated via molecular mechanics generalized Born surface area (MM-GBSA) and Poisson-Boltzmann surface area (MM-PBSA) methods using AMBER 16 as a molecular dynamics (MD) simulation package. Based on potential findings, we report that selected candidates are more likely to be used as DPP-4 inhibitors or as starting leads for the development of novel and potent DPP-4 inhibitors.     

基于药效团模型的DHODH抑制剂构效关系研究

利用药效团模型研究二氢乳清酸脱氢酶(Dihydroorotate dehydrogenase,DHODH)抑制剂的构效关系,为DHODH抑制剂的虚拟筛选提供新的方法.以31个具有DHODH抑制活性的化合物为训练集化合物,半数抑制浓度(IC50)范围为7~63000 nmol/L,利用Catalyst/HypoGen算法构建DHODH抑制剂药效团模型,通过对训练集化合物多个构象进行叠合,提取药效团特征及三维空间限制构建药效团模型.利用基于CatScramble的交叉验证方法及评价模型对已知活性化合物的活性预测能力,确定较优药效团模型.模型包含1个氢键受体、3个疏水中心,表征了受体配体相互作用时可能发生的氢键相互作用、疏水相互作用和π-π相互作用,4个药效特征在三维空间的排列概括了DHODH抑制剂产生活性的结构特点.所得较优模型对训练集化合物及测试集化合物的计算活性值与实验活性值的相关系数分别为0.8405和0.8788.利用药效团模型对来源于微生物的系列化合物进行虚拟筛选,筛选出59个预测活性较好的化合物,可作为进一步药物研发的候选化合物.     

Pharmacophore and receptor models for neurokinin receptors

Three neurokinin (NK) antagonist pharmacophore models (Models 1-3) accounting for hydrogen bonding groups in the 'head' and 'tail' of NK receptor ligands have been developed by use of a new procedure for treatment of hydrogen bonds during superimposition. Instead of modelling the hydrogen bond acceptor vector in the strict direction of the lone pair, an angle is allowed between the hydrogen bond acceptor direction and the ideal lone pair direction. This approach adds flexibility to hydrogen bond directions and produces more realistic RMS values. By using this approach, two novel pharmacophore models were derived (Models 2 and 3) and a hydrogen bond acceptor was added to a previously published NK2 pharmacophore model [Poulsen et al., J. Comput.-Aided Mol. Design, 16 (2002) 273] (Model 1). Model 2 as well as Model 3 are described by seven pharmacophore elements: three hydrophobic groups, three hydrogen bond acceptors and a hydrogen bond donor. Model 1 contains the same hydrophobic groups and hydrogen bond donor as Models 2 and 3, but only one hydrogen bond acceptor. The hydrogen bond acceptors and donor are represented as vectors. Two of the hydrophobic groups are always aromatic rings whereas the other hydrophobic group can be either aromatic or aliphatic. In Model 1 the antagonists bind in an extended conformation with two aromatic rings in a parallel displaced and tilted conformation. Model 2 has the same two aromatic rings in a parallel displaced conformation whereas Model 3 has the rings in an edge to face conformation. The pharmacophore models were evaluated using both a structure (NK receptor homology models) and a ligand based approach. By use of exhaustive conformational analysis (MMFFs force field and the GB/SA hydration model) and least-squares molecular superimposition studies, 21 non-peptide antagonists from several structurally diverse classes were fitted to the pharmacophore models. More antagonists could be fitted to Model 2 with a low RMS and a low conformational energy penalty than to Models 1 and 3. Pharmacophore Model 2 was also able to explain the NK1, NK2 and NK3 subtype selectivity of the compounds fitted to the model. Three NK 7TM receptor models were constructed, one for each receptor subtype. The location of the antagonist binding site in the three NK receptor models is identical. Compounds fitted to pharmacophore Model 2 could be docked into the NK1, NK2 and NK3 receptor models after adjustment of the conformation of the flexible linker connecting the head and tail. Models I and 3 are not compatible with the receptor models.     

表皮生长因子受体酪氨酸激酶抑制剂的药效团研究

根据一系列表皮生长因子受体酪氨酸激酶抑制剂的三维定量构效关系研究，得 到了该类抑制剂的药效团，研究结果与Novartis的药效团模型相当类似．药效团包 括一个氢键受体，一个氢键给体，一个疏水区和一个带有氯或溴原子药效团对于研 究表皮生长因子受体酪氨酸激酶抑制剂结构与活性的关系具有重要的意义．通过三 维数据库搜索可能会得到新的先导化合物．     

The Pharmacophore Hypothesis of Novel Inhibitors for Aurora A Kinase

A three-dimensional pharmacophore model was developed from a series of inhibitors of Aurora A kinase to discover new potent anti-cancer agents using the HypoGen module in the Catalyst software. The pharmacophore model was developed based on the structure of 20 currently available inhibitors, which were carefully selected from the literature. The best hypothesis （Hypo 1） was defined by four features： one hydrogen-bond donor and three hy- drophobic points, with the best correlation coefficient of 0.909, the lowest rms deviation of 1.563, and the highest cost difference of 99.075. The Hypo 1 was then validated by a test set consisting of 24 compounds and by a cross-validation of 95% confidence level through randomizing the data using the CatScramble program, which suggested that a predictive pharmacophore model had been successfully obtained.     

基于药效团的BRD4靶向小分子抑制剂虚拟筛选

BRD4靶点和多种肿瘤密切相关,是具有良好成药性的热门靶点。本文选取活性较好且结构差异较大的BRD4小分子抑制剂作为训练集分子,基于配体小分子共同特征(HipHop)方法使用Discovery Studio 3.0分子模拟软件构建了药效团。药效团通过测试集验证、ROC曲线验证(SE(sensitivity)=0.93765、SP(specificity)=0.89500、(AUC)=0.956),结果表明构建得到的药效团具有较强的可靠性和较高的可信度。药效团模型含有1个芳环中心、1个疏水基团、2个氢键受体四个药效特征元素。此药效团被用于ZINC数据库进行虚拟筛选,共筛选了861203个分子,命中率为0.782%。再对筛选得到的分子经过分子对接、ADMET成药性预测、构象分析并讨论分子-蛋白相互作用模式,最终得到了21个有潜力的BRD4小分子抑制剂。     

氧肟酸类组蛋白去乙酰化酶抑制剂的药效团模型研究

基于24个目前已知的氧肟酸类组蛋白去乙酰化酶抑制剂,我们运用Catalyst软件建立了一个三维药效团模型。其中,最好的药效团模型1,包含了四个化学特征（一个氢键供体,一个芳环和两个疏水基）,相关系数达到0.946,并由另外20个化合物进行了测试验证。我们第一次特征性描述了组蛋白去乙酰化酶的帽子（CAP）部分。我们的研究结果对于设计全新组蛋白去乙酰化酶抑制剂具有很好的指导作用。