Spectroscopic and density functional theory studies of the blue-copper site in M121SeM and C112SeC azurin: Cu-Se versus Cu-S bonding

Theoretical insights into the reduction of Azurin metal site with unnatural amino acid substitutions

Copper-containing proteins play crucial roles in biological systems. Azurin is a copper-containing protein which has a Type 1 copper site that facilitates electron transfer in the cytochrome chain. Previous research has highlighted the significant impact of mutations in the axial Met121 of the copper site on the reduction potential. However, the mechanism of this regulation has not been fully established. In this study, we employed theoretical modeling to investigate the reduction of the Type 1 copper site, focusing on how unnatural amino acid substitutions at Met121 influence its behavior. Our findings demonstrated a strong linear correlation between electrostatic interactions and the reduction potential of the copper site, which indicates that the perturbation of the reduction potential is primarily influenced by electrostatic interactions between the metal ion and the ligating atom. Furthermore, we found that CF/π and CF…H interactions could induce subtle changes in geometry and hence impact the electronic properties of the systems under study. In addition, our calculations suggest the coordination mode and ion-ligand distance could significantly impact the reduction potential of a copper site. Overall, this study offers valuable insights into the structural and electronic properties of the Type 1 copper site, which could potentially guide the design of future artificial catalysts.

Long-Range Charge Delocalization Mediates the Ultrafast Ligand-to-Metal Charge Transfer Dynamics at the Cu2+-Active Site in Azurin

The blue color metalloprotein in azurin has traditionally been attributed to the intense cysteine-to-Cu²⁺ ligand-to-metal charge transfer transition centered at 628 nm. Although resonance Raman measurements of the Cu²⁺ active site have implied that the LMCT transition electronically couples to the protein scaffold well beyond its primary metal-ligand coordination shell, the structural extent of this electronic coupling and visualization of the protein-mediated charge transfer dynamics have remained elusive. Here, using femtosecond broadband transient absorption and impulsive Raman spectroscopy, we provide direct evidence for a rapid relaxation between two distinct charge transfer states, having different spatial delocalization, within ∼300 fs followed by recombination of charges in subpicosecond time scales. We invoke the formation of a protein-centered radical cation, possibly Trp48 or a Phe residue, within 100 fs substantiating the long-range electronic coupling for the first time beyond the traditional copper active site. The Raman spectra of the excited CT state show the presence of protein-centric vibrations along with the vibrational modes assigned to the copper active site. Our results demonstrate a large delocalization length scale of the initially populated CT state, thereby highlighting the possibility of exploiting azurin photochemistry for energy conversion techniques.

Converged Structural and Spectroscopic Properties for Refined QM/MM Models of Azurin

Blue copper proteins continue to challenge experiment and theory with their electronic structure and spectroscopic properties that respond sensitively to the coordination environment of the copper ion. In this work, we report state-of-the art electronic structure studies for geometric and spectroscopic properties of the archetypal "Type I" copper protein azurin in its Cu(II) state. A hybrid quantum mechanics/molecular mechanics (QM/MM) approach is used, employing both density functional theory (DFT) and coupled cluster with singles, doubles, and perturbative triples (CCSD(T)) methods for the QM region, the latter method making use of the domain-based local pair natural orbital (DLPNO) approach. Models of increasing QM size are employed to investigate the convergence of critical geometric parameters. It is shown that convergence is slow and that a large QM region is critical for reproducing the short experimental Cu-SCys112 distance. The study of structural convergence is followed by investigation of spectroscopic parameters using both DFT and DLPNO-CC methods and comparing these to the experimental spectrum using simulations. The results allow us to examine for the first time the distribution of spin densities and hyperfine coupling constants at the coupled cluster level, leading us to revisit the experimental assignment of the 33S hyperfine splitting. The wavefunction-based approach to obtain spin-dependent properties of open-shell systems demonstrated here for the case of azurin is transferable and applicable to a large array of bioinorganic systems.

Structural Basis for a Quadratic Relationship between Electronic Absorption and Electronic Paramagnetic Resonance Parameters of Type 1 Copper Proteins

Type 1 copper (T1Cu) proteins play important roles in electron transfer in biology, largely due to the unique structure of the T1Cu center, which is reflected by its spectroscopic properties. Previous reports have suggested a correlation between a high ratio of electronic absorbance at ∼450 nm to that at ∼600 nm (R = A₄₅₀/A₆₀₀) and a large copper(II) hyperfine coupling in the z direction (A_z) in electron paramagnetic resonance (EPR). However, this correlation does not have a clear physical meaning, nor does it hold for many proteins with a perturbed T1Cu center. To address this issue, a new parameter of R' [A₄₅₀/(A₄₅₀ + A₆₀₀)] with a better physical meaning of a fractional S_Cys pseudo-σ to Cu(II) charge transfer transition intensity is defined and a quadratic relationship between R' and A_z is found on the basis of a comprehensive analysis of ultraviolet-visible absorption, EPR, and structural parameters of T1Cu proteins. We are able to find good correlations between R' and the displacement of copper from the trigonal plane defined by the His₂Cys ligands and the angle between the N_His1-Cu-N_His2 plane and the S_Cys-Cu-axial ligand plane, providing a structural basis for the observed correlation. These findings and analyses provide a new framework for a deeper understanding of the spectroscopic and electronic properties of T1Cu proteins, which may allow better design and applications of this important class of proteins for redox and electron transfer functions.

Understanding and Modulating Metalloenzymes with Unnatural Amino Acids, Non-Native Metal Ions, and Non-Native Metallocofactors

Metalloproteins set the gold standard for performing important functions, including catalyzing demanding reactions under mild conditions. Designing artificial metalloenzymes (ArMs) to catalyze abiological reactions has been a major endeavor for many years, but most ArM activities are far below those of native enzymes, making them unsuitable for most pratical applications. A critical step to advance the field is to fundamentally understand what it takes to not only confer but also fine-tune ArM activities so they match those of native enzymes. Indeed, only once we can freely modulate ArM activity to rival (or surpass!) natural enzymes can the potential of ArMs be fully realized. A key to unlocking ArM potential is the observation that one metal primary coordination sphere can display a range of functions and levels of activity, leading to the realization that secondary coordination sphere (SCS) interactions are critically important. However, SCS interactions are numerous, long-range, and weak, making them very difficult to reproduce in ArMs. Furthermore, natural enzymes are tied to a small set of biologically available functional moieties from canonical amino acids and physiologically available metal ions and metallocofactors, severely limiting the chemical space available to probe and tune ArMs. In this Account, we summarize the use of unnatural amino acids (UAAs) and non-native metal ions and metallocofactors by our group and our collaborators to probe and modulate ArM functions. We incorporated isostructural UAAs in a type 1 copper (T1Cu) protein azurin to provide conclusive evidence that axial ligand hydrophobicity is a major determinant of T1Cu redunction potential ( E°'). Closely related work from other groups are also discussed. We also probed the role of protein backbone interactions that cannot be altered by standard mutagenesis by replacing the peptide bond with an ester linkage. We used insight gained from these studies to tune the E°' of azurin across the entire physiological range, the broadest range ever achieved in a single metalloprotein. Introducing UAA analogues of Tyr into ArM models of heme-copper oxidase (HCO) revealed a linear relationship between p K_a, E°', and activity. We also substituted non-native hemes and non-native metal ions for their native equivalents in these models to resolve several issues that were intractable in native HCOs and the closely related nitric oxide reductases, such as their roles in modulating substrate affinity, electron transfer rate, and activity. We incorporated abiological cofactors such as ferrocene and Mn(salen) into azurin and myoglobin, respectively, to stabilize these inorganic and organometallic compounds in water, confer abiological functions, tune their E°' and activity through SCS interactions, and show that the approach to metallocofactor anchoring and orientation can tune enantioselectivity and alter function. Replacing Cu in azurin with non-native Fe or Ni can impart novel activities, such as superoxide reduction and C-C bond formation. While progress was made, we have identified only a small fraction of the interactions that can be generally applied to ArMs to fine-tune their functions. Because SCS interactions are subtle and heavily interconnected, it has been difficult to characterize their effects quantitatively. It is vital to develop spectroscopic and computational techniques to detect and quantify their effects in both resting states and catalytic intermediates.

Ab initio electronic structure calculations of entire blue copper azurins

We present a theoretical study of the blue-copper azurin extracted from Pseudomonas aeruginosa and several of its single amino acid mutants.

Design of Heteronuclear Metalloenzymes

Heteronuclear metalloenzymes catalyze some of the most fundamentally interesting and practically useful reactions in nature. However, the presence of two or more metal ions in close proximity in these enzymes makes them more difficult to prepare and study than homonuclear metalloenzymes. To meet these challenges, heteronuclear metal centers have been designed into small and stable proteins with rigid scaffolds to understand how these heteronuclear centers are constructed and the mechanism of their function. This chapter describes methods for designing heterobinuclear metal centers in a protein scaffold by giving specific examples of a few heme-nonheme bimetallic centers engineered in myoglobin and cytochrome c peroxidase. We provide step-by-step procedures on how to choose the protein scaffold, design a heterobinuclear metal center in the protein scaffold computationally, incorporate metal ions into the protein, and characterize the resulting metalloproteins, both structurally and functionally. Finally, we discuss how an initial design can be further improved by rationally tuning its secondary coordination sphere, electron/proton transfer rates, and the substrate affinity.

The mechanochemistry of copper reports on the directionality of unfolding in model cupredoxin proteins

Understanding the directionality and sequence of protein unfolding is crucial to elucidate the underlying folding free energy landscape. An extra layer of complexity is added in metalloproteins, where a metal cofactor participates in the correct, functional fold of the protein. However, the precise mechanisms by which organometallic interactions are dynamically broken and reformed on (un)folding are largely unknown. Here we use single molecule force spectroscopy AFM combined with protein engineering and MD simulations to study the individual unfolding pathways of the blue-copper proteins azurin and plastocyanin. Using the nanomechanical properties of the native copper centre as a structurally embedded molecular reporter, we demonstrate that both proteins unfold via two independent, competing pathways. Our results provide experimental evidence of a novel kinetic partitioning scenario whereby the protein can stochastically unfold through two distinct main transition states placed at the N and C termini that dictate the direction in which unfolding occurs.

Using solution- and solid-state S K-edge X-ray absorption spectroscopy with density functional theory to evaluate M-S bonding for MS4(2-) (M = Cr, Mo, W) dianions

Herein, we have evaluated relative changes in M-S electronic structure and orbital mixing in Group 6 MS4(2-) dianions using solid- and solution-phase S K-edge X-ray absorption spectroscopy (XAS; M = Mo, W), as well as density functional theory (DFT; M = Cr, Mo, W) and time-dependent density functional theory (TDDFT) calculations. To facilitate comparison with solution measurements (conducted in acetonitrile), theoretical models included gas-phase calculations as well as those that incorporated an acetonitrile dielectric, the latter of which provided better agreement with experiment. Two pre-edge features arising from S 1s → e* and t electron excitations were observed in the S K-edge XAS spectra and were reasonably assigned as (1)A1 → (1)T2 transitions. For MoS4(2-), both solution-phase pre-edge peak intensities were consistent with results from the solid-state spectra. For WS4(2-), solution- and solid-state pre-edge peak intensities for transitions involving e* were equivalent, while transitions involving the t orbitals were less intense in solution. Experimental and computational results have been presented in comparison to recent analyses of MO4(2-) dianions, which allowed M-S and M-O orbital mixing to be evaluated as the principle quantum number (n) for the metal valence d orbitals increased (3d, 4d, 5d). Overall, the M-E (E = O, S) analyses revealed distinct trends in orbital mixing. For example, as the Group 6 triad was descended, e* (π*) orbital mixing remained constant in the M-S bonds, but increased appreciably for M-O interactions. For the t orbitals (σ* + π*), mixing decreased slightly for M-S bonding and increased only slightly for the M-O interactions. These results suggested that the metal and ligand valence orbital energies and radial extensions delicately influenced the orbital compositions for isoelectronic ME4(2-) (E = O, S) dianions.

Importance of polarization effect in the study of metalloproteins: application of polarized protein specific charge scheme in predicting the reduction potential of azurin

Molecular dynamics (MD) simulation is commonly used in the study of protein dynamics, and in recent years, the extension of MD simulation to the study of metalloproteins is gaining much interest. Choice of force field is crucial in MD studies, and the inclusion of metal centers complicates the process of accurately describing the electrostatic environment that surrounds the redox centre. Herein, we would like to explore the importance of including electrostatic contribution from both protein and solvent in the study of metalloproteins. MD simulations with the implementation of thermodynamic integration will be conducted to model the reduction process of azurin from Pseudomonas aeruginosa. Three charge schemes will be used to derive the partial charges of azurin. These charge schemes differ in terms of the amount of immediate environment, respective to copper, considered during charge fitting, which ranges from the inclusion of copper and residues in the first coordination sphere during density functional theory charge fitting to the comprehensive inclusion of protein and solvent effect surrounding the metal centre using polarized protein-specific charge scheme. From the simulations conducted, the relative reduction potential of the mutated azurins respective to that of wild-type azurin (ΔEcal) were calculated and compared with experimental values. The ΔEcal approached experimental value with increasing consideration of environmental effect hence substantiating the importance of polarization effect in the study of metalloproteins. This study also attests the practicality of polarized protein-specific charge as a computational tool capable of incorporating both protein environment and solvent effect into MD simulations.

Positive and negative ion formation in deep-core excited molecules: S 1s excitation in dimethyl sulfoxide

The photo-fragmentation of the dimethyl sulfoxide (DMSO) molecule was studied using synchrotron radiation and a magnetic mass spectrometer. The total cationic yield spectrum was recorded in the photon energy region around the sulfur K edge. The sulfur composition of the highest occupied molecular orbital's and lowest unoccupied molecular orbital's in the DMSO molecule has been obtained using both ab initio and density functional theory methods. Partial cation and anion-yield measurements were obtained in the same energy range. An intense resonance is observed at 2475.4 eV. Sulfur atomic ions present a richer structure around this resonant feature, as compared to other fragment ions. The yield curves are similar for most of the other ionic species, which we interpret as due to cascade Auger processes leading to multiply charged species which then undergo Coulomb explosion. The anions S(-), C(-), and O(-) are observed for the first time in deep-core-level excitation of DMSO.

Spectroscopic and DFT studies of second-sphere variants of the type 1 copper site in azurin: covalent and nonlocal electrostatic contributions to reduction potentials

The reduction potentials (E(0)) of type 1 (T1) or blue copper (BC) sites in proteins and enzymes with identical first coordination spheres around the redox active copper ion can vary by ~400 mV. Here, we use a combination of low-temperature electronic absorption and magnetic circular dichroism, electron paramagnetic resonance, resonance Raman, and S K-edge X-ray absorption spectroscopies to investigate a series of second-sphere variants--F114P, N47S, and F114N in Pseudomonas aeruginosa azurin--which modulate hydrogen bonding to and protein-derived dipoles nearby the Cu-S(Cys) bond. Density functional theory calculations correlated to the experimental data allow for the fractionation of the contributions to tuning E(0) into covalent and nonlocal electrostatic components. These are found to be significant, comparable in magnitude, and additive for active H-bonds, while passive H-bonds are mostly nonlocal electrostatic in nature. For dipoles, these terms can be additive to or oppose one another. This study provides a methodology for uncoupling covalency from nonlocal electrostatics, which, when coupled to X-ray crystallographic data, distinguishes specific local interactions from more long-range protein/active interactions, while affording further insight into the second-sphere mechanisms available to the protein to tune the E(0) of electron-transfer sites in biology.

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.

Analysis of resonance Raman data on the blue copper site in pseudoazurin: excited state π and σ charge transfer distortions and their relation to ground state reorganization energy

The short Cu(2+)-S(Met) bond in pseudoazurin (PAz) results in the presence of two relatively intense S(p)(π) and S(p)(σ) charge transfer (CT) transitions. This has enabled resonance Raman (rR) data to be obtained for each excited state. The rR data show very different intensity distribution patterns for the vibrations in the 300-500 cm(-1) region. Time-dependent density functional theory (TDDFT) calculations have been used to determine that the change in intensity distribution between the S(p)(π) and S(p)(σ) excited states reflects the differential enhancement of S(Cys) backbone modes with Cu-S(Cys)-C(β) out-of-plane (oop) and in-plane (ip) bend character in their respective potential energy distributions (PEDs). The rR excited state distortions have been related to ground state reorganization energies (λ s) and predict that, in addition to M-L stretches, the Cu-S(Cys)-C(β) oop bend needs to be considered. DFT calculations predict a large distortion in the Cu-S(Cys)-C(β) oop bending coordinate upon reduction of a blue copper (BC) site; however, this distortion is not present in the X-ray crystal structures of reduced BC sites. The lack of Cu-S(Cys)-C(β) oop distortion upon reduction corresponds to a previously unconsidered constraint on the thiolate ligand orientation in the reduced state of BC proteins and can be considered as a contribution to the entatic/rack nature of BC sites.

Crucial role of copper in detection of metal-coordinating odorants

Odorant receptors (ORs) in olfactory sensory neurons (OSNs) mediate detection of volatile odorants. Divalent sulfur compounds, such as thiols and thioethers, are extremely potent odorants. We identify a mouse OR, MOR244-3, robustly responding to (methylthio)methanethiol (MeSCH(2)SH; MTMT) in heterologous cells. Found specifically in male mouse urine, strong-smelling MTMT [human threshold 100 parts per billion (ppb)] is a semiochemical that attracts female mice. Nonadjacent thiol and thioether groups in MTMT suggest involvement of a chelated metal complex in MOR244-3 activation. Metal ion involvement in thiol-OR interactions was previously proposed, but whether these ions change thiol-mediated OR activation remained unknown. We show that copper ion among all metal ions tested is required for robust activation of MOR244-3 toward ppb levels of MTMT, structurally related sulfur compounds, and other metal-coordinating odorants (e.g., strong-smelling trans-cyclooctene) among >125 compounds tested. Copper chelator (tetraethylenepentamine, TEPA) addition abolishes the response of MOR244-3 to MTMT. Histidine 105, located in the third transmembrane domain near the extracellular side, is proposed to serve as a copper-coordinating residue mediating interaction with the MTMT-copper complex. Electrophysiological recordings of the OSNs in the septal organ, abundantly expressing MOR244-3, revealed neurons responding to MTMT. Addition of copper ion and chelator TEPA respectively enhanced and reduced the response of some MTMT-responding neurons, demonstrating the physiological relevance of copper ion in olfaction. In a behavioral context, an olfactory discrimination assay showed that mice injected with TEPA failed to discriminate MTMT. This report establishes the role of metal ions in mammalian odor detection by ORs.

Valence-to-core X-ray emission spectroscopy: a sensitive probe of the nature of a bound ligand

The sensitivity of iron Kβ X-ray emission spectroscopy (XES) to the nature of the bound ligands (σ-donating, π-donating, and π-accepting) has been explored. A combination of experiment and theory has been employed, with a DFT approach being utilized to elucidate ligand effects on the spectra and to assign the spectral intensity mechanisms. Knowledge of the various contributions to the spectra allows for a deeper understanding of spectral features and demonstrates the sensitivity of this method to the identity of the interacting ligands. The potential of XES for identifying intermediate species in nonheme iron enzymes is highlighted.

The essential role of the Cu(II) state of Sco in the maturation of the Cu(A) center of cytochrome oxidase: evidence from H135Met and H135SeM variants of the Bacillus subtilis Sco

Sco is a red copper protein that plays an essential yet poorly understood role in the metalation of the Cu(A) center of cytochrome oxidase, and is stable in both the Cu(I) and Cu(II) forms. To determine which oxidation state is important for function, we constructed His135 to Met or selenomethionine (SeM) variants that were designed to stabilize the Cu(I) over the Cu(II) state. H135M was unable to complement a scoΔ strain of Bacillus subtilis, indicating that the His to Met substitution abrogated cytochrome oxidase maturation. The Cu(I) binding affinities of H135M and H135SeM were comparable to that of the WT and 100-fold tighter than that of the H135A variant. The coordination chemistry of the H135M and H135SeM variants was studied by UV/vis, EPR, and XAS spectroscopy in both the Cu(I) and the Cu(II) forms. Both oxidation states bound copper via the S atoms of C45, C49 and M135. In particular, EXAFS data collected at both the Cu and the Se edges of the H135SeM derivative provided unambiguous evidence for selenomethionine coordination. Whereas the coordination chemistry and copper binding affinity of the Cu(I) state closely resembled that of the WT protein, the Cu(II) state was unstable, undergoing autoreduction to Cu(I). H135M also reacted faster with H(2)O(2) than WT Sco. These data, when coupled with the complete elimination of function in the H135M variant, imply that the Cu(I) state cannot be the sole determinant of function; the Cu(II) state must be involved in function at some stage of the reaction cycle.

Magnetic circular dichroism spectroscopy as a probe of the structures of the metal sites in metalloproteins

Magnetic circular dichroism (MCD) is a powerful probe of the electronic and geometric structures of metal centres in metalloproteins. MCD has provided significant insight into the nature of the axial donors at haem centres and, more recently, sophisticated methods for the analysis of MCD spectra have had a major impact on the study of the electronic structures of the ground states of a range of Cu, non-haem iron and Mo-containing active sites. This detail, together with data from other complimentary spectroscopies, has played a major role in defining the chemistry underpinning the catalysis achieved by these metal centres.

Density functional theory for transition metals and transition metal chemistry

We introduce density functional theory and review recent progress in its application to transition metal chemistry. Topics covered include local, meta, hybrid, hybrid meta, and range-separated functionals, band theory, software, validation tests, and applications to spin states, magnetic exchange coupling, spectra, structure, reactivity, and catalysis, including molecules, clusters, nanoparticles, surfaces, and solids.

Probing the role of the proximal heme ligand in cytochrome P450cam by recombinant incorporation of selenocysteine

The unique monooxygenase activity of cytochrome P450cam has been attributed to coordination of a cysteine thiolate to the heme cofactor. To investigate this interaction, we replaced cysteine with the more electron-donating selenocysteine. Good yields of the selenoenzyme were obtained by bacterial expression of an engineered gene containing the requisite UGA codon for selenocysteine and a simplified yet functional selenocysteine insertion sequence (SECIS). The sulfur-to-selenium substitution subtly modulates the structural, electronic, and catalytic properties of the enzyme. Catalytic activity decreases only 2-fold, whereas substrate oxidation becomes partially uncoupled from electron transfer, implying a more complex role for the axial ligand than generally assumed.

Copper-Carbon Bonds in Mechanistic and Structural Probing of Proteins as well as in Situations where Copper is a Catalytic or Receptor Site

While copper-carbon bonds are well appreciated in organometallic synthetic chemistry, such occurrences are less known in biological settings. By far, the greatest incidence of copper-carbon moieties is in bioinorganic research aimed at probing copper protein active site structure and mechanism; for example, carbon monoxide (CO) binding as a surrogate for O2. Using infrared (IR) spectroscopy, CO coordination to cuprous sites has proven to be an extremely useful tool for determining active site copper ligation (e.g., donor atom number and type). The coupled (hemocyanin, tyrosinase, catechol oxidase) and non-coupled (peptidylglycine α-hydroxylating monooxygenase, dopamine β-monooxygenase) binuclear copper proteins as well as the heme-copper oxidases (HCOs) have been studied extensively via this method. In addition, environmental changes within the vicinity of the active site have been determined based on shifts in the CO stretching frequencies, such as for copper amine oxidases, nitrite reductases and again in the binuclear proteins and HCOs. In many situations, spectroscopic monitoring has provided kinetic and thermodynamic data on CuI-CO formation and CO dissociation from copper(I); recently, processes occurring on a femtosecond timescale have been reported. Copper-cyano moieties have also been useful for obtaining insights into the active site structure and mechanisms of copper-zinc superoxide dismutase, azurin, nitrous oxide reductase, and multi-copper oxidases. Cyanide is a good ligand for both copper(I) and copper(II), therefore multiple physical-spectroscopic techniques can be applied. A more obvious occurrence of a “Cu-C” moiety was recently described for a CO dehydrogenase which contains a novel molybdenum-copper catalytic site. A bacterial copper chaperone (CusF) was recently established to have a novel d-π interaction comprised of copper(I) with the arene containing side-chain of a tryptophan amino acid residue. Meanwhile, good evidence exists that a plant receptor site (ETR1) utilizes copper(I) to sense ethylene, a growth hormone. A copper olfactory receptor has also been suggested. All of the above mentioned occurrences or uses of carbon-containing substrates and/or probes are reviewed and discussed within the framework of copper proteins and other relevant systems.

Expressed protein ligation for metalloprotein design and engineering

Metalloproteins contain highly specialized metal-binding sites that are designed to accept specific metal ions to maintain correct function. Although many of the sites have been modified with success, the relative paucity of functional group availability within proteinogenic amino acids can sometimes leave open questions about specific functions of the metal binding ligands. Attaining a more thorough analysis of individual amino acid function within metalloproteins has been realized using expressed protein ligation (EPL). Here we describe our recent efforts using EPL to incorporate nonproteinogenic cysteine and methionine analogues into the type 1 copper site found in Pseudomonas aeruginosa azurin.