Matrix calculus-based approach to orthogonal polynomial sequences