Stability and robustness of a generalized pump-leak model for epithelial cell and lumen volume regulation
Epithelial cells regulate ion concentrations and volume through coordinated membrane pumps, ion channels, and paracellular pathways which can be modeled by classical single-compartment pump-leak equations (PLEs). Many epithelial functions, however, depend on the interaction between a cell and an enclosed luminal space, a geometry that cannot be captured by classical PLEs. To address this, we develop a two-compartment model consisting of an intracellular compartment coupled to a luminal compartment through the apical membrane, with both compartments interfacing an infinite extracellular bath and connected to it through the basolateral membrane and a paracellular pathway. Building on the five-dimensional single-cell PLEs, we formulate a ten-dimensional PLE system for this geometry and derive analytical equilibria and steady-state formulas for both the passive system and the Na+/K+-ATPase (NKA) driven active system. We characterize how these equilibria depend on physiologically relevant parameters, analyze local stability across wide parameter ranges, and apply global sensitivity and robustness methods to identify the principal determinants of ion and volume homeostasis. The model reveals fundamental differences between basolateral and apical placement of the NKA, including the onset of luminal volume blow-up when apical potassium recycling is insufficient. More broadly, this framework provides a mathematically tractable and physiologically grounded foundation for studying epithelial transport and for predicting conditions under which pump localization and conductance changes lead to stable function or pathological lumen expansion.