NIAID/SMB Workshop on Multiscale Modeling of Infectious and Immune-Mediated Diseases.

On July 19th, 2023, the National Institute of Allergy and Infectious Diseases co-organized a workshop with the Society of Mathematical Biology, with the authors of this paper as the organizing committee. The workshop, "Bridging multiscale modeling and practical clinical applications in infectious diseases" sought to create an environment for mathematical modelers, statisticians, and infectious disease researchers and clinicians to exchange ideas and perspectives.

SteadyCellPhenotype: a web-based tool for the modeling of biological networks with ternary logic.

SUMMARY: We introduce SteadyCellPhenotype, a browser-based interface for the analysis of ternary biological networks. It includes tools for deterministically finding all steady states of a network, as well as the simulation and visualization of trajectories with publication quality graphics. Simulations allow us to approximate the size of the basin for attractors and deterministic simulations of trajectories nearby specified points allow us to explore the behavior of the system in that neighborhood. AVAILABILITY AND IMPLEMENTATION: https://github.com/knappa/steadycellphenotype MIT License.

COVID-19-associated pulmonary aspergillosis in immunocompetent patients: a virtual patient cohort study.

The opportunistic fungus Aspergillus fumigatus infects the lungs of immunocompromised hosts, including patients undergoing chemotherapy or organ transplantation. More recently however, immunocompetent patients with severe SARS-CoV2 have been reported to be affected by COVID-19 Associated Pulmonary Aspergillosis (CAPA), in the absence of the conventional risk factors for invasive aspergillosis. This paper explores the hypothesis that contributing causes are the destruction of the lung epithelium permitting colonization by opportunistic pathogens. At the same time, the exhaustion of the immune system, characterized by cytokine storms, apoptosis, and depletion of leukocytes may hinder the response to A. fumigatus infection. The combination of these factors may explain the onset of invasive aspergillosis in immunocompetent patients. We used a previously published computational model of the innate immune response to infection with Aspergillus fumigatus. Variation of model parameters was used to create a virtual patient population. A simulation study of this virtual patient population to test potential causes for co-infection in immunocompetent patients. The two most important factors determining the likelihood of CAPA were the inherent virulence of the fungus and the effectiveness of the neutrophil population, as measured by granule half-life and ability to kill fungal cells. Varying these parameters across the virtual patient population generated a realistic distribution of CAPA phenotypes observed in the literature. Computational models are an effective tool for hypothesis generation. Varying model parameters can be used to create a virtual patient population for identifying candidate mechanisms for phenomena observed in actual patient populations.

Aspergillus Utilizes Extracellular Heme as an Iron Source During Invasive Pneumonia, Driving Infection Severity.

BACKGROUND: Depriving microbes of iron is critical to host defense. Hemeproteins, the largest source of iron within vertebrates, are abundant in infected tissues in aspergillosis due to hemorrhage, but Aspergillus species have been thought to lack heme import mechanisms. We hypothesized that heme provides iron to Aspergillus during invasive pneumonia, thereby worsening the outcomes of the infection. METHODS: We assessed the effect of heme on fungal phenotype in various in vitro conditions and in a neutropenic mouse model of invasive pulmonary aspergillosis. RESULTS: In mice with neutropenic invasive aspergillosis, we found a progressive and compartmentalized increase in lung heme iron. Fungal cells cultured under low iron conditions took up heme, resulting in increased fungal iron content, resolution of iron starvation, increased conidiation, and enhanced resistance to oxidative stress. Intrapulmonary administration of heme to mice with neutropenic invasive aspergillosis resulted in markedly increased lung fungal burden, lung injury, and mortality, whereas administration of heme analogs or heme with killed Aspergillus did not. Finally, infection caused by fungal germlings cultured in the presence of heme resulted in a more severe infection. CONCLUSIONS: Invasive aspergillosis induces local hemolysis in infected tissues, thereby supplying heme iron to the fungus, leading to lethal infection.

Mathematical modeling of the Candida albicans yeast to hyphal transition reveals novel control strategies.

Candida albicans, an opportunistic fungal pathogen, is a significant cause of human infections, particularly in immunocompromised individuals. Phenotypic plasticity between two morphological phenotypes, yeast and hyphae, is a key mechanism by which C. albicans can thrive in many microenvironments and cause disease in the host. Understanding the decision points and key driver genes controlling this important transition and how these genes respond to different environmental signals is critical to understanding how C. albicans causes infections in the host. Here we build and analyze a Boolean dynamical model of the C. albicans yeast to hyphal transition, integrating multiple environmental factors and regulatory mechanisms. We validate the model by a systematic comparison to prior experiments, which led to agreement in 17 out of 22 cases. The discrepancies motivate alternative hypotheses that are testable by follow-up experiments. Analysis of this model revealed two time-constrained windows of opportunity that must be met for the complete transition from the yeast to hyphal phenotype, as well as control strategies that can robustly prevent this transition. We experimentally validate two of these control predictions in C. albicans strains lacking the transcription factor UME6 and the histone deacetylase HDA1, respectively. This model will serve as a strong base from which to develop a systems biology understanding of C. albicans morphogenesis.

Topological data analysis in biomedicine: A review.

Significant technological advances made in recent years have shepherded a dramatic increase in utilization of digital technologies for biomedicine- everything from the widespread use of electronic health records to improved medical imaging capabilities and the rising ubiquity of genomic sequencing contribute to a "digitization" of biomedical research and clinical care. With this shift toward computerized tools comes a dramatic increase in the amount of available data, and current tools for data analysis capable of extracting meaningful knowledge from this wealth of information have yet to catch up. This article seeks to provide an overview of emerging mathematical methods with the potential to improve the abilities of clinicians and researchers to analyze biomedical data, but may be hindered from doing so by a lack of conceptual accessibility and awareness in the life sciences research community. In particular, we focus on topological data analysis (TDA), a set of methods grounded in the mathematical field of algebraic topology that seeks to describe and harness features related to the "shape" of data. We aim to make such techniques more approachable to non-mathematicians by providing a conceptual discussion of their theoretical foundations followed by a survey of their published applications to scientific research. Finally, we discuss the limitations of these methods and suggest potential avenues for future work integrating mathematical tools into clinical care and biomedical informatics.

The contribution of microRNA-mediated regulation to short- and long-term gene expression predictability.

MicroRNAs are a class of short, noncoding RNAs which are essential for the coordination and timing of cell differentiation and embryonic development. However, despite their guiding role in development, microRNAs are dysregulated in many pathologies, including nearly all cases of cancer. While both development and oncogenesis can be thought of as extremes of phenotypic plasticity, they characteristically manifest on much different time scales: one taking place over a matter of weeks, the other typically requiring decades. Because microRNAs are believed to support this plasticity, a critically important question is how microRNAs affect phenotypic stability on different time scales, and what dynamical characteristics shift the balance between these two roles. To address this question, we extend a well-established mathematical model of transcriptional gene regulation to include translational regulation by microRNAs, and examine their effects on both short- and long-term gene expression predictability. Our findings show that microRNAs greatly improve short-term predictability for earlier, developmental phenotypes while causing a small decrease in long-term predictability, and that these effects are difficult to separate. In addition to providing a theoretical explanation for this seemingly duplicitous behavior, we describe some of the properties which determine the cost-benefit balance between short-term stabilization and long-term destabilization.

Systems biology of ferroptosis: A modeling approach.

Ferroptosis is a recently discovered form of iron-dependent regulated cell death (RCD) that occurs via peroxidation of phospholipids containing polyunsaturated fatty acid (PUFA) moieties. Activating this form of cell death is an emerging strategy in cancer treatment. Because multiple pathways and molecular species contribute to the ferroptotic process, predicting which tumors will be sensitive to ferroptosis is a challenge. We thus develop a mathematical model of several critical pathways to ferroptosis in order to perform a systems-level analysis of the process. We show that sensitivity to ferroptosis depends on the activity of multiple upstream cascades, including PUFA incorporation into the phospholipid membrane, and the balance between levels of pro-oxidant factors (reactive oxygen species, lipoxogynases) and antioxidant factors (GPX4). We perform a systems-level analysis of ferroptosis sensitivity as an outcome of five input variables (ACSL4, SCD1, ferroportin, transferrin receptor, and p53) and organize the resulting simulations into 'high' and 'low' ferroptosis sensitivity groups. We make a novel prediction corresponding to the combinatorial requirements of ferroptosis sensitivity to SCD1 and ACSL4 activity. To validate our prediction, we model the ferroptotic response of an ovarian cancer stem cell line following single- and double-knockdown of SCD1 and ACSL4. We find that the experimental outcomes are consistent with our simulated predictions. This work suggests that a systems-level approach is beneficial for understanding the complex combined effects of ferroptotic input, and in predicting cancer susceptibility to ferroptosis.

Fostering bioinformatics education through skill development of professors: Big Genomic Data Skills Training for Professors.

Bioinformatics has become an indispensable part of life science over the past 2 decades. However, bioinformatics education is not well integrated at the undergraduate level, especially in liberal arts colleges and regional universities in the United States. One significant obstacle pointed out by the Network for Integrating Bioinformatics into Life Sciences Education is the lack of faculty in the bioinformatics area. Most current life science professors did not acquire bioinformatics analysis skills during their own training. Consequently, a great number of undergraduate and graduate students do not get the chance to learn bioinformatics or computational biology skills within a structured curriculum during their education. To address this gap, we developed a module-based, week-long short course to train small college and regional university professors with essential bioinformatics skills. The bioinformatics modules were built to be adapted by the professor-trainees afterward and used in their own classes. All the course materials can be accessed at https://github.com/TheJacksonLaboratory/JAXBD2K-ShortCourse.

An important role for periplasmic storage in Pseudomonas aeruginosa copper homeostasis revealed by a combined experimental and computational modeling study.

Biological systems require precise copper homeostasis enabling metallation of cuproproteins while preventing metal toxicity. In bacteria, sensing, transport, and storage molecules act in coordination to fulfill these roles. However, there is not yet a kinetic schema explaining the system integration. Here, we report a model emerging from experimental and computational approaches that describes the dynamics of copper distribution in Pseudomonas aeruginosa. Based on copper uptake experiments, a minimal kinetic model describes well the copper distribution in the wild-type bacteria but is unable to explain the behavior of the mutant strain lacking CopA1, a key Cu+ efflux ATPase. The model was expanded through an iterative hypothesis-driven approach, arriving to a mechanism that considers the induction of compartmental pools and the parallel function of CopA and Cus efflux systems. Model simulations support the presence of a periplasmic copper storage with a crucial role under dyshomeostasis conditions in P. aeruginosa. Importantly, the model predicts not only the interplay of periplasmic and cytoplasmic pools but also the existence of a threshold in the concentration of external copper beyond which cells lose their ability to control copper levels.

Control of Intracellular Molecular Networks Using Algebraic Methods.

Many problems in biology and medicine have a control component. Often, the goal might be to modify intracellular networks, such as gene regulatory networks or signaling networks, in order for cells to achieve a certain phenotype, what happens in cancer. If the network is represented by a mathematical model for which mathematical control approaches are available, such as systems of ordinary differential equations, then this problem might be solved systematically. Such approaches are available for some other model types, such as Boolean networks, where structure-based approaches have been developed, as well as stable motif techniques. However, increasingly many published discrete models are mixed-state or multistate, that is, some or all variables have more than two states, and thus the development of control strategies for multistate networks is needed. This paper presents a control approach broadly applicable to general multistate models based on encoding them as polynomial dynamical systems over a finite algebraic state set, and using computational algebra for finding appropriate intervention strategies. To demonstrate the feasibility and applicability of this method, we apply it to a recently developed multistate intracellular model of E2F-mediated bladder cancerous growth and to a model linking intracellular iron metabolism and oncogenic pathways. The control strategies identified for these published models are novel in some cases and represent new hypotheses, or are supported by the literature in others as potential drug targets. Our Macaulay2 scripts to find control strategies are publicly available through GitHub at https://github.com/luissv7/multistatepdscontrol.

Activated Oncogenic Pathway Modifies Iron Network in Breast Epithelial Cells: A Dynamic Modeling Perspective.

Dysregulation of iron metabolism in cancer is well documented and it has been suggested that there is interdependence between excess iron and increased cancer incidence and progression. In an effort to better understand the linkages between iron metabolism and breast cancer, a predictive mathematical model of an expanded iron homeostasis pathway was constructed that includes species involved in iron utilization, oxidative stress response and oncogenic pathways. The model leads to three predictions. The first is that overexpression of iron regulatory protein 2 (IRP2) recapitulates many aspects of the alterations in free iron and iron-related proteins in cancer cells without affecting the oxidative stress response or the oncogenic pathways included in the model. This prediction was validated by experimentation. The second prediction is that iron-related proteins are dramatically affected by mitochondrial ferritin overexpression. This prediction was validated by results in the pertinent literature not used for model construction. The third prediction is that oncogenic Ras pathways contribute to altered iron homeostasis in cancer cells. This prediction was validated by a combination of simulation experiments of Ras overexpression and catalase knockout in conjunction with the literature. The model successfully captures key aspects of iron metabolism in breast cancer cells and provides a framework upon which more detailed models can be built.

Costimulation Endows Immunotherapeutic CD8 T Cells with IL-36 Responsiveness during Aerobic Glycolysis.

CD134- and CD137-primed CD8 T cells mount powerful effector responses upon recall, but even without recall these dual-costimulated T cells respond to signal 3 cytokines such as IL-12. We searched for alternative signal 3 receptor pathways and found the IL-1 family member IL-36R. Although IL-36 alone did not stimulate effector CD8 T cells, in combination with IL-12, or more surprisingly IL-2, it induced striking and rapid TCR-independent IFN-γ synthesis. To understand how signal 3 responses functioned in dual-costimulated T cells we showed that IL-2 induced IL-36R gene expression in a JAK/STAT-dependent manner. These data help delineate a sequential stimulation process where IL-2 conditioning must precede IL-36 for IFN-γ synthesis. Importantly, this responsive state was transient and functioned only in effector T cells capable of aerobic glycolysis. Specifically, as the effector T cells metabolized glucose and consumed O2, they also retained potential to respond through IL-36R. This suggests that T cells use innate receptor pathways such as the IL-36R/axis when programmed for aerobic glycolysis. To explore a function for IL-36R in vivo, we showed that dual costimulation therapy reduced B16 melanoma tumor growth while increasing IL-36R gene expression. In summary, cytokine therapy to eliminate tumors may target effector T cells, even outside of TCR specificity, as long as the effectors are in the correct metabolic state.

AlgoRun: a Docker-based packaging system for platform-agnostic implemented algorithms.

MOTIVATION: There is a growing need in bioinformatics for easy-to-use software implementations of algorithms that are usable across platforms. At the same time, reproducibility of computational results is critical and often a challenge due to source code changes over time and dependencies. RESULTS: The approach introduced in this paper addresses both of these needs with AlgoRun, a dedicated packaging system for implemented algorithms, using Docker technology. Implemented algorithms, packaged with AlgoRun, can be executed through a user-friendly interface directly from a web browser or via a standardized RESTful web API to allow easy integration into more complex workflows. The packaged algorithm includes the entire software execution environment, thereby eliminating the common problem of software dependencies and the irreproducibility of computations over time. AlgoRun-packaged algorithms can be published on http://algorun.org, a centralized searchable directory to find existing AlgoRun-packaged algorithms. AVAILABILITY AND IMPLEMENTATION: AlgoRun is available at http://algorun.org and the source code under GPL license is available at https://github.com/algorun CONTACT: laubenbacher@uchc.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Optimal harvesting for a predator-prey agent-based model using difference equations.

In this paper, a method known as Pareto optimization is applied in the solution of a multi-objective optimization problem. The system in question is an agent-based model (ABM) wherein global dynamics emerge from local interactions. A system of discrete mathematical equations is formulated in order to capture the dynamics of the ABM; while the original model is built up analytically from the rules of the model, the paper shows how minor changes to the ABM rule set can have a substantial effect on model dynamics. To address this issue, we introduce parameters into the equation model that track such changes. The equation model is amenable to mathematical theory­we show how stability analysis can be performed and validated using ABM data. We then reduce the equation model to a simpler version and implement changes to allow controls from the ABM to be tested using the equations. Cohen's weighted κ is proposed as a measure of similarity between the equation model and the ABM, particularly with respect to the optimization problem. The reduced equation model is used to solve a multi-objective optimization problem via a technique known as Pareto optimization, a heuristic evolutionary algorithm. Results show that the equation model is a good fit for ABM data; Pareto optimization provides a suite of solutions to the multi-objective optimization problem that can be implemented directly in the ABM.

Steady state analysis of Boolean molecular network models via model reduction and computational algebra.

BACKGROUND: A key problem in the analysis of mathematical models of molecular networks is the determination of their steady states. The present paper addresses this problem for Boolean network models, an increasingly popular modeling paradigm for networks lacking detailed kinetic information. For small models, the problem can be solved by exhaustive enumeration of all state transitions. But for larger models this is not feasible, since the size of the phase space grows exponentially with the dimension of the network. The dimension of published models is growing to over 100, so that efficient methods for steady state determination are essential. Several methods have been proposed for large networks, some of them heuristic. While these methods represent a substantial improvement in scalability over exhaustive enumeration, the problem for large networks is still unsolved in general. RESULTS: This paper presents an algorithm that consists of two main parts. The first is a graph theoretic reduction of the wiring diagram of the network, while preserving all information about steady states. The second part formulates the determination of all steady states of a Boolean network as a problem of finding all solutions to a system of polynomial equations over the finite number system with two elements. This problem can be solved with existing computer algebra software. This algorithm compares favorably with several existing algorithms for steady state determination. One advantage is that it is not heuristic or reliant on sampling, but rather determines algorithmically and exactly all steady states of a Boolean network. The code for the algorithm, as well as the test suite of benchmark networks, is available upon request from the corresponding author. CONCLUSIONS: The algorithm presented in this paper reliably determines all steady states of sparse Boolean networks with up to 1000 nodes. The algorithm is effective at analyzing virtually all published models even those of moderate connectivity. The problem for large Boolean networks with high average connectivity remains an open problem.

A systems biology approach to iron metabolism.

Iron is critical to the survival of almost all living organisms. However, inappropriately low or high levels of iron are detrimental and contribute to a wide range of diseases. Recent advances in the study of iron metabolism have revealed multiple intricate pathways that are essential to the maintenance of iron homeostasis. Further, iron regulation involves processes at several scales, ranging from the subcellular to the organismal. This complexity makes a systems biology approach crucial, with its enabling technology of computational models based on a mathematical description of regulatory systems. Systems biology may represent a new strategy for understanding imbalances in iron metabolism and their underlying causes.

Control of Medical Digital Twins with Artificial Neural Networks.

The objective of personalized medicine is to tailor interventions to an individual patient's unique characteristics. A key technology for this purpose involves medical digital twins, computational models of human biology that can be personalized and dynamically updated to incorporate patient-specific data collected over time. Certain aspects of human biology, such as the immune system, are not easily captured with physics-based models, such as differential equations. Instead, they are often multi-scale, stochastic, and hybrid. This poses a challenge to existing model-based control and optimization approaches that cannot be readily applied to such models. Recent advances in automatic differentiation and neural-network control methods hold promise in addressing complex control problems. However, the application of these approaches to biomedical systems is still in its early stages. This work introduces dynamics-informed neural-network controllers as an alternative approach to control of medical digital twins. As a first use case for this method, the focus is on agent-based models, a versatile and increasingly common modeling platform in biomedicine. The effectiveness of the proposed neural-network control method is illustrated and benchmarked against other methods with two widely-used agent-based model types. The relevance of the method introduced here extends beyond medical digital twins to other complex dynamical systems.