The Asn 38−Cys 84 H-Bond in Plastocyanin

Conformational Analyses of Physiological Binary and Ternary Copper(II) Complexes with l-Asparagine and l-Histidine; Study of Tridentate Binding of Copper(II) in Aqueous Solution

This study explores the structural properties and energy landscapes of the physiologically important bis(l-asparaginato)copper(II) [Cu(l-Asn)2] and (l-histidinato)(l-asparaginato)copper(II) [Cu(l-His)(l-Asn)]. The conformational analyses in the gas phase and implicitly modeled water medium, and magnetic parameters of electron paramagnetic resonance spectra were attained using density functional theory calculations. The apical CuII coordination and hydrogen bonding were analyzed. Predicted lower-energy structures enabled the confirmation and, for apical bonding, also the refinement of structural proposals from literature. Available experimental results were indecisive regarding the amido-group binding in the CuII equatorial plane in solutions, but the examination of the relative stability of Cu(l-Asn)2 conformers in 30 binding modes confirms the glycine-like mode as the most stable one. Previously reported experimental results for Cu(l-His)(l-Asn) were interpreted for l-His to have a tridentate histamine-like mode. However, the aqueous conformers with l-His in the glycinato mode are also predicted to have low energies, which does not contradict the tridentate l-His binding. The predicted magnetic parameters of conformers with an apical oxygen atom (intramolecular or from a water molecule) can reproduce the experimental data. An extent of conformational flexibility and abundance of l-His-containing ternary copper(II) amino acid complexes under physiological conditions may be related.

Glutamate Ligation in the Ni(II)- and Co(II)-Responsive Escherichia coli Transcriptional Regulator, RcnR

Escherichia coli RcnR (resistance to cobalt and nickel regulator, EcRcnR) is a metal-responsive repressor of the genes encoding the Ni(II) and Co(II) exporter proteins RcnAB by binding to P_RcnAB. The DNA binding affinity is weakened when the cognate ions Ni(II) and Co(II) bind to EcRcnR in a six-coordinate site that features a (N/O)₅S ligand donor-atom set in distinct sites: while both metal ions are bound by the N terminus, Cys35, and His64, Co(II) is additionally bound by His3. On the other hand, the noncognate Zn(II) and Cu(I) ions feature a lower coordination number, have a solvent-accessible binding site, and coordinate protein ligands that do not include the N-terminal amine. A molecular model of apo-EcRcnR suggested potential roles for Glu34 and Glu63 in binding Ni(II) and Co(II) to EcRcnR. The roles of Glu34 and Glu63 in metal binding, metal selectivity, and function were therefore investigated using a structure/function approach. X-ray absorption spectroscopy was used to assess the structural changes in the Ni(II), Co(II), and Zn(II) binding sites of Glu → Ala and Glu → Cys variants at both positions. The effect of these structural alterations on the regulation of P_rcnA by EcRcnR in response to metal binding was explored using LacZ reporter assays. These combined studies indicate that while Glu63 is a ligand for both metal ions, Glu34 is a ligand for Co(II) but possibly not for Ni(II). The Glu34 variants affect the structure of the cognate metal sites, but they have no effect on the transcriptional response. In contrast, the Glu63 variants affect both the structure and transcriptional response, although they do not completely abolish the function of EcRcnR. The structure of the Zn(II) site is not significantly perturbed by any of the glutamic acid variations. The spectroscopic and functional data obtained on the mutants were used to calculate models of the metal-site structures of EcRcnR bound to Ni(II), Co(II), and Zn(II). The results are interpreted in terms of a switch mechanism, in which a subset of the metal-binding ligands is responsible for the allosteric response required for DNA release.

Multiconfigurational g tensor calculations as a probe for the covalency of the copper-ligand bonds in copper(II) complexes: [CuCl4]2-, [Cu(NH3)4]2+, and plastocyanin

Calculations of the g tensor of three copper(II) complexes [Cu(NH3)4]2+, [CuCl4]2-, and plastocyanin are presented. Two different sum-over-states-based approaches are considered, making use of the multistate CASPT2 method for excitation energies and PMCAS (perturbation modified CAS) wave functions for the computation of the angular momentum and spin-orbit coupling matrix elements. Test calculations on [Cu(NH3)4]2+ and [CuCl4]2- point to the need of including in the MS-CASPT2 treatment the specific charge-transfer state with an electron excited out of the bonding counterpart of the ground-state SOMO. The computed g shifts for these two molecules present a considerable improvement with respect to the results obtained from our previous g tensor calculations based instead on CASSCF/CASPT2. This is shown to be related to an improved description of the covalency of the Cu-L bonds. For the calculations on plastocyanin, different models are used, taken from a recent (QM/MM) DFT study by Sinnecker and Neese. The effect of the surrounding protein is taken into account by surrounding the central cluster either with a dielectric continuum (epsilon = 4) or with a set of point charges. The second approach is found to be indispensable for an accurate description of environmental effects. With this approach, the calculated g values compare to within 30 ppt with the experimental data of plastocyanin.

NMR hyperfine shifts in blue copper proteins: a quantum chemical investigation

We present the results of the first quantum chemical investigations of 1H NMR hyperfine shifts in the blue copper proteins (BCPs): amicyanin, azurin, pseudoazurin, plastocyanin, stellacyanin, and rusticyanin. We find that very large structural models that incorporate extensive hydrogen bond networks, as well as geometry optimization, are required to reproduce the experimental NMR hyperfine shift results, the best theory vs experiment predictions having R2 = 0.94, a slope = 1.01, and a SD = 40.5 ppm (or approximately 4.7% of the overall approximately 860 ppm shift range). We also find interesting correlations between the hyperfine shifts and the bond and ring critical point properties computed using atoms-in-molecules theory, in addition to finding that hyperfine shifts can be well-predicted by using an empirical model, based on the geometry-optimized structures, which in the future should be of use in structure refinement.

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.

Basic requirements for a metal-binding site in a protein: the influence of loop shortening on the cupredoxin azurin

The main active-site loop of the copper-binding protein azurin (a cupredoxin) has been shortened from C(112)TFPGH(117)SALM(121) to C(112)TPH(115)PFM(118) (the native loop from the cupredoxin amicyanin) and also to C(112)TPH(115)PM(117). The Cu(II) site structure is almost unaffected by shortening, as is that of the Cu(I) center at alkaline pH in the variant with the C(112)TPH(115)PM(117) loop sequence. Subtle spectroscopic differences due to alterations in the spin density distribution at the Cu(II) site can be attributed mainly to changes in the hydrogen-bonding pattern. Electron transfer is almost unaffected by the introduction of the C(112)TPH(115)PFM(118) loop, but removal of the Phe residue has a sizable effect on reactivity, probably because of diminished homodimer formation. At mildly acidic pH values, the His-115 ligand protonates and dissociates from the cuprous ion, an effect that has a dramatic influence on the reactivity of cupredoxins. These studies demonstrate that the amicyanin loop adopts a conformation identical to that found in the native protein when introduced into azurin, that a shorter than naturally occurring C-terminal active-site loop can support a functional T1 copper site, that CTPHPM is the minimal loop length required for binding this ubiquitous electron transfer center, and that the length and sequence of a metal-binding loop regulates a range of structural and functional features of the active site of a metalloprotein.