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

Histidine oxidation in lytic polysaccharide monooxygenase

The lytic polysaccharide monooxygenases (LPMOs) comprise a super-family of copper enzymes that boost the depolymerisation of polysaccharides by oxidatively disrupting the glycosidic bonds connecting the sugar units. Industrial use of LPMOs for cellulose depolymerisation has already begun but is still far from reaching its full potential. One issue is that the LPMOs self-oxidise and thereby deactivate. The mechanism of this self-oxidation is unknown, but histidine residues coordinating to the copper atom are the most susceptible. An unusual methyl modification of the NE2 atom in one of the coordinating histidine residues has been proposed to have a protective role. Furthermore, substrate binding is also known to reduce oxidative damage. We here for the first time investigate the mechanism of histidine oxidation with combined quantum and molecular mechanical (QM/MM) calculations, with outset in intermediates previously shown to form from a reaction with peroxide and a reduced LPMO. We show that an intermediate with a [Cu-O]⁺ moiety is sufficiently potent to oxidise the nearest C-H bond on both histidine residues, but methylation of the NE2 atom of His-1 increases the reaction barrier of this reaction. The substrate further increases the activation barrier. We also investigate a [Cu-OH]²⁺ intermediate with a deprotonated tyrosine radical. This intermediate was previously proposed to have a protective role, and we also find it to have higher barriers than the corresponding a [Cu-O]⁺ intermediate.

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.

Challenges in Protein QM/MM Simulations with Intra-Backbone Link Atoms

Hybrid quantum mechanical/molecular mechanical (QM/MM) simulations fuel discoveries in many fields of science including computational biochemistry and enzymology. Development of more convenient tools leads to an increase in the number of works in which mechanical insights into enzymes' mode of operation are obtained. Most commonly, these tools feature hydrogen-capping (link atom) approach to provide coupling between QM and MM subsystems across a covalent bond. Extensive studies were conducted to provide a solid foundation for the correctness of such an approach when a bond to a nonpolar MM atom is considered. However, not every task may be accomplished this way. Certain scenarios of using QM/MM in computational enzymology encourage or even necessitate the incorporation of backbone atoms into the QM region. Two out of three backbone atoms are polar, and in QM/MM with electrostatic embedding, a neighboring link atom will be hyperpolarized. Several schemes to mitigate this effect were previously proposed alongside a rigorous assessment of quantitative effects on model systems. However, it was not clear whether they may translate into qualitatively different results and how link atom hyperpolarization may manifest itself in a real-life enzymological scenario. Here, we show that the consequences of such an artifact may be severe and may completely overturn the conclusions drawn from the simulations. Our case advocates for the use of charge redistribution schemes whenever intra-backbone QM/MM boundaries are considered. Moreover, we addressed how different boundary types and charge redistribution schemes influence backbone dynamics. We showed that the results are heavily dependent on which boundary MM terms are retained, with charge alteration being of secondary importance. In the worst case, only three intra-backbone boundaries may be used with relative confidence in the adequacy of resulting simulations, irrespective of the hyperpolarization mitigation scheme. Thus, advances in the field are certainly needed to fuel new discoveries. As of now, we believe that issues raised in this work might encourage authors in the field to report what boundaries, boundary MM terms, and charge redistribution schemes they are using, so their results may be correctly interpreted.

Changes in active-site geometry on X-ray photoreduction of a lytic polysaccharide monooxygenase active-site copper and saccharide binding

The recently discovered lytic polysaccharide monooxygenases (LPMOs) are Cu-containing enzymes capable of degrading polysaccharide substrates oxidatively. The generally accepted first step in the LPMO reaction is the reduction of the active-site metal ion from Cu²⁺ to Cu⁺. Here we have used a systematic diffraction data collection method to monitor structural changes in two AA9 LPMOs, one from Lentinus similis (LsAA9_A) and one from Thermoascus auranti-acus (TaAA9_A), as the active-site Cu is photoreduced in the X-ray beam. For LsAA9_A, the protein produced in two different recombinant systems was crystallized to probe the effect of post-translational modifications and different crystallization conditions on the active site and metal photoreduction. We can recommend that crystallographic studies of AA9 LPMOs wishing to address the Cu²⁺ form use a total X-ray dose below 3 × 10⁴ Gy, while the Cu⁺ form can be attained using 1 × 10⁶ Gy. In all cases, we observe the transition from a hexa-coordinated Cu site with two solvent-facing ligands to a T-shaped geometry with no exogenous ligands, and a clear increase of the θ₂ parameter and a decrease of the θ₃ parameter by averages of 9.2° and 8.4°, respectively, but also a slight increase in θ_T. Thus, the θ₂ and θ₃ parameters are helpful diagnostics for the oxidation state of the metal in a His-brace protein. On binding of cello-oligosaccharides to LsAA9_A, regardless of the production source, the θ_T parameter increases, making the Cu site less planar, while the active-site Tyr-Cu distance decreases reproducibly for the Cu²⁺ form. Thus, the θ_T increase found on copper reduction may bring LsAA9_A closer to an oligosaccharide-bound state and contribute to the observed higher affinity of reduced LsAA9_A for cellulosic substrates.

Refining details of the structural and electronic properties of the CuB site in pMMO enzyme through sequential molecular dynamics/CPKS-EPR calculations

This work investigated the structural and electronic properties of the copper mononuclear site of the PmoB part of the pMMO enzyme at the molecular level. We propose that the Cu_B catalytic site in the soluble portion of pMMO at room temperature and under physiological conditions is a mononuclear copper complex in a distorted octahedral arrangement with the residues His³³, His¹³⁷, and His¹³⁹ on the equatorial base and two water molecules on the axial axis. Our view was based on the molecular dynamics results and DFT calculations of the electronic paramagnetic resonance parameters and comparisons with experimental EPR data. This new proposed model for the Cu_B site brings additional support concerning the recent experimental evidence, which pointed out that a saturated coordination sphere of the copper ion in the Cu_B center is an essential factor that makes it less efficient than the Cu_C site in the methane oxidation. Therefore, according to the Cu_B site model proposed here, an additional step involving a displacement of at least one water molecule of the copper coordination sphere by the O₂ molecule prior to its activation must be necessary. This scenario is less likely to occur in the Cu_C center once this one is buried in the alpha-helices, which are part of the pMMO structure bound to the membrane wall, and consequently located in a less solvent-exposed region. In addition, we also present a simple and efficient sequential S-MD/CPKS protocol to compute EPR parameters that can, in principle, be expanded for the study of other copper-containing proteins.

Decoding the Ambiguous Electron Paramagnetic Resonance Signals in the Lytic Polysaccharide Monooxygenase from Photorhabdus luminescens

Understanding the structure and function of lytic polysaccharide monooxygenases (LPMOs), copper enzymes that degrade recalcitrant polysaccharides, requires the reliable atomistic interpretation of electron paramagnetic resonance (EPR) data on the Cu(II) active site. Among various LPMO families, the chitin-active PlAA10 shows an intriguing phenomenology with distinct EPR signals, a major rhombic and a minor axial signal. Here, we combine experimental and computational investigations to uncover the structural identity of these signals. X-band EPR spectra recorded at different pH values demonstrate pH-dependent population inversion: the major rhombic signal at pH 6.5 becomes minor at pH 8.5, where the axial signal dominates. This suggests that a protonation change is involved in the interconversion. Precise structural interpretations are pursued with quantum chemical calculations. Given that accurate calculations of Cu g-tensors remain challenging for quantum chemistry, we first address this problem via a thorough calibration study. This enables us to define a density functional that achieves accurate and reliable prediction of g-tensors, giving confidence in our evaluation of PlAA10 LPMO models. Large models were considered that include all parts of the protein matrix surrounding the Cu site, along with the characteristic second-sphere features of PlAA10. The results uniquely identify the rhombic signal with a five-coordinate Cu ion bearing two water molecules in addition to three N-donor ligands. The axial signal is attributed to a four-coordinate Cu ion where only one of the waters remains bound, as hydroxy. Alternatives that involve decoordination of the histidine brace amino group are unlikely based on energetics and spectroscopy. These results provide a reliable spectroscopy-consistent view on the plasticity of the resting state in PlAA10 LPMO as a foundation for further elucidating structure-property relationships and the formation of catalytically competent species. Our strategy is generally applicable to the study of EPR parameters of mononuclear copper-containing metalloenzymes.

EPR Spectroscopy of Cu(II) Complexes: Prediction of g-Tensors Using Double-Hybrid Density Functional Theory

Computational electron paramagnetic resonance (EPR) spectroscopy is an important field of applied quantum chemistry that contributes greatly to connecting spectroscopic observations with the fundamental description of electronic structure for open-shell molecules. However, not all EPR parameters can be predicted accurately and reliably for all chemical systems. Among transition metal ions, Cu(II) centers in inorganic chemistry and biology, and their associated EPR properties such as hyperfine coupling and g-tensors, pose exceptional difficulties for all levels of quantum chemistry. In the present work, we approach the problem of Cu(II) g-tensor calculations using double-hybrid density functional theory (DHDFT). Using a reference set of 18 structurally and spectroscopically characterized Cu(II) complexes, we evaluate a wide range of modern double-hybrid density functionals (DHDFs) that have not been applied previously to this problem. Our results suggest that the current generation of DHDFs consistently and systematically outperform other computational approaches. The B2GP-PLYP and PBE0-DH functionals are singled out as the best DHDFs on average for the prediction of Cu(II) g-tensors. The performance of the different functionals is discussed and suggestions are made for practical applications and future methodological developments.

Mapping the protonation states of the histidine brace in an AA10 lytic polysaccharide monooxygenase using CW-EPR spectroscopy and DFT calculations

The active site of the polysaccharide-degrading lytic polysaccharide monooxygenase (LPMO) enzyme features a single copper ion coordinated by a histidine brace. The primary coordination sphere of the copper contains several ligating atoms which are bonded to ionisable protons (e.g. OH₂, NH₂), the pK_as of which are unknown. Using a combination of CW-EPR X-band spectroscopy over a range of pH values and DFT calculations, we show that the active site of a chitin-active AA10 LPMO can exist in three different protonation states (pK_a1 = 8.7, pK_a2 ∼ 11.5), representing the ionisation of the coordinating groups. The middle pH species (fully formed at pH ∼ 10.5) is proposed to be Cu(II)(His)₂(OH)₂ (N₂O₂ coordination) with a decoordinated R-NH₃⁺ group at the amino terminus. This species also sees a rotation of the SOMO equatorial plane from the canonical histidine brace plane, whereby the nominal Cu d(x² - y²)-orbital has rotated some 45° along the His-Cu(II)-His axis, driven by the elongation and decoordination of the amino group. The highest pH species (>12) is proposed to exist as a Cu(II)-azanide, in which the NH₂ of the amino terminus has been deprotonated. The high pH means that this species is unlikely to be biologically relevant in the catalytic cycle of AA10 LPMOs.

Mechanism of DOPA radical generation and transfer in metal-free class Ie ribonucleotide reductase based on density functional theory

•

The mechanism of DOPA radical generation, transfer and regeneration is revealed.

•

The superoxide O₂^•− should be protonated to HO₂^• prior to oxidizing Tyr126 to DOPA radical.

•

The protonation of Asp88 is the prerequisite for the DOPA radical generation and radical transfer.

•

Lys213 is a key residue for the transfer of the DOPA radical.

Quantum mechanical/molecular mechanical (QM/MM) calculations were carried out to investigate the mechanisms of the generation, transfer, and regeneration of the DOPA radical for metal-free class Ie ribonucleotide reductase. The crystal structure of MfR2 (Nature, 2018, 563, 416–420) was adopted for the calculations. The QM/MM calculations have revealed several key points that are vital for understanding the mechanisms. The superoxide O₂^•− provided by the flavoprotein NrdI cannot directly oxidize the residue Tyr126 to the DOPA radical. It should be protonated to HO₂^•. The calculation results suggest that the covalent modification of Tyr126 and the DOPA radical generation can be carried out with no involvement of metal cofactors. This addresses the concerns of the articles (Nature, 2018, 563, 416–420; PNAS, 2018, 115, 10022–10027). Another concern from the articles is that how the DOPA radical is transferred from the radical trap. The DFT calculations have demonstrated that Lys213 is a key residue for the radical transfer from the DOPA radical. The ε-amino group of Lys213 is used not only as a bridge for the electron transfer but also as a proton donor. It can provide a proton to DOPA126 via a water molecule, and thus the radical transfer from DOPA126 to Trp52 is facilitated. It has also been revealed that the protonation of Asp88 is the prerequisite for the DOPA radical generation and the radical transfer in class Ie. Once the radical is quenched, it can be regenerated via the oxidations by superoxide O₂^•− and hydroperoxyl radical HO₂^•.