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.     

Theoretical study of flow,pressure and solute concentration in double membrane systems

A model system consisting of two rigidly held membranes in series was investigated through the application of the Kedem and Katchalsky thermodynamic single membrane flow equations. This analysis results in predictions of the steady state flow properties as well as values for the solute concentration and pressure of the internal compartment when the system is under the influence of a constant solute concentration or hydrostatic pressure gradient. It is demonstrated that although the flow properties and internal compartment pressure are complicated functions of the membrane permeability coefficients and driving gradient across the system, the relationships are greatly simplified by the explicit appearance of the internal compartment steady state solute concentration in the equations. It is shown that the steady state volume flow rate depends on the absolute value of the solute concentration in the external compartments, as well as the solute concentration gradient across the system. The properties of non-linear dependence of volume flow on concentration gradient, and rectification of volume flow are discussed and shown to be independent properties of the system. For the system under the influence of a solute concentration gradient, the internal compartment pressure can be greater or less than the ambient pressure, and depends mainly on the order in which the membranes are encountered by the volume flow. These properties are qualitatively correlated with certain available experimental observations in biological systems.     

Steady-state voltages, ion fluxes, and volume regulation in syncytial tissues.

Equations are developed that describe the steady-state relationships among ion fluxes, solute fluxes, water flow, voltage, concentration of solute, and hydrostatic pressure in a spherically symmetrical syncytial tissue. Each cell of the syncytium is assumed to have membrane channels for Na, K, and Cl, a membrane pump for Na/K, and some concentration of intracellular protein of net negative charge. However, the surface cells and inner cells of the tissue are assumed to have different distributions of membrane transport properties, hence there is a radial circulation of fluxes and a radial distribution of forces. Some reasonable approximations are made that allow analytic solutions of the nonlinear differential equations. These solutions are used to analyze data from the frog lens and are shown to account for the known steady-state properties of this tissue. Moreover, these solutions are used to make predictions on other steady-state properties, which have not been directly measured, and graphical results on the circulation of water, ions and solute through the frog lens are presented.     

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.     

A theoretical study on the sucrose gap technique as applied to multicellular muscle preparations. I. Saline-sucrose interdiffusion.

Voltage-clamp analysis of membrane currents in multicellular muscle preparations by means of the sucrose gap method is complicated by diffusion of saline and sucrose in the interstitial fluid spaces. This paper is the first part of a theoretical study made to analyze electrical events related to this diffusion process. Concentration profiles of ions and sucrose (both axial and radial) were computed by solving diffusion equations with boundary conditions appropriate for the different types of preparations and experimental arrangements used. In addition to steady-state solutions, analytical expressions were derived that describe the time-course with which concentration profiles become established after a stepwise change of the solute concentration in one of the compartments of the sucrose gap apparatus. The model accounts for the presence of an endothelial surface layer, or endocardium, which acts as an external diffusion barrier and is important in determining concentration gradients of solutes within heart cell preparations. Results of numerical computations dealing with several cases of experimental interest are presented.     

The Time Dependence of Single File Diffusion

The single file diffusion of particles through a narrow pore membrane separating two media is treated as a stochastic birth and death process. A set of differential-difference equations is derived to describe the probability of finding n particles in the pore at any time whose source is the left-hand medium. Explicit time-dependent solutions for an arbitary number of sites are obtained. These can be used to calculate both one-way and net flux as a function of time. Parameters are estimated from steady state permeability data, and the results of some numerical calculations are presented to illustrate the time required to approach a steady state. In many cases, significant time delays can occur.     

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.     

Dynamic loading of deformable porous media can induce active solute transport

Active solute transport mediated by molecular motors across porous membranes is a well-recognized mechanism for transport across the cell membrane. In contrast, active transport mediated by mechanical loading of porous media is a non-intuitive mechanism that has only been predicted recently from theory, but not yet observed experimentally. This study uses agarose hydrogel and dextran molecules as a model experimental system to explore this mechanism. Results show that dynamic loading can enhance the uptake of dextran by a factor greater than 15 over passive diffusion, for certain combinations of gel concentration and dextran molecular weight. Upon cessation of loading, the concentration reverts back to that achieved under passive diffusion. Thus, active solute transport in porous media can indeed be mediated by cyclical mechanical loading.     

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.     

Effectiveness Factors and Conversion in a Biocatalytic Membrane Reactor

Analytical expressions of the effectiveness factor of a biocatalytic membrane reactor, and its asymptote as the Thiele modulus becomes large, are presented. The evaluation of the effectiveness factor is based on the solution of the governing equations for solute transport in the two regions of the reactor, i.e. the lumen and the matrix (with the biofilm immobilized in the matrix). The lumen solution accounts for both axial diffusion and radial convective flow, while the matrix solution is based on Robin-type boundary conditions. The effectiveness factor is shown to be a function of the Thiele modulus, the partition coefficient, the Sherwood number, the Peclet number, and membrane thickness. Three regions of Thiele moduli are defined in the effectiveness factor graphs. These correspond with reaction rate limited, internal-diffusion limited, and external mass transfer limited solute transport. Radial convective flows were shown to only improve the effectiveness factor in the region of internal diffusion limitation. The assumption of first order kinetics is shown to be applicable only in the Thiele modulus regions of internal and external mass transfer limitation. An iteration scheme is also presented for estimating the effectiveness factor when the solute fractional conversion is known. The model is validated with experimental data from a membrane gradostat reactor immobilised with Phanerochaete chrysosporium for the production of lignin and manganese peroxidases. The developed model and experimental data allow for the determination of the Thiele modulus at which the effectiveness factor and fractional conversion are optimal.     

A Numerical Study of the Hydrodynamic Stable Concentration Boundary Layers in a Membrane System Under Microgravitational Conditions

On the basis of the classic formula of the concentration Rayleigh number and the Kedem–Katchalsky equation for diffusive membrane transport, we derived the equations of sixteenth order which show the dependence of the thicknesses of the concentration boundary layers on the difference of the solution concentrations, the concentration Rayleigh number, the solute permeability coefficient of the membrane and the diffusion coefficients in the solution, the kinematic viscosity of the solution, the density of solutions, the temperature and gravitational acceleration. The obtained equation has numerical solutions in the first, third and fourth quadrant of a co-ordinate system. However, only two solutions from the first quadrant of the co-ordinate system have physical meaning. Confining ourselves to the set of solutions with physical meaning only, the thicknesses of concentration boundary layers for different parameters occurring in the obtained equation were calculated numerically.     

Dynamic light scattering from collagen solutions. II. Photon correlation study of the depolarized light.

The depolarized forward-scattered light from solutions of rat tail collagen has been studied by photon correlation spectroscopy. The measured autocorrelation function is seen to decay on two widely different time scales. The decay time for the fast component is consistent with the rotational diffusion of rodlike collagen monomers. The slowly decaying autocorrelation component is attributed to large nonspecific aggregates of collagen. A substantial fraction of the collagen is in this aggregated form. Extrapolation of the faster decay times to zero concentration yields a value of θ = 1082 ± 30 sec^?1 for the rotational diffusion coefficient of the collagen monomer.     

Transport across homoporous and heteroporous membranes in nonideal, nondilute solutions. I. Inequality of reflection coefficients for volume flow and solute flow.

The Kirkwood formulation of the Stefan-Maxwell equations is used to develop the transport equations for a membrane bounded by nonideal, nondilute solutions. The reflection coefficients for volume flow and solute flow are not equal but are related by a simple expression that depends on the concentration of the bounding solutions. The ratio of the two coefficients is independent of heteroporous membrane structure and the thickness of adjacent boundary layers. Experimental measurements of these reflection coefficients for sucrose transport across Cuprophan verify this relationship; this indicates that the Onsager reciprocal relation, which is assumed by the theory, holds for nonideal, nondilute solutions. The two reflection coefficients may be made operationally identical by a simple redefination of the osmotic driving force.     

Electroosmosis in Membranes: Effects of Unstirred Layers and Transport Numbers: I. Theory

When a current is passed through a membrane system, differences in transport numbers between the membrane and the adjacent solutions will, in general, result in depletion and enhancement of concentrations at the membrane-solution interfaces. This will be balanced by diffusion back into the bulk solution, diffusion of solute back across the membrane itself, and osmosis resulting from these local concentration gradients. The two main results of such a phenomenon are (1) that there is a current-induced volume flow, which may be mistaken for electroosmosis, and (2) that there will generally develop transient changes in potential difference (PD) across membranes during and after the passage of current through them.     

Coupled water transport in standing gradient models of the lateral intercellular space.

A standing gradient model of the lateral intercellular space is presented which includes a basement membrane of finite solute permeability. The solution to the model equations is estimated analytically using the "isotonic convection approximation" of Segel. In the case of solute pumps uniformly distributed along the length of the channel, the achievement of isotonic transport depends only on the water permeability of the cell membranes. The ability of the model to transport water against an adverse osmotic gradient is the sum of two terms: The first term is simply that for a well-stirred compartment model and reflects basement membrane solute permeability. The second term measures the added strength due to diffusion limitation within the interspace. It is observed, however, that the ability for uphill water transport due to diffusion limitation is diminished by high cell membrane water permeability. For physiologically relevant parameters, it appears that the high water permeability required for isotonic transport renders the contribution of the standing gradient relatively ineffective in transport against an osmotic gradient. Finally, when the model transports both isotonically and against a gradient, it is shown that substantial intraepithelial solute polarization effects are unavoidable. Thus, the measured epithelial water permeability will grossly underestimate the water permeability of the cell membranes. The accuracy of the analytic approximation is demonstrated by numerical solution of the complete model equations.     

The development of osmotic flow through an unstirred layer

We investigate the errors involved in estimating the osmotic permeability of a semi-permeable membrane, from the measured osmotic flow and the difference in concentration of osmotically active solute across it, without taking account of the unstirred layer in the solution next to the membrane. In the problem solved, this layer is represented as a region of thickness δ at the far side of which a solute concentration C_b is imposed for time . The initial diffusion of solute towards the membrane causes the concentration at the membrane C_m to rise, generating an osmotic flow of water, J, whose convective effect opposes the diffusion. The problem is made non-linear by the dependence of J on C_m. Ultimately a steady state is set up, in which C_m is less than C_b. The solution is shown to depend on a single parameter β, equal to (L_pRT) δ C_b/D, where L_pRT is the osmotic permeability of the membrane and D is the diffusivity of solute. Solution of the steady state leads to a prediction of C_m/C_b as a function of β, and analysis of the decay of transient terms leads to a prediction of the decay time π, also as a function of β. Numerical data for membranes with a wide range of osmotic permeabilites, and for a reasonable range of solute, i.e. sucrose, concentrations, suggest that values of β can range from 0.001 or below to 7.5 or above. The former value implies negligible error in neglecting the unstirred layer, while the latter implies a 79%. error. For β = 0.1 and for δ = 2 × 10⁻⁴ m, π is predicted to be around 74 s. This decreases as β increases (for fixed δ); for values of β above about 27, the decay of transients is no longer monotonic but takes the form of damped oscillations.     

Passive Asymmetric Transport through Biological Membranes

The magnitude of passive diffusional solute transfer through artificial membranes is usually considered to be independent of the direction of the concentration gradient driving force. It can be shown, however, that a composite membrane, having as one component a membrane with a chemical reaction-facilitated diffusion transport mechanism, can result in an asymmetrical flux. An asymmetric flux caused by this type of structural heterogeneity may be one mechanism contributing to the asymmetric properties of biological membranes. Similar vectorial fluxes can be generated in interfacial solute transfer through membranes if hydrodynamic boundary layers occur at the membrane interface and reversible chemical reactions with the permeant species are involved in either phase.     

Ion and water fluxes in the ileum of rats

Studies have been carried out on the movement of salt and water across the small intestine of the rat. Segments of the ileum of anesthetized rats have been perfused in vivo with unbuffered NaCl solutions or isotonic solutions of NaCl and mannitol. Kinetic analysis of movements of Na²⁴ and Cl³⁶ has permitted determination of the efflux and influx of Na and Cl. Net water absorption has been measured using hemoglobin as a reference substance. Water was found to move freely in response to gradients of osmotic pressure. Net water flux from isotonic solutions with varying NaCl concentration was directly dependent on net solute flux. The amount of water absorbed was equivalent to the amount required to maintain the absorbed solute at isotonic concentration. These results have been interpreted as indicating that water movement is a passive process depending on gradients of water activity and on the rate of absorption of solute. The effluxes of Na and Cl are linear functions of concentration in the lumen, but both ions are actively transported by the ileum according to the criterion of Ussing (Acta Physiol. Scand., 1949, 19, 43). The electrical potential difference between the lumen and plasma has been interpreted as a diffusion potential slightly modified by the excess of active Cl flux over active Na flux. The physical properties of the epithelial membrane indicate that it is equivalent to a membrane having negatively charged uniform right circular pores of 36 Å radius occupying 0.001 per cent of the surface area.     

A multiscale theoretical model for diffusive mass transfer in cellular biological media

An integrated methodology is developed for the theoretical analysis of solute transport and reaction in cellular biological media, such as tissues, microbial flocs, and biofilms. First, the method of local spatial averaging with a weight function is used to establish the equation which describes solute conservation at the cellular biological medium scale, starting with a continuum-based formulation of solute transport at finer spatial scales. Second, an effective-medium model is developed for the self-consistent calculation of the local diffusion coefficient in the cellular biological medium, including the effects of the structural heterogeneity of the extra-cellular space and the reversible adsorption to extra-cellular polymers. The final expression for the local effective diffusion coefficient is: D(Abeta)=lambda(beta)D(Aupsilon), where D(Aupsilon) is the diffusion coefficient in water, and lambda(beta) is a function of the composition and fundamental geometric and physicochemical system properties, including the size of solute molecules, the size of extra-cellular polymer fibers, and the mass permeability of the cell membrane. Furthermore, the analysis sheds some light on the function of the extra-cellular hydrogel as a diffusive barrier to solute molecules approaching the cell membrane, and its implications on the transport of chemotherapeutic agents within a cellular biological medium. Finally, the model predicts the qualitative trend as well as the quantitative variability of a large number of published experimental data on the diffusion coefficient of oxygen in cell-entrapping gels, microbial flocs, biofilms, and mammalian tissues.