Rate theory models for ion transport through rigid pores. II. Time-dependent analysis of the single-file model with special reference to oscillatory behavior

This article deals with the time-dependent evolution of the single-file movement of ions through channels of both biological and artificial membranes. The single-file transport process may exhibit not only the usual relaxation behaviour but also oscillatory behaviour as a steady state is approached after an initial perturbation. A necessary condition for the occurrence of oscillations is that the system acts sufficiently far from equilibrium. The occurrence of oscillations is due to the interactions within the transport system which are taken into account by the single-file model; these are the electrostatic repulsion between the ions being transported, and the competition of the ions for the free binding sites within the pore. Information about the strength of the interactions can be obtained by measuring the damping of the transport observables (e.g. the electric current): The stronger the inter-ionic repulsion, the more apparent the oscillatory behaviour will become. Furthermore, the damping is influenced by the microscopic structure of the transport system (i.e. the energy profile of the pores). With an increasing degree of microscopicity, i.e. with a decreasing number of binding sites and an increasingly irregular pore profile, the oscillations become more damped. However, a considerable oscillatory behaviour can only be predicted for pores with both a sufficiently regular structure and a sufficiently large number of binding sites. For this class of pores, however, the measurement of the damping represents an appropriate method of gaining information which could exceed that obtainable from the usual methods of measuring stationary quantities (e.g. stationary conductance). Moreover, our goal is to explain theoretically how the oscillatory behaviour can be interpreted in terms of the order inherent in the ionic movement, which is determined by both the external and internal forces and the microscopic properties of the system.     

General continuum analysis of transport through pores. II. Nonuniform pores.

A general continuum derivation of the nonelectrolyte (Js) and volume (Jv) flux through a pore whose cross section is a function of axial position (nonuniform) is given. In general, the flux equations cannot be reduced to the same form as for a uniform pore and it is not possible to characterize the pore kinetics by three constants as in the uniform pore case. However, it is shown that under certain conditions, the nonuniform pore equations can be approximated by the uniform pore form and can be characterized by three constants (omega, sigma, Lp). The only condition needed to reduce the Jv equation to the uniform form is that the solution be dilute. The deviation of the Js equation from the uniform form is characterized by an asymmetrical function of Jv whose maximum value is estimated. It is shown that the maximum posible fractional deviation of the Js equation from the uniform form is given by the parameter: 0:5sigmaJv/omegaRT. Since this parameter is less then 0.15 for most membrane studies, the nonuniform Js equation can usually be approximated by the uniform pore form. The general results are illustrated by explicit calculations on several models of nonuniform pores. It is shown, for example, that the "equivalent pore radius" defined in the usual way is a function of the experimental parameter that is measured and is not unique.     

A hydrodynamic theory of ion conductance through Ohmic pores

A method of calculating the size of membrane pores lacking strong ionic selectivity is presented. By treating the flow of ions through a small channel as a hydrodynamic phenomenon, the electrical conductance becomes a function of the ratio of ion radius to channel radius. Thus when both the channel conductance and the ion size are known, the radius of the channel may be estimated. The method gives good agreement among radii predicted from conductances of four different alkali cations in alamethicin pores.     

General continuum analysis of transport through pores. I. Proof of Onsager's reciprocity postulate for uniform pore.

The nonelectrolyte (Js) and volume (Jv) flux across a membrane is usually described in terms of two equations derived from the theory of irreversible thermodynamics: (see article) where delta c and delta P are the concentration and pressure difference; omega and Lp are the diffuse and hydraulic permeability; and sigma s and sigma v are the reflection coefficients. If Onsager's reciprocity postulate is assumed, it can be shown that signa s and sigma v are equal. This is an important assumption because it allows one to apply the continuum theory relationship between sigma s and the pore radius to experimental measurements of sigma v. In this paper, general continuum expressions for both the Jv (a new result) and Js equation will be derived and the equality of sigma s and sigma v proved. The proof uses only general hydrodynamic results and does not require explicit solutions for the drag coefficients or, for example, the assumption that the solute is in the center of the pore. The proof applys to arbitrarily shaped solutes and any pore whose shape is independent of axial position (uniform). In addition, new expressions for the functional dependence of omega and sigma on the pore radius are derived (including the effect of the particle lying off the pore axis). These expressions differ slightly from earlier results and are probably more accurate.     

Two-carrier models for mediated transport. I. Theoretical analysis of several two-carrier models.

Several possible models of two sequential and two simultaneous carriers of different affinities are theoretically analysed. Following the analysis we suggest for each model an experimental procedure capable of testing and rejecting the model.     

Two-carrier models for mediated transport I. Theoretical analysis of several two-carrier models

Several possible models of two sequential and two simultaneous carriers of different affinities are theoretically analysed. Following the analysis we suggest for each model an experimental procedure capable of testing and rejecting the model.     

A unified approach to ion transport through membranes. I. Membrane-soluble ions.

Theoretical membrane potential transient produced by applying a current step to nerve cells has been derived based on the compartment neuron model and also on the equivalent cylinder model developed by W. Rall. It is expressed as a sum of exponential functions as

$\text{∑}\text{i=0}\text{n?1}\text{E}{}_{\mathrm{i}}\text{[1?}\text{exp}\text{(}\text{t}\text{τ}{}_{\mathrm{i}}\text{)]}$

where n is the number of compartments. The ratio of the amplitudes of the first and the second largest exponential functions, ($\text{E}{}_1\text{E}{}_0$), was found to be proportional to that of their respective time constants, ($\text{τ}{}_1\text{τ}{}_0$), in these neuron models. The constant of proportionality is given in a form that depends on the number of compartments as $\text{E}{}_1\text{E}{}_0\text{= (1 +}\text{cos}\text{π}\text{n}\text{)}\text{τ}{}_1\text{τ}{}_0$. This theoretical result is discussed in the light of recent experimental results in cat red nucleus neurons.     

Semiconductor theory of ion transport in thin lipid membranes: I. Potential and field distributions

In the theory as presented in this paper and the following one, we shall attempt to apply the semiconductor principles and methods to the study of ion transport in thin lipid membranes. Detailed formulations are given on the potential energy barriers at the interfaces, voltage drops in the polar and non-polar regions, and potential and field distributions in the diffuse double layer and within a charged membrane. These results will be used mainly as the boundary conditions for the solution of ion flow as to be given in the following paper. The analysis clearly indicates that the ion transport is interface-limited and is profoundly influenced by the presence of surface charges. An explanation of Na⁺ extrusion in nerve membrane is given based on the field distribution analysis. The theory also suggests that the “membrane potential” depends mainly on surface charges but not necessarily on ion permeation through the membrane.     

A lattice relaxation algorithm for three-dimensional Poisson-Nernst-Planck theory with application to ion transport through the gramicidin A channel.

A lattice relaxation algorithm is developed to solve the Poisson-Nernst-Planck (PNP) equations for ion transport through arbitrary three-dimensional volumes. Calculations of systems characterized by simple parallel plate and cylindrical pore geometries are presented in order to calibrate the accuracy of the method. A study of ion transport through gramicidin A dimer is carried out within this PNP framework. Good agreement with experimental measurements is obtained. Strengths and weaknesses of the PNP approach are discussed.     

Tests of continuum theories as models of ion channels. I. Poisson-Boltzmann theory versus Brownian dynamics

Continuum theories of electrolytes are widely used to describe physical processes in various biological systems. Although these are well-established theories in macroscopic situations, it is not clear from the outset that they should work in small systems whose dimensions are comparable to or smaller than the Debye length. Here, we test the validity of the mean-field approximation in Poisson-Boltzmann theory by comparing its predictions with those of Brownian dynamics simulations. For this purpose we use spherical and cylindrical boundaries and a catenary shape similar to that of the acetylcholine receptor channel. The interior region filled with electrolyte is assumed to have a high dielectric constant, and the exterior region representing protein a low one. Comparisons of the force on a test ion obtained with the two methods show that the shielding effect due to counterions is overestimated in Poisson-Boltzmann theory when the ion is within a Debye length of the boundary. As the ion gets closer to the boundary, the discrepancy in force grows rapidly. The implication for membrane channels, whose radii are typically smaller than the Debye length, is that Poisson-Boltzmann theory cannot be used to obtain reliable estimates of the electrostatic potential energy and force on an ion in the channel environment.     

Simple models for the analysis of binding protein-dependent transport systems.

Mathematical modeling was used to evaluate experimental data for bacterial binding protein-dependent transport systems. Two simple models were considered in which ligand-free periplasmic binding protein interacts with the membrane-bound components of transport. In one, this interaction was viewed as a competition with the ligand-bound binding protein, whereas in the other, it was considered to be a consequence of the complexes formed during the transport process itself. Two sets of kinetic parameters were derived for each model that fit the available experimental results for the maltose system. By contrast, a model that omitted the interaction of ligand-free binding protein did not fit the experimental data. Some applications of the successful models for the interpretation of existing mutant data are illustrated, as well as the possibilities of using mutant data to test the original models and sets of kinetic parameters. Practical suggestions are given for further experimental design.     

Statistical-mechanical theory of passive transport through semipermeable membranes.

The first general multicomponent equations for transport through semipermeable membranes are derived from basic statistical-mechanical principles. The procedure follows that used earlier for open membranes, but semipermeability is modelled mathematically by the introduction of external forces on the impermeant species. Gases are treated first in order to clarify the problems involved, but the final results apply to general nonideal solutions of any concentration. The mixed-solvent effect is treated rigorously, and a mixed-solvent osmotic pressure is defined. A useful specific identification of so-called osmotic flow is given, along with a demonstration that such an identification cannot be unique. Results are obtained both for discontinuous membrane models, and for a continuous model.     

Structure and dynamics of ion transport through gramicidin A.

Molecular dynamics calculations in which all atoms were allowed to move were performed on a water-filled ion channel of the polypeptide dimer gramicidin A (approximately 600 atoms total) in the head-to-head Urry model conformation. Comparisons were made among nine simulations in which four different ions (lithium, sodium, potassium, and cesium) were each placed at two different locations in the channel as well as a reference simulation with only water present. Each simulation lasted for 5 ps and was carried out at approximately 300 K. The structure and dynamics of the peptide and interior waters were found to depend strongly on the ion tested and upon its location along the pore. Speculations on the solution and diffusion of ions in gramicidin are offered based on the observations in our model that smaller ions tended to lie off axis and to distort the positions of the carbonyl oxygens more to achieve proper solvation and that the monomer-monomer junction was more distortable than the center of the monomer. With the potential energy surface used, the unique properties of the linear chain of interior water molecules were found to be important for optimal solvation of the various ions. Strongly correlated motions persisting over 25 A among the waters in the interior single-file column were observed.     

Low Reynolds number steady state flow through a branching network of rigid vessels: I. A mixture theory

A mixture theory has been used to formulate a theory of blood perfusion. By means of a formal averaging procedure the discrete network of microvessels is transformed into a continuum. During this procedure, the distinction between arterioles, capillaries and venules is preserved by means of an arteriovenous parameter. In this paper, two equations are derived for the case of low Reynolds number steady-state flow through a rigid vessel network: the extended Darcy equation and the continuity equation. A verification of the theory is presented, on the basis of a network analysis.