Regularized relativistic corrections for polyelectronic and polyatomic systems with explicitly correlated Gaussians