Dynamical hybrid modeling of human metabolism

Human metabolism plays a key role in disease pathogenesis and drug action. Half a century of biochemical literature leveraged by the advent of genomics allowed the emergence of computational modeling techniques and the in silico analysis of complex biological systems. In particular, Constraint-Based Reconstruction and Analysis (COBRA) methods address the complexity of metabolism through building tissue-specific networks in their steady state. It is known that biological systems respond to perturbations induced by pathogens, drugs or malignant processes by shifting their activity to safeguard key metabolic functions. Extending the modeling framework to consider the dynamics of these complex systems will bring simulations closer to observable human phenotypes. In this thesis, I combined physiologically-based pharmacokinetic (PBPK) models with genome-scale metabolic models (GSMMs) to form hybrid genome-scale dynamical models that provide a hypothesis-free framework to study the perturbations induced by one or more perturbagen on human tissues. On a first stage, these methodologies were applied to decipher the absorption of levodopa and amino acids by the intestinal epithelium and allowed to derive a model-based diet for Parkinson's Disease patients. In the next phase, we extended the study to 605 drugs in order to predict the occurrence of gastrointestinal side effects through a machine learning classifier, using a combination of gene expression and metabolic reactions set as features. Finally, the approach upscaled to several tissues, specifically to investigate the genesis of metabolic symptoms in type 1 diabetes and to suggest key metabolic players underlying within and between-individual variability to insulin action. Taken as whole, the integration of two modeling techniques constrained by expert biological knowledge and heterogeneous data types will be a step forward in achieving convergence in human biology

DRAGON: Determining Regulatory Associations using Graphical models on multi-Omic Networks

The increasing quantity of multi-omic data, such as methylomic and transcriptomic profiles collected on the same specimen or even on the same cell, provides a unique opportunity to explore the complex interactions that define cell phenotype and govern cellular responses to perturbations. We propose a network approach based on Gaussian Graphical Models (GGMs) that facilitates the joint analysis of paired omics data. This method, called DRAGON (Determining Regulatory Associations using Graphical models on multi-Omic Networks), calibrates its parameters to achieve an optimal trade-off between the network’s complexity and estimation accuracy, while explicitly accounting for the characteristics of each of the assessed omics ‘layers.’ In simulation studies, we show that DRAGON adapts to edge density and feature size differences between omics layers, improving model inference and edge recovery compared to state-of-the-art methods. We further demonstrate in an analysis of joint transcriptome - methylome data from TCGA breast cancer specimens that DRAGON can identify key molecular mechanisms such as gene regulation via promoter methylation. In particular, we identify Transcription Factor AP-2 Beta (TFAP2B) as a potential multi-omic biomarker for basal-type breast cancer. DRAGON is available as open-source code in Python through the Network Zoo package (netZooPy v0.8; netzoo.github.io)

Toward a Generalizable Framework of Disturbance Ecology Through Crowdsourced Science

© 2021 Graham, Averill, Bond-Lamberty, Knelman, Krause, Peralta, Shade, Smith, Cheng, Fanin, Freund, Garcia, Gibbons, Van Goethem, Guebila, Kemppinen, Nowicki, Pausas, Reed, Rocca, Sengupta, Sihi, Simonin, Słowiński, Spawn, Sutherland, Tonkin, Wisnoski, Zipper and Contributor Consortium.Disturbances fundamentally alter ecosystem functions, yet predicting their impacts remains a key scientific challenge. While the study of disturbances is ubiquitous across many ecological disciplines, there is no agreed-upon, cross-disciplinary foundation for discussing or quantifying the complexity of disturbances, and no consistent terminology or methodologies exist. This inconsistency presents an increasingly urgent challenge due to accelerating global change and the threat of interacting disturbances that can destabilize ecosystem responses. By harvesting the expertise of an interdisciplinary cohort of contributors spanning 42 institutions across 15 countries, we identified an essential limitation in disturbance ecology: the word ‘disturbance’ is used interchangeably to refer to both the events that cause, and the consequences of, ecological change, despite fundamental distinctions between the two meanings. In response, we developed a generalizable framework of ecosystem disturbances, providing a well-defined lexicon for understanding disturbances across perspectives and scales. The framework results from ideas that resonate across multiple scientific disciplines and provides a baseline standard to compare disturbances across fields. This framework can be supplemented by discipline-specific variables to provide maximum benefit to both inter- and intra-disciplinary research. To support future syntheses and meta-analyses of disturbance research, we also encourage researchers to be explicit in how they define disturbance drivers and impacts, and we recommend minimum reporting standards that are applicable regardless of scale. Finally, we discuss the primary factors we considered when developing a baseline framework and propose four future directions to advance our interdisciplinary understanding of disturbances and their social-ecological impacts: integrating across ecological scales, understanding disturbance interactions, establishing baselines and trajectories, and developing process-based models and ecological forecasting initiatives. Our experience through this process motivates us to encourage the wider scientific community to continue to explore new approaches for leveraging Open Science principles in generating creative and multidisciplinary ideas.This research was supported by the U.S. Department of Energy (DOE), Office of Biological and Environmental Research (BER), as part of Subsurface Biogeochemical Research Program’s Scientific Focus Area (SFA) at the Pacific Northwest National Laboratory (PNNL). PNNL is operated for DOE by Battelle under contract DE-AC06-76RLO 1830

Global network of computational biology communities: ISCB's regional student groups breaking barriers [version 1; peer review: Not peer reviewed]

Regional Student Groups (RSGs) of the International Society for Computational Biology Student Council (ISCB-SC) have been instrumental to connect computational biologists globally and to create more awareness about bioinformatics education. This article highlights the initiatives carried out by the RSGs both nationally and internationally to strengthen the present and future of the bioinformatics community. Moreover, we discuss the future directions the organization will take and the challenges to advance further in the ISCB-SC main mission: “Nurture the new generation of computational biologists”.Fil: Shome, Sayane. University of Iowa; Estados UnidosFil: Parra, Rodrigo Gonzalo. European Molecular Biology Laboratory; Alemania. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Fatima, Nazeefa. Uppsala Universitet; SueciaFil: Monzon, Alexander Miguel. Università di Padova; ItaliaFil: Cuypers, Bart. Universiteit Antwerp; BélgicaFil: Moosa, Yumna. University of KwaZulu Natal; SudáfricaFil: Da Rocha Coimbra, Nilson. Universidade Federal de Minas Gerais; BrasilFil: Assis, Juliana. Universidade Federal de Minas Gerais; BrasilFil: Giner Delgado, Carla. Universitat Autònoma de Barcelona; EspañaFil: Dönertaş, Handan Melike. European Molecular Biology Laboratory. European Bioinformatics Institute; Reino UnidoFil: Cuesta Astroz, Yesid. Universidad de Antioquia; Colombia. Universidad Ces. Facultad de Medicina.; ColombiaFil: Saarunya, Geetha. University of South Carolina; Estados UnidosFil: Allali, Imane. Universite Mohammed V. Rabat; Otros paises de África. University of Cape Town; SudáfricaFil: Gupta, Shruti. Jawaharlal Nehru University; IndiaFil: Srivastava, Ambuj. Indian Institute of Technology Madras; IndiaFil: Kalsan, Manisha. Jawaharlal Nehru University; IndiaFil: Valdivia, Catalina. Universidad Andrés Bello; ChileFil: Olguín Orellana, Gabriel José. Universidad de Talca; ChileFil: Papadimitriou, Sofia. Vrije Unviversiteit Brussel; Bélgica. Université Libre de Bruxelles; BélgicaFil: Parisi, Daniele. Katholikie Universiteit Leuven; BélgicaFil: Kristensen, Nikolaj Pagh. Technical University of Denmark; DinamarcaFil: Rib, Leonor. Universidad de Copenhagen; DinamarcaFil: Guebila, Marouen Ben. University of Luxembourg; LuxemburgoFil: Bauer, Eugen. University of Luxembourg; LuxemburgoFil: Zaffaroni, Gaia. University of Luxembourg; LuxemburgoFil: Bekkar, Amel. Universite de Lausanne; SuizaFil: Ashano, Efejiro. APIN Public Health Initiatives; NigeriaFil: Paladin, Lisanna. Università di Padova; ItaliaFil: Necci, Marco. Università di Padova; ItaliaFil: Moreyra, Nicolás Nahuel. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Ecología, Genética y Evolución de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Ecología, Genética y Evolución de Buenos Aires; Argentin

The Network Zoo: a multilingual package for the inference and analysis of gene regulatory networks

Inference and analysis of gene regulatory networks (GRNs) require software that integrates multi-omic data from various sources. The Network Zoo (netZoo; netzoo.github.io) is a collection of open-source methods to infer GRNs, conduct differential network analyses, estimate community structure, and explore the transitions between biological states. The netZoo builds on our ongoing development of network methods, harmonizing the implementations in various computing languages and between methods to allow better integration of these tools into analytical pipelines. We demonstrate the utility using multi-omic data from the Cancer Cell Line Encyclopedia. We will continue to expand the netZoo to incorporate additional methods

Crowdsourced mapping of unexplored target space of kinase inhibitors

Despite decades of intensive search for compounds that modulate the activity of particular protein targets, a large proportion of the human kinome remains as yet undrugged. Effective approaches are therefore required to map the massive space of unexplored compound-kinase interactions for novel and potent activities. Here, we carry out a crowdsourced benchmarking of predictive algorithms for kinase inhibitor potencies across multiple kinase families tested on unpublished bioactivity data. We find the top-performing predictions are based on various models, including kernel learning, gradient boosting and deep learning, and their ensemble leads to a predictive accuracy exceeding that of single-dose kinase activity assays. We design experiments based on the model predictions and identify unexpected activities even for under-studied kinases, thereby accelerating experimental mapping efforts. The open-source prediction algorithms together with the bioactivities between 95 compounds and 295 kinases provide a resource for benchmarking prediction algorithms and for extending the druggable kinome. The IDG-DREAM Challenge carried out crowdsourced benchmarking of predictive algorithms for kinase inhibitor activities on unpublished data. This study provides a resource to compare emerging algorithms and prioritize new kinase activities to accelerate drug discovery and repurposing efforts

Predicting gastrointestinal drug effects using contextualized metabolic models.

Gastrointestinal side effects are among the most common classes of adverse reactions associated with orally absorbed drugs. These effects decrease patient compliance with the treatment and induce undesirable physiological effects. The prediction of drug action on the gut wall based on in vitro data solely can improve the safety of marketed drugs and first-in-human trials of new chemical entities. We used publicly available data of drug-induced gene expression changes to build drug-specific small intestine epithelial cell metabolic models. The combination of measured in vitro gene expression and in silico predicted metabolic rates in the gut wall was used as features for a multilabel support vector machine to predict the occurrence of side effects. We showed that combining local gut wall-specific metabolism with gene expression performs better than gene expression alone, which indicates the role of small intestine metabolism in the development of adverse reactions. Furthermore, we reclassified FDA-labeled drugs with respect to their genetic and metabolic profiles to show hidden similarities between seemingly different drugs. The linkage of xenobiotics to their transcriptomic and metabolic profiles could take pharmacology far beyond the usual indication-based classifications

Un modèle d'optimisation nutritionelle pour les Parkinsoniens sous levodopa.

Levodopa has been the gold standard for Parkinson’s disease treatment for more than 40 years. Its bioavailability is hindered by dietary amino acids, leading to fluctuations in the motor response particularly in late-stage (stage 3 and 4 on Hoehn and Yahr scale) patients. The routine dietary intervention consists of low-protein (<0.8 g/kg) diets or the redistribution of daily protein allowance to the last meal. Computational modeling was used to examine the fluctuation of gastrointestinal levodopa absorption under consideration of the diet by (i) identifying the group of patients that could benefit from dietary interventions, (ii) comparing existing diet recommendations for their impact on levodopa bioavailability, and (iii) suggesting a mechanism-based dietary intervention. We developed a multiscale computational model consisting of an ordinary differential equations-based advanced compartmentalized absorption and transit (ACAT) gut model and metabolic genome-scale small intestine epithelial cell model. We used this model to investigate complex spatiotemporal relationship between dietary amino acids and levodopa absorption. Our model predicted an improvement in bioavailability, as reflected by blood concentrations of levodopa with protein redistribution diet by 34% compared with a low-protein diet and by 11% compared with the ante cibum (a.c.) administration. These results are consistent with the reported better outcome in late-stage patients. A systematic analysis of the effect of different amino acids in the diet suggested that a serine-rich diet could improve the bioavailability by 22% compared with the a.c. administration. In addition, the slower gastric emptying rate in PD patients exacerbates the loss of levodopa due to competition. Optimizing dietary recommendations in quantity, composition, and intake time holds the promise to improve levodopa efficiency and patient’s quality of life based on highly detailed, mechanistic models of gut physiology endowed with improved extrapolative properties, thus paving the way for precision medical treatment

Regulatory network of PD1 signaling is associated with prognosis in glioblastoma multiforme

Glioblastoma is an aggressive cancer of the brain and spine. While analysis of glioblastoma ‘omics data has somewhat improved our understanding of the disease, it has not led to direct improvement in patient survival. Cancer survival is often characterized by differences in gene expression, but the mechanisms that drive these differences are generally unknown. We therefore set out to model the regulatory mechanisms associated with glioblastoma survival. We inferred individual patient gene regulatory networks using data from two different expression platforms from The Cancer Genome Atlas. We performed comparative network analysis between patients with long- and short-term survival. Seven pathways were identified as associated with survival, all of them involved in immune signaling; differential regulation of PD1 signaling was validated to correspond with outcome in an independent dataset from the German Glioma Network. In this pathway, transcriptional repression of genes for which treatment options are available was lost in short-term survivors; this was independent of mutational burden and only weakly associated with T-cell infiltration. Collectively, these results provide a new way to stratify patients with glioblastoma that uses network features as biomarkers to predict survival. They also identify new potential therapeutic interventions, underscoring the value of analyzing gene regulatory networks in individual patients with cancer