This work presents an extension of the projector operator embedding scheme of Manby and coworkers J. Chem. Theory Comput., 8, 2564 (2012) in a multicomponent framework. Here, a molecular system containing electrons and other types of quantum species is divided into a wavefunction (WF) subsystem of interest and a density functional theory (DFT) environment. The WF-in-DFT partition decreases the computational costs by partially truncating the WF subsystem basis set at the cost of introducing a controllable embedding error. To explore the applicability of the MC extension, third-order propagator-in-DFT calculations were performed for positron-anion complexes for alkoxides and carboxylates with carbon chains of different sizes. For these systems, it was found that selecting a WF subsystem with the positron and only the oxygen atoms caused an error of 0.1 eV or lower in the positron-binding energies, while reducing between 33% to 55% the basis set size. The reduction of computational costs achieved with the embedding scheme allowed us to improve the molecular positron-binding energy predictions by performing complete basis set limit extrapolations. Combining the WF-in-DFT embedding and the complete basis set extrapolation, positronium aliphatic alkoxides were predicted to be energetically stable by 0.3 eV with respect to Ps emission. Similarly, positronium carboxylates, both aromatic and aliphatic, were predicted to be stable by 1.3 eV.