A rapid and sensitive high-throughput screening method to identify compounds targeting protein–nucleic acids interactions

DNA-binding and RNA-binding proteins are usually considered ‘undruggable’ partly due to the lack of an efficient method to identify inhibitors from existing small molecule repositories. Here we report a rapid and sensitive high-throughput screening approach to identify compounds targeting protein–nucleic acids interactions based on protein–DNA or protein–RNA interaction enzyme-linked immunosorbent assays (PDI-ELISA or PRI-ELISA). We validated the PDI-ELISA method using the mammalian high-mobility-group protein AT-hook 2 (HMGA2) as the protein of interest and netropsin as the inhibitor of HMGA2–DNA interactions. With this method we successfully identified several inhibitors and an activator for HMGA2–DNA interactions from a collection of 29 DNA-binding compounds. Guided by this screening excise, we showed that netropsin, the specific inhibitor of HMGA2–DNA interactions, strongly inhibited the differentiation of the mouse pre-adipocyte 3T3-L1 cells into adipocytes, most likely through a mechanism by which the inhibition is through preventing the binding of HMGA2 to the target DNA sequences. This method should be broadly applicable to identify compounds or proteins modulating many DNA-binding or RNA-binding proteins.     

Protein–DNA binding: complexities and multi-protein codes

Binding of proteins to particular DNA sites across the genome is a primary determinant of specificity in genome maintenance and gene regulation. DNA-binding specificity is encoded at multiple levels, from the detailed biophysical interactions between proteins and DNA, to the assembly of multi-protein complexes. At each level, variation in the mechanisms used to achieve specificity has led to difficulties in constructing and applying simple models of DNA binding. We review the complexities in protein–DNA binding found at multiple levels and discuss how they confound the idea of simple recognition codes. We discuss the impact of new high-throughput technologies for the characterization of protein–DNA binding, and how these technologies are uncovering new complexities in protein–DNA recognition. Finally, we review the concept of multi-protein recognition codes in which new DNA-binding specificities are achieved by the assembly of multi-protein complexes.     

From Nonspecific DNA–Protein Encounter Complexes to the Prediction of DNA–Protein Interactions

DNA–protein interactions are involved in many essential biological activities. Because there is no simple mapping code between DNA base pairs and protein amino acids, the prediction of DNA–protein interactions is a challenging problem. Here, we present a novel computational approach for predicting DNA-binding protein residues and DNA–protein interaction modes without knowing its specific DNA target sequence. Given the structure of a DNA-binding protein, the method first generates an ensemble of complex structures obtained by rigid-body docking with a nonspecific canonical B-DNA. Representative models are subsequently selected through clustering and ranking by their DNA–protein interfacial energy. Analysis of these encounter complex models suggests that the recognition sites for specific DNA binding are usually favorable interaction sites for the nonspecific DNA probe and that nonspecific DNA–protein interaction modes exhibit some similarity to specific DNA–protein binding modes. Although the method requires as input the knowledge that the protein binds DNA, in benchmark tests, it achieves better performance in identifying DNA-binding sites than three previously established methods, which are based on sophisticated machine-learning techniques. We further apply our method to protein structures predicted through modeling and demonstrate that our method performs satisfactorily on protein models whose root-mean-square Cα deviation from native is up to 5 Å from their native structures. This study provides valuable structural insights into how a specific DNA-binding protein interacts with a nonspecific DNA sequence. The similarity between the specific DNA–protein interaction mode and nonspecific interaction modes may reflect an important sampling step in search of its specific DNA targets by a DNA-binding protein.     

High-resolution profiling of homing endonuclease binding and catalytic specificity using yeast surface display

Experimental analysis and manipulation of protein–DNA interactions pose unique biophysical challenges arising from the structural and chemical homogeneity of DNA polymers. We report the use of yeast surface display for analytical and selection-based applications for the interaction between a LAGLIDADG homing endonuclease and its DNA target. Quantitative flow cytometry using oligonucleotide substrates facilitated a complete profiling of specificity, both for DNA-binding and catalysis, with single base pair resolution. These analyses revealed a comprehensive segregation of binding specificity and affinity to one half of the pseudo-dimeric interaction, while the entire interface contributed specificity at the level of catalysis. A single round of targeted mutagenesis with tandem affinity and catalytic selection steps provided mechanistic insights to the origins of binding and catalytic specificity. These methods represent a dynamic new approach for interrogating specificity in protein–DNA interactions.     

Mechanism of high-mobility group protein B enhancement of progesterone receptor sequence-specific DNA binding

The DNA-binding domain (DBD) of progesterone receptor (PR) is bipartite containing a zinc module core that interacts with progesterone response elements (PRE), and a short flexible carboxyl terminal extension (CTE) that interacts with the minor groove flanking the PRE. The chromosomal high-mobility group B proteins (HMGB), defined as DNA architectural proteins capable of bending DNA, also function as auxiliary factors that increase the DNA-binding affinity of PR and other steroid receptors by mechanisms that are not well defined. Here we show that the CTE of PR contains a specific binding site for HMGB that is required for stimulation of PR-PRE binding, whereas the DNA architectural properties of HMGB are dispensable. Specific PRE DNA inhibited HMGB binding to the CTE, indicating that DNA and HMGB–CTE interactions are mutually exclusive. Exogenous CTE peptide increased PR-binding affinity for PRE as did deletion of the CTE. In a PR-binding site selection assay, A/T sequences flanking the PRE were enriched by HMGB, indicating that PR DNA-binding specificity is also altered by HMGB. We conclude that a transient HMGB–CTE interaction alters a repressive conformation of the flexible CTE enabling it to bind to preferred sequences flanking the PRE.     

Highly conserved amino acids in Pax and Ets proteins are required for DNA binding and ternary complex assembly

Combinatorial association of DNA-binding proteins on composite binding sites enhances their nucleotide sequence specificity and functional synergy. As a paradigm for these interactions, Pax-5 (BSAP) assembles ternary complexes with Ets proteins on the B cell-specific mb-1 promoter through interactions between their respective DNA-binding domains. Pax-5 recruits Ets-1 to bind the promoter, but not the closely related Ets protein SAP1a. Here we show that, while several different mutations increase binding of SAP1a to an optimized Ets binding site, only conversion of Val68 to an acidic amino acid facilitates ternary complex assembly with Pax-5 on the mb-1 promoter. This suggests that enhanced DNA binding by SAP1a is not sufficient for recruitment by Pax-5, but instead involves protein–protein interactions mediated by the acidic side chain. Recruitment of Ets proteins by Pax-5 requires Gln22 within the N-terminal β-hairpin motif of its paired domain. The β-hairpin also participates in recognition of a subset of Pax-5-binding sites. Thus, Pax-5 incorporates protein–protein interaction and DNA recognition functions in a single motif. The Caenorhabditis elegans Pax protein EGL-38 also binds specifically to the mb-1 promoter and recruits murine Ets-1 or the C.elegans Ets protein T08H4.3, but not the related LIN-1 protein. Together, our results define specific amino acid requirements for Pax–Ets ternary complex assembly and show that the mechanism is conserved between evolutionarily related proteins of diverse animal species. Moreover, the data suggest that interactions between Pax and Ets proteins are an important mechanism that regulates fundamental biological processes in worms and humans.     

Cooperative binding of DNA and CBFbeta to the Runt domain of the CBFalpha studied via MD simulations

The Runt domain (RD) is the DNA-binding region of the Runx genes. A related protein, known as core binding factor β (CBFβ) also binds to the RD to enhance RD–DNA interaction by 6- to 10-fold. Here, we report results from molecular dynamics (MD) simulations of RD alone, as a dimer in complexes with DNA and CBFβ and in a ternary complex with DNA and CBFβ. Consistent with the experimental findings, in the presence of CBFβ the estimated free energy of binding of RD to the DNA is more favorable, which is shown to be due to more favorable intermolecular interactions and desolvation contributions. Also contributing to the enhanced binding are favorable intramolecular interactions between the ‘wing’ residues (RD residues 139–145) and the ‘wing1’ residues (RD residues 104–116). The simulation studies also indicate that the RD–CBFβ binding is more favorable in the presence of DNA due to a more favorable RD–CBFβ interaction energy. In addition, it is predicted that long-range interactions involving ionic residues contribute to binding cooperativity. Results from the MD calculations are used to interpret a variety of experimental mutagenesis data. A novel role for RD Glu116 to the RD–CBFβ interaction is predicted.     

NAR Breakthrough Article: A systematic survey of the Cys2His2 zinc finger DNA-binding landscape

Cys₂His₂ zinc fingers (C2H2-ZFs) comprise the largest class of metazoan DNA-binding domains. Despite this domain''s well-defined DNA-recognition interface, and its successful use in the design of chimeric proteins capable of targeting genomic regions of interest, much remains unknown about its DNA-binding landscape. To help bridge this gap in fundamental knowledge and to provide a resource for design-oriented applications, we screened large synthetic protein libraries to select binding C2H2-ZF domains for each possible three base pair target. The resulting data consist of >160 000 unique domain–DNA interactions and comprise the most comprehensive investigation of C2H2-ZF DNA-binding interactions to date. An integrated analysis of these independent screens yielded DNA-binding profiles for tens of thousands of domains and led to the successful design and prediction of C2H2-ZF DNA-binding specificities. Computational analyses uncovered important aspects of C2H2-ZF domain–DNA interactions, including the roles of within-finger context and domain position on base recognition. We observed the existence of numerous distinct binding strategies for each possible three base pair target and an apparent balance between affinity and specificity of binding. In sum, our comprehensive data help elucidate the complex binding landscape of C2H2-ZF domains and provide a foundation for efforts to determine, predict and engineer their DNA-binding specificities.     

Protein–DNA footprinting by endcapped duplex oligodeoxyribonucleotides

Oligodeoxyribonucleotides (5′-phosphorylated) of varying lengths were capped using a polyamide linker to form thermodynamically stable, endcapped DNA duplexes containing 8–14 bp. We have employed these endcapped DNA duplexes as tools to determine the DNA footprint of T4 DNA ligase. By high-performance liquid chromatography and PAGE analysis of the ligation mixtures of the endcapped DNA duplexes, we have found that by varying the lengths and the position of the nick, we can determine the minimal DNA-binding site as well as the mode of binding (symmetrical or asymmetrical binding) by the enzyme. The results of the study revealed that a 11 bp endcapped duplex was the shortest duplex effectively ligated. Dependence of ligation efficiency on nick position demonstrates that T4 DNA ligase bound asymmetrically to its DNA substrate. The use of a set of thermodynamically stable endcapped deoxyribonucleoside duplexes as a tool to elucidate the DNA footprint provides an efficient strategy for footprinting, which avoids ambiguities associated with chemical and biochemical footprinting methods.     

Mechanistic studies on DNA damage by minor groove binding copper-phenanthroline conjugates

Copper–phenanthroline complexes oxidatively damage and cleave nucleic acids. Copper bis-phenanthroline and copper complexes of mono- and bis-phenanthroline conjugates are used as research tools for studying nucleic acid structure and binding interactions. The mechanism of DNA oxidation and cleavage by these complexes was examined using two copper–phenanthroline conjugates of the sequence-specific binding molecule, distamycin. The complexes contained either one or two phenanthroline units that were bonded to the DNA-binding domain through a linker via the 3-position of the copper ligand. A duplex containing independently generated 2-deoxyribonolactone facilitated kinetic analysis of DNA cleavage. Oxidation rate constants were highly dependent upon the ligand environment but rate constants describing elimination of the alkali-labile 2-deoxyribonolactone intermediate were not. Rate constants describing DNA cleavage induced by each molecule were 11–54 times larger than the respective oxidation rate constants. The experiments indicate that DNA cleavage resulting from β-elimination of 2-deoxyribonolactone by copper–phenanthroline complexes is a general mechanism utilized by this family of molecules. In addition, the experiments confirm that DNA damage mediated by mono- and bis-phenanthroline copper complexes proceeds through distinct species, albeit with similar outcomes.     

Binding of hairpin pyrrole and imidazole polyamides to DNA: relationship between torsion angle and association rate constants

N-methylpyrrole (Py)-N-methylimidazole (Im) polyamides are small organic molecules that bind to DNA with sequence specificity and can be used as synthetic DNA-binding ligands. In this study, five hairpin eight-ring Py–Im polyamides 1–5 with different number of Im rings were synthesized, and their binding behaviour was investigated with surface plasmon resonance assay. It was found that association rate (k_a) of the Py–Im polyamides with their target DNA decreased with the number of Im in the Py–Im polyamides. The structures of four-ring Py–Im polyamides derived from density functional theory revealed that the dihedral angle of the Py amide carbonyl is 14∼18°, whereas that of the Im is significantly smaller. As the minor groove of DNA has a helical structure, planar Py–Im polyamides need to change their conformation to fit it upon binding to the minor groove. The data explain that an increase in planarity of Py–Im polyamide induced by the incorporation of Im reduces the association rate of Py–Im polyamides. This fundamental knowledge of the binding of Py–Im polyamides to DNA will facilitate the design of hairpin Py–Im polyamides as synthetic DNA-binding modules.     

Thermodynamic, kinetic and structural basis for recognition and repair of abasic sites in DNA by apurinic/apyrimidinic endonuclease from human placenta

X-ray analysis of enzyme–DNA interactions is very informative in revealing molecular contacts, but provides neither quantitative estimates of the relative importance of these contacts nor information on the relative contributions of specific and nonspecific interactions to the total affinity of enzymes for specific DNA. A stepwise increase in the ligand complexity approach is used to estimate the relative contributions of virtually every nucleotide unit of synthetic DNA containing abasic sites to its affinity for apurinic/apyrimidinic endonuclease (APE1) from human placenta. It was found that APE1 interacts with 9–10 nt units or base pairs of single-stranded and double-stranded ribooligonucleotides and deoxyribooligonucleotides of different lengths and sequences, mainly through weak additive contacts with internucleotide phosphate groups. Such nonspecific interactions of APE1 with nearly every nucleotide within its DNA-binding cleft provides up to seven orders of magnitude (ΔG° ~ −8.7 to −9.0 kcal/mol) of the enzyme affinity for any DNA substrate. In contrast, interactions with the abasic site together with other specific APE1–DNA interactions provide only one order of magnitude (ΔG° ~ −1.1 to −1.5 kcal/mol) of the total affinity of APE1 for specific DNA. We conclude that the enzyme's specificity for abasic sites in DNA is mostly due to a great increase (six to seven orders of magnitude) in the reaction rate with specific DNA, with formation of the Michaelis complex contributing to the substrate preference only marginally.     

Structure of the DNA-binding domain of NgTRF1 reveals unique features of plant telomere-binding proteins

Telomeres are protein–DNA elements that are located at the ends of linear eukaryotic chromosomes. In concert with various telomere-binding proteins, they play an essential role in genome stability. We determined the structure of the DNA-binding domain of NgTRF1, a double-stranded telomere-binding protein of tobacco, using multidimensional NMR spectroscopy and X-ray crystallography. The DNA-binding domain of NgTRF1 contained the Myb-like domain and C-terminal Myb-extension that is characteristic of plant double-stranded telomere-binding proteins. It encompassed amino acids 561–681 (NgTRF1^561–681), and was composed of 4 α-helices. We also determined the structure of NgTRF1^561–681 bound to plant telomeric DNA. We identified several amino acid residues that interacted directly with DNA, and confirmed their role in the binding of NgTRF1 to telomere using site-directed mutagenesis. Based on a structural comparison of the DNA-binding domains of NgTRF1 and human TRF1 (hTRF1), NgTRF1 has both common and unique DNA-binding properties. Interaction of Myb-like domain with telomeric sequences is almost identical in NgTRF1^561–681 with the DNA-binding domain of hTRF1. The interaction of Arg-638 with the telomeric DNA, which is unique in NgTRF1^561–681, may provide the structural explanation for the specificity of NgTRF1 to the plant telomere sequences, (TTTAGGG)_n.     

Delineation of the high-affinity single-stranded telomeric DNA-binding domain of Saccharomyces cerevisiae Cdc13

Cdc13 is an essential protein from Saccharomyces cerevisiae that caps telomeres by protecting the C-rich telomeric DNA strand from degradation and facilitates telomeric DNA replication by telomerase. In vitro, Cdc13 binds TG-rich single-stranded telomeric DNA with high affinity and specificity. A previously identified domain of Cdc13 encompassing amino acids 451–694 (the 451–694 DBD) retains the single-stranded DNA-binding properties of the full-length protein; however, this domain contains a large unfolded region identified in heteronuclear NMR experiments. Trypsin digestion and MALDI mass spectrometry were used to identify the minimal DNA-binding domain (the 497–694 DBD) necessary and sufficient for full DNA-binding activity. This domain was completely folded, and the N-terminal unfolded region removed was shown to be dispensable for function. Using affinity photocrosslinking to site-specifically modified telomeric single-stranded DNA, the 497–694 DBD was shown to contact the entire 11mer required for high-affinity binding. Intriguingly, both domains bound single-stranded telomeric DNA with much greater affinity than the full-length protein. The full-length protein exhibited the same rate of dissociation as both domains, however, indicating that the full-length protein contains a region that inhibits association with single-stranded telomeric DNA.     

Amino acid–base interactions: a three-dimensional analysis of protein–DNA interactions at an atomic level

To assess whether there are universal rules that govern amino acid–base recognition, we investigate hydrogen bonds, van der Waals contacts and water-mediated bonds in 129 protein–DNA complex structures. DNA–backbone interactions are the most numerous, providing stability rather than specificity. For base interactions, there are significant base–amino acid type correlations, which can be rationalised by considering the stereochemistry of protein side chains and the base edges exposed in the DNA structure. Nearly two-thirds of the direct read-out of DNA sequences involves complex networks of hydrogen bonds, which enhance specificity. Two-thirds of all protein–DNA interactions comprise van der Waals contacts, compared to about one-sixth each of hydrogen and water-mediated bonds. This highlights the central importance of these contacts for complex formation, which have previously been relegated to a secondary role. Although common, water-mediated bonds are usually non-specific, acting as space-fillers at the protein–DNA interface. In conclusion, the majority of amino acid–base interactions observed follow general principles that apply across all protein–DNA complexes, although there are individual exceptions. Therefore, we distinguish between interactions whose specificities are ‘universal’ and ‘context-dependent’. An interactive Web-based atlas of side chain–base contacts provides access to the collected data, including analyses and visualisation of the three-dimensional geometry of the interactions.