Development and applications of in-vitro and in-silico models of the cardiovascular system to study the effects of mechanical circulatory support.

Cardiovascular diseases (CVDs) are the leading cause of mortality globally. With ongoing interest in CVDs treatment, preclinical models for drug/therapeutic development that allow for fast iterative research are needed. Owing to the inherent complexity of the cardiovascular system, current in-vitro models of the cardiovascular system fail to replicate many of the physiological aspects of the cardiovascular system. In this dissertation, the main concern is with heart failure (HF). In advanced HF, patients may receive Left Ventricular Assist Devices (LVADs) as a bridge to transplant or destination therapy. However, LVADs have many limitations, including inability to adapt to varying tissue demand conditions, risk of ventricular suction, and diminished arterial pulsatility. To address these issues, this dissertation aims to use and develop computer, cellular, and tissue models of the cardiovascular system. 1) Use an in-silico model of the cardiovascular system to develop a novel control algorithm for LVADs. The control system was rigorously tested and showed adequate perfusion during rest and exercise, protect against ventricular suction under reduced heart preload, and augment arterial pulsatility through pulse modulation without requiring sensor implantation or model-based estimations. 2) While pulsatility augmentation was feasible through the developed control algorithm, the pulse waveform that could normalize the vascular phenotype is unknown. To address this, an endothelial cell-smooth muscle cell microfluidic coculture model was developed to recreate the physiological mechanical stimulants in the vascular wall. The results demonstrated different effects of pulsatile shear stress and stretch on endothelial cells and may indicate that a pulse pressure of at least 30 mmHg is needed to maintain normal endothelial morphology. 3) In order to study the effects of mechanical unloading on the native ventricle, a novel cardiac tissue culture model (CTCM) was developed. CTCM provided physiological electromechanical and humoral stimulation with 25% preload stretch and thyroid and glucocorticoid treatment maintained the cardiac phenotype for 12 days. The device was thoroughly characterized and tested. Results demonstrated improved viability, energy utilization, fibrotic remodeling, and structural integrity compared to available culture systems. The system was also used to reproduce ventricular volume-overload and the results demonstrated hypertrophic and fibrotic remodeling, typical of volume-overload pathology

From Benchtop to Beside: Patient-specific Outcomes Explained by Invitro Experiment

Study: Recent analyses show that females have higher early postoperative (PO) mortality and right ventricular failure (RVF) than males after left ventricular assist device (LVAD) implantation; and that this association is partially mediated by smaller LV size in females. Benchtop experiments allow us to investigate patient-specific (PS) characteristics in a reproducible way given the fact that the PS anatomy and physiology is mimicked accurately. With multiple heart models of varying LV size, we can directly study the individual effects of titrating the LVAD speed and the resulting bi-ventricular volumes, shedding light on the interplay between LV and RV as well as resulting inter-ventricular septum (IVS) positions, which may cause the different outcomes pertaining to sex. Methods: In vitro, we studied the impact of the heart size to IVS position using two smaller and two larger sized PS silicone heart phantoms derived from clinical CT images (Fig. 1A). With ultrasound crystals that were integrated on a placeholder inflow cannula, the IVS position was measured during LV and RV volume changes (dV) mimicking varying ventricular loading states (Fig. 1B). Figure 1 A Two small (blue) and two large PS heart phantoms (orange) on B benchtop. C Median septum curvature results. LVEDD/LVV/RVV: LV enddiastolic diameter/LV and RV volume. Results: Going from small to large dV, at zero curvature, the septum starts to shift towards the left; for smaller hearts at dV = -40 mL and for larger hearts at dV = -50 mL (Fig. 1C). This result indicates that smaller hearts are more prone to an IVS shift to the left than larger hearts. We conclude that smaller LV size may therefore mediate increased early PO LVAD mortality and RVF observed in females compared to males. Novel 3D silicone printing technology enables us to study accurate, PS heart models across a heterogeneous patient population. PS relationships can be studied simultaneously to clinical assessments and support the decision-making prior to LVAD implantation

Feasibility study of intelligent LVAD control for optimal heart failure therapy.

Background: Left ventricular assist devices (LVAD) are operated at constant speeds (rpm), consequently, pump flow is passively determined by the pressure difference between the LV and aorta. Since the diastolic pressure gradient (~70 mmHg) is much larger than the systolic gradient (~10 mmHg), the majority of pump flow occurs during systole. This limitation results in sub-optimal LV volume unloading, LV washing, and diminished vascular pulsatility that may be associated with increased risk for clinically-significant adverse events, including stroke, bleeding, arteriovenous malformations, and aortic insufficiency. To address these clinical adverse events, an intelligent control strategy using pump speed modulation was developed to provide dynamic LV unloading during the cardiac cycle to produce near-physiologic pulsatile flow delivery similar to that of the native heart. Materials and Methods: The objective of this study was to integrate a novel algorithm to dynamically control Medtronic HVAD pump speed and demonstrate proofof-concept by characterizing hemodynamic performance in a mock flow loop primed with a blood analog solution (glycerol-saline, 3 cP) and tuned to simulate class IV heart failure (HF). The intelligent LVAD control was operated a varying pump speeds (Dspeed = 0, 1000, 1500, 2000, 2500 rpm) and systolic durations (30%, 35%, and 40%); systolic duration correlates to the time spent at either the high or low pump speed setting. The intelligent LVAD control strategy modulates pump speed within a cardiac cycle triggered from an R-wave of an EKG waveform set to 80 BPM. This pump speed modulation control strategy allows for pulsatile operation of a continuous flow LVAD within a single cardiac cycle. Hemodynamic waveforms (LV pressure-volume, aortic pressure-flow, and pump flow) and intrinsic pump parameters (speed and current) were recorded and analyzed for each test condition. We hypothesize that pump speed modulation may be configured for optimal volume unloading (rest), vascular pulsatility (reloading), and/or washing. Results and Discussion: The intelligent LVAD control system successfully demonstrated the ability to rapidly increase and decrease HVAD pump speed within a single cardiac cycle to provide asynchronous, synchronous co-pulsation, and synchronous counter-pulsation profiles for all systolic durations (30, 35, 40%) and Drpm tested (D1000, D1500, D2000, D2500). Asynchronous support was achieved when pump speed increase (or decrease) was independent of the cardiac cycle, co-pulsation support was achieved when increase in pump speed was timed with beginning of systole corresponding with ventricular contraction (systole), and counter-pulsation support was when increase in pump speed was timed with the end of systole corresponding with ventricular filling (diastole). Ideally, the intelligent control would increase (or decrease) the HVAD pump speed instantaneously upon R-wave detection; however, two distinct time delays were observed: (1) a time delay from detection of the R-wave trigger and increase (or decrease) of pump speed for systolic durations of 35% and 40% (being 45 ± 3.0 ms and 82 ± 3.0 ms respectively and (2) a delay in LVAD flow when pump speed was increased which is hypothesized to be from the blood analog solution’s fluid inertia. Left ventricular stroke volume decreased for all LVAD pump speed modulation operating conditions compared to baseline (HF with LVAD off) indicating that the intelligent control strategy was able to reduce LV volume with increasing HVAD support. The highest flow was achieved with the HVAD operated at a fixed speed of 4000 rpm; however, co-pulsation pump speed modulation at the largest pump speed differential (low = 1500, high = 4000, Drpm = 2500, and systolic duration 30%) resulted in a mean pump speed 3,300 ± 1,200 rpm. By comparison, the forward flow at fixed pump speed of 4,000 rpm was 4.8 L/min compared to a mean co-pulsation rpm was 4.5 L/min. Additionally, all operating settings for the intelligent control during pulsatile function produced an average forward flow through the aortic valve, while in contrast at higher fixed speeds (3,500 and 4,000 rpm) the mean aortic flow was negative. Pulse pressure (DP) decreased with increasing mean pump speed (rpm) for all operating modes (fixed, asynchronous, co-pulsation, counter-pulsation). When operating at the same mean pump speed (rpm) copulsation has increased hemodynamic benefit for pulsatility when compared to counterpulsation and fixed speed at the same mean pump (rpm). Conclusion: The results of this study show the ability of the intelligent HVAD control strategy to increase and decrease pump speed within a single cardiac cycle. This study showed that asynchronous modulation with phases of co-pulsation can generate near physiologic pulse pressure and vascular pulsatility when compared to counterpulsation support, while counter-pulsation can generate greater ventricular volume unloading and diastolic augmentation when compared to co-pulsation. Furthermore, the clinical impact of this study is that through speed modulation adverse events of continuous flow LVADs may be reduced such as incidences of bleeding associated with decreased pulsatility and a decrease in the risk of thrombus formation from poor washing around the aortic valve

Preload-based Starling-like control of rotary blood pumps: an in-vitro evaluation

Due to a shortage of donor hearts, rotary left ventricular assist devices (LVADs) are used to provide mechanical circulatory support. To address the preload insensitivity of the constant speed controller (CSC) used in conventional LVADs, we developed a preload-based Starling-like controller (SLC). The SLC emulates the Starling law of the heart to maintain mean pump flow ([Formula: see text]) with respect to mean left ventricular end diastolic pressure (PLVEDm) as the feedback signal. The SLC and CSC were compared using a mock circulation loop to assess their capacity to increase cardiac output during mild exercise while avoiding ventricular suction (marked by a negative PLVEDm) and maintaining circulatory stability during blood loss and severe reductions in left ventricular contractility (LVC). The root mean squared hemodynamic deviation (RMSHD) metric was used to assess the clinical acceptability of each controller based on pre-defined hemodynamic limits. We also compared the in-silico results from our previously published paper with our in-vitro outcomes. In the exercise simulation, the SLC increased [Formula: see text] by 37%, compared to only 17% with the CSC. During blood loss, the SLC maintained a better safety margin against left ventricular suction with PLVEDm of 2.7 mmHg compared to -0.1 mmHg for CSC. A transition to reduced LVC resulted in decreased mean arterial pressure (MAP) and [Formula: see text] with CSC, whilst the SLC maintained MAP and [Formula: see text]. The results were associated with a much lower RMSHD value with SLC (70.3%) compared to CSC (225.5%), demonstrating improved capacity of the SLC to compensate for the varying cardiac demand during profound circulatory changes. In-vitro and in-silico results demonstrated similar trends to the simulated changes in patient state however the magnitude of hemodynamic changes were different, thus justifying the progression to in-vitro evaluation.Mahdi Mansouri, Shaun D. Gregory, Robert F. Salamonsen, Nigel H. Lovell, Michael C. Stevens, Jo P. Pauls, Rini Akmeliawati, Einly Li

Design and development of a pulsatile axial flow blood pump as a left ventricular assist device

This thesis was submitted for the award of Doctor of Philosophy and was awarded by Brunel University LondonEach year all over the world, Millions of patients from infants to adults are diagnosed with heart failure. A limited number of donor hearts available for these patients results in a tremendous demand of mechanical circulatory support (MCS) system, either in the form of total artificial heart (TAH) or a ventricular assist device (VAD). Physiologically MCS are expected to provide heart; a time to rest and potential recovery by unloading the ventricle, while maintaining the adequate peripheral as well as coronary circulation. Existing ventricular assist devices (VAD) have employed either displacement type pulsatile flow pumping systems or continuous flow type centrifugal/rotodynamic pumps systems. Displacement type devices produce a pulsatile outflow, which has significant benefits on vital organ function and end organ recovery. Continuous flow devices are small and can be placed within body using minimal invasive procedures, in addition they reduces infection as well as mechanical failure related complications. Despite availability of success stories for both types of pumping systems, the selection of the either of them is an ongoing debate. This thesis aims to merge the advantages of displacement pumps (pulsatile flow) and axial-flow pumps (continuous flow) into a novel left vertical assist device (LVAD), by designing a novel minimal invasive, miniature axial-flow pump producing pulsating outflow for the patients having early heart failure and myocardial infarction as a Bridge-To-Recovery (BTR) or Bridge-To-Decision (BTD) device. The design of VAD, the experimental setup and dedicated control system were developed for the in vitro evaluation of pulsatile flow. Computational fluid dynamics (CFD) had been employed for the detail investigation of pulsatile flow. In addition, CFD was also applied to optimize the pulse generation for low haemolysis levels. Outcome of the study produces comprehensive understanding for the generation of pulsatile flow using an axial flow pump. Further, it provides the means of generating a controlled pulse that can regulate flow rate for varying heart rate within low haemolysis levels

Towards patient-specific modelling as a pre-operative planning strategy and follow up assessment for the treatment of advanced heart failure with rotary blood pumps

Background: Ventricular Assist Devices (VADs) insertion is an established treatment for patients with end-stage heart failure waiting for a heart transplant or in need for long-term circulatory support (destination therapy). Rotary blood pumps (RBP) are the most popular devices in view of their size and performance. Pre-operative planning strategy for the insertion of a left ventricular assist device (LVAD) requires a timely discussion at a Multi-Disciplinary Team Meeting (MDT). Clinical-decision making is based according to the needs of the patient and must be processed without delays. Nevertheless, thrombus formation remains a feared complication which affects outcome. VADs operate in a flow regime which is difficult to simulate: the transitional region at the boundary of laminar and turbulent flow (low Reynolds number). Different methods have been used but the best approach remains debatable. Computational Fluid Dynamics (CFD) is an attractive and invaluable tool for the study of the interactions between VADs and the cardiovascular system. The aim of this thesis is three-fold: a) to investigate the use of pressure-volume analysis in a clinical setting through the review of six heart failure patients previously discussed at a MDT meeting with a view to predict or guide further management; b) to review the theory behind modelling approaches to VADs and their interactions with the cardiovascular system for better understanding of their clinical use. Then, an overview of computational fluid dynamics (CFD) is considered as a prelude to its application to the analysis of VADs performance. Additionally, the development of a simplified model of centrifugal pump will be used in initial simulations as preliminary analysis; c) to examine an example of a proof-of-concept pilot patient-specific model of an axial flow pump (HeartMate II) as pre-operative planning strategy in a patient-specific model with a view to identify potential critical areas that may affect pump function and outcome in a clinical setting. Material and Methods: 3D reconstruction from CT-scan images of patients who underwent the insertion of rotary blood pumps, namely HeartWare HVAD and HeartMate II. Ansys Fluent has been used for CFD analysis based on the fundamental governing equations of motion. Blood has been modelled as incompressible, Newtonian fluid with density = 1060 and viscosity = 0.0035 kg/m-s. The laminar and SST models have been used for comparison purposes. The rotational motion of the impeller has been implemented using the moving reference frame (MRF) approach. The sliding mesh method has also been used to account for unsteady interaction between stationary and moving part. The no-slip condition has been applied to all walls, which were assumed to be rigid. Boundary conditions consisting of velocity inlet and pressure outlet of the pump based on different settings and constant rotational speed for the impeller. Pressure-velocity coupling has been based on the coupled scheme. Spatial discretisation consisted of the “least square cell based” gradient for velocity and “PRESTO” or second order for pressure. Second order upwind has been set for the momentum, turbulent kinetic energy and specific dissipation rate. First order implicit has been set for transient formulation. The pseudo transient algorithm (steady state), the high order relaxation term and the warped-face gradient correction have been used to add an unsteady term to the solution equations with the aim to improve stability and enhance convergence. Specific settings have been considered for comparison purposes. Results: Pressure-volume simulation analysis in six advanced heart failure patients showed that an integrated model of the cardiovascular system based on lumped-parameter representation, modified time-varying elastance and pressure-volume analysis of ventricular function seems a feasible and suitable approach yielding a sufficiently accurate quantitative analysis in real time, therefore applicable within the time-constraints of a clinical setting. Lumped-parameter models consist of simultaneous ordinary differential equations complemented by an algebraic balance equation and are suitable for examination of global distribution of pressure, flow and volume over a range of physiological conditions with inclusion of the interaction between modelled components. Higher level lumped-parameter modelling is needed to address the interaction between the circulation and other systems based on a compromise between complexity and ability to set the required parameters to personalise an integrated lumped-parameter model for a patient-specific approach. CARDIOSIM© fulfils these requirements and does address the systems interaction with its modular approach and assembly of models with varying degree of complexity although 0-D and 1-D coupling may be required for the evaluation of long-term VAD support. The challenge remains the ability to predict outcome over a longer period of time. The preliminary CFD simulations with the HeartWare HVAD centrifugal pump demonstrated that it is possible to obtain an accurate analysis in a timely manner to complement the clinical review process. The simulations with the pilot patient-specific model of the HeartMate II axial flow pump revealed that a complex 3D reconstruction is feasible in a timely manner and can be used to generate sufficiently accurate results to be used in the context of a MDT meeting for the purposes of clinical decision-making. Overall, these three studies demonstrate that the time frame of the simulations was within hours which may fit the time constraints of the clinical environment in the context of a MDT meeting. More specifically, it was shown that the laminar model may be used for an initial evaluation of the flow development within the pump. Nonetheless, the k- model offers higher accuracy if the timeline of the clinical setting allows for a longer simulation. Conclusion: This thesis aimed at the understanding of the use of computational modelling as a pre-operative planning strategy and follow up assessment for the treatment of advanced heart failure with rotary blood pumps. The novelty lays in the use of both pressure-volume simulation analysis and 3D flow dynamics studies in VADs with a view to treatment optimisation and outcome prediction within the time constraints of a clinical setting in the context of a MDT meeting. The clinical significance and the contribution to the field is a more targeted approach for different groups of patients and a more quantitative evaluation in the clinical decision process based on a pro-active co-operation between clinicians and scientists reducing the potential for “guess work”. The results of this thesis are a proof-of-concept as a prelude to a potential future implementation of patient-specific modelling within a clinical setting on a daily basis demonstrating a clear clinical significance and contribution to the field. The proposed approach does not consider modelling and simulation as a substitute for clinical experience but an additional tool to guide therapeutic intervention and complement the clinical decision process in which the clinician remains the ultimate decision-maker. Such an approach may well add a different dimension to the problem of heart failure with potential for high return in terms of patient’s outcome and long-term surveillance. The same principles would be applicable to other cardiovascular problems in line with the current concept of “Team Approach” such as the Heart Team, the Structural Heart Team or the Aortic Team. The present work has taken this concept closer to clinical delivery and has highlighted its potential but further work remains to be done in refining the technique.Background: Ventricular Assist Devices (VADs) insertion is an established treatment for patients with end-stage heart failure waiting for a heart transplant or in need for long-term circulatory support (destination therapy). Rotary blood pumps (RBP) are the most popular devices in view of their size and performance. Pre-operative planning strategy for the insertion of a left ventricular assist device (LVAD) requires a timely discussion at a Multi-Disciplinary Team Meeting (MDT). Clinical-decision making is based according to the needs of the patient and must be processed without delays. Nevertheless, thrombus formation remains a feared complication which affects outcome. VADs operate in a flow regime which is difficult to simulate: the transitional region at the boundary of laminar and turbulent flow (low Reynolds number). Different methods have been used but the best approach remains debatable. Computational Fluid Dynamics (CFD) is an attractive and invaluable tool for the study of the interactions between VADs and the cardiovascular system. The aim of this thesis is three-fold: a) to investigate the use of pressure-volume analysis in a clinical setting through the review of six heart failure patients previously discussed at a MDT meeting with a view to predict or guide further management; b) to review the theory behind modelling approaches to VADs and their interactions with the cardiovascular system for better understanding of their clinical use. Then, an overview of computational fluid dynamics (CFD) is considered as a prelude to its application to the analysis of VADs performance. Additionally, the development of a simplified model of centrifugal pump will be used in initial simulations as preliminary analysis; c) to examine an example of a proof-of-concept pilot patient-specific model of an axial flow pump (HeartMate II) as pre-operative planning strategy in a patient-specific model with a view to identify potential critical areas that may affect pump function and outcome in a clinical setting. Material and Methods: 3D reconstruction from CT-scan images of patients who underwent the insertion of rotary blood pumps, namely HeartWare HVAD and HeartMate II. Ansys Fluent has been used for CFD analysis based on the fundamental governing equations of motion. Blood has been modelled as incompressible, Newtonian fluid with density = 1060 and viscosity = 0.0035 kg/m-s. The laminar and SST models have been used for comparison purposes. The rotational motion of the impeller has been implemented using the moving reference frame (MRF) approach. The sliding mesh method has also been used to account for unsteady interaction between stationary and moving part. The no-slip condition has been applied to all walls, which were assumed to be rigid. Boundary conditions consisting of velocity inlet and pressure outlet of the pump based on different settings and constant rotational speed for the impeller. Pressure-velocity coupling has been based on the coupled scheme. Spatial discretisation consisted of the “least square cell based” gradient for velocity and “PRESTO” or second order for pressure. Second order upwind has been set for the momentum, turbulent kinetic energy and specific dissipation rate. First order implicit has been set for transient formulation. The pseudo transient algorithm (steady state), the high order relaxation term and the warped-face gradient correction have been used to add an unsteady term to the solution equations with the aim to improve stability and enhance convergence. Specific settings have been considered for comparison purposes. Results: Pressure-volume simulation analysis in six advanced heart failure patients showed that an integrated model of the cardiovascular system based on lumped-parameter representation, modified time-varying elastance and pressure-volume analysis of ventricular function seems a feasible and suitable approach yielding a sufficiently accurate quantitative analysis in real time, therefore applicable within the time-constraints of a clinical setting. Lumped-parameter models consist of simultaneous ordinary differential equations complemented by an algebraic balance equation and are suitable for examination of global distribution of pressure, flow and volume over a range of physiological conditions with inclusion of the interaction between modelled components. Higher level lumped-parameter modelling is needed to address the interaction between the circulation and other systems based on a compromise between complexity and ability to set the required parameters to personalise an integrated lumped-parameter model for a patient-specific approach. CARDIOSIM© fulfils these requirements and does address the systems interaction with its modular approach and assembly of models with varying degree of complexity although 0-D and 1-D coupling may be required for the evaluation of long-term VAD support. The challenge remains the ability to predict outcome over a longer period of time. The preliminary CFD simulations with the HeartWare HVAD centrifugal pump demonstrated that it is possible to obtain an accurate analysis in a timely manner to complement the clinical review process. The simulations with the pilot patient-specific model of the HeartMate II axial flow pump revealed that a complex 3D reconstruction is feasible in a timely manner and can be used to generate sufficiently accurate results to be used in the context of a MDT meeting for the purposes of clinical decision-making. Overall, these three studies demonstrate that the time frame of the simulations was within hours which may fit the time constraints of the clinical environment in the context of a MDT meeting. More specifically, it was shown that the laminar model may be used for an initial evaluation of the flow development within the pump. Nonetheless, the k- model offers higher accuracy if the timeline of the clinical setting allows for a longer simulation. Conclusion: This thesis aimed at the understanding of the use of computational modelling as a pre-operative planning strategy and follow up assessment for the treatment of advanced heart failure with rotary blood pumps. The novelty lays in the use of both pressure-volume simulation analysis and 3D flow dynamics studies in VADs with a view to treatment optimisation and outcome prediction within the time constraints of a clinical setting in the context of a MDT meeting. The clinical significance and the contribution to the field is a more targeted approach for different groups of patients and a more quantitative evaluation in the clinical decision process based on a pro-active co-operation between clinicians and scientists reducing the potential for “guess work”. The results of this thesis are a proof-of-concept as a prelude to a potential future implementation of patient-specific modelling within a clinical setting on a daily basis demonstrating a clear clinical significance and contribution to the field. The proposed approach does not consider modelling and simulation as a substitute for clinical experience but an additional tool to guide therapeutic intervention and complement the clinical decision process in which the clinician remains the ultimate decision-maker. Such an approach may well add a different dimension to the problem of heart failure with potential for high return in terms of patient’s outcome and long-term surveillance. The same principles would be applicable to other cardiovascular problems in line with the current concept of “Team Approach” such as the Heart Team, the Structural Heart Team or the Aortic Team. The present work has taken this concept closer to clinical delivery and has highlighted its potential but further work remains to be done in refining the technique