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.

Association of car ownership and physical activity across the spectrum of human development : Modeling the Epidemiologic Transition Study (METS).

BACKGROUND: Variations in physical activity (PA) across nations may be driven by socioeconomic position. As national incomes increase, car ownership becomes within reach of more individuals. This report characterizes associations between car ownership and PA in African-origin populations across 5 sites at different levels of economic development and with different transportation infrastructures: US, Seychelles, Jamaica, South Africa, and Ghana. METHODS: Twenty-five hundred adults, ages 25-45, were enrolled in the study. A total of 2,101 subjects had valid accelerometer-based PA measures (reported as average daily duration of moderate to vigorous PA, MVPA) and complete socioeconomic information. Our primary exposure of interest was whether the household owned a car. We adjusted for socioeconomic position using household income and ownership of common goods. RESULTS: Overall, PA levels did not vary largely between sites, with highest levels in South Africa, lowest in the US. Across all sites, greater PA was consistently associated with male gender, fewer years of education, manual occupations, lower income, and owning fewer material goods. We found heterogeneity across sites in car ownership: after adjustment for confounders, car owners in the US had 24.3 fewer minutes of MVPA compared to non-car owners in the US (20.7 vs. 45.1 minutes/day of MVPA); in the non-US sites, car-owners had an average of 9.7 fewer minutes of MVPA than non-car owners (24.9 vs. 34.6 minutes/day of MVPA). CONCLUSIONS: PA levels are similar across all study sites except Jamaica, despite very different levels of socioeconomic development. Not owning a car in the US is associated with especially high levels of MVPA. As car ownership becomes prevalent in the developing world, strategies to promote alternative forms of active transit may become important.

The dynamics of interfaces: Seven authors in search of encounters across levels of description of an event involving a mother, father, and baby

The interfaces between the intrapsychic, interactional, and intergenerational domains are a new frontier. As a pilot, we exposed ourselves to a complex but controllable situation as viewed by people whose main interest is in one of the three interfaces; we also fully integrated the subjects in the team, to learn about their subjective perspectives and to provide them with an enriching experience. We started with a brief "triadification" sequence (i.e., moving from a "two plus one" to a "three together" family organization). Considering this sequence as representing at a micro level many larger family transitions, we proceeded with a microanalytic interview, a psychodynamic investigation, and a family interview. As expected, larger patterns of correspondences are emerging. Central questions under debate are: What are the most appropriate units at each level of description and what are their articulations between these levels? What is the status of "triadification"? Les interfaces entre les domaines intrapsychiques, interactionnels et intergénérationnels représentent une nouvelle frontiére. A titre exploratoire, nous nous sommes exposés à une situation complexe mais contrǒlable ainsi que le voient ceux dont I'intérět principal se porte sur l'une de ces trois interfaces. Nous avons aussi entièrement intégré les sujets dans l'équipe, de facon à comprendre leur perspective subjective et à leur offrir une expérience enrichissante. Nous avons commencé avec une brève séquence de "triadification," c'est-à-dire passer d'une organisation familiale "deux plus un" à Ltne organisation familiale "trois (add sentenc)ensemble." Considérant cette séquence comme representative à un niveau microscopique de transitions familiales bien plus larges, nous avons procedé à l'entretien microanalytique, à une enquěte psychodynamique et à un entretien familial. Comme prévu, de grands patterns de correspondances émergent. Les questions essentielles sur lesquelles portent le débat sont: quelles les unités les plus appropiées à chaque niveau de description et quelles sont les articulations entre ces niveaux? Quel est le statut de la "triadification"?

Ab initio modeling and molecular dynamics simulation of the alpha 1b-adrenergic receptor activation.

This work describes the ab initio procedure employed to build an activation model for the alpha 1b-adrenergic receptor (alpha 1b-AR). The first version of the model was progressively modified and complicated by means of a many-step iterative procedure characterized by the employment of experimental validations of the model in each upgrading step. A combined simulated (molecular dynamics) and experimental mutagenesis approach was used to determine the structural and dynamic features characterizing the inactive and active states of alpha 1b-AR. The latest version of the model has been successfully challenged with respect to its ability to interpret and predict the functional properties of a large number of mutants. The iterative approach employed to describe alpha 1b-AR activation in terms of molecular structure and dynamics allows further complications of the model to allow prediction and interpretation of an ever-increasing number of experimental data.

Modeling sequencing errors by combining Hidden Markov models.

Among the largest resources for biological sequence data is the large amount of expressed sequence tags (ESTs) available in public and proprietary databases. ESTs provide information on transcripts but for technical reasons they often contain sequencing errors. Therefore, when analyzing EST sequences computationally, such errors must be taken into account. Earlier attempts to model error prone coding regions have shown good performance in detecting and predicting these while correcting sequencing errors using codon usage frequencies. In the research presented here, we improve the detection of translation start and stop sites by integrating a more complex mRNA model with codon usage bias based error correction into one hidden Markov model (HMM), thus generalizing this error correction approach to more complex HMMs. We show that our method maintains the performance in detecting coding sequences.

A protocol guided by transpulmonary thermodilution and lactate levels for resuscitation of patients with severe burns.

Over-resuscitation is deleterious in many critically ill conditions, including major burns. For more than 15 years, several strategies to reduce fluid administration in burns during the initial resuscitation phase have been proposed, but no single or simple parameter has shown superiority. Fluid administration guided by invasive hemodynamic parameters usually resulted in over-resuscitation. As reported in the previous issue of Critical Care, Sánchez-Sánchez and colleagues analyzed the performance of a 'permissive hypovolemia' protocol guided by invasive hemodynamic parameters (PiCCO, Pulsion Medical Systems, Munich, Germany) and vital signs in a prospective cohort over a 3-year period. The authors' results confirm that resuscitation can be achieved with below-normal levels of preload but at the price of a fluid administration greater than predicted by the Parkland formula (2 to 4 mL/kg per% burn). The classic approach based on an adapted Parkland equation may still be the simplest until further studies identify the optimal bundle of resuscitation goals.

D-amphetamine fails to increase extracellular dopamine levels in mice lacking alpha 1b-adrenergic receptors: relationship between functional and nonfunctional dopamine release.

It was found recently that locomotor and rewarding effects of psychostimulants and opiates were dramatically decreased or suppressed in mice lacking alpha1b-adrenergic receptors [alpha1b-adrenergic receptor knock-outs (alpha1bAR-KOs)] (Drouin et al., 2002). Here we show that blunted locomotor responses induced by 3 and 6 mg/kg d-amphetamine in alpha1bAR-KO mice [-84 and -74%, respectively, when compared with wild-type (WT) mice] are correlated with an absence of d-amphetamine-induced increase in extracellular dopamine (DA) levels in the nucleus accumbens of alpha1bAR-KO mice. Moreover, basal extracellular DA levels in the nucleus accumbens are lower in alpha1bAR-KO than in WT littermates (-28%; p < 0.001). In rats however, prazosin, an alpha1-adrenergic antagonist, decreases d-amphetamine-induced locomotor hyperactivity without affecting extracellular DA levels in the nucleus accumbens, a finding related to the presence of an important nonfunctional release of DA (Darracq et al., 1998). We show here that local d-amphetamine releases nonfunctional DA with the same affinity but a more than threefold lower amplitude in C57BL6/J mice than in Sprague Dawley rats. Altogether, this suggests that a trans-synaptic mechanism amplifies functional DA into nonfunctional DA release. Our data confirm the presence of a powerful coupling between noradrenergic and dopaminergic neurons through the stimulation of alpha1b-adrenergic receptors and indicate that nonfunctional DA release is critical in the interpretation of changes in extracellular DA levels. These results suggest that alpha1b-adrenergic receptors may be important therapeutic pharmacological targets not only in addiction but also in psychosis because most neuroleptics possess anti-alpha1-adrenergic properties.

HCV infection after liver transplantation is associated with lower levels of alloreactive CD4+CD25+CD45RO+IL-7R+high+ T cells

Aim: Expression of IL-7R discriminates alloreactive CD4 T cells (Foxp3 negative), from IL-7Rlow regulatory CD4 T cells (Foxp3 positive). Chronic hepatitis C virus infection (HCV) reduces expression of IL-7R on T cells thus promoting persistence of infection. The aim of this study was to analyze the effect of HCV infection on the expression of IL-7R of activated CD4+ T cells in liver transplant patients. Patients and methods: We analyzed PBMC from liver transplant recipients for the expression of CD4, CD25, FoxP3, IL-7R (24 HCV negative and 29 HCV-chronically infected). We compared these data with non-transplanted individuals (52 HCV-chronically infected patients and 38 healthy donors). Results: In HCV-infected liver transplant recipients, levels of CD4+CD25+CD45RO+IL-7R+ T cells were significantly reduced (10.5+/-0.9%) when compared to non-HCV-infected liver transplant recipients (17.6+/-1.4%) (P<0.001), while both groups (HCV-infected and negative transplant recipients) had significantly higher levels than healthy individuals (6.6+/-0.9%) (P<0.0001). After successful antiviral therapy (sustained antiviral response), 6 HCV-infected transplant recipients showed an increase of CD4+CD25+CD45RO+IL-7R+ T cells, reaching levels similar to that of non-HCVinfected recipients (10.73+/-2.63% prior therapy versus 21.7+/-6.3% after clearance of HCV). (P<0.05) In 4 non-responders (i.e. HCVRNA remaining present in serum), levels of CD4+CD25+CD45RO+IL-7R+ T cells remained unmodified during and after antiviral treatment (11.8+/- 3.3% versus 11.3+/-3.3% respectively). Conclusions: Overall, these data indicate that CD4+CD25+CD45RO+IL-7R+ T cells appear to be modulated by chronic HCV infection after liver transplantation. Whether lower levels of alloreactive T cells in HCV-infected liver transplant recipients are associated with a tolerogenic profile remains to be studied.

Relationships between imatinib plasma levels and efficacy

Purpose: While imatinib has revolutionized the treatment of chronic myeloid leukaemia (CML) and gastrointestinal stromal tumors (GIST), its pharmacokinetic-pharmacodynamic relationships have been poorly studied. This study aimed to explore the issue in oncologic patients, and to evaluate the specific influence of the target genotype in a GIST subpopulation. Patients and methods: Data from 59 patients (321 plasma samples) were collected during a previous pharmacokinetic study. Based on a population model purposely developed, individual post-hoc Bayesian estimates of pharmacokinetic parameters were derived, and used to estimate drug exposure (AUC; area under curve). Free fraction parameters were deduced from a model incorporating plasma alpha1-acid glycoprotein levels. Associations between AUC (or clearance) and therapeutic response (coded on a 3-point scale), or tolerability (4-point scale), were explored by ordered logistic regression. Influence of KIT genotype on response was also assessed in GIST patients. Results: Total and free drug exposure correlated with the number of side effects (p < 0.005). A relationship with response was not evident in the whole patient set (with good-responders tending to receive lower doses and bad-responders higher doses). In GIST patients however, higher free drug exposure predicted better responses. A strong association was notably observed in patients harboring an exon 9 mutation or a wild type KIT, known to decrease tumor sensitivity towards imatinib (p < 0.005). Conclusions: Our results are arguments to further evaluate the potential benefit of a therapeutic monitoring program for imatinib. Our data also suggest that stratification by genotype will be important in future trials.

SUCLG2 identified as both a determinator of CSF Aβ1-42 levels and an attenuator of cognitive decline in Alzheimer's disease.

Cerebrospinal fluid amyloid-beta 1-42 (Aβ1-42) and phosphorylated Tau at position 181 (pTau181) are biomarkers of Alzheimer's disease (AD). We performed an analysis and meta-analysis of genome-wide association study data on Aβ1-42 and pTau181 in AD dementia patients followed by independent replication. An association was found between Aβ1-42 level and a single-nucleotide polymorphism in SUCLG2 (rs62256378) (P = 2.5×10(-12)). An interaction between APOE genotype and rs62256378 was detected (P = 9.5 × 10(-5)), with the strongest effect being observed in APOE-ε4 noncarriers. Clinically, rs62256378 was associated with rate of cognitive decline in AD dementia patients (P = 3.1 × 10(-3)). Functional microglia experiments showed that SUCLG2 was involved in clearance of Aβ1-42.

An alteration in the levels of populations of CD4+ Treg is in part responsible for altered cytokine production by cells of aged mice which follows injection with a fetal liver extract.

We have shown previously that a fetal sheep liver extract (FSLE) containing significant quantities of fetal ovine gamma globin chain (Hbgamma) and LPS injected into aged (>20 months) mice could reverse the altered polarization (increased IL-4 and IL-10 with decreased IL-2 and IFNgamma) in cytokine production seen from ConA stimulated lymphoid cells of those mice. The mechanism(s) behind this change in cytokine production were not previously investigated. We report below that aged mice show a >60% decline in numbers and suppressive function of both CD4(+)CD25(+)Foxp3(+) Treg and so-called Tr3 (CD4(+)TGFbeta(+)), and that their number/function is restored to levels seen in control (8-week-old) mice by FSLE. In addition, on a per cell basis, CD4(+)CD25(-)Treg from aged mice were >4-fold more effective in suppression of proliferation and IL-2 production from ConA-activated lymphoid cells of a pool of CD4(+)CD25(-)T cells from 8-week-old mice than similar cells from young animals, and this suppression by CD25(-)T cells was also ameliorated following FSLE treatment. Infusion of anti-TGFbeta and anti-IL-10 antibodies in vivo altered Treg development following FSLE treatment, and attenuated FSLE-induced alterations in cytokine production profiles.

Association between conventional risk factors and different ultrasound-based markers of atherosclerosis at carotid and femoral levels in a middle-aged population.

Ultrasound detection of sub-clinical atherosclerosis (ATS) may help identify individuals at high cardiovascular risk. Most studies evaluated intima-media thickness (IMT) at carotid level. We compared the relationships between main cardiovascular risk factors (CVRF) and five indicators of ATS (IMT, mean and maximal plaque thickness, mean and maximal plaque area) at both carotid and femoral levels. Ultrasound was performed on 496 participants aged 45-64 years randomly selected from the general population of the Republic of Seychelles. 73.4 % participants had ≥ 1 plaque (IMT thickening ≥ 1.2 mm) at carotid level and 67.5 % at femoral level. Variance (adjusted R2) contributed by age, sex and CVRF (smoking, LDL-cholesterol, HDL-cholesterol, blood pressure, diabetes) in predicting any of the ATS markers was larger at femoral than carotid level. At both carotid and femoral levels, the association between CVRF and ATS was stronger based on plaque-based markers than IMT. Our findings show that the associations between CVRF and ATS markers were stronger at femoral than carotid level, and with plaque-based markers rather than IMT. Pending comparison of these markers using harder cardiovascular endpoints, our findings suggest that markers based on plaque morphology assessed at femoral artery level might be useful cardiovascular risk predictors.

Body mass index, cigarette smoking, and alcohol consumption and cancers of the oral cavity, pharynx, and larynx: modeling odds ratios in pooled case-control data.

Odds ratios for head and neck cancer increase with greater cigarette and alcohol use and lower body mass index (BMI; weight (kg)/height(2) (m(2))). Using data from the International Head and Neck Cancer Epidemiology Consortium, the authors conducted a formal analysis of BMI as a modifier of smoking- and alcohol-related effects. Analysis of never and current smokers included 6,333 cases, while analysis of never drinkers and consumers of < or =10 drinks/day included 8,452 cases. There were 8,000 or more controls, depending on the analysis. Odds ratios for all sites increased with lower BMI, greater smoking, and greater drinking. In polytomous regression, odds ratios for BMI (P = 0.65), smoking (P = 0.52), and drinking (P = 0.73) were homogeneous for oral cavity and pharyngeal cancers. Odds ratios for BMI and drinking were greater for oral cavity/pharyngeal cancer (P < 0.01), while smoking odds ratios were greater for laryngeal cancer (P < 0.01). Lower BMI enhanced smoking- and drinking-related odds ratios for oral cavity/pharyngeal cancer (P < 0.01), while BMI did not modify smoking and drinking odds ratios for laryngeal cancer. The increased odds ratios for all sites with low BMI may suggest related carcinogenic mechanisms; however, BMI modification of smoking and drinking odds ratios for cancer of the oral cavity/pharynx but not larynx cancer suggests additional factors specific to oral cavity/pharynx cancer.