In the conventional GW method, the three-point vertex function (Γ) is approximated to unity (Γ ∼ 1). Here, we developed an all-electron first-principles GWΓ method beyond a conventional GW method by considering a first-order three-point vertex function (Γ(1) = 1 + iGGW) in a one-electron self-energy operator. We applied the GWΓ method to simulate the binding energies (BEs) of B1s, C1s, N1s, O1s, and F1s for 19 small-sized molecules. Contrary to the one-shot GW method [or G0W0(LDA)], which underestimates the experimentally determined absolute BEs by about 3.7eV for B1s, 5.1eV for C1s, 6.9eV for N1s, 7.8eV for O1s, and 5.8eV for F1s, the GWΓ method successfully reduces these errors by approximately 1-2eV for all the elements studied here. Notably, the first-order three-point vertex corrections are more significant for heavier elements, following the order of F > O > N > C > B1s. Finally, the computational cost analysis revealed that one term in the GWΓ one-electron self-energy operator, despite being computationally intensive, contributes negligibly (<0.1eV) to the C1s, N1s, O1s, and F1s.