Multi-scale interaction of 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 their 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 flow, and proceed to propose, then validate (in practical, qualitative terms, at least) a multiply coupled, multi-scale modeling strategy, which, eventually, aims at a quantitative, patient-specific understanding of the coupling between the flow and the endothelial state

Linking quantitative radiology to molecular mechanism for improved vascular disease therapy selection and follow-up

Objective: Therapeutic advancements in atherosclerotic cardiovascular disease have improved the prevention of ischemic stroke and myocardial infarction. However, diagnostic methods for atherosclerotic plaque phenotyping to aid individualized therapy are lacking. In this thesis, we aimed to elucidate plaque biology through the analysis of computed-tomography angiography (CTA) with sufficient sensitivity and specificity to capture the differentiated drivers of the disease. We then aimed to use such data to calibrate a systems biology model of atherosclerosis with adequate granularity to be clinically relevant. Such development may be possible with computational modeling, but given, the multifactorial biology of atherosclerosis, modeling must be based on complete biological networks that capture protein-protein interactions estimated to drive disease progression. Approach and Results: We employed machine intelligence using CTA paired with a molecular assay to determine cohort-level associations and individual patient predictions. Examples of predicted transcripts included ion transporters, cytokine receptors, and a number of microRNAs. Pathway analyses elucidated enrichment of several biological processes relevant to atherosclerosis and plaque pathophysiology. The ability of the models to predict plaque gene expression from CTAs was demonstrated using sequestered patients with transcriptomes of corresponding lesions. We further performed a case study exploring the relationship between biomechanical quantities and plaque morphology, indicating the ability to determine stress and strain from tissue characteristics. Further, we used a uniquely constituted plaque proteomic dataset to create a comprehensive systems biology disease model, which was finally used to simulate responses to different drug categories in individual patients. Individual patient response was simulated for intensive lipid-lowering, anti-inflammatory drugs, anti-diabetic, and combination therapy. Plaque tissue was collected from 18 patients with 6735 proteins at two locations per patient. 113 pathways were identified and included in the systems biology model of endothelial cells, vascular smooth muscle cells, macrophages, lymphocytes, and the integrated intima, altogether spanning 4411 proteins, demonstrating a range of 39-96% plaque instability. Simulations of drug responses varied in patients with initially unstable lesions from high (20%, on combination therapy) to marginal improvement, whereas patients with initially stable plaques showed generally less improvement, but importantly, variation across patients. Conclusion: The results of this thesis show that atherosclerotic plaque phenotyping by multi-scale image analysis of conventional CTA can elucidate the molecular signatures that reflect atherosclerosis. We further showed that calibrated system biology models may be used to simulate drug response in terms of atherosclerotic plaque instability at the individual level, providing a potential strategy for improved personalized management of patients with cardiovascular disease. These results hold promise for optimized and personalized therapy in the prevention of myocardial infarction and ischemic stroke, which warrants further investigations in larger cohorts

A Fluid Structure Interaction Model Of Intracoronary Atherosclerotic Plaque Rupture

Plaque rupture with superimposed thrombosis is the primary cause of acute coronary syndromes of unstable angina, myocardial infarction and sudden death. Although intensive studies in the past decade have shed light on the mechanism that causes unstable atheroma, none has directly addressed the clinical observation that most myocardial infarction (MI) patients have moderate stenoses (less than 50%). Considering the important role the arterial wall compliance and pulsitile blood flow play in atheroma rupture, fluid-structure interaction (FSI) phenomenon has been of interest in recent studies. In this thesis, the impact is investigated numerically of coupled blood flow and structural dynamics on coronary plaque rupture. The objective is to determine a unique index that can be used to characterize plaque rupture potential. The FSI index, developed in this study for the first time derives from the theory of buckling of thin-walled cylinder subjected to radial pressure. Several FSI indices are first defined by normalizing the predicted hemodynamic endothelial shear stress by the structural stresses, specifically, by the maximum principal stress (giving the ratio ), and the Von Mises stress (giving the ratio ). The predicted at the location of maximum (i.e { }) denoted , is then chosen to characterize plaque rupture through systematic investigation of a variety of plaque characteristics and simulated patient conditions. The conditions investigated include varying stenosis levels ranging from 20% to 70%, blood pressure drop ranging from 3125 Pa/m to 9375 Pa/m, fibrous cap thickness ranging from to , lipid pool location ranging from the leading to the trailing edge of plaque, lipid pool volume relative to stenosis volume ranging from 24% to 80%, Calcium volume relative to stenosis volume ranging from 24% to 80% and arterial remodeling. The predicted varies with the stenosis severity and indicates that the plaques investigated are prone to rupture at approximately 40-45% stenosis levels. It predicts that high pressure significantly lowers the threshold stenosis rate for plaque rupture. In addition, the plaque potential to rupture increases for relatively thin fibrous cap, lipid core located near the leading plaque shoulder, and dramatically for relative lipid pool volume above 60%. However, calcium deposit has marginal effect on plaque rupture. Overall, the predicted results are consistent with clinical observations, indicating that the has the potential to characterize plaque rupture when properly established. In the appendix, the unsteady flow in a collapsible tube model of a diseased artery is solved analytically. The novelty of our approach is that the set of governing equations is reduced to a single integro-differential equation in the transient state. The equation was solved using the finite difference method to obtain the pressure and compliant wall behavior. The analytical approach is less computer-intensive than solving the full set of governing equations. The predicted membrane deflection is quite large at low inlet velocity, suggesting possible approach to breakdown in equilibrium. As the transmural pressure increases with wall deflection, bulges appear at the ends of the membrane indicating critical stage of stability, consistent with previous studies. An increase in wall thickness reduces the wall deflection and ultimately results in its collapse. The collapse is due to breakdown in the balance of wall governing equation. An increase in internal pressure is required to maintain membrane stability

Evolution and rupture of vulnerable plaques: a review of mechanical effects

Atherosclerosis occurs as a result of the buildup and infiltration of lipid streaks in artery walls, leading to plaques. Understanding the development of atherosclerosis and plaque vulnerability is of critical importance, since plaque rupture can result in heart attack or stroke. Plaques can be divided into two distinct types: those that rupture (vulnerable) and those that are less likely to rupture (stable). In the last few decades, researchers have been interested in studying the influence of the mechanical effects (blood shear stress, pressure forces, and structural stress) on the plaque formation and rupture processes. In the literature, physiological experimental studies are limited by the complexity of in vivo experiments to study such effects, whereas the numerical approach often uses simplified models compared with realistic conditions, so that no general agreement of the mechanisms responsible for plaque formation has yet been reached. In addition, in a large number of cases, the presence of plaques in arteries is asymptomatic. The prediction of plaque rupture remains a complex question to elucidate, not only because of the interaction of numerous phenomena involved in this process (biological, chemical, and mechanical) but also because of the large time scale on which plaques develop. The purpose of the present article is to review the current mechanical models used to describe the blood flow in arteries in the presence of plaques, as well as reviewing the literature treating the influence of mechanical effects on plaque formation, development, and rupture. Finally, some directions of research, including those being undertaken by the authors, are described

Design and Fatigue Study of Intravascular Coronary Stent using Finite Element Analysis

A thesis presented to the faculty of the College of Science & Technology at Morehead State University in partial fulfillment of the requirements for the Degree of Master of Science by Jared May on March 15, 2012

Cap inflammation leads to higher plaque cap strain and lower cap stress: An MRI-PET/CT-based FSI modeling approach.

Plaque rupture may be triggered by extreme stress/strain conditions. Inflammation is also implicated and can be imaged using novel imaging techniques. The impact of cap inflammation on plaque stress/strain and flow shear stress were investigated. A patient-specific MRI-PET/CT-based modeling approach was used to develop 3D fluid-structure interaction models and investigate the impact of inflammation on plaque stress/strain conditions for better plaque assessment. 18FDG-PET/CT and MRI data were acquired from 4 male patients (average age: 66) to assess plaque characteristics and inflammation. Material stiffness for the fibrous cap was adjusted lower to reflect cap weakening causing by inflammation. Setting stiffness ratio (SR) to be 1.0 (fibrous tissue) for baseline, results for SR=0.5, 0.25, and 0.1 were obtained. Thin cap and hypertension were also considered. Combining results from the 4 patients, mean cap stress from 729 cap nodes was lowered by 25.2% as SR went from 1.0 to 0.1. Mean cap strain value for SR=0.1 was 0.313, 114% higher than that from SR=1.0 model. The thin cap SR=0.1 model had 40% mean cap stress decrease and 81% cap strain increase compared with SR=1.0 model. The hypertension SR=0.1 model had 19.5% cap stress decrease and 98.6% cap strain increase compared with SR=1.0 model. Differences of flow shear stress with 4 different SR values were limited (<10%). Cap inflammation may lead to large cap strain conditions when combined with thin cap and hypertension. Inflammation also led to lower cap stress. This shows the influence of inflammation on stress/strain calculations which are closely related to plaque assessment.This work was supported in part by NIH grants NIH/NIBIB R01 EB004759, NIH/NHLBI R01 HL071021, and National Natural Sciences Foundation of China grant 11672001, 11171030

In Vivo MRI-Based Three-Dimensional Fluid-Structure Interaction Models and Mechanical Image Analysis for Human Carotid Atherosclerotic Plaques

Introduction. Atherosclerotic plaque rupture may occur without warning leading to severe clinical events such as heart attack and stroke. The mechanisms causing plaque rupture are not well understood. It is hypothesized that mechanical forces may play an important role in the plaque rupture process and that image-based computational mechanical analysis may provide useful information for more accurate plaque vulnerability assessment. The objectives of this dissertation are: a) develop in vivo magnetic resonance imaging (MRI)-based 3D computational models with fluid-structure Interactions (FSI) for human atherosclerotic carotid plaques; b) perform mechanical analysis using 3D FSI models to identify critical stress/strain conditions which may be used for possible plaque rupture predictions. Data, Model, and Methods. Histological, ex vivo/ in vivo MRI data of human carotid plaques were provided by the University of Washington Medical School and Washington University Medical School. Blood flow was assumed to be laminar, Newtonian, viscous and incompressible. The Navier-Stokes equations with arbitrary Lagrangian-Eulerian (ALE) formulation were used as the governing equations for the flow model. The vessel and plaque components were assumed to be hyperelastic, isotropic, nearly-incompressible and homogeneous. The nonlinear Mooney-Rivlin model was used to describe the nonlinear properties of the materials with parameter values chosen to match available experimental data. The fully-coupled FSI models were solved by a commercial finite element software ADINA to obtain full 3D flow and stress/strain distributions for analysis. Validation of the computational models and Adina software were provided by comparing computational solutions with analytic solutions and experimental data. Several novel methods were introduced to address some fundamental issues for construction of in vivo MRI-based 3D FSI models: a) an automated MRI segmentation technique using a Bayes theorem with normal probability distribution was implemented to obtain plaque geometry with enclosed components; b) a pre-shrink process was introduced to shrink the in vivo MRI geometry to obtain the no-load shape of the plaque; c) a Volume Component-Fitting Method was introduced to generate a 3D computational mesh for the plaque model with deformable complex geometry, FSI and inclusions; d) a method using MRI data obtained under in vitro pressurized conditions was introduced to determine vessel material properties. Results. The effects of material properties on flow and wall stress/strain behaviors were evaluated. The results indicate that a 100% stiffness increase may decrease maximal values of maximum principal stress (Stress-P1) and maximum principal strain (Strain-P1) by about 20% and 40%, respectively; flow Maximum-Shear-Stress (FMSS) and flow velocity did not show noticeable changes. By comparing ex vivo and in vivo data of 10 plaque samples, the average axial (25%) and inner circumferential (7.9%) shrinkages of the plaques between loaded and unloaded state were obtained. Effects of the shrink-stretch process on plaque stress/strain distributions were demonstrated based on six adjusted 3D FSI models with different shrinkages. Stress-P1 and Strain-P1 increased 349.8% and 249% respectively with 33% axial stretch. The effects of a lipid-rich necrotic core and fibrous cap thickness on structure/flow behaviors were investigated. The mean values of wall Stress-P1 and Strain-P1 from lipid nodes from a ruptured plaque were significantly higher than those from a non-ruptured plaque (112.3 kPa, 0.235 & 80.1 kPa, 0.185), which was 40.2% and 26.8% higher, respectively (p\u3c0.001). High stress/strain concentrations were found at the thin fibrous cap regions. These results indicate that high stress concentrations and thin fibrous cap thickness might be critical indicators for plaque vulnerability. Conclusion. In vivo image-based 3D FSI models and mechanical image analysis may have the potential to provide quantitative risk indicators for plaque vulnerability assessment

Carotid plaque stress analysis by fluid structure interaction based on in-vivo MRI: Implications to plaque vulnerability assessment

This thesis was submitted for the degree of Doctor of Philosophy and awarded by Brunel University, 2010.Stroke is one of the leading causes of death in the world, resulting mostly from the sudden rupture of atherosclerotic plaques. From a biomechanical view, plaque rupture can be considered as a mechanical failure caused by extremely high plaque stress. In this PhD project, we are aiming to predict 3D plaque stress based on in-vivo MRI by using fluid structure interaction (FSI) method, and provide information for plaque rupture risk assessment. Fluid structure interaction was implemented with ANSYS 11.0, followed by a parameter study on fibrous cap thickness and lipid core size with realistic carotid plaque geometry. Twenty patients with carotid plaques imaged by in-vivo MRI were provided in the project. A framework of reconstructing 3D plaque geometry from in-vivo multispectral MRI was designed. The followed reproducibility study on plaque geometry reconstruction procedure and its effect on plaque stress analysis filled the gap in the literature on imaging based plaque stress modeling. The results demonstrated that current MRI technology can provide sufficient information for plaque structure characterization; however stress analysis result is highly affected by MRI resolution and quality. The application of FSI stress analysis to 4 patients with different plaque burdens has showed that the whole procedure from plaque geometry reconstruction to FSI stress analysis was applicable. In the study, plaque geometries from three patients with recent transient ischemic attack were reconstructed by repairing ruptured fibrous cap. The well correlated relationship between local stress concentrations and plaque rupture sites indicated that extremely high plaque stress could be a factor responsible for plaque rupture. Based on the 20 reconstructed carotid plaques from two groups (symptomatic and asymptomatic), fully coupled fluid structure interaction was performed. It was found that there is a significant difference between symptomatic and asymptomatic patients in plaque stress levels, indicating plaque stress could be used as one of the factors for plaque vulnerability assessment. A corresponding plaque morphological feature study showed that plaque stress is significantly affected by fibrous cap thickness, lipid core size and fibrous cap surface irregularities (curvedness). A procedure was proposed for predicting plaque stress by using fibrous cap thickness and curvedness, which requires much less computational time, and has the potential for clinical routine application. The effects of residual stress on plaque stress analysis and arterial wall material property characterization by using in-vivo MRI data were also discussed for patient specific modeling. As the further development, histological study of plaque sample has been combined with conventional plaque stress analysis by assigning material properties to each computational element, based on the data from histological analysis. This method could bridge the gap between biochemistry and biomechanical study of atherosclerosis plaques. In conclusion, extreme stress distributions in the plaque region can be predicted by modern numerical methods, and used for plaque rupture risk assessment, which will be helpful in clinical practice. The combination of plaque MR imaging analysis, computational modelling, and clinical study/ validation would advance our understandings of plaque rupture, prediction of future rupture, and establish new procedures for patient diagnose, management, and treatment.Financial Support was obtained from British Heart Foundation, Brunel Institute for Bioengineering and Brunel Graduate School