Comparative genomics usually involves managing the functional aspects of genomes, by simply comparing gene-by-gene functions. Following this approach, Mushegian and Koonin proposed a hypothetical minimal genome, Minimal Gene Set (MGS), aiming for a possible oldest ancestor genome. They obtained MGS by comparing the genomes of two simple bacteria and eliminating duplicated or functionally identical genes. The authors raised the fundamental question of whether a hypothetical organism possessing MGS is able to live or not. We attacked this viability problem specifying in silico the metabolic pathways of the MGS-based prokaryote. We then performed a dynamic simulation of cellular metabolic activities in order to check whether the MGS-prokaryote reaches some equilibrium state and produces the necessary biomass. We assumed these two conditions to be necessary for a living organism. Our simulations clearly show that the MGS does not express an organism that is able to live. We then iteratively proceeded with functional replacements in order to obtain a genome composition that gives rise to equilibrium. We ruled out 76 of the original 254 genes in the MGS, because they resulted in duplication from a functional point of view. We also added seven genes not present in the MGS. These genes encode for enzymes involved in critical nodes of the metabolic network. These modifications led to a genome composed of 187 elements expressing a virtually living organism, Virtual Cell (ViCe), that exhibits homeostatic capabilities and produces biomass. Moreover, the steady-state distribution of the concentrations of virtual metabolites that resulted was similar to that experimentally measured in bacteria. We conclude then that ViCe is able to “live in silico.”
References
[1]
Glass J, Assad-Garcia N, Alperovich N, Yooseph S, Lewis MR, et al. (2006) Essential genes of a minimal bacterium. Proc Natl Acad Sci 103: 425–430.
[2]
Itaya M (1995) An estimation of minimal genome size required for life. FEBS Lett 362: 257–260.
[3]
Mushegian AR, Koonin EV (1996) A minimal gene set for cellular life derived by comparison of complete bacterial genome. Proc Natl Acad Sci 93: 10268–10273.
[4]
Cho MK, Magnus D, Caplan AL, McGee D, Ethics of Genomics Group (1999) Ethical considerations in synthesizing a minimal genome. Science 286: 2087–2090.
[5]
Huynen H (2000) Constructing a minimal genome. Trends Genetics 16: 116–116.
[6]
Priami C, Regev A, Silverman W, Shapiro E (2001) Application of a stochastic passing-name calculus to representation and simulation of molecular processes. Inf Proc Lett 80: 25–31.
[7]
Regev A, Silverman W, Shapiro E (2001) Representation and simulation of biochemical processes using the π-calculus process algebra. Proceedings of the Pacific Symposium on Biocomputing (PSB2001). Big Island, Hawaii, United States: pp. 459–470.
[8]
Curti M, Degano P, Baldari CT (2003) Causal π-calculus for biochemical modeling. Lect Notes Comput Sci 2003: 21–33.
[9]
Curti M, Degano P, Priami C, Baldari CT (2004) Modeling biochemical pathways through enhanced π-calculus. Theor Comput Sci 325: 111–140.
[10]
Milner R, Parrow J, Walker D (1992) A calculus of mobile processes (i and ii). Information Computation 100: 1–72.
[11]
Degano P, Priami C (2001) Enhanced Operational Semantics: A tool for describing and analysing concurrent systems. ACM Comput Surveys 33: 135–176.
[12]
Hucka M, Finney A, Sauro HM, Bolouri H, Doyle JC, et al. (2003) The systems biology markup language (sbml): A medium for representation and exchange of biochemical models. Bioinformatics 19: 524–531.
[13]
Paun G (2001) Membrane computing. An introduction. Berlin: Springer-Verlag.
[14]
Reddy VN, Mavrovouniotis ML, Liebman MN (1993) Petri Net representations in metabolic pathways. In: Proceedings of the Fourteenth International Conference on Intelligent Systems in Molecular Biology; July 1993; Bethesda, Maryland, United States. Association for the Advancement of Artificial Intelligence. ISMB 1: 328–336.
[15]
Prandi D, Priami C, Quaglia P (2005) Process calculi in a biological context. Bull EATCS 85: 53–69.
[16]
Chiarugi D, Curti M, Degano P, Marangoni R (2005) ViCe: A VIrtual CEll. Lect Notes Comput Sci 3082: 207–220.
[17]
Wandrey C, Buchholz A, Takors R (2001) Quantification of intracellular metabolites in k12 using liquid chromatographic-electrospray ionization tandem mass spectrometric techniques. Anal Biochem 295: 129–137.