Large-scale fermentations (»100 m³) often encounter concentration gradients which may significantly affect microbial activities and production performance. Reliably investigating such scenarios in silico would allow to optimize bioproduction. But related simulations are very rare in particular for large bubble columns. Here, we pioneer the spatially resolved investigation of a 600 m³ bubble column operating for Escherichia coli based l-phenylalanine fed-batch production. Microbial kinetics are derived from experimental data. Advanced Euler-Lagrange (EL) computational fluid dynamics (CFD) simulations are applied to track individual bubble dynamics that result from a recently developed bubble breakage model. Thereon, the complex nonlinear characteristics of hydrodynamics, mass transfer, and microbial activities are simulated for large scale and compared with real data. As a key characteristic, zones for upriser, downcomer, and circulation cells were identified that dominate mixing and mass transfer. This results in complex gradients of glucose, dissolved oxygen, and microbial rates dividing the bioreactor into sections. Consequently, alternate feed designs are evaluated splitting real feed rates in two feeds at different locations. The opposite reversed installation of feed spots and spargers improved the product synthesis by 6.24% while alternate scenarios increased the growth rate by 11.05%. The results demonstrate how sophisticated, spatially resolved simulations of hydrodynamics, mass transfer, and microbial kinetics help to optimize bioreactors in silico.