Search

Recombinant retroviruses are effective vectors able to integrate transgenes into the target cell's genome to achieve longer-term expression. This study investigates the effect of cell lysis products, a common cell culture by-product, on the transduction of suspension cells by gammaretroviral vectors. Cell lysates derived from human and murine suspension cell lines significantly increased the transduction of human TF-1 and K-562 cell lines by gibbon ape leukemia virus-pseudotyped retroviral vectors without altering tropism. The transduction efficiency of TF-1 cells increased as a function of lysate concentration and decreased with increasing target cell concentrations. This was adequately predicted using a saturation equation based on the lysed-to-target cell concentration ratio, R, where: Fold increase = 1+Fold_(Max) (R/(K_(L)+R)). Lysate completely masked the effects of fibronectin when the two were added in combination. With protamine sulfate, the transduction efficiency was increased by lysate to 58% from 20% for protamine sulfate alone. Overall, the presence of cell lysate significantly influenced the outcome of the transduction process, either alone or in the presence of protamine sulfate or fibronectin., ((c) 2009 Wiley Periodicals, Inc.)

A quantitative understanding of the process of retrovirus-mediated gene transfer into mammalian cells should assist the design and optimization of transduction protocols. We present a mathematical model of the process that incorporates the essential rate-limiting transduction steps including diffusion, convection and decay of viral vectors, their binding at the cell surface and entry into the cell cytoplasm, reverse transcription of uncoated RNA to form DNA intermediates, transport of the latter through the cytosol to the cell nucleus and, finally, nuclear import and integration of the delivered DNA into the target cell genome. Cell and virus population balances are used to account for the kinetics of multiple vector infections which influence the transduction efficiency and govern the integrated copy number. The mathematical model is validated using gibbon ape leukemia virus envelope pseudotyped retroviral vectors and K562 target cells. Viral intermediate complexes derived from the internalized retroviral vectors are found to remain stable inside the K562 cells and the cytoplasmic trafficking time is consistent with the time scale for retrovirus uncoating, reverse transcription and transport to the cell nucleus. The model predictions of transduction efficiency and integrated copy number agree well with experimental data for both static (i.e., standard gravity) and centrifugation-based gene transfer protocols. The formulation of the model can also be applied to transduction protocols involving lenti- or foamy-viruses and so should prove to be useful for the optimization of several types of gene transfer processes.

The expansion of stem cell numbers while retaining their developmental properties is a bioprocess challenge. We compared the growth rates and embryoid body (EB) formation yields of R1 and EFC murine embryonic stem cells (mESC) cultured in two basal media (DMEM or DMEM:F12) with additions of 1.7-15% fetal bovine serum (FBS) or serum replacer (KOSR). Whereas the basal medium or KOSR dose did not have a significant effect on growth rate for either cell line, increasing doses of KOSR had a significant negative effect on the EB yield of EFC cells. Use of DMEM:F12 and increasing doses of FBS independently and significantly increased the growth rate for both cell lines. DMEM:F12 also significantly increased EB yields for both cell lines. The results show that use of DMEM:F12 and several-fold lower than conventional concentrations of KOSR can efficiently support maintenance of mESC and that KOSR should be dose as well as lot optimized.

A cell retention device that provides reliable high-separation efficiency with minimal negative effects on the cell culture is essential for robust perfusion culture processes. External separation devices generally expose cells to periodic variations in temperature, most commonly temperatures below 37 degrees C, while the cells are outside the bioreactor. To examine this phenomenon, aliquots of approximately 5% of a CHO cell culture were exposed to 60 s cyclic variations of temperature simulating an acoustic separator environment. It was found that, for average exposure temperatures between 31.5 and 38.5 degrees C, there were no significant impacts on the rates of growth, glucose consumption, or t-PA production, defining an acceptable range of operating temperatures. These results were subsequently confirmed in perfusion culture experiments for average exposure temperatures between 31.6 and 38.1 degrees C. A 2(5-1) central composite factorial design experiment was then performed to systematically evaluate the effects of different operating variables on the inlet and outlet temperatures of a 10L acoustic separator. The power input, ambient temperature, as well as the perfusion and recycle flow rates significantly influenced the temperature, while the cell concentration did not. An empirical model was developed that predicted the temperature changes between the inlet and the outlet of the acoustic separator within +/-0.5 degrees C. A series of perfusion experiments determined the ranges of the significant operational settings that maintained the acoustic separator inlet and outlet temperatures within the acceptable range. For example, these objectives were always met by using the manufacturer-recommended operational settings as long as the recirculation flow rate was maintained above 15 L day(-1) and the ambient temperature was near 22 degrees C.

Controlled shear affinity filtration (CSAF) is a novel integrated processing technology that positions a rotor directly above an affinity membrane chromatography column to permit protein capture and purification directly from cell culture. The conical rotor is intended to provide a uniform and tunable shear stress at the membrane surface that inhibits membrane fouling and cell cake formation by providing a hydrodynamic force away from and a drag force parallel to the membrane surface. Computational fluid dynamics (CFD) simulations are used to show that the rotor in the original CSAF device (Vogel et al., 2002) does not provide uniform shear stress at the membrane surface. This results in the need to operate the system at unnecessarily high rotor speeds to reach a required shear stress of at least 0.17 Pa at every radial position of the membrane surface, compromising the scale-up of the technology. Results from CFD simulations are compared with particle image velocimetry (PIV) experiments and a numerical solution for low Reynolds number conditions to confirm that our CFD model accurately describes the hydrodynamics in the rotor chamber of the CSAF device over a range of rotor velocities, filtrate fluxes, and (both laminar and turbulent) retentate flows. CFD simulations were then carried out in combination with a root-finding method to optimize the shape of the CSAF rotor. The optimized rotor geometry produces a nearly constant shear stress of 0.17 Pa at a rotational velocity of 250 rpm, 60% lower than the original CSAF design. This permits the optimized CSAF device to be scaled up to a maximum rotor diameter 2.5 times larger than is permissible in the original device, thereby providing more than a sixfold increase in volumetric throughput., (Copyright 2006 Wiley Periodicals, Inc.)

Acoustic cell filters operate at high separation efficiencies with minimal fouling and have provided a practical alternative for up to 200 L/d perfusion cultures. However, the operation of cell retention systems depends on several settings that should be adjusted depending on the cell concentration and perfusion rate. The impact of operating variables on the separation efficiency performance of a 10-L acoustic separator was characterized using a factorial design of experiments. For the recirculation mode of separator operation, bioreactor cell concentration, perfusion rate, power input, stop time and recirculation ratio were studied using a fractional factorial 2(5-1) design, augmented with axial and center point runs. One complete replicate of the experiment was carried out, consisting of 32 more runs, at 8 runs per day. Separation efficiency was the primary response and it was fitted by a second-order model using restricted maximum likelihood estimation. By backward elimination, the model equation for both experiments was reduced to 14 significant terms. The response surface model for the separation efficiency was tested using additional independent data to check the accuracy of its predictions, to explore robust operation ranges and to optimize separator performance. A recirculation ratio of 1.5 and a stop time of 2 s improved the separator performance over a wide range of separator operation. At power input of 5 W the broad range of robust high SE performance (95% or higher) was raised to over 8 L/d. The reproducible model testing results over a total period of 3 months illustrate both the stable separator performance and the applicability of the model developed to long-term perfusion cultures., (Copyright (c) 2005 Wiley Periodicals, Inc.)

There is an expanding need for predictive mathematical models to accelerate the optimization of cell therapy culture processes. Here we demonstrate the ability of simple mathematical models to describe quantitatively the cytokine growth-rate dependence of two human hematopoietic cell lines, TF-1 and MO7e. These cells are immortal but depend on either interleukin-3 (IL-3) or granulocyte-macrophage colony stimulating factor (GM-CSF) for their continued survival and maximal proliferation. They are also responsive to interleukin-6 (IL-6) and exhibit saturation kinetics when these cytokines are limiting. A Monod-type relationship consistently failed to fit measured cytokine dose-proliferation response curves while a Hill-type relationship showed a good fit. Cytokine interactions were first modeled by modifying the Hill-function to include an interaction parameter, gamma. This model did not indicate either synergistic or even additive effects between IL-3 and GM-CSF. Based on the reported competition between IL-3 and GM-CSF for their common receptor (beta(c)) subunit, a competitive model was also developed. This model had no new parameters beyond those obtained from single cytokine cultures and provided improved prediction of the growth rates for both cell lines exposed to combinations of IL-3 and GM-CSF over a wide range of concentrations. As expected, the competitive model failed to fit the data for IL-6 in combination with either IL-3 or GM-CSF, since IL-6 signaling does not involve the beta(c) chain of the IL-3/GM-CSF receptors. Interestingly, the cell-specific rates of GM-CSF uptake and cell proliferation were found to be uncoupled processes. Taken together, these results illustrate the utility of appropriately designed empirical models to describe the proliferative responses of hematopoietic cells to cytokine stimulation., (Copyright 2004 Wiley Periodicals, Inc.)

The distribution and transport of fluid, ions, and other solutes (plasma proteins and glucose) are described in a mathematical model of unresuscitated hemorrhage. The model is based on balances of each material in both the circulation and its red blood cells, as well as in a whole-body tissue compartment along with its cells. Exchange between these four compartments occurs by a number of different mechanisms. The hemorrhage model has as its basis a validated model, due to Gyenge et al., of fluid and solute exchange in the whole body of a standard human. Hypothetical but physiologically based features such as glucose and small ion releases along with cell membrane changes are incorporated into the hemorrhage model to describe the system behavior, particularly during larger hemorrhages. Moderate (10%-30% blood volume loss) and large (> 30% blood loss) hemorrhage dynamics are simulated and compared with available data. The model predictions compare well with the available information for both types of hemorrhages and provide a reasonable description of the progression of a large hemorrhage from the compensatory phase through vascular collapse.

This study is concerned with the formulation of a 'kidney module' linked to the plasma compartment of a larger mathematical model previously developed. Combined, these models can be used to predict, amongst other things, fluid and small ion excretion rates by the kidney; information that should prove useful in evaluating values and trends related to whole-body fluid balance for different clinical conditions to establish fluid administration protocols and for educational purposes. The renal module assumes first-order, negative-feedback responses of the kidney to changes in plasma volume and/or plasma sodium content from their normal physiological set points. Direct hormonal influences are not explicitly formulated in this empiric model. The model also considers that the renal excretion rates of small ions other than sodium are proportional to the excretion rate of sodium. As part of the model development two aspects are emphasized (1): the estimation of parameters related to the renal elimination of fluid and small ions, and (2) model validation via comparisons between the model predictions and selected experimental data. For validation, model predictions of the renal dynamics are compared with new experimental data for two cases: plasma overload resulting from external fluid infusion (e.g. infusions of iso-osmolar solutions and/or hypertonic/hyperoncotic saline solutions), and untreated hypo volemic conditions that result from the external loss of blood. The present study demonstrates that the empiric kidney module presented above can provide good short-term predictions with respect to all renal outputs considered here. Physiological implications of the model are also presented., (Copyright Acta Anaesthesiologica Scandinavica 47 (2003))

Increasing worldwide demand for mammalian cell production capacity will likely be partially satisfied by a greater use of higher volumetric productivity perfusion processes. An important additional component of any perfusion system is the cell retention device that can be based on filtration, sedimentation, and/or acoustic technologies. A common concern with these systems is that pumping and transient exposure to suboptimal medium conditions may damage the cells or influence the product quality. A novel air-backflush mode of operating an acoustic cell separator was developed in which an injection of bioreactor air downstream of the separator periodically returned the captured cells to the reactor, allowing separation to resume within 20 s. This mode of operation eliminated the need to pump the cells and allows the selection of a residence time in the separator depending on the sensitivity of the cell line. The air-backflush mode of operating a 10L acoustic separator was systematically tested at 10(7) cells/mL to define reliable ranges of operation. Consistent separation performance was obtained for wide ranges of cooling airflow rates from 0 to 15 L/min and for backflush frequencies between 10 and 40 h(-1). The separator performance was optimized at a perfusion rate of 10 L/day to obtain a maximum separation efficiency of 92 +/- 0.3%. This was achieved by increasing the power setting to 8 W and using duty cycle stop and run times of 4.5 and 45 s, respectively. Acoustic cell separation with air backflush was successfully applied over a 110 day CHO cell perfusion culture at 10(7) cells/mL and 95% viability.

In order to understand better the complex, dynamic behaviour of the redistribution and exchange of fluid and solutes administered to normal individuals or to those with acute hypovolemia, mathematical models are used in addition to direct experimental investigation. Initial validation of a model developed by our group involved data from animal experiments (Gyenge, C.C., Bowen, B.D., Reed, R.K. & Bert, J.L. 1999b. Am J Physiol 277 (Heart Circ Physiol 46), H1228-H1240). For a first validation involving humans, we compare the results of simulations with a wide range of different types of data from two experimental studies. These studies involved administration of normal saline or hypertonic saline with Dextran to both normal and 10% haemorrhaged subjects. We compared simulations with data including the dynamic changes in plasma and interstitial fluid volumes VPL and VIT respectively, plasma and interstitial colloid osmotic pressures PiPL and PiIT respectively, haematocrit (Hct), plasma solute concentrations and transcapillary flow rates. The model predictions were overall in very good agreement with the wide range of experimental results considered. Based on the conditions investigated, the model was also validated for humans. We used the model both to investigate mechanisms associated with the redistribution and transport of fluid and solutes administered following a mild haemorrhage and to speculate on the relationship between the timing and amount of fluid infusions and subsequent blood volume expansion.

A compartmental model of short-term whole body fluid, protein, and ion distribution and transport is formulated. The model comprises four compartments: a vascular and an interstitial compartment, each with an embedded cellular compartment. The present paper discusses the assumptions on which the model is based and describes the equations that make up the model. Fluid and protein transport parameters from a previously validated model as well as ionic exchange parameters from the literature or from statistical estimation [see companion paper: C. C. Gyenge, B. D. Bowen, R. K. Reed, and J. L. Bert. Am. J. Physiol. 277 (Heart Circ. Physiol. 46): H1228-H1240, 1999] are used in formulating the model. The dynamic model has the ability to simulate 1) transport across the capillary membrane of fluid, proteins, and small ions and their distribution between the vascular and interstitial compartments; 2) the changes in extracellular osmolarity; 3) the distribution and transport of water and ions associated with each of the cellular compartments; 4) the cellular transmembrane potential; and 5) the changes of volume in the four fluid compartments. The validation and testing of the proposed model against available experimental data are presented in the companion paper.

A mathematical model of short-term whole body fluid, protein, and ion distribution and transport developed earlier [see companion paper: C. C. Gyenge, B. D. Bowen, R. K. Reed, and J. L. Bert. Am. J. Physiol. 277 (Heart Circ. Physiol. 46): H1215-H1227, 1999] is validated using experimental data available in the literature. The model was tested against data measured for the following three types of experimental infusions: 1) hyperosmolar saline solutions with an osmolarity in the range of 2,000-2,400 mosmol/l, 2) saline solutions with an osmolarity of approximately 270 mosmol/l and composition comparable with Ringer solution, and 3) an isosmotic NaCl solution with an osmolarity of approximately 300 mosmol/l. Good agreement between the model predictions and the experimental data was obtained with respect to the trends and magnitudes of fluid shifts between the intra- and extracellular compartments, extracellular ion and protein contents, and hematocrit values. The model is also able to yield information about inaccessible or difficult-to-measure system variables such as intracellular ion contents, cellular volumes, and fluid fluxes across the vascular capillary membrane, data that can be used to help interpret the behavior of the system.

Within the spectrum of current applications for cell culture technologies, efficient large-scale mammalian cell production processes are typically carried out in stirred fed-batch or perfusion bioreactors. The specific aspects of each individual process that can be considered when determining the method of choice are presented. A major challenge for perfusion reactor design and operation is the reliability of the cell retention device. Current retention systems include cross-flow membrane filters, spin-filters, inclined settlers, continuous centrifuges and ultrasonic separators. The relative merits and limitations of these technologies for cell retention and their suitability for large-scale perfusion are discussed.

A dynamic compartmental model is developed to describe the redistribution of fluid and albumin between the circulation and the intact and injured interstitia following burn injury in humans. Transcapillary fluid and albumin exchange is described by a coupled Starling mechanism, while the effect of the burn is represented by time-dependent perturbations to all three compartments. The unknown model parameters are determined for two groups of patients, having less than and greater than 25% total body surface area burns, by statistical fitting of model predictions to patient data from two sources. The parameters include the perturbations to the fluid filtration coefficients in uninjured and injured tissue, GkF,Tl and GkF,BT, respectively, the relaxation coefficient, r, which describes the exponential decay of the perturbations, and the exudation factor, EXFAC, which relates the protein concentration in the exudate to that in the injured tissue. Perturbations to other parameters, including the membrane permeability-surface area product and the albumin reflection coefficient in the injured and uninjured tissues, are determined based on interrelationships with GkF,Tl and GkF,BT. The values of GkF,BT, when corrected for tissue destruction and decreased post-injury perfusion, are in reasonable agreement with the limited experimental data available from the literature. The model and its parameters are further validated by comparing the simulated patient responses to the clinical data used in the parameter estimation as well as to data available from two additional sources.

A compartmental model consisting of the circulation, a general interstitium, and the lymphatics, is formulated to describe the transport and distribution of fluid and plasma proteins (albumin) in the human microvascular exchange system. Transcapillary mass exchange is assumed to occur via a coupled Starling mechanism. Unknown or poorly quantified model parameters are estimated by statistical fitting of simulation predictions to five different sets of experimental data. The data consist of steady-state and transient plasma and interstitial volumes and colloid osmotic pressures measured under laboratory or clinical conditions for normal humans and for patients with nephrotic syndrome or mild heart disease. In all cases, it is assumed that the system response to perturbations imposed either artificially or through illness is due to changes in the Starling driving forces. The three best-fit parameters were found to be normal capillary hydrostatic pressure, Pc,o = 11.0 mm Hg; albumin reflection coefficient, sigma = 0.99; and lymph flow sensitivity, LS = 43.1 ml/mm Hg.hr. Three other parameters, which were unknown but related to the estimated parameters through steady-state mass balance equations, were determined to be fluid filtration coefficient, KF = 121.1 ml/mm Hg.hr; albumin permeability-surface area product, PS = 73.0 ml/hr; and normal lymph flow, JL,o = 75.7 ml/hr. The fully described model was validated by comparisons between (1) simulation predictions and data used in parameter estimation, (2) estimated transport parameters and available literature values, and (3) model predictions and an additional set of experimental data.

The recovery of serum-free medium proteins from poly-sulfone hollow fiber bioreactors (HFBRs) was investigated. More than 99% of the initial transferrin was adsorbed to the hydrophobic hollow fibers within 2 h of HFBR operation. A methodology to minimize transferrin adsorption by pre-adsorption of bovine serum albumin (BSA) was developed. BSA adsorption on suspended cut fibers was virtually complete within 1 h. BSA-coated fibers adsorbed only 5% of the transferrin within 10 days, whereas uncoated cut fibers adsorbed more than 99% of the transferrin within 1 h. An improved HFBR startup procedure, using a BSA-coating step before inoculation, resulted in substantially higher transferrin recovery. Additional factors influenced extracapillary space (ECS) transferrin concentrations. Pronounced downstream polarization of transferrin was observed in the ECS. In addition, the 30-kDa nominal molecular weight cutoff ultrafiltration membranes rapidly leaked transferrin from the ECS to the lumen.

A mathematical model is formulated and used to describe the distribution and transport of fluid and albumin in the human circulation, interstitium and lymphatics. Two transcapillary mass exchange mechanisms are investigated: a homoporous 'coupled Starling model', in which transcapillary albumin diffusion and convection occur within the same pathway, and a heteroporous 'plasma leak model', in which variations in structure and pressure are permitted along the length of the capillary. Parameters used in the transport models are determined based on statistical fitting of simulation predictions to experimental data from normal humans and nephrotic patients. The data consists of interstitial fluid volumes and interstitial colloid osmotic pressures as functions of plasma colloid osmotic pressure. Model validation is carried out based on comparison of (i) simulation predictions with experimental data used in parameter estimation, (ii) estimated transport parameters with experimentally determined values, and (iii) simulation predictions with a set of dynamic data from an albumin infusion study. While both models with their best-fit parameter estimates provide a good representation of experimental data, the drawbacks of the plasma leak model are three-fold: it requires more estimated parameters than the coupled Starling model, little experimental information exists with which to compare these parameters and, with the best fit values obtained, the plasma leak mechanism becomes insignificant. The model that employs a Starling-type exchange mechanism will therefore be favoured in future applications.

In this paper, a mathematical model of interstitial transport and microvascular exchange within a rigid mesenteric tissue segment is employed to simulate the transient exchange of fluid and plasma proteins following two systemic disturbances: hypoproteinemia and venous congestion. In each case, the model system behavior is studied as a function of interstitial plasma protein transport mechanisms and mesothelial transport properties. Plasma protein washout was generally predicted in cases of hypoproteinemia. However, following venous congestion, the transient change in interstitial plasma protein content also depended on the relative sieving properties of the filtering and draining boundaries. When these boundaries display similar sieving characteristics, the interstitial plasma protein content increases following the disturbance. Such behavior may have some bearing on transient exchange in the hepatic microcirculation during venous congestion.

The present work is a continuation of studies concerned with mathematical modelling and simulation of microvascular fluid and protein exchange following burn injuries [Bert et al.: Circulatory Shock 28: 199-219, 1989: Bowen et al.: Circulatory Shock 28: 221-233, 1989]. The model has been extended to include the effects of different types of fluid resuscitation on the circulatory and microvascular exchange systems. The model and a statistical fitting procedure were used to find the ranges of fitting parameter values that best describe the changes in interstitial fluid volume and protein mass as well as transcapillary protein extravasation for three sets of experiments (no resuscitation, resuscitation with Ringer's or resuscitation with plasma). Typical changes in mass exchange related parameters postburn that resulted in simulation predictions which were a good fit to the experimental data include: an increase in the large pore pathway for protein of 100 times in the injured skin and 5 times in non-injured skin and skeletal muscle, an increase in fluid filtration coefficients in injured skin of 10 times and an instantaneous decrease of 50% in the area available for exchange in injured skin at the time of the burn.

A simplified version of the model of interstitial transport developed earlier (D. G. Taylor, J. L. Bert, and B. D. Bowen, 1990, Microvasc. Res. 39, 253-278) is used to investigate microvascular exchange of fluid and a single "aggregate" plasma protein species in mesenteric tissue. The interstitium is approximated by a rigid, rectangular, porous slab displaying two fluid pathways, only one of which is available to plasma proteins. The model is used to explore the effects of the interstitial plasma protein diffusivity, the tissue hydraulic conductivity, the restricted convection of plasma proteins, and the mesothelial transport characteristics on the steady-state distribution and transport of plasma proteins and flow of fluid in the tissue. The simulations predict significant convective plasma protein transport and complex fluid flow patterns within the interstitium. These flow patterns can produce local regions of high fluid and plasma protein exchange along the mesothelium which might be erroneously identified as "leaky sites."

A generalized model is developed to describe the transport of fluid and plasma proteins or other macromolecules within the interstitium. To account for the effects of plasma protein exclusion and interstitial swelling, the interstitium is treated as a multiphase deformable porous medium. Fluid flow is assumed proportional to the gradient in fluid chemical potential and therefore depends not only on the local hydrostatic pressure but also on the local plasma protein concentrations through appropriate colloid osmotic pressure relationships. Plasma protein transport is assumed to occur by restricted convection, molecular diffusion, and convective dispersion. In a companion paper (D. G. Taylor, J. L. Bert, and B. D. Bowen, 1990, Microvasc. Res. 39, 279-306) a simplified version of the model is used to analyze steady-state fluid and plasma protein exchange within mesentery.

A dynamic mathematical model is formulated and used to describe the distribution and transport of fluid and plasma proteins between the circulation, interstitial space of skin and muscle, and the lymphatics in the rat. Two descriptions of transcapillary exchange are investigated: a homoporous "Starling model" and a heteroporous "plasma leak model." Parameters used in the two hypothetical transport mechanisms are determined based on statistical fitting procedures between simulation predictions and selected experimental data. These data consist of interstitial fluid volume and colloid osmotic pressure measurements as a function of venous pressure for muscle and interstitial colloid osmotic pressure vs. venous pressure for skin. The values determined for the transport parameters compare well with data in the literature. The fully determined model is used to simulate steady-state conditions of hypoproteinemia, overhydration, and dehydration, as well as the dynamic response to changes in venous pressure and intravascularly administered protein tracers. Comparisons between the simulation predictions and experimental data for these various perturbations are made. The plasma leak model appears to provide a better description of microvascular exchange.

A mathematical model of microvascular exchange in the rat following a burn injury was developed by extending an existing model of normal microvascular exchange to include perturbations characteristic of burn injuries without fluid resuscitation. The changes anticipated for small (10% body surface area) and large (40% body surface area) burns are incorporated systematically into the model until there is no improvement in the statistical fit of the simulation predictions with the experimental data of Lund and Reed (Circulatory Shock 20:91-104, 1986). The "best fit" perturbations for the small burn include the experimentally measured changes in mean arterial pressure and injured tissue pressure as well as changes to plasma protein and fluid transport coefficients in the injured tissue. The larger burn "best fit" simulation required changes to the plasma protein transport coefficients in the intact tissues as well as all of the changes listed above. The simulation results are compared with the available experimental information on burn injuries as well as with the specific data of Lund and Reed (Circulatory Shock 20:91-104, 1986).

The present work uses and extends a dynamic mathematical model [J. L. Bert, B. D. Bowen, and R. K. Reed. Am. J. Physiol. 254 (Heart Circ. Physiol. 23): H384-H399, 1988] to investigate microvascular exchange and interstitial fluid volume regulation in the rat. Alternative concepts of transcapillary exchange as well as other parametric changes were incorporated into the model. In all cases, predictions resulting from these changes did not describe the available experimental information as well as the original model. A sensitivity analysis of the model showed the microvascular exchange system to be well regulated near its normal steady-state conditions through passive readjustment of the forces participating in the volume regulation. The transient rates of fluid and protein exchange were studied in order to determine the mechanisms inherent in the model that lead to fluid volume regulation during episodes of increased venous pressure and hypoproteinemia. In addition to interstitial compliance, lymph flow characteristics, and washdown of interstitial proteins, it was found that the magnitude and direction of reabsorption played an important role in the regulation process. Edema was always associated with a permanent reversal of the reabsorptive flow.

The present work investigates the implications of the predictions of a dynamic mathematical model of microvascular exchange following a nonresuscitated burn injury in a rat (Bert et al.: Circulatory Shock 28:199-219, 1989). Transport coefficients, transmicrovascular pressures, and the resultant fluid and protein fluxes were examined in order to assess their quantitative importance to the dynamic behavior of small (10% body surface area) and large (40% body surface area) burns. Edema accumulation in the injured tissue is dependent not only on events occurring in that tissue but is influenced strongly by interaction with the plasma and the noninjured tissue compartments.