Highly-accurate numerical methods that can efficiently handle problems with interfaces and/or problems in domains with complex geometry are crucial for the resolution of different temporal and spatial scales in many problems from physics and biology. In this paper we continue the work started in [8], and we use modest one-dimensional parabolic problems as the initial step towards the development of high-order accurate methods based on the Difference Potentials approach. The designed methods are well-suited for variable coefficient parabolic models in heterogeneous media and/or models with non-matching interfaces and with non-matching grids. Numerical experiments are provided to illustrate high-order accuracy and efficiency of the developed schemes. While the method and analysis are simpler in the one-dimensional settings, they illustrate and test several important ideas and capabilities of the developed approach.