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.

Application of contrast media in post-mortem imaging (CT and MRI).

The application of contrast media in post-mortem radiology differs from clinical approaches in living patients. Post-mortem changes in the vascular system and the absence of blood flow lead to specific problems that have to be considered for the performance of post-mortem angiography. In addition, interpreting the images is challenging due to technique-related and post-mortem artefacts that have to be known and that are specific for each applied technique. Although the idea of injecting contrast media is old, classic methods are not simply transferable to modern radiological techniques in forensic medicine, as they are mostly dedicated to single-organ studies or applicable only shortly after death. With the introduction of modern imaging techniques, such as post-mortem computed tomography (PMCT) and post-mortem magnetic resonance (PMMR), to forensic death investigations, intensive research started to explore their advantages and limitations compared to conventional autopsy. PMCT has already become a routine investigation in several centres, and different techniques have been developed to better visualise the vascular system and organ parenchyma in PMCT. In contrast, the use of PMMR is still limited due to practical issues, and research is now starting in the field of PMMR angiography. This article gives an overview of the problems in post-mortem contrast media application, the various classic and modern techniques, and the issues to consider by using different media.

Single centre experience of the application of self navigated 3D whole heart cardiovascular magnetic resonance for the assessment of cardiac anatomy in congenital heart disease.

BACKGROUND: For free-breathing cardiovascular magnetic resonance (CMR), the self-navigation technique recently emerged, which is expected to deliver high-quality data with a high success rate. The purpose of this study was to test the hypothesis that self-navigated 3D-CMR enables the reliable assessment of cardiovascular anatomy in patients with congenital heart disease (CHD) and to define factors that affect image quality. METHODS: CHD patients ≥2 years-old and referred for CMR for initial assessment or for a follow-up study were included to undergo a free-breathing self-navigated 3D CMR at 1.5T. Performance criteria were: correct description of cardiac segmental anatomy, overall image quality, coronary artery visibility, and reproducibility of great vessels diameter measurements. Factors associated with insufficient image quality were identified using multivariate logistic regression. RESULTS: Self-navigated CMR was performed in 105 patients (55% male, 23 ± 12y). Correct segmental description was achieved in 93% and 96% for observer 1 and 2, respectively. Diagnostic quality was obtained in 90% of examinations, and it increased to 94% if contrast-enhanced. Left anterior descending, circumflex, and right coronary arteries were visualized in 93%, 87% and 98%, respectively. Younger age, higher heart rate, lower ejection fraction, and lack of contrast medium were independently associated with reduced image quality. However, a similar rate of diagnostic image quality was obtained in children and adults. CONCLUSION: In patients with CHD, self-navigated free-breathing CMR provides high-resolution 3D visualization of the heart and great vessels with excellent robustness.

Early versus delayed application of Thomas splints in patients with isolated femur shaft fractures: The benefits quantified.

AIMS: To investigate and quantify the clinical benefits of early versus delayed application of Thomas splints in patients with isolated femur shaft fractures. MATERIALS AND METHODS: Level IV retrospective clinical and radiological analysis of patients presenting from January to December 2012 at a Level 1 Trauma Unit. All skeletally mature patients with isolated femur shaft fractures independently of their mechanism of injury were included. Exclusion criteria were: ipsilateral fracture of the lower limb, neck and supracondylar femur fractures, periprosthetic and incomplete fractures. Their clinical records were analysed for blood transfusion requirements, pulmonary complications, surgery time, duration of hospital stay and analgesic requirements. RESULTS: A total of 106 patients met our inclusion criteria. There were 74 males and 32 females. Fifty seven (54%) patients were in the 'early splinted' group and 49 patients (46%) were in the 'delayed splinted' group (P>0.05). The need for blood transfusion was significantly reduced in the 'early splinted' group (P=0.04). There was a significantly higher rate of pulmonary complications in the 'delayed splinted' group (P=0.008). All other parameters were similar between the two groups. CONCLUSION: The early application of Thomas splints for isolated femur fractures in non-polytraumatised patients has a clinically and statistically significant benefit of reducing the need for blood transfusions and the incidence of pulmonary complications.

Clinical Application of Prognostic Gene Expression Signature in Fusion Gene-Negative Rhabdomyosarcoma: A Report from the Children's Oncology Group.

PURPOSE: Pediatric rhabdomyosarcoma (RMS) has two common histologic subtypes: embryonal (ERMS) and alveolar (ARMS). PAX-FOXO1 fusion gene status is a more reliable prognostic marker than alveolar histology, whereas fusion gene-negative (FN) ARMS patients are clinically similar to ERMS patients. A five-gene expression signature (MG5) previously identified two diverse risk groups within the fusion gene-negative RMS (FN-RMS) patients, but this has not been independently validated. The goal of this study was to test whether expression of the MG5 metagene, measured using a technical platform that can be applied to routine pathology material, would correlate with outcome in a new cohort of patients with FN-RMS. EXPERIMENTAL DESIGN: Cases were taken from the Children's Oncology Group (COG) D9803 study of children with intermediate-risk RMS, and gene expression profiling for the MG5 genes was performed using the nCounter assay. The MG5 score was correlated with clinical and pathologic characteristics as well as overall and event-free survival. RESULTS: MG5 standardized score showed no significant association with any of the available clinicopathologic variables. The MG5 signature score showed a significant correlation with overall (N = 57; HR, 7.3; 95% CI, 1.9-27.0; P = 0.003) and failure-free survival (N = 57; HR, 6.1; 95% CI, 1.9-19.7; P = 0.002). CONCLUSIONS: This represents the first, validated molecular prognostic signature for children with FN-RMS who otherwise have intermediate-risk disease. The capacity to measure the expression of a small number of genes in routine pathology material and apply a simple mathematical formula to calculate the MG5 metagene score provides a clear path toward better risk stratification in future prospective clinical trials. Clin Cancer Res; 21(20); 4733-9. ©2015 AACR.

Application of intravoxel incoherent motion perfusion imaging to shoulder muscles after a lift-off test of varying duration.

Intravoxel incoherent motion (IVIM) MRI is a method to extract microvascular blood flow information out of diffusion-weighted images acquired at multiple b-values. We hypothesized that IVIM can identify the muscles selectively involved in a specific task, by measuring changes in activity-induced local muscular perfusion after exercise. We tested this hypothesis using a widely used clinical maneuver, the lift-off test, which is known to assess specifically the subscapularis muscle functional integrity. Twelve shoulders from six healthy male volunteers were imaged at 3 T, at rest, as well as after a lift-off test hold against resistance for 30 s, 1 and 2 min respectively, in three independent sessions. IVIM parameters, consisting of perfusion fraction (f), diffusion coefficient (D), pseudo-diffusion coefficient D* and blood flow-related fD*, were estimated within outlined muscles of the rotator cuff and the deltoid bundles. The mean values at rest and after the lift-off tests were compared in each muscle using a one-way ANOVA. A statistically significant increase in fD* was measured in the subscapularis, after a lift-off test of any duration, as well as in D. A fD* increase was the most marked (30 s, +103%; 1 min, +130%; 2 min, +156%) and was gradual with the duration of the test (in 10(-3) mm(2) /s: rest, 1.41 ± 0.50; 30 s, 2.86 ± 1.17; 1 min, 3.23 ± 1.22; 2 min, 3.60 ± 1.21). A significant increase in fD* and D was also visible in the posterior bundle of the deltoid. No significant change was consistently visible in the other investigated muscles of the rotator cuff and the other bundles of the deltoid. In conclusion, IVIM fD* allows the demonstration of a task-related microvascular perfusion increase after a specific task and suggests a direct relationship between microvascular perfusion and the duration of the effort. It is a promising method to investigate non-invasively skeletal muscle physiology and clinical perfusion-related muscular disorders.

Characterization and application of two RANK-specific antibodies with different biological activities.

Antibodies play an important role in therapy and investigative biomedical research. The TNF-family member Receptor Activator of NF-κB (RANK) is known for its role in bone homeostasis and is increasingly recognized as a central player in immune regulation and epithelial cell activation. However, the study of RANK biology has been hampered by missing or insufficient characterization of high affinity tools that recognize RANK. Here, we present a careful description and comparison of two antibodies, RANK-02 obtained by phage display (Newa, 2014 [1]) and R12-31 generated by immunization (Kamijo, 2006 [2]). We found that both antibodies recognized mouse RANK with high affinity, while RANK-02 and R12-31 recognized human RANK with high and lower affinities, respectively. Using a cell apoptosis assay based on stimulation of a RANK:Fas fusion protein, and a cellular NF-κB signaling assay, we showed that R12-31 was agonist for both species. R12-31 interfered little or not at all with the binding of RANKL to RANK, in contrast to RANK-02 that efficiently prevented this interaction. Depending on the assay and species, RANK-02 was either a weak agonist or a partial antagonist of RANK. Both antibodies recognized human Langerhans cells, previously shown to express RANK, while dermal dendritic cells were poorly labeled. In vivo R12-31 agonist activity was demonstrated by its ability to induce the formation of intestinal villous microfold cells in mice. This characterization of two monoclonal antibodies should now allow better evaluation of their application as therapeutic reagents and investigative tools.