With computer, we can use the branch and bound strategy to prove the inequality.
Proof: The inequality is written as
$$F(a, b, c, d) = \mathrm{e}^{ba\ln a} + \mathrm{e}^{cb\ln b} + \mathrm{e}^{dc\ln c} + \mathrm{e}^{a d\ln d} > \pi.$$
We first introduce some auxiliary results (Facts 1 through 2). The proofs are given later.
Fact 1: Given four real numbers $0 \le A < B, \ 0 \le C < D$ with $B-A \le \frac{1}{8}, \ D-C\le \frac{1}{8}$. Let $x\in [A, B]$ and $y\in [C, D]$. For convenience, we set $0\ln 0 = 0$. We have
i) If $\frac{1}{\mathrm{e}} \in [A, B]$, then $yx\ln x \ge -D \mathrm{e}^{-1}$;
ii) If $\frac{1}{\mathrm{e}} \notin [A, B]$ and $ A < 1$, then $yx\ln x \ge D \min(A\ln A, B\ln B)$;
iii) If $A \ge 1$, then $yx\ln x \ge C A \ln A$.
Fact 2: If $a\in [0, 1], \ b \in [0, 2], \ c \in [0, 1], \ d\ge 0$ and $a+b+c+d=4$, then $F(a,b,c,d) > \pi$.
Now let us proceed. WLOG, assume that $d= \max(a,b,c,d)$. Clearly $b \le 2$.
If $a > 1$, we have
$$\mathrm{e}^{ba\ln a} + \mathrm{e}^{cb\ln b} + \mathrm{e}^{dc\ln c} + \mathrm{e}^{a d\ln d} \ge 2\mathrm{e}^{\frac{1}{2}(cb\ln b + dc\ln c)} + 2
\ge 2\mathrm{e}^{\frac{1}{2}(-\frac{3}{\mathrm{e}})} + 2 > \pi$$
since $x\ln x \ge -\frac{1}{\mathrm{e}}$ for all $x > 0$, and $x\mapsto \mathrm{e}^x$ is convex.
If $c > 1$, we have
$$\mathrm{e}^{ba\ln a} + \mathrm{e}^{cb\ln b} + \mathrm{e}^{dc\ln c} + \mathrm{e}^{a d\ln d} \ge 2\mathrm{e}^{\frac{1}{2}(ba\ln a + cb\ln b)} + 2
\ge 2\mathrm{e}^{\frac{1}{2}(-\frac{b+c}{\mathrm{e}})} + 2
\ge 2\mathrm{e}^{\frac{1}{2}(-\frac{8/3}{\mathrm{e}})} + 2 > \pi.$$
If $a \le 1$ and $c\le 1$, from Fact 2, the desired result follows.
We are done.
$\phantom{2}$
Proof of Fact 1: Let $f(x) = x\ln x$. Then $f(x)$ is strictly decreasing on $[0, \frac{1}{\mathrm{e}})$ and strictly increasing on $(\frac{1}{\mathrm{e}}, \infty)$.
Also $f(x)$ achieves its global minimum at $x = \frac{1}{\mathrm{e}}$. It is easy to prove the desired results.
Proof of Fact 2: We use the branch and bound strategy.
Let
$$\Omega = \{(a,b,c,d): \ a\in [0, 1], \ b \in [0, 2], \ c \in [0, 1], \ d\ge 0, \ a+b+c+d=4\}.$$
For $i=0, 1, \cdots, 127; \ j=0, 1, \cdots, 255; \ k = 0, 1, \cdots, 127$, let
\begin{align}
S(i,j,k) &= \Big\{(a,b,c,d): \
a \in \Big[\frac{i}{128}, \ \frac{i+1}{128}\Big],
\quad b \in \Big[\frac{j}{128}, \ \frac{j+1}{128}\Big], \\
&\qquad\quad c \in \Big[\frac{k}{128}, \ \frac{k+1}{128}\Big], \quad
d \in \Big[4 - \frac{3 + i + j + k}{128}, \ 4 - \frac{i + j + k}{128}\Big]\Big\}.
\end{align}
We have $\Omega \subset \bigcup_{i,j,k} S(i,j,k)$.
From Fact 1, we evaluate a lower bound $L(i,j,k)$ of $F(a,b,c,d)$ on $S(i,j,k)$, i.e.,
$$F(a,b,c,d) \ge L(i,j,k), \quad \forall (a,b,c,d) \in S(i,j,k).$$
Thus, we have
$$ F(a,b,c,d) \ge \min_{i,j,k} L(i,j,k), \quad \forall (a,b,c,d)\in \Omega.$$
With computer, we obtain $\min_{i,j,k} L(i,j,k) > \pi$. The desired result follows.
This completes the proof of Fact 2.
3.1605859174508652189…
, usingNMinimize[{(a^2)^(a^2 b^2) + (b^2)^(b^2 c^2) + (c^2)^(c^2 d^2) + \ (d^2)^(d^2 a^2), a != 0 && b != 0 && c != 0 && d != 0 && a^2 + b^2 + c^2 + d^2 == 4}, {a, b, c, d}, WorkingPrecision -> 100]
$\endgroup$