Dataset Information

Accurate Quantum Mechanical/Molecular Mechanical Calculations of Reduction Potentials in Azurin Variants.

ABSTRACT: Understanding the regulation mechanism and molecular determinants of the reduction potential of metalloprotein is a major challenge. An ab initio quantum mechanical/molecular mechanical (QM/MM) method combining the minimum free energy path (MFEP) and fractional number of electron (FNE) approaches has been developed in our group to simulate the redox processes of large systems. The FNE scheme provides an efficient unique description for the redox process, while the MFEP method provides improved conformational sampling on complex environments such as protein in the QM/MM calculations. The reduction potentials of wild-type and seven mutants of azurin, a type 1 copper metalloprotein, were simulated with the QM/MM-MFEP+FNE approach in this paper. A range of 350 mV for the variations of the reduction potentials of these azurin proteins was reproduced faithfully with relative errors around 20 mV. The correlation between structural interactions and reduction potentials observed in simulations provides in-depth insight into the regulation of reduction potentials, which potentially can also be very useful to the engineering of metalloprotein-based electrocatalysts in artificial photosynthesis. The excellent accuracy and efficiency of the QM/MM-MFEP+FNE approach demonstrate the potential for simulations of many electron transfer processes in condensed phases and biochemical systems.

SUBMITTER: Shen L

PROVIDER: S-EPMC6195766 | biostudies-literature | 2018 Sep

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Accurate Quantum Mechanical/Molecular Mechanical Calculations of Reduction Potentials in Azurin Variants.

Shen Lin L Zeng Xiancheng X Hu Hao H Hu Xiangqian X Yang Weitao W

Journal of chemical theory and computation 20180810 9

Understanding the regulation mechanism and molecular determinants of the reduction potential of metalloprotein is a major challenge. An ab initio quantum mechanical/molecular mechanical (QM/MM) method combining the minimum free energy path (MFEP) and fractional number of electron (FNE) approaches has been developed in our group to simulate the redox processes of large systems. The FNE scheme provides an efficient unique description for the redox process, while the MFEP method provides improved c ...[more]

PMID: 30040901

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.

Dataset Information

Accurate Quantum Mechanical/Molecular Mechanical Calculations of Reduction Potentials in Azurin Variants.

Publications

Accurate Quantum Mechanical/Molecular Mechanical Calculations of Reduction Potentials in Azurin Variants.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets