A model for the modulation of microvessel permeability by junction strands

To investigate the effect of junction strands on microvessel permeability, we extend the previous analytical model developed by Fu et al. (1994, J. Biomech. Eng., 116, pp. 502-513), for the interendothelial cleft to include multiple junction strands in the cleft and an interface between the surface glycocalyx layer and the cleft entrance. Based on the electron microscopic observations by Adamson et al. (1998, Am. J. Physiol., 274(43), pp. H1885-H1894), that elevation of intracellular cAMP levels would increase number of tight junction strands, this two-junction-strand and two-pore model can successfully account for the experimental data for the decreased permeability to water, small and intermediate-sized solutes by cAMP.     

An electrodiffusion-filtration model for effects of endothelial surface glycocalyx on microvessel permeability to macromolecules

Endothelial surface glycocalyx plays an important role in the regulation of microvessel permeability by possibly changing its charge and configuration. To investigate the mechanisms by which surface properties of the endothelial cells control the changes in microvessel permeability, we extended the electrodiffusion model developed by Fu et al. [Am. J. Physiol. 284, H1240-1250 (2003)], which is for the interendothelial cleft with a negatively charged surface glycocalyx layer, to include the filtration due to hydrostatic and oncotic pressures across the microvessel wall as well as the electrical potential across the glycocalyx layer On the basis of the hypotheses proposed by Curry [Microcirculation 1(1): 11-26 (1994)], the predictions from this electrodiffusion-filtration model provide a good agreement with experimental data for permeability of negatively charged a-lactalbumin summarized in Curry [Microcirculation 1(1), 11-26 (1994)] under various conditions. In addition, we applied this new model to describe the transport of negatively charged macromolecules, bovine serum albumin (BSA), across venular microvessels in frog mesentery. According to the model, the convective component of the albumin transport is greatly diminished by the presence of a negatively charged glycocalyx under both normal and increased permeability conditions.     

An electrodiffusion model for effects of surface glycocalyx layer on microvessel permeability

To investigate the charge effect of the endothelial surface glycocalyx on microvessel permeability, we extended the three-dimensional model developed by Fu et al. (J Biomech Eng 116: 502-513, 1994) for the interendothelial cleft to include a negatively charged glycocalyx layer at the entrance of the cleft. Both electrostatic and steric exclusions on charged solutes were considered within the glycocalyx layer and at the interfaces. Four charge-density profiles were assumed for the glycocalyx layer. Our model indicates that the overall solute permeability across the microvessel wall including the surface glycocalyx layer and the cleft region is independent of the charge-density profiles as long as they have the same maximum value and the same total charge. On the basis of experimental data, this model predicts that the charge density would be 25-35 meq/l in the glycolcalyx of frog mesenteric capillaries. An intriguing prediction of this model is that when the concentrations of cations and anions are unequal in the lumen due to the presence of negatively charged proteins, the negatively charged glycocalyx would provide more resistance to positively charged solutes than to negatively charged ones.     

Model for the chemotactic response of a bacterial population.

We present a mathematical model for the motion of a bacterial population in prescribed attractant or repellent gradients. The model is suggested by the observations of Mesibov et al. (1973, J. Gen. Physiol. 62:203) and Brown and Berg (1974, Proc. Natl. Acad. Sci. U.S.A. 71:1388) who found that the sensitivity of the chemotactic response depends on the concentration of attractant. Predictions of the theory are in general agreement with the experiments of Dahlquist et al. (1972, Nat. New Biol. 236:120) and of Mesibov et al. on populations of motile bacteria in fixed attractant gradients. Additional tests of the model are proposed.     

Effects of heterogeneities on the partitioning of airway and tissue properties in normal mice.

We measured the mechanical properties of the respiratory system of C57BL/6 mice using the optimal ventilation waveform method in closed- and open-chest conditions at different positive end-expiratory pressures. The tissue damping (G), tissue elastance (H), airway resistance (Raw), and hysteresivity were obtained by fitting the impedance data to three different models: a constant-phase model by Hantos et al. (Hantos Z, Daroczy B, Suki B, Nagy S, Fredberg JJ. J Appl Physiol 72: 168-178, 1992), a heterogeneous Raw model by Suki et al. (Suki B, Yuan H, Zhang Q, Lutchen KR. J Appl Physiol 82: 1349-1359, 1997), and a heterogeneous H model by Ito et al. (Ito S, Ingenito EP, Arold SP, Parameswaran H, Tgavalekos NT, Lutchen KR, Suki B. J Appl Physiol 97: 204-212, 2004). Both in the closed- and open-chest conditions, G and hysteresivity were the lowest and Raw the highest in the heterogeneous Raw model, and G and H were the largest in the heterogeneous H model. Values of G, Raw, and hysteresivity were significantly higher in the closed-chest than in the open-chest condition. However, H was not affected by the conditions. When the tidal volume of the optimal ventilation waveform was decreased from 8 to 4 ml/kg in the closed-chest condition, G and hysteresivity significantly increased, but there were smaller changes in H or Raw. In summary, values of the obtained mechanical properties varied among these models, primarily due to heterogeneity. Moreover, the mechanical parameters were significantly affected by the chest wall and tidal volume in mice. Contribution of the chest wall and heterogeneity to the mechanical properties should be carefully considered in physiological studies in which partitioning of airway and tissue properties are attempted.     

Starling forces that oppose filtration after tissue oncotic pressure is increased

We tested the hypothesis that the effective oncotic force that opposes fluid filtration across the microvessel wall is the local oncotic pressure difference across the endothelial surface glycocalyx and not the global difference between the plasma and tissue. In single frog mesenteric microvessels perfused and superfused with solutions containing 50 mg/ml albumin, the effective oncotic pressure exerted across the microvessel wall was not significantly different from that measured when the perfusate alone contained albumin at 50 mg/ml. Measurements were made during transient and steady-state filtration at capillary pressures between 10 and 35 cmH(2)O. A cellular-level model of coupled water and solute flows in the interendothelial cleft showed water flux through small breaks in the junctional strand limited back diffusion of albumin into the protected space on the tissue side of the glycocalyx. Thus oncotic forces opposing filtration are larger than those estimated from blood-to-tissue protein concentration differences, and transcapillary fluid flux is smaller than estimated from global differences in oncotic and hydrostatic pressures.     

Calculation of the reflection coefficient from measurements of endogenous vascular indicators

The solvent drag reflection coefficient (sigma) for total proteins can be estimated by comparing the relative degrees of concentration of erythrocytes and plasma proteins that occur during fluid filtration in an isolated perfused organ. In this analysis, we evaluated the accuracy of equations proposed by Pilati and Maron [Am. J. Physiol. 247 (Heart Circ. Physiol. 16): H1-H7, 1984] and Wolf et al. [Am. J. Physiol. 253 (Heart Circ. Physiol. 22): H194-H204, 1987] to calculate sigma from these concentration changes. We calculated sigma with each equation using data generated from a mathematical model of fluid and solute flux in membranes with known sigma's. We found that the equation of Wolf et al. provided the closest approximation to the true sigma over the entire range of filtration fractions tested (0.1-0.6), with the differences between the two equations increasing with filtration fraction. At low filtration fractions, the difference in sigma obtained using either approach was found to be inconsequential. At larger filtration fractions, a closer approximation of the true sigma can be obtained using the equation of Wolf et al.     

Microstructural analysis of deformation-induced hypoxic damage in skeletal muscle

Deep pressure ulcers are caused by sustained mechanical loading and involve skeletal muscle tissue injury. The exact underlying mechanisms are unclear, and the prevalence is high. Our hypothesis is that the aetiology is dominated by cellular deformation (Bouten et al. in Ann Biomed Eng 29:153-163, 2001; Breuls et al. in Ann Biomed Eng 31:1357-1364, 2003; Stekelenburg et al. in J App Physiol 100(6):1946-1954, 2006) and deformation-induced ischaemia. The experimental observation that mechanical compression induced a pattern of interspersed healthy and dead cells in skeletal muscle (Stekelenburg et al. in J App Physiol 100(6):1946-1954, 2006) strongly suggests to take into account the muscle microstructure in studying damage development. The present paper describes a computational model for deformation-induced hypoxic damage in skeletal muscle tissue. Dead cells stop consuming oxygen and are assumed to decrease in stiffness due to loss of structure. The questions addressed are if these two consequences of cell death influence the development of cell injury in the remaining cells. The results show that weakening of dead cells indeed affects the damage accumulation in other cells. Further, the fact that cells stop consuming oxygen after they have died, delays cell death of other cells.     

20-hydroxyeicosatetraenoic acid (20-HETE): structural determinants for renal vasoconstriction

The effects of natural and synthetic eicosanoids on the diameter of rat interlobular arteries studied in vitro were compared to that of the potent, endogenous vasoconstrictor 20-HETE. Vasoconstrictor activity was optimum for chain lengths of 20-22 carbons with at least one olefin or epoxide between located between C(13)-C(15) and an oxygen substituent at C(20)-C(22). The presence of delta (Zou et al. Am. J. Physiol. 1996, 270, R228; Gebremedhin, D. et al. Am. J. Physiol. 1998, 507, 771)-, delta (Carroll et al. Am. J. Physiol. 1996, 271, R863; Vazquez et al. Life Sci. 1995, 56, 1455)-, or delta (Imig et al. Hypertension 2000, 35, 307; Lopez et al. Amer. J. Physiol. 2001, 281, F420)-olefins had no influence on the vasoconstrictor response whereas the introduction of a C(7)-thiomethylene enhanced potency. A sulfonamide or alcohol, but not a lactone, could replace the C(1)-carboxylate. These data were used to construct a putative binding domain map of the 20-HETE receptor consisting of: (i) a comparatively open, hydrophilic binding site accommodating the C(1)-functionality; (ii) a hydrophobic trough spanning the olefins; (iii) a shallow pocket containing a critical pi-pi binding site in the vicinity of the pi (Ito et al. Am. J. Physiol. 1998, 274, F395; Quigley, R.; Baum, M.; Reddy, K. M.; Griener, J. C.; Falck, J. R. Am. J. Physiol. 2000, 278, F949)-olefin; and (iv) an oxyphilic binding site proximate to the omega-terminus.     

A one-dimensional model for the propagation of transient pressure waves through the lung

The propagation of pressure waves in the lung has been investigated by many authors concerned with respiratory physiology, ultrasound medical techniques or thoracic impact injuries. In most of the theoretical studies, the lung has been modeled as an isotropic and homogeneous medium, and by using Hooke's constitutive law (see e.g. Ganesan et al. Respir. Physiol. 110 (1997) 19; Jahed et al. J. Appl. Physiol. 66 (1989) 2675; Grimal et al. C.R. Acad. Sci., Paris 329 (IIb) (2001) 655-662), or more elaborated material laws (see, e.g. Bush and Challener (Proceedings of the International Research Council on Biokinetics Impacts (IRCOBI), Bergish-gladbach, 1988); Stuhmiller et al. J. Trauma 28 (1988) S132; Yang and Wang, Finite element modeling of the human thorax. Web page: http://wwwils.nlm.nih.gov/research/visible/vhpconf98/AUTHORS/YANG/YANG.HTM.). The hypothesis of homogeneous medium may be inappropriate for certain problems. Because of its foam-like structure, the behavior of the lung-even if the air and the soft tissue are assumed to behave like linearly elastic materials-is susceptible to be frequency dependent. In the present study, the lung is viewed as a one-dimensional stack of air and soft tissue layers; wave propagation in such a stack can be investigated in an equivalent mass-spring chain (El-Raheb (J. Acoust. Soc. Am. 94 (1993) 172; Int. J. Solids Struct. 34 (1997) 2969), where the masses and springs, respectively, represent the alveolar walls and alveolar gas. Results are presented in the time and frequency domains. The frequency dependence (cutoff frequency, variations in phase velocity) of the lung model is found to be highly dependent on the mean alveolar size. We found that short pulses induced by high velocity impacts (bullet stopped by a bulletproof jacket) can be highly distorted during the propagation. The pressure differential between two alveoli is discussed as a possible injury criterion.     

The epithelial Na(+) channel (ENaC) is related to the hypertonicity-induced Na(+) conductance in rat hepatocytes

The epithelial Na(+) channel (ENaC) is composed of the subunits alpha, beta, and gamma [Canessa et al., Nature 367 (1994) 463-467] and typically exhibits a high affinity to amiloride [Canessa et al., Nature 361 (1993) 467-470]. When expressed in Xenopus oocytes, conflicting results were reported concerning the osmo-sensitivity of the channel [Ji et al., Am. J. Physiol. 275 (1998) C1182-C1190; Hawayda and Subramanyam, J. Gen. Physiol. 112 (1998) 97-111; Rossier, J. Gen. Physiol. 112 (1998) 95-96]. Rat hepatocytes were the first system in which amiloride-sensitive sodium currents in response to hypertonic stress were reported [Wehner et al., J. Gen. Physiol. 105 (1995) 507-535; Wehner et al., Physiologist 40 (1997) A-4]. Moreover, all three ENaC subunits are expressed in these cells [B?hmer et al., Cell. Physiol. Biochem. 10 (2000) 187-194]. Here, we injected specific antisense oligonucleotides directed against alpha-rENaC into single rat hepatocytes in confluent primary culture and found an inhibition of hypertonicity-induced Na(+) currents by 70%. This is the first direct evidence for a role of the ENaC in cell volume regulation.     

Structural mechanisms of acute VEGF effect on microvessel permeability

To investigate the ultrastructural mechanisms of acute microvessel hyperpermeability by vascular endothelial growth factor (VEGF), we combined a mathematical model (J Biomech Eng 116: 502-513, 1994) with experimental data of the effect of VEGF on microvessel hydraulic conductivity (L(p)) and permeability of various-sized solutes. We examined the effect of VEGF on microvessel permeability to a small solute (sodium fluorescein, Stokes radius 0.45 nm), an intermediate solute (alpha-lactalbumin, Stokes radius 2.01 nm), and a large solute [albumin (BSA), Stokes radius 3.5 nm]. Exposure to 1 nM VEGF transiently increased apparent permeability to 2.3, 3.3, and 6.2 times their baseline values for sodium fluorescein, alpha-lactalbumin, and BSA, respectively, within 30 s, and all returned to control within 2 min. On the basis of L(p) (DO Bates and FE Curry. Am J Physiol Heart Circ Physiol 271: H2520-H2528, 1996) and permeability data, the prediction from the model suggested that the most likely structural changes in the interendothelial cleft induced by VEGF would be a approximately 2.5-fold increase in its opening width and partial degradation of the surface glycocalyx.     

Studies on some factors affecting solasodine contents in tissue cultures of Solanum nigrum

Tissue cultures of Solanum nigrum L. were initiated from leaf explants on a solid medium containing inorganic salts [Murashige and Skoog (1962), Physiol. Plant. 15: 473–497], vitamins [Gamborg et al. (1968) Exp. Cell Res. 50:151–158], 3% sucrose and combinations of indoleacetic acid and benzyladenine. Solasodine content was determined in differentiated and undifferentiated (callus) tissues by a colorimetric technique and thin layer chromatography. Indoleacetic acid and sucrose in the medium markedly stimulated the production of solasodine in the tissue cultures. In the cultures grown in darkness the differentiated tissues produced significantly more (anywhere from 1.5 to 10 times) solasodine than the callus in several media. When sucrose concentration was increased to 4, 6 and 10% level in the medium which contained 10 μ M benzyladenine as the sole growth regulator, a significant increase of solasodine production in cultures was found.     

A model of the single atrial cell: relation between calcium current and calcium release

The hypothesis that calcium release from the sarcoplasmic reticulum in cardiac muscle is induced by rises in free cytosolic calcium (Fabiato 1983, Am. J. Physiol 245) allows the possibility that the release could be at least partly regenerative. There would then be a non-linear relation between calcium current and calcium release. We have investigated this possibility in a single-cell version of the rabbit-atrial model developed by Hilgemann & Noble (1987, Proc. R. Soc. Lond. B 230). The model predicts different voltage ranges of activation for calcium-dependent processes (like the sodium-calcium exchange current, contraction or Fura-2 signals) and the calcium current, in agreement with the experimental results obtained by Earm et al. (1990, Proc. R. Soc. Lond. B 240) on exchange current tails, Cannell et al. (1987, Science, Wash. 238) by using Fura-2 signals, and Fedida et al. (1987, J. Physiol., Lond. 385) and Talo et al. (1988, Biology of isolated adult cardiac myocytes) by using contraction. However, when the Fura-2 concentration is sufficiently high (greater than 200 microM) the activation ranges become very similar as the buffering properties of Fura-2 are sufficient to remove the regenerative effect. It is therefore important to allow for the buffering properties of calcium indicators when investigating the correlation between calcium current and calcium release.     

Effect of variable parenchymal expansion on gas mixing

Using implanted radiopaque markers, Hubmayr et al. (J. Appl. Physiol. 54: 1048-1056, 1983) and Olson et al. (J. Appl. Physiol. 57: 1710-1714, 1984) detected a variability in the volume changes of regions defined by the markers in intact and excised dog lungs, respectively. In dogs lying prone and in excised lobes, there is virtually no large-scale spatial organization of the variability. We interpret these data as evidence of an intrinsic heterogeneity of parenchymal expansion. The effect of variability of parenchymal expansion on gas mixing is calculated. From a statistical model, we infer that the variability of volume changes observed by Olson et al. is a result of an underlying variability with a larger magnitude at a smaller scale and that the variability at the smaller scale is large enough to explain the inefficiency of mixing observed in single-breath oxygen tests on excised dog lobes.     

Solubility of Freon 22 in human blood and lung tissue

The solubility of Freon 22 in human blood and lung tissue was determined using the chromatographic method of Wagner et al. (J. Appl. Physiol. 36: 600-605, 1974). In normal human blood, the mean Bunsen coefficient of solubility (alpha B) was 0.804 cm3 STPD.cm-3.ATA-1 at 37 degrees C. It increased with hematocrit (Hct) according to the equation alpha B = 0.274 Hct + 0.691. Tissue homogenates were prepared from macroscopically normal lung pieces obtained at thoracotomy from eight patients undergoing resection for lung carcinoma. The Bunsen solubility coefficients were 0.537 +/- 0.068 and 0.635 +/- 0.091 in washed and unwashed lung, respectively. These values can be used in the determination of both cardiac output and pulmonary tissue volume in humans by use of the rebreathing technique.     

Modeling of respiratory system impedances in dogs

Mechanical impedances between 4 and 64 Hz of the respiratory system in dogs have been reported (A.C. Jackson et al. J. Appl. Physiol. 57: 34-39, 1984) previously by this laboratory. It was observed that resistance (the real part of impedance) decreased slightly with frequency between 4 and 22 Hz then increased considerably with frequency above 22 Hz. In the current study, these impedance data were analyzed using nonlinear regression analysis incorporating several different lumped linear element models. The five-element model of Eyles and Pimmel (IEEE Trans. Biomed. Eng. 28: 313-317, 1981) could only fit data where resistance decreased with frequency. However, when the model was applied to these data the returned parameter estimates were not physiologically realistic. Over the entire frequency range, a significantly improved fit was obtained with the six-element model of DuBois et al. (J. Appl. Physiol. 8: 587-594, 1956), since it could follow the predominate frequency-dependent characteristic that was the increase in resistance. The resulting parameter estimates suggested that the shunt compliance represents alveolar gas compressibility, the central branch represents airways, and the peripheral branch represents lung and chest wall tissues. This six-element model could not fit, with the same set of parameter values, both the frequency-dependent decrease in Rrs and the frequency-dependent increase in resistance. A nine-element model recently proposed by Peslin et al. (J. Appl. Physiol. 39: 523-534, 1975) was capable of fitting both the frequency-dependent decrease and the frequency-dependent increase in resistance. However, the data only between 4 and 64 Hz was not sufficient to consistently determine unique values for all nine parameters.