ML-Space: hybrid spatial Gillespie and Brownian motion simulation at multiple levels, and a rule-based description language

Computer simulations of biological cells as well-stirred systems are well established but neglect the spatial distribution of key actors. In this thesis, a simulation algorithm "ML-Space" for spatial models with dynamic hierarchies is presented. It combines stochastic spatial algorithms in discretized space with individual particles moving in continuous space that have spatial extensions and can contain other particles. For formal descriptions of the systems to be simulated spatially, ML-Space provides a rule-based specification language.Computersimulationen mikrobiologischer Prozesse, bei denen eine homogene Verteilung der Akteure einer Zelle angenommen wird, sind gut etabliert. In dieser Arbeit wird ein räumlicher Simulationsalgorithmus "ML-Space" für Mehrebenenmodelle vorgestellt, der auch dynamische Hierarchien abdeckt. Er vereint stochastische räumliche Algorithmen in diskretisiertem Raum mit individuellen Partikeln mit kontinuierlichen Koordinaten, die andere Partikel enthalten können. Zur formalen Beschreibung der räumlich zu simulierenden Systeme bietet ML-Space eine regelbasierte Modellierungssprache

Quantification of growth asymmetries in developing epithelia

Many developmental processes of multicellular organisms involve the patterning and growth of two-dimensional tissues, so called epithelia. We have quantified the growth of the wing imaginal disk, which is the precursor of the adult wing, of the fruit fly Drosophila melanogaster. We find that growth follows a simple rule with exponentially decreasing area growth rate. Anisotropies of growth can be precisely determined by comparing experimental results to a continuum theory. Growth anisotropies are to good approximation constant in space and time. They are weak in wild-type wing disks but threefold increased in GFP-Dpp disks in which the morphogen Dpp is overexpressed. Our findings indicate that morphogens such as Dpp control tissue shape via oriented cell divisions that generate anisotropic growt

Incorporating chemical signalling factors into cell-based models of growing epithelial tissues

In this paper we present a comprehensive computational framework within which the effects of chemical signalling factors on growing epithelial tissues can be studied. The method incorporates a vertex-based cell model, in conjunction with a solver for the governing chemical equations. The vertex model provides a natural mesh for the finite element method (FEM), with node movements determined by force laws. The arbitrary Lagrangian–Eulerian formulation is adopted to account for domain movement between iterations. The effects of cell proliferation and junctional rearrangements on the mesh are also examined. By implementing refinements of the mesh we show that the finite element (FE) approximation converges towards an accurate numerical solution. The potential utility of the system is demonstrated in the context of Decapentaplegic (Dpp), a morphogen which plays a crucial role in development of the Drosophila imaginal wing disc. Despite the presence of a Dpp gradient, growth is uniform across the wing disc. We make the growth rate of cells dependent on Dpp concentration and show that the number of proliferation events increases in regions of high concentration. This allows hypotheses regarding mechanisms of growth control to be rigorously tested. The method we describe may be adapted to a range of potential application areas, and to other cell-based models with designated node movements, to accurately probe the role of morphogens in epithelial tissues

Bottom mixed layer oxygen dynamics in the Celtic Sea

The seasonally stratified continental shelf seas are highly productive, economically important environments which are under considerable pressure from human activity. Global dissolved oxygen concentrations have shown rapid reductions in response to anthropogenic forcing since at least the middle of the twentieth century. Oxygen consumption is at the same time linked to the cycling of atmospheric carbon, with oxygen being a proxy for carbon remineralisation and the release of CO2. In the seasonally stratified seas the bottom mixed layer (BML) is partially isolated from the atmosphere and is thus controlled by interplay between oxygen consumption processes, vertical and horizontal advection. Oxygen consumption rates can be both spatially and temporally dynamic, but these dynamics are often missed with incubation based techniques. Here we adopt a Bayesian approach to determining total BML oxygen consumption rates from a high resolution oxygen time-series. This incorporates both our knowledge and our uncertainty of the various processes which control the oxygen inventory. Total BML rates integrate both processes in the water column and at the sediment interface. These observations span the stratified period of the Celtic Sea and across both sandy and muddy sediment types. We show how horizontal advection, tidal forcing and vertical mixing together control the bottom mixed layer oxygen concentrations at various times over the stratified period. Our muddy-sand site shows cyclic spring-neap mediated changes in oxygen consumption driven by the frequent resuspension or ventilation of the seabed. We see evidence for prolonged periods of increased vertical mixing which provide the ventilation necessary to support the high rates of consumption observed

Colloquium: Mechanical formalisms for tissue dynamics

The understanding of morphogenesis in living organisms has been renewed by tremendous progressin experimental techniques that provide access to cell-scale, quantitative information both on theshapes of cells within tissues and on the genes being expressed. This information suggests that ourunderstanding of the respective contributions of gene expression and mechanics, and of their crucialentanglement, will soon leap forward. Biomechanics increasingly benefits from models, which assistthe design and interpretation of experiments, point out the main ingredients and assumptions, andultimately lead to predictions. The newly accessible local information thus calls for a reflectionon how to select suitable classes of mechanical models. We review both mechanical ingredientssuggested by the current knowledge of tissue behaviour, and modelling methods that can helpgenerate a rheological diagram or a constitutive equation. We distinguish cell scale ("intra-cell")and tissue scale ("inter-cell") contributions. We recall the mathematical framework developpedfor continuum materials and explain how to transform a constitutive equation into a set of partialdifferential equations amenable to numerical resolution. We show that when plastic behaviour isrelevant, the dissipation function formalism appears appropriate to generate constitutive equations;its variational nature facilitates numerical implementation, and we discuss adaptations needed in thecase of large deformations. The present article gathers theoretical methods that can readily enhancethe significance of the data to be extracted from recent or future high throughput biomechanicalexperiments.Comment: 33 pages, 20 figures. This version (26 Sept. 2015) contains a few corrections to the published version, all in Appendix D.2 devoted to large deformation

Regge Field Theory in zero transverse dimensions: loops versus "net" diagrams

Toy models of interacting Pomerons with triple and quaternary Pomeron vertices in zero transverse dimension are investigated. Numerical solutions for eigenvalues and eigenfunctions of the corresponding Hamiltonians are obtained, providing the quantum solution for the scattering amplitude in both models. The equations of motion for the Lagrangians of the theories are also considered and the classical solutions of the equations are found. Full two-point Green functions ("effective" Pomeron propagator) and amplitude of diffractive dissociation process are calculated in the framework of RFT-0 approach. The importance of the loops contribution in the amplitude at different values of the model parameters is discussed as well as the difference between the models with and without quaternary Pomeron vertex.Comment: 34 pages, 36 figure

Dynamics and Mechanical Stability of the Developing Dorsoventral Organizer of the Wing Imaginal Disc

Shaping the primordia during development relies on forces and mechanisms able to control cell segregation. In the imaginal discs of Drosophila the cellular populations that will give rise to the dorsal and ventral parts on the wing blade are segregated and do not intermingle. A cellular population that becomes specified by the boundary of the dorsal and ventral cellular domains, the so-called organizer, controls this process. In this paper we study the dynamics and stability of the dorsal-ventral organizer of the wing imaginal disc of Drosophila as cell proliferation advances. Our approach is based on a vertex model to perform in silico experiments that are fully dynamical and take into account the available experimental data such as: cell packing properties, orientation of the cellular divisions, response upon membrane ablation, and robustness to mechanical perturbations induced by fast growing clones. Our results shed light on the complex interplay between the cytoskeleton mechanics, the cell cycle, the cell growth, and the cellular interactions in order to shape the dorsal-ventral organizer as a robust source of positional information and a lineage controller. Specifically, we elucidate the necessary and sufficient ingredients that enforce its functionality: distinctive mechanical properties, including increased tension, longer cell cycle duration, and a cleavage criterion that satisfies the Hertwig rule. Our results provide novel insights into the developmental mechanisms that drive the dynamics of the DV organizer and set a definition of the so-called Notch fence model in quantitative terms