An in silico hybrid continuum-/agent-based procedure to modelling cancer development: Interrogating the interplay amongst glioma invasion, vascularity and necrosis

This paper develops a three-dimensional in silico hybrid model of cancer, which describes the multi-variate phenotypic behaviour of tumour and host cells. The model encompasses the role of cell migration and adhesion, the influence of the extracellular matrix, the effects of oxygen and nutrient availability, and the signalling triggered by chemical cues and growth factors. The proposed in silico hybrid modelling framework combines successfully the advantages of continuum-based and discrete methods, namely the finite element and agent-based method respectively. The framework is thus used to realistically model cancer mechano-biology in a multiscale fashion while maintaining the resolution power of each method in a computationally cost-effective manner. The model is tailored to simulate glioma progression, and is subsequently used to interrogate the balance between the host cells and small sized gliomas, while the go-or-grow phenotype characteristic in glioblastomas is also investigated. Also, cell-cell and cell-matrix interactions are examined with respect to their effect in (macroscopic) tumour growth, brain tissue perfusion and tumour necrosis. Finally, we use the in silico framework to assess differences between low-grade and high-grade glioma growth, demonstrating significant differences in the distribution of cancer as well as host cells, in accordance with reported experimental findings.

Functionalized Anatomical Models for Computational Life Sciences

The advent of detailed computational anatomical models has opened new avenues for computational life sciences (CLS). To date, static models representing the anatomical environment have been used in many applications but are insufficient when the dynamics of the body prevents separation of anatomical geometrical variability from physics and physiology. Obvious examples include the assessment of thermal risks in magnetic resonance imaging and planning for radiofrequency and acoustic cancer treatment, where posture and physiology-related changes in shape (e.g., breathing) or tissue behavior (e.g., thermoregulation) affect the impact. Advanced functionalized anatomical models can overcome these limitations and dramatically broaden the applicability of CLS in basic research, the development of novel devices/therapies, and the assessment of their safety and efficacy. Various forms of functionalization are discussed in this paper: (i) shape parametrization (e.g., heartbeat, population variability), (ii) physical property distributions (e.g., image-based inhomogeneity), (iii) physiological dynamics (e.g., tissue and organ behavior), and (iv) integration of simulation/measurement data (e.g., exposure conditions, “validation evidence” supporting model tuning and validation). Although current model functionalization may only represent a small part of the physiology, it already facilitates the next level of realism by (i) driving consistency among anatomy and different functionalization layers and highlighting dependencies, (ii) enabling third-party use of validated functionalization layers as established simulation tools, and (iii) therefore facilitating their application as building blocks in network or multi-scale computational models. Integration in functionalized anatomical models thus leverages and potentiates the value of sub-models and simulation/measurement data toward ever-increasing simulation realism. In our o²S²PARC platform, we propose to expand the concept of functionalized anatomical models to establish an integration and sharing service for heterogeneous computational models, ranging from the molecular to the organ level. The objective of o²S²PARC is to integrate all models developed within the National Institutes of Health SPARC initiative in a unified anatomical and computational environment, to study the role of the peripheral nervous system in controlling organ physiology. The functionalization concept, as outlined for the o²S²PARC platform, could form the basis for many other application areas of CLS. The relationship to other ongoing initiatives, such as the Physiome Project, is also presented.

Statistical ensemble of gene regulatory networks of macrophage differentiation

BACKGROUND

Macrophages cover a major role in the immune system, being the most plastic cell yielding several key immune functions.

METHODS

Here we derived a minimalistic gene regulatory network model for the differentiation of macrophages into the two phenotypes M1 (pro-) and M2 (anti-inflammatory).

RESULTS

To test the model, we simulated a large number of such networks as in a statistical ensemble. In other words, to enable the inter-cellular crosstalk required to obtain an immune activation in which the macrophage plays its role, the simulated networks are not taken in isolation but combined with other cellular agents, thus setting up a discrete minimalistic model of the immune system at the microscopic/intracellular (i.e., genetic regulation) and mesoscopic/intercellular scale.

CONCLUSIONS

We show that within the mesoscopic level description of cellular interaction and cooperation, the gene regulatory logic is coherent and contributes to the overall dynamics of the ensembles that shows, statistically, the expected behaviour.

Multiscale modelling in immunology: a review

One of the greatest challenges in biomedicine is to get a unified view of observations made from the molecular up to the organism scale. Towards this goal, multiscale models have been highly instrumental in contexts such as the cardiovascular field, angiogenesis, neurosciences and tumour biology. More recently, such models are becoming an increasingly important resource to address immunological questions as well. Systematic mining of the literature in multiscale modelling led us to identify three main fields of immunological applications: host-virus interactions, inflammatory diseases and their treatment and development of multiscale simulation platforms for immunological research and for educational purposes. Here, we review the current developments in these directions, which illustrate that multiscale models can consistently integrate immunological data generated at several scales, and can be used to describe and optimize therapeutic treatments of complex immune diseases.

Multiscale digital Arabidopsis predicts individual organ and whole-organism growth

Understanding how dynamic molecular networks affect whole-organism physiology, analogous to mapping genotype to phenotype, remains a key challenge in biology. Quantitative models that represent processes at multiple scales and link understanding from several research domains can help to tackle this problem. Such integrated models are more common in crop science and ecophysiology than in the research communities that elucidate molecular networks. Several laboratories have modeled particular aspects of growth in Arabidopsis thaliana, but it was unclear whether these existing models could productively be combined. We test this approach by constructing a multiscale model of Arabidopsis rosette growth. Four existing models were integrated with minimal parameter modification (leaf water content and one flowering parameter used measured data). The resulting framework model links genetic regulation and biochemical dynamics to events at the organ and whole-plant levels, helping to understand the combined effects of endogenous and environmental regulators on Arabidopsis growth. The framework model was validated and tested with metabolic, physiological, and biomass data from two laboratories, for five photoperiods, three accessions, and a transgenic line, highlighting the plasticity of plant growth strategies. The model was extended to include stochastic development. Model simulations gave insight into the developmental control of leaf production and provided a quantitative explanation for the pleiotropic developmental phenotype caused by overexpression of miR156, which was an open question. Modular, multiscale models, assembling knowledge from systems biology to ecophysiology, will help to understand and to engineer plant behavior from the genome to the field.

Modeling biology spanning different scales: an open challenge

It is coming nowadays more clear that in order to obtain a unified description of the different mechanisms governing the behavior and causality relations among the various parts of a living system, the development of comprehensive computational and mathematical models at different space and time scales is required. This is one of the most formidable challenges of modern biology characterized by the availability of huge amount of high throughput measurements. In this paper we draw attention to the importance of multiscale modeling in the framework of studies of biological systems in general and of the immune system in particular.

A vision and strategy for the virtual physiological human: 2012 update

European funding under Framework 7 (FP7) for the virtual physiological human (VPH) project has been in place now for 5 years. The VPH Network of Excellence (NoE) has been set up to help develop common standards, open source software, freely accessible data and model repositories, and various training and dissemination activities for the project. It is also working to coordinate the many clinically targeted projects that have been funded under the FP7 calls. An initial vision for the VPH was defined by the FP6 STEP project in 2006. In 2010, we wrote an assessment of the accomplishments of the first two years of the VPH in which we considered the biomedical science, healthcare and information and communications technology challenges facing the project (Hunter et al. 2010 Phil. Trans. R. Soc. A 368, 2595-2614 (doi:10.1098/rsta.2010.0048)). We proposed that a not-for-profit professional umbrella organization, the VPH Institute, should be established as a means of sustaining the VPH vision beyond the time-frame of the NoE. Here, we update and extend this assessment and in particular address the following issues raised in response to Hunter et al.: (i) a vision for the VPH updated in the light of progress made so far, (ii) biomedical science and healthcare challenges that the VPH initiative can address while also providing innovation opportunities for the European industry, and (iii) external changes needed in regulatory policy and business models to realize the full potential that the VPH has to offer to industry, clinics and society generally.

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.

3D Imaging of vascular networks for biophysical modeling of perfusion distribution within the heart

One of the main determinants of perfusion distribution within an organ is the structure of its vascular network. Past studies were based on angiography or corrosion casting and lacked quantitative three dimensional, 3D, representation. Based on branching rules and other properties derived from such imaging, 3D vascular tree models were generated which were rather useful for generating and testing hypotheses on perfusion distribution in organs. Progress in advanced computational models for prediction of perfusion distribution has raised the need for more realistic representations of vascular trees with higher resolution. This paper presents an overview of the different methods developed over time for imaging and modeling the structure of vascular networks and perfusion distribution, with a focus on the heart. The strengths and limitations of these different techniques are discussed. Episcopic fluorescent imaging using a cryomicrotome is presently being developed in different laboratories. This technique is discussed in more detail, since it provides high-resolution 3D structural information that is important for the development and validation of biophysical models but also for studying the adaptations of vascular networks to diseases. An added advantage of this method being is the ability to measure local tissue perfusion. Clinically, indices for patient-specific coronary stenosis evaluation derived from vascular networks have been proposed and high-resolution noninvasive methods for perfusion distribution are in development. All these techniques depend on a proper representation of the relevant vascular network structures.

Multiscale integration of -omic, imaging, and clinical data in biomedical informatics

This paper reviews challenges and opportunities in multiscale data integration for biomedical informatics. Biomedical data can come from different biological origins, data acquisition technologies, and clinical applications. Integrating such data across multiple scales (e.g., molecular, cellular/tissue, and patient) can lead to more informed decisions for personalized, predictive, and preventive medicine. However, data heterogeneity, community standards in data acquisition, and computational complexity are big challenges for such decision making. This review describes genomic and proteomic (i.e., molecular), histopathological imaging (i.e., cellular/tissue), and clinical (i.e., patient) data; it includes case studies for single-scale (e.g., combining genomic or histopathological image data), multiscale (e.g., combining histopathological image and clinical data), and multiscale and multiplatform (e.g., the Human Protein Atlas and The Cancer Genome Atlas) data integration. Numerous opportunities exist in biomedical informatics research focusing on integration of multiscale and multiplatform data.

Minimum Information about a Cardiac Electrophysiology Experiment (MICEE): standardised reporting for model reproducibility, interoperability, and data sharing

Cardiac experimental electrophysiology is in need of a well-defined Minimum Information Standard for recording, annotating, and reporting experimental data. As a step towards establishing this, we present a draft standard, called Minimum Information about a Cardiac Electrophysiology Experiment (MICEE). The ultimate goal is to develop a useful tool for cardiac electrophysiologists which facilitates and improves dissemination of the minimum information necessary for reproduction of cardiac electrophysiology research, allowing for easier comparison and utilisation of findings by others. It is hoped that this will enhance the integration of individual results into experimental, computational, and conceptual models. In its present form, this draft is intended for assessment and development by the research community. We invite the reader to join this effort, and, if deemed productive, implement the Minimum Information about a Cardiac Electrophysiology Experiment standard in their own work.

Theoretical models for coronary vascular biomechanics: progress & challenges

A key aim of the cardiac Physiome Project is to develop theoretical models to simulate the functional behaviour of the heart under physiological and pathophysiological conditions. Heart function is critically dependent on the delivery of an adequate blood supply to the myocardium via the coronary vasculature. Key to this critical function of the coronary vasculature is system dynamics that emerge via the interactions of the numerous constituent components at a range of spatial and temporal scales. Here, we focus on several components for which theoretical approaches can be applied, including vascular structure and mechanics, blood flow and mass transport, flow regulation, angiogenesis and vascular remodelling, and vascular cellular mechanics. For each component, we summarise the current state of the art in model development, and discuss areas requiring further research. We highlight the major challenges associated with integrating the component models to develop a computational tool that can ultimately be used to simulate the responses of the coronary vascular system to changing demands and to diseases and therapies.