-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathantiporter.html
283 lines (239 loc) · 15.5 KB
/
antiporter.html
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
<!DOCTYPE html>
<html>
<head>
<title>Steady State</title>
<style type="text/css">
body {
background-color: #FFFFFF;
font-family: Verdana, sans-serif;
font-size: 12 px
}
</style>
<script type="text/x-mathjax-config">
MathJax.Hub.Config({tex2jax: {inlineMath: [['$','$'], ['\\(','\\)']]}});
</script>
<script type="text/x-mathjax-config">
MathJax.Hub.Config({
TeX: {
equationNumbers: { autoNumber: "AMS" },
extensions: ["AMSmath.js", "AMSsymbols.js"]
}
});
</script>
<script type="text/javascript"
src="https://c328740.ssl.cf1.rackcdn.com/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML">
</script>
<script type="text/javascript" src="functions.js"></script>
</head>
<body>
<div style="display:none">
$
\newcommand{\conc}[1]{[\mathrm{#1}]}
\newcommand{\conceq}[1]{[\mathrm{#1}]^{\mathrm{eq}}}
\newcommand{\fidl}{F^{\mathrm{idl}}}
\newcommand{\kcat}{k_{\mathrm{cat}}}
\newcommand{\kmmon}{\kon^{\mathrm{ES}}}
\newcommand{\kmmoff}{\koff^{\mathrm{ES}}}
\newcommand{\kconf}{k_{\mathrm{conf}}}
\newcommand{\koff}{k_{\mathrm{off}}}
\newcommand{\kon}{k_{\mathrm{on}}}
\newcommand{\nai}{N_A^i}
\newcommand{\nao}{N_A^o}
\newcommand{\nbi}{N_B^i}
\newcommand{\nbo}{N_B^o}
\newcommand{\ss}{\mathrm{SS}}
$
</div>
<h2>
Active transport in an alternating-access antiporter model
</h2>
<a href="javascript:changeTo('antiporter','free')">Free energy</a> stored in a non-equilibrium gradient of concentrations across a membrane can be used by an antiporter ("exchanger") to pump another molecule against its gradient.
<ul>
<li> Membrane-embedded transporter proteins can function as machines because binding events are coupled to conformational changes. </li>
<li> The cycle below is a simplified, schematic model of a transporter for which binding events are coupled to "eversion" - transition between inward-facing and outward-facing conformations. </li>
<ul>
<li> The alpha helices in membrane proteins achieve eversion in roughly the same manner as you might re-adjust the orientations of a group of irregular pens and pencils held in your hand. </li>
</ul>
</ul>
<p style="text-align:center"><img src="images/antiporter.gif" alt="Antiporter cycle" /></p>
<h3>
Basics of the antiporter cycle
</h3>
<ul>
<li> The antiporter schematized above can perform active transport using the gradient of one molecule or ion (e.g., A) to pump the other (B) in the opposite direction.
<a href="javascript:changeTo('antiporter','free')">Free energy</a> is supplied by the gradient(s) -- <em> no other source of free energy, such as ATP, is needed.</em>
</li>
<li> If the ratio of outside to inside concentrations of one molecule (e.g., A) is
higher than the other, the
<a href="javascript:changeTo('antiporter','cycles')">cycle</a> will be driven in a
particular direction (e.g., clockwise). </li>
<li> A transporter such as the antiporter shown can be considered a "passive element": it does not supply any energy itself, but only transduces free energy stored in gradients. The transduction occurs due to the coordination of binding and conformational changes.
<li> Because it is a passive element, the antiporter does not "know" which molecule's gradient is meant to drive the transport of the other. Thus, if the outside/inside concentration ratio is greater for B, then the cycle will run counter-clockwise and A will be pumped by B.
<li> The model as shown is a simplified version that omits some possible
states and connections among states - e.g., both binding sites are assumed
to always face the same direction.
<a href="javascript:changeTo('antiporter','cycleLogicAdvanced')">Additional states and connections</a> could reduce the efficiency of pumping. See
also discussion in Hill's book.
The presence of additional states or connections would have to be verified on a
system-by-system basis. </li>
</ul>
<h3>
Understanding how the cycle is driven
</h3>
<ul>
<li> A simple but powerful qualitative analysis starts by considering the condition of <a href="javascript:changeTo('atpTransport','equil')">equilibrium</a>.
Equilibrium describes a (hypothetical) very large set or ensemble of identical systems of which - on average - an equal number are executing the forward and reverse of every process.
Thus, for example, between steps 1 and 2, if there are $N_{12}$ systems in which extracellular ("outside") A binds every second, there is an equal number ($N_{21} = N_{12}$) for which A unbinds to the outside.
</li>
<li> In the perfectly balanced state of equilibrium, there will be equal numbers of systems performing full clockwise and full counter-clockwise cycles per unit time.
Thus, although individual cycles may pump A and B, there is no net pumping.
</li>
<li> Because the balance of equilibrium is so perfect, however, it can be disturbed at almost any point in the cycle.
For example, adding an excess of A outside will cause more 1-to-2 transitions, which will cause 2-to-3 transitions and so on, leading to clockwise cycling.
It is the binding process that provides directionality.
<em>Binding events can be considered the "handles" used to drive a cycle in a given direction.</em>
</li>
<li> Importantly, the cycle can also be driven by imbalancing the detailed equilibrium at <em>any</em> point in the cycle:</li>
<ul>
<li> Excess B added inside will also drive the cycle clockwise.</li>
<li> Excess A added inside will drive the cycle <em>counter-</em>clockwise.</li>
<li> Excess B will also drive the cycle counter-clockwise.</li>
</ul>
</ul>
</ul>
<h3>
A quantitative model: Simulation and analysis
</h3>
<ul>
<li> We will employ a chemical-kinetics model, which is formulated solely in terms of state populations and rate constants for transitions among states, assumed to obey <a href="javascript:changeTo('atpTransport','massAction')">mass-action</a> behavior.</li>
<li> To keep the model simple, we will assume all rates are the same:
<ul>
<li> $\kon$ will be the on-rate (see <a href="javascript:notation()">Notation</a>) for all A and B binding -- for transitions from state 1 to 2, from 1 to 6, from 3 to 4, and from 5 to 4.</li>
<li> $\koff$ will be the off-rate for all unbinding of A or B -- reversals of the transitions noted above.
<li> $\kconf$ will be the rate for all conformational transitions -- both directions between states 2 and 3, and between 5 and 6.
<li> We will assume that the outside and inside volumes are the same -- so that the numbers of A and B molecules tell us the concentrations in a simple way.
</ul>
<li> The model is now fully specified. For example, the differential equation governing the population of state 1 (denoted [1]) is gievn by
\begin{equation*}
\frac{ d \conc{1} }{ dt } =
- \conc{1}\conc{A}_{\mathrm{out}} \, \kon - \conc{1}\conc{B}_{\mathrm{out}} \, \kon
+ \conc{2} \, \koff + \conc{6} \, \koff
\end{equation*}
<ul>
<li> In words, the equation means that the population of the unbound state 1 decreases due to binding of A or B (from the outside) with rate $\kon$ and increases due to unbinding from either of states 2 or 6 with rate $\koff$.
</ul>
</li>
</ul>
<h3>
Simulating the model
</h3>
If we start a simulation with all A molecules on the outside and an equal number of B molecules solely on the inside, the system equilibrates to equal concentrations on both sides -- for both species (A and B).
<p style="text-align:center"><img src="images/antiporter-equal-A-B-justA.png" alt="Simulation data" height="225" /></p>
<p style="text-align:center"><img src="images/antiporter-equal-A-B-justB.png" alt="Simulation data" height="225" /></p>
However, if we start a simulation with B molecules equally distributed between outside and inside, and the A molecule all outside as before, there is a very different outcome.
The greater outside concentration of A tends to drive A molecules inward, but doing so requires "pumping" B molecules to the outside.
Because the B molecules "resist" having different outside and inside concentrations, the final equilibrium is one where the <emph>ratio</emph> of outside/inside concentrations equalizes between the two species.
<p style="text-align:center"><img src="images/antiporter-unequal-A-B.png" alt="Simulation data" height="225" /></p>
These simulations were performed using <a href="http://bionetgen.org/" target="_blank">BioNetGen</a>, a rule-based platform for kinetic modeling.
The brief source code for the model (a .bngl file) can be downloaded by right-clicking <a href="pubFiles/antiporter.bngl">here</a>.
<h3>
Analyzing the model
</h3>
Perhaps the simplest and most important analysis of the model is deriving and understanding the <a href="javascript:changeTo('antiporter','equil')">equilibrium</a> behavior.
In equilibrium, every process and its reverse will be in exact (detailed) balance -- the flow in both directions will match.
Hence, we need to write down balance equations for each set of arrows (connectings states 1 and 2, 2 and 3, ...).
Using $\conceq{X}$ to refer to the equilibrium concentration of component X, we write down the equilibrium equations for the three processes on the right side of the cycle.
\begin{equation}
\conceq{1} \, \conceq{A}_{\mathrm{out}} \, \kon
= \conceq{2} \, \koff
\end{equation}
\begin{equation}
\conceq{2} \, \kconf = \conceq{3} \, \kconf
\end{equation}
\begin{equation}
\conceq{3} \, \conceq{B}_{\mathrm{in}} \, \kon
= \conceq{4} \, \koff
\end{equation}
<p>
The equations for the left side mirror these.
</p>
<p>
Solving the full set of equations (by going around the cycle and eliminating one of the numbered-state concentrations at a time) yields the result
</p>
\begin{equation}
\label{equiltot}
\frac{ \conceq{A}_{\mathrm{out}} }{ \conceq{A}_{\mathrm{in}} }
= \frac{ \conceq{B}_{\mathrm{out}} }{ \conceq{B}_{\mathrm{in}} }
\end{equation}
<p>
As we saw in the simulations, equilibrium is reached when the outside-to-inside concentration ratios match for A and B.
In fact, because our transporter is a passive element, this result is independent of the rates chosen -- even if the rates are not symmetric (but be careful that the set of rates in a <a href="javascript:changeTo('antiporter','cycles')">cycle</a> has to be physically consistent).
</p>
Note that, perhaps surprisingly, the equilibrium point will depend on the initial condition.
This is true because our simple antiporter model does not allow "slippage" (where only an A or a B is transported in some cycles).
If there were some mechanism for A and/or B to switch compartments independently of the other, then the equilibrium point would not only have equal concentration ratios, but the ratios would be 1 -- equal inside and outside concentrations.
<a href="javascript:changeTo('antiporter','cycleLogicAdvanced')">A more advanced discussion of the antiporter cycle</a> describes further possibilities and the relation to a symporter cycle.
<h3>
Thermodynamic analysis
</h3>
In contrast to the kinetic analysis above, which includes a description of specific states of the antiporter, we can perform a thermodynamic analysis of the same system.
Using an <a href="javascript:changeTo('antiporter','idealGas')">ideal gas</a> description, in analogy to what was done for <a href="javascript:changeTo('antiporter','binding')">binding</a> and for understanding <a href="javascript:changeTo('antiporter','freeGradient')">free energy storage in a gradient</a>, arguably leads to a simpler description than we encountered using kinetics.
The thermodynamic analysis will directly lead to the equilibrium condition \eqref{equiltot}.
<p style="text-align:center"><img src="images/antiporter-compartments.gif" height="275" /></p>
The thermodynamic model consistent with the mass-action picture above consists of ideal gases for each of the components: A molecules inside ($i$) and out ($o$), as well as B molecules inside and out.
The combined free energy $F$ for the whole system is the sum of the free energies for each of these ideal gases.
\begin{align}
F = & \fidl(\nai, V_i) + \fidl(\nbi, V_i) \nonumber \\
& + \; \fidl(\nao, V_o) + \fidl(\nbo, V_o) ,
\label{ftot}
\end{align}
where we have omitted the temperaure dependence that will not be pertinent here.
The total free energy is the simple sum of the component free energies because there are no (energetic) interactions between the components.
In an ideal gas, by definition, there are no interactions among molecules.
<p>
To obtain concrete results, we must specify one specific state of our system - something like an initial condition.
After all, the 1:1 exchange of A and B means that not every possible set of populations ($\nai, \nao, \nbi, \nbo$) is possible.
For simplicity, we shall assume that our reference state has all A molecules inside ($\nai = N_A, \nao = 0$) and all B molecules outside ($\nbo = N_B, \nbi = 0$).
</p>
Further, because we are considering a 1:1 antiporter (where one A always exchanges for one B), there is really only one variable that can change.
To see this, we can express all the species numbers in terms of $\nao$, the number of A molecules outside, and the <i>constant</i> species totals $N_A$ and $N_B$:
<ul>
<li> because an A molecule must be either inside or out, we have $\nai = N_A - \nao$; </li>
<li> because of the 1:1 exchange, $\nbi = \nao$; and </li>
<li> because a B molecule must be either inside or out, $\nbo = N_B - \nbi = N_B - \nao$. </li>
</ul>
Using these relations and noting that $\fidl(N,V,T) = N k_B T \ln [ N / (V/\lambda^3)]$, we can now re-write the total free energy from \eqref{ftot} as
\begin{align}
\frac{F}{k_B T} =
& \left( N_A - \nao \right) \ln \left( \frac{ N_A - \nao }{V_i \left/ \lambda_a^3 \right.} \right) \nonumber \\
&+ \; \nao \ln \left( \frac{ \nao }{V_o \left/ \lambda_a^3 \right.} \right)
+ \nao \ln \left( \frac{ \nao }{V_i \left/ \lambda_b^3 \right.} \right) \nonumber \\
&+ \; \left( N_B - \nao \right) \ln \left( \frac{ N_B - \nao }{V_o \left/ \lambda_b^3 \right.} \right)
\end{align}
We find the most probable state (which is equilibrium in the large $N$ limit) by setting to zero the derivative of $F$ with respect to $\nao$.
After some algebra, the result is
\begin{equation}
\label{equilthermo}
\frac{ \left. \nao \right/ V_o }{ \left. \nai \right/ V_i }
= \frac{ \left. \nbo \right/ V_o }{ \left. \nbi \right/ V_i } ,
\end{equation}
which is equivalent to \eqref{equiltot} because $N_X/V$ is the concentration of species X.
<p>
Importantly, <i>the thermodynamic analysis does not require knowledge of the system's substates, but only the 1:1 stoichiometry</i>.
</p>
<h3>
More on antiporters
</h3>
The stoichiometry of transport has very interesting effects: see the discussion of <a href="javascript:changeTo('antiporter','transportStoich')">2:1 exchange</a>.
Also, all the analysis above assumed perfect 1:1 exchange.
In reality <a href="javascript:changeTo('antiporter','cycleLogicAdvanced')">slippage</a> could be important, and is a critial general phenomenon in complex molecular machines.
<h3>
Exercises
</h3>
<ol>
<li> Derive \eqref{equiltot}. </li>
<li> Derive \eqref{equiltot} using arbitrary rates $\{ k_{12}, k_{21}, k_{23}, \ldots \}$. In doing so, you will want to be sure your <a href="javascript:changeTo('antiporter','cycles')">cycle is "thermodynamically consistent" with equilibrium principles</a>. In other words, you will need to determine a constraint, which can be done by assuming the system starts at the state of <a href="javascript:changeTo('antiporter','freeGradient')">minimum global free energy in which each species has equal concentrations on both sides</a>. </li>
<li> Derive \eqref{equilthermo}. </li>
</ol>
</body>
</html>