COVID-19 and cholesterol biosynthesis: Towards innovative decision support systems.

With COVID-19 becoming endemic, there is a continuing need to find biomarkers characterizing the disease and aiding in patient stratification. We studied the relation between COVID-19 and cholesterol biosynthesis by comparing 10 intermediates of cholesterol biosynthesis during the hospitalization of 164 patients (admission, disease deterioration, discharge) admitted to the University Medical Center of Ljubljana. The concentrations of zymosterol, 24-dehydrolathosterol, desmosterol, and zymostenol were significantly altered in COVID-19 patients. We further developed a predictive model for disease severity based on clinical parameters alone and their combination with a subset of sterols. Our machine learning models applying 8 clinical parameters predicted disease severity with excellent accuracy (AUC = 0.96), showing substantial improvement over current clinical risk scores. After including sterols, model performance remained better than COVID-GRAM. This is the first study to examine cholesterol biosynthesis during COVID-19 and shows that a subset of cholesterol-related sterols is associated with the severity of COVID-19.

Integrative computational modeling to unravel novel potential biomarkers in hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is a major health problem around the world. The management of this disease is complicated by the lack of noninvasive diagnostic tools and the few treatment options available. Better clinical outcomes can be achieved if HCC is detected early, but unfortunately, clinical signs appear when the disease is in its late stages. We aim to identify novel genes that can be targeted for the diagnosis and therapy of HCC. We performed a meta-analysis of transcriptomics data to identify differentially expressed genes and applied network analysis to identify hub genes. Fatty acid metabolism, complement and coagulation cascade, chemical carcinogenesis and retinol metabolism were identified as key pathways in HCC. Furthermore, we integrated transcriptomics data into a reference human genome-scale metabolic model to identify key reactions and subsystems relevant in HCC. We conclude that fatty acid activation, purine metabolism, vitamin D, and E metabolism are key processes in the development of HCC and therefore need to be further explored for the development of new therapies. We provide the first evidence that GABRP, HBG1 and DAK (TKFC) genes are important in HCC in humans and warrant further studies.

A computational design of a programmable biological processor.

Basic synthetic information processing structures, such as logic gates, oscillators and flip-flops, have already been implemented in living organisms. Current implementations of these structures have yet to be extended to more complex processing structures that would constitute a biological computer. We make a step forward towards the construction of a biological computer. We describe a model-based computational design of a biological processor that uses transcription and translation resources of the host cell to perform its operations. The proposed processor is composed of an instruction memory containing a biological program, a program counter that is used to address this memory, and a biological oscillator that triggers the execution of the next instruction in the memory. We additionally describe the implementation of a biological compiler that compiles a sequence of human-readable instructions into ordinary differential equation-based models, which can be used to simulate and analyse the dynamics of the processor. The proposed implementation presents the first programmable biological processor that exploits cellular resources to execute the specified instructions. We demonstrate the application of the described processor on a set of simple yet scalable biological programs. Biological descriptions of these programs can be produced manually or automatically using the provided compiler.

Review and assessment of Boolean approaches for inference of gene regulatory networks.

Boolean descriptions of gene regulatory networks can provide an insight into interactions between genes. Boolean networks hold predictive power, are easy to understand, and can be used to simulate the observed networks in different scenarios. We review fundamental and state-of-the-art methods for inference of Boolean networks. We introduce a methodology for a straightforward evaluation of Boolean inference approaches based on the generation of evaluation datasets, application of selected inference methods, and evaluation of performance measures to guide the selection of the best method for a given inference problem. We demonstrate this procedure on inference methods REVEAL (REVerse Engineering ALgorithm), Best-Fit Extension, MIBNI (Mutual Information-based Boolean Network Inference), GABNI (Genetic Algorithm-based Boolean Network Inference) and ATEN (AND/OR Tree ENsemble algorithm), which infers Boolean descriptions of gene regulatory networks from discretised time series data. Boolean inference approaches tend to perform better in terms of dynamic accuracy, and slightly worse in terms of structural correctness. We believe that the proposed methodology and provided guidelines will help researchers to develop Boolean inference approaches with a good predictive capability while maintaining structural correctness and biological relevance.

Integration of omics data to generate and analyse COVID-19 specific genome-scale metabolic models.

COVID-19 presents a complex disease that needs to be addressed using systems medicine approaches that include genome-scale metabolic models (GEMs). Previous studies have used a single model extraction method (MEM) and/or a single transcriptomic dataset to reconstruct context-specific models, which proved to be insufficient for the broader biological contexts. We have applied four MEMs in combination with five COVID-19 datasets. Models produced by GIMME were separated by infection, while tINIT preserved the biological variability in the data and enabled the best prediction of the enrichment of metabolic subsystems. Vitamin D3 metabolism was predicted to be down-regulated in one dataset by GIMME, and in all by tINIT. Models generated by tINIT and GIMME predicted downregulation of retinol metabolism in different datasets, while downregulated cholesterol metabolism was predicted only by tINIT-generated models. Predictions are in line with the observations in COVID-19 patients. Our data indicated that GIMME and tINIT models provided the most biologically relevant results and should have a larger emphasis in further analyses. Particularly tINIT models identified the metabolic pathways that are a part of the host response and are potential antiviral targets. The code and the results of the analyses are available to download from https://github.com/CompBioLj/COVID_GEMs_and_MEMs.

Guided extraction of genome-scale metabolic models for the integration and analysis of omics data.

Omics data can be integrated into a reference model using various model extraction methods (MEMs) to yield context-specific genome-scale metabolic models (GEMs). How to chose the appropriate MEM, thresholding rule and threshold remains a challenge. We integrated mouse transcriptomic data from a Cyp51 knockout mice diet experiment (GSE58271) using five MEMs (GIMME, iMAT, FASTCORE, INIT an tINIT) in a combination with a recently published mouse GEM iMM1865. Except for INIT and tINIT, the size of extracted models varied with the MEM used (t-test: p-value < 0.001). The Jaccard index of iMAT models ranged from 0.27 to 1.0. Out of the three factors under study in the experiment (diet, gender and genotype), gender explained most of the variability ( > 90%) in PC1 for FASTCORE. In iMAT, each of the three factors explained less than 40% of the variability within PC1, PC2 and PC3. Among all the MEMs, FASTCORE captured the most of the true variability in the data by clustering samples by gender. Our results show that for the efficient use of MEMs in the context of omics data integration and analysis, one should apply various MEMs, thresholding rules, and thresholding values to select the MEM and its configuration that best captures the true variability in the data. This selection can be guided by the methodology as proposed and used in this paper. Moreover, we describe certain approaches that can be used to analyse the results obtained with the selected MEM and to put these results in a biological context.

Common Transcriptional Program of Liver Fibrosis in Mouse Genetic Models and Humans.

Multifactorial metabolic diseases, such as non-alcoholic fatty liver disease, are a major burden to modern societies, and frequently present with no clearly defined molecular biomarkers. Herein we used system medicine approaches to decipher signatures of liver fibrosis in mouse models with malfunction in genes from unrelated biological pathways: cholesterol synthesis-Cyp51, notch signaling-Rbpj, nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) signaling-Ikbkg, and unknown lysosomal pathway-Glmp. Enrichment analyses of Kyoto Encyclopedia of Genes and Genomes (KEGG), Reactome and TRANScription FACtor (TRANSFAC) databases complemented with genome-scale metabolic modeling revealed fibrotic signatures highly similar to liver pathologies in humans. The diverse genetic models of liver fibrosis exposed a common transcriptional program with activated estrogen receptor alpha (ERα) signaling, and a network of interactions between regulators of lipid metabolism and transcription factors from cancer pathways and the immune system. The novel hallmarks of fibrosis are downregulated lipid pathways, including fatty acid, bile acid, and steroid hormone metabolism. Moreover, distinct metabolic subtypes of liver fibrosis were proposed, supported by unique enrichment of transcription factors based on the type of insult, disease stage, or potentially, also sex. The discovered novel features of multifactorial liver fibrotic pathologies could aid also in improved stratification of other fibrosis related pathologies.

Field-programmable biological circuits and configurable (bio)logic blocks for distributed biological computing.

Synthetic biology applications often require engineered computing structures, which can be programmed to process the information in a given way. However, programming of these structures usually requires significant amount of trial-and-error genetic engineering. This process is to some degree analogous to the design of application-specific integrated circuits (ASIC) in the domain of digital electronic circuits, which often require complex and time-consuming workflows to obtain a desired response. We describe a design of programmable biological circuits that can be configured without additional genetic engineering. Their configuration can be changed in vivo, i.e. during the execution of their biological program, simply with an introduction of programming inputs. These, e.g., increase the degradation rates of selected proteins that store the current configuration of the circuit. Programming can be thus performed in the field as in the case of field-programmable gate array (FPGA) circuits, which present an attractive alternative of ASICs in digital electronics. We describe a basic programmable unit, which we denote configurable (bio)logical block (CBLB) inspired by the architecture of configurable logic blocks (CLBs), basic functional units within the FPGA circuits. The design of a CBLB is based on distributed cellular computing modules, which makes its biological implementation easier to achieve. We establish a computational model of a CBLB and analyse its response with a given set of biologically feasible parameter values. Furthermore, we show that the proposed CBLB design exhibits correct behaviour for a vast range of kinetic parameter values, different population ratios, and as well preserves this response in stochastic simulations.

Correction to: Computational analysis of viable parameter regions in models of synthetic biological systems.

[This corrects the article DOI: 10.1186/s13036-019-0205-0.].

Computational analysis of viable parameter regions in models of synthetic biological systems.

BACKGROUND: Gene regulatory networks with different topological and/or dynamical properties might exhibit similar behavior. System that is less perceptive for the perturbations of its internal and external factors should be preferred. Methods for sensitivity and robustness assessment have already been developed and can be roughly divided into local and global approaches. Local methods focus only on the local area around nominal parameter values. This can be problematic when parameters exhibits the desired behavior over a large range of parameter perturbations or when parameter values are unknown. Global methods, on the other hand, investigate the whole space of parameter values and mostly rely on different sampling techniques. This can be computationally inefficient. To address these shortcomings 'glocal' approaches were developed that apply global and local approaches in an effective and rigorous manner. RESULTS: Herein, we present a computational approach for 'glocal' analysis of viable parameter regions in biological models. The methodology is based on the exploration of high-dimensional viable parameter spaces with global and local sampling, clustering and dimensionality reduction techniques. The proposed methodology allows us to efficiently investigate the viable parameter space regions, evaluate the regions which exhibit the largest robustness, and to gather new insights regarding the size and connectivity of the viable parameter regions. We evaluate the proposed methodology on three different synthetic gene regulatory network models, i.e. the repressilator model, the model of the AC-DC circuit and the model of the edge-triggered master-slave D flip-flop. CONCLUSIONS: The proposed methodology provides a rigorous assessment of the shape and size of viable parameter regions based on (1) the mathematical description of the biological system of interest, (2) constraints that define feasible parameter regions and (3) cost function that defines the desired or observed behavior of the system. These insights can be used to assess the robustness of biological systems, even in the case when parameter values are unknown and more importantly, even when there are multiple poorly connected viable parameter regions in the solution space. Moreover, the methodology can be efficiently applied to the analysis of biological systems that exhibit multiple modes of the targeted behavior.

Initial state perturbations as a validation method for data-driven fuzzy models of cellular networks.

BACKGROUND: Data-driven methods that automatically learn relations between attributes from given data are a popular tool for building mathematical models in computational biology. Since measurements are prone to errors, approaches dealing with uncertain data are especially suitable for this task. Fuzzy models are one such approach, but they contain a large amount of parameters and are thus susceptible to over-fitting. Validation methods that help detect over-fitting are therefore needed to eliminate inaccurate models. RESULTS: We propose a method to enlarge the validation datasets on which a fuzzy dynamic model of a cellular network can be tested. We apply our method to two data-driven dynamic models of the MAPK signalling pathway and two models of the mammalian circadian clock. We show that random initial state perturbations can drastically increase the mean error of predictions of an inaccurate computational model, while keeping errors of predictions of accurate models small. CONCLUSIONS: With the improvement of validation methods, fuzzy models are becoming more accurate and are thus likely to gain new applications. This field of research is promising not only because fuzzy models can cope with uncertainty, but also because their run time is short compared to conventional modelling methods that are nowadays used in systems biology.

Computational Modelling of Liver Metabolism and its Applications in Research and the Clinics.

Computational models of liver metabolism are gaining an increasing importance within the research community. Moreover, their first clinical applications have been reported in recent years in the context of personalised and systems medicine. Herein, we survey selected experimental models together with the computational modelling approaches that are used to describe the metabolic processes of the liver in silico. We also review the recent developments in the large-scale hepatic computational models where we focus on object-oriented models as a part of our research. The object-oriented modelling approach is beneficial in efforts to describe the interactions between the tissues, such as how metabolism of the liver interacts with metabolism of other tissues via blood. Importantly, this modelling approach can account as well for transcriptional and post-translational regulation of metabolic reactions which is a difficult task to achieve. The current and potential clinical applications of large-scale hepatic models are also discussed. We conclude with the future perspectives within the systems and translational medicine research community.

LiverSex Computational Model: Sexual Aspects in Hepatic Metabolism and Abnormalities.

The liver is to date the best example of a sexually dimorphic non-reproductive organ. Over 1,000 genes are differentially expressed between sexes indicating that female and male livers are two metabolically distinct organs. The spectrum of liver diseases is broad and is usually prevalent in one or the other sex, with different contributing genetic and environmental factors. It is thus difficult to predict individual's disease outcomes and treatment options. Systems approaches including mathematical modeling can aid importantly in understanding the multifactorial liver disease etiology leading toward tailored diagnostics, prognostics and therapy. The currently established computational models of hepatic metabolism that have proven to be essential for understanding of non-alcoholic fatty liver disease (NAFLD) and hepatocellular carcinoma (HCC) are limited to the description of gender-independent response or reflect solely the response of the males. Herein we present LiverSex, the first sex-based multi-tissue and multi-level liver metabolic computational model. The model was constructed based on in silico liver model SteatoNet and the object-oriented modeling. The crucial factor in adaptation of liver metabolism to the sex is the inclusion of estrogen and androgen receptor responses to respective hormones and the link to sex-differences in growth hormone release. The model was extensively validated on literature data and experimental data obtained from wild type C57BL/6 mice fed with regular chow and western diet. These experimental results show extensive sex-dependent changes and could not be reproduced in silico with the uniform model SteatoNet. LiverSex represents the first large-scale liver metabolic model, which allows a detailed insight into the sex-dependent complex liver pathologies, and how the genetic and environmental factors interact with the sex in disease appearance and progression. We used the model to identify the most important sex-dependent metabolic pathways, which are involved in accumulation of triglycerides representing initial steps of NAFLD. We identified PGC1A, PPARα, FXR, and LXR as regulatory factors that could become important in sex-dependent personalized treatment of NAFLD.

Grohar: Automated Visualization of Genome-Scale Metabolic Models and Their Pathways.

Genome-scale metabolic models (GEMs) have become a powerful tool for the investigation of the entire metabolism of the organism in silico. These models are, however, often extremely hard to reconstruct and also difficult to apply to the selected problem. Visualization of the GEM allows us to easier comprehend the model, to perform its graphical analysis, to find and correct the faulty relations, to identify the parts of the system with a designated function, etc. Even though several approaches for the automatic visualization of GEMs have been proposed, metabolic maps are still manually drawn or at least require large amount of manual curation. We present Grohar, a computational tool for automatic identification and visualization of GEM (sub)networks and their metabolic fluxes. These (sub)networks can be specified directly by listing the metabolites of interest or indirectly by providing reference metabolic pathways from different sources, such as KEGG, SBML, or Matlab file. These pathways are identified within the GEM using three different pathway alignment algorithms. Grohar also supports the visualization of the model adjustments (e.g., activation or inhibition of metabolic reactions) after perturbations are induced.

Computational modelling of genome-scale metabolic networks and its application to CHO cell cultures.

Genome-scale metabolic models (GEMs) have become increasingly important in recent years. Currently, GEMs are the most accurate in silico representation of the genotype-phenotype link. They allow us to study complex networks from the systems perspective. Their application may drastically reduce the amount of experimental and clinical work, improve diagnostic tools and increase our understanding of complex biological phenomena. GEMs have also demonstrated high potential for the optimisation of bio-based production of recombinant proteins. Herein, we review the basic concepts, methods, resources and software tools used for the reconstruction and application of GEMs. We overview the evolution of the modelling efforts devoted to the metabolism of Chinese Hamster Ovary (CHO) cells. We present a case study on CHO cell metabolism under different amino acid depletions. This leads us to the identification of the most influential as well as essential amino acids in selected CHO cell lines.

Large-scale computational models of liver metabolism: How far from the clinics?

Understanding the dynamics of human liver metabolism is fundamental for effective diagnosis and treatment of liver diseases. This knowledge can be obtained with systems biology/medicine approaches that account for the complexity of hepatic responses and their systemic consequences in other organs. Computational modeling can reveal hidden principles of the system by classification of individual components, analyzing their interactions and simulating the effects that are difficult to investigate experimentally. Herein, we review the state-of-the-art computational models that describe liver dynamics from metabolic, gene regulatory, and signal transduction perspectives. We focus especially on large-scale liver models described either by genome scale metabolic networks or an object-oriented approach. We also discuss the benefits and limitations of each modeling approach and their value for clinical applications in diagnosis, therapy, and prevention of liver diseases as well as precision medicine in hepatology. (Hepatology 2017;66:1323-1334).

Classical Mechanics Approach Applied to Analysis of Genetic Oscillators.

Biological oscillators present a fundamental part of several regulatory mechanisms that control the response of various biological systems. Several analytical approaches for their analysis have been reported recently. They are, however, limited to only specific oscillator topologies and/or to giving only qualitative answers, i.e., is the dynamics of an oscillator given the parameter space oscillatory or not. Here, we present a general analytical approach that can be applied to the analysis of biological oscillators. It relies on the projection of biological systems to classical mechanics systems. The approach is able to provide us with relatively accurate results in the meaning of type of behavior system reflects (i.e., oscillatory or not) and periods of potential oscillations without the necessity to conduct expensive numerical simulations. We demonstrate and verify the proposed approach on three different implementations of amplified negative feedback oscillator.

Computational Framework for Modeling Multiple Noncooperative Transcription Factor Binding and Its Application to the Analysis of Nuclear Factor Kappa B Oscillatory Response.

Recent studies have shown that regulation of many important genes is achieved with multiple transcription factor (TF) binding sites with low or no cooperativity. Additionally, noncooperative binding sites are gaining more and more importance in the field of synthetic biology. Here, we introduce a computational framework that can be applied to dynamical modeling and analysis of gene regulatory networks with multiple noncooperative TF binding sites. We propose two computational methods to be used within the framework, that is, average promoter state approximation and expression profiles based modeling. We demonstrate the application of the proposed framework on the analysis of nuclear factor kappa B (NF-κB) oscillatory response. We show that different promoter expression hypotheses in a combination with the number of TF binding sites drastically affect the dynamics of the observed system and should not be ignored in the process of quantitative dynamical modeling, as is usually the case in existent state-of-the-art computational analyses.

Fuzzy Logic as a Computational Tool for Quantitative Modelling of Biological Systems with Uncertain Kinetic Data.

Quantitative modelling of biological systems has become an indispensable computational approach in the design of novel and analysis of existing biological systems. However, kinetic data that describe the system's dynamics need to be known in order to obtain relevant results with the conventional modelling techniques. These data are often hard or even impossible to obtain. Here, we present a quantitative fuzzy logic modelling approach that is able to cope with unknown kinetic data and thus produce relevant results even though kinetic data are incomplete or only vaguely defined. Moreover, the approach can be used in the combination with the existing state-of-the-art quantitative modelling techniques only in certain parts of the system, i.e., where kinetic data are missing. The case study of the approach proposed here is performed on the model of three-gene repressilator.

Stochastic simulation algorithm for gene regulatory networks with multiple binding sites.

Promoters with multiple binding sites present a regulatory mechanism of several natural biological systems. It has been shown that such systems reflect a higher stability in comparison to the systems with small numbers of binding sites. Regulatory mechanisms with multiple binding sites are therefore used more frequently in artificially designed biological systems in recent years. While the number of possible promoter states increases exponentially with the number of binding sites, it is extremely hard to model such systems accurately. Here we present an adaptation of stochastic simulation algorithm for accurate modeling of gene regulatory networks with multiple binding sites. Small computational complexity of adapted algorithm allows us to model any feasible number of binding sites per promoter. The approach introduced in this work is demonstrated on the model of switching mechanism in Epstein-Barr virus, where 20 binding sites are observed on one of the promoters. We show that the presented approach is easy to adapt to any biological systems based on the regulatory mechanisms with multiple binding sites in order to obtain and analyze their behavior.