Closed-loop investigation of the dynamical properties of cardiomyocyte electrophysiology by means of Dynamic clamp and Mathematical modelling

The research field of my PhD concerns mathematical modeling and numerical simulation, applied to the cardiac electrophysiology analysis at a single cell level. This is possible thanks to the development of mathematical descriptions of single cellular components, ionic channels, pumps, exchangers and subcellular compartments. Due to the difficulties of vivo experiments on human cells, most of the measurements are acquired in vitro using animal models (e.g. guinea pig, dog, rabbit). Moreover, to study the cardiac action potential and all its features, it is necessary to acquire more specific knowledge about single ionic currents that contribute to the cardiac activity. Electrophysiological models of the heart have become very accurate in recent years giving rise to extremely complicated systems of differential equations. Although describing the behavior of cardiac cells quite well, the models are computationally demanding for numerical simulations and are very difficult to analyze from a mathematical (dynamical-systems) viewpoint. Simplified mathematical models that capture the underlying dynamics to a certain extent are therefore frequently used. The results presented in this thesis have confirmed that a close integration of computational modeling and experimental recordings in real myocytes, as performed by dynamic clamp, is a useful tool in enhancing our understanding of various components of normal cardiac electrophysiology, but also arrhythmogenic mechanisms in a pathological condition, especially when fully integrated with experimental data.

Modelling the dynamic response of rockfall protection barriers

Mountainous areas are prone to natural hazards like rockfalls. Among the many countermeasures, rockfall protection barriers represent an effective solution to mitigate the risk. They are metallic structures designed to intercept rocks falling from unstable slopes, thus dissipating the energy deriving from the impact. This study aims at providing a better understanding of the response of several rockfall barrier types, through the development of rather sophisticated three-dimensional numerical finite elements models which take into account for the highly dynamic and non-linear conditions of such events. The models are built considering the actual geometrical and mechanical properties of real systems. Particular attention is given to the connecting details between the structural components and to their interactions. The importance of the work lies in being able to support a wide experimental activity with appropriate numerical modelling. The data of several full-scale tests carried out on barrier prototypes, as well as on their structural components, are combined with results of numerical simulations. Though the models are designed with relatively simple solutions in order to obtain a low computational cost of the simulations, they are able to reproduce with great accuracy the test results, thus validating the reliability of the numerical strategy proposed for the design of these structures. The developed models have shown to be readily applied to predict the barrier performance under different possible scenarios, by varying the initial configuration of the structures and/or of the impact conditions. Furthermore, the numerical models enable to optimize the design of these structures and to evaluate the benefit of possible solutions. Finally it is shown they can be also used as a valuable supporting tool for the operators within a rockfall risk assessment procedure, to gain crucial understanding of the performance of existing barriers in working conditions.

Mechanical properties and structure of gel systems

Understanding the origins of the mechanical properties and its correlation withrnthe microstructure of gel systems is of great scientific and industrial interest. Inrngeneral, colloidal gels can be classified into chemical and physical gels, accordingrnto the life time of the network bonds. The characteristic di↵erences in gelationrndynamics can be observed with rheological measurements.rnAs a model system, a mixture of sodium silicate and low concentration sulfuric acidrnwas used. Nano-sized silica particles grow and aggregate to a system-spanning gelrnnetwork. The influence of the finite solubility of silica at high pH on the gelationrnwas studied with classical and piezo rheometer. The storage modulus of therngel grew logarithmically with time with two distinct growth laws. A relaxationrnat low frequency was observed in the frequency dependent measurements. I attributernthese two behaviors as a sign of structural rearrangements due to the finiternsolubility of silica at high pH. The reaction equilibrium between formation andrndissolution of bonds leads to a finite life time of the bonds and behavior similar tornphysical gel. The frequency dependence was more pronounced for lower water concentrations,rnhigher temperatures and shorter reaction times. With two relaxationrnmodels, I deduced characteristic relaxation times from the experimental data. Besidesrnrheology, the evolution of silica gels at high pH on di↵erent length scales wasrnstudied by NMR and dynamic light scattering. The results revealed that the primaryrnparticles existed already in sodium silicate and aggregated after the mixingrnof reactants due to a chemical reaction. Throughout the aggregation process thernsystem was in its chemical reaction equilibrium. Applying large oscillatory shearrnstrain to the gel allowed for modifying the gel modulus. The e↵ect of shear andrnshear history on the rheological properties of the gel were investigated. The storagernmodulus of the final gel increased with increasing strain. This behavior can be explained with (i) shear-induced aggregate compaction and (ii) combination ofrnbreakage and new formation of bonds.rnIn comparison with the physical gel-like behavior of the silica gel at high pH, typicalrnchemical gel features were exhibited by other gels formed from various chemicalrnreactions. Influences of the chemical structure modification on the gelation wererninvestigated with the piezo-rheometer. The external stimuli can be applied to tunernthe mechanical properties of the gel systems.

The dynamics of ultra-thin polymer films in different confinements studied by resonance enhanced dynamic light scattering

Die Kontroverse über den Glasübergang im Nanometerbereich, z. B. die Glas¬über¬gangs-temperatur Tg von dünnen Polymerfilmen, ist nicht vollständig abgeschlossen. Das dynamische Verhalten auf der Nanoskala ist stark von den einschränkenden Bedingungen abhängig, die auf die Probe wirken. Dünne Polymerfilme sind ideale Systeme um die Dynamik von Polymerketten unter der Einwirkung von Randbedingungen zu untersuchen, wie ich sie in dieser Arbeit variiert habe, um Einblick in dieses Problem zu erhalten.rnrnResonanzverstärkte dynamische Lichtstreuung ist eine Methode, frei von z.B. Fluoreszenzmarkern, die genutzt werden kann um in dünnen Polymerfilmen dynamische Phänomene

Synthesis and application of stimuli responsive tunable systems based on physical interactions

Novel single step synthetic procedure for hydrophobically modified alkali soluble latexes (HASE) via a miniemulsion-analogous method is presented. This facile method simplifies the copolymerization of the monomers with basically “opposite” character in terms of their hydrophilic/hydrophobic nature, which represent one of the main challenges in water based systems. Considered systems do not represent classical miniemulsions due to a high content of water soluble monomers. However, the polymerization mechanism was found to be rather similar to miniemulsion polymerization process.rnThe influence of the different factors on the system stability has been investigated. The copolymerization behavior studies typically showed strong composition drifts during copolymerization. It was found that the copolymer composition drift can be suppressed via changing the initial monomer ratio.rnThe neutralization behavior of the obtained HASE systems was investigated via potentiometric titration. The rheological behavior of the obtained systems as a function of the different parameters, such as pH, composition (ultrahydrophobe content) and additive type and content has also been investigated.rnDetailed investigation of the storage and loss moduli, damping factor and the crossover frequencies of the samples showed that at the initial stages of the neutralization the systems show microgel-like behavior.rnThe dependence of the rheological properties on the content and the type of the ultrahydrophobe showed that the tuning of the mechanical properties can be easily achieved via minor (few percent) but significant changes in the content of the latter. Besides, changing the hydrophobicity of the ultrahydrophobe via increasing the carbon chain length represents another simple method for achieving the same results.rnThe influence of amphiphilic additives (especially alcohols) on the rheological behavior of the obtained systems has been studied. An analogy was made between micellation of surfactants and the formation of hydrophobic domains between hydrophobic groups of the polymer side chain.rnDilution induced viscosity reduction was investigated in different systems, without or with different amounts or types of the amphiphilic additive. Possibility of the controlled response to dilution was explored. It was concluded that the sensitivity towards dilution can be reduced, and in extreme cases even the increase of the dynamic modulus can be observed, which is of high importance for the setting behavior of the adhesive material.rnIn the last part of this work, the adhesive behavior of the obtained HASE systems was investigated on different substrates (polypropylene and glass) for the standard labeling paper. Wet tack and setting behavior was studied and the trends for possible applications have been evaluated.rnThe novel synthetic procedure, investigation of rheological properties and the possibility of the tuning via additives, investigated in this work create a firm background for the development of the HASE based adhesives as well as rheology modifiers with vast variety of possible applications due to ease of tuning the mechanical and rheological properties of the systems.

Stochastic foundations of dynamic trade and labor market models

This dissertation consists of three self-contained papers that are related to two main topics. In particular, the first and third studies focus on labor market modeling, whereas the second essay presents a dynamic international trade setup.rnrnIn Chapter "Expenses on Labor Market Reforms during Transitional Dynamics", we investigate the arising costs of a potential labor market reform from a government point of view. To analyze various effects of unemployment benefits system changes, this chapter develops a dynamic model with heterogeneous employed and unemployed workers.rn rnIn Chapter "Endogenous Markup Distributions", we study how markup distributions adjust when a closed economy opens up. In order to perform this analysis, we first present a closed-economy general-equilibrium industry dynamics model, where firms enter and exit markets, and then extend our analysis to the open-economy case.rn rnIn Chapter "Unemployment in the OECD - Pure Chance or Institutions?", we examine effects of aggregate shocks on the distribution of the unemployment rates in OECD member countries.rn rnIn all three chapters we model systems that behave randomly and operate on stochastic processes. We therefore exploit stochastic calculus that establishes clear methodological links between the chapters.

Dynamic computer simulations of electrophoresis: a versatile research and teaching tool

Software is available, which simulates all basic electrophoretic systems, including moving boundary electrophoresis, zone electrophoresis, ITP, IEF and EKC, and their combinations under almost exactly the same conditions used in the laboratory. These dynamic models are based upon equations derived from the transport concepts such as electromigration, diffusion, electroosmosis and imposed hydrodynamic buffer flow that are applied to user-specified initial distributions of analytes and electrolytes. They are able to predict the evolution of electrolyte systems together with associated properties such as pH and conductivity profiles and are as such the most versatile tool to explore the fundamentals of electrokinetic separations and analyses. In addition to revealing the detailed mechanisms of fundamental phenomena that occur in electrophoretic separations, dynamic simulations are useful for educational purposes. This review includes a list of current high-resolution simulators, information on how a simulation is performed, simulation examples for zone electrophoresis, ITP, IEF and EKC and a comprehensive discussion of the applications and achievements.

Development of Dynamic Constraint Models for a Model Based Transient Calibration Process

Model based calibration has gained popularity in recent years as a method to optimize increasingly complex engine systems. However virtually all model based techniques are applied to steady state calibration. Transient calibration is by and large an emerging technology. An important piece of any transient calibration process is the ability to constrain the optimizer to treat the problem as a dynamic one and not as a quasi-static process. The optimized air-handling parameters corresponding to any instant of time must be achievable in a transient sense; this in turn depends on the trajectory of the same parameters over previous time instances. In this work dynamic constraint models have been proposed to translate commanded to actually achieved air-handling parameters. These models enable the optimization to be realistic in a transient sense. The air handling system has been treated as a linear second order system with PD control. Parameters for this second order system have been extracted from real transient data. The model has been shown to be the best choice relative to a list of appropriate candidates such as neural networks and first order models. The selected second order model was used in conjunction with transient emission models to predict emissions over the FTP cycle. It has been shown that emission predictions based on air-handing parameters predicted by the dynamic constraint model do not differ significantly from corresponding emissions based on measured air-handling parameters.

Investigating the Effects of Various Crowd Characteristics on the Dynamic Properties of an Occupied Structure

One of the challenges for structural engineers during design is considering how the structure will respond to crowd-induced dynamic loading. It has been shown that human occupants of a structure do not simply add mass to the system when considering the overall dynamic response of the system, but interact with it and may induce changes of the dynamic properties from those of the empty structure. This study presents an investigation into the human-structure interaction based on several crowd characteristics and their effect on the dynamic properties of an empty structure. The dynamic properties including frequency, damping, and mode shapes were estimated for a single test structure by means of experimental modal analysis techniques. The same techniques were utilized to estimate the dynamic properties when the test structure was occupied by a crowd with different combinations of size, posture, and distribution. The goal of this study is to isolate the occupant characteristics in order to determine the significance of each to be considered when designing new structures to avoid crowd serviceability issues. The results are presented and summarized based on the level of influence of each characteristic. The posture that produces the most significant effects based on the scope of this research is standing with bent knees with a maximum decrease in frequency of the first mode of the empty structure by 32 percent atthe highest mass ratio. The associated damping also increased 36 times the damping of the empty structure. In addition to the analysis of the experimental data, finite element models and a two degree-of-freedom model were created. These models were used to gain an understanding of the test structure, model a crowd as an equivalent mass, and also to develop a single degree-of-freedom (SDOF) model to best represent a crowd of occupants based on the experimental results. The SDOF models created had an averagefrequency of 5.0 Hz, within the range presented in existing biomechanics research, and combined SDOF systems of the test structure and crowd were able to reproduce the frequency and damping ratios associated with experimental tests. Results of this study confirmed the existence of human-structure interaction andthe inability to simply model a crowd as only additional mass. The two degree-offreedom model determined was able to predict the change in natural frequency and damping ratio for a structure occupied by multiple group sizes in a single posture. These results and model are the preliminary steps in the development of an appropriate methodfor modeling a crowd in combination with a more complex FE model of the empty structure.

A Systems-Level Approach to Human Epileptic Seizures

Epileptic seizures are due to the pathological collective activity of large cellular assemblies. A better understanding of this collective activity is integral to the development of novel diagnostic and therapeutic procedures. In contrast to reductionist analyses, which focus solely on small-scale characteristics of ictogenesis, here we follow a systems-level approach, which combines both small-scale and larger-scale analyses. Peri-ictal dynamics of epileptic networks are assessed by studying correlation within and between different spatial scales of intracranial electroencephalographic recordings (iEEG) of a heterogeneous group of patients suffering from pharmaco-resistant epilepsy. Epileptiform activity as recorded by a single iEEG electrode is determined objectively by the signal derivative and then subjected to a multivariate analysis of correlation between all iEEG channels. We find that during seizure, synchrony increases on the smallest and largest spatial scales probed by iEEG. In addition, a dynamic reorganization of spatial correlation is observed on intermediate scales, which persists after seizure termination. It is proposed that this reorganization may indicate a balancing mechanism that decreases high local correlation. Our findings are consistent with the hypothesis that during epileptic seizures hypercorrelated and therefore functionally segregated brain areas are re-integrated into more collective brain dynamics. In addition, except for a special sub-group, a highly significant association is found between the location of ictal iEEG activity and the location of areas of relative decrease of localised EEG correlation. The latter could serve as a clinically important quantitative marker of the seizure onset zone (SOZ).

A global ”imaging” view on systems approaches in immunology

The immune system exhibits an enormous complexity. High throughput methods such as the "-omic'' technologies generate vast amounts of data that facilitate dissection of immunological processes at ever finer resolution. Using high-resolution data-driven systems analysis, causal relationships between complex molecular processes and particular immunological phenotypes can be constructed. However, processes in tissues, organs, and the organism itself (so-called higher level processes) also control and regulate the molecular (lower level) processes. Reverse systems engineering approaches, which focus on the examination of the structure, dynamics and control of the immune system, can help to understand the construction principles of the immune system. Such integrative mechanistic models can properly describe, explain, and predict the behavior of the immune system in health and disease by combining both higher and lower level processes. Moving from molecular and cellular levels to a multiscale systems understanding requires the development of methodologies that integrate data from different biological levels into multiscale mechanistic models. In particular, 3D imaging techniques and 4D modeling of the spatiotemporal dynamics of immune processes within lymphoid tissues are central for such integrative approaches. Both dynamic and global organ imaging technologies will be instrumental in facilitating comprehensive multiscale systems immunology analyses as discussed in this review.

The new insight into dynamic crossover in glass forming liquids from the apparent enthalpy analysis

One of the most intriguing phenomena in glass forming systems is the dynamic crossover (T(B)), occurring well above the glass temperature (T(g)). So far, it was estimated mainly from the linearized derivative analysis of the primary relaxation time τ(T) or viscosity η(T) experimental data, originally proposed by Stickel et al. [J. Chem. Phys. 104, 2043 (1996); J. Chem. Phys. 107, 1086 (1997)]. However, this formal procedure is based on the general validity of the Vogel-Fulcher-Tammann equation, which has been strongly questioned recently [T. Hecksher et al. Nature Phys. 4, 737 (2008); P. Lunkenheimer et al. Phys. Rev. E 81, 051504 (2010); J. C. Martinez-Garcia et al. J. Chem. Phys. 134, 024512 (2011)]. We present a qualitatively new way to identify the dynamic crossover based on the apparent enthalpy space (H(a)(') = dlnτ/d(1/T)) analysis via a new plot lnH(a)(') vs. 1∕T supported by the Savitzky-Golay filtering procedure for getting an insight into the noise-distorted high order derivatives. It is shown that depending on the ratio between the "virtual" fragility in the high temperature dynamic domain (m(high)) and the "real" fragility at T(g) (the low temperature dynamic domain, m = m(low)) glass formers can be splitted into two groups related to f < 1 and f > 1, (f = m(high)∕m(low)). The link of this phenomenon to the ratio between the apparent enthalpy and activation energy as well as the behavior of the configurational entropy is indicated.

Predicting asthma control and exacerbations: chronic asthma as a complex dynamic model

PURPOSE OF REVIEW: Predicting asthma episodes is notoriously difficult but has potentially significant consequences for the individual, as well as for healthcare services. The purpose of this review is to describe recent insights into the prediction of acute asthma episodes in relation to classical clinical, functional or inflammatory variables, as well as present a new concept for evaluating asthma as a dynamically regulated homeokinetic system. RECENT FINDINGS: Risk prediction for asthma episodes or relapse has been attempted using clinical scoring systems, considerations of environmental factors and lung function, as well as inflammatory and immunological markers in induced sputum or exhaled air, and these are summarized here. We have recently proposed that newer mathematical methods derived from statistical physics may be used to understand the complexity of asthma as a homeokinetic, dynamic system consisting of a network comprising multiple components, and also to assess the risk for future asthma episodes based on fluctuation analysis of long time series of lung function. SUMMARY: Apart from the classical analysis of risk factor and functional parameters, this new approach may be used to assess asthma control and treatment effects in the individual as well as in future research trials.

Theoretical considerations to the verification of dynamic multileaf collimated fields with a SLIC-type portal image detector

The verification possibilities of dynamically collimated treatment beams with a scanning liquid ionization chamber electronic portal image device (SLIC-EPID) are investigated. The ion concentration in the liquid of a SLIC-EPID and therefore the read-out signal is determined by two parameters of a differential equation describing the creation and recombination of the ions. Due to the form of this equation, the portal image detector describes a nonlinear dynamic system with memory. In this work, the parameters of the differential equation were experimentally determined for the particular chamber in use and for an incident open 6 MV photon beam. The mathematical description of the ion concentration was then used to predict portal images of intensity-modulated photon beams produced by a dynamic delivery technique, the sliding window approach. Due to the nature of the differential equation, a mathematical condition for 'reliable leaf motion verification' in the sliding window technique can be formulated. It is shown that the time constants for both formation and decay of the equilibrium concentration in the chamber is in the order of seconds. In order to guarantee reliable leaf motion verification, these time constants impose a constraint on the rapidity of the image-read out for a given maximum leaf speed. For a leaf speed of 2 cm s(-1), a minimum image acquisition frequency of about 2 Hz is required. Current SLIC-EPID systems are usually too slow since they need about a second to acquire a portal image. However, if the condition is fulfilled, the memory property of the system can be used to reconstruct the leaf motion. It is shown that a simple edge detecting algorithm can be employed to determine the leaf positions. The method is also very robust against image noise.