Development and application of computational methodologies in qualitative modeling

La présente thèse s'intitule "Développent et Application des Méthodologies Computationnelles pour la Modélisation Qualitative". Elle comprend tous les différents projets que j'ai entrepris en tant que doctorante. Plutôt qu'une mise en oeuvre systématique d'un cadre défini a priori, cette thèse devrait être considérée comme une exploration des méthodes qui peuvent nous aider à déduire le plan de processus regulatoires et de signalisation. Cette exploration a été mue par des questions biologiques concrètes, plutôt que par des investigations théoriques. Bien que tous les projets aient inclus des systèmes divergents (réseaux régulateurs de gènes du cycle cellulaire, réseaux de signalisation de cellules pulmonaires) ainsi que des organismes (levure à fission, levure bourgeonnante, rat, humain), nos objectifs étaient complémentaires et cohérents. Le projet principal de la thèse est la modélisation du réseau de l'initiation de septation (SIN) du S.pombe. La cytokinèse dans la levure à fission est contrôlée par le SIN, un réseau signalant de protéines kinases qui utilise le corps à pôle-fuseau comme échafaudage. Afin de décrire le comportement qualitatif du système et prédire des comportements mutants inconnus, nous avons décidé d'adopter l'approche de la modélisation booléenne. Dans cette thèse, nous présentons la construction d'un modèle booléen étendu du SIN, comprenant la plupart des composantes et des régulateurs du SIN en tant que noeuds individuels et testable expérimentalement. Ce modèle utilise des niveaux d'activité du CDK comme noeuds de contrôle pour la simulation d'évènements du SIN à différents stades du cycle cellulaire. Ce modèle a été optimisé en utilisant des expériences d'un seul "knock-out" avec des effets phénotypiques connus comme set d'entraînement. Il a permis de prédire correctement un set d'évaluation de "knock-out" doubles. De plus, le modèle a fait des prédictions in silico qui ont été validées in vivo, permettant d'obtenir de nouvelles idées de la régulation et l'organisation hiérarchique du SIN. Un autre projet concernant le cycle cellulaire qui fait partie de cette thèse a été la construction d'un modèle qualitatif et minimal de la réciprocité des cyclines dans la S.cerevisiae. Les protéines Clb dans la levure bourgeonnante présentent une activation et une dégradation caractéristique et séquentielle durant le cycle cellulaire, qu'on appelle communément les vagues des Clbs. Cet évènement est coordonné avec la courbe d'activation inverse du Sic1, qui a un rôle inhibitoire dans le système. Pour l'identification des modèles qualitatifs minimaux qui peuvent expliquer ce phénomène, nous avons sélectionné des expériences bien définies et construit tous les modèles minimaux possibles qui, une fois simulés, reproduisent les résultats attendus. Les modèles ont été filtrés en utilisant des simulations ODE qualitatives et standardisées; seules celles qui reproduisaient le phénotype des vagues ont été gardées. L'ensemble des modèles minimaux peut être utilisé pour suggérer des relations regulatoires entre les molécules participant qui peuvent ensuite être testées expérimentalement. Enfin, durant mon doctorat, j'ai participé au SBV Improver Challenge. Le but était de déduire des réseaux spécifiques à des espèces (humain et rat) en utilisant des données de phosphoprotéines, d'expressions des gènes et des cytokines, ainsi qu'un réseau de référence, qui était mis à disposition comme donnée préalable. Notre solution pour ce concours a pris la troisième place. L'approche utilisée est expliquée en détail dans le dernier chapitre de la thèse. -- The present dissertation is entitled "Development and Application of Computational Methodologies in Qualitative Modeling". It encompasses the diverse projects that were undertaken during my time as a PhD student. Instead of a systematic implementation of a framework defined a priori, this thesis should be considered as an exploration of the methods that can help us infer the blueprint of regulatory and signaling processes. This exploration was driven by concrete biological questions, rather than theoretical investigation. Even though the projects involved divergent systems (gene regulatory networks of cell cycle, signaling networks in lung cells), as well as organisms (fission yeast, budding yeast, rat, human), our goals were complementary and coherent. The main project of the thesis is the modeling of the Septation Initiation Network (SIN) in S.pombe. Cytokinesis in fission yeast is controlled by the SIN, a protein kinase signaling network that uses the spindle pole body as scaffold. In order to describe the qualitative behavior of the system and predict unknown mutant behaviors we decided to adopt a Boolean modeling approach. In this thesis, we report the construction of an extended, Boolean model of the SIN, comprising most SIN components and regulators as individual, experimentally testable nodes. The model uses CDK activity levels as control nodes for the simulation of SIN related events in different stages of the cell cycle. The model was optimized using single knock-out experiments of known phenotypic effect as a training set, and was able to correctly predict a double knock-out test set. Moreover, the model has made in silico predictions that have been validated in vivo, providing new insights into the regulation and hierarchical organization of the SIN. Another cell cycle related project that is part of this thesis was to create a qualitative, minimal model of cyclin interplay in S.cerevisiae. CLB proteins in budding yeast present a characteristic, sequential activation and decay during the cell cycle, commonly referred to as Clb waves. This event is coordinated with the inverse activation curve of Sic1, which has an inhibitory role in the system. To generate minimal qualitative models that can explain this phenomenon, we selected well-defined experiments and constructed all possible minimal models that, when simulated, reproduce the expected results. The models were filtered using standardized qualitative ODE simulations; only the ones reproducing the wave-like phenotype were kept. The set of minimal models can be used to suggest regulatory relations among the participating molecules, which will subsequently be tested experimentally. Finally, during my PhD I participated in the SBV Improver Challenge. The goal was to infer species-specific (human and rat) networks, using phosphoprotein, gene expression and cytokine data and a reference network provided as prior knowledge. Our solution to the challenge was selected as in the final chapter of the thesis.

Modelling of iron-filled magneto-active polymers with a dispersed chain-like microstructure

Magneto-active polymers are a class of smart materials commonly manufactured by mixing micron-sized iron particles in a rubber-like matrix. When cured in the presence of an externally applied magnetic field, the iron particles arrange themselves into chain-like structures that lend an overall anisotropy to the material. It has been observed through electron micrographs and X-ray tomographs that these chains are not always perfect in structure, and may have dispersion due to the conditions present during manufacturing or some undesirable material properties. We model the response of these materials to coupled magneto-mechanical loading in this paper using a probability based structure tensor that accounts for this imperfect anisotropy. The response of the matrix material is decoupled from the chain phase, though still being connected through kinematic constraints. The latter is based on the definition of a 'chain deformation gradient' and a 'chain magnetic field'. We conclude with numerical examples that demonstrate the effect of chain dispersion on the response of the material to magnetoelastic loading.

Twenty years of research into Chlamydia-like organisms: a revolution in our understanding of the biology and pathogenicity of members of the phylum Chlamydiae.

Chlamydiae are obligate intracellular bacteria that share a unique but remarkably conserved biphasic developmental cycle that relies on a eukaryotic host cell for survival. Although the phylum was originally thought to only contain one family, the Chlamydiaceae, a total of nine families are now recognized. These so-called Chlamydia-like organisms (CLOs) are also referred to as 'environmental chlamydiae', as many were initially isolated from environmental sources. However, these organisms are also emerging pathogens, as many, such as Parachlamydia sp., Simkania sp. and Waddlia sp., have been associated with human disease, and others, such as Piscichlamydia sp. and Parilichlamydia sp., have been documented in association with diseases in animals. Their strict intracellular nature and the requirement for cell culture have been a confounding factor in characterizing the biology and pathogenicity of CLOs. Nevertheless, the genomes of seven CLO species have now been sequenced, providing new information on their potential ability to adapt to a wide range of hosts. As new isolation and diagnostic methods advance, we are able to further explore the richness of this phylum with further research likely to help define the true pathogenic potential of the CLOs while also providing insight into the origins of the 'traditional' chlamydiae.

Effects of a post-weaning cafeteria diet in young rats: metabolic syndrome, reduced activity and low anxiety-like behaviour.

Among adolescents, overweight, obesity and metabolic syndrome are rapidly increasing in recent years as a consequence of unhealthy palatable diets. Animal models of diet-induced obesity have been developed, but little is known about the behavioural patterns produced by the consumption of such diets. The aim of the present study was to determine the behavioural and biochemical effects of a cafeteria diet fed to juvenile male and female rats, as well as to evaluate the possible recovery from these effects by administering standard feeding during the last week of the study. Two groups of male and female rats were fed with either a standard chow diet (ST) or a cafeteria (CAF) diet from weaning and for 8 weeks. A third group of males (CAF withdrawal) was fed with the CAF diet for 7 weeks and the ST in the 8th week. Both males and females developed metabolic syndrome as a consequence of the CAF feeding, showing overweight, higher adiposity and liver weight, increased plasma levels of glucose, insulin and triglycerides, as well as insulin resistance, in comparison with their respective controls. The CAF diet reduced motor activity in all behavioural tests, enhanced exploration, reduced anxiety-like behaviour and increased social interaction; this last effect was more pronounced in females than in males. When compared to animals only fed with a CAF diet, CAF withdrawal increased anxiety in the open field, slightly decreased body weight, and completely recovered the liver weight, insulin sensitivity and the standard levels of glucose, insulin and triglycerides in plasma. In conclusion, a CAF diet fed to young animals for 8 weeks induced obesity and metabolic syndrome, and produced robust behavioural changes in young adult rats, whereas CAF withdrawal in the last week modestly increased anxiety, reversed the metabolic alterations and partially reduced overweight.

N=1 SQCD-like theories with N_f massive flavors from AdS/CFT and beta functions

We study new supergravity solutions related to large-N c N=1 supersymmetric gauge field theories with a large number N f of massive flavors. We use a recently proposed framework based on configurations with N c color D5 branes and a distribution of N f flavor D5 branes, governed by a function N f S(r). Although the system admits many solutions, under plausible physical assumptions the relevant solution is uniquely determined for each value of x ≡ N f /N c . In the IR region, the solution smoothly approaches the deformed Maldacena-Núñez solution. In the UV region it approaches a linear dilaton solution. For x < 2 the gauge coupling β g function computed holographically is negative definite, in the UV approaching the NSVZ β function with anomalous dimension γ 0 = −1/2 (approaching − 3/(32π 2)(2N c  − N f )g 3)), and with β g  → −∞ in the IR. For x = 2, β g has a UV fixed point at strong coupling, suggesting the existence of an IR fixed point at a lower value of the coupling. We argue that the solutions with x > 2 describe a"Seiberg dual" picture where N f  − 2N c flips sign.

Protein oligomers studied by solid-state NMR the case of the full-length nucleoid-associated protein histone-like nucleoid structuring protein

Members of the histone-like nucleoid structuring protein (H-NS) family play roles both as architectural proteins and as modulators of gene expression in Gram-negative bacteria. The H-NS protein participates in modulatory processes that respond to environmental changes in osmolarity, pH, or temperature. H-NS oligomerization is essential for its activity. Structural models of different truncated forms are available. However, high-resolution structural details of full-length H-NS and its DNA-bound state have largely remained elusive. We report on progress in characterizing the biologically active H-NS oligomers with solid-state NMR. We compared uniformly ((13)C,(15)N)-labeled ssNMR preparations of the isolated N-terminal region (H-NS 1-47) and full-length H-NS (H-NS 1-137). In both cases, we obtained ssNMR spectra of good quality and characteristic of well-folded proteins. Analysis of the results of 2D and 3D (13)C-(13)C and (15)N-(13)C correlation experiments conducted at high magnetic field led to assignments of residues located in different topological regions of the free full-length H-NS. These findings confirm that the structure of the N-terminal dimerization domain is conserved in the oligomeric full-length protein. Small changes in the dimerization interface suggested by localized chemical shift variations between solution and solid-state spectra may be relevant for DNA recoginition.

Toll-like receptor 5 deficiency exacerbates cardiac injury and inflammation induced by myocardial ischaemia-reperfusion in the mouse.

Myocardial ischaemia-reperfusion (MIR) triggers a sterile inflammatory response important for myocardial healing, but which may also contribute to adverse ventricular remodelling. Such inflammation is initiated by molecular danger signals released by damaged myocardium, which induce innate immune responses by activating toll-like receptors (TLRs). Detrimental roles have been recently reported for TLR2, TLR3 and TLR4. The role of other TLRs is unknown. We therefore evaluated the role of TLR5, expressed at high level in the heart, in the development of myocardial damage and inflammation acutely triggered by MIR. TLR5-/- and wild-type (WT) mice were exposed to MIR (30 min ischaemia, 2 h reperfusion). We measured infarct size, markers of cardiac oxidative stress, myocardial phosphorylation state of mitogen-activated protein (MAP) kinases and AKT, expression levels of chemokines and cytokines in the heart and plasma, as well as cardiac function by echography and conductance volumetry. TLR5-deficient mice had normal cardiac morphology and function under physiological conditions. After MIR, the absence of TLR5 promoted an increase in infarct size and myocardial oxidative stress. Lack of TLR5 fostered p38 phosphorylation, reduced AKT phosphorylation and markedly increased the expression of inflammatory cytokines, whereas it precipitated acute LV (left ventricle) dysfunction. Therefore, contrary to the detrimental roles of TLR2, TLR3 and TLR4 in the infarcted heart, TLR5 is important to limit myocardial damage, inflammation and functional compromise after MIR.

Patterns of nucleotide diversity at the regions encompassing the Drosophila insulin-like peptide (dilp) genes: demography vs positive selection in Drosophila melanogaster.

In Drosophila, the insulin-signaling pathway controls some life history traits, such as fertility and lifespan, and it is considered to be the main metabolic pathway involved in establishing adult body size. Several observations concerning variation in body size in the Drosophila genus are suggestive of its adaptive character. Genes encoding proteins in this pathway are, therefore, good candidates to have experienced adaptive changes and to reveal the footprint of positive selection. The Drosophila insulin-like peptides (DILPs) are the ligands that trigger the insulin-signaling cascade. In Drosophila melanogaster, there are several peptides that are structurally similar to the single mammalian insulin peptide. The footprint of recent adaptive changes on nucleotide variation can be unveiled through the analysis of polymorphism and divergence. With this aim, we have surveyed nucleotide sequence variation at the dilp1-7 genes in a natural population of D. melanogaster. The comparison of polymorphism in D. melanogaster and divergence from D. simulans at different functional classes of the dilp genes provided no evidence of adaptive protein evolution after the split of the D. melanogaster and D. simulans lineages. However, our survey of polymorphism at the dilp gene regions of D. melanogaster has provided some evidence for the action of positive selection at or near these genes. The regions encompassing the dilp1-4 genes and the dilp6 gene stand out as likely affected by recent adaptive events.

Sustained CTL activation by murine pulmonary epithelial cells promotes the development of COPD-like disease

Chronic obstructive pulmonary disease (COPD) is a lethal progressive lung disease culminating in permanent airway obstruction and alveolar enlargement. Previous studies suggest CTL involvement in COPD progression; however, their precise role remains unknown. Here, we investigated whether the CTL activation receptor NK cell group 2D (NKG2D) contributes to the development of COPD. Using primary murine lung epithelium isolated from mice chronically exposed to cigarette smoke and cultured epithelial cells exposed to cigarette smoke extract in vitro, we demonstrated induced expression of the NKG2D ligand retinoic acid early tran - script 1 (RAET1)as well as NKG2D-mediated cytotoxicity. Furthermore, a genetic model of inducible RAET1 expression on mouse pulmonary epithelial cells yielded a severe emphysematous phenotype characterized by epithelial apoptosis and increased CTL activation, which was reversed by blocking NKG2D activation. We also assessed whether NKG2D ligand expression corresponded with pulmonary disease in human patients by staining airway and peripheral lung tissues from never smokers, smokers with normal lung function, and current and former smokers with COPD. NKG2D ligand expression was independent of NKG2D receptor expression in COPD patients, demonstrating that ligand expression is the limiting factor in CTL activation. These results demonstrate that aberrant, persistent NKG2D ligand expression in the pulmonary epithelium contributes to the development of COPD pathologies.

Innate Immune Sensing of Fusarium culmorum by Mouse Dendritic Cells.

Chronic inhalation of grain dust is associated with asthma and chronic bronchitis in grain worker populations. Exposure to fungal particles was postulated to be an important etiologic agent of these pathologies. Fusarium species frequently colonize grain and straw and produce a wide array of mycotoxins that impact human health, necessitating an evaluation of risk exposure by inhalation of Fusarium and its consequences on immune responses. Data showed that Fusarium culmorum is a frequent constituent of aerosols sampled during wheat harvesting in the Vaud region of Switzerland. The aim of this study was to examine cytokine/chemokine responses and innate immune sensing of F. culmorum in bone-marrow-derived dendritic cells and macrophages. Overall, dendritic cells and macrophages responded to F. culmorum spores but not to its secreted components (i.e., mycotoxins) by releasing large amounts of macrophage inflammatory protein (MIP)-1α, MIP-1β, MIP-2, monocyte chemoattractant protein (MCP)-1, RANTES, and interleukin (IL)-12p40, intermediate amounts of tumor necrosis factor (TNF), IL-6, IL-12p70, IL-33, granulocyte colony-stimulating factor (G-CSF), and interferon gamma-induced protein (IP-10), but no detectable amounts of IL-4 and IL-10, a pattern of mediators compatible with generation of Th1 or Th17 antifungal protective immune responses rather than with Th2-dependent allergic responses. The sensing of F. culmorum spores by dendritic cells required dectin-1, the main pattern recognition receptor involved in β-glucans detection, but likely not MyD88 and TRIF-dependent Toll-like receptors. Taken together, our results indicate that F. culmorum stimulates potently innate immune cells in a dectin-1-dependent manner, suggesting that inhalation of F. culmorum from grain dust may promote immune-related airway diseases in exposed worker populations.

The development of decision limits for the GH-2000 detection methodology using additional insulin-like growth factor-I and amino-terminal pro-peptide of type III collagen assays.

The GH-2000 and GH-2004 projects have developed a method for detecting GH misuse based on measuring insulin-like growth factor-I (IGF-I) and the amino-terminal pro-peptide of type III collagen (P-III-NP). The objectives were to analyze more samples from elite athletes to improve the reliability of the decision limit estimates, to evaluate whether the existing decision limits needed revision, and to validate further non-radioisotopic assays for these markers. The study included 998 male and 931 female elite athletes. Blood samples were collected according to World Anti-Doping Agency (WADA) guidelines at various sporting events including the 2011 International Association of Athletics Federations (IAAF) World Athletics Championships in Daegu, South Korea. IGF-I was measured by the Immunotech A15729 IGF-I IRMA, the Immunodiagnostic Systems iSYS IGF-I assay and a recently developed mass spectrometry (LC-MS/MS) method. P-III-NP was measured by the Cisbio RIA-gnost P-III-P, Orion UniQ? PIIINP RIA and Siemens ADVIA Centaur P-III-NP assays. The GH-2000 score decision limits were developed using existing statistical techniques. Decision limits were determined using a specificity of 99.99% and an allowance for uncertainty because of the finite sample size. The revised Immunotech IGF-I - Orion P-III-NP assay combination decision limit did not change significantly following the addition of the new samples. The new decision limits are applied to currently available non-radioisotopic assays to measure IGF-I and P-III-NP in elite athletes, which should allow wider flexibility to implement the GH-2000 marker test for GH misuse while providing some resilience against manufacturer withdrawal or change of assays. Copyright © 2015 John Wiley & Sons, Ltd.

Quantitative proteomics in the characterization of T helper lymphocyte differentiation

The term proteome is used to define the complete set of proteins expressed in cells or tissues of an organism at a certain timepoint. Respectively, proteomics is used to describe the methods, which are used to study such proteomes. These methods include chromatographic and electrophoretic techniques for protein or peptide fractionation, mass spectrometry for their identification, and use of computational methods to assist the complicated data analysis. A primary aim in this Ph.D. thesis was to set-up, optimize, and develop proteomics methods for analysing proteins extracted from T-helper (Th) lymphocytes. First, high-throughput LC-MS/MS and ICAT labeling methods were set-up and optimized for analysing the microsomal fraction proteins extracted from Th lymphocytes. Later, iTRAQ method was optimized to study cytokine regulated protein expression in the nuclei of Th lymphocytes. High-throughput LC-MS/MS analyses, like ICAT and iTRAQ, produce large quantities of data and robust software and data analysis pipelines are needed. Therefore, different software programs used for analysing such data were evaluated. Moreover, a pre-filtering algorithm was developed to classify good-quality and bad-quality spectra prior to the database searches. Th-lymphocytes can differentiate into Th1 or Th2 cells based on surrounding antigens, co-stimulatory molecules, and cytokines. Both subsets have individual cytokine secretion profiles and specific functions. Th1 cells participate in the cellular immunity against intracellular pathogens, while Th2 cells have important role in the humoral immunity against extracellular parasites. An abnormal response of Th1 and Th2 cells and imbalance between the subsets are charasteristic of several diseases. Th1 specific reactions and cytokines have been detected in autoimmune diseases, while Th2 specific response and cytokine profile is common in allergy and asthma. In this Ph. D. thesis mass spectrometry-based proteomics was used to study the effects of Th1 and Th2 promoting cytokines IL-12 and IL-4 on the proteome of Th lymphocytes. Characterization of microsomal fraction proteome extracted from IL-12 treated lymphobasts and IL-4 stimulated cord blood CD4+ cells resulted in finding of cytokine regulated proteins. Galectin-1 and CD7 were down-regulated in IL-12 treated cells, while IL-4 stimulation decreased the expression of STAT1, MXA, GIMAP1, and GIMAP4. Interestingly, the transcription of both GIMAP genes was up-regulated in Th1 polarized cells and down-regulated in Th2 promoting conditions.

TH17 cells promote microbial killing and innate immune sensing of DNA via interleukin 26.

Interleukin 17-producing helper T cells (TH17 cells) have a major role in protection against infections and in mediating autoimmune diseases, yet the mechanisms involved are incompletely understood. We found that interleukin 26 (IL-26), a human TH17 cell-derived cytokine, is a cationic amphipathic protein that kills extracellular bacteria via membrane-pore formation. Furthermore, TH17 cell-derived IL-26 formed complexes with bacterial DNA and self-DNA released by dying bacteria and host cells. The resulting IL-26-DNA complexes triggered the production of type I interferon by plasmacytoid dendritic cells via activation of Toll-like receptor 9, but independently of the IL-26 receptor. These findings provide insights into the potent antimicrobial and proinflammatory function of TH17 cells by showing that IL-26 is a natural human antimicrobial that promotes immune sensing of bacterial and host cell death.

A frequent hypofunctional IRAK2 variant is associated with reduced spontaneous hepatitis C virus clearance.

UNLABELLED: Patients carrying very rare loss-of-function mutations in interleukin-1 receptor-associated kinase 4 (IRAK4), a critical signaling mediator in Toll-like receptor signaling, are severely immunodeficient, highlighting the paramount role of IRAK kinases in innate immunity. We discovered a comparatively frequent coding variant of the enigmatic human IRAK2, L392V (rs3844283), which is found homozygously in ∼15% of Caucasians, to be associated with a reduced ability to induce interferon-alpha in primary human plasmacytoid dendritic cells in response to hepatitis C virus (HCV). Cytokine production in response to purified Toll-like receptor agonists was also impaired. Additionally, rs3844283 was epidemiologically associated with a chronic course of HCV infection in two independent HCV cohorts and emerged as an independent predictor of chronic HCV disease. Mechanistically, IRAK2 L392V showed intact binding to, but impaired ubiquitination of, tumor necrosis factor receptor-associated factor 6, a vital step in signal transduction. CONCLUSION: Our study highlights IRAK2 and its genetic variants as critical factors and potentially novel biomarkers for human antiviral innate immunity. (Hepatology 2015;62:1375-1387).

Papillary cystoadenoma lymphomatosum (Warthin-like) of minor salivary glands

Papillary cystadenoma lymphomatosum is a benign salivary gland tumor most frequently located in the parotid gland (Warthin"s tumor). Its presentation in other major, or in minor, salivary glands is rare. Clinically, it manifests as a slow growing tumor, fluctuant on palpation due to its cystic morphology. The treatment of choice is complete excision with wide tumor-free margins. We present a 73-year-old female patient with an asymptomatic tumor of 8 years evolution in the right posterior area of the hard palate. We performed surgical excision and a biopsy, which was reported as papillary cystadenoma lymphomatosum. During the post-operative examination carried out after 3 weeks, it was observed that the lesion had recurred. The lesion was re-operated, performing the excision with CO2 laser and including the periosteum to ensure complete resection of the tumor. At 10 months follow-up, there was no recurrence of the lesion. This article includes a review of this condition and discusses its most important clinical and pathologic features and therapeutic approaches.