Important signal transduction pathways originate on the plasma membrane, where microdomains may transiently entrap diffusing receptors. This results in a non-random distribution of receptors even in the resting state, which can be visualized as "clusters" by high resolution imaging methods. Here, we explore how spatial in-homogeneities in the plasma membrane might influence the dimerization and phosphorylation status of ErbB2 and ErbB3, two receptor tyrosine kinases that preferentially heterodimerize and are often co-expressed in cancer. This theoretical study is based upon spatial stochastic simulations of the two-dimensional membrane landscape, where variables include differential distributions and overlap of transient confinement zones ("domains") for the two receptor species. The in silico model is parameterized and validated using data from single particle tracking experiments. We report key differences in signaling output based on the degree of overlap between domains and the relative retention of receptors in such domains, expressed as escape probability. Results predict that a high overlap of domains, which favors transient co-confinement of both receptor species, will enhance the rate of hetero-interactions. Where domains do not overlap, simulations confirm expectations that homo-interactions are favored. Since ErbB3 is uniquely dependent on ErbB2 interactions for activation of its catalytic activity, variations in domain overlap or escape probability markedly alter the predicted patterns and time course of ErbB3 and ErbB2 phosphorylation. Taken together, these results implicate membrane domain organization as an important modulator of signal initiation, motivating the design of novel experimental approaches to measure these important parameters across a wider range of receptor systems.
Keywords: ErbB receptors; ErbB2; ErbB3; membrane domains; spatial stochastic modeling.