[step:Derive the exponential moment estimate for bounded Lipschitz functions]Let
\begin{align*}
g:X &\to \mathbb{R}
\end{align*}
be a bounded $1$-Lipschitz Borel function, and define its $\mu$-mean by
\begin{align*}
m_g := \int_X g\,d\mu(x).
\end{align*}
Fix $\lambda \geq 0$. Define the normalizing constant
\begin{align*}
Z_\lambda := \int_X e^{\lambda g}\,d\mu(x),
\end{align*}
and define the exponentially tilted probability measure $\nu_\lambda$ on the Borel $\sigma$-algebra of $X$ by
\begin{align*}
\nu_\lambda(A) := \frac{1}{Z_\lambda}\int_A e^{\lambda g}\,d\mu(x)
\end{align*}
for every Borel set $A \subset X$. Since $g$ is bounded, $0 < Z_\lambda < \infty$, and $\nu_\lambda \ll \mu$. Moreover, if $M_\lambda := \exp(\lambda\|g\|_\infty)/Z_\lambda$, then $d\nu_\lambda/d\mu \leq M_\lambda$, so
\begin{align*}
\int_X d(x,x_0)\,d\nu_\lambda(x) \leq M_\lambda \int_X d(x,x_0)\,d\mu(x) < \infty.
\end{align*}
Thus $\nu_\lambda \in \mathcal{P}_1(X)$, and the hypothesis $T_1(C)$ applies to $\nu_\lambda$. Its relative entropy is
\begin{align*}
H(\nu_\lambda \mid \mu)
= \int_X \log\left(\frac{d\nu_\lambda}{d\mu}\right)\,d\nu_\lambda(x)
= \lambda \int_X g\,d\nu_\lambda(x) - \log Z_\lambda.
\end{align*}
We first record the coupling estimate used for $W_1$. For every coupling $\pi$ of $\nu_\lambda$ and $\mu$, the function $g$ being $1$-Lipschitz gives
\begin{align*}
g(x)-g(y) \leq d(x,y)
\end{align*}
for every $(x,y) \in X \times X$. Integrating with respect to $\pi$ yields
\begin{align*}
\int_X g\,d\nu_\lambda(x)-\int_X g\,d\mu(y)
= \int_{X \times X} (g(x)-g(y))\,d\pi(x,y)
\leq \int_{X \times X} d(x,y)\,d\pi(x,y).
\end{align*}
Applying the same argument to the $1$-Lipschitz function $-g:X\to\mathbb{R}$ gives
\begin{align*}
\int_X g\,d\mu(y)-\int_X g\,d\nu_\lambda(x) \leq \int_{X \times X} d(x,y)\,d\pi(x,y).
\end{align*}
Taking the infimum over all couplings $\pi$ in both inequalities gives
\begin{align*}
\left|\int_X g\,d\nu_\lambda(x)-m_g\right| \leq W_1(\nu_\lambda,\mu).
\end{align*}
Set
\begin{align*}
a_\lambda := \int_X g\,d\nu_\lambda(x)-m_g.
\end{align*}
By the previous inequality and $T_1(C)$,
\begin{align*}
|a_\lambda| \leq \sqrt{2C H(\nu_\lambda \mid \mu)}.
\end{align*}
Since $H(\nu_\lambda \mid \mu)=\lambda a_\lambda+\lambda m_g-\log Z_\lambda$, this gives
\begin{align*}
a_\lambda^2 \leq 2C\left(\lambda a_\lambda+\lambda m_g-\log Z_\lambda\right).
\end{align*}
Therefore
\begin{align*}
\log Z_\lambda-\lambda m_g \leq \lambda a_\lambda-\frac{a_\lambda^2}{2C}.
\end{align*}
For every real number $a$, completing the square gives
\begin{align*}
\lambda a-\frac{a^2}{2C} \leq \frac{C\lambda^2}{2}.
\end{align*}
Applying this with $a=a_\lambda$ gives
\begin{align*}
\log \int_X e^{\lambda(g-m_g)}\,d\mu(x)
= \log Z_\lambda-\lambda m_g
\leq \frac{C\lambda^2}{2}.
\end{align*}
Exponentiating,
\begin{align*}
\int_X e^{\lambda(g-m_g)}\,d\mu(x) \leq \exp\left(\frac{C\lambda^2}{2}\right).
\end{align*}[/step]