Dissection of integrated readout reveals the structural thermodynamics of DNA selection by transcription factors.

Nucleobases such as inosine have been extensively utilized to map direct contacts by proteins in the DNA groove. Their deployment as targeted probes of dynamics and hydration, which are dominant thermodynamic drivers of affinity and specificity, has been limited by a paucity of suitable experimental models. We report a joint crystallographic, thermodynamic, and computational study of the bidentate complex of the arginine side chain with a Watson-Crick guanine (Arg×GC), a highly specific configuration adopted by major transcription factors throughout the eukaryotic branches in the Tree of Life. Using the ETS-family factor PU.1 as a high-resolution structural framework, inosine substitution for guanine resulted in a sharp dissection of conformational dynamics and hydration and elucidated their role in the DNA specificity of PU.1. Our work suggests an under-exploited utility of modified nucleobases in untangling the structural thermodynamics of interactions, such as the Arg×GC motif, where direct and indirect readout are tightly integrated.

DNA selection by the master transcription factor PU.1.

The master transcriptional regulator PU.1/Spi-1 engages DNA sites with affinities spanning multiple orders of magnitude. To elucidate this remarkable plasticity, we have characterized 22 high-resolution co-crystallographic PU.1/DNA complexes across the addressable affinity range in myeloid gene transactivation. Over a purine-rich core (such as 5'-GGAA-3') flanked by variable sequences, affinity is negotiated by direct readout on the 5' flank via a critical glutamine (Q226) sidechain and by indirect readout on the 3' flank by sequence-dependent helical flexibility. Direct readout by Q226 dynamically specifies PU.1's characteristic preference for purines and explains the pathogenic mutation Q226E in Waldenström macroglobulinemia. The structures also reveal how disruption of Q226 mediates strand-specific inhibition by DNA methylation and the recognition of non-canonical sites, including the authentic binding sequence at the CD11b promoter. A re-synthesis of phylogenetic and structural data on the ETS family, considering the centrality of Q226 in PU.1, unifies the model of DNA selection by ETS proteins.

Dissecting knowledge, guessing, and blunder in multiple choice assessments.

Multiple choice results are inherently probabilistic outcomes, as correct responses reflect a combination of knowledge and guessing, while incorrect responses additionally reflect blunder, a confidently committed mistake. To objectively resolve knowledge from responses in an MC test structure, we evaluated probabilistic models that explicitly account for guessing, knowledge and blunder using eight assessments (>9,000 responses) from an undergraduate biotechnology curriculum. A Bayesian implementation of the models, aimed at assessing their robustness to prior beliefs in examinee knowledge, showed that explicit estimators of knowledge are markedly sensitive to prior beliefs with scores as sole input. To overcome this limitation, we examined self-ranked confidence as a proxy knowledge indicator. For our test set, three levels of confidence resolved test performance. Responses rated as least confident were correct more frequently than expected from random selection, reflecting partial knowledge, but were balanced by blunder among the most confident responses. By translating evidence-based guessing and blunder rates to pass marks that statistically qualify a desired level of examinee knowledge, our approach finds practical utility in test analysis and design.

Detecting recurrent passenger mutations in melanoma by targeted UV damage sequencing.

Sequencing of melanomas has identified hundreds of recurrent mutations in both coding and non-coding DNA. These include a number of well-characterized oncogenic driver mutations, such as coding mutations in the BRAF and NRAS oncogenes, and non-coding mutations in the promoter of telomerase reverse transcriptase (TERT). However, the molecular etiology and significance of most of these mutations is unknown. Here, we use a new method known as CPD-capture-seq to map UV-induced cyclobutane pyrimidine dimers (CPDs) with high sequencing depth and single nucleotide resolution at sites of recurrent mutations in melanoma. Our data reveal that many previously identified drivers and other recurrent mutations in melanoma occur at CPD hotspots in UV-irradiated melanocytes, often associated with an overlapping binding site of an E26 transformation-specific (ETS) transcription factor. In contrast, recurrent mutations in the promoters of a number of known or suspected cancer genes are not associated with elevated CPD levels. Our data indicate that a subset of recurrent protein-coding mutations are also likely caused by ETS-induced CPD hotspots. This analysis indicates that ETS proteins profoundly shape the mutation landscape of melanoma and reveals a method for distinguishing potential driver mutations from passenger mutations whose recurrence is due to elevated UV damage.

Molecular mechanism of UV damage modulation in nucleosomes.

Exposure to ultraviolet (UV) light causes the formation of mutagenic cyclobutane pyrimidine dimers (CPDs) in cellular DNA. Previous studies have revealed that CPD formation in nucleosomes, the building blocks of chromatin, shows a striking ∼10 base pair (bp) periodic pattern. CPD formation is suppressed at positions where the DNA minor groove faces toward the histone octamer (minor-in) and elevated CPD formation at positions where the minor groove faces away from the histone octamer (minor-out). However, the molecular mechanism underlying this nucleosome photofootprint is unclear. Here, we analyzed ∼180 high-resolution nucleosome structures to characterize whether differences in DNA mobility or conformation are responsible for the CPD modulation in nucleosomes. Our results indicate that differences in DNA mobility cannot explain CPD modulation in nucleosome. Instead, we find that the sharp DNA bending around the histone octamer results in DNA conformations with structural parameters more susceptible to UV damage formation at minor-out positions and more resistant to CPD formation at minor-in positions. This analysis reveals the molecular mechanism responsible for periodic modulation of CPD formation and UV mutagenesis in nucleosomal DNA.

Drug design and DNA structural research inspired by the Neidle laboratory: DNA minor groove binding and transcription factor inhibition by thiophene diamidines.

The understanding of sequence-specific DNA minor groove interactions has recently made major steps forward and as a result, the goal of development of compounds that target the minor groove is an active research area. In an effort to develop biologically active minor groove agents, we are preparing and exploring the DNA interactions of diverse diamidine derivatives with a 5'-GAATTC-3' binding site using a powerful array of methods including, biosensor-SPR methods, and X-ray crystallography. The benzimidazole-thiophene module provides an excellent minor groove recognition component. A central thiophene in a benzimidazole-thiophene-phenyl aromatic system provides essentially optimum curvature for matching the shape of the minor groove. Comparison of that structure to one with the benzimidazole replaced with an indole shows that the two structures are very similar, but have some interesting and important differences in electrostatic potential maps, the DNA minor groove binding structure based on x-ray crystallographic analysis, and inhibition of the major groove binding PU.1 transcription factor complex. The binding KD for both compounds is under 10 nM and both form amidine H-bonds to DNA bases. They both have bifurcated H-bonds from the benzimidazole or indole groups to bases at the center of the -AATT- binding site. Analysis of the comparative results provides an excellent understanding of how thiophene compounds recognize the minor groove and can act as transcription factor inhibitors.

Self-Consistent Parameterization of DNA Residues for the Non-Polarizable AMBER Force Fields.

Fixed-charge (non-polarizable) forcefields are accurate and computationally efficient tools for modeling the molecular dynamics of nucleic acid polymers, particularly DNA, well into the µs timescale. The continued utility of these forcefields depends in part on expanding the residue set in step with advancing nucleic acid chemistry and biology. A key step in parameterizing new residues is charge derivation which is self-consistent with the existing residues. As atomic charges are derived by fitting against molecular electrostatic potentials, appropriate structural models are critical. Benchmarking against the existing charge set used in current AMBER nucleic acid forcefields, we report that quantum mechanical models of deoxynucleosides, even at a high level of theory, are not optimal structures for charge derivation. Instead, structures from molecular mechanics minimization yield charges with up to 6-fold lower RMS deviation from the published values, due to the choice of such an approach in the derivation of the original charge set. We present a contemporary protocol for rendering self-consistent charges as well as optimized charges for a panel of nine non-canonical residues that will permit comparison with literature as well as studying the dynamics of novel DNA polymers.

High-resolution mapping demonstrates inhibition of DNA excision repair by transcription factors.

DNA base damage arises frequently in living cells and needs to be removed by base excision repair (BER) to prevent mutagenesis and genome instability. Both the formation and repair of base damage occur in chromatin and are conceivably affected by DNA-binding proteins such as transcription factors (TFs). However, to what extent TF binding affects base damage distribution and BER in cells is unclear. Here, we used a genome-wide damage mapping method, N-methylpurine-sequencing (NMP-seq), and characterized alkylation damage distribution and BER at TF binding sites in yeast cells treated with the alkylating agent methyl methanesulfonate (MMS). Our data show that alkylation damage formation was mainly suppressed at the binding sites of yeast TFs ARS binding factor 1 (Abf1) and rDNA enhancer binding protein 1 (Reb1), but individual hotspots with elevated damage levels were also found. Additionally, Abf1 and Reb1 binding strongly inhibits BER in vivo and in vitro, causing slow repair both within the core motif and its adjacent DNA. Repair of ultraviolet (UV) damage by nucleotide excision repair (NER) was also inhibited by TF binding. Interestingly, TF binding inhibits a larger DNA region for NER relative to BER. The observed effects are caused by the TF-DNA interaction, because damage formation and BER can be restored by depletion of Abf1 or Reb1 protein from the nucleus. Thus, our data reveal that TF binding significantly modulates alkylation base damage formation and inhibits repair by the BER pathway. The interplay between base damage formation and BER may play an important role in affecting mutation frequency in gene regulatory regions.

The Non-continuum Nature of Eukaryotic Transcriptional Regulation.

Eukaryotic transcription factors are versatile mediators of specificity in gene regulation. This versatility is achieved through mutual specification by context-specific DNA binding on the one hand, and identity-specific protein-protein partnerships on the other. This interactivity, known as combinatorial control, enables a repertoire of complex transcriptional outputs that are qualitatively disjoint, or non-continuum, with respect to binding affinity. This feature contrasts starkly with prokaryotic gene regulators, whose activities in general vary quantitatively in step with binding affinity. Biophysical studies on prokaryotic model systems and more recent investigations on transcription factors highlight an important role for folded state dynamics and molecular hydration in protein/DNA recognition. Analysis of molecular models of combinatorial control and recent literature in low-affinity gene regulation suggest that transcription factors harbor unique conformational dynamics that are inaccessible or unused by prokaryotic DNA-binding proteins. Thus, understanding the intrinsic dynamics involved in DNA binding and co-regulator recruitment appears to be a key to understanding how transcription factors mediate non-continuum outcomes in eukaryotic gene expression, and how such capability might have evolved from ancient, structurally conserved counterparts.

CTCF binding modulates UV damage formation to promote mutation hot spots in melanoma.

Somatic mutations in DNA-binding sites for CCCTC-binding factor (CTCF) are significantly elevated in many cancers. Prior analysis has suggested that elevated mutation rates at CTCF-binding sites in skin cancers are a consequence of the CTCF-cohesin complex inhibiting repair of UV damage. Here, we show that CTCF binding modulates the formation of UV damage to induce mutation hot spots. Analysis of genome-wide CPD-seq data in UV-irradiated human cells indicates that formation of UV-induced cyclobutane pyrimidine dimers (CPDs) is primarily suppressed by CTCF binding but elevated at specific locations within the CTCF motif. Locations of CPD hot spots in the CTCF-binding motif coincide with mutation hot spots in melanoma. A similar pattern of damage formation is observed at CTCF-binding sites in vitro, indicating that UV damage modulation is a direct consequence of CTCF binding. We show that CTCF interacts with binding sites containing UV damage and inhibits repair by a model repair enzyme in vitro. Structural analysis and molecular dynamic simulations reveal the molecular mechanism for how CTCF binding modulates CPD formation.

DNA-facilitated target search by nucleoproteins: Extension of a biosensor-surface plasmon resonance method.

To extend the value of biosensor-SPR in the characterization of DNA recognition by nucleoproteins, we report a comparative analysis of DNA-facilitated target search by two ETS-family transcription factors: Elk1 and ETV6. ETS domains represent an attractive system for developing biosensor-based techniques due to a broad range of physicochemical properties encoded within a highly conserved DNA-binding motif. Building on a biosensor approach in which the protein is quantitatively sequestered and presented to immobilized cognate DNA as nonspecific complexes, we assessed the impact of intrinsic cognate and nonspecific affinities on long-range (intersegmental) target search. The equilibrium constants of DNA-facilitated binding were sensitive to the intrinsic binding properties of the proteins such that their relative specificity for cognate DNA were reinforced when binding occurred by transfer vs. without nonspecific DNA. Direct measurement of association and dissociation kinetics revealed ionic features of the activated complex that evidenced DNA-facilitated dissociation, even though Elk1 and ETV6 harbor only a single DNA-binding surface. At salt concentrations that masked the effects of nonspecific pre-binding at equilibrium, the dissociation kinetics of cognate binding were nevertheless distinct from conditions under which nonspecific DNA was absent. These results further strengthen the significance of long-range DNA-facilitated translocation in the physiologic environment.

Salt bridge dynamics in protein/DNA recognition: a comparative analysis of Elk1 and ETV6.

Electrostatic protein/DNA interactions arise from the neutralization of the DNA phosphodiester backbone as well as coupled exchanges by charged protein residues as salt bridges or with mobile ions. Much focus has been and continues to be paid to interfacial ion pairs with DNA. The role of extra-interfacial ionic interactions, particularly as dynamic drivers of DNA sequence selectivity, remain poorly known. The ETS family of transcription factors represents an attractive model for addressing this knowledge gap given their diverse ionic composition in primary structures that fold to a tightly conserved DNA-binding motif. To probe the importance of extra-interfacial salt bridges in DNA recognition, we compared the salt-dependent binding by Elk1 with ETV6, two ETS homologs differing markedly in ionic composition. While both proteins exhibit salt-dependent binding with cognate DNA that corresponds to interfacial phosphate contacts, their nonspecific binding diverges from cognate binding as well as each other. Molecular dynamics simulations in explicit solvent, which generated ionic interactions in agreement with the experimental binding data, revealed distinct salt-bridge dynamics in the nonspecific complexes formed by the two proteins. Impaired DNA contact by ETV6 resulted in fewer backbone contacts in the nonspecific complex, while Elk1 exhibited a redistribution of extra-interfacial salt bridges via residues that are non-conserved between the two ETS relatives. Thus, primary structure variation in ionic residues can encode highly differentiated specificity mechanisms in a highly conserved DNA-binding motif.

Constrained chromatin accessibility in PU.1-mutated agammaglobulinemia patients.

The pioneer transcription factor (TF) PU.1 controls hematopoietic cell fate by decompacting stem cell heterochromatin and allowing nonpioneer TFs to enter otherwise inaccessible genomic sites. PU.1 deficiency fatally arrests lymphopoiesis and myelopoiesis in mice, but human congenital PU.1 disorders have not previously been described. We studied six unrelated agammaglobulinemic patients, each harboring a heterozygous mutation (four de novo, two unphased) of SPI1, the gene encoding PU.1. Affected patients lacked circulating B cells and possessed few conventional dendritic cells. Introducing disease-similar SPI1 mutations into human hematopoietic stem and progenitor cells impaired early in vitro B cell and myeloid cell differentiation. Patient SPI1 mutations encoded destabilized PU.1 proteins unable to nuclear localize or bind target DNA. In PU.1-haploinsufficient pro-B cell lines, euchromatin was less accessible to nonpioneer TFs critical for B cell development, and gene expression patterns associated with the pro- to pre-B cell transition were undermined. Our findings molecularly describe a novel form of agammaglobulinemia and underscore PU.1's critical, dose-dependent role as a hematopoietic euchromatin gatekeeper.

CTCF puts a new twist on UV damage and repair in skin cancer.

Somatic mutations in skin cancers are highly enriched at binding sites for CCCTC-binding factor (CTCF). We have discovered that CTCF binding alters the DNA structure to render it more susceptible to UV damage. Elevated UV damage formation at CTCF binding sites, in conjunction with subsequent repair inhibition, promotes UV mutagenesis.

DNA mismatches reveal conformational penalties in protein-DNA recognition.

Transcription factors recognize specific genomic sequences to regulate complex gene-expression programs. Although it is well-established that transcription factors bind to specific DNA sequences using a combination of base readout and shape recognition, some fundamental aspects of protein-DNA binding remain poorly understood1,2. Many DNA-binding proteins induce changes in the structure of the DNA outside the intrinsic B-DNA envelope. However, how the energetic cost that is associated with distorting the DNA contributes to recognition has proven difficult to study, because the distorted DNA exists in low abundance in the unbound ensemble3-9. Here we use a high-throughput assay that we term SaMBA (saturation mismatch-binding assay) to investigate the role of DNA conformational penalties in transcription factor-DNA recognition. In SaMBA, mismatched base pairs are introduced to pre-induce structural distortions in the DNA that are much larger than those induced by changes in the Watson-Crick sequence. Notably, approximately 10% of mismatches increased transcription factor binding, and for each of the 22 transcription factors that were examined, at least one mismatch was found that increased the binding affinity. Mismatches also converted non-specific sites into high-affinity sites, and high-affinity sites into 'super sites' that exhibit stronger affinity than any known canonical binding site. Determination of high-resolution X-ray structures, combined with nuclear magnetic resonance measurements and structural analyses, showed that many of the DNA mismatches that increase binding induce distortions that are similar to those induced by protein binding-thus prepaying some of the energetic cost incurred from deforming the DNA. Our work indicates that conformational penalties are a major determinant of protein-DNA recognition, and reveals mechanisms by which mismatches can recruit transcription factors and thus modulate replication and repair activities in the cell10,11.

Dissecting Dynamic and Hydration Contributions to Sequence-Dependent DNA Minor Groove Recognition.

Sequence selectivity is a critical attribute of DNA-binding ligands and underlines the need for detailed molecular descriptions of binding in representative sequence contexts. We investigated the binding and volumetric properties of DB1976, a model bis(benzimidazole)-selenophene diamidine compound with emerging therapeutic potential in acute myeloid leukemia, debilitating fibroses, and obesity-related liver dysfunction. To sample the scope of cognate DB1976 target sites, we evaluated three dodecameric duplexes spanning >103-fold in binding affinity. The attendant changes in partial molar volumes varied substantially, but not in step with binding affinity, suggesting distinct modes of interactions in these complexes. Specifically, whereas optimal binding was associated with loss of hydration water, low-affinity binding released more hydration water. Explicit-atom molecular dynamics simulations showed that minor groove binding perturbed the conformational dynamics and hydration at the termini and interior of the DNA in a sequence-dependent manner. The impact of these distinct local dynamics on hydration was experimentally validated by domain-specific interrogation of hydration with salt, which probed the charged axial surfaces of oligomeric DNA preferentially over the uncharged termini. Minor groove recognition by DB1976, therefore, generates dynamically distinct domains that can make favorable contributions to hydration release in both high- and low-affinity binding. Because ligand binding at internal sites of DNA oligomers modulates dynamics at the termini, the results suggest both short- and long-range dynamic effects along the DNA target that can influence their effectiveness as low-MW competitors of protein binding.

Intrinsic disorder controls two functionally distinct dimers of the master transcription factor PU.1.

Transcription factors comprise a major reservoir of conformational disorder in the eukaryotic proteome. The hematopoietic master regulator PU.1 presents a well-defined model of the most common configuration of intrinsically disordered regions (IDRs) in transcription factors. We report that the structured DNA binding domain (DBD) of PU.1 regulates gene expression via antagonistic dimeric states that are reciprocally controlled by cognate DNA on the one hand and by its proximal anionic IDR on the other. The two conformers are mediated by distinct regions of the DBD without structured contributions from the tethered IDRs. Unlike DNA-bound complexes, the unbound dimer is markedly destabilized. Dimerization without DNA is promoted by progressive phosphomimetic substitutions of IDR residues that are phosphorylated in immune activation and stimulated by anionic crowding agents. These results suggest a previously unidentified, nonstructural role for charged IDRs in conformational control by mitigating electrostatic penalties that would mask the interactions of highly cationic DBDs.

Modulating DNA by polyamides to regulate transcription factor PU.1-DNA binding interactions.

Hairpin polyamides are synthetic small molecules that bind DNA minor groove sequence-selectively and, in many sequences, induce widening of the minor groove and compression of the major groove. The structural distortion of DNA caused by polyamides has enhanced our understanding of the regulation of DNA-binding proteins via polyamides. Polyamides have DNA binding affinities that are comparable to those proteins, therefore, can potentially be used as therapeutic agents to treat diseases caused by aberrant gene expression. In fact, many diseases are characterized by over- or under-expressed genes. PU.1 is a transcription factor that regulates many immune system genes. Aberrant expression of PU.1 has been associated with the development of acute myeloid leukemia (AML). We have, therefore, designed and synthesized ten hairpin polyamides to investigate their capacity in controlling the PU.1-DNA interaction. Our results showed that nine of the polyamides disrupt PU.1-DNA binding and the inhibition capacity strongly correlates with binding affinity. One molecule, FH1024, was observed forming a FH1024-PU.1-DNA ternary complex instead of inhibiting PU.1-DNA binding. This is the first report of a small molecule that is potentially a weak agonist that recruits PU.1 to DNA. This finding sheds light on the design of polyamides that exhibit novel regulatory mechanisms on protein-DNA binding.

Mechanism of cognate sequence discrimination by the ETS-family transcription factor ETS-1.

Functional evidence increasingly implicates low-affinity DNA recognition by transcription factors as a general mechanism for the spatiotemporal control of developmental genes. Although the DNA sequence requirements for affinity are well-defined, the dynamic mechanisms that execute cognate recognition are much less resolved. To address this gap, here we examined ETS1, a paradigm developmental transcription factor, as a model for which cognate discrimination remains enigmatic. Using molecular dynamics simulations, we interrogated the DNA-binding domain of murine ETS1 alone and when bound to high-and low-affinity cognate sites or to nonspecific DNA. The results of our analyses revealed collective backbone and side-chain motions that distinguished cognate versus nonspecific as well as high- versus low-affinity cognate DNA binding. Combined with binding experiments with site-directed ETS1 mutants, the molecular dynamics data disclosed a triad of residues that respond specifically to low-affinity cognate DNA. We found that a DNA-contacting residue (Gln-336) specifically recognizes low-affinity DNA and triggers the loss of a distal salt bridge (Glu-343/Arg-378) via a large side-chain motion that compromises the hydrophobic packing of two core helices. As an intact Glu-343/Arg-378 bridge is the default state in unbound ETS1 and maintained in high-affinity and nonspecific complexes, the low-affinity complex represents a unique conformational adaptation to the suboptimization of developmental enhancers.

Quantifying length-dependent DNA end-binding by nucleoproteins.

The ends of nucleic acids oligomers alter the statistics of interior nonspecific ligand binding and act as binding sites with altered properties. While the former aspect of "end effects" has received much theoretical attention in the literature, the physical nature of end-binding, and hence its potential impact on a wide range of studies with oligomers, remains poorly known. Here, we report for the first time end-binding to DNA using a model helix-turn-helix motif, the DNA-binding domain of ETV6, as a function of DNA sequence length. Spectral analysis of ETV6 intrinsic tryptophan fluorescence by singular value decomposition showed that end-binding to nonspecific fragments was negligible at >0.2 kbp and accumulated to 8% of total binding to 23-bp oligomers. The affinity for end-binding was insensitive to salt but tracked the affinity of interior binding, suggesting translocation from interior sites rather than free solution as its mechanism. As the presence of a cognate site in the 23-bp oligomer suppressed end-binding, neglect of end-binding to the short cognate DNA does not introduce significant error. However, the same applies to nonspecific DNA only if longer fragments (>0.2 kbp) are used.