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: ΔΔG_{linkers, otherstructuralfactors}?=?ΔG_zif268???(ΔG_F1+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.     

Expanding the DNA-recognition repertoire for zinc finger proteins beyond 20 amino acids

The design of DNA binding domains based on the Cys2His2 zinc finger motif has proven to be a successful strategy for the specific recognition of novel DNA sequences. Although considerable effort has been devoted to the generation of zinc finger proteins with widely varying DNA-binding preferences, only a limited number of potential DNA binding sites have been targeted with a high degree of specificity. These restrictions on zinc finger design appear to be a consequence of the limited repertoire of side-chain lengths and functionalities available with the 20 proteinogenic amino acids. To demonstrate that these limitations can be overcome through the use of "unnatural" amino acids, expressed protein ligation was employed to incorporate the amino acid citrulline into a single position within a three-zinc finger protein. As anticipated, the resulting semisynthetic protein specifically recognizes adenine in the appropriate position of its binding site.     

Interaction identification of Zif268 and TATAZF proteins with GC‐/AT‐rich DNA sequence: A theoretical study

Molecular dynamics (MD) simulations for Zif268 (a zinc‐finger‐protein binding specifically to the GC‐rich DNA)‐d(A₁G₂C₃G₄T₅G₆G₇G₈C₉A₁₀C₁₁)₂ and TATA_ZF (a zinc‐finger‐protein recognizing the AT‐rich DNA)‐d(A₁C₂G₃C₄T₅A₆T₇A₈A₉A₁₀A₁₁G₁₂G₁₃)₂ complexes have been performed for investigating the DNA binding affinities and specific recognitions of zinc fingers to GC‐rich and AT‐rich DNA sequences. The binding free energies for the two systems have been further analyzed by using the molecular mechanics Poisson‐Boltzmann surface area (MM‐PBSA) method. The calculations of the binding free energies reveal that the affinity energy of Zif268‐DNA complex is larger than that of TATA_ZF‐DNA one. The affinity between the zinc‐finger‐protein and DNA is mainly driven by more favorable van‐der‐Waals and nonpolar/solvation interactions in both complexes. However, the affinity energy difference of the two binding systems is mainly caused by the difference of van‐der‐Waals interactions and entropy components. The decomposition analysis of MM‐PBSA free energies on each residue of the proteins predicts that the interactions between the residues with the positive charges and DNA favor the binding process; while the interactions between the residues with the negative charges and DNA behave in the opposite way. The interhydrogen‐bonds at the protein‐DNA interface and the induced intrafinger hydrogen bonds between the residues of protein for the Zif268‐DNA complex have been identified at some key contact sites. However, only the interhydrogen‐bonds between the residues of protein and DNA for TATA_ZF‐DNA complex have been found. The interactions of hydrogen‐bonds, electrostatistics and van‐der‐Waals type at some new contact sites have been identified. Moreover, the recognition characteristics of the two studied zinc‐finger‐proteins have also been discussed. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011     

Can Co(II) or Cd(II) substitute for Zn(II) in zinc fingers?

Zinc finger domains consist of sequences of amino acids containing cysteine and histidine residues tetrahedrally coordinated to a zinc ion. The role of zinc in a DNA binding finger was considered purely structural due to the absence of redox chemistry in zinc. However, whether other metals e.g. Co(II) or Cd(II) can substitute Zn(II) is not settled. For an answer the detailed interaction of Co(II) and Cd(II) with cysteine methylester and histidine methylester has been investigated as a model for the zinc core in zinc fingers. The study was extended to different temperatures to evaluate the thermodynamic parameters associated with these interactions. The results suggest that zinc has a unique role.     

Control of thiolate nucleophilicity and specificity in zinc metalloproteins by hydrogen bonding: lessons from model compound studies

A single hydrogen bond between an amide N-H and a thiolate sulfur in model complexes designed to mimic the binding site of zinc thiolate proteins, is shown to reduce the reactivity of the thiolate toward electrophiles by up to 2 orders of magnitude. In addition a single such bond is also sufficient to achieve nearly 100% regiospecificity of reaction between a strong, and hence inherently indiscriminate, alkylating agent like trimethyl oxonium tetrafluoroborate and a single sulfur in a dithiolate construct. The importance of these results in understanding how two systems such as the zinc fingers of the GATA family and the Escherichia coli DNA repair protein Ada which share the same pseudotetrahedral structure and tetrascysteinyl ligation around the zinc can fulfill such widely divergent (structural vs reactive) roles and how specificity of reaction in multithiolate-containing systems can be achieved is discussed.     

H-bonding interactions and control of thiolate nucleophilicity and specificity in model complexes of zinc metalloproteins

It is shown in model complexes designed to mimic the binding site of zinc-thiolate proteins that a single hydrogen bond between an amide N-H and a Zn-coordinated thiolate reduces its reactivity toward electrophiles by up to 2 orders of magnitude. In addition, we show that a single N-H...S hydrogen bond is sufficient to achieve near 100% regiospecificity of reaction between a strong, and hence inherently indiscriminate, alkylating agent like trimethyloxonium tetraflouroborate and a single sulfur in a dithiolate construct. The importance of these results in understanding how systems such as the zinc fingers of the GATA family and the E. coli DNA repair protein Ada, which share the same pseudotetrahedral structure and tetracysteinyl ligation around the zinc, can fulfill such widely divergent (structural vs reactive) roles and how specificity of reaction in such multi-thiolate containing systems can be achieved is discussed.     

A bacterial one-hybrid selection system for interrogating zinc finger-DNA interactions

We have developed two bacterial one-hybrid systems for interrogating and selecting zinc finger-DNA interactions. Our systems utilize two plasmids: a zinc finger-plasmid containing the gene for the zinc finger fused to a fragment of the alpha subunit of RNA polymerase and a reporter plasmid where the zinc finger-binding site is located upstream of a reporter gene-either the gene encoding the green fluorescent protein (GFP) or chloramphenicol acetyltransferase (CAT). Binding of the zinc finger domain to the target binding site results in a 10-fold increase in chloramphenicol resistance with the CAT reporter and an 8- to 22-fold increase in total cell fluorescence with the GFP reporter. The CAT reporter allows for sequence specific zinc fingers to be isolated in a single selection step whereas the GFP reporter enables quantitative evaluation of libraries using flow cytometry and theoretically allows for both negative and positive selection. Both systems have been used to select for zinc fingers that have affinity for the motif 5'-GGGGCAGAA-3' from a library of approximately 2 x 10(5) variants. The systems have been engineered to report on zinc finger-DNA binding with dissociation constants less than about 1 microM in order to be most applicable for evaluating binding specificity in an in vivo setting.     

A general strategy to determine a target DNA sequence of a short peptide: application to a d-peptide.

Short peptides could potentially provide a novel element to read-out DNA sequences from the major groove. However, it is difficult to determine sequence-preference of de novo designed monomeric short peptides. Because DNS-binding affinity and specificity of short peptides are usually much lower than those of native DNA-binding proteins, determining the sequence-preference of short peptides by conventional methods utilized to deduce the target sequence of proteins often produces an unclear outcome. We report here a general strategy to defining the sequence-preference of a DNA-binding short peptide by using the heterodimers. A GCN4 basic region peptide tethers a low-affinity DNA-binding peptide adjacent to a GCN4 binding sequence through the cyclodextrin-adamantane association, thereby increasing local concentration of the low-affinity peptide on degenerated DNA sequences. An increase of the local concentration allows one to select a preferential sequence for the low-affinity DNA binding peptide. The method successfully identified specific sequences of short peptides derived from native DNA-binding proteins. The usefulness of this approach has been demonstrated by identifying preferred DNA targets for a peptide composed only of d-amino acids. The method is potentially applicable not only to artificial peptides, but also to other synthethic ligands.     

Cysteine and histidine shuffling: mixing and matching cysteine and histidine residues in zinc finger proteins to afford different folds and function

Zinc finger proteins utilize zinc for structural purposes: zinc binds to a combination of cysteine and histidine ligands in a tetrahedral coordination geometry facilitating protein folding and function. While much is known about the classical zinc finger proteins, which utilize a Cys(2)His(2) ligand set to coordinate zinc and fold into an anti-parallel beta sheet/alpha helical fold, there are thirteen other families of 'non-classical' zinc finger proteins for which relationships between metal coordination and protein structure/function are less defined. This 'Perspective' article focuses on two classes of these non-classical zinc finger proteins: Cys(3)His type zinc finger proteins and Cys(2)His(2)Cys type zinc finger proteins. These proteins bind zinc in a tetrahedral geometry, like the classical zinc finger proteins, yet they adopt completely different folds and target different oligonucleotides. Our current understanding of the relationships between ligand set, metal ion, fold and function for these non-classical zinc fingers is discussed.     

Coordination properties of zinc finger peptides revisited: ligand competition studies reveal higher affinities for zinc and cobalt

Zinc fingers are ubiquitous small protein domains which have a Zn(Cys)(4-x)(His)(x) site. They possess great diversity in their structure and amino acid composition. Using a family of six peptides, it was possible to assess the influence of hydrophobic amino acids on the metal-peptide affinities and on the rates of metal association and dissociation. A model of a treble-clef zinc finger, a model of the zinc finger site of a redox-switch protein, and four variants of the classical ββα zinc finger were used. They differ in their coordination set, their sequence length, and their hydrophobic amino acid content. The speciation, metal binding constants, and structure of these peptides have been investigated as a function of pH. The zinc binding constants of peptides, which adopt a well-defined structure, were found to be around 10(15) at pH 7.0. The rates of zinc exchange between EDTA and the peptides were also assessed. We evidenced that the packing of hydrophobic amino acids into a well-defined hydrophobic core can have a drastic influence on both the binding constant and the kinetics of metal exchange. Notably, well-packed hydrophobic amino acids can increase the stability constant by 4 orders of magnitude. The half-life of zinc exchange was also seen to vary significantly depending on the sequence of the zinc finger. The possible causes for this behavior are discussed. This work will help in understanding the dynamics of zinc exchange in zinc-containing proteins.     

Classical Cys2His2 zinc finger peptides are rapidly oxidized by either H2O2 or O2 irrespective of metal coordination

ZIF268, a member of the classical zinc finger protein family, contains three Cys(2)His(2) zinc binding domains that together recognize the DNA sequence 5'-AGCGTGGGCGT-3'. These domains can be fused to an endonuclease to make a chimeric protein to target and cleave specific DNA sequences. A peptide corresponding to these domains, named ZIF268-3D, has been prepared to determine if the zinc finger domain itself can promote DNA cleavage when a redox active metal ion, Fe(II), is coordinated. The UV-vis absorption spectrum of Fe(II)-ZIF268-3D is indicative of Fe(II) coordination. Using fluorescence anisotropy, we demonstrate that Fe(II)-ZIF268-3D binds selectively to its target DNA in the same manner as Zn(II)-ZIF268-3D. In the presence of added oxidant, H(2)O(2) or O(2), DNA cleavage is not observed by Fe(II)-ZIF268-3D. Instead, the peptide itself is rapidly oxidized. Similarly, Zn(II)-ZIF268-3D and apo-ZIF268-3D are rapidly oxidized by H(2)O(2) or O(2), and we propose that ZIF268-3D is highly susceptible to oxidation.     

Semi-Synthesis and Analysis of Chemically Modified Zif268 Zinc-Finger Domains

Total synthesis of proteins can be challenging despite assembling techniques, such as native chemical ligation (NCL) and expressed protein ligation (EPL). Especially, the combination of recombinant protein expression and chemically addressable solid-phase peptide synthesis (SPPS) is well suited for the redesign of native protein structures. Incorporation of analytical probes and artificial amino acids into full-length natural protein domains, such as the sequence-specific DNA binding zinc-finger motifs, are of interest combining selective DNA recognition and artificial function. The semi-synthesis of the natural 90 amino acid long sequence of the zinc-finger domain of Zif268 is described including various chemically modified constructs. Our approach offers the possibility to exchange any amino acid within the third zinc finger. The realized modifications of the natural sequence include point mutations, attachment of a fluorophore, and the exchange of amino acids at different positions in the zinc finger by artificial amino acids to create additional metal binding sites. The individual constructs were analyzed by circular dichroism (CD) spectroscopy with respect to the integrity of the zinc-finger fold and DNA binding.     

Analysis and prediction of hydrogen bonding in protein-DNA complexes using parallel processors

A number of essential biological functions are controlled by proteins that bind to specific sequences in genomic DNA. In this article we present a simplified model for analyzing DNA-protein interactions mediated exclusively by hydrogen bonds. Based on this model, an optimized algorithm for geometric pattern recognition was developed. The large number of local energy minima are efficiently screened by using a geometric approach to pattern matching based on a square-well potential. The second part of the algorithm represents a closed form solution for minimization based on a quadratic potential. A Monte Carlo method applied to a modified Lennard-Jones potential is used as a third step to rank DNA sequences in terms of pattern matching. Using protein structures derived from four DNA-protein complexes with three-dimensional coordinates established by X-ray diffraction analysis, all possible DNA sequences to which these proteins could bind were ranked in terms of binding energies. The algorithm predicts the correct DNA sequence when at least two hydrogen bonds per base pair are involved in binding to the protein, providing a partial solution to the three-dimensional docking problem. This study lays a framework for future refinements of the algorithm in which the number of assumptions made in the present analysis are reduced. © 1996 by John Wiley & Sons, Inc.