Yes, there are such functions. Take a very narrow region $D$ containing the positive ray,
with nice boundary and such that $D$ intersects every any horizontal other than the real line by a bounded interval.
Let $g$ be a conformal map of $D'$ onto the right half-plane, where
$D'\subset D$ is another similar region. Then with appropriate choice of $D$ and $D'$
$$f(z)=\int_{\partial D} \frac{e^{g(\zeta)}}{\zeta-z}d\zeta$$

will converge for $z$ outside $D$ and the function $f$ will be bounded outside $D$ and extend to an entire function (by a deformation of the contour).

For the details, see for example  MR2753600 or  MR0545054.

EDIT. The method is very flexible of course. Taking $D$ to be a half-strip and $g(z)=e^z$ one obtains a Mittag-Leffler functiuon. Replacing it by $f(z+4i)$ you obtain a function that is bounded on every line from the origin.  But to obtain a function as you ask, a half-strip $D$ will not work, so the function is less elementary.