Charged spherically symmetric NUT and Taub-NUT black hole solutions in quadratic $f(R)$ gravity