Mapping hole hopping escape routes in proteins

Elucidating the Role of O2 Uncoupling for the Adaptation of Bacterial Biodegradation Reactions Catalyzed by Rieske Oxygenases

Oxygenation of aromatic and aliphatic hydrocarbons by Rieske oxygenases is the initial step of various biodegradation pathways for environmental organic contaminants. Microorganisms carrying Rieske oxygenases are able to quickly adapt their substrate spectra to alternative carbon and energy sources that are structurally related to the original target substrate, yet the molecular events responsible for this rapid adaptation are not well understood. Here, we evaluated the hypothesis that reactive oxygen species (ROS) generated by unproductive activation of O2, the so-called O2 uncoupling, in the presence of the alternative substrate exert a selective pressure on the bacterium for increasing the oxygenation efficiency of Rieske oxygenases. To that end, we studied wild-type 2-nitrotoluene dioxygenase from Acidovorax sp. strain JS42 and five enzyme variants that have evolved from adaptive laboratory evolution experiments with 3- and 4-nitrotoluene as alternative growth substrates. The enzyme variants showed a substantially increased oxygenation efficiency toward the new target substrates concomitant with a reduction of ROS production, while mechanisms and kinetics of enzymatic O2 activation remained unchanged. Structural analyses and docking studies suggest that amino acid substitutions in enzyme variants occurred at residues lining both substrate and O2 transport tunnels, enabling tighter binding of the target substrates in the active site. Increased oxygenation efficiencies measured in vitro for the various enzyme (variant)-substrate combinations correlated linearly with in vivo changes in growth rates for evolved Acidovorax strains expressing the variants. Our data suggest that the selective pressure from oxidative stress toward more efficient oxygenation by Rieske oxygenases was most notable when O2 uncoupling exceeded 60%.

An Investigation of the Influence of Tyrosine Local Interactions on Electron Hopping in a Model Protein

Multi-step electron transfer reactions are important to the function of many cellular systems. The ways in which such systems have evolved to direct electrons along specific pathways are largely understood, but less so are the ways in which the reduction-oxidation potentials of individual redox sites are controlled. We prepared a series of three new artificial variants of Pseudomonas aeruginosa azurin where a tyrosine (Tyr109) is situated between the native Cu ion and a Ru(II) photosensitizer tethered to a histidine (His107). Arginine, glutamine, or methionine were introduced as position 122, which is near to Tyr109. We investigated the rate of CuI oxidation by a flash-quench generated Ru(III) oxidant over pH values from 5 to 9. While the identity of the residue at position 122 affects some of the physical properties of Tyr109, the rates of CuI oxidation are only weakly dependent on the identity of the residue at 122. The results highlight that more work is still needed to understand how non-covalent interactions of redox active groups are affected in redox proteins.

New insights into controlling radical migration pathways in heme enzymes gained from the study of a dye-decolorising peroxidase

In heme enzymes, such as members of the dye-decolorising peroxidase (DyP) family, the formation of the highly oxidising catalytic Fe(iv)-oxo intermediates following reaction with hydrogen peroxide can lead to free radical migration (hole hopping) from the heme to form cationic tyrosine and/or tryptophan radicals. These species are highly oxidising (∼1 V vs. NHE) and under certain circumstances can catalyse the oxidation of organic substrates. Factors that govern which specific tyrosine or tryptophan the free radical migrates to in heme enzymes are not well understood, although in the case of tyrosyl radical formation the nearby proximity of a proton acceptor is a recognised facilitating factor. By using an A-type member of the DyP family (DtpAa) as an exemplar, we combine protein engineering, X-ray crystallography, hole-hopping calculations, EPR spectroscopy and kinetic modelling to provide compelling new insights into the control of radical migration pathways following reaction of the heme with hydrogen peroxide. We demonstrate that the presence of a tryptophan/tyrosine dyad motif displaying a T-shaped orientation of aromatic rings on the proximal side of the heme dominates the radical migration landscape in wild-type DtpAa and continues to do so following the rational engineering into DtpAa of a previously identified radical migration pathway in an A-type homolog on the distal side of the heme. Only on disrupting the proximal dyad, through removal of an oxygen atom, does the radical migration pathway then switch to the engineered distal pathway to form the desired tyrosyl radical. Implications for protein design and biocatalysis are discussed.

Long-Range Conductivity in Proteins Mediated by Aromatic Residues

Single-molecule measurements show that many proteins, lacking any redox cofactors, nonetheless exhibit electrical conductance on the order of a nanosiemen over 10 nm distances, implying that electrons can transit an entire protein in less than a nanosecond when subject to a potential difference of less than 1 V. This is puzzling because, for fast transport (i.e., a free energy barrier of zero), the hopping rate is determined by the reorganization energy of approximately 0.8 eV, and this sets the time scale of a single hop to at least 1 μs. Furthermore, the Fermi energies of typical metal electrodes are far removed from the energies required for sequential oxidation and reduction of the aromatic residues of the protein, which should further reduce the hopping current. Here, we combine all-atom molecular dynamics (MD) simulations of non-redox-active proteins (consensus tetratricopeptide repeats) with an electron transfer theory to demonstrate a molecular mechanism that can account for the unexpectedly fast electron transport. According to our MD simulations, the reorganization energy produced by the energy shift on charging (the Stokes shift) is close to the conventional value of 0.8 eV. However, the non-ergodic sampling of molecular configurations by the protein results in reaction-reorganization energies, extracted directly from the distribution of the electrostatic energy fluctuations, that are only ∼0.2 eV, which is small enough to enable long-range conductivity, without invoking quantum coherent transport. Using the MD values of the reorganization energies, we calculate a current decay with distance that is in agreement with experiment.

Control of charge transport in electronically active systems towards integrated biomolecular circuits (IbC)

The miniaturization of traditional silicon-based electronics will soon reach its limitation as quantum tunneling and heat become serious problems at the several-nanometer scale. Crafting integrated circuits via self-assembly of electronically active molecules using a "bottom-up" paradigm provides a potential solution to these technological challenges. In particular, integrated biomolecular circuits (IbC) offer promising advantages to achieve this goal, as nature offers countless examples of functionalities entailed by self-assembly and examples of controlling charge transport at the molecular level within the self-assembled structures. To this end, the review summarizes the progress in understanding how charge transport is regulated in biosystems and the key redox-active amino acids that enable the charge transport. In addition, charge transport mechanisms at different length scales are also reviewed, offering key insights for controlling charge transport in IbC in the future.

Study and design of amino acid-based radical enzymes using unnatural amino acids

Radical enzymes harness the power of reactive radical species by placing them in a protein scaffold, and they are capable of catalysing many important reactions. New native radical enzymes, especially those with amino acid-based radicals, in the category of non-heme iron enzymes (including ribonucleotide reductases), heme enzymes, copper enzymes, and FAD-radical enzymes have been discovered and characterized. We discussed recent research efforts to discover new native amino acid-based radical enzymes, and to study the roles of radicals in processes such as enzyme catalysis and electron transfer. Furthermore, design of radical enzymes in a small and simple scaffold not only allows us to study the radical in a well-controlled system and test our understanding of the native enzymes, but also allows us to create powerful enzymes. In the study and design of amino acid-based radical enzymes, the use of unnatural amino acids allows precise control of pKa values and reduction potentials of the residue, as well as probing the location of the radical through spectroscopic methods, making it a powerful research tool. Our understanding of amino acid-based radical enzymes will allow us to tailor them to create powerful catalysts and better therapeutics.

Adenine Fine-Tunes DNA Photolyase's Repair Mechanism

The comparative study of DNA repair by mesophilic and extremophilic photolyases helps us understand the evolution of these enzymes and their role in preserving life on our changing planet. The mechanism of repair of cyclobutane pyrimidine dimer lesions in DNA by electron transfer from the flavin adenine dinucleotide cofactor is the subject of intense interest. The role of adenine in mediating this process remains unresolved. Using microsecond molecular dynamics simulations, we find that adenine mediates the electron transfer in both mesophile and extremophile DNA photolyases through a similar mechanism. In fact, in all photolyases studied, the molecular conformations with the largest electronic couplings between the enzyme cofactor and DNA show the presence of adenine in 10-20% of the strongest-coupling tunneling pathways between the atoms of the electron donor and acceptor. Our theoretical analysis finds that adenine serves the critical role of fine-tuning rather than maximizing the donor-acceptor coupling within the range appropriate for the repair function.

Profiling charge transport: A new computational approach

To maintain life, charge transfer processes must be efficient to allow electrons to migrate across distances as large as 30-50 Å within a timescale from picoseconds to milliseconds, and the free-energy cost should not exceed one electron volt. By employing local ionization and local affinity energies, we calculated the pathway for electron and electron-hole transport, respectively. The pathway is then used to calculate both the driving force and the activation energy. The electronic coupling is calculated using configuration interaction procedure. When the charge acceptor is not known, as in oxidative stress, the charge transport terminals are found using Monte-Carlo simulation. These parameters were used to calculate the rate described by Marcus theory. Our approach has been elaborately explained using the famous androstane example and then applied to two proteins: electron transport in azurin protein and hole-hopping migration route from the heme center of cytochrome c peroxidase to its surface. This model gives an effective method to calculate the charge transport pathway and the free-energy profile within 0.1 eV from the experimental measurements and electronic coupling within 3 meV.

Theory of Protein Charge Transfer: Electron Transfer between Tryptophan Residue and Active Site of Azurin

One reaction step in the conductivity relay of azurin, electron transfer between the Cu-based active site and the tryptophan residue, is studied theoretically and by classical molecular dynamics simulations. Oxidation of tryptophan results in electrowetting of this residue. This structural change makes the free energy surfaces of electron transfer nonparabolic as described by the Q-model of electron transfer. We analyze the medium dynamical effect on protein electron transfer produced by coupled Stokes-shift dynamics and the dynamics of the donor-acceptor distance modulating electron tunneling. The equilibrium donor-acceptor distance falls in the plateau region of the rate constant, where it is determined by the protein-water dynamics, and the probability of electron tunneling does not affect the rate. The crossover distance found here puts most intraprotein electron-transfer reactions under the umbrella of dynamical control. The crossover between the medium-controlled and tunneling-controlled kinetics is combined with the effect of the protein-water medium on the activation barrier to formulate principles of tunability of protein-based charge-transfer chains. The main principle in optimizing the activation barrier is the departure from the Gaussian-Gibbsian statistics of fluctuations promoting activated transitions. This is achieved either by incomplete (nonergodic) sampling, breaking the link between the Stokes-shift and variance reorganization energies, or through wetting-induced structural changes of the enzyme's active site.

Cofactor Dynamics Couples the Protein Surface to the Heme in Cytochrome c, Facilitating Electron Transfer

Electron transport through biomolecules and in biological transport networks is of great importance to bioenergetics and biocatalysis. More generally, it is of crucial importance to understand how the pathways that connect buried metallocofactors to other cofactors, and to protein surfaces, affect the biological chemistry of metalloproteins. In terms of electron transfer (ET), the strongest coupling pathways usually comprise covalent and hydrogen bonded networks, with a limited number of through-space contacts. Herein, we set out to determine the relative roles of hydrogen bonds involved in ET via an established heme-to-surface tunneling pathway in cytochrome (cyt) c (i.e., heme-W59-D60-E61-N62). A series of cyt c variants were produced where a ruthenium tris(diimine) photooxidant was placed at position 62 via covalent modification of the N62C residue. Surprisingly, variants where the H-bonding residues W59 and D60 were replaced (i.e., W59F and D60A) showed no change in ET rate from the ferrous heme to Ru(III). In contrast, changing the composition of an alternative tunneling pathway (i.e., heme-M64-N63-C62) with the M64L substitution shows a factor of 2 decrease in the rate of heme-to-Ru ET. This pathway involves a through-space tunneling step between the heme and M64 residue, and such steps are usually disfavored. To rationalize why the heme-M64-N63-C62 is preferred, molecular dynamics (MD) simulations and Pathways analysis were employed. These simulations show that the change in heme-Ru ET rates is attributed to different conformations with compressed donor-acceptor distances, by ∼2 Å in pathway distance, in the M64-containing protein as compared to the M64L protein. The change in distance is correlated with changes in the electronic coupling that are in accord with the experimentally observed heme-Ru ET rates. Remarkably, the M64L variation at the core of the protein translates to changes in cofactor dynamics at the protein surface. The surface changes identified by MD simulations include dynamic anion-π and dipole-dipole interactions. These interactions influence the strength of tunneling pathways and ET rates by facilitating decreases in through-space tunneling distances in key coupling pathways.

Insights into the Thermodynamics and Kinetics of Amino-Acid Radicals in Proteins

Some oxidoreductase enzymes use redox-active tyrosine, tryptophan, cysteine, and/or glycine residues as one-electron, high-potential redox (radical) cofactors. Amino-acid radical cofactors typically perform one of four tasks-they work in concert with a metallocofactor to carry out a multielectron redox process, serve as storage sites for oxidizing equivalents, activate the substrate molecules, or move oxidizing equivalents over long distances. It is challenging to experimentally resolve the thermodynamic and kinetic redox properties of a single-amino-acid residue. The inherently reactive and highly oxidizing properties of amino-acid radicals increase the experimental barriers further still. This review describes a family of stable and well-structured model proteins that was made specifically to study tyrosine and tryptophan oxidation-reduction. The so-called α₃X model protein system was combined with very-high-potential protein film voltammetry, transient absorption spectroscopy, and theoretical methods to gain a comprehensive description of the thermodynamic and kinetic properties of protein tyrosine and tryptophan radicals.

Probing the Role of Cysteine Thiyl Radicals in Biology: Eminently Dangerous, Difficult to Scavenge

Thiyl radicals are exceptionally interesting reactive sulfur species (RSS), but rather rarely considered in a biological or medical context. We here review the reactivity of protein thiyl radicals in aqueous and lipid phases and provide an overview of their most relevant reaction partners in biological systems. We deduce that polyunsaturated fatty acids (PUFAs) are their preferred reaction substrates in lipid phases, whereas protein side chains arguably prevail in aqueous phases. In both cellular compartments, a single, dominating thiyl radical-specific antioxidant does not seem to exist. This conclusion is rationalized by the high reaction rate constants of thiyl radicals with several highly concentrated substrates in the cell, precluding effective interception by antioxidants, especially in lipid bilayers. The intractable reactivity of thiyl radicals may account for a series of long-standing, but still startling biochemical observations surrounding the amino acid cysteine: (i) its global underrepresentation on protein surfaces, (ii) its selective avoidance in aerobic lipid bilayers, especially the inner mitochondrial membrane, (iii) the inverse correlation between cysteine usage and longevity in animals, (iv) the mitochondrial synthesis and translational incorporation of cysteine persulfide, and potentially (v) the ex post introduction of selenocysteine into the genetic code.

On the roles of methionine and the importance of its microenvironments in redox metalloproteins

The amino acid residue methionine (Met) is commonly thought of as a ligand in redox metalloproteins, for example in cytochromes c and in blue copper proteins. However, the roles of Met can go beyond a simple ligand. The thioether functional group of Met allows it to be considered as a hydrophobic residue as well as one that is capable of weak dipolar interactions. In addition, the lone pairs on sulphur allow Met to interact with other groups, inluding the aforementioned metal ions. Because of its properties, Met can play diverse roles in metal coordination, fine tuning of redox reactions, or supporting protein structures. These roles are strongly influenced by the nature of the surrounding medium. Herein, we describe several common interactions between Met and surrounding aromatic amino acids and how they affect the physical properties of both copper and iron metalloproteins. While the importance of interactions between Met and other groups is established in biological systems, less is known about their roles in redox metalloproteins and our view is that this is an area that is ready for greater attention.

Clustering of Aromatic Amino Acid Residues around Methionine in Proteins

Short-range, non-covalent interactions between amino acid residues determine protein structures and contribute to protein functions in diverse ways. The interactions of the thioether of methionine with the aromatic rings of tyrosine, tryptophan, and/or phenylalanine has long been discussed and such interactions are favorable on the order of 1–3 kcal mol⁻¹. Here, we carry out a new bioinformatics survey of known protein structures where we assay the propensity of three aromatic residues to localize around the [-CH₂-S-CH₃] of methionine. We term these groups “3-bridge clusters”. A dataset consisting of 33,819 proteins with less than 90% sequence identity was analyzed and such clusters were found in 4093 structures (or 12% of the non-redundant dataset). All sub-classes of enzymes were represented. A 3D coordinate analysis shows that most aromatic groups localize near the CH₂ and CH₃ of methionine. Quantum chemical calculations support that the 3-bridge clusters involve a network of interactions that involve the Met-S, Met-CH₂, Met-CH₃, and the π systems of nearby aromatic amino acid residues. Selected examples of proposed functions of 3-bridge clusters are discussed.

Functional and protective hole hopping in metalloenzymes

Electrons can tunnel through proteins in microseconds with a modest release of free energy over distances in the 15 to 20 Å range. To span greater distances, or to move faster, multiple charge transfers (hops) are required. When one of the reactants is a strong oxidant, it is convenient to consider the movement of a positively charged "hole" in a direction opposite to that of the electron. Hole hopping along chains of tryptophan (Trp) and tyrosine (Tyr) residues is a critical function in several metalloenzymes that generate high-potential intermediates by reactions with O2 or H2O2, or by activation with visible light. Examination of the protein structural database revealed that Tyr/Trp chains are common protein structural elements, particularly among enzymes that react with O2 and H2O2. In many cases these chains may serve a protective role in metalloenzymes by deactivating high-potential reactive intermediates formed in uncoupled catalytic turnover.

Multi-state harmonic models with globally shared bath for nonadiabatic dynamics in the condensed phase

Model Hamiltonians constructed from quantum chemistry calculations and molecular dynamics simulations are widely used for simulating nonadiabatic dynamics in the condensed phase. The most popular two-state spin-boson model could be built by mapping the all-atom anharmonic Hamiltonian onto a two-level system bilinearly coupled to a harmonic bath using the energy gap time correlation function. However, for more than two states, there lacks a general strategy to construct multi-state harmonic (MSH) models since the energy gaps between different pairs of electronic states are not entirely independent and need to be considered consistently. In this paper, we extend the previously proposed approach for building three-state harmonic models for photoinduced charge transfer to the arbitrary number of electronic states with a globally shared bath and the system-bath couplings are scaled differently according to the reorganization energies between each pair of states. We demonstrate the MSH model construction for an organic photovoltaic carotenoid-porphyrin-C₆₀ molecular triad dissolved in explicit tetrahydrofuran solvent. Nonadiabatic dynamics was simulated using mixed quantum-classical techniques, including the linearized semiclassical and symmetrical quasiclassical dynamics with the mapping Hamiltonians, mean-field Ehrenfest, and mixed quantum-classical Liouville dynamics in two-state, three-state, and four-state harmonic models of the triad system. The MSH models are shown to provide a general and flexible framework for simulating nonadiabatic dynamics in complex systems.

Photoinduced hole hopping through tryptophans in proteins

Hole hopping through tryptophan/tyrosine chains enables rapid unidirectional charge transport over long distances. We have elucidated structural and dynamical factors controlling hopping speed and efficiency in two modified azurin constructs that include a rhenium(I) sensitizer, Re(His)(CO)₃(dmp)⁺, and one or two tryptophans (W₁, W₂). Experimental kinetics investigations showed that the two closely spaced (3 to 4 Å) intervening tryptophans dramatically accelerated long-range electron transfer (ET) from Cu^I to the photoexcited sensitizer. In our theoretical work, we found that time-dependent density-functional theory (TDDFT) quantum mechanics/molecular mechanics/molecular dynamics (QM/MM/MD) trajectories of low-lying triplet excited states of Re^I(His)(CO)₃(dmp)⁺-W₁(-W₂) exhibited crossings between sensitizer-localized (*Re) and charge-separated [Re^I(His)(CO)₃(dmp^•-)/(W₁ ^•+ or W₂ ^•+)] (CS1 or CS2) states. Our analysis revealed that the distances, angles, and mutual orientations of ET-active cofactors fluctuate in a relatively narrow range in which the cofactors are strongly coupled, enabling adiabatic ET. Water-dominated electrostatic field fluctuations bring *Re and CS1 states to a crossing where *Re(CO)₃(dmp)⁺←W₁ ET occurs, and CS1 becomes the lowest triplet state. ET is promoted by solvation dynamics around *Re(CO)₃(dmp)⁺(W₁); and CS1 is stabilized by Re(dmp^•-)/W₁ ^•+ electron/hole interaction and enhanced W₁ ^•+ solvation. The second hop, W₁ ^•+←W₂, is facilitated by water fluctuations near the W₁/W₂ unit, taking place when the electrostatic potential at W₂ drops well below that at W₁ ^•+ Insufficient solvation and reorganization around W₂ make W₁ ^•+←W₂ ET endergonic, shifting the equilibrium toward W₁ ^•+ and decreasing the charge-separation yield. We suggest that multiscale TDDFT/MM/MD is a suitable technique to model the simultaneous evolution of photogenerated excited-state manifolds.

An interesting possibility of forming special hole stepping stones with high-stacking aromatic rings in proteins: three-π five-electron and four-π seven-electron resonance bindings

Long-range hole transfer of proteins plays an important role in many biological processes of living organisms. Therefore, it is highly useful to examine the possible hole stepping stones, which can facilitate hole transfer in proteins. However, the structures of stepping stones are diverse because of the complexity of the protein structures. In the present work, we proposed a series of special stepping stones, which are instantaneously formed by three and four packing aromatic side chains of amino acids to capture a hole, corresponding to three-π five-electron (π:π∴π↔π∴π:π) and four-π seven-electron (π:π∴π:π↔π:π:π∴π) resonance bindings with appropriate binding energies. The aromatic amino acids include histidine (His), phenylalanine (Phe), tyrosine (Tyr) and tryptophan (Trp). The formations of these special stepping stones can effectively reduce the local ionization potential of the high π-stacking region to efficiently capture the migration hole. The quick formations and separations of them promote the efficient hole transfer in proteins. More interestingly, we revealed that a hole cannot delocalize over infinite aromatic rings along the high π-π packing structure at the same time and the micro-surroundings of proteins can modulate the formations of π:π∴π↔π∴π:π and π:π∴π:π↔π:π:π∴π bindings. These results may contribute a new avenue to better understand the potential hole transfer pathway in proteins.

Why Do Most Aromatics Fail to Support Hole Hopping in the Cytochrome c Peroxidase-Cytochrome c Complex?

Electron transport through aromatic species (especially tryptophan and tyrosine) plays a central role in water splitting, redox signaling, oxidative damage protection, and bioenergetics. The cytochrome c peroxidase (CcP)-cytochrome c (Cc) complex (CcP:Cc) is used widely to study interprotein electron transfer (ET) mechanisms. Tryptophan 191 (Trp191) of CcP supports hole hopping charge recombination in the CcP:Cc complex. Experimental studies find that when Trp191 is substituted by tyrosine, phenylalanine, or redox-active aniline derivatives bound in the W191G cavity, enzymatic activity and charge recombination rates both decrease. Theoretical analysis of these CcP:Cc complexes finds that the ET kinetics depend strongly on the chemistry of the modified Trp site. The computed electronic couplings in the W191F and W191G species are orders of magnitude smaller than in the native protein, due largely to the absence of a hopping intermediate and the large tunneling distance. Small molecules bound in the W191G cavity are weakly coupled electronically to the Cc heme, and the structural disorder of the guest molecule in the binding pocket may contribute further to the lack of enzymatic activity. The couplings in W191Y are not substantially weakened compared to the native species, but the redox potential difference for tyrosine vs tryptophan oxidation accounts for the slower rate in the Tyr mutant. Thus, theoretical analysis explains why only the native Trp supports rapid hole hopping in the CcP:Cc complex. Favorable free energies and electronic couplings are essential for establishing an efficient hole hopping relay in this protein-protein complex.

Oxalate decarboxylase uses electron hole hopping for catalysis

The hexameric low-pH stress response enzyme oxalate decarboxylase catalyzes the decarboxylation of the oxalate mono-anion in the soil bacterium Bacillus subtilis. A single protein subunit contains two Mn-binding cupin domains, and catalysis depends on Mn(III) at the N-terminal site. The present study suggests a mechanistic function for the C-terminal Mn as an electron hole donor for the N-terminal Mn. The resulting spatial separation of the radical intermediates directs the chemistry toward decarboxylation of the substrate. A π-stacked tryptophan pair (W96/W274) links two neighboring protein subunits together, thus reducing the Mn-to-Mn distance from 25.9 Å (intrasubunit) to 21.5 Å (intersubunit). Here, we used theoretical analysis of electron hole-hopping paths through redox-active sites in the enzyme combined with site-directed mutagenesis and X-ray crystallography to demonstrate that this tryptophan pair supports effective electron hole hopping between the C-terminal Mn of one subunit and the N-terminal Mn of the other subunit through two short hops of ∼8.5 Å. Replacement of W96, W274, or both with phenylalanine led to a large reduction in catalytic efficiency, whereas replacement with tyrosine led to recovery of most of this activity. W96F and W96Y mutants share the wildtype tertiary structure. Two additional hole-hopping networks were identified leading from the Mn ions to the protein surface, potentially protecting the enzyme from high Mn oxidation states during turnover. Our findings strongly suggest that multistep hole-hopping transport between the two Mn ions is required for enzymatic function, adding to the growing examples of proteins that employ aromatic residues as hopping stations.

Amino-acid interactions with the Au(111) surface: adsorption, band alignment, and interfacial electronic coupling

The charge transport properties of biological molecules like peptides and proteins are intensively studied for the great flexibility, redox-state variability, long-range efficiency, and biocompatibility of potential bioelectronic applications. Yet, the electronic interactions of biomolecules with solid metal surfaces, determining the conductivities of the biomolecular junctions, are hard to predict and usually unavailable. Here, we present accurate adsorption structures and energies, electronic band alignment, and interfacial electronic coupling data for all 20 natural amino acids computed using the DFT+Σ scheme based on the vdW-DF and OT-RSH functionals. For comparison, data obtained using the popular PBE functional are provided as well. Tryptophan, compared to other amino acids, is shown to be distinctly exceptional in terms of the electronic properties related to charge transport. Its high adsorption energy, frontier-orbital levels aligned relatively close to the Fermi energy of gold and strong interfacial electronic coupling make it an ideal candidate for facilitating charge transfer on such heterogeneous interfaces. Although the amino acids in peptides and proteins are affected by the structural interactions hindering their contact with the surface, knowledge of the single-molecule surface interactions is necessary for a detailed understanding of such structural effects and tuning of potential applications.

Peripheral Methionine Residues Impact Flavin Photoreduction and Protonation in an Engineered LOV Domain Light Sensor

Proton-coupled electron transfer reactions play critical roles in many aspects of sensory phototransduction. In the case of flavoprotein light sensors, reductive quenching of flavin excited states initiates chemical and conformational changes that ultimately transmit light signals to downstream targets. These reactions generally require neighboring aromatic residues and proton-donating side chains for rapid and coordinated electron and proton transfer to flavin. Although photoreduction of flavoproteins can produce either the anionic (ASQ) or neutral semiquinone (NSQ), the factors that favor one over the other are not well understood. Here we employ a biologically active variant of the light-oxygen-voltage (LOV) domain protein VVD devoid of the adduct-forming Cys residue (VVD-III) to probe the mechanism of flavin photoreduction and protonation. A series of isosteric and conservative residue replacements studied by rate measurements, fluorescence quantum yields, FTIR difference spectroscopy, and molecular dynamics simulations indicate that tyrosine residues facilitate charge recombination reactions that limit sustained flavin reduction, whereas methionine residues facilitate radical propagation and quenching and also gate solvent access for flavin protonation. Replacement of a single surface Met residue with Leu favors formation of the ASQ over the NSQ and desensitizes photoreduction to oxidants. In contrast, increasing site hydrophilicity by Gln substitution promotes rapid NSQ formation and weakens the influence of the redox environment. Overall, the photoreactivity of VVD-III can be understood in terms of redundant electron donors, internal hole quenching, and coupled proton transfer reactions that all depend upon protein conformation, dynamics, and solvent penetration.

Nickel-Sulfonate Mode of Substrate Binding for Forward and Reverse Reactions of Methyl-SCoM Reductase Suggest a Radical Mechanism Involving Long-Range Electron Transfer

Methyl-coenzyme M reductase (MCR) catalyzes both the synthesis and the anaerobic oxidation of methane (AOM). Its catalytic site contains Ni at the core of cofactor F430. The Ni ion, in its low-valent Ni(I) state, lights the fuse leading to homolysis of the C-S bond of methyl-coenzyme M (methyl-SCoM) to generate a methyl radical, which abstracts a hydrogen atom from coenzyme B (HSCoB) to generate methane and the mixed disulfide CoMSSCoB. Direct reversal of this reaction activates methane to initiate anaerobic methane oxidation. On the basis of the crystal structures, which reveal a Ni-thiol interaction between Ni(II)-MCR and inhibitor CoMSH, a Ni(I)-thioether complex with substrate methyl-SCoM has been transposed to canonical MCR mechanisms. Similarly, a Ni(I)-disulfide with CoMSSCoB is proposed for the reverse reaction. However, this Ni(I)-sulfur interaction poses a conundrum for the proposed hydrogen-atom abstraction reaction because the >6 Å distance between the thiol group of SCoB and the thiol of SCoM observed in the structures appears to be too long for such a reaction. The spectroscopic, kinetic, structural, and computational studies described here establish that both methyl-SCoM and CoMSSCoB bind to the active Ni(I) state of MCR through their sulfonate groups, forming a hexacoordinate Ni(I)-N/O complex, not Ni(I)-S. These studies rule out direct Ni(I)-sulfur interactions in both substrate-bound states. As a solution to the mechanistic conundrum, we propose that both the forward and the reverse MCR reactions emanate through long-range electron transfer from the Ni(I)-sulfonate complexes with methyl-SCoM and CoMSSCoB, respectively.

Correlation between Charge Transport and Base Excision Repair in the MutY-DNA Glycosylase

Experimental evidence suggests that DNA-mediated redox signaling between high-potential [Fe₄S₄] proteins is relevant to DNA replication and repair processes, and protein-mediated charge transfer (CT) between [Fe₄S₄] clusters and nucleic acids is a fundamental process of the signaling and repair mechanisms. We analyzed the dominant CT pathways in the base excision repair glycosylase MutY using molecular dynamics simulations and hole hopping pathway analysis. We find that the adenine nucleobase of the mismatched A·oxoG DNA base pair facilitates [Fe₄S₄]-DNA CT prior to adenine excision by MutY. We also find that the R153L mutation in MutY (linked to colorectal adenomatous polyposis) influences the dominant [Fe4S4]-DNA CT pathways and appreciably decreases their effective CT rates.

Mutation effects on charge transport through the p58c iron-sulfur protein

Growing experimental evidence indicates that iron–sulfur proteins play key roles in DNA repair and replication. In particular, charge transport between [Fe₄S₄] clusters, mediated by proteins and DNA, may convey signals to coordinate enzyme action. Human primase is a well studied [Fe₄S₄] protein, and its p58c domain (which contains an [Fe₄S₄] cluster) plays a role in the initiation of DNA replication. The Y345C mutation in p58c is linked to gastric tumors and may influence the protein-mediated charge transport. The complexity of protein–DNA systems, and the intricate electronic structure of [Fe₄S₄] clusters, have impeded progress into understanding functional charge transport in these systems. In this study, we built force fields to describe the high potential [Fe₄S₄] cluster in both oxidation states. The parameterization is compatible with AMBER force fields and enabled well-balanced molecular dynamics simulations of the p58c–RNA/DNA complex relevant to the initiation of DNA replication. Using the molecular mechanics Poisson–Boltzmann and surface area solvation method on the molecular dynamics trajectories, we find that the p58c mutation induces a modest change in the p58c–duplex binding free energy in agreement with recent experiments. Through kinetic modeling and analysis, we identify key features of the main charge transport pathways in p58c. In particular, we find that the Y345C mutation partially changes the composition and frequency of the most efficient (and potentially relevant to the biological function) charge transport pathways between the [Fe₄S₄] cluster and the duplex. Moreover, our approach sets the stage for a deeper understanding of functional charge transfer in [Fe₄S₄] protein–DNA complexes.

Functional electron transfer between the [Fe₄S₄] cluster and the nucleic acid is impacted by a Y345C mutation in the p58c subunit of human primase.

Formation of a Copper(II)-Tyrosyl Complex at the Active Site of Lytic Polysaccharide Monooxygenases Following Oxidation by H2O2

Hydrogen peroxide is a cosubstrate for the oxidative cleavage of saccharidic substrates by copper-containing lytic polysaccharide monooxygenases (LPMOs). The rate of reaction of LPMOs with hydrogen peroxide is high, but it is accompanied by rapid inactivation of the enzymes, presumably through protein oxidation. Herein, we use UV-vis, CD, XAS, EPR, VT/VH-MCD, and resonance Raman spectroscopies, augmented with mass spectrometry and DFT calculations, to show that the product of reaction of an AA9 LPMO with H2O2 at higher pHs is a singlet Cu(II)-tyrosyl radical species, which is inactive for the oxidation of saccharidic substrates. The Cu(II)-tyrosyl radical center entails the formation of significant Cu(II)-(●OTyr) overlap, which in turn requires that the plane of the d(x2-y2) SOMO of the Cu(II) is orientated toward the tyrosyl radical. We propose from the Marcus cross-relation that the active site tyrosine is part of a "hole-hopping" charge-transfer mechanism formed of a pathway of conserved tyrosine and tryptophan residues, which can protect the protein active site from inactivation during uncoupled turnover.

Discovery, activity and characterisation of an AA10 lytic polysaccharide oxygenase from the shipworm symbiont Teredinibacter turnerae

BACKGROUND

The quest for novel enzymes for cellulosic biomass-degradation has recently been focussed on lytic polysaccharide monooxygenases (LPMOs/PMOs), Cu-containing proteins that catalyse the oxidative degradation of otherwise recalcitrant polysaccharides using O₂ or H₂O₂ as a co-substrate.

RESULTS

Although classical saprotrophic fungi and bacteria have been a rich source of lytic polysaccharide monooxygenases (LPMOs), we were interested to see if LPMOs from less evident bio-environments could be discovered and assessed for their cellulolytic activity in a biofuel context. In this regard, the marine shipworm Lyrodus pedicellatus represents an interesting source of new enzymes, since it must digest wood particles ingested during its natural tunnel boring behaviour and plays host to a symbiotic bacterium, Teredinibacter turnerae, the genome of which has revealed a multitude of enzymes dedicated to biomass deconstruction. Here, we show that T. turnerae encodes a cellulose-active AA10 LPMO. The 3D structure, at 1.4 Å resolution, along with its EPR spectrum is distinct from other AA10 polysaccharide monooxygenases insofar as it displays a "histidine-brace" catalytic apparatus with changes to the surrounding coordination sphere of the copper. Furthermore, TtAA10A possesses a second, surface accessible, Cu site 14 Å from the classical catalytic centre. Activity measurements show that the LPMO oxidises cellulose and thereby significantly augments the rate of degradation of cellulosic biomass by classical glycoside hydrolases.

CONCLUSION

Shipworms are wood-boring marine molluscs that can live on a diet of lignocellulose. Bacterial symbionts of shipworms provide many of the enzymes needed for wood digestion. The shipworm symbiont T. turnerae produces one of the few LPMOs yet described from the marine environment, notably adding to the capability of shipworms to digest recalcitrant polysaccharides.