The decoupled sites representation (DSR) for one type of ligand allows to regard complex overall titration curves as sum of classical Henderson-Hasselbalch (HH) titration curves. In this work we transfer this theoretical approach to molecules with different types of interacting ligands (e.g. protons and electrons), prove the existence of decoupled systems for n 1 and one binding sites for two different ligands, and point out some difficulties and limits of this transfer. A major difference to the DSR for one type of ligand is the loss of uniqueness of the decoupled system. However, all decoupled systems share a unique set of microstate probabilities and each decoupled system corresponds to a certain permutation of these microstate probabilities. Moreover, we show that the titration curve of a certain binding site in the original system can be regarded as linear combination of the titration curves of the individual sites of the decoupled system if the weights of the linear combination are substituted by functions in the activity of the second ligand. In the underlying model with only pairwise interaction, an important observation of our theoretical investigation is the following: Even though the binding sites of ligand L 1 may not interact directly, they can show secondary interaction due to the interaction with the second type of ligand. This means, if the activity of the second ligand is fixed and we regard the 1-dimensional titration curve of an individual binding site for ligand L 1 depending on its activity, we may observe a strong deviation from the classical HH shape in spite of non-interacting sites for ligand L 1 .