Let $B(a,b;x):=\int_0^x t^{a-1}(1-t)^{b-1}dt$ denote, as usual, the value of the incomplete Beta function. We have to show that $$ n\binom{n}{k}B(k+1,n-k+1,\tfrac k{n+1}) \le 1/2 \quad\text{for $k=0,\dots,n$}. \tag{1} $$ As suggested by Brendan McKay, it should be easy to use asymptotic analysis to check that $(1)$ holds for any real $c>1/2$ if $n$ is large enough. Thus, we can use the reverse induction in $n$. That is, it is enough to show that, if $(1)$ holds for some natural $n+1\ge2$ in place of $n$ and some real $c>0$ in place of $1/2$, then it holds as written, but with the same real $c>0$ in place of $1/2$. The key here is the simple identity \begin{equation} B(a,b;x)=B(a+1,b;x)+B(a,b+1;x), \end{equation} which yields, by the reverse induction, $$ B(k+1,n-k+1,\tfrac k{n+1})=B(k+2,n-k+1,\tfrac k{n+1})+B(k+1,n-k+2,\tfrac k{n+1}) $$ $$ =B(k+2,n-k+1,\tfrac{k+1}{n+2})+B(k+1,n-k+2,\tfrac k{n+2})+I_1-I_2 $$ $$ \le\frac c{n+1}/\binom{n+1}{k+1} +\frac c{n+1}/\binom{n+1}{k}+I_1-I_2 =\frac{c(n+2)n}{(n+1)^2}\frac1{n\binom nk}+I_1-I_2 <\frac c{n\binom nk}+I_1-I_2, $$ where \begin{equation} I_1:=\int_{t_0}^{t_1}f(t)\,\frac{dt}t,\quad I_2:=\int_{t_1}^{t_2}f(t)\,\frac{dt}{1-t}, \end{equation} \begin{equation} f(t):=t^{k+1}(1-t)^{n-k+1},\quad t_0:=\tfrac k{n+2}, \quad t_1:=\tfrac k{n+1}, \quad t_2:=\tfrac{k+1}{n+2}, \end{equation} so that $0\le t_0<t_1<t_2<1$ and the function $f$ is increasing on the interval $(t_0,t_2)$. It remains to check the inequality $I_1\le I_2$. This is easy: since $f$ is increasing, \begin{equation} I_1-I_2\le \int_{t_0}^{t_1}f(t_1)\,\frac{dt}t-\int_{t_1}^{t_2}f(t_1)\,\frac{dt}{1-t}=0. \end{equation} $\qed$