Non-asymptotic convergence bounds for Sinkhorn iterates and their gradients: a coupling approach