One answer - subgaussian variables generalize this property.
Let $\mu = \sum_i p_i z_i$, then the distribution is considered $\sigma^2$-subgaussian if for all $\lambda \in \mathbb{R}$,
$$ \log\left(\sum_i p_i e^{\lambda(z_i - \mu)}\right) \leq \frac{\lambda^2 \sigma^2}{2} $$
i.e.
$$ \log\left(e^{-\lambda \mu} \sum_i p_i e^{\lambda z_i}\right) \leq \frac{\lambda^2 \sigma^2}{2} $$
i.e.
$$ \log\left(\sum_i p_i e^{\lambda z_i}\right) - \lambda \mu \leq \frac{\lambda^2 \sigma^2}{2} . $$
Your expression is the case $\lambda=1$.
(In a sense this is not an answer, just a definition for when and how your inequality can be satisfied, but hopefully useful.)