The role of zinc finger linkers in zinc finger protein binding to DNA.

Zinc finger proteins (ZFP) play important roles in cellular processes. The DNA binding region of ZFP consists of 3 zinc finger DNA binding domains connected by amino acid linkers, the sequence TGQKP connects ZF1 and ZF2, and TGEKP connects ZF2 with ZF3. Linkers act to tune the zinc finger protein in the right position to bind its DNA target, the type of amino acid residues and length of linkers reflect on ZF1-ZF2-ZF3 interactions and contribute to the search and recognition process of ZF protein to its DNA target. Linker mutations and the affinity of the resulting mutants to specific and nonspecific DNA targets were studied by MD simulations and MM_GB(PB)SA. The affinity of mutants to DNA varied with type and position of amino acid residue. Mutation of K in TGQKP resulted in loss in affinity due to the loss of positive K interaction with phosphates, mutation of G showed loss in affinity to DNA, WT protein and all linker mutants showed loss in affinity to a nonspecific DNA target, this finding confirms previous reports which interpreted this loss in affinity as due to ZF1 having an anchoring role, and ZF3 playing an explorer role in the binding mechanism. The change in ZFP-DNA affinity with linker mutations is discussed in view of protein structure and role of linker residues in binding.

Role of protein structure and the role of individual fingers in zinc finger protein-DNA recognition: a molecular dynamics simulation study and free energy calculations.

Molecular dynamics and MM_GBSA energy calculations on various zinc finger proteins containing three and four fingers bound to their target DNA gave insights into the role of each finger in the DNA binding process as part of the protein structure. The wild type Zif 268 (PDB code: 1AAY) gave a ΔG value of - 76.1 (14) kcal/mol. Zinc fingers ZF1, ZF2 and ZF3 were mutated in one experiment and in another experiment one finger was cut and the rest of the protein was studied for binding. The ΔΔG values for the Zinc Finger protein with both ZF1 and ZF2 mutated was + 80 kcal/mol, while mutating only ZF1 the ΔΔG value was + 52 kcal/mol (relative to the wild type). Cutting ZF3 and studying the protein consisting only of ZF1 linked to ZF2 gave a ΔΔG value of + 68 kcal/mol. Upon cutting ZF1, the resulting ZF2 linked to ZF3 protein gave a ΔΔG value of + 41 kcal/mol. The above results shed light on the importance of each finger in the binding process, especially the role of ZF1 as the anchoring finger followed in importance by ZF2 and ZF3. The energy difference between the binding of the wild type protein Zif268 (1AAY) and that for individual finger binding to DNA according to the formula: ΔΔGlinkers, otherstructuralfactors = ΔGzif268 - (ΔGF1+F2+F3) gave a value = - 44.5 kcal/mol. This stabilization can be attributed to the contribution of linkers and other structural factors in the intact protein in the DNA binding process. DNA binding energies of variant proteins of the wild type Zif268 which differ in their ZF1 amino acid sequence gave evidence of a good relationship between binding energy and recognition and specificity, this finding confirms the reported vital role of ZF1 in the ZF protein scanning and anchoring to the target DNA sequence. The role of hydrogen bonds in both specific and nonspecific amino acid-DNA contacts is discussed in relation to mutations. The binding energies of variant Zinc Finger proteins confirmed the role of ZF1 in the recognition, specificity and anchoring of the zinc finger protein to DNA.

Zinc finger protein binding to DNA: an energy perspective using molecular dynamics simulation and free energy calculations on mutants of both zinc finger domains and their specific DNA bases.

Energy calculations based on MM-GBSA were employed to study various zinc finger protein (ZF) motifs binding to DNA. Mutants of both the DNA bound to their specific amino acids were studied. Calculated energies gave evidence for a relationship between binding energy and affinity of ZF motifs to their sites on DNA. ΔG values were -15.82(12), -3.66(12), and -12.14(11.6) kcal/mol for finger one, finger two, and finger three, respectively. The mutations in the DNA bases reduced the value of the negative energies of binding (maximum value for ΔΔG = 42Kcal/mol for F1 when GCG mutated to GGG, and ΔΔG = 22 kcal/mol for F2, the loss in total energy of binding originated in the loss in electrostatic energies upon mutation (r = .98). The mutations in key amino acids in the ZF motif in positions-1, 2, 3, and 6 showed reduced binding energies to DNA with correlation coefficients between total free energy and electrostatic was .99 and with Van der Waal was .93. Results agree with experimentally found selectivity which showed that Arginine in position-1 is specific to G, while Aspartic acid (D) in position 2 plays a complicated role in binding. There is a correlation between the MD calculated free energies of binding and those obtained experimentally for prepared ZF motifs bound to triplet bases in other reports (), our results may help in the design of ZF motifs based on the established recognition codes based on energies and contributing energies to the total energy.

Binding of the Zn2+ ion to ferric uptake regulation protein from E. coli and the competition with Fe2+ binding: a molecular modeling study of the effect on DNA binding and conformational changes of Fur.

The three dimensional structure of Ferric uptake regulation protein dimer from E. coli, determined by molecular modeling, was docked on a DNA fragment (iron box) and Zn(2+) ions were added in two steps. The first step involved the binding of one Zn2+ ion to what is known as the zinc site which consists of the residues Cys 92, Cys 95, Asp 137, Asp141, Arg139, Glu 140, His 145 and His 143 with an average metal-Nitrogen distance of 2.5 A and metal-oxygen distance of 3.1-3.2 A. The second Zn2+ ion is bound to the iron activating site formed from the residues Ile 50, His 71, Asn 72, Gly 97, Asp 105 and Ala 109. The binding of the second Zn2+ ion strengthened the binding of the first ion as indicated by the shortening of the zinc-residue distances. Fe2+, when added to the complex consisting of 2Zn2+/Fur dimer/DNA, replaced the Zn2+ ion in the zinc site and when a second Fe2+ was added, it replaced the second zinc ion in the iron activating site. The binding of both zinc and iron ions induced a similar change in Fur conformations, but shifted residues closer to DNA in a different manner. This is discussed along with a possible role for the Zn2+ ion in the Fur dimer binding of DNA in its repressor activity.

Iron(II) triggered conformational changes in Escherichia coli fur upon DNA binding: a study using molecular modeling.

In order to identify the Fur dimerization domain, a three-dimensional structure of the ferric uptake regulation protein from Escherichia coli (Fur EC) was determined using homology modeling and energy minimization. The Fur monomer consists of turn- helix -turn motif on the N-terminal domain, followed by another helix-turn-helix-turn motif, and two beta-strands separated by a turn which forms the wing. The C-terminal domain, separated by a long coil from the N-terminal, and consisting of two anti parallel beta strands, and a turn-helix-turn-helix-turn motif. Residues in central domain were found to aid the dimer formation, residues 45-70 as evident in the calculated distances; this region is rich in hydrophobic residues. Most interactions occur between residues Val55, Leu53, Gln52, Glu49 and Tyr56 with closest contacts occurring at residues 49-56. These residues are part of an alpha-helix (alpha(4)) near the N-terminal. Upon raising the Fe(2+) concentration the binding of Fur dimer to DNA was enhanced, this was evident when, the Fur EC dimer was docked onto DNA "iron box" (it was found to bind the AT-rich region) and upon addition of Fe(2+) the helices near the N-terminal bound to the major groove of the DNA. Addition of high Fe(2+) concentration triggered further conformational changes in the Fur dimer as was measured by distances between the two subunits, Fe(2+) mediated the Fur binding to DNA by attaching itself to the DNA. At the same time DNA changed conformation as was evident in the distortion in the backbone and the shrinking of major groove distance from 11.4 to 9.3A. Two major Fe(2+) sites were observed on the C-terminal domain: site 1, the traditional Zn site, the cavity contains the residues Cys92, Cys95, Asp137, Asp141, Arg139, Glu 140, His 145 and His 143 at distances range from 1.3 to 2.2A. Site 2 enclave consists of His71, Ile50, Asn72, Gly97, Asp105 and Ala109 at very close proximity to Fe(2+). The closest contacts between Fur dimer and DNA at the AT-rich region were at residues Ala11, Gly12, Leu13, Pro18 and Arg19 mostly hydrophobic residues near the N-terminal domain. Close contacts repeated at His87, His88 and Arg112, and a third region near the C-terminal at Asn137, Arg 139, Glu140, Asn141, His143, Asn141 and His145. Fur dimer has three major contact regions with DNA, the first on the N-terminal domain, a second smaller region at His87, His88 and Arg112 mediated by Fe(2+) ions, and a third region on the C-terminal domain consisting mainly of hydrophobic contacts and mediated by Fe(2+) ions at high concentration.