Reaction of Cysteine and Selenocysteine with a Model Au(I)-N-Heterocyclic Carbene Complex: A QM and QM/MM MD Comparative Study
- 1Department of Chemistry, Faculty of Natural Sciences, Masaryk University, Kamenice 5, Brno 62500, Czech Republic
- 2Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, Prague 2 121 16, Czech Republic
Abstract
The substitution reaction of a model Au(I)-N-heterocyclic carbene (NHC) complex with cysteine and selenocysteine represents a key elementary step in the inhibition mechanism of thioredoxin reductase by gold-based metallodrugs. We present a systematic quantum-mechanical (QM) and combined QM/MM investigation of this reaction with the aim of assessing the consistency and reliability of activation and reaction free energies obtained at different theoretical levels. The reaction was first studied at the QM level using several density functional theory (DFT) functionals and basis sets and subsequently within a QM/MM framework, employing multiple QM descriptions of the reactive region and explicit molecular-mechanics (MM) treatment of the solvent (water). Comparison of QM and QM/MM results reveals how the inclusion of the surrounding environment affects both activation barriers and reaction energies. Several protonation states of cysteine and selenocysteine were considered, and a pH-dependent model was constructed using a grand-canonical transformation of the free energy. Combining QM/MM reaction energetics with pK a values yielded pH-dependent Gibbs–Alberty free energy profiles with a coordination-preference switch at pH ≈ 8. In acidic and neutral solutions the Au(I) complex coordinates with selenocysteine, while in a basic environment, coordination with cysteine is energetically more stable.
Related
- Synthesis, characterization and biological evaluation of newly designed silver(I)-N-heterocyclic carbene (NHC) complexes
- Synthesis, Cytotoxicity, and Mechanistic Evaluation of Thiazolium-Based Gold(I) Bis‑N‑Heterocyclic Carbene Complexes
- Toward Chemoselective Modification of Selenocysteine-Containing Biomolecules Via Doubly Activated Cyclopropanes
- Comparative analysis of adenosine deaminase (ADA) and polymerase chain reaction (PCR/GeneXpert) for diagnosing pleural and peritoneal tuberculosis
- An auto-generated ADP-ribosyl-cysteine degron drives DTX2-dependent turnover of PARP7
BibTeX
@article{riha2026reaction,
title = {{Reaction of Cysteine and Selenocysteine with a Model Au(I)-N-Heterocyclic Carbene Complex: A QM and QM/MM MD Comparative Study}},
author = {Milan Říha and Markéta Munzarová and Jaroslav V. Burda},
journal = {Inorganic Chemistry},
year = {2026},
month = oct,
doi = {10.1021/acs.inorgchem.6c04543},
url = {https://doi.org/10.1021/acs.inorgchem.6c04543},
note = {CC BY 4.0},
}
RIS
TY - JOUR
TI - Reaction of Cysteine and Selenocysteine with a Model Au(I)-N-Heterocyclic Carbene Complex: A QM and QM/MM MD Comparative Study
AU - Milan Říha
AU - Markéta Munzarová
AU - Jaroslav V. Burda
PY - 2026
DA - 2026/10/15
JO - Inorganic Chemistry
DO - 10.1021/acs.inorgchem.6c04543
UR - https://doi.org/10.1021/acs.inorgchem.6c04543
AB - The substitution reaction of a model Au(I)-N-heterocyclic carbene (NHC) complex with cysteine and selenocysteine represents a key elementary step in the inhibition mechanism of thioredoxin reductase by gold-based metallodrugs. We present a systematic quantum-mechanical (QM) and combined QM/MM investigation of this reaction with the aim of assessing the consistency and reliability of activation and reaction free energies obtained at different theoretical levels. The reaction was first studied at the QM level using several density functional theory (DFT) functionals and basis sets and subsequently within a QM/MM framework, employing multiple QM descriptions of the reactive region and explicit molecular-mechanics (MM) treatment of the solvent (water). Comparison of QM and QM/MM results reveals how the inclusion of the surrounding environment affects both activation barriers and reaction energies. Several protonation states of cysteine and selenocysteine were considered, and a pH-dependent model was constructed using a grand-canonical transformation of the free energy. Combining QM/MM reaction energetics with pK a values yielded pH-dependent Gibbs–Alberty free energy profiles with a coordination-preference switch at pH ≈ 8. In acidic and neutral solutions the Au(I) complex coordinates with selenocysteine, while in a basic environment, coordination with cysteine is energetically more stable.
SN - 0020-1669
SN - 1520-510X
ER -