High-affinity DNA Research Articles

In silico post-SELEX screening and experimental characterizations for acquisition of high affinity DNA aptamers against carcinoembryonic antigen.

DNA aptamers against carcinoembryonic antigen (CEA) have been identified through the systematic evolution of ligands by exponential enrichment (SELEX) technique, but their affinity needs to be improved. In this study, an in silico approach was firstly used to screen the mutation sequences of a reported DNA aptamer (the parent aptamer, denoted as P) against CEA. The affinities of several high-score DNA mutants were determined by the biolayer interferometry technique. Finally, the newly obtained aptamers were verified in an aptasensor application. For the in silico approach, Mfold and RNA Composer were combined to generate the 3D RNA structures of the DNA mutants. The RNA structures were then modified to 3D DNA structures with the Write program. The docking model and binding ability of the 3D DNA structures with CEA were simulated and predicted with the ZDOCK program. Two mutation sequences (P-ATG and GAC-P) exhibited significantly higher ZDOCK scores than P. The dissociation constant of P-ATG and GAC-P to CEA was determined to be 4.62 and 3.93 nM respectively, obviously superior to that of P (6.95 nM). The detection limit of the P-ATG and GAC-P based aptasensors was 1.5 and 1.2 ng mL−1, respectively, markedly better than that based on P (3.4 ng mL−1). The consistency between the in silico and the experimental results indicates that the developed in silico post-SELEX screening approach is feasible for improving DNA aptamers. The P-ATG and GAC-P aptamers found in this study could be used for future CEA aptasensor design and fabrication, promisingly applicable for highly sensitive CEA detection and early cancer diagnosis.

Genetic Alphabet Expansion Provides Versatile Specificities and Activities of Unnatural-Base DNA Aptamers Targeting Cancer Cells

The potential of genetic alphabet expansion technologies using artificial extra base pairs (unnatural base pairs) has been rapidly expanding and increasing. We present that the hydrophobic unnatural base, 7-(2-thienyl)imidazo[4,5-b]pyridine (Ds), which acts as a fifth letter in a DNA library, provides a series of high-affinity DNA aptamers with versatile binding specificities and activities to cancer cells. These Ds-containing DNA aptamers were generated by a method called cell-ExSELEX to target three breast cancer cell lines: MCF7, MDA-MB-231, and T-47D. Aptamer 14A-MCF7, which targets MCF7 cells, specifically binds to MCF7 cells, but not other cancer cell lines. Aptamer 07-MB231, which targets MDA-MB-231 cells, binds to a series of metastatic bone and lung cancer cell lines. Aptamer 05-MB231 targets MDA-MB-231 cells, but it also binds to all of the cancer and leukemia cell lines that we examined. None of these aptamers bind to normal cell lines, such as MCF10A and HUVEC. In addition, aptamers 14A-MCF7 and 05-MB231 are internalized within the cancer cells, and aptamer 05-MB231 possesses anti-proliferative properties against most cancer cell lines that we examined. These aptamers and the generation method are broadly applicable to cancer cell imaging, biomarker discovery, cancer cell profiling, anti-cancer therapies, and drug delivery systems.

Computational studies on G-quadruplex DNA-stabilizing property of novel Wittig-based Schiff-Base ligands and their copper(II) complexes

A series of mono imine (C=N) group that contained Wittig-based Schiff-Base ligands was optimized using the DFT-based computational method and Gaussian 09 program package. These optimized molecules were docked with Quadruplex DNA (PDB ID: 1KF1) and duplex DNA (PDB ID: 1BNA) using AutoDock Vina program along with the reference molecules. Schiff-Base ligands derived from fused aromatic rings contained amines and precursor aldehyde (PA-5 both Z and E isomers) showed lower binding energy for G-quadruplex DNA among all and N-5 category both Z and E isomer Schiff-Base ligands have shown high selectivity for G-quadruplex DNA over duplex DNA which is a very important phenomenon to develop the G-quadruplex DNA stabilizers. For a few Schiff-Base molecules like Ligand-6 (1-{[2-Hydroxy-5-(2-pyren-1-yl-vinyl)-benzylidene]-amino}-naphthalen-2-ol) of N-5 category both Z and E isomers with groove binding and end stacking modes, molecular dynamic calculations were carried out. The study revealed that Ligand-6 of N-5 category E isomer with groove binding mode has higher stabilizing effect on G-quadruplex DNA in spite of having the higher binding energy value. Among Schiff-Base copper(II) complexes, Complex-3 (E-(1-{[2-Hydroxy-5-(2-pyren-1-yl-vinyl)-benzylidene]-amino}-naphthalen-2-ol)Cu) is having high binding affinity for G-quadruplex DNA as compared to others.

Novel Amphiphilic G-Quadruplex Binding Synthetic Derivative of TMPyP4 and Its Effect on Cancer Cell Proliferation and Apoptosis Induction.

Porphyrins are well-known anticancer agents because of their high binding affinity for G-quadruplex DNA and excellent photophysical properties. Several studies carried out using TMPyP4 established it as an efficient chemotherapeutic and a photodynamic therapeutic (PDT) agent, but its use as a lead molecule has been restricted because of its high level of binding to double-stranded DNA (dsDNA), which may have side effects on normal cells and tissues. To minimize its interaction with dsDNA and to enhance internalization into cells, an analogue of TMPyP4 (5Me) was synthesized. Its selectivity for G-quadruplex DNA over dsDNA was evaluated by spectroscopic methods, and its role in stabilizing G-quadruplex DNA was assessed by fluorescence lifetime and thermal melting experiments. Biophysical studies indicated that 5Me interacts well with G-quadruplex DNA. In vitro cytotoxicity experiments with tumor cell lines (PANC-1, B16F10, and MDA MB 231) have revealed that 5Me can inhibit the growth of cancer cells comparable to TMPyP4. MTT and apoptotic assays demonstrated the ability of 5Me to specifically affect cancer cells over normal cells. Cell cycle analysis showed that 5Me, like TMPyP4, induces G2/M phase cell cycle arrest. In addition, 5Me is more effectively taken up by both cancer and normal cells than TMPyP4. In addition, we have noticed that 5Me is more efficient than TMPyP4 in inhibiting the growth of the cancer cells after irradiation with light (600-720 nm, 20 J/cm2, 50 mW/cm2). By and large, these experimental results indicate that 5Me can be an efficient chemotherapeutic as well as a PDT agent.

Specificity landscapes unmask submaximal binding site preferences of transcription factors

We have developed Differential Specificity and Energy Landscape (DiSEL) analysis to comprehensively compare DNA-protein interactomes (DPIs) obtained by high-throughput experimental platforms and cutting edge computational methods. While high-affinity DNA binding sites are identified by most methods, DiSEL uncovered nuanced sequence preferences displayed by homologous transcription factors. Pairwise analysis of 726 DPIs uncovered homolog-specific differences at moderate- to low-affinity binding sites (submaximal sites). DiSEL analysis of variants of 41 transcription factors revealed that many disease-causing mutations result in allele-specific changes in binding site preferences. We focused on a set of highly homologous factors that have different biological roles but "read" DNA using identical amino acid side chains. Rather than direct readout, our results indicate that DNA noncontacting side chains allosterically contribute to sculpt distinct sequence preferences among closely related members of transcription factor families.

The ATPase activity of E. coli RecA prevents accumulation of toxic complexes formed by erroneous binding to undamaged double stranded DNA.

The Escherichia coli RecA protein catalyzes the central step of homologous recombination using its homology search and strand exchange activity. RecA is a DNA-dependent ATPase, but its homology search and strand exchange activities are largely independent of its ATPase activity. ATP hydrolysis converts a high affinity DNA binding form, RecA-ATP, to a low affinity form RecA-ADP, thereby supporting an ATP hydrolysis-dependent dynamic cycle of DNA binding and dissociation. We provide evidence for a novel function of RecA’s dynamic behavior; RecA’s ATPase activity prevents accumulation of toxic complexes caused by direct binding of RecA to undamaged regions of dsDNA. We show that a mutant form of RecA, RecA-K250N, previously shown to be toxic to E. coli, is a loss-of-function ATPase-defective mutant. We use a new method for detecting RecA complexes involving nucleoid surface spreading and immunostaining. The method allows detection of damage-induced RecA foci; STED microscopy revealed these to typically be between 50 and 200 nm in length. RecA-K250N, and other toxic variants of RecA, form spontaneous DNA-bound complexes that are independent of replication and of accessory proteins required to load RecA onto tracts of ssDNA in vivo, supporting the hypothesis that RecA’s expenditure of ATP serves an error correction function.

Non-SELEX isolation of DNA aptamers for the homogeneous-phase fluorescence anisotropy sensing of tau Proteins

Herein, we report for the first time the isolation of DNA aptamers directed against the whole tau protein, an important Alzheimer's disease (AD) biomarker. Non-SELEX approach based on the capillary electrophoresis partitioning technique was employed to isolate a high-affinity DNA sequence pool towards the target in only three rounds and one working day. High-throughput sequencing was next performed and the recognition ability of five selected aptamers was preliminary evaluated by surface plasmon resonance using the protein target immobilized on the chip. Finally, the analytical potential of the most affine aptamer was demonstrated through the design of a homogeneous-phase fluorescence anisotropy assay. This DNA aptamer was found to be able to recognize not only the whole τ-441 but also the τ-381, τ-352, τ-383 isoforms. The sensing platform allowed the determination of these four targets with a detection limit of 28 nM, 3.2 nM, 6.3 nM and 22 nM, respectively.

Generation and application of DNA aptamers against HspX for accurate diagnosis of tuberculous meningitis

Tuberculous meningitis (TBM) is the most severe manifestation of tuberculosis and its diagnosis remains a challenge even today due to the lack of an adequate test. HspX antigen of Mycobacterium tuberculosis was previously established as a reliable diagnostic biomarker for TBM in an ELISA test format using anti-HspX polyclonal antibodies. Towards overcoming the limitations of batch-to-batch variation and challenges of scalability in antibody generation, we utilized Systematic Evolution of Ligands by EXponential enrichment (SELEX) to develop high affinity DNA aptamers against HspX as an alternative diagnostic reagent. Post-SELEX optimization of the best-performing aptamer candidate, H63, established its derivative H63 SL-2 M6 to be superior to its parent. Aptamer H63 SL-2 M6 displayed a specific and high affinity interaction with HspX (Kd ∼9.0 × 10−8 M). In an Aptamer Linked Immobilized Sorbent Assay (ALISA), H63 SL-2 M6 significantly differentiated between cerebrospinal fluid specimens from TBM and non-TBM subjects (n = 87, ***p < 0.0001) with ∼100% sensitivity and ∼91% specificity. Notably, ALISA exhibited comparable performance with previously reported antibody-based ELISA and qPCR. Altogether, our findings establish the utility of HspX aptamer for the reliable diagnosis of TBM and pave the way for developing an aptamer-based point-of-care test for TBM.

Novel polypyridyl ruthenium complexes acting as high affinity DNA intercalators, potent transcription inhibitors and antitumor reagents

Six novel polypyridyl ruthenium complexes with (E)-2-styryl-1H- imidazo[4,5-f][1,10]phenanthroline ligand and its analogues have been designed to enhance the DNA intercalation ability of their model compound [Ru(bpy)2(pip)]2+ (bpy = 2,2′-bipyridine, pip = 2-phenyl-1H-imidazo[4,5-f][1,10]phenanthroline). As shown in the optimized geometry of the complexes, the introduction of styryl group not only extended the conjugated area of the intercalative ligand, but also retained the excellent planarity. These two merits have been proven to be beneficial for their DNA intercalation, thus greatly improved their inhibition activity towards DNA transcription by RNA polymerase and DNA topoisomerase, two enzymes closely related to both DNA and tumor cell growth. The relationships between the substituent group structures and the biological activities have also been investigated from energetic and electronic aspects by quantum chemistry calculations. Results from cell cytotoxicity and apoptosis assay testified that the styryl substituted ruthenium complexes possessed higher antitumor activity than [Ru(bpy)2(pip)]2+, as expected. As quantified in the MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide) assay, the tumor cell death is caused mostly through apoptosis for Ru2 and Ru3, while non-apoptotic processes for Ru1, Ru4 and Ru5. In vitro fluorescence evaluation revealed that all complexes located mainly in cytoplasm, but the three complexes with high antiproliferative activity could enter nucleus. All complexes have shown apparent lower cytotoxicity towards normal human colon epithelial cell CCD-841-CON than the examined tumor cell lines.

Protein Cofactors Are Essential for High-Affinity DNA Binding by the Nuclear Factor κB RelA Subunit.

Transcription activator proteins typically contain two functional domains: a DNA binding domain (DBD) that binds to DNA with sequence specificity and an activation domain (AD) whose established function is to recruit RNA polymerase. In this report, we show that purified recombinant nuclear factor κB (NF-κB) RelA dimers bind specific κB DNA sites with an affinity significantly lower than that of the same dimers from nuclear extracts of activated cells, suggesting that additional nuclear cofactors might facilitate DNA binding by the RelA dimers. Additionally, recombinant RelA binds DNA with relatively low affinity at a physiological salt concentration in vitro. The addition of p53 or RPS3 (ribosomal protein S3) increases RelA:DNA binding affinity 2- to >50-fold depending on the protein and ionic conditions. These cofactor proteins do not form stable ternary complexes, suggesting that they stabilize the RelA:DNA complex through dynamic interactions. Surprisingly, the RelA-DBD alone fails to bind DNA under the same solution conditions even in the presence of cofactors, suggesting an important role of the RelA-AD in DNA binding. Reduced RelA:DNA binding at a physiological ionic strength suggests that multiple cofactors might be acting simultaneously to mitigate the electrolyte effect and stabilize the RelA:DNA complex in vivo. Overall, our observations suggest that the RelA-AD and multiple cofactor proteins function cooperatively to prime the RelA-DBD and stabilize the RelA:DNA complex in cells. Our study provides a mechanism for nuclear cofactor proteins in NF-κB-dependent gene regulation.

A highly sensitive aptamer-based HIV reverse transcriptase detection assay

Although many new assays for HIV have been developed, several labs still use simple and reliable radioactivity-based reverse transcriptase (RT) nucleotide incorporation assays for detection and quantification. We describe here a new assay for detection and quantitation of HIV RT activity that is based on a high affinity DNA aptamer to RT. The aptamer is sequestered on 96-well plates where it can bind to RT and other constituents can be removed by extensive washing. Since the aptamer mimics a primer-template, upon radiolabeled nucleotide addition, bound RT molecules can extend the aptamer and the radioactive signal can be detected by standard methods. In addition to being procedurally simple, the assay demonstrated high sensitivity (detection limits for RT and virions were ≤6400 molecules (∼4 × 10−8 units) and ∼100–300 virions, respectively) and was essentially linear over a range of at least 104. Both wild type and drug-resistant forms of HIV-1 RT were detectable as was HIV-2 RT, although there were some modest differences in sensitivity.

Determining if DNA Stained with a Cyanine Dye Can Be Digested with Restriction Enzymes

Visualization of DNA for fluorescence microscopy utilizes a variety of dyes such as cyanine dyes. These dyes are utilized due to their high affinity and sensitivity for DNA. In order to determine if the DNA molecules are full length after the completion of the experiment, a method is required to determine if the stained molecules are full length by digesting DNA with restriction enzymes. However, stained DNA may inhibit the enzymes, so a method is needed to determine what enzymes one could use for fluorochrome stained DNA. In this method, DNA is stained with a cyanine dye overnight to allow the dye and DNA to equilibrate. Next, stained DNA is digested with a restriction enzyme, loaded into a gel and electrophoresed. The experimental DNA digest bands are compared to an in silico digest to determine the restriction enzyme activity. If there is the same number of bands as expected, then the reaction is complete. More bands than expected indicate partial digestion and less bands indicate incomplete digestion. The advantage of this method is its simplicity and it uses equipment that a scientist would need for a restriction enzyme assay and gel electrophoresis. A limitation of this method is that the enzymes available to most scientists are commercially available enzymes; however, any restriction enzymes could be used.

Determining if DNA Stained with a Cyanine Dye Can Be Digested with Restriction Enzymes.

Visualization of DNA for fluorescence microscopy utilizes a variety of dyes such as cyanine dyes. These dyes are utilized due to their high affinity and sensitivity for DNA. In order to determine if the DNA molecules are full length after the completion of the experiment, a method is required to determine if the stained molecules are full length by digesting DNA with restriction enzymes. However, stained DNA may inhibit the enzymes, so a method is needed to determine what enzymes one could use for fluorochrome stained DNA. In this method, DNA is stained with a cyanine dye overnight to allow the dye and DNA to equilibrate. Next, stained DNA is digested with a restriction enzyme, loaded into a gel and electrophoresed. The experimental DNA digest bands are compared to an in silico digest to determine the restriction enzyme activity. If there is the same number of bands as expected, then the reaction is complete. More bands than expected indicate partial digestion and less bands indicate incomplete digestion. The advantage of this method is its simplicity and it uses equipment that a scientist would need for a restriction enzyme assay and gel electrophoresis. A limitation of this method is that the enzymes available to most scientists are commercially available enzymes; however, any restriction enzymes could be used.

The WYL Domain of the PIF1 Helicase from the Thermophilic Bacterium Thermotoga elfii is an Accessory Single-Stranded DNA Binding Module.

PIF1 family helicases are conserved from bacteria to man. With the exception of the well-studied yeast PIF1 helicases (e.g., ScPif1 and ScRrm3), however, very little is known about how these enzymes help maintain genome stability. Indeed, we lack a basic understanding of the protein domains found N- and C-terminal to the characteristic central PIF1 helicase domain in these proteins. Here, using chimeric constructs, we show that the ScPif1 and ScRrm3 helicase domains are interchangeable and that the N-terminus of ScRrm3 is important for its function in vivo. This suggests that PIF1 family helicases evolved functional modules fused to a generic motor domain. To investigate this hypothesis, we characterized the biochemical activities of the PIF1 helicase from the thermophilic bacterium Thermotoga elfii (TePif1), which contains a C-terminal WYL domain of unknown function. Like helicases from other thermophiles, recombinant TePif1 was easily prepared, thermostable in vitro, and displayed activities similar to its eukaryotic homologues. We also found that the WYL domain was necessary for high-affinity single-stranded DNA (ssDNA) binding and affected both ATPase and helicase activities. Deleting the WYL domain from TePif1 or mutating conserved residues in the predicted ssDNA binding site uncoupled ATPase activity and DNA unwinding, leading to higher rates of ATP hydrolysis but less efficient DNA helicase activity. Our findings suggest that the domains of unknown function found in eukaryotic PIF1 helicases may also confer functional specificity and additional activities to these enzymes, which should be investigated in future work.

BET-seq: Binding energy topographies revealed by microfluidics and high-throughput sequencing.

Biophysical models of transcriptional regulation rely on energetic measurements of the binding affinities between transcription factors (TFs) and target DNA binding sites. Historically, assays capable of measuring TF-DNA binding affinities have been relatively low-throughput (measuring ~103 sequences in parallel) and have required significant specialized equipment, limiting their use to a handful of laboratories. Recently, we developed an experimental assay and analysis pipeline that allows measurement of binding energies between a single TF and up to 106 DNA species in a single experiment (Binding Energy Topography by sequencing, or BET-seq) (Le et al., 2018). BET-seq employs the Mechanically Induced Trapping of Molecular Interactions (MITOMI) platform to purify DNA bound to a TF at equilibrium followed by high coverage sequencing to reveal relative differences in binding energy for each sequence. While we have previously used BET-seq to refine the binding affinity landscapes surrounding high-affinity DNA consensus target sites, we anticipate this technique will be applied in future work toward measuring a wide variety of TF-DNA landscapes. Here, we provide detailed instructions and general considerations for DNA library design, performing BET-seq assays, and analyzing the resulting data.

The C-Terminal Zinc Fingers of ZBTB38 are Novel Selective Readers of DNA Methylation

Methyl-CpG binding proteins play an essential role in translating DNA methylation marks into a downstream transcriptional response, which has implications for both normal cell function as well as disease. Although for many of these proteins, a detailed mechanistic understanding for how this cellular process is mediated remains to be determined. ZBTB38 is an under-characterized member of the zinc finger (ZF) family of methyl-CpG binding proteins. Functional knowledge has been gained for its conserved methylated DNA binding N-terminal ZF region; however, a specific role for the C-terminal set of five ZFs remains to be elucidated. Here we demonstrate for the first time that a subset of the C-terminal ZBTB38 ZFs exhibit high-affinity DNA interactions and that preferential targeting of the consensus DNA site is methyl specific. Utilizing a hybrid approach, a model for the C-terminal ZBTB38 ZFs in complex with its cognate DNA target is proposed, providing insight into a possible novel mode of methylated DNA recognition. Furthermore, it is shown that the C-terminal ZFs of ZBTB38 can directly occupy promoters harboring the newly identified sequence motif in cell in a methyl-dependent manner and, depending on the gene context, contribute to modulating transcriptional response. Combined, these findings provide evidence for a key and novel physiological function for the C-terminal ZF domain of ZBTB38.

A peculiar IclR family transcription factor regulates para-hydroxybenzoate catabolism in Streptomyces coelicolor.

In Streptomyces coelicolor, we identified a para-hydroxybenzoate (PHB) hydroxylase, encoded by gene pobA (SCO3084), which is responsible for conversion of PHB into PCA (protocatechuic acid), a substrate of the β-ketoadipate pathway which yields intermediates of the Krebs cycle. We also found that the transcription of pobA is induced by PHB and is negatively regulated by the product of SCO3209, which we named PobR. The product of this gene is highly unusual in that it is the apparent fusion of two IclR family transcription factors. Bioinformatic analyses, in vivo transcriptional assays, electrophoretic mobility shift assays (EMSAs), DNase I footprinting, and isothermal calorimetry (ITC) were used to elucidate the regulatory mechanism of PobR. We found that PobR loses its high affinity for DNA (i.e., the pobA operator) in the presence of PHB, the inducer of pobA transcription. PHB binds to PobR with a KD of 5.8 μM. Size-exclusion chromatography revealed that PobR is a dimer in the absence of PHB and a monomer in the presence of PHB. The crystal structure of PobR in complex with PHB showed that only one of the two IclR ligand binding domains was occupied, and defined how the N-terminal ligand binding domain engages the effector ligand.

Amiloride-enhanced gene transfection of octa-arginine functionalized calcium phosphate nanoparticles.

Nanoparticles represent promising gene delivery systems in biomedicine to facilitate prolonged gene expression with low toxicity compared to viral vectors. Specifically, nanoparticles of calcium phosphate (nCaP), the main inorganic component of human bone, exhibit high biocompatibility and good biodegradability and have been reported to have high affinity for protein or DNA, having thus been used as gene transfer vectors. On the other hand, Octa-arginine (R8), which has a high permeability to cell membrane, has been reported to improve intracellular delivery systems. Here, we present an optimized method for nCaP-mediated gene delivery using an octa-arginine (R8)-functionalized nCaP vector containing a marker or functional gene construct. nCaP particle size was between 220–580 nm in diameter and all R8-functionalized nCaPs carried a positive charge. R8 concentration significantly improved nCaP transfection efficiency with high cell compatibility in human mesenchymal stem cells (hMSC) and human osteoblasts (hOB) in particular, suggesting nCaPs as a good option for non-viral vector gene delivery. Furthermore, pre-treatment with different endocytosis inhibitors identified that the endocytic pathway differed among cell lines and functionalized nanoparticles, with amiloride increasing transfection efficiency of R8-functionalized nCaPs in hMSC and hOB.

Site-specific aptamer inhibitors of Thermus RNA polymerase

Bacterial RNA polymerase (RNAP) is an RNA-synthesizing molecular machine and a target for antibiotics. In transcription, RNAP can interact with DNA sequence-specifically, during promoter recognition by the σ-containing holoenzyme, or nonspecifically, during productive RNA elongation by the core RNAP. We describe high-affinity single-stranded DNA aptamers that are specifically recognized by the core RNAP from Thermus aquaticus. The aptamers interact with distinct epitopes inside the RNAP main channel, including the rifamycin pocket, and sense the binding of other RNAP ligands such as rifamycin and the σA subunit. The aptamers inhibit RNAP activity and can thus be used for functional studies of transcription and development of novel RNAP inhibitors.

Coop-Seq Analysis Demonstrates that Sox2 Evokes Latent Specificities in the DNA Recognition by Pax6

Sox2 and Pax6 co-regulate genes in neural lineages and the lens by forming a ternary complex likely facilitated allosterically through DNA. We used the quantitative and scalable cooperativity-by-sequencing (Coop-seq) approach to interrogate Sox2/Pax6 dimerization on a DNA library where five positions of the Pax6 half-site were randomized yielding 1024 cooperativity factors. Consensus positions normally required for the high-affinity DNA binding by Pax6 need to be mutated for effective dimerization with Sox2. Out of the five randomized bases, a 5′ thymidine is present in most of the top ranking elements. However, this thymidine maps to a region outside of the Pax half site and is not expected to directly interact with Pax6 in known binding modes suggesting structural reconfigurations. Re-analysis of ChIP-seq data identified several genomic regions where the cooperativity promoting sequence pattern is co-bound by Sox2 and Pax6. A highly conserved Sox2/Pax6 bound site near the Sprouty2 locus was verified to promote cooperative dimerization designating Sprouty2 as a potential target reliant on Sox2/Pax6 cooperativity in several neural cell types. Collectively, the functional interplay of Sox2 and Pax6 demands the relaxation of high-affinity binding sites and is enabled by alternative DNA sequences. We conclude that this binding mode evolved to warrant that a subset of target genes is only regulated in the presence of suitable partner factors.