Modeling CH <sub>3</sub> SOH–aromatic complexes to probe cysteine sulfenic acid–aromatic interactions in proteins

Cysteine sulfenic acid (CysOH), formed on oxidation of Cys residues, is an intermediate in the catalytic cycles of numerous antioxidant enzymes and participates in oxidative-stress sensing and redox signaling. Proteins control CysOH reactivity in part by its interactions with aromatic residues. To characterize such interactions, we performed extensive <i>ab initio</i> quantum mechanical calculations with MP2(full)/6-311++G(d,p) on complexes of CH<sub>3</sub>SOH as a CysOH model with side-chain models for Phe (toluene), Trp (3-methylindole), Tyr/Tyr<sup>-</sup> (4-methylphenol/4-methylphenolate) and His/HisH<sup>+</sup> (4-methylimidazole/4-methylimidazolium) residues. The gas-phase global minima conformers extracted from the 67 aromatic complexes found exhibit binding energies of ∼-5 to -25 kcal mol<sup>-1</sup>. In the neutral CH<sub>3</sub>SOH-aromatics, the center oxidized is the stronger H-bond donor, which varies with the geometry of the complex as does the ionization potential (IPV). While CH<sub>3</sub>SOH (IPV = 9.20 eV) is exclusively oxidized when complexed to 4-methylimidazolium (IPV = 14.64 eV), the phenol ring is oxidized in all CH<sub>3</sub>SOH complexes with 4-methylphenolate (IPV = 3.31 eV). To perform molecular dynamics (MD) simulations of the aqueous complexes, a potential model was optimized for CH<sub>3</sub>SOH and calibrated for its interactions with the aromatic ions. The MD simulations reveal that in bulk water the S atom preferentially adopts en-face or intermediate binding geometry with binding free energies of ∼-0.6, -2.5 and -5 kcal mol<sup>-1</sup> for the neutral, imidazolium and phenolate complexes, respectively. Overall, the gas-phase and aqueous CH<sub>3</sub>SOH complexes are 40-170% more stable and 0-40% less stable, respectively, than their CH<sub>3</sub>SH counterparts. Exceptionally, aqueous 4-methylphenolate binds CH<sub>3</sub>SOH ∼50% more tightly than CH<sub>3</sub>SH due to strong σ-type O-H⋯O<sub>ar</sub> H-bond bonding. Examination of a subset of CysOH-aromatics from the Protein Data Bank highlight their role in CysOH formation and stabilization in proteins.

Authors

Institutions

Publication Details

Journal
Physical Chemistry Chemical Physics
Published
2025-01-01
DOI
https://doi.org/10.1039/d5cp02976g
Citations
1
Primary Topic
Free Radicals and Antioxidants
Type
article
Field-Weighted Citation Impact
0.94
Controls
|||
ALL TIME
JAN
FEB
MAR
APR
MAY
JUN
JUL
AUG
SEP
OCT
article

Modeling CH 3 SOH–aromatic complexes to probe cysteine sulfenic acid–aromatic interactions in proteins

Esam A. Orabi
1 citations
Physical Chemistry Chemical Physics
Free Radicals and Antioxidants
0.94
article

Modeling CH 3 SOH–aromatic complexes to probe cysteine sulfenic acid–aromatic interactions in proteins

Esam A. Orabi
article en
1 citations

Abstract

Cysteine sulfenic acid (CysOH), formed on oxidation of Cys residues, is an intermediate in the catalytic cycles of numerous antioxidant enzymes and participates in oxidative-stress sensing and redox signaling. Proteins control CysOH reactivity in part by its interactions with aromatic residues. To characterize such interactions, we performed extensive <i>ab initio</i> quantum mechanical calculations with MP2(full)/6-311++G(d,p) on complexes of CH<sub>3</sub>SOH as a CysOH model with side-chain models for Phe (toluene), Trp (3-methylindole), Tyr/Tyr<sup>-</sup> (4-methylphenol/4-methylphenolate) and His/HisH<sup>+</sup> (4-methylimidazole/4-methylimidazolium) residues. The gas-phase global minima conformers extracted from the 67 aromatic complexes found exhibit binding energies of ∼-5 to -25 kcal mol<sup>-1</sup>. In the neutral CH<sub>3</sub>SOH-aromatics, the center oxidized is the stronger H-bond donor, which varies with the geometry of the complex as does the ionization potential (IPV). While CH<sub>3</sub>SOH (IPV = 9.20 eV) is exclusively oxidized when complexed to 4-methylimidazolium (IPV = 14.64 eV), the phenol ring is oxidized in all CH<sub>3</sub>SOH complexes with 4-methylphenolate (IPV = 3.31 eV). To perform molecular dynamics (MD) simulations of the aqueous complexes, a potential model was optimized for CH<sub>3</sub>SOH and calibrated for its interactions with the aromatic ions. The MD simulations reveal that in bulk water the S atom preferentially adopts en-face or intermediate binding geometry with binding free energies of ∼-0.6, -2.5 and -5 kcal mol<sup>-1</sup> for the neutral, imidazolium and phenolate complexes, respectively. Overall, the gas-phase and aqueous CH<sub>3</sub>SOH complexes are 40-170% more stable and 0-40% less stable, respectively, than their CH<sub>3</sub>SH counterparts. Exceptionally, aqueous 4-methylphenolate binds CH<sub>3</sub>SOH ∼50% more tightly than CH<sub>3</sub>SH due to strong σ-type O-H⋯O<sub>ar</sub> H-bond bonding. Examination of a subset of CysOH-aromatics from the Protein Data Bank highlight their role in CysOH formation and stabilization in proteins.

Physical Chemistry Chemical PhysicsVol. 27(39)
National Heart, Lung, and Blood Institute (US), PROTEO (CA), Concordia University (CA)
Openalex Percentile: Top 21%
Free Radicals and Antioxidants
0.94
AI Navigator

Ask Laika to Summarize, Analyze, and Connect papers live on the map.

Summarize Papers & Methodologies

Extract key findings, datasets, and comparative methods across publications.

Benchmark Rankings & Visual Analytics

Rank top research institutions, authors, funders, topics, and journals by Field-Weighted Citation Impact (FWCI) and paper volume with instant charts.

Connect Distant Disciplines

Bridge topological clusters on the map to find hidden collaborative intersections.