## Abstract

Voltage-gated ion channels are crucial for electrical activity and chemical signaling in a variety of cell types. Structure-activity studies involving electrophysiological characterization of mutants are widely used and allow us to quickly realize the energetic effects of a mutation by measuring macroscopic currents and fitting the observed voltage dependence of conductance to a Boltzmann equation. However, such an approach is somewhat limiting, principally because of the inherent assumption that the channel activation is a two-state process. In this analysis, we show that the area delineated by the gating charge displacement curve and its ordinate axis is related to the free energy of activation of a voltage-gated ion channel. We derive a parameter, the median voltage of charge transfer (V_{m}), which is proportional to this area, and prove that the chemical component of free energy change of a system can be obtained from the knowledge of V_{m} and the maximum number of charges transferred. Our method is not constrained by the number or connectivity of intermediate states and is applicable to instances in which the observed responses show a multiphasic behavior. We consider various models of ion channel gating with voltage-dependent steps, latent charge movement, inactivation, etc. and discuss the applicability of this approach in each case. Notably, our method estimates a net free energy change of approximately −14 kcal/mol associated with the full-scale activation of the Shaker potassium channel, in contrast to −2 to −3 kcal/mol estimated from a single Boltzmann fit. Our estimate of the net free energy change in the system is consistent with those derived from detailed kinetic models (Zagotta et al. 1994. *J. Gen. Physiol.* doi:10.1085/jgp.103.2.321). The median voltage method can reliably quantify the magnitude of free energy change associated with activation of a voltage-dependent system from macroscopic equilibrium measurements. This will be particularly useful in scanning mutagenesis experiments.

## INTRODUCTION

Voltage-gated ion channels are integral membrane proteins that, upon sensing a change in the membrane electric field, open a passage for ions to flux through the membrane. They are represented in all major kingdoms of life and are crucial for both electrical and chemical signaling pathways in higher organisms (Hille, 2001). Several inherited diseases such as arrhythmias and epilepsies have been shown to be correlated to mutations in these proteins, thereby underscoring their physiological importance (Lehmann-Horn and Jurkat-Rott, 1999). To understand the mechanisms of ion channel gating and function, it is necessary to obtain accurate estimates of energetic effects of site-specific mutations. One widely used approach involves measurements of macroscopic ionic currents for a series of voltage steps from which one can derive the relative fraction of open channel at various potentials (P_{O}-V). These responses show a sigmoidal voltage dependence and are typically characterized by fitting to a single Boltzmann equation. For such a curve, the chemical free energy difference between the open and the closed state is defined by *z*, the Boltzmann slope, and V_{1/2}, which is the voltage that elicits the half-maximal response. This free energy difference is also referred to as the free energy of channel opening at zero voltage.

However, the voltage-gated ion channels are known to transit through several intermediates to reach their final open states (Cole and Moore, 1960; Vandenberg and Bezanilla, 1991; Bezanilla et al., 1994; Zagotta et al., 1994a; Schoppa and Sigworth, 1998). The slope factor, which for a two-state process is the charge translocated, in the case of a multistate process depends not only on the charge displaced during activation but also on stabilities of the intermediate states (McCormack et al., 1991; Sigworth, 1994; Bezanilla, 2000; Villalba-Galea et al., 2008). Thus, mutations that affect only the energies of the intermediate states would generate erroneous values of the net free energy difference between the initial closed and the final open states when derived from Boltzmann fits.

Furthermore, when the central pore opens, the channel molecule has already undergone significant conformational transitions as indicated by separation in the charge-voltage (QV) and conductance-voltage (GV) curves (Armstrong and Bezanilla, 1973; Perozo et al., 1993). Therefore, the free energy associated with the activation process cannot be fully captured from conductance-voltage relationships. Single Boltzmann fits of GV curves for Shaker potassium ion channel give an estimate of −2 to −3 kcal/mol in the chemical free energy difference between the open and closed state (Monks et al., 1999; Li-Smerin et al., 2000a; DeCaen et al., 2008). This estimate, which corresponds to a single hydrogen bond or salt bridge interaction, is surprising and seems to suggest that a single interaction site could tip the balance between the two states.

Finally, calculating the free energy change using single Boltzmann fits becomes even more problematic when mutations result in discernibly multiphasic P_{O}-V curves (Li-Smerin et al., 2000b; Yifrach and MacKinnon, 2002). Boltzmann fits to the P_{O}-V curve also neglect the energetic effect of nonunity values of the maximum open probability, which has been observed in voltage-dependent ion channels (Schoppa et al., 1992). These observations prompted us to carefully analyze the free energy principles associated with a voltage-dependent process at a theoretical level. Our objective was to derive a more general, physically consistent measure of free energy change that is easily obtainable from experimental observables.

In this paper, we analyzed the fundamental basis of free energy change in a voltage-dependent system from the standpoint of classical statistical mechanics. The pioneering work of Wyman and Gill (1990) on theoretical thermodynamics has greatly furthered our understanding of ligand binding equilibria. Of specific interest to us was the concept of median ligand activity. The elegant derivation of Wyman (1964, 1967) shows that from the knowledge of the maximum number of ligands that bind a macromolecule and the median ligand activity, one can obtain an accurate estimate of the mean ligand binding affinity. Here, we adopted a similar approach to analyze the free energy relationships of voltage-dependent ion channels. We found that a comparable parameter, the median voltage of activation, as calculated from the QV curve, can be related to the net chemical free energy change associated with voltage-dependent activation of these proteins. Specifically, this procedure yielded a value of −14.6 kcal/mol associated with full-scale channel activation of the Shaker potassium channels and −16.1 kcal/mol for the voltage-gated sodium channels. We have discussed several characteristic models of ion channel activation such as those with inactivation, cooperativity, voltage-independent steps, latent charge movement between multiple open states, etc. and illustrate the validity and robustness of our proposed method.

## MATERIALS AND METHODS

### Oocyte expression

The cDNA voltage-gated Shaker potassium channel, with a deletion in the N-terminal region (residues 6–46) and with the W434F mutation, was linearized with NotI (New England Biolabs, Inc.) and in vitro transcribed to generate the mRNA using the mMESSAGE mMACHINE T7 kit (Life Technologies). Stage V *Xenopus laevis* oocytes were injected with 50.6 nl mRNA (at a concentration of 0.1 µg/µl). After injection, the oocytes were kept at 18°C in a solution containing 100 mM NaCl, 2 mM KCl, 1.8 mM CaCl_{2}, 1 mM MgCl_{2}, 5 mM Hepes, 0.1 mM DTT, and 0.2 mM EDTA, supplemented with 100 µg/ml gentamicin and 100 mg/ml bovine serum albumin. Measurements were performed 1–2 d after injection.

cDNAs of both the α and β subunit of the rNa_{V}1.4 were transcribed as described in the previous paragraph. Equimolar ratios of the α and β subunit mRNAs were coinjected into *Xenopus* oocytes to a final volume of 50 nl. Injected oocytes were preserved as in the previous paragraph, and measurements were performed 3–5 d after injection.

### Gating current measurements

The gating current measurements were performed on a cut-open oocyte voltage clamp set-up (CA-1B; Dagan Corporation) as described previously (Muroi et al., 2010; Lacroix and Bezanilla, 2011). For the potassium channel gating current measurement, the external solution was 115 mM NMG-MES (*N*-methyl-d-glucamine methanesulfonate), 2 mM Ca-MES, and 10 mM Hepes, pH 7.4. For the sodium channel gating current measurement, the external solution was 115 mM Na-MES, 2 mM Ca-MES, and 10 mM Hepes, pH 7.4. In the latter case, all ionic currents were blocked by the application of 10 µM tetrodotoxin to the external and middle chambers. For both channels, the internal solution was 115 mM NMG-MES, 2 mM EGTA, and 10 mM Hepes, pH 7.4. The recording pipette resistance was 0.3–0.5 MΩ. Analogue signals were sampled at 250 kHz with a Digidata 1440 interface (Molecular Devices) and low-pass filtered at 10 kHz. The capacitive transient currents were subtracted online using the P/4 method with a subtraction holding potential of −120 mV for the potassium channels and 50 mV for the sodium channels. Gating currents were obtained by applying a depolarizing pulse (50 ms for potassium and 20 ms for sodium channels) to voltages from −120 to 10 mV (at 5-mV intervals) for the potassium channels and −160 to 30 mV (at 10-mV intervals) for the sodium channels. The holding potential was −90 mV, and a 50-ms-long pre- and postpulse at −130 mV was used.

### Analysis and simulations

Each of the QV curves is a mean obtained from five oocytes. The curve was fitted to various equations as described in the Table 1. Computation of the area was performed using the trapezoid method. In brief, for a fractional gating charge displacement versus voltage (QV) curve, with *n* points, the area between the curve and the ordinate (Q) axis was calculated as* _{i}* is the

*i*

^{th}point on the QV curve and Q

*is the fraction of charge that is transferred at voltage V*

_{i}*. Numerical simulations of the different models were performed using MATLAB 2008b (MathWorks).*

_{i}## RESULTS

### Theory

#### Free energy change in a two-state process is obtained from a Boltzmann fit.

Consider a voltage-dependent ion channel that can exist in two states, closed (C) and open (O). On increasing the voltage, the open state of the channel becomes increasingly populated. The associated free energy difference between the closed and open states can be written as *q* is the gating charge translocated when the ion channel activates and is responsible for an electrical component in the net free energy change associated with activation. The equilibrium constant for such a transition at any voltage will be_{O} and P_{C} are the fractional occupancy of the open and closed states, respectively, and *K*^{0} is the voltage-independent component of the equilibrium constant, which is related to *K*^{0} as exp(−*qFV*_{1/2}), where V_{1/2} (Stevens, 1978)is the voltage at which half of the channels are open, P_{O} can be reexpressed as

By fitting an experimentally derived P_{O}-V curve to Eq. 1, one can derive the parameters of the Boltzmann fit, V_{1/2}, and *q* and thereby estimate *q*FV_{1/2}. A similar equation can be derived for the gating charge displacement versus V relation. If a channel can access only two states, the normalized QV and P_{O}-V curves will superimpose on each other. Also, the Boltzmann slope of the curves will equal the exact number of gating charges displaced during channel activation. However, voltage-gated ion channels are known to transit through a large number of intermediates during activation; the experimentally derived P_{O}-V and QV curves do not superpose, with the latter almost always preceding the P_{O}-V curve on the voltage axis. Also the Boltzmann slope of either curve is much lower than the actual amount of gating charge displaced during channel activation (Schoppa et al., 1992; Aggarwal and MacKinnon, 1996; Seoh et al., 1996). Thus, it is likely that free energy estimates of ion channel activation via a Boltzmann fit to activation response are inaccurate.

#### Area under the activation curve is a measure of the free energy change of activation.

Consider that the voltage-gated ion channel exists in *n* different conformational states, with each state *i* associated with a gating charge (or valence), given by*z _{j}* of the channel protein in state

*i*, which are informed principally by the charges on the side chains of specific ionizable residues of the protein. λ

*is a state-dependent parameter representing the fraction of the membrane electric field sensed by the*

_{ji}*j*

^{th}charge of the protein in the state

*i*(Stevens, 1978; Roux, 1997).

*q*, as described by Eq. 2, essentially governs the voltage dependence of the energy of a given conformational state, E

_{i}*, i.e.,*

_{i}The mean gating charge of the entire ensemble of conformational states at any given voltage is (Sigworth, 1994)* _{i}*(V) is the equilibrium occupancy of the

*i*

^{th}state at voltage V, which is assumed to follow a Boltzmann distribution. Thus,

*k*T and Z is the partition function of the system,

_{B}Changing the electrical potential (V) will alter the distribution of the channel states and thus the mean ensemble free energy. At a constant temperature, differentiating Eq. 5 with respect to voltage gives_{1} to V_{2} can be given by

The gating charge displacement measured when the voltage is switched from a reference voltage (V_{ref}) to a voltage V is essentially the difference _{ref} is usually a hyperpolarizing voltage when all protein charges are retracted to their initial resting configuration, and, without any loss of generality, _{1} as V_{ref} and using V instead of V_{2}, Eq. 8 can be integrated by parts, giving_{max}, which makes Eq. 9_{max}V) scales linearly with voltage, and the second component (

#### The median voltage of activation.

According to the definition of median ligand activity proposed by Wyman (1967), we define the median voltage of activation, V_{m}, as the voltage at which_{m} be the gating charge displacement at V_{m}.

By definition, _{m} (Eq. 11) and the fact that

Fig. 1 B shows the QV curve intersected by the V = V_{m} axis, with the two shaded regions being equal in area. The areas on the right and left of the V_{m} axis correspond to the integrals on the right and left side of Eq. 11, respectively. The area of the rectangle shaded in Fig. 1 C is equal to the sum of the two areas shaded in Fig. 1 A, taking the sign of the latter areas into consideration.

#### A physical interpretation of Δ G − C and V_{m}.

Thus far, we have shown that the area under an experimentally measured QV curve can be separated into two components, a linear and a saturating (converging) component, both being dependent on the total charge moved between the initial and final states of the system. The saturating component,

According to Eq. 3, the energy of each state of the system can be written as *i*^{th} state. The Boltzmann weight of each of these states can be expressed as *K _{i}* exp(

*q*β), where

_{i}V*K*= exp(−

_{i}*E*β). We assume that at sufficiently hyperpolarizing voltages, the system occupies just one (reference) state for which

_{i}*q*= 0 and

_{r}*i*. The partition function of the system is now

Now we turn to Eq. 11 and rewrite it as*k _{B}*T and rewriting Q

_{max}Vβ as ln exp(Q

_{max}Vβ), where (

*k*T = 1/β), Eq. 18 becomes

_{B}_{max}V

_{m}, is thus the measure of chemical free energy difference between the initial resting state of the system and the final activated state of the system, i.e., the free energy difference between the initial and final states of the system in the absence of any electric field (0 mV). Second, at V

_{m}, the Boltzmann weight of the fully activated state of the system,

_{m}, the state or states where all charge has moved (fully activated) and that where none has moved (fully resting state) are equally populated.

The relations established are independent of the pathway connecting the initial and final states of the system and the number of possible states of the system. Measurement of the QV can thus directly be related to the change in the free energy without elaborate model fitting procedures or assumptions about the nature of the conformational transitions (such as the two-state process assumption required by Boltzmann fits). The relations remain valid even when the QV curves show biphasic or asymmetric behavior. It requires the knowledge of Q_{max}, which, at least for the voltage-gated potassium (K_{V}) channels, is known with reasonable certainty (Schoppa et al., 1992; Aggarwal and MacKinnon, 1996; Seoh et al., 1996). Thus, at a first level of approximation, the measurement of the median voltage of activation would facilitate a straightforward and accurate calculation of the free energy changes associated with the voltage-dependent activation of any system.

### Accuracy of Q_{max}V_{m} as a measure of free energy change

To compare the free energy values derived from Q_{max}V_{m} and the Boltzmann fits for real systems, we experimentally recorded the gating currents in the inactivation-removed Shaker W434F potassium channel construct (Sh-IR W434F) and in the voltage-gated sodium channel (Na_{V}1.4), with the ionic currents blocked with tetrodotoxin, and evaluated their normalized QV curves (Fig. 2). V_{1/2} for each was obtained by fitting the Boltzmann equation (Eq. 1), whereas V_{m} was obtained by calculating the area between the curve and the Q axis (see Materials and methods). Previously published values of Q_{max} for each channel (Schoppa et al., 1992; Hirschberg et al., 1995; Aggarwal and MacKinnon, 1996; Seoh et al., 1996) were used for the free energy calculations using the median method. We find that the free energy values computed via the median method were four to eight times the values computed via the Boltzmann method (Table 1).

Occasionally for potassium channels, the sigmoid activation curves have been fitted to a fourth-power Boltzmann equation (Zagotta et al., 1994b; Smith-Maxwell et al., 1998b; Gonzalez et al., 2000, 2001), assuming that identical conformational changes occur in a single step in each of the four identical voltage-sensors, which “move” independent of the other. This fit provides a much closer value of the free energy change to that obtained from the median method for Sh-IR W434F, but for Na_{V}1.4 the difference between the two estimates continues to be large (Table 1). Fitting the QV curve of Na_{V}1.4 to the product of four nonidentical Boltzmann terms (which takes into consideration the four nonidentical domains of the protein) would be ill constrained. We fitted two other model dependent sigmoidal functions (Perozo et al., 1994; Rodríguez et al., 1998), both of which resulted in free energy estimates significantly different from the median method for Sh-IR W434F as well as Na_{V}1.4 (Table 1).

For the Shaker potassium channel, elaborate multistep kinetic models have been proposed by various groups (Bezanilla et al., 1994; Zagotta et al., 1994a; Schoppa and Sigworth, 1998) based on single channel conductances, gating currents, and macroscopic ionic currents. Here, we consider the model proposed by Zagotta et al. (1994a; ZHA model).

(SCHEME 1)According to this model, each voltage sensor independently activates in two discrete steps (with equilibrium constants *K*_{1} and *K*_{2}), and when all the voltage-sensors are activated, a concerted transition gates the ion channel pore open (equilibrium constant *K*_{L}). All the equilibrium constants have an exponential voltage dependence: *i* = 1, 2, or L). The net chemical free energy change in going from the reference state to the final open state, *k _{B}* and the inclusion of Faraday constant [F] in the numerical free energy calculations.) Using the parameters reported, we find that the free energy estimate from the ZHA model parameters and the median method are nearly identical (Table 2). Next, we considered the ILT mutant of the Shaker potassium channel (Smith-Maxwell et al., 1998a; Ledwell and Aldrich, 1999), which is characteristic in that it results in a QV curve with a secondary phase, widely separated from the primary (Fig. 2), and thus cannot be fitted to a Boltzmann curve. However, the value of chemical free energy change for this mutant computed using the median method from the numerically simulated QV curve is again almost the same as that obtained from the kinetic model parameters. The fitting exercise described herein thus shows that the median method is a physically consistent, experimentally feasible method to accurately estimate voltage-dependent energetics in the system.

### Free energy estimation using the median method in the case where certain steps are voltage independent

Next we considered models of ion channel activation with voltage-independent steps to illustrate how our inferred energetics can be modified by consideration of such transitions. Let us assume that a voltage-independent transition occurs in the middle of a linear transition pathway, involving six states with the final state being the open state O_{f} (Fig. 3 A, scheme I). Because all of the steps beyond it are voltage driven, as we increase the voltage, the system will be pushed further and further to the right until it finally saturates in O_{f}. For this system, the chemical free energy difference between the initial and final states, _{max}V_{m} as described previously (see Theory section A physical interpretation of _{m}).

To test this numerically, we generated a large number of QV curves, based on this model, for different values of *z*_{2} and *z*_{2} alters the maximum gating charge transferred in the process of activation. From each simulated curve, we computed the median voltage of activation, V_{m}, by integrating the area between the curve and the ordinate (Q) axis (_{max}FV_{m} against _{max}. The surface plot of *z*_{2} and log *z*_{2}. However, V_{m} shows a logarithmic dependence on *z*_{2} (Fig. 3 D). _{max} in the system, V_{m} is linearly related to _{max} is altered without changing any of the equilibrium constants, _{m} will be altered because Q_{max}FV_{m} =

Next, we modified scheme I (Fig. 3 A) by making the last (rather than the middle) step of the pathway voltage independent (Fig. 4 A, scheme II). Now at saturating depolarization, both the final states, C_{5} and O_{f}, will be occupied, their relative occupancies being determined by the voltage-independent equilibrium constant of the last transition (_{max}V_{m} reports the change in chemical free energy in taking the system from its initial reference state to a saturating condition, which is an ensemble rather than a single state. This brings us to an important point: there are two nonequivalent definitions of free energy change in the system. One is the change in the free energy of the ensemble, whereas the other is the free energy difference between the unique final and initial states of the system (_{max}V_{m}, whereas to estimate the latter, we rewrite Eq. 21:_{max} and V_{m}, would let us estimate the net free energy change associated with transferring the channel from its initial state C_{1} to its final state O_{f}.

We illustrate this point through numerical simulations. We calculated the V_{m} for several QV curves, based on scheme II (Fig. 4 A), using different values of _{max}FV_{m} against _{max}FV_{m} and

Next, we consider scheme III, which is a six-state linear model with two final open states, the transition between them being voltage independent (Fig. 4 A). For this case, the Q_{max}V_{m} measure is equal to_{max}V_{m} will be the free energy change of the ensemble but not the chemical free energy difference between the first and last states (

### Effects of cooperativity and ion channel inactivation

In the models considered thus far, the median measure of the chemical free energy change of the system correlates well with the model-based estimates. To further check and validate these correlations, we performed numerical simulations on an allosteric model of ion channel activation. Such models have been used with success to understand the voltage-dependent gating of the BK channels (Horrigan and Aldrich, 1999, 2002). We used an MWC model (Monod et al., 1965) for our simulations here.

Consider scheme IV, a 10-state MWC model of ion channel activation as shown in Fig. 5 A. The channel is assumed to comprise four identical voltage-sensing modules and a single pore domain. *K*_{V} is the equilibrium constant of activation of each of the four voltage-sensing modules, whereas *K*_{P} is the equilibrium constant of the pore-opening transition. Both of the equilibrium constants are assumed to have an exponential voltage dependence: *i* = V, P) where *z _{i}* reflects the voltage dependence of each transition and

_{1}at hyperpolarizing voltages and end in state O

_{f}at highly depolarizing voltages. The chemical free energy difference between the states C

_{1}and O

_{f},

We generated several QV curves for scheme IV (Fig. 5 A) with different values of *z*_{V} and θ when *z*_{P} is not 0. The free energy change calculated from median voltage of activation in each case was compared with _{max} for scheme IV (Fig. 5 A) is 4*z*_{V} + *z*_{P}. The plot of Q_{max}FV_{m} against *z*_{P} ≠ 0), the two free energy measures are equal, dependent on θ but independent of Q_{max}. This illustrates that the presence of multiple open states or latent charge movement between the open states does not limit the applicability of the median method to estimate the chemical free energy change.

When the pore-opening process is voltage independent, at saturating voltages, both states C_{f} and O_{f} will be populated, their relative proportions being determined by the magnitude of the equilibrium constant, _{0} and O_{0} will be populated, their relative proportions being determined by _{0} to be our reference state, we see that

For illustration purposes, we compared numerically simulated values of Q_{max}FV_{m} and *z*_{P} = 0) at different values of the model parameters. When _{max}FV_{m}. Alternatively, when both θ and

We finally consider the case of channels that undergo inactivation on sustained depolarization (Fig. 6 A, scheme V). The system can be modeled by considering that there are three specialized units, V, P, and I, each capable of existing in two conformations. Conformational change of V (V_{R} to V_{A}) governs the initial voltage-sensing steps of the channel, that of P (P_{C} to P_{O}) is the gate opening transition, and that of I (I_{R} to I_{A}) controls the inactivation process. Scheme IV (Fig. 6 A) has 12 parameters. This cubic model of an inactivating channel can be represented as a nested coupled model as depicted in Fig. 6 B. Each of three structural units has an intrinsic activation constant described by a voltage-independent component, *z _{i}* (where

*i*is V, P, or inact). Additionally there are three coupling parameters, θ

_{VP}, θ

_{PI}, and θ

_{IV}, which are voltage independent and depict the pairwise interactions between the structural units. Each of the 12 parameters of the cubic model (Fig. 6 A) can be uniquely related to the parameters of Fig. 6 B using the principle of microscopic reversibility.

If all of the equilibrium constants in this model have a positive voltage dependence, at highly depolarizing voltages, the system will saturate in the state V_{A}P_{O}I_{A}, whereas at very low voltages, only the initial reference state, V_{R}P_{C}I_{R}, will be populated. Q_{max}V_{m} (V_{m} being obtained from the QV curve) will thus be the chemical free energy difference between the initial and final states of the system.

Next, assume that the inactivation process is voltage independent. As a result, all of the equilibrium constants between two connected states differing in the conformation of I (Fig. 6 A) will be voltage independent, i.e., *K*_{9}, *K*_{10}, *K*_{11}, and *K*_{12} do not change with voltage. Because all other transitions are voltage dependent, both states V_{A}P_{O}I_{R} and V_{A}P_{O}I_{A} will be populated at depolarized potentials, and at very low voltages, both V_{R}P_{C}I_{R} and V_{R}P_{C}I_{A} will exist. This situation is similar to our description of the MWC model, with voltage-independent pore opening transitions, in that multiple states are populated at both hyperpolarized and depolarized voltages.

Thus, Eq. 18 here transforms to_{max}V_{m} is the change in the chemical free energy of the ensemble but not the chemical free energy difference between a unique final and initial states of the system, _{max}) in this system, at depolarizing voltages, is _{max}V_{m} is small when I_{max} is close to unity and continues to increase with decreasing values of I_{max}. These points are illustrated through simulations performed using scheme V (Fig. 6 A). Q_{max}FV_{m} is plotted against _{max} is ∼1 (Fig. 6 D). Thus, our overall proposition that Q_{max}FV_{m} is the chemical free energy of the ensemble remains valid even in inactivating channels.

## DISCUSSION

In summary, we have derived a general expression for the change in the chemical free energy change associated with the activation of a voltage-dependent ion channel, _{max} is the maximum amount of gating charge displacement, F is the Faraday constant, and V_{m} is the median voltage of activation estimated from a QV curve. V_{m} is mathematically defined by Eq. 11. It is the voltage at which the area bounded by the QV curve, V_{m} axis, and Q = 0 axis becomes equal to that bounded by the QV curve, V_{m} axis, and the Q = Q_{max} axis, as depicted pictorially in Fig. 1 B. We have considered various classes of models with cooperativity, voltage-independent transitions, multiple open states, latent charge movement, etc. and discussed the measures of free energy change in each case. We demonstrate that the V_{m} measure of free energy change is independent of the nature of the transition pathway and of the symmetry relations within the protein (i.e., whether parts of them are identical or not). In some instances where initial or terminal steps of the transition pathway are voltage independent, certain correction factors need to be added to the V_{m} measure to calculate the free energy difference between the first and last state of the channel. General relations between the median estimates of free energy changes and the chemical free energy differences between the first and last states of the system in six broad classes of ion channel gating models are listed in Table 3.

Although these relations are fairly general, their application to understand the energetic relationships in ion channels should be performed with careful deliberation of the assumptions on the basis of which the relations are derived and the behavior of the specific system of interest. The median estimate of free energy is no substitute for full quantitative models of channel gating. However, it might serve as a useful check/constraint of the parameters for detail models of ion channel gating. Also, this measured free energy change includes free energy contributions from the intrinsic stabilities of the different structural units as well as interactions between them. Parsing out these individual components via measurement of QV curves alone is not straightforward and might require the application of site-specific measurements along with more detailed analyses such as the recently proposed χ-value analysis (Chowdhury and Chanda, 2010).

The median voltage of charge movement (V_{m}) is a more accurate and physically meaningful free energy correlate than V_{1/2}. The two parameters will be identical when the QV curve is symmetric, but the difference between them can be easily appreciated in cases in which the QV curve is nonsymmetric and/or shows multiple phases (Perozo et al., 1994; Ledwell and Aldrich, 1999; Tao et al., 2010; Lacroix and Bezanilla, 2011). We estimate that the chemical free energy change for the Shaker potassium channels is approximately −14 kcal/mol, which suggests that most electrical energy is being efficiently used to drive conformational change in the ion channel.

Our derivations establish QV as the direct free energy correlate of the free energy change in a voltage-dependent system. In many instances, P_{O}-V curves have been used to obtain empirical estimates of free energy of activation of voltage-dependent ion channels (Monks et al., 1999; Li-Smerin et al., 2000a; Hackos et al., 2002; Yifrach and MacKinnon, 2002; Sadovsky and Yifrach, 2007). Such energy measures are substantially different from those proposed in this study. Operationally, a P_{O}-V–based estimate of the free energy change can be interpreted as the free energy difference between the ensemble of closed states of the system at the threshold voltage of pore opening and the final ensemble of open states. This also implies that a P_{O}-V curve does not fully capture the activation energetics of a voltage-dependent ion channel that activates in multiple steps.

On a final note, it must be mentioned that, in a case in which the Q_{max} of a channel is not measured, the energetic effect of a mutation on the activation energetics can be enumerated via the V_{m} of its QV curve alone, under the assumption that Q_{max} is unaltered by the mutation. If the QV curve of the mutant channel is symmetric, V_{1/2} might be used instead of V_{m} (both will be identical). This implies that considering the Boltzmann slope in the energy terms introduces a theoretical anomaly in the expression for the change in free energy in cases beyond the realm of a two-state approximation. Therefore when Q_{max} is not measured, an empirical estimate of the energetic effect of a mutation on the voltage-gated ion channel can be assessed simply via comparison of the V_{m} value of its QV curve with that of a control (Soler-Llavina et al., 2006; Gagnon and Bezanilla, 2010; Muroi et al., 2010; Xu et al., 2010). The obvious caveats to this approach are cases in which mutations affect Q_{max}. In such scenarios, to quantify the magnitude of energetic perturbation, apart from the QV curve, Q_{max} for each charge-altering mutation needs to be measured separately.

## Acknowledgments

We thank Deborah L. Capes for providing the gating current data for the rNa_{V}1.4 channels and Dr. Meyer Jackson and Dr. Marcel Goldschen for their valuable comments.

This project was supported by funds from the National Institutes of Health (grant GM084140) and the Shaw Scientist Award to B. Chanda.

Christopher Miller served as editor.

- Submitted: 15 September 2011
- Accepted: 22 November 2011

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/).