Transport across homoporous and heteroporous membranes in nonideal, nondilute solutions. II. Inequality of phenomenological and tracer solute permeabilities.

The phenomenological solute permeability (omega p) of a membrane measures the flux of solute across it when the concentrations of the solutions on the two sides of the membrane differ. The relationship between omega p and the the conventionally measured tracer permeability (omega T) is examined for homoporous and heteroporous (parallel path) membranes in nonideal, nondilute solutions and in the presence of boundary layers. In general, omega p and omega T are not equal; therefore, predictions of transmembrane solute flux based on omega T are always subject to error. For a homoporous membrane, the two permeabilities become equal as the solutions become ideal and dilute. For heteroporous membranes, omega p is always greater than omega T. An upper bound on omega p- omega T is derived to provide an estimate of the maximum error in predicted solute flux. This bound is also used to show that the difference between omega P and omega T demonstrated earlier for the sucrose-Cuprophan system can be explained if the membrane is heteroporous. The expressions for omega P developed here support the use of a modified osmotic driving force to describe membrane transport in nonideal, nondilute solutions.     

Osmotic Transport across Cell Membranes in Nondilute Solutions: A New Nondilute Solute Transport Equation

The fundamental physical mechanisms of water and solute transport across cell membranes have long been studied in the field of cell membrane biophysics. Cryobiology is a discipline that requires an understanding of osmotic transport across cell membranes under nondilute solution conditions, yet many of the currently-used transport formalisms make limiting dilute solution assumptions. While dilute solution assumptions are often appropriate under physiological conditions, they are rarely appropriate in cryobiology. The first objective of this article is to review commonly-used transport equations, and the explicit and implicit assumptions made when using the two-parameter and the Kedem-Katchalsky formalisms. The second objective of this article is to describe a set of transport equations that do not make the previous dilute solution or near-equilibrium assumptions. Specifically, a new nondilute solute transport equation is presented. Such nondilute equations are applicable to many fields including cryobiology where dilute solution conditions are not often met. An illustrative example is provided. Utilizing suitable transport equations that fit for two permeability coefficients, fits were as good as with the previous three-parameter model (which includes the reflection coefficient, σ). There is less unexpected concentration dependence with the nondilute transport equations, suggesting that some of the unexpected concentration dependence of permeability is due to the use of inappropriate transport equations.     

A Biomechanical Triphasic Approach to the Transport of Nondilute Solutions in Articular Cartilage

Biomechanical models for biological tissues such as articular cartilage generally contain an ideal, dilute solution assumption. In this article, a biomechanical triphasic model of cartilage is described that includes nondilute treatment of concentrated solutions such as those applied in vitrification of biological tissues. The chemical potential equations of the triphasic model are modified and the transport equations are adjusted for the volume fraction and frictional coefficients of the solutes that are not negligible in such solutions. Four transport parameters, i.e., water permeability, solute permeability, diffusion coefficient of solute in solvent within the cartilage, and the cartilage stiffness modulus, are defined as four degrees of freedom for the model. Water and solute transport in cartilage were simulated using the model and predictions of average concentration increase and cartilage weight were fit to experimental data to obtain the values of the four transport parameters. As far as we know, this is the first study to formulate the solvent and solute transport equations of nondilute solutions in the cartilage matrix. It is shown that the values obtained for the transport parameters are within the ranges reported in the available literature, which confirms the proposed model approach.     

Dilute Solution Approximation and Generalization of the Reflection Coefficient Method of Describing Volume and Solute Flows

Kedem and Katchalsky introduced an approximation for dilute solutions which requires that the quantity (Δπ/Δπ_i)ø_i be much less than one. Zelman attempted to generalize the reflection coefficient concept to apply to solutions of multiple solutes, both penetrable and impenetrable, of concentrations sufficiently high for the approximation not to work. By simple algebraic manipulation, Zelman introduced a pair of new reflection coefficients, and a third new parameter γ which he misleadingly calls the “deviation from the dilute solution approximation.” It is shown here that the original Kedem-Katchalsky form for the flow equations can be preserved in such a way that no new coefficients need be introduced and an explicit statement of the effect of the dilute solution approximation can be made. There is an option of using a new set of conjugate driving forces for the solute flows or, alternatively, incorporating the nondilute solution correction in the coefficients in a clear way.     

Equations for membrane transport. Experimental and theoretical tests of the frictional model.

Frictional models for membrane transport are tested experimentally and theoretically for the simple case of a solution consisting of a mixture of two perfect gases and a membrane consisting of a porous graphite septum. Serious disagreement is found, which is traced to a missing viscous term. Kinetic theory is then used as a guide in formulating a corrected set of transport equations, and in giving a physical interpretation to the frictional coefficients. Sieving effects are found to be attributable to entrance effects rather than to true frictional effects within the body of the membrane. The results are shown to be compatible with nonequilibrium thermodynamics. Some correlations and predictions are made of the behavior of various transport coefficients for general solutions.     

Thermodynamic Model Equations for Heterogeneous Multicomponent Non-Ionic Solution Transport in a Multimembrane System

Non-equilibrium thermodynamic model equations for non-ionic and heterogeneous n-component solution transport in a m-membrane system are presented. This model is based on two equations. The first one describes the volume transport of the solution and the second the transport of the solute. Definitions of the hydraulic permeability, reflection and diffusive permeability coefficients of the m-membrane system and relations between the coefficients of the m-membrane system and the respective membranes of the system are also given. The validity of this model for binary and ternary solutions was verified, using a double-membrane cell with a horizontally mounted membrane. In the cell, volume and solute fluxes were measured as a function of concentration and gravitational configuration.     

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.     

The osmotic flow of an electrolyte through a charged porous membrane

A pore model in which the pore wall has a continuous distribution of electrical charge is used to investigate the osmotic flow through a charged permeable membrane separating electrolyte solutions of unequal concentrations. The pore is treated as a long, circular, cylindrical duct. The analysis is based on a continuum formulation in which a dilute electrolyte solution is described by the coupled Nernst-Planck/Poisson creeping flow equations. Account is taken of the significant size of the electrolyte ions (assumed to be rigid spheres) when compared with the diameter of the membrane pores. Analytical solutions for the ion concentrations, hydrostatic pressure and electrostatic potential in the electrolyte solutions are given and an intra-pore flow solution is derived. A mathematical expression for the osmotic reflection coefficient as a function of the solute ion: pore diameter ratio λ and the solute fluxes is obtained. Approximate solutions are quoted which relate the solute fluxes and the solution electrostatic potentials at the membrane surfaces to the bulk solution concentrations, the membrane pore charge and pore geometry. The osmotic reflection coefficient is thus determined as a function of these parameters.     

Reflection coefficients in pulsatile flow through converging junctions and the pressure distribution in a simple loop

Analytical expressions for the reflection coefficients in pulsatile flow through converging junctions are derived by two independent methods and are used to study the effects of wave reflections on the pressure distribution in a simple vascular loop. A simulated physiological situation is used as an example in which the loop is formed by the combination of a bypass and a bypassed vessel, the relative diameter of the latter being varied in order to simulate a narrowing. The results demonstrate how, in the case of a converging junction, the effects of wave reflections on the pressure distribution in one vessel depend on conditions within the vessel itself as well as in the other. The new reflection coefficients take into account this interdependence of flow in the two vessels forming a converging junction, and are shown to be consistent with reflection coefficients commonly used in diverging junctions.     

Existence and uniqueness of solutions in general multisolute renal flow problems

This paper considers systems of differential equations that describe flows in renal networks. The flow geometry is of the type that occurs in modelling the renal medulla. The unknowns in the system include the flow rate, the hydrostatic pressure, and the concentrations of the various solutes. Existence and uniqueness of solutions of the appropriate boundary value problems are established, in the case of small permeability coefficients and transport rates, or large diffusion coefficients and small resistance to flow constants.Work supported in part by NIH Grants 5-R01-AM28617 and 7-R01-DK38817Work supported in part by NIH Grant 5-R01-AM20373     

Membrane transport of electrolyte ions and time-dependent membrane potential

Equations are derived for the transport of a symmetrical electrolyte, consisting of cations and anions of equal valency, through a neutral membrane that separates two solutions of finite volume under quasi-steady-state conditions. The time-dependent membrane potential produced by the flow of ions is taken into account. Deviation of the time course of the solute concentrations from that of neutral solutes is found to be determined by the permeability ratio of cations and anions (when this ratio equals unity, the derived membrane transport equations reduce to those for neutral substances). Simple approximate expressions for the solute concentrations and of the membrane potential as functions of time are proposed, which are in excellent agreement with the exact numerical results.     

The effect of membrane heterogeneity on the predictability of fluxes,with application to the cornea

Many phenomenological treatments of biological membrane transport are based on the assumption that the membrane consists of a single class of passive transport paths; i.e. that the membrane is simple. Simplicity is assumed both in the measurement of the membrane's transport coefficients and in the use of these coefficients to predict membrane fluxes. Such a procedure will in general lead to an error in the prediction of solute flux across parallel arrays. The error depends on the distribution of the reflection coefficients of the parallel paths and upper and lower bounds on it are given in terms of conventionally measured transport coefficients. The frictional representation of electrolyte transport across membranes possessing metabolic pumps is generalized to take this structure effect into account. The flux error resulting from the neglect of membrane heteroreflectivity is essentially the same for nonelectrolytes and electrolytes, irrespective of whether phenomenological or frictional membrane transport properties are used. It is shown that some information about transport structure can be obtained from global measurements made without regard for the organization of pathways across the membrane.The values of the measured transport coefficients of the corneal epithelium and endothelium imply that these cell layers are heteroreflective. Analysis of corneal transport, taking structure effects into account, shows that the corneal thickness may be intrinsically insensitive to tear tonicity, by a mechanism which may be of more general homeostatic significance.     

Generalization of the Spiegler-Kedem-Katchalsky frictional model equations of the transmembrane transport for multicomponent non-electrolyte solutions.

The Spiegler-Kedem-Katchalsky frictional model equations of the transmembrane transport for systems containing n-component, non-ionic solutions is presented. The frictional interpretation of the phenomenological coefficients of membrane and the expressions connecting the practical coefficients (Lp, sigma i, omega ij) with frictional coefficients (fij) are presented.     

Electrokinetic membrane processes in relation to properties excitable tissues. II. Some theoretical considerations

A quantitative theory is presented for the behavior of a membrane-electrolyte system subject to an electric current flow (the "membrane oscillator"). If the membrane is porous, carries "fixed charges," and separates electrolyte solutions of different conductances, it can be the site of repetitive oscillatory changes in the membrane potential, the membrane resistance, and the hydrostatic pressure difference across the membrane. These events are accompanied by a pulsating transport of bulk solutions. The theory assumes the superposition of electrochemical and hydrostatic gradients and centers round the kinetics of resistance changes within the membrane, as caused by effects from diffusion and electro-osmotic fluid streaming. The results are laid down in a set of five simple, basic expressions, which can be transformed into a pair of non-linear differential equations yielding oscillatory solutions. A graphical integration method is also outlined (Appendix II). The agreement between the theory and previous experimental observations is satisfactory. The applied electrokinetic concepts may have importance in relation to analyses of the behavior of living excitable cells or tissues.     

Network thermodynamic analysis and stimulation of isotonic solute-coupled volume flow in leaky epithelia: an example of the use of network theory to provide the qualitative aspects of a complex system and its verification by stimulation

A network thermodynamic model was developed to provide insights into the nature of isotonic solute-coupled volume flow in "leaky" epithelia, where the transepithelial volume flow is assumed to be primarily through the cellular pathway. The coupled flows of solute and volume at each membrane in this four membrane model are described by the practical phenomenological equations as developed by Kedem & Katchalsky (1958). The model contains one permeable non-electrolyte solute (s) and a fixed amount of an impermeable non-electrolyte (i) inside the cell. The cell is assumed to be capable of volume regulation under the steady-state experimental conditions simulated. A solute-pump, located in the basolateral membrane, uses feedback regulation to adjust Cs in the cell in order to maintain cell volume at or near control levels in all simulations. Model behavior is, in general, very consistent with experimental observations with respect to tonicity and magnitude of volume flow over a wide range of experimental conditions. Examination of the parameter space suggests the following important features when isotonic solute-coupled volume flow moves primarily through the cellular pathway: (1) the apical membrane reflection coefficient must be less than that of the basolateral membrane; (2) the basement membrane reflection coefficient must be small; (3) the apical membrane solute permeability and reflection coefficient are the two most "sensitive" parameters and need to vary in an inverse manner in order to maintain isotonicity when both solute and volume flows increase; and (4) relationships (1) and (3) above imply the need for at least two separate solute pathways in the apical membrane, one that is shared with volume flow and one that is not.     

Global flow equations for membrane transport from local equations of motion—II. The case of a single nonelectrolyte solute plus water

The global flow equations of nonequilibrium thermodynamics for a single nonelectrolyte solute and water passing through a membrane are obtained by solving the local equations of motion. The method follows that developed for the general,n-solute case in the previous paper (Mikulecky, 1978). It is easily seen in this simple case that the passage from local interactions, formulated as position dependent frictional interactions in the equations of motion, to ghe global result involves a loss of any simple way of identifying particulars about local information. Two particular cases are analyzed in further detail: the case of no interaction within the pore and the case of constant interaction for both solute and solvent across the pore. In the former case, Onsager reciprocity survives in the global result if a self-consistent definition of the partial viscosity coefficients is used, while in the latter case, reciprocity is lost. Since, in many biologically interesting cases, the presence of interaction of the type considered here is likely to occur, the reciprocity condition should not automatically be assumed to hold.     

Convective paracellular solute flux. A source of ion-ion interaction in the epithelial transport equations

An electrolyte model of an epithelium (a cell and a tight junction in parallel, both in series with a lateral interspace basement membrane) is analyzed using the formalism of nonequilibrium thermodynamics. It is shown that if the parallel structures are heteroporous (i.e., reflection coefficients for two ion species differ between the components), then a cross-term will appear in the overall transport equations of the epithelium. Formally, this cross-term represents an ion-ion interaction. With respect to the rat proximal tubule, data indicating epithelial ionic reflection coefficients less than unity, together with the assumption of no transcellular solvent drag, imply the presence of convective paracellular solute flux. This means that a model applicable to a heteroporous structure must be used to represent the tubule, and, in particular, the cross-terms for ion-ion interaction must also be evaluated in permeability determinations. A series of calculations is presented that permits the estimation of the Na-Cl interaction for rat proximal tubule from available experimental data. One consequence of tubule heteroporosity is that an electrical potential may be substantially less effective than an equivalent concentration gradient in driving reabsorptive ion fluxes.