Hybrid Modelling of Titan's Interaction with the Magnetosphere of Saturn

In this dissertation we study the interaction between Saturn's moon Titan and the magnetospheric plasma and magnetic field. The method of research is a three-dimensional computer simulation model, that is used to simulate this interaction. The simulation model used is a hybrid model. Hybrid models enable individual tracking or tracing of ions and also take into account the particle motion in the propagation of the electromagnetic fields. The hybrid model has been developed at the Finnish Meteorological Institute. This thesis gives a general description of the effects that the solar wind has on Earth and other planets of our solar system. Planetary satellites can also have similar interactions with the solar wind but also with the plasma flows of planetary magnetospheres. Titan is clearly the largest among the satellites of Saturn and also the only known satellite with a dense atmosphere. It is the atmosphere that makes Titan's plasma interaction with the magnetosphere of Saturn so unique. Nevertheless, comparisons with the plasma interactions of other solar system bodies are valuable. Detecting charged plasma particles requires in situ measurements obtainable through scientific spacecraft. The Cassini mission has been one of the most remarkable international efforts in space science. Since 2004 the measurements and images obtained from instruments onboard the Cassini spacecraft have increased the scientific knowledge of Saturn as well as its satellites and magnetosphere in a way no one was probably able to predict. The current level of science on Titan is practically unthinkable without the Cassini mission. Many of the observations by Cassini instrument teams have influenced this research both the direct measurements of Titan as well as observations of its plasma environment. The theoretical principles of the hybrid modelling approach are presented in connection to the broader context of plasma simulations. The developed hybrid model is described in detail: e.g. the way the equations of the hybrid model are solved is shown explicitly. Several simulation techniques, such as the grid structure and various boundary conditions, are discussed in detail as well. The testing and monitoring of simulation runs is presented as an essential routine when running sophisticated and complex models. Several significant improvements of the model, that are in preparation, are also discussed. A main part of this dissertation are four scientific articles based on the results of the Titan model. The Titan model developed during the course of the Ph.D. research has been shown to be an important tool to understand Titan's plasma interaction. One reason for this is that the structures of the magnetic field around Titan are very much three-dimensional. The simulation results give a general picture of the magnetic fields in the vicinity of Titan. The magnetic fine structure of Titan's wake as seen in the simulations seems connected to Alfvén waves an important wave mode in space plasmas. The particle escape from Titan is also a major part of these studies. Our simulations show a bending or turning of Titan's ionotail that we have shown to be a direct result of the basic principles in plasma physics. Furthermore, the ion flux from the magnetosphere of Saturn into Titan's upper atmosphere has been studied. The modelled ion flux has asymmetries that would likely have a large impact in the heating in different parts of Titan's upper atmosphere.

Characterization of Ion-Exchange Fibers for Controlled Drug Delivery

Increasing attention has been focused on methods that deliver pharmacologically active compounds (e.g. drugs, peptides and proteins) in a controlled fashion, so that constant, sustained, site-specific or pulsatile action can be attained. Ion-exchange resins have been widely studied in medical and pharmaceutical applications, including controlled drug delivery, leading to commercialisation of some resin based formulations. Ion-exchangers provide an efficient means to adjust and control drug delivery, as the electrostatic interactions enable precise control of the ion-exchange process and, thus, a more uniform and accurate control of drug release compared to systems that are based only on physical interactions. Unlike the resins, only few studies have been reported on ion-exchange fibers in drug delivery. However, the ion-exchange fibers have many advantageous properties compared to the conventional ion-exchange resins, such as more efficient compound loading into and release from the ion-exchanger, easier incorporation of drug-sized compounds, enhanced control of the ion-exchange process, better mechanical, chemical and thermal stability, and good formulation properties, which make the fibers attractive materials for controlled drug delivery systems. In this study, the factors affecting the nature and strength of the binding/loading of drug-sized model compounds into the ion-exchange fibers was evaluated comprehensively and, moreover, the controllability of subsequent drug release/delivery from the fibers was assessed by modifying the conditions of external solutions. Also the feasibility of ion-exchange fibers for simultaneous delivery of two drugs in combination was studied by dual loading. Donnan theory and theoretical modelling were applied to gain mechanistic understanding on these factors. The experimental results imply that incorporation of model compounds into the ion-exchange fibers was attained mainly as a result of ionic bonding, with additional contribution of non-specific interactions. Increasing the ion-exchange capacity of the fiber or decreasing the valence of loaded compounds increased the molar loading, while more efficient release of the compounds was observed consistently at conditions where the valence or concentration of the extracting counter-ion was increased. Donnan theory was capable of fully interpreting the ion-exchange equilibria and the theoretical modelling supported precisely the experimental observations. The physico-chemical characteristics (lipophilicity, hydrogen bonding ability) of the model compounds and the framework of the fibrous ion-exchanger influenced the affinity of the drugs towards the fibers and may, thus, affect both drug loading and release. It was concluded that precisely controlled drug delivery may be tailored for each compound, in particularly, by choosing a suitable ion-exchange fiber and optimizing the delivery system to take into account the external conditions, also when delivering two drugs simultaneously.

Algorithms for 13C metabolic flux analysis

The metabolism of an organism consists of a network of biochemical reactions that transform small molecules, or metabolites, into others in order to produce energy and building blocks for essential macromolecules. The goal of metabolic flux analysis is to uncover the rates, or the fluxes, of those biochemical reactions. In a steady state, the sum of the fluxes that produce an internal metabolite is equal to the sum of the fluxes that consume the same molecule. Thus the steady state imposes linear balance constraints to the fluxes. In general, the balance constraints imposed by the steady state are not sufficient to uncover all the fluxes of a metabolic network. The fluxes through cycles and alternative pathways between the same source and target metabolites remain unknown. More information about the fluxes can be obtained from isotopic labelling experiments, where a cell population is fed with labelled nutrients, such as glucose that contains 13C atoms. Labels are then transferred by biochemical reactions to other metabolites. The relative abundances of different labelling patterns in internal metabolites depend on the fluxes of pathways producing them. Thus, the relative abundances of different labelling patterns contain information about the fluxes that cannot be uncovered from the balance constraints derived from the steady state. The field of research that estimates the fluxes utilizing the measured constraints to the relative abundances of different labelling patterns induced by 13C labelled nutrients is called 13C metabolic flux analysis. There exist two approaches of 13C metabolic flux analysis. In the optimization approach, a non-linear optimization task, where candidate fluxes are iteratively generated until they fit to the measured abundances of different labelling patterns, is constructed. In the direct approach, linear balance constraints given by the steady state are augmented with linear constraints derived from the abundances of different labelling patterns of metabolites. Thus, mathematically involved non-linear optimization methods that can get stuck to the local optima can be avoided. On the other hand, the direct approach may require more measurement data than the optimization approach to obtain the same flux information. Furthermore, the optimization framework can easily be applied regardless of the labelling measurement technology and with all network topologies. In this thesis we present a formal computational framework for direct 13C metabolic flux analysis. The aim of our study is to construct as many linear constraints to the fluxes from the 13C labelling measurements using only computational methods that avoid non-linear techniques and are independent from the type of measurement data, the labelling of external nutrients and the topology of the metabolic network. The presented framework is the first representative of the direct approach for 13C metabolic flux analysis that is free from restricting assumptions made about these parameters.In our framework, measurement data is first propagated from the measured metabolites to other metabolites. The propagation is facilitated by the flow analysis of metabolite fragments in the network. Then new linear constraints to the fluxes are derived from the propagated data by applying the techniques of linear algebra.Based on the results of the fragment flow analysis, we also present an experiment planning method that selects sets of metabolites whose relative abundances of different labelling patterns are most useful for 13C metabolic flux analysis. Furthermore, we give computational tools to process raw 13C labelling data produced by tandem mass spectrometry to a form suitable for 13C metabolic flux analysis.

Ion-regulatory proreins in neuronal development and communication

Brain function is critically dependent on the ionic homeostasis in both the extra- and intracellular compartment. The regulation of brain extracellular ionic composition mainly relies on active transport at blood brain and at blood cerebrospinal fluid interfaces whereas intracellular ion regulation is based on plasmalemmal transporters of neurons and glia. In addition, the latter mechanisms can generate physiologically as well as pathophysiologically significant extracellular ion transients. In this work I have studied molecular mechanisms and development of ion regulation and how these factors alter neuronal excitability and affect synaptic and non-synaptic transmission with a particular emphasis on intracellular pH and chloride (Cl-) regulation. Why is the regulation of acid-base equivalents (H+ and HCO3-) and Cl- of such interest and importance? First of all, GABAA-receptors are permeable to both HCO3- and Cl-. In the adult mammalian central nervous system (CNS) fast postsynaptic inhibition relies on GABAA-receptor mediated transmission. Today, excitatory effects of GABAA-receptors, both in mature neurons and during the early development, have been recognized and the significance of the dual actions of GABA on neuronal communication has become an interesting field of research. The transmembrane gradients of Cl- and HCO3- determine the reversal potential of GABAA-receptor mediated postsynaptic potentials and hence, the function of pH and Cl- regulatory proteins have profound consequences on GABAergic signaling and neuronal excitability. Secondly, perturbations in pH can cause a variety of changes in cellular function, many of them resulting from the interaction of protons with ionizable side chains of proteins. pH-mediated alterations of protein conformation in e.g. ion channels, transporters, and enzymes can powerfully modulate neurotransmission. In the context of pH homeostasis, the enzyme carbonic anhydrase (CA) needs to be taken into account in parallel with ion transporters: for CO2/HCO3- buffering to act in a fast manner, CO2 (de)hydration must be catalyzed by this enzyme. The acid-base equivalents that serve as substrates in the CO2 dehydration-hydration reaction are also engaged in many carrier and channel mediated ion movements. In such processes, CA activity is in key position to modulate transmembrane solute fluxes and their consequences. The bicarbonate transporters (BTs; SLC4) and the electroneutral cation-chloride cotransporters (CCCs; SLC12) belong the to large gene family of solute carriers (SLCs). In my work I have studied the physiological roles of the K+-Cl- cotransporter KCC2 (Slc12a5) and the Na+-driven Cl--HCO3- exchanger NCBE (Slc4a10) and the roles of these two ion transporters in the modualtion of neuronal communication and excitability in the rodent hippocampus. I have also examined the cellular localization and molecular basis of intracellular CA that has been shown to be essential for the generation of prolonged GABAergic excitation in the mature hippocampus. The results in my Thesis provide direct evidence for the view that the postnatal up-regulation of KCC2 accounts for the developmental shift from depolarizing to hyperpolarizing postsynaptic EGABA-A responses in rat hippocampal pyramidal neurons. The results also indicate that after KCC2 expression the developmental onset of excitatory GABAergic transmission upon intense GABAA-receptor stimulation depend on the expression of intrapyramidal CA, identified as the CA isoform VII. Studies on mice with targeted Slc4a10 gene disruption revealed an important role for NCBE in neuronal pH regulation and in pH-dependent modulation of neuronal excitability. Furthermore, this ion transporter is involved in the basolateral Na+ and HCO3- uptake in choroid plexus epithelial cells, and is thus likely to contribute to cerebrospinal fluid production.

Ion channel gene variants predisposing to severe human ventricular arrhythmias

Congenital long QT syndrome (LQTS) with an estimated prevalence of 1:2000-1:10 000 manifests with prolonged QT interval on electrocardiogram and risk for ventricular arrhythmias and sudden death. Several ion channel genes and hundreds of mutations in these genes have been identified to underlie the disorder. In Finland, four LQTS founder mutations of potassium channel genes account for up to 40-70% of genetic spectrum of LQTS. Acquired LQTS has similar clinical manifestations, but often arises from usage of QT-prolonging medication or electrolyte disturbances. A prolonged QT interval is associated with increased morbidity and mortality not only in clinical LQTS but also in patients with ischemic heart disease and in the general population. The principal aim of this study was to estimate the actual prevalence of LQTS founder mutations in Finland and to calculate their effect on QT interval in the Finnish background population. Using a large population-based sample of over 6000 Finnish individuals from the Health 2000 Survey, we identified LQTS founder mutations KCNQ1 G589D (n=8), KCNQ1 IVS7-2A>G (n=1), KCNH2 L552S (n=2), and KCNH2 R176W (n=16) in 27 study participants. This resulted in a weighted prevalence estimate of 0.4% for LQTS in Finland. Using a linear regression model, the founder mutations resulted in a 22- to 50-ms prolongation of the age-, sex-, and heart rate-adjusted QT interval. Collectively, these data suggest that one of 250 individuals in Finland may be genetically predisposed to ventricular arrhythmias arising from the four LQTS founder mutations. A KCNE1 D85N minor allele with a frequency of 1.4% was associated with a 10-ms prolongation in adjusted QT interval and could thus identify individuals at increased risk of ventricular arrhythmias at the population level. In addition, the previously reported associations of KCNH2 K897T, KCNH2 rs3807375, and NOS1AP rs2880058 with QT interval duration were confirmed in the present study. In a separate study, LQTS founder mutations were identified in a subgroup of acquired LQTS, providing further evidence that congenital LQTS gene mutations may underlie acquired LQTS. Catecholaminergic polymorphic ventricular tachycardia (CPVT) is characterized by exercise-induced ventricular arrhythmias in a structurally normal heart and results from defects in the cardiac Ca2+ signaling proteins, mainly ryanodine receptor type 2 (RyR2). In a patient population of typical CPVT, RyR2 mutations were identifiable in 25% (4/16) of patients, implying that noncoding variants or other genes are involved in CPVT pathogenesis. A 1.1 kb RyR2 exon 3 deletion was identified in two patients independently, suggesting that this region may provide a new target for RyR2-related molecular genetic studies. Two novel RyR2 mutations showing a gain-of-function defect in vitro were identified in three victims of sudden cardiac death. Extended pedigree analyses revealed some surviving mutation carriers with mild structural abnormalities of the heart and resting ventricular arrhythmias suggesting that not all RyR2 mutations lead to a typical CPVT phenotype, underscoring the relevance of tailored risk stratification of a RyR2 mutation carrier.

Applications of Gauge/Gravity Duality to QCD and Heavy Ion Collisions

The description of quarks and gluons, using the theory of quantum chromodynamics (QCD), has been known for a long time. Nevertheless, many fundamental questions in QCD remain unanswered. This is mainly due to problems in solving the theory at low energies, where the theory is strongly interacting. AdS/CFT is a duality between a specific string theory and a conformal field theory. Duality provides new tools to solve the conformal field theory in the strong coupling regime. There is also some evidence that using the duality, one can get at least qualitative understanding of how QCD behaves at strong coupling. In this thesis, we try to address some issues related to QCD and heavy ion collisions, applying the duality in various ways.