Converged Structural and Spectroscopic Properties for Refined QM/MM Models of Azurin

HYSCORE and QM/MM Studies of Second Sphere Variants of the Type 1 Copper Site in Azurin: Influence of Mutations on the Hyperfine Couplings of Remote Nitrogens

Secondary coordination sphere (SCS) interactions have been shown to play important roles in tuning reduction potentials and electron transfer (ET) properties of the Type 1 copper proteins, but the precise roles of these interactions are not fully understood. In this work, we examined the influence of F114P, F114N, and N47S mutations in the SCS on the electronic structure of the T1 copper center in azurin (Az) by studying the hyperfine couplings of (i) histidine remote Nε nitrogens and (ii) the amide Np using the two-dimensional (2D) pulsed electron paramagnetic resonance (EPR) technique HYSCORE (hyperfine sublevel correlation) combined with quantum mechanics/molecular mechanics (QM/MM) and DLPNO-CCSD calculations. Our data show that some components of hyperfine tensor and isotropic coupling in N47SAz and F114PAz (but not F114NAz) deviate by up to ∼±20% from WTAz, indicating that these mutations significantly influence the spin density distribution between the CuII site and coordinating ligands. Furthermore, our calculations support the assignment of Np to the backbone amide of residue 47 (both in Asn and Ser variants). Since the spin density distributions play an important role in tuning the covalency of the Cu-Scys bond of Type 1 copper center that has been shown to be crucial in controlling the reduction potentials, this study provides additional insights into the electron spin factor in tuning the reduction potentials and ET properties.

Primary and Secondary Coordination Sphere Effects on the Structure and Function of S-Nitrosylating Azurin

Much progress has been made in understanding the roles of the secondary coordination sphere (SCS) in tuning redox potentials of metalloproteins. In contrast, the impact of SCS on reactivity is much less understood. A primary example is how copper proteins can promote S-nitrosylation (SNO), which is one of the most important dynamic post-translational modifications, and is crucial in regulating nitric oxide storage and transportation. Specifically, the factors that instill CuII with S-nitrosylating capabilities and modulate activity are not well understood. To address this issue, we investigated the influence of the primary and secondary coordination sphere on CuII-catalyzed S-nitrosylation by developing a series of azurin variants with varying catalytic capabilities. We have employed a multidimensional approach involving electronic absorption, S and Cu K-edge XAS, EPR, and resonance Raman spectroscopies together with QM/MM computational analysis to examine the relationships between structure and molecular mechanism in this reaction. Our findings have revealed that kinetic competency is correlated with three balancing factors, namely Cu-S bond strength, Cu spin localization, and relative S(ps) vs S(pp) contributions to the ground state. Together, these results support a reaction pathway that proceeds through the attack of the Cu-S bond rather than electrophilic addition to CuII or radical attack of SCys. The insights gained from this work provide not only a deeper understanding of SNO in biology but also a basis for designing artificial and tunable SNO enzymes to regulate NO and prevent diseases due to SNO dysregulation.

Dissecting the Molecular Origin of g-Tensor Heterogeneity and Strain in Nitroxide Radicals in Water: Electron Paramagnetic Resonance Experiment versus Theory

Nitroxides are common EPR sensors of microenvironmental properties such as polarity, numbers of H-bonds, pH, and so forth. Their solvation in an aqueous environment is facilitated by their high propensity to form H-bonds with the surrounding water molecules. Their g- and A-tensor elements are key parameters to extracting the properties of their microenvironment. In particular, the gxx value of nitroxides is rich in information. It is known to be characterized by discrete values representing nitroxide populations previously assigned to have different H-bonds with the surrounding waters. Additionally, there is a large g-strain, that is, a broadening of g-values associated with it, which is generally correlated with environmental and structural micro-heterogeneities. The g-strain is responsible for the frequency dependence of the apparent line width of the EPR spectra, which becomes evident at high field/frequency. Here, we address the molecular origin of the gxx heterogeneity and of the g-strain of a nitroxide moiety (HMI: 2,2,3,4,5,5-hexamethylimidazolidin-1-oxyl, C9H19N2O) in water. To treat the solvation effect on the g-strain, we combined a multi-frequency experimental approach with ab initio molecular dynamics simulations for structural sampling and quantum chemical EPR property calculations at the highest realistically affordable level, including an explicitly micro-solvated HMI ensemble and the embedded cluster reference interaction site model. We could clearly identify the distinct populations of the H-bonded nitroxides responsible for the gxx heterogeneity experimentally observed, and we dissected the role of the solvation shell, H-bond formation, and structural deformation of the nitroxide in the creation of the g-strain associated with each nitroxide subensemble. Two contributions to the g-strain were identified in this study. The first contribution depends on the number of hydrogen bonds formed between the nitroxide and the solvent because this has a large and well-understood effect on the gxx-shift. This contribution can only be resolved at high resonance frequencies, where it leads to distinct peaks in the gxx region. The second contribution arises from configurational fluctuations of the nitroxide that necessarily lead to g-shift heterogeneity. These contributions cannot be resolved experimentally as distinct resonances but add to the line broadening. They can be quantitatively analyzed by studying the apparent line width as a function of microwave frequency. Interestingly, both theory and experiment confirm that this contribution is independent of the number of H-bonds. Perhaps even more surprisingly, the theoretical analysis suggests that the configurational fluctuation broadening is not induced by the solvent but is inherently present even in the gas phase. Moreover, the calculations predict that this broadening decreases upon solvation of the nitroxide.

Hybrid Quantum Mechanical/Molecular Mechanical Methods For Studying Energy Transduction in Biomolecular Machines

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods have become indispensable tools for the study of biomolecules. In this article, we briefly review the basic methodological details of QM/MM approaches and discuss their applications to various energy transduction problems in biomolecular machines, such as long-range proton transports, fast electron transfers, and mechanochemical coupling. We highlight the particular importance for these applications of balancing computational efficiency and accuracy. Using several recent examples, we illustrate the value and limitations of QM/MM methodologies for both ground and excited states, as well as strategies for calibrating them in specific applications. We conclude with brief comments on several areas that can benefit from further efforts to make QM/MM analyses more quantitative and applicable to increasingly complex biological problems.

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.

Factors That Determine the Variation of Equilibrium and Kinetic Properties of QM/MM Enzyme Simulations: QM Region, Conformation, and Boundary Condition

To analyze the impact of various technical details on the results of quantum mechanical (QM)/molecular mechanical (MM) enzyme simulations, including the QM region size, catechol-O-methyltransferase (COMT) is studied as a model system using an approximate QM/MM method (DFTB3/CHARMM). The results show that key equilibrium and kinetic properties for methyl transfer in COMT exhibit limited variations with respect to the size of the QM region, which ranges from ∼100 to ∼500 atoms in this study. With extensive sampling, local and global structural characteristics of the enzyme are largely conserved across the studied QM regions, while the nature of the transition state (e.g., secondary kinetic isotope effect) and reaction exergonicity are largely maintained. Deviations in the free energy profile with different QM region sizes are similar in magnitude to those observed with changes in other simulation protocols, such as different initial enzyme conformations and boundary conditions. Electronic structural properties, such as the covariance matrix of residual charge fluctuations, appear to exhibit rather long-range correlations, especially when the peptide backbone is included in the QM region; this observation holds when a range-separated DFT approach is used as the QM region, suggesting that delocalization error is unlikely the origin. Overall, the analyses suggest that multiple simulation details determine the results of QM/MM enzyme simulations with comparable contributions.

Harder, better, faster, stronger: Large-scale QM and QM/MM for predictive modeling in enzymes and proteins

Computational prediction of enzyme mechanism and protein function requires accurate physics-based models and suitable sampling. We discuss recent advances in large-scale quantum mechanical (QM) modeling of biochemical systems that have reduced the cost of high-accuracy models. Tradeoffs between sampling and accuracy have motivated modeling with molecular mechanics (MM) in a multiscale QM/MM or iterative approach. Limitations to both conventional density-functional theory and classical MM force fields remain for describing noncovalent interactions in comparison to experiment or wavefunction theory. Because predictions of enzyme action (i.e. electrostatics), free energy barriers, and mechanisms are sensitive to the protocol and embedding method in QM/MM, convergence tests and systematic methods for quantifying QM-level interactions are a needed, active area of development.