Proteins catalyze crucial reactions via unstable, high-energy chemical intermediates. In the absence of physiological substrates, activated redox cofactors become ticking time bombs, capable of producing oxidative damage to the protein. In PNAS, Gray and Winkler (1) propose that chains of tryptophan (Trp) and tyrosine (Tyr) residues may serve as escape routes for potentially damaging, highly oxidizing electron holes (oxidizing equivalents) that are generated in enzymatic reactions. The aromatic residues are suggested to provide charge-shuttling pathways that lead out to the protein surface from buried active sites. Once outside the protein, the holes can be safely scavenged by cellular reductants.Since at least the 1960s, aromatic amino acids were postulated to play a special role in biological redox chemistry as electron hole carriers (2). Circadian rhythm photochemistry, photosynthetic water splitting, nucleic acid biosynthesis, and cell signaling are now well understood to engage redox-active Trp and Tyr residues as essential parts of their function (3). Could Tyr and Trp also play a protective role during enzymatic reactions (such as the oxidation of aliphatic C-H bonds) that require high-energy intermediates (4, 5)? The answer will be dictated by the set of time scales at play: Hole hopping from the activated redox cofactor to the first Trp/Tyr must be slower than catalysis (to avoid disabling the enzyme's valuable catalytic function) but faster than oxidative protein damage. Gray and Winkler's (1) query of the Protein Data Bank (e.g., they searched for Trp/Tyr pairs within a 5-Å edge-to-edge cutoff distance) provides a critical starting point for more detailed investigations of these hole-transfer chains.Electron-transfer (ET) rates in biochemical reactions are usually proportional to the square of the electronic coupling interaction energy (V, Fig. 1B) between the electron donor and acceptor, facilitated by the protein medium, or by direct van der Waals (vdW) contact between redox species. V decays rapidly with distance, and the steepness of this decay depends on the structure of the intervening medium (9). ET pathways and relative donor/acceptor orientations influence wave function overlap, and thus determine V and the ET rate. The relative Trp and Tyr orientations fine-tune the hole-hopping rates. Two orientations are common among protein aromatics (Fig. 1C), the so-called parallel-displaced (PD) and the T-shape (T) orientations, where ring interactions capture quadrupolar electrostatic stabilization (10). PD geometries may use strong stacking interactions (akin to sigma bonding between p-orbitals) to maximize V; however, the strength of such interactions depends on orbital symmetry, as shown in Fig. 1D for indole pairs. Indeed, the greater than 10-fold stronger V for Tyr90-Trp231 in Fig. 2 arises from its PD geometry, compared with the T geometry of the other pairs.Within a subset of orientations (e.g., T or PD), additional geometric effects between Fig. 1. The electronic coupling interaction energy (...