Mathematical models for cellular aggregation: the chemotactic instability and clustering formation

In this thesis we present a mathematical formulation of the interaction between microorganisms such as bacteria or amoebae and chemicals, often produced by the organisms themselves. This interaction is called chemotaxis and leads to cellular aggregation. We derive some models to describe chemotaxis. The first is the pioneristic Keller-Segel parabolic-parabolic model and it is derived by two different frameworks: a macroscopic perspective and a microscopic perspective, in which we start with a stochastic differential equation and we perform a mean-field approximation. This parabolic model may be generalized by the introduction of a degenerate diffusion parameter, which depends on the density itself via a power law. Then we derive a model for chemotaxis based on Cattaneo's law of heat propagation with finite speed, which is a hyperbolic model. The last model proposed here is a hydrodynamic model, which takes into account the inertia of the system by a friction force. In the limit of strong friction, the model reduces to the parabolic model, whereas in the limit of weak friction, we recover a hyperbolic model. Finally, we analyze the instability condition, which is the condition that leads to aggregation, and we describe the different kinds of aggregates we may obtain: the parabolic models lead to clusters or peaks whereas the hyperbolic models lead to the formation of network patterns or filaments. Moreover, we discuss the analogy between bacterial colonies and self gravitating systems by comparing the chemotactic collapse and the gravitational collapse (Jeans instability).

Electrolyte disorders in the human ventricle. A simulation study.

Ventricular cells are immersed in a bath of electrolytes and these ions are essential for a healthy heart and a regular rhythm. Maintaining physiological concentration of them is fundamental for reducing arrhythmias and risk of sudden cardiac death, especially in haemodialysis patients and in the heart diseases treatments. Models of electrically activity of the heart based on mathematical formulation are a part of the efforts to improve the understanding and prediction of heart behaviour. Modern models incorporate the extensive and ever increasing amounts of experimental data in incorporating biophysically detailed mechanisms to allow the detailed study of molecular and subcellular mechanisms of heart disease. The goal of this project was to simulate the effects of changes in potassium and calcium concentrations in the extracellular space between experimental data and and a description incorpored into two modern biophysically detailed models (Grandi et al. Model; O’Hara Rudy Model). Moreover the task was to analyze the changes in the ventricular electrical activity, in particular by studying the modifications on the simulated electrocardiographic signal. We used the cellular information obtained by the heart models in order to build a 1D tissue description. The fibre is composed by 165 cells, it is divided in four groups to differentiate the cell types that compound human ventricular tissue. The main results are the following: Grandi et al. (GBP) model is not even able to reproduce the correct action potential profile in hyperkalemia. Data from hospitalized patients indicates that the action potential duration (APD) should be shorter than physiological state but in this model we have the opposite. From the potassium point of view the results obtained by using O’Hara model (ORD) are in agreement with experimental data for the single cell action potential in hypokalemia and hyperkalemia, most of the currents follow the data from literature. In the 1D simulations we were able to reproduce ECGs signal in most the potassium concentrations we selected for this study and we collected data that can help physician in understanding what happens in ventricular cells during electrolyte disorder. However the model fails in the conduction of the stimulus under hyperkalemic conditions. The model emphasized the ECG modifications when the K+ is slightly more than physiological value. In the calcium setting using the ORD model we found an APD shortening in hypocalcaemia and an APD lengthening in hypercalcaemia, i.e. the opposite to experimental observation. This wrong behaviour is kept in one dimensional simulations bringing a longer QT interval in the ECG under higher [Ca2+]o conditions and vice versa. In conclusion it has highlighted that the actual ventricular models present in literature, even if they are useful in the original form, they need an improvement in the sensitivity of these two important electrolytes. We suggest an use of the GBP model with modifications introduced by Carro et al. who understood that the failure of this model is related to the Shannon et al. model (a rabbit model) from which the GBP model was built. The ORD model should be modified in the Ca2+ - dependent IcaL and in the influence of the Iks in the action potential for letting it him produce a correct action potential under different calcium concentrations. In the 1D tissue maybe a heterogeneity setting of intra and extracellular conductances for the different cell types should improve a reproduction of the ECG signal.

Parallelization of the algorithm WHAM with NVIDIA CUDA.

The aim of my thesis is to parallelize the Weighting Histogram Analysis Method (WHAM), which is a popular algorithm used to calculate the Free Energy of a molucular system in Molecular Dynamics simulations. WHAM works in post processing in cooperation with another algorithm called Umbrella Sampling. Umbrella Sampling has the purpose to add a biasing in the potential energy of the system in order to force the system to sample a specific region in the configurational space. Several N independent simulations are performed in order to sample all the region of interest. Subsequently, the WHAM algorithm is used to estimate the original system energy starting from the N atomic trajectories. The parallelization of WHAM has been performed through CUDA, a language that allows to work in GPUs of NVIDIA graphic cards, which have a parallel achitecture. The parallel implementation may sensibly speed up the WHAM execution compared to previous serial CPU imlementations. However, the WHAM CPU code presents some temporal criticalities to very high numbers of interactions. The algorithm has been written in C++ and executed in UNIX systems provided with NVIDIA graphic cards. The results were satisfying obtaining an increase of performances when the model was executed on graphics cards with compute capability greater. Nonetheless, the GPUs used to test the algorithm is quite old and not designated for scientific calculations. It is likely that a further performance increase will be obtained if the algorithm would be executed in clusters of GPU at high level of computational efficiency. The thesis is organized in the following way: I will first describe the mathematical formulation of Umbrella Sampling and WHAM algorithm with their apllications in the study of ionic channels and in Molecular Docking (Chapter 1); then, I will present the CUDA architectures used to implement the model (Chapter 2); and finally, the results obtained on model systems will be presented (Chapter 3).

Algoritmi esatti per il Job Shop Scheduling: approcci Mathematical Programming

In questa tesi ci occuperemo di fornire un modello MIP di base e di alcune sue varianti, realizzate allo scopo di comprenderne il comportamento ed eventualmente migliorarne l’efficienza. Le diverse varianti sono state costruite agendo in particolar modo sulla definizione di alcuni vincoli, oppure sui bound delle variabili, oppure ancora nell’obbligare il risolutore a focalizzarsi su determinate decisioni o specifiche variabili. Sono stati testati alcuni dei problemi tipici presenti in letteratura e i diversi risultati sono stati opportunamente valutati e confrontati. Tra i riferimenti per tale confronto sono stati considerati anche i risultati ottenibili tramite un modello Constraint Programming, che notoriamente produce risultati apprezzabili in ambito di schedulazione. Un ulteriore scopo della tesi è, infatti, comparare i due approcci Mathematical Programming e Constraint Programming, identificandone quindi i pregi e gli svantaggi e provandone la trasferibilità al modello raffrontato.

Mathematical Modeling to Investigate Antiarrhythmic Drug Side Effects: Rate-Dependence Role in Ionic Currents and Action Potentials Shape in the O’Hara Model

Sudden cardiac death due to ventricular arrhythmia is one of the leading causes of mortality in the world. In the last decades, it has proven that anti-arrhythmic drugs, which prolong the refractory period by means of prolongation of the cardiac action potential duration (APD), play a good role in preventing of relevant human arrhythmias. However, it has long been observed that the “class III antiarrhythmic effect” diminish at faster heart rates and that this phenomenon represent a big weakness, since it is the precise situation when arrhythmias are most prone to occur. It is well known that mathematical modeling is a useful tool for investigating cardiac cell behavior. In the last 60 years, a multitude of cardiac models has been created; from the pioneering work of Hodgkin and Huxley (1952), who first described the ionic currents of the squid giant axon quantitatively, mathematical modeling has made great strides. The O’Hara model, that I employed in this research work, is one of the modern computational models of ventricular myocyte, a new generation began in 1991 with ventricular cell model by Noble et al. Successful of these models is that you can generate novel predictions, suggest experiments and provide a quantitative understanding of underlying mechanism. Obviously, the drawback is that they remain simple models, they don’t represent the real system. The overall goal of this research is to give an additional tool, through mathematical modeling, to understand the behavior of the main ionic currents involved during the action potential (AP), especially underlining the differences between slower and faster heart rates. In particular to evaluate the rate-dependence role on the action potential duration, to implement a new method for interpreting ionic currents behavior after a perturbation effect and to verify the validity of the work proposed by Antonio Zaza using an injected current as a perturbing effect.

A mathematical introduction to Kerr black holes

Lo scopo della tesi è descrivere i buchi neri di Kerr. Dopo aver introdotto tutti gli strumenti matematici necessari quali tensori, vettori di Killing e geodetiche, enunceremo la metrica di Kerr, il teorema no-hair e il frame-dragging. In seguito, a partire dalla metrica di Kerr, calcoleremo e descriveremo le ergosfere, gli orizzonti degli eventi e il moto dei fotoni nel piano equatoriale.

Synthesis, characterization and formulation of a cathode active material: Copper nitroprusside

Nowadays, rechargeable Li-ion batteries play an important role in portable consumer devices. Formulation of such batteries is improvable by researching new cathodic materials that present higher performances of cyclability and negligible efficiency loss over cycles. Goal of this work was to investigate a new cathodic material, copper nitroprusside, which presents a porous 3D framework. Synthesis was carried out by a low-cost and scalable co-precipitation method. Subsequently, the product was characterized by means of different techniques, such as TGA, XRF, CHN elemental analysis, XRD, Mössbauer spectroscopy and cyclic voltammetry. Electrochemical tests were finally performed both in coin cells and by using in situ cells: on one hand, coin cells allowed different formulations to be easily tested, on the other operando cycling led a deeper insight to insertion process and both chemical and physical changes. Results of several tests highlighted a non-reversible electrochemical behavior of the material and a rapid capacity fading over time. Moreover, operando techniques report that amorphisation occurs during the discharge.

Keller-Segel system, chemotaxis and applications to a mathematical model for Alzheimer's disease

In questa tesi viene presentato il modello di Keller-Segel per la chemiotassi, un sistema di tipo parabolico-ellittico che appare nella descrizione di molti fenomeni in ambito biologico e medico. Viene mostrata l'esistenza globale della soluzione debole del modello, per dati iniziali sufficientemente piccoli in dimensione N>2. La scelta di dati iniziali abbastanza grandi invece può causare il blow-up della soluzione e viene mostrato sotto quali condizioni questo si verifica. Infine il modello della chemiotassi è stato applicato per descrivere una fase della malattia di Alzheimer ed è stata effettuata un'analisi di stabilità del sistema.