Pharmacophore-based discovery of 3,4-disubstituted pyrrolidines as a novel class of monoamine transporter inhibitors

3,4-Disubstituted pyrrolidines were discovered as a novel class of monoamine transporter inhibitors through 3-D database pharmacophore searching using a new pharmacophore model. The most potent analogue 12 has Ki values of 0.084 microM in [3H]mazindol binding, 0.20, 0.23, and 0.031 microM in inhibition of dopamine (DA), serotonin (SER), and norepinephrine (NE) reuptake, respectively. Functional antagonism testing in vitro showed that 11 and 12 are weak cocaine antagonists.     

Pharmacophore modeling of human adenosine receptor A2A antagonists

Three-dimensional pharmacophore models of human adenosine receptor A_2A antagonists were developed based on 23 diverse compounds selected from a large number of A_2A antagonists. The best pharmacophore model, Hypo1, contained five features: one hydrogen bond donor , three hydrophobic points and one ring aromatic. Its correlation coefficient, root mean square deviation, and cost difference values were 0.955, 0.921 and 84.4, respectively, suggested that the Hypo1 model was reasonable and reliable. This model was validated by three methods: a test set of 106 diverse compounds, a simulated virtual screening, and superimposition with the crystal structure of A_2A receptor. The results showed that Hypo1 was not only in agreement with the A_2A crystal structure and literature reports, but also well identified active A_2A antagonists from the virtual database. This methodology provides helpful information and a robust tool for the discovery of potent A_2A antagonists.     

Molecular basis for lysophosphatidic acid receptor antagonist selectivity

Recent characterization of lysophosphatidic acid (LPA) receptors has made possible studies elucidating the structure-activity relationships (SAR) for agonist activity at individual receptors. Additionally, the availability of these receptors has allowed the identification of antagonists of LPA-induced effects. Two receptor-subtype selective LPA receptor antagonists, one selective for the LPA1/EDG2 receptor (a benzyl-4-oxybenzyl N-acyl ethanolamide phosphate, NAEPA, derivative) and the other selective for the LPA3/EDG7 receptor (diacylglycerol pyrophosphate, DGPP, 8:0), have recently been reported. The receptor SAR for both agonists and antagonists are reviewed, and the molecular basis for the difference between agonism and antagonism as well as for receptor-subtype antagonist selectivity identified by molecular modeling is described. The implications of the newly available receptor-subtype selective antagonists are also discussed.     

A novel potential therapeutic avenue for autism: design, synthesis and pharmacophore generation of SSRIs with dual action

Autism symptoms are currently modulated by Selective Serotonin Reuptake Inhibitors (SSRIs). SSRIs slow onset of action limits their efficiency. The established synergistic activity of SSRIs and 5HT(1B/1D) autoreceptors antagonists motivated us to incorporate SSRIs and 5HT(1B/1D) antagonists in one 'hybrid' molecule. A library of virtual 'hybrid' molecules was designed using the tethering technique. A pharmacophore model was generated derived from 16 structurally diverse SSRIs (K(i)=0.013-5000 nM) and used as 3D query. Compounds with fit values (≥2) were chosen for synthesis and subsequent in vitro biological evaluation. Our pharmacophore model is a promising milestone to a class of SSRIs with dual action.     

Pharmacophore-based design, synthesis, biological evaluation, and 3D-QSAR studies of aryl-piperazines as alpha(1)-adrenoceptor antagonists

Phenyl-piperazines were designed and synthesized based on pharmacophore for uro-selective alpha(1)-adrenoceptor antagonists and 3D chemical database searching. Within this series, three compounds, 2, 3, and 13, showed similar or better alpha(1)-AR antagonistic activity compared with prazosin. The 3D-QSAR study of these compounds may provide useful information for the development of novel aryl-piperazines as uro-selective alpha(1)-adrenoceptor antagonists, which can be used for the treatment of BPH with fewer side effects.     

3D-Quantitative structure-activity relationship and docking studies of the tachykinin NK3 receptor

The tachykinin NK(3) receptor (NK(3)R) is a novel drug target for schizophrenia and drug abuse. Since few non-peptide antagonists of this G protein-coupled receptor are available, we have initiated this study to gain a better understanding of the structure-activity relationships of NK(3) antagonist compounds. We developed a 3D comparative molecular similarity index analysis (CoMSIA) model that gave cross-validated PLS values with q(2) >0.5 which were validated using a test set. We also describe the development of a homology model of the NK(3)R. The model was then used to develop a pharmacophore for docked ligands. This pharmacophore showed two aromatic, two hydrogen donor and one acceptor/aromatic points. These data will be useful for future structure-based drug discovery of ligands for the NK(3)R.     

Discovery of potent inhibitors of pseudomonal quorum sensing via pharmacophore modeling and in silico screening

HipHop-Refine was employed to derive a binding hypothesis for pseudomonal quorum sensing (QS) antagonists. The model was employed as 3D search query to screen the National Cancer Institute (NCI) database. One of the hits illustrated nanomolar QS inhibitory activity. The fact that this compound contained tetravalent lead (Pb) prompted us to evaluate the activities of phenyl mercuric nitrate and thimerosal, both fit the binding pharmacophore. The two mercurials illustrated nanomolar to low micromolar IC50 inhibitory values against pseudomonal QS. The three compounds represent a new class of QS inhibitors.     

Pharmacophore modeling and virtual screening to identify potential RET kinase inhibitors

Chemical features based 3D pharmacophore model for REarranged during Transfection (RET) tyrosine kinase were developed by using a training set of 26 structurally diverse known RET inhibitors. The best pharmacophore hypothesis, which identified inhibitors with an associated correlation coefficient of 0.90 between their experimental and estimated anti-RET values, contained one hydrogen-bond acceptor, one hydrogen-bond donor, one hydrophobic, and one ring aromatic features. The model was further validated by a testing set, Fischer’s randomization test, and goodness of hit (GH) test. We applied this pharmacophore model to screen NCI database for potential RET inhibitors. The hits were docked to RET with GOLD and CDOCKER after filtering by Lipinski’s rules. Ultimately, 24 molecules were selected as potential RET inhibitors for further investigation.     

Pharmacophore identification of alpha(1A)-adrenoceptor antagonists

A chemical feature based pharmacophore model was developed for alpha(1A)-adrenoceptor antagonists by HypoGen module implemented in catalyst software package. The best scoring pharmacophore hypothesis, Hypo1, consisted of four important chemical features (one positive ion, one hydrogen-bond donor, one aromatic ring, and one hydrophobic group). The results of our study provide a valuable tool in designing new leads with desired biological activity by virtual screening.     

Identification of a phosphothionate analogue of lysophosphatidic acid (LPA) as a selective agonist of the LPA3 receptor

Lysophosphatidic acid (LPA) is a bioactive lysophospholipid mediator that acts through G protein-coupled receptors. Most cell lines in culture express one or more LPA receptors, making it difficult to assign a response to specific LPA receptors. Dissection of the signaling properties of LPA has been hampered by lack of LPA receptor subtype-specific agonists and antagonists. The present study characterizes an ester-linked thiophosphate derivative (1-oleoyl-2-O-methyl-rac-glycerophosphothionate, OMPT) of LPA. OMPT is a functional LPA analogue with potent mitogenic activity in fibroblasts. In contrast to LPA, OMPT does not couple to the pheromone response through the LPA(1) receptor in yeast cells. OMPT induces intracellular calcium increases efficiently in LPA(3) receptor-expressing Sf9 cells but poorly in LPA(2) receptor-expressing cells. Guanosine 5'-O-(3-[(35)S]thio)triphosphate binding assays in mammalian cells showed that LPA exhibits agonistic activity on all three LPA receptor subtypes, whereas OMPT has a potent agonistic effect only on the LPA(3) receptor. In transiently transfected HEK293 cells, OMPT stimulates mitogen-activated protein kinases through the LPA(3) but not the LPA(1) or LPA(2) receptors. Furthermore, OMPT-induced intracellular calcium mobilization in mammalian cells is efficiently inhibited by the LPA(1)/LPA(3) receptor-selective antagonist VPC12249. These results establish that OMPT is an LPA(3)-selective agonist. OMPT binding to the LPA(3) receptor in mammalian cells is sufficient to elicit multiple responses, including activation of G proteins, calcium mobilization, and activation of mitogen-activated protein kinases. Thus OMPT offers a powerful probe for the dissection of LPA signaling events in complex mammalian systems.     

Structural Determinants of the Transient Receptor Potential 1 (TRPV1) Channel Activation by Phospholipid Analogs

The transient receptor potential vanilloid 1 (TRPV1) ion channel is a polymodal protein that responds to various stimuli, including capsaicin (the pungent compound found in chili peppers), extracellular acid, and basic intracellular pH, temperatures close to 42 °C, and several lipids. Lysophosphatidic acid (LPA), an endogenous lipid widely associated with neuropathic pain, is an agonist of the TRPV1 channel found in primary afferent nociceptors and is activated by other noxious stimuli. Agonists or antagonists of lipid and other chemical natures are known to possess specific structural requirements for producing functional effects on their targets. To better understand how LPA and other lipid analogs might interact and affect the function of TRPV1, we set out to determine the structural features of these lipids that result in the activation of TRPV1. By changing the acyl chain length, saturation, and headgroup of these LPA analogs, we established strict requirements for activation of TRPV1. Among the natural LPA analogs, we found that only LPA 18:1, alkylglycerophosphate 18:1, and cyclic phosphatidic acid 18:1, all with a monounsaturated C18 hydrocarbon chain activate TRPV1, whereas polyunsaturated and saturated analogs do not. Thus, TRPV1 shows a more restricted ligand specificity compared with LPA G-protein-coupled receptors. We synthesized fatty alcohol phosphates and thiophosphates and found that many of them with a single double bond in position Δ9, 10, or 11 and Δ9 cyclopropyl group can activate TRPV1 with efficacy similar to capsaicin. Finally, we developed a pharmacophore and proposed a mechanistic model for how these lipids could induce a conformational change that activates TRPV1.     

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.     

Pharmacophore models generation by catalyst and phase consensus-based virtual screening protocol against PI3Kα inhibitors

Phosphoinositide 3-kinases (PI3Ks) family has emerged as promising targets for novel therapeutic agents against neoplastic diseases. Pharmacophore and 3D-quantitative structure–activity relationship modelling were applied to study the structure–activity relationship of PI3K inhibitors. The best HypoGen pharmacophore hypothesis Hypo1 with a correlation coefficient of 0.961 consists of one hydrogen-bond acceptor, one hydrogen-bond donor and two hydrophobic features, whereas the best phase hypothesis AADRRR.378 with favourable statistics (q² = 0.7368, r² = 0.9863) has two hydrogen-bond acceptors, one hydrogen-bond donor and three ring aromatic features. Multiple methods, such as Fischer validation, molecular docking and mapping of test set molecules, were carried out to validate these pharmacophore models. Furthermore, a comparative molecular similarity indices analysis candidate hypothesis model was generated as a supplement of pharmacophore hypothesis. Detailed protein–ligand binding information obtained by Glide was utilised in compound optimisation and virtual screening. A molecular database of 133 known inhibitors and 6179 decoys was built for a screening test to quantitatively analyse various hypotheses and scoring parameters. Finally, we designed a workflow integrating HypoGen pharmacophore searching, phase pharmacophore searching and molecular docking for screening the database. With an improved criterion of enrichment factor (EF = 17.43) and ROC curve (AUC = 0.946), this workflow would provide us an original method for novel PI3K inhibitors.     

Synthesis and pharmacological evaluation of second-generation phosphatidic acid derivatives as lysophosphatidic acid receptor ligands

Short-chain phosphatidic acid derivatives, dioctanoyl glycerol pyrophosphate (DGPP 8:0, 1) and phosphatidic acid 8:0 (PA 8:0, 2), were previously identified as subtype-selective LPA(1) and LPA(3) receptor antagonists. Recently, we reported that the replacement of the phosphate headgroup by thiophosphate in a series of fatty alcohol phosphates (FAP) improves agonist as well as antagonist activities at LPA GPCR. Here, we report the synthesis of stereoisomers of PA 8:0 analogs and their biological evaluation at LPA GPCR, PPARgamma, and ATX. The results indicate that LPA receptors stereoselectively interact with glycerol backbone modified ligands. We observed entirely stereospecific responses by dioctyl PA 8:0 compounds, in which (R)-isomers were found to be agonists and (S)-isomers were antagonists of LPA GPCR. From this series, we identified compound 13b as the most potent LPA(3) receptor subtype-selective agonist (EC(50)=3 nM), and 8b as a potent and selective LPA(3) receptor antagonist (K(i)=5 nM) and inhibitor of ATX (IC(50)=600 nM). Serinediamide phosphate 19b was identified as an LPA(3) receptor specific antagonist with no effect on LPA(1), LPA(2), and PPARgamma.     

In silico investigation of interactions between human cannabinoid receptor-1 and its antagonists

Cannabinoid receptor-1 (CB(1)) is widely expressed in the central nervous system and plays a vital role in regulating food intake and energy expenditure. CB(1) antagonists such as Rimonabant have been used in clinic to inhibit food intake, and therefore reduce body weight in obese animals and humans. To investigate the binding modes of CB(1) antagonists to the receptor, both receptor- and ligand-based methods were implemented in this study. At first, a pharmacophore model was generated based on 31 diverse CB(1) antagonists collected from literature. A test set validation and a simulated virtual screening evaluation were then performed to verify the reliability and discriminating ability of the pharmacophore. Meanwhile, the homology model of CB(1) receptor was constructed based on the crystal structure of human β (2) adrenergic receptor (β (2)-AR). Several classical antagonists were then docked into the optimized homology model with induced fit docking method. A hydrogen bond between the antagonists and Lys192 on the third transmembrane helix of the receptor was formed in the docking study, which has proven to be critical for receptor-ligand interaction by biological experiments. The structure obtained from induced fit docking was then confirmed to be a reliable model for molecular docking from the result of the simulated virtual screening. The consistency between the pharmacophore and the homology structure further proved the previous observation. The built receptor structure and antagonists' pharmacophore should be useful for the understanding of inhibitory mechanism and development of novel CB(1) antagonists.     

Stereochemical properties of lysophosphatidic acid receptor activation and metabolism

Ligand recognition by G protein-coupled receptors (GPCR), as well as substrate recognition by enzymes, almost always shows a preference for a naturally occurring enantiomer over the unnatural one. Recognition of lysophosphatidic acid (LPA) by its receptors is an exception, as both the natural L (R) and unnatural D (S) stereoisomers of LPA are equally active in bioassays. In contrast to the enantiomers of LPA, analogs of N-acyl-serine phosphoric acid (NASPA) and N-acyl-ethanolamine phosphoric acid (NAEPA), which contain a serine and an ethanolamine backbone, respectively, in place of glycerol, are recognized in a stereoselective manner. This stereoselective interaction may lead to the development of receptor subtype-selective antagonists. In the present study, we review the stereochemical aspects of LPA pharmacology and describe the chemical synthesis of pure LPA enantiomers together with their ligand-binding properties toward the LPA1, LPA2, and LPA3 receptors and their metabolism by lipid phosphate phosphatase 1 (LPP1). Finally, we evaluate the concept of stereopharmacology in developing novel ligands for LPA receptors.     

Initial structure-activity relationships of lysophosphatidic acid receptor antagonists: discovery of a high-affinity LPA1/LPA3 receptor antagonist

A recently reported dual LPA1/LPA3 receptor antagonist (VPC12249, 1) has been modified herein so as to optimize potency and selectivity at LPA receptors. Compounds containing variation in the acyl lipid chain and linker region have been synthesized and screened for activity at individual LPA receptors. LPA1-selective (14b) and LPA3-selective (10g,m) compounds of modest potency have been discovered. Additionally, 2-pyridyl derivative 10t exhibits a Ki value of 18 nM at the LPA1 receptor and is significantly more potent than 1 at the LPA3 receptor. This paper describes the synthetic methods, biological evaluation, and structure-activity relationships (SARs) of LPA receptor antagonists.     

The utilization of a unified pharmacophore query in the discovery of new antagonists of the adenosine receptor family

Pharmacophore queries from previously known potent selective A3 antagonists were generated by Chem-X. These queries were used to search a pharmacophore database of diverse compounds (CNS-Set). In vitro assays of 186 'hits' yielded over 30 active compounds, for four adenosine receptor subtypes. This search strategy may also be applicable to the discovery of new ligands via receptor homology data.     

Structure-based drug design identifies novel LPA3 antagonists

Compound 5 ([5-(3-nitrophenoxy)-1,3-dioxo-1,3-dihydro-2-isoindol-2-yl]acetic acid) was identified as a weak selective LPA₃ antagonist (IC₅₀ = 4504 nM) in a virtual screening effort to optimize a dual LPA_{2 and 3} antagonist. Structure-based drug design techniques were used to prioritize similarity search matches of compound 5. This strategy rapidly identified 10 novel antagonists. The two most efficacious compounds identified inhibit activation of the LPA₃ receptor by 200 nM LPA with IC₅₀ values of 752 nM and 2992 nM. These compounds additionally define changes to our previously reported pharmacophore that will improve its ability to identify more potent and selective LPA₃ receptor antagonists. The results of the combined computational and experimental screening are reported.