Electron tunneling routes for the electron transfer from the bacteriopheophytin anion to the primary quinone in the bacterial photosynthetic reaction center of Rhodobactor sphaeroides are investigated by a combined method of molecular dynamics simulations for the protein conformation fluctuation and quantum chemical calculations for the electronic states of the donor, acceptor, and protein medium. The analysis of the tunneling route is made by mapping interatomic electron tunneling currents for each protein conformation. We found that there are two dominant routes mainly passing through TrpM252 (Trp route) or mainly passing through MetM218 (Met route). Actual electron tunneling pathways alternate between the two routes, depending on the protein conformation which varies with time. When either the Trp route or the Met route dominates, the electron tunneling matrix element |TDA| becomes large. When both the Trp route and the Met route dominate, |TDA| becomes very small due to the destructive interference of the electron tunneling currents between the two routes. We found that a linear relationship exists between the value of |TDA| and the inverse of the degree of destructive interference Q for a wide range of values (ca. 3−103 for Q). A similar relationship was also found previously for electron transfer in ruthenium-modified azurins (Kawatsu, T.; Kakitani, T.; Yamato, T. J. Phys. Chem. B 2002, 106, 11356), suggesting that this relationship holds true in general. From these results, we are led to the conclusion that |TDA| cannot exceed a maximum value at Q = 1, even if much variation of |TDA| happens due to the fluctuation of protein conformation. We also conclude that the property of the electron transfer alternates between constructive and destructive interference, due to the fluctuation of protein conformation. It is impossible to keep a system in either constructive or destructive interference because thermal fluctuation of protein conformation takes place.