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.     

Pharmacophore Modeling and in Silico Screening Studies to Design Potential KDR Kinase Inhibitors

A novel ligand‐based pharmacophore model for KDR kinase was generated on the basis of chemical features of 30 KDR kinase inhibitors. This pharmacophore model consists of one hydrogen‐bond acceptor, one hydrogen‐bond donor and two hydrophobic groups. Several methods have been used to validate the model, suggesting that it can serve as a reliable tool for virtual screening to facilitate the discovery of novel KDR inhibitors. The model was then used as database search query from the National Cancer Institute (NCI) database for the rational design to identify new hit compound.     

Multiple receptor-ligand based pharmacophore modeling and molecular docking to screen the selective inhibitors of matrix metalloproteinase-9 from natural products

Matrix metalloproteinase-9 (MMP-9) is an attractive target for cancer therapy. In this study, the pharmacophore model of MMP-9 inhibitors is built based on the experimental binding structures of multiple receptor-ligand complexes. It is found that the pharmacophore model consists of six chemical features, including two hydrogen bond acceptors, one hydrogen bond donor, one ring aromatic regions, and two hydrophobic (HY) features. Among them, the two HY features are especially important because they can enter the S1′ pocket of MMP-9 which determines the selectivity of MMP-9 inhibitors. The reliability of pharmacophore model is validated based on the two different decoy sets and relevant experimental data. The virtual screening, combining pharmacophore model with molecular docking, is performed to identify the selective MMP-9 inhibitors from a database of natural products. The four novel MMP-9 inhibitors of natural products, NP-000686, NP-001752, NP-014331, and NP-015905, are found; one of them, NP-000686, is used to perform the experiment of in vitro bioassay inhibiting MMP-9, and the IC₅₀ value was estimated to be only 13.4 µM, showing the strongly inhibitory activity of NP-000686 against MMP-9, which suggests that our screening results should be reliable. The binding modes of screened inhibitors with MMP-9 active sites were discussed. In addition, the ADMET properties and physicochemical properties of screened four compounds were assessed. The found MMP-9 inhibitors of natural products could serve as the lead compounds for designing the new MMP-9 inhibitors by carrying out structural modifications in the future.     

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.     

Virtual screening studies to design potent CDK2-cyclin A inhibitors

The cell division cycle is controlled by cyclin-dependent kinases (CDK), which consist of a catalytic subunit (CDK1-CDK8) and a regulatory subunit (cyclin A-H). Pharmacophore analysis indicates that the best inhibitor model consists of (1) two hydrogen bond acceptors, (2) one hydrogen bond donor, and (3) one hydrophobic feature. The HypoRefine pharmacophore model gave an enrichment factor of 1.31 and goodness of fit score of 0.76. Docking studies were carried out to explore the structural requirements for the CDK2-cyclin A inhibitors and to construct highly predictive models for the design of new inhibitors. Docking studies demonstrate the important role of hydrogen bond and hydrophobic interactions in determining the inhibitor-receptor binding affinity. The validated pharmacophore model is further used for retrieving the most active hits/lead from a virtual library of molecules. Subsequently, docking studies were performed on the hits, and novel series of potent leads were suggested based on the interaction energy between CDK2-cyclin A and the putative inhibitors.     

3D-QSAR studies on indole and 7-azoindole derivatives as ROCK-2 inhibitors: An integrative computational approach

Rho Kinases (ROCK) has been found to regulate a wide range of fundamental cell functions such as contraction, motility, proliferation, and apoptosis. Recent experiments have defined new functions of ROCKs in cells, including centrosome positioning and cell-size regulation, which might contribute to various physiological and pathological states. In this study, we have performed pharmacophore modeling and 3D QSAR studies on a series of 36 indoles and 7-azoindoles derivatives as ROCK2 inhibitors to elucidate the structural variations with their inhibitory activities. Ligand based CoMFA and CoMSIA models were generated based on three different alignment methods such as systematic search, simulated annealing and pharmacophore. A total of 15 CoMFA models and 27 CoMSIA were generated using different alignments. One model from each alignment is selected based on the statistical values. Contour maps of the selected models were compared, analysed and reported. The 3D QSAR study revealed that electro positive group linked to the methoxy-benzene ring position of the structure will enhance the biological activity and bulkier substitutions are preferred in the methyl dihydroindole region. Also, it is found that the hydrogen bond donor substituted at the R₁ position enhances the inhibitory activity. In future, this study would give proper guidelines to further enhance the activity of novel inhibitors for ROCK2.     

黄酮类醛糖还原酶抑制剂的抑制机理研究

为了研究黄酮类醛糖还原酶抑制剂的抑制机理, 选择了31个黄酮类化合物作为训练集, 使用Catalyst软件包构建了此类抑制剂的药效团模型. 并专门针对黄酮类化合物定制了氢键给体和受体模型, 效果优于使用Catalyst内预定义的模型. 最终的药效团模型由两个氢键给体和一个氢键受体组成, 对训练集具有较好预测能力(Correl=0.9013). 此外, 使用InsightII/Affinity对6个黄酮类化合物进行了分子对接研究. 综合药效团模型和分子对接研究的结果, 发现黄酮类化合物的抑制活性主要源于黄酮骨架上的C4’或C3’位的羟基与醛糖还原酶活性口袋中的TYR48、VAL47、GLN49和C7位的羟基与HIS110, TRP111所形成的两组氢键.     

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.     

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.     

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

BRD4靶点和多种肿瘤密切相关,是具有良好成药性的热门靶点.本文选取活性较好且结构差异较大的BRD4小分子抑制剂作为训练集分子,基于配体小分子共同特征(HipHop)方法使用Discovery Studio 3.0分子模拟软件构建了药效团.药效团通过测试集验证、ROC曲线验证(SE(sensitivity)=0.937...     

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

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

3D-pharmacophore models for selective A2A and A2B adenosine receptor antagonists

Three-dimensional pharmacophore models were generated for A2A and A2B adenosine receptors (ARs) based on highly selective A2A and A2B antagonists using the Catalyst program. The best pharmacophore model for selective A2A antagonists (Hypo-A2A) was obtained through a careful validation process. Four features contained in Hypo-A2A (one ring aromatic feature (R), one positively ionizable feature (P), one hydrogen bond acceptor lipid feature (L), and one hydrophobic feature (H)) seem to be essential for antagonists in terms of binding activity and A2A AR selectivity. The best pharmacophore model for selective A2B antagonists (Hypo-A2B) was elaborated by modifying the Catalyst common features (HipHop) hypotheses generated from the selective A2B antagonists training set. Hypo-A2B also consists of four features: one ring aromatic feature (R), one hydrophobic aliphatic feature (Z), and two hydrogen bond acceptor lipid features (L). All features play an important role in A2B AR binding affinity and are essential for A2B selectivity. Both A2A and A2B pharmacophore models have been validated toward a wide set of test molecules containing structurally diverse selective antagonists of all AR subtypes. They are capable of identifying correspondingly high potent antagonists and differentiating antagonists between subtypes. The results of our study will act as a valuable tool for retrieving structurally diverse compounds with desired biological activities and designing novel selective adenosine receptor ligands.     

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.     

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

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

Construction of the Pharmacophore Model of Glycogen Synthase Kinase-3 Inhibitors

A three dimensional pharmacophore model has been generated for glycogen synthase kinase-3 （GSK-3） inhibitots. A dataset consisting of 89 compounds was selected on the basis of the information content of the structures and activity data as required by the CATALYST system. The optimum model with four features （one hydrogen-bond acceptor unit, one ring aromatic unit, and two hydrophobic aromatic units） was selected with a good correlation coefficient （0.95）. This model is able to predict the activity of other known GSK-3 inhibitors not included in the model generation, and can be further used to identify structurally diverse compounds with desired biological activity by virtual screening.     

Mangiferin, an anti-HIV-1 agent targeting protease and effective against resistant strains

The anti-HIV-1 activity of mangiferin was evaluated. Mangiferin can inhibit HIV-1(Ⅲ)(B) induced syncytium formation at non-cytotoxic concentrations, with a 50% effective concentration (EC??) at 16.90 μM and a therapeutic index (TI) above 140. Mangiferin also showed good activities in other laboratory-derived strains, clinically isolated strains and resistant HIV-1 strains. Mechanism studies revealed that mangiferin might inhibit the HIV-1 protease, but is still effective against HIV peptidic protease inhibitor resistant strains. A combination of docking and pharmacophore methods clarified possible binding modes of mangiferin in the HIV-1 protease. The pharmacophore model of mangiferin consists of two hydrogen bond donors and two hydrogen bond acceptors. Compared to pharmacophore features found in commercially available drugs, three pharmacophoric elements matched well and one novel pharmacophore element was observed. Moreover, molecular docking analysis demonstrated that the pharmacophoric elements play important roles in binding HIV-1 protease. Mangiferin is a novel nonpeptidic protease inhibitor with an original structure that represents an effective drug development strategy for combating drug resistance.     

Identification of potential HIV-1 integrase strand transfer inhibitors: In silico virtual screening and QM/MM docking studies

HIV-1 integrase (IN) is a retroviral enzyme that catalyses integration of the reverse-transcribed viral DNA into the host genome, which is necessary for efficient viral replication. In this study, we have performed an in silico virtual screening for the identification of potential HIV-1 IN strand transfer (ST) inhibitors. Pharmacophore modelling and atom-based 3D-QSAR studies were carried out for a series of compounds belonging to 3-Hydroxypyrimidine-2,4-diones. Based on the ligand-based pharmacophore model, we obtained a five-point pharmacophore with two hydrogen bond acceptors (A), one hydrogen bond donor (D), one hydrophobic group (H) and one aromatic ring (R) as pharmacophoric features. The pharmacophore hypothesis AADHR was used as a 3D query in a sequential virtual screening study to filter small molecule databases Maybridge, ChemBridge and Asinex. Hits matching with pharmacophore hypothesis AADHR were retrieved and passed progressively through Lipinski’s rule of five filtering, molecular docking and hierarchical clustering. The five compounds with best hits with novel and diverse chemotypes were subjected to QM/MM docking, which showed improved docking accuracy. We further performed molecular dynamics simulation and found three compounds that form stable interactions with key residues. These compounds could be used as a leads for further drug development and rational design of HIV-1 IN inhibitors.     

Identification of new dual spleen tyrosine kinase (Syk) and phosphoionositide-3-kinase δ (PI3Kδ) inhibitors using ligand and structure-based integrated ideal pharmacophore models

Owing to the complex pathophysiology of autoimmune disorders, it is very challenging to develop successful treatment strategies. Single-target agents are not desired therapeutics for such multi-factorial disorders. Considering the current need for the treatment of complex autoimmune disorders, dual inhibitors of Syk and PI3Kδ have been designed using ligand and structure-based molecular modelling strategies. In the present work, structure and ligand-based pharmacophore modelling was implemented for a varied set of Syk and PI3Kδ inhibitors. Ligand-based pharmacophore models (LBPMs) were developed for two kinases: ADPR.14 (r²_train = 0.809) for Syk, comprising one hydrogen bond acceptor, one hydrogen bond donor, one positive ionisable and one ring aromatic feature, and for PI3Kδ: AAARR.45 (r²_train = 0.942) consisting of three hydrogen bond acceptor and two ring aromatic features. The generated e-pharmacophore models revealed an additional ring aromatic and hydrophobic feature important for Syk and PI3Kδ inhibition, respectively. Subsequently, LBPMs were modified resulting in APDRR.14 hypothesis for Syk inhibitors and AAAHRR.45 hypothesis for PI3Kδ inhibitors employed for virtual screening. Thus, the combination of ligand and structure-based pharmacophore modelling helped in developing ideal pharmacophore models that may be an efficient tool for the designing of novel dual inhibitors of Syk and PI3Kδ.     

Structure-based design and optimization of antihypertensive peptides to obtain high inhibitory potency against both renin and angiotensin I-converting enzyme

The human renin–angiotensin system (RAS) plays an essential role in regulating blood pressure and systemic vascular resistance. Renin and angiotensin I-converting enzyme (ACE) are two key enzymes in RAS and have long been recognized as attractive antihypertensive targets. Here, a synthetic strategy was proposed integrating quantitative structure–activity relationship (QSAR), molecular dynamics (MD) simulation and binding free energy analysis to discover novel dual renin and ACE peptidic inhibitors. With the strategy a number of candidates were generated virtually, from which eight promising peptides were selected and synthesized for biological assay. Consequently, three peptides (RYLP, YTAWVP and YRAWVL) were successfully identified to have satisfactory inhibitory profile against both renin and ACE with IC₅₀ values of <1 mM and <10 μM, respectively. Structural analysis and energetic dissection revealed different binding modes of peptide to renin and ACE; a peptide only inserts its C-terminus into the active site of ACE, whereas the whole peptide packs tightly against renin. In addition, when limited to structural diversity it is hard to reconcile the renin and ACE inhibitory activities of short peptides such as dipeptides. These findings can be used to guide peptide optimization with improved biological activity.     

Comparative pharmacophore modeling of human adenosine receptor A1 and A3 antagonists

Adenosine receptors are promising therapeutic targets in drug discovery. In this study, three-dimensional pharmacophore models of human adenosine receptor A1 and A3 antagonists were developed based on 26 and 23 diverse compounds, respectively. The best A1 pharmacophore model (A 1 _Hopy1) consists of four features: one hydrogen bond donor, one hydrophobic point and two ring aromatics, while the best A 3 pharmacophore model (A3 _Hopy1) also has four features: one hydrogen bond acceptor, one hydrophobic point and two ring aromatics. The correlation coefficients were 0.840 for A 1 test set with 146 diverse compounds and 0.827 for A3 test set with 238 diverse compounds. In the simulated virtual screening experiments, high enrichment factors of 6.51 and 6.90 were obtained for A 1 _Hopy1 and A3 _Hopy1 models, respectively. Moreover, two models also showed high subtype-selectivity in the simulated virtual screening experiments. These results could be helpful for the discovery of novel potent and selective A 1 and A3 antagonists.