Adhesion between cells, diffusion of growth factors, and elasticity of the AER produce the paddle shape of the chick limb

This paper has been withdrawnComment: This paper has been withdraw

Computer Simulation of Cellular Patterning Within the Drosophila Pupal Eye

We present a computer simulation and associated experimental validation of assembly of glial-like support cells into the interweaving hexagonal lattice that spans the Drosophila pupal eye. This process of cell movements organizes the ommatidial array into a functional pattern. Unlike earlier simulations that focused on the arrangements of cells within individual ommatidia, here we examine the local movements that lead to large-scale organization of the emerging eye field. Simulations based on our experimental observations of cell adhesion, cell death, and cell movement successfully patterned a tracing of an emerging wild-type pupal eye. Surprisingly, altering cell adhesion had only a mild effect on patterning, contradicting our previous hypothesis that the patterning was primarily the result of preferential adhesion between IRM-class surface proteins. Instead, our simulations highlighted the importance of programmed cell death (PCD) as well as a previously unappreciated variable: the expansion of cells' apical surface areas, which promoted rearrangement of neighboring cells. We tested this prediction experimentally by preventing expansion in the apical area of individual cells: patterning was disrupted in a manner predicted by our simulations. Our work demonstrates the value of combining computer simulation with in vivo experiments to uncover novel mechanisms that are perpetuated throughout the eye field. It also demonstrates the utility of the Glazier–Graner–Hogeweg model (GGH) for modeling the links between local cellular interactions and emergent properties of developing epithelia as well as predicting unanticipated results in vivo

3D Multi-Cell Simulation of Tumor Growth and Angiogenesis

We present a 3D multi-cell simulation of a generic simplification of vascular tumor growth which can be easily extended and adapted to describe more specific vascular tumor types and host tissues. Initially, tumor cells proliferate as they take up the oxygen which the pre-existing vasculature supplies. The tumor grows exponentially. When the oxygen level drops below a threshold, the tumor cells become hypoxic and start secreting pro-angiogenic factors. At this stage, the tumor reaches a maximum diameter characteristic of an avascular tumor spheroid. The endothelial cells in the pre-existing vasculature respond to the pro-angiogenic factors both by chemotaxing towards higher concentrations of pro-angiogenic factors and by forming new blood vessels via angiogenesis. The tumor-induced vasculature increases the growth rate of the resulting vascularized solid tumor compared to an avascular tumor, allowing the tumor to grow beyond the spheroid in these linear-growth phases. First, in the linear-spherical phase of growth, the tumor remains spherical while its volume increases. Second, in the linear-cylindrical phase of growth the tumor elongates into a cylinder. Finally, in the linear-sheet phase of growth, tumor growth accelerates as the tumor changes from cylindrical to paddle-shaped. Substantial periods during which the tumor grows slowly or not at all separate the exponential from the linear-spherical and the linear-spherical from the linear-cylindrical growth phases. In contrast to other simulations in which avascular tumors remain spherical, our simulated avascular tumors form cylinders following the blood vessels, leading to a different distribution of hypoxic cells within the tumor. Our simulations cover time periods which are long enough to produce a range of biologically reasonable complex morphologies, allowing us to study how tumor-induced angiogenesis affects the growth rate, size and morphology of simulated tumors

SBML Level 3: an extensible format for the exchange and reuse of biological models

Abstract Systems biology has experienced dramatic growth in the number, size, and complexity of computational models. To reproduce simulation results and reuse models, researchers must exchange unambiguous model descriptions. We review the latest edition of the Systems Biology Markup Language (SBML), a format designed for this purpose. A community of modelers and software authors developed SBML Level 3 over the past decade. Its modular form consists of a core suited to representing reaction‐based models and packages that extend the core with features suited to other model types including constraint‐based models, reaction‐diffusion models, logical network models, and rule‐based models. The format leverages two decades of SBML and a rich software ecosystem that transformed how systems biologists build and interact with models. More recently, the rise of multiscale models of whole cells and organs, and new data sources such as single‐cell measurements and live imaging, has precipitated new ways of integrating data with models. We provide our perspectives on the challenges presented by these developments and how SBML Level 3 provides the foundation needed to support this evolution

Adhesion Failures Determine the Pattern of Choroidal Neovascularization in the Eye: A Computer Simulation Study

Choroidal neovascularization (CNV) of the macular area of the retina is the major cause of severe vision loss in adults. In CNV, after choriocapillaries initially penetrate Bruch's membrane (BrM), invading vessels may regress or expand (CNV initiation). Next, during Early and Late CNV, the expanding vasculature usually spreads in one of three distinct patterns: in a layer between BrM and the retinal pigment epithelium (sub-RPE or Type 1 CNV), in a layer between the RPE and the photoreceptors (sub-retinal or Type 2 CNV) or in both loci simultaneously (combined pattern or Type 3 CNV). While most studies hypothesize that CNV primarily results from growth-factor effects or holes in BrM, our three-dimensional simulations of multi-cell model of the normal and pathological maculae recapitulate the three growth patterns, under the hypothesis that CNV results from combinations of impairment of: 1) RPE-RPE epithelial junctional adhesion, 2) Adhesion of the RPE basement membrane complex to BrM (RPE-BrM adhesion), and 3) Adhesion of the RPE to the photoreceptor outer segments (RPE-POS adhesion). Our key findings are that when an endothelial tip cell penetrates BrM: 1) RPE with normal epithelial junctions, basal attachment to BrM and apical attachment to POS resists CNV. 2) Small holes in BrM do not, by themselves, initiate CNV. 3) RPE with normal epithelial junctions and normal apical RPE-POS adhesion, but weak adhesion to BrM (e.g. due to lipid accumulation in BrM) results in Early sub-RPE CNV. 4) Normal adhesion of RBaM to BrM, but reduced apical RPE-POS or epithelial RPE-RPE adhesion (e.g. due to inflammation) results in Early sub-retinal CNV. 5) Simultaneous reduction in RPE-RPE epithelial binding and RPE-BrM adhesion results in either sub-RPE or sub-retinal CNV which often progresses to combined pattern CNV. These findings suggest that defects in adhesion dominate CNV initiation and progression

SBML Level 3: an extensible format for the exchange and reuse of biological models

Systems biology has experienced dramatic growth in the number, size, and complexity of computational models. To reproduce simulation results and reuse models, researchers must exchange unambiguous model descriptions. We review the latest edition of the Systems Biology Markup Language (SBML), a format designed for this purpose. A community of modelers and software authors developed SBML Level 3 over the past decade. Its modular form consists of a core suited to representing reaction-based models and packages that extend the core with features suited to other model types including constraint-based models, reaction-diffusion models, logical network models, and rule-based models. The format leverages two decades of SBML and a rich software ecosystem that transformed how systems biologists build and interact with models. More recently, the rise of multiscale models of whole cells and organs, and new data sources such as single-cell measurements and live imaging, has precipitated new ways of integrating data with models. We provide our perspectives on the challenges presented by these developments and how SBML Level 3 provides the foundation needed to support this evolution

Bifurcation analysis of regulatory modules in cell biology

Das Kernstueck der vorliegenden Arbeit ist die Betonung von kleinen Modulen als Schluesselkomponenten von biologischen Netzwerken. Unter den zahlreichen moeglichen Modulen scheinen besondere diejenigen interessant zu sein, welche die Rueckkopplungen realisieren und in regulatorischen Einheiten auftreten. Prozesse wie Genregulation, Differentiation oder Homeostasis benoetigen haeufig Autoregulation. Auf Grund dessen ist die detaillierte Kenntnis der dynamischen Eigenschaften von kleinen Modulen von groesserem Interesse. Es werden zwei biologische Systeme analysiert. Das erste beschaftigt sich mit dem Zellzyklus, das zweite Beispiel kommt aus der Immunologie und betrifft die Aktivierung von T-Zellen. Beide Modelle, d.h. ihre zugrundeliegende Netzwerke, lassen sich in Untereinheiten mit wohldefinierten Funktionen zerlegen. Diese Module entscheiden ueber das Verhalten des gesamten Netzwerkes. Mit anderen Worten, die von den Modulen getroffenen Entscheidungen, werden von dem gesamten System uebernommen. Bei der Analyse des Modells zum Zellzyklus wurde eine interessante Eigenschaft von gekoppelten Modulen deutlich, die wir dann getrennt behandelt haben. Seriell geschaltete Module mit positiver Rueckkopplung liefern ueberraschende Konstruktionsmoeglichkeiten fuer Systeme mit mehreren stabilen Gleichgewichtslagen. Obwohl nicht alle hier aufgestellten Hypothesen derzeit experimentell ueberpruefbar sind, es kann eine wichtige Aussage getroffen werden. Uebereinstimmende Strukturen und Mechanismen, die in verschiedenen biologischen Systemen vorkommen, bieten uns die Moeglichkeit einer Klassifizierung von biologischen Systemen bezueglich ihrer strukturellen Aehnlichkeiten.The thesis emphasizes the importance of small modules as key components of biological networks. Especially, those which perform positive feedbacks seem to be involved in a number of regulatory units. Processes like gene regulation, differentiation and homeostasis often require autoregulation. Therefore, detailed knowledge of dynamics of small modules becomes nowadays an important subject of study. We analyze two biological systems: one regarding cell cycle regulation and one immunological example related to T-cell activation. Their underlying networks can be dissected into subunits with well defined functions. These modules decide about the behavior of the global network. In other words, they have decision taking function, which is inherited by the whole system. Stimulated by the cell cycle model and its interesting dynamics resulting from coupled modules, we analyzed the switching issue separately. Serial coupling of positive feedback circuits provides astonishing possibilities to construct systems with multiple stable steady states. Even though, in current stage, no exact experimental proof of all hypotheses is possible, one important observation can be made. Common structures and mechanisms found in different biological systems allow to classify biological systems with respect to their structural similarities

Modeling and Simulation Tools: From Systems Biology to Systems Medicine

Modeling is an integral component of modern biology. In this chapter we look into the role of the model, as it pertains to Systems Medicine, and the software that is required to instantiate and run it. We do this by comparing the development, implementation, and characteristics of tools that have been developed to work with two divergent methodologies: Systems Biology and Pharmacometrics. From the Systems Biology perspective we consider the concept of “Software as a Medical Device” and what this may imply for the migration of research-oriented, simulation software into the domain of human health

Systems oncology:towards patient-specific treatment regimes informed by multiscale mathematical modelling

The multiscale complexity of cancer as a disease necessitates a corresponding multiscale modelling approach to produce truly predictive mathematical models capable of improving existing treatment protocols. To capture all the dynamics of solid tumour growth and its progression, mathematical modellers need to couple biological processes occurring at various spatial and temporal scales (from genes to tissues). Because effectiveness of cancer therapy is considerably affected by intracellular and extracellular heterogeneities as well as by the dynamical changes in the tissue microenvironment, any model attempt to optimise existing protocols must consider these factors ultimately leading to improved multimodal treatment regimes. By improving existing and building new mathematical models of cancer, modellers can play important role in preventing the use of potentially sub-optimal treatment combinations. In this paper, we analyse a multiscale computational mathematical model for cancer growth and spread, incorporating the multiple effects of radiation therapy and chemotherapy in the patient survival probability and implement the model using two different cell based modelling techniques. We show that the insights provided by such multiscale modelling approaches can ultimately help in designing optimal patient-specific multi-modality treatment protocols that may increase patients quality of life.PostprintPeer reviewe