In this report, we study the photoactivation process in Escherichia coli DNA photolyase, involving long-range electron transport along a conserved chain of Trp residues between the protein surface and the flavin adenine dinucleotide (FAD) cofactor. Fully coupled nonadiabatic (Ehrenfest) quantum mechanics/molecular mechanics (QM/MM) simulations allow us to follow the time evolution of charge distributions over the natural time scale of multiple charge transfer events and conduct rigorous statistical analysis. Charge transfer rates in excellent agreement with experimental data are obtained without the need for any system-specific parametrization. The simulations are shown to provide a more detailed picture of electron transfer than a classical analysis of Marcus parameters. The protein and solvent both strongly influence the localization and transport properties of a positive charge, but the directionality of the process is mainly caused by solvent polarization. The time scales of charge movement, delocalization, protein relaxation and solvent reorganization overlap and lead to nonequilibrium reaction conditions. All these contributions are explicitly considered and fully resolved in the model used and provide an intricate picture of multistep biochemical electron transfer in a flexible, heterogeneous environment.