Agent-Based Modeling in Molecular Systems Biology

A three-dimensional lattice-free agent-based model of intracellular ice formation and propagation and intercellular mechanics in liver tissues

A successful cryopreservation of tissues and organs is crucial for medical procedures and drug development acceleration. However, there are only a few instances of successful tissue cryopreservation. One of the main obstacles to successful cryopreservation is intracellular ice damage. Understanding how ice spreads can accelerate protocol development and enable model-based decision-making. Previous models of intracellular ice formation in individual cells have been extended to one-cell-wide arrays to establish the theory of intercellular ice propagation in tissues. The current lattice-based ice propagation models do not account for intercellular forces resulting from cell solidification, which could lead to mechanical disruption of tissue structures during freezing. Moreover, these models have not been expanded to include more realistic tissue architectures. In this article, we discuss the development and validation of a stochastic model for the formation and propagation of ice in small tissues using lattice-free agent-based model. We have improved the existing model by incorporating the mechanical effects of water crystallization within cells. Using information from previous research, we have also created a new model that accounts for ice growth in tissue slabs, spheroids and hepatocyte discs. Our model demonstrates that individual cell freezing can have mechanical consequences and is consistent with earlier findings.

A perfectly imperfect engine: Utilizing the digital twin paradigm in pulmonary hypertension

Pulmonary hypertension (PH) is a severe medical condition with a number of treatment options, the majority of which are introduced without consideration of the underlying mechanisms driving it within an individual and thus a lack of tailored approach to treatment. The one exception is a patient presenting with apparent pulmonary arterial hypertension and shown to have vaso-responsive disease, whose clinical course and prognosis is significantly improved by high dose calcium channel blockers. PH is however characterized by a relative abundance of available data from patient cohorts, ranging from molecular data characterizing gene and protein expression in different tissues to physiological data at the organ level and clinical information. Integrating available data with mechanistic information at the different scales into computational models suggests an approach to a more personalized treatment of the disease using model-based optimization of interventions for individual patients. That is, constructing digital twins of the disease, customized to a patient, promises to be a key technology for personalized medicine, with the aim of optimizing use of existing treatments and developing novel interventions, such as new drugs. This article presents a perspective on this approach in the context of a review of existing computational models for different aspects of the disease, and it lays out a roadmap for a path to realizing it.

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.

Multiscale transport and 4D time-lapse imaging in precision-cut liver slices (PCLS)

Background

Monitoring cellular processes across different levels of complexity, from the cellular to the tissue scale, is important for understanding tissue structure and function. However, it is challenging to monitor and estimate these structural and dynamic interactions within three-dimensional (3D) tissue models.

Objective

The aim of this study was to design a method for imaging, tracking, and quantifying 3D changes in cell morphology (shape and size) within liver tissue, specifically a precision-cut liver slice (PCLS). A PCLS is a 3D model of the liver that allows the study of the structure and function of liver cells in their native microenvironment.

Methods

Here, we present a method for imaging liver tissue during anisosmotic exposure in a multispectral four-dimensional manner. Three metrics of tissue morphology were measured to quantify the effects of osmotic stress on liver tissue. We estimated the changes in the volume of whole precision cut liver slices, quantified the changes in nuclei position, and calculated the changes in volumetric responses of tissue-embedded cells.

Results

During equilibration with cell-membrane-permeating and non-permeating solutes, the whole tissue experiences shrinkage and expansion. As nuclei showed a change in position and directional displacement under osmotic stress, we demonstrate that nuclei could be used as a probe to measure local osmotic and mechanical stress. Moreover, we demonstrate that cells change their volume within tissue slices as a result of osmotic perturbation and that this change in volume is dependent on the position of the cell within the tissue and the duration of the exposure.

Conclusion

The results of this study have implications for a better understanding of multiscale transport, mechanobiology, and triggered biological responses within complex biological structures.

3D animation as a tool for integrative modeling of dynamic molecular mechanisms

As the scientific community accumulates diverse data describing how molecular mechanisms occur, creating and sharing visual models that integrate the richness of this information has become increasingly important to help us explore, refine, and communicate our hypotheses. Three-dimensional (3D) animation is a powerful tool to capture dynamic hypotheses that are otherwise difficult or impossible to visualize using traditional 2D illustration techniques. This perspective discusses the current and future roles that 3D animation can play in the research sphere.

Antibody-dependent enhancement and neutralization against CVB4 investigated in vitro and in silico through an agent-based model

The infection with coxsackievirus B4 (CVB4) can be enhanced in vitro by antibodies directed against the viral capsid protein VP4. In peripheral blood mononuclear cells, antibody-dependent enhancement (ADE) of CVB4 infection leads to the production of interferon alpha (IFN-α). To investigate ADE of CVB4-induced production of IFN-α, an agent-based model was constructed with enhancing and neutralizing antibodies. The model recapitulates viral neutralization and ADE in silico. The enhancing and neutralizing activities of serum samples were evaluated in vitro to confront the model predictions with experimental results. Increasing the incubation time of CVB4 with serum samples improves virus neutralization in silico as well as in vitro. It also results in ADE at lower antibody numbers in silico, which is confirmed in vitro with IFN-α production at lower serum concentrations. Furthermore, incubation of CVB4 with serum at a low temperature does not induce IFN-α production in vitro. Thus, taken together our results suggest that enhancing antibodies bind cryptic epitopes, more accessible with longer incubation time and at higher temperature due to changes in capsid conformation, consistent with previous results indicating that enhancing antibodies are anti-VP4 antibodies.

A multiscale model of the role of microenvironmental factors in cell segregation and heterogeneity in breast cancer development

We analyzed a quantitative multiscale model that describes the epigenetic dynamics during the growth and evolution of an avascular tumor. A gene regulatory network (GRN) formed by a set of ten genes that are believed to play an important role in breast cancer development was kinetically coupled to the microenvironmental agents: glucose, estrogens, and oxygen. The dynamics of spontaneous mutations was described by a Yule-Furry master equation whose solution represents the probability that a given cell in the tissue undergoes a certain number of mutations at a given time. We assumed that the mutation rate is modified by a spatial gradient of nutrients. The tumor mass was simulated by means of cellular automata supplemented with a set of reaction diffusion equations that described the transport of microenvironmental agents. By analyzing the epigenetic state space described by the GRN dynamics, we found three attractors that were identified with cellular epigenetic states: normal, precancer and cancer. For two-dimensional (2D) and three-dimensional (3D) tumors we calculated the spatial distribution of the following quantities: (i) number of mutations, (ii) mutation of each gene and, (iii) phenotypes. Using estrogen as the principal microenvironmental agent that regulates cell proliferation process, we obtained tumor shapes for different values of estrogen consumption and supply rates. It was found that he majority of mutations occurred in cells that were located close to the 2D tumor perimeter or close to the 3D tumor surface. Also, it was found that the occurrence of different phenotypes in the tumor are controlled by estrogen concentration levels since they can change the individual cell threshold and gene expression levels. All results were consistently observed for 2D and 3D tumors.

Fluorescent biosensor imaging meets deterministic mathematical modelling: quantitative investigation of signalling compartmentalization

Cells execute specific responses to diverse environmental cues by encoding information in distinctly compartmentalized biochemical signalling reactions. Genetically encoded fluorescent biosensors enable the spatial and temporal monitoring of signalling events in live cells. Temporal and spatiotemporal computational models can be used to interpret biosensor experiments in complex biochemical networks and to explore hypotheses that are difficult to test experimentally. In this review, we first provide brief discussions of the experimental toolkit of fluorescent biosensors as well as computational basics with a focus on temporal and spatiotemporal deterministic models. We then describe how we used this combined approach to identify and investigate a protein kinase A (PKA) - cAMP - Ca2+ oscillatory circuit in MIN6 β cells, a mouse pancreatic β cell system. We describe the application of this combined approach to interrogate how this oscillatory circuit is differentially regulated in a nano-compartment formed at the plasma membrane by the scaffolding protein A kinase anchoring protein 79/150. We leveraged both temporal and spatiotemporal deterministic models to identify the key regulators of this oscillatory circuit, which we confirmed with further experiments. The powerful approach of combining live-cell biosensor imaging with quantitative modelling, as discussed here, should find widespread use in the investigation of spatiotemporal regulation of cell signalling.

Validation of agent-based models of surface receptor oligomerisation

Receptor dimerisation and higher order oligomerisation regulates signalling by a wide variety of transmembrane receptors. We discuss how agent-based modelling (ABM) combined with advanced microscopy and structural studies can provide new insights into the regulation of clustering, including spatial considerations, revealing novel targets for therapeutic intervention.

Verifiable biology

The formalization of biological systems using computational modelling approaches as an alternative to mathematical-based methods has recently received much interest because computational models provide a deeper mechanistic understanding of biological systems. In particular, formal verification, complementary approach to standard computational techniques such as simulation, is used to validate the system correctness and obtain critical information about system behaviour. In this study, we survey the most frequently used computational modelling approaches and formal verification techniques for computational biology. We compare a number of verification tools and software suites used to analyse biological systems and biochemical networks, and to verify a wide range of biological properties. For users who have no expertise in formal verification, we present a novel methodology that allows them to easily apply formal verification techniques to analyse their biological or biochemical system of interest.

Improving the hole picture: towards a consensus on the mechanism of nuclear transport

Nuclear pore complexes (NPCs) mediate the exchange of materials between the nucleoplasm and cytoplasm, playing a key role in the separation of nucleic acids and proteins into their required compartments. The static structure of the NPC is relatively well defined by recent cryo-EM and other studies. The functional roles of dynamic components in the pore of the NPC, phenylalanyl-glycyl (FG) repeat rich nucleoporins, is less clear because of our limited understanding of highly dynamic protein systems. These proteins form a 'restrained concentrate' which interacts with and concentrates nuclear transport factors (NTRs) to provide facilitated nucleocytoplasmic transport of cargoes. Very rapid on- and off-rates among FG repeats and NTRs supports extremely fast facilitated transport, close to the rate of macromolecular diffusion in cytoplasm, while complexes without specific interactions are entropically excluded, though details on several aspects of the transport mechanism and FG repeat behaviors remain to be resolved. However, as discussed here, new technical approaches combined with more advanced modeling methods will likely provide an improved dynamic description of NPC transport, potentially at the atomic level in the near future. Such advances are likely to be of major benefit in comprehending the roles the malfunctioning NPC plays in cancer, ageing, viral diseases, and neurodegeneration.

DNA storage-from natural biology to synthetic biology

Natural DNA storage allows cellular differentiation, evolution, the growth of our children and controls all our ecosystems. Here, we discuss the fundamental aspects of DNA storage and recent advances in this field, with special emphasis on natural processes and solutions that can be exploited. We point out new ways of efficient DNA and nucleotide storage that are inspired by nature. Within a few years DNA-based information storage may become an attractive and natural complementation to current electronic data storage systems. We discuss rapid and directed access (e.g. DNA elements such as promotors, enhancers), regulatory signals and modulation (e.g. lncRNA) as well as integrated high-density storage and processing modules (e.g. chromosomal territories). There is pragmatic DNA storage for use in biotechnology and human genetics. We examine DNA storage as an approach for synthetic biology (e.g. light-controlled nucleotide processing enzymes). The natural polymers of DNA and RNA offer much for direct storage operations (read-in, read-out, access control). The inbuilt parallelism (many molecules at many places working at the same time) is important for fast processing of information. Using biology concepts from chromosomal storage, nucleic acid processing as well as polymer material sciences such as electronical effects in enzymes, graphene, nanocellulose up to DNA macramé , DNA wires and DNA-based aptamer field effect transistors will open up new applications gradually replacing classical information storage methods in ever more areas over time (decades).

Mathematical modeling of antihypertensive therapy

Hypertension is a multifactorial disease arising from complex pathophysiological pathways. Individual characteristics of patients result in different responses to various classes of antihypertensive medications. Therefore, evaluating the efficacy of therapy based on in silico predictions is an important task. This study is a continuation of research on the modular agent-based model of the cardiovascular and renal systems (presented in the previously published article). In the current work, we included in the model equations simulating the response to antihypertensive therapies with different mechanisms of action. For this, we used the pharmacodynamic effects of the angiotensin II receptor blocker losartan, the calcium channel blocker amlodipine, the angiotensin-converting enzyme inhibitor enalapril, the direct renin inhibitor aliskiren, the thiazide diuretic hydrochlorothiazide, and the β-blocker bisoprolol. We fitted therapy parameters based on known clinical trials for all considered medications, and then tested the model's ability to show reasonable dynamics (expected by clinical observations) after treatment with individual drugs and their dual combinations in a group of virtual patients with hypertension. The extended model paves the way for the next step in personalized medicine that is adapting the model parameters to a real patient and predicting his response to antihypertensive therapy. The model is implemented in the BioUML software and is available at https://gitlab.sirius-web.org/virtual-patient/antihypertensive-treatment-modeling.

The dawn of interoperating spatial models in cell biology

Spatial simulations are becoming an increasingly ubiquitous component in the cycle of discovery, experimentation, and communication across the sciences. In cell biology, many researchers share a vision of developing multiscale models that recapitulate observable behaviors spanning from atoms to cells to tissues. For this dream to become a reality, however, simulation technologies must provide a means for integration and interoperability as they advance. Already, the field has developed numerous methods that span scales of length, time, and complexity to create an extensive body of effective simulation approaches, and although these approaches rarely interoperate, they collectively cover a large spectrum of knowledge that future models may handle in a more unified manner. Here, we discuss the importance of making the data, workflows, and outputs of spatial simulations shareable and interoperable; and how democratization could encourage diverse biologists to participate more easily in developing models to advance our understanding of biological systems.

Agent-based modeling predicts RAC1 is critical for ovarian cancer metastasis

Experimental and computational studies pinpoint rate-limiting step(s) in metastasis governed by Rac1. Using ovarian cancer cell and animal models, Rac1 expression was manipulated, and quantitative measurements of cell-cell and cell-substrate adhesion, cell invasion, mesothelial clearance, and peritoneal tumor growth discriminated the tumor behaviors most highly influenced by Rac1. The experimental data were used to parameterize an agent-based computational model simulating peritoneal niche colonization, intravasation, and hematogenous metastasis to distant organs. Increased ovarian cancer cell survival afforded by the more rapid adhesion and intravasation upon Rac1 overexpression is predicted to increase the numbers of and the rates at which tumor cells are disseminated to distant sites. Surprisingly, crowding of cancer cells along the blood vessel was found to decrease the numbers of cells reaching a distant niche irrespective of Rac1 overexpression or knockdown, suggesting that sites for tumor cell intravasation are rate limiting and become accessible if cells intravasate rapidly or are displaced due to diminished viability. Modeling predictions were confirmed through animal studies of Rac1-dependent metastasis to the lung. Collectively, the experimental and modeling approaches identify cell adhesion, rapid intravasation, and survival in the blood as parameters in the ovarian metastatic cascade that are most critically dependent on Rac1.

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.

Scar formation from the perspective of complexity science: a new look at the biological system as a whole

A burn wound is a complex systemic disease at multiple levels. Current knowledge of scar formation after burn injury has come from traditional biological and clinical studies. These are normally focused on just a small part of the entire process, which has limited our ability to sufficiently understand the underlying mechanisms and to predict systems behaviour. Scar formation after burn injury is a result of a complex biological system-wound healing. It is a part of a larger whole. In this self-organising system, many components form networks of interactions with each other. These networks of interactions are typically non-linear and change their states dynamically, responding to the environment and showing emergent long-term behaviour. How molecular and cellular data relate to clinical phenomena, especially regarding effective therapies of burn wounds to achieve minimal scarring, is difficult to unravel and comprehend. Complexity science can help bridge this gap by integrating small parts into a larger whole, such that relevant biological mechanisms and data are combined in a computational model to better understand the complexity of the entire biological system. A better understanding of the complex biological system of post-burn scar formation could bring research and treatment regimens to the next level. The aim of this review/position paper is to create more awareness of complexity in scar formation after burn injury by describing the basic principles of complexity science and its potential for burn care professionals.

Modeling cholesterol metabolism and atherosclerosis

Atherosclerotic cardiovascular disease (ASCVD) is the leading cause of morbidity and mortality among Western populations. Many risk factors have been identified for ASCVD; however, elevated low-density lipoprotein cholesterol (LDL-C) remains the gold standard. Cholesterol metabolism at the cellular and whole-body level is maintained by an array of interacting components. These regulatory mechanisms have complex behavior. Likewise, the mechanisms which underpin atherogenesis are nontrivial and multifaceted. To help overcome the challenge of investigating these processes mathematical modeling, which is a core constituent of the systems biology paradigm has played a pivotal role in deciphering their dynamics. In so doing models have revealed new insights about the key drivers of ASCVD. The aim of this review is fourfold; to provide an overview of cholesterol metabolism and atherosclerosis, to briefly introduce mathematical approaches used in this field, to critically discuss models of cholesterol metabolism and atherosclerosis, and to highlight areas where mathematical modeling could help to investigate in the future. This article is categorized under: Cardiovascular Diseases > Computational Models.

Characterizing Binding Interactions That Are Essential for Selective Transport through the Nuclear Pore Complex

Specific macromolecules are rapidly transported across the nuclear envelope via the nuclear pore complex (NPC). The selective transport process is facilitated when nuclear transport receptors (NTRs) weakly and transiently bind to intrinsically disordered constituents of the NPC, FG Nups. These two types of proteins help maintain the selective NPC barrier. To interrogate their binding interactions in vitro, we deployed an NPC barrier mimic. We created the stationary phase by covalently attaching fragments of a yeast FG Nup called Nsp1 to glass coverslips. We used a tunable mobile phase containing NTR, nuclear transport factor 2 (NTF2). In the stationary phase, three main factors affected binding: the number of FG repeats, the charge of fragments, and the fragment density. We also identified three main factors affecting binding in the mobile phase: the avidity of the NTF2 variant for Nsp1, the presence of nonspecific proteins, and the presence of additional NTRs. We used both experimentally determined binding parameters and molecular dynamics simulations of Nsp1FG fragments to create an agent-based model. The results suggest that NTF2 binding is negatively cooperative and dependent on the density of Nsp1FG molecules. Our results demonstrate the strengths of combining experimental and physical modeling approaches to study NPC-mediated transport.

Magnesium ions regulate mesenchymal stem cells population and osteogenic differentiation: A fuzzy agent-based modeling approach

Mesenchymal stem cells (MSCs) are proliferative and multipotent cells that play a key role in the bone regeneration process. Empirical data have repeatedly shown the bioregulatory importance of magnesium (Mg) ions in MSC growth and osteogenesis. In this study, we propose an agent-based model to predict the spatiotemporal dynamics of the MSC population and osteogenic differentiation in response to Mg2+ ions. A fuzzy-logic controller was designed to govern the decision-making process of cells by predicting four cellular processes of proliferation, differentiation, migration, and mortality in response to several important bioregulatory factors such as Mg2+ ions, pH, BMP2, and TGF-β1. The model was calibrated using the empirical data obtained from three sets of cell culture experiments. The model successfully reproduced the empirical observations regarding live cell count, viability, DNA content, and the differentiation-related markers of alkaline phosphate (ALP) and osteocalcin (OC). The simulation results, in agreement with the empirical data, showed that Mg2+ ions within 3-6 mM concentration have the highest stimulation effect on cell population growth. The model also correctly reproduced the stimulatory effect of Mg2+ ions on ALP and its inhibitory effect on OC as the early and late differentiation markers, respectively. Besides, the numerical simulation shed light on the innate cellular differences of the cells cultured in different experiments in terms of the proliferative capacity as well as sensitivity to Mg2+ ions. The proposed model can be adopted in the study of the osteogenesis around Mg-based implants where ions released due to degradation interact with local cells and regulate bone regeneration.

Interdisciplinary Approaches to COVID-19

The coronavirus disease 2019 (COVID-19) pandemic has been a significant concern worldwide. The pandemic has demonstrated that public health issues are not merely a health concern but also affect society as a whole. In this chapter, we address the importance of bringing together the world's scientists to find appropriate solutions for controlling and managing the COVID-19 pandemic. Interdisciplinary cooperation, through modern scientific methods, could help to handle the consequences of the pandemic and to avoid the recurrence of future pandemics.

Viral Pandemics of the Last Four Decades: Pathophysiology, Health Impacts and Perspectives

The twenty-first century has witnessed some of the deadliest viral pandemics with far-reaching consequences. These include the Human Immunodeficiency Virus (HIV) (1981), Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV) (2002), Influenza A virus subtype H1N1 (A/H1N1) (2009), Middle East Respiratory Syndrome Coronavirus (MERS-CoV) (2012) and Ebola virus (2013) and the Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) (2019-present). Age- and gender-based characterizations suggest that SARS-CoV-2 resembles SARS-CoV and MERS-CoV with regard tohigher fatality rates in males, and in the older population with comorbidities. The invasion-mechanism of SARS-CoV-2 and SARS-CoV, involves binding of its spike protein with angiotensin-converting enzyme 2 (ACE2) receptors; MERS-CoV utilizes dipeptidyl peptidase 4 (DPP4), whereas H1N1 influenza is equipped with hemagglutinin protein. The viral infections-mediated immunomodulation, and progressive inflammatory state may affect the functions of several other organs. Although no effective commercial vaccine is available for any of the viruses, those against SARS-CoV-2 are being developed at an unprecedented speed. Until now, only Pfizer/BioNTech's vaccine has received temporary authorization from the UK Medicines and Healthcare products Regulatory Agency. Given the frequent emergence of viral pandemics in the 21st century, proper understanding of their characteristics and modes of action are essential to address the immediate and long-term health consequences.

Cancer systems immunology

Tumor immunology is undergoing a renaissance due to the recent profound clinical successes of tumor immunotherapy. These advances have coincided with an exponential growth in the development of -omics technologies. Armed with these technologies and their associated computational and modeling toolsets, systems biologists have turned their attention to tumor immunology in an effort to understand the precise nature and consequences of interactions between tumors and the immune system. Such interactions are inherently multivariate, spanning multiple time and size scales, cell types, and organ systems, rendering systems biology approaches particularly amenable to their interrogation. While in its infancy, the field of 'Cancer Systems Immunology' has already influenced our understanding of tumor immunology and immunotherapy. As the field matures, studies will move beyond descriptive characterizations toward functional investigations of the emergent behavior that govern tumor-immune responses. Thus, Cancer Systems Immunology holds incredible promise to advance our ability to fight this disease.

Integrative Approaches to Cancer Immunotherapy

Cancer immunotherapy aims to arm patients with cancer-fighting immunity. Many new cancer-specific immunotherapeutic drugs have gained approval in the past several years, demonstrating immunotherapy's efficacy and promise as an anticancer modality. Despite these successes, several outstanding questions remain for cancer immunotherapy, including how to make immunotherapy more efficacious in a broader range of cancer types and patients, and how to predict which patients will respond or not respond to therapy. We present a case for integrative systems approaches that will answer these questions. This involves applying mechanistic and statistical modeling, establishing consistent and widely adopted experimental tools to generate systems-level data, and creating sustained mechanisms of support. If implemented, these approaches will lead to major advances in cancer treatment.

The urgent need for integrated science to fight COVID-19 pandemic and beyond

The COVID-19 pandemic has become the leading societal concern. The pandemic has shown that the public health concern is not only a medical problem, but also affects society as a whole; so, it has also become the leading scientific concern. We discuss in this treatise the importance of bringing the world's scientists together to find effective solutions for controlling the pandemic. By applying novel research frameworks, interdisciplinary collaboration promises to manage the pandemic's consequences and prevent recurrences of similar pandemics.

Trends in mathematical modeling of host-pathogen interactions

Pathogenic microorganisms entail enormous problems for humans, livestock, and crop plants. A better understanding of the different infection strategies of the pathogens enables us to derive optimal treatments to mitigate infectious diseases or develop vaccinations preventing the occurrence of infections altogether. In this review, we highlight the current trends in mathematical modeling approaches and related methods used for understanding host-pathogen interactions. Since these interactions can be described on vastly different temporal and spatial scales as well as abstraction levels, a variety of computational and mathematical approaches are presented. Particular emphasis is placed on dynamic optimization, game theory, and spatial modeling, as they are attracting more and more interest in systems biology. Furthermore, these approaches are often combined to illuminate the complexities of the interactions between pathogens and their host. We also discuss the phenomena of molecular mimicry and crypsis as well as the interplay between defense and counter defense. As a conclusion, we provide an overview of method characteristics to assist non-experts in their decision for modeling approaches and interdisciplinary understanding.