We investigate the isospin symmetry breaking effects within a recently derived Nambu-Jona-Lasinio related model by fitting the measured pseudoscalar meson masses and weak decay couplings f π , f K . Our model contains the next to leading order terms in the 1=N c expansion of the effective multiquark Lagrangian, including the ones that break the chiral symmetry explicitly. We show the important phenomenological role of these interactions: (1) They lead to an accurate fit of the low-lying pseudoscalar nonet characteristics. (2) They account for a very good agreement of the current quark masses with the present PDG values. (3) They reduce by 40% the ratio ϵ=ϵ 0 of the π 0 − η and π 0 − η 0 mixing angles, as compared to the case that contemplates explicit breaking only in the leading order, bringing it in consonance with the quoted values in the literature. The conventional NJL-type models fail in the joint description of these parameters.