Parametric Sensitivity Analysis of the Most Recent Computational Models of Rabbit Cardiac Pacemaking

The cellular basis of cardiac pacemaking activity, and specifically the quantitative contributions of particular mechanisms, is still debated. Reliable computational models of sinoatrial nodal (SAN) cells may provide mechanistic insights, but competing models are built from different data sets and with different underlying assumptions. To understand quantitative differences between alternative models, we performed thorough parameter sensitivity analyses of the SAN models of Maltsev & Lakatta (2009) and Severi et al (2012). Model parameters were randomized to generate a population of cell models with different properties, simulations performed with each set of random parameters generated 14 quantitative outputs that characterized cellular activity, and regression methods were used to analyze the population behavior. Clear differences between the two models were observed at every step of the analysis. Specifically: (1) SR Ca2+ pump activity had a greater effect on SAN cell cycle length (CL) in the Maltsev model; (2) conversely, parameters describing the funny current (If) had a greater effect on CL in the Severi model; (3) changes in rapid delayed rectifier conductance (GKr) had opposite effects on action potential amplitude in the two models; (4) within the population, a greater percentage of model cells failed to exhibit action potentials in the Maltsev model (27%) compared with the Severi model (7%), implying greater robustness in the latter; (5) confirming this initial impression, bifurcation analyses indicated that smaller relative changes in GKr or Na+-K+ pump activity led to failed action potentials in the Maltsev model. Overall, the results suggest experimental tests that can distinguish between models and alternative hypotheses, and the analysis offers strategies for developing anti-arrhythmic pharmaceuticals by predicting their effect on the pacemaking activity.

The effects of hypocalcemia on spatial alternans and ventricular fibrillation studied with the optical mapping technique

The heart is a wonderful but complex organ: it uses electrochemical mechanisms in order to produce mechanical energy to pump the blood throughout the body and allow the life of humans and animals. This organ can be subject to several diseases and sudden cardiac death (SCD) is the most catastrophic manifestation of these diseases, responsible for the death of a large number of people throughout the world. It is estimated that 325000 Americans annually die for SCD. SCD most commonly occurs as a result of reentrant tachyarrhythmias (ventricular tachycardia (VT) and ventricular fibrillation (VF)) and the identification of those patients at higher risk for the development of SCD has been a difficult clinical challenge. Nowadays, a particular electrocardiogram (ECG) abnormality, “T-wave alternans” (TWA), is considered a precursor of lethal cardiac arrhythmias and sudden death, a sensitive indicator of risk for SCD. TWA is defined as a beat-to-beat alternation in the shape, amplitude, or timing of the T-wave on the ECG, indicative of the underlying repolarization of cardiac cells [5]. In other words TWA is the macroscopic effect of subcellular and celluar mechanisms involving ionic kinetics and the consequent depolarization and repolarization of the myocytes. Experimental activities have shown that TWA on the ECG is a manifestation of an underlying alternation of long and short action potential durations (APDs), the so called APD-alternans, of cardiac myocytes in the myocardium. Understanding the mechanism of APDs-alternans is the first step for preventing them to occur. In order to investigate these mechanisms it’s very important to understand that the biological systems are complex systems and their macroscopic properties arise from the nonlinear interactions among the parts. The whole is greater than the sum of the parts, and it cannot be understood only by studying the single parts. In this sense the heart is a complex nonlinear system and its way of working follows nonlinear dynamics; alternans also, they are a manifestation of a phenomenon typical in nonlinear dynamical systems, called “period-dubling bifurcation”. Over the past decade, it has been demonstrated that electrical alternans in cardiac tissue is an important marker for the development of ventricular fibrillation and a significant predictor for mortality. It has been observed that acute exposure to low concentration of calcium does not decrease the magnitude of alternans and sustained ventricular Fibrillation (VF) is still easily induced under these condition. However with prolonged exposure to low concentration of calcium, alternans disappears, but VF is still inducible. This work is based on this observation and tries to make it clearer. The aim of this thesis is investigate the effect of hypocalcemia spatial alternans and VF doing experiments with canine hearts and perfusing them with a solution with physiological ionic concentration and with a solution with low calcium concentration (hypocalcemia); in order to investigate the so called memory effect, the experimental activity was modified during the way. The experiments were performed with the optical mapping technique, using voltage-sensitive dye, and a custom made Java code was used in post-processing. Finding the Nolasco and Dahlen’s criterion [8] inadequate for the prediction of alternans, and takin into account the experimental results, another criterion, which consider the memory effect, has been implemented. The implementation of this criterion could be the first step in the creation of a method, AP-based, discriminating who is at risk if developing VF. This work is divided into four chapters: the first is a brief presentation of the physiology of the heart; the second is a review of the major theories and discovers in the study of cardiac dynamics; the third chapter presents an overview on the experimental activity and the optical mapping technique; the forth chapter contains the presentation of the results and the conclusions.

Physical characterization and numerical modeling of a WEST-like tokamak divertor plasma

Il presente elaborato è incentrato sulla modellizzazione del plasma di bordo nei dispositivi per la produzione di energia da fusione nucleare noti come tokamak. La tecnologia che nel corso di tutta la seconda metà del XX secolo fino ad oggi è stata sviluppata a questo fine deve necessariamente scontrarsi con alcuni limiti. Nei tokamak il confinamento del plasma è di tipo magnetico e vincola le particelle a muoversi di moto elicoidale all'interno del vessel, tuttavia il confinamento non risulta perfetto e parte dell'energia si scarica sulle pareti della camera, rischiando pertanto di fondere i materiali. Alcune strategie possono essere messe in atto per limitare questo problema, per esempio agendo sulla geometria del tokamak, oppure sulla fisica, inducendo nel plasma una data concentrazione di impurezze che ionizzino irraggiando parte dell'energia di plasma. Proprio tale meccanismo di perdita è stato simulato in un modello monodimensionale di plasma monofluido di bordo. I risultati del codice numerico relativo al modello dimostrano che per concentrazioni di impurezze crescenti è possibile diminuire in modo significativo flusso di calore e temperatura al divertore. Per di più risulta possibile controllare la posizione del fronte di irraggiamento per mezzo di parametri di controllo del plasma quali la pressione. Si osserva inoltre l'insorgere del cosiddetto fenomeno di biforcazione alle basse temperature di divertore, fenomeno in cui il plasma si comporta in modo instabile a causa di fenomeni fisici tipici delle basse energie ("detachment") e a seguito del quale può improvvisamente spegnersi (disruzione). Infine lo stesso modello è stato migliorato inserendo l'ipotesi di plasma bifluido. Anche per gli ioni viene osservato il fenomeno di biforcazione. I risultati numerici evidenziano le dinamiche dello scambio energetico fra le specie gettando le basi di una progettazione efficiente della chimica del plasma finalizzata al raffreddamento del divertore.