Cited by 1,300Gene Regulatory Network AnalysisMicrobial Metabolic Engineering and BioproductionBioinformatics and Genomic Networks

A Whole-Cell Computational Model Predicts Phenotype from Genotype

Open access
Karr, Jonathan R., Sanghvi, Jayodita C., Macklin, Derek N., Gutschow, Miriam V., Jacobs, Jared M., Bolival, Benjamin, Assad-Garcia, Nacyra, Glass, John I. +1 more·2012·Cell

Abstract

SUMMARY 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 rates. 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.

Integrity signals

Registry

Crossref · 10.1016/j.cell.2012.05.044

Cites

47 of 47 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Claims the paper makes

5 of 28 cited sources checked · 5 supported · 0 flagged

  • Findingcitation

    Predicted high-occupancy RNA polymerase regions correspond to highly transcribed rRNAs and tRNAs.

    “Consistent with a recent experimental study by Vora et al., 2009, the predicted high-occupancy RNA polymerase regions correspond to highly transcribed rRNAs and tRNAs.”
    Quote in the paperNot a substantial claim1 of 1 sources checked
    • cites Vora et al., 2009 source
  • Interpretationown data

    The results support the assertion that large-scale modeling can aid biological discovery.

    “These results support the assertion that large-scale modeling can be used to support biological discovery (Kitano, 2002; Brenner, 2010).”
    Quote in the paperSubstantial claim1 of 2 sources checked
    • cites Kitano, 2002 source
    • cites Brenner, 2010 source
  • Findingown data

    Model calculations matched the observed doubling time.

    “We found that the model calculations were consistent with the observed doubling time (Figure 2A and 2B), cellular chemical composition (Figure 2C), replication of major cell mass fractions (Figure 2D), and gene expression (R2 = 0.68; Figure S1A).”
    Quote in the paperSubstantial claim
  • Findingown data

    The authors built a whole-cell computational model of the M. genitalium life cycle that incorporates all molecular components and their interactions.

    “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.”
    Quote in the paperSubstantial claim
    • The authors state this as an assumption: The cellular process sub-models are approximately independent of one another on short time scales of less than one second.
  • Findingown data

    The model predicts an inverse relation between the length of DNA replication initiation and replication rate.

    “including in vivo rates of protein-DNA association and an inverse relationship between the durations of DNA replication initiation and replication rates.”
    Quote in the paperNot a substantial claim
  • Classified as probably not a substantial claim. 7
  • May combine two claims. 26
  • Borderline: one of the tests for a substantial claim is unclear. 1
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisecitation

    Later methods needing fewer parameters include Boolean network and constraint-based modeling.

    “Subsequently, alternative approaches were developed that require fewer parameters, including Boolean network modeling (Davidson et al., 2002) and constraint-based modeling (Orth et al., 2010; Thiele et al., 2009).”
    Quote in the paperNot a substantial claim2 of 3 sources checked
    • cites Davidson et al., 2002 source
    • cites Orth et al., 2010 source
    • cites Thiele et al., 2009 source
  • May combine two claims.

    Findingown data

    Predicted metabolite concentrations are within an order of magnitude of E. coli measurements for 100% of metabolites in one dataset and 70% in a newer high-throughput study.

    “the predicted metabolite concentrations are within an order of magnitude of concentrations measured in Escherichia coli for 100% of the metabolites in one compilation of data (Sundararaj et al., 2004) and for 70% in a more recent high-throughput study (Bennett et al., 2009; Figure 2F).”
    Quote in the paperSubstantial claim1 of 2 sources checked
    • cites Sundararaj et al., 2004 source
    • cites Bennett et al., 2009 source
  • Borderline on S3 (load-bearing), S4 (based): unclear from the text. May combine two claims.

    Findingnone given

    The model predicts over 30,000 protein-protein collisions per cell cycle on average, displacing 0.93 proteins per second.

    “Our model predicts that over 30,000 collisions occur on average per cell cycle, leading to the displacement of 0.93 proteins per second.”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Findingown data

    The most debilitating disruptions involved metabolic genes and prevented production of any major cell mass components.

    “The most debilitating disruptions involved metabolic genes and resulted in the inability to produce any of the major cell mass components.”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Findingown data

    Most collisions are caused by RNA polymerase (84%) and DNA polymerase (8%), and most commonly displace SMC proteins (70%) or single-stranded binding proteins (6%).

    “the majority of collisions are caused by RNA polymerase (84%) and DNA polymerase (8%), most commonly resulting in the displacement of Structural Maintenance of Chromosome (SMC) proteins (70%), or single-stranded binding proteins (6%) (Figure 3E and Table S2F).”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Findingown data

    The model predicts that 50% of the chromosome is bound by at least one protein within 6 min and 90% within 20 min.

    “The model further predicts that the chromosome is explored very rapidly, with 50% of the chromosome having been bound by at least one protein within the first 6 min of the cell cycle, and 90% within the first 20 min (Figure 3B).”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    Reducing Tdk kcat in the model lowered the predicted thyA strain growth rate but also affected wild type growth rate.

    “Reducing the Tdk kcat in the model did indeed reduce the predicted growth rate of the thyA strain, but it also affected the wild type growth rate (Figure 7F).”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    Growth rates of twelve single-gene disruption strains were measured experimentally; ten were correctly predicted viable and two incorrectly predicted nonviable.

    “we experimentally measured the growth rates of twelve single-gene disruption strains – ten of which were correctly predicted to be viable and two of which that were incorrectly predicted to be nonviable – for comparison to our models' predictions (Figure 7A).”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    Stochastic transcription and translation create variability in new DnaA monomer numbers over time and in DnaA binding and unbinding events.

    “we observed that the stochastic aspect of the transcription and translation sub-models creates variability in the number of new DnaA monomers produced over time, as well as the DnaA binding and unbinding events themselves.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    Predicted DNA polymerase chromosomal occupancy is significantly lower and biased toward terC.

    “In contrast, the predicted DNA polymerase chromosomal occupancy is significantly lower and biased toward the terC (see below for further discussion).”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Findingown data

    For Tdk and Pdp, original values came from distantly related organisms, while newly predicted kcat values agree with reports from more closely related species.

    “For Tdk and Pdp, these values corresponded to distantly related organisms; however, the newly predicted kcat values are consistent with reports from more closely related species (Figure 7H).”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Findingcitation

    The model predicts 90% of genes are expressed within the first 143 min, with transcription lagging RNA polymerase exploration because non-specific RNA polymerase-DNA interactions contribute to RNA polymerase diffusion.

    “On average, this results in expression of 90% of genes within the first 143 min (Figure 3C), with transcription lagging RNA polymerase exploration due to the significant contribution of non-specific RNA polymerase-DNA interactions to RNA polymerase diffusion (Harada et al., 1999).”
    Quote in the paperSubstantial claim0 of 1 sources checked
    • cites Harada et al., 1999 source
  • Classified as probably not a substantial claim.

    Study design and conductnone given

    The model rested on a detailed M. genitalium reconstruction assembled from more than 900 primary sources, reviews, books, and databases.

    “The whole-cell model was based on a detailed reconstruction of M. genitalium developed from over 900 primary sources, reviews, books, and databases.”
    Quote in the paperSubstantial claim
    • The authors state this as an assumption: The cellular process sub-models are approximately independent of one another on short time scales of less than one second.
  • May combine two claims.

    Study design and conductnone given

    The model synthesizes over 900 publications and contains more than 1,900 experimentally observed parameters.

    “Our model is based on a synthesis of over 900 publications and includes more than 1,900 experimentally observed parameters.”
    Quote in the paperSubstantial claim
    • The authors state this as an assumption: The cellular process sub-models are approximately independent of one another on short time scales of less than one second.
  • May combine two claims.

    Premisereasoning

    The assumptions behind Boolean and constraint-based methods do not hold for all cellular processes and conditions, so a whole-cell model based only on either is impractical.

    “However, the underlying assumptions of these methods do not apply to all cellular processes and conditions, and building a whole-cell model entirely based on either method is therefore impractical.”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing). May combine two claims.

    Premisecitation

    M. genitalium is a human urogenital parasite with 525 genes in its genome.

    “Here, we present a “whole-cell” model of the bacterium Mycoplasma genitalium, a human urogenital parasite whose genome contains 525 genes (Fraser et al., 1995).”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites Fraser et al., 1995 source
  • Classified as probably not a substantial claim.

    Premisenone given

    Before recent advances, knowledge of individual molecules and their interactions was insufficient to model any organism completely.

    “First, until recently not enough has been known about the individual molecules and their interactions to completely model any one organism.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Interpretationown data

    In the lpdA, deoD, and thyA cases, resolving model-experiment discrepancies led to further analysis that resolved the discrepancy and gave new insight into M. genitalium biology.

    “In each of these three cases (lpdA, deoD, thyA), identifying a discrepancy between model predictions and experimental measurements led to further analysis which resolved the discrepancy, and also provided novel insight into M. genitalium biology (Figure 7I).”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim: may not pass S4 (based). May combine two claims.

    Premisenone given

    The thyA gene product catalyzes dTMP production and can be complemented by the tdk gene product.

    “The thyA gene product catalyzes dTMP production and can be complemented by the tdk gene product.”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Premisecitation

    The lpdA gene product is part of the pyruvate dehydrogenase complex, which transfers electrons to NAD as part of the overall pyruvate dehydrogenase reaction.

    “The lpdA gene product is part of the pyruvate dehydrogenase complex, which catalyzes the transfer of electrons to NAD as a subset of the overall pyruvate dehydrogenase chemical reaction (de Kok et al., 1998).”
    Quote in the paperSubstantial claim0 of 1 sources checked
    • cites de Kok et al., 1998 source
  • May combine two claims.

    Interpretationown data

    Replication initially proceeds quickly because dNTPs are abundant; after the dNTP pool is exhausted, replication slows to the rate of dNTP synthesis.

    “Initially replication proceeds quickly, due to the free dNTP content in the cell (Figure 4B, bottom). When DNA polymerase initially binds to the replication origin, dNTPs are abundant and replication can proceed unimpeded. When the dNTP pool is exhausted, however, the rate of replication slows to the rate of dNTP synthesis.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Premisenone given

    DeoD converts deoxyadenosine to adenine and D-ribose-1-phosphate; the pdp gene product can also produce these products from deoxyuridine.

    “DeoD catalyzes the conversion of deoxyadenosine to adenine and D-ribose-1-phosphate; these products can also be produced by the pdp gene product from deoxyuridine.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Interpretationown data

    The whole-cell model can predict instantaneous chromosomal occupancy and the temporal dynamics and interactions of every DNA-binding protein at genomic scale and single-cell resolution.

    “the whole-cell model can predict both the instantaneous protein chromosomal occupancy as well as the temporal dynamics and interactions of every DNA-binding protein at the genomic scale at single-cell resolution.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Premisecitation

    The nox locus lies in a sub-operon containing two other pyruvate dehydrogenase genes and is coexpressed with pdhA.

    “Moreover, the nox locus falls in a sub-operon that contains two other pyruvate dehydrogenase genes and has been shown to be coexpressed with pdhA (Guell et al., 2009) (Figure 7D)”
    Quote in the paperSubstantial claim0 of 1 sources checked
    • cites Guell et al., 2009 source
  • Classified as probably not a substantial claim.

    Premisecitation

    Many M. genitalium genes are multi-functional.

    “Since previous studies have shown that many M. genitalium genes are multi-functional (Pollack et al., 2002; Cordwell et al., 1997)”
    Quote in the paperSubstantial claim0 of 2 sources checked
    • cites Pollack et al., 2002 source
    • cites Cordwell et al., 1997 source
  • May combine two claims.

    Interpretationown data

    Disrupting lpdA had a severe but non-critical effect on cell growth.

    “Taken together, the data suggested that disrupting the lpdA gene had a severe, but non-critical impact on cell growth.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Interpretationreasoning

    Cells needing extra time to initiate replication accumulate a large dNTP surplus, which speeds replication.

    “This occurs because cells that require extra time to initiate replication also build up a large dNTP surplus, leading to faster replication.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Interpretationreasoning

    The authors developed a comprehensive whole-cell model that accounts for all annotated M. genitalium gene functions and explains emergent behaviors through molecular interactions.

    “We have developed a comprehensive, whole-cell model that accounts for all of the annotated gene functions identified in M. genitalium and explains a variety of emergent behaviors in terms of molecular interactions.”
    Quote in the paperSubstantial claim
    • The authors state this as an assumption: The cellular process sub-models are approximately independent of one another on short time scales of less than one second.
  • Classified as probably not a substantial claim. May combine two claims.

    Premisenone given

    M. genitalium is difficult to study experimentally because of resistance to most antibiotics, lack of a chemically defined medium, and a cell size requiring advanced microscopy.

    “Resistance to most antibiotics, the lack of a chemically defined medium, and a cell size that requires advanced microscopy techniques for visualization, all greatly limit the range of experimental techniques available to study this organism.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Interpretationown data

    The model enables rapid identification of novel gene functions and specific cellular parameters.

    “Our model accurately recapitulates a broad set of experimental data, provides insight into several biological processes for which experimental assessment is not readily feasible, and enables the rapid identification of novel gene functions as well as specific cellular parameters.”
    Quote in the paperSubstantial claim
30 of 30
  • The paper relies on this without support.

    Biological systems are modular, so cells can be modeled by dividing them into functional processes, modeling each process independently on a short time scale, and integrating the sub-models at longer time scales.

    “Because biological systems are modular, cells can be modeled by (1) dividing cells into functional processes, (2) independently modeling each process on a short time scale, and (3) integrating process sub-models at longer time scales.”
  • The authors state this as an assumption.

    The cellular process sub-models are approximately independent of one another on short time scales of less than one second.

    “We began with the assumption that the sub-models are approximately independent on short time scales (less than one second in this work).”
  • The paper relies on this without support.

    Data and parameters obtained from other organisms are applicable enough to M. genitalium to build and validate the model.

    “As a result, much of the data used to build and validate the model was obtained from other organisms.”
3 of 3

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

  • Existing cell models are often rule-based and combine assumptions about the underlying biological mechanisms with parameters fit from observational data.

    “Existing cell models are often rule-based and combine assumptions about the underlying biological mechanisms with parameters fit from observational data.”

    From How to build the virtual cell with artificial intelligence: Priorities and opportunities (2024) · cites this paper as [8]

    The Introduction discusses Boolean and constraint-based cell models and their underlying assumptions. Results, under “Model training and parameter reconciliation,” reports more than 1,900 experimentally observed parameters in the authors’ model. The paper does not establish that existing cell models are often rule-based or that their parameters are generally fitted to observational data; most parameters in this model were implemented as originally reported. (LLM, conf=0.82)

    Experimental — accuracy not yet measured
  • Genome-wide models have been developed to represent other bacterial organisms, including Escherichia coli.

    “Since this pioneering work, genome-wide models have been developed to represent other bacterial organisms, including Escherichia coli.8,20–22”

    From How to build the virtual cell with artificial intelligence: Priorities and opportunities (2024) · cites this paper as [8]

    The Abstract identifies this paper as a whole-cell model of Mycoplasma genitalium, supporting its role as the pioneering work referenced in the sentence. The Discussion says that modeling E. coli would be needed for comprehensive validation; it does not report a developed genome-wide E. coli model or establish the later developments attributed to the grouped citations. Page numbers are unavailable in our copy. (LLM, conf=0.84)

    Experimental — accuracy not yet measured
  • The first whole-cell model was developed in 2012, representing all 482 genes and molecular functions known for Mycobacterium genitalium.

    “The first whole-cell model was developed in 2012, representing all 482 genes and molecular functions known for an organism: the bacteria Mycobacterium genitalium.8”

    From How to build the virtual cell with artificial intelligence: Priorities and opportunities (2024) · cites this paper as [8]

    The Experimental Procedures (Reconstruction; page unavailable) describe a whole-cell model of M. genitalium and functional annotation of each gene. But the Results (single-gene disruption phenotypes; page unavailable) describe 525 possible single-gene disruptions, not 482, and the organism is Mycoplasma genitalium, not Mycobacterium genitalium. (LLM, conf=0.90)

    Experimental — accuracy not yet measured
  • No one organism has a complete set of experimentally determined kinetic parameters, chemical composition, and physical characteristics, so the model presented here, as well as other whole-cell models, has depended on the transferability of quantities among related organisms.

    “Unfortunately, no one organism has a complete set of experimentally determined kinetic parameters, chemical composition, and physical characteristics, so the model we present here, as well as other whole-cell models, has depended on the transferability of quantities among related organisms.1”

    From Bringing the genetically minimal cell to life on a computer in 4D (2026) · cites this paper as [3]

    Split into 2 atoms — 1 confirmed, 1 partially supported, 0 disputed, 0 misattributed, 0 unverifiable. Expand the claim to see per-atom verdicts.

    Experimental — accuracy not yet measured
  • Existing cell models generally rely on explicitly defined mathematical or computational approaches, such as differential equations, stochastic simulations, or agent-based models.

    “They generally rely on explicitly defined mathematical or computational approaches, such as differential equations,7–9 stochastic simulations,10,11 or agent-based models.12,13”

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

Each contradiction between two passages listed here was confirmed by an independent check.

  • Shorter doubling time conflicts with the reported growth impairment

    Minor

    A doubling time 40% lower than wild type indicates faster growth, whereas the authors interpret the result as a severe impairment of growth. The direction of the reported change conflicts with that conclusion.

  • Tdk adjustment is attributed to different disruption strains

    Minor

    The stated hypothesis concerns the tdk disruption strain, but the next sentence says the adjustment confirmed it by reducing growth of the thyA strain. The strain named in the predicted outcome does not match the one in the reported outcome.

  • Sub-model time scale is both below and equal to one second

    Minor

    The results describe the short time scale used in this work as less than one second, while the methods say each process was modeled independently on a one-second time scale.

3 of 3

Cites

Flagged first, then by citations

Cited by

Sorted by citations

Authors
Karr, J. R., Sanghvi, J. C., Macklin, D. N., Gutschow, M. V., Jacobs, J. M., Bolival, B., Assad-Garcia, N., Glass, J. I., & Covert, M. W.
Affiliations
  • Stanford University (US)
  • J. Craig Venter Institute (US)
Corresponding author
Covert, Markus W.
Venue
Cell
Publisher
Elsevier BV
Published
2012-07
Volume
150
Issue
2
Pages
389-401
ISSN
0092-8674 (ISSN-L)
OpenAlex
W2081781433
Cited by
1,300
References
47
Citations by year
2026: 84 · 2025: 59 · 2024: 48 · 2023: 52 · 2022: 51 · 2021: 91 · 2020: 75 · 2019: 101 · 2018: 105 · 2017: 104 · 2016: 105 · 2015: 137 · 2014: 119 · 2013: 126 · 2012: 41
FWCI
54.33
Type
journal article
Language
en
Open access
Locations
  • Cell (publishedVersion)
  • PubMed (publishedVersion)
  • Europe PMC (PubMed Central) (acceptedVersion)
  • PubMed Central (acceptedVersion)
Funders
  • National Institutes of Health (DP1OD006413)
  • U.S. National Library of Medicine (DP1LM011510)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Gene Regulatory Network Analysis
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Microbial Metabolic Engineering and Bioproduction
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Bioinformatics and Genomic Networks
Keywords
Mycoplasma genitalium, genotype-phenotype correlation, life cycle modeling, DNA replication initiation, gene expression regulation, kinetic parameter estimation
MeSH terms
Bacterial Proteins, Cell Cycle, Computer Simulation, DNA-Binding Proteins, Models, Biological, Phenotype, Mycoplasma genitalium, Molecular Sequence Annotation

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Claims read by deepseek-v4.1-flash under definition version 3 (extract-v3) on 2026-09-28.

Classified by jev-1.13-20260917 (gold-guideline-v0.2) on 2026-09-28.

Cited sources are checked by gpt-6-sol.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review