G-Protein coupled receptors: models, mutagenesis, and drug design

Probing the orientation of inhibitor and epoxy-eicosatrienoic acid binding in the active site of soluble epoxide hydrolase

Soluble epoxide hydrolase (sEH) is an important therapeutic target of many diseases, such as chronic obstructive pulmonary disease (COPD) and diabetic neuropathic pain. It acts by hydrolyzing and thus regulating specific bioactive long chain polyunsaturated fatty acid epoxides (lcPUFA), like epoxyeicosatrienoic acids (EETs). To better predict which epoxides could be hydrolyzed by sEH, one needs to dissect the important factors and structural requirements that govern the binding of the substrates to sEH. This knowledge allows further exploration of the physiological role played by sEH. Unfortunately, a crystal structure of sEH with a substrate bound has not yet been reported. In this report, new photoaffinity mimics of a sEH inhibitor and EET regioisomers were prepared and used in combination with peptide sequencing and computational modeling, to identify the binding orientation of different regioisomers and enantiomers of EETs into the catalytic cavity of sEH. Results indicate that the stereochemistry of the epoxide plays a crucial role in dictating the binding orientation of the substrate.

Design, synthesis and pharmacological evaluation of indolylsulfonamide amines as potent and selective 5-HT6 receptor antagonists

A series of N'-[3-(indole-1-sulfonyl) aryl]-N,N-dimethyl ethane-1,2-diamines and N'-[3-(indole-1-sulfonyl) aryl]-N,N-dimethyl propane-1,3-diamines was designed and synthesized as 5-HT6 receptor ligands. These compounds, when screened in a functional reporter gene-based assay, displayed potent antagonistic activity with Kb values in the range of 1.8-60 nM. The lead compound 9y has shown good ADME surrogate properties, acceptable pharmacokinetic profile and is active in animal models of cognition like novel object recognition test and Morris water maze. It was selected for detailed profiling.

Genomics, Other “Omic” Technologies, Personalized Medicine, and Additional Biotechnology-Related Techniques

The products resulting for biotechnologies continue to grow at an exponential rate, and the expectations are that an even greater percentage of drug development will be in the area of the biologics. In 2011, worldwide there were over 800 new biotech drugs and treatments in development including 23 antisense, 64 cell therapy, 50 gene therapy, 300 monoclonal antibodies, 78 recombinant proteins, and 298 vaccines (PhRMA 2012). Pharmaceutical biotechnology techniques are at the core of most methodologies used today for drug discovery and development of both biologics and small molecules. While recombinant DNA technology and hybridoma techniques were the major methods utilized in pharmaceutical biotechnology through most of its historical timeline, our ever-widening understanding of human cellular function and disease processes and a wealth of additional and innovative biotechnologies have been, and will continue to be, developed in order to harvest the information found in the human genome. These technological advances will provide a better understanding of the relationship between genetics and biological function, unravel the underlying causes of disease, explore the association of genomic variation and drug response, enhance pharmaceutical research, and fuel the discovery and development of new and novel biopharmaceuticals. These revolutionary technologies and additional biotechnology-related techniques are improving the very competitive and costly process of drug development of new medicinal agents, diagnostics, and medical devices. Some of the technologies and techniques described in this chapter are both well established and commonly used applications of biotechnology producing potential therapeutic products now in development including clinical trials. New techniques are emerging at a rapid and unprecedented pace and their full impact on the future of molecular medicine has yet to be imagined.

Classification of 5-HT(1A) receptor agonists and antagonists using GA-SVM method

AIM

To construct a reliable computational model for the classification of agonists and antagonists of 5-HT(1A) receptor.

METHODS

Support vector machine (SVM), a well-known machine learning method, was employed to build a prediction model, and genetic algorithm (GA) was used to select the most relevant descriptors and to optimize two important parameters, C and r of the SVM model. The overall dataset used in this study comprised 284 ligands of the 5-HT(1A) receptor with diverse structures reported in the literatures.

RESULTS

A SVM model was successfully developed that could be used to predict the probability of a ligand being an agonist or antagonist of the 5-HT(1A) receptor. The predictive accuracy for training and test sets was 0.942 and 0.865, respectively. For compounds with probability estimate higher than 0.7, the predictive accuracy of the model for training and test sets was 0.954 and 0.927, respectively. To further validate our model, the receiver operating characteristic (ROC) curve was plotted, and the Area-Under-the-ROC- Curve (AUC) value was calculated to be 0.883 for training set and 0.906 for test set.

CONCLUSION

A reliable SVM model was successfully developed that could effectively distinguish agonists and antagonists among the ligands of the 5-HT(1A) receptor. To our knowledge, this is the first effort for the classification of 5-HT(1A) receptor agonists and antagonists based on a diverse dataset. This method may be used to classify the ligands of other members of the GPCR family.

Comparison analysis of primary ligand-binding sites in seven-helix membrane proteins

Seven-helix transmembrane proteins, including the G-protein-coupled receptors (GPCRs), mediate a broad range of fundamental cellular activities through binding to a wide range of ligands. Understanding the structural basis for the ligand-binding selectivity of these proteins is of significance to their structure-based drug design. Comparison analysis of proteins' ligand-binding sites provides a useful way to study their structure-activity relationships. Various computational methods have been developed for the binding-site comparison of soluble proteins. In this work, we applied this approach to the analysis of the primary ligand-binding sites of 92 seven-helix transmembrane proteins. Results of the studies confirmed that the binding site of bacterial rhodopsins is indeed different from all GPCRs. In the latter group, further comparison of the binding sites indicated a group of residues that could be responsible for ligand-binding selectivity and important for structure-based drug design. Furthermore, unexpected binding-site dissimilarities were observed among adrenergic and adenosine receptors, suggesting that the percentage of the overall sequence identity between a target protein and a template protein alone is not sufficient for selecting the best template for homology modeling of seven-helix membrane proteins. These results provided novel insight into the structural basis of ligand-binding selectivity of seven-helix membrane proteins and are of practical use to the computational modeling of these proteins.

Comparative 3D QSAR study on β(1)-, β(2)-, and β(3)-adrenoceptor agonists

A quantitative structure–activity relationship study of tryptamine-based derivatives of β₁-, β₂-, and β₃-adrenoceptor agonists was conducted using comparative molecular field analysis (CoMFA). Correlation coefficients (cross-validated r²) of 0.578, 0.595, and 0.558 were obtained for the three subtypes, respectively, in three different CoMFA models. All three CoMFA models have different steric and electrostatic contributions, implying different requirements inside the binding cavity. The CoMFA coefficient contour plots of the three models and comparisons among these plots provide clues regarding the main chemical features responsible for the biological activity variations and also result in predictions which correlate very well with the observed biological activity. Based on the analysis, a summary regeospecific description of the requirements for improving β-adrenoceptor subtype selectivity is given.

The beta-D-glucose scaffold as a beta-turn mimetic

Activity and selectivity are typically the first considerations when designing a drug. However, absorption, distribution, metabolism, excretion, and toxicity (ADMET) are equally important considerations. Peptides can provide a combination of potent binding and exquisite selectivity, as evidenced by their pervasive use as enzymes, hormones, and signaling agents within living systems. In particular, peptidic turn motifs are key elements of molecular recognition. They may be found at the exposed surfaces of globular proteins, where they are available for binding interactions with other peptides and small molecules. However, despite these advantages, peptides often make poor drugs. The amide backbone is subject to rapid enzymatic proteolysis, resulting in short half-lives. Furthermore, the ability of the amide backbone to hydrogen bond with water restricts its ability to cross membranes and, consequentially, results in poor oral bioavailability. Accordingly, the development of nonpeptidic scaffolds that mimic peptidic turn motifs represents a promising means of converting peptidic agents into more drugable molecules. In this Account, we describe the design and synthesis of beta-turn mimetics that use a beta-D-glucose scaffold, the first use of a sugar scaffold for this purpose. Somatostatin (SRIF) is a small protein (14 amino acid residues) human hormone; a shorter (6 amino acid residues) synthetic peptide, L-363,301, is a fully peptidal agonist. These two cyclic peptides share the beta-turn motif comprising Phe(7)-Trp(8)-Lys(9)-Thr(10) (d-Trp(8) in the case of L-363,301), of which the tryptophan and lysine residues in the i + 1 and i + 2 positions, respectively, are critical for binding. In 1988, we initiated a program that tested and validated the then-novel proposition that the beta-D-glucose scaffold can mimic the beta-turn in L-363,301. The beta-D-glucose scaffold proved to be an attractive mimic of a beta-turn in part because it permits the convenient attachment of amino acid side chains via facile etherification reactions, rather than carbon-carbon bond formations; it is also an inexpensive starting material with well-defined stereochemistry. From the beginning, biological assays were used alongside physical measurements to assess the relevance of the design. Our first two synthetic targets, compounds 6 and 7, bound the SRIF receptors on benchmark (AtT-20) cells, albeit weakly, consistent with the objective of the design. Subsequently, a better ligand (8) and two congeners were found to be agonists at the SRIF receptors, providing convincing evidence that the peptide backbone is not required for receptor binding or signal transduction. The unexpectedly high level of receptor affinity of selected analogs, as well as the fortuitous discovery that our peptidomimetics were active against several chemically distinct receptors, led us to hypothesize that these monosaccharides could access multiple potential binding modes. Our later studies of this sugar scaffold confirmed this property, which we termed pseudosymmetry, whereby multiple similar but nonidentical motifs are displayed within a single analog. We propose the presence of pseudosymmetry to be an element of privilege and an advantage for lead discovery.

Molecular modeling of the peptide agonist-binding site in a neurokinin-2 receptor

The neurokinin-2 receptor is a member of the rhodopsin family of G-protein coupled receptors, which represents one of the most relevant target families in small-molecule drug design. NK-2 receptors have been implicated in playing a pathophysiological role in asthma. Activation of the NK-2 receptor by its endogenous peptide agonist, tachykinins, is associated with diverse biological responses like bronchoconstriction, vasodepression, and regulation of endocrine functions. Agonist binding to the receptor is a crucial event in initiating signaling, and therefore characterization of the structural features of the agonists can reveal the molecular basis of receptor activation and help in rational design of novel therapeutics. In this study a molecular model for the interaction of the primary ligand NKA with its G-protein coupled receptor neurokinin-2 receptor has been developed. A three-dimensional model for the NK-2 receptor has been generated by homology modeling using rhodopsin as a template. A knowledge based docking of the NMR derived bioactive conformation of NKA to the receptor has been performed utilizing the ligand binding data obtained from the photoaffinity labeling and site-directed mutagenesis studies. The molecular model for the NKA/NK-2 receptor complex thus obtained sheds light on the topographical features of the binding pocket of the receptor and provides atomic insight into the biochemical data currently available for the receptor. The results of the receptor modeling studies have been used to discuss the molecular determinants for NK-2 receptor selectivity.

The cholesterol-complexing agent digitonin modulates ligand binding of the bovine hippocampal serotonin1Areceptor

The serotonin(1A) (5-HT(1A)) receptor is an important member of the superfamily of seven transmembrane domain G-protein-coupled receptors. We have examined the modulatory role of cholesterol on the ligand binding of the bovine hippocampal 5-HT(1A) receptor by cholesterol complexation in native membranes using digitonin. Complexation of cholesterol from bovine hippocampal membranes using digitonin results in a concentration-dependent reduction in specific binding of the agonist 8-OH-DPAT and antagonist p-MPPF to 5-HT(1A) receptors. The corresponding changes in membrane order were monitored by analysis of fluorescence polarization data of the membrane depth-specific probes, DPH and TMA-DPH. Taken together, our results point out the important role of membrane cholesterol in maintaining the function of the 5-HT(1A) receptor. An important aspect of these results is that non-availability of free cholesterol in the membrane due to complexation with digitonin rather than physical depletion is sufficient to significantly reduce the 5-HT(1A) receptor function. These results provide a comprehensive understanding of the effects of the sterol-complexing agent digitonin in particular, and the role of membrane cholesterol in general, on the 5-HT(1A) receptor function.

Dopamine D3 receptor ligands—Recent advances in the control of subtype selectivity and intrinsic activity

Various pharmacological studies have implicated the dopamine D(3) receptor as an interesting therapeutic target in the treatment of different neurological disorders. Because of these putative therapeutic applications, D(3) receptor ligands with diverse intrinsic activities have been an active field of research in recent years. Separation of purely D(3)-mediated drug effects from effects produced by interactions with similar biogenic amine receptors allows to verify the therapeutic impact of D(3) receptors and to reduce possible side-effects caused by "promiscuous" receptor interactions. The requirement to gain control of receptor selectivity and in particular subtype selectivity has been a challenging task in rational drug discovery for quite a few years. In this review, recently developed structural classes of D(3) ligands are discussed, which cover a broad spectrum of intrinsic activities and show interesting selectivities.

Location of the hydrophobic pocket in the binding site of fentanyl analogs in the µ-opioid receptor

Fentanyl is a highly potent and clinically widely used narcotic analgesic. The synthesis of its analogs remains a challenge in an attempt to develop highly selective ?-opioid receptor agonists with specific pharmacological properties. In this paper, the use of flexible molecular docking of several specific fentanyl analogs to the ?-opioid receptor model, in order to test the hypothesis that the hydrophobic pocket accommodates alkyl groups at position 3 of the fentanyl skeleton, is described. The stereoisomers of the following compounds were studied: cis- and trans-3-methylfentanyl, 3,3-dimethylfentanyl, cis- and trans-3-ethylfentanyl, cis- and trans-3-propylfentanyl, cis-3-isopropylfentanyl and cis-3-benzylfentanyl. The optimal position and orientation of these fentanyl analogs in the binding pocket of the ?-receptor, explaining their enantiospecific potency, were determined. It was found that the 3-alkyl group of cis-3R,4S and trans-3S,4S stereoisomers of all the active compounds occupies the hydrophobic pocket between TM5, TM6 and TM7, made up of the amino acids Trp318 (TM7), Ile322 (TM7), Ile301 (TM6) and Phe237 (TM5). However, the fact that this hydrophobic pocket can also accommodate the bulky 3-alkyl substituents of the two inactive compounds: cis-3-isopropylfentanyl, and cis-3-benzylfentanyl, indicates that this hydrophobic pocket in the employed receptor model is probably too large. .

Approaches to the molecular modeling of 7-transmembrane helical receptors

7-Transmembrane helical receptors (7TMs) represent the single most important class of target for drug therapy; therefore, a great deal of effort has gone into computational studies of their structures. Historically, these were based on low resolution electron diffraction data, together with the use of computational methods such as multiple sequence alignments, distance geometry, and molecular mechanics calculations. In the year 2000 the situation changed when the first crystal structure of a 7TM, was published. It was then possible to use the homology modeling techniques to generate more accurate models of these proteins. This unit reviews the modeling of 7TMs and describes in detail how homology modeling can be used to build a structure of the 5-HT2a receptor. Special attention is given to the initial sequence alignment, the most important step in the process. Use of automatic alignment programs often produces incorrect results, and manual intervention is necessary before proceeding further.

Regulation of G proteins by human 5-HT1a receptor TM3/i2 and TM5/i3 loop peptides

A bioactive synthetic 11 amino acid peptide probe (P11) was constructed according to the published sequence of the human 5HT1a receptor. The probe was used to enhance understanding of cytoplasmic loop 2/G protein coupling and activation. Additionally, two peptides (P8, P9) from the cytoplasmic loop 3 region were synthesized and studied. These probes were tested in a model system of human 5HT1a receptor stably expressed in Chinese Hamster Ovary cells. In agonist inhibition studies, P11 was active in all three receptor preparations tested: whole cells, membrane bound, and solubilized. In analyses of the membrane bound receptor system, P11 demonstrated uncompetitive inhibition characteristics. When forskolin-stimulated cAMP levels were measured, P11 was inactive in this negatively coupled system. Utilizing a [35S]gamma-S-GTP incorporation assay, P11 was unable to stimulate G protein incorporation of GTP. While P8 and P9 were also broadly active as non-competitive agonist inhibitors, their characteristics differed in the signal transduction system. P8 and P9 did not significantly change forskolin-stimulated cAMP levels. However, P8 increased [35S]gamma-S-GTP incorporation, while P9 decreased incorporation. Thus, P11, a synthetic peptide from the TM3/i2 region of the receptor, provides suggestive evidence that this receptor region is involved in G protein coupling but not activation. On the other hand, P8 and P9 activities suggest that the TM5/i3 region is involved in both coupling to and regulation of G protein activity. The current evidence from these cytoplasmic loop regions is discussed in the overall context of an emerging model for human 5HT1a receptor-G protein interactions.

The structural evolution of dopamine D3 receptor ligands: structure-activity relationships and selected neuropharmacological aspects

"Evolution consists largely of molecular tinkering."-Following the famous concept of the molecular geneticist and medicine Nobel laureate François Jacob, in this review we describe the structural evolution of dopamine D3 receptor ligands from the natural agonist dopamine (DA) to highly potent and subtype selective new agents by bioisosteric tinkering with well-established and privileged or novel and fancy chemical functionalities and scaffolds. Some of the more than 200 ligands presented herein have already achieved therapeutic or scientific value up to now, some will most likely achieve it in the future. Hence, great importance is not only attached to the relationship between structure and activity of the ligands, but also to their utility as pharmacological tools in animal models or as therapeutics in patients with neurological diseases or other disorders.

A new series of pyridinyl-alkynes as antagonists of the metabotropic glutamate receptor 5 (mGluR5)

Synthesis and some structure-activity relationships for a new series of propargyl ethers as mGluR5 antagonists are reported.

Application of the Goldilocks effect to the design of potent and selective inhibitors of phenylethanolamine N-methyltransferase: balancing pKa and steric effects in the optimization of 3-methyl-1,2,3,4-tetrahydroisoquinoline inhibitors by beta-fluorination

3-Methyl-1,2,3,4-tetrahydroisoquinolines (3-methyl-THIQs) are potent inhibitors of phenylethanolamine N-methyltransferase (PNMT), but are not selective due to significant affinity for the alpha(2)-adrenoceptor. Fluorination of the methyl group lowers the pK(a) of the THIQ amine from 9.53 (CH(3)) to 7.88 (CH(2)F), 6.42 (CHF(2)), and 4.88 (CF(3)). This decrease in pK(a) results in a reduction in affinity for the alpha(2)-adrenoceptor. However, increased fluorination also results in a reduction in PNMT inhibitory potency, apparently due to steric and electrostatic factors. Biochemical evaluation of a series of 3-fluoromethyl-THIQs and 3-trifluoromethyl-THIQs showed that the former were highly potent inhibitors of PNMT, but were often nonselective due to significant affinity for the alpha(2)-adrenoceptor, while the latter were devoid of alpha(2)-adrenoceptor affinity, but also lost potency at PNMT. 3-Difluoromethyl-7-substituted-THIQs have the proper balance of both steric and pK(a) properties and thus have enhanced selectivity versus the corresponding 3-fluoromethyl-7-substituted-THIQs and enhanced PNMT inhibitory potency versus the corresponding 3-trifluoromethyl-7-substituted-THIQs. Using the "Goldilocks Effect" analogy, the 3-fluoromethyl-THIQs are too potent (too hot) at the alpha(2)-adrenoceptor and the 3-trifluoromethyl-THIQs are not potent enough (too cold) at PNMT, but the 3-difluoromethyl-THIQs are just right. They are both potent inhibitors of PNMT and highly selective due to low affinity for the alpha(2)-adrenoceptor. This seems to be the first successful use of the beta-fluorination of aliphatic amines to impart selectivity to a pharmacological agent while maintaining potency at the site of interest.

Steric interactions and the activity of fentanyl analogs at the μ-opioid receptor

Fentanyl is a highly potent and clinically widely used narcotic analgesic. The synthesis of its analogs remains a challenge in the attempt to develop highly selective mu-opioid receptor agonists with specific pharmacological properties. In this paper, the use of flexible molecular docking in a study of the formation of complexes between a series of active fentanyl analogs and the mu-opioid receptor is described. The optimal position and orientation of fourteen fentanyl analogs in the binding pocket of the mu-receptor were determined. The major receptor amino acids and the ligand functional groups participating in the complex formation were identified. Stereochemical effects on the potency and binding are explained. The proposed model of ligand-receptor binding is in agreement with point mutation experiments explaining the role of the amino acids: Asp147, Tyr148, Asn230, His297, Trp318, His319, Cys321, and Tyr326 in the complex formation. In addition, the following amino acids were identified as being important for ligand binding or receptor activation: Ile322, Gly325, Val300, Met203, Leu200, Val143, and Ile144.

Focused library design in GPCR projects on the example of 5-HT(2c) agonists: comparison of structure-based virtual screening with ligand-based search methods

The aim of this study was to investigate the usefulness of structure-based virtual screening (VS) for focused library design in G protein-coupled receptors (GPCR) projects on the example of 5-HT(2c) agonists. We compared the performance of structure-based VS against two different homology models using FRED for docking and ScreenScore, FlexX, and PMF for rescoring with the results of 12 ligand-based similarity searches using four different query compounds and three different similarity metrics (Daylight, FTree, Phacir). The result of the similarity search showed much variation, from an enrichment factor up to 3.2 to worse than random, whereas the structure-based VS gave a more stable result with a constant enrichment factor around 2. Additionally, actives retrieved by the structure-based approach were more diverse than the actives among the top scorers of the similarity searches. Based on these results, we suggest basing a focused library design for a GPCR project on a combination of a ligand-based similarity search and structure-based docking.

A physicogenetic method to assign ligand-binding relationships between 7TM receptors

A computational protocol has been devised to relate 7TM receptor proteins (GPCRs) with respect to physicochemical features of the core ligand-binding site as defined from the crystal structure of bovine rhodopsin. The identification of such receptors that already are associated with ligand information (e.g., small molecule ligands with mutagenesis or SAR data) is used to support structure-guided drug design of novel ligands. A case targeting the newly identified prostaglandin D2 receptor CRTH2 serves as a primary example to illustrate the procedure.

Improved approach for proteochemometrics modeling: application to organic compound--amine G protein-coupled receptor interactions

MOTIVATION

Proteochemometrics is a novel technology for the analysis of interactions of series of proteins with series of ligands. We have here customized it for analysis of large datasets and evaluated it for the modeling of the interaction of psychoactive organic amines with all the five known families of amine G protein-coupled receptors (GPCRs).

RESULTS

The model exploited data for the binding of 22 compounds to 31 amine GPCRs, correlating chemical descriptions and cross-descriptions of compounds and receptors to binding affinity using a novel strategy. A highly valid model (q2 = 0.76) was obtained which was further validated by external predictions using data for 10 other entirely independent compounds, yielding the high q2ext = 0.67. Interpretation of the model reveals molecular interactions that govern psychoactive organic amines overall affinity for amine GPCRs, as well as their selectivity for particular amine GPCRs. The new modeling procedure allows us to obtain fully interpretable proteochemometrics models using essentially unlimited number of ligand and protein descriptors.

New 5-Hydroxytryptamine1A Receptor Ligands Containing a Norbornene Nucleus: Synthesis and in Vitro Pharmacological Evaluation

New arylpiperazine derivatives were prepared to identify highly selective and potent ligands for the 5-hydroxytryptamine 1A (5-HT(1A)) receptor as potential pharmacological tools in studies of central nervous system (CNS) disorders. The combination of structural elements (heterocyclic nucleus, oxyalkyl chain, and arylpiperazine) known to introduce 5-HT(1A) receptor affinity and the proper selection of substituents led to compounds with higher receptor specificity and affinity. In binding studies, several molecules showed affinity in the nanomolar and subnanomolar ranges at 5-HT(1A) and moderate to no affinity for other relevant receptors (5-HT(2A), 5-HT(2C), D(1), D(2), alpha(1), and alpha(2)). The 4-[3-[4-(o-methoxyphenyl)piperazin-1-yl]propoxy]-4-aza-tricyclo[5.2.1.02,6]dec-8-ene-3,5-dione, with K(i) = 0.021 nM, was the most active and selective derivative for the 5-HT(1A) receptor with respect to other serotonin receptors, whereas the most selective derivative for dopaminergic and adrenergic receptors was a CF(3)-substituted arylpiperazine. As a general trend, compounds with a piperazinylpropoxy chain showed a preferential affinity for the 5-HT(1A) receptor, suggesting that the alkyl chain length represents a critical structural feature in determining 5-HT(1A) receptor affinity and selectivity, as confirmed by the molecular modeling invoked for explaining the differential binding affinities of the new arylpiperazines.

A Three-Dimensional Pharmacophore Model for 5-Hydroxytryptamine6 (5-HT6) Receptor Antagonists

Forty-five structurally diverse 5-hydroxytryptamine(6) receptor (5-HT(6)R) antagonists were selected to develop a 3D pharmacophore model with the Catalyst software. The structural features for antagonism at this receptor are a positive ionizable atom interacting with Asp(3.32), a hydrogen bond acceptor group interacting with Ser(5.43) and Asn(6.55), a hydrophobic site interacting with residues in a hydrophobic pocket between transmembranes 3, 4, and 5, and an aromatic-ring hydrophobic site interacting with Phe(6.52).

The serotonin1A receptor: a representative member of the serotonin receptor family

1. Serotonin is an intrinsically fluorescent biogenic amine that acts as a neurotransmitter and is found in a wide variety of sites in the central and peripheral nervous system. Serotonergic signaling appears to play a key role in the generation and modulation of various cognitive and behavioral functions. 2. Serotonin exerts its diverse actions by binding to distinct cell surface receptors which have been classified into many groups. The serotonin1A (5-HT1A) receptor is the most extensively studied of the serotonin receptors and belongs to the large family of seven transmembrane domain G-protein coupled receptors. 3. The tissue and sub-cellular distribution, structural characteristics, signaling of the serotonin1A receptor and its interaction with G-proteins are discussed. 4. The pharmacology of serotonin1A receptors is reviewed in terms of binding of agonists and antagonists and sensitivity of their binding to guanine nucleotides. 5. Membrane biology of 5-HT1A receptors is presented using the bovine hippocampal serotonin1A receptor as a model system. The ligand binding activity and G-protein coupling of the receptor is modulated by membrane cholesterol thereby indicating the requirement of cholesterol in maintaining the receptor organization and function. This, along with the reported detergent resistance characteristics of the receptor, raises important questions on the role of membrane lipids and domains in the function of this receptor.

Binding properties of dipropyltryptamine at the human 5-HT1a receptor

Dipropyltryptamine (DPT) is a synthetic indolealkylamine first characterized in the 1960s. Largely forgotten since the discovery of multiple serotonin receptor subtypes, some of the properties of DPT at the cloned human 5-HT1a receptor are described here. When [3H]8-OH-DPAT is bound to the receptor, DPT inhibits the interaction with an IC50 of 0.1 micromol/l. This interaction is shown to be competitive when double-reciprocal plots of the DPT/agonist interaction are analyzed. DPT's effects in the signal transduction system are complex. While DPT alone (0.1-1,000 micromol/l) activates Gi when both cAMP and gamma-S-GTP incorporation are measured, in the presence of 5-HT (0.1-10 micromol/l), DPT blocks the agonist effect. In combination, the findings suggest that DPT is a moderate affinity partial agonist at the human 5-HT1a receptor. These results provide evidence that DPT has potential as a versatile experimental tool at 5-HT1a receptors.

Benzimidazole derivatives. Part 5: design and synthesis of new benzimidazole-arylpiperazine derivatives acting as mixed 5-HT1A/5-HT3 ligands

A series of new mixed benzimidazole-arylpiperazine derivatives were designed by incorporating in general structure III the pharmacophoric elements of 5-HT(1A) and 5-HT(3) receptors. Compounds 1-11 were synthesized and evaluated for binding affinity at both serotoninergic receptors, all of them exhibiting high 5-HT(3)R affinity (K(i)=10-62nM), and derivatives with an o-alkoxy group in the arylpiperazine ring showing nanomolar affinity for the 5-HT(1A)R (K(i)=18-150nM). Additionally, all the synthesized compounds were selective over alpha(1)-adrenergic and dopamine D(2) receptors (K(i)>1000-10,000nM). Compound 3 was selected for further pharmacological characterization due to its interesting binding profile as mixed 5-HT(1A)/5-HT(3) ligand with high affinity for both receptors (5-HT(1A): K(i)=18.0nM, 5-HT(3): K(i)=27.2nM). In vitro and in vivo findings suggest that this compound acts as a partial agonist at 5-HT(1A)Rs and as a 5-HT(3)R antagonist. This novel mixed 5-HT(1A)/5-HT(3) ligand was also effective in preventing the cognitive deficits induced by muscarinic receptor blockade in a passive avoidance learning test, suggesting a potential interest in the treatment of cognitive dysfunction.

Synthesis and Structure−Activity Relationships of a New Model of Arylpiperazines. 8.1Computational Simulation of Ligand−Receptor Interaction of 5-HT1AR Agonists with Selectivity over α1-Adrenoceptors

We have designed and synthesized a new series of arylpiperazines V exhibiting high 5-HT(1A)R affinity and selectivity over alpha(1)-adrenoceptors. The new selective 5-HT(1A)R ligands contain a hydantoin (m = 0) or diketopiperazine (m = 1) moiety and an arylpiperazine moiety separated by one methylene unit (n = 1). The aryl substituent of the piperazine moiety (Ar) consists of different benzofused rings mimicking the favorable voluminous substituents at ortho and meta positions predicted by 3D-QSAR analysis in the previously reported series I. In particular, (S)-2-[[4-(naphth-1-yl)piperazin-1-yl]methyl]-1,4-dioxoperhydropyrrolo[1,2-a]pyrazine [(S)-9, CSP-2503] (5-HT(1A), K(i) = 4.1 nM; alpha(1), K(i) > 1000 nM) has been pharmacologically characterized as a 5-HT(1A)R agonist at somatodendritic and postsynaptic sites, endowed with anxiolytic properties. Ligand (S)-9 is predicted, in computer simulations, to bind Asp(3.32) in TMH 3, Thr(5.39) and Ser(5.42) in TMH 5, and Trp(6.48) in TMH 6. We propose that agonists modify, by means of an explicit hydrogen bond, the conformation of Trp(6.48) from pointing toward TMH 7, in the inactive gauche+ conformation, to pointing toward the ligand binding site, in the active trans conformation.

Atypical and typical antipsychotic drug interactions with the dopamine D2 receptor

A model of the dopamine D2 receptor was used to study the receptor interactions of dopamine, the typical antipsychotics haloperidol and loxapine, and the atypical antipsychotics clozapine and melperone. The atypical antipsychotics interacted with the halogen atom of the ring system in the direction of the transmembrane helices (TMHs) 2, 3 and 7, while the typical had the corresponding halogen atom in the direction of TMH5. Molecular dynamics simulations indicated that the average helical displacement upon binding increased in the order: typical < atypical < dopamine. Upon binding, the atypical induced larger displacements into TMH5 than did the typical. The typical had stronger non-bonded interactions with the receptor than had the atypical, which is in agreement with the experimental observation that the atypical antipsychotic drugs dissociate faster from the receptor than the typical antipsychotic drugs.

Amino acids 270 to 510 of the severe acute respiratory syndrome coronavirus spike protein are required for interaction with receptor

A novel coronavirus, severe acute respiratory syndrome coronavirus (SARS-CoV), has recently been identified as the causative agent of severe acute respiratory syndrome (SARS). SARS-CoV appears similar to other coronaviruses in both virion structure and genome organization. It is known for other coronaviruses that the spike (S) glycoprotein is required for both viral attachment to permissive cells and for fusion of the viral envelope with the host cell membrane. Here we describe the construction and expression of a soluble codon-optimized SARS-CoV S glycoprotein comprising the first 1,190 amino acids of the native S glycoprotein (S(1190)). The codon-optimized and native S glycoproteins exhibit similar molecular weight as determined by Western blot analysis, indicating that synthetic S glycoprotein is modified correctly in a mammalian expression system. S(1190) binds to the surface of Vero E6 cells, a cell permissive to infection, as demonstrated by fluorescence-activated cell sorter analysis, suggesting that S(1190) maintains the biologic activity present in native S glycoprotein. This interaction is blocked with serum obtained from recovering SARS patients, indicating that the binding is specific. In an effort to map the ligand-binding domain of the SARS-CoV S glycoprotein, carboxy- and amino-terminal truncations of the S(1190) glycoprotein were constructed. Amino acids 270 to 510 were the minimal receptor-binding region of the SARS-CoV S glycoprotein as determined by flow cytometry. We speculate that amino acids 1 to 510 of the SARS-CoV S glycoprotein represent a unique domain containing the receptor-binding site (amino acids 270 to 510), analogous to the S1 subunit of other coronavirus S glycoproteins.

Dopamine D4 Ligands and Models of Receptor Activation: 2-(4-Pyridin-2-ylpiperazin-1-ylmethyl)-1H-benzimidazole and Related Heteroarylmethylarylpiperazines Exhibit a Substituent Effect Responsible for Additional Efficacy Tuning

A series of subtype selective dopamine D(4) receptor ligands from the hetroarylmethylphenylpiperazine class have been discovered that exhibit a remarkable structure-activity relationship (SAR), revealing a substituent effect in which regiosubstitution on the terminal arylpiperazine ring can modulate functional or intrinsic activity. Other structure-dependent efficacy studies in the dopamine D(4) field have suggested a critical interaction of the heteroarylmethyl moiety with specific protein microdomains in controlling intrinsic activity. Our studies indicate that for some binding orientations, the phenylpiperazine moiety also plays a key role in determining efficacy. These data also implicate a kinetic or efficiency term, contained within measured functional affinities for agonists, which support a sequential binding and conformational stabilization model for receptor activation. The structural similarity between partial agonist and antagonist, within this subset of ligands, and lack of bioisosterism for this substituent effect are key phenomena for these hypotheses.

Synthesis and structure–activity relationships of a new model of arylpiperazines. Part 7: Study of the influence of lipophilic factors at the terminal amide fragment on 5-HT1A affinity/selectivity

The influence of lipophilic factors at the amide fragment of a new series of (+/-)-7a-alkyl-2-[4-(4-arylpiperazin-1-yl)butyl]-1,3-dioxoperhydropyrrolo[1,2-c]imidazoles 2 and of (+/-)-7a-alkyl-2-[(4-arylpiperazin-1-yl)methyl]-1,3-dioxoperhydropyrrolo[1,2-c]imidazoles 3 has been studied. Variations of logP have been carried out by introducing different hydrocarbonated substituents (R(1)) at the position 7a of the bicyclohydantoin, namely the non-pharmacophoric part. All the new compounds exhibit high potency for the 5-HT(1A) receptor; however, affinities for the alpha(1) receptor are high for compounds 2a-l while compounds 3a-f are selective over this adrenergic receptor. On the other hand, differences in logP do not notably affect the K(i) values for the above receptors.

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.

Differences in the central nervous system distribution and pharmacology of the mouse 5-hydroxytryptamine-6 receptor compared with rat and human receptors investigated by radioligand binding, site-directed mutagenesis, and molecular modeling

There is increasing evidence for a role of 5-hydroxytrypta-mine-6 (5-HT6) receptors in cognitive function. In the rat and human brain, 5-HT6 receptors are widely expressed and highly enriched in the basal ganglia. However, in the mouse brain, only very low levels of 5-HT6 receptor mRNA and receptor protein, measured by TaqMan reverse transcriptase-polymerase chain reaction and selective radioligand binding, could be detected, with no evidence of enrichment in the basal ganglia. The mouse receptor was cloned and transiently expressed in human embryonic kidney 293 cells to characterize its pharmacological profile. Despite significant sequence homology between human, rat, and mouse 5-HT6 receptors, the pharmacological profile of the mouse receptor was significantly different from the rat and human receptors. Four amino acid residues, conserved in rat and human and divergent in mouse receptors, were identified, and various mutant receptors were generated and their pharmacologies studied. Residues 188 (tyrosine in mouse, phenylalanine in rat and human) in transmembrane region 5 and 290 (serine in mouse, asparagine in rat and human) in transmembrane region 6 were identified as key amino acids responsible for the different pharmacological profiles. Molecular modeling of the receptor and docking of selective and nonselective compounds was undertaken to elucidate the ligand receptor interactions. The binding pocket was predicted to be different in the mouse compared with rat and human 5-HT6 receptors, and the models were in excellent agreement with the observed mutation results and have been used extensively in the design of further selective 5-HT6 antagonists.

Molecular modelling of fentanyl analogs

Fentanyl is a highly potent and clinically widely used narcotic analgesic. A large number of its analogs have been synthesized, some of which (sufentanil and alfentanyl) are also in clinical use. Theoretical studies, in recent years, afforded a better understanding of the structure-activity relationships of this class of opiates and allowed insight into the molecular mechanism of the interactions of fentanyl analogs with their receptors. An overview of the current computational techniques for modeling fentanyl analogs, their receptors and ligand-receptor interactions is presented in this paper.

Homology modeling of the transmembrane domain of the human calcium sensing receptor and localization of an allosteric binding site

A homology model for the human calcium sensing receptor (hCaR) transmembrane domain utilizing bovine rhodopsin (bRho) structural information was derived and tested by docking the allosteric antagonist, NPS 2143, followed by mutagenesis of predicted contact sites. Mutation of residues Phe-668 (helix II), Arg-680, or Phe-684 (helix III) to Ala (or Val or Leu) and Glu-837 (helix VII) to Ile (or Gln) reduced the inhibitory effects of NPS 2143 on [Ca2+]i responses. The calcimimetic NPS R-568 increases the potency of Ca2+ in functional assays of CaR. Mutations at Phe-668, Phe-684, or Glu-837 attenuated the effects of this compound, but mutations at Arg-680 had no effect. In all cases, mutant CaRs responded normally to Ca2+ or phenylalanine, which act at distinct site(s). Discrimination by the Arg-680 mutant is consistent with the structural differences between NPS 2143, which contains an alkyl bridge hydroxyl group, and NPS R-568, which does not. The homology model of the CaR transmembrane domain robustly accounts for binding of both an allosteric antagonist and agonist, which share a common site, and provides a basis for the development of more specific and/or potent allosteric modulators of CaR. These studies suggest that the bRho backbone can be used as a starting point for homology modeling of even distantly related G protein-coupled receptors and provide a rational framework for investigation of the contributions of the transmembrane domain to CaR function.

Three-dimensional models for beta-adrenergic receptor complexes with agonists and antagonists

Molecular modeling methods have been used to construct three-dimensional models for agonist and antagonist complexes with beta-adrenergic receptors. The recent rhodopsin crystal structure was used as a template in standard homology modeling methods. The rhodopsin-based homology models were assessed for agreement with experimental results for beta-adrenergic receptors, and compared with receptor models developed using de novo modeling techniques. While the de novo and homology-derived receptor models are generally quite similar, there are some localized structural differences that impact the putative ligand-binding site significantly. The de novo receptor models appear to provide much better agreement with experimental data, particularly for receptor complexes with agonist ligands. The de novo receptor models also yield some interesting and testable hypotheses for the structural basis of beta-adrenergic receptor subtype ligand selectivity.

Design and synthesis of new benzimidazole-arylpiperazine derivatives acting as mixed 5-HT1A/5-HT3 ligands

A series of new benzimidazole-arylpiperazine derivatives III were designed, synthesized and evaluated for binding affinity at serotoninergic 5-HT(1A) and 5-HT(3) receptors. Compound IIIc was identified as a novel mixed 5-HT(1A)/5-HT(3) ligand with high affinity for both serotonin receptors and excellent selectivity over alpha(1)-adrenergic and dopamine D(2) receptors. This compound was characterized as a partial agonist at 5-HT(1A)Rs and a 5-HT(3)R antagonist, and was effective in preventing the cognitive deficits induced by muscarinic receptor blockade in a passive avoidance learning test.