Understanding the molecular recognition process of nucleobases is one of the greatest challenges for both computational chemistry and biophysics fields. In fact, our results point out that it is a hard task to take into account the hydrophobic interactions, such as pi-pi and T-stacking interactions, by theoretical calculations using conventional force fields due to quantum effects of hyperconjugation and electronic correlation. In this line, our findings put in evidence that simple modifications in the Lennard-Jones potential can improve theoretical predictions in scenarios where hydrophobic interactions can drive the molecular recognition.