Spectroscopic and Geometric Variations in Perturbed Blue Copper Centers: Electronic Structures of Stellacyanin and Cucumber Basic Protein

Reparameterization of Non-Bonded Parameters for Copper Ions in Plastocyanin: An Adaptive Force Matching Study

Molecular mechanics rely on existing experimental and theoretical inputs to confidently calculate the trajectories of molecular systems. These calculations, however, are often hindered by missing force field parameters. A notable subject of this problem is metal centers of proteins. This study parameterized, through an adaptive force matching (AFM) workflow, the copper cofactor of plastocyanin in its two oxidation states. New 12-6 Lennard-Jones (LJ) parameters and atomic partial charges were generated to complete the non-bonded description of the copper site. Our models show uniform distorted tetrahedral structures for reduced plastocyanin, Cu(I), and oxidized plastocyanin, Cu(II). These structures align with the QM/MM MD results and existing crystallography studies. TD-DFT calculations, meanwhile, showed that conformations with elongated axial Cu-SMet and shortened equatorial Cu-SCys bonds retain the experimental UV-Vis profile of blue copper (BC) proteins, thus signifying the importance of Cu-S interactions on BC proteins' unique spectroscopic properties.

Biochemical Characterization of the Copper Nitrite Reductase from Neisseria gonorrhoeae

The copper-containing nitrite reductase from Neisseria gonorrhoeae has been shown to play a critical role in the infection mechanism of this microorganism by producing NO and abolishing epithelial exfoliation. This enzyme is a trimer with a type 1 copper center per subunit and a type 2 copper center in the subunits interface, with the latter being the catalytic site. The two centers were characterized for the first time by EPR and CD spectroscopy, showing that the type 1 copper center has a high rhombicity due to its lower symmetry and more tetragonal structure, while the type 2 copper center has the usual properties, but with a smaller hyperfine coupling constant (A// = 10.5 mT). The thermostability of the enzyme was analyzed by differential scanning calorimetry, which shows a single endothermic transition in the thermogram, with a maximum at 94 °C, while the CD spectra in the visible region indicate the presence of the type 1 copper center up to 80 °C. The reoxidation of the N. gonorrhoeae copper-containing nitrite reductase in the presence of nitrite were analyzed by visible spectroscopy and showed a pH dependence, being higher at pH 5.5-6.0. The high thermostability of this enzyme may be important to maintaining a high activity in the extracellular space and to making it less susceptible to denaturation and proteolysis, contributing to the proliferation of N. gonorrhoeae.

Orchestrating copper binding: structure and variations on the cupredoxin fold

A large number of copper binding proteins coordinate metal ions using a shared three-dimensional fold called the cupredoxin domain. This domain was originally identified in Type 1 "blue copper" centers but has since proven to be a common domain architecture within an increasingly large and diverse group of copper binding domains. The cupredoxin fold has a number of qualities that make it ideal for coordinating Cu ions for purposes including electron transfer, enzyme catalysis, assembly of other copper sites, and copper sequestration. The structural core does not undergo major conformational changes upon metal binding, but variations within the coordination environment of the metal site confer a range of Cu-binding affinities, reduction potentials, and spectroscopic properties. Here, we discuss these proteins from a structural perspective, examining how variations within the overall cupredoxin fold and metal binding sites are linked to distinct spectroscopic properties and biological functions. Expanding far beyond the blue copper proteins, cupredoxin domains are used by a growing number of proteins and enzymes as a means of binding copper ions, with many more likely remaining to be identified.

Spiers Memorial Lecture: Activating Metal Sites for Biological Electron Transfer

Metal sites in biology often exhibit unique spectroscopic features that reflect novel geometric and electronic structures imposed by the protein that are key to reactivity. The Blue copper active site...

Electrostatic regulation of blue copper sites

In the last 50 years, the blue copper proteins became central targets of investigation. Extensive experiments focused on the Cu coordination to probe the effect of local perturbations on its properties. We found that local electric fields, generated by charged residues evolutionarily placed throughout the protein edifice, mainly second sphere, but also more remotely, constitute an additional significant factor regulating blue copper proteins. These fields are not random, but exhibit a highly specific directionality, negative with respect to the and vectors in the Cu first shell. The field magnitude contributes to fine-tuning of the geometric and electronic properties of Cu sites in individual blue copper proteins. Specifically, the local electric fields evidently control the Cu–S_Met bond distance, Cu(ii)–S_Cys bond covalency, and the energies of the frontier molecular orbitals, which, in turn, govern the Cu(ii/i) reduction potential and the relative absorption intensities at 450 nm and 600 nm.

Intramolecular electric fields in blue copper proteins are oriented in a fixed way to modulate properties of their copper sites: they control the first-shell copper interactions to influence geometric, spectroscopic, and redox behavior.

A [3Cu:2S] cluster provides insight into the assembly and function of the CuZ site of nitrous oxide reductase

Nitrous oxide reductase (N₂OR) is the only known enzyme reducing environmentally critical nitrous oxide (N₂O) to dinitrogen (N₂) as the final step of bacterial denitrification. The assembly process of its unique catalytic [4Cu:2S] cluster Cu_Z remains scarcely understood. Here we report on a mutagenesis study of all seven histidine ligands coordinating this copper center, followed by spectroscopic and structural characterization and based on an established, functional expression system for Pseudomonas stutzeri N₂OR in Escherichia coli. While no copper ion was found in the Cu_Z binding site of variants H129A, H130A, H178A, H326A, H433A and H494A, the H382A variant carried a catalytically inactive [3Cu:2S] center, in which one sulfur ligand, S_Z2, had relocated to form a weak hydrogen bond to the sidechain of the nearby lysine residue K454. This link provides sufficient stability to avoid the loss of the sulfide anion. The UV-vis spectra of this cluster are strikingly similar to those of the active enzyme, implying that the flexibility of S_Z2 may have been observed before, but not recognized. The sulfide shift changes the metal coordination in Cu_Z and is thus of high mechanistic interest.

Variants of all seven histidine ligands of the [4Cu:2S] active site of nitrous oxide reductase mostly result in loss of the metal site. However, a H382A variant retains a [3Cu:2S] cluster that hints towards a structural flexibility also present in the intact site.

Structural Basis for a Quadratic Relationship between Electronic Absorption and Electronic Paramagnetic Resonance Parameters of Type 1 Copper Proteins

Type 1 copper (T1Cu) proteins play important roles in electron transfer in biology, largely due to the unique structure of the T1Cu center, which is reflected by its spectroscopic properties. Previous reports have suggested a correlation between a high ratio of electronic absorbance at ∼450 nm to that at ∼600 nm (R = A₄₅₀/A₆₀₀) and a large copper(II) hyperfine coupling in the z direction (A_z) in electron paramagnetic resonance (EPR). However, this correlation does not have a clear physical meaning, nor does it hold for many proteins with a perturbed T1Cu center. To address this issue, a new parameter of R' [A₄₅₀/(A₄₅₀ + A₆₀₀)] with a better physical meaning of a fractional S_Cys pseudo-σ to Cu(II) charge transfer transition intensity is defined and a quadratic relationship between R' and A_z is found on the basis of a comprehensive analysis of ultraviolet-visible absorption, EPR, and structural parameters of T1Cu proteins. We are able to find good correlations between R' and the displacement of copper from the trigonal plane defined by the His₂Cys ligands and the angle between the N_His1-Cu-N_His2 plane and the S_Cys-Cu-axial ligand plane, providing a structural basis for the observed correlation. These findings and analyses provide a new framework for a deeper understanding of the spectroscopic and electronic properties of T1Cu proteins, which may allow better design and applications of this important class of proteins for redox and electron transfer functions.

CuA-based chimeric T1 copper sites allow for independent modulation of reorganization energy and reduction potential

Attaining rational modulation of thermodynamic and kinetic redox parameters of metalloproteins is a key milestone towards the (re)design of proteins with new or improved redox functions. Here we report that implantation of ligand loops from natural T1 proteins into the scaffold of a CuA protein leads to a series of distorted T1-like sites that allow for independent modulation of reduction potentials (E°') and electron transfer reorganization energies (λ). On the one hand E°' values could be fine-tuned over 120 mV without affecting λ. On the other, λ values could be modulated by more than a factor of two while affecting E°' only by a few millivolts. These results are in sharp contrast to previous studies that used T1 cupredoxin folds, thus highlighting the importance of the protein scaffold in determining such parameters.

Effect of NaOH concentration on antibacterial activities of Cu nanoparticles and the antibacterial mechanism

A series of Cu nanoparticles (NPs) have been prepared by a facile hydrothermal method at 180 °C using different concentrations of NaOH solutions and characterized by XRD, SEM, TEM and FT-IR spectra. Their antibacterial activities were assessed by means of Gram-positive S. aureus and Gram-negative E. coli bacteria, where various dosages (3, 5, 7, 10 mg) of the antibacterial agents were applied, and compared with that of the commercial CuSO₄ salt. The antibacterial mechanism was explored based on series of control experiments. The results show that the NaOH concentration affects the crystallinity, crystal size and surface hydroxyl content of the Cu NPs, which significantly influence the antibacterial activities. Compared to the commercial CuSO₄ salt, the four Cu samples prepared using no <4 mol L^-1 of NaOH display excellent antibacterial activities with low concentrations of copper leachates, which is great beneficial to the practical applications. The experimental results support that the highly reactive and soluble copper species in the antibacterial system of the Cu NPs is a Cu (II)-peptide complex, but not free Cu²⁺ ions.

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.

A QM/MM study of the nature of the entatic state in plastocyanin

Plastocyanin is a copper containing protein that is involved in the electron transfer process in photosynthetic organisms. The active site of plastocyanin is described as an entatic state whereby its structure represents a compromise between the structures favored by the oxidized and reduced forms. In this study, the nature of the entatic state is investigated through density functional theory‐based hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations. The strain energy is computed to be 12.8 kcal/mol and 14.5 kcal/mol for the oxidized and reduced forms of the protein, indicating that the active site has an intermediate structure. It is shown that the energy gap between the oxidized and reduced forms varies significantly with the fluctuations in the structure of the active site at room temperature. An accurate determination of the reorganization energy requires averaging over conformation and a large region of the protein around the active site to be treated at the quantum mechanical level. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Activating Metal Sites for Biological Electron Transfer

This review focuses on the unique spectroscopic features of the blue copper active sites. These reflect a novel electronic structure that activates the site for rapid long-range electron transfer in its biological function. The role of the protein in determining the geometric and electronic structure of this site is defined, as is its contribution to function. This has been referred to as the entatic/rack-induced state. These concepts are then extended to cytochrome c, which is also determined to be in an entatic state.

Modular Artificial Cupredoxins

Cupredoxins are electron-transfer proteins that have active sites containing a mononuclear Cu center with an unusual trigonal monopyramidal structure (Type 1 Cu). A single Cu-Scys bond is present within the trigonal plane that is responsible for its unique physical properties. We demonstrate that a cysteine-containing variant of streptavidin (Sav) can serve as a protein host to model the structure and properties of Type 1 Cu sites. A series of artificial Cu proteins are described that rely on Sav and a series of biotinylated synthetic Cu complexes. Optical and EPR measurements highlight the presence of a Cu-Scys bond, and XRD analysis provides structural evidence. We further provide evidence that changes in the linker between the biotin and Cu complex within the synthetic constructs allows for small changes in the placement of Cu centers within Sav that have dramatic effects on the structural and physical properties of the resulting artificial metalloproteins. These findings highlight the utility of the biotin-Sav technology as an approach for simulating active sites of metalloproteins.

Electron transfer and reaction mechanism of laccases

Laccases are part of the family of multicopper oxidases (MCOs), which couple the oxidation of substrates to the four electron reduction of O2 to H2O. MCOs contain a minimum of four Cu's divided into Type 1 (T1), Type 2 (T2), and binuclear Type 3 (T3) Cu sites that are distinguished based on unique spectroscopic features. Substrate oxidation occurs near the T1, and electrons are transferred approximately 13 Å through the protein via the Cys-His pathway to the T2/T3 trinuclear copper cluster (TNC), where dioxygen reduction occurs. This review outlines the electron transfer (ET) process in laccases, and the mechanism of O2 reduction as elucidated through spectroscopic, kinetic, and computational data. Marcus theory is used to describe the relevant factors which impact ET rates including the driving force, reorganization energy, and electronic coupling matrix element. Then, the mechanism of O2 reaction is detailed with particular focus on the intermediates formed during the two 2e(-) reduction steps. The first 2e(-) step forms the peroxide intermediate, followed by the second 2e(-) step to form the native intermediate, which has been shown to be the catalytically relevant fully oxidized form of the enzyme.

Anisotropic covalency contributions to superexchange pathways in type one copper active sites

Type one (T1) Cu sites deliver electrons to catalytic Cu active sites: the mononuclear type two (T2) Cu site in nitrite reductases (NiRs) and the trinuclear Cu cluster in the multicopper oxidases (MCOs). The T1 Cu and the remote catalytic sites are connected via a Cys-His intramolecular electron-transfer (ET) bridge, which contains two potential ET pathways: P1 through the protein backbone and P2 through the H-bond between the Cys and the His. The high covalency of the T1 Cu-S(Cys) bond is shown here to activate the T1 Cu site for hole superexchange via occupied valence orbitals of the bridge. This covalency-activated electronic coupling (H(DA)) facilitates long-range ET through both pathways. These pathways can be selectively activated depending on the geometric and electronic structure of the T1 Cu site and thus the anisotropic covalency of the T1 Cu-S(Cys) bond. In NiRs, blue (π-type) T1 sites utilize P1 and green (σ-type) T1 sites utilize P2, with P2 being more efficient. Comparing the MCOs to NiRs, the second-sphere environment changes the conformation of the Cys-His pathway, which selectively activates HDA for superexchange by blue π sites for efficient turnover in catalysis. These studies show that a given protein bridge, here Cys-His, provides different superexchange pathways and electronic couplings depending on the anisotropic covalencies of the donor and acceptor metal sites.

Metalloproteins diversified: the auracyanins are a family of cupredoxins that stretch the spectral and redox limits of blue copper proteins

The metal sites of electron transfer proteins are tuned for function. The type 1 copper site is one of the most utilized metal sites in electron transfer reactions. This site can be tuned by the protein environment from +80 mV to +680 mV in typical type 1 sites. Accompanying this huge variation in midpoint potentials are large changes in electronic structure, resulting in proteins that are blue, green, or even red. Here, we report a family of blue copper proteins, the auracyanins, from the filamentous anoxygenic phototroph Chloroflexus aurantiacus that display the entire known spectral and redox variations known in the type 1 copper site. C. aurantiacus encodes four auracyanins, labeled A-D. The midpoint potentials vary from +83 mV (auracyanin D) to +423 mV (auracyanin C). The electronic structures vary from classical blue copper UV-vis absorption spectra (auracyanin B) to highly perturbed spectra (auracyanins C and D). The spectrum of auracyanin C is temperature-dependent. The expansion and divergent nature of the auracyanins is a previously unseen phenomenon.

Computational study of the structure and electronic circular dichroism spectroscopy of blue copper proteins

The calculation of the electronic circular dichroism (CD) spectra of the oxidized form of the blue copper proteins plastocyanin and cucumber basic protein and the relationship between the observed spectral features and the structure of the active site of the protein is investigated. Excitation energies and transition strengths are computed using multireference configuration interaction, and it is shown that computed spectra based on coordinates from the crystal structure or a single structure optimized in quantum mechanics/molecular mechanics (QM/MM) or ligand field molecular mechanics (LFMM) are qualitatively incorrect. In particular, the rotational strength of the ligand to metal charge transfer band is predicted to be too small or have the incorrect sign. By considering calculations on active site models with modified structures, it is shown that the intensity of this band is sensitive to the nonplanarity of the histidine and cysteine ligands coordinated to copper. Calculation of the ultraviolet absorption and CD spectra based upon averaging over many structures drawn from a LFMM molecular dynamics simulation are in good agreement with experiment, and superior to analogous calculations based upon structures from a classical molecular dynamics simulation. This provides evidence that the LFMM force field provides an accurate description of the molecular dynamics of these proteins.

Comparative analysis of the performance of commonly available density functionals in the determination of geometrical parameters for copper complexes

In this study, a set of 50 transition-metal complexes of Cu(I) and Cu(II), were used in the evaluation of 18 density functionals in geometry determination. In addition, 14 different basis sets were considered, including four commonly used Pople's all-electron basis sets; four basis sets including popular types of effective-core potentials: Los Alamos, Steven-Basch-Krauss, and Stuttgart-Dresden; and six triple-ζ basis sets. The results illustrate the performance of different methodological alternatives for the treatment of geometrical properties in relevant copper complexes, pointing out Double-Hybrid (DH) and Long-range Correction (LC) Generalized Gradient Approximation (GGA) methods as better descriptors of the geometry of the evaluated systems. These however, are associated with a computational cost several times higher than some of the other methods employed, such as the M06 functional, which has also demonstrated a comparable performance. Regarding the basis sets, 6-31+G(d) and 6-31+G(d,p) were the best performing approaches. In addition, the results show that the use of effective-core potentials has a limited impact, in terms of the accuracy in the determination of metal-ligand bond-lengths and angles in our dataset of copper complexes. Hence, these could become a good alternative for the geometrical description of these systems, particularly CEP-121G and SDD basis sets, if one is considering larger copper complexes where the computational cost could be an issue.

Analysis of resonance Raman data on the blue copper site in pseudoazurin: excited state π and σ charge transfer distortions and their relation to ground state reorganization energy

The short Cu(2+)-S(Met) bond in pseudoazurin (PAz) results in the presence of two relatively intense S(p)(π) and S(p)(σ) charge transfer (CT) transitions. This has enabled resonance Raman (rR) data to be obtained for each excited state. The rR data show very different intensity distribution patterns for the vibrations in the 300-500 cm(-1) region. Time-dependent density functional theory (TDDFT) calculations have been used to determine that the change in intensity distribution between the S(p)(π) and S(p)(σ) excited states reflects the differential enhancement of S(Cys) backbone modes with Cu-S(Cys)-C(β) out-of-plane (oop) and in-plane (ip) bend character in their respective potential energy distributions (PEDs). The rR excited state distortions have been related to ground state reorganization energies (λ s) and predict that, in addition to M-L stretches, the Cu-S(Cys)-C(β) oop bend needs to be considered. DFT calculations predict a large distortion in the Cu-S(Cys)-C(β) oop bending coordinate upon reduction of a blue copper (BC) site; however, this distortion is not present in the X-ray crystal structures of reduced BC sites. The lack of Cu-S(Cys)-C(β) oop distortion upon reduction corresponds to a previously unconsidered constraint on the thiolate ligand orientation in the reduced state of BC proteins and can be considered as a contribution to the entatic/rack nature of BC sites.

Cytotoxic origin of copper(II) oxide nanoparticles: comparative studies with micron-sized particles, leachate, and metal salts

The work investigates the source of toxicity of copper oxide nanoparticles (CuO NPs) with respect to its leaching characteristic and speciation. Complexation-mediated leaching of CuO NPs by amino acids was identified as the source of toxicity toward Escherichia coli, the model microorganism used in the current study. The leached copper-peptide complex induces a multiple-fold increase in intracellular reactive oxygen species generation and reduces the fractions of viable cells, resulting in the overall inhibition of biomass growth. The cytotoxicity of the complex leachate is however different from that of equivalent soluble copper salts (nitrates and sulfates). A pH-dependent copper speciation during the addition of copper salts gives rise to uncoordinated copper ions, which in turn result in greater toxicity and cell lysis, the latter of which was not observed for CuO NPs even at comparable pH. Since leaching did not occur with micrometer-sized CuO, no cytotoxicty effect was observed, thus highlighting the prominence of materials toxicity at the nanoscale.

Folding and unfolding in the blue copper protein rusticyanin: role of the oxidation state

The unfolding process of the blue copper protein rusticyanin has been studied from the structural and the thermodynamic points of view at two pH values (pH 2.5 and 7.0). When Rc unfolds, copper ion remains bound to the polypeptide chain. Nuclear magnetic resonance data suggest that three of the copper ligands in the folded state are bound to the metal ion in the unfolded form, while the other native ligand is detached. These structural changes are reflected in the redox potentials of the protein in both folded and unfolded forms. The affinities of the copper ion in both redox states have been also determined at the two specified pH values. The results indicate that the presence of two histidine ligands in the folded protein can compensate the change in the net charge that the copper ion receives from their ligands, while, in the unfolded protein, charges of aminoacids are completely transferred to the copper ion, altering decisively the relative stability of its two-redox states.

Thermodynamic equilibrium between blue and green copper sites and the role of the protein in controlling function

A combination of spectroscopies and density functional theory calculations indicate that there are large temperature-dependent absorption spectral changes present in green nitrite reductases (NiRs) due to a thermodynamic equilibrium between a green and a blue type 1 (T1) copper site. The axial methionine (Met) ligand is unconstrained in the oxidized NiRs, which results in an enthalpically favored (DeltaH approximately 4.6 kcal/mol) Met-bound green copper site at low temperatures, and an entropically favored (TDeltaS approximately 4.5 kcal/mol, at room temperature) Met-elongated blue copper site at elevated temperatures. In contrast to the NiRs, the classic blue copper sites in plastocyanin and azurin show no temperature-dependent behavior, indicating that a single species is present at all temperatures. For these blue copper proteins, the polypeptide matrix opposes the gain in entropy that would be associated with the loss of the weak axial Met ligand at physiological temperatures by constraining its coordination to copper. The potential energy surfaces of Met binding indicate that it stabilizes the oxidized state more than the reduced state. This provides a mechanism to tune down the reduction potential of blue copper sites by >200 mV.

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.

Periodic trends within a series of five-coordinate thiolate-ligated [MII(SMe2N4(tren))]+ (M = Mn, Fe, Co, Ni, Cu, Zn) complexes, including a rare example of a stable CuII-thiolate

A series of five-coordinate thiolate-ligated complexes [M(II)(tren)N4S(Me2)]+ (M = Mn, Fe, Co, Ni, Cu, Zn; tren = tris(2-aminoethyl)amine) are reported, and their structural, electronic, and magnetic properties are compared. Isolation of dimeric [Ni(II)(SN4(tren)-RS(dang))]2 ("dang"= dangling, uncoordinated thiolate supported by H bonds), using the less bulky [(tren)N4S](1-) ligand, pointed to the need for gem-dimethyls adjacent to the sulfur to sterically prevent dimerization. All of the gem-dimethyl derivatized complexes are monomeric and, with the exception of [Ni(II)(S(Me2)N4(tren)]+, are isostructural and adopt a tetragonally distorted trigonal bipyramidal geometry favored by ligand constraints. The nickel complex uniquely adopts an approximately ideal square pyramidal geometry and resembles the active site of Ni-superoxide dismutase (Ni-SOD). Even in coordinating solvents such as MeCN, only five-coordinate structures are observed. The MII-S thiolate bonds systematically decrease in length across the series (Mn-S > Fe-S > Co-S > Ni-S approximately Cu-S < Zn-S) with exceptions occurring upon the occupation of sigma* orbitals. The copper complex, [Cu(II)(S(Me2)N4(tren)]+, represents a rare example of a stable CuII-thiolate, and models the perturbed "green" copper site of nitrite reductase. In contrast to the intensely colored, low-spin Fe(III)-thiolates, the M(II)-thiolates described herein are colorless to moderately colored and high-spin (in cases where more than one spin-state is possible), reflecting the poorer energy match between the metal d- and sulfur orbitals upon reduction of the metal ion. As the d-orbitals drop in energy proceeding across the across the series M(2+) (M= Mn, Fe, Co, Ni, Cu), the sulfur-to-metal charge-transfer transition moves into the visible region, and the redox potentials cathodically shift. The reduced M(+1) oxidation state is only accessible with copper, and the more oxidized M(+4) oxidation state is only accessible for manganese.

Spectroscopic characterization of a high-potential lipo-cupredoxin found in Streptomyces coelicolor

For many streptomycetes, a distinct dependence on the "bioavailability" of copper ions for their morphological development has been reported. Analysis of the Streptomyces coelicolor genome reveals a number of gene products encoding for putative copper-binding proteins. One of these appears as an unusual copper-binding protein with a lipoprotein signal sequence and a cupredoxin-like domain harboring a putative Type-1 copper-binding motif. Cloning of this gene from S. coelicolor and subsequent heterologous expression in Escherichia coli has allowed for a thorough spectroscopic interrogation of this putative copper-binding protein. Optical and electron paramagnetic resonance spectroscopies have confirmed the presence of a "classic" Type-1 copper site with the axial ligand to the copper a methionine. Paramagnetic NMR spectroscopy on both the native Cu(II) form and Co(II)-substituted protein has yielded active-site structural information, which on comparison with that of other cupredoxin active sites reveals metal-ligand interactions most similar to the "classic" Type-1 copper site found in the amicyanin family of cupredoxins. Despite this high structural similarity, the Cu(II)/(I) midpoint potential of the S. coelicolor protein is an unprecedented +605 mV vs normal hydrogen electrode at neutral pH (amicyanin approximately +250 mV), with no active-site protonation of the N-terminal His ligand observed. Suggestions for the physiological role/function of this high-potential cupredoxin are discussed.

Comparison of intra- vs intermolecular long-range electron transfer in crystals of ruthenium-modified azurin

Selective metal-ion incorporation and ligand substitution are employed to control whether electrons tunnel over intra- or intermolecular separations in crystals of P. aeruginosa azurin modified with Ru-polypyridine complexes. Cu(1+)-to-Ru3+ electron transfer (ET) across a specific protein-protein interface in the crystal lattice has a time constant 5-10 times longer than ET between the same donor and acceptor within a single protein (tauET = 5 vs 0.5-1.0 micros). Slower intermolecular ET agrees well with a longer distance between redox centers across the inter-protein (18.9 A) compared to the intra-protein separation (17.0 A) and indicates that the closest donor/acceptor pair dominates crystal ET. Lowering the crystal pH accelerates inter-protein ET (tauET = 1.0 micros) but not intra-protein ET. Faster inter-protein ET likely results from a pH-induced peptide bond flip that perturbs hydrogen bonding in the path between Ru and Cu centers on adjacent molecules.

Active site comparison of CoII blue and green nitrite reductases

Copper-containing nitrite reductases (NiRs) possess type 1 (T1) and type 2 (T2) copper sites and can be either green or blue in color owing to differences at their T1 centers. The active sites of a green and a blue NiR were studied by utilizing their T1CuI/T2CoII and T1CoII/T2CoII-substituted forms. The UV/Vis spectra of these derivatives highlight the similarity of the T2 centers in these enzymes and that T1 site differences are also present in the CoII forms. The paramagnetic NMR spectra of T1CuI/T2CoII enzymes allow hyperfine shifted resonances from the three T2 His ligands to be assigned: these exhibit remarkably similar positions in the spectra of both NiRs, emphasizing the homology of the T2 centers. The addition of nitrite results in subtle alterations in the paramagnetic NMR spectra of the T1CuI/T2CoII forms at pH<7, which indicate a geometry change upon the binding of substrate. Shifted resonances from all of the T1 site ligands have been assigned and the CoII--N(His) interactions are alike, whereas the CbetaH proton resonances of the Cys ligand exhibit subtle chemical shift differences in the blue and green NiRs. The strength of the axial CoII--S(Met) interaction is similar in the two NiRs studied, but the altered conformation of the side chain of this ligand results in a dramatically different chemical shift pattern for the CgammaH protons. This indicates an alteration in the bonding of the axial ligand in these derivatives, which could be influential in the CuII proteins.

Comprehensive Molecular Mechanics Model for Oxidized Type I Copper Proteins: Active Site Structures, Strain Energies, and Entatic Bulging

The ligand field molecular mechanics (LFMM) model has been applied to the oxidized Type 1 copper center. In conjunction with the AMBER94 force field implemented in DommiMOE, the ligand field extension of the molecular operating environment (MOE), LFMM parameters for Cu-N(imidazole), Cu-S(thiolate), Cu-S(thioether), and Cu-O(carbonyl) interactions were developed on the basis of experimental and theoretical data for homoleptic model systems. Subsequent LFMM optimizations of the active site model complex [Cu(imidazole)2(SMe)(SMe2]+ agree with high level quantum results both structurally and energetically. Stable trigonal and tetragonal structures are located with the latter about 1.5 kcal mol-1 lower in energy. Fully optimized unconstrained structures were computed for 24 complete proteins containing T1 centers spanning four-coordinate, plastocyanin-like CuN2SS' and stellacyanin-like CuN2SO sites, plus the five-coordinate CuN2SS'O sites of the azurins. The initial structures were based on PDB coordinates augmented by a 10 A layer of water molecules. Agreement between theory and experiment is well within the experimental uncertainties. Moreover, the LFMM results for plastocyanin (Pc), cucumber basic protein (CBP) and azurin (Az) are at least as good as previously reported QM/MM structures and are achieved several orders of magnitude faster. The LFMM calculations suggest the protein provides an entatic strain of about 10 kcal mol-1. However, when combined with the intrinsic 'plasticity' of d9 Cu(II), different starting protein/solvent configurations can have a significant effect on the final optimized structure. This 'entatic bulging' results in relatively large fluctuations in the calculated metal-ligand bond lengths. For example, simply on the basis of 25 different starting configurations of the solvent molecules, the optimized Cu-S(thiolate) bond lengths in Pc vary by 0.04 A while the Cu-S(thioether) distance spans over 0.3 A. These variations are the same order of magnitude as the differences often quoted to correlate the spectroscopic properties from a set of proteins. Isolated optimizations starting from PDB coordinates (or indeed, the PDB structures themselves) may only accidentally correlate with spectroscopic measurements. The present calculations support the work of Warshel who contends that adequate configurational averaging is necessary to make proper contact with experimental properties measured in solution. The LFMM is both sufficiently accurate and fast to make this practical.

Reduction potential tuning of the blue copper center in Pseudomonas aeruginosa azurin by the axial methionine as probed by unnatural amino acids

The conserved axial ligand methionine 121 from Pseudomonas aeruginosa azurin (Az) has been replaced by isostructural unnatural amino acid analogues, oxomethionine (OxM), difluoromethionine (DFM), trifluoromethionine (TFM), selenomethionine (SeM), and norleucine (Nle) using expressed protein ligation. The replacements resulted in < 6 nm shifts in the S(Cys)-Cu charge transfer (CT) band in the electronic absorption spectra and < 8 gauss changes in the copper hyperfine coupling constants (AII) in the X-band electron paramagnetic resonance spectra, suggesting that isostructural replacement of Met resulted in minimal structural perturbation of the copper center. The slight blue shifts of the CT band follow the trend of stronger electronegativity of the ligands. This trend is supported by 19F NMR studies of the fluorinated methionine analogues. However, the order of AII differs, suggesting additional factors influencing AII. In contrast to the small changes in the UV-vis and EPR spectra, a large variation of > 227 mV in reduction potential was observed for the series of variants reported here. Additionally, a linear correlation was established between the reduction potentials and hydrophobicity of the variants. Extension of this analysis to other type 1 copper-containing proteins reveals a linear correlation between change in hydrophobicity and change in reduction potential, independent of the protein scaffold, experimental conditions, measurement techniques, and steric modifications. This analysis has also revealed for the first time high and low potential states for type 1 centers, and the difference may be attributable to destabilization of the protein fold by disruption of hydrophobic or hydrogen bonding interactions that stabilize the type 1 center.