The topological and dynamical features of small scales are studied in the context of decaying magnetohydrodynamic turbulent flows using direct numerical simulations. Joint probability density functions (PDFs) of the invariants of gradient quantities related to the velocity and the magnetic fields demonstrate that structures and dynamics at the time of maximum dissipation depend on the large scale initial conditions at the examined Reynolds numbers. This is evident in particular from the fact that each flow has a different shape for the joint PDF of the invariants of the velocity gradient in contrast to the universal teardrop shape of hydrodynamic turbulence. The general picture that emerges from the analysis of the invariants is that regions of high vorticity are correlated with regions of high strain rate S also in contrast to hydrodynamic turbulent flows. Magnetic strain dominated regions are also well correlated with region of high current density j. Viscous dissipation (\documentclass[12pt]{minimal}\begin{document}${\propto } \bm S^2$\end{document}∝S2) as well as Ohmic dissipation (\documentclass[12pt]{minimal}\begin{document}${\propto } \bm j^2$\end{document}∝j2) resides in regions where strain and rotation are locally almost in balance. The structures related to the velocity gradient possess different characteristics than those associated with the magnetic field gradient with the latter being locally more quasi-two dimensional.