Computational Modelling of Cardiac Electrophysiology: from Cell to Bedside

Heart diseases are the leading cause of death worldwide, both for men and women. However, the ionic mechanisms underlying many cardiac arrhythmias and genetic disorders are not completely understood, thus leading to a limited efficacy of the current available therapies and leaving many open questions for cardiac electrophysiologists. On the other hand, experimental data availability is still a great issue in this field: most of the experiments are performed in vitro and/or using animal models (e.g. rabbit, dog and mouse), even when the final aim is to better understand the electrical behaviour of in vivo human heart either in physiological or pathological conditions. Computational modelling constitutes a primary tool in cardiac electrophysiology: in silico simulations, based on the available experimental data, may help to understand the electrical properties of the heart and the ionic mechanisms underlying a specific phenomenon. Once validated, mathematical models can be used for making predictions and testing hypotheses, thus suggesting potential therapeutic targets. This PhD thesis aims to apply computational cardiac modelling of human single cell action potential (AP) to three clinical scenarios, in order to gain new insights into the ionic mechanisms involved in the electrophysiological changes observed in vitro and/or in vivo. The first context is blood electrolyte variations, which may occur in patients due to different pathologies and/or therapies. In particular, we focused on extracellular Ca2+ and its effect on the AP duration (APD). The second context is haemodialysis (HD) therapy: in addition to blood electrolyte variations, patients undergo a lot of other different changes during HD, e.g. heart rate, cell volume, pH, and sympatho-vagal balance. The third context is human hypertrophic cardiomyopathy (HCM), a genetic disorder characterised by an increased arrhythmic risk, and still lacking a specific pharmacological treatment.

Modelling, simulation and characterization of epithelial cell culture biochip

A microfluidic Organ-on-Chip has been developed for monitoring the epithelial cells monolayer. Equivalent circuit Model was used to determine the electrical properties from the impedance spectra of the epithelial cells monolayer. Black platinum on platinum electrodes was electrochemically deposited onto the surface of electrodes to reduce the influence of the electrical double layer on the impedance measurements. Measurements of impedance with an Impedance Analyzer were done to validate the equivalent circuit model and the decrease of the double layer effect. A Lock-in Amplifier was designed to measure the impedance.

Tomographic 3D reconstruction through linear inversion for breast microwave imaging

In the present thesis we address the problem of detecting and localizing a small spherical target with characteristic electrical properties inside a volume of cylindrical shape, representing female breast, with MWI. One of the main works of this project is to properly extend the existing linear inversion algorithm from planar slice to volume reconstruction; results obtained, under the same conditions and experimental setup are reported for the two different approaches. Preliminar comparison and performance analysis of the reconstruction algorithms is performed via numerical simulations in a software-created environment: a single dipole antenna is used for illuminating the virtual breast phantom from different positions and, for each position, the corresponding scattered field value is registered. Collected data are then exploited in order to reconstruct the investigation domain, along with the scatterer position, in the form of image called pseudospectrum. During this process the tumor is modeled as a dielectric sphere of small radius and, for electromagnetic scattering purposes, it's treated as a point-like source. To improve the performance of reconstruction technique, we repeat the acquisition for a number of frequencies in a given range: the different pseudospectra, reconstructed from single frequency data, are incoherently combined with MUltiple SIgnal Classification (MUSIC) method which returns an overall enhanced image. We exploit multi-frequency approach to test the performance of 3D linear inversion reconstruction algorithm while varying the source position inside the phantom and the height of antenna plane. Analysis results and reconstructed images are then reported. Finally, we perform 3D reconstruction from experimental data gathered with the acquisition system in the microwave laboratory at DIFA, University of Bologna for a recently developed breast-phantom prototype; obtained pseudospectrum and performance analysis for the real model are reported.

Regulation of Cx45 hemichannels mediated by extracellular and intracellular calcium

Connexin45 (Cx45) hemichannels (HCs) open in the absence of Ca(2+) and close in its presence. To elucidate the underlying mechanisms, we examined the role of extra- and intracellular Ca(2+) on the electrical properties of HCs. Experiments were performed on HeLa cells expressing Cx45 using electrical (voltage clamp) and optical (Ca(2+) imaging) methods. HCs exhibit a time- and voltage-dependent current (I(hc)), activating with depolarization and inactivating with hyperpolarization. Elevation of [Ca(2+)](o) from 20 nM to 2 μM reversibly decreases I(hc), decelerates its rate of activation, and accelerates its deactivation. Our data suggest that [Ca(2+)](o) modifies the channel properties by adhering to anionic sites in the channel lumen and/or its outer vestibule. In this way, it blocks the channel pore and reversibly lowers I(hc) and modifies its kinetics. Rapid lowering of [Ca(2+)](o) from 2 mM to 20 nM, achieved early during a depolarizing pulse, led to an outward I(hc) that developed with virtually no delay and grew exponentially in time paralleled by unaffected [Ca(2+)](i). A step increase of [Ca(2+)](i) evoked by photorelease of Ca(2+) early during a depolarizing pulse led to a transient decrease of I(hc) superimposed on a growing outward I(hc); a step decrease of [Ca(2+)](i) elicited by photoactivation of a Ca(2+) scavenger provoked a transient increase in I(hc). Hence, it is tempting to assume that Ca(2+) exerts a direct effect on Cx45 hemichannels.

Subconductance states of Cx30 gap junction channels: data from transfected HeLa cells versus data from a mathematical model

Human HeLa cells expressing mouse connexin30 were used to study the electrical properties of gap junction channel substates. Experiments were performed on cell pairs using a dual voltage-clamp method. Single-channel currents revealed discrete levels attributable to a main state, a residual state, and five substates interposed, suggesting the operation of six subgates provided by the six connexins of a gap junction hemichannel. Substate conductances, gamma(j,substate), were unevenly distributed between the main-state and the residual-state conductance (gamma(j,main state) = 141 pS, gamma(j,residual state) = 21 pS). Activation of the first subgate reduced the channel conductance by approximately 30%, and activation of subsequent subgates resulted in conductance decrements of 10-15% each. Current transitions between the states were fast (<2 ms). Substate events were usually demarcated by transitions from and back to the main state; transitions among substates were rare. Hence, subgates are recruited simultaneously rather than sequentially. The incidence of substate events was larger at larger gradients of V(j). Frequency and duration of substate events increased with increasing number of synchronously activated subgates. Our mathematical model, which describes the operation of gap junction channels, was expanded to include channel substates. Based on the established V(j)-sensitivity of gamma(j,main state) and gamma(j,residual state), the simulation yielded unique functions gamma(j,substate) = f(V(j)) for each substate. Hence, the spacing of subconductance levels between the channel main state and residual state were uneven and characteristic for each V(j).

Calcium indicator loading of neurons using single-cell electroporation

Studies of subcellular Ca(2+) signaling rely on methods for labeling cells with fluorescent Ca(2+) indicator dyes. In this study, we demonstrate the use of single-cell electroporation for Ca(2+) indicator loading of individual neurons and small neuronal networks in rat neocortex in vitro and in vivo. Brief voltage pulses were delivered through glass pipettes positioned close to target cells. This approach resulted in reliable and rapid (within seconds) loading of somata and subsequent complete labeling of dendritic and axonal arborizations. By using simultaneous whole-cell recordings in brain slices, we directly addressed the effect of electroporation on neurons. Cell viability was high (about 85%) with recovery from the membrane permeabilization occurring within a minute. Electrical properties of recovered cells were indistinguishable before and after electroporation. In addition, Ca(2+) transients with normal appearance could be evoked in dendrites, spines, and axonal boutons of electroporated cells. Using negative-stains of somata, targeted single-cell electroporation was equally applicable in vivo. We conclude that electroporation is a simple approach that permits Ca(2+) indicator loading of multiple cells with low background staining within a short amount of time, which makes it especially well suited for functional imaging of subcellular Ca(2+) dynamics in small neuronal networks.

Immobilizing bacteriorhodopsin on a single electron transistor

As awareness of potential human and environmental impacts from toxins has increased, so has the development of innovative sensors. Bacteriorhodopsin (bR) is a light activated proton pump contained in the purple membrane (PM) of the bacteria Halobacterium salinarum. Bacteriorhodopsin is a robust protein which can function in both wet and dry states and can withstand extreme environmental conditions. A single electron transistor(SET) is a nano-scale device that exploits the quantum mechanical properties of electrons to switch on and off. SETs have tremendous potential in practical applications due to their size, ultra low power requirements, and electrometer-like sensitivity. The main goal of this research was to create a bionanohybrid device by integrating bR with a SET device. This was achieved by a multidisciplinary approach. The SET devices were created by a combination of sputtering, photolithography, and focused ion beam machining. The bionanomaterial bacteriorhodopsin was created through oxidative fermentation and a series of transmembrane purification processes. The bR was then integrated with the SET by electrophoretic deposition, creating a bionanohybrid device. The bionanohybrid device was then characterized using a semiconductor parametric analyzer. Characterization demonstrated that the bR modulated the operational characteristics of the SET when bR was activated with light within its absorbance spectrum. To effectively integrate bacteriorhodopsin with microelectromechanical systems (MEMS) and nanoelectromechanical systems (NEMS), it is critical to know the electrical properties of the material and to understand how it will affect the functionality of the device. Tests were performed on dried films of bR to determine if there is a relationship between inductance, capacitance, and resistance (LCR) measurements and orientation, light-on/off, frequency, and time. The results indicated that the LCR measurements of the bR depended on the thickness and area of the film, but not on the orientation, as with other biological materials such as muscle. However, there was a transient LCR response for both oriented and unoriented bR which depended on light intensity. From the impedance measurements an empirical model was suggested for the bionanohybrid device. The empirical model is based on the dominant electrical characteristics of the bR which were the parallel capacitance and resistance. The empirical model suggests that it is possible to integrate bR with a SET without influencing its functional characteristics.

NANOSCALE ELECTROCHEMISTRY BY IN-SITU TRANSMISSION ELECTRON MICROSCOPY

Nanoscale research in energy storage has recently focused on investigating the properties of nanostructures in order to increase energy density, power rate, and capacity. To better understand the intrinsic properties of nanomaterials, a new and advanced in situ system was designed that allows atomic scale observation of materials under external fields. A special holder equipped with a scanning tunneling microscopy (STM) probe inside a transmission electron microscopy (TEM) system was used to perform the in situ studies on mechanical, electrical, and electrochemical properties of nanomaterials. The nanostructures of titanium dioxide (TiO2) nanotubes are characterized by electron imaging, diffraction, and chemical analysis techniques inside TEM. TiO2 nanotube is one of the candidates as anode materials for lithium ion batteries. It is necessary to study their morphological, mechanical, electrical, and electrochemical properties at atomic level. The synthesis of TiO2 nanotubes showed that the aspect ratio of TiO2 could be controlled by processing parameters, such as anodization time and voltage. Ammonium hydroxide (NH4OH) treated TiO2 nanotubes showed unexpected instability. Observation revealed the nanotubes were disintegrated into nanoparticles and the tubular morphology was vanished after annealing. The nitrogen compounds incorporated in surface defects weaken the nanotube and result in the collapse of nanotube into nanoparticles during phase transformation. Next, the electrical and mechanical properties of TiO2 nanotubes were studied by in situ TEM system. Phase transformation of anatase TiO2 nanotubes into rutile nanoparticles was studied by in situ Joule heating. The results showed that single anatase TiO2 nanotubes broke into ultrafine small anatase nanoparticles. On further increasing the bias, the nanoclusters of anatase particles became prone to a solid state reaction and were grown into stable large rutile nanoparticles. The relationship between mechanical and electrical properties of TiO2 nanotubes was also investigated. Initially, both anatase and amorphous TiO2 nanotubes were characterized by using I-V test to demonstrate the semiconductor properties. The observation of mechanical bending on TiO2 nanotubes revealed that the conductivity would increase when bending deformation happened. The defects on the nanotubes created by deformation helped electron transportation to increase the conductivity. Lastly, the electrochemical properties of amorphous TiO2 nanotubes were characterized by in situ TEM system. The direct chemical and imaging evidence of lithium-induced atomic ordering in amorphous TiO2 nanotubes was studied. The results indicated that the lithiation started with the valance reduction of Ti4+ to Ti3+ leading to a LixTiO2 intercalation compound. The continued intercalation of Li ions in TiO2 nanotubes triggered an amorphous to crystalline phase transformation. The crystals were formed as nano islands and identified to be Li2Ti2O4 with cubic structure (a = 8.375 Å). This phase transformation is associated with local inhomogeneities in Li distribution. Based on these observations, a new reaction mechanism is proposed to explain the first cycle lithiation behavior in amorphous TiO2 nanotubes.

Modeling Calcium Concentration and Biochemical Reactions

Almost all regions of the brain receive one or more neuromodulatory inputs, and disrupting these inputs produces deficits in neuronal function. Neuromodulators act through intracellular second messenger pathways to influence the electrical properties of neurons, integration of synaptic inputs, spatio-temporal firing dynamics of neuronal networks, and, ultimately, systems behavior. Second messengers pathways consist of series of bimolecular reactions, enzymatic reactions, and diffusion. Calcium is the second messenger molecule with the most effectors, and thus is highly regulated by buffers, pumps and intracellular stores. Computational modeling provides an innovative, yet practical method to evaluate the spatial extent, time course and interaction among second messenger pathways, and the interaction of second messengers with neuron electrical properties. These processes occur both in compartments where the number of molecules are large enough to describe reactions deterministically (e.g. cell body), and in compartments where the number of molecules is small enough that reactions occur stochastically (e.g. spines). – In this tutorial, I explain how to develop models of second messenger pathways and calcium dynamics. The first part of the tutorial explains the equations used to model bimolecular reactions, enzyme reactions, calcium release channels, calcium pumps and diffusion. The second part explains some of the GENESIS, Kinetikit and Chemesis objects that implement the appropriate equations. In depth explanation of calcium and second messenger models is provided by reviewing code, both in XPP, Chemesis and Kinetikit, that implements simple models of calcium dynamics and second messenger cascades.

Functional Characterization and Comparison of Intercellular Communication in Stem Cell-Derived Cardiomyocytes

One novel treatment strategy for the diseased heart focuses on the use of pluripotent stem cell-derived cardiomyocytes (SC-CMs) to overcome the heart's innate deficiency for self-repair. However, targeted application of SC-CMs requires in-depth characterization of their true cardiogenic potential in terms of excitability and intercellular coupling at cellular level and in multicellular preparations. In this study, we elucidated the electrical characteristics of single SC-CMs and intercellular coupling quality of cell pairs, and concomitantly compared them with well-characterized murine native neonatal and immortalized HL-1 cardiomyocytes. Firstly, we investigated the electrical properties and Ca2+ signaling mechanisms specific to cardiac contraction in single SC-CMs. Despite heterogeneity of the new cardiac cell population, their electrophysiological activity and Ca2+ handling were similar to native cells. Secondly, we investigated the capability of paired SC-CMs to form an adequate subunit of a functional syncytium and analyzed gap junctions and signal transmission by dye transfer in cell pairs. We discovered significantly diminished coupling in SC-CMs compared with native cells, which could not be enhanced by a coculture approach combining SC-CMs and primary CMs. Moreover, quantitative and structural analysis of gap junctions presented significantly reduced connexin expression levels compared with native CMs. Strong dependence of intercellular coupling on gap junction density was further confirmed by computational simulations. These novel findings demonstrate that despite the cardiogenic electrophysiological profile, SC-CMs present significant limitations in intercellular communication. Inadequate coupling may severely impair functional integration and signal transmission, which needs to be carefully considered for the prospective use of SC-CMs in cardiac repair.

CNT and PDCs: A fruitful association? Study of a polycarbosilane–MWCNT composite

The effect of MWCNT introduction in a polycarbosilane based ceramic on its electrical properties is presented. The electrical conductivity of two MWCNT powders was measured under dynamic compaction up to 20 MPa when it reached 3–5 S/cm. The compaction behavior was also analyzed and modeled. A composite was then realized using allylhydridopolycarbosilane SMP10® and divinylbenzene as matrix. Intact 10 mm MWCNT-SiC ceramic discs samples with 2 wt.% filler load were produced pressure-less via liquid route despite the linear shrinkage of about 30%. Nanotubes microstructure and distribution in the matrix were confirmed after pyrolysis with TEM and SEM analysis. Anyhow similar electrical conductivity values after pyrolysis between the loaded and unloaded samples were measured. The microstructure analysis via XRD and TEM revealed that the percolative carbon network formed through the use of divinylbenzene improves the electric conductivity more than that of MWCNT addition and also simplifies the whole process.

Characterization of HFO2 Capacitors

Hafnium oxide (HfOn) is a promising dielectric for future microelectronic applications. Hf02 thin films (10-75nm) were deposited on Pt/Si02/Si substrates by pulsed DC magnetron reactive sputtering. Top electrodes of Pt were formed by e-beam evapo- ration through an aperture mask on the samples to create MIM (Metal-Insulator-Metal) capacitors. Various processing conditions (Arloz ratio, DC power and deposition rate) and post-deposition annealing conditions (time and temperature) were investigated. The structure of the Hf02 films was characterized by X-ray diffraction (XRD) and the roughness was measured by a profilometer. The electrical properties were characterized in terms of their relative permittivity (E,(T) and ~,.(f)) and leakage behavior (I-V, I-T and I- time). The electrical measurements were performed over a temperature range from -5 to 200°C. For the samples with best experimental results, the relative permittivity of HfOa was found to be -- 27 after anneal and increased by 0.027%/"C with increasing temperature over the measured temperature range. At 25"C, leakage current density was below lop8 ~ l c m ' at 1 volt. The leakage current increased with temperature above a specific threshold temperature below which the leakage current didn't change much. The leakage current increased with voltage. At voltages below lvolt, it's ohmic; at higher voltages, it follows Schottky model. The breakdown field is - 1 . 8 2 ~ lo6 Vlcm. The optical bandgap was measured with samples deposited on quartz substrates to be 5.4eV after anneal.

Gridded bathymetry from multibeam echosounder EM122 data of the cruise MSM35

The aim of Maria S. MERIAN voyage MSM35 was to conduct marine controlled source electromagnetic (CSEM) measurements in the Danube Delta, Black Sea, to identify suitable gas hydrate deposits for a future MEBO drilling test site for methane production and CO2 sequestration in hydrate form. Marine CSEM is a geophysical exploration method to derive the electrical properties, i.e. resistivity of the seafloor. Gas hydrates and free gas are electrically insulating and replace conductive pore fluid.