Multiscale modeling of lymphatic drainage from tissues using homogenization theory

Lymphatic capillary drainage of interstitial fluid under both steady-state and inflammatory conditions is important for tissue fluid balance, cancer metastasis, and immunity. Lymphatic drainage function is critically coupled to the fluid mechanical properties of the interstitium, yet this coupling is poorly understood. Here we sought to effectively model the lymphatic-interstitial fluid coupling and ask why the lymphatic capillary network often appears with roughly a hexagonal architecture. We use homogenization method, which allows tissue-scale lymph flow to be integrated with the microstructural details of the lymphatic capillaries, thus gaining insight into the functionality of lymphatic anatomy. We first describe flow in lymphatic capillaries using the Navier-Stokes equations and flow through the interstitium using Darcy's law. We then use multiscale homogenization to derive macroscale equations describing lymphatic drainage, with the mouse tail skin as a basis. We find that the limiting resistance for fluid drainage is that from the interstitium into the capillaries rather than within the capillaries. We also find that between hexagonal, square, and parallel tube configurations of lymphatic capillary networks, the hexagonal structure is the most efficient architecture for coupled interstitial and capillary fluid transport; that is, it clears the most interstitial fluid for a given network density and baseline interstitial fluid pressure. Thus, using homogenization theory, one can assess how vessel microstructure influences the macroscale fluid drainage by the lymphatics and demonstrate why the hexagonal network of dermal lymphatic capillaries is optimal for interstitial tissue fluid clearance.     

Functional recovery of fluid drainage precedes lymphangiogenesis in acute murine foreleg lymphedema

Secondary lymphedema in humans is a common consequence of axillary lymph node dissection (ALND) to treat breast cancer. It is commonly hypothesized that lymphatic growth is required to increase fluid drainage and ameliorate lymphedema. Although there is a pronounced alteration in the balance of interstitial forces regulating fluid transport that sustains the chronic form of lymphedema, it is presently unknown whether changes occur to the balance of interstitial forces during acute lymphedema that may play a role in the recovery of fluid drainage. Here, we compared the relative importance of lymphangiogenesis of lymphatic vessels and interstitial flows for restoring fluid drainage and resolving acute lymphedema in the mouse foreleg after ALND. We found that removal of the axillary lymph nodes reduced lymph drainage in the foreleg at days 0 and 5 postsurgery, with fluid tracer spreading interstitially through subcutaneous tissues. Interstitial fluid drainage returned to normal by day 10, whereas functional regrowth of lymphatic vessels was first detected by indocyanine green fluorescence lymphography at day 15, demonstrating that the recovery of interstitial fluid drainage preceded the regrowth of lymphatic vessels. This was confirmed by the administration of VEGF receptor-3-neutralizing antibodies, which completely blocks lymphatic regrowth. It was found that the recovery of interstitial fluid drainage and the natural resolution of acute lymphedema produced by ALND were not hindered by VEGF receptor-3 neutralization, demonstrating that interstitial fluid drainage recovery and the resolution of acute lymphedema are lymphangiogenesis independent. The data highlight the central role of the interstitial environment in adapting to lymphatic injury to increase fluid drainage.     

Influence of interstitial fluid dynamics on growth and therapy of angiogenic tumor. Analysis by mathematical model

We have developed a spatially distributed mathematical model of angiogenic tumor growth in tissue with account of interstitial fluid dynamics and bevacizumab monotherapy. In this model the process of neovascularization is initiated by tumor cells in a state of metabolic stress, vascular endothelial growth factor (VEGF) being its main mediator. The model takes into consideration the convection flows arising in dense tissue due to active proliferation and migration of tumor cells as well as interstitial fluid inflow from blood vascular system, its outflow through lymphatic system and redistribution in the area of tumor growth. The work considers the diffusive approximation of interstitial fluid dynamics in tumor and normal tissue. Numerical study of the model showed that in absence of therapy a peritumoral edema is formed due to the increase of interstitial fluid inflow from angiogenic capillaries. In the case of rapid interstitial fluid outflow through lymphatic system and its fast transport from necrotic zone to normal tissue the regimes of full growth stop are observed in case of low-invasive tumor. Under bevacizumab monotherapy the peritumoral edema vanishes and low-invasive tumor may not only decelerate its growth, but also start shrinking for a large range of parameters.     

Fibroblast alignment under interstitial fluid flow using a novel 3-D tissue culture model

Interstitial flow is an important component of the microcirculation and interstitial environment, yet its effects on cell organization and tissue architecture are poorly understood, in part due to the lack of in vitro models. To examine the effects of interstitial flow on cell morphology and matrix remodeling, we developed a tissue culture model that physically supports soft tissue cultures and allows microscopic visualization of cells within the three-dimensional matrix. In addition, pressure-flow relationships can be continuously monitored to evaluate the bulk hydraulic resistance as an indicator of changes in the overall matrix integrity. We observed that cells such as human dermal fibroblasts aligned perpendicular to the direction of interstitial flow. In contrast, fibroblasts in static three-dimensional controls remained randomly oriented, whereas cells subjected to fluid shear as a two-dimensional monolayer regressed. Also, the dynamic measurements of hydraulic conductivity suggest reorganization toward a steady state. These primary findings help establish the importance of interstitial flow on the biology of tissue organization and interstitial fluid balance.     

Modeling of interstitial fluid movement in soft tissue under negative pressure – relevance to treatment of tissue swelling

Exact physiological mechanisms behind the potential positive treatment effects of pathological tissue swelling (edema), such as increased interstitial fluid flow, are poorly understood. Finite-element model was created and the model response was matched with the deformation data from the negative pressure (suction) measurements in human (N = 11) forearm. Two experimental suction protocols were simulated to evaluate their impact on interstitial fluid flow in soft tissues. Simulated continuous suction was up to 27 times more efficient in fluid transportation compared to the cyclic suction. The continuous suction that transports the interstitial fluid effectively may help to decrease soft tissue edema.     

Edemagenic gain and interstitial fluid volume regulation

Under physiological conditions, interstitial fluid volume is tightly regulated by balancing microvascular filtration and lymphatic return to the central venous circulation. Even though microvascular filtration and lymphatic return are governed by conservation of mass, their interaction can result in exceedingly complex behavior. Without making simplifying assumptions, investigators must solve the fluid balance equations numerically, which limits the generality of the results. We thus made critical simplifying assumptions to develop a simple solution to the standard fluid balance equations that is expressed as an algebraic formula. Using a classical approach to describe systems with negative feedback, we formulated our solution as a "gain" relating the change in interstitial fluid volume to a change in effective microvascular driving pressure. The resulting "edemagenic gain" is a function of microvascular filtration coefficient (K(f)), effective lymphatic resistance (R(L)), and interstitial compliance (C). This formulation suggests two types of gain: "multivariate" dependent on C, R(L), and K(f), and "compliance-dominated" approximately equal to C. The latter forms a basis of a novel method to estimate C without measuring interstitial fluid pressure. Data from ovine experiments illustrate how edemagenic gain is altered with pulmonary edema induced by venous hypertension, histamine, and endotoxin. Reformulation of the classical equations governing fluid balance in terms of edemagenic gain thus yields new insight into the factors affecting an organ's susceptibility to edema.     

The resolution of lymphedema by interstitial flow in the mouse tail skin

Lymphangiogenesis is considered a promising approach for increasing fluid drainage during secondary lymphedema. However, organization of lymphatics into functional capillaries may be dependent upon interstitial flow (IF). The present study was undertaken to determine the importance of lymphangiogenesis for lymphedema resolution. We created a lymphatic obstruction that produces lymphedema in mouse tail skin. The relatively scar-free skin regeneration that occurred across the obstruction allowed the progression of lymphangiogenesis to be observed and compared with the evolution of lymphedema. The role of vascular endothelial growth factor-C (VEGF-C)/VEGF receptor (VEGFR)-3 signaling in lymphedema resolution was investigated by exogenous administration of VEGF-C or neutralizing antibodies against VEGFR-3. VEGF-C protein improved lymphedema at 15 days [reducing dermal thickness from 742 +/- 105 to 559 +/- 141 microm with 95% confidence intervals (CIs), P < 0.05] without increasing lymphatic capillary coverage (11.6 +/- 6.4% following VEGF-C treatment relative to 9.6 +/- 6.2% with 95% CIs, P > 0.50). Blocking VEGFR-3 signaling did not inhibit lymphedema resolution at 25 days (dermal thickness of 462 +/- 127 microm following VEGFR-3 inhibition relative to 502 +/- 87 microm with 95% CIs) or inhibit IF, although VEGFR-3 blocking prevented lymphangiogenesis (reducing lymphatic coverage to 0.2 +/- 0.7% relative to 8.7 +/- 7.3% with 95% CIs, P < 0.005). A second mouse tail lymphedema model was employed to investigate the ability of VEGF-C to increase fluid drainage across a scar. We found that neither neutralization of VEGFR-3 nor administration of VEGF-C affected the course of skin swelling over 25 days. These findings suggest that resolution of lymphedema in the mouse tail skin may be more dependent upon IF and regeneration of the extracellular matrix across the obstruction than lymphatic capillary regeneration.     

The structure of lymphatic capillaries in lymph formation.

The lymphatic vascular system consists of endothelial lined vessels which begin as blind-end tubes or saccules that are located within the connective tissue areas. This system serves as a one-way drainage apparatus for the removal of diffusible substances as well as plasma proteins that escape the blood capillaries. If permitted to accumulate, these escaped components would deplete the circulatory system of its plasma colloids and disrupt the balance of forces responsible for the control of fluid movement and the exchange of gases and fluids across the blood vascular wall. The lymphatic capillaries are strategically placed and anatomically constructed to permit a continuous and rapid removal of the transient interstitial fluids, plasma proteins, and cells from the interstitium. Structurally the lymphatic capillaries consist of a continuous endothelium that is extremely attenuated over major aspects of its diameter, except in the perinuclear region which bulges into the lumen. These vessels lack a continuous basal lamina and maintain a close relationship with the adjoining interstitium by way of anchoring filaments. The adjacent cells are extensively overlapped and lack adhesion devices in many areas. When electron-opaque tracers are injected intravenously (i.e., horseradish peroxidase and ferritin), subsequent electron microscopic examination of tissues reveals the presence of tracer particles within the interstitium and the lymphatic capillary lumen. These particles gain access into the lymphatic capillaries via two major pathways: 1) the intercellular clefts of patent junctions and 2) plasmalemmal vesicles (pinocytotic vesicles). Another salient feature of the lymphatic endothelial cell includes the presence of numerous cytoplasmic filaments, which are similar in morphology to the actin filaments observed in a variety of cell types. The ultrastructural features of the lymphatic capillaries are discussed in relation to their role in the removal of interstitial fluids and particulate matter, and in the formation of lymph.     

Clinical assessment of human lymph flow using removal rate constants of interstitial macromolecules: a critical review of lymphoscintigraphy

Edema is a common clinical problem, and the daily avoidance of edema depends critically on the lymphatic system, which clears leaked plasma proteins and fluid from the interstitial compartment. There is often confusion as to the difference between chronic edema and lymphedema. Lymphedema is by definition primarily a disease of impaired lymphatic drainage and lymph flow, and progress in lymphedema research, currently an increasingly active field, requires a clinically viable method for the quantitative assessment of lymph drainage rate in patients. Measurement of the rate of clearance of a new protein marker, radiolabelled human immunoglobulin, from skin, subcutis, and muscle provides a way of measuring human lymph flow quantitatively and is the only viable clinical method currently available. Considerable strides have been made over the last 5-10 years in evaluating the method and its pitfalls, including potential complications such as vascular clearance, peripheral lymphovenous communications and label dissociation. The review assesses critically, for the first time, the evidence relating to the method: its pitfalls; human lymph flow in various healthy and oedematous tissues; and how this is altered in hyperfiltration edemas, inflammation, vasoconstriction and various primary and secondary human lymphedemas.     

Mechanoinduction of lymph vessel expansion

In the mammalian embryo, few mechanical signals have been identified to influence organ development and function. Here, we report that an increase in the volume of interstitial or extracellular fluid mechanically induces growth of an organ system, that is, the lymphatic vasculature. We first demonstrate that lymph vessel expansion in the developing mouse embryo correlates with a peak in interstitial fluid pressure and lymphatic endothelial cell (LEC) elongation. In 'loss-of-fluid' experiments, we then show that aspiration of interstitial fluid reduces the length of LECs, decreases tyrosine phosphorylation of vascular endothelial growth factor receptor-3 (VEGFR3), and inhibits LEC proliferation. Conversely, in 'gain-of-fluid' experiments, increasing the amount of interstitial fluid elongates the LECs, and increases both VEGFR3 phosphorylation and LEC proliferation. Finally, we provide genetic evidence that β1 integrins are required for the proliferative response of LECs to both fluid accumulation and cell stretching and, therefore, are necessary for lymphatic vessel expansion and fluid drainage. Thus, we propose a new and physiologically relevant mode of VEGFR3 activation, which is based on mechanotransduction and is essential for normal development and fluid homeostasis in a mammalian embryo.     

Fluid Flow along Venous Adventitia in Rabbits: Is It a Potential Drainage System Complementary to Vascular Circulations?

Background

Our previous research and other studies with radiotracers showed evidence of a centripetal drainage pathway, separate from blood or lymphatic vessels, that can be visualized when a small amount of low molecular weight tracer is injected subcutaneously into a given region on skin of humans. In order to further characterize this interesting biological phenomenon, animal experiments are designed to elucidate histological and physiologic characteristics of these visualized pathways.

Methods

Multiple tracers are injected subcutaneously into an acupuncture point of KI3 to visualize centripetal pathways by magnetic resonance imaging or fluorescein photography in 85 healthy rabbits. The pathways are compared with venography and indirect lymphangiography. Fluid flow through the pathways is observed by methods of altering their hydrated state, hydrolyzing by different collagenases, and histology is elucidated by optical, fluorescein and electron microscopy.

Results

Histological and magnetic imaging examinations of these visualized pathways show they consist of perivenous loose connective tissues. As evidenced by examinations of tracers’ uptake, they appear to function as a draining pathway for free interstitial fluid. Fluorescein sodium from KI3 is found in the pathways of hind limbs and segments of the small intestines, partial pulmonary veins and results in pericardial effusion, suggesting systematical involvement of this perivenous pathway. The hydraulic conductivity of these pathways can be compromised by the collapse of their fiber-rich beds hydrolyzed by either of collagenase type I, III, IV or V.

Conclusions

The identification of pathways comprising perivenous loose connective tissues with a high hydraulic conductivity draining interstitial fluid in hind limbs of a mammal suggests a potential drainage system complementary to vascular circulations. These findings may provide new insights into a systematically distributed collagenous connective tissue with a circulatory function and their potential relevance to the nature of acupuncture meridians.     

Isolation of interstitial fluid from rat mammary tumors by a centrifugation method

Access to interstitial fluid is of fundamental importance to understand tumor transcapillary fluid balance, including the distribution of probes and therapeutic agents. Tumors were induced by gavage of 9,10-dimethyl-1,2-benzanthracene to rats, and fluid was isolated after anesthesia by exposing tissue to consecutive centrifugations from 27 to 6,800 g. The observed (51)Cr-EDTA (extracellular tracer) tissue fluid-to-plasma ratio obtained from whole tumor or from superficial tumor tissue by centrifugation at 27-424 g was not significantly different from 1.0 (0.92-0.99), suggesting an extracellular origin only. However, fluid collected from excised central tumor parts had a significantly lower ratio (0.66-0.77) for all imposed G forces, suggesting dilution by fluid deriving from a space unavailable for (51)Cr-EDTA. The colloid osmotic pressure in tumor fluid was generally higher than in fluid isolated from the subcutis, attributable to less selective capillaries and impaired lymphatic drainage in tumors. HPLC analysis of tumor fluid showed that low-molecular-weight macromolecules not present in arterial plasma were present in tumor fluid obtained by centrifugation and in venous blood draining the tumor, most likely representing proteins derived from tumor cells. We conclude that low-speed centrifugation may be a simple and reliable method to isolate interstitial fluid from tumors.     

A model of a radially expanding and contracting lymphangion

The lymphatic system is an extensive vascular network featuring valves and contractile walls that pump interstitial fluid and plasma proteins back to the main circulation. Immune function also relies on the lymphatic system's ability to transport white blood cells. Failure to drain and pump this excess fluid results in edema characterized by fluid retention and swelling of limbs. It is, therefore, important to understand the mechanisms of fluid transport and pumping of lymphatic vessels. Unfortunately, there are very few studies in this area, most of which assume Poiseuille flow conditions. In vivo observations reveal that these vessels contract strongly, with diameter changes of the order of magnitude of the diameter itself over a cycle that lasts typically 2-3s. The radial velocity of the contracting vessel is on the order of the axial fluid velocity, suggesting that modeling flow in these vessels with a Poiseuille model is inappropriate. In this paper, we describe a model of a radially expanding and contracting lymphatic vessel and investigate the validity of assuming Poiseuille flow to estimate wall shear stress, which is presumably important for lymphatic endothelial cell mechanotransduction. Three different wall motions, periodic sinusoidal, skewed sinusoidal and physiologic wall motions, were investigated with steady and unsteady parabolic inlet velocities. Despite high radial velocities resulting from the wall motion, wall shear stress values were within 4% of quasi-static Poiseuille values. Therefore, Poiseuille flow is valid for the estimation of wall shear stress for the majority of the lymphangion contractile cycle.     

A model for mechanics of primary lymphatic valves

Recent experimental evidence indicates that lymphatics have two valve systems, a set of primary valves in the wall of the endothelial cells of initial lymphatics and a secondary valve system in the lumen of the lymphatics. While the intralymphatic secondary valves are well described, no analysis of the primary valves is available. We propose a model for primary lymphatics valves at the junctions between lymphatic endothelial cells. The model consists of two overlapping endothelial extensions at a cell junction in the initial lymphatics. One cell extension is firmly attached to the adjacent connective tissue while the other cell extension is not attached to the interstitial collagen. It is free to bend into the lumen of the lymphatic when the lymphatic pressure falls below the adjacent interstitial fluid pressure. Thereby the cell junction opens a gap permitting entry of interstitial fluid into the lymphatic lumen. When the lymphatic fluid pressure rises above the adjacent interstitial fluid pressure, the endothelial extensions contact each other and the junction is closed preventing fluid reflow into the interstitial space. The model illustrates the mechanics of valve action and provides the first time a rational analysis of the mechanisms underlying fluid collection in the initial lymphatics and lymph transport in the microcirculation.     

Implications of Transvascular Fluid Exchange in Nonlinear,Biphasic Analyses of Flow-Controlled Infusion in Brain

A nonlinear, coupled biphasic-mass transport model that includes transvascular fluid exchange is proposed for flow-controlled infusions in brain tissue. The model accounts for geometric and material nonlinearities, a hydraulic conductivity dependent on deformation, and transvascular fluid exchange according to Starling’s law. The governing equations were implemented in a custom-written code assuming spherical symmetry and using an updated Lagrangian finite-element algorithm. Results of the model indicate that, using normal physiological values of vascular permeability, transvascular fluid exchange has negligible effects on tissue deformation, fluid pressure, and transport of the infused agent. As vascular permeability may be increased artificially through methods such as administering nitric oxide, a parametric study was conducted to determine how increased vascular permeability affects flow-controlled infusion. Increased vascular permeability reduced both tissue deformation and fluid pressure, possibly reducing damage to tissue adjacent to the infusion catheter. Furthermore, the loss of fluid to the vasculature resulted in a significantly increased interstitial fluid concentration but a modestly increased tissue concentration. From a clinical point of view, this increase in concentration could be beneficial if limited to levels below which toxicity would not occur. However, the modestly increased tissue concentration may make the increase in interstitial fluid concentration difficult to assess in vivo using co-infused radiolabeled agents.     

Toll-like receptor deficiency worsens inflammation and lymphedema after lymphatic injury

Mechanisms regulating lymphedema pathogenesis remain unknown. Recently, we have shown that lymphatic fluid stasis increases endogenous danger signal expression, and these molecules influence lymphatic repair (Zampbell JC, et al. Am J Physiol Cell Physiol 300: C1107-C1121, 2011). Endogenous danger signals activate Toll-like receptors (TLR) 2, 4, and 9 and induce homeostatic or harmful responses, depending on physiological context. The purpose of this study was to determine the role of TLRs in regulating tissue responses to lymphatic fluid stasis. A surgical model of lymphedema was used in which wild-type or TLR2, 4, or 9 knockout (KO) mice underwent tail lymphatic excision. Six weeks postoperatively, TLR KOs demonstrated markedly increased tail edema compared with wild-type animals (50-200% increase; P < 0.01), and this effect was most pronounced in TLR4 KOs (P < 0.01). TLR deficiency resulted in decreased interstitial and lymphatic transport, abnormal lymphatic architecture, and fewer capillary lymphatics (40-50% decrease; P < 0.001). Lymphedematous tissues of TLR KOs demonstrated increased leukocyte infiltration (P < 0.001 for TLR4 KOs), including higher numbers of infiltrating CD3+ cells (P < 0.05, TLR4 and TLR9 KO), yet decreased infiltrating F4/80+ macrophages (P < 0.05, all groups). Furthermore, analysis of isolated macrophages revealed twofold reductions in VEGF-C (P < 0.01) and LYVE-1 (P < 0.05) mRNA from TLR2-deficient animals. Finally, TLR deficiency was associated with increased collagen type I deposition and increased transforming growth factor-β1 expression (P < 0.01, TLR4 and TLR9 KO), contributing to dermal fibrosis. In conclusion, TLR deficiency worsens tissue responses to lymphatic fluid stasis and is associated with decreased lymphangiogenesis, increased fibrosis, and reduced macrophage infiltration. These findings suggest a role for innate immune responses, including TLR signaling, in lymphatic repair and lymphedema pathogenesis.     

A lumped parameter model of mechanically mediated acute and long-term adaptations of contractility and geometry in lymphatics for characterization of lymphedema

A primary purpose of the lymphatic system is to transport fluid from peripheral tissues to the central venous system in order to maintain tissue–fluid balance. Failure to perform this task results in lymphedema marked by swelling of the affected limb as well as geometric remodeling and reduced contractility of the affected lymphatic vessels. The mechanical environment has been implicated in the regulation of lymphatic contractility, but it is unknown how changes in the mechanical environment are related to loss of contractile function and remodeling of the tissue. The purpose of this paper was to introduce a new theoretical framework for acute and long-term adaptations of lymphatic vessels to changes in mechanical loading. This theoretical framework combines a simplified version of a published lumped parameter model for lymphangion function and lymph transport, a published microstructurally motivated constitutive model for the active and passive mechanical behavior of isolated rat thoracic ducts, and novel models for acute mechanically mediated vasoreactive adaptations and long-term volumetric growth to simulate changes in muscle contractility and geometry of a single isolated rat thoracic duct in response to a sustained elevation in afterload. The illustrative examples highlight the potential role of the mechanical environment in the acute maintenance of contractility and long-term geometric remodeling, presumably aimed at meeting fluid flow demands while also maintaining mechanical homeostasis. Results demonstrate that contractility may adapt in response to shear stress to meet fluid flow demands and show that pressure-induced long-term geometric remodeling may attenuate these adaptations and reduce fluid flow. The modeling framework and illustrative simulations help suggest relevant experiments that are necessary to accurately quantify and predict the acute and long-term adaptations of lymphangions to altered mechanical loading.     

Evaluation of gravimetric techniques to estimate the microvascular filtration coefficient

Microvascular permeability to water is characterized by the microvascular filtration coefficient (K(f)). Conventional gravimetric techniques to estimate K(f) rely on data obtained from either transient or steady-state increases in organ weight in response to increases in microvascular pressure. Both techniques result in considerably different estimates and neither account for interstitial fluid storage and lymphatic return. We therefore developed a theoretical framework to evaluate K(f) estimation techniques by 1) comparing conventional techniques to a novel technique that includes effects of interstitial fluid storage and lymphatic return, 2) evaluating the ability of conventional techniques to reproduce K(f) from simulated gravimetric data generated by a realistic interstitial fluid balance model, 3) analyzing new data collected from rat intestine, and 4) analyzing previously reported data. These approaches revealed that the steady-state gravimetric technique yields estimates that are not directly related to K(f) and are in some cases directly proportional to interstitial compliance. However, the transient gravimetric technique yields accurate estimates in some organs, because the typical experimental duration minimizes the effects of interstitial fluid storage and lymphatic return. Furthermore, our analytical framework reveals that the supposed requirement of tying off all draining lymphatic vessels for the transient technique is unnecessary. Finally, our numerical simulations indicate that our comprehensive technique accurately reproduces the value of K(f) in all organs, is not confounded by interstitial storage and lymphatic return, and provides corroboration of the estimate from the transient technique.     

Model study of the effects of interactions between systemic and peripheral circulation on interstitial fluid balance

Interstitial fluid balance is severely altered in microgravity, but the mechanisms underlying the fluid shift from lower to upper body are still partially unclear. A lumped parameter model of the arterial tree with active and non linear modulation of peripheral resistances and capillary fluid exchange was adopted to simulate the response of microcirculation to pulsatility and edema. Results suggest that myogenic regulation not only impinges on arteriolar radius, but it also indirectly affects interstitial fluid balance. Non linear dynamics of blood pressure (BP) and flow in capillary beds are influenced by systemic pulsatility, hinting that local activity is involved in the response to peripheral edema as well.