Premium
Prediction of Reduction Potentials of Copper Proteins with Continuum Electrostatics and Density Functional Theory
Author(s) -
Fowler Nicholas J.,
Blanford Christopher F.,
Warwicker Jim,
de Visser Sam P.
Publication year - 2017
Publication title -
chemistry – a european journal
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.687
H-Index - 242
eISSN - 1521-3765
pISSN - 0947-6539
DOI - 10.1002/chem.201702901
Subject(s) - azurin , electrostatics , reduction (mathematics) , density functional theory , chemistry , mutant , static electricity , copper , electron transfer , molecular dynamics , work (physics) , chemical physics , copper protein , biological system , computational chemistry , physics , thermodynamics , mathematics , biology , quantum mechanics , biochemistry , geometry , organic chemistry , gene
Blue copper proteins, such as azurin, show dramatic changes in Cu 2+ /Cu + reduction potential upon mutation over the full physiological range. Hence, they have important functions in electron transfer and oxidation chemistry and have applications in industrial biotechnology. The details of what determines these reduction potential changes upon mutation are still unclear. Moreover, it has been difficult to model and predict the reduction potential of azurin mutants and currently no unique procedure or workflow pattern exists. Furthermore, high‐level computational methods can be accurate but are too time consuming for practical use. In this work, a novel approach for calculating reduction potentials of azurin mutants is shown, based on a combination of continuum electrostatics, density functional theory and empirical hydrophobicity factors. Our method accurately reproduces experimental reduction potential changes of 30 mutants with respect to wildtype within experimental error and highlights the factors contributing to the reduction potential change. Finally, reduction potentials are predicted for a series of 124 new mutants that have not yet been investigated experimentally. Several mutants are identified that are located well over 10 Å from the copper center that change the reduction potential by more than 85 mV. The work shows that secondary coordination sphere mutations mostly lead to long‐range electrostatic changes and hence can be modeled accurately with continuum electrostatics.