A variational model for anisotropic and naturally twisted ribbons