We know that
$$
1+\tau_iL_i=\frac{1}{P(t_i,t_{i+1})}
$$
and, under the risk-neutral measure,
$$
P(t,T)=P(0,T)\exp\left(\int_0^tr(s)\,ds+\int_0^t\sigma(s,T)\,dW_s-\frac{1}{2}\int_0^t\sigma^2(s,T)\,ds\right)\,.\quad\quad\quad (1)
$$
Using $P(t_i,t_i)=1$ this implies
\begin{align}
P(t_i,t_{i+1})&=\frac{P(0,t_{i+1})}{P(0,t_i)}\exp\Bigg(\int_0^{t_i}\sigma(s,t_{i+1})-\sigma(s,t_i)\,dW_s\\
&\quad\quad-\frac{1}{2}\int_0^{t_i}\sigma^2(s,t_{i+1})-\sigma^2(s,t_i)\,ds \Bigg)\,.
\end{align}
Therefore,
\begin{align}
\prod_{i=0}^{n-1}1+\tau_iL_i&=\frac{P(0,t_0)}{P(0,t_n)}\exp\Bigg(\sum_{i=0}^{n-1}\int_0^{t_i}\sigma(s,t_i)-\sigma(s,t_{i+1})\,dW_s\\
&\quad\quad-\sum_{i=0}^{n-1}\frac{1}{2}\int_0^{t_i}\sigma^2(s,t_i)-\sigma^2(s,t_{i+1})\,ds \Bigg)\,.
\end{align}
Further, (1) implies with $P(T,T)=1$ that
$$
\exp\left(-\int_0^{T}r(s)\,ds\right)=P(0,T)\exp\left(\int_0^T\sigma(s,T)\,dW_s-\frac{1}{2}\int_0^T\sigma^2(s,T)\,ds\right)\,.
$$
So we have to calculate the expectation of
\begin{align}
\exp\left(-\int_0^{T}r(s)\,ds\right)\left(\prod_{i=0}^{n-1}1+\tau_iL_i\right)
\end{align}
which is
$$
\frac{P(0,t_0)P(0,T)}{P(0,t_n)}\exp\left(-\sum_{i=0}^{n-1}\frac{1}{2}\int_0^{t_i}\sigma^2(s,t_i)-\sigma^2(s,t_{i+1})\,ds-\frac{1}{2}\int_0^T\sigma^2(s,T)\,ds\right)
$$
times the expectation of the lognormal variable
$$
e^Y:=\exp\left(\int_0^T\sigma(s,T)\,dW_s+\sum_{i=0}^{n-1}\int_0^{t_i}\sigma(s,t_i)-\sigma(s,t_{i+1})\,dW_s\right)\,.
$$
To calculate that expectation you need to know only the variance of $Y$ which is
a sum of integrals of the form
\begin{align}
&\int_0^{t_i}(\sigma(s,t_i)-\sigma(t_{i+1}))(\sigma(s,t_j)-\sigma(t_{j+1}))\,ds\quad\quad i\le j\\
&\int_0^{t_i}(\sigma(s,t_i)-\sigma(t_{i+1}))\sigma(s,T)\,ds\,.
\end{align}
This is as explicit as it gets. Even
using the time dependent HW-form of $\sigma(t,T)\,,$ that is,
$$
\sigma(t,T)=\int_t^T\sigma(s)\exp\left(-\int_s^T\lambda(u)\,du\right)\,ds
$$
does not help much, unless you assume constant mean reversion $\lambda$
and constant volatility $\sigma$ of the short rate.