Current fluctuations in discrete transport systems far from equilibrium. Breakdown of the fluctuation dissipation theorem

A recently developed theoretical approach to transport fluctuations around stable steady states in discrete biological transport systems is used in order to investigate general fluctuation properties at nonequilibrium. An expression for the complex frequency dependent admittance at nonequilibrium is derived by calculation of the linear current response of the transport systems to small disturbances in the applied external voltage. It is shown that the Nyquist or fluctuation dissipation theorem, by which at equilibrium the macroscopic admittance or linear response can be expressed in terms of fluctuation properties of the system, breaks down at nonequilibrium. The spectral density of current fluctuations is decomposed into one term containing the macroscopic admittance and a second term which is bilinear in current. This second term is generated by microscopic disturbances, which cannot be excited by external macroscopic perturbations. At special examples it is demonstrated that this second term is decisive for the occurrence of excess noise e.g. the 1/f(2)-Iorentzian noise generated by the opening and closing of nerve channels in biological membranes.     

Irreversible thermodynamic analysis of electrophoretic light scattering experiments

A fluctuation theory for electrolyte solutions is developed based on the coupling between the equations of nonequilibrium thermodynamics and the Poisson equation. The resulting fluctuation theory is applied to the analysis of electrophoretic light scattering. It is shown that in a binary electrolyte solution (two ionic species), the Doppler shift is not determined by the electrical mobility of either ion, but depends instead on the rate of change of transference number with salt concentration. In addition the ionic relaxation time is shown to be proportional to the conductivity of the solution.     

Amendments to the theory underlying Ussing chamber data of chloride ion secretion after bacterial enterotoxin exposure

Bacterial enterotoxins may cause life-threatening diarrhoeal fluid loss in part because they stimulate enterocytes to secrete fluid into the small intestine as well as preventing normal fluid uptake. Abnormal chloride ion secretion is believed to provide the osmotic driving force for the inappropriate fluid movement. Evidence for enhanced chloride secretion consists of isotopic flux measurements in Ussing chambers, the standard apparatus for permeation studies. Flux from the lumen of the intestine is assumed to be determined solely by absorptive processes and flux towards the lumen solely by secretory processes. Bacterial enterotoxin increased flux towards the lumen is taken as an evidence of enhanced secretion. Examination of the flux equation solutions shows that the existing theoretical treatment of the Ussing chamber consists of the super-imposition of two contradictory unidirectional models. In contrast, the present analysis shows that a measured 'unidirectional' flux contains information both about absorptive and secretory processes, regardless of which flux is measured. Reciprocity is predicted for the fluxes, as decreases in the absorptive processes will cause increases in apparent secretory flux. Data from the literature show that mucosal-to-serosal chloride ion flux in rabbit ileum after exposure to secretagogues correlates inversely and highly significantly (r=0.74, n=17, p<0.001) with increases in serosal-to-mucosal chloride ion flux. As a category of evidence, flux data do not provide conclusive evidence of enhanced chloride secretion after exposure to enterotoxins, since an apparently enhanced serosal-to-mucosal flux would also be noted after inhibition of the mucosal-to-serosal flux. As interruption of absorptive processes can be misinterpreted as enhanced secretion in the Ussing chamber, this is a serious deficiency in the evidence for direct enterotoxin enhancement of the intestinal chloride ion channel as a basis for diarrhoeal disease.     

Electrostatic radius of the gramicidin channel determined from voltage dependence of H+ ion conductance.

The results of Decker and Levitt (1987) suggest that the conductance of H+ ion through the gramicidin channel is limited primarily by diffusion in the bulk solution at the channel mouth. It is assumed in this paper that the H+ conductance is 100% diffusion limited. This means that all the factors that influence the H+ flux are external to the channel and are presumed to be known. In particular, the diffusion coefficient of H+ in this region is assumed to be equal to the bulk solution value and the only force acting on the ion is that due to the applied voltage. A model of the H+ flux is derived, based on the Nernst-Planck equation. It has three adjustable parameters: the electrostatic radius, the capture distance, and the radius of the H+ ion. The acceptable range of the parameters was determined by comparing the predictions of the model with the experimental measurements of the H+ conductance at pH 3.75. The best fit was obtained for an electrostatic radius in the range 2.3-2.7 A. This is in good agreement with earlier predictions (2.5 A) based on the assumption that the dielectric constant of the channel water is equal to that of bulk water. The addition of 1 M choline Cl- (an impermeant) increases the H+ current at low voltage and decreases it at high voltage. The increase can be explained by the small surface charge that results from the separation of charge produced by exclusion of the large choline cation (relative to Cl-) from the membrane surface. The decrease at high voltages can be accounted for by the change in the profile of the applied potential produced by the increase in ionic strength.     

Conformational model for ion permeation in membrane channels: a comparison with multi-ion models and applications to calcium channel permeability.

The permeation properties of ion channels existing in several conductive states were analyzed. Each state was represented by the one-ion model. A special emphasis was placed on features, assumed to be indicative of a multi-ion mode of channel occupancy such as a deviation of concentration dependence of channel conductance from the Michaelis-Menten equation, an anomalous mole fraction effect, a strong voltage dependence of ion block and coupling of unidirectional fluxes (anomalous Ussing flux ratio). The conformational model was shown to have all these properties. The ion permeation through voltage-sensitive calcium channels fulfilled all the characteristics of the model proposed.     

Single channel gating events in tracer flux experiments I. Theory

Tracer ion flux measurements are a commonly used method for studying ion transport through membranes of cellular systems, where the rate of ion flow is determined by gating processes which control the opening and closing of transmembrane channels. Due to recent advances in the theoretical analysis of tracer flux from or into closed membrane structures (CMS), the mechanism of gating reactions can, in principle, be derived from flux data. A physically well founded analysis is presented for the dependence of the total tracer ion content of a collection of CMS on the gating processes. For functionally uncoupled gating units a mean single channel flux contribution [equation, see text] can be defined, where k is the intrinsic single channel flux coefficient, t the time over which flux is measured, and p(tau,t) is the probability that a given channel was open for a total period tau during t. This quantity reflects the mean time course of the tracer content due to flux through a single channel. Expressions for are derived that explicitly take into account a distribution in the lifetime of open channels. On the basis of the results, kinetic and thermodynamic parameters of multiphasic gating reactions can be determined from the time course of the overall tracer content in a colleciion of CMS.     

Stochastic properties of ion channel openings and bursts in a membrane patch that contains two channels: evidence concerning the number of channels present when a record containing only single openings is observed.

If a single ion channel record is observed in which two ion channels are never simultaneously open, then it is often of interest to know whether the observations indeed arose from the activity of only one ion channel. This question can be answered if it is possible to calculate the distribution of the duration of runs of single openings in a membrane patch that contains two active channels. If the observed run of single openings is much longer than that expected for a patch with two channels it is likely that only one channel was active. An approximate method is presented for calculating the distribution of the duration of runs of single openings in a patch with two active channels; this method has the advantage that it can be calculated from observable quantities, and requires no knowledge of the details of the ion-channel mechanism or its rate constants. The accuracy of this approximation is tested by exact calculations of the properties of runs of single openings, and of single bursts, for two specific mechanisms and a large range of rate constants. The approximation is good in all cases in which openings occur singly, or in closely spaced bursts. If, as is common in practice, openings occur in clusters that are separated by long shut periods, then overlap of clusters from two different channels may be detected, if no double opening is produced, as a period in the middle of a cluster in which the probability of being open doubles. The results derived here can be applied to such a period to test whether it results from the simultaneous activity of two channels, rather than from a change in the properties of a single channel.     

Application of the ensemble nonequilibrium response spectroscopy to shaker potassium channel gating

Standard electrophysiology techniques study relaxation transients in voltage-gated ion channels generated by discrete voltage steps. The nonequilibrium response spectroscopy involves analyzing responses to fluctuating potentials. We apply the ensemble NRS method to gating kinetics of Shaker potassium ion channels. We evaluate various proposed Markov models of channel gating from the nonequilibrium response viewpoint. These new NRS protocols can be used to test otherwise indistinguishable models or improve estimates for parameters of channel kinetics models.     

A ratiometric imaging method for mapping ion flux densities

A heterogeneous distribution of ion channels on the cell surface is a prerequisite for several cellular functions. Thus, there has been considerable interest in methods allowing the mapping of ion channel distributions. Here we report on a novel ratiometric imaging technique appropriate to measure spatially resolved ion flux signals by using ion sensitive dyes. However, given that certain relevant cell properties like the surface to volume ratio may exhibit significant spatial heterogeneities, the local influx signal cannot be interpreted as a measure of the local open channel concentration or flux density. To overcome this problem, we suggest an internal normalization procedure, which, in analogy to, but clearly distinct from, well-established ratioing techniques, eliminates effects which would otherwise obscure the desired result. Ratioing is performed on flux signals from a given cell, triggered by two different, subsequent stimuli. If the two stimuli address different ion channels, the flux density distribution caused by two channel types can be determined relative to each other. In cases where one of the stimuli triggers a spatially homogeneous flux signal, ratioing yields an ion flux density map for a given channel type. Thus distribution patterns of ion channels active during a given stimulus may be derived.     

Flux ratio theorems for nonlinear membrane transport under nonstationary conditions

We extend flux ratio theorems concerning ratios of unidirectional flux transients passed (in complementary experiments) through a medium of spatially inhomogeneous transport properties pertaining to diffusion, migration and temporary trapping of the transported substance. Any nonlinearity in the transport equations leads to a breakdown of the Ussing flux ratio theorem pertaining to all times. An integrated flux ratio theorem is proved for the case when the nonlinearity is in the kinetics of trapping, as when trapping sites can be saturated. The new theorem is shown to fail when the nonlinearity is due to a concentration-dependence of the diffusion coefficient, as in facilitated transport. The nature of a nonlinearity in membrane transport can therefore be elucidated experimentally by the use of the integrated flux ratio.     

Kudos to Reviewers for the JGP: You Make Our Science Better

Rate-equilibrium free energy relationship (REFER) analysis provides information on transition-state structures and has been applied to reveal the temporal sequence in which the different regions of an ion channel protein move during a closed–open conformational transition. To date, the theory used to interpret REFER relationships has been developed only for equilibrium mechanisms. Gating of most ion channels is an equilibrium process, but recently several ion channels have been identified to have retained nonequilibrium traits in their gating cycles, inherited from transporter-like ancestors. So far it has not been examined to what extent REFER analysis is applicable to such systems. By deriving the REFER relationships for a simple nonequilibrium mechanism, this paper addresses whether an equilibrium mechanism can be distinguished from a nonequilibrium one by the characteristics of their REFER plots, and whether information on the transition-state structures can be obtained from REFER plots for gating mechanisms that are known to be nonequilibrium cycles. The results show that REFER plots do not carry information on the equilibrium nature of the underlying gating mechanism. Both equilibrium and nonequilibrium mechanisms can result in linear or nonlinear REFER plots, and complementarity of REFER slopes for opening and closing transitions is a trivial feature true for any mechanism. Additionally, REFER analysis provides limited information about the transition-state structures for gating schemes that are known to be nonequilibrium cycles.     

Strong electrolyte continuum theory solution for equilibrium profiles, diffusion limitation, and conductance in charged ion channels.

The solution for the ion flux through a membrane channel that incorporates the electrolyte nature of the aqueous solution is a difficult theoretical problem that, until now, has not been properly formulated. The difficulty arises from the complicated electrostatic problem presented by a high dielectric aqueous channel piercing a low dielectric lipid membrane. The problem is greatly simplified by assuming that the ratio of the dielectric constant of the water to that of the lipid is infinite. It is shown that this is a good approximation for most channels of biological interest. This assumption allows one to derive simple analytical expressions for the Born image potential and the potential from a fixed charge in the channel, and it leads to a differential equation for the potential from the background electrolyte. This leads to a rigorous solution for the ion flux or the equilibrium potential based on a combination of the Nernst-Planck equation and strong electrolyte theory (i.e., Gouy-Chapman or Debye-Huckel). This approach is illustrated by solving the system of equations for the specific case of a large channel containing fixed negative charges. The following characteristics of this channels are discussed: anion and mono- and divalent cation conductance, saturation of current with increasing concentration, current-voltage relationship, influence of location and valence of fixed charge, and interaction between ions. The qualitative behavior of this channel is similar to that of the acetylcholine receptor channel.     

Complexity and demographic stability in population models

This article is concerned with relating the stability of a population, as defined by the rate of decay of fluctuations induced by demographic stochasticity, with its heterogeneity in age-specific birth and death rates. We invoke the theory of large deviations to establish a fluctuation theorem: The demographic stability of a population is positively correlated with evolutionary entropy, a measure of the variability in the age of reproducing individuals in the population. This theorem is exploited to predict certain correlations between ecological constraints and evolutionary trends in demographic stability, namely, (i) bounded growth constraints--a uni-directional increase in stability, (ii) unbounded growth constraints (large population size)--a uni-directional decrease in stability, (iii) unbounded growth constraints (small population size)--random, non-directional change in stability. These principles relating ecological constraints with trends in demographic stability are shown to be far reaching generalizations of the tenets derived from classical studies of stability in an evolutionary context. These results thus provide a new conceptual framework for explaining patterns of variation in population numbers observed in natural populations.     

Electrodiffusion of ions approaching the mouth of a conducting membrane channel.

The movement of ions in the aqueous medium as they approach the mouth (radius a) of a conducting membrane channel is analyzed. Starting with the Nernst-Planck and Poisson equations, we derive a nonlinear integrodifferential equation for the electric potential, phi(r), a less than or equal to r less than infinity. The formulation allows deviations from charge neutrality and dependence of phi(r) on ion flux. A numerical solution is obtained by converting the equation to an integral equation that is solved by an iterative method for an assumed mouth potential, combined with a shooting method to adjust the mouth potential until the numerical solution agrees with an asymptotic expansion of the potential at r-a much greater than lambda (lambda = Debye length). Approximate analytic solutions are obtained by assuming charge neutrality (Läuger, 1976) and by linearizing. The linear approximation agrees with the exact solution under most physiological conditions, but the charge-neutrality solution is only valid for r much greater than lambda and thus cannot be used unless a much greater than lambda. Families of curves of ion flux vs. potential drop across the electrolyte, phi(infinity)-phi (a), and of permeant ion density at the channel mouth, n1(a), vs. flux are obtained for different values of a/lambda and S = a d phi/dr(a). If a much greater than lambda and S = O, the maximum flux (which is approached when n1(a)----0) is reduced by 50% compared to the value predicted by the charge-neutrality solution. Access resistance is shown to be a factor a/[2 (a + lambda)] times the published formula (Hille, 1968), which was derived without including deviations from charge neutrality and ion density gradients and hence does not apply when there is no counter-ion current. The results are applied to an idealized diffusion-limited channel with symmetric electrolytes. For S = O, the current/voltage curves saturate at a value dependent on a/lambda; for S greater than O, they increase linearly for large voltage.     

Constant fields and constant gradients in open ionic channels.

Ions enter cells through pores in proteins that are holes in dielectrics. The energy of interaction between ion and charge induced on the dielectric is many kT, and so the dielectric properties of channel and pore are important. We describe ionic movement by (three-dimensional) Nemst-Planck equations (including flux and net charge). Potential is described by Poisson's equation in the pore and Laplace's equation in the channel wall, allowing induced but not permanent charge. Asymptotic expansions are constructed exploiting the long narrow shape of the pore and the relatively high dielectric constant of the pore's contents. The resulting one-dimensional equations can be integrated numerically; they can be analyzed when channels are short or long (compared with the Debye length). Traditional constant field equations are derived if the induced charge is small, e.g., if the channel is short or if the total concentration gradient is zero. A constant gradient of concentration is derived if the channel is long. Plots directly comparable to experiments are given of current vs voltage, reversal potential vs. concentration, and slope conductance vs. concentration. This dielectric theory can easily be tested: its parameters can be determined by traditional constant field measurements. The dielectric theory then predicts current-voltage relations quite different from constant field, usually more linear, when gradients of total concentration are imposed. Numerical analysis shows that the interaction of ion and channel can be described by a mean potential if, but only if, the induced charge is negligible, that is to say, the electric field is spatially constant.     

Stationary-state distributions in membrane channels. I. Conservative flux systems.

The population distribution of cations at binding sites in membrane channels is determined for stationary states characterized by a net ion flux. For conservative systems, the net flux conserves the total intra-channel cation population for an ensemble of channels. The requisite matrix formalism is developed and illustrated for some homogeneous channels with N single-ion occupancy states with chemical and electrochemical transmembrane gradients. The lowest eigen-value approximation, which is used effectively for systems with no net cation flux, is applied to these stationary-state systems for comparison with the exact solutions. The stationary states are characterized by special thermodynamic functions for entropy, internal energy, and free energy defined by using information theory. The changes in these parameters for the transition from the equilibrium to the stationary state depend on the ratio of the next flux to the transition velocity for nearest-neighbor transitions within the channel. The free energy changes with the square of this ratio while the internal energy changes linearly.     

Propagation and amplification of microwave radiation in a plasma channel created in gas by a high-power femtosecond UV laser pulse

The time evolution of a nonequilibrium plasma channel created in a noble gas by a high-power femtosecond KrF laser pulse is investigated. It is shown that such a channel possesses specific electrodynamic properties and can be used as a waveguide for efficient transportation and amplification of microwave pulses. The propagation of microwave radiation in a plasma waveguide is analyzed by self-consistently solving (i) the Boltzmann kinetic equation for the electron energy distribution function at different spatial points and (ii) the wave equation in the parabolic approximation for a microwave pulse transported along the plasma channel.     

Application of the fluctuation theorem to motor proteins: from F<Subscript>1</Subscript>-ATPase to axonal cargo transport by kinesin and dynein

The fluctuation theorem is a representative theorem in non-equilibrium statistical physics actively studied in the 1990s. Relating to entropy production in non-equilibrium states, the theorem has been used to estimate the driving power of motor proteins from fluctuation in their motion. In this review, usage of the fluctuation theorem in experiments on motor proteins is illustrated for biologists, especially those who study mechanobiology, in which force measurement is a central issue. We first introduce the application of the fluctuation theorem in measuring the rotary torque of the rotary motor protein F₁-ATPase. Next, as an extension of this application, a recent trial estimating the force generated during cargo transport in vivo by the microtubule motors kinesin and dynein is introduced. Elucidation of the physical mechanism of such transport is important, especially for neurons, in which deficits in cargo transport are deeply related to neuronal diseases. Finally, perspectives on the fluctuation theorem as a new technique in the field of neuroscience are discussed.     

Perspectives on the landscape and flux theory for describing emergent behaviors of the biological systems

We give a review on the landscape theory of the equilibrium biological systems and landscape-flux theory of the nonequilibrium biological systems as the global driving force. The emergences of the behaviors, the associated thermodynamics in terms of the entropy and free energy and dynamics in terms of the rate and paths have been quantitatively demonstrated. The hierarchical organization structures have been discussed. The biological applications ranging from protein folding, biomolecular recognition, specificity, biomolecular evolution and design for equilibrium systems as well as cell cycle, differentiation and development, cancer, neural networks and brain function, and evolution for nonequilibrium systems, cross-scale studies of genome structural dynamics and experimental quantifications/verifications of the landscape and flux are illustrated. Together, this gives an overall global physical and quantitative picture in terms of the landscape and flux for the behaviors, dynamics and functions of biological systems.