MOTIVATION: Maintenance of the self-renewal state in human embryonic stem cells (hESCs) is the foremost critical step for regenerative therapy applications. The insulin-mediated PI3K/AKT pathway is well appreciated as being the central pathway supporting hESC self-renewal; however, the regulatory interactions in the pathway that maintain cell state are not yet known. Identification of these regulatory pathway components will be critical for designing targeted interventions to facilitate a completely defined platform for hESC propagation and differentiation. Here, we have developed a systems analysis approach to identify regulatory components that control PI3K/AKT pathway in self-renewing hESCs. RESULTS: A detailed mathematical model was adopted to explain the complex regulatory interactions in the PI3K/AKT pathway. We evaluated globally sensitive processes of the pathway in a computationally efficient manner by replacing the detailed model by a surrogate meta-model. Our mathematical analysis, supported by experimental validation, reveals that negative regulators of the molecules IRS1 and PIP3 primarily govern the steady state of the pathway in hESCs. Among the regulators, negative feedback via IRS1 reduces the sensitivity of various reactions associated with direct trunk of the pathway and also constraints the propagation of parameter uncertainty to the levels of post receptor signaling molecules. Furthermore, our results suggest that inhibition of negative feedback can significantly increase p-AKT levels and thereby, better support hESC self-renewal. Our integrated mathematical modeling and experimental workflow demonstrates the significant advantage of computationally efficient meta-model approaches to detect sensitive targets from signaling pathways. AVAILABILITY AND IMPLEMENTATION: FORTRAN codes for the PI3K/AKT pathway and the RS-HDMR implementation are available from the authors upon request.

MOTIVATION: Primary purpose of modeling gene regulatory networks for developmental process is to reveal pathways governing the cellular differentiation to specific phenotypes. Knowledge of differentiation network will enable generation of desired cell fates by careful alteration of the governing network by adequate manipulation of cellular environment. RESULTS: We have developed a novel integer programming-based approach to reconstruct the underlying regulatory architecture of differentiating embryonic stem cells from discrete temporal gene expression data. The network reconstruction problem is formulated using inherent features of biological networks: (i) that of cascade architecture which enables treatment of the entire complex network as a set of interconnected modules and (ii) that of sparsity of interconnection between the transcription factors. The developed framework is applied to the system of embryonic stem cells differentiating towards pancreatic lineage. Experimentally determined expression profile dynamics of relevant transcription factors serve as the input to the network identification algorithm. The developed formulation accurately captures many of the known regulatory modes involved in pancreatic differentiation. The predictive capacity of the model is tested by simulating an in silico potential pathway of subsequent differentiation. The predicted pathway is experimentally verified by concurrent differentiation experiments. Experimental results agree well with model predictions, thereby illustrating the predictive accuracy of the proposed algorithm. CONTACT: [email protected] SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.