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.

The effects of positive end-expiratory pressure on respiratory system mechanics and hemodynamics in postoperative cardiac surgery patients

We prospectively evaluated the effects of positive end-expiratory pressure (PEEP) on the respiratory mechanical properties and hemodynamics of 10 postoperative adult cardiac patients undergoing mechanical ventilation while still anesthetized and paralyzed. The respiratory mechanics was evaluated by the inflation inspiratory occlusion method and hemodynamics by conventional methods. Each patient was randomized to a different level of PEEP (5, 10 and 15 cmH2O), while zero end-expiratory pressure (ZEEP) was established as control. PEEP of 15-min duration was applied at 20-min intervals. The frequency dependence of resistance and the viscoelastic properties and elastance of the respiratory system were evaluated together with hemodynamic and respiratory indexes. We observed a significant decrease in total airway resistance (13.12 ± 0.79 cmH2O l-1 s-1 at ZEEP, 11.94 ± 0.55 cmH2O l-1 s-1 (P<0.0197) at 5 cmH2O of PEEP, 11.42 ± 0.71 cmH2O l-1 s-1 (P<0.0255) at 10 cmH2O of PEEP, and 10.32 ± 0.57 cmH2O l-1 s-1 (P<0.0002) at 15 cmH2O of PEEP). The elastance (Ers; cmH2O/l) was not significantly modified by PEEP from zero (23.49 ± 1.21) to 5 cmH2O (21.89 ± 0.70). However, a significant decrease (P<0.0003) at 10 cmH2O PEEP (18.86 ± 1.13), as well as (P<0.0001) at 15 cmH2O (18.41 ± 0.82) was observed after PEEP application. Volume dependence of viscoelastic properties showed a slight but not significant tendency to increase with PEEP. The significant decreases in cardiac index (l min-1 m-2) due to PEEP increments (3.90 ± 0.22 at ZEEP, 3.43 ± 0.17 (P<0.0260) at 5 cmH2O of PEEP, 3.31 ± 0.22 (P<0.0260) at 10 cmH2O of PEEP, and 3.10 ± 0.22 (P<0.0113) at 15 cmH2O of PEEP) were compensated for by an increase in arterial oxygen content owing to shunt fraction reduction (%) from 22.26 ± 2.28 at ZEEP to 11.66 ± 1.24 at PEEP of 15 cmH2O (P<0.0007). We conclude that increments in PEEP resulted in a reduction of both airway resistance and respiratory elastance. These results could reflect improvement in respiratory mechanics. However, due to possible hemodynamic instability, PEEP should be carefully applied to postoperative cardiac patients.

Gap junctions in cells of the immune system: structure, regulation and possible functional roles

Gap junction channels are sites of cytoplasmic communication between contacting cells. In vertebrates, they consist of protein subunits denoted connexins (Cxs) which are encoded by a gene family. According to their Cx composition, gap junction channels show different gating and permeability properties that define which ions and small molecules permeate them. Differences in Cx primary sequences suggest that channels composed of different Cxs are regulated differentially by intracellular pathways under specific physiological conditions. Functional roles of gap junction channels could be defined by the relative importance of permeant substances, resulting in coordination of electrical and/or metabolic cellular responses. Cells of the native and specific immune systems establish transient homo- and heterocellular contacts at various steps of the immune response. Morphological and functional studies reported during the last three decades have revealed that many intercellular contacts between cells in the immune response present gap junctions or "gap junction-like" structures. Partial characterization of the molecular composition of some of these plasma membrane structures and regulatory mechanisms that control them have been published recently. Studies designed to elucidate their physiological roles suggest that they might permit coordination of cellular events which favor the effective and timely response of the immune system.

Role of sympathetic nervous system and neuropeptides in obesity hypertension

Obesity is the most common cause of human essential hypertension in most industrialized countries. Although the precise mechanisms of obesity hypertension are not fully understood, considerable evidence suggests that excess renal sodium reabsorption and a hypertensive shift of pressure natriuresis play a major role. Sympathetic activation appears to mediate at least part of the obesity-induced sodium retention and hypertension since adrenergic blockade or renal denervation markedly attenuates these changes. Recent observations suggest that leptin and its multiple interactions with neuropeptides in the hypothalamus may link excess weight gain with increased sympathetic activity. Leptin is produced mainly in adipocytes and is believed to regulate energy balance by acting on the hypothalamus to reduce food intake and to increase energy expenditure via sympathetic activation. Short-term administration of leptin into the cerebral ventricles increases renal sympathetic activity, and long-term leptin infusion at rates that mimic plasma concentrations found in obesity raises arterial pressure and heart rate via adrenergic activation in non-obese rodents. Transgenic mice overexpressing leptin also develop hypertension. Acute studies suggest that the renal sympathetic effects of leptin may depend on interactions with other neurochemical pathways in the hypothalamus, including the melanocortin-4 receptor (MC4-R). However, the role of this pathway in mediating the long-term effects of leptin on blood pressure is unclear. Also, it is uncertain whether there is resistance to the chronic renal sympathetic and blood pressure effects of leptin in obese subjects. In addition, leptin also has other cardiovascular and renal actions, such as stimulation of nitric oxide formation and improvement of insulin sensitivity, which may tend to reduce blood pressure in some conditions. Although the role of these mechanisms in human obesity has not been elucidated, this remains a fruitful area for further investigation, especially in view of the current "epidemic" of obesity in most industrialized countries.

Towards understanding the kallikrein-kinin system: insights from measurement of kinin peptides

The kallikrein-kinin system is complex, with several bioactive peptides that are formed in many different compartments. Kinin peptides are implicated in many physiological and pathological processes including the regulation of blood pressure and sodium homeostasis, inflammatory processes, and the cardioprotective effects of preconditioning. We established a methodology for the measurement of individual kinin peptides in order to study the function of the kallikrein-kinin system. The levels of kinin peptides in tissues were higher than in blood, confirming the primary tissue localization of the kallikrein-kinin system. Moreover, the separate measurement of bradykinin and kallidin peptides in man demonstrated the differential regulation of the plasma and tissue kallikrein-kinin systems, respectively. Kinin peptide levels were increased in the heart of rats with myocardial infarction, in tissues of diabetic and spontaneously hypertensive rats, and in urine of patients with interstitial cystitis, suggesting a role for kinin peptides in the pathogenesis of these conditions. By contrast, blood levels of kallidin, but not bradykinin, peptides were suppressed in patients with severe cardiac failure, suggesting that the activity of the tissue kallikrein-kinin system may be suppressed in this condition. Both angiotensin converting enzyme (ACE) and neutral endopeptidase (NEP) inhibitors increased bradykinin peptide levels. ACE and NEP inhibitors had different effects on kinin peptide levels in blood, urine, and tissues, which may be accounted for by the differential contributions of ACE and NEP to kinin peptide metabolism in the multiple compartments in which kinin peptide generation occurs. Measurement of the levels of individual kinin peptides has given important information about the operation of the kallikrein-kinin system and its role in physiology and disease states.

CFD modeling of dispersion water feed in wastewater cleaning application

Fluid particle breakup and coalescence are important phenomena in a number of industrial flow systems. This study deals with a gas-liquid bubbly flow in one wastewater cleaning application. Three-dimensional geometric model of a dispersion water system was created in ANSYS CFD meshing software. Then, numerical study of the system was carried out by means of unsteady simulations performed in ANSYS FLUENT CFD software. Single-phase water flow case was setup to calculate the entire flow field using the RNG k-epsilon turbulence model based on the Reynolds-averaged Navier-Stokes (RANS) equations. Bubbly flow case was based on a computational fluid dynamics - population balance model (CFD-PBM) coupled approach. Bubble breakup and coalescence were considered to determine the evolution of the bubble size distribution. Obtained results are considered as steps toward optimization of the cleaning process and will be analyzed in order to make the process more efficient.

A study of the interaction between Helicobacter pylori and components of the human fibrinolytic system

The interaction of plasminogen, tissue plasminogen activator (t-PA) and urokinase with a clinical strain of Helicobacter pylori was studied. Plasminogen bound to the surface of H. pylori cells in a concentration-dependent manner and could be activated to the enzymatic form, plasmin, by t-PA. Affinity chromatography assays revealed a plasminogen-binding protein of 58.9 kDa in water extracts of surface proteins. Surface-associated plasmin activity, detected with the chromogenic substrate CBS 00.65, was observed only when plasminogen and an exogenous activator were added to the cell suspension. The two physiologic plasminogen activators, t-PA and urokinase, were also shown to bind to and remain active on the surface of bacterial cells. epsilon-Aminocaproic acid caused partial inhibition of t-PA binding, suggesting that the kringle 2 structure of this activator is involved in the interaction with surface receptors. The activation of plasminogen by t-PA, but not urokinase, strongly depended on the presence of cells and a 25-fold enhancer effect on the initial velocity of activation by t-PA compared to urokinase was established. Furthermore, a relationship between cell concentration and the initial velocity of activation was demonstrated. These findings support the concept that plasminogen activation by t-PA on the bacterial surface is a surface-dependent reaction which offers catalytic advantages.

Adaptive Markov chain Monte Carlo and Bayesian filtering for state space models

This thesis is concerned with the state and parameter estimation in state space models. The estimation of states and parameters is an important task when mathematical modeling is applied to many different application areas such as the global positioning systems, target tracking, navigation, brain imaging, spread of infectious diseases, biological processes, telecommunications, audio signal processing, stochastic optimal control, machine learning, and physical systems. In Bayesian settings, the estimation of states or parameters amounts to computation of the posterior probability density function. Except for a very restricted number of models, it is impossible to compute this density function in a closed form. Hence, we need approximation methods. A state estimation problem involves estimating the states (latent variables) that are not directly observed in the output of the system. In this thesis, we use the Kalman filter, extended Kalman filter, Gauss–Hermite filters, and particle filters to estimate the states based on available measurements. Among these filters, particle filters are numerical methods for approximating the filtering distributions of non-linear non-Gaussian state space models via Monte Carlo. The performance of a particle filter heavily depends on the chosen importance distribution. For instance, inappropriate choice of the importance distribution can lead to the failure of convergence of the particle filter algorithm. In this thesis, we analyze the theoretical Lᵖ particle filter convergence with general importance distributions, where p ≥2 is an integer. A parameter estimation problem is considered with inferring the model parameters from measurements. For high-dimensional complex models, estimation of parameters can be done by Markov chain Monte Carlo (MCMC) methods. In its operation, the MCMC method requires the unnormalized posterior distribution of the parameters and a proposal distribution. In this thesis, we show how the posterior density function of the parameters of a state space model can be computed by filtering based methods, where the states are integrated out. This type of computation is then applied to estimate parameters of stochastic differential equations. Furthermore, we compute the partial derivatives of the log-posterior density function and use the hybrid Monte Carlo and scaled conjugate gradient methods to infer the parameters of stochastic differential equations. The computational efficiency of MCMC methods is highly depend on the chosen proposal distribution. A commonly used proposal distribution is Gaussian. In this kind of proposal, the covariance matrix must be well tuned. To tune it, adaptive MCMC methods can be used. In this thesis, we propose a new way of updating the covariance matrix using the variational Bayesian adaptive Kalman filter algorithm.

Evidence of muscarinic acetylcholine receptors in the retinal centrifugal system of the chick

In this study we characterize the presence of muscarinic acetylcholine receptors (mAChR) in the isthmo-optic nucleus (ION) of chicks by immunohistochemistry with the M35 antibody. Some M35-immunoreactive fibers were observed emerging from the retinal optic nerve insertion, suggesting that they could be centrifugal fibers. Indeed, intraocular injections of cholera toxin B (CTb), a retrograde tracer, and double-labeling with M35 and CTb in the ION confirmed this hypothesis. The presence of M35-immunoreactive cells and the possible mAChR expression in ION and ectopic neuron cells in the chick brain strongly suggest the existence of such a cholinergic system in this nucleus and that acetylcholine release from amacrine cells may mediate interactions between retinal cells and ION terminals.

The neuroimmune-endocrine axis: pathophysiological implications for the central nervous system cytokines and hypothalamus-pituitary-adrenal hormone dynamics

Cytokines are molecules that were initially discovered in the immune system as mediators of communication between various types of immune cells. However, it soon became evident that cytokines exert profound effects on key functions of the central nervous system, such as food intake, fever, neuroendocrine regulation, long-term potentiation, and behavior. In the 80's and 90's our group and others discovered that the genes encoding various cytokines and their receptors are expressed in vascular, glial, and neuronal structures of the adult brain. Most cytokines act through cell surface receptors that have one transmembrane domain and which transduce a signal through the JAK/STAT pathway. Of particular physiological and pathophysiological relevance is the fact that cytokines are potent regulators of hypothalamic neuropeptidergic systems that maintain neuroendocrine homeostasis and which regulate the body's response to stress. The mechanisms by which cytokine signaling affects the function of stress-related neuroendocrine systems are reviewed in this article.

A comparative study of aging of the elastic fiber system of the diaphragm and the rectus abdominis muscles in rats

In the present study the age-related changes of the striated muscle elastic fiber system were investigated in the diaphragm and rectus abdominis muscles of 1-, 4-, 8- and 18-month-old rats. The activation patterns of these muscles differ in that the diaphragm is regularly mobilized tens of times every minute during the entire life of the animal whereas the rectus abdominis, although mobilized in respiration, is much less and more irregularly activated. The elastic fibers were stained by the Verhoeff technique for mature elastic fibers. Weigert stain was used to stain mature and elaunin elastic fibers, and Weigert-oxone to stain mature, elaunin and oxytalan elastic fibers. The density of mature and elaunin elastic fibers showed a progressive increase with age, whereas the amount of oxytalan elastic fibers decreased in both diaphragm and rectus abdominis muscles and their muscular fascias. These age-related quantitative and structural changes of the elastic fiber system may reduce the viscoelastic properties of the diaphragm and rectus abdominis muscles, which may compromise the transmission of tensile muscle strength to the tendons and may affect maximum total strength.

Determination of macromolecular exchange and PO2 in the microcirculation: a simple system for in vivo fluorescence and phosphorescence videomicroscopy

We have developed a system with two epi-illumination sources, a DC-regulated lamp for transillumination and mechanical switches for rapid shift of illumination and detection of defined areas (250-750 µm²) by fluorescence and phosphorescence videomicroscopy. The system permits investigation of standard microvascular parameters, vascular permeability as well as intra- and extravascular PO2 by phosphorescence quenching of Pd-meso-tetra (4-carboxyphenyl) porphine (PORPH). A Pechan prism was used to position a defined region over the photomultiplier and TV camera. In order to validate the system for in vivo use, in vitro tests were performed with probes at concentrations that can be found in microvascular studies. Extensive in vitro evaluations were performed by filling glass capillaries with solutions of various concentrations of FITC-dextran (diluted in blood and in saline) mixed with different amounts of PORPH. Fluorescence intensity and phosphorescence decay were determined for each mixture. FITC-dextran solutions without PORPH and PORPH solutions without FITC-dextran were used as references. Phosphorescence decay curves were relatively unaffected by the presence of FITC-dextran at all concentrations tested (0.1 µg/ml to 5 mg/ml). Likewise, fluorescence determinations were performed in the presence of PORPH (0.05 to 0.5 mg/ml). The system was successfully used to study macromolecular extravasation and PO2 in the rat mesentery circulation under controlled conditions and during ischemia-reperfusion.

Influence of acute renal failure on the mononuclear phagocytic system

Several studies show the ability of macrophages to remove particles injected into the bloodstream. This function seems to be increased in the presence of acute renal failure. The objective of the present study was to assess the phagocytic function of the main organs (spleen, liver and lung) of the mononuclear phagocytic system in renal and postrenal failures. Fifteen rats (250-350 g) were divided into three groups (N = 5): group I - control; group II - ligature of both ureters, and group III - bilateral nephrectomy. On the third postoperative day, all animals received an iv injection of 1 ml/kg 99mTc sulfur colloid. Blood samples were collected for the assessment of plasma urea, creatinine, sodium, and potassium concentrations and arterial gasometry. Samples of liver, spleen, lung and blood clots were obtained and radioactivity was measured. Samples of liver, spleen, lung and kidney were prepared for routine histopathological analysis. Plasma urea, creatinine and potassium concentrations in groups II and III were higher than in group I (P<0.05). Plasma sodium concentrations in groups II and III were lower than in group I (P<0.05). Compensated metabolic acidosis was observed in the presence of postrenal failure. Group II animals showed a lower level of radioactivity in the spleen (0.98) and lung (2.63), and a higher level in the liver (105.51) than control. Group III animals showed a lower level of radioactivity in the spleen (11.94) and a higher level in the liver (61.80), lung (11.30) and blood clot (5.13) than control. In groups II and III liver steatosis and bronchopneumonia were observed. Renal and postrenal failures seem to interfere with blood clearance by the mononuclear phagocytic system.

Frequent occurrence of nicotinic acetylcholine receptors in GABAergic neurons of the chick visual system

Double-labeling immunohistochemical methods were used to investigate the occurrence of the alpha8 and alpha5 nicotinic receptor subunits in presumptive GABAergic neurons of the chick nervous system. Nicotinic receptor immunoreactivity was often found in cells exhibiting GABA-like immunoreactivity, especially in the visual system. The alpha8 subunit appeared to be present in presumptive GABAergic cells of the ventral lateral geniculate nucleus, nucleus of the basal optic root of the accessory optic system, and the optic tectum, among several other structures. The alpha5 subunit was also found in GABA-positive neurons, as observed in the lentiform nucleus of the mesencephalon and other pretectal nuclei. The numbers of alpha8- and alpha5-positive neurons that were also GABA-positive represented high percentages of the total number of neurons containing nicotinic receptor labeling in several brain areas, which indicates that most of the alpha8 and alpha5 nicotinic receptor subunits are present in GABAergic cells. Taken together with data from other studies, our results indicate an important role of the nicotinic acetylcholine receptors in the functional organization of GABAergic circuits in the visual system.

Effect of race, genetic population structure, and genetic models in two-locus association studies: clustering of functional renin-angiotensin system gene variants in hypertension association studies

Previous genetic association studies have overlooked the potential for biased results when analyzing different population structures in ethnically diverse populations. The purpose of the present study was to quantify this bias in two-locus association studies conducted on an admixtured urban population. We studied the genetic structure distribution of angiotensin-converting enzyme insertion/deletion (ACE I/D) and angiotensinogen methionine/threonine (M/T) polymorphisms in 382 subjects from three subgroups in a highly admixtured urban population. Group I included 150 white subjects; group II, 142 mulatto subjects, and group III, 90 black subjects. We conducted sample size simulation studies using these data in different genetic models of gene action and interaction and used genetic distance calculation algorithms to help determine the population structure for the studied loci. Our results showed a statistically different population structure distribution of both ACE I/D (P = 0.02, OR = 1.56, 95% CI = 1.05-2.33 for the D allele, white versus black subgroup) and angiotensinogen M/T polymorphism (P = 0.007, OR = 1.71, 95% CI = 1.14-2.58 for the T allele, white versus black subgroup). Different sample sizes are predicted to be determinant of the power to detect a given genotypic association with a particular phenotype when conducting two-locus association studies in admixtured populations. In addition, the postulated genetic model is also a major determinant of the power to detect any association in a given sample size. The present simulation study helped to demonstrate the complex interrelation among ethnicity, power of the association, and the postulated genetic model of action of a particular allele in the context of clustering studies. This information is essential for the correct planning and interpretation of future association studies conducted on this population.