Reduction potential tuning of the blue copper center in Pseudomonas aeruginosa azurin by the axial methionine as probed by unnatural amino acids

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.

Facile incorporation of non-canonical heme ligands in myoglobin through chemical protein synthesis

The incorporation of non-canonical amino acids (ncAAs) into the metal coordination environments of proteins has endowed metalloproteins with enhanced properties and novel activities, particularly in hemoproteins. In this work, we disclose a scalable synthetic strategy that enables the production of myoglobin (Mb) variants with non-canonical heme ligands, i.e., HoCys and f4Tyr. The ncAA-containing Mb* variants (with H64V/V68A mutations) were obtained through two consecutive native chemical ligations and a subsequent desulfurization step, with overall isolated yield up to 28.6 % in over 10-milligram scales. After refolding and heme b cofactor reconstitution, the synthetic Mb* variants showed typical electronic absorption bands. When subjected to the catalysis of the cyclopropanation of styrene, both synthetic variants, however, were not as competent as the His-ligated Mb*. We envisioned that the synthetic method reported herein would be useful for incorporating a variety of ncAAs with diverse structures and properties into Mb for varied purposes.

Azole reagents enabled ligation of peptide acyl pyrazoles for chemical protein synthesis

Native chemical ligation (NCL) has been playing an increasingly important role in chemical protein synthesis (CPS). More efficient ligation methods that circumvent the requirement of a peptidyl thioester and thiol additive-which allow the following desulfurization or refolding in one pot-are urgently needed for the synthesis of more complex protein targets and in large quantities. Herein, we discover that the weak acyl donor peptidyl N-acyl pyrazole can be activated by azole reagents like 3-methylpyrazole or imidazole to facilitate its ligation directly with an N-terminal cysteine peptide. As it requires no thioester or thiol additive, this ligation strategy can be conveniently combined with metal-free desulfurization (MFD) or oxidative protein folding to allow various one-pot protocols. The utility and generality of the strategy are showcased by the total synthesis of ubiquitin via an N-to-C sequential ligation-MFD strategy, the semi-synthesis of the copper protein azurin, and the efficient assembly of a sulfated hirudin variant and the cyclotide kalata B1, all in a one-pot fashion.

Chemical protein synthesis enabled engineering of saccharide oxidative cleavage activity in artificial metalloenzymes

Chemical protein (semi-)synthesis is a powerful technique allowing the incorporation of unnatural functionalities at any desired protein site. Herein we describe a facile one-pot semi-synthetic strategy for the construction of a type 2 copper center in the active site of azurin, which is achieved by substitution of Met121 with unnatural amino acid residues bearing a strong ligand N,N-bis(pyridylmethyl)amine (DPA) to mimic the function of typical histidine brace-bearing copper monooxygenases, such as lytic polysaccharide monooxygenases (LPMOs) involved in polysaccharide breakdown. The semi-synthetic proteins were routinely obtained in over 10-mg scales to allow for spectroscopic measurements (UV-Vis, CD, and EPR), which provides structural evidences for the CuII-DPA-modified azurins. 4-nitrophenyl-β-D-glucopyranoside (PNPG) was used as a model substrate for the H2O2-driven oxidative cleavage reaction facilitated by semi-synthetic azurins, and the CuII-6 complex showed a highest activity (TTN 253). Interestingly, our semi-synthetic azurins were able to tolerate high H2O2 concentrations (up to 4000-fold of the enzyme), making them promising for practical applications. Collectively, we establish that chemical protein synthesis can be exploited as a reliable technology in affording large quantities of artificial metalloproteins to facilitate the transformation of challenging chemical reactions.

Integration of molecular and computational approaches paints a holistic portrait of obscure metabolisms

Microorganisms are essential drivers of earth's geochemical cycles. However, the significance of elemental redox cycling mediated by microorganisms is often underestimated beyond the most well-studied nutrient cycles. Phosphite, (per)chlorate, and iodate are each considered esoteric substrates metabolized by microorganisms. However, recent investigations have indicated that these metabolisms are widespread and ubiquitous, affirming a need to continue studying the underlying microbiology to understand their biogeochemical effects and their interface with each other and our biosphere. This review focuses on combining canonical techniques of culturing microorganisms with modern omic approaches to further our understanding of obscure metabolic pathways and elucidate their importance in global biogeochemical cycles. Using these approaches, marker genes of interest have already been identified for phosphite, (per)chlorate, and iodate using traditional microbial physiology and genetics. Subsequently, their presence was queried to reveal the distribution of metabolic pathways in the environment using publicly available databases. In conjunction with each other, computational and experimental techniques provide a more comprehensive understanding of the location of these microorganisms, their underlying biochemistry and genetics, and how they tie into our planet's geochemical cycles.

A study of ab initio folding of chignolins using replica-exchange molecular dynamics simulations

More and more studies have confirmed the importance of polarization effects in hydrogen bonding interactions in protein folding simulations. In this paper, a recently developed charge update scheme termed polarized structure-specific backbone charge (PSBC) model was applied to the folding of 10-residue chignolin. A comparison between simulations performed using PSBC and a nonpolarizable (AMBER99SB) force field demonstrably showed the importance of the electrostatic polarization effect in the folding of the short β-hairpin peptide by a series of analyses such as DSSP, free-energy landscape, hydrogen bond occupancy, and melting curve. The PSBC model was further validated by folding two other chignolin variants.

The structure of plastocyanin tunes the midpoint potential by restricting axial ligation of the reduced copper ion

Blue copper proteins are models for illustrating how proteins tune metal properties. Nevertheless, the mechanisms by which the protein controls the metal site remain to be fully elucidated. A hindrance is that the closed shell Cu(I) site is inaccessible to most spectroscopic analyses. Carbon deuterium (C-D) bonds used as vibrational probes afford nonperturbative, selective characterization of the key cysteine and methionine copper ligands in both redox states. The structural integrity of Nostoc plastocyanin was perturbed by disrupting potential hydrogen bonds between loops of the cupredoxin fold via mutagenesis (S9A, N33A, N34A), variably raising the midpoint potential. The C-D vibrations show little change to suggest substantial alteration to the Cu(II) coordination in the oxidized state or in the Cu(I) interaction with the cysteine ligand. They rather indicate, along with visible and NMR spectroscopy, that the methionine ligand distinctly interacts more strongly with the Cu(I) ion, in line with the increases in midpoint potential. Here we show that the protein structure determines the redox properties by restricting the interaction between the methionine ligand and Cu(I) in the reduced state.

Electron Tunneling in Biology: When Does it Matter?

Electrons can tunnel between cofactor molecules positioned along biological electron transport chains up to a distance of ≃ 20 Å on the millisecond time scale of enzymatic turnover. This tunneling range determines the design of biological energy chains facilitating the cross-membrane transport of electrons. Tunneling distance and cofactors' redox potentials become the main physical parameters affecting the rate of electron transport. In addition, universal charge-transport properties are assigned to all proteins, making protein identity, flexibility, and dynamics insignificant. This paradigm is challenged by dynamical models of electron transfer, showing that the electron hopping rate is constant within the crossover distance R* ≃ 12 Å, followed with an exponential falloff at longer distances. If this hypothesis is fully confirmed, natural and man-made energy chains for electron transport should be best designed by placing redox cofactors near the crossover distance R*. Protein flexibility and dynamics affect the magnitude of the maximum hopping rate within the crossover distance. Changes in protein flexibility between forward and backward transitions contribute to vectorial charge transport. For biological energy chains, charge transport through proteins is not defined by universal parameters, and protein identity matters.

Tuning the Type 1 Reduction Potential of Multicopper Oxidases: Uncoupling the Effects of Electrostatics and H-Bonding to Histidine Ligands

In multicopper oxidases (MCOs), the type 1 (T1) Cu accepts electrons from the substrate and transfers these to the trinuclear Cu cluster (TNC) where O₂ is reduced to H₂O. The T1 potential in MCOs varies from 340 to 780 mV, a range not explained by the existing literature. This study focused on the ∼350 mV difference in potential of the T1 center in Fet3p and Trametes versicolor laccase (TvL) that have the same 2His1Cys ligand set. A range of spectroscopies performed on the oxidized and reduced T1 sites in these MCOs shows that they have equivalent geometric and electronic structures. However, the two His ligands of the T1 Cu in Fet3p are H-bonded to carboxylate residues, while in TvL they are H-bonded to noncharged groups. Electron spin echo envelope modulation spectroscopy shows that there are significant differences in the second-sphere H-bonding interactions in the two T1 centers. Redox titrations on type 2-depleted derivatives of Fet3p and its D409A and E185A variants reveal that the two carboxylates (D409 and E185) lower the T1 potential by 110 and 255-285 mV, respectively. Density functional theory calculations uncouple the effects of the charge of the carboxylates and their difference in H-bonding interactions with the His ligands on the T1 potential, indicating 90-150 mV for anionic charge and ∼100 mV for a strong H-bond. Finally, this study provides an explanation for the generally low potentials of metallooxidases relative to the wide range of potentials of the organic oxidases in terms of different oxidized states of their TNCs involved in catalytic turnover.

Effect of Water Deuteration on Protein Electron Transfer

Traditional theories of long-range protein electron transfer describe the reaction rate in terms of the tunneling distance and the reaction free energy. They do not recognize two physical effects: (i) local wetting of the active site by hydration water and (ii) protein identity affecting the rate through dynamics and flexibility. We find, by molecular dynamics simulations, a significant, ∼25 times, slowing down of the rate of protein electron transfer upon deuteration. H/D substitution changes the rate constant pre-exponential factor in the regime of electron transfer controlled by medium dynamics. Switching from light to heavy water increases the effective medium relaxation time. The effect is caused by both a global change in the flexibility of the protein backbone and locally stronger hydrogen bonds to charged residues.

Theory of Protein Charge Transfer: Electron Transfer between Tryptophan Residue and Active Site of Azurin

One reaction step in the conductivity relay of azurin, electron transfer between the Cu-based active site and the tryptophan residue, is studied theoretically and by classical molecular dynamics simulations. Oxidation of tryptophan results in electrowetting of this residue. This structural change makes the free energy surfaces of electron transfer nonparabolic as described by the Q-model of electron transfer. We analyze the medium dynamical effect on protein electron transfer produced by coupled Stokes-shift dynamics and the dynamics of the donor-acceptor distance modulating electron tunneling. The equilibrium donor-acceptor distance falls in the plateau region of the rate constant, where it is determined by the protein-water dynamics, and the probability of electron tunneling does not affect the rate. The crossover distance found here puts most intraprotein electron-transfer reactions under the umbrella of dynamical control. The crossover between the medium-controlled and tunneling-controlled kinetics is combined with the effect of the protein-water medium on the activation barrier to formulate principles of tunability of protein-based charge-transfer chains. The main principle in optimizing the activation barrier is the departure from the Gaussian-Gibbsian statistics of fluctuations promoting activated transitions. This is achieved either by incomplete (nonergodic) sampling, breaking the link between the Stokes-shift and variance reorganization energies, or through wetting-induced structural changes of the enzyme's active site.

Designing Artificial Metalloenzymes by Tuning of the Environment beyond the Primary Coordination Sphere

Metalloenzymes catalyze a variety of reactions using a limited number of natural amino acids and metallocofactors. Therefore, the environment beyond the primary coordination sphere must play an important role in both conferring and tuning their phenomenal catalytic properties, enabling active sites with otherwise similar primary coordination environments to perform a diverse array of biological functions. However, since the interactions beyond the primary coordination sphere are numerous and weak, it has been difficult to pinpoint structural features responsible for the tuning of activities of native enzymes. Designing artificial metalloenzymes (ArMs) offers an excellent basis to elucidate the roles of these interactions and to further develop practical biological catalysts. In this review, we highlight how the secondary coordination spheres of ArMs influence metal binding and catalysis, with particular focus on the use of native protein scaffolds as templates for the design of ArMs by either rational design aided by computational modeling, directed evolution, or a combination of both approaches. In describing successes in designing heme, nonheme Fe, and Cu metalloenzymes, heteronuclear metalloenzymes containing heme, and those ArMs containing other metal centers (including those with non-native metal ions and metallocofactors), we have summarized insights gained on how careful controls of the interactions in the secondary coordination sphere, including hydrophobic and hydrogen bonding interactions, allow the generation and tuning of these respective systems to approach, rival, and, in a few cases, exceed those of native enzymes. We have also provided an outlook on the remaining challenges in the field and future directions that will allow for a deeper understanding of the secondary coordination sphere a deeper understanding of the secondary coordintion sphere to be gained, and in turn to guide the design of a broader and more efficient variety of ArMs.

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

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

Altered coordination in a blue copper protein upon association with redox partner revealed by carbon-deuterium vibrational probes

Proteins tune the reactivity of metal sites; less understood is the impact of association with a redox partner. We demonstrate the utility of carbon-deuterium labels for selective analysis of delicate...

Delivery of Various Cargos into Cancer Cells and Tissues via Cell-Penetrating Peptides: A Review of the Last Decade

Cell-penetrating peptides (CPPs), also known as protein transduction domains, are a class of diverse amino acid sequences with the ability to cross cellular membranes. CPPs can deliver several bioactive cargos, including proteins, peptides, nucleic acids and chemotherapeutics, into cells. Ever since their discovery, synthetic and natural CPPs have been utilized in therapeutics delivery, gene editing and cell imaging in fundamental research and clinical experiments. Over the years, CPPs have gained significant attention due to their low cytotoxicity and high transduction efficacy. In the last decade, multiple investigations demonstrated the potential of CPPs as carriers for the delivery of therapeutics to treat various types of cancer. Besides their remarkable efficacy owing to fast and efficient delivery, a crucial benefit of CPP-based cancer treatments is delivering anticancer agents selectively, rather than mediating toxicities toward normal tissues. To obtain a higher therapeutic index and to improve cell and tissue selectivity, CPP-cargo constructions can also be complexed with other agents such as nanocarriers and liposomes to obtain encouraging outcomes. This review summarizes various types of CPPs conjugated to anticancer cargos. Furthermore, we present a brief history of CPP utilization as delivery systems for anticancer agents in the last decade and evaluate several reports on the applications of CPPs in basic research and preclinical studies.

Contributions of methionine to recognition of trimethyllysine in aromatic cage of PHD domains: implications of polarizability, hydrophobicity, and charge on binding

Recognition of trimethyllysine (Kme3) by reader proteins is an important regulator of gene expression. This recognition event is mediated by an aromatic cage made up of 2-4 aromatic residues in the reader proteins that bind Kme3 via cation-π interactions. A small subset of reader proteins contain a methionine (Met) residue in place of an aromatic sidechain in the binding pocket. The unique role of sulfur in molecular recognition has been demonstrated in a number of noncovalent interactions recently, including interactions of thiols, thioethers, and sulfoxides with aromatic rings. However, the interaction of a thioether with an ammonium ion has not previously been investigated and the role of Met in binding Kme3 has not yet been explored. Herein, we systematically vary the Met in two reader proteins, DIDO1 and TAF3, and the ligand, Kme3 or its neutral analog tert-butyl norleucine (tBuNle), to determine the role of Met in the recognition of the cationic Kme3. Our studies demonstrate that Met contributes to binding via dispersion forces, with about an equal contribution to binding Kme3 and tBuNle, indicating that electrostatic interactions do not play a role. During the course of these studies, we also discovered that DIDO1 exhibits equivalent binding to tBuNle and Kme3 through a change in the mechanism of binding.

Noncanonical Heme Ligands Steer Carbene Transfer Reactivity in an Artificial Metalloenzyme*

Changing the primary metal coordination sphere is a powerful strategy for tuning metalloprotein properties. Here we used amber stop codon suppression with engineered pyrrolysyl-tRNA synthetases, including two newly evolved enzymes, to replace the proximal histidine in myoglobin with N_δ -methylhistidine, 5-thiazoylalanine, 4-thiazoylalanine and 3-(3-thienyl)alanine. In addition to tuning the heme redox potential over a >200 mV range, these noncanonical ligands modulate the protein's carbene transfer activity with ethyl diazoacetate. Variants with increased reduction potential proved superior for cyclopropanation and N-H insertion, whereas variants with reduced E^o values gave higher S-H insertion activity. Given the functional importance of histidine in many enzymes, these genetically encoded analogues could be valuable tools for probing mechanism and enabling new chemistries.

Chemical modifications of proteins and their applications in metalloenzyme studies

Protein chemical modifications are important tools for elucidating chemical and biological functions of proteins. Several strategies have been developed to implement these modifications, including enzymatic tailoring reactions, unnatural amino acid incorporation using the expanded genetic codes, and recognition-driven transformations. These technologies have been applied in metalloenzyme studies, specifically in dissecting their mechanisms, improving their enzymatic activities, and creating artificial enzymes with non-natural activities. Herein, we summarize some of the recent efforts in these areas with an emphasis on a few metalloenzyme case studies.

Enantiopure 5-CF3-Proline: Synthesis, Incorporation in Peptides, and Tuning of the Peptide Bond Geometry

The straightforward synthesis of enantiopure 5-(R)-and 5-(S)-trifluoromethylproline is reported. The key steps are a Ruppert-Prakash reagent addition on l-pyroglutamic esters followed by an elimination reaction and a selective reduction. The solution-phase and solid-phase incorporation of this unprotected enantiopure fluorinated amino acid in a short peptide chain was demonstrated. Compared to proline, the CF₃ group provides a decrease of the trans to cis amide bond isomerization energy and an increase of the cis conformer population.

A Binuclear CuA Center Designed in an All α-Helical Protein Scaffold

The primary and secondary coordination spheres of metal binding sites in metalloproteins have been investigated extensively, leading to the creation of high-performing functional metalloproteins; however, the impact of the overall structure of the protein scaffold on the unique properties of metalloproteins has rarely been studied. A primary example is the binuclear Cu_A center, an electron transfer cupredoxin domain of photosynthetic and respiratory complexes and, recently, a protein coregulated with particulate methane and ammonia monooxygenases. The redox potential, Cu-Cu spectroscopic features, and a valence delocalized state of Cu_A are difficult to reproduce in synthetic models, and every artificial protein Cu_A center to-date has used a modified cupredoxin. Here, we present a fully functional Cu_A center designed in a structurally nonhomologous protein, cytochrome c peroxidase (CcP), by only two mutations (Cu_ACcP). We demonstrate with UV-visible absorption, resonance Raman, and magnetic circular dichroism spectroscopy that Cu_ACcP is valence delocalized. Continuous wave and pulsed (HYSCORE) X-band EPR show it has a highly compact g_z area and small A_z hyperfine principal value with g and A tensors that resemble axially perturbed Cu_A. Stopped-flow kinetics found that Cu_A formation proceeds through a single T2Cu intermediate. The reduction potential of Cu_ACcP is comparable to native Cu_A and can transfer electrons to a physiological redox partner. We built a structural model of the designed Cu binding site from extended X-ray absorption fine structure spectroscopy and validated it by mutation of coordinating Cys and His residues, revealing that a triad of residues (R48C, W51C, and His52) rigidly arranged on one α-helix is responsible for chelating the first Cu(II) and that His175 stabilizes the binuclear complex by rearrangement of the CcP heme-coordinating helix. This design is a demonstration that a highly conserved protein fold is not uniquely necessary to induce certain characteristic physical and chemical properties in a metal redox center.

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.

Enantiopure α-Trifluoromethylated Aziridine-2-carboxylic Acid (α-TfmAzy): Synthesis and Peptide Coupling

A straightforward synthesis of enantiopure α-trifluoromethyl aziridine-2-carboxylic acid (α-TfmAzy) is reported from a trifluoropyruvate derived enantiopure oxazolidine. A key Strecker-type synthetic step and a late cyanide basic hydrolysis gave the target compounds in six steps and 41% yield. A final peptide coupling was performed to demonstrate the usefulness of this highly constrained fluorinated unnatural amino acid.

Tuning Structure and Dynamics of Blue Copper Azurin Junctions via Single Amino-Acid Mutations

In the growing field of biomolecular electronics, blue-copper Azurin stands out as one of the most widely studied protein in single-molecule contacts. Interestingly, despite the paramount importance of the structure/dynamics of molecular contacts in their transport properties, these factors remain largely unexplored from the theoretical point of view in the context of single Azurin junctions. Here we address this issue using all-atom Molecular Dynamics (MD) of Pseudomonas Aeruginosa Azurin adsorbed to a Au(111) substrate. In particular, we focus on the structure and dynamics of the free/adsorbed protein and how these properties are altered upon single-point mutations. The results revealed that wild-type Azurin adsorbs on Au(111) along two well defined configurations: one tethered via cysteine groups and the other via the hydrophobic pocket surrounding the Cu 2 + . Surprisingly, our simulations revealed that single amino-acid mutations gave rise to a quenching of protein vibrations ultimately resulting in its overall stiffening. Given the role of amino-acid vibrations and reorientation in the dehydration process at the protein-water-substrate interface, we suggest that this might have an effect on the adsorption process of the mutant, giving rise to new adsorption configurations.

Correlation between the T1 copper reduction potential and catalytic activity of a small laccase

Laccases are multicopper enzymes that catalyze oxidation of electron-rich substrates coupled to reduction of molecular oxygen to water. Since the Type 1 copper (T1 Cu) is the site where electrons are withdrawn from the substrate, it is assumed that the reduction potential of this copper correlates with enzyme activity. Herein, we studied the correlation of the T1 Cu reduction potential and the enzymatic activity of the small two-domain laccase Ssl1 from Streptomyces sviceus. For a systematic approach, we aimed to minimize any effects other than the reduction potential difference. To this end, we constructed a series of Ssl1 mutants with reduction potentials varying from <290 to 560 mV. Along with the hydrophobicity of the axial ligand of the T1 Cu also structural changes in the substrate binding site and additional hydrogen bonding increased the reduction potential. Enzyme activity experiments demonstrated that the T1 Cu reduction potential has a different effect on oxidation of different substrates. Whereas there was no obvious correlation between the T1 Cu reduction potential and kinetic parameters for the oxidation of syringaldazine (with a reduction potential of 390 mV), a good correlation was observed between the T1 Cu reduction potential and the conversion of substituted phenols with reduction potentials between 660 and 820 mV. This correlation was pronounced for the Ssl1 variants with reduction potentials above 470 mV, which demonstrated increased activities also during the oxidation of two dyes, alizarin red S and indigo carmine.

Redox Tuning via Ligand-Induced Geometric Distortions at a YMn3O4 Cubane Model of the Biological Oxygen Evolving Complex

The function of proteins involved in electron transfer is dependent on cofactors attaining the necessary reduction potentials. We establish a mode of cluster redox tuning through steric pressure on a synthetic model related to Photosystem II. Resembling the cuboidal [CaMn₃O₄] subsite of the biological oxygen evolving complex (OEC), [Mn₄O₄] and [YMn₃O₄] complexes featuring ligands of different basicity and chelating properties were characterized by cyclic voltammetry. In the absence of ligand-induced distortions, increasing the basicity of the ligands results in a decrease of cluster reduction potential. Contraction of Y-oxo/Y–Mn distances by 0.1/0.15 Å enforced by a chelating ligand results in an increase of cluster reduction potential even in the presence of strongly basic donors. Related protein-induced changes in Ca-oxo/Ca–Mn distances may have similar effects in tuning the redox potential of the OEC through entatic states and may explain the cation size dependence on the progression of the S-state cycle.

The redox properties of [YMn₃O₄] and [Mn₄O₄] complexes featuring bridging ligands of different basicity and chelating properties are reported. In the absence of ligand-induced geometric distortions, increasing the basicity of the ligands results in a decrease of cluster reduction potential. A chelating ligand results in contractions of Y-oxo distances by ∼0.1 Å, which correlates with an increase of cluster reduction potential even in the presence of strongly basic donors.

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.

What Are We Missing by Not Measuring the Net Charge of Proteins?

The net electrostatic charge (Z) of a folded protein in solution represents a bird's eye view of its surface potentials-including contributions from tightly bound metal, solvent, buffer, and cosolvent ions-and remains one of its most enigmatic properties. Few tools are available to the average biochemist to rapidly and accurately measure Z at pH≠pI. Tools that have been developed more recently seem to go unnoticed. Most scientists are content with this void and estimate the net charge of a protein from its amino acid sequence, using textbook values of pK_a . Thus, Z remains unmeasured for nearly all folded proteins at pH≠pI. When marveling at all that has been learned from accurately measuring the other fundamental property of a protein-its mass-one wonders: what are we missing by not measuring the net charge of folded, solvated proteins? A few big questions immediately emerge in bioinorganic chemistry. When a single electron is transferred to a metalloprotein, does the net charge of the protein change by approximately one elementary unit of charge or does charge regulation dominate, that is, do the pK_a values of most ionizable residues (or just a few residues) adjust in response to (or in concert with) electron transfer? Would the free energy of charge regulation (ΔΔG_z ) account for most of the outer sphere reorganization energy associated with electron transfer? Or would ΔΔG_z contribute more to the redox potential? And what about metal binding itself? When an apo-metalloprotein, bearing minimal net negative charge (e.g., Z=-2.0) binds one or more metal cations, is the net charge abolished or inverted to positive? Or do metalloproteins regulate net charge when coordinating metal ions? The author's group has recently dusted off a relatively obscure tool-the "protein charge ladder"-and used it to begin to answer these basic questions.

Accurate Quantum Mechanical/Molecular Mechanical Calculations of Reduction Potentials in Azurin Variants

Understanding the regulation mechanism and molecular determinants of the reduction potential of metalloprotein is a major challenge. An ab initio quantum mechanical/molecular mechanical (QM/MM) method combining the minimum free energy path (MFEP) and fractional number of electron (FNE) approaches has been developed in our group to simulate the redox processes of large systems. The FNE scheme provides an efficient unique description for the redox process, while the MFEP method provides improved conformational sampling on complex environments such as protein in the QM/MM calculations. The reduction potentials of wild-type and seven mutants of azurin, a type 1 copper metalloprotein, were simulated with the QM/MM-MFEP+FNE approach in this paper. A range of 350 mV for the variations of the reduction potentials of these azurin proteins was reproduced faithfully with relative errors around 20 mV. The correlation between structural interactions and reduction potentials observed in simulations provides in-depth insight into the regulation of reduction potentials, which potentially can also be very useful to the engineering of metalloprotein-based electrocatalysts in artificial photosynthesis. The excellent accuracy and efficiency of the QM/MM-MFEP+FNE approach demonstrate the potential for simulations of many electron transfer processes in condensed phases and biochemical systems.

Engineered Metalloenzymes with Non-Canonical Coordination Environments

Nature employs a limited number of genetically encoded, metal-coordinating residues to create metalloenzymes with diverse structures and functions. Engineered components of the cellular translation machinery can now be exploited to encode non-canonical ligands with user-defined electronic and structural properties. This ability to install "chemically programmed" ligands into proteins can provide powerful chemical probes of metalloenzyme mechanism and presents excellent opportunities to create metalloprotein catalysts with augmented properties and novel activities. In this Concept article, we provide an overview of several recent studies describing the creation of engineered metalloenzymes with interesting catalytic properties, and reveal how characterization of these systems has advanced our understanding of nature's bioinorganic mechanisms. We also highlight how powerful laboratory evolution protocols can be readily adapted to allow optimization of metalloenzymes with non-canonical ligands. This approach combines beneficial features of small molecule and protein catalysis by allowing the installation of a greater variety of local metal coordination environments into evolvable protein scaffolds, and holds great promise for the future creation of powerful metalloprotein catalysts for a host of synthetically valuable transformations.

Synthesis and physical chemical properties of 2-amino-4-(trifluoromethoxy)butanoic acid - a CF3O-containing analogue of natural lipophilic amino acids

2-Amino-2-(trifluoromethoxy)butanoic acid (O-trifluoromethyl homoserine) was synthesized as a racemate and in both enantiomeric forms. The measured pK_a and log D values establish the compound as a promising analogue of natural aliphatic amino acids.

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.

Prediction of Reduction Potentials of Copper Proteins with Continuum Electrostatics and Density Functional Theory

Blue copper proteins, such as azurin, show dramatic changes in Cu²⁺/Cu⁺ reduction potential upon mutation over the full physiological range. Hence, they have important functions in electron transfer and oxidation chemistry and have applications in industrial biotechnology. The details of what determines these reduction potential changes upon mutation are still unclear. Moreover, it has been difficult to model and predict the reduction potential of azurin mutants and currently no unique procedure or workflow pattern exists. Furthermore, high‐level computational methods can be accurate but are too time consuming for practical use. In this work, a novel approach for calculating reduction potentials of azurin mutants is shown, based on a combination of continuum electrostatics, density functional theory and empirical hydrophobicity factors. Our method accurately reproduces experimental reduction potential changes of 30 mutants with respect to wildtype within experimental error and highlights the factors contributing to the reduction potential change. Finally, reduction potentials are predicted for a series of 124 new mutants that have not yet been investigated experimentally. Several mutants are identified that are located well over 10 Å from the copper center that change the reduction potential by more than 85 mV. The work shows that secondary coordination sphere mutations mostly lead to long‐range electrostatic changes and hence can be modeled accurately with continuum electrostatics.

Stabilization of protein structure through π-π interaction in the second coordination sphere of pseudoazurin

Noncovalent, weak interactions in the second coordination sphere of the copper active site of Pseudoazurin (PAz) from Achromobacter cycloclastes were examined using a series of Met16X variants. In this study, the differences in protein stability due to the changes in the nature of the 16th amino acid (Met, Phe, Val, Ile) were investigated by electrospray ionization mass spectrometry (ESI-MS) and far-UV circular dichroism (CD) as a result of acid denaturation. The percentage of native states (folded holo forms) of Met16Phe variants was estimated to be 75% at pH 2.9 although the wild-type (WT), Met16Val and Met16Ile PAz, became completely unfolded. The high stability under acidic conditions is correlated with the result of the active site being stabilized by the aromatic substitution of the Met16 residue. The π-π interaction in the second coordination sphere makes a significant contribution to the stability of active site and the protein matrix.

Tuning of Enthalpic/Entropic Parameters of a Protein Redox Center through Manipulation of the Electronic Partition Function

Manipulation of the partition function (Q) of the redox center Cu_A from cytochrome c oxidase is attained by tuning the accessibility of a low lying alternative electronic ground state and by perturbation of the electrostatic potential through point mutations, loop engineering and pH variation. We report clear correlations of the entropic and enthalpic contributions to redox potentials with Q and with the identity and hydrophobicity of the weak axial ligand, respectively.

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.

Improving artificial metalloenzymes' activity by optimizing electron transfer

This feature article discusses the strategies to optimize electron transfer efficiency, towards enhancing the activity of artificial metalloenzymes.

Probing the role of the backbone carbonyl interaction with the CuA center in azurin by replacing the peptide bond with an ester linkage

The role of a backbone carbonyl interaction with an engineered CuA center in azurin was investigated by developing a method of synthesis and incorporation of a depsipeptide where one of the amide bonds in azurin is replaced by an ester bond using expressed protein ligation. Studies by electronic absorption and electron paramagnetic resonance spectroscopic techniques indicate that, while the substitution does not significantly alter the geometry of the site, it weakens the axial interaction to the CuA center and strengthens the Cu-Cu bond, as evidenced by the blue shift of the near-IR absorption that has been assigned to the Cu-Cu ψ → ψ* transition. Interestingly, the changes in the electronic structure from the replacement did not result in a change in the reduction potential of the CuA center, suggesting that the diamond core structure of Cu2SCys2 is resistant to variations in axial interactions.