Personalized surgical planning for coronary bypass graft configurations using patient-specific computational modeling to avoid flow competition in arterial grafts

Objectives Flow competition between coronary artery bypass grafts (CABG) and native coronary arteries is a significant problem affecting arterial graft patency. The objectives of this study were to compare the predictive hemodynamic flow resulting from various total arterial grafting configurations and to evaluate whether the use of computational fluid dynamics (CFD) models capable of predicting flow can assist surgeons to make better decisions for individual patients by avoiding poorly functioning grafts. Methods Sixteen cardiac surgeons declared their preferred CABG configuration using bilateral internal mammary and radial arteries for each of 5 patients who had differing degrees of severe triple vessel coronary disease. Surgeons selected both a preferred 'aortic' strategy, with at least one graft arising from the ascending aorta, and a preferred “anaortic” strategy which could be performed as a “no-aortic touch” operation. CT coronary angiograms of the 5 patients were coupled to CFD models using a novel flow solver “COMCAB.” Twelve different CABG configurations were compared for each patient of which 4 were “aortic” and 8 were “anaortic.” Surgeons then selected their preferred grafting configurations after being shown predictive hemodynamic metrics including functional assessment of stenoses (instantaneous wave-free ratio; fractional flow reserve), transit time flowmetry graft parameters (mean graft flow; pulsatility index) and myocardial perfusion. Results A total of 87.5% (7/8) of “anaortic” configurations compared to 25% (1/4) of “aortic” configurations led to unsatisfactory grafts in at least 1 of the 5 patients (P = 0.038). The use of the computational models led to a significant decrease in the selection of unsatisfactory grafting configurations when surgeons employed “anaortic” (21.25% (17/80) vs. 1.25% (1/80), P < 0.001) but not “aortic” techniques (5% (4/80) vs. 0% (0/80), P = 0.64). Similarly, there was an increase in the selection of ideal configurations for “anaortic” (6.25% (5/80) vs. 28.75% (23/80), P < 0.001) but not “aortic” techniques (65% (52/80) vs. 61.25% (49/80), P = 0.74). Furthermore, surgeons who planned to use more than one unique “anaortic” configuration across all 5 patients increased (12.5% (2/16) vs. 87.5% (14/16), P<0.001). Conclusions “COMCAB” is a promising tool to improve personalized surgical planning particularly for CABG configurations involving composite or sequential grafts which are used more frequently in anaortic operations.

Objectives: Flow competition between coronary artery bypass grafts (CABG) and native coronary arteries is a significant problem a ecting arterial graft patency. The objectives of this study were to compare the predictive hemodynamic flow resulting from various total arterial grafting configurations and to evaluate whether the use of computational fluid dynamics (CFD) models capable of predicting flow can assist surgeons to make better decisions for individual patients by avoiding poorly functioning grafts.
Methods: Sixteen cardiac surgeons declared their preferred CABG configuration using bilateral internal mammary and radial arteries for each of patients who had di ering degrees of severe triple vessel coronary disease. Surgeons selected both a preferred 'aortic' strategy, with at least one graft arising from the ascending aorta, and a preferred "anaortic" strategy which could be performed as a "no-aortic touch" operation. CT coronary angiograms of the patients were coupled to CFD models using a novel flow solver "COMCAB." Twelve di erent CABG configurations were compared for each patient of which were "aortic" and were "anaortic." Surgeons then selected their preferred grafting configurations after being shown predictive hemodynamic metrics including functional assessment of stenoses (instantaneous wave-free ratio; fractional flow reserve), transit time flowmetry graft parameters (mean graft flow; pulsatility index) and myocardial perfusion.

. Introduction
Total arterial grafting (TAG), which involves the exclusive use of arterial conduits such as bilateral internal mammary arteries (BIMA) and radial arteries (RA) for coronary artery bypass grafting (CABG), has been associated with improved long-term clinical outcomes (1). Total arterial, anaortic, off-pump coronary artery bypass grafting (OPCABG) is an operative approach that avoids aortic manipulation with no graft attached to the ascending aorta and thus this strategy has been endorsed for its reduction in stroke risk and other complications associated with cardiopulmonary bypass (2). However, anaortic grafting configurations typically require more composite and sequential grafts rather than separate grafts.
Composite BIMA grafting to left coronary targets with <70% diameter stenosis leads to a high rate of graft occlusion or constriction due to competitive flow (3). For a radial artery graft anastomosed to the right coronary artery (RCA) territory, the percent diameter stenosis should be at least 80-90% because the radial artery is more prone to spasm (4,5). Functional assessment of coronary stenoses using fractional flow reserve (FFR) < 0.80 or instantaneous wavefree ratio (iFR) < 0.90 (6), and functional assessment of bypass grafts using transit-time flowmetry (TTFM) with mean graft flow (MGF) ≥ 15 ml/min and pulsatility index (PI) < 5 can avoid situations of competitive flow leading to poor graft patency following arterial grafting (7).
There is no consensus among surgeons regarding the optimal TAG configuration using BIMA for an individual patient and it has even been argued that configuration is not important (8). However, others have maintained that graft configuration is significant (9) and that suboptimal judgement in the arrangement of grafts can lead to steal of flow between grafts and native coronary arteries particularly when composite or sequential grafts are used with unbalanced native coronary stenoses (10). In this study, the importance of graft configuration was investigated using predictive patient-specific hemodynamic computational modeling and the impact of such predictive information on surgical planning was evaluated.
. Materials and methods

. . Selection of patient cases
Institutional ethics review was obtained to conduct this study. Five patient cases for inclusion in the study were identified by screening patients who had undergone both a CT coronary angiogram (CTCA) and invasive coronary angiogram. All patients were required to have severe triple-vessel coronary artery disease with >90% diameter stenosis in the RCA territory and >75% stenosis in the left anterior descending (LAD) and circumflex (CIRC) territories such that surgeons would not deliberately avoid using a RA graft on a distal target due to concerns of competitive flow. Patients 1 and 2 both had a ramus intermedius artery ( Figures 1A, B), Patients 3 and 4 had additional less severe stenoses ( Figures 1C, D), and Patient 5 had an in-stent restenosis in the LAD ( Figure 1E). The differing degrees of severity and distribution of stenoses evident from the coronary tree diagrams of the five patients are summarized in Table 1.

. . Selection of grafting configurations
Twelve grafting configurations utilizing BIMA and RA, used by surgeons internationally, were examined for each of the 5 patient-specific diseased coronary circulations. Four of the grafting configurations involved aorto-coronary grafts ("aortic" configurations). One configuration using BIMA exclusively, recycles the free RIMA for use to the RCA by anastomosing it from the aorta while using a LIMA/RIMA composite graft for the leftsided vessels (11) (configuration A) (Figure 2A). Other configurations have three separate inflows with separate in situ use of both internal mammary arteries (IMA) and the RA off the aorta. The RIMA can be used in situ via the transverse sinus to the CIRC branches (12) (configuration B) ( Figure 2B), in situ to the distal RCA or proximal posterior descending artery (PDA) (13) (configuration C) ( Figure 2C), or, if used for the LAD, then the LIMA is anastomosed to the CIRC branches (14) (configuration D) ( Figure 2D).
The other eight grafting configurations involved no proximal graft anastomosis to the aorta ("anaortic" configurations). A composite Y-graft can be constructed using the RA off the in situ LIMA to LAD and if the in situ RIMA cannot reach the RCA target then it is lengthened by a small I-graft with the remaining RA (15) (configuration E) ( Figure 2E). An alternative configuration involves the RIMA in the Y-composite graft which spares the LIMA to be used as an individual graft (16) (configuration F) ( Figure 2F). Certain surgeons perform a double-Y graft to revascularize all three territories and this is based on a single inflow (17) (configuration G) ( Figure 2G). Another composite approach used by surgeons, based on a single inflow requiring use of the BIMA only, is a configuration that uses the free RIMA as a sequential graft (11, 18) (configuration H) ( Figure 2H).
The appeal in using both IMA as in situ grafts, has led surgeons to using the RA as an I-graft to lengthen the RIMA and configure it .
/fcvm. .  as a sequential graft to anastomose the CIRC and RCA territories. The in situ RIMA-RA sequential can be used in an anticlockwise orientation where the first anastomosis is to the RCA and the final anastomosis is to the CIRC branches (19) (configuration I) ( Figure 2I). The orientation of the RIMA-RA sequential in the clockwise version is such that it is first brought through the transverse sinus to anastomose to the CIRC branches and then terminates on the RCA (20) (configuration J) ( Figure 2J). Another configuration uses an I-graft constructed with the RIMA and RA but uses it as a free Y-graft from the LIMA to CIRC (21) (configuration K) ( Figure 2K). Finally, an uncommon "bail-out" configuration involves the use of a jump graft where the radial artery conduit is anastomosed from the distal LAD to the PDA (22) (configuration L) ( Figure 2L).

. . Computational model analysis of grafting configurations
Patient-specific 1D-0D computational fluid dynamics models were created for each of the 5 patients' theoretical non-diseased coronary circulation, their stenotic circulation and their 12 grafted circulations and these 70 networks were solved using the novel software "COMCAB" created by the authors for this purpose. This involved a manual segmentation of the native coronary artery geometry from the CTCA for each patient and the mapping of each vessel onto the 1D domain. Thereafter, coronary artery side branches were added to account for a physiological loss in pressure from proximal to distal along the coronary arteries as such branches were not visualized on the CTCA. Terminal vessels were coupled to 0D lumped parameter models. Lumped parameter stenoses were added to create the stenotic network for each patient and 1D grafts added to create the different grafting network topologies for the grafted circulations. The networks were solved for blood flow and pressure using the Richtmyer two-step Lax-Wendroff 1D numerical method (23,24) with prescription of appropriate boundary conditions. The chosen meshgrid size, x, was set at x ≈ 0.1cm and the timestep, t, was chosen to be the largest value possible whilst ensuring stability of the numerical scheme with t ≈ 2.95 × 10 −5 . The methodology for the creation of the CFD models has previously been described in detail in the literature (25). The vessel segment .
/fcvm. . The computational model inputs included a generic aortic root pressure waveform of 120/77 mmHg at 65 beats per minute (cardiac period 0.917 s), a generic left ventricular (LV) pressure waveform, and cardiac output (CO) of 5 L/min. Approximately 4.5% of CO was assigned to myocardial blood flow with this flow distributed to the three main coronary territories: LAD, CIRC, and RCA, as dictated by distribution of vessels on the CTCA. Total arterial compliance was estimated at 1.15 cm 3 /mmHg from the aortic pressure waveform and a 3-element Windkessel resistance-capacitance-resistance (RCR) model was applied at the outlets of the terminal vessels to represent the distal microcirculation (25). iFR was chosen as the metric for functional stenosis severity as it is measured at rest like TTFM graft measurements, not at hyperemia which is used for FFR. However, as iFR is a recent concept an equivalent FFR was also calculated: (26). Creation of the predictive models took ∼6 h for each patient. Running all simulations using high performance computing on Xeon Broadwell CPUs (2.1 GHz) with embarrassing parallelisation, took up to 75 min. MGF and PI were calculated at the distal end of each graft ( Figure 3). Myocardial territory perfusion was calculated for the theoretical non-diseased circulations and the diseased circulations along with the improvement in myocardial perfusion following restoration of blood flow by the grafted coronary circulations. Validation of the computational models' predictions was performed by comparing the calculated iFR, MGF and PI with in vivo measurements available in the literature.
Grafting configurations were classified as unsatisfactory, satisfactory, or ideal, based on the graft performance indices. An unsatisfactory grafting configuration was defined as having either MGF < 15 ml/min or PI > 5 in any graft. A satisfactory grafting configuration was defined as MGF ≥ 15 ml/min in all grafts with 3 < PI < 5 for a graft to a LV target and an ideal grafting configuration was defined as MGF ≥ 15 ml/min and PI < 3 in all grafts (27).

. . Selection of participant surgeons
Sixteen cardiac surgeons were recruited for participation in the study from five centers. The power calculation assumed that ∼50% of the selected grafting configurations would be deemed unsatisfactory with standard surgeon decision-making but close to 0% after the provision of computer modeling predictions. Thus, the study would require a sample size of 16 to achieve two-tailed statistical significance with an alpha level of 0.05 and a power of 0.8. No dropouts were expected due to the administration of the survey in one sitting. Surgeon experience was noted by recording the volume of CABG operations performed in their career and current practice of total arterial revascularisation using BIMA.

. . Survey administration procedures
The coronary tree diagrams, as representations of each patient's diseased coronary circulation, were presented to the surgeons . /fcvm. .    ( Figure 1). Surgeons were informed that the only available conduits would be the LIMA, RIMA and one RA with no contraindications for their use. All patients had an overall preserved LV function with all distal grafting target sites being suitable for grafting. Conduit lengths were adequate for an in situ RIMA to reach the obtuse marginal (OM) target and for an in situ RIMA to reach a distal RCA/PDA. Each of the surgeons selected one preferred "aortic" configuration and one preferred "anaortic" configuration for each of the five patients (standard decision-making). They were then asked to indicate their overall grafting preference if there were no constraints on their choice such as a porcelain aorta. The surgeons were shown the 12 grafting configurations used by surgeons internationally and asked to rank their top five selections as well as indicate any configurations they would not use. They were then presented with the predictive hemodynamic graft flow information from the computer model and asked for their preferred selections before, and after, the provision of regional myocardial perfusion data (computer-informed decision-making) (Tables 2-6).

. . Outcomes of interest
The outcomes of interest were the selected number of unsatisfactory graft configurations (primary outcome measure), the number of ideal graft configurations, and the number of unique grafting configurations among the five patients for both "aortic" and "anaortic" configurations. Additional exploratory measures included the number of configuration rankings that changed with additional regional myocardial perfusion data and the number of configurations that individual surgeons indicated they would not initially use but later would use after viewing the computational models' hemodynamic predictions.

. . Statistical analysis
Differences in outcome measures were compared as continuous variables for flows with paired t-tests and as proportions for discrete variables with chi-squared tests. If parametric assumptions were not satisfied, then the Mann-Whitney U-test for independent comparisons was performed or the Wilcoxon signed rank-test for dependent comparisons. For the primary outcome measure, McNemar's test was performed to determine the influence of the computational model in changing individual surgeons' selection of unsatisfactory grafting configurations. Outcomes were considered statistically significant for P < 0.05 and the statistical tests were twotailed with power 0.8 to test hypotheses. A learning gain attributable the computational model was set at 30% of surgeons improving their unsatisfactory decisions, which is considered effective for an educational intervention (28).

. . Validation of computational fluid dynamics model predictions
Seven out of the 9 stenoses (78%) in this study that were between 71 and 90% diameter stenosis were functionally significant ( Figure 4A). This proportion correlated with the findings of a large clinical study in which 80% of such lesions were found to be functionally significant (29). The MGF and PI predicted by the computational models compared well with available in vivo TTFM measurements from other clinical studies, noting that the measurement of PI depends on location along the graft and is graft-specific (Table 7). Furthermore, the correlation between MGF and iFR (R = −0.455) ( Figure 4B) as well as PI and iFR (R = 0.522) for 15 separate arterial grafts in this study ( Figure 4C) was consistent with an in vivo clinical study involving 25 arterial grafts where the correlation coefficients were R = −0.460 and 0.563, respectively (30).

. . Patient-specific grafting configuration predictions
For Patient 1, configurations G and H had MGF < 15 ml/min in the composite LIMA to LAD graft mainly due to lower flows down this limb where the LIMA was the single inflow source for grafts to all three targets ( Table 2). For Patient 2, configurations A, F, G, I and K had MGF < 15 ml/min and higher PI because composite grafts to the OM1 were subject to steal of flow, which was . /fcvm. .  exacerbated by flow competition from a native OM1 75% stenosis with iFR 0.86 (Table 3). For Patients 3 and 4 only configuration L was unsatisfactory due to 90% and 95% stenosis, respectively, in the RCA that compromised the MGF and PI in the RA to PDA jump grafts (Tables 4, 5). For Patient 5, configurations A, E, G, and H were unsatisfactory with MGF < 15 ml/min in the LIMA to LAD graft due to a steal of blood flow down the other limb of the composite Y-graft, which was accentuated by the native competitive flow in the LAD due to a stenosis of 90% with an iFR 0.93 (Table 6). For Patient 5, configuration L was unsatisfactory due to a reversal of flow in the RA to PDA jump graft from a 90% RCA stenosis combined with limited blood flow supplied by the upstream LIMA to LAD (Table 6).

. . Standard surgeon decision-making
There was a wide range of cardiac surgical experience among the participating surgeons (Table 8). The most common "aortic" grafting strategy selected with standard decision-making was configuration C (62.5%), whereas the most common "anaortic" strategy was configuration E (71.25%) and no surgeons in this study chose configurations G or L with very few selecting configurations I or K ( Table 9). As few surgeons selected configuration A [15% (12/80)] which was the only "aortic" configuration that led to unsatisfactory grafting in two patients, standard decision-making led to very few unsatisfactory "aortic" grafting configurations [5% (4/80)] (Table 9, Figure 5). However, as most surgeons selected "anaortic" configuration E [71.25% (57/80)] which led to unsatisfactory grafting in one patient as well as "anaortic" configuration H [13.75% (11/80)] with unsatisfactory grafting in two patients the rate of unsatisfactory "anaortic" selections with standard decision-making was higher [21.25% (17/80)] (Table 9, Figure 5). An "aortic" configuration was preferred over an "anaortic" configuration, 77.5% (62/80) of the time. Half (8/16) of the surgeons chose at least one "anaortic" strategy among the five patients in preference to an "aortic" configuration.

. . Impact of computational models on surgeon decision-making
As standard surgeon decision-making led to the selection of fewer unsatisfactory "aortic" grafting configurations compared with "anaortic" configurations, the integration of the computational model-generated predictions by the cardiac surgeons led to a significant decrease in the selection of unsatisfactory grafting configurations for "anaortic"    (Table 9). Similarly, there was an increase in the selection of ideal configurations for "anaortic" [6.25% (5/80) vs. 28.75% (23/80), P < 0.001] but not "aortic" techniques [65% (52/80) vs. 61.25% (49/80), P = 0.743] ( Table 9). For the primary outcome measure, the number of surgeons that changed from choosing at least one unsatisfactory grafting configuration to having no unsatisfactory grafting configurations was significant for "anaortic" configurations (13/16 = 81.25%, P < 0.001) but not for "aortic" configurations (2/16 = 12.5%, P = 0.480) ( Table 10). The computational model predictions also led surgeons to make more patient-specific grafting selections for the "anaortic" configurations. The number of surgeons that changed from choosing no unique grafting configurations to having at least one unique grafting configuration across the five patients was significant (12/16 = 75%, P = 0.002) ( Table 10). Half (8/16) of computer-informed surgeons used a configuration that they earlier stated they would not use and 56.25% (9/16) changed the order of their ranking preferences with the addition of myocardial perfusion data.

. . Role for predictive computational flow modeling in surgical planning
After the degree of functional coronary stenosis, graft configuration has the most significant influence on graft and native coronary artery flows (10). Despite the tendency of surgeons for using a "one-size fits all" approach, hemodynamic predictions from this study suggest that surgeons need to tailor composite and sequential grafting configurations for each individual patient. Anaortic configurations based on a smaller number of separate inflows are more prone to steal of flow and competitive flow affecting upstream segments (17).
"Anaortic" configuration J, with its sequential arrangement, was favorable for Patients 1, 2, 3, and 4 as the RCA stenoses were equal to or greater than the CIRC stenoses (31). In Patient 5, although the OM vessel stenosis was at 99% and the RCA at 90% (iFR 0.73), the flows in the graft segment to PDA were still adequate at 15.29 ml/min. A higher iFR for the RCA stenosis in this patient could compromise configuration J as being universally satisfactory. This observation highlights that the complexity of determining unsatisfactory grafting configurations is potentially beyond the capabilities of using simple heuristics and may require the quantitative predictive value of computational modeling to determine the complex interplay between graft and host vessels and distal run-off.
Standard surgeon decision-making culminated in four composite Y-graft "aortic" selections [5% (4/80)] and 17 "anaortic" selections [21.25% (17/80)] which, on the basis of predictive hemodynamics, could lead to unsatisfactory grafts with poor patency. This result reveals a significant clinical opportunity for improvement in graft selection and function for both surgeons and patients given that options exist to avoid these configurations. In a clinical study of 120 patients undergoing complex composite grafting procedures using BIMA, the graft patency at a mean follow-up of 29.9 ± 33.1 months for arterial grafts was between 80 and 98.7% (32). Although grafts may appear to be patent initially, those with poor flow have been found to be occluded within 1 year (33, 34).
The computational model was effective in reducing the proportion of surgeons choosing unsatisfactory "anaortic" configurations from 87.5 to 6.25% (P < 0.001), far exceeding the learning change set a priori at 30% (28). One surgeon continued to pursue an unsatisfactory grafting configuration for one patient which may be explained by cognitive biases when evaluating computational models (35). The computational model influenced 50% of surgeons to select grafting configurations that they stated they would not usually use and led to more patient-specific tailoring of . /fcvm. .

Relationships between metrics of functional stenoses and grafts. (A)
The majority of stenoses > % diameter were functionally significant at rest (iFR < . ). Two stenoses between and % were also functionally significant on account of their increased lesion length. (B) Correlation between mean graft flow and iFR (R = − .
) and (C) correlation between pulsatility index and iFR (R = . ) using the arterial grafts across patients using separate arterial grafts from configuration C, as an example. iFR, instantaneous wave-free ratio; R , coe cient of determination.
"anaortic" configurations, rather than a "one-size fits all" approach. This change in decision-making is understandable as the biophysics of the coronary circulation is complex, often difficult to completely measure and highly patient-specific. Thus, the standard decisionmaking for a CABG configuration is made in an environment of significant uncertainty (36). This study demonstrates that, when faced with uncertainty, the additional quantified flow information provided by predictive computational modeling will change surgical planning especially for configurations involving composite grafts.
The provision of additional myocardial perfusion data led 56.25% (9/16) of surgeons to change their graft configuration rankings. Indeed, surgeons tend to be poor at predicting  (25). Values are presented as mean ± standard deviation. MGF, mean graft flow; PI, pulsatility index; LIMA, left internal mammary artery; RIMA, right internal mammary artery; RA, radial artery; LAD, left anterior descending coronary artery; CIRC, circumflex coronary artery; RCA, right coronary artery; PDA, posterior descending coronary artery. qualitative effects on perfusion post CABG (37) and are not accustomed to assimilating data providing quantification of myocardial perfusion.

. . Limitations and future directions
There were a number of limitations in this study relating to how surgeons were engaged. Using only five patients did allow 12 grafting configurations to be interrogated on each patient. While relevant predictive hemodynamic data were presented to all 16 surgeons, other TTFM parameters which "COMCAB" is capable of measuring such as diastolic filling percentage and backward flow were not provided, to avoid information overload. This would not affect results as this study addressed situations of flow competition rather than technical anastomotic errors and hence providing MGF and PI were considered sufficient (38). Despite this, the study still typically took up to 45 min for each surgeon to complete. As this virtual surgical planning study omitted the real-life execution of CABG, the authenticity of surgeon responses may also be questioned. Only four  variability is consistent with 20% TAG utilization in Europe, up to 80% in some Australian centers (39) and 3% overall BIMA usage on multi-institutional analysis (40). A further subgroup analysis is planned as an extension to this study to investigate the effect of surgeon experience on engagement with information provided by the computational modeling. Future studies could also be streamlined by including more surgeons who perform anaortic total arterial OPCABG revascularisation using BIMA or even CABG involving more composite and sequential grafts. The limitations arising from computational modeling have been described previously (25). With more invasive clinical patient data available, the idealized generic parameters of aortic root pressure, LV pressure, heart rate, CO, total arterial compliance and systemic arterial branch dimensions could be made more patient-specific. Despite this, the differences in patient-specific coronary artery disease resulted in differing outcomes from different grafting configurations between patients. The 7 h duration required to create the computer models included manual processing and therefore further automation of methods would be beneficial. Although hemodynamic predictions generally agreed with in vivo clinical data from other studies, further validation of the computational models could also be performed in the same patient data set by measuring in vivo pre-operative and post-operative iFR, intraoperative TTFM, as well as post-operative myocardial perfusion imaging. The present study investigated 12 different TAG configurations using BIMA and RA for severe triple vessel coronary disease for each of five patients and is therefore the most comprehensive study of its kind in the literature, modeling 60 virtual grafting configurations (41-43). The computational methodology developed has the potential to investigate a wide range of other CABG configurations and hemodynamic scenarios. It is uncertain whether composite grafts can sustain adequate perfusion at both rest and hyperemia (44,45). Future models could also account for variations in flow and incorporate the effects of cardiopulmonary bypass, anesthesia, exercise, hyperemia, coronary autoregulation and increased graft diameters over time. They could also investigate patients with saphenous vein grafts, poorer ejection fractions, or those with LV hypertrophy and microvascular coronary disease where the use of anaortic configurations with composite Y-grafts has been clinically cautioned (44).
Patient-specific computational modeling provides important predictive hemodynamic flow information that cardiac surgeons can incorporate into their decision-making when planning graft configurations for individual patients. Surgeons take heed of computer model predictions as their current practice involves significant uncertainty regarding native coronary and bypass graft flows achieved following CABG. Based on these results, "COMCAB" has the potential to be a promising clinical decision-support software tool for personalized surgical planning for CABG configurations.
. /fcvm. . This is particularly important in avoiding situations of flow competition affecting bypass graft patency arising from the use of composite and sequential grafts such as those used more frequently in "anaortic" grafting techniques.

Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.

Ethics statement
The studies involving human participants were reviewed and approved by the Health and Disability Ethics Committee (HDEC) New Zealand (

Author contributions
KC and NS conceptualized the computational modeling. AP improved the code efficiency and visualization of the computational models. KC, PM, and NS conceptualized the study involving the cardiac surgeons. KC, AP, SW, PM, and NS edited the manuscript for intellectual content. All authors gave approval for the final version of the manuscript to be published.