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

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

Pharmacophore Mode lingand Virtual Screening to Design the Potential Influenza Virus Endonuclease Inhibitors

Influenza virus endonuclease is an attractive target for antiviral therapy in the treatment of influenza infection. The purpos e of this study is to design a novel antiviral agent with improved biological activities against the influenza virus endonuclease. In this study, chemical feature‐based 3D pharmacophore models were developed from 41 known influenza virus endonuclease inhibitors. The best quantitative pharmacohore model (Hypo 1), which consists of two hydrogen‐bond acceptors and two hydrophobic features, yields the highest correlation coefficient (R = 0.886). Hypo 1 was further validated by the cross validation method and the test set compounds. The application of this model for predicting the activities of 11 known influenza virus endonuclease inhibitors in the test set shows great success. The correlation coefficient of 0.942 and a cross validation of 95;% confidence level prove that this model is reliable in identifying structurally diverse compounds for influenza virus endonuclease inhibition. The most active compound (compound 1) from the training set was docked into the active site of the influenza virus endonuclease as an additional verification that the pharmacophore model is accurate. The docked conformation showed important hydrogen bond interactions between the compound and two amino acids, Lys 134 and Lys 137. After validation, this model was used to screen the NCI chemical database to identify new influenza virus endonuclease inhibitors. Our study shows that the to pranking compound out of the 10 newly identified compounds using fit value ranking has an estimated activity of 0.049 μM. These newly identified lead compounds can be further experimentally validated using in vitro techniques.     

5-羟色胺(5-HT)和去甲肾上腺素(NE)重摄取抑制剂药效团模型的建立及比较分析

近年来的研究表明, 去甲肾上腺素(NE)能系统和5-羟色胺(5-HT)能系统可能共同参与了抑郁症的发病机制. 采用Catalyst软件的Hypogen方法, 利用22个不同结构类型的5-HT重摄取抑制剂和19个不同结构类型和活性的NE重摄取抑制剂分别建立了5-HT药效团模型和NE药效团模型, 它们的相关系数分别为0.935, 0.844, 这表明所得到的模型能较好地表征重摄取抑制剂化合物的特征; 此外, 我们还选择了四种不同活性的预测集分别对所建立的药效团模型进行检验, 结果表明所建立的药效团模型具有较好的预测能力. 对这两个药效团模型进行了比较分析, 其结果可以为设计高活性的双重5-HT和NE重摄取抑制剂提供依据.     

A Novel Pharmacophore Model Derived from a Class of Capsid Protein Enterovirus 71 Inhibitors

Capsid protein enterovirus 71 (EV71) is one of the major viruses that cause the severe encephalitis and thus result in a high mortality in children less than 5 years of age.In an effort to discover new potent inhibitors against EV71,a novel three-dimensional pharmacophore model was developed on 24 inhibitors with different molecular structures and bioactivities.The best hypothesis (Hypo1) has a high predictive power and consists of four features,namely,one hydrophobic point (HY) and three hydrogen-bond acceptors (HA).Two key features of the best Hypo1,HY1 and HA3 match well with an important narrow hydrophobic canyon and with the surface of LYS274 in the target EV71 active site,respectively.The more versatile feature,HA1,is firstly found to be very influential on these compounds' bioactivities,which may interact with the other side of the active site in the EV71 receptor.The application of the model is successful in predicting the activities of 30 known EV71 inhibitors with a correlation coefficient of 0.831.Furthermore,Hypo1 demonstrates a superior screening capability for retrieving inhibitors from the database with a high enrichment factor of 70.This study provides some important clues in search for more potent inhibitors against EV71 infection.     

Modeling ligand-receptor interaction for some MHC class II HLA-DR4 peptide mimetic inhibitors using several molecular docking and 3D QSAR techniques

The ligand-receptor interaction between some peptidomimetic inhibitors and a class II MHC peptide presenting molecule, the HLA-DR4 receptor, was modeled using some three-dimensional (3D) quantitative structure-activity relationship (QSAR) methods such as the Comparative Molecular Field Analysis (CoMFA), Comparative Molecular Similarity Indices Analysis (CoMSIA), and a pharmacophore building method, the Catalyst program. The structures of these peptidomimetic inhibitors were generated theoretically, and the conformations used in the 3D QSAR studies were defined by docking them into the known structure of HLA-DR4 receptor through the GOLD, GLIDE Rigidly, GLIDE Flexible, and Xscore programs. Some of the parameters used in these docking programs were selected by docking an X-ray ligand into the receptor and comparing the root-means-square difference (RMSD) computed between the coordinates of the X-ray and docked structure. However, the goodness of a docking result for docking a series of peptidomimetic inhibitors into the HLA-DR4 receptor was judged by comparing the Spearman's rank correlation coefficient computed between each docking result and the activity data taken from the literature. The best CoMFA and CoMSIA models were constructed using the aligned structures of the best docking result. The CoMSIA was conducted in a stepwise manner to identify some important molecular features that were further employed in a pharmacophore building process by the Catalyst program. It was found that most inhibitors of the training set were accurately predicted by the best pharmacophore model, the Hypo1 hypothesis constructed. The deviation or conflict found between the actual and predicted activities of some inhibitors of both the training and the test sets were also investigated by mapping the Hypo1 hypothesis onto the corresponding structures of the inhibitors.     

Pharmacophore modeling coupled with molecular dynamic simulation approach to identify new leads for meprin-β metalloprotease

Human meprin beta metalloprotease, a small subgroup of the astacin family, is a potent drug target for the treatment of several disorders such as fibrosis, neurodegenerative disease in particular Alzheimer and inflammatory bowel diseases. In this study, a ligand-based pharmacophore approach has been used for the selection of potentially active compounds to understand the inhibitory activities of meprin-β by using the sulfonamide scaffold based inhibitors. Using this dataset, a pharmacophore model (Hypo1) was selected on the basis of a highest correlation coefficient (0.959), lowest total cost (105.89) and lowest root mean square deviation (1.31 Å) values. All the pharmacophore hypotheses generated from the candidate inhibitors comprised four features: two hydrogen-bond acceptor, one hydrogen-bond donor and one zinc binder feature. The best validated pharmacophore model (Hypo1) was used for virtual screening of compounds from several databases. The selective hit compounds were filtered by drug likeness property, acceptable ADMET profile, molecular docking and DFT study. Molecular dynamic simulations with the final 10 hit compounds revealed that a large number of non-covalent interactions were formed with the active site and specificity sub-pockets of the meprin beta metalloprotease. This study assists in the development of the new lead molecules as well as gives a better understanding of their interaction with meprin-β.     

Evaluation of potential inhibitors of squalene synthase based on virtual screening and in vitro studies

Squalene synthase (SQS) is a potential target for hyperlipidemia treatment. To identify novel chemical scaffolds of SQS inhibitors, we generated 3D-QSAR pharmacophore models using HypoGen. The best quantitative pharmacophore model, Hypo 1, was selected for virtual screening using two chemical databases, Specs and Traditional Chinese Medicine database (TCM). The best-mapped hit compounds were then subjected to filtering by Lipinski^’s rule of five and docking studies to refine the hits. Finally, five compounds were selected from the top-ranked hit compounds for SQS inhibitory assay in vitro. Three of these compounds could inhibit SQS in vitro, and should be further evaluated pre-clinically as a treatment for hyperlipidemia.     

Identification of Influenza PAN Endonuclease Inhibitors via 3D-QSAR Modeling and Docking-Based Virtual Screening

Structural and biochemical studies elucidate that PA_N may contribute to the host protein shutdown observed during influenza A infection. Thus, inhibition of the endonuclease activity of viral RdRP is an attractive approach for novel antiviral therapy. In order to envisage structurally diverse novel compounds with better efficacy as PA_N endonuclease inhibitors, a ligand-based-pharmacophore model was developed using 3D-QSAR pharmacophore generation (HypoGen algorithm) methodology in Discovery Studio. As the training set, 25 compounds were taken to generate a significant pharmacophore model. The selected pharmacophore Hypo1 was further validated by 12 compounds in the test set and was used as a query model for further screening of 1916 compounds containing 71 HIV-1 integrase inhibitors, 37 antibacterial inhibitors, 131 antiviral inhibitors and other 1677 approved drugs by the FDA. Then, six compounds (Hit01–Hit06) with estimated activity values less than 10 μM were subjected to ADMET study and toxicity assessment. Only one potential inhibitory ‘hit’ molecule (Hit01, raltegravir’s derivative) was further scrutinized by molecular docking analysis on the active site of PA_N endonuclease (PDB ID: 6E6W). Hit01 was utilized for designing novel potential PA_N endonuclease inhibitors through lead optimization, and then compounds were screened by pharmacophore Hypo1 and docking studies. Six raltegravir’s derivatives with significant estimated activity values and docking scores were obtained. Further, these results certainly do not confirm or indicate the seven compounds (Hit01, Hit07, Hit08, Hit09, Hit10, Hit11 and Hit12) have antiviral activity, and extensive wet-laboratory experimentation is needed to transmute these compounds into clinical drugs.     

基于吡咯烷与正丁烷类衍生物CCR5拮抗剂的药效团模型构建

吡咯烷与正丁烷类CCR5(化学趋化因子受体5)拮抗剂可通过抑制人类免疫缺陷病毒(HIV-1)包膜蛋白与CCR5的相互作用而阻断病毒进入细胞. 本文使用已知拮抗剂结构和活性信息构建了一个三维药效团模型. 按照Catalyst/HypoGen模块的要求, 选择了25个结构和活性均具备差异性的分子作为药效团产生的训练集. 其中训练集分子以IC50值表示的生物活性值跨度为0.06到10000 nmol·L-1. 最好的药效团模型(Hypo 1)由两个正离子化特征以及三个疏水特征组成, 训练集预测相关系数为0.924, 均方根偏差为1.068. 模型用于预测由74个分子组成的测试集化合物活性, 结果表明模型可以提供较好的活性预测结果并用于新的拮抗剂的设计.     

Pharmacophore Model Generation Based on Pyrrolidine- and Butane-derived CCR5 Antagonists

A three-dimensional pharmacophore model was developed for a considerable number of pyrrolidine-based and butane-based chemokine (C-C motif) receptor 5 (CCR5) antagonists, which can block the entry of human immunodeficiency virus type 1 (HIV-1) by inhibiting the interaction of HIV-1 envelope protein and CCR5. The pharmacophore model was generated using a training set consisting of 25 carefully selected antagonists with the diverse molecular architecture and bioactivity, as required by the Catalyst/HypoGen program. The activity of the training set molecules expressed in IC₅₀ (half-inhibitory concentration) covered from 0.06 to 10000 nmol·L^–1. The most predictive pharmacophore model (Hypo 1), consisting of two positive ionizable points and three hydrophobic groups, had a correlation of 0.924 and a root mean square of 1.068, and a cost difference of 63.67 bits between the null cost and the total cost. The model was applied in predicting the activity of 74 compounds as a test set. The results indicated that the model was able to provide clear guidelines and accurate activity prediction for novel antagonist design.     

Design of novel FLT-3 inhibitors based on dual-layer 3D-QSAR model and fragment-based compounds in silico

FMS-like tyrosine kinase 3 (FLT-3) is strongly correlated with acute myeloid leukemia, but no FLT-3-inhibitor cocomplex structure is available to assist the design of therapeutic inhibitors. Hence, we propose a dual-layer 3D-QSAR model for FLT-3 that integrates the pharmacophore, CoMFA, and CoMSIA. We then coupled the model with the fragment-based design strategy to identify novel FLT-3 inhibitors. In the first layer, the previously established model, Hypo02, was evaluated in terms of its correlation coefficient (r), RMS, cost difference, and configuration cost, with values of 0.930, 1.24, 106.45, and 16.44, respectively. Moreover, Fischer's cross-validation test of data generated by Hypo02 yielded a 98% confidence level, and the validation of the testing set yielded a best r value of 0.87. The features of Hypo02 were separated into two parts and then used to screen the MiniMaybridge fragment compound database. Nine novel FLT-3 inhibitors were generated in this layer. In the second layer, Hypo02 was subjected to an alignment rule to generate CoMFA- and CoMSIA-based models, for which the partial least-squares validation method was utilized. The values of q(2), r(2), and predictive r(2) were 0.58, 0.98, and 0.76, respectively, derived from the CoMFA model with steric and electrostatic fields. The CoMSIA model with five different fields yielded values of 0.54, 0.97, and 0.76 for q(2), r(2), and predictive r(2), respectively. The CoMFA and CoMSIA models were used to constrain 3D structures of the nine novel FLT-3 inhibitors. This dual-layer 3D-QSAR model constitutes a valuable tool to easily and quickly screen and optimize novel potential FLT-3 inhibitors for the treatment of acute myeloid leukemia.     

Pharmacophore Modeling and Virtual Screening of Novel Inhibitors for c-Kit Kinase

The stem cell factor receptor (c‐Kit) has been known to play critical roles in regulating numerous aspects of cellular behavior including cell growth, differentiation, migration and metabolism. In this investigation, a three‐dimensional pharmacophore model of c‐Kit inhibitors has been established by using the HypoGen algorithms implemented in the catalyst software package. The best quantitative pharmacophore model, hypothesis 1, which has the highest correlation coefficient (0.989), consists of one hydrogen bond acceptor, two hydrogen bond donors and one hydrophobic feature. To our knowledge, this is the first report on the pharmacophore modeling study of c‐Kit inhibitors. The best hypothesis, hypothesis 1, was used to screen molecular structural databases, including Specs and China Natural Products Database for potential lead compounds. The hit compounds were subsequently subjected to filtering by Lipinski's rules and docking study to refine the retrieved hits and as a result to reduce the rate of false positive. Finally 28 compounds were purchased or synthesized for further in vitro assay against several human tumour cell lines including A549, MCF‐7, HepG2 and PC‐3, in which c‐Kit is overexpressed. Two compounds show very low micromolar inhibition potency against the PC‐3 and HepG2 cell lines respectively. And they were selected for further modification and testing.     

Accounting of ligand–receptor interactions to explore and design novel architecture for PTP 1B inhibition: a legitimate approach

Computer‐aided drug design was performed on a diverse set of 103 biphenyl derivatives that demonstrated antidiabetic activity by restraining the protein tyrosine phosphatase 1B (PTP 1B) receptor. A four‐point pharmacophore hypothesis using the PHASE module of Schrödinger suite with one hydrogen bond acceptor (A) and three aromatic rings (R) as pharmacophoric features was generated. The hypothesis, ARRR.2, considered the best hypothesis in the present study is characterized by survival score (3.553), cross‐validated r² (Q²) (0.722), regression coefficient (0.949), Pearson R (0.867), and F value (492.6). The developed pharmacophore model was externally validated by predicting the activity of test set molecules. Docking algorithm combined with the drug–receptor binding free energetic and pharmacokinetic drug profile envisaged a novel concept, which may provide structural insights for the development of potential PTP 1B inhibitors. The study also provided a valid rapport between pharmacophore drug mapping, atom‐based three‐dimensional quantitative structure–activity relationship, molecular docking, sitemap, molecular simulations, and pharmacokinetic prediction approaches demonstrating the trends in activity. The results of these ligand–receptor relationship studies may account to design a legitimate template for the development and optimization of highly selective and potent PTP 1B inhibitors. Copyright © 2012 John Wiley & Sons, Ltd.     

Development of predictive pharmacophore model for in silico screening,and 3D QSAR CoMFA and CoMSIA studies for lead optimization,for designing of potent tumor necrosis factor alpha converting enzyme inhibitors

A chemical feature-based pharmacophore model was developed for Tumor Necrosis Factor-α converting enzyme (TACE) inhibitors. A five point pharmacophore model having two hydrogen bond acceptors (A), one hydrogen bond donor (D) and two aromatic rings (R) with discrete geometries as pharmacophoric features was developed. The pharmacophore model so generated was then utilized for in silico screening of a database. The pharmacophore model so developed was validated by using four compounds having proven TACE inhibitory activity which were grafted into the database. These compounds mapped well onto the five listed pharmacophoric features. This validated pharmacophore model was also used for alignment of molecules in CoMFA and CoMSIA analysis. The contour maps of the CoMFA/CoMSIA models were utilized to provide structural insight for activity improvement of potential novel TACE inhibitors. The pharmacophore model so developed could be used for in silico screening of any commercial/in house database for identification of TACE inhibiting lead compounds, and the leads so identified could be optimized using the developed CoMSIA model. The present work highlights the tremendous potential of the two mutually complementary ligand-based drug designing techniques (i.e. pharmacophore mapping and 3D-QSAR analysis) using TACE inhibitors as prototype biologically active molecules.     

基于药效团模型的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个预测活性较好的化合物,可作为进一步药物研发的候选化合物.     

Homology modeling,docking and structure-based pharmacophore of inhibitors of DNA methyltransferase

DNA methyltransferase 1 (DNMT1) is an emerging epigenetic target for the treatment of cancer and other diseases. To date, several inhibitors from different structural classes have been published. In this work, we report a comprehensive molecular modeling study of 14 established DNTM1 inhibitors with a herein developed homology model of the catalytic domain of human DNTM1. The geometry of the homology model was in agreement with the proposed mechanism of DNA methylation. Docking results revealed that all inhibitors studied in this work have hydrogen bond interactions with a glutamic acid and arginine residues that play a central role in the mechanism of cytosine DNA methylation. The binding models of compounds such as curcumin and parthenolide suggest that these natural products are covalent blockers of the catalytic site. A pharmacophore model was also developed for all DNMT1 inhibitors considered in this work using the most favorable binding conformations and energetic terms of the docked poses. To the best of our knowledge, this is the first pharmacophore model proposed for compounds with inhibitory activity of DNMT1. The results presented in this work represent a conceptual advance for understanding the protein–ligand interactions and mechanism of action of DNMT1 inhibitors. The insights obtained in this work can be used for the structure-based design and virtual screening for novel inhibitors targeting DNMT1.     

具有全新机理的DNA旋转酶抑制剂的筛选及抑菌活性

利用具有新机制的抗耐药菌DNA旋转酶抑制剂GSK299423与DNA旋转酶的晶体复合物(PDB code:2XCS)构建基于配体-受体复合物的药效团模型, 诱骗集(Decoy set)验证结果表明该药效团模型具有较强的活性识别能力. 将药效团模型与分子对接相结合用于筛选化合物库, 通过抑菌活性测定, 获得了具有抗多药耐药菌活性的DNA旋转酶抑制剂LTH02.     

A ligand-based molecular modeling study on some matrix metalloproteinase-1 inhibitors using several 3D QSAR techniques

Some three-dimensional quantitative structure-activity relationship (3D-QSAR) models have been constructed using the comparative molecular field analysis (CoMFA) and comparative molecular similarity indices (CoMSIA) for a series of 84 proline-based plus 12 structurally more diversified nonproline matrix metalloproteinase inhibitors. The structures of these inhibitors were built from a structure template extracted from the crystal structure of stromelysin. The structures built were divided into the training and test sets for both the CoMFA and CoMSIA analyses for each being composed of 60 and 24 inhibitors, respectively. The structures in the training set were aligned using some alignment rules derived from the analysis of the Ligplot program on a recent crystal structure of ligand-collagenase-1 complex. Some stepwise CoMSIA's were performed on the aligned training set on which the best CoMFA result was obtained. The best CoMSIA model was identified from the stepwise results, and the corresponding pharmacophore features were used for the construction of a pharmacophore hypothesis by the Catalyst 4.9 program. The training set was extended to include 11 structurally more diversified and nonproline inhibitors. To construct a pharmacophore hypothesis, the conformation of 60 structurally aligned proline-based inhibitors was fixed, while that of the 11 structurally more diversified nonproline inhibitors was allowed to vary during the hypothesis construction process. It was found that the predicted activities by the top hypothesis constructed for both the training and test sets were as good in statistics as those predicted by the best CoMSIA model from which the hypothesis was derived. The top hypothesis was mapped onto the structures of several highly active inhibitors selected from both the training and test sets. The goodness of mapping on each inhibitor was found to be correlated well with the activity of each inhibitor.     

Identification of critical chemical features for Aurora kinase-B inhibitors using Hip-Hop, virtual screening and molecular docking

This study was performed to find the selective chemical features for Aurora kinase-B inhibitors using the potent methods like Hip-Hop, virtual screening, homology modeling, molecular dynamics and docking. The best hypothesis, Hypo1 was validated toward a wide range of test set containing the selective inhibitors of Aurora kinase-B. Homology modeling and molecular dynamics studies were carried out to perform the molecular docking studies. The best hypothesis Hypo1 was used as a 3D query to screen the chemical databases. The screened molecules from the databases were sorted based on ADME and drug like properties. The selective hit compounds were docked and the hydrogen bond interactions with the critical amino acids present in Aurora kinase-B were compared with the chemical features present in the Hypo1. Finally, we suggest that the chemical features present in the Hypo1 are vital for a molecule to inhibit the Aurora kinase-B activity.