Systems biology最新文献

Understanding the large reaction networks found in biological systems is a daunting task. One approach is to divide a network into more manageable smaller modules, thus simplifying the problem. This is a common strategy used in engineering. However, the process of identifying biological modules is still in its infancy and very little is understood about the range and capabilities of motif structures found in biological modules. In order to delineate these modules, a library of functional motifs has been generated via in silico evolution techniques. On the basis of their functional forms, networks were evolved from four broad areas: oscillators, bistable switches, homeostatic systems and frequency filters. Some of these motifs were constructed from simple mass action kinetics, others were based on Michaelis-Menten kinetics as found in protein/protein networks and the remainder were based on Hill equations as found in gene/protein interaction networks. The purpose of the study is to explore the capabilities of different network architectures and the rich variety of functional forms that can be generated. Ultimately, the library may be used to delineate functional motifs in real biological networks.

The transient kinetic behaviour of an open single enzyme, single substrate reaction is examined. The reaction follows the Van Slyke-Cullen mechanism, a spacial case of the Michaelis-Menten reaction. The analysis is performed both with and without applying the quasi-steady-state approximation. The analysis of the full system shows conditions for biochemical pathway coupling, which yield sustained oscillatory behaviour in the enzyme reaction. The reduced model does not demonstrate this behaviour. The results have important implications in the analysis of open biochemical reactions and the modelling of metabolic systems.

Large, complex data sets that are generated from microarray experiments, create a need for systematic analysis techniques to unravel the underlying connectivity of gene regulatory networks. A modular approach, previously proposed by Kholodenko and co-workers, helps to scale down the network complexity into more computationally manageable entities called modules. A functional module includes a gene's mRNA, promoter and resulting products, thus encompassing a large set of interacting states. The essential elements of this approach are described in detail for a three-gene model network and later extended to a ten-gene model network, demonstrating scalability. The network architecture is identified by analysing in silico steady-state changes in the activities of only the module outputs, communicating intermediates, that result from specific perturbations applied to the network modules one at a time. These steady-state changes form the system response matrix, which is used to compute the network connectivity or network interaction map. By employing a known biochemical network, the accuracy of the modular approach and its sensitivity to key assumptions are evaluated.

To unravel the complex in vivo regulatory interdependences of biochemical networks, experiments with the living organism are absolutely necessary. Stimulus response experiments (SREs) have become increasingly popular in recent years. The response of metabolite concentrations from all major parts of the central metabolism is monitored over time by modem analytical methods, producing several thousand data points. SREs are applied to determine enzyme kinetic parameters and to find unknown enzyme regulatory mechanisms. Owing to the complex regulatory structure of metabolic networks and the amount of measured data, the evaluation of an SRE has to be extensively supported by modelling. If the enzyme regulatory mechanisms are part of the investigation, a large number of models with different enzyme kinetics have to be tested for their ability to reproduce the observed behaviour. In this contribution, a systematic model-building process for data-driven exploratory modelling is introduced with the aim of discovering essential features of the biological system. The process is based on data pre-processing, correlation-based hypothesis generation, automatic model family generation, large-scale model selection and statistical analysis of the best-fitting models followed by an extraction of common features. It is illustrated by the example of the aromatic amino acid synthesis pathway in Escherichia coli.

The concepts and methods of systems biology are extended to neuropharmacology in order to test and design drugs for the treatment of neurological and psychiatric disorders. Computational modelling by integrating compartmental neural modelling techniques and detailed kinetic descriptions of pharmacological modulation of transmitter-receptor interaction is offered as a method to test the electrophysiological and behavioural effects of putative drugs. Even more, an inverse method is suggested as a method for controlling a neural system to realise a prescribed temporal pattern. In particular, as an application of the proposed new methodology, a computational platform is offered to analyse the generation and pharmacological modulation of theta rhythm related to anxiety.

Gene regulatory, signal transduction and metabolic networks are major areas of interest in the newly emerging field of systems biology. In living cells, stochastic dynamics play an important role; however, the kinetic parameters of biochemical reactions necessary for modelling these processes are often not accessible directly through experiments. The problem of estimating stochastic reaction constants from molecule count data measured, with error, at discrete time points is considered. For modelling the system, a hidden Markov process is used, where the hidden states are the true molecule counts, and the transitions between those states correspond to reaction events following collisions of molecules. Two different algorithms are proposed for estimating the unknown model parameters. The first is an approximate maximum likelihood method that gives good estimates of the reaction parameters in systems with few possible reactions in each sampling interval. The second algorithm, treating the data as exact measurements, approximates the number of reactions in each sampling interval by solving a simple linear equation. Maximising the likelihood based on these approximations can provide good results, even in complex reaction systems.

The ability to build genetic circuits with a reproducible response to external stimuli depends on the experimental and theoretical methods used in the process. A theoretical formalism that describes the response of a nonlinear stochastic genetic network to the external stimuli (input signals), is proposed. Two applications are studied in detail: the design of a logic pulse and the interference of three signal generators in the E2F1 regulatory element. The gene interactions are presented using molecular diagrams that have a precise mathematical structure and retain the biological meaning of the processes.

Biological systems that have been experimentally verified to be robust to significant changes in their environments require mathematical models that are themselves robust. In this context, a necessary condition for model robustness is that the model dynamics should not be sensitive to small variations in the model's parameters. Robustness analysis problems of this type have been extensively studied in the field of robust control theory and have been found to be very difficult to solve in general. The authors describe how some tools from robust control theory and nonlinear optimisation can be used to analyse the robustness of a recently proposed model of the molecular network underlying adenosine 3',5'-cyclic monophosphate (cAMP) oscillations observed in fields of chemotactic Dictyostelium cells. The network model, which consists of a system of seven coupled nonlinear differential equations, accurately reproduces the spontaneous oscillations in cAMP observed during the early development of D. discoideum. The analysis by the authors reveals, however, that very small variations in the model parameters can effectively destroy the required oscillatory dynamics. A biological interpretation of the analysis results is that correct functioning of a particular positive feedback loop in the proposed model is crucial to maintaining the required oscillatory dynamics.

Compared evaluation of different methods is presented for estimating missing values in microarray data: weighted K-nearest neighbours imputation (KNNimpute), regression-based methods such as local least squares imputation (LLSimpute) and partial least squares imputation (PLSimpute) and Bayesian principal component analysis (BPCA). The influence in prediction accuracy of some factors, such as methods' parameters, type of data relationships used in the estimation process (i.e. row-wise, column-wise or both), missing rate and pattern and type of experiment [time series (TS), non-time series (NTS) or mixed (MIX) experiments] is elucidated. Improvements based on the iterative use of data (iterative LLS and PLS imputation--ILLSimpute and IPLSimpute), the need to perform initial imputations (modified PLS and Helland PLS imputation--MPLSimpute and HPLSimpute) and the type of relationships employed (KNNarray, LLSarray, HPLSarray and alternating PLS--APLSimpute) are proposed. Overall, it is shown that data set properties (type of experiment, missing rate and pattern) affect the data similarity structure, therefore influencing the methods' performance. LLSimpute and ILLSimpute are preferable in the presence of data with a stronger similarity structure (TS and MIX experiments), whereas PLS-based methods (MPLSimpute, IPLSimpute and APLSimpute) are preferable when estimating NTS missing data.