Heterodimer-heterotetramer formation mediates enhanced sensor activity in a biophysical model for BMP signaling

Numerous stages of organismal development rely on the cellular interpretation of gradients of secreted morphogens including members of the Bone Morphogenetic Protein (BMP) family through transmembrane receptors. Early gradients of BMPs drive dorsal/ventral patterning throughout the animal kingdom in both vertebrates and invertebrates. Growing evidence in Drosophila, zebrafish, murine and other systems suggests that BMP ligand heterodimers are the primary BMP signaling ligand, even in systems in which mixtures of BMP homodimers and heterodimers are present. Signaling by heterodimers occurs through a hetero-tetrameric receptor complex comprising of two distinct type one BMP receptors and two type II receptors. To understand the system dynamics and determine whether kinetic assembly of heterodimer-heterotetramer BMP complexes is favored, as compared to other plausible BMP ligand-receptor configurations, we developed a kinetic model for BMP tetramer formation based on current measurements for binding rates and affinities. We find that contrary to a common hypothesis, heterodimer-heterotetramer formation is not kinetically favored over the formation of homodimer-tetramer complexes under physiological conditions of receptor and ligand concentrations and therefore other mechanisms, potentially including differential kinase activities of the formed heterotetramer complexes, must be the cause of heterodimer-heterotetramer signaling primacy. Further, although BMP complex assembly favors homodimer and homomeric complex formation over a wide range of parameters, ignoring these signals and instead relying on the heterodimer improves the range of morphogen interpretation in a broad set of conditions, suggesting a performance advantage for heterodimer signaling in patterning multiple cell types in a gradient. TGF-beta signaling is an important cell signaling system through which cells respond to external information. In the TGF-beta system, signaling is initiated when a ligand dimer pair binds to a receptor tetramer. Interestingly, in numerous developmental contexts, TGF-beta signaling has a greater response to heterodimeric ligands (dimers of multiple ligands), as compared to homomeric ligands (dimers made of two molecules of a single ligand. However, neither the cause of heterodimer signaling primacy, nor the systemic effects of heterodimeric vs homomeric signaling are understood. We use a biophysically-informed computational modeling approach to investigate the system dynamics of heterodimer-heterotetramer BMP signaling, to understand the cause and consequence of the requirement for BMP2/7-mediated signaling in dorsoventral patterning in zebrafish development. Using our model, we demonstrate that BMP heterodimer signaling complex formation is not kinetically favored over homodimer signaling complexes, suggesting subfunctionalization of BMP receptors may be required to explain heterodimer signaling. Additionally, we show that heterodimer signaling provides a performance advantage via increased range of morphogen interpretation. Our findings provide insight into the systems principles involved in developmental signaling.