[step:Derive the moment bound for bounded centered Lipschitz functions]Let $\mathcal{B}(E)$ denote the Borel $\sigma$-algebra generated by the metric topology on $E$. By the meaning of the hypothesis $T_1(C)$, $\mu$ is a probability measure on $(E,\mathcal{B}(E))$ with finite first moment and satisfies
\begin{align*}
W_1(\nu,\mu) \leq \sqrt{2C\,\operatorname{Ent}(\nu\mid\mu)}
\end{align*}
for every probability measure $\nu$ on $(E,\mathcal{B}(E))$ such that $\nu\ll\mu$ and $\operatorname{Ent}(\nu\mid\mu)<\infty$. Here $\operatorname{Ent}(\nu\mid\mu)$ denotes the relative entropy
\begin{align*}
\operatorname{Ent}(\nu\mid\mu):=\int_E \log\left(\frac{d\nu}{d\mu}(x)\right)\,d\nu(x),
\end{align*}
$\Pi(\nu,\mu)$ denotes the set of probability measures on $(E\times E,\mathcal{B}(E)\otimes\mathcal{B}(E))$ whose first marginal is $\nu$ and whose second marginal is $\mu$, and $W_1(\nu,\mu)$ denotes the first Wasserstein distance
\begin{align*}
W_1(\nu,\mu):=\inf_{\pi\in\Pi(\nu,\mu)}\int_{E\times E} d(x,y)\,d\pi(x,y).
\end{align*}
Let $G:E \to \mathbb{R}$ be a bounded Borel measurable $1$-Lipschitz function satisfying
\begin{align*}
\int_E G(x)\,d\mu(x) = 0.
\end{align*}
For $\lambda \geq 0$, define the logarithmic moment generating function $\Lambda:[0,\infty) \to \mathbb{R}$ by
\begin{align*}
\Lambda(\lambda) := \log\left(\int_E e^{\lambda G(x)}\,d\mu(x)\right).
\end{align*}
Since $G$ is bounded, $\Lambda$ is continuously differentiable and
\begin{align*}
\Lambda'(\lambda) = \frac{\int_E G(x)e^{\lambda G(x)}\,d\mu(x)}{\int_E e^{\lambda G(x)}\,d\mu(x)}.
\end{align*}
Let $\nu$ be a probability measure on $(E,\mathcal{B}(E))$ such that $\nu \ll \mu$, $\operatorname{Ent}(\nu\mid\mu)<\infty$, and $\nu$ has finite first moment. Because $G$ is bounded, the integrals $\int_E G(x)\,d\nu(x)$ and $\int_E G(y)\,d\mu(y)$ are finite. Because $\nu$ and $\mu$ have finite first moments, $W_1(\nu,\mu)$ is finite and the coupling costs below are meaningful. Since $G$ is centered and $1$-Lipschitz, every coupling $\pi \in \Pi(\nu,\mu)$ satisfies
\begin{align*}
\int_E G(x)\,d\nu(x) - \int_E G(y)\,d\mu(y)
= \int_{E \times E} \bigl(G(x)-G(y)\bigr)\,d\pi(x,y)
\leq \int_{E \times E} d(x,y)\,d\pi(x,y).
\end{align*}
Taking the infimum over $\pi \in \Pi(\nu,\mu)$ gives
\begin{align*}
\int_E G(x)\,d\nu(x) \leq W_1(\nu,\mu).
\end{align*}
The $T_1(C)$ hypothesis applies to such $\nu$, so
\begin{align*}
\int_E G(x)\,d\nu(x) \leq \sqrt{2C\,\operatorname{Ent}(\nu\mid\mu)}.
\end{align*}
We now prove the entropy variational step needed for the logarithmic moment bound. For $\lambda \geq 0$, define the tilted probability measure $\nu_\lambda$ on $(E,\mathcal{B}(E))$ by
\begin{align*}
\frac{d\nu_\lambda}{d\mu}(x) := \frac{e^{\lambda G(x)}}{\int_E e^{\lambda G(y)}\,d\mu(y)}.
\end{align*}
Since $G$ is bounded, $\nu_\lambda \ll \mu$, the density is bounded above and below by positive constants, and $\operatorname{Ent}(\nu_\lambda\mid\mu)<\infty$. To verify the finite-first-moment condition, fix a point $x_0\in E$ and let $M_\lambda:=\sup_{x\in E} d\nu_\lambda/d\mu(x)<\infty$. The finite-first-moment assumption on $\mu$ gives
\begin{align*}
\int_E d(x,x_0)\,d\nu_\lambda(x) \leq M_\lambda\int_E d(x,x_0)\,d\mu(x) < \infty.
\end{align*}
Thus $\nu_\lambda$ is admissible for the $T_1(C)$ hypothesis. Its entropy is
\begin{align*}
\operatorname{Ent}(\nu_\lambda\mid\mu)
= \int_E \log\left(\frac{d\nu_\lambda}{d\mu}(x)\right)\,d\nu_\lambda(x)
= \lambda\int_E G(x)\,d\nu_\lambda(x)-\Lambda(\lambda).
\end{align*}
Hence
\begin{align*}
\Lambda(\lambda)=\lambda\int_E G(x)\,d\nu_\lambda(x)-\operatorname{Ent}(\nu_\lambda\mid\mu).
\end{align*}
Applying the preceding transport-entropy estimate to $\nu_\lambda$ gives
\begin{align*}
\Lambda(\lambda) \leq \lambda\sqrt{2C\,\operatorname{Ent}(\nu_\lambda\mid\mu)}-\operatorname{Ent}(\nu_\lambda\mid\mu).
\end{align*}
For the scalar $r_\lambda:=\operatorname{Ent}(\nu_\lambda\mid\mu)\geq0$, completing the square gives
\begin{align*}
\lambda\sqrt{2Cr_\lambda}-r_\lambda
= \frac{C\lambda^2}{2}-\left(\sqrt{r_\lambda}-\lambda\sqrt{\frac{C}{2}}\right)^2
\leq \frac{C\lambda^2}{2}.
\end{align*}
Therefore
\begin{align*}
\Lambda(\lambda) \leq \frac{C\lambda^2}{2}.
\end{align*}
Thus, for every $\lambda \geq 0$,
\begin{align*}
\int_E e^{\lambda G(x)}\,d\mu(x) \leq \exp\left(\frac{C\lambda^2}{2}\right).
\end{align*}[/step]