Understanding Rubredoxin Redox Potentials: Role of H-Bonds on Model Complexes

Long-range charge transfer mechanism of the III2IV2 mycobacterial supercomplex

Aerobic life is powered by membrane-bound redox enzymes that shuttle electrons to oxygen and transfer protons across a biological membrane. Structural studies suggest that these energy-transducing enzymes operate as higher-order supercomplexes, but their functional role remains poorly understood and highly debated. Here we resolve the functional dynamics of the 0.7 MDa III2IV2 obligate supercomplex from Mycobacterium smegmatis, a close relative of M. tuberculosis, the causative agent of tuberculosis. By combining computational, biochemical, and high-resolution (2.3 Å) cryo-electron microscopy experiments, we show how the mycobacterial supercomplex catalyses long-range charge transport from its menaquinol oxidation site to the binuclear active site for oxygen reduction. Our data reveal proton and electron pathways responsible for the charge transfer reactions, mechanistic principles of the quinone catalysis, and how unique molecular adaptations, water molecules, and lipid interactions enable the proton-coupled electron transfer (PCET) reactions. Our combined findings provide a mechanistic blueprint of mycobacterial supercomplexes and a basis for developing drugs against pathogenic bacteria.

Theoretical Study on Redox Potential Control of Iron-Sulfur Cluster by Hydrogen Bonds: A Possibility of Redox Potential Programming

The effect of hydrogen bonds around the active site of Anabaena [2Fe-2S] ferredoxin (Fd) on a vertical ionization potential of the reduced state (IP(red)) is examined based on the density functional theory (DFT) calculations. The results indicate that a single hydrogen bond increases the relative stability of the reduced state, and shifts IP(red) to a reductive side by 0.31–0.33 eV, regardless of the attached sulfur atoms. In addition, the IP(red) value can be changed by the number of hydrogen bonds around the active site. The results also suggest that the redox potential of [2Fe-2S] Fd is controlled by the number of hydrogen bonds because IP(red) is considered to be a major factor in the redox potential. Furthermore, there is a possibility that the redox potentials of artificial iron-sulfur clusters can be finely controlled by the number of the hydrogen bonds attached to the sulfur atoms of the cluster.

Terminal Electron-Proton Transfer Dynamics in the Quinone Reduction of Respiratory Complex I

Complex I functions as a redox-driven proton pump in aerobic respiratory chains. By reducing quinone (Q), complex I employs the free energy released in the process to thermodynamically drive proton pumping across its membrane domain. The initial Q reduction step plays a central role in activating the proton pumping machinery. In order to probe the energetics, dynamics, and molecular mechanism for the proton-coupled electron transfer process linked to the Q reduction, we employ here multiscale quantum and classical molecular simulations. We identify that both ubiquinone (UQ) and menaquinone (MQ) can form stacking and hydrogen-bonded interactions with the conserved Q-binding-site residue His-38 and that conformational changes between these binding modes modulate the Q redox potentials and the rate of electron transfer (eT) from the terminal N2 iron-sulfur center. We further observe that, while the transient formation of semiquinone is not proton-coupled, the second eT process couples to a semiconcerted proton uptake from conserved tyrosine (Tyr-87) and histidine (His-38) residues within the active site. Our calculations indicate that both UQ and MQ have low redox potentials around -260 and -230 mV, respectively, in the Q-binding site, respectively, suggesting that release of the Q toward the membrane is coupled to an energy transduction step that could thermodynamically drive proton pumping in complex I.

Prediction of Reduction Potentials of Copper Proteins with Continuum Electrostatics and Density Functional Theory

Blue copper proteins, such as azurin, show dramatic changes in Cu²⁺/Cu⁺ reduction potential upon mutation over the full physiological range. Hence, they have important functions in electron transfer and oxidation chemistry and have applications in industrial biotechnology. The details of what determines these reduction potential changes upon mutation are still unclear. Moreover, it has been difficult to model and predict the reduction potential of azurin mutants and currently no unique procedure or workflow pattern exists. Furthermore, high‐level computational methods can be accurate but are too time consuming for practical use. In this work, a novel approach for calculating reduction potentials of azurin mutants is shown, based on a combination of continuum electrostatics, density functional theory and empirical hydrophobicity factors. Our method accurately reproduces experimental reduction potential changes of 30 mutants with respect to wildtype within experimental error and highlights the factors contributing to the reduction potential change. Finally, reduction potentials are predicted for a series of 124 new mutants that have not yet been investigated experimentally. Several mutants are identified that are located well over 10 Å from the copper center that change the reduction potential by more than 85 mV. The work shows that secondary coordination sphere mutations mostly lead to long‐range electrostatic changes and hence can be modeled accurately with continuum electrostatics.

Oxidative Unfolding of the Rubredoxin Domain and the Natively Disordered N-terminal Region Regulate the Catalytic Activity of Mycobacterium tuberculosis Protein Kinase G

Mycobacterium tuberculosis escapes killing in human macrophages by secreting protein kinase G (PknG). PknG intercepts host signaling to prevent fusion of the phagosome engulfing the mycobacteria with the lysosome and, thus, their degradation. The N-terminal NORS (no regulatory secondary structure) region of PknG (approximately residues 1-75) has been shown to play a role in PknG regulation by (auto)phosphorylation, whereas the following rubredoxin-like metal-binding motif (RD, residues ∼74-147) has been shown to interact tightly with the subsequent catalytic domain (approximately residues 148-420) to mediate its redox regulation. Deletions or mutations in NORS or the redox-sensitive RD significantly decrease PknG survival function. Based on combined NMR spectroscopy, in vitro kinase assay, and molecular dynamics simulation data, we provide novel insights into the regulatory roles of the N-terminal regions. The NORS region is indeed natively disordered and rather dynamic. Consistent with most earlier data, autophosphorylation occurs in our assays only when the NORS region is present and, thus, in the NORS region. Phosphorylation of it results only in local conformational changes and does not induce interactions with the subsequent RD. Although the reduced, metal-bound RD makes tight interactions with the following catalytic domain in the published crystal structures, it can also fold in its absence. Our data further suggest that oxidation-induced unfolding of the RD regulates substrate access to the catalytic domain and, thereby, PknG function under different redox conditions, e.g. when exposed to increased levels of reactive oxidative species in host macrophages.

Mysteries of metals in metalloenzymes

Natural metalloenzymes are often the most proficient catalysts in terms of their activity, selectivity, and ability to operate at mild conditions. However, metalloenzymes are occasionally surprising in their selection of catalytic metals, and in their responses to metal substitution. Indeed, from the isolated standpoint of producing the best catalyst, a chemist designing from first-principles would likely choose a different metal. For example, some enzymes employ a redox active metal where a simple Lewis acid is needed. Such are several hydrolases. In other cases, substitution of a non-native metal leads to radical improvements in reactivity. For example, histone deacetylase 8 naturally operates with Zn(2+) in the active site but becomes much more active with Fe(2+). For β-lactamases, the replacement of the native Zn(2+) with Ni(2+) was suggested to lead to higher activity as predicted computationally. There are also intriguing cases, such as Fe(2+)- and Mn(2+)-dependent ribonucleotide reductases and W(4+)- and Mo(4+)-dependent DMSO reductases, where organisms manage to circumvent the scarcity of one metal (e.g., Fe(2+)) by creating protein structures that utilize another metal (e.g., Mn(2+)) for the catalysis of the same reaction. Naturally, even though both metal forms are active, one of the metals is preferred in every-day life, and the other metal variant remains dormant until an emergency strikes in the cell. These examples lead to certain questions. When are catalytic metals selected purely for electronic or structural reasons, implying that enzymatic catalysis is optimized to its maximum? When are metal selections a manifestation of competing evolutionary pressures, where choices are dictated not just by catalytic efficiency but also by other factors in the cell? In other words, how can enzymes be improved as catalysts merely through the use of common biological building blocks available to cells? Addressing these questions is highly relevant to the enzyme design community, where the goal is to prepare maximally efficient quasi-natural enzymes for the catalysis of reactions that interest humankind. Due to competing evolutionary pressures, many natural enzymes may not have evolved to be ideal catalysts and can be improved for the isolated purpose of catalysis in vitro when the competing factors are removed. The goal of this Account is not to cover all the possible stories but rather to highlight how variable enzymatic catalysis can be. We want to bring up possible factors affecting the evolution of enzyme structure, and the large- and intermediate-scale structural and electronic effects that metals can induce in the protein, and most importantly, the opportunities for optimization of these enzymes for catalysis in vitro.

Computational electrochemistry: prediction of liquid-phase reduction potentials

This article reviews recent developments and applications in the area of computational electrochemistry. Our focus is on predicting the reduction potentials of electron transfer and other electrochemical reactions and half-reactions in both aqueous and nonaqueous solutions. Topics covered include various computational protocols that combine quantum mechanical electronic structure methods (such as density functional theory) with implicit-solvent models, explicit-solvent protocols that employ Monte Carlo or molecular dynamics simulations (for example, Car-Parrinello molecular dynamics using the grand canonical ensemble formalism), and the Marcus theory of electronic charge transfer. We also review computational approaches based on empirical relationships between molecular and electronic structure and electron transfer reactivity. The scope of the implicit-solvent protocols is emphasized, and the present status of the theory and future directions are outlined.

Pharmacophore modeling and 3D quantitative structure-activity relationship analysis of febrifugine analogues as potent antimalarial agent

Febrifugine and its derivatives are effective against Plasmodium falciparum. Using PHASE algorithm, a five-point pharmacophore model with two hydrogen bond acceptor (A), one positively ionizable (P) and two aromatic rings (R), was developed to derive a predictive ligand-based statistically significant 3D-quantitative structure-activity relationship (QSAR) model (r² = 0.972, SD = 0.3, F = 173.4, Q² = 0.712, RMSE = 0.3, Person-R = 0.94, and r²_pred = 0.8) to explicate the structural attributes crucial for antimalarial activity. The developed pharmacophore model and 3D QSAR model can be a substantial tool for virtual screening and related antimalarial drug discovery research.

Parameters for molecular dynamics simulations of iron-sulfur proteins

Iron-sulfur proteins involved in electron transfer reactions have finely tuned redox potentials, which allow them to be highly efficient and specific. Factors such as metal center solvent exposure, interaction with charged residues, or hydrogen bonds between the ligand residues and amide backbone groups have all been pointed out to cause such specific redox potentials. Here, we derived parameters compatible with the AMBER force field for the metal centers of iron-sulfur proteins and applied them in the molecular dynamics simulations of three iron-sulfur proteins. We used density-functional theory (DFT) calculations and Seminario's method for the parameterization. Parameter validation was obtained by matching structures and normal frequencies at the quantum mechanics and molecular mechanics levels of theory. Having guaranteed a correct representation of the protein coordination spheres, the amide H-bonds and the water exposure to the ligands were analyzed. Our results for the pattern of interactions with the metal centers are consistent to those obtained by nuclear magnetic resonance spectroscopy (NMR) experiments and DFT calculations, allowing the application of molecular dynamics to the study of those proteins.

Hybrid dynamics simulation engine for metalloproteins

Quality computational description of metalloproteins is a great challenge due to the vast span of time- and lengthscales characteristic of their existence. We present an efficient new method that allows for robust characterization of metalloproteins. It combines quantum mechanical (QM) description of the metal-containing active site, and extensive dynamics of the protein captured by discrete molecular dynamics (DMD) (QM/DMD). DMD samples the entire protein, including the backbone, and most of the active site, except for the immediate coordination region of the metal. QM operates on the part of the protein of electronic and chemical significance, which may include tens to hundreds of atoms. The breathing quantum-classical boundary provides a continuous mutual feedback between the two machineries. We test QM/DMD using the Fe-containing electron transporter protein, rubredoxin, and its three mutants as a model. QM/DMD can provide a reliable balanced description of metalloproteins' structure, dynamics, and electronic structure in a reasonable amount of time. As an illustration of QM/DMD capabilities, we then predict the structure of the Ca(2+) form of the enzyme catechol O-methyl transferase, which, unlike the native Mg(2+) form, is catalytically inactive. The Mg(2+) site is ochtahedral, but the Ca(2+) is 7-coordinate and features the misalignment of the reacting parts of the system. The change is facilitated by the backbone adjustment. QM/DMD is ideal and fast for providing this level of structural insight.