On the relative stability of tetragonal and trigonal Cu(II) complexes with relevance to the blue copper proteins

 The role of the cysteine thiolate ligand for the unusual copper coordination geometry in the blue copper proteins has been studied by comparing the electronic structure, geometry, and energetics of a number of small Cu(II) complexes. The geometries have been optimised with the density functional B3LYP method, and energies have been calculated by multiconfigurational second-order perturbation theory (the CASPT2 method). Most small inorganic Cu(II) complexes assume a tetragonal geometry, where four ligands make σ bonds to a Cu 3d orbital. If a ligand lone-pair orbital instead forms a π bond to the copper ion, it formally occupies two ligand positions in a square coordination, and the structure becomes trigonal. Large, soft, and polarisable ligands, such as SH^– and SeH^–, give rise to covalent copper-ligand bonds and structures close to a tetrahedron, which might be trigonal or tetragonal with approximately the same stability. On the other hand, small and hard ligands, such as NH₃, OH₂, and OH^–, give ionic bonds and flattened tetragonal structures. It is shown that axial type 1 (blue) copper proteins have a trigonal structure with a π bond to the cysteine sulphur atom, whereas rhombic type 1 and type 2 proteins have a tetragonal structure with σ bonds to all strong ligands. The soft cysteine ligand is essential for the stabilisation of a structure that is close to a tetrahedron (either trigonal or tetragonal), which ensures a low reorganisation energy during electron transfer. Received: 9 July 1997 / 26 November 1997     

Control of metalloprotein reduction potential: compensation phenomena in the reduction thermodynamics of blue copper proteins

The reduction thermodynamics (Delta H degrees '(rc) and Delta S degrees '(rc)) for native Paracoccus versutus amicyanin, for Alcaligenes faecalis S-6 pseudoazurin, and for the G45P, M64E, and K27C variants of Pseudomonas aeruginosa azurin were measured electrochemically. Comparison with the data available for other native and mutated blue copper proteins indicates that the features of metal coordination and the electrostatic potential due to the protein matrix and the solvent control the reduction enthalpy in a straightforward way. However, the effects on the reduction potential are rather unpredictable owing to the entropic contribution to E degrees ', which is mainly determined by solvent reorganization effects. Analysis of all the Delta H degrees '(rc) and Delta S degrees '(rc) values available for this protein class indicates that enthalpy-entropy compensation occurs in the reduction thermodynamics of wt cupredoxins from different sources, as well as for mutants of the same species. The findings indicate that the reduction enthalpies and entropies for these species are strongly affected by reduction-induced reorganization of solvent molecules within the solvation sphere of the protein. The absence of a perfect enthalpy-entropy compensation is due to the fact that while the differences between reduction entropies are dominated by solvent reorganization effects, those between reduction enthalpies are significantly controlled by intrinsic molecular factors related to the selective stabilization of the reduced form by coordination features of the copper site and electrostatic effects at the interface with the protein matrix.     

Structural control of electron-transfer properties in metalloproteins

Summary The factors that control long-range electron transfer between two redox centers in a protein are summarized. Rack-induced bonding in blue copper proteins is described. The protein conformation forces the Cu(II) ion into a distorted geometry, lying at least 70 kJ mol^–1 above the preferred square-planar geometry in energy. The distortion has the effect that the structural change associated with electron transfer is minimal and thus the reorganization energy small. Variations in back bonding are suggested to modulate the reduction potentials of blue proteins without any change in the energy of the charge-transfer transitions. In proton pumps there must be a structural control of the electron transfer rates (electron gating) and model studies suggest that this is best achieved by variations in the reorganization energy.     

Protein strain in blue copper proteins studied by free energy perturbations.

Free energy perturbations have been performed on two blue copper proteins, plastocyanin and nitrite reductase. By changing the copper coordination geometry, force constants, and charges, we have estimated the maximum energy with which the proteins may distort the copper coordination sphere. By comparing this energy with the quantum chemical energy cost for the same perturbation on the isolated copper complex, various hypotheses about protein strain have been tested. The calculations show that the protein can only modify the copper-methionine bond length by a modest amount of energy-<5 kJ/mol-and they lend no support to the suggestion that the quite appreciable difference in the copper coordination geometry encountered in the two proteins is a result of the proteins enforcing different Cu-methionine bond lengths. On the contrary, this bond is very flexible, and neither the geometry nor the electronic structure change appreciably when the bond length is changed. Moreover, the proteins are rather indifferent to the length of this bond. Instead, the Cu(II) coordination geometries in the two proteins represent two distinct minima on the potential surface of the copper ligand sphere, characterized by different electronic structures, a tetragonal, mainly sigma-bonded, structure in nitrite reductase and a trigonal, pi-bonded, structure in plastocyanin. In vacuum, the structures have almost the same energy, and they are stabilized in the proteins by a combination of geometric and electrostatic interactions. Plastocyanin favors the bond lengths and electrostatics of the trigonal structure, whereas in nitrite reductase, the angles are the main discriminating factor. Proteins 1999;36:157-174.     

An outer-sphere hydrogen-bond network constrains copper coordination in blue proteins

In azurins and other blue copper proteins with relatively low reduction potentials (E(0) [Cu(II)/Cu(I)]<400 mV vs. normal hydrogen electrode), the folded polypeptide framework constrains both copper(II) and copper(I) in such a way as to tune the reduction potentials to values that differ greatly from those for most copper complexes. Largely conserved networks of hydrogen bonds organize and lock the rest of the folded protein structure to a loop that contains three of the ligands to copper. Changes in hydrogen bonds that allow copper(I) to revert more closely to its preferred geometry [relative to the copper(II) geometry] accordingly lead to an increase in E(0). This paper reports mutations in the ligand loop of amicyanin from P. denitrificans that relax the constraints on ligation for copper(I) and significantly raise E(0) for these mutants (for example 415+/-4 mV) relative to that of the native amicyanin (265+/-4 mV). These mutations also shift the pK(a) of a ligand histidine to below 5 relative to 7.0 in the wild type.     

Energy saving electron pathways in proteins

This paper is a contribution to the discussion of whether the general architecture of electron transfer sites in blue copper proteins is mainly a result of the structural preferences of the metal ion or is induced by the protein. Although the site is probably stable only when protected by the protein, there appears to be no strain from the latter on the structure in the vicinity of the copper atom. For an operative redox site it is further required that the geometry of the site is acceptable for both oxidation states, to avoid high reorganization energy. The site must also be connected to the outer world by suitable tunneling pathways. The blue copper sites appear to fulfill these requirements, but it is difficult to assess the role of evolutionary pressure to form electron transfer proteins in general.     

Quantum chemical calculations of the reorganization energy of blue-copper proteins.

The inner-sphere reorganization energy for several copper complexes related to the active site in blue-copper protein has been calculated with the density functional B3LYP method. The best model of the blue-copper proteins, Cu(Im)2(SCH3)(S(CH3)2)(0/+), has a self-exchange inner-sphere reorganization energy of 62 kJ/mol, which is at least 120 kJ/mol lower than for Cu(H2O)4(+/2+). This lowering of the reorganization energy is caused by the soft ligands in the blue-copper site, especially the cysteine thiolate and the methionine thioether groups. Soft ligands both make the potential surfaces of the complexes flatter and give rise to oxidized structures that are quite close to a tetrahedron (rather than tetragonal). Approximately half of the reorganization energy originates from changes in the copper-ligand bond lengths and half of this contribution comes from the Cu-S(Cys) bond. A tetragonal site, which is present in the rhombic type 1 blue-copper proteins, has a slightly higher (16 kJ/mol) inner-sphere reorganization energy than a trigonal site, present in the axial type 1 copper proteins. A site with the methionine ligand replaced by an amide group, as in stellacyanin, has an even higher reorganization energy, about 90 kJ/mol.     

The use of osmolytes to facilitate protein NMR spectroscopy

Summary A method of stabilizing folded proteins is described, which allows NMR studies under conditions where a protein would normally be unfolded. This enables stable proteins to be examined at elevated temperatures, or spectra recorded on samples that are insufficiently stable under normal conditions. Up to two molar perdeuterated glycine, a potent osmolyte, can be added to aqueous protein NMR samples without altering the folded three-dimensional structure or function of the protein. However, the stability of the folded form is dramatically increased. This is illustrated for the protein lysozyme at high temperature (348 K) where the structural integrity is destroyed in standard aqueous solution, but is retained in the osmolyte solution. We hope that the technique will be of value to those studying by NMR the structural biology of protein fragments and mutants, which are often of reduced stability compared with the original proteins.To whom correspondence should be addressed.     

Catalysis and electron transfer in protein crystals: the binary and ternary complexes of methylamine dehydrogenase with electron acceptors

Polarized absorption microspectrophotometry has been used to detect catalysis and intermolecular electron transfer in single crystals of two multiprotein complexes: (1) the binary complex between Paracoccus denitrificans methylamine dehydrogenase, which contains tryptophan-tryptophylquinone (TTQ) as a cofactor, and its redox partner, the blue copper protein amicyanin; (2) the ternary complex between the same two proteins and cytochrome c-551i. Continuous wave electron paramagnetic resonance has been used to compare the state of copper in polycrystalline powders of the two systems. While catalysis and intermolecular electron transfer from reduced TTQ to copper are too fast to be accessible to our measurements, heme reduction occurs over a period of several minutes. The observed rate constant is about four orders of magnitude lower than in solution. The analysis of the temperature dependence of this apparent constant provides values for the parameters H(AB), related to electronic coupling between the two centers, and lambda, the reorganizational energy, that are compatible with electron transfer being the rate-determining step. From these parameters and the known distance between copper and heme, it is possible to calculate the parameter beta, which depends on the nature of the intervening medium, obtaining a value typical of electron transfer across a protein matrix. These findings suggest that the ternary complex in solution might achieve a higher efficiency than the rigid crystal structure thanks to an as yet unidentified role of protein dynamics.     

The medium reorganization energy for the charge transfer reactions in proteins

A low static dielectric permittivity of proteins causes the low reorganization energies for the charge transfer reactions inside them. This reorganization energy does not depend on the pre-existing intraprotein electric field. The charge transferred inside the protein interacts with its aqueous surroundings; for many globular proteins, the effect of this surroundings on the reorganization energy is comparable with the effect of reorganization of the protein itself while for the charge transfer in the middle of membrane the aqueous phase plays a minor role. Reorganization energy depends strongly on the system considered, and hence there is no sense to speak on the "protein reorganization energy" as some permanent characteristic parameter. We employed a simple algorithm for calculation of the medium reorganization energy using the numerical solution of the Poisson-Boltzmann equation. Namely, the reaction field energy was computed in two versions - all media having optical dielectric permittivity, and all the media with the static one; the difference of these two quantities gives the reorganization energy. We have calculated reorganization energies for electron transfer in cytochrome c, various ammine-ruthenated cytochromes c, azurin, ferredoxin, cytochrome c oxidase, complex of methylamine dehydrogenase with amicyanin, and for proton transfer in α-chymotrypsin. It is shown that calculation of the medium reorganization energy can be a useful tool in analysis of the mechanisms of the charge transfer reactions in proteins.     

Reduction potentials of blue and purple copper proteins in their unfolded states: a closer look at rack-induced coordination

 Cyclic voltammetry has been used to determine the reduction potentials of blue (Pseudomonas aeruginosa azurin) and purple (Thermus thermophilus Cu_A domain) copper proteins unfolded by guanidine hydrochloride. These Cu(II/I) potentials [456 (azurin); 453 (Cu_A) mV vs., NHE] are higher than those of the folded proteins. The downshift of the potential in the folded state can be accounted for by assuming that rack-induced axial coordination stabilizes Cu(II) relative to Cu(I) in a protein-encapsulated active site. Received: 3 March 1998 / Accepted: 6 April 1998     

Reduction potential variations in azurin through secondary coordination sphere phenylalanine incorporations

Recent evidence has shown that the properties of metal binding sites can be tuned by more than the ligands in the primary coordination sphere. We investigated the incorporation of four phenylalanine residues into the secondary coordination sphere of the small soluble blue copper protein azurin. The locations for placement of these residues in azurin were based on the structure of the highly hydrophobic blue copper protein rusticyanin, which is known to have a significantly higher reduction potential than azurin. Using site-directed mutagenesis, these residues in close proximity to the copper binding site were mutated to large hydrophobic phenylalanine residues individually and in combination. We also added the Met121Leu mutation on top of the Phe mutations to construct a total of 13 variants. We found little change in the UV-visible absorption and EPR data for these proteins, however modest increases in reduction potential were observed with increases by as much as 30 mV per Phe residue. Furthermore, we observed the increases in potential to be additive.     

Crystal structure of a novel red copper protein from Nitrosomonas europaea

Nitrosocyanin (NC) is a mononuclear red copper protein isolated from the ammonia oxidizing bacterium Nitrosomonas europaea. Although NC exhibits some sequence homology to classic blue copper proteins, its spectroscopic and electrochemical properties are drastically different. The 1.65 A resolution crystal structure of oxidized NC reveals an unprecedented trimer of single domain cupredoxins. Each copper center is partially covered by an unusual extended beta-hairpin structure from an adjacent monomer. The copper ion is coordinated by His 98, His 103, Cys 95, a single side chain oxygen of Glu 60, and a solvent molecule. In the 2.3 A resolution structure of reduced NC, His 98 shifts away from the copper ion, and the solvent molecule is not observed. The arrangement of these ligands renders the coordination geometry of the NC red copper center distinct from that of blue copper centers. In particular, the red copper center has a higher coordination number and lacks the long Cu-S(Met) and short Cu-S(Cys) bond distances characteristic of blue copper. Moreover, the red copper center is square pyramidal whereas blue copper is typically distorted tetrahedral. Analysis of the NC structure provides insight into possible functions of this new type of biological copper center.     

Partly folded states of bovine carbonic anhydrase interact with zwitterionic and anionic lipid membranes

The acidic, partly folded states of bovine carbonic anhydrase II (BCAII) were used as an experimental system to study the interactions of partly denatured proteins with lipid membranes. The pH dependence of their interactions with palmitoyloleoyl phosphatidylcholine (POPC) and palmitoyloleoyl phosphatidylglycerol (POPG) membranes was studied. A filtration binding assay shows that acidic partly folded states of BCAII bind to POPC membranes. Fluorescence emission spectra from Trp residues of the bound protein are slightly shifted to shorter wavelength and can be quenched by a water-soluble quencher of fluorescence, indicating that the binding occurs without deep penetration of Trp residues into the membrane. The content of beta-structures of the protein in solution, as revealed by FT-IR spectroscopy, decreases in the partly folded states and the binding to POPC membrane occurs without further changes of secondary structure. In the presence of 0.1 M NaCl, a partly folded state self-aggregates and does not bind to POPC membrane. At acidic pH, BCAII binds to POPG membranes both at high and low ionic strength. The binding to the anionic lipid occurs with protein self-aggregation within the lipid-protein complexes and with changes in the secondary structure; large blue shifts in the fluorescence emission spectra and the decrease in the exposure to water-soluble acrylamide quencher of Trp fluorescence strongly suggest that BCAII penetrates the hydrocarbon domain in the POPG-protein complexes.     

Stellacyanin. Studies of the metal-binding site using x-ray absorption spectroscopy.

Stellacyanin is a mucoprotein of molecular weight approximately 20,000 containing one copper atom in a blue or type I site. The metal ion can exist in both the Cu(II) and Cu(I) redox states. The metal binding site in plastocyanin, another blue copper protein, contains one cysteinyl, one methionyl, and two imidazoyl residues (Colman et al. 1978. Nature [Lond.]. 272:319-324.), but an exactly analogous site cannot exist in stellacyanin as it lacks methionine. The copper coordination in stellacyanin has been studied by x-ray edge absorption and extended x-ray absorption fine structure (EXAFS) analysis. A new, very conservative data analysis procedure has been introduced, which suggests that the there are two nitrogen atoms in the first coordination shell of the oxidized [Cu(II)] protein and one in the reduced [Cu(I)] protein; these N atoms have normal Cu--N distances: 1.95-2.05 A. In both redox states there are either one or two sulfur atoms coordinating the copper, the exact number being indeterminable from the present data. In the oxidized state the Cu--S distance is intermediate between the short bond found in plastocyanin and those found in near tetragonal copper model compounds. Above -140 degree C, radiation damage of the protein occurs. At room temperature the oxidized proteins is modified in the x-ray beam at a rate of 0.25%/s.     

Conformation of peptide fragments of proteins in aqueous solution: implications for initiation of protein folding

Applications of sensitive new technologies, in particular, two-dimensional NMR spectroscopy, have allowed detection of folded structures in short peptide fragments of proteins in aqueous solution under conditions where native proteins fold. These structures are in rapid dynamic exchange with unfolded states. These observations provide evidence in support of models for protein folding which postulate localized regions of folded structure as initiation sites for the folding process. Since these initiation processes are expected to be rapid, such models are consistent with kinetic evidence that the rate-determining steps of protein folding occur late in the process and probably involve rearrangement of incorrectly folded intermediates.     

Electron spin relaxation of CuA and cytochrome a in cytochrome c oxidase. Comparison to heme, copper, and sulfur radical complexes

The method of continuous saturation has been used to measure the electron spin relaxation parameter T1T2 at temperatures between 10 and 50 K for a variety of S = 1/2 species including: CuA and cytochrome a of cytochrome c oxidase, the type 1 copper in several blue copper proteins, the type 2 copper in laccase, inorganic Cu(II) complexes, sulfur radicals, and low spin heme proteins. The temperature dependence and the magnitude of T1T2 for all of the species examined are accounted for by assuming that the Van Vleck Raman process dominates the electron spin-lattice relaxation. Over the entire temperature range examined, the relaxation of the type 1 coppers in six to seven times faster than that of type 2 copper, inorganic copper, and sulfur radicals, in spite of the similar g-anisotropies of these species. This result may indicate that the coupling of the phonon bath to the spin center is more effective in type 1 coppers than in the other complexes studied. The relaxation of CuA of cytochrome oxidase exhibits an unusual temperature dependence relative to the other copper complexes studied, suggesting that the protein environment of this center is different from that of the other copper centers studied and/or that CuA is influenced by a magnetic dipolar interaction with another, faster-relaxing paramagnetic site in the enzyme. A comparison of the saturation characteristics of the CuA EPR signal in native and partially reduced CO complexes of the enzyme also suggests the existence of such an interaction. The implications of these results with respect to the disposition of the metal centers in cytochrome oxidase are discussed.     

Conformational properties of azurin in solution as determined from resolution-enhanced Fourier-transform infrared spectra

Infrared spectra of the blue copper protein azurin and of apoazurin from P. fluorescens were obtained in aqueous solution. Using resolution enhancement procedures, a number of component bands were identified in the region of the amide I mode, and these bands were assigned to various components of protein secondary structure. A quantitative analysis of these infrared spectra indicates that the secondary structure of P. fluorescens azurin in solution is very similar to those determined previously by X-ray diffraction for the crystals of azurins from other bacterial species. The major components of this structure are beta strands and turns. Infrared spectra also evince a remarkable thermal stability of the native azurin. A significant unfolding of the protein could only be detected at temperatures above approximately 76 degrees C. While the secondary structure of apoazurin is practically indistinguishable from that of the native protein at room temperature, the thermal stability of the apo form is significantly reduced.     

Residual structure in disordered peptides and unfolded proteins from multivariate analysis and ab initio simulation of Raman optical activity data

Vibrational Raman optical activity (ROA), measured as a small difference in the intensity of Raman scattering from chiral molecules in right- and left-circularly polarized incident light, or as the intensity of a small circularly polarized component in the scattered light, is a powerful probe of the aqueous solution structure of proteins. The large number of structure-sensitive bands in protein ROA spectra makes multivariate analysis techniques such as nonlinear mapping (NLM) especially favorable for determining structural relationships between different proteins. We have previously used NLM to map a large dataset of peptide, protein, and virus ROA spectra into a readily visualizable two-dimensional space in which points close to or distant from each other, respectively, represent similar or dissimilar structures. As well as folded proteins, our dataset contains ROA spectra from many natively unfolded proteins, proteins containing both folded and unfolded domains, denatured partially structured molten globule and reduced protein states, together with folded proteins containing little or no alpha-helix or beta-sheet. In this article, the relative positions of these systems in the NLM plot are used to obtain information about any residual structure that they may contain. The striking differences between the structural propensities of proteins that are unfolded in their native states and those that are unfolded due to denaturation may be responsible for their often very different behavior, especially with regard to aggregation. An ab initio simulation of the Raman and ROA spectra of an alanine oligopeptide in the poly(L-proline) II-helical conformation confirms previous suggestions that this conformation is a significant structural element in disordered peptides and natively unfolded proteins. The use of ROA to identify and characterize proteins containing significant amounts of unfolded structure will, inter alia, be valuable in structural genomics/proteomics since unfolded sequences often inhibit crystallization.     

Comparative analysis of the human and chicken prion protein copper binding regions at pH 6.5

Recent experimental evidence supports the hypothesis that prion proteins (PrPs) are involved in the Cu(II) metabolism. Moreover, the copper binding region has been implicated in transmissible spongiform encephalopathies, which are caused by the infectious isoform of prion proteins (PrP(Sc)). In contrast to mammalian PrP, avian prion proteins have a considerably different N-terminal copper binding region and, most interestingly, are not able to undergo the conversion process into an infectious isoform. Therefore, we applied x-ray absorption spectroscopy to analyze in detail the Cu(II) geometry of selected synthetic human PrP Cu(II) octapeptide complexes in comparison with the corresponding chicken PrP hexapeptide complexes at pH 6.5, which mimics the conditions in the endocytic compartments of neuronal cells. Our results revealed that structure and coordination of the human PrP copper binding sites are highly conserved in the pH 6.5-7.4 range, indicating that the reported pH dependence of copper binding to PrP becomes significant at lower pH values. Furthermore, the different chicken PrP hexarepeat motifs display homologous Cu(II) coordination at sub-stoichiometric copper concentrations. Regarding the fully cation-saturated prion proteins, however, a reduced copper coordination capability is supposed for the chicken prion protein based on the observation that chicken PrP is not able to form an intra-repeat Cu(II) binding site. These results provide new insights into the prion protein structure-function relationship and the conversion process of PrP.