Here's a proof of the statement for $f=0$, so that $X=W$ is a Wiener process. (I will give a more completeThe proof of the case with general $f$ further belowis a bit more involved, onceand I have written it outgive this further below). I'll base the proof on the following simple result. Here, I am using $B_\epsilon=\left\{\omega\in\mathcal{W}^n\colon\sup_{t\in[0,1]}\lVert\omega(t)\rVert\le\epsilon\right\}$ for the $\epsilon$-ball in Wiener space.
Here's a proof of the general case, based on the following lemma which generalizes Lemma 1 above (I have aa proof of this, and will update the answer once I have written it outis given further below).
Lemma 2: Suppose that $X$ satisfies an SDE of the form $dX=g(X_t,t)dt+dW$, $X_0=0$, where $g\colon\mathbb{R}^n\times[0,1]\to\mathbb{R}^n$. is (jointly) continuous and $g(x,t)$ has (jointly) continuous first and second order derivatives with respect to $x$. Then, for any continuous path $\phi\colon[0,1]\to\mathbb{R}^d$, $$ \limsup_{\epsilon\to0}\frac{\mathbb{P}(X\in B_\epsilon+\phi)}{\mathbb{P}(W\in B_\epsilon)}\le\exp\left(\frac12\sup_t\lVert Df(\phi(t),t)\rVert+c\int_0^1\lVert D^2g(\phi(t),t)\rVert dt\right)$$$$ \limsup_{\epsilon\to0}\frac{\mathbb{P}(X\in B_\epsilon+\phi)}{\mathbb{P}(W\in B_\epsilon)}\le K$$ where $c$$K$ is some constant (depending on $n$, but not onincreasing function of $g$ or the path$\lVert Dg(\phi,\cdot)\rVert$ and $\phi$)$\lVert D^2g(\phi,\cdot)\rVert$.
Let's use this to prove the result in the question. For any smooth $\gamma\colon[0,1]\to\mathbb{R}^n$, define $Y$ by $Y_0=0$ and $dY=f(Y+\gamma)dt+dW$. We see that $X-\gamma$ satisfies the same SDE as $Y$, but with $W-\gamma$ in place of $W$. Applying a Girsanov transformation as above, $$ \frac{\mathbb{P}(X\in B_\epsilon+\phi)}{\mathbb{P}(Y\in B_\epsilon+\phi-\gamma)}=\frac{\mathbb{P}(X-\gamma\in B_\epsilon+\phi-\gamma)}{\mathbb{P}(Y\in B_\epsilon+\phi-\gamma)}=\mathbb{E}\left[\exp\left(-\int_0^1\left(\dot\gamma dW_t+\frac12\dot\gamma^2dt\right)\right)\Bigg\vert Y\in B_\epsilon+\phi-\gamma\right]=\mathbb{E}\left[\exp\left(-\int_0^1\left(\dot\gamma dY_t-\dot\gamma f(Y_t+\gamma)dt+\frac12\dot\gamma^2dt\right)\right)\Bigg\vert Y\in B_\epsilon+\phi-\gamma\right]. $$ Taking limits as $\epsilon$ goes to 0 and applying Lemma 2, $$ \limsup_{\epsilon\to0} \frac{\mathbb{P}(X\in B_\epsilon+\phi)}{\mathbb{P}(W\in B_\epsilon)}\le K\exp\left(\int_0^1\left(\frac12\dot\gamma^2-\dot\gamma\dot\phi-\dot\gamma f(\phi)\right)dt\right). $$ As above, we use integration by parts to take the limit, and $\dot\phi$ is the derivative in the sense of distributions. The term $K$ is the right hand side of the inequality in Lemma 2, evaluated with $g(x,t)=f(x+\gamma_t,t)$ evaluated along the path $\phi-\gamma$. This is the same as the right hand side of the inequality evaluated with $g(x,t)=f(x)$ along the path $\phi$, so is independent of $\gamma$. Letting $\gamma$ be smooth approximations to $\phi$, the right hand side of the above inequality can be made arbitrarily close to 0, giving the required result.
Proof of Lemma 2:
Note that the SDE can be expressed as $W_t=X_t-\int_0^tg(X_s,s)ds$. Let us define the path $\psi_t=\phi_t-\int_0^tg(\phi_s,s)ds$. Expanding $g(X_s,s)$ as a power series to second order in $X-\phi$ gives $$ W_t-\psi_t=X_t-\phi_t-\int_0^t\left(A_s(X_s-\phi_s)+u_s\right)ds, $$ where $A_s=Dg(\phi_s,s)$ and $u_s$ is of second order in $X_s-\phi_s$. Choosing a fixed $L > \lVert D^2g(\phi_t,t)\rVert/2$ then $\lVert u_s\rVert\le L\lVert X_s-\phi_s\rVert^2\le L\epsilon^2$ for $X\in B_\epsilon+\phi$ and small enough $\epsilon$. Now, define $S$ to be set of all paths $\omega$ of the form $$ \omega_t=y_t-\int_0^t(A_sy_s+v_s)ds $$ for some $y\in B_\epsilon$ and $v\in B_{L\epsilon^2}$. From the above, we have $W\in S+\psi$ whenever $X\in B_\epsilon+\phi$. As $S$ is convex and symmetric, using log-concavity as in the proof of Lemma 1, $\mathbb{P}(W\in S+\psi)\le\mathbb{P}(W\in S)$.
Now, define a process $Y$ by $$ dY_t=dW_t+AY_t,\ \ Y_0=0. $$ Suppose that $\omega\in S$, $y\in B_\epsilon$ and $v\in B_{L\epsilon^2}$ are as above and $W_t=\omega_t$ then, we have $$ d\lVert Y_t-y_t\rVert/dt\le\lVert A(Y_t-y_t)-v_t\rVert\le a\lVert Y_t-y_t\rVert+L\epsilon^2 $$ where $a=\sup_t\lVert A_t\rVert$. From this we can bound $Y_t-y_t$ by $L\epsilon^2a^{-1}(e^{at}-1)$, which I will write as $M\epsilon^2$. So, $Y\in B_{\epsilon+M\epsilon^2}$, and we have $\mathbb{P}(W\in S)\le\mathbb{P}(Y\in B_{\epsilon+M\epsilon^2})$.
Next, define the process $$ dZ = dW - aZdt,\ \ Z_0=0. $$ This has drift $-a\lVert Z\rVert$ in the radial direction, whereas $Y$ has radial drift bounded below by $-a\lVert Y\rVert$. By comparing these processes we have $\mathbb{P}(Y\in B_{\epsilon+M\epsilon^2})\le\mathbb{P}(Z\in B_{\epsilon+M\epsilon^2})$.
Now, applying the Onsager--Machlup functional to the simple case with $f(x)=-ax$ and $\phi=0$ gives $$ \lim_{\epsilon\to0}\frac{\mathbb{P}(Z\in B_{\epsilon+M\epsilon^2})}{\mathbb{P}(W\in B_{\epsilon+M\epsilon^2})}=\exp(a). $$ We can also use the standard result that $\mathbb{P}(W\in B_\epsilon)\sim\exp(-c/\epsilon)$ in the limit $\epsilon\to0$ (where $c$ depends only on the dimension $n$) to get $$ \lim_{\epsilon\to0}\frac{\mathbb{P}(W\in B_{\epsilon+M\epsilon^2})}{\mathbb{P}(W\in B_{\epsilon})}=\exp(cM). $$
Putting the above results together gives the following sequence of inequalities, which should be understood in the asymptotic limit $\epsilon\to0$. $$ \begin{align} \frac{\mathbb{P}(X\in B_\epsilon+\phi)}{\mathbb{P}(W\in B_\epsilon)} &\le \frac{\mathbb{P}(W\in S+\psi)}{\mathbb{P}(W\in B_\epsilon)} \le \frac{\mathbb{P}(W\in S)}{\mathbb{P}(W\in B_\epsilon)}\\ &\le \frac{\mathbb{P}(Y\in B_{\epsilon+M\epsilon^2})}{\mathbb{P}(W\in B_\epsilon)} \le \frac{\mathbb{P}(Z\in B_{\epsilon+M\epsilon^2})}{\mathbb{P}(W\in B_\epsilon)}\\ &\le \exp(a)\frac{\mathbb{P}(W\in B_{\epsilon+M\epsilon^2})}{\mathbb{P}(W\in B_\epsilon)} \le\exp(a+cM). \end{align} $$ QED