Computational modeling reveals molecular details of epidermal growth factor binding

An integrated mechanistic and data-driven computational model predicts cell responses to high- and low-affinity EGFR ligands

The biophysical properties of ligand binding heavily influence the ability of receptors to specify cell fates. Understanding the rules by which ligand binding kinetics impact cell phenotype is challenging, however, because of the coupled information transfers that occur from receptors to downstream signaling effectors and from effectors to phenotypes. Here, we address that issue by developing an integrated mechanistic and data-driven computational modeling platform to predict cell responses to different ligands for the epidermal growth factor receptor (EGFR). Experimental data for model training and validation were generated using MCF7 human breast cancer cells treated with the high- and low-affinity ligands epidermal growth factor (EGF) and epiregulin (EREG), respectively. The integrated model captures the unintuitive, concentration-dependent abilities of EGF and EREG to drive signals and phenotypes differently, even at similar levels of receptor occupancy. For example, the model correctly predicts the dominance of EREG over EGF in driving a cell differentiation phenotype through AKT signaling at intermediate and saturating ligand concentrations and the ability of EGF and EREG to drive a broadly concentration-sensitive migration phenotype through cooperative ERK and AKT signaling. Parameter sensitivity analysis identifies EGFR endocytosis, which is differentially regulated by EGF and EREG, as one of the most important determinants of the alternative phenotypes driven by different ligands. The integrated model provides a new platform to predict how phenotypes are controlled by the earliest biophysical rate processes in signal transduction and may eventually be leveraged to understand receptor signaling system performance depends on cell context.

One-sentence summary

Integrated kinetic and data-driven EGFR signaling model identifies the specific signaling mechanisms that dictate cell responses to EGFR activation by different ligands.

Applications of personalised signalling network models in precision oncology

As our ability to provide in-depth, patient-specific characterisation of the molecular alterations within tumours rapidly improves, it is becoming apparent that new approaches will be required to leverage the power of this data and derive the full benefit for each individual patient. Systems biology approaches are beginning to emerge within this field as a potential method of incorporating large volumes of network level data and distilling a coherent, clinically-relevant prediction of drug response. However, the initial promise of this developing field is yet to be realised. Here we argue that in order to develop these precise models of individual drug response and tailor treatment accordingly, we will need to develop mathematical models capable of capturing both the dynamic nature of drug-response signalling networks and key patient-specific information such as mutation status or expression profiles. We also review the modelling approaches commonly utilised within this field, and outline recent examples of their use in furthering the application of systems biology for a precision medicine approach to cancer treatment.

Cooperativity transitions driven by higher-order oligomer formations in ligand-induced receptor dimerization

While cooperativity in ligand-induced receptor dimerization has been linked with receptor-receptor couplings via minimal representations of physical observables, effects arising from higher-order oligomer, e.g., trimer and tetramer, formations of unobserved receptors have received less attention. Here we propose a dimerization model of ligand-induced receptors in multivalent form representing physical observables under basis vectors of various aggregated receptor states. Our simulations of multivalent models not only reject Wofsy-Goldstein parameter conditions for cooperativity, but show that higher-order oligomer formations can shift cooperativity from positive to negative.

Computational Systems Biology in Cancer: Modeling Methods and Applications

In recent years it has become clear that carcinogenesis is a complex process, both at the molecular and cellular levels. Understanding the origins, growth and spread of cancer, therefore requires an integrated or system-wide approach. Computational systems biology is an emerging sub-discipline in systems biology that utilizes the wealth of data from genomic, proteomic and metabolomic studies to build computer simulations of intra and intercellular processes. Several useful descriptive and predictive models of the origin, growth and spread of cancers have been developed in an effort to better understand the disease and potential therapeutic approaches. In this review we describe and assess the practical and theoretical underpinnings of commonly-used modeling approaches, including ordinary and partial differential equations, petri nets, cellular automata, agent based models and hybrid systems. A number of computer-based formalisms have been implemented to improve the accessibility of the various approaches to researchers whose primary interest lies outside of model development. We discuss several of these and describe how they have led to novel insights into tumor genesis, growth, apoptosis, vascularization and therapy.

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.

Influences of excluded volume of molecules on signaling processes on the biomembrane

We investigate the influences of the excluded volume of molecules on biochemical reaction processes on 2-dimensional surfaces using a model of signal transduction processes on biomembranes. We perform simulations of the 2-dimensional cell-based model, which describes the reactions and diffusion of the receptors, signaling proteins, target proteins, and crowders on the cell membrane. The signaling proteins are activated by receptors, and these activated signaling proteins activate target proteins that bind autonomously from the cytoplasm to the membrane, and unbind from the membrane if activated. If the target proteins bind frequently, the volume fraction of molecules on the membrane becomes so large that the excluded volume of the molecules for the reaction and diffusion dynamics cannot be negligible. We find that such excluded volume effects of the molecules induce non-trivial variations of the signal flow, defined as the activation frequency of target proteins, as follows. With an increase in the binding rate of target proteins, the signal flow varies by monotonically increasing; increasing then decreasing in a bell-shaped curve; or increasing, decreasing, then increasing in an S-shaped curve. We further demonstrate that the excluded volume of molecules influences the hierarchical molecular distributions throughout the reaction processes. In particular, when the system exhibits a large signal flow, the signaling proteins tend to surround the receptors to form receptor-signaling protein clusters, and the target proteins tend to become distributed around such clusters. To explain these phenomena, we analyze the stochastic model of the local motions of molecules around the receptor.

The impact of high density receptor clusters on VEGF signaling

Vascular endothelial growth factor (VEGF) signaling is involved in the process of blood vessel development and maintenance. Signaling is initiated by binding of the bivalent VEGF ligand to the membrane-bound receptors (VEGFR), which in turn stimulates receptor dimerization. Herein, we discuss experimental evidence that VEGF receptors localize in caveloae and other regions of the plasma membrane, and for other receptors, it has been shown that receptor clustering has an impact on dimerization and thus also on signaling. Overall, receptor clustering is part of a complex ecosystem of interactions and how receptor clustering impacts dimerization is not well understood. To address these questions, we have formulated the simplest possible model. We have postulated the existence of a single high affinity region in the cell membrane, which acts as a transient trap for receptors. We have defined an ODE model by introducing high- and low-density receptor variables and introduce the corresponding reactions from a realistic model of VEGF signal initiation. Finally, we use the model to investigate the relation between the degree of VEGFR concentration, ligand availability, and signaling. In conclusion, our simulation results provide a deeper understanding of the role of receptor clustering in cell signaling.

Mathematical simulation of membrane protein clustering for efficient signal transduction

Initiation and propagation of cell signaling depend on productive interactions among signaling proteins at the plasma membrane. These diffusion-limited interactions can be influenced by features of the membrane that introduce barriers, such as cytoskeletal corrals, or microdomains that transiently confine both transmembrane receptors and membrane-tethered peripheral proteins. Membrane topographical features can lead to clustering of receptors and other membrane components, even under very dynamic conditions. This review considers the experimental and mathematical evidence that protein clustering impacts cell signaling in complex ways. Simulation approaches used to consider these stochastic processes are discussed.

The shuttling scaffold model for prevention of yeast pheromone pathway misactivation

The molecular scaffold in the yeast pheromone pathway, Ste5, shuttles continuously between the nucleus and the cytoplasm. Ste5 undergoes oligomerization reaction in the nucleus. Upon pheromone stimulation, the Ste5 dimer is rapidly exported out of the nucleus and recruited to the plasma membrane for pathway activation. This clever device on part of the yeast cell is thought to prevent pathway misactivation at high enough levels of Ste5 in the absence of pheromone. We have built a spatiotemporal model of signaling in this pathway to describe its regulation. Our present work underscores the importance of spatial modeling of cell signaling networks to understand their control and functioning.

Comparing the epidermal growth factor interaction with four different cell lines: intriguing effects imply strong dependency of cellular context

The interaction of the epidermal growth factor (EGF) with its receptor (EGFR) is known to be complex, and the common over-expression of EGF receptor family members in a multitude of tumors makes it important to decipher this interaction and the following signaling pathways. We have investigated the affinity and kinetics of (125)I-EGF binding to EGFR in four human tumor cell lines, each using four culturing conditions, in real time by use of LigandTracer®.Highly repeatable and precise measurements show that the overall apparent affinity of the (125)I-EGF - EGFR interaction is greatly dependent on cell line at normal culturing conditions, ranging from K(D) ≈ 200 pM on SKBR3 cells to K(D)≈8 nM on A431 cells. The (125)I-EGF - EGFR binding curves (irrespective of cell line) have strong signs of multiple simultaneous interactions. Furthermore, for the cell lines A431 and SKOV3, gefitinib treatment increases the (125)I-EGF - EGFR affinity, in particular when the cells are starved. The (125)I-EGF - EGFR interaction on cell line U343 is sensitive to starvation while as on SKBR3 it is insensitive to gefitinib and starvation.The intriguing pattern of the binding characteristics proves that the cellular context is important when deciphering how EGF interacts with EGFR. From a general perspective, care is advisable when generalizing ligand-receptor interaction results across multiple cell-lines.

Monte Carlo simulations of plasma membrane corral-induced EGFR clustering

Experimental evidence suggests that the cell membrane is a highly organized structure that is compartmentalized by the underlying membrane cytoskeleton (MSK). The interaction between the cell membrane and the cytoskeleton led to the "picket-fence" model, which was proposed to explain certain aspects of membrane compartmentalization. This model assumes that the MSK hinders and confines the motion of receptors and lipids to compartments in the membrane. However, the impact of the MSK on receptor clustering, aggregation, and downstream signaling remains unclear. For example, some evidence suggests that the MSK enhances dimerization, while other evidence suggests decreased dimerization and signaling. Herein, we use computational Monte Carlo simulations to examine the effects of MSK density and receptor concentration on receptor dimerization and clustering. Preliminary results suggest that the MSK may have the potential to induce receptor clustering, which is a function of both picket-fence density and receptor concentration.

Quantitative understanding of cell signaling: the importance of membrane organization

Systems biology modeling of signal transduction pathways traditionally employs ordinary differential equations, deterministic models based on the assumptions of spatial homogeneity. However, this can be a poor approximation for certain aspects of signal transduction, especially its initial steps: the cell membrane exhibits significant spatial organization, with diffusion rates approximately two orders of magnitude slower than those in the cytosol. Thus, to unravel the complexities of signaling pathways, quantitative models must consider spatial organization as an important feature of cell signaling. Furthermore, spatial separation limits the number of molecules that can physically interact, requiring stochastic simulation methods that account for individual molecules. Herein, we discuss the need for mathematical models and experiments that appreciate the importance of spatial organization in the membrane.

Adaptive coarse-grained Monte Carlo simulation of reaction and diffusion dynamics in heterogeneous plasma membranes

Background

An adaptive coarse-grained (kinetic) Monte Carlo (ACGMC) simulation framework is applied to reaction and diffusion dynamics in inhomogeneous domains. The presented model is relevant to the diffusion and dimerization dynamics of epidermal growth factor receptor (EGFR) in the presence of plasma membrane heterogeneity and specifically receptor clustering. We perform simulations representing EGFR cluster dissipation in heterogeneous plasma membranes consisting of higher density clusters of receptors surrounded by low population areas using the ACGMC method. We further investigate the effect of key parameters on the cluster lifetime.

Results

Coarse-graining of dimerization, rather than of diffusion, may lead to computational error. It is shown that the ACGMC method is an effective technique to minimize error in diffusion-reaction processes and is superior to the microscopic kinetic Monte Carlo simulation in terms of computational cost while retaining accuracy. The low computational cost enables sensitivity analysis calculations. Sensitivity analysis indicates that it may be possible to retain clusters of receptors over the time scale of minutes under suitable conditions and the cluster lifetime may depend on both receptor density and cluster size.

Conclusions

The ACGMC method is an ideal platform to resolve large length and time scales in heterogeneous biological systems well beyond the plasma membrane and the EGFR system studied here. Our results demonstrate that cluster size must be considered in conjunction with receptor density, as they synergistically affect EGFR cluster lifetime. Further, the cluster lifetime being of the order of several seconds suggests that any mechanisms responsible for EGFR aggregation must operate on shorter timescales (at most a fraction of a second), to overcome dissipation and produce stable clusters observed experimentally.

Coupled stochastic spatial and non-spatial simulations of ErbB1 signaling pathways demonstrate the importance of spatial organization in signal transduction

Background

The ErbB family of receptors activates intracellular signaling pathways that control cellular proliferation, growth, differentiation and apoptosis. Given these central roles, it is not surprising that overexpression of the ErbB receptors is often associated with carcinogenesis. Therefore, extensive laboratory studies have been devoted to understanding the signaling events associated with ErbB activation.

Methodology/Principal Findings

Systems biology has contributed significantly to our current understanding of ErbB signaling networks. However, although computational models have grown in complexity over the years, little work has been done to consider the spatial-temporal dynamics of receptor interactions and to evaluate how spatial organization of membrane receptors influences signaling transduction. Herein, we explore the impact of spatial organization of the epidermal growth factor receptor (ErbB1/EGFR) on the initiation of downstream signaling. We describe the development of an algorithm that couples a spatial stochastic model of membrane receptors with a nonspatial stochastic model of the reactions and interactions in the cytosol. This novel algorithm provides a computationally efficient method to evaluate the effects of spatial heterogeneity on the coupling of receptors to cytosolic signaling partners.

Conclusions/Significance

Mathematical models of signal transduction rarely consider the contributions of spatial organization due to high computational costs. A hybrid stochastic approach simplifies analyses of the spatio-temporal aspects of cell signaling and, as an example, demonstrates that receptor clustering contributes significantly to the efficiency of signal propagation from ligand-engaged growth factor receptors.

Therapeutically targeting ErbB3: a key node in ligand-induced activation of the ErbB receptor-PI3K axis

The signaling network downstream of the ErbB family of receptors has been extensively targeted by cancer therapeutics; however, understanding the relative importance of the different components of the ErbB network is nontrivial. To explore the optimal way to therapeutically inhibit combinatorial, ligand-induced activation of the ErbB-phosphatidylinositol 3-kinase (PI3K) axis, we built a computational model of the ErbB signaling network that describes the most effective ErbB ligands, as well as known and previously unidentified ErbB inhibitors. Sensitivity analysis identified ErbB3 as the key node in response to ligands that can bind either ErbB3 or EGFR (epidermal growth factor receptor). We describe MM-121, a human monoclonal antibody that halts the growth of tumor xenografts in mice and, consistent with model-simulated inhibitor data, potently inhibits ErbB3 phosphorylation in a manner distinct from that of other ErbB-targeted therapies. MM-121, a previously unidentified anticancer therapeutic designed using a systems approach, promises to benefit patients with combinatorial, ligand-induced activation of the ErbB signaling network that are not effectively treated by current therapies targeting overexpressed or mutated oncogenes.

Computational systems biology in cancer: modeling methods and applications

In recent years it has become clear that carcinogenesis is a complex process, both at the molecular and cellular levels. Understanding the origins, growth and spread of cancer, therefore requires an integrated or system-wide approach. Computational systems biology is an emerging sub-discipline in systems biology that utilizes the wealth of data from genomic, proteomic and metabolomic studies to build computer simulations of intra and intercellular processes. Several useful descriptive and predictive models of the origin, growth and spread of cancers have been developed in an effort to better understand the disease and potential therapeutic approaches. In this review we describe and assess the practical and theoretical underpinnings of commonly-used modeling approaches, including ordinary and partial differential equations, petri nets, cellular automata, agent based models and hybrid systems. A number of computer-based formalisms have been implemented to improve the accessibility of the various approaches to researchers whose primary interest lies outside of model development. We discuss several of these and describe how they have led to novel insights into tumor genesis, growth, apoptosis, vascularization and therapy.

Dimerization of VEGF receptors and implications for signal transduction: a computational study

Vascular endothelial growth factor (VEGF) is a potent cytokine involved in the induction of neovascularization. Secreted as a cysteine-linked dimer, it has two binding sites at opposite poles through which it may bind VEGF receptors (VEGFRs), receptor tyrosine kinases found on the surface of endothelial and other cells. The binding of a VEGF molecule to two VEGFR molecules induces transphosphorylation of the intracellular domains of the receptors, leading to signal transduction. The dominant mechanism of receptor dimerization is not clear: the receptors may be present in an inactive pre-dimerized form, VEGF binding first to one of the receptors, the second receptor then ideally located for dimerization; or VEGF may bind receptor monomers on the cell surface, which then diffuse and bind to available unligated receptor monomers to complete the activation. Both processes take place and one or other may dominate on different cell types. We demonstrate the impact of dimerization mechanism on the binding of VEGF to the cell surface and on the formation of active signaling receptor complexes. We describe two methods to determine which process dominates, based on binding and phosphorylation assays. The presence of two VEGF receptor populations, VEGFR1 and VEGFR2, can result in receptor heterodimer formation. Our simulations predict that heterodimers will comprise 10-50% of the active, signaling VEGF receptor complexes, and that heterodimers will form at the expense of homodimers of VEGFR1 when VEGFR2 populations are larger. These results have significant implications for VEGF signal transduction and interpretation of experimental studies. These results may be applicable to other ligand-receptor pairs, in particular PDGF.

The EGF receptor family: spearheading a merger of signaling and therapeutics

The ErbB receptor tyrosine kinases evolved as key regulatory entities enabling the extracellular milieu to communicate with the intracellular machinery to bring forth the appropriate biological response in an ever-changing environment. Since its discovery, many aspects of the ErbB family have been deciphered, with emphasis on aberration of signaling in human diseases. However, only now, with the availability of the atomic coordinates of these receptors, can we construct a comprehensive model of the mechanisms underlying ligand-induced receptor dimerization and subsequent tyrosine kinase activation. Furthermore, the recent introduction of new high-throughput screening methodologies, combined with the materialization of a systems biology perspective, reveals an overwhelming network complexity, enabling robust signaling and evolvability. This knowledge is likely to impact our view of diseases as system perturbations and resistance to ErbB-targeted therapeutics as manifestations of robustness.

Single-molecule analysis of epidermal growth factor binding on the surface of living cells

Global cellular responses induced by epidermal growth factor (EGF) receptor (EGFR) occur immediately with a less than 1% occupancy among tens of thousands of EGFR molecules on single cell surface. Activation of EGFR requires the formation of a signaling dimer of EGFR bound with a single ligand to each molecule. How sufficient numbers of signaling dimers are formed at such low occupancy rate is still not known. Here, we have analyzed the kinetics of EGF binding and the formation of the signaling dimer using single-molecule imaging and mathematical modeling. A small number of EGFR on the cell surface formed dimeric binding sites, which bound EGF two orders of magnitude faster than the monomeric binding sites. There was a positive cooperative binding of EGF to the dimeric binding sites through a newly discovered kinetic intermediate. These two mechanisms facilitate the formation of signaling dimers of EGF/EGFR complexes.

Spatial modeling of dimerization reaction dynamics in the plasma membrane: Monte Carlo vs. continuum differential equations

Bimolecular reactions in the plasma membrane, such as receptor dimerization, are a key signaling step for many signaling systems. For receptors to dimerize, they must first diffuse until a collision happens, upon which a dimerization reaction may occur. Therefore, study of the dynamics of cell signaling on the membrane may require the use of a spatial modeling framework. Despite the availability of spatial simulation methods, e.g., stochastic spatial Monte Carlo (MC) simulation and partial differential equation (PDE) based approaches, many biological models invoke well-mixed assumptions without completely evaluating the importance of spatial organization. Whether one is to utilize a spatial or non-spatial simulation framework is therefore an important decision. In order to evaluate the importance of spatial effects a priori, i.e., without performing simulations, we have assessed the applicability of a dimensionless number, known as second Damköhler number (Da), defined here as the ratio of time scales of collision and reaction, for 2-dimensional bimolecular reactions. Our study shows that dimerization reactions in the plasma membrane with Da approximately >0.1 (tested in the receptor density range of 10(2)-10(5)/microm(2)) require spatial modeling. We also evaluated the effective reaction rate constants of MC and simple deterministic PDEs. Our simulations show that the effective reaction rate constant decreases with time due to time dependent changes in the spatial distribution of receptors. As a result, the effective reaction rate constant of simple PDEs can differ from that of MC by up to two orders of magnitude. Furthermore, we show that the fluctuations in the number of copies of signaling proteins (noise) may also depend on the diffusion properties of the system. Finally, we used the spatial MC model to explore the effect of plasma membrane heterogeneities, such as receptor localization and reduced diffusivity, on the dimerization rate. Interestingly, our simulations show that localization of epidermal growth factor receptor (EGFR) can cause the diffusion limited dimerization rate to be up to two orders of magnitude higher at higher average receptor densities reported for cancer cells, as compared to a normal cell.