Brian E. Carlson

Theoretical model of blood flow autoregulation: roles of myogenic, shear-dependent, and metabolic responses

The autoregulation of blood flow, the maintenance of almost constant blood flow in the face of variations in arterial pressure, is characteristic of many tissue types. Here, contributions to the autoregulation of pressure-dependent, shear stress-dependent, and metabolic vasoactive responses are analyzed using a theoretical model. Seven segments, connected in series, represent classes of vessels: arteries, large arterioles, small arterioles, capillaries, small venules, large venules, and veins. The large and small arterioles respond actively to local changes in pressure and wall shear stress and to the downstream metabolic state communicated via conducted responses. All other segments are considered fixed resistances. The myogenic, shear-dependent, and metabolic responses of the arteriolar segments are represented by a theoretical model based on experimental data from isolated vessels. To assess autoregulation, the predicted flow at an arterial pressure of 130 mmHg is compared with that at 80 mmHg. If the degree of vascular smooth muscle activation is held constant at 0.5, there is a fivefold increase in blood flow. When myogenic variation of tone is included, flow increases by a factor of 1.66 over the same pressure range, indicating weak autoregulation. The inclusion of both myogenic and shear-dependent responses results in an increase in flow by a factor of 2.43. A further addition of the metabolic response produces strong autoregulation with flow increasing by a factor of 1.18 and gives results consistent with experimental observation. The model results indicate that the combined effects of myogenic and metabolic regulation overcome the vasodilatory effect of the shear response and lead to the autoregulation of blood flow.

Theoretical model of metabolic blood flow regulation: roles of ATP release by red blood cells and conducted responses

A proposed mechanism for metabolic flow regulation involves the saturation-dependent release of ATP by red blood cells, which triggers an upstream conducted response signal and arteriolar vasodilation. To analyze this mechanism, a theoretical model is used to simulate the variation of oxygen and ATP levels along a flow pathway of seven representative segments, including two vasoactive arteriolar segments. The conducted response signal is defined by integrating the ATP concentration along the vascular pathway, assuming exponential decay of the signal in the upstream direction with a length constant of approximately 1 cm. Arteriolar tone depends on the conducted metabolic signal and on local wall shear stress and wall tension. Arteriolar diameters are calculated based on vascular smooth muscle mechanics. The model predicts that conducted responses stimulated by ATP release in venules and propagated to arterioles can account for increases in perfusion in response to increased oxygen demand that are consistent with experimental findings at low to moderate oxygen consumption rates. Myogenic and shear-dependent responses are found to act in opposition to this mechanism of metabolic flow regulation.

A theoretical model for the myogenic response based on the length-tension characteristics of vascular smooth muscle

Objective: A theoretical model is developed to describe the myogenic response of resistance vessels to changes in intravascular pressure, based on a consideration of the active and passive length–tension characteristics of vascular smooth muscle (VSM). The dependence of model parameters on vessel diameter is examined.

Candidate Urinary Biomarker Discovery in Ureteropelvic Junction Obstruction: A Proteomic Approach

PURPOSE Ureteropelvic junction obstruction may either worsen and require surgery, improve or remain stable. It may take upward of 3 years for the natural history to unfold. Urinary proteome analysis using capillary electrophoresis mass spectrometry has been shown to differentiate between normal infants and those with ureteropelvic junction obstruction. We sought to confirm these findings using liquid chromatography/nano-spray mass spectrometry to examine the urinary proteome in patients with unilateral grade IV ureteropelvic junction obstruction compared to age matched healthy infants. MATERIALS AND METHODS Urine specimens were obtained from 21 healthy infants with normal maternal/fetal ultrasound and 25 infants with grade IV unilateral ureteropelvic junction obstruction. Specimens were prepared using standard methods and subjected to liquid chromatography/tandem mass spectrometry analysis. Normalized data were annotated using the IPA(R) knowledge platform. RESULTS There were 31 proteins significantly different in their level of abundance at 1 to 6 months, and 18 at 7 to 12 months compared to age matched controls. These proteins clustered into major functional networks. All of the biomarkers previously reported in clinical studies of ureteropelvic junction obstruction were observed with the notable exception of transforming growth factor-beta1. CONCLUSIONS These results confirm the presence of significant differences in the urinary proteome in unilateral ureteropelvic junction obstruction compared to age matched normal individuals. This study adds new information about levels of abundance of specific proteins and peptides in ureteropelvic junction obstruction, which may allow for better classification of disease subgroups and help to establish improved indications for the early selection of surgical candidates based on urinary protein biomarkers.

Multiscale modeling and data integration in the Virtual Physiological Rat Project

It has become increasingly evident that the descriptions of many complex diseases are only possible by taking into account multiple influences at different physiological scales. To do this with computational models often requires the integration of several models that have overlapping scales (genes to molecules, molecules to cells, cells to tissues). The Virtual Physiological Rat (VPR) Project, a National Institute of General Medical Sciences (NIGMS) funded National Center of Systems Biology, is tasked with mechanistically describing several complex diseases and is therefore identifying methods to facilitate the process of model integration across physiological scales. In addition, the VPR has a considerable experimental component and the resultant data must be integrated into these composite multiscale models and made available to the research community. A perspective of the current state of the art in model integration and sharing along with archiving of experimental data will be presented here in the context of multiscale physiological models. It was found that current ontological, model and data repository resources and integrative software tools are sufficient to create composite models from separate existing models and the example composite model developed here exhibits emergent behavior not predicted by the separate models.

Integration of multi-scale biosimulation models via light-weight semantics.

Currently, biosimulation researchers use a variety of computational environments and languages to model biological processes. Ideally, researchers should be able to semiautomatically merge models to more effectively build larger, multi-scale models. However, current modeling methods do not capture the underlying semantics of these models sufficiently to support this type of model construction. In this paper, we both propose a general approach to solve this problem, and we provide a specific example that demonstrates the benefits of our methodology. In particular, we describe three biosimulation models: (1) a cardio-vascular fluid dynamics model, (2) a model of heart rate regulation via baroreceptor control, and (3) a sub-cellular-level model of the arteriolar smooth muscle. Within a light-weight ontological framework, we leverage reference ontologies to match concepts across models. The light-weight ontology then helps us combine our three models into a merged model that can answer questions beyond the scope of any single model.

Effect of P2X4 and P2X7 receptor antagonism on the pressure diuresis relationship in rats

Reduced glomerular filtration, hypertension and renal microvascular injury are hallmarks of chronic kidney disease, which has a global prevalence of ~10%. We have shown previously that the Fischer (F344) rat has lower GFR than the Lewis rat, and is more susceptible to renal injury induced by hypertension. In the early stages this injury is limited to the pre-glomerular vasculature. We hypothesized that poor renal hemodynamic function and vulnerability to vascular injury are causally linked and genetically determined. In the present study, normotensive F344 rats had a blunted pressure diuresis relationship, compared with Lewis rats. A kidney microarray was then interrogated using the Endeavour enrichment tool to rank candidate genes for impaired blood pressure control. Two novel candidate genes, P2rx7 and P2rx4, were identified, having a 7− and 3− fold increased expression in F344 rats. Immunohistochemistry localized P2X4 and P2X7 receptor expression to the endothelium of the pre-glomerular vasculature. Expression of both receptors was also found in the renal tubule; however there was no difference in expression profile between strains. Brilliant Blue G (BBG), a relatively selective P2X7 antagonist suitable for use in vivo, was administered to both rat strains. In Lewis rats, BBG had no effect on blood pressure, but increased renal vascular resistance, consistent with inhibition of some basal vasodilatory tone. In F344 rats BBG caused a significant reduction in blood pressure and a decrease in renal vascular resistance, suggesting that P2X7 receptor activation may enhance vasoconstrictor tone in this rat strain. BBG also reduced the pressure diuresis threshold in F344 rats, but did not alter its slope. These preliminary findings suggest a physiological and potential pathophysiological role for P2X7 in controlling renal and/or systemic vascular function, which could in turn affect susceptibility to hypertension-related kidney damage.

A Reappraisal of How to Build Modular, Reusable Models of Biological Systems

61Department of Bioengineering, University of Washington, Seattle, Washington, United States of America, 2Auckland Bioengineering Institute, University of Auckland,Auckland, New Zealand, 3Department of Physiology, Medical College of Wisconsin, Milwaukee, Wisconsin, United States of America, 4Department of Molecular andIntegrative Physiology, University of Michigan, Ann Arbor, Michigan, United States of America, 5Department of Physiology and Biophysics, University of Washington,Seattle, Washington, United States of America, 6Department of Biomedical Informatics and Medical Education, University of Washington, Seattle, Washington, UnitedStates of America

Mechanical control of cation channels in the myogenic response

Microcirculatory vessel response to changes in pressure, known as the myogenic response, is a key component of a tissues ability to regulate blood flow. Experimental studies have not clearly elucidated the mechanical signal in the vessel wall governing steady-state reduction in vessel diameter upon an increase in intraluminal pressure. In this study, a multiscale computational model is constructed from established models of vessel wall mechanics, vascular smooth muscle (VSM) force generation, and VSM Ca(2+) handling and electrophysiology to compare the plausibility of vessel wall stress or strain as an effective mechanical signal controlling steady-state vascular contraction in the myogenic response. It is shown that, at the scale of a resistance vessel, wall stress, and not stretch (strain), is the likely physiological signal controlling the steady-state myogenic response. The model is then used to test nine candidate VSM stress-controlled channel variants by fitting two separate sets of steady-state myogenic response data. The channel variants include nonselective cation (NSC), supplementary Ca(2+) and Na(+), L-type Ca(2+), and large conductance Ca(2+)-activated K(+) channels. The nine variants are tested in turn, and model fits suggest that stress control of Ca(2+) or Na(+) influx through NSC, supplementary Ca(2+) or Na(+), or L-type Ca(2+) channels is sufficient to produce observed steady-state diameter changes with pressure. However, simulations of steady-state VSM membrane potential, cytosolic Ca(2+), and Na(+) with pressure show only that Na(+) influx through NSC channel also generates known trends with increasing pressure, indicating that stress-controlled Na(+) influx through NSC is sufficient to generate the myogenic response.

A computational analysis of the long-term regulation of arterial pressure

The asserted dominant role of the kidneys in the chronic regulation of blood pressure and in the etiology of hypertension has been debated since the 1970s. At the center of the theory is the observation that the acute relationships between arterial pressure and urine production—the acute pressure-diuresis and pressure-natriuresis curves—physiologically adapt to perturbations in pressure and/or changes in the rate of salt and volume intake. These adaptations, modulated by various interacting neurohumoral mechanisms, result in chronic relationships between water and salt excretion and pressure that are much steeper than the acute relationships. While the view that renal function is the dominant controller of arterial pressure has been supported by computer models of the cardiovascular system known as the “Guyton-Coleman model”, no unambiguous description of a computer model capturing chronic adaptation of acute renal function in blood pressure control has been presented. Here, such a model is developed with the goals of: 1. representing the relevant mechanisms in an identifiable mathematical model; 2. identifying model parameters using appropriate data; 3. validating model predictions in comparison to data; and 4. probing hypotheses regarding the long-term control of arterial pressure and the etiology of primary hypertension. The developed model reveals: long-term control of arterial blood pressure is primarily through the baroreflex arc and the renin-angiotensin system; and arterial stiffening provides a sufficient explanation for the etiology of primary hypertension associated with ageing. Furthermore, the model provides the first consistent explanation of the physiological response to chronic stimulation of the baroreflex.