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 (Vmem) 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. Using in silico simulations, we highlight the capacity for Vmem to alter steady-state concentrations of key signalling molecules inside and out of cells. We characterize fundamental feedbacks where Vmem both controls, and is in turn regulated by, biochemical signals and thereby demonstrate Vmem homeostatic control, Vmem memory and Vmem controlled state switching. BIGR networks demonstrating hysteresis are identified as a mechanisms through which more complex patterns of stable Vmem 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.