Dear colleagues, We thank you for your comments on our manuscript [1] and for raising the important question of automation [2]. The publications by Sepp et al. (2019) [3] and Liu et al. (2024) [4] focus on one type of PBPK models, denominated ‘for biologics’. In these models, each compartment (tissue/organ) is divided into vascular, endothelial endosomal and interstitial (sub-)compartments, and lymphatic flows are also represented [3, 4]. In these examples, the pharmacokinetics of proteins are represented. Notably, the protein distribution within tissues includes processes of passive transport across two types of pores via diffusion or fluid convection, processes of pinocytosis, binding, recycling and degradation [4]. This type of PBPK models is complex and their authors successfully used specific tools for automated code generation. In particular, Liu et al. (2024) used ‘mathematical sets’ available in the Ubiquity package, which facilitates assembling model components in R-language [4]. Sepp et al. (2019) used the MATLAB script PBPKassembler.m for automated code generation, where files in Simbiology (MATLAB) and Excel formats were combined [3]. Briefly, automation can help build models, which is highly valuable notably for complex models. Having a robust, thoroughly validated platform for automated code generation will indeed contribute to accessibility and reproducibility. Also, automation of model building seems to be a natural way to address the case of large models, in which manual scripting will be prone to introducing errors. Smaller models or those with fewer equations do not necessarily require automation. From our experience, writing the script manually provides a deep understanding of the underlying model mechanics. Manually written scripts also allow to give a clear view of how equations and parameters are applied, which we consider especially valuable for teaching and knowledge transfer. To conclude, as observed in other niches of systems biology, automation brings significant benefits to the modelling field, e.g. there is a plethora of tools that reconstruct and benchmark genome-scale metabolic models. As the correspondence authors argue, there is more scope for automation in the current practices of PBPK modelling. To facilitate this, we have opted to begin the journey of standardization via open collaboration through ELIXIR. We are looking forward to continuing to engage the research community towards automatic construction, validation and deposition of PBPK models. Thank you. Best regards, Elena Domínguez-Romero, Stanislav Mazurenko, Martin Scheringer, Vítor Martins dos Santos, Chris Evelo, Mihail Anton, John M. Hancock, Anže Županič, and Maria Suarez-Diez. None declared.

Systems biology aims to understand living organisms through mathematically modeling their behaviors at different organizational levels, ranging from molecules to populations. Modeling involves several steps, from determining the model purpose to developing the mathematical model, implementing it computationally, simulating the model's behavior, evaluating, and refining the model. Importantly, model simulation results must be reproducible, ensuring that other researchers can obtain the same results after writing the code de novo and/or using different software tools. Guidelines to increase model reproducibility have been published. However, reproducibility remains a major challenge in this field. In this paper, we tackle this challenge for physiologically-based pharmacokinetic (PBPK) models, which represent the pharmacokinetics of chemicals following exposure in humans or animals. We summarize recommendations for PBPK model reporting that should apply during model development and implementation, in order to ensure model reproducibility and comprehensibility. We make a proposal aiming to harmonize abbreviations used in PBPK models. To illustrate these recommendations, we present an original and reproducible PBPK model code in MATLAB, alongside an example of MATLAB code converted to Systems Biology Markup Language format using MOCCASIN. As directions for future improvement, more tools to convert computational PBPK models from different software platforms into standard formats would increase the interoperability of these models. The application of other systems biology standards to PBPK models is encouraged. This work is the result of an interdisciplinary collaboration involving the ELIXIR systems biology community. More interdisciplinary collaborations like this would facilitate further harmonization and application of good modeling practices in different systems biology fields.

Harnessing the power of microbial consortia is integral to a diverse range of sectors, from healthcare to biotechnology to environmental remediation. To fully realize this potential, it is critical to understand the mechanisms behind the interactions that structure microbial consortia and determine their functions. Constraint-based reconstruction and analysis (COBRA) approaches, employing genome-scale metabolic models (GEMs), have emerged as the state-of-the-art tool to simulate the behavior of microbial communities from their constituent genomes. In the last decade, many tools have been developed that use COBRA approaches to simulate multi-species consortia, under either steady-state, dynamic, or spatiotemporally varying scenarios. Yet, these tools have not been systematically evaluated regarding their software quality, most suitable application, and predictive power. Hence, it is uncertain which tools users should apply to their system and what are the most urgent directions that developers should take in the future to improve existing capacities. This study conducted a systematic evaluation of COBRA-based tools for microbial communities using datasets from two-member communities as test cases. First, we performed a qualitative assessment in which we evaluated 24 published tools based on a list of FAIR (Findability, Accessibility, Interoperability, and Reusability) features essential for software quality. Next, we quantitatively tested the predictions in a subset of 14 of these tools against experimental data from three different case studies: a) syngas fermentation by C. autoethanogenum and C. kluyveri for the static tools, b) glucose/xylose fermentation with engineered E. coli and S. cerevisiae for the dynamic tools, and c) a Petri dish of E. coli and S. enterica for tools incorporating spatiotemporal variation. Our results show varying performance levels of the best qualitatively assessed tools when examining the different categories of tools. The differences in the mathematical formulation of the approaches and their relation to the results were also discussed. Ultimately, we provide recommendations for refining future GEM microbial modeling tools.

Atlantic salmon (Salmo salar) is the most valuable farmed fish globally and there is much interest in optimizing its genetics and rearing conditions for growth and feed efficiency. Marine feed ingredients must be replaced to meet global demand, with challenges for fish health and sustainability. Metabolic models can address this by connecting genomes to metabolism, which converts nutrients in the feed to energy and biomass, but such models are currently not available for major aquaculture species such as salmon. We present SALARECON, a model focusing on energy, amino acid, and nucleotide metabolism that links the Atlantic salmon genome to metabolic fluxes and growth. It performs well in standardized tests and captures expected metabolic (in)capabilities. We show that it can explain observed hypoxic growth in terms of metabolic fluxes and apply it to aquaculture by simulating growth with commercial feed ingredients. Predicted limiting amino acids and feed efficiencies agree with data, and the model suggests that marine feed efficiency can be achieved by supplementing a few amino acids to plant- and insect-based feeds. SALARECON is a high-quality model that makes it possible to simulate Atlantic salmon metabolism and growth. It can be used to explain Atlantic salmon physiology and address key challenges in aquaculture such as development of sustainable feeds.

BACKGROUND: Several computational methods have been developed that integrate transcriptomics data with genome-scale metabolic reconstructions to increase accuracy of inferences of intracellular metabolic flux distributions. Even though existing methods use transcript abundances as a proxy for enzyme activity, each method uses a different hypothesis and assumptions. Most methods implicitly assume a proportionality between transcript levels and flux through the corresponding function, although these proportionality constant(s) are often not explicitly mentioned nor discussed in any of the published methods. E-Flux is one such method and, in this algorithm, flux bounds are related to expression data, so that reactions associated with highly expressed genes are allowed to carry higher flux values. RESULTS: Here, we extended E-Flux and systematically evaluated the impact of an assumed proportionality constant on model predictions. We used data from published experiments with Escherichia coli and Saccharomyces cerevisiae and we compared the predictions of the algorithm to measured extracellular and intracellular fluxes. CONCLUSION: We showed that detailed modelling using a proportionality constant can greatly impact the outcome of the analysis. This increases accuracy and allows for extraction of better physiological information.

Summary: To unlock the full potential of genome data and to enhance data interoperability and reusability of genome annotations we have developed SAPP, a Semantic Annotation Platform with Provenance. SAPP is designed as an infrastructure supporting FAIR de novo computational genomics but can also be used to process and analyze existing genome annotations. SAPP automatically predicts, tracks and stores structural and functional annotations and associated dataset- and element-wise provenance in a Linked Data format, thereby enabling information mining and retrieval with Semantic Web technologies. This greatly reduces the administrative burden of handling multiple analysis tools and versions thereof and facilitates multi-level large scale comparative analysis. Availability and implementation: SAPP is written in JAVA and freely available at https://gitlab.com/sapp and runs on Unix-like operating systems. The documentation, examples and a tutorial are available at https://sapp.gitlab.io. Contact: [email protected] or [email protected].

BACKGROUND: Systems biology takes a holistic approach by handling biomolecules and their interactions as big systems. Network based approach has emerged as a natural way to model these systems with the idea of representing biomolecules as nodes and their interactions as edges. Very often the input data come from various sorts of omics analyses. Those resulting networks sometimes describe a wide range of aspects, for example different experiment conditions, species, tissue types, stimulating factors, mutants, or simply distinct interaction features of the same network produced by different algorithms. For these scenarios, synchronous visualization of more than one distinct network is an excellent mean to explore all the relevant networks efficiently. In addition, complementary analysis methods are needed and they should work in a workflow manner in order to gain maximal biological insights. RESULTS: In order to address the aforementioned needs, we have developed a Synchronous Network Data Integration (SyNDI) framework. This framework contains SyncVis, a Cytoscape application for user-friendly synchronous and simultaneous visualization of multiple biological networks, and it is seamlessly integrated with other bioinformatics tools via the Galaxy platform. We demonstrated the functionality and usability of the framework with three biological examples - we analyzed the distinct connectivity of plasma metabolites in networks associated with high or low latent cardiovascular disease risk; deeper insights were obtained from a few similar inflammatory response pathways in Staphylococcus aureus infection common to human and mouse; and regulatory motifs which have not been reported associated with transcriptional adaptations of Mycobacterium tuberculosis were identified. CONCLUSIONS: Our SyNDI framework couples synchronous network visualization seamlessly with additional bioinformatics tools. The user can easily tailor the framework for his/her needs by adding new tools and datasets to the Galaxy platform.