Journal of Biophysics

Volume 2013 (2013), Article ID 683256, 14 pages

http://dx.doi.org/10.1155/2013/683256

## Transport Reversal during Heteroexchange: A Kinetic Study

^{1}Department of Physics, UPR, San Juan, PR 00931, USA^{2}Department of Biochemistry, UCC, Bayamon, PR 00960, USA^{3}Department of Physiology, UCC, Bayamon, PR 00960, USA

Received 21 May 2013; Revised 18 August 2013; Accepted 18 August 2013

Academic Editor: Kuo-Chen Chou

Copyright © 2013 V. Makarov et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

It is known that secondary transporters, which utilize transmembrane ionic gradients to drive their substrates up a concentration gradient, can reverse the uptake and instead release their substrates. Unfortunately, the *Michaelis-Menten* kinetic scheme, which is popular in transporter studies, does not include transporter reversal, and it completely neglects the possibility of equilibrium between the substrate concentrations on both sides of the membrane. We have developed a complex two-substrate kinetic model that includes transport reversal. This model allows us to construct analytical formulas allowing the calculation of a “heteroexchange” and “transacceleration” using standard Michaelis coefficients for respective substrates. This approach can help to understand how glial and other cells accumulate substrates without synthesis and are able to release such substrates and gliotransmitters.

#### 1. Introduction

Unlike “primary” or ATP dependent transporters that create the major ionic gradients of K/Na/H and Cl/CO_{2} ions across cellular membranes harnessing the energy reserved in ATP, the “secondary transporters” utilize the energy available from transmembrane ionic and/or pH gradients and membrane potential to drive their substrates up a steep concentration gradient. Transporters on neurons and astrocytes clearing neurotransmitters from the synaptic cleft and extracellular space mainly belong to different “secondary transporters” families. Recently, it has been shown that astrocytes and other glial cells accumulate monoamines [1] and polyamines [2, 3] while lacking the enzymes for their synthesis [1, 4–6]. One among many known representatives of the “secondary transporters” that utilize the transmembrane ionic gradients and membrane potential is the family of organic cation transporters (OCT). These transporters take up different mono- and polyamines [7], and cells expressing such transporters also release these substrates using possibly two pathways: (i) large pores and (ii) transport reversal. Here we analyze one of transport reversal mechanisms.

*Energy Calculations*. Experimentally, it has been shown that secondary transporters can reverse their uptake releasing their substrates instead [8–10]. Energy based calculations were introduced to analyze the conditions for substrate release or uptake for this kind of transporter [11, 12]. It was established that substrate transport depends on the energy balance of coupled transport of the substrate and simultaneously transported ions (see Appendix A). Most secondary transporters could be reversed by membrane potential and by changes in the principal ion gradients and substrate concentrations. Experimentally, the reversal was shown for the glutamate transporters (for the review see [13]), GABA (reviewed by [12]), and for glial organic cation transporters [14, 15]. Being reversed, electrogenic transporters usually change the direction of the net ion flow. We summarize the energy balance study, introduced by Rudnick [11] in Appendix A. This analysis only studies one substrate uptake/release by a secondary transporter.

*Michaelis-Menten Scheme*. The kinetic concept based on the Michaelis-Menten scheme proved very useful for transporter mediated substrate uptake and inhibition [16, 17]. This kinetic model predicts saturability and specificity of secondary transporters in many cases, and atypical transport kinetics can be explained by multiple binding sites [18]. We have summarized this classic concept in Appendix B. Unfortunately, as one can see, the Michaelis-Menten model does not include transporter reversal, and it completely neglects the reversal constant (see Appendix B). A more complex transporter kinetic model is needed to predict quantitatively at least the following well-established experimental observations.(1)It has been shown that one transporter substrate can release another one already accumulated inside the cell. Sometimes this is called “heteroexchange.” For example, dopamine, tyramine, and amphetamine, which are substrates for the neuronal dopamine transporter (DAT), can release the substrate named N-methyl-4-phenylpyridinium (MPP) through DAT [19], with releasing ability of these substances correlated with the elicited coupled transport current. Also, it was shown that L-glutamate and its transportable analogs (substrates for EAATs) specifically release L-aspartate (another EAAT substrate) through this transporter and can be blocked by nontransportable analogs [20].(2)A special term was coined for the release of the (tracer) substrate by the same substrate, a process named “transacceleration.” While the phenomenon is not kinetically different from the “heteroexchange” described in the previous paragraph, it is well established experimentally (see, e.g., [21]). As new transporter models arise (e.g., a channel-transporter model [22]), it might be important to get this phenomenon explained by a purely thermodynamic model, not by using kinematic assumptions.

Here we present kinetic algorithms that more accurately explain the behavior of a secondary transporter pumping two substrates simultaneously; it predicts transporter reversal by the application of an additional second substrate to the transporter already in equilibrium with the first substrate, “transacceleration” and other interactions.

#### 2. Results

We modified the Michaelis-Menten kinetic model to include two different transportable substrates and also additional elementary steps, characterized by their kinetic coefficients, which are necessary for the transporter not only to uptake but also to release substrates. The model is presented in Appendix C by relations 1–8. This model can be considered as a system of kinetic equations describing the dynamics of the model –. A general solution for this scheme is difficult to obtain analytically. But some particularly interesting cases can be resolved (see Appendices C.1, C.2, and C.3), and we are presenting them below.

##### 2.1. Equilibrium Conditions for Both Substrates (See Appendix C.1)

Practically, the initial concentrations of substrates and are considered known (i.e., and ), and then substrate concentration can be measured in the outside solution ( and in our notation for this section of Appendix C). In this way, we tried to reduce all equations to measurable parameters.

It follows from relationships (Appendix C) that at fixed concentration of (initial concentration of first substrate ) and variable concentration of (different initial concentrations of ), the equilibrium concentration of (the substrate outside) increases with increasing (the effect of substrate *releasing* from the cell), and similarly, the equilibrium concentration of decreases with decreasing (effect of substrate transport inside of the cell). The same behavior follows from for the equilibrium concentration of at fixed concentration of and variable concentration of . These respective dependencies are shown in Figure 1.

*Conclusions of Appendix C.1*(1)Effect of substrate being *released* in case of competition in the two-substrate system can be observed, if at equilibrium condition most of the transporter is coupled by both substrates of interest: .(2)Efficiency of the substrate *releasing* process is dependent on equilibrium constant values describing processes of substrate-transporter intermediate complex formation.(3)Correct sign of the square root term in relations is defined by the conditions of
(4)Relationships can be used for analysis of equilibrium substrate concentration dependence on initial substrate concentrations.(5)The relationship of
can be used to determine parameters and , if concentrations of and can be simultaneously measured as functions of and , the initial concentrations of the first and second substrates.

##### 2.2. Two-Substrate System Dynamics at the Initial Time

Transporter velocities (transport rates) can be determined if (similar to Michaelis-Menten scheme) there is no equilibrium between substrate concentrations inside and outside and processes of the type can be neglected. We also assume there are the initial conditions where and are added to the external solution, thus and . In that case (see Appendix C.2), (1) are analogous to the Michaelis-Menten formulation for a two-substrate system. If , we obtain the exact Michaelis-Menten formula for the first substrate velocity, and if , we obtain the exact Michaelis-Menten formula for the second substrate velocity. Note also that the term can be interpreted as , and as .(2)The constants , , , and can be determined experimentally similar to the standard procedures used in the Michaelis formulation. There are some important equations:(i)if and : (ii)If and :

##### 2.3. Effect of the Equilibrium Reverse Bias for a First Substrate When a Second One is Added to the System

If previously the equilibrium was established for a first substrate between outside concentration of the substrate and the inside concentration, the addition of a second substrate will produce a reverse bias (equilibrium shift). In the beginning, at initial time, some of the transporter molecules in the outside bind to the second substrate while inside there is still no second substrate. That means the availability of outside transporter for a first substrate becomes reduced. Thus equilibrium for a first substrate starts to break down; that is, the velocity of first substrate transport to outside (release) becomes bigger than its transport to the inside. At initial times during the start of the process and far from equilibrium for a second transporter, (8) allows the calculation of the velocity of the first substrate release due to transport reversal (see Appendix C.3): where is Michaelis constant for a second substrate and is initial concentration of a second substrate, where Thus, finally we have Taking into consideration the relation (8) we have calculated the dependence of the velocity of first substrate release on the concentration of a second substrate, at initial times after it was added to the system. Functional dependence (8) is represented in Figure 2.

It can be seen from (8) and Figure 2 that the velocity of reversed transport (release) of a first substrate is 0 if , because there is equilibrium between the velocities of inward and outward flow of the first substrate through the transporter. Thus, the “net” velocity, is equal to zero. Also, from (8) and as seen in Figure 2, with increase of a second substrate concentration, when , the velocity of a first substrate release becomes saturated and can be calculated as Release velocity depends on the first substrate concentration , and at given value of the value of is a constant. Thus, (8) for the velocity at half maximal value at a certain concentration of second substrate can be written as and because of this equation it can be calculated as There is similarity between the formula of velocity of transporter reversal due to second substrate addition and Michaelis-like formulas for the velocity of substrate uptake.

The formula that predicts the velocity of substrate uptake (see Appendix C or Appendix B) can be written as where the maximum velocity is represented by Thus, for the half maximal velocity, Thus we can write To say in plain words, the Michaelis constant for a second substrate can be determined in two ways: (i) from the standard Michaelis formulas at transport velocity measurements for the second substrate, or (ii) from the release velocity measurements of a first substrate, from our formula, where a second substrate produces release of the first one.

In the most important case, if (8) then can be interpreted as the release force for a first substrate after the addition of a second one. In the case of , the release force can be written as that is, in this case the release force for a first substrate after the addition of a second one has a linear dependence on the first substrate concentration.

#### 3. Discussion and Conclusions

We have studied the extended kinetic model for a secondary transporter simultaneously dealing with two substrates, which includes direct (outside-in) and reverse transport (inside-out). The model was solved in different equilibrium conditions (see Appendices C.1, C.2, and C.3). We have shown that when both substrates are in equilibrium, addition of one of them leads to reequilibrium and release of the second substrate (Appendix C.1). This was emphasized in Appendix C.3, when the system was studied for conditions where a first substrate is in equilibrium (inside-outside concentrations) and a second one is just added and is far from equilibrium. This situation is of a special interest as it has been studied experimentally [19, 20]. Also, this is what probably happens when methamphetamine, ephedrine, or other similar substances induce dopamine (and other monoamine) release from monoamine neurons primarily via membrane transporters, reversing the dopamine transporter (DAT), norepinephrine transporter (NET), and/or serotonin transporter (SERT) [23–27] and also reversing VMAT vesicular transport [28]. In addition, it has been recently shown that astrocytes and other glial cells accumulate polyamines [2, 3] while lacking the enzymes for their synthesis [4–6], and OCT type of transporters (that are expressed in glia) take up different polyamines [7]. Polyamines are released in brain from glial cells, but the mechanisms of such release are unknown [29].

Actually, as we understand now from formula (8) it can be ANY transportable substrate. This formula allows us to classify experimental measurements of a “heteroexchange” related substrate release for substrate-transporter pairs, using standard Michaelis coefficients.

The special term for the release of the (tracer) substrate by the same substrate, a process named “transacceleration,” can be explained by changes in equilibrium according to formula (8). There is no fundamental thermodynamic difference if the system has two chemically distinct substrates for the same transporter or there are radiolabelled and unlabelled chemically similar substrates. Thus, a new added substrate produces the release of a similar tracer substrate (labelled, e.g., with radioactive isotope) by equilibrium shift as shown in Appendix C.3.

We also have shown that if we assume both substrates are far away from equilibrium, and transporter reversal can be neglected (Appendix C.2), the formulas for the uptake velocity of both substrates become the same as in the Michaelis-Menten scheme (see Appendix C.2, and ), with the respective inhibitory coefficients.

We suggest that formula (8) will be especially useful in the study of polyspecific transporters with known multiple substrates, such as the organic cation transporters (OCT) that participate in the transport of different monoamines [30], as well as polyamines [7].

#### Appendices

#### A. Classical Energy Based Calculations

Similar to the analysis of ion channels, substrate flux through transporters can be determined by the transmembrane electrochemical potential () which is the sum of the electrical potential () and chemical potential (). For a single molecule , the driving force is quantified as where = the valence of , = Faraday’s constant, = membrane potential, = universal gas constant, and = temperature. It should be noted that when (i.e., if is at equilibrium), this equation reduces to the Nernst equation.

For the glial glutamate transporter GLT1, for example, all of the cotransported ions are coupled to each other as they cross the membrane (i.e., they are not independent). Therefore, the total electrochemical driving force for GLT1 is the sum of the linked contributions from each co-transported ion. Because one thermodynamic reaction cycle for GLT-1 involves coupled translocation of three sodium ions, one proton and one negative glutamate molecule and countertransport of one potassium ion are quantified as
In equilibrium, when and knowing that the K^{+} term must be negative as it is going in the opposite direction, this equation reduces to
which defines the reversal potential for the transporter. It should be noted that this latter equation has a very similar form to the Goldman-Hodgkin-Katz equation for ion channels. The only difference is that ion fluxes are coupled unlike the fluxes for ion channels. The calculated driving force for a transporter can be viewed in the same way as the driving force for an ion channel; thus, there will be no net substrate flux when membrane potential is equal to the reversal potential [12, 13]. This last equation can be rearranged also like this, clearly showing the substrate gradient produced by the transporter:

#### B. Derivation of Michaelis-Menten Equation for the Transporter

Kinetic concept takes into account only steady-state velocities of transport that can be divided in adhesion, transport, and release of substrate on other side of the membrane. Let the transporter and the substrate first form the complex in the outer membrane, and then substrate is transported to the inner membrane and released. The reversal is not taken into consideration. This can be written as follows: Formation of complex is proportional to the concentration of substrate and the free transporter: when the complex transports the substrate and releases it inside the cell proportionally to the concentration of , and also some complex just releases substrate again without transporting it. Thus, the removal of from the system is After complex formation the quantity of free transporter left is At steady state (at equilibrium), formation and removal of became the same: As velocity of substrate transport is proportional to and , we can write As in abundance of the substrate all velocity depends only on transporter , maximal velocity of the transport () is and the constant , we can write The formula allows determination of from the experimental curve: , which means becomes equal to if () = 1, and it happens when .

So, when the speed of transport is saturated and becomes it is simple to find half of that is equal to . Michaelis constant are used as a measure of substrate affinity to the transporter. Please, note that does not include the constant, reflecting unidirectionality of transporters in the Michaelis-Menten approach.

*Competitive Inhibition from the Michaelis-Menten Point of View.* Transporters can have another substrate . It can bind to the transporter, whether it is transported or not
The constant of dissociation of this complex can be written as follows:
On the other hand, from the equation
if we replace from (B.7) with the last formula, we will get
The transporter can be in free form or can be occupied by inhibitor or by substrate; therefore, one can write
This means the velocity of the transport of the main substrate will be

#### C. Kinetic Scheme with Two-Substrate Uptake-Release by the Same Transporter

In a general case, the simplest representation of the two-substrate scheme can be determined by set of elementary steps, characterized by their kinetic coefficients:(1); —formation of the intermediate complex “first substrate-transporter” () outside the cell;(2)— () dissociation outside the cell;(3)—dissociation of () complex inside of cell, with substrate released inside;(4)—formation of the first substrate complex with transporter inside the cell;(5)—formation of the intermediate complex “second substrate-transporter” () outside the cell;(6)—dissociation of () complex outside of cell;(7)—dissociation of () complex inside of cell;(8)—formation of the second substrate complex with transporter inside the cell.

Let us introduce notations of Series of kinetics equations describing dynamics of the system of interest can be represented as follows: Summarizing (C.3), (C.4), and (C.7), we will obtain As well as we can write Therefore, the equation system represented above can be rewritten as follows:

##### C.1. Equilibrium Conditions for Both Substrates

In equilibrium conditions, we can write where are equilibrium constants of processes of interest. Taking into account relations (C.12) for equilibrium constants, we can form combinations of such constants, which may be represented as follows: Let us consider the case, when , - all transporter is coupled to complexes of () and (). In this case, we can write with a fair approximation the relations of From the last relationship, the concentration can be represented as follows: Applying the same approach to equations related to the second substrate and taking into account the last relationship for concentration of we can write If relationship of is satisfied, we can write and finally we can write where Let us consider an additional relationship, which can be formed by multiplication of all equilibrium constants, that is, by multiplication of relationships of In this case, we will obtain relationship of where with respect to relationship (C.18), we may assume that the relationships of are satisfied, and we may finally write Taking into account relationship of we can write and if we will substitute the last value to relationship (C.24b), we will obtain where Thus, Solution is physically reasonable, if the relation of is satisfied. Since parameter , it can easy be shown that condition (C.31) is satisfied at all reasonable values of and parameters. It also follows from relationships (C.30) that at fixed concentration of and variable concentration of , equilibrium concentration of increases with increase of (effect of substrate realizing from the cell) and decreasing of equilibrium concentration of with decreasing of (effect of substrate transport inside of the cell). The same behavior follows from relationships of interest for equilibrium concentration of at fixed concentration of and variable concentration of . The respective dependences are shown in Figure 1.

*Conclusions of Appendix C.1*(1)Effect of substrate realizing in case of transportation competition for two substrate system can be observable, if in equilibrium condition most of the transporter is coupled by both substrates of interest: .(2)Efficiency of substrate realizing process is dependent on equilibrium constant values describing processes of substrate-transporter intermediate complex formation.(3)Correct sign in front of square root in relations (C.30) is defined by the conditions of
(4)Relationships (C.30) can be used for analysis of equilibrium substrate concentration dependence on initial substrate concentrations.(5)The relationship of
can be used to determine parameters of and , if concentration of and can simultaneously be measured as function of and initial concentrations of first and second substrates.

*Procedure of Experimental Measurements of ** Relationship *(1)Measurements of equilibrium concentrations set of , versus and .(i)Let us choose a fixed set of concentrations from the interval []. For a fixed concentration , let us measure the dependence of concentrations and versus []. The same procedure must be repeated for another fixed concentration , the third fixed concentration , and so forth.(ii)Using the equation
it can be rewritten as
We can experimentally measure , , and , with unknown parameters: . Applying the least-squares method, we can write
Thus,
Main system determinants can be represented as follows:
Therefore,
Earlier we had determined the parameters of interest as
That is,
As , , , :
we can write

##### C.2. Two-Substrate System Dynamics at the Initial Time

In this case we will not take into account processes (4) and (8): In this case, the series of kinetics equations can be represented as follows: For , , we may use quasi-stationary approximation: Thus, Thus, For the initial time after the start of all processes, we can assume that and ; that is, Hence we can write where are Michaelis constants for the first and second substrate, respectively, and —is a factor of transporter inhibition.

*Conclusions for Substrate Dynamics Section*(1)Equations
are the analogs to the Michaelis formula for a two-substrate system. If , we get the Michaelis formula for the first substrate, and if , we get the formula for the second one.(2)Constants , , , can be determined experimentally:(iii):
(iv):

##### C.3. Effect of the Equilibrium Reverse Bias for a First Substrate When a Second One Is Added to the System

If previously the equilibrium was established for a first substrate between the outside concentration of the substrate and the inside concentration, the addition of a second substrate will produce a reverse bias (equilibrium shift). To study this, we first can write the equation representing the velocity of second substrate transport: where In the beginning, at initial time, some of the transporter molecules in the outside bind to the second substrate while inside there is still no second substrate. This means the availability of outside transporter for a first substrate becomes reduced. Thus, equilibrium for a first substrate starts to break down; that is, the velocity of first substrate transport to outside (release) becomes bigger than its transport to the inside. Our model allows evaluation of the release of a first substrate with simple approximation.

Unlike that in the Appendix C.1, where all the transport molecules were available for both substrates, let us approximate that the amount of transporter binding the first substrate outside is reduced by the second substrate also binding to it. Thus, using this rough approximation, the available outside transporter is reduced to the amount: where while inside (the cell) the amount of transporter available for the first substrate still remains . Consequently, the velocity of the first substrate reverse flow can be represented as where and are some effective parameters that can be only analytically solved in more complex model that will implicate intermediate transporter-substrate complexes and are out of the scope of this study. Here we can only get formulas for the initial times during the start of the process and far from equilibrium for a second transporter. In (C.61) we had balanced the flow disparity between the first substrate in-flow (uptake) and out-flow (release), where and are the equilibrium concentrations of the first substrate outside and inside (the cell). Thus, these concentrations can be approximately found using the approximation that still the system is close to the equilibrium for a first substrate: Thus, For initial time conditions with relatively good precision, we can accept that as where Thus First substrate reverse flow (the release speed) velocity dependence on the external second substrate concentration is represented on Figure 2.

One can see that the reverse flow of the first substrate is absent if