Geant4 Hadronic Cascade Models and CMS Data Analysis : Computational Challenges in the LHC era

This work belongs to the field of computational high-energy physics (HEP). The key methods used in this thesis work to meet the challenges raised by the Large Hadron Collider (LHC) era experiments are object-orientation with software engineering, Monte Carlo simulation, the computer technology of clusters, and artificial neural networks. The first aspect discussed is the development of hadronic cascade models, used for the accurate simulation of medium-energy hadron-nucleus reactions, up to 10 GeV. These models are typically needed in hadronic calorimeter studies and in the estimation of radiation backgrounds. Various applications outside HEP include the medical field (such as hadron treatment simulations), space science (satellite shielding), and nuclear physics (spallation studies). Validation results are presented for several significant improvements released in Geant4 simulation tool, and the significance of the new models for computing in the Large Hadron Collider era is estimated. In particular, we estimate the ability of the Bertini cascade to simulate Compact Muon Solenoid (CMS) hadron calorimeter HCAL. LHC test beam activity has a tightly coupled cycle of simulation-to-data analysis. Typically, a Geant4 computer experiment is used to understand test beam measurements. Thus an another aspect of this thesis is a description of studies related to developing new CMS H2 test beam data analysis tools and performing data analysis on the basis of CMS Monte Carlo events. These events have been simulated in detail using Geant4 physics models, full CMS detector description, and event reconstruction. Using the ROOT data analysis framework we have developed an offline ANN-based approach to tag b-jets associated with heavy neutral Higgs particles, and we show that this kind of NN methodology can be successfully used to separate the Higgs signal from the background in the CMS experiment.

Morphogenesis of Mammalian Endoplasmic Reticulum and Golgi Apparatus throughout the Cell Cycle

The endoplasmic reticulum (ER) and the Golgi apparatus are organelles that produce, modify and transport proteins and lipids and regulate Ca2+ environment within cells. Structurally they are composed of sheets and tubules. Sheets may take various forms: intact, fenestrated, single or stacked. The ER, including the nuclear envelope, is a single continuous network, while the Golgi shows only some level of connectivity. It is often unclear, how different morphologies correspond to particular functions. Previous studies indicate that the structures of the ER and Golgi are dynamic and regulated by fusion and fission events, cytoskeleton, rate of protein synthesis and secretion, and specific structural proteins. For example, many structural proteins shaping tubular ER have been identified, but sheet formation is much more unclear. In this study, we used light and electron microscopy to study morphological changes of the ER and Golgi in mammalian cells. The proportion, type, location and dynamics of ER sheets and tubules were found to vary in a cell type or cell cycle stage dependent manner. During interphase, ER and Golgi structures were demonstrated to be regulated by p37, a cofactor of the fusion factor p97, and microtubules, which also affected the localization of the organelles. Like previously shown for the Golgi, the ER displayed a tendency for fenestration and tubulation during mitosis. However, this shape change did not result in ER fragmentation as happens to Golgi, but a continuous network was retained. The activity of p97/p37 was found to be important for the reassembly of both organelles after mitosis. In EM images, ER sheet membranes appear rough, since they contain attached ribosomes, whereas tubular membranes appear smooth. Our studies revealed that structural changes of the ER towards fenestrated and tubular direction correlate with loss of ER-bound ribosomes and vice versa. High and low curvature ER membranes have a low and high density of ribosomes, respectively. To conclude, both ER and Golgi architecture depend on fusion activity of p97/p37. ER morphogenesis, particularly of the sheet shape, is intimately linked to the density of membrane bound ribosomes.

Methods in general model localization

The aim of this study was to evaluate and test methods which could improve local estimates of a general model fitted to a large area. In the first three studies, the intention was to divide the study area into sub-areas that were as homogeneous as possible according to the residuals of the general model, and in the fourth study, the localization was based on the local neighbourhood. According to spatial autocorrelation (SA), points closer together in space are more likely to be similar than those that are farther apart. Local indicators of SA (LISAs) test the similarity of data clusters. A LISA was calculated for every observation in the dataset, and together with the spatial position and residual of the global model, the data were segmented using two different methods: classification and regression trees (CART) and the multiresolution segmentation algorithm (MS) of the eCognition software. The general model was then re-fitted (localized) to the formed sub-areas. In kriging, the SA is modelled with a variogram, and the spatial correlation is a function of the distance (and direction) between the observation and the point of calculation. A general trend is corrected with the residual information of the neighbourhood, whose size is controlled by the number of the nearest neighbours. Nearness is measured as Euclidian distance. With all methods, the root mean square errors (RMSEs) were lower, but with the methods that segmented the study area, the deviance in single localized RMSEs was wide. Therefore, an element capable of controlling the division or localization should be included in the segmentation-localization process. Kriging, on the other hand, provided stable estimates when the number of neighbours was sufficient (over 30), thus offering the best potential for further studies. Even CART could be combined with kriging or non-parametric methods, such as most similar neighbours (MSN).

The function of NR3B and NR4A orphan nuclear receptors in osteoblasts

Nuclear receptors (NRs) comprise a large family of proteins that mediate the effects of small lipophilic molecules such as steroid hormones. In addition, there are a group of NRs which lack identified natural ligands and are referred as orphan NRs. In this thesis, the function of two such orphan NR families, the NR3B (ERRα, ERRβ and ERRγ) and the NR4A family (NGFI-B, Nurr1 and Nor1), was studied. NR3B and NR4A receptors regulate many biological processes such as energy metabolism and carcinogenesis. In addition, NR3B and NR4A receptors are expressed in bone. Therefore, the signaling and function of NR3B and NR4A orphan nuclear receptors was studied specifically in osteoblasts. NR4A receptors were found to be regulated by NR3B receptors and the Wnt/β-catenin signaling pathway as ERRα, ERRγ and β-catenin repressed the transcriptional activity of NR4A receptors in U2-OS cells. NGFI-B was found to repress the transcriptional activity of ERRγ in HeLa cells. The phytoestrogen equol was identified as a new agonist for ERRγ and ERRβ in PC-3, U2-OS, and SaOS-2 cells. Equol increased the transcriptional activity of ERRγ by increasing ERRγ co-activator binding and by inducing a conformational change in the ligand binding pocket of ERRγ. The growth inhibitory effect of equol on PC-3 prostate cancer cells was decreased by blocking ERRγ expression by siRNA. Therefore, ERRγ could mediate some of the beneficial health effects of equol. The Wnt/β-catenin signaling pathway is important for the differentiation and function of osteoblasts. NR3B and NR4A receptors were found to repress the transcriptional activity mediated by β-catenin in U2-OS cells. The mesenchymal stem cells (MSCs) isolated from ERRα knockout (KO) mice showed diminished proliferation and osteoblastic differentiation compared to the wild-type cells. The overexpression of ERRα in osteoblastic MC3T3-E1 cell line increased their mineralization. Bone sialoprotein (BSP) was shown to be a direct target gene for ERRα and ERRγ as the BSP promoter was activated by ERRα or ERRγ and PGC-1α in HeLa cells. The adipogenic differentiation of ERRα KO MSCs was also decreased and they expressed less adipogenic marker genes. In conclusion, the studies described in this thesis demonstrated that the transcriptional activity of NR3B and NR4A receptors can be regulated by other orphan NRs and signaling pathways in osteoblasts. NR3B receptors can also be regulated by ligands and a new agonist, equol, was identified for ERRβ and ERRγ. New roles for NR3B and NR4A were also identified as they were shown to converge with the Wnt signaling pathway in osteoblasts, ERRγ was shown to mediate the growth inhibitory effect of equol in prostate cancer cells, and ERRα was shown to regulate positively MSC proliferation, osteoblastic differentiation and adipogenesis.

mRNA-binding protein HuR in breast carcinogenesis

Breast cancer is the most common cancer among women. Although its prognosis has improved nowadays, methods to predict the progression of the disease or to treat it are not comprehensive. This thesis work was initiated to elucidate in breast carcinogenesis the role of HuR, a ubiquitously expressed mRNA-binding protein that regulates gene expression posttranscriptionally. HuR is predominantly nuclear, but it shuttles between the nucleus and the cytoplasm, and this nucleocytoplasmic translocation is important for its function as a RNA-stabilizing and translational regulator. HuR has been associated with diverse cellular processes, for example carcinogenesis. The specific aims of my thesis work were to study the prognostic value of HuR in breast cancer and to clarify the mechanisms by which HuR contributes to breast carcinogenesis. My ultimate goal is, by better understanding the role of HuR in breast carcinogenesis, to aid in the discovery of novel targets for cancer therapies. HuR expression and localization was studied in paraffin-embedded preinvasive (atypical ductal hyperplasia, ADH, and ductal carcinoma in situ, DCIS) specimens as well in sporadic and familial breast cancer specimens. Our results show that cytoplasmic HuR expression was already elevated in ADH and remained elevated in DCIS as well as in cancer specimens. Clinicopathological analysis showed that cytoplasmic HuR expression associated with the more aggressive form of the disease in DCIS, and in cancer specimens it proved an independent marker for poor prognosis. Importantly, cytoplasmic HuR expression was significantly associated with poor outcome in the subgroups of small (2 cm) and axillary lymph node-negative breast cancers. HuR proved to be the first mRNA stability protein the expression of which is associated in breast cancer with poor outcome. To explore the mechanisms of HuR in breast carcinogenesis, lentiviral constructs were developed to inhibit and to overexpress the HuR expression in a breast epithelial cell line (184B5Me). Our results suggest that HuR mediates breast carcinogenesis by participating in processes important in cell transformation, in programmed cell death, and in cell invasion. Global gene expression analysis shows that HuR regulates genes participating in diverse cellular processes, and affects several pathways important in cancer development. In addition, we identified two novel target transcripts (connective tissue growth factor, CTGF, and Ras oncogene family member 31, RAB31) for HuR. In conclusion, because cytoplasmic HuR expression in breast cancer can predict the outcome of the disease it could serve in clinics as a prognostic marker. HuR accumulates in the cytoplasm even at its non-invasive stage (ADH and DCIS) of the carcinogenic process and supports functions essential in cell alteration. These data suggest that HuR contributes to carcinogenesis of the breast epithelium.

Hyperfine effects in the nuclear magnetic resonance of paramagnetic molecules

Paramagnetic, or open-shell, systems are often encountered in the context of metalloproteins, and they are also an essential part of molecular magnets. Nuclear magnetic resonance (NMR) spectroscopy is a powerful tool for chemical structure elucidation, but for paramagnetic molecules it is substantially more complicated than in the diamagnetic case. Before the present work, the theory of NMR of paramagnetic molecules was limited to spin-1/2 systems and it did not include relativistic corrections to the hyperfine effects. It also was not systematically expandable. --- The theory was first expanded by including hyperfine contributions up to the fourth power in the fine structure constant α. It was then reformulated and its scope widened to allow any spin state in any spatial symmetry. This involved including zero-field splitting effects. In both stages the theory was implemented into a separate analysis program. The different levels of theory were tested by demonstrative density functional calculations on molecules selected to showcase the relative strength of new NMR shielding terms. The theory was also tested in a joint experimental and computational effort to confirm assignment of 11 B signals. The new terms were found to be significant and comparable with the terms in the earlier levels of theory. The leading-order magnetic-field dependence of shielding in paramagnetic systems was formulated. The theory is now systematically expandable, allowing for higher-order field dependence and relativistic contributions. The prevailing experimental view of pseudocontact shift was found to be significantly incomplete, as it only includes specific geometric dependence, which is not present in most of the new terms introduced here. The computational uncertainty in density functional calculations of the Fermi contact hyperfine constant and zero-field splitting tensor sets a limit for quantitative prediction of paramagnetic shielding for now.

Functional expression and subcellular localization of the Cl- cotransporters KCC2 and NKCC1 in rodent hippocampal and neocortical neurons

The work presented here has focused on the role of cation-chloride cotransporters (CCCs) in (1) the regulation of intracellular chloride concentration within postsynaptic neurons and (2) on the consequent effects on the actions of the neurotransmitter gamma-aminobutyric acid (GABA) mediated by GABAA receptors (GABAARs) during development and in pathophysiological conditions such as epilepsy. In addition, (3) we found that a member of the CCC family, the K-Cl cotransporter isoform 2 (KCC2), has a structural role in the development of dendritic spines during the differentiation of pyramidal neurons. Despite the large number of publications dedicated to regulation of intracellular Cl-, our understanding of the underlying mechanisms is not complete. Experiments on GABA actions under resting steady-state have shown that the effect of GABA shifts from depolarizing to hyperpolarizing during maturation of cortical neurons. However, it remains unclear, whether conclusions from these steady-state measurements can be extrapolated to the highly dynamic situation within an intact and active neuronal network. Indeed, GABAergic signaling in active neuronal networks results in a continuous Cl- load, which must be constantly removed by efficient Cl- extrusion mechanisms. Therefore, it seems plausible to suggest that key parameters are the efficacy and subcellular distribution of Cl- transporters rather than the polarity of steady-state GABA actions. A further related question is: what are the mechanisms of Cl- regulation and homeostasis during pathophysiological conditions such as epilepsy in adults and neonates? Here I present results that were obtained by means of a newly developed method of measurements of the efficacy of a K-Cl cotransport. In Study I, the developmental profile of KCC2 functionality during development was analyzed both in dissociated neuronal cultures and in acute hippocampal slices. A novel method of photolysis of caged GABA in combination with Cl- loading to the somata was used in this study to assess the extrusion efficacy of KCC2. We demonstrated that these two preparations exhibit a different temporal profile of functional KCC2 upregulation. In Study II, we reported an observation of highly distorted dendritic spines in neurons cultured from KCC2-/- embryos. During their development in the culture dish, KCC2-lacking neurons failed to develop mature, mushroom-shaped dendritic spines but instead maintained an immature phenotype of long, branching and extremely motile protrusions. It was shown that the role of KCC2 in spine maturation is not based on its transport activity, but is mediated by interactions with cytoskeletal proteins. Another important player in Cl- regulation, NKCC1 and its role in the induction and maintenance of native Cl- gradients between the axon initial segment (AIS) and soma was the subject of Study III. There we demonstrated that this transporter mediates accumulation of Cl- in the axon initial segment of neocortical and hippocampal principal neurons. The results suggest that the reversal potential of the GABAA response triggered by distinct populations of interneurons show large subcellular variations. Finally, a novel mechanism of fast post-translational upregulation of the membrane-inserted, functionally active KCC2 pool during in-vivo neonatal seizures and epileptiform-like activity in vitro was identified and characterized in Study IV. The seizure-induced KCC2 upregulation may act as an intrinsic antiepileptogenic mechanism.

The many ways to delimit species: hairs, genes and surface chemistry

Species identification forms the basis for understanding the diversity of the living world, but it is also a prerequisite for understanding many evolutionary patterns and processes. The most promising approach for correctly delimiting and identifying species is to integrate many types of information in the same study. Our aim was to test how cuticular hydro- carbons, traditional morphometrics, genetic polymorphisms in nuclear markers (allozymes and DNA microsatellites) and DNA barcoding (partial mitochondrial COI gene) perform in delimiting species. As an example, we used two closely related Formica ants, F. fusca and F. lemani, sampled from a sympatric population in the northern part of their distribu- tion. Morphological characters vary and overlap in different parts of their distribution areas, but cuticular hydrocarbons include a strong taxonomic signal and our aim is to test the degree to which morphological and genetic data correspond to the chemical data. In the morphological analysis, species were best separated by the combined number of hairs on pro- notum and mesonotum, but individual workers overlapped in hair numbers, as previously noted by several authors. Nests of the two species were separated but not clustered according to species in a Principal Component Analysis made on nuclear genetic data. However, model-based Bayesian clustering resulted in perfect separation of the species and gave no indication of hybridization. Furthermore, F. lemani and F. fusca did not share any mitochondrial haplotypes, and the species were perfectly separated in a phylogenetic tree. We conclude that F. fusca and F. lemani are valid species that can be separated in our study area relatively well with all methods employed. However, the unusually small genetic differen- tiation in nuclear markers (FST = 0.12) shows that they are closely related, and occasional hybridization between F. fusca and F. lemani cannot be ruled out.

Moonlighting Proteins of Lactobacillus crispatus : Extracellular Localization, Cell Wall Anchoring and Interactions with the Host

Moonlighting functions have been described for several proteins previously thought to localize exclusively in the cytoplasm of bacterial or eukaryotic cells. Moonlighting proteins usually perform conserved functions, e. g. in glycolysis or as chaperonins, and their traditional and moonlighting function(s) usually localize to different cell compartments. The most characterized moonlighting proteins in Grampositive bacteria are the glycolytic enzymes enolase and glyceraldehyde-3-phosphate dehydrogenase (GAPDH), which function in bacteria-host interactions, e. g. as adhesins or plasminogen receptors. Research on bacterial moonlighting proteins has focused on Gram-positive bacterial pathogens, where many of their functions have been associated with bacterial virulence. In this thesis work I show that also species of the genus Lactobacillus have moonlighting proteins that carry out functions earlier associated with bacterial virulence only. I identified enolase, GAPDH, glutamine synthetase (GS), and glucose-6-phosphate isomerase (GPI) as moonlighting proteins of Lactobacillus crispatus strain ST1 and demonstrated that they are associated with cell surface and easily released from the cell surface into incubation buffer. I also showed that these lactobacillar proteins moonlight either as adhesins with affinity for basement membrane and extracellular matrix proteins or as plasminogen receptors. The mechanisms of surface translocation and anchoring of bacterial moonlighting proteins have remained enigmatic. In this work, the surface localization of enolase, GAPDH, GS and GPI was shown to depend on environmental factors. The members of the genus Lactobacillus are fermentative organisms that lower the ambient pH by producing lactic acid. At acidic pH enolase, GAPDH, GS and GPI were associated with the cell surface, whereas at neutral pH they were released into the buffer. The release did not involve de novo protein synthesis. I showed that purified recombinant His6-enolase, His6-GAPDH, His6-GS and His6-GPI reassociate with cell wall and bind in vitro to lipoteichoic acids at acidic pH. The in-vitro binding of these proteins localizes to cell division septa and cell poles. I also show that the release of moonlighting proteins is enhanced in the presence of cathelicidin LL- 37, which is an antimicrobial peptide and a central part of the innate immunity defence. I found that the LL-37-induced detachment of moonlighting proteins from cell surface is associated with cell wall permeabilization by LL-37. The results in this thesis work are compatible with the hypothesis that the moonlighting proteins of L. crispatus associate to the cell wall via electrostatic or ionic interactions and that they are released into surroundings in stress conditions. Their surface translocation is, at least in part, a result from their release from dead or permeabilized cells and subsequent reassociation onto the cell wall. The results of this thesis show that lactobacillar cells rapidly change their surface architecture in response to environmental factors and that these changes influence bacterial interactions with the host.