Emergent Threshold Dynamics in Unified Metabolic-Regulatory Models: Achieving Spontaneous Homeostasis Through Bidirectional ATP-Pathway Coupling

Abstract

Background: For decades, computational biology has failed to create unified models where metabolic state and regulatory control are bidirectionally coupled: metabolic models optimize flux but cannot represent dynamic regulation, while regulatory models treat ATP as a fixed parameter rather than a dynamic variable affected by pathway activity. This fundamental limitation prevents computational recapitulation of emergent threshold behaviors-spontaneous homeostasis, adaptive reorganization, pathway switching-observed in living organisms. The challenge requires formalisms where (1) metabolic state governs regulatory decisions AND (2) regulatory choices consume metabolic resources, producing emergent dynamics from feedback rather than programming. Methods: We introduce Signal Hierarchical Petri Nets, extending Hybrid Petri Nets with bidirectional metabolic-regulatory coupling through energy-dependent layer organization. Unlike classical approaches, ATP is simultaneously a regulatory signal (governing pathway availability through quantitative thresholds) and a material substrate (consumed by pathway activity). When ATP depletes below 1000 M, high-cost pathways automatically become unavailable; pathway activity consuming ATP creates feedback affecting subsequent pathway accessibility. This bidirectional coupling enables emergent threshold behaviors impossible in classical formalisms. We demonstrate the paradigm through macrocyclic peptide transport across 53 metabolic conditions, where drug accumulation depends on ATP-governed pathway reorganization. Results: The formalism produces three emergent behaviors never achieved in unified metabolic-regulatory models. (1) Spontaneous homeostasis without programming: Despite 113-fold permeability variation from N-methylation, ATP-replete cells maintain constant drug accumulation (CV=0.066%)-homeostatic compensation emerges from ATP-consumption feedback, not explicit control logic. (2) Threshold-triggered reorganization: ATP depletion to 300 M triggers 8533-fold active-to-passive transport shifts with paradoxical 141% accumulation increase from efflux collapse. (3) Tissue-specific dynamics from identical parameters: Tumor hypoxia (ATP=1200 M) versus normal tissue (ATP=5000 M) produces 6.62-fold selectivity differences from differential pathway accessibility-same model, different emergent outcomes. Computational predictions achieve r=0.911 correlation with experimental cyclosporin permeability (n=32). Conclusions: Signal Hierarchical Petri Nets represent the first computational formalism achieving emergent threshold dynamics through bidirectional metabolic-regulatory coupling. The paradigm enables in silico recapitulation of adaptive cellular behaviors previously impossible to model, with applications extending beyond drug transport to any biological system where metabolic state governs regulatory reorganization: cancer metabolism, ischemia, synthetic biology, and aging research.