Effect of spin–orbit coupling on an ab initio Gutzwiller approach for Pu