A Holistic Evolutionary and 3D Pharmacophore Modelling Study Provides Insights into the Metabolism,Function, and Substrate Selectivity of the Human Monocarboxylate Transporter 4 (hMCT4)

Monocarboxylate transporters (MCTs) are of great research interest for their role in cancer cell metabolism and their potential ability to transport pharmacologically relevant compounds across the membrane. Each member of the MCT family could potentially provide novel therapeutic approaches to various diseases. The major differences among MCTs are related to each of their specific metabolic roles, their relative substrate and inhibitor affinities, the regulation of their expression, their intracellular localization, and their tissue distribution. MCT4 is the main mediator for the efflux of L-lactate produced in the cell. Thus, MCT4 maintains the glycolytic phenotype of the cancer cell by supplying the molecular resources for tumor cell proliferation and promotes the acidification of the extracellular microenvironment from the co-transport of protons. A promising therapeutic strategy in anti-cancer drug design is the selective inhibition of MCT4 for the glycolytic suppression of solid tumors. A small number of studies indicate molecules for dual inhibition of MCT1 and MCT4; however, no selective inhibitor with high-affinity for MCT4 has been identified. In this study, we attempt to approach the structural characteristics of MCT4 through an in silico pipeline for molecular modelling and pharmacophore elucidation towards the identification of specific inhibitors as a novel anti-cancer strategy.     

乙酰胆碱M3受体拮抗剂的从头设计

M3受体是人体内分布广泛的一种重要的毒蕈碱乙酰胆碱受体亚型,与膀胱过度活动症、心律失常、呼吸道及支气管疾病密切相关.设计与合成M3受体拮抗剂对于开发治疗相关疾病的药物具有重要理论意义和应用价值。本文采用同源模建方法,构建了人M3受体蛋白的三维结构,基于M3受体拮抗剂的结构构建了其药效团模型,其中较理想的模型含有6个药效团特征元素。利用药效团模型进行了虚拟筛选,虚拟筛选的数据库是本实验室构建的虚拟化合物库,发现了10个新型的对M3受体有拮抗作用的化合物。最后对这10个化合物进行了分子对接,根据对接结果,发现3个化合物对接能量及结合方式比较合理,为M3受体拮抗剂的合成研究及活性测定奠定了基础。     

流感病毒神经氨酸酶抑制剂的理性筛选

介绍了理性筛选流感病毒神经氨酸酶抑制剂的全过程，共分4个阶段：1)化合物数据库类药性处理；2)建立神经氨酸酶抑制剂三维药效团并对目标数据库进行构象搜索；3)分子对接及对接后分析；4)神经氨酸酶抑制模型的建立及待测化合物的活性检测。活性检测后发现4个活性化合物，其中Ic。为10。M的化合物1个．Ic，。为10^-6M的化合物2个，IC50为10^-7M的化合物1个。应用理性筛选方法，从化合物数据库中挑选出部分化合物进行神经氨酸酶抑制活性的筛选，减少了药物筛选的盲目性，提高了药物发现的机率。     

3D QSAR Studies,Pharmacophore Modeling and Virtual Screening on a Series of Steroidal Aromatase Inhibitors

Aromatase inhibitors are the most important targets in treatment of estrogen-dependent cancers. In order to search for potent steroidal aromatase inhibitors (SAIs) with lower side effects and overcome cellular resistance, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were performed on a series of SAIs to build 3D QSAR models. The reliable and predictive CoMFA and CoMSIA models were obtained with statistical results (CoMFA: q² = 0.636, r²_ncv = 0.988, r²_pred = 0.658; CoMSIA: q² = 0.843, r²_ncv = 0.989, r²_pred = 0.601). This 3D QSAR approach provides significant insights that can be used to develop novel and potent SAIs. In addition, Genetic algorithm with linear assignment of hypermolecular alignment of database (GALAHAD) was used to derive 3D pharmacophore models. The selected pharmacophore model contains two acceptor atoms and four hydrophobic centers, which was used as a 3D query for virtual screening against NCI2000 database. Six hit compounds were obtained and their biological activities were further predicted by the CoMFA and CoMSIA models, which are expected to design potent and novel SAIs.     

Discovery of Dual ETA/ETB Receptor Antagonists from Traditional Chinese Herbs through in Silico and in Vitro Screening

Endothelin-1 receptors (ETAR and ETBR) act as a pivotal regulator in the biological effects of ET-1 and represent a potential drug target for the treatment of multiple cardiovascular diseases. The purpose of the study is to discover dual ETA/ETB receptor antagonists from traditional Chinese herbs. Ligand- and structure-based virtual screening was performed to screen an in-house database of traditional Chinese herbs, followed by a series of in vitro bioassay evaluation. Aristolochic acid A (AAA) was first confirmed to be a dual ETA/ETB receptor antagonist based intracellular calcium influx assay and impedance-based assay. Dose-response curves showed that AAA can block both ETAR and ETBR with IC₅₀ of 7.91 and 7.40 μM, respectively. Target specificity and cytotoxicity bioassay proved that AAA is a selective dual ETA/ETB receptor antagonist and has no significant cytotoxicity on HEK293/ETAR and HEK293/ETBR cells within 24 h. It is a feasible and effective approach to discover bioactive compounds from traditional Chinese herbs using in silico screening combined with in vitro bioassay evaluation. The structural characteristic of AAA for its activity was especially interpreted, which could provide valuable reference for the further structural modification of AAA.     

In silico and in vitro screening to identify structurally diverse non-azole CYP51 inhibitors as potent antifungal agent

The problem of resistance to azole class of antifungals is a serious cause of concern to the medical fraternity and thus there is an urgent need to identify non-azole scaffolds with high affinity for lanosterol 14α-demethylase (CYP51). In view of this we have attempted to identify novel non-azole CYP51 inhibitors through the application of pharmacophore based virtual screening and in vitro evaluation. A rigorously validated pharmacophore model comprising of 2 hydrogen bond acceptor and 2 hydrophobic features has been developed and used to mine NCI database. Out of 265 retrieved hits, NSC 1215 and 1520 have been chosen on the basis of Lipinski’s rule of five, fit and estimated values. Both the hits were docked into the active site of CYP51. In view of high fit value and CDocker score, NSC 1215 and 1520 have been subjected to in vitro microbiological assay. The result reveals that NSC 1215 and 1520 are active against Candida albicans, Candida parapsilosis, Candida tropicalis, and Aspergillus niger. In addition to this the absorption characteristics of both the hits have also been determined using the rat sac technique and permeation in order of NSC 1520 > NSC 1215 has been observed.     

In silico screening and study of novel ERK2 inhibitors using 3D QSAR,docking and molecular dynamics

ERK2 is a dual specificity protein kinase, part of the Ras/Raf/MEK/ERK signal transduction cascade. It forms an interesting target for inhibition based on its relationship with cell proliferation and oncogenesis. A 3D QSAR pharmacophore model (Hypo1) with high correlation (r = 0.938) was developed for ERK2 ATP site on the basis of experimentally known inhibitors. The model included three hydrogen bonds, and one hydrophobic site. Assessment of Hypo1 through Fisher randomization, cost analysis, leave one out method and decoy test suggested that the model can reliably detect ERK2 inhibitors. Hypo1 has been used for virtual screening of potential inhibitors from ZINC, Drug Bank, NCI, Maybridge and Chembank databases. Using Hypo1 as a query, databases have been interrogated for compounds who meet the pharmacophore features. The resulting hit compounds were subject to docking and analysis. Docking and molecular dynamics analysis showed that in order to achieve a higher potency compounds have to interact with catalytic site, glycine rich loop, Hinge region, Gatekeeper region and ATP site entrance residues. We also identified catalytic site and Glycine rich loop as important regions to bind by molecules for better potency and selectivity.     

一个组织型纤溶性溶解酶原激活剂(T-pa)非共价类抑制剂的药效模型

为了研究组织型纤溶性溶解酶原激活剂（T－pa)的非共价类抑制剂的三维定量构效关系，本文采用分子模拟软件Catalyst4.0(molecular simulation company)构建T-pa的非共价类抑制剂的三维药效模型。14个非共价类的抑制剂及其体外活性数据被用于构建此药效模型。此药效模型有3个疏水区特性，一个正离子化区特性和氢键供体区特性，且结构与活性相关系数为r=0.962。与T－pa的晶体结构比较，此药物模型在化学状态和疏水性上与其能很好的相匹配。且利用此是找到的抑制剂（bbzi14)的活性构象，与晶体结构中此抑制剂的结合构象基本一致。通过此药效集团模型反映的抑制剂与受体的相互作用模式，能为发现新型亲和配基和抑制剂提供有用的启示。     

Insight into the Structural Determinants of Imidazole Scaffold-Based Derivatives as TNF-α Release Inhibitors by in Silico Explorations

Presently, 151 widely-diverse pyridinylimidazole-based compounds that show inhibitory activities at the TNF-α release were investigated. By using the distance comparison technique (DISCOtech), comparative molecular field analysis (CoMFA), and comparative molecular similarity index analysis (CoMSIA) methods, the pharmacophore models and the three-dimensional quantitative structure-activity relationships (3D-QSAR) of the compounds were explored. The proposed pharmacophore model, including two hydrophobic sites, two aromatic centers, two H-bond donor atoms, two H-bond acceptor atoms, and two H-bond donor sites characterizes the necessary structural features of TNF-α release inhibitors. Both the resultant CoMFA and CoMSIA models exhibited satisfactory predictability (with Q² (cross-validated correlation coefficient) = 0.557, R²_ncv (non-cross-validated correlation coefficient) = 0.740, R²_pre (predicted correlation coefficient) = 0.749 and Q² = 0.598, R²_ncv = 0.767, R²_pre = 0.860, respectively). Good consistency was observed between the 3D-QSAR models and the pharmacophore model that the hydrophobic interaction and hydrogen bonds play crucial roles in the mechanism of actions. The corresponding contour maps generated by these models provide more diverse information about the key intermolecular interactions of inhibitors with the surrounding environment. All these models have extended the understanding of imidazole-based compounds in the structure-activity relationship, and are useful for rational design and screening of novel 2-thioimidazole-based TNF-α release inhibitors.     

Selective Pharmacophore Models of Dopamine D1 and D2 Full Agonists Based on Extended Pharmacophore Features

This study is focused on the identification of structural features that determine the selectivity of dopamine receptor agonists toward D₁ and D₂ receptors. Selective pharmacophore models were developed for both receptors. The models were built by using projected pharmacophoric features that represent the main agonist interaction sites in the receptor (the Ser residues in TM5 and the Asp in TM3), a directional aromatic feature in the ligand, a feature with large positional tolerance representing the positively charged nitrogen in the ligand, and sets of excluded volumes reflecting the shapes of the receptors. The sets of D₁ and D₂ ligands used for modeling were carefully selected from published sources and consist of structurally diverse, conformationally rigid full agonists as active ligands together with structurally related inactives. The robustness of the models in discriminating actives from inactives was tested against four ensembles of conformations generated by using different established methods and different force fields. The reasons for the selectivity can be attributed to both geometrical differences in the arrangement of the features, e.g., different tilt angels of the π system, as well as shape differences covered by the different sets of excluded volumes. This work provides useful information for the design of new D₁ and D₂ agonists and also for comparative homology modeling of D₁ and D₂ receptors. The approach is general and could therefore be applied to other ligand–protein interactions for which no experimental protein structure is available.