I briefly review results obtained within the variational Hamiltonian approach to Yang-Mills theory in Coulomb gauge and confront them with recent lattice data. The variational approach is extended to non-Gaussian wave functionals including three- and four-gluon kernels in the exponential of the vacuum wave functional and used to calculate the three-gluon vertex. A new functional renormalization group flow equation for Hamiltonian Yang--Mills theory in Coulomb gauge is solved for the gluon and ghost propagator under the assumption of ghost dominance. The results are compared to those obtained in the variational approach.