Perspectives on computational modeling of biological systems and the significance of the SysMod community.

Motivation: In recent years, applying computational modeling to systems biology has caused a substantial surge in both discovery and practical applications and a significant shift in our understanding of the complexity inherent in biological systems.

Results: In this perspective article, we briefly overview computational modeling in biology, highlighting recent advancements such as multi-scale modeling due to the omics revolution, single-cell technology, and integration of artificial intelligence and machine learning approaches. We also discuss the primary challenges faced: integration, standardization, model complexity, scalability, and interdisciplinary collaboration.

LACE 2.0: an interactive R tool for the inference and visualization of longitudinal cancer evolution.

Background: Longitudinal single-cell sequencing experiments of patient-derived models are increasingly employed to investigate cancer evolution. In this context, robust computational methods are needed to properly exploit the mutational profiles of single cells generated via variant calling, in order to reconstruct the evolutionary history of a tumor and characterize the impact of therapeutic strategies, such as the administration of drugs. To this end, we have recently developed the LACE framework for the Longitudinal Analysis of Cancer Evolution.

Combining denoising of RNA-seq data and flux balance analysis for cluster analysis of single cells.

Background: Sophisticated methods to properly pre-process and analyze the increasing collection of single-cell RNA sequencing (scRNA-seq) data are increasingly being developed. On the contrary, the best practices to integrate these data into metabolic networks, aiming at describing metabolic phenotypes within a heterogeneous cell population, have been poorly investigated. In this regard, a critical factor is the presence of false zero values in reactions essential for a fundamental metabolic function, such as biomass or energy production.

Machine Learning for Efficient Prediction of Protein Redox Potential: The Flavoproteins Case.

Determining the redox potentials of protein cofactors and how they are influenced by their molecular neighborhoods is essential for basic research and many biotechnological applications, from biosensors and biocatalysis to bioremediation and bioelectronics. The laborious determination of redox potential with current experimental technologies pushes forward the need for computational approaches that can reliably predict it. Although current computational approaches based on quantum and molecular mechanics are accurate, their large computational costs hinder their usage.

GPRuler: Metabolic gene-protein-reaction rules automatic reconstruction.

Metabolic network models are increasingly being used in health care and industry. As a consequence, many tools have been released to automate their reconstruction process de novo. In order to enable gene deletion simulations and integration of gene expression data, these networks must include gene-protein-reaction (GPR) rules, which describe with a Boolean logic relationships between the gene products (e.

PMCE: efficient inference of expressive models of cancer evolution with high prognostic power.

Motivation: Driver (epi)genomic alterations underlie the positive selection of cancer subpopulations, which promotes drug resistance and relapse. Even though substantial heterogeneity is witnessed in most cancer types, mutation accumulation patterns can be regularly found and can be exploited to reconstruct predictive models of cancer evolution. Yet, available methods can not infer logical formulas connecting events to represent alternative evolutionary routes or convergent evolution.

On the Use of Topological Features of Metabolic Networks for the Classification of Cancer Samples.

Background: The increasing availability of omics data collected from patients affected by severe pathologies, such as cancer, is fostering the development of data science methods for their analysis.

Introduction: The combination of data integration and machine learning approaches can provide new powerful instruments to tackle the complexity of cancer development and deliver effective diagnostic and prognostic strategies.

Methods: We explore the possibility of exploiting the topological properties of sample-specific metabolic networks as features in a supervised classification task.

Accelerated global sensitivity analysis of genome-wide constraint-based metabolic models.

Background: Genome-wide reconstructions of metabolism opened the way to thorough investigations of cell metabolism for health care and industrial purposes. However, the predictions offered by Flux Balance Analysis (FBA) can be strongly affected by the choice of flux boundaries, with particular regard to the flux of reactions that sink nutrients into the system. To mitigate possible errors introduced by a poor selection of such boundaries, a rational approach suggests to focus the modeling efforts on the pivotal ones.

MaREA4Galaxy: Metabolic reaction enrichment analysis and visualization of RNA-seq data within Galaxy.

We present MaREA4Galaxy, a user-friendly tool that allows a user to characterize and to graphically compare groups of samples with different transcriptional regulation of metabolism, as estimated from cross-sectional RNA-seq data. The tool is available as plug-in for the widely-used Galaxy platform for comparative genomics and bioinformatics analyses. MaREA4Galaxy combines three modules.

Systems metabolomics: from metabolomic snapshots to design principles.

Metabolomics is a rapidly expanding technology that finds increasing application in a variety of fields, form metabolic disorders to cancer, from nutrition and wellness to design and optimization of cell factories. The integration of metabolic snapshots with metabolic fluxes, physiological readouts, metabolic models, and knowledge-informed Artificial Intelligence tools, is required to obtain a system-level understanding of metabolism. The emerging power of multi-omic approaches and the development of integrated experimental and computational tools, able to dissect metabolic features at cellular and subcellular resolution, provide unprecedented opportunities for understanding design principles of metabolic (dis)regulation and for the development of precision therapies in multifactorial diseases, such as cancer and neurodegenerative diseases.

Single-cell Digital Twins for Cancer Preclinical Investigation.

Laboratory models derived from clinical samples represent a solid platform in preclinical research for drug testing and investigation of disease mechanisms. The integration of these laboratory models with their digital counterparts (i.e.

Integration of single-cell RNA-seq data into population models to characterize cancer metabolism.

Metabolic reprogramming is a general feature of cancer cells. Regrettably, the comprehensive quantification of metabolites in biological specimens does not promptly translate into knowledge on the utilization of metabolic pathways. By estimating fluxes across metabolic pathways, computational models hold the promise to bridge this gap between data and biological functionality.

Integration of transcriptomic data and metabolic networks in cancer samples reveals highly significant prognostic power.

Effective stratification of cancer patients on the basis of their molecular make-up is a key open challenge. Given the altered and heterogenous nature of cancer metabolism, we here propose to use the overall expression of central carbon metabolism as biomarker to characterize groups of patients with important characteristics, such as response to ad-hoc therapeutic strategies and survival expectancy. To this end, we here introduce the data integration framework named Metabolic Reaction Enrichment Analysis (MaREA), which strives to characterize the metabolic deregulations that distinguish cancer phenotypes, by projecting RNA-seq data onto metabolic networks, without requiring metabolic measurements.

Emerging ensembles of kinetic parameters to characterize observed metabolic phenotypes.

Background: Determining the value of kinetic constants for a metabolic system in the exact physiological conditions is an extremely hard task. However, this kind of information is of pivotal relevance to effectively simulate a biological phenomenon as complex as metabolism.

Results: To overcome this issue, we propose to investigate emerging properties of ensembles of sets of kinetic constants leading to the biological readout observed in different experimental conditions.

A metabolic core model elucidates how enhanced utilization of glucose and glutamine, with enhanced glutamine-dependent lactate production, promotes cancer cell growth: The WarburQ effect.

Cancer cells share several metabolic traits, including aerobic production of lactate from glucose (Warburg effect), extensive glutamine utilization and impaired mitochondrial electron flow. It is still unclear how these metabolic rearrangements, which may involve different molecular events in different cells, contribute to a selective advantage for cancer cell proliferation. To ascertain which metabolic pathways are used to convert glucose and glutamine to balanced energy and biomass production, we performed systematic constraint-based simulations of a model of human central metabolism.

popFBA: tackling intratumour heterogeneity with Flux Balance Analysis.

Motivation: Intratumour heterogeneity poses many challenges to the treatment of cancer. Unfortunately, the transcriptional and metabolic information retrieved by currently available computational and experimental techniques portrays the average behaviour of intermixed and heterogeneous cell subpopulations within a given tumour. Emerging single-cell genomic analyses are nonetheless unable to characterize the interactions among cancer subpopulations.

Divergent in vitro/in vivo responses to drug treatments of highly aggressive NIH-Ras cancer cells: a PET imaging and metabolomics-mass-spectrometry study.

Oncogenic K-ras is capable to control tumor growth and progression by rewiring cancer metabolism. In vitro NIH-Ras cells convert glucose to lactate and use glutamine to sustain anabolic processes, but their in vivo environmental adaptation and multiple metabolic pathways activation ability is poorly understood. Here, we show that NIH-Ras cancer cells and tumors are able to coordinate nutrient utilization to support aggressive cell proliferation and survival.

Zooming-in on cancer metabolic rewiring with tissue specific constraint-based models.

The metabolic rearrangements occurring in cancer cells can be effectively investigated with a Systems Biology approach supported by metabolic network modeling. We here present tissue-specific constraint-based core models for three different types of tumors (liver, breast and lung) that serve this purpose. The core models were extracted and manually curated from the corresponding genome-scale metabolic models in the Human Metabolic Atlas database with a focus on the pathways that are known to play a key role in cancer growth and proliferation.

CABeRNET: a Cytoscape app for augmented Boolean models of gene regulatory NETworks.

Background: Dynamical models of gene regulatory networks (GRNs) are highly effective in describing complex biological phenomena and processes, such as cell differentiation and cancer development. Yet, the topological and functional characterization of real GRNs is often still partial and an exhaustive picture of their functioning is missing.

Results: We here introduce CABERNET, a Cytoscape app for the generation, simulation and analysis of Boolean models of GRNs, specifically focused on their augmentation when a only partial topological and functional characterization of the network is available.

Lipid-based nanoparticles containing cationic derivatives of PTA (1,3,5-triaza-7-phosphaadamantane) as innovative vehicle for Pt complexes: Production, characterization and in vitro studies.

The aliphatic phosphine PTA (1,3,5-triaza-7-phosphaadamantane) is a promising ligand for metal complexes designed and developed as innovative inorganic drugs. In the present paper, an N-alkylated derivative of PTA, (PTAC16H33)X (X=I, C1, or X=PF6, C2) and its platinum coordination complex cis-[PtCl2(PTAC16H33)2](PF6)2, C3, were considered as components of cationic lipid nanoparticles (CLNs). Particularly, CLN1, CLN2 and CLN3 were obtained by adding derivatives C1, C2 or C3 during nanoparticles preparation, while CLN2-Pt were obtained by treating preformed CLN2 with Pt(II).

Growth and division in a dynamic protocell model.

In this paper a new model of growing and dividing protocells is described, whose main features are (i) a lipid container that grows according to the composition of the molecular milieu (ii) a set of "genetic memory molecules" (GMMs) that undergo catalytic reactions in the internal aqueous phase and (iii) a set of stochastic kinetic equations for the GMMs. The mass exchange between the external environment and the internal phase is described by simulating a semipermeable membrane and a flow driven by the differences in chemical potentials, thereby avoiding to resort to sometimes misleading simplifications, e.g.

Computational strategies for a system-level understanding of metabolism.

Cell metabolism is the biochemical machinery that provides energy and building blocks to sustain life. Understanding its fine regulation is of pivotal relevance in several fields, from metabolic engineering applications to the treatment of metabolic disorders and cancer. Sophisticated computational approaches are needed to unravel the complexity of metabolism.