We have studied electron transfers (ET) between electron donors and acceptors, taking as illustrative example the case of ET in methionine enkephalin. Recent pulse and gamma radiolysis experiments suggested that an ultrafast ET takes place from the C-terminal tyrosine residue to the N-terminal, oxidized, methionine residue. According to standard theoretical frameworks like the Marcus theory, ET can be decomposed into two successive steps: i) the achievement through thermal fluctuations, of a set of nuclear coordinates associated with degeneracy of the two electronic states, ii) the electron tunneling from the donor molecular orbital to the acceptor molecular orbital. Here, we focus on the analysis of the time-dependent electronic dynamics during the tunneling event. This is done by extending the approaches based on the topological analyses of stationary electronic density and of the electron localization function (ELF) to the time-dependent domain. Furthermore, we analyzed isosurfaces of the divergence of the current density, showing the paths that are followed by the tunneling electron from the donor to the acceptor. We show how these functions can be calculated with constrained density functional theory. Beyond this work, the topological tools used here can open up new opportunities for the electronic description in the time-dependent domain.