Interpolation by hypercyclic functions for differential operators

We prove that, given a sequence of points in a complex domain Ω without accumulation points, there are functions having prescribed values at the points of the sequence and, simultaneously, having dense orbit in the space of holomorphic functions on Ω. The orbit is taken with respect to any fixed nonscalar differential operator generated by an entire function of subexponential type, thereby extending a recent result about MacLane-hypercyclicity due to Costakis, Vlachou and Niess.