Bioelectric gene and reaction networks: computational modelling of genetic, biochemical and bioelectrical dynamics in pattern regulation

Gene regulatory networks (GRNs) describe interactions between gene products and transcription factors that control gene expression. In combination with reaction–diffusion models, GRNs have enhanced comprehension of biological pattern formation. However, although it is well known that biological systems exploit an interplay of genetic and physical mechanisms, instructive factors such as transmembrane potential (V mem ) have not been integrated into full GRN models. Here we extend regulatory networks to include bioelectric signalling, developing a novel synthesis: the bioelectricity-integrated gene and reaction (BIGR) network. Usingin silico simulations, we highlight the capacity forV mem to alter steady-state concentrations of key signalling molecules inside and out of cells. We characterize fundamental feedbacks whereV mem both controls, and is in turn regulated by, biochemical signals and thereby demonstrateV mem homeostatic control,V mem memory andV mem controlled state switching. BIGR networks demonstrating hysteresis are identified as a mechanisms through which more complex patterns of stableV mem spots and stripes, along with correlated concentration patterns, can spontaneously emerge. As further proof of principle, we present and analyse a BIGR network model that mechanistically explains key aspects of the remarkable regenerative powers of creatures such as planarian flatworms. The functional properties of BIGR networks generate the first testable, quantitative hypotheses for biophysical mechanisms underlying the stability and adaptive regulation of anatomical bioelectric pattern.