Separating the Role of Protein Restraints and Local Metal-Site Interaction Chemistry in the Thermodynamics of a Zinc Finger Protein

Abstract

We express the effective Hamiltonian of an ion-binding site in a protein as a combination of the Hamiltonian of the ion-bound site in vacuum and the restraints of the protein on the site. The protein restraints are described by the quadratic elastic network model. The Hamiltonian of the ion-bound site in vacuum is approximated as a generalized Hessian around the minimum energy configuration. The resultant of the two quadratic Hamiltonians is cast into a pure quadratic form. In the canonical ensemble, the quadratic nature of the resultant Hamiltonian allows us to express analytically the excess free energy, enthalpy, and entropy of ion binding to the protein. The analytical expressions allow us to separate the roles of the dynamic restraints imposed by the protein on the binding site and the temperature-independent chemical effects in metal-ligand coordination. For the consensus zinc-finger peptide, relative to the aqueous phase, the calculated free energy of exchanging Zn2+ with Fe2+, Co2+, Ni2+, and Cd2+ are in agreement with experiments. The predicted excess enthalpy of ion exchange between Zn2+ and Co2+ also agrees with the available experimental estimate. The free energy of applying the protein restraints reveals that relative to Zn2+, the Co2+, and Cd2+-site clusters are more destabilized by the protein restraints. This leads to an experimentally testable hypothesis that a tetrahedral metal binding site with minimal protein restraints will be less selective for Zn2+ over Co2+ and Cd2+ compared to a zinc finger peptide. No appreciable change is expected for Fe2+ and Ni2+. The framework presented here may prove useful in protein engineering to tune metal selectivity.