Preconditioning with selective autoretroperfusion: In vivo and in silico evidence of washout hypothesis.

Introduction: Prompt reperfusion of coronary artery after acute myocardial infarction (AMI) is crucial for minimizing heart injury. The myocardium, however, may experience additional injury due to the flow restoration itself (reperfusion injury, RI). The purpose of this study was to demonstrate that short preconditioning (10 min) with selective autoretroperfusion (SARP) ameliorates RI, based on a washout hypothesis. Methods: AMI was induced in 23 pigs (3 groups) by occluding the left anterior descending (LAD) artery. In SARP-b (SARP balloon inflated) and SARP-nb (SARP balloon deflated) groups, arterial blood was retroperfused for 10 min via the great cardiac vein before releasing the arterial occlusion. A mathematical model of coronary circulation was used to simulate the SARP process and evaluate the potential washout effect. Results: SARP restored left ventricular function during LAD occlusion. Ejection fraction in the SARP-b group returned to baseline levels, compared to SARP-nb and control groups. Infarct area was significantly larger in the control group than in the SARP-b and SARP-nb groups. End-systolic wall thickness was preserved in the SARP-b compared to the SARP-nb and control groups. Analyte values (pH, lactate, glucose, and others), measured every 2 min during retroperfusion, suggest a "washout" effect as one important mechanism of action of SARP in reducing infarct size. With SARP, the values progressively approached baseline levels. The mathematical model also confirmed a possible washout effect of tracers. Discussion: RI can be ameliorated by delaying restoration of arterial flow for a brief period of time while pretreating the infarction with SARP to restore homeostasis via a washout mechanism.

Mechanisms of coronary sinus reducer for treatment of myocardial ischemia: in silico study.

The coronary sinus reducer (CSR) is an emerging medical device for treating patients with refractory angina, often associated with myocardial ischemia. Patients implanted with CSR have shown positive outcomes, but the underlying mechanisms are unclear. This study sought to understand the mechanisms of CSR by investigating its effects on coronary microcirculation hemodynamics that may help explain the therapy's efficacy. We applied a validated computer model of the coronary microcirculation to investigate how CSR affects hemodynamics under different degrees of coronary artery stenosis. With moderate coronary stenosis, an increase in capillary transit time (CTT) [up to 69% with near-complete coronary sinus (CS) occlusion] is the key change associated with CSR. Because capillaries in the microcirculation can still receive oxygenated blood from the upstream artery with moderate stenosis, the increase in CTT allows more time for the exchange of gases and nutrients, aiding tissue oxygenation. With severe coronary stenosis; however, the redistribution of blood draining from the nonischemic region to the ischemic region (up to 96% with near-complete CS occlusion) and the reduction in capillary flow heterogeneity are the key changes associated with CSR. Because blood draining from the nonischemic region is not completely devoid of O2, the redistribution of blood to the capillaries in the ischemic region by CSR is beneficial especially when little or no oxygenated blood reaches these capillaries. This simulation study provides insights into the mechanisms of CSR in improving clinical symptoms. The mechanisms differ with the severity of the upstream stenosis.NEW & NOTEWORTHY Emerging coronary venous retroperfusion treatments, particularly coronary sinus reducer (CSR) for refractory angina linked to myocardial ischemia, show promise; however, their mechanisms of action are not well understood. We find that CSR's effectiveness varies with the severity of coronary stenosis. In moderate stenosis, CSR improves tissue oxygenation by increasing capillary transit time, whereas in severe stenosis, it redistributes blood from nonischemic to ischemic regions and reduces capillary flow heterogeneity.

Simulation of coronary capillary transit time based on full vascular model of the heart.

Capillary transit time (CTT) is a fundamental determinant of gas exchange between blood and tissues in the heart and other organs. Despite advances in experimental techniques, it remains difficult to measure coronary CTT in vivo. Here, we developed a novel computational framework that couples coronary microcirculation with cardiac mechanics in a closed-loop system that enables prediction of hemodynamics in the entire coronary network, including arteries, veins, and capillaries. We also developed a novel "particle-tracking" approach for computing CTT where "virtual tracers" are individually tracked as they traverse the capillary network. Model predictions compare well with blood pressure and flow rate distributions in the arterial network reported in previous studies. Model predictions of transit times in the capillaries (1.21 ± 1.5 s) and entire coronary network (11.8 ± 1.8 s) also agree with measurements. We show that, with increasing coronary artery stenosis (as quantified by fractional flow reserve, FFR), intravascular pressure and flow rate downstream are reduced but remain non-stationary even at 100 % stenosis because some flow (∼3 %) is redistributed from the non-occluded to the occluded territories. Importantly, the model predicts that occlusion of a large artery results in higher CTT. For moderate stenosis (FFR > 0.6), the increase in CTT (from 1.21 s without stenosis to 2.23 s at FFR=0.6) is caused by a decrease in capillary flow rate. In severe stenosis (FFR = 0.1), the increase in CTT to 14.2 s is due to both a decrease in flow rate and an increase in path length taken by "virtual tracers" in the capillary network.

Pre-arterialization of coronary veins prior to retroperfusion of ischemic myocardium: percutaneous closure device.

Background: Chronic coronary retroperfusion to treat myocardial ischemia has previously failed due to edema and hemorrhage of coronary veins suddenly exposed to arterial pressures. The objective of this study was to selectively adapt the coronary veins to become arterialized prior to coronary venous retroperfusion to avoid vascular edema and hemorrhage. Methods and results: In 32 animals (Group I = 19 and Group II = 13), the left anterior descending (LAD) artery was occluded using an ameroid occlusion model. In Group I, the great cardiac vein was blocked with suture ligation (Group IA = 11) or with occlusion device (Group IB = 8) to arterialize the venous system within 2 weeks at intermediate pressure (between arterial and venous levels) before a coronary venous bypass graft (CVBG) was implemented through a left internal mammary artery (LIMA) anastomosis. Group II only received the LAD artery occlusion and served as control. Serial echocardiograms showed recovery of left ventricular (LV) function with this adaptation-arterialization approach, with an increase in ejection fraction (EF) in Group I from 38% ± 5% after coronary occlusion to 53% ± 7% eight weeks after CVBG, whereas in Group II the EF never recovered (41% ± 2%-33% ± 7%). The remodeling of the venous system not only allowed restoration of myocardial function when CVBG was implemented but possibly promoted a novel form of "collateralization" between the native arterioles and the newly arterialized venules, which revascularized the ischemic myocardium. Conclusions: These findings form a potential rationale for a venous arterialization-revascularization treatment for the refractory angina and the "no-option" patients using a hybrid percutaneous (closure device for arterialization)/surgical approach (CVBG) to revascularize the myocardium.

An in silico study of the effects of left ventricular assist device on right ventricular function and inter-ventricular interaction.

BACKGROUND: Left ventricular assist device (LVAD) is associated with a high incidence of right ventricular (RV) failure, which is hypothesized to be caused by the occurring inter-ventricular interactions when the LV is unloaded. Factors contributing to these interactions are unknown. METHODS: We used computer modeling to investigate the impact of the HeartMate 3 LVAD on RV functions. The model was first calibrated against pressure-volume (PV) loops associated with a heart failure (HF) patient and validated against measurements of inter-ventricular interactions in animal experiments. The model was then applied to investigate the effects of LVAD on (1) RV chamber contractility indexed by V 60 derived from its end-systolic PV relationship, and (2) RV diastolic function indexed by V 20 derived from its end-diastolic PV relationship. We also investigated how septal wall thickness and regional contractility affect the impact of LVAD on RV function. RESULTS: The impact of LVAD on RV chamber contractility is small at a pump speed lower than 4k rpm. At a higher pump speed between 4k and 9k rpm, however, RV chamber contractility is reduced (by ~3% at 6k rpm and ~10% at 9k rpm). The reduction of RV chamber contractility is greater with a thinner septal wall or with a lower myocardial contractility at the LV free wall, septum, or RV free wall. CONCLUSION: RV chamber contractility is reduced at a pump speed higher than 4k rpm, and this reduction is greater with a thinner septal wall or lower regional myocardial contractility. Findings here may have clinical implications in identifying LVAD patients who may suffer from RV failure.

Mechanism of exercise intolerance in heart diseases predicted by a computer model of myocardial demand-supply feedback system.

BACKGROUND AND OBJECTIVE: The myocardial demand-supply feedback system plays an important role in augmenting blood supply in response to exercise-induced increased myocardial demand. During this feedback process, the myocardium and coronary blood flow interact bidirectionally at many different levels. METHODS: To investigate these interactions, a novel computational framework that considers the closed myocardial demand-supply feedback system was developed. In the framework coupling the systemic circulation of the left ventricle and coronary perfusion with regulation, myocardial work affects coronary perfusion via flow regulation mechanisms (e.g., metabolic regulation) and myocardial-vessel interactions, whereas coronary perfusion affects myocardial contractility in a closed feedback system. The framework was calibrated based on the measurements from healthy subjects under graded exercise conditions, and then was applied to simulate the effects of graded exercise on myocardial demand-supply under different physiological and pathological conditions. RESULTS: We found that the framework can recapitulate key features found during exercise in clinical and animal studies. We showed that myocardial blood flow is increased but maximum hyperemia is reduced during exercise, which led to a reduction in coronary flow reserve. For coronary stenosis and myocardial inefficiency, the model predicts that an increase in heart rate is necessary to maintain the baseline cardiac output. Correspondingly, the resting coronary flow reserve is exhausted and the range of heart rate before exhaustion of coronary flow reserve is reduced. In the presence of metabolic regulation dysfunction, the model predicts that the metabolic vasodilator signal is higher at rest, saturates faster during exercise, and as a result, causes quicker exhaustion of coronary flow reserve. CONCLUSIONS: Model predictions showed that the coronary flow reserve deteriorates faster during graded exercise, which in turn, suggests a decrease in exercise tolerance for patients with stenosis, myocardial inefficiency and metabolic flow regulation dysfunction. The findings in this study may have clinical implications in diagnosing cardiovascular diseases.

Optimization of cardiac resynchronization therapy based on a cardiac electromechanics-perfusion computational model.

Cardiac resynchronization therapy (CRT) is an established treatment for left bundle branch block (LBBB) resulting in mechanical dyssynchrony. Approximately 1/3 of patients with CRT, however, are non-responders. To understand factors affecting CRT response, an electromechanics-perfusion computational model based on animal-specific left ventricular (LV) geometry and coronary vascular networks located in the septum and LV free wall is developed. The model considers contractility-flow and preload-activation time relationships, and is calibrated to simultaneously match the experimental measurements in terms of the LV pressure, volume waveforms and total coronary flow in the left anterior descending and left circumflex territories from 2 swine models under right atrium and right ventricular pacing. The model is then applied to investigate the responses of CRT indexed by peak LV pressure and (dP/dt)max at multiple pacing sites with different degrees of perfusion in the LV free wall. Without the presence of ischemia, the model predicts that basal-lateral endocardial region is the optimal pacing site that can best improve (dP/dt)max by 20%, and is associated with the shortest activation time. In the presence of ischemia, a non-ischemic region becomes the optimal pacing site when coronary flow in the ischemic region fell below 30% of its original value. Pacing at the ischemic region produces little response at that perfusion level. The optimal pacing site is associated with one that optimizes the LV activation time. These findings suggest that CRT response is affected by both pacing site and coronary perfusion, which may have clinical implication in improving CRT responder rates.

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.

Role of coronary flow regulation and cardiac-coronary coupling in mechanical dyssynchrony associated with right ventricular pacing.

Mechanical dyssynchrony (MD) affects left ventricular (LV) mechanics and coronary perfusion. To understand the multifactorial effects of MD, we developed a computational model that bidirectionally couples the systemic circulation with the LV and coronary perfusion with flow regulation. In the model, coronary flow in the left anterior descending (LAD) and left circumflex (LCX) arteries affects the corresponding regional contractility based on a prescribed linear LV contractility-coronary flow relationship. The model is calibrated with experimental measurements of LV pressure and volume, as well as LAD and LCX flow rate waveforms acquired under regulated and fully dilated conditions from a swine under right atrial (RA) pacing. The calibrated model is applied to simulate MD. The model can simultaneously reproduce the reduction in mean LV pressure (39.3%), regulated flow (LAD: 7.9%; LCX: 1.9%), LAD passive flow (21.6%), and increase in LCX passive flow (15.9%). These changes are associated with right ventricular pacing compared with RA pacing measured in the same swine only when LV contractility is affected by flow alterations with a slope of 1.4 mmHg/mL2 in a contractility-flow relationship. In sensitivity analyses, the model predicts that coronary flow reserve (CFR) decreases and increases in the LAD and LCX with increasing delay in LV free wall contraction. These findings suggest that asynchronous activation associated with MD impacts 1) the loading conditions that further affect the coronary flow, which may explain some of the changes in CFR, and 2) the coronary flow that reduces global contractility, which contributes to the reduction in LV pressure.NEW & NOTEWORTHY A computational model that couples the systemic circulation of the left ventricular (LV) and coronary perfusion with flow regulation is developed to study the effects of mechanical dyssynchrony. The delayed contraction in the LV free wall with respect to the septum has a significant effect on LV function and coronary flow reserve.

Application of feed forward and recurrent neural networks in simulation of left ventricular mechanics.

An understanding of left ventricle (LV) mechanics is fundamental for designing better preventive, diagnostic, and treatment strategies for improved heart function. Because of the costs of clinical and experimental studies to treat and understand heart function, respectively, in-silico models play an important role. Finite element (FE) models, which have been used to create in-silico LV models for different cardiac health and disease conditions, as well as cardiac device design, are time-consuming and require powerful computational resources, which limits their use when real-time results are needed. As an alternative, we sought to use deep learning (DL) for LV in-silico modeling. We used 80 four-chamber heart FE models for feed forward, as well as recurrent neural network (RNN) with long short-term memory (LSTM) models for LV pressure and volume. We used 120 LV-only FE models for training LV stress predictions. The active material properties of the myocardium and time were features for the LV pressure and volume training, and passive material properties and element centroid coordinates were features of the LV stress prediction models. For six test FE models, the DL error for LV volume was 1.599 ± 1.227 ml, and the error for pressure was 1.257 ± 0.488 mmHg; for 20 LV FE test examples, the mean absolute errors were, respectively, 0.179 ± 0.050 for myofiber, 0.049 ± 0.017 for cross-fiber, and 0.039 ± 0.011 kPa for shear stress. After training, the DL runtime was in the order of seconds whereas equivalent FE runtime was in the order of several hours (pressure and volume) or 20 min (stress). We conclude that using DL, LV in-silico simulations can be provided for applications requiring real-time results.

Effects of Mechanical Dyssynchrony on Coronary Flow: Insights From a Computational Model of Coupled Coronary Perfusion With Systemic Circulation.

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 dP LV /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.

Intra-myocardial alginate hydrogel injection acts as a left ventricular mid-wall constraint in swine.

Despite positive initial outcomes emerging from preclinical and early clinical investigation of alginate hydrogel injection therapy as a treatment for heart failure, the lack of knowledge about the mechanism of action remains a major shortcoming that limits the efficacy of treatment design. To identify the mechanism of action, we examined previously unobtainable measurements of cardiac function from in vivo, ex vivo, and in silico states of clinically relevant heart failure (HF) in large animals. High-resolution ex vivo magnetic resonance imaging and histological data were used along with state-of-the-art subject-specific computational model simulations. Ex vivo data were incorporated in detailed geometric computational models for swine hearts in health (n = 5), ischemic HF (n = 5), and ischemic HF treated with alginate hydrogel injection therapy (n = 5). Hydrogel injection therapy mitigated elongation of sarcomere lengths (1.68 ± 0.10µm [treated] vs. 1.78 ± 0.15µm [untreated], p<0.001). Systolic contractility in treated animals improved substantially (ejection fraction = 43.9 ± 2.8% [treated] vs. 34.7 ± 2.7% [untreated], p<0.01). The in silico models realistically simulated in vivo function with >99% accuracy and predicted small myofiber strain in the vicinity of the solidified hydrogel that was sustained for up to 13 mm away from the implant. These findings suggest that the solidified alginate hydrogel material acts as an LV mid-wall constraint that significantly reduces adverse LV remodeling compared to untreated HF controls without causing negative secondary outcomes to cardiac function. STATEMENT OF SIGNIFICANCE: Heart failure is considered a growing epidemic and hence an important health problem in the US and worldwide. Its high prevalence (5.8 million and 23 million, respectively) is expected to increase by 25% in the US alone by 2030. Heart failure is associated with high morbidity and mortality, has a 5-year mortality rate of 50%, and contributes considerably to the overall cost of health care ($53.1 billion in the US by 2030). Despite positive initial outcomes emerging from preclinical and early clinical investigation of alginate hydrogel injection therapy as a treatment for heart failure, the lack of knowledge concerning the mechanism of action remains a major shortcoming that limits the efficacy of treatment design. To understand the mechanism of action, we combined high-resolution ex vivo magnetic resonance imaging and histological data in swine with state-of-the-art subject-specific computational model simulations. The in silico models realistically simulated in vivo function with >99% accuracy and predicted small myofiber strain in the vicinity of the solidified hydrogel that was sustained for up to 13 mm away from the implant. These findings suggest that the solidified alginate hydrogel material acts as a left ventricular mid-wall constraint that significantly reduces adverse LV remodeling compared to untreated heart failure controls without causing negative secondary outcomes to cardiac function. Moreover, if the hydrogel can be delivered percutaneously rather than via the currently used open-chest procedure, this therapy may become routine for heart failure treatment. A minimally invasive procedure would be in the best interest of this patient population; i.e., one that cannot tolerate general anesthesia and surgery, and it would be significantly more cost-effective than surgery.

Selective Autoretroperfusion Provides Substantial Cardioprotection in Swine: Incremental Improvements With Mild Hypothermia.

Mild hypothermia (MH) and retroperfusion are 2 techniques proposed to reduce infarct size due to myocardial infarction. The authors evaluated the effects of focal MH combined with selective coronary venous autoretroperfusion (SARP) as an acute cardioprotective modality before percutaneous coronary intervention (PCI) in a swine model of left ventricular myocardial infarction. Significant reduction in infarct size with preservation of cardiac function and cardiomyocyte viability were achieved. The authors propose that SARP alone or in combination with MH may provide a clinically relevant percutaneous short-term option of cardiac support to high-risk patients undergoing PCI.

Method for Calibration of Left Ventricle Material Properties Using Three-Dimensional Echocardiography Endocardial Strains.

We sought to calibrate mechanical properties of left ventricle (LV) based on three-dimensional (3D) speckle tracking echocardiographic imaging data recorded from 16 segments defined by American Heart Association (AHA). The in vivo data were used to create finite element (FE) LV and biventricular (BV) models. The orientation of the fibers in the LV model was rule based, but diffusion tensor magnetic resonance imaging (MRI) data were used for the fiber directions in the BV model. A nonlinear fiber-reinforced constitutive equation was used to describe the passive behavior of the myocardium, whereas the active tension was described by a model based on tissue contraction (Tmax). isight was used for optimization, which used abaqus as the forward solver (Simulia, Providence, RI). The calibration of passive properties based on the end diastolic pressure volume relation (EDPVR) curve resulted in relatively good agreement (mean error = -0.04 ml). The difference between the experimental and computational strains decreased after segmental strain metrics, rather than global metrics, were used for calibration: for the LV model, the mean difference reduced from 0.129 to 0.046 (circumferential) and from 0.076 to 0.059 (longitudinal); for the BV model, the mean difference nearly did not change in the circumferential direction (0.061) but reduced in the longitudinal direction from 0.076 to 0.055. The calibration of mechanical properties for myocardium can be improved using segmental strain metrics. The importance of realistic fiber orientation and geometry for modeling of the LV was shown.

Prediction of Left Ventricular Mechanics Using Machine Learning.

The goal of this paper was to provide a real-time left ventricular (LV) mechanics simulator using machine learning (ML). Finite element (FE) simulations were conducted for the LV with different material properties to obtain a training set. A hyperelastic fiber-reinforced material model was used to describe the passive behavior of the myocardium during diastole. The active behavior of the heart resulting from myofiber contractions was added to the passive tissue during systole. The active and passive properties govern the LV constitutive equation. These mechanical properties were altered using optimal Latin hypercube design of experiments to obtain training FE models with varied active properties (volume and pressure predictions) and varied passive properties (stress predictions). For prediction of LV pressures, we used eXtreme Gradient Boosting (XGboost) and Cubist, and XGBoost was used for predictions of LV pressures, volumes as well as LV stresses. The LV pressure and volume results obtained from ML were similar to FE computations. The ML results could capture the shape of LV pressure as well as LV pressure-volume loops. The results predicted by Cubist were smoother than those from XGBoost. The mean absolute errors were as follows: XGBoost volume: 1.734 ± 0.584 ml, XGBoost pressure: 1.544 ± 0.298 mmHg, Cubist volume: 1.495 ± 0.260 ml, Cubist pressure: 1.623 ± 0.191 mmHg, myofiber stress: 0.334 ± 0.228 kPa, cross myofiber stress: 0.075 ± 0.024 kPa, and shear stress: 0.050 ± 0.032 kPa. The simulation results show ML can predict LV mechanics much faster than the FE method. The ML model can be used as a tool to predict LV behavior. Training of our ML model based on a large group of subjects can improve its predictability for real world applications.

Intramyocardial Injections to De-Stiffen the Heart: A Subject-Specific in Silico Approach.

We hypothesized that minimally invasive injections of a softening agent at strategic locations in stiff myocardium could de-stiffen the left ventricle (LV) globally. Physics-based finite element models of the LV were created from LV echocardiography images and pressures recorded during experiments in four swine. Results confirmed animal models of LV softening by systemic agents. Regional de-stiffening of myocardium led to global de-stiffening of LV. The mathematical set up was used to design LV global de-stiffening by regional softening of myocardium. At an end diastolic pressure of 23 mmHg, when 8 ml of the free wall was covered by intramyocardial injections, end diastolic volume (EDV) increased by 15.0%, whereas an increase up to 11 ml due to intramyocardial injections in the septum and free wall led to a 26.0% increase in EDV. Although the endocardial intramyocardial injections occupied a lower LV wall volume, they led to an EDV (44 ml) that was equal compared to intramyocardial injections in the mid-wall (44 ml) and larger compared to intramyocardial injections in the epicardium (41 ml). Using an in silico set up, sites of regional myocardium de-stiffening could be planned in order to globally soften overly stiff LV in heart failure with preserved ejection fraction. This novel treatment is built on subject-specific data. Hypothesis-testing of these simulation findings in animal models is warranted.

Construction and Validation of Subject-Specific Biventricular Finite-Element Models of Healthy and Failing Swine Hearts From High-Resolution DT-MRI.

Predictive computational modeling has revolutionized classical engineering disciplines and is in the process of transforming cardiovascular research. This is particularly relevant for investigating emergent therapies for heart failure, which remains a leading cause of death globally. The creation of subject-specific biventricular computational cardiac models has been a long-term endeavor within the biomedical engineering community. Using high resolution (0.3 × 0.3 × 0.8 mm) ex vivo data, we constructed a precise fully subject-specific biventricular finite-element model of healthy and failing swine hearts. Each model includes fully subject-specific geometries, myofiber architecture and, in the case of the failing heart, fibrotic tissue distribution. Passive and active material properties are prescribed using hyperelastic strain energy functions that define a nearly incompressible, orthotropic material capable of contractile function. These materials were calibrated using a sophisticated multistep approach to match orthotropic tri-axial shear data as well as subject-specific hemodynamic ventricular targets for pressure and volume to ensure realistic cardiac function. Each mechanically beating heart is coupled with a lumped-parameter representation of the circulatory system, allowing for a closed-loop definition of cardiovascular flow. The circulatory model incorporates unidirectional fluid exchanges driven by pressure gradients of the model, which in turn are driven by the mechanically beating heart. This creates a computationally meaningful representation of the dynamic beating of the heart coupled with the circulatory system. Each model was calibrated using subject-specific experimental data and compared with independent in vivo strain data obtained from echocardiography. Our methods produced highly detailed representations of swine hearts that function mechanically in a remarkably similar manner to the in vivo subject-specific strains on a global and regional comparison. The degree of subject-specificity included in the models represents a milestone for modeling efforts that captures realism of the whole heart. This study establishes a foundation for future computational studies that can apply these validated methods to advance cardiac mechanics research.

Efficacy of intramyocardial injection of Algisyl-LVR for the treatment of ischemic heart failure in swine.

BACKGROUND: Progressive thinning and dilation of the LV due to ischemic heart failure (IHF) increases wall stress and myocardial oxygen consumption. Injectable biopolymers implanted in the myocardial wall have been used to increase wall thickness to reduce chamber volume, decrease wall stress, and improve cardiac function. We sought to evaluate the efficacy of a biopolymer (Algisyl-LVR) to prevent left ventricular (LV) remodeling in a swine model of IHF. METHODS: IHF was induced in 11 swine by occluding the marginal obtuse branches of the left circumflex artery. Eight weeks later, Algisyl-LVR was injected into the LV myocardial free wall in five of the 11 animals. Echocardiographic examinations were done every 2weeks for 16weeks. RESULTS: Within eight weeks of treatment, the ejection fraction increased from 30.5%±7.7% to 42.4%±3.5% (treated group) vs. 37.3%±3.8% to 34.3%±2.9% (control), p<0.01. Stroke volume increased from 18.5±9.3mL to 41.3±13.3mL (treated group) vs. 25.4±2.3mL to 31.4±5.3mL (control), p<0.05. Wall thickness in end-diastole of the infarcted region changed from 0.69±0.06cm to 0.81±0.13cm (treated group) vs. 0.73±0.09cm to 0.68±0.11cm (control), p<0.05. Sphericity index remained almost unchanged after treatment, although differences were found at the end of the study between both groups (p<0.001). Average myofiber stress changed from 16.3±5.8kPa to 10.2±4.0kPa (treated group) vs. 15.2±4.8kPa to 17.9±5.6kPa (control), p<0.05. CONCLUSIONS: Algisyl-LVR is an effective strategy that serves as a micro-LV assist device to reduce stress and hence prevent or reverse maladaptive cardiac remodeling caused by IHF in swine.

Compensatory Enlargement of Ossabaw Miniature Swine Coronary Arteries in Diffuse Atherosclerosis.

Studies in human and non-human primates have confirmed the compensatory enlargement or positive remodeling (Glagov phenomenon) of coronary vessels in the presence of focal stenosis. To our knowledge, this is the first study to document arterial enlargement in a metabolic syndrome animal model with diffuse coronary artery disease (DCAD) in the absence of severe focal stenosis. Two different groups of Ossabaw miniature pigs were fed a high fat atherogenic diet for 4 months (Group I) and 12 months (Group II), respectively. Group I (6 pigs) underwent contrast enhanced computed tomographic angiography (CCTA) and intravascular ultrasound (IVUS) at baseline and after 4 months of high fat diet, whereas Group II (7 pigs) underwent only IVUS at 12 months of high fat diet. IVUS measurements of the left anterior descending (LAD), left circumflex (LCX) and right coronary (RCA) arteries in Group I showed an average increase in their lumen cross-sectional areas (CSA) of 25.8%, 11.4%, and 43.4%, respectively, as compared to baseline. The lumen CSA values of LAD in Group II were found to be between the baseline and 4 months values in Group I. IVUS and CCTA measurements showed a similar trend and positive correlation. Fractional flow reserve (FFR) was 0.91±0.07 at baseline and 0.93±0.05 at 4 months with only 2.2%, 1.6% and 1% stenosis in the LAD, LCX and RCA, respectively. The relation between percent stenosis and lumen CSA shows a classical Glagov phenomenon in this animal model of DCAD.

Longitudinal hemodynamic measurements in swine heart failure using a fully implantable telemetry system.

Chronic monitoring of heart rate, blood pressure, and flow in conscious free-roaming large animals can offer considerable opportunity to understand the progression of cardiovascular diseases and can test new diagnostics and therapeutics. The objective of this study was to demonstrate the feasibility of chronic, simultaneous measurement of several hemodynamic parameters (left ventricular pressure, systemic pressure, blood flow velocity, and heart rate) using a totally implantable multichannel telemetry system in swine heart failure models. Two solid-state blood pressure sensors were inserted in the left ventricle and the descending aorta for pressure measurements. Two Doppler probes were placed around the left anterior descending (LAD) and the brachiocephalic arteries for blood flow velocity measurements. Electrocardiographic (ECG) electrodes were attached to the surface of the left ventricle to monitor heart rate. The telemeter body was implanted in the right side of the abdomen under the skin for approximately 4 to 6 weeks. The animals were subjected to various heart failure models, including volume overload (A-V fistula, n = 3), pressure overload (aortic banding, n = 2) and dilated cardiomyopathy (pacing-induced tachycardia, n = 3). Longitudinal changes in hemodynamics were monitored during the progression of the disease. In the pacing-induced tachycardia animals, the systemic blood pressure progressively decreased within the first 2 weeks and returned to baseline levels thereafter. In the aortic banding animals, the pressure progressively increased during the development of the disease. The pressure in the A-V fistula animals only showed a small increase during the first week and remained stable thereafter. The results demonstrated the ability of this telemetry system of long-term, simultaneous monitoring of blood flow, pressure and heart rate in heart failure models, which may offer significant utility for understanding cardiovascular disease progression and treatment.