A framework for different levels of integration of computational models into web-based virtual patients.

BACKGROUND: Virtual patients are increasingly common tools used in health care education to foster learning of clinical reasoning skills. One potential way to expand their functionality is to augment virtual patients' interactivity by enriching them with computational models of physiological and pathological processes. OBJECTIVE: The primary goal of this paper was to propose a conceptual framework for the integration of computational models within virtual patients, with particular focus on (1) characteristics to be addressed while preparing the integration, (2) the extent of the integration, (3) strategies to achieve integration, and (4) methods for evaluating the feasibility of integration. An additional goal was to pilot the first investigation of changing framework variables on altering perceptions of integration. METHODS: The framework was constructed using an iterative process informed by Soft System Methodology. The Virtual Physiological Human (VPH) initiative has been used as a source of new computational models. The technical challenges associated with development of virtual patients enhanced by computational models are discussed from the perspectives of a number of different stakeholders. Concrete design and evaluation steps are discussed in the context of an exemplar virtual patient employing the results of the VPH ARCH project, as well as improvements for future iterations. RESULTS: The proposed framework consists of four main elements. The first element is a list of feasibility features characterizing the integration process from three perspectives: the computational modelling researcher, the health care educationalist, and the virtual patient system developer. The second element included three integration levels: basic, where a single set of simulation outcomes is generated for specific nodes in the activity graph; intermediate, involving pre-generation of simulation datasets over a range of input parameters; advanced, including dynamic solution of the model. The third element is the description of four integration strategies, and the last element consisted of evaluation profiles specifying the relevant feasibility features and acceptance thresholds for specific purposes. The group of experts who evaluated the virtual patient exemplar found higher integration more interesting, but at the same time they were more concerned with the validity of the result. The observed differences were not statistically significant. CONCLUSIONS: This paper outlines a framework for the integration of computational models into virtual patients. The opportunities and challenges of model exploitation are discussed from a number of user perspectives, considering different levels of model integration. The long-term aim for future research is to isolate the most crucial factors in the framework and to determine their influence on the integration outcome.

Evaluation of models of sequestration flow in coronary arteries-Physiology versus anatomy?

BACKGROUND: Myocardial ischaemia results from insufficient coronary blood flow. Computed virtual fractional flow reserve (vFFR) allows quantification of proportional flow loss without the need for invasive pressure-wire testing. In the current study, we describe a novel, conductivity model of side branch flow, referred to as 'leak'. This leak model is a function of taper and local pressure, the latter of which may change radically when focal disease is present. This builds upon previous techniques, which either ignore side branch flow, or rely purely on anatomical factors. This study aimed to describe a new, conductivity model of side branch flow and compare this with established anatomical models. METHODS AND RESULTS: The novel technique was used to quantify vFFR, distal absolute flow (Qd) and microvascular resistance (CMVR) in 325 idealised 1D models of coronary arteries, modelled from invasive clinical data. Outputs were compared to an established anatomical model of flow. The conductivity model correlated and agreed with the reference model for vFFR (r = 0.895, p < 0.0001; ＋0.02, 95% CI 0.00 to ＋ 0.22), Qd (r = 0.959, p < 0.0001; -5.2 mL/min, 95% CI -52.2 to ＋13.0) and CMVR (r = 0.624, p < 0.0001; ＋50 Woods Units, 95% CI -325 to ＋2549). CONCLUSION: Agreement between the two techniques was closest for vFFR, with greater proportional differences seen for Qd and CMVR. The conductivity function assumes vessel taper was optimised for the healthy state and that CMVR was not affected by local disease. The latter may be addressed with further refinement of the technique or inferred from complementary image data. The conductivity technique may represent a refinement of current techniques for modelling coronary side-branch flow. Further work is needed to validate the technique against invasive clinical data.

Sex differences in coronary microvascular resistance measured by a computational fluid dynamics model.

Background: Increased coronary microvascular resistance (CMVR) is associated with coronary microvascular dysfunction (CMD). Although CMD is more common in women, sex-specific differences in CMVR have not been demonstrated previously. Aim: To compare CMVR between men and women being investigated for chest pain. Methods and results: We used a computational fluid dynamics (CFD) model of human coronary physiology to calculate absolute CMVR based on invasive coronary angiographic images and pressures in 203 coronary arteries from 144 individual patients. CMVR was significantly higher in women than men (860 [650-1,205] vs. 680 [520-865] WU, Z = -2.24, p = 0.025). None of the other major subgroup comparisons yielded any differences in CMVR. Conclusion: CMVR was significantly higher in women compared with men. These sex-specific differences may help to explain the increased prevalence of CMD in women.

In Vitro Modelling for Bulging Sinus Effects of an Expanded Polytetrafluoroethylene Valved Conduit Based on High-Speed 3D Leaflet Evaluation.

The study aimed to develop a pulmonary circulatory system capable of high-speed 3D reconstruction of valve leaflets to elucidate the local hemodynamic characteristics in the valved conduits with bulging sinuses. Then a simultaneous measurement system for leaflet structure and pressure and flow characteristics was designed to obtain valve leaflet dynamic behaviour with different conduit structures. An image preprocessing method was established to obtain the three leaflets behaviour simultaneously for one sequence with two leaflets images from each pair of three high-speed cameras. Firstly, the multi-digital image correlation analyses were performed, and then the valve leaflet structure was measured under the static condition with fixed opening angles in the water-filled visualization chamber and the pulsatile flow tests simulating paediatric pulmonary flow conditions in the different types of conduit structures; with or without bulging sinuses. The results showed the maximum 3D reconstruction error to be around 0.06 mm. In the steady flow test, the evaluation of opening angles under the different flow rates conditions was achieved. In the pulsatile flow test, each leaflet's opening and closing behaviours were successfully reconstructed simultaneously at the high-frequency recording rate of 960fps. Therefore, the system developed in this study confirms the design evaluation method of an ePTFE valved conduit behaviour with leaflet structures interacting with local fluid dynamics in the vicinity of valves. Clinical Relevance- The system reveals the bulging sinus effects on ePTFE valve leaflet motion by the 3D reconstruction using multi-camera high-speed sequential imaging in vitro.

The Complementary Value of Absolute Coronary Flow in the Assessment of Patients with Ischaemic Heart Disease (the COMPAC-Flow Study).

Fractional flow reserve (FFR) is the current gold-standard invasive assessment of coronary artery disease (CAD). FFR reports coronary blood flow (CBF) as a fraction of a hypothetical and unknown normal value. Although used routinely to diagnose CAD and guide treatment, how accurately FFR predicts actual CBF changes remains unknown. Here we compared fractional CBF with the absolute CBF (aCBF in mL/min), measured with a computational method during standard angiography and pressure-wire assessment, on 203 diseased arteries (143 patients). We found a substantial correlation between the two measurements (r 0.89, Cohen's Kappa 0.71). Concordance between fractional and absolute CBF reduction was high when FFR was >0.80 (91%), but reduced when FFR was ≤0.80 (81%), 0.70-0.80 (68%) and, particularly 0.75-0.80 (62%). Discordance was associated with coronary microvascular resistance, vessel diameter and mass of myocardium subtended, all factors to which FFR is agnostic. Assessment of aCBF complements FFR, and may be valuable to assess CBF, particularly in cases within the FFR 'grey-zone'.

The importance of three dimensional coronary artery reconstruction accuracy when computing virtual fractional flow reserve from invasive angiography.

Three dimensional (3D) coronary anatomy, reconstructed from coronary angiography (CA), is now being used as the basis to compute 'virtual' fractional flow reserve (vFFR), and thereby guide treatment decisions in patients with coronary artery disease (CAD). Reconstruction accuracy is therefore important. Yet the methods required remain poorly validated. Furthermore, the magnitude of vFFR error arising from reconstruction is unkown. We aimed to validate a method for 3D CA reconstruction and determine the effect this had upon the accuracy of vFFR. Clinically realistic coronary phantom models were created comprosing seven standard stenoses in aluminium and 15 patient-based 3D-printed, imaged with CA, three times, according to standard clinical protocols, yielding 66 datasets. Each was reconstructed using epipolar line projection and intersection. All reconstructions were compared against the real phantom models in terms of minimal lumen diameter, centreline and surface similarity. 3D-printed reconstructions (n = 45) and the reference files from which they were printed underwent vFFR computation, and the results were compared. The average error in reconstructing minimum lumen diameter (MLD) was 0.05 (± 0.03 mm) which was < 1% (95% CI 0.13-1.61%) compared with caliper measurement. Overall surface similarity was excellent (Hausdorff distance 0.65 mm). Errors in 3D CA reconstruction accounted for an error in vFFR of ± 0.06 (Bland Altman 95% limits of agreement). Errors arising from the epipolar line projection method used to reconstruct 3D coronary anatomy from CA are small but contribute to clinically relevant errors when used to compute vFFR.

A novel method for measuring absolute coronary blood flow and microvascular resistance in patients with ischaemic heart disease.

AIMS: Ischaemic heart disease is the reduction of myocardial blood flow, caused by epicardial and/or microvascular disease. Both are common and prognostically important conditions, with distinct guideline-indicated management. Fractional flow reserve (FFR) is the current gold-standard assessment of epicardial coronary disease but is only a surrogate of flow and only predicts percentage flow changes. It cannot assess absolute (volumetric) flow or microvascular disease. The aim of this study was to develop and validate a novel method that predicts absolute coronary blood flow and microvascular resistance (MVR) in the catheter laboratory. METHODS AND RESULTS: A computational fluid dynamics (CFD) model was used to predict absolute coronary flow (QCFD) and coronary MVR using data from routine invasive angiography and pressure-wire assessment. QCFD was validated in an in vitro flow circuit which incorporated patient-specific, three-dimensional printed coronary arteries; and then in vivo, in patients with coronary disease. In vitro, QCFD agreed closely with the experimental flow over all flow rates [bias +2.08 mL/min; 95% confidence interval (error range) -4.7 to +8.8 mL/min; R2 = 0.999, P < 0.001; variability coefficient <1%]. In vivo, QCFD and MVR were successfully computed in all 40 patients under baseline and hyperaemic conditions, from which coronary flow reserve (CFR) was also calculated. QCFD-derived CFR correlated closely with pressure-derived CFR (R2 = 0.92, P < 0.001). This novel method was significantly more accurate than Doppler-wire-derived flow both in vitro (±6.7 vs. ±34 mL/min) and in vivo (±0.9 vs. ±24.4 mmHg). CONCLUSIONS: Absolute coronary flow and MVR can be determined alongside FFR, in absolute units, during routine catheter laboratory assessment, without the need for additional catheters, wires or drug infusions. Using this novel method, epicardial and microvascular disease can be discriminated and quantified. This comprehensive coronary physiological assessment may enable a new level of patient stratification and management.

Modeling Approach for An Aortic Dissection with Endovascular Stenting.

Repair of dissected aorta requires remodeling the structure of the media. Modeling approaches specific to endovascular stenting for aortic dissection have been reported. We created a goat model of descending thoracic aortic dissection and reproduced its morphological characteristics in a mock circulatory system. The purpose of this study was to examine a newly developed aortic stent which was capable of installing to the aortic dissected lesion for biomedical hemodynamics point of view. In this study, we examined the changes in hemodynamics of dissected lesions and the amelioration by endovascular stent intervention. Firstly, we performed animal experiments with the dissected aorta and examined the effects of stenting on volumetric changes in the false lumen. Secondly, we made several types of 3-D stereolithographic dissected aortic models with silicone rubber membrane between the false and the true lumens. Then, the hemodynamic characteristics in each model were evaluated in the pulsatile flow conditions in a mock circulatory system. These modelling approaches enabled the quantitative examination of post-therapeutic effects of stenting followed by elucidating of hemodynamic changes in the vicinity of stents, which may follow the management of clinical amelioration of interventional treatment with aortic stenting.Clinical Relevance- This study represents a modelling approach of the dissected aorta for endovascular intervention using stenting followed by the examination of false lumen volumetric changes resulting in the deterioration of pressure increase in diseased lesions.

Measurement of in vitro cardiac deformation by means of 3D digital image correlation and ultrasound 2D speckle-tracking echocardiography.

Ultrasound-based 2D speckle-tracking echocardiography (US-2D-STE) is increasingly used to assess the functionality of the heart. In particular, the analysis of cardiac strain plays an important role in the identification of several cardiovascular diseases. However, this imaging technique presents some limitations associated with its operating principle that result in low accuracy and reproducibility of the measurement. In this study, an experimental framework for multimodal strain imaging in an in vitro porcine heart was developed. Specifically, the aim of this work was to analyse displacement and strain in the heart by means of 3D digital image correlation (3D-DIC) and US-2D-STE. Over a single cardiac cycle, displacement values obtained from the two techniques were in strong correlation, although systematically larger displacements were observed with 3D-DIC. Notwithstanding an absolute comparison of the strain measurements was not possible to achieve between the two methods, maximum principal strain directions computed with 3D-DIC were consistent with the longitudinal and circumferential strain distribution measured with US-2D-STE. 3D-DIC confirmed its high repeatability in quantifying displacement and strain over multiple cardiac cycles, unlike US-2D-STE which is affected by accumulated errors over time (i.e. drift). To conclude, this study demonstrates the potential of 3D-DIC to perform dynamic measurement of displacement and strain during heart deformations and supports future applications of this method in ex vivo beating heart platforms, which replicate more fully the complex contraction of the heart.

Full-field analysis of epicardial strain in an in vitro porcine heart platform.

The quantitative assessment of cardiac strain is increasingly performed to provide valuable insights on heart function. Currently, the most frequently used technique in the clinic is ultrasound-based speckle tracking echocardiography (STE). However, verification and validation of this modality are still under investigation and further reference measurements are required to support this activity. The aim of this work was to enable these reference measurements using a dynamic beating heart simulator to ensure reproducible, controlled, and realistic haemodynamic conditions and to validate the reliability of optical-based three-dimensional digital image correlation (3D-DIC) for a dynamic full-field analysis of epicardial strain. Specifically, performance assessment of 3D-DIC was carried out by evaluating the accuracy and repeatability of the strain measurements across multiple cardiac cycles in a single heart and between five hearts. Moreover, the ability of this optical method to differentiate strain variations when different haemodynamic conditions were imposed in the same heart was examined. Strain measurements were successfully accomplished in a region of the lateral left ventricle surface. Results were highly repeatable over heartbeats and across hearts (intraclass correlation coefficient = 0.99), whilst strain magnitude was significantly different between hearts, due to change in anatomy and wall thickness. Within an individual heart, strain variations between different haemodynamic scenarios were greater than the estimated error of the measurement technique. This study demonstrated the feasibility of applying 3D-DIC in a dynamic passive heart simulator. Most importantly, non-contact measurements were obtained at a high spatial resolution (~ 1.5 mm) allowing resolution of local variation of strain on the epicardial surface during ventricular filling. The experimental framework developed in this paper provides detailed measurement of cardiac strains under controlled conditions, as a reference for validation of clinical cardiac strain imaging modalities.

The role of venous valves in pressure shielding.

BACKGROUND: It is widely accepted that venous valves play an important role in reducing the pressure applied to the veins under dynamic load conditions, such as the act of standing up. This understanding is, however, qualitative and not quantitative. The purpose of this paper is to quantify the pressure shielding effect and its variation with a number of system parameters. METHODS: A one-dimensional mathematical model of a collapsible tube, with the facility to introduce valves at any position, was used. The model has been exercised to compute transient pressure and flow distributions along the vein under the action of an imposed gravity field (standing up). RESULTS: A quantitative evaluation of the effect of a valve, or valves, on the shielding of the vein from peak transient pressure effects was undertaken. The model used reported that a valve decreased the dynamic pressures applied to a vein when gravity is applied by a considerable amount. CONCLUSION: The model has the potential to increase understanding of dynamic physical effects in venous physiology, and ultimately might be used as part of an interventional planning tool.

Measurement of in vitro and in vivo stent geometry and deformation by means of 3D imaging and stereo-photogrammetry.

PURPOSE: To quantify variability of in vitro and in vivo measurement of 3D device geometry using 3D and biplanar imaging. METHODS: Comparison of stent reconstruction is reported for in vitro coronary stent deployment (using micro-CT and optical stereo-photogrammetry) and in vivo pulmonary valve stent deformation (using 4DCT and biplanar fluoroscopy). Coronary stent strut length and inter-strut angle were compared in the fully deployed configuration. Local (inter-strut angle) and global (dog-boning ratio) measures of stent deformation were reported during stent deployment. Pulmonary valve stent geometry was assessed throughout the cardiac cycle by reconstruction of stent geometry and measurement of stent diameter. RESULTS: Good agreement was obtained between methods for assessment of coronary stent geometry with maximum disagreement of +/- 0.03 mm (length) and +/- 3 degrees (angle). The stent underwent large, non-uniform, local deformations during balloon inflation, which did not always correlate with changes in stent diameter. Three-dimensional reconstruction of the pulmonary valve stent was feasible for all frames of the fluoroscopy and for 4DCT images, with good correlation between the diameters calculated from the two methods. The largest compression of the stent during the cardiac cycle was 6.98% measured from fluoroscopy and 7.92% from 4DCT, both in the most distal ring. CONCLUSIONS: Quantitative assessment of stent geometry reconstructed from biplanar imaging methods in vitro and in vivo has shown good agreement with geometry reconstructed from 3D techniques. As a result of their short image acquisition time, biplanar methods may have significant advantages in the measurement of dynamic 3D stent deformation.

Uncertainty assessment of imaging techniques for the 3D reconstruction of stent geometry.

This paper presents a quantitative assessment of uncertainty for the 3D reconstruction of stents. This study investigates a CP stent (Numed, USA) used in congenital heart disease applications with a focus on the variance in measurements of stent geometry. The stent was mounted on a model of patient implantation site geometry, reconstructed from magnetic resonance images, and imaged using micro-computed tomography (CT), conventional CT, biplane fluoroscopy and optical stereo-photogrammetry. Image data were post-processed to retrieve the 3D stent geometry. Stent strut length, separation angle and cell asymmetry were derived and repeatability was assessed for each technique along with variation in relation to µCT data, assumed to represent the gold standard. The results demonstrate the performance of biplanar reconstruction methods is comparable with volumetric CT scans in evaluating 3D stent geometry. Uncertainty on the evaluation of strut length, separation angle and cell asymmetry using biplanar fluoroscopy is of the order ±0.2mm, 3° and 0.03, respectively. These results support the use of biplanar fluoroscopy for in vivo measurement of 3D stent geometry and provide quantitative assessment of uncertainty in the measurement of geometric parameters.

From histology and imaging data to models for in-stent restenosis.

The implantation of stents has been used to treat coronary artery stenosis for several decades. Although stenting is successful in restoring the vessel lumen and is a minimally invasive approach, the long-term outcomes are often compromised by in-stent restenosis (ISR). Animal models have provided insights into the pathophysiology of ISR and are widely used to evaluate candidate drug inhibitors of ISR. Such biological models allow the response of the vessel to stent implantation to be studied without the variation of lesion characteristics encountered in patient studies.This paper describes the development of complementary in silico models employed to improve the understanding of the biological response to stenting using a porcine model of restenosis. This includes experimental quantification using microCT imaging and histology and the use of this data to establish numerical models of restenosis. Comparison of in silico results with histology is used to examine the relationship between spatial localization of fluid and solid mechanics stimuli immediately post-stenting. Multi-scale simulation methods are employed to study the evolution of neointimal growth over time and the variation in the extent of neointimal hyperplasia within the stented region. Interpretation of model results through direct comparison with the biological response contributes to more detailed understanding of the pathophysiology of ISR, and suggests the focus for follow-up studies.In conclusion we outline the challenges which remain to both complete our understanding of the mechanisms responsible for restenosis and translate these models to applications in stent design and treatment planning at both population-based and patient-specific levels.

Balloon folding affects the symmetry of stent deployment: experimental and computational evidence.

The level of restenosis following coronary artery stenting may be related to the deployed stent geometry. This study investigated the influence of two balloon folding patterns (;C' and ;S' shaped) on stent deployment. In vitro stent expansion showed ;S' shape folding produced more uniform expansion than ;C' shape folding. A numerical contact model (NCM) was developed to study the detail of load transfer between balloon and stent. Finite element analysis of the Palmaz-Schatz 204C stent provided a composite non-linear material model for the NCM. Agreement between the predicted final stent geometry and experimental results was strongly dependent on the frictional coefficient between the stent and balloon. We conclude that non-uniform contact may contribute to the asymmetry of deployed stents reported clinically.

Influence of intermittent compression cuff design on calf deformation: computational results.

The intermittent compression of the calf with an external pressure cuff for the prevention of deep vein thrombosis (DVT) is a well established treatment for surgical patients. The exact mechanisms by which DVT is prevented are poorly understood. This study presents a finite element model of calf cross section, based on MR images of calf geometry, to examine the variation in calf deformation during compression with four different cuff types. Cuff pressure distribution is modelled using interface pressures obtained in a volunteer study. The model has been validated against gross calf deformation obtained from MR images of the compressed calf. This validation has illustrated the importance of out-of-plane boundary conditions, material properties and the variation in cuff loading in the axial direction. In the future this model may have merit in determining optimum pressure loading regimes for Intermittent Pneumatic Compression (IPC) cuff design.

Influence of intermittent compression cuff design on interface pressure and calf deformation: experimental results.

Intermittent pneumatic compression (IPC) is widely used for deep vein thrombosis (DVT) prophylaxis. The technique involves periodic inflation of a compression cuff around a limb, which acts to simulate the muscle pump mechanism, encouraging venous blood flow. However, there is uncertainty regarding the relationship between compression, vascular effects and clinical outcomes. This study investigates calf compression provided by four IPC cuffs with different air bladder configurations. Interface pressure between the cuff and the skin surface is measured and magnetic resonance (MR) images are obtained showing the calf cross section before and during compression. The data will be used to inform numerical simulations of IPC, leading to increased understanding of the implications of cuff design in relation to IPC and DVT prophylaxis.