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.

Deep Learning Automates the Quantitative Analysis of Individual Cells in Live-Cell Imaging Experiments.

Live-cell imaging has opened an exciting window into the role cellular heterogeneity plays in dynamic, living systems. A major critical challenge for this class of experiments is the problem of image segmentation, or determining which parts of a microscope image correspond to which individual cells. Current approaches require many hours of manual curation and depend on approaches that are difficult to share between labs. They are also unable to robustly segment the cytoplasms of mammalian cells. Here, we show that deep convolutional neural networks, a supervised machine learning method, can solve this challenge for multiple cell types across the domains of life. We demonstrate that this approach can robustly segment fluorescent images of cell nuclei as well as phase images of the cytoplasms of individual bacterial and mammalian cells from phase contrast images without the need for a fluorescent cytoplasmic marker. These networks also enable the simultaneous segmentation and identification of different mammalian cell types grown in co-culture. A quantitative comparison with prior methods demonstrates that convolutional neural networks have improved accuracy and lead to a significant reduction in curation time. We relay our experience in designing and optimizing deep convolutional neural networks for this task and outline several design rules that we found led to robust performance. We conclude that deep convolutional neural networks are an accurate method that require less curation time, are generalizable to a multiplicity of cell types, from bacteria to mammalian cells, and expand live-cell imaging capabilities to include multi-cell type systems.

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.

Competing pathways control host resistance to virus via tRNA modification and programmed ribosomal frameshifting.

Viral infection depends on a complex interplay between host and viral factors. Here, we link host susceptibility to viral infection to a network encompassing sulfur metabolism, tRNA modification, competitive binding, and programmed ribosomal frameshifting (PRF). We first demonstrate that the iron-sulfur cluster biosynthesis pathway in Escherichia coli exerts a protective effect during lambda phage infection, while a tRNA thiolation pathway enhances viral infection. We show that tRNA(Lys) uridine 34 modification inhibits PRF to influence the ratio of lambda phage proteins gpG and gpGT. Computational modeling and experiments suggest that the role of the iron-sulfur cluster biosynthesis pathway in infection is indirect, via competitive binding of the shared sulfur donor IscS. Based on the universality of many key components of this network, in both the host and the virus, we anticipate that these findings may have broad relevance to understanding other infections, including viral infection of humans.

Simultaneous cross-evaluation of heterogeneous E. coli datasets via mechanistic simulation.

The extensive heterogeneity of biological data poses challenges to analysis and interpretation. Construction of a large-scale mechanistic model of Escherichia coli enabled us to integrate and cross-evaluate a massive, heterogeneous dataset based on measurements reported by various groups over decades. We identified inconsistencies with functional consequences across the data, including that the total output of the ribosomes and RNA polymerases described by data are not sufficient for a cell to reproduce measured doubling times, that measured metabolic parameters are neither fully compatible with each other nor with overall growth, and that essential proteins are absent during the cell cycle-and the cell is robust to this absence. Finally, considering these data as a whole leads to successful predictions of new experimental outcomes, in this case protein half-lives.

Live-cell measurements of kinase activity in single cells using translocation reporters.

Although kinases are important regulators of many cellular processes, measuring their activity in live cells remains challenging. We have developed kinase translocation reporters (KTRs), which enable multiplexed measurements of the dynamics of kinase activity at a single-cell level. These KTRs are composed of an engineered construct in which a kinase substrate is fused to a bipartite nuclear localization signal (bNLS) and nuclear export signal (NES), as well as to a fluorescent protein for microscopy-based detection of its localization. The negative charge introduced by phosphorylation of the substrate is used to directly modulate nuclear import and export, thereby regulating the reporter's distribution between the cytoplasm and nucleus. The relative cytoplasmic versus nuclear fluorescence of the KTR construct (the C/N ratio) is used as a proxy for the kinase activity in living, single cells. Multiple KTRs can be studied in the same cell by fusing them to different fluorescent proteins. Here, we present a protocol to execute and analyze live-cell microscopy experiments using KTRs. We describe strategies for development of new KTRs and procedures for lentiviral expression of KTRs in a cell line of choice. Cells are then plated in a 96-well plate, from which multichannel fluorescent images are acquired with automated time-lapse microscopy. We provide detailed guidance for a computational analysis and parameterization pipeline. The entire procedure, from virus production to data analysis, can be completed in ∼10 d.

Measuring Signaling and RNA-Seq in the Same Cell Links Gene Expression to Dynamic Patterns of NF-κB Activation.

Signaling proteins display remarkable cell-to-cell heterogeneity in their dynamic responses to stimuli, but the consequences of this heterogeneity remain largely unknown. For instance, the contribution of the dynamics of the innate immune transcription factor nuclear factor κB (NF-κB) to gene expression output is disputed. Here we explore these questions by integrating live-cell imaging approaches with single-cell sequencing technologies. We used this approach to measure both the dynamics of lipopolysaccharide-induced NF-κB activation and the global transcriptional response in the same individual cell. Our results identify multiple, distinct cytokine expression patterns that are correlated with NF-κB activation dynamics, establishing a functional role for NF-κB dynamics in determining cellular phenotypes. Applications of this approach to other model systems and single-cell sequencing technologies have significant potential for discovery, as it is now possible to trace cellular behavior from the initial stimulus, through the signaling pathways, down to genome-wide changes in gene expression, all inside of a single cell.

The future of whole-cell modeling.

Integrated whole-cell modeling is poised to make a dramatic impact on molecular and systems biology, bioengineering, and medicine--once certain obstacles are overcome. From our group's experience building a whole-cell model of Mycoplasma genitalium, we identified several significant challenges to building models of more complex cells. Here we review and discuss these challenges in seven areas: first, experimental interrogation; second, data curation; third, model building and integration; fourth, accelerated computation; fifth, analysis and visualization; sixth, model validation; and seventh, collaboration and community development. Surmounting these challenges will require the cooperation of an interdisciplinary group of researchers to create increasingly sophisticated whole-cell models and make data, models, and simulations more accessible to the wider community.