I can't see a good way to answer your question without some more data.
I tried fitting a quadratic and cubic to your distance time data, and both give pretty good fits, but the fitted velocity at the 400m point is around 8m/sec so it's a big extrapolation down to 2m/sec. In fact the quadratic and cubic fits diverge quite strongly beyond about 450m so I wouldn't trust either of them.
I think the only good way to estimate the time to 2m/sec is if you know a specific mathematical model underlies the motion. If this is coursework there may be some hints about what model they expect you to use. If it's a real life example I think you'd be bold to assume any particular model. As I recall, resistance tends to go as $v^2$ at high speeds and switch to $v$ at low speeds, with an additional constant (i.e. independant of velocity) resistance due to friction in the drive train. But the details will depend on the car.
If you want the fits (fitted using Excel) they are:
$$ s = -0.1148t^2 + 15.596t + 1.4563 $$
$$ s = 0.0017t^3 - 0.201t^2 + 16.64t + 0.0021 $$
Later: I'm probably taking this further than the data supports, but if you assume the deceleration is proportional to $v^2$ you get $v \propto 1/t$ and hence $s \propto ln(t)$. If I fit a function of this form to your data I find the following gives a pretty good fit:
$$ s = 577.5 \space ln(t + 34.1) + 2038 $$
The point of this fit is that it's physically motivated while the fits above were purely numerical. Anyhow, differentiating this to get the velocity gives:
$$ v = \frac{577.5}{t + 34.1} $$
So if you're prepared to assume the deceleration stays proportional to $v^2$, the car reaches 2 m/sec after 255 seconds.
Response to comment: I can only speak as an (ex) experimental physicist, but faced with some data the first thing I would do is graph it and see what the curve looks like. Experience often immediately tells you what the functional form is e.g. whether it's a parabola, hyperbola, exponential etc. However in this case nothing obvious sprang to mind.
The next step is to fit a few random functions and see if anything looks plausible. Excel has a polynomial fit built in. Graph your data as a scatter plot then right click the curve in the graph and select "add trendline". Select the type of fit, then on the options tab click "include formula" or you'll just get the line and not the fitted equation.
In this case the polynomial fit didn't shed any light, so you have to get out the pencil and paper and try to work out a physically reasonable description. Wikipedia will tell you that at high speed air resistance goes as $v^2$, so a reasonable guess is that the retarding force is proportional to $v^2$. I used my polynomial fits to calculate $v$ and $a$ from your data. With only five points the results were pretty noisy, but it did look as if a 2-fold change in velocity did cause a 4-fold change in acceleration so the $v^2$ dependance looked plausible.
So the equation of motion is something like:
$$ \frac{dv}{dt} = -A \space v^2 $$
for some constant, $A$, and the solution to this is:
$$ v = \frac{B}{t} $$
for some other constant $B$, i.e. velocity falls of as 1/time. This function for $v$ says that $v$ takes infinite time to fall to zero, but we expect it only applies at high speeds so that's OK. Since $v = ds/dt$ we integrate the expression for $v$ to give the logarithmic dependance of distance on time.
The equation I used:
$$ s = a \space ln(t + b) + c $$
is just the most general form of the equation for $s$. I went back to excel and graphed the data along with an extra column for the log equation, then I manually adjusted the values of $a$, $b$ and $c$ to get a rough fit.
At this point I used an Excel add-on called the Solver. I won't go into how to use this because it's quite fiddly and there are lots of articles on how to use it out there in Googleland. The solver can adjust variables you specify to get the best fit of your function to data, and that gave me the equation I posted.
On last comment, extrapolating data is always dangerous because you have to assume that whatever equation you use applies outside the range of data you've collected. The velocity in your data falls from about 16 m/sec at $t = 0$ to about 8 m/sec at the last point in your data, and it's a large extrapolation to try and extend this to 2 m/sec. I doubt that the $v^2$ dependance of deceleration would apply down to 2 m/sec so my calculated value for the time taken to decelerate to 2 m/sec is probably much too high.
Best Answer
If we apply a force $F$ to a mass $m$ and a friction force (drag) $F_d$ also acts on it the force diagram becomes:
With $a$ the acceleration the object experiences, the equation of motion becomes:
$F=ma+F_d$.
As the mass moves towards the right, say for an infinitesimal distance $dx$, an infinitesimal amount of work $dW$ is performed on $m$ by $F$:
$dW=(ma+F_d)dx$.
If we divide both sides with $dt$ then $\frac{dW}{dt}=P$ with $P$ the power, constant in this case. So we have:
$P=\frac{dW}{dt}=(ma+F_d)\frac{dx}{dt}$ and by definition $\frac{dx}{dt}=v$, so:
$P=(ma+F_d)v$.
Now the question becomes, what is $F_d$?
We know that in general $F_d \propto v^n$, where $v$ is the velocity and $n$ is some exponent. For instance in the case Navier Stokes drag (viscous drag of a fluid on a spherical object), $n=1$.
For the case of air drag the exponent $n=2$ is generally assumed.
Let's however explore the case of $n=1$, so that $F_d=kv$, with $k$ a proportionality constant, so we get:
$P=(ma+kv)v$.
With $a=\frac{dv}{dt}$, we get:
$P=mv\frac{dv}{dt}+kv^2$, a differential equation that can be separated by variables to yield:
$m \frac{v}{P-kv^2}dv=dt$.
This can be integrated between $t=0, v=0$ and $t, v$ and yields after reworking:
$\Large{v=\sqrt{\frac{P}{k}(1-e^{-\frac{2kt}{m}})}}$.
For $t \to +\infty$ the exponential term $e^{-\frac{2kt}{m}} \to 0$, so that the terminal velocity $v_t$ is achieved for $t = +\infty$ and is given by:
$\large{v_t=\sqrt{\frac{P}{k}}}$.
Since as the terminal speed is only achieved for $t = \infty$, it will be achieved also only for $x = \infty$.
The general shape of the $(v,t)$ function is as follows, with $v_t$ being reached only asymptotically: