Towards robust network based complex systems : from evolutionary cellular automata to biological models

Abstract Sitting between your past and your future doesn't mean you are in the present. Dakota Skye Complex systems science is an interdisciplinary field grouping under the same umbrella dynamical phenomena from social, natural or mathematical sciences. The emergence of a higher order organization or behavior, transcending that expected of the linear addition of the parts, is a key factor shared by all these systems. Most complex systems can be modeled as networks that represent the interactions amongst the system's components. In addition to the actual nature of the part's interactions, the intrinsic topological structure of underlying network is believed to play a crucial role in the remarkable emergent behaviors exhibited by the systems. Moreover, the topology is also a key a factor to explain the extraordinary flexibility and resilience to perturbations when applied to transmission and diffusion phenomena. In this work, we study the effect of different network structures on the performance and on the fault tolerance of systems in two different contexts. In the first part, we study cellular automata, which are a simple paradigm for distributed computation. Cellular automata are made of basic Boolean computational units, the cells; relying on simple rules and information from- the surrounding cells to perform a global task. The limited visibility of the cells can be modeled as a network, where interactions amongst cells are governed by an underlying structure, usually a regular one. In order to increase the performance of cellular automata, we chose to change its topology. We applied computational principles inspired by Darwinian evolution, called evolutionary algorithms, to alter the system's topological structure starting from either a regular or a random one. The outcome is remarkable, as the resulting topologies find themselves sharing properties of both regular and random network, and display similitudes Watts-Strogtz's small-world network found in social systems. Moreover, the performance and tolerance to probabilistic faults of our small-world like cellular automata surpasses that of regular ones. In the second part, we use the context of biological genetic regulatory networks and, in particular, Kauffman's random Boolean networks model. In some ways, this model is close to cellular automata, although is not expected to perform any task. Instead, it simulates the time-evolution of genetic regulation within living organisms under strict conditions. The original model, though very attractive by it's simplicity, suffered from important shortcomings unveiled by the recent advances in genetics and biology. We propose to use these new discoveries to improve the original model. Firstly, we have used artificial topologies believed to be closer to that of gene regulatory networks. We have also studied actual biological organisms, and used parts of their genetic regulatory networks in our models. Secondly, we have addressed the improbable full synchronicity of the event taking place on. Boolean networks and proposed a more biologically plausible cascading scheme. Finally, we tackled the actual Boolean functions of the model, i.e. the specifics of how genes activate according to the activity of upstream genes, and presented a new update function that takes into account the actual promoting and repressing effects of one gene on another. Our improved models demonstrate the expected, biologically sound, behavior of previous GRN model, yet with superior resistance to perturbations. We believe they are one step closer to the biological reality.

Multi-scale systems analysis of plant development : beyond cellular signaling and tissue morphodynamics

Les plantes sont essentielles pour les sociétés humaines. Notre alimentation quotidienne, les matériaux de constructions et les sources énergétiques dérivent de la biomasse végétale. En revanche, la compréhension des multiples aspects développementaux des plantes est encore peu exploitée et représente un sujet de recherche majeur pour la science. L'émergence des technologies à haut débit pour le séquençage de génome à grande échelle ou l'imagerie de haute résolution permet à présent de produire des quantités énormes d'information. L'analyse informatique est une façon d'intégrer ces données et de réduire la complexité apparente vers une échelle d'abstraction appropriée, dont la finalité est de fournir des perspectives de recherches ciblées. Ceci représente la raison première de cette thèse. En d'autres termes, nous appliquons des méthodes descriptives et prédictives combinées à des simulations numériques afin d'apporter des solutions originales à des problèmes relatifs à la morphogénèse à l'échelle de la cellule et de l'organe. Nous nous sommes fixés parmi les objectifs principaux de cette thèse d'élucider de quelle manière l'interaction croisée des phytohormones auxine et brassinosteroïdes (BRs) détermine la croissance de la cellule dans la racine du méristème apical d'Arabidopsis thaliana, l'organisme modèle de référence pour les études moléculaires en plantes. Pour reconstruire le réseau de signalement cellulaire, nous avons extrait de la littérature les informations pertinentes concernant les relations entre les protéines impliquées dans la transduction des signaux hormonaux. Le réseau a ensuite été modélisé en utilisant un formalisme logique et qualitatif pour pallier l'absence de données quantitatives. Tout d'abord, Les résultats ont permis de confirmer que l'auxine et les BRs agissent en synergie pour contrôler la croissance de la cellule, puis, d'expliquer des observations phénotypiques paradoxales et au final, de mettre à jour une interaction clef entre deux protéines dans la maintenance du méristème de la racine. Une étude ultérieure chez la plante modèle Brachypodium dystachion (Brachypo- dium) a révélé l'ajustement du réseau d'interaction croisée entre auxine et éthylène par rapport à Arabidopsis. Chez ce dernier, interférer avec la biosynthèse de l'auxine mène à la formation d'une racine courte. Néanmoins, nous avons isolé chez Brachypodium un mutant hypomorphique dans la biosynthèse de l'auxine qui affiche une racine plus longue. Nous avons alors conduit une analyse morphométrique qui a confirmé que des cellules plus anisotropique (plus fines et longues) sont à l'origine de ce phénotype racinaire. Des analyses plus approfondies ont démontré que la différence phénotypique entre Brachypodium et Arabidopsis s'explique par une inversion de la fonction régulatrice dans la relation entre le réseau de signalisation par l'éthylène et la biosynthèse de l'auxine. L'analyse morphométrique utilisée dans l'étude précédente exploite le pipeline de traitement d'image de notre méthode d'histologie quantitative. Pendant la croissance secondaire, la symétrie bilatérale de l'hypocotyle est remplacée par une symétrie radiale et une organisation concentrique des tissus constitutifs. Ces tissus sont initialement composés d'une douzaine de cellules mais peuvent aisément atteindre des dizaines de milliers dans les derniers stades du développement. Cette échelle dépasse largement le seuil d'investigation par les moyens dits 'traditionnels' comme l'imagerie directe de tissus en profondeur. L'étude de ce système pendant cette phase de développement ne peut se faire qu'en réalisant des coupes fines de l'organe, ce qui empêche une compréhension des phénomènes cellulaires dynamiques sous-jacents. Nous y avons remédié en proposant une stratégie originale nommée, histologie quantitative. De fait, nous avons extrait l'information contenue dans des images de très haute résolution de sections transverses d'hypocotyles en utilisant un pipeline d'analyse et de segmentation d'image à grande échelle. Nous l'avons ensuite combiné avec un algorithme de reconnaissance automatique des cellules. Cet outil nous a permis de réaliser une description quantitative de la progression de la croissance secondaire révélant des schémas développementales non-apparents avec une inspection visuelle classique. La formation de pôle de phloèmes en structure répétée et espacée entre eux d'une longueur constante illustre les bénéfices de notre approche. Par ailleurs, l'exploitation approfondie de ces résultats a montré un changement de croissance anisotropique des cellules du cambium et du phloème qui semble en phase avec l'expansion du xylème. Combinant des outils génétiques et de la modélisation biomécanique, nous avons démontré que seule la croissance plus rapide des tissus internes peut produire une réorientation de l'axe de croissance anisotropique des tissus périphériques. Cette prédiction a été confirmée par le calcul du ratio des taux de croissance du xylème et du phloème au cours de développement secondaire ; des ratios élevés sont effectivement observés et concomitant à l'établissement progressif et tangentiel du cambium. Ces résultats suggèrent un mécanisme d'auto-organisation établi par un gradient de division méristématique qui génèrent une distribution de contraintes mécaniques. Ceci réoriente la croissance anisotropique des tissus périphériques pour supporter la croissance secondaire. - Plants are essential for human society, because our daily food, construction materials and sustainable energy are derived from plant biomass. Yet, despite this importance, the multiple developmental aspects of plants are still poorly understood and represent a major challenge for science. With the emergence of high throughput devices for genome sequencing and high-resolution imaging, data has never been so easy to collect, generating huge amounts of information. Computational analysis is one way to integrate those data and to decrease the apparent complexity towards an appropriate scale of abstraction with the aim to eventually provide new answers and direct further research perspectives. This is the motivation behind this thesis work, i.e. the application of descriptive and predictive analytics combined with computational modeling to answer problems that revolve around morphogenesis at the subcellular and organ scale. One of the goals of this thesis is to elucidate how the auxin-brassinosteroid phytohormone interaction determines the cell growth in the root apical meristem of Arabidopsis thaliana (Arabidopsis), the plant model of reference for molecular studies. The pertinent information about signaling protein relationships was obtained through the literature to reconstruct the entire hormonal crosstalk. Due to a lack of quantitative information, we employed a qualitative modeling formalism. This work permitted to confirm the synergistic effect of the hormonal crosstalk on cell elongation, to explain some of our paradoxical mutant phenotypes and to predict a novel interaction between the BREVIS RADIX (BRX) protein and the transcription factor MONOPTEROS (MP),which turned out to be critical for the maintenance of the root meristem. On the same subcellular scale, another study in the monocot model Brachypodium dystachion (Brachypodium) revealed an alternative wiring of auxin-ethylene crosstalk as compared to Arabidopsis. In the latter, increasing interference with auxin biosynthesis results in progressively shorter roots. By contrast, a hypomorphic Brachypodium mutant isolated in this study in an enzyme of the auxin biosynthesis pathway displayed a dramatically longer seminal root. Our morphometric analysis confirmed that more anisotropic cells (thinner and longer) are principally responsible for the mutant root phenotype. Further characterization pointed towards an inverted regulatory logic in the relation between ethylene signaling and auxin biosynthesis in Brachypodium as compared to Arabidopsis, which explains the phenotypic discrepancy. Finally, the morphometric analysis of hypocotyl secondary growth that we applied in this study was performed with the image-processing pipeline of our quantitative histology method. During its secondary growth, the hypocotyl reorganizes its primary bilateral symmetry to a radial symmetry of highly specialized tissues comprising several thousand cells, starting with a few dozens. However, such a scale only permits observations in thin cross-sections, severely hampering a comprehensive analysis of the morphodynamics involved. Our quantitative histology strategy overcomes this limitation. We acquired hypocotyl cross-sections from tiled high-resolution images and extracted their information content using custom high-throughput image processing and segmentation. Coupled with an automated cell type recognition algorithm, it allows precise quantitative characterization of vascular development and reveals developmental patterns that were not evident from visual inspection, for example the steady interspace distance of the phloem poles. Further analyses indicated a change in growth anisotropy of cambial and phloem cells, which appeared in phase with the expansion of xylem. Combining genetic tools and computational modeling, we showed that the reorientation of growth anisotropy axis of peripheral tissue layers only occurs when the growth rate of central tissue is higher than the peripheral one. This was confirmed by the calculation of the ratio of the growth rate xylem to phloem throughout secondary growth. High ratios are indeed observed and concomitant with the homogenization of cambium anisotropy. These results suggest a self-organization mechanism, promoted by a gradient of division in the cambium that generates a pattern of mechanical constraints. This, in turn, reorients the growth anisotropy of peripheral tissues to sustain the secondary growth.

Comparative study and identification of potent eukaryotic transcriptional repressors in gene switch systems.

In mammalian cells, proper gene regulation is achieved by the complex interplay of transcription factors that activate or repress gene expression by binding to the regulatory regions of target promoters. While transcriptional activators have been extensively characterised and classified into functional groups, relatively little is known about the comparative strength and cell type-specificity of transcriptional repressors. Here, we have compared the ability of a series of eukaryotic repression domains to silence basal and activated transcription. A series of the most potent repression domains was further tested in the context of a gene therapy gene-switch system in various cell types. The results indicate that the analysed repression domains exert varying silencing activities in different promoter contexts. Furthermore, their potential for gene silencing varies also depending on the cellular context. When multimerised within one chimeric repressor protein, particular combinations of repressor domains were found to display synergistic repressing effects and efficient repression in a panel of cell lines. This approach thus allowed the identification of transcriptional repressors that are both potent and versatile in terms of cellular specificity as a basis for gene switch systems.

Cell Therapy Assistance in Reconstructive Surgery for Musculoskeletal Tissues Following Burn and Trauma: Swiss Cellular Transplantation Platform

The sterol compositions of three oceanic jellyfish have been determined using gas chromatographic mass spectrometric techniques involving the use of two separate gas chromatographic column systems. The components in overlapping peaks have been identified by comparison of the mass spectra of peaks in the two column systems using subtractive techniques. A mid-water animal, Periphylla periphylla, was found to contain a very complex and unusual sterol profile including rare 5alpha-stanols, whereas two other oceanic jellyfish Pelagia noctiluca and Atolla wyvillei contained similar mixtures of delta5 sterols to those previously isolated from coastal species.

Interaction of alpha-melanotropin with central dopamine systems: role of hormonal state and molecular structure

The tubero-infundibular and nigrostriatal DA neurone systems of rats respond to systemic (i.p.) injection of alpha-MSH (2-100 microgram/kg). The response of the tubero-infundibular (arcuate) DA neurones, an increase in cellular fluorescence intensity which can be interpreted as a sign of increased neuronal activity, is essentially the same in males, estrogen-progesterone-pretreated ovariectomized females and hypophysectomized males, whereas the type of response elicited by alpha-MSH in the nigral DA neurones depends upon the hormonal state of the animal. Differences between the two DA neurone groups exist also with regard to the effects of peptide fragments containing the two active sites of the alpha-MSH molecule. Results of lesion experiments in the lower brainstem (area postrema) and of blockade of muscarinic mechanisms by atropine further point to differences in the mechanisms underlying the peptide effects on the two neurone systems. The reaction of the tubero-infundibular DA system (which controls the pars intermedia of the pituitary) can be considered to reflect the activation of a feedback mechanism on MSH secretion, while the functional counterpart of the changes observed in the nigral system remains unknown at the present time.

Identification of GPx8 as a novel cellular substrate of the hepatitis C virus NS3-4A protease

Background: The hepatitis C virus (HCV) NS3-4A protease is not only an essential component of the viral replication complex and a prime target for antiviral intervention but also a key player in the persistence and pathogenesis of HCV. It cleaves and thereby inactivates two crucial adaptor proteins in viral RNA sensing and innate immunity (MAVS and TRIF) as well as a phosphatase involved in growth factor signaling (TC-PTP). The aim of this study was to identify novel cellular substrates of the NS3-4A protease and to investigate their role in the life cycle and pathogenesis of HCV. Methods: Cell lines inducibly expressing the NS3-4A protease were analyzed in basal as well as interferon- α -stimulated states by stable isotopic labeling using amino acids in cell culture (SILAC) coupled with protein separation and mass spectrometry. Candidates fulfilling strin- gent criteria for potential substrates or products of the NS3-4A protease were further investigated in different experimental sys- tems as well as in liver biopsies from patients with chronic hep- atitis C. Results: SILAC coupled with protein separation and mass spectrometry yielded > 5000 proteins of which 21 can- didates were selected for further analyses. These allowed us to identify GPx8, a membrane-associated peroxidase involved in disulfide bond formation in the endoplasmic reticulum, as a novel cellular substrate of the HCV NS3-4A protease. Cleavage occurs at cysteine in position 11, removing the cytosolic tip of GPx8, and was observed in different experimental systems as well as in liver biopsies from patients with chronic hepatitis C. Further functional studies, involving overexpression and RNA silencing, revealed that GPx8 is a proviral factor involved in viral particle production but not in HCV entry or RNA replica- tion. Conclusions: GPx8 is a proviral host factor cleaved by the HCV NS3-4A protease. Studies investigating the consequences of cleavage for GPx8 function are underway. The identification of novel cellular substrates of the HCV NS3-4A protease should yield new insights into the HCV life cycle and the pathogenesis of hepatitis C and may reveal novel angles for therapeutic inter- vention.

Identification of GPx8 as a novel cellular substrate of the hepatitis C virus NS3-4A protease

Background: The hepatitis C virus (HCV) NS3-4A protease is not only an essential component of the viral replication complex and a prime target for a ntiviral intervention but also a key player i n the persistence and pathogenesis of HCV. It cleaves and thereby inactivates two crucial adaptor proteins in viral RNA sensing and innate immunity (MAVS and TRIF) as well as a phosphatase involved in growth factor signaling (TCPTP). T he aim of this study was to identify novel cellular substrates o f the N S3-4A protease and to investigate their role in the replication and pathogenesis of HCV. Methods: Cell lines inducibly expressing t he NS3-4A protease were analyzed in basal as well as interferon-α-stimulated states by stable isotopic l abeling using amino acids in cell culture (SILAC) coupled with protein separation and mass spectrometry. Candidates fulfilling stringent criteria for potential substrates or products of the NS3-4A protease were further i nvestigated in different experimental systems as well a s in liver biopsies from patients with chronic hepatitis C. Results: SILAC coupled with protein separation and mass spectrometry yielded > 5000 proteins of which 18 candidates were selected for further analyses. These allowed us to identify GPx8, a membrane-associated peroxidase involved in disulfide bond formation in the endoplasmic reticulum, as a n ovel cellular substrate of the H CV NS3-4A protease. Cleavage occurs at cysteine in position 11, removing the cytosolic tip of GPx8, and was observed in different experimental systems as well as in liver biopsies from patients with chronic hepatitis C. Further functional studies, involving overexpression and RNA silencing, revealed that GPx8 is a p roviral factor involved in viral particle production but not in HCV entry or HCV RNA replication. Conclusions: GPx8 is a proviral host factor cleaved by the HCV NS3-4A protease. Studies investigating the consequences of GPx8 cleavage for protein function are underway. The identification of novel cellular substrates o f the HCV N S3-4A protease should yield new insights i nto the HCV life cycle and the pathogenesis of hepatitis C and may reveal novel targets for antiviral intervention.

Cellular delivery of pyrenyl-arene ruthenium complexes by a water-soluble arene ruthenium metalla-cage.

Three pyrenyl-arene ruthenium complexes (M(1)-M(3)) of the general formula [Ru(η(6)-arene-pyrenyl)Cl(2)(pta)] (pta = 1,3,5-triaza-7-phosphaadamantane) have been synthesised and characterised. Prior to the coordination to ruthenium, pyrene was connected to the arene ligand via an alkane chain containing different functional groups: ester (L(1)), ether (L(2)) and amide (L(3)), respectively. Furthermore, the pyrenyl moieties of the M(n) complexes were encapsulated within the hydrophobic cavity of the water soluble metalla-cage, [Ru(6)(η(6)-p-cymene)(6)(tpt)(2)(donq)(3)](6+) (tpt = 2,4,6-tri-(pyridin-4-yl)-1,3,5-triazine; donq = 5,8-dioxydo-1,4-naphthoquinonato), while the arene ruthenium end was pointing out of the cage, thus giving rise to the corresponding host-guest systems [M(n)⊂Ru(6)(η(6)-p-cymene)(6)(tpt)(2)(donq)(3)](6+) ([M(n)⊂cage](6+)). The antitumor activity of the pyrenyl-arene ruthenium complexes (M(n)) and the corresponding host-guest systems [M(n)⊂cage][CF(3)SO(3)](6) were evaluated in vitro in different types of human cancer cell lines (A549, A2780, A2780cisR, Me300 and HeLa). Complex M(2), which contains an ether group within the alkane chain, demonstrated at least a 10 times higher cytotoxicity than the reference compound [Ru(η(6)-p-cymene)Cl(2)(pta)] (RAPTA-C). All host-guest systems [M(n)⊂cage](6+) showed good anticancer activity with IC(50) values ranging from 2 to 8 μM after 72 h exposure. The fluorescence of the pyrenyl moiety allowed the monitoring of the cellular uptake and revealed an increase of uptake by a factor two of the M(2) complex when encapsulated in the metalla-cage [Ru(6)(η(6)-p-cymene)(6)(tpt)(2)(donq)(3)](6+).

Behavioural queueing : cellular automata and a laboratory experiment

Dans cette thèse, nous étudions les aspects comportementaux d'agents qui interagissent dans des systèmes de files d'attente à l'aide de modèles de simulation et de méthodologies expérimentales. Chaque période les clients doivent choisir un prestataire de servivce. L'objectif est d'analyser l'impact des décisions des clients et des prestataires sur la formation des files d'attente. Dans un premier cas nous considérons des clients ayant un certain degré d'aversion au risque. Sur la base de leur perception de l'attente moyenne et de la variabilité de cette attente, ils forment une estimation de la limite supérieure de l'attente chez chacun des prestataires. Chaque période, ils choisissent le prestataire pour lequel cette estimation est la plus basse. Nos résultats indiquent qu'il n'y a pas de relation monotone entre le degré d'aversion au risque et la performance globale. En effet, une population de clients ayant un degré d'aversion au risque intermédiaire encoure généralement une attente moyenne plus élevée qu'une population d'agents indifférents au risque ou très averses au risque. Ensuite, nous incorporons les décisions des prestataires en leur permettant d'ajuster leur capacité de service sur la base de leur perception de la fréquence moyenne d'arrivées. Les résultats montrent que le comportement des clients et les décisions des prestataires présentent une forte "dépendance au sentier". En outre, nous montrons que les décisions des prestataires font converger l'attente moyenne pondérée vers l'attente de référence du marché. Finalement, une expérience de laboratoire dans laquelle des sujets jouent le rôle de prestataire de service nous a permis de conclure que les délais d'installation et de démantèlement de capacité affectent de manière significative la performance et les décisions des sujets. En particulier, les décisions du prestataire, sont influencées par ses commandes en carnet, sa capacité de service actuellement disponible et les décisions d'ajustement de capacité qu'il a prises, mais pas encore implémentées. - Queuing is a fact of life that we witness daily. We all have had the experience of waiting in line for some reason and we also know that it is an annoying situation. As the adage says "time is money"; this is perhaps the best way of stating what queuing problems mean for customers. Human beings are not very tolerant, but they are even less so when having to wait in line for service. Banks, roads, post offices and restaurants are just some examples where people must wait for service. Studies of queuing phenomena have typically addressed the optimisation of performance measures (e.g. average waiting time, queue length and server utilisation rates) and the analysis of equilibrium solutions. The individual behaviour of the agents involved in queueing systems and their decision making process have received little attention. Although this work has been useful to improve the efficiency of many queueing systems, or to design new processes in social and physical systems, it has only provided us with a limited ability to explain the behaviour observed in many real queues. In this dissertation we differ from this traditional research by analysing how the agents involved in the system make decisions instead of focusing on optimising performance measures or analysing an equilibrium solution. This dissertation builds on and extends the framework proposed by van Ackere and Larsen (2004) and van Ackere et al. (2010). We focus on studying behavioural aspects in queueing systems and incorporate this still underdeveloped framework into the operations management field. In the first chapter of this thesis we provide a general introduction to the area, as well as an overview of the results. In Chapters 2 and 3, we use Cellular Automata (CA) to model service systems where captive interacting customers must decide each period which facility to join for service. They base this decision on their expectations of sojourn times. Each period, customers use new information (their most recent experience and that of their best performing neighbour) to form expectations of sojourn time at the different facilities. Customers update their expectations using an adaptive expectations process to combine their memory and their new information. We label "conservative" those customers who give more weight to their memory than to the xiv Summary new information. In contrast, when they give more weight to new information, we call them "reactive". In Chapter 2, we consider customers with different degree of risk-aversion who take into account uncertainty. They choose which facility to join based on an estimated upper-bound of the sojourn time which they compute using their perceptions of the average sojourn time and the level of uncertainty. We assume the same exogenous service capacity for all facilities, which remains constant throughout. We first analyse the collective behaviour generated by the customers' decisions. We show that the system achieves low weighted average sojourn times when the collective behaviour results in neighbourhoods of customers loyal to a facility and the customers are approximately equally split among all facilities. The lowest weighted average sojourn time is achieved when exactly the same number of customers patronises each facility, implying that they do not wish to switch facility. In this case, the system has achieved the Nash equilibrium. We show that there is a non-monotonic relationship between the degree of risk-aversion and system performance. Customers with an intermediate degree of riskaversion typically achieve higher sojourn times; in particular they rarely achieve the Nash equilibrium. Risk-neutral customers have the highest probability of achieving the Nash Equilibrium. Chapter 3 considers a service system similar to the previous one but with risk-neutral customers, and relaxes the assumption of exogenous service rates. In this sense, we model a queueing system with endogenous service rates by enabling managers to adjust the service capacity of the facilities. We assume that managers do so based on their perceptions of the arrival rates and use the same principle of adaptive expectations to model these perceptions. We consider service systems in which the managers' decisions take time to be implemented. Managers are characterised by a profile which is determined by the speed at which they update their perceptions, the speed at which they take decisions, and how coherent they are when accounting for their previous decisions still to be implemented when taking their next decision. We find that the managers' decisions exhibit a strong path-dependence: owing to the initial conditions of the model, the facilities of managers with identical profiles can evolve completely differently. In some cases the system becomes "locked-in" into a monopoly or duopoly situation. The competition between managers causes the weighted average sojourn time of the system to converge to the exogenous benchmark value which they use to estimate their desired capacity. Concerning the managers' profile, we found that the more conservative Summary xv a manager is regarding new information, the larger the market share his facility achieves. Additionally, the faster he takes decisions, the higher the probability that he achieves a monopoly position. In Chapter 4 we consider a one-server queueing system with non-captive customers. We carry out an experiment aimed at analysing the way human subjects, taking on the role of the manager, take decisions in a laboratory regarding the capacity of a service facility. We adapt the model proposed by van Ackere et al (2010). This model relaxes the assumption of a captive market and allows current customers to decide whether or not to use the facility. Additionally the facility also has potential customers who currently do not patronise it, but might consider doing so in the future. We identify three groups of subjects whose decisions cause similar behavioural patterns. These groups are labelled: gradual investors, lumpy investors, and random investor. Using an autocorrelation analysis of the subjects' decisions, we illustrate that these decisions are positively correlated to the decisions taken one period early. Subsequently we formulate a heuristic to model the decision rule considered by subjects in the laboratory. We found that this decision rule fits very well for those subjects who gradually adjust capacity, but it does not capture the behaviour of the subjects of the other two groups. In Chapter 5 we summarise the results and provide suggestions for further work. Our main contribution is the use of simulation and experimental methodologies to explain the collective behaviour generated by customers' and managers' decisions in queueing systems as well as the analysis of the individual behaviour of these agents. In this way, we differ from the typical literature related to queueing systems which focuses on optimising performance measures and the analysis of equilibrium solutions. Our work can be seen as a first step towards understanding the interaction between customer behaviour and the capacity adjustment process in queueing systems. This framework is still in its early stages and accordingly there is a large potential for further work that spans several research topics. Interesting extensions to this work include incorporating other characteristics of queueing systems which affect the customers' experience (e.g. balking, reneging and jockeying); providing customers and managers with additional information to take their decisions (e.g. service price, quality, customers' profile); analysing different decision rules and studying other characteristics which determine the profile of customers and managers.

Human fetal bone cells in delivery systems for bone engineering.

The aim of this study was to culture human fetal bone cells (dedicated cell banks of fetal bone derived from 14 week gestation femurs) within both hyaluronic acid gel and collagen foam, to compare the biocompatibility of both matrices as potential delivery systems for bone engineering and particularly for oral application. Fetal bone cell banks were prepared from one organ donation and cells were cultured for up to 4 weeks within hyaluronic acid (Mesolis(®)) and collagen foams (TissueFleece(®)). Cell survival and differentiation were assessed by cell proliferation assays and histology of frozen sections stained with Giemsa, von Kossa and ALP at 1, 2 and 4 weeks of culture. Within both materials, fetal bone cells could proliferate in three-dimensional structure at ∼70% capacity compared to monolayer culture. In addition, these cells were positive for ALP and von Kossa staining, indicating cellular differentiation and matrix production. Collagen foam provides a better structure for fetal bone cell delivery if cavity filling is necessary and hydrogels would permit an injectable technique for difficult to treat areas. In all, there was high biocompatibility, cellular differentiation and matrix deposition seen in both matrices by fetal bone cells, allowing for easy cell delivery for bone stimulation in vivo. Copyright © 2011 John Wiley & Sons, Ltd.

Quantitative proteomics identifies the membrane-associated peroxidase GPx8 as a cellular substrate of the hepatitis C virus NS3-4A protease.

The hepatitis C virus (HCV) NS3-4A protease is not only an essential component of the viral replication complex and a prime target for antiviral intervention but also a key player in the persistence and pathogenesis of HCV. It cleaves and thereby inactivates two crucial adaptor proteins in viral RNA sensing and innate immunity, mitochondrial antiviral signaling protein (MAVS) and TRIF, a phosphatase involved in growth factor signaling, T-cell protein tyrosine phosphatase (TC-PTP), and the E3 ubiquitin ligase component UV-damaged DNA-binding protein 1 (DDB1). Here we explored quantitative proteomics to identify novel cellular substrates of the NS3-4A protease. Cell lines inducibly expressing the NS3-4A protease were analyzed by stable isotopic labeling using amino acids in cell culture (SILAC) coupled with protein separation and mass spectrometry. This approach identified the membrane-associated peroxidase GPx8 as a bona fide cellular substrate of the HCV NS3-4A protease. Cleavage by NS3-4A occurs at Cys 11, removing the cytosolic tip of GPx8, and was observed in different experimental systems as well as in liver biopsies from patients with chronic HCV. Overexpression and RNA silencing studies revealed that GPx8 is involved in viral particle production but not in HCV entry or RNA replication. Conclusion: We provide proof-of-concept for the use of quantitative proteomics to identify cellular substrates of a viral protease and describe GPx8 as a novel proviral host factor targeted by the HCV NS3-4A protease. (Hepatology 2014;59:423-433).

Systems Genetics of Metabolism: The Use of the BXD Murine Reference Panel for Multiscalar Integration of Traits.

Metabolic homeostasis is achieved by complex molecular and cellular networks that differ significantly among individuals and are difficult to model with genetically engineered lines of mice optimized to study single gene function. Here, we systematically acquired metabolic phenotypes by using the EUMODIC EMPReSS protocols across a large panel of isogenic but diverse strains of mice (BXD type) to study the genetic control of metabolism. We generated and analyzed 140 classical phenotypes and deposited these in an open-access web service for systems genetics (www.genenetwork.org). Heritability, influence of sex, and genetic modifiers of traits were examined singly and jointly by using quantitative-trait locus (QTL) and expression QTL-mapping methods. Traits and networks were linked to loci encompassing both known variants and novel candidate genes, including alkaline phosphatase (ALPL), here linked to hypophosphatasia. The assembled and curated phenotypes provide key resources and exemplars that can be used to dissect complex metabolic traits and disorders.

A qualitative continuous model of cellular auxin and brassinosteroid signaling and their crosstalk.

Motivation: Hormone pathway interactions are crucial in shaping plant development, such as synergism between the auxin and brassinosteroid pathways in cell elongation. Both hormone pathways have been characterized in detail, revealing several feedback loops. The complexity of this network, combined with a shortage of kinetic data, renders its quantitative analysis virtually impossible at present.Results: As a first step towards overcoming these obstacles, we analyzed the network using a Boolean logic approach to build models of auxin and brassinosteroid signaling, and their interaction. To compare these discrete dynamic models across conditions, we transformed them into qualitative continuous systems, which predict network component states more accurately and can accommodate kinetic data as they become available. To this end, we developed an extension for the SQUAD software, allowing semi-quantitative analysis of network states. Contrasting the developmental output depending on cell type-specific modulators enabled us to identify a most parsimonious model, which explains initially paradoxical mutant phenotypes and revealed a novel physiological feature.

The blood-brain barrier: cellular basis.

Perfusion experiments with horseradish peroxidase have established that the morphological substrate of the blood-brain barrier is represented by microvascular endothelial cells. They are characterized by complexly arranged tight junctions and a very low rate of transcytotic vesicular transport. They express transport enzymes, carrier systems and brain endothelial cell-specific molecules of unknown function not expressed by any other endothelial cell population. These blood-brain barrier properties are not intrinsic to these cells but are inducible by the surrounding brain tissue. Type I astrocytes injected into the anterior eye chamber of the rat or onto the chick chorioallantoic membrane are able to induce a host-derived angiogenesis and some blood-brain barrier properties in endothelial cells of non-neural origin. Recently we have shown that this cellular interaction is due to the secretion of a soluble astrocyte derived factor(s). Astrocytes are also implicated in the maintenance, functional regulation and the repair of the blood-brain barrier. Complex interactions between other constituents of the microenvironment surrounding the endothelial cells, such as the basement membrane, pericytes, nerve endings, microglial cells and the extracellular fluid, take place and are required for the proper functioning of the blood-brain barrier, which in addition is regionally different as reflected by endothelial cell heterogeneity.