The medium reorganization energy for the charge transfer reactions in proteins

Interdomain Interactions Modulate the Active Site Dynamics of Human Inducible Nitric Oxide Synthase

Nitric oxide synthase (NOS) is a homodimeric flavohemoprotein responsible for catalyzing the oxidation of l-arginine (l-Arg) to citrulline and nitric oxide. Electrons are supplied for the reaction via interdomain electron transfer between an N-terminal heme-containing oxygenase domain and a FMN-containing (sub)domain of a C-terminal reductase domain. Extensive attention has focused on elucidating how conformational dynamics regulate electron transfer between the domains. Here we investigate the impact of the interdomain FMN-heme interaction on the heme active site dynamics of inducible NOS (iNOS). Steady state linear and time-resolved two-dimensional infrared (2D IR) spectroscopy was applied to probe a CO ligand at the heme within the oxygenase domain for full-length and truncated or mutated constructs of human iNOS. Whereas the linear IR spectra of the CO ligand were identical among the constructs, 2D IR spectroscopy revealed variation in the frequency dynamics. The wild-type constructs that can properly form the FMN/oxygenase docked state due to the presence of both the FMN and oxygenase domains showed slower dynamics than the oxygenase domain alone. Introduction of the mutation (E546N) predicted to perturb electrostatic interactions between the domains resulted in measured dynamics intermediate between those for the full-length and individual oxygenase domain, consistent with perturbation to the docked/undocked equilibrium. These results indicate that docking of the FMN domain to the oxygenase domain not only brings the FMN cofactor within electron transfer distance of the heme domain but also modulates the dynamics sensed by the CO ligand within the active site in a way expected to promote efficient electron transfer.

Methemoglobin formation in mutant hemoglobin α chains: electron transfer parameters and rates

Hemoglobin-mediated transport of dioxygen (O₂) critically depends on the stability of the reduced (Fe²⁺) form of the heme cofactors. Some protein mutations stabilize the oxidized (Fe³⁺) state (methemoglobin, Hb M), causing methemoglobinemia, and can be lethal above 30%. The majority of the analyses of factors influencing Hb oxidation are retrospective and give insights only for inner-sphere mutations of heme (His58, His87). Herein, we report the first all-atom molecular dynamics simulations on both redox states and calculations of the Marcus electron transfer (ET) parameters for the α chain Hb oxidation and reduction rates for Hb M. The Hb wild-type (WT) and most of the studied α chain variants maintain globin structure except the Hb M Iwate (H87Y). The mutants forming Hb M tend to have lower redox potentials and thus stabilize the oxidized (Fe³⁺) state (in particular, the Hb Miyagi variant with K61E mutation). Solvent reorganization (λ_solv 73-96%) makes major contributions to reorganization free energy, whereas protein reorganization (λ_prot) accounts for 27-30% except for the Miyagi and J-Buda variants (λ_prot ∼4%). Analysis of heme-solvent H-bonding interactions among variants provide insights into the role of Lys61 residue in stabilizing the Fe²⁺ state. Semiclassical Marcus ET theory-based calculations predict experimental k_ET for the Cyt b5-Hb complex and provide insights into relative reduction rates for Hb M in Hb variants. Thus, our methodology provides a rationale for the effect of mutations on the structure, stability, and Hb oxidation reduction rates and has potential for identification of mutations that result in methemoglobinemia.

Active-Site Environmental Factors Customize the Photophysics of Photoenzymatic Old Yellow Enzymes

The development of non-natural photoenzymatic systems has reinvigorated the study of photoinduced electron transfer (ET) within protein active sites, providing new and unique platforms for understanding how biological environments affect photochemical processes. In this work, we use ultrafast spectroscopy to compare the photoinduced electron transfer in known photoenzymes. 12-Oxophytodienoate reductase 1 (OPR1) is compared to Old Yellow Enzyme 1 (OYE1) and morphinone reductase (MR). The latter enzymes are structurally homologous to OPR1. We find that slight differences in the amino acid composition of the active sites of these proteins determine their distinct electron-transfer dynamics. Our work suggests that the inside of a protein active site is a complex/heterogeneous dielectric network where genetically programmed heterogeneity near the site of biological ET can significantly affect the presence and lifetime of various intermediate states. Our work motivates additional tunability of Old Yellow Enzyme active-site reorganization energy and electron-transfer energetics that could be leveraged for photoenzymatic redox approaches.

Kinetic Isotope Effects Quantify pH-Sensitive Water Dynamics at the Pt Electrode Interface

The pH-dependent kinetics of the hydrogen oxidation and evolution reactions (HERs and HORs) remain a fundamental conundrum in modern electrochemistry. Recent efforts have focused on the impact of the interfacial water network on the reaction kinetics. In this work, we quantify the importance of interfacial water dynamics on the overall hydrogen reaction kinetics with kinetic isotope effect (KIE) voltammetry experiments on single-crystal Pt(111) and Pt(110). Our results find a surface-sensitive KIE for both the HER and the HOR that is measurable in base but not in acid. Remarkably, the HOR in KOD on Pt(111) yields a KIE of up to 3.4 at moderate overpotentials, greater than any expected secondary KIE values, yet the HOR in DClO₄ yields no measurable KIE. These results provide direct evidence that solvent dynamics play a crucial role in the alkaline but not in the acidic hydrogen reactions, thus reinforcing the importance of "beyond adsorption" phenomena in modern electrocatalysis.

Towards direct enzyme wiring: a theoretical investigation of charge carrier transfer mechanisms between glucose oxidase and organic semiconductors

In this work, a general theoretical and numerical approach based on semiconductor theory, which could be applied to a study of direct enzyme wiring, has been discussed. Marcus-Hush theory was applied to evaluate the potential transfer of charge carriers (holes and electrons) between glucose oxidase (GOx) and organic semiconductors. Two mechanisms of multistep hopping of charge to/from the oxidised/reduced flavin-based moiety through residues of aromatic amino acids located in GOx and long range charge direct tunnelling from the cofactor to the organic semiconductor surface have been proposed and evaluated. It was determined that the hole-hopping mechanism is possible and proceeds at a low ionization potential of the organic semiconductor. The calculations reveal that hopping of electrons is blocked, but direct electron tunnelling between the cofactor and the organic semiconductor is still probable. The most optimal conditions and tunable characteristics of GOx-based biosensors such as the ionization potential, electron affinity of organic semiconductors and distance between the enzyme and surface were estimated for the first time.

Reorganization energies of the electron transfer reactions involving quinones in the reaction center of Rhodobacter sphaeroides

In framework of the continuum electrostatics theory, the reorganization energies of the electron transfers Q_A^--Q_B (fast phase), Bph^--Q_A, P⁺-Q_A^-, and P⁺-Q_B^- in the photosynthetic bacterial reaction center have been calculated. The calculations were based on the static dielectric permittivity spatial distribution derived from the data on the electrogenesis, with the corresponding characteristic times relatively close to the reaction times of Q_A^--Q_B (fast phase) and Bph^--Q_A but much shorter than those times of the latter two recombination reactions. The calculated reorganization energies were reasonably close to the experimental estimates for Q_A^--Q_B (fast phase) and Bph^--Q_A but substantially lower than those of P⁺-Q_A^- and P⁺-Q_B^-. A higher effective dielectric permittivity contributes to this effect, but the dominant contribution is most probably made by a non-dielectric relaxation, especially for the P⁺-Q_B^- recombination influenced by the proton transfer. This situation calls for reconsidering of the current electron transfer rate estimates.

Simulating biological charge transfer: Continuum dielectric theory or molecular dynamics?

We study the thermodynamic parameters of Marcus's theory of charge transfer, the driving forces and the reorganization energies, using two widely applied approaches to bioenergetic problems that seem to be radically different: continuum dielectric theory via a numerical solution of Poisson's equation, and the thermodynamic integration approach based upon classical Newtonian molecular dynamics, as perfomed by Na et al., PCCP 19, 18,938 (2017). With application to a nitrite reductase NrfHA protein heterodimer, we obtain an excellent agreement between the respective driving forces with an r.m.s. deviation of 1.7 kcal/mol, and a lower limit to the reorganization energies. The computational methods turn out to be mutually supportive: molecular dynamics can be used to determine the parameters of a dielectric theory computation, which on the other hand can be used to properly rescale the reorganization energies and partition them into aqueous and protein contributions. In addition, we use the electrostatic approach to study the influence of Ca²⁺ ions on the free energy landscape of charge transfer.

Photochemical Generation of a Tryptophan Radical within the Subunit Interface of Ribonucleotide Reductase

The Escherichia coli class Ia ribonucleotide reductase (RNR) achieves forward and reverse proton-coupled electron transfer (PCET) over a pathway of redox active amino acids (β-Y122 ⇌ β-Y356 ⇌ α-Y731 ⇌ α-Y730 ⇌ α-C439) spanning ∼35 Å and two subunits every time it turns over. We have developed photoRNRs that allow radical transport to be phototriggered at tyrosine (Y) or fluorotyrosine (FnY) residues along the PCET pathway. We now report a new photoRNR in which photooxidation of a tryptophan (W) residue replacing Y356 within the α/β subunit interface proceeds by a stepwise ET/PT (electron transfer then proton transfer) mechanism and provides an orthogonal spectroscopic handle with respect to radical pathway residues Y731 and Y730 in α. This construct displays an ∼3-fold enhancement in photochemical yield of W(•) relative to F3Y(•) and a ∼7-fold enhancement relative to Y(•). Photogeneration of the W(•) radical occurs with a rate constant of (4.4 ± 0.2) × 10(5) s(-1), which obeys a Marcus correlation for radical generation at the RNR subunit interface. Despite the fact that the Y → W variant displays no enzymatic activity in the absence of light, photogeneration of W(•) within the subunit interface results in 20% activity for turnover relative to wild-type RNR under the same conditions.

Single-Molecule Conductance of Viologen-Cucurbit[8]uril Host-Guest Complexes

The local molecular environment is a critical factor which should be taken into account when measuring single-molecule electrical properties in condensed media or in the design of future molecular electronic or single molecule sensing devices. Supramolecular interactions can be used to control the local environment in molecular assemblies and have been used to create microenvironments, for instance, for chemical reactions. Here, we use supramolecular interactions to create microenvironments which influence the electrical conductance of single molecule wires. Cucurbit[8]uril (CB[8]) with a large hydrophobic cavity was used to host the viologen (bipyridinium) molecular wires forming a 1:1 supramolecular complex. Significant increases in the viologen wire single molecule conductances are observed when it is threaded into CB[8] due to large changes of the molecular microenvironment. The results were interpreted within the framework of a Marcus-type model for electron transfer as arising from a reduction in outer-sphere reorganization energy when the viologen is confined within the hydrophobic CB[8] cavity.

An Empirical Polarizable Force Field Based on the Classical Drude Oscillator Model: Development History and Recent Applications

Molecular mechanics force fields that explicitly account for induced polarization represent the next generation of physical models for molecular dynamics simulations. Several methods exist for modeling induced polarization, and here we review the classical Drude oscillator model, in which electronic degrees of freedom are modeled by charged particles attached to the nuclei of their core atoms by harmonic springs. We describe the latest developments in Drude force field parametrization and application, primarily in the last 15 years. Emphasis is placed on the Drude-2013 polarizable force field for proteins, DNA, lipids, and carbohydrates. We discuss its parametrization protocol, development history, and recent simulations of biologically interesting systems, highlighting specific studies in which induced polarization plays a critical role in reproducing experimental observables and understanding physical behavior. As the Drude oscillator model is computationally tractable and available in a wide range of simulation packages, it is anticipated that use of these more complex physical models will lead to new and important discoveries of the physical forces driving a range of chemical and biological phenomena.

Charge-Transfer Dynamics at the α/β Subunit Interface of a Photochemical Ribonucleotide Reductase

Ribonucleotide reductase (RNR) catalyzes the conversion of ribonucleotides to deoxyribonucleotides to provide the monomeric building blocks for DNA replication and repair. Nucleotide reduction occurs by way of multistep proton-coupled electron transfer (PCET) over a pathway of redox active amino acids spanning ∼35 Å and two subunits (α2 and β2). Despite the fact that PCET in RNR is rapid, slow conformational changes mask examination of the kinetics of these steps. As such, we have pioneered methodology in which site-specific incorporation of a [Re(I)] photooxidant on the surface of the β2 subunit (photoβ2) allows photochemical oxidation of the adjacent PCET pathway residue β-Y356 and time-resolved spectroscopic observation of the ensuing reactivity. A series of photoβ2s capable of performing photoinitiated substrate turnover have been prepared in which four different fluorotyrosines (FnYs) are incorporated in place of β-Y356. The FnYs are deprotonated under biological conditions, undergo oxidation by electron transfer (ET), and provide a means by which to vary the ET driving force (ΔG°) with minimal additional perturbations across the series. We have used these features to map the correlation between ΔG° and kET both with and without the fully assembled photoRNR complex. The photooxidation of FnY356 within the α/β subunit interface occurs within the Marcus inverted region with a reorganization energy of λ ≈ 1 eV. We also observe enhanced electronic coupling between donor and acceptor (HDA) in the presence of an intact PCET pathway. Additionally, we have investigated the dynamics of proton transfer (PT) by a variety of methods including dependencies on solvent isotopic composition, buffer concentration, and pH. We present evidence for the role of α2 in facilitating PT during β-Y356 photooxidation; PT occurs by way of readily exchangeable positions and within a relatively "tight" subunit interface. These findings show that RNR controls ET by lowering λ, raising HDA, and directing PT both within and between individual polypeptide subunits.

Converting the bis-FeIV state of the diheme enzyme MauG to Compound I decreases the reorganization energy for electron transfer

The electron transfer (ET) properties of two types of high-valent hemes were studied within the same protein matrix; the bis-Fe(IV) state of MauG and the Compound I state of Y294H MauG. The latter is formed as a consequence of mutation of the tyrosine which forms the distal axial ligand of the six-coordinate heme that allows it to stabilize Fe(IV) in the absence of an external ligand. The rates of the ET reaction of each high-valent species with the type I copper protein, amicyanin, were determined at different temperatures and analysed by ET theory. The reaction with bis-Fe(IV) wild-type (WT) MauG exhibited a reorganization energy (λ) that was 0.39 eV greater than that for the reaction of Compound I Y295H MauG. It is concluded that the delocalization of charge over the two hemes in the bis-Fe(IV) state is responsible for the larger λ, relative to the Compound I state in which the Fe(V) equivalent is isolated on one heme. Although the increase in λ decreases the rate of ET, the delocalization of charge decreases the ET distance to its natural substrate protein, thus increasing the ET rate. This describes how proteins can balance different ET properties of complex redox cofactors to optimize each system for its particular ET or catalytic reaction.

Recent Advances in Polarizable Force Fields for Macromolecules: Microsecond Simulations of Proteins Using the Classical Drude Oscillator Model

In this Perspective, we summarize recent efforts to include the explicit treatment of induced electronic polarization in biomolecular force fields. Methods used to treat polarizability, including the induced dipole, fluctuating charge, and classical Drude oscillator models, are presented, including recent advances in force fields using those methods. This is followed by recent results obtained with the Drude model, including microsecond molecular dynamics (MD) simulations of multiple proteins in explicit solvent. Results show significant variability of backbone and side-chain dipole moments as a function of environment, including significant changes during individual simulations. Dipole moments of water in the vicinity of the proteins reveal small but systematic changes, with the direction of the changes dependent on the environment. Analyses of the full proteins show that the polarizable Drude model leads to larger values of the dielectric constant of the protein interior, especially in the case of hydrophobic regions. These results indicate that the inclusion of explicit electronic polarizability leads to significant differences in the physical forces affecting the structure and dynamics of proteins, which can be investigated in a computationally tractable fashion in the context of the Drude model.

The mechanism of ubihydroquinone oxidation at the Qo-site of the cytochrome bc1 complex

1. Recent results suggest that the major flux is carried by a monomeric function, not by an intermonomer electron flow. 2. The bifurcated reaction at the Qo-site involves sequential partial processes, - a rate limiting first electron transfer generating a semiquinone (SQ) intermediate, and a rapid second electron transfer in which the SQ is oxidized by the low potential chain. 3. The rate constant for the first step in a strongly endergonic, proton-first-then-electron mechanism, is given by a Marcus-Brønsted treatment in which a rapid electron transfer is convoluted with a weak occupancy of the proton configuration needed for electron transfer. 4. A rapid second electron transfer pulls the overall reaction over. Mutation of Glu-295 of cyt b shows it to be a key player. 5. In more crippled mutants, electron transfer is severely inhibited and the bell-shaped pH dependence of wildtype is replaced by a dependence on a single pK at ~8.5 favoring electron transfer. Loss of a pK ~6.5 is explained by a change in the rate limiting step from the first to the second electron transfer; the pK ~8.5 may reflect dissociation of QH. 6. A rate constant (<10(3)s(-1)) for oxidation of SQ in the distal domain by heme bL has been determined, which precludes mechanisms for normal flux in which SQ is constrained there. 7. Glu-295 catalyzes proton exit through H(+) transfer from QH, and rotational displacement to deliver the H(+) to exit channel(s). This opens a volume into which Q(-) can move closer to the heme to speed electron transfer. 8. A kinetic model accounts well for the observations, but leaves open the question of gating mechanisms. For the first step we suggest a molecular "escapement"; for the second a molecular ballet choreographed through coulombic interactions. This article is part of a Special Issue entitled: Respiratory complex III and related bc complexes.

Molecular mechanisms for generating transmembrane proton gradients

Membrane proteins use the energy of light or high energy substrates to build a transmembrane proton gradient through a series of reactions leading to proton release into the lower pH compartment (P-side) and proton uptake from the higher pH compartment (N-side). This review considers how the proton affinity of the substrates, cofactors and amino acids are modified in four proteins to drive proton transfers. Bacterial reaction centers (RCs) and photosystem II (PSII) carry out redox chemistry with the species to be oxidized on the P-side while reduction occurs on the N-side of the membrane. Terminal redox cofactors are used which have pKas that are strongly dependent on their redox state, so that protons are lost on oxidation and gained on reduction. Bacteriorhodopsin is a true proton pump. Light activation triggers trans to cis isomerization of a bound retinal. Strong electrostatic interactions within clusters of amino acids are modified by the conformational changes initiated by retinal motion leading to changes in proton affinity, driving transmembrane proton transfer. Cytochrome c oxidase (CcO) catalyzes the reduction of O2 to water. The protons needed for chemistry are bound from the N-side. The reduction chemistry also drives proton pumping from N- to P-side. Overall, in CcO the uptake of 4 electrons to reduce O2 transports 8 charges across the membrane, with each reduction fully coupled to removal of two protons from the N-side, the delivery of one for chemistry and transport of the other to the P-side.