Starting with a substitution tiling, we demonstrate a method for constructing infinitely many new substitution tilings. Each of these new tilings is derived from a graph iterated function system and the tiles have fractal boundary. We show that each of the new tilings is mutually locally derivable to the original tiling. Thus, at the tiling space level, the new substitution rules are expressing geometric and combinatorial, rather than topological, features of the original. Our method is easy to apply to particular substitution tilings, permits experimentation, and can be used to construct border-forcing substitution rules. For a large class of examples we show that the combinatorial dual tiling has a realization as a substitution tiling. Since the boundaries of our new tilings are fractal we are led to compute their fractal dimension. As an application of our techniques we show how to compute the Čech cohomology of a (not necessarily border-forcing) tiling using a graph iterated function system of a fractal tiling.