This is the first report of a computational study of the bioluminescence of ligand-bound Photinus pyralis luciferase. A time-dependent PBE0/molecular mechanics approach was used to study the interaction between excited-state oxyluciferin (Keto-(-1)) and neighboring active site molecules. The results of these calculations demonstrated that the most important intermolecular interactions are: blue-shifting ionic interactions, red-shifting π-π stacking, and red/blue shifting hydrogen bonding. Subsequent molecular dynamics simulations further supported these conclusions.