Multiscale mathematical modeling to support drug development

Structural and practical identifiability analysis in bioengineering: a beginner's guide

Advancements in digital technology have brought modelling to the forefront in many disciplines from healthcare to architecture. Mathematical models, often represented using parametrised sets of ordinary differential equations, can be used to characterise different processes. To infer possible estimates for the unknown parameters, these models are usually calibrated using associated experimental data. Structural and practical identifiability analyses are a key component that should be assessed prior to parameter estimation. This is because identifiability analyses can provide insights as to whether or not a parameter can take on single, multiple, or even infinitely or countably many values which will ultimately have an impact on the reliability of the parameter estimates. Also, identifiability analyses can help to determine whether the data collected are sufficient or of good enough quality to truly estimate the parameters or if more data or even reparameterization of the model is necessary to proceed with the parameter estimation process. Thus, such analyses also provide an important role in terms of model design (structural identifiability analysis) and the collection of experimental data (practical identifiability analysis). Despite the popularity of using data to estimate the values of unknown parameters, structural and practical identifiability analyses of these models are often overlooked. Possible reasons for non-consideration of application of such analyses may be lack of awareness, accessibility, and usability issues, especially for more complicated models and methods of analysis. The aim of this study is to introduce and perform both structural and practical identifiability analyses in an accessible and informative manner via application to well established and commonly accepted bioengineering models. This will help to improve awareness of the importance of this stage of the modelling process and provide bioengineering researchers with an understanding of how to utilise the insights gained from such analyses in future model development.

Attenuation of Aging-Related Oxidative Stress Pathways by Phytonutrients: A Computational Systems Biology Analysis

Aging results from gradual accumulation of damage to the cellular functions caused by biochemical processes such as oxidative stress, inflammation-driven prolonged cellular senescence state, immune system malfunction, psychological stress, and epigenetic changes due to exposure to environmental toxins. Plant-derived bioactive molecules have been shown to ameliorate the damage from oxidative stress. This research seeks to uncover the mechanisms of action of how phytochemicals from fruit/berry/vegetable (FBV) juice powder mitigate oxidative stress. The study uses a computational systems biology approach to (1) identify biomolecular pathways of oxidative stress; (2) identify phytochemicals from FBV juice powder and their specific action on oxidative stress mechanisms; and (3) quantitatively estimate the effects of FBV juice powder bioactive compounds on oxidative stress. The compounds in FBV affected two oxidative stress molecular pathways: (1) reactive oxygen species (ROS) production and (2) antioxidant enzyme production. Six bioactive compounds including cyanidin, delphinidin, ellagic acid, kaempherol, malvidin, and rutin in FBV significantly lowered production of ROS and increased the production of antioxidant enzymes such as catalase, heme oxygenase-1, superoxide dismutase, and glutathione peroxidase. FBV juice powder provides a combination of bioactive compounds that attenuate aging by affecting multiple pathways of oxidative stress.

Mechanistic Understanding of D-Glucaric Acid to Support Liver Detoxification Essential to Muscle Health Using a Computational Systems Biology Approach

Liver and muscle health are intimately connected. Nutritional strategies that support liver detoxification are beneficial to muscle recovery. Computational-in silico-molecular systems' biology analysis of supplementation of calcium and potassium glucarate salts and their metabolite D-glucaric acid (GA) reveals their positive effect on mitigation of liver detoxification via four specific molecular pathways: (1) ROS production, (2) deconjugation, (3) apoptosis of hepatocytes, and (4) β-glucuronidase synthesis. GA improves liver detoxification by downregulating hepatocyte apoptosis, reducing glucuronide deconjugates levels, reducing ROS production, and inhibiting β-Glucuronidase enzyme that reduces re-absorption of toxins in hepatocytes. Results from this in silico study provide an integrative molecular mechanistic systems explanation for the mitigation of liver toxicity by GA.

Molecular Systems Architecture of Interactome in the Acute Myeloid Leukemia Microenvironment

A molecular systems architecture is presented for acute myeloid leukemia (AML) to provide a framework for organizing the complexity of biomolecular interactions. AML is a multifactorial disease resulting from impaired differentiation and increased proliferation of hematopoietic precursor cells involving genetic mutations, signaling pathways related to the cancer cell genetics, and molecular interactions between the cancer cell and the tumor microenvironment, including endothelial cells, fibroblasts, myeloid-derived suppressor cells, bone marrow stromal cells, and immune cells (e.g., T-regs, T-helper 1 cells, T-helper 17 cells, T-effector cells, natural killer cells, and dendritic cells). This molecular systems architecture provides a layered understanding of intra- and inter-cellular interactions in the AML cancer cell and the cells in the stromal microenvironment. The molecular systems architecture may be utilized for target identification and the discovery of single and combination therapeutics and strategies to treat AML.

Bioactive compounds in green tea may improve transplant tolerance: A computational systems biology analysis

BACKGROUND

Green tea (Camellia sinensis) has bioactive compounds that have been shown to possess nutritive effects on various biomolecular processes such as immunomodulation. This research explores the immunomodulatory effects of green tea in reducing transplant rejection.

METHOD

The study employs computational systems biology: 1) to identify biomolecular mechanisms of immunomodulation in transplant rejection; 2) to identify the bioactive compounds of green tea and their specific effects on mechanisms of immunomodulation in transplant rejection; and, 3) to predict the quantitative effects of those bioactive compounds on immunomodulation in transplant rejection.

RESULTS

Three bioactive compounds of green tea - epicatechin (EC), gallic acid (GA), and epigallocatechin gallate (EGCG), were identified for their potential effects on immunomodulation of transplant rejection. Of the three, EGCG was the only one determined to enhance anti-inflammatory activity by: 1) upregulating synthesis of HO-1 that is known to promote Treg and Th2 phenotypes associated with enabling transplant tolerance; and, 2) downregulating pro-inflammatory cytokines IL-2, IL-17, IFN-γ, TNF-α, NO, IL-6, and IL-1β that are known to promote Th1 and Th17 phenotypes associated with transplant rejection.

CONCLUSIONS

To the best of our knowledge, this study provides the first molecular mechanistic understanding the clinical nutritive value of green tea, specifically the bioactive compound EGCG, in enabling transplant tolerance.

Integrating Multiscale Modeling with Drug Effects for Cancer Treatment

In this paper, we review multiscale modeling for cancer treatment with the incorporation of drug effects from an applied system's pharmacology perspective. Both the classical pharmacology and systems biology are inherently quantitative; however, systems biology focuses more on networks and multi factorial controls over biological processes rather than on drugs and targets in isolation, whereas systems pharmacology has a strong focus on studying drugs with regard to the pharmacokinetic (PK) and pharmacodynamic (PD) relations accompanying drug interactions with multiscale physiology as well as the prediction of dosage-exposure responses and economic potentials of drugs. Thus, it requires multiscale methods to address the need for integrating models from the molecular levels to the cellular, tissue, and organism levels. It is a common belief that tumorigenesis and tumor growth can be best understood and tackled by employing and integrating a multifaceted approach that includes in vivo and in vitro experiments, in silico models, multiscale tumor modeling, continuous/discrete modeling, agent-based modeling, and multiscale modeling with PK/PD drug effect inputs. We provide an example application of multiscale modeling employing stochastic hybrid system for a colon cancer cell line HCT-116 with the application of Lapatinib drug. It is observed that the simulation results are similar to those observed from the setup of the wet-lab experiments at the Translational Genomics Research Institute.

In silico modeling of shear-stress-induced nitric oxide production in endothelial cells through systems biology

Nitric oxide (NO) produced by vascular endothelial cells is a potent vasodilator and an antiinflammatory mediator. Regulating production of endothelial-derived NO is a complex undertaking, involving multiple signaling and genetic pathways that are activated by diverse humoral and biomechanical stimuli. To gain a thorough understanding of the rich diversity of responses observed experimentally, it is necessary to account for an ensemble of these pathways acting simultaneously. In this article, we have assembled four quantitative molecular pathways previously proposed for shear-stress-induced NO production. In these pathways, endothelial NO synthase is activated 1), via calcium release, 2), via phosphorylation reactions, and 3), via enhanced protein expression. To these activation pathways, we have added a fourth, a pathway describing actual NO production from endothelial NO synthase and its various protein partners. These pathways were combined and simulated using CytoSolve, a computational environment for combining independent pathway calculations. The integrated model is able to describe the experimentally observed change in NO production with time after the application of fluid shear stress. This model can also be used to predict the specific effects on the system after interventional pharmacological or genetic changes. Importantly, this model reflects the up-to-date understanding of the NO system, providing a platform upon which information can be aggregated in an additive way.

Computational modeling of chemo-electro-mechanical coupling: a novel implicit monolithic finite element approach

Computational modeling of the human heart allows us to predict how chemical, electrical, and mechanical fields interact throughout a cardiac cycle. Pharmacological treatment of cardiac disease has advanced significantly over the past decades, yet it remains unclear how the local biochemistry of an individual heart cell translates into global cardiac function. Here, we propose a novel, unified strategy to simulate excitable biological systems across three biological scales. To discretize the governing chemical, electrical, and mechanical equations in space, we propose a monolithic finite element scheme. We apply a highly efficient and inherently modular global-local split, in which the deformation and the transmembrane potential are introduced globally as nodal degrees of freedom, whereas the chemical state variables are treated locally as internal variables. To ensure unconditional algorithmic stability, we apply an implicit backward Euler finite difference scheme to discretize the resulting system in time. To increase algorithmic robustness and guarantee optimal quadratic convergence, we suggest an incremental iterative Newton-Raphson scheme. The proposed algorithm allows us to simulate the interaction of chemical, electrical, and mechanical fields during a representative cardiac cycle on a patient-specific geometry, robust and stable, with calculation times on the order of 4 days on a standard desktop computer.