Detailed information about the influences of the cooperative motion of water and methanol molecules on practical solid-liquid heterogeneous photocatalysis reactions is critical for our understanding of photocatalytic reactions. The present work addresses this issue by applying operando nuclear magnetic resonance (NMR) spectroscopy, in conjunction with density functional theory (DFT) calculations and ab initio molecular dynamics (AIMD) simulations, to investigate the dynamic behaviors of heterogeneous photocatalytic systems with different molar ratios of water to methanol on rutile-TiO2 photocatalyst. The results demonstrate that methanol and water molecules are involved in the cooperative motions, and the cooperation often takes the form of methanol-water clusters that govern the number of methanol molecules reaching to the active sites of the photocatalyst per unit time, as confirmed by the diffusion coefficients of the methanol molecule calculated in the binary methanol-water solutions. Nuclear Overhauser effect spectroscopy experiments reveal that the clusters are formed by the hydrogen bonding between the -OH groups of CH3OH and H2O. The formation of such methanol-water clusters is likely from an energetic standpoint in low-concentration methanol, which eventually determines the yields of methanol reforming products.