Cross-evaluation of E. coli's operon structures via a whole-cell model suggests alternative cellular benefits for low- versus high-expressing operons.

Many bacteria use operons to coregulate genes, but it remains unclear how operons benefit bacteria. We integrated E. coli's 788 polycistronic operons and 1,231 transcription units into an existing whole-cell model and found inconsistencies between the proposed operon structures and the RNA-seq read counts that the model was parameterized from. We resolved these inconsistencies through iterative, model-guided corrections to both datasets, including the correction of RNA-seq counts of short genes that were misreported as zero by existing alignment algorithms. The resulting model suggested two main modes by which operons benefit bacteria. For 86% of low-expression operons, adding operons increased the co-expression probabilities of their constituent proteins, whereas for 92% of high-expression operons, adding operons resulted in more stable expression ratios between the proteins. These simulations underscored the need for further experimental work on how operons reduce noise and synchronize both the expression timing and the quantity of constituent genes. A record of this paper's transparent peer review process is included in the supplemental information.

A forecast for large-scale, predictive biology: Lessons from meteorology.

Quantitative systems biology, in which predictive mathematical models are constructed to guide the design of experiments and predict experimental outcomes, is at an exciting transition point, where the foundational scientific principles are becoming established, but the impact is not yet global. The next steps necessary for mathematical modeling to transform biological research and applications, in the same way it has already transformed other fields, is not completely clear. The purpose of this perspective is to forecast possible answers to this question-what needs to happen next-by drawing on the experience gained in another field, specifically meteorology. We review here a number of lessons learned in weather prediction that are directly relevant to biological systems modeling, and that we believe can enable the same kinds of global impact in our field as atmospheric modeling makes today.

Oncogenic mutant RAS signaling activity is rescaled by the ERK/MAPK pathway.

Activating mutations in RAS are present in ~ 30% of human tumors, and the resulting aberrations in ERK/MAPK signaling play a central role in oncogenesis. However, the form of these signaling changes is uncertain, with activating RAS mutants linked to both increased and decreased ERK activation in vivo. Rationally targeting the kinase activity of this pathway requires clarification of the quantitative effects of RAS mutations. Here, we use live-cell imaging in cells expressing only one RAS isoform to quantify ERK activity with a new level of accuracy. We find that despite large differences in their biochemical activity, mutant KRAS isoforms within cells have similar ranges of ERK output. We identify roles for pathway-level effects, including variation in feedback strength and feedforward modulation of phosphatase activity, that act to rescale pathway sensitivity, ultimately resisting changes in the dynamic range of ERK activity while preserving responsiveness to growth factor stimuli. Our results reconcile seemingly inconsistent reports within the literature and imply that the signaling changes induced by RAS mutations early in oncogenesis are subtle.

Systems-Level Properties of EGFR-RAS-ERK Signaling Amplify Local Signals to Generate Dynamic Gene Expression Heterogeneity.

Intratumoral heterogeneity is associated with aggressive tumor behavior, therapy resistance, and poor patient outcomes. Such heterogeneity is thought to be dynamic, shifting over periods of minutes to hours in response to signaling inputs from the tumor microenvironment. However, models of this process have been inferred from indirect or post-hoc measurements of cell state, leaving the temporal details of signaling-driven heterogeneity undefined. Here, we developed a live-cell model system in which microenvironment-driven signaling dynamics can be directly observed and linked to variation in gene expression. Our analysis reveals that paracrine signaling between two cell types is sufficient to drive continual diversification of gene expression programs. This diversification emerges from systems-level properties of the EGFR-RAS-ERK signaling cascade, including intracellular amplification of amphiregulin-mediated paracrine signals and differential kinetic filtering by target genes including Fra-1, c-Myc, and Egr1. Our data enable more precise modeling of paracrine-driven transcriptional variation as a generator of gene expression heterogeneity. A record of this paper's transparent peer review process is included in the Supplemental Information.

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.

Neuronal ER-plasma membrane junctions organized by Kv2-VAP pairing recruit Nir proteins and affect phosphoinositide homeostasis.

The association of plasma membrane (PM)-localized voltage-gated potassium (Kv2) channels with endoplasmic reticulum (ER)-localized vesicle-associated membrane protein-associated proteins VAPA and VAPB defines ER-PM junctions in mammalian brain neurons. Here, we used proteomics to identify proteins associated with Kv2/VAP-containing ER-PM junctions. We found that the VAP-interacting membrane-associated phosphatidylinositol (PtdIns) transfer proteins PYK2 N-terminal domain-interacting receptor 2 (Nir2) and Nir3 specifically associate with Kv2.1 complexes. When coexpressed with Kv2.1 and VAPA in HEK293T cells, Nir2 colocalized with cell-surface-conducting and -nonconducting Kv2.1 isoforms. This was enhanced by muscarinic-mediated PtdIns(4,5)P2 hydrolysis, leading to dynamic recruitment of Nir2 to Kv2.1 clusters. In cultured rat hippocampal neurons, exogenously expressed Nir2 did not strongly colocalize with Kv2.1, unless exogenous VAPA was also expressed, supporting the notion that VAPA mediates the spatial association of Kv2.1 and Nir2. Immunolabeling signals of endogenous Kv2.1, Nir2, and VAP puncta were spatially correlated in cultured neurons. Fluorescence-recovery-after-photobleaching experiments revealed that Kv2.1, VAPA, and Nir2 have comparable turnover rates at ER-PM junctions, suggesting that they form complexes at these sites. Exogenous Kv2.1 expression in HEK293T cells resulted in significant differences in the kinetics of PtdIns(4,5)P2 recovery following repetitive muscarinic stimulation, with no apparent impact on resting PtdIns(4,5)P2 or PtdIns(4)P levels. Finally, the brains of Kv2.1-knockout mice had altered composition of PtdIns lipids, suggesting a crucial role for native Kv2.1-containing ER-PM junctions in regulating PtdIns lipid metabolism in brain neurons. These results suggest that ER-PM junctions formed by Kv2 channel-VAP pairing regulate PtdIns lipid homeostasis via VAP-associated PtdIns transfer proteins.

Linear Integration of ERK Activity Predominates over Persistence Detection in Fra-1 Regulation.

ERK signaling regulates the expression of target genes, but it is unclear how ERK activity dynamics are interpreted. Here, we investigate this question using simultaneous, live, single-cell imaging of two ERK activity reporters and expression of Fra-1, a target gene controlling epithelial cell identity. We find that Fra-1 is expressed in proportion to the amplitude and duration of ERK activity. In contrast to previous "persistence detector" and "selective filter" models in which Fra-1 expression only occurs when ERK activity persists beyond a threshold duration, our observations demonstrate that the network regulating Fra-1 expression integrates total ERK activity and responds to it linearly. However, exploration of a generalized mathematical model of the Fra-1 coherent feedforward loop demonstrates that it can perform either linear integration or persistence detection, depending on the basal mRNA production rate and protein production delays. Our data indicate that significant basal expression and short delays cause Fra-1 to respond linearly to integrated ERK activity.

Akt regulation of glycolysis mediates bioenergetic stability in epithelial cells.

Cells use multiple feedback controls to regulate metabolism in response to nutrient and signaling inputs. However, feedback creates the potential for unstable network responses. We examined how concentrations of key metabolites and signaling pathways interact to maintain homeostasis in proliferating human cells, using fluorescent reporters for AMPK activity, Akt activity, and cytosolic NADH/NAD+ redox. Across various conditions, including glycolytic or mitochondrial inhibition or cell proliferation, we observed distinct patterns of AMPK activity, including both stable adaptation and highly dynamic behaviors such as periodic oscillations and irregular fluctuations that indicate a failure to reach a steady state. Fluctuations in AMPK activity, Akt activity, and cytosolic NADH/NAD+ redox state were temporally linked in individual cells adapting to metabolic perturbations. By monitoring single-cell dynamics in each of these contexts, we identified PI3K/Akt regulation of glycolysis as a multifaceted modulator of single-cell metabolic dynamics that is required to maintain metabolic stability in proliferating cells.

Single-Cell Imaging of ERK Signaling Using Fluorescent Biosensors.

Single-cell analysis of the mitogen-activated protein kinase (MAPK) extracellular signal-regulated kinase (ERK) provides a means to perform highly detailed kinetic studies, assess heterogeneity between cells, and distinguish the subcellular localization of ERK activity. We describe here the methods needed to perform such measurements in a cell type of the investigator's choosing. We discuss the selection of appropriate reporters and provide detailed methods for stably introducing reporters, collecting live-cell data, and automatically extracting quantitative information from individual cells.

Cell division orientation in animals.

Cell division orientation during animal development can serve to correctly organize and shape tissues, create cellular diversity or both. The underlying cellular mechanism is regulated spindle orientation. Depending on the developmental context, extrinsic signals or intrinsic cues control the correct orientation of the mitotic spindle. Cell geometry has been known to be another determinant of spindle orientation and recent results have shed new light on the link between cellular shape and cell division orientation. The importance of controlling spindle orientation is manifested in neurodevelopmental defects such as microcephaly, tumor initiation as well as defects in tissue architecture and cell fate misspecification. Here, we summarize the role of oriented cell division during animal development and also outline the cellular and molecular mechanisms in selected invertebrate and vertebrate systems.