First of all, periodic boundary conditions look like a pair of boundary conditions, e.g. like $$u(-1,t)=u(1,t), \quad u_x(-1,t)=u_x(1,t).$$
Start by separating variables, finding the eigenvalues $\lambda_n$, eigenfunctions $X_n(x)$, and temporal solutions $T_n(t)$.
Superimpose those to get something of the form $u(x,t)=\sum_n c_n X_n(x)T_n(t)$.
To obtain the coefficients in the series solution above, evaluate $u(x,t)$ at $t=0$, apply the initial data, and recognize this as an appropriate Fourier series. Use your knowledge of Fourier series and/or orthogonality on $[-1,1]$ to deduce formulas for the coefficients.
Edit based on progress in the comments:
The superimposed solution should look like $$u(x,t)=\sum_{n=1}^\infty [a_n\cos(n\pi x)+b_n\sin(n\pi x)]e^{10(n\pi)^2 t}.$$
Then, the initial condition $u(x,0)=f(x)$, $-1<x<1$, leads to
$$u(x,0)=\underbrace{\sum_{n=1}^\infty a_n\cos(n\pi x)+b_n\sin(n\pi x)=f(x)}, \quad -1<x<1,$$
which can be found as standard Fourier coefficients.