Mechanistic insights on the reduction of glutathione disulfide by protein disulfide isomerase

Abstract

We explore the enzymatic mechanism of the reduction of glutathione disulfide (GSSG) by the reduced a domain of human protein disulfide isomerase (hPDI) with atomistic resolution. We use classical molecular dynamics and hybrid quantum mechanics/molecular mechanics calculations at the mPW1N/6-311+G(2d,2p):FF99SB//mPW1N/6-31G(d): FF99SB level. The reaction proceeds in two stages: (i) a thiol-disulfide exchange through nucleophilic attack of the Cys53-thiolate to the GSSG-disulfide followed by the deprotonation of Cys56-thiol by Glu47-carboxylate and (ii) a second thiol-disulfide exchange between the Cys56-thiolate and the mixed disulfide intermediate formed in the first step. The Gibbs activation energy for the first stage was 18.7 kcal·mol, and for the second stage, it was 7.2 kcal·mol ). Our results also suggest that the catalysis by protein disulfide isomerase (PDI) and thiol-disulfide exchange is mostly enthalpy-driven (entropy changes below 2 kcal·mol −1 at all stages of the reaction). Hydrogen bonds formed between the backbone of His55 and Cys56 and the Cys56-thiol result in an increase in the Gibbs energy barrier of the first thiol-disulfide exchange. The solvent plays a key role in stabilizing the leaving glutathione thiolate formed. This role is not exclusively electrostatic, because an explicit inclusion of several water molecules at the density-functional theory level is a requisite to form the mixed disulfide intermediate. In the intramolecular oxidation of PDI, a transition state is only observed if hydrogen bond donors are nearby the mixed disulfide intermediate, which emphasizes that the thermochemistry of thiol-disulfide exchange in PDI is influenced by the presence of hydrogen bond donors.PDI | GSH/GSSG buffer | thiol-disulfide exchange | protein folding | ONIOM

The Catalytic Mechanism of Human Transketolase

“…Only vibrational temperatures larger than 100 K have been accounted for, during the calculation of the entropy contribution, according to a validated procedure, [49] already successfully used in other works. [50][51][52] The final energy profiles were obtained adding the Zero Point Energy, the empirical D3-dispersion correction [47] and entropy contributions, as reported in Table S1.…”

Human Fatty Acid Synthase: A Computational Study of the Transfer of the Acyl Moieties from MAT to the ACP Domain

“…Vibrational frequency calculations at 25 °C were performed for each optimized structure to validate their nature as minima (absence of imaginary frequencies) or TS (single imaginary frequency). Only vibrational temperatures larger than 120 K (≈100 cm −1 ) were considered for the calculation of enthalpic and entropic corrections, as validated before . Single‐point (SP) QM/MM energy calculations were performed for all stationary points at the B3LYP/6‐311+G(2d,2p)‐D3:AMBER level.…”

Protocol for Computational Enzymatic Reactivity Based on Geometry Optimisation

“…In those cases, an extensive sampling of the conformational space is of central importance. However, if a good starting X‐ray geometry of the enzyme is available, our experience fortunately shows that dynamics are not crucial to unravel the correct enzymatic mechanistic pathway . Several methods were developed, which can tackle this.…”