A whole-cell computational model predicts phenotype from genotype.

Understanding how complex phenotypes arise from individual molecules and their interactions is a primary challenge in biology that computational approaches are poised to tackle. We report a whole-cell computational model of the life cycle of the human pathogen Mycoplasma genitalium that includes all of its molecular components and their interactions. An integrative approach to modeling that combines diverse mathematics enabled the simultaneous inclusion of fundamentally different cellular processes and experimental measurements. Our whole-cell model accounts for all annotated gene functions and was validated against a broad range of data. The model provides insights into many previously unobserved cellular behaviors, including in vivo rates of protein-DNA association and an inverse relationship between the durations of DNA replication initiation and replication. In addition, experimental analysis directed by model predictions identified previously undetected kinetic parameters and biological functions. We conclude that comprehensive whole-cell models can be used to facilitate biological discovery.

Accelerated discovery via a whole-cell model.

To test the promise of whole-cell modeling to facilitate scientific inquiry, we compared growth rates simulated in a whole-cell model with experimental measurements for all viable single-gene disruption Mycoplasma genitalium strains. Discrepancies between simulations and experiments led to predictions about kinetic parameters of specific enzymes that we subsequently validated. These findings represent, to our knowledge, the first application of whole-cell modeling to accelerate biological discovery.

WholeCellKB: model organism databases for comprehensive whole-cell models.

Whole-cell models promise to greatly facilitate the analysis of complex biological behaviors. Whole-cell model development requires comprehensive model organism databases. WholeCellKB (http://wholecellkb.stanford.edu) is an open-source web-based software program for constructing model organism databases. WholeCellKB provides an extensive and fully customizable data model that fully describes individual species including the structure and function of each gene, protein, reaction and pathway. We used WholeCellKB to create WholeCellKB-MG, a comprehensive database of the Gram-positive bacterium Mycoplasma genitalium using over 900 sources. WholeCellKB-MG is extensively cross-referenced to existing resources including BioCyc, KEGG and UniProt. WholeCellKB-MG is freely accessible through a web-based user interface as well as through a RESTful web service.

A forward-genetic screen and dynamic analysis of lambda phage host-dependencies reveals an extensive interaction network and a new anti-viral strategy.

Latently infecting viruses are an important class of virus that plays a key role in viral evolution and human health. Here we report a genome-scale forward-genetics screen for host-dependencies of the latently-infecting bacteriophage lambda. This screen identified 57 Escherichia coli (E. coli) genes--over half of which have not been previously associated with infection--that when knocked out inhibited lambda phage's ability to replicate. Our results demonstrate a highly integrated network between lambda and its host, in striking contrast to the results from a similar screen using the lytic-only infecting T7 virus. We then measured the growth of E. coli under normal and infected conditions, using wild-type and knockout strains deficient in one of the identified host genes, and found that genes from the same pathway often exhibited similar growth dynamics. This observation, combined with further computational and experimental analysis, led us to identify a previously unannotated gene, yneJ, as a novel regulator of lamB gene expression. A surprising result of this work was the identification of two highly conserved pathways involved in tRNA thiolation-one pathway is required for efficient lambda replication, while the other has anti-viral properties inhibiting lambda replication. Based on our data, it appears that 2-thiouridine modification of tRNAGlu, tRNAGln, and tRNALys is particularly important for the efficient production of infectious lambda phage particles.

Crystal structure, stability and in vitro RNAi activity of oligoribonucleotides containing the ribo-difluorotoluyl nucleotide: insights into substrate requirements by the human RISC Ago2 enzyme.

Short interfering RNA (siRNA) duplexes are currently being evaluated as antisense agents for gene silencing. Chemical modification of siRNAs is widely expected to be required for therapeutic applications in order to improve delivery, biostability and pharmacokinetic properties. Beyond potential improvements in the efficacy of oligoribonucleotides, chemical modification may also provide insight into the mechanism of mRNA downregulation mediated by the RNA-protein effector complexes (RNA-induced silencing complex or RISC). We have studied the in vitro activity in HeLa cells of siRNA duplexes against firefly luciferase with substitutions in the guide strand of U for the apolar ribo-2,4-difluorotoluyl nucleotide (rF) [Xia, J. et al. (2006) ACS Chem. Biol., 1, 176-183] as well as of C for rF. Whereas an internal rF:A pair adjacent to the Ago2 ('slicer' enzyme) cleavage site did not affect silencing relative to the native siRNA duplex, the rF:G pair and other mismatches such as A:G or A:A were not tolerated. The crystal structure at atomic resolution determined for an RNA dodecamer duplex with rF opposite G manifests only minor deviations between the geometries of rF:G and the native U:G wobble pair. This is in contrast to the previously found, significant deviations between the geometries of rF:A and U:A pairs. Comparison between the structures of the RNA duplex containing rF:G and a new structure of an RNA with A:G mismatches with the structures of standard Watson-Crick pairs in canonical duplex RNA leads to the conclusion that local widening of the duplex formed by the siRNA guide strand and the targeted region of mRNA is the most likely reason for the intolerance of human Ago2 (hAgo2), the RISC endonuclease, toward internal mismatch pairs involving native or chemically modified RNA. Contrary to the influence of shape, the thermodynamic stabilities of siRNA duplexes with single rF:A, A:A, G:A or C:A (instead of U:A) or rF:G pairs (instead of C:G) show no obvious correlation with their activities. However, incorporation of three rF:A pairs into an siRNA duplex leads to loss of activity. Our structural and stability data also shed light on the role of organic fluorine as a hydrogen bond acceptor. Accordingly, UV melting (T(M)) data, osmotic stress measurements, X-ray crystallography at atomic resolution and the results of semi-empirical calculations are all consistent with the existence of weak hydrogen bonds between fluorine and the H-N1(G) amino group in rF:G pairs of the investigated RNA dodecamers.

Mutator/Hypermutable fetal/juvenile metakaryotic stem cells and human colorectal carcinogenesis.

Adult age-specific colorectal cancer incidence rates increase exponentially from maturity, reach a maximum, then decline in extreme old age. Armitage and Doll (1) postulated that the exponential increase resulted from "n" mutations occurring throughout adult life in normal "cells at risk" that initiated the growth of a preneoplastic colony in which subsequent "m" mutations promoted one of the preneoplastic "cells at risk" to form a lethal neoplasia. We have reported cytologic evidence that these "cells at risk" are fetal/juvenile organogenic, then preneoplastic metakaryotic stem cells. Metakaryotic cells display stem-like behaviors of both symmetric and asymmetric nuclear divisions and peculiarities such as bell shaped nuclei and amitotic nuclear fission that distinguish them from embryonic, eukaryotic stem cells. Analyses of mutant colony sizes and numbers in adult lung epithelia supported the inferences that the metakaryotic organogenic stem cells are constitutively mutator/hypermutable and that their contributions to cancer initiation are limited to the fetal/juvenile period. We have amended the two-stage model of Armitage and Doll and incorporated these several inferences in a computer program CancerFit v.5.0. We compared the expectations of the amended model to adult (15-104 years) age-specific colon cancer rates for European-American males born 1890-99 and observed remarkable concordance. When estimates of normal colonic fetal/juvenile APC and OAT gene mutation rates (∼2-5 × 10(-5) per stem cell doubling) and preneoplastic colonic gene loss rates (∼8 × 10(-3)) were applied, the model was in accordance only for the values of n = 2 and m = 4 or 5.

A noisy paracrine signal determines the cellular NF-kappaB response to lipopolysaccharide.

Nearly identical cells can exhibit substantially different responses to the same stimulus. We monitored the nuclear localization dynamics of nuclear factor kappaB (NF-kappaB) in single cells stimulated with tumor necrosis factor-alpha (TNF-alpha) and lipopolysaccharide (LPS). Cells stimulated with TNF-alpha have quantitative differences in NF-kappaB nuclear localization, whereas LPS-stimulated cells can be clustered into transient or persistent responders, representing two qualitatively different groups based on the NF-kappaB response. These distinct behaviors can be linked to a secondary paracrine signal secreted at low concentrations, such that not all cells undergo a second round of NF-kappaB activation. From our single-cell data, we built a computational model that captures cell variability, as well as population behaviors. Our findings show that mammalian cells can create "noisy" environments to produce diversified responses to stimuli.