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

Towards a partial order graph for interactive pharmacophore exploration: extraction of pharmacophores activity delta

This paper presents a novel approach called Pharmacophore Activity Delta for extracting outstanding pharmacophores from a chemogenomic dataset, with a specific focus on a kinase target known as BCR-ABL. The method involves constructing a Hasse diagram, referred to as the pharmacophore network, by utilizing the subgraph partial order as an initial step, leading to the identification of pharmacophores for further evaluation. A pharmacophore is classified as a 'Pharmacophore Activity Delta' if its capability to effectively discriminate between active vs inactive molecules significantly deviates (by at least δ standard deviations) from the mean capability of its related pharmacophores. Among the 1479 molecules associated to BCR-ABL binding data, 130 Pharmacophore Activity Delta were identified. The pharmacophore network reveals distinct regions associated with active and inactive molecules. The study includes a discussion on representative key areas linked to different pharmacophores, emphasizing structure-activity relationships.

In silico modeling studies of N9-substituted harmine derivatives as potential anticancer agents: combination of ligand-based and structure-based approaches

A computational study was carried out to develop quantitative-structure activity relationship (QSAR), pharmacophore, molecular docking and molecular dynamics simulations of a series of N9-substituted harmine derivatives in order to investigate the structural factors involved in the cytotoxic activity and thus design new active derivatives. A valid 3 D-QSAR (R²= 0.89, q²=0.67, R²pred = 0.72) and 2 D-QSAR (R²= 0.81, q²=0.69, R²pred = 0.76) models were obtained correlating the cytotoxic activity with hydrophobic and hydrogen bond acceptor (HBA) features for 3 D-QSAR and SlogP and a_acc descriptors for 2 D-QSAR. Analysis of the selected descriptors for both models highlighted that lipophilicity and hydrogen bonding acceptor atoms remain the crucial properties and those on which cytotoxic activity depends. Also, these findings are in agreement with the characteristics of the generated pharmacophore. Furthermore, molecular docking revealed that the binding energy (-9.74 kcal/mol) and inhibition constant (0.071 µmol) correlate with the activity of the most active compound that forms hydrophobic interactions and two hydrogen bonds with the the dual specificity tyrosine phosphorylation regulated kinase 1 A (DYRK1A). The molecular dynamics simulations revealed that the protein-ligand equilibrium is stable after 100000 fs of trajectories. Based on these results, we designed new N⁹-substituted harmine derivatives with improved properties: predicted cytotoxic activities, estimated binding energies, estimated inhibition constants and interaction modes with amino acid residues of DYRK1A, compared to the best compound in the studied dataset. Additionally, these newly designed inhibitors showed promising results in the preliminary in silico Absorption, Distribution, Metabolism, Excretion and Toxicity (ADMET) evaluations.Communicated by Ramaswamy H. Sarma.

Virtual Screening Identifies Chebulagic Acid as an Inhibitor of the M2(S31N) Viral Ion Channel and Influenza A Virus

The increasing prevalence of drug-resistant influenza viruses emphasizes the need for new antiviral countermeasures. The M2 protein of influenza A is a proton-gated, proton-selective ion channel, which is essential for influenza replication and an established antiviral target. However, all currently circulating influenza A virus strains are now resistant to licensed M2-targeting adamantane drugs, primarily due to the widespread prevalence of an M2 variant encoding a serine to asparagine 31 mutation (S31N). To identify new chemical leads that may target M2(S31N), we performed a virtual screen of molecules from two natural product libraries and identified chebulagic acid as a candidate M2(S31N) inhibitor and influenza antiviral. Chebulagic acid selectively restores growth of M2(S31N)-expressing yeast. Molecular modeling also suggests that chebulagic acid hydrolysis fragments preferentially interact with the highly-conserved histidine residue within the pore of M2(S31N) but not adamantane-sensitive M2(S31). In contrast, chebulagic acid inhibits in vitro influenza A replication regardless of M2 sequence, suggesting that it also acts on other influenza targets. Taken together, results implicate chebulagic acid and/or its hydrolysis fragments as new chemical leads for M2(S31N) and influenza-directed antiviral development.

Screening of the Pan-African natural product library identifies ixoratannin A-2 and boldine as novel HIV-1 inhibitors

The continued burden of HIV in resource-limited regions such as parts of sub-Saharan Africa, combined with adverse effects and potential risks of resistance to existing antiretroviral therapies, emphasize the need to identify new HIV inhibitors. Here we performed a virtual screen of molecules from the pan-African Natural Product Library, the largest collection of medicinal plant-derived pure compounds on the African continent. We identified eight molecules with structural similarity to reported interactors of Vpu, an HIV-1 accessory protein with reported ion channel activity. Using in vitro HIV-1 replication assays with a CD4+ T cell line and peripheral blood mononuclear cells, we confirmed antiviral activity and minimal cytotoxicity for two compounds, ixoratannin A-2 and boldine. Notably, ixoratannin A-2 retained inhibitory activity against recombinant HIV-1 strains encoding patient-derived mutations that confer resistance to protease, non-nucleoside reverse transcriptase, or integrase inhibitors. Moreover, ixoratannin A-2 was less effective at inhibiting replication of HIV-1 lacking Vpu, supporting this protein as a possible direct or indirect target. In contrast, boldine was less effective against a protease inhibitor-resistant HIV-1 strain. Both ixoratannin A-2 and boldine also inhibited in vitro replication of hepatitis C virus (HCV). However, BIT-225, a previously-reported Vpu inhibitor, demonstrated antiviral activity but also cytotoxicity in HIV-1 and HCV replication assays. Our work identifies pure compounds derived from African plants with potential novel activities against viruses that disproportionately afflict resource-limited regions of the world.

Modulation of exteroceptive suppression periods in human jaw-closing muscles by subanaesthetic dose of propofol

Exteroceptive suppression (ES) periods in human jaw-closing muscles can be conditioned by a wide range of somatosensory stimuli and cognitive states. The aim of this study was to examine the effects of subanaesthetic doses of midazolam, ketamine and propofol on the short latency (ES1) and long latency (ES2) reflex in the jaw-closing muscles. First, we tried to evaluate the various methodological criteria for ES recording. We then examined the effect of subanaesthetic doses of midazolam (0·035 mg kg(-1)), ketamine (0·30 mg kg(-1)) and propofol (0·35 mg kg(-1)) on these reflexes of recording left masseter and temporalis muscle. ES duration did not differ greatly in the present study, recorded with the correct adjustment of stimulating and recording conditions. None of the subanaesthetic doses of the agents influenced ES1, and no significant effects on ES2 were observed with midazolam and ketamine. However, significant inhibitory change was observed in ES2 with propofol. ES2 is thought to be mediated by afferents, which descend in the spinal trigeminal tract and connect with a polysynaptic chain of excitatory interneurones located in the lateral reticular formation. Our observations indicate that propofol is uniquely effective not only through involvement of the gamma-aminobutyric acid type A receptor, but also through a range of other effects.

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.

ErG: 2D Pharmacophore Descriptions for Scaffold Hopping

An extended reduced graph approach (ErG) is presented that uses pharmacophore-type node descriptions to encode the relevant molecular properties. The basic idea of the method can be described as a hybrid approach of reduced graphs (Gillet et al. J. Chem. Inf. Comput. Sci. 2003, 43, 338-345) and binding property pairs (Kearsley et al. J. Chem. Inf. Comput. Sci. 1996, 36, 118-127). However, specific extension modifications to correctly describe the pharmacophoric properties, size, and shape of the molecules under study result in a very stable and good performance as compared to DAYLIGHT fingerprints (DFP). This is exemplified for 11 activity classes of the MDL Drug Data Report database, for which ErG performs as well or better than DFP in 10 cases. On the basis of the example data sets, the ability of ErG to switch from one chemotype to another (often referred to as "scaffold hopping") is highlighted. Additionally, possible pitfalls of reduced graph approaches as well as suitable solutions are discussed with the help of example structures. Overall, it is shown that ErG is a widely applicable method capable of identifying structurally diverse actives for a given active search query. This diversity is achieved by a high degree of molecular abstraction, which in turn results in a low dimensional descriptor vector that allows very low computation times for similarity searches.

Ligand-Based Molecular Modeling Study on a Chemically Diverse Series of Cholecystokinin-B/Gastrin Receptor Antagonists: Generation of Predictive Model

Pharmacophore hypotheses were developed for six structurally diverse series of cholecystokinin-B/gastrin receptor (CCK-BR) antagonists. A training set consisting of 33 compounds was carefully selected. The activity spread of the training set molecules was from 0.1 to 2100 nM. The most predictive pharmacophore model (hypothesis 1), consisting of four features, namely, two hydrogen bond donors, one hydrophobic aliphatic, and one hydrophobic aromatic feature, had a correlation (r) of 0.884 and a root-mean-square deviation of 1.1526, and the cost difference between null cost and fixed cost was 81.5 bits. The model was validated on a test set consisting of six different series of 27 structurally diverse compounds and performed well in classifying active and inactive molecules correctly. This validation approach provides confidence in the utility of the predictive pharmacophore model developed in this work as a 3D query tool in the virtual screening of drug-like molecules to retrieve new chemical entities as potent CCK-BR antagonists. The model can also be used to predict the biological activities of compounds prior to their costly and time-consuming synthesis.

Molecular design based on 3D pharmacophores. Applications to 5-HT7 receptors

A definition of a pharmacophore for the 5-HT7 antagonists was carried out by searching the common chemical features of selective antagonists from the literature. A molecular design is described by analyzing the differences between this new pharmacophore and three other 3D serotonin pharmacophores previously described. This comparison led to the synthesis of a new series of potent 5-HT7 antagonists.

Comparison of Conformational Analysis Techniques To Generate Pharmacophore Hypotheses Using Catalyst

Generation of reliable pharmacophore models is a key strategy in drug design. The quality of a pharmacophore model is known to depend on several factors, with the quality of the conformer sets used perhaps being one of the most important. The goal of this study was to compare different conformational analysis methods to determine if one was superior to the others for pharmacophore generation using Catalyst/HypoGen. The five methods selected were Catalyst/Fast, Catalyst/Best, Omega, Chem-X and MacroModel. Data sets for which Catalysts models had previously been published were selected using defined quality measures. Hypotheses were generated for each of the data sets and the performance of the different conformational analysis methods was compared using both quantitative (cost and correlation coefficients) and qualitative measures (by comparing the hypotheses in terms of the features present and their spatial relationships). Two main conclusions emerged from the study. First, it was not always possible to replicate the literature results. The reasons for these failures are explored in detail, and a template for use in publications that apply the Catalyst methodology is proposed. Second, the faster rule-based methods for conformational analysis give pharmacophore models that are just as good as, and in some cases better than, the models generated using the slower, more rigorous approaches.

POT-DMC: A Virtual Screening Method for the Identification of Potent Hits

A method for ligand-based virtual screening (LBVS), dynamic mapping of consensus positions (DMC), has been extended to take different potency levels of template compounds into account. This potency scaling technique is designed to tune search calculations toward the detection of increasingly potent hits. LBVS analysis of three different compound classes confirmed the ability of potency-scaled DMC (POT-DMC) to identify active database compounds with higher potency than conventional calculations.

Optimization of the pharmacophore model for 5-HT7R antagonism. Design and synthesis of new naphtholactam and naphthosultam derivatives

We present in this study an optimization of a preliminary pharmacophore model for 5-HT(7)R antagonism, with the incorporation of recently reported ligands and using an efficient procedure with the CATALYST program. The model consists of five features: a positive ionizable atom (PI), a H-bonding acceptor group (HBA), and three hydrophobic regions (HYD). This model has been supported by the design, synthesis, and biological evaluation of new naphtholactam and naphthosultam derivatives of general structure I (39-72). A systematic structure-affinity relationship (SAFIR) study on these analogues has allowed us to confirm that the model incorporates the essential structural features for 5-HT(7)R antagonism. In addition, computational simulation of the complex between compound 56 and a rhodopsin-based 3D model of the 5-HT(7)R transmembrane domain has permitted us to define the molecular details of the ligand-receptor interaction and gives additional support to the proposed pharmacophore model for 5-HT(7)R antagonism: (i) the HBA feature of the pharmacophore model binds Ser(5.42) and Thr(5.43), (ii) the HYD1 feature interacts with Phe(6.52), (iii) the PI feature forms an ionic interaction with Asp(3.32), and (iv) the HYD3 (AR) feature interacts with a set of aromatic residues (Phe(3.28), Tyr(7.43)). These results provide the tools for the design and synthesis of new ligands with predetermined affinities and pharmacological properties.

Arylguanidine and arylbiguanide binding at 5-HT3 serotonin receptors: a QSAR study

For a series of monosubstituted arylguanidines, 5-HT3 receptor affinity was found generally related to the electron withdrawing nature of the substituent at the aryl 3-position and the lipophilicity of the 4-position substituent. A broader examination of 35 arylguanidines and arylbiguanides revealed that affinity could be described by molecular polarizability, a Chi index term (8chiP), and the sum of all (-Cl) E-State values (SsCl) in the molecule.

Synthesis and pharmacological characterization of a new benzoxazole derivative as a potent 5-HT3 receptor agonist

N-(2-Benzoxazol-2-yl-ethyl)-guanidine hydrochloride (10) was synthesized and pharmacologically tested. This compound showed high affinity for the 5-HT(3) receptor (K(i)=0.77 nM) and potently triggered the von Bezold-Jarisch reflex (BJR) in rats with an ED(50)=0.52 microg/kg iv and intrinsic activity next to 1 (i.a.=0.94). This stimulant effect was abolished by pretreatment with the 5-HT(3) receptor antagonist granisetron and was subject to a rapid and pronounced tachyphylaxis, due to desensitization of the peripheric cardiac 5-HT(3) receptor. Consequently, 10 acts as an in vivo 5-HT(3) antagonist inhibiting the BJR responses evoked by submaximal doses of 5-HT with an ID(50)=5.8 microg/kg iv.

Prediction of 5-HT3 receptor agonist-binding residues using homology modeling

5-HT(3) receptors demonstrate significant structural and functional homology to other members of the Cys-loop ligand-gated ion channel superfamily. The extracellular domains of these receptors share similar sequence homology (approximately 20%) with Limnaea acetylcholine binding protein, for which an x-ray crystal structure is available. We used this structure as a template for computer-based homology modeling of the 5-HT(3) receptor extracellular domain. AutoDock software was used to dock 5-HT into the putative 5-HT(3) receptor ligand-binding site, resulting in seven alternative energetically favorable models. Residues located no more than 5 A from the docked 5-HT were identified for each model; of these, 12 were found to be common to all seven models with five others present in only certain models. Some docking models reflected the cation-pi interaction previously demonstrated for W183, and data from these and other studies were used to define our preferred models.

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.

Functional group interactions of a 5-HT3R antagonist

BACKGROUND

Lerisetron, a competitive serotonin type 3 receptor (5-HT3R) antagonist, contains five functional groups capable of interacting with amino acids in the 5-HT3R binding site. Site directed mutagenesis studies of the 5-HT3AR have revealed several amino acids that are thought to form part of the binding domain of this receptor. The specific functional groups on the ligand that interact with these amino acids are, however, unknown. Using synthetic analogs of lerisetron as molecular probes in combination with site directed mutagenesis, we have identified some of these interactions and have proposed a model of the lerisetron binding site.

RESULTS

Two analogs of lerisetron were synthesized to probe 5-HT3R functional group interactions with this compound. Analog 1 lacks the N1 benzyl group of lerisetron and analog 2 contains oxygen in place of the distal piperazine nitrogen. Both analogs show significantly decreased binding affinity to wildtype 5-HT3ASRs. Mutations at W89, R91, Y142 and Y152 produced significant decreases in binding compared to wildtype receptors. Binding affinities of analogs 1 and 2 were altered only by mutations at W89, and Y152.

CONCLUSIONS

Based on the data obtained for lerisetron and analogs 1 and 2, we have proposed a tentative model of the lerisetron binding pocket of the 5-HT3ASR. According to this model, The N-benzyl group interacts in a weak interaction with R91 while the benzimidazole group interacts with W89. Our data support an interaction of the distal amino nitrogen with Y142 and Y152.

Molecular design based on 3D-pharmacophore. Application to 5-HT subtypes receptors

A first definition of a pharmacophore for the serotonin reuptake inhibitors was carried out by considering a three-dimensional model which correlates the chemical structures of series of reuptake inhibitors with their biological affinities. A molecular design was described by analyzing two different 3D serotonin pharmacophores. This successful approach enabled us to consider the design of new serotonin ligands by the same method.

Pharmacophore mapping of a series of 2,4-diamino-5-deazapteridine inhibitors of Mycobacterium avium complex dihydrofolate reductase

Pharmacophore hypotheses were developed for a series of 2,4-diamino-5-deazapteridine inhibitors of Mycobacterium avium complex (MAC) and human dihydrofolate reductase (hDHFR). Training sets consisting of 20 inhibitors were selected in each case on the basis of the information content of the structures and activity data as required by the HypoGen program in the Catalyst software. In the case of MAC DHFR inhibitors, the best pharmacophore in terms of statistics and predictive value consisted of four features: two hydrogen bond acceptors (HA), one hydrophobic (HY) feature, and one ring aromatic (RA) feature. The selected pharmacophore hypothesis yielded an rms deviation of 0.730 and a correlation coefficient of 0.967 with a cost difference (null cost minus total cost) of approximately 52. The pharmacophore was validated on a large set of test inhibitors. For the test series, a classification scheme was used to distinguish highly active from moderately active and inactive compounds on the basis of activity ranges. This classification scheme is more practical than actual estimated values because these values have no meaning for compounds yet to be tested except that they indicate whether the compounds will be active or inactive in a biological assay. For the training set, the success rate for predicting active and inactive compounds was 100%. For the test set, the success rate in predicting active compounds was greater than 92% while about 7% of the inactive compounds were predicted to be active. This successful prediction was further validated on three structurally diverse compounds active against MAC DHFR. Two compounds mapped well onto three of the four features of the pharmacophore. The third compound was mapped to all four features of the pharmacophore. This validation study provided confidence for the usefulness of the selected pharmacophore model to identify compounds with diverse structures from a database search. Comparison of pharmacophores for inhibitors of human and MAC DHFR is expected to reveal fundamental differences between these two pharmacophores that may be effectively exploited to identify and design compounds with high selectivity for MAC DHFR.

5-HT(3)R binding of lerisetron: an interdisciplinary approach to drug-Receptor interactions

The design, synthesis, and use of lerisetron-based molecular probes to investigate the 5-HT(3)R binding site are described. A SAR study, which involved distance and electronic parameter modifications of lerisetron's N-benzyl group, resulted in the discovery of a partial agonist.

The binding of arylguanidines at 5-HT(3) serotonin receptors: a structure-affinity investigation

The 5-HT(3) receptor binding affinities of nine pairs of aryl-substituted arylguanidines and arylbiguanides were examined and the results suggest the likelihood that both classes of agents utilize common receptor binding features. The effects of structural modification were also examined using CoMFA. 1-(3,4,5-Trichlorophenyl)guanidine (5-HT(3) K(i)=0.7 nM) was identified as a very high-affinity arylguanidine. The structures of the high-affinity arylguanidines are inconsistent with current 5-HT(3) pharmacophore models.

Association of two 3D QSAR analyses. application to the study of partial agonist serotonin-3 ligands

CATALYST and COMFA, two software packages for 3D QSAR studies, were associated to correlate the three-dimensional structures of 75 serotonin 5-HT3 ligands to their biological affinities. The conformational analysis and the influence of chemical function-based alignments (the basis of this association) on final results are discussed in this publication. These two analyses allow for precisely quantitating the weights of significant chemical groups or functions on the biological affinities.

A feature based pharmacophore for Candida albicans MyristoylCoA: protein N-myristoyltransferase inhibitors

A three-dimensional pharmacophore model has been generated for Candida albicans MyristoylCoA: protein N-myristoyltransferase (NMT) inhibitors, using the software program CATALYST. The in vitro NMT inhibitory activity of a series of peptidic inhibitors was used for pharmacophore generation. The effect of altering the control parameters and feature selection was studied to arrive at the pharmacophore model. The selection of the best hypothesis model was based on the total cost, predictive ability, difference in the cost from the null hypothesis and alignment of the training set compounds on to the hypothesis. The pharmacophore model selected has four features; one hydrophobic, two hydrogen bond acceptor and one positive ionisable function. Groups identified as necessary by scanning alanine mutagenesis studies of the peptidic substrate of C. albicans NMT, have been identified as pharmacophore features. Comparison of the ligand binding with the enzyme in the crystal structure of NMT and that proposed by the phamacophore is consistent. The pharmacophore thus generated can be used as a template for designing non-peptidic inhibitors of NMT.