Strategies To Design Selective Histone Deacetylase Inhibitors

This review classifies drug-design strategies successfully implemented in the development of histone deacetylase (HDAC) inhibitors, which have many applications including cancer treatment. Our focus is on especially demanded selective HDAC inhibitors and their structure-activity relationships in relation to corresponding protein structures. The main part of the paper is divided into six subsections each narrating how optimization of one of six structural features can influence inhibitor selectivity. It starts with the impact of the zinc binding group on selectivity, continues with the optimization of the linker placed in the substrate binding tunnel as well as the adjustment of the cap group interacting with the surface of the protein, and ends with the addition of groups targeting class-specific sub-pockets: the side-pocket-, lower-pocket- and foot-pocket-targeting groups. The review is rounded off with a conclusion and an outlook on the future of HDAC inhibitor design.     

Molecular Dynamics-Derived Pharmacophore Model Explaining the Nonselective Aspect of KV10.1 Pore Blockers

The K_V10.1 voltage-gated potassium channel is highly expressed in 70% of tumors, and thus represents a promising target for anticancer drug discovery. However, only a few ligands are known to inhibit K_V10.1, and almost all also inhibit the very similar cardiac hERG channel, which can lead to undesirable side-effects. In the absence of the structure of the K_V10.1–inhibitor complex, there remains the need for new strategies to identify selective K_V10.1 inhibitors and to understand the binding modes of the known K_V10.1 inhibitors. To investigate these binding modes in the central cavity of K_V10.1, a unique approach was used that allows derivation and analysis of ligand–protein interactions from molecular dynamics trajectories through pharmacophore modeling. The final molecular dynamics-derived structure-based pharmacophore model for the simulated K_V10.1–ligand complexes describes the necessary pharmacophore features for K_V10.1 inhibition and is highly similar to the previously reported ligand-based hERG pharmacophore model used to explain the nonselectivity of K_V10.1 pore blockers. Moreover, analysis of the molecular dynamics trajectories revealed disruption of the π–π network of aromatic residues F359, Y464, and F468 of K_V10.1, which has been reported to be important for binding of various ligands for both K_V10.1 and hERG channels. These data indicate that targeting the K_V10.1 channel pore is also likely to result in undesired hERG inhibition, and other potential binding sites should be explored to develop true K_V10.1-selective inhibitors as new anticancer agents.     

Discovery and Pharmacophore Mapping of a Low-Nanomolar Inhibitor of P. falciparum Growth

The treatment of malaria, the most common parasitic disease worldwide and the third deadliest infection after HIV and tuberculosis, is currently compromised by the dramatic increase and diffusion of drug resistance among the various species of Plasmodium, especially P. falciparum (Pf). In this view, the development of new antiplasmodial agents that are able to act via innovative mechanisms of action, is crucial to ensure efficacious antimalarial treatments. In one of our previous communications, we described a novel class of compounds endowed with high antiplasmodial activity, characterized by a pharmacophore never described before as antiplasmodial and identified by their 4,4’-oxybisbenzoyl amide cores. Here, through a detailed structure-activity relationship (SAR) study, we thoroughly investigated the chemical features of the reported scaffolds and successfully built a novel antiplasmodial agent active on both chloroquine (CQ)-sensitive and CQ-resistant Pf strains in the low nanomolar range, without displaying cross-resistance. Moreover, we conducted an in silico pharmacophore mapping.     

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.     

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.     

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.     

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

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

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.     

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

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

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

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