ABSTRACT
Background and Objective The main indications for right ventricular assist device (RVAD) support are right heart failure after implantation of a left ventricular assist device (LVAD) or early graft failure following heart transplantation. About 30-40% of patients will need RVAD support after LVAD implantation. Pulmonary hypertension is also an indication for right heart assistance. Several types of RVAD generating pulsatile or continuous flow are available on the market. These assist devices can be connected to the cardiovascular system in different ways. We sought to analyse the effects induced by different RVAD connections when right ventricular elastance is modified using a numerical simulator. The analysis was based on the behaviour of both left and right ventricular and atrial loops in the pressure-volume plane.
Methods New modules of the cardiovascular network and a right ventricular centrifugal pump were implemented in CARDIOSIM© software simulator platform. The numerical pump model generated continuous flow when connected in series or parallel to the right ventricle. When the RVAD was connected in series (parallel), the pump removed blood from the right ventricle (atrium) and ejected it into the pulmonary artery. In our study, we analysed the effects induced by RVAD support on left/right ventricular/atrial loops when right ventricular elastance slope (EesRIGHT) changed from 0.3 to 0.8 mmHg/ml with the pump connected either in series or parallel. The effect of low and high rotational pump speed was also addressed.
Results Percentage changes up to 60% were observed for left ventricular pressure-volume area and external work during in-parallel RVAD support at 4000 rpm with EesRIGHT = 0.3mmHg/ml. The same pump setting and connection type led to percentage variation up to 20% for left ventricular ESV and up to 25% for left ventricular EDV with EesRIGHT = 0.3mmHg/ml. Again the same pump setting and connection generated up to 50% change in left atrial pressure-volume loop area (PVLAL-A) and only 3% change in right atrial pressure-volume loop area (PVLAR-A) when EesRIGHT = 0.3mmHg/ml. Percentage variation was lower when EesRIGHT was increased up to 0.8 mmHg/ml.
Conclusion Early recognition of right ventricular failure followed by aggressive treatment is desirable to achieve a more favourable outcome. RVAD support remains an option for advanced right ventricular failure although onset of major adverse events may preclude its use.
INTRODUCTION
The main indications for right ventricular assist device (RVAD) support are right heart failure after implantation of a left ventricular assist device (LVAD) or early graft failure following heart transplantation. Measurement of right heart function is often overlooked during daily clinical practice. The study by [1] is an attempt to better understand RV pathophysiology in the setting of pulmonary hypertension, which represents the prevalent pathological condition for RVAD indication, whether to support a failing left ventricle or a primarily involved RV [2].
However, the lack of data does not allow an in-depth analysis of right ventricular and atrial behaviour during RVAD support. An attempt in this direction can be made using a numerical simulator of the cardiovascular system with a view to reproduce pathological conditions requiring right ventricular assistance.
The aim of our study was the analysis of the haemodynamic and energetic variables of both ventricles and atria in different cardiovascular conditions reproduced by changing right ventricular contractility from 0.3 to 0.8 mmHg/ml during RVAD support with different rotational speeds.
The first step of this work required the implementation of two new modules within CARDIOSIM© software simulation platform, which would reproduce the characteristics of the left circulatory network and a continuous flow centrifugal pump (RVAD) to be connected either in series or parallel to the right ventricle. The new modules were based on a 0-D (lumped-parameter) numerical model including input and output cannula of the RVAD.
We simulated RVAD support following in series and parallel connection driven by different rotational speeds in a heart failure setting. Subsequently, the right ventricular elastance was increased from 0.3 to 0.8 mmHg/ml in a stepwise manner. During each setting, we focused our attention on the following hemodynamic and energetic variables:
✓ Right and left ventricular end-systolic volume (ESVR-V and ESVL-V);
✓ Right and left ventricular end-diastolic volume (EDVR-V and EDVL-V);
✓ Stroke volume (SV);
✓ Right and left ventricular external work (EWR-V and EWL-V) and pressure-volume area (PVAR-V and PVAL-V);
✓ Right and left atrial end-systolic volume (ESVR-A and ESVL-A);
✓ Right and left atrial end-diastolic volume (EDVR-A and EDVL-A);
✓ Right and left atrial pressure-volume loop area (PVLAR-A and PVLAL-A);
✓ Cardiac output (CO);
✓ Systolic, diastolic and mean systemic aortic pressure (AoP);
✓ Systolic, diastolic and mean pulmonary arterial pressure (PAP);
✓ Pulmonary capillary wedge pressure (PCWP);
✓ Right and left atrial pressure (RAP and LAP).
MATERIAL AND METHODS
The heart and circulatory numerical network
CARDIOSIM© software platform has been previously described [3]-[8]. The modules of the software simulator are: systemic and pulmonary arterial section, systemic and pulmonary venous section and the coronary circulation. Native left and right ventricles, atria and septum reproduce the entire cardiac activity; they are implemented in a single module (Fig. 1) using the time-varying elastance concept [4],[5]. The ventricular, atrial and septal activity is synchronized with the electrocardiographic (ECG) signal [5]. Using the modified time-varying elastance theory, the inter-ventricular septum (IVS) interaction and the instantaneous left and right ventricular pressure can be reproduced by [5],[10]:
In the same way, the inter-atrial septum (IAS) interaction and the left and right instantaneous atrial pressure can be reproduced by: The symbols used in Eqs. (1) and (2) have been listed in Table 1. These features allow the simulation of inter-ventricular and intra-ventricular dyssynchrony [10].
Specific modules of the coronary circulation (Fig. 1) are also available in CARDIOSIM© platform [3].
For the purposes of our study, we have assembled the cardiovascular network with the new module of the systemic circulation whilst the behavior of the heart is modelled as described in [3],[4]. The systemic venous section [3],[6] and the entire pulmonary circulation [7],[8],[11],[12] are modelled as described in current literature. We have selected the module presented in [13] for the coronary circulation. The tricuspid, mitral, pulmonary, and aortic valves are modelled with an ideal diode: when the pressure across the valve is positive, the valve opens and allows the flow of blood; when the pressure is less than or equal to zero, the valve is closes and the flow of blood is zero [3]-[5],[14]
New lumped-parameter model of the systemic circulation
The new module is described in Fig. 1 using resistance, inertance and compliance (RLC) elements. The systemic network consists of the following compartments: ascending and descending aorta, carotid artery, peripheral arteries and superior and inferior vena cava (Fig. 1a).
The ascending (descending) aorta is modelled using two resistances RAA and RVaa (RDA and RVda), inertance LAA (LDA) and compliance CAA (CDA). A single resistance (RCA) reproduces the behavior of the carotid district. The peripheral arterial circulation is reproduced with resistances RA and RVa and with compliance CA (Fig. 1b). The superior vena cava compartment consists of resistances RSVC_I and RSVC_II, inertance LSVC and compliance CSVC. The systemic venous network is implemented using compliance CVS and resistances RVS and RVvs (Fig. 1a). Finally, the inferior vena cava district is modelled with resistances RIVC, RIVC_I and RIVC_II, inertance LIVC and compliance CIVC (Fig. 1a). The resistances RVaa, RVda, RVa, and RVvs account for viscous losses of the vessel wall. Pt is the mean intrathoracic pressure.
Right ventricular assist device (RVAD)
A 0-D numerical model described in [7] was used to implement a centrifugal pump that reproduced the behavior of the right ventricular assist device. The pump can be connected in series or parallel to the right ventricle with two cannulae modelled using RLC elements (Fig. 1c). When the RVAD removes blood from the right atrium (parallel connection – SW1=ON and SW2=OFF in Fig. 1c), the flow through the inlet cannula is: RAP is the right atrial pressure, LiCANN, RiCANN and CiCANN are the inertance, resistance and compliance of the inlet cannula, respectively (Fig 1c). QiCANN (QPUMP) is the flow through the inlet cannula (generated by the centrifugal pump). ΔP is the pressure on the pump head.
When the RVAD removes blood from the right ventricle (in series connection - SW1=OFF and SW2=ON in Fig. 1c), Eq. 3 becomes: RVP is the right ventricular pressure.
When the RVAD ejects blood, the flow through the outlet cannula is: MPAP is the mean pulmonary artery pressure (Fig. 1c), RoCANN and CoCANN are the inertance, resistance and compliance of the outlet cannula, respectively. QoCANN is the flow through the outlet cannula.
Simulation Protocol
For the purposes of our simulations, we considered values for right and left ventricular elastance that could reproduce a realistic diseased heart. According to the available literature, normal values for right ventricular elastance fluctuate around 1 mmHg/ml [15],[16] and range from 1.6 to 5 mmHg/ml for left ventricular elastance [17],[18]. Therefore, we considered EesLEFT = 0.7 mmHg/ml for the left ventricle and EesRIGHT = 0.3 mmHg/ml for the right ventricle as reference values for a failing heart. Our simulation approach consisted of three steps. In the first step after setting HR= 90 bpm, the slope of left ventricular End-Systolic Pressure-Volume Relationship (ESPVR) EesLEFT=0.7 mmHg/ml and the slope of right ESPVR EesRIGHT=0.3 mmHg/ml [8], the simulator generated the following values: cardiac output (CO) 4.51 l/min, aortic systolic (diastolic) pressure 82.1 (60.4) mmHg, mean right atrial pressure 23.3 mmHg, pulmonary systolic (diastolic) pressure 51.2 (31.6) mmHg and pulmonary capillary wedge pressure 20.7 mmHg. The mean pressure (flow) value was calculated as the mean value of all blood pressure (flow) measurements during a cardiac cycle.
In the second step, the slope of right ESPVR EesRIGHT=0.3 mmHg/ml was set to 0.4, 0.5, 0.6, 0.7 and 0.8 mmHg/ml [1],[19] and for each value the parameters described above were measured.
In the third step, RVAD support was applied both in series and parallel mode with rotational speed of 2000, 2500, 3000, 3500 and 4000 rpm. In the fourth step, the slope of right ESPVR was changed from 0.3 to 0.8 mmHg/ml (0.1 mmHg/ml stepwise increase) during RVAD support connected in series and parallel mode. The measured parameters were those described above.
Considering that we did not measure data in patients undergoing RVAD support to allow us to reproduce their haemodynamic conditions with our numerical simulator, and given that literature data are largely incomplete, we decided against a direct data comparison. Therefore, we considered percentage variation to evaluate the trend of the effects induced by RVAD support on the haemodynamic and energetic variables in line with other available simulation work.
RESULTS
Figure 2 shows the effects induced by different values for ESPVR slope (EesRIGHT) on left and right atrial pressure-volume loop area (left upper panel). The percentage changes calculated with respect to the reference value for EesRIGHT=0.3 mmHg/ml have been listed. The red bars show the percentage change with respect to the reference value for PVLAL-A increase at high EesRIGHT values. On the contrary, the PVLAR-A (yellow bars) decreases when the percentage change is referred to high EesRIGHT values. These results show that when EesRIGHT increases from 0.3 to 08 mmHg/ml, the area of the pressure-volume loop of the right atrium decreases whilst the area of the pressure-volume loop of the left atrium increases (Fig. 3). The left and right lower panel shows the left and right atrial pressure-volume loop obtained with different EesRIGHT values. The black pressure-volume loop was obtained setting EesRIGHT=0.3 mmHg/ml, while the blue pressure-volume loop was obtained setting EesRIGHT=0.6 mmHg/ml. The left (left side) and the right (right side) ventricular pressure-volume loops are placed in the upper panels. When EesRIGHT increases, EWL_V and PVAL-V (left upper panel) increase leading to right-sided shift with increased left ESV and EDV. The effect induced by different EesRIGHT values on left and right ventricular EW (EWL-V and EWR-V) are reported in Fig. 2 (left lower panel). The effect induced on EDV and ESV of the left (right) atrium is reported in the upper (lower) right panel (Fig. 2). The effect induced on right atrial ESV is more evident than the one produced on right atrial EDV. Figure 3 (right upper panel) shows that an increase in EesRIGHT leads to an increase in right ventricular external work and PVAR-V with left-sided shift and decrease of both right ESV and EDV. Figure 4 shows the effect induced by RVAD support driven in parallel connection (left panels) on left ventricular PVAL-V (left upper panel) and EWL-V (left lower panel). The effect induced on PVAL-V and EWL-V when the RVAD is driven in series is available in the right panels. The data were measured for different values of EesRIGHT (0.3; 0.5 and 0.8 mmHg/ml) and with increasing RVAD rotational speed (2000; 2500; 3000; 3500 and 4000 rpm). For each rotational speed and EesRIGHT values, the percentage changes calculated with respect to the reference value measured in pathological conditions have been listed (Fig. 4). The highest percentage changes in PVAL-V and EWL-V (50% or more for both variables) were recorded in a diseased condition with right ventricular elastance set to EesRIGHT=0.3 mmHg/ml during in-parallel RVAD assistance with pump rotational speed at 4000 rpm. Figure 5 shows the effect induced by in-parallel RVAD support on left ventricular ESVL-V (left upper panel) and EDVL-V (left lower panel). The effect induced on ESVL-V (right upper panel) and EDVL-V (right lower panel) by in-series RVAD support is also available for comparison purposes. When the rotational speed of the pump was set to 2000 rpm, the percentage changes for each value of right ventricular ESPVR slope were no more than 3%. The effects induced by different rotational speed of the RVAD on left and right atrial pressure-volume loop area are available in Fig. 6 for three different values of right ESPVR slope. In addition, for each value of the slope, the effect induced by the pump has been calculated in percentage terms with respect to the value obtained without assistance. In the case of the right atrium, we observed a higher percentage reduction of the loop area (from 2% to 9%) during in-series assistance (left lower panel). In the case of parallel connection (left upper panel), the percentage reduction of the pressure-volume loop area of the right atrium becomes more evident at high values for both the slope and the pump rotational speed. Right heart assistance produces significant effects on the left atrium loop area (right panels). Percentage variations between 20% and 50% are observed when the RVAD is connected in parallel to the right ventricle (right upper panel). Whether in-series or in-parallel assistance is considered, the most significant percentage variations are observed when the right ventricular ESPVR slope is set to EesRIGHT=0.3 mmHg/ml. The effects induced by RVAD support on the right ventricular flow output and on the total cardiac output are showed in Fig. 7. RVAD parallel assistance generated a lower reduction in right ventricular output compared to in series assistance for each EesRIGHT value and each pump rotational speed. More specifically, when the pump speed was set to 4000 rpm, parallel assistance induced a reduction between 20% and 45% compared to the baseline value while the in series assistance induced a reduction between 38% and 78%. The upper panels in Fig. 7 show that parallel RVAD assistance generates a higher increase in total cardiac output compared to in series assistance at pump speeds higher than 3500 rpm (for each EesRIGHT value).
The left upper (lower) panel in Fig. 8 shows the effects induced on left ventricular (atrial) loop when different types of RVAD assistance were applied to pathological conditions (black loops) reproduced by setting EesRIGHT to 0.3 mmHg/ml. The assistance was applied in parallel (red loops) and in series (blue loops) mode with rotational speed of 3000 rpm.
The right upper (lower) panel shows the effects induced on right ventricular (atrial) loops. When in-parallel assistance was applied, a right-sided shift in the left and right ventricular loop (red lines) was observed. Although in-series RVAD assistance did not cause changes in right end-systolic ventricular volume, it leaded to a reduction in right end-diastolic ventricular loop. The type of assistance did not cause relevant changes in the right atrial loop (right lower panel).
DISCUSSION
The left ventricle (LV) is coupled to the low-compliance, high-resistance peripheral arterial circulation and is more adaptable to changes in pressure than volume. In contrast, the right ventricle (RV) is coupled to the high-compliance, low-resistance pulmonary circulation and is more adaptable to changes in volume than pressure. The right ventricle consists of a free wall containing a wrap-around circumferential muscle at its base and a septum made of oblique helical fibres crossing each other at 60° angles. This is consistent with the helical ventricular myocardial band concept, which defines two interconnected muscle bands: a basal loop with transverse fibres surrounding the left and right ventricles and an apical loop made of a right- and left-handed helix forming an apical vortex [20],[21]. The wrap-around transverse fibres constrict or compress leading to a bellows motion responsible for 20% of right ventricular output whilst the oblique fibres are responsible for shortening and lengthening, which contribute to 80% of right ventricular systolic function [22]. The crista supra-ventricularis shares muscle fibers with the inter-ventricular septum and the free wall playing a key anatomical and functional role [23]. A reduction in longitudinal contraction and an increase in transverse shortening are observed following cardiopulmonary bypass and pericardiotomy [24]. This is quite an important aspect to bear in mind and may be addressed initially with pulmonary vasodilators [25]. The relationship between structure and function plays a key role in clinical decision-making, which must be based on detailed knowledge of normality and recognise how disease can be addressed to restore normality [22]. The important contribution of right ventricular function has been neglected for a long time due to previous observations and assumptions. The onset of right ventricular dysfunction should trigger the search for the main underlying cause in relation to pressure overload, volume overload or primary myocardial disease [26]. Right heart failure (RHF) is difficult to manage because of its complex geometry and a lack of specific treatments aimed at stabilisation and recovery of right ventricular function. Nevertheless, right ventricular dysfunction remains associated with poor clinical outcome regardless of the underlying disease mechanism [27].
A simulation approach overcomes ethical issues and the risk of offering an ineffective or potentially dangerous therapeutic option. At the same time, it may help focus on the specific problem to address. Our starting point was to develop a failing right heart, which would require support at a subsequent stage. The easier way to do it was to act on the ESPVR slope of both ventricles. A range between 0.35 and 0.74 mmHg/ml is observed in patients with pulmonary hypertension [1] with a cut-off of 0.8 for £eS/£a ratio as the onset of right ventricular maladaptation. Our initial aim was to observe the effects of RVAD support with either in series or parallel connection following stepwise variation of right ventricular end-systolic elastance in patients with increased right ventricular afterload. Therefore, the right ventricular end-systolic elastance considered in the present study ranged between 0.3 and 0.8 mmHg/ml as per previously reported values observed in clinical practice [1],[28],[29]. According to [19], increased native cardiac output is observed in the presence of left ventricular systolic impairment when right ventricular end-systolic elastance increases from 0.1 to 1.0 mmHg/ml. An increased native cardiac output is still observed during VA ECMO support following stepwise increase in right ventricular end-systolic elastance but to a lesser degree. A left-to-right ventricular septal shift is observed during diastole following a stepwise increase in right ventricular end-systolic elastance both with and without VA ECMO support. Our considerations are based in the context of pure RVAD support with either in series or parallel connection.
PVA and EW of the left ventricle gain benefit when the pump speed of the RVAD is at least 3500 rpm. The highest effect is obtained when EesRIGHT is 0.3 mmHg/ml, which is consistent with significant RV dysfunction (left panels in Fig. 4). The pathological range considered would suggest that early recognition and aggressive RVAD support is advisable where a lower degree of assistance is required generating enough benefit for a less compromised right ventricle with more potential for recovery. This is an important point to consider in the context of LVAD support when the right heart shows signs of failure requiring attention.
The highest beneficial effect is obtained when the RVAD is connected in parallel to the right ventricle with up to nearly 35% increase in total cardiac output (right upper panel in Fig. 7) and lower reduction of right ventricular output compared to in series RVAD connection. No significant effect is observed on the right atrium regardless the type of RVAD connection to the right ventricle. Instead, in-parallel RVAD connection has a more beneficial effect on the left atrium. Again, RVAD support has a more beneficial effect on ESV and EDV of the left ventricle when connected in parallel to the right ventricle (left panels in Fig. 5). The role of the inter-ventricular septum is critical in this context.
Our aim was to observe the effect of pure RVAD assistance at different stages of right ventricular dysfunction to determine appropriate timing for intervention. Our target was early graft failure secondary to right heart dysfunction following orthotopic heart transplant and right ventricular failure following LVAD insertion in an apparently preserved right heart function preoperatively. We have focused our attention mainly on EesRIGHT neglecting Vo. We have not considered the progressive increase in afterload. Despite these limitation, our preliminary findings support the concept of early intervention in the presence of a failing right heart regardless of its aetiology. This simulation study confirms what had been previously advocated but not always put into practice [30]. A more liberal right ventricular support may be the way forward [31] when considering different support strategies for a failing right ventricle [32]. Late onset of right ventricular failure remains associated with worse survival and higher cumulative incidence of major adverse events [33].
CONCLUSION
Although RVAD support may be effective in advanced right heart failure, intervention at an earlier stage may require a lesser degree of assistance and be more beneficial for the left ventricle. Early recognition of right ventricular failure followed by aggressive treatment is desirable to achieve a more favourable outcome. RVAD support remains an option for advanced right ventricular failure although the onset of major adverse events may preclude its use. In-parallel RVAD connection to the right ventricle seems a more suitable option.
Data Availability
All data produced in the present study are available upon reasonable request to the authors
FUNDING AND CONFLICT OF INTEREST
None
ACKNOWLEDGEMENTS
None
Footnotes
beatrice.delazzari{at}gmail.com
aiacovoni{at}asst-pg23.it
capoccia{at}doctors.org.uk
silvia.papa{at}uniroma1.it
roberto.badagliacca{at}uniroma1.it
dario.vizza{at}gmail.com
claudio.delazzari{at}ifc.cnr.it