Geometric and electronic structure changes in the copper (Cu) centers in bilirubin oxidase (BOD) upon a four-electron reduction were investigated by quantum mechanics/molecular mechanics (QM/MM) calculations. For the QM region, the unrestricted density functional theory (UDFT) method was adopted for the open-shell system. We found new candidates of the native intermediate (NI, intermediate II) and the resting oxidized (RO) states, i.e., NIH+ and RO₀. Elongations of the Cu-Cu atomic distances for the trinuclear Cu center (TNC) and very small structural changes around the type I Cu (T1Cu) were calculated as the results of a four-electron reduction. The QM/MM optimized structures are in good agreement with recent high-resolution X-ray structures. As the structural change in the TNC upon reduction was revealed to be the change in the size of the triangle spanned by the three Cu atoms of TNC, we introduced a new index (l) to characterize the specific structural change. Not only the wild-type, but also the M467Q, which mutates the amino acid residue coordinating T1Cu, were precisely analyzed in terms of their molecular orbital levels, and the optimized redox potential of T1Cu was theoretically reconfirmed.
Keywords: Cu oxidation state; SONO; X-ray reduction; catalytic intermediate; molecular orbital analysis; protonation.