MOTIVATION: Biological reality can in silico only be comprehensively represented in multi-scaled models. To this end, cell behavioural models addressing the multi-cellular level have to be semantically linked with mechanistic molecular models. These requirements have to be met by flexible software workflows solving the issues of different time scales, inter-model variable referencing and flexible sub-model embedding. RESULTS: We developed a novel software workflow (EPISIM) for the semantic integration of Systems Biology Markup Language (SBML)-based quantitative models in multi-scaled tissue models and simulations. This workflow allows to import and access SBML-based models. SBML model species, reactions and parameters are semantically integrated in cell behavioural models (CBM) represented by graphical process diagrams. By this, cellular states like proliferation and differentiation can be flexibly linked to gene-regulatory or biochemical reaction networks. For a multi-scale agent-based tissue simulation executable code is automatically generated where different time scales of imported SBML models and CBM have been mapped. We demonstrate the capabilities of the novel software workflow by integrating Tyson's cell cycle model in our model of human epidermal tissue homeostasis. Finally, we show the semantic interplay of the different biological scales during tissue simulation. AVAILABILITY: The EPISIM platform is available as binary executables for Windows, Linux and Mac OS X at http://www.tiga.uni-hd.de. Supplementary data are available at http://www.tiga.uni-hd.de/supplements/SemSBMLIntegration.html. CONTACT: [email protected].

MOTIVATION: For a mechanistic understanding of skin and its response to an induced perturbation, systems biology is gaining increasing attention. Unfortunately, quantitative and spatial expression data for skin, like for most other tissues, are almost not available. RESULTS: Integrating organotypic skin cultures, whole-slide scanning and subsequent image processing provides bioinformatics with a novel source of spatial expression data. We here used this approach to quantitatively describe the effect of treating organotypic skin cultures with sodium dodecyl sulphate in a non-corrosive concentration. We first measured the differentiation-related spatial expression gradient of Heat-Shock-Protein 27 in a time series of up to 24 h. Secondly, a multi-dimensional tissue classifier for predicting skin irritation was developed based on abstract features of these profiles. We obtained a high specificity of 0.94 and a sensitivity of 0.92 compared with manual classification. Our results demonstrate that the integration of tissue cultures, whole-slide scanning and image processing is well suited for both the standardized data acquisition for systems biological tissue models and a highly robust classification of tissue responses.

MOTIVATION: For the efficient application of multi-agent systems to spatial and functional modeling of tissues flexible and intuitive modeling tools are needed, which allow the graphical specification of cellular behavior in a tissue context without presuming specialized programming skills. RESULTS: We developed a graphical modeling system for multi-agent based simulation of tissue homeostasis. An editor allows the intuitive and hierarchically structured specification of cellular behavior. The models are then automatically compiled into highly efficient source code and dynamically linked to an interactive graphical simulation environment. The system allows the quantitative analysis of the morphological and functional tissue properties emerging from the cell behavioral model. We demonstrate the relevance of the approach using a recently published model of epidermal homeostasis as well as a series of cell-cycle models. AVAILABILITY: The complete software is available in binary executables for MS-Windows and Linux at tiga.uni-hd.de.

BACKGROUND: The epidermal physiology results from a complex regulated homeostasis of keratinocyte proliferation, differentiation and death and is tightly regulated by a specific protein expression during cellular maturation. Cellular in silico models are considered a promising and inevitable tool for the understanding of this complex system. Hence, we need to incorporate the information of the differentiation dependent protein expression in cell based systems biological models of tissue homeostasis. Such methods require measuring tissue differentiation quantitatively while correlating it with biomarker expression intensities. RESULTS: Differentiation of a keratinocyte is characterized by its continuously changing morphology concomitant with its movement from the basal layer to the surface, leading to a decreased average nuclei density throughout the tissue. Based thereon, we designed and evaluated three different mathematical measures (nuclei based, distance based, and joint approach) for quantifying differentiation in epidermal keratinocytes. We integrated them with an immunofluorescent staining and image analysis method for tissue sections, automatically quantifying epidermal differentiation and measuring the corresponding expression of biomarkers. When studying five well-known differentiation related biomarkers in an epidermal neck sample only the resulting biomarker profiles incorporating the relative distance information of cells to the tissue borders (distance based and joint approach) provided a high-resolution view on the whole process of keratinocyte differentiation. By contrast, the inverse nuclei density approach led to an increased resolution at early but heavily decreased resolution at late differentiation. This effect results from the heavy non-linear decay of DAPI intensity per area, probably caused by cytoplasmic growth and chromatin decondensation. In the joint approach this effect could be compensated again by incorporating distance information. CONCLUSION: We suppose that key mechanisms regulating tissue homeostasis probably depend more on distance information rather than on nuclei reorganization. Concluding, the distance approach appears well suited for comprehensively observing keratinocyte differentiation.

MOTIVATION: For systems biology of complex stratified epithelia like human epidermis, it will be of particular importance to reconstruct the spatiotemporal gene and protein networks regulating keratinocyte differentiation and homeostasis. RESULTS: Inside the epidermis, the differentiation state of individual keratinocytes is correlated with their respective distance from the connective tissue. We here present a novel method to profile this correlation for multiple epithelial protein biomarkers in the form of quantitative spatial profiles. Profiles were computed by applying image processing algorithms to histological sections stained with tri-color indirect immunofluorescence. From the quantitative spatial profiles, reflecting the spatiotemporal changes of protein expression during cellular differentiation, graphs of protein networks were reconstructed. CONCLUSION: Spatiotemporal networks can be used as a means for comparing and interpreting quantitative spatial protein expression profiles obtained from different tissue samples. In combination with automated microscopes, our new method supports the large-scale systems biological analysis of stratified epithelial tissues.

An investigation of the chaotic behavior in heart rate variability (HRV) signals of a normal person and a heart transplant recipient is presented. A statistical approach of chaos identification in time series is described and applied to HRV signals. This method compares the short-term predictability for a given time series to an ensemble of random data which has the same Fourier spectrum as the original time series. The forecasting error is computed as a statistic for performing statistical hypothesis testing. The results suggest that the HRV signal of the transplant recipient recorded three months after the transplanting shows the same dynamical behavior as that of the HRV signal for a normal person.