This is an alternative proof to show that $\beta \le \sqrt n \max(\sigma_1,\dotsc,\sigma_n,\beta_1,\dotsc, \beta_n)$.
In the proof, we assume for simplicity that all $\mu_i=0$, as Bernstein's condition is only concerned with central moments. We use a shorthand $\beta_\max \triangleq \max \{\beta_1,\ldots,\beta_n, \sigma_1,\dotsc,\sigma_n\}$.
For $k=3$, the result can be checked in a straightforward manner. For $n\ge 2$ and $k \ge 4$, we have:
\begin{align*}
&\left| \mathbb E{\left(X_1+\dotsb+X_n\right)^k} \right| \\
&= \left| \sum_{\substack{i_1+\dotsb+i_n=k \\ i_1,\dotsc,i_n \ge 0}} \frac{k!}{i_1!\cdot \dots \cdot i_n!} \mathbb E{X_1^{i_1}} \cdot \dots \cdot \mathbb E{X_n^{i_n}}\right| \\
&\le \sum_{\substack{i_1+\dotsb+i_n=k \\ i_1,\dotsc,i_n \ge 0}} \frac{k!}{i_1!\cdot \dots \cdot i_n!} \prod_{j=1}^n \left|\mathbb E{X_j^{i_j}}\right|
\end{align*}
Note that $\mathbb E{X_j^0} = 1$ and $\mathbb E{X_j}=\mu_j=0$. We continue:
\begin{align*}
&\le \sum_{j=1}^n \left|\mathbb E{X_j^k}\right| + \sum_{\substack{i_1+\dotsb+i_n=k \\ i_1,\dotsc,i_n \neq 1,k-1,k}} \frac{k!}{i_1!\cdot \dots \cdot i_n!} \prod_{\substack{j=1\\i_j \neq 0}}^n \left|\mathbb E{X_j^{i_j}}\right| \\
&\le \sum_{j=1}^n \frac 12 k! \sigma_j^2 \beta_j^{k-2} \\
&\qquad + \sum_{\substack{i_1+\dotsb+i_n=k \\ i_1,\dotsc,i_n \neq 1,k-1,k}} \frac{k!}{i_1!\cdot \dots \cdot i_n!} \prod_{\substack{j=1\\i_j \neq 0}}^n \frac 12 i_j! \sigma_j^2 \beta_j^{i_j-2} \\
&\le \frac 12 k! \sigma^2 \beta_{\mathrm{max}}^{k-2} + k! \sum_{\substack{i_1+\dotsb+i_n=k \\ i_1,\dotsc,i_n \neq 1,k-1,k}} \frac 1{2^{w_H(i_1,\dotsc,i_n)}} \beta_{\mathrm{max}}^{k-w_H(i_1,\dotsc,i_n)} \prod_{\substack{j=1\\i_j \neq 0}}^n \sigma_j,
\end{align*}
where $w_H(i_1,\dotsc,i_n)$ denotes the number of non-zeros among $i_1,\dotsc,i_n$.
Observe that if $J = \{j \in [n] \mid i_j \neq 0\}$, where $|J| = w \ge 2$, then for any $j_1,j_2 \in J$, $j_1 < j_2$, it is true that
$$
\prod_{j \in J} \sigma_j \le \sigma_{j_1} \sigma_{j_2} \beta_{\mathrm{max}}^{w-2}.
$$
Summing up over the choice of ordered pairs $(j_1,j_2) \in J^2$, $j_1 < j_2$, we get:
$$
\prod_{j \in J} \sigma_j \le \frac{\beta_{\mathrm{max}}^{w-2}}{\binom w2} \sum_{\substack{j_1,j_2 \in J \\ j_1 < j_2}} \sigma_{j_1} \sigma_{j_2}
$$
We continue:
\begin{align*}
&\le \frac 12 k! \sigma^2 \beta_{\mathrm{max}}^{k-2}
+ k! \sum_{\substack{i_1+\dotsb+i_n=k \\ i_1,\dotsc,i_n \neq 1,k-1,k}} \frac 1{2^{w_H(i_1,\dotsc,i_n)}} \beta_{\mathrm{max}}^{k-w_H(i_1,\dotsc,i_n)} \prod_{\substack{j=1\\i_j \neq 0}}^n \sigma_j \\
&\le \frac 12 k! \sigma^2 \beta_{\mathrm{max}}^{k-2} \\
&\qquad + \frac 12 k! \beta_{\mathrm{max}}^{k-2} \sum_{\substack{i_1+\dotsb+i_n=k \\ i_1,\dotsc,i_n \neq 1,k-1,k}} \frac 1{2^{w_H(i_1,\dotsc,i_n)-1} \binom{w_H(i_1,\dotsc,w_n)}2} \sum_{\substack{1 \le j_1 < j_2 \le n \\ i_{j_1},i_{j_2} \neq 0}} \sigma_{j_1} \sigma_{j_2}
\end{align*}
Consider the following polynomial in $\binom n2$ variables $s_{j_1j_2}$, $j_1,j_2 \in [n]$, $j_1 < j_2$:
\begin{align}
&f(s_{12},\dotsc,s_{n-1,n}) \\
&\qquad= \sum_{\substack{i_1+\dotsb+i_n=k \\ i_1,\dotsc,i_n \neq 1,k-1,k}} \frac 1{2^{w_\mathrm{H}(i_1,\dotsc,i_n)-1} \binom{w_\mathrm{H}(i_1,\dotsc,i_n)}2} \sum_{\substack{1 \le j_1 < j_2 \le n \\ i_{j_1},i_{j_2} \neq 0}} s_{j_1j_2}
\tag{1}
\end{align}
It has degree $1$ (linear) and symmetric in its variables. Any degree-$1$ symmetric polynomial in $\binom n2$ variables has the form
\begin{equation}\tag{2}
f(s_{12},\dotsc,s_{n-1,n}) = a \frac{s_{12} + s_{13} + \dotsb + s_{n-1,n}}{\binom n2}.
\end{equation}
Setting all $s_{j_1j_2} = 1$, and comparing (1) with (2), we obtain that
$$
a = \sum_{\substack{i_1+\dotsb+i_n=k \\ i_1,\dotsc,i_n \neq 1,k-1,k}} \frac 1{2^{w_\mathrm{H}(i_1,\dotsc,i_n)-1}}.
$$
Therefore,
\begin{align*}
&\frac 12 k! \sigma^2 \beta_{\mathrm{max}}^{k-2} \\
&\qquad + \frac 12 k! \beta_{\mathrm{max}}^{k-2} \sum_{\substack{i_1+\dotsb+i_n=k \\ i_1,\dotsc,i_n \neq 1,k-1,k}} \frac 1{2^{w_\mathrm{H}(i_1,\dotsc,i_n)-1} \binom{w_\mathrm{H}(i_1,\dotsc,i_n)}2} \sum_{\substack{1 \le j_1 < j_2 \le n \\ i_{j_1},i_{j_2} \neq 0}} \sigma_{j_1} \sigma_{j_2} \\
&=\frac 12 k! \sigma^2 \beta_{\mathrm{max}}^{k-2} \\
&\qquad + \frac 12 k! \frac{\sum_{1 \le j_1 < j_2 \le n} \sigma_{j_1} \sigma_{j_2}}{\binom n2} \beta_{\mathrm{max}}^{k-2} \sum_{\substack{i_1+\dotsb+i_n=k \\ i_1,\dotsc,i_n \neq 1,k-1,k}} \frac 1{2^{w_\mathrm{H}(i_1,\dotsc,i_n)-1}}
\end{align*}
For a fixed $2 \le w \le \lfloor \frac k2 \rfloor$, the number of integer solutions of
\begin{gather*}
i_1+\dotsb+i_n=k, \\
i_1,\dotsc,i_n \neq 1,k-1,k, \\
w=w_\mathrm{H}(i_1,\dotsc,i_n),
\end{gather*}
is $\binom nw$ times the number of integer solutions of (everywhere we assume that $\binom xy = 0$ if $x < y$)
\begin{gather*}
x_1 + \dotsb + x_w = k, \\
2 \le x_1, \dotsc, x_w \le k-2,
\end{gather*}
Since $w \ge 2$, the constraint $\le k-2$ is redundant and can be dropped. Next, by substitution $x_i = y_i+1$, the equation is equivalent to the following one:
\begin{gather*}
y_1 + \dotsb + y_w = k-w, \\
y_1, \dotsc, y_w \ge 1,
\end{gather*}
which has $\binom{k-w-1}{w-1}$ integer solutions (the number of $w$-compositions).
Altogether, we continue:
\begin{align*}
&\frac 12 k! \sigma^2 \beta_\max^{k-2} \\
&\qquad + \frac 12 k! \frac{\sum_{1 \le j_1 < j_2 \le n} \sigma_{j_1} \sigma_{j_2}}{\binom n2} \beta_\max^{k-2} \sum_{\substack{i_1+\dotsb+i_n=k \\ i_1,\dotsc,i_n \neq 1,k-1,k}} \frac 1{2^{w_\mathrm{H}(i_1,\dotsc,i_n)-1}} \\
&\le \frac 12 k! \sigma^2 \beta_{\mathrm{max}}^{k-2} + \frac 12 k! \left(\sum_{1 \le j_1 < j_2 \le n} \sigma_{j_1} \sigma_{j_2}\right) A_{nk} \beta_{\mathrm{max}}^{k-2},
\end{align*}
where
$$
A_{nk} \triangleq \frac{1}{\binom n2}\sum_{w=2}^{\lfloor \frac k2 \rfloor} \frac 1{2^{w-1}} \binom nw \binom{k-w-1}{w-1}
$$
Now we will require that
\begin{align*}
&\frac 12 k! \sigma^2 \beta_{\mathrm{max}}^{k-2} + \frac 12 k! \left(\sum_{1 \le j_1 < j_2 \le n} \sigma_{j_1} \sigma_{j_2}\right) A_{nk} \beta_{\mathrm{max}}^{k-2} \\
&\qquad \le \frac 12 k! \sigma^2 \left(\alpha_{nk} \sqrt n \beta_{\mathrm{max}}\right)^{k-2},
\end{align*}
where $\alpha_{nk} > 0$ is some expression we want to bound. The last inequality is equivalent to
\begin{equation}
\sigma^2 \left( \alpha_{nk}^{k-2} n^{\frac k2 - 1} - 1 \right)
- A_{nk} \sum_{1 \le i<j \le n} \sigma_i \sigma_j \ge 0.\tag{3}
\end{equation}
Observe that
$$
\sigma^2 = \frac{\sum_{1 \le i<j \le n} (\sigma_i - \sigma_j)^2 + 2\sigma_i\sigma_j}{n-1}
$$
Thus, (3) is equivalent to
$$
\sum_{1 \le i<j \le n} (\sigma_i - \sigma_j)^2
+ \left(2 - \frac{(n-1)A_{nk}}{\alpha_{nk}^{k-2} n^{\frac k2 - 1} - 1} \right) \sum_{1 \le i<j \le n} \sigma_i \sigma_j \ge 0.
$$
Since all standard deviations are non-negative, it is sufficient to require that
\begin{gather*}
2 - \frac{(n-1)A_{nk}}{\alpha_{nk}^{k-2} n^{\frac k2 - 1} - 1} = 0, \\
\alpha_{nk}^{k-2} = \frac{1}{n^{\frac k2 - 1}}
+ \frac{n-1}{2 n^{\frac k2 - 1}} A_{nk}.
\end{gather*}
Now, we want to find an upper bound for all $n \ge 2$ and $k \ge 4$ on $\alpha_{nk}^{k-2}$:
\begin{align*}
&\frac{1}{n^{\frac k2 - 1}}
+ \frac{n-1}{2 n^{\frac k2 - 1}} A_{nk} \\
&=\frac{1}{n^{\frac k2 - 1}}
+ \frac{n-1}{2 n^{\frac k2 - 1}} \frac{1}{\binom n2}\sum_{w=2}^{\lfloor \frac k2 \rfloor} \frac 1{2^{w-1}} \binom nw \binom{k-w-1}{w-1} \tag{4}\\
&=\frac{1}{n^{\frac k2 - 1}}
+ \frac{1}{n^{\frac k2}} \sum_{w=2}^{\lfloor \frac k2 \rfloor} \frac 1{2^{w-1}} \binom nw \binom{k-w-1}{w-1} \\
&\le\frac{1}{n^{\lfloor\frac k2\rfloor - 1}}
+ \frac{1}{n^{\lfloor\frac k2\rfloor}} \sum_{w=2}^{\lfloor \frac k2 \rfloor} \frac 1{2^{w-1}} \binom nw \binom{2 \lfloor \frac k2 \rfloor-w}{w-1} \\
&= \left[ \ell \triangleq \lfloor \frac k2 \rfloor \ge 2 \right] \\
&= \frac{1}{n^{\ell-1}} + \frac{1}{n^\ell} \sum_{w=2}^{\ell} \frac{1}{2^{w-1}} \binom nw \binom{2\ell-w}{w-1} \\
&=\left[\text{for $w \le \ell$, we have } \frac{\binom nw}{n^\ell} \le \frac{1}{w!n^{\ell-w}}\right] \\
&\le \frac{1}{n^{\ell-1}} + \sum_{w=2}^{\ell} \frac{1}{2^{w-1} w! n^{\ell-w}} \binom{2\ell-w}{w-1}
\end{align*}
The last expression for a fixed $\ell \ge 2$ has the following form:
$$
a_0 + \frac{a_1}{n} + \frac{a_2}{n^2} + \dotsb + \frac{a_{\ell-1}}{n^{\ell-1}},
$$
where $a_0,a_1,\dotsc,a_{\ell-1}$ are positive constants. Thus, for a fixed $\ell \ge 2$, it is decreasing in $n$, and
\begin{align*}
&\frac{1}{n^{\ell-1}} + \sum_{w=2}^{\ell} \frac{1}{2^{w-1} w! n^{\ell-w}} \binom{2\ell-w}{w-1} \\
&\le \left.\left( \frac{1}{n^{\ell-1}} + \sum_{w=2}^{\ell} \frac{1}{2^{w-1} w! n^{\ell-w}} \binom{2\ell-w}{w-1} \right)\right|_{n=2} \\
&= \frac{1}{2^{\ell-1}} + \frac{1}{2^{\ell-1}}\sum_{w=2}^{\ell} \frac{1}{w!} \binom{2\ell-w}{w-1}
\end{align*}
Since $\frac{1}{w!} \binom{2l-w}{w-1} \le 2^{l-w}$, we have:
\begin{align*}
&\frac{1}{2^{\ell-1}} + \frac{1}{2^{\ell-1}}\sum_{w=2}^{\ell} \frac{1}{w!} \binom{2\ell-w}{w-1} \\
&\qquad\le \frac{1}{2^{\ell-1}} + \frac{1}{2^{\ell-1}}\sum_{w=2}^{\ell} 2^{\ell-w} = 1.
\end{align*}
Therefore, $\alpha_{nk}^{k-2} \le 1$, and thus, $\alpha_{nk} \le 1$. We can take $\alpha_{nk} = 1$ and thus obtain that for any $n \ge 2$ and $k \ge 4$,
\begin{align*}
\left| \mathbb E{\left(X_1+\dotsb+X_n\right)^k} \right| \le \frac 12 k! (\sqrt n \beta_{\mathrm{max}})^{k-2},
\end{align*}
and $X_1 + \dotsb + X_n$ satisfies the Bernstein's condition with the parameter $\sqrt n \beta_{\mathrm{max}}$.
If instead of bounding, we calculate (4) with Mathematica for
$n$ and
$k$ up to
$100$, we get that
$\alpha_{nk} \le \frac{\sqrt[4]{7}}2 \approx 0.81$. This achieved for
$n=2$ and
$k=6$ and seems to decrease for larger
$n$ and
$k$. So perhaps the constant in front of
$\sqrt n \beta_\max$ can be further improved.