Systematic review and meta-analysis of Murray's law in the coronary arterial circulation.

Murray's law has been viewed as a fundamental law of physiology. Relating blood flow ([Formula: see text]) to vessel diameter (D) ([Formula: see text]·âˆ·D3), it dictates minimum lumen area (MLA) targets for coronary bifurcation percutaneous coronary intervention (PCI). The cubic exponent (3.0), however, has long been disputed, with alternative theoretical derivations, arguing this should be closer to 2.33 (7/3). The aim of this meta-analysis was to quantify the optimum flow-diameter exponent in human and mammalian coronary arteries. We conducted a systematic review and meta-analysis of all articles quantifying an optimum flow-diameter exponent for mammalian coronary arteries within the Cochrane library, PubMed Medline, Scopus, and Embase databases on 20 March 2023. A random-effects meta-analysis was used to determine a pooled flow-diameter exponent. Risk of bias was assessed with the National Institutes of Health (NIH) quality assessment tool, funnel plots, and Egger regression. From a total of 4,772 articles, 18 were suitable for meta-analysis. Studies included data from 1,070 unique coronary trees, taken from 372 humans and 112 animals. The pooled flow diameter exponent across both epicardial and transmural arteries was 2.39 (95% confidence interval: 2.24-2.54; I2 = 99%). The pooled exponent of 2.39 showed very close agreement with the theoretical exponent of 2.33 (7/3) reported by Kassab and colleagues. This exponent may provide a more accurate description of coronary morphometric scaling in human and mammalian coronary arteries, as compared with Murray's original law. This has important implications for the assessment, diagnosis, and interventional treatment of coronary artery disease.

Validation of a novel numerical model to predict regionalized blood flow in the coronary arteries.

Aims: Ischaemic heart disease results from insufficient coronary blood flow. Direct measurement of absolute flow (mL/min) is feasible, but has not entered routine clinical practice in most catheterization laboratories. Interventional cardiologists, therefore, rely on surrogate markers of flow. Recently, we described a computational fluid dynamics (CFD) method for predicting flow that differentiates inlet, side branch, and outlet flows during angiography. In the current study, we evaluate a new method that regionalizes flow along the length of the artery. Methods and results: Three-dimensional coronary anatomy was reconstructed from angiograms from 20 patients with chronic coronary syndrome. All flows were computed using CFD by applying the pressure gradient to the reconstructed geometry. Side branch flow was modelled as a porous wall boundary. Side branch flow magnitude was based on morphometric scaling laws with two models: a homogeneous model with flow loss along the entire arterial length; and a regionalized model with flow proportional to local taper. Flow results were validated against invasive measurements of flow by continuous infusion thermodilution (Coroventis™, Abbott). Both methods quantified flow relative to the invasive measures: homogeneous (r 0.47, P 0.006; zero bias; 95% CI -168 to +168 mL/min); regionalized method (r 0.43, P 0.013; zero bias; 95% CI -175 to +175 mL/min). Conclusion: During angiography and pressure wire assessment, coronary flow can now be regionalized and differentiated at the inlet, outlet, and side branches. The effect of epicardial disease on agreement suggests the model may be best targeted at cases with a stenosis close to side branches.

The EurValve model execution environment.

The goal of this paper is to present a dedicated high-performance computing (HPC) infrastructure which is used in the development of a so-called reduced-order model (ROM) for simulating the outcomes of interventional procedures which are contemplated in the treatment of valvular heart conditions. Following a brief introduction to the problem, the paper presents the design of a model execution environment, in which representative cases can be simulated and the parameters of the ROM fine-tuned to enable subsequent deployment of a decision support system without further need for HPC. The presentation of the system is followed by information concerning its use in processing specific patient cases in the context of the EurValve international collaboration.

Estimation of valvular resistance of segmented aortic valves using computational fluid dynamics.

Aortic valve stenosis is associated with an elevated left ventricular pressure and transaortic pressure drop. Clinicians routinely use Doppler ultrasound to quantify aortic valve stenosis severity by estimating this pressure drop from blood velocity. However, this method approximates the peak pressure drop, and is unable to quantify the partial pressure recovery distal to the valve. As pressure drops are flow dependent, it remains difficult to assess the true significance of a stenosis for low-flow low-gradient patients. Recent advances in segmentation techniques enable patient-specific Computational Fluid Dynamics (CFD) simulations of flow through the aortic valve. In this work a simulation framework is presented and used to analyze data of 18 patients. The ventricle and valve are reconstructed from 4D Computed Tomography imaging data. Ventricular motion is extracted from the medical images and used to model ventricular contraction and corresponding blood flow through the valve. Simplifications of the framework are assessed by introducing two simplified CFD models: a truncated time-dependent and a steady-state model. Model simplifications are justified for cases where the simulated pressure drop is above 10 mmHg. Furthermore, we propose a valve resistance index to quantify stenosis severity from simulation results. This index is compared to established metrics for clinical decision making, i.e. blood velocity and valve area. It is found that velocity measurements alone do not adequately reflect stenosis severity. This work demonstrates that combining 4D imaging data and CFD has the potential to provide a physiologically relevant diagnostic metric to quantify aortic valve stenosis severity.

A cell-based mechanical model of coronary artery tunica media.

A three-dimensional cell-based mechanical model of coronary artery tunica media is proposed. The model is composed of spherical cells forming a hexagonal close-packed lattice. Tissue anisotropy is taken into account by varying interaction forces with the direction of intercellular connection. Several cell-centre interaction potentials for repulsion and attraction are considered, including the Hertz contact model and its neo-Hookean extension, the Johnson-Kendall-Roberts model of adhesive contact, and a wormlike chain model. The model is validated against data from in vitro uni-axial tension tests performed on dissected strips of tunica media. The wormlike chain potential in combination with the neo-Hookean Hertz contact model produces stress-stretch curves which represent the experimental data very well.

MRI model-based non-invasive differential diagnosis in pulmonary hypertension.

Pulmonary hypertension(PH) is a disorder characterised by increased mean pulmonary arterial pressure. Currently, the diagnosis of PH relies upon measurements taken during invasive right heart catheterisation (RHC). This paper describes a process to derive diagnostic parameters using only non-invasive methods based upon MRI imaging alone. Simultaneous measurements of main pulmonary artery (MPA) anatomy and flow are interpreted by 0D and 1D mathematical models, in order to infer the physiological status of the pulmonary circulation. Results are reported for 35 subjects, 27 of whom were patients clinically investigated for PH and eight of whom were healthy volunteers. The patients were divided into 3 sub-groups according to the severity of the disease state, one of which represented a negative diagnosis (NoPH), depending on the results of the clinical investigation, which included RHC and complementary MR imaging. Diagnostic indices are derived from two independent mathematical models, one based on the 1D wave equation and one based on an RCR Windkessel model. Using the first model it is shown that there is an increase in the ratio of the power in the reflected wave to that in the incident wave (Wpb/Wptotal) according to the classification of the disease state. Similarly, the second model shows an increase in the distal resistance with the disease status. The results of this pilot study demonstrate that there are statistically significant differences in the parameters derived from the proposed models depending on disease status, and thus suggest the potential for development of a non-invasive, image-based diagnostic test for pulmonary hypertension.

Measuring three-dimensional strain distribution in tendon.

Tendons are tough fibrous tissues that facilitate skeletal movement by transferring muscular force to bone. Studies into the effects of mechanical stress on tendons have shown that these can either accelerate healing or cause tendon injuries depending on the load applied. It is known that local strain magnitude and direction play an important role in tendon remodelling and also failure, and different techniques to study strain distribution have been proposed. Image registration and processing techniques are among the recently employed methods. In this study, a novel three-dimensional image processing technique using the Sheffield Image Registration Toolkit is introduced to study local strain and displacement distribution in tendon. The results show that the local normal strain values in the loading axis are smaller than the global applied load, and fibre sliding was detected as a dominant mechanism for transferring the applied load within tendon. However, results from different samples suggest three distinct modes of deformation during loading, as some show only parallel sliding of fibres in respect to the loading axis, whereas others are twisted or deflected in directions transverse to the loading axis. The proposed 3D image registration method is essential for analysing this out-of-plane movement, which cannot be detected using a standard 2D method.

Generating anatomical models of the heart and the aorta from medical images for personalized physiological simulations.

The anatomy and motion of the heart and the aorta are essential for patient-specific simulations of cardiac electrophysiology, wall mechanics and hemodynamics. Within the European integrated project euHeart, algorithms have been developed that allow to efficiently generate patient-specific anatomical models from medical images from multiple imaging modalities. These models, for instance, account for myocardial deformation, cardiac wall motion, and patient-specific tissue information like myocardial scar location. Furthermore, integration of algorithms for anatomy extraction and physiological simulations has been brought forward. Physiological simulations are linked closer to anatomical models by encoding tissue properties, like the muscle fibers, into segmentation meshes. Biophysical constraints are also utilized in combination with image analysis to assess tissue properties. Both examples show directions of how physiological simulations could provide new challenges and stimuli for image analysis research in the future.

Multi-scale interaction of particulate flow and the artery wall.

We discuss, from the perspective of basic science, the physical and biological processes which underlie atherosclerotic (plaque) initiation at the vascular endothelium, identifying the widely separated spatial and temporal scales which participate. We draw on current, related models of vessel wall evolution, paying particular attention to the role of particulate flow (blood is not a continuum fluid), and proceed to propose, then validate all the key components in a multiply-coupled, multi-scale modeling strategy (in qualitative terms only, note). Eventually, this strategy should lead to a quantitative, patient-specific understanding of the coupling between particulate flow and the endothelial state.

@neurIST complex information processing toolchain for the integrated management of cerebral aneurysms.

Cerebral aneurysms are a multi-factorial disease with severe consequences. A core part of the European project @neurIST was the physical characterization of aneurysms to find candidate risk factors associated with aneurysm rupture. The project investigated measures based on morphological, haemodynamic and aneurysm wall structure analyses for more than 300 cases of ruptured and unruptured aneurysms, extracting descriptors suitable for statistical studies. This paper deals with the unique challenges associated with this task, and the implemented solutions. The consistency of results required by the subsequent statistical analyses, given the heterogeneous image data sources and multiple human operators, was met by a highly automated toolchain combined with training. A testimonial of the successful automation is the positive evaluation of the toolchain by over 260 clinicians during various hands-on workshops. The specification of the analyses required thorough investigations of modelling and processing choices, discussed in a detailed analysis protocol. Finally, an abstract data model governing the management of the simulation-related data provides a framework for data provenance and supports future use of data and toolchain. This is achieved by enabling the easy modification of the modelling approaches and solution details through abstract problem descriptions, removing the need of repetition of manual processing work.

Adaptation and development of software simulation methodologies for cardiovascular engineering: present and future challenges from an end-user perspective.

This paper describes the use of diverse software tools in cardiovascular applications. These tools were primarily developed in the field of engineering and the applications presented push the boundaries of the software to address events related to venous and arterial valve closure, exploration of dynamic boundary conditions or the inclusion of multi-scale boundary conditions from protein to organ levels. The future of cardiovascular research and the challenges that modellers and clinicians face from validation to clinical uptake are discussed from an end-user perspective.

Automated tracking of migrating cells in phase-contrast video microscopy sequences using image registration.

Analysis of in vitro cell motility is a useful tool for assessing cellular response to a range of factors. However, the majority of cell-tracking systems available are designed primarily for use with fluorescently labelled images. In this paper, five commonly used tracking systems are examined for their performance compared with the use of a novel in-house cell-tracking system based on the principles of image registration and optical flow. Image registration is a tool commonly used in medical imaging to correct for the effects of patient motion during imaging procedures and works well on low-contrast images, such as those found in bright-field and phase-contrast microscopy. The five cell-tracking systems examined were Retrac, a manual tracking system used as the gold standard; CellTrack, a recently released freely downloadable software system that uses a combination of tracking methods; ImageJ, which is a freely available piece of software with a plug-in for automated tracking (MTrack2) and Imaris and Volocity, both commercially available automated tracking systems. All systems were used to track migration of human epithelial cells over ten frames of a phase-contrast time-lapse microscopy sequence. This showed that the in-house image-registration system was the most effective of those tested when tracking non-dividing epithelial cells in low-contrast images, with a successful tracking rate of 95%. The performance of the tracking systems was also evaluated by tracking fluorescently labelled epithelial cells imaged with both phase-contrast and confocal microscopy techniques. The results showed that using fluorescence microscopy instead of phase contrast does improve the tracking efficiency for each of the tested systems. For the in-house software, this improvement was relatively small (<5% difference in tracking success rate), whereas much greater improvements in performance were seen when using fluorescence microscopy with Volocity and ImageJ.

A validated model of calf compression and deep vessel collapse during external cuff inflation.

This paper presents a validated model of calf compression with an external pressure cuff as used for deep vein thrombosis. Magnetic resonance (MR) images of calf geometry were used to generate subject-specific finite-element (FE) models of the calf cross section. Ultrasound images of deep vessel collapse obtained through a water-filled cuff were used to validate model behavior. Calf/cuff pressure interface measurements were applied to the FE model and the resulting tissue deformation was compared with MR image in normal volunteers (three females, four males, age range 20-55) using two distinct cuffs. MR observations and the model results showed good qualitative agreement. A similar reduction in cross-sectional area of the posterior tibial veins was obtained under both symmetric compression (89%) and asymmetric compression (81%), but greater compression of the anterior tibial veins was achieved with symmetric compression. The need to account for the effective compressibility of the calf tissue suggests that external measurements of the calf tissue deformation will not accurately predict deep vessel collapse. These results have implications for the modification of venous haemodynamics by such systems and could help to improve cuff design.

Modelling the interaction of haemodynamics and the artery wall: current status and future prospects.

Clinical research has historically focused on the two main strategies of in vivo and in vitro experimentation. The concept of applying scientific theory to direct clinical applications is relatively recent. In this paper we focus on the interaction of wall shear stress with the endothelium and discuss how 'state of the art' computer modelling techniques can provide valuable data to aid understanding. Such data may be used to inform experiment and further, may help identify the key features of this complex system. Current emphasis is on coupling haemodynamics with models of biological phenomena to test hypotheses or predict the likely outcome of a disease or an intervention. New technologies to enable the integration of models of different types, levels of complexity and scales, are being developed. As will be discussed, the ultimate goal is the translation of this technology to the clinical arena.

The application of multiscale modelling to the process of development and prevention of stenosis in a stented coronary artery.

The inherent complexity of biomedical systems is well recognized; they are multiscale, multiscience systems, bridging a wide range of temporal and spatial scales. While the importance of multiscale modelling in this context is increasingly recognized, there is little underpinning literature on the methodology and generic description of the process. The COAST (complex autonoma simulation technique) project aims to address this by developing a multiscale, multiscience framework, coined complex autonoma (CxA), based on a hierarchical aggregation of coupled cellular automata (CA) and agent-based models (ABMs). The key tenet of COAST is that a multiscale system can be decomposed into N single-scale CA or ABMs that mutually interact across the scales. Decomposition is facilitated by building a scale separation map on which each single-scale system is represented according to its spatial and temporal characteristics. Processes having well-separated scales are thus easily identified as the components of the multiscale model. This paper focuses on methodology, introduces the concept of the CxA and demonstrates its use in the generation of a multiscale model of the physical and biological processes implicated in a challenging and clinically relevant problem, namely coronary artery in-stent restenosis.

The EuroPhysiome, STEP and a roadmap for the virtual physiological human.

Biomedical science and its allied disciplines are entering a new era in which computational methods and technologies are poised to play a prevalent role in supporting collaborative investigation of the human body. Within Europe, this has its focus in the virtual physiological human (VPH), which is an evolving entity that has emerged from the EuroPhysiome initiative and the strategy for the EuroPhysiome (STEP) consortium. The VPH is intended to be a solution to common infrastructure needs for physiome projects across the globe, providing a unifying architecture that facilitates integration and prediction, ultimately creating a framework capable of describing Homo sapiens in silico. The routine reliance of the biomedical industry, biomedical research and clinical practice on information technology (IT) highlights the importance of a tailor-made and robust IT infrastructure, but numerous challenges need to be addressed if the VPH is to become a mature technological reality. Appropriate investment will reap considerable rewards, since it is anticipated that the VPH will influence all sectors of society, with implications predominantly for improved healthcare, improved competitiveness in industry and greater understanding of (patho)physiological processes. This paper considers issues pertinent to the development of the VPH, highlighted by the work of the STEP consortium.

Reproducibility of haemodynamical simulations in a subject-specific stented aneurysm model--a report on the Virtual Intracranial Stenting Challenge 2007.

This paper presents the results of the Virtual Intracranial Stenting Challenge (VISC) 2007, an international initiative whose aim was to establish the reproducibility of state-of-the-art haemodynamical simulation techniques in subject-specific stented models of intracranial aneurysms (IAs). IAs are pathological dilatations of the cerebral artery walls, which are associated with high mortality and morbidity rates due to subarachnoid haemorrhage following rupture. The deployment of a stent as flow diverter has recently been indicated as a promising treatment option, which has the potential to protect the aneurysm by reducing the action of haemodynamical forces and facilitating aneurysm thrombosis. The direct assessment of changes in aneurysm haemodynamics after stent deployment is hampered by limitations in existing imaging techniques and currently requires resorting to numerical simulations. Numerical simulations also have the potential to assist in the personalized selection of an optimal stent design prior to intervention. However, from the current literature it is difficult to assess the level of technological advancement and the reproducibility of haemodynamical predictions in stented patient-specific models. The VISC 2007 initiative engaged in the development of a multicentre-controlled benchmark to analyse differences induced by diverse grid generation and computational fluid dynamics (CFD) technologies. The challenge also represented an opportunity to provide a survey of available technologies currently adopted by international teams from both academic and industrial institutions for constructing computational models of stented aneurysms. The results demonstrate the ability of current strategies in consistently quantifying the performance of three commercial intracranial stents, and contribute to reinforce the confidence in haemodynamical simulation, thus taking a step forward towards the introduction of simulation tools to support diagnostics and interventional planning.

Efficient computational fluid dynamics mesh generation by image registration.

Most implementations of computational fluid dynamics (CFD) solutions require a discretisation or meshing of the solution domain. The production from a medical image of a computationally efficient mesh representing the structures of interest can be time consuming and labour-intensive, and remains a major bottleneck in the clinical application of CFD. This paper presents a method for deriving a patient-specific mesh from a medical image. The method uses volumetric registration of a pseudo-image, produced from an idealised template mesh, with the medical image. The registration algorithm used is robust and computationally efficient. The accuracy of the new algorithm is measured in terms of the distance between a registered surface and a known surface, for image data derived from casts of the lumen of two different vessels. The true surface is identified by laser profiling. The average distance between the surface points measured by the laser profiler and the surface of the mapped mesh is better than 0.2 mm. For the images analysed, the new algorithm is shown to be 2-3 times more accurate than a standard published algorithm based on maximising normalised mutual information. Computation times are approximately 18 times faster for the new algorithm than the standard algorithm. Examples of the use of the algorithm on two clinical examples are also given. The registration methodology lends itself immediately to the construction of dynamic mesh models in which vessel wall motion is obtained directly using registration.

Application and validation of the lattice Boltzmann method for modelling flow-related clotting.

The purpose of this paper is to present a simple clotting model, based on residence time and shear stress distribution, that can simulate the deposition over time of enzyme-activated milk in an in vitro system. Results for the model are compared with experiments exhibiting clot deposition in the region of a sharp-edged stenosis. The milk experiments have been shown to be a valuable analogue for the experimental representation of flow-induced blood clotting, particularly in the context of separation of hydrodynamic from biochemical factors. The facility to predict the flow-induced clotting of the blood analogue, in which the chemistry reduces to what is effectively a zeroth order reaction, gives confidence in this physics-based approach to simulation of the final part of the coagulation cascade. This type of study is a necessary precursor to the development of a complex, multi-factorial, biochemical model of the process of thrombosis. In addition to the clotting simulations, comparisons are reported between the computed flow patterns prior to clot deposition and flow visualisation studies. Excellent agreement of hydrodynamic parameters is reported for a Reynolds number of 100, and qualitative agreement is seen for the complex, disturbed flow occurring at a physiologically relevant Reynolds number of 550. The explicit, time-stepping lattice Boltzmann approach may have particular merit for the transitional flow at this higher Reynolds number.

Characterisation of the haemodynamics of the superior mesenteric artery.

In contrast to its prevalence in the surrounding vasculature, occurrence of primary atherosclerotic disease in the superior mesenteric artery (SMA) is rare (Glagov et al., 1988. Hemodynamics and atherosclerosis, Insights and perspectives gained from studies of human arteries. Archives of Pathology and Laboratory Medicine 112(10), 1018-1031; Hansen et al., 2004. Mesenteric artery disease in the elderly. Journal of Vascular Surgery 40(1), 45-52). We hypothesise that this sparing might be attributed to more favourable haemodynamic characteristics in the SMA than in other vessels locally. Dynamic magnetic resonance imaging (MRI) images established that the SMA is highly mobile (Jeays, 2006. Investigation of blood flow in the superior mesenteric artery and its potential influence on atheroma and gut ischaemia. Ph.D. Thesis, University of Sheffield), and thus that an analysis based on rigid geometry might be inappropriate. This paper describes an efficient methodology for the construction of a patient-specific, time-dependent model of an arterial segment and reports the results of a haemodynamic characterisation of the SMA for one individual. A transient computational fluid dynamic (CFD) model was constructed by morphing a parametric mesh constructed from simple geometric primitives. This process has the merit that it is easy to control the element size distribution mapped onto the original geometric primitives. It is robust in operation, and is ideally suited to the generation of dynamic CFD meshes of arterial systems that are free from major pathology. Flow boundary conditions were determined based on phase contrast MRI velocity measurements. Comparative studies with rigid walls and with moving walls, based on the transient data, indicated that, despite the significant motion of the SMA (radial dilation of the order of 10% and translation of the order of the radius), the maximum (spatially and temporally-resolved) wall shear stresses changed by no more than 21.6% of a global norm, and the average change was less than 2.1%.