We present a continuous modelling framework for simulating the dynamics of metabolic-regulatory networks (MRNs), designed to overcome the scalability limitations of traditional hybrid models. Hybrid approaches, often based on Boolean logic to represent regulatory interactions, become computationally intractable as the number of regulatory proteins increases, due to an exponential growth in discrete modes and transitions. To address this, our framework replaces discrete logic with smooth Hill functions, enabling the approximation of switch-like regulatory behaviour without introducing combinatorial complexity. This continuous formulation maintains the biological interpretability of hybrid models while greatly enhancing computational efficiency. Parameter estimation, a common bottleneck in continuous models, is simplified in our approach by requiring fewer kinetic parameters than typical hybrid models. We further employ sparse-based system identification, a data-driven technique that efficiently infers network dynamics by selecting a minimal set of nonlinear terms. This method avoids exhaustive search procedures and yields interpretable kinetic models. Applied to MRNs, our framework demonstrates the ability to capture essential regulatory mechanisms with reduced complexity and improved scalability.