17Increased technological methods have enabled the investigation of biology at nanoscale levels. 18 Nevertheless, such systems necessitate the use of computational methods to comprehend the complex 19 interactions occurring. Traditionally, dynamics of metabolic systems are described by ordinary differential 20 equations producing a deterministic result which neglects the intrinsic heterogeneity of biological systems. 21 More recently, stochastic modeling approaches have gained popularity with the capacity to provide more 22 realistic outcomes. Yet, solving stochastic algorithms tend to be computationally intensive processes. 23 Employing the queueing theory, an approach commonly used to evaluate telecommunication networks, 24 reduces the computational power required to generate simulated results, while simultaneously reducing 25 expansion of errors inherent to classical deterministic approaches. Herein, we present the application of 26 queueing theory to efficiently simulate stochastic metabolic networks. For the current model, we utilize 27 glycolysis to demonstrate the power of the proposed modeling methods, and we describe simulation and 28 pharmacological inhibition in glycolysis to further exemplify modeling capabilities.
30
Author Summary
31Computational biology is increasingly used to understand biological occurances and complex 32 dynamics. Biological modeling, in general, aims to represent a biological system with computational 33 approaches, as realistically and accurate as current methods allow. Metabolomics and metabolic systems 34 have emerged as an important aspect of cellular biology, allowing a more sentive view for understanding 35 the complex interactions occurring intracellularly as a result of normal or perturbed (or diseased) states. To 36 understand metabolic changes, many researchers have commonly used Ordianary Differential Equations to 37 produce in silico models of the in vitro system of interest. While these have been beneficial to date, 38 continuing to advance computational methods of analyzing such systems is of interest. Stochastic models 39 that include randomness have been known to produce more reaslistic results, yet the difficulty and intesive 40 time component urges additional methods and techniques to be developed. In the present research, we 3 41 propose using queueing networks as a technique to model complex metabolic systems, doing such with a 42 model of glycolysis, a core metabolic pathway. 43 44 58 are a representation of reality, aiming to accurately represent the system of study. Inclusion of all cellular 59 components indirectly or directly involved are considered far too complex to model. Consequently, 60 simplifications and assumptions must be made and often the perceived non-pivotal details, such as 61 stochasticity, omitted. Nevertheless, the accuracy and competence of the model is dependent on these 62 assumptions and simplifications. 63 Many approaches may be taken to model the dynamics of metabolic systems; importantly, the 64 categoriza...