Outline
Comptes Rendus

Towards a systems biology approach of G protein-coupled receptor signalling: Challenges and expectations
Comptes Rendus. Biologies, Integrative Biology and Modelling in Agronomical Sciences, Volume 332 (2009) no. 11, pp. 947-957.

Abstracts

G protein-coupled receptors (GPCRs) control all the main physiological functions and are targeted by more than 50% of therapeutics. Our perception of GPCRs signalling has grown increasingly complex since it is now accepted that they activate large signalling networks which are integrating the information fluxes into appropriate biological responses. These concepts lead the way to the development of pathway-selective agonists (or antagonists) with fewer side effects. Systems biology approaches focused on GPCR-mediated signalling would help dealing with the huge complexity of these mechanisms therefore speeding-up the discovery of new drug classes. In this review, we present the various technical and conceptual possibilities allowing a systems approach of GPCR-mediated signalling. The main remaining limitations are also discussed.

Les récepteurs couplés aux protéines G (RCPGs) contrôlent toutes les grandes fonctions physiologiques et sont la cible de plus de 50% des médicaments. Notre perception de la signalisation des RCPGs s'est considérablement complexifiée puisqu'il est maintenant admis qu'ils activent de larges réseaux de signalisation capables d'intégrer les flux d'information en des réponses biologiques appropriées. Ces concepts ouvrent la voie au développement d'agonistes (ou d'antagonistes) sélectifs de voies de signalisation qui présenteraient moins d'effets indésirables. Une démarche de biologie systémique appliquée à la signalisation des RCPGs aiderait à appréhender la complexité des mécanismes de signalisation et accélérerait ainsi l'émergence de nouvelles classes de médicaments. Dans cette revue, nous présentons les différentes possibilités techniques et conceptuelles permettant la mise en place d'une approche systémique de la signalisation des RCPGs. Les principaux défis auxquels nous restons confrontés sont également discutés.

Metadata
Published online:
DOI: 10.1016/j.crvi.2009.09.002
Keywords: GPCR, Signalling, Pharmacology, Systems biology, Bioinformatics, Modelling
Mots-clés : RCPG, Signalisation, Pharmacologie, Biologie systémique, Bioinformatique, Modélisation

Domitille Heitzler 1, 2, 3, 4; Pascale Crépieux 1, 2, 3, 4; Anne Poupon 1, 2, 3, 4; Frédérique Clément 5; François Fages 6; Eric Reiter 1, 2, 3, 4

1 BIOS group, INRA, UMR85, Unité physiologie de la reproduction et des comportements, 37380 Nouzilly, France
2 CNRS, UMR6175, 37380 Nouzilly, France
3 Université François-Rabelais, 37041 Tours, France
4 Haras Nationaux, 37380 Nouzilly, France
5 Sisyphe team, INRIA Paris-Rocquencourt, 78153 Le Chesnay, France
6 Contraintes team, INRIA Paris-Rocquencourt, 78153 Le Chesnay, France
@article{CRBIOL_2009__332_11_947_0,
     author = {Domitille Heitzler and Pascale Cr\'epieux and Anne Poupon and Fr\'ed\'erique Cl\'ement and Fran\c{c}ois Fages and Eric Reiter},
     title = {Towards a systems biology approach of {G} protein-coupled receptor signalling: {Challenges} and expectations},
     journal = {Comptes Rendus. Biologies},
     pages = {947--957},
     publisher = {Elsevier},
     volume = {332},
     number = {11},
     year = {2009},
     doi = {10.1016/j.crvi.2009.09.002},
     language = {en},
}
TY  - JOUR
AU  - Domitille Heitzler
AU  - Pascale Crépieux
AU  - Anne Poupon
AU  - Frédérique Clément
AU  - François Fages
AU  - Eric Reiter
TI  - Towards a systems biology approach of G protein-coupled receptor signalling: Challenges and expectations
JO  - Comptes Rendus. Biologies
PY  - 2009
SP  - 947
EP  - 957
VL  - 332
IS  - 11
PB  - Elsevier
DO  - 10.1016/j.crvi.2009.09.002
LA  - en
ID  - CRBIOL_2009__332_11_947_0
ER  - 
%0 Journal Article
%A Domitille Heitzler
%A Pascale Crépieux
%A Anne Poupon
%A Frédérique Clément
%A François Fages
%A Eric Reiter
%T Towards a systems biology approach of G protein-coupled receptor signalling: Challenges and expectations
%J Comptes Rendus. Biologies
%D 2009
%P 947-957
%V 332
%N 11
%I Elsevier
%R 10.1016/j.crvi.2009.09.002
%G en
%F CRBIOL_2009__332_11_947_0
Domitille Heitzler; Pascale Crépieux; Anne Poupon; Frédérique Clément; François Fages; Eric Reiter. Towards a systems biology approach of G protein-coupled receptor signalling: Challenges and expectations. Comptes Rendus. Biologies, Integrative Biology and Modelling in Agronomical Sciences, Volume 332 (2009) no. 11, pp. 947-957. doi : 10.1016/j.crvi.2009.09.002. https://comptes-rendus.academie-sciences.fr/biologies/articles/10.1016/j.crvi.2009.09.002/

Original version of the full text

1 Introduction

An impressive amount of detailed information has been gathered over the past decades on how external stimuli activate plasma membrane receptors, how they translate to the activation of linear downstream signalling cascades and eventually affect cell fate. Recently, the advent of highly sensitive proteomic methods has produced maps of protein interactions and led to the reconstruction of biochemical networks [1]. As a consequence, it is now widely accepted that signalling pathways are organized as coordinated communication networks in which multi-protein complexes process and integrate the signal fluxes. Now, the challenge in cell signalling is to understand the behaviour of these intertwined communication networks in order to decipher the cellular language [2]. G protein-coupled receptors (GPCR) represent the largest class of membrane receptors. They are capable of binding a wide diversity of molecules that regulate most physiological processes and are involved in a plethora of diseases. Noteworthy, GPCRs have long been preferential targets of therapeutic research and development and they currently account for up to 50% of marketed drugs [3].

2 The growing complexity of GPCR-induced signalling

Classically, upon ligand binding, GPCRs undergo a conformational change that leads to heterotrimeric G protein recruitment and activation, followed by the generation of diffusible second messengers such as cAMP (cyclic Adenosine Mono-Phosphate), calcium or phosphoinositides. However, it is increasingly recognized that GPCRs trigger multiple signalling pathways which lead to the formation of signalling networks [4] (Fig. 1). For instance, some GPCRs have the ability to couple to multiple G protein subtypes [5] and many GPCRs directly interact with non-G protein signalling effectors through specific protein–protein interaction domains [6]. But quite remarkably, outside of heterotrimeric G proteins, only two protein families are able to specifically interact with the majority of GPCRs in their activated conformation: G protein-coupled receptor kinases (GRKs) and β-arrestins [7]. Historically, GRKs and β-arrestins have been associated with the desensitization and internalization/recycling of most GPCRs [8]. However, recently, GPCRs have also been demonstrated to elicit signals, independently of heterotrimeric G protein coupling, through interaction with β-arrestins 1 and 2. Indeed, β-arrestins have been shown to act as multifunctional scaffolds and activators for a growing number of signalling proteins including ERK, p38, JNK, I-κB, Akt and RhoA [7,9–13]. Moreover, a recent proteomic study has reported as many as 337 protein interactions involving β-arrestins [14], strongly suggesting that they play a central role in the ability of GPCRs to activate very complex signalling networks. In addition, GRKs have also been reported to elicit signalling responses on their own right through protein/protein interactions [13]. Indeed, GRKs interact with a variety of proteins involved in signalling and trafficking such as Gαq, Gβγ, PI3Kγ, clathrin, GIT (G protein-coupled receptor kinase-interacting protein) and caveolin [15]. Phosphorylation of Raf kinase inhibitor protein (RKIP) by PKC displaces it from Raf and increases its association with GRK2 [16]. In addition, the physical interaction between GRK2 and Akt leads to the inhibition of Akt activity [17]. Finally, GRK2 and MEK1 have been found in the same multimolecular complex and this interaction is correlated with an inhibition of MEK activity [18].

Fig. 1

GPCR-induced signalling networks. Schematic representation of the different transduction mechanisms leading to the activation of signalling pathways upon activation of a GPCR and how they define a complex signalling network. Cross-activations and cross-inhibitions between the G protein-dependent and -independent sub-networks are depicted as green and red arrows respectively. Gs, Gq, Gi and G12: different alpha subunits of heterotrimeric G protein; Gβγ: G beta and gamma subunits of heterotrimeric G proteins; other GIPs: other GPCR interacting proteins; β-arrs: β-arrestin 1 and 2. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.)

Adding to this complexity is the fact that GPCR-induced signals can be spatially and temporally encoded. Signalling networks actively modulate the transmitted signals: negative feedback allows pathways to adapt or desensitize to persistent stimuli whereas cross inhibition is used to avoid crosstalk between pathways [19,20]. In addition to transmit qualitative information (e.g. the presence or absence of a stimulus), signalling pathways must also convey quantitative information about the strength of the stimulus (i.e. ligand concentration). It has been recently shown that signalling pathways can take advantage of their non-linear nature to convert stimulus intensity into signal duration [21]. Modulation of signal duration increases the range of stimulus concentrations for which dose-dependent responses are possible as dose-dependent responses are still possible after apparent saturation of the receptors. Another well documented example of spatial and temporal encoding in GPCR-induced signalling pathways is the dual activation mechanism of ERK by G protein and β-arrestins [22–26]. G protein-mediated ERK activation is rapid, transient and translocates to the nucleus. In contrast, the ERK activated via β-arrestins are slower in onset (∼5–10 min to reach maximum), very persistent (t1/2>1 hour) and are sequestered in the cytosol. Such spatial and temporal differences in GPCR-induced signals substantially increase the complexity of signalling systems, hence their processing power. Collectively, these features highlight the importance of considering the dynamic properties of signalling pathways when characterizing their behaviour.

3 Pathway-selective ligands for GPCR: A new era in drug discovery?

Interestingly, this emerging conceptual framework opens new research avenues for the development of therapeutics [27]. There is increasing evidence that some GPCR targeting drugs can selectively modulate a subset of the signalling events triggered by the full agonist. These effects have been given various names including “stimulus-trafficking”, “biased agonism”, “collateral efficacy” or “functional selectivity” [28]. Moreover, several GPCR ligands have already been reported to selectively activate or inhibit β-arrestin signalling [29–31]. Consideration of these new concepts might lead to the development of therapeutics with more selective actions, hence less side-effects. For instance, one drug, carvedilol, a β adrenergic receptor antagonist, has proven particularly effective in the treatment of heart failure. Interestingly, of 16 clinically relevant β adrenergic receptor antagonists, carvedilol displays a unique ability to stimulate β-arrestin-mediated signalling while preventing receptor coupling to Gs [31].

4 What added value can systems biology provide?

GPCRs' complex signalling mechanisms probably lead to context-adapted cellular responses relying on emerging system-level properties that cannot be predicted from the individual components of the induced networks. Therefore, it would be of paramount interest to provide a conceptual framework for deciphering and possibly predicting how an extra-cellular signal that activates a GPCR translates into a given biological or pathological response. This “global” level of analysis of GPCRs' biology is in its infancy. Pioneer work carried out in yeast has recently shown the value of systems biology for elucidating complex signalling mechanisms triggered by GPCRs [32,33]. Aspects of cell signalling and the mechanisms (i.e. feedback and feed-forward regulations) that regulate pathway activity triggered by GPCRs have been studied in yeast, and nicely illustrate how mathematical modelling can be used to understand the logic of various pathway architectures. In mammals, a system-level grasp of GPCR-mediated signalling networks would be a significant asset to rationalize and speed-up the discovery of new “pathway-selective” drugs. Indeed, the rate of new drug discovery using standard approaches, based essentially on heterotrimeric G protein-dependent activities, such as second messenger accumulation, has been slowing down despite increased investments by the pharmaceutical industry.

It has been proposed that computational modelling offers a powerful tool for examining GPCR pathways [34]. Such models can be used to better understand hypothesized mechanisms, run virtual (in silico) experiments, interpret data, suggest new drug targets, motivate experiments, and offer new explanations for observed phenomena. In the remaining part of the present review, we will identify and discuss the different challenges that the scientists in the field will be confronted with in their efforts to establish a systems biology approach to GPCR signalling.

Understanding the whole cell as an integrated system has been a goal for almost a century [35]. As in other fields of cell biology, deciphering GPCR signalling by at least descriptive, at most analytical methods has occupied the last 4 decades. Consequently, an enormous amount of data is published every year on some GPCR signalling aspects. However, these studies are highly heterogeneous in terms of the nature of the GPCR studied, the signalling pathway studied, the cellular system used, etc. As a consequence, it has been tremendously difficult to transform this huge amount of information into general concepts. An important effort of standardization and data sharing between laboratories in the field is needed. In addition, publicly accessible databases gathering and distributing standardized raw data related to GPCR signalling, especially dynamic data, would undoubtedly stimulate the modelling of GPCR-induced signalling networks.

5 Challenges in high throughput generation of signalling data

More recently, feeding system-level analyses with relevant high quality biological data have become possible thanks to new experimental techniques that allow large scale accumulation of unbiased signalling data (Fig. 2). Speeding up the production of new data sets and enhancing their quality is not only essential to feed model-building but also to allow testing of key model findings. Signalling events are often propagated within the cell by post-translational modifications involving protein–protein interactions and enzymatic activities. Noteworthy, reversible protein phosphorylation is centrally involved in signal transmission within cells. The comprehensive and quantitative analysis of the protein phosphorylation patterns in different cellular backgrounds is therefore critical to reach a system level analysis of cell signalling. Lately, high-throughput has spread to molecular biology and biochemistry, giving access to most items of information necessary for the comprehension of organisms' behaviour [36]. In particular, breakthroughs have been achieved in the isolation of phosphorylated peptides from complex samples, as well as in their analysis by mass spectrometry coupled to computational methods. Using such approaches, thousands of phosphopeptides and phosphorylation sites can now be identified in a single sample [37–39]. However, despite their unparalleled analytical power, mass spectrometry-based approaches present a static snapshot of cellular events; they do not allow the acquisition of dynamical phosphorylation data at large-scale and with high throughput. This limitation remains a major hurdle towards the development of powerful systems biology approaches in the field of cellular signalling. Over the last decade, the possibility to analyze the proteome using protein microarray-based methods has emerged into proteomics research, diagnostics, and drug discovery. In particular, automated spotting of concentrated and complex protein extracts permits their analysis with phosphospecific antibodies. This method, referred to as Reverse-Phase Protein Array (RPPA), uses very small quantities of biological material, which allows a wide sample collection with a high number of different antibodies to be screened [40–43]. Since it has a very high throughput, detailed kinetic experiments can be systematically carried out and analyzed. Moreover, the simultaneous quantification of thousands of samples achieved with RPPA drastically reduces data heterogeneity and variability which traditionally hamper modelling. Therefore, RPPA potentially represents a very attractive approach to capture the subtle and highly dynamic nature of phosphorylation cascades out of very large sample collections.

Fig. 2

Schematic representation of three high-throughput methods particularly adapted to the model-building process in the context of GPCR-induced signalling. A. Reverse phase protein array (RPPA) relies on the ordered micro-arraying of cellular lysates. Specific signals are detected using a primary antibody-fluorescently-labelled secondary antibody revelation system. This technology is easy to automate and when coupled with the use of a repertoire of phospho-specific antibodies raised against signalling intermediates, it allows the analysis of signalling networks at an unprecedented throughput. B. Fluorescence resonance energy transfer (FRET) biosensors are genetically-encoded fusion proteins linking two fluorescent proteins with a peptide sensitive to either phosphorylation by a specific kinase or second messenger binding. Once altered by a signalling event, the linker undergoes a conformational change that modifies the relative positions of the two fluorescent proteins hence their resonance. Cells transfected with the sensor plasmid are analyzed using a microscope. The emission wavelength of the first fluorescent protein (generally CFP) corresponds to the excitation wavelength of the second fluorescent protein (generally YFP). Therefore, FRET can be quantified in living cells with exquisite time-resolution by using an excitation wavelength specific of the first fluorescent protein while recording the emission at the wavelength specific of the second fluorescent protein. C. Transfected cell array (TCA) consists of parallelized micro-spots containing a transfection reagent, a plasmid and/or a siRNA. Each spot is different and a very large number of them can be arrayed on one slide allowing large-scale screening. Cells are seeded directly on the slide so that they attach and grow on it while being specifically transfected by the genetic material underneath. The effect of each transfection can be analyzed for each spot individually based on morphology criteria, immunocytochemistry or reporter-gene activity.

Recently developed imaging approaches that use fluorescent sensors of signalling activities combine unmatched time and spatial resolution. Genetically encoded fluorescence resonance energy transfer (FRET)-based reporters have been used in living cells to monitor the spatiotemporal patterns of diffusible second messengers, kinase activities and GPCR activation [44–46]. Protein–protein interactions can also be analyzed in real time in living cells using either FRET or BRET (bioluminescence resonance energy transfer) [46–48]. When used in multi-well plate format, both FRET and BRET-based methods ensure the production of huge amounts of high content dynamic data which are very well suited to feed systems biology approaches.

If the ability to measure signalling outputs in a dynamic way is a key-point, so is the aptitude to control the inputs (i.e. agonist stimulation) both in time and concentrations. Chemical signalling in/between cells is crucial as it determines the cellular response, and is characterized by various time-scales ranging from a few milliseconds to several minutes [44]. However, with conventional experiments, mimicking in vivo conditions remains a challenge. The implementation of microfluidic devices might help suppress the limitations found with conventional approaches and allows multiplexed analysis with various stimulation patterns to be performed. Indeed, microfluidics offers the possibility of not only tightly controlling and modulating cell culture conditions but also of applying well-defined (in intensity, space and time) chemical stimulation (receptor ligands, inhibitors, etc.) [49,50]. This should notably facilitate an iterative dialog between computational modelling and “wet lab” experiments (i.e. prediction versus validation).

An important aspect when trying to decipher and model signalling networks is the ability to specifically apply perturbations and to measure the interactions of signalling pathways considered in the recent past as isolated entities (Fig. 3). In addition to classical approaches (e.g. kinase inhibitors, dominant negative constructs, etc.), interfering RNAs offer the unique opportunity to easily and specifically achieve gene knock-downs. Moreover, genome-wide siRNA screening are now available either in multi-well liquid phase or in transfected cell array format [51–53]. Large-scale siRNA screening may very well be an experimental breakthrough facilitating the edification of highly complex signalling networks.

Fig. 3

GPCR-induced signalling systems can be specifically perturbed using different strategies: A. chemical inhibitors and toxins can be used to specifically block some of the proximal transduction mechanisms as well as some of the nodes of the downstream network; B. interfering RNA provide the unique possibility to knock-down any single gene putatively involved in the signalling network; C. when available, the use of signalling pathway-selective “biased ligands” provides a very interesting way to imbalance the proximal transduction mechanism allowing to measure the consequences on downstream network; D. alternatively, signalling pathway-selective mutant GPCR can be engineered; E. an interesting and potentially powerful way to favour the much needed iterative interactions between model-based predictions and experimental validations consists in applying tightly controlled agonist stimulation patterns to cells. Of course, this approach implies the use of adapted microfluidic devices; F. finally, genetically-encoded modulators can also be serviceable for instance to modify reaction stoichiometry by overexpressing a protein. Dominant negative constructs can also block a specific reaction or interaction.

6 Challenges in bioinformatics

Knowledge and data management is instrumental to systems biology. Indeed, high-throughput approaches generate huge quantities of heterogeneous data that cannot be handled with classical labbooks, or even flat files. Therefore, computational methods, standards and tools must be developed and used to tackle this problem. Beyond the handling problem, the design, validation and refinement of high fidelity models require the use of all available experimental and non-experimental data. To reach this objective, data produced by experimental biology and bioinformatics have to be accessible, in compatible formats, and easy to correlate. Experimental data management systems, called Laboratory Information Management Systems (LIMS) should allow biologists to capture all the data relative to an experiment (conditions, protocol, results). So far, the availability of a LIMS to manage all the pieces of information required has been the major stumbling-block of the design of models in systems biology.

Once the data and knowledge are available and organized, the first modelling step consists in the construction of a detailed representation of the studied system. For intracellular signalling networks, the reaction/interaction graph is commonly used. Its construction is achieved manually, in a hypothesis-driven manner by the expert, using prior knowledge of the system and chosen experimental data. Sophisticated data analysis and visualization tools, such as clustering, help in this task. However, when the system is very large, choices have to be made by the expert, and the resulting model is biased, incomplete and often fails to reveal emerging features from the studied network. The development of methods to automatically infer influence graphs directly from the ensemble of data would allow the production of unbiased, thus potentially very innovative models from large data sets. It is certainly a promising research area but to date, only a few and limited attempts have been reported [54].

7 Challenges in mathematical modelling

Currently, relatively detailed mechanistic models of GPCR-induced signalling network can be formalized using the systems biology Markup Language (SBML) for representing the elementary interactions [55]. In this initial formalization step, signalling networks can be either directly written in SBML or drawn using CellDesigner which provides an intuitive SBML-based graphical modelling environment [56,57]. Of note is the Biochemical Abstract Machine BIOCHAM [58] which can be used to compute the influence graph between molecular species from the reaction/interaction graph [59], formalize the (observed) biological properties of the system by temporal logic formulae, and automatically verify their satisfiability by model-checking algorithms. From the influence graph, temporal logics and model-checking algorithms have proven useful to express biological properties of complex biochemical systems and automatically verify if they are appropriate. This approach has allowed the analysis of reachability and temporal logic properties of signalling systems under various conditions [60].

Considering the highly dynamic nature both in time and space of GPCR-induced signalling pathways, a dynamic modelling approach is interesting and should yield more predictive power than static approaches. The aim of a dynamic model is to reproduce time-course experiments and provide intrinsic concentrations for molecular species that are experimentally unreachable. There are several dynamical modelling approaches such as ordinary differential equations (ODEs, population view), Petri nets (discrete and independent mechanisms), and pi-calculus (stochastic approach) [61,62]. For all those methods, kinetics parameters (activation rates, probability of transition) are needed to simulate the changes in molecule concentrations over time. Frequently, too many parameters are unknown and the initial large reaction/interaction graphs have to be reduced. Stability analysis of the model (study of steady states, bifurcation diagrams, etc.) can help to restrict the parameter space to biologically relevant areas and reveal important differences in reaction speeds therefore providing an accurate way to reduce the model. Despite its usefulness, stability analysis remains too limited to ensure sufficient model reduction on its own. In addition, the qualitative properties of the influence graph can help develop a reduced dynamic model amenable to numerical simulations and parameter optimization with respect to quantitative data [63]. This model reduction approach remains however mostly empirical and even if the main graphical properties of the network (i.e. positive and negative circuits in the influence graph, reachability properties in the reaction graph) are preserved, the simplification of the structure of the network may suppress some delays and may limit the repertoire of dynamic behaviours permitted by the original model.

Even though stability analysis can lead to the determination of some unknown kinetic parameters in function of the others, the number of remaining unknown parameters is often large. Thus, another central difficulty concerns the non-linear optimization techniques that are needed to infer the unknown kinetic parameter values from the experimental data obtained under various conditions. Several techniques of parameter learning by data/property fitting can be used, such as gradient-based methods, Monte Carlo methods, and the Covariance Matrix Adaptation Evolution Strategy CMAES [64]. The latter two approaches minimize the error by repeated random sampling. Once the unknown parameters have been optimized, simulations can be performed and provide all component quantities over time. Simulations can be performed with different stimulation patterns or with in silico perturbations (modifications of total protein amounts, suppression of pathways, modulation of kinetics parameters, etc.). Thanks to these perturbations, the robustness of the system can be appreciated. Indeed, robustness is classically defined as the error between the perturbed and the initial simulations. Interestingly, the formalization of the expected dynamic properties of the system in temporal logic with numerical constraints also makes it possible to quantify the robustness of the system with respect to some important properties [65]. This opens the way for integrating robustness criteria in the process of building the model.

8 Necessity of an iterative dialog between experimentation and modelling

The value of a computational model also needs to be assessed by its ability to fit all the available experimental observations made in either control or perturbed conditions. Next, it can be used to make in silico predictions. Simulations present the decisive advantage of providing the temporal evolution of all the constitutive molecular species involved, including those that are experimentally unreachable. In addition, it is easy to systematically perturb the signalling system in silico, for instance by suppressing one molecule or one reaction; the predicted effects can then be verified experimentally. Interesting predictions often result from this process. Another way to evaluate the predictive power of a signalling network model is to modify the agonist input to the system. Typically, different agonist concentrations or patterns (i.e. positive or negative gradients, pulsatile mode, etc.) can be applied in silico.

After initial experimental data gathering, a model is designed, parameterized and simulated. Predictions are made based on this initial model, and validation experiments are undertaken. In this initial phase, it often happens that predictions are not experimentally validated. In this case, the model has to be modified accordingly, and submitted to experimental validation again. This iterative process between modelling and experimentation is referred to as model refinement. According to this workflow, any prediction leads to either hypothesis validation or model refinement (Fig. 4). Thanks to this virtuous circle, the decoding of the intimate functioning of signalling networks' should move forward. Understanding how signalling pathways encode and transmit quantitative information about the external environment not only deepens our understanding of these systems, but will also lead to major achievements both in agriculture and medicine: models will point to molecular targets to be optimized in normal animal physiology and will help restore the proper function of pathways that have become deregulated in disease.

Fig. 4

Proposed workflow to model signalling systems. Data are produced in the laboratory, using low and high throughput methods. Already existing relevant and carefully standardized data can be extracted from publicly accessible databases (note that, once exploited, locally produced data have the vocation to be deposited in such databases). Literature mining is also essential. Data and knowledge from these three sources are combined, organized, analyzed and made available for modelling using various bioinformatics tools integrated into a comprehensive LIMS (laboratory information management system). Influence graph is constructed and its qualitative analysis can lead to new hypotheses. In parallel, this initial model is reduced and parameters are optimized. The dynamical model can be analyzed and its predictive power evaluated. An iterative process then is engaged between experimentation and modelling which leads to novel findings and model refinement.

9 Conclusion

The acquisition of reliable data has long been the primary challenge in the field of GPCR signalling. Now, with the advent of high throughput technologies potentially able to generate an unprecedented amount of dynamic signalling data, and to take into account the growing complexity of GPCR signalling, the challenges are being re-centred towards more theoretical and conceptual aspects. State of the art bioinformatics and mathematics can already help better handle the complexity associated with the large signalling networks and with the tremendous development of systems biology, rapid improvements can be expected.

Acknowledgements

The authors thank members of the BIOS team for their advices and support. D.H. was funded by a fellowship from the INRA (ASC). This work was supported by the INRA AIP AgroBI and by the large scale project REGATE (REgulation of the GonAdoTropE axis).


References

[1] C. Preisinger; A. von Kriegsheim; D. Matallanas; W. Kolch Proteomics and phosphoproteomics for the mapping of cellular signalling networks, Proteomics, Volume 8 (2008), pp. 4402-4415

[2] D.B. Kell Theodor Bucher Lecture: Metabolomics, modelling and machine learning in systems biology – Towards an understanding of the languages of cells, FEBS J., Volume 273 (2006), pp. 873-894 (Delivered on 3 July 2005 at the 30th FEBS Congress and the 9th IUBMB Conference in Budapest)

[3] P. Ma; R. Zemmel Value of novelty?, Nat. Rev. Drug Discov., Volume 1 (2002), pp. 571-572

[4] L.M. Luttrell Transmembrane signaling by G protein-coupled receptors, Methods Mol. Biol., Volume 332 (2006), pp. 3-49

[5] R.J. Lefkowitz; K.L. Pierce; L.M. Luttrell Dancing with different partners: Protein kinase a phosphorylation of seven membrane-spanning receptors regulates their G protein-coupling specificity, Mol. Pharmacol., Volume 62 (2002), pp. 971-974

[6] J. Bockaert; A. Dumuis; L. Fagni; P. Marin GPCR-GIP networks: A first step in the discovery of new therapeutic drugs?, Curr. Opin. Drug Discov. Devel., Volume 7 (2004), pp. 649-657

[7] R.J. Lefkowitz; S.K. Shenoy Transduction of receptor signals by beta-arrestins, Science, Volume 308 (2005), pp. 512-517

[8] R.J. Lefkowitz G protein-coupled receptors. III. New roles for receptor kinases and beta-arrestins in receptor signaling and desensitization, J. Biol. Chem., Volume 273 (1998), pp. 18677-18680

[9] D.S. Witherow; T.R. Garrison; W.E. Miller; R.J. Lefkowitz Beta-arrestin inhibits NF-kappaB activity by means of its interaction with the NF-kappaB inhibitor IkappaBalpha, Proc. Natl. Acad. Sci. USA, Volume 101 (2004), pp. 8603-8607

[10] W.G. Barnes; E. Reiter; J.D. Violin; X.R. Ren; G. Milligan; R.J. Lefkowitz {beta}-Arrestin 1 and G{alpha}q/11 coordinately activate RhoA and stress fiber formation following receptor stimulation, J. Biol. Chem., Volume 280 (2005), pp. 8041-8050

[11] J.M. Beaulieu; T.D. Sotnikova; S. Marion; R.J. Lefkowitz; R.R. Gainetdinov; M.G. Caron An Akt/beta-arrestin 2/PP2A signaling complex mediates dopaminergic neurotransmission and behavior, Cell, Volume 122 (2005), pp. 261-273

[12] D.L. Hunton; W.G. Barnes; J. Kim; X.R. Ren; J.D. Violin; E. Reiter; G. Milligan; D.D. Patel; R.J. Lefkowitz Beta-arrestin 2-dependent angiotensin II type 1A receptor-mediated pathway of chemotaxis, Mol. Pharmacol., Volume 67 (2005), pp. 1229-1236

[13] E. Reiter; R.J. Lefkowitz GRKs and beta-arrestins: Roles in receptor silencing, trafficking and signaling, Trends Endocrinol. Metab., Volume 17 (2006), pp. 159-165

[14] K. Xiao; D.B. McClatchy; A.K. Shukla; Y. Zhao; M. Chen; S.K. Shenoy; J.R. Yates; R.J. Lefkowitz Functional specialization of beta-arrestin interactions revealed by proteomic analysis, Proc. Natl. Acad. Sci. USA, Volume 104 (2007), pp. 12011-12016

[15] P. Penela; C. Ribas; F. Mayor Mechanisms of regulation of the expression and function of G protein-coupled receptor kinases, Cell. Signal., Volume 15 (2003), pp. 973-981

[16] K. Lorenz; M.J. Lohse; U. Quitterer Protein kinase C switches the Raf kinase inhibitor from Raf-1 to GRK-2, Nature, Volume 426 (2003), pp. 574-579

[17] S. Liu; R.T. Premont; C.D. Kontos; S. Zhu; D.C. Rockey A crucial role for GRK2 in regulation of endothelial cell nitric oxide synthase function in portal hypertension, Nat. Med., Volume 11 (2005), pp. 952-958

[18] M.C. Jimenez-Sainz; C. Murga; A. Kavelaars; M. Jurado-Pueyo; B.F. Kraastak; C.J. Heijnen; F. Mayor; A.M. Aragay G protein-coupled receptor kinase 2 negatively regulates chemokine signaling at a level downstream from G protein subunits, Mol. Biol. Cell, Volume 17 (2006), pp. 25-31

[19] N.L. Komarova; X. Zou; Q. Nie; L. Bardwell A theoretical framework for specificity in cell signaling, Mol. Syst. Biol., Volume 1 (2005), p. 2005.0023

[20] M.N. McClean; A. Mody; J.R. Broach; S. Ramanathan Cross-talk and decision making in MAP kinase pathways, Nat. Genet., Volume 39 (2007), pp. 409-414

[21] M. Behar; N. Hao; H.G. Dohlman; T.C. Elston Dose-to-duration encoding and signaling beyond saturation in intracellular signaling networks, PLoS Comput. Biol., Volume 4 (2008), p. e1000197

[22] X.R. Ren; E. Reiter; S. Ahn; J. Kim; W. Chen; R.J. Lefkowitz Different G protein-coupled receptor kinases govern G protein and {beta}-arrestin-mediated signaling of V2 vasopressin receptor, Proc. Natl. Acad. Sci. USA, Volume 102 (2005), pp. 1448-1453

[23] S.K. Shenoy; M.T. Drake; C.D. Nelson; D.A. Houtz; K. Xiao; S. Madabushi; E. Reiter; R.T. Premont; O. Lichtarge; R.J. Lefkowitz Beta-arrestin-dependent, G protein-independent ERK1/2 activation by the beta2 adrenergic receptor, J. Biol. Chem., Volume 281 (2006), pp. 1261-1273

[24] S. Ahn; S.K. Shenoy; H. Wei; R.J. Lefkowitz Differential kinetic and spatial patterns of beta-arrestin and G protein-mediated ERK activation by the angiotensin II receptor, J. Biol. Chem., Volume 279 (2004), pp. 35518-35525

[25] H. Wei; S. Ahn; S.K. Shenoy; S.S. Karnik; L. Hunyady; L.M. Luttrell; R.J. Lefkowitz Independent beta-arrestin 2 and G protein-mediated pathways for angiotensin II activation of extracellular signal-regulated kinases 1 and 2, Proc. Natl. Acad. Sci. USA, Volume 100 (2003), pp. 10782-10787

[26] E. Kara; P. Crepieux; C. Gauthier; N. Martinat; V. Piketty; F. Guillou; E. Reiter A phosphorylation cluster of five serine and threonine residues in the C-terminus of the follicle-stimulating hormone receptor is important for desensitization but not for beta-arrestin-mediated ERK activation, Mol. Endocrinol., Volume 20 (2006), pp. 3014-3026

[27] J.D. Violin; R.J. Lefkowitz Beta-arrestin-biased ligands at seven-transmembrane receptors, Trends Pharmacol. Sci., Volume 28 (2007), pp. 416-422

[28] T. Kenakin Collateral efficacy in drug discovery: Taking advantage of the good (allosteric) nature of 7TM receptors, Trends Pharmacol. Sci., Volume 28 (2007), pp. 407-415

[29] D. Gesty-Palmer; M. Chen; E. Reiter; S. Ahn; C.D. Nelson; S. Wang; A.E. Eckhardt; C.L. Cowan; R.F. Spurney; L.M. Luttrell; R.J. Lefkowitz Distinct beta-arrestin- and G protein-dependent pathways for parathyroid hormone receptor-stimulated ERK1/2 activation, J. Biol. Chem., Volume 281 (2006), pp. 10856-10864

[30] M.T. Drake; J.D. Violin; E.J. Whalen; J.W. Wisler; S.K. Shenoy; R.J. Lefkowitz Beta-arrestin-biased agonism at the beta2-adrenergic receptor, J. Biol. Chem., Volume 283 (2008), pp. 5669-5676

[31] J.W. Wisler; S.M. DeWire; E.J. Whalen; J.D. Violin; M.T. Drake; S. Ahn; S.K. Shenoy; R.J. Lefkowitz A unique mechanism of beta-blocker action: Carvedilol stimulates beta-arrestin signaling, Proc. Natl. Acad. Sci. USA, Volume 104 (2007), pp. 16657-16662

[32] N. Hao; M. Behar; T.C. Elston; H.G. Dohlman Systems biology analysis of G protein and MAP kinase signaling in yeast, Oncogene, Volume 26 (2007), pp. 3254-3266

[33] N. Hao; M. Behar; S.C. Parnell; M.P. Torres; C.H. Borchers; T.C. Elston; H.G. Dohlman A systems-biology analysis of feedback inhibition in the Sho1 osmotic-stress-response pathway, Curr. Biol., Volume 17 (2007), pp. 659-667

[34] J.J. Linderman Modeling of G-protein-coupled receptor signaling pathways, J. Biol. Chem., Volume 284 (2009), pp. 5427-5431

[35] L. Von Bertalanffy The theory of open systems in physics and biology, Science, Volume 111 (1950), pp. 23-29

[36] U. Stelzl; U. Worm; M. Lalowski; C. Haenig; F.H. Brembeck; H. Goehler; M. Stroedicke; M. Zenkner; A. Schoenherr; S. Koeppen; J. Timm; S. Mintzlaff; C. Abraham; N. Bock; S. Kietzmann; A. Goedde; E. Toksoz; A. Droege; S. Krobitsch; B. Korn; W. Birchmeier; H. Lehrach; E.E. Wanker A human protein–protein interaction network: A resource for annotating the proteome, Cell, Volume 122 (2005), pp. 957-968

[37] R. Linding; L.J. Jensen; G.J. Ostheimer; M.A. van Vugt; C. Jorgensen; I.M. Miron; F. Diella; K. Colwill; L. Taylor; K. Elder; P. Metalnikov; V. Nguyen; A. Pasculescu; J. Jin; J.G. Park; L.D. Samson; J.R. Woodgett; R.B. Russell; P. Bork; M.B. Yaffe; T. Pawson Systematic discovery of in vivo phosphorylation networks, Cell, Volume 129 (2007), pp. 1415-1426

[38] E.S. Witze; W.M. Old; K.A. Resing; N.G. Ahn Mapping protein post-translational modifications with mass spectrometry, Nat. Methods, Volume 4 (2007), pp. 798-806

[39] J.V. Olsen; B. Blagoev; F. Gnad; B. Macek; C. Kumar; P. Mortensen; M. Mann Global, in vivo, and site-specific phosphorylation dynamics in signaling networks, Cell, Volume 127 (2006), pp. 635-648

[40] B. Spurrier; S. Ramalingam; S. Nishizuka Reverse-phase protein lysate microarrays for cell signaling analysis, Nat. Protoc., Volume 3 (2008), pp. 1796-1808

[41] S.M. Chan; J. Ermann; L. Su; C.G. Fathman; P.J. Utz Protein microarrays for multiplex analysis of signal transduction pathways, Nat. Med., Volume 10 (2004), pp. 1390-1396

[42] D.S. Gembitsky; K. Lawlor; A. Jacovina; M. Yaneva; P. Tempst A prototype antibody microarray platform to monitor changes in protein tyrosine phosphorylation, Mol. Cell. Proteomics, Volume 3 (2004), pp. 1102-1118

[43] K.M. Sheehan; V.S. Calvert; E.W. Kay; Y. Lu; D. Fishman; V. Espina; J. Aquino; R. Speer; R. Araujo; G.B. Mills; L.A. Liotta; E.F. Petricoin; J.D. Wulfkuhle Use of reverse phase protein microarrays and reference standard development for molecular network analysis of metastatic ovarian carcinoma, Mol. Cell. Proteomics, Volume 4 (2005), pp. 346-355

[44] M.J. Lohse; V.O. Nikolaev; P. Hein; C. Hoffmann; J.P. Vilardaga; M. Bunemann Optical techniques to analyze real-time activation and signaling of G-protein-coupled receptors, Trends Pharmacol. Sci., Volume 29 (2008), pp. 159-165

[45] M.D. Allen; L.M. DiPilato; B. Ananthanarayanan; R.H. Newman; Q. Ni; J. Zhang Dynamic visualization of signaling activities in living cells, Sci. Signal., Volume 1 (2008) (pt6)

[46] B.N. Giepmans; S.R. Adams; M.H. Ellisman; R.Y. Tsien The fluorescent toolbox for assessing protein location and function, Science, Volume 312 (2006), pp. 217-224

[47] P.G. Charest; S. Terrillon; M. Bouvier Monitoring agonist-promoted conformational changes of beta-arrestin in living cells by intramolecular BRET, EMBO Rep., Volume 6 (2005), pp. 334-340

[48] C. Gales; R.V. Rebois; M. Hogue; P. Trieu; A. Breit; T.E. Hebert; M. Bouvier Real-time monitoring of receptor and G-protein interactions in living cells, Nat. Methods, Volume 2 (2005), pp. 177-184

[49] F. Azizi; C.H. Mastrangelo Generation of dynamic chemical signals with pulse code modulators, Lab Chip, Volume 8 (2008), pp. 907-912

[50] P.J. Hung; P.J. Lee; P. Sabounchi; R. Lin; L.P. Lee Continuous perfusion microfluidic cell culture array for high-throughput cell-based assays, Biotechnol. Bioeng., Volume 89 (2005), pp. 1-8

[51] H. Erfle; B. Neumann; P. Rogers; J. Bulkescher; J. Ellenberg; R. Pepperkok Work flow for multiplexing siRNA assays by solid-phase reverse transfection in multiwell plates, J. Biomol. Screen, Volume 13 (2008), pp. 575-580

[52] V. Starkuviene; R. Pepperkok; H. Erfle Transfected cell microarrays: An efficient tool for high-throughput functional analysis, Expert Rev. Proteomics, Volume 4 (2007), pp. 479-489

[53] D. Castel; M.A. Debily; A. Pitaval; X. Gidrol Cell microarray for functional exploration of genomes, Methods Mol. Biol., Volume 381 (2007), pp. 375-384

[54] S. Nelander; W. Wang; B. Nilsson; Q.B. She; C. Pratilas; N. Rosen; P. Gennemark; C. Sander Models from experiments: Combinatorial drug perturbations of cancer cells, Mol. Syst. Biol., Volume 4 (2008), p. 216

[55] M. Hucka; A. Finney; H.M. Sauro; H. Bolouri; J.C. Doyle; H. Kitano; A.P. Arkin; B.J. Bornstein; D. Bray; A. Cornish-Bowden; A.A. Cuellar; S. Dronov; E.D. Gilles; M. Ginkel; V. Gor Goryanin; W.J. Hedley; T.C. Hodgman; J.H. Hofmeyr; P.J. Hunter; N.S. Juty; J.L. Kasberger; A. Kremling; U. Kummer; N. Le Novere; L.M. Loew; D. Lucio; P. Mendes; E. Minch; E.D. Mjolsness; Y. Nakayama; M.R. Nelson; P.F. Nielsen; T. Sakurada; J.C. Schaff; B.E. Shapiro; T.S. Shimizu; H.D. Spence; J. Stelling; K. Takahashi; M. Tomita; J. Wagner; J. Wang The systems biology markup language (SBML): A medium for representation and exchange of biochemical network models, Bioinformatics, Volume 19 (2003), pp. 524-531

[56] H. Kitano; A. Funahashi; Y. Matsuoka; K. Oda Using process diagrams for the graphical representation of biological networks, Nat. Biotechnol., Volume 23 (2005), pp. 961-966

[57] A. Funahashi, Y. Matsuoka, A. Jouraku, M. Morohashi, N. Kikuchi, H. Kitano, CellDesigner 3.5: A versatile modeling tool for biochemical networks, in: Proceedings of the IEEE, 2008, pp. 1254–1265

[58] L. Calzone; F. Fages; S. Soliman BIOCHAM: An environment for modeling biological systems and formalizing experimental knowledge, Bioinformatics, Volume 22 (2006), pp. 1805-1807

[59] F. Fages; S. Soliman From reaction models to influence graphs and back: A theorem, Formal Methods in Systems Biology, LNBI, vol. 5054, Springer-Verlag, 2008, pp. 90-102

[60] F. Fages; S. Soliman Formal cell biology in BIOCHAM, International School on Formal Methods for the Design of Computer, Communication and Software Systems: Computational Systems Biology, Bertinoro, Italy, LNCS, vol. 5016, Springer-Verlag, 2008, pp. 54-80

[61] D. Gilbert; H. Fuss; X. Gu; R. Orton; S. Robinson; V. Vyshemirsky; M.J. Kurth; C.S. Downes; W. Dubitzky Computational methodologies for modelling, analysis and simulation of signalling networks, Brief Bioinform., Volume 7 (2006), pp. 339-353

[62] R.J. Orton; O.E. Sturm; V. Vyshemirsky; M. Calder; D.R. Gilbert; W. Kolch Computational modelling of the receptor-tyrosine-kinase-activated MAPK pathway, Biochem. J., Volume 392 (2005), pp. 249-261

[63] A. Rizk; G. Batt; F. Fages; S. Soliman On a continuous degree of satisfaction of temporal logic formulae with applications to systems biology, Computational Methods in Systems Biology CMSB'08, LNBI, vol. 5307, Springer-Verlag, 2008, pp. 251-268

[64] N. Hansen; A. Ostermeier Completely derandomized self-adaptation in evolution strategies, Evol. Comput., Volume 9 (2001), pp. 159-195

[65] A. Rizk; G. Batt; F. Fages; S. Soliman A general computational method for robustness analysis with applications to synthetic gene networks, Bioinformatics, Volume 25 (2009), pp. 169-178


Comments - Policy