Active site loop dictates the thermodynamics of reduction and ligand protonation in cupredoxins

Molecular dynamics simulations of apocupredoxins: insights into the formation and stabilization of copper sites under entatic control

Cupredoxins perform copper-mediated long-range electron transfer (ET) in biological systems. Their copper-binding sites have evolved to force copper ions into ET-competent systems with decreased reorganization energy, increased reduction potential, and a distinct electronic structure compared with those of non-ET-competent copper complexes. The entatic or rack-induced state hypothesis explains these special properties in terms of the strain that the protein matrix exerts on the metal ions. This idea is supported by X-ray structures of apocupredoxins displaying "closed" arrangements of the copper ligands like those observed in the holoproteins; however, it implies completely buried copper-binding atoms, conflicting with the notion that they must be exposed for copper loading. On the other hand, a recent work based on NMR showed that the copper-binding regions of apocupredoxins are flexible in solution. We have explored five cupredoxins in their "closed" apo forms through molecular dynamics simulations. We observed that prearranged ligand conformations are not stable as the X-ray data suggest, although they do form part of the dynamic landscape of the apoproteins. This translates into variable flexibility of the copper-binding regions within a rigid fold, accompanied by fluctuations of the hydrogen bonds around the copper ligands. Major conformations with solvent-exposed copper-binding atoms could allow initial binding of the copper ions. An eventual subsequent incursion to the closed state would result in binding of the remaining ligands, trapping the closed conformation thanks to the additional binding energy and the fastening of noncovalent interactions that make up the rack.

How the dynamics of the metal-binding loop region controls the acid transition in cupredoxins

Many reduced cupredoxins undergo a pH-dependent structural rearrangement, triggered by protonation of the His ligand belonging to the C-terminal hydrophobic loop, usually termed the acid transition. At variance with several members of the cupredoxin family, the acid transition is not observed for azurin (AZ). We have addressed this issue by performing molecular dynamics simulations of AZ and four mutants, in which the C-terminal loop has been replaced with those of other cupredoxins or with polyalanine loops. All of the loop mutants undergo the acid transition in the pH range of 4.4-5.5. The main differences between AZ and its loop mutants are the average value of the active site solvent accessible surface area and the extent of its fluctuations with time, together with an altered structure of the water layer around the copper center. Using functional mode analysis, we found that these variations arise from changes in nonbonding interactions in the second coordination sphere of the copper center, resulting from the loop mutation. Our results strengthen the view that the dynamics at the site relevant for function and its surroundings are crucial for protein activity and for metal-containing electron transferases.

The Active Site Loop Modulates the Reorganization Energy of Blue Copper Proteins by Controlling the Dynamic Interplay with Solvent

Understanding the factors governing the rate of electron transfer processes in proteins is crucial not only to a deeper understanding of redox processes in living organisms but also for the design of efficient devices featuring biological molecules. Here, molecular dynamics simulations performed on native azurin and four chimeric cupredoxins allow for the calculation of the reorganization energy and of structure-related quantities that were used to clarify the molecular determinants to the dynamics/function relationship in blue copper proteins. We find that the dynamics of the small, metal-binding loop region controls the outer-sphere reorganization energy not only by determining the exposure of the active site to solvent but also through the modulation of the redox-dependent rearrangement of the whole protein scaffold and of the surrounding water molecules.

pH and solvent H/D isotope effects on the thermodynamics and kinetics of electron transfer for electrode-immobilized native and urea-unfolded stellacyanin

The thermodynamics of Cu(II) to Cu(I) reduction and the kinetics of the electron transfer (ET) process for Rhus vernicifera stellacyanin (STC) immobilized on a decane-1-thiol coated gold electrode have been measured through cyclic voltammetry at varying pH and temperature, in the presence of urea and in D(2)O. Immobilized STC undergoes a limited conformational change that mainly results in an enhanced exposure of one or both copper binding histidines to solvent which slightly stabilizes the cupric state and increases histidine basicity. The large immobilization-induced increase in the pK(a) for the acid transition (from 4.5 to 6.3) makes this electrode-SAM-protein construct an attractive candidate as a biomolecular ET switch operating near neutral pH in molecular electronics. Such a potential interest is increased by the robustness of this interface against chemical unfolding as it undergoes only moderate changes in the reduction thermodynamics and in the ET rate in the presence of up to 8 M urea. The sensitivity of these parameters to solvent H/D isotope effects testifies to the role of protein solvation as effector of the thermodynamics and kinetics of ET.