Morphogenesis is a precise and robust dynamic process during metazoan embryogenesis, consisting of both cell proliferation and cell migration. Despite the fact that much is known about specific regulations at molecular level, how cell proliferation and migration together drive the morphogenesis at cellular and organismic levels is not well understood. Using Caenorhabditis elegans as the model animal, we present a phase field model to compute early embryonic morphogenesis within a confined eggshell. With physical information about cell division obtained from three-dimensional time-lapse cellular imaging experiments, the model can precisely reproduce the early morphogenesis process as seen in vivo, including time evolution of location and morphology of each cell. Furthermore, the model can be used to reveal key cell-cell attractions critical to the development of C. elegans embryo. Our work demonstrates how genetic programming and physical forces collaborate to drive morphogenesis and provides a predictive model to decipher the underlying mechanism.

Metabolic networks undergo gene expression regulation in response to external nutrient signals. In microbes, the synthesis of enzymes that are used to transport and catabolize less preferred carbon sources is repressed in the presence of a preferred carbon source. For most microbes, glucose is a preferred carbon source, and it has long been believed that as long as glucose is present in the environment, the expression of genes related to the metabolism of alternative carbon sources is shut down, due to catabolite repression. However, recent studies have shown that the induction of the galactose (GAL) metabolic network does not solely depend on the exhaustion of glucose. Instead, the GAL genes respond to the external concentration ratio of galactose to glucose, a phenomenon of unknown mechanism that we termed ratio-sensing. Using mathematical modeling, we found that ratio-sensing is a general phenomenon that can arise from competition between two carbon sources for shared transporters, between transcription factors for binding to communal regulatory sequences of the target genes, or a combination of the aforementioned two levels of competition. We analyzed how the parameters describing the competitive interaction influenced ratio-sensing behaviors in each scenario and found that the concatenation of both layers of signal integration could expand the dynamical range of ratio-sensing. Finally, we investigated the influence of circuit topology on ratio-sensing and found that incorporating negative auto-regulation and/or coherent feedforward loop motifs to the basic signal integration unit could tune the sensitivity of the response to the external nutrient signals. Our study not only deepened our understanding of how ratio-sensing is achieved in yeast GAL metabolic regulation, but also elucidated design principles for ratio-sensing signal processing that can be used in other biological settings, such as being introduced into circuit designs for synthetic biology applications.

MOTIVATION: Epistatic Miniarray Profiles (EMAP) has enabled the mapping of large-scale genetic interaction networks; however, the quantitative information gained from EMAP cannot be fully exploited since the data are usually interpreted as a discrete network based on an arbitrary hard threshold. To address such limitations, we adopted a mixture modeling procedure to construct a probabilistic genetic interaction network and then implemented a Bayesian approach to identify densely interacting modules in the probabilistic network. RESULTS: Mixture modeling has been demonstrated as an effective soft-threshold technique of EMAP measures. The Bayesian approach was applied to an EMAP dataset studying the early secretory pathway in Saccharomyces cerevisiae. Twenty-seven modules were identified, and 14 of those were enriched by gold standard functional gene sets. We also conducted a detailed comparison with state-of-the-art algorithms, hierarchical cluster and Markov clustering. The experimental results show that the Bayesian approach outperforms others in efficiently recovering biologically significant modules.

Drug molecules not only interact with specific targets, but also alter the state and function of the associated biological network. How to design drugs and evaluate their functions at the systems level becomes a key issue in highly efficient and low-side-effect drug design. The arachidonic acid metabolic network is the network that produces inflammatory mediators, in which several enzymes, including cyclooxygenase-2 (COX-2), have been used as targets for anti-inflammatory drugs. However, neither the century-old nonsteriodal anti-inflammatory drugs nor the recently revocatory Vioxx have provided completely successful anti-inflammatory treatment. To gain more insights into the anti-inflammatory drug design, the authors have studied the dynamic properties of arachidonic acid (AA) metabolic network in human polymorphous leukocytes. Metabolic flux, exogenous AA effects, and drug efficacy have been analyzed using ordinary differential equations. The flux balance in the AA network was found to be important for efficient and safe drug design. When only the 5-lipoxygenase (5-LOX) inhibitor was used, the flux of the COX-2 pathway was increased significantly, showing that a single functional inhibitor cannot effectively control the production of inflammatory mediators. When both COX-2 and 5-LOX were blocked, the production of inflammatory mediators could be completely shut off. The authors have also investigated the differences between a dual-functional COX-2 and 5-LOX inhibitor and a mixture of these two types of inhibitors. Their work provides an example for the integration of systems biology and drug discovery.

MOTIVATION: The use of gene microchips has enabled a rapid accumulation of gene-expression data. One of the major challenges of analyzing this data is the diversity, in both size and signal strength, of the various modules in the gene regulatory networks of organisms. RESULTS: Based on the iterative signature algorithm [Bergmann,S., Ihmels,J. and Barkai,N. (2002) Phys. Rev. E 67, 031902], we present an algorithm-the progressive iterative signature algorithm (PISA)-that, by sequentially eliminating modules, allows unsupervised identification of both large and small regulatory modules. We applied PISA to a large set of yeast gene-expression data, and, using the Gene Ontology database as a reference, found that the algorithm is much better able to identify regulatory modules than methods based on high-throughput transcription-factor binding experiments or on comparative genomics.