|
Advertisement | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 280, Issue 41, 34813-34822, October 14, 2005
A Conserved Glutamate Residue in Transmembrane Helix 10 Influences Substrate Specificity of Rabbit OCT2 (SLC22A2)*
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
In human kidney, OCT2 appears to be the predominant OC transporter. Expression of mRNA for OCT2 far exceeds that for OCT1, and immunocytochemistry clearly shows a basolateral expression for OCT2 and little or no presence of OCT1 (9). There are, however, clear species differences in this profile, with substantial expression of both OCT1 and OCT2 evident in both rat and rabbit kidney (10, 11). Notably, the relative functional expression of OCT1 and OCT2 has only been established in the mouse (8) and rabbit (12), with clear evidence indicating that both OCT1 and OCT2 can play a quantitatively significant role in OC secretion. Nevertheless, OCT2, which interacts with a variety of OCs, including many clinically used drugs (e.g. procainamide and cimetidine), hormones (e.g. norepinephrine and dopamine), and toxic substances (e.g. 1-methyl-4-phenylpyridinium (MPP+)), evidently plays a key role in all species in governing the entry of OCs from the blood into the renal tubule, thereby controlling the first step in renal secretion of OCs.
The orthologs of OCT1 and OCT2 have been cloned from human (13, 14), rat (15, 16), mouse (17, 18), and rabbit (19, 20), and in each species, these homologs share
70% amino acid identity and >90% sequence similarity. This high degree of shared structure is consistent with the shared functional characteristics of these homologous transporters. Indeed, there is strong evidence supporting the contention that these transporters arose by gene duplication (21). However, despite the similar affinities for some substrates generally observed between OCT1 and OCT2 homologs in any of these species (e.g. for TEA) (11), these transporters are typically distinguished from one another by markedly different affinities for selected substrates (11, 22). These differences in selectivity presumably reflect differences in the identity of a subset of the amino acid residues that distinguish OCT1 from OCT2 homologs. Here, we have tested the hypothesis that shifts in selectivity of OCT2 toward a more OCT1-like phenotype will result from mutation of one or more of the residues that distinguish these two homologous proteins. We identified one such residue in OCT2 (Glu447 in the rabbit ortholog) that exerted a marked influence on substrate selectivity. This observation has been interpreted in light of a proposed homology model of the structure of rabbit OCT2 and a recently described model of the structure of rat OCT1 (23), both of which are based upon recent developments in the understanding of the structure of proteins in the major facilitator superfamily (MFS) of transporters.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
Cell Culture and TransfectionChinese hamster ovary (CHO) cells were grown at 37 °C in a humidified atmosphere (5% CO2) in plastic culture flasks. The medium was Kaighn's modified Ham's F-12 medium supplemented with 10% fetal calf serum. Cultures were split every 3 days. Cells were transfected via electroporation (ECM 630, BTX, San Diego, CA) with a plasmid containing the cDNA sequence for the desired construct. Briefly, cells were transfected with 10 µg of DNA at 260 V (time constant of
25 ms) and seeded onto 12-well plates at 320,000 cells/well. Uptake was typically measured 48 h after transfection, after the cells had reached confluence.
Measurement of TransportCHO cells were incubated at room temperature (
25 °C) in Waymouth buffer (135 mM NaCl, 13 mM HEPESNaOH, pH 7.4, 28 mM D-glucose, 5 mM KCl, 1.2 mM MgCl2, 2.5 mM CaCl2, and 0.8 mM MgSO4) to which labeled substrate and appropriate test agents were added. Uptake was stopped by rinsing the cells three times with 2 ml of ice-cold Waymouth buffer. The cells were then solubilized in 400 µl of 0.5 N NaOH with 1% (v/v) SDS. (The extract was neutralized with 200 µl of 1 N HCl.) Accumulated radioactivity was determined by liquid scintillation spectrometry. Rates of uptake are expressed as moles/cm2 of nominal cell surface of the confluent monolayer.
ImmunocytochemistryCHO cells were electroporated with plasmid DNA containing the various epitope-tagged sequences and seeded onto coverslips. Immunocytochemistry was generally performed on a confluent monolayer 24 h after plating. Cells were fixed in ice-cold 100% methanol for 20 min, washed with phosphate-buffered saline (PBS; 137 mM NaCl, 2.7 mM KCl, 8.0 mM Na2HPO4, and 1.5 mM KH2PO4,pH8;all washes were done in triplicate), and incubated for 1 h with anti-V5 antibody (Invitrogen) diluted 1:500 in PBS. The cells were then washed with PBS and incubated for 1 h in the dark with fluorescein isothiocyanate-conjugated goat anti-mouse antibody (Molecular Probes, Inc.) diluted 1:1000 in PBS.
Western BlottingCrude CHO membranes were prepared as described (24), and protein concentration was determined by the Bradford protein assay (Bio-Rad). Thirty micrograms of total membrane protein was separated on a 10% Tris/glycine/SDS-polyacrylamide gel and electrophoretically transferred to a nitrocellulose membrane (25). The membrane was blocked for 1 h with 6% milk in PBS at room temperature, followed by incubation overnight with anti-V5 antibody (1:5000 dilution) at 4 °C. After extensive washing with 0.05% PBS/Tween, the membrane was incubated with horseradish peroxidase-conjugated goat anti-mouse IgG (Promega Corp.). The secondary antibody was detected with the SuperSignal West Femto maximum sensitivity substrate kit (Pierce).
Homology Modeling of the Three-dimensional Structure of rbOCT2 The rbOCT2 sequence was used as a "probe" to search homologous sequences (PSI-BLAST, NCBI Database) and sequence-based structural relatives (3D-PSSM) (26). GlpT (Protein Data Bank code 1PW4) and LacY (code 1PV6) were returned as the top hits in both search algorithms. Notably, although the sequence identity (between rbOCT2 and GlpT or LacY) was low (<15%), the sequence similarity was
50%. The transmembrane
-helices of rbOCT2 were predicted using TMHMM (27). The putative transmembrane-spanning helices (TMHs) of rbOCT2 were then aligned with the respective
-helices of LacY and GlpT using ClustalW. In parallel, we also computed the hydropathy plots and helical wheel plots. These suggested that TMH1 and TMH12 were not aligned correctly with respect to the crystal structures. Significantly, the highly conserved motifs that fall between TMH2 and TMH3 and between TMH7 and TMH8 of MFS proteins (28) were well aligned between the rbOCT2 and GlpT sequences. Based on these factors, we performed some deletions and modifications to fine-tune the sequence alignment. The final sequence alignment was used for target-template homology modeling, which was performed using ICM 3.0 (MolSoft L.L.C., La Jolla, CA). The model was refined in SYBYL 6.9 (Tripos, Inc., St. Louis, MO) using the Powell method (29), and the minimization routine was run for 7000 cycles. The model was validated using PRO-CHECK (30), WHATIF (31), and PROVE (32).
StatisticsAll data are expressed as means ± S.E., with calculations of standard errors based on the number of separate experiments. Statistical comparisons were performed using either a t test or one-way analysis of variance, followed by a Student-Newman-Keuls post hoc test (ProStat 3.81, Poly Software International, Pearl River, NY).
| RESULTS |
|---|
|
|
|---|
Effect of Mutations on the Kinetics of Substrate InteractionWe elected to use the rabbit ortholog of OCT2 to assess the influence of the selected mutations on transport activity because rbOCT2 generally displays a greater degree of homolog-specific selectivity over rbOCT1 than noted for the human orthologs of these transporters (35). In particular, the >20-fold difference in affinity for cimetidine that distinguishes rbOCT2 from rbOCT1 (11, 20) was anticipated as being a valuable tool for determining the bases of such differences in selectivity. We also reasoned that structural information gained from the study of rbOCT orthologs is likely to be relevant to understanding the structure of OCTs in other species as well. Fig. 1 compares the uptake of [3H]TEA by wild-type rbOCT2 with that supported by the N353L, R403I, and E447Q single, double, and triple mutants of rbOCT2. Each mutant displayed significant transport of TEA. Uptake into transiently transfected CHO cells was nearly linear for at least 5 min (data not shown), and 5-min uptakes were generally used to generate initial rates of transport for kinetic analyses.
OCT1 and OCT2 have a similar apparent affinity for TEA (Kt = 141 and 165 µM, respectively; p > 0.05) (TABLE ONE; see also Ref. 22), and we hypothesized that replacing OCT2-specific residues with those found in OCT1 would have comparatively little effect on the affinity of the mutant transporter for TEA. As shown in Fig. 2 and TABLE ONE, this proved to be the case: although there was a general trend for a slight increase in apparent affinity of some of the mutant transporters for TEA (i.e. a decrease in IC50), it was not significant. We did expect, however, that one or more of the residue replacements would result in a decrease in affinity for cimetidine and for the fluorescent cation N,N,N-trimethyl-2-(methyl(7-nitro-2,1,3-benzoxadiazol-4-yl)amino)ethanaminium iodide (NBD-MTMA) (36), i.e. a shift toward a more "OCT1-like" selectivity. This proved to be the case: although neither the N353L nor R403I mutation by itself resulted in a significant change in the IC50 for inhibition of TEA uptake by cimetidine and NBD-MTMA, every protein that included the E447Q mutation displayed a significant decrease in apparent affinity for cimetidine and NBD-MTMA (Fig. 2). The IC50 for cimetidine inhibition of TEA uptake by the E447Q mutant was increased >15-fold (from 2.1 to 33.4 µM; p < 0.05), and that for NBDMTMA was increased 2.5-fold (from 11.6 to 28.3 µM; p < 0.05). The inclusion of both the E447Q and N353I conversions produced a further decrease in apparent affinity for both cimetidine (IC50 increase from 33.4 to 104 µM) and NBD-MTMA (from 28.0 to 78.9 µM) (p < 0.05). Interestingly, the affinity of the N353I/R403I/E447Q triple mutant for both cimetidine and NBD-MTMA did not differ significantly from that displayed by OCT1.
|
|
|
|
|
|
|
3.5 Å. Furthermore, the LacY structure was solved with (Protein Data Bank code 1PV7) and without (code 1PV6) a ligand in the substrate-binding site. There is now strong evidence to support the contention that all members of the MFS share a common structural fold with similar topological organization of their
-helices (40). With that in mind, we developed a homology model of the three-dimensional structure of rbOCT2 using the GlpT structure as the template. Because of the comparatively low degree of similarity between the full-length sequences of GlpT and rbOCT2 (13.8% amino acid identity and a total of 41.5% similarity), aspects of the alignment were, of necessity, based on structural insight (see supplemental Fig. A). An 11-residue motif between TMH2 and TMH3 and between TMH8 and TMH9 that is highly conserved in MFS members (28) was evident in the alignment of GlpT and rbOCT2, and this enhances the confidence in the alignment of the adjacent helices (i.e. helices 2/3 and 7/9). However, there is less confidence associated with the alignment adjacent to the "long loops," i.e. the long cytoplasmic loops found in both GlpT and the OCTs between TMH6 and TMH7 and the long external loop found between TMH1 and TMH2 in the OCTs (which is not found in other MFS transporters). To facilitate modeling, these loops were edited out of the sequences, but the rationale for selection of precise regions for elimination introduced a degree of arbitrariness to the resulting alignment of the adjacent helices. The use of TMH prediction algorithms provided a degree of independent assessment of the general accuracy of the model. It should be noted that the computationally derived position of TMH1 of both OCT1 and OCT2 routinely placed the valine at position 40 within helix 1. However, the homologous position in the closely related protein OAT1 (Asn39) has been shown to be a site of glycosylation (41). Consequently, we adjusted the alignment to ensure that Val40 was outside of the membrane (supplemental Fig. A).
Fig. 6 shows the resulting three-dimensional homology model of rbOCT2 based on the GlpT template. The root mean square deviation for the superposition of the refined and minimized rbOCT2 model based on the C-
backbone with the crystal structure of GlpT was 2.3 Å. Clearly evident is the pseudo 2-fold symmetry that arises from the helical organization of the N- and C-terminal halves of the protein and that appears to be a common structural feature of MFS proteins (40). The water-filled "pocket" that exists between these two domains is suspected to contain the sites involved in substrate binding (as found for both GlpT (38) and LacY (37)). This pocket is formed by the relative positions of TMH1, TMH2, TMH4, and TMH5 (in the N-terminal half of the protein) and TMH7, TMH8, TMH10, and TMH11 (in the C-terminal half of the protein). Consequently, it is reasonable to hypothesize that selected residues found in these helices may reside in or adjacent to the water-filled cavity that composes the pocket and could strongly influence the binding of substrates to rbOCT2. Lending support to this hypothesis was the observation that Glu447 proved to reside within TMH10 and was directed toward the putative binding pocket. Additional support for the structure presented in Fig. 6 was derived from the observation that residues shown in other studies to influence binding were also found to be directed toward the putative binding region of OCT2, including residues in TMH4 (Trp217, Tyr221, and Thr225) (23), TMH11 (Asp474) (42), and two additional residues in TMH10 (Ile442 and Leu446) (34). Interestingly, Asp474, the first residue implicated in the binding of substrates to OCTs (i.e. Asp475 in rat (r) OCT1) (42), was found at a position within the putative binding pocket that was virtually identical to the positions of the residues associated with the binding of sugar to LacY (37) and the postulated binding of glycerol phosphate to GlpT (38).
Docking AnalysisOne of the principal reasons for developing a three-dimensional structure of OCT2 is development of a predictive model of transporter/substrate interactions. Consequently, it was of particular interest to compare the predicted docking energies of substrates within the putative binding region of the rbOCT2 homology model with experimentally determined kinetic parameters for substrate/inhibitor interactions with rbOCT2. We generated three distinct docking surfaces, each derived from appropriate "receptor description files," for use with the FlexX docking program (BioSolvIT GmbH). The first two were based upon the location of ligand within the LacY binding site. The rbOCT2 model was superimposed over the LacY/ligand crystal structure (Protein Data Bank code 1PV7), and the rbOCT2 residues that were within either 5 or 10 Å (first and second models, respectively) of the ligand (
-D-galactopyranosyl 1-thio-
-D-galactopyranoside) were identified and used to prepare receptor description files. The third model was based upon a docking surface that contained only those residues identified in this and previous (23, 34, 42) studies as influencing substrate binding to OCTs (i.e. the residues shown in Fig. 6, A and B). We then used FlexX to generate G scores for the binding interaction of each of 14 different compounds selected from a set of compounds examined in a previous study on the molecular determinants of substrate interaction with OCT2 (35). The selected compounds (TABLE TWO) provided a range of three-dimensional structures and apparent affinities for OCT2. FlexX generates four parameters related to ligand docking: PMF score, G score, D score, and F score. Scatter diagrams of these scores versus the experimentally determined IC50 values for the test compounds showed that the measured IC50 values most closely correlated with the G score (values listed in TABLE TWO). The G score is related to the binding free energy (Gb) by the empirically determined relationship
Gb = (G score)/4 (reflecting the uniform weighting of the four Flex parameters in our analysis) (43, 44). Apparent dissociation constants (KD) for ligand transport interactions were then derived from the relationship KD = exp(
Gb/RT).
|
-D-galactopyranosyl 1-thio-
-D-galactopyranoside ligand. To determine whether the inclusion of these residues in the 5 Å docking surface of rbOCT2 would improve its calculated interaction with the test battery of OC ligands, we constructed a new docking surface consisting of the "5 Å residues" plus the additional "specific" residues noted above. Although the predicted KD values obtained using this docking surface (TABLE TWO) correlated with the measured IC50 values for the test OCs, the correlation (r2 = 0.544) was not as good as that obtained using the 5 Å surface (data not shown). Fig. 6C shows a structural model of the 5 Å docking surface and the predicted lowest energy configuration of TEA within the associated cleft.
|
-D-galactopyranosyl 1-thio-
-D-galactopyranoside ligand within the putative binding pocket. As shown in Fig. 7D, there was no correlation between the measured IC50 values and KD values predicted for docking to either GlpT or LacY. | DISCUSSION |
|---|
|
|
|---|
We hypothesized that the homolog-specific selectivity of OCT2, as exemplified by its high affinity for cimetidine and NBD-MTMA compared with rbOCT1 (11), is a consequence of the influence of a comparatively small number of amino acid residues within the sequences of these two homologous proteins that are conserved across all orthologs of OCT1 as one amino acid, but are conserved as a different amino acid in all OCT2 orthologs. There are 24 such residues in OCT1/OCT2 (supplemental Table A). To test this hypothesis, we focused on the effect on rbOCT2 selectivity of three residues, i.e. N353L, R403I, and E447Q. These sites were targeted because the replacement residue for each (i.e. the residue expressed in OCT1) has markedly different physicochemical properties compared with that in OCT2. In addition, preliminary studies on the activity of OCT1/OCT2 chimeric constructs suggested that the structural factors that influence homolog-specific selectivity reside predominantly in the C-terminal half of these transporters.3 Consistent with our hypothesis, none of the mutants (single, double, or triple) displayed a significant change in apparent affinity for TEA, a substrate for which OCT1 and OCT2 have a similar affinity (Fig. 2). Selected residue replacements did, however, exert a significant effect on the apparent affinity of rbOCT2 for cimetidine and NBD-MTMA. All the replacement combinations that included E447Q had a decreased affinity (i.e. increased IC50) for these selective substrates (TABLE ONE). Interestingly, although neither N353L nor R403I by itself had an effect on the interaction with either compound, the inclusion of N353L with E447Q exerted a synergistic effect on the decrease in affinity for cimetidine and NBD-MTMA (Fig. 2). These data suggest that Glu447 and Asn353 play significant roles in defining a binding surface within rbOCT2 that is responsible for its characteristic homolog-specific selectivity.
Replacing Glu447 with residues other than glutamine also caused altered binding profiles in the resulting rbOCT2 mutants. Replacing the anionic glutamate residue with either lysine or arginine completely eliminated transport activity; the resulting mutants did not support uptake of TEA, cimetidine, MPP+, or PAH despite the fact that the mutant proteins appeared to be expressed in the plasma membrane of transfected CHO cells (Fig. 3). The E447A mutant displayed greatly reduced uptake of both TEA and MPP+, although it showed accumulation of the fluorescent cation NBD-MTMA. The E447L mutant supported nearly normal uptake of MPP+ (and of NBD-MTMA), but only very low rates of TEA and cimetidine uptake. Notably, although the low rate of TEA transport was associated with a marked decrease in apparent affinity of the E447L mutant for TEA, cimetidine continued to be a high affinity inhibitor of E447L transport activity (of MPP+) (Fig. 5) despite the fact that cimetidine was itself transported only poorly (Fig. 4).
The evidence presented here implicating Glu447 as a key contributor to the binding properties of rbOCT2 is consistent with the recent observation by Gorboulev et al. (34) that the homologous residue in rOCT1, i.e. Gln448, plays a key role in the low affinity of rOCT1 for corticosterone compared with rOCT2. Two additional residues in the putative TMH10 in rOCT1, i.e. Ala443 and Leu447, were also shown to influence the affinity of that transporter for corticosterone, with replacement of these residues with those found in OCT2 resulting in a marked increase in affinity of the OCT1 mutant for corticosterone. In addition, mutation of Leu447 and Gln448 was associated with differential changes in affinity for rOCT1 substrates, e.g. although the L447Y/Q448E double mutant showed no change in affinity for TEA (consistent with our observations (Fig. 2)), its affinity for MPP+was substantially increased. Significantly, Gorboulev et al. (34) reported that corticosterone binding to one of their mutants exerted an allosteric effect on the binding of other cationic substrates, leading to the conclusion that the binding surface of OCTs is large enough to support the simultaneous binding of multiple ligands, which can then exert short-range allosteric interactions with the binding region.
The differential influence of selected residues on substrate interactions has been observed in several additional studies. In their assessment of the influence on transport of amino acid residues in TMH4 of rOCT1, Popp et al. (23) reported that replacement of Trp218 with tyrosine and Tyr222 with leucine resulted in significant decreases in Kt values for both TEA and MPP+, whereas the Y222F and T226A mutants decreased the Kt values for only TEA or MPP+, respectively. In addition, in their earlier study on the influence of acidic residues on the transport activity of rOCT1, Gorboulev et al. (42) found that the D475E mutant displayed a marked increase in apparent affinity for TEA, with no change in interaction with MPP+. Lending further support to the view that Asp475 plays a key role in defining a binding surface in OCTs is the observation by Wolff et al. (45) that the homologous position in the sequence of all members of the OAT family (which are also members of the OCT transport family (21)) is filled with a cationic residue. Although the R478D mutant of the flounder ortholog of OAT1 still supported transport of the monovalent anion PAH (albeit at a reduced rate compared with the wild-type protein), this mutation eliminated interaction with the dicarboxylate glutarate. In the rat ortholog of OAT1, similar differential effects on interaction with different substrates have also been reported in studies that used site-directed methods to probe the influence of selected residues on the transport activity of rOAT3. The homologous replacement in rOAT3 (i.e. R454D) showed a profound change in interaction with PAH (33), and mutations of several hydrophobic residues in TMH7 (W334A, F335A, Y341A, and Y342Q) and one in TMH8 (F362S) resulted in a marked decrease in transport of the hydrophilic substrates PAH and cimetidine, with comparatively little effect on transport of the more hydrophobic substrate estrone sulfate (46), leading to the conclusion that the structure of OAT3 includes a general binding domain with no single binding site.
The homology model of the three-dimensional structure of rbOCT2 provides a means to interpret, in the context of the potential mechanism of substrate binding and translocation, the effects of specific residues on transport activity as determined in site-directed studies. The development of the rbOCT2 model (Fig. 6) was made possible by the recent observation that the high resolution crystal structures of two other MFS transporters, LacY (37) and GlpT (38), revealed a common organization of TMHs. Despite the low degree of sequence identity of these two transporters (
13%), the observed structural homology supported the contention that the structure of other MFS transporters may be deduced using the methods of homology modeling (40). This approach has been subsequently used to develop structural models of human GLUT1 (47, 48) and rOCT1 (23) based upon the template structures of GlpT and LacY, respectively. Not surprisingly, the GlpT-based model of the rbOCT2 structure reported here displays a number of key similarities to the LacY-based model of the rOCT1 structure developed by Popp et al. (23). Most significantly, all the residues that were identified in rOCT1 as exerting influence on ligand binding (on both TMH4 and TMH10) and that were observed to be directed toward the putative pore (or "cleft") region of the rOCT1 model structure were also observed to face the pore region of the rbOCT2 model described here (Fig. 6). There are, however, differences in the proposed models as well. First, the rbOCT2 and rOCT1 models were based upon different template structures (i.e. GlpT and LacY, respectively). Although these structures share a remarkably similar fold, they nevertheless show distinct differences in the absolute positions of residues that are otherwise homologous in their positions within the two sequences (the C-
root mean square deviation for superimposition of the two sequences is on the order of 3.5 Å), and these differences in position of homologous residues clearly carry over to the homology models of rbOCT2 and rOCT1. Second, the alignment of rbOCT2 with GlpT used to generate the OCT2 structural model differs slightly from the alignment of rOCT1 with LacY used by Popp et al. (23). Interestingly, these differences are particularly evident in TMH4 and TMH10, where the alignments are shifted by four residues or approximately a full helical turn. Consequently, both of the resulting homology models direct the same set of residues toward the putative binding cleft, but these residues reside at different levels relative to the "center point" of the membrane. The differences in alignments reflect different decisions concerning the length of selected loops found between the TMHs. Given the ambiguities associated with constructing alignments of sequences that share comparatively little amino acid identity, it is premature to conclude whether either alignment is correct.
The similarity of the protein fold of LacY and GlpT and the identification in both structures of similar ligand-binding regions (i.e. within the cleft between the N- and C-terminal halves of the proteins) suggested that the homologous region in rbOCT2 might be involved in binding interactions associated with transporter/substrate interaction. This hypothesis was supported by the strong correlation between the binding energies calculated for a battery of test compounds and their measured IC50 values for inhibition of OCT2 transport activity (Fig. 7). Notably, no correlation was noted between the measured IC50 values and the binding energies calculated from the structures of GlpT and LacY, suggesting that the rank order of substrate/inhibitor docking calculated to occur within the modeled OCT2 reflects the specific amino acid composition of the docking surface within the putative binding region of the protein.
The relationship between predicted and measured binding interactions of substrate with the homology model of the rbOCT2 structure must be interpreted cautiously. Despite the evident correlation between these parameters for the test set of substrates (Fig. 7), not all predictions based upon the postulated docking surface of the OCT2 homology model proved to be as accurate. For example, although the calculated KD for cimetidine of the 5 Å docking surface that included the E447Q mutation (following energy minimization of the E447Q-containing mutant protein) was increased slightly (i.e. from the wild-type value of 10 µM to the "mutant" value of 34 µM), the KD for TEA was increased much more (from 9 µM to >2 mM). Similarly, the KD profile of the E447L mutant showed some consistencies with the observed behavior of the mutant protein, i.e. maintenance of reasonably high affinity for MPP+and a marked decrease in affinity for TEA; the profile also included a profound decrease in affinity for cimetidine, which contrasted sharply with the observed profile of behavior (Fig. 5). These deviations from predicted behaviors serve to underscore the current limitations of homology modeling with respect to providing a quantitative predictive model of substrate interactions. Even comparatively slight differences in the spatial distribution of amino acid residues within a docking region can, as we have seen, result in marked differences in predicted binding interactions. As noted earlier, the alignment of target and template sequences of proteins that share a low degree of homology offers many opportunities for small shifts in the position of residues within the homology model. In addition, the absolute spatial position of OCT2 residues within the hypothetical docking region reflects the use of the GlpT template structure; those positions would have been different had LacY been used as the template (and indeed, they are different in the homology model of rOCT1 (23) that did use LacY). Existing docking algorithms make simplifying assumptions to facilitate their calculations; for example, the FlexX algorithm assumes that the docking surface is "rigid," and that may not be the case, particularly if the docking region of OCTs is viewed as a comparatively large surface with multiple binding domains. Finally, it should be emphasized that all homology models that arise from the deduced structures of MFS transporters reflect an overall orientation in which the binding surface is directed toward the cytoplasmic face of the membrane; strong evidence supports the contention that OCTs can display markedly different affinities for substrate when oriented toward the cytoplasmic versus external aspect of the membrane (49).
Nevertheless, the correlation between predicted and measured binding interactions with a docking surface based upon the homology model of rbOCT2 and the fact that the correlation is based upon the specific sequence information contained in OCT2 (versus the GlpT template structure) support the conclusion that homology modeling offers a powerful tool for probing structure-function interactions of OCTs. Furthermore, in silico tests of the influence of selected residues on binding interactions also offer the opportunity to generate specific testable hypotheses. For example, the homology model suggests that replacement of the glutamate at position 447 with aspartate, i.e. a change in residue volume with no change in charge, is predicted to decrease the affinity of the transporter for TEA, cimetidine, and NBD, with little effect on MPP+ binding, a prediction readily testable in future site-directed studies.
In conclusion, we have identified amino acid residues in OCT2 that play a role in defining the homolog-specific selectivity that distinguishes the OC transporters OCT1 and OCT2. Mutations of one of these (Glu447) had a profound effect on both selectivity and apparent turnover of the transporter. The results have been interpreted in the context of a homology model of the three-dimensional structure of OCT2. A docking surface that included this and other amino acid residues known to influence substrate interactions with OCTs produced calculated binding energies that correlated well with IC50 values for substrate/inhibitor interactions with OCT2. The model of the OCT2 structure offers the opportunity to design experiments to test the influence of specific residues on both binding and translocation properties of this transporter and to develop predictive models of substrate/inhibitor interaction.
| FOOTNOTES |
|---|
The on-line version of this article (available at http://www.jbc.org) contains supplemental Table A and Fig. A. ![]()
1 To whom correspondence should be addressed: Dept. of Physiology, College of Medicine, University of Arizona, Tucson, AZ 85724. Tel.: 520-626-4253; Fax: 520-626-2383; E-mail: shwright{at}u.arziona.edu.
2 The abbreviations used are: OCs, organic cations; OCT, organic cation transporter; TEA, tetraethylammonium; MPP+, 1-methyl-4-phenylpyridinium; MFS, major facilitator superfamily; rb, rabbit; CHO, Chinese hamster ovary; PBS, phosphate-buffered saline; TMH, transmembrane-spanning helix; NBD-MTMA, N,N,N-trimethyl-2-(methyl(7-nitro-2,1,3-benzoxadiazol-4-yl)amino)ethanaminium iodide; PAH, p-aminohippurate; OAT, organic anion transporter; r, rat. ![]()
3 X. Zhang and S. H. Wright, unpublished data. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
O. Zolk, T. F. Solbach, J. Konig, and M. F. Fromm Functional Characterization of the Human Organic Cation Transporter 2 Variant p.270Ala>Ser Drug Metab. Dispos., June 1, 2009; 37(6): 1312 - 1318. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Gorbunov, V. Gorboulev, N. Shatskaya, T. Mueller, E. Bamberg, T. Friedrich, and H. Koepsell High-Affinity Cation Binding to Organic Cation Transporter 1 Induces Movement of Helix 11 and Blocks Transport after Mutations in a Modeled Interaction Domain between Two Helices Mol. Pharmacol., January 1, 2008; 73(1): 50 - 61. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Sturm, V. Gorboulev, D. Gorbunov, T. Keller, C. Volk, B. M. Schmitt, P. Schlachtbauer, G. Ciarimboli, and H. Koepsell Identification of cysteines in rat organic cation transporters rOCT1 (C322, C451) and rOCT2 (C451) critical for transport activity and substrate affinity Am J Physiol Renal Physiol, September 1, 2007; 293(3): F767 - F779. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Chen, S. Zhang, M. Sorani, and K. M. Giacomini Transport of Paraquat by Human Organic Cation Transporters and Multidrug and Toxic Compound Extrusion Family J. Pharmacol. Exp. Ther., August 1, 2007; 322(2): 695 - 700. [Abstract] [Full Text] [PDF] |
||||
![]() |
J.-i. Asaka, T. Terada, M. Tsuda, T. Katsura, and K.-i. Inui Identification of Essential Histidine and Cysteine Residues of the H+/Organic Cation Antiporter Multidrug and Toxin Extrusion (MATE) Mol. Pharmacol., June 1, 2007; 71(6): 1487 - 1493. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. N. Rizwan, W. Krick, and G. Burckhardt The Chloride Dependence of the Human Organic Anion Transporter 1 (hOAT1) Is Blunted by Mutation of a Single Amino Acid J. Biol. Chem., May 4, 2007; 282(18): 13402 - 13409. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. M. Pelis, Y. Dangprapai, T. M. Wunz, and S. H. Wright Inorganic mercury interacts with cysteine residues (C451 and C474) of hOCT2 to reduce its transport activity Am J Physiol Renal Physiol, May 1, 2007; 292(5): F1583 - F1591. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Hong, F. Zhou, K. Lee, and G. You The Putative Transmembrane Segment 7 of Human Organic Anion Transporter hOAT1 Dictates Transporter Substrate Binding and Stability J. Pharmacol. Exp. Ther., March 1, 2007; 320(3): 1209 - 1215. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Zhou, L. Xia, K. Engel, and J. Wang Molecular Determinants of Substrate Selectivity of a Novel Organic Cation Transporter (PMAT) in the SLC29 Family J. Biol. Chem., February 2, 2007; 282(5): 3188 - 3195. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. L. Perry, N. Dembla-Rajpal, L. A. Hall, and J. B. Pritchard A Three-dimensional Model of Human Organic Anion Transporter 1: AROMATIC AMINO ACIDS REQUIRED FOR SUBSTRATE TRANSPORT J. Biol. Chem., December 8, 2006; 281(49): 38071 - 38079. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. M. Pelis, X. Zhang, Y. Dangprapai, and S. H. Wright Cysteine Accessibility in the Hydrophilic Cleft of Human Organic Cation Transporter 2 J. Biol. Chem., November 17, 2006; 281(46): 35272 - 35280. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. M. Pelis, W. M. Suhre, and S. H. Wright Functional influence of N-glycosylation in OCT2-mediated tetraethylammonium transport Am J Physiol Renal Physiol, May 1, 2006; 290(5): F1118 - F1126. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| All ASBMB Journals | Molecular and Cellular Proteomics |
| Journal of Lipid Research | ASBMB Today |