Theoretical investigation of intra- and inter-cellular spatiotemporal calcium patterns in microcirculation

Microcirculatory vessels are lined by endothelial cells (ECs) which are surrounded by a single or multiple layer of smooth muscle cells (SMCs). Spontaneous and agonist induced spatiotemporal calcium (Ca2+) events are generated in ECs and SMCs, and regulated by complex bi-directional signaling between the two layers which ultimately determines the vessel tone. The contractile state of microcirculatory vessels is an important factor in the determination of vascular resistance, blood flow and blood pressure. This dissertation presents theoretical insights into some of the important and currently unresolved phenomena in microvascular tone regulation. Compartmental and continuum models of isolated EC and SMC, coupled EC-SMC and a multi-cellular vessel segment with deterministic and stochastic descriptions of the cellular components were developed, and the intra- and inter-cellular spatiotemporal Ca2+ mobilization was examined.^ Coupled EC-SMC model simulations captured the experimentally observed localized subcellular EC Ca2+ events arising from the opening of EC transient receptor vanilloid 4 (TRPV4) channels and inositol triphosphate receptors (IP3Rs). These localized EC Ca2+ events result in endothelium-derived hyperpolarization (EDH) and Nitric Oxide (NO) production which transmit to the adjacent SMCs to ultimately result in vasodilation. The model examined the effect of heterogeneous distribution of cellular components and channel gating kinetics in determination of the amplitude and spread of the Ca2+ events. The simulations suggested the necessity of co-localization of certain cellular components for modulation of EDH and NO responses. Isolated EC and SMC models captured intracellular Ca2+ wave like activity and predicted the necessity of non-uniform distribution of cellular components for the generation of Ca2+ waves. The simulations also suggested the role of membrane potential dynamics in regulating Ca2+ wave velocity. The multi-cellular vessel segment model examined the underlying mechanisms for the intercellular synchronization of spontaneous oscillatory Ca2+ waves in individual SMC. ^ From local subcellular events to integrated macro-scale behavior at the vessel level, the developed multi-scale models captured basic features of vascular Ca2+ signaling and provide insights for their physiological relevance. The models provide a theoretical framework for assisting investigations on the regulation of vascular tone in health and disease.^

Resum de la tasca realitzada al Center Biomedical Engineering (CBE) del Massachussets Institute of Technology (MIT). Juliol-Agost 2005

Estudi elaborat a partir d’una estada al Center Biomedical Engineering (CBE) del Massachussets Institute of Technology (MIT), durant els mesos de juliol i agost del 2005. S’investiga una metodologia amb l’objectiu d’obtenir biomaterials que puguin actuar de bastida en la interfície os/cartílag, afavorint la diferenciació i creixement cel·lular de cartílag ossificat que pugui actuar d’unió entre l’articulació i l’os. S’experimenta una metodologia per a establir quins són els péptids afavoridors de la formació de teixit ossi utilitzats en materials d’hidroxiapatita. Es conclou que la tecnologia desenvolupada permet disposar d’una plataforma per assajar l’estudi del signaling sobre cèl·lules embrionàries, que permeti desenvolupar materials amb més capacitat diferenciadora.

Hand-in-hand advances in biomedical engineering and sensorimotor restoration.

BACKGROUND: Living in a multisensory world entails the continuous sensory processing of environmental information in order to enact appropriate motor routines. The interaction between our body and our brain is the crucial factor for achieving such sensorimotor integration ability. Several clinical conditions dramatically affect the constant body-brain exchange, but the latest developments in biomedical engineering provide promising solutions for overcoming this communication breakdown. NEW METHOD: The ultimate technological developments succeeded in transforming neuronal electrical activity into computational input for robotic devices, giving birth to the era of the so-called brain-machine interfaces. Combining rehabilitation robotics and experimental neuroscience the rise of brain-machine interfaces into clinical protocols provided the technological solution for bypassing the neural disconnection and restore sensorimotor function. RESULTS: Based on these advances, the recovery of sensorimotor functionality is progressively becoming a concrete reality. However, despite the success of several recent techniques, some open issues still need to be addressed. COMPARISON WITH EXISTING METHOD(S): Typical interventions for sensorimotor deficits include pharmaceutical treatments and manual/robotic assistance in passive movements. These procedures achieve symptoms relief but their applicability to more severe disconnection pathologies is limited (e.g. spinal cord injury or amputation). CONCLUSIONS: Here we review how state-of-the-art solutions in biomedical engineering are continuously increasing expectances in sensorimotor rehabilitation, as well as the current challenges especially with regards to the translation of the signals from brain-machine interfaces into sensory feedback and the incorporation of brain-machine interfaces into daily activities.

Biomedical engineering for healthy ageing. Predictive tools for falls

Falls are common and burdensome accidents among the elderly. About one third of the population aged 65 years or more experience at least one fall each year. Fall risk assessment is believed to be beneficial for fall prevention. This thesis is about prognostic tools for falls for community-dwelling older adults. We provide an overview of the state of the art. We then take different approaches: we propose a theoretical probabilistic model to investigate some properties of prognostic tools for falls; we present a tool whose parameters were derived from data of the literature; we train and test a data-driven prognostic tool. Finally, we present some preliminary results on prediction of falls through features extracted from wearable inertial sensors. Heterogeneity in validation results are expected from theoretical considerations and are observed from empirical data. Differences in studies design hinder comparability and collaborative research. According to the multifactorial etiology of falls, assessment on multiple risk factors is needed in order to achieve good predictive accuracy.

Rapid prototyping for biomedical engineering: current capabilities and Challenges

A new set of manufacturing technologies has emerged in the past decades to address market requirements in a customized way and to provide support for research tasks that require prototypes. These new techniques and technologies are usually referred to as rapid prototyping and manufacturing technologies, and they allow prototypes to be produced in a wide range of materials with remarkable precision in a couple of hours. Although they have been rapidly incorporated into product development methodologies, they are still under development, and their applications in bioengineering are continuously evolving. Rapid prototyping and manufacturing technologies can be of assistance in every stage of the development process of novel biodevices, to address various problems that can arise in the devices' interactions with biological systems and the fact that the design decisions must be tested carefully. This review focuses on the main fields of application for rapid prototyping in biomedical engineering and health sciences, as well as on the most remarkable challenges and research trends.

Fluid dynamics and platelet deposition in a model arterial stenosis

The objective of this study was to gain further understanding and elucidation of the fluid dynamic factors and flow-induced mechanisms of the thrombogenic process of platelet deposition onto, and possible subsequent embolization from, the walls of an arterial stenosis. This has been accomplished by measurement of the axial dependence of platelet deposition within a modeled arterial stenosis for a transitional flow and a completely laminar flow field. The stenotic region of the model was collagen-coated to simulate a damaged endothelial lining of an artery. Fluid dynamics within a stenosis was studied using qualitative flow visualization, and was further compared to the in vitro platelet deposition studies. Normalized platelet density (NPD) measurements indicate decreased levels of NPD in the high shear throat region of the stenosis for a Reynolds number of 300 and a drastic increase in NPD at the throat for a Reynolds number of 175. This study provides further understanding of the flow dynamic effects on thrombus development within a stenosis. ^

Influence of Outlet Boundary Conditions on Cerebrovascular Aneurysm Hemodynamics

Computational fluid dynamic (CFD) studies of blood flow in cerebrovascular aneurysms have potential to improve patient treatment planning by enabling clinicians and engineers to model patient-specific geometries and compute predictors and risks prior to neurovascular intervention. However, the use of patient-specific computational models in clinical settings is unfeasible due to their complexity, computationally intensive and time-consuming nature. An important factor contributing to this challenge is the choice of outlet boundary conditions, which often involves a trade-off between physiological accuracy, patient-specificity, simplicity and speed. In this study, we analyze how resistance and impedance outlet boundary conditions affect blood flow velocities, wall shear stresses and pressure distributions in a patient-specific model of a cerebrovascular aneurysm. We also use geometrical manipulation techniques to obtain a model of the patient’s vasculature prior to aneurysm development, and study how forces and stresses may have been involved in the initiation of aneurysm growth. Our CFD results show that the nature of the prescribed outlet boundary conditions is not as important as the relative distributions of blood flow through each outlet branch. As long as the appropriate parameters are chosen to keep these flow distributions consistent with physiology, resistance boundary conditions, which are simpler, easier to use and more practical than their impedance counterparts, are sufficient to study aneurysm pathophysiology, since they predict very similar wall shear stresses, time-averaged wall shear stresses, time-averaged pressures, and blood flow patterns and velocities. The only situations where the use of impedance boundary conditions should be prioritized is if pressure waveforms are being analyzed, or if local pressure distributions are being evaluated at specific time points, especially at peak systole, where the use of resistance boundary conditions leads to unnaturally large pressure pulses. In addition, we show that in this specific patient, the region of the blood vessel where the neck of the aneurysm developed was subject to abnormally high wall shear stresses, and that regions surrounding blebs on the aneurysmal surface were subject to low, oscillatory wall shear stresses. Computational models using resistance outlet boundary conditions may be suitable to study patient-specific aneurysm progression in a clinical setting, although several other challenges must be addressed before these tools can be applied clinically.

Effects of Load History on Ovine Spinal Range of Motion

Loading of spinal motion segment units alters biomechanical properties by modifying flexibility and range of motion. This study utilizes angular displacement due to an applied bending moment to assess biomechanical function during high-magnitude and prolonged compressive loading of ovine lumbar motion segments. High compressive loads, representative of physiological lifestyle and occupational behaviors, appear to limit fluid recovery of the intervertebral disc, thereby modifying spinal flexibility and increasing spinal instability. Intermittent extensions, or backwards bending movements, may provide a protective effect against the load-induced spinal instability. This study contributes a greater understanding of the effects of load history on the function and health of the lumbar spine. Findings may inform future efforts investigating adjustments in spinal posture to preserve or promote the recovery of lumbar spinal biomechanics.

Movement effects on the flow physics and nutrient delivery in engineered valvular tissues

Mechanical conditioning has been shown to promote tissue formation in a wide variety of tissue engineering efforts. However the underlying mechanisms by which external mechanical stimuli regulate cells and tissues are not known. This is particularly relevant in the area of heart valve tissue engineering (HVTE) owing to the intense hemodynamic environments that surround native valves. Some studies suggest that oscillatory shear stress (OSS) caused by steady flow and scaffold flexure play a critical role in engineered tissue formation derived from bone marrow derived stem cells (BMSCs). In addition, scaffold flexure may enhance nutrient (e.g. oxygen, glucose) transport. In this study, we computationally quantified the i) magnitude of fluid-induced shear stresses; ii) the extent of temporal fluid oscillations in the flow field using the oscillatory shear index (OSI) parameter, and iii) glucose and oxygen mass transport profiles. Noting that sample cyclic flexure induces a high degree of oscillatory shear stress (OSS), we incorporated moving boundary computational fluid dynamic simulations of samples housed within a bioreactor to consider the effects of: 1) no flow, no flexure (control group), 2) steady flow-alone, 3) cyclic flexure-alone and 4) combined steady flow and cyclic flexure environments. We also coupled a diffusion and convention mass transport equation to the simulated system. We found that the coexistence of both OSS and appreciable shear stress magnitudes, described by the newly introduced parameter OSI-:τ: explained the high levels of engineered collagen previously observed from combining cyclic flexure and steady flow states. On the other hand, each of these metrics on its own showed no association. This finding suggests that cyclic flexure and steady flow synergistically promote engineered heart valve tissue production via OSS, so long as the oscillations are accompanied by a critical magnitude of shear stress. In addition, our simulations showed that mass transport of glucose and oxygen is enhanced by sample movement at low sample porosities, but did not play a role in highly porous scaffolds. Preliminary in-house in vitro experiments showed that cell proliferation and phenotype is enhanced in OSI-:τ: environments.^

Translation of biomechanical concepts in bone tissue engineering: from animal study to revision knee arthroplasty.

Bone defects in revision knee arthroplasty are often located in load-bearing regions. The goal of this study was to determine whether a physiologic load could be used as an in situ osteogenic signal to the scaffolds filling the bone defects. In order to answer this question, we proposed a novel translation procedure having four steps: (1) determining the mechanical stimulus using finite element method, (2) designing an animal study to measure bone formation spatially and temporally using micro-CT imaging in the scaffold subjected to the estimated mechanical stimulus, (3) identifying bone formation parameters for the loaded and non-loaded cases appearing in a recently developed mathematical model for bone formation in the scaffold and (4) estimating the stiffness and the bone formation in the bone-scaffold construct. With this procedure, we estimated that after 3 years mechanical stimulation increases the bone volume fraction and the stiffness of scaffold by 1.5- and 2.7-fold, respectively, compared to a non-loaded situation.

Orthotropic active strain models for the numerical simulation of cardiac biomechanics

An active strain formulation for orthotropic constitutive laws arising in cardiac mechanics modeling is introduced and studied. The passive mechanical properties of the tissue are described by the Holzapfel-Ogden relation. In the active strain formulation, the Euler-Lagrange equations for minimizing the total energy are written in terms of active and passive deformation factors, where the active part is assumed to depend, at the cell level, on the electrodynamics and on the specific orientation of the cardiac cells. The well-posedness of the linear system derived from a generic Newton iteration of the original problem is analyzed and different mechanical activation functions are considered. In addition, the active strain formulation is compared with the classical active stress formulation from both numerical and modeling perspectives. Taylor-Hood and MINI finite elements are employed to discretize the mechanical problem. The results of several numerical experiments show that the proposed formulation is mathematically consistent and is able to represent the main key features of the phenomenon, while allowing savings in computational costs.

Development of a prototype for measuring plantar forces

Studies on the distribution of plantar pressure between the sole of the foot and the ground were developed before the 19th century. Currently, the most often employed plantar pressure measurement systems are Pedar® and FScan®, which have restrictions such as operational difficulty and high cost. In the present study, a device was constructed from two pressure plates capable of measuring plantar forces in discreet areas of the feet at a low cost, using strain-gages attached to sixteen strategic points of the mechanical elements. Sixteen prismatic beams were soldered to each frame, for which the free extremity of each beam represented a specific point of the foot. Two strain gauges were attached to each beam - one near the upper fixed extremity and the other near the lower fixed extremity. Using a Wheatstone bridge electric circuit, the gauges were used to measure the force acting on the extremity of the beam. Precision and accuracy of the prototype was about 10%. In some measurements, accuracy was 2%. The low precision and accuracy were mainly due to the restrictions of the available equipment, which only permitted four measurements at a time. Thus, it was necessary for participants to stand on the plates four separate times, which signified possible changes in the position of the feet on the pressure plates. Despite some limitations, the aim was achieved. The prototype has been used in some studies and represents a contribution to biomechanics, demonstrating the viability of using strain gauges.

Mechanical behavior of isolated intervertebral disc

The study on several components of intervertebral joints is essential to understand the spine's degenerative mechanisms and to assess the best method for their treatment. For such study it is necessary to know the mechanical properties of the isolated intervertebral disc (ID) mechanical properties and, it is necessary to evaluate its stresses and strains. In order to assess the ID displacements, a fine, U-shaped blade was developed, over which two extensometers connected in a Wheatstone bridge were placed. The device was then tested on porcine spine ID, where compression loads were applied and the extremities displacements of the blade coupled to the intervertebral disc were measured. Stress/strain diagram, both on the compression and on the decompression phases, evidencing the non-linear nature of such relationship. With the experiment, it was possible to obtain approximate values of the longitudinal elasticity module (E) of the disc material and of the Poisson coefficient (n ). After several tests, E results are compatible with those obtained by others studies, with very simple and low-cost device. This experiments can be used for obtained others mechanical properties of isolated ID with precision and accuracy.