We have coupled hybrid quantum mechanics (density functional theory; Car-Parrinello)/molecular mechanics molecular dynamics simulations to a grand-canonical scheme, to calculate the in situ redox potential of the Cu 2Ű Ű e Ű 3 Cu Ű half reaction in azurin from Pseudomonas aeruginosa. An accurate description at atomistic level of the environment surrounding the metal-binding site and finite-temperature fluctuations of the protein structure are both essential for a correct quantitative description of the electronic properties of this system. We report a redox potential shift with respect to copper in water of 0.2 eV (experimental 0.16 eV) and a reorganization free energy â«Ű⏠0.76 eV (experimental 0.6 -0.8 eV). The electrostatic field of the protein plays a crucial role in fine tuning the redox potential and determining the structure of the solvent. The inner-sphere contribution to the reorganization energy is negligible. The overall small value is mainly due to solvent rearrangement at the protein surface.density functional theory Í electron transfer Í molecular dynamics Í reorganization energy E lectron transfer (ET) processes are ubiquitous chemical reactions that occur in a variety of essential biological functions, such as immune response, respiration, and photosynthesis (1-10). Among the different families of ET proteins, single-copper cupredoxins, also known as blue copper proteins, are particularly appealing for theoretical studies because of their relatively small size and the great availability of experimental measurements in the literature. Blue copper proteins exchange electrons among themselves or with other redox proteins, such as cytochrome c551 or nitrite reductase, using a protein-bound Cu metal ion, that can exist in the Cu(II) or Cu(I) redox states (11,12). The copper ion forms a type 1 Cu-binding site ( Fig. 1), which is characterized by a bright blue color, a narrow hyperfine splitting in the electron paramagnetic resonance (EPR) spectra, a high reduction potential (13-15), and a strong structural similarity between oxidized and reduced states (13-17). In fact, in both oxidation states, the copper ion is coordinated by a cysteine (Cys) thiolate group and two histidine (His) nitrogen atoms in a trigonal planar conformation. The coordination polyhedron is completed by one axial ligand, typically a methionine (Met) thioether group. In azurin, a backbone amide oxygen of a glycine (Gly) constitutes an additional axial ligand (18). The structural similarity between the two redox states provides an advantage for the functionality of these ET proteins, as the reorganization free energy () for the redox process is small [0.6-0.8 eV (17, 19)], allowing a high ET rate (20-23). This is not the case for solvated copper ions and synthetic Cu complexes, which tend to have large structural changes concurrent with any variation in their oxidation state (24). The origin of the low value of is still unclear. The initial hypothesis implied that the rigidity of the protein would force Cu(II) to be bound in a geomet...