Many isotopes of Np, Pu, Am and Cm around N = 126 shell still have not been produced in the laboratory. This work with the aims to investigate the cross sections and yields of the neutron-deficient nuclei Np, Pu, Am and Cm produced in the proton-induced spallations of transuranium element. The isospin-dependent quantum molecular dynamics (IQMD) model is applied to study the dynamical process of reaction, and the subsequent decay process is simulated by the GEMINI++ model. The IQMD-GEMINI++ model is applied to calculate the cross section, kinetic energy and angle distribution of the isotopic productions around N=126. The Lindhand Scharff and Schiott theory is applied to calculate the energy loss of different heavy nuclei in target material. Comparison between the data and calculations shows that the IQMD-GEMINI++ model can reproduce the production cross section of the neutron-deficient nuclei in spallation within 1 order of magnitude approximately. The maximum cross section of undiscovered isotope of Np, Pu, Am and Cm is about $10^{-5} $ mb, while the kinetic energy of productions are all less than 16 MeV. The angle distribution shows that the emission direction of production is in the forward downward angle mostly. The range of the production in the target is within the range of $10^{-7}$ to $10^{-5} $ cm. This range is effective target thickness for online identification of the undiscovered isotope. Based on the effective thickness of the target and assuming a intensity 120 $\mu$A of the proton beam, the yields of the undiscovered neutron-deficient nuclei are calculated. Productions of the undiscovered isotope of Np, Pu, Am and Cm by the proton-induced spallation of transuranium element are feasible. However, the experimental techniques for on-line identification of the neutron-deficient nuclei produced in the proton-induced spallation should be developed.