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.

Biological impact assessment of nanomaterial used in nanomedicine. Introduction to the NanoTEST project.

Therapeutic nanoparticles (NPs) are used in nanomedicine as drug carriers or imaging agents, providing increased selectivity/specificity for diseased tissues. The first NPs in nanomedicine were developed for increasing the efficacy of known drugs displaying dose-limiting toxicity and poor bioavailability and for enhancing disease detection. Nanotechnologies have gained much interest owing to their huge potential for applications in industry and medicine. It is necessary to ensure and control the biocompatibility of the components of therapeutic NPs to guarantee that intrinsic toxicity does not overtake the benefits. In addition to monitoring their toxicity in vitro, in vivo and in silico, it is also necessary to understand their distribution in the human body, their biodegradation and excretion routes and dispersion in the environment. Therefore, a deep understanding of their interactions with living tissues and of their possible effects in the human (and animal) body is required for the safe use of nanoparticulate formulations. Obtaining this information was the main aim of the NanoTEST project, and the goals of the reports collected together in this special issue are to summarise the observations and results obtained by the participating research teams and to provide methodological tools for evaluating the biological impact of NPs.

Comparison of Immunogenicity in Rhesus Macaques of Transmitted-Founder, HIV-1 Group M Consensus, and Trivalent Mosaic Envelope Vaccines Formulated as a DNA Prime, NYVAC, and Envelope Protein Boost.

An effective human immunodeficiency virus type 1 (HIV-1) vaccine must induce protective antibody responses, as well as CD4(+) and CD8(+) T cell responses, that can be effective despite extraordinary diversity of HIV-1. The consensus and mosaic immunogens are complete but artificial proteins, computationally designed to elicit immune responses with improved cross-reactive breadth, to attempt to overcome the challenge of global HIV diversity. In this study, we have compared the immunogenicity of a transmitted-founder (T/F) B clade Env (B.1059), a global group M consensus Env (Con-S), and a global trivalent mosaic Env protein in rhesus macaques. These antigens were delivered using a DNA prime-recombinant NYVAC (rNYVAC) vector and Env protein boost vaccination strategy. While Con-S Env was a single sequence, mosaic immunogens were a set of three Envs optimized to include the most common forms of potential T cell epitopes. Both Con-S and mosaic sequences retained common amino acids encompassed by both antibody and T cell epitopes and were central to globally circulating strains. Mosaics and Con-S Envs expressed as full-length proteins bound well to a number of neutralizing antibodies with discontinuous epitopes. Also, both consensus and mosaic immunogens induced significantly higher gamma interferon (IFN-γ) enzyme-linked immunosorbent spot assay (ELISpot) responses than B.1059 immunogen. Immunization with these proteins, particularly Con-S, also induced significantly higher neutralizing antibodies to viruses than B.1059 Env, primarily to tier 1 viruses. Both Con-S and mosaics stimulated more potent CD8-T cell responses against heterologous Envs than did B.1059. Both antibody and cellular data from this study strengthen the concept of using in silico-designed centralized immunogens for global HIV-1 vaccine development strategies. IMPORTANCE: There is an increasing appreciation for the importance of vaccine-induced anti-Env antibody responses for preventing HIV-1 acquisition. This nonhuman primate study demonstrates that in silico-designed global HIV-1 immunogens, designed for a human clinical trial, are capable of eliciting not only T lymphocyte responses but also potent anti-Env antibody responses.

Comparative genomics suggests primary homothallism of Pneumocystis species.

UNLABELLED: Pneumocystis species are fungal parasites of mammal lungs showing host specificity. Pneumocystis jirovecii colonizes humans and causes severe pneumonia in immunosuppressed individuals. In the absence of in vitro cultures, the life cycle of these fungi remains poorly known. Sexual reproduction probably occurs, but the system of this process and the mating type (MAT) genes involved are not characterized. In the present study, we used comparative genomics to investigate the issue in P. jirovecii and Pneumocystis carinii, the species infecting rats, as well as in their relative Taphrina deformans. We searched sex-related genes using 103 sequences from the relative Schizosaccharomyces pombe as queries. Genes homologous to several sex-related role categories were identified in all species investigated, further supporting sexuality in these organisms. Extensive in silico searches identified only three putative MAT genes in each species investigated (matMc, matMi, and matPi). In P. jirovecii, these genes clustered on the same contig, proving their contiguity in the genome. This organization seems compatible neither with heterothallism, because two different MAT loci on separate DNA molecules would have been detected, nor with secondary homothallism, because the latter involves generally more MAT genes. Consistently, we did not detect cis-acting sequences for mating type switching in secondary homothallism, and PCR revealed identical MAT genes in P. jirovecii isolates from six patients. A strong synteny of the genomic region surrounding the putative MAT genes exists between the two Pneumocystis species. Our results suggest the hypothesis that primary homothallism is the system of reproduction of Pneumocystis species and T. deformans. IMPORTANCE: Sexual reproduction among fungi can involve a single partner (homothallism) or two compatible partners (heterothallism). We investigated the issue in three pathogenic fungal relatives: Pneumocystis jirovecii, which causes severe pneumonia in immunocompromised humans; Pneumocystis carinii, which infects rats; and the plant pathogen Taphrina deformans. The nature, the number, and the organization within the genome of the genes involved in sexual reproduction were determined. The three species appeared to harbor a single genomic region gathering only three genes involved in sexual differentiation, an organization which is compatible with sexual reproduction involving a single partner. These findings illuminate the strategy adopted by fungal pathogens to infect their hosts.

Predicting drug metabolism: experiment and/or computation?

Drug metabolism can produce metabolites with physicochemical and pharmacological properties that differ substantially from those of the parent drug, and consequently has important implications for both drug safety and efficacy. To reduce the risk of costly clinical-stage attrition due to the metabolic characteristics of drug candidates, there is a need for efficient and reliable ways to predict drug metabolism in vitro, in silico and in vivo. In this Perspective, we provide an overview of the state of the art of experimental and computational approaches for investigating drug metabolism. We highlight the scope and limitations of these methods, and indicate strategies to harvest the synergies that result from combining measurement and prediction of drug metabolism.

Angiogenic activity of breast cancer patients' monocytes reverted by combined use of systems modeling and experimental approaches.

Angiogenesis plays a key role in tumor growth and cancer progression. TIE-2-expressing monocytes (TEM) have been reported to critically account for tumor vascularization and growth in mouse tumor experimental models, but the molecular basis of their pro-angiogenic activity are largely unknown. Moreover, differences in the pro-angiogenic activity between blood circulating and tumor infiltrated TEM in human patients has not been established to date, hindering the identification of specific targets for therapeutic intervention. In this work, we investigated these differences and the phenotypic reversal of breast tumor pro-angiogenic TEM to a weak pro-angiogenic phenotype by combining Boolean modelling and experimental approaches. Firstly, we show that in breast cancer patients the pro-angiogenic activity of TEM increased drastically from blood to tumor, suggesting that the tumor microenvironment shapes the highly pro-angiogenic phenotype of TEM. Secondly, we predicted in silico all minimal perturbations transitioning the highly pro-angiogenic phenotype of tumor TEM to the weak pro-angiogenic phenotype of blood TEM and vice versa. In silico predicted perturbations were validated experimentally using patient TEM. In addition, gene expression profiling of TEM transitioned to a weak pro-angiogenic phenotype confirmed that TEM are plastic cells and can be reverted to immunological potent monocytes. Finally, the relapse-free survival analysis showed a statistically significant difference between patients with tumors with high and low expression values for genes encoding transitioning proteins detected in silico and validated on patient TEM. In conclusion, the inferred TEM regulatory network accurately captured experimental TEM behavior and highlighted crosstalk between specific angiogenic and inflammatory signaling pathways of outstanding importance to control their pro-angiogenic activity. Results showed the successful in vitro reversion of such an activity by perturbation of in silico predicted target genes in tumor derived TEM, and indicated that targeting tumor TEM plasticity may constitute a novel valid therapeutic strategy in breast cancer.

The SwissLipids knowledgebase for lipid biology.

MOTIVATION: Lipids are a large and diverse group of biological molecules with roles in membrane formation, energy storage and signaling. Cellular lipidomes may contain tens of thousands of structures, a staggering degree of complexity whose significance is not yet fully understood. High-throughput mass spectrometry-based platforms provide a means to study this complexity, but the interpretation of lipidomic data and its integration with prior knowledge of lipid biology suffers from a lack of appropriate tools to manage the data and extract knowledge from it. RESULTS: To facilitate the description and exploration of lipidomic data and its integration with prior biological knowledge, we have developed a knowledge resource for lipids and their biology-SwissLipids. SwissLipids provides curated knowledge of lipid structures and metabolism which is used to generate an in silico library of feasible lipid structures. These are arranged in a hierarchical classification that links mass spectrometry analytical outputs to all possible lipid structures, metabolic reactions and enzymes. SwissLipids provides a reference namespace for lipidomic data publication, data exploration and hypothesis generation. The current version of SwissLipids includes over 244 000 known and theoretically possible lipid structures, over 800 proteins, and curated links to published knowledge from over 620 peer-reviewed publications. We are continually updating the SwissLipids hierarchy with new lipid categories and new expert curated knowledge. AVAILABILITY: SwissLipids is freely available at http://www.swisslipids.org/. CONTACT: alan.bridge@isb-sib.ch SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Congenital hypogonadotropic hypogonadism with split hand/foot malformation: a clinical entity with a high frequency of FGFR1 mutations.

PURPOSE: Congenital hypogonadotropic hypogonadism (CHH) and split hand/foot malformation (SHFM) are two rare genetic conditions. Here we report a clinical entity comprising the two. METHODS: We identified patients with CHH and SHFM through international collaboration. Probands and available family members underwent phenotyping and screening for FGFR1 mutations. The impact of identified mutations was assessed by sequence- and structure-based predictions and/or functional assays. RESULTS: We identified eight probands with CHH with (n = 3; Kallmann syndrome) or without anosmia (n = 5) and SHFM, seven of whom (88%) harbor FGFR1 mutations. Of these seven, one individual is homozygous for p.V429E and six individuals are heterozygous for p.G348R, p.G485R, p.Q594*, p.E670A, p.V688L, or p.L712P. All mutations were predicted by in silico analysis to cause loss of function. Probands with FGFR1 mutations have severe gonadotropin-releasing hormone deficiency (absent puberty and/or cryptorchidism and/or micropenis). SHFM in both hands and feet was observed only in the patient with the homozygous p.V429E mutation; V429 maps to the fibroblast growth factor receptor substrate 2α binding domain of FGFR1, and functional studies of the p.V429E mutation demonstrated that it decreased recruitment and phosphorylation of fibroblast growth factor receptor substrate 2α to FGFR1, thereby resulting in reduced mitogen-activated protein kinase signaling. CONCLUSION: FGFR1 should be prioritized for genetic testing in patients with CHH and SHFM because the likelihood of a mutation increases from 10% in the general CHH population to 88% in these patients.Genet Med 17 8, 651-659.