Computational Methods for Reconstruction and Analysis of Genome-Scale Metabolic Networks

Metabolism is the cellular subsystem responsible for generation of energy from nutrients and production of building blocks for larger macromolecules. Computational and statistical modeling of metabolism is vital to many disciplines including bioengineering, the study of diseases, drug target identification, and understanding the evolution of metabolism. In this thesis, we propose efficient computational methods for metabolic modeling. The techniques presented are targeted particularly at the analysis of large metabolic models encompassing the whole metabolism of one or several organisms. We concentrate on three major themes of metabolic modeling: metabolic pathway analysis, metabolic reconstruction and the study of evolution of metabolism. In the first part of this thesis, we study metabolic pathway analysis. We propose a novel modeling framework called gapless modeling to study biochemically viable metabolic networks and pathways. In addition, we investigate the utilization of atom-level information on metabolism to improve the quality of pathway analyses. We describe efficient algorithms for discovering both gapless and atom-level metabolic pathways, and conduct experiments with large-scale metabolic networks. The presented gapless approach offers a compromise in terms of complexity and feasibility between the previous graph-theoretic and stoichiometric approaches to metabolic modeling. Gapless pathway analysis shows that microbial metabolic networks are not as robust to random damage as suggested by previous studies. Furthermore the amino acid biosynthesis pathways of the fungal species Trichoderma reesei discovered from atom-level data are shown to closely correspond to those of Saccharomyces cerevisiae. In the second part, we propose computational methods for metabolic reconstruction in the gapless modeling framework. We study the task of reconstructing a metabolic network that does not suffer from connectivity problems. Such problems often limit the usability of reconstructed models, and typically require a significant amount of manual postprocessing. We formulate gapless metabolic reconstruction as an optimization problem and propose an efficient divide-and-conquer strategy to solve it with real-world instances. We also describe computational techniques for solving problems stemming from ambiguities in metabolite naming. These techniques have been implemented in a web-based sofware ReMatch intended for reconstruction of models for 13C metabolic flux analysis. In the third part, we extend our scope from single to multiple metabolic networks and propose an algorithm for inferring gapless metabolic networks of ancestral species from phylogenetic data. Experimenting with 16 fungal species, we show that the method is able to generate results that are easily interpretable and that provide hypotheses about the evolution of metabolism.

The Networks of Learning in Technological Innovation : The Emergence of Collaboration Across Fields of Expertise

This study addresses four issues concerning technological product innovations. First, the nature of the very early phases or "embryonic stages" of technological innovation is addressed. Second, this study analyzes why and by what means people initiate innovation processes outside the technological community and the field of expertise of the established industry. In other words, this study addresses the initiation of innovation that occurs without the expertise of established organizations, such as technology firms, professional societies and research institutes operating in the technological field under consideration. Third, the significance of interorganizational learning processes for technological innovation is dealt with. Fourth, this consideration is supplemented by considering how network collaboration and learning change when formalized product development work and the commercialization of innovation advance. These issues are addressed through the empirical analysis of the following three product innovations: Benecol margarine, the Nordic Mobile Telephone system (NMT) and the ProWellness Diabetes Management System (PDMS). This study utilizes the theoretical insights of cultural-historical activity theory on the development of human activities and learning. Activity-theoretical conceptualizations are used in the critical assessment and advancement of the concept of networks of learning. This concept was originally proposed by the research group of organizational scientist Walter Powell. A network of learning refers to the interorganizational collaboration that pools resources, ideas and know-how without market-based or hierarchical relations. The concept of an activity system is used in defining the nodes of the networks of learning. Network collaboration and learning are analyzed with regard to the shared object of development work. According to this study, enduring dilemmas and tensions in activity explain the participants' motives for carrying out actions that lead to novel product concepts in the early phases of technological innovation. These actions comprise the initiation of development work outside the relevant fields of expertise and collaboration and learning across fields of expertise in the absence of market-based or hierarchical relations. These networks of learning are fragile and impermanent. This study suggests that the significance of networks of learning across fields of expertise becomes more and more crucial for innovation activities.

Human brain networks of auditory attention and working memory

This thesis examines brain networks involved in auditory attention and auditory working memory using measures of task performance, brain activity, and neuroanatomical connectivity. Auditory orienting and maintenance of attention were compared with visual orienting and maintenance of attention, and top-down controlled attention was compared to bottom-up triggered attention in audition. Moreover, the effects of cognitive load on performance and brain activity were studied using an auditory working memory task. Corbetta and Shulman s (2002) model of visual attention suggests that what is known as the dorsal attention system (intraparietal sulcus/superior parietal lobule, IPS/SPL and frontal eye field, FEF) is involved in the control of top-down controlled attention, whereas what is known as the ventral attention system (temporo-parietal junction, TPJ and areas of the inferior/middle frontal gyrus, IFG/MFG) is involved in bottom-up triggered attention. The present results show that top-down controlled auditory attention also activates IPS/SPL and FEF. Furthermore, in audition, TPJ and IFG/MFG were activated not only by bottom-up triggered attention, but also by top-down controlled attention. In addition, the posterior cerebellum and thalamus were activated by top-down controlled attention shifts and the ventromedial prefrontal cortex (VMPFC) was activated by to-be-ignored, but attention-catching salient changes in auditory input streams. VMPFC may be involved in the evaluation of environmental events causing the bottom-up triggered engagement of attention. Auditory working memory activated a brain network that largely overlapped with the one activated by top-down controlled attention. The present results also provide further evidence of the role of the cerebellum in cognitive processing: During auditory working memory tasks, both activity in the posterior cerebellum (the crus I/II) and reaction speed increased when the cognitive load increased. Based on the present results and earlier theories on the role of the cerebellum in cognitive processing, the function of the posterior cerebellum in cognitive tasks may be related to the optimization of response speed.

Xanthine oxidoreductase in essential hypertension and metabolic syndrome : Experimental studies on rodent models

Hypertension, obesity, dyslipidemia and dysglycemia constitute metabolic syndrome, a major public health concern, which is associated with cardiovascular mortality. High dietary salt (NaCl) is the most important dietary risk factor for elevated blood pressure. The kidney has a major role in salt-sensitive hypertension and is vulnerable to harmful effects of increased blood pressure. Elevated serum urate is a common finding in these disorders. While dysregulation of urate excretion is associated with cardiovascular diseases, present studies aimed to clarify the role of xanthine oxidoreductase (XOR), i.e. xanthine dehydrogenase (XDH) and its post-translational isoform xanthine oxidase (XO), in cardiovascular diseases. XOR yields urate from hypoxanthine and xanthine. Low oxygen levels upregulate XOR in addition to other factors. In present studies higher renal XOR activity was found in hypertension-prone rats than in the controls. Furthermore, NaCl intake increased renal XOR dose-dependently. To clarify whether XOR has any causal role in hypertension, rats were kept on NaCl diets for different periods of time, with or without a XOR inhibitor, allopurinol. While allopurinol did not alleviate hypertension, it prevented left ventricular and renal hypertrophy. Nitric oxide synthases (NOS) produce nitric oxide (NO), which mediates vasodilatation. A paucity of NO, produced by NOS inhibition, aggravated hypertension and induced renal XOR, whereas NO generating drug, alleviated salt-induced hypertension without changes in renal XOR. Zucker fa/fa rat is an animal model of metabolic syndrome. These rats developed substantial obesity and modest hypertension and showed increased hepatic and renal XOR activities. XOR was modified by diet and antihypertensive treatment. Cyclosporine (CsA) is a fungal peptide and one of the first-line immunosuppressive drugs used in the management of organ transplantation. Nephrotoxicity ensue high doses resulting in hypertension and limit CsA use. CsA increased renal XO substantially in salt-sensitive rats on a high NaCl diet, indicating a possible role for this reactive oxygen species generating isoform in CsA nephrotoxicity. Renal hypoxia, common to these rodent models of hypertension and obesity, is one of the plausible XOR inducing factors. Although XOR inhibition did not prevent hypertension, present experimental data indicate that XOR plays a role in the pathology of salt-induced cardiac and renal hypertrophy.

Serum uric acid and metabolic risk factors in three ethnic groups: Asian Indians and Creoles in Mauritius and Chinese in Qingdao, China

In humans with a loss of uricase the final oxidation product of purine catabolism is uric acid (UA). The prevalence of hyperuricemia has been increasing around the world accompanied by a rapid increase in obesity and diabetes. Since hyperuricemia was first described as being associated with hyperglycemia and hypertension by Kylin in 1923, there has been a growing interest in the association between elevated UA and other metabolic abnormalities of hyperglycemia, abdominal obesity, dyslipidemia, and hypertension. The direction of causality between hyperuricemia and metabolic disorders, however, is unceartain. The association of UA with metabolic abnormalities still needs to be delineated in population samples. Our overall aims were to study the prevalence of hyperuricemia and the metabolic factors clustering with hyperuricemia, to explore the dynamical changes in blood UA levels with the deterioration in glucose metabolism and to estimate the predictive capability of UA in the development of diabetes. Four population-based surveys for diabetes and other non-communicable diseases were conducted in 1987, 1992, and 1998 in Mauritius, and in 2001-2002 in Qingdao, China. The Qingdao study comprised 1 288 Chinese men and 2 344 women between 20-74, and the Mauritius study consisted of 3 784 Mauritian Indian and Mauritian Creole men and 4 442 women between 25-74. In Mauritius, re-exams were made in 1992 and/or 1998 for 1 941 men (1 409 Indians and 532 Creoles) and 2 318 non pregnant women (1 645 Indians and 673 Creoles), free of diabetes, cardiovascular diseases, and gout at baseline examinations in 1987 or 1992, using the same study protocol. The questionnaire was designed to collect demographic details, physical examinations and standard 75g oral glucose tolerance tests were performed in all cohorts. Fasting blood UA and lipid profiles were also determined. The age-standardized prevalence in Chinese living in Qingdao was 25.3% for hyperuricemia (defined as fasting serum UA > 420 μmol/l in men and > 360 μmol/l in women) and 0.36% for gout in adults between 20-74. Hyperuricemia was more prevalent in men than in women. One standard deviation increase in UA concentration was associated with the clustering of metabolic risk factors for both men and women in three ethnic groups. Waist circumference, body mass index, and serum triglycerides appeared to be independently associated with hyperuricemia in both sexes and in all ethnic groups except in Chinese women, in whom triglycerides, high-density lipoprotein cholesterol, and total cholesterol were associated with hyperuricemia. Serum UA increased with increasing fasting plasma glucose levels up to a value of 7.0 mmol/l, but significantly decreased thereafter in mainland Chinese. An inverse relationship occurred between 2-h plasma glucose and serum UA when 2-h plasma glucose higher than 8.0 mmol/l. In the prospective study in Mauritius, 337 (17.4%) men and 379 (16.4%) women developed diabetes during the follow-up. Elevated UA levels at baseline increased 1.14-fold in risk of incident diabetes in Indian men and 1.37-fold in Creole men, but no significant risk was observed in women. In conclusion, the prevalence of hyperuricemia was high in Chinese in Qingdao, blood UA was associated with the clustering of metabolic risk factors in Mauritian Indian, Mauritian Creole, and Chinese living in Qingdao, and a high baseline UA level independently predicted the development of diabetes in Mauritian men. The clinical use of UA as a marker of hyperglycemia and other metabolic disorders needs to be further studied. Keywords: Uric acid, Hyperuricemia, Risk factors, Type 2 Diabetes, Incidence, Mauritius, Chinese

IP Mobility in Wireless Operator Networks

Wireless network access is gaining increased heterogeneity in terms of the types of IP capable access technologies. The access network heterogeneity is an outcome of incremental and evolutionary approach of building new infrastructure. The recent success of multi-radio terminals drives both building a new infrastructure and implicit deployment of heterogeneous access networks. Typically there is no economical reason to replace the existing infrastructure when building a new one. The gradual migration phase usually takes several years. IP-based mobility across different access networks may involve both horizontal and vertical handovers. Depending on the networking environment, the mobile terminal may be attached to the network through multiple access technologies. Consequently, the terminal may send and receive packets through multiple networks simultaneously. This dissertation addresses the introduction of IP Mobility paradigm into the existing mobile operator network infrastructure that have not originally been designed for multi-access and IP Mobility. We propose a model for the future wireless networking and roaming architecture that does not require revolutionary technology changes and can be deployed without unnecessary complexity. The model proposes a clear separation of operator roles: (i) access operator, (ii) service operator, and (iii) inter-connection and roaming provider. The separation allows each type of an operator to have their own development path and business models without artificial bindings with each other. We also propose minimum requirements for the new model. We present the state of the art of IP Mobility. We also present results of standardization efforts in IP-based wireless architectures. Finally, we present experimentation results of IP-level mobility in various wireless operator deployments.

A Decentralized Session Management Framework for Heterogeneous Ad-Hoc and Fixed Networks

Wireless technologies are continuously evolving. Second generation cellular networks have gained worldwide acceptance. Wireless LANs are commonly deployed in corporations or university campuses, and their diffusion in public hotspots is growing. Third generation cellular systems are yet to affirm everywhere; still, there is an impressive amount of research ongoing for deploying beyond 3G systems. These new wireless technologies combine the characteristics of WLAN based and cellular networks to provide increased bandwidth. The common direction where all the efforts in wireless technologies are headed is towards an IP-based communication. Telephony services have been the killer application for cellular systems; their evolution to packet-switched networks is a natural path. Effective IP telephony signaling protocols, such as the Session Initiation Protocol (SIP) and the H 323 protocol are needed to establish IP-based telephony sessions. However, IP telephony is just one service example of IP-based communication. IP-based multimedia sessions are expected to become popular and offer a wider range of communication capabilities than pure telephony. In order to conjoin the advances of the future wireless technologies with the potential of IP-based multimedia communication, the next step would be to obtain ubiquitous communication capabilities. According to this vision, people must be able to communicate also when no support from an infrastructured network is available, needed or desired. In order to achieve ubiquitous communication, end devices must integrate all the capabilities necessary for IP-based distributed and decentralized communication. Such capabilities are currently missing. For example, it is not possible to utilize native IP telephony signaling protocols in a totally decentralized way. This dissertation presents a solution for deploying the SIP protocol in a decentralized fashion without support of infrastructure servers. The proposed solution is mainly designed to fit the needs of decentralized mobile environments, and can be applied to small scale ad-hoc networks or also bigger networks with hundreds of nodes. A framework allowing discovery of SIP users in ad-hoc networks and the establishment of SIP sessions among them, in a fully distributed and secure way, is described and evaluated. Security support allows ad-hoc users to authenticate the sender of a message, and to verify the integrity of a received message. The distributed session management framework has been extended in order to achieve interoperability with the Internet, and the native Internet applications. With limited extensions to the SIP protocol, we have designed and experimentally validated a SIP gateway allowing SIP signaling between ad-hoc networks with private addressing space and native SIP applications in the Internet. The design is completed by an application level relay that permits instant messaging sessions to be established in heterogeneous environments. The resulting framework constitutes a flexible and effective approach for the pervasive deployment of real time applications.

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.

TCP Performance in Heterogeneous Wireless Networks

The TCP protocol is used by most Internet applications today, including the recent mobile wireless terminals that use TCP for their World-Wide Web, E-mail and other traffic. The recent wireless network technologies, such as GPRS, are known to cause delay spikes in packet transfer. This causes unnecessary TCP retransmission timeouts. This dissertation proposes a mechanism, Forward RTO-Recovery (F-RTO) for detecting the unnecessary TCP retransmission timeouts and thus allow TCP to take appropriate follow-up actions. We analyze a Linux F-RTO implementation in various network scenarios and investigate different alternatives to the basic algorithm. The second part of this dissertation is focused on quickly adapting the TCP's transmission rate when the underlying link characteristics change suddenly. This can happen, for example, due to vertical hand-offs between GPRS and WLAN wireless technologies. We investigate the Quick-Start algorithm that, in collaboration with the network routers, aims to quickly probe the available bandwidth on a network path, and allow TCP's congestion control algorithms to use that information. By extensive simulations we study the different router algorithms and parameters for Quick-Start, and discuss the challenges Quick-Start faces in the current Internet. We also study the performance of Quick-Start when applied to vertical hand-offs between different wireless link technologies.