On the theory of ion transport across the nerve membrane. VI. Free energy and activation free energies of conformational change

A Microscopic Capacitor Model of Voltage Coupling in Membrane Proteins: Gating Charge Fluctuations in Ci-VSD

The voltage sensitivity of membrane proteins is reflected in the response of the voltage sensing domains (VSDs) to the changes in membrane potential. This response is implicated in the displacement of positively charged residues, associated with the conformational changes of VSDs. The displaced charges generate nonlinear (i.e., voltage-dependent) capacitance current called the gating current (and its corresponding gating charge), which is a key experimental quantity that characterizes voltage activation in VSMP. However, the relevant theoretical/computational approaches, aimed to correlate the structural information on VSMP to electrophysiological measurements, have been rather limited, posing a broad challenge in computer simulations of VSMP. Concomitant with the development of our coarse-graining (CG) model of voltage coupling, we apply our theoretical framework for the treatments of voltage effects in membrane proteins to modeling the VSMP activation, taking the VSDs (Ci-VSD) derived from the Ciona intestinalis voltage sensitive phosphatase (Ci-VSP) as a model system. Our CG model reproduces the observed gating charge of Ci-VSD activation in several different perspectives. In particular, a new closed-form expression of the gating charge is evaluated in both nonequilibrium and equilibrium ways, while considering the fluctuation-dissipation relation that connects a nonequilibrium measurement of the gating charge to an equilibrium measurement of charge fluctuations (i.e., the voltage-independent linear component of membrane capacitance). In turn, the expression uncovers a novel link that connects an equilibrium measurement of the voltage-independent linear capacitance to a nonequilibrium measurement of the voltage-dependent nonlinear capacitance (whose integral over voltage is equal to the gating charge). In addition, our CG model yields capacitor-like voltage dependent free energy parabolas, resulting in the free energy difference and the free energy barrier for the Ci-VSD activation at "zero" (depolarization) membrane potential. Significantly, the resultant voltage dependent energetics enables a direct evaluation of capacitance-voltage relationship (C-V curve) as well as charge-voltage relationship (Q-V curve) that is in a good agreement with the observed measurement of Ci-VSD voltage activation. Importantly, an extension of our kinetic/thermodynamic model of voltage dependent activation in VSMP allows for novel derivations of voltage-dependent rate constants, whose parameters are expressed by the intrinsic properties of VSMP. These novel closed-form expressions offer a physicochemical foundation for the semiempirical Eyring-type voltage dependent rate equations that have been the cornerstone for the phenomenological (kinetic) descriptions of gating and membrane currents in the mechanistic study of ion channels and transporters. Our extended theoretical framework developed in the present study has potential implications on the roles played by gating charge fluctuations for the spike generations in nerve cells within the framework of the Hodgkin-Huxley-type model.

A geometric comparison of single chain multi-state models of ion channel gating

Multi-state models of ion channel gating have been used extensively, but choosing optimally small yet sufficiently complex models to describe particular experimental data remains a difficult task. In order to provide some insight into appropriate model selection, this paper presents some basic results about the behavior of solutions of multi-state models, particularly those arranged in a chain formation. Some properties of the eigenvalues and eigenvectors of constant-rate multi-state models are presented. A geometric description of a three-state chain is given and, in particular, differences between a chain equivalent to an Hodgkin-Huxley model and a chain with identical rates are analyzed. One distinguishing feature between these two types of systems is that decay from the open state in the Hodgkin-Huxley model is dominated by the most negative eigenvalue while the identical rate chain displays a mix of modes over all eigenvalues.

Thermodynamically equivalent silicon models of voltage-dependent ion channels

We model ion channels in silicon by exploiting similarities between the thermodynamic principles that govern ion channels and those that govern transistors. Using just eight transistors, we replicate--for the first time in silicon--the sigmoidal voltage dependence of activation (or inactivation) and the bell-shaped voltage-dependence of its time constant. We derive equations describing the dynamics of our silicon analog and explore its flexibility by varying various parameters. In addition, we validate the design by implementing a channel with a single activation variable. The design's compactness allows tens of thousands of copies to be built on a single chip, facilitating the study of biologically realistic models of neural computation at the network level in silicon.

An improved non-linear thermodynamic model of voltage-dependent ionic currents

Thermodynamic models of ionic currents are used to deduce exact functional form of rate constant. It is first assumed that free energy depends linearly on voltage. A major criticism of this approach is that time constant can reach arbitrarily small values resulting in an aberrant behavior. Recently, non-linear effects of electric field on the free energy were considered to solve this problem based on T-type calcium channel. In this study, we show that the current model approximately captures the voltage-dependence of the time constant only in a specific range of voltage and does not provide saturation of the time constant outside this range. Then, we propose an improved non-linear thermodynamic model and illustrate its applicability based on T-type calcium channel.

Nonlinear thermodynamic models of voltage-dependent currents

Hodgkin and Huxley provided the first quantitative description of voltage-dependent currents and adjusted their model to experimental data using empirical functions of voltage. A physically plausible formalism was proposed later by assuming that transition rates depend exponentially on a free-energy barrier, by analogy with the theory of reaction rates. It was also assumed that the free energy depends linearly on voltage. This thermodynamic formalism can accurately describe many processes, but the resulting time constants can be arbitrarily fast, which may also lead to aberrant behavior. We considered here a physically plausible solution to this problem by including nonlinear effects of the electrical field on the free energy. We show that including effects such as mechanical constraints, inherent to the structure of the ion channel protein, leads to more accurate thermodynamic models. These models can account for voltage-dependent transitions that are rate-limited in a given voltage range, without invoking additional states. We illustrate their applicability to fit experimental data by considering the case of the T-type calcium current in thalamic neurons.

Recognition and processing of randomly fluctuating electric signals by Na,K-ATPase

Previous work has shown that Na,K-ATPase of human erythrocytes can extract free energy from sinusoidal electric fields to pump cations up their respective concentration gradients. Because regularly oscillating waveform is not a feature of the transmembrane electric potential of cells, questions have been raised whether these observed effects are biologically relevant. Here we show that a random-telegraph fluctuating electric field (RTF) consisting of alternating square electric pulses with random lifetimes can also stimulate the Rb(+)-pumping mode of the Na,K-ATPase. The net RTF-stimulated, ouabain-sensitive Rb+ pumping was monitored with 86Rb+. The tracer-measured, Rb+ influx exhibited frequency and amplitude dependencies that peaked at the mean frequency of 1.0 kHz and amplitude of 20 V/cm. At 4 degrees C, the maximal pumping activity under these optimal conditions was 28 Rb+/RBC-hr, which is approximately 50% higher than that obtained with the sinusoidal electric field. These findings indicate that Na,K-ATPase can recognize an electric signal, either regularly oscillatory or randomly fluctuating, for energy coupling, with high fidelity. The use of RTF for activation also allowed a quantitative theoretical analysis of kinetics of a membrane transport model of any complexity according to the theory of electroconformational coupling (ECC) by the diagram methods. A four-state ECC model was shown to produce the amplitude and the frequency windows of the Rb(+)-pumping if the free energy of interaction of the transporter with the membrane potential was to include a nonlinear quadratic term. Kinetic constants for the ECC model have been derived. These results indicate that the ECC is a plausible mechanism for the recognition and processing of electric signals by proteins of the cell membrane.

Voltage-noise-induced transitions in electrically excitable membranes

A quantitative study of the steady-state behavior of the sodium and potassium conductance for the Hodgkin-Huxley axon under the influence of an externally driven voltage noise is reported. The dichotomous Markov noise (random telegraph signal) considered allows for an exact evaluation of the stationary probability density of the conductances. Phase diagrams are constructed to represent the response of the system as a function of the amplitude and the correlation time of the noise. The results obtained for the Hodgkin-Huxley axon are compared with some molecular models used in the literature.

Biological effects of oscillating electric fields: role of voltage-sensitive ion channels

An alternating component of potential across the membrane of an excitable cell may change the membrane conductance by interacting with the voltage-sensing charged groups of the protein macromolecules that form voltage-sensitive ion channels. Because the probability than a voltage sensor is in a given state is a highly nonlinear function of the applied electric field, the average occupancy of a particular state will change in an oscillating electric field of sufficient magnitude. This "rectification" at the level of the voltage sensors could results in conformational changes (gating) that would modify channel conductance. A simplified two-state model is examined where the relaxation time of the voltage sensor is assumed to be considerably faster than the fastest changes of ionic conductance. Significant changes in the occupancy of voltage sensor states in response to an applied oscillating electric field are predicted by the model.

On the physico-chemical basis of voltage-dependent molecular gating mechanisms in biological membranes

The possible nature and theoretical treatment of electric field-induced molecular processes in a membrane are examined. Special attention is given to fairly fast switching phenomena as reflected by asymmetry currents as well as ionic gating in squid axon and similar systems. The apparent charge displacement associated with the underlying mechanisms is argued to be brought about by conformational transitions of integral macromolecular structures. Under these circumstances, voltage changes can actaully control the functional state of membranes by direct interference with conformational equilibria. A basic model is quantitatively discussed and shown to account for certain observed asymmetry currents. Effects due to temperature, pressure, or chemical interactions can be readily described. It is indicated how more complicated voltage-dependent membrane processes may be approached along these lines.

Dipole moment, enthalpy, and entropy changes of Hodgkin-Huxley type kinetic units

Dipole moment, enthalpy, and entropy changes were calculated for hypothetical structural units which control the opening and closing of ionic channels in axon membranes. The changes of these thermodynamic functions were calculated both for activation (transition to intermediate complex) and for the structural transformation as a whole. The calculations are based on the experimentally determined Q10 values and the empirical formulae for the rate constants (alpha's and beta's) as functions of membrane potentials in Hodgkin-Huxley type models. From the calculated thermodynamic functions we suggest that the specific structural units of the axon membranes are probably of macromolecular (possible protein-like) dimensions with large dipole moments (hundreds of debyes). The calculated dipole moment changes of a single structural unit indicate that in many cases these dipole moments saturate at strong depolarizations or hyperpolarizations. The transitions in structural units show substantial activation enthalpies and entropies but the net enthalpy and entropy changes are practically negligible for the transition as a whole, i.e. the structural units presumably undergo displacements. While the calculated dipole moment changes associated with structural transformations in Loligo and Myxicola show similar potential dependencies, those for Rana usually show a different behavior. The relevance of the dipole moment changes to gating currents is discussed.