Hydrogen Bonding in the Blue-Copper Site. Resonance Raman Study

The structure of plastocyanin tunes the midpoint potential by restricting axial ligation of the reduced copper ion

Blue copper proteins are models for illustrating how proteins tune metal properties. Nevertheless, the mechanisms by which the protein controls the metal site remain to be fully elucidated. A hindrance is that the closed shell Cu(I) site is inaccessible to most spectroscopic analyses. Carbon deuterium (C-D) bonds used as vibrational probes afford nonperturbative, selective characterization of the key cysteine and methionine copper ligands in both redox states. The structural integrity of Nostoc plastocyanin was perturbed by disrupting potential hydrogen bonds between loops of the cupredoxin fold via mutagenesis (S9A, N33A, N34A), variably raising the midpoint potential. The C-D vibrations show little change to suggest substantial alteration to the Cu(II) coordination in the oxidized state or in the Cu(I) interaction with the cysteine ligand. They rather indicate, along with visible and NMR spectroscopy, that the methionine ligand distinctly interacts more strongly with the Cu(I) ion, in line with the increases in midpoint potential. Here we show that the protein structure determines the redox properties by restricting the interaction between the methionine ligand and Cu(I) in the reduced state.

Systematic elucidation of the second coordination sphere effect on the structure and properties of a blue copper protein, pseudoazurin

The rational structural and computational studies of a blue copper protein, pseudoazurin (PAz), and its Met16X (X = Phe, Leu, Val, Ile) variants gave clear functional meanings of the noncovalent interaction (NCI) through the second coordination sphere. The high-resolution X-ray crystal structures of Met16X PAz demonstrated that the active site geometry is significantly affected by the substitution of Met16, which is located within the NCI distance from the His81 imidazole ring at the copper active site. The computational chemistry calculations based on the crystal structure analyses confirmed that the NCI of S-π/CH-π (wild-type), π-π (Met16Phe), double CH-π (Met16Leu), and single CH-π (Met16Val and Met16Ile). The estimated interaction energies for the NCI demonstrated that the fine-tuning of the protein stability and Cu site properties form the second coordination sphere of PAz.

Engineering a bifunctional copper site in the cupredoxin fold by loop-directed mutagenesis

Copper sites in proteins are designed to perform either electron transfer or redox catalysis. Type 1 and CuA sites are electron transfer hubs bound to a rigid protein fold that prevents binding of exogenous ligands and side reactions. Here we report the engineering of two Type 1 sites by loop-directed mutagenesis within a CuA scaffold with unique electronic structures and functional features. A copper-thioether axial bond shorter than the copper-thiolate bond is responsible for the electronic structure features, in contrast to all other natural or chimeric sites where the copper thiolate bond is short. These sites display highly unusual features, such as: (1) a high reduction potential despite a strong interaction with the axial ligand, which we attribute to changes in the hydrogen bond network and (2) the ability to bind exogenous ligands such as imidazole and azide. This strategy widens the possibility of using natural protein scaffolds with functional features not present in nature.

QM/MM calculations with DFT for taking into account protein effects on the EPR and optical spectra of metalloproteins. Plastocyanin as a case study

A detailed study of the influence of the surrounding protein on magnetic and optical spectra of metalloproteins is presented using the quantum-mechanical/molecular mechanical (QM/MM) approach. The well-studied type I copper site in plastocyanin in the cupric oxidation state is taken as a test case because its spectroscopic properties have been extensively studied and are well understood. The calculations have been performed using nonrelativistic and scalar relativistic (at the level of the zeroth order regular approximation, ZORA) calculations (B3LYP functional). Linear response theory has been used to calculate first- and second-order properties, namely the EPR g-tensor, the central metal hyperfine couplings (HFCs), the HFCs of the directly coordinating ligands, as well as superhyperfine couplings (1H, 14N) from remote nuclei, transition energies, and oscillator strengths. Two different model systems have been defined that do not and do include important amino acids from the second coordination sphere, respectively. For comparison, calculations have been carried out in the gas phase and in a dielectric continuum (conductor like screening model, COSMO) with a dielectric constant of four. The best results were obtained at the scalar relativistic ZORA level for the largest model in conjunction with explicit modeling of the protein environment through the QM/MM procedure, which is also considered to be the highest level of theory used in this work. The protein effects beyond the second coordination sphere were found to be quite substantial (up to 30% changes on some properties), and were found to require an explicit treatment of the protein beyond the second coordination sphere. In addition, the embedding water cage was found to have a nonnegligible influence on the calculated spectroscopic data, which is of the same order as the influence of the protein backbone charges. However, while qualitatively satisfactory, the errors in the calculated spectroscopic parameters are still substantial, and can all be traced back to the fact that the linear-response of the presently available functionals is "too stiff" with respect to the external perturbations at least for the model systems studied here. Ligand field-based approaches are used to correct for systematic errors in the DFT procedures. As a consequence, we propose a new breakdown of the copper hyperfine interaction into Fermi-contact, spin-dipolar and spin-orbit contributions.

The role of hydrogen bonding at the active site of a cupredoxin: the Phe114Pro azurin variant

The Phe114Pro mutation to the cupredoxin azurin (AZ) leads to a number of structural changes at the active site attributed to deletion of one of the hydrogen bonds to the Cys112 ligand, removal of the bulky phenyl group from the hydrophobic patch of the protein, and steric interactions made by the introduced Pro. The remaining hydrogen bond between the coordinating thiolate and the backbone amide of Asn47 is strengthened. At the type-1 copper site, the Cu(II)-O(Gly45) axial interaction decreases, while the metal moves out of the plane formed by the equatorial His46, Cys112, and His117 ligands, shortening the bond to the axially coordinating Met121. The resulting distorted tetrahedral geometry is distinct from the trigonal bipyramidal arrangement in the wild-type (WT) protein. The unique position of the main S(Cys) --> Cu(II) ligand-to-metal charge-transfer transition in AZ (628 nm) has shifted in the Phe114Pro variant to a value that is more typical for cupredoxins (599 nm). This probably occurs because of the removal of the Phe114-Cys112 hydrogen bond. The Phe114Pro mutation results in a 90 mV decrease in the reduction potential of AZ, and removal of the second hydrogen bond to the Cys ligand seems to be the major cause of this change. The C-terminal His117 ligand does not protonate in the reduced Phe114Pro AZ variant, which suggests that none of the structural features altered by the mutation are responsible for the absence of this effect in the WT protein. Upon reduction, the copper displaces further from the equatorial ligand plane and the Cu-S(Met121) bond length decreases. These changes are larger than those seen in the WT protein and contribute to the order of magnitude decrease in the intrinsic electron-transfer capabilities of the Phe114Pro variant.