In many bacteria the flavoenzyme thymidylate synthase ThyX produces the DNA nucleotide deoxythymidine monophosphate from dUMP, using methylenetetrahydrofolate as carbon donor and NADPH as hydride donor. Because all three substrates bind in close proximity to the catalytic flavin adenine dinucleotide group, substantial flexibility of the ThyX active site has been hypothesized. Using femtosecond time-resolved fluorescence spectroscopy, we have studied the conformational heterogeneity and the conformational interconversion dynamics in real time in ThyX from the hyperthermophilic bacterium Thermotoga maritima. The dynamics of electron transfer to excited flavin adenine dinucleotide from a neighboring tyrosine residue are used as a sensitive probe of the functional dynamics of the active site. The fluorescence decay spanned a full three orders of magnitude, demonstrating a very wide range of conformations. In particular, at physiological temperatures, multiple angstrom cofactor-residue displacements occur on the picoseconds timescale. These experimental findings are supported by molecular dynamics simulations. Binding of the dUMP substrate abolishes this flexibility and stabilizes the active site in a configuration where dUMP closely interacts with the flavin cofactor and very efficiently quenches fluorescence itself. Our results indicate a dynamic selected-fit mechanism where binding of the first substrate dUMP at high temperature stabilizes the enzyme in a configuration favorable for interaction with the second substrate NADPH, and more generally have important implications for the role of active site flexibility in enzymes interacting with multiple poly-atom substrates and products. Moreover, our data provide the basis for exploring the effect of inhibitor molecules on the active site dynamics of ThyX and other multisubstrate flavoenzymes.protein dynamics | flavoprotein | ultrafast fluorescence spectroscopy | quenching C onfigurational flexibility is essential for enzyme function during catalysis. Binding of one or more substrates, accommodation of the transition state where the actual reaction takes place, relaxation to the product state, and release of the product (s) is possible because different configurations of the enzyme are continuously sampled, by thermal or reaction-driven motions, on timescales ranging from femtoseconds to microseconds. The configurational space sampled in a certain time range will depend on the local protein flexibility/energy landscape and the temperature (1). During these configurational changes, distances between constituents of the enzyme complex change. It has been recognized that the fastest (femtosecond to picosecond) localized motions exist alongside the slower motions that occur on the (typically millisecond) timescale of catalysis (2-4).Various experimental techniques allow monitoring changes in interactions resulting from configurational changes. Long-range micro/millisecond domain motions in flexible proteins have been studied by NMR (5) and FRET techniques (6, 7). Molecula...