A critical review on modelling formalisms and simulation tools in computational biosystems

Integration of different kinds of biological processes is an ultimate goal for whole-cell modelling. We briefly review modelling formalisms that have been used in Systems Biology and identify the criteria that must be addressed by an integrating framework capable of modelling, analysing and simulating different biological networks. Aware that no formalism can fit all purposes we realize Petri nets as a suitable model for Metabolic Engineering and take a deeper perspective on the role of this formalism as an integrating framework for regulatory and metabolic networks.Research supported by PhD grant SFRH/BD/35215/2007 from the Fundacao para a Ciencia e a Tecnologia (FCT) and the MIT-Portugal program

TbUNC119 and Its Binding Protein Complex Are Essential for Propagation, Motility, and Morphogenesis of Trypanosoma brucei Procyclic Form Cells

Flagellum-mediated motility of Trypanosoma brucei is considered to be essential for the parasite to complete stage development in the tsetse fly vector, while the mechanism by which flagellum-mediated motility is controlled are not fully understood. We thus compared T. brucei whole gene products (amino acid sequence) with Caenorhabditis elegans UNC (uncoordinated) proteins, in order to find uncharacterized motility-related T. brucei genes. Through in silico analysis, we found 88 gene products which were highly similar to C. elegans UNC proteins and categorized them as TbCEUN (T. brucei gene products which have high similarity to C. elegans UNC proteins). Approximately two thirds of the 88 TbCEUN gene products were kinesin-related molecules. A gene product highly similar to C. elegans UNC119 protein was designated as TbUNC119. RNAi-mediated depletion of TbUNC119 showed no apparent phenotype. However, knock-down analysis of both TbUNC119 and its binding protein (TbUNC119BP) which was found by yeast two-hybrid analysis showed characteristic phenotypes, including reduced motility, morphological change (extended cell shape), and cellular apoptosis. Based on the observed phenotypes, possible function of the TbUNC119 and TbUNC119BP is discussed

Systems Biology Modeling Reveals a Possible Mechanism of the Tumor Cell Death upon Oncogene Inactivation in EGFR Addicted Cancers

Despite many evidences supporting the concept of “oncogene addiction” and many hypotheses rationalizing it, there is still a lack of detailed understanding to the precise molecular mechanism underlying oncogene addiction. In this account, we developed a mathematic model of epidermal growth factor receptor (EGFR) associated signaling network, which involves EGFR-driving proliferation/pro-survival signaling pathways Ras/extracellular-signal-regulated kinase (ERK) and phosphoinositol-3 kinase (PI3K)/AKT, and pro-apoptotic signaling pathway apoptosis signal-regulating kinase 1 (ASK1)/p38. In the setting of sustained EGFR activation, the simulation results show a persistent high level of proliferation/pro-survival effectors phospho-ERK and phospho-AKT, and a basal level of pro-apoptotic effector phospho-p38. The potential of p38 activation (apoptotic potential) due to the elevated level of reactive oxygen species (ROS) is largely suppressed by the negative crosstalk between PI3K/AKT and ASK1/p38 pathways. Upon acute EGFR inactivation, the survival signals decay rapidly, followed by a fast increase of the apoptotic signal due to the release of apoptotic potential. Overall, our systems biology modeling together with experimental validations reveals that inhibition of survival signals and concomitant release of apoptotic potential jointly contribute to the tumor cell death following the inhibition of addicted oncogene in EGFR addicted cancers

Hematopoietic prostaglandin D synthase (HPGDS): A high stability, Val187Ile isoenzyme common among African Americans and its relationship to risk for colorectal cancer

Intestinal tumors in ApcMin/+ mice are suppressed by over-production of HPGDS, which is a glutathione transferase that forms prostaglandin D2 (PGD2). We characterized naturally occurring HPGDS isoenzymes, to see if HPGDS variation is associated with human colorectal cancer risk. We used DNA heteroduplex analysis and sequencing to identify HPGDS variants among healthy individuals. HPGDS isoenzymes were produced in bacteria, and their catalytic activities were tested. To determine in vivo effects, we conducted pooled case-control analyses to assess whether there is an association of the isoenzyme with colorectal cancer. Roughly 8% of African Americans and 2% of Caucasians had a highly stable Val187lle isoenzyme (with isoleucine instead of valine at position 187). At 37 °C, the wild-type enzyme lost 15% of its activity in one hour, whereas the Val187Ile form remained >95% active. At 50 °C, the half life of native HPGDS was 9 minutes, compared to 42 minutes for Val187Ile. The odds ratio for colorectal cancer among African Americans with Val187Ile was 1.10 (95% CI, 0.75–1.62; 533 cases, 795 controls). Thus, the Val187Ile HPGDS isoenzyme common among African Americans is not associated with colorectal cancer risk. Other approaches will be needed to establish a role for HPGDS in occurrence of human intestinal tumors, as indicated by a mouse model

Principled simulation of cell proliferation dynamics using the CoSMoS approach

A collaboration between cancer biologists and academic software engineers has been exploring the development of an agent-based simulator to inform and support work on the dynamics of cell proliferation in the study of prostate disorders. The research has influenced and been informed by the CoSMoS project. This paper presents the simulation project (which is not yet complete). We reflect on the reality of following CoSMoS principles; we describe the domain exploration and show how software modelling approaches (here, Petri nets, state diagrams) can be used to express both biological and software models. We explore fitness for purpose and consider ways to present a fitness argument. We consider issues in choosing simulation media and mapping from domain models through to code. The implementation emphasis is on traceability to support reuse and extension of the simulator, as well as demonstrable fitness for purpose. Initial work on calibration is presented. We discuss the calibration results, that both support and challenge the design and assumptions captured in the domain modelling and development activities