Contributions of conserved serine residues to the interactions of ligands with dopamine D2 receptors

Evaluating the conformational space of the active site of D2 dopamine receptor. Scope and limitations of the standard docking methods

We report here for the first time the potential energy surfaces (PES) of phenyletilamine (PEA) and meta-tyramine (m-OH-PEA) at the D₂ dopamine receptor (D2DR) binding site. PESs not only allow us to observe all the critical points of the surface (minimums, maximums, and transition states), but also to note the ease or difficulty that each local minima have for their conformational inter-conversions and therefore know the conformational flexibility that these ligands have in their active sites. Taking advantage of possessing this valuable information, we analyze how accurate a standard docking study is in these cases. Our results indicate that although we have to be careful in how to carry out this type of study and to consider performing some extra-simulations, docking calculations can be satisfactory. In order to analyze in detail the different molecular interactions that are stabilizing the different ligand-receptor (L-R) complexes, we carried out quantum theory of atoms in molecules (QTAIM) computations and NMR shielding calculations. Although some of these techniques are a bit tedious and require more computational time, our results demonstrate the importance of performing computational simulations using different types of combined techniques (docking/MD/hybrid QM-MM/QTAIM and NMR shielding calculations) in order to obtain more accurate results. Our results allow us to understand in details the molecular interactions stabilizing and destabilizing the different L-R complexes reported here. Thus, the different activities observed for dopamine (DA), m-OH-PEA, and PEA can be clearly explained at molecular level.

Evidence for Two Modes of Binding of the Negative Allosteric Modulator SB269,652 to the Dopamine D2 Receptor

SB269,652 has been described as the first negative allosteric modulator (NAM) of the dopamine D₂ receptor (D₂R), however, the binding mode and allosteric mechanism of action of this ligand remain incompletely understood. SB269,652 comprises an orthosteric, primary pharmacophore and a secondary (or allosteric) pharmacophore joined by a hydrophilic cyclohexyl linker and is known to form corresponding interactions with the orthosteric binding site (OBS) and the secondary binding pocket (SBP) in the D₂R. Here, we observed a surprisingly low potency of SB269,652 to negatively modulate the D₂R-mediated activation of G protein-coupled inward-rectifier potassium channels (GIRK) and decided to perform a more detailed investigation of the interaction between dopamine and SB269,652. The results indicated that the SB269,652 inhibitory potency is increased 6.6-fold upon ligand pre-incubation, compared to the simultaneous co-application with dopamine. Mutagenesis experiments implicated both S193 in the OBS and E95 in the SBP in the effect of pre-application. The present findings extend previous knowledge about how SB269,652 competes with dopamine at the D₂R and may be useful for the development of novel D₂R ligands, such as antipsychotic drug candidates.

Mechanistic insights into dopaminergic and serotonergic neurotransmission - concerted interactions with helices 5 and 6 drive the functional outcome

Brain functions rely on neurotransmitters that mediate communication between billions of neurons. Disruption of this communication can result in a plethora of psychiatric and neurological disorders. In this work, we combine molecular dynamics simulations, live-cell biosensor and electrophysiological assays to investigate the action of the neurotransmitter dopamine at the dopaminergic D2 receptor (D2R). The study of dopamine and closely related chemical probes reveals how neurotransmitter binding translates into the activation of distinct subsets of D2R effectors (i.e.: Gi2, GoB, Gz and β-arrestin 2). Ligand interactions with key residues in TM5 (S5.42) and TM6 (H6.55) in the D2R binding pocket yield a dopamine-like coupling signature, whereas exclusive TM5 interaction is typically linked to preferential G protein coupling (in particular GoB) over β-arrestin. Further experiments for serotonin receptors indicate that the reported molecular mechanism is shared by other monoaminergic neurotransmitter receptors. Ultimately, our study highlights how sequence variation in position 6.55 is used by nature to fine-tune β-arrestin recruitment and in turn receptor signaling and internalization of neurotransmitter receptors.

Dopamine D2 Receptor Agonist Binding Kinetics-Role of a Conserved Serine Residue

The forward (k_on) and reverse (k_off) rate constants of drug–target interactions have important implications for therapeutic efficacy. Hence, time-resolved assays capable of measuring these binding rate constants may be informative to drug discovery efforts. Here, we used an ion channel activation assay to estimate the k_ons and k_offs of four dopamine D₂ receptor (D₂R) agonists; dopamine (DA), p-tyramine, (R)- and (S)-5-OH-dipropylaminotetralin (DPAT). We further probed the role of the conserved serine S193^5.42 by mutagenesis, taking advantage of the preferential interaction of (S)-, but not (R)-5-OH-DPAT with this residue. Results suggested similar k_offs for the two 5-OH-DPAT enantiomers at wild-type (WT) D₂R, both being slower than the k_offs of DA and p-tyramine. Conversely, the k_on of (S)-5-OH-DPAT was estimated to be higher than that of (R)-5-OH-DPAT, in agreement with the higher potency of the (S)-enantiomer. Furthermore, S193^5.42A mutation lowered the k_on of (S)-5-OH-DPAT and reduced the potency difference between the two 5-OH-DPAT enantiomers. Kinetic K_ds derived from the k_off and k_on estimates correlated well with EC₅₀ values for all four compounds across four orders of magnitude, strengthening the notion that our assay captured meaningful information about binding kinetics. The approach presented here may thus prove valuable for characterizing D₂R agonist candidate drugs.

Structure of the dopamine D2 receptor in complex with the antipsychotic drug spiperone

In addition to the serotonin 5-HT_2A receptor (5-HT_2AR), the dopamine D₂ receptor (D₂R) is a key therapeutic target of antipsychotics for the treatment of schizophrenia. The inactive state structures of D₂R have been described in complex with the inverse agonists risperidone (D₂R_ris) and haloperidol (D₂R_hal). Here we describe the structure of human D₂R in complex with spiperone (D₂R_spi). In D₂R_spi, the conformation of the extracellular loop (ECL) 2, which composes the ligand-binding pocket, was substantially different from those in D₂R_ris and D₂R_hal, demonstrating that ECL2 in D₂R is highly dynamic. Moreover, D₂R_spi exhibited an extended binding pocket to accommodate spiperone’s phenyl ring, which probably contributes to the selectivity of spiperone to D₂R and 5-HT_2AR. Together with D₂R_ris and D₂R_hal, the structural information of D₂R_spi should be of value for designing novel antipsychotics with improved safety and efficacy.

The dopamine D₂ receptor (D₂R) is a GPCR and an important drug target for schizophrenia treatment. Here, the authors present the crystal structure of human D₂R in complex with the antipsychotic drug spiperone, which is of interest for designing antipsychotics with improved receptor selectivity.

Combined MD/QTAIM techniques to evaluate ligand-receptor interactions. Scope and limitations

In medicinal chemistry, it is extremely important to evaluate, as accurately as possible, the molecular interactions involved in the formation of different ligand-receptor (L-R) complexes. Evaluating the different molecular interactions by quantum mechanics calculations is not a simple task, since formation of an L-R complex is a dynamic process. In this case, the use of combined techniques of molecular dynamics (MD) and quantum calculations is one the best possible approaches. In this work we report a comparative study using combined MD and QTAIM (Quantum Theory of Atoms In Molecules) calculations for five biological systems with different levels of structural complexity. We have studied Acetylcholinesterase (AChE), D2 Dopamine Receptor (D2DR), beta Secretase (BACE1), Dihydrofolate Reductase (DHFR) and Sphingosine Kinase 1 (SphK1). In these molecular targets, we have analyzed different ligands with diverse structural characteristics. The inhibitory activities of most of them have been previously measured in our laboratory. Our results indicate that QTAIM calculations can be extremely useful for in silico studies. It is possible to obtain very accurate information about the strength of the molecular interactions that stabilize the formation of the different L-R complexes. Better correlations can be obtained between theoretical and experimental data by using QTAIM calculations, allowing us to discriminate among ligands with similar affinities. QTAIM analysis gives fairly accurate information for weak interactions which are not well described by MD simulations. QTAIM study also allowed us to evaluate and determine which parts of the ligand need to be modified in order to increase its interactions with the molecular target. In this study we have discussed the importance of combined MD/QTAIM calculations for this type of simulations, showing their scopes and limitations.

Crystal structure of dopamine receptor D4 bound to the subtype selective ligand, L745870

Multiple subtypes of dopamine receptors within the GPCR superfamily regulate neurological processes through various downstream signaling pathways. A crucial question about the dopamine receptor family is what structural features determine the subtype-selectivity of potential drugs. Here, we report the 3.5-angstrom crystal structure of mouse dopamine receptor D4 (DRD4) complexed with a subtype-selective antagonist, L745870. Our structure reveals a secondary binding pocket extended from the orthosteric ligand-binding pocket to a DRD4-specific crevice located between transmembrane helices 2 and 3. Additional mutagenesis studies suggest that the antagonist L745870 prevents DRD4 activation by blocking the relative movement between transmembrane helices 2 and 3. These results expand our knowledge of the molecular basis for the physiological functions of DRD4 and assist new drug design.

Aminergic GPCR-Ligand Interactions: A Chemical and Structural Map of Receptor Mutation Data

The aminergic family of G protein-coupled receptors (GPCRs) plays an important role in various diseases and represents a major drug discovery target class. Structure determination of all major aminergic subfamilies has enabled structure-based ligand design for these receptors. Site-directed mutagenesis data provides an invaluable complementary source of information for elucidating the structural determinants of binding of different ligand chemotypes. The current study provides a comparative analysis of 6692 mutation data points on 34 aminergic GPCR subtypes, covering the chemical space of 540 unique ligands from mutagenesis experiments and information from experimentally determined structures of 52 distinct aminergic receptor-ligand complexes. The integrated analysis enables detailed investigation of structural receptor-ligand interactions and assessment of the transferability of combined binding mode and mutation data across ligand chemotypes and receptor subtypes. An overview is provided of the possibilities and limitations of using mutation data to guide the design of novel aminergic receptor ligands.

The structural determinants of the bitopic binding mode of a negative allosteric modulator of the dopamine D2 receptor

SB269652 is a negative allosteric modulator of the dopamine D₂ receptor (D₂R) yet possesses structural similarity to ligands with a competitive mode of interaction. In this study, we aimed to understand the ligand-receptor interactions that confer its allosteric action. We combined site-directed mutagenesis with molecular dynamics simulations using both SB269652 and derivatives from our previous structure activity studies. We identify residues within the conserved orthosteric binding site (OBS) and a secondary binding pocket (SBP) that determine affinity and cooperativity. Our results indicate that interaction with the SBP is a requirement for allosteric pharmacology, but that both competitive and allosteric derivatives of SB269652 can display sensitivity to the mutation of a glutamate residue (E95^2.65) within the SBP. Our findings provide the molecular basis for the differences in affinity between SB269652 derivatives, and reveal how changes to interactions made by the primary pharmacophore of SB269652 in the orthosteric pocket can confer changes in the interactions made by the secondary pharmacophore in the SBP. Our insights provide a structure-activity framework towards rational optimization of bitopic ligands for D₂R with tailored competitive versus allosteric properties.

Critical analysis of the successes and failures of homology models of G protein-coupled receptors

We present a critical assessment of the performance of our homology model refinement method for G protein-coupled receptors (GPCRs), called LITICon that led to top ranking structures in a recent structure prediction assessment GPCRDOCK2010. GPCRs form the largest class of drug targets for which only a few crystal structures are currently available. Therefore, accurate homology models are essential for drug design in these receptors. We submitted five models each for human chemokine CXCR4 (bound to small molecule IT1t and peptide CVX15) and dopamine D3DR (bound to small molecule eticlopride) before the crystal structures were published. Our models in both CXCR4/IT1t and D3/eticlopride assessments were ranked first and second, respectively, by ligand RMSD to the crystal structures. For both receptors, we developed two types of protein models: homology models based on known GPCR crystal structures, and ab initio models based on the prediction method MembStruk. The homology-based models compared better to the crystal structures than the ab initio models. However, a robust refinement procedure for obtaining high accuracy structures is needed. We demonstrate that optimization of the helical tilt, rotation, and translation is vital for GPCR homology model refinement. As a proof of concept, our in-house refinement program LITiCon captured the distinct orientation of TM2 in CXCR4, which differs from that of adrenoreceptors. These findings would be critical for refining GPCR homology models in future.

Development of 7TM receptor-ligand complex models using ligand-biased, semi-empirical helix-bundle repacking in torsion space: application to the agonist interaction of the human dopamine D2 receptor

Prediction of 3D structures of membrane proteins, and of G-protein coupled receptors (GPCRs) in particular, is motivated by their importance in biological systems and the difficulties associated with experimental structure determination. In the present study, a novel method for the prediction of 3D structures of the membrane-embedded region of helical membrane proteins is presented. A large pool of candidate models are produced by repacking of the helices of a homology model using Monte Carlo sampling in torsion space, followed by ranking based on their geometric and ligand-binding properties. The trajectory is directed by weak initial restraints to orient helices towards the original model to improve computation efficiency, and by a ligand to guide the receptor towards a chosen conformational state. The method was validated by construction of the β₁ adrenergic receptor model in complex with (S)-cyanopindolol using bovine rhodopsin as template. In addition, models of the dopamine D₂ receptor were produced with the selective and rigid agonist (R)-N-propylapomorphine ((R)-NPA) present. A second quality assessment was implemented by evaluating the results from docking of a library of 29 ligands with known activity, which further discriminated between receptor models. Agonist binding and recognition by the dopamine D₂ receptor is interpreted using the 3D structure model resulting from the approach. This method has a potential for modeling of all types of helical transmembrane proteins for which a structural template with sequence homology sufficient for homology modeling is not available or is in an incorrect conformational state, but for which sufficient empirical information is accessible.

Cloning and characterization of the dopamine like receptor in the oyster Crassostrea angulata: expression during the ovarian cycle

We cloned and characterized a complete cDNA encoding a dopamine receptor (DAR) named Ca-DA1R from Fujian oyster, Crassostrea angulata. The 2843 bp long cDNA sequence includes a 916-bp 5'-UTR, the 1197 bp ORF which encodes a putative protein of 399 amino acids, and a 729 bp 3'-UTR. The Ca-DA1R sequence possesses typical characteristics of a D1 receptor: two main features being a short third intracellular loop and a long inner COOH-terminal tail domain. Using a real-time PCR approach, expression profiles of Ca-DA1R were analyzed in adult tissues and during the four stages of ovarian development. Ca-DA1R was expressed ubiquitously, although transcript levels varied between tissues, with higher mRNA levels detected in the ovary, labial palps and mantle. During the four stages of ovarian development, Ca-DA1R mRNA expression level was higher in the proliferation stage than in the other three stages during the ovary cycle. In situ hybridization results reveal that the Ca-DA1R mRNA is mainly expressed in the epithelium of the gonoducts. These observations suggest that Ca-DA1R binding of DA probably plays an important role in early ovarian development and via regulating oocyte locomotion cooperates with the 5-HT receptor system during the ovarian cycle in C. angulata.

Novel bivalent ligands for D2/D3 dopamine receptors: Significant co-operative gain in D2 affinity and potency

This report describes development of a series of novel bivalent molecules with a pharmacophore derived from the D2/D3 agonist 5-OH-DPAT. Spacer length in the bivalent compounds had a pronounced influence on affinity for D2 receptors. A 23-fold increase of D2 affinity was observed at a spacer length of 9 or 10 (compounds 11d and 14b) compared to monovalent 5-OH-DPAT (Ki; 2.5 and 2.0 vs. 59 nM for 11d and 14b vs. 5-OH-DPAT, respectively). Functional potency of 11d and 14b indicated a 24- and 94-fold increase in potency at the D2 receptor compared to 5-OH-DPAT (EC50; 1.7 and 0.44 vs. 41 nM for 11d and 14b vs. 5-OH-DPAT, respectively). These are the most potent bivalent agonists for D2 receptor known to date. This synergism is consonant with cooperative interaction at the two orthosteric binding sites in the homodimeric receptor.

Receptor conformations involved in dopamine D(2L) receptor functional selectivity induced by selected transmembrane-5 serine mutations

Although functional selectivity is now widely accepted, the molecular basis is poorly understood. We have studied how aspects of transmembrane region 5 (TM5) of the dopamine D(2L) receptor interacts with three rationally selected rigid ligands (dihydrexidine, dinapsoline, and dinoxyline) and the reference compounds dopamine and quinpirole. As was expected from homology modeling, mutation of three TM5 serine residues to alanine (S5.42A, S5.43A, S5.46A) had little effect on antagonist affinity. All three mutations decreased the affinity of the agonist ligands to different degrees, S5.46A being somewhat less affected. Four functions [adenylate cyclase (AC), extracellular signal-regulated kinase 1/2 phosphorylation (MAPK), arachidonic acid release (AA), and guanosine 5'-O-(3-thio)triphosphate binding (GTPγS)] were assessed. The intrinsic activity (IA) of quinpirole was unaffected by any of the mutations, whereas S5.42A and S5.46A mutations abolished the activity of dopamine and the three rigid ligands, although dihydrexidine retained IA at MAPK function only with S5.42A. Remarkably, S5.43A did not markedly affect IA for AC and MAPK for any of the ligands and eliminated AA activity for dinapsoline and dihydrexidine but not dinoxyline. These data suggest that this mutation did not disrupt the overall conformation or signaling ability of the mutant receptors but differentially affected ligand activation. Computational studies indicate that these D(2) agonists stabilize multiple receptor conformations. This has led to models showing the stabilized conformations and interhelical and receptor-ligand contacts corresponding to the different activation pathways stabilized by various agonists. These data provide a basis for understanding D(2L) functional selectivity and rationally discovering functionally selective D(2) drugs.

Investigation of D₂ receptor-agonist interactions using a combination of pharmacophore and receptor homology modeling

A combined modeling approach was used to identify structural factors that underlie the structure–activity relationships (SARs) of full dopamine D₂ receptor agonists and structurally similar inactive compounds. A 3D structural model of the dopamine D₂ receptor was constructed, with the agonist (−)-(R)-2-OH-NPA present in the binding site during the modeling procedure. The 3D model was evaluated and compared with our previously published D₂ agonist pharmacophore model. The comparison revealed an inconsistency between the projected hydrogen bonding feature (Ser-TM5) in the pharmacophore model and the TM5 region in the structure model. A new refined pharmacophore model was developed, guided by the shape of the binding site in the receptor model and with less emphasis on TM5 interactions. The combination of receptor and pharmacophore modeling also identified the importance of His393^6.55 for agonist binding. This convergent 3D pharmacophore and protein structure modeling strategy is considered to be general and can be highly useful in less well-characterized systems to explore ligand–receptor interactions. The strategy has the potential to identify weaknesses in the individual models and thereby provides an opportunity to improve the discriminating predictivity of both pharmacophore searches and structure-based virtual screens.

Agonist-specific voltage sensitivity at the dopamine D2S receptor--molecular determinants and relevance to therapeutic ligands

Voltage sensitivity has been demonstrated for some GPCRs. At the dopamine D(2S) receptor, this voltage sensitivity is agonist-specific; some agonists, including dopamine, exhibit decreased potency at depolarized potentials, whereas others are not significantly affected. In the present study, we examined some of the receptor-agonist interactions contributing to these differences, and investigated how dopamine D(2S) receptor voltage sensitivity affects clinically used dopamine agonists. GIRK channel activation in voltage-clamped Xenopus oocytes was used as readout of receptor activation. Structurally distinct agonists and complementary site-directed mutagenesis of the receptor's binding site were used to investigate the role of agonist-receptor interactions. We also confirmed that the depolarization-induced decrease of dopamine potency in GIRK activation is correlated by decreased binding of radiolabeled dopamine, and by decreased potency in G protein activation. In the mutagenesis experiments, a conserved serine residue as well as the conserved aspartate in the receptor's binding site were found to be important for voltage sensitive potency of dopamine. Furthermore, the voltage sensitivity of the receptor had distinct effects on different therapeutic D(2) agonists. Depolarization decreased the potency of several compounds, whereas for others, efficacy was reduced. For some agonists, both potency and efficacy were diminished, whereas for others still, neither parameter was significantly altered. The present work identifies some of the ligand-receptor interactions which determine agonist-specific effects of voltage at the dopamine D(2S) receptor. The observed differences between therapeutic agonists might be clinically relevant, and make them potential tools for investigating the roles of dopamine D(2) receptor voltage sensitivity in native tissue.

Interaction of novel hybrid compounds with the D3 dopamine receptor: Site-directed mutagenesis and homology modeling studies

The dopamine D3 receptor has been implicated as a potential target for drug development in various complex psychiatric disorders including psychosis, drug dependence, and Parkinson's disease. In our overall goal to develop molecules with preferential affinity at D3 receptors, we undertook a hybrid drug development approach by combining a known dopamine agonist moiety with a substituted piperazine fragment. In the present study, three compounds produced this way with preferential D3 agonist activity, were tested at D3 receptors with mutations in the agonist binding pocket of three residues known to be important for agonist binding activity. At S192A and T369V, the hybrid agonist compounds produced an interaction profile in [(3)H]spiperone binding assays similar to that of the parent 5-OH-DPAT and 7-OH-DPAT molecules. The loss of affinity at the S192A mutant was most prominent for 5-OH-DPAT and its corresponding hybrid compound D237. D110N did not show any radioligand binding. Homology modeling indicated that 7-OH-DPAT-derived D315 uniquely shares H-bonding with Tyr365 which produced favorable interaction and no loss of H-bonding in the S192A mutant, suggesting that agonist activity may not be solely controlled by residues in the binding pocket.

Dopamine D2 receptors as treatment targets in schizophrenia

The antipsychotic effectiveness of chlorpromazine and haloperidol started a search for their therapeutic targets. The antipsychotic receptor target turned out to be a dopamine receptor, now cloned as the dopamine D2 receptor. The D2 receptor is the common target for antipsychotics. Antipsychotic clinical doses correlate with their affinities for this receptor. Therapeutic doses of antipsychotics occupy 60 to 80% of brain D2 receptors in patients, but aripiprazole occupies up to 90%. While antipsychotics may take up to six hours to occupy D2 receptors, much clinical improvement occurs within a few days. The receptor has high- and low-affinity states. The D2High state is functional for dopamine-like agonists such as aripiprazole. Most individuals with schizophrenia are supersensitive to dopamine. Animal models of psychosis show that a variety of risk factors, genetic and nongenetic, are associated with behavioral supersensitivity to dopamine, reflected in elevated levels of dopamine D2High receptors. Although antipsychotics such as haloperidol alleviate psychosis and reverse the elevation of D2High receptors, long-term use of traditional antipsychotics can further enhance dopamine supersensitivity in patients. Therefore, switching from a traditional antipsychotic to an agonist antipsychotic such as aripiprazole can result in the emergence of psychotic signs and symptoms. Clozapine and quetiapine do not elicit parkinsonism and rarely result in tardive dyskinesia because they are released from D2 within 12 to 24 hours. Traditional antipsychotics remain attached to D2 receptors for days, preventing relapse, but allowing accumulation that can lead to tardive dyskinesia. Future goals include imaging D2High receptors and desensitizing them in early-stage psychosis.

The dopamine D4 receptor: biochemical and signalling properties

Dopamine is an important neurotransmitter that regulates several key functions in the brain, such as motor output, motivation and reward, learning and memory, and endocrine regulation. Dopamine does not mediate fast synaptic transmission, but rather modulates it by triggering slow-acting effects through the activation of dopamine receptors, which belong to the G-protein-coupled receptor superfamily. Besides activating different effectors through G-protein coupling, dopamine receptors also signal through interaction with a variety of proteins, collectively termed dopamine receptor-interacting proteins. We focus on the dopamine D4 receptor, which contains an important polymorphism in its third intracellular loop. This polymorphism has been the subject of numerous studies investigating links with several brain disorders, such as attention-deficit hyperactivity disorder and schizophrenia. We provide an overview of the structure, signalling properties and regulation of dopamine D4 receptors, and briefly discuss their physiological and pathophysiological role in the brain.

Transmembrane segment five serines of the D4 dopamine receptor uniquely influence the interactions of dopamine, norepinephrine, and Ro10-4548

Conserved serines of transmembrane segment (TM) five (TM5) are critical for the interactions of endogenous catecholamines with alpha(1)- and alpha(2)-adrenergic, beta(2)-adrenergic, and D1, D2, and D3 dopamine receptors. The unique high-affinity interaction of the D4 dopamine receptor subtype with both norepinephrine and dopamine, and the fact that TM5 serine interactions have never been studied for this receptor subtype, led us to investigate the interactions of ligands with D4 receptor TM5 serines. Serine-to-alanine mutations at positions 5.42 and 5.46 drastically decreased affinities of dopamine and norepinephrine for the D4 receptor. The D4-S5.43A receptor mutant had substantially reduced affinity for norepinephrine, but a modest loss of affinity for dopamine. In functional assays of cAMP accumulation, norephinephrine was unable to activate any of the mutant receptors, even though the agonist quinpirole displayed wild-type functional properties for all of them. Dopamine was unable to activate the S5.46A mutant and had reduced potency for the S5.43A mutant and reduced potency and efficacy for the S5.42A mutant. In contrast, Ro10-4548 [RAC-2'-2-hydroxy-3-4-(4-hydroxy-2-methoxyphenyl)-1-piperazinyl-propoxy-acetanilide], a catechol-like antagonist of the wild-type receptor unexpectedly functions as an agonist of the S5.43A mutant. Other noncatechol ligands had similar properties for mutant and wild-type receptors. This is the first example of a dopamine receptor point mutation selectively changing the receptor's interaction with a specific antagonist to that of an agonist, and together with other data, provides evidence, supported by molecular modeling, that catecholamine-type agonism is induced by different ligand-specific configurations of intermolecular H-bonds with the TM5 conserved serines.

Synthesis and evaluation of a set of 4-phenylpiperidines and 4-phenylpiperazines as D2 receptor ligands and the discovery of the dopaminergic stabilizer 4-[3-(methylsulfonyl)phenyl]-1-propylpiperidine (huntexil, pridopidine, ACR16)

Modification of the partial dopamine type 2 receptor (D(2)) agonist 3-(1-benzylpiperidin-4-yl)phenol (9a) generated a series of novel functional D(2) antagonists with fast-off kinetic properties. A representative of this series, pridopidine (4-[3-(methylsulfonyl)phenyl]-1-propylpiperidine; ACR16, 12b), bound competitively with low affinity to D(2) in vitro, without displaying properties essential for interaction with D(2) in the inactive state, thereby allowing receptors to rapidly regain responsiveness. In vivo, neurochemical effects of 12b were similar to those of D(2) antagonists, and in a model of locomotor hyperactivity, 12b dose-dependently reduced activity. In contrast to classic D(2) antagonists, 12b increased spontaneous locomotor activity in partly habituated animals. The "agonist-like" kinetic profile of 12b, combined with its lack of intrinsic activity, induces a functional state-dependent D(2) antagonism that can vary with local, real-time dopamine concentration fluctuations around distinct receptor populations. These properties may contribute to its unique "dopaminergic stabilizer" characteristics, differentiating 12b from D(2) antagonists and partial D(2) agonists.

Selective pharmacophore models of dopamine D(1) and D(2) full agonists based on extended pharmacophore features

This study is focused on the identification of structural features that determine the selectivity of dopamine receptor agonists toward D(1) and D(2) receptors. Selective pharmacophore models were developed for both receptors. The models were built by using projected pharmacophoric features that represent the main agonist interaction sites in the receptor (the Ser residues in TM5 and the Asp in TM3), a directional aromatic feature in the ligand, a feature with large positional tolerance representing the positively charged nitrogen in the ligand, and sets of excluded volumes reflecting the shapes of the receptors. The sets of D(1) and D(2) ligands used for modeling were carefully selected from published sources and consist of structurally diverse, conformationally rigid full agonists as active ligands together with structurally related inactives. The robustness of the models in discriminating actives from inactives was tested against four ensembles of conformations generated by using different established methods and different force fields. The reasons for the selectivity can be attributed to both geometrical differences in the arrangement of the features, e.g., different tilt angels of the pi system, as well as shape differences covered by the different sets of excluded volumes. This work provides useful information for the design of new D(1) and D(2) agonists and also for comparative homology modeling of D(1) and D(2) receptors. The approach is general and could therefore be applied to other ligand-protein interactions for which no experimental protein structure is available.

Investigation of a dopamine receptor in Schistosoma mansoni: functional studies and immunolocalization

A dopamine receptor (SmD2) was cloned from adult Schistosoma mansoni. The receptor has the classical heptahelical topology of class A (rhodopsin-like) G protein-coupled receptors (GPCR) and shares sequence homology with D2-like receptors from other species. The full length SmD2 cDNA was expressed in the yeast Saccharomyces cerevisiae and mammalian HEK293 cells. Functional assays in both expression systems revealed that SmD2 was responsive to dopamine in a dose-dependent manner, whereas other structurally related amines had no effect. Activation of SmD2 in mammalian cells caused an elevation in intracellular cAMP but not calcium, suggesting that the receptor coupled to Gs and the stimulation of adenylate cyclase. Pharmacological studies showed that the S. mansoni dopamine receptor was inhibited by apomorphine, a classical dopamine agonist, as well as known dopaminergic antagonists, including chlorpromazine, spiperone and haloperidol. SmD2 immunoreactivity was detected in membrane protein fractions of S. mansoni cercaria, in vitro transformed schistosomula and adult parasites, using a specific peptide antibody. When tested by confocal immunofluorescence, SmD2 was detected in the subtegumental somatic musculature and acetabulum of all larval stages tested. In the adults, SmD2 was enriched in the somatic muscles and, to a lesser extent, the muscular lining of the caecum. The results suggest that SmD2 is an important component of the neuromuscular system in schistosomes.

Potential usefulness of D2R reporter gene imaging by IBF as gene therapy monitoring for cerebellar neurodegenerative diseases

We investigated a gene expression imaging method to examine the level of therapeutic gene expression in the cerebellum. Using a human immunodeficiency virus derived lentivial vector, we expressed the dopamine D(2) receptor (D(2)R) as a reporter protein to mouse cerebellar Purkinje cells. Biodistribution and ex vivo autoradiography studies were performed by giving [(125)I]5-iodo-7-N-[(1-ethyl-2-pyrrolidinyl)methyl]carboxamide-2,3-dihydrobenzofuran ([(125)I]IBF) (1.85 MBq), as a radioactive D(2)R ligand, to model mice expressing the D(2)R with an HA tag (HA-D(2)R) in the cerebellum. In this study, [(125)I]IBF was bound to the D(2)R expressed in the cerebellum of the model mice selectively. Immunostaining was performed to confirm the HA-D(2)R expression in the cerebellum of the model mice. A significant correlation (r=0.900, P<0.001) between areas that expressed HA-D(2)R by immunostaining and areas in which [(125)I]IBF accumulated by the ex vivo autoradiograms was found. These results indicated that radioiodinated IBF is useful as a reporter probe to detect D(2)R reporter gene expression, which can be used for monitoring therapeutic gene expression in the cerebellum.

Structurally constrained hybrid derivatives containing octahydrobenzo[g or f]quinoline moieties for dopamine D2 and D3 receptors: binding characterization at D2/D3 receptors and elucidation of a pharmacophore model

A series of structurally constrained analogues based on hybrid compounds containing octahydrobenzo[g or f]quinoline moieties were designed, synthesized, and characterized for their binding to dopamine D2 and D3 receptors expressed in HEK-293 cells. Among the newly developed constrained molecules, trans-octahydrobenzo[f]quinolin-7-ol (8) exhibited the highest affinity for D2 and D3 receptors, the (-)-isomer being the eutomer. Interestingly, this hybrid constrained version 8 showed significant affinity over the corresponding nonhybrid version 1 (representing a constrained version of the aminotetralin structure only) when assayed under same conditions (K(i) of 49.1 and 14.9 nM for 8 vs 380 and 96.0 nM for 1 at D2 and D3, respectively). Similar results were found with other lead hybrid compounds, indicating a contribution of the piperazine moiety in the observed enhanced affinity. On the basis of the data of new lead constrained derivatives and other lead hybrid derivatives developed by us, a unique pharmacophore model was proposed consisting of three pharmacophoric centers, two with aromatic/hydrophobic and one with cationic features.

Voltage-sensitivity at the human dopamine D2S receptor is agonist-specific

Recently, we and others have shown that agonist potencies at some, but not all, G protein-coupled receptors are voltage-sensitive. Several of those studies employed electrophysiology assays in Xenopus oocytes with G protein-coupled potassium channels as a readout. Using this assay, we have now obtained evidence that voltage-sensitivity at the dopamine D(2S) receptor is agonist-specific. Whereas the potency of dopamine at the D(2S) receptor is decreased by depolarization, the potencies of beta-phenethylamine, p- and m-tyramine are voltage-insensitive. Furthermore, both monohydroxylated and non-hydroxylated N,N-dipropyl-2-aminotetralin compounds are voltage-sensitive. Differential activation of G protein subtypes or differential ratios between effector and active G protein do not underlie this agonist-selective voltage-sensitivity. This is the first demonstration of voltage-sensitive and voltage-insensitive behaviour of different agonists acting via the same receptor.

Binding of serotonin and N1-benzenesulfonyltryptamine-related analogs at human 5-HT6 serotonin receptors: receptor modeling studies

A population of 100 graphics models of the human 5-HT6 serotonin receptor was constructed based on the structure of bovine rhodopsin. The endogenous tryptamine-based agonist serotonin (5-HT; 1) and the benzenesulfonyl-containing tryptamine-derived 5-HT6 receptor antagonist MS-245 (4a) were automatically docked with each of the 100 receptor models using a genetic algorithm approach. Similar studies were conducted with the more selective 5-HT6 receptor agonist EMDT (5) and optical isomers of EMDT-related analog 8, as well as with optical isomers of MS-245 (4a)-related and benzenesulfonyl-containing pyrrolidine 6 and aminotetralin 7. Although associated with the same general aromatic/hydrophobic binding cluster, 5-HT (1) and MS-245 (4a) were found to preferentially bind with distinct receptor conformations, and did so with different binding orientations (i.e., poses). A 5-HT pose/model was found to be common to EMDT (5) and its analogs, whereas that identified for MS-245 (4a) was found common to benzenesulfonyl-containing compounds. Specific amino acid residues were identified that can participate in binding, and evaluation of a sulfenamide analog of MS-245 indicates for the first time that the presence of the sulfonyl oxygen atoms enhances receptor affinity. The results indicate that the presence or absence of an N1-benzenesulfonyl group is a major determinant of the manner in which tryptamine-related agents bind at 5-HT6 serotonin receptors.

Dopamine D1 receptor agonist and D2 receptor antagonist effects of the natural product (-)-stepholidine: molecular modeling and dynamics simulations

(-)-Stepholidine (SPD), an active ingredient of the Chinese herb Stephania, is the first compound found to have dual function as a dopamine receptor D1 agonist and D2 antagonist. Insights into dynamical behaviors of D1 and D2 receptors and their interaction modes with SPD are crucial in understanding the structural and functional characteristics of dopamine receptors. In this study a computational approach, integrating protein structure prediction, automated molecular docking, and molecular dynamics simulations were employed to investigate the dual action mechanism of SPD on the D1 and D2 receptors, with the eventual aim to develop new drugs for treating diseases affecting the central nervous system such as schizophrenia. The dynamics simulations revealed the surface features of the electrostatic potentials and the conformational "open-closed" process of the binding entrances of two dopamine receptors. Potential binding conformations of D1 and D2 receptors were obtained, and the D1-SPD and D2-SPD complexes were generated, which are in good agreement with most of experimental data. The D1-SPD structure shows that the K-167_EL-2-E-302_EL-3 (EL-2: extracellular loop 2; EL-3: extracellular loop 3) salt bridge plays an important role for both the conformational change of the extracellular domain and the binding of SPD. Based on our modeling and simulations, we proposed a mechanism of the dual action of SPD and a subsequent signal transduction model. Further mutagenesis and biophysical experiments are needed to test and improve our proposed dual action mechanism of SPD and signal transduction model.

Structural determinants of pharmacological specificity between D(1) and D(2) dopamine receptors

To test the hypothesis that pharmacological differentiation between D(1) and D(2) dopamine receptors results from interactions of selective ligands with nonconserved residues lining the binding pocket, we mutated amino acid residues in the D(2) receptor to the corresponding aligned residues in the D(1) receptor and vice versa and expressed the receptors in human embryonic kidney 293 cells. Determinations of the affinity of the 14 mutant D(2) receptors and 11 mutant D(1) receptors for D(1)- and D(2)-selective antagonists, and rhodopsin-based homology models of the two receptors, identified two residues whose direct interactions with certain ligands probably contribute to ligand selectivity. The D(1) receptor mutant W99(3.28)F showed dramatically increased affinity for several D(2)-selective antagonists, particularly spiperone (225-fold), whereas the D(2) receptor mutant Y417(7.43)W had greatly decreased affinity for benzamide ligands such as raclopride (200-fold) and sulpiride (125-fold). The binding of the D(1)-selective ligand R-(+)-7-chloro-8-hydroxy-3-methyl-1-phenyl-2,3,4,5-tetrahydro-1H-3-benzazepine (SCH23390) was unaffected, indicating that SCH23390 makes little contact with these ancillary pocket residues. Mutation of A/V(5.39) caused modest but consistent and reciprocal changes in affinity of the receptors for D(1) and D(2)-selective ligands, perhaps reflecting altered packing of the interface of helices 5 and 6. We also obtained some evidence that residues in the second extracellular loop contribute to ligand binding. We conclude that additional determinants of D(1)/D(2) receptor-selective binding are located either in that loop or in the transmembrane helices but, like residue 5.39, indirectly influence the interactions of selective ligands with conserved residues by altering the shape of the primary and ancillary binding pockets.

Characterization of a novel D2-like dopamine receptor with a truncated splice variant and a D1-like dopamine receptor unique to invertebrates from Caenorhabditis elegans

We have cloned two novel Caenorhabditis elegans dopamine receptors, DOP-3 and DOP-4. DOP-3 shows high sequence homology with other D2-like dopamine receptors. As a result of alternative splicing, a truncated splice variant of DOP-3, DOP-3nf, was produced. Because of the in-frame insertion of a stop codon in the third intracellular loop, DOP-3nf lacks the sixth and seventh transmembrane domains that are found in the full-length DOP-3 receptor. Reporter gene assay showed that DOP-3 attenuates forskolin-stimulated cAMP formation in response to dopamine stimulation, whereas DOP-3nf does not. When DOP-3 was coexpressed with DOP-3nf, the ability to inhibit forskolin-stimulated cAMP formation was reduced. DOP-4 shows high sequence homology with D1-like dopamine receptors unique to invertebrates, which are distinct from mammalian D1-like dopamine receptors. Reporter gene assay showed that DOP-4 stimulates cAMP accumulation in response to dopamine stimulation. These two receptors provide new opportunities to understand dopaminergic signaling at the molecular level.

Dopamine receptor microdomains involved in molecular recognition and the regulation of drug affinity and function

A cationic protonatable amine moiety on dopaminergic ligands forms a high affinity reinforced ionic bond with an anionic aspartic acid at position 3.32 of dopamine receptors. When present, catechol hydroxyls of the ligands form hydrogen bonds with serines at position 5.42, 5.43, and 5.46, and this network of hydrogen bonds serves to orient ligands in the binding-site crevice and increase their binding affinity. A steric clash between aromatic moieties of the ligands and aromatic amino acids of the receptor (e.g., H6.55, F6.52 or F6.51 and W6.48) is likely to be propagated in domino-like fashion along the length of TM6, which is believed to trigger activation of the receptor. Specifically, it is the change in the conformation of W6.48 from an orientation perpendicular to the plane of the lipid membrane to one that is parallel that is believed to result in activation. Molecular determinants that mediate the D4/D2-selectivity of many extremely D4-selective 1,4-DAP ligands, include a nonconserved cluster of bulky amino acids at the TM2/TM3 interface (positions 2.61, 3.28 and 3.29).

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.

PREDICT modeling and in-silico screening for G-protein coupled receptors

G-protein coupled receptors (GPCRs) are a major group of drug targets for which only one x-ray structure is known (the nondrugable rhodopsin), limiting the application of structure-based drug discovery to GPCRs. In this paper we present the details of PREDICT, a new algorithmic approach for modeling the 3D structure of GPCRs without relying on homology to rhodopsin. PREDICT, which focuses on the transmembrane domain of GPCRs, starts from the primary sequence of the receptor, simultaneously optimizing multiple 'decoy' conformations of the protein in order to find its most stable structure, culminating in a virtual receptor-ligand complex. In this paper we present a comprehensive analysis of three PREDICT models for the dopamine D2, neurokinin NK1, and neuropeptide Y Y1 receptors. A shorter discussion of the CCR3 receptor model is also included. All models were found to be in good agreement with a large body of experimental data. The quality of the PREDICT models, at least for drug discovery purposes, was evaluated by their successful utilization in in-silico screening. Virtual screening using all three PREDICT models yielded enrichment factors 9-fold to 44-fold better than random screening. Namely, the PREDICT models can be used to identify active small-molecule ligands embedded in large compound libraries with an efficiency comparable to that obtained using crystal structures for non-GPCR targets.

Noninvasive imaging of reporter gene expression in living subjects

The development of noninvasive imaging technologies designed specifically for use with small animals has provided new paradigms for cancer research. Traditional molecular biology techniques are being melded with noninvasive imaging technologies to develop a new research domain, "molecular imaging." One of the most exciting advances in this research area is the adaptation and application of conventional reporter-gene imaging techniques, used extensively by cell and molecular biologists, to living animals. Using these new assays, investigators can image noninvasively, repeatedly, and quantitatively the location, magnitude, and duration of reporter-gene expression in living animals. This review will describe the instrumentation used for noninvasive imaging of reporter genes, the reporter genes developed for noninvasive imaging with radio-nuclide-based assays such as positron emission tomography, and the reporter genes used for optically based noninvasive assays using sensitive charged-coupled device cameras. Applications of noninvasive, whole-animal imaging to gene therapy for cancer, to cell-based therapy for cancer, to lymphocyte activation, to cancer progression and dissemination in engrafted models, to tumor initiation, promotion and metastasis in conditional murine models of cancer induction, and to the noninvasive monitoring of tumor responses to a variety of therapies are described. New developments in multimodality molecular imaging are discussed, and the potential utility of noninvasive reporter gene expression in the diagnosis and management of human cancer is presented.

Ligand-Binding Modes in Cationic Biogenic Amine Receptors

The binding site in G protein-coupled cationic biogenic amine receptors is formed in the cleft of the seven transmembrane segments. Upon binding the ligand, the receptors are activated or inactivated through the conformational changes of the transmembrane segments. G protein-coupled receptors bind four functionally distinct ligands; inverse agonists, antagonists, partial agonists, and full agonists. Hence, putative structural models for biogenic amine receptors corresponding to the ligand function (inverse agonist-, antagonist-, partial agonist-, and full agonist-bound receptor models) were built by using photointermediate models in the rhodopsin photocascade (M. Ishiguro et al. ChemBioChem. 2004, 5, 298-310). The ligand-receptor recognition of each was examined by modeling receptor-ligand complexes with functional ligands. The complex models suggested that each functional ligand binds the corresponding receptor structure and that ligand-specific interactions contribute to stabilization of the corresponding receptor structure.

Coupling of the human A1 adenosine receptor to different heterotrimeric G proteins: evidence for agonist-specific G protein activation

The present study investigates the effect of varying ligand structure on the ability of agonists to activate guanine nucleotide-binding proteins of the Gi, Gs and Gq families via the A(1) adenosine receptor. In CHO cells expressing this receptor, inhibition or potentiation of forskolin-stimulated cAMP accumulation was used as an end point to measure the activation of Gi and, in Pertussis toxin (PTX)-treated cells, Gs, respectively. Stimulation of inositol phosphate accumulation in PTX-treated cells was used as an index of Gq activation. CPA (N(6)-cyclopentyladenosine), NECA (5'-N-ethyl-carboxyamidoadenosine) and eight analogues of these ligands presented a range of guanine nucleotide-binding protein (G-protein)-activating profiles. Some ligands could only activate Gi (e.g. 2'deoxyCPA), some primarily Gi and Gs (and only weakly Gq) (e.g. 3'deoxyCPA), highlighting the importance of the ribose hydroxyls in agonist activation of multiple G proteins. CHA (N(6)-cyclohexyladenosine) activated Gi, Gs and Gq, but was more efficacious than CPA in activating Gs. The NECA analogues 5'-N-cyclopropyl-carboxamidoadenosine, 5'-N-cyclobutyl-carboxamidoadenosine and 5'-N-cyclopentyl-carboxamidoadenosine (CPeCA) also activated all three G proteins, although their ability to activate Gs and Gq (relative to CPA) was reduced with increasing substituent size, such that CPeCA produced only a small stimulation (at 100 microM) at Gq, but was a full agonist, relative to CPA, at Gi and Gs. This study suggests that the A(1) adenosine receptor can adopt agonist-specific conformations, arising from small changes in ligand structure, which lead to the differential activation of Gi, Gs and Gq.

Non-invasive imaging of reporter genes

Non-invasive, quantitative and repetitive imaging of biological processes in living animals is rapidly changing the way in which many experiments in models of human disease and normal physiological processes are conducted. This review summarizes the newest molecular imaging approaches to analyzing reporter gene expression, with particular emphasis on pre-clinical cancer research. Alternative modes of imaging are summarized, followed by descriptions of the major reporter gene systems now used for radionuclide imaging in vivo of gene expression. Several somatic delivery paradigms for co-ordinate expression of therapeutic and imaging genes are presented, and our own emphasis on the dopamine D2 receptor and Herpes Simplex Virus Type I thymidine kinase reporter genes are detailed.

Sequence analysis reveals how G protein-coupled receptors transduce the signal to the G protein

Sequence entropy-variability plots based on alignments of very large numbers of sequences-can indicate the location in proteins of the main active site and modulator sites. In the previous article in this issue, we applied this observation to a series of well-studied proteins and concluded that it was possible to detect most of the residues with a known functional role. Here, we apply the method to rhodopsin-like G protein-coupled receptors. Our conclusion is that G protein binding is the main evolutionary constraint on these receptors, and that other ligands, such as agonists, act as modulators. The activation of the receptors can be described as a simple, two-step process, and the residues involved in signal transduction can be identified.

Micro-PET imaging and small animal models of disease

Positron emission tomography (PET) has been used clinically to measure enzyme reactions, ligand-receptor interactions, cellular metabolism and cell proliferation. Until recently, however, PET has not been suitable for small animal models because of resolution limitations. Development of micro-PET instrumentation for small animal imaging and the availability of positron-emitting tracers has made this technology accessible for the non-invasive, quantitative and repetitive imaging of biological function in living animals. The development of new probes and positron-imaging based reporter genes has extended micro-PET applications to investigations of metabolism, enzyme activity, receptor-ligand interactions, protein-protein interactions, gene expression, adoptive cell therapy and somatic gene therapy. Because small animal PET is immediately extrapolatable to the clinic, laboratory advances should rapidly be translated to clinical practice.