Conventionally, most amino acid substitutions at important protein positions are expected to abolish function. However, in several soluble-globular proteins, we identified a class of non-conserved positions for which various substitutions produced progressive functional changes; we consider these evolutionary “rheostats”. Here, we report a strong rheostat position in the integral membrane protein, Na+/taurocholate cotransporting polypeptide (NTCP), at the site of a pharmacologically-relevant polymorphism (S267F). Functional studies were performed for all 20 substitutions (“S267X”) with three substrates (taurocholate, estrone-3-sulfate and rosuvastatin). The S267X set showed strong rheostatic effects on overall transport, and individual substitutions showed varied effects on transport kinetics (Km and Vmax) and substrate specificity. To assess protein stability, we measured surface expression and used the Rosetta software suite to model structure and stability changes of S267X. Although buried near the substrate binding site, S267X substitutions were easily accommodated in the NTCP structure model. Across the modest range of changes, calculated stabilities correlated with surface-expression differences, but neither parameter correlated with altered transport. Thus, substitutions at rheostat position 267 had wide-ranging effects on the phenotype of this integral membrane protein. We further propose that polymorphic positions in other proteins might be locations of rheostat positions.

AbstractConventionally, most amino acid substitutions at important protein positions are expected to abolish function. However, in several soluble-globular proteins, we identified a class of non-conserved positions for which various substitutions produced progressive functional changes; we consider these evolutionary “rheostats”. Here, we report a strong rheostat position in the integral membrane protein, Na+/taurocholate cotransporting polypeptide (NTCP), at the site of a pharmacologically-relevant polymorphism (S267F). Functional studies were performed for all 20 substitutions (“S267X”) with three substrates (taurocholate, estrone-3-sulfate and rosuvastatin). The S267X set showed strong rheostatic effects on overall transport, and individual substitutions showed varied effects on transport kinetics (Km and Vmax). However, the outcomes were substrate dependent, indicating altered specificity. To assess protein stability, we measured surface expression and used the Rosetta software suite to model structure and stability changes of S267X. Although buried near the substrate binding site, S267X substitutions were easily accommodated in the NTCP structure model. Across the modest range of changes, calculated stabilities correlated with surface-expression differences, but neither parameter correlated with altered transport. Thus, substitutions at rheostat position 267 had wide-ranging effects on the phenotype of this integral membrane protein. We further propose that polymorphic positions in other proteins might be locations of rheostat positions.

Na+/taurocholate cotransporting polypeptide (NTCP) is important for the enterohepatic circulation of bile acids, which has been suggested to contribute to the long serum elimination half-lives of perfluoroalkyl substances in humans. We demonstrated that some perfluoroalkyl sulfonates are transported by NTCP; however, little was known about carboxylates. The purpose of this study was to determine if perfluoroalkyl carboxylates would interact with NTCP and potentially act as substrates. Sodium-dependent transport of [3H]-taurocholate was measured in human embryonic kidney cells (HEK293) stably expressing NTCP in the absence or presence of perfluoroalkyl carboxylates with varying chain lengths. PFCAs with 8 (PFOA), 9 (PFNA), and 10 (PFDA) carbons were the strongest inhibitors. Inhibition kinetics demonstrated competitive inhibition and indicated that PFNA was the strongest inhibitor followed by PFDA and PFOA. All three compounds are transported by NTCP, and kinetics experiments revealed that PFOA had the highest affinity for NTCP with a Km value of 1.8 ± 0.4 mM. The Km value PFNA was estimated to be 5.3 ± 3.5 mM and the value for PFDA could not be determined due to limited solubility. In conclusion, our results suggest that, in addition to sulfonates, perfluorinated carboxylates are substrates of NTCP and have the potential to interact with NTCP-mediated transport.

Pravastatin (PVA) can be isomerized to its inactive metabolite 3'α-iso-pravastatin acid (3αPVA) under acidic pH conditions. Previous studies reported inter-individual differences in circulating concentrations of PVA and 3αPVA. This study investigated the functional consequences of PVA isomerization on OATP1B1-mediated transport. We characterized 3αPVA inhibition of OATP1B1-mediated PVA uptake into HEK293 cells expressing the four different OATP1B1 proteins (*1a, *1b, *5 and *15). 3αPVA inhibited OATP1B1-mediated PVA uptake in all four OATP1B1 gene products but with lower IC 50 /K i values for OATP1B1*5 and *15 than for the reference proteins (*1a and *1b). PVA and 3αPVA were transported by all four OATP1B1 proteins. Kinetic experiments revealed that maximal transport rates (V max values) for OATP1B1 variants *5 and *15 were lower than for *1a and *1b for both substrates. Apparent affinities for 3αPVA transport were similar for all four variants. However, the apparent affinity of OATP1B1*5 for 3αPVA was higher (lower K m value) than for PVA. These data confirm that PVA conversion to 3αPVA can have functional consequences on PVA uptake and impacts OATP1B1 variants more than the reference protein, thus highlighting another source variation that must be taken into consideration when optimizing the PVA dose-exposure relationship for patients. Significance Statement 3'α-iso-pravastatin acid inhibits pravastatin uptake for all OATP1B1 protein types, however, the IC 50 values were significantly lower in OATP1B1*5 and *15 transfected cells. This suggests that a lower concentration of 3'α-iso-pravastatin is needed to disrupt OATP1B1-mediated pravastatin uptake, secondary to decreased cell surface expression of functional OATP1B1 in variant expressing cells. These data will refine previous pharmacokinetic models that are utilized to characterize pravastatin inter-individual variability with an ultimate goal of maximizing efficacy at the lowest possible risk for toxicity.

In the Na+/taurocholate cotransporting polypeptide (NTCP), the clinically relevant S267F polymorphism occurs at a “rheostat position”. That is, amino acid substitutions at this position (“S267X”) lead to a wide range of functional outcomes. This result was particularly striking because molecular models predicted the S267X side chains are buried, and thus, usually expected to be less tolerant of substitutions. To assess whether structural tolerance to buried substitutions is widespread in NTCP, here we used Rosetta to model all 19 potential substitutions at another 13 buried positions. Again, only subtle changes in the calculated stabilities and structures were predicted. Calculations were experimentally validated for 19 variants at codon 271 (“N271X”). Results showed near wildtype expression and rheostatic modulation of substrate transport, implicating N271 as a rheostat position. Notably, each N271X substitution showed a similar effect on the transport of three different substrates and thus did not alter substrate specificity. This differs from S267X, which altered both transport kinetics and specificity. As both transport and specificity may change during protein evolution, the recognition of such rheostat positions may be important for evolutionary studies. We further propose that the presence of rheostat positions is facilitated by local plasticity within the protein structure. Finally, we note that identifying rheostat positions may advance efforts to predict new biomedically relevant missense variants in NTCP and other membrane transport proteins.