We develop a new approach, based on quantization methods, to study higher symmetries of invariant differential operators. We focus here on conformally invariant powers of the Laplacian over a conformally flat manifold, and recover results of Eastwood, Leistner, Gover and Silhan. In particular, conformally equivariant quantization establishes a crystal clear correspondence between hamiltonian symmetries of the null geodesic flow and the algebra of higher symmetries of the conformal Laplacian. Resorting to symplectic reduction, this leads to a quantization of the minimal nilpotent coadjoint orbit of the conformal group and allows to identify the latter algebra of symmetries in terms of the Joseph ideal. By the way, we obtain a tangential star-product for a family of coadjoint orbits of the conformal group.