The main focus of enzymology is on the enzyme rates, substrate structures, and reactivity, whereas the role of solvent dynamics in mediating the biological reaction is often left aside owing to its complex molecular behavior. We used integrated X-ray-and terahertz-based time-resolved spectroscopic tools to study proteinwater dynamics during proteolysis of collagen-like substrates by a matrix metalloproteinase. We show equilibration of structural kinetic transitions in the millisecond timescale during degradation of the two model substrates collagen and gelatin, which have different supersecondary structure and flexibility. Unexpectedly, the detected changes in collective enzyme-substrate-water-coupled motions persisted well beyond steady state for both substrates while displaying substrate-specific behaviors. Molecular dynamics simulations further showed that a hydration funnel (i.e., a gradient in retardation of hydrogen bond (HB) dynamics toward the active site) is substrate-dependent, exhibiting a steeper gradient for the more complex enzyme-collagen system. The long-lasting changes in protein-water dynamics reflect a collection of local energetic equilibrium states specifically formed during substrate conversion. Thus, the observed long-lasting water dynamics contribute to the net enzyme reactivity, impacting substrate binding, positional catalysis, and product release.solvation dynamics | enzyme catalysis | metalloenzymes I t has now been a century since the Michaelis-Menten (MM) steady-state theory provided a highly satisfactory description of the kinetic behavior of many enzymes (1, 2). The simplest kinetic mechanism with a single substrate assumes that an enzyme E combines with a substrate S to form an ES complex (known as the Michaelis complex), which undergoes an irreversible reaction to form the product P and the original enzyme. When the steady state is reached, it is assumed that the reaction rate is constant and follows the MM equation:This theory has been proven to be true for a wide variety of isolated enzymatic reactions in vitro (3) as well as more complicated reaction schemes involving multiple ligands and/or multisubunit enzymes (4). Recently, however, studies of enzyme dynamics together with single-molecule experiments have questioned the generality of existing kinetics steady-state laws for fluctuating enzyme systems (5-7). Steady-state kinetics were shown to depend on conformational transition rates of enzyme-substrate association, catalysis, and product release (6,8). In addition, the kinetics of simple or isolated reactions can differ from the kinetics of network reactions in a crowded environment as present in the cell (9, 10).Moreover, within the classical description, the contribution of the solvent in effecting catalysis is almost completely neglected, despite its potential role in mediating enzyme-substrate interactions (11). Importantly, the solvent plays an active role in various protein reactions (12-17), mediates molecular recognition (18-20), or acts as an adhesive to facilitat...