Valence topological Charge-Transfer (CT) indices are applied to the calculation of pH at the isoelectric point (pI). The model is generalized for molecules with heteroatoms. The ability of the indices for the description of molecular charge distribution is established by comparing them with the pI of 21 amino acids. Linear correlation models are obtained. The CT indices improve multivariable regression equations for pI. The variance decreases by 95%. No superimposition of the corresponding G(k)-J(k) and G(k)(V)-J(k)(V) pairs is observed in most fits, which diminishes the risk of collinearity. The inclusion of heteroatoms in pi-electron system is beneficial for the description of pI, the because of either the role of the additional p orbitals provided by heteroatom or the role of steric factors in pi-electron conjugation. The use of CT/valence CT indices {G(k), J(k), G(k)(V), J(k)(V)} gives limited results for modelling amino acid pI. The inclusion of the number of acidic/basic groups improves the models. The effect is specially noticeable for amino acids with more than two functional groups. The fitting line for 21 amino acids is used to estimate lysozyme pI replacing (1 + Delta n/n(T)) with (M + Delta n)/n(T). The results for lysozyme fragments can estimate the pI of the whole protein with 1-13% error.