Here's a proof.
Using Lagrange's identity leads to
$$\sum_{cyc}\frac{x}{y+z}=\frac32+\frac12\sum_{cyc}\frac{(x-y)^2}{(x+z)(y+z)} \qquad (*)$$
We get this by letting $a = \sqrt{x+y}$, $b = \sqrt{y+z}$, $c = \sqrt{z+x}$,
$d = \frac{1}{\sqrt{x+y}}$, $e = \frac{1}{\sqrt{y+z}}$, $f = \frac{1}{\sqrt{z+x}}$, and by Lagrange's identity,
$$(a^2+b^2+c^2)(d^2+e^2+f^2) = (ad+be+cf)^2 + (ae-bd)^2+(af-cd)^2+(bf-ce)^2$$
which gives
$$2 (x+y+z)(\sum_{cyc}\frac{1}{y+z}) = 9 + \sum_{cyc} \left(\frac{\sqrt{x+y}}{\sqrt{y+z}} - \frac{\sqrt{y+z}}{\sqrt{x+y}}\right)^2$$
or
$$3 +\sum_{cyc}\frac{x}{y+z} = \frac92 + \frac12\sum_{cyc} \frac{\left({x+y} - (y+z)\right)^2}{({y+z})(x+y)} $$
which is the desired equation $(*)$.
Squaring both sides of the inequality then gives
$$
\tag{1} \frac16 \left[ \sum_{cyc}\frac{(x-y)^2}{(x+z)(y+z)} \right ]^2 + \sum_{cyc}\frac{(x-y)^2}{(x+z)(y+z)} \geq \frac{(y-z)^2}{xy + yz + xz}
$$
We will follow two paths for separate cases.
Path 1:
Omitting the square term it suffices to prove
$$\sum_{cyc}\frac{(x-y)^2}{(x+z)(y+z)} \geqslant \frac{(y-z)^2}{xy + yz + xz}$$
Clearing denominators, we obtain
$$
(x-y)^2 (x+y) + (y-z)^2 (y+z) + (z-x)^2 (z+x) \geq \frac{(y-z)^2}{xy + yz + xz} (x+y) (y+z) (z+x)
$$
Using $(x+y)(y+z)(x+z) = (x+y+z)(xy + yz+ xz) - x y z$
it suffices to show
$$
(x-y)^2 (x+y) + (y-z)^2 (y+z) + (z-x)^2 (z+x) \geq {(y-z)^2} (x+y+z)
$$
or
$$
(x-y)^2 (x+y) + (z-x)^2 (z+x) \geq (y-z)^2 x
$$
Since
$$
(y-z)^2 = (y-x + x- z)^2 = (y-x)^2 + (x- z)^2 + 2 (y-x)(x-z)
$$
this translates into
$$
(y-x)^2 y + (x-z)^2 z \geq 2 x(y-x)(x-z)
$$
For the two cases $y\geq x ; z\geq x $ and $y\leq x ; z\leq x $ the RHS $\leq 0$ so we are done. For the other two cases, by symmetry, it remains to show the case $y> x ; z < x $.
Rearranging terms, we can also write
$$
(y-x)^3 - (x-z)^3 +x ((y - x) + (z-x))^2 \geq 0
$$
This holds true at least for $(y-x)^3 \geq (x-z)^3$ or $y+z\geq 2 x$.
So the proof is complete other than for the case $y+z < 2x$ and [ $y> x ; z < x $ or $z> x ; y < x $ ].
Path 2.
For the remaining case $y+z < 2 x $ and [$y> x ; z < x$ or $z> x ; y < x $] we will follow a different path. Again, by symmetry, we must inspect only $y+z < 2x$ and $y> x > z$.
A remark up front: In the following, some high order polynomials of one variable have to be inspected. MATLAB is used, also for plotting behaviours of these polynomials. There is no case in spending effort for further analytical work on polynomials where their behaviour is obvious. Still, what follows contains some "ugly" parts.
In the squared version (1) of the inequality, we can use a further inequality which has been proved here:
$$\sum_{cyc}\frac{(x-y)^2}{(x+z)(y+z)} \geqslant \frac{27}{8} \frac{(y-z)^2}{(x+y+z)^2}$$
So it suffices to prove
$$
\frac16 \left[ \frac{27}{8} \frac{(y-z)^2}{(x+y+z)^2} \right ]^2 + \sum_{cyc}\frac{(x-y)^2}{(x+z)(y+z)} \geq \frac{(y-z)^2}{xy + yz + xz}
$$
Some numerical inspection shows immediately that the first term cannot be ommitted.
Clearing some denominators gives
$$
\tag{2} \frac16 \left[ \frac{27}{8} \right ]^2 \frac{(y-z)^4 (xy + yz + xz) (x+y)(y+z) (z+x)}{(x+y+z)^4} + \\
(xy + yz + xz) \sum_{cyc} {(x-y)^2}{(x+y)} - (y-z)^2 (x+y)(y+z) (z+x) \geqslant 0 \quad
$$
By homogeneity, we set $y=1+z$.
The condition $y+z < 2x$ then translates into $1+2z < 2x$, hence we further set $x = z + (1 +q)/2$ where $0\leq q \leq 1$ since also $x = z + (1 +q)/2 < y = 1 +z$.
Inserting $y=1+z$ and $x = z + (1 +q)/2$ into (2) is straightforward, the result is lengthy (not displayed here).
Let us start with focussing on the first term in (2), calling that fraction $F$:
$$
F= \frac{(y-z)^4 (xy + yz + xz) (x+y)(y+z) (z+x)}{(x+y+z)^4}
$$
With the setting $y=1+z$ this can be simplified to
$$
F = \frac{(xy + yz + xz) (x+y)(y+z) (z+x)}{(x+y+z)^4}
$$
Since by the settings $y=1+z$ and $x = z + (1 +q)/2$ both x and y are linear in z, the numerator of $F$ is of fifth order in $z$, whereas the denominator is of fourth order in $z$. In leading order, the whole term is therefore of first order in $z$ and will therefore rise with $z$ for large enough $z$. This motivates to show that indeed
$$F(q,z) \geq F(q,z=0) = 2 \frac{(q+1)^2}{(q+3)^2}$$
for all $z$ and $q$.
Showing that directly requires a condition
$$
G = (xy + yz + xz) (x+y)(y+z) (z+x) (q+3)^2 - 2 (q+1)^2 (x+y+z)^4 \geq 0
$$
Inserting $y=1+z$ and $x = z + (1 +q)/2$ into $G$, and expanding the brackets, gives a very lengthy expression which however contains only positive terms, so the condition is proved immediately:
$$
G = (z(2q^6z + 2q^6 + 22q^5z^2 + 51q^5z + 22q^5 + 80q^4z^3 + 358q^4z^2 + 357q^4z + 100q^4 + 96q^3z^4 + 960q^3z^3 + 1828q^3z^2 + 1206q^3z + 260q^3 + 864q^2z^4 + 3672q^2z^3 + 4788q^2z^2 + 2484q^2z + 450q^2 + 2592qz^4 + 7344qz^3 + 7398qz^2 + 3159qz + 486q + 2592z^4 + 5832z^3 + 4806z^2 + 1701z + 216))/4 \geq 0
$$
Hence it suffices, instead of (2), to prove the following:
$$
\tag{3} \frac16 \left[ \frac{27}{8} \right ]^2 2 \frac{(q+1)^2}{(q+3)^2} +
(xy + yz + xz) \sum_{cyc} {(x-y)^2}{(x+y)} \\ - (x+y)(y+z) (z+x) \geqslant 0 \quad
$$
After insertion of $y=1+z$ and $x = z + (1 +q)/2$, the factors $(x-y)^2$ in the cyclic sum will not be functions of $z$. Hence,
the LHS is a third order expression in $z$ with leading (in $z^3$ ) term $( 1 + 3 q^2) z^3$, so for large enough $z$ it is rising with $z$. A remarkable feature of this expression is that for the considered range $0\leq q \leq 1$ it is actually monotonously rising for all $z$. To see this, consider whether there are points with zero slope. The first derivative of the expression with respect to $z$ is
$$
q^4/4 + (7q^3z)/2 + (7q^3)/4 + 9q^2z^2 + 9q^2z + (5q^2)/4 - (7qz)/2 - (7q)/4 + 3z^2 + 3z + 1/2
$$
Equating this to zero gives
$$
z_{1} =
-(2((13q^6)/4 + (17q^4)/2 + (133q^2)/4 + 3)^{(1/2)} - 7q + 18q^2 + 7q^3 + 6)/(36q^2 + 12)\\
z_{2} =
-(-2((13q^6)/4 + (17q^4)/2 + (133q^2)/4 + 3)^{(1/2)} - 7q + 18q^2 + 7q^3 + 6)/(36q^2 + 12)
$$
One now shows that in the range $0\leq q \leq 1$, there are only negative solutions $z_{1,2}$. So we will have no zero slopes for $z\geq 0$. Since the polynomials are (roots of) sixth order in q, we investigate the following figures:
So monotonicity (rising with $z$) is established for all $q$.
Hence, to show the inequality it suffices to inspect (3) at the smallest $z=0$. This gives
$$
((q + 1)(8q^6 + 88q^5 + 336q^4 + 432q^3 - 216q^2 - 405q + 243))/(64(q + 3)^3) \geqslant 0 \quad
$$
or
$$
\tag{4} 8q^6 + 88q^5 + 336q^4 + 432q^3 - 216q^2 - 405q + 243 \geqslant 0 \quad
$$
An even stronger requirement is
$$
h(q) =
432q^3 - 216q^2 - 405q + 243 \geq 0
$$
In the considered range $0\leq q \leq 1$, $h(q)$ has a minimum which is obtained by taking the first derivative,
$$
1296 q^2 - 432 q - 405
$$
and equating to zero, which gives $q = 3/4$, and the above $h(q)$ then gives
$$
h(q = 3/4) = 0
$$
This establishes the inequality. $ \qquad \Box$
Best Answer
Thus, $$(a^4+b^4+c^4)(3(a^2+b^2+c^2)-(a+b+c)^2)-25\prod_{cyc}(a-b)^2=$$ $$=((c+u)^4+(c+v)^4+c^4)\sum_{cyc}(a-b)^2-25\prod_{cyc}(a-b)^2\geq$$ $$=(u^4+v^4)(u^2+v^2+(u-v)^2)-25(u-v)^2u^2v^2=$$ $$=u^3v^3((4t^2-2)(4t-2)-25(2t-2))=u^3v^3(16t^3-8t^2-58t+54)=$$ $$=u^3v^3(t(4t-5)^2+32t^2-83t+54)\geq0$$ because $$83^2-4\cdot32\cdot54=-23<0.$$
The number $48$ is a best constant.