Functional renormalization group for fermion lattice models in three dimensions: Application to the Hubbard model on the cubic lattice