To model dynamic expression patterns in somitogenesis we developed a Java-application for simulating gene regulatory networks in many cells in parallel and visualising the results using the Java3D API, thus simulating the collective behaviour of many thousand cells. According to the 'clock-and-wave-front' model mesodermal segmentation of vertebrate embryos is regulated by a 'segmentation clock', which oscillates with a period of about 2h in mice, and a 'wave front' moving back with the growing caudal end of the presomitic mesoderm. The clock is realised through cycling expression of genes such as Hes1 and Hes7, whose gene products repress the transcription of their encoding genes in a negative feedback loop. By coupling the decay of the Hes1 mRNA to a gradient with the same features and mechanism of formation as the mesodermal Fgf8 gradient we can simulate typical features of the dynamic expression pattern of Hes1 in the presomitic mesoderm. Furthermore, our program is able to synchronise Hes1 oscillations in thousands of cells through simulated Delta-Notch signalling interactions.

The authors propose piecewise deterministic Markov processes as an alternative approach to model gene regulatory networks. A hybrid simulation algorithm is presented and discussed, and several standard regulatory modules are analysed by numerical means. It is shown that despite of the partial simplification of the mesoscopic nature of regulatory networks such processes are suitable to reveal the intrinsic noise effects because of the low copy numbers of genes.

Somitogenesis describes the segmentation of vertebrate embryonic bodies, which is thought to be induced by ultradian clocks (i.e., clocks with relatively short cycles compared to circadian clocks). One candidate for such a clock is the bHLH factor Hes1, forming dimers which repress the transcription of its own encoding gene. Most models for such small autoregulative networks are based on delay equations where a Hill function represents the regulation of transcription. The aim of the present paper is to estimate the Hill coefficient in the switch of an Hes1 oscillator and to suggest a more detailed model of the autoregulative network. The promoter of Hes1 consists of three to four binding sites for Hes1 dimers. Using the sparse data from literature, we find, in contrast to other statements in literature, that there is not much evidence for synergistic binding in the regulatory region of Hes1, and that the Hill coefficient is about three. As a model for the negative feedback loop, we use a Goodwin system and find sustained oscillations for systems with a large enough number of linear differential equations. By a suitable variation of the number of equations, we provide a rational lower bound for the Hill coefficient for such a system. Our results suggest that there exist additional nonlinear processes outside of the regulatory region of Hes1.

The segmentation of the vertebrate body is laid down during early embryogenesis. The formation of signaling gradients, the periodic expression of genes of the Notch-, Fgf- and Wnt-pathways and their interplay in the unsegmented presomitic mesoderm (PSM) precedes the rhythmic budding of nascent somites at its anterior end, which later develops into epithelialized structures, the somites. Although many in silico models describing partial aspects of somitogenesis already exist, simulations of a complete causal chain from gene expression in the growth zone via the interaction of multiple cells to segmentation are rare. Here, we present an enhanced gene regulatory network (GRN) for mice in a simulation program that models the growing PSM by many virtual cells and integrates WNT3A and FGF8 gradient formation, periodic gene expression and Delta/Notch signaling. Assuming Hes7 as core of the somitogenesis clock and LFNG as modulator, we postulate a negative feedback of HES7 on Dll1 leading to an oscillating Dll1 expression as seen in vivo. Furthermore, we are able to simulate the experimentally observed wave of activated NOTCH (NICD) as a result of the interactions in the GRN. We esteem our model as robust for a wide range of parameter values with the Hes7 mRNA and protein decays exerting a strong influence on the core oscillator. Moreover, our model predicts interference between Hes1 and HES7 oscillators when their intrinsic frequencies differ. In conclusion, we have built a comprehensive model of somitogenesis with HES7 as core oscillator that is able to reproduce many experimentally observed data in mice.