The hybrid quantum mechanics/molecular mechanics (QM/MM) approach, which combines the accuracy of QM methods with the efficiency of MM methods, is widely used in the study of complex systems. However, past QM/MM implementations often neglect or face challenges in addressing nuclear quantum effects, despite their crucial role in many key chemical and biological processes. Recently, our group developed the constrained nuclear-electronic orbital (CNEO) theory, a cost-efficient approach that accurately addresses nuclear quantum effects, especially quantum nuclear delocalization effects. In this work, we integrate CNEO with the QM/MM approach through the electrostatic embedding scheme and apply the resulting CNEO QM/MM to two hydrogen-bonded complexes. We find that both solvation effects and nuclear quantum effects significantly impact hydrogen bond structures and dynamics. Notably, in the glutamic acid–glutamate complex, which mimics a common low barrier hydrogen bond in biological systems, CNEO QM/MM accurately predicts nearly equal proton sharing between the two residues. With an accurate description of both quantum nuclear delocalization effects and environmental effects, CNEO QM/MM is a promising new approach for simulating complex chemical and biological systems.