The virtual cell--a candidate co-ordinator for 'middle-out' modelling of biological systems

Towards early diagnosis of Alzheimer's disease: advances in immune-related blood biomarkers and computational approaches

Alzheimer's disease has an increasing prevalence in the population world-wide, yet current diagnostic methods based on recommended biomarkers are only available in specialized clinics. Due to these circumstances, Alzheimer's disease is usually diagnosed late, which contrasts with the currently available treatment options that are only effective for patients at an early stage. Blood-based biomarkers could fill in the gap of easily accessible and low-cost methods for early diagnosis of the disease. In particular, immune-based blood-biomarkers might be a promising option, given the recently discovered cross-talk of immune cells of the central nervous system with those in the peripheral immune system. Here, we give a background on recent advances in research on brain-immune system cross-talk in Alzheimer's disease and review machine learning approaches, which can combine multiple biomarkers with further information (e.g. age, sex, APOE genotype) into predictive models supporting an earlier diagnosis. In addition, mechanistic modeling approaches, such as agent-based modeling open the possibility to model and analyze cell dynamics over time. This review aims to provide an overview of the current state of immune-system related blood-based biomarkers and their potential for the early diagnosis of Alzheimer's disease.

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.

Development and Analysis of Multiscale Models for Tuberculosis: From Molecules to Populations

Although infectious disease dynamics are often analyzed at the macro-scale, increasing numbers of drug-resistant infections highlight the importance of within-host modeling that simultaneously solves across multiple scales to effectively respond to epidemics. We review multiscale modeling approaches for complex, interconnected biological systems and discuss critical steps involved in building, analyzing, and applying such models within the discipline of model credibility. We also present our two tools: CaliPro, for calibrating multiscale models (MSMs) to datasets, and tunable resolution, for fine- and coarse-graining sub-models while retaining insights. We include as an example our work simulating infection with Mycobacterium tuberculosis to demonstrate modeling choices and how predictions are made to generate new insights and test interventions. We discuss some of the current challenges of incorporating novel datasets, rigorously training computational biologists, and increasing the reach of MSMs. We also offer several promising future research directions of incorporating within-host dynamics into applications ranging from combinatorial treatment to epidemic response.

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.

Rules of Engagement: A Guide to Developing Agent-Based Models

Agent-based models (ABM), also called individual-based models, first appeared several decades ago with the promise of nearly real-time simulations of active, autonomous individuals such as animals or objects. The goal of ABMs is to represent a population of individuals (agents) interacting with one another and their environment. Because of their flexible framework, ABMs have been widely applied to study systems in engineering, economics, ecology, and biology. This chapter is intended to guide the users in the development of an agent-based model by discussing conceptual issues, implementation, and pitfalls of ABMs from first principles. As a case study, we consider an ABM of the multi-scale dynamics of cellular interactions in a microbial community. We develop a lattice-free agent-based model of individual cells whose actions of growth, movement, and division are influenced by both their individual processes (cell cycle) and their contact with other cells (adhesion and contact inhibition).

Modeling Cancer Dynamics

Cancer is a complex disease that continues to pose formidable challenges to therapeutic interventions. An increased understanding of cancer complexity and in particular tumor growth dynamics is critical to the development of more effective therapies. In this respect, an agent-based model of tumor growth is explored with the consideration of cancer reprograming of metabolism and the immune response, to seek insight about the coupling between these two key determinants of tumor growth dynamics. Ultimately, this exploration is intended to inform the development of therapies that can induce a more effective immune response despite the metabolic constraints of the tumor microenvironment.

Modelling cell surface dynamics and cell-cell interactions using Cell Studio: a three-dimensional visualization tool based on gaming technology

Predictive modelling of complex biological systems and biophysical interactions requires the inclusion of multiple nano- and micro-scale events. In many scenarios, however, numerical solutions alone do not necessarily enhance the understanding of the system. Instead, this work explores the use of an agent-based model with visualization capabilities to elucidate interactions between single cells. We present a model of juxtacrine signalling, using Cell Studio, an agent-based modelling system, based on gaming and three-dimensional visualization tools. The main advantages of the system are its ability to apply any cell geometry and to dynamically visualize the diffusion and interactions of the molecules within the cells in real time. These provide an excellent tool for obtaining insight about different biological scenarios, as the user may view the dynamics of a system and observe its emergent behaviour as it unfolds. The agent-based model was validated against the results of a mean-field model of Notch receptors and ligands in two neighbouring cells. The conversion to an agent-based model is described in detail. To demonstrate the advantages of the model, we further created a filopodium-mediated signalling model. Our model revealed that diffusion and endocytosis alone are insufficient to produce significant signalling in a filopodia scenario. This is due to the bottleneck at the cell-filopodium contact region and the long distance to the end of the filopodium. However, allowing active transport of ligands into filopodia enhances the signalling significantly compared with a face-to-face scenario. We conclude that the agent-based approach can provide insights into mechanisms underlying cell signalling. The open-source model can be found in the Internet hosting service GitHub.

An agent-based model of the Notch signaling pathway elucidates three levels of complexity in the determination of developmental patterning

BACKGROUND

The Notch signaling pathway is involved in cell fate decision and developmental patterning in diverse organisms. A receptor molecule, Notch (N), and a ligand molecule (in this case Delta or Dl) are the central molecules in this pathway. In early Drosophila embryos, these molecules determine neural vs. skin fates in a reproducible rosette pattern.

RESULTS

We have created an agent-based model (ABM) that simulates the molecular components for this signaling pathway as agents acting within a spatial representation of a cell. The model captures the changing levels of these components, their transition from one state to another, and their movement from the nucleus to the cell membrane and back to the nucleus again. The model introduces stochastic variation into the system using a random generator within the Netlogo programming environment. The model uses these representations to understand the biological systems at three levels: individual cell fate, the interactions between cells, and the formation of pattern across the system. Using a set of assessment tools, we show that the current model accurately reproduces the rosette pattern of neurons and skin cells in the system over a wide set of parameters. Oscillations in the level of the N agent eventually stabilize cell fate into this pattern. We found that the dynamic timing and the availability of the N and Dl agents in neighboring cells are central to the formation of a correct and stable pattern. A feedback loop to the production of both components is necessary for a correct and stable pattern.

CONCLUSIONS

The signaling pathways within and between cells in our model interact in real time to create a spatially correct field of neurons and skin cells. This model predicts that cells with high N and low Dl drive the formation of the pattern. This model also be used to elucidate general rules of biological self-patterning and decision-making.

Systems biology primer: the basic methods and approaches

Systems biology is an integrative discipline connecting the molecular components within a single biological scale and also among different scales (e.g. cells, tissues and organ systems) to physiological functions and organismal phenotypes through quantitative reasoning, computational models and high-throughput experimental technologies. Systems biology uses a wide range of quantitative experimental and computational methodologies to decode information flow from genes, proteins and other subcellular components of signaling, regulatory and functional pathways to control cell, tissue, organ and organismal level functions. The computational methods used in systems biology provide systems-level insights to understand interactions and dynamics at various scales, within cells, tissues, organs and organisms. In recent years, the systems biology framework has enabled research in quantitative and systems pharmacology and precision medicine for complex diseases. Here, we present a brief overview of current experimental and computational methods used in systems biology.

Cell studio: A platform for interactive, 3D graphical simulation of immunological processes

The field of computer modeling and simulation of biological systems is rapidly advancing, backed by significant progress in the fields of experimentation techniques, computer hardware, and programming software. The result of a simulation may be delivered in several ways, from numerical results, through graphs of the simulated run, to a visualization of the simulation. The vision of an in-silico experiment mimicking an in-vitro or in-vivo experiment as it is viewed under a microscope is appealing but technically demanding and computationally intensive. Here, we report “Cell Studio,” a generic, hybrid platform to simulate an immune microenvironment with biological and biophysical rules. We use game engines—generic programs for game creation which offer ready-made assets and tools—to create a visualized, interactive 3D simulation. We also utilize a scalable architecture that delegates the computational load to a server. The user may view the simulation, move the “camera” around, stop, fast-forward, and rewind it and inject soluble molecules into the extracellular medium at any point in time. During simulation, graphs are created in real time for a broad view of system-wide processes. The model is parametrized using a user-friendly Graphical User Interface (GUI). We show a simple validation simulation and compare its results with those from a “classical” simulation, validated against a “wet” experiment. We believe that interactive, real-time 3D visualization may aid in generating insights from the model and encourage intuition about the immunological scenario.

A cell-based mechanical model of coronary artery tunica media

A three-dimensional cell-based mechanical model of coronary artery tunica media is proposed. The model is composed of spherical cells forming a hexagonal close-packed lattice. Tissue anisotropy is taken into account by varying interaction forces with the direction of intercellular connection. Several cell-centre interaction potentials for repulsion and attraction are considered, including the Hertz contact model and its neo-Hookean extension, the Johnson-Kendall-Roberts model of adhesive contact, and a wormlike chain model. The model is validated against data from in vitro uni-axial tension tests performed on dissected strips of tunica media. The wormlike chain potential in combination with the neo-Hookean Hertz contact model produces stress-stretch curves which represent the experimental data very well.

An agent-based model of anoikis in the colon crypt displays novel emergent behaviour consistent with biological observations

Colorectal cancer (CRC) is a major cause of cancer mortality. Colon crypts are multi-cellular flask-shaped invaginations of the colonic epithelium, with stem cells at their base which support the continual turnover of the epithelium with loss of cells by anoikis from the flat mucosa. Mutations in these stem cells can become embedded in the crypts, a process that is strongly implicated in CRC initiation. We describe a computational model which includes novel features, including an accurate representation of the geometry of the crypt mouth. Model simulations yield previously unseen emergent phenomena, such as localization of cell death to a small region of the crypt mouth which corresponds with that observed in vivo. A mechanism emerges in the model for regulation of crypt cellularity in response to changes in either cell proliferation rates or membrane adhesion strengths. We show that cell shape assumptions influence this behaviour, with cylinders recapitulating biology better than spheres. Potential applications of the model include determination of roles of mutations in neoplasia and exploring factors for altered crypt morphodynamics.

Comparing individual-based approaches to modelling the self-organization of multicellular tissues

The coordinated behaviour of populations of cells plays a central role in tissue growth and renewal. Cells react to their microenvironment by modulating processes such as movement, growth and proliferation, and signalling. Alongside experimental studies, computational models offer a useful means by which to investigate these processes. To this end a variety of cell-based modelling approaches have been developed, ranging from lattice-based cellular automata to lattice-free models that treat cells as point-like particles or extended shapes. However, it remains unclear how these approaches compare when applied to the same biological problem, and what differences in behaviour are due to different model assumptions and abstractions. Here, we exploit the availability of an implementation of five popular cell-based modelling approaches within a consistent computational framework, Chaste (http://www.cs.ox.ac.uk/chaste). This framework allows one to easily change constitutive assumptions within these models. In each case we provide full details of all technical aspects of our model implementations. We compare model implementations using four case studies, chosen to reflect the key cellular processes of proliferation, adhesion, and short- and long-range signalling. These case studies demonstrate the applicability of each model and provide a guide for model usage.

Towards the virtual artery: a multiscale model for vascular physiology at the physics-chemistry-biology interface

This discussion paper introduces the concept of the Virtual Artery as a multiscale model for arterial physiology and pathologies at the physics-chemistry-biology (PCB) interface. The cellular level is identified as the mesoscopic level, and we argue that by coupling cell-based models with other relevant models on the macro- and microscale, a versatile model of arterial health and disease can be composed. We review the necessary ingredients, both models of arteries at many different scales, as well as generic methods to compose multiscale models. Next, we discuss how this can be combined into the virtual artery. Finally, we argue that the concept of models at the PCB interface could or perhaps should become a powerful paradigm, not only as in our case for studying physiology, but also for many other systems that have such PCB interfaces.This article is part of the themed issue 'Multiscale modelling at the physics-chemistry-biology interface'.

A hybrid agent-based model of the developing mammary terminal end bud

Mammary gland ductal elongation is spearheaded by terminal end buds (TEBs), where populations of highly proliferative cells are maintained throughout post-pubertal organogenesis in virgin mice until the mammary fat pad is filled by a mature ductal tree. We have developed a hybrid multiscale agent-based model to study how cellular differentiation pathways, cellular proliferation capacity, and endocrine and paracrine signaling play a role during development of the mammary gland. A simplified cellular phenotypic hierarchy that includes stem, progenitor, and fully differentiated cells within the TEB was implemented. Model analysis finds that mammary gland development was highly sensitive to proliferation events within the TEB, with progenitors likely undergoing 2-3 proliferation cycles before transitioning to a non-proliferative phenotype, and this result is in agreement with our previous experimental work. Endocrine and paracrine signaling were found to provide reliable ductal elongation rate regulation, while variations in the probability a new daughter cell will be of a proliferative phenotype were seen to have minimal effects on ductal elongation rates. Moreover, the distribution of cellular phenotypes within the TEB was highly heterogeneous, demonstrating significant allowable plasticity in possible phenotypic distributions while maintaining biologically relevant growth behavior. Finally, simulation results indicate ductal elongation rates due to cellular proliferation within the TEB may have a greater sensitivity to upstream endocrine signaling than endothelial to stromal paracrine signaling within the TEB. This model provides a useful tool to gain quantitative insights into cellular population dynamics and the effects of endocrine and paracrine signaling within the pubertal terminal end bud.

Multi-scale Modeling of the Cardiovascular System: Disease Development, Progression, and Clinical Intervention

Cardiovascular diseases (CVDs) are the leading cause of death in the western world. With the current development of clinical diagnostics to more accurately measure the extent and specifics of CVDs, a laudable goal is a better understanding of the structure-function relation in the cardiovascular system. Much of this fundamental understanding comes from the development and study of models that integrate biology, medicine, imaging, and biomechanics. Information from these models provides guidance for developing diagnostics, and implementation of these diagnostics to the clinical setting, in turn, provides data for refining the models. In this review, we introduce multi-scale and multi-physical models for understanding disease development, progression, and designing clinical interventions. We begin with multi-scale models of cardiac electrophysiology and mechanics for diagnosis, clinical decision support, personalized and precision medicine in cardiology with examples in arrhythmia and heart failure. We then introduce computational models of vasculature mechanics and associated mechanical forces for understanding vascular disease progression, designing clinical interventions, and elucidating mechanisms that underlie diverse vascular conditions. We conclude with a discussion of barriers that must be overcome to provide enhanced insights, predictions, and decisions in pre-clinical and clinical applications.

Sources of variability in cytosolic calcium transients triggered by stimulation of homogeneous uro-epithelial cell monolayers

Epithelial tissue structure is the emergent outcome of the interactions between large numbers of individual cells. Experimental cell biology offers an important tool to unravel these complex interactions, but current methods of analysis tend to be limited to mean field approaches or representation by selected subsets of cells. This may result in bias towards cells that respond in a particular way and/or neglect local, context-specific cell responses. Here, an automated algorithm was applied to examine in detail the individual calcium transients evoked in genetically homogeneous, but asynchronous populations of cultured non-immortalized normal human urothelial cells when subjected to either the global application of an external agonist or a localized scratch wound. The recorded calcium transients were classified automatically according to a set of defined metrics and distinct sub-populations of cells that responded in qualitatively different ways were observed. The nature of this variability in the homogeneous cell population was apportioned to two sources: intrinsic variation in individual cell responses and extrinsic variability due to context-specific factors of the environment, such as spatial heterogeneity. Statistically significant variation in the features of the calcium transients evoked by scratch wounding according to proximity to the wound edge was identified. The manifestation of distinct sub-populations of cells is considered central to the coordination of population-level response resulting in wound closure.

Introduction of a Framework for Dynamic Knowledge Representation of the Control Structure of Transplant Immunology: Employing the Power of Abstraction with a Solid Organ Transplant Agent-Based Model

Agent-based modeling has been used to characterize the nested control loops and non-linear dynamics associated with inflammatory and immune responses, particularly as a means of visualizing putative mechanistic hypotheses. This process is termed dynamic knowledge representation and serves a critical role in facilitating the ability to test and potentially falsify hypotheses in the current data- and hypothesis-rich biomedical research environment. Importantly, dynamic computational modeling aids in identifying useful abstractions, a fundamental scientific principle that pervades the physical sciences. Recognizing the critical scientific role of abstraction provides an intellectual and methodological counterweight to the tendency in biology to emphasize comprehensive description as the primary manifestation of biological knowledge. Transplant immunology represents yet another example of the challenge of identifying sufficient understanding of the inflammatory/immune response in order to develop and refine clinically effective interventions. Advances in immunosuppressive therapies have greatly improved solid organ transplant (SOT) outcomes, most notably by reducing and treating acute rejection. The end goal of these transplant immune strategies is to facilitate effective control of the balance between regulatory T cells and the effector/cytotoxic T-cell populations in order to generate, and ideally maintain, a tolerant phenotype. Characterizing the dynamics of immune cell populations and the interactive feedback loops that lead to graft rejection or tolerance is extremely challenging, but is necessary if rational modulation to induce transplant tolerance is to be accomplished. Herein is presented the solid organ agent-based model (SOTABM) as an initial example of an agent-based model (ABM) that abstractly reproduces the cellular and molecular components of the immune response to SOT. Despite its abstract nature, the SOTABM is able to qualitatively reproduce acute rejection and the suppression of acute rejection by immunosuppression to generate transplant tolerance. The SOTABM is intended as an initial example of how ABMs can be used to dynamically represent mechanistic knowledge concerning transplant immunology in a scalable and expandable form and can thus potentially serve as useful adjuncts to the investigation and development of control strategies to induce transplant tolerance.

Strategies for efficient numerical implementation of hybrid multi-scale agent-based models to describe biological systems

Biologically related processes operate across multiple spatiotemporal scales. For computational modeling methodologies to mimic this biological complexity, individual scale models must be linked in ways that allow for dynamic exchange of information across scales. A powerful methodology is to combine a discrete modeling approach, agent-based models (ABMs), with continuum models to form hybrid models. Hybrid multi-scale ABMs have been used to simulate emergent responses of biological systems. Here, we review two aspects of hybrid multi-scale ABMs: linking individual scale models and efficiently solving the resulting model. We discuss the computational choices associated with aspects of linking individual scale models while simultaneously maintaining model tractability. We demonstrate implementations of existing numerical methods in the context of hybrid multi-scale ABMs. Using an example model describing Mycobacterium tuberculosis infection, we show relative computational speeds of various combinations of numerical methods. Efficient linking and solution of hybrid multi-scale ABMs is key to model portability, modularity, and their use in understanding biological phenomena at a systems level.

Integrated PK-PD and agent-based modeling in oncology

Mathematical modeling has become a valuable tool that strives to complement conventional biomedical research modalities in order to predict experimental outcome, generate new medical hypotheses, and optimize clinical therapies. Two specific approaches, pharmacokinetic-pharmacodynamic (PK-PD) modeling, and agent-based modeling (ABM), have been widely applied in cancer research. While they have made important contributions on their own (e.g., PK-PD in examining chemotherapy drug efficacy and resistance, and ABM in describing and predicting tumor growth and metastasis), only a few groups have started to combine both approaches together in an effort to gain more insights into the details of drug dynamics and the resulting impact on tumor growth. In this review, we focus our discussion on some of the most recent modeling studies building on a combined PK-PD and ABM approach that have generated experimentally testable hypotheses. Some future directions are also discussed.

Semantically linking in silico cancer models

Multiscale models are commonplace in cancer modeling, where individual models acting on different biological scales are combined within a single, cohesive modeling framework. However, model composition gives rise to challenges in understanding interfaces and interactions between them. Based on specific domain expertise, typically these computational models are developed by separate research groups using different methodologies, programming languages, and parameters. This paper introduces a graph-based model for semantically linking computational cancer models via domain graphs that can help us better understand and explore combinations of models spanning multiple biological scales. We take the data model encoded by TumorML, an XML-based markup language for storing cancer models in online repositories, and transpose its model description elements into a graph-based representation. By taking such an approach, we can link domain models, such as controlled vocabularies, taxonomic schemes, and ontologies, with cancer model descriptions to better understand and explore relationships between models. The union of these graphs creates a connected property graph that links cancer models by categorizations, by computational compatibility, and by semantic interoperability, yielding a framework in which opportunities for exploration and discovery of combinations of models become possible.

Systems biology approaches for studying the pathogenesis of non-alcoholic fatty liver disease

Non-alcoholic fatty liver disease (NAFLD) is a progressive disease of increasing public health concern. In western populations the disease has an estimated prevalence of 20%-40%, rising to 70%-90% in obese and type II diabetic individuals. Simplistically, NAFLD is the macroscopic accumulation of lipid in the liver, and is viewed as the hepatic manifestation of the metabolic syndrome. However, the molecular mechanisms mediating both the initial development of steatosis and its progression through non-alcoholic steatohepatitis to debilitating and potentially fatal fibrosis and cirrhosis are only partially understood. Despite increased research in this field, the development of non-invasive clinical diagnostic tools and the discovery of novel therapeutic targets has been frustratingly slow. We note that, to date, NAFLD research has been dominated by in vivo experiments in animal models and human clinical studies. Systems biology tools and novel computational simulation techniques allow the study of large-scale metabolic networks and the impact of their dysregulation on health. Here we review current systems biology tools and discuss the benefits to their application to the study of NAFLD. We propose that a systems approach utilising novel in silico modelling and simulation techniques is key to a more comprehensive, better targeted NAFLD research strategy. Such an approach will accelerate the progress of research and vital translation into clinic.

Agent-based model of fecal microbial transplant effect on bile acid metabolism on suppressing Clostridium difficile infection: an example of agent-based modeling of intestinal bacterial infection

Agent-based modeling is a computational modeling method that represents system-level behavior as arising from multiple interactions between the multiple components that make up a system. Biological systems are thus readily described using agent-based models (ABMs), as multi-cellular organisms can be viewed as populations of interacting cells, and microbial systems manifest as colonies of individual microbes. Intersections between these two domains underlie an increasing number of pathophysiological processes, and the intestinal tract represents one of the most significant locations for these inter-domain interactions, so much so that it can be considered an internal ecology of varying robustness and function. Intestinal infections represent significant disturbances of this internal ecology, and one of the most clinically relevant intestinal infections is Clostridium difficile infection (CDI). CDI is precipitated by the use of broad-spectrum antibiotics, involves the depletion of commensal microbiota, and alterations in bile acid composition in the intestinal lumen. We present an example ABM of CDI (the C. difficile Infection ABM, or CDIABM) to examine fundamental dynamics of the pathogenesis of CDI and its response to treatment with anti-CDI antibiotics and a newer treatment therapy, fecal microbial transplant. The CDIABM focuses on one specific mechanism of potential CDI suppression: commensal modulation of bile acid composition. Even given its abstraction, the CDIABM reproduces essential dynamics of CDI and its response to therapy, and identifies a paradoxical zone of behavior that provides insight into the role of intestinal nutritional status and the efficacy of anti-CDI therapies. It is hoped that this use case example of the CDIABM can demonstrate the usefulness of both agent-based modeling and the application of abstract functional representation as the biomedical community seeks to address the challenges of increasingly complex diseases with the goal of personalized medicine.

Simulating cancer growth with multiscale agent-based modeling

There have been many techniques developed in recent years to in silico model a variety of cancer behaviors. Agent-based modeling is a specific discrete-based hybrid modeling approach that allows simulating the role of diversity in cell populations as well as within each individual cell; it has therefore become a powerful modeling method widely used by computational cancer researchers. Many aspects of tumor morphology including phenotype-changing mutations, the adaptation to microenvironment, the process of angiogenesis, the influence of extracellular matrix, reactions to chemotherapy or surgical intervention, the effects of oxygen and nutrient availability, and metastasis and invasion of healthy tissues have been incorporated and investigated in agent-based models. In this review, we introduce some of the most recent agent-based models that have provided insight into the understanding of cancer growth and invasion, spanning multiple biological scales in time and space, and we further describe several experimentally testable hypotheses generated by those models. We also discuss some of the current challenges of multiscale agent-based cancer models.

Genome scale modeling in systems biology: algorithms and resources

In recent years, in silico studies and trial simulations have complemented experimental procedures. A model is a description of a system, and a system is any collection of interrelated objects; an object, moreover, is some elemental unit upon which observations can be made but whose internal structure either does not exist or is ignored. Therefore, any network analysis approach is critical for successful quantitative modeling of biological systems. This review highlights some of most popular and important modeling algorithms, tools, and emerging standards for representing, simulating and analyzing cellular networks in five sections. Also, we try to show these concepts by means of simple example and proper images and graphs. Overall, systems biology aims for a holistic description and understanding of biological processes by an integration of analytical experimental approaches along with synthetic computational models. In fact, biological networks have been developed as a platform for integrating information from high to low-throughput experiments for the analysis of biological systems. We provide an overview of all processes used in modeling and simulating biological networks in such a way that they can become easily understandable for researchers with both biological and mathematical backgrounds. Consequently, given the complexity of generated experimental data and cellular networks, it is no surprise that researchers have turned to computer simulation and the development of more theory-based approaches to augment and assist in the development of a fully quantitative understanding of cellular dynamics.

Multiple types of data are required to identify the mechanisms influencing the spatial expansion of melanoma cell colonies

BACKGROUND

The expansion of cell colonies is driven by a delicate balance of several mechanisms including cell motility, cell-to-cell adhesion and cell proliferation. New approaches that can be used to independently identify and quantify the role of each mechanism will help us understand how each mechanism contributes to the expansion process. Standard mathematical modelling approaches to describe such cell colony expansion typically neglect cell-to-cell adhesion, despite the fact that cell-to-cell adhesion is thought to play an important role.

RESULTS

We use a combined experimental and mathematical modelling approach to determine the cell diffusivity, D, cell-to-cell adhesion strength, q, and cell proliferation rate, λ, in an expanding colony of MM127 melanoma cells. Using a circular barrier assay, we extract several types of experimental data and use a mathematical model to independently estimate D, q and λ. In our first set of experiments, we suppress cell proliferation and analyse three different types of data to estimate D and q. We find that standard types of data, such as the area enclosed by the leading edge of the expanding colony and more detailed cell density profiles throughout the expanding colony, does not provide sufficient information to uniquely identify D and q. We find that additional data relating to the degree of cell-to-cell clustering is required to provide independent estimates of q, and in turn D. In our second set of experiments, where proliferation is not suppressed, we use data describing temporal changes in cell density to determine the cell proliferation rate. In summary, we find that our experiments are best described using the range D=161-243μm2 hour-1, q=0.3-0.5 (low to moderate strength) and λ=0.0305-0.0398 hour-1, and with these parameters we can accurately predict the temporal variations in the spatial extent and cell density profile throughout the expanding melanoma cell colony.

CONCLUSIONS

Our systematic approach to identify the cell diffusivity, cell-to-cell adhesion strength and cell proliferation rate highlights the importance of integrating multiple types of data to accurately quantify the factors influencing the spatial expansion of melanoma cell colonies.

Autism and EMF? Plausibility of a pathophysiological link part II

Autism spectrum conditions (ASCs) are defined behaviorally, but they also involve multileveled disturbances of underlying biology that find striking parallels in the physiological impacts of electromagnetic frequency and radiofrequency radiation exposures (EMF/RFR). Part I (Vol 776) of this paper reviewed the critical contributions pathophysiology may make to the etiology, pathogenesis and ongoing generation of behaviors currently defined as being core features of ASCs. We reviewed pathophysiological damage to core cellular processes that are associated both with ASCs and with biological effects of EMF/RFR exposures that contribute to chronically disrupted homeostasis. Many studies of people with ASCs have identified oxidative stress and evidence of free radical damage, cellular stress proteins, and deficiencies of antioxidants such as glutathione. Elevated intracellular calcium in ASCs may be due to genetics or may be downstream of inflammation or environmental exposures. Cell membrane lipids may be peroxidized, mitochondria may be dysfunctional, and various kinds of immune system disturbances are common. Brain oxidative stress and inflammation as well as measures consistent with blood-brain barrier and brain perfusion compromise have been documented. Part II of this paper documents how behaviors in ASCs may emerge from alterations of electrophysiological oscillatory synchronization, how EMF/RFR could contribute to these by de-tuning the organism, and policy implications of these vulnerabilities. It details evidence for mitochondrial dysfunction, immune system dysregulation, neuroinflammation and brain blood flow alterations, altered electrophysiology, disruption of electromagnetic signaling, synchrony, and sensory processing, de-tuning of the brain and organism, with autistic behaviors as emergent properties emanating from this pathophysiology. Changes in brain and autonomic nervous system electrophysiological function and sensory processing predominate, seizures are common, and sleep disruption is close to universal. All of these phenomena also occur with EMF/RFR exposure that can add to system overload ('allostatic load') in ASCs by increasing risk, and can worsen challenging biological problems and symptoms; conversely, reducing exposure might ameliorate symptoms of ASCs by reducing obstruction of physiological repair. Various vital but vulnerable mechanisms such as calcium channels may be disrupted by environmental agents, various genes associated with autism or the interaction of both. With dramatic increases in reported ASCs that are coincident in time with the deployment of wireless technologies, we need aggressive investigation of potential ASC-EMF/RFR links. The evidence is sufficient to warrant new public exposure standards benchmarked to low-intensity (non-thermal) exposure levels now known to be biologically disruptive, and strong, interim precautionary practices are advocated.

Predictive modelling of complex agronomic and biological systems

Biological systems are tremendously complex in their functioning and regulation. Studying the multifaceted behaviour and describing the performance of such complexity has challenged the scientific community for years. The reduction of real-world intricacy into simple descriptive models has therefore convinced many researchers of the usefulness of introducing mathematics into biological sciences. Predictive modelling takes such an approach another step further in that it takes advantage of existing knowledge to project the performance of a system in alternating scenarios. The ever growing amounts of available data generated by assessing biological systems at increasingly higher detail provide unique opportunities for future modelling and experiment design. Here we aim to provide an overview of the progress made in modelling over time and the currently prevalent approaches for iterative modelling cycles in modern biology. We will further argue for the importance of versatility in modelling approaches, including parameter estimation, model reduction and network reconstruction. Finally, we will discuss the difficulties in overcoming the mathematical interpretation of in vivo complexity and address some of the future challenges lying ahead.

Simulating the decentralized processes of the human immune system in a virtual anatomy model

Background

Many physiological processes within the human body can be perceived and modeled as large systems of interacting particles or swarming agents. The complex processes of the human immune system prove to be challenging to capture and illustrate without proper reference to the spacial distribution of immune-related organs and systems. Our work focuses on physical aspects of immune system processes, which we implement through swarms of agents. This is our first prototype for integrating different immune processes into one comprehensive virtual physiology simulation.

Results

Using agent-based methodology and a 3-dimensional modeling and visualization environment (LINDSAY Composer), we present an agent-based simulation of the decentralized processes in the human immune system. The agents in our model - such as immune cells, viruses and cytokines - interact through simulated physics in two different, compartmentalized and decentralized 3-dimensional environments namely, (1) within the tissue and (2) inside a lymph node. While the two environments are separated and perform their computations asynchronously, an abstract form of communication is allowed in order to replicate the exchange, transportation and interaction of immune system agents between these sites. The distribution of simulated processes, that can communicate across multiple, local CPUs or through a network of machines, provides a starting point to build decentralized systems that replicate larger-scale processes within the human body, thus creating integrated simulations with other physiological systems, such as the circulatory, endocrine, or nervous system. Ultimately, this system integration across scales is our goal for the LINDSAY Virtual Human project.

Conclusions

Our current immune system simulations extend our previous work on agent-based simulations by introducing advanced visualizations within the context of a virtual human anatomy model. We also demonstrate how to distribute a collection of connected simulations over a network of computers. As a future endeavour, we plan to use parameter tuning techniques on our model to further enhance its biological credibility. We consider these in silico experiments and their associated modeling and optimization techniques as essential components in further enhancing our capabilities of simulating a whole-body, decentralized immune system, to be used both for medical education and research as well as for virtual studies in immunoinformatics.

Longitudinal measures of proteostasis in live neurons: features that determine fate in models of neurodegenerative disease

Protein misfolding and proteostasis decline is a common feature of many neurodegenerative diseases. However, modeling the complexity of proteostasis and the global cellular consequences of its disruption is a challenge, particularly in live neurons. Although conventional approaches, based on population measures and single "snapshots", can identify cellular changes during neurodegeneration, they fail to determine if these cellular events drive cell death or act as adaptive responses. Alternatively, a "systems" cell biology approach known as longitudinal survival analysis enables single neurons to be followed over the course of neurodegeneration. By capturing the dynamics of misfolded proteins and the multiple cellular events that occur along the way, the relationship of these events to each other and their importance and role during cell death can be determined. Quantitative models of proteostasis dysfunction may yield unique insight and novel therapeutic strategies for neurodegenerative disease.

Agent-based models of cellular systems

Software agents are particularly suitable for engineering models and simulations of cellular systems. In a very natural and intuitive manner, individual software components are therein delegated to reproduce "in silico" the behavior of individual components of alive systems at a given level of resolution. Individuals' actions and interactions among individuals allow complex collective behavior to emerge. In this chapter we first introduce the readers to software agents and multi-agent systems, reviewing the evolution of agent-based modeling of biomolecular systems in the last decade. We then describe the main tools, platforms, and methodologies available for programming societies of agents, possibly profiting also of toolkits that do not require advanced programming skills.

The modulatory effect of cell–cell contact on the tumourigenic potential of pre-malignant epithelial cells: a computational exploration

Malignant development cannot be attributed alone to genetic changes in a single cell, but occurs as a result of the complex interplay between the failure of cellular regulation mechanisms and the presence of a permissive microenvironment. Although E-cadherin is classified as a 'metastasis suppressor' owing to its role in intercellular adhesion, the observation that it may be downregulated at a premalignant stage is indicative of additional roles in neoplastic development. We have used an agent-based computational model to explore the emergent behaviour resulting from the interaction of single and subpopulations of E-cadherin-compromised cells with unaffected normal epithelial cells within a monolayer environment. We have extended this to investigate the importance of local tissue perturbations in the form of scratch-wounding, or ablation of randomly-dispersed normal cells, on the growth of a single cell exhibiting E-cadherin loss. Our results suggest that the microenvironment with respect to localized cell density and normal/E-cadherin-compromised neighbours is crucial in determining whether an abnormal individual cell proliferates or remains dormant within the monolayer. These predictions raise important questions relating to the propensity for individual mutations to give rise to disease, and future experimental exploration of these will enhance our understanding of a complex, multifactorial pathological process.

A structural systems biology approach for quantifying the systemic consequences of missense mutations in proteins

Gauging the systemic effects of non-synonymous single nucleotide polymorphisms (nsSNPs) is an important topic in the pursuit of personalized medicine. However, it is a non-trivial task to understand how a change at the protein structure level eventually affects a cell's behavior. This is because complex information at both the protein and pathway level has to be integrated. Given that the idea of integrating both protein and pathway dynamics to estimate the systemic impact of missense mutations in proteins remains predominantly unexplored, we investigate the practicality of such an approach by formulating mathematical models and comparing them with experimental data to study missense mutations. We present two case studies: (1) interpreting systemic perturbation for mutations within the cell cycle control mechanisms (G2 to mitosis transition) for yeast; (2) phenotypic classification of neuron-related human diseases associated with mutations within the mitogen-activated protein kinase (MAPK) pathway. We show that the application of simplified mathematical models is feasible for understanding the effects of small sequence changes on cellular behavior. Furthermore, we show that the systemic impact of missense mutations can be effectively quantified as a combination of protein stability change and pathway perturbation.

In silico modelling of tumour margin diffusion and infiltration: review of current status

As a result of advanced treatment techniques, requiring precise target definitions, a need for more accurate delineation of the Clinical Target Volume (CTV) has arisen. Mathematical modelling is found to be a powerful tool to provide fairly accurate predictions for the Microscopic Extension (ME) of a tumour to be incorporated in a CTV. In general terms, biomathematical models based on a sequence of observations or development of a hypothesis assume some links between biological mechanisms involved in cancer development and progression to provide quantitative or qualitative measures of tumour behaviour as well as tumour response to treatment. Generally, two approaches are taken: deterministic and stochastic modelling. In this paper, recent mathematical models, including deterministic and stochastic methods, are reviewed and critically compared. It is concluded that stochastic models are more promising to provide a realistic description of cancer tumour behaviour due to being intrinsically probabilistic as well as discrete, which enables incorporation of patient-specific biomedical data such as tumour heterogeneity and anatomical boundaries.

Spatial simulations in systems biology: from molecules to cells

Cells are highly organized objects containing millions of molecules. Each biomolecule has a specific shape in order to interact with others in the complex machinery. Spatial dynamics emerge in this system on length and time scales which can not yet be modeled with full atomic detail. This review gives an overview of methods which can be used to simulate the complete cell at least with molecular detail, especially Brownian dynamics simulations. Such simulations require correct implementation of the diffusion-controlled reaction scheme occurring on this level. Implementations and applications of spatial simulations are presented, and finally it is discussed how the atomic level can be included for instance in multi-scale simulation methods.

Computational modelling of maternal interactions with spermatozoa: potentials and prospects

Understanding the complex interactions between gametes, embryos and the maternal tract is required knowledge for combating infertility and developing new methods of contraception. Here we present some main aspects of spermatozoa interactions with the mammalian oviduct before fertilisation and discuss how computational modelling can be used as an invaluable aid to experimental investigation in this field. A complete predictive computational model of gamete and embryo interactions with the female reproductive tract is a long way off. However, the enormity of this task should not discourage us from working towards it. Computational modelling allows us to investigate aspects of maternal communication with gametes and embryos, which are financially, ethically or practically difficult to look at experimentally. In silico models of maternal communication with gametes and embryos can be used as tools to complement in vivo experiments, in the same way as in vitro and in situ models.

Systems biology of fungal infection

Elucidation of pathogenicity mechanisms of the most important human-pathogenic fungi, Aspergillus fumigatus and Candida albicans, has gained great interest in the light of the steadily increasing number of cases of invasive fungal infections. A key feature of these infections is the interaction of the different fungal morphotypes with epithelial and immune effector cells in the human host. Because of the high level of complexity, it is necessary to describe and understand invasive fungal infection by taking a systems biological approach, i.e., by a comprehensive quantitative analysis of the non-linear and selective interactions of a large number of functionally diverse, and frequently multifunctional, sets of elements, e.g., genes, proteins, metabolites, which produce coherent and emergent behaviors in time and space. The recent advances in systems biology will now make it possible to uncover the structure and dynamics of molecular and cellular cause-effect relationships within these pathogenic interactions. We review current efforts to integrate omics and image-based data of host-pathogen interactions into network and spatio-temporal models. The modeling will help to elucidate pathogenicity mechanisms and to identify diagnostic biomarkers and potential drug targets for therapy and could thus pave the way for novel intervention strategies based on novel antifungal drugs and cell therapy.

Multiscale cancer modeling

Simulating cancer behavior across multiple biological scales in space and time, i.e., multiscale cancer modeling, is increasingly being recognized as a powerful tool to refine hypotheses, focus experiments, and enable more accurate predictions. A growing number of examples illustrate the value of this approach in providing quantitative insights in the initiation, progression, and treatment of cancer. In this review, we introduce the most recent and important multiscale cancer modeling works that have successfully established a mechanistic link between different biological scales. Biophysical, biochemical, and biomechanical factors are considered in these models. We also discuss innovative, cutting-edge modeling methods that are moving predictive multiscale cancer modeling toward clinical application. Furthermore, because the development of multiscale cancer models requires a new level of collaboration among scientists from a variety of fields such as biology, medicine, physics, mathematics, engineering, and computer science, an innovative Web-based infrastructure is needed to support this growing community.

Modeling formalisms in Systems Biology

Systems Biology has taken advantage of computational tools and high-throughput experimental data to model several biological processes. These include signaling, gene regulatory, and metabolic networks. However, most of these models are specific to each kind of network. Their interconnection demands a whole-cell modeling framework for a complete understanding of cellular systems. We describe the features required by an integrated framework for modeling, analyzing and simulating biological processes, and review several modeling formalisms that have been used in Systems Biology including Boolean networks, Bayesian networks, Petri nets, process algebras, constraint-based models, differential equations, rule-based models, interacting state machines, cellular automata, and agent-based models. We compare the features provided by different formalisms, and discuss recent approaches in the integration of these formalisms, as well as possible directions for the future.

From systems biology to systems biomedicine

Systems Biology is about combining theory, technology, and targeted experiments in a way that drives not only data accumulation but knowledge as well. The challenge in Systems Biomedicine is to furthermore translate mechanistic insights in biological systems to clinical application, with the central aim of improving patients' quality of life. The challenge is to find theoretically well-chosen models for the contextually correct and intelligible representation of multi-scale biological systems. In this review, we discuss the current state of Systems Biology, highlight the emergence of Systems Biomedicine, and highlight some of the topics and views that we think are important for the efficient application of Systems Theory in Biomedicine.

Systems biology of kidney diseases

Kidney diseases manifest in progressive loss of renal function, which ultimately leads to complete kidney failure. The mechanisms underlying the origins and progression of kidney diseases are not fully understood. Multiple factors involved in the pathogenesis of kidney diseases have made the traditional candidate gene approach of limited value toward full understanding of the molecular mechanisms of these diseases. A systems biology approach that integrates computational modeling with large-scale data gathering of the molecular changes could be useful in identifying the multiple interacting genes and their products that drive kidney diseases. Advances in biotechnology now make it possible to gather large data sets to characterize the role of the genome, epigenome, transcriptome, proteome, and metabolome in kidney diseases. When combined with computational analyses, these experimental approaches will provide a comprehensive understanding of the underlying biological processes. Multiscale analysis that connects the molecular interactions and cell biology of different kidney cells to renal physiology and pathology can be utilized to identify modules of biological and clinical importance that are perturbed in disease processes. This integration of experimental approaches and computational modeling is expected to generate new knowledge that can help to identify marker sets to guide the diagnosis, monitor disease progression, and identify new therapeutic targets.

Plant systems biology: network matters

Systems biology is all about networks. A recent trend has been to associate systems biology exclusively with the study of gene regulatory or protein-interaction networks. However, systems biology approaches can be applied at many other scales, from the subatomic to the ecosystem scales. In this review, we describe studies at the sub-cellular, tissue, whole plant and crop scales and highlight how these studies can be related to systems biology. We discuss the properties of system approaches at each scale as well as their current limits, and pinpoint in each case advances unique to the considered scale but representing potential for the other scales. We conclude by examining plant models bridging different scales and considering the future prospects of plant systems biology.

Discovering Molecular Targets in Cancer with Multiscale Modeling

Multiscale modeling is increasingly being recognized as a promising research area in computational cancer systems biology. Here, exemplified by two pioneering studies, we attempt to explain why and how such a multiscale approach paired with an innovative cross-scale analytical technique can be useful in identifying high-value molecular therapeutic targets. This novel, integrated approach has the potential to offer a more effective in silico framework for target discovery and represents an important technical step towards systems medicine.