Theoretical Biology and Biophysics Group, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico, United States of America.
Systems Biology Ireland, University College Dublin, Belfield, Dublin, Ireland.
PLoS Comput Biol. 2019 Jan 17;15(1):e1006706. doi: 10.1371/journal.pcbi.1006706. eCollection 2019 Jan.
Receptor tyrosine kinases (RTKs) typically contain multiple autophosphorylation sites in their cytoplasmic domains. Once activated, these autophosphorylation sites can recruit downstream signaling proteins containing Src homology 2 (SH2) and phosphotyrosine-binding (PTB) domains, which recognize phosphotyrosine-containing short linear motifs (SLiMs). These domains and SLiMs have polyspecific or promiscuous binding activities. Thus, multiple signaling proteins may compete for binding to a common SLiM and vice versa. To investigate the effects of competition on RTK signaling, we used a rule-based modeling approach to develop and analyze models for ligand-induced recruitment of SH2/PTB domain-containing proteins to autophosphorylation sites in the insulin-like growth factor 1 (IGF1) receptor (IGF1R). Models were parameterized using published datasets reporting protein copy numbers and site-specific binding affinities. Simulations were facilitated by a novel application of model restructuration, to reduce redundancy in rule-derived equations. We compare predictions obtained via numerical simulation of the model to those obtained through simple prediction methods, such as through an analytical approximation, or ranking by copy number and/or KD value, and find that the simple methods are unable to recapitulate the predictions of numerical simulations. We created 45 cell line-specific models that demonstrate how early events in IGF1R signaling depend on the protein abundance profile of a cell. Simulations, facilitated by model restructuration, identified pairs of IGF1R binding partners that are recruited in anti-correlated and correlated fashions, despite no inclusion of cooperativity in our models. This work shows that the outcome of competition depends on the physicochemical parameters that characterize pairwise interactions, as well as network properties, including network connectivity and the relative abundances of competitors.
受体酪氨酸激酶(RTKs)通常在其细胞质结构域中含有多个自身磷酸化位点。一旦被激活,这些自身磷酸化位点可以募集含有Src 同源 2(SH2)和磷酸酪氨酸结合(PTB)结构域的下游信号蛋白,这些结构域识别含有磷酸酪氨酸的短线性基序(SLiM)。这些结构域和 SLiM 具有多特异性或混杂的结合活性。因此,多个信号蛋白可能会竞争结合共同的 SLiM,反之亦然。为了研究竞争对 RTK 信号的影响,我们使用基于规则的建模方法,开发和分析胰岛素样生长因子 1(IGF1)受体(IGF1R)中配体诱导的 SH2/PTB 结构域蛋白募集到自身磷酸化位点的模型。使用报告蛋白质拷贝数和特定位点结合亲和力的已发表数据集对模型进行参数化。通过模型重构的新应用简化了模拟,减少了规则衍生方程中的冗余。我们将通过模型数值模拟获得的预测与通过简单预测方法获得的预测进行比较,例如通过分析近似值或通过拷贝数和/或 KD 值进行排序,结果发现简单方法无法再现数值模拟的预测。我们创建了 45 个特定于细胞系的模型,展示了 IGF1R 信号的早期事件如何取决于细胞的蛋白质丰度谱。通过模型重构进行的模拟确定了以反相关和相关方式募集的 IGF1R 结合伴侣对,尽管我们的模型中没有包含协同作用。这项工作表明,竞争的结果取决于表征成对相互作用的物理化学参数,以及网络特性,包括网络连通性和竞争者的相对丰度。