Classical zinc fingers domains (ZFs) bind Zn(II) ion by a pair of cysteine and histidine residues to adopt a characteristic and stable ββα fold containing a small hydrophobic core. As a component of transcription factors, they recognize specific DNA sequences to transcript particular genes. The loss of Zn(II) disrupts the unique structure and function of the whole protein. It has been shown that the saturation of ZFs under cellular conditions is strictly related to their affinity for Zn(II). High affinity warrants their constant saturation, while medium affinity results in their transient structurization depending on cellular zinc availability. Therefore, there must be factors hidden in the sequence and structure of ZFs that impact Zn(II)-to-protein affinities to control their function. Using molecular dynamics simulations and experimental spectroscopic and calorimetric approaches, we showed that particular non-conserved residues derived from ZF sequences impact hydrogen bond formation. Our in silico and in vitro studies show that non-conserved residues can alter metal-coupled folding mechanisms and overall ZF stability. Furthermore, we show that Zn(II) binding to ZFs can also be entropically driven. This preference does not correlate either with Zn(II) binding site or with the extent of the secondary structure but is strictly related to a reservoir of interactions within the second coordination shell, which may loosen or tighten up the structure. Our findings shed new light on how the functionality of ZFs is modulated by non-coordinating residues diversity under cellular conditions. Moreover, they can be helpful for systematic backbone alteration of native ZF ββα scaffold to create artificial foldamers and proteins with improved stability.
Keywords: DNA; energetics; hydrogen bond; isothermal titration calorimetry; metal binding affinity; metal-coupled folding; molecular dynamics.