Modeling cerebrovascular responses to assess the impact of the collateral circulation following middle cerebral artery occlusion.

OBJECTIVE: An improved understanding of the role of the leptomeningeal collateral circulation in blood flow compensation following middle cerebral artery (MCA) occlusion can contribute to more effective treatment development for ischemic stroke. The present study introduces a model of the cerebral circulation to predict cerebral blood flow and tissue oxygenation following MCA occlusion. METHODS: The model incorporates flow regulation mechanisms based on changes in pressure, shear stress, and metabolic demand. Oxygen saturation in cerebral vessels and tissue is calculated using a Krogh cylinder model. The model is used to assess the effects of changes in oxygen demand and arterial pressure on cerebral blood flow and oxygenation after MCA occlusion. RESULTS: An increase from five to 11 leptomeningeal collateral vessels was shown to increase the oxygen saturation in the region distal to the occlusion by nearly 100%. Post-occlusion, the model also predicted a loss of autoregulation and a decrease in flow to the ischemic territory as oxygen demand was increased; these results were consistent with data from experiments that induced cerebral ischemia. CONCLUSIONS: This study highlights the importance of leptomeningeal collaterals following MCA occlusion and reinforces the idea that lower oxygen demand and higher arterial pressure improve conditions of flow and oxygenation.

Modeling acute and chronic vascular responses to a major arterial occlusion.

OBJECTIVE: To incorporate chronic vascular adaptations into a mathematical model of the rat hindlimb to simulate flow restoration following total occlusion of the femoral artery. METHODS: A vascular wall mechanics model is used to simulate acute and chronic vascular adaptations in the collateral arteries and collateral-dependent arterioles of the rat hindlimb. On an acute timeframe, the vascular tone of collateral arteries and distal arterioles is determined by responses to pressure, shear stress, and metabolic demand. On a chronic timeframe, sustained dilation of arteries and arterioles induces outward vessel remodeling represented by increased passive vessel diameter (arteriogenesis), and low venous oxygen saturation levels induce the growth of new capillaries represented by increased capillary number (angiogenesis). RESULTS: The model predicts that flow compensation to an occlusion is enhanced primarily by arteriogenesis of the collateral arteries on a chronic time frame. Blood flow autoregulation is predicted to be disrupted and to occur for higher pressure values following femoral arterial occlusion. CONCLUSIONS: Structural adaptation of the vasculature allows for increased blood flow to the collateral-dependent region after occlusion. Although flow is still below pre-occlusion levels, model predictions indicate that interventions which enhance collateral arteriogenesis would have the greatest potential for restoring flow.

Modeling acute blood flow responses to a major arterial occlusion.

OBJECTIVE: The development of earlier and less invasive treatments for peripheral arterial disease requires a more complete understanding of vascular responses following a major arterial occlusion. A mechanistic model of the vasculature of the rat hindlimb is developed to predict acute (immediate) changes in vessel diameters and smooth muscle tone following femoral arterial occlusion. METHODS: Vascular responses of collateral arteries and distal arterioles to changes in pressure, shear stress, and metabolism are assessed before and after occlusion. The effects of exercise are also simulated and compared with venous flow measurements from WKY rats. RESULTS: The model identifies collateral arteries as the primary contributors to flow compensation following occlusion. Increasing the number of capillaries has minimal effect on blood flow while increasing the number of collateral arteries significantly increases flow, since the primary site of resistance shifts upstream to the collateral arteries following occlusion. Despite significant collateral dilation, calf flow remains below pre-occlusion levels and the deficit becomes more severe with increased activity. CONCLUSIONS: Although unable to compensate fully for an occlusion, the model demonstrates the importance of the shear response in collateral arteries and the metabolic response in the distal microcirculation in acute adaptations to a major arterial occlusion.

Two-dimensional Finite Element Model of Breast Cancer Cell Motion Through a Microfluidic Channel.

A two-dimensional model for red blood cell motion is adapted to consider the dynamics of breast cancer cells in a microfluidic channel. Adjusting parameters to make the membrane stiffer, as is the case with breast cancer cells compared with red blood cells, allows the model to produce reasonable estimates of breast cancer cell trajectories through the channel. In addition, the model produces estimates of quantities not as easily obtained from experiment such as velocity and stress field information throughout the fluid and on the cell membrane. This includes locations of maximum stress along the membrane wall. A sensitivity analysis shows that the model is capable of producing useful insights into various systems involving breast cancer cells. Current results suggest that dynamics taking place when cells are near other objects are most sensitive to membrane and cytoplasm elasticity, dynamics taking place when cells are not near other objects are most sensitive to cytoplasm viscosity, and dynamics are significantly affected by low membrane bending elasticity. These results suggest that continued calibration and application of this model can yield useful predictions in other similar systems.

The effect of the endothelial surface layer on cell-cell interactions in microvessel bifurcations.

Red blood cells (RBCs) carry oxygen and make up 40-45% of blood by volume in large vessels down to 10% or less in smaller capillaries. Because of their finite size and large volume fraction, they are heterogeneously distributed throughout the body. This is partially because RBCs are distributed or partitioned nonuniformly at diverging vessel bifurcations where blood flows from one vessel into two. Despite its increased recognition as an important player in the microvasculature, few studies have explored how the endothelial surface layer (ESL; a vessel wall coating) may affect partitioning and RBC dynamics at diverging vessel bifurcations. Here, we use a mathematical and computational model to consider how altering ESL properties, as can occur in pathological scenarios, change RBC partitioning, deformation, and penetration of the ESL. The two-dimensional finite element model considers pairs of cells, represented by interconnected viscoelastic elements, passing through an ESL-lined diverging vessel bifurcation. The properties of the ESL include the hydraulic resistivity and an osmotic pressure difference modeling how easily fluid flows through the ESL and how easily the ESL is structurally compressed, respectively. We find that cell-cell interaction leads to more uniform partitioning and greatly enhances the effects of ESL properties, especially for deformation and penetration. This includes the trend that increased hydraulic resistivity leads to more uniform partitioning, increased deformation, and decreased penetration. It also includes the trend that decreased osmotic pressure increases penetration.

A three-dimensional mathematical and computational model of necrotizing enterocolitis.

Necrotizing enterocolitis (NEC) is a severe disease that affects the gastrointestinal (GI) tract of premature infants. Different areas of NEC research have often been isolated from one another and progress on the role of the inflammatory response in NEC, on the dynamics of epithelial layer healing, and on the positive effects of breast feeding have not been synthesized to produce a more integrated understanding of the pathogenesis of NEC. We seek to synthesize these areas of research by creating a mathematical model that incorporates the current knowledge on these aspects. Unlike previous models that are based on ordinary differential equations, our mathematical model takes into account not only transient effects but also spatial effects. A system of nonlinear transient partial differential equations is solved numerically using cell-centered finite differences and an explicit Euler method. The model is used to track the evolution of a prescribed initial injured area in the intestinal wall. It is able to produce pathophysiologically realistic results; decreasing the initial severity of the injury in the system and introducing breast feeding to the system both lead to healthier overall simulations, and only a small fraction of epithelial injuries lead to full-blown NEC. In addition, in the model, changing the initial shape of the injured area can significantly alter the overall outcome of a simulation. This finding suggests that taking into account spatial effects may be important in assessing the outcome for a given NEC patient. This model can provide a platform with which to test competing hypotheses regarding pathological mechanisms of inflammation in NEC, suggest experimental approaches by which to clarify pathogenic drivers of NEC, and may be used to derive potential intervention strategies.

Dependence of red blood cell dynamics in microvessel bifurcations on the endothelial surface layer's resistance to flow and compression.

Red blood cells (RBCs) make up 40-45% of blood and play an important role in oxygen transport. That transport depends on the RBC distribution throughout the body, which is highly heterogeneous. That distribution, in turn, depends on how RBCs are distributed or partitioned at diverging vessel bifurcations where blood flows from one vessel into two. Several studies have used mathematical modeling to consider RBC partitioning at such bifurcations in order to produce useful insights. These studies, however, assume that the vessel wall is a flat impenetrable homogeneous surface. While this is a good first approximation, especially for larger vessels, the vessel wall is typically coated by a flexible, porous endothelial glycocalyx or endothelial surface layer (ESL) that is on the order of 0.5-1 µm thick. To better understand the possible effects of this layer on RBC partitioning, a diverging capillary bifurcation is analyzed using a flexible, two-dimensional model. In addition, the model is also used to investigate RBC deformation and RBC penetration of the ESL region when ESL properties are varied. The RBC is represented using interconnected viscoelastic elements. Stokes flow equations (viscous flow) model the surrounding fluid. The flow in the ESL is modeled using the Brinkman approximation for porous media with a corresponding hydraulic resistivity. The ESL's resistance to compression is modeled using an osmotic pressure difference. One cell passes through the bifurcation at a time, so there are no cell-cell interactions. A range of physiologically relevant hydraulic resistivities and osmotic pressure differences are explored. Decreasing hydraulic resistivity and/or decreasing osmotic pressure differences (ESL resistance to compression) produced four behaviors: (1) RBC partitioning nonuniformity increased slightly; (2) RBC deformation decreased; (3) RBC velocity decreased relative to blood flow velocity; and (4) RBCs penetrated more deeply into the ESL. Decreasing the ESL's resistance to flow and/or compression to pathological levels could lead to more frequent cell adhesion and clotting as well as impaired vascular regulation due to weaker ATP and nitric oxide release. Potential mechanisms that can contribute to these behaviors are also discussed.

Kalman filter parameter estimation for a nonlinear diffusion model of epithelial cell migration using stochastic collocation and the Karhunen-Loeve expansion.

Several Kalman filter algorithms are presented for data assimilation and parameter estimation for a nonlinear diffusion model of epithelial cell migration. These include the ensemble Kalman filter with Monte Carlo sampling and a stochastic collocation (SC) Kalman filter with structured sampling. Further, two types of noise are considered -uncorrelated noise resulting in one stochastic dimension for each element of the spatial grid and correlated noise parameterized by the Karhunen-Loeve (KL) expansion resulting in one stochastic dimension for each KL term. The efficiency and accuracy of the four methods are investigated for two cases with synthetic data with and without noise, as well as data from a laboratory experiment. While it is observed that all algorithms perform reasonably well in matching the target solution and estimating the diffusion coefficient and the growth rate, it is illustrated that the algorithms that employ SC and KL expansion are computationally more efficient, as they require fewer ensemble members for comparable accuracy. In the case of SC methods, this is due to improved approximation in stochastic space compared to Monte Carlo sampling. In the case of KL methods, the parameterization of the noise results in a stochastic space of smaller dimension. The most efficient method is the one combining SC and KL expansion.

Simulated Red Blood Cell Motion in Microvessel Bifurcations: Effects of Cell-Cell Interactions on Cell Partitioning.

Partitioning of red blood cell (RBC) fluxes between the branches of a diverging microvessel bifurcation is generally not proportional to the flow rates, as RBCs preferentially enter the higher-flow branch. A two-dimensional model for RBC motion and deformation is used to investigate the effects of cell-cell mechanical interactions on RBC partitioning in bifurcations. The RBC membrane and cytoplasm are represented by sets of viscoelastic elements immersed in a low Reynolds number flow. Several types of two-cell interactions that can affect partitioning are found. In the most frequent interactions, a `trade-off' occurs, in which a cell entering one branch causes a following cell to enter the other branch. Other types of interactions include `herding,' where the leading cell is caused to enter the same branch as the following cell, and `following,' where the trailing cell is caused to enter the same branch as the leading cell. The combined effect of these cell-cell interactions is a tendency towards more uniform partitioning, which results from the trade-off effect but is reduced by the herding and following effects. With increasing hematocrit, the frequency of interactions increases, and more uniform partitioning results. This prediction is consistent with experimental observations on how hematocrit affects RBC partitioning.

Simulated two-dimensional red blood cell motion, deformation, and partitioning in microvessel bifurcations.

Movement, deformation, and partitioning of mammalian red blood cells (RBCs) in diverging microvessel bifurcations are simulated using a two-dimensional, flexible-particle model. A set of viscoelastic elements represents the RBC membrane and the cytoplasm. Motion of isolated cells is considered, neglecting cell-to-cell interactions. Center-of-mass trajectories deviate from background flow streamlines due to migration of flexible cells towards the mother vessel centerline upstream of the bifurcation and due to flow perturbations caused by cell obstruction in the neighborhood of the bifurcation. RBC partitioning in the bifurcation is predicted by determining the RBC fraction entering each branch, for a given partition of total flow and for a given upstream distribution of RBCs. Typically, RBCs preferentially enter the higher-flow branch, leading to unequal discharge hematocrits in the downstream branches. This effect is increased by migration toward the centerline but decreased by the effects of obstruction. It is stronger for flexible cells than for rigid circular particles of corresponding size, and decreases with increasing parent vessel diameter. For unequally sized daughter vessels, partitioning is asymmetric, with RBCs tending to enter the smaller vessel. Partitioning is not significantly affected by branching angles. Model predictions are consistent with previous experimental results.