A novel method based on diffusion NMR for the epitope mapping of ligand binding is presented. The intermolecular NOE builds up during a long diffusion period and creates a deviation from the linearity. The ligand proton nearest the protein generates the strongest NOE from protein during the diffusion period and has the largest deviation. Therefore, this diffusion artifact can be used to characterize the ligand binding epitope. The concept was investigated using dihydrofolate reductase (DHFR) and its ligand trimethoprim (TMP), and the epitope map of TMP on DHFR generated with this method is in excellent agreement with the structural and dynamic studies by crystallography and NMR, as well as the medicinal chemistry results.