ORIGINAL RESEARCH article
Sec. Computational Physiology and Medicine
Effects of Mechanical Dyssynchrony on Coronary Flow: Insights From a Computational Model of Coupled Coronary Perfusion With Systemic Circulation
- 1Department of Mechanical Engineering, Michigan State University, East Lansing, MI, United States
- 2California Medical Innovation Institute, San Diego, CA, United States
Mechanical dyssynchrony affects left ventricular (LV) mechanics and coronary perfusion. Due to the confounding effects of their bi-directional interactions, the mechanisms behind these changes are difficult to isolate from experimental and clinical studies alone. Here, we develop and calibrate a closed-loop computational model that couples the systemic circulation, LV mechanics, and coronary perfusion. The model is applied to simulate the impact of mechanical dyssynchrony on coronary flow in the left anterior descending artery (LAD) and left circumflex artery (LCX) territories caused by regional alterations in perfusion pressure and intramyocardial pressure (IMP). We also investigate the effects of regional coronary flow alterations on regional LV contractility in mechanical dyssynchrony based on prescribed contractility-flow relationships without considering autoregulation. The model predicts that LCX and LAD flows are reduced by 7.2%, and increased by 17.1%, respectively, in mechanical dyssynchrony with a systolic dyssynchrony index of 10% when the LAD's IMP is synchronous with the arterial pressure. The LAD flow is reduced by 11.6% only when its IMP is delayed with respect to the arterial pressure by 0.07 s. When contractility is sensitive to coronary flow, mechanical dyssynchrony can affect global LV mechanics, IMPs and contractility that in turn, further affect the coronary flow in a feedback loop that results in a substantial reduction of dPLV/dt, indicative of ischemia. Taken together, these findings imply that regional IMPs play a significant role in affecting regional coronary flows in mechanical dyssynchrony and the changes in regional coronary flow may produce ischemia when contractility is sensitive to the changes in coronary flow.
Heart failure (HF) is a complex syndrome associated with significant morbidity, mortality and socioeconomic burden. Many HF patients are presented with mechanical dyssynchrony (Claridge et al., 2015) that is a strong predictor of cardiovascular mortality (Saporito and Fisher, 1987; Modin et al., 2018). Cardiac resynchronization therapy (CRT), on the other hand, has emerged as a powerful treatment for mechanical dyssynchrony and HF (Moss et al., 2009). About 30% of patients, however, still do not improve after therapy (non-responders) and the percentage of CRT non-responders has remained constant over the past decade (Gorcsan, 2011). The relatively high-rate of CRT non-responders has imposed a significant burden on the clinical course and outcome of patients (Corbisiero et al., 2016).
Studies have suggested that CRT response is sensitive to coronary perfusion (Trimble et al., 2008; Svendsen et al., 2012) and mechanical dyssynchrony elicits changes in coronary perfusion (Claridge et al., 2015). During normal electrical activation of the left ventricle (LV) in which contraction of the interventricular septum and LV free wall (LVFW) are synchronized, there is sufficient blood flow to balance the regional myocardial demand and supply (Borlaug and Kass, 2008). In left bundle branch block (LBBB), however, there is a delay in the electrical activation of the LVFW with respect to the septum that can produce acute changes in the LV function and coronary blood flow to the left anterior descending artery (LAD) (septal hypoperfusion) and left circumflex artery (LCX) territories in the LV (Fang et al., 2013; Claridge et al., 2015). Right ventricular (RV) pacing can also elicit mechanical dyssynchrony in which the septum and LVFW contract asynchronously (Heyndrickx et al., 1985), which also produces several abnormalities resembling that of LBBB [e.g., delayed LV relaxation, reduced cardiac output (Zile et al., 1987; Xiao et al., 1992), and depressed LV function (Shefer et al., 1987; Kolettis et al., 2000)]. Similar to LBBB, RV pacing can produce perfusion defects in the form of a reduction in LAD flow. This reduction has been found in patients without coronary artery disease (Matsumura et al., 1987) as well as in an animal model of RV pacing (Bailer, 1981). One study in LBBB patients with normal coronary arteries had an abnormal increase in lactate production during a high-rate of pacing, suggesting that perfusion changes associated with mechanical dyssynchrony may induce ischemia (Breithardt and Breithardt, 2012). It is, however, unclear if this ischemia, plausibly induced by changes in perfusion during mechanical dyssynchrony can reduce myocardial contractility, which in turn, causes further changes to the coronary flow that may exacerbate ischemia. Due to the confounding effects of the two-way interactions between LV mechanics and coronary flow during mechanical dyssynchrony, it is difficult to isolate the mechanisms behind these changes from experimental and clinical measurements alone.
To overcome the limitations associated with experiments and clinical studies of mechanical dyssynchrony (e.g., difficulty in decoupling the confounding effects between LV mechanics and coronary flow), computational models have been developed. Specifically, finite element (FE) models and lumped parameter models [e.g., TriSeg/CircAdap model (Arts et al., 2005)] of ventricular mechanics (Trayanova et al., 2011; Kerckhoffs et al., 2012; Sermesant et al., 2012; Walmsley et al., 2015; Lee et al., 2018; Arumugam et al., 2019) have been developed and applied to investigate global and regional LV mechanics (i.e., stresses and strains) and energetics (i.e., work done) in mechanical dyssynchrony. The impact of mechanical dyssynchrony on coronary flow, however, has not been incorporated, and correspondingly, interrogated in these models. This is in part because of the complexity and the need for a large number of parameters necessary to model the coupling between coronary perfusion and ventricular mechanics in a theoretical framework (Beard and Bassingthwaighte, 2000; Smith et al., 2002; Huo et al., 2009; Namani et al., 2018a). Moreover, most computational models of coronary perfusion (Beard and Bassingthwaighte, 2000; Smith et al., 2002; Huo et al., 2009; Namani et al., 2018a) are not coupled in a closed-loop manner with the systemic circulation and the heart. As a result, it is difficult to apply these models to directly investigate the bi-directional interaction mechanisms affecting the LV mechanics-coronary perfusion relationship, which is not completely understood in both normal and pathophysiological settings. This is especially significant since changes in the LV mechanics can affect coronary perfusion (via intramyocardial and perfusion pressures) and changes in coronary flow can affect LV mechanics (via changes in myocardial contractility) in a closed-loop feedback manner. Computational models that account for bi-directional interactions between the LV and the coronary network can, therefore, help discriminate the roles and contributions of the various confounding factors affecting coronary flow in mechanical dyssynchrony.
Motivated by the limitations of current models and the need to better understand the impact of mechanical dyssynchrony on coronary flow, we developed a novel computational framework that couples the systemic circulation and the LV's contraction-relaxation mechanics with the coronary flow dynamics in both the LAD and LCX territories. The framework is calibrated using experimental measurements of the coronary flow waveforms in the LAD and the LCX, pressure waveforms in the LV, as well as the LV volume waveform from a swine model under normal conditions. The calibrated model is then extended to model mechanical dyssynchrony with and without considering the effects of regional coronary flow on regional LV contractility (ischemia). A sensitivity analysis is performed by varying key parameters in the modeling framework [such as time delays in the regional intramyocardial pressure (IMP) waveforms] to investigate how they affect changes in the LAD and the LCX flows, as well as the global LV function in mechanical dyssynchrony. In addition, we investigate the role of the regional contractility—coronary flow relationship in producing regional ischemia that affects the (global and regional) LV mechanics, IMP, and contractility in mechanical dyssynchrony.
All experiments were performed following the national and local ethical guidelines, including the Institute of Laboratory Animal Research guidelines, the Public Health Service policy, and the Animal Welfare Act, and an approved California Medical Innovations Institute IACUC protocol regarding the use of animals in research. Yorkshire domestic swine (n = 3), of either sex, were fasted overnight and sedated with TKX (Telazol 10 mg/kg, Ketamine 5 mg/kg, Xylazine 5 mg/kg; IM). The surgical plane was maintained with 1–2% Isoflurane and 100% O2 to keep PCO2 at 35–40 mmHg. Limb and precordial ECG leads were attached to the animals and cardiac electrical activity was monitored on a Physio-Control Lifepak-12 defibrillator. Two introducer sheaths were percutaneously inserted into the jugular vein for administration of fluids and placement of a pacing lead, and to access to the right side of the heart. An additional introducer sheath was percutaneously inserted into the right femoral artery to access the coronary arteries and the LV. Heparin 100 IU/kg was administered with IV fluids before further instrumentation of the animal, which was then supplemented as needed. Measurements were recorded with right atrial (RA) pacing at 100 beats per minute.
Left Ventricular Pressure
After proper calibration, a Millar pressure catheter (Ventri Cath 507), which was connected to an MPVS Ultra PV loop system (Millar, Houston, Texas, USA), was advanced into the LV using a guide catheter. Pressure waveforms were recorded using a data acquisition system (LabChart Pro, ADInstruments, Colorado Springs, CO, USA).
Transesophageal echocardiograms were obtained simultaneously at the time of LV pressure measurements, using an EPIQ 7C ultrasound system (Philips, Andover, MA) with an X8-2t transducer. Two- and three-dimensional echocardiographic images in the long axis view were acquired with the animal placed in the supine position. The images were analyzed offline using QLAB 10.8 (Philips, Andover, MA) to determine ejection fraction (EF), LV volume waveform, cardiac output, and stroke volume.
Coronary Flow Rate
The chest was opened through a midline sternotomy and an incision was made in the pericardium with the creation of a sling to support the heart. The LAD and LCX arteries were carefully dissected free from their surrounding tissue and a flow probe connected to a flow meter (Transonic, Ithaca, NY, USA) was placed around the arteries to measure their flow rates during RA pacing.
Cardiac-Coronary Modeling Framework
The closed-loop lumped parameter modeling framework consists of the LV, systemic and coronary circulations to describe the bi-directional interactions between LV mechanics and coronary perfusion. A schematic of the modeling framework is given in Figure 1A, showing the connection between the LV and systemic circulations with coronary network.
Figure 1. (A) Schematic of the cardiac-coronary modeling framework: the LV and the coronary tree (LAD and LCX) are connected to a closed-loop lumped model with their pressures (P), volumes (V), flow rate (Q), compliances (C), and resistances (R) where ao, art, ven, and mv denote the aorta, peripheral arteries, peripheral veins, and mitral valves, respectively; (B) A representation of the model to separate the LV into LAD and LCX compartments with their pressures (P), volumes (V), maximal chamber elastances (Ees), and time-varying elastances (e); (C) Time-varying elastance function in the LAD and LCX, eLAD and eLCX in mechanical dyssynchrony. Note, the LAD and LCX networks are graphically depicted by the 3D coronary tree in the lower left figure.
The systemic circulation consists of four key compartments: LV, left atrium (LA), peripheral arterial and venous networks that are each represented by their corresponding electrical analogs (Figure 1A) (Shavik et al., 2017, 2018). Total blood volume is conserved by relating the rate of volume change in each storage compartment in the circulatory system to the inflow and outflow rates as
In Equations (1)–(4), VLA, VLV, Vart, and Vven are the volumes of the four compartments with the subscripts denoting the LA, LV, peripheral arteries, and peripheral veins, respectively. In Equations (3) and (4), qLAD and qLCX are the coronary flow rates associated with the LAD and LCX territories, respectively. Flow rates in the venous, mitral valve, aortic, and peripheral arteries are denoted by qven, qmv, qao, and qper, which are each determined by their segment's resistance and the pressure difference between the two connecting storage compartments as
where PLV, Part, Pven, and PLA are the cavity pressures, and Rao, Rper, Rven, and Rmv are the resistances of the aorta, peripheral arteries, peripheral veins, and mitral valves, respectively. We note that the flows across the heart valves are unidirectional (assuming no valvular regurgitation) (Shavik et al., 2019). Pressures in the storage compartments representing the peripheral arterial and venous networks are given as
where (Vart0, Vven0) are the prescribed (arteries and veins) resting volumes, and (Cart, Cven) are the prescribed total (arteries and veins) compliances. The pumping characteristics of the LV and LA are represented by a time-varying elastance model (Santamore and Burkhoff, 1991; Witzenburg and Holmes, 2018). For the LA, instantaneous pressure, PLA is related to the instantaneous volume, VLA by a time-varying elastance function, eLA(t) as follows
In Equation (11), Pes and Ped are the end-systolic and end-diastolic pressures, respectively, Ees,LA is the maximal chamber elastance of the LA, VLA0 is the volume with zero end-systolic pressure, and both ALA and BLA are parameters defining the end-diastolic pressure volume relationship of the LA. In Equation (14), Tmax,LA is the time to end systole and τLA is the relaxation time constant.
In order to simulate asynchronous activation of the septum associated with the LAD territory and the LVFW associated with the LCX territory, we divided the LV cavity volume (VLV) into two compartments with volumes (VLV,LAD) and (VLV,LCX) based on the approach by Sunagawa et al. (1983) (i.e., VLV = VLV, LAD + VLV,LCX) (Figure 1B) and assumed the pressures in the compartments associated with the LAD and LCX to be equal to the LV cavity pressure (PLV). In Sunagawa's study on ischemia (Sunagawa et al., 1983), LV pressure was derived by prescribing different Ees in the normal and ischemic regions. Here, based on the same assumption (i.e., PLV = PLAD = PLCX), different time-varying elastance functions (eLAD, eLCX) and different maximal chamber elastances (Ees,LAD, Ees,LCX) are prescribed for the LAD and LCX territories. To simulate mechanical dyssynchrony, we prescribed a positive time delay (Δt > 0) in the time-varying elastance function of the LCX with respect to the LAD [i.e., eLCX = e(t + Δt) and eLAD = e(t)] so that the contraction of LCX territory happens later than that in the LAD territory with respect to end of diastole (Figure 1C). The resultant LV cavity pressure is given as
where F = VLV,LAD/VLV is a prescribed value, and VLV0 is the LV volume at zero end-systolic pressure (refer to Appendix A for derivation). We note that the zero end-systolic pressure corresponds to the volume intercept of the end-systolic pressure volume relationship (ESPVR), which we have assumed to be linear (Davidson et al., 2015). In Equation (15), Pes, Ped, e(t), and e(t + Δt) are obtained from Equations (12)–(14) by replacing the subscript LA with LV.
Coronary network morphometry of the LAD and LCX territories are each described using a previous coronary arterial model (Namani et al., 2018a). Morphometry of the arterial tree in that model is pruned from a previously reconstructed coronary arterial network and consists of 400 vessels (from order 6 to order 0) with 195 bifurcations, 3 trifurcations, and 79 terminal vessels (Kaimovitz et al., 2010). The microvascular tree, an input to the flow analysis, was generated from the global statistics of the measured coronary vessel diameters (Kassab et al., 1993; Kaimovitz et al., 2005, 2010). Because there are many capillary networks at the distal end of the coronary arterial tree, the terminal vessels of the reconstructed tree are truncated with either pre-capillary arterioles of vessel order 1 and/or capillary arterioles of vessel order 0. Larger arterioles are designated as vessels of orders 1, 2, 3, …, 6 in the direction of increasing diameter (Kassab et al., 1993). To maintain physiological flow dispersion levels, the diameter of each vessel order was reassigned within the range of the measurement in the reconstructed microvascular tree as described previously (Namani et al., 2018b). To take into account the differences in the morphometry (e.g., diameter and length of the vessels) and flow in the LAD vs. the LCX, we prescribed different proximal resistances RLAD and RLCX in the corresponding segments connecting the coronary networks to the systemic circulation (Figure 1A) (Kassab et al., 1993). To couple the coronary perfusion (both the LAD and the LCX networks) with the systemic circulation, the arterial pressure (Part) and venous pressure (Pven) are employed as the inlet and outlet pressure boundary conditions of the coronary arterial tree, based on the fact that (1) most of the pressure drop occurs across the coronary arterial microvascular networks (Marcus et al., 1990) and (2) the systemic venous pressure (~10 mm Hg) is of the same magnitude as the RA pressure that the networks drain into under physiological conditions. We also note that the venous pressure is imposed at the outlet pressure in the “intramyocardial pump” model introduced by Spaan et al. (1981). The model is solved for coronary flows in the LAD (qLAD) and the LCX (qLCX) trees, which are fed into the systemic circulation model in Equations (3) and (4). In the flow analysis, it is assumed that the diameter is homogeneous within each vessel segment but different from vessel to vessel. The non-linear three-element Windkessel electrical representation that has been investigated previously (Namani et al., 2018a) is employed in this work for coronary network analysis. Modeling the flow in a trifurcation is analogous to that in a bifurcation (Namani et al., 2018a), with the difference being that an additional vessel needs to be accounted for in the mass conservation. For a single bifurcation (j = 3) or trifurcation network (j = 4), mass conservation is satisfied at the nodal position as
where Qi is the flow in each vessel (Figure 2A). For each single vessel i, using mass conservation at Windkessel's internal (“mid”) node, the ordinary differential equation in terms of the flow is
where Pmid is the unknown pressure in the geometric center of the vessel, PT is the intramyocardial pressure with details found in the next section, the resistances R1 and R2 (assumed to be the same) and the capacitor C are obtained by
Figure 2. (A) Single bifurcation three-vessel network; (B) Electrical non-linear analog of a single vessel.
In Equation (18), L is the length and D(t) is the diameter of single coronary vessel, respectively, and μ denotes the effective blood viscosity, which is assumed to be a constant to simplify the flow analysis (Pries et al., 1994) (Figure 2B). The capacitance of an intramyocardial arteriole is non-linear and defined by Equation (19) where the partial derivative means the change in lumen volume (V) with respect to the transvascular pressure when the other variables affecting the lumen volume are not varied (Jacobs et al., 2008). The vessel diameter is a sigmoidal function of transvascular pressure as
where Ap and Bp are the asymptotical highest and lowest radii, respectively (under the highest and lowest transvascular pressure), φp is the transvascular pressure corresponding to the mean of radii Ap and Bp, and Cp is the passive response bandwidth. For the parameters used in Equation (20), refer to Namani et al. (2018a). As the resistance and capacitor in the network is not a constant and depend on the flow conditions, the resulting system of ordinary differential equations is non-linear.
Myocardium Vessel Interaction
In the model, intramyocardial pressure in each vessel (i.e., extravascular forces on the coronary vessels) of the LAD and LCX arterial trees are prescribed to be homogeneous and have the value of IMPLAD and IMPLCX, respectively, that are defined as follows. Because the distribution of IMP in vivo is difficult to measure (Heineman and Grayson, 1985; Mihailescu and Abel, 1994), some studies have assumed that the IMP equals the LV pressure at the endocardium and decreases linearly to zero at the epicardium (Downey and Kirk, 1975; Bellamy, 1982). On the other hand, it has been suggested that the LV pressure is not the sole determinant of IMP but myocardial contraction (associated with changing myocardial elasticity and myocyte shortening) also affects the coronary flow myocardium-vessel interaction mechanisms (Westerhof, 1990; Westerhof et al., 2006): (1) time-varying elasticity (VE), (2) myocardial shortening-induced intracellular pressure (SIP), and (3) ventricular cavity-induced extracellular pressure (CEP) (Algranati et al., 2010). Correspondingly, the intramyocardial pressure is prescribed as
where α, β, and γ are prescribed parameters, and SSR is the myocardial stretch. These parameters are prescribed based on a sensitivity analysis (refer to Appendix E) to (1) best match the coronary flow rate waveform and (2) be consistent with experiments. Specifically, α is constrained between 0 and 1 based on previous modeling studies (Arts and Reneman, 1985; Bruinsma et al., 1988) that assume CEP to vary linearly from the cavity pressure PLV at the endocardium to zero at the epicardium. The parameter γ is chosen to produce a peak SIP that is about 20% of IMP based on a previous study (Mynard et al., 2014). For the SIP component in Equation (21), SSR is defined as the ratio of VLV to the LV's end-diastolic volume (Algranati et al., 2010). The sum of peak SIP and VE produced by any combination of α, β and γ should be <30% of IMP based on the experimental data by Rabbany et al. (1989), which shows that the peak IMP is 10–30% greater than the peak PLV. Based on these conditions and sensitivity analysis, we arrived at the values α = 0.8, β = 5, and γ = 20. Here, we note that the different time-varying elastance functions (eLAD, eLCX) prescribed to simulate mechanical dyssynchrony will result in different IMPs on the LAD and the LCX networks.
Simulations for Mechanical Dyssynchrony and Ischemia
We consider three types of simulations: (1) control (CON) based on measurements from three RA pacing swine models, (2) isolated mechanical dyssynchrony (MD), and (3) mechanical dyssynchrony + ischemia (MD + IS). In the control simulations (without mechanical dyssynchrony) (n = 3), model parameters were calibrated to fit the coronary flow (LAD, LCX), LV volume and pressure, and arterial pressure waveforms measured in the animals. Specifically, the peak LV pressure and end-diastolic LV pressure were matched by manually adjusting the aorta and mitral valve resistances, (Rao and Rmv), and the passive parameters defining the end-diastolic pressure volume relationship (ALV and BLV). The arterial compliance (Cart), the peripheral arterial resistance (Rper), the arterial resting volume (Vart0), and the initial value for the arterial volume (Vart) were manually adjusted to match both the peak and the end systolic arterial pressures. The LV volume waveform was matched by manually adjusting the resting LV volume (VLV0), the initial value of the LV volume, (VLV) and the venous vessel compliance (Vven). For the control simulations, Ees,LAD was set equal to Ees,LCX and Δt was set to zero. Therefore, the LAD and LCX networks are exposed to the same IMP, and flows in the LAD and LCX were matched by altering their corresponding inlet resistances (RLAD and RLCX). Based on the initial guess of the volume in each compartment (Vven, Vart, VLV, VLA), the simulation is run until the LV pressure, volume, and coronary flow rate waveforms are at steady-state time-periodic conditions. Reference parameters used in previous studies on lumped-parameter models (Shavik et al., 2019; Namani et al., 2020) are selected and manually varied after each run until the waveforms are in agreement with the experimental data. We note that parameters for the coronary model (section Coronary Circulation) were obtained from a previous study and fixed here whereas parameters associated with IMP (α, β and γ) were obtained based on a sensitivity study as discussed earlier. Only parameters associated with the lumped parameter model of the systemic circulation and the proximal resistances (RLAD and RLCX) of the coronary arterial tree are calibrated based on the experimental measurements. The parameters from the calibrated model for the three swine model are listed in Appendix B (Table B1). We note that VLV0 (volume intercept of the prescribed linear ESPVR) is negative in Table B1 because the time-varying elastance model does not take into account the non-linearity of ESPVR at small VLV. We further note that a negative VLV0 has been found experimentally and numerically (Crottogini et al., 1987; Pironet et al., 2013; Davidson et al., 2015). Parameters that are fixed for all simulations are listed with their value in Table 1.
The calibrated model is then applied to study the independent effect of mechanical dyssynchrony (MD) and its combined effect with ischemia (MD + IS) on LV function. Since there is a delayed activation in the LCX in the MD simulations (n = 3), we prescribed a systolic dyssynchrony index (SDI > 0), which is defined by the time delay Δt (Δt = duration of one cardiac cycle × SDI) to model the delayed contraction (see Equation 15) of the LCX with respect to the LAD territory. We note that SDI is a common metric used to assess mechanical dyssynchrony in the clinic (Delgado et al., 2008; Tani et al., 2012). Additionally, we also prescribed a different time delay for both the elastance functions [i.e., and in the IMP function in Equation (21) to simulate the effects of time delay on the contraction of the LAD territory with respect to PLV (whose elastance function in Equation 15 is not altered).
In the MD + IS simulations (n = 3), we include the bi-directional interaction between the coronary flow and the LV contractility. To model this interaction, we prescribed a relationship between contractility and the coronary flow separately in the LAD and LCX territories based on experimental observations (Figure 3) (Gallagher et al., 1984; Ross, 1991; Heusch and Schulz, 1999). This relationship consists of two regimes, an ischemic and a non-ischemic one that are separated by the transitional flows ∑ QLAD, n and ∑ QLCX, n in the LAD and the LCX, respectively. The transitional flows ∑ QLAD, n and ∑ QLCX, n are prescribed, respectively, by the corresponding total coronary flow (over a cardiac cycle) measured in the normal swine model (“n” denotes to normal). In the ischemic regime, maximal chamber elastances of the LAD (Ees,LAD) and LCX (Ees,LCX) territories vary linearly with their corresponding total coronary flows (over a cardiac cycle) ∑ QLAD and ∑ QLCX, respectively. The ischemic regime is controlled by the slope k, which has different values in the LAD and LCX (Figure 3). In the non-ischemic regime, the maximal chamber elastances remains constant with respect to their corresponding coronary flows (Heusch, 2019).
Figure 3. The prescribed contractility-coronary flow relationship in both the ischemic and non-ischemic myocardial regimes. The parameter “k” represents the slope of this relationship in the ischemic regime and “n” denotes normal condition.
This linear relationship between contractility and coronary flow in the ischemic regime is constructed based on ample and robust experimental evidence of “perfusion-contraction matching” in the ischemic myocardium as described by Schulz et al. (1989) and Ross (1991). Considering the prescribed contractility-flow relationship, coronary flow alterations caused by mechanical dyssynchrony affect regional contractility (Ees,LAD, Ees,LCX) that further affects both IMPs in the LAD and the LCX networks, and the coronary networks' inlet pressures, which changes the coronary flow again in a closed-loop feedback. To elucidate the complex effects quantitatively, we consider two cases with different values of the contractility-flow gradient k.
We note that the parameters given in Table 1 and the animal specific parameters in Table B1 are fixed for the three simulation cases. All simulations are run with a heart rate at 100 beats per minute for several cardiac cycles until the pressure-volume loop reaches a steady state.
Experimental measurements and model predicted results of the coronary flow in LAD (∑ QLAD) and LCX (∑ QLCX), LV end-diastolic volume (LV EDV), and peak LV pressure (peak LVP) from three swine are summarized in Table 2. The LV EDV and peak LVP are within the normal range as compared to a previous work (Lionetti et al., 2013). Coronary flow measured in the LAD is larger than that in the LCX in the three swine models. This is consistent with previous studies (Seiler et al., 1992; Miguel, 1999), where the epicardial vessel in the LAD network has been observed to have a larger diameter than that in the LCX network. This allows for more flow in the LAD artery as volumetric flow is proportional to the fourth power of the vessel's diameter.
The arterial pressure, LV pressure, LAD, and LCX flow rates, and LV volume predicted by the simulations are in agreement with the experimental measurements (see Table 2 and a representative case in Figure 4). The average error between the model predictions and measurements for the three animals for ∑ QLAD, ∑ QLCX, LV EDV, and peak LVP are 0, 0, 0.4, and 0.2%, respectively. While there are some slight discrepancies in the model predicted time-course (e.g., in the LV volume waveform and the magnitude of the first up-shoot during systole in the LAD coronary flow rate waveform) when compared to the experimental measurements, the calibrated model, nevertheless, can reproduce key features of the measured LV volume, pressure and coronary flow rate waveforms (Figures 4D,E). Specifically, the model predicted a reduced coronary flow rate in the LCX and LAD networks during early systole (t = 0–0.06 s) because as IMP starts to increase, Part has yet to increase during this period of time (Figures 4A,F). The flow rate then oscillates (t = 0.06–0.33 s) as the perfusion pressure, which is largely driven by Part, rises and competes with IMP to control the flow. With the decay of IMP during diastole (t > 0.33 s), the perfusion pressure dominates, and the flow rate increases. Three components of IMP associated with the myocardium-vessel interaction mechanism show that the contribution of each mechanism is different (i.e., mean CEP, ~76%, mean SIP, ~18% and mean VE, ~6%) (Figure 4C). We note that since the LCX territory (LVFW) and LAD territory (LV septum) contract synchronously in the control simulations, they both have the same IMP waveforms (i.e., IMPLAD = IMPLCX) due to their identical VE and SIP components.
Figure 4. A representative simulation is compared with data from the control group. Experimentally (Exp.) and numerically (Num.) predicted waveforms of (A) LV pressure; (B) LV volume; (C) IMP and its three components (Num. only); (D) qLAD; (E) qLCX; and (F) arterial pressure (Num. only). “S” and “D” denote systole and diastole, respectively.
Isolated Mechanical Dyssynchrony
From the simulations for isolated mechanical dyssynchrony (SDI > 0), the waveforms for the LV pressure, LV volume, LCX flow rate, and arterial volume are all shifted to the right compared to the control simulations (SDI = 0) as a result of a time delay, Δt in the contraction of the LCX territory with respect to the LAD territory (Figures 5A,B,E,F). The LAD flow rate waveform is, however, less shifted (Figure 5D). The model predicted a larger EDV because the delayed contraction causes a delay in the rise of the LV pressure, which in turn results in a delay of the onset of end-diastole. Because of the larger EDV, peak PLV is increased (based on Equation 15). Mean value of PLV is, however, reduced compared to that of the control group because of the delay in the rise of PLV during systole (Figure 5A). As a result, the peak value of the CEP component increases but its mean value decreases (Figure 5C). On the other hand, the VELAD waveform remains synchronous with the CEP waveform while the VELCX waveform is dyssynchronous but with the same magnitude. Correspondingly, the peak IMPLAD is greater while the peak IMPLCX is smaller than the corresponding values in the control simulations (Figures 4C, 5C). Because of the wave shift of IMPs and perfusion pressure that is largely driven by Part (during t = 0−0.12s), there is a corresponding wave shift in terms of qLAD and qLCX in which the first up shoot during systole is delayed by ~0.07s (Figures 5D,E) Due to the greater increase of Part than IMPLAD during t = 0.12 − 0.22 s, coronary flow in the LAD increases in the isolated mechanical dyssynchrony simulations. While due to the decay of peak IMPLCX resulted from the time delay, qLCX decreases in late systole and early diastole, resulting in a reduction of the LCX coronary flow in the isolated mechanical dyssynchrony simulations.
Figure 5. A representative simulation of the control case (Δt = 0) and isolated mechanical dyssynchrony case (t ≠ 0 with SDI = 10%). Comparison of (A) LV pressure; (B) LV volume; (C) IMP and its three components at the LAD and the LCX for only in isolated mechanical dyssynchrony; (D) qLAD; (E) qLCX; and (F) arterial pressure. Dotted and dash lines denote a control and an isolated dyssynchrony with both systole (S) and diastole (D).
These effects of mechanical dyssynchrony on Part, PLV, qLAD, qLCX, IMPLAD, and IMPLCX are all consistent in the animal-specific simulations of the three swine models. Coronary flow and IMP in the LAD and LCX, as well as the mean and peak values of PLV and Part all vary monotonically with increasing degree of mechanical dyssynchrony as SDI increases from 5 to 15% (Figure 6A). Specifically, ∑ QLAD, peak PLV, peak Part, and peak IMPLAD increase with SDI while all the other quantities decrease with SDI. We note that peak IMPLAD increases while the mean IMPLAD decreases with increasing SDI as explained earlier. Therefore, the disparity between LAD and LCX flows increases with increasing SDI. For an SDI of 15%, ∑ QLAD is increased by about 28% whereas ∑ QLCX is decreased by about 12% compared to control simulation (Figure 6A). Given that the coronary flow rate is a competition between perfusion pressure and IMP, the percentage change in mean Part (~1% at SDI = 15%) relative to that of the IMPs (~3% in LAD and ~6% LCX at SDI = 15%) at a fixed SDI in mechanical dyssynchrony implies that the changes in coronary flow rates are affected more by the changes in regional IMPs than the perfusion pressure.
Figure 6. Effects of SDI on total coronary flow, mean arterial and LV pressures, and IMPLAD and IMPLCX over three swine. (A) Percentage difference of the quantities with respect to the control simulation for different SDI. (B) Percentage difference of the quantities with respect to the control simulation for different with SDI = 10%.
With a time delay prescribed to the VE and CEP components of both the IMPLAD and IMPLCX in Equation (21), the simulation shows that at an SDI = 10%, LAD flow is reduced (Figure 6B). This is because IMPLAD is delayed with respect to PLV with , which reduces the flow. With Δ, there is more delay in IMPLCX, which impedes more flow during early diastole. Increasing the time delay also reduces the peak IMPLAD as VELAD and CEP becomes asynchronous. With s, ∑ QLAD and ∑ QLCX are decreased by 12 and 20% with respect to the control simulations (Figure 6B).
Mechanical Dyssynchrony With Ischemia
In the mechanical dyssynchrony with ischemia (MD + IS) simulations for an SDI = 10%, both Ees,LAD and Ees,LCX are prescribed as a function of ∑ QLAD and ∑ QLCX, respectively, based on a contractility-flow relationship (Figure 3). The gradients in ischemic regions (kLAD, kLCX) were obtained by taking the ratio of the maximal chamber elastance under control condition (Ees = 3.18) with the experimentally measured coronary flow (∑ QLAD, n, ∑ QLCX, n) for each swine (Table 2). The values of (kLAD, kLCX) for the three swine models are (0.18, 0.34), (0.27, 0.38), and (0.17, 0.22) mmHg/ml2. Coupling of the contractility with coronary flow for these values of k produced only a small decrease in ∑ QLAD (0.21%) and ∑ QLCX (4.1%) when compared to the isolated mechanical dyssynchrony case (Figure 7A). For these values of k, coupling of the contractility with coronary flow, which changes during mechanical dyssynchrony, did not produce substantial changes in peak PLV, Part, and IMPLAD. The coupling produced only a slightly reduced IMPLCX (4.1%) compared to that in the isolated mechanical dyssynchrony case (Figure 7B). The IMPLCX is reduced because of a reduction in ∑ QLCX (4.1%) in mechanical dyssynchrony that decreases Ees,LCX (contractility in LCX), which in turn reduces the VE component of IMPLCX.
Figure 7. A comparison between the simulations for the control case, the isolated mechanical dyssynchrony case, and the mechanical dyssynchrony (SDI = 10%) + ischemia case for k = 0.34 and 10k = 3.40 mmHg/ml2 for the following variables: (A) Coronary flow in the LAD and LCX with k; (B) Arterial Pressure, LV pressure waveforms and IMPs with k; (C) Coronary flow in the LAD and LCX with 10k; (D) Arterial Pressure, LV pressure waveforms and IMPs with 10k; and (E) Percentage differences over three swine using different slopes (k) in the contractility-flow relationship in ischemic region of LCX in terms of flow, dPLV/dt and mean arterial and LV pressures between the isolated mechanical dyssynchrony and mechanical dyssynchrony + ischemia cases.
Increasing the gradient k in the ischemic regime of the contractility-flow relationship for LCX by 10 times produced pronounced changes in the hemodynamics (Figures 7C,D). A steeper contractility-flow relationship produced a substantial reduction in Ees,LCX, which decreases the overall LV contractility and led to a substantial reduction in PLV and Part (Figure 7D) that in turn significantly decreases flows to the LAD and LCX territory (Figure 7C). The decrease in LAD and LCX flows further reduced the overall LV contractility until the pressure waveforms reached a steady periodic state. Averaging the effects of incorporating the different contractility-flow relationships in the three swine models (Figure 7E), we find that the resultant maximum rate of LV pressure rise (positive and negative dPLV/dt) which is an index of contractility is reduced by 45–55% in mechanical dyssynchrony with ischemia when k is large.
The key findings of this study are: (1) an increase in time delay Δt (SDI) in the IMPLCX waveform with respect to IMPLAD, LV, and arterial pressure waveforms produces a decrease in LCX flow and an increase in LAD flow; (2) a time delay in both IMPLAD and IMPLCX waveform with respect to the LV and arterial pressure waveforms causes both LAD and LCX coronary flow to decrease; (3) prescribing a direct relationship between regional contractility and its corresponding coronary flow (Figure 3) can affect global LV mechanics, hemodynamics and function, which in turn further affects the coronary flow in a feedback loop in mechanical dyssynchrony, especially when contractility is highly sensitive to flow (i.e., large gradient k in the ischemic regime of the relationship). These findings imply that during mechanical dyssynchrony, regional IMP plays a significant role in affecting regional coronary flows, which can then affect global LV function in a feedback loop that produces ischemia if contractility is sensitive to the coronary flow changes. Correspondingly, correcting for mechanical dyssynchrony by CRT may help restore regional IMP waveform and coronary flow that improve both regional and global LV contractility.
We developed a closed-loop cardiac-coronary modeling framework consisting of the systemic and the coronary circulations to interrogate the bi-directional interactions between LV mechanics and coronary flow in mechanical dyssynchrony. The modeling framework was calibrated using the LV pressure and volume, LAD and LCX flow rate waveforms that were measured experimentally in three (normal) RA pacing swine models. We then applied the calibrated framework to investigate the isolated effects of mechanical dyssynchrony on coronary flow, as well as the effects of mechanical dyssynchrony when regional coronary flow and contractility are coupled. We show that regional coronary flows in the LCX and LAD are altered by asynchronous LV contraction due to changes in the arterial pressure, LV pressure, IMPLAD and IMPLCX waveforms as discussed below.
We show that a delay or phase shift in the IMPLCX waveform with respect to the IMPLAD waveform that is synchronous with the LV and arterial pressures produces a reduction in LCX flow, but an increase in LAD flow. To produce a reduction in LAD flow (in addition to a reduction in LCX flow) requires a delay or phase shift in IMPLAD with respect to the LV and arterial pressures. These model predictions are consistent with experimental studies, which have shown that in RV pacing or LBBB (Ono et al., 1992; Skalidis et al., 1999, 2001; Kyriakides et al., 2007), regional myocardial blood flow in the septum is significantly reduced compared to the reduction in LVFW after LBBB induction. Studies on canine model with pacing at different ventricular sites (Delhaas et al., 1994; Amitzur et al., 1995) also show a significant reduction in regional blood flow in the early activated regions compared to the late activated regions. Specifically, Amitzur et al. (1995) shows that RV pacing causes a reduction in LAD flow without any change in the LCX flow when perfusion pressure is fixed (17–45%) or not fixed (~10%) as compared to RA pacing (representing normal activation). Under maximal vasodilation using adenosine, however, they reported no changes in LAD flow with RV pacing. In contrast to that study, a few clinical studies (Skalidis et al., 2001; Itoh et al., 2015) show that flow (as measured by flow velocity time integral) is reduced at maximal vasodilation in both the LCX and LAD arteries in LBBB patients compared to when they are treated with CRT (Itoh et al., 2015) or to control subjects (Skalidis et al., 2001). This discrepancy between the experimental and clinical studies may be attributed to anatomical differences in the coronary vasculature between humans and canine (Vernooy et al., 2007). Our model prediction that both LAD and LCX flow changes with mechanical dyssynchrony is consistent with clinical studies (Skalidis et al., 2001; Itoh et al., 2015).
Because mechanical dyssynchrony produces regional changes in mechanics (e.g., regional myocardial fiber strain, work done and wall tension) leading to a redistribution of myocardial oxygen demand (Prinzen et al., 1999) [e.g., RV pacing is found to reduce and increase mechanical work in the septum and the LVFW, respectively (Prinzen et al., 1990, 1999)], some studies have hypothesized that the regional alteration in coronary flow found in mechanical dyssynchrony is caused by regional changes in metabolic demand (Prinzen et al., 1990). These alterations, however, persist even under hyperemia or maximal vasodilation in the absence of autoregulation mechanisms (i.e., metabolic, myogenic, and shear regulations) (Prinzen et al., 1990) as described in the above studies. Indeed, a study (Itoh et al., 2015) has shown that both LAD and LCX flows are moderately reduced by about 14 and 6% in asynchronous activation under hyperemia, which are within the range of reductions predicted by our model when IMPLAD is delayed by about (Figure 6B). These findings further suggest that regional changes in IMP and perfusion pressure play a significant role on the regional changes in coronary flow during mechanical dyssynchrony. The mechanisms of IMP are attributed to the contraction of the myocyte (by the changing of myocardial tissue elasticity and shortening induced pressure) and the cavity pressure (Algranati et al., 2010). Some studies (Skalidis et al., 1999, 2001; Kyriakides et al., 2007) suggest that asynchronous activation due to RV pacing or LBBB causes an increase in IMP in the interventricular septum that reduces myocardial perfusion but the effects on IMPLCX has not been measured. Here, in our model prediction, the regional variation of LV mechanics in mechanical dyssynchrony produces regional changes in IMP (IMPLAD increases and IMPLCX decreases) that in turn alters the coronary flow. Taken together with the clinical and experimental studies mentioned above, our finding underscores the significance of the regional IMP in regulating the regional coronary flow during mechanical dyssynchrony.
Because myocardial tissue contractility depends on coronary flow (Ross, 1991), we hypothesized that regional variation in coronary flow can in turn affect the regional contractility, especially when coronary flow is impeded. To interrogate this effect, we prescribed (for the first time) a regional contractility-flow relationship (Ees − ∑ Q) in the modeling framework based on experimental observations (Ross, 1991). The model predicts that the reduction in LCX coronary flow during mechanical dyssynchrony led to a decrease in Ees,LCX that affects LV and arterial pressures as well as IMPLAD and IMPLCX, causing further effects on the coronary flow. Since the contractility is highly sensitive to the coronary flow in the ischemic regime (i.e., large gradient k) of the contractility-flow relationship, a small reduction in the regional coronary flow in mechanical dyssynchrony leads to a large decrease in the corresponding regional and overall LV contractility as indicated by a reduction in dPLV/dt, a metric associated with myocardial contractility (Theroux et al., 1974). The reduction in overall LV contractility in turn reduces the perfusion pressure and regional IMP that further reduces coronary flow in a feedback loop. These results suggest that mechanical dyssynchrony can produce ischemia (i.e., a reduction in global LV contractility due to a reduction in regional coronary flow), especially when contractility-flow relationship has a steep gradient. Existing evidence of ischemia caused solely by asynchronous activation from previous studies are conflicting (Kuhn et al., 1982; Ono et al., 1992). While an experimental study has reported no change in lactate production with the induction of LBBB in a canine model at resting heart rate (Ono et al., 1992), a clinical study measuring lactate production in the myocardium in 95 patients with no coronary disease and a normal coronary angiogram has suggested that mechanical dyssynchrony may induce myocardial ischemia at a high-rate of pacing (Kuhn et al., 1982). The latter study is supported by the findings that LBBB patients often complain of exertional dyspnea despite normal left ventricular dimensions and function at rest, which might imply the induction of ischemia (Breithardt and Breithardt, 2012). Although we have not taken into account coronary autoregulation due to changes in metabolic demand that should be present in these studies (especially at higher heart rate), our model suggests that the local changes in coronary flow may produce global changes in LV function that is indicative of ischemia (i.e., reduced dPLV/dt). This is especially so when the myocardial contractility is very sensitive to flow conditions, which may occur when metabolic regulation is exhausted.
While modeling of coronary perfusion using the poroelastic theory in a finite element framework (Chapelle et al., 2010; Lee et al., 2016) and 1D Navier-Stokes equation (Smith et al., 2002) can take into account more details associated with coronary flow (e.g., myocardial tissue microstructure, geometrical variation in the tissue, wave reflection, etc.), lumped models of the coronary vasculature (such as this work) (Wang et al., 1989; Duanmu et al., 2018; Namani et al., 2020) sacrifice these details for computational efficiency in simulating mechanical and hemodynamic changes associated with heart diseases. To the best of our knowledge, we have, for the first time, bi-directionally coupled the coronary circulation with LV mechanics in a closed-loop manner to model the effects of mechanical dyssynchrony on both the LV's perfusion and its mechanics. Bi-directional interactions between the regional LV mechanics and the coronary vasculature in this novel modeling framework occur both directly (IMP affects flow and flow affects contractility) and indirectly (LV mechanics affects perfusion pressure that in turn affects flow). By dividing the LV into two compartments associated with the LAD and LCX regions based on the idea from Sunagawa et al. (1983), we were able to simulate the corresponding regional changes in coronary flow associated with asynchronous activation of the perfused territories. By prescribing a contractility-flow relationship in the LAD and the LCX regions, we were also able to simulate, for the first time, how those changes affect the regional LV contractility that in turn affects the LV hemodynamics globally in a feedback loop. The calibrated flow rate waveform in the LAD predicted by our model (0D) (Figure 4A) is also comparable with the predictions of other published coronary flow models (Liang et al., 2007; Mynard et al., 2014; Duanmu et al., 2018), including models based on the 1D Navier-Stokes equation (refer to Appendix F for the comparison). Given these capabilities, our model offers a valuable and efficient approach to understand the mechanisms that affect coronary flow in mechanical dyssynchrony that maybe useful for applying CRT in heart failure. In particular when taken together with the consistent observation that LVFW mechanical work is increased in mechanical dyssynchrony, our model predictions imply the presence of a stronger compensatory control mechanism to increase (or at least maintain) coronary blood flow to the LVFW, further suggesting that defects in autoregulation may play a central role in discriminating CRT responder from non-responder.
There are some limitations associated with this study. First, we did not consider the autoregulation mechanism in coronary flow analysis. As such, the flow predicted by the model is only driven by perfusion pressure and IMP under passive vessel's condition. Coronary autoregulation has been shown to maintain blood flow relatively constant across a wide range of perfusion pressure (Johnson, 1986; Namani et al., 2018a). This mechanism can be incorporated into the existing framework in future studies to investigate the role of autoregulation in asynchronous activation. Second, 0D lumped model used in this study is able to predict the main features of coronary flow rate waveforms but unable to capture the wave oscillations. This issue can be addressed by coupling with a 1D model. Third, LV mechanics is simplified using a time-varying elastance model by a rule-of-mixture based on the approach by Sunagawa et al. (1983). Although this approach is able to represent key characteristics of the LV, it does not take into account detailed mechanical properties (e.g., anisotropy), microstructure of the LV (e.g., myofiber orientation) and the length-dependent effects in the cardiac muscles. These limitations may account for some of the discrepancies between model predicted and experimentally measured LV pressure and volume waveforms in terms of the time course. Nevertheless, we have shown that the model predictions are consistent with the experimental results in reproducing key features of the waveforms. This limitation can be overcome by using a finite element model for describing LV mechanics (Arumugam et al., 2019). Fourth, we have not included the RV in the modeling framework. Although the RV pressure is small compared to the LV, it may affect the cavity pressure component of the IMP in the LAD territory, which was assumed to be the same as in the LCX territory here. To address this issue, we have performed a sensitivity analysis with a different CEP component in the LAD that considers the effects of RV pressure (refer to Appendix G for sensitivity analysis). The analysis show that the conclusion is not altered when we consider the effects of RV pressure in the LAD IMP. The effects of RV on IMP can be investigated in greater detail in future studies using finite element modeling. Fifth, because the goal here is not to infer and interpret any changes in the model parameters arising from the calibration of measurements associated with the experimental/clinical data, we did not perform a rigorous sensitivity analysis of all the parameters and uncertainty quantification analysis. Such analyses will be performed in future studies. Last, model predictions of RV pacing were compared with experimental data from previous published studies as mentioned in the discussion, but not with animal-specific measurements. This can be addressed in future studies with animal-specific measurements made under RV and RA pacing.
In summary, we have developed and calibrated a closed-loop cardiac-coronary modeling framework that takes into account the bi-directional interaction between LV mechanics, systemic circulation and coronary perfusion. We applied the modeling framework to simulate mechanical dyssynchrony and show that LV mechanics and coronary flow are affected by asynchronous contraction and interacts bi-directionally. We show that the changes in coronary flow are sensitive to a phase shift in the regional IMP waveform that may arise from asynchronous activation. We further show that, depending on the regional contractility-flow relationship, the altered regional coronary flow can influence the global LV hemodynamics and regional LV mechanics that in turn further affects the coronary flow to induce myocardial ischemia.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
The animal study was reviewed and approved by the Institutional Animal Care and Use Committee.
LF, RN, GK, and LL contributed to the conception and design of the study. JC contributed to the conduction experimental procedure. LF, RN, JC, GK, and LL contributed to the analysis and interpretation of data, the drafting of the manuscript, revising it critically for important intellectual content, and final approval of the manuscript. All authors have significantly contributed to the submitted work.
This work was supported by NIH (R01 HL134841).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2020.00915/full#supplementary-material
Amitzur, G., Manor, D., Pressman, A., Adam, D., Hammerman, H., Shofti, R., et al. (1995). Modulation of the arterial coronary blood flow by asynchronous activation with ventricular pacing. Pacing Clin. Electrophysiol. 18, 697–710. doi: 10.1111/j.1540-8159.1995.tb04664.x
Arts, T., Delhaas, T., Bovendeerd, P., Verbeek, X., and Prinzen, F. W. (2005). Adaptation to mechanical load determines shape and properties of heart and circulation: the circadapt model. Am. J. Physiol. Heart Circ. Physiol. 288, 1943–1954. doi: 10.1152/ajpheart.00444.2004
Bellamy, R. F. (1982). Comments on “diastolic-systolic coronary flow differences are caused by intramyocardial pump action in the anesthetized dog”. Circ. Res. 50, 747–748. doi: 10.1161/01.RES.50.5.747
Bruinsma, P., Arts, T., Dankelman, J., and Spaan, J. A. E. (1988). Model of the coronary circulation based on pressure dependence of coronary resistance and compliance. Basic Res. Cardiol. 83, 510–524. doi: 10.1007/BF01906680
Chapelle, D., Gerbeau, J. F., Sainte-Marie, J., and Vignon-Clementel, I. E. (2010). A poroelastic model valid in large strains with applications to perfusion in cardiac modeling. Comput. Mech. 46, 91–101. doi: 10.1007/s00466-009-0452-x
Claridge, S., Chen, Z., Jackson, T., Sammut, E., Sohal, M., Behar, J., et al. (2015). Current concepts relating coronary flow, myocardial perfusion and metabolism in left bundle branch block and cardiac resynchronisation therapy. Int. J. Cardiol. 181, 65–72. doi: 10.1016/j.ijcard.2014.11.194
Corbisiero, R., Buck, D. C., Muller, D., Bharmi, R., Dalal, N., and Kazemian, P. (2016). What is the cost of non-response to cardiac resynchronization therapy? Hospitalizations and healthcare utilization in the CRT-D population. J. Interv. Card. Electrophysiol. 47, 189–195. doi: 10.1007/s10840-016-0180-z
Crottogini, A. J., Willshaw, P., Barra, J. G., Armentano, R., Cabrera Fischer, E. I., and Pichel, R. H. (1987). Inconsistency of the slope and the volume intercept of the end-systolic pressure-volume relationship as individual indexes of inotropic state in conscious dogs: presentation of an index combining both variables. Circulation 76, 1115–1126. doi: 10.1161/01.CIR.76.5.1115
Davidson, S. M., Kannangara, D. O., Pretty, C. G., Kamoi, S., Pironet, A., Desaive, T., et al. (2015). Modelling of the nonlinear end-systolic pressure-volume relation and volume-at-zero-pressure in porcine experiments. Proc. Annu. Int. Conf. IEEE Eng. Med. Biol. Soc. EMBS 2015, 6544–6547. doi: 10.1109/EMBC.2015.7319892
Delgado, V., Sitges, M., Vidal, B., Silva, E., Azqueta, M., Tolosana, J. M., et al. (2008). Assessment of left ventricular dyssynchrony by real-time three-dimensional echocardiography. Rev. Española Cardiol. 61, 825–834. doi: 10.1157/13124993
Delhaas, T., Arts, T., Prinzen, F. W., and Reneman, R. S. (1994). Regional fibre stress-fibre strain area as an estimate of regional blood flow and oxygen demand in the canine heart. J. Physiol. 477, 481–496. doi: 10.1113/jphysiol.1994.sp020209
Fang, F., Chan, J. Y. S., Lee, A. P. W., Sung, S. H., Luo, X. X., Jiang, X., et al. (2013). Improved coronary artery blood flow following the correction of systolic dyssynchrony with cardiac resynchronization therapy. Int. J. Cardiol. 167, 2167–2171. doi: 10.1016/j.ijcard.2012.05.094
Gallagher, K. P., Matsuzaki, M., and Koziol, J. A. (1984). Regional myocardial perfusion and wall thickening during ischemia in conscious dogs. Am. J. Physiol. Heart Circ. Physiol. 247, H727–H738. doi: 10.1152/ajpheart.1984.247.5.H727
Heineman, F. W., and Grayson, J. (1985). Transmural distribution of intramyocardial pressure measured by micropipette technique. Am. J. Physiol. Heart Circ. Physiol. 249, H1216–H223. doi: 10.1152/ajpheart.1985.249.6.H1216
Heusch, G. (2019). Myocardial ischemia: lack of coronary blood flow, myocardial oxygen supply-demand imbalance, or what? Am. J. Physiol. Heart Circ. Physiol. 316, H1439–H1446. doi: 10.1152/ajpheart.00139.2019
Heyndrickx, G. R., Vilaine, J. P., Knight, D. R., and Vatner, S. F. (1985). Effects of altered site of electrical activation on myocardial performance during inotropic stimulation. Circulation 71, 1010–1016. doi: 10.1161/01.CIR.71.5.1010
Huo, Y., Kaimovitz, B., Lanir, Y., Wischgoll, T., Hoffman, J. I. E., and Kassab, G. S. (2009). Biophysical model of the spatial heterogeneity of myocardial flow. Biophys. J. 96, 4035–4043. doi: 10.1016/j.bpj.2009.02.047
Itoh, M., Shinke, T., Yoshida, A., Kozuki, A., Takei, A., Fukuzawa, K., et al. (2015). Reduction in coronary microvascular resistance through cardiac resynchronization and its impact on chronic reverse remodelling of left ventricle in patients with non-ischaemic cardiomyopathy. Europace 17, 1407–1414. doi: 10.1093/europace/euu361
Kaimovitz, B., Lanir, Y., and Kassab, G. S. (2005). Large-scale 3-D geometric reconstruction of the porcine coronary arterial vasculature based on detailed anatomical data. Ann. Biomed. Eng. 33, 1517–1535. doi: 10.1007/s10439-005-7544-3
Kaimovitz, B., Lanir, Y., and Kassab, G. S. (2010). A full 3-D reconstruction of the entire porcine coronary vasculature. Am. J. Physiol. Heart Circ. Physiol. 299, 1064–1076. doi: 10.1152/ajpheart.00151.2010
Kerckhoffs, R. C. P., Omens, J. H., and McCulloch, A. D. (2012). Mechanical discoordination increases continuously after the onset of left bundle branch block despite constant electrical dyssynchrony in a computational model of cardiac electromechanics and growth. Europace 14, 65–72. doi: 10.1093/europace/eus274
Kolettis, T. M., Kyriakides, Z. S., Tsiapras, D., Popov, T., Paraskevaides, I. A., and Kremastinos, D. T. (2000). Improved left ventricular relaxation during short-term right ventricular outflow tract compared to apical pacing. Chest 117, 60–64. doi: 10.1378/chest.117.1.60
Kuhn, H., Losse, B., and Hart, W. (1982). “Studies in patients with abnormalelectrocardiogram of unknown etiology,” in Detection of Ischemic Myocardium With Exercise, eds F. Loogen and L. Seipel (Berlin; Heidelberg: Springer), 122–124. doi: 10.1007/978-3-642-68387-9_15
Lee, A. W. C., Costa, C. M., Strocchi, M., Rinaldi, C. A., and Niederer, S. A. (2018). Computational modeling for cardiac resynchronization therapy. J. Cardiovasc. Transl. Res. 11, 92–108. doi: 10.1007/s12265-017-9779-4
Lee, J., Nordsletten, D., Cookson, A., Rivolo, S., and Smith, N. (2016). In silico coronary wave intensity analysis: application of an integrated one-dimensional and poromechanical model of cardiac perfusion. Biomech. Model. Mechanobiol. 15, 1535–1555. doi: 10.1007/s10237-016-0782-5
Liang, F., Taniguchi, H., and Liu, H. (2007). A multi-scale computational method applied to the quantitative evaluation of the left ventricular function. Comput. Biol. Med. 37, 700–715. doi: 10.1016/j.compbiomed.2006.06.011
Lionetti, V., Romano, S. L., Bianchi, G., Bernini, F., Dushpanova, A., Mascia, G., et al. (2013). Impact of acute changes of left ventricular contractility on the transvalvular impedance: validation study by pressure-volume loop analysis in healthy pigs. PLoS ONE 8:e0080591. doi: 10.1371/journal.pone.0080591
Marcus, M. L., Chilian, W. M., Kanatsuka, H., Dellsperger, K. C., Eastham, C. L., and Lamping, K. G. (1990). Understanding the coronary circulation through studies at the microvascular level. Circulation 82, 1–7. doi: 10.1161/01.CIR.82.1.1
Matsumura, K., Ozaki, M., Yamagishi, T., Ishine, K., Furutani, Y., Nagano, H., et al. (1987). Thallium myocardial perfusion defects and wall motion abnormalities during right ventricular pacing in patients with normal coronary arteries. Kakuigaku 24, 955–962.
Mihailescu, L. S., and Abel, F. L. (1994). Intramyocardial pressure gradients in working and nonworking isolated cat hearts. Am. J. Physiol. Heart Circ. Physiol. 266, H1233–H1241. doi: 10.1152/ajpheart.1994.266.3.H1233
Modin, D., Biering-Sørensen, S. R., Møgelvang, R., Jensen, J. S., and Biering-Sørensen, T. (2018). Prognostic importance of left ventricular mechanical dyssynchrony in predicting cardiovascular death in the general population. Circ. Cardiovasc. Imaging 11:e007528. doi: 10.1161/CIRCIMAGING.117.007528
Moss, A. J., Hall, W. J., Cannom, D. S., Klein, H., Brown, M. W., Daubert, J. P., et al. (2009). Cardiac-resynchronization therapy for the prevention of heart-failure events. N. Engl. J. Med. 361, 1329–1338. doi: 10.1056/NEJMoa0906431
Mynard, J. P., Penny, D. J., and Smolich, J. J. (2014). Scalability and in vivo validation of a multiscale numerical model of the left coronary circulation. Am. J. Physiol. Heart Circ. Physiol. 306, 517–528. doi: 10.1152/ajpheart.00603.2013
Namani, R., Kassab, G. S., and Lanir, Y. (2018a). Integrative model of coronary flow in anatomically based vasculature under myogenic, shear, and metabolic regulation. J. Gen. Physiol. 150, 145–168. doi: 10.1085/jgp.201711795
Namani, R., Kassab, G. S., and Lanir, Y. (2018b). Morphometric reconstruction of coronary vasculature incorporating uniformity of flow dispersion. Front. Physiol. 9:1069. doi: 10.3389/fphys.2018.01069
Namani, R., Lee, L. C., Lanir, Y., Kaimovitz, B., Shavik, S. M., and Kassab, G. S. (2020). Effects of myocardial function and systemic circulation on regional coronary perfusion. J. Appl. Physiol. 128, 1106–1122. doi: 10.1152/japplphysiol.00450.2019
Ono, S., Nohara, R., Kambara, H., Okuda, K., and Kawai, C. (1992). Regional myocardial perfusion and glucose metabolism in experimental left bundle branch block. Circulation 85, 1125–1131. doi: 10.1161/01.CIR.85.3.1125
Pironet, A., Desaive, T., Kosta, S., Lucas, A., Paeme, S., Collet, A., et al. (2013). A multi-scale cardiovascular system model can account for the load-dependence of the end-systolic pressure-volume relationship. Biomed. Eng. Online 12:8. doi: 10.1186/1475-925X-12-8
Pries, A. R., Secomb, T. W., Geßner, T., Sperandio, M. B., Gross, J. F., and Gaehtgens, P. (1994). Resistance to blood flow in microvessels in vivo. Circ. Res. 75, 904–915. doi: 10.1161/01.RES.75.5.904
Prinzen, F. W., Augustijn, C. H., Arts, T., Allessie, M. A., and Reneman, R. S. (1990). Redistribution of myocardial fiber strain and blood flow by asynchronous activation. Am. J. Physiol. Heart Circ. Physiol. 259, H300–H308. doi: 10.1152/ajpheart.1990.259.2.H300
Prinzen, F. W., Hunter, W. C., Wyman, B. T., and McVeigh, E. R. (1999). Mapping of regional myocardial strain and work during ventricular pacing: experimental study using magnetic resonance imaging tagging. J. Am. Coll. Cardiol. 33, 1735–1742. doi: 10.1016/S0735-1097(99)00068-6
Rabbany, S. Y., Kresh, J. Y., and Noordergraaf, A. (1989). Intramyocardial pressure: interaction of myocardial fluid pressure and fiber stress. Am. J. Physiol. Hear. Circ. Physiol. 257, H357–H364. doi: 10.1152/ajpheart.1989.257.2.H357
Santamore, W. P., and Burkhoff, D. (1991). Hemodynamic consequences of ventricular interaction as assessed by model analysis. Am. J. Physiol. Heart Circ. Physiol. 260, H146–H157. doi: 10.1152/ajpheart.1991.260.1.H146
Saporito, M., and Fisher, L. D. (1987). Bundle branch block in patients with chronic coronary artery disease: angiographic correlates and prognostic significance. J. Am. Coll. Cardiol. 10, 73–80. doi: 10.1016/S0735-1097(87)80162-6
Schulz, R., Miyazaki, S., Miller, M., Thaulow, E., Heusch, G., Ross, J., et al. (1989). Consequences of regional inotropic stimulation of ischemic myocardium on regional myocardial blood flow and function in anesthetized swine. Circ. Res. 64, 1116–1126. doi: 10.1161/01.RES.64.6.1116
Seiler, C., Kirkeeide, R. L., and Gould, K. L. (1992). Basic structure-function relations of the epicardial coronary vascular tree: Basis of quantitative coronary arteriography for diffuse coronary artery disease. Circulation 85, 1987–2003. doi: 10.1161/01.CIR.85.6.1987
Sermesant, M., Chabiniok, R., Chinchapatnam, P., Mansi, T., Billet, F., Moireau, P., et al. (2012). Patient-specific electromechanical models of the heart for the prediction of pacing acute effects in CRT: a preliminary clinical validation. Med. Image Anal. 16, 201–215. doi: 10.1016/j.media.2011.07.003
Shavik, S. M., Jiang, Z., Baek, S., and Lee, L. C. (2018). High spatial resolution multi-organ finite element modeling of ventricular-arterial coupling. Front. Physiol. 9:119. doi: 10.3389/fphys.2018.00119
Shavik, S. M., Wall, S. T., Sundnes, J., Burkhoff, D., and Lee, L. C. (2017). Organ-level validation of a cross-bridge cycling descriptor in a left ventricular finite element model: effects of ventricular loading on myocardial strains. Physiol. Rep. 5, 1–14. doi: 10.14814/phy2.13392
Shavik, S. M., Zhong, L., Zhao, X., and Lee, L. C. (2019). In-silico assessment of the effects of right ventricular assist device on pulmonary arterial hypertension using an image based biventricular modeling framework. Mech. Res. Commun. 97, 101–111. doi: 10.1016/j.mechrescom.2019.04.008
Shefer, A., Rozenman, Y., David, Y., Ben, F.lugelman, M. Y., Gotsman, M. S., and Lewis, B. S. (1987). Left ventricular function during physiological cardiac pacing: relation to rate, pacing mode, and underlying myocardial disease. Pacing Clin. Electrophysiol. 10, 315–325. doi: 10.1111/j.1540-8159.1987.tb05971.x
Skalidis, E. I., Kochiadakis, G. E., Koukouraki, S. I., Chrysostomakis, S. I., Igoumenidis, N. E., Karkavitsas, N. S., et al. (2001). Myocardial perfusion in patients with permanent ventricular pacing and normal coronary arteries. J. Am. Coll. Cardiol. 37, 124–129. doi: 10.1016/S0735-1097(00)01096-2
Skalidis, E. I., Kochiadakis, G. E., Koukouraki, S. I., Parthenakis, F. I., Karkavitsas, N. S., and Vardas, P. E. (1999). Phasic coronary flow pattern and flow reserve in patients with left bundle branch block and normal coronary arteries. J. Am. Coll. Cardiol. 33, 1338–1346. doi: 10.1016/S0735-1097(98)00698-6
Spaan, J. A. E., Breuls, N. P., and Laird, J. D. (1981). Diastolic-systolic coronary flow differences are caused by intramyocardial pump action in the anesthetized dog. Circ. Res. 49, 584–593. doi: 10.1161/01.RES.49.3.584
Sunagawa, K., Maughan, W. L., and Sagawa, K. (1983). Effect of regional ischemia on the left ventricular end-systolic pressure-volume relationship of isolated canine hearts. Circ. Res. 52, 170–178. doi: 10.1161/01.RES.52.2.170
Svendsen, M., Prinzen, F. W., Das, M. K., Berwick, Z., Rybka, M., Tune, J. D., et al. (2012). Bi-ventricular pacing improves pump function only with adequate myocardial perfusion in canine hearts with pseudo-left bundle branch block. Exp. Biol. Med. 237, 644–651. doi: 10.1258/ebm.2012.012023
Tani, T., Sumida, T., Tanabe, K., Ehara, N., Yamaguchi, K., Kawai, J., et al. (2012). Left ventricular systolic dyssynchrony index by three-dimensional echocardiography in patients with decreased left ventricular function: Comparison with tissue doppler echocardiography. Echocardiography 29, 346–352. doi: 10.1111/j.1540-8175.2011.01577.x
Theroux, P., Franklin, D., Ross, J., and Kemper, W. S. (1974). Regional myocardial function during acute coronary artery occlusion and its modification by pharmacologic agents in the dog. Circ. Res. 35, 896–908. doi: 10.1161/01.RES.35.6.896
Trimble, M. A., Borges-Neto, S., Velazquez, E. J., Chen, J., Shaw, L. K., Pagnanelli, R., et al. (2008). Emerging role of myocardial perfusion imaging to evaluate patients for cardiac resynchronization therapy. Am. J. Cardiol. 102, 211–217. doi: 10.1016/j.amjcard.2008.03.043
Vernooy, K., Cornelussen, R. N. M., Verbeek, X. A. A. M., Vanagt, W. Y. R., Van Hunnik, A., Kuiper, M., et al. (2007). Cardiac resynchronization therapy cures dyssynchronopathy in canine left bundle-branch block hearts. Eur. Heart J. 28, 2148–2155. doi: 10.1093/eurheartj/ehm207
Walmsley, J., Arts, T., Derval, N., Bordachar, P., Cochet, H., Ploux, S., et al. (2015). Fast simulation of mechanical heterogeneity in the electrically asynchronous heart using the multipatch module. PLoS Comput. Biol. 11, 1–23. doi: 10.1371/journal.pcbi.1004284
Witzenburg, C. M., and Holmes, J. W. (2018). Predicting the time course of ventricular dilation and thickening using a rapid compartmental model. J. Cardiovasc. Transl. Res. 11, 109–122. doi: 10.1007/s12265-018-9793-1
Xiao, H. B., Brecker, S. J. D., and Gibson, D. G. (1992). Effects of abnormal activation on the time course of the left ventricular pressure pulse in dilated cardiomyopathy. Heart 68, 403–407. doi: 10.1136/hrt.68.10.403
Keywords: systemic circulation, coronary perfusion, mechanical dyssynchrony, ischemia, bi-directional interactions
Citation: Fan L, Namani R, Choy JS, Kassab GS and Lee LC (2020) Effects of Mechanical Dyssynchrony on Coronary Flow: Insights From a Computational Model of Coupled Coronary Perfusion With Systemic Circulation. Front. Physiol. 11:915. doi: 10.3389/fphys.2020.00915
Received: 30 April 2020; Accepted: 08 July 2020;
Published: 14 August 2020.
Edited by:Sanjay Ram Kharche, University of Western Ontario, Canada
Reviewed by:Hao Gao, University of Glasgow, United Kingdom
Stephen Michael Moore, IBM Research, Australia
Copyright © 2020 Fan, Namani, Choy, Kassab and Lee. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lei Fan, email@example.com