Initial stress in biomechanical models of atherosclerotic plaques

Rupture of atherosclerotic plaques is the underlying cause for the majority of acute strokes and myocardial infarctions. Rupture of the plaque occurs when the stress in the plaque exceeds the strength of the material locally. Biomechanical stress analyses are commonly based on pressurized geometries, in most cases measured by in-vivo MRI. The geometry is therefore not stress-free. The aim of this study is to identify the effect of neglecting the initial stress state on the plaque stress distribution. Fifty 2D histological sections (7 patients, 9 diseased coronary artery segments), perfusion fixed at 100 mmHg, were segmented and finite element models were created. The Backward Incremental method was applied to determine the initial stress state and the zero-pressure state. Peak plaque and cap stresses were compared with and without initial stress. The effect of initial stress on the peak stress was related to the minimum cap thickness, maximum necrotic core thickness, and necrotic core angle. When accounting for initial stress, the general relations between geometrical features and peak cap stress remain intact. However, on a patient-specific basis, accounting for initial stress has a different effect on the absolute cap stress for each plaque. Incorporating initial stress may therefore improve the accuracy of future stress based rupture risk analyses for atherosclerotic plaques.     

The effects of plaque morphology and material properties on peak cap stress in human coronary arteries

Heart attacks are often caused by rupture of caps of atherosclerotic plaques in coronary arteries. Cap rupture occurs when cap stress exceeds cap strength. We investigated the effects of plaque morphology and material properties on cap stress. Histological data from 77 coronary lesions were obtained and segmented. In these patient-specific cross sections, peak cap stresses were computed by using finite element analyses. The finite element analyses were 2D, assumed isotropic material behavior, and ignored residual stresses. To represent the wide spread in material properties, we applied soft and stiff material models for the intima. Measures of geometric plaque features for all lesions were determined and their relations to peak cap stress were examined using regression analyses. Patient-specific geometrical plaque features greatly influence peak cap stresses. Especially, local irregularities in lumen and necrotic core shape as well as a thin intima layer near the shoulder of the plaque induce local stress maxima. For stiff models, cap stress increased with decreasing cap thickness and increasing lumen radius (R = 0.79). For soft models, this relationship changed: increasing lumen radius and increasing lumen curvature were associated with increased cap stress (R = 0.66). The results of this study imply that not only accurate assessment of plaque geometry, but also of intima properties is essential for cap stress analyses in atherosclerotic plaques in human coronary arteries.     

Advanced human carotid plaque progression correlates positively with flow shear stress using follow-up scan data: An in vivo MRI multi-patient 3D FSI study

Although it has been well-accepted that atherosclerosis initiation and early progression correlate negatively with flow wall shear stresses (FSS), increasing evidence suggests mechanisms governing advanced plaque progression are not well understood. Fourteen patients were scanned 2–4 times at 18 month intervals using a histologically validated multi-contrast magnetic resonance imaging (MRI) protocol to acquire carotid plaque progression data. Thirty-two scan pairs (baseline and follow-up scans) were formed with slices matched for model construction and analysis. 3D fluid–structure interaction (FSI) models were constructed and plaque wall stress (PWS) and flow shear stress (FSS) were obtained from all matching lumen data points (400–1000 per plaque; 100 points per matched slice) to quantify correlations with plaque progression measured by vessel wall thickness increase (WTI). Using FSS and PWS data from follow-up scan, 21 out of 32 scan pairs showed a significant positive correlation between WTI and FSS (positive/negative/no significance ratio=21/8/3), and 26 out of 32 scan pairs showed a significant negative correlation between WTI and PWS (positive/negative/no significance ratio=2/26/4). The mean FSS value of lipid core nodes (n=5294) from all 47 plaque models was 63.5 dyn/cm², which was 45% higher than that from all normal vessel nodes (n=27553, p<0.00001). The results from this intensive FSI study indicate that flow shear stress from follow-up scan correlates positively with advanced plaque progression which is different from what has been observed in plaque initiation and early-stage progression. It should be noted that the correlation results do not automatically lead to any causality conclusions.     

Anisotropic mechanical properties of tissue components in human atherosclerotic plaques

Knowledge of the biomechanical properties of human atherosclerotic plaques is of essential importance for developing more insights in the pathophysiology of the cardiovascular system and for better predicting the outcome of interventional treatments such as balloon angioplasty. Available data are mainly based on uniaxial tests, and most of the studies investigate the mechanical response of fibrous plaque caps only. However, stress distributions during, for example, balloon angioplasty are strongly influenced by all components of atherosclerotic lesions. A total number of 107 samples from nine human high-grade stenotic iliac arteries were tested; associated anamnesis of donors reported. Magnetic resonance imaging was employed to test the usability of the harvested arteries. Histological analyses has served to characterize the different tissue types. Prepared strips of 7 different tissue types underwent cyclic quasistatic uniaxial tension tests in axial and circumferential directions; ultimate tensile stresses and stretches were documented. Experimental data of individual samples indicated anisotropic and highly nonlinear tissue properties as well as considerable interspecimen differences. The calcification showed, however a linear property, with about the same stiffness as observed for the adventitia in high stress regions. The stress and stretch values at calcification fracture are smaller (179 +/- 56 kPa and 1.02 +/- 0.005) than for each of the other tissue components. Of all intimal tissues investigated, the lowest fracture stress occurred in the circumferential direction of the fibrous cap (254.8 +/- 79.8 kPa at stretch 1.182 +/- 0.1). The adventitia demonstrated the highest and the nondiseased media the lowest mechanical strength on average.     

Tensile and compressive properties of fresh human carotid atherosclerotic plaques

Accurate characterisation of the mechanical properties of human atherosclerotic plaque is important for our understanding of the role of vascular mechanics in the development and treatment of atherosclerosis. The majority of previous studies investigating the mechanical properties of human plaque are based on tests of plaque tissue removed following autopsy. This study aims to characterise the mechanical behaviour of fresh human carotid plaques removed during endarterectomy and tested within 2 h. A total of 50 radial compressive and 17 circumferential tensile uniaxial tests were performed on samples taken from 14 carotid plaques. The clinical classification of each plaque, as determined by duplex ultrasound is also reported. Plaques were classified as calcified, mixed or echolucent. Experimental data indicated that plaques were highly inhomogeneous; with variations seen in the mechanical properties of plaque obtained from individual donors and between donors. The mean behaviour of samples for each classification indicated that calcified plaques had the stiffest response, while echolucent plaques were the least stiff. Results also indicated that there may be a difference in behaviour of samples taken from different anatomical locations (common, internal and external carotid), however the large variability indicates that more testing is needed to reach significant conclusions. This work represents a step towards a better understanding of the in vivo mechanical behaviour of human atherosclerotic plaque.     

Effect of a lipid pool on stress/strain distributions in stenotic arteries: 3-D fluid-structure interactions (FSI) models

Nonlinear 3-D models with fluid-structure interactions (FSI) based on in vitro experiments are introduced and solved by ADINA to perform flow and stress/strain analysis for stenotic arteries with lipid cores. Navier-Stokes equations are used as the governing equations for the fluid. Hyperelastic Mooney-Rivlin models are used for both the arteries and lipid cores. Our results indicate that critical plaque stress/strain conditions are affected considerably by stenosis severity, eccentricity, lipid pool size, shape and position, plaque cap thickness, axial stretch, pressure, and fluid-structure interactions, and may be used for possible plaque rupture predictions.     

Simulation of human atherosclerotic femoral plaque tissue: the influence of plaque material model on numerical results

Background

Due to the limited number of experimental studies that mechanically characterise human atherosclerotic plaque tissue from the femoral arteries, a recent trend has emerged in current literature whereby one set of material data based on aortic plaque tissue is employed to numerically represent diseased femoral artery tissue. This study aims to generate novel vessel-appropriate material models for femoral plaque tissue and assess the influence of using material models based on experimental data generated from aortic plaque testing to represent diseased femoral arterial tissue.

Methods

Novel material models based on experimental data generated from testing of atherosclerotic femoral artery tissue are developed and a computational analysis of the revascularisation of a quarter model idealised diseased femoral artery from a 90% diameter stenosis to a 10% diameter stenosis is performed using these novel material models. The simulation is also performed using material models based on experimental data obtained from aortic plaque testing in order to examine the effect of employing vessel appropriate material models versus those currently employed in literature to represent femoral plaque tissue.

Results

Simulations that employ material models based on atherosclerotic aortic tissue exhibit much higher maximum principal stresses within the plaque than simulations that employ material models based on atherosclerotic femoral tissue. Specifically, employing a material model based on calcified aortic tissue, instead of one based on heavily calcified femoral tissue, to represent diseased femoral arterial vessels results in a 487 fold increase in maximum principal stress within the plaque at a depth of 0.8 mm from the lumen.

Conclusions

Large differences are induced on numerical results as a consequence of employing material models based on aortic plaque, in place of material models based on femoral plaque, to represent a diseased femoral vessel. Due to these large discrepancies, future studies should seek to employ vessel-appropriate material models to simulate the response of diseased femoral tissue in order to obtain the most accurate numerical results.

     

Using in vivo Cine and 3D multi-contrast MRI to determine human atherosclerotic carotid artery material properties and circumferential shrinkage rate and their impact on stress/strain predictions

In vivo magnetic resonance image (MRI)-based computational models have been introduced to calculate atherosclerotic plaque stress and strain conditions for possible rupture predictions. However, patient-specific vessel material properties are lacking in those models, which affects the accuracy of their stress/strain predictions. A noninvasive approach of combining in vivo Cine MRI, multicontrast 3D MRI, and computational modeling was introduced to quantify patient-specific carotid artery material properties and the circumferential shrinkage rate between vessel in vivo and zero-pressure geometries. In vivo Cine and 3D multicontrast MRI carotid plaque data were acquired from 12 patients after informed consent. For each patient, one nearly-circular slice and an iterative procedure were used to quantify parameter values in the modified Mooney-Rivlin model for the vessel and the vessel circumferential shrinkage rate. A sample artery slice with and without a lipid core and three material parameter sets representing stiff, median, and soft materials from our patient data were used to demonstrate the effect of material stiffness and circumferential shrinkage process on stress/strain predictions. Parameter values of the Mooney-Rivlin models for the 12 patients were quantified. The effective Young's modulus (YM, unit: kPa) values varied from 137 (soft), 431 (median), to 1435 (stiff), and corresponding circumferential shrinkages were 32%, 12.6%, and 6%, respectively. Using the sample slice without the lipid core, the maximum plaque stress values (unit: kPa) from the soft and median materials were 153.3 and 96.2, which are 67.7% and 5% higher than that (91.4) from the stiff material, while the maximum plaque strain values from the soft and median materials were 0.71 and 0.293, which are about 700% and 230% higher than that (0.089) from the stiff material, respectively. Without circumferential shrinkages, the maximum plaque stress values (unit: kPa) from the soft, median, and stiff models were inflated to 330.7, 159.2, and 103.6, which were 116%, 65%, and 13% higher than those from models with proper shrinkage. The effective Young's modulus from the 12 human carotid arteries studied varied from 137 kPa to 1435 kPa. The vessel circumferential shrinkage to the zero-pressure condition varied from 6% to 32%. The inclusion of proper shrinkage in models based on in vivo geometry is necessary to avoid over-estimating the stresses and strains by up 100%. Material stiffness had a greater impact on strain (up to 700%) than on stress (up to 70%) predictions. Accurate patient-specific material properties and circumferential shrinkage could considerably improve the accuracy of in vivo MRI-based computational stress/strain predictions.     

Quantifying dynamic mechanical properties of human placenta tissue using optimization techniques with specimen-specific finite-element models

Motor-vehicle crashes are the leading cause of fetal deaths resulting from maternal trauma in the United States, and placental abruption is the most common cause of these deaths. To minimize this injury, new assessment tools, such as crash-test dummies and computational models of pregnant women, are needed to evaluate vehicle restraint systems with respect to reducing the risk of placental abruption. Developing these models requires accurate material properties for tissues in the pregnant abdomen under dynamic loading conditions that can occur in crashes. A method has been developed for determining dynamic material properties of human soft tissues that combines results from uniaxial tensile tests, specimen-specific finite-element models based on laser scans that accurately capture non-uniform tissue-specimen geometry, and optimization techniques. The current study applies this method to characterizing material properties of placental tissue. For 21 placenta specimens tested at a strain rate of 12/s, the mean failure strain is 0.472±0.097 and the mean failure stress is 34.80±12.62 kPa. A first-order Ogden material model with ground-state shear modulus (μ) of 23.97±5.52 kPa and exponent (α₁) of 3.66±1.90 best fits the test results. The new method provides a nearly 40% error reduction (p<0.001) compared to traditional curve-fitting methods by considering detailed specimen geometry, loading conditions, and dynamic effects from high-speed loading. The proposed method can be applied to determine mechanical properties of other soft biological tissues.     

Higher critical plaque wall stress in patients who died of coronary artery disease compared with those who died of other causes: A 3D FSI study based on ex vivo MRI of coronary plaques

Mechanical forces play an important role in the rupture of vulnerable plaques. This process is often associated with cardiovascular syndromes, such as heart attack and stroke. In this study, magnetic resonance imaging (MRI)-based models were used to investigate the association between plaque wall stress (PWS) and coronary artery disease (CAD).     

Profilin-1 is expressed in human atherosclerotic plaques and induces atherogenic effects on vascular smooth muscle cells

Background

Profilin-1 is an ubiquitous actin binding protein. Under pathological conditions such as diabetes, profilin-1 levels are increased in the vascular endothelium. We recently demonstrated that profilin-1 overexpression triggers indicators of endothelial dysfunction downstream of LDL signaling, and that attenuated expression of profilin-1 confers protection from atherosclerosis in vivo.

Methodology

Here we monitored profilin-1 expression in human atherosclerotic plaques by immunofluorescent staining. The effects of recombinant profilin-1 on atherogenic signaling pathways and cellular responses such as DNA synthesis (BrdU-incorporation) and chemotaxis (modified Boyden-chamber) were evaluated in cultured rat aortic and human coronary vascular smooth muscle cells (VSMCs). Furthermore, the correlation between profilin-1 serum levels and the degree of atherosclerosis was assessed in humans.

Principal Findings

In coronary arteries from patients with coronary heart disease, we found markedly enhanced profilin expression in atherosclerotic plaques compared to the normal vessel wall. Stimulation of rat aortic and human coronary VSMCs with recombinant profilin-1 (10⁻⁶ M) in vitro led to activation of intracellular signaling cascades such as phosphorylation of Erk1/2, p70^S6 kinase and PI3K/Akt within 10 minutes. Furthermore, profilin-1 concentration-dependently induced DNA-synthesis and migration of both rat and human VSMCs, respectively. Inhibition of PI3K (Wortmannin, {type:entrez-nucleotide,attrs:{text:LY294002,term_id:1257998346,term_text:LY294002}}LY294002) or Src-family kinases (SU6656, PP2), but not PLCγ ({type:entrez-nucleotide,attrs:{text:U73122,term_id:4098075,term_text:U73122}}U73122), completely abolished profilin-induced cell cycle progression, whereas PI3K inhibition partially reduced the chemotactic response. Finally, we found that profilin-1 serum levels were significantly elevated in patients with severe atherosclerosis in humans (p<0.001 vs. no atherosclerosis or control group).

Conclusions

Profilin-1 expression is significantly enhanced in human atherosclerotic plaques compared to the normal vessel wall, and the serum levels of profilin-1 correlate with the degree of atherosclerosis in humans. The atherogenic effects exerted by profilin-1 on VSMCs suggest an auto-/paracrine role within the plaque. These data indicate that profilin-1 might critically contribute to atherogenesis and may represent a novel therapeutic target.     

Mechanical characterization of 3D printed mimic of human artery affected by atherosclerotic plaque through numerical and experimental methods

Biomechanics and Modeling in Mechanobiology - This paper proposes a novel experimental investigation based on 3D printing to validate numerical models for biomechanics simulations. Soft elastomeric...     

Computational evaluation of oxygen and shear stress distributions in 3D perfusion culture systems: macro-scale and micro-structured models

We present a combined macro-scale/micro-scale computational approach to quantify oxygen transport and flow-mediated shear stress to human chondrocytes cultured in three-dimensional scaffolds in a perfusion bioreactor system. A macro-scale model was developed to assess the influence of the bioreactor design and to identify the proper boundary conditions for the micro-scale model. The micro-scale model based on a micro-computed tomography (microCT) reconstruction of a poly(ethylene glycol terephthalate)/poly(butylene terephthalate) (PEGT/PBT) foam scaffold, was developed to assess the influence of the scaffold micro-architecture on local shear stress and oxygen levels within the scaffold pores. Experiments were performed to derive specific oxygen consumption rates for constructs perfused under flow rates of 0.3 and 0.03 ml min(-1). While macro-scale and micro-scale models predicted similar average oxygen levels at different depths within the scaffold, microCT models revealed small local oxygen variations within the scaffold micro-architecture. The combined macro-scale/micro-scale approach indicated that 0.3 ml min(-1), which subjected 95% of the cells to less than 6.3 mPa shear, would maintain the oxygen supply throughout the scaffold above anoxic levels (>1%), with 99.5% of the scaffold supplied with 8-2% O(2). Alternatively, at 0.03 ml min(-1), the macro-scale model predicted 6% of the cells would be supplied with 0.5-1% O(2), although this region of cells was confined to the periphery of the scaffold. Together with local variations predicted by the micro-scale model, the simulations underline that in the current model system, reducing the flow below 0.03 ml min(-1) would likely have dire consequences on cell viability to pronounced regions within the engineered construct. The presented approach provides a sensitive tool to aid efficient bioreactor optimization and scaffold design.     

Construction of 3D human distal femoral surface models using a 3D statistical deformable model

Construction of 3D geometric surface models of human knee joint is always a challenge in biomedical engineering. This study introduced an improved statistical shape model (SSM) method that only uses 2D images of a joint to predict the 3D joint surface model. The SSM was constructed using 40 distal femur models of human knees. In this paper, a series validation and parametric analysis suggested that more than 25 distal femur models are needed to construct the SSM; each distal femur should be described using at least 3000 nodes in space; and two 2D fluoroscopic images taken in 45° directions should be used for the 3D surface shape prediction. Using this SSM method, ten independent distal femurs from 10 independent living subjects were predicted using their 2D plane fluoroscopic images. The predicted models were compared to their native 3D distal femur models constructed using their 3D MR images. The results demonstrated that using two fluoroscopic images of the knee, the overall difference between the predicted distal femur surface and the MR image-based surface was 0.16±1.16 mm. These data indicated that the SSM method could be a powerful method for construction of 3D surface geometries of the distal femur.     

Accumulation of glycosphingolipids in human atherosclerotic plaque and unaffected aorta tissues

We have measured the levels of glycosphingolipids and the activityof glycosphingolipid glycosyltransferases in human aortic intimaand media from patients who died of atherosclerosis. The effectsof lactosylceramide (LacCer) and glucosylceramide (GlcCer) fromplaque intima on smooth muscle cell proliferation were assessed.When the GIcCer data was expressed as (pg GlcCer/mg cholesteroland/mg total phospholipid, a 28-fold and 7-told increase inplaque intima compared to normal intima was observed. Similarly,the level of LacCer was elevated 5-fold and 4-fold, respectively,compared to unaffected intima. The activity of UDPGicCer: ceramideß1     

The influence of computational strategy on prediction of mechanical stress in carotid atherosclerotic plaques: Comparison of 2D structure-only, 3D structure-only,one-way and fully coupled fluid-structure interaction analyses

Background

Compositional and morphological features of carotid atherosclerotic plaques provide complementary information to luminal stenosis in predicting clinical presentations. However, they alone cannot predict cerebrovascular risk. Mechanical stress within the plaque induced by cyclical changes in blood pressure has potential to assess plaque vulnerability. Various modeling strategies have been employed to predict stress, including 2D and 3D structure-only, 3D one-way and fully coupled fluid-structure interaction (FSI) simulations. However, differences in stress predictions using different strategies have not been assessed.

Methods

Maximum principal stress (Stress-P1) within 8 human carotid atherosclerotic plaques was calculated based on geometry reconstructed from in vivo computerized tomography and high resolution, multi-sequence magnetic resonance images. Stress-P1 within the diseased region predicted by 2D and 3D structure-only, and 3D one-way FSI simulations were compared to 3D fully coupled FSI analysis.

Results

Compared to 3D fully coupled FSI, 2D structure-only simulation significantly overestimated stress level (94.1 kPa [65.2, 117.3] vs. 85.5 kPa [64.4, 113.6]; median [inter-quartile range], p=0.0004). However, when slices around the bifurcation region were excluded, stresses predicted by 2D structure-only simulations showed a good correlation (R²=0.69) with values obtained from 3D fully coupled FSI analysis. 3D structure-only model produced a small yet statistically significant stress overestimation compared to 3D fully coupled FSI (86.8 kPa [66.3, 115.8] vs. 85.5 kPa [64.4, 113.6]; p<0.0001). In contrast, one-way FSI underestimated stress compared to 3D fully coupled FSI (78.8 kPa [61.1, 100.4] vs. 85.5 kPa [64.4, 113.7]; p<0.0001).

Conclusions

A 3D structure-only model seems to be a computationally inexpensive yet reasonably accurate approximation for stress within carotid atherosclerotic plaques with mild to moderate luminal stenosis as compared to fully coupled FSI analysis.     

Identification and quantitation of unique fatty acid oxidation products in human atherosclerotic plaque using high-performance liquid chromatography

Oxidation of lipoproteins, particularly low-density lipoprotein, is thought to play a major role in the development of atherosclerosis. We set out to identify and quantitate the major fatty acid oxidation products in human atherosclerotic plaque obtained from individuals undergoing carotid endarterectomy. Oxidized lipids were extracted from plaque homogenate under conditions to prevent artifactual oxidation. Identification and quantitation was performed using HPLC and GC-MS. High levels of hydroxyoctadecanoic acids (0.51 +/- 0.17 ng/microg of linoleic acid), 15-hydroxyeicosatetranoic acid (HETE) (0.66 +/- 0.24 ng/microg of arachidonic acid), and 11-HETE (0.84 +/- 0.24 ng/microg of arachidonic acid) were detected in all atherosclerotic plaques (n = 10). Low levels of 9-oxo-octadecanoic acid (oxoODE) (0.04 +/- 0.01 ng/microg of linoleic acid), were present in all samples, while 13-oxoODE (0.01 +/- 0.008 ng/microg of linoleic acid) was present in only 4 of the 10 plaque samples. Of interest was the identification of two previously unidentified compounds in atherosclerotic plaque, 11-oxo-eicosatetranoic acid in 9 of the 10 samples and 5,6-dihydroxyeicosatetranoic acid in 3 samples. Chiral analysis revealed that all the major compounds identified in this study are of a nonenzymatic origin. This study is the first to provide a convenient HPLC method to quantify all the products of both linoleic acid and arachidonic acid oxidation in human atherosclerotic plaque. The quantitation of lipid peroxidation products in plaque may be important given the potential biological activity of these compounds and their possible relationship to plaque pathogenesis and instability.     

Altered dermatan sulfate structure and reduced heparin cofactor II-stimulating activity of biglycan and decorin from human atherosclerotic plaque

Biglycan and decorin are small dermatan sulfate-containing proteoglycans in the extracellular matrix of the artery wall. The dermatan sulfate chains are known to stimulate thrombin inhibition by heparin cofactor II (HCII), a plasma proteinase inhibitor that has been detected within the artery wall. The purpose of this study was to analyze the HCII-stimulatory activity of biglycan and decorin isolated from normal human aorta and atherosclerotic lesions type II through VI and to correlate activity with dermatan sulfate chain composition and structure. Biglycan and decorin from plaque exhibited a 24-75% and 38-79% loss of activity, respectively, in thrombin-HCII inhibition assays relative to proteoglycan from normal aorta. A significant negative linear relationship was observed between lesion severity and HCII stimulatory activity (r = 0.79, biglycan; r = 0.63, decorin; p < 0.05). Biglycan, but not decorin, from atherosclerotic plaque contained significantly reduced amounts of iduronic acid and disulfated disaccharides DeltaDi-2,4S and DeltaDi-4,6S relative to proteoglycan from normal artery. Affinity coelectrophoresis analysis of a subset of samples demonstrated that increased interaction of proteoglycan with HCII in agarose gels paralleled increased activity in thrombin-HCII inhibition assays. In conclusion, both biglycan and decorin from atherosclerotic plaque possessed reduced activity with HCII, but only biglycan demonstrated a correlation between activity and specific glycosaminoglycan structural features. Loss of the ability of biglycan and decorin in atherosclerotic lesions to regulate thrombin activity through HCII may be critical in the progression of the disease.     

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.     

Influence of non-Newtonian properties of blood on the wall shear stress in human atherosclerotic right coronary arteries

The objective of this work is to investigate the effect of non-Newtonian properties of blood on the wall shear stress (WSS) in atherosclerotic coronary arteries using both Newtonian and non-Newtonian models. Numerical simulations were performed to examine how the spatial and temporal WSS distributions are influenced by the stenosis size, blood viscosity, and flow rate. The computational results demonstrated that blood viscosity properties had considerable effect on the magnitude of the WSS, especially where disturbed flow was observed. The WSS distribution is highly non-uniform both temporally and spatially, especially in the stenotic region. The maximum WSS occurred at the proximal side of the stenosis, near the outer wall in the curved artery with no stenosis. The lumen area near the inner wall distal to the stenosis region experienced a lower WSS during the entire cardiac cycle. Among the factors of stenosis size, blood viscosity, and flow rate, the size of the stenosis has the most significant effect on the spatial and temporal WSS distributions qualitatively and quantitatively.