The vector operator $\vec\nabla^2\{\cdot\}$ can be expressed (and is defined in general coordinates) in $\mathbb{R^3}$ as
$$\begin{align}
\vec\nabla^2\{\cdot\}&=\text{Grad}(\text{Div} \{\cdot\})-\text{Curl}\left(\text{Curl}\{\cdot\}\right)\\\\
&=\nabla \left(\nabla\cdot\{\cdot\}\right)-\nabla\times\left(\nabla\times\{\cdot\}\right)
\end{align}$$
In Cartesian coordinates in $\mathbb{R}^N$, the vector Laplacian is given by
$$\begin{align}
\vec\nabla^2\{\cdot\}&=\sum_{i=1}^N\frac{\partial^2 \{\cdot\}}{\partial x_i^2}\\\\
\end{align}$$
So, we have for a vector field $\vec A$, in $\mathbb{R^N}$
$$\begin{align}
\nabla^2 \vec A(\vec r)&=\sum_{i=1}^N\sum_{j=1}^N\hat x_j \frac{\partial^2A_j(\vec r)}{\partial x_i^2} \\\\
&=\sum_{j=1}^N \hat x_j \nabla^2 A_j(\vec r)
\end{align}$$
For the specific problem, $\vec A(\vec r)=\vec E(\vec r)$, and $\vec E(\vec r)=\vec E_0 f(\vec k\cdot \vec r-\omega t)$, with $\vec r\in \mathbb{R}^3$. Therefore, the vector Laplacian applied to $\vec E(\vec r)$ simplifies to
$$\begin{align}
\vec \nabla^2\vec E(\vec r)&=\nabla \left(\nabla\cdot\vec E(\vec r)\right)-\nabla\times\left(\nabla\times\vec E(\vec r)\right)\\\\
&=\vec k(\vec k\cdot \vec E_0)f''(\vec k\cdot\vec r-\omega t)-\vec k\times(\vec k\times \vec E_0)f''(\vec k\cdot\vec r-\omega t)\\\\
&=|\vec k|^2f''(\vec k\cdot\vec r-\omega t)\vec E_0
\end{align}$$