Maximum likelihood estimation for three-dimensional cubic contingency tables under modified quasi-symmetry