Transition Flux Formula for the Electronic Coupling Matrix Element

Respiratory complex II acting as a homeostatic regulatory sensor

The succinate-ubiquinone oxidoreductase (SQR) complex connects two of the cell's most vital energy-producing metabolic processes: the tricarboxylic acid cycle and the electron transport chain. Hence, the SQR complex is essential in cell metabolism, and its malfunction leads to the progression of multiple metabolic disorders and other diseases, such as cancer. In the current study, we calculated the electron tunneling (ET) pathways between the different redox systems in the SQR complex, including the SQR ligands and the distant heme b redox center, using the broken-symmetry semi-empirical ZINDO method. Interestingly, we discovered a water channel running from the mitochondrial matrix, filling the space between Fe3S4 and heme b redox centers. To investigate the physiological function of the water channel, we performed extensive molecular dynamics (MD) simulations of the membrane-embedded SQR complex in small and large water boxes, representing regular (MDA) and extended (MDB) volume states, respectively. We found that under regular volume conditions (MDA), the ET reaction is conducted through both the iron-sulfur cluster chain (i.e., pathway A) and through heme b (i.e., pathway B). Hence, the SQR complex encompasses an internal interferometer similar to the Mach-Zender interferometer, such that the tunneling electron experiences a self-interference effect through pathways A and B, enhancing the SQR complex's overall ET thermodynamics and favoring the forward ET direction of oxidizing succinate to fumarate and reducing ubiquinone to ubiquinol. On the other hand, we found that under extended volume conditions (MDB), the internal water channel of the SQR complex "senses" the expansion in the mitochondrial volume, pushing the heme b and Fe4S3 redox centers apart and hence lowering the SQR equilibrium constant to almost unity. Therefore, the SQR complex could be driven to work in the reverse direction, catalyzing the production of ubiquinone molecules essential for the physiological function of respiratory complexes I and III and restoring the inner-mitochondrial membrane potential, which leads to restoring the function of the H-K anti-porter, pumping K+ outward from the matrix and restoring the regular mitochondrial volume.

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.

Fragment-orbital tunneling currents and electronic couplings for analysis of molecular charge-transfer systems

In theoretical charge-transfer research, calculation of the electronic coupling element is crucial for examining the degree of the electronic donor-acceptor interaction. The tunneling current (TC), representing the magnitudes and directions of electron flow, provides a way of evaluating electronic couplings, along with the ability of visualizing how electrons flow in systems. Here, we applied the TC theory to π-conjugated organic dimer systems, in the form of our fragment-orbital tunneling current (FOTC) method, which uses the frontier molecular-orbitals of system fragments as diabatic states. For a comprehensive test of FOTC, we assessed how reasonable the computed electronic couplings and the corresponding TC densities are for the hole- and electron-transfer databases HAB11 and HAB7. FOTC gave 12.5% mean relative unsigned error with regard to the high-level ab initio reference. The shown performance is comparable with that of fragment-orbital density functional theory, which gave the same error by 20.6% or 13.9% depending on the formulation. In the test of a set of nucleobase π stacks, we showed that the original TC expression is also applicable to nondegenerate cases under the condition that the overlap between the charge distributions of diabatic states is small enough to offset the energy difference. Lastly, we carried out visual analysis on the FOTC densities of thiophene dimers with different intermolecular alignments. The result depicts an intimate topological connection between the system geometry and electron flow. Our work provides quantitative and qualitative grounds for FOTC, showing it to be a versatile tool in characterization of molecular charge-transfer systems.

Electron tunneling in proteins program

We developed a unique integrated software package (called Electron Tunneling in Proteins Program or ETP) which provides an environment with different capabilities such as tunneling current calculation, semi-empirical quantum mechanical calculation, and molecular modeling simulation for calculation and analysis of electron transfer reactions in proteins. ETP program is developed as a cross-platform client-server program in which all the different calculations are conducted at the server side while only the client terminal displays the resulting calculation outputs in the different supported representations. ETP program is integrated with a set of well-known computational software packages including Gaussian, BALLVIEW, Dowser, pKip, and APBS. In addition, ETP program supports various visualization methods for the tunneling calculation results that assist in a more comprehensive understanding of the tunneling process. © 2016 Wiley Periodicals, Inc.

Novel Inhibitors for a Novel Binding Site in Respiratory Complex III

A new binding site and potential novel inhibitors of the respiratory complex III are described. The site is located at the opposite side of the enzyme with respect to ubiquinol binding site (Qo site), and distinctly different from both Qo and Qi sites (hence designated as Non-Q binding site, NQ). NQ site binding pocket extends up close to Phe90 residue, an internal switch (LH switch) that regulates electron transfer between heme bL and heme bH of the low potential redox chain. Docking studies and molecular dynamics simulations of different molecules to the NQ site revealed potential ligands which exhibit a novel inhibitory effect for bc1 complex by switching the LH switch to "off" conformation, thereby shutting down electron transfer in the low potential redox chain. Moreover, the novel inhibitors have lower binding affinity for both Qo and Qi sites, and hence do not interfere with binding of the natural ligands to those sites. The inhibitory activity of those novel ligands in bc1 complex is suggested to promote the production of reactive oxygen species (ROS) at the Qo site. Hence those ligands are potential candidates for designing new "mitocan" drugs.