Search

Hydrogen peroxide is a natural oxidant that can oxidize protein thiols (RSH) via sulfenic acid (RSOH) and sulfinic acid (RSO2H) to sulfonic acid (RSO3H). In this paper, we study the complete anionic and neutral oxidation pathway from thiol to sulfonic acid. Reaction barriers and reaction free energies for all three oxidation steps are computed, both for the isolated substrates and for the substrates in the presence of different model ligands (CH4, H2O, NH3) mimicking the enzymatic environment. We found for all three barriers that the anionic thiolate is more reactive than the neutral thiol. However, the assistance of the environment in the neutral pathway in a solvent-assisted proton-exchange (SAPE) mechanism can lower the reaction barrier noticeably. Polar ligands can decrease the reaction barriers, whereas apolar ligands do not influence the barrier heights. The same holds for the reaction energies: they decrease (become more negative) in the presence of polar ligands whereas apolar ligands do not have an influence. The consistently negative consecutive reaction energies for the oxidation in the anionic pathway when going from thiolate over sulfenic and sulfinic acid to sulfonic acid are in agreement with biological reversibility.

Peroxiredoxins catalyze the reduction of peroxides, a process of vital importance to survive oxidative stress. A nucleophilic cysteine, also known as the peroxidatic cysteine, is responsible for this catalytic process. We used the Mycobacterium tuberculosis alkyl hydroperoxide reductase E (MtAhpE) as a model to investigate the effect of the chemical environment on the specificity of the reaction. Using an integrative structural (R116A – PDB 4XIH; F37H – PDB 5C04), kinetic and computational approach, we explain the mutational effects of key residues in its environment. This study shows that the active site residues are specifically oriented to create an environment which selectively favours a reaction with peroxides.

In many established methods, identification of hydrogen bonds (H-bonds) is primarily based on pairwise comparison of distances between atoms. These methods often give rise to systematic errors when sulfur is involved. A more accurate method is the non-covalent interaction index, which determines the strength of the H-bonds based on the associated electron density and its gradient. We applied the NCI index on the active site of a single-cysteine peroxiredoxin. We found a different sulfur hydrogen-bonding network to that typically found by established methods and we propose a more accurate equation for determining sulfur H-bonds based on geometrical criteria. This new algorithm will be implemented in the next release of the widely-used CHARMM program (version 41b) and will be particularly useful for analyzing water molecule-mediated H-bonds involving different atom types. Furthermore, based on the identification of the weakest sulfur-water H-bond, the location of hydrogen peroxide for the nucleophilic attack by the cysteine sulfur can be predicted. In general, current methods to determine H-bonds will need to be reevaluated, thereby leading to better understanding of the catalytic mechanisms in which sulfur chemistry is involved.

Recently, it was found that mutations in the binding cavity of drug-metabolizing Cytochrome P450 BM3 mutants can result in major changes in regioselectivity in testosterone (TES) hydroxylation. In the current work, we report the intrinsic reactivity of TES' C-H bonds and our attempts to rationalize experimentally observed changes in TES hydroxylation using a protein structure-based in silico approach, by setting up and employing a combined Molecular Dynamics (MD) and ligand docking approach to account for the flexibility and plasticity of BM3 mutants. Using this approach, about 100,000 TES binding poses were obtained per mutant. The predicted regioselectivity in TES hydroxylation by the mutants was found to be in disagreement with experiment. As revealed in a detailed structural analysis of the obtained docking poses, this disagreement is due to limitations in correctly scoring hydrogen-bonding and steric interactions with specific active-site residues, which could explain the experimentally observed trends in regioselectivity in TES hydroxylation. © 2012 Bentham Science Publishers.

Inversion of stereoselectivity: screening of a minimal mutant library revealed a cytochrome P450 BM3 variant M01 A82W S72I capable of producing 16 α-OH-testosterone. Remarkably, a single active site mutation S72I in M01 A82W inverted the stereoselectivity of hydroxylation from 16 β to 16 α. Introduction of S72I mutation in another 16 β-OH-selective variant M11 V87I, also resulted in similar inversion of stereoselectivity.

The conjugation of reactive drug metabolites to GSH is considered an important detoxification mechanism that can be spontaneous and/or mediated by glutathione S-transferases (GSTs). In case GSTs play an important role in GSH conjugation, genetically determined deficiencies in GSTs may be a risk factor for adverse drug reactions (ADRs) resulting from reactive drug metabolites. So far, the role of GSTs in the detoxification of reactive intermediates of clozapine, a drug-causing idiosyncratic drug reactions (IDRs), has not been studied. In the present study, we studied the ability of four recombinant human GSTs (hGST A1-1, hGST M1-1, hGST P1-1, and hGST T1-1) to catalyze the GSH conjugation of reactive metabolites of clozapine, formed in vitro by human and rat liver microsomes and drug-metabolizing P450 BM3 mutant, P450 102A1M11H. Consistent with previous studies, in the absence of GSTs, three GSH conjugates were identified derived from the nitrenium ion of clozapine. In the presence of three of the GSTs, hGST P1-1, hGST M1-1, and hGST A1-1, total GSH conjugation was strongly increased in all bioactivation systems tested. The highest activity was observed with hGST P1-1, whereas hGST M1-1 and hGST A1-1 showed slightly lower activity. Polymorphic hGST T1-1 did not show any activity in catalyzing GSH conjugation of reactive clozapine metabolites. Interestingly, the addition of hGSTs resulted in major changes in the regioselectivity of GSH conjugation of the reactive clozapine metabolite, possibly due to the different active site geometries of hGSTs. Two GSH conjugates found were completely dependent on the presence of hGSTs. Chlorine substitution of the clozapine nitrenium ion, which so far was only observed in in vivo studies, appeared to be the major pathway of hGST P1-1-catalyzed GSH conjugation, whereas hGST A1-1 and hGST M1-1 also showed significant activity. The second GSH conjugate, previously also only found in in vivo studies, was also formed by hGST P1-1 and to a small extent by hGST A1-1. These results demonstrate that human GSTs may play a significant role in the inactivation of reactive intermediates of clozapine. Therefore, further studies are required to investigate whether genetic polymorphisms of hGST P1-1 and hGST M1-1 contribute to the interindividual differences in susceptibility to clozapine-induced adverse drug reactions. © 2010 American Chemical Society.

Protein thiol/sulfenic acid oxidation potentials provide a tool to select specific oxidation agents, but are experimentally difficult to obtain. Here, insights into the thiol sulfenylation thermodynamics are obtained from model calculations on small systems and from a quantum mechanics/molecular mechanics (QM/MM) analysis on human 2-Cys peroxiredoxin thioredoxin peroxidase B (Tpx-B). To study thiol sulfenylation in Tpx-B, our recently developed computational method to determine reduction potentials relatively compared to a reference system and based on reaction energies reduction potential from electronic energies is updated. Tpx-B forms a sulfenic acid (R-SO−) on one of its active site cysteines during reactive oxygen scavenging. The observed effect of the conserved active site residues is consistent with the observed hydrogen bond interactions in the QM/MM optimized Tpx-B structures and with free energy calculations on small model systems. The ligand effect could be linked to the complexation energies of ligand L with CH3S− and CH3SO−. Compared to QM only calculations on Tpx-B’s active site, the QM/MM calculations give an improved understanding of sulfenylation thermodynamics by showing that other residues from the protein environment other than the active site residues can play an important role.