Predicting Reduction Potentials of Blue Copper Proteins Using Quantum Mechanical Calculations.

Area of Science:

• Biophysical Chemistry
• Computational Chemistry
• Protein Science

Background:

• Blue copper proteins are vital metalloproteins involved in electron transfer.
• Accurate prediction of their redox potentials is crucial for understanding biological function.
• Computational methods offer a powerful tool for studying these properties.

Purpose of the Study:

• To systematically evaluate 64 computational methods for calculating blue copper protein redox potentials.
• To identify optimal computational parameters for predicting both relative and absolute redox potentials.
• To assess the accuracy and limitations of quantum mechanics/molecular mechanics (QM/MM) approaches.

Main Methods:

• Calculated redox potentials for 12 blue copper protein sites.
• Systematically varied quantum mechanics (QM) system size, dielectric constants, density functionals, and basis sets.
• Utilized QM/MM for structure optimization and QM-cluster calculations in a continuum solvent.

Main Results:

• An intermediate QM system size (∼70 atoms) with TPSS functional and SV(P) basis set yielded the best relative potentials (MAD = 0.09 V).
• Larger QM systems (∼340 atoms) with B3LYP functional and larger basis sets improved absolute potential accuracy (MSE = -0.27 V).
• Blue copper proteins showed higher accuracy than iron-sulfur clusters due to simpler coordination.

Conclusions:

• QM-cluster calculations in continuum solvent effectively balance accuracy and computational cost.
• Method selection is critical and depends on whether relative or absolute potentials are targeted.
• Accurate modeling of redox-active sites in proteins remains challenging but achievable with optimized computational strategies.