Dynamics of the firefly luciferase-oxyluciferin complex in its electronic ground and excited states are studied using various theoretical approaches. By mimicking the physiological conditions with realistic models of the chromophore oxyluciferin, the enzyme luciferase, and solvating water molecules and by performing real time simulations with a molecular dynamics technique on the model surfaces, we reveal that the local chromophore-surrounding interaction patterns differ rather severely in the two states. Because of the presence of protein, the solvation dynamics of water around the chromophore is also peculiar and shows widely different time scales on the two terminal oxygen atoms. In addition, simulations of the emission with the quantum-mechanics/molecular-mechanics approach show a close relationship between the emission color variation and the environmental dynamics, mostly through electrostatic effects from the chromophore-surrounding interaction. We also discuss the importance of considering the time scales of the luminescence and the dynamics of the interaction. © 2011 American Chemical Society