Let $E\subseteq\mathbb R^d$ be a convex set, $\beta\geq 0$ be a given real number and $f:E\to\mathbb R$ be a convex and differentiable function satisfying: $$f(y)\leq f(x)+\nabla f(x)^\top (y-x) +\frac{\beta}{2}\|x-y\|_2^2, \quad \forall x,y\in E.$$ Show that $\nabla f$ is $\beta$-Lipschitz, i.e., $$\|\nabla f(x)-\nabla f(y)\|_2\leq\beta\|x-y\|_2, \quad \forall x,y\in E.$$
Edit: I am interested in the case when $E$ can be any convex subset of $\mathbb R^d$. When $E=\mathbb R^d$, the claim can be proven as proposed by the answer below. However, in many convex optimisation textbooks, the claim is stated for general domains $E$ without proof.