Crystal structure analysis of amicyanin and apoamicyanin from Paracoccus denitrificans at 2.0 A and 1.8 A resolution

Substitution of the sole tryptophan of the cupredoxin, amicyanin, with 5-hydroxytryptophan alters fluorescence properties and energy transfer to the type 1 copper site

Amicyanin is a type 1 copper protein with a single tryptophan residue. Using genetic code expansion, the tryptophan was selectively replaced with the unnatural amino acid, 5-hydroxytryptophan (5-HTP). The 5-HTP substituted amicyanin exhibited absorbance at 300-320 nm, characteristic of 5-HTP and not seen in native amicyanin. The fluorescence emission maximum in 5-HTP substituted amicyanin is redshifted from 318 nm in native amicyanin to 331 nm and to 348 nm in the unfolded protein. The fluorescence quantum yield of 5-HTP substituted amicyanin mutant was much less than that of native amicyanin. Differences in intrinsic fluorescence are explained by differences in the excited states of tryptophan versus 5-HTP and the intraprotein environment. The substitution of tryptophan with 5-HTP did not affect the visible absorbance and redox potential of the copper, which is 10 Å away. In amicyanin and other cupredoxins, an unexplained quenching of the intrinsic fluorescence by the bound copper is observed. However, the fluorescence of 5-HTP substituted amicyanin is not quenched by the copper. It is shown that the mechanism of quenching in native amicyanin is Förster, or fluorescence, resonance energy transfer (FRET). This does not occur in 5-HTP substituted amicyanin because the fluorescence quantum yield is significantly lower and the red-shift of fluorescence emission maximum decreases overlap with the near UV absorbance of copper. Characterization of the distinct fluorescence properties of 5-HTP relative to tryptophan in amicyanin provides a basis for spectroscopic interrogation of the protein microenvironment using 5-HTP, and long-distance interactions with transition metals.

Metalloprotein catalysis: structural and mechanistic insights into oxidoreductases from neutron protein crystallography

Metalloproteins catalyze a range of reactions, with enhanced chemical functionality due to their metal cofactor. The reaction mechanisms of metalloproteins have been experimentally characterized by spectroscopy, macromolecular crystallography and cryo-electron microscopy. An important caveat in structural studies of metalloproteins remains the artefacts that can be introduced by radiation damage. Photoreduction, radiolysis and ionization deriving from the electromagnetic beam used to probe the structure complicate structural and mechanistic interpretation. Neutron protein diffraction remains the only structural probe that leaves protein samples devoid of radiation damage, even when data are collected at room temperature. Additionally, neutron protein crystallography provides information on the positions of light atoms such as hydrogen and deuterium, allowing the characterization of protonation states and hydrogen-bonding networks. Neutron protein crystallography has further been used in conjunction with experimental and computational techniques to gain insight into the structures and reaction mechanisms of several transition-state metal oxidoreductases with iron, copper and manganese cofactors. Here, the contribution of neutron protein crystallography towards elucidating the reaction mechanism of metalloproteins is reviewed.

Copper-binding energetics of amicyanin in different folding states

Amicyanin is a type I copper protein that mediates electron transfer between methylamine dehydrogenase and cytochrome c-_551i for energy production in Paracoccus denitrificans. Although the Met98 axial ligand of amicyanin has been shown to dictate metal selectivity and specificity during protein folding, the mechanism involved in copper-mediated amicyanin folding is unknown. Here, we kinetically and spectroscopically described reaction steps for incorporating copper into fully and less folded apo-amicyanin and established thermodynamic parameters for two amicyanin folding states. The rate constant for the incorporation of copper into fully folded apo-amicyanin at 25 °C was almost 1.5-fold lower than that for the initial phase of copper addition to the less folded apo-amicyanin. However, the rate constant was 10-fold higher than that of the second phase of copper addition to less folded apo-amicyanin at 25 °C. When overall binding energetic parameters (ΔH° and ΔS°) for the incorporation of copper into fully folded apo-amicyanin were measured by the van't Hoff method and isothermal titration calorimetry, the values were more positive than those determined for less folded apo-amicyanin. This indicates that during amicyanin biogenesis, copper rapidly binds to an unfolded apo-amicyanin active site, inducing protein folding and favorably influencing subsequent organization of copper ligands.

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.

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.

Converting the bis-FeIV state of the diheme enzyme MauG to Compound I decreases the reorganization energy for electron transfer

The electron transfer (ET) properties of two types of high-valent hemes were studied within the same protein matrix; the bis-Fe(IV) state of MauG and the Compound I state of Y294H MauG. The latter is formed as a consequence of mutation of the tyrosine which forms the distal axial ligand of the six-coordinate heme that allows it to stabilize Fe(IV) in the absence of an external ligand. The rates of the ET reaction of each high-valent species with the type I copper protein, amicyanin, were determined at different temperatures and analysed by ET theory. The reaction with bis-Fe(IV) wild-type (WT) MauG exhibited a reorganization energy (λ) that was 0.39 eV greater than that for the reaction of Compound I Y295H MauG. It is concluded that the delocalization of charge over the two hemes in the bis-Fe(IV) state is responsible for the larger λ, relative to the Compound I state in which the Fe(V) equivalent is isolated on one heme. Although the increase in λ decreases the rate of ET, the delocalization of charge decreases the ET distance to its natural substrate protein, thus increasing the ET rate. This describes how proteins can balance different ET properties of complex redox cofactors to optimize each system for its particular ET or catalytic reaction.

Characterization of the free energy dependence of an interprotein electron transfer reaction by variation of pH and site-directed mutagenesis

The interprotein electron transfer (ET) reactions of the cupredoxin amicyanin, which mediates ET from the tryptophan tryptophylquinone (TTQ) cofactor of methylamine dehydrogenase to cytochrome c-551i have been extensively studied. However, it was not possible to perform certain key experiments in that native system. This study examines the ET reaction from reduced amicyanin to an alternative electron acceptor, the diheme protein MauG. It was possible to vary the ΔG° for this ET reaction by simply changing pH to determine the dependence of kET on ΔG°. A P94A mutation of amicyanin significantly altered its oxidation-reduction midpoint potential value. It was not possible to study the ET from reduced P94A amicyanin to cytochrome c-551i in the native system because that reaction was kinetically coupled. However, the reaction from reduced P94A amicyanin to MauG was a true ET reaction and it was possible to determine values of reorganization energy (λ) and electronic coupling for the reactions of this variant as well as native amicyanin. Comparison of the λ values associated with the ET reactions between amicyanin and the TTQ of methylamine dehydrogenase, the diheme center of MauG and the single heme of cytochrome c-551i, provides insight into the factors that dictate the λ values for the respective reactions. These results demonstrate how study of ET reactions with alternative redox partner proteins can complement and enhance our understanding of the reactions with the natural redox partners, and further our understanding of mechanisms of protein ET reactions.

Evaluating membrane affinity by integrating protein orientations

Energetic interactions of a protein with lipid bilayers determine its propensity to reside in the membrane. Here we seek to evaluate the membrane interactions for EMAPII, a protein found to be released from the cell by unknown mechanisms, as well as several other proteins. Using a knowledge-based coarse-grained membrane potential, we calculate the free energy profiles for these proteins by integrating out the orientation degrees of freedom. Due to the invariance of energy under in-plane rotations about the membrane normal, the orientation space can be reduced to two dimensions and mapped onto the surface of a unit sphere, thus making visualization, sampling and integration more convenient. The integrated free energy profiles determine the relative probabilities along the membrane normal for the proteins regardless of their orientations, and display distinctive characteristics for membrane proteins and water-soluble proteins. The membrane interactions for EMAPII exhibit typical features of a water-soluble protein, with a high energetic barrier to enter or cross the membrane. Our results thus suggest that similar to the non-classical export of FGF1, the release of EMAPII would involve more complicated mechanisms than simple passive diffusion across the membrane.

Mechanisms for control of biological electron transfer reactions

Electron transfer (ET) through and between proteins is a fundamental biological process. The rates and mechanisms of these ET reactions are controlled by the proteins in which the redox centers that donate and accept electrons reside. The protein influences the magnitudes of the ET parameters, the electronic coupling and reorganization energy that are associated with the ET reaction. The protein can regulate the rates of the ET reaction by requiring reaction steps to optimize the system for ET, leading to kinetic mechanisms of gated or coupled ET. Amino acid residues in the segment of the protein through which long range ET occurs can also modulate the ET rate by serving as staging points for hopping mechanisms of ET. Specific examples are presented to illustrate these mechanisms by which proteins control rates of ET reactions.

Nature's palette: the search for natural blue colorants

The food and beverage industry is seeking to broaden the palette of naturally derived colorants. Although considerable effort has been devoted to the search for new blue colorants in fruits and vegetables, less attention has been directed toward blue compounds from other sources such as bacteria and fungi. The current work reviews known organic blue compounds from natural plant, animal, fungal, and microbial sources. The scarcity of blue-colored metabolites in the natural world relative to metabolites of other colors is discussed, and structural trends common among natural blue compounds are identified. These compounds are grouped into seven structural classes and evaluated for their potential as new color additives.

A simple method to engineer a protein-derived redox cofactor for catalysis

The 6×-Histidine tag which is commonly used for purification of recombinant proteins was converted to a catalytic redox-active center by incorporation of Co(2+). Two examples of the biological activity of this engineered protein-derived cofactor are presented. After inactivation of the natural diheme cofactor of MauG, it was shown that the Co(2+)-loaded 6×His-tag could substitute for the hemes in the H2O2-driven catalysis of tryptophan tryptophylquinone biosynthesis. To further demonstrate that the Co(2+)-loaded 6×His-tag could mediate long range electron transfer, it was shown that addition of H2O2 to the Co(2+)-loaded 6×His-tagged Cu(1+) amicyanin oxidizes the copper site which is 20Å away. These results provide proof of principle for this simple method by which to introduce a catalytic redox-active site into proteins for potential applications in research and biotechnology.

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.

DJ-1 is a copper chaperone acting on SOD1 activation

Lack of oxidative stress control is a common and often prime feature observed in many neurodegenerative diseases. Both DJ-1 and SOD1, proteins involved in familial Parkinson disease and amyotrophic lateral sclerosis, respectively, play a protective role against oxidative stress. Impaired activity and modified expression of both proteins have been observed in different neurodegenerative diseases. A potential cooperative action of DJ-1 and SOD1 in the same oxidative stress response pathway may be suggested based on a copper-mediated interaction between the two proteins reported here. To investigate the mechanisms underlying the antioxidative function of DJ-1 in relation to SOD1 activity, we investigated the ability of DJ-1 to bind copper ions. We structurally characterized a novel copper binding site involving Cys-106, and we investigated, using different techniques, the kinetics of DJ-1 binding to copper ions. The copper transfer between the two proteins was also examined using both fluorescence spectroscopy and specific biochemical assays for SOD1 activity. The structural and functional analysis of the novel DJ-1 copper binding site led us to identify a putative role for DJ-1 as a copper chaperone. Alteration of the coordination geometry of the copper ion in DJ-1 may be correlated to the physiological role of the protein, to a potential failure in metal transfer to SOD1, and to successive implications in neurodegenerative etiopathogenesis.

Molecular dynamics simulations of apocupredoxins: insights into the formation and stabilization of copper sites under entatic control

Cupredoxins perform copper-mediated long-range electron transfer (ET) in biological systems. Their copper-binding sites have evolved to force copper ions into ET-competent systems with decreased reorganization energy, increased reduction potential, and a distinct electronic structure compared with those of non-ET-competent copper complexes. The entatic or rack-induced state hypothesis explains these special properties in terms of the strain that the protein matrix exerts on the metal ions. This idea is supported by X-ray structures of apocupredoxins displaying "closed" arrangements of the copper ligands like those observed in the holoproteins; however, it implies completely buried copper-binding atoms, conflicting with the notion that they must be exposed for copper loading. On the other hand, a recent work based on NMR showed that the copper-binding regions of apocupredoxins are flexible in solution. We have explored five cupredoxins in their "closed" apo forms through molecular dynamics simulations. We observed that prearranged ligand conformations are not stable as the X-ray data suggest, although they do form part of the dynamic landscape of the apoproteins. This translates into variable flexibility of the copper-binding regions within a rigid fold, accompanied by fluctuations of the hydrogen bonds around the copper ligands. Major conformations with solvent-exposed copper-binding atoms could allow initial binding of the copper ions. An eventual subsequent incursion to the closed state would result in binding of the remaining ligands, trapping the closed conformation thanks to the additional binding energy and the fastening of noncovalent interactions that make up the rack.

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.

The X-ray absorption spectroscopic model of the copper(II) imidazole complex ion in liquid aqueous solution: a strongly solvated square pyramid

Cu K-edge extended X-ray absorption fine structure (EXAFS) and Minuit X-ray absorption near-edge structure (MXAN) analyses were combined to evaluate the structure of the copper(II) imidazole complex ion in liquid aqueous solution. Both methods converged to the same square-pyramidal inner coordination sphere [Cu(Im)(4)L(ax)](2+) (L(ax) indeterminate) with four equatorial nitrogen atoms at EXAFS, 2.02 ± 0.01 Å, and MXAN, 1.99 ± 0.03 Å. A short-axial N/O scatterer (L(ax)) was found at 2.12 ± 0.02 Å (EXAFS) or 2.14 ± 0.06 Å (MXAN). A second but very weak axial Cu-N/O interaction was found at 2.9 ± 0.1 Å (EXAFS) or 3.0 ± 0.1 Å (MXAN). In the MXAN fits, only a square-pyramidal structural model successfully reproduced the doubled maximum of the rising K-edge X-ray absorption spectrum, specifically excluding an octahedral model. Both EXAFS and MXAN also found eight outlying oxygen scatterers at 4.2 ± 0.3 Å that contributed significant intensity over the entire spectral energy range. Two prominent rising K-edge shoulders at 8987.1 and 8990.5 eV were found to reflect multiple scattering from the 3.0 Å axial scatterer and the imidazole rings, respectively. In the MXAN fits, the imidazole rings took in-plane rotationally staggered positions about copper. The combined (EXAFS and MXAN) model for the unconstrained cupric imidazole complex ion in liquid aqueous solution is an axially elongated square-pyramidal core, with a weak nonbonded interaction at the second axial coordination position and a solvation shell of eight nearest-neighbor water molecules. This core square-pyramidal motif has persisted through [Cu(H(2)O)(5)](2+), [Cu(NH(3))(4)(NH(3),H(2)O)](2+), (1, 2) and now [Cu(Im)(4)L(ax))](2+) and appears to be the geometry preferred by unconstrained aqueous-phase copper(II) complex ions.

Tryptophan tryptophylquinone biosynthesis: a radical approach to posttranslational modification

Protein-derived cofactors are formed by irreversible covalent posttranslational modification of amino acid residues. An example is tryptophan tryptophylquinone (TTQ) found in the enzyme methylamine dehydrogenase (MADH). TTQ biosynthesis requires the cross-linking of the indole rings of two Trp residues and the insertion of two oxygen atoms onto adjacent carbons of one of the indole rings. The diheme enzyme MauG catalyzes the completion of TTQ within a precursor protein of MADH. The preMADH substrate contains a single hydroxyl group on one of the tryptophans and no crosslink. MauG catalyzes a six-electron oxidation that completes TTQ assembly and generates fully active MADH. These oxidation reactions proceed via a high valent bis-Fe(IV) state in which one heme is present as Fe(IV)=O and the other is Fe(IV) with both axial heme ligands provided by amino acid side chains. The crystal structure of MauG in complex with preMADH revealed that catalysis does not involve direct contact between the hemes of MauG and the protein substrate. Rather it is accomplished through long-range electron transfer, which presumably generates radical intermediates. Kinetic, spectrophotometric, and site-directed mutagenesis studies are beginning to elucidate how the MauG protein controls the reactivity of the hemes and mediates the long range electron/radical transfer required for catalysis. This article is part of a Special Issue entitled: Radical SAM enzymes and Radical Enzymology.

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.

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.

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 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.

EfeO-cupredoxins: major new members of the cupredoxin superfamily with roles in bacterial iron transport

The EfeUOB system of Escherichia coli is a tripartite, low pH, ferrous iron transporter. It resembles the high-affinity iron transporter (Ftr1p-Fet3p) of yeast in that EfeU is homologous to Ftr1p, an integral-membrane iron-permease. However, EfeUOB lacks an equivalent of the Fet3p component--the multicopper oxidase with three cupredoxin-like domains. EfeO and EfeB are periplasmic but their precise roles are unclear. EfeO consists primarily of a C-terminal peptidase-M75 domain with a conserved 'HxxE' motif potentially involved in metal binding. The smaller N-terminal domain (EfeO-N) is predicted to be cupredoxin (Cup) like, suggesting a previously unrecognised similarity between EfeO and Fet3p. Our structural modelling of the E. coli EfeO Cup domain identifies two potential metal-binding sites. Site I is predicted to bind Cu(2+) using three conserved residues (C41 and 103, and E66) and M101. Of these, only one (C103) is conserved in classical cupredoxins where it also acts as a Cu ligand. Site II most probably binds Fe(3+) and consists of four well conserved surface Glu residues. Phylogenetic analysis indicates that the EfeO-Cup domains form a novel Cup family, designated the 'EfeO-Cup' family. Structural modelling of two other representative EfeO-Cup domains indicates that different subfamilies employ distinct ligand sets at their proposed metal-binding sites. The ~100 efeO homologues in the bacterial sequence databases are all associated with various iron-transport related genes indicating a common role for EfeO-Cup proteins in iron transport, supporting a new copper-iron connection in biology.

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.

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.

Respiratory transformation of nitrous oxide (N2O) to dinitrogen by Bacteria and Archaea

N2O is a potent greenhouse gas and stratospheric reactant that has been steadily on the rise since the beginning of industrialization. It is an obligatory inorganic metabolite of denitrifying bacteria, and some production of N2O is also found in nitrifying and methanotrophic bacteria. We focus this review on the respiratory aspect of N2O transformation catalysed by the multicopper enzyme nitrous oxide reductase (N2OR) that provides the bacterial cell with an electron sink for anaerobic growth. Two types of Cu centres discovered in N2OR were both novel structures among the Cu proteins: the mixed-valent dinuclear Cu(A) species at the electron entry site of the enzyme, and the tetranuclear Cu(Z) centre as the first catalytically active Cu-sulfur complex known. Several accessory proteins function as Cu chaperone and ABC transporter systems for the biogenesis of the catalytic centre. We describe here the paradigm of Z-type N2OR, whose characteristics have been studied in most detail in the genera Pseudomonas and Paracoccus. Sequenced bacterial genomes now provide an invaluable additional source of information. New strains harbouring nos genes and capability of N2O utilization are being uncovered. This reveals previously unknown relationships and allows pattern recognition and predictions. The core nos genes, nosZDFYL, share a common phylogeny. Most principal taxonomic lineages follow the same biochemical and genetic pattern and share the Z-type enzyme. A modified N2OR is found in Wolinella succinogenes, and circumstantial evidence also indicates for certain Archaea another type of N2OR. The current picture supports the view of evolution of N2O respiration prior to the separation of the domains Bacteria and Archaea. Lateral nos gene transfer from an epsilon-proteobacterium as donor is suggested for Magnetospirillum magnetotacticum and Dechloromonas aromatica. In a few cases, nos gene clusters are plasmid borne. Inorganic N2O metabolism is associated with a diversity of physiological traits and biochemically challenging metabolic modes or habitats, including halorespiration, diazotrophy, symbiosis, pathogenicity, psychrophily, thermophily, extreme halophily and the marine habitat down to the greatest depth. Components for N2O respiration cover topologically the periplasm and the inner and outer membranes. The Sec and Tat translocons share the task of exporting Nos components to their functional sites. Electron donation to N2OR follows pathways with modifications depending on the host organism. A short chronology of the field is also presented.

The ligand geometry of copper determines the stability of amicyanin

Solution differential scanning calorimetry (DSC) of oxidized amicyanin, a Type I copper protein, at pH 7.5 reveals two thermal transitions. The major transition at 67.7 degrees C corresponds to the disruption of the Cys(92) thiolate to Cu(II) charge transfer as evidenced by a corresponding temperature-dependent loss of amicyanin visible absorbance. A minor transition at 75.5 degrees C describes the further irreversible protein unfolding. Reduced amicyanin exhibits a pH-dependent change of the copper ligand geometry. At pH 8.5 where the Type I tetrahedral geometry is maintained, DSC reveals two thermal transitions with T(m) values similar to that of oxidized amicyanin. At pH 6.2 where the Cu(I) coordination is trigonal planar, reduced amicyanin exhibits a single thermal transition with a lower T(m) of 64.0 degrees C. Apoamicyanin, from which copper has been removed, also exhibits a single thermal transition but with a much lower T(m) of 51.8 degrees C. Thus, the thermal stability of amicyanin is dictated both by the presence or absence of copper and its ligand geometry, but not its redox state. The physiological relevance of these data is discussed.

The effect of resolution-dependent global shape modifications on rigid-body protein-protein docking

Docking unbound molecules presents a challenge in the case where no prior biological or bioinformatic knowledge exists. This is mainly due to differences between the structures of the molecules when in a complex and in the free state. Presumably, these differences interfere with the ability of protein-protein docking algorithms, which rely on a dominant shape descriptor, to identify the correct solution and rank it higher than false solutions. In this study we verify the notion that small discords in the molecular fit can be eliminated by using appropriately designed low-resolution shape descriptors, thereby improving the docking results. We exploit the inherent gradual resolution dependency of Fourier transforms and formulate a resolution-dependent shape descriptor by truncating selected Fourier transform terms. Thus, different levels of shape modification are attained, affecting the degree of detail in the depiction of the molecular surface. We applied the modified descriptor to a selection of 23 protein-protein systems, using the unbound structures where possible. The docking results obtained with the new geometric descriptor were considerably superior to former results, improving the ranks of nearly correct solutions for 17 systems. Unification of the results of scans in which different resolutions were employed further improved the ranks of nearly correct solutions to less than 100 for 12 of the 23 systems and less than 300 for 20 systems. The new geometric descriptor can be combined with other descriptors, which typify the electrostatic or hydrophobic character of the molecular surface, and with external experimental or bioinformatic data.

The linked conservation of structure and function in a family of high diversity: the monomeric cupredoxins

The monomeric cupredoxins are a highly divergent family of copper binding electron transport proteins that function in photosynthesis and respiration. To determine how function and structure are conserved in the context of large sequence differences, we have carried out a detailed analysis of the cupredoxins of known structure and their sequence homologs. The common structure of the cupredoxins is formed by a sandwich of two beta sheets which support a copper binding site. The structure of the deeply buried core is intimately coupled to the binding site on the surface of the protein; in each protein the conserved regions form one continuous substructure that extends from the surface active site and through the center of the molecule. Residues around the active site are conserved for functional reasons, while those deeper in the structure will be conserved for structural reasons. Together the two sets support each other.

Hydrophobic complementarity in protein-protein docking

Formation of hydrophobic contacts across a newly formed interface is energetically favorable. Based on this observation we developed a geometric-hydrophobic docking algorithm that estimates quantitatively the hydrophobic complementarity at protein-protein interfaces. Each molecule to be docked is represented as a grid of complex numbers, storing information regarding the shape of the molecule in the real part and information regarding the hydropathy of the surface in the imaginary part. The grid representations are correlated using fast Fourier transformations. The algorithm is used to compare the extent of hydrophobic complementarity in oligomers (represented by D2 tetramers) and in hetero-dimers of soluble proteins (complexes). We also test the implication of hydrophobic complementarity in distinguishing correct from false docking solutions. We find that hydrophobic complementarity at the interface exists in oligomers and in complexes, and in both groups the extent of such complementarity depends on the size of the interface. Thus, the non-polar portions of large interfaces are more often juxtaposed than non-polar portions of small interfaces. Next we find that hydrophobic complementarity helps to point out correct docking solutions. In oligomers it significantly improves the ranks of nearly correct reassembled and modeled tetramers. Combining geometric, electrostatic and hydrophobic complementarity for complexes gives excellent results, ranking a nearly correct solution < 10 for 5 of 23 tested systems, < 100 for 8 systems and < 1000 for 19 systems.

Amicyanin metal-site structure and interaction with MADH: PAC and NMR spectroscopy of Ag-, Cd-, and Cu-amicyanin

To investigate the structural control mechanisms in the metal site of amicyanin when interacting with MADH, redox-inactive Ag(+)- and Cd(2+)-substituted amicyanins were studied with perturbed angular correlations of gamma-rays (PAC) spectroscopy. PAC experiments on (111m)Cd-substituted amicyanin revealed two different metal-site structures, which are very likely in dynamic exchange on a ~5 ns timescale. Only one structure binds to MADH. The dissociation constants, K(d), are 9+/-2 microM with MADH(red) and 38+/-11 microM with MADH(ox), indicating that the Cd-amicyanin binding affinity is regulated by the MADH redox state. PAC experiments on (111)Ag-substituted amicyanin also showed two different forms of Ag-amicyanin, probably reflecting relaxation from Ag to Cd geometry. No binding of Ag-amicyanin to MADH could be observed with PAC, suggesting that the K(d) is larger than 43 microM, based on the 95% confidence limit. NMR revealed large chemical shift differences between native copper amicyanin and both metal-substituted forms. Affected residues are found up to 15 A away from the metal ion. The Ag(+)- and Cd(2+)-substituted amicyanins demonstrate no change in coordination as a function of pH, contrary to Cu(+)-amicyanin which shows protonation of the copper ligand His96 with p K(a)=6.8. It is concluded that, contrary to other blue copper proteins, Ag(+)-amicyanin is not a close mimic of Cu(+)-amicyanin, and that structural changes in the metal site have large effects on the affinity for the redox partner.

Blue copper proteins: a comparative analysis of their molecular interaction properties

Blue copper proteins are type-I copper-containing redox proteins whose role is to shuttle electrons from an electron donor to an electron acceptor in bacteria and plants. A large amount of experimental data is available on blue copper proteins; however, their functional characterization is hindered by the complexity of redox processes in biological systems. We describe here the application of a semiquantitative method based on a comparative analysis of molecular interaction fields to gain insights into the recognition properties of blue copper proteins. Molecular electrostatic and hydrophobic potentials were computed and compared for a set of 33 experimentally-determined structures of proteins from seven blue copper subfamilies, and the results were quantified by means of similarity indices. The analysis provides a classification of the blue copper proteins and shows that (I) comparison of the molecular electrostatic potentials provides useful information complementary to that highlighted by sequence analysis; (2) similarities in recognition properties can be detected for proteins belonging to different subfamilies, such as amicyanins and pseudoazurins, that may be isofunctional proteins; (3) dissimilarities in interaction properties, consistent with experimentally different binding specificities, may be observed between proteins belonging to the same subfamily, such as cyanobacterial and eukaryotic plastocyanins; (4) proteins with low sequence identity, such as azurins and pseudoazurins, can have sufficient similarity to bind to similar electron donors and acceptors while having different binding specificity profiles.

Tyr(30) of amicyanin is not critical for electron transfer to cytochrome c-551i: implications for predicting electron transfer pathways

A Pathways analysis of the methylamine dehydrogenase-amicyanin-cytochrome c-551i protein electron transfer (ET) complex predicts two sets of ET pathways of comparable efficiency from the type I copper of amicyanin to the heme of cytochrome c-551i. In one pathway, the electron exits copper via the Cys(92) copper ligand, and in the other, it exits via the Met(98) copper ligand. If the Pathways algorithm is modified to include contributions from the anisotropy of metal-ligand coupling, independent of differences in copper-ligand bond length, then the pathways via Cys(92) are predicted to be at least 100-fold more strongly coupled than the pathways via any of the other copper ligands. All of the favored pathways via Cys(92) include a through-space jump from Cys(92) to the side chain of Tyr(30). To determine whether or not the pathways via Cys(92) are preferentially used for ET, Tyr(30) was changed to other amino acid residues by site-directed mutagenesis. Some mutant proteins were very unstable suggesting a role for Tyr(30) in stabilizing the protein structure. Y30F and Y30I mutant amicyanins could be isolated and analyzed. For the Y30I mutant, the modified Pathways analysis which favors ET via Cys(92) predicts a decrease in ET rate of at least two orders of magnitude, whereas the standard Pathways analysis predicts no change in ET rate since ET via Met(98) is not affected. Experimentally, the ET rates of the Y30I and Y30F mutants were indistinguishable from that of wild-type amicyanin. Likely explanations for these observations are discussed as are their implications for predicting pathways for ET reactions of metalloproteins.