A new feature in our MiMiC-based QM/MM implementation compared to the previous implementation in the CPMD program by Laio et al. is the ability to use an open-ended multipole expansion for use in the long-range coupling. Using higher orders in the multipole expansion will in many cases produce a more accurate electrostatic potential thus allowing the long-range coupling to be used at shorter distances. Compared to the previous implementation, this reduces the computational cost of the QM-MM interactions even further without compromising the accuracy. To investigate how much the short-range cutoff distances can be decreased while still retaining a high accuracy, we performed single-point calculations on a series of snapshots using varying orders of the multipole expansion and short-range cutoff distances, and compared key properties to calculations where the short-range coupling was used for the entire system. Specifically, we examined the total energy, forces, and, as a measure of how well the electron density is reproduced, the dipole moment of the QM subsystem. For the GB1 system, we investigate the Hirshfeld partial charge of the zinc ion instead of the dipole moment.