Development of Clinical Decision Support Systems based on Mathematical Models of Physiological Systems

In the last years of research, I focused my studies on different physiological problems. Together with my supervisors, I developed/improved different mathematical models in order to create valid tools useful for a better understanding of important clinical issues. The aim of all this work is to develop tools for learning and understanding cardiac and cerebrovascular physiology as well as pathology, generating research questions and developing clinical decision support systems useful for intensive care unit patients. I. ICP-model Designed for Medical Education We developed a comprehensive cerebral blood flow and intracranial pressure model to simulate and study the complex interactions in cerebrovascular dynamics caused by multiple simultaneous alterations, including normal and abnormal functional states of auto-regulation of the brain. Individual published equations (derived from prior animal and human studies) were implemented into a comprehensive simulation program. Included in the normal physiological modelling was: intracranial pressure, cerebral blood flow, blood pressure, and carbon dioxide (CO2) partial pressure. We also added external and pathological perturbations, such as head up position and intracranial haemorrhage. The model performed clinically realistically given inputs of published traumatized patients, and cases encountered by clinicians. The pulsatile nature of the output graphics was easy for clinicians to interpret. The manoeuvres simulated include changes of basic physiological inputs (e.g. blood pressure, central venous pressure, CO2 tension, head up position, and respiratory effects on vascular pressures) as well as pathological inputs (e.g. acute intracranial bleeding, and obstruction of cerebrospinal outflow). Based on the results, we believe the model would be useful to teach complex relationships of brain haemodynamics and study clinical research questions such as the optimal head-up position, the effects of intracranial haemorrhage on cerebral haemodynamics, as well as the best CO2 concentration to reach the optimal compromise between intracranial pressure and perfusion. We believe this model would be useful for both beginners and advanced learners. It could be used by practicing clinicians to model individual patients (entering the effects of needed clinical manipulations, and then running the model to test for optimal combinations of therapeutic manoeuvres). II. A Heterogeneous Cerebrovascular Mathematical Model Cerebrovascular pathologies are extremely complex, due to the multitude of factors acting simultaneously on cerebral haemodynamics. In this work, the mathematical model of cerebral haemodynamics and intracranial pressure dynamics, described in the point I, is extended to account for heterogeneity in cerebral blood flow. The model includes the Circle of Willis, six regional districts independently regulated by autoregulation and CO2 reactivity, distal cortical anastomoses, venous circulation, the cerebrospinal fluid circulation, and the intracranial pressure-volume relationship. Results agree with data in the literature and highlight the existence of a monotonic relationship between transient hyperemic response and the autoregulation gain. During unilateral internal carotid artery stenosis, local blood flow regulation is progressively lost in the ipsilateral territory with the presence of a steal phenomenon, while the anterior communicating artery plays the major role to redistribute the available blood flow. Conversely, distal collateral circulation plays a major role during unilateral occlusion of the middle cerebral artery. In conclusion, the model is able to reproduce several different pathological conditions characterized by heterogeneity in cerebrovascular haemodynamics and can not only explain generalized results in terms of physiological mechanisms involved, but also, by individualizing parameters, may represent a valuable tool to help with difficult clinical decisions. III. Effect of Cushing Response on Systemic Arterial Pressure. During cerebral hypoxic conditions, the sympathetic system causes an increase in arterial pressure (Cushing response), creating a link between the cerebral and the systemic circulation. This work investigates the complex relationships among cerebrovascular dynamics, intracranial pressure, Cushing response, and short-term systemic regulation, during plateau waves, by means of an original mathematical model. The model incorporates the pulsating heart, the pulmonary circulation and the systemic circulation, with an accurate description of the cerebral circulation and the intracranial pressure dynamics (same model as in the first paragraph). Various regulatory mechanisms are included: cerebral autoregulation, local blood flow control by oxygen (O2) and/or CO2 changes, sympathetic and vagal regulation of cardiovascular parameters by several reflex mechanisms (chemoreceptors, lung-stretch receptors, baroreceptors). The Cushing response has been described assuming a dramatic increase in sympathetic activity to vessels during a fall in brain O2 delivery. With this assumption, the model is able to simulate the cardiovascular effects experimentally observed when intracranial pressure is artificially elevated and maintained at constant level (arterial pressure increase and bradicardia). According to the model, these effects arise from the interaction between the Cushing response and the baroreflex response (secondary to arterial pressure increase). Then, patients with severe head injury have been simulated by reducing intracranial compliance and cerebrospinal fluid reabsorption. With these changes, oscillations with plateau waves developed. In these conditions, model results indicate that the Cushing response may have both positive effects, reducing the duration of the plateau phase via an increase in cerebral perfusion pressure, and negative effects, increasing the intracranial pressure plateau level, with a risk of greater compression of the cerebral vessels. This model may be of value to assist clinicians in finding the balance between clinical benefits of the Cushing response and its shortcomings. IV. Comprehensive Cardiopulmonary Simulation Model for the Analysis of Hypercapnic Respiratory Failure We developed a new comprehensive cardiopulmonary model that takes into account the mutual interactions between the cardiovascular and the respiratory systems along with their short-term regulatory mechanisms. The model includes the heart, systemic and pulmonary circulations, lung mechanics, gas exchange and transport equations, and cardio-ventilatory control. Results show good agreement with published patient data in case of normoxic and hyperoxic hypercapnia simulations. In particular, simulations predict a moderate increase in mean systemic arterial pressure and heart rate, with almost no change in cardiac output, paralleled by a relevant increase in minute ventilation, tidal volume and respiratory rate. The model can represent a valid tool for clinical practice and medical research, providing an alternative way to experience-based clinical decisions. In conclusion, models are not only capable of summarizing current knowledge, but also identifying missing knowledge. In the former case they can serve as training aids for teaching the operation of complex systems, especially if the model can be used to demonstrate the outcome of experiments. In the latter case they generate experiments to be performed to gather the missing data.

Ruolo dei glicosaminoglicani su struttura e funzione dei crimps nei tendini e nei legamenti

The biomechanical roles of both tendons and ligaments are fulfilled by extracellular matrix of these tissues. In particular, tension is mainly transmitted and resisted by fibrous proteins (collagen, elastin), whereas compressive load is absorbed by water-soluble glycosaminoglycans (GAGs). GAGs spanning the interfibrillar spaces and interacting with fibrils also seem to play a part in transmitting and resisting tensile stresses. Apart from different functional roles and collagen array, tendons and ligaments share the same basic structure showing periodic undulations of collagen fibers or crimps. Each crimp is composed of many knots of each single fibril or fibrillar crimps. Fibrillar and fiber crimps act as shock absorbers during the initial elongation of both tendons and ligaments and assist the elastic recoil of fibrils and fibers when the tensile stress is removed. The aim of this thesis was to evaluate whether GAGs directly affect the 3D microstructural integrity of fibrillar crimp and fiber crimps in both tendons and ligaments. Achilles tendons and medial collateral ligaments of the knee from eight female Sprague-Dawley rats (90 days old) were digested with chondroitinase ABC to remove GAGs and observed under a scanning electron microscope (SEM). In addition, isolated fibrils from these tissues obtained by mechanical homogenization were analyzed by a transmission electron microscope (TEM). Both samples digested with chondroitinase ABC or mechanically disrupted still showed crimps and fibrillar crimps comparable to tissues with a normal GAGs content. All fibrils in the fibrillar crimp region always twisted leftwards, thus changing their running plane, and then sharply bent, changing their course on a new plane. These data suggest that GAGs do not affect structural integrity or fibrillar crimps functions that seem mainly related to the local fibril leftward twisting and the alternating handedness of collagen from a molecular to a supramolecular level.

Biomechanics of sprint running: a methodological contribution

Sports biomechanics describes human movement from a performance enhancement and an injury reduction perspective. In this respect, the purpose of sports scientists is to support coaches and physicians with reliable information about athletes’ technique. The lack of methods allowing for in-field athlete evaluation as well as for accurate joint force estimates represents, to date, the main limitation to this purpose. The investigations illustrated in the present thesis aimed at providing a contribution towards the development of the above mentioned methods. Two complementary approaches were adopted: a Low Resolution Approach – related to performance assessment – where the use of wearable inertial measurement units is exploited during different phases of sprint running, and a High Resolution Approach – related to joint kinetics estimate for injury prevention – where subject-specific, non-rigid constraints for knee joint kinematic modelling used in multi-body optimization techniques are defined. Results obtained using the Low Resolution Approach indicated that, due to their portability and inexpensiveness, inertial measurement systems are a valid alternative to laboratory-based instrumentation for in-field performance evaluation of sprint running. Using acceleration and angular velocity data, the following quantities were estimated: trunk inclination and angular velocity, instantaneous horizontal velocity and displacement of a point approximating the centre of mass, and stride and support phase durations. As concerns the High Resolution Approach, results indicated that the length of the anterior cruciate and lateral collateral ligaments decreased, while that of the deep bundle of the medial collateral ligament increased significantly during flexion. Variations of the posterior cruciate and the superficial bundle of the medial collateral ligament lengths were concealed by the experimental indeterminacy. A mathematical model was provided that allowed the estimate of subject-specific ligament lengths as a function of knee flexion and that can be integrated in a multi-body optimization procedure.

Cáncer y tratamiento oncológico: Representaciones sociales de la población general, el paciente oncológico y miembros del equipo de salud

Se realizaron tres estudios cualitativos que tuvieron como propósito conocer las representaciones que ha construido la población general, los pacientes oncológicos y los profesionales de la salud, sobre el cáncer, la quimioterapia y el trasplante de médula ósea y realizar un análisis sobre las semejanzas y diferencias entre ellos. Se realizó en la ciudad de Bogotá (Colombia) con 55 personas: 20 pacientes con cáncer en proceso de trasplante de médula ósea, 20 personas no diagnosticadas con cáncer y 15 personas que trabajan en la atención de pacientes con cáncer. Se realizó una entrevista en profundidad con todos los participantes y asociaciones libres, clásicas y por sustitución sobre las palabras “cáncer”, “quimioterapia” y “trasplante de médula”. Los datos conseguidos se analizaron a la luz de la Teoría de las Representaciones Sociales (TRS). El análisis de la información siguió la técnica de análisis cualitativo de contenido para encontrar significados simbólicos y construir, denominar y definir categorías. Para los tres grupos el cáncer es una enfermedad terrible, que puede llevar a la muerte. El personal de salud y la población general creen que la enfermedad genera terror, angustia y miedo. Los pacientes tienen conciencia de la gravedad y del temor consecuente por una enfermedad que lo cambia todo, produce sufrimiento, dolor, obliga a depender de alguien y puede conducir a la muerte. El personal de salud considera que los pacientes lo pueden vivir como un castigo y la población general que puede ser la consecuencia de estilos de vida poco saludables. Para todos, la quimioterapia es un tratamiento para la enfermedad, que por un lado presenta efectos colaterales difíciles y visibles y que producen sentimientos negativos de temor y de angustia y al mismo tiempo constituye una opción y posibilidad de curación. El Trasplante de Médula Ósea representa para todos una oportunidad.