I do not like this proof, for reasons that will become obvious.
By shifting the summation index and using a binomial ID in the denominator, we can write the LHS (left-hand side) of the proposed identity as
$$D_n:=\sum_{k=1}^{n-1} \frac{ (k-1)\binom{n-k+1}{k} + (k+1)\binom{n-k}{k} }{n\binom{n-1}{k}} $$
For the denominator, use the ID (beta integral identity)
$$\frac{1}{n\binom{n-1}{k}}=\int_0^1t^k(1-t)^{n-1-k} dt = 2\int_0^{\pi/2}
\cos^{2k+1}(t) \ \sin^{2n-2k-1}(t) dt $$
This is done so that the sum can be solved in explicit form in terms of the Chebyshev polynomials of first and second kind. e.g.,
$$\sum_{k=0}^n (-1)^kx^{2k}\binom{n-k}{k} = U_n\big(\frac{1}{2x}\big)x^n$$
With one integration by parts and simplification, we can write an equivalent expression, with $\omega_n(t)=(n\cos^{2n-1}(t)-(n+1)\cos^{2n+1}(t))\sin{(t)},$
$$ (eq.1) \quad \frac{1}{2}D_n = \int_{0}^{\pi/2} \omega_n(t)(i\tan{t})^{n+1}U_{n+1}\big(\frac{1}{2i\ \tan{t}} \big) dt \quad+ $$
$$+ \quad \int_{0}^{\pi/2} \omega_n(t)(i\tan{t})^{n}U_{n}\big(\frac{1}{2i\ \tan{t}} \big) dt
-\frac{2}{n+1}\int_{0}^{\pi/2} \omega_n(t)(i\tan{t})^{n+1}T_{n+1}\big(\frac{1}{2i\ \tan{t}} \big) dt$$
At this point I tried to use recurrence formulas for the Chebyshev polynomials to prove $D_{n+1} = D_n,$ but I was unsuccessful. Therefore, I did the following: create a generating function by operating on both sides with $\sum_{n=1}^\infty x^n/n $. Now by explicit calculation, $D_1 = 0.$ Therefore we want to verify
$$ (eq.2)\quad -(x+\log{(1-x)})/2 = \sum_{n=1}^\infty x^n/n \ RHS(*) $$
To do this, you need six explicit sums, which are not difficult to derive:
$$\sum_{n=1}^\infty \frac{x^n}{n} U_{n}(y) = \frac{y}{\sqrt{1-y^2}} \text{arctan}\big(x\frac{\sqrt{1-y^2}}{1-xy}\big)-\frac{1}{2}\log{(1+x^2-2xy)} $$
$$\sum_{n=1}^\infty \frac{x^n}{n} U_{n+1}(y) = \frac{2y^2-1}{\sqrt{1-y^2}} \text{arctan}\big(x\frac{\sqrt{1-y^2}}{1-xy}\big)-y\log{(1+x^2-2xy)} $$
$$\sum_{n=1}^\infty x^n U_{n}(y) = \frac{1}{1-2xy+x^2}-1 $$
$$\sum_{n=1}^\infty x^n U_{n+1}(y) =\frac{1}{x}\Big( \frac{1}{1-2xy+x^2}-1 -2xy\Big)$$
$$\sum_{n=1}^\infty \frac{x^n}{n+1} T_{n+1}(y) = \frac{1}{x}\Big(-\frac{1}{2}\log{(1+x^2-2xy)} - xy\Big) $$
$$\sum_{n=1}^\infty \frac{x^n}{n} T_{n+1}(y) = -\sqrt{1-y^2} \text{arctan}\big(x\frac{\sqrt{1-y^2}}{1-xy}\big)-\frac{y}{2}\log{(1+x^2-2xy)} $$
With these closed forms, the indefinite integrals are solvable in closed form by Mathematica. There is considerable simplification needed to finally get the form on the LHS of (eq.2).
I detest this proof because it does not give a hint at all about how the identity was discovered, and the huge amounts of simplification that must be performed with a CAS and human intervention.