This site needs JavaScript to work properly. Please enable it to take advantage of the complete set of features!
Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation

Save citation to file

Add to Collections

Name must be less than 100 characters
Unable to load your collection due to an error
Please try again

Add to My Bibliography

Unable to load your delegates due to an error
Please try again

Your saved search

Would you like email updates of new search results?
Saved Search Alert Radio Buttons
()

Create a file for external citation management software

Your RSS Feed

Review
. 2014 Jan-Feb;6(1):13-36.
doi: 10.1002/wsbm.1245. Epub 2013 Sep 30.

Rule-based modeling: a computational approach for studying biomolecular site dynamics in cell signaling systems

Affiliations

Affiliation

  • 1 Department of Chemistry and Chemical Biology, Cornell University, Ithaca, NY, USA.
Review

Rule-based modeling: a computational approach for studying biomolecular site dynamics in cell signaling systems

Lily A Chylek et al. Wiley Interdiscip Rev Syst Biol Med. 2014 Jan-Feb.
. 2014 Jan-Feb;6(1):13-36.
doi: 10.1002/wsbm.1245. Epub 2013 Sep 30.

Affiliation

  • 1 Department of Chemistry and Chemical Biology, Cornell University, Ithaca, NY, USA.

Abstract

Rule-based modeling was developed to address the limitations of traditional approaches for modeling chemical kinetics in cell signaling systems. These systems consist of multiple interacting biomolecules (e.g., proteins), which themselves consist of multiple parts (e.g., domains, linear motifs, and sites of phosphorylation). Consequently, biomolecules that mediate information processing generally have the potential to interact in multiple ways, with the number of possible complexes and posttranslational modification states tending to grow exponentially with the number of binary interactions considered. As a result, only large reaction networks capture all possible consequences of the molecular interactions that occur in a cell signaling system, which is problematic because traditional modeling approaches for chemical kinetics (e.g., ordinary differential equations) require explicit network specification. This problem is circumvented through representation of interactions in terms of local rules. With this approach, network specification is implicit and model specification is concise. Concise representation results in a coarse graining of chemical kinetics, which is introduced because all reactions implied by a rule inherit the rate law associated with that rule. Coarse graining can be appropriate if interactions are modular, and the coarseness of a model can be adjusted as needed. Rules can be specified using specialized model-specification languages, and recently developed tools designed for specification of rule-based models allow one to leverage powerful software engineering capabilities. A rule-based model comprises a set of rules, which can be processed by general-purpose simulation and analysis tools to achieve different objectives (e.g., to perform either a deterministic or stochastic simulation).

PubMed Disclaimer

Figures

Figure 1

Figure 1

Illustration of proteins and interactions…

Figure 1

Illustration of proteins and interactions involved in EGFR signaling. (A) In this graph,…

Figure 1
Illustration of proteins and interactions involved in EGFR signaling. (A) In this graph, nodes correspond to proteins and edges correspond to direct physical interactions. This type of graph is commonly used to visualize protein interaction networks. Numbers next to edges refer to descriptions of interactions given in Table 1. (B) In this graph, arrows represent positive and negative influences. Note that there is a (negative) arrow connecting PTPN11 and RASA1 even though these proteins do not directly interact (cf. Panels A and B) because PTPN11 is responsible for dephosphorylation of a site in EGFR that interacts with RASA1 (Table 1). Interesting regulatory circuits embedded within the diagram of Panel B are highlighted in the panels at bottom. (C) EGFR generates competing positive and negative signals for KRAS activation. (D) An incoherent type 1 feed-forward loop (FFL) motif. (E) A coherent type 1 FFL motif. (F) A positive feedback loop.
Figure 2

Figure 2

The edges/arrows of a typical…

Figure 2

The edges/arrows of a typical pathway diagram hide significant combinatorial complexity. (A) The…

Figure 2
The edges/arrows of a typical pathway diagram hide significant combinatorial complexity. (A) The interactions and phosphorylation states of the three proteins considered in this diagram (EGFR and two adaptor proteins, GRB2 and SHC1) imply 11 distinct chemical species. The horizontal bars serve to label the 11 species and to delineate their compositions. Note that EGFR is taken to be phosphorylated at two sites of EGFR autophosphorylation: Y1092 (or both Y1092 and Y1110, which are important docking sites of GRB2) and Y1172 (or both Y1172 and Y1197, which are important docking sites of SHC1). GRB2 and SHC1 interact with EGFR via domains that recognize phosphotyrosines: the SH2 domain in GRB2 and the PTB domain in SHC1. The diagram shown here does not comprehensively depict the known interactions and phosphorylation states of these proteins; for example, EGFR-mediated phosphorylation of SHC1 and interaction between GRB2 and phosphorylated SHC1 are not considered. (B) The 11 chemical species are connected in a reaction network encompassing 24 unidirectional reactions. Reactions 1–6 are phosphorylation reactions, Reactions 7–12 are dephosphorylation reactions, Reactions 13–18 are bimolecular association reactions, and Reactions 19–24 are dissociation reactions. For simplicity, all reactions are represented as single-step transformations and catalysts are considered implicitly. For more information about the interactions considered here, see Table 1.
Figure 3

Figure 3

A protein or protein complex…

Figure 3

A protein or protein complex can be represented at different levels of resolution.

Figure 3
A protein or protein complex can be represented at different levels of resolution. (A) An atomic model of EGFR (amino acid residues 25–982). The model is based on experimentally determined structures of three domains of EGFR (PDB accession codes 3NJP, 2M20, and 2JIV). The model does not include the cytoplasmic tail of EGFR. (B) An atomic model of an EGF-induced oligomer of EGFR, a cyclic side-by-side dimer of dimers. EGF is not shown. The model was constructed with preservation of experimentally determined structural interfaces between domains and with preservation of the chemical integrity of individual domains, which were treated as rigid bodies. Consecutive domains are within the allowed lengths of the connecting loops. The top view highlights ectodomain-ectodomain interactions, and the bottom view highlights cytoplasmic domain-cytoplasmic domain interactions. (C) Domain architecture of EGFR. This diagram provides a schematic representation of EGFR and its component parts and includes a depiction of five sites of autophosphorylation (cf. Figs. 2 and 5). (D) BNGL-encoded representation of EGF. This encoding indicates that the molecule type EGF contains one functional component, rec. More formally, this line of code introduces a graph that has the color EGF and one vertex labeled rec. (E) BNGL-encoded representation of EGFR (cf. Panel A). This encoding indicates that the molecule type EGFR has four functional components: lig, back, cd, and Y. The lig component represents the ligand-binding site, which comprises domains I and III of the EGFR ectodomain. The back component represents domain II of the EGFR ectodomain, which is responsible for a self-interaction. The cd component, which is responsible for another self-interaction, represents the cytoplasmic domain of EGFR and encompasses the juxtramembrane region and kinase domain. The cd component is taken to have two possible internal states (or more formally, vertex attributes): closed (c) and open (o). In the o state, the cd component is able to interact with a second copy of itself (also in the o state). The Y component represents an autophosphorylation site, which is taken to have two possible internal states. These states, u and p, represent different phosphorylation states, unphosphorylated and phosphorylated. (F) BNGL-encoded representation of a cyclic EGFR tetramer (cf. Panel B). As indicated, this complex is held together through alternating back-back and cd-cd interactions, which are abstractions of the interactions illustrated at atomic resolution in Panel B. Atomic models were visualized using the UCSF Chimera package. The model of Panel B represents only one plausible structure.
Figure 4

Figure 4

A rule-based model for IGF1…

Figure 4

A rule-based model for IGF1 interaction with IGF1R. (A) Two reactions from the…

Figure 4
A rule-based model for IGF1 interaction with IGF1R. (A) Two reactions from the model of Kiselyov et al.. Reversible capture of a free ligand (IGF1) by an inactive, unbound receptor (IGF1R) is illustrated at left, and reversible ligand-mediated crosslinking of receptor subunits containing Sites 1 and 2 (or S1 and S2) is illustrated at right. A receptor is dimeric and each monomer contains two ligand-binding sites. Common shading indicates sites that are found in the same receptor monomer. (B) Two rules from the model of Kiselyov et al. drawn according to the graphical conventions of Faeder et al.. The rule at top (bottom) implies the left (right) reaction in Panel A. (C) Excerpts from an executable encoding of the model of Kiselyov et al. in BNGL. In the listing shown here, definitions of graphs for IGF1 and IGF1R (i.e., molecule types) are given at top, definitions of seed species are given next, definition of a pattern used in a local function is given in the observables block, and definitions of rules are given at bottom. Rule 1 is the text encoding of the top rule in Panel B, and Rule 3 is the text encoding of the bottom rule in Panel B. Rule 2 is similar to Rule 1. Rules 3 and 4 together define crosslinking of Sites 1 and 2 when a ligand is bound at Site 1. Rule 5 defines crosslinking of Sites 1 and 2 when a ligand is bound at Site 2.
Figure 5

Figure 5

Example of a contact map…

Figure 5

Example of a contact map for the model of Kiselyov et al. (Fig.…

Figure 5
Example of a contact map for the model of Kiselyov et al. (Fig. 4). The map shown here is essentially the same as the one that is generated automatically by RuleBender, from the model specification that is partially shown in Fig. 4C.
Figure 6

Figure 6

Example of an extended contact…

Figure 6

Example of an extended contact map for selected proteins and protein-protein interactions involved…

Figure 6
Example of an extended contact map for selected proteins and protein-protein interactions involved in EGFR signaling. The same proteins and interactions are considered in Fig. 1. Boxes with rounded corners represent selected proteins, domains, and linear motifs. Small square boxes attached to vertical lines represent sites of phosphorylation. Lines that begin and end with an arrowhead represent noncovalent direct-binding interactions. Lines that begin at a box representing a catalytic subunit and end with a circle point to substrates of an enzyme (a kinase or phosphatase). An open circle indicates a post-translational modification (phosphorylation); a circle with a line through it indicates reversal of a modification (dephosphorylation). See Table 1 for more information about the interactions represented in this map.
Figure 7

Figure 7

Illustration of how PySB can…

Figure 7

Illustration of how PySB can simplify the specification of related but distinct rules.

Figure 7
Illustration of how PySB can simplify the specification of related but distinct rules. (A) PySB code that specifies two phosphorylation rules for each of seven sites of IGF1R autophosphorylation. IGF1R is a receptor tyrosine kinase. The rules all have the same form, differing only with respect to substrate and symmetry of ligand binding. It is assumed that an IGF1-crosslinked receptor can mediate autophosphorylation of sites in both of its monomer subunits. (B) Examples of the 14 rules defined by the PySB code in Panel A, in BNGL format.

References

    1. Bhalla US, Iyengar R. Emergent properties of networks of biological signaling pathways. Science. 1999;283:381–387. - PubMed
    1. Kitano H. Systems biology: a brief overview. Science. 2002;295:1662–1664. - PubMed
    1. Nurse P. Life, logic and information. Nature. 2008;454:424–426. - PubMed
    1. Hlavacek WS, Faeder JR, Blinov ML, Posner RG, Hucka M, Fontana W. Rules for modeling signal-transduction systems. Sci STKE. 2006:re6. - PubMed
    1. Kholodenko B, Yaffe MB, Kolch W. Computational approaches for analyzing information flow in biological networks. Sci Signal. 2012;5:re1. - PubMed

Publication types

LinkOut - more resources

Cite
Morty Proxy This is a proxified and sanitized view of the page, visit original site.