Simulating non-small cell lung cancer with a multiscale agent-based model

Background The epidermal growth factor receptor (EGFR) is frequently overexpressed in many cancers, including non-small cell lung cancer (NSCLC). In silcio modeling is considered to be an increasingly promising tool to add useful insights into the dynamics of the EGFR signal transduction pathway. However, most of the previous modeling work focused on the molecular or the cellular level only, neglecting the crucial feedback between these scales as well as the interaction with the heterogeneous biochemical microenvironment. Results We developed a multiscale model for investigating expansion dynamics of NSCLC within a two-dimensional in silico microenvironment. At the molecular level, a specific EGFR-ERK intracellular signal transduction pathway was implemented. Dynamical alterations of these molecules were used to trigger phenotypic changes at the cellular level. Examining the relationship between extrinsic ligand concentrations, intrinsic molecular profiles and microscopic patterns, the results confirmed that increasing the amount of available growth factor leads to a spatially more aggressive cancer system. Moreover, for the cell closest to nutrient abundance, a phase-transition emerges where a minimal increase in extrinsic ligand abolishes the proliferative phenotype altogether. Conclusions Our in silico results indicate that, in NSCLC, in the presence of a strong extrinsic chemotactic stimulus, and depending on the cell's location, downstream EGFR-ERK signaling may be processed more efficiently, thereby yielding a migration-dominant cell phenotype and overall, an accelerated spatio-temporal expansion rate.Comment: 37 pages, 7 figure

Simulating Brain Tumor Heterogeneity with a Multiscale Agent-Based Model: Linking Molecular Signatures, Phenotypes and Expansion Rate

We have extended our previously developed 3D multi-scale agent-based brain tumor model to simulate cancer heterogeneity and to analyze its impact across the scales of interest. While our algorithm continues to employ an epidermal growth factor receptor (EGFR) gene-protein interaction network to determine the cells' phenotype, it now adds an explicit treatment of tumor cell adhesion related to the model's biochemical microenvironment. We simulate a simplified tumor progression pathway that leads to the emergence of five distinct glioma cell clones with different EGFR density and cell 'search precisions'. The in silico results show that microscopic tumor heterogeneity can impact the tumor system's multicellular growth patterns. Our findings further confirm that EGFR density results in the more aggressive clonal populations switching earlier from proliferation-dominated to a more migratory phenotype. Moreover, analyzing the dynamic molecular profile that triggers the phenotypic switch between proliferation and migration, our in silico oncogenomics data display spatial and temporal diversity in documenting the regional impact of tumorigenesis, and thus support the added value of multi-site and repeated assessments in vitro and in vivo. Potential implications from this in silico work for experimental and computational studies are discussed.Comment: 37 pages, 10 figure

Identification of Critical Molecular Components in a Multiscale Cancer Model Based on the Integration of Monte Carlo, Resampling, and ANOVA

To date, parameters defining biological properties in multiscale disease models are commonly obtained from a variety of sources. It is thus important to examine the influence of parameter perturbations on system behavior, rather than to limit the model to a specific set of parameters. Such sensitivity analysis can be used to investigate how changes in input parameters affect model outputs. However, multiscale cancer models require special attention because they generally take longer to run than does a series of signaling pathway analysis tasks. In this article, we propose a global sensitivity analysis method based on the integration of Monte Carlo, resampling, and analysis of variance. This method provides solutions to (1) how to render the large number of parameter variation combinations computationally manageable, and (2) how to effectively quantify the sampling distribution of the sensitivity index to address the inherent computational intensity issue. We exemplify the feasibility of this method using a two-dimensional molecular-microscopic agent-based model previously developed for simulating non-small cell lung cancer; in this model, an epidermal growth factor (EGF)-induced, EGF receptor-mediated signaling pathway was implemented at the molecular level. Here, the cross-scale effects of molecular parameters on two tumor growth evaluation measures, i.e., tumor volume and expansion rate, at the microscopic level are assessed. Analysis finds that ERK, a downstream molecule of the EGF receptor signaling pathway, has the most important impact on regulating both measures. The potential to apply this method to therapeutic target discovery is discussed

Diversification, Relatedness, And Firm Performance: Empirical Evidence From China

The relationship between diversification, relatedness and performance has long been a controversial issue in mainstream strategic management research. Research in this area, however, has focused primarily on developed countries. This study argues that the conclusions drawn from developed countries may not apply to developing countries. In an investigation of 227 publicly-listed companies in China, this study found that: 1) firm scale significantly contributes to the improvement of economic performance; 2) relatedness correlates negatively with firm performance, and 3) the relationship between diversification and performance fits the intermediate model. This study also provided evidence to support the argument that differences do exist in the rationales between firms in developed and developing countries

Dynamic Targeting in Cancer Treatment

With the advent of personalized medicine, design and development of anti-cancer drugs that are specifically targeted to individual or sets of genes or proteins has been an active research area in both academia and industry. The underlying motivation for this approach is to interfere with several pathological crosstalk pathways in order to inhibit or at the very least control the proliferation of cancer cells. However, after initially conferring beneficial effects, if sub-lethal, these artificial perturbations in cell function pathways can inadvertently activate drug-induced up- and down-regulation of feedback loops, resulting in dynamic changes over time in the molecular network structure and potentially causing drug resistance as seen in clinics. Hence, the targets or their combined signatures should also change in accordance with the evolution of the network (reflected by changes to the structure and/or functional output of the network) over the course of treatment. This suggests the need for a “dynamic targeting” strategy aimed at optimizing tumor control by interfering with different molecular targets, at varying stages. Understanding the dynamic changes of this complex network under various perturbed conditions due to drug treatment is extremely challenging under experimental conditions let alone in clinical settings. However, mathematical modeling can facilitate studying these effects at the network level and beyond, and also accelerate comparison of the impact of different dosage regimens and therapeutic modalities prior to sizeable investment in risky and expensive clinical trials. A dynamic targeting strategy based on the use of mathematical modeling can be a new, exciting research avenue in the discovery and development of therapeutic drugs

Differential neuroproteomic and systems biology analysis of spinal cord injury

Acute spinal cord injury (SCI) is a devastating condition with many consequences and no known effective treatment. Although it is quite easy to diagnose traumatic SCI, the assessment of injury severity and projection of disease progression or recovery are often challenging, as no consensus biomarkers have been clearly identified. Here rats were subjected to experimental moderate or severe thoracic SCI. At 24h and 7d postinjury, spinal cord segment caudal to injury center versus sham samples was harvested and subjected to differential proteomic analysis. Cationic/anionic-exchange chromatography, followed by 1D polyacrylamide gel electrophoresis, was used to reduce protein complexity. A reverse phase liquid chromatography-tandem mass spectrometry proteomic platform was then utilized to identify proteome changes associated with SCI. Twenty-two and 22 proteins were up-regulated at 24 h and 7 day after SCI, respectively; whereas 19 and 16 proteins are down-regulated at 24 h and 7 day after SCI, respectively, when compared with sham control. A subset of 12 proteins were identified as candidate SCI biomarkers - TF (Transferrin), FASN (Fatty acid synthase), NME1 (Nucleoside diphosphate kinase 1), STMN1 (Stathmin 1), EEF2 (Eukaryotic translation elongation factor 2), CTSD (Cathepsin D), ANXA1 (Annexin A1), ANXA2 (Annexin A2), PGM1 (Phosphoglucomutase 1), PEA15 (Phosphoprotein enriched in astrocytes 15), GOT2 (Glutamic-oxaloacetic transaminase 2), and TPI-1 (Triosephosphate isomerase 1), data are available via ProteomeXchange with identifier PXD003473. In addition, Transferrin, Cathepsin D, and TPI-1 and PEA15 were further verified in rat spinal cord tissue and/or CSF samples after SCI and in human CSF samples from moderate/severe SCI patients. Lastly, a systems biology approach was utilized to determine the critical biochemical pathways and interactome in the pathogenesis of SCI. Thus, SCI candidate biomarkers identified can be used to correlate with disease progression or to identify potential SCI therapeutic targets

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

Radiation fluxes in a business district of Shanghai, China

Radiative fluxes are key drivers of surface-atmosphere heat exchanges in cities. Here the first year-long (December 2012 – November 2013) measurements of the full radiation balance for a dense urban site in Shanghai are presented, collected with a net radiometer CNR4 mounted 80 m above ground. Clear sky incoming shortwave radiation (K↓) (median daytime maxima) ranges from 575 W m-2 in winter to 875 W m-2 in spring, with cloud cover reducing the daily maxima by about 160 W m-2. The median incoming longwave radiation daytime maxima is 305 and 468 W m-2 in winter and summer, respectively, with increases of 30 and 15 W m-2 for cloudy conditions. The effect of air quality is evident: ‘haze’ conditions decrease hourly median K↓ by 11.3%. The midday (11:00 -13:00 LST) clear sky surface albedo (α) is 0.128, 0.141, 0.143 and 0.129 for winter, spring, summer and autumn, respectively. α varies with solar elevation and azimuth angle due to heterogeneity of the urban surface. In winter, shadows play an important role in decreasing α in the late afternoon. For the site, the bulk α is 0.14. The NARP/SUEWS land surface model reproduces the radiation components at this site well, a promising result for applications elsewhere. These observations help to fill the gap of long-term radiation measurements in East Asian and low-latitude cities quantifying the effects of season, cloud cover and air quality