The bioluminescence emission from fireflies is an astounding tool to mark and view cells. However, the bioluminescent mechanism is not completely deciphered, limiting the comprehension of key processes. We use a theoretical approach to study for the first time the arrival of a dioxygen molecule inside the fireflies protein and one path of the formation of the dioxetanone ring, the high-energy intermediate precursor of the bioluminescence. To describe this reaction step, a joint approach combining classical molecular dynamics (MD) simulations and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations is used. The formation of the dioxetanone ring has been studied for both singlet and triplet states with the help of MS-CASPT2 calculations. We also emphasize the role played by the proteinic environment in the formation of the dioxetanone ring. The results obtained shed some light on an important reaction step and give new insights concerning the bioluminescence in fireflies.