Modeling to link regional myocardial work, metabolism and blood flows

Transmural Distribution of Coronary Perfusion and Myocardial Work Density Due to Alterations in Ventricular Loading, Geometry and Contractility

Myocardial supply changes to accommodate the variation of myocardial demand across the heart wall to maintain normal cardiac function. A computational framework that couples the systemic circulation of a left ventricular (LV) finite element model and coronary perfusion in a closed loop is developed to investigate the transmural distribution of the myocardial demand (work density) and supply (perfusion) ratio. Calibrated and validated against measurements of LV mechanics and coronary perfusion, the model is applied to investigate changes in the transmural distribution of passive coronary perfusion, myocardial work density, and their ratio in response to changes in LV contractility, preload, afterload, wall thickness, and cavity volume. The model predicts the following: (1) Total passive coronary flow varies from a minimum value at the endocardium to a maximum value at the epicardium transmurally that is consistent with the transmural distribution of IMP; (2) Total passive coronary flow at different transmural locations is increased with an increase in either contractility, afterload, or preload of the LV, whereas is reduced with an increase in wall thickness or cavity volume; (3) Myocardial work density at different transmural locations is increased transmurally with an increase in either contractility, afterload, preload or cavity volume of the LV, but is reduced with an increase in wall thickness; (4) Myocardial work density-perfusion mismatch ratio at different transmural locations is increased with an increase in contractility, preload, wall thickness or cavity volume of the LV, and the ratio is higher at the endocardium than the epicardium. These results suggest that an increase in either contractility, preload, wall thickness, or cavity volume of the LV can increase the vulnerability of the subendocardial region to ischemia.

Impact of baseline calibration on semiquantitative assessment of myocardial perfusion reserve by adenosine stress MRI

In this study, we sought to investigate the impact of baseline calibration, which is used in quantitative cardiac MRI perfusion analysis to correct for surface coil inhomogeneity and noise, on myocardial perfusion reserve index (MPRI) and its contribution to previously reported paradoxical low MPRI < 1.0 in patients with unobstructed coronary arteries. Semiquantitative perfusion analysis was performed in 20 patients with unobstructed coronary arteries undergoing stress/rest perfusion CMR and in ten patients undergoing paired rest perfusion CMR. The following baseline calibration settings were compared: (1) baseline division, (2) baseline subtraction and (3) no baseline calibration. In uncalibrated analysis, we observed ~ 20% segmental dispersion of signal intensity (SI)-over-time curves. Both baseline subtraction and baseline division reduced relative dispersion of t0-SI (p < 0.001), but only baseline division corrected for dispersion of peak-SI and maximum upslope also (p < 0.001). In the assessment of perfusion indices, however, baseline division resulted in paradoxical low MPRI (1.01 ± 0.23 vs. 1.63 ± 0.38, p < 0.001) and rest perfusion index (RPI 0.54 ± 0.07 vs. 0.94 ± 0.12, p < 0.001), respectively. This was due to a reversed ratio of blood-pool and myocardial baseline-SI before the second perfusion study caused by circulating contrast agent from the first injection. In conclusion, baseline division reliably corrects for inhomogeneity of the surface coil sensitivity profile facilitating comparisons of regional myocardial perfusion during hyperemia or at rest. However, in the assessment of MPRI, baseline division can lead to paradoxical low results (even MPRI < 1.0 in patients with unobstructed coronary arteries) potentially mimicking severely impaired perfusion reserve. Thus, in the assessment of MPRI we propose to waive baseline calibration.

Morphometric Reconstruction of Coronary Vasculature Incorporating Uniformity of Flow Dispersion

Experimental limitations in measurements of coronary flow in the beating heart have led to the development of in silico models of reconstructed coronary trees. Previous coronary reconstructions relied primarily on anatomical data, including statistical morphometry (e.g., diameters, length, connectivity, longitudinal position). Such reconstructions are non-unique, however, often leading to unrealistic predicted flow features. Thus, it is necessary to impose physiological flow constraints to ensure realistic tree reconstruction. Since a vessel flow depends on its diameter to fourth power, diameters are the logical candidates to guide vascular reconstructions to achieve realistic flows. Here, a diameter assignment method was developed where each vessel diameter was determined depending on its downstream tree size, aimed to reduce flow dispersion to within measured range. Since the coronary micro-vessels are responsible for a major portion of the flow resistance, the auto regulated coronary flow was analyzed in a morphometry-based reconstructed 400 vessel arterial microvascular sub-tree spanning vessel orders 1–6. Diameters in this subtree were re-assigned based on the flow criteria. The results revealed that diameter re-assignment, while adhering to measured morphometry, significantly reduced the flow dispersion to realistic levels while adhering to measured morphometry. The resulting network flow has longitudinal pressure distribution, flow fractal nature, and near-neighboring flow autocorrelation, which agree with measured coronary flow characteristics. Collectively, these results suggest that a realistic coronary tree reconstruction should impose not only morphometric data but also flow considerations. The work is of broad significance in providing a novel computational framework in the field of coronary microcirculation. It is essential for the study of coronary circulation by model simulation, based on a realistic network structure.

Coronary microcirculatory pathophysiology: can we afford it to remain a black box?

Coronary microvascular networks play the key role in determining blood flow distribution in the heart. Matching local blood supply to tissue metabolic demand entails continuous adaptation of coronary vessels via regulation of smooth muscle tone and structural dilated vessel diameter. The importance of coronary microcirculation for relevant pathological conditions including angina in patients with normal or near-normal coronary angiograms [microvascular angina (MVA)] and heart failure with preserved ejection fraction (HFpEF) is increasingly recognized. For MVA, clinical studies have shown a prevalence of up to 40% in patients with suspected coronary artery disease and a relevant impact on adverse cardiovascular events including cardiac death, stroke, and heart failure. Despite a continuously increasing number of corresponding clinical studies, the knowledge on pathophysiological cause–effect relations involving coronary microcirculation is, however, still very limited. A number of pathophysiological hypotheses for MVA and HFpEF have been suggested but are not established to a degree, which would allow definition of nosological entities, stratification of affected patients, or development of effective therapeutic strategies. This may be related to a steep decline in experimental (animal) pathophysiological studies in this area during the last 15 years. Since technology to experimentally investigate microvascular pathophysiology in the beating heart is increasingly, in principle, available, a concerted effort to build ‘coronary microcirculatory observatories’ to close this gap and to accelerate clinical progress in this area is suggested.

Fractal regional myocardial blood flows pattern according to metabolism, not vascular anatomy

Regional myocardial blood flows are markedly heterogeneous. Fractal analysis shows strong near-neighbor correlation. In experiments to distinguish control by vascular anatomy vs. local vasomotion, coronary flows were increased in open-chest dogs by stimulating myocardial metabolism (catecholamines + atropine) with and without adenosine. During control states mean left ventricular (LV) myocardial blood flows (microspheres) were 0.5-1 ml·g(-1)·min(-1) and increased to 2-3 ml·g(-1)·min(-1) with catecholamine infusion and to ∼4 ml·g(-1)·min(-1) with adenosine (Ado). Flow heterogeneity was similar in all states: relative dispersion (RD = SD/mean) was ∼25%, using LV pieces 0.1-0.2% of total. During catecholamine infusion local flows increased in proportion to the mean flows in 45% of the LV, "tracking" closely (increased proportionately to mean flow), while ∼40% trended toward the mean. Near-neighbor regional flows remained strongly spatially correlated, with fractal dimension D near 1.2 (Hurst coefficient 0.8). The spatial patterns remain similar at varied levels of metabolic stimulation inferring metabolic dominance. In contrast, adenosine vasodilation increased flows eightfold times control while destroying correlation with the control state. The Ado-induced spatial patterns differed from control but were self-consistent, inferring that with full vasodilation the relaxed arterial anatomy dominates the distribution. We conclude that vascular anatomy governs flow distributions during adenosine vasodilation but that metabolic vasoregulation dominates in normal physiological states.

A Green's function method for simulation of time-dependent solute transport and reaction in realistic microvascular geometries

A novel theoretical method is presented for simulating the spatially resolved convective and diffusive transport of reacting solutes between microvascular networks and the surrounding tissues. The method allows for efficient computational solution of problems involving convection and non-linear binding of solutes in blood flowing through microvascular networks with realistic 3D geometries, coupled with transvascular exchange and diffusion and reaction in the surrounding tissue space. The method is based on a Green's function approach, in which the solute concentration distribution in the tissue is expressed as a sum of fields generated by time-varying distributions of discrete sources and sinks. As an example of the application of the method, the washout of an inert diffusible tracer substance from a tissue region perfused by a network of microvessels is simulated, showing its dependence on the solute's transvascular permeability and tissue diffusivity. Exponential decay of the washout concentration is predicted, with rate constants that are about 10-30% lower than the rate constants for a tissue cylinder model with the same vessel length, vessel surface area and blood flow rate per tissue volume.

Simple accurate mathematical models of blood HbO2 and HbCO2 dissociation curves at varied physiological conditions: evaluation and comparison with other models

PURPOSE

Equations for blood oxyhemoglobin (HbO2) and carbaminohemoglobin (HbCO2) dissociation curves that incorporate nonlinear biochemical interactions of oxygen and carbon dioxide with hemoglobin (Hb), covering a wide range of physiological conditions, are crucial for a number of practical applications. These include the development of physiologically-based computational models of alveolar-blood and blood-tissue O2–CO2 transport, exchange, and metabolism, and the analysis of clinical and in vitro data.

METHODS AND RESULTS

To this end, we have revisited, simplified, and extended our previous models of blood HbO2 and HbCO2 dissociation curves (Dash and Bassingthwaighte, Ann Biomed Eng 38:1683–1701, 2010), validated wherever possible by available experimental data, so that the models now accurately fit the low HbO2 saturation (SHbO2) range over a wide range of values of PCO2, pH, 2,3-DPG, and temperature. Our new equations incorporate a novel PO2-dependent variable cooperativity hypothesis for the binding of O2 to Hb, and a new equation for P50 of O2 that provides accurate shifts in the HbO2 and HbCO2 dissociation curves over a wide range of physiological conditions. The accuracy and efficiency of these equations in computing PO2 and PCO2 from the SHbO2 and SHbCO2 levels using simple iterative numerical schemes that give rapid convergence is a significant advantage over alternative SHbO2 and SHbCO2 models.

CONCLUSION

The new SHbO2 and SHbCO2 models have significant computational modeling implications as they provide high accuracy under non-physiological conditions, such as ischemia and reperfusion, extremes in gas concentrations, high altitudes, and extreme temperatures.

Simple accurate mathematical models of blood HbO2 and HbCO2 dissociation curves at varied physiological conditions: evaluation and comparison with other models

PURPOSE

Equations for blood oxyhemoglobin (HbO2) and carbaminohemoglobin (HbCO2) dissociation curves that incorporate nonlinear biochemical interactions of oxygen and carbon dioxide with hemoglobin (Hb), covering a wide range of physiological conditions, are crucial for a number of practical applications. These include the development of physiologically-based computational models of alveolar-blood and blood-tissue O2–CO2 transport, exchange, and metabolism, and the analysis of clinical and in vitro data.

METHODS AND RESULTS

To this end, we have revisited, simplified, and extended our previous models of blood HbO2 and HbCO2 dissociation curves (Dash and Bassingthwaighte, Ann Biomed Eng 38:1683–1701, 2010), validated wherever possible by available experimental data, so that the models now accurately fit the low HbO2 saturation (SHbO2) range over a wide range of values of PCO2, pH, 2,3-DPG, and temperature. Our new equations incorporate a novel PO2-dependent variable cooperativity hypothesis for the binding of O2 to Hb, and a new equation for P50 of O2 that provides accurate shifts in the HbO2 and HbCO2 dissociation curves over a wide range of physiological conditions. The accuracy and efficiency of these equations in computing PO2 and PCO2 from the SHbO2 and SHbCO2 levels using simple iterative numerical schemes that give rapid convergence is a significant advantage over alternative SHbO2 and SHbCO2 models.

CONCLUSION

The new SHbO2 and SHbCO2 models have significant computational modeling implications as they provide high accuracy under non-physiological conditions, such as ischemia and reperfusion, extremes in gas concentrations, high altitudes, and extreme temperatures.

Heart-on-a-chip based on stem cell biology

Heart diseases are one of the main causes of death around the world. The great challenge for scientists is to develop new therapeutic methods for these types of ailments. Stem cells (SCs) therapy could be one of a promising technique used for renewal of cardiac cells and treatment of heart diseases. Conventional in vitro techniques utilized for investigation of heart regeneration do not mimic natural cardiac physiology. Lab-on-a-chip systems may be the solution which could allow the creation of a heart muscle model, enabling the growth of cardiac cells in conditions similar to in vivo conditions. Microsystems can be also used for differentiation of stem cells into heart cells, successfully. It will help better understand of proliferation and regeneration ability of these cells. In this review, we present Heart-on-a-chip systems based on cardiac cell culture and stem cell biology. This review begins with the description of the physiological environment and the functions of the heart. Next, we shortly described conventional techniques of stem cells differentiation into the cardiac cells. This review is mostly focused on describing Lab-on-a-chip systems for cardiac tissue engineering. Therefore, in the next part of this article, the microsystems for both cardiac cell culture and SCs differentiation into cardiac cells are described. The section about SCs differentiation into the heart cells is divided in sections describing biochemical, physical and mechanical stimulations. Finally, we outline present challenges and future research concerning Heart-on-a-chip based on stem cell biology.

Applications of computational models to better understand microvascular remodelling: a focus on biomechanical integration across scales

Microvascular network remodelling is a common denominator for multiple pathologies and involves both angiogenesis, defined as the sprouting of new capillaries, and network patterning associated with the organization and connectivity of existing vessels. Much of what we know about microvascular remodelling at the network, cellular and molecular scales has been derived from reductionist biological experiments, yet what happens when the experiments provide incomplete (or only qualitative) information? This review will emphasize the value of applying computational approaches to advance our understanding of the underlying mechanisms and effects of microvascular remodelling. Examples of individual computational models applied to each of the scales will highlight the potential of answering specific questions that cannot be answered using typical biological experimentation alone. Looking into the future, we will also identify the needs and challenges associated with integrating computational models across scales.

Tuneable resolution as a systems biology approach for multi-scale, multi-compartment computational models

The use of multi-scale mathematical and computational models to study complex biological processes is becoming increasingly productive. Multi-scale models span a range of spatial and/or temporal scales and can encompass multi-compartment (e.g., multi-organ) models. Modeling advances are enabling virtual experiments to explore and answer questions that are problematic to address in the wet-lab. Wet-lab experimental technologies now allow scientists to observe, measure, record, and analyze experiments focusing on different system aspects at a variety of biological scales. We need the technical ability to mirror that same flexibility in virtual experiments using multi-scale models. Here we present a new approach, tuneable resolution, which can begin providing that flexibility. Tuneable resolution involves fine- or coarse-graining existing multi-scale models at the user's discretion, allowing adjustment of the level of resolution specific to a question, an experiment, or a scale of interest. Tuneable resolution expands options for revising and validating mechanistic multi-scale models, can extend the longevity of multi-scale models, and may increase computational efficiency. The tuneable resolution approach can be applied to many model types, including differential equation, agent-based, and hybrid models. We demonstrate our tuneable resolution ideas with examples relevant to infectious disease modeling, illustrating key principles at work. WIREs Syst Biol Med 2014, 6:225–245. doi:10.1002/wsbm.1270

How to cite this article:WIREs Syst Biol Med 2014, 6:289–309. doi:10.1002/wsbm.1270

Left ventricular finite element model bounded by a systemic circulation model

A series of models were developed in which a circulatory system model was coupled to an existing series of finite element (FE) models of the left ventricle (LV). The circulatory models were used to provide realistic boundary conditions for the LV models. This was developed for the JSim analysis package and was composed of a systemic arterial, capillary, and venous system in a closed loop with a varying elastance LV and left atria to provide the driving pressures and flows matching those of the FE model. Three coupled models were developed, a normal LV under normotensive aortic loading (116/80 mm Hg), a mild hypertension (137/89 mm Hg) model, and a moderate hypertension model (165/100 mm Hg). The initial step in the modeling analysis was that the circulation was optimized to the end-diastolic pressure and volume values of the LV model. The cardiac FE models were then optimized to the systolic pressure/volume characteristics of the steady-state JSim circulatory model solution. Comparison of the stress predictions for the three models indicated that the mild hypertensive case produced a 21% increase in the average fiber stress levels, and the moderate hypertension case had a 36% increase in average stress. The circulatory work increased by 18% and 43% over that of the control for the mild and moderate hypertensive cases, respectively.

Are the different patterns of stress-induced (Takotsubo) cardiomyopathy explained by regional mechanical overload and demand: supply mismatch in selected ventricular regions?

Takotsubo cardiomyopathy (TCM) or stress-induced cardiomyopathy is an increasingly recognized syndrome characterized by severe regional left ventricular dysfunction in the absence of an explanatory coronary lesion. TCM may lead to lethal complications but is completely reversible if the patient survives the acute phase. The pathogenesis of TCM and the mechanism behind this remarkable recovery are unknown. Plasma levels of catecholamine are elevated in many TCM patients and exogenously administered catecholamine induces TCM-like cardiac dysfunction in both humans and rats. A catecholamine excess increases myocardial metabolic demand by increasing the force of contraction as well as the heart rate, and also alters cardiac depolarization patterns. We propose that an altered spatiotemporal pattern of cardiac contraction and excessive force of contraction may lead to a redistribution of wall stresses in the left ventricle. This redistribution of wall stress causes regional mechanical overload of regions where wall tension becomes disproportionately great and renders these cardiomyocytes "metabolically insufficient". In other words, these cardiomyocytes experience a demand: supply mismatch on the basis of excessive metabolic demand. In order to prevent the death of these cardiomyocytes and to prevent excessive wall tension from developing in neighboring regions, a protective metabolic shutdown occurs in the affected cardiomyocytes. This metabolic shutdown, i.e., acute down regulation of non-vital cellular functions, serves to protect the affected regions from necrosis and explains the apparently complete recovery observed in TCM. We propose that this phenomenon may share important characteristics with phenomena such as ischemic conditioning, stunning and hibernation. In this manuscript, we discuss our hypothesis in the context of available knowledge and discuss important experiments that would help to corroborate or refute the hypothesis.