Although the crystal structure of Vibrio harveyi luciferase has been elucidated, the binding sites for the flavin mononucleotide and fatty aldehyde substrates are still unknown. The determined location of the phosphate-binding site close to Arg 107 on the ␣ subunit of luciferase is supported here by point mutagenesis. This information, together with previous structure-activity data for the length of the linker connecting the phosphate group to the isoalloxazine ring represent important characteristics of the luciferasebound conformation of the flavin mononucleotide. A model of the luciferase-flavin complex is developed here using flexible docking supplemented by these structural constraints. The location of the phosphate moiety was used as the anchor in a flexible docking procedure performed by conformation search by using the Monte Carlo minimization approach. The resulting databases of energy-ranked feasible conformations of the luciferase complexes with flavin mononucleotide, -phosphopentylflavin, -phosphobutylflavin, and -phosphopropylflavin were filtered according to the structure-activity profile of these analogs. A unique model was sought not only on energetic criteria but also on the geometric requirement that the isoalloxazine ring of the active flavin analogs must assume a common orientation in the luciferase-binding site, an orientation that is also inaccessible to the inactive flavin analog. The resulting model of the bacterial luciferase-flavin mononucleotide complex is consistent with the experimental data available in the literature. Specifically, the isoalloxazine ring of the flavin mononucleotide interacts with the Ala 74-Ala 75 cis-peptide bond as well as with the Cys 106 side chain in the ␣ subunit of luciferase. The model of the binary complex reveals a distinct cavity suitable for aldehyde binding adjacent to the isoalloxazine ring and flanked by other key residues (His 44 and Trp 250) implicated in the active site.Keywords: Bacterial luciferase; flavin mononucleotide; flexible docking; structural constraints Bacterial luciferase, a heterodimeric protein constituted of two homologous subunits (␣ and ), catalyzes the oxidation of reduced flavin mononucleotide and fatty aldehyde in a process that reduces molecular oxygen to water and releases energy in the form of light (∼490 nm). Both the ␣ (LuxA) and  (LuxB) subunits fold into (/␣) 8 barrels in the crystal structures of bacterial luciferase from Vibrio harveyi (Fisher et al. 1995(Fisher et al. , 1996. The ␣ subunit is primarily responsible for kinetic properties; however, the presence of the  subunit is essential for high catalytic efficiency (Cline and Hastings 1972;Li et al. 1993). The active site thus is believed to be formed by residues in the ␣ subunit with the phosphate moiety of flavin mononucleotide anchored at an electron dense inorganic phosphate site detected on LuxA (Fisher et al. 1995). Another critical factor controlling the biolumi-