We compute the one-loop corrected effective Lagrangian for the neutralino-neutralino-neutral Higgs interactions ${\ensuremath{\chi}}_{\ensuremath{\ell}}^{0}{\ensuremath{\chi}}_{k}^{0}{H}_{m}^{0}$. The analysis completes the previous analyses where similar corrections were computed for the $\overline{f}f{H}_{m}^{0}$ couplings, where $f$ stands for standard model quarks and leptons and for the chargino-chargino-neutral Higgs couplings ${\ensuremath{\chi}}_{l}^{+}{\ensuremath{\chi}}_{k}^{\ensuremath{-}}{H}_{m}^{0}$ within the minimal supersymmetric standard model. The effective one-loop Lagrangian is then applied to the computation of the neutral Higgs decays. The sizes of the supersymmetric loop corrections of the neutral Higgs decay widths into ${\ensuremath{\chi}}_{\ensuremath{\ell}}^{0}{\ensuremath{\chi}}_{k}^{0}$ ($\ensuremath{\ell}=1$, 2, 3, 4; $k=1$, 2, 3, 4) are investigated and the supersymmetric loop correction is found to be in the range of $\ensuremath{\sim}10%$ in significant regions of the parameter space. By including the loop corrections of the other decay channels $\overline{b}b$, $\overline{t}t$, $\overline{\ensuremath{\tau}}\ensuremath{\tau}$, $\overline{c}c$, and ${\ensuremath{\chi}}_{i}^{\ensuremath{-}}{\ensuremath{\chi}}_{j}^{+}$ ($i=1$, 2; $j=1$, 2), the corrections to branching ratios for ${H}_{m}^{0}\ensuremath{\rightarrow}{\ensuremath{\chi}}_{\ensuremath{\ell}}^{0}{\ensuremath{\chi}}_{k}^{0}$ can reach as high as 50%. The effects of $CP$ phases on the branching ratio are also investigated. A discussion of the implications of the analysis for colliders is given.