Crystal structure analysis of oxidized Pseudomonas aeruginosa azurin at pH 5.5 and pH 9.0. A pH-induced conformational transition involves a peptide bond flip

Inner- and outer-sphere metal coordination in blue copper proteins

Blue copper proteins (BCPs) comprise classic cases of Nature's profound control over the electronic structures and chemical reactivity of transition metal ions. Early studies of BCPs focused on their inner coordination spheres, that is, residues that directly coordinate Cu. Equally important are the electronic and geometric perturbations to these ligands provided by the outer coordination sphere. In this tribute to Hans Freeman, we review investigations that have advanced the understanding of how inner-sphere and outer-sphere coordination affects biological Cu properties.

Flexibility of the metal-binding region in apo-cupredoxins

Protein-mediated electron transfer is an essential event in many biochemical processes. Efficient electron transfer requires the reorganization energy of the redox event to be minimized, which is ensured by the presence of rigid donor and acceptor sites. Electron transfer copper sites are present in the ubiquitous cupredoxin fold, able to bind one or two copper ions. The low reorganization energy in these metal centers has been accounted for by assuming that the protein scaffold creates an entatic/rack-induced state, which gives rise to a rigid environment by means of a preformed metal chelating site. However, this notion is incompatible with the need for an exposed metal-binding site and protein-protein interactions enabling metallochaperone-mediated assembly of the copper site. Here we report an NMR study that reveals a high degree of structural heterogeneity in the metal-binding region of the nonmetallated Cu(A)-binding cupredoxin domain, arising from microsecond to second dynamics that are quenched upon metal binding. We also report similar dynamic features in apo-azurin, a paradigmatic blue copper protein, suggesting a general behavior. These findings reveal that the entatic/rack-induced state, governing the features of the metal center in the copper-loaded protein, does not require a preformed metal-binding site. Instead, metal binding is a major contributor to the rigidity of electron transfer copper centers. These results reconcile the seemingly contradictory requirements of a rigid, occluded center for electron transfer, and an accessible, dynamic site required for in vivo copper uptake.

Outer-sphere contributions to the electronic structure of type zero copper proteins

Bioinorganic canon states that active-site thiolate coordination promotes rapid electron transfer (ET) to and from type 1 copper proteins. In recent work, we have found that copper ET sites in proteins also can be constructed without thiolate ligation (called "type zero" sites). Here we report multifrequency electron paramagnetic resonance (EPR), magnetic circular dichroism (MCD), and nuclear magnetic resonance (NMR) spectroscopic data together with density functional theory (DFT) and spectroscopy-oriented configuration interaction (SORCI) calculations for type zero Pseudomonas aeruginosa azurin variants. Wild-type (type 1) and type zero copper centers experience virtually identical ligand fields. Moreover, O-donor covalency is enhanced in type zero centers relative that in the C112D (type 2) protein. At the same time, N-donor covalency is reduced in a similar fashion to type 1 centers. QM/MM and SORCI calculations show that the electronic structures of type zero and type 2 are intimately linked to the orientation and coordination mode of the carboxylate ligand, which in turn is influenced by outer-sphere hydrogen bonding.

Biological applications of hybrid quantum mechanics/molecular mechanics calculation

Since in most cases biological macromolecular systems including solvent water molecules are remarkably large, the computational costs of performing ab initio calculations for the entire structures are prohibitive. Accordingly, QM calculations that are jointed with MM calculations are crucial to evaluate the long-range electrostatic interactions, which significantly affect the electronic structures of biological macromolecules. A UNIX-shell-based interface program connecting the quantum mechanics (QMs) and molecular mechanics (MMs) calculation engines, GAMESS and AMBER, was developed in our lab. The system was applied to a metalloenzyme, azurin, and PU.1-DNA complex; thereby, the significance of the environmental effects on the electronic structures of the site of interest was elucidated. Subsequently, hybrid QM/MM molecular dynamics (MD) simulation using the calculation system was employed for investigation of mechanisms of hydrolysis (editing reaction) in leucyl-tRNA synthetase complexed with the misaminoacylated tRNA(Leu), and a novel mechanism of the enzymatic reaction was revealed. Thus, our interface program can play a critical role as a powerful tool for state-of-the-art sophisticated hybrid ab initio QM/MM MD simulations of large systems, such as biological macromolecules.

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.

Spin delocalization over type zero copper

Hard-ligand, high-potential copper sites have been characterized in double mutants of Pseudomonas aeruginosa azurin (C112D/M121X (X = L, F, I)). These sites feature a small A(zz)(Cu) splitting in the EPR spectrum together with enhanced electron transfer activity. Due to these unique properties, these constructs have been called "type zero" copper sites. In contrast, the single mutant, C112D, features a large A(zz)(Cu) value characteristic of the typical type 2 Cu(II). In general, A(zz)(Cu) comprises contributions from Fermi contact, spin dipolar, and orbital dipolar terms. In order to understand the origin of the low A(zz)(Cu) value of type zero Cu(II), we explored in detail its degree of covalency, as manifested by spin delocalization over its ligands, which affects A(zz)(Cu) through the Fermi contact and spin dipolar contributions. This was achieved by the application of several complementary EPR hyperfine spectroscopic techniques at X- and W-band (∼9.5 and 95 GHz, respectively) frequencies to map the ligand hyperfine couplings. Our results show that spin delocalization over the ligands in type zero Cu(II) is different from that of type 2 Cu(II) in the single C112D mutant. The (14)N hyperfine couplings of the coordinated histidine nitrogens are smaller by about 25-40%, whereas that of the (13)C carboxylate of D112 is about 50% larger. From this comparison, we concluded that the spin delocalization of type zero copper over its ligands is not dramatically larger than in type 2 C112D. Therefore, the reduced A(zz)(Cu) value of type zero Cu(II) is largely attributable to an increased orbital dipolar contribution that is related to its larger g(zz) value, as a consequence of the distorted tetrahedral geometry. The increased spin delocalization over the D112 carboxylate in type zero mutants compared to type 2 C112D suggests that electron transfer paths involving this residue are enhanced.

Role of metal in folding and stability of copper proteins in vitro

Metal coordination is required for function of many proteins. For biosynthesis of proteins coordinating a metal, the question arises if the metal binds before, during or after folding of the polypeptide. Moreover, when the metal is bound to the protein, how does its coordination affect biophysical properties such as stability and dynamics? Understanding how metals are utilized by proteins in cells on a molecular level requires accurate descriptions of the thermodynamic and kinetic parameters involved in protein-metal complexes. Copper is one of the essential transition metals found in the active sites of many key proteins. To avoid toxicity of free copper ions, living systems have developed elaborate copper-transport systems that involve dedicated proteins that facilitate efficient and specific delivery of copper to target proteins. This review describes in vitro and in silico biophysical work assessing the role of copper in folding and stability of copper-binding proteins. Examples of proteins discussed are: a blue-copper protein (Pseudomonas aeruginosa azurin), members of copper-transport systems (bacterial CopZ, human Atox1 and ATP7B domains) and multi-copper ferroxidases (yeast Fet3p and human ceruloplasmin). The consequences of interactions between copper proteins and platinum-complexes are also discussed. This article is part of a Special Issue entitled: Cell Biology of Metals.

Predicting folding free energy changes upon single point mutations

MOTIVATION

The folding free energy is an important characteristic of proteins stability and is directly related to protein's wild-type function. The changes of protein's stability due to naturally occurring mutations, missense mutations, are typically causing diseases. Single point mutations made in vitro are frequently used to assess the contribution of given amino acid to the stability of the protein. In both cases, it is desirable to predict the change of the folding free energy upon single point mutations in order to either provide insights of the molecular mechanism of the change or to design new experimental studies.

RESULTS

We report an approach that predicts the free energy change upon single point mutation by utilizing the 3D structure of the wild-type protein. It is based on variation of the molecular mechanics Generalized Born (MMGB) method, scaled with optimized parameters (sMMGB) and utilizing specific model of unfolded state. The corresponding mutations are built in silico and the predictions are tested against large dataset of 1109 mutations with experimentally measured changes of the folding free energy. Benchmarking resulted in root mean square deviation = 1.78 kcal/mol and slope of the linear regression fit between the experimental data and the calculations was 1.04. The sMMGB is compared with other leading methods of predicting folding free energy changes upon single mutations and results discussed with respect to various parameters.

AVAILABILITY

All the pdb files we used in this article can be downloaded from http://compbio.clemson.edu/downloadDir/mentaldisorders/sMMGB_pdb.rar.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Electron transfer reactivity of type zero Pseudomonas aeruginosa azurin

Type zero copper is a hard-ligand analogue of the classical type 1 or blue site in copper proteins that function as electron transfer (ET) agents in photosynthesis and other biological processes. The EPR spectroscopic features of type zero Cu(II) are very similar to those of blue copper, although lacking the deep blue color, due to the absence of thiolate ligation. We have measured the rates of intramolecular ET from the pulse radiolytically generated C3-C26 disulfide radical anion to the Cu(II) in both type zero C112D/M121L and type 2 C112D Pseudomonas aeruginosa azurins in pH 7.0 aqueous solutions between 8 and 45 °C. We also have obtained rate/temperature (10-30 °C) profiles for ET reactions between these mutants and the wild-type azurin. Analysis of the rates and activation parameters for both intramolecular and intermolecular ET reactions indicates that the type zero copper reorganization energy falls in a range (0.9-1.1 eV) slightly above that for type 1 (0.7-0.8 eV), but substantially smaller than that for type 2 (>2 eV), consistent with XAS and EXAFS data that reveal minimal type zero site reorientation during redox cycling.

Mass spectrometric characterization of oligomers in Pseudomonas aeruginosa azurin solutions

We have employed laser-induced liquid bead ion desorption mass spectroscopy (LILBID MS) to study the solution behavior of Pseudomonas aeruginosa azurin as well as two mutants and corresponding Re-labeled derivatives containing a Re(CO)(3)(4,7-dimethyl-1,10-phenanthroline)(+) chromophore appended to a surface histidine. LILBID spectra show broad oligomer distributions whose particular patterns depend on the solution composition (pure H(2)O, 20-30 mM NaCl, 20 and 50 mM NaP(i) or NH(4)P(i) at pH = 7). The distribution maximum shifts to smaller oligomers upon decreasing the azurin concentration and increasing the buffer concentration. Oligomerization is less extensive for native azurin than its mutants. The oligomerization propensities of unlabeled and Re-labeled proteins are generally comparable, and only Re126 shows some preference for the dimer that persists even in highly diluted solutions. Peak shifts to higher masses and broadening in 20-50 mM NaP(i) confirm strong azurin association with buffer ions and solvation. We have found that LILBID MS reveals the solution behavior of weakly bound nonspecific protein oligomers, clearly distinguishing individual components of the oligomer distribution. Independently, average data on oligomerization and the dependence on solution composition were obtained by time-resolved anisotropy of the Re-label photoluminescence that confirmed relatively long rotation correlation times, 6-30 ns, depending on Re-azurin and solution composition. Labeling proteins with Re-chromophores that have long-lived phosphorescence extends the time scale of anisotropy measurements to hundreds of nanoseconds, thereby opening the way for investigations of large oligomers with long rotation times.

Cupredoxins--a study of how proteins may evolve to use metals for bioenergetic processes

Cupredoxins are small proteins that contain type I copper centers, which are ubiquitous in nature. They function as electron transfer shuttles between proteins. This review of the structure and properties of native cupredoxins, and those modified by site-directed mutagenesis, illustrates how these proteins may have evolved to specifically bind copper, develop recognition sites for specific redox partners, tune redox potential for a particular function, and allow for efficient electron transfer through the protein matrix. This is relevant to the general understanding of the roles of metals in energy metabolism, respiration and photosynthesis.

Outer-sphere effects on reduction potentials of copper sites in proteins: the curious case of high potential type 2 C112D/M121E Pseudomonas aeruginosa azurin

Redox and spectroscopic (electronic absorption, multifrequency electron paramagnetic resonance (EPR), and X-ray absorption) properties together with X-ray crystal structures are reported for the type 2 Cu(II) C112D/M121E variant of Pseudomonas aeruginosa azurin. The results suggest that Cu(II) is constrained from interaction with the proximal glutamate; this structural frustration implies a "rack" mechanism for the 290 mV (vs NHE) reduction potential measured at neutral pH. At high pH (∼9), hydrogen bonding in the outer coordination sphere is perturbed to allow axial glutamate ligation to Cu(II), with a decrease in potential to 119 mV. These results highlight the role played by outer-sphere interactions, and the structural constraints they impose, in determining the redox behavior of transition metal protein cofactors.

Transforming a blue copper into a red copper protein: engineering cysteine and homocysteine into the axial position of azurin using site-directed mutagenesis and expressed protein ligation

Interactions of the axial ligand with its blue copper center are known to be important in tuning spectroscopic and redox properties of cupredoxins. While conversion of the blue copper center with a weak axial ligand to a green copper center containing a medium strength axial ligand has been demonstrated in cupredoxins, converting the blue copper center to a red copper center with a strong axial ligand has not been reported. Here we show that replacing Met121 in azurin from Pseudomonas aeruginosa with Cys caused an increased ratio (R(L)) of absorption at 447 nm over that at 621 nm. Whereas no axial Cu-S(Cys121) interaction in Met121Cys was detectable by extended X-ray absorption fine structure (EXAFS) spectroscopy at pH 5, similar to what was observed in native azurin with Met121 as the axial ligand, the Cu-S(Cys121) interaction at 2.74 A is clearly visible at higher pH. Despite the higher R(L) and stronger axial Cys121 interaction with Cu(II) ion, the Met121Cys variant remains largely a type 1 copper protein at low pH (with hyperfine coupling constant A( parallel) = 54 x 10(-4) cm(-1) at pH 4 and 5), or distorted type 1 or green copper protein at high pH (A(parallel) = 87 x 10(-4) cm(-1) at pH 8 and 9), attributable to the relatively long distance between the axial ligand and copper and the constraint placed by the protein scaffold. To shorten the distance between axial ligand and copper, we replaced Met121 with a nonproteinogenic amino acid homocysteine that contains an extra methylene group, resulting in a variant whose spectra (R(L)= 1.5, and A(parallel) = 180 x 10(-4) cm(-1)) and Cu-S(Cys) distance (2.22 A) are very similar to those of the red copper protein nitrosocyanin. Replacing Met121 with Cys or homocysteine resulted in lowering of the reduction potential from 222 mV in the native azurin to 95 +/- 3 mV for Met121Cys azurin and 113 +/- 6 mV for Met121Hcy azurin at pH 7. The results strongly support the "coupled distortion" model that helps explain axial ligand tuning of spectroscopic properties in cupredoxins, and demonstrate the power of using unnatural amino acids to address critical chemical biological questions.

Recent advances in jointed quantum mechanics and molecular mechanics calculations of biological macromolecules: schemes and applications coupled to ab initio calculations

We review the recent research on the functional mechanisms of biological macromolecules using theoretical methodologies coupled to ab initio quantum mechanical (QM) treatments of reaction centers in proteins and nucleic acids. Since in most cases such biological molecules are large, the computational costs of performing ab initio calculations for the entire structures are prohibitive. Instead, simulations that are jointed with molecular mechanics (MM) calculations are crucial to evaluate the long-range electrostatic interactions, which significantly affect the electronic structures of biological macromolecules. Thus, we focus our attention on the methodologies/schemes and applications of jointed QM/MM calculations, and discuss the critical issues to be elucidated in biological macromolecular systems.

Structural characteristics of the hydrophobic patch of azurin and its interaction with p53: a site-directed spin labeling study

Site-directed spin labeling (SDSL) is a powerful tool for monitoring protein structure, dynamics and conformational changes. In this study, the domain-specific properties of azurin and its interaction with p53 were studied using this technique. Mutations of six residues, that are located in the hydrophobic patch of azurin, were prepared and spin labeled. Spectra of the six azurin mutants in solution showed that spin labeled residues 45 and 63 are in a very restricted environment, residues 59 and 65 are in a spacious environment and have free movement, and residues 49 and 51 are located in a relatively closed pocket. Polarity experiments confirmed these results. The changes observed in the spectra of spin labeled azurin upon interaction with p53 indicate that the hydrophobic patch is involved in this interaction. Our results provide valuable insight into the topographic structure of the hydrophobic domain of azurin, as well as direct evidence of its interaction with p53 in solution via the hydrophobic patch. Cytotoxicity studies of azurin mutants showed that residues along the hydrophobic patch are important for its cytotoxicity.

Protein dynamics and pressure: what can high pressure tell us about protein structural flexibility?

After a brief overview of NMR and X-ray crystallography studies on protein flexibility under pressure, we summarize the effects of hydrostatic pressure on the native fold of azurin from Pseudomonas aeruginosa as inferred from the variation of the intrinsic phosphorescence lifetime and the acrylamide bimolecular quenching rate constants of the buried Trp residue. The pressure/temperature response of the globular fold and modulation of its dynamical structure is analyzed both in terms of a reduction of internal cavities and of the hydration of the polypeptide. The study of the effect of two single point cavity forming mutations, F110S and I7S, on the unfolding volume change (ΔV(0)) of azurin and on the internal dynamics of the protein fold under pressure demonstrate that the creation of an internal cavity will enhance the plasticity and lower the stability of the globular structure. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

The 1.4 A resolution structure of Paracoccus pantotrophus pseudoazurin

Pseudoazurins are small type 1 copper proteins that are involved in the flow of electrons between various electron donors and acceptors in the bacterial periplasm, mostly under denitrifying conditions. The previously determined structure of Paracoccus pantotrophus pseudoazurin in the oxidized form was improved to a nominal resolution of 1.4 A, with R and R(free) values of 0.188 and 0.206, respectively. This high-resolution structure makes it possible to analyze the interactions between the monomers and the solvent structure in detail. Analysis of the high-resolution structure revealed the structural regions that are responsible for monomer-monomer recognition during dimer formation and for protein-protein interaction and that are important for partner recognition. The pseudoazurin structure was compared with other structures of various type 1 copper proteins and these were grouped into families according to similarities in their secondary structure; this may be useful in the annotation of copper proteins in newly sequenced genomes and in the identification of novel copper proteins.

A combined atomic force microscopy imaging and docking study to investigate the complex between p53 DNA binding domain and Azurin

The tumor suppressor p53 interacts with the redox copper protein Azurin (AZ) forming a complex which is of some relevance in biomedicine and cancer therapy. To obtain information on the spatial organization of this complex when it is immobilized on a substrate, we have used tapping mode-atomic force microscopy (TM-AFM) imaging combined with computational docking. The vertical dimension and the bearing volume of the DNA binding domain (DBD) of p53, anchored to functionalized gold substrate through exposed lysine residues, alone and after deposing AZ, have been measured by TM-AFM. By a computational docking approach, a three-dimensional model for the DBD of p53, before and after addition of AZ, have been predicted. Then we have calculated the possible arrangements of these biomolecular systems on gold substrate by finding a good agreement with the related experimental distribution of the height. The potentiality of the approach combining TM-AFM imaging and computational docking for the study of biomolecular complexes immobilized on substrates is briefly discussed.

A Proton ENDOR Study of Azurin

As part of our ongoing project that aims at the optimum characterization of the electronic structure of the blue-copper site of azurin from Pseudomonas aeruginosa, we present the complete hyperfine tensors of the protons bound to the Cbeta atom of the copper-bound cysteine 112. These tensors have been obtained from a 95 GHz pulsed electron-nuclear double resonance study of a single crystal of the protein.

Electron Flow through Proteins

Electron transfers in photosynthesis and respiration commonly occur between metal-containing cofactors that are separated by large molecular distances. Employing laser flash-quench triggering methods, we have shown that 20-Å, coupling-limited Fe(II) to Ru(III) and Cu(I) to Ru(III) electron tunneling in Ru-modified cytochromes and blue copper proteins can occur on the microsecond timescale both in solutions and crystals. Redox equivalents can be transferred even longer distances by multistep tunneling, often called hopping, through intervening amino acid side chains. Our work has established that 20-Å hole hopping through an intervening tryptophan is two orders of magnitude faster than single-step electron tunneling in a Re-modified blue copper protein.

Modeling the interaction between the N-terminal domain of the tumor suppressor p53 and azurin

It is known that the half life of the tumor suppressor p53 can be increased by the interaction with the bacterial protein azurin, resulting in an enhanced anti-tumoral activity. The understanding of the molecular mechanisms on the basis of this phenomenon can open the way to new anti-cancer strategies. Some experimental works have given evidence of an interaction between p53 and azurin (AZ); however the binding regions of the proteins are still unknown. Recently, fluorescence studies have shown that p53 partakes in the binding with the bacterial protein by its N-terminal (NT) domain. Here we have used a computational method to get insight into this interacting mode. The model that we propose for the best complex between AZ and p53 has been obtained from a rigid-body docking, coupled with a molecular dynamics (MD) simulation, a free energy calculation, and validated by mutagenesis analysis. We have found a high degree of geometric fit between the two proteins that are kept together by several hydrophobic interactions and numerous hydrogen bonds. Interestingly, it has emerged that AZ binds essentially to the helices H(I) and H(III) of the p53 NT domain, which are also interacting regions for the foremost inhibitor of p53, MDM2.

Metal-binding loop length and not sequence dictates structure

The C-terminal copper-binding loop in the beta-barrel fold of the cupredoxin azurin has been replaced with a range of sequences containing alanine, glycine, and valine residues to assess the importance of amino acid composition and the length of this region. The introduction of 2 and 4 alanines between the coordinating Cys, His, and Met results in loop structures matching those in naturally occurring proteins with the same loop lengths. A loop with 4 alanines between the Cys and His and 3 between the His and Met ligands has a structure identical to that of the WT protein, whose loop is the same length. Loop structure is dictated by length and not sequence allowing the properties of the main surface patch for interactions with partners, to which the loop is a major contributor, to be optimized. Loops with 2 amino acids between the ligands using glycine, alanine, and valine residues have been compared. An empirical relationship is found between copper site protection by the loop and reduction potential. A loop adorned with 4 methyl groups is sufficient to protect the copper ion, enabling most sequences to adequately perform this task. The mutant with 3 alanine residues between the ligands forms a strand-swapped dimer in the crystal structure, an arrangement that has not, to our knowledge, been seen previously for this family of proteins. Cupredoxins function as redox shuttles and are required to be monomeric; therefore, none have evolved with a metal-binding loop of this length.

Active site loop dictates the thermodynamics of reduction and ligand protonation in cupredoxins

The thermodynamics of reduction and His ligand protonation have been determined for a range of loop-contraction variants of the electron transferring type 1 copper protein azurin (AZ). For AZPC, in which the native C-terminal loop containing the Cys, His and Met ligands has been replaced with the shorter sequence from plastocyanin (PC) and AZAMI, in which the even shorter amicyanin (AMI) loop has been inserted, the thermodynamics of reduction match those of the protein whose loop has been introduced which are different to the values for AZ. The enthalpic contribution to His ligand protonation, which is not observed in AZ, is similar in AZAMI and AMI. The thermodynamics of this process in AZPC are more dissimilar to those for PC. In the case of AZAMI-F, a variant possessing the (non natural) minimal loop that can bind a type 1 copper site, the reduction thermodynamics are intermediate between those of AZPC and AZAMI, whilst the thermodynamic data for His ligand protonation are very similar to those for AMI. The results for AZAMI and AZPC are primarily due to protein based enthalpic effects related to the interaction of the metal with permanent protein dipoles from the loop, and to the decreased loop length which favors His ligand protonation in the cuprous proteins. Entropic factors related to loop flexibility have little influence because of constraints imposed by metal coordination and the fact that the introduced loops pack well against the AZ scaffold. Thus, the host scaffold in general plays a minor thermodynamic role in both processes, although for AZAMI-F differences in the first and second coordination spheres influence the thermodynamics of reduction.