Hyperbolic heat conduction equation with temperature dependent thermal properties is considered. The thermal conductivity, specific heat and density are assumed to be functions of temperature. The equation is cast into a non-dimensional form suitable for perturbation analysis. By employing a newly developed approximate symmetry theory, the approximate symmetries of the equation are calculated for the case of small variations in thermal properties. Various similarity solutions corresponding to the symmetries of first order equations are presented. For second order equations, the method of constructing approximate symmetries and similarity solutions are discussed. A linear functional variation is assumed for the thermal properties and a similarity solution is constructed using one of the first order solutions as an example.