## Abstract

Prokaryotic and eukaryotic Na^{+}-driven transporters couple the movement of one or more Na^{+} ions down their electrochemical gradient to the active transport of a variety of solutes. When more than one Na^{+} is involved, Na^{+}-binding data are usually analyzed using the Hill equation with a non-integer exponent *n*. The results of this analysis are an overall *K _{d}*-like constant equal to the concentration of ligand that produces half saturation and

*n*, a measure of cooperativity. This information is usually insufficient to provide the basis for mechanistic models. In the case of transport using two Na

^{+}ions, an

*n*< 2 indicates that molecules with only one of the two sites occupied are present at low saturation. Here, we propose a new way of analyzing Na

^{+}-binding data for the case of two Na

^{+}ions that, by taking into account binding to individual sites, provides far more information than can be obtained by using the Hill equation with a non-integer coefficient: it yields pairs of possible values for the Na

^{+}affinities of the individual sites that can only vary within narrowly bounded ranges. To illustrate the advantages of the method, we present experimental scintillation proximity assay (SPA) data on binding of Na

^{+}to the Na

^{+}/I

^{−}symporter (NIS). SPA is a method widely used to study the binding of Na

^{+}to Na

^{+}-driven transporters. NIS is the key plasma membrane protein that mediates active I

^{−}transport in the thyroid gland, the first step in the biosynthesis of the thyroid hormones, of which iodine is an essential constituent. NIS activity is electrogenic, with a 2:1 Na

^{+}/I

^{−}transport stoichiometry. The formalism proposed here is general and can be used to analyze data on other proteins with two binding sites for the same substrate.

## INTRODUCTION

Membrane transport proteins serve a key function in living cells by mediating the transport of many different solutes, such as ions, nutrients, and neurotransmitters across biological membranes. Numerous pathological conditions—cystinuria (Fotiadis et al., 2013), glucose–galactose malabsorption (Wright, 2013), hypothyroidism caused by iodide transport defects (Levy et al., 1998; De la Vieja et al., 2007; Paroder-Belenitsky et al., 2011; Portulano et al., 2014)—result from impairments in the function of molecules of this class. Membrane transporters also play a crucial role in the delivery of chemotherapeutic agents and internal radiation (Klutz et al., 2011; Ho et al., 2013) to cancer cells.

Many prokaryotic and eukaryotic transporters exploit the electrochemical gradient of Na^{+} to translocate their substrates across the membrane. Of those, a significant fraction binds and transports more than one Na^{+} ion per transport cycle (Caplan et al., 2008; Shi et al., 2008; Abramson and Wright, 2009). Binding assays using a radiotracer in scintillation proximity assays (SPAs) have become a routine procedure for investigating the function of detergent-solubilized transport proteins and, critically, for determining the substrate-binding constants of purified transporters (Pirch et al., 2002; Berry and Price-Jones, 2005; Quick and Javitch, 2007; Harder and Fotiadis, 2012; Khafizov et al., 2012; Zhou et al., 2014).

Conventionally, when more than one Na^{+} is involved, SPA data are analyzed using the Hill equation (Hill, 1910, 1913; Weiss, 1997; Goutelle et al., 2008) with a non-integer exponent *n*:*f*_{n} is the fraction of transporter molecules that are occupied by Na^{+}, and *K _{d}* is an apparent microscopic dissociation constant. Cases with

*n*> 1 are instances of positive cooperativity, and the resulting curves for

*f*

_{n}versus [Na

^{+}] are sigmoidal; the

*K*value yields the [Na

_{d}^{+}] at half-saturation. This equation, however, does not describe a true equilibrium situation: there is no real binding reaction that can lead to this expression, as it would have a fractional stoichiometric coefficient. In the case of two binding sites, for example, a fractional Hill coefficient indicates that proteins with only one site occupied represent a significant fraction of the molecules in equilibrium, especially at low [Na

^{+}].

Here, we propose a novel way of analyzing SPA data obtained from experiments with transporters that bind Na^{+} at two sites. The method is exemplified by analyzing experimental SPA data from the Na^{+}/I^{−} symporter (NIS). This approach, which takes into account binding to individual sites, provides information not accessible with the traditional Hill equation analysis.

## MATERIALS AND METHODS

### Protein purification

293-H cells grown in suspension were transduced with a lentiviral vector containing a rat NIS cDNA construct with an HA tag at the N terminus and a His_{8} tag and a streptavidin-binding protein (IBA GmbH) tag at the C terminus. I^{−} transport levels were comparable to those in FRTL-5 cells, a line of rat thyroid–derived cells that express NIS endogenously (Dai et al., 1996). A membrane fraction prepared from 293 cells was solubilized with 40 mM *n*-dodecyl-β-d-maltopyranoside (DDM; Anatrace) and purified via Strep-Tactin affinity chromatography. The S353A/T354A NIS mutant was expressed and purified using the same procedures.

### Binding experiments

Binding experiments were performed using the SPA technique (Quick and Javitch, 2007). Affinity-purified NIS (WT or mutant) was bound to Cu^{2+} chelate YSi scintillation SPA beads (PerkinElmer) via the C-terminal His_{8} tag. 250 µg SPA beads was used per assay, in a volume of 100 µl. Experiments consisted of competing bound ^{22}Na^{+} with cold Na^{+} at the specified concentrations. The assay buffer consisted of 100–600 mM Tris/Mes, pH 7.5, 20% glycerol, 0–500 mM NaCl (equimolar replacement with Tris/Mes), 1 mM TCEP, and 0.1% DDM with 250 ng of affinity-purified NIS. [^{22}Na]Cl (1,017 mCi/mg; PerkinElmer) and the indicated concentrations of cold NaCl were added to the bead solution simultaneously with the protein. Nonspecific binding was determined in the presence of 800 mM imidazole. All binding assays were performed in white-walled, clear-bottomed 96-well plates and quantitated in a photomultiplier (Wallac; PerkinElmer) tube counter (MicroBeta; PerkinElmer) after overnight equilibration.

### Analysis of NIS SPA data

All readings were background-subtracted and then normalized using the following relation:*f _{x}* is the fraction of NIS molecules for which the bound

^{22}Na

^{+}has been displaced by cold Na

^{+}at a [Na

^{+}] = X; and C

_{0}, C

_{X}, and C

_{500}are the counts per minute at cold [Na

^{+}] = 0, [Na

^{+}] = X, and [Na

^{+}] = 500 mM.

## RESULTS AND DISCUSSION

### Non-integer Hill coefficients

The equation proposed by Hill (Hill, 1910, 1913; Dahlquist, 1978; Weiss, 1997; Goutelle et al., 2008) can be written as*K _{a}* is the equilibrium constant and

*n*is a whole number > 1” (Hill, 1913), and Y is the fraction of protein molecules that have

*n*ligands bound.

This reaction and its associated equilibrium constant imply that the protein (P) has *n* binding sites, and the ligand (L) must bind to the *n* sites simultaneously: binding of L to a single site of the protein does not take place. The exponent in this equation is called the Hill coefficient.

In many cases, a better fit to the experimental data can be achieved by using a non-integer Hill coefficient, yielding the equation,*n* > 1 indicates positive cooperativity; an *n* < 1 indicates negative cooperativity.

A more general case would be one in which the protein has two nonequivalent sites and L can bind to either site first and then bind to the other. The two sites will, in general, have different affinities for the ligand, and the affinity for binding to a second site may be enhanced (positive cooperativity) or reduced (negative cooperativity) when the other site is already occupied. The reaction scheme for the general case involving two sites, A and B, is represented in Fig. 1, where *K _{a,A}* and

*K*are the association constants for the binding of L to the empty protein,

_{a,B}*K*is the association constant for binding L to site B when site A is occupied, and

_{a,A−B}*K*is the association constant for binding L to site A when site B is occupied. Note that only three of the four constants are needed because the upper and the lower paths lead to the same species, and therefore

_{a,B−A}*K*·

_{a,A}*K*=

_{a,A−B}*K*·

_{a,B}*K*. Introducing ϕ, a unitless enhancement/reduction factor, we can write

_{a,B−A}*K*=

_{a,A−B}*K*· ϕ and

_{a,B}*K*=

_{a,B−A}*K*· ϕ. For positive cooperativity, ϕ > 1, and ϕ < 1 corresponds to negative cooperativity (ϕ = 1 means independent sites).

_{a,A}The fraction of doubly occupied molecules can be written as

The normalized signal from an experimental measurement sig(L) is given by*sr* is the ratio of the signal when one L is bound to that when two L’s are bound.

The values of *K _{p}* and

*K*can be obtained by a least-squares fit of normalized experimental data. Eqs. 3 and 4 can then be used to solve for

_{x}*K*and

_{a,A}*K*as a function of ϕ. This yields a quadratic equation with solutions

_{a,B}*K*and the minus sign for

_{a,A}*K*. Because

_{a,B}*K*and

_{a,A}*K*must be real numbers, the expression within the square root “

_{a,B}*sq*”

Three cases have to be considered: ϕ = 1, ϕ > 1, and ϕ < 1. The case ϕ = 1 corresponds to two independent sites (no cooperativity). If *sq* = 0, *K _{a,A}* =

*K*= 1/2

_{a,B}*K*and corresponds to two independent identical sites; otherwise,

_{p}*K*>

_{a,A}*K*.

_{a,B}The case ϕ > 1 corresponds to positive cooperativity. The condition * _{min,p}* ≤ ϕ < ∞. In contrast, the case ϕ < 1 corresponds to negative cooperativity. The condition

*≤ ϕ < 1.*

_{min,p}In both cases, when ϕ = ϕ* _{min}*,

*sq*= 0; therefore,

In the case of positive cooperativity, there is another limiting case. It corresponds to the situation where ϕ is very large, making *K _{a,A}* approach

*K*and

_{p}*K*approach 0 (from Eq. 7a). This case corresponds to absolute sequential binding (i.e., the ligand does not bind to site B unless site A is already occupied).

_{a,B}After the values of *K _{p}* and

*K*are obtained by fitting the experimental data, combinations of values of

_{x}*K*and

_{a,A}*K*can be obtained as a function of the allowed values of ϕ. Although the information obtained is only a range of possible values for the affinities of sites A and B of the empty protein for the ligand,

_{a,B}*K*=

_{a,B}*K*–

_{p}*K*.

_{a,A}This analysis provides more information than just the single value that is obtained from the Hill equation using non-integer exponents. Furthermore, the values of *K _{a,A}*,

*K*, and

_{a,B}*K*, is similar to the

_{x}^{−1/2}*K*obtained from the Hill equation.) One issue that may seem puzzling is that

_{d}*K*approaches zero (and

_{a.B}*K*approaches ∞) as

_{d,B}*K*approaches

_{a,A}*K*. Why does it still correspond to a two-site situation? The answer is simple. Both

_{p}*K*and ϕ are a function of (

_{a,B}*K*−

_{p}*K*);

_{a,A}*K*is directly proportional and ϕ is inversely proportional, such that

_{a,B}*K*approaches zero,

_{a,B}*K*· ϕ remains bound, and the second site has a finite dissociation constant when site A is occupied. Because

_{a,B}*K*varies only within a factor of 2,

_{a,A}*K*· ϕ also varies within a factor of 2.

_{a,B}Because the equations were derived as a function of [L], the concentration of free ligand, they work directly in cases in which [L] is very close to L_{added}, the concentration of ligand added at each experimental point. SPA, in which the amount of protein is much less than the amount of ligand, and fluorescence measurements performed with concentration of protein [P] << *K _{d}* satisfy this condition.

In the case of SPA experiments with ^{22}Na (L = Na^{+}), the transporter with two Na^{+} ions bound and the two species with one Na^{+} bound all produce a signal. Because the signal from the molecules with one Na^{+} is half that from the fully occupied transporter (*sr* = 0.5 in Eq. 6), the normalized SPA signal is proportional to the sum of the fraction of molecules with two Na^{+} ions bound and half the fraction of molecules with one Na^{+} bound:_{max}, which can also be fitted to the data. The value of the fitted SPA_{max} should remain close to 1.0.

### Sodium binding by WT NIS SPA data

NIS is the key plasma membrane protein that mediates active I^{−} transport in the thyroid gland, the first step in the biosynthesis of the thyroid hormones, of which iodine is an essential constituent (Dai et al., 1996). NIS mutations have been linked to iodide transport defects, which lead to congenital hypothyroidism, deficit in the central nervous system development, and mental retardation if patients are not treated soon after birth (Dohán et al., 2002; De la Vieja et al., 2004, 2005; Paroder-Belenitsky et al., 2011; Li et al., 2013; Paroder et al., 2013; Portulano et al., 2014). NIS couples the inward “uphill” translocation of I^{−} against its electrochemical gradient to the inward “downhill” transport of Na^{+} down its electrochemical gradient. NIS activity is electrogenic, with a 2:1 Na^{+}/I^{−} transport stoichiometry (Eskandari et al., 1997). It is able to accumulate iodide from serum concentrations orders of magnitude below its *K _{d}*. We have recently shown that at physiological Na

^{+}concentrations, ∼79% of all NIS molecules have two Na

^{+}ions bound, and are therefore ready to bind and transport I

^{−}(Nicola et al., 2014).

SPA experiments were performed to determine the affinity of Na^{+} binding by NIS. NIS, like many other transporters in its family, has two independent binding sites for Na^{+}; both Na^{+} ions that bind to NIS are used to transport I^{−}. The *f _{x}* values, obtained and normalized as detailed in Materials and methods, were fitted to Eq. 8 as a function of [Na

^{+}] to obtain the

*K*and

_{p}*K*values (Fig. 2 A). The agreement between the calculated values (Fig. 2 A, continuous curve) and the experimental data is excellent

_{x}*f*values in the range of 0–1.0;

_{x}*m*is the number of experimental data points, and

*p*is the number of parameters.

The extremes of the ranges provide examples of the combinations of affinities compatible with the data (Fig. 3, top and bottom). If *K _{d,A}* =

*K*= 225 mM (

_{d,B}*K*= 1/

_{d,X}*K*) (Fig. 3, top), binding to a second site becomes highly favorable, with

_{a,X}*K*/ϕ = 4.5 mM (Fig. 3, bottom), a 50-fold enhancement. At the other extreme, when

_{d,B}*K*approaches 112.5 mM (i.e.,

_{d,A}*K*approaches infinity),

_{d,B}*K*/ϕ becomes 9.2 mM, corresponding to a very large enhancement. The actual values for NIS must fall within this range, somewhere between these two extremes.

_{d,B}This analysis shows that the values of the dissociation constants for empty NIS, *K _{d,A}*, and

*K*(Fig. 3, top, left side) are compatible at one extreme, with both sites having the same affinity (

_{d,B}*K*=

_{d,A}*K*= 225 mM), and as the affinity of one of the sites (for example, site A) increases, that of the other must decrease to remain compatible with the experimental data. As

_{d,B}*K*decreases, the value of

_{d,A}*K*increases dramatically, approaching infinite as

_{d,B}*K*approaches 112.5 mM (= 1/

_{d,A}*K*). However, because of the cooperativity,

_{p}*K*/ϕ becomes highly favorable when site A is occupied (range of 4.6–9.2 mM; Fig. 3, bottom, left side). These values are very similar to those obtained by a statistical thermodynamic analysis of kinetic data (

_{d,B}*K*= 137 mM and

_{d,A}*K*= 22.4 mM) (Nicola et al., 2014).

_{d,B}It should be noted that the dissociation constant for the first site (*K _{d,A}*) and the enhanced dissociation constant for the second site (

*K*/ϕ) are always known within a factor of 2. Furthermore, if additional information becomes available—for example, from binding data for affinity mutants—the values of all three constants (

_{d,B}*K*,

_{d,A}*K*, and

_{d,B}*K*/ϕ) can be obtained.

_{d,B}The same data fitted to Eq. 1 (Fig. 2 B) gave a *K _{d}* value of 31.9 mM and an

*n*of 1.74. This

*K*is highly similar to the geometric mean of

_{d}*K*and

_{d,A}*K*/ϕ (= 32.1 mM). Although the two equations yield fits of similar quality (see below), there is much more information in the values in Fig. 3 than in the two numbers obtained from fitting the data to the Hill equation.

_{d,B}### Sodium binding by S353A/T354A NIS SPA data

The S353A/T354A NIS mutant was studied as a possible example of a NIS molecule that may bind a single Na^{+} ion. Residues S353 and T354 were studied because the NIS mutant T354P causes iodide transport defect (Levy et al., 1998), and were identified as residues that are part of the coordination of one of the Na^{+} ions (De la Vieja et al., 2007), like the corresponding residues in the bacterial leucine transporter LeuT (Yamashita et al., 2005). Alanine substitutions at NIS positions 353 and 354 abolish I^{−} transport (De la Vieja et al., 2007). To determine whether or not the lack of transport activity is caused by a lack of Na^{+} binding, we generated the double mutant S353A/T354A and measured, using SPA, the binding of Na^{+} to this NIS mutant. Surprisingly, the mutant binds two Na^{+} ions, but the binding curve (Fig. 4 A) fails to show the positive cooperativity between the sites observed for the WT (Fig. 2 A).

The SPA data were fitted and analyzed using three different formalisms: (1) fit to Eq. 8 followed by the analysis proposed above (Figs. 4 A and 5); (2) fit to Hill equation (Eq. 1) (Fig. 4 B); and (3) fit assuming two independent sites (Fig. 4 C). The three formalisms fit the experimental data extremely well. Interestingly, the fit to Eq. 1 yields a single *K _{d}* = 20.6 mM and a Hill coefficient of 0.42. There is little information in these values besides a general

*K*and the indication of negative cooperativity (

_{d,50}*n*≤ 1). Furthermore, the fact that the value of

*n*is smaller than the minimum value (1/N = 0.5) expected for negative cooperativity (Abeliovich, 2005) suggests that Eq. 1 is not an appropriate model for ligand binding to this NIS mutant.

Analysis using the equations presented above yields a range of *K _{d}* values for the first site (

*K*) between 0.99 and 1.98 mM (Fig. 5, top). For most of this range, the dissociation constant for binding Na

_{d,A}^{+}to site B when site A is not occupied remains small (

*K*= 4.75 mM at

_{d,B}*K*= 0.88 mM

_{a,A}^{−1};

*K*= 6.23 mM at

_{d,B}*K*= 0.85 mM

_{a,A}^{−1}; Fig. 5), but it increases rapidly as

*K*approaches

_{a,A}*K*. Negative cooperativity is manifested in binding to site B when site A is occupied (

_{p}*K*/ϕ; Fig. 5, bottom, cyan). In this situation,

_{d,B}*K*/ϕ, the dissociation constant for site B when site A is occupied, varies between 43.84 and 87.68 mM (Fig. 5, bottom, cyan). The values on the right side of the graph, 0.99 and 87.68 mM, are highly similar to those obtained by fitting two independent sites (0.90 and 80.48 mM). The similarity between these sets of values reflects the fact that, under these conditions, the two models describe similar situations: in both cases, site B starts binding after site A is substantially occupied. The analysis based on the formalism introduced here, however, describes other situations that are equally compatible with the data from SPA experiments. As mentioned before, analysis with the Hill equation provides much less information.

_{d,B}### Equivalence of fitting by the two equations

It could be argued that the similarity between the quality of the fit using the Hill equation with a non-integer coefficient and that using the equations proposed here (i.e., similar rms_{resid} values) is a coincidence because of the particular values of the constants and *n* in the case of NIS. To test this possibility, we computed points with the Hill equation (Eq. 1) using values for *n* ranging from 1.7 to 1.9 and a *K _{d}^{n}* value of 0.001. The points in each individual curve were fitted by least squares using Eq. 8 (Fig. 6). The equation fits the calculated points with an rms

_{resid}of <0.004, indicating that the two equations can fit the same data equally well.

### Analysis of data obtained by other experimental methods

The formalism presented here can be used to analyze data obtained by other methods, for example, fluorescence (Reyes et al., 2013). In this case, *f* versus [L] data can be fitted by least squares to a version of Eq. 8 modified to take into account the fact that the signal with one L bound may not be half the signal when two L’s are bound. The modified equation has the form_{1} is the fraction of the signal that occurs in response to binding of a single L (*sr* in Eq. 6). The value of S_{1} can be estimated as part of the least-squares fit. The results obtained can be analyzed as described above. As mentioned earlier, these equations can be used when [L] ≈ L_{added}. In the case of fluorescence, this condition is met when [P] << *K _{d}*.

### Summary of the approach

The equations for binding of a ligand to an allosteric macromolecule with two sites can be derived using three parameters, *K _{a,A}, K_{a,B}*, and ϕ (the enhancement of the affinity for binding the ligand to one of the sites that results from binding the ligand to the other site) (Hines et al., 2014), where

*K*=

_{a,A−B}*K*· ϕ and

_{a,B}*K*=

_{a,B−A}*K*· ϕ (Fig. 1). Fitting the experimental data to the three parameters in Eqs. 2a and 2b results in many combinations of parameters compatible with the data (Hines et al., 2014). The simple variable change shown in Eqs. 3 (

_{a,A}*K*=

_{p}*K*+

_{a,A}*K*) and 4 (

_{a,B}*K*=

_{x}*K*·

_{a,A}*K*· ϕ) results in equations with only two parameters (

_{a,B}*K*and

_{p}*K*; Eqs. 5a and 5b) that fit the data much more stably. Once the values of

_{x}*K*and

_{p}*K*are obtained by least squares or other methods,

_{x}*K*and

_{a,A}*K*can be obtained as the two roots of a quadratic equation (Eq. 7a). Choosing A (arbitrarily) as the site with the higher affinity (the plus sign in Eq. 7a), the solutions are bound by the conditions

_{a,B}*K*approaches zero,

_{a,B}*K*,

_{a,A}*K*, and

_{a,B}*K*· ϕ (

_{a,B}*K*,

_{d,A}*K*, and

_{d,B}*K*/ϕ), in which

_{d,B}*K*and

_{d,A}*K*/ϕ are known within a factor of 2. This information, along with the fact that the values correspond to a well-defined binding reaction, makes the values obtained by this procedure ideally suited for establishing correlations with structural and computational results. Because

_{d,B}*K*and

_{a,A}*K*· ϕ are known within a factor of 2, the free energies of binding ΔG

_{a,B}_{a,A}and ΔG

_{a,A−B}, calculated as −RT

*ln*K, are known within ∼0.4 Kcal/mol (= RT

*ln*2). Additional information can be obtained from the values of ϕ. The free energy of the interaction between the two sites (Δg

_{A,B}) can be calculated as −RT

*ln*(ϕ).

### Conclusions

In summary, we present a method for analyzing binding data from a protein that binds the same ligand at two sites. In the case of a symmetric homodimer, the analysis provides values for the two binding constants: that for binding to the empty protein at either of the two sites (*K _{d,A}* and

*K*), and that for binding to the second site when the first site is occupied (

_{d,B}*K*/ϕ). In the case in which the two sites on the protein are different, the approach provides bounded pairs of possible values for the same constants. This information is missed when using the Hill equation with a non-integer coefficient. We tested the usefulness of the equations we propose by analyzing SPA data for NIS. In the case of WT NIS, the Hill equation fit gives a

_{d,B}*K*of 31.8 mM and an

_{d}*n*of 1.74 (Fig. 2 B). The approach presented here indicates that Na

^{+}binds to the first site with an affinity between 112.5 and 225 mM, and to the second site—once the first site is occupied—with an affinity between 9.2 and 4.6 mM. Although the square root of the product of the constants (32.2 mM) is highly similar to the Hill equation

*K*, significant additional information is provided by the values obtained using the new approach. Thus, as opposed to the Hill equation, which provides

_{d}*n*, an indication of cooperativity, and an overall

*K*-like constant, the approach presented here provides narrowly bounded ranges of pairs of the values of all

_{d}*K*s necessary to fully describe binding of the ligand. This formalism can also be used to analyze data of proteins that bind the same substrate at two binding sites.

_{d}## Acknowledgments

We thank Drs. D. Leahy and A. Lau for critical reading of the manuscript and insightful discussions.

This study was supported by National Institutes of Health (grant DK-41544 to N. Carrasco).

The authors declare no competing financial interests.

Eduardo Rios served as editor.

## Footnotes

- Abbreviations used in this paper:
- NIS
- Na
^{+}/I^{−}symporter - SPA
- scintillation proximity assay

- Submitted: 15 January 2015
- Accepted: 17 April 2015

This article is distributed under the terms of an Attribution–Noncommercial–Share Alike–No Mirror Sites license for the first six months after the publication date (see http://www.rupress.org/terms). After six months it is available under a Creative Commons License (Attribution–Noncommercial–Share Alike 3.0 Unported license, as described at http://creativecommons.org/licenses/by-nc-sa/3.0/).