Optimal settings for large DFT calculations

We are looking into large molecule containing Co atom with ligands (49 atoms). With def2-TZVP basis, it comes to 1057 basis functions. We are interested in using wB97X-D and wB97M-V functionals.
What performance-enhancing options can we use for such calculation?