A mixed meshless local Petrov–Galerkin (MLPG) collocation method is proposed for solving the two-dimensional boundary value problem of heterogeneous structures. The heterogeneous structures are defined by partitioning the total material domain into subdomains with different linear-elastic isotropic properties which define homogeneous materials. The discretization and approximation of unknown field variables is done for each homogeneous material independently, therein the interface of the homogeneous materials is discretized with overlapping nodes. For the approximation, the moving least square method with the imposed interpolation condition is utilized. The solution for the entire heterogeneous structure is obtained by enforcing displacement continuity and traction reciprocity conditions at the nodes representing the interface boundary. The accuracy and numerical efficiency of the proposed mixed MLPG collocation method is demonstrated by numerical examples. The obtained results are compared with a standard fully displacement (primal) meshless approach as well as with available analytical and numerical solutions. Excellent agreement of the solutions is exhibited and a more robust, superior and stable modeling of material discontinuity is achieved using the mixed method.