Back

Spatiotemporal 4D Whole-cell Modeling of a Minimal Autotroph Reveals Central Carbon Metabolism Regulated Locally by Protein Megacomplexes via Post-translational Modifications under Light Disturbance

Johnson, C. G.; Chan, A.; Rozum, J.; George, A.; Parvate, A. D.; Kim, D. N.; Feng, S.; Bohutskyi, P.; Johnson, Z.; Sadler, N.; Garcia, M.; Li, X.; Trejo, J.; Wu, R.; Sineath, W.; Anderson, L. N.; Evans, J. E.; Mehta, A. P.; Qian, W.-J.; Luthey-Schulten, Z.; Cheung, M. S.

2026-01-23 systems biology
10.64898/2026.01.21.700212 bioRxiv
Show abstract

Photosynthetic microorganisms rely on multiple pathways in central carbon metabolism to adapt to fluctuating light and energy availability across diel cycles. Mechanistic insight into the regulatory dynamics of this adaptation requires integrating processes spanning disparate timescales, from rapid redox-dependent post-translational modifications (PTMs) to slower changes in protein expression and metabolic pathway usage. To address this complexity beyond genome-based inference and traditional modeling, we develop a whole-cell four-dimensional (3D + time) model of the marine cyanobacterium Prochlorococcus marinus MED4 that explicitly represents the spatial organization of enzymatic and molecular processes in central carbon metabolism under light perturbation. We employ a perturbation-based research design to experimentally generate time-series, multi-omics measurements that provide molecular descriptors and cryo-ET images as constraints for this dynamic 4D framework. The integration of experiments and modeling across defined light regimes enables quantitative validation of system-level responses and forecasting under distinct light disturbances. We test the hypothesis that light-dependent redox PTMs regulating the structural assembly of a protein megacomplex, the "dark complex," modulate metabolic flux at a conserved regulatory node of the Calvin-Benson cycle (CBC) in cyanobacteria. Our model shows that subcellular spatial organization buffers rapid light-induced changes in thylakoid reaction rates, which are followed by redox-PTM-mediated sequestration or release of CBC enzymes in the dark complex, ultimately impacting carbon fixation dynamics within carboxysomes. Comparison with an equivalently parameterized well-mixed stochastic model demonstrates that post-translational regulation not only buffers transcriptional noise and diffusion-driven fluctuations but also stabilizes phenotypic outcomes, underscoring the importance of spatial heterogeneity in phenotypic robustness. This ability to probe adaptive, spatiotemporally resolved mechanisms in photosynthetic machinery and central carbon metabolism addresses a critical gap in genotype-to-phenotype inference and expands modeling and design capabilities for understudied or genetically intractable autotrophs such as P. marinus MED4. Significance StatementThis work advances 4D whole-cell modeling by presenting the first spatiotemporal simulation of a photosynthetic autotroph using the Lattice Microbes platform. Using Prochlorococcus marinus MED4, we show that subcellular spatial organization of organelles, diffusion constraints, and redox regulation collectively shape central carbon metabolism across orders of magnitude in space and time. Through a perturbation-based strategy that generates multi-omics data sets over time, we construct and validate a spatially and temporally resolved model of MED4, constrained by high-resolution (10 nm) cryo-electron tomography. Our results highlight the importance of localized biochemical reactions and redox-dependent post-translational modification of enzymes in regulating carbon fixation in a noisy environment under light disturbance. This study establishes a spatiotemporal, whole-cell physiology modeling framework as a transformative tool for uncovering multiscale regulatory responses to environmental gradients.

Matching journals

The top 4 journals account for 50% of the predicted probability mass.

50% of probability mass above

"Similar papers" are the closest papers from that journal in the model's embedding space. They show what the match is built on, but the ranking comes mostly from a classifier over the whole training set, not from these examples alone.