Molecular design based on 3D-pharmacophore. Application to 5-HT4 receptor

A definition of a pharmacophore for the 5-HT4 antagonist was carried out by considering a three-dimensional model which correlates the chemical structures of series of antagonists with their biological affinities. A molecular design is described by analyzing the differences between two 3D serotonin pharmacophores. This successful structural modification demonstrates the efficiency of this approach to design new serotonin ligands.     

New serotonin 5-HT(6) ligands from common feature pharmacophore hypotheses

Serotonin 5-HT6 receptor antagonists are thought to play an important role in the treatment of psychiatry, Alzheimer's disease, and probably obesity. To find novel and potent 5-HT6 antagonists and to provide a new idea for drug design, we used a ligand-based pharmacophore to perform the virtual screening of a commercially available database. A three-dimensional common feature pharmacophore model was developed by using the HipHop program provided in Catalyst software and was used as a query for screening the database. A recursive partitioning (RP) model which can separate active and inactive compounds was used as a filtering system. Finally a sequential virtual screening procedure (SQSP) was conducted, wherein both the common feature pharmacophore and the RP model were used in succession to improve the results. Some of the hits were selected based on druglikeness, ADME properties, structural diversity, and synthetic accessibility for real biological evaluation. The best hit compound showed a significant IC50 value of 9.6 nM and can be used as a lead for further drug development.     

5-HT3受体拮抗剂药效团模型的构建

以31个来源于MDDR数据库中具有抑制鼠Bezold-Jarisch反射作用的5-HT3受体拮抗剂作为训练集化合物, 构建5-HT3受体拮抗剂药效团模型. 训练集化合物具备结构多样性, 来源于相同药理模型, 活性值ED50范围为0.05~320 μg/kg i.v.. 利用Catalyst计算5-HT3受体拮抗剂的最优药效团由一个氢键受体、一个疏水基团、一个正电离子化基团、一个芳香环特征和6个排除体积组成; Fixed cost值、Null cost 值、Δcost值和Configuration cost值分别为112.6, 172.0, 59.4和7.248. 训练集化合物活性的计算值与实测值相关系数为0.9031, 偏差值为0.8976, 基于Fischer的交叉验证结果表明药效团模型具有较高的置信度, 所得药效团对训练集化合物活性值的预测结果显示有较好的预测能力, 可用于数据库搜索指导发现新的具有该活性的先导化合物, 也可用于中药或天然产物药物研究开发.     

Molecular modeling studies focused on 5-HT7 versus 5-HT1A selectivity. Discovery of novel phenylpyrrole derivatives with high affinity for 5-HT7 receptors

The present study discusses the well-known 5-HT7/5-HT1A selectivity issue through a new series of phenylpyrrole derivatives. The first hits emerged from a virtual screening performed on a chemolibrary. Further study led to an optimization of a preliminary 5-HT7 pharmacophore model. The importance of each pharmacophoric feature is confirmed, but these characteristics have to be coupled to geometric constraints in order to achieve a 5-HT7 selectivity. Indeed, 5-HT1A affinity probably arises from extended conformations, whereas a bent one appears to be best suited for 5-HT7 selectivity.     

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.     

人类腺苷受体A3亚型拮抗剂的构效关系分析

构建人类腺苷受体A3亚型药效团模型和三维蛋白结构模型用于作用模式研究.以18个来源于文献具有腺苷受体A3亚型拮抗活性的化合物作为训练集,使用HypoGen方法构建药效团模型.通过同源模建和分子动力学模拟构建了人类腺苷受体A3亚型的三维蛋白模型,并利用PROCHECK方法评估该模型的合理性,对所得的结构使用分子对接程序进行作用模式分析,药效团模型和同源模建结果相互匹配较好.使用新药效团模型对MDL药物数据库(MDDR)中包含的约120000个化合物进行虚拟筛选,得到了8个候选化合物,用于进一步的生物学评价和活性测定.本工作对于人类腺苷受体A3亚型拮抗剂的设计和抗哮喘药物的研发具有一定的理论指导和应用价值.     

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重摄取抑制剂提供依据.     

On the value of homology models for virtual screening: discovering hCXCR3 antagonists by pharmacophore-based and structure-based approaches

Human chemokine receptor CXCR3 (hCXCR3) antagonists have potential therapeutic applications as antivirus, antitumor, and anti-inflammatory agents. A novel virtual screening protocol, which combines pharmacophore-based and structure-based approaches, was proposed. A three-dimensional QSAR pharmacophore model and a structure-based docking model were built to virtually screen for hCXCR3 antagonists. The hCXCR3 antagonist binding site was constructed by homology modeling and molecular dynamics (MD) simulation. By combining the structure-based and ligand-based screenings results, 95% of the compounds satisfied either pharmacophore or docking score criteria and would be chosen as hits if the union of the two searches was taken. The false negative rates were 15% for the pharmacophore model, 14% for the homology model, and 5% for the combined model. Therefore, the consistency of the pharmacophore model and the structural binding model is 219/273 = 80%. The hit rate for the virtual screening protocol is 273/286 = 95%. This work demonstrated that the quality of both the pharmacophore model and homology model can be measured by the consistency of the two models, and the false negatives in virtual screening can be reduced by combining two virtual screening approaches.     

Benzimidazole derivatives. 4. The recognition of the voluminous substituent attached to the basic amino group of 5-HT4 receptor antagonists

We report a structure-affinity analysis of an important element in the pharmacophore model for the recognition of 5-HT4 receptor antagonists: the voluminous substituent attached to the basic nitrogen of the ligand. We have designed, synthesized and pharmacologically evaluated a series of benzimidazole derivatives 1 containing a common molecular skeleton formed by N-[(4-piperidyl)methyl]-6-chlorobenzimidazole-4-carboxamide and four different substituents (R = butyl, 2-[(methylsulfonyl)amino]ethyl, 5-[(phenylacetyl)amino]pentyl, and 5-[(benzylsulfonyl)amino]pentyl). These compounds possess binding affinities in the nM range (Ki = 0.11-1.50 nM). Moreover, a ligand that contains a hydrogen atom attached to the basic nitrogen (R = H; Ki = 150 nM) is used as a control for structure-affinity relationships.     

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.     

α1A-亚型肾上腺素受体拮抗剂3D药效团模型的研究

运用Catalyst软件以34个α_1A-AR拮抗剂分子为训练集, 构建了包含一个氢键受体、一个正电中心和一个芳环中心的三元素药效团模型, 线性回归相关系数为0.89. 经13个分子组成的测试集验证该药效团模型具有较好的活性预测能力, 为寻找新的α_1A-AR拮抗剂分子提供了理论基础.     

基于吡咯烷与正丁烷类衍生物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.     

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.     

新型5-HT重摄取抑制剂的设计、合成及活性评价

在对已知各种结构类型的5-HT重摄取抑制剂分子结构全面分析的基础上, 建立了SSRIs药效团模型. 基于该模型应用UNITY程序对NCI-3D和Maybridge-3D数据库进行三维结构的限制性查询, 在获得的命中结构的信息指导下, 设计合成了3种全新结构类型的化合物, 并完成了初步的药理活性评价. 这些化合物均显示出不同程度的5-HT重摄取抑制活性, 其中5个化合物显示高抑制活性. 哌嗪取代的二苯脒类化合物的结构新颖, 较好地符合5-HT重摄取抑制剂药效团模型, 与SSRIs类化合物三维定量构效关系研究得到的CoMFA模型有较好的适配性.     

An improved synthesis of piperazino-piperidine based CCR5 antagonists with flexible variation on pharmacophore sites

An improved and efficient synthetic route towards piperidino-piperazine based CCR5 antagonists was developed. The new approach was flexible for introducing various substituents in the pharmacophore sites via Grignard reagent addition and reductive amination. l-Amino acids were used as a chiral pool to introduce and then induce the desired stereochemistries, meanwhile rendering the variable substitution. The efficient construction of the piperazino-piperidine nucleus was achieved in a highly convergent manner with a key building block of N¹-Boc-4-substituent-4-aminopiperidine, exhibiting significant advantages in terms of concise synthetic route and environmental-friendly reagents over the previously described stepwise synthesis, in which a modified Strecker reaction was involved with highly toxic reagents such as diethylaluminum cyanide.     

Definition of a pharmacophore for partial agonists of serotonin 5-HT3 receptors

A definition of a partial agonists serotonin 5-HT3 pharmacophore was carried out by considering a three-dimensional model which correlates the chemical structures of series of piperazinopyrrolothienopyrazines, piperazinopyridopyrrolopyrazines, piperazinopyrroloquinolaxines, piperazinopyridopyrroloquinoxalines, aminoalkyloximinopyrroloindoles, aminoalkyloximinothienopyrrolizines, and aminoalkyloximinopyrrolizines with the biological affinities. The model is formed by five features corresponding to two hydrogen bond acceptors, one aromatic ring, one hydrophobic group, and one positive ionizable site (quaternary ammonium ions). The nature of the features and the distances between them explain the partial agonist activities of these compounds.     

Structure-activity relationship studies on CXCR4 antagonists having cyclic pentapeptide scaffolds

Structure-activity relationship studies on CXCR4 antagonists, which were previously found by using cyclic pentapeptide libraries, were performed to optimize side-chain functional groups, involving conformationally constrained analogues. In addition, a new lead of cyclic pentapeptides with the introduction of a novel pharmacophore was developed.     

A first step towards the understanding of the 5-HT(3) receptor subunit heterogeneity from a computational point of view

The functional serotonin type-3 receptor (5-HT(3)-R), which is the target of many neuroactive drugs, is known to be a homopentamer made of five identical subunits A (5-HT(3A)-R) or a binary heteropentamer made of subunits A and B (5-HT(3A/B)-R) with a still debated arrangement and stoichiometry. This complex picture has been recently further complicated by the discovery of additional 5-HT(3)-R subunits, C, D, and E, which, similarly to the B subunit, are apparently able to form functional receptors only if co-expressed with subunit A. Being the binding site for both serotonin and antagonists (i.e. drugs) located at the extracellular interface between two adjacent subunits, the large variability of the 5-HT(3)-R composition becomes a crucial issue, since it can originate many different interfaces providing non-equivalent ligand binding sites and complicating the pharmacological modulation. Here, the different 5-HT(3)-R interfaces are analysed, on the bases of the structural conformations of previously built 3D homology models and of the known subunit sequences, by addressing their physicochemical characterization. The results confirm the presence of an aromatic cluster located in the core of the A-A interface as a key determinant for having an interface both stable and functional. This is used as a discriminant to make hypotheses about the capability of all the other possible interfaces constituted by the known 5-HT(3)-R sequences A, B, C, D, and E to build active receptors.