Computing the Stochastic Complexity of Simple Probabilistic Graphical Models

Minimum Description Length (MDL) is an information-theoretic principle that can be used for model selection and other statistical inference tasks. There are various ways to use the principle in practice. One theoretically valid way is to use the normalized maximum likelihood (NML) criterion. Due to computational difficulties, this approach has not been used very often. This thesis presents efficient floating-point algorithms that make it possible to compute the NML for multinomial, Naive Bayes and Bayesian forest models. None of the presented algorithms rely on asymptotic analysis and with the first two model classes we also discuss how to compute exact rational number solutions.

A Probabilistic Approach to the Primary Visual Cortex

What can the statistical structure of natural images teach us about the human brain? Even though the visual cortex is one of the most studied parts of the brain, surprisingly little is known about how exactly images are processed to leave us with a coherent percept of the world around us, so we can recognize a friend or drive on a crowded street without any effort. By constructing probabilistic models of natural images, the goal of this thesis is to understand the structure of the stimulus that is the raison d etre for the visual system. Following the hypothesis that the optimal processing has to be matched to the structure of that stimulus, we attempt to derive computational principles, features that the visual system should compute, and properties that cells in the visual system should have. Starting from machine learning techniques such as principal component analysis and independent component analysis we construct a variety of sta- tistical models to discover structure in natural images that can be linked to receptive field properties of neurons in primary visual cortex such as simple and complex cells. We show that by representing images with phase invariant, complex cell-like units, a better statistical description of the vi- sual environment is obtained than with linear simple cell units, and that complex cell pooling can be learned by estimating both layers of a two-layer model of natural images. We investigate how a simplified model of the processing in the retina, where adaptation and contrast normalization take place, is connected to the nat- ural stimulus statistics. Analyzing the effect that retinal gain control has on later cortical processing, we propose a novel method to perform gain control in a data-driven way. Finally we show how models like those pre- sented here can be extended to capture whole visual scenes rather than just small image patches. By using a Markov random field approach we can model images of arbitrary size, while still being able to estimate the model parameters from the data.

ORP1L, a new Rab7 effector and regulator of late endosome functions

Oxysterol binding protein (OSBP) homologues have been found in eukaryotic organisms ranging from yeast to humans. These evolutionary conserved proteins have in common the presence of an OSBP-related domain (ORD) which contains the fully conserved EQVSHHPP sequence motif. The ORD forms a barrel structure that binds sterols in its interior. Other domains and sequence elements found in OSBP-homologues include pleckstrin homology domains, ankyrin repeats and two phenylalanines in an acidic tract (FFAT) motifs, which target the proteins to distinct subcellular compartments. OSBP homologues have been implicated in a wide range of intracellular processes, including vesicle trafficking, lipid metabolism and cell signaling, but little is known about the functional mechanisms of these proteins. The human family of OSBP homologues consists of twelve OSBP-related proteins (ORP). This thesis work is focused on one of the family members, ORP1, of which two variants were found to be expressed tissue-specifically in humans. The shorter variant, ORP1S contains an ORD only. The N-terminally extended variant, ORP1L, comprises a pleckstrin homology domain and three ankyrin repeats in addition to the ORD. The two ORP1 variants differ in intracellular localization. ORP1S is cytosolic, while the ankyrin repeat region of ORP1L targets the protein to late endosomes/lysosomes. This part of ORP1L also has profound effects on late endosomal morphology, inducing perinuclear clustering of late endosomes. A central aim of this study was to identify molecular interactions of ORP1L on late endosomes. The morphological changes of late endosomes induced by overexpressed ORP1L implies involvement of small Rab GTPases, regulators of organelle motility, tethering, docking and/or fusion, in generation of the phenotype. A direct interaction was demonstrated between ORP1L and active Rab7. ORP1L prolongs the active state of Rab7 by stabilizing its GTP-bound form. The clustering of late endosomes/lysosomes was also shown to be linked to the minus end-directed microtubule-based dynein-dynactin motor complex through the ankyrin repeat region of ORP1L. ORP1L, Rab7 and the Rab7-interacting lysosomal protein (RILP) were found to be part of the same effector complex recruiting the dynein-dynactin complex to late endosomes, thereby promoting minus end-directed movement. The proteins were found to be physically close to each other on late endosomes and RILP was found to stabilize the ORP1L-Rab7 interaction. It is possible that ORP1L and RILP bind to each other through their C-terminal and N-terminal regions, respectively, when they are bridged by Rab7. With the results of this study we have been able to place a member of the uncharacterized OSBP-family, ORP1L, in the endocytic pathway, where it regulates motility and possibly fusion of late endosomes through interaction with the small GTPase Rab7.

The biological functions of mouse twinfilin isoforms

The actin cytoskeleton is essential for many cellular processes, including motility, morphogenesis, endocytosis and signal transduction. Actin can exist in monomeric (G-actin) or filamentous (F-actin) form. Actin filaments are considered to be the functional form of actin, generating the protrusive forces characteristic for the actin cytoskeleton. The structure and dynamics of the actin filament and monomer pools are regulated by a large number of actin-binding proteins in eukaryotic cells. Twinfilin is an evolutionarily conserved small actin monomer binding protein. Twinfilin is composed of two ADF/cofilin-like domains, separated by a short linker and followed by a C-terminal tail. Twinfilin forms a stable, high affinity complex with ADP-G-actin, inhibits the nucleotide exchange on actin monomers, and prevents their assembly into filament ends. Twinfilin was originally identified from yeast and has since then been found from all organisms studied except plants. Not much was known about the role of twinfilin in the actin dynamics in mammalian cells before this study. We set out to unravel the mysteries still covering twinfilins functions using biochemistry, cell biology, and genetics. We identified and characterized two mouse isoforms for the previously identified mouse twinfilin-1. The new isoforms, twinfilin-2a and -2b, are generated from the same gene through alternative promoter usage. The three isoforms have distinctive expression patterns, but are similar biochemically. Twinfilin-1 is the major isoform during development and is expressed in high levels in almost all tissues examined. Twinfilin-2a is also expressed almost ubiquitously, but at lower levels. Twinfilin-2b turned out to be a muscle-specific isoform, with very high expression in heart and skeletal muscle. It seems all mouse tissues express at least two twinfilin isoforms, indicating that twinfilins are important regulators of actin dynamics in all cell and tissue types. A knockout mouse line was generated for twinfilin-2a. The mice homozygous for this knockout were viable and developed normally, indicating that twinfilin-2a is dispensable for mouse development. However, it is important to note that twinfilin-2a shows similar expression pattern to twinfilin-1, suggesting that these proteins play redundant roles in mice. All mouse isoforms were shown to be able to sequester actin filaments and have higher affinity for ADP-G-actin than ATP-G-actin. They are also able to directly interact with heterodimeric capping protein and PI(4,5)P2 similar to yeast twinfilin. In this study we also uncovered a novel function for mouse twinfilins; capping actin filament barbed ends. All mouse twinfilin isoforms were shown to possess this function, while yeast and Drosophila twinfilin were not able to cap filament barbed ends. Twinfilins localize to the cytoplasm but also to actin-rich regions in mammalian cells. The subcellular localizations of the isoforms are regulated differently, indicating that even though twinfilins biochemical functions in vitro are very similar, in vivo they can play different roles through different regulatory pathways. Together, this study show that twinfilins regulate actin filament assembly both by sequestering actin monomers and by capping filament barbed ends, and that mammals have three biochemically similar twinfilin isoforms with partially overlapping expression patterns.

Neuromagnetic studies on cortical somatosensory functions in infants and children : Normal development and effect of early brain lesions

Until recently, objective investigation of the functional development of the human brain in vivo was challenged by the lack of noninvasive research methods. Consequently, fairly little is known about cortical processing of sensory information even in healthy infants and children. Furthermore, mechanisms by which early brain insults affect brain development and function are poorly understood. In this thesis, we used magnetoencephalography (MEG) to investigate development of cortical somatosensory functions in healthy infants, very premature infants at risk for neurological disorders, and adolescents with hemiplegic cerebral palsy (CP). In newborns, stimulation of the hand activated both the contralateral primary (SIc) and secondary somatosensory cortices (SIIc). The activation patterns differed from those of adults, however. Some of the earliest SIc responses, constantly present in adults, were completely lacking in newborns and the effect of sleep stage on SIIc responses differed. These discrepancies between newborns and adults reflect the still developmental stage of the newborns’ somatosensory system. Its further maturation was demonstrated by a systematic transformation of the SIc response pattern with age. The main early adult­like components were present by age two. In very preterm infants, at term age, the SIc and SIIc were activated at similar latencies as in healthy fullterm newborns, but the SIc activity was weaker in the preterm group. The SIIc response was absent in four out of the six infants with brain lesions of the underlying hemisphere. Determining the prognostic value of this finding remains a subject for future studies, however. In the CP adolescents with pure subcortical lesions, contrasting their unilateral symptoms, the SIc responses of both hemispheres differed from those of controls: For example the distance between SIc representation areas for digits II and V was shorter bilaterally. In four of the five CP patients with cortico­subcortical brain lesions, no normal early SIc responses were evoked by stimulation of the palsied hand. The varying differences in neuronal functions, underlying the common clinical symptoms, call for investigation of more precisely designed rehabilitation strategies resting on knowledge about individual functional alterations in the sensorimotor networks.

Characterization of LRRTM and NGR gene families: Expression and functions

Some leucine-rich repeat (LRR) -containing membrane proteins are known regulators of neuronal growth and synapse formation. In this work I characterize two gene families encoding neuronal LRR membrane proteins, namely the LRRTM (leucine-rich repeat, transmembrane neuronal) and NGR (Nogo-66 receptor) families. I studied LRRTM and NGR family member's mRNA tissue distribution by RT-PCR and by in situ hybridization. Subcellular localization of LRRTM1 protein was studied in neurons and in non-neuronal cells. I discovered that LRRTM and NGR family mRNAs are predominantly expressed in the nervous system, and that each gene possesses a specific expression pattern. I also established that LRRTM and NGR family mRNAs are expressed by neurons, and not by glial cells. Within neurons, LRRTM1 protein is not transported to the plasma membrane; rather it localizes to endoplasmic reticulum. Nogo-A (RTN4), MAG, and OMgp are myelin-associated proteins that bind to NgR1 to limit axonal regeneration after central nervous system injury. To better understand the functions of NgR2 and NgR3, and to explore the possible redundancy in the signaling of myelin inhibitors of neurite growth, I mapped the interactions between NgR family and the known and candidate NgR1 ligands. I identified high-affinity interactions between RTN2-66, RTN3-66 and NgR1. I also demonstrate that Rtn3 mRNA is expressed in the same glial cell population of mouse spinal cord white matter as Nogo-A mRNA, and thus it could have a role in myelin inhibition of axonal growth. To understand how NgR1 interacts with multiple structurally divergent ligands, I aimed first to map in more detail the nature of Nogo-A:NgR1 interactions, and then to systematically map the binding sites of multiple myelin ligands in NgR1 by using a library of NgR1 expression constructs encoding proteins with one or multiple surface residues mutated to alanine. My analysis of the Nogo-A:NgR1 -interactions revealed a novel interaction site between the proteins, suggesting a trivalent Nogo-A:NgR1-interaction. Our analysis also defined a central binding region on the concave side of NgR1's LRR domain that is required for the binding of all known ligands, and a surrounding region critical for binding MAG and OMgp. To better understand the biological role of LRRTMs, I generated Lrrtm1 and Lrrtm3 knock out mice. I show here that reporter genes expressed from the targeted loci can be used for maping the neuronal connections of Lrrtm1 and Lrrtm3 expressing neurons in finer detail. With regard to LRRTM1's role in humans, we found a strong association between a 70 kb-spanning haplotype in the proposed promoter region of LRRTM1 gene and two possibly related phenotypes: left-handedness and schizophrenia. Interestingly, the responsible haplotype was linked to phenotypic variability only when paternally inherited. In summary, I identified two families of neuronal receptor-like proteins, and mapped their expression and certain protein-protein interactions. The identification of a central binding region in NgR1 shared by multiple ligands may facilitate the design and development of small molecule therapeutics blocking binding of all NgR1 ligands. Additionally, the genetic association data suggests that allelic variation upstream of LRRTM1 may play a role in the development of left-right brain asymmetry in humans. Lrrtm1 and Lrrtm3 knock out mice developed as a part of this study will likely be useful for schizophrenia and Alzheimer s disease research.

Augustine and the functions of concupiscence

This study analyses Augustine s concept of concupiscentia, or evil desire (together with two cognate terms, libido and cupiditas) in the context of his entire oeuvre. By the aid of systematic analysis, the concept and its development is explored in four distinct ways. It is claimed that Augustine used the concept of concupiscentia for several theological purposes, and the task of the study is to represent these distinct functions, and their connections to Augustine s general theological and philosophical convictions. The study opens with a survey on terminology. A general overview of the occurrences of the negatively connoted words for desire in Latin literature precedes a corresponding examination of Augustine s own works. In this introductory chapter it is shown that, despite certain preferences in the uses of the words, a sufficient degree of synonymy reigns so as to allow an analysis of the concept without tightly discriminating between the terms. The theological functions of concupiscentia with its distinct contexts are analysed in separate chapters. The function of concupiscentia as a divine punishment is explored first (Ch 3). It is seen how Augustine links together concupiscentia and ideas about divine justice, and finally suggests that in the inordinate, psychologically experienced sexual desire, the original theological disobedience of Adam and Eve can be perceived. Augustine was criticized for this solution already in his own times, and the analysis of the function of concupiscentia as a divine punishment ends in a discussion on the critical response of punitive concupiscentia by Julian of Aeclanum. Augustine also attached to concupiscentia another central theological function by viewing evil desire as an inward originating cause for all external evil actions. In the study, this function is analysed by surveying two formally distinct images of evil desire, i.e. as the root (radix) of all evil, and as a threefold (triplex) matrix of evil actions (Ch 4). Both of these images were based on a single verse of the Bible (1 Jn 2, 16 and 1 Tim 6, 10). This function of concupiscentia was formed both parallel to, and in answer to, Manichaean insights into concupiscentia. Being familiar with the traditional philosophical discussions on the nature and therapy of emotions, Augustine situated concupiscentia also into this context. It is acknowledged that these philosophical traditions had an obvious impact into his way of explaining psychological processes in connection with concupiscentia. Not only did Augustine implicitly receive and exploit these traditions, but he also explicitly moulded and criticized them in connection with concupiscentia. Eventually, Augustine conceives the philosophical traditions of emotions as partly useful but also partly inadequate to deal with concupiscentia (Ch 5). The role of concupiscentia in connection to divine grace and Christian renewal is analysed in the final chapter of the study. Augustine s gradual development in internalizing the effects of concupiscentia also into the life of a baptized Christian are elucidated, as are the strong limitations and mitigations Augustine makes to the concept when attaching it into the life under grace (sub gratia). A crucial part in the development of this function is played by Augustine s changing interpretation of Rom 7, and the way concupiscentia appears in Augustine s readings of this text is therefore also analysed. As a result of the analysis of these four distinct functions and contexts of concupiscentia, it is concluded that Augustine s concept of concupiscentia is fairly tightly and coherently connected to his views of central theological importance. Especially the functions of concupiscentia as a punishment and the function of concupiscentia in Christian renewal were both tightly interwoven into Augustine s view of God s being and God s grace. The study shows the importance of reading Augustine s discussions on evil desire with a constant awareness of their role in their larger context, that is, of their function in each situation. The study warns against too simplistic and unifying readings of Augustine s concupiscentia, emphasizing the need to acknowledge both the necessitating, sinful aspects of concupiscentia, and the domesticated features of concupiscentia during Christian renewal.

The Structure and Functions of Hantavirus Nucleocapsid Protein

Hantaviruses, members of the genus Hantavirus in the Bunyaviridae family, are enveloped single-stranded RNA viruses with tri-segmented genome of negative polarity. In humans, hantaviruses cause two diseases, hemorrhagic fever with renal syndrome (HFRS) and hantavirus pulmonary syndrome (HPS), which vary in severity depending on the causative agent. Each hantavirus is carried by a specific rodent host and is transmitted to humans through excreta of infected rodents. The genome of hantaviruses encodes four structural proteins: the nucleocapsid protein (N), the glycoproteins (Gn and Gc), and the polymerase (L) and also the nonstructural protein (NSs). This thesis deals with the functional characterization of hantavirus N protein with regard to its structure. Structural studies of the N protein have progressed slowly and the crystal structure of the whole protein is still not available, therefore biochemical assays coupled with bioinformatical modeling proved essential for studying N protein structure and functions. Presumably, during RNA encapsidation, the N protein first forms intermediate trimers and then oligomers. First, we investigated the role of N-terminal domain in the N protein oligomerization. The results suggested that the N-terminal region of the N protein forms a coiled-coil, in which two antiparallel alpha helices interact via their hydrophobic seams. Hydrophobic residues L4, I11, L18, L25 and V32 in the first helix and L44, V51, L58 and L65 in the second helix were crucial for stabilizing the structure. The results were consistent with the head-to-head, tail-to-tail model for hantavirus N protein trimerization. We demonstrated that an intact coiled-coil structure of the N terminus is crucial for the oligomerization capacity of the N protein. We also added new details to the head-to-head, tail-to-tail model of trimerization by suggesting that the initial step is based on interaction(s) between intact intra-molecular coiled-coils of the monomers. We further analyzed the importance of charged aa residues located within the coiled-coil for the N protein oligomerization. To predict the interacting surfaces of the monomers we used an upgraded in silico model of the coiled-coil domain that was docked into a trimer. Next the predicted target residues were mutated. The results obtained using the mammalian two-hybrid assay suggested that conserved charged aa residues within the coiled-coil make a substantial contribution to the N protein oligomerization. This contribution probably involves the formation of interacting surfaces of the N monomers and also stabilization of the coiled-coil via intramolecular ionic bridging. We proposed that the tips of the coiled-coils are the first to come into direct contact and thus initiate tight packing of the three monomers into a compact structure. This was in agreement with the previous results showing that an increase in ionic strength abolished the interaction between N protein molecules. We also showed that residues having the strongest effect on the N protein oligomerization are not scattered randomly throughout the coiled-coil 3D model structure, but form clusters. Next we found evidence for the hantaviral N protein interaction with the cytoplasmic tail of the glycoprotein Gn. In order to study this interaction we used the GST pull-down assay in combination with mutagenesis technique. The results demonstrated that intact, properly folded zinc fingers of the Gn protein cytoplasmic tail as well as the middle domain of the N protein (that includes aa residues 80 248 and supposedly carries the RNA-binding domain) are essential for the interaction. Since hantaviruses do not have a matrix protein that mediates the packaging of the viral RNA in other negatve stranded viruses (NSRV), hantaviral RNPs should be involved in a direct interaction with the intraviral domains of the envelope-embedded glycoproteins. By showing the N-Gn interaction we provided the evidence for one of the crucial steps in the virus replication at which RNPs are directed to the site of the virus assembly. Finally we started analysis of the N protein RNA-binding region, which is supposedly located in the middle domain of the N protein molecule. We developed a model for the initial step of RNA-binding by the hantaviral N protein. We hypothesized that the hantaviral N protein possesses two secondary structure elements that initiate the RNA encapsidation. The results suggest that amino acid residues (172-176) presumably act as a hook to catch vRNA and that the positively charged interaction surface (aa residues 144-160) enhances the initial N-RNA interacation. In conclusion, we elucidated new functions of hantavirus N protein. Using in silico modeling we predicted the domain structure of the protein and using experimental techniques showed that each domain is responsible for executing certain function(s). We showed that intact N terminal coiled-coil domain is crucial for oligomerization and charged residues located on its surface form a interaction surface for the N monomers. The middle domain is essential for interaction with the cytoplasmic tail of the Gn protein and RNA binding.