You can consider getting the NUTS regions via [{giscoR}](https://dieghernan.github.io/giscoR/) package - it serves formally valid NUTS regions.

    nuts <- giscoR::gisco_get_nuts(year = '2016')
    all(sf::st_is_valid(nuts))
    [1] TRUE