The present work deals with dynamical system analysis of Holographic Dark Energy cosmological model with different infra-red (IR)-cutoff. By suitable transformation of variables the Einstein field equations are converted to an autonomous system. The critical points are determined and the stability of the equilibrium points are examined by Center Manifold Theory and Liapunov function method. Possible bifurcation scenarios have also been explained. 98.80.Cq, 95.35.+d, 95.36.+x, 0.Ud, 6402.1.60.Fr,