Free energy of charge transfer and intraprotein electric field: method of calculation depends on the charge state of protein at a given structure

Semi-continuum electrostatic calculations of redox potentials in photosystem I

The midpoint redox potentials (E(m)) of all cofactors in photosystem I from Synechococcus elongatus as well as of the iron-sulfur (Fe(4)S(4)) clusters in two soluble ferredoxins from Azotobacter vinelandii and Clostridium acidiurici were calculated within the framework of a semi-continuum dielectric approach. The widely used treatment of proteins as uniform media with single dielectric permittivity is oversimplified, particularly, because permanent charges are considered both as a source for intraprotein electric field and as a part of dielectric polarizability. Our approach overcomes this inconsistency by using two dielectric constants: optical epsilon(o)=2.5 for permanent charges pre-existing in crystal structure, and static epsilon(s) for newly formed charges. We also take into account a substantial dielectric heterogeneity of photosystem I revealed by photoelectric measurements and a liquid junction potential correction for E(m) values of relevant redox cofactors measured in aprotic solvents. We show that calculations based on a single permittivity have the discrepancy with experimental data larger than 0.7 V, whereas E(m) values calculated within our approach fall in the range of experimental estimates. The electrostatic analysis combined with quantum chemistry calculations shows that (i) the energy decrease upon chlorophyll dimerization is essential for the downhill mode of primary charge separation between the special pair P(700) and the primary acceptor A(0); (ii) the primary donor is apparently P(700) but not a pair of accessory chlorophylls; (iii) the electron transfer from the A branch quinone Q(A) to the iron-sulfur cluster F(X) is most probably downhill, whereas that from the B branch quinone Q(B) to F(X) is essentially downhill.

Catalysis of electron transfer during activation of O2 by the flavoprotein glucose oxidase

Two prototropic forms of glucose oxidase undergo aerobic oxidation reactions that convert FADH(-) to FAD and form H(2)O(2) as a product. Limiting rate constants of k(cat)K(M)(O(2)) = (5.7 +/- 1.8) x 10(2) M(-1).s(-1) and k(cat)K(M)(O(2)) = (1.5 +/- 0.3) x 10(6) M(-1).s(-1) are observed at high and low pH, respectively. Reactions exhibit oxygen-18 kinetic isotope effects but no solvent kinetic isotope effects, consistent with mechanisms of rate-limiting electron transfer from flavin to O(2). Site-directed mutagenesis studies reveal that the pH dependence of the rates is caused by protonation of a highly conserved histidine in the active site. Temperature studies (283-323 K) indicate that protonation of His-516 results in a reduction of the activation energy barrier by 6.0 kcal.mol(-1) (0.26 eV). Within the context of Marcus theory, catalysis of electron transfer is attributed to a 19-kcal.mol(-1) (0.82 eV) decrease in the reorganization energy and a much smaller 2.2-kcal.mol(-1) (0.095 eV) enhancement of the reaction driving force. An explanation is advanced that is based on changes in outer-sphere reorganization as a function of pH. The active site is optimized at low pH, but not at high pH or in the H516A mutant where rates resemble the uncatalyzed reaction in solution.

Macromolecular electrostatics: continuum models and their growing pains

Theoretical understanding of macromolecular electrostatics has advanced substantially over the past year. Continuum models have given promising results for calculating protein-ligand binding free energy differences, as well as pK(a)s and redox properties, particularly with explicit treatment of multiple conformers. Generalized Born and other techniques have led to the first molecular dynamics simulations of proteins and RNA with continuum solvent. Continuum and microscopic descriptions of dielectric relaxation have been critically compared.