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.

Deux approches de l'influence du discours économique sur les phénomènes sociaux

The influence of theoretical discourse on the corresponding field of investigation is an important feature of social science : social scientists shape the world by describing it. This phenomenon has been studied in several ways. Economic sociology has recently focused on the fact that the economy is embedded in the economics through which it is scrutinized. Great names of economic theory have focused on the effects of economists' discourse on agents' behaviours. This article aims to bring out the distinction between these two kinds of contributions which deal with the same object. Finally, it explains the analytical distinction by a difference in the initial problematic. Following trends in the French school of the economics of convention, the social influence of theoretical discourse is analyzed on the strategic side (in connection with economics) and on the interpretative side (in connection with sociology).

Use of DNA profiles for investigation using a simulated national DNA database: Part I. Partial SGM Plus profiles.

In traditional criminal investigation, uncertainties are often dealt with using a combination of common sense, practical considerations and experience, but rarely with tailored statistical models. For example, in some countries, in order to search for a given profile in the national DNA database, it must have allelic information for six or more of the ten SGM Plus loci for a simple trace. If the profile does not have this amount of information then it cannot be searched in the national DNA database (NDNAD). This requirement (of a result at six or more loci) is not based on a statistical approach, but rather on the feeling that six or more would be sufficient. A statistical approach, however, could be more rigorous and objective and would take into consideration factors such as the probability of adventitious matches relative to the actual database size and/or investigator's requirements in a sensible way. Therefore, this research was undertaken to establish scientific foundations pertaining to the use of partial SGM Plus loci profiles (or similar) for investigation.

Clinical and genetic investigation of a large Tunisian family with complete achromatopsia: identification of a new nonsense mutation in GNAT2 gene.

Complete achromatopsia is a rare autosomal recessive disease associated with CNGA3, CNGB3, GNAT2 and PDE6C mutations. This retinal disorder is characterized by complete loss of color discrimination due to the absence or alteration of the cones function. The purpose of the present study was the clinical and the genetic characterization of achromatopsia in a large consanguineous Tunisian family. Ophthalmic evaluation included a full clinical examination, color vision testing and electroretinography. Linkage analysis using microsatellite markers flanking CNGA3, CNGB3, GNAT2 and PDE6C genes was performed. Mutations were screened by direct sequencing. A total of 12 individuals were diagnosed with congenital complete achromatopsia. They are members of six nuclear consanguineous families belonging to the same large consanguineous family. Linkage analysis revealed linkage to GNAT2. Mutational screening of GNAT2 revealed three intronic variations c.119-69G>C, c.161+66A>T and c.875-31G>C that co-segregated with a novel mutation p.R313X. An identical GNAT2 haplotype segregating with this mutation was identified, indicating a founder mutation. All patients were homozygous for the p.R313X mutation. This is the first report of the clinical and genetic investigation of complete achromatopsia in North Africa and the largest family with recessive achromatopsia involving GNAT2; thus, providing a unique opportunity for genotype-phenotype correlation for this extremely rare condition.

Subcellular investigation of photosynthesis-driven carbon assimilation in the symbiotic reef coral Pocillopora damicornis.

Reef-building corals form essential, mutualistic endosymbiotic associations with photosynthetic Symbiodinium dinoflagellates, providing their animal host partner with photosynthetically derived nutrients that allow the coral to thrive in oligotrophic waters. However, little is known about the dynamics of these nutritional interactions at the (sub)cellular level. Here, we visualize with submicrometer spatial resolution the carbon and nitrogen fluxes in the intact coral-dinoflagellate association from the reef coral Pocillopora damicornis by combining nanoscale secondary ion mass spectrometry (NanoSIMS) and transmission electron microscopy with pulse-chase isotopic labeling using [(13)C]bicarbonate and [(15)N]nitrate. This allows us to observe that (i) through light-driven photosynthesis, dinoflagellates rapidly assimilate inorganic bicarbonate and nitrate, temporarily storing carbon within lipid droplets and starch granules for remobilization in nighttime, along with carbon and nitrogen incorporation into other subcellular compartments for dinoflagellate growth and maintenance, (ii) carbon-containing photosynthates are translocated to all four coral tissue layers, where they accumulate after only 15 min in coral lipid droplets from the oral gastroderm and within 6 h in glycogen granules from the oral epiderm, and (iii) the translocation of nitrogen-containing photosynthates is delayed by 3 h. IMPORTANCE: Our results provide detailed in situ subcellular visualization of the fate of photosynthesis-derived carbon and nitrogen in the coral-dinoflagellate endosymbiosis. We directly demonstrate that lipid droplets and glycogen granules in the coral tissue are sinks for translocated carbon photosynthates by dinoflagellates and confirm their key role in the trophic interactions within the coral-dinoflagellate association.

Native-like, long synthetic peptides as components of sub-unit vaccines: practical and theoretical considerations for their use in humans.

Vaccines have been used as a successful tool in medicine by way of controlling many major diseases. In spite of this, vaccines today represent only a handful of all infectious diseases. Therefore, there is a pressing demand for improvements of existing vaccines with particular reference to higher efficacy and undisputed safety profiles. To this effect, as an alternative to available vaccine technologies, there has been a drive to develop vaccine candidate polypeptides by chemical synthesis. In our laboratory, we have recently developed a technology to manufacture long synthetic peptides of up to 130 residues, which are correctly folded and biologically active. This paper discusses the advantages of the molecularly defined, long synthetic peptide approach in the context of vaccine design, development and use in human vaccination.

Theoretical study on receptor-G protein recognition: new insights into the mechanisms of the alpha1b-adrenergic recptor activation

This work compares the structural/dynamics features of the wild-type alb-adrenergic receptor (AR) with those of the D142A active mutant and the agonist-bound state. The two active receptor forms were compared in their isolated states as well as in their ability to form homodimers and to recognize the G alpha q beta 1 gamma 2 heterotrimer. The analysis of the isolated structures revealed that, although the mutation- and agonist-induced active states of the alpha 1b-AR are different, they, however, share several structural peculiarities including (a) the release of some constraining interactions found in the wild-type receptor and (b) the opening of a cytosolic crevice formed by the second and third intracellular loops and the cytosolic extensions of helices 5 and 6. Accordingly, also their tendency to form homodimers shows commonalties and differences. In fact, in both the active receptor forms, helix 6 plays a crucial role in mediating homodimerization. However, the homodimeric models result from different interhelical assemblies. On the same line of evidence, in both of the active receptor forms, the cytosolic opened crevice recognizes similar domains on the G protein. However, the docking solutions are differently populated and the receptor-G protein preorientation models suggest that the final complexes should be characterized by different interaction patterns.

Investigation of field and diffusion time dependence of the diffusion-weighted signal at ultrahigh magnetic fields.

Over the last decade, there has been a significant increase in the number of high-magnetic-field MRI magnets. However, the exact effect of a high magnetic field strength (B0 ) on diffusion-weighted MR signals is not yet fully understood. The goal of this study was to investigate the influence of different high magnetic field strengths (9.4 T and 14.1 T) and diffusion times (9, 11, 13, 15, 17 and 24 ms) on the diffusion-weighted signal in rat brain white matter. At a short diffusion time (9 ms), fractional anisotropy values were found to be lower at 14.1 T than at 9.4 T, but this difference disappeared at longer diffusion times. A simple two-pool model was used to explain these findings. The model describes the white matter as a first hindered compartment (often associated with the extra-axonal space), characterized by a faster orthogonal diffusion and a lower fractional anisotropy, and a second restricted compartment (often associated with the intra-axonal space), characterized by a slower orthogonal diffusion (i.e. orthogonal to the axon direction) and a higher fractional anisotropy. Apparent T2 relaxation time measurements of the hindered and restricted pools were performed. The shortening of the pseudo-T2 value from the restricted compartment with B0 is likely to be more pronounced than the apparent T2 changes in the hindered compartment. This study suggests that the observed differences in diffusion tensor imaging parameters between the two magnetic field strengths at short diffusion time may be related to differences in the apparent T2 values between the pools. Copyright © 2013 John Wiley & Sons, Ltd.

Investigation of Postmortem Gases: A Forensic Imaging and Analytical Chemistry Approach

Background: Distinguishing postmortem gas accumulations in the body due to natural decomposition and other phenomena such as gas embolism can prove a difficult task using purely Multi-Detector Computed Tomography (MDCT). The Radiological Alteration Index (RAI) was created with the intention to be able to identify bodies undergoing the putrefaction process based on the quantity of gas detected within the body. The flaw in this approach is the inability to absolutely determine putrefaction as the origin of gas volumes in cases of moderate alteration. The aim of the current study is to identify percentage compositions of O2, N2, CO2 and the presence of gases such as H2 and H2S within these sampling sites in order to resolve this complication. Materials and methods: All cases investigated in our University Center of Legal Medicine are undergoing a Post-Mortem Computed Tomography (PMCT)-scan before external examination or autopsy as a routine investigation. In the obtained images, areas of gas were characterized as 0, I, II or III based on the amount of gas present according to the RAI (1). The criteria for these characterizations were dependent of the site of gas, for example thoracic and abdominal cavities were graded as I (1 - 3cm gas), II (3 - 5cm gas) and III (>5cm gas). Cases showing gaseous sites with grade II or III were selected for this study. The sampling was performed under CT-guidance to target the regions to be punctured. Luer-lock PTFE syringes equipped with a three-way valve and needles were used to sample the gas directly (2). Gaseous samples were then analysed using gas chromatography coupled to a thermal conductivity detector (GC-TCD). The components present in the samples were expressed as a percentage of the overall gas present. Results: Up to now, we have investigated more than 40 cases using our standardized procedure for sampling and analysis of gas. O2, N2 and CO2 were present in most samples. The following distributions were found to correlate to gas origins of gas embolism/scuba diving accidents, trauma and putrefaction: ? Putrefaction → O2 = 1 - 5%; CO2 > 15%; N2 = 10 - 70%; H2 / H2S / CH4 variable presence ? Gas embolism/Scuba diving accidents → O2 and N2= varying percentages; CO2 > 20% ? Trauma → O2 = small percentage; CO2 < 15%; N2 > 65% H2 and H2S indicated levels of putrefaction along with methane which can also gauge environmental conditions or conditions of body storage/burial. Many cases showing large RAI values (advanced alteration) did reveal a radiological diagnosis which was in concordance with the interpretation of the gas composition. However, in certain cases (gas embolism, scuba divers) radiological interpretation was not possible and only chemical gas analysis was found to lead to the correct diagnosis, meaning that it provided complementary information to the radiological diagnosis. Conclusion: Investigation of postmortem gases is a useful tool to determine origin of gas generation which can aid the diagnosis of the cause of death. Levels of gas can provide information on stage of putrefaction and help to perform essential medico-legal diagnosis such as vital gas embolism.

Terrestrial Lidar investigation of the 2004 rockslide along Petit Champlain Street, Québec City (Quebec, Canada)

Au cours des 240 dernières années, 53 mouvements de versant se sont produits le long du promontoire de Québec, causant la mort de 88 personnes principalement lors de chutes de blocs. En octobre 2004, un petit éboulement a atteint la route dans une zone proche de l'éboulement de 1889 qui a tué 35 personnes et blessé 30 autres. Une image 3D a été créée par l'utilisation d'un scanner Lidar terrestre (SLT). Les sept familles de joints identifiées sont en accord avec les mesures effectuées dans de précédentes études. L'imagerie SLT a aussi permit d'estimer les volumes des instabilités passées et d'en analyser le mécanisme : un glissement rocheux qui affecte des blocs débités en parallélépipèdes par d'autres familles de joints. De plus la zone étudiée montre qu'elle est favorable aux chutes de blocs.

A collaborative approach for incorporating forensic case data into crime investigation using criminal intelligence analysis and visualisation

There is an increasing awareness that the articulation of forensic science and criminal investigation is critical to the resolution of crimes. However, models and methods to support an effective collaboration between these partners are still poorly expressed or even lacking. Three propositions are borrowed from crime intelligence methods in order to bridge this gap: (a) the general intelligence process, (b) the analyses of investigative problems along principal perspectives: entities and their relationships, time and space, quantitative aspects and (c) visualisation methods as a mode of expression of a problem in these dimensions. Indeed, in a collaborative framework, different kinds of visualisations integrating forensic case data can play a central role for supporting decisions. Among them, link-charts are scrutinised for their abilities to structure and ease the analysis of a case by describing how relevant entities are connected. However, designing an informative chart that does not bias the reasoning process is not straightforward. Using visualisation as a catalyser for a collaborative approach integrating forensic data thus calls for better specifications.