The reaction dynamics of excited-state intramolecular proton transfer (ESIPT) tautomerism in 10-hydroxybenzoquinoline (HBQ), 1-hydroxyanthraquinone (HAQ), methyl salicylate (MS) and 4-methyl-2,6-diformyl phenol (MFOH) has been investigated by quantum mechanics/molecular mechanics–molecular dynamics (QM/MM–MD) simulation with four different solvents H₂O, CH₃OH, CH₂Cl₂ and CHCl₃. The ESIPT of HBQ proceeds 20 ± 8 fs in water and methanol solvents; whereas, for other two solvents, CH₂Cl₂ and CHCl₃, proton transfer (PT) occurs within 62 ± 6 fs and 76 ± 6 fs, respectively. Similarly, the ESIPT of HAQ proceeds 21 ± 6 fs in H₂O, 33 ± 6 fs in CH₃OH, 51 ± 6 fs in CH₂Cl₂ and 61 ± 6 fs in CHCl₃. The ESIPT of MS proceeds in 23 ± 6 fs in H₂O, 38 ± 6 fs in CH₃OH, 60 ± 6 fs in CH₂Cl₂ and 73 ± 6 fs in CHCl₃. The ESIPT of MFOH proceeds 20 ± 6, 22 ± 6, 35 ± 6 and 58 ± 6 fs in H₂O, CH₃OH, CH₂Cl₂ and CHCl₃ solvent, respectively. Our QM/MM–MD simulations show that the rate of PT is inversely proportional to the polarity of the solvents. In addition, intramolecular hydrogen bonding capacity between the hydroxyl oxygen (proton donor) and the benzoquinolinic nitrogen (proton acceptor) plays an important role for the ESIPT reaction.