Full Title: Spatio-temporal BMP signaling model
This simulation reproduces Fig. 3B of the referenced paper. Parameter values are as published in Suppl. Section "Mathematical modeling" except for the size of the mBMP expressing clone, see Note 1 below. The results illustrate that different signal transduction kinetics can determine specific activity ranges of pSmad5 and pSmad2 in response to the same gradient of mBMP4, the latter being established by a localized source and diffusion, decay in a 2D domain.
A radial coordinate "x" is used to optionally plot the concentration profiles. Also optionally, time courses of mean concentrations over the domain can be plotted. Both plotting options are included but temporarily inactivated in the model file.
Reference: The model is described in the peer-reviewed publication "Integration of Nodal and BMP Signaling by Mutual Signaling Effector Antagonism" by Gary Huiming Soh, Autumn Penecilla Pomreinke, Patrick MÃ¼ller. Cell Reports 31, 107487, 2020, https://doi.org/10.1016/j.celrep.2020.03.051
Author of Morpheus model: peter.brusch.dd@gmail.com
Date of Morpheus model: November 21, 2020
Software: Morpheus (open-source), download from https://morpheus.gitlab.io
Units are: [space]=micrometer, [time]=second
Note: 1. This Morpheus model uses r=50 (micrometers) for the radius of the (blue) source region in the center of the (black) simulation domain with radius 300 micrometers. This corresponds to the geometry used in Fig. 3B of the referenced paper.
2. The resulting state at 7200sec is a transient state as the mBMP4 field is initialized as zero, a steady state with higher concentration values is approached at ten times longer duration.
3. The color scales for the two displayed fields are chosen such that red corresponds to maximum pSmad5 of 5000 and green to maximum pSmad2 of 5.