Variational Monte Carlo for spin-orbit interacting systems