We report on our on-going project to calculate the nucleon decay matrix elements with domain-wall fermions. Operator mixing is discussed employing a non-perturbative renormalization. Bare matrix elements of all the possible decay modes induced by the dimension-six operators are calculated with the direct method, which are compared with the indirect calculation using chiral perturbation theory.