The thermal lubrication of an entangled polymeric liquid in wall-driven shear flows between parallel plates is investigated by using a multiscale hybrid method coupling molecular dynamics and the hydrodynamics (i.e., the synchronized molecular dynamics method). The temperature of the polymeric liquid rapidly increases due to viscous heating once the drive force exceeds a certain threshold value. The rheological properties of the polymeric liquid drastically change at around the critical drive force. In the weak viscous-heating regime, the conformation of polymer chains is dominated by the local shear flow so that the anisotropy of the bond orientation tensor grows as the drive force increases. However, in the large viscous-heating regime, the conformation dynamics is dominated by the thermal agitation of polymer chains so that the bond orientation tensor recovers more uniform and random structures as the drive force increases, even though the local shear flows are further enhanced. Remarkably, these counter-intuitive transitional behaviors give an interesting re-entrant transition in the stress--optical relation, where a linear formalism in the stress--optical relation approximately holds even though each of the macroscopic quantities behaves nonlinearly. The robustness of the linear stress--optical relation is also confirmed in the spatiotemporal evolution at the hydrodynamic level.