Ab Initio Quantum Mechanical Study of Metal Substitution in Analogues of Rubredoxin: Implications for Redox Potential Control

Quantitative Account of the Bonding Properties of a Rubredoxin Model Complex [Fe(SCH3)4]q, q = -2, -1, +2, +3

Iron-sulfur clusters play important roles in biology as parts of electron-transfer chains and catalytic cofactors. Here, we report a detailed computational analysis of a structural model of the simplest natural iron-sulfur cluster of rubredoxin and its cationic counterparts. Specifically, we investigated adiabatic reduction energies, dissociation energies, and bonding properties of the low-lying electronic states of the complexes [Fe(SCH₃)₄]^2-/1-/2+/3+ using multireference (CASSCF, MRCISD), and coupled cluster [CCSD(T)] methodologies. We show that the nature of the Fe-S chemical bond and the magnitude of the ionization potentials in the anionic and cationic [Fe(SCH₃)₄] complexes offer a physical rationale for the relative stabilization, structure, and speciation of these complexes. Anionic and cationic complexes present different types of chemical bonds: prevalently ionic in [Fe(SCH₃)₄]^2-/1- complexes and covalent in [Fe(SCH₃)₄]^2+/3+ complexes. The ionic bonds result in an energy gain for the transition [Fe(SCH₃)₄]^2- → [Fe(SCH₃)₄]^- (i.e., Fe^II → Fe^III) of 1.5 eV, while the covalent bonds result in an energy loss for the transition [Fe(SCH₃)₄]²⁺ → [Fe(SCH₃)₄]³⁺ of 16.6 eV, almost half of the ionization potential of Fe²⁺. The ionic versus covalent bond character influences the Fe-S bond strength and length, that is, ionic Fe-S bonds are longer than covalent ones by about 0.2 Å (for Fe^II) and 0.04 Å (for Fe^II). Finally, the average Fe-S heterolytic bond strength is 6.7 eV (Fe^II) and 14.6 eV (Fe^III) at the RCCSD(T) level of theory.

Theoretical study of DNA damage recognition via electron transfer from the [4Fe-4S] complex of MutY

The mechanism of site-specific recognition of DNA by proteins has been a long-standing issue. The DNA glycosylase MutY, for instance, must find the rare 8-oxoguanine-adenine mismatches among the large number of basepairs in the DNA. This protein has a [4Fe-4S] cluster, which is highly conserved in species as diverse as Escherichia Coli and Homo sapiens. The mixed-valent nature of this cluster suggests that charge transfer may play a role in MutY's function. We have studied the energetics of the charge transfer in Bacillus stearothermophilus MutY-DNA complex using multiscale calculation including density functional theory and molecular dynamics. The [4Fe-4S] cluster in MutY is found to undergo 2+ to 3+ oxidation when coupling to DNA through hole transfer, especially when MutY is near an oxoguanine modified base (oxoG). Employing the Marcus theory for electron transfer, we find near optimal Frank-Condon factors for electron transfer from MutY to oxoguanine modified base. MutY has modest selectivity for oxoguanine over guanine due to the difference in oxidation potential. The tunneling matrix element is significantly reduced with the mutation R149W, whereas the mutation L154F reduces the tunneling matrix element as well as the Frank-Condon factor. Both L154F and R149W mutations are known to dramatically reduce or eliminate repair efficiency. We suggest a scenario where the charge transfer leads to a stabilization of the specific binding conformation, which is likely the recognition mode, thus enabling it to find the damaged site efficiently.

A billion-fold range in acidity for the solvent-exposed amides of Pyrococcus furiosus rubredoxin

The exchange rates of the static solvent-accessible amide hydrogens of Pyrococcus furiosus rubredoxin range from near the diffusion-limited rate to a billion-fold slower for the non-hydrogen-bonded Val 38 (eubacterial numbering). Hydrogen exchange directly monitors the kinetic acidity of the peptide nitrogen. Electrostatic solvation free energies were calculated by Poisson-Boltzmann methods for the individual peptide anions that form during the hydroxide-catalyzed exchange reaction to examine how well the predicted thermodynamic acidities match the experimentally determined kinetic acidities. With the exception of the Ile 12 amide, the differential exchange rate constant for each solvent-exposed amide proton that is not hydrogen bonded to a backbone carbonyl can be predicted within a factor of 6 (10 (0.78)) root-mean-square deviation (rmsd) using the CHARMM22 electrostatic parameter set and an internal dielectric value of 3. Under equivalent conditions, the PARSE parameter set yields a larger rmsd value of 1.28 pH units, while the AMBER parm99 parameter set resulted in a considerably poorer correlation. Either increasing the internal dielectric value to 4 or reducing it to a value of 2 significantly degrades the quality of the prediction. Assigning the excess charge of the peptide anion equally between the peptide nitrogen and the carbonyl oxygen also reduces the correlation to the experimental data. These continuum electrostatic calculations were further analyzed to characterize the specific structural elements that appear to be responsible for the wide range of peptide acidities observed for these solvent-exposed amides. The striking heterogeneity in the potential at sites along the protein-solvent interface should prove germane to the ongoing challenge of quantifying the contribution that electrostatic interactions make to the catalytic acceleration achieved by enzymes.

Spatial distribution of dielectric shielding in the interior of Pyrococcus furiosus rubredoxin as sampled in the subnanosecond timeframe by hydrogen exchange

Experimental pK values of ionizable sidechains provide the most direct test for models representing dielectric shielding within the interior of a protein. However, only the strongly shifted pK values are particularly useful for discriminating among models. NMR titration studies have usually found only one or two such shifted pK values in each protein, so that the fitting of the experimental data to a uniform internal dielectric (epsilon(int)) model is not well constrained. The observed variation among proteins for such epsilon(int) estimates may reflect nonuniformity of dielectric shielding within each protein interior or qualitative differences between individual proteins. The differential amide kinetic acidities for a series of metal-substituted rubredoxins are shown to be consistent with Poisson-Boltzmann predictions of dielectric shielding that is relatively uniform for all of the amides that are sensitive to the metal charge, a region which corresponds to roughly 1/3 of the internal volume. The effective epsilon(int) values near 6 that are found in this study are significantly lower than many such estimates derived from sidechain pK measurements. The differing timeframes in which dielectric relaxation can respond to the highly transient peptide anion as compared to the longer lived states of the charged sidechains offers an explanation for the lower apparent dielectric constant deduced from these measurements.

Structure and Redox Properties of the Protein, Rubredoxin, and Its Ligand and Metal Mutants Studied by Electronic Structure Calculation

The changes in the structural and electronic properties accompanying metal ionization of the iron-containing protein, rubredoxin, and of some ligand and metal mutants, have been explored using density functional theory (DFT) calculations of the metal atom coordinated to the four immediate residues. Both isolated and embedded cluster studies have been carried out, the latter using the hybrid quantum mechanics/molecular mechanics (QM/MM) approach. The replacement of a cysteine by a serine residue has a considerable effect on both the electronic and geometric structure of the core, which can be qualitatively understood on the basis of the isolated cluster studies. The modulation of these properties caused by the protein environment is quite accurately described by the QM/MM calculations. The predicted core geometries are in good accord with both X-ray and EXAFS data, and the changes in the redox potentials are predicted, at least semiquantitatively, by considering only the core part of the protein.

Calculation of the Redox Potential of the Protein Azurin and Some Mutants

Azurin from Pseudomonas aeruginosa is a small 128-residue, copper-containing protein. Its redox potential can be modified by mutating the protein. Free-energy calculations based on classical molecular-dynamics simulations of the protein and from mutants in aqueous solution at different pH values were used to compute relative redox potentials. The precision of the free-energy calculations with the lambda coupling-parameter approach is evaluated as function of the number and sequence of lambda values, the sampling time and initial conditions. It is found that the precision is critically dependent on the relaxation of hydrogen-bonding networks when changing the atomic-charge distribution due to a change of redox state or pH value. The errors in the free energies range from 1 to 10 k(B)T, depending on the type of process. Only qualitative estimates of the change in redox potential by protein mutation can be obtained.