Deviations from statistical binding, that is cooperativity, in self-assembled polynuclear complexes partly result from intermetallic interactions DeltaE(M,M), whose magnitudes in solution depend on a balance between electrostatic repulsion and solvation energies. These two factors have been reconciled in a simple point-charge model, which suggests severe and counter-intuitive deviations from predictions based solely on the Coulomb law when considering the variation of DeltaE(M,M) with metallic charge and intermetallic separation in linear polynuclear helicates. To demonstrate this intriguing behaviour, the ten microscopic interactions that define the thermodynamic formation constants of some twenty-nine homometallic and heterometallic polynuclear triple-stranded helicates obtained from the coordination of the segmental ligands L1-L11 with Zn(2+) (a spherical d-block cation) and Lu(3+) (a spherical 4f-block cation), have been extracted by using the site binding model. As predicted, but in contrast with the simplistic coulombic approach, the apparent intramolecular intermetallic interactions in solution are found to be i) more repulsive at long distance (${{ m{Delta }}E_{{ m{1 - 4}}};{{ m{Lu}}{ m{,Lu}}} }$>${{ m{Delta }}E_{{ m{1 - 2}}};{{ m{Lu}}{ m{,Lu}}} }$), ii) of larger magnitude when Zn(2+) replaces Lu(3+) (${{ m{Delta }}E_{{ m{1 - 2}}};{{ m{Zn}}{ m{,Lu}}} }$>${{ m{Delta }}E_{{ m{1 - 2}}};{{ m{Lu}}{ m{,Lu}}} }$) and iii) attractive between two triply charged cations held at some specific distance (${{ m{Delta }}E_{{ m{1 - 3}}};{{ m{Lu}}{ m{,Lu}}} }$<0). The consequences of these trends are discussed for the design of polynuclear complexes in solution.