|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 282, Issue 33, 23841-23853, August 17, 2007
Structural Variation Governs Substrate Specificity for Organic Anion Transporter (OAT) HomologsPOTENTIAL REMOTE SENSING BY OAT FAMILY MEMBERS*
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
-lactam antibiotics, and diuretics, reviewed in Ref. 4); however, there has not been a comprehensive examination of substrate preferences within a single study or clear mechanistic insight into the structural requirements for transport by OATs.
Oat6 is a close phylogenetic relation of Oat1, first identified as NKT (5, 7, 8).3 Unique among the OATs, it appears absent from the kidney but is instead expressed in olfactory mucosa (9, 10). In addition to different tissue localization, substrate discrimination between the two transporters (10) may reflect their distinctive physiological functions. Both Oat1 and Oat6 interact (though with different affinity) with organic anions implicated in urinary odor-type specification (i.e. those contributing to the urinary odors and thus helping mice distinguish different individuals) (3, 10). This raises the possibility that odorants secreted into the urine via one OAT might be absorbed and/or sensed in the olfactory mucosa via another OAT (10). Furthermore, Oat6 or similar transporters might mediate nasal absorption of some of the organic anionic drugs that are also substrates for other SLC22 family transporters. This could allow for intranasal administration as an alternative route of drug delivery that might bypass the blood-brain barrier, suggesting the possibility of selective drug targeting for enhanced brain delivery (11–13).
To characterize the multispecificity of OATs, two kinds of substrate-transporter interaction should be analyzed separately as follows: 1) interactions leading to the substrate association with the transporter-binding site, and 2) interactions of the bound substrate within the binding site that affect transporter conformational dynamics and substrate translocation. Previous studies by Ullrich (14), performed on whole kidney proximal tubules prior to the cloning and identification of the key transporters expressed at renal as well as extrarenal sites, assumed specific interactions between hydrophobic moieties of the substrate with particular hydrophobic residues within the transporter channel; furthermore, optimal substrates contained two carboxylates spaced 6–7 Å apart. Aside from this work, which is focused on the whole kidney, no comprehensive studies have been carried out to clarify chemical space requirements for substrate specificity. Given the tremendous pharmacological relevance of the OATs, as well as their role in handling key metabolites in many different cells and tissues, it is important to define the structural basis for substrate preferences.
To address this issue, we have performed a detailed wet lab and computational characterization of the substrate binding preferences and transport function of Oat6 in comparison to those of Oat1, testing over 40 compounds for their interaction with and/or transport by these OATs. We first analyzed correlations of the binding preferences of these compounds with their one-dimensional property activities (hydrophobicity, net charge, and molecular mass). Next, we extended these analyses by determining the three-dimensional quantitative structure-activity relationships (QSAR) of Oat1 and Oat6 binding using comparative molecular field analysis (CoMFA). In three-dimensional QSAR, differences in observed biological properties are correlated with determinations made by placing ligands in a three-dimensional field box and evaluating the electrostatic (Coulombic interactions) and steric (van der Waals interactions) fields at regularly spaced grid points. This approach has been used with much success on other transporters (15) and in fact can be used to predict the affinity of a ligand for a transporter.
Our experimental data together with computational analyses strongly support several novel conclusions regarding the structural determinants of OAT substrate selectivity, helping to explain the difference in Oat1 versus Oat6 substrate preferences and presumably reflecting specific differences in the transporter binding pockets. We find substrate potency (pKi) to be correlated with hydrophobicity (logP) mainly in the case of Oat6, suggesting a hydrophobicity-driven association of organic anions with the transporter-binding site. In contrast, substrate binding to Oat1 involves roughly equal contributions from hydrophobic and electrostatic interactions. Differences in substrate selectivity between Oat6 and Oat1 are influenced by ligand size and charge, rather than hydrophobicity, with smaller mono-anions manifesting a higher affinity for Oat6 and di-anions for Oat1. In addition, we determined the maximum uptake rates (Vmax) of six radiolabeled substrates for both transporters. Little correlation was found between Ki and Vmax, suggesting that for these transporters different substrate characteristics influence binding to the transporter versus actual translocation across the membrane. The data suggest ways of targeting drugs to specific OATs and also a mechanism for remote sensing, via endogenous organic anions, between OAT-expressing tissue compartments.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
Unlabeled organic anions (potential substrates/inhibitors) were obtained from Sigma as free acids or sodium salts, and their stock solutions (100 mM) were prepared and adjusted to pH 7.4 as described previously (10). LogP values (octanol/water partition coefficients) of the organic anions were calculated using Molinspiration software.
Uptake in Xenopus Oocytes—Capped RNA was synthesized from linearized plasmid DNA (mOat6, Image clone ID 6309674; mOat1, Image clone ID 4163278) using mMessage mMachine in vitro transcription kit (Ambion, Inc., Austin, TX). Xenopus oocyte uptake assays were performed as described previously (10). Briefly, oocytes were isolated and maintained in Barth's buffer (88 mM NaCl, 1 mM KCl, 0.33 mM Ca(NO3)2, 0.41 mM CaCl2, 0.82 mM MgSO4, 2.4 mM NaHCO3, and 10 mM HEPES, pH 7.4) supplemented with 5% fetal horse serum, 0.05 mg/ml gentamycin sulfate, and 2.5 mM sodium pyruvate. The day after oocytes were harvested, cRNA solution (mOat6 or mOat1, 1 µg/µl) was injected into oocytes (23 nl/oocyte) using the Nanoliter 2000 nanoinjector (World Precision Instruments, Sarasota, FL). Two days after injection, oocytes were washed in serum-free Barth's buffer, and experimental groups of 20–30 oocytes each were placed in wells of a 24-well plate with 1 ml of Barth's buffer containing 1 µCi of a 3H-labeled tracer ion and an unlabeled organic anion (tracer uptake inhibitor), with no inhibitor in a control group. After 1 h of incubation at room temperature, oocytes were washed three times with ice-cold Barth's buffer, and each experimental group was divided into four samples of 5–8 oocytes, and radioactivity was measured by scintillation counting.
Transport activity was calculated for [3H]-labeled tracers as tracer clearance from the incubation medium (CL = Vtransport/S), by dividing the amount of tracer absorbed per oocyte per unit time (Vtransport, cpm/oocyte/h) by the tracer concentration in the incubation medium (S, cpm/µl). In all experiments, the background (non-OAT-mediated) tracer uptake was measured in uninjected oocytes (in preliminary experiments, uptake in water-injected and uninjected oocytes was not found to be significantly different). This background uptake (probably combining nonsaturable tracer diffusion and endogenous transport expressed by the oocytes) was subtracted from the uptake measured in injected oocytes to calculate the OAT-mediated component of uptake. In a typical experiment, the control (uninhibited) transporter-mediated tracer clearance was 0.35 to 0.9 µl/oocyte/h for [3H]ES uptake in Oat6-injected oocytes and 0.7–2.0 µl/oocyte/h for [3H]PAH uptake in Oat1-injected oocytes. This indicates that the percentage of tracer that is transported in 1 h did not exceed 6% (assuming maximum experimental group size of 30 oocytes and maximum tracer clearance of 2.0 µl/oocyte/h), so that the tracer concentration was practically constant during the incubation time. In all inhibition experiments, the uninhibited OAT-mediated clearance exceeded the background clearance of the same tracer in uninjected oocytes by at least 10-fold.
Calculations and Statistics—In each experimental group, tracer clearance was calculated as mean ± S.E. of quadruplicate samples (with exception of the control group that included eight samples). In uptake inhibition experiments, the OAT-mediated clearance (i.e. difference between cRNA-injected and uninjected oocytes) in the presence of an inhibiting organic anion was expressed as a percentage of the mean OAT-mediated clearance in the control group.
To determine potencies of organic anions (tested as tracer uptake inhibitors), tracer uptake was measured in the presence of increasing inhibitor concentration. For each organic anion, a series of 3–4 concentrations in successive 10-fold increments was tested. The inhibition data were considered sufficient for curve-fitting when the inhibitor concentration points spanned the interval comprising 50% inhibition. The tracer uptake versus log[inhibitor] data were fit in the one-site competition equation incorporated in Prism 4.0 software (GraphPad Inc., San Diego, CA), and log(IC50) was calculated as mean ± S.E. The IC50 value was calculated as 10mean and standard error for IC50 as S.E.IC50 = 10mean – 10mean – S.E.. Analysis of the uptake inhibition of 3H-labeled tracer substrate (237 nM [3H]PAH in Oat1-injected oocytes or 17 nM [3H]ES in Oat6-injected oocytes) by "cold" PAH and ES yielded Kd values of 9.4 ± 1.1 and 58 ± 10 µM, respectively (10). These results indicate that the concentrations of both tracers in our uptake inhibition experiments were far below their Kd values, so that analyzing inhibition data using the Cheng-Prusoff equation (Ki = IC50/(1 + [substrate]/Kd)) for all organic anions yielded Ki values not distinguishable from respective IC50 values (the difference being always significantly less than the standard error for IC50; data not shown). Based on these results, we used the calculated values of IC50 as inhibition constant (Ki) values and calculated organic anion potencies as pKi =–log(Ki).
The maximum transporter-mediated uptake rate was determined for six 3H-labeled substrates in Oat6- and Oat1-injected oocytes as Vmax = CL x (S + Ki) (see under "Results" for details). Standard error for Vmax was calculated as S.E.Vmax = (S.E.CL2 x (S + Ki)2 + S.E.Ki2 x CL2)1/2.
Molecular Modeling and Alignment—Three-dimensional structure building was performed using SYBYL 7.1 (Sybyl). Energy minimizations were performed using the Tripos force field (17) and Gasteiger-Hückel charges with distance-dependent dielectrics and the conjugate gradient method with a convergence criterion of 0.001 kcal/mol. The most important requirement for CoMFA (18) is the structural overlap of the molecules to be analyzed to a suitable template, which is assumed to be a "bioactive conformation." The three-dimensional coordinates of methotrexate were obtained from the Protein Data Bank (access code 1RX3) (19) and used as a template for the superimposition of all other molecules using the "Flexible Superposition" option in the FlexS (20) module of Sybyl.
CoMFA Three-dimensional QSAR Models—CoMFA explains the gradual changes in observed biological properties by evaluating the electrostatic (Coulombic interactions) and steric (van der Waals interactions) fields at regularly spaced grid points surrounding a set of mutually aligned ligands for a specific target protein. A statistical algorithm, partial least square (PLS), was used to correlate the field descriptors with biological activities (i.e. pKi). The standard CoMFA settings were applied for developing Oat1, Oat6, and selectivity models (i.e. pKi(Oat6) – pKi(Oat1); see Table 3). The global model contains 28 molecules in training set (19 mono-anions and 9 di-anions) and 7 molecules in test set (5 mono-anions and 2 di-anions). The pKi values of 2-hydroxy-3-methyl butyrate, 2-methyl butyrate, 3-hydroxy butyrate, and 2-ethyl hexanoate were tested using racemic mixtures; hence, their activities were predicted as separate enantiomers. The aligned molecules were placed inside a three-dimensional cubic lattice box with grid spacing of 2.0 Å in x, y, and z directions. The CoMFA steric and electrostatic descriptors were calculated at each grid point using an sp3 hybridized carbon atom probe with +1 charge at 1.52 Å van der Waals radius. The cutoff value for steric and electrostatic fields was set to +30 kcal/mol. To derive CoMFA models, pKi values were used as the dependent variables, and CoMFA (steric and electrostatics) and logP values were used as the independent variables in PLS analysis. The predictive ability of the model was evaluated using leave-one out (LOO) cross-validation studies (r2cv) and the predictive residual sum of squares (PRESS) of the test set. The optimum number of components from r2cv was used for the final PLS model and calculating the conventional r2 of the model. The r2cv was calculated using Equation 1,
![]() |
|
| RESULTS |
|---|
|
|
|---|
Fig. 1 shows the side-by-side comparison of the inhibition of Oat6- and Oat1-mediated uptake by different organic anions tested at a single concentration (in the range of 1 µM to 1 mM, depending on the organic anion potency). The majority of the organic anions in the study significantly inhibited radioactive tracer uptake through Oat6 and/or Oat1, when tested at up to 1 mM concentration. Of all organic anions tested, 16 displayed significantly stronger inhibition of Oat6, 15 inhibited Oat1 to a greater extent, and 10 showed no apparent preference for either transporter.
To determine the potency of organic anions that exhibited submillimolar activity upon at least one of the two transporters, we examined the concentration-dependent inhibition of tracer uptake. Organic anions producing <50% inhibition of both Oat1- and Oat6-mediated transport at 1 mM (ascorbate, maleate, tartrate, taurocholate, cholate, and gluconate) (Fig. 1) were excluded from the potency assessment experiments. With the remaining 35 anions, concentration-dependent inhibition of tracer uptake was measured, and the potencies (apparent pKi values, see under "Experimental Procedures" for details) were calculated for both transporters. Fig. 2 shows representative examples of the concentration-dependent inhibition of tracer uptake by organic anions selective for Oat6 (benzoate and pyruvate) or for Oat1 (PAH and fluorescein), and by those displaying similar potency for the two transporters (probenecid and penicillin G). The summarized results, along with the structures of the organic anions, their molecular masses, and logP values (octanol/water partition coefficient, indicating hydrophobicity), are presented in Table 1.
|
|
|
Physicochemical and Structural Determinants of the Oat6/Oat1 Affinity and Selectivity of Organic Anions—The data in Table 1 were analyzed to determine the one-dimensional quantitative property-activity relationship between the potency of organic anions for the OATs and the most common molecular features that impact the interaction of a small molecule with protein-binding sites (molecular mass, hydrophobicity, and net charge).
The potencies of the organic anions (pKi values) plotted versus their molecular masses are presented in Fig. 4, A and B, as determined for Oat6 and Oat1, respectively. Although no correlation between organic anion potency and molecular mass was detected for mono-anions on Oat6 or for dianions with Oat1, a weak correlation was observed for di-anions on Oat6 (Fig. 4A) and a strong correlation for mono-anions with Oat1 (Fig. 4B) (in both cases, smaller anions tend to manifest lower potency). To verify the influence of the molecular mass of organic anions on their selectivity between the two transporters, we also analyzed the potency difference between the transporters, pKi(Oat6) – pKi(Oat1). This approach compensates for the structural features of organic anions that contribute equally to their interaction with both Oat1 and Oat6 and reveals the features distinguishing between the two transporters. Single-charged organic anions, but not di-anions, showed a very strong correlation between the potency difference and molecular mass (p < 0.0001; Fig. 4C), with a distinct Oat6 preference for mono-anions with a molecular mass below or around 150.
We also investigated whether organic anion hydrophobicity (logP value) correlated with differences in potency between Oat6 and Oat1. In contrast to the data on the effect of organic anion molecular mass on selectivity between Oat6 and Oat1 (Fig. 4, A–C), the organic anion logP value apparently affected the potency for both Oat6 (strong potencylogP correlation; Fig. 4D) and Oat1 (weak correlation; Fig. 4E) but not the selectivity between the transporters (no correlation between pKi(Oat6) – pKi(Oat1) and logP; Fig. 4F). Interestingly, di-anions displayed a consistently lower potency for Oat6 than mono-anions with comparable logP (Fig. 4D), whereas for Oat1, di-anions tended to be more potent than mono-anions (Fig. 4E). Computational analysis later strikingly confirmed these differences in substrate preferences (see global and anion-specific CoMFA models below).
|
Uptake of Radiolabeled Tracers in Oat6- and Oat1-injected Oocytes—To directly assess organic anion transport via Oat6 and Oat1, we measured the rate of transporter-mediated uptake for a subset of six organic anions (3H-labeled) varying in their potencies and transporter binding preferences (Table 1) as follows: PAH, ES, ochratoxin A, prostaglandin E2, methotrexate, and ibuprofen. The entire set of measurements was performed in a single experiment using Oat1-injected, Oat6-injected, and uninjected oocytes from the same batch to improve the likelihood that transporter expression level would not contribute to uptake differences between substrates. The results (Fig. 5A) revealed significant Oat6-mediated uptake of [3H]ES and Oat1-mediated uptake of [3H]PAH and [3H]ochratoxin A that considerably exceeded (typically by more than 10-fold) the background uptake of these compounds in uninjected oocytes. Lesser transporter-mediated uptake was measured for [3H]ochratoxin A and [3H]prostaglandin E2 in Oat6-injected oocytes and for [3H]prostaglandin E2 and [3H]methotrexate in Oat1-injected oocytes. The high background uptake observed for the most hydrophobic of these six organic anions, ibuprofen, is likely because of its nontransporter-mediated diffusion across the membrane.
Transport Efficacy (Maximum Translocation Rate) for Different Oat6/Oat1 Substrates—In our experiments, uptake of 3H-labeled tracers was measured as the initial rate on the linear part of the uptake kinetics (10). Under these conditions, the dependence of transporter-mediated uptake rate (Vtransport) on substrate concentration (S) can be described as Vtransport = Vmax x S/(S + Km), where Km is the apparent equilibrium constant of the substrate-transporter complex, and Vmax is the maximum uptake rate.
Oat6- and Oat1-mediated clearance of six 3H-labeled substrates (ES, PAH, prostaglandin E2, ibuprofen, ochratoxin A, and methotrexate) (Fig. 5B) was determined by subtracting the background clearance (in uninjected oocytes) from the total clearance observed in Oat6- and Oat1-injected oocytes (Fig. 5A). Based on obtained values of CL = Vtransport/S, the maximum uptake rates, Vmax, can be calculated as Vmax = (Vtransport/S) x (S + Km), provided that the respective Km values are known. It is reasonable to assume that the Ki values found for these six organic anions in the tracer uptake inhibition experiments (Table 1) should represent an approximation of the Km values. Accordingly, the maximum uptake rate for each substrate was calculated as Vmax = (Vtransport/S) x (S + Ki) = CL x (S + Ki). The discussion of Vmax that follows should therefore be understood to be qualified by the above assumption.
These calculations (Table 2) revealed that the apparent Vmax varies greatly among the organic anions tested (Fig. 5C), and for both Oat1 and Oat6, there is little correlation between organic anion potency (Ki) and efficacy (Vmax). For Oat6, ES is the most efficacious substrate, followed by ochratoxin A and PAH. With Oat1, the highest Vmax value was calculated for methotrexate followed by ochratoxin A, PAH, and ES.
|
CoMFA, Global Models—To further delineate the role of hydrophobic (steric) interactions in substrate affinity for Oat1 and Oat6, we applied CoMFA to analyze the relative contributions of steric and electrostatic molecular properties in three-dimensional space. Furthermore, logP was used as an independent variable to determine its relative contribution to steric properties, if any. The global CoMFA includes both mono- and di-anions. Six different CoMFA models were built. Models 1–3 used CoMFA as independent and pKi as dependent variables to develop models for Oat1, Oat6, and Oat1/Oat6 selectivity, respectively. Models 4–6 included pKi as dependent and CoMFA and logP as independent variables (Table 3). The predictive ability of all models was evaluated based on PRESS, where a lower value is indicative of a better model.
|
50% steric and 50% electrostatic molecular field contributions (Table 3). Although the inclusion of logP as an independent variable (models 4–6) produces robust CoMFA models, the relative contribution to the models is primarily drawn from steric interaction fields. This is not surprising, because logP is an intrinsically steric property. Thus, we have extended and superseded the prior one-dimensional logP models (Fig. 4) by analyzing both steric and electrostatic interactions that govern Oat1 and Oat6 inhibition.
|
Correlation between calculated and experimental pKi values (Fig. 7, A–C) indicates that all models were linear with at least a 5–6 log unit spread in data points; randomly selected test set molecules were evenly distributed across the range in activity data for all three models. Further inspection of residual values for both training and test set data (Fig. 7D) revealed that the mean residual values for Oat1 and Oat6 training models were –0.04 and –0.01, respectively; the Oat6/1 selectivity training model had the highest mean residual value (0.19). As expected for studies involving extrapolation of data, test set mean residual values were markedly higher for Oat1 (–0.43) and Oat6/1 selectivity model (0.32); however, Oat6 test set predictions were identical to the performance of the training set estimates (–0.01), indicating a significantly predictive model.
Comparative Molecular Field Analysis, Anion-specific Models—Based on the global CoMFA model, it is clear that logP is an important physicochemical property governing affinity for Oat6. In line with one-dimensional property-activity relationship analyses (Fig. 4), individual CoMFA models were built for mono- and di-anions to investigate the relationship between specific anions and logP for OAT binding. The entire data set was divided into mono- and di-anions, and using pKi as dependent variable and CoMFA and logP as independent variables, three individual models were developed for Oat1, Oat6, and Oat6/Oat1 selectivity (Table 4). The predictive ability of the Oat1 mono- and di-anion models is comparable as indicated by similar r2 values. However, the predictive ability of the Oat6 mono-anion model (r2 = 0.691) is poor compared with that of the di-anion model (r2 = 0.833). Overall, the individual mono- and di-anion models clearly indicate that equal contributions of steric and electrostatic factors govern affinity for Oat1, whereas logP contributions predominate affinity for Oat6.
|
| DISCUSSION |
|---|
|
|
|---|
Clearly defining the similarities and differences in substrate structural preferences between these two transporters then could have important ramifications not only for drug design (e.g. inhibition to prolong half-life or synthesis of drugs that go through a specific transporter), but for understanding how tissues communicate via small molecule within the body and between individuals. Our results indicate very different charge and hydrophobicity preferences for these two transporters. Because they are such close phylogenetic relatives, our data suggest that although the general three-dimensional structures of the transporters might be very similar, the binding pockets have crucial differences.
Characterization of the Affinity of Organic Anions for Oat6 and Oat1—To investigate the molecular features influencing organic anion selectivity between the transporters, we used the potencies of 35 organic anions determined for both Oat6 and Oat1 (Table 1). Our results for Oat1 are consistent with the Km and Ki data available elsewhere for many of these organic anions (reviewed in Refs. 4, 14). Our study includes four of the total of five organic anions assayed in other studies (28) (ES, probenecid, salicylate, and penicillin G), and our results show general agreement.
Comparative analysis of organic anion potencies (pKi values, see Table 1) reveals the influence of molecular features such as anion net charge, molecular mass, and hydrophobicity on affinity and selectivity for Oat6 and Oat1. The most obvious observation is a clear-cut preference of single-charged organic anions for Oat6 and of double-charged ones for Oat1 (Figs. 3 and 4).
Examination of the organic anion potency correlation with molecular mass also shows a difference in the binding preferences of the two transporters. With single-charged organic anions, no correlation was observed for Oat6 (Fig. 4A), whereas Oat1 tends to "dislike" smaller ligands (Fig. 4B). This tendency becomes more pronounced when the potency difference, pKi(Oat6) – pKi(Oat1), is analyzed (Fig. 4C). Thus, Oat6 would seem to represent a better fit, compared with Oat1, for binding of low molecular mass odorant anions. Interestingly, we did not see a similar correlation for double-charged organic anions (Fig. 4C).
|
Compared with Oat6, the binding of organic anions to Oat1 displays a much weaker correlation with logP, suggesting a significant role for structure-specific molecular interactions in the ligand-transporter binding energy. Moreover, although a second negatively charged group consistently lowers organic anion potency for Oat6 (Fig. 4D), an apparently positive contribution of a second anionic group to potency is noted for Oat1 (Fig. 4E). This overall pattern suggests a difference between the two transporters in anion-binding modes as follows: primarily hydrophobicity-driven binding for Oat6, with relatively little involvement of specific molecular interactions (the present study), versus largely structure-dependent ligand association with Oat1 this study and Ref. 14. Fig. 4F clearly demonstrates that the ligand hydrophobicity equally contributes to the binding energy with both transporters; no correlation is detected between ligand hydrophobicity and potency difference, pKi(Oat6) – pKi(Oat1), which is calculated to factor out nontransporter-specific aspects of ligand-transporter interactions and uncover the role of transporter-specific interactions.
In the case of a more homogeneous collection of organic anions (aliphatic mono- and di-carboxylates) lacking additional groups (other than carboxyls) capable of specific interactions with hydrogen bond-forming or aromatic groups of the transporter, it is noteworthy that a strong potency-hydrophobicity correlation is observed for inhibition of proximal tubular PAH transport (which is largely mediated by Oat1 (3)) (Fig. 4G, based on the data of Ullrich and co-workers (23, 24)); a very strong correlation observed for mono-carboxylates (R2 = 0.917, p < 0.0001) supports the idea that the single carboxyl group in organic anions varying in size (C3 to C10) and molecular structure (linear, branched, and unsaturated) interacts uniformly with the binding site of the transporter, with the overall hydrophobicity of the ligand controlling its equilibrium between the aqueous medium and the transporter binding pocket. With dicarboxylates, a weaker pKi – logP correlation is noted (R2 = 0.484, p < 0.01), which is consistent with specific interactions of both carboxyl groups being essential, so that their relative position affects affinity. Interestingly, no significant correlation was obtained when analyzing a larger data set from the same group of investigators (29–33), which included structurally diverse single charged organic anions containing additional groups (aromatic, hydrogen bond-forming) that are potentially capable of specific interactions with the transporter-binding site (data not shown).
CoMFA—To gain further understanding of the binding mode of organic anions to Oat1 and Oat6, we used a computational approach to model in vitro inhibition data. Although a one-dimensional property-activity relationship is a useful approach to analyze the influence of physicochemical properties on biological activity, it does not take into consideration the three-dimensional molecular orientation and the spatial distribution of electrostatic and steric interactions. To address these limitations, we extended our analyses using the CoMFA three-dimensional QSAR algorithm. Using an iterative process, predictive QSAR models were derived for the individual OAT transporters as well as a model focusing on Oat6/1 selectivity (Table 3; Fig. 6). The latter model is expected to predict compounds that specifically bind to Oat6. As a further illustration, Fig. 7 displays the residual plots of divergence between predicted and actual activity values for the Oat1, Oat6, and Oat6/1 models, respectively. The CoMFA models feature robust r2cv values, indicative of internally consistent models.
|
|
OATs belong to the major facilitator superfamily of transporters currently postulated to function via a single binding site-alternate access mechanism (6, 16). Based on this transport mechanism, the inwardly directed transporter-mediated uptake of a substrate anion can be described as involving the following steps (Fig. 8): 1) substrate binding to the transporter in the conformation open to extracellular medium (Co); 2) switching of the substrate-transporter complex to the conformation open to intracellular medium (Co
Ci); 3) substrate dissociation into the intracellular medium; and 4) switching of the "empty" transporter back to the Co conformation. In this scheme, the equilibrium constant of the substrate-transporter complex determines the substrate concentration profile (Km value, see Fig. 8, step 1), and the rates of conformational transition of the substrate-transporter complex (step 2) and of the "empty" transporter (step 4) control the maximum uptake rate (Vmax).
Substantial differences among the Vmax values found for the substrates tested in this work suggest that the rate of the conformational switch of the transporter-substrate complex (step 2, substrate-dependent), rather than that of empty transporter (step 4, substrate-independent), mostly determines the substrate uptake rate. It should be noted that the empty transporter might actually mean a complex of the transporter with one of the anionic metabolites that are intracellular OAT substrates. In this case, an intracellular substrate (such as
-ketoglutarate for Oat1 (4)) would be transported from the oocytes in exchange for the extracellular substrate taken up. This notion, however, does not affect the conclusion (based on the present view of single binding site-alternating access transport model) about the Co
Ci conformational switch of the transporter-substrate complex (Fig. 8, step 2) being the rate-limiting step in substrate uptake.
In contrast to the observed correlation between the potencies of organic anions for OATs (pKi values, describing the substrate affinity for the transporter-binding site) and their general molecular characteristics (molecular size, charge, and hydrophobicity) (Fig. 4), little correlation was found for the efficacies of organic anions (Vmax values, reflecting the conformational dynamics of the transporter-substrate complex). It should also be emphasized that our results (Table 2) revealed no apparent relation between the Vmax values for Oat6 and Oat1 or between the Vmax and Ki values for either of the transporters. For example, ES and ochratoxin A interact with Oat6 with approximately equal potency, but the former is translocated across the membrane with much higher efficacy. As a counter example, the potency of PAH for Oat1 is much higher than that of ES, but when bound both appear equally efficacious in inducing the transporter conformation switch.
Thus, there appear to be very distinct requirements for the two steps of OAT-mediated transport; substrate binding to the transporter (Fig. 8, step 1) is largely dependent on general physicochemical characteristics of the substrate (charge, size, and hydrophobicity) (14), whereas actual translocation across the membrane, driven by switching the conformation of the substrate-transporter complex (step 2), is apparently determined by distinctive molecular features of the substrate. Consequently, when describing the facilitated transport of a substrate, two kinds of the substrate-transporter interaction should be analyzed separately as follows: 1) interactions leading to the substrate association with the transporter-binding site, and 2) interactions of the bound substrate within the binding site that affect transporter conformational dynamics. Our results indicate marked differences between these two kinds of interaction.
| FOOTNOTES |
|---|
The on-line version of this article (available at http://www.jbc.org) contains supplemental Table 1. ![]()
1 To whom correspondence should be addressed: University of California, San Diego, 9500 Gilman Dr., La Jolla, CA 92093.
2 The abbreviations used are: OAT, organic anion transporter; CoMFA, comparative molecular field analysis; ES, estrone sulfate; PAH, p-aminohippurate; QSAR, quantitative structure-activity relationship; PLS, partial least square; CoMFA, comparative molecular field analysis. ![]()
3 C. E. Lopez-Nieto, GenBankTM accession number MMU528. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
W. Wu, M. E. Baker, S. A. Eraly, K. T. Bush, and S. K. Nigam Analysis of a large cluster of SLC22 transporter genes, including novel USTs, reveals species-specific amplification of subsets of family members Physiol Genomics, July 1, 2009; 38(2): 116 - 124. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. Vallon, S. A. Eraly, W. R. Wikoff, T. Rieg, G. Kaler, D. M. Truong, S.-Y. Ahn, N. R. Mahapatra, S. K. Mahata, J. A. Gangoiti, et al. Organic Anion Transporter 3 Contributes to the Regulation of Blood Pressure J. Am. Soc. Nephrol., September 1, 2008; 19(9): 1732 - 1740. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. M. Truong, G. Kaler, A. Khandelwal, P. W. Swaan, and S. K. Nigam Multi-level Analysis of Organic Anion Transporters 1, 3, and 6 Reveals Major Differences in Structural Determinants of Antiviral Discrimination J. Biol. Chem., March 28, 2008; 283(13): 8654 - 8663. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. A. Zager, A. C. M. Johnson, M. Naito, and K. Bomsztyk Maleate nephrotoxicity: mechanisms of injury and correlates with ischemic/hypoxic tubular cell death Am J Physiol Renal Physiol, January 1, 2008; 294(1): F187 - F197. [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 |