Simulation study of autoregulation responses of peripheral circulation to systemic pulsatility
© Aletti et al; licensee BioMed Central Ltd. 2009
Received: 01 March 2009
Accepted: 24 July 2009
Published: 24 July 2009
This simulation study investigated potential modulations of total peripheral resistance (TPR), due to distributed peripheral vascular activity, by means of a lumped model of the arterial tree and a non linear model of microcirculation, inclusive of local controls of blood flow and tissue-capillary fluid exchange.
Numerical simulations of circulation were carried out to compute TPR under different conditions of blood flow pulsatility, and to extract the pressure-flow characteristics of the cardiovascular system. Simulations showed that TPR seen by the large arteries was increased in absence of pulsatility, while it decreased with an augmented harmonic content. This is a typically non linear effect due to the contribution of active, non linear autoregulation of the peripheral microvascular beds, which also generated a nonlinear relationship between arterial blood pressure and cardiac output.
This simulation study, though focused on a simple effect attaining TPR modulation due to pulsatility, suggests that non-linear autoregulation mechanisms cannot be overlooked while studying the integrated behavior of the global cardiovascular system, including the arterial tree and the peripheral vascular bed.
The study of the interactions between the control of hemodynamics in large arteries and in the peripheral districts of the arterial tree is only marginally addressed in the literature. In the investigation of global control of circulation, cardiovascular (CV) variability (CVV) research has mainly focused on heart rate (HR) variability (HRV), baroreflex and sympathetic control of vascular tone [1, 2]. Several studies concerning peripheral autoregulation have instead shown that a local activity, independent from neural drives, is able to regulate blood flow to tissues: for instance, modulations of peripheral resistances in response to mechanical and metabolic stimuli have been repeatedly observed and modeled [3–8].
Although some studies have taken local mechanisms into account when discussing the genesis and sustenance of slow arterial pressure waves [9–11], to our knowledge no systematic effort for a comprehensive interpretation of the mutual interactions between global control of systemic circulation and local autoregulation has been proposed yet. It may be ascribed to the experimental limitations affecting the in vivo study of microcirculation and the consequently difficult validation of the hypotheses regarding these phenomena, or to the uncertainty about which CV variables carry information relative to the activity of the peripheries of the arterial tree. Neural and humoral control mechanisms of blood flow and pressure act in parallel on HR, peripheral resistances, arterial compliances, venous return, and heart contractility. Therefore, the unfolding of the specific roles of each control represents a challenging task still largely overlooked. Consequently, our focus was not on the arterial tree nor on the peripheral autoregulation, per se, rather on the integration of a complex system.
It is well known that TPR is a major determinant of ABP variability ; hence, in the hypothesis of an active response of the peripheral autoregulation to pressure wave features, a mutual interaction between hemodynamics in large vessels and local control of blood flow follows.
To our knowledge, only two works by Aljuri and Cohen [13, 14] tried to identify autoregulation induced components of TPR oscillations and, in doing so, they introduced a critical interpretation of TPR as seen in a classical, Windkessel like context. In such a perspective, TPR is assumed as a constant constitutive parameter of the arterial tree [15, 16]. However, TPR results from the summation of all peripheral resistances of the system and should therefore reflect the local modulations, which are distributed in the peripheries of the arterial tree.
In this paper, we adopt a model for the simulation of hemodynamics in large vessels and autoregulation of microcirculation, in order to formulate hypotheses about the effect of peripheral responses on systemic hemodynamics, which can be hardly investigated through experiments. In particular, the modulation of TPR by local dynamics was simulated in function of the pulsatility of circulation, identified as a powerful drive for autoregulation, with a special regard to the comparison with non pulsatile circulation, which is the typical clinical condition characterizing heart surgery requiring non pulsatile cardiopulmonary bypass (NP-CPB).
Global, neural control mechanisms were not included in this model because our goal consisted in observing the potential effect of local non linear activity on systemic hemodynamics, which implied to fix the tone of sympathetic outflow and temporarily exclude global control mechanisms, though physiologically important.
Model of the arterial tree and microcirculation
Each segment of the model was described in terms of ordinary differential equations, representing mass and energy conservation .
Peripheral resistance characterization represents the novelty with respect to classical models, since they are not constant resistive elements, as considered either with lumped parameter approaches [15, 16] or structured tree approaches . On the contrary, our model includes the analytical description of the passive elastic relaxation of peripheral vessel walls, of the active mechanisms of local regulation of blood flow, and of the filtration of fluids through the capillary walls.
where R is the resistance of the considered element, μ the blood dynamic viscosity, l the length of the vessel, r the vessel radius and N the number of vessels in parallel.
The number of arterioles, capillaries and venules varies from network to network while geometrical parameters of single vessels are the same for all districts, in accordance with Lanzarone et al. . Fahareus effect was taken into account for small diameter vessels, according to the in vivo formulation of Pries et al .
The capacitor Cp represents the compliance of terminal arteries, which cuts off the high frequencies typical of blood flow in large vessels, in order to obtain low pulsatility, typical of blood circulation in microvessels. Values of the compliances were computed on the basis of the time constant RC of each network .
The flow generator Qf describes the filtration processes through the capillary wall. The pressure generator pv lumps the venous system downstream peripheral circulation and its value was set to 8 mmHg in all the districts.
Myogenic regulation of pre-capillary resistance
The myogenic regulation controls the arteriolar radius, modulating it in function of the pressure drop on arterioles [3, 8, 17]. Its action consists in maintaining the flow to a district of the system by regulating the arteriolar tone, and translates into an active regulation characterized by a positive feedback (typical of autoregulation), of vascular tone, lumen and resistance, following changes in pressure. That is to say, a constriction is elicited by a raise in pressure: this raises arteriolar resistance too, causing in turn a further increase in pressure.
where pc is the blood pressure at the inlet of the peripheral network, Ra is the arteriolar resistance and Qa is the blood flow entering the peripheral network.
The resulting expression for the radius is the sum of rpass and ract. All the constants in eqs.3–8 are in accordance with Lanzarone et al. .
Metabolic regulation of post-capillary resistance
The metabolic regulation modulates the peripheral vascular tone evaluating the specific O2 consumption by the district and comparing it to the actual oxygen need of the tissue. In case of a deficit of oxygen, the metabolic regulation induces vasodilation, so that the resistance opposed to the flow in the network can be reduced. The implementation of the metabolic control in the model followed [5, 17] in attributing its effects only to postcapillary vessels.
where Qa is the flow entering the peripheral network through the resistance Ra, Cons the physiological total oxygen consumption of the district, assumed constant in time, and CAO2 the oxygen concentration in the arterial blood entering the district .
The possible role of arteriolar oxygenation , though not denied, was not included in this model.
Fluid balance across capillary walls
The microcirculation model also describes the interactions between the vascular compartment and the interstitial space. The filtration through capillary walls is modeled through a flow generator Qf (fig. 2), representing the filtration flow.
where J is the filtration flow through the capillary wall, and Lp stands for the permeability index of capillary membranes, Δp(x) is the hydraulic pressure gradient between the blood pressure in capillaries and the interstitial pressure outside microvessels in function of the variable x, which accounts for the length of the microvessel. Δπ(x) is the osmotic pressure gradient between the blood pressure in capillaries and the interstitial pressure in function of the variable x.
The capillary radius is subject to the effects of the filtration process through capillary walls . Therefore, this aspect of the analytical model of peripheral districts actually introduces one more controller of the peripheral conductance, specifically of the capillary conductance.
Simulations: parameters of blood and circulation, inputs to the system, numerical integration
Numerical simulations of circulation were carried out to estimate TPR, to derive the pressure-flow characteristic curve of the arterial tree in presence and absence of local controls of blood flow, and to investigate local dynamics in response to systemic pulsatility, which explain TPR responses.
As regards blood parameters, physiological values were imposed for temperature (37°C), hematocrit (Ht) (45%) and concentration of blood proteins (7.5 g/100 ml of blood). Based on these values, blood density, dynamic viscosity and oncotic pressure were determined according to .
The arterial tree was solicited with a set of different flow inputs, imposed at the inlet of the system, in order to investigate the different responses in terms of TPR and to analyze how these systemic differences are due to local activity.
The inputs were either pulsatile flow waves and a non pulsatile one. The system was solicited by different heart flow waveforms at a fixed CO (5 l/min). A steady state pressure input was conversely applied in order to characterize pressure/flow relationships, with and without peripheral control.
where Ts is the systole duration in s, and HR, which was set at 1.25 Hz (75 bpm).
The idea behind the adoption of this set of input waves was to look for a relation between different degrees of pulsatility and the overall response of the arterial tree, reflected in mean TPR changes.
CO and HR were fixed so that the effect of pulsatility on the system could be evaluated under equal amplitude of stroke volume in any cardiac cycle.
MAP was computed as the mean pressure in the segment 1a of the model. A more accurate definition of TPR includes the difference between MAP and central venous pressure (CVP) in the numerator of eq.16. Still, since our model does not describe the venous side of circulation but only pictures the arterial side of circulation, TPR definition did not take CVP into account.
TPR of the arterial tree under different cardiac output waveforms
Pressure-flow characteristic of the arterial tree
The pressure-flow characteristic curve of the arterial tree was linear in absence of the effect of local controls. This basically implied a constant value of TPR, as typically assumed by classical lumped parameter models, in spite of changes in the steady value of pressure applied at the inlet of the arterial tree. On the other hand, the curve obtained with active local control showed a non linear behavior and lower steepness. The intersection between the two curves represents the physiological working point of the system, where the metabolic demand is balanced by CO.
Local responses to systemic pulsatility
Mean values of pressures and resistances in a peripheral network as a function of pulsatility and controls.
The sensitivity of TPR to aortic flow shapes and harmonic content (fig. 4) appeared to depend on the superposition of the active peripheral responses in all districts.
In our simulated results, the propagation of pressure waves along the arterial tree and the degree of residual pulsatility at the inlet of peripheries was affected by the flow waveform at the inlet of the arterial tree. As a consequence, a constant TPR, as hypothesized by a linear model, should be a theoretical concept bound to specific conditions. Changes of an apparent TPR are present and supposed to derive from the complex peripheral modulation, which is responsive not only to mean pressure and flow values, but also to other dynamical factors such as the shape of the systemic input waveform. Any type of pulsatile wave was able to induce a lower TPR than without systemic pulsatility.
In other terms, the global system under periodical regime sets around mean values (alias, DC components) which are influenced by the forced oscillations and their harmonic content, which is a typically non-linear behavior.
Translated into the field of NP-CPB, this may support theories and experimental work proposing the use of pulsatile devices in heart-lung machines, in order to reduce and possibly minimize the invasiveness of the procedures. Several studies showed that the presence of a pulsatile flow CPB proves to be beneficial in terms of better district and organ perfusion, compared to the continuous flow CPB , although results in this sense are not conclusive and still rather controversial. Recent studies confirmed a different behavior, underlining that the better perfusion can be obtained only in the presence of an adequate reproduction of physiological pulsatility, in terms of left ventricle flow waveform [24, 25]. In some cases , experimental results have also specifically shown a better preservation of microcirculation, which could be a validating element for our simulated results. However, the problem of vascular responses to a pulsatile or non pulsatile flow CPB, which was also studied from different perspectives with respect to classical assumptions on CV oscillations [26–28], is still widely debated and, even neglecting the potential effects which pulsatility may have on baroceptors stimulation or other global regulatory mechanisms, understanding the response of local activity may shed light on the combined interactions of pulsatility, hemodilution, tissue-capillary fluid exchange and vasomotion.
The pressure-flow characteristic (fig. 5) showed an expected behavior of controls stabilizing flow against mean pressure changes and largely reducing slope, coherently with the modeling description of autoregulation in .
This result was also considered a validation for our control equations: as previously pointed out, local controls of blood flow always aim at preserving blood flow to tissues by varying vessels resistance in response to changes in blood pressure. Thus, it appears that autoregulation itself, considering an unchanging sympathetic tone on vessels, may produce modulations of peripheral resistances. Therefore, TPR should not be considered constant as in the hydraulic adaptation of Ohm's law (eq.15), while, from the point of view of autoregulation mechanisms, blood flow is. TPR and ABP modulations mutually affect one each other in virtue of the controls of blood flow and arterial pressure.
Moreover, such a result appeared consistent with experimental observations , that demonstrated that autoregulation does modulate regional and systemic vascular resistance in the absence of reflex neurohormonal control of the circulation.
Since our simulations were carried out under the same value of mean CO, the metabolic needs of all peripheral districts were considered matched by their flows. Thus, the myogenic control appeared to be the main responsible for local responses sensitive to flow dynamics, modulating peripheral resistances in response to pressure fluctuations. The time courses of pc, pinm and ppI (fig. 6) clearly showed that the variable ppI (eq.2), which drives the myogenic regulation, basically had the same dynamic characterization of the arterial pressure at the inlet of the peripheral network because pressure at the inlet of the capillary bed was fundamentally non pulsatile, as an effect of the role of the capacitor Cp. However, the lumped representation adopted did not permit to exclude residual pulsatility in the initial segment of the capillary bed, according to compliance/resistance distribution.
The simulated response of the myogenic control to systemic pulsatility (table 1) showed a lower arteriolar resistance in order to keep the physiological perfusion of peripheries unaltered: hence, the precapillary pressure drop reduced as well and the working pressures in the whole peripheral district are lower than under a non pulsatile regimen of perfusion. The lower value of resistance obtained in the periphery p16s in presence of pulsatility confirmed that the macroscopic effect shown by TPR in fig. 4 is a reflection of phenomena occurring at the peripheral level.
Summarizing, the core of the discussion of the results of the simulated studies that were presented here can be synthesized in a few key points:
1. Results of our simulations hinted that the myogenic mechanism and the peripheral compliance, which buffers the systemic pulsatility at the inlet of any peripheral network, may represent the link between systemic and peripheral circulation.
2. The combined effect of non linearity due to controls and to Hagen-Poiseuille law (which, although not the object of this discussion, has to be considered as a further element of non linearity because it states that the radius of a vessel and its resistance are related through the inverse fourth power of the radius) results in a lower peripheral resistance under pulsatility.
3. Pulsatility is able to elicit rhythmical oscillations of the arteriolar walls, due to the action of the myogenic mechanisms, whose net effect consist in a decrease of arteriolar resistances.
The rationale of this study identified pulsatility as a strong drive in peripheral responses. The idea that autoregulation is responsive to cardiac output (CO) dynamic properties  and is characterized by a positive feedback response to variations of pressure, could also add material, in a different perspective, to the discussion regarding the concept of "positive feedback" proposed and theorized by Malliani , referring to positive feedback autonomic responses, and constituting a distinct contribution to arterial control from negative feedback reflexes.
The presented simulated results may also encourage efforts in experimental microcirculation physiology towards the integration of results deriving from the analysis of isolated microvessels and the investigation of the whole circulatory effects. Both experimental and modeling developments may also cast new light on the contribution of other dynamic components, not addressed here, but which could be easily included in the proposed description of periphery to system interactions. For instance, locally induced vasomotion , associated with intracellular Ca2+ modulation , may potentially contribute to the modulation of apparent TPR in a non-linear portrait, thus affecting important systemic features such as the heart workload.
Nonetheless, it would be very important to envision possible fields of applications for simulation studies of very significant physiological mechanisms that result deeply affected by relevant pathological conditions or induced altered conditions. In this paper, heart surgery requiring CPB was mentioned, and it certainly represents a major condition for which this type of modeling may prove helpful.
This study was supported by the Italian Space Agency (ASI, Agenzia Spaziale Italiana) under the DCMC Project and by the HeartCycle ICT-216695 Grant of the European Union.
- Akselrod S, Gordon D, Ubel FA, Shannon DC, Berger AC, Cohen RJ: Power spectrum analysis of heart rate fluctuation: a quantitative probe of beat-to-beat cardiovascular control. Science. 1981, 213: 220-222. 10.1126/science.6166045.View ArticleADSGoogle Scholar
- Pagani M, Lombardi F, Guzzetti S, Rimoldi O, Furlan R, Pizzinelli P, Sandrone G, Malfatto G, Dell'Orto S, Piccaluga E, Turiel M, Baselli G, Cerutti S, Malliani A: Power spectral analysis of heart rate and arterial pressure variabilities as a marker of sympatho-vagal interaction in man and conscious dog. Circ Res. 1986, 59: 178-193.View ArticleGoogle Scholar
- Davis MJ: Myogenic response gradient in an arteriolar network. Am J Physiol. 1993, 264: H2168-2179.Google Scholar
- Kim MH, Harris NR, Korzick DH, Tarbell JM: Control of the arteriolar myogenic response by transvascular fluid filtration. Microvasc Res. 2004, 68: 30-37. 10.1016/j.mvr.2004.03.002.View ArticleGoogle Scholar
- Groebe K: Precapillary servo control of blood pressure and postcapillary adjustment of flow to tissue metabolic status; a new paradigm for local perfusion regulation. Circulation. 1996, 94: 1876-1885.View ArticleGoogle Scholar
- Cornelissen AJM, Dankelman J, Van Bavel E, Spaan JAE: Balance between myogenic, flow-dependent, and metabolic flow control in coronary arterial tree: a model study. Am J Physiol Heart Circ Physiol. 2002, 282: H2224-2237.View ArticleGoogle Scholar
- Cavalcanti S, Ursino M: Chaotic oscillations in microvessel arterial networks. Ann Biomed Eng. 1996, 24: 37-47. 10.1007/BF02770993.View ArticleGoogle Scholar
- Ursino M, Fabbri G, Belardinelli E: A mathematical analysis of vasomotion in the peripheral vascular bed. Cardioscience. 1992, 3: 13-25.Google Scholar
- Baselli G, Porta A, Pagani M: Coupling arterial Windkessel with peripheral vasomotion: modeling the effects on low-frequency oscillations. IEEE Trans Biomed Eng. 2006, 53: 53-64. 10.1109/TBME.2005.859787.View ArticleGoogle Scholar
- Julien C: The enigma of Mayer waves: facts and models. Cardiovasc Res. 2006, 70: 12-21. 10.1016/j.cardiores.2005.11.008.View ArticleGoogle Scholar
- Killip T: Oscillation of blood flow and vascular resistance during Mayer waves. Circ Res. 1962, 11: 987-993.View ArticleGoogle Scholar
- O'Leary DS, Woodbury DJ: Role of cardiac output in mediating arterial blood pressure oscillations. Am J Physiol. 1996, 271: 641-646.Google Scholar
- Aljuri N, Cohen RJ: Theoretical considerations in the dynamic closed-loop baroreflex and autoregulatory control of total peripheral resistance. Am J Physiol Heart Circ Physiol. 2004, 287: H2252-2273. 10.1152/ajpheart.00489.2003.View ArticleGoogle Scholar
- Aljuri N, Marini R, Cohen RJ: Test of dynamic closed-loop baroreflex and autoregulatory control of total peripheral resistance in intact and conscious sheep. Am J Physiol Heart Circ Physiol. 2004, 287: H2274-2286. 10.1152/ajpheart.00490.2003.View ArticleGoogle Scholar
- Burattini R, Gnudi G: Computer identification of models for the arterial tree input impedance: comparison between two new models and first experimental results. Med Biol Eng Comput. 1982, 20: 134-144. 10.1007/BF02441348.View ArticleGoogle Scholar
- Westerhof N, Bosman F, De Vries CJ, Noordergraaf A: Analogue studies of the human systemic arterial tree. J Biomech. 1969, 2: 121-143. 10.1016/0021-9290(69)90024-4.View ArticleGoogle Scholar
- Lanzarone E, Liani P, Baselli G, Costantino ML: Model of arterial tree and peripheral control for the study of physiological and assisted circulation. Med Eng Phys. 2007, 29: 542-555. 10.1016/j.medengphy.2006.08.004.View ArticleGoogle Scholar
- Olufsen M, Peskin CS, Kim WY, Pedersen EM, Nadim A, Larsen J: Numerical simulation and experimental validation of blood flow in arteries with structured-tree outflow conditions. Ann Biomed Eng. 2000, 28: 1281-1299. 10.1114/1.1326031.View ArticleGoogle Scholar
- Pries AR, Secomb TW, Gaehtgens P: Biophysical aspects of blood flow in the microvasculature. Cardiovasc Res. 1996, 32: 654-667.View ArticleGoogle Scholar
- Bird RB, Stewart WE, Lightfoot EN: Transport phenomena. Revised Second Edition ed. 2007, John Wiley & Sons, New York, USAGoogle Scholar
- Swanson WM, Clark RE: Cardiovascular system simulation requirements, design and performance. J Bioeng. 1997, 1: 121-133.Google Scholar
- Katz LN, Feil HS: Clinical observation on the dynamics of ventricular systole; I: Auricular fibrillation. Arch Intern Med. 1923, 32: 672-692.View ArticleGoogle Scholar
- Fumero R, Montevecchi FM, Scuri S, Carrara B, Gamba A, Parenzan L: Clinical experience with a new pulsatile pump for infant and pediatric cardiopulmonary bypass. Int J Artif Organs. 1989, 12: 314-320.Google Scholar
- Ündar A, Rosenberg G, Myers JL: Major factors in the controversy of pulsatile versus nonpulsatile flow during acute and chronic cardiac support. ASAIO J. 2005, 51: 173-175. 10.1097/01.MAT.0000161944.20233.40.View ArticleGoogle Scholar
- Ji B, Ündar A: An evaluation of the benefits of pulsatile versus nonpulsatile perfusion during cardiopulmonary bypass procedures in pediatric and adult cardiac patients. ASAIO J. 2006, 52: 357-361. 10.1097/01.mat.0000225266.80021.9b.View ArticleGoogle Scholar
- Podgoreanu MV, Stout RG, El-Moalem HE, Silverman DG: Synchronous rhythmical vasomotion in the human cutaneous microvasculature during nonpulsatile cardiopulmonary bypass. Anesthesiology. 2002, 97: 1110-1117. 10.1097/00000542-200211000-00012.View ArticleGoogle Scholar
- Silverman DG, Stout RG, Lee FA, Ferneini EM: Detection and characterization of cholinergic oscillatory control in the forehead microvasculature in response to systemic alpha-agonist infusion in healthy volunteers. Microvasc Res. 2001, 61: 144-147. 10.1006/mvre.2000.2283.View ArticleGoogle Scholar
- Silverman DG, Stout RG: Distinction between atropine-sensitive control of microvascular and cardiac oscillatory activity. Microvasc Res. 2002, 63: 196-208. 10.1006/mvre.2001.2385.View ArticleGoogle Scholar
- Burattini R, Borgdorff P, Westerhof N: The baroreflex is counteracted by autoregulation, thereby preventing circulatory instability. Exp Physiol. 2004, 89: 397-405. 10.1113/expphysiol.2003.027094.View ArticleGoogle Scholar
- Létienne R, Barrès C, Cerutti C, Julien C: Short-term haemodynamic variability in the conscious areflexic rat. J Physiol. 1998, 506: 263-274. 10.1111/j.1469-7793.1998.263bx.x.View ArticleGoogle Scholar
- Malliani A: Principles of cardiovascular neural regulation in health and disease. 2000, Kluwer Academic, Boston, USAView ArticleGoogle Scholar
- Nilsson H, Aalkjaer C: Vasomotion: mechanisms and physiological importance. Mol Interv. 2003, 3: 79-89. 10.1124/mi.3.2.79.View ArticleGoogle Scholar
- Peng H, Matchkov V, Ivarsen A, Aalkjaer C, Nilsson H: Hypothesis for the initiation of vasomotion. Circ Res. 2001, 88: 810-815. 10.1161/hh0801.089603.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.