Just for sake of sharing, I add my elementary idea: For any $f\in C^\infty_0(\mathbb{R})$ (for more general function it is a common practice to define via distribution $\langle (-\Delta)^sf,g\rangle:=\langle f,(-\Delta)^s g\rangle$ for all $g\in C_0^\infty(\mathbb{R})$ with some given conditions on $f$), we have
$$
\begin{align*}
&(-\Delta)^sf(0)\\
&=c(s)P.V.\int_{-\infty}^\infty\frac{f(0)-f(x)}{|x|^{1+2s}}\,dx\\
&=c(s)P.V.\int_{B_\epsilon(0)}\frac{f(0)-f(x)}{|x|^{1+2s}}\,dx\\
&+c(s)\int_{B_\epsilon^c(0)}\frac{f(0)-f(x)}{|x|^{1+2s}}\,dx\\
&=c(s)P.V.\int_{B_\epsilon(0)}-f'(0)\textrm{sign}(x)|x|^{-2s}-\frac{f''(0)}{2}|x|^{1-2s}+O(|x|^{2-2s})\,dx+c(s)C_{\epsilon,f}O(1)\\
&=-f''(0)\frac{c(s)\epsilon^{2-2s}}{2-2s}+c(s)O(\epsilon^{4-2s})+c(s)C_{\epsilon,f}O(1)\\
&\xrightarrow{s\rightarrow 1} -f''(0),
\end{align*}
$$
because $c(s)$ has order one zero at $s=1$, i.e. $c(s)\sim 2(1-s)$, and the second and third term just vanishs in the limit of $s\rightarrow 1-$. I did not verify that $\lim_{s\rightarrow 1-}\frac{c(s)\epsilon^{2-2s}}{2-2s}=1$ with the precise definition of $c(s)$, but I guess it should be true if $c_s$ is choosen properly. So we just have the classical Laplacian as we hoped.
EDIT: For the completeness I add the interpretation for $(-\Delta)^{s}\xrightarrow{s\rightarrow 0+}Id$ in the limit $s\rightarrow 0+$ as well. In this case all terms vanish except the one multiplied with $f(0)$.
$$
\begin{align*}
&(-\Delta)^sf(0)=c(s)P.V.\int_{B_\epsilon(0)}\frac{f(0)-f(x)}{|x|^{1+2s}}\,dx\\
&+c(s)\int_{B_\epsilon^c(0)}\frac{f(0)-f(x)}{|x|^{1+2s}}\,dx\\
&=c(s)C_{\epsilon,f}O(1)+f(0)\frac{c(s)\epsilon^{-2s}}{2s}\\
&\xrightarrow{s\rightarrow 0+}f(0),
\end{align*}
$$
because $c(s)\sim 2s$ near $s\sim 0$. Note that we have used $\int_{B_\epsilon^c(0)}\frac{f(x)}{|x|^{1+2s}}=C_{\epsilon,f}O(1)$.