Mathematical and numerical analysis of the time-dependent Maxwell--Schr\"{o}dinger Equations in the Coulomb gauge