Solution hardening in body-centered cubic quaternary alloys interpreted using Suzuki's kink-solute interaction model