Due in large measure to the explosive progress in molecular biology, biology has become arguably the most exciting scientific field. The first half of the 21st century is sometimes referred to as the 'era of biology', analogous to the first half of the 20th century, which was considered to be the 'era of physics'. Yet, biology is facing a crisis--or is it an opportunity--reminiscent of the state of biology in pre-double-helix time. The principal challenge facing systems biology is complexity. According to Hood, 'Systems biology defines and analyses the interrelationships of all of the elements in a functioning system in order to understand how the system works.' With 30000+ genes in the human genome the study of all relationships simultaneously becomes a formidably complex problem. Hanahan and Weinberg raised the question as to whether progress will consist of 'adding further layers of complexity to a scientific literature that is already complex almost beyond measure' or whether the progress will lead to a 'science with a conceptual structure and logical coherence that rivals that of chemistry or physics.' At the core of the challenge is the need for a new approach, a shift from reductionism to a holistic perspective. However, more than just a pronouncement of a new approach is needed. We suggest that what is needed is to provide a conceptual framework for systems biology research. We propose that the concept of a complex system, i.e. a system of systems as defined in mathematical general systems theory (MGST), is central to provide such a framework. We further argue that for a deeper understanding in systems biology investigations should go beyond building numerical mathematical or computer models--important as they are. Biological phenomena cannot be predicted with the level of numerical precision as in classical physics. Explanations in terms of how the categories of systems are organised to function in ever changing conditions are more revealing. Non-numerical mathematical tools are appropriate for the task. Such a categorical perspective led us to propose that the core of understanding in systems biology depends on the search for organising principles rather than solely on construction of predictive descriptions (i.e. models) that exactly outline the evolution of systems in space and time. The search for organising principles requires an identification/discovery of new concepts and hypotheses. Some of them, such as coordination motifs for transcriptional regulatory networks and bounded autonomy of levccels in a hierarchy, are outlined in this article. Experimental designs are outlined to help verify the applicability of the interaction balance principle of coordination to transcriptional and posttranscriptional networks.

In the present chapter we discuss methodologies for the modelling, calibration and validation of cellular signalling pathway dynamics. The discussion begins with the typical range of techniques for modelling that might be employed to go from the chemical kinetics to a mathematical model of biochemical pathways. In particular, we consider the decision-making processes involved in selecting the right mechanism and level of detail of representation of the biochemical interactions. These include the choice between (i) deterministic and stochastic chemical kinetics representations, (ii) discrete and continuous time models and (iii) representing continuous and discrete state processes. We then discuss the task of calibrating the models using information available in web-based databases. For situations in which the data are not available from existing sources we discuss model calibration based upon measured data and system identification methods. Such methods, together with mathematical modelling databases and computational tools, are often available in standard packages. We therefore make explicit mention of a range of popular and useful sites. As an example of the whole modelling and calibration process, we discuss a study of the cross-talk between the IL-1 (interleukin-1)-stimulated NF-kappaB (nuclear factor kappaB) pathway and the TGF-beta (transforming growth factor beta)-stimulated Smad2 pathway.

Few of >150 published cell cycle modeling efforts use significant levels of data for tuning and validation. This reflects the difficultly to generate correlated quantitative data, and it points out a critical uncertainty in modeling efforts. To develop a data-driven model of cell cycle regulation, we used contiguous, dynamic measurements over two time scales (minutes and hours) calculated from static multiparametric cytometry data. The approach provided expression profiles of cyclin A2, cyclin B1, and phospho-S10-histone H3. The model was built by integrating and modifying two previously published models such that the model outputs for cyclins A and B fit cyclin expression measurements and the activation of B cyclin/Cdk1 coincided with phosphorylation of histone H3. The model depends on Cdh1-regulated cyclin degradation during G1, regulation of B cyclin/Cdk1 activity by cyclin A/Cdk via Wee1, and transcriptional control of the mitotic cyclins that reflects some of the current literature. We introduced autocatalytic transcription of E2F, E2F regulated transcription of cyclin B, Cdc20/Cdh1 mediated E2F degradation, enhanced transcription of mitotic cyclins during late S/early G2 phase, and the sustained synthesis of cyclin B during mitosis. These features produced a model with good correlation between state variable output and real measurements. Since the method of data generation is extensible, this model can be continually modified based on new correlated, quantitative data.

In this work, we search for coordination as an organizing principle in a complex signaling system using a multilevel hierarchical paradigm. The objective is to explain the underlying mechanism of Interferon (IFN γ ) induced JAK-STAT (specifically JAK1/JAK2-STAT1) pathway behavior. Starting with a mathematical model of the pathway from the literature, we modularize the system using biological knowledge via principles of biochemical cohesion, biological significance, and functionality. The modularized system is then used as a basis for in silico inhibition, knockdown/ deletion and perturbation experiments to discover a coordination mechanism. Our analysis shows that a module representing the SOCS1 complex can be identified as the coordinator. Analysis of the coordinator can then be used for the selection of biological experiments for the discovery of 'soft' molecular drug targets, that could lead to the development of improved therapeutics. The coordinator identified is also being investigated to determine its relationship to pathological conditions.

BackgroundAn imprecise quantitative sense for the oscillating levels of proteins and their modifications, interactions, and translocations as a function of the cell cycle is fundamentally important for a cartoon/narrative understanding for how the cell cycle works. Mathematical modeling of the same cartoon/narrative models would be greatly enhanced by an open-ended methodology providing precise quantification of many proteins and their modifications, etc. Here we present methodology that fulfills these features.MethodologyMultiparametric flow cytometry was performed on Molt4 cells to measure cyclins A2 and B1, phospho-S10-histone H3, DNA content, and light scatter (cell size). The resulting 5 dimensional data were analyzed as a series of bivariate plots to isolate the data as segments of an N-dimensional “worm” through the data space. Sequential, unidirectional regions of the data were used to assemble expression profiles for each parameter as a function of cell frequency.ResultsAnalysis of synthesized data in which the true values where known validated the approach. Triplicate experiments demonstrated exceptional reproducibility. Comparison of three triplicate experiments stained by two methods (single cyclin or dual cyclin measurements with common DNA and phospho-histone H3 measurements) supported the feasibility of combining an unlimited number of epitopes through this methodology. The sequential degradations of cyclin A2 followed by cyclin B1 followed by de-phosphorylation of histone H3 were precisely mapped. Finally, a two phase expression rate during interphase for each cyclin was robustly identified.ConclusionsVery precise, correlated expression profiles for important cell cycle regulating and regulated proteins and their modifications can be produced, limited only by the number of available high-quality antibodies. These profiles can be assembled into large information libraries for calibration and validation of mathematical models.