Dataset Information

Mechanistic implications of the cysteine-nicotinamide adduct in aldehyde dehydrogenase based on quantum mechanical/molecular mechanical simulations.

ABSTRACT: Recent computer simulations of the cysteine nucleophilic attack on propanal in human mitochondrial aldehyde dehydrogenase (ALDH2) yielded an unexpected result: the chemically reasonable formation of a dead-end cysteine-cofactor adduct when NAD+ was in the "hydride transfer" position. More recently, this adduct found independent crystallographic support in work on formyltetrahydrofolate dehydrogenase, work which further found evidence of the same adduct on re-examination of deposited electron densities of ALDH2. Although the experimental data showed that this adduct was reversible, several mechanistic questions arise from the fact that it forms at all. Here, we present results from further quantum mechanical/molecular mechanical (QM/MM) simulations toward understanding the mechanistic implications of adduct formation. These simulations revealed formation of the oxyanion thiohemiacetal intermediate only when the nicotinamide ring of NAD+ is oriented away from the active site, contrary to prior arguments. In contrast, and in seeming paradox, when NAD is oriented to receive the hydride, disassociation of the oxyanion intermediate to form the dead-end adduct is more thermodynamically favored than maintaining the oxyanion intermediate necessary for catalysis to proceed. However, this disassociation to the adduct could be avoided through proton transfer from the enzyme to the intermediate. Our results continue to indicate that the unlikely source of this proton is the Cys302 main chain amide.

SUBMITTER: Wymore T

PROVIDER: S-EPMC2529467 | biostudies-literature | 2007 Aug

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Mechanistic implications of the cysteine-nicotinamide adduct in aldehyde dehydrogenase based on quantum mechanical/molecular mechanical simulations.

Wymore Troy T Deerfield David W DW Hempel John J

Biochemistry 20070727 33

Recent computer simulations of the cysteine nucleophilic attack on propanal in human mitochondrial aldehyde dehydrogenase (ALDH2) yielded an unexpected result: the chemically reasonable formation of a dead-end cysteine-cofactor adduct when NAD+ was in the "hydride transfer" position. More recently, this adduct found independent crystallographic support in work on formyltetrahydrofolate dehydrogenase, work which further found evidence of the same adduct on re-examination of deposited electron den ...[more]

PMID: 17655326

Similar Datasets

Project description:Flavin adenine dinucleotide (FAD) is a common cofactor in redox proteins, and its reduction potentials are controlled by the protein environment. This regulation is mainly responsible for the versatile catalytic functions of flavoenzymes. In this article, we report computations of the reduction potentials of FAD in medium-chain acyl-CoA dehydrogenase (MCAD) and cholesterol oxidase (CHOX). In addition, the reduction potentials of lumiflavin in aqueous solution have also been computed. Using molecular dynamics and free-energy perturbation techniques, we obtained the free-energy changes for two-electron/two-proton as well as one-electron/one-proton addition steps. We employed a combined quantum mechanical and molecular mechanical (QM/MM) potential, in which the flavin ring was represented by the self-consistent-charge density functional tight-binding (SCC-DFTB) method, while the rest of the enzyme-solvent system was treated by classical force fields. The computed two-electron/two-proton reduction potentials for lumiflavin and the two enzyme-bound FADs are in reasonable agreement with experimental data. The calculations also yielded the pKa values for the one-electron reduced semiquinone (FH*) and the fully reduced hydroquinone (FH2) forms. The pKa of the FAD semiquinone in CHOX was found to be around 4, which is 4 units lower than that in the enzyme-free state and 2 units lower than that in MCAD; this supports the notion that oxidases have a greater ability than dehydrogenases to stabilize anionic semiquinones. In MCAD, the flavin ring interacts with four hydrophobic residues and has a significantly bent structure, even in the oxidized state. The present study shows that this bending of the flavin imparts a significant destabilization (approximately 5 kcal/mol) to the oxidized state. The reduction potential of lumiflavin was also computed using DFT (M06-L and B3LYP functionals with 6-31+G(d,p) basis set) with the SM6 continuum solvation model, and the results are in good agreement with results from explicit free-energy simulations, which supports the conclusion that the SCC-DFTB/MM computation is reasonably accurate for both 1e(-)/1H+ and 2e(-)/2H+ reduction processes. These results suggest that the first coupled electron-proton addition is stepwise for both the free and the two enzyme-bound flavins. In contrast, the second coupled electron-proton addition is also stepwise for the free flavin but is likely to be concerted when the flavin is bound to either the dehydrogenase or the oxidase enzyme.

Project description:DNA bases can adopt energetically unfavorable tautomeric forms that enable the formation of Watson-Crick-like (WC-like) mispairs, which have been proposed to give rise to spontaneous mutations in DNA and misincorporation errors in DNA replication and translation. Previous NMR and computational studies have indicated that the population of WC-like guanine-thymine (G-T) mispairs depends on the environment, such as the local nucleic acid sequence and solvation. To investigate these environmental effects, herein G-T mispair tautomerization processes are studied computationally in aqueous solution, in A-form and B-form DNA duplexes, and within the active site of a DNA polymerase λ variant. The wobble G-T (wG-T), WC-like G-T*, and WC-like G*-T forms are considered, where * indicates the enol tautomer of the base. The minimum free energy paths for the tautomerization from the wG-T to the WC-like G-T* and from the WC-like G-T* to the WC-like G*-T are computed with mixed quantum mechanical/molecular mechanical (QM/MM) free energy simulations. The reaction free energies and free energy barriers are found to be significantly influenced by the environment. The wG-T→G-T* tautomerization is predicted to be endoergic in aqueous solution and the DNA duplexes but slightly exoergic in the polymerase, with Arg517 and Asn513 providing electrostatic stabilization of G-T*. The G-T*→G*-T tautomerization is also predicted to be slightly more thermodynamically favorable in the polymerase relative to these DNA duplexes. These simulations are consistent with an experimentally driven kinetic misincorporation model suggesting that G-T mispair tautomerization occurs in the ajar polymerase conformation or concertedly with the transition from the ajar to the closed polymerase conformation. Furthermore, the order of the associated two proton transfer reactions is predicted to be different in the polymerase than in aqueous solution and the DNA duplexes. These studies highlight the impact of the environment on the thermodynamics, kinetics, and fundamental mechanisms of G-T mispair tautomerization, which plays a role in a wide range of biochemically important processes.

Dataset Information

Mechanistic implications of the cysteine-nicotinamide adduct in aldehyde dehydrogenase based on quantum mechanical/molecular mechanical simulations.

Publications

Mechanistic implications of the cysteine-nicotinamide adduct in aldehyde dehydrogenase based on quantum mechanical/molecular mechanical simulations.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets