X-ray Structure of the Cupredoxin Amicyanin, fromParacoccus denitrificans, Refined at 1.31 Å Resolution

Beyond the coupled distortion model: structural analysis of the single domain cupredoxin AcoP, a green mononuclear copper centre with original features

Cupredoxins are widely occurring copper-binding proteins with a typical Greek-key beta barrel fold. They are generally described as electron carriers that rely on a T1 copper centre coordinated by four ligands provided by the folded polypeptide. The discovery of novel cupredoxins demonstrates the high diversity of this family, with variations in terms of copper-binding ligands, copper centre geometry, redox potential, as well as biological function. AcoP is a periplasmic cupredoxin belonging to the iron respiratory chain of the acidophilic bacterium Acidithiobacillus ferrooxidans. AcoP presents original features, including high resistance to acidic pH and a constrained green-type copper centre of high redox potential. To understand the unique properties of AcoP, we undertook structural and biophysical characterization of wild-type AcoP and of two Cu-ligand mutants (H166A and M171A). The crystallographic structures, including native reduced AcoP at 1.65 Å resolution, unveil a typical cupredoxin fold. The presence of extended loops, never observed in previously characterized cupredoxins, might account for the interaction of AcoP with physiological partners. The Cu-ligand distances, determined by both X-ray diffraction and EXAFS, show that the AcoP metal centre seems to present both T1 and T1.5 features, in turn suggesting that AcoP might not fit well to the coupled distortion model. The crystal structures of two AcoP mutants confirm that the active centre of AcoP is highly constrained. Comparative analysis with other cupredoxins of known structures, suggests that in AcoP the second coordination sphere might be an important determinant of active centre rigidity due to the presence of an extensive hydrogen bond network. Finally, we show that other cupredoxins do not perfectly follow the coupled distortion model as well, raising the suspicion that further alternative models to describe copper centre geometries need to be developed, while the importance of rack-induced contributions should not be underestimated.

On the roles of methionine and the importance of its microenvironments in redox metalloproteins

The amino acid residue methionine (Met) is commonly thought of as a ligand in redox metalloproteins, for example in cytochromes c and in blue copper proteins. However, the roles of Met can go beyond a simple ligand. The thioether functional group of Met allows it to be considered as a hydrophobic residue as well as one that is capable of weak dipolar interactions. In addition, the lone pairs on sulphur allow Met to interact with other groups, inluding the aforementioned metal ions. Because of its properties, Met can play diverse roles in metal coordination, fine tuning of redox reactions, or supporting protein structures. These roles are strongly influenced by the nature of the surrounding medium. Herein, we describe several common interactions between Met and surrounding aromatic amino acids and how they affect the physical properties of both copper and iron metalloproteins. While the importance of interactions between Met and other groups is established in biological systems, less is known about their roles in redox metalloproteins and our view is that this is an area that is ready for greater attention.

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.

Motif-Aware PRALINE: Improving the alignment of motif regions

Protein or DNA motifs are sequence regions which possess biological importance. These regions are often highly conserved among homologous sequences. The generation of multiple sequence alignments (MSAs) with a correct alignment of the conserved sequence motifs is still difficult to achieve, due to the fact that the contribution of these typically short fragments is overshadowed by the rest of the sequence. Here we extended the PRALINE multiple sequence alignment program with a novel motif-aware MSA algorithm in order to address this shortcoming. This method can incorporate explicit information about the presence of externally provided sequence motifs, which is then used in the dynamic programming step by boosting the amino acid substitution matrix towards the motif. The strength of the boost is controlled by a parameter, α. Using a benchmark set of alignments we confirm that a good compromise can be found that improves the matching of motif regions while not significantly reducing the overall alignment quality. By estimating α on an unrelated set of reference alignments we find there is indeed a strong conservation signal for motifs. A number of typical but difficult MSA use cases are explored to exemplify the problems in correctly aligning functional sequence motifs and how the motif-aware alignment method can be employed to alleviate these problems.

The most important functional parts of proteins are often small—but very specific—sequence motifs. Moreover, these motifs tend to be strongly conserved during evolution due to their functional role. Nevertheless, when trying to align protein sequences of the same family, it is often very difficult to align such motifs using standard multiple sequence alignment methods. Aligning functional residues correctly is essential to detect motif conservation, which can be used to filter out spuriously occurring motifs. Additionally, many downstream analyses, such as phylogenetics, are strongly reliant on alignment quality. We have developed a sequence alignment program named Motif-Aware PRALINE (MA-PRALINE) that incorporates information about motifs explicitly. Motifs are provided to MA-PRALINE in the PROSITE pattern syntax; it then scans the input sequences for instances of the pattern and provides a score bonus to matching sequence positions. Our method provides a reproducible alternative to editing alignments by hand in order to account for motif conservation, which is a tedious and error-prone process. We will show that MA-PRALINE allows the alignment of motif-rich regions to be fine-tuned while not degrading the rest of the alignment. MA-PRALINE is available on GitHub as open source software; this allows it to be easily tailored to similar problems. We apply MA-PRALINE on the HIV-1 envelope glycoprotein (gp120) to get an improved alignment of the N-terminal glycosylation motifs. The presence of these motifs is essential for the virus in evading the immune response of the host.

Geometrical Description of Protein Structural Motifs

We present a geometrical method that can identify secondary structural motifs in proteins via angular correlations. The method uses crystal structure coordinates to calculate angular and radial signatures of each residue relative to an external reference point as the number of nearest-neighbor residues increases. We apply our approach to the blue copper protein amicyanin using the copper cofactor as the external reference point. We define a signature termed Δβ which describes the change in angular correlation as the span of nearest neighbor residues increases. We find that three turn regions of amicyanin harbor residues with Δβ near zero, while residues in other secondary structures have Δβ greater than zero: for β-strands, Δβ changes gradually residue by residue along the strand. Extension of our analysis to other blue copper proteins demonstrated that the noted structural trends are general. Importantly, our geometrical description of the folded protein accounts for all forces holding the structure together. Through this analysis, we identified some of the turns in amicyanin as symmetrical anchor points.

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.

Relaxation of structural constraints during Amicyanin unfolding

We study the thermal unfolding of amicyanin by quantifying the resiliency of the native state to structural perturbations. Three signatures characterizing stages of unfolding are identified. The first signature, lateral extension of the polypeptide chain, is calculated directly from the reported crystallographic data. Two other signatures, the radial displacement of each residue from Cu(II) and the angular spread in the chain as the protein unfolds, are calculated using crystallographic data in concert with a geometrical model we introduced previously (J.J. Kozak, H. B. Gray, R. A. Garza-López, J. Inorg. Biochem. 155(2016) 44-55). Particular attention is paid to the resiliency of the two beta sheets in amicyanin. The resiliency of residues in the near neighborhood of the Cu center to destabilization provides information on the persistence of the entatic state. Similarly, examining the resiliency of residues intercalated between structured regions (beta sheets, the alpha helix) provides a basis for identifying a "hydrophobic core." A principal focus of our study is to compare results obtained using our geometrical model with the experimental results (C. La Rosa, D. Milardi, D. M. Grasso, M. P. Verbeet, G. W. Canters, L. Sportelli, R. Guzzi, Eur. Biophy. J.30(8),(2002) 559-570) on the denaturation of amicyanin, and we show that our results support a classical model proposed by these authors.

The sole tryptophan of amicyanin enhances its thermal stability but does not influence the electronic properties of the type 1 copper site

The cupredoxin amicyanin possesses a single tryptophan residue, Trp45. Its fluorescence is quenched when copper is bound even though it is separated by 10.1Å. Mutation of Trp45 to Ala, Phe, Leu and Lys resulted in undetectable protein expression. A W45Y amicyanin variant was isolated. The W45Y mutation did not alter the spectroscopic properties or intrinsic redox potential of amicyanin, but increased the pKa value for the pH-dependent redox potential by 0.5 units. This is due to a hydrogen-bond involving the His95 copper ligand which is present in reduced W45Y amicyanin but not in native amicyanin. The W45Y mutation significantly decreased the thermal stability of amicyanin, as determined by changes in the visible absorbance of oxidized amicyanin and in the circular dichroism spectra for oxidized, reduced and apo forms of amicyanin. Comparison of the crystal structures suggests that the decreased stability of W45Y amicyanin may be attributed to the loss of a strong interior hydrogen bond between Trp45 and Tyr90 in native amicyanin which links two of the β-sheets that comprise the overall structure of amicyanin. Thus, Trp45 is critical for stabilizing the structure of amicyanin but it does not influence the electronic properties of the copper which quenches its fluorescence.

Cross-linking mass spectrometry and mutagenesis confirm the functional importance of surface interactions between CYP3A4 and holo/apo cytochrome b(5)

Cytochrome b(5) (cyt b(5)) is one of the key components in the microsomal cytochrome P450 monooxygenase system. Consensus has not been reached about the underlying mechanism of cyt b(5) modulation of CYP catalysis. Both cyt b(5) and apo b(5) are reported to stimulate the activity of several P450 isoforms. In this study, the surface interactions of both holo and apo b(5) with CYP3A4 were investigated and compared for the first time. Chemical cross-linking coupled with mass spectrometric analysis was used to identify the potential electrostatic interactions between the protein surfaces. Subsequently, the models of interaction of holo/apo b(5) with CYP3A4 were built using the identified interacting sites as constraints. Both cyt b(5) and apo b(5) were predicted to bind to the same groove on CYP3A4 with close contacts to the B-B' loop of CYP3A4, a substrate recognition site. Mutagenesis studies further confirmed that the interacting sites on CYP3A4 (Lys96, Lys127, and Lys421) are functionally important. Mutation of these residues reduced or abolished cyt b(5) binding affinity. The critical role of Arg446 on CYP3A4 in binding to cyt b(5) and/or cytochrome P450 reductase was also discovered. The results indicated that electrostatic interactions on the interface of the two proteins are functionally important. The results indicate that apo b(5) can dock with CYP3A4 in a manner analogous to that of holo b(5), so electron transfer from cyt b(5) is not required for its effects.

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.

Purification and characterization of azurin from the methylamine-utilizing obligate methylotrophMethylobacillus flagellatusKT

Methylamine dehydrogenase (MADH) and azurin were purified from the periplasmic fraction of the methylamine-grown obligate methylotroph Methylobacillus flagellatus KT. The molecular mass of the purified azurin was 16.3 kDa, as measured by SDS–PAGE, or 13 920 Da as determined by MALDI–TOF mass spectrometry. Azurin of M. flagellatus KT contained 1 copper atom per molecule and had an absorption maximum at 620 nm in the oxidized state. The redox potential of azurin measured at pH 7.0 by square-wave voltammetry was +275 mV versus normal hydrogen electrode. MADH reduced azurin in the presence of methylamine, indicating that this cupredoxin is likely to be the physiological electron acceptor for MADH in the electron transport chain of the methylotroph. A scheme of electron transport functioning in M. flagellatus KТ during methylamine oxidation is proposed.

Replacement of the axial copper ligand methionine with lysine in amicyanin converts it to a zinc-binding protein that no longer binds copper

The mutation of the axial ligand of the type I copper protein amicyanin from Met to Lys results in a protein that is spectroscopically invisible and redox inactive. M98K amicyanin acts as a competitive inhibitor in the reaction of native amicyanin with methylamine dehydrogenase indicating that the M98K mutation has not affected the affinity for its natural electron donor. The crystal structure of M98K amicyanin reveals that its overall structure is very similar to native amicyanin but that the type I binding site is occupied by zinc. Anomalous difference Fourier maps calculated using the data collected around the absorption edges of copper and zinc confirm the presence of Zn(2+) at the type I site. The Lys98 NZ donates a hydrogen bond to a well-ordered water molecule at the type I site which enhances the ability of Lys98 to provide a ligand for Zn(2+). Attempts to reconstitute M98K apoamicyanin with copper resulted in precipitation of the protein. The fact that the M98K mutation generated such a selective zinc-binding protein was surprising as ligation of zinc by Lys is rare and this ligand set is unique for zinc.

Proline 96 of the copper ligand loop of amicyanin regulates electron transfer from methylamine dehydrogenase by positioning other residues at the protein-protein interface

Amicyanin is a type 1 copper protein that serves as an electron acceptor for methylamine dehydrogenase (MADH). The site of interaction with MADH is a "hydrophobic patch" of amino acid residues including those that comprise a "ligand loop" that provides three of the four copper ligands. Three prolines are present in this region. Pro94 of the ligand loop was previously shown to strongly influence the redox potential of amicyanin but not affinity for MADH or mechanism of electron transfer (ET). In this study Pro96 of the ligand loop was mutated. P96A and P96G mutations did not affect the spectroscopic or redox properties of amicyanin but increased the K(d) for complex formation with MADH and altered the kinetic mechanism for the interprotein ET reaction. Values of reorganization energy (λ) and electronic coupling (H(AB)) for the ET reaction with MADH were both increased by the mutation, indicating that the true ET reaction observed with native amicyanin was now gated by or coupled to a reconfiguration of the proteins within the complex. The crystal structure of P96G amicyanin was very similar to that of native amicyanin, but notably, in addition to the change in Pro96, the side chains of residues Phe97 and Arg99 were oriented differently. These two residues were previously shown to make contacts with MADH that were important for stabilizing the amicyanin-MADH complex. The values of K(d), λ, and H(AB) for the reactions of the Pro96 mutants with MADH are remarkably similar to those obtained previously for P52G amicyanin. Mutation of this proline, also in the hydrophobic patch, caused reorientation of the side chain of Met51, another reside that interacted with MADH and caused a change in the kinetic mechanism of ET from MADH. These results show that proline residues near the copper site play key roles in positioning other amino acid residues at the amicyanin-MADH interface not only for specific binding to the redox protein partner but also to optimize the orientation of proteins for interprotein ET.

A joint x-ray and neutron study on amicyanin reveals the role of protein dynamics in electron transfer

The joint x-ray/neutron diffraction model of the Type I copper protein, amicyanin from Paracoccus denitrificans was determined at 1.8 A resolution. The protein was crystallized using reagents prepared in D(2)O. About 86% of the amide hydrogen atoms are either partially or fully exchanged, which correlates well with the atomic depth of the amide nitrogen atom and the secondary structure type, but with notable exceptions. Each of the four residues that provide copper ligands is partially deuterated. The model reveals the dynamic nature of the protein, especially around the copper-binding site. A detailed analysis of the presence of deuterated water molecules near the exchange sites indicates that amide hydrogen exchange is primarily due to the flexibility of the protein. Analysis of the electron transfer path through the protein shows that residues in that region are highly dynamic, as judged by hydrogen/deuterium exchange. This could increase the rate of electron transfer by transiently shortening through-space jumps in pathways or by increasing the atomic packing density. Analysis of C-HX bonding reveals previously undefined roles of these relatively weak H bonds, which, when present in sufficient number can collectively influence the structure, redox, and electron transfer properties of amicyanin.

Defining the role of the axial ligand of the type 1 copper site in amicyanin by replacement of methionine with leucine

The effects of replacing the axial methionine ligand of the type 1 copper site with leucine on the structure and function of amicyanin have been characterized. The crystal structures of the oxidized and reduced forms of the protein reveal that the copper site is now tricoordinate with no axial ligand, and that the copper coordination distances for the two ligands provided by histidines are significantly increased. Despite these structural changes, the absorption and EPR spectra of M98L amicyanin are only slightly altered and still consistent with that of a typical type 1 site. The oxidation-reduction midpoint potential (E(m)) value becomes 127 mV more positive as a consequence of the M98L mutation, most likely because of the increased hydrophobicity of the copper site. The most dramatic effect of the mutation was on the electron transfer (ET) reaction from reduced M98L amicyanin to cytochrome c(551i) within the protein ET complex. The rate decreased 435-fold, which was much more than expected from the change in E(m). Examination of the temperature dependence of the ET rate (k(ET)) revealed that the mutation caused a 13.6-fold decrease in the electronic coupling (H(AB)) for the reaction. A similar decrease was predicted from a comparative analysis of the crystal structures of reduced M98L and native amicyanins. The most direct route of ET for this reaction is through the Met98 ligand. Inspection of the structures suggests that the major determinant of the large decrease in the experimentally determined values of H(AB) and k(ET) is the increased distance from the copper to the protein within the type 1 site of M98L amicyanin.

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.

Molecular dynamics of amicyanin reveals a conserved dynamical core for blue copper proteins

Molecular dynamics simulation has been carried out for the blue copper protein amicyanin from two different sources, Paracoccus denitrificans and Paraccocus versutus, to investigate the structural and dynamical properties common to the two molecules and to identify prominent features shared with proteins of the same family, the monomeric cupredoxins. The two amicyanins have almost identical secondary and tertiary structure. In the simulation, they differ for the number of hydrogen bonds in the main chain and the conformation of some beta-strands. However, they strictly maintain the arrangement of the portions of the beta-barrel that are conserved in the folding architecture of the blue copper proteins. Paracoccus versutus amicyanin equilibrates more rapidly, shows lower atomic deviation values, and is less rigid with respect to Paracoccus denitrificans amicyanin. Principal component analysis reveals that the conformational subspaces corresponding to eigenvectors with the same index for each of the two molecules are not necessarily equivalent. Nevertheless, a core scaffold with constrained dynamics exist for both amicyanins. In addition, two fairly flexible regions that are located on the opposite side with respect to the interaction sites with the partner molecules in the redox process have been evidenced in the protein structure. This description of amicyanin, with a few mobile regions remote from the active site and a rigid scaffold including most of the protein beta-barrel, has a close similarity with that of azurin and plastocyanin, two other cupredoxins previously investigated in simulation. Furthermore, similarities in the distribution of the atomic fluctuations indicate that amicyanin, azurin, and plastocyanin possess common dynamical features, in spite of differences in their structure. On the basis of these findings, we suggest that topological constraints imposed by the folding in correspondence of protein regions that are the most conserved determine the protein dynamics of the cupredoxin family. The dynamical properties of the cupredoxins might be controlled for functional advantages that include the binding mechanism with the biological partners and the collective inner motions of the protein matrix required for the electron transfer, whereas long-range conformational changes in the redox reaction should be excluded.

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.

The role of ligand-containing loops at copper sites in proteins

Many approaches are being used to engineer metalloproteins, with most of these informed by, and aiming to further elucidate, the basic structural requirements for biological metal centers. Cupredoxins are type 1 (T1) copper-containing electron transfer (ET) proteins with a -barrel fold that is thought to constrain metal site structure. The T1 copper ion is bound by ligands mainly originating from a single active site loop whose length and structure varies. This Highlight article will focus on protein engineering studies which have investigated the role of the metal-binding loop for active site integrity and functionality. Scaffold differences are present within the cupredoxin family and their influence has also been assessed. Given the widespread occurrence of -barrel domains in nature, and the array of metal sites in proteins composed of loop regions, the studies described on this model system have implications for a variety of metalloproteins.

The axial ligand and extent of protein folding determine whether Zn or Cu binds to amicyanin

M98Q amicyanin is isolated with zinc bound to its type 1 copper-binding site. The influence of the axial ligand of the type 1 copper site on metal specificity is strongest prior to the completion of protein folding and adoption of the final type 1 site geometry. The preference for zinc over copper correlated with the selectivity of apoamicyanin in vitro in the partially folded, rather than the completely folded state. These results suggest that metal incorporation in vivo occurs during protein folding in the periplasm and not to a preformed type 1 site.

Correlation of rhombic distortion of the type 1 copper site of M98Q amicyanin with increased electron transfer reorganization energy

Mutation of the axial Met ligand of the type 1 copper site of amicyanin to Ala or Gln yielded M98A amicyanin, which exhibits typical axial type 1 ligation geometry but with a water molecule providing the axial ligand, and M98Q amicyanin, which exhibits significant rhombic distortion of the type 1 site (Carrell, C. J., Ma, J. K., Antholine, W. E., Hosler, J. P., Mathews, F. S., and Davidson, V. L. (2007) Biochemistry 46, 1900-1912). Despite the change of the axial ligand, the M98Q and M98A mutations had little effect on the redox potential of copper. The true electron transfer (ET) reactions from O-quinol methylamine dehydrogenase to oxidized native and mutant amicyanins revealed that the M98A mutation had little effect on kET, but the M98Q mutation reduced kET 45-fold. Thermodynamic analysis of the latter showed that the decrease in kET was due to an increase of 0.4 eV in the reorganization energy (lambda) associated with the ET reaction to M98Q amicyanin. No change in the experimentally determined electronic coupling or ET distance was observed, confirming that the mutation had not altered the rate-determining step for ET and that this was still a true ET reaction. The basis for the increased lambda is not the nature of the atom that provides the axial ligand because each uses an oxygen from Gln in M98Q amicyanin and from water in M98A amicyanin. Comparisons of the distance of the axial copper ligand from the equatorial plane that is formed by the other three copper ligands in isomorphous crystals of native and mutant amicyanins at atomic resolution indicate an increase in distance from 0.20 A in the native to 0.42 A in M98Q amicyanin and a slight decrease in distance for M98A amicyanin. This correlates with the rhombic distortion caused by the M98Q mutation that is clearly evident in the EPR and visible absorption spectra of the protein and suggests that the extent of rhombicity of the type 1 copper site influences the magnitude of lambda.

Engineering copper sites in proteins: loops confer native structures and properties to chimeric cupredoxins

The ligand-containing loops of two copper-binding electron-transfer proteins (cupredoxins) have been swapped. In the azurin (AZ) variant in which the plastocyanin (PC) sequence is introduced (AZPC), the loop adopts a conformation identical to that in PC. The reduction potential of AZPC is raised as compared to AZ and matches that of PC. In the previously published AZAMI variant (AMI = amicyanin), the shorter introduced loop adopts the same conformation as in AMI, and the reduction potential is lowered to equal that of AMI (Yanagisawa, S.; Dennison, C. J. Am. Chem. Soc. 2004, 126, 15711-15719. Li, C.; et al. Proc. Natl. Acad. Sci. U.S.A. 2006, 103, 7258-7263). Thus, the loop structure plays an important role in tuning the reduction potential of a type 1 copper site with contributions from protein dipoles in this region probably the most important feature. The structure of the loop also seems to be a major factor in controlling dissociation and protonation of the C-terminal His ligand, which can act as a switch to regulate electron-transfer reactivity. The PCAZ variant (PC with the AZ loop) possesses an active site, which is different from those of both PC and AZ, and it is assumed that the introduced loop does not adopt a structure as in AZ. This contributes to the observed instability of PCAZ and highlights that loop-scaffold interactions are important for stabilizing the active site of a cupredoxin.

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

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

Conformer generation under restraints

Conformational sampling by direct optimization of an all-atom energy function is ineffective and inefficient because of the ruggedness of the energy landscape. Discrete sampling schemes represent an attractive alternative for generating ensembles of conformers consistent with spatial restraints derived from empirical data. Conformational sampling is becoming increasingly important for structure prediction as the bottleneck in accurate prediction shifts from energy functions to the methods used to find low-energy conformers. Experimental structure determination remains a perennial challenge as investigators tackle larger macromolecular systems, and begin to incorporate more complete descriptions of uncertainty, heterogeneity and dynamics into their models. Computational approaches that combine dense, discrete sampling with all-atom energy evaluation and refinement may help to overcome the remaining barriers to solving these problems.

Heterogeneity and Inaccuracy in Protein Structures Solved by X-Ray Crystallography

Proteins are dynamic molecules, exhibiting structural heterogeneity in the form of anisotropic motion and discrete conformational substates, often of functional importance. In protein structure determination by X-ray crystallography, the observed diffraction pattern results from the scattering of X-rays by an ensemble of heterogeneous molecules, ordered and oriented by packing in a crystal lattice. The majority of proteins diffract to resolutions where heterogeneity is difficult to identify and model, and are therefore approximated by a single, average conformation with isotropic variance. Here we show that disregarding structural heterogeneity introduces degeneracy into the structure determination process, as many single, isotropic models exist that explain the diffraction data equally well. The large differences among these models imply that the accuracy of crystallographic structures has been widely overestimated. Further, it suggests that analyses that depend on small differences in the relative positions of atoms may be flawed.

An engineered CuA Amicyanin capable of intermolecular electron transfer reactions

The type I copper center of amicyanin was replaced with a binuclear CuA center. To create this model CuA protein, a portion of the amino acid sequence that contains three of the ligands to the native type I copper center of Paracoccus denitrificans amicyanin was replaced with the corresponding portion of sequence that provides five ligands for the CuA center of cytochrome c oxidase from P. denitrificans. UV-visible and electron paramagnetic resonance spectroscopy confirm that the engineered protein as isolated possesses the mixed-valence Cu1.5Cu1.5 (purple) CuA center. Comparison of the spectroscopic properties of this CuA amicyanin with those of the CuA centers of other natural and engineered CuA proteins suggests that the spectroscopic features may be dictated more by the protein host than the sequence of the CuA loop. Novel reactions for a simple CuA model protein are also described. In contrast to other natural and engineered CuA proteins, the fully reduced CuA amicyanin may be reoxidized by molecular oxygen to the mixed-valence state. It is also shown that CuA amicyanin can serve as an electron donor and an electron acceptor for other redox proteins. The mixed-valence form accepts electrons from cytochromes c-551i and c-550 from P. denitrificans. The fully reduced form donates electrons to native and P94F amicyanin. The function as either an electron donor or acceptor is consistent with the measured redox potential of CuA amicyanin of +273 mV. These data indicate that this CuA amicyanin will be a particularly useful model protein for structure-function studies of reactivity and the electron transfer properties of the CuA redox center.

Metal-ligand interplay in blue copper proteins studied by 1H NMR spectroscopy: Cu(II)-pseudoazurin and Cu(II)-rusticyanin

The blue copper proteins (BCPs), pseudoazurin from Achromobacter cycloclastes and rusticyanin from Thiobacillus ferrooxidans, have been investigated by (1)H NMR at a magnetic field of 18.8 T. Hyperfine shifts of the protons belonging to the coordinated ligands have been identified by exchange spectroscopy, including the indirect detection for those resonances that cannot be directly observed (the beta-CH(2) of the Cys ligand, and the NH amide hydrogen bonded to the S(gamma)(Cys) atom). These data reveal that the Cu(II)-Cys interaction in pseudoazurin and rusticyanin is weakened compared to that in classic blue sites (plastocyanin and azurin). This weakening is not induced by a stronger interaction with the axial ligand, as found in stellacyanin, but might be determined by the protein folding around the metal site. The average chemical shift of the beta-CH(2) Cys ligand in all BCPs can be correlated to geometric factors of the metal site (the Cu-S(gamma)(Cys) distance and the angle between the CuN(His)N(His) plane and the Cu-S(gamma)(Cys) vector). It is concluded that the degree of tetragonal distortion is not necessarily related to the strength of the Cu(II)-S(gamma)(Cys) bond. The copper-His interaction is similar in all BCPs, even for the solvent-exposed His ligand. It is proposed that the copper xy magnetic axes in blue sites are determined by subtle geometrical differences, particularly the orientation of the His ligands. Finally, the observed chemical shifts for beta-CH(2) Cys and Ser NH protons in rusticyanin suggest that a less negative charge at the sulfur atom could contribute to the high redox potential (680 mV) of this protein.

Crystal structures of the Met148Leu and Ser86Asp mutants of rusticyanin from Thiobacillus ferrooxidans: insights into the structural relationship with the cupredoxins and the multi copper proteins

The crystal structures of the Met148Leu and Ser86Asp mutants of rusticyanin are presented at 1.82 and 1.65 A resolution, respectively. Both of these structures have two molecules in the asymmetric unit compared to the one present in the crystal form of the native protein. This provides an opportunity to investigate intramolecular electron transfer pathways in rusticyanin. The redox potential of the Met148Leu mutant ( approximately 800 mV) is elevated compared to that of the native protein ( approximately 670 mV at pH 3.2) while that of the Ser86Asp mutant ( approximately 623 mV at pH 3.2) is decreased. The effect of the Ser86Asp mutation on the hydrogen bonding near the type 1 Cu site is discussed and hence its role in determining acid stability is examined. The type 1 Cu site of Met148Leu mimics the structural and biochemical characteristics of those found in domain II of ceruloplasmin and fungal laccase. Moreover, the native rusticyanin's cupredoxin core and the type 1 Cu site closely resemble those found in ascorbate oxidase and nitrite reductase. Structure based phylogenetic trees have been re-examined in view of the additional structural data on rusticyanin and fungal laccase. We confirm that rusticyanin is in the same class as nitrite reductase domain 2, laccase domain 3 and ceruloplasmin domains 2, 4 and 6.

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.

Electronic characterization of the oxidized state of the blue copper protein rusticyanin by 1H NMR: is the axial methionine the dominant influence for the high redox potential?

The oxidized state of rusticyanin, the blue copper protein with the highest redox potential in its class, has been investigated through (1)H nuclear magnetic resonance applied to its cobalt(II) derivative. The assignment of the protons belonging to the coordinated residues has been performed. Many other amino acids situated in the vicinity of the metal ion, including six hydrophobic residues (isoleucine140 and five phenylalanines) have also been identified. The orientation of the main axes of the magnetic susceptibility tensor for the cobalt(II)-rusticyanin as well as its axial, Deltachi(ax), and rhombic, Deltachi(rh), magnetic susceptibility anisotropy components have been determined. A comparison of the present results with those previously obtained for cobalt(II)azurin [Donaire, A., Salgado, J., Moratal, J. M. (1998) Biochemistry 37, 8659-8673] allows us to provide further insights into the reasons for the high redox potential of this protein. According to our results, the interaction between the metal ion and the thioether Sdelta of the axial methionine is not as influential as the strong destabilizing effect that the hydrophobic residues close to the metal ion undergo in the oxidized state.