We study the multiple Eisenstein series introduced by Gangl, Kaneko and Zagier. We give a proof of (restricted) finite double shuffle relations for multiple Eisenstein series by developing an explicit connection between the Fourier expansion of multiple Eisenstein series and the Goncharov coproduct on Hopf algebras of iterated integrals.