Multiscale Modelling of Lytic Polysaccharide Monooxygenases

Impact of the Copper Second Coordination Sphere on Catalytic Performance and Substrate Specificity of a Bacterial Lytic Polysaccharide Monooxygenase

Lytic polysaccharide monooxygenases (LPMOs) catalyze the oxidative cleavage of glycosidic bonds in recalcitrant polysaccharides, such as cellulose and chitin, using a single copper cofactor bound in a conserved histidine brace with a more variable second coordination sphere. Cellulose-active LPMOs in the fungal AA9 family and in a subset of bacterial AA10 enzymes contain a His-Gln-Tyr second sphere motif, whereas other cellulose-active AA10s have an Arg-Glu-Phe motif. To shine a light on the impact of this variation, we generated single, double, and triple mutations changing the His216-Gln219-Tyr221 motif in cellulose- and chitin-oxidizing MaAA10B toward Arg-Glu-Phe. These mutations generally reduced enzyme performance due to rapid inactivation under turnover conditions, showing that catalytic fine-tuning of the histidine brace is complex and that the roles of these second sphere residues are strongly interconnected. Studies of copper reactivity showed remarkable effects, such as an increase in oxidase activity following the Q219E mutation and a strong dependence of this effect on the presence of Tyr at position 221. In reductant-driven reactions, differences in oxidase activity, which lead to different levels of in situ generated H2O2, correlated with differences in polysaccharide-degrading ability. The single Q219E mutant displayed a marked increase in activity on chitin in both reductant-driven reactions and reactions fueled by exogenously added H2O2. Thus, it seems that the evolution of substrate specificity in LPMOs involves both the extended substrate-binding surface and the second coordination sphere.

A frontier-orbital view of the initial steps of lytic polysaccharide monooxygenase reactions

Lytic polysaccharide monooxygenases (LPMOs) are copper enzymes that oxidatively cleave the strong C-H bonds in recalcitrant polysaccharide substrates, thereby playing a crucial role in biomass degradation. Recently, LPMOs have also been shown to be important for several pathogens. It is well established that the Cu(II) resting state of LPMOs is inactive, and the electronic structure of the active site needs to be altered to transform the enzyme into an active form. Whether this transformation occurs due to substrate binding or due to a unique priming reduction has remained speculative. Starting from four different crystal structures of the LPMO LsAA9A with well-defined oxidation states, we use a frontier molecular orbital approach to elucidate the initial steps of the LPMO reaction. We give an explanation for the requirement of the unique priming reduction and analyse electronic structure changes upon substrate binding. We further investigate how the presence of the substrate could facilitate an electron transfer from the copper active site to an H2O2 co-substrate. Our findings could help to control experimental LPMO reactions.

Understanding the initial events of the oxidative damage and protection mechanisms of the AA9 lytic polysaccharide monooxygenase family

Lytic polysaccharide monooxygenase (LPMO) is a new class of oxidoreductases that boosts polysaccharide degradation employing a copper active site. This boost may facilitate the cost-efficient production of biofuels and high-value chemicals from polysaccharides such as lignocellulose. Unfortunately, self-oxidation of the active site inactivates LPMOs. Other oxidoreductases employ hole-hopping mechanisms as protection against oxidative damage, but little is generally known about the details of these mechanisms. Herein, we employ highly accurate theoretical models based on density functional theory (DFT) molecular mechanics (MM) hybrids to understand the initial steps in LPMOs' protective measures against self-oxidation; we identify several intermediates recently proposed from experiment, and quantify which are important for protective hole-hopping pathways. Investigations on two different LPMOs show consistently that a tyrosine residue close to copper is crucial for protection: this explains recent experiments, showing that LPMOs without this tyrosine are more susceptible to self-oxidation.

Kβ X-ray Emission Spectroscopy of Cu(I)-Lytic Polysaccharide Monooxygenase: Direct Observation of the Frontier Molecular Orbital for H2O2 Activation

Lytic polysaccharide monooxygenases (LPMOs) catalyze the degradation of recalcitrant carbohydrate polysaccharide substrates. These enzymes are characterized by a mononuclear Cu(I) active site with a three-coordinate T-shaped "His-brace" configuration including the N-terminal histidine and its amine group as ligands. This study explicitly investigates the electronic structure of the d10 Cu(I) active site in a LPMO using Kβ X-ray emission spectroscopy (XES). The lack of inversion symmetry in the His-brace site enables the 3d/p mixing required for intensity in the Kβ valence-to-core (VtC) XES spectrum of Cu(I)-LPMO. These Kβ XES data are correlated to density functional theory (DFT) calculations to define the bonding, and in particular, the frontier molecular orbital (FMO) of the Cu(I) site. These experimentally validated DFT calculations are used to evaluate the reaction coordinate for homolytic cleavage of the H2O2 O-O bond and understand the contribution of this FMO to the low barrier of this reaction and how the geometric and electronic structure of the Cu(I)-LPMO site is activated for rapid reactivity with H2O2.

Molecular Mechanism of Substrate Oxidation in Lytic Polysaccharide Monooxygenases: Insight from Theoretical Investigations

Lytic polysaccharide monooxygenases (LPMOs) are copper enzymes that today comprise a large enzyme superfamily, grouped into the distinct members AA9-AA17 (with AA12 exempted). The LPMOs have the potential to facilitate the upcycling of biomass waste products by boosting the breakdown of cellulose and other recalcitrant polysaccharides. The cellulose biopolymer is the main component of biomass waste and thus comprises a large, unexploited resource. The LPMOs work through a catalytic, oxidative reaction whose mechanism is still controversial. For instance, the nature of the intermediate performing the oxidative reaction is an open question, and the same holds for the employed co-substrate. Here we review theoretical investigations addressing these questions. The applied theoretical methods are usually based on quantum mechanics (QM), often combined with molecular mechanics (QM/MM). We discuss advantages and disadvantages of the employed theoretical methods and comment on the interplay between theoretical and experimental results.

Distal Hydrophobic Loop Modulates the Copper Active Site and Reaction of AA13 Polysaccharide Monooxygenases

Polysaccharide monooxygenases (PMOs) use a type-2 copper center to activate O₂ for the selective hydroxylation of one of the two C-H bonds of glycosidic linkages. Our electron paramagnetic resonance (EPR) analysis and molecular dynamics (MD) simulations suggest the unprecedented dynamic roles of the loop containing the residue G89 (G89 loop) on the active site structure and reaction cycle of starch-active PMOs (AA13 PMOs). In the Cu(II) state, the G89 loop could switch between an "open" and "closed" conformation, which is associated with the binding and dissociation of an aqueous ligand in the distal site, respectively. The conformation of the G89 loop influences the positioning of the copper center on the preferred substrate of AA13 PMOs. The dissociation of the distal ligand results in the bending of the T-shaped core of the Cu(II) active site, which could help facilitate its reduction to the active Cu(I) state. In the Cu(I) state, the G89 loop is in the "closed" conformation with a confined copper center, which could allow for efficient O₂ binding. In addition, the G89 loop remains in the "closed" conformation in the Cu(II)-superoxo intermediate, which could prevent off-pathway superoxide release via exchange with the distal aqueous ligand. Finally, at the end of the reaction cycle, aqueous ligand binding to the distal site could switch the G89 loop to the "open" conformation and facilitate product release.

The crystal structure of CbpD clarifies substrate-specificity motifs in chitin-active lytic polysaccharide monooxygenases

The 3 Å resolution crystal structure of the Pseudomonas aeruginosa virulence factor CbpD both supports and challenges the current model of how lytic polysaccharide monooxygenases bind chitin and raises interesting possibilities about how type 2 secretion-system substrates may interact with the secretion machinery. This structure also demonstrates the utility of new, AI-powered, protein structure-prediction algorithms in making challenging structural targets tractable.

Pseudomonas aeruginosa secretes diverse proteins via its type 2 secretion system, including a 39 kDa chitin-binding protein, CbpD. CbpD has recently been shown to be a lytic polysaccharide monooxygenase active on chitin and to contribute substantially to virulence. To date, no structure of this virulence factor has been reported. Its first two domains are homologous to those found in the crystal structure of Vibrio cholerae GbpA, while the third domain is homologous to the NMR structure of the CBM73 domain of Cellvibrio japonicusCjLPMO10A. Here, the 3.0 Å resolution crystal structure of CbpD solved by molecular replacement is reported, which required ab initio models of each CbpD domain generated by the artificial intelligence deep-learning structure-prediction algorithm RoseTTAFold. The structure of CbpD confirms some previously reported substrate-specificity motifs among LPMOAA10s, while challenging the predictive power of others. Additionally, the structure of CbpD shows that post-translational modifications occur on the chitin-binding surface. Moreover, the structure raises interesting possibilities about how type 2 secretion-system substrates may interact with the secretion machinery and demonstrates the utility of new artificial intelligence protein structure-prediction algorithms in making challenging structural targets tractable.

Preliminary results of neutron and X-ray diffraction data collection on a lytic polysaccharide monooxygenase under reduced and acidic conditions

Lytic polysaccharide monooxygenases (LPMOs) are copper-center enzymes that are involved in the oxidative cleavage of the glycosidic bond in crystalline cellulose and other polysaccharides. The LPMO reaction is initiated by the addition of a reductant and oxygen to ultimately form an unknown activated copper-oxygen species that is responsible for polysaccharide-substrate H-atom abstraction. Given the sensitivity of metalloproteins to radiation damage, neutron protein crystallography provides a nondestructive technique for structural characterization while also informing on the positions of H atoms. Neutron cryo-crystallography permits the trapping of catalytic intermediates, thereby providing insight into the protonation states and chemical nature of otherwise short-lived species in the reaction mechanism. To characterize the reaction-mechanism intermediates of LPMO9D from Neurospora crassa, a cryo-neutron diffraction data set was collected from an ascorbate-reduced crystal. A second neutron diffraction data set was collected at room temperature from an LPMO9D crystal exposed to low-pH conditions to probe the protonation states of ionizable groups involved in catalysis under acidic conditions.

Estimating the accuracy of calculated electron paramagnetic resonance hyperfine couplings for a lytic polysaccharide monooxygenase

Lytic polysaccharide monooxygenases (LPMOs) are enzymes that bind polysaccharides followed by an (oxidative) disruption of the polysaccharide surface, thereby boosting depolymerization. The binding process between the LPMO catalytic domain and polysaccharide is key to the mechanism and establishing structure-function relationships for this binding is therefore crucial. The hyperfine coupling constants (HFCs) from EPR spectroscopy have proven useful for this purpose. Unfortunately, EPR does not provide direct structural data and therefore the experimental EPR parameters have to be supported with parameters calculated with density functional theory. Yet, calculated HFCs are extremely sensitive to the employed computational setup. Using the LPMO Ls(AA9)A catalytic domain, we here quantify the importance of several choices in the computational setup, ranging from the use of specialized basis, the underlying structures, and the employed exchange-correlation functional. We show that specialized basis sets are an absolute necessity, and also that care has to be taken in the optimization of the underlying structure: only by allowing large parts of the protein around the active site to structurally relax could we obtain results that uniformly reproduced experimental trends. We compare our results to previously published X-ray structures and experimental HFCs for Ls(AA9)A as well as to recent experimental/theoretical results for another (AA10) family of LPMOs.

The role of the active site tyrosine in the mechanism of lytic polysaccharide monooxygenase

Catalytic breakdown of polysaccharides can be achieved more efficiently by means of the enzymes lytic polysaccharide monooxygenases (LPMOs). However, the LPMO mechanism has remained controversial, preventing full exploitation of their potential. One of the controversies has centered around an active site tyrosine, present in most LPMO classes. Recent investigations have for the first time obtained direct (spectroscopic) evidence for the possibility of chemical modification of this tyrosine. However, the spectroscopic features obtained in the different investigations are remarkably different, with absorption maximum at 420 and 490 nm, respectively. In this paper we use density functional theory (DFT) in a QM/MM formulation to reconcile these (apparently) conflicting results. By modeling the spectroscopy as well as the underlying reaction mechanism we can show how formation of two isomers (both involving deprotonation of tyrosine) explains the difference in the observed spectroscopic features. Both isomers have a [TyrO–Cu–OH]⁺ moiety with the OH in either the cis- or trans-position to a deprotonated tyrosine. Although the cis-[TyrO–Cu–OH]⁺ moiety is well positioned for oxidation of the substrate, preliminary calculations with the substrate reveal that the reactivity is at best moderate, making a protective role of tyrosine more likely.

With QM/MM, we investigate the mechanism of tyrosine deprotonation in lytic polysaccharide monooxygenases. Our results support deprotonation and our calculated UV-vis spectra show that two isomers must be formed to match recent experiments.

Is density functional theory accurate for lytic polysaccharide monooxygenase enzymes?

The lytic polysaccharide monooxygenase (LPMO) enzymes boost polysaccharide depolymerization through oxidative chemistry, which has fueled the hope for more energy-efficient production of biofuel. We have recently proposed a mechanism for the oxidation of the polysaccharide substrate (E. D. Hedegård and U. Ryde, Chem. Sci., 2018, 9, 3866-3880). In this mechanism, intermediates with superoxide, oxyl, as well as hydroxyl (i.e. [CuO2]+, [CuO]+ and [CuOH]2+) cores were involved. These complexes can have both singlet and triplet spin states, and both spin-states may be important for how LPMOs function during catalytic turnover. Previous calculations on LPMOs have exclusively been based on density functional theory (DFT). However, different DFT functionals are known to display large differences for spin-state splittings in transition-metal complexes, and this has also been an issue for LPMOs. In this paper, we study the accuracy of DFT for spin-state splittings in superoxide, oxyl, and hydroxyl intermediates involved in LPMO turnover. As reference we employ multiconfigurational perturbation theory (CASPT2).

Dimethylanilinic N-Oxides and Their Oxygen Surrogacy Role in the Formation of a Putative High-Valent Copper-Oxygen Species

The reaction of p-cyano-N,N-dimethylaniline N-oxide, an O-atom donor, with different copper(I) complexes (at room temperature and in acetone) indicates the formation via O-atom transfer of a high-valent copper oxyl species, CuII-O•, a putative key intermediate in the catalytic cycle of copper-containing monooxygenases. The formation of p-cyano-N-hydroxymethyl-N-methylaniline and p-cyano-N-methylaniline as the main products of the reaction highlight the capability of this species to hydroxylate strong C-H bonds (bond dissociation energy ∼ 90 kcal/mol). A plausible mechanism for the reactivity of this catalytic system is proposed.

Molecular Mechanisms of Oxygen Activation and Hydrogen Peroxide Formation in Lytic Polysaccharide Monooxygenases

Lytic polysaccharide monooxygenases (LPMOs) are copper-dependent enzymes for the degradation of recalcitrant polysaccharides such as chitin and cellulose. Unlike classical hydrolytic enzymes (cellulases), LPMOs catalyze the cleavage of the glycosidic bond via an oxidative mechanism using oxygen and a reductant. The full enzymatic molecular mechanisms, starting from the initial electron transfer from a reductant to oxygen activation and hydrogen peroxide formation, are not yet understood. Using quantum mechanics/molecular mechanics (QM/MM) metadynamics simulations, we have uncovered the oxygen activation mechanisms by LPMO in the presence of ascorbic acid, one of the most-used reductants in LPMOs assays. Our simulations capture the sequential formation of Cu(II)-O₂ ^- and Cu(II)-OOH^- intermediates via facile H atom abstraction from ascorbate. By investigating all the possible reaction pathways from the Cu(II)-OOH^- intermediate, we ruled out Cu(II)-O^• - formation via direct O-O cleavage of Cu(II)-OOH^-. Meanwhile, we identified a possible pathway in which the proximal O atom of Cu(II)-OOH^- abstracts a hydrogen atom from ascorbate, leading to Cu(I) and H₂O₂. The in-situ-generated H₂O₂ either converts to LPMO-Cu(II)-O^• - via a homolytic reaction, or diffuses into the bulk water in an uncoupled pathway. The competition of these two pathways is strongly dependent on the binding of the carbohydrate substrate, which plays a role in barricading the in-situ-generated H₂O₂ molecule, preventing its diffusion from the active site into the bulk water. Based on the present results, we propose a catalytic cycle of LPMOs that is consistent with the experimental information available. In particular, it explains the enigmatic substrate dependence of the reactivity of the LPMO with H₂O₂.

Oxidoreductases and Reactive Oxygen Species in Conversion of Lignocellulosic Biomass

Biomass constitutes an appealing alternative to fossil resources for the production of materials and energy. The abundance and attractiveness of vegetal biomass come along with challenges pertaining to the intricacy of its structure, evolved during billions of years to face and resist abiotic and biotic attacks. To achieve the daunting goal of plant cell wall decomposition, microorganisms have developed many (enzymatic) strategies, from which we seek inspiration to develop biotechnological processes. A major breakthrough in the field has been the discovery of enzymes today known as lytic polysaccharide monooxygenases (LPMOs), which, by catalyzing the oxidative cleavage of recalcitrant polysaccharides, allow canonical hydrolytic enzymes to depolymerize the biomass more efficiently. Very recently, it has been shown that LPMOs are not classical monooxygenases in that they can also use hydrogen peroxide (H2O2) as an oxidant. This discovery calls for a revision of our understanding of how lignocellulolytic enzymes are connected since H2O2 is produced and used by several of them. The first part of this review is dedicated to the LPMO paradigm, describing knowns, unknowns, and uncertainties. We then present different lignocellulolytic redox systems, enzymatic or not, that depend on fluxes of reactive oxygen species (ROS). Based on an assessment of these putatively interconnected systems, we suggest that fine-tuning of H2O2 levels and proximity between sites of H2O2 production and consumption are important for fungal biomass conversion. In the last part of this review, we discuss how our evolving understanding of redox processes involved in biomass depolymerization may translate into industrial applications.

Reactivity of O2 versus H2O2 with polysaccharide monooxygenases

Enzymatic conversion of polysaccharides into lower-molecular-weight, soluble oligosaccharides is dependent on the action of hydrolytic and oxidative enzymes. Polysaccharide monooxygenases (PMOs) use an oxidative mechanism to break the glycosidic bond of polymeric carbohydrates, thereby disrupting the crystalline packing and creating new chain ends for hydrolases to depolymerize and degrade recalcitrant polysaccharides. PMOs contain a mononuclear Cu(II) center that is directly involved in C-H bond hydroxylation. Molecular oxygen was the accepted cosubstrate utilized by this family of enzymes until a recent report indicated reactivity was dependent on H2O2 Reported here is a detailed analysis of PMO reactivity with H2O2 and O2, in conjunction with high-resolution MS measurements. The cosubstrate utilized by the enzyme is dependent on the assay conditions. PMOs will directly reduce O2 in the coupled hydroxylation of substrate (monooxygenase activity) and will also utilize H2O2 (peroxygenase activity) produced from the uncoupled reduction of O2 Both cosubstrates require Cu reduction to Cu(I), but the reaction with H2O2 leads to nonspecific oxidation of the polysaccharide that is consistent with the generation of a hydroxyl radical-based mechanism in Fenton-like chemistry, while the O2 reaction leads to regioselective substrate oxidation using an enzyme-bound Cu/O2 reactive intermediate. Moreover, H2O2 does not influence the ability of secretome from Neurospora crassa to degrade Avicel, providing evidence that molecular oxygen is a physiologically relevant cosubstrate for PMOs.

A Random-Sequential Kinetic Mechanism for Polysaccharide Monooxygenases

Polysaccharide monooxygenases (PMOs) are mononuclear copper enzymes that catalyze the hydroxylation of polysaccharides leading to the scission of the glycosidic bond. The mechanism, in which PMOs utilize molecular oxygen to oxidize the polysaccharide substrate, still remains largely unknown. Here, steady-state kinetics assays were used to probe the mechanism of oxygen-dependent cellohexaose oxidation catalyzed by MtPMO9E. Kinetic analysis indicated that both k_cat/ K_M(O2) and k_cat/ K_M(Glc6) were dependent on the concentration of the second substrate. Inhibition studies using carbon monoxide were also carried out. In addition, K_D values for Glc6 were determined for the Cu(I) and Cu(II) forms of the enzyme. Taken together, PMOs follow a random-sequential kinetic mechanism to form a ternary ES-O₂ complex. The optimal pH for MtPMO9E turnover was determined to be between pH 6.00 and pH 7.00. Furthermore, the kinetic parameters k_cat, k_cat/ K_M(O2), and k_cat/ K_M(Glc6) demonstrate a decrease in PMO activity at a low pH and provide equivalent kinetic p K_a's of 5.10. This points to the protonation of a general base required for turnover. These results provide a basis for the initial chemical steps in the mechanism of PMOs.

Molecular mechanism of lytic polysaccharide monooxygenases

The lytic polysaccharide monooxygenases (LPMOs) are copper metalloenzymes that can enhance polysaccharide depolymerization through an oxidative mechanism and hence boost generation of biofuel from e.g. cellulose. By employing density functional theory in a combination of quantum mechanics and molecular mechanics (QM/MM), we report a complete description of the molecular mechanism of LPMOs.

The lytic polysaccharide monooxygenases (LPMOs) are copper metalloenzymes that can enhance polysaccharide depolymerization through an oxidative mechanism and hence boost generation of biofuel from e.g. cellulose. By employing density functional theory in a combination of quantum mechanics and molecular mechanics (QM/MM), we report a complete description of the molecular mechanism of LPMOs. The QM/MM scheme allows us to describe all reaction steps with a detailed protein environment and we show that this is necessary. Several active species capable of abstracting a hydrogen from the substrate have been proposed previously and starting from recent crystallographic work on a substrate–LPMO complex, we investigate previously suggested paths as well as new ones. We describe the generation of the reactive intermediates, the abstraction of a hydrogen atom from the polysaccharide substrate, as well as the final recombination step in which OH is transferred back to the substrate. We show that a superoxo [CuO₂]⁺ complex can be protonated by a nearby histidine residue (suggested by recent mutagenesis studies and crystallographic work) and, provided an electron source is available, leads to formation of an oxyl-complex after cleavage of the O–O bond and dissociation of water. The oxyl complex either reacts with the substrate or is further protonated to a hydroxyl complex. Both the oxyl and hydroxyl complexes are also readily generated from a reaction with H₂O₂, which was recently suggested to be the true co-substrate, rather than O₂. The C–H abstraction by the oxyl and hydroxy complexes is overall favorable with activation barriers of 69 and 94 kJ mol^–1, compared to the much higher barrier (156 kJ mol^–1) obtained for the copper–superoxo species. We obtain good structural agreement for intermediates for which structural data are available and the estimated reaction energies agree with experimental rate constants. Thus, our suggested mechanism is the most complete to date and concur with available experimental evidence.

Oxygen Activation by Cu LPMOs in Recalcitrant Carbohydrate Polysaccharide Conversion to Monomer Sugars

Natural carbohydrate polymers such as starch, cellulose, and chitin provide renewable alternatives to fossil fuels as a source for fuels and materials. As such, there is considerable interest in their conversion for industrial purposes, which is evidenced by the established and emerging markets for products derived from these natural polymers. In many cases, this is achieved via industrial processes that use enzymes to break down carbohydrates to monomer sugars. One of the major challenges facing large-scale industrial applications utilizing natural carbohydrate polymers is rooted in the fact that naturally occurring forms of starch, cellulose, and chitin can have tightly packed organizations of polymer chains with low hydration levels, giving rise to crystalline structures that are highly recalcitrant to enzymatic degradation. The topic of this review is oxidative cleavage of carbohydrate polymers by lytic polysaccharide mono-oxygenases (LPMOs). LPMOs are copper-dependent enzymes (EC 1.14.99.53-56) that, with glycoside hydrolases, participate in the degradation of recalcitrant carbohydrate polymers. Their activity and structural underpinnings provide insights into biological mechanisms of polysaccharide degradation.

Targeting the reactive intermediate in polysaccharide monooxygenases

Lytic polysaccharide monooxygenases (LPMOs) are copper metalloenzymes that can enhance polysaccharide depolymerization through an oxidative mechanism, making them interesting for the production of biofuel from cellulose. However, the details of this activation are unknown; in particular, the nature of the intermediate that attacks the glycoside C-H bond in the polysaccharide is not known, and a number of different species have been suggested. The homolytic bond-dissociation energy (BDE) has often been used as a descriptor for the bond-activation power, especially for inorganic model complexes. We have employed quantum-chemical cluster calculations to estimate the BDE for a number of possible LPMO intermediates to bridge the gap between model complexes and the actual LPMO active site. The calculated BDEs suggest that the reactive intermediate is either a Cu(II)-oxyl, a Cu(III)-oxyl, or a Cu(III)-hydroxide, which indicate that O-O bond breaking occurs before the C-H activation step.