Orchestration of ErbB3 signaling through heterointeractions and homointeractions.

Members of the ErbB family of receptor tyrosine kinases are capable of both homointeractions and heterointeractions. Because each receptor has a unique set of binding sites for downstream signaling partners and differential catalytic activity, subtle shifts in their combinatorial interplay may have a large effect on signaling outcomes. The overexpression and mutation of ErbB family members are common in numerous human cancers and shift the balance of activation within the signaling network. Here we report the development of a spatial stochastic model that addresses the dynamics of ErbB3 homodimerization and heterodimerization with ErbB2. The model is based on experimental measures for diffusion, dimer off-rates, kinase activity, and dephosphorylation. We also report computational analysis of ErbB3 mutations, generating the prediction that activating mutations in the intracellular and extracellular domains may be subdivided into classes with distinct underlying mechanisms. We show experimental evidence for an ErbB3 gain-of-function point mutation located in the C-lobe asymmetric dimerization interface, which shows enhanced phosphorylation at low ligand dose associated with increased kinase activity.

Analytical solution of steady-state equations for chemical reaction networks with bilinear rate laws.

True steady states are a rare occurrence in living organisms, yet their knowledge is essential for quasi-steady-state approximations, multistability analysis, and other important tools in the investigation of chemical reaction networks (CRN) used to describe molecular processes on the cellular level. Here, we present an approach that can provide closed form steady-state solutions to complex systems, resulting from CRN with binary reactions and mass-action rate laws. We map the nonlinear algebraic problem of finding steady states onto a linear problem in a higher-dimensional space. We show that the linearized version of the steady-state equations obeys the linear conservation laws of the original CRN. We identify two classes of problems for which complete, minimally parameterized solutions may be obtained using only the machinery of linear systems and a judicious choice of the variables used as free parameters. We exemplify our method, providing explicit formulae, on CRN describing signal initiation of two important types of RTK receptor-ligand systems, VEGF and EGF-ErbB1.