Abstract
Bifurcating electron transferring flavoproteins (Bf-ETFs) tune chemically identical flavins to two contrasting roles. To understand how, we used hybrid quantum mechanical molecular mechanical calculations to characterize noncovalent interactions applied to each flavin by the protein. Our computations replicated the differences between the reactivities of the flavins: the electron transferring flavin ((ET)flavin) was calculated to stabilize anionic semiquinone (ASQ) as needed to execute its single-electron transfers, whereas the Bf flavin ((Bf)flavin) was found to disfavor the ASQ state more than does free flavin and to be less susceptible to reduction. The stability of (ET)flavin ASQ was attributed in part to H-bond donation to the flavin O2 from a nearby His side chain, via comparison of models employing different tautomers of His. This H-bond between O2 and the ET site was uniquely strong in the ASQ state, whereas reduction of (ET)flavin to the anionic hydroquinone (AHQ) was associated with side chain reorientation, backbone displacement, and reorganization of its H-bond network including a Tyr from the other domain and subunit of the ETF. The Bf site was less responsive overall, but formation of the (Bf)flavin AHQ allowed a nearby Arg side chain to adopt an alternative rotamer that can H-bond to the (Bf)flavin O4. This would stabilize the anionic (Bf)flavin and rationalize effects of mutation at this position. Thus, our computations provide insights on states and conformations that have not been possible to characterize experimentally, offering explanations for observed residue conservation and raising possibilities that can now be tested.