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.     

Engineering of a redox protein for DNA-directed assembly

Engineered enzyme conjugate of the small laccase enzyme from Streptomyces coelicolor and zinc finger DNA binding domain from Zif268 is demonstrated to bind double stranded DNA in a site specific manner while retaining enzymatic activity.     

Automatic polymerase chain reaction product detection system for food safety monitoring using zinc finger protein fused to luciferase

An automatic polymerase chain reaction (PCR) product detection system for food safety monitoring using zinc finger (ZF) protein fused to luciferase was developed. ZF protein fused to luciferase specifically binds to target double stranded DNA sequence and has luciferase enzymatic activity. Therefore, PCR products that comprise ZF protein recognition sequence can be detected by measuring the luciferase activity of the fusion protein. We previously reported that PCR products from Legionella pneumophila and Escherichia coli (E. coli) O157 genomic DNA were detected by Zif268, a natural ZF protein, fused to luciferase. In this study, Zif268–luciferase was applied to detect the presence of Salmonella and coliforms. Moreover, an artificial zinc finger protein (B2) fused to luciferase was constructed for a Norovirus detection system. In the luciferase activity detection assay, several bound/free separation process is required. Therefore, an analyzer that automatically performed the bound/free separation process was developed to detect PCR products using the ZF–luciferase fusion protein. By means of the automatic analyzer with ZF–luciferase fusion protein, target pathogenic genomes were specifically detected in the presence of other pathogenic genomes. Moreover, we succeeded in the detection of 10 copies of E. coli BL21 without extraction of genomic DNA by the automatic analyzer and E. coli was detected with a logarithmic dependency in the range of 1.0 × 10 to 1.0 × 10⁶ copies.     

Sugar edge/sugar edge base pairs in RNA: stabilities and structures from quantum chemical calculations

Cis and trans sugar edge/sugar edge (SE/SE) binding patterns are essential building units of RNAs. For example, SE/SE interactions form the A-minor motifs, the most important tertiary interaction type in functional RNAs. This study provides an in-depth structure and stability analysis for these two base pair families. Gas-phase-optimized geometries are reported for 12 cis and 7 trans SE/SE base pairs and contrasted to their X-ray counterparts. Interaction energies are computed at the RIMP2 level of theory using the density-functional-theory-optimized geometries. There is a good overall agreement between the optimized and X-ray geometries of the cis SE/SE base pairs. In contrast, only three of the seven trans SE/SE binding patterns could be optimized without a significant distortion of the X-ray geometry. Note, however, that many SE/SE base pairs participate in broader networks of interactions; thus it is not surprising to see some of them to deviate from the X-ray geometry in a complete isolation. Computed interaction energies reveal that all 12 known cis SE/SE binding patterns are very stable. Among the trans SE/SE binding patterns, only the rG/rG, rG/rC, and rA/rG base pairs are sufficiently stable in the crystal geometry. Prediction has been made for some structures not yet detected by crystallography, namely, cis rC/rC, rG/rC, rG/rU, and rU/rU and trans rG/rA base pairs. Interestingly, the new cis SE/SE binding patterns are not necessarily isosteric with the remaining 12 members of this family. The trans rG/rA base pair represents a viable option for base pairing in RNA to be identified by future X-ray studies. In a complete lack of structural information, prediction of other unknown members of the trans SE/SE family was not attempted. Analysis of the interaction energies shows a very large electron correlation component of the interaction energy, pointing at the elevated role of dispersion energy as compared to other types of base pairs. This likely is profitable for stabilization of SE/SE binding patterns in polar environments and could be one of the reasons why the A-minor motif is the leading type of tertiary interactions in RNAs.     

Targeting retroviral Zn finger-DNA interactions: a small-molecule approach using the electrophilic nature of trans-platinum-nucleobase compounds

Noncovalent interactions are ubiquitous in ternary systems involving metal ions, DNA/RNA, and proteins and represent a structural motif for design of selective inhibitors of biological function. This contribution shows that small molecules containing platinated purine nucleobases mimic the natural DNA(RNA)-tryptophan recognition interaction of zinc finger peptides, specifically the C-terminal finger of HIV NCp7 protein. Interaction with platinum results in Zn ejection from the peptide accompanied by loss of tertiary structure. Targeting the NCp7-DNA interaction for drug design represents a conceptual advance over electrophiles designed for chemical attack on the zinc finger alone. These results demonstrate examples of a new platinum structural class targeting specific biological processes, distinct from the bifunctional DNA-DNA binding of cytotoxic agents like cisplatin. The results confirm the validity of a chemical biological approach for metallodrug design for selective ternary DNA(RNA)-protein interactions.     

Site-selective RNA cleavage by DNA bearing a base pair-mimic nucleoside

We have synthesized the deoxyadenosine derivative tethering a phenyl group (X), which mimics the Watson-Crick A/T base pair. The RNA/DNA hybrid duplexes containing X in the middle of the DNA sequence showed a similar thermal stability regardless of the ribonucleotide species (A, G, C, or U) opposite to X, probably because of the phenyl group stacking inside of the duplex accompanied by the opposite ribonucleotide base flipped in an extrahelical position. The RNA strand hybridized with the DNA strand bearing X was cleaved on the 3'-side of the ribonucleotide opposite to X in the presence of MgCl2, and the RNA sequence to be cleaved was not restricted. The site-specific RNA hydrolysis suggests that the DNA strand bearing X has the advantage of the site-selective base flipping in the target sequence and the development of a "universal deoxyribozyme" to exclusively cleave a target RNA sequence.     

The structure of a mixed LNA/DNA:RNA duplex is driven by conformational coupling between LNA and deoxyribose residues as determined from 13C relaxation measurements

A study of the internal dynamics of an LNA/DNA:RNA duplex has been performed to further characterize the conformational changes associated with the incorporation of locked nucleic acid (LNA) nucleotides in a DNA:RNA duplex. In general, it was demonstrated that the LNA/DNA:RNA duplex has a very high degree of order compared to dsDNA and dsRNA duplexes. The order parameters of the aromatic carbon atoms in the LNA/DNA strand are uniformly high, whereas a sharp drop in the degree of order was seen in the RNA strand in the beginning of the AUAU stretch in the middle of the strand. This can be related to a return to normal dsRNA dynamics for the central A:U base pair. The high order of the heteroduplex is consistent with preorganization of the chimera strand for an A-form duplex conformation. These results partly explain the dramatic increase in T(m) of the chimeric heteroduplex over dsDNA and DNA:RNA hybrids of the same sequence.     

Synthesis,incorporation into triplex-forming oligonucleotide,and binding properties of a novel 2'-deoxy-C-nucleoside featuring a 6-(thiazolyl-5)benzimidazole nucleobase

[reaction: see text] 6-(Thiazolyl-5)benzimidazole (B(t)()) was designed as a novel nucleobase for the specific recognition of an inverted A.T base pair in a triple helix motif. It was successfully incorporated into an 18-mer triplex-forming oligonucleotide (TFO) using the 2'-deoxy-C-nucleoside phosphoramidite 16. The triple helix binding properties of the modified TFO were examined by means of thermal denaturation experiments targeting an oligopyrimidine.oligopurine 26-mer DNA duplex containing an A.T base pair inversion.     

DNA小沟结合药物DB818的分子动力学模拟和结合自由能分析

采用分子动力学模拟了DNA小沟结合芳香二脒药物DB818形成的复合物. 通过5 ns的模拟研究表明: DB818药物分子可紧密结合在DNA的AATTC小沟区域, 和双螺旋d[CGCGAATTCGCG]2形成稳定的复合物. 由于噻吩硫原子的弱电负性, 使DB818能够以更大的伸展程度与DNA的小沟结合, 形成更强的结合力. DB818苯并咪唑的氮原子能够与DNA 7位和19位T碱基上的氧原子形成两个稳定的氢键, 同时, DB818末端氨基氮原子分别与DNA 的20位T碱基的氧原子和9位C碱基的氧原子形成两个氢键. 另外, 运用MM_PBSA方法计算了DB293-DNA和DB818-DNA复合物的结合自由能, 计算结合能与实验值能较好的吻合, 通过比较其结合自由能, 从热力学能量角度说明了DB818有较大的熵值与较小的焓值贡献, 从而与DNA小沟结合的结合力比DB293强. 本文在分子水平上提供了DB818直接与双螺旋DNA相互作用的结构及复合物的动态变化情况, 为设计出更高生物活性的DNA小沟结合剂提供一定的理论依据.     

Hg(II) ion specifically binds with T:T mismatched base pair in duplex DNA

Metal-mediated base pair formation, resulting from the interaction between metal ions and artificial bases in oligonucleotides, has been developed for its potential application in nanotechnology. We have recently found that the T:T mismatched base pair binds with Hg(II) ions to generate a novel metal-mediated base pair in duplex DNA. The thermal stability of the duplex with the T-Hg-T base pair was comparable to that of the corresponding T:A or A:T. The novel T-Hg-T base pair involving the natural base thymine is more convenient than the metal-mediated base pairs involving artificial bases due to the lack of time-consuming synthesis. Here, we examine the specificity and thermodynamic properties of the binding between Hg(II) ions and the T:T mismatched base pair. Only the melting temperature of the duplex with T:T and not of the perfectly matched or other mismatched base pairs was found to specifically increase in the presence of Hg(II) ions. Hg(II) specifically bound with the T:T mismatched base pair at a molar ratio of 1:1 with a binding constant of 10(6) M(-1), which is significantly higher than that for nonspecific metal ion-DNA interactions. Furthermore, the higher-order structure of the duplex was not significantly distorted by the Hg(II) ion binding. Our results support the idea that the T-Hg-T base pair could eventually lead to progress in potential applications of metal-mediated base pairs in nanotechnology.     

Modelling CreA protein–DNA recognition determinants. A molecular dynamics study of fully charged CreA–DNA model in water

CreA, the negative regulator mediating carbon catabolism repression in Aspergillus nidulans, is a protein that contains a DNA-binding domain comprising two zinc finger motifs. A 3D model for the CreA–G4 (5′-GCGGGGGCGT-3′) complex is constructed on the basis of the structure of the Zif268–DNA crystal complex [Science 252 (1991) 809] and using similarity analysis and computer assisted modelling techniques. The CreA–G4 model was then subjected to a set of molecular dynamics (MD) studies. Based on our previous nano second long Zif268–DNA MD simulation, a 170 ps long trajectory was deemed sufficient to test possible DNA–protein interactions. A screening of the static model and of the trajectory was performed for protein amino acids, nucleotide bases, phosphates backbone and water molecules mediating protein–DNA contacts. Energy, root mean square deviation (RMSD), principal inertial moment and distances were analysed. For this time span, the stability shown in fluctuation patterns reveals the presence of a complex behaving in a manner similar to Zif268–DNA.

An unambiguous characterisation of the amino acids involved in DNA-binding was obtained. These results could contribute towards the establishment of a code of protein–DNA recognition for this class of DNA-binding motifs.     

碱基对在DNA双螺旋链上分离的自由能计算

DNA是大部分生物包括病毒的基因载体。DNA双螺旋链通过A＝T和G≡C两种碱基对编码实现对遗传信息的存储。碱基对中的相互作用对DNA双螺旋链的稳定性起到重要作用,直接关系到基因的复制和转录。当前研究中,我们构建了四组不同结构的DNA双螺旋链,进行了总共4.3 μs的分子动力学模拟。通过伞形取样技术计算了DNA双螺旋链中碱基对分离的自由能曲线,并从分子尺度细节和相互作用能对自由能曲线进行解析。在碱基对G≡C的自由能曲线(PMF-PGC)上观察到三个峰,通过监测氢键数目的变化发现分别对应于G≡C三个氢键的断裂;而在A＝T的自由能曲线(PMF-PAT)上只出现一个峰,说明A＝T的两个氢键在分离过程中几乎同时断裂。PMF-PGC的总能垒比PMF-PAT高,主要是因为G≡C比A＝T多一个氢键,更稳定。两条曲线的后段自由能仍然升高,而此时碱基对的氢键已断裂,这是DNA链骨架刚性所导致。我们还研究了碱基对稳定性受相邻碱基对的影响,发现邻近G≡C碱基对会增强A＝T的稳定性, C≡G会削弱A＝T的稳定性, T＝A对A＝T的影响较小。     

Efficient incorporation of a copper hydroxypyridone base pair in DNA

Recently, we reported the first artificial nucleoside for alternative DNA base pairing through metal complexation (J. Org. Chem. 1999, 64, 5002-5003). In this regard, we report here the synthesis of a hydroxypyridone-bearing nucleoside and the incorporation of a neutral Cu(2+)-mediated base pair of hydroxypyridone nucleobases (H-Cu-H) in a DNA duplex. When the hydroxypyridone bases are incorporated into the middle of a 15 nucleotide duplex, the duplex displays high thermal stabilization in the presence of equimolar Cu(2+) ions in comparison with a duplex containing an A-T pair in place of the H-H pair. Monitoring temperature dependence of UV-absorption changes verified that a Cu(2+)-mediated base pair is stoichiometrically formed inside the duplex and dissociates upon thermal denaturation at elevated temperature. In addition, EPR and CD studies suggested that the radical site of a Cu(2+) center is formed within the right-handed double-strand structure of the oligonucleotide. The present strategy could be developed for controlled and periodic spacing of neutral metallobase pairs along the helix axis of DNA.     

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     

Altered sequence specificity identified from a library of DNA-binding small molecules

BACKGROUND: The ability to target specific DNA sequences using small molecules has major implications for basic research and medicine. Previous studies revealed that a bis-intercalating molecule containing two 1,4,5,8-napthalenetetracarboxylic diimides separated by a lysine-tris-glycine linker binds to DNA cooperatively, in pairs, with a preference for G + C-rich sequences. Here we investigate the binding properties of a library of bis-intercalating molecules that have partially randomized peptide linkers. RESULTS: A library of bis-intercalating derivatives with varied peptide linkers was screened for sequence specificity using DNase I footprinting on a 231 base pair (bp) restriction fragment. The library mixtures produced footprints that were generally similar to the parent bis-intercalator, which bound within a 15 bp G + C-rich repeat above 125 nM. Nevertheless, subtle differences in cleavage enhancement bands followed by library deconvolution revealed a derivative with novel specificity. A lysine-tris-beta-alanine derivative was found to bind preferentially within a 19 bp palindrome, without substantial loss of affinity. CONCLUSIONS: Synthetically simple changes in the bis-intercalating compounds can produce derivatives with novel sequence specificity. The large size and symmetrical nature of the preferred binding sites suggest that cooperativity may be retained despite modified sequence specificity. Such findings, combined with structural data, could be used to develop versatile DNA ligands of modest molecular weight that target relatively long DNA sequences in a selective manner.     

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.     

Cyclohexyl "base pairs" stabilize duplexes and intensify pyrene fluorescence by shielding it from natural base pairs

In this study, we investigated the stability and structure of artificial base pairs that contain cyclohexyl rings. The introduction of a single pair of isopropylcyclohexanes into the middle of DNA slightly destabilized the duplex. Interestingly, as the number of the "base pairs" increased, the duplex was remarkably stabilized. A duplex with six base pairs was even more stable than one containing six A-T pairs. Thermodynamic analysis revealed that changes in entropy and not enthalpy contributed to duplex stability, demonstrating that hydrophobic interactions between isopropyl groups facilitated the base pairing, and thus stabilized the duplex. NOESY of a duplex containing an isopropylcyclohexane-methylcyclohexane pair unambiguously demonstrated its "pairing" in the duplex because distinct NOEs between the protons of cyclohexyl moieties and imino protons of both of the neighboring natural base pairs were observed. CD spectra of duplexes tethering cyclohexyl moieties also showed a positive-negative couplet that is characteristic of the B-form DNA duplex. Taken together, these results showed that cyclohexyl moieties formed base pairs in the DNA duplex without severely disturbing the helical structure of natural DNA. Next, we introduced cyclohexyl base pairs between pyrene and nucleobases as an "insulator" that suppresses electron transfer between them. We found a massive increase in the quantum yield of pyrene due to the efficient shielding of pyrene from nucleobases. The cyclohexyl base pairs reported here have the potential to prepare highly fluorescent labeling agents by multiplying fluorophores and insulators alternately into DNA duplexes.