Metal-Binding Loop Length Is a Determinant of the pKa of a Histidine Ligand at a Type 1 Copper Site

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.

DFTB3 Parametrization for Copper: The Importance of Orbital Angular Momentum Dependence of Hubbard Parameters

We report the parametrization of a density functional tight binding method (DFTB3) for copper in a spin-polarized formulation. The parametrization is consistent with the framework of 3OB for main group elements (ONCHPS) and can be readily used for biological applications that involve copper proteins/peptides. The key to our parametrization is to introduce orbital angular momentum dependence of the Hubbard parameter and its charge derivative, thus allowing the 3d and 4s orbitals to adopt different sizes and responses to the change of charge state. The parametrization has been tested by applying to a fairly broad set of molecules of biological relevance, and the properties of interest include optimized geometries, ligand binding energies, and ligand proton affinities. Compared to the reference QM level (B3LYP/aug-cc-pVTZ, which is shown here to be similar to the B97-1 and CCSD(T) results, in terms of many properties of interest for a set of small copper containing molecules), our parametrization generally gives reliable structural properties for both Cu(I) and Cu(II) compounds, although several exceptions are also noted. For energetics, the results are more accurate for neutral ligands than for charged ligands, likely reflecting the minimal basis limitation of DFTB3; the results generally outperform NDDO based methods such as PM6 and even PBE with the 6-31+G(d,p) basis. For all ligand types, single-point B3LYP calculations at DFTB3 geometries give results very close (∼1-2 kcal/mol) to the reference B3LYP values, highlighting the consistency between DFTB3 and B3LYP structures. Possible further developments of the DFTB3 model for a better treatment of transition-metal ions are also discussed. In the current form, our first generation of DFTB3 copper model is expected to be particularly valuable as a method that drives sampling in systems that feature a dynamical copper binding site.

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.

A Euclidean Perspective on the Unfolding of Azurin: Spatial Correlations

We investigate the stability to structural perturbation of Pseudomonas aeruginosa azurin using a previously developed geometric model. Our analysis considers Ru(2,2',6',2″-terpyridine)(1,10-phenanthroline)(His83)-labeled wild-type azurin and five variants with mutations to Cu-ligating residues. We find that in the early stages of unfolding, the β-strands exhibit the most structural stability. The conserved residues comprising the hydrophobic core are dislocated only after nearly complete unfolding of the β-barrel. Attachment of the Ru-complex at His83 does not destabilize the protein fold, despite causing some degree of structural rearrangement. Notably, replacing the Cys112 and/or Met121 Cu ligands does not affect the conformational integrity of the protein. Notably, these results are in accord with experimental evidence, as well as molecular dynamics simulations of the denaturation of azurin.

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.

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.