The explored solutions described some different solutions as, Lump soliton, a solitary wave and exponential solutions. These solutions are investigated through some new Lie infinitesimals for the (3 + 1) dimensional variable coefficients Kudryashov-Sinelshchikov (VCKS). We used the fourth prolongation to investigate fifteen cases of Lie vectors. In each case, there is an infinite number of possibilities of vectors due to the unknown arbitrary functions and the variable coefficients for the considered model. We selected one case and examined the commutative product between multi unknown Lie infinitesimals for the (3 + 1) dimensional (VCKS) equation and this complicated process resulted from some new Lie vectors. The commutative product generates a system of nonlinear ODEs which had been solved manually. Through three stages of Lie symmetry reduction using the equivalence transformation, (VCKS) equation is reduced to solvable nonlinear ODEs using various combinations of optimal Lie vectors. By solving these ODEs, we investigate new analytical solutions for these ODEs. Back substituting to the original variables generates new solutions for (VCKS). Some selected solutions are illustrated through three-dimensional plots.