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

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.

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.

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.

Detecting transitions in protein dynamics using a recurrence quantification analysis based bootstrap method

Background

Proteins undergo conformational transitions over different time scales. These transitions are closely intertwined with the protein’s function. Numerous standard techniques such as principal component analysis are used to detect these transitions in molecular dynamics simulations. In this work, we add a new method that has the ability to detect transitions in dynamics based on the recurrences in the dynamical system. It combines bootstrapping and recurrence quantification analysis. We start from the assumption that a protein has a “baseline” recurrence structure over a given period of time. Any statistically significant deviation from this recurrence structure, as inferred from complexity measures provided by recurrence quantification analysis, is considered a transition in the dynamics of the protein.

Results

We apply this technique to a 132 ns long molecular dynamics simulation of the β-Lactamase Inhibitory Protein BLIP. We are able to detect conformational transitions in the nanosecond range in the recurrence dynamics of the BLIP protein during the simulation. The results compare favorably to those extracted using the principal component analysis technique.

Conclusions

The recurrence quantification analysis based bootstrap technique is able to detect transitions between different dynamics states for a protein over different time scales. It is not limited to linear dynamics regimes, and can be generalized to any time scale. It also has the potential to be used to cluster frames in molecular dynamics trajectories according to the nature of their recurrence dynamics. One shortcoming for this method is the need to have large enough time windows to insure good statistical quality for the recurrence complexity measures needed to detect the transitions.

Effect of circular permutation on the structure and function of type 1 blue copper center in azurin

Type 1 copper (T1Cu) proteins are electron transfer (ET) proteins involved in many important biological processes. While the effects of changing primary and secondary coordination spheres in the T1Cu ET function have been extensively studied, few report has explored the effect of the overall protein structural perturbation on active site configuration or reduction potential of the protein, even though the protein scaffold has been proposed to play a critical role in enforcing the entatic or "rack-induced" state for ET functions. We herein report circular permutation of azurin by linking the N- and C-termini and creating new termini in the loops between 1st and 2nd β strands or between 3rd and 4th β strands. Characterization by electronic absorption, electron paramagnetic spectroscopies, as well as crystallography and cyclic voltammetry revealed that, while the overall structure and the primary coordination sphere of the circular permutated azurins remain the same as those of native azurin, their reduction potentials increased by 18 and 124 mV over that of WTAz. Such increases in reduction potentials can be attributed to subtle differences in the hydrogen-bonding network in secondary coordination sphere around the T1Cu center.

Engineering the metal-binding loop at a type 1 copper center by circular permutation

Circular permutation of the cupredoxin azurin creates a break on the metal binding loop, highlighting the loop's flexibility.

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.

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

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

How the dynamics of the metal-binding loop region controls the acid transition in cupredoxins

Many reduced cupredoxins undergo a pH-dependent structural rearrangement, triggered by protonation of the His ligand belonging to the C-terminal hydrophobic loop, usually termed the acid transition. At variance with several members of the cupredoxin family, the acid transition is not observed for azurin (AZ). We have addressed this issue by performing molecular dynamics simulations of AZ and four mutants, in which the C-terminal loop has been replaced with those of other cupredoxins or with polyalanine loops. All of the loop mutants undergo the acid transition in the pH range of 4.4-5.5. The main differences between AZ and its loop mutants are the average value of the active site solvent accessible surface area and the extent of its fluctuations with time, together with an altered structure of the water layer around the copper center. Using functional mode analysis, we found that these variations arise from changes in nonbonding interactions in the second coordination sphere of the copper center, resulting from the loop mutation. Our results strengthen the view that the dynamics at the site relevant for function and its surroundings are crucial for protein activity and for metal-containing electron transferases.

HHM motif at the CuH-site of peptidylglycine monooxygenase is a pH-dependent conformational switch

Peptidylglycine monooxygenase is a copper-containing enzyme that catalyzes the amidation of neuropeptides hormones, the first step of which is the conversion of a glycine-extended pro-peptide to its α-hydroxyglcine intermediate. The enzyme contains two mononuclear Cu centers termed CuM (ligated to imidazole nitrogens of H242, H244 and the thioether S of M314) and CuH (ligated to imidazole nitrogens of H107, H108, and H172) with a Cu-Cu separation of 11 Å. During catalysis, the M site binds oxygen and substrate, and the H site donates the second electron required for hydroxylation. The WT enzyme shows maximum catalytic activity at pH 5.8 and undergoes loss of activity at lower pHs due to a protonation event with a pKA of 4.6. Low pH also causes a unique structural transition in which a new S ligand coordinates to copper with an identical pKA, manifest by a large increase in Cu-S intensity in the X- ray absorption spectroscopy. In previous work (Bauman, A. T., Broers, B. A., Kline, C. D., and Blackburn, N. J. (2011) Biochemistry 50, 10819-10828), we tentatively assigned the new Cu-S interaction to binding of M109 to the H-site (part of an HHM conserved motif common to all but one member of the family). Here we follow up on these findings via studies on the catalytic activity, pH-activity profiles, and spectroscopic (electron paramagnetic resonance, XAS, and Fourier transform infrared) properties of a number of H-site variants, including H107A, H108A, H172A, and M109I. Our results establish that M109 is indeed the coordinating ligand and confirm the prediction that the low pH structural transition with associated loss of activity is abrogated when the M109 thioether is absent. The histidine mutants show more complex behavior, but the almost complete lack of activity in all three variants coupled with only minor differences in their spectroscopic properties suggests that unique structural elements at H are critical for functionality. The data suggest a more general utility for the HHM motif as a copper- and pH-dependent conformational switch.

The Active Site Loop Modulates the Reorganization Energy of Blue Copper Proteins by Controlling the Dynamic Interplay with Solvent

Understanding the factors governing the rate of electron transfer processes in proteins is crucial not only to a deeper understanding of redox processes in living organisms but also for the design of efficient devices featuring biological molecules. Here, molecular dynamics simulations performed on native azurin and four chimeric cupredoxins allow for the calculation of the reorganization energy and of structure-related quantities that were used to clarify the molecular determinants to the dynamics/function relationship in blue copper proteins. We find that the dynamics of the small, metal-binding loop region controls the outer-sphere reorganization energy not only by determining the exposure of the active site to solvent but also through the modulation of the redox-dependent rearrangement of the whole protein scaffold and of the surrounding water molecules.

The influence of protein folding on the copper affinities of trafficking and target sites

The relative influence of protein unfolding on the Cu(I) affinity of trafficking and target sites for copper has been determined. For the copper metallochaperone Atx1 from Synechocystis PCC 6803 (a cyanobacterium), Saccharomyces cerevisiae and humans unfolding in urea results in a decrease in the Cu(I) affinity from (4-5) × 10(17) M(-1) to (1-3) × 10(16) M(-1) at pH 7. The affinities of the unfolded Atx1s are similar to those for CXXC-containing peptides. Partial unfolding, due to the loop 5 His61Lys mutation in Synechocystis Atx1, gives rise to a more limited decrease in Cu(I) affinity. For the copper target protein plastocyanin from Synechocystis, chemical unfolding results in the Cu(I) affinity decreasing by 5-orders of magnitude. This differential influence of protein unfolding on Cu(I) affinity is due to a more complex copper site structure in the target protein, including numerous interactions of non-coordinating residues with ligating amino acids. This second-coordination sphere is much simpler in the Atx1s with the main interaction provided by the loop 5 residue that tunes the Cu(I) affinity by altering the pK(a) of the C-terminal Cys ligand of the CXXC motif. This interaction and others are absent in the unfolded Atx1s and the two Cys ligands have pK(a) values reminiscent of free thiols (>8) resulting in lowered Cu(I) affinities at pH 7. Residues close to the active site of the thiol-disulfide oxidoreductase thioredoxin appear to lower the Cu(I) affinity of its CXXC motif to 3.1 × 10(15) M(-1) at pH 7, presumably to prevent copper binding in vivo. The structure of a copper site, including the number and relative position of ligands in the primary structure and the complexity of the second-coordination sphere, results in dramatically different effects of unfolding on Cu(I) affinity that has important implications for copper homeostasis.

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.

The influence of active site loop mutations on the thermal stability of azurin from Pseudomonas aeruginosa

The copper site and overall structures of azurin (AZ) variants in which the amicyanin (AMI) and plastocyanin (PC) metal binding loops have been introduced, AZAMI and AZPC, respectively, are similar to that of AZ, whereas the loop conformations resemble those in the native proteins. To assess the influence of these loop mutations on stability, the thermal unfolding of AZAMI and AZPC has been investigated by differential scanning calorimetry, absorption and fluorescence spectroscopy. The calorimetric profiles of both variants exhibit a complex shape consisting of two endothermic peaks and an exothermic peak. The temperature of the maximum heat of absorption for the single endothermic peak is 82.7°C for AZ, whereas for AZAMI and AZPC the most intense endothermic peaks are at 74.9 and 68.1°C comparable to values for AMI and PC, respectively. Denaturation investigated using the temperature dependence of the absorbance at ∼600nm and Trp emission, also demonstrates decreased stability for both loop mutants. The thermal transition between the native and the denaturated states is irreversible, scan rate dependent and consistent with the two-state irreversible model. The structure of the active-site loop has a dramatic effect on the kinetic stability and the unfolding pathway of cupredoxins.

The importance of slow motions for protein functional loops

Loops in proteins that connect secondary structures such as alpha-helix and beta-sheet, are often on the surface and may play a critical role in some functions of a protein. The mobility of loops is central for the motional freedom and flexibility requirements of active-site loops and may play a critical role for some functions. The structures and behaviors of loops have not been studied much in the context of the whole structure and its overall motions, especially how these might be coupled. Here we investigate loop motions by using coarse-grained structures (C(α) atoms only) to solve the motions of the system by applying Lagrange equations with elastic network models to learn about which loops move in an independent fashion and which move in coordination with domain motions, faster and slower, respectively. The normal modes of the system are calculated using eigen-decomposition of the stiffness matrix. The contribution of individual modes and groups of modes is investigated for their effects on all residues in each loop by using Fourier analyses. Our results indicate overall that the motions of functional sets of loops behave in similar ways as the whole structure. But overall only a relatively few loops move in coordination with the dominant slow modes of motion, and these are often closely related to function.

A copper-methionine interaction controls the pH-dependent activation of peptidylglycine monooxygenase

The pH dependence of native peptidylglycine monooxygenase (PHM) and its M314H variant has been studied in detail. For wild-type (WT) PHM, the intensity of the Cu-S interaction visible in the Cu(I) extended X-ray absorption fine structure (EXAFS) data is inversely proportional to catalytic activity over the pH range of 3-8. A previous model based on more limited data was interpreted in terms of two protein conformations involving an inactive Met-on form and an active flexible Met-off form [Bauman, A. T., et al. (2006) Biochemistry 45, 11140-11150] that derived its catalytic activity from the ability to couple into vibrational modes critical for proton tunneling. The new studies comparing the WT and M314H variant have led to the evolution of this model, in which the Met-on form has been found to be derived from coordination of an additional Met residue, rather than a more rigid conformer of M314 as previously proposed. The catalytic activity of the mutant decreased by 96% because of effects on both k(cat) and K(M), but it displayed the same activity-pH profile with a maximum around pH 6. At pH 8, the reduced Cu(I) form gave spectra that could be simulated by replacement of the Cu(M) Cu-S(Met) interaction with a Cu-N/O interaction, but the data did not unambiguously assign the ligand to the imidazole side chain of H314. At pH 3.5, the EXAFS still showed the presence of a strong Cu-S interaction, establishing that the Met-on form observed at low pH in WT cannot be due to a strengthening of the Cu(M)-methionine interaction but must arise from a different Cu-S interaction. Therefore, lowering the pH causes a conformational change at one of the Cu centers that brings a new S donor residue into a favorable orientation for coordination to copper and generates an inactive form. Cys coordination is unlikely because all Cys residues in PHM are engaged in disulfide cross-links. Sequence comparison with the PHM homologues tyramine β-monooxygenase and dopamine β-monooxygenase suggests that M109 (adjacent to H site ligands H107 and H108) is the most likely candidate. A model is presented in which H108 is protonated with a pK(a) of 4.6 to generate the inactive low-pH form with Cu(H) coordinated by M109, H107, and H172.

Models to Approximate the Motions of Protein Loops

We approximate the loop motions of various proteins by using a coarse-grained model and the theory of rubberlike elasticity of polymer chains. The loops are considered as chains where only the first and the last residues thereof are tethered by their connections to the main structure; while within the loop, the loop residues are connected only to their sequence neighbors. We applied these approximate models to five proteins. Our approximation shows that the loop motions can usually be computed locally which shows these motions are robust and not random. But most interestingly, the new method presented here can be used to compute the likely motions of loops that are missing in the structures.

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.

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

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

The crystal structure of auracyanin A at 1.85 A resolution: the structures and functions of auracyanins A and B, two almost identical "blue" copper proteins, in the photosynthetic bacterium Chloroflexus aurantiacus

Auracyanins A and B are two closely similar "blue" copper proteins produced by the filamentous anoxygenic phototrophic bacterium Chloroflexus aurantiacus. Both proteins have a water-soluble 140-residue globular domain, which is preceded in the sequence by an N-terminal tail. The globular domains of auracyanins A and B have sequences that are 38% identical. The sequences of the N-terminal tails, on the other hand, are distinctly different, suggesting that auracyanins A and B occupy different membrane sites and have different functions. The crystal structure of auracyanin A has been solved and refined at 1.85 A resolution. The polypeptide fold is similar to that of auracyanin B (Bond et al. in J Mol Biol 306:47-67, 2001), but the distribution of charged and polar residues on the molecular surface is different. The Cu-site dimensions of the two auracyanins are identical. This is unexpected, since auracyanin A has a shorter polypeptide loop between two of the Cu-binding residues, and the two proteins have significantly different EPR, UV-visible and resonance Raman spectra. The genes for the globular domains of auracyanins A and B have been cloned in a bacterial expression system, enabling purification of large quantities of protein. It is shown that auracyanin A is expressed only when C. aurantiacus cells are grown in light, whereas auracyanin B is expressed under dark as well as light conditions. The inference is that auracyanin A has a function in photosynthesis, and that auracyanin B has a function in aerobic respiration.

Conformational exchange in pseudoazurin: different kinds of microsecond to millisecond dynamics characterized by their pH and buffer dependence using 15N NMR relaxation

The dynamics of the reduced form of the blue copper protein pseudoazurin from Alcaligenes faecalis S-6 was investigated using (15)N relaxation measurements with a focus on the dynamics of the micro- to millisecond time scale. Different types of conformational exchange processes are observed in the protein on this time scale. At low pH, the protonation of the C-terminal copper-ligated histidine, His81, is observed. A comparison of the exchange rates in the presence and absence of added buffers shows that the protonation is the rate-limiting step at low buffer concentrations. This finding agrees with previous observations for other blue copper proteins, e.g., amicyanin and plastocyanin. However, in contrast to plastocyanin but similar to amicyanin, a second conformational exchange between different conformations of the protonated copper site is observed at low pH, most likely triggered by the protonation of His81. This process has been further characterized using CPMG dispersion methods and is found to occur with a rate of a few thousands per second. Finally, micro- to millisecond motions are observed in one of the loop regions and in the alpha-helical regions. These motions are unaffected by pH and are unrelated to the conformational changes in the active site of pseudoazurin.

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.