Spontaneous oscillations in a model for active control of microvessel diameters

Vasomotion in Retinal Arterioles Is Modified by Exercise and Flicker Stimulation

PURPOSE

Vasomotion is spontaneous oscillations in the diameter of resistance vessels with derived effects on blood flow, and it has been proposed that disturbances in vasomotion may be involved in retinal vascular disease. The purpose of this study was to investigate whether retinal vasomotion shows regional variation and is modified by activated autoregulation.

METHODS

Video recordings of the diameter of retinal arterioles previously obtained from 55 normal persons were subjected to Fourier analysis to characterize the frequencies and propagation of spontaneous diameter changes in retinal arterioles. The analyses were performed on peripapillary temporal retinal arterioles, on arteriolar branches toward the macular area and the retinal periphery, and were performed during rest, during an increase in the arterial blood pressure induced by isometric exercise, and during increased retinal metabolism induced by flickering light.

RESULTS

There was no propagation of diameter changes along the studied vascular segments. Isometric exercise constricted the arterioles significantly by (mean ± SD) 1.76% ± 3.56% (P = 0.02) and increased the power of diameter oscillations at very low frequencies (0.1-1.4 c/min). Flicker stimulation dilated the arterioles significantly by (mean ± SD) 5.10% ± 2.91% (P < 0.0001) and reduced the power of diameter oscillations at all but the very low frequencies (P < 0.006 for all comparisons). Flicker-induced dilation and changes in hydraulic conductance were lower in peripheral than in macular arterioles.

CONCLUSIONS

Retinal vasomotion in normal persons increases during increased arterial blood pressure and decreases during flicker stimulation. The findings may act as a basis for the study of vasomotion in retinal vascular disease.

Functional implications of microvascular heterogeneity for oxygen uptake and utilization

In the vascular system, an extensive network structure provides convective and diffusive transport of oxygen to tissue. In the microcirculation, parameters describing network structure, blood flow, and oxygen transport are highly heterogeneous. This heterogeneity can strongly affect oxygen supply and organ function, including reduced oxygen uptake in the lung and decreased oxygen delivery to tissue. The causes of heterogeneity can be classified as extrinsic or intrinsic. Extrinsic heterogeneity refers to variations in oxygen demand in the systemic circulation or oxygen supply in the lungs. Intrinsic heterogeneity refers to structural heterogeneity due to stochastic growth of blood vessels and variability in flow pathways due to geometric constraints, and resulting variations in blood flow and hematocrit. Mechanisms have evolved to compensate for heterogeneity and thereby improve oxygen uptake in the lung and delivery to tissue. These mechanisms, which involve long-term structural adaptation and short-term flow regulation, depend on upstream responses conducted along vessel walls, and work to redistribute flow and maintain blood and tissue oxygenation. Mathematically, the variance of a functional quantity such as oxygen delivery that depends on two or more heterogeneous variables can be reduced if one of the underlying variables is controlled by an appropriate compensatory mechanism. Ineffective regulatory mechanisms can result in poor oxygen delivery even in the presence of adequate overall tissue perfusion. Restoration of endothelial function, and specifically conducted responses, should be considered when addressing tissue hypoxemia and organ failure in clinical settings.

Using a Human Circulation Mathematical Model to Simulate the Effects of Hemodialysis and Therapeutic Hypothermia

Background: We developed a hemodynamic mathematical model of human circulation coupled to a virtual hemodialyzer. The model was used to explore mechanisms underlying our clinical observations involving hemodialysis. Methods: The model consists of whole body human circulation, baroreflex feedback control, and a hemodialyzer. Four model populations encompassing baseline, dialysed, therapeutic hypothermia treated, and simultaneous dialysed with hypothermia were generated. In all populations atrial fibrillation and renal failure as co-morbidities, and exercise as a treatment were simulated. Clinically relevant measurables were used to quantify the effects of each in silico experiment. Sensitivity analysis was used to uncover the most relevant parameters. Results: Relative to baseline, the modelled dialysis increased the population mean diastolic blood pressure by 5%, large vessel wall shear stress by 6%, and heart rate by 20%. Therapeutic hypothermia increased systolic blood pressure by 3%, reduced large vessel shear stress by 15%, and did not affect heart rate. Therapeutic hypothermia reduced wall shear stress by 15% in the aorta and 6% in the kidneys, suggesting a potential anti-inflammatory benefit. Therapeutic hypothermia reduced cardiac output under atrial fibrillation by 12% and under renal failure by 20%. Therapeutic hypothermia and exercise did not affect dialyser function, but increased water removal by approximately 40%. Conclusions: This study illuminates some mechanisms of the action of therapeutic hypothermia. It also suggests clinical measurables that may be used as surrogates to diagnose underlying diseases such as atrial fibrillation.

Cyclic Hypoxia: An Update on Its Characteristics, Methods to Measure It and Biological Implications in Cancer

Regions of hypoxia occur in most if not all solid cancers. Although the presence of tumor hypoxia is a common occurrence, the levels of hypoxia and proportion of the tumor that are hypoxic vary significantly. Importantly, even within tumors, oxygen levels fluctuate due to changes in red blood cell flux, vascular remodeling and thermoregulation. Together, this leads to cyclic or intermittent hypoxia. Tumor hypoxia predicts for poor patient outcome, in part due to increased resistance to all standard therapies. However, it is less clear how cyclic hypoxia impacts therapy response. Here, we discuss the causes of cyclic hypoxia and, importantly, which imaging modalities are best suited to detecting cyclic vs. chronic hypoxia. In addition, we provide a comparison of the biological response to chronic and cyclic hypoxia, including how the levels of reactive oxygen species and HIF-1 are likely impacted. Together, we highlight the importance of remembering that tumor hypoxia is not a static condition and that the fluctuations in oxygen levels have significant biological consequences.

Effects of impaired microvascular flow regulation on metabolism-perfusion matching and organ function

Impaired tissue oxygen delivery is a major cause of organ damage and failure in critically ill patients, which can occur even when systemic parameters, including cardiac output and arterial hemoglobin saturation, are close to normal. This review addresses oxygen transport mechanisms at the microcirculatory scale, and how hypoxia may occur in spite of adequate convective oxygen supply. The structure of the microcirculation is intrinsically heterogeneous, with wide variations in vessel diameters and flow pathway lengths, and consequently also in blood flow rates and oxygen levels. The dynamic processes of structural adaptation and flow regulation continually adjust microvessel diameters to compensate for heterogeneity, redistributing flow according to metabolic needs to ensure adequate tissue oxygenation. A key role in flow regulation is played by conducted responses, which are generated and propagated by endothelial cells and signal upstream arterioles to dilate in response to local hypoxia. Several pathophysiological conditions can impair local flow regulation, causing hypoxia and tissue damage leading to organ failure. Therapeutic measures targeted to systemic parameters may not address or may even worsen tissue oxygenation at the microvascular level. Restoration of tissue oxygenation in critically ill patients may depend on restoration of endothelial cell function, including conducted responses.

Rude mechanicals in brain haemodynamics: non-neural actors that influence blood flow

Fluctuations in blood oxygenation and flow are widely used to infer brain activity during resting-state functional magnetic resonance imaging (fMRI). However, there are strong systemic and vascular contributions to resting-state signals that are unrelated to ongoing neural activity. Importantly, these non-neural contributions to haemodynamic signals (or 'rude mechanicals') can be as large as or larger than the neurally evoked components. Here, we review the two broad classes of drivers of these signals. One is systemic and is tied to fluctuations in external drivers such as heart rate and breathing, and the robust autoregulatory mechanisms that try to maintain a constant milieu in the brain. The other class comprises local, active fluctuations that appear to be intrinsic to vascular tissue and are likely similar to active local fluctuations seen in vasculature all over the body. In this review, we describe these non-neural fluctuations and some of the tools developed to correct for them when interpreting fMRI recordings. However, we also emphasize the links between these vascular fluctuations and brain physiology and point to ways in which fMRI measurements can be used to exploit such links to gain valuable information about neurovascular health and about internal brain states. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

Electrical Propagation of Vasodilatory Signals in Capillary Networks

We have developed a computational model to study electrical propagation of vasodilatory signals and arteriolar regulation of blood flow depending on the oxygen tension and agonist distribution in the capillary network. The involving key parameters of endothelial cell-to-cell electrical conductivity and plasma membrane area per unit volume were calibrated with the experimental data on an isolated endothelial tube of mouse skeletal feeding arteries. We have estimated the oxygen saturation parameters in terms of erythrocyte ATP release from the data of a left anterior descending coronary blood perfusion of dog. Regarding the acetylcholine-induced upstream conduction, our model shows that spatially uniform superfusion of acetylcholine attenuates the electrical signal propagation, and blocking calcium-activated potassium channels suppresses that attenuation. On the other hand, a local infusion of acetylcholine induces enhanced electrical propagation that corresponds to physiological relevance. Integrating the electrophysiology of endothelial tube and the electrophysiology/mechanics of a lumped arteriole, we show mechanistically that endothelial purinergic oxygen sensing of ATP released from erythrocytes and local infusion of acetylcholine are individually effective to induce vasodilatory signals to regulate blood flow in arterioles. We have recapitulated the upstream vasomotion in arterioles from the elevated oxygen tension in the downstream capillary domain. This study is a foundation for characterizing effective pharmaceutical strategies for ascending vasodilation and oxygenation.

Spatial and temporal patterns of nitric oxide diffusion and degradation drive emergent cerebrovascular dynamics

Nitric oxide (NO) is a gaseous signaling molecule that plays an important role in neurovascular coupling. NO produced by neurons diffuses into the smooth muscle surrounding cerebral arterioles, driving vasodilation. However, the rate of NO degradation in hemoglobin is orders of magnitude higher than in brain tissue, though how this might impact NO signaling dynamics is not completely understood. We used simulations to investigate how the spatial and temporal patterns of NO generation and degradation impacted dilation of a penetrating arteriole in cortex. We found that the spatial location of NO production and the size of the vessel both played an important role in determining its responsiveness to NO. The much higher rate of NO degradation and scavenging of NO in the blood relative to the tissue drove emergent vascular dynamics. Large vasodilation events could be followed by post-stimulus constrictions driven by the increased degradation of NO by the blood, and vasomotion-like 0.1-0.3 Hz oscillations could also be generated. We found that these dynamics could be enhanced by elevation of free hemoglobin in the plasma, which occurs in diseases such as malaria and sickle cell anemia, or following blood transfusions. Finally, we show that changes in blood flow during hypoxia or hyperoxia could be explained by altered NO degradation in the parenchyma. Our simulations suggest that many common vascular dynamics may be emergent phenomena generated by NO degradation by the blood or parenchyma.

Correlation-and-Spectral Analysis of Time Series for Microcirculatory Parameters in Vessels of Symmetrical Organs

The correlation relationships between microcirculatory parameters were studied in vessels of symmetrical organs in humans. Parameters of microcirculation in healthy volunteers were measured by the method of laser Doppler flowmetry (LDF). LDF-signal transducers were placed symmetrically on the lower parts of the right and left shoulders (3 cm above the elbow bend). Correlations between various components of the amplitude-frequency spectrum for blood flow fluctuations (myogenic, neurogenic, respiratory, and cardiac) were evaluated on symmetrical sides. A significant correlation was found for the mean value of vascular perfusion on the right side with the neurogenic, myogenic, respiratory, and cardiac components on the same side. However, the mean value of vascular perfusion on the light side was interrelated only with the cardiac component. A positive correlation was revealed between the neurogenic and myogenic components on the left side. However, the neurogenic component on the right side was positively related to the myogenic component on the opposite side. Asymmetry of correlation relationships was also observed for the respiratory and cardiac components of study spectrum. Our results illustrate the specific regulation of microcirculatory blood flow in paired organs, which is associated with the existence of functional asymmetry. The physiological mechanisms for this asymmetry require further experimental and clinical studies.

A Multicellular Vascular Model of the Renal Myogenic Response

The myogenic response is a key autoregulatory mechanism in the mammalian kidney. Triggered by blood pressure perturbations, it is well established that the myogenic response is initiated in the renal afferent arteriole and mediated by alterations in muscle tone and vascular diameter that counterbalance hemodynamic perturbations. The entire process involves several subcellular, cellular, and vascular mechanisms whose interactions remain poorly understood. Here, we model and investigate the myogenic response of a multicellular segment of an afferent arteriole. Extending existing work, we focus on providing an accurate—but still computationally tractable—representation of the coupling among the involved levels. For individual muscle cells, we include detailed Ca2+ signaling, transmembrane transport of ions, kinetics of myosin light chain phosphorylation, and contraction mechanics. Intercellular interactions are mediated by gap junctions between muscle or endothelial cells. Additional interactions are mediated by hemodynamics. Simulations of time-independent pressure changes reveal regular vasoresponses throughout the model segment and stabilization of a physiological range of blood pressures (80–180 mmHg) in agreement with other modeling and experimental studies that assess steady autoregulation. Simulations of time-dependent perturbations reveal irregular vasoresponses and complex dynamics that may contribute to the complexity of dynamic autoregulation observed in vivo. The ability of the developed model to represent the myogenic response in a multiscale and realistic fashion, under feasible computational load, suggests that it can be incorporated as a key component into larger models of integrated renal hemodynamic regulation.

Bifurcation in Blood Oscillatory Rhythms for Patients with Ischemic Stroke: A Small Scale Clinical Trial using Laser Doppler Flowmetry and Computational Modeling of Vasomotion

We describe application of spectral analysis of laser Doppler flowmetry (LDF) signals to investigation of cerebrovascular haemodynamics in patients with post-acute ischemic stroke (AIS) and cerebrovascular insufficiency. LDF was performed from 3 to 7 days after the onset of AIS on forehead in the right and left supraorbital regions in patients. Analysis of LDF signals showed that perfusion in the microvasculature in AIS patients was lower than that in patients with cerebrovascular insufficiency. As a result of wavelet analysis of the LDF signals we obtained activation of the vasomotion in the frequency range of myogenic oscillation of 0.1 Hz and predominantly nutritive regime microcirculation after systemic thrombolytic therapy of the AIS patients. In case of significant stroke size, myogenic activity, and nutritive pattern microhaemodynamics were reduced, in some cases non-nutritive pattern and/or venular stasis was revealed. Wavelet analysis of the LDF signals also showed asymmetry in wavelet spectra of the LDF signals obtained in stroke-affected and unaffected hemispheres in the AIS patients. A mechanism underlying the observed asymmetry was analyzed by computational modeling of vasomotion developed in Arciero and Secomb (2012). We applied this model to describe relaxation oscillation of arteriole diameter which is forced by myogenic oscillation induced by synchronous calcium oscillation in vascular smooth muscle cells. Calculation showed that vasomotion frequency spectrum at the low-frequency range (0.01 Hz) is reciprocally modulated by myogenic oscillation (0.1 Hz) that correlates with experimental observation of inter-hemispheric variation in the LDF spectrum.

Bifurcation study of blood flow control in the kidney

Renal blood flow is maintained within a narrow window by a set of intrinsic autoregulatory mechanisms. Here, a mathematical model of renal hemodynamics control in the rat kidney is used to understand the interactions between two major renal autoregulatory mechanisms: the myogenic response and tubuloglomerular feedback. A bifurcation analysis of the model equations is performed to assess the effects of the delay and sensitivity of the feedback system and the time constants governing the response of vessel diameter and smooth muscle tone. The results of the bifurcation analysis are verified using numerical simulations of the full nonlinear model. Both the analytical and numerical results predict the generation of limit cycle oscillations under certain physiologically relevant conditions, as observed in vivo.

Functional sympatholysis and sympathetic escape in a theoretical model for blood flow regulation

A mathematical simulation of flow regulation in vascular networks is used to investigate the interaction between arteriolar vasoconstriction due to sympathetic nerve activity (SNA) and vasodilation due to increased oxygen demand. A network with 13 vessel segments in series is used, each segment representing a different size range of arterioles or venules. The network includes five actively regulating arteriolar segments with time-dependent diameters influenced by shear stress, wall tension, metabolic regulation, and SNA. Metabolic signals are assumed to be propagated upstream along vessel walls via a conducted response. The model exhibits functional sympatholysis, in which sympathetic vasoconstriction is partially abrogated by increases in metabolic demand, and sympathetic escape, in which SNA elicits an initial vasoconstriction followed by vasodilation. In accordance with experimental observations, these phenomena are more prominent in small arterioles than in larger arterioles when SNA is assumed to act equally on arterioles of all sizes. The results imply that a mechanism based on the competing effects on arteriolar tone of SNA and conducted metabolic signals can account for several observed characteristics of functional sympatholysis, including the different responses of large and small arterioles.

Capillary recruitment in a theoretical model for blood flow regulation in heterogeneous microvessel networks

In striated muscle, the number of capillaries containing moving red blood cells increases with increasing metabolic demand. This phenomenon, termed capillary recruitment, has long been recognized, but its mechanism has been unclear. Here, a theoretical model for metabolic blood flow regulation in a heterogeneous network is used to test the hypothesis that capillary recruitment occurs as a result of active control of arteriolar diameters, combined with unequal partition of hematocrit at diverging microvascular bifurcations. The network structure is derived from published observations of hamster cremaster muscle in control and dilated states. The model for modulation of arteriolar diameters includes length-tension characteristics of vascular smooth muscle and responses of smooth muscle tone to myogenic, shear-dependent, and metabolic stimuli. Blood flow is simulated including nonuniform hematocrit distribution. Convective and diffusive oxygen transport in the network is simulated. Oxygen-dependent metabolic signals are assumed to be conducted upstream from distal vessels to arterioles. With increasing oxygen demand, arterioles dilate, blood flow increases, and the numbers of flowing arterioles and capillaries, as defined by red blood cell flux above a small threshold value, increase. Unequal hematocrit partition at diverging bifurcations contributes to recruitment and enhances tissue oxygenation. The results imply that capillary recruitment, as observed in the hamster cremaster preparations, can occur as a consequence of local control of arteriolar tone and the resulting nonuniform changes in red blood cell fluxes, and provide an explanation for observations of sequential recruitment of individual capillaries in response to modulation of terminal arteriolar diameter.

Theoretical comparison of wall-derived and erythrocyte-derived mechanisms for metabolic flow regulation in heterogeneous microvascular networks

The objective of this study is to compare the effectiveness of metabolic signals derived from erythrocytes and derived from the vessel wall for regulating blood flow in heterogeneous microvascular networks. A theoretical model is used to simulate blood flow, mass transport, and vascular responses. The model accounts for myogenic, shear-dependent, and metabolic flow regulation. Metabolic signals are assumed to be propagated upstream along vessel walls via a conducted response. Arteriolar tone is assumed to depend on the conducted metabolic signal as well as local wall shear stress and wall tension, and arteriolar diameters are calculated based on vascular smooth muscle mechanics. The model shows that under certain conditions metabolic regulation based on wall-derived signals can be more effective in matching perfusion to local oxygen demand relative to regulation based on erythrocyte-derived signals, resulting in higher extraction and lower oxygen deficit. The lower effectiveness of the erythrocyte-derived signal is shown to result in part from the unequal partition of hematocrit at diverging bifurcations, such that low-flow vessels tend to receive a reduced hematocrit and thereby experience a reduced erythrocyte-derived metabolic signal. The model simulations predict that metabolic signals independent of erythrocytes may play an important role in local metabolic regulation of vascular tone and flow distribution in heterogeneous microvessel networks.