In this paper, nonlinear boundary value problems for fourth-order fractional integro-differential equations are solved by modified Laplace transform variational iteration method (MLT-VIM). This modification is designed to improve the accuracy of the approximate solutions obtained by LT-VIM. The suitable initial approximation will be selected based on a rule or standard to adjust the determination of the initial guess that satisfies the boundary conditions. Furthermore, the residual error will be cancelled at several points of the appropriate interval. The convergence conditions of the proposed method are investigated. Some illustrative examples of nonlinear boundary value problems for fourth-order fractional integro-differential equations are presented in order to demonstrate the effectiveness of the proposed method. Comparisons are made between the results obtained by MLT-VIM and LT-VIM depending on the exact solutions, revealing that the MLT-HPM contributes to reduce the amount of computational work to obtain the first-order approximate solution which greatly accelerates the convergence of the solution. Whereas LT-HPM needs more iterations to obtain a suitable approximate solution, which results in an increase in the computational workload.