A multidisciplinary survey of modeling techniques for biochemical networks

All processes of life are dominated by networks of interacting biochemical components. The purpose of modeling these networks is manifold. From a theoretical point of view it allows the exploration of network structures and dynamics, to find emergent properties or to explain the organization and evolution of networks. From a practical point of view, in silico experiments can be performed that would be very expensive or impossible to achieve in the laboratory, such as hypothesis-testing with regard to knockout experiments or overexpression, or checking the validity of a proposed molecular mechanism. The literature on modeling biochemical networks is growing rapidly and the motivations behind different modeling techniques are sometimes quite distant from each other. To clarify the current context, we present a systematic overview of the different philosophies to model biochemical networks. We put particular emphasis on three main domains which have been playing a major role in the past, namely: mathematics with ordinary and partial differential equations, statistics with stochastic simulation algorithms, Bayesian networks and Markov chains, and the field of computer science with process calculi, term rewriting systems and state based systems. For each school, we evaluate advantages and disadvantages such as the granularity of representation, scalability, accessibility or availability of analysis tools. Following this, we describe how one can combine some of those techniques and thus take advantages of several techniques through the use of bridging tools. Finally, we propose a next step for modeling biochemical networks by using artificial chemistries and evolutionary computation. This work was funded by ESIGNET (Evolving Cell Signaling Networks in Silico), an European Integrated Project in the EU FP6 NEST Initiative (contract no. 12789)

Mathematical model of a telomerase transcriptional regulatory network developed by cell-based screening: analysis of inhibitor effects and telomerase expression mechanisms

Cancer cells depend on transcription of telomerase reverse transcriptase (TERT). Many transcription factors affect TERT, though regulation occurs in context of a broader network. Network effects on telomerase regulation have not been investigated, though deeper understanding of TERT transcription requires a systems view. However, control over individual interactions in complex networks is not easily achievable. Mathematical modelling provides an attractive approach for analysis of complex systems and some models may prove useful in systems pharmacology approaches to drug discovery. In this report, we used transfection screening to test interactions among 14 TERT regulatory transcription factors and their respective promoters in ovarian cancer cells. The results were used to generate a network model of TERT transcription and to implement a dynamic Boolean model whose steady states were analysed. Modelled effects of signal transduction inhibitors successfully predicted TERT repression by Src-family inhibitor SU6656 and lack of repression by ERK inhibitor FR180204, results confirmed by RT-QPCR analysis of endogenous TERT expression in treated cells. Modelled effects of GSK3 inhibitor 6-bromoindirubin-3′-oxime (BIO) predicted unstable TERT repression dependent on noise and expression of JUN, corresponding with observations from a previous study. MYC expression is critical in TERT activation in the model, consistent with its well known function in endogenous TERT regulation. Loss of MYC caused complete TERT suppression in our model, substantially rescued only by co-suppression of AR. Interestingly expression was easily rescued under modelled Ets-factor gain of function, as occurs in TERT promoter mutation. RNAi targeting AR, JUN, MXD1, SP3, or TP53, showed that AR suppression does rescue endogenous TERT expression following MYC knockdown in these cells and SP3 or TP53 siRNA also cause partial recovery. The model therefore successfully predicted several aspects of TERT regulation including previously unknown mechanisms. An extrapolation suggests that a dominant stimulatory system may programme TERT for transcriptional stability

UNDERSTANDING CONDITIONAL MODES OF ACTIONS IN CHEMICAL-INDUCED TOXICITY USING RULE MODELS

It is estimated that 115 million animals are used in experimental testing each year. Hence, shifting efforts toward alternative methods for toxicity assessment is essential. However, slow regulatory acceptance of new approaches is governed by knowledge gaps in toxicity modes of action. In this thesis, I describe these challenges and the use of in vitro screening as an alternative of animal testing. I also discuss common data-based methods to derive hypotheses about toxicity modes of actions, and the associated limitations in capturing multiple biological perturbations. I applied novel data-based workflows, using rule models, to prioritize in vitro assays predictive of toxicity as well as to detect significant polypharmacology profiles. I explain how constraints were applied to rule-based models to inform meaningful mechanistic interpretation for two toxicity endpoints: rat hepatotoxicity and acute toxicity. I compared assays selected, by rules, for predicting hepatotoxicity with endpoints used in in vitro models from commercial sources. An overlap was observed including cytochrome activity, mitochondrial toxicity and immunological responses. However, nuclear receptor activity, identified in rules, is not currently covered in commercial setups. I also demonstrate that endocrine disruption endpoints extrapolate better into in vivo toxicity when a set of specific conditions are met, such as physicochemical properties associated with good bioavailability. Next, I examined synergistic interactions between conditions in rules describing acute toxicity. I gained novel insights into how specific stressors potentiate the perturbation by known key events, such as acetylcholinesterase inhibition and neuro-signalling disruption. I show that examining polypharmacology profiles is particularly important at low bioactive potencies. Further, the overall predictive performance of rules describing acute toxicity was tested against a benchmark Random Forest model in a conformal prediction framework. Irrespective to the data type used in the training, the models were prone to bias over compounds promiscuity, by which high promiscuous compounds were more likely to be predicted as toxic. Overall, the studies conducted in this thesis provide novel insights into molecular mechanisms of toxicity, namely hepatotoxicity and acute toxicity, and with regards to chemical properties and polypharmacology. This knowledge can be used to improve the utility and design of alternative methods for toxicity, and hence, accelerate the regulatory acceptance.Islamic Development Bank Cambridge Trust Fun

Application of Petri net based analysis techniques to signal transduction pathways

BACKGROUND: Signal transduction pathways are usually modelled using classical quantitative methods, which are based on ordinary differential equations (ODEs). However, some difficulties are inherent in this approach. On the one hand, the kinetic parameters involved are often unknown and have to be estimated. With increasing size and complexity of signal transduction pathways, the estimation of missing kinetic data is not possible. On the other hand, ODEs based models do not support any explicit insights into possible (signal-) flows within the network. Moreover, a huge amount of qualitative data is available due to high-throughput techniques. In order to get information on the systems behaviour, qualitative analysis techniques have been developed. Applications of the known qualitative analysis methods concern mainly metabolic networks. Petri net theory provides a variety of established analysis techniques, which are also applicable to signal transduction models. In this context special properties have to be considered and new dedicated techniques have to be designed. METHODS: We apply Petri net theory to model and analyse signal transduction pathways first qualitatively before continuing with quantitative analyses. This paper demonstrates how to build systematically a discrete model, which reflects provably the qualitative biological behaviour without any knowledge of kinetic parameters. The mating pheromone response pathway in Saccharomyces cerevisiae serves as case study. RESULTS: We propose an approach for model validation of signal transduction pathways based on the network structure only. For this purpose, we introduce the new notion of feasible t-invariants, which represent minimal self-contained subnets being active under a given input situation. Each of these subnets stands for a signal flow in the system. We define maximal common transition sets (MCT-sets), which can be used for t-invariant examination and net decomposition into smallest biologically meaningful functional units. CONCLUSION: The paper demonstrates how Petri net analysis techniques can promote a deeper understanding of signal transduction pathways. The new concepts of feasible t-invariants and MCT-sets have been proven to be useful for model validation and the interpretation of the biological system behaviour. Whereas MCT-sets provide a decomposition of the net into disjunctive subnets, feasible t-invariants describe subnets, which generally overlap. This work contributes to qualitative modelling and to the analysis of large biological networks by their fully automatic decomposition into biologically meaningful modules

Bioinformatics and mathematical modelling in the study of receptor-receptor interactions and receptor oligomerization: focus on adenosine receptors.

none8sìThe concept of intra-membrane receptor-receptor interactions (RRIs) between different types of G protein-coupled receptors (GPCRs) and evidence for their existence was introduced by Agnati and Fuxe in 1980/81 through the biochemical analysis of the effects of neuropeptides on the binding characteristics of monoamine receptors in membrane preparations from discrete brain regions and functional studies of the interactions between neuropeptides and monoamines in the control of specific functions such as motor control and arterial blood pressure control in animal models. Whether GPCRs can form high-order structures is still a topic of an intense debate. Increasing evidence, however, suggests that the hypothesis of the existence of high-order receptor oligomers is correct. A fundamental consequence of the view describing GPCRs as interacting structures, with the likely formation at the plasma membrane of receptor aggregates of multiple receptors (Receptor Mosaics) is that it is no longer possible to describe signal transduction simply as the result of the binding of the chemical signal to its receptor, but rather as the result of a filtering/integration of chemical signals by the Receptor Mosaics (RMs) and membrane-associated proteins. Thus, in parallel with experimental research, significant efforts were spent in bioinformatics and mathematical modelling. We review here the main approaches that have been used to assess the interaction interfaces allowing the assembly of GPCRs and to shed some light on the integrative functions emerging from the complex behaviour of these RMs. Particular attention was paid to the RMs generated by adenosine A(2A), dopamine D-2, cannabinoid CB1, and metabotropic glutamate mGlu(5) receptors (A(2A). D-2, CB1, and mGlu(5), respectively), and a possible approach to model the interplay between the D-2-A(2A)-CB1 and D-2-A(2A)-mGlu(5) trimers is proposed. This article is part of a Special Issue entitled: "Adenosine Receptors". (C) 2010 Elsevier B.V. All rights reserved.openD. GUIDOLIN; F. CIRUELA; S. GENEDANI; M. GUESCINI; C. TORTORELLA; G. ALBERTIN; K. FUXE; L.F. AGNATID., Guidolin; F., Ciruela; S., Genedani; Guescini, Michele; C., Tortorella; G., Albertin; K., Fuxe; L. F., Agnat

pdCSM-GPCR: predicting potent GPCR ligands with graph-based signatures.

Funder: Newton Fund RCUK-CONFAPFunder: Victorian Government’s Operational Infrastructure Support ProgramMOTIVATION: G protein-coupled receptors (GPCRs) can selectively bind to many types of ligands, ranging from light-sensitive compounds, ions, hormones, pheromones and neurotransmitters, modulating cell physiology. Considering their role in many essential cellular processes, they are one of the most targeted protein families, with over a third of all approved drugs modulating GPCR signalling. Despite this, the large diversity of receptors and their multipass transmembrane architectures make the identification and development of novel specific, and safe GPCR ligands a challenge. While computational approaches have the potential to assist GPCR drug development, they have presented limited performance and generalization capabilities. Here, we explored the use of graph-based signatures to develop pdCSM-GPCR, a method capable of rapidly and accurately screening potential GPCR ligands. RESULTS: Bioactivity data (IC50, EC50, Ki and Kd) for individual GPCRs were curated. After curation, we used the data for developing predictive models for 36 major GPCR targets, across 4 classes (A, B, C and F). Our models compose the most comprehensive computational resource for GPCR bioactivity prediction to date. Across stratified 10-fold cross-validation and blind tests, our approach achieved Pearson's correlations of up to 0.89, significantly outperforming previous methods. Interpreting our results, we identified common important features of potent GPCRs ligands, which tend to have bicyclic rings, leading to higher levels of aromaticity. We believe pdCSM-GPCR will be an invaluable tool to assist screening efforts, enriching compound libraries and ranking candidates for further experimental validation. AVAILABILITY AND IMPLEMENTATION: pdCSM-GPCR predictive models and datasets used have been made available via a freely accessible and easy-to-use web server at http://biosig.unimelb.edu.au/pdcsm_gpcr/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics Advances online