|
|
|
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Department of Medicinal Chemistry and Molecular Pharmacology, School of Pharmacy and Pharmaceutical Sciences, Purdue University, West Lafayette, Indiana
Received June 22, 2007; accepted August 22, 2007
| Abstract |
|---|
|
|
|---|
As an example of this approach, the in silico-activated homology model of the h5-HT2A receptor developed in our laboratory (Chambers and Nichols, 2002
) has been used to predict several ligand-receptor interactions. Our results thus far have provided qualitative support for the model and have identified directions for further investigation of the structure-activity relationships of agonist ligands (Parrish et al., 2005
; Braden et al., 2006
; McLean et al., 2006a
,b
). Although such homology receptor models must be viewed with caution, they can often be validated by receptor mutagenesis experiments, along with complementary changes in ligand structure.
To refine our homology model and provide empirical support for the receptor functional topography it defines, we extended results derived from previous mutagenesis studies of 5-HT2A receptors. Agonist ligands for this receptor are of two chemical types: tryptamines and phenylalkylamines, the former being analogs of the endogenous ligand serotonin (5-hydroxytryptamine; 5-HT) and the latter comprising phenethylamines and phenylisopropylamines, all typically with polar substituents attached to the aromatic ring. A polar oxygen atom in the 4- or 5-position of tryptamines generally confers high affinity and activity at this receptor (McKenna and Towers, 1984
; Nichols and Glennon, 1984
). Tricyclic ergolines such as LSD can be considered special cases of rigidified tryptamines and are particularly interesting because they lack polar aromatic ring substituents. Extensive studies of the phenylalkylamine pharmacophore have revealed optimal polar ring substitution when methoxy groups are located at the 2- and 5-positions (Nichols, 1981
; Glennon et al., 1986
, 1989; Shulgin and Shulgin, 1991
; Nichols, 1994
, 1997
; Glennon, 1999
). Although most previous pharmacological investigations of 5-HT2A receptor agonists have been carried out in rats or in rat-derived preparations, we believe that the mechanistic effects and possible clinical utility of these compounds must be based on an understanding of the human receptor.
Before the advent of computationally derived homology models, putative binding sites had been identified using site-directed mutagenesis data. A strong consensus exists for the small molecule agonist binding site of these receptors to be nestled within the transmembrane domains (TM) of the seven heptahelical bundle, with contacts primarily between the ligand and residues in helices 3, 5, and 6 (van Rhee and Jacobsen, 1996
; Roth et al., 1998
).
The present study focused on serine residues at positions 5.43 [Ser5.43(239)] and 5.46]Ser5.46(242)] in TM5 of the 5-HT2A receptor [see Ballesteros and Weinstein (1995
) for numbering scheme]. These two residues are located one and two turns, respectively, above the highly conserved proline residue Pro5.50(246) in TM5, toward the extracellular face of the receptor. In other homologous GPCRs with a polar residue at position 5.43, mutagenesis to a nonpolar residue dramatically reduced the affinity and activity of nearly all agonists tested (Strader et al., 1989
; Ho et al., 1992
; Pollock et al., 1992
; Wess et al., 1992
) or increased the affinity of antagonists or was involved in species selectivity (Link et al., 1992
).
Mutation of a polar residue at position 5.46 in nonserotonin aminergic GPCRs also decreased affinity of agonists (Strader et al., 1989
; Wang et al., 1991
; Pollock et al., 1992
), affected the selectivity of agonists or antagonists (Mansour et al., 1992
), and/or decreased the affinity of antagonists (Gantz et al., 1992
; Leurs et al., 1994
). Functional activity was generally unaffected (van Rhee and Jacobsen, 1996
).
|
The data from the present study support the hypothesis that Ser5.43(239) serves as the most critical residue engaged by substituted tryptamine and phenylalkylamine agonists in the process of receptor binding and activation. We find that most oxygen-substituted tryptamines interact with both Ser5.43(239) and Ser5.46(242), whereas unsubstituted tryptamines probably adopt alternate binding orientations. The potent hallucinogen LSD seems to engage only Ser5.46(242) directly, whereas the phenylalkylamines engage only Ser5.43(239).
| Materials and Methods |
|---|
|
|
|---|
Cell Culture Methods. The construction of human embryonic kidney 293 human embryonic kidney cells with stable expression of wild-type h5-HT2A receptors at high (Hh2Ahi; 8000 fmol/mg) and moderate (Hh2Alo; 1600 fmol/mg) expression levels has been described previously (Parrish et al., 2005
; Braden et al., 2006
). All cell types were maintained in complete Dulbecco's modified Eagle's medium with 30 µg/ml Zeocin (Invitrogen) as described previously (Parrish et al., 2005
; Braden et al., 2006
). The Hh2Ahi cell line was used for competition binding assays and the Hh2Alo cell line was used for functional assays.
Establishing Mutant Human 5-HT2A Receptor Cell Lines. Vector construction and site-directed mutagenesis using the QuikChange kit (Stratagene, La Jolla, CA) method were performed as described previously (Braden et al., 2006
) with the following sense primers and corresponding antisense primers (Integrated DNA Technologies, Coralville, IA): S5.43(239)A, CTTTGTCCTGATCGGCTCTTTTGTGTCATTTTTCATTCCC; S5.46(242)A, CCTGATCGGCTCTTTTGTGTCATTTTTCATTCCCTTAACC. Mutant inserts verified by primer-directed sequencing (Retrogen, San Diego, CA) were then subcloned into the pBudCE4 vector (Invitrogen, Carlsbad, CA). Human embryonic kidney 293 cells were transfected, colonies selected, and receptor expression verified as described previously (Parrish et al., 2005
). The Hh2A/S5.43(239)A and Hh2A/S5.46(242)A cell lines were chosen for moderate expression (3800 and 3200 fmol/mg, respectively).
Radioligand Binding Assays. Membrane preparations, saturation isotherm, and competition binding assays were performed as described in detail previously (Chambers et al., 2002
; Marona-Lewicka et al., 2002
). Saturation isotherm binding assays used 0.25 to 10 nM [3H]ketanserin or 0.125 to 5.0 nM (±)-[125I]DOI. Ligands were tested in competition binding assays for their ability to displace 0.25 nM [125I]DOI. Binding assays at wild-type receptors used the high-expressing cell line (Hh2Ahi; 8000 fmol/mg); those at mutant receptors used the cell lines described above.
Inositol phosphate accumulation assays. Compounds were tested for their ability to stimulate hydrolysis of radiolabeled phosphatidyl inositides (PI) by measurement of radiolabeled inositol phosphate accumulation as described previously in detail (Marona-Lewicka et al., 2002
; Kurrasch-Orbaugh et al., 2003
). Each assay plate was normalized to wells stimulated with water (0%) and a concentration of serotonin chosen to be maximally stimulating (100%). PI accumulation assays at wild-type receptors used the lower expressing cell line (Hh2Alo; 1600 fmol/mg); assays at mutant receptors used the cell lines described above.
Computational Modeling/Virtual Docking. Ligand structures were virtually docked into an in silico-activated h5-HT2AR homology model constructed from the published X-ray crystal structure of bovine rhodopsin after undergoing rigid-body molecular dynamics subsequent to isomerization of the bound 11-cis-retinal to all-trans-retinal (Chambers and Nichols, 2002
). Local energy-minimized ensemble structures were obtained as described previously (Parrish et al., 2005
; Braden et al., 2006
). In brief, energy-minimized structures were virtually docked using the GOLD software package (Cambridge Crystallographic Data Center, Cambridge, UK). Docking algorithms for phenylalkylamines and ergolines were performed without constraints. Docking algorithms for tryptamines were performed with the following constraints: 1) a distance constraint of 1.5 to 3.5Å between the ligand amine (nonindole) nitrogen and the carbonyl carbon of the side chain of Asp3.32(155); 2) a protein hydrogen bond constraint for the side-chain oxygen of Ser5.46(242); and 3) a protein hydrogen bond constraint for the side-chain OH hydrogen of Ser5.43(239). The latter two constraints favorably weight the scoring of docking orientations that place any polar ligand atom in proximity to the protein atoms that could form a possible hydrogen bond without constraining a particular ligand atom. Ligand-receptor ensemble structures were each obtained by merging the highest ranked docking output ligand orientation structures with the input h5-HT2A homology model structure using the SYBYL software package (Tripos, St. Louis, MO), followed by energy minimization, molecular dynamics, and a final energy minimization simulation. Aggregates for molecular dynamics and minimization simulations were defined as residues more than 6 Å from the ligand as well as the backbone atoms. The molecular dynamics and minimization simulations were performed with constraints between Asp3.32(155) and the ligand aliphatic amine, and between any polar residue within 2.5 Å of a polar ligand group, including Ser3.36(159), Thr3.37(160), Ser5.43(239), or Ser5.46(242). Final energy minimization simulations also used these constraints.
Data Analysis. Prism software (GraphPad Software Inc., San Diego, CA) was used to calculate nonlinear regression curves based on the Cheng-Prusoff equation for a one-site model to obtain Ki (affinity) values for radioligand displacement and variable slope sigmoidal dose-response curves for EC50 (potency) and intrinsic activity from PI hydrolysis. This software package was also used to perform two-way ANOVA calculations with Bonferroni post-tests comparing wild-type and mutant pKi, pEC50, and intrinsic activity values. Values obtained from the mutant receptors were considered statistically distinguishable from wild-type if the models generated p < 0.05. Change in the standard Gibbs free energy (
G°) of binding due to the mutation was calculated from the Ki values at 25°C as follows:
(
G°) =
G°mutant –
G°wild-type = RT ln(Kmutant/kwild-type), where R is the gas constant, and T is absolute temperature. Changes in binding affinity and EC50 were transformed to normalize the scale by taking the difference of the log10 value (
pKi and
EC50, respectively) as follows:
pKi = pKi-mutant–pKi-WT = –logKi-mutant – (–logKi-WT) and
pEC50 = pEC50-mutant – pEC50-WT = –logEC50-mutant – (–logEC50-WT). Changes in intrinsic activity (
Int.Act.) were calculated as follows:
Int.Act. = Int.Act.mutant – Int.Act.WT.
Virtual docking figures were generated using pyMol (DeLano Scientific, San Carlos, CA; http://www.pymol.org). Amino acid residues are numbered with their position relative to the most highly conserved residue of that transmembrane region and their sequence position in the h5-HT2A receptor (Ballesteros and Weinstein, 1995
).
| Results |
|---|
|
|
|---|
|
Virtual Docking of Tryptamines, Ergolines, and Phenylalkylamines to an in Silico-Activated h5-HT2A Receptor Homology Model Indicated That Ser5.43(239) Interacted with Polar Substituents on the Aromatic Moiety of the Ligand, Whereas Ser5.46(242) Interacted with the Indole Nitrogen of Tryptamines and Ergolines, but Did Not Interact with Phenylalkylamines. Fig. 2 illustrates representative binding poses for virtual docking and subsequent energy minimization simulations of tryptamines, ergolines, and phenylalkylamines in an in silico-activated homology model of the h5-HT2A receptor. Virtual docking of LSD or the phenylalkylamines consistently oriented the ligand amine near Asp3.32(155) and the polar ligand substituents near polar residues in TM3 and TM5. Initial unconstrained docking orientations for tryptamines, however, were highly variable. Few of these orientations positioned the ligand so that the strong ionic interaction could occur between the protonated amine of the ligand and the highly conserved acidic residue at position 3.32, Asp3.32(155), believed to be the counter-ion in all biogenic amine-binding GPCRs (Ho et al., 1992
; Mansour et al., 1992
; Kristiansen et al., 2000
; Abekawa et al., 2003
). Applying a distance constraint between the side-chain nitrogen of the ligand and the carbonyl carbon of the side chain of Asp3.32(155) gave more reasonable and less variable highest-ranked output structures from docking simulations, although many structures were still oriented out of the putative small agonist binding pocket. Thus, additional docking constraints for tryptamines were introduced to encourage any hydrogen bonding between the ligand and either Ser5.43(239) and Ser5.46(242) but without forcing a particular receptor-ligand interaction.
|
The S5.43(239)A Mutation Detrimentally Affected the Binding and Activity of Tryptamines and Phenylalkylamines Predicted to Interact with This Residue. Table 2 and Fig. 3 present the results of competition binding assays with wild-type, S5.43(239)A mutant, and S5.46(242)A mutant h5-HT2A receptors. To aid visual interpretation, a loss in binding affinity is defined as a negative
(
G°). Binding of tryptamines with a polar aromatic substituent, such as 5-MeO-T, 5-HT, 5-MeO-DMT, psilocin, and 1-Me-5-HT, was attenuated by the S5.43(239)A mutation to a degree consistent with the loss of a hydrogen bond (0.5 to 1.5 kcal/mol; Fersht, 1988
), although 1-iPr-5-MeO-T was surprisingly unaffected. Both 5-HT and 5-MeO-DMT suffered the greatest losses of affinity, approximately 14- and 11-fold (
1.5 kcal/mol), respectively. The affinities of the phenylisopropylamines DOI, DOM, and DOH were detrimentally affected by the S5.43(239)A mutation, with 
G° well within the range given above for hydrogen bonds. The phenethylamines 2CI and 2CH also were detrimentally affected, but the 
G° values were less than 0.5 kcal/mol, and not statistically discernible from wild-type.
|
|
Table 3 and Fig. 4 present measures of functional activity of the test compounds using PI hydrolysis assays at wild-type and mutant receptors. In general, the trends observed in affinity changes for polar-substituted tryptamines were maintained in shifts of functional potency, as defined by
pEC50. All compounds were affected by this mutation to a certain degree, although in some cases the effect was very slight (e.g., 5-MeO-T). Again, 5-MeO-DMT and 5-HT were the polar-substituted tryptamines most affected, their potencies being reduced by approximately 20- and 35-fold, respectively. Psilocin was the only tryptamine to show a statistically discernible change in intrinsic activity (Table 3;
33% compared with the wild type) as well as a
30-fold decrease in potency.
|
|
By contrast, LSD lacks a polar substituent in the area indicated by virtual docking experiments to be near Ser5.43(239) and, not surprisingly, this mutation had no effect on its affinity and only a weak effect on its potency. Likewise, tryptamine and DMT also lack a polar substituent in this region; their binding was relatively unaffected and potency was only weakly affected. Although affinity of 5-Me-T also was unaffected, we were surprised to observe a dramatic 65-fold decrease in potency.
All phenylalkylamines tested showed significant decreases in potency (
pEC50) at the S5.43(239)A mutant receptor compared with wild type. DOM and DOH were the most markedly affected ligands, with decreases of
60- and
220-fold, respectively. Of all the phenylalkylamines tested at all the mutant receptors in this study, we observed a decrease in efficacy (
Int.Act.
–35%) only with 2CI at the S5.43(239)A receptor.
The S5.46(242)A Mutation in the h5-Ht2A Receptor Attenuated the Affinity and Activity of Some Tryptamines Predicted to Interact with This Residue, but Enhanced the Binding and Activity of Tryptamines with Alkyl Substitution on the Indole (N1) Nitrogen and Did Not Affect Phenylalkylamines. As Fig. 3 illustrates, the S5.46(242)A mutation attenuated the affinity only of LSD, 5-HT, and 5-MeO-DMT to a degree consistent with the loss of a hydrogen bond (Fersht, 1988
). Changes to the binding energetics of psilocin and 5-MeO-T approach the lower limit of this energy range (0.5–1.5 kcal/mol). The affinity of other tryptamines, particularly those lacking a polar ring substituent, was not significantly affected by the S5.46(242)A mutation. By contrast, the tryptamines with alkyl substitution on the indole (N1) nitrogen, 1-Me-5-HT and 1-iPr-5-MeO-T, showed 12- and 21-fold increased affinity, respectively.
Similar trends were generally observed in functional potency, as illustrated in Fig. 4. With the exception of 5-MeO-T, the potency of all the tryptamines with polar aromatic substituents, as well as LSD, was weakly to moderately decreased. No shifts in potency were observed in any of the ring-unsubstituted tryptamines 5-Me-T, tryptamine, or DMT. Again, we observed dramatic 45- and 51-fold increases in potency for the N(1)-alkyl analogs 1-Me-5-HT and 1-iPr-5-MeO-T, respectively.
The S5.46(242)A mutation did not affect the binding of any of the phenylalkylamines tested to a degree statistically discernible from wild type. At the S5.46(242)A mutant receptor, DOM, DOH, and 2CH suffered generally weak losses in potency (2-, 6-, and 5-fold, respectively), whereas the potency of DOI was slightly increased (
3-fold). This result is consistent with our observation that there is no significant difference between the affinity of phenethylamine agonists at the human and rat 5-HT2A receptors, which differ only in that the human receptor has a serine at this position, whereas the rat receptor has an alanine. The S5.46(242)A mutation did not lead to a significant change in efficacy for any of the phenylalkylamines tested.
A Complementary Alteration of Phenylalkylamine Structure Could Produce Binding and Activity Changes at the WT h5-HT2A Receptor Similar to Those of the Unmodified Ligand at the S5.43(239)A Mutation and Could Rescue Some Binding and Activity at the S5.43(239)A Mutant. As Table 4 indicates, removal of a polar group from the aromatic ring of the phenylisopropylamine DOM to yield 2-Et-DOM, 5-Et-DOM, or 5-H-DOM gave changes in affinities at the wild-type receptor consistent with the energetics of a lost hydrogen bond (Fersht, 1988
). These shifts are similar in magnitude to those observed when comparing the polar-substituted ligand homologues in wild-type versus S5.43(239)A mutant receptors [e.g., the effect of the S5.43(239)A mutation on DOI or DOM in Fig. 3]. Removal of the polar group at the 5-position of DOM to give 5-Et-DOM did not detrimentally affect binding at the S5.43(239)A mutant receptor. This mutant was still sensitive to the complete removal of the 5-methoxy to give 5-H-DOM (–
G° = 1.2). It was not, however, as sensitive as the wild-type receptor (–
G° = 2.0).
|
Table 4 also illustrates the difference in functional potency between wild-type and S5.43(239)A mutant receptors when a polar group is removed from phenylisopropylamines. Similar to the binding energies, the functional potencies at h5-HT2A wild-type receptors are sensitive to the removal of a polar group either at the 2- or 5-position of phenylalkylamines. By contrast, the S5.43(239)A mutant receptor is generally less sensitive or even insensitive to the loss of a polar group, but only at the ligand 5-position. This insensitivity is clearly evident in the comparison of DOM and 5-Et-DOM. Although we did not expect 2-Et-DOM and DOM to have significantly different potencies in the S5.43(239)A mutant receptor, the loss of sensitivity to substituting the 2-oxygen (2-Et-DOM;
pEC50 = 1.2 versus 0.7) is not as great as seen with either substitution (5-Et-DOM;
pEC50 = 1.4 versus 0.0), or removal (5-H-DOM
pEC50 = 2.3 versus 0.6) of the 5-oxygen substituent.
| Discussion |
|---|
|
|
|---|
A nondisruptive mutation to an alanine was chosen for each residue because it was anticipated that this change would abolish specific ligand-receptor interactions without affecting global receptor structure (Fersht et al., 1987
). Indeed, as Table 1 indicates, only a 2-fold shift in KD was observed for [3H]ketanserin at S5.43(239)A receptors (1.1 versus 2.2 nM) and a slight enhancement at S5.46(242)A receptors (1.1 versus 0.71 nM). A shift of
3-fold in the KD for [125I]DOI was observed at S5.43(239)A receptors (0.78 versus 2.19 nM) and approximately 2-fold at S5.46(242)A receptors (0.78 versus 1.7 nM). The mutant S5.43(239)A and S5.46(242)A receptors still were able to induce the release of radiolabeled inositol phosphates in a dose-dependent manner that was approximately comparable in magnitude with that of the wild-type receptor (data not shown), although functional potencies were slightly decreased for all chemical classes tested at the S5.43(239)A receptor. Although these findings indicate that the S5.43(239)A mutation may have slightly altered G-protein coupling, we do not believe it produces a significant global disruption of receptor structure. Thus, our initial hypothesis was that any differences in ligand affinity or activity between the wild-type and S5.43(239)A or S5.46(242)A mutant receptors would mainly result from a loss of direct interaction between the ligand and those serine residues.
The clearest results were obtained with the potent hallucinogen LSD, which was predicted not to interact with Ser5.43(239). The S5.43(239)A mutation had no effect at all on affinity or intrinsic activity and only a minor 3-fold effect on potency (Figs. 3 and 4). By contrast, the S5.46(242)A mutation caused a 4-fold shift in affinity (Table 2), resulting in a change of the standard free energy of binding (–
G°) of –0.8 kcal/mol (Fig. 4), within the –0.5 to –1.5 kcal/mol range determined for a hydrogen bond (Fersht, 1988
). This mutation similarly caused a shift of approximately 6-fold in potency for LSD, although it had no effect on intrinsic activity (Fig. 4). Our previous virtual docking studies (Chambers and Nichols 2002
) placed LSD into the agonist binding domain such that there was no interaction with Ser5.43(239), whereas the indole N(1)H engaged Ser5.46(242), located approximately one turn below Ser5.43(239). The present data are consistent with this predicted binding pose. The fact that LSD is such a rigid molecule indicates that the bound ligand-ensemble has little conformational freedom to involve Ser5.43(239) in binding. The slight loss of potency in the S5.43(239)A mutant does suggest, however, that even when a ligand does not directly engage this residue, it may still be important for secondary intramolecular interactions that create an "efficient" activated receptor state.
A second compound where there was no difference in affinity between the wild-type and mutant receptors is 5-Me-T, an example of a "mutated" ligand chosen to mimic the loss of interaction between Ser5.43(239) and the tryptamine 5-position oxygen. We initially envisioned that addition of the 5-methyl group to tryptamine might allow a favorable van der Waals interaction between the methyl group and A5.43(239). Indeed, although there is no difference in affinity of 5-Me-T in the wild-type and S5.43(239)A receptors (Fig. 3), this ligand does have higher affinity at the S5.43(239)A receptor compared with both 5-HT and tryptamine (Table 2).
We were surprised to find that affinity and potency at the S5.46(242)A mutant receptor indicate that 5-Me-T does not engage this residue (Figs. 3 and 4). This finding may indicate that this ligand does not bind in an orientation similar to other tryptamines. Further confounding interpretation is the fact that although this ligand has the same intrinsic activity at the wild-type and either of the two mutant receptors, there is a loss of potency in the S5.43(239)A mutant of more than 100-fold (Fig. 4). As was the case for LSD, we believe this result suggests that there may be additional indirect effects of the S5.43(239)A mutation that affect the ability of the ligand-receptor ensemble to adopt or maintain an "activated" state.
The use of "mutated" ligands and shifts in sensitivity indicated by Table 4 yields the most compelling evidence for the interaction of Ser5.43(239) with the phenylalkylamines, particularly the phenylisopropylamine subclass. At wild-type receptors, the data provide strong support for the necessity of the polar oxygen at the ligand 5-position, reflected in the sensitivity to its removal from DOM in the 5-ethyl or 5-hydrogen homologues 5-Et-DOM and 5-H-DOM, respectively. Moreover, the binding energetics and functional potencies at the S5.43(239)A mutant receptor indicate insensitivity to the absence of a 5-oxygen atom in 5-Et-DOM, the most direct 5-desoxy homolog of DOM. This mutant receptor is not completely insensitive to this change, as seen with 5-H-DOM, and also has reduced sensitivity to 2-Et-DOM. The magnitude of sensitivity loss is greatest, however, with the 5-position homologues. It should not be surprising that 5-H-DOM has lower affinity and potency than DOM in the S5.43(239)A mutant receptor, in that there are likely entropic and/or solvation factors with the removal of so many ligand atoms. A diagram of this interaction is provided in Fig. 5.
|
As further illustrated in Fig. 3, the classic phenylalkylamines affected by the S5.43(239)A mutation were limited to the phenylisopropylamine subclass (i.e., DOI, DOM, and DOH), whereas the S5.46(242)A mutation had no effect on the binding of all but one of the phenylalkylamines tested. This trend is maintained in the functional data, as illustrated in Fig. 4.
Compared with phenethylamines, the phenylisopropylamines DOI and DOM have additional steric bulk on the
carbon that can alter their conformation. We hypothesized that the S5.43(239)A mutation would eliminate the favorable interaction with the 5-methoxy of these ligands. Not surprisingly, we observed a 3-fold and nearly 8-fold loss of affinity at the mutant receptor for DOI and DOM, respectively, compared with the wild type (Table 2). This difference can be related to the nature of the nonpolar 4-substituent, which is known to affect both affinity and potency (Nichols and Glennon, 1984
). That is, DOI has nearly 10-fold higher affinity than DOM in the wild-type receptor. Thus, the 4-iodo substituent adds binding energy, probably through a Van der Waals interaction, so that the affinity of DOI is less affected by the S5.43(239)A mutation than DOM. A parallel to this reasoning in seen in the comparison of the potencies of these two ligands, as illustrated in Fig. 4. The potency of DOM in the S5.43(239)A mutant receptor is reduced more than 200-fold, whereas the potency of DOI is reduced only approximately 10-fold.
The phenethylamines 2CI and 2CH also were predicted to interact with Ser5.43(239), yet the changes in energetics of binding are not consistent with the range of energy for a hydrogen bond (Fig. 3) (Fersht, 1988
). One observes, however, that even though an interaction with the phenethylamines 2CI and 2CH is not apparent in the affinity assays, a moderate effect on the functional potency of these compounds is observed at the S5.43(239)A receptor, surpassing that of DOI. In addition, 2CI was the only phenylalkylamine tested that had significant and dramatically decreased intrinsic activity in the S5.43(239)A mutant (Table 3). In general, and as hypothesized, it seems that Ser5.43(239) is a key interaction site for most phenylalkylamines, whereas Ser5.46(242) does not interact with phenylalkylamine ligands in a manner that affects binding or potency at the h5-HT2A receptor.
The interaction of all the tryptamines with residue Ser5.46(242) is not so clearly evident. Indeed, it would seem that only oxygen-substituted tryptamines are detrimentally affected to a degree that would suggest loss of a hydrogen bond (Figs. 3 and 4). It seems possible that the strong polar interaction with Ser5.43(239) serves to orient the ligand, placing the indole N(1)H in proximity to Ser5.46(242), where a much weaker hydrogen bond is formed. In the absence of a polar ring substituent, the indole N(1)H interaction with Ser5.46(242) may be too weak, forcing those tryptamines to adopt a different binding orientation. Consistent with that argument, the ring-unsubstituted tryptamine and DMT, as well as the 5-methyl substituted 5-Me-T, were unaffected by the S5.46(242)A mutation. Moreover, the two polar substituted tryptamines with N(1)-alkylation, 1-Me-5-HT and 1-iPr-5-MeO-T, had dramatic increases in affinity and potency at the S5.46(242)A receptor compared with wild-type (Figs. 3 and 4). Additional support for an interaction between the tryptamine N(1)-substituent and either Ser5.46(242) or Ala5.46(242) in human or rat 5-HT2A receptors, respectively, was provided by the studies of Johnson et al. (1994
).
Another earlier investigation developed de novo models of the rat 5-HT2A receptor, where both Ser5.43(239) and Ala5.46(242) were inaccessible to binding and were instead projected toward TM4, or had Ser5.43(239) engaged with the tryptamine indole N(1) hydrogen (Shapiro et al., 2000
). Our data are not consistent with such an orientation of TM5 in the human 5-HT2A receptor, although we agree with the conclusion of that study that not all tryptamines may be binding in the same orientation. Instead, we believe our data support the hypothesis that Ser5.43(239) is directly engaged when binding tryptamines with polar substituents at either the 4- or 5-position. Although the ergoline LSD lacks a polar substituent to interact with Ser5.43(239), it has other structural features that may orient it within the binding site, so that it can engage Ser5.46(242). Even when a direct ligand-receptor interaction cannot occur, our data suggest that residue Ser5.43(239) may be involved in additional indirect interactions essential for potency at the 5-HT2A receptor. These indirect interactions are difficult to model and indicate that our model may not be fully accounting for all possible "active" receptor conformations. Moreover, it may be possible, as Strader et al. (1989
) suggested with the
-adrenergic receptor, that some ligands are able to bind to the mutant receptors in orientations that are competitive but not able to affect receptor activation.
Polar residues have been shown to be critical for agonist activity in a variety of monoamine GPCRs, as well as muscarinic acetylcholine receptors; representative TM5 sequences are illustrated in Fig. 6. In particular, in catecholamine receptors, the serine(s) at position 5.43 (and/or Ser5.42) probably engages the "meta" OH of the catechol moiety, a key structural feature of catecholamine agonists recognized for many years (e.g., Strader et al., 1989
).
|
Our data also indicate that not all tryptamines interact with the receptor in the same way, a conclusion reached by previous studies based on interactions with residue Ser5.43(239) in the rat 5-HT2A receptor (Johnson et al., 1997
; Shapiro et al., 2000
). The present results extend previous investigations into the role of these residues in the h5-HT2A receptor. As with investigations by our and other laboratories into the aromatic contacts of TM6 (Choudhary et al., 1993
, 1995
; Roth et al., 1997
; Braden et al., 2006
), it seems that the functional topography of the receptor binding site is such that a contact with TM5 is necessary for high affinity and functional potency of agonists. These results further support the topography and utility of our in silico-activated homology model of the h5-HT2A receptor, and we anticipate that the present findings can be extended to the understanding of other family A monoamine-binding GPCRs.
| Acknowledgements |
|---|
| Footnotes |
|---|
Article, publication date, and citation information can be found at http://molpharm.aspetjournals.org.
ABBREVIATIONS: 1-iPr-5-MeO-T, 1-isopropyl-5-methoxytryptamine; 1-Me-5-HT, 1-methyl-5-hydroxytryptamine (1-methylserotonin); 2CH, 2,5-dimethoxyphenethylamine; 2CI, 2,5-dimethoxy-4-iodophenethylamine; 2-Et-DOM, 2-ethyl-5-methoxy-4-methylphenyl-isopropylamine; 5-Et-DOM, 5-ethyl-2-methoxy-4-methylphenylisopropylamine; 5-H-DOM, 2-methoxy-4-methylphenylisopropylamine; 5-HT, 5-hydroxtryptamine (serotonin); 5-Me-T, 5-methyltryptamine; 5-MeO-DMT, 5-methoxy-N,N-dimethyltryptamine; 5-MeO-T, 5-methoxytryptamine; DMT, N,N-dimethyltryptamine; DOH, 2,5-dimethoxyphenyl-isopropylamine; DOI, 4-iodo-2,5-dimethoxyphenylisopropylamine; DOM, 2,5-dimethoxy-4-methylphenylisopropylamine; PI, phosphoinositide; TM, transmembrane; GPCR, G protein-coupled receptor; LSD, lysergic acid diethylamide; ANOVA, analysis of variance; WT, wild type.
Address correspondence to: Dr. David E. Nichols, Dept. of Medicinal Chemistry and Molecular Pharmacology, School of Pharmacy and Pharmaceutical Sciences, 575 Stadium Mall Drive, Purdue University, West Lafayette, IN 47907-2091. E-mail: drdave{at}pharmacy.purdue.edu
| References |
|---|
|
|
|---|
Baker D and Sali A (2001) Protein structure prediction and structural genomics. Science 294: 93–96.
Ballesteros JA and Weinstein H (1995) Integrated methods for the construction of three-dimensional models and computational probing of structure-function relations in G protein-coupled receptor. Methods Neurosci 25: 366–428.
Braden MR, Parrish JC, Naylor JC, and Nichols DE (2006) Molecular interaction of serotonin 5 HT2A receptor residues Phe339(6.51) and Phe340(6.52) with superpotent N-benzyl phenethylamine agonists. Mol Pharmacol 70: 1956–1964.
Chambers JJ, Kurrasch-Orbaugh DM, and Nichols DE (2002) Translocation of the 5-alkoxy substituent of 2,5-dialkoxyarylalkylamines to the 6-position: effects on 5-HT(2A/2C) receptor affinity. Bioorg Med Chem Lett 12: 1997–1999.[CrossRef][Medline]
Chambers JJ and Nichols DE (2002) A homology-based model of the human 5-HT2A receptor derived from an in silico activated G-protein coupled receptor. J Comput Aided Mol Des 16: 511–520.[CrossRef][Medline]
Choudhary MS, Craigo S, and Roth BL (1993) A single point mutation (Phe340–>LEU340) of a conserved phenylalanine abolishes 4-[125I]iodo-(2,5-dimethoxy)-phenylisopropylamine and [3H]mesulergine but not [3H]ketanserin binding to 5-hydroxytryptamine2 receptors. Mol Pharmacol 43: 755–761.[Abstract]
Choudhary MS, Sachs N, Uluer A, Glennon RA, Westkaemper RB, and Roth BL (1995) Differential ergoline and ergopeptine binding to 5-hydroxytryptamine2A receptors: ergolines require an aromatic residue at position 340 for high affinity binding. Mol Pharmacol 47: 450–457.[Abstract]
Crozier PS, Stevens MJ, and Woolf TB (2007) How a small change in retinal leads to G-protein activation: initial events suggested by molecular dynamics calculations. Proteins 66: 559–574.[CrossRef][Medline]
Fersht AR (1988) Relationships between apparent binding energies measured in site-directed mutagenesis experiments and energetics of binding and catalysis. Biochemistry 27: 1577–1580.[CrossRef][Medline]
Fersht AR, Leatherbarrow RJ, and Wells TN (1987) Structure-activity relationships in engineered proteins: analysis of use of binding energy by linear free energy relationships. Biochemistry 26: 6030–6038.[CrossRef][Medline]
Gantz I, Delvalle J, Wang LD, Tashiro T, Munzert G, Guo YJ, Konda Y, and Yamada T (1992) Molecular basis for the interaction of histamine with the histamine H2 receptor. J Biol Chem 267: 20840–20843.
Glennon RA (1989) Stimulus properties of hallucinogenic phenalkylamines and related designer drugs: formulation of structure-activity relationships. NIDA Res Monogr 94: 43–67.[Medline]
Glennon RA (1999) Arylalkylamine drugs of abuse: an overview of drug discrimination studies. Pharmacol Biochem Behav 64: 251–256.[CrossRef][Medline]
Glennon RA, Titeler M, and Young R (1986) Structure-activity relationships and mechanism of action of hallucinogenic agents based on drug discrimination and radioligand binding studies. Psychopharmacol Bull 22: 953–958.[Medline]
Ho BY, Karschin A, Branchek T, Davidson N, and Lester HA (1992) The role of conserved aspartate and serine residues in ligand binding and in function of the 5-HT1A receptor: a site-directed mutation study. FEBS Lett 312: 259–262.[CrossRef][Medline]
Johnson MP, Audia JE, Nissen JS, and Nelson DL (1993) N(1)-substituted ergolines and tryptamines show species differences for the agonist-labeled 5-HT2 receptor. Eur J Pharmacol 239: 111–118.[CrossRef][Medline]
Johnson MP, Loncharich RJ, Baez M, and Nelson DL (1994) Species variations in transmembrane region v of the 5-hydroxytryptamine type 2A receptor alter the structure-activity relationship of certain ergolines and tryptamines. Mol Pharmacol 45: 277–286.[Abstract]
Johnson MP, Wainscott DB, Lucaites VL, Baez M, and Nelson DL (1997) Mutations of transmembrane IV and V serines indicate that all tryptamines do not bind to the rat 5-HT2A receptor in the same manner. Brain Res Mol Brain Res 49: 1–6.[Medline]
Kao HT, Adham N, Olsen MA, Weinshank RL, Branchek TA, and Hartig PR (1992) Site-directed mutagenesis of a single residue changes the binding properties of the serotonin 5-HT2 receptor from a human to a rat pharmacology. FEBS Lett 307: 324–328.[CrossRef][Medline]
Kristiansen K, Kroeze WK, Willins DL, Gelber EI, Savage JE, Glennon RA, and Roth BL (2000) A highly conserved aspartic acid (Asp-155) anchors the terminal amine moiety of tryptamines and is involved in membrane targeting of the 5-HT(2A) serotonin receptor but does not participate in activation via a "salt-bridge disruption" mechanism. J Pharmacol Exp Ther 293: 735–746.
Kurrasch-Orbaugh DM, Watts VJ, Barker EL, and Nichols DE (2003) Serotonin 5-hydroxytryptamine 2A receptor-coupled phospholipase C and phospholipase A2 signaling pathways have different receptor reserves. J Pharmacol Exp Ther 304: 229–237.
Leurs R, Smit MJ, Tensen CP, Ter Laak AM, and Timmerman H (1994) Site-directed mutagenesis of the histamine H1-receptor reveals a selective interaction of asparagine207 with subclasses of H1-receptor agonists. Biochem Biophys Res Commun 201: 295–301.[CrossRef][Medline]
Li J, Edwards PC, Burghammer M, Villa C, and Schertler GF (2004) Structure of bovine rhodopsin in a trigonal crystal form. J Mol Biol 343: 1409–1438.[CrossRef][Medline]
Link R, Daunt D, Barsh G, Chruscinski A, and Kobilka B (1992) Cloning of two mouse genes encoding
2-adrenergic receptor subtypes and identification of a single amino acid in the mouse
2-C10 homolog responsible for an interspecies variation in antagonist binding. Mol Pharmacol 42: 16–27.[Abstract]
Mansour A, Meng F, Meador-Woodruff JH, Taylor LP, Civelli O, and Akil H (1992) Site-directed mutagenesis of the human dopamine D2 receptor. Eur J Pharmacol 227: 205–214.[CrossRef][Medline]
Marona-Lewicka D, Kurrasch-Orbaugh DM, Selken JR, Cumbay MG, Lisnicchia JG, and Nichols DE (2002) Re-evaluation of lisuride pharmacology: 5-hydroxytryptamine1A receptor-mediated behavioral effects overlap its other properties in rats. Psychopharmacology (Berl) 164: 93–107.[CrossRef][Medline]
McKenna DJ, and Towers GH (1984) Biochemistry and pharmacology of tryptamines and beta-carbolines. A minireview. J Psychoactive Drugs 16: 347–358.[Medline]
McLean TH, Chambers JJ, Parrish JC, Braden MR, Marona-Lewicka D, Kurrasch-Orbaugh D, and Nichols DE (2006a) C-(4,5,6-trimethoxyindan-1-yl)methanamine: a mescaline analogue designed using a homology model of the 5-HT2A receptor. J Med Chem 49: 4269–4274.[CrossRef][Medline]
McLean TH, Parrish JC, Braden MR, Marona-Lewicka D, Gallardo-Godoy A, and Nichols DE (2006b) 1-Aminomethylbenzocycloalkanes: conformationally restricted hallucinogenic phenethylamine analogues as functionally selective 5-HT2A receptor agonists. J Med Chem 49: 5794–5803.[CrossRef][Medline]
Nichols DE (1981) Structure-activity relationships of phenethylamine hallucinogens. J Pharm Sci 70: 839–849.[CrossRef][Medline]
Nichols DE (1994) Medicinal chemistry and structure-activity relationships, in Amphetamine and Its Analogs (Cho AK and Segal DS eds) pp 3–41, Academic Press, New York.
Nichols DE (1997) Role of serotoninergic neurons and 5-HT receptors in the action of hallucinogens, in Serotoninergic Neurons and 5-HT Receptors in the CNS (Baumgarten HG and Gothert M eds) pp 563–585, Springer-Verlag, Berlin Heidelberg
Nichols DE (2004) Hallucinogens. Pharmacol Ther 101: 131–181.[CrossRef][Medline]
Nichols DE and Glennon RA (1984) Medicinal chemistry and structure-activity relationships of hallucinogens, in Hallucinogens: Neurochemical, Behavioral, and Clinical Perspectives (Jacobs BL ed) pp 95–142, Raven Press, New York
Palczewski K, Kumasaka T, Hori T, Behnke CA, Motoshima H, Fox BA, Le T, I, Teller DC, Okada T, et al. (2000) Crystal structure of rhodopsin: a G protein-coupled receptor. Science 289: 739–745.
Parrish JC, Braden MR, Gundy E, and Nichols DE (2005) Differential phospholipase C activation by phenylalkylamine serotonin 5-HT 2A receptor agonists. J Neurochem 95: 1575–1584.[CrossRef][Medline]
Pollock NJ, Manelli AM, Hutchins CW, Steffey ME, MacKenzie RG, and Frail DE (1992) Serine mutations in transmembrane V of the dopamine D1 receptor affect ligand interactions and receptor activation. J Biol Chem 267: 17780–17786.
Roth BL, Shoham M, Choudhary MS, and Khan N (1997) Identification of conserved aromatic residues essential for agonist binding and second messenger production at 5-hydroxytryptamine2a receptors. Mol Pharmacol 52: 259–266.
Roth BL, Willins DL, Kristiansen K, and Kroeze WK (1998) 5-hydroxytryptamine2-family receptors (5-hydroxytryptamine2a, 5-hydroxytryptamine2b, 5-hydroxytryptamine2c): where structure meets function. Pharmacol Ther 79: 231–257.[CrossRef][Medline]
Shapiro DA, Kristiansen K, Kroeze WK, and Roth BL (2000) Differential modes of agonist binding to 5-hydroxytryptamine(2A) serotonin receptors revealed by mutation and molecular modeling of conserved residues in transmembrane region 5. Mol Pharmacol 58: 877–886.
Shulgin A and Shulgin A (1991) PIHKAL: A Chemical Love Story. Transform Press, Berkeley, CA.
Strader CD, Candelore MR, Hill WS, Sigal IS, and Dixon RA (1989) Identification of two serine residues involved in agonist activation of the beta-adrenergic receptor. J Biol Chem 264: 13572–13578.
van Rhee AM and Jacobsen KA (1996) Molecular architecture of G protein-coupled receptors. Drug Dev Res 37: 1–38.[CrossRef]
Wang CD, Buck MA, and Fraser CM (1991) Site-directed mutagenesis of
2A-adrenergic receptors: identification of amino acids involved in ligand binding and receptor activation by agonists. Mol Pharmacol 40: 168–179.[Abstract]
Wess J, Maggio R, Palmer JR, and Vogel Z (1992) Role of conserved threonine and tyrosine residues in acetylcholine binding and muscarinic receptor activation. A study with M3 muscarinic receptor point mutants. J Biol Chem 267: 19313–19319.
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||