Cell signaling in space and time: where proteins come together and when they're apart

Attributes of short linear motifs

Traditionally, protein-protein interactions were thought to be mediated by large, structured domains. However, it has become clear that the interactome comprises a wide range of binding interfaces with varying degrees of flexibility, ranging from rigid globular domains to disordered regions that natively lack structure. Enrichment for disorder in highly connected hub proteins and its correlation with organism complexity hint at the functional importance of disordered regions. Nevertheless, they have not yet been extensively characterised. Shifting the attention from globular domains to disordered regions of the proteome might bring us closer to elucidating the dense and complex connectivity of the interactome. An important class of disordered interfaces are the compact mono-partite, short linear motifs (SLiMs, or eukaryotic linear motifs (ELMs)). They are evolutionarily plastic and interact with relatively low affinity due to the limited number of residues that make direct contact with the binding partner. These features confer to SLiMs the ability to evolve convergently and mediate transient interactions, which is imperative to network evolution and to maintain robust cell signalling, respectively. The ability to discriminate biologically relevant SLiMs by means of different attributes will improve our understanding of the complexity of the interactome and aid development of bioinformatics tools for motif discovery. In this paper, the curated instances currently available in the Eukaryotic Linear Motif (ELM) database are analysed to provide a clear overview of the defining attributes of SLiMs. These analyses suggest that functional SLiMs have higher levels of conservation than their surrounding residues, frequently evolve convergently, preferentially occur in disordered regions and often form a secondary structure when bound to their interaction partner. These results advocate searching for small groupings of residues in disordered regions with higher relative conservation and a propensity to form the secondary structure. Finally, the most interesting conclusions are examined in regard to their functional consequences.

Synthetic mammalian gene networks as a blueprint for the design of interactive biohybrid materials

Synthetic biology aims at the rational design and construction of devices, systems and organisms with desired functionality based on modular well-characterized biological building blocks. Based on first proof-of-concept studies in bacteria a decade ago, synthetic biology strategies have rapidly entered mammalian cell technology providing novel therapeutic solutions. Here we review how biological building blocks can be rewired to interactive regulatory genetic networks in mammalian cells and how these networks can be transformed into open- and closed-loop control configurations for autonomously managing disease phenotypes. In the second part of this tutorial review we describe how the regulatory biological sensors and switches can be transferred from mammalian cell synthetic biology to materials sciences in order to develop interactive biohybrid materials with similar (therapeutic) functionality as their synthetic biological archetypes. We develop a perspective of how the convergence of synthetic biology with materials sciences might contribute to the development of truly interactive and adaptive materials for autonomous operation in a complex environment.

AKAP220 protein organizes signaling elements that impact cell migration

Cell movement requires the coordinated reception, integration, and processing of intracellular signals. We have discovered that the protein kinase A anchoring protein AKAP220 interacts with the cytoskeletal scaffolding protein IQGAP1 to influence cell motility. AKAP220/IQGAP1 networks receive and integrate calcium and cAMP second messenger signals and position signaling enzymes near their intended substrates at leading edges of migrating cells. IQGAP1 supports calcium/calmodulin-dependent association of factors that modulate microtubule dynamics. AKAP220 suppresses GSK-3β and positions this kinase to allow recruitment of the plus-end microtubule tracking protein CLASP2. Gene silencing of AKAP220 alters the rate of microtubule polymerization and the lateral tracking of growing microtubules and retards cell migration in metastatic human cancer cells. This reveals an unappreciated role for this anchored kinase/microtubule effector protein network in the propagation of cell motility.

Expanding the genetic code of Escherichia coli with phosphoserine

O-Phosphoserine (Sep), the most abundant phosphoamino acid in the eukaryotic phosphoproteome, is not encoded in the genetic code, but synthesized posttranslationally. Here, we present an engineered system for specific cotranslational Sep incorporation (directed by UAG) into any desired position in a protein by an Escherichia coli strain that harbors a Sep-accepting transfer RNA (tRNA(Sep)), its cognate Sep-tRNA synthetase (SepRS), and an engineered EF-Tu (EF-Sep). Expanding the genetic code rested on reengineering EF-Tu to relax its quality-control function and permit Sep-tRNA(Sep) binding. To test our system, we synthesized the activated form of human mitogen-activated ERK activating kinase 1 (MEK1) with either one or two Sep residues cotranslationally inserted in their canonical positions (Sep(218), Sep(222)). This system has general utility in protein engineering, molecular biology, and disease research.

Regulation of G protein-coupled receptor function by Na+/H+ exchange regulatory factors

Many G protein-coupled receptors (GPCR) exert patterns of cell-specific signaling and function. Mounting evidence now supports the view that cytoplasmic adapter proteins contribute critically to this behavior. Adapter proteins recognize highly conserved motifs such as those for Src homology 3 (SH3), phosphotyrosine-binding (PTB), and postsynaptic density 95/discs-large/zona occludens (PDZ) docking sequences in candidate GPCRs. Here we review the behavior of the Na+/H+ exchange regulatory factor (NHERF) family of PDZ adapter proteins on GPCR signalling, trafficking, and function. Structural determinants of NHERF proteins that allow them to recognize targeted GPCRs are considered. NHERF1 and NHERF2 are capable also of modifying the assembled complex of accessory proteins such as β-arrestins, which have been implicated in regulating GPCR signaling. In addition, NHERF1 and NHERF2 modulate GPCR signaling by altering the G protein to which the receptor binds or affect other regulatory proteins that affect GTPase activity, protein kinase A, phospholipase C, or modify downstream signaling events. Small molecules targeting the site of NHERF1-GPCR interaction are being developed and may become important and selective drug candidates.

Differential inhibitor sensitivity between human kinases VRK1 and VRK2

Human vaccinia-related kinases (VRK1 and VRK2) are atypical active Ser-Thr kinases implicated in control of cell cycle entry, apoptosis and autophagy, and affect signalling by mitogen activated protein kinases (MAPK). The specific structural differences in VRK catalytic sites make them suitable candidates for development of specific inhibitors. In this work we have determined the sensitivity of VRK1 and VRK2 to kinase inhibitors, currently used in biological assays or in preclinical studies, in order to discriminate between the two proteins as well as with respect to the vaccinia virus B1R kinase. Both VRK proteins and vaccinia B1R are poorly inhibited by inhibitors of different types targeting Src, MEK1, B-Raf, JNK, p38, CK1, ATM, CHK1/2 and DNA-PK, and most of them have no effect even at 100 µM. Despite their low sensitivity, some of these inhibitors in the low micromolar range are able to discriminate between VRK1, VRK2 and B1R. VRK1 is more sensitive to staurosporine, RO-31-8220 and TDZD8. VRK2 is more sensitive to roscovitine, RO 31–8220, Cdk1 inhibitor, AZD7762, and IC261. Vaccinia virus B1R is more sensitive to staurosporine, KU55933, and RO 31–8220, but not to IC261. Thus, the three kinases present a different pattern of sensitivity to kinase inhibitors. This differential response to known inhibitors can provide a structural framework for VRK1 or VRK2 specific inhibitors with low or no cross-inhibition. The development of highly specific VRK1 inhibitors might be of potential clinical use in those cancers where these kinases identify a clinical subtype with a poorer prognosis, as is the case of VRK1 in breast cancer.

Distinct involvement of the Gab1 and Grb2 adaptor proteins in signal transduction by the related receptor tyrosine kinases RON and MET

Although the signal transduction mechanisms of the receptor tyrosine kinase MET are well defined, less is known about its close relative RON. MET initiates intracellular signaling by autophosphorylation on specific cytoplasmic tyrosines that form docking sites for the adaptor proteins Grb2 and Gab1. Grb2 binds directly and is essential for all of the biological activities of MET. Gab1 docks either directly or indirectly via Grb2 and controls only a subset of MET functions. Because MET and RON possess similar adaptor binding sites, it was anticipated that their adaptor interactions would be conserved. Here we show that in contrast to MET, RON relies primarily on Gab1 for signal transmission. Surprisingly, disruption of the Grb2 docking site of RON or Grb2 depletion augments activity, whereas enhancement of Grb2 binding attenuates Gab1 recruitment and signaling. Hence, RON and MET differ in their adaptor interactions; furthermore, Grb2 performs a novel antagonistic role in the context of RON signaling.

Signals of LOV1: a computer simulation study on the wildtype LOV1-domain of Chlamydomonas reinhardtii and its mutants

Phototropins are photoreceptors regulating the blue-light response in plants and bacteria. They consist of two LOV (light oxygen voltage sensitive) domains each containing a non-covalently bound flavin-mononucleotide (FMN) chromophore, which are connected to a serine/threonine-kinase. Upon illumination, the LOV-domains undergo conformational changes, triggering a signal cascade in the organism through kinase activation. Here, we present results from molecular dynamics simulations in which we investigate the signal transduction pathway of the wildtype LOV1-domain of Chlamydomonas reinhardtii and a methyl-mercaptan (MM) adduct of its Cys57Gly-mutant at the molecular level. In particular, we analyzed the effect of covalent-bond formation between the reactive cysteine Cys57 and the FMN-reaction center, as well as the subsequent charge redistribution, on the spatio-dynamical behavior of the LOV1-domain. We compare the calculation results with experimental data and demonstrate that these adduct state characteristics have an important influence on the response of this photosensor. The light-induced changes implicate primarily an alteration of the surface charge distribution through rearrangement of the highly flexible Cα-, Dα- and Eα-helices including the Glu51-Lys91-salt bridge on the hydrophilic side of the protein domain and a β-sheet tightening process via coupling of the Aβ- and Bβ-strands. Our findings confirm the aptitude of the LOV1-domain to function as a dimerization partner, allowing the green alga to adapt its reproduction and growth speed to the environmental conditions.

Host-pathogen interactions: cheating the host by making new connections

Dynamic signaling networks are required to perform complex cellular processes. Structural and functional data now indicate the intriguing possibility that extracellular bacterial pathogens use catalytic scaffolds to assemble unique supramolecular signaling networks that effectively subvert key cellular processes in the host.

NMR analysis of a stress response metabolic signaling network

We previously hypothesized that Staphylococcus epidermidis senses a diverse set of environmental and nutritional factors associated with biofilm formation through a modulation in the activity of the tricarboxylic acid (TCA) cycle. Herein, we report our further investigation of the impact of additional environmental stress factors on TCA cycle activity and provide a detailed description of our NMR methodology. S. epidermidis wild-type strain 1457 was treated with stressors that are associated with biofilm formation, a sublethal dose of tetracycline, 5% NaCl, 2% glucose, and autoinducer-2 (AI-2). As controls and to integrate our current data with our previous study, 4% ethanol stress and iron-limitation were also used. Consistent with our prior observations, the effect of many environmental stress factors on the S. epidermidis metabolome was essentially identical to the effect of TCA cycle inactivation in the aconitase mutant strain 1457-acnA::tetM. A detailed quantitative analysis of metabolite concentration changes using 2D (1)H-(13)C HSQC and (1)H-(1)H TOCSY spectra identified a network of 37 metabolites uniformly affected by the stressors and TCA cycle inactivation. We postulate that the TCA cycle acts as the central pathway in a metabolic signaling network.

The cell in an era of systems biology

The increasing use of high-throughput technologies and computational modeling is revealing new levels of biological function and organization. How are these features of systems biology influencing our view of the cell?

Integrating cardiac PIP3 and cAMP signaling through a PKA anchoring function of p110γ

Adrenergic stimulation of the heart engages cAMP and phosphoinositide second messenger signaling cascades. Cardiac phosphoinositide 3-kinase p110γ participates in these processes by sustaining β-adrenergic receptor internalization through its catalytic function and by controlling phosphodiesterase 3B (PDE3B) activity via an unknown kinase-independent mechanism. We have discovered that p110γ anchors protein kinase A (PKA) through a site in its N-terminal region. Anchored PKA activates PDE3B to enhance cAMP degradation and phosphorylates p110γ to inhibit PIP(3) production. This provides local feedback control of PIP(3) and cAMP signaling events. In congestive heart failure, p110γ is upregulated and escapes PKA-mediated inhibition, contributing to a reduction in β-adrenergic receptor density. Pharmacological inhibition of p110γ normalizes β-adrenergic receptor density and improves contractility in failing hearts.

Protein pathway activation mapping of brain metastasis from lung and breast cancers reveals organ type specific drug target activation

Brain metastases are the most common fatal complication of systemic cancer, especially of lung (40-50%) and breast (20-30%) cancers. In this era of personalized therapy, there is a critical need to uncover the signaling architecture of brain metastases; however, little is known about what signaling pathways are activated in the context of the brain microenvironment. Using a unique study set of 42 brain metastases from patients with breast or nonsmall cell lung cancer (NSCLC), the phosphorylation/activation states of 128 key signaling proteins involved in cancer signaling were measured in laser capture microdissected tumor epithelium using reverse phase protein microarray (RPMA) technology. Distinct pathway activation subgroups from both breast and lung metastases were underpinned by, among others, ERBB2, AKT, mTOR, EGFR, SMAD, and ERK-p38 signaling. Breast cancer metastases showed significantly (p < 0.05) higher activation of the c-ERBB2/IGFR-AKT pathway network compared to NSCLC metastases, whereas NSCLC metastases to the brain exhibited higher relative levels of many members of the EGFR-ERK signaling network. Protein pathway activation mapping using RPMA revealed both the heterogeneity of signaling networks in brain metastases that would require a prior stratification to targeted therapies as well as the requirement of direct analysis of the metastatic lesion.

A roadmap to generate renewable protein binders to the human proteome

Despite the wealth of commercially available antibodies to human proteins, research is often hindered by their inconsistent validation, their poor performance and the inadequate coverage of the proteome. These issues could be addressed by systematic, genome-wide efforts to generate and validate renewable protein binders. We report a multicenter study to assess the potential of hybridoma and phage-display technologies in a coordinated large-scale antibody generation and validation effort. We produced over 1,000 antibodies targeting 20 SH2 domain proteins and evaluated them for potency and specificity by enzyme-linked immunosorbent assay (ELISA), protein microarray and surface plasmon resonance (SPR). We also tested selected antibodies in immunoprecipitation, immunoblotting and immunofluorescence assays. Our results show that high-affinity, high-specificity renewable antibodies generated by different technologies can be produced quickly and efficiently. We believe that this work serves as a foundation and template for future larger-scale studies to create renewable protein binders.

Regulation of human EGF receptor by lipids

The human epidermal growth factor receptor (EGFR) is a key representative of tyrosine kinase receptors, ubiquitous actors in cell signaling, proliferation, differentiation, and migration. Although the receptor is well-studied, a central issue remains: How does the compositional diversity and functional diversity of the surrounding membrane modulate receptor function? Reconstituting human EGFR into proteoliposomes of well-defined and controlled lipid compositions represents a minimal synthetic approach to systematically address this question. We show that lipid composition has little effect on ligand-binding properties of the EGFR but rather exerts a profound regulatory effect on kinase domain activation. Here, the ganglioside GM3 but not other related lipids strongly inhibited the autophosphorylation of the EGFR kinase domain. This inhibitory action of GM3 was only seen in liposomes compositionally poised to phase separate into coexisting liquid domains. The inhibition by GM3 was released by either removing the neuraminic acid of the GM3 headgroup or by mutating a membrane proximal lysine of EGFR (K642G). Our results demonstrate that GM3 exhibits the potential to regulate the allosteric structural transition from inactive to a signaling EGFR dimer, by preventing the autophosphorylation of the intracellular kinase domain in response to ligand binding.

Scaffold proteins: hubs for controlling the flow of cellular information

The spatial and temporal organization of molecules within a cell is critical for coordinating the many distinct activities carried out by the cell. In an increasing number of biological signaling processes, scaffold proteins have been found to play a central role in physically assembling the relevant molecular components. Although most scaffolds use a simple tethering mechanism to increase the efficiency of interaction between individual partner molecules, these proteins can also exert complex allosteric control over their partners and are themselves the target of regulation. Scaffold proteins offer a simple, flexible strategy for regulating selectivity in pathways, shaping output behaviors, and achieving new responses from preexisting signaling components. As a result, scaffold proteins have been exploited by evolution, pathogens, and cellular engineers to reshape cellular behavior.

Control of host cell phosphorylation by legionella pneumophila

Phosphorylation is one of the most frequent modifications in intracellular signaling and is implicated in many processes ranging from transcriptional control to signal transduction in innate immunity. Many pathogens modulate host cell phosphorylation pathways to promote growth and establish an infectious disease. The intracellular pathogen Legionella pneumophila targets and exploits the host phosphorylation system throughout the infection cycle as part of its strategy to establish an environment beneficial for replication. Key to this manipulation is the L. pneumophila Icm/Dot type IV secretion system, which translocates bacterial proteins into the host cytosol that can act directly on phosphorylation cascades. This review will focus on the different stages of L. pneumophila infection, in which host kinases and phosphatases contribute to infection of the host cell and promote intracellular survival of the pathogen. This includes the involvement of phosphatidylinositol 3-kinases during phagocytosis as well as the role of phosphoinositide metabolism during the establishment of the replication vacuole. Furthermore, L. pneumophila infection modulates the NF-κB and mitogen-activated protein kinase pathways, two signaling pathways that are central to the host innate immune response and involved in regulation of host cell survival. Therefore, L. pneumophila infection manipulates host cell signal transduction by phosphorylation at multiple levels.

Anchored protein kinase A recruitment of active Rac GTPase

Protein kinase A-anchoring proteins (AKAPs) influence fundamental cellular processes by directing the cAMP-dependent protein kinase (PKA) toward its intended substrates. In this report we describe the identification and characterization of a ternary complex of AKAP220, the PKA holoenzyme, and the IQ domain GTPase-activating protein 2 isoform (IQGAP2) that is enriched at cortical regions of the cell. Formation of an IQGAP2-AKAP220 core complex initiates a subsequent phase of protein recruitment that includes the small GTPase Rac. Biochemical and molecular biology approaches reveal that PKA phosphorylation of Thr-716 on IQGAP2 enhances association with the active form of the Rac GTPase. Cell-based experiments indicate that overexpression of an IQGAP2 phosphomimetic mutant (IQGAP2 T716D) enhances the formation of actin-rich membrane ruffles at the periphery of HEK 293 cells. In contrast, expression of a nonphosphorylatable IQGAP2 T716A mutant or gene silencing of AKAP220 suppresses formation of membrane ruffles. These findings imply that IQGAP2 and AKAP220 act synergistically to sustain PKA-mediated recruitment of effectors such as Rac GTPases that impact the actin cytoskeleton.

Cross-talk between carboxypeptidase M and the kinin B1 receptor mediates a new mode of G protein-coupled receptor signaling

G protein-coupled receptor (GPCR) signaling is affected by formation of GPCR homo- or heterodimers, but GPCR regulation by other cell surface proteins is not well understood. We reported that the kinin B1 receptor (B1R) heterodimerizes with membrane carboxypeptidase M (CPM), facilitating receptor signaling via CPM-mediated conversion of bradykinin or kallidin to des-Arg kinin B1R agonists. Here, we found that a catalytically inactive CPM mutant that still binds substrate (CPM-E264Q) also facilitates efficient B1R signaling by B2 receptor agonists bradykinin or kallidin. This response required co-expression of B1R and CPM-E264Q in the same cell, was disrupted by antibody that dissociates CPM from B1R, and was not found with a CPM-E264Q-B1R fusion protein. An additional mutation that reduced the affinity of CPM for C-terminal Arg and increased the affinity for C-terminal Lys inhibited the B1R response to bradykinin (with C-terminal Arg) but generated a response to Lys(9)-bradykinin. CPM-E264Q-mediated activation of B1Rs by bradykinin resulted in increased intramolecular fluorescence resonance energy transfer (FRET) in a B1R FRET construct, similar to that generated directly by a B1R agonist. In cytokine-treated human lung microvascular endothelial cells, disruption of B1R-CPM heterodimers inhibited B1R-dependent NO production stimulated by bradykinin and blocked the increased endothelial permeability caused by treatment with bradykinin and pyrogallol (a superoxide generator). Thus, CPM and B1Rs on cell membranes form a critical complex that potentiates B1R signaling. Kinin peptide binding to CPM causes a conformational change in the B1R leading to intracellular signaling and reveals a new mode of GPCR activation by a cell surface peptidase.

Plugging PKA into ERK scaffolds

Cancers often arise in part through derangements in protein kinase signaling. A striking example of this is the finding that approximately 30% of human tumors have mutations in Ras or B-Raf, leading to aberrant ERK kinase activation. Kinase signaling networks are often organized by scaffolding and anchoring proteins that help shape the dynamics of signal processing. AKAP-Lbc associates with the ERK scaffold protein KSR-1 to organize a growth factor and cAMP responsive signaling network. AKAP-Lbc also directs PKA phosphorylation of KSR-1 on a critical residue to ensure maximal signaling efficiency.

Exploring the spatial and temporal organization of a cell's proteome

To increase our current understanding of cellular processes, such as cell signaling and division, knowledge is needed about the spatial and temporal organization of the proteome at different organizational levels. These levels cover a wide range of length and time scales: from the atomic structures of macromolecules for inferring their molecular function, to the quantitative description of their abundance, and spatial distribution in the cell. Emerging new experimental technologies are greatly increasing the availability of such spatial information on the molecular organization in living cells. This review addresses three fields that have significantly contributed to our understanding of the proteome's spatial and temporal organization: first, methods for the structure determination of individual macromolecular assemblies, specifically the fitting of atomic structures into density maps generated from electron microscopy techniques; second, research that visualizes the spatial distributions of these complexes within the cellular context using cryo electron tomography techniques combined with computational image processing; and third, methods for the spatial modeling of the dynamic organization of the proteome, specifically those methods for simulating reaction and diffusion of proteins and complexes in crowded intracellular fluids. The long-term goal is to integrate the varied data about a proteome's organization into a spatially explicit, predictive model of cellular processes.

Local control of β-adrenergic stimulation: Effects on ventricular myocyte electrophysiology and Ca(2+)-transient

Local signaling domains and numerous interacting molecular pathways and substrates contribute to the whole-cell response of myocytes during β-adrenergic stimulation (βARS). We aimed to elucidate the quantitative contribution of substrates and their local signaling environments during βARS to the canine epicardial ventricular myocyte electrophysiology and calcium transient (CaT). We present a computational compartmental model of βARS and its electrophysiological effects. Novel aspects of the model include localized signaling domains, incorporation of β1 and β2 receptor isoforms, a detailed population-based approach to integrate the βAR and Ca(2+)/Calmodulin kinase (CaMKII) signaling pathways and their effects on a wide range of substrates that affect whole-cell electrophysiology and CaT. The model identifies major roles for phosphodiesterases, adenylyl cyclases, PKA and restricted diffusion in the control of local cAMP levels and shows that activation of specific cAMP domains by different receptor isoforms allows for specific control of action potential and CaT properties. In addition, the model predicts increased CaMKII activity during βARS due to rate-dependent accumulation and increased Ca(2+) cycling. CaMKII inhibition, reduced compartmentation, and selective blockade of β1AR is predicted to reduce the occurrence of delayed afterdepolarizations during βARS. Finally, the relative contribution of each PKA substrate to whole-cell electrophysiology is quantified by comparing simulations with and without phosphorylation of each target. In conclusion, this model enhances our understanding of localized βAR signaling and its whole-cell effects in ventricular myocytes by incorporating receptor isoforms, multiple pathways and a detailed representation of multiple-target phosphorylation; it provides a basis for further studies of βARS under pathological conditions.

Kinase associated-1 domains drive MARK/PAR1 kinases to membrane targets by binding acidic phospholipids

Phospholipid-binding modules such as PH, C1, and C2 domains play crucial roles in location-dependent regulation of many protein kinases. Here, we identify the KA1 domain (kinase associated-1 domain), found at the C terminus of yeast septin-associated kinases (Kcc4p, Gin4p, and Hsl1p) and human MARK/PAR1 kinases, as a membrane association domain that binds acidic phospholipids. Membrane localization of isolated KA1 domains depends on phosphatidylserine. Using X-ray crystallography, we identified a structurally conserved binding site for anionic phospholipids in KA1 domains from Kcc4p and MARK1. Mutating this site impairs membrane association of both KA1 domains and intact proteins and reveals the importance of phosphatidylserine for bud neck localization of yeast Kcc4p. Our data suggest that KA1 domains contribute to "coincidence detection," allowing kinases to bind other regulators (such as septins) only at the membrane surface. These findings have important implications for understanding MARK/PAR1 kinases, which are implicated in Alzheimer's disease, cancer, and autism.

Evolving specificity from variability for protein interaction domains

An important question in modular domain-peptide interactions, which play crucial roles in many biological processes, is how the diverse specificities exhibited by different members of a domain family are encoded in a common scaffold. Analysis of the Src homology (SH) 2 family has revealed that its specificity is determined, in large part, by the configuration of surface loops that regulate ligand access to binding pockets. In a distinct manner, SH3 domains employ loops for ligand recognition. The PDZ domain, in contrast, achieves specificity by co-evolution of binding-site residues. Thus, the conformational and sequence variability afforded by surface loops and binding sites provides a general mechanism by which to encode the wide spectrum of specificities observed for modular protein interaction domains.

The nerve growth factor reduces APOBEC3G synthesis and enhances HIV-1 transcription and replication in human primary macrophages

Macrophages infected with HIV-1 sustain viral replication for long periods of time, functioning as viral reservoirs. Therefore, recognition of factors that maintain macrophage survival and influence HIV-1 replication is critical to understanding the mechanisms that regulate the HIV-1-replicative cycle. Because HIV-1-infected macrophages release the nerve growth factor (NGF), and NGF neutralization reduces viral production, we further analyzed how this molecule affects HIV-1 replication. In the present study, we show that NGF stimulates HIV-1 replication in primary macrophages by signaling through its high-affinity receptor Tropomyosin-related Kinase A (TrKA), and with the involvement of reticular calcium, protein kinase C, extracellular signal-regulated kinase, p38 kinase, and nuclear factor-κB. NGF-induced enhancement of HIV-1 replication occurred during the late events of the HIV-1-replicative cycle, with a concomitant increase in viral transcription and production. In addition, NGF reduced the synthesis of the cellular HIV-1 restriction factor APOBEC3G and also overrode its interferon-γ-induced up-regulation, allowing the production of a well-fitted virus. Because NGF-TrKA signaling is a crucial event for macrophage survival, it is possible that NGF-induced HIV-1 replication plays a role in the maintenance of HIV-1 reservoirs. Our study may contribute to the understanding of the immunopathogenesis of HIV-1 infection and provide insights about approaches aimed at limiting viral replication in HIV-1 reservoirs.

PP2B/calcineurin-mediated desensitization of TRPV1 does not require AKAP150

Activation of protein kinases and phosphatases at the plasma membrane often initiates agonist-dependent signalling events. In sensory neurons, AKAP150 (A-kinase-anchoring protein 150) orientates PKA (protein kinase A), PKC (protein kinase C) and the Ca2+/calmodulin-dependent PP2B (protein phosphatase 2B, also known as calcineurin) towards membrane-associated substrates. Recent evidence indicates that AKAP150-anchored PKA and PKC phosphorylate and sensitize the TRPV1 (transient receptor potential subfamily V type 1 channel, also known as the capsaicin receptor). In the present study, we explore the hypothesis that an AKAP150-associated pool of PP2B catalyses the dephosphorylation and desensitization of TRPV1. Biochemical, electrophysiological and cell-based experiments indicate that PP2B associates with AKAP150 and TRPV1 in cultured TG (trigeminal ganglia) neurons. Gene silencing of AKAP150 reduces basal phosphorylation of TRPV1. However, functional studies in neurons isolated from AKAP150-/- mice indicate that the anchoring protein is not required for pharmacological desensitization of TRPV1. Behavioural analysis of AKAP150-/- mice further support this notion, demonstrating that agonist-stimulated desensitization of TRPV1 is sensitive to PP2B inhibition and does not rely on AKAP150. These findings allow us to conclude that pharmacological desensitization of TRPV1 by PP2B may involve additional regulatory components.

The integration of a Stat3 specific peptide aptamer into the thioredoxin scaffold protein strongly enhances its inhibitory potency

We are characterizing peptides which are able to interact with functional domains of oncoproteins and thus inhibit their activity. The yeast two-hybrid system was used to derive a peptide sequence which specifically interacts with the dimerization domain of the transcription factor Stat3. The activated form of Stat3 is required for the survival of many transformed cells and Stat3 inhibition can cause tumor cell death. The genetic selection of specific peptide sequences from random peptide libraries requires the integration into a scaffold protein and the expression in yeast cells. The scaffold protein, a variant of the human thioredoxin protein, has previously been optimized and also allows for effective bacterial expression of the recombinant protein and the cellular uptake of the purified, recombinant protein. We investigated the contributions of the scaffold protein to the inhibitory properties of rS3-PA. For this purpose we compared rS3-PA in which the ligand peptide is embedded within the thioredoxin scaffold protein with a minimal Stat3-interacting peptide sequence. sS3-P45 is a synthetic peptide of 45 amino acids in length and consists only of the Stat3-binding sequence of 20 amino acids, a protein transduction domain (PTD) and a Flag-tag. Both, the recombinant rS3-PA of 19.3 kDa and the synthetic sS3-P45 of 5.1 kDa, were taken up into the cytoplasm of cells by the PTD-mediated transduction process, inhibited Stat3 target gene expression and caused the death of Stat3-dependent tumor cells. Stat3-independent normal cells were unaffected. rS3-PA effectively inhibited Stat3 function at 2 μM, however, sS3-P45 was required at a concentration of 100 μM to exert the same effects. The more potent action of rS3-PA is most probably due to a conformational stabilization of the Stat3-interacting peptide in the context of the scaffold protein.

Proteins move! Protein dynamics and long-range allostery in cell signaling

An emerging point of view in protein chemistry is that proteins are not the static objects that are displayed in textbooks but are instead dynamic actors. Protein dynamics plays a fundamental role in many diseases, and spans a large hierarchy of timescales, from picoseconds to milliseconds or even longer. Nanoscale protein domain motion on length scales comparable to protein dimensions is key to understanding how signals are relayed through multiple protein-protein interactions. A canonical example is how the scaffolding proteins NHERF1 and ezrin work in coordination to assemble crucial membrane complexes. As membrane-cytoskeleton scaffolding proteins, these provide excellent prototypes for understanding how regulatory signals are relayed through protein-protein interactions between the membrane and the cytoskeleton. Here, we review recent progress in understanding the structure and dynamics of the interaction. We describe recent novel applications of neutron spin echo spectroscopy to reveal the dynamic propagation of allosteric signals by nanoscale protein motion, and present a guide to the future study of dynamics and its application to the cure of disease.

The assembly of a GTPase-kinase signalling complex by a bacterial catalytic scaffold

The fidelity and specificity of information flow within a cell is controlled by scaffolding proteins that assemble and link enzymes into signalling circuits. These circuits can be inhibited by bacterial effector proteins that post-translationally modify individual pathway components. However, there is emerging evidence that pathogens directly organize higher-order signalling networks through enzyme scaffolding, and the identity of the effectors and their mechanisms of action are poorly understood. Here we identify the enterohaemorrhagic Escherichia coli O157:H7 type III effector EspG as a regulator of endomembrane trafficking using a functional screen, and report ADP-ribosylation factor (ARF) GTPases and p21-activated kinases (PAKs) as its relevant host substrates. The 2.5 Å crystal structure of EspG in complex with ARF6 shows how EspG blocks GTPase-activating-protein-assisted GTP hydrolysis, revealing a potent mechanism of GTPase signalling inhibition at organelle membranes. In addition, the 2.8 Å crystal structure of EspG in complex with the autoinhibitory Iα3-helix of PAK2 defines a previously unknown catalytic site in EspG and provides an allosteric mechanism of kinase activation by a bacterial effector. Unexpectedly, ARF and PAKs are organized on adjacent surfaces of EspG, indicating its role as a 'catalytic scaffold' that effectively reprograms cellular events through the functional assembly of GTPase-kinase signalling complex.

The large isoforms of A-kinase anchoring protein 18 mediate the phosphorylation of inhibitor-1 by protein kinase A and the inhibition of protein phosphatase 1 activity

Inhibitor-1 (I-1) is phosphorylated on threonine residue 35 (Thr35) by the cAMP-dependent protein kinase (PKA), inducing the potent inhibition of the serine-threonine-specific protein phosphatase 1 (PP1). We now report that the formation of a signaling complex containing PKA and I-1 by the A-kinase anchoring protein 18 (AKAP18) facilitates this regulation in cells. AKAP18 directly bound I-1, and AKAP18/I-1 complexes were isolated from both rat heart extract and transfected heterologous cells. It is noteworthy that prevention of PKA binding to the AKAP18 scaffold decreased I-1 phosphorylation by 48% in cells. Moreover, the I-1 target PP1 was also associated with AKAP18 complexes. The cAMP-mediated inhibition of phosphatase activity was contingent on PKA binding to the scaffold. These observations reveal an additional level of complexity in PP1 regulation because of its association with AKAP18 multimolecular signaling complexes and suggest that targeting of AKAP18 complexes may be an alternative method to alter phosphatase activity and modulate specific substrate dephosphorylation.

Comparison of strong cation exchange and SDS-PAGE fractionation for analysis of multiprotein complexes

Affinity purification of protein complexes followed by identification using liquid chromatography/mass spectrometry (LC-MS/MS) is a robust method to study the fundamental process of protein interaction. Although affinity isolation reduces the complexity of the sample, fractionation prior to LC-MS/MS analysis is still necessary to maximize protein coverage. In this study, we compared the protein coverage obtained via LC-MS/MS analysis of protein complexes prefractionated using two commonly employed methods, SDS-PAGE and strong cation exchange chromatography (SCX). The two complexes analyzed focused on the nuclear proteins Bmi-1 and GATA3 that were expressed within the cells at low and high levels, respectively. Prefractionation of the complexes at the peptide level using SCX consistently resulted in the identification of approximately 3-fold more proteins compared to separation at the protein level using SDS-PAGE. The increase in the number of identified proteins was especially pronounced for the Bmi-1 complex, where the target protein was expressed at a low level. The data show that prefractionation of affinity isolated protein complexes using SCX prior to LC-MS/MS analysis significantly increases the number of identified proteins and individual protein coverage, particularly for target proteins expressed at low levels.

Stress induced subcellular distribution of ALG-2, RBM22 and hSlu7

ALG-2 is a highly conserved calcium binding protein in the cytoplasm which belongs to the family of penta-EF hand proteins. Recently, we showed that ALG-2 is interacting with RBM22, a highly conserved spliceosomal nuclear protein (Montaville et al. Biochim. Biophys. Acta 1763, 1335, 2006; Krebs, Biochim. Biophys. Acta 1793, 979, 2009). In NIH 3T3 cells expressing both proteins a significant amount of ALG-2mRFP is translocated to the nucleus due to the interaction with RBM22-EGFP. hSlu7, another spliceosomal nuclear protein, known to interact with RBM22 in yeast, has been shown to translocate to the cytoplasm under cellular stress conditions. Here we provide evidence that the 2 spliceosomal proteins differ significantly in their subcellular distributions under stress conditions, and that RBM22 enhances the cytoplasmic translocation of hSlu7 under stress, especially a stress induced by thapsigargin. On the other hand, in NIH 3T3 cells expressing RBM22-EGFP and ALG-2-mRFP, ALG-2 remains translocated into the nucleus under both stress conditions, i.e. heat shock or treatment with thapsigargin. We could further demonstrate that these stress conditions had a different influence on the splicing pattern of XBP-1, a marker for the unfolded protein response indicating that ER stress may play a role in stress-induced translocation of spliceosomal proteins. The article is part of a Special Issue entitled: 11th European Symposium on Calcium.

Alpha-dystrobrevin-1 recruits alpha-catulin to the alpha1D-adrenergic receptor/dystrophin-associated protein complex signalosome

α(1D)-Adrenergic receptors (ARs) are key regulators of cardiovascular system function that increase blood pressure and promote vascular remodeling. Unfortunately, little information exists about the signaling pathways used by this important G protein-coupled receptor (GPCR). We recently discovered that α(1D)-ARs form a "signalosome" with multiple members of the dystrophin-associated protein complex (DAPC) to become functionally expressed at the plasma membrane and bind ligands. However, the molecular mechanism by which the DAPC imparts functionality to the α(1D)-AR signalosome remains a mystery. To test the hypothesis that previously unidentified molecules are recruited to the α(1D)-AR signalosome, we performed an extensive proteomic analysis on each member of the DAPC. Bioinformatic analysis of our proteomic data sets detected a common interacting protein of relatively unknown function, α-catulin. Coimmunoprecipitation and blot overlay assays indicate that α-catulin is directly recruited to the α(1D)-AR signalosome by the C-terminal domain of α-dystrobrevin-1 and not the closely related splice variant α-dystrobrevin-2. Proteomic and biochemical analysis revealed that α-catulin supersensitizes α(1D)-AR functional responses by recruiting effector molecules to the signalosome. Taken together, our study implicates α-catulin as a unique regulator of GPCR signaling and represents a unique expansion of the intricate and continually evolving array of GPCR signaling networks.

Emerging roles for the FSH receptor adapter protein APPL1 and overlap of a putative 14-3-3τ interaction domain with a canonical G-protein interaction site

The interaction of cytoplasmic proteins with intracellular domains of membrane receptors can occur at several opportunities, including: during biosynthesis, while in membrane residency and during internalization and recycling following ligand binding. Since the initial discovery that it interacts with the FSH receptor (FSHR) together with additional members of a potential signaling complex, APPL1 has been shown to interact with a variety of membrane receptors. Recent subcellular localizations of APPL1 place it in dynamic and varied venues in the cell, including at the cell membrane, the nucleus and the early endosomes. Another adapter protein family the 14-3-3 proteins, are largely recognized as binding to phosphorylation sites but recent work demonstrated that in the case of FSHR, the 14-3-3 site overlaps with the canonical G-protein binding site. G-proteins appear to sample the environment and exchange between the membrane and intracellular locales and this binding could be mediated by or modulated by receptor interactions at the 14-3-3 binding site. Observations that multiple proteins can interact with cytoplasmic domains of GPCRs leads to the inescapable conclusion that either the interactions occur via orderly replacement or exchange, or that receptors are simultaneously occupied by a variety of adapters and effectors or even that oligomers of dimeric GPCRs provide for platforms that can simultaneously interact with effectors and adaptors.

AKAP-Lbc enhances cyclic AMP control of the ERK1/2 cascade

Mitogen-activated protein kinase (MAPK) cascades propagate a variety of cellular activities. Processive relay of signals through RAF-MEK-ERK modulates cell growth and proliferation. Signalling through this ERK cascade is frequently amplified in cancers, and drugs such as sorafenib (which is prescribed to treat renal and hepatic carcinomas) and PLX4720 (which targets melanomas) inhibit RAF kinases. Natural factors that influence ERK1/2 signalling include the second messenger cyclic AMP. However, the mechanisms underlying this cascade have been difficult to elucidate. We demonstrate that the A-kinase-anchoring protein AKAP-Lbc and the scaffolding protein kinase suppressor of Ras (KSR-1) form the core of a signalling network that efficiently relay signals from RAF, through MEK, and on to ERK1/2. AKAP-Lbc functions as an enhancer of ERK signalling by securing RAF in the vicinity of MEK1 and synchronizing protein kinase A (PKA)-mediated phosphorylation of Ser 838 on KSR-1. This offers mechanistic insight into cAMP-responsive control of ERK signalling events.

Mechanism of signal transduction of the LOV2-Jα photosensor from Avena sativa

Fusion proteins containing blue-light-activable protein domains possess great potential as molecular switches in cell signalling. This has recently been impressively demonstrated by connecting the light oxygen voltage LOV2-Jα-protein domain of A. sativa (AsLOV2-Jα) with the Rac1-GTPase, responsible for regulating the morphology and motility of metazoan cells. However, a target-oriented development of fusion proteins in conjunction with this photosensor is still very challenging, because a detailed understanding of its signal transduction pathway on a molecular level is still lacking. Here, we show through molecular dynamics simulation that, after formation of the cysteinyl-flavin mononucleotide (FMN) adduct, the signalling pathway begins with a rotational reorientation of the residue glutamine 1029 adjacent to the FMN chromophore, transmitting stress through the Iβ strand towards the LOV2-Jα interface. This then results in the breakage of two H-bonds, namely, glutamic acid 1034-Gln995 and aspartic acid (Asp) 1056-Gln1013, at opposite sides of the interface between the Jα helix and the LOV2 domain, ultimately leading to a disruption of Jα helix from the LOV2 core.

The role of oligomerization and cooperative regulation in protein function: the case of tryptophan synthase

The oligomerization/co-localization of protein complexes and their cooperative regulation in protein function is a key feature in many biological systems. The synergistic regulation in different subunits often enhances the functional properties of the multi-enzyme complex. The present study used molecular dynamics and Brownian dynamics simulations to study the effects of allostery, oligomerization and intermediate channeling on enhancing the protein function of tryptophan synthase (TRPS). TRPS uses a set of α/β–dimeric units to catalyze the last two steps of L-tryptophan biosynthesis, and the rate is remarkably slower in the isolated monomers. Our work shows that without their binding partner, the isolated monomers are stable and more rigid. The substrates can form fairly stable interactions with the protein in both forms when the protein reaches the final ligand–bound conformations. Our simulations also revealed that the α/β–dimeric unit stabilizes the substrate–protein conformation in the ligand binding process, which lowers the conformation transition barrier and helps the protein conformations shift from an open/inactive form to a closed/active form. Brownian dynamics simulations with a coarse-grained model illustrate how protein conformations affect substrate channeling. The results highlight the complex roles of protein oligomerization and the fine balance between rigidity and dynamics in protein function.

Conformational changes of enzymes are often related to regulating and creating an optimal environment for efficient chemistry. An increasing number of evidences also indicate that oligomerization/co-localization of proteins contributes to the efficiency of metabolic pathways. Although static structures have been available for many multi-enzyme complexes, their efficiency is also governed by the synergistic regulation between the multi-units. Our study applies molecular dynamics and Brownian dynamics simulations to the model system, the tryptophan synthase complex. The multi-enzyme complex is a bienzyme nanomachine and its catalytic activity is intimately related to the allosteric signaling and the metabolite transfer between its α– and β–subunits connected by a 25-Å long channel. Our studies suggest that the binding partner is crucial for the ligand binding processes. Although the isolated monomers are stable in the ligand–free state and can form stable interaction if the substrate is in the final bound conformation, it has higher energy barrier when ligand binds to the active site. We also show that the channel does not always exist, but it may be blocked before the enzyme reaches its final bound conformation. The results highlight the importance of forming protein complexes and the cooperative changes during different states.

Chemotherapy and signaling: How can targeted therapies supercharge cytotoxic agents?

In recent years, oncologists have begun to conclude that chemotherapy has reached a plateau of efficacy as a primary treatment modality, even if toxicity can be effectively controlled. Emerging specific inhibitors of signaling and metabolic pathways (i.e., targeted agents) contrast with traditional chemotherapy drugs in that the latter primarily interfere with the DNA biosynthesis and the cell replication machinery. In an attempt to improve on the efficacy, combination of targeted drugs with conventional chemotherapeutics has become a routine way of testing multiple new agents in early phase clinical trials. This review discusses the recent advances including integrative systematic biology and RNAi approaches to counteract the chemotherapy resistance and to buttress the selectivity, efficacy and personalization of anti-cancer drug therapy.

Chk-ing in and Chk-ing out: kinase compartmentalization comes to checkpoint control

Checkpoints are the sentinels of cell-cycle progression. In this issue of Molecular Cell, Yaffe and colleagues (Reinhardt et al., 2010) show that spatial and temporal resolution of Chk1 and MK2, checkpoint kinases with identical substrate specificity, are necessary to signal different aspects of DNA damage signaling.

Cell cycle arrest induced by the bacterial adenylate cyclase toxins from Bacillus anthracis and Bordetella pertussis

Bacillus anthracis oedema toxin (ET) and Bordetella pertussis adenylate cyclase toxin (ACT) enter host cells and produce cAMP. To understand the cellular consequences, we exposed J774 cells to these toxins at ng ml(-1) (pM) concentrations, then followed cell number and changes in cell signalling pathways. Under these conditions, both toxins produce a concentration-dependent inhibition of cell proliferation without cytotoxicity. ET and ACT increase the proportion of cells in G(1) /G(0) and reduce S phase, such that a single addition of ET or ACT inhibits cell division for 3-6 days. Treatment with ET or ACT produces striking changes in proteins controlling cell cycle, including virtual elimination of phosphorylated ERK 1/2 and Cyclin D1 and increases in phospho-CREB and p27(Kip1) . Importantly, PD98059, a MEK inhibitor, elicits a comparable reduction in Cyclin D1 to that produced by the toxins and blocks proliferation. These data show that non-lethal concentrations of ET and ACT impose a prolonged block on the proliferation of J774 cells by impairment of the progression from G(1) /G(0) to S phase in a process involving cAMP-mediated increases in phospho-CREB and p27(Kip1) and reductions in phospho-ERK 1/2 and Cyclin D1. This phenomenon represents a new mechanism by which these toxins affect host cells.

Organizing signal transduction through A-kinase anchoring proteins (AKAPs)

A fundamental role for protein-protein interactions in the organization of signal transduction pathways is evident. Anchoring, scaffolding and adapter proteins function to enhance the precision and directionality of these signaling events by bringing enzymes together. The cAMP signaling pathway is organized by A-kinase anchoring proteins. This family of proteins assembles enzyme complexes containing the cAMP-dependent protein kinase, phosphoprotein phosphatases, phosphodiesterases and other signaling effectors to optimize cellular responses to cAMP and other second messengers. Selected A-kinase anchoring protein signaling complexes are highlighted in this minireview.

Azirine ligation: fast and selective protein conjugation via photoinduced azirine-alkene cycloaddition

We report a new bioorthogonal ligation reaction between p-nitrodiphenylazirine and dimethyl fumarate. This photoinduced azirine-alkene cycloaddition provides a rapid (~2 min) and highly selective route to protein conjugation at neutral pH and room temperature in biological medium.

TAZ interacts with zonula occludens-1 and -2 proteins in a PDZ-1 dependent manner

The transcriptional coactivator TAZ recognizes L/PPxY motifs in transcription factors like Runx1/2 through its WW domain. We show that the first PDZ domain of zona occludens-1 (ZO-1) and 2 (ZO-2) interacts with the carboxy-terminal PDZ binding motif of TAZ. Deletion of this motif abrogates binding. ZO-2 colocalizes with TAZ in the nucleus of MDCK cells and ZO-2 expression alters TAZ localization in human embryonic kidney cells. Luciferase assays demonstrate ZO-2 inhibition of TAZ-mediated transactivation. We propose that zonula occludens is a negative regulator of TAZ and suggest that selected tight junction proteins control nuclear translocation and activity of TAZ.

Topographic studies of the GroEL-GroES chaperonin complex by chemical cross-linking using diformyl ethynylbenzene: the power of high resolution electron transfer dissociation for determination of both peptide sequences and their attachment sites

Many essential cellular processes depend upon the self-assembly of stable multiprotein entities. The architectures of the vast majority of these protein machines remain unknown because these structures are difficult to obtain by biophysical techniques alone. However, recent progress in defining the architecture of protein complexes has resulted from integrating information from all available biochemical and biophysical sources to generate computational models. Chemical cross-linking is a technique that holds exceptional promise toward achieving this goal by providing distance constraints that reflect the topography of protein complexes. Combined with the available structural data, these constraints can yield three-dimensional models of higher order molecular machines. However, thus far the utility of cross-linking has been thwarted by insufficient yields of cross-linked products and tandem mass spectrometry methods that are unable to unambiguously establish the identity of the covalently labeled peptides and their sites of modification. We report the cross-linking of amino moieties by 1,3-diformyl-5-ethynylbenzene (DEB) with analysis by high resolution electron transfer dissociation. This new reagent coupled with this new energy deposition technique addresses these obstacles by generating cross-linked peptides containing two additional sites of protonation relative to conventional cross-linking reagents. In addition to excellent coverage of sequence ions by electron transfer dissociation, DEB cross-linking produces gas-phase precursor ions in the 4+, 5+, or 6+ charge states that are readily segregated from unmodified and dead-end modified peptides using charge-dependent precursor selection of only quadruply and higher charge state ions. Furthermore, electron transfer induces dissociation of the DEB-peptide bonds to yield diagnostic ion signals that reveal the "molecular ions" of the unmodified peptides. We demonstrate the power of this strategy by cross-linking analysis of the 21-protein, ADP-bound GroEL-GroES chaperonin complex. Twenty-five unique sites of cross-linking were determined.

Dopamine reveals neural circuit mechanisms of fly memory

A goal of memory research is to understand how changing the weight of specific synapses in neural circuits in the brain leads to an appropriate learned behavioral response. Finding the relevant synapses should allow investigators to probe the underlying physiological and molecular operations that encode memories and permit their retrieval. In this review I discuss recent work in Drosophila that implicates specific subsets of dopaminergic (DA) neurons in aversive reinforcement and appetitive motivation. The zonal architecture of these DA neurons is likely to reveal the functional organization of aversive and appetitive memory in the mushroom bodies. Combinations of fly DA neurons might code negative and positive value, consistent with a motivational systems role as proposed in mammals.

RuleMonkey: software for stochastic simulation of rule-based models

Background

The system-level dynamics of many molecular interactions, particularly protein-protein interactions, can be conveniently represented using reaction rules, which can be specified using model-specification languages, such as the BioNetGen language (BNGL). A set of rules implicitly defines a (bio)chemical reaction network. The reaction network implied by a set of rules is often very large, and as a result, generation of the network implied by rules tends to be computationally expensive. Moreover, the cost of many commonly used methods for simulating network dynamics is a function of network size. Together these factors have limited application of the rule-based modeling approach. Recently, several methods for simulating rule-based models have been developed that avoid the expensive step of network generation. The cost of these "network-free" simulation methods is independent of the number of reactions implied by rules. Software implementing such methods is now needed for the simulation and analysis of rule-based models of biochemical systems.

Results

Here, we present a software tool called RuleMonkey, which implements a network-free method for simulation of rule-based models that is similar to Gillespie's method. The method is suitable for rule-based models that can be encoded in BNGL, including models with rules that have global application conditions, such as rules for intramolecular association reactions. In addition, the method is rejection free, unlike other network-free methods that introduce null events, i.e., steps in the simulation procedure that do not change the state of the reaction system being simulated. We verify that RuleMonkey produces correct simulation results, and we compare its performance against DYNSTOC, another BNGL-compliant tool for network-free simulation of rule-based models. We also compare RuleMonkey against problem-specific codes implementing network-free simulation methods.

Conclusions

RuleMonkey enables the simulation of rule-based models for which the underlying reaction networks are large. It is typically faster than DYNSTOC for benchmark problems that we have examined. RuleMonkey is freely available as a stand-alone application http://public.tgen.org/rulemonkey. It is also available as a simulation engine within GetBonNie, a web-based environment for building, analyzing and sharing rule-based models.