Method and Apparatus for Soft Tissue Material Parameter Estimation Using Tissue Tagged Magnetic Resonance Imaging

Validating MRI-Derived Myocardial Stiffness Estimates Using In Vitro Synthetic Heart Models

Impaired cardiac filling in response to increased passive myocardial stiffness contributes to the pathophysiology of heart failure. By leveraging cardiac MRI data and ventricular pressure measurements, we can estimate in vivo passive myocardial stiffness using personalized inverse finite element models. While it is well-known that this approach is subject to uncertainties, only few studies quantify the accuracy of these stiffness estimates. This lack of validation is, at least in part, due to the absence of ground truth in vivo passive myocardial stiffness values. Here, using 3D printing, we created soft, homogenous, isotropic, hyperelastic heart phantoms of varying geometry and stiffness and simulate diastolic filling by incorporating the phantoms into an MRI-compatible left ventricular inflation system. We estimate phantom stiffness from MRI and pressure data using inverse finite element analyses based on a Neo-Hookean model. We demonstrate that our identified softest and stiffest values of 215.7 and 512.3 kPa agree well with the ground truth of 226.2 and 526.4 kPa. Overall, our estimated stiffnesses revealed a good agreement with the ground truth ([Formula: see text] error) across all models. Our results suggest that MRI-driven computational constitutive modeling can accurately estimate synthetic heart material stiffnesses in the range of 200-500 kPa.

Robust and efficient fixed-point algorithm for the inverse elastostatic problem to identify myocardial passive material parameters and the unloaded reference configuration

Image-based computational models of the heart represent a powerful tool to shed new light on the mechanisms underlying physiological and pathological conditions in cardiac function and to improve diagnosis and therapy planning. However, in order to enable the clinical translation of such models, it is crucial to develop personalized models that are able to reproduce the physiological reality of a given patient. There have been numerous contributions in experimental and computational biomechanics to characterize the passive behavior of the myocardium. However, most of these studies suffer from severe limitations and are not applicable to high-resolution geometries. In this work, we present a novel methodology to perform an automated identification of in vivo properties of passive cardiac biomechanics. The highly-efficient algorithm fits material parameters against the shape of a patient-specific approximation of the end-diastolic pressure-volume relation (EDPVR). Simultaneously, an unloaded reference configuration is generated, where a novel line search strategy to improve convergence and robustness is implemented. Only clinical image data or previously generated meshes at one time point during diastole and one measured data point of the EDPVR are required as an input. The proposed method can be straightforwardly coupled to existing finite element (FE) software packages and is applicable to different constitutive laws and FE formulations. Sensitivity analysis demonstrates that the algorithm is robust with respect to initial input parameters.

An Implementation of Patient-Specific Biventricular Mechanics Simulations With a Deep Learning and Computational Pipeline

Parameterised patient-specific models of the heart enable quantitative analysis of cardiac function as well as estimation of regional stress and intrinsic tissue stiffness. However, the development of personalised models and subsequent simulations have often required lengthy manual setup, from image labelling through to generating the finite element model and assigning boundary conditions. Recently, rapid patient-specific finite element modelling has been made possible through the use of machine learning techniques. In this paper, utilising multiple neural networks for image labelling and detection of valve landmarks, together with streamlined data integration, a pipeline for generating patient-specific biventricular models is applied to clinically-acquired data from a diverse cohort of individuals, including hypertrophic and dilated cardiomyopathy patients and healthy volunteers. Valve motion from tracked landmarks as well as cavity volumes measured from labelled images are used to drive realistic motion and estimate passive tissue stiffness values. The neural networks are shown to accurately label cardiac regions and features for these diverse morphologies. Furthermore, differences in global intrinsic parameters, such as tissue anisotropy and normalised active tension, between groups illustrate respective underlying changes in tissue composition and/or structure as a result of pathology. This study shows the successful application of a generic pipeline for biventricular modelling, incorporating artificial intelligence solutions, within a diverse cohort.

Surface deformation tracking and modelling of soft materials

Many computer vision algorithms have been presented to track surface deformations, but few have provided a direct comparison of measurements with other stereoscopic approaches and physics-based models. We have previously developed a phase-based cross-correlation algorithm to track dense distributions of displacements over three-dimensional surfaces. In the present work, we compare this algorithm with one that uses an independent tracking system, derived from an array of fluorescent microspheres. A smooth bicubic Hermite mesh was fitted to deformations obtained from the phase-based cross-correlation data. This mesh was then used to estimate the microsphere locations, which were compared to stereo reconstructions of the microsphere positions. The method was applied to a 35 mm × 35 mm × 35 mm soft silicone gel cube under indentation, with three square bands of microspheres placed around the indenter tip. At an indentation depth of 4.5 mm, the root-mean-square (RMS) differences between the reconstructed positions of the microspheres and their identified positions for the inner, middle, and outer bands were 60 µm, 20 µm, and 19 µm, respectively. The usefulness of the strain-tracking data for physics-based finite element modelling of large deformation mechanics was then demonstrated by estimating a neo-Hookean stiffness parameter for the gel. At the optimal constitutive parameter estimate, the RMS difference between the measured microsphere positions and their finite element model-predicted locations was 143 µm.

Cardiovascular magnetic resonance elastography: A review

Cardiovascular diseases are the leading cause of death worldwide. These cardiovascular diseases are associated with mechanical changes in the myocardium and aorta. It is known that stiffness is altered in many diseases, including the spectrum of ischemia, diastolic dysfunction, hypertension and hypertrophic cardiomyopathy. In addition, the stiffness of the aortic wall is altered in multiple diseases, including hypertension, coronary artery disease and aortic aneurysm formation. For example, in diastolic dysfunction in which the ejection fraction is preserved, stiffness can potentially be an important biomarker. Similarly, in aortic aneurysms, stiffness can provide valuable information with regard to rupture potential. A number of studies have addressed invasive and non-invasive approaches to test and measure the mechanical properties of the myocardium and aorta. One of the non-invasive approaches is magnetic resonance elastography (MRE). MRE is a phase-contrast magnetic resonance imaging technique that measures tissue stiffness non-invasively. This review article highlights the technical details and application of MRE in the quantification of myocardial and aortic stiffness in different disease states.

Left Ventricular Diastolic Myocardial Stiffness and End-Diastolic Myofibre Stress in Human Heart Failure Using Personalised Biomechanical Analysis

Understanding the aetiology of heart failure with preserved (HFpEF) and reduced (HFrEF) ejection fraction requires knowledge of biomechanical factors such as diastolic myocardial stiffness and stress. Cine CMR images and intra-ventricular pressure recordings were acquired in 8 HFrEF, 11 HFpEF and 5 control subjects. Diastolic myocardial stiffness was estimated using biomechanical models and found to be greater in HFrEF (6.4 ± 1.2 kPa) than HFpEF (2.7 ± 0.6 kPa, p < 0.05) and also greater than control (1.2 ± 0.4 kPa, p < 0.005). End-diastolic mid-ventricular myofibre stress derived from the personalised biomechanics model was higher in HFrEF (2.9 ± 0.3 kPa) than control (0.9 ± 0.3 kPa, p < 0.01). Chamber stiffness, measured from the slope of the diastolic pressure-volume relationship, is determined by the intrinsic tissue properties as well as the size and shape of the heart, and was unable to distinguish between any of the three groups (p > 0.05). Personalised biomechanical analysis may provide more specific information about myocardial mechanical behaviour than global chamber indices, which are confounded by variations in ventricular geometry.

High-resolution data assimilation of cardiac mechanics applied to a dyssynchronous ventricle

Computational models of cardiac mechanics, personalized to a patient, offer access to mechanical information above and beyond direct medical imaging. Additionally, such models can be used to optimize and plan therapies in-silico, thereby reducing risks and improving patient outcome. Model personalization has traditionally been achieved by data assimilation, which is the tuning or optimization of model parameters to match patient observations. Current data assimilation procedures for cardiac mechanics are limited in their ability to efficiently handle high-dimensional parameters. This restricts parameter spatial resolution, and thereby the ability of a personalized model to account for heterogeneities that are often present in a diseased or injured heart. In this paper, we address this limitation by proposing an adjoint gradient-based data assimilation method that can efficiently handle high-dimensional parameters. We test this procedure on a synthetic data set and provide a clinical example with a dyssynchronous left ventricle with highly irregular motion. Our results show that the method efficiently handles a high-dimensional optimization parameter and produces an excellent agreement for personalized models to both synthetic and clinical data.

Improved identifiability of myocardial material parameters by an energy-based cost function

Myocardial stiffness is a valuable clinical biomarker for the monitoring and stratification of heart failure (HF). Cardiac finite element models provide a biomechanical framework for the assessment of stiffness through the determination of the myocardial constitutive model parameters. The reported parameter intercorrelations in popular constitutive relations, however, obstruct the unique estimation of material parameters and limit the reliable translation of this stiffness metric to clinical practice. Focusing on the role of the cost function (CF) in parameter identifiability, we investigate the performance of a set of geometric indices (based on displacements, strains, cavity volume, wall thickness and apicobasal dimension of the ventricle) and a novel CF derived from energy conservation. Our results, with a commonly used transversely isotropic material model (proposed by Guccione et al.), demonstrate that a single geometry-based CF is unable to uniquely constrain the parameter space. The energy-based CF, conversely, isolates one of the parameters and in conjunction with one of the geometric metrics provides a unique estimation of the parameter set. This gives rise to a new methodology for estimating myocardial material parameters based on the combination of deformation and energetics analysis. The accuracy of the pipeline is demonstrated in silico, and its robustness in vivo, in a total of 8 clinical data sets (7 HF and one control). The mean identified parameters of the Guccione material law were [Formula: see text] and [Formula: see text] ([Formula: see text], [Formula: see text], [Formula: see text]) for the HF cases and [Formula: see text] and [Formula: see text] ([Formula: see text], [Formula: see text], [Formula: see text]) for the healthy case.

Improved identifiability of myocardial material parameters by an energy-based cost function

Myocardial stiffness is a valuable clinical biomarker for the monitoring and stratification of heart failure (HF). Cardiac finite element models provide a biomechanical framework for the assessment of stiffness through the determination of the myocardial constitutive model parameters. The reported parameter intercorrelations in popular constitutive relations, however, obstruct the unique estimation of material parameters and limit the reliable translation of this stiffness metric to clinical practice. Focusing on the role of the cost function (CF) in parameter identifiability, we investigate the performance of a set of geometric indices (based on displacements, strains, cavity volume, wall thickness and apicobasal dimension of the ventricle) and a novel CF derived from energy conservation. Our results, with a commonly used transversely isotropic material model (proposed by Guccione et al.), demonstrate that a single geometry-based CF is unable to uniquely constrain the parameter space. The energy-based CF, conversely, isolates one of the parameters and in conjunction with one of the geometric metrics provides a unique estimation of the parameter set. This gives rise to a new methodology for estimating myocardial material parameters based on the combination of deformation and energetics analysis. The accuracy of the pipeline is demonstrated in silico, and its robustness in vivo, in a total of 8 clinical data sets (7 HF and one control). The mean identified parameters of the Guccione material law were [Formula: see text] and [Formula: see text] ([Formula: see text], [Formula: see text], [Formula: see text]) for the HF cases and [Formula: see text] and [Formula: see text] ([Formula: see text], [Formula: see text], [Formula: see text]) for the healthy case.

Estimation of passive and active properties in the human heart using 3D tagged MRI

Advances in medical imaging and image processing are paving the way for personalised cardiac biomechanical modelling. Models provide the capacity to relate kinematics to dynamics and-through patient-specific modelling-derived material parameters to underlying cardiac muscle pathologies. However, for clinical utility to be achieved, model-based analyses mandate robust model selection and parameterisation. In this paper, we introduce a patient-specific biomechanical model for the left ventricle aiming to balance model fidelity with parameter identifiability. Using non-invasive data and common clinical surrogates, we illustrate unique identifiability of passive and active parameters over the full cardiac cycle. Identifiability and accuracy of the estimates in the presence of controlled noise are verified with a number of in silico datasets. Unique parametrisation is then obtained for three datasets acquired in vivo. The model predictions show good agreement with the data extracted from the images providing a pipeline for personalised biomechanical analysis.

A self-taught artificial agent for multi-physics computational model personalization

Personalization is the process of fitting a model to patient data, a critical step towards application of multi-physics computational models in clinical practice. Designing robust personalization algorithms is often a tedious, time-consuming, model- and data-specific process. We propose to use artificial intelligence concepts to learn this task, inspired by how human experts manually perform it. The problem is reformulated in terms of reinforcement learning. In an off-line phase, Vito, our self-taught artificial agent, learns a representative decision process model through exploration of the computational model: it learns how the model behaves under change of parameters. The agent then automatically learns an optimal strategy for on-line personalization. The algorithm is model-independent; applying it to a new model requires only adjusting few hyper-parameters of the agent and defining the observations to match. The full knowledge of the model itself is not required. Vito was tested in a synthetic scenario, showing that it could learn how to optimize cost functions generically. Then Vito was applied to the inverse problem of cardiac electrophysiology and the personalization of a whole-body circulation model. The obtained results suggested that Vito could achieve equivalent, if not better goodness of fit than standard methods, while being more robust (up to 11% higher success rates) and with faster (up to seven times) convergence rate. Our artificial intelligence approach could thus make personalization algorithms generalizable and self-adaptable to any patient and any model.

Multiphysics and multiscale modelling, data-model fusion and integration of organ physiology in the clinic: ventricular cardiac mechanics

With heart and cardiovascular diseases continually challenging healthcare systems worldwide, translating basic research on cardiac (patho)physiology into clinical care is essential. Exacerbating this already extensive challenge is the complexity of the heart, relying on its hierarchical structure and function to maintain cardiovascular flow. Computational modelling has been proposed and actively pursued as a tool for accelerating research and translation. Allowing exploration of the relationships between physics, multiscale mechanisms and function, computational modelling provides a platform for improving our understanding of the heart. Further integration of experimental and clinical data through data assimilation and parameter estimation techniques is bringing computational models closer to use in routine clinical practice. This article reviews developments in computational cardiac modelling and how their integration with medical imaging data is providing new pathways for translational cardiac modelling.

Adjoint multi-start-based estimation of cardiac hyperelastic material parameters using shear data

Cardiac muscle tissue during relaxation is commonly modeled as a hyperelastic material with strongly nonlinear and anisotropic stress response. Adapting the behavior of such a model to experimental or patient data gives rise to a parameter estimation problem which involves a significant number of parameters. Gradient-based optimization algorithms provide a way to solve such nonlinear parameter estimation problems with relatively few iterations, but require the gradient of the objective functional with respect to the model parameters. This gradient has traditionally been obtained using finite differences, the calculation of which scales linearly with the number of model parameters, and introduces a differencing error. By using an automatically derived adjoint equation, we are able to calculate this gradient more efficiently, and with minimal implementation effort. We test this adjoint framework on a least squares fitting problem involving data from simple shear tests on cardiac tissue samples. A second challenge which arises in gradient-based optimization is the dependency of the algorithm on a suitable initial guess. We show how a multi-start procedure can alleviate this dependency. Finally, we provide estimates for the material parameters of the Holzapfel and Ogden strain energy law using finite element models together with experimental shear data.

Estimation of passive and active properties in the human heart using 3D tagged MRI

Advances in medical imaging and image processing are paving the way for personalised cardiac biomechanical modelling. Models provide the capacity to relate kinematics to dynamics and-through patient-specific modelling-derived material parameters to underlying cardiac muscle pathologies. However, for clinical utility to be achieved, model-based analyses mandate robust model selection and parameterisation. In this paper, we introduce a patient-specific biomechanical model for the left ventricle aiming to balance model fidelity with parameter identifiability. Using non-invasive data and common clinical surrogates, we illustrate unique identifiability of passive and active parameters over the full cardiac cycle. Identifiability and accuracy of the estimates in the presence of controlled noise are verified with a number of in silico datasets. Unique parametrisation is then obtained for three datasets acquired in vivo. The model predictions show good agreement with the data extracted from the images providing a pipeline for personalised biomechanical analysis.

Clinical Diagnostic Biomarkers from the Personalization of Computational Models of Cardiac Physiology

Computational modelling of the heart is rapidly advancing to the point of clinical utility. However, the difficulty of parameterizing and validating models from clinical data indicates that the routine application of truly predictive models remains a significant challenge. We argue there is significant value in an intermediate step towards prediction. This step is the use of biophysically based models to extract clinically useful information from existing patient data. Specifically in this paper we review methodologies for applying modelling frameworks for this goal in the areas of quantifying cardiac anatomy, estimating myocardial stiffness and optimizing measurements of coronary perfusion. Using these indicative examples of the general overarching approach, we finally discuss the value, ongoing challenges and future potential for applying biophysically based modelling in the clinical context.

Parameter estimation in a Holzapfel-Ogden law for healthy myocardium

A central problem in biomechanical studies of personalized human left ventricular (LV) modelling is to estimate material properties from in vivo clinical measurements. In this work we evaluate the passive myocardial mechanical properties inversely from the in vivo LV chamber pressure-volume and strain data. The LV myocardium is described using a structure-based orthotropic Holzapfel-Ogden constitutive law with eight parameters. In the first part of the paper we demonstrate how to use a multi-step non-linear least-squares optimization procedure to inversely estimate the parameters from the pressure-volume and strain data obtained from a synthetic LV model in diastole. In the second part, we show that to apply this procedure to clinical situations with limited in vivo data, additional constraints are required in the optimization procedure. Our study, based on three different healthy volunteers, demonstrates that the parameters of the Holzapfel-Ogden law could be extracted from pressure-volume and strain data with a suitable multi-step optimization procedure. Although the uniqueness of the solution cannot be addressed using our approaches, the material response is shown to be robustly determined.

Analysis of passive cardiac constitutive laws for parameter estimation using 3D tagged MRI

An unresolved issue in patient-specific models of cardiac mechanics is the choice of an appropriate constitutive law, able to accurately capture the passive behavior of the myocardium, while still having uniquely identifiable parameters tunable from available clinical data. In this paper, we aim to facilitate this choice by examining the practical identifiability and model fidelity of constitutive laws often used in cardiac mechanics. Our analysis focuses on the use of novel 3D tagged MRI, providing detailed displacement information in three dimensions. The practical identifiability of each law is examined by generating synthetic 3D tags from in silico simulations, allowing mapping of the objective function landscape over parameter space and comparison of minimizing parameter values with original ground truth values. Model fidelity was tested by comparing these laws with the more complex transversely isotropic Guccione law, by characterizing their passive end-diastolic pressure–volume relation behavior, as well as by considering the in vivo case of a healthy volunteer. These results show that a reduced form of the Holzapfel–Ogden law provides the best balance between identifiability and model fidelity across the tests considered.

Estimating passive mechanical properties in a myocardial infarction using MRI and finite element simulations

Myocardial infarction (MI) triggers a series of maladaptive events that lead to structural and functional changes in the left ventricle. It is crucial to better understand the progression of adverse remodeling, in order to develop effective treatment. In addition, being able to assess changes in vivo would be a powerful tool in the clinic. The goal of the current study is to quantify the in vivo material properties of infarcted and remote myocardium 1 week after MI, as well as the orientation of collagen fibers in the infarct. This will be accomplished by using a combination of magnetic resonance imaging (MRI), catheterization, finite element modeling, and numerical optimization to analyze a porcine model ([Formula: see text]) of posterolateral myocardial infarction. Specifically, properties will be determined by minimizing the difference between in vivo strains and volume calculated from MRI and finite element model predicted strains and volume. The results indicate that the infarct region is stiffer than the remote region and that the infarct collagen fibers become more circumferentially oriented 1 week post-MI. These findings are consistent with previous studies, which employed ex vivo techniques. The proposed methodology will ultimately provide a means of predicting remote and infarct mechanical properties in vivo at any time point post-MI.

Estimating passive mechanical properties in a myocardial infarction using MRI and finite element simulations

Myocardial infarction (MI) triggers a series of maladaptive events that lead to structural and functional changes in the left ventricle. It is crucial to better understand the progression of adverse remodeling, in order to develop effective treatment. In addition, being able to assess changes in vivo would be a powerful tool in the clinic. The goal of the current study is to quantify the in vivo material properties of infarcted and remote myocardium 1 week after MI, as well as the orientation of collagen fibers in the infarct. This will be accomplished by using a combination of magnetic resonance imaging (MRI), catheterization, finite element modeling, and numerical optimization to analyze a porcine model ([Formula: see text]) of posterolateral myocardial infarction. Specifically, properties will be determined by minimizing the difference between in vivo strains and volume calculated from MRI and finite element model predicted strains and volume. The results indicate that the infarct region is stiffer than the remote region and that the infarct collagen fibers become more circumferentially oriented 1 week post-MI. These findings are consistent with previous studies, which employed ex vivo techniques. The proposed methodology will ultimately provide a means of predicting remote and infarct mechanical properties in vivo at any time point post-MI.

Improved measurement of brain deformation during mild head acceleration using a novel tagged MRI sequence

In vivo measurements of human brain deformation during mild acceleration are needed to help validate computational models of traumatic brain injury and to understand the factors that govern the mechanical response of the brain. Tagged magnetic resonance imaging is a powerful, noninvasive technique to track tissue motion in vivo which has been used to quantify brain deformation in live human subjects. However, these prior studies required from 72 to 144 head rotations to generate deformation data for a single image slice, precluding its use to investigate the entire brain in a single subject. Here, a novel method is introduced that significantly reduces temporal variability in the acquisition and improves the accuracy of displacement estimates. Optimization of the acquisition parameters in a gelatin phantom and three human subjects leads to a reduction in the number of rotations from 72 to 144 to as few as 8 for a single image slice. The ability to estimate accurate, well-resolved, fields of displacement and strain in far fewer repetitions will enable comprehensive studies of acceleration-induced deformation throughout the human brain in vivo.

A computational pipeline for quantification of mouse myocardial stiffness parameters

The mouse is an important model for theoretical-experimental cardiac research, and biophysically based whole organ models of the mouse heart are now within reach. However, the passive material properties of mouse myocardium have not been much studied. We present an experimental setup and associated computational pipeline to quantify these stiffness properties. A mouse heart was excised and the left ventricle experimentally inflated from 0 to 1.44kPa in eleven steps, and the resulting deformation was estimated by echocardiography and speckle tracking. An in silico counterpart to this experiment was built using finite element methods and data on ventricular tissue microstructure from diffusion tensor MRI. This model assumed a hyperelastic, transversely isotropic material law to describe the force-deformation relationship, and was simulated for many parameter scenarios, covering the relevant range of parameter space. To identify well-fitting parameter scenarios, we compared experimental and simulated outcomes across the whole range of pressures, based partly on gross phenotypes (volume, elastic energy, and short- and long-axis diameter), and partly on node positions in the geometrical mesh. This identified a narrow region of experimentally compatible values of the material parameters. Estimation turned out to be more precise when based on changes in gross phenotypes, compared to the prevailing practice of using displacements of the material points. We conclude that the presented experimental setup and computational pipeline is a viable method that deserves wider application.

Cardiovascular magnetic resonance: deeper insights through bioengineering

Heart disease is the main cause of morbidity and mortality worldwide, with coronary artery disease, diabetes, and obesity being major contributing factors. Cardiovascular magnetic resonance (CMR) can provide a wealth of quantitative information on the performance of the heart, without risk to the patient. Quantitative analyses of these data can substantially augment the diagnostic quality of CMR examinations and can lead to more effective characterization of disease and quantification of treatment benefit. This review provides an overview of the current state of the art in CMR with particular regard to the quantification of motion, both microscopic and macroscopic, and the application of bioengineering analysis for the evaluation of cardiac mechanics. We discuss the current clinical practice and the likely advances in the next 5-10 years, as well as the ways in which clinical examinations can be augmented by bioengineering analysis of strain, compliance, and stress.

The estimation of patient-specific cardiac diastolic functions from clinical measurements

An unresolved issue in patients with diastolic dysfunction is that the estimation of myocardial stiffness cannot be decoupled from diastolic residual active tension (AT) because of the impaired ventricular relaxation during diastole. To address this problem, this paper presents a method for estimating diastolic mechanical parameters of the left ventricle (LV) from cine and tagged MRI measurements and LV cavity pressure recordings, separating the passive myocardial constitutive properties and diastolic residual AT. Dynamic C1-continuous meshes are automatically built from the anatomy and deformation captured from dynamic MRI sequences. Diastolic deformation is simulated using a mechanical model that combines passive and active material properties. The problem of non-uniqueness of constitutive parameter estimation using the well known Guccione law is characterized by reformulation of this law. Using this reformulated form, and by constraining the constitutive parameters to be constant across time points during diastole, we separate the effects of passive constitutive properties and the residual AT during diastolic relaxation. Finally, the method is applied to two clinical cases and one control, demonstrating that increased residual AT during diastole provides a potential novel index for delineating healthy and pathological cases.

Efficient computational methods for strongly coupled cardiac electromechanics

Strongly coupled cardiac electromechanical models can further our understanding of the relative importance of feedback mechanisms in the heart, but computational challenges currently remain a major obstacle, which limit their widespread use. To address this issue, we present a set of efficient computational methods including an efficient adaptive cell model integration scheme and a solution method for the monodomain equations that maintains high conduction velocity for time steps greater than 0.1 ms. We also present a novel method for increasing the efficiency of simulating electromechanical coupling, which shows a significant reduction in computational cost of the mechanical component on a personalized left ventricular geometry with an active contraction cell model reparametrized for human cells.

A computationally efficient formal optimization of regional myocardial contractility in a sheep with left ventricular aneurysm

A non-invasive method for estimating regional myocardial contractility in vivo would be of great value in the design and evaluation of new surgical and medical strategies to treat and/or prevent infarction-induced heart failure. As a first step towards developing such a method, an explicit finite element (FE) model-based formal optimization of regional myocardial contractility in a sheep with left ventricular (LV) aneurysm was performed using tagged magnetic resonance (MR) images and cardiac catheterization pressures. From the tagged MR images, 3-dimensional (3D) myocardial strains, LV volumes and geometry for the animal-specific 3D FE model of the LV were calculated, while the LV pressures provided physiological loading conditions. Active material parameters (T(max_B) and T(max_R)) in the non-infarcted myocardium adjacent to the aneurysm (borderzone) and in myocardium remote from the aneurysm were estimated by minimizing the errors between FE model-predicted and measured systolic strains and LV volumes using the successive response surface method for optimization. The significant depression in optimized T(max_B) relative to T(max_R) was confirmed by direct ex vivo force measurements from skinned fiber preparations. The optimized values of T(max_B) and T(max_R) were not overly sensitive to the passive material parameters specified. The computation time of less than 5 hours associated with our proposed method for estimating regional myocardial contractility in vivo makes it a potentially very useful clinical tool.

A finite element study of invariant-based orthotropic constitutive equations in the context of myocardial material parameter estimation

A previous study investigated a number of invariant-based orthotropic and transversely isotropic constitutive equations for their suitability to fit three-dimensional simple shear mechanics data of passive myocardial tissue. The study was based on the assumption of a homogeneous deformation. Here, we extend the previous study by performing an inverse finite element material parameter estimation. This ensures a more realistic deformation state and material parameter estimates. The constitutive relations were compared on the basis of (i) 'goodness of fit': how well they fit a set of six shear deformation tests and (ii) 'variability': how well determined the material parameters are over the range of experiments. These criteria were utilised to discuss the advantages and disadvantages of the constitutive relations. It was found that a specific form of the polyconvex type as well as the exponential Fung-type equations were most suitable for modelling the orthotropic behaviour of myocardium under simple shear.

Coupling multi-physics models to cardiac mechanics

We outline and review the mathematical framework for representing mechanical deformation and contraction of the cardiac ventricles, and how this behaviour integrates with other processes crucial for understanding and modelling heart function. Building on general conservation principles of space, mass and momentum, we introduce an arbitrary Eulerian-Lagrangian framework governing the behaviour of both fluid and solid components. Exploiting the natural alignment of cardiac mechanical properties with the tissue microstructure, finite deformation measures and myocardial constitutive relations are referred to embedded structural axes. Coupling approaches for solving this large deformation mechanics framework with three dimensional fluid flow, coronary hemodynamics and electrical activation are described. We also discuss the potential of cardiac mechanics modelling for clinical applications.

Modelling passive diastolic mechanics with quantitative MRI of cardiac structure and function

The majority of patients with clinically diagnosed heart failure have normal systolic pump function and are commonly categorized as suffering from diastolic heart failure. The left ventricle (LV) remodels its structure and function to adapt to pathophysiological changes in geometry and loading conditions, which in turn can alter the passive ventricular mechanics. In order to better understand passive ventricular mechanics, a LV finite element (FE) model was customized to geometric data segmented from in vivo tagged magnetic resonance images (MRI) data and myofibre orientation derived from ex vivo diffusion tensor MRI (DTMRI) of a canine heart using nonlinear finite element fitting techniques. MRI tissue tagging enables quantitative evaluation of cardiac mechanical function with high spatial and temporal resolution, whilst the direction of maximum water diffusion in each voxel of a DTMRI directly corresponds to the local myocardial fibre orientation. Due to differences in myocardial geometry between in vivo and ex vivo imaging, myofibre orientations were mapped into the geometric FE model using host mesh fitting (a free form deformation technique). Pressure recordings, temporally synchronized to the tagging data, were used as the loading constraints to simulate the LV deformation during diastole. Simulation of diastolic LV mechanics allowed us to estimate the stiffness of the passive LV myocardium based on kinematic data obtained from tagged MRI. Integrated physiological modelling of this kind will allow more insight into mechanics of the LV on an individualized basis, thereby improving our understanding of the underlying structural basis of mechanical dysfunction under pathological conditions.

Myocardial material parameter estimation: a comparison of invariant based orthotropic constitutive equations

This study investigated a number of invariant based orthotropic and transversely isotropic constitutive equations for their suitability to fit three-dimensional simple shear mechanics data of passive myocardial tissue. A number of orthotropic laws based on Green strain components and one microstructurally based law have previously been investigated to fit experimental measurements of stress-strain behaviour. Here we extend this investigation to include several recently proposed functional forms, i.e. invariant based orthotropic and transversely isotropic constitutive relations. These laws were compared on the basis of (i) 'goodness of fit': how well they fit a set of six shear deformation tests, (ii) 'variability': how well determined the material parameters are over the range of experiments. These criteria were utilised to discuss the advantages and disadvantages of the constitutive laws. It was found that a specific form of the polyconvex type as well as the exponential Fung-type law from the previous study were most suitable for modelling the orthotropic behaviour of myocardium under simple shear.

Passive ventricular mechanics modelling using MRI of structure and function

Patients suffering from dilated cardiomyopathy or myocardial infarction can develop left ventricular (LV) diastolic impairment. The LV remodels its structure and function to adapt to pathophysiological changes in geometry and loading conditions and this remodeling process can alter the passive ventricular mechanics. In order to better understand passive ventricular mechanics, a LV finite element model was developed to incorporate physiological and mechanical information derived from in vivo magnetic resonance imaging (MRI) tissue tagging, in vivo LV cavity pressure recording and ex vivo diffusion tensor MRI (DTMRI) of a canine heart. MRI tissue tagging enables quantitative evaluation of cardiac mechanical function with high spatial and temporal resolution, whilst the direction of maximum water diffusion (the primary eigenvector) in each voxel of a DTMRI directly correlates with the myocardial fibre orientation. This model was customized to the geometry of the canine LV during diastasis by fitting the segmented epicardial and endocardial surface data from tagged MRI using nonlinear finite element fitting techniques. Myofibre orientations, extracted from DTMRI of the same heart, were incorporated into this geometric model using a free form deformation methodology. Pressure recordings, temporally synchronized to the tissue tagging MRI data, were used to simulate the LV deformation during diastole. Simulation of the diastolic LV mechanics allowed us to estimate the stiffness of the passive LV myocardium based on kinematic data obtained from tagged MRI. This integrated physiological model will allow more insight into the regional passive diastolic mechanics of the LV on an individualized basis, thereby improving our understanding of the underlying structural basis of mechanical dysfunction in pathological conditions.

Estimating material parameters of human skin in vivo

An accurate mathematical representation of the mechanical behaviour of human skin is essential when simulating deformations occurring in the skin during body movements or clinical procedures. In this study constitutive stress-strain relationships based on experimental data from human skin in vivo were obtained. A series of multiaxial loading experiments were performed on the forearms of four age- and gender matched subjects. The tissue geometry, together with recorded displacements and boundary forces, were combined in an analysis using finite element modelling. A non-linear optimization technique was developed to estimate values for the material parameters of a previously published constitutive law, describing the stress-strain relationship as a non-linear anisotropic membrane. Ten sets of material parameters where estimated from the experiments, showing considerable differences in mechanical behaviour both between individual subjects as well as mirrored body locations on a single subject. The accuracy of applications that simulate large deformations of human skin could be improved by using the parameters found from the in vivo experiments as described in this study.

Estimation of cardiac hyperelastic material properties from MRI tissue tagging and diffusion tensor imaging

The passive material properties of myocardium are important in the understanding of diastolic cardiac dysfunction. We determined hyperelastic myocardial material parameters in four isolated arrested pig hearts undergoing passive inflation of the left ventricle. Using geometry from MRI, recorded boundary conditions, muscle fiber architecture from diffusion tensor imaging, and deformation from tissue tagging, finite element models were constructed to solve the finite elasticity stress estimation problem. The constitutive parameters of a hyperelastic transversely isotropic material law were determined by minimizing the difference between the predicted and imaged deformation field. The optimized parameters were in a similar range as those reported by previous studies, showing increased passive stiffness in the muscle fiber direction. The average RMS error was 0.92 mm, similar to the image resolution of 0.80 mm. Optimization of hyperelastic models of myocardial mechanics can thus be performed to extract meaningful biophysical parameters from MRI data.

A biomechanical model of mammographic compressions

A number of biomechanical models have been proposed to improve nonrigid registration techniques for multimodal breast image alignment. A deformable breast model may also be useful for overcoming difficulties in interpreting 2D X-ray projections (mammograms) of 3D volumes (breast tissues). If a deformable model could accurately predict the shape changes that breasts undergo during mammography, then the model could serve to localize suspicious masses (visible in mammograms) in the unloaded state, or in any other deformed state required for further investigations (such as biopsy or other medical imaging modalities). In this paper, we present a validation study that was conducted in order to develop a biomechanical model based on the well-established theory of continuum mechanics (finite elasticity theory with contact mechanics) and demonstrate its use for this application. Experimental studies using gel phantoms were conducted to test the accuracy in predicting mammographic-like deformations. The material properties of the gel phantom were estimated using a nonlinear optimization process, which minimized the errors between the experimental and the model-predicted surface data by adjusting the parameter associated with the neo-Hookean constitutive relation. Two compressions (the equivalent of cranio-caudal and medio-lateral mammograms) were performed on the phantom, and the corresponding deformations were recorded using a MRI scanner. Finite element simulations were performed to mimic the experiments using the estimated material properties with appropriate boundary conditions. The simulation results matched the experimental recordings of the deformed phantom, with a sub-millimeter root-mean-square error for each compression state. Having now validated our finite element model of breast compression, the next stage is to apply the model to clinical images.

A Computationally Efficient Optimization Kernel for Material Parameter Estimation Procedures

Estimating material parameters is an important part in the study of soft tissue mechanics. Computational time can easily run to days, especially when all available experimental data are taken into account. The material parameter estimation procedure is examplified on a set of homogeneous simple shear experiments to estimate the orthotropic constitutive parameters of myocardium. The modification consists of changing the traditional least-squares approach to a weighted least-squares. This objective function resembles a L2-norm type integral which is approximated using Gaussian quadrature. This reduces the computational time of the material parameter estimation by two orders of magnitude.

Dobutamine cardiovascular magnetic resonance: A review

Dobutamine cardiovascular magnetic resonance (DCMR) is useful for identifying myocardial ischemia and viability in patients with known or suspected coronary artery disease (CAD). This article reviews the performance and utility of DCMR, its association with dobutamine stress echocardiography (DSE), and areas of active investigative research.