The myogenic response in isolated rat cerebrovascular arteries: smooth muscle cell model

A chemo-mechanical constitutive model for muscle activation in bat wing skins

Birds, bats and insects have evolved unique wing structures to achieve a wide range of flight capabilities. Insects have relatively stiff and passive wings, birds have a complex and hierarchical feathered structure and bats have an articulated skeletal system integrated with a highly stretchable skin. The compliant skin of the wing distinguishes bats from all other flying animals and contributes to bats' remarkable, highly manoeuvrable flight performance and high energetic efficiency. The structural and functional complexity of the bat wing skin is one of the least understood although important elements of the bat flight anatomy. The wing skin has two unusual features: a discrete array of very soft elastin fibres and a discrete array of skeletal muscle fibres. The latter is intriguing because skeletal muscle is typically attached to bone, so the arrangement of intramembranous muscle in soft skin raises questions about its role in flight. In this paper, we develop a multi-scale chemo-mechanical constitutive model for bat wing skin. The chemo-mechanical model links cross-bridge cycling to a structure-based continuum model that describes the active viscoelastic behaviour of the soft anisotropic skin tissue. Continuum models at the tissue length-scale are valuable as they are easily implemented in commercial finite element codes to solve problems involving complex geometries, loading and boundary conditions. The constitutive model presented in this paper will be used in detailed finite element simulations to improve our understanding of the mechanics of bat flight in the context of wing kinematics and aerodynamic performance.

A computational model predicts sex-specific responses to calcium channel blockers in mammalian mesenteric vascular smooth muscle

The function of the smooth muscle cells lining the walls of mammalian systemic arteries and arterioles is to regulate the diameter of the vessels to control blood flow and blood pressure. Here, we describe an in silico model, which we call the 'Hernandez-Hernandez model', of electrical and Ca2+ signaling in arterial myocytes based on new experimental data indicating sex-specific differences in male and female arterial myocytes from murine resistance arteries. The model suggests the fundamental ionic mechanisms underlying membrane potential and intracellular Ca2+ signaling during the development of myogenic tone in arterial blood vessels. Although experimental data suggest that KV1.5 channel currents have similar amplitudes, kinetics, and voltage dependencies in male and female myocytes, simulations suggest that the KV1.5 current is the dominant current regulating membrane potential in male myocytes. In female cells, which have larger KV2.1 channel expression and longer time constants for activation than male myocytes, predictions from simulated female myocytes suggest that KV2.1 plays a primary role in the control of membrane potential. Over the physiological range of membrane potentials, the gating of a small number of voltage-gated K+ channels and L-type Ca2+ channels are predicted to drive sex-specific differences in intracellular Ca2+ and excitability. We also show that in an idealized computational model of a vessel, female arterial smooth muscle exhibits heightened sensitivity to commonly used Ca2+ channel blockers compared to male. In summary, we present a new model framework to investigate the potential sex-specific impact of antihypertensive drugs.

A new model for evaluating pressure-induced vascular tone in small cerebral arteries

The capacity of small cerebral arteries (SCAs) to adapt to pressure fluctuations has a fundamental physiological role and appears to be relevant in different pathological conditions. Here, we present a new computational model for quantifying the link, and its contributors, between luminal pressure and vascular tone generation in SCAs. This is assembled by combining a chemical sub-model, representing pressure-induced smooth muscle cell (SMC) signalling, with a mechanical sub-model for the tone generation and its transduction at tissue level. The devised model can accurately reproduce the impact of luminal pressure on different cytoplasmic components involved in myogenic signalling, both in the control case and when combined with some specific pharmacological interventions. Furthermore, the model is also able to capture and predict experimentally recorded pressure-outer diameter relationships obtained for vessels under control conditions, both in a Ca2 + -free bath and under drug inhibition. The modularity of the proposed framework allows the integration of new components for the study of a broad range of processes involved in the vascular function.

A computational model predicts sex-specific responses to calcium channel blockers in mammalian mesenteric vascular smooth muscle

The function of the smooth muscle cells lining the walls of mammalian systemic arteries and arterioles is to regulate the diameter of the vessels to control blood flow and blood pressure. Here, we describe an in-silico model, which we call the "Hernandez-Hernandez model", of electrical and C a 2+ signaling in arterial myocytes based on new experimental data indicating sex-specific differences in male and female arterial myocytes from murine resistance arteries. The model suggests the fundamental ionic mechanisms underlying membrane potential and intracellular C a 2+ signaling during the development of myogenic tone in arterial blood vessels. Although experimental data suggest that KV1.5 channel currents have similar amplitudes, kinetics, and voltage dependencies in male and female myocytes, simulations suggest that the KV1.5 current is the dominant current regulating membrane potential in male myocytes. In female cells, which have larger KV2.1 channel expression and longer time constants for activation than male myocytes, predictions from simulated female myocytes suggest that KV2.1 plays a primary role in the control of membrane potential. Over the physiological range of membrane potentials, the gating of a small number of voltage-gated K+ channels and L-type C a 2+ channels are predicted to drive sex-specific differences in intracellular C a 2+ and excitability. We also show that in an idealized computational model of a vessel, female arterial smooth muscle exhibits heightened sensitivity to commonly used C a 2+ channel blockers compared to male. In summary, we present a new model framework to investigate the potential sex-specific impact of anti-hypertensive drugs.

A network-based model of dynamic cerebral autoregulation

Cerebrovascular diseases continue to be one of the leading causes of morbidity and mortality in humans. Abnormalities in dynamic cerebral autoregulation (dCA) have been implicated in many of these disease conditions. Accurate models are therefore needed to better understand the complex pathophysiology behind impaired dCA. We thus present here a simple framework for modelling a vessel-driven network model of dCA in the microvasculature, as opposed to the conventional compartmental modelling approach. Network models incorporate the actual connectivity and anatomy of the vasculature, thereby allowing us to include and trace changes in the calibre and morphology of individual vessels, investigate the spatial specificity and heterogeneity of the various control mechanisms to help disentangle their contributions, and link the model parameters to the actual network physiology. The proposed control feedback mechanisms are incorporated at the level of the individual vessel, and the dynamic pressure and flow fields are solved for here within a simple vessel network. In response to an upstream pressure drop, the network is found to be able to recover cerebral blood flow (CBF) while exhibiting the characteristic autoregulatory behaviour in terms of changes in vessel calibre and the biphasic flow response. We assess the feasibility of our formulation in larger networks by comparing the simulation results to those obtained using a one-dimensional (1D) model of CBF applied to the same microvasculature network and find that our model results are in very good agreement with the 1D solution, while significantly reducing the computational cost, thus enabling more detailed models of network behaviour to be adopted in the future. Accurate and computationally feasible models of dCA that are more representative of the vasculature can help increase the translatability of haemodynamic models into the clinical environment, which would help develop more informed treatment guidelines for patients with cerebrovascular diseases.

Contractility-induced self-organization of smooth muscle cells: from multilayer cell sheets to dynamic three-dimensional clusters

Smooth muscle cells (SMCs) are mural cells that play a vital contractile function in many tissues. Abnormalities in SMC organization are associated with many diseases including atherosclerosis, asthma, and uterine fibroids. Various studies have reported that SMCs cultured on flat surfaces can spontaneously form three-dimensional clusters whose organization resembles that encountered in some of these pathological settings. Remarkably, how these structures form remains unknown. Here we combine in vitro experiments and physical modeling to show that three-dimensional clusters initiate when cellular contractile forces induce a hole in a flat SMC sheet, a process that can be modeled as the brittle fracture of a viscoelastic material. The subsequent evolution of the nascent cluster can be modeled as an active dewetting process with cluster shape evolution driven by a balance between cluster surface tension, arising from both cell contractility and adhesion, and cluster viscous dissipation. The description of the physical mechanisms governing the spontaneous emergence of these intriguing three-dimensional clusters may offer insight into SMC-related disorders.

Chemo-mechanical modeling of smooth muscle cell activation for the simulation of arterial walls under changing blood pressure

In this paper, a novel chemo-mechanical model is proposed for the description of the stretch-dependent chemical processes known as Bayliss effect and their impact on the active contraction in vascular smooth muscle. These processes are responsible for the adaptive reaction of arterial walls to changing blood pressure by which the blood vessels actively support the heart in providing sufficient blood supply for varying demands in the supplied tissues. The model is designed to describe two different stretch-dependent mechanisms observed in smooth muscle cells (SMCs): a calcium-dependent and a calcium-independent contraction. For the first one, stretch of the SMCs leads to an inlet of calcium ions which activates the myosin light chain kinase (MLCK). The increased activity of MLCK triggers the contractile units of the cells resulting in the contraction on a comparatively short time scale. For the calcium-independent contraction mechanism, stretch-dependent receptors of the cell membrane stimulate an intracellular reaction leading to an inhibition of the antagonist of MLCK, the myosin light chain phosphatase resulting in a contraction on a comparatively long time scale. An algorithmic framework for the implementation of the model in finite element programs is derived. Based thereon, it is shown that the proposed approach agrees well with experimental data. Furthermore, the individual aspects of the model are analyzed in numerical simulations of idealized arteries subject to internal pressure waves with changing intensities. The simulations show that the proposed model is able to describe the experimentally observed contraction of the artery as a reaction to increased internal pressure, which can be considered a crucial aspect of the regulatory mechanism of muscular arteries.

Modeling Reactive Hyperemia to Better Understand and Assess Microvascular Function: A Review of Techniques

Reactive hyperemia is a well-established technique for the non-invasive evaluation of the peripheral microcirculatory function, measured as the magnitude of limb re-perfusion after a brief period of ischemia. Despite widespread adoption by researchers and clinicians alike, many uncertainties remain surrounding interpretation, compounded by patient-specific confounding factors (such as blood pressure or the metabolic rate of the ischemic limb). Mathematical modeling can accelerate our understanding of the physiology underlying the reactive hyperemia response and guide in the estimation of quantities which are difficult to measure experimentally. In this work, we aim to provide a comprehensive guide for mathematical modeling techniques that can be used for describing the key phenomena involved in the reactive hyperemia response, alongside their limitations and advantages. The reported methodologies can be used for investigating specific reactive hyperemia aspects alone, or can be combined into a computational framework to be used in (pre-)clinical settings.

Modeling and experimental approaches for elucidating multi-scale uterine smooth muscle electro- and mechano-physiology: A review

The uterus provides protection and nourishment (via its blood supply) to a developing fetus, and contracts to deliver the baby at an appropriate time, thereby having a critical contribution to the life of every human. However, despite this vital role, it is an under-investigated organ, and gaps remain in our understanding of how contractions are initiated or coordinated. The uterus is a smooth muscle organ that undergoes variations in its contractile function in response to hormonal fluctuations, the extreme instance of this being during pregnancy and labor. Researchers typically use various approaches to studying this organ, such as experiments on uterine muscle cells, tissue samples, or the intact organ, or the employment of mathematical models to simulate the electrical, mechanical and ionic activity. The complexity exhibited in the coordinated contractions of the uterus remains a challenge to understand, requiring coordinated solutions from different research fields. This review investigates differences in the underlying physiology between human and common animal models utilized in experiments, and the experimental interventions and computational models used to assess uterine function. We look to a future of hybrid experimental interventions and modeling techniques that could be employed to improve the understanding of the mechanisms enabling the healthy function of the uterus.

A multiscale sliding filament model of lymphatic muscle pumping

The lymphatics maintain fluid balance by returning interstitial fluid to veins via contraction/compression of vessel segments with check valves. Disruption of lymphatic pumping can result in a condition called lymphedema with interstitial fluid accumulation. Lymphedema treatments are often ineffective, which is partially attributable to insufficient understanding of specialized lymphatic muscle lining the vessels. This muscle exhibits cardiac-like phasic contractions and smooth muscle-like tonic contractions to generate and regulate flow. To understand the relationship between this sub-cellular contractile machinery and organ-level pumping, we have developed a multiscale computational model of phasic and tonic contractions in lymphatic muscle and coupled it to a lymphangion pumping model. Our model uses the sliding filament model (Huxley in Prog Biophys Biophys Chem 7:255-318, 1957) and its adaptation for smooth muscle (Mijailovich in Biophys J 79(5):2667-2681, 2000). Multiple structural arrangements of contractile components and viscoelastic elements were trialed but only one provided physiologic results. We then coupled this model with our previous lumped parameter model of the lymphangion to relate results to experiments. We show that the model produces similar pressure, diameter, and flow tracings to experiments on rat mesenteric lymphatics. This model provides the first estimates of lymphatic muscle contraction energetics and the ability to assess the potential effects of sub-cellular level phenomena such as calcium oscillations on lymphangion outflow. The maximum efficiency value predicted (40%) is at the upper end of estimates for other muscle types. Spontaneous calcium oscillations during diastole were found to increase outflow up to approximately 50% in the range of frequencies and amplitudes tested.

A multiscale model of cerebral autoregulation

The mechanism of cerebral autoregulation ensures a continuous and sufficient blood supply to the brain to maintain normal function in the presence of changes in blood pressure. Impaired cerebral autoregulation is implicated in a range of brain diseases. We thus present here a multiscale model of cerebral autoregulation to provide a more detailed basis for a better understanding of the mechanisms behind impaired autoregulation. This model is built around a model of single arteriole, which includes a model of Nitric Oxide (NO) transport, the myogenic response, and a 4-state kinetic model coupled to a mechanical model of the vessel wall. In particular, the NO component of the model is added here to better understand the interaction mode between NO and the myogenic response, since the role of NO, the recognized effective vasodilator, is poorly understood in this context. This vessel model is then integrated within a model of the full-brain vasculature. The model is validated using a range of experimental data from the literature, both steady-state and dynamic. The model is able to predict the response of the arteriole to changes in both driving pressure and baseline pressure, indicating that the model captures well the balance between the myogenic and metabolic mechanisms. We next plan to examine the ways in which impaired autoregulation is manifested in different patient groups, potentially leading to improved therapy.

A Chemomechanobiological Model of the Long-Term Healing Response of Arterial Tissue to a Clamping Injury

Vascular clamping often causes injury to arterial tissue, leading to a cascade of cellular and extracellular events. A reliable in silico prediction of these processes following vascular injury could help us to increase our understanding thereof, and eventually optimize surgical techniques or drug delivery to minimize the amount of long-term damage. However, the complexity and interdependency of these events make translation into constitutive laws and their numerical implementation particularly challenging. We introduce a finite element simulation of arterial clamping taking into account acute endothelial denudation, damage to extracellular matrix, and smooth muscle cell loss. The model captures how this causes tissue inflammation and deviation from mechanical homeostasis, both triggering vascular remodeling. A number of cellular processes are modeled, aiming at restoring this homeostasis, i.e., smooth muscle cell phenotype switching, proliferation, migration, and the production of extracellular matrix. We calibrated these damage and remodeling laws by comparing our numerical results to in vivo experimental data of clamping and healing experiments. In these same experiments, the functional integrity of the tissue was assessed through myograph tests, which were also reproduced in the present study through a novel model for vasodilator and -constrictor dependent smooth muscle contraction. The simulation results show a good agreement with the in vivo experiments. The computational model was then also used to simulate healing beyond the duration of the experiments in order to exploit the benefits of computational model predictions. These results showed a significant sensitivity to model parameters related to smooth muscle cell phenotypes, highlighting the pressing need to further elucidate the biological processes of smooth muscle cell phenotypic switching in the future.

Computational Study of Growth and Remodeling in Ascending Thoracic Aortic Aneurysms Considering Variations of Smooth Muscle Cell Basal Tone

In this paper, we investigate the progression of Ascending Thoracic Aortic Aneurysms (ATAA) using a computational model of Growth and Remodeling (G&R) taking into account the composite (elastin, four collagen fiber families and Smooth Muscle Cells—SMCs) and multi-layered (media and adventitia) nature of the aorta. The G&R model, which is based on the homogenized Constrained Mixture theory, is implemented as a UMAT in the Abaqus finite-element package. Each component of the mixture is assigned a strain energy density function: nearly-incompressible neo-Hookean for elastin and Fung-type for collagen and SMCs. Active SMCs tension is additionally considered, through a length-tension relationship having a classic inverted parabola shape, in order to investigate its effects on the progression of ATAA in a patient-specific model. A sensitivity analysis is performed to evaluate the potential impact of variations in the parameters of the length-tension relationships. These variations reflect in variations of SMCs normal tone during ATAA progression, with active stress contributions ranging between 30% (best case scenario) and 0% (worst case scenario) of the total wall circumferential stress. Low SMCs active stress in the worst case scenarios, in fact, affect the rates of collagen deposition by which the elastin loss is gradually compensated by collagen deposition in the simulated ATAA progression, resulting eventually in larger aneurysm diameters. The types of length-tension relationships leading to a drop of SMCs active stress in our simulations reveal a critical condition which could also result in SMCs apoptosis.

The Neurovascular Unit: Focus on the Regulation of Arterial Smooth Muscle Cells

The neurovascular unit is a physiological unit present in the brain, which is constituted by elements of the nervous system (neurons and astrocytes) and the vascular system (endothelial and mural cells). This unit is responsible for the homeostasis and regulation of cerebral blood flow. There are two major types of mural cells in the brain, pericytes and smooth muscle cells. At the arterial level, smooth muscle cells are the main components that wrap around the outside of cerebral blood vessels and the major contributors to basal tone maintenance, blood pressure and blood flow distribution. They present several mechanisms by which they regulate both vasodilation and vasoconstriction of cerebral blood vessels and their regulation becomes even more important in situations of injury or pathology. In this review, we discuss the main regulatory mechanisms of brain smooth muscle cells and their contributions to the correct brain homeostasis.

A myofibre model for the study of uterine excitation-contraction dynamics

As the uterus remodels in preparation for delivery, the excitability and contractility of the uterine smooth muscle layer, the myometrium, increase drastically. But when remodelling proceeds abnormally it can contribute to preterm birth, slow progress of labour, and failure to initiate labour. Remodelling increases intercellular coupling and cellular excitability, which are the main targets of pharmaceutical treatments for uterine contraction disorders. However, the way in which electrical propagation and force development depend on intercellular coupling and cellular excitability is not fully understood. Using a computational myofibre model we study the dependency of electrical propagation and force development on intercellular coupling and cellular excitability. This model reveals that intercellular coupling determines the conduction velocity. Moreover, our model shows that intercellular coupling alone does not regulate force development. Further, cellular excitability controls whether conduction across the cells is blocked. Lastly, our model describes how cellular excitability regulates force development. Our results bridge cellular factors, targeted by drugs to regulate uterine contractions, and tissue level electromechanical properties, which are responsible for delivery. They are a step forward towards understanding uterine excitation-contraction dynamics and developing safer and more efficient pharmaceutical treatments for uterine contraction disorders.

A mathematical model for persistent post-CSD vasoconstriction

Cortical spreading depression (CSD) is the propagation of a relatively slow wave in cortical brain tissue that is linked to a number of pathological conditions such as stroke and migraine. Most of the existing literature investigates the dynamics of short term phenomena such as the depolarization and repolarization of membrane potentials or large ion shifts. Here, we focus on the clinically-relevant hour-long state of neurovascular malfunction in the wake of CSDs. This dysfunctional state involves widespread vasoconstriction and a general disruption of neurovascular coupling. We demonstrate, using a mathematical model, that dissolution of calcium that has aggregated within the mitochondria of vascular smooth muscle cells can drive an hour-long disruption. We model the rate of calcium clearance as well as the dynamical implications on overall blood flow. Based on reaction stoichiometry, we quantify a possible impact of calcium phosphate dissolution on the maintenance of F0F1-ATP synthase activity.

Overview of mathematical modeling of myocardial blood flow regulation

The oxygen consumption by the heart and its extraction from the coronary arterial blood are the highest among all organs. Any increase in oxygen demand due to a change in heart metabolic activity requires an increase in coronary blood flow. This functional requirement of adjustment of coronary blood flow is mediated by coronary flow regulation to meet the oxygen demand without any discomfort, even under strenuous exercise conditions. The goal of this article is to provide an overview of the theoretical and computational models of coronary flow regulation and to reveal insights into the functioning of a complex physiological system that affects the perfusion requirements of the myocardium. Models for three major control mechanisms of myogenic, flow, and metabolic control are presented. These explain how the flow regulation mechanisms operating over multiple spatial scales from the precapillaries to the large coronary arteries yield the myocardial perfusion characteristics of flow reserve, autoregulation, flow dispersion, and self-similarity. The review not only introduces concepts of coronary blood flow regulation but also presents state-of-the-art advances and their potential to impact the assessment of coronary microvascular dysfunction (CMD), cardiac-coronary coupling in metabolic diseases, and therapies for angina and heart failure. Experimentalists and modelers not trained in these models will have exposure through this review such that the nonintuitive and highly nonlinear behavior of coronary physiology can be understood from a different perspective. This survey highlights knowledge gaps, key challenges, future research directions, and novel paradigms in the modeling of coronary flow regulation.

Electro-Mechanical Ionic Channel Modeling for Uterine Contractions and Oxytocin Effect during Pregnancy

Uterine contractions during normal pregnancy and preterm birth are an important physiological activity. Although the cause of preterm labor is usually unknown, preterm birth creates very serious health concerns in many cases. Therefore, understanding normal birth and predicting preterm birth can help both newborn babies and their families. In our previous work, we developed a multiscale dynamic electrophysiology model of uterine contractions. In this paper, we mainly focus on the cellular level and use electromyography (EMG) and cell force generation methods to construct a new ionic channel model and a corresponding mechanical force model. Specifically, the ionic channel model takes into consideration the knowledge of individual ionic channels, which include the electrochemical and bioelectrical characteristics of individual myocytes. We develop a new sodium channel and a new potassium channel based on the experimental data from the human myometrium and the average correlations are 0.9946 and 0.9945, respectively. The model is able to generate the single spike, plateau type and bursting type of action potentials. Moreover, we incorporate the effect of oxytocin on changing the properties of the L-type and T-type calcium channels and further influencing the output action potentials. In addition, we develop a mechanical force model based on the new ionic channel model that describes the detailed ionic dynamics. Our model produces cellular mechanical force that propagates to the tissue level. We illustrate the relationship between the cellular mechanical force and the intracellular ionic dynamics and discuss the relationship between the application of oxytocin and the output mechanical force. We also propose a simplified version of the model to enable large scale simulations using sensitivity analysis method. Our results show that the model is able to reproduce the bioelectrical and electromechanical characteristics of uterine contractions during pregnancy.

Biological Relativity Requires Circular Causality but Not Symmetry of Causation: So, Where, What and When Are the Boundaries?

Since the Principle of Biological Relativity was formulated and developed there have been many implementations in a wide range of biological fields. The purpose of this article is to assess the status of the applications of the principle and to clarify some misunderstandings. The principle requires circular causality between levels of organization. But the forms of causality are also necessarily different. They contribute in asymmetric ways. Upward causation can be represented by the differential or similar equations describing the mechanics of lower level processes. Downward causation is then best represented as determining initial and boundary conditions. The questions tackled in this article are: (1) where and when do these boundaries exist? and (2) how do they convey the influences between levels? We show that not all boundary conditions arise from higher-level organization. It is important to distinguish those that do from those that don't. Both forms play functional roles in organisms, particularly in their responses to novel challenges. The forms of causation also change according to the levels concerned. These principles are illustrated with specific examples.

Biomechanical relevance of the microstructure in artery walls with a focus on passive and active components

The microstructure of arteries, consisting, in particular, of collagen, elastin, and vascular smooth muscle cells, plays a very significant role in their biomechanical response during a cardiac cycle. In this article, we highlight the microstructure and the contributions of each of its components to the overall mechanical behavior. We also describe the changes of the microstructure that occur as a result of abdominal aortic aneurysms and disease, such as atherosclerosis. We also focus on how the passive and active constituents are incorporated into a mathematical model without going into detail of the mathematical formulation. We conclude by mentioning open problems toward a better characterization of the biomechanical aspects of arteries that will be beneficial for a better understanding of cardiovascular pathophysiology.

A computational model of excitation and contraction in uterine myocytes from the pregnant rat

Aberrant uterine myometrial activities in humans are major health issues. However, the cellular and tissue mechanism(s) that maintain the uterine myometrium at rest during gestation, and that initiate and maintain long-lasting uterine contractions during delivery are incompletely understood. In this study we construct a computational model for describing the electrical activity (simple and complex action potentials), intracellular calcium dynamics and mechanical contractions of isolated uterine myocytes from the pregnant rat. The model reproduces variant types of action potentials - from spikes with a smooth plateau, to spikes with an oscillatory plateau, to bursts of spikes - that are seen during late gestation under different physiological conditions. The effects of the hormones oestradiol (via reductions in calcium and potassium selective channel conductance), oxytocin (via an increase in intracellular calcium release) and the tocolytic nifedipine (via a block of L-type calcium channels currents) on action potentials and contractions are also reproduced, which quantitatively match to experimental data. All of these results validated the cell model development. In conclusion, the developed model provides a computational platform for further investigations of the ionic mechanism underlying the genesis and control of electrical and mechanical activities in the rat uterine myocytes.

Calcium and electrical dynamics in lymphatic endothelium

KEY POINTS

Endothelial cell function in resistance arteries integrates Ca2+ signalling with hyperpolarization to promote relaxation of smooth muscle cells and increase tissue blood flow. Whether complementary signalling occurs in lymphatic endothelium is unknown. Intracellular calcium and membrane potential were evaluated in endothelial cell tubes freshly isolated from mouse collecting lymphatic vessels of the popliteal fossa. Resting membrane potential measured using intracellular microelectrodes averaged ∼-70 mV. Stimulation of lymphatic endothelium by acetylcholine or a TRPV4 channel agonist increased intracellular Ca2+ with robust depolarization. Findings from Trpv4-/- mice and with computational modelling suggest that the initial mobilization of intracellular Ca2+ leads to influx of Ca2+ and Na+ through TRPV4 channels to evoke depolarization. Lymphatic endothelial cells lack the Ca2+ -activated K+ channels present in arterial endothelium to generate endothelium-derived hyperpolarization. Absence of this signalling pathway with effective depolarization may promote rapid conduction of contraction along lymphatic muscle during lymph propulsion.

ABSTRACT

Subsequent to a rise in intracellular Ca2+ ([Ca2+ ]i ), hyperpolarization of the endothelium coordinates vascular smooth muscle relaxation along resistance arteries during blood flow control. In the lymphatic vasculature, collecting vessels generate rapid contractions coordinated along lymphangions to propel lymph, but the underlying signalling pathways are unknown. We tested the hypothesis that lymphatic endothelial cells (LECs) exhibit Ca2+ and electrical signalling properties that facilitate lymph propulsion. To study electrical and intracellular Ca2+ signalling dynamics in lymphatic endothelium, we excised collecting lymphatic vessels from the popliteal fossa of mice and removed their muscle cells to isolate intact LEC tubes (LECTs). Intracellular recording revealed a resting membrane potential of ∼-70 mV. Acetylcholine (ACh) increased [Ca2+ ]i with a time course similar to that observed in endothelium of resistance arteries (i.e. rapid initial peak with a sustained 'plateau'). In striking contrast to the endothelium-derived hyperpolarization (EDH) characteristic of arteries, LECs depolarized (>15 mV) to either ACh or TRPV4 channel activation. This depolarization was facilitated by the absence of Ca2+ -activated K+ (KCa ) channels as confirmed with PCR, persisted in the absence of extracellular Ca2+ , was abolished by LaCl3 and was attenuated ∼70% in LECTs from Trpv4-/- mice. Computational modelling of ion fluxes in LECs indicated that omitting K+ channels supports our experimental results. These findings reveal novel signalling events in LECs, which are devoid of the KCa activity abundant in arterial endothelium. Absence of EDH with effective depolarization of LECs may promote the rapid conduction of contraction waves along lymphatic muscle during lymph propulsion.

An effective model of cerebrovascular pressure reactivity and blood flow autoregulation

Understanding cerebral blood flow dynamics is crucial for the care of patients at risk of poor cerebral perfusion. We describe an effective model of cerebral hemodynamics designed to reveal important macroscopic features of cerebral blood flow without having to resolve the detailed microvasculature of the brain. Based on principles of fluid and elastic dynamics and vascular pressure-reactivity, the model quantifies the physical means by which the vasculature executes autoregulatory reflexes. We demonstrate that the frequency response of the proposed model matches experimental measurements and explains the influence of mechanical factors on the autoregulatory performance. Analysis of the model indicates the existence of an optimal mean arterial pressure which minimizes the sensitivity of the flow to changes in perfusion pressure across the frequency spectrum of physiological oscillations. We highlight the simplicity of the model and its potential to improve monitoring of brain perfusion via real-time computational simulations of cerebro- and cardio-vascular interventions.

Importance of Altered Levels of SERCA, IP3R, and RyR in Vascular Smooth Muscle Cell

Calcium cycling between the sarcoplasmic reticulum (SR) and the cytosol via the sarco-/endoplasmic reticulum Ca-ATPase (SERCA) pump, inositol-1,4,5-triphosphate receptor (IP₃R), and Ryanodine receptor (RyR), plays a major role in agonist-induced intracellular calcium ([Ca²⁺]_cyt) dynamics in vascular smooth muscle cells (VSMC). Levels of these calcium handling proteins in SR get altered under disease conditions. We have developed a mathematical model to understand the significance of altered levels of SERCA, IP₃R, and RyR on the intracellular calcium dynamics of VSMC and to understand how variation in protein levels that arise due to diabetes contribute to different VSMC behavior and thus vascular disease. SR is modeled as a single continuous entity with homogeneous intra-SR calcium. Model results show that agonist-induced intracellular calcium dynamics can be modified by changing the levels of SERCA, IP₃R, and/or RyR. Lowering SERCA level will enable intracellular calcium oscillations at low agonist concentrations whereas lowered levels of IP₃R and RyR need higher agonist concentration for intracellular calcium oscillations. This research suggests that reduced SERCA level is the main factor responsible for the reduced intracellular calcium transients and contractility in VSMCs.

Identification of the mechanical parameters for the human uterus in vivo using intrauterine pressure measurements

There are limited experimental data to characterize the mechanical response of human myometrium. A method is presented in this work to identify mechanical parameters describing the active response of human myometrium from the in vivo intrauterine pressure measurements. The human uterine contraction during labor is simulated by implementing a coupled model in a finite element scheme, and the intrauterine pressure is evaluated as the outcome. A meta-model is developed to approximate the finite element simulation response with a lower computational cost and used to identify model parameters through fitting its prediction to the in vivo measurements. Copyright © 2016 John Wiley & Sons, Ltd.

A combination of experimental and numerical methods to investigate the role of strain rate on the mechanical properties and collagen fiber orientations of the healthy and atherosclerotic human coronary arteries

Atherosclerosis enables to alter not only the microstructural but also the physical properties of the arterial walls by plaque forming. Few studies so far have been conducted to calculate the isotropic or anisotropic mechanical properties of the healthy and atherosclerotic human coronary arteries. To date there is a paucity of knowledge on the mechanical response of the arteries under different strain rates. Therefore, the objective of the concurrent research was to comprehend whether the alteration in the strain rates of the human atherosclerotic arteries in comparison with the healthy ones contribute to the biomechanical behaviors. To do this, healthy and atherosclerotic human coronary arteries were removed from 18 individuals during autopsy. Histological analyses by both an expert histopathologist and an imaged-based recognizer software were performed to figure out the average angle of collagen fibers in the healthy and atherosclerotic arterial walls. Thereafter, the samples were subjected to 3 diverse strain rates, i.e., 5, 20, and 50 mm/min, until the material failure occurs. The stress-strain diagrams of the arterial tissues were calculated in order to capture their linear elastic and nonlinear hyperelastic mechanical properties. In addition, Artificial Neural Networks (ANNs) was employed to predict the alteration of mean angle of collagen fibers during load bearing up to failure. The findings suggest that strain rate has a significant (p < 0.05) role in the linear elastic and nonlinear hyperelastic mechanical properties as well as the mean angle of collagen fibers of the atherosclerotic arteries, whereas no specific impact on the healthy ones. Furthermore, the mean angle of collagen fibers during the load bearing up to the failure at each strain rate was well predicted by the proposed ANNs code.

On the roles of vascular smooth muscle contraction in cerebral blood flow autoregulation - a modeling perspective

We here review existing models of vascular smooth muscle cell, endothelial cell and cell-cell communication, which have been developed to better understand vascular tone and blood flow autoregulation. In particular, we discuss models that intended to explain modulation of myogenic tone by intraluminal pressure in resistance arterioles. Modeling efforts in the recent past have witnessed a shift from empirical models to models with mechanistic details that underscore different physical aspects of vascular regulation. Future models should synthesize mechanistic interactions in a hierarchy, from molecular regulation of ion channels to whole organ blood flow control.

Predicted effects of nitric oxide and superoxide on the vasoactivity of the afferent arteriole

We expanded a published mathematical model of an afferent arteriole smooth muscle cell in rat kidney (Edwards A, Layton, AT. Am J Physiol Renal Physiol 306: F34-F48, 2014) to understand how nitric oxide (NO) and superoxide (O(2)(-)) modulate the arteriolar diameter and its myogenic response. The present model includes the kinetics of NO and O(2)(-) formation, diffusion, and reaction. Also included are the effects of NO and its second messenger cGMP on cellular Ca²⁺ uptake and efflux, Ca²⁺-activated K⁺ currents, and myosin light chain phosphatase activity. The model considers as well pressure-induced increases in O(2)(-) production, O(2)(-)-mediated regulation of L-type Ca²⁺ channel conductance, and increased O(2)(-) production in spontaneous hypertensive rats (SHR). Our results indicate that elevated O(2)(-) production in SHR is sufficient to account for observed differences between normotensive and hypertensive rats in the response of the afferent arteriole to NO synthase inhibition, Tempol, and angiotensin II at baseline perfusion pressures. In vitro, whether the myogenic response is stronger in SHR remains uncertain. Our model predicts that if mechanosensitive cation channels are not modulated by O(2)(-), then fractional changes in diameter induced by pressure elevations should be smaller in SHR than in normotensive rats. Our results also suggest that most NO diffuses out of the smooth muscle cell without being consumed, whereas most O(2)(-) is scavenged, by NO and superoxide dismutase. Moreover, the predicted effects of superoxide on arteriolar constriction are not predominantly due to its scavenging of NO.

Simulating uterine contraction by using an electro-chemo-mechanical model

Contractions of uterine smooth muscle cells consist of a chain of physiological processes. These contractions provide the required force to expel the fetus from the uterus. The inclusion of these physiological processes is, therefore, imperative when studying uterine contractions. In this study, an electro-chemo-mechanical model to replicate the excitation, activation, and contraction of uterine smooth muscle cells is developed. The presented modeling strategy enables efficient integration of knowledge about physiological processes at the cellular level to the organ level. The model is implemented in a three-dimensional finite element setting to simulate uterus contraction during labor in response to electrical discharges generated by pacemaker cells and propagated within the myometrium via gap junctions. Important clinical factors, such as uterine electrical activity and intrauterine pressure, are predicted using this simulation. The predictions are in agreement with clinically measured data reported in the literature. A parameter study is also carried out to investigate the impact of physiologically related parameters on the uterine contractility.

Conducted vasoreactivity: the dynamical point of view

Conducted vasodilation is part of the physiological response to increasing metabolic demand of the tissue. Similar responses can be elicited by focal electrical or chemical stimulation. Some evidence suggests an endothelial pathway for nondecremental transmission of hyperpolarizing pulses. However, the underlying mechanisms are debated. Here, we focus on dynamical aspects of the problem hypothesizing the existence of a bistability-powered mechanism for regenerative pulse transmission along the endothelium. Bistability implies that the cell can have two different stable resting potentials and can switch between those states following an appropriate stimulus. Bistability is possible if the current–voltage curve is N shaped instead of monotonically increasing. Specifically, the presence of an inwardly rectifying potassium current may provide the endothelial cell with such properties. We provide a theoretical analysis as well as numerical simulations of both single- and multiunit bistable systems mimicking endothelial cells to investigate the self-consistence and stability of the proposed mechanism. We find that the individual cell may switch readily between two stable potentials. An array of coupled cells, however, as found in the vascular wall, requires a certain adaptation of the membrane currents after a switch, in order to switch back. Although the formulation is generic, we suggest a combination of specific membrane currents that could underlie the phenomenon.

Bio-Chemo-Mechanical Models of Vascular Mechanics

Models of vascular mechanics are necessary to predict the response of an artery under a variety of loads, for complex geometries, and in pathological adaptation. Classic constitutive models for arteries are phenomenological and the fitted parameters are not associated with physical components of the wall. Recently, microstructurally-linked models have been developed that associate structural information about the wall components with tissue-level mechanics. Microstructurally-linked models are useful for correlating changes in specific components with pathological outcomes, so that targeted treatments may be developed to prevent or reverse the physical changes. However, most treatments, and many causes, of vascular disease have chemical components. Chemical signaling within cells, between cells, and between cells and matrix constituents affects the biology and mechanics of the arterial wall in the short- and long-term. Hence, bio-chemo-mechanical models that include chemical signaling are critical for robust models of vascular mechanics. This review summarizes bio-mechanical and bio-chemo-mechanical models with a focus on large elastic arteries. We provide applications of these models and challenges for future work.

Investigating the role of smooth muscle cells in large elastic arteries: a finite element analysis

Physiological loading in large elastic arteries is considered to be mainly carried by the passive components of the media but it is not known how much the contraction of the smooth muscle cells is actually involved in the load carrying. Smooth muscle contraction is considered to occur in a relatively slow time domain but the contraction is able to produce significant tension. In the present work the role of smooth muscle contraction in large elastic arteries is investigated by analyzing how changes in the intracellular calcium, and thereby the active tone of smooth muscle cells, influence the deformation and stress behavior; different intracellular calcium functions and medial wall thicknesses with cycling internal pressure are studied. In particular, a recently proposed mechanochemical model (Murtada et al., 2012. J. Theor. Biol. 297, 176-186), which links intracellular calcium with mechanical contraction and an anisotropic model representing the elastin/collagen composite, was implemented into a 3D finite element framework. Details of the implementation procedure are described and a verification of the model implementation is provided by means of the isometric contraction/relaxation analysis of a medial strip at optimal muscle length. In addition, numerically obtained pressure-radius relationships of arterial rings modeled with one and two layers are analyzed with different geometries and at different calcium levels; a comparison with the Laplace equation is provided. Finally, a two-layer arterial ring is loaded with a realistic pressure wave and with various intracellular calcium functions (different amplitudes and mean values) and medial wall thicknesses; residual stresses are considered. The finite element results show that changes in the calcium amplitudes hardly have an influence on the current inner ring radius and the circumferential stress. However, an increase in the mean intracellular calcium value and the medial wall thickness leads to a clear influence on the deformation and the stress behavior.

Changes in biomechanical properties of the coronary artery wall contribute to maintained contractile responses to endothelin-1 in atherosclerosis

AIMS

Our aim was to determine whether alterations in biomechanical properties of human diseased compared to normal coronary artery contribute to changes in artery responsiveness to endothelin-1 in atherosclerosis.

MAIN METHODS

Concentration-response curves were constructed to endothelin-1 in normal and diseased coronary artery. The passive mechanical properties of arteries were determined using tensile ring tests from which finite element models of passive mechanical properties of both groups were created. Finite element modelling of artery endothelin-1 responses was then performed.

KEY FINDINGS

Maximum responses to endothelin-1 were significantly attenuated in diseased (27±3 mN, n=55) compared to normal (38±2 mN, n=68) artery, although this remained over 70% of control. There was no difference in potency (pD2 control=8.03±0.06; pD2 diseased=7.98±0.06). Finite element modelling of tensile ring tests resulted in hyperelastic shear modulus μ=2004±410 Pa and hardening exponent α=22.8±2.2 for normal wall and μ=2464±1075 Pa and α=38.3±6.7 for plaque tissue and distensibility of diseased vessels was decreased. Finite element modelling of active properties of both groups resulted in higher muscle contractile strain (represented by thermal reactivity) of the atherosclerotic artery model than the normal artery model. The models suggest that a change in muscle response to endothelin-1 occurs in atherosclerotic artery to increase its distensibility towards that seen in normal artery.

SIGNIFICANCE

Our data suggest that an adaptation occurs in medial smooth muscle of atherosclerotic coronary artery to maintain distensibility of the vessel wall in the presence of endothelin-1. This may contribute to the vasospastic effect of locally increased endothelin-1 production that is reported in this condition.

A continuum model for excitation-contraction of smooth muscle under finite deformations

The main focus in most of the continuum based muscle models is the mechanics of muscle contraction while other physiological processes governing muscle contraction, e.g., cell membrane excitation and activation, are ignored. These latter processes are essential to initiate contraction and to determine the amount of generated force, and by excluding them, the developed model cannot replicate the true behavior of the muscle in question. The aim of this study is to establish a thermodynamically and physiologically consistent framework which allows us to model smooth muscle contraction by including cell membrane excitability and kinetics of myosin phosphorylation, along with dynamics of smooth muscle contraction. The model accounts for these processes through a set of coupled dissipative constitutive equations derived by applying first principles. To show the performance of the derived model, it is evaluated for two different cases: a chemo-mechanical study of pig taenia coli cells where the excitation process is excluded, and an electro-chemo-mechanical study of rat myometrium. The results show that the model is able to replicate important aspects of the smooth muscle excitation-contraction process.

Calcium dynamics underlying the myogenic response of the renal afferent arteriole

The renal afferent arteriole reacts to an elevation in blood pressure with an increase in muscle tone and a decrease in luminal diameter. This effect, known as the myogenic response, is believed to stabilize glomerular filtration and to protect the glomerulus from systolic blood pressure increases, especially in hypertension. To study the mechanisms underlying the myogenic response, we developed a mathematical model of intracellular Ca(2+) signaling in an afferent arteriole smooth muscle cell. The model represents detailed transmembrane ionic transport, intracellular Ca(2+) dynamics, the kinetics of myosin light chain phosphorylation, and the mechanical behavior of the cell. It assumes that the myogenic response is initiated by pressure-induced changes in the activity of nonselective cation channels. Our model predicts spontaneous vasomotion at physiological luminal pressures and KCl- and diltiazem-induced diameter changes comparable to experimental findings. The time-periodic oscillations stem from the dynamic exchange of Ca(2+) between the cytosol and the sarcoplasmic reticulum, coupled to the stimulation of Ca(2+)-activated potassium (KCa) and chloride (ClCa) channels, and the modulation of voltage-activated L-type channels; blocking sarco/endoplasmic reticulum Ca(2+) pumps, ryanodine receptors (RyR), KCa, ClCa, or L-type channels abolishes these oscillations. Our results indicate that the profile of the myogenic response is also strongly dependent on the conductance of ClCa and L-type channels, as well as the activity of plasmalemmal Ca(2+) pumps. Furthermore, inhibition of KCa is not necessary to induce myogenic contraction. Lastly, our model suggests that the kinetic behavior of L-type channels results in myogenic kinetics that are substantially faster during constriction than during dilation, consistent with in vitro observations (Loutzenhiser R, Bidani A, Chilton L. Circ. Res. 90: 1316-1324, 2002).

Microstructural constitutive model of active coronary media

Although vascular smooth muscle cells (VSMCs) are pivotal in physiology and pathology, there is a lack of detailed morphological data on these cells. The objective of this study was to determine dimensions (width and length) and orientation of swine coronary VSMCs and to develop a microstructural constitutive model of active media. The dimensions, spatial aspect ratio and orientation angle of VSMCs measured at zero-stress state were found to follow continuous normal (or bimodal normal) distributions. The VSMCs aligned off circumferential direction of blood vessels with symmetrical polar angles 18.7° ± 10.9°, and the local VSMC deformation was affine with tissue-level deformation. A microstructure-based active constitutive model was developed to predict the biaxial vasoactivity of coronary media, based on experimental measurements of geometrical and deformation features of VSMCs. The results revealed that the axial active response of blood vessels is associated with multi-axial contraction as well as oblique VSMC arrangement. The present morphological database is essential for developing accurate structural models and is seminal for understanding the biomechanics of muscular vessels.

A three-dimensional chemo-mechanical continuum model for smooth muscle contraction

Based on two fields, namely the placement and the calcium concentration, a chemo-mechanically coupled three-dimensional model, describing the contractile behaviour of smooth muscles, is presented by means of a strain energy function. The strain energy function (Schmitz and Böl, 2011) is additively decomposed into a passive part, relating to elastin and collagen, and an active calcium-driven part related to the chemical contraction of the smooth muscle cells. For the description of the calcium phase the four state cross-bridge model of Hai and Murphy (Hai and Murphy, 1988) has been implemented into the finite element method. Beside three-dimensional illustrative boundary-value problems demonstrating the features of the presented modelling concept, simulations on an idealised artery document the applicability of the model to more realistic geometries.

Extracellular matrix and the mechanics of large artery development

The large, elastic arteries, as their name suggests, provide elastic distention and recoil during the cardiac cycle in vertebrate animals. The arteries are distended from the pressure of ejecting blood during the active contraction of the left ventricle (LV) during systole and recoil to their original dimensions during relaxation of the LV during diastole. The cyclic distension occurs with minimal energy loss, due to the elastic properties of one of the major structural extracellular matrix (ECM) components, elastin. The maximum distension is limited to prevent damage to the artery by another major ECM component, collagen. The mix of ECM components in the wall largely determines the passive mechanical behavior of the arteries and the subsequent load on the heart during systole. While much research has focused on initial artery formation, there has been less attention on the continuing development of the artery to produce the mature composite wall complete with endothelial cells (ECs), smooth muscle cells (SMCs), and the necessary mix of ECM components for proper cardiovascular function. This review focuses on the physiology of large artery development, including SMC differentiation and ECM production. The effects of hemodynamic forces and ECM deposition on the evolving arterial structure and function are discussed. Human diseases and mouse models with genetic mutations in ECM proteins that affect large artery development are summarized. A review of constitutive models and growth and remodeling theories is presented, along with future directions to improve understanding of ECM and the mechanics of large artery development.

Modeling Ca2+ signaling in the microcirculation: intercellular communication and vasoreactivity

A network of intracellular signaling pathways and complex intercellular interactions regulate calcium mobilization in vascular cells, arteriolar tone, and blood flow. Different endothelium-derived vasoreactive factors have been identified and the importance of myoendothelial communication in vasoreactivity is now well appreciated. The ability of many vascular networks to conduct signals upstream also is established. This phenomenon is critical for both short-term changes in blood perfusion as well as long-term adaptations of a vascular network. In addition, in a phenomenon termed vasomotion, arterioles often exhibit spontaneous oscillations in diameter. This is thought to improve tissue oxygenation and enhance blood flow. Experimentation has begun to reveal important aspects of the regulatory machinery and the significance of these phenomena for the regulation of local perfusion and oxygenation. Mathematical modeling can assist in elucidating the complex signaling mechanisms that participate in these phenomena. This review highlights some of the important experimental studies and relevant mathematical models that provide the current understanding of these mechanisms in vasoreactivity.

Cross-bridge apparent rate constants of human gallbladder smooth muscle

This paper studies human gallbladder (GB) smooth muscle contractions. A two-state cross-bridge model was used to estimate the apparent attachment and detachment rate constants, as well as increased Ca2+ concentration from the peak active stress during the isometric contraction. The active stress was computed from a mechanical model based entirely on non-invasive routine ultrasound scans. In the two-state cross-bridge model, the two apparent rate constants, representing the total attached/detached cross-bridges, respectively, were estimated using active stress prediction for 51 subjects undergoing cholecystokinin-provocation test, together with estimates from the four-state cross-bridge model for a swine carotid, bovine tracheal and guinea pig GB smooth muscles. The study suggests that the apparent rate constants should be patient-specific, i.e. patients with a lower stress level are characterized by smaller apparent rate constants. In other words, the diseased GB may need to develop fast cycling cross-bridges to compensate in the emptying process. This is a first step towards more quantitative and non-invasive measures of GB pain, and may provide useful insight in understanding GB motility and developing effective drug treatments.

Time course of isotonic shortening and the underlying contraction mechanism in airway smooth muscle

Although the structure of the contractile unit in smooth muscle is poorly understood, some of the mechanical properties of the muscle suggest that a sliding-filament mechanism, similar to that in striated muscle, is also operative in smooth muscle. To test the applicability of this mechanism to smooth muscle function, we have constructed a mathematical model based on a hypothetical structure of the smooth muscle contractile unit: a side-polar myosin filament sandwiched by actin filaments, each attached to the equivalent of a Z disk. Model prediction of isotonic shortening as a function of time was compared with data from experiments using ovine tracheal smooth muscle. After equilibration and establishment of in situ length, the muscle was stimulated with ACh (100 μM) until force reached a plateau. The muscle was then allowed to shorten isotonically against various loads. From the experimental records, length-force and force-velocity relationships were obtained. Integration of the hyperbolic force-velocity relationship and the linear length-force relationship yielded an exponential function that approximated the time course of isotonic shortening generated by the modeled sliding-filament mechanism. However, to obtain an accurate fit, it was necessary to incorporate a viscoelastic element in series with the sliding-filament mechanism. The results suggest that a large portion of the shortening is due to filament sliding associated with muscle activation and that a small portion is due to continued deformation associated with an element that shows viscoelastic or power-law creep after a step change in force.

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 tissue's 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.

Calcium dynamics and signaling in vascular regulation: computational models

Calcium is a universal signaling molecule with a central role in a number of vascular functions including in the regulation of tone and blood flow. Experimentation has provided insights into signaling pathways that lead to or affected by Ca(2+) mobilization in the vasculature. Mathematical modeling offers a systematic approach to the analysis of these mechanisms and can serve as a tool for data interpretation and for guiding new experimental studies. Comprehensive models of calcium dynamics are well advanced for some systems such as the heart. This review summarizes the progress that has been made in modeling Ca(2+) dynamics and signaling in vascular cells. Model simulations show how Ca(2+) signaling emerges as a result of complex, nonlinear interactions that cannot be properly analyzed using only a reductionist's approach. A strategy of integrative modeling in the vasculature is outlined that will allow linking macroscale pathophysiological responses to the underlying cellular mechanisms.

Theoretical models for coronary vascular biomechanics: progress & challenges

A key aim of the cardiac Physiome Project is to develop theoretical models to simulate the functional behaviour of the heart under physiological and pathophysiological conditions. Heart function is critically dependent on the delivery of an adequate blood supply to the myocardium via the coronary vasculature. Key to this critical function of the coronary vasculature is system dynamics that emerge via the interactions of the numerous constituent components at a range of spatial and temporal scales. Here, we focus on several components for which theoretical approaches can be applied, including vascular structure and mechanics, blood flow and mass transport, flow regulation, angiogenesis and vascular remodelling, and vascular cellular mechanics. For each component, we summarise the current state of the art in model development, and discuss areas requiring further research. We highlight the major challenges associated with integrating the component models to develop a computational tool that can ultimately be used to simulate the responses of the coronary vascular system to changing demands and to diseases and therapies.

Influence of dispersion in myosin filament orientation and anisotropic filament contractions in smooth muscle

A new constitutive model for the biomechanical behaviour of smooth muscle tissue is proposed. The active muscle contraction is accomplished by the relative sliding between actin and myosin filaments, comprising contractile units in the smooth muscle cells. The orientation of the myosin filaments, and thereby the contractile units, are taken to exhibit a statistical dispersion around a preferred direction. The number of activated cross-bridges between the actin and myosin filaments governs the contractile force generated by the muscle and also the contraction speed. A strain-energy function is used to describe the mechanical behaviour of the smooth muscle tissue. Besides the active contractile apparatus, the mechanical model also incorporates a passive elastic part. The constitutive model was compared to histological and isometric tensile test results for smooth muscle tissue from swine carotid artery. In order to be able to predict the active stress at different muscle lengths, a filament dispersion significantly larger than the one observed experimentally was required. Furthermore, a comparison of the predicted active stress for a case of uniaxially oriented myosin filaments and a case of filaments with a dispersion based on the experimental histological data shows that the difference in generated stress is noticeable but limited. Thus, the results suggest that myosin filament dispersion alone cannot explain the increase in active muscle stress with increasing muscle stretch.

Cellular mechanisms underlying nitric oxide-induced vasodilation of descending vasa recta

It has been observed that vasoactivity of explanted descending vasa recta (DVR) is modulated by intrinsic nitric oxide (NO) and superoxide (O(2)(-)) production (Cao C, Edwards A, Sendeski M, Lee-Kwon W, Cui L, Cai CY, Patzak A, Pallone TL. Am J Physiol Renal Physiol 299: F1056-F1064, 2010). To elucidate the cellular mechanisms by which NO, O(2)(-) and hydrogen peroxide (H(2)O(2)) modulate DVR pericyte cytosolic Ca(2+) concentration ([Ca](cyt)) and vasoactivity, we expanded our mathematical model of Ca(2+) signaling in pericytes. We incorporated simulations of the pathways that translate an increase in [Ca](cyt) to the activation of myosin light chain (MLC) kinase and cell contraction, as well as the kinetics of NO and reactive oxygen species formation and their effects on [Ca](cyt) and MLC phosphorylation. The model reproduced experimentally observed trends of DVR vasoactivity that accompany exposure to N(ω)-nitro-L-arginine methyl ester, 8-Br-cGMP, Tempol, and H(2)O(2). Our results suggest that under resting conditions, NO-induced activation of cGMP maintains low levels of [Ca](cyt) and MLC phosphorylation to minimize basal tone. This results from stimulation of Ca(2+) uptake from the cytosol into the SR via SERCA pumps, Ca(2+) efflux into the extracellular space via plasma membrane Ca(2+) pumps, and MLC phosphatase (MLCP) activity. We predict that basal concentrations of O(2)(-) and H(2)O(2) have negligible effects on Ca(2+) signaling and MLC phosphorylation. At concentrations above 1 nM, O(2)(-) is predicted to modulate [Ca(cyt)] and MCLP activity mostly by reducing NO bioavailability. The DVR vasoconstriction that is induced by high concentrations of H(2)O(2) can be explained by H(2)O(2)-mediated downregulation of MLCP and SERCA activity. We conclude that intrinsic generation of NO by the DVR wall may be sufficient to inhibit vasoconstriction by maintaining suppression of MLC phosphorylation.