The equation of state (EOS) of the isospin asymmetric nucleonic matter (IASM) is studied in the Lowest Order Constrained Variational (LOCV) formalism, using the AV18 two-body potential with a density-dependent three-body interaction (TBI), which was proposed by Lagaris and Pandharipande (LP) in 1981. In the LOCV framework, an adjustment mechanism is introduced for finding the appropriate parameters of the TBI. It is shown that, assuming the TBI, the corresponding LOCV symmetric nuclear (pure neutron) matter results fairly agree with those of Fermi Hypernetted Chain (FHNC). Also, the contribution of the averaged three-body cluster (TBC) energy is evaluated, at different proton to neutron ratios R. Employing the TBI, up to ρ=0.5 fm−3, the TBC energy ranges from −6 to 6 MeV per particle which considerably changes the IASMEOS, especially at low densities. Moreover, employing the TBI, the LOCVIASM state-averaged correlation functions (effective potentials) and the nucleon–nucleon distribution functions (NNDF) are reported. It is illustrated that, by increasing the density (ratio R), the impact of the TBI becomes much more evident. Finally, it should be noted that, the IASMFHNC and Monte Carlo (MC) data are not available for comparison.