Benzodiazepine-induced hyperphagia: development and assessment of a 3D pharmacophore by computational methods

Benzodiazepine receptor (BDZR) ligands are structurally diverse compounds that bind to specific binding sites on GABA(A) receptors and allosterically modulate the effect of GABA on chloride ion flux. The binding of BDZR ligands to this receptor system results in activity at multiple behavioral endpoints, including anxiolytic, sedative, anticonvulsant, and hyperphagic effects. In the work presented here, a computational procedure developed in our laboratory has been used to obtain a 3D pharmacophore for ligand recognition of the GABA(A)/BDZRs initiating the hyperphagic response. To accomplish this goal, 17 structurally diverse compounds, previously assessed in our laboratory for activity at the hyperphagic endpoint, were used. The result is a four-component 3D pharmacophore. It consists of two proton acceptor atoms, the centroid of an aromatic ring and the centroid of a hydrophobic moiety in a common geometric arrangement in all compounds with activity at this endpoint. This 3D pharmacophore was then assessed and successfully validated using three different tests. First, two BDZR ligands, which were included as negative controls in the set of seventeen compounds used for the pharmacophore development, did not fit the pharmacophore. Second, some benzodiazepine ligands known to have activity at the hyperphagia endpoint, but not included in the pharmacophore development, were used as positive controls and were found to fit the pharmacophore. Finally, using the 3D pharmacophore developed in the present work to search 3D databases, over 50 classical benzodiazepines were found. Among them, were benzodiazepine ligands known to have an effect at the hyperphagic endpoint. In addition, the novel compounds also found in this search are promising therapeutic agents that could beneficially affect feeding behavior.     

Pharmacophore mapping and in silico screening to identify new potent leads for A(2A) adenosine receptor as antagonists

A(2A) adenosine receptor (AR) antagonists play an important role in neurodegenerative diseases like Parkinson's disease. A 3D-QSAR study of A(2A) AR antagonists, was taken up to design best pharmacophore model. The pharmacophoric features (ADHRR) containing a hydrogen bond acceptor (A), a hydrogen bond donor (D), a hydrophobic group (H) and two aromatic rings (R), is projected as the best predictive pharmacophore model. The QSAR model was further treated as a template for in silico search of databases to identify new scaffolds. The binding patterns of the leads with A(2A) AR are analysed using docking studies and novel potent ligands of A(2A) AR are projected.     

Benzodiazepine-Induced Hyperphagia: Development and Assessment of a 3D Pharmacophore By Computational Methods

Abstract

Benzodiazepine receptor (BDZR) ligands are structurally diverse compounds that bind to specific binding sites on GABAA receptors and allosterically modulate the effect of GABA on chloride ion flux. The binding of BDZR ligands to this receptor system results in activity at multiple behavioral endpoints, including anxiolytic, sedative, anticonvulsant, and hyperphagic effects. In the work presented here, a computational procedure developed in our laboratory has been used to obtain a 3D pharmacophore for ligand recognition of the GABAA/BDZRS initiating the hyperphagic response. To accomplish this goal, 17 structurally diverse compounds, previously assessed in our laboratory for activity at the hyperphagic endpoint, were used. The result is a four-component 3D pharmacophore. It consists of two proton acceptor atoms, the centroid of an aromatic ring and the centroid of a hydrophobic moiety in a common geometric arrangement in all compounds with activity at this endpoint. This 3D pharmacophore was then assessed and successfully validated using three different tests. First, two BDZR ligands, which were included as negative controls in the set of seventeen compounds used for the pharmacophore development, did not fit the pharmacophore. Second, some benzodiazepine ligands known to have activity at the hyperphagia endpoint, but not included in the pharmacophore development, were used as positive controls and were found to fit the pharmacophore. Finally, using the 3D pharmacophore developed in the present work to search 3D databases, over 50 classical benzodiazepines were found. Among them, were benzodiazepine ligands known to have an effect at the hyperphagic endpoint. In addition, the novel compounds also found in this search are promising therapeutic agents that could beneficially affect feeding behavior.     

A three-dimensional model of the delta-opioid pharmacophore: comparative molecular modeling of peptide and nonpeptide ligands

A comparative molecular modeling study of delta-opioid ligands was performed under the assumption that potent peptide and nonpeptide agonists may have common three-dimensional (3D) arrangement of pharmacophore groups upon binding to the delta-receptor. Low-energy conformations of the agonists 7-spiroindanyloxymorphone (SIOM) and 2-methyl-4a-alpha-(3-hydroxyphenyl)-1,2,3,4,4a,5,12, 12a-alpha-octahydro-quinolino[2,3,3-g]isoquinoline (TAN-67), and a partial agonist oxomorphindole (OMI) were determined by high-temperature molecular dynamics (MD). A good spatial overlap was found for the pharmacophore groups of SIOM, TAN-67, and OMI, including the basic nitrogen, phenol hydroxyl, and two aromatic ring. Based on this overlap we proposed a 3D pharmacophore model for nonpeptide delta-opioid agonists with a distance of 7.0 +/- 1.3 A between the two aromatic rings and of 8.2 +/- 1.0 A between the nitrogen and phenyl ring. The potent and highly delta-opioid receptor selective agonist [(2S,3R)-TMT(1)]DPDPE, which shares global backbone constraints of the 14-membered disulfide cycle and a strong preference for the trans rotamer of the TMT(1) side chain, was chosen as a peptide template of the delta-opioid pharmacophore. Extensive MD simulations at 300 K with the AMBER force field were performed for [(2S,3R)-TMT(1)]DPDPE and the less potent [(2S, 3S)-TMT(1)]DPDPE analogue. Multiple MD trajectories were collected for each peptide starting from the x-ray structures of DPDPE and [L-Ala(3)]DPDPE and from models proposed in the literature. Low-energy MD conformations were filtered by the nonpeptide pharmacophore query and then directly superimposed with SIOM, OMI, and TAN-67. Two conformers of [(2S,3R)-TMT(1)]DPDPE that showed the best overlap with the nonpeptide pharmacophore (rms deviation 相似文献   

Development of a 3D pharmacophore for nonspecific ligand recognition of alpha1, alpha2, alpha3, alpha5, and alpha6 containing GABA(A)/benzodiazepine receptors

Transfected cells containing GABA(A) benzodiazepine receptors (BDZRs) have been utilized to systematically determine the affinity of ligands at alpha1, alpha2, alpha3, alpha5 and alpha6 subtypes in combination with beta2 and gamma2. All but a few of the ligands thus far studied have relatively high affinities for each of these alpha subtype receptors. Thus, these ligands must contain common stereochemical properties favorable for recognition by each of the subtype combinations. In the present work, such a common three-dimensional (3D) pharmacophore for recognition of alpha1, alpha2, alpha3, alpha5 and alpha6 containing GABA(A)/BDZRs types of receptors has been developed and assessed, using as a database receptor affinities measured in transfected cells for 27 diverse compounds. The 3D-recognition pharmacophore developed consists of three proton accepting groups, a hydrophobic group, and the centroid of an aromatic ring found in a common geometric arrangement in the 19 nonselective ligands used. Three tests were made to assess this pharmacophore: (i) Four low affinity compounds were used as negative controls, (ii) Four high affinity compounds, excluded from the pharmacophore development, were used as compounds for pharmacophore validation, (iii) The 3D pharmacophore was used to search 3D databases. The results of each of these types of assessments provided robust validation of the 3D pharmacophore. This 3D pharmacophore can now be used to discover novel nonselective ligands that could be activation selective at different behavioral end points. Additionally, it may serve as a guide in the design of more selective ligands, by determining if candidate ligands proposed for synthesis conform to this pharmacophore and selecting those that do not for further experimental assessment.     

Development and assessment of a 3D pharmacophore for ligand recognition of BDZR/GABAA receptors initiating the anxiolytic response

Benzodiazepine receptor (BDZR) ligands are structurally diverse compounds that bind to specific binding sites on GABAA receptors and allosterically modulate the effect of GABA on chloride flux. The binding of BDZR ligands to this receptor system results in activity at multiple behavioral end points including anxiolytic, sedative, hyperphagic, anticonvulsant and hyperthermic effects. In the work presented here, 17 structurally diverse BDZR ligands of the receptors initiating the anxiolytic response have been studied using a systematic computational procedure developed in our laboratory. Using this procedure, a five component 3D recognition pharmacophore was obtained consisting of two proton acceptors, a hydrophobic group, an aromatic electron accepting ring and a ring containing polar moieties, all found in a common geometric arrangement in the 15 compounds with an effect at the anxiolytic end point and absent in two control compounds. The 3D pharmacophore developed was validated by searching 3D databases and finding known BDZR ligands active at the anxiolytic end point, including 1,4-BDZ derivatives, imidazo BDZ and beta-carboline ligands.     

Pharmacophore mapping and in silico screening to identify new potent leads for A2A adenosine receptor as antagonists

A_2A adenosine receptor (AR) antagonists play an important role in neurodegenerative diseases like Parkinson’s disease. A 3D-QSAR study of A_2A AR antagonists, was taken up to design best pharmacophore model. The pharmacophoric features (ADHRR) containing a hydrogen bond acceptor (A), a hydrogen bond donor (D), a hydrophobic group (H) and two aromatic rings (R), is projected as the best predictive pharmacophore model. The QSAR model was further treated as a template for in silico search of databases to identify new scaffolds. The binding patterns of the leads with A_2A AR are analysed using docking studies and novel potent ligands of A_2A AR are projected.     

Pharmacophore modeling, 3D-QSAR,and in silico ADME prediction of N-pyridyl and pyrimidine benzamides as potent antiepileptic agents

Biological mechanism attributing mutations in KCNQ2/Q3 results in benign familial neonatal epilepsy (BFNE), a rare form of epilepsy and thus neglected. It offers a potential target for antiepileptic drug discovery. In the present work, a pharmacophore-based 3D-QSAR model was generated for a series of N-pyridyl and pyrimidine benzamides possessing KCNQ2/Q3 opening activity. The pharmacophore model generated contains one hydrogen bond donor (D), one hydrophobic (H), and two aromatic rings (R). They are the crucial molecular write-up detailing predicted binding efficacy of high affinity and low affinity ligands for KCNQ2/Q3 opening activity. Furthermore, it has been validated by using a biological correlation between pharmacophore hypothesis-based 3D-QSAR variables and functional fingerprints of openers responsible for the receptor binding and also by docking of these benzamides into the validated homology model. Excellent statistical computational tools of QSAR model such as good correlation coefficient (R^2?>^?0.80), higher F value (F?>?39), and excellent predictive power (Q²?>?0.7) with low standard deviation (SD <0.3) strongly suggest that the developed model could be used for prediction of antiepileptic activity of newer analogs. A preliminary pharmacokinetic profile of these derivatives was also performed on the basis of QikProp predictions.     

Synthetic and computer assisted analysis of the pharmacophore for agonists at benzodiazepine receptors.

In order to employ rational drug design in the discovery of selective benzodiazepine receptor agonists and inverse agonists, pharmacophore/receptor models for both these activities must first be established. Recently, a pharmacophore for the inverse agonist site has been formulated employing the most recent receptor mapping techniques (22). The continuation of this approach to the pharmacophore for agonist ligands has permitted a definition of this site independently of the inverse agonist model. The agonist pharmacophore/receptor contains two hydrogen bond donating sites of interaction (H1 and H2) located about 6.5 A from each other, as well as three areas of lipophilic interaction (L1-L3). The areas L1 and L2 are critical for agonist activity; moreover, some ligands also require an interaction in a third lipophilic area termed L3. This is in agreement with previous work (12-23). In addition, an area of negative steric interaction (S1) between the ligand and receptor-binding protein is defined. In regard to the pharmacophore, it was established that the alignment rule for agonist beta-carbolines is different from that which elicits inverse agonist activity. Consideration of the pharmacophore has resulted in the synthesis of a new beta-carboline 16 which elicits agonist activity. This ligand 16 not only satisfied the requirements of the pharmacophore, but more importantly it elicited both anticonvulsant and anxiolytic activity, but was devoid of the myorelaxant/ataxic properties associated with the benzodiazepines.     

Pharmacophore modelling and atom-based 3D-QSAR studies on N-methyl pyrimidones as HIV-1 integrase inhibitors

Pharmacophore modelling and atom-based 3D-QSAR studies were carried out for a series of compounds belonging to N-methyl pyrimidones as HIV-1 integrase inhibitors. Based on the ligand-based pharmacophore model, we got 5-point pharmacophore model AADDR, with two hydrogen bond acceptors (A), two hydrogen bond donors (D) and one aromatic ring (R). The generated pharmacophore-based alignment was used to derive a predictive atom-based 3D-QSAR model for the training set (r(2)?=?0.92, SD?=?0.16, F?=?84.8, N?=?40) and for test set (Q(2)?=?0.71, RMSE?=?0.06, Pearson R?=?0.90, N?=?10). From these results, AADDR pharmacophore feature was selected as best common pharmacophore hypothesis, and atom-based 3D-QSAR results also support the outcome by means of favourable and unfavourable regions of hydrophobic and electron-withdrawing groups for the most potent compound 30. These results can be useful for further design of new and potent HIV-1 IN inhibitors.     

3D-QSAR studies of 2,2-diphenylpropionates to aid discovery of novel potent muscarinic antagonists

Muscarinic acetylcholine receptors (mAChRs) consisting of five known subtypes, are widely distributed in both central and peripheral nervous systems for regulation of a variety of critical functions. The present theoretical study describes correlations between experimental and calculated molecular properties of 15 α-substituted 2,2-diphenylpropionate antimuscarinics using quantum chemical and pharmacophore generation methods to characterize the drug mAChR properties and design new therapeutics. The calculated stereoelectronic properties, such as total energies, bond distances, valence angles, torsion angles, HOMO–LUMO energies, reactivity indices, vibrational frequencies of ether and carbonyl moieties, and nitrogen atom proton affinity were found to be well correlated when compared with experimentally determined inhibition constants from the literature using three muscarinic receptor assays: [³H]NMS receptor binding, α-amylase release from rat pancreas, and guinea pig ileum contraction. In silico predicted toxicity on rat oral LD₅₀ values correlated well with the [³H]NMS binding in N4TG1 cells and α-amylase release assays, but not the ileum contraction assay. Next, to explore the functional requirements for potent activity of the compounds, we developed a preliminary 3D pharmacophore model using the in silico techniques. The resulting model contained a hydrogen bond acceptor site on the carbonyl oxygen atom and a ring aromatic feature on one of the two aromatic rings in these compounds. This model was used as a template to search an in-house database for novel analogs. We found compounds equal in inhibition potency to atropine and, importantly, six not reported before as antimuscarinics. These results demonstrate that this QSAR approach not only provides a basis for understanding the molecular mechanism of action but a pharmacophore to aid in the discovery and design of novel potent muscarinic antagonists.     

Electrostatic potential surfaces of 5-HT(3)R agonists suggest accessory cation-pi site adjacent to agonist binding domain

Electrostatic potential surface mapping of various aromatic ring systems contained in 5-HT(3)R agonists indicate that some agonists contain an aromatic moiety capable of a favorable cation-pi interaction next to the e-face of pyridine (or its bioisostere). A pharmacophore model has been proposed based on superimposition of two distinct 'aryl' interactions.     

Acidic elements in histamine H3 receptor antagonists

Antagonists of the human histamine H₃ receptor (hH₃R) often contain a second basic moiety, which is well known to boost affinity on this histamine receptor subtype. Here, we prepared compounds with acidic moieties of different pK_a values to figure out that the hH₃R tolerates these functionalities when added to a common pharmacophore blueprint. Depending on the acidic, electronic and steric features the designed ligands showed hH₃R affinities in the nanomolar concentration range. Additionally, selected ligands were tested but failed as dual acting hH₃R/hPPAR (human peroxisome proliferator-activated receptor) ligands.     

Polycyclic dinitriles: a novel class of potent GABAergic insecticides provides a new radioligand, [3H]BIDN

The polycyclic dinitriles are a potent class of insecticides which are non-competitive GABA (γ-aminobutyric a acid) antagonists acting at the convulsant site. Comparison with other classes of GABA convulsant site ligands using molecular modelling has shown significant structural similarities. We have developed a pharmacophore model which unifies this class and some previous classes of GABA convulsants. Key pharmacophore elements are a polarizable functionality separated by a fixed distance from two H-bond accepting elements. This model is based on information from X-ray crystal structures and Sybyl using the Tripos force field. Using this pharmacophore model, numerous structural modifications were explored to enhance understanding of structure-activity relationships at the GABA receptor convulsant site of insects and mammals. A radiolabelled bicyclic dinitrile, [³H]BIDN ([³H]3,3-bis-trifluoromethyl-bicyclo[2,2,1] heptane-2,2-dicarbonitrile), was prepared from this area of chemistry and was used as a probe for the interaction of polycyclic dinitriles at the target site.     

Hierarchical Bayesian Modeling of Pharmacophores in Bioinformatics

Summary One of the key ingredients in drug discovery is the derivation of conceptual templates called pharmacophores. A pharmacophore model characterizes the physicochemical properties common to all active molecules, called ligands, bound to a particular protein receptor, together with their relative spatial arrangement. Motivated by this important application, we develop a Bayesian hierarchical model for the derivation of pharmacophore templates from multiple configurations of point sets, partially labeled by the atom type of each point. The model is implemented through a multistage template hunting algorithm that produces a series of templates that capture the geometrical relationship of atoms matched across multiple configurations. Chemical information is incorporated by distinguishing between atoms of different elements, whereby different elements are less likely to be matched than atoms of the same element. We illustrate our method through examples of deriving templates from sets of ligands that all bind structurally related protein active sites and show that the model is able to retrieve the key pharmacophore features in two test cases.     

In silico modelling and molecular dynamics simulation studies of thiazolidine based PTP1B inhibitors

Protein tyrosine phosphatase 1B (PTP1B) has been identified as a negative regulator of insulin and leptin signalling pathway; hence, it can be considered as a new therapeutic target of intervention for the treatment of type 2 diabetes. Inhibition of this molecular target takes care of both diabetes and obesity, i.e. diabestiy. In order to get more information on identification and optimization of lead, pharmacophore modelling, atom-based 3D QSAR, docking and molecular dynamics studies were carried out on a set of ligands containing thiazolidine scaffold. A six-point pharmacophore model consisting of three hydrogen bond acceptor (A), one negative ionic (N) and two aromatic rings (R) with discrete geometries as pharmacophoric features were developed for a predictive 3D QSAR model. The probable binding conformation of the ligands within the active site was studied through molecular docking. The molecular interactions and the structural features responsible for PTP1B inhibition and selectivity were further supplemented by molecular dynamics simulation study for a time scale of 30 ns. The present investigation has identified some of the indispensible structural features of thiazolidine analogues which can further be explored to optimize PTP1B inhibitors.     

Exploration of new scaffolds as potential MAO-A inhibitors using pharmacophore and 3D-QSAR based in silico screening

Monoamine oxidase-A (MAO-A) inhibitors are of particular importance in the treatment of depressive disorders. Herein described is pharmacophore generation and atom-based 3D-QSAR analysis of previously reported pyrrole based MAO-A inhibitors in order to get insight into their structural requirements responsible for high affinity. The best pharmacophore model generated consisted of four features DHHR: a hydrogen bond donor (D), two hydrophobic groups (H) and an aromatic ring (R). Based on model generated, a statistically valid 3D-QSAR with good predictability was developed. Derived pharmacophore was used as a query to search Zinc ‘clean drug-like’ database. Hits retrieved were passed progressively through filters like fitness score, predicted activity and docking scores. The survived hits present new scaffolds with a potential for MAO-A inhibition.     

Naphtho[2,1-b][1,5] and [1,2-f][1,4]oxazocines as selective NK1 antagonists

Previously we reported on the synthesis and properties of a series of highly potent piperidinyl 2-subsituted-3-cyano-1-naphthamide NK1 antagonists that includes 3 and 4. Here we report our efforts to alleviate a troublesome atropisomeric property of those derivatives by introduction of a tethering bridge that, in addition, could be used to lock the resulting cyclic derivatives in a purported NK1 pharmacophore conformation. Using 3 as a starting point, the naphtho[2,1-b][1,5]oxazocine, 17, was found to contain the optimal ring tether size (8) for retaining NK1 activity, was more NK1 versus NK2 selective, and reduced the number of atropisomers from four to two. Cyclic derivatives 29 and 32, which exist as essentially single atropisomers in the purported pharmacophore conformation, were prepared in the closely related naphtho[1,2-f][1,4]oxazocine series as part of an effort to use mono methyl substitution of the tethering bridge as a conformation stabilizing factor. Both 29 and 32 were found to be less active as NK1 antagonists than the non-methylated parent 28 possibly due to methyl group destabilization of receptor interaction. We discuss the above findings in the context of a previously proposed NK1 pharmacophore model and present a further refinement of that model.     

Interaction of Met297 in the seventh transmembrane segment of the tachykinin NK2 receptor with neurokinin A

We report the use of thiol chemistry to define specific and reversible disulfide interactions of Cys-substituted NK2 receptor mutants with analogues of neurokinin A (NKA) containing single cysteine substitutions. The NKA analogues were N-biotinylated to facilitate the rapid detection of covalent analogue-receptor interactions utilizing streptavidin reactivity. N-biotinyl-[Tyr1,Cys9]NKA, N-biotinyl-[Tyr1,Cys10]NKA were both found to reversibly disulfide bond to the NK2 receptor mutant Met297 --> Cys. This is consistent with the improved affinities of these particular analogues for the Met297 --> Cys receptor as compared with those for the wild-type and Met297 --> Leu receptors. In our three-dimensional model, Met297 occupies the equivalent position in helix 7 to the retinal binding Lys296 in rhodopsin. Binding of the NK2 receptor antagonist [3H]SR 48968 and of 125I-NKA was used to characterize additional receptor mutants. It seems that the aromatic residues Trp99 (helix 3), His198 (helix 5), Tyr266, His267, and Phe270 play an important role in NKA binding as structural determinants. The existence of overlapping SR 48968 and NKA binding sites is also evident. These data suggest that the peptide binding site of the NK2R is at least in part formed by residues buried deep within the transmembrane bundle and that this intramembranous binding domain may correspond to the binding sites for substantially smaller endogenous GPCR ligands.     

SR 142801, a tachykinin NK(3) receptor antagonist,prevents beta(2)-adrenoceptor agonist-induced hyperresponsiveness to neurokinin A in guinea-pig isolated trachea

We investigated whether fenoterol was able to enhance contractile responsiveness to neurokinin A (NKA) on the guinea-pig isolated trachea. We then studied the effects of two inhibitors of nuclear factor kappa B (NFkappaB), gliotoxin and pyrrolidine dithiocarbamate, and of the tachykinin NK(1), NK(2) and NK(3) receptor antagonists, SR 140333, SR 48968 and SR 142801 and determined whether tachykinin receptor gene expression was up-regulated in the trachea after exposure to fenoterol. Fenoterol (0.1 microM, 15 h, 21 degrees C) induced an increased contractile response to NKA (mean of difference in maximal tension between control and fenoterol +/- S.E.M; +0.47 +/- 0.14 g, n = 26, P < 0.01). This hyperresponsiveness was strongly reduced by co-incubation with gliotoxin (0.1 microg/ml) or pyrrolidine dithiocarbamate (0.1 mM) and abolished by SR 140333 (0.1 microM) and SR 142801 (0.1 microM). SR 48968 (0.1 microM) diminished the tracheal contractility to NKA but failed to reduce the hyperreactivity induced by fenoterol. Tachykinin NK(1) receptor (NK(1)R), NK(2) receptor (NK(2)R) and NK(3) receptor (NK(3)R) gene expression was analyzed by semiquantitative RT-PCR. Compared to control tissues, NK(1)R and NK(2)R mRNA expression was increased by about 1.6-fold and 1.4-fold, respectively, in tissues treated with fenoterol. We were unable to detect the presence of NK(3)R mRNA in the guinea-pig trachea. In conclusion, fenoterol induces tracheal hyperresponsiveness to NKA and an up-regulation of NK(1)R and NK(2)R gene expression. The hyperresponsiveness implicates the NFkappaB pathway and is abolished by tachykinin NK(1) (SR 140333) and NK(3) (SR 142801) receptor antagonists.