An assortment of curves for fitting chemistry examples is presented in these Colby College class notes. Of particular application is the sigmoid response curve with variable "slope" for the central part of the curve:
$$ f(x) = \frac{a}{1 + e^{bx - c} } + d $$
[This is similar to the suggested logistic function proposed in the first Answer, but has four rather than three parameters allowing for the inflection points (changing sign of second derivative) to appear in a location other than at the origin.]
Such functions are inherently monotone (decreasing when $b\gt 0$ and increasing when $b\lt 0$), so they will "automatically" possess that property if you fit your data to one of these models, avoiding the tendency of polynomial curve fits to oscillations seen in your graph.
A disadvantage of this nonlinear parameterization is that computing it requires an iterative process, rather than by the direct solvers available for a linear least squares fitting as you seem to have used for the fifth degree polynomial.
But nonlinear least squares fitting is often not difficult when reasonable initial estimates for the parameters are available. There are online solvers (at a glance Nonlinear Least Squares Regression suggested in the class note above seems to be intact, although a second site mentioned there returns 404:Not found). Most spreadsheets have either built-in or add-in solvers that would be capable of doing the fit. But you may find it edifying to do the fit "manually" so that you understand the roles played by the parameters.
I would start with the two extremes where the curve levels off. That is, for arguments well below the "drop-off" portion of the curve, the data appears to approach $f(0) = 360$ more or less. From this one might estimate:
$$ a + d \approx 360 $$
while the data well above the "drop-off" suggests $d \approx 130$. So we can get initial estimates for parameters $a,d$ easily.
Visually the inflection point (change from concave down to concave up) appears to occur at about $x=3.5$. Solving $f''(x) = 0$ gives us $bx = c$, thus eliminating one more parameter estimate. Finally the first derivative of $f(x)$ at the inflection point can be estimated visually from the slope exhibited by the data (although this seems rather steep, circa $600$, so less firmly estimated by sight).
The iterative procedure (nonlinear least squares regression) will seek to adjust the parameters $a,b,c,d$ to minimize the least squares measure of error:
$$ \sum_{i=1}^n (y_i - f(x_i))^2 $$
where summation is taken over the set of data $(x_i,y_i)$ from your titration experiments.