eng International Academy of Ecology and Environmental Sciences Network Biology 2220-8879 2022-12-1 12 4 142 159 2 article Md. Jasim Uddin 1 2 S. M. Sohel Rana 1 2 Department of Mathematics, University of Dhaka, Dhaka-1000, Bangladesh Glycolysis model has been considered by Caputo fractional derivative. We give the topological classifications of fixed points of this model. Then, we show analytically that under certain parametric conditions fractional order glycolysis model underlies a Neimark-Sacker (NS) bifurcation and Flip bifurctaion. By using central manifold and bifurcation theory, we confirm the existence and direction of both NS and Flip bifurcations. To reinforce our analytical findings, we perform numerical simulations that include bifurcations, phase portraits, periodic orbits, invariant closed cycles, abrupt emergence of chaos and abrupt elimination of chaos. At the end, OGY method is applied to eliminate chaotic trajectories of the system. http://www.iaees.org/publications/journals/nb/articles/2022-12(4)/chaotic-dynamics-of-a-fractional-order-glycolysis-model.pdf fractional order glycolysis model Flip and Neimark-Sacker (NS) bifurcations Maximum Lyapunov Exponents (MLEs) Fractal Dimensions (FDs) chaos control