Advanced PE methods are significantly more computationally intensive than non-polarizable QM/MM methods due to a quadratic scaling with the number of polarizable sites. Additionally, the necessity to frequently recalculate electric fields emanating from the quantum region, and to update the Fock matrix to reflect the effects of induced dipoles, further escalates the computational requirements. In this work, a number of approximations have been followed to reduce the computational cost of PE for this large system. These approximations cover MM fields, and approximations on the QM/MM coupling.