To study the network dynamics of the riboflavin biosynthesis pathway and to identify potential bottlenecks in the system, an ordinary differential equation-based model was constructed using available literature data for production strains. The results confirmed that the RibA protein is rate limiting in the pathway. Under the conditions investigated, we determined a potential limiting order of the remaining enzymes under increased RibA concentration (> 0.102 mM) and therefore higher riboflavin production (> 0.045 mmol g (CDW) (-1) h(-1) and 0.0035 mM s(-1), respectively). The reductase activity of RibG and lumazine synthase (RibH) might be the next most limiting steps. The computational minimization of the enzyme concentrations of the pathway suggested the need for a greater RibH concentration (0.251 mM) compared with the other enzymes (RibG: 0.188 mM, RibB: 0.023 mM).