Internal switches modulating electron tunneling currents in respiratory complex III

Respiratory Complex III: A Bioengine with a Ligand-Triggered Electron-Tunneling Gating Mechanism

Respiratory complex III (a.k.a., the bc1 complex) plays a key role in the electron transport chain in aerobic cells. The bc1 complex exhibits multiple unique electron tunneling (ET) processes, such as ET-bifurcation at the Qo site and movement of the Rieske domain. Moreover, we previously discovered that electron tunneling in the low potential arm of the bc1 complex is regulated by a key phenylalanine residue (Phe90). The main goal of the current work is to study the dynamics of the key Phe90 residue in the electron tunneling reaction between heme bL and heme bH as a function of the occupancy of the Qo and Qi binding sites in the bc1 complex. We simulated the molecular dynamics of four model systems of respiratory complex III with different ligands bound at the Qo and Qi binding sites. In addition, we calculated the electron tunneling rate constants between heme bL and heme bH along the simulated molecular dynamics trajectories. The binding of aromatic ligands at the Qo site induces a conformational cascade that properly positions the Phe90 residue, reducing the through-space ET distance from ∼7 to ∼5.5 Å and thus enhancing the electron transfer rate between the heme bL and the heme bH redox pair. Also, the binding of aromatic ligands at the Qi site induces conformational changes that stabilize the Phe90 conformational variation from ∼1.5 to ∼0.5 Å. Hence, our molecular dynamics simulation results show an on-demand two-step conformational connection between the occupancy of the Qo and Qi binding sites and the conformational dynamics of the Phe90 residue. Additionally, our dynamic electron tunneling results confirm our previously reported findings that the Phe90 residue acts as an electron-tunneling gate or switch, controlling the electron transfer rate between the heme bL and heme bH redox systems.

Tunneling and Nonadiabatic Effects on a Proton-Coupled Electron Transfer Model for the Qo Site in Cytochrome bc1

Cytochrome bc₁ is a fundamental enzyme for cellular respiration and photosynthesis. This dimeric protein complex catalyzes a proton-coupled electron transfer (PCET) from the reduced coenzyme-Q substrate (Q) to a bimetallic iron-sulfur cluster in the Q_o active site. Herein, we combine molecular dynamics simulations of the complete cytochrome bc₁ protein with electronic-structure calculations of truncated models and a semiclassical tunneling theory to investigate the electron-proton adiabaticity of the initial reaction catalyzed in the Q_o site. After sampling possible orientations between the Q substrate and a histidine side chain that functions as hydrogen acceptor, we find that a truncated model composed by ubiquinol-methyl and imidazole-iron(III)-sulfide captures the expected changes in oxidation and spin states of the electron donor and acceptor. Diabatic electronic surfaces obtained for this model with multiconfigurational wave function calculations demonstrate that this reaction is electronic nonadiabatic, and proton tunneling is faster than mixing of electronic configurations. These results indicate the formalism that should be used to calculate vibronic couplings and kinetic parameters for the initial reaction in the Q_o site of cytochrome bc₁. This framework for molecular simulation may also be applied to investigate other PCET reactions in the Q-cycle or in various metalloproteins that catalyze proton translocation coupled to redox processes.