Cells execute their functions through dynamic operations of biological networks. Dynamic networks delineate the operation of biological networks in terms of temporal changes of abundances or activities of nodes (proteins and RNAs), as well as formation of new edges and disappearance of existing edges over time. Global genomic and proteomic technologies can be used to decode dynamic networks. However, using these experimental methods, it is still challenging to identify temporal transition of nodes and edges. Thus, several computational methods for estimating dynamic topological and functional characteristics of networks have been introduced. In this review, we summarize concepts and applications of these computational methods for inferring dynamic networks and further summarize methods for estimating spatial transition of biological networks.

MOTIVATION: Time-evolving differential protein-protein interaction (PPI) networks are essential to understand serial activation of differentially regulated (up- or downregulated) cellular processes (DRPs) and their interplays over time. Despite developments in the network inference, current methods are still limited in identifying temporal transition of structures of PPI networks, DRPs associated with the structural transition and the interplays among the DRPs over time. RESULTS: Here, we present a probabilistic model for estimating Time-Evolving differential PPI networks with MultiPle Information (TEMPI). This model describes probabilistic relationships among network structures, time-course gene expression data and Gene Ontology biological processes (GOBPs). By maximizing the likelihood of the probabilistic model, TEMPI estimates jointly the time-evolving differential PPI networks (TDNs) describing temporal transition of PPI network structures together with serial activation of DRPs associated with transiting networks. This joint estimation enables us to interpret the TDNs in terms of temporal transition of the DRPs. To demonstrate the utility of TEMPI, we applied it to two time-course datasets. TEMPI identified the TDNs that correctly delineated temporal transition of DRPs and time-dependent associations between the DRPs. These TDNs provide hypotheses for mechanisms underlying serial activation of key DRPs and their temporal associations. AVAILABILITY AND IMPLEMENTATION: Source code and sample data files are available at http://sbm.postech.ac.kr/tempi/sources.zip. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

MOTIVATION: Systems biology attempts to describe complex systems behaviors in terms of dynamic operations of biological networks. However, there is lack of tools that can effectively decode complex network dynamics over multiple conditions. RESULTS: We present principal network analysis (PNA) that can automatically capture major dynamic activation patterns over multiple conditions and then generate protein and metabolic subnetworks for the captured patterns. We first demonstrated the utility of this method by applying it to a synthetic dataset. The results showed that PNA correctly captured the subnetworks representing dynamics in the data. We further applied PNA to two time-course gene expression profiles collected from (i) MCF7 cells after treatments of HRG at multiple doses and (ii) brain samples of four strains of mice infected with two prion strains. The resulting subnetworks and their interactions revealed network dynamics associated with HRG dose-dependent regulation of cell proliferation and differentiation and early PrPSc accumulation during prion infection. AVAILABILITY: The web-based software is available at: http://sbm.postech.ac.kr/pna.

SUMMARY: We describe an integrative software platform, Prequips, for comparative proteomics-based systems biology analysis that: (i) integrates all information generated from mass spectrometry (MS)-based proteomics as well as from basic proteomics data analysis tools, (ii) visualizes such information for various proteomic analyses via graphical interfaces and (iii) links peptide and protein abundances to external tools often used in systems biology studies. AVAILABILITY: http://prequips.sourceforge.net

UNLABELLED: MS-BID (MS Biomarker Discovery Platform) is an integrative computational pipeline for biomarker discovery using LC-MS-based comparative proteomic analysis. This platform consists of several computational tools for: (i) detecting peptides in the collected patterns; (ii) matching detected peptides across a number of LC-MS datasets and (iii) selecting discriminatory peptides between classes of samples. AVAILABILITY: MS-BID source codes, binaries and documentations are freely available under LGPL from http://tools.proteomecenter.org/msBID.php.

MOTIVATIONS: Tissue engineering constitutes an important field with its potential of addressing the current shortage in organ availability. To successfully develop tissue-engineered organs, it is crucial to understand how to maintain the cells under conditions that maximize their ability to perform their physiological roles, regardless of the environment, whether the cells are part of an extracorporeal system, such as the bioartificial liver assist device, or an implantable tissue-engineered device. Our goals are to (1) provide insight into how cells will behave when confronted with changes in its environment and (2) determine the optimal environmental factors to achieve a desired level of cellular function. RESULTS: Diverse sets of environmental factors were used to systematically perturb the metabolic behavior associated with pre-conditioning and plasma supplementation. To probe metabolic state of hepatocytes, metabolic flux analysis was used to obtain the metabolic profile. We applied a multi-block partial least square (MPLS) model to relate environmental factors and fluxes to levels of intracellular lipids and urea synthesis. The MPLS model identified: (1) the most influential environmental factors and (2) how the metabolic pathways are altered by these factors. Finally, we inverted the MPLS model to determine the concentrations and types of environmental factors required to obtain the most economical solution for achieving optimal levels of cellular function for practical situations.

MOTIVATIONS: Classification of biological samples for diagnostic purposes is a difficult task because of the many decisions involved on the number, type and functional manipulations of the input variables. This study presents a generally applicable strategy for systematic formulation of optimal diagnostic indexes. To this end, we develop a novel set of computational tools by integrating regression optimization, stepwise variable selection and cross-validation algorithms. RESULTS: The proposed discrimination methodology was applied to plasma and tissue (liver) metabolic profiling data describing the time progression of liver dysfunction in a rat model of acute hepatic failure generated by d-galactosamine (GalN) injection. From the plasma data, our methodology identified seven (out of a total of 23) metabolites, and the corresponding transform functions, as the best inputs to the optimal diagnostic index. This index showed better time resolution and increased noise robustness compared with an existing metabolic index, Fischer's BCAA/AAA molar ratio, as well as indexes generated using other commonly used discriminant analysis tools. Comparison of plasma and liver indexes found two consensus metabolites, lactate and glucose, which implicate glycolysis and/or gluconeogenesis in mediating the metabolic effects of GalN.

MOTIVATION: Transcriptional profiling using microarrays can reveal important information about cellular and tissue expression phenotypes, but these measurements are costly and time consuming. Additionally, tissue sample availability poses further constraints on the number of arrays that can be analyzed in connection with a particular disease or state of interest. It is therefore important to provide a method for the determination of the minimum number of microarrays required to separate, with statistical reliability, distinct disease states or other physiological differences. RESULTS: Power analysis was applied to estimate the minimum sample size required for two-class and multi-class discrimination. The power analysis algorithm calculates the appropriate sample size for discrimination of phenotypic subtypes in a reduced dimensional space obtained by Fisher discriminant analysis (FDA). This approach was tested by applying the algorithm to existing data sets for estimation of the minimum sample size required for drawing certain conclusions on multi-class distinction with statistical reliability. It was confirmed that when the minimum number of samples estimated from power analysis is used, group means in the FDA discrimination space are statistically different. CONTACT: [email protected]

MOTIVATION: The increasing use of DNA microarrays to probe cell physiology requires methods for visualizing different expression phenotypes and explicitly connecting individual genes to discriminating expression features. Such methods should be robust and maintain biological interpretability. RESULTS: We propose a method for the mapping of the physiological state of cells and tissues from multidimensional expression data such as those obtained with DNA microarrays. The method uses Fisher discriminant analysis to create a linear projection of gene expression measurements that maximizes the separation of different sample classes. Relative to other typical classification methods, this method provides insights into the discriminating characteristics of expression measurements in terms of the contribution of individual genes to the definition of distinct physiological states. This projection method also facilitates visualization of classification results in a reduced dimensional space. Examples from four different cases demonstrate the ability of the method to produce well-separated groups in the projection space and to identify important genes for defining physiological states. The method can be augmented to also include data from the proteomic and metabolic phenotypes and can be useful in disease diagnosis, drug screening and bioprocessing applications.