This does not answer the question. I don't prove anything here but just present the result of some simulations that suggest the answer to the question is yes.
My experience make me being careful about my intuition when it is about to think about geometrical construction in more then $3$ dimensions. So before investigating more this question I wanted to check if there is no counter example I could get with Mathematica.
I generated some random simplices with a random point inside and counted the number of balls covering the inside point. For all the tests I did, no counter example appeared. Moreover, the distribution of the number of balls in my simulations seem to be close to a normal distribution.
For example, I generated $1000$ random $3$-dimensional simplices as the convex hull of $4$ random points uniformly distributed in the unit cube. For each simplex I generated $1000$ random inside points with baricenter coordinate uniformly distributed.
Below this is the histogram of the number of balls covering the inside points considered.
We observe that it starts at $3$, which means we did not find counter example to the conjecture.
dimension$=3$; number of edges$={4 \choose 2}=6$.

Below I give more histograms obtained from simulations in higher dimension. For each I generated $100$ random simplices and and tested $100$ random inside points.
dimension$=4$; number of edges$={5 \choose 2}=10$.

dimension$=5$; number of edges$={6 \choose 2}=15$.

dimension$=6$; number of edges$={7 \choose 2}=21$.

dimension$=7$; number of edges$={8 \choose 2}=28$.

dimension$=8$; number of edges$={9 \choose 2}=36$.

dimension$=20$; number of edges$={20 \choose 2}=190$.

The code I used with Mathematica (it can certainly be optimized if you wanna run biger simulations)
> RandomSimplex[dim_] := RandomReal[{0, 1}, {dim, dim + 1}]
> RandomBaricenter[dimension_] := Module[{baricenter},
> baricenter = RandomReal[{0, 1}, dimension + 1];
> baricenter = baricenter/Sum[baricenter[[i]], {i, 1, dimension + 1}];(*normalisation*)
> Return[baricenter]]
> Distance[A_, B_] := Module[{dim},
> dim = Dimensions[A][[1]];
> Return[Sqrt[Sum[(B[[i]] - A[[i]])^2, {i, 1, dim}]]];]
> NumberClosePoints[simplex_, point_] := Module[{dim, number, i, j, M, A, B},
> dim = Dimensions[simplex][[1]];
> number = 0;
> For[i = 1, i <= dim, i++,
> For[j = i + 1, j <= dim + 1, j++,
> A = Transpose[simplex][[i]];
> B = Transpose[simplex][[j]];
> M = (A + B)/2;
> If[Distance[M, A] > Distance[M, point], number++,]]];
> Return[number]]
and to run the tests and produce and histogram:
> dimension = 3;
> nBar = 1000;(*number of barycenter tested*)
> nSimplices = 1000;(*number of simplices tested*)
> list = {};
> For[i = 1, i <= nBar, i++,
> bar = RandomBaricenter[dimension];
> For[j = 1, j <= nSimplices, j++,
> simplex = RandomSimplex[dimension];
> point = simplex.bar;
> n = NumberClosePoints[simplex, point];
> list = Append[list, n];
> If[n < dimension, Print[{bar, simplex, n}],]];]
> Histogram[list]