Direct Numerical Simulation (DNS) of decaying isotropic 3D magnetohydrodynamic (MHD) turbulence based on the 10243-modes in a periodic box is used to study the statistical properties of turbulence. In this paper, the presence of intermittency in MHD turbulence is investigated through the analysis of the Probability Distribution Function (PDF) for Elsässer fields and total energy fluctuations. We observe that the PDFs of the Elsässer fields fluctuations display a strong non-Gaussian behavior at small scale, which can be ascribed to multifractality feature, while the PDFs of the total energy fluctuations have the same shape over all observed scales and are monofractal. The PDFs have stretched exponential tail and satisfy the function P(|δX|) ∼ exp(−A|δX|μ). Numerically, we extract the exponent μ and find that it is constant for monofractal behavior as the length scale varies. To check the notion of self-similarity in the respective fluctuation, we apply the compensated structure functions.