Computational modeling of stented coronary arteries

The application of computational fluid dynamics (CFD) and finite element analysis (FEA) has been growing rapidly in the various fields of science and technology. One of the areas of interest is in biomedical engineering. The altered hemodynamics inside the blood vessels plays a key role in the development of the arterial disease called atherosclerosis, which is the major cause of human death worldwide. Atherosclerosis is often treated with the stenting procedure to restore the normal blood flow. A stent is a tubular, flexible structure, usually made of metals, which is driven and expanded in the blocked arteries. Despite the success rate of the stenting procedure, it is often associated with the restenosis (re-narrowing of the artery) process. The presence of non-biological device in the artery causes inflammation or re-growth of atherosclerotic lesions in the treated vessels. Several factors including the design of stents, type of stent expansion, expansion pressure, morphology and composition of vessel wall influence the restenosis process. Therefore, the role of computational studies is crucial in the investigation and optimisation of the factors that influence post-stenting complications. This thesis focuses on the stent-vessel wall interactions followed by the blood flow in the post-stenting stage of stenosed human coronary artery. Hemodynamic and mechanical stresses were analysed in three separate stent-plaque-artery models. Plaque was modeled as a multi-layer (fibrous cap (FC), necrotic core (NC), and fibrosis (F)) and the arterial wall as a single layer domain. CFD/FEA simulations were performed using commercial software packages in several models mimicking the various stages and morphologies of atherosclerosis. The tissue prolapse (TP) of stented vessel wall, the distribution of von Mises stress (VMS) inside various layers of vessel wall, and the wall shear stress (WSS) along the luminal surface of the deformed vessel wall were measured and evaluated. The results revealed the role of the stenosis size, thickness of each layer of atherosclerotic wall, thickness of stent strut, pressure applied for stenosis expansion, and the flow condition in the distribution of stresses. The thicknesses of FC, and NC and the total thickness of plaque are critical in controlling the stresses inside the tissue. A small change in morphology of artery wall can significantly affect the distribution of stresses. In particular, FC is the most sensitive layer to TP and stresses, which could determine plaque’s vulnerability to rupture. The WSS is highly influenced by the deflection of artery, which in turn is dependent on the structural composition of arterial wall layers. Together with the stenosis size, their roles could play a decisive role in controlling the low values of WSS (<0.5 Pa) prone to restenosis. Moreover, the time dependent flow altered the percentage of luminal area with WSS values less than 0.5 Pa at different time instants. The non- Newtonian viscosity model of the blood properties significantly affects the prediction of WSS magnitude. The outcomes of this investigation will help to better understand the roles of the individual layers of atherosclerotic vessels and their risk to provoke restenosis at the post-stenting stage. As a consequence, the implementation of such an approach to assess the post-stented stresses will assist the engineers and clinicians in optimizing the stenting techniques to minimize the occurrence of restenosis.

Flexible Multibody Simulation Approach in the Dynamic Analysis of Bone Strains Activity

The objective of this study is to show that bone strains due to dynamic mechanical loading during physical activity can be analysed using the flexible multibody simulation approach. Strains within the bone tissue play a major role in bone (re)modeling. Based on previous studies, it has been shown that dynamic loading seems to be more important for bone (re)modeling than static loading. The finite element method has been used previously to assess bone strains. However, the finite element method may be limited to static analysis of bone strains due to the expensive computation required for dynamic analysis, especially for a biomechanical system consisting of several bodies. Further, in vivo implementation of strain gauges on the surfaces of bone has been used previously in order to quantify the mechanical loading environment of the skeleton. However, in vivo strain measurement requires invasive methodology, which is challenging and limited to certain regions of superficial bones only, such as the anterior surface of the tibia. In this study, an alternative numerical approach to analyzing in vivo strains, based on the flexible multibody simulation approach, is proposed. In order to investigate the reliability of the proposed approach, three 3-dimensional musculoskeletal models where the right tibia is assumed to be flexible, are used as demonstration examples. The models are employed in a forward dynamics simulation in order to predict the tibial strains during walking on a level exercise. The flexible tibial model is developed using the actual geometry of the subject’s tibia, which is obtained from 3 dimensional reconstruction of Magnetic Resonance Images. Inverse dynamics simulation based on motion capture data obtained from walking at a constant velocity is used to calculate the desired contraction trajectory for each muscle. In the forward dynamics simulation, a proportional derivative servo controller is used to calculate each muscle force required to reproduce the motion, based on the desired muscle contraction trajectory obtained from the inverse dynamics simulation. Experimental measurements are used to verify the models and check the accuracy of the models in replicating the realistic mechanical loading environment measured from the walking test. The predicted strain results by the models show consistency with literature-based in vivo strain measurements. In conclusion, the non-invasive flexible multibody simulation approach may be used as a surrogate for experimental bone strain measurement, and thus be of use in detailed strain estimation of bones in different applications. Consequently, the information obtained from the present approach might be useful in clinical applications, including optimizing implant design and devising exercises to prevent bone fragility, accelerate fracture healing and reduce osteoporotic bone loss.

Modeling the transport phenomena within arterial wall: porous media approach

The transport of macromolecules, such as low-density lipoprotein (LDL), and their accumulation in the layers of the arterial wall play a critical role in the creation and development of atherosclerosis. Atherosclerosis is a disease of large arteries e.g., the aorta, coronary, carotid, and other proximal arteries that involves a distinctive accumulation of LDL and other lipid-bearing materials in the arterial wall. Over time, plaque hardens and narrows the arteries. The flow of oxygen-rich blood to organs and other parts of the body is reduced. This can lead to serious problems, including heart attack, stroke, or even death. It has been proven that the accumulation of macromolecules in the arterial wall depends not only on the ease with which materials enter the wall, but also on the hindrance to the passage of materials out of the wall posed by underlying layers. Therefore, attention was drawn to the fact that the wall structure of large arteries is different than other vessels which are disease-resistant. Atherosclerosis tends to be localized in regions of curvature and branching in arteries where fluid shear stress (shear rate) and other fluid mechanical characteristics deviate from their normal spatial and temporal distribution patterns in straight vessels. On the other hand, the smooth muscle cells (SMCs) residing in the media layer of the arterial wall respond to mechanical stimuli, such as shear stress. Shear stress may affect SMC proliferation and migration from the media layer to intima. This occurs in atherosclerosis and intimal hyperplasia. The study of blood flow and other body fluids and of heat transport through the arterial wall is one of the advanced applications of porous media in recent years. The arterial wall may be modeled in both macroscopic (as a continuous porous medium) and microscopic scales (as a heterogeneous porous medium). In the present study, the governing equations of mass, heat and momentum transport have been solved for different species and interstitial fluid within the arterial wall by means of computational fluid dynamics (CFD). Simulation models are based on the finite element (FE) and finite volume (FV) methods. The wall structure has been modeled by assuming the wall layers as porous media with different properties. In order to study the heat transport through human tissues, the simulations have been carried out for a non-homogeneous model of porous media. The tissue is composed of blood vessels, cells, and an interstitium. The interstitium consists of interstitial fluid and extracellular fibers. Numerical simulations are performed in a two-dimensional (2D) model to realize the effect of the shape and configuration of the discrete phase on the convective and conductive features of heat transfer, e.g. the interstitium of biological tissues. On the other hand, the governing equations of momentum and mass transport have been solved in the heterogeneous porous media model of the media layer, which has a major role in the transport and accumulation of solutes across the arterial wall. The transport of Adenosine 5´-triphosphate (ATP) is simulated across the media layer as a benchmark to observe how SMCs affect on the species mass transport. In addition, the transport of interstitial fluid has been simulated while the deformation of the media layer (due to high blood pressure) and its constituents such as SMCs are also involved in the model. In this context, the effect of pressure variation on shear stress is investigated over SMCs induced by the interstitial flow both in 2D and three-dimensional (3D) geometries for the media layer. The influence of hypertension (high pressure) on the transport of lowdensity lipoprotein (LDL) through deformable arterial wall layers is also studied. This is due to the pressure-driven convective flow across the arterial wall. The intima and media layers are assumed as homogeneous porous media. The results of the present study reveal that ATP concentration over the surface of SMCs and within the bulk of the media layer is significantly dependent on the distribution of cells. Moreover, the shear stress magnitude and distribution over the SMC surface are affected by transmural pressure and the deformation of the media layer of the aorta wall. This work reflects the fact that the second or even subsequent layers of SMCs may bear shear stresses of the same order of magnitude as the first layer does if cells are arranged in an arbitrary manner. This study has brought new insights into the simulation of the arterial wall, as the previous simplifications have been ignored. The configurations of SMCs used here with elliptic cross sections of SMCs closely resemble the physiological conditions of cells. Moreover, the deformation of SMCs with high transmural pressure which follows the media layer compaction has been studied for the first time. On the other hand, results demonstrate that LDL concentration through the intima and media layers changes significantly as wall layers compress with transmural pressure. It was also noticed that the fraction of leaky junctions across the endothelial cells and the area fraction of fenestral pores over the internal elastic lamina affect the LDL distribution dramatically through the thoracic aorta wall. The simulation techniques introduced in this work can also trigger new ideas for simulating porous media involved in any biomedical, biomechanical, chemical, and environmental engineering applications.

Kokeellisen moodianalyysin hyödyntäminen sääriluun lujuuden määrittämisessä

Diplomityö tehtiin Lappeenrannan teknillisen yliopiston konetekniikan laitokselle. Diplomityö on osa teknillisen yliopiston biomekaanista tutkimusta, jonka tarkoituksena on mallintaa ihmisen tuki- ja liikuntaelimistön toimintaa. Työssä pyrittiin selvittämään, voitaisiinko sääriluuhun kohdistetun mekaanisen herätteen aiheuttamaa värähtelyvastetta analysoimalla saada tietoa luun ominaistaajuuksista ja lujuudesta. Tietoa voitaisiin käyttää esimerkiksi ostoporoosiriskin arvioinnissa sekä ihmiskehon osien toimintaa kuvaavien simulointimallien verifioinnissa. Mittauslaitteistona käytettiin Brüel & Kjær-moodianalyysilaitteistoa. Laitteistokokonaisuuteen kuuluivat herätevasara, elektromagneettinen täristin, voima-anturi, kaksi kiihtyvyysmitta-anturia sekä PulseLab 2.0 –ohjelmistolla varustettu PC-laitteisto. Tulosten jatkoanalyysi suoritettiin MathWorks yhtiön MatLab v 4.0 -ohjelmistolla. Työssä esitellyn mittaustavan ja -laitteiston todettiin soveltuvan sääriluun värähtelyvasteen mittaamiseen. Mittaustulokset eri mittauskertojen välillä samalla henkilöllä ovat yhtenevät. Tutkimuksen tulosten perusteella ei voida osoittaa luun värähtelyvasteen ja lujuuden välistä suoraa korrelaatiota.

CT and PET/CT Hybrid Imaging of Coronary Artery Disease

Coronary artery disease (CAD) is a chronic process that evolves over decades and may culminate in myocardial infarction (MI). While invasive coronary angiography (ICA) is still considered the gold standard of imaging CAD, non-invasive assessment of both the vascular anatomy and myocardial perfusion has become an intriguing alternative. In particular, computed tomography (CT) and positron emission tomography (PET) form an attractive combination for such studies. Increased radiation dose is, however, a concern. Our aim in the current thesis was to test novel CT and PET techniques alone and in hybrid setting in the detection and assessment of CAD in clinical patients. Along with diagnostic accuracy, methods for the reduction of the radiation dose was an important target. The study investigating the coronary arteries of patients with atrial fibrillation (AF) showed that CAD may be an important etiology of AF because a high prevalence of CAD was demonstrated within AF patients. In patients with suspected CAD, we demonstrated that a sequential, prospectively ECG-triggered CT technique was applicable to nearly 9/10 clinical patients and the radiation dose was over 60% lower than with spiral CT. To detect the functional significance of obstructive CAD, a novel software for perfusion quantification, CarimasTM, showed high reproducibility with 15O-labelled water in PET, supporting feasibility and good clinical accuracy. In a larger cohort of 107 patients with moderate 30-70% pre-test probability of CAD, hybrid PET/CT was shown to be a powerful diagnostic method in the assessment of CAD with diagnostic accuracy comparable to that of invasive angiography and fractional flow reserve (FFR) measurements. A hybrid study may be performed with a reasonable radiation dose in a vast majority of the cases, improving the performance of stand-alone PET and CT angiography, particularly when the absolute quantification of the perfusion is employed. These results can be applied into clinical practice and will be useful for daily clinical diagnosis of CAD.

Application of Computational Fluid Dynamics in Modeling Blood Flow in Human Thoracic Aorta

The aim of this study was to simulate blood flow in thoracic human aorta and understand the role of flow dynamics in the initialization and localization of atherosclerotic plaque in human thoracic aorta. The blood flow dynamics in idealized and realistic models of human thoracic aorta were numerically simulated in three idealized and two realistic thoracic aorta models. The idealized models of thoracic aorta were reconstructed with measurements available from literature, and the realistic models of thoracic aorta were constructed by image processing Computed Tomographic (CT) images. The CT images were made available by South Karelia Central Hospital in Lappeenranta. The reconstruction of thoracic aorta consisted of operations, such as contrast adjustment, image segmentations, and 3D surface rendering. Additional design operations were performed to make the aorta model compatible for the numerical method based computer code. The image processing and design operations were performed with specialized medical image processing software. Pulsatile pressure and velocity boundary conditions were deployed as inlet boundary conditions. The blood flow was assumed homogeneous and incompressible. The blood was assumed to be a Newtonian fluid. The simulations with idealized models of thoracic aorta were carried out with Finite Element Method based computer code, while the simulations with realistic models of thoracic aorta were carried out with Finite Volume Method based computer code. Simulations were carried out for four cardiac cycles. The distribution of flow, pressure and Wall Shear Stress (WSS) observed during the fourth cardiac cycle were extensively analyzed. The aim of carrying out the simulations with idealized model was to get an estimate of flow dynamics in a realistic aorta model. The motive behind the choice of three aorta models with distinct features was to understand the dependence of flow dynamics on aorta anatomy. Highly disturbed and nonuniform distribution of velocity and WSS was observed in aortic arch, near brachiocephalic, left common artery, and left subclavian artery. On the other hand, the WSS profiles at the roots of branches show significant differences with geometry variation of aorta and branches. The comparison of instantaneous WSS profiles revealed that the model with straight branching arteries had relatively lower WSS compared to that in the aorta model with curved branches. In addition to this, significant differences were observed in the spatial and temporal profiles of WSS, flow, and pressure. The study with idealized model was extended to study blood flow in thoracic aorta under the effects of hypertension and hypotension. One of the idealized aorta models was modified along with the boundary conditions to mimic the thoracic aorta under the effects of hypertension and hypotension. The results of simulations with realistic models extracted from CT scans demonstrated more realistic flow dynamics than that in the idealized models. During systole, the velocity in ascending aorta was skewed towards the outer wall of aortic arch. The flow develops secondary flow patterns as it moves downstream towards aortic arch. Unlike idealized models, the distribution of flow was nonplanar and heavily guided by the artery anatomy. Flow cavitation was observed in the aorta model which was imaged giving longer branches. This could not be properly observed in the model with imaging containing a shorter length for aortic branches. The flow circulation was also observed in the inner wall of the aortic arch. However, during the diastole, the flow profiles were almost flat and regular due the acceleration of flow at the inlet. The flow profiles were weakly turbulent during the flow reversal. The complex flow patterns caused a non-uniform distribution of WSS. High WSS was distributed at the junction of branches and aortic arch. Low WSS was distributed at the proximal part of the junction, while intermedium WSS was distributed in the distal part of the junction. The pulsatile nature of the inflow caused oscillating WSS at the branch entry region and inner curvature of aortic arch. Based on the WSS distribution in the realistic model, one of the aorta models was altered to induce artificial atherosclerotic plaque at the branch entry region and inner curvature of aortic arch. Atherosclerotic plaque causing 50% blockage of lumen was introduced in brachiocephalic artery, common carotid artery, left subclavian artery, and aortic arch. The aim of this part of the study was first to study the effect of stenosis on flow and WSS distribution, understand the effect of shape of atherosclerotic plaque on flow and WSS distribution, and finally to investigate the effect of lumen blockage severity on flow and WSS distributions. The results revealed that the distribution of WSS is significantly affected by plaque with mere 50% stenosis. The asymmetric shape of stenosis causes higher WSS in branching arteries than in the cases with symmetric plaque. The flow dynamics within thoracic aorta models has been extensively studied and reported here. The effects of pressure and arterial anatomy on the flow dynamic were investigated. The distribution of complex flow and WSS is correlated with the localization of atherosclerosis. With the available results we can conclude that the thoracic aorta, with complex anatomy is the most vulnerable artery for the localization and development of atherosclerosis. The flow dynamics and arterial anatomy play a role in the localization of atherosclerosis. The patient specific image based models can be used to diagnose the locations in the aorta vulnerable to the development of arterial diseases such as atherosclerosis.

Generation and Characterization of Knockout Mice and Analysis of Patient Material Demonstrates a Role for Tissue Inhibitor of Metalloproteinases 4 (Timp4) in the Cardiovascular System and Tumor Metastasis

This thesis focuses on tissue inhibitor of metalloproteinases 4 (TIMP4) which is the newest member of a small gene and protein family of four closely related endogenous inhibitors of extracellular matrix (ECM) degrading enzymes. Existing data on TIMP4 suggested that it exhibits a more restricted expression pattern than the other TIMPs with high expression levels in heart, brain, ovary and skeletal muscle. These observations and the fact that the ECM is of special importance to provide the cardiovascular system with structural strength combined with elasticity and distensibility, prompted the present molecular biologic investigation on TIMP4. In the first part of the study the murine Timp4 gene was cloned and characterized in detail. The structure of murine Timp4 genomic locus resembles that in other species and of the other Timps. The highest Timp4 expression was detected in heart, ovary and brain. As the expression pattern of Timp4 gives only limited information about its role in physiology and pathology, Timp4 knockout mice were generated next. The analysis of Timp4 knockout mice revealed that Timp4 deficiency has no obvious effect on the development, growth or fertility of mice. Therefore, Timp4 deficient mice were challenged using available cardiovascular models, i.e. experimental cardiac pressure overload and myocardial infarction. In the former model, Timp4 deficiency was found to be compensated by Timp2 overexpression, whereas in the myocardial infarct model, Timp4 deficiency resulted in increased mortality due to increased susceptibility for cardiac rupture. In the wound healing model, Timp4 deficiency was shown to result in transient retardation of re-epithelialization of cutaneous wounds. Melanoma tumor growth was similar in Timp4 deficient and control mice. Despite of this, lung metastasis of melanoma cells was significantly increased in Timp4 null mice. In an attempt to translate the current findings to patient material, TIMP4 expression was studied in human specimens representing different inflammatory cardiovascular pathologies, i.e. giant cell arteritis, atherosclerotic coronary arteries and heart allografts exhibiting signs of chronic rejection. The results showed that cardiovascular expression of TIMP4 is elevated particularly in areas exhibiting inflammation. The results of the present studies suggest that TIMP4 has a special role in the regulation of tissue repair processes in the heart, and also in healing wounds and metastases. Furthermore, evidence is provided suggesting the usefulness of TIMP4 as a novel systemic marker for vascular inflammation.

The effects of maternal hyperglycemia on human umbilical vascular gene expression and neonatal rat lung development

Background: Maternal diabetes affects many fetal organ systems, including the vasculature and the lungs. The offspring of diabetic mothers have respiratory adaptation problems after birth. The mechanisms are multifactorial and the effects are prolonged during the postnatal period. An increasing incidence of diabetic pregnancies accentuates the importance of identifying the pathological mechanisms, which cause the metabolic and genetic changes that occur in offspring, born to diabetic mothers. Aims and methods: The aim of this thesis was to determine changes both in human umbilical cord exposed to maternal type 1 diabetes and in neonatal rat lungs after streptozotocin-induced maternal hyperglycemia, during pregnancy. Rat lungs were used as a model for the potential disease mechanisms. Gene expression alterations were determined in human umbilical cords at birth and in rat pup lungs at two week of age. During the first two postnatal weeks, rat lung development was studied morphologically and histologically. Further, the effect of postnatal hyperoxia on hyperglycemia-primed rat lungs was investigated at one week of age to mimic the clinical situation of supplemental oxygen treatment. Results: In the umbilical cord, maternal diabetes had a major negative effect on the expression of genes involved in blood vessel development. The genes regulating vascular tone were also affected. In neonatal rat lungs, intrauterine hyperglycemia had a prolonged effect on gene expression during late alveolarization. The most affected pathway was the upregulation of extracellular matrix proteins. Newborn rat lungs exposed to intrauterine hyperglycemia had thinner saccular walls without changes in airspace size, a smaller relative lung weight and lung total tissue area, and increased cellular apoptosis and proliferation compared to control lungs, possibly reflecting an aberrant maturational adaptation. At one and two weeks of age, cell proliferation and secondary crest formation were accelerated in hyperglycemia-exposed lungs. Postnatal hyperoxic exposure, alone caused arrested alveolarization with thin-walled and enlarged alveoli. In contrast, the dual exposure of intrauterine hyperglycemia and postnatal hyperoxia resulted in the phenotype of thick septa together with arrested alveolarization and decreased number of small pulmonary arteries. Conclusions: Maternal diabetic environment seems to alter the umbilical cord gene expression profile of the regulation of vascular development and function. Fetal hyperglycemia may additionally affect the genetic regulation of the postnatal lung development and may actually induce prolonged structural alterations in neonatal lungs together with a modifying effect on the deleterious pulmonary exposure of postnatal hyperoxia. This, combined with the novel human umbilical cord gene data could serve as stepping stones for future therapies to curb developmental aberrations.

Intracranial Nimodipine Implant: Feasibility and Implications for the Treatment of Subarachnoid Hemorrhage – A Pre-Clinical Study

Intracranial aneurysmal subarachnoid hemorrhage (aSAH) is a life-threatening condition requiring immediate neurocritical care. A ruptured aneurysm must be isolated from arterial circulation to prevent rebleeding. Open surgical clipping of the neck of the aneurysm or intra-arterial filling of the aneurysm sack with platinum coils are major treatment strategies in an acute phase. About 40% of the patients suffering from aSAH die within a year of the bleeding despite the intensive treatment. After aSAH, the patient may develop a serious complication called vasospasm. Major risk for the vasospasm takes place at days 5–14 after the primary bleeding. In vasospasm, cerebral arteries contract uncontrollably causing brain ischemia that may lead to death. Nimodipine (NDP) is used to treat of vasospasm and it is administrated intravenously or orally every four hours for 21 days. NDP treatment has been scientifically proven to improve patients’ clinical outcome. The therapeutic effect of L-type calcium channel blocker NDP is due to the ability to dilate cerebral arteries. In addition to vasodilatation, recent research has shown the pleiotropic effect of NDP such as inhibition of neuronal apoptosis and inhibition of microthrombi formation. Indeed, NDP inhibits cortical spreading ischemia. Knowledge of the pathophysiology of the vasospasm has evolved in recent years to a complex entity of early brain injury, secondary injuries and cortical spreading ischemia, instead of being pure intracranial vessel spasm. High NDP levels are beneficial since they protect neurons and inhibit the cortical spreading ischemia. One of the drawbacks of the intravenous or oral administration of NPD is systemic hypotension, which is harmful particularly when the brain is injured. Maximizing the beneficial effects and avoiding systemic hypotension of NDP, we developed a sustained release biodegradable NDP implant that was surgically positioned in the basal cistern of animal models (dog and pig). Higher concentrations were achieved locally and lower concentrations systemically. Using this treatment approach in humans, it may be possible to reduce incidence of harmful hypotension and potentiate beneficial effects of NDP on neurons. Intracellular calcium regulation has a pivotal role in brain plasticity. NDP blocks L-type calcium channels in neurons, substantially decreasing intracellular calcium levels. Thus, we were interested in how NDP affects brain plasticity and tested the hypothesis in a mouse model. We found that NDP activates Brain-derived neurotrophic factor (BDNF) receptor TrkB and its downstream signaling in a reminiscent of antidepressant drugs. In contrast to antidepressant drugs, NDP activates Akt, a major survival-promoting factor. Our group’s previous findings demonstrate that long-term antidepressant treatment reactivates developmental-type of plasticity mechanisms in the adult brain, which allows the remodeling of neuronal networks if combined with appropriate rehabilitation. It seems that NDP has antidepressant-like properties and it is able to induce neuronal plasticity. In general, drug induced neuronal plasticity has a huge potential in neurorehabilitation and more studies are warranted.