Relationship between Hartree-Fock and Hückel coefficients
Let start by definition the Hartree-Fock-Roothaan equation:

$$\sum_\nu F_{\mu\nu}C_{\nu i}=\varepsilon_i\sum_\nu S_{\mu\nu}C_{\nu i}$$

Where the Fock matrix elements are giving as:

$$F_{\mu \nu} = \left&lt; \phi_\mu \left| \hat{f} \right|\phi_\nu\right&gt;$$

And the overlap matrix is given as:

$$S_{\mu \nu} = \left&lt;  \left. \phi_\mu \right|\phi_\nu\right&gt;$$

Here $\phi_\nu$ is our atomic orbitals (i.e. what you get from your basis set).  

When we do a Hartree-Fock calculation we want to find a set of orthonormal orbitals, also known as molecular orbitals. The molecular orbitals are constructed as:

$$\psi_i=\sum_{\mu=1}^KC_{\nu i}\phi_\nu$$

For the molecular orbitals to be orthonormal we can write the following condition:

$$\left&lt; \left. \psi_i \right|\psi_j\right&gt;=\delta_{ij}$$

This means that the orbitals have zero overlap with one-another and an overlap of $1$ with them self. If we insert the expansion of the molecular orbitals in atomic orbitals in the above equation we get:

$$\sum_{\mu=1}^K\sum_{\nu=1}^LC_{\nu i}^\dagger C_{\mu j}\left&lt; \left.\phi_\nu \right|\phi_\mu\right&gt;=\delta_{ij}$$

I.e. we can see that when we do a Hartree-Fock calculation we have the constraint on our molecular orbital coefficients that have to transform our atomic orbital overlap to a matrix with ones on the diagonal and zero elsewhere. This can also be written more compact in matrix form:

$$\mathbb{C^\dagger SC=1}$$

Which is a condition equivalent to the one you know from Hückel theory. At last let me make a breif comment about basis set, since you also talked about those in your question. In general we use contracted Gaussian type basis sets, just like STO3G, that you mentioned. Our atomic orbitals can thus be written as:

$$\phi_\nu = \sum_{k=1}^M N_k d_k P(r)\cdot e^{-\alpha_k\cdot r^2}$$

Here $P(r)$ is a function related to the angular moment quantum number.The way the contraction coefficients $d_k$ are choose implies that we normalize our Gaussian functions, i.e:

$$N_k^2 \int_{-\infty}^\infty \left( P(r)\cdot e^{-\alpha_k\cdot r^2} \right)^2 dr = 1$$

We can now write $\left&lt; \left.\phi_\nu \right|\phi_\nu\right&gt; = 1$ as:

$$\overline{N}\sum_{k=1}^M N_k d_k\sum_{l=1}^T N_l d_l\left&lt; \left.P(r)\cdot e^{-\alpha_k\cdot r^2} \right| P(r)\cdot e^{-\alpha_l\cdot r^2}\right&gt; = 1$$

Where $\overline{N}$ is the normalization we need. As we can see there is no easy condition we can put on the contraction coefficients, in order to avoid evaluation $\overline{N}$. It can be shown (I will not show it here) that $\overline{N}$ only depends on the total angular moment and on the Gaussian exponents, thus it is possible to choose the contraction coefficients in such a way, that we do not need to evaluate $\overline{N}$.