Multiscale Modeling of Lytic Polysaccharide Monooxygenases in Oxidized and Reduced Redox States
- 1Tunneling Group, Biotechnology Centre, Silesian University of Technology, Gliwice 44-100, Poland
- 2Molecular Simulations and Design Group, Max Planck Institute for Dynamics of Complex Technical Systems, Magdeburg 39106, Germany
Abstract
Lytic polysaccharide monooxygenases (LPMOs) are copper-dependent oxidative enzymes whose catalytic activity depends on redox cycling of a single copper ion in the active site, coordinated by the conserved histidine brace motif. Although the static coordination geometry of the active site has been characterized crystallographically, the dynamic response of the histidine brace to Cu(I)/Cu(II) cycling and the role of second-sphere residues in modulating this response remain poorly understood at the atomistic level. Here, we address this gap using a multiscale computational approach combining molecular dynamics (MD) simulations, QM/MM optimization, and EPR parameter calculations, applied to two representative LPMOs: NcAA9C (AA9, axial tyrosine) and BaAA10A (AA10, axial phenylalanine). The structural response to the change in oxidation state is seen almost entirely in the exogenous ligands rather than in the histidine brace itself, with the oxidized form showing markedly higher and longer-lived water coordination. In BaAA10A, the equatorial water is anchored by a second-sphere glutamate with no counterpart in NcAA9C. QM/MM electrostatic perturbation analysis identifies Asp76 (NcAA9C) and Glu68 (BaAA10A) as the dominant stabilizing residues of the catalytic core, acting through water-mediated hydrogen-bond bridges that are more pronounced in the Cu(II) state. Computed EPR parameters respond to the displacement of the copper site from a planar geometry, and the equatorial water position was shown to have the dominant effect, differentiating the two LPMOs.
Related
- Multiscale electro-thermal modeling of the Branly effect in oxidized copper granular media
- How Metal-to-Linker Ratio and Concentration Govern the Nucleation Mechanism of MOF-5-inspired Crystals: Insights from Multiscale Modeling
- From Materials to Devices: Multiscale Modeling of Colloidal Quantum Dot Optoelectronics
- Quantum Multiscale Modeling: A Hierarchy of Algorithms for Complex Chemical Systems
- Advances in Multiscale Computational Modeling of 2D Nanomaterials for Solar Cells: From Interface Design to Fabrication and Stability
BibTeX
@article{szlepershendy2026multiscale,
title = {{Multiscale Modeling of Lytic Polysaccharide Monooxygenases in Oxidized and Reduced Redox States}},
author = {Katarzyna Szleper-Shendy and Andrés M. Escorcia and Ravi Kumar and Artur Góra and Matthias Stein},
journal = {The Journal of Physical Chemistry. B},
year = {2026},
month = oct,
doi = {10.1021/acs.jpcb.6c04364},
url = {https://doi.org/10.1021/acs.jpcb.6c04364},
note = {CC BY 4.0},
}
RIS
TY - JOUR
TI - Multiscale Modeling of Lytic Polysaccharide Monooxygenases in Oxidized and Reduced Redox States
AU - Katarzyna Szleper-Shendy
AU - Andrés M. Escorcia
AU - Ravi Kumar
AU - Artur Góra
AU - Matthias Stein
PY - 2026
DA - 2026/10/15
JO - The Journal of Physical Chemistry. B
DO - 10.1021/acs.jpcb.6c04364
UR - https://doi.org/10.1021/acs.jpcb.6c04364
AB - Lytic polysaccharide monooxygenases (LPMOs) are copper-dependent oxidative enzymes whose catalytic activity depends on redox cycling of a single copper ion in the active site, coordinated by the conserved histidine brace motif. Although the static coordination geometry of the active site has been characterized crystallographically, the dynamic response of the histidine brace to Cu(I)/Cu(II) cycling and the role of second-sphere residues in modulating this response remain poorly understood at the atomistic level. Here, we address this gap using a multiscale computational approach combining molecular dynamics (MD) simulations, QM/MM optimization, and EPR parameter calculations, applied to two representative LPMOs: NcAA9C (AA9, axial tyrosine) and BaAA10A (AA10, axial phenylalanine). The structural response to the change in oxidation state is seen almost entirely in the exogenous ligands rather than in the histidine brace itself, with the oxidized form showing markedly higher and longer-lived water coordination. In BaAA10A, the equatorial water is anchored by a second-sphere glutamate with no counterpart in NcAA9C. QM/MM electrostatic perturbation analysis identifies Asp76 (NcAA9C) and Glu68 (BaAA10A) as the dominant stabilizing residues of the catalytic core, acting through water-mediated hydrogen-bond bridges that are more pronounced in the Cu(II) state. Computed EPR parameters respond to the displacement of the copper site from a planar geometry, and the equatorial water position was shown to have the dominant effect, differentiating the two LPMOs.
SN - 1520-6106
SN - 1520-5207
ER -