Julio Vera

A systems biology approach to analyse amplification in the JAK2-STAT5 signalling pathway

BackgroundThe amplification of signals, defined as an increase in the intensity of a signal through networks of intracellular reactions, is considered one of the essential properties in many cell signalling pathways. Despite of the apparent importance of signal amplification, there have been few attempts to formalise this concept.ResultsIn this work we investigate the amplification and responsiveness of the JAK2-STAT5 pathway using a kinetic model. The recruitment of EpoR to the plasma membrane, activation by Epo, and deactivation of the EpoR/JAK2 complex are considered as well as the activation and nucleocytoplasmic shuttling of STAT5. Using qualitative biological knowledge, we first establish the structure of a general power-law model. We then generate a family of models from which we select suitable candidates. The parameter values of the model are estimated from experimental quantitative time-course data. The final model, whether it is conventional model with fixed predefined integer kinetic orders or a model with variable non-integer kinetic orders, is selected on the basis of a good agreement between simulations and the experimental data. The model is used to analyse the responsiveness and amplification properties of the pathway with sustained, transient, and oscillatory stimulation.ConclusionThe selected kinetic model predicts that the system acts as an amplifier with maximum amplification and sensitivity for input signals whose intensity match physiological values for Epo concentration and with duration in the range of one to 100 minutes. The response of the system reaches saturation for more intense and longer stimulation with Epo. We hypothesise that these properties of the system directly relate to the saturation of Epo receptor activation, its low recruitment to the plasma membrane and intense deactivation as predicted by the model.

MicroRNA-regulated networks: the perfect storm for classical molecular biology, the ideal scenario for systems biology.

MicroRNAs (miRNAs) are involved in many regulatory pathways some of which are complex networks enriched in regulatory motifs like positive or negative feedback loops or coherent and incoherent feedforward loops. Their complexity makes the understanding of their regulation difficult and the interpretation of experimental data cumbersome. In this book chapter we claim that systems biology is the appropriate approach to investigate the regulation of these miRNA-regulated networks. Systems biology is an interdisciplinary approach by which biomedical questions on biochemical networks are addressed by integrating experiments with mathematical modelling and simulation. We here introduce the foundations of the systems biology approach, the basic theoretical and computational tools used to perform model-based analyses of miRNA-regulated networks and review the scientific literature in systems biology of miRNA regulation, with a focus on cancer.

Systems biologists seek fuller integration of systems biology approaches in new cancer research programs

Systems biology takes an interdisciplinary approach to the systematic study of complex interactions in biological systems. This approach seeks to decipher the emergent behaviors of complex systems rather than focusing only on their constituent properties. As an increasing number of examples illustrate the value of systems biology approaches to understand the initiation, progression, and treatment of cancer, systems biologists from across Europe and the United States hope for changes in the way their field is currently perceived among cancer researchers. In a recent EU-US workshop, supported by the European Commission, the German Federal Ministry for Education and Research, and the National Cancer Institute of the NIH, the participants discussed the strengths, weaknesses, hurdles, and opportunities in cancer systems biology.

Kinetic Modeling-Based Detection of Genetic Signatures that Provide Chemoresistance via the E2F1-p73/DNp73-miR-205 Network

Drug resistance is a major cause of deaths from cancer. E2F1 is a transcription factor involved in cell proliferation, apoptosis. and metastasis through an intricate regulatory network, which includes other transcription factors like p73 and cancer-related microRNAs like miR-205. To investigate the emergence of drug resistance, we developed a methodology that integrates experimental data with a network biology and kinetic modeling. Using a regulatory map developed to summarize knowledge on E2F1 and its interplay with p73/DNp73 and miR-205 in cancer drug responses, we derived a kinetic model that represents the network response to certain genotoxic and cytostatic anticancer drugs. By perturbing the model parameters, we simulated heterogeneous cell configurations referred to as in silico cell lines. These were used to detect genetic signatures characteristic for single or double drug resistance. We identified a signature composed of high E2F1 and low miR-205 expression that promotes resistance to genotoxic drugs. In this signature, downregulation of miR-205, can be mediated by an imbalance in the p73/DNp73 ratio or by dysregulation of other cancer-related regulators of miR-205 expression such as TGFβ-1 or TWIST1. In addition, we found that a genetic signature composed of high E2F1, low miR-205, and high ERBB3 can render tumor cells insensitive to both cytostatic and genotoxic drugs. Our model simulations also suggested that conventional genotoxic drug treatment favors selection of chemoresistant cells in genetically heterogeneous tumors, in a manner requiring dysregulation of incoherent feedforward loops that involve E2F1, p73/DNp73, and miR-205.

Optimization of biochemical systems through mathematical programming: Methods and applications

In this work we present a general (mono and multiobjective) optimization framework for the technological improvement of biochemical systems. The starting point of the method is a mathematical model in ordinary differential equations (ODEs) of the investigated system, based on qualitative biological knowledge and quantitative experimental data. In the method we take advantage of the special structural features of a family of ODEs called power-law models to reduce the computational complexity of the optimization program. In this way, the genetic manipulation of a biochemical system to meet a certain biotechnological goal can be expressed as an optimization program with some desired properties such as linearity or convexity. The general method of optimization is presented and discussed in its linear and geometric programming versions. We furthermore illustrate the use of the method by several real case studies. We conclude that the technological improvement of microorganisms can be afforded using the combination of mathematical modelling and optimization. The systematic nature of this approach facilitates the redesign of biochemical systems and makes this a predictive exercise rather than a trial-and-error procedure.

Understanding microRNA-mediated gene regulatory networks through mathematical modelling.

The discovery of microRNAs (miRNAs) has added a new player to the regulation of gene expression. With the increasing number of molecular species involved in gene regulatory networks, it is hard to obtain an intuitive understanding of network dynamics. Mathematical modelling can help dissecting the role of miRNAs in gene regulatory networks, and we shall here review the most recent developments that utilise different mathematical modelling approaches to provide quantitative insights into the function of miRNAs in the regulation of gene expression. Key miRNA regulation features that have been elucidated via modelling include: (i) the role of miRNA-mediated feedback and feedforward loops in fine-tuning of gene expression; (ii) the miRNA–target interaction properties determining the effectiveness of miRNA-mediated gene repression; and (iii) the competition for shared miRNAs leading to the cross-regulation of genes. However, there is still lack of mechanistic understanding of many other properties of miRNA regulation like unconventional miRNA–target interactions, miRNA regulation at different sub-cellular locations and functional miRNA variant, which will need future modelling efforts to deal with. This review provides an overview of recent developments and challenges in this field.

Systems biology of JAK-STAT signalling in human malignancies.

Originally implicated in the regulation of survival, proliferation and differentiation of haematopoietic cells, the JAK-STAT pathway has also been linked to developmental processes, growth control and maintenance of homeostasis in a variety of other cells and tissues. Although it remains a complex system, its relative simplicity and the availability of molecular data makes it particularly attractive for modelling approaches. In this review, we will focus on JAK-STAT signalling in the context of cancer and present efforts to investigate signalling dynamics with the help of mathematical models. We describe the modelling workflow that realises a systems biology approach and give an example for interferon-γ signalling in pancreatic stellate cells.

Cooperative gene regulation by microRNA pairs and their identification using a computational workflow

MicroRNAs (miRNAs) are an integral part of gene regulation at the post-transcriptional level. Recently, it has been shown that pairs of miRNAs can repress the translation of a target mRNA in a cooperative manner, which leads to an enhanced effectiveness and specificity in target repression. However, it remains unclear which miRNA pairs can synergize and which genes are target of cooperative miRNA regulation. In this paper, we present a computational workflow for the prediction and analysis of cooperating miRNAs and their mutual target genes, which we refer to as RNA triplexes. The workflow integrates methods of miRNA target prediction; triplex structure analysis; molecular dynamics simulations and mathematical modeling for a reliable prediction of functional RNA triplexes and target repression efficiency. In a case study we analyzed the human genome and identified several thousand targets of cooperative gene regulation. Our results suggest that miRNA cooperativity is a frequent mechanism for an enhanced target repression by pairs of miRNAs facilitating distinctive and fine-tuned target gene expression patterns. Human RNA triplexes predicted and characterized in this study are organized in a web resource at www.sbi.uni-rostock.de/triplexrna/.

A system biology approach to understand functional activity of cell communication systems.

Systems Biology is the quantitative analysis of dynamic interactions among several components of a biochemical system, aiming at an understanding of the behavior of the system as a whole. From an experimental perspective, systems biology is a suitable tool to support the biologist in the process of hypotheses generation and the efficient design of experiments. In this chapter, we discuss the elements of a systems biology methodology based on the interaction between experimental biologists and theoreticians. We, furthermore, show the use of such a methodology in a case study, analyzing receptor and transcription factor modulation affecting the responsiveness of the JAK2/STAT5 pathway.

Dynamic properties of a delayed protein cross talk model.

In this paper we investigate how the inclusion of time delay alters the dynamical properties of the Jacob-Monod model, describing the control of the beta-galactosidase synthesis by the lac repressor protein in E. coli. The consequences of a time delay on the dynamics of this system are analysed using Hopfs theorem and Lyapunov-Andronovs theory applied to the original mathematical model and to an approximated version. Our analytical calculations predict that time delay acts as a key bifurcation parameter. This is confirmed by numerical simulations. A critical value of time delay, which depends on the values of the model parameters, is analytically established. Around this critical value, the properties of the system change drastically, allowing under certain conditions the emergence of stable limit cycles, that is self-sustained oscillations. In addition, the features of the end product repression play an essential role in the characterisation of these limit cycles: if cooperativity is considered in the end product repression, time delay higher than the mentioned critical value induce differentiated responses during the oscillations, provoking cycles of all-or-nothing response in the concentration of the species.