Non-normal distributions with zero skewness and zero excess kurtosis? Mostly theoretical question. Are there any examples of non-normal distributions that has first four moment equal to those of normal? Could they exist in theory?
 A: Yes, examples with skewness and excess kurtosis both zero are relatively easy to construct. (Indeed examples (a) to (d) below also have Pearson mean-median skewness 0)
(a) For example, in this answer an example is given by taking a 50-50 mixture of a gamma variate, (which I call $X$), and the negative of a second one, which has a density that looks like this:

Clearly the result is symmetric and not normal. The scale parameter is unimportant here, so we can make it 1. Careful choice of the shape parameter of the gamma yields the required kurtosis:

*

*The variance of this double-gamma ($Y$) is easy to work out in terms of the gamma variate it's based on: $\text{Var}(Y)=E(X^2)=\text{Var}(X)+E(X)^2=\alpha+\alpha^2$.


*The fourth central moment of the variable $Y$ is the same as $E(X^4)$, which for a gamma($\alpha$) is $\alpha(\alpha+1)(\alpha+2)(\alpha+3)$
As a result the kurtosis is $\frac{\alpha(\alpha+1)(\alpha+2)(\alpha+3)}{\alpha^2(\alpha+1)^2}=\frac{(\alpha+2)(\alpha+3)}{\alpha(\alpha+1)}$. This is $3$ when $(\alpha+2)(\alpha+3)=3\alpha(\alpha+1)$, which happens when $\alpha=(\sqrt{13}+1)/2\approx 2.303$.

(b) We could also create an example as a scale mixture of two uniforms. Let $U_1\sim U(-1,1)$ and let $U_2\sim U(-a,a)$, and let $M=\frac12 U_1+\frac12 U_2$. Clearly by considering that $M$ is symmetric and has finite range, we must have $E(M)=0$; the skewness will also be 0 and central moments and raw moments will be the same.
$\text{Var}(M)=E(M^2)=\frac12\text{Var}(U1)+\frac12\text{Var}(U_2)=\frac16[1+a^2]$.
Similarly, $E(M^4)=\frac{1}{10} (1+a^4)$ and so
the kurtosis is $\frac{\frac{1}{10} (1+a^4)}{[\frac16 (1+a^2)]^2}=3.6\frac{1+a^4}{(1+a^2)^2}$
If we choose $a=\sqrt{5+\sqrt{24}}\approx 3.1463$, then kurtosis is 3, and the density looks like this:


(c) here's a fun example. Let $X_i\stackrel{_\text{iid}}{\sim}\text{Pois}(\lambda)$, for $i=1,2$.
Let $Y$ be a 50-50 mixture of $\sqrt{X_1}$ and $-\sqrt{X_2}$:

by symmetry $E(Y)=0$ (we also need $E(|Y|)$ to be finite but given $E(X_1)$ is finite, we have that)
$Var(Y)=E(Y^2)=E(X_1)=\lambda$
by symmetry (and the fact that the absolute 3rd moment exists) skew=0
4th moment: $E(Y^4) = E(X_1^2) = \lambda+\lambda^2$
kurtosis = $\frac{\lambda+\lambda^2}{\lambda^2}= 1+1/\lambda$
so when $\lambda=\frac12$, kurtosis is 3. This is the case illustrated above.

(d) all my examples so far have been symmetric, since symmetric answers are easier to create -- but asymmetric solutions are also possible. Here's a discrete example.

(e) Now, here's an asymmetric continuous family. It will perhaps be the most surprising for some readers, so I'll describe it in detail. I'll begin by describing a discrete example and then build a continuous example from it (indeed I could have started with the one in (d), and it would have been simpler to play with, but I didn't, so we also have another discrete example for free).
$\:\,$ (i) At $x=-2,1$ and $m=\frac12 (5+\sqrt{33})$ ($\approx 5.3723$) place probabilities of $p_{-2}= \frac{1}{36}(7+\sqrt{33})$, $p_1=\frac{1}{36}(17+\sqrt{33})$, and $p_m=\frac{1}{36}(12-2\sqrt{33})$ (approximately 35.402%, 63.179% and 1.419%), respectively. This asymmetric three-point discrete distribution has zero skewness and zero excess kurtosis (as with all the above examples, it also has mean zero, which simplifies the calculations).
$\:$ (ii) Now, let's make a continuous mixture. Centered at each of the ordinates above (-2,1,m), place a Gaussian kernel with common standard deviation $\sigma$, and probability-weight given by the probabilities above (i.e. $w=(p_{-2},p_1,p_m)$). Phrased another way, take a mixture of three Gaussians with means at $-2,1$ and $m$ each with standard deviation $\sigma$ in the proportions $(p_{-2},p_1,p_m)$ respectively. For any choice of $\sigma$ the resulting continuous distribution has skewness 0 and excess kurtosis 0.
Here's one example (here the common $\sigma$ for the normal components is 1.25):

(The marks below the density show the locations of the centers of the Gaussian components.)

As you see, none of these examples look particularly "normal". It would be a simple matter to make any number of discrete, continuous or mixed variables with the same properties. While most of my examples were constructed as mixtures, there's nothing special about mixtures, other than they're often a convenient way to make distributions with properties the way you want, a bit like building things with Lego.
This answer gives some additional details on kurtosis that should make some of the considerations involved in constructing other examples a little clearer.

You could match more moments in similar fashion, though it requires more effort to do so. However, because the MGF of the normal exists, you can't match all integer moments of a normal with some non-normal distribution, since that would mean their MGFs match, implying the second distribution was normal as well.
A: Good points are made by Glen_b. I would only add consideration of the Dirac Delta function as additional grist for the mill. As Wikipedia notes, "The DDF is a generalized function, or distribution, on the real number line that is zero everywhere except at zero, with an integral of one over the entire real line" with the consequence that all higher moments of the DDF are zero. 
Paul Dirac applies it to quantum mechanics in his 1931 book The Principles of Quantum Mechanics but it's origins date back to Fourier, Lesbesgue, Cauchy and others. The DDF also has physical analogues in modeling the distribution, e.g., of the crack of a bat hitting a baseball. 
