An advanced constitutive model incorporating two specific aspects of Ni-base superalloy deformation behaviour is proposed. Several deformation mechanisms are active in these two-phase materials. In the matrix phase, cube slip plays an important role in the orientation dependence of the material. Moreover, inelastic deformation of the precipitate phase leads to non-Schmid effects in the material response. Macroscopic cube slip is modelled here by incorporating a zig-zag cross slip mechanism into the constitutive relations for the matrix phase. A cross slip factor is proposed that quantifies the amount of cross slip and consequently represents the orientation dependence of the cube slip. Further, a detailed precipitate phase constitutive model is proposed, which enables the simulation of non-Schmid effects, like the tension–compression asymmetry. The cross slip mechanism and the associated splitting of partial dislocations in the γ′-phase, which are responsible for the anomalous yield behaviour, are incorporated in the model. The proposed formulations are implemented in a recently developed crystal plasticity framework for single crystal Ni-base superalloys and a consistent set of model parameters for the commercial alloy CMSX-4 is determined. The model is shown to reasonably predict the material tensile response and creep behaviour for a range of temperatures and stress or strain rate levels. The incorporation of the cross slip mechanisms in the matrix and precipitate results in an adequate simulation of the material orientation dependence and the experimentally determined tension–compression asymmetry.