On the Orthogonality of Orbitals in Subsystem Kohn–Sham Density Functional Theory