Designed azurins show lower reorganization free energies for intraprotein electron transfer

HYSCORE and QM/MM Studies of Second Sphere Variants of the Type 1 Copper Site in Azurin: Influence of Mutations on the Hyperfine Couplings of Remote Nitrogens

Secondary coordination sphere (SCS) interactions have been shown to play important roles in tuning reduction potentials and electron transfer (ET) properties of the Type 1 copper proteins, but the precise roles of these interactions are not fully understood. In this work, we examined the influence of F114P, F114N, and N47S mutations in the SCS on the electronic structure of the T1 copper center in azurin (Az) by studying the hyperfine couplings of (i) histidine remote Nε nitrogens and (ii) the amide Np using the two-dimensional (2D) pulsed electron paramagnetic resonance (EPR) technique HYSCORE (hyperfine sublevel correlation) combined with quantum mechanics/molecular mechanics (QM/MM) and DLPNO-CCSD calculations. Our data show that some components of hyperfine tensor and isotropic coupling in N47SAz and F114PAz (but not F114NAz) deviate by up to ∼±20% from WTAz, indicating that these mutations significantly influence the spin density distribution between the CuII site and coordinating ligands. Furthermore, our calculations support the assignment of Np to the backbone amide of residue 47 (both in Asn and Ser variants). Since the spin density distributions play an important role in tuning the covalency of the Cu-Scys bond of Type 1 copper center that has been shown to be crucial in controlling the reduction potentials, this study provides additional insights into the electron spin factor in tuning the reduction potentials and ET properties.

Contribution Factors of the First Kind Calculated for the Marcus Electron-Transfer Rate and Their Applications

In this study, we applied the concept of the "contribution factor of the first kind (CFFK)" to the original electron-transfer (ET) rate theory proposed by Marcus. Mathematical derivations provided simple and convenient formulas for estimating the relative contributions of ten physical and chemical parameters involved in the Marcus ET rate formula: (1) the maximum strength of the electronic coupling energy between two molecules, (2) the exponential decay rate of the electronic coupling energy versus the distance between both molecules, (3) the distance between both molecules, (4) the equilibrium distance between both molecules, (5) the Gibbs free energy, (6) reorganization free energy in the prefactor of the Marcus ET rate equation, (7) reorganization free energy in the denominator of the exponential term, (8) reorganization free energy in the argument of the exponential term, (9) Boltzmann constant times absolute temperature in the prefactor of the rate equation, and (10) Boltzmann constant times absolute temperature in the denominator of the exponential term. We applied our theories to (i) ET reactions at bacterial photosynthesis reaction centers, PSI and PSII, and soluble ferredoxins (Fd); (ii) intraprotein ET reactions for designed azurin mutants; and (iii) ET reactions in flavodoxin (Fld). The formulas and calculations suggest that the theory behind the CFFK is useful for quantitatively identifying major and minor physical and chemical factors and corresponding trade-offs, all of which affect the magnitude of the Marcus ET rate.

Pyruvate dehydrogenase operates as an intramolecular nitroxyl generator during macrophage metabolic reprogramming

M1 macrophages enter a glycolytic state when endogenous nitric oxide (NO) reprograms mitochondrial metabolism by limiting aconitase 2 and pyruvate dehydrogenase (PDH) activity. Here, we provide evidence that NO targets the PDH complex by using lipoate to generate nitroxyl (HNO). PDH E2-associated lipoate is modified in NO-rich macrophages while the PDH E3 enzyme, also known as dihydrolipoamide dehydrogenase (DLD), is irreversibly inhibited. Mechanistically, we show that lipoate facilitates NO-mediated production of HNO, which interacts with thiols forming irreversible modifications including sulfinamide. In addition, we reveal a macrophage signature of proteins with reduction-resistant modifications, including in DLD, and identify potential HNO targets. Consistently, DLD enzyme is modified in an HNO-dependent manner at Cys477 and Cys484, and molecular modeling and mutagenesis show these modifications impair the formation of DLD homodimers. In conclusion, our work demonstrates that HNO is produced physiologically. Moreover, the production of HNO is dependent on the lipoate-rich PDH complex facilitating irreversible modifications that are critical to NO-dependent metabolic rewiring.

Why charging Li-air batteries with current low-voltage mediators is slow and singlet oxygen does not explain degradation

Although Li-air rechargeable batteries offer higher energy densities than lithium-ion batteries, the insulating Li₂O₂ formed during discharge hinders rapid, efficient re-charging. Redox mediators are used to facilitate Li₂O₂ oxidation; however, fast kinetics at a low charging voltage are necessary for practical applications and are yet to be achieved. We investigate the mechanism of Li₂O₂ oxidation by redox mediators. The rate-limiting step is the outer-sphere one-electron oxidation of Li₂O₂ to LiO₂, which follows Marcus theory. The second step is dominated by LiO₂ disproportionation, forming mostly triplet-state O₂. The yield of singlet-state O₂ depends on the redox potential of the mediator in a way that does not correlate with electrolyte degradation, in contrast to earlier views. Our mechanistic understanding explains why current low-voltage mediators (<+3.3 V) fail to deliver high rates (the maximum rate is at +3.74 V) and suggests important mediator design strategies to deliver sufficiently high rates for fast charging at potentials closer to the thermodynamic potential of Li₂O₂ oxidation (+2.96 V).

Applicability of perturbed matrix method for charge transfer studies at bio/metallic interfaces: a case of azurin

As the field of nanoelectronics based on biomolecules such as peptides and proteins rapidly grows, there is a need for robust computational methods able to reliably predict charge transfer properties at bio/metallic interfaces. Traditionally, hybrid quantum-mechanical/molecular-mechanical techniques are employed for systems where the electron hopping transfer mechanism is applicable to determine physical parameters controlling the thermodynamics and kinetics of charge transfer processes. However, these approaches are limited by a relatively high computational cost when extensive sampling of a configurational space is required, like in the case of soft biomatter. For these applications, semi-empirical approaches such as the perturbed matrix method (PMM) have been developed and successfully used to study charge-transfer processes in biomolecules. Here, we explore the performance of PMM on prototypical redox-active protein azurin in various environments, from solution to vacuum interfaces with gold surfaces and protein junction. We systematically benchmarked the robustness and convergence of the method with respect to the quantum-centre size, size of the Hamiltonian, number of samples, and level of theory. We show that PMM can adequately capture all the trends associated with the structural and electronic changes related to azurin oxidation at bio/metallic interfaces.

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.

Reorganization Free Energy of Copper Proteins in Solution, in Vacuum, and on Metal Surfaces

Metalloproteins, known to efficiently transfer electronic charge in biological systems, recently found their utilization in nanobiotechnological devices where the protein is placed into direct contact with metal surfaces. The feasibility of oxidation/reduction of the protein redox sites is affected by the reorganization free energies, one of the key parameters determining the transfer rates. While their values have been measured and computed for proteins in their native environments, i.e., in aqueous solution, the reorganization free energies of dry proteins or proteins adsorbed to metal surfaces remain unknown. Here, we investigate the redox properties of blue copper protein azurin, a prototypical redox-active metalloprotein previously probed by various experimental techniques both in solution and on metal/vacuum interfaces. We used a hybrid QM/MM computational technique based on DFT to explore protein dynamics, flexibility, and corresponding reorganization free energies in aqueous solution, vacuum, and on vacuum gold interfaces. Somewhat surprisingly, the reorganization free energy only slightly decreases when azurin is dried because the loss of the hydration shell leads to larger flexibility of the protein near its redox site. At the vacuum gold surfaces, the energetics of the structure relaxation depends on the adsorption geometry, however, significant reduction of the reorganization free energy was not observed. These findings have important consequences for the charge transport mechanism in vacuum devices, showing that the free energy barriers for protein oxidation remain significant even under ultra-high vacuum conditions.

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.

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.

Mechanical Deformation and Electronic Structure of a Blue Copper Azurin in a Solid-State Junction

Protein-based electronics is an emerging field which has attracted considerable attention over the past decade. Here, we present a theoretical study of the formation and electronic structure of a metal-protein-metal junction based on the blue-copper azurin from pseudomonas aeruginosa. We focus on the case in which the protein is adsorbed on a gold surface and is contacted, at the opposite side, to an STM (Scanning Tunneling Microscopy) tip by spontaneous attachment. This has been simulated through a combination of molecular dynamics and density functional theory. We find that the attachment to the tip induces structural changes in the protein which, however, do not affect the overall electronic properties of the protein. Indeed, only changes in certain residues are observed, whereas the electronic structure of the Cu-centered complex remains unaltered, as does the total density of states of the whole protein.

Hole Hopping Across a Protein-Protein Interface

We have investigated photoinduced hole hopping in a Pseudomonas aeruginosa azurin mutant Re126WWCuI, where two adjacent tryptophan residues (W124 and W122) are inserted between the CuI center and a Re photosensitizer coordinated to a H126 imidazole (Re = ReI(H126)(CO)3(dmp)+, dmp = 4,7-dimethyl-1,10-phenanthroline). Optical excitation of this mutant in aqueous media (≤40 μM) triggers 70 ns electron transport over 23 Å, yielding a long-lived (120 μs) ReI(H126)(CO)3(dmp•-)WWCuII product. The Re126FWCuI mutant (F124, W122) is not redox-active under these conditions. Upon increasing the concentration to 0.2-2 mM, {Re126WWCuI}2 and {Re126FWCuI}2 are formed with the dmp ligand of the Re photooxidant of one molecule in close contact (3.8 Å) with the W122' indole on the neighboring chain. In addition, {Re126WWCuI}2 contains an interfacial tryptophan quadruplex of four indoles (3.3-3.7 Å apart). In both mutants, dimerization opens an intermolecular W122' → //*Re ET channel (// denotes the protein interface, *Re is the optically excited sensitizer). Excited-state relaxation and ET occur together in two steps (time constants of ∼600 ps and ∼8 ns) that lead to a charge-separated state containing a Re(H126)(CO)3(dmp•-)//(W122•+)' unit; then (CuI)' is oxidized intramolecularly (60-90 ns) by (W122•+)', forming ReI(H126)(CO)3(dmp•-)WWCuI//(CuII)'. The photocycle is closed by ∼1.6 μs ReI(H126)(CO)3(dmp•-) → //(CuII)' back ET that occurs over 12 Å, in contrast to the 23 Å, 120 μs step in Re126WWCuI. Importantly, dimerization makes Re126FWCuI photoreactive and, as in the case of {Re126WWCuI}2, channels the photoproduced "hole" to the molecule that was not initially photoexcited, thereby shortening the lifetime of ReI(H126)(CO)3(dmp•-)//CuII. Although two adjacent W124 and W122 indoles dramatically enhance CuI → *Re intramolecular multistep ET, the tryptophan quadruplex in {Re126WWCuI}2 does not accelerate intermolecular electron transport; instead, it acts as a hole storage and crossover unit between inter- and intramolecular ET pathways. Irradiation of {Re126WWCuII}2 or {Re126FWCuII}2 also triggers intermolecular W122' → //*Re ET, and the Re(H126)(CO)3(dmp•-)//(W122•+)' charge-separated state decays to the ground state by ∼50 ns ReI(H126)(CO)3(dmp•-)+ → //(W122•+)' intermolecular charge recombination. Our findings shed light on the factors that control interfacial hole/electron hopping in protein complexes and on the role of aromatic amino acids in accelerating long-range electron transport.

The Periodic Table's Impact on Bioinorganic Chemistry and Biology's Selective Use of Metal Ions

Despite the availability of a vast variety of metal ions in the periodic table, biology uses only a selective few metal ions. Most of the redox active metals used belong to the first row of transition metals in the periodic table and include Fe, Co, Ni, Mn and Cu. On the other hand, Ca, Zn and Mg are the most commonly used redox inactive metals in biology. In this chapter, we discuss the periodic table's impact on bio-inorganic chemistry, by exploring reasons behind this selective choice of metals biology. A special focus is placed on the chemical and functional reasons why one metal ion is preferred over another one. We discuss the implications of metal choice in various biological processes including catalysis, electron transfer, redox sensing and signaling. We find that bioavailability of metal ions along with their redox potentials, coordination flexibility, valency and ligand affinity determine the specificity of metals for biological processes. Understanding the implications underlying the selective choice of metals of the periodic table in these biological processes can help design more efficient catalysts, more precise biosensors and more effective drugs.

Fluctuating hydrogen-bond networks govern anomalous electron transfer kinetics in a blue copper protein

We combine experimental and computational methods to address the anomalous kinetics of long-range electron transfer (ET) in mutants of Pseudomonas aeruginosa azurin. ET rates and driving forces for wild type (WT) and three N47X mutants (X = L, S, and D) of Ru(2,2'-bipyridine)₂ (imidazole)(His83) azurin are reported. An enhanced ET rate for the N47L mutant suggests either an increase of the donor-acceptor (DA) electronic coupling or a decrease in the reorganization energy for the reaction. The underlying atomistic features are investigated using a recently developed nonadiabatic molecular dynamics method to simulate ET in each of the azurin mutants, revealing unexpected aspects of DA electronic coupling. In particular, WT azurin and all studied mutants exhibit more DA compression during ET (>2 Å) than previously recognized. Moreover, it is found that DA compression involves an extended network of hydrogen bonds, the fluctuations of which gate the ET reaction, such that DA compression is facilitated by transiently rupturing hydrogen bonds. It is found that the N47L mutant intrinsically disrupts this hydrogen-bond network, enabling particularly facile DA compression. This work, which reveals the surprisingly fluctional nature of ET in azurin, suggests that hydrogen-bond networks can modulate the efficiency of long-range biological ET.

Protein bioelectronics: a review of what we do and do not know

We review the status of protein-based molecular electronics. First, we define and discuss fundamental concepts of electron transfer and transport in and across proteins and proposed mechanisms for these processes. We then describe the immobilization of proteins to solid-state surfaces in both nanoscale and macroscopic approaches, and highlight how different methodologies can alter protein electronic properties. Because immobilizing proteins while retaining biological activity is crucial to the successful development of bioelectronic devices, we discuss this process at length. We briefly discuss computational predictions and their connection to experimental results. We then summarize how the biological activity of immobilized proteins is beneficial for bioelectronic devices, and how conductance measurements can shed light on protein properties. Finally, we consider how the research to date could influence the development of future bioelectronic devices.

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.

Structure, Dynamics, and Electron Transfer of Azurin Bound to a Gold Electrode

Blue copper redox protein azurin (AZ) constitutes an ideal active element for building bionano-optoelectronic devices based on the intriguing interplay among its electron transfer (ET), vibrational, and optical properties. A full comprehension of its dynamical and functional behavior is required for efficient applications. Here, AZ bound to gold electrode via its disulfide bridge was investigated by a molecular dynamics simulation approach taking into account for gold electron polarization which provides a more realistic description of the protein-gold interaction. Upon binding to gold, AZ undergoes slight changes in its secondary structure with the preservation of the copper-containing active site structure. Binding of AZ to gold promotes new collective motions, with respect to free AZ, as evidenced by essential dynamics. Analysis of the ET from the AZ copper ion to the gold substrate, performed by the Pathways model, put into evidence the main residues and structural motifs of AZ involved in the ET paths. During the dynamical evolution of the bionanosystem, transient contacts between some lateral protein atoms and the gold substrate occurred; concomitantly, the opening of additional ET channels with much higher rates was registered. These results provide new and detailed insights on the dynamics and ET properties of the AZ-gold system, by also helping to rationalize some imaging and conductive experimental evidences and also to design new bionanodevices with tailored features.

Extending the essential dynamics analysis to investigate molecular properties: application to the redox potential of proteins

Here, a methodology is proposed to investigate the collective fluctuation modes of an arbitrary set of observables, maximally contributing to the fluctuation of another functionally relevant observable. The methodology, based on the analysis of fully classical molecular dynamics (MD) simulations, exploits the essential dynamics (ED) method, originally developed to analyse the collective motions in proteins. We apply this methodology to identify the residues that are more relevant for determining the reduction potential (E(0)) of a redox-active protein. To this aim, the fluctuation modes of the single-residue electrostatic potentials mostly contributing to the fluctuations of the total electrostatic potential (the main determinant of E(0)) are investigated for wild-type azurin and two of its mutants with a higher E(0). By comparing the results here obtained with a previous study on the same systems [Zanetti-Polzi et al., Org. Biomol. Chem., 2015, 13, 11003] we show that the proposed methodology is able to identify the key sites that determine E(0). This information can be used for a general deeper understanding of the molecular mechanisms on the basis of the redox properties of the proteins under investigation, as well as for the rational design of mutants with a higher or lower E(0). From the results of the present analysis we propose a new azurin mutant that, according to our calculations, shows a further increase of E(0).

Asn47 and Phe114 modulate the inner sphere reorganization energies of type zero copper proteins

The geometric structures and electron transfer properties of type 1 Cu proteins are reasonably understood at the molecular level (E. I. Solomon and R. G. Hadt, Coord. Chem. Rev., 2011, 255, 774-789, J. J. Warren, K. M. Lancaster, J. H. Richards and H. B. Gray, J. Inorg. Biochem., 2012, 115, 119-126). Much understanding of type 1 copper electron transfer reactivity has come from site directed mutagenesis studies. For example, artificial "type zero" Cu-centres constructed in cupredoxin-azurin have showcased the capacity of outer-sphere hydrogen bonding networks to enhance Cu II/I electron transfer reactivity. In this paper, we have elaborated on earlier kinetics and electronic structural studies of type zero Cu by calculating the inner sphere reorganization energies of type 1, type 2, and type zero Cu proteins using density functional theory (DFT). Although the choice of density functionals for copper systems is not straightforward, we have benchmarked the density functionals against the recently reported ESI-PES data for two synthetic copper models (S. Niu, D.-L. Huang, P. D. Dau, H.-T. Liu, L.-S. Wang and T. J. Ichiye, Chem. Theory Comput., 2014, 10, 1283). For the Cu proteins, our calculations predict that changes in the coordination number upon metal reduction lead to large inner sphere reorganization energies for type 2 Cu sites, whereas retention in the coordination number is observed for type zero Cu sites. These variations in the coordination number are modulated by the outer-sphere coordinating residues Asn47 and Phe114, which are involved in hydrogen bonding with the Asp112 side chain.

Design and fine-tuning redox potentials of metalloproteins involved in electron transfer in bioenergetics

Redox potentials are a major contributor in controlling the electron transfer (ET) rates and thus regulating the ET processes in the bioenergetics. To maximize the efficiency of the ET process, one needs to master the art of tuning the redox potential, especially in metalloproteins, as they represent major classes of ET proteins. In this review, we first describe the importance of tuning the redox potential of ET centers and its role in regulating the ET in bioenergetic processes including photosynthesis and respiration. The main focus of this review is to summarize recent work in designing the ET centers, namely cupredoxins, cytochromes, and iron-sulfur proteins, and examples in design of protein networks involved these ET centers. We then discuss the factors that affect redox potentials of these ET centers including metal ion, the ligands to metal center and interactions beyond the primary ligand, especially non-covalent secondary coordination sphere interactions. We provide examples of strategies to fine-tune the redox potential using both natural and unnatural amino acids and native and nonnative cofactors. Several case studies are used to illustrate recent successes in this area. Outlooks for future endeavors are also provided. This article is part of a Special Issue entitled Biodesign for Bioenergetics--the design and engineering of electronic transfer cofactors, proteins and protein networks, edited by Ronald L. Koder and J.L. Ross Anderson.

A dynamical approach to non-adiabatic electron transfers at the bio-inorganic interface

A methodology is proposed to investigate the role of the energy fluctuations, determined by the dynamical evolution of a system, and the role of non-adiabaticity in affecting the kinetic rate of electron transfer reactions at the bio-inorganic interface.

A few key residues determine the high redox potential shift in azurin mutants

The wide range of variability of the reduction potential (E(0)) of blue-copper proteins has been the subject of a large number of studies in the past several years. In particular, a series of azurin mutants have been recently rationally designed tuning E(0) over a very broad range (700 mV) without significantly altering the redox-active site [Marshall et al., Nature, 2009, 462, 113]. This clearly suggests that interactions outside the primary coordination sphere are relevant to determine E(0) in cupredoxins. However, the molecular determinants of the redox potential variability are still undisclosed. Here, by means of atomistic molecular dynamics simulations and hybrid quantum/classical calculations, the mechanisms that determine the E(0) shift of two azurin mutants with high potential shifts are unravelled. The reduction potentials of native azurin and of the mutants are calculated obtaining results in good agreement with the experiments. The analysis of the simulations reveals that only a small number of residues (including non-mutated ones) are relevant in determining the experimentally observed E(0) variation via site-specific, but diverse, mechanisms. These findings open the path to the rational design of new azurin mutants with different E(0).

Sulfonamido tripods: tuning redox potentials via ligand modifications

A series of Fe^II-OH₂ complexes were synthesized with ligands based on the tetradentate sulfonamido tripod N,N',N"-[2,2',2"-nitrilotris(ethane-2,1-diyl)]-tris-({R-Ph}-sulfonamido). These complexes differ by the substituent on the aryl rings and were fully characterized, including their molecular structures via X-ray diffraction methods. All the complexes were five-coordinate with trigonal bipyramidal geometry. A linear correlation was observed between the electronic effects of each ligand, given by the Hammett constants of the para-substituents, and the potential of the Fe^II/Fe^III redox couple, which were determined using cyclic voltammetry. It was found that the range of redox potentials for the complexes spanned approximately 160 mV.

Long-Range Electron Transfer in Engineered Azurins Exhibits Marcus Inverted Region Behavior

The Marcus theory of electron transfer (ET) predicts that while the ET rate constants increase with rising driving force until it equals a reaction's reorganization energy, at higher driving force the ET rate decreases, having reached the Marcus inverted region. While experimental evidence of the inverted region has been reported for organic and inorganic ET reactions as well as for proteins conjugated with ancillary redox moieties, evidence of the inverted region in a "protein-only" system has remained elusive. We herein provide such evidence in a series of nonderivatized proteins. These results may facilitate the design of ET centers for future applications such as advanced energy conversions.