CP violation in neutrino interactions is described by three phases contained in Pontecorvo-Maki-Nakagawa-Sakata mixing matrix (U PMNS ). We argue that the phenomenologically consistent result of the Dirac CP violation can be obtained if U PMNS is constructed along bipair neutrino mixing scheme, namely, requiring that |U 12 | = |U 32 | and |U 22 | = |U 23 | (case 1) and |U 12 | = |U 22 | and |U 32 | = |U 33 | (case 2), where U ij stands for the i × j matrix element of U PMNS . As a result, the solar, atmospheric and reactor neutrino mixing angles θ 12 , θ 23 and θ 13 , respectively, are correlated to satisfy cos 2θ 12 = sin 2 θ 23 − tan 2 θ 13 (case 1) or cos 2θ 12 = cos 2 θ 23 − tan 2 θ 13 (case 2). Furthermore, if Dirac CP violation is observed to be maximal, θ 23 is determined by θ 13 to be: sin 2 θ 23 ≈ √ 2 − 1 cos 2 θ 13 + √ 2 sin 2 θ 13 (case 1) or cos 2 θ 23 ≈ √ 2 − 1 cos 2 θ 13 + √ 2 sin 2 θ 13 (case 2). For the case of non-maximal Dirac CP violation, we perform numerical computation to show relations between the CP-violating Dirac phase and the mixing angles.