Increase of Redox Potential during the Evolution of Enzymes Degrading Recalcitrant Lignin

Fluctuation Relations to Calculate Protein Redox Potentials from Molecular Dynamics Simulations

The tunable design of protein redox potentials promises to open a range of applications in biotechnology and catalysis. Here, we introduce a method to calculate redox potential changes by combining fluctuation relations with molecular dynamics simulations. It involves the simulation of reduced and oxidized states, followed by the instantaneous conversion between them. Energy differences introduced by the perturbations are obtained using the Kubo-Onsager approach. Using a detailed fluctuation relation coupled with Bayesian inference, these are postprocessed into estimates for the redox potentials in an efficient manner. This new method, denoted MD + CB, is tested on a de novo four-helix bundle heme protein (the m4D2 "maquette") and five designed mutants, including some mutants characterized experimentally in this work. The MD + CB approach is found to perform reliably, giving redox potential shifts with reasonably good correlation (0.85) to the experimental values for the mutants. The MD + CB approach also compares well with redox potential shift predictions using a continuum electrostatic method. The estimation method employed within the MD + CB approach is straightforwardly transferable to standard equilibrium MD simulations and holds promise for redox protein engineering and design applications.

Recent Theoretical Insights into the Oxidative Degradation of Biopolymers and Plastics by Metalloenzymes

Molecular modeling techniques have become indispensable in many fields of molecular sciences in which the details related to mechanisms and reactivity need to be studied at an atomistic level. This review article provides a collection of computational modeling works on a topic of enormous interest and urgent relevance: the properties of metalloenzymes involved in the degradation and valorization of natural biopolymers and synthetic plastics on the basis of both circular biofuel production and bioremediation strategies. In particular, we will focus on lytic polysaccharide monooxygenase, laccases, and various heme peroxidases involved in the processing of polysaccharides, lignins, rubbers, and some synthetic polymers. Special attention will be dedicated to the interaction between these enzymes and their substrate studied at different levels of theory, starting from classical molecular docking and molecular dynamics techniques up to techniques based on quantum chemistry.

Ancestral sequence reconstruction as a tool to study the evolution of wood decaying fungi

The study of evolution is limited by the techniques available to do so. Aside from the use of the fossil record, molecular phylogenetics can provide a detailed characterization of evolutionary histories using genes, genomes and proteins. However, these tools provide scarce biochemical information of the organisms and systems of interest and are therefore very limited when they come to explain protein evolution. In the past decade, this limitation has been overcome by the development of ancestral sequence reconstruction (ASR) methods. ASR allows the subsequent resurrection in the laboratory of inferred proteins from now extinct organisms, becoming an outstanding tool to study enzyme evolution. Here we review the recent advances in ASR methods and their application to study fungal evolution, with special focus on wood-decay fungi as essential organisms in the global carbon cycling.

Lignin peroxidase-catalyzed direct oxidation of trace organic pollutants through a long-range electron transfer mechanism: Using propranolol as an example

In this work, lignin peroxidase (LiP) was extracted for the in vitro degradation of a persistent compound (propranolol, PPN). The results showed that 94.2% of PPN was degraded at 30 U L^-1 LiP activity and 10 mg L^-1 PPN. The PPN degradation rate increased from 33.5% to 94.2% when the veratryl alcohol (VA) concentration varied from 0 to 180 µM, but decreased to 73.1% with further VA addition. This phenomenon confirmed that VA was indispensable, however, it also acted as a competitive inhibitor of PPN oxidation. Computational analysis revealed that the Trp171…iron porphyrin (TRP-FeP) path was responsible for specific substrate (e.g., VA) transformation, and another long-range electron transfer (LRET) path through His-Asp…FeP (HSP-FeP) was discovered for non-specific substrate (e.g., PPN) degradation. These two electron-transfer routes shared one catalytic center, and VA protected the enzyme from H₂O₂-dependent inactivation. The HSP-FeP path transformed PPN through single electron transfer or H abstraction mechanisms. In addition, hydroxyl radicals generated in the LiP/H₂O₂ system were involved in the hydroxylation of the PPN intermediates. Possible degradation pathways were deduced using these degradation mechanisms and mass-spectrometry analysis. The multipath degradation mechanism endowed LiP with a remarkable capacity for removing various recalcitrant pollutants in environmental remediation.

Deep Insights into the Specific Evolution of Fungal Hybrid B Heme Peroxidases

Simple Summary

Fungi are well equipped to cope with oxidative stress and the reactive oxygen species that are, in the case of phytopathogens, produced mainly by the plant host for defence purposes. Peroxidases represent the major line of evolution for rapid decomposition of harmful peroxides in all aerobically metabolising organisms. In all the sequenced fungal genomes, many divergent genes coding for various peroxidases have been discovered, and Hybrid B heme peroxidases represent a distinctive mode of fungal-gene evolution within a large peroxidase–catalase superfamily that ranges from bacteria to plants.

Abstract

In this study, we focus on a detailed bioinformatics analysis of hyBpox genes, mainly within the genomes of Sclerotiniaceae (Ascomycota, Leotiomycetes), which is a specifically evolved fungal family of necrotrophic host generalists and saprophytic or biotrophic host specialists. Members of the genus Sclerotium produce only sclerotia and no fruiting bodies or spores. Thus, their physiological role for peroxidases remains open. A representative species, S. cepivorum, is a dangerous plant pathogen causing white rot in Allium species, particularly in onions, leeks, and garlic. On a worldwide basis, the white rot caused by this soil-borne fungus is apparently the most serious threat to Allium-crop production. We have also found very similar peroxidase sequences in the related fungus S. sclerotiorum, although with minor yet important modifications in the architecture of its active centre. The presence of ScephyBpox1-specific mRNA was confirmed by transcriptomic analysis. The presence of Hybrid B peroxidase at the protein level as the sole extracellular peroxidase of this fungus was confirmed in the secretome of S. cepivorum through detailed proteomic analyses. This prompted us to systematically search for all available genes coding for Hybrid B heme peroxidases in the whole fungal family of Sclerotiniaceae. We present here a reconstruction of their molecular phylogeny and analyse the unique aspects of their conserved-sequence features and structural folds in corresponding ancestral sequences.

Agaricales Mushroom Lignin Peroxidase: From Structure-Function to Degradative Capabilities

Lignin biodegradation has been extensively studied in white-rot fungi, which largely belong to order Polyporales. Among the enzymes that wood-rotting polypores secrete, lignin peroxidases (LiPs) have been labeled as the most efficient. Here, we characterize a similar enzyme (ApeLiP) from a fungus of the order Agaricales (with ~13,000 described species), the soil-inhabiting mushroom Agrocybe pediades. X-ray crystallography revealed that ApeLiP is structurally related to Polyporales LiPs, with a conserved heme-pocket and a solvent-exposed tryptophan. Its biochemical characterization shows that ApeLiP can oxidize both phenolic and non-phenolic lignin model-compounds, as well as different dyes. Moreover, using stopped-flow rapid spectrophotometry and 2D-NMR, we demonstrate that ApeLiP can also act on real lignin. Characterization of a variant lacking the above tryptophan residue shows that this is the oxidation site for lignin and other high redox-potential substrates, and also plays a role in phenolic substrate oxidation. The reduction potentials of the catalytic-cycle intermediates were estimated by stopped-flow in equilibrium reactions, showing similar activation by H₂O₂, but a lower potential for the rate-limiting step (compound-II reduction) compared to other LiPs. Unexpectedly, ApeLiP was stable from acidic to basic pH, a relevant feature for application considering its different optima for oxidation of phenolic and nonphenolic compounds.

Mutational and structural analysis of an ancestral fungal dye-decolorizing peroxidase

Dye-decolorizing peroxidases (DyPs) constitute a superfamily of heme-containing peroxidases that are related neither to animal nor to plant peroxidase families. These are divided into four classes (types A, B, C, and D) based on sequence features. The active site of DyPs contains two highly conserved distal ligands, an aspartate and an arginine, the roles of which are still controversial. These ligands have mainly been studied in class A-C bacterial DyPs, largely because no effective recombinant expression systems have been developed for the fungal (D-type) DyPs. In this work, we employ ancestral sequence reconstruction (ASR) to resurrect a D-type DyP ancestor, AncDyPD-b1. Expression of AncDyPD-b1 in Escherichia coli results in large amounts of a heme-containing soluble protein and allows for the first mutagenesis study on the two distal ligands of a fungal DyP. UV-Vis and resonance Raman (RR) spectroscopic analyses, in combination with steady-state kinetics and the crystal structure, reveal fine pH-dependent details about the heme active site structure and show that both the aspartate (D222) and the arginine (R390) are crucial for hydrogen peroxide reduction. Moreover, the data indicate that these two residues play important but mechanistically different roles on the intraprotein long-range electron transfer process. DATABASE: Structural data are available in the PDB database under the accession number 7ANV.

Experimental and theoretical insights into the effects of pH on catalysis of bond-cleavage by the lignin peroxidase isozyme H8 from Phanerochaete chrysosporium

BACKGROUND

Lignin peroxidases catalyze a variety of reactions, resulting in cleavage of both β-O-4' ether bonds and C-C bonds in lignin, both of which are essential for depolymerizing lignin into fragments amendable to biological or chemical upgrading to valuable products. Studies of the specificity of lignin peroxidases to catalyze these various reactions and the role reaction conditions such as pH play have been limited by the lack of assays that allow quantification of specific bond-breaking events. The subsequent theoretical understanding of the underlying mechanisms by which pH modulates the activity of lignin peroxidases remains nascent. Here, we report on combined experimental and theoretical studies of the effect of pH on the enzyme-catalyzed cleavage of β-O-4' ether bonds and of C-C bonds by a lignin peroxidase isozyme H8 from Phanerochaete chrysosporium and an acid stabilized variant of the same enzyme.

RESULTS

Using a nanostructure initiator mass spectrometry assay that provides quantification of bond breaking in a phenolic model lignin dimer we found that catalysis of degradation of the dimer to products by an acid-stabilized variant of lignin peroxidase isozyme H8 increased from 38.4% at pH 5 to 92.5% at pH 2.6. At pH 2.6, the observed product distribution resulted from 65.5% β-O-4' ether bond cleavage, 27.0% Cα-C1 carbon bond cleavage, and 3.6% Cα-oxidation as by-product. Using ab initio molecular dynamic simulations and climbing-image Nudge Elastic Band based transition state searches, we suggest the effect of lower pH is via protonation of aliphatic hydroxyl groups under which extremely acidic conditions resulted in lower energetic barriers for bond-cleavages, particularly β-O-4' bonds.

CONCLUSION

These coupled experimental results and theoretical explanations suggest pH is a key driving force for selective and efficient lignin peroxidase isozyme H8 catalyzed depolymerization of the phenolic lignin dimer and further suggest that engineering of lignin peroxidase isozyme H8 and other enzymes involved in lignin depolymerization should include targeting stability at low pH.

Setting the stage for evolution of a new enzyme

The evolution of novel enzymes has fueled the diversification of life on earth for billions of years. Insights into events that set the stage for the evolution of a new enzyme can be obtained from ancestral reconstruction and laboratory evolution. Ancestral reconstruction can reveal the emergence of a promiscuous activity in a pre-existing protein and the impact of subsequent mutations that enhance a new activity. Laboratory evolution provides a more holistic view by revealing mutations elsewhere in the genome that indirectly enhance the level of a newly important enzymatic activity. This review will highlight recent studies that probe the early stages of the evolution of a new enzyme from these complementary points of view.

Comparing Ligninolytic Capabilities of Bacterial and Fungal Dye-Decolorizing Peroxidases and Class-II Peroxidase-Catalases

We aim to clarify the ligninolytic capabilities of dye-decolorizing peroxidases (DyPs) from bacteria and fungi, compared to fungal lignin peroxidase (LiP) and versatile peroxidase (VP). With this purpose, DyPs from Amycolatopsis sp., Thermomonospora curvata, and Auricularia auricula-judae, VP from Pleurotus eryngii, and LiP from Phanerochaete chrysosporium were produced, and their kinetic constants and reduction potentials determined. Sharp differences were found in the oxidation of nonphenolic simple (veratryl alcohol, VA) and dimeric (veratrylglycerol-β- guaiacyl ether, VGE) lignin model compounds, with LiP showing the highest catalytic efficiencies (around 15 and 200 s⁻¹·mM⁻¹ for VGE and VA, respectively), while the efficiency of the A. auricula-judae DyP was 1–3 orders of magnitude lower, and no activity was detected with the bacterial DyPs. VP and LiP also showed the highest reduction potential (1.28–1.33 V) in the rate-limiting step of the catalytic cycle (i.e., compound-II reduction to resting enzyme), estimated by stopped-flow measurements at the equilibrium, while the T. curvata DyP showed the lowest value (1.23 V). We conclude that, when using realistic enzyme doses, only fungal LiP and VP, and in much lower extent fungal DyP, oxidize nonphenolic aromatics and, therefore, have the capability to act on the main moiety of the native lignin macromolecule.

Functional annotation and comparative modeling of ligninolytic enzymes from Trametes villosa (SW.) Kreisel for biotechnological applications

Functional annotation of Trametes villosa genome was performed to search Class II peroxidase proteins in this white-rot fungus, which can be valuable for several biotechnological processes. After sequence identification and manual curation, five proteins were selected to build 3 D models by comparative modeling. Analysis of sequential and structural sequences from selected targets revealed the presence of two putative Lignin Peroxidase and three putative Manganese Peroxidase on this fungal genome. All 3 D models had a similar folding pattern from selected 3 D structure templates. After minimization and validation steps, the best 3 D models were subjected to docking studies and molecular dynamics to identify structural requirements and the interactions required for molecular recognition. Two reliable 3 D models of Class II peroxidases, with typical catalytic site and architecture, and its protein sequences are indicated to recombinant production in biotechnological applications, such as bioenergy.Communicated by Ramaswamy H. Sarma.

Peroxidase evolution in white-rot fungi follows wood lignin evolution in plants

We analyze the evolution of ligninolytic peroxidases from wood-rotting fungi using conifer and angiosperm lignin as representatives of 2 steps of lignin evolution. By enzyme resurrection, we show that during fungal evolution, these enzymes improved their activity and switched their degradative preferences with the rise of a surface tryptophan conferring on them the ability to oxidize nonphenolic lignin. We calibrated the peroxidase phylogeny and determined that this residue appeared coincident with angiosperm diversification, characterized by the synthesis of a more complex and less phenolic lignin due to the general incorporation of a new unit in its structure. This way, we show that fungal evolution followed that of lignin synthesis, pointing to a coevolution between fungal saprotrophs and their plant hosts.

A comparison of sequenced Agaricomycotina genomes suggests that efficient degradation of wood lignin was associated with the appearance of secreted peroxidases with a solvent-exposed catalytic tryptophan. This hypothesis is experimentally demonstrated here by resurrecting ancestral fungal peroxidases, after sequence reconstruction from genomes of extant white-rot Polyporales, and evaluating their oxidative attack on the lignin polymer by state-of-the-art analytical techniques. Rapid stopped-flow estimation of the transient-state constants for the 2 successive one-electron transfers from lignin to the peroxide-activated enzyme (k_2app and k_3app) showed a progressive increase during peroxidase evolution (up to 50-fold higher values for the rate-limiting k_3app). The above agreed with 2-dimensional NMR analyses during steady-state treatments of hardwood lignin, showing that its degradation (estimated from the normalized aromatic signals of lignin units compared with a control) and syringyl-to-guaiacyl ratio increased with the enzyme evolutionary distance from the first peroxidase ancestor. More interestingly, the stopped-flow estimations of electron transfer rates also showed how the most recent peroxidase ancestors that already incorporated the exposed tryptophan into their molecular structure (as well as the extant lignin peroxidase) were comparatively more efficient at oxidizing hardwood (angiosperm) lignin, while the most ancestral “tryptophanless” enzymes were more efficient at abstracting electrons from softwood (conifer) lignin. A time calibration of the ancestry of Polyporales peroxidases localized the appearance of the first peroxidase with a solvent-exposed catalytic tryptophan to 194 ± 70 Mya, coincident with the diversification of angiosperm plants characterized by the appearance of dimethoxylated syringyl lignin units.

Different fungal peroxidases oxidize nitrophenols at a surface catalytic tryptophan

Dye-decolorizing peroxidase (DyP) from Auricularia auricula-judae and versatile peroxidase (VP) from Pleurotus eryngii oxidize the three mononitrophenol isomers. Both enzymes have been overexpressed in Escherichia coli and in vitro activated. Despite their very different three-dimensional structures, the nitrophenol oxidation site is located at a solvent-exposed aromatic residue in both DyP (Trp377) and VP (Trp164), as revealed by liquid chromatography coupled to mass spectrometry and kinetic analyses of nitrophenol oxidation by the native enzymes and their tryptophan-less variants (the latter showing 10-60 fold lower catalytic efficiencies).