In this work, we develop such an integration using the QM/MM electrostatic embedding scheme. By studying two bimolecular complex systems, one of which is of strong biological relavance, we show that CNEO QM/MM outperforms conventional QM/MM in describing hydrogen bonds and hydrogen-bond dynamics in the condensed phase, aligning well with experimental evidence.