The heterogeneous multi-scale method (HMM) provides an efficient framework for constructing stable algorithms for tackling a large class of problems with multiple scales. In this work, we apply the HMM framework to interfaces, which may represent the boundary of a tumor or the growth of a crystal, in heterogeneous media. The goal is to capture the macroscale interface dynamics and is achieved through a macroscale solver, in this case the level set method, that obtains essential, unknown information through numerical experiments in the microscale rather than constitutive relations which may be difficult to derive. We present results revealing the details and demonstrating the effectiveness and generality of the approach.