An overview of the interaction between charged macromolecules in aqueous solution is presented. The starting point is the dielectric continuum model and the Debye-Hückel equation. The usefulness of the simple theory is emphasized in particular for biological macromolecules, whose net charge or surface charge density often is low. With more highly charged macromolecules or aggregates it may be necessary to go beyond the simple Debye-Hückel theory and invoke the non-linear Poisson-Boltzmann equation or even to approach an exact solution using Monte Carlo simulations or similar techniques. The latter approach becomes indispensable when studying systems with divalent or multivalent (counter)-ions. The long range character of the electrostatic interactions means that charged systems of varying geometry -spheres, planes, cylinders... -often have many properties in common. Another consequence is that the detailed charge distribution on a macromolecule is less important. Many biological macromolecules contain titratable groups, which means that the net charge will vary as a consequence of solution conditions. This gives an extra attractive contribution to the interaction between two macromolecules, which might be particularly important close to their respective isoelectric points. The treatment of flexible polyelectrolytes/polyampholytes requires some extra efforts in order to handle the increasingly complex geometry. A theoretical consequence is that the number of parameters -chain length, charge density, polydispersity etc -prohibits the presentation of a simple unified picture. An additional experimental, and theoretical, difficulty in this context is the slow approach towards equilibrium, in particular with high molecular weight polymers. A few generic situations where polyelectrolytes can act both as stabilizers and coagulants can, however, be demonstrated using simulation techniques.