$$f(x)=x\left(\sin^2(x)+\frac1{1+x^2}\right)$$
- There is only one value of $x$ such that $f(x)=0$
- $f'(0)=1\ne0$
- $f$ oscillates between values arbitrarily close to $0$ and arbitrarily large values.
The last property guarantees that $f^{-1}((-\varepsilon,\varepsilon))$ is made of a union of infinitely many disjoint open intervals, hence it's not homeomorphic to $(-\varepsilon,\varepsilon)$, for any $\varepsilon>0$.

With the example in my comment above, $f(x)=\dfrac{x}{1+x^2}$, the constraints where only partially fulfilled: for a large enough $\varepsilon>0$, $f^{-1}((-\varepsilon,\varepsilon))=\Bbb R$, because that $f$ is bounded. Here it doesn't happen.