Matrix Orthogonal Polynomials: A Riemann--Hilbert approach