Predator-prey models with Michaelis-Menten-Holling type ratio- dependent functional response exhibit very rich and complex dynamical behavior, such as the existence of degenerate equilibria, appearance of limit cycles and heteroclinic loops, and the coexistence of two attractive equilibria. In this paper, we study heteroclinic bifurcations of such a predator-prey model. We first calculate the higher order Melnikov functions by transforming the model into a Hamiltonian system and then provide an algorithm for computing higher order approximations of the heteroclinic bifurcation curves.