It looks like there is a mistake somewhere.
If this were true,
we would have
$$xf(x) = -f'(x)F(x)^2 - 2f(x)^2F(x) \quad, \forall x \in \mathbb{R}$$
Let's see if $f(x) = c e^{-x^2/2}$ satisfies this (with $c = \frac{1}{\sqrt{2\pi}}$). Substituting $x=0$ gives:
$$ 0 = -f'(0)F(0)^2 - 2f(0)^2F(0) $$
But $F(0)=1/2$, $f'(0)=0$, $f(0)=c$, so we get
$$ 0 = -c^2$$
So $c=0$, a contradiction.
I think the given information is stated incorrectly. It should be stated as
$$ \int_{-\infty}^b x f(x)dx = -f(b) \quad, \forall b \in \mathbb{R} $$
This yields an ODE that is solved by a Gaussian. This modified information is like saying that the mean of $Y_b$ is equal to $-f(b)/F(b)$ (rather than $-f(b)F(b)$) for all $b$.