Multiscale modeling of dorsoventral patterning in Drosophila.

The role of mathematical models of signaling networks is showcased by examples from Drosophila development. Three models of consecutive stages in dorsoventral patterning are presented. We begin with a compartmental model of intracellular reactions that generates a gradient of nuclear-localized Dorsal, exhibiting constant shape and dynamic amplitude. A simple thermodynamic model of equilibrium binding explains how a spatially uniform transcription factor, Zelda, can act in combination with a graded factor, Dorsal, to cooperatively regulate gene expression borders. Finally, we formulate a dynamic and stochastic model that predicts spatiotemporal patterns of Sog expression based on known patterns of its transcription factor, Dorsal. The future of coupling multifarious models across multiple temporal and spatial scales is discussed.

Distinguishing graded and ultrasensitive signalling cascade kinetics by the shape of morphogen gradients in Drosophila.

Recently, signalling gradients in cascades of two-state reaction-diffusion systems were described as a model for understanding key biochemical mechanisms that underlie development and differentiation processes in the Drosophila embryo. Diffusion-trapping at the exterior of the cell membrane triggers the mitogen-activated protein kinase (MAPK) cascade to relay an appropriate signal from the membrane to the inner part of the cytosol, whereupon another diffusion-trapping mechanism involving the nucleus reads out this signal to trigger appropriate changes in gene expression. Proposed mathematical models exhibit equilibrium distributions consistent with experimental measurements of key spatial gradients in these processes. A significant property of the formulation is that the signal is assumed to be relayed from one system to the next in a linear fashion. However, the MAPK cascade often exhibits nonlinear dose-response properties and the final remark of Berezhkovskii et al. (2009) is that this assumption remains an important property to be tested experimentally, perhaps via a new quantitative assay across multiple genetic backgrounds. In anticipation of the need to be able to sensibly interpret data from such experiments, here we provide a complementary analysis that recovers existing formulae as a special case but is also capable of handling nonlinear functional forms. Predictions of linear and nonlinear signal relays and, in particular, graded and ultrasensitive MAPK kinetics, are compared.

Stochastic chemical kinetics and the total quasi-steady-state assumption: application to the stochastic simulation algorithm and chemical master equation.

Recently the application of the quasi-steady-state approximation (QSSA) to the stochastic simulation algorithm (SSA) was suggested for the purpose of speeding up stochastic simulations of chemical systems that involve both relatively fast and slow chemical reactions [Rao and Arkin, J. Chem. Phys. 118, 4999 (2003)] and further work has led to the nested and slow-scale SSA. Improved numerical efficiency is obtained by respecting the vastly different time scales characterizing the system and then by advancing only the slow reactions exactly, based on a suitable approximation to the fast reactions. We considerably extend these works by applying the QSSA to numerical methods for the direct solution of the chemical master equation (CME) and, in particular, to the finite state projection algorithm [Munsky and Khammash, J. Chem. Phys. 124, 044104 (2006)], in conjunction with Krylov methods. In addition, we point out some important connections to the literature on the (deterministic) total QSSA (tQSSA) and place the stochastic analogue of the QSSA within the more general framework of aggregation of Markov processes. We demonstrate the new methods on four examples: Michaelis-Menten enzyme kinetics, double phosphorylation, the Goldbeter-Koshland switch, and the mitogen activated protein kinase cascade. Overall, we report dramatic improvements by applying the tQSSA to the CME solver.