Osmotic flow equations for leaky porous membranes

A basic set of equations describing the flows of volume (Jv) and solute (Js) across a leaky porous membrane, coupled to the differences of osmotic and hydrostatic pressures d pi and dP has been derived by using general frictional theory. Denoting the mean pore concentration of solute by c*s and the hydraulic and diffusive conductances by Lp and Ps/RT the equations take the form Jv = LpdP + sigma sLp d pi Js = c*s(1 - sigma f)Jv + Ps d pi/RT sigma s = theta (1 - DsVs/DwVw - Ds/Dos) sigma f = 1 - theta DsVs/DwVw - Ds/Dos in which Dw and Ds are the diffusion coefficients for water and solute in the pore and Dos that for free solution. The relation between the reflection coefficients sigma s and sigma f for osmosis and ultrafiltration is then given by sigma s = sigma f - (1- theta)(1 - Ds/Dos), where theta is the diffusive-driven:pressure-driven flow ratio. These equations follow from the fact that in leaky pores osmosis occurs by diffusion alone and that there cannot be any Onsager symmetry leading to sigma s = sigma f. Symmetry holds in the limits where either the pore is small, when sigma s = sigma f = 1, or where the pore is large when sigma s = sigma f = 0.     

Osmotic flow in membrane pores of molecular size

Water transfer by osmosis through pores occurs either by viscous flow or diffusion depending on whether the driving osmolyte is able to enter the pore. Analysis of osmotic permeabilities (P _os )measured in antibiotic and cellular pore systems supports this distinction, showing that P _os approaches either the viscous value (P _f ) or the diffusive value (P _d )depending on the size of the osmolyte in relation to the pore radius. Macroscopic hydrodynamics and diffusion theory, when used with drag and steric coefficients within an appropriate osmotic model, apply with remarkable accuracy to channels of molecular dimensions where water molecules cannot pass each other, without the need to postulate any special flow regimes. It becomes apparent that the true viscous to diffusive flow ratio, P _f /P _d , can be separated from the effects of tracer filing by osmotic measurements alone. It does not monotonically decrease with the pore radius but rises steeply at the smaller radii which would apply to pores in cell membranes. Consequently, the application of the theory to osmotic and diffusive flow data for the red cell predicts a pore radius of 0.2 nm in agreement with other recent measurements on isolated components of the system, showing that the viscous-diffusive distinction applies even in molecular pores.     

Mechanism of osmotic flow in a periodic fiber array

The classic analysis by Anderson and Malone (Biophys J 14: 957-982, 1974) of the osmotic flow across membranes with long circular cylindrical pores is extended to a fiber matrix layer wherein the confining boundaries are the fibers themselves. The equivalent of the well-known result for the reflection coefficient sigma0 = (1 - phi)2, where phi is the partition coefficient, is derived for a periodic fiber array of hexagonally ordered core proteins. The boundary value problem for the potential energy function describing the solute distribution surrounding each fiber is solved by defining an equivalent fluid annulus in which the pressures and osmotic forces are determined. This model is of special interest in the osmotic flow of water across a capillary wall, where recent experimental studies suggest that the endothelial glycocalyx is a quasiperiodic fiber array that serves as the primary molecular sieve for plasma proteins. Results for the reflection coefficient are presented in terms of two dimensionless numbers, alpha = a/R and beta = b/R, where a and b are the solute and fiber radii, respectively, and R is the outer radius of the fluid annulus. In general, the results differ substantially from the classic expression for a circular pore because of the large difference in the shape of the boundary along which the osmotic force is generated. However, as in circular pore theory, one finds that the reflection coefficients for osmosis and filtration are the same.     

Osmotic Flow of Water across Permeable Cellulose Membranes

Direct measurements have been made of the net volume flow through cellulose membranes, due to a difference in concentration of solute across the membrane. The aqueous solutions used included solutes ranging in size from deuterated water to bovine serum albumin. For the semipermeable membrane (impermeable to the solute) the volume flow produced by the osmotic gradient is equal to the flow produced by the hydrostatic pressure RT ΔC, as given by the van't Hoff relationship. In the case in which the membrane is permeable to the solute, the net volume flow is reduced, as predicted by the theory of Staverman, based on the thermodynamics of the steady state. A means of establishing the amount of this reduction is given, depending on the size of the solute molecule and the effective pore radius of the membrane. With the help of these results, a hypothetical biological membrane moving water by osmotic and hydrostatic pressure gradients is discussed.     

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.     

Single Water Channels of Aquaporin-1 Do not Obey the Kedem-Katchalsky Equations

The Kedem-Katchalsky (KK) equations are often used to obtain information about the osmotic properties and conductance of channels to water. Using human red cell membranes, in which the osmotic flow is dominated by Aquaporin-1, we show here that compared to NaCl the reflexion coefficient of the channel for methylurea, when corrected for solute volume exchange and for the water permeability of the lipid membrane, is 0.54. The channels are impermeable to these two solutes which would seem to rule out flow interaction and require a reflexion coefficient close to 1.0 for both. Thus, two solutes can give very different osmotic flow rates through a semi-permeable pore, a result at variance with both classical theory and the KK formulation. The use of KK equations to analyze osmotic volume changes, which results in a single hybrid reflexion coefficient for each solute, may explain the discrepancy in the literature between such results and those where the equations have not been employed. Osmotic reflexion coefficients substantially different from 1.0 cannot be ascribed to the participation of other 'hidden' parallel aqueous channels consistently with known properties of the membrane. Furthermore, we show that this difference cannot be due to second-order effects, such as a solute-specific interaction with water in only part of the channel, because the osmosis is linear with driving force down to zero solute concentration, a finding which also rules out the involvement of unstirred-layer effects. Reflexion coefficients smaller than 1.0 do not necessitate water-solute flow interaction in permeable aqueous channels; rather, the osmotic behaviour of impermeable molecular-sized pores can be explained by differences in the fundamental nature of water flow in regions either accessible or inaccessible to solute, created by a varying cross-section of the channel.     

Osmotic flow caused by polyelectrolytes

Osmotic flow of water caused by high concentrations of anionic polyelectrolytes across semipermeable membranes, permeable only to solvent and simple electrolyte, has been measured in a newly designed flow cell. The flow cell features small solution and solvent compartments and an efficient stirring mechanism. We have demonstrated that, while the osmotic pressure of the anionic polyelectrolytes is determined primarily by micro-counterions, the osmotic flow is determined by solution-dependent properties as embodied in the hydrodynamic frictional coefficient which is determined by the polymer backbone segment of the polyelectrolyte. The variation of the osmotic permeability coefficient, L(p)(o), with concentration and osmotic pressure closely correlated with the concentration dependence of this frictional coefficient. These studies confirm previous work that the kinetics of osmotic flow across a membrane impermeable to the osmotically active solute is primarily determined by the diffusive mobility of the solute.     

Permeability studies on red cell membranes of dog, cat, and beef

Water permeability coefficients of dog, cat, and beef red cell membranes have been measured under an osmotic pressure gradient. The measurements employed a rapid reaction stop flow apparatus with which cell shrinking was measured under a relative osmotic pressure gradient of 1.25 to 1.64 times the isosmolar concentration. For the dog red cell the osmotic permeability coefficient is 0.36 cm⁴/(sec, osmol). The water permeability coefficient for the dog red cell under a diffusion gradient was also measured (rate constant = 0.10/msec). The ratio between the two permeabilities was used to calculate an equivalent pore radius of 5.9 A. This value agrees welt with an equivalent pore radius of 6.2 A obtained from reflection coefficients of nonelectrolyte water-soluble molecules, and is consistent with data on the permeability of the dog red cell membrane to glucose. These data provide evidence supporting the existence of equivalent pores in single biological membranes.     

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.     

Fluid mechanical and electrostatic study on the osmotic flow through cylindrical pores

When two solutions differing in solute concentration are separated by a porous membrane, the osmotic pressure will generate a net volume flux of the suspending fluid across the membrane; this is termed osmotic flow. We consider the osmotic flow across a membrane with circular cylindrical pores when the solute and the pore walls are electrically charged, and the suspending fluid is an electrolytic solution containing small cations and anions. Under the condition in which the radius of the pores and that of the solute molecules greatly exceed those of the solvent as well as the ions, a fluid mechanical and electrostatic theory is introduced to describe the osmotic flow in the presence of electric charge. The interaction energy, including the electrostatic interaction between the solute and the pore wall, plays a key role in determining the osmotic flow. We examine the electrostatic effect on the osmotic flow and discuss the difference in the interaction energy determined from the nonlinear Poisson-Boltzmann equation and from its linearized equation (the Debye-Hückel equation).     

A transmural pressure gradient induces mechanical and biological adaptive responses in endothelial cells

A sudden increase in the transmural pressure gradient across endothelial monolayers reduces hydraulic conductivity (L(p)), a phenomenon known as the sealing effect. To further characterize this endothelial adaptive response, we measured bovine aortic endothelial cell (BAEC) permeability to albumin and 70-kDa dextran, L(p), and the solvent-drag reflection coefficients (sigma) during the sealing process. The diffusional permeability coefficients for albumin (1.33 +/- 0.18 x 10(-6) cm/s) and dextran (0.60 +/- 0.16 x 10(-6) cm/s) were measured before pressure application. The effective permeabilities (measured when solvent drag contributes to solute transport) of albumin and dextran (P(ealb) and P(edex)) were measured after the application of a 10 cmH(2)O pressure gradient; during the first 2 h of pressure application, P(ealb), P(edex), and L(p) were significantly reduced by 2.0 +/- 0.3-, 2.1 +/- 0.3-, and 3.7 +/- 0.3-fold, respectively. Immunostaining of the tight junction (TJ) protein zonula occludens-1 (ZO-1) was significantly increased at cell-cell contacts after the application of transmural pressure. Cytochalasin D treatment significantly elevated transport but did not inhibit the adaptive response, whereas colchicine treatment had no effect on diffusive permeability but inhibited the adaptive response. Neither cytoskeletal inhibitor altered sigma despite significantly elevating both L(p) and effective permeability. Our data suggest that BAECs actively adapt to elevated transmural pressure by mobilizing ZO-1 to intercellular junctions via microtubules. A mechanical (passive) component of the sealing effect appears to reduce the size of a small pore system that allows the transport of water but not dextran or albumin. Furthermore, the structures of the TJ determine transport rates but do not define the selectivity of the monolayer to solutes (sigma).     

Osmosis and solute—Solvent drag

In 1903, George Hulett explained how solute alters water in an aqueous solution to lower the vapor pressure of its water. Hulett also explained how the same altered water causes osmosis and osmotic pressure when the solution is separated from liquid water by a membrane permeable to the water only. Hulett recognized that the solute molecules diffuse toward all boundaries of the solution containing the solute. Solute diffusion is stopped at all boundaries, at an open-unopposed surface of the solution, at a semipermeable membrane, at a container wall, or at the boundary of a solid or gaseous inclusion surrounded by solution but not dissolved in it. At each boundary of the solution, the solute molecules are reflected, they change momentum, and the change of momentum of all reflected molecules is a pressure, a solute pressure (i.e., a force on a unit area of reflecting boundary). When a boundary of the solution is open and unopposed, the solute pressure alters the internal tension in the force bonding the water in its liquid phase, namely, the hydrogen bond. All altered properties of the water in the solution are explained by the altered internal tension of the water in the solution. We acclaim Hulett's explanation of osmosis, osmotic pressure, and lowering of the vapor pressure of water in an aqueous solution. His explanation is self-evident. It is the necessary, sufficient, and inescapable explanation of all altered properties of the water in the solution relative to the same property of pure liquid water at the same externally applied pressure and the same temperature. We extend Hulett's explanation of osmosis to include the osmotic effects of solute diffusing through solvent and dragging on the solvent through which it diffuses. Therein lies the explanations of (1) the extravasation from and return of interstitial fluid to capillaries, (2) the return of luminal fluid in the proximal and distal convoluted tubules of a kidney nephron to their peritubular capillaries, (3) the return of interstitial fluid to the vasa recta, (4) return of aqueous humor to the episcleral veins, and (5) flow of phloem from source to sink in higher plants and many more examples of fluid transport and fluid exchange in animal and plant physiology. When a membrane is permeable to water only and when it separates differing aqueous solutions, the flow of water is from the solution with the lower osmotic pressure to the solution with the higher osmotic pressure.     

Measurement of grouped intracellular solute osmotic virial coefficients

Models of cellular osmotic behaviour depend on thermodynamic solution theories to calculate chemical potentials in the solutions inside and outside the cell. These solutions are generally thermodynamically non-ideal under cryobiological conditions. The molality-based Elliott et al. form of the multi-solute osmotic virial equation is a solution theory which has been demonstrated to provide accurate predictions in cryobiological solutions, accounting for the non-ideality of these solutions using solute-specific thermodynamic parameters called osmotic virial coefficients. However, this solution theory requires as inputs the exact concentration of every solute in the solution being modeled, which poses a problem for the cytoplasm, where such detailed information is rarely available. This problem can be overcome by using a grouped solute approach for modeling the cytoplasm, where all the non-permeating intracellular solutes are treated as a single non-permeating “grouped” intracellular solute. We have recently shown (Zielinski et al., J Physical Chemistry B, 2017) that such a grouped solute approach is theoretically valid when used with the Elliott et al. model, and Ross-Rodriguez et al. (Biopreservation and Biobanking, 2012) have previously developed a method for measuring the cell type-specific osmotic virial coefficients of the grouped intracellular solute. However, the Ross-Rodriguez et al. method suffers from a lack of precision, which—as we demonstrate in this work—can severely impact the accuracy of osmotic model predictions under certain conditions. Thus, we herein develop a novel method for measuring grouped intracellular solute osmotic virial coefficients which yields more precise values than the existing method and then apply this new method to measure these coefficients for human umbilical vein endothelial cells.     

The Kinetics of Osmotic Transport through Pores of Molecular Dimensions

This paper presents a theoretical analysis of the kinetics of osmotic transport across a semipermeable membrane. There is a thermodynamic connection between the rate of flow under a hydrostatic pressure difference and the rate of flow due to a difference in solute concentration on the two sides. One might therefore attempt to calculate the osmotic transport coefficient by applying Poiseuille's equation to the flow produced by a difference in hydrostatic pressure. Such a procedure is, however, inappropriate if the pores in the membrane are too small to allow molecules to “overtake.” It then becomes necessary to perform a statistical calculation of the transport coefficient, and such a calculation is described in this paper. The resulting expression for the number of solvent molecules passing through a pore per second is J = m D₁ δn₁/l² where m is the number of solvent molecules in the pore, l is the length of the pore, D₁ is the self-diffusion coefficient of the solute, and δn₁ the difference in solvent mole fraction on the two sides of the membrane. This equation is used for estimating the number of pores per unit area of the squid axon membrane; the result is 6 × 10⁹ pores/cm².     

Some compartmental models of the root: steady-state behavior

Plots of the pressure difference (DeltaP) applied to plant roots vs. the resulting volume flow rate (Q(v)) often exhibit an anomalous offset that has been difficult to explain. The present analysis suggests that solute build-up in two- and three-compartment models of the root cannot account for this offset. The Ginsburg-Newman three-compartment model explains the offset in terms of differing reflection coefficients for the membranes bounding the intermediate compartment. This model appears more promising, but it predicts a minimum in the plot of xylem-sap osmotic pressure vs. Q(v)which is not observed in practice. Fiscus hypothesized that an internal asymmetric distribution of non-mobile solutes is responsible for the offset. In the present study, this hypothesis is incorporated into a four-compartment model of the root that is conceptually related to the three-compartment model of Miller. But according to the four-compartment model, the asymmetric solute distribution does not arise because of solvent drag. Rather the anomalous offset is associated with a concentration gradient of photoassimilates (the non-mobile solutes) that exists in the absence of volume flow, and which drives the diffusive transport of these solutes from the stele to the cortex via endodermal plasmodesmata. This model is consistent with the existence of radial symplastic osmotic-pressure gradients, and it appears to have greater explanatory power than the Ginsburg-Newman model. In particular, it suggests explanations for diurnal variations in DeltaP-Q(v)curves, as well as the effects of changing external solute concentrations. It also shows how the overall root reflection coefficient can be less than unity, even when the cell membranes are effectively ideally semipermeable, and there is negligible extracellular transport of water and solutes. The model makes a number of experimentally testable predictions.     

Reflection coefficients of homopore membranes: effect of molecular size and configuration

Osmotic water flow through membranes with uniform defined pores was measured for a variety of macromolecular solutes. Water flow increased linearly with applied hydrostatic pressure, allowing the effective osmotic pressure of the solutes to be estimated by extrapolation. Reflection coefficients for each solute-membrane combination were calculated and correlated with the ratio of solute size to pore size. For the same mean molecular size, proteins were found to have larger reflection coefficients than dextrans. Molecular rigidity may play a role in this difference in behavior.     

NaCl reflection coefficients in proximal tubule apical and basolateral membrane vesicles. Measurement by induced osmosis and solvent drag.

Two independent methods, induced osmosis and solvent drag, were used to determine the reflection coefficients for NaCl (sigma NaCl) in brush border and basolateral membrane vesicles isolated from rabbit proximal tubule. In the induced osmosis method, vesicles loaded with sucrose were subjected to varying inward NaCl gradients in a stopped-flow apparatus. sigma NaCl was determined from the osmolality of the NaCl solution required to cause no initial osmotic water flux as measured by light scattering (null point). By this method sigma NaCl was greater than 0.92 for both apical and basolateral membranes with best estimates of 1.0. sigma NaCl was determined by the solvent drag method using the Cl-sensitive fluorescent indicator, 6-methoxy-N-[3-sulfopropyl]quinolinium (SPQ), to detect the drag of Cl into vesicles by inward osmotic water movement caused by an outward osmotic gradient. sigma NaCl was determined by comparing experimental data with theoretical curves generated using the coupled flux equations of Kedem and Katchalsky. By this method we found that sigma NaCl was greater than 0.96 for apical and greater than 0.98 for basolateral membrane vesicles, with best estimates of 1.0 for both membranes. These results demonstrate that sigma NaCl for proximal tubule apical and basolateral membranes are near unity. Taken together with previous results, these data suggest that proximal tubule water channels are long narrow pores that exclude NaCl.     

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.