Wavefunction-based methods are probably not applicable to such a big system, so DFT would be the first option. Hopefully, some insights might be obtained by using QM-MM or DFT-MM techniques, too.