Mathematical and Computational Modeling in Complex Biological Systems

Advancing the quantitative understanding of adverse outcome pathways: current status, methodologies, and future directions

An adverse outcome pathway (AOP) framework maps the sequence of events leading to adverse outcomes from chemical exposures, providing a mechanistic understanding often absent in traditional methods. The quantitative AOP (qAOP) advances AOP by integrating quantitative data and mathematical modeling, thereby providing a more precise comprehension of relationships between molecular initiating events, key events, and adverse outcomes. This review critically examines three primary methodologies: systems toxicology, regression modeling, and Bayesian network modeling, highlighting their strengths, limitations, and specific data requirements within toxicology. Through an analysis of current methodologies and challenges, this review emphasizes the integration of experimental and computational approaches to elucidate key event relationships and proposes strategies for overcoming limitations through standardized protocols and advanced computational tools. By outlining future research directions and the potential of qAOPs to transform chemical risk assessment, this review aims to contribute to the advancement of regulatory science and the protection of public health and the environment.

Exploring the molecular landscape of osteosarcoma through PTTG family genes using a detailed multi-level methodology

Background

Osteosarcoma (OS) poses a significant clinical challenge, necessitating a comprehensive exploration of its molecular underpinnings.

Methods

This study explored the roles of PTTG family genes (PTTG1, PTTG2, and PTTG3P) in OS, employing a multifaceted approach encompassing molecular experiments, including OS cell lines culturing, RT-qPCR, bisulfite and Whole Exome Sequencing (WES) and in silico experiments, including The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets-based validation, overall survival, gene enrichment, functional assays, and molecular docking analyses.

Results

Our findings reveal a consistent up-regulation of PTTG genes in OS cell lines, supported by RT-qPCR experiments and corroborated across various publically available expression datasets databases. Importantly, ROC curve analyses highlight their potential as diagnostic markers. Moving beyond expression profiles, we unveil the epigenetic landscape by demonstrating significant hypomethylation of CpG islands associated with PTTG genes in OS. The negative correlation between methylation status and mRNA expression emphasizes the regulatory role of promoter methylation in PTTG gene expression. Contrary to expectations, genetic mutations in PTTG genes are rare in OS, with only benign mutations observed. Moreover, functional assays also confirmed the oncogenic roles of the PTTG gene in the development of OS. Lastly, we also revealed that Calcitriol is the most appropriate drug that can be utilized to treat OS in the context of PTTG genes.

Conclusion

The identification of PTTG genes as potential diagnostic markers and their association with epigenetic alterations opens new avenues for understanding OS pathogenesis and developing targeted therapies. As we navigate the complex landscape of OS, this study contributes essential insights that may pave the way for improved diagnostic and therapeutic strategies in its management.

Aflibercept Off-Target Effects in Diabetic Macular Edema: An In Silico Modeling Approach

Intravitreal aflibercept injection (IAI) is a treatment for diabetic macular edema (DME), but its mechanism of action (MoA) has not been completely elucidated. Here, we aimed to explore IAI's MoA and its multi-target nature in DME pathophysiology with an in silico (computer simulation) disease model. We used the Therapeutic Performance Mapping System (Anaxomics Biotech property) to generate mathematical models based on the available scientific knowledge at the time of the study, describing the relationship between the modulation of vascular endothelial growth factor receptors (VEGFRs) by IAI and DME pathophysiological processes. We also undertook an enrichment analysis to explore the processes modulated by IAI, visualized the effectors' predicted protein activity, and specifically evaluated the role of VEGFR1 pathway inhibition on DME treatment. The models simulated the potential pathophysiology of DME and the likely IAI's MoA by inhibiting VEGFR1 and VEGFR2 signaling. The action of IAI through both signaling pathways modulated the identified pathophysiological processes associated with DME, with the strongest effects in angiogenesis, blood-retinal barrier alteration and permeability, and inflammation. VEGFR1 inhibition was essential to modulate inflammatory protein effectors. Given the role of VEGFR1 signaling on the modulation of inflammatory-related pathways, IAI may offer therapeutic advantages for DME through sustained VEGFR1 pathway inhibition.

Incorporating temporal information during feature engineering bolsters emulation of spatio-temporal emergence

MOTIVATION

Emergent biological dynamics derive from the evolution of lower-level spatial and temporal processes. A long-standing challenge for scientists and engineers is identifying simple low-level rules that give rise to complex higher-level dynamics. High-resolution biological data acquisition enables this identification and has evolved at a rapid pace for both experimental and computational approaches. Simultaneously harnessing the resolution and managing the expense of emerging technologies-e.g. live cell imaging, scRNAseq, agent-based models-requires a deeper understanding of how spatial and temporal axes impact biological systems. Effective emulation is a promising solution to manage the expense of increasingly complex high-resolution computational models. In this research, we focus on the emulation of a tumor microenvironment agent-based model to examine the relationship between spatial and temporal environment features, and emergent tumor properties.

RESULTS

Despite significant feature engineering, we find limited predictive capacity of tumor properties from initial system representations. However, incorporating temporal information derived from intermediate simulation states dramatically improves the predictive performance of machine learning models. We train a deep-learning emulator on intermediate simulation states and observe promising enhancements over emulators trained solely on initial conditions. Our results underscore the importance of incorporating temporal information in the evaluation of spatio-temporal emergent behavior. Nevertheless, the emulators exhibit inconsistent performance, suggesting that the underlying model characterizes unique cell populations dynamics that are not easily replaced.

AVAILABILITY AND IMPLEMENTATION

All source codes for the agent-based model, emulation, and analyses are publicly available at the corresponding DOIs: 10.5281/zenodo.10622155, 10.5281/zenodo.10611675, 10.5281/zenodo.10621244, respectively.

Hair and Nail-On-Chip for Bioinspired Microfluidic Device Fabrication and Biomarker Detection

The advent of biosensors has tremendously increased our potential of identifying and solving important problems in various domains, ranging from food safety and environmental analysis, to healthcare and medicine. However, one of the most prominent drawbacks of these technologies, especially in the biomedical field, is to employ conventional samples, such as blood, urine, tissue extracts and other body fluids for analysis, which suffer from the drawbacks of invasiveness, discomfort, and high costs encountered in transportation and storage, thereby hindering these products to be applied for point-of-care testing that has garnered substantial attention in recent years. Therefore, through this review, we emphasize for the first time, the applications of switching over to noninvasive sampling techniques involving hair and nails that not only circumvent most of the aforementioned limitations, but also serve as interesting alternatives in understanding the human physiology involving minimal costs, equipment and human interference when combined with rapidly advancing technologies, such as microfluidics and organ-on-a-chip to achieve miniaturization on an unprecedented scale. The coalescence between these two fields has not only led to the fabrication of novel microdevices involving hair and nails, but also function as robust biosensors for the detection of biomarkers, chemicals, metabolites and nucleic acids through noninvasive sampling. Finally, we have also elucidated a plethora of futuristic innovations that could be incorporated in such devices, such as expanding their applications in nail and hair-based drug delivery, their potential in serving as next-generation wearable sensors and integrating these devices with machine-learning for enhanced automation and decentralization.

Modelling the Tumour Microenvironment, but What Exactly Do We Mean by "Model"?

The Oxford English Dictionary includes 17 definitions for the word "model" as a noun and another 11 as a verb. Therefore, context is necessary to understand the meaning of the word model. For instance, "model railways" refer to replicas of railways and trains at a smaller scale and a "model student" refers to an exemplary individual. In some cases, a specific context, like cancer research, may not be sufficient to provide one specific meaning for model. Even if the context is narrowed, specifically, to research related to the tumour microenvironment, "model" can be understood in a wide variety of ways, from an animal model to a mathematical expression. This paper presents a review of different "models" of the tumour microenvironment, as grouped by different definitions of the word into four categories: model organisms, in vitro models, mathematical models and computational models. Then, the frequencies of different meanings of the word "model" related to the tumour microenvironment are measured from numbers of entries in the MEDLINE database of the United States National Library of Medicine at the National Institutes of Health. The frequencies of the main components of the microenvironment and the organ-related cancers modelled are also assessed quantitatively with specific keywords. Whilst animal models, particularly xenografts and mouse models, are the most commonly used "models", the number of these entries has been slowly decreasing. Mathematical models, as well as prognostic and risk models, follow in frequency, and these have been growing in use.

Comparison of rule- and ordinary differential equation-based dynamic model of DARPP-32 signalling network

Dynamic modelling has considerably improved our understanding of complex molecular mechanisms. Ordinary differential equations (ODEs) are the most detailed and popular approach to modelling the dynamics of molecular systems. However, their application in signalling networks, characterised by multi-state molecular complexes, can be prohibitive. Contemporary modelling methods, such as rule- based (RB) modelling, have addressed these issues. The advantages of RB modelling over ODEs have been presented and discussed in numerous reviews. In this study, we conduct a direct comparison of the time courses of a molecular system founded on the same reaction network but encoded in the two frameworks. To make such a comparison, a set of reactions that underlie an ODE model was manually encoded in the Kappa language, one of the RB implementations. A comparison of the models was performed at the level of model specification and dynamics, acquired through model simulations. In line with previous reports, we confirm that the Kappa model recapitulates the general dynamics of its ODE counterpart with minor differences. These occur when molecules have multiple sites binding the same interactor. Furthermore, activation of these molecules in the RB model is slower than in the ODE one. As reported for other molecular systems, we find that, also for the DARPP-32 reaction network, the RB representation offers a more expressive and flexible syntax that facilitates access to fine details of the model, easing model reuse. In parallel with these analyses, we report a refactored model of the DARPP-32 interaction network that can serve as a canvas for the development of more complex dynamic models to study this important molecular system.

Gene Regulatory Identification Based on the Novel Hybrid Time-Delayed Method

Gene regulatory network (GRN) inference with biology data is a difficult and serious issue in the field of system biology. In order to detect the direct associations of GRN more accurately, a novel two-step GRN inference technique based on the time-delayed correlation coefficient (TDCC) and time-delayed complex-valued S-system model (TDCVSS) is proposed. First, a TDCC algorithm is utilized to construct an initial network. Second, a TDCVSS model is utilized to prune the network topology in order to delete false-positive regulatory relationships for each target gene. The complex-valued restricted additive tree and complex-valued differential evolution are proposed to approximate the optimal TDCVSS model. Finally, the overall network could be inferred by integrating the regulations of all target genes. Two real gene expression datasets from E. coli and S. cerevisiae gene networks are utilized to evaluate the performances of our proposed two-step GRN inference algorithm. The results demonstrated that the proposed algorithm could infer GRN more correct than classical methods and time-delayed methods.

Probing the Immune System Dynamics of the COVID-19 Disease for Vaccine Designing and Drug Repurposing Using Bioinformatics Tools

The pathogenesis of COVID-19 is complicated by immune dysfunction. The impact of immune-based therapy in COVID-19 patients has been well documented, with some notable studies on the use of anti-cytokine medicines. However, the complexity of disease phenotypes, patient heterogeneity and the varying quality of evidence from immunotherapy studies provide problems in clinical decision-making. This review seeks to aid therapeutic decision-making by giving an overview of the immunological responses against COVID-19 disease that may contribute to the severity of the disease. We have extensively discussed theranostic methods for COVID-19 detection. With advancements in technology, bioinformatics has taken studies to a higher level. The paper also discusses the application of bioinformatics and machine learning tools for the diagnosis, vaccine design and drug repurposing against SARS-CoV-2.

Control of Cholesterol Metabolism Using a Systems Approach

Simple Summary

Cholesterol is the main sterol in mammals that is essential for healthy cell functionining. It plays a key role in metabolic regulation and signaling, it is a precursor molecule of bile acids, oxysterols, and all steroid hormones. It also contributes to the structural makeup of the membranes. Its homeostasis is tightly controlled since it can harm the body if it is allowed to reach abnormal blood concentrations. One of the diseases associated with elevated cholesterol levels being the major cause of morbidities and mortalities worldwide, is atherosclerosis. In this study, we have developed a model of the cholesterol metabolism taking into account local inflammation and oxidative stress. The aim was to investigate the impact of the interplay of those processes and cholesterol metabolism disturbances on the atherosclerosis development and progression. We have also analyzed the effect of combining different classes of drugs targeting selected components of cholesterol metabolism.

Abstract

Cholesterol is an essential component of mammalian cells and is involved in many fundamental physiological processes; hence, its homeostasis in the body is tightly controlled, and any disturbance has serious consequences. Disruption of the cellular metabolism of cholesterol, accompanied by inflammation and oxidative stress, promotes the formation of atherosclerotic plaques and, consequently, is one of the leading causes of death in the Western world. Therefore, new drugs to regulate disturbed cholesterol metabolism are used and developed, which help to control cholesterol homeostasis but still do not entirely cure atherosclerosis. In this study, a Petri net-based model of human cholesterol metabolism affected by a local inflammation and oxidative stress, has been created and analyzed. The use of knockout of selected pathways allowed us to observe and study the effect of various combinations of commonly used drugs on atherosclerosis. The analysis results led to the conclusion that combination therapy, targeting multiple pathways, may be a fundamental concept in the development of more effective strategies for the treatment and prevention of atherosclerosis.

The future of burn care from a complexity science perspective

Healthcare is undergoing a profound technological and digital transformation and has become increasingly complex. It is important for burns professionals and researchers to adapt to these developments which may require new ways of thinking and subsequent new strategies. As Einstein has put it: 'We must learn to see the world anew'. The relatively new scientific discipline "Complexity science" can give more direction to this and is the metaphorical open door that should not go unnoticed in view of the burn care of the future. Complexity sciences studies 'why the whole is more than the sum of the parts'. It studies how multiple separate components interact with each other and their environment and how these interactions lead to 'behavior of the system'. Biological systems are always part of smaller and larger systems and exhibit the behavior of adaptivity, hence the name complex adaptive systems. From the perspective of complexity science, a severe burn injury is an extreme disruption of the 'human body system'. But this disruption also applies to the systems at the organ and cellular level. All these systems follow principles of complex systems. Awareness of the scaling process at multilevel helps to understand and manage the complex situation when dealing with severe burn cases. The aim of this paper is to create awareness of the concept of complexity and to demonstrate the value and possibilities of complexity science methods and tools for the future of burn care through examples from preclinical, clinical, and organizational perspective in burn care.

BioDynaMo: a modular platform for high-performance agent-based simulation

MOTIVATION

Agent-based modeling is an indispensable tool for studying complex biological systems. However, existing simulation platforms do not always take full advantage of modern hardware and often have a field-specific software design.

RESULTS

We present a novel simulation platform called BioDynaMo that alleviates both of these problems. BioDynaMo features a modular and high-performance simulation engine. We demonstrate that BioDynaMo can be used to simulate use cases in: neuroscience, oncology and epidemiology. For each use case, we validate our findings with experimental data or an analytical solution. Our performance results show that BioDynaMo performs up to three orders of magnitude faster than the state-of-the-art baselines. This improvement makes it feasible to simulate each use case with one billion agents on a single server, showcasing the potential BioDynaMo has for computational biology research.

AVAILABILITY AND IMPLEMENTATION

BioDynaMo is an open-source project under the Apache 2.0 license and is available at www.biodynamo.org. Instructions to reproduce the results are available in the supplementary information.

SUPPLEMENTARY INFORMATION

Available at https://doi.org/10.5281/zenodo.5121618.

Developing an Embedding, Koopman and Autoencoder Technologies-Based Multi-Omics Time Series Predictive Model (EKATP) for Systems Biology research

It is very important for systems biologists to predict the state of the multi-omics time series for disease occurrence and health detection. However, it is difficult to make the prediction due to the high-dimensional, nonlinear and noisy characteristics of the multi-omics time series data. For this reason, this study innovatively proposes an Embedding, Koopman and Autoencoder technologies-based multi-omics time series predictive model (EKATP) to predict the future state of a high-dimensional nonlinear multi-omics time series. We evaluate this EKATP by using a genomics time series with chaotic behavior, a proteomics time series with oscillating behavior and a metabolomics time series with flow behavior. The computational experiments demonstrate that our proposed EKATP can substantially improve the accuracy, robustness and generalizability to predict the future state of a time series for multi-omics data.

Nets-within-nets for modeling emergent patterns in ontogenetic processes

Ontogenesis is the development of an organism from its earliest stage to maturity, including homeostasis maintenance throughout adulthood despite environmental perturbations. Almost all cells of a multicellular organism share the same genomic information. Nevertheless, phenotypic diversity and complex supra-cellular architectures emerge at every level, starting from tissues and organs. This is possible thanks to a robust and dynamic interplay of regulative mechanisms. To study ontogenesis, it is necessary to consider different levels of regulation, both genetic and epigenetic. Each cell undergoes a specific path across a landscape of possible regulative states affecting both its structure and its functions during development. This paper proposes using the Nets-Within-Nets formalism, which combines Petri Nets' simplicity with the capability to represent and simulate the interplay between different layers of regulation connected by non-trivial and context-dependent hierarchical relations. In particular, this work introduces a modeling strategy based on Nets-Within-Nets that can model several critical processes involved in ontogenesis. Moreover, it presents a case study focusing on the first phase of Vulval Precursor Cells specification in C.Elegans. The case study shows that the proposed model can simulate the emergent morphogenetic pattern corresponding to the observed developmental outcome of that phase, in both the physiological case and different mutations. The model presented in the results section is available online at https://github.com/sysbio-polito/NWN_CElegans_VPC_model/.

Reverse engineering gene regulatory network based on complex-valued ordinary differential equation model

BACKGROUND

The growing researches of molecular biology reveal that complex life phenomena have the ability to demonstrating various types of interactions in the level of genomics. To establish the interactions between genes or proteins and understand the intrinsic mechanisms of biological systems have become an urgent need and study hotspot.

RESULTS

In order to forecast gene expression data and identify more accurate gene regulatory network, complex-valued version of ordinary differential equation (CVODE) is proposed in this paper. In order to optimize CVODE model, a complex-valued hybrid evolutionary method based on Grammar-guided genetic programming and complex-valued firefly algorithm is presented.

CONCLUSIONS

When tested on three real gene expression datasets from E. coli and Human Cell, the experiment results suggest that CVODE model could improve 20-50% prediction accuracy of gene expression data, which could also infer more true-positive regulatory relationships and less false-positive regulations than ordinary differential equation.

Multiscale modeling in disease

Multiscale computational modeling aims to connect the complex networks of effects at different length and/or time scales. For example, these networks often include intracellular molecular signaling, crosstalk, and other interactions between neighboring cell populations, and higher levels of emergent phenomena across different regions of tissues and among collections of tissues or organs interacting with each other in the whole body. Recent applications of multiscale modeling across intracellular, cellular, and/or tissue levels are highlighted here. These models incorporated the roles of biochemical and biomechanical modulation in processes that are implicated in the mechanisms of several diseases including fibrosis, joint and bone diseases, respiratory infectious diseases, and cancers.

In Silico Factorial Screening and Optimization of Chitosan Based Gel for Urapidil Loaded Microparticle using Reduced Factorial Design

OBJECTIVE

Literature study revealed the poor mechanical strength of chitosan-based microparticles. Our research aimed at developing sufficient strength of microparticle with a suitable concentration of chitosan and non-ionic surfactants such as poloxamer-188 (pluronic). It also aimed to develop and study the effect of variables for prepared microparticles utilizing insilico screening methodology, such as reduced factorial design, followed by optimization.

METHODS

Preliminary trial batches were prepared with variable concentration of chitosan and poloxamer-188 utilizing cross-linked ion-gelation technique. A 20% w/v sodium citrate solution was used as a cross-linking solution. The resolution-IV of 24-1 reduced factorial design was selected to screen the possible and significant independent variables or factors in the dosage form design. A total number of eight runs were suggested by statistical software and responses were recorded. The responses such as spreadability, pH, viscosity and percentage of drug released at 12 h were considered in the screening study. Based on the result, selected factors were included in the optimization technique, including graphical and numerical methods.

RESULTS

The signified factors based on reduced two-level factorial screening design with randomized subtype, were identified by Half-normal and Pareto chart. Mathematical fitting and analysis were performed by the factorial equation during the optimization process. The validation and fitting of models were suggested and evaluated by p-value, adjusted R2, and predicted R2 values. The significant and non-significant terms were evaluated, followed by finding the optimal concentration and region with yellow color highlighted in an overlay plot. Based on the data obtained by the overlay study, the final formulation batch was prepared and the observed value was found to be pretty much nearer as compared to predicted values. Drug-polymer interaction study included attenuated total reflectance, differential scanning calorimetry, and X-Ray diffraction study.

CONCLUSION

The principal of the study design was based on finding the prefixed set parameter values utilizing the concept of in-silico screening technique and optimization with a minimal number of trials and study expenses. It concluded that Poloxamer-188 (0.94%), chitosan (2.38%), swelling time (1.81 h), and parts of chitosan (78.51%) in a formulation batch would fulfill the predetermined parameter with specific values.

Multi-scale modeling for systematically understanding the key roles of microglia in AD development

Alzheimer's disease (AD) is the leading cause of age-related dementia, affecting over 5 million people in the United States. Unfortunately, current therapies are largely palliative and several potential drug candidates have failed in late-stage clinical trials. Studies suggest that microglia-mediated neuroinflammation might be responsible for the failures of various therapies. Microglia contribute to Aβ clearance in the early stage of neurodegeneration and may contribute to AD development at the late stage by releasing pro-inflammatory cytokines. However, the activation profile and phenotypic changes of microglia during the development of AD are poorly understood. To systematically understand the key role of microglia in AD progression and predict the optimal therapeutic strategy in silico, we developed a 3D multi-scale model of AD (MSMAD) by integrating multi-level experimental data, to manipulate the neurodegeneration in a simulated system. Based on our analysis, we revealed that how TREM2-related signal transduction leads to an imbalance in the activation of different microglia phenotypes, thereby promoting AD development. Our MSMAD model also provides an optimal therapeutic strategy for improving the outcome of AD treatment.

The Optimal Tetralogy of Fallot Repair Using Generative Adversarial Networks

Background

Tetralogy of Fallot (TOF) is a type of congenital cardiac disease with pulmonary artery (PA) stenosis being the most common defect. Repair surgery needs an appropriate patch to enlarge the narrowed artery from the right ventricular (RV) to the PA.

Methods

In this work, we proposed a generative adversarial networks (GANs) based method to optimize the patch size, shape, and location. Firstly, we built the 3D PA of patients by segmentation from cardiac computed tomography angiography. After that, normal and stenotic areas of each PA were detected and labeled into two sub-images groups. Then a GAN was trained based on these sub-images. Finally, an optimal prediction model was utilized to repair the PA with patch augmentation in the new patient.

Results

The fivefold cross-validation (CV) was performed for optimal patch prediction based on GANs in the repair of TOF and the CV accuracy was 93.33%, followed by the clinical outcome. This showed that the GAN model has a significant advantage in finding the best balance point of patch optimization.

Conclusion

This approach has the potential to reduce the intraoperative misjudgment rate, thereby providing a detailed surgical plan in patients with TOF.

Mapping the Metabolic Networks of Tumor Cells and Cancer-Associated Fibroblasts

Metabolism is considered to be the core of all cellular activity. Thus, extensive studies of metabolic processes are ongoing in various fields of biology, including cancer research. Cancer cells are known to adapt their metabolism to sustain high proliferation rates and survive in unfavorable environments with low oxygen and nutrient concentrations. Hence, targeting cancer cell metabolism is a promising therapeutic strategy in cancer research. However, cancers consist not only of genetically altered tumor cells but are interwoven with endothelial cells, immune cells and fibroblasts, which together with the extracellular matrix (ECM) constitute the tumor microenvironment (TME). Cancer-associated fibroblasts (CAFs), which are linked to poor prognosis in different cancer types, are one important component of the TME. CAFs play a significant role in reprogramming the metabolic landscape of tumor cells, but how, and in what manner, this interaction takes place remains rather unclear. This review aims to highlight the metabolic landscape of tumor cells and CAFs, including their recently identified subtypes, in different tumor types. In addition, we discuss various in vitro and in vivo metabolic techniques as well as different in silico computational tools that can be used to identify and characterize CAF-tumor cell interactions. Finally, we provide our view on how mapping the complex metabolic networks of stromal-tumor metabolism will help in finding novel metabolic targets for cancer treatment.

Understanding the Determination of Meat Quality Using Biochemical Characteristics of the Muscle: Stress at Slaughter and Other Missing Keys

Despite increasingly detailed knowledge of the biochemical processes involved in the determination of meat quality traits, robust models, using biochemical characteristics of the muscle to predict future meat quality, lack. The neglecting of various aspects of the model paradigm may explain this. First, preslaughter stress has a major impact on meat quality and varies according to slaughter context and individuals. Yet, it is rarely taken into account in meat quality models. Second, phenotypic similarity does not imply similarity in the underlying biological causes, and several models may be needed to explain a given phenotype. Finally, the implications of the complexity of biological systems are discussed: a homeostatic equilibrium can be reached in countless ways, involving thousands of interacting processes and molecules at different levels of the organism, changing over time and differing between animals. Consequently, even a robust model may explain a significant part, but not all of the variability between individuals.

Tumor cell invasiveness in the initial stages of bladder cancer development - A computational study

Bladder cancer is one of the most common types of cancer, being the sixth more frequent in men, and one with higher recurrence rates and overall treatment costs. We introduce an agent-based computational model of the urothelium, adopting a Cellular Potts Model (CPM) approach to describe both a healthy urothelium and the development of bladder cancer. We focus on the identification of the conditions in which cancer cells cross, by mechanical means, the basement membrane and invade the bladder lamina propria. When within the urothelium the tumor grows in a very constrained environment. These tight conditions imply that the urothelium layer where the tumor initiates greatly determines tumor growth and invasiveness. Moreover, we demonstrate how specific mechanical properties of the cancer cells, as their stiffness or the adhesion to neighboring cells, heavily modulate the critical initial moments of tumor development. We propose that these characteristics should be considered as therapeutic targets to control tumor growth.

Computational study on electromechanics of electroactive hydrogels for cartilage-tissue repair

BACKGROUND AND OBJECTIVE

The self-repair capability of articular cartilage is limited because of non-vascularization and low turnover of its extracellular matrix. Regenerating hyaline cartilage remains a significant clinical challenge as most non-surgical and surgical treatments provide only mid-term relief. Eventually, further pain and mobility loss occur for many patients in the long run due to further joint deterioration. Repair of articular cartilage tissue using electroactive scaffolds and biophysical stimuli like electrical and osmotic stimulation may have the potential to heal cartilage defects occurring due to trauma, osteoarthritis, or sport-related injuries. Therefore, the focus of the current study is to present a computational model of electroactive hydrogels for the cartilage-tissue repair as a first step towards an optimized experimental design.

METHODS

The multiphysics transport model that mainly includes the Poisson-Nernst-Planck equations and the mechanical equation is used to find the electrical stimulation response of the polyelectrolyte hydrogels. Based upon this, a numerical model on electromechanics of electroactive hydrogels seeded with chondrocytes is presented employing the open-source software FEniCS, which is a Python library for finite-element analysis.

RESULTS

We analyzed the ionic concentrations and electric potential in a hydrogel sample and the cell culture medium, the osmotic pressure created due to ionic concentration variations and the resulting hydrogel displacement. The proposed mathematical model was validated with examples from literature.

CONCLUSIONS

The presented model for the electrical and osmotic stimulation of a hydrogel sample can serve as a useful tool for the development and analysis of a cartilaginous scaffold employing electrical stimulation. By analyzing various parameters, we pave the way for future research on a finer scale using open-source software.

Network Analysis of Gut Microbiome and Metabolome to Discover Microbiota-Linked Biomarkers in Patients Affected by Non-Small Cell Lung Cancer

Several studies in recent times have linked gut microbiome (GM) diversity to the pathogenesis of cancer and its role in disease progression through immune response, inflammation and metabolism modulation. This study focused on the use of network analysis and weighted gene co-expression network analysis (WGCNA) to identify the biological interaction between the gut ecosystem and its metabolites that could impact the immunotherapy response in non-small cell lung cancer (NSCLC) patients undergoing second-line treatment with anti-PD1. Metabolomic data were merged with operational taxonomic units (OTUs) from 16S RNA-targeted metagenomics and classified by chemometric models. The traits considered for the analyses were: (i) condition: disease or control (CTRLs), and (ii) treatment: responder (R) or non-responder (NR). Network analysis indicated that indole and its derivatives, aldehydes and alcohols could play a signaling role in GM functionality. WGCNA generated, instead, strong correlations between short-chain fatty acids (SCFAs) and a healthy GM. Furthermore, commensal bacteria such as Akkermansia muciniphila, Rikenellaceae, Bacteroides, Peptostreptococcaceae, Mogibacteriaceae and Clostridiaceae were found to be more abundant in CTRLs than in NSCLC patients. Our preliminary study demonstrates that the discovery of microbiota-linked biomarkers could provide an indication on the road towards personalized management of NSCLC patients.

Advances in Astrocyte Computational Models: From Metabolic Reconstructions to Multi-omic Approaches

The growing importance of astrocytes in the field of neuroscience has led to a greater number of computational models devoted to the study of astrocytic functions and their metabolic interactions with neurons. The modeling of these interactions demands a combined understanding of brain physiology and the development of computational frameworks based on genomic-scale reconstructions, system biology, and dynamic models. These computational approaches have helped to highlight the neuroprotective mechanisms triggered by astrocytes and other glial cells, both under normal conditions and during neurodegenerative processes. In the present review, we evaluate some of the most relevant models of astrocyte metabolism, including genome-scale reconstructions and astrocyte-neuron interactions developed in the last few years. Additionally, we discuss novel strategies from the multi-omics perspective and computational models of other glial cell types that will increase our knowledge in brain metabolism and its association with neurodegenerative diseases.

A New ODE-Based Model for Tumor Cells and Immune System Competition

Changes in diet are heavily associated with high mortality rates in several types of cancer. In this paper, a new mathematical model of tumor cells growth is established to dynamically demonstrate the effects of abnormal cell progression on the cells affected by the tumor in terms of the immune system’s functionality and normal cells’ dynamic growth. This model is called the normal-tumor-immune-unhealthy diet model (NTIUNHDM) and governed by a system of ordinary differential equations. In the NTIUNHDM, there are three main populations normal cells, tumor cell and immune cells. The model is discussed analytically and numerically by utilizing a fourth-order Runge–Kutta method. The dynamic behavior of the NTIUNHDM is discussed by analyzing the stability of the system at various equilibrium points and the Mathematica software is used to simulate the model. From analysis and simulation of the NTIUNHDM, it can be deduced that instability of the response stage, due to a weak immune system, is classified as one of the main reasons for the coexistence of abnormal cells and normal cells. Additionally, it is obvious that the NTIUNHDM has only one stable case when abnormal cells begin progressing into early stages of tumor cells such that the immune cells are generated once. Thus, early boosting of the immune system might contribute to reducing the risk of cancer.

A Role of Inflammation and Immunity in Essential Hypertension-Modeled and Analyzed Using Petri Nets

Recent studies have shown that the innate and adaptive immune system, together with low-grade inflammation, may play an important role in essential hypertension. In this work, to verify the importance of selected factors for the development of essential hypertension, we created a Petri net-based model and analyzed it. The analysis was based mainly on t-invariants, knockouts of selected fragments of the net and its simulations. The blockade of the renin-angiotensin (RAA) system revealed that the most significant effect on the emergence of essential hypertension has RAA activation. This blockade affects: (1) the formation of angiotensin II, (2) inflammatory process (by influencing C-reactive protein (CRP)), (3) the initiation of blood coagulation, (4) bradykinin generation via the kallikrein-kinin system, (5) activation of lymphocytes in hypertension, (6) the participation of TNF alpha in the activation of the acute phase response, and (7) activation of NADPH oxidase-a key enzyme of oxidative stress. On the other hand, we found that the blockade of the activation of the RAA system may not eliminate hypertension that can occur due to disturbances associated with the osmotically independent binding of Na in the interstitium. Moreover, we revealed that inflammation alone is not enough to trigger primary hypertension, but it can coexist with it. We believe that our research may contribute to a better understanding of the pathology of hypertension. It can help identify potential subprocesses, which blocking will allow better control of essential hypertension.

From tumour perfusion to drug delivery and clinical translation of in silico cancer models

In silico cancer models have demonstrated great potential as a tool to improve drug design, optimise the delivery of drugs to target sites in the host tissue and, hence, improve therapeutic efficacy and patient outcome. However, there are significant barriers to the successful translation of in silico technology from bench to bedside. More precisely, the specification of unknown model parameters, the necessity for models to adequately reflect in vivo conditions, and the limited amount of pertinent validation data to evaluate models' accuracy and assess their reliability, pose major obstacles in the path towards their clinical translation. This review aims to capture the state-of-the-art in in silico cancer modelling of vascularised solid tumour growth, and identify the important advances and barriers to success of these models in clinical oncology. Particular emphasis has been put on continuum-based models of cancer since they - amongst the class of mechanistic spatio-temporal modelling approaches - are well-established in simulating transport phenomena and the biomechanics of tissues, and have demonstrated potential for clinical translation. Three important avenues in in silico modelling are considered in this contribution: first, since systemic therapy is a major cancer treatment approach, we start with an overview of the tumour perfusion and angiogenesis in silico models. Next, we present the state-of-the-art in silico work encompassing the delivery of chemotherapeutic agents to cancer nanomedicines through the bloodstream, and then review continuum-based modelling approaches that demonstrate great promise for successful clinical translation. We conclude with a discussion of what we view to be the key challenges and opportunities for in silico modelling in personalised and precision medicine.

Parallelisation strategies for agent based simulation of immune systems

BACKGROUND

In recent years, the study of immune response behaviour using bottom up approach, Agent Based Modeling (ABM), has attracted considerable efforts. The ABM approach is a very common technique in the biological domain due to high demand for a large scale analysis tools for the collection and interpretation of information to solve biological problems. Simulating massive multi-agent systems (i.e. simulations containing a large number of agents/entities) requires major computational effort which is only achievable through the use of parallel computing approaches.

RESULTS

This paper explores different approaches to parallelising the key component of biological and immune system models within an ABM model: pairwise interactions. The focus of this paper is on the performance and algorithmic design choices of cell interactions in continuous and discrete space where agents/entities are competing to interact with one another within a parallel environment.

CONCLUSIONS

Our performance results demonstrate the applicability of these methods to a broader class of biological systems exhibiting typical cell to cell interactions. The advantage and disadvantage of each implementation is discussed showing each can be used as the basis for developing complete immune system models on parallel hardware.

Assembling systems biology, embryo development and teratogenesis: What do we know so far and where to go next?

The recognition of molecular mechanisms of a teratogen can provide insights to understand its embryopathy, and later to plan strategies for the prevention of new exposures. In this context, experimental research is the most invested approach. Despite its relevance, these assays require financial and time investment. Hence, the evaluation of such mechanisms through systems biology rise as an alternative for this conventional methodology. Systems biology is an integrative field that connects experimental and computational analyses, assembling interaction networks between genes, proteins, and even teratogens. It is a valid strategy to generate new hypotheses, that can later be confirmed in experimental assays. Here, we present a literature review of the application of systems biology in embryo development and teratogenesis studies. We provide a glance at the data available in public databases, and evaluate common mechanisms between different teratogens. Finally, we discuss the advantages of using this strategy in future teratogenesis researches.

Systems biology: perspectives on multiscale modeling in research on endocrine-related cancers

Drawing on concepts from experimental biology, computer science, informatics, mathematics and statistics, systems biologists integrate data across diverse platforms and scales of time and space to create computational and mathematical models of the integrative, holistic functions of living systems. Endocrine-related cancers are well suited to study from a systems perspective because of the signaling complexities arising from the roles of growth factors, hormones and their receptors as critical regulators of cancer cell biology and from the interactions among cancer cells, normal cells and signaling molecules in the tumor microenvironment. Moreover, growth factors, hormones and their receptors are often effective targets for therapeutic intervention, such as estrogen biosynthesis, estrogen receptors or HER2 in breast cancer and androgen receptors in prostate cancer. Given the complexity underlying the molecular control networks in these cancers, a simple, intuitive understanding of how endocrine-related cancers respond to therapeutic protocols has proved incomplete and unsatisfactory. Systems biology offers an alternative paradigm for understanding these cancers and their treatment. To correctly interpret the results of systems-based studies requires some knowledge of how in silico models are built, and how they are used to describe a system and to predict the effects of perturbations on system function. In this review, we provide a general perspective on the field of cancer systems biology, and we explore some of the advantages, limitations and pitfalls associated with using predictive multiscale modeling to study endocrine-related cancers.

Sensitivity analysis for reproducible candidate values of model parameters in signaling hub model

Mathematical models for signaling pathways are helpful for understanding molecular mechanism in the pathways and predicting dynamic behavior of the signal activity. To analyze the robustness of such models, local sensitivity analysis has been implemented. However, such analysis primarily focuses on only a certain parameter set, even though diverse parameter sets that can recapitulate experiments may exist. In this study, we performed sensitivity analysis that investigates the features in a system considering the reproducible and multiple candidate values of the model parameters to experiments. The results showed that although different reproducible model parameter values have absolute differences with respect to sensitivity strengths, specific trends of some relative sensitivity strengths exist between reactions regardless of parameter values. It is suggested that (i) network structure considerably influences the relative sensitivity strength and (ii) one might be able to predict relative sensitivity strengths specified in the parameter sets employing only one of the reproducible parameter sets.

Evaluation of anti-EGFR-iRGD recombinant protein with GOLD nanoparticles: synergistic effect on antitumor efficiency using optimized deep neural networks

NP screening through a deep learning approach against Anti-EGFR and validation through docking with AuNP. Biochemical pathway and simulation of AuNP with Anti-EGFR and further implementation in biological circuits.

How scientific literature analysis yields innovative therapeutic hypothesis through integrative iterations

It is becoming generally accepted that the current diagnostic system often guarantees, rather than diminishes, disease heterogeneity. In effects, syndrome-dominated conceptual thinking has become a barrier to understanding the biological causes of complex, multifactorial diseases characterized by clinical and therapeutic heterogeneity. Furthermore, not only is the flood of currently available medical and biological information highly heterogeneous, it is also often conflicting. Together with the entire absence of functional models of pathogenesis and pathological evolution of complex diseases, this leads to a situation where illness activity cannot be coherently approached and where therapeutic developments become highly problematic. Acquisition of the necessary knowledge can be obtained, in parts, using in silico models produced through analytical approaches and processes collectively known as `Systems Biology'. However, without analytical approaches that specifically incorporate the facts that all that is called `information' is not necessarily useful nor utilisable and that all information should be considered as a priori suspect, modelling attempts will fail because of the much too numerous conflicting and, although correct in molecular terms, physiologically invalid reports. In the present essay, we suggest means whereby this body of problems could be functionally attacked and describe new analytical approaches that have demonstrated their efficacy in alleviating these difficulties.

Data-Driven Methods for Advancing Precision Oncology

PURPOSE OF REVIEW

This article discusses the advances, methods, challenges, and future directions of data-driven methods in advancing precision oncology for biomedical research, drug discovery, clinical research, and practice.

RECENT FINDINGS

Precision oncology provides individually tailored cancer treatment by considering an individual's genetic makeup, clinical, environmental, social, and lifestyle information. Challenges include voluminous, heterogeneous, and disparate data generated by different technologies with multiple modalities such as Omics, electronic health records, clinical registries and repositories, medical imaging, demographics, wearables, and sensors. Statistical and machine learning methods have been continuously adapting to the ever-increasing size and complexity of data. Precision Oncology supportive analytics have improved turnaround time in biomarker discovery and time-to-application of new and repurposed drugs. Precision oncology additionally seeks to identify target patient populations based on genomic alterations that are sensitive or resistant to conventional or experimental treatments. Predictive models have been developed for cancer progression and survivorship, drug sensitivity and resistance, and identification of the most suitable combination treatments for individual patient scenarios. In the future, clinical decision support systems need to be revamped to better incorporate knowledge from precision oncology, thus enabling clinical practitioners to provide precision cancer care.

SUMMARY

Open Omics datasets, machine learning algorithms, and predictive models have enabled the advancement of precision oncology. Clinical decision support systems with integrated electronic health record and Omics data are needed to provide data-driven recommendations to assist clinicians in disease prevention, early identification, and individualized treatment. Additionally, as cancer is a constantly evolving disorder, clinical decision systems will need to be continually updated based on more recent knowledge and datasets.