I improve the result see (Elegant) proof of : $x \log_2\frac{1}{x}+(1-x) \log_2\frac{1}{1-x} \geq 1- (1-\frac{x}{1-x})^2$ :
What is the best constant $C$ such that $x\in(0,0.5]$:
$\frac{1}{\ln\left(2\right)}\left(x\ln\left(\frac{1}{x}\right)+\left(1-x\right)\ln\left(\frac{1}{1-x}\right)\right)+\frac{\ln\left(1-2\left(1-x\right)x\right)}{\ln\left(2\right)}\leq -\left(1-\left(1-\frac{x}{1-x}\right)^{2}+\frac{\ln\left(1-f\left(\frac{1}{2}\right)\left(\frac{1}{2^{C-1}}-x^{C}\right)x^{C}\right)}{\ln\left(2\right)}\right)$
Where :
$$f\left(x\right)=\frac{1}{2}\left(\frac{1}{2^{C-1}}-x^{C}\right)^{-1}x^{-C}$$
It seems there is a problem with the derivatives around zero .So to find $C$ it might be useful to have an asymptotic approximation around this limit of domain .
You can find an attempt of mine in the linked question and other answer which could be useful
Currently $C\simeq 0.76131$.
How to find the $C$ such that we have the inequality ?
Sides notes :
If $C$ as optimal constant does'nt admit a closed form we can attempt to show it for $C=e^{-3/11}$
We can also on the other hand use :
$$x\ln x+(1-x)\ln 1-x \geq\ln \left(1-\sqrt{x(1-x)}\right)$$
Wich can be solved without derivative here New bound for Am-Gm of 2 variables
After some hours it seems we have the equality :
$$\lim_{x\to 0}\frac{\ln\left(C\right)}{f\left(\frac{1}{2}\right)}-\frac{\ln\left(kC\right)}{\ln\left(\frac{k}{d}\right)\cdot f\left(x\right)}-\ln\left(d\right)=\lim_{x\to 0}-\frac{\ln\left(kC\right)}{\ln\left(2\right)\cdot f\left(x\right)}-\frac{\ln\left(kC\right)}{\ln\left(2-\frac{k}{1+d}\right)\cdot f\left(x\right)}=0$$
Where :
$$\ln\left(\frac{k+1}{-\left(1-f\left(\frac{1}{2}\right)\right)}\right)=2,d=0.827...$$
Edit :12/01/2023
Let $d+c=k$ as above then define :
$$h\left(x\right)=\left(\frac{\ln\left(C\right)}{f\left(\frac{1}{2}\right)}-\frac{\ln\left(pC\right)}{\ln\left(\frac{p}{d}\right)\cdot f\left(x\right)}-\ln\left(d\right)\right)$$
Where $p$ is the Catalan's constant plus one or $p=1.915...$
Now define :
$$r\left(x\right)=\sum_{n=1}^{\infty}\frac{1}{e^{x^{2}n}+1+\ln\left(1+e\right)}$$
Then we have :
$$r(c)\simeq h(c)$$
See also the derivation due to Wolfram Alpha . See also the q-PolygammaFunction