Redefining the bacteriophage mv4 site-specific recombination system and the sequence specificity of its attB and core-attP sites

Through their involvement in the integration and excision of a large number of mobile genetic elements, such as phages and integrative and conjugative elements (ICEs), site-specific recombination systems based on heterobivalent tyrosine recombinases play a major role in genome dynamics and evolution. However, despite hundreds of these systems having been identified in genome databases, very few have been described in detail, with none from phages that infect Bacillota (formerly Firmicutes). In this study, we reanalyzed the recombination module of Lactobacillus delbrueckii subsp. bulgaricus phage mv4, previously considered atypical compared with classical systems. Our results reveal that mv4 integrase is a 369 aa protein with all the structural hallmarks of recombinases from the Tn916 family and that it cooperatively interacts with its recombination sites. Using randomized DNA libraries, NGS sequencing, and other molecular approaches, we show that the 21-bp core-attP and attB sites have structural similarities to classical systems only if considering the nucleotide degeneracy, with two 7-bp inverted regions corresponding to mv4Int core-binding sites surrounding a 7-bp strand-exchange region. We also examined the different compositional constraints in the core-binding regions, which define the sequence space of permissible recombination sites.

Pleiotropic Effects of PhaR Regulator in Bradyrhizobium diazoefficiens Microaerobic Metabolism

Bradyrhizobium diazoefficiens can live inside soybean root nodules and in free-living conditions. In both states, when oxygen levels decrease, cells adjust their protein pools by gene transcription modulation. PhaR is a transcription factor involved in polyhydroxyalkanoate (PHA) metabolism but also plays a role in the microaerobic network of this bacterium. To deeply uncover the function of PhaR, we applied a multipronged approach, including the expression profile of a phaR mutant at the transcriptional and protein levels under microaerobic conditions, and the identification of direct targets and of proteins associated with PHA granules. Our results confirmed a pleiotropic function of PhaR, affecting several phenotypes, in addition to PHA cycle control. These include growth deficiency, regulation of carbon and nitrogen allocation, and bacterial motility. Interestingly, PhaR may also modulate the microoxic-responsive regulatory network by activating the expression of fixK2 and repressing nifA, both encoding two transcription factors relevant for microaerobic regulation. At the molecular level, two PhaR-binding motifs were predicted and direct control mediated by PhaR determined by protein-interaction assays revealed seven new direct targets for PhaR. Finally, among the proteins associated with PHA granules, we found PhaR, phasins, and other proteins, confirming a dual function of PhaR in microoxia.

The Impact of Lens Epithelium-Derived Growth Factor p75 Dimerization on Its Tethering Function

The transcriptional co-activator lens epithelium-derived growth factor/p75 (LEDGF/p75) plays an important role in the biology of the cell and in several human diseases, including MLL-rearranged acute leukemia, autoimmunity, and HIV-1 infection. In both health and disease, LEDGF/p75 functions as a chromatin tether that interacts with proteins such as MLL1 and HIV-1 integrase via its integrase-binding domain (IBD) and with chromatin through its N-terminal PWWP domain. Recently, dimerization of LEDGF/p75 was shown, mediated by a network of electrostatic contacts between amino acids from the IBD and the C-terminal α6-helix. Here, we investigated the functional impact of LEDGF/p75 variants on the dimerization using biochemical and cellular interaction assays. The data demonstrate that the C-terminal α6-helix folds back in cis on the IBD of monomeric LEDGF/p75. We discovered that the presence of DNA stimulates LEDGF/p75 dimerization. LEDGF/p75 dimerization enhances binding to MLL1 but not to HIV-1 integrase, a finding that was observed in vitro and validated in cell culture. Whereas HIV-1 replication was not dependent on LEDGF/p75 dimerization, colony formation of MLLr-dependent human leukemic THP-1 cells was. In conclusion, our data indicate that intricate changes in the quaternary structure of LEDGF/p75 modulate its tethering function.

ULDNA: integrating unsupervised multi-source language models with LSTM-attention network for high-accuracy protein-DNA binding site prediction

Efficient and accurate recognition of protein-DNA interactions is vital for understanding the molecular mechanisms of related biological processes and further guiding drug discovery. Although the current experimental protocols are the most precise way to determine protein-DNA binding sites, they tend to be labor-intensive and time-consuming. There is an immediate need to design efficient computational approaches for predicting DNA-binding sites. Here, we proposed ULDNA, a new deep-learning model, to deduce DNA-binding sites from protein sequences. This model leverages an LSTM-attention architecture, embedded with three unsupervised language models that are pre-trained on large-scale sequences from multiple database sources. To prove its effectiveness, ULDNA was tested on 229 protein chains with experimental annotation of DNA-binding sites. Results from computational experiments revealed that ULDNA significantly improves the accuracy of DNA-binding site prediction in comparison with 17 state-of-the-art methods. In-depth data analyses showed that the major strength of ULDNA stems from employing three transformer language models. Specifically, these language models capture complementary feature embeddings with evolution diversity, in which the complex DNA-binding patterns are buried. Meanwhile, the specially crafted LSTM-attention network effectively decodes evolution diversity-based embeddings as DNA-binding results at the residue level. Our findings demonstrated a new pipeline for predicting DNA-binding sites on a large scale with high accuracy from protein sequence alone.

Protein-DNA binding sites prediction based on pre-trained protein language model and contrastive learning

Protein-DNA interaction is critical for life activities such as replication, transcription and splicing. Identifying protein-DNA binding residues is essential for modeling their interaction and downstream studies. However, developing accurate and efficient computational methods for this task remains challenging. Improvements in this area have the potential to drive novel applications in biotechnology and drug design. In this study, we propose a novel approach called Contrastive Learning And Pre-trained Encoder (CLAPE), which combines a pre-trained protein language model and the contrastive learning method to predict DNA binding residues. We trained the CLAPE-DB model on the protein-DNA binding sites dataset and evaluated the model performance and generalization ability through various experiments. The results showed that the area under ROC curve values of the CLAPE-DB model on the two benchmark datasets reached 0.871 and 0.881, respectively, indicating superior performance compared to other existing models. CLAPE-DB showed better generalization ability and was specific to DNA-binding sites. In addition, we trained CLAPE on different protein-ligand binding sites datasets, demonstrating that CLAPE is a general framework for binding sites prediction. To facilitate the scientific community, the benchmark datasets and codes are freely available at https://github.com/YAndrewL/clape.

Cooperative action of separate interaction domains promotes high-affinity DNA binding of Arabidopsis thaliana ARF transcription factors

The signaling molecule auxin coordinates many growth and development processes in plants, mainly through modulating gene expression. Transcriptional response is mediated by the family of auxin response factors (ARF). Monomers of this family recognize a DNA motif and can homodimerize through their DNA-binding domain (DBD), enabling cooperative binding to an inverted binding site. Most ARFs further contain a C-terminal PB1 domain that is capable of homotypic interactions and mediating interactions with Aux/IAA repressors. Given the dual role of the PB1 domain, and the ability of both DBD and PB1 domain to mediate dimerization, a key question is how these domains contribute to DNA-binding specificity and affinity. So far, ARF-ARF and ARF-DNA interactions have mostly been approached using qualitative methods that do not provide a quantitative and dynamic view on the binding equilibria. Here, we utilize a DNA binding assay based on single-molecule Förster resonance energy transfer (smFRET) to study the affinity and kinetics of the interaction of several Arabidopsis thaliana ARFs with an IR7 auxin-responsive element (AuxRE). We show that both DBD and PB1 domains of AtARF2 contribute toward DNA binding, and we identify ARF dimer stability as a key parameter in defining binding affinity and kinetics across AtARFs. Lastly, we derived an analytical solution for a four-state cyclic model that explains both the kinetics and the affinity of the interaction between AtARF2 and IR7. Our work demonstrates that the affinity of ARFs toward composite DNA response elements is defined by dimerization equilibrium, identifying this as a key element in ARF-mediated transcriptional activity.

Structural and Molecular Kinetic Features of Activities of DNA Polymerases

DNA polymerases catalyze DNA synthesis during the replication, repair, and recombination of DNA. Based on phylogenetic analysis and primary protein sequences, DNA polymerases have been categorized into seven families: A, B, C, D, X, Y, and RT. This review presents generalized data on the catalytic mechanism of action of DNA polymerases. The structural features of different DNA polymerase families are described in detail. The discussion highlights the kinetics and conformational dynamics of DNA polymerases from all known polymerase families during DNA synthesis.

Residues located in the primase domain of the bacteriophage T7 primase-helicase are essential for loading the hexameric complex onto DNA

The T7 primase-helicase plays a pivotal role in the replication of T7 DNA. Using affinity isolation of peptide–nucleic acid crosslinks and mass spectrometry, we identify protein regions in the primase-helicase and T7 DNA polymerase that form contacts with the RNA primer and DNA template. The contacts between nucleic acids and the primase domain of the primase-helicase are centered in the RNA polymerase subdomain of the primase domain, in a cleft between the N-terminal subdomain and the topoisomerase-primase fold. We demonstrate that residues along a beta sheet in the N-terminal subdomain that contacts the RNA primer are essential for phage growth and primase activity in vitro. Surprisingly, we found mutations in the primase domain that had a dramatic effect on the helicase. Substitution of a residue conserved in other DnaG-like enzymes, R84A, abrogates both primase and helicase enzymatic activities of the T7 primase-helicase. Alterations in this residue also decrease binding of the primase-helicase to ssDNA. However, mass photometry measurements show that these mutations do not interfere with the ability of the protein to form the active hexamer.

Fine-Tuning Modulation of Oxidation-Mediated Posttranslational Control of Bradyrhizobium diazoefficiens FixK2 Transcription Factor

FixK2 is a CRP/FNR-type transcription factor that plays a central role in a sophisticated regulatory network for the anoxic, microoxic and symbiotic lifestyles of the soybean endosymbiont Bradyrhizobium diazoefficiens. Aside from the balanced expression of the fixK2 gene under microoxic conditions (induced by the two-component regulatory system FixLJ and negatively auto-repressed), FixK2 activity is posttranslationally controlled by proteolysis, and by the oxidation of a singular cysteine residue (C183) near its DNA-binding domain. To simulate the permanent oxidation of FixK2, we replaced C183 for aspartic acid. Purified C183D FixK2 protein showed both low DNA binding and in vitro transcriptional activation from the promoter of the fixNOQP operon, required for respiration under symbiosis. However, in a B. diazoefficiens strain coding for C183D FixK2, expression of a fixNOQP'-'lacZ fusion was similar to that in the wild type, when both strains were grown microoxically. The C183D FixK2 encoding strain also showed a wild-type phenotype in symbiosis with soybeans, and increased fixK2 gene expression levels and FixK2 protein abundance in cells. These two latter observations, together with the global transcriptional profile of the microoxically cultured C183D FixK2 encoding strain, suggest the existence of a finely tuned regulatory strategy to counterbalance the oxidation-mediated inactivation of FixK2 in vivo.

Development of a FOXM1-DBD Binding Assay for High-Throughput Screening Using TR-FRET Assay

Electrophoretic mobility shift assay (EMSA) technology has been widely employed for the analysis of transcription factors such as Forkhead box protein M1 (FOXM1). However, the application of high-throughput screening (HTS) in performing, such analyses are limited as it uses time consuming electrophoresis procedure and radioisotopes. In this study, we developed a FOXM1-DNA binding domain (DBD) binding assay based on time-resolved fluorescence energy transfer (TR-FRET) that enables HTS for the inhibitors of FOXM1-DNA interaction. This assay was robust, highly reproducible and could be easily miniaturized into 384-well plate format. The signal-to-background (S/B) ratio and Z' factor were calculated as 7.46 and 0.74, respectively, via a series of optimization of the assay conditions. A pilot library screening of 1019 natural compounds was performed using the FOXM1-DBD binding assay. Five hit compounds, namely, AC1LXM, BRN5, gangaleoidin, leoidin, and roemerine were identified as the inhibitors of FOXM1. In a cell viability assay, it was demonstrated that cell proliferation of FOXM1 overexpressed cell lines was suppressed in cell lines such as MDA-MB-231 and MCF-7 by five hit compounds. These results indicate that developed FOXM1-DBD binding assay can be applied to highly efficiency HTS of compound libraries.

A feature-based approach to predict hot spots in protein-DNA binding interfaces

DNA-binding hot spot residues of proteins are dominant and fundamental interface residues that contribute most of the binding free energy of protein-DNA interfaces. As experimental methods for identifying hot spots are expensive and time consuming, computational approaches are urgently required in predicting hot spots on a large scale. In this work, we systematically assessed a wide variety of 114 features from a combination of the protein sequence, structure, network and solvent accessible information and their combinations along with various feature selection strategies for hot spot prediction. We then trained and compared four commonly used machine learning models, namely, support vector machine (SVM), random forest, Naïve Bayes and k-nearest neighbor, for the identification of hot spots using 10-fold cross-validation and the independent test set. Our results show that (1) features based on the solvent accessible surface area have significant effect on hot spot prediction; (2) different but complementary features generally enhance the prediction performance; and (3) SVM outperforms other machine learning methods on both training and independent test sets. In an effort to improve predictive performance, we developed a feature-based method, namely, PrPDH (Prediction of Protein-DNA binding Hot spots), for the prediction of hot spots in protein-DNA binding interfaces using SVM based on the selected 10 optimal features. Comparative results on benchmark data sets indicate that our predictor is able to achieve generally better performance in predicting hot spots compared to the state-of-the-art predictors. A user-friendly web server for PrPDH is well established and is freely available at http://bioinfo.ahu.edu.cn:8080/PrPDH.

Characterization of the specific DNA-binding properties of Tnp26, the transposase of insertion sequence IS26

The bacterial insertion sequence (IS) IS26 mobilizes and disseminates antibiotic resistance genes. It differs from bacterial IS that have been studied to date as it exclusively forms cointegrates via either a copy-in (replicative) or a recently discovered targeted conservative mode. To investigate how the Tnp26 transposase recognizes the 14-bp terminal inverted repeats (TIRs) that bound the IS, amino acids in two domains in the N-terminal (amino acids M1-P56) region were replaced. These changes substantially reduced cointegration in both modes. Tnp26 was purified as a maltose-binding fusion protein and shown to bind specifically to dsDNA fragments that included an IS26 TIR. However, Tnp26 with an R49A or a W50A substitution in helix 3 of a predicted trihelical helix-turn-helix domain (amino acids I13-R53) or an F4A or F9A substitution replacing the conserved amino acids in a unique disordered N-terminal domain (amino acids M1-D12) did not bind. The N-terminal M1-P56 fragment also bound to the TIR but only at substantially higher concentrations, indicating that other parts of Tnp26 enhance the binding affinity. The binding site was confined to the internal part of the TIR, and a G to T nucleotide substitution in the TGT at positions 6 to 8 of the TIR that is conserved in most IS26 family members abolished binding of both Tnp26 (M1-M234) and Tnp26 M1-P56 fragment. These findings indicate that the helix-turn-helix and disordered domains of Tnp26 play a role in Tnp26-TIR complex formation. Both domains are conserved in all members of the IS26 family.

DeepDISE: DNA Binding Site Prediction Using a Deep Learning Method

It is essential for future research to develop a new, reliable prediction method of DNA binding sites because DNA binding sites on DNA-binding proteins provide critical clues about protein function and drug discovery. However, the current prediction methods of DNA binding sites have relatively poor accuracy. Using 3D coordinates and the atom-type of surface protein atom as the input, we trained and tested a deep learning model to predict how likely a voxel on the protein surface is to be a DNA-binding site. Based on three different evaluation datasets, the results show that our model not only outperforms several previous methods on two commonly used datasets, but also demonstrates its robust performance to be consistent among the three datasets. The visualized prediction outcomes show that the binding sites are also mostly located in correct regions. We successfully built a deep learning model to predict the DNA binding sites on target proteins. It demonstrates that 3D protein structures plus atom-type information on protein surfaces can be used to predict the potential binding sites on a protein. This approach should be further extended to develop the binding sites of other important biological molecules.

Structural basis for the DNA-binding activity of human ARID4B Tudor domain

Human ARID4A and ARID4B are homologous proteins that are important in controlling gene expression and epigenetic regulation but have distinct functions. Previous studies have shown that the N-terminal domain of ARID4A is an unusual interdigitated double Tudor domain with DNA-binding activity. However, how the Tudor domain of ARID4B differs from that of ARID4A remains unknown. Here, we found that the ARID4B Tudor domain has significantly weaker DNA affinity than the ARID4A Tudor domain despite sharing more than 80% sequence identity. Structure determination and DNA titration analysis indicated that the ARID4B Tudor domain is also an interdigitated double Tudor domain with a DNA-binding surface similar to ARID4A. We identified a residue close to the DNA-binding site of the Tudor domain that differs between ARID4A and ARID4B. The Leu50 in ARID4A is Glu50 in ARID4B, and the latter forms salt bridges with two lysine residues at the DNA-binding surface. This causes a decrease in the strength of positive charge, thus reducing DNA-binding affinity while significantly increasing protein stability. We also found that a C-terminal extension region enhances the DNA-binding affinity of the ARID4B Tudor domain. This C-terminal extension is disordered and contains a positively charged RGR motif, providing an additional DNA-binding site. Finally, sequence and phylogenetic analyses indicated that the residue differences and the presence of the RGR extension region are conserved. These results provide new insight into the functional differences between ARID4A and ARID4B proteins, as well as elucidating the function of the disordered regions in these proteins.

Single-molecule multiplexed profiling of protein-DNA complexes using magnetic tweezers

Epigenetics, such as the dynamic interplay between DNA methylation and demethylation, play diverse roles in critical cellular events. Enzymatic activity at CpG sites, where cytosines are methylated or demethylated, is known to be influenced by the density of CpGs, methylation states, and the flanking sequences of a CpG site. However, how the relevant enzymes are recruited to and recognize their target DNA is less clear. Moreover, although DNA-binding epigenetic enzymes are ideal targets for therapeutic intervention, these targets have been rarely exploited. Single-molecule techniques offer excellent capabilities to probe site-specific protein-DNA interactions and unravel the dynamics. Here, we develop a single-molecule approach that allows multiplexed profiling of protein-DNA complexes using magnetic tweezers. When a DNA hairpin with multiple binding sites is unzipping, strand separation pauses at the positions bound by a protein. We can thus measure site-specific binding probabilities and dissociation time directly. Taking the TET1 CXXC domain as an example, we show that TET1 CXXC binds multiple CpG motifs with various flanking nucleotides or different methylation patterns in an AT-rich DNA. We are able to establish for the first time, at nanometer resolution, that TET1 CXXC prefers G/C flanked CpG motif over C/G, A/T, or T/A flanked ones. CpG methylation strengthens TET1 CXXC recruitment but has little effect on dissociation time. Finally, we demonstrate that TET1 CXXC can distinguish five CpG clusters in a CpG island with crowded binding motifs. We anticipate that the feasibility of single-molecule multiplexed profiling assays will contribute to the understanding of protein-DNA interactions.

Poly(ADP-ribose) polymerase 1 regulates mitochondrial DNA repair in an NAD-dependent manner

Mitochondrial DNA is located in organelle that house essential metabolic reactions and contains high reactive oxygen species. Therefore, mitochondrial DNA suffers more oxidative damage than its nuclear counterpart. Formation of a repair enzyme complex is beneficial to DNA repair. Recent studies have shown that mitochondrial DNA polymerase (Pol γ) and poly(ADP-ribose) polymerase 1 (PARP1) were found in the same complex along with other mitochondrial DNA repair enzymes, and mitochondrial PARP1 level is correlated with mtDNA integrity. However, the molecular basis for the functional connection between Pol γ and PARP1 has not yet been elucidated because cellular functions of PARP1 in DNA repair are intertwined with metabolism via NAD+ (nicotinamide adenosine dinucleotide), the substrate of PARP1, and a metabolic cofactor. To dissect the direct effect of PARP1 on mtDNA from the secondary perturbation of metabolism, we report here biochemical studies that recapitulated Pol γ PARylation observed in cells and showed that PARP1 regulates Pol γ activity during DNA repair in a metabolic cofactor NAD⁺ (nicotinamide adenosine dinucleotide)-dependent manner. In the absence of NAD⁺, PARP1 completely inhibits Pol γ, while increasing NAD⁺ levels to a physiological concentration that enables Pol γ to resume maximum repair activity. Because cellular NAD+ levels are linked to metabolism and to ATP production via oxidative phosphorylation, our results suggest that mtDNA damage repair is coupled to cellular metabolic state and the integrity of the respiratory chain.

Structure of HIRAN domain of human HLTF bound to duplex DNA provides structural basis for DNA unwinding to initiate replication fork regression

Replication fork regression is a mechanism to rescue a stalled fork by various replication stresses, such as DNA lesions. Helicase-like transcription factor, a SNF2 translocase, plays a central role in the fork regression and its N-terminal domain, HIRAN (HIP116 and Rad5 N-terminal), binds the 3'-hydroxy group of single-stranded DNA. Furthermore, HIRAN is supposed to bind double-stranded DNA (dsDNA) and involved in strand separation in the fork regression, whereas structural basis for mechanisms underlying dsDNA binding and strand separation by HIRAN are still unclear. Here, we report the crystal structure of HIRAN bound to duplex DNA. The structure reveals that HIRAN binds the 3'-hydroxy group of DNA and unexpectedly unwinds three nucleobases of the duplex. Phe-142 is involved in the dsDNA binding and the strand separation. In addition, the structure unravels the mechanism underlying sequence-independent recognition for purine bases by HIRAN, where the N-glycosidic bond adopts syn conformation. Our findings indicate direct involvement of HIRAN in the fork regression by separating of the daughter strand from the parental template.

Architecture of DNA elements mediating ARF transcription factor binding and auxin-responsive gene expression in Arabidopsis

The hormone auxin controls many aspects of the plant life cycle by regulating the expression of thousands of genes. The transcriptional output of the nuclear auxin signaling pathway is determined by the activity of AUXIN RESPONSE transcription FACTORs (ARFs), through their binding to cis-regulatory elements in auxin-responsive genes. Crystal structures, in vitro, and heterologous studies have fueled a model in which ARF dimers bind with high affinity to distinctly spaced repeats of canonical AuxRE motifs. However, the relevance of this "caliper" model, and the mechanisms underlying the binding affinities in vivo, have remained elusive. Here we biochemically and functionally interrogate modes of ARF-DNA interaction. We show that a single additional hydrogen bond in Arabidopsis ARF1 confers high-affinity binding to individual DNA sites. We demonstrate the importance of AuxRE cooperativity within repeats in the Arabidopsis TMO5 and IAA11 promoters in vivo. Meta-analysis of transcriptomes further reveals strong genome-wide association of auxin response with both inverted (IR) and direct (DR) AuxRE repeats, which we experimentally validated. The association of these elements with auxin-induced up-regulation (DR and IR) or down-regulation (IR) was correlated with differential binding affinities of A-class and B-class ARFs, respectively, suggesting a mechanistic basis for the distinct activity of these repeats. Our results support the relevance of high-affinity binding of ARF transcription factors to uniquely spaced DNA elements in vivo, and suggest that differential binding affinities of ARF subfamilies underlie diversity in cis-element function.

High-speed atomic force microscopy directly visualizes conformational dynamics of the HIV Vif protein in complex with three host proteins

Vif (viral infectivity factor) is a protein that is essential for the replication of the HIV-1 virus. The key function of Vif is to disrupt the antiviral activity of host APOBEC3 (apolipoprotein B mRNA-editing enzyme catalytic subunit 3) proteins, which mutate viral nucleic acids. Inside the cell, Vif binds to the host cell proteins Elongin-C, Elongin-B, and core-binding factor subunit β, forming a four-protein complex called VCBC. The structure of VCBC-Cullin5 has recently been solved by X-ray crystallography, and, using molecular dynamics simulations, the dynamics of VCBC have been characterized. Here, we applied time-lapse high-speed atomic force microscopy to visualize the conformational changes of the VCBC complex. We determined the three most favorable conformations of this complex, which we identified as the triangle, dumbbell, and globular structures. Moreover, we characterized the dynamics of each of these structures. Our data revealed the very dynamic behavior of all of them, with the triangle and dumbbell structures being the most dynamic. These findings provide insight into the structure and dynamics of the VCBC complex and may support efforts to improve HIV treatment, because Vif is essential for virus survival in the cell.

Structure of the HLTF HIRAN domain and its functional implications in regression of a stalled replication fork

HLTF (helicase-like transcription factor) is a yeast RAD5 homolog that is found in mammals. HLTF has E3 ubiquitin ligase and DNA helicase activities, and is a pivotal protein in template-switched DNA synthesis that allows DNA replication to continue even in the presence of DNA damage by utilizing a newly synthesized undamaged strand as a template. In addition, HLTF has a DNA-binding domain termed HIRAN (HIP116 and RAD5 N-terminal). HIRAN has been hypothesized to play a role in DNA binding; however, the structural basis of its role in DNA binding has remained unclear. In the past five years, several crystal structures of HIRAN have been reported. These structures revealed new insights into the molecular mechanism underlying DNA binding by HIRAN. Here, the structural information on HIRAN is summarized and the function of HIRAN in recognizing the 3'-terminus of the daughter strand at a stalled replication fork and the implications for its involvement in fork regression are discussed.

Mechanistic insights into heme-mediated transcriptional regulation via a bacterial manganese-binding iron regulator, iron response regulator (Irr)

The transcription factor iron response regulator (Irr) is a key regulator of iron homeostasis in the nitrogen-fixating bacterium Bradyrhizobium japonicum Irr acts by binding to target genes, including the iron control element (ICE), and is degraded in response to heme binding. Here, we examined this binding activity using fluorescence anisotropy with a 6-carboxyfluorescein-labeled ICE-like oligomer (FAM-ICE). In the presence of Mn²⁺, Irr addition increased the fluorescence anisotropy, corresponding to formation of the Irr-ICE complex. The addition of EDTA to the Irr-ICE complex reduced fluorescence anisotropy, but fluorescence was recovered after Mn²⁺ addition, indicating that Mn²⁺ binding is a prerequisite for complex formation. Binding activity toward ICE was lost upon introduction of substitutions in a His-cluster region of Irr, revealing that Mn²⁺ binds to this region. We observed that the His-cluster region is also the heme binding site; results from fluorescence anisotropy and electrophoretic mobility shift analyses disclosed that the addition of a half-equivalent of heme dissociates Irr from ICE, likely because of Mn²⁺ release due to heme binding. We hypothesized that heme binding to another heme binding site, Cys-29, would also inhibit the formation of the Irr-ICE complex because it is proximal to the ICE binding site, which was supported by the loss of ICE binding activity in a Cys-29-mutated Irr. These results indicate that Irr requires Mn²⁺ binding to form the Irr-ICE complex and that the addition of heme dissociates Irr from ICE by replacing Mn²⁺ with heme or by heme binding to Cys-29.

Hydrophobic Amino Acids as Universal Elements of Protein-Induced DNA Structure Deformation

Interaction with the DNA minor groove is a significant contributor to specific sequence recognition in selected families of DNA-binding proteins. Based on a statistical analysis of 3D structures of protein–DNA complexes, we propose that distortion of the DNA minor groove resulting from interactions with hydrophobic amino acid residues is a universal element of protein–DNA recognition. We provide evidence to support this by associating each DNA minor groove-binding amino acid residue with the local dimensions of the DNA double helix using a novel algorithm. The widened DNA minor grooves are associated with high GC content. However, some AT-rich sequences contacted by hydrophobic amino acids (e.g., phenylalanine) display extreme values of minor groove width as well. For a number of hydrophobic amino acids, distinct secondary structure preferences could be identified for residues interacting with the widened DNA minor groove. These results hold even after discarding the most populous families of minor groove-binding proteins.

Biophysical Screens Identify Fragments That Bind to the Viral DNA-Binding Proteins EBNA1 and LANA

The human gamma-herpesviruses Epstein-Barr virus (EBV) (HHV-4) and Kaposi's sarcoma-associated herpesvirus (KSHV) (HHV-8) are responsible for a number of diseases, including various types of cancer. Epstein-Barr nuclear antigen 1 (EBNA1) from EBV and latency-associated nuclear antigen (LANA) from KSHV are viral-encoded DNA-binding proteins that are essential for the replication and maintenance of their respective viral genomes during latent, oncogenic infection. As such, EBNA1 and LANA are attractive targets for the development of small-molecule inhibitors. To this end, we performed a biophysical screen of EBNA1 and LANA using a fragment library by saturation transfer difference (STD)-NMR spectroscopy and surface plasmon resonance (SPR). We identified and validated a number of unique fragment hits that bind to EBNA1 or LANA. We also determined the high-resolution crystal structure of one fragment bound to EBNA1. Results from this screening cascade provide new chemical starting points for the further development of potent inhibitors for this class of viral proteins.

The C-terminal tails of the mitochondrial transcription factors Mtf1 and TFB2M are part of an autoinhibitory mechanism that regulates DNA binding

The structurally homologous Mtf1 and TFB2M proteins serve as transcription initiation factors of mitochondrial RNA polymerases in Saccharomyces cerevisiae and humans, respectively. These transcription factors directly interact with the nontemplate strand of the transcription bubble to drive promoter melting. Given the key roles of Mtf1 and TFB2M in promoter-specific transcription initiation, it can be expected that the DNA binding activity of the mitochondrial transcription factors is regulated to prevent DNA binding at inappropriate times. However, little information is available on how mitochondrial DNA transcription is regulated. While studying C-terminal (C-tail) deletion mutants of Mtf1 and TFB2M, we stumbled upon a finding that suggested that the flexible C-tail region of these factors autoregulates their DNA binding activity. Quantitative DNA binding studies with fluorescence anisotropy-based titrations revealed that Mtf1 with an intact C-tail has no affinity for DNA but deletion of the C-tail greatly increases Mtf1's DNA binding affinity. Similar observations were made with TFB2M, although autoinhibition by the C-tail of TFB2M was not as complete as in Mtf1. Analysis of available TFB2M structures disclosed that the C-tail engages in intramolecular interactions with the DNA binding groove in the free factor, which, we propose, inhibits its DNA binding activity. Further experiments showed that RNA polymerase relieves this autoinhibition by interacting with the C-tail and engaging it in complex formation. In conclusion, our biochemical and structural analyses reveal autoinhibitory and activation mechanisms of mitochondrial transcription factors that regulate their DNA binding activities and aid in specific assembly of transcription initiation complexes.

Atomic force microscopy-based characterization of the interaction of PriA helicase with stalled DNA replication forks

In bacteria, the restart of stalled DNA replication forks requires the DNA helicase PriA. PriA can recognize and remodel abandoned DNA replication forks, unwind DNA in the 3'-to-5' direction, and facilitate the loading of the helicase DnaB onto the DNA to restart replication. Single-stranded DNA-binding protein (SSB) is typically present at the abandoned forks, but it is unclear how SSB and PriA interact, although it has been shown that the two proteins interact both physically and functionally. Here, we used atomic force microscopy to visualize the interaction of PriA with DNA substrates with or without SSB. These experiments were done in the absence of ATP to delineate the substrate recognition pattern of PriA before its ATP-catalyzed DNA-unwinding reaction. These analyses revealed that in the absence of SSB, PriA binds preferentially to a fork substrate with a gap in the leading strand. Such a preference has not been observed for 5'- and 3'-tailed duplexes, suggesting that it is the fork structure that plays an essential role in PriA's selection of DNA substrates. Furthermore, we found that in the absence of SSB, PriA binds exclusively to the fork regions of the DNA substrates. In contrast, fork-bound SSB loads PriA onto the duplex DNA arms of forks, suggesting a remodeling of PriA by SSB. We also demonstrate that the remodeling of PriA requires a functional C-terminal domain of SSB. In summary, our atomic force microscopy analyses reveal key details in the interactions between PriA and stalled DNA replication forks with or without SSB.

Activating transcription factor 4 (ATF4) promotes skeletal muscle atrophy by forming a heterodimer with the transcriptional regulator C/EBPβ

Skeletal muscle atrophy is a highly-prevalent and debilitating condition that remains poorly understood at the molecular level. Previous work found that aging, fasting, and immobilization promote skeletal muscle atrophy via expression of activating transcription factor 4 (ATF4) in skeletal muscle fibers. However, the direct biochemical mechanism by which ATF4 promotes muscle atrophy is unknown. ATF4 is a member of the basic leucine zipper transcription factor (bZIP) superfamily. Because bZIP transcription factors are obligate dimers, and because ATF4 is unable to form highly-stable homodimers, we hypothesized that ATF4 may promote muscle atrophy by forming a heterodimer with another bZIP family member. To test this hypothesis, we biochemically isolated skeletal muscle proteins that associate with the dimerization- and DNA-binding domain of ATF4 (the bZIP domain) in mouse skeletal muscle fibers in vivo Interestingly, we found that ATF4 forms at least five distinct heterodimeric bZIP transcription factors in skeletal muscle fibers. Furthermore, one of these heterodimers, composed of ATF4 and CCAAT enhancer-binding protein β (C/EBPβ), mediates muscle atrophy. Within skeletal muscle fibers, the ATF4-C/EBPβ heterodimer interacts with a previously unrecognized and evolutionarily conserved ATF-C/EBP composite site in exon 4 of the Gadd45a gene. This three-way interaction between ATF4, C/EBPβ, and the ATF-C/EBP composite site activates the Gadd45a gene, which encodes a critical mediator of muscle atrophy. Together, these results identify a biochemical mechanism by which ATF4 induces skeletal muscle atrophy, providing molecular-level insights into the etiology of skeletal muscle atrophy.

Structural basis of binding of homodimers of the nuclear receptor NR4A2 to selective Nur-responsive DNA elements

Proteins of nuclear receptor subfamily 4 group A (NR4A), including NR4A1/NGFI-B, NR4A2/Nurr1, and NR4A3/NOR-1, are nuclear transcription factors that play important roles in metabolism, apoptosis, and proliferation. NR4A proteins recognize DNA response elements as monomers or dimers to regulate the transcription of a variety of genes involved in multiple biological processes. In this study, we determined two crystal structures of the NR4A2 DNA-binding domain (NR4A2-DBD) bound to two Nur-responsive elements: an inverted repeat and an everted repeat at 2.6-2.8 Å resolution. The structures revealed that two NR4A2-DBD molecules bind independently to the everted repeat, whereas two other NR4A2-DBD molecules form a novel dimer interface on the inverted repeat. Moreover, substitution of the interfacial residue valine 298 to lysine as well as mutation of DNA bases involved in the interactions abolished the dimerization. Overall, our structural, biochemical, and bioinformatics analyses provide a molecular basis for the binding of the NR4A2 protein dimers to NurREs and advance our understanding of the dimerization specificity of nuclear receptors.

Understanding the Effect of Structural Diversity in WRKY Transcription Factors on DNA Binding Efficiency through Molecular Dynamics Simulation

A precise understanding of the molecular mechanism involved in stress conditions has great importance for crop improvement. Biomolecules, such as WRKY proteins, which are the largest transcription factor family that is widely distributed in higher plants, plays a significant role in plant defense response against various biotic and abiotic stressors. In the present study, an extensive homology-based three-dimensional model construction and subsequent interaction study of WRKY DNA-binding domain (DBD) in CcWRKY1 (Type I), CcWRKY51 (Type II), and CcWRKY70 (Type III) belonging to pigeonpea, a highly tolerant crop species, was performed. Evaluation of the generated protein models was done to check their reliability and accuracy based on the quantitative and qualitative parameters. The final model was subjected to investigate the comparative binding analysis of different types of WRKY–DBD with DNA-W-box (a cis-acting element) by protein–DNA docking and molecular dynamics (MD) simulation. The DNA binding specificity with WRKY variants was scrutinized through protein–DNA interaction using the HADDOCK server. The stability, as well as conformational changes of protein–DNA complex, was investigated through molecular dynamics (MD) simulations for 100 ns using GROMACS. Additionally, the comparative stability and dynamic behavior of each residue of the WRKY–DBD type were analyzed in terms of root mean square deviation (RMSD), root mean square fluctuation (RMSF)values of the backbone atoms for each frame taking the minimized structure as a reference. The details of DNA binding activity of three different types of WRKY–DBD provided here will be helpful to better understand the regulation of WRKY gene family members in plants.

A Novel G-Quadruplex Binding Protein in Yeast-Slx9

G-quadruplex (G4) structures are highly stable four-stranded DNA and RNA secondary structures held together by non-canonical guanine base pairs. G4 sequence motifs are enriched at specific sites in eukaryotic genomes, suggesting regulatory functions of G4 structures during different biological processes. Considering the high thermodynamic stability of G4 structures, various proteins are necessary for G4 structure formation and unwinding. In a yeast one-hybrid screen, we identified Slx9 as a novel G4-binding protein. We confirmed that Slx9 binds to G4 DNA structures in vitro. Despite these findings, Slx9 binds only insignificantly to G-rich/G4 regions in Saccharomyces cerevisiae as demonstrated by genome-wide ChIP-seq analysis. However, Slx9 binding to G4s is significantly increased in the absence of Sgs1, a RecQ helicase that regulates G4 structures. Different genetic and molecular analyses allowed us to propose a model in which Slx9 recognizes and protects stabilized G4 structures in vivo.

Structural insights into methylated DNA recognition by the C-terminal zinc fingers of the DNA reader protein ZBTB38

Methyl-CpG-binding proteins (MBPs) are selective readers of DNA methylation that play an essential role in mediating cellular transcription processes in both normal and diseased cells. This physiological function of MBPs has generated significant interest in understanding the mechanisms by which these proteins read and interpret DNA methylation signals. Zinc finger and BTB domain-containing 38 (ZBTB38) represents one member of the zinc finger (ZF) family of MBPs. We recently demonstrated that the C-terminal ZFs of ZBTB38 exhibit methyl-selective DNA binding within the ((A/G)TmCG(G/A)(mC/T)(G/A)) context both in vitro and within cells. Here we report the crystal structure of the first four C-terminal ZBTB38 ZFs (ZFs 6-9) in complex with the previously identified methylated consensus sequence at 1.75 Å resolution. From the structure, methyl-selective binding is preferentially localized at the 5' mCpG site of the bound DNA, which is facilitated through a series of base-specific interactions from residues within the α-helices of ZF7 and ZF8. ZF6 and ZF9 primarily stabilize ZF7 and ZF8 to facilitate the core base-specific interactions. Further structural and biochemical analyses, including solution NMR spectroscopy and electrophoretic mobility gel shift assays, revealed that the C-terminal ZFs of ZBTB38 utilize an alternative mode of mCpG recognition from the ZF MBPs structurally evaluated to date. Combined, these findings provide insight into the mechanism by which this ZF domain of ZBTB38 selectively recognizes methylated CpG sites and expands our understanding of how ZF-containing proteins can interpret this essential epigenetic mark.

Transcription Factor ANAC074 Binds to NRS1, NRS2, or MybSt1 Element in Addition to the NACRS to Regulate Gene Expression

NAC (NAM, ATAF1/2, and CUC2) transcription factors play important roles in many biological processes, and mainly bind to the NACRS with core sequences "CACG" or "CATGTG" to regulate gene expression. However, whether NAC proteins can bind to other motifs without these core sequences remains unknown. In this study, we employed a Transcription Factor-Centered Yeast one Hybrid (TF-Centered Y1H) screen to study the motifs recognized by ANAC074. In addition to the NACRS core cis-element, we identified that ANAC074 could bind to MybSt1, NRS1, and NRS2. Y1H and GUS assays showed that ANAC074 could bind the promoters of ethylene responsive genes and stress responsive genes via the NRS1, NRS2, or MybSt1 element. ChIP study further confirmed that the bindings of ANAC074 to MybSt1, NRS1, and NRS2 actually occurred in Arabidopsis. Furthermore, ten NAC proteins from different NAC subfamilies in Arabidopsis thaliana were selected and confirmed to bind to the MybSt1, NRS1, and NRS2 motifs, indicating that they are recognized commonly by NACs. These findings will help us to further reveal the functions of NAC proteins.

Electrostatic repulsion causes anticooperative DNA binding between tumor suppressor ETS transcription factors and JUN-FOS at composite DNA sites

Many different transcription factors (TFs) regulate gene expression in a combinatorial fashion, often by binding in close proximity to each other on composite cis-regulatory DNA elements. Here, we investigated how ETS TFs bind with the AP1 TFs JUN-FOS at composite DNA-binding sites. DNA-binding ability with JUN-FOS correlated with the phenotype of ETS proteins in prostate cancer. We found that the oncogenic ETS-related gene (ERG) and ETS variant (ETV) 1/4/5 subfamilies co-occupy ETS-AP1 sites with JUN-FOS in vitro, whereas JUN-FOS robustly inhibited DNA binding by the tumor suppressors ETS homologous factor (EHF) and SAM pointed domain-containing ETS TF (SPDEF). EHF bound ETS-AP1 DNA with tighter affinity than ERG in the absence of JUN-FOS, possibly enabling EHF to compete with ERG and JUN-FOS for binding to ETS-AP1 sites. Genome-wide mapping of EHF- and ERG-binding sites in prostate epithelial cells revealed that EHF is preferentially excluded from closely spaced ETS-AP1 DNA sequences. Structural modeling and mutational analyses indicated that adjacent positively charged surfaces from EHF and JUN-FOS use electrostatic repulsion to disfavor simultaneous DNA binding. Conservation of positive residues on the JUN-FOS interface identified E74-like ETS TF 1 (ELF1) as an additional ETS TF exhibiting anticooperative DNA binding with JUN-FOS, and we found that ELF1 is frequently down-regulated in prostate cancer. In summary, divergent electrostatic features of ETS TFs at their JUN-FOS interface enable distinct binding events at ETS-AP1 DNA sites, which may drive specific targeting of ETS TFs to facilitate distinct transcriptional programs.

Recognition mechanism of Wilms' tumour suppressor protein and DNA triplets: insights from molecular dynamics simulation and free energy analysis

The Wilms' tumour suppressor protein (WT1) plays a multifaceted role in human cancer processes. Mutations on its DNA recognition domain could lead to Denys-Drash syndrome, and alternate splicing results in insertion of the tripeptide Lys-Thr-Ser (KTS) between the third and fourth zinc fingers (ZFs), leading to changes in the DNA-binding function. However, detailed recognition mechanisms of the WT1-DNA complex have not been explored. To clarify the mutational effects upon WT1 towards DNA binding at the atomic level, molecular dynamics simulations and the molecular mechanics/Poisson Boltzmann surface area (MM/PBSA) method were employed. The simulation results indicate that mutations in ZF domains (E427Q and Q369H) may weaken the binding affinity, and the statistical analyses of the hydrogen bonds and hydrophobic interactions show that eight residues (Lys351, Arg366, Arg375, Arg376, Lys399, Arg403, Arg424 and Arg430) have a significant influence on recognition and binding to DNA. Insertion of the tripeptide KTS could form an immobilized hydrogen-bonding network with Arg403, affecting the flexibility and angle of the linker between ZF3 and ZF4, thus influencing the recognition between the protein and the DNA triplet at its 5' terminus. These results represent the first step towards a thorough characterization of the WT1 recognition mechanisms, providing a better understanding of the structure-function relationship of WT1 and its mutants.

Evaluation of Sox2 binding affinities for distinct DNA patterns using steered molecular dynamics simulation

Transcription factors (TFs) are gene expression regulators that bind to DNA in a sequence‐specific manner and determine the functional characteristics of the gene. It is worthwhile to study the unique characteristics of such specific TF‐binding pattern in DNA. Sox2 recognizes a 6‐ to 7‐base pair consensus DNA sequence; the central four bases of the binding site are highly conserved, whereas the two to three flanking bases are variable. Here, we attempted to analyze the binding affinity and specificity of the Sox2 protein for distinct DNA sequence patterns via steered molecular dynamics, in which a pulling force is employed to dissociate Sox2 from Sox2–DNA during simulation to study the behavior of a complex under nonequilibrium conditions. The simulation results revealed that the first two stacking bases of the binding pattern have an exclusive impact on the binding affinity, with the corresponding mutant complexes showing greater binding and longer dissociation time than the experimental complexes do. In contrast, mutation of the conserved bases tends to reduce the affinity, and mutation of the complete conserved region disrupts the binding. It might pave the way to identify the most likely binding pattern recognized by Sox2 based on the affinity of each configuration. The α2‐helix of Sox2 was found to be the key player in the Sox2–DNA association. The characterization of Sox2's binding patterns for the target genes in the genome helps in understanding of its regulatory functions.

Predicting Variation of DNA Shape Preferences in Protein-DNA Interaction in Cancer Cells with a New Biophysical Model

DNA shape readout is an important mechanism of transcription factor target site recognition, in addition to the sequence readout. Several machine learning-based models of transcription factor-DNA interactions, considering DNA shape features, have been developed in recent years. Here, we present a new biophysical model of protein-DNA interactions by integrating the DNA shape properties. It is based on the neighbor dinucleotide dependency model BayesPI2, where new parameters are restricted to a subspace spanned by the dinucleotide form of DNA shape features. This allows a biophysical interpretation of the new parameters as a position-dependent preference towards specific DNA shape features. Using the new model, we explore the variation of DNA shape preferences in several transcription factors across various cancer cell lines and cellular conditions. The results reveal that there are DNA shape variations at FOXA1 (Forkhead Box Protein A1) binding sites in steroid-treated MCF7 cells. The new biophysical model is useful for elucidating the finer details of transcription factor-DNA interaction, as well as for predicting cancer mutation effects in the future.

MarR family transcription factors: dynamic variations on a common scaffold

Members of the multiple antibiotic resistance regulator (MarR) family of transcription factors are critical for bacterial cells to respond to chemical signals and to convert such signals into changes in gene activity. Obligate dimers belonging to the winged helix-turn-helix protein family, they are critical for regulation of a variety of functions, including degradation of organic compounds and control of virulence gene expression. The conventional regulatory paradigm is based on a genomic locus in which the gene encoding the MarR protein is divergently oriented from a gene under its control; MarR binding to the intergenic region controls expression of both genes by changing the interaction of RNA polymerase with gene promoters. MarR protein oxidation or binding of a small molecule ligand adversely affects DNA binding, resulting in altered expression of the divergent genes. The generality of this simple paradigm, including the regulation of Escherichia coli MarR by direct binding of antibiotics, has been challenged by reports published in recent years. In addition, structural and biochemical analyses of ligand binding to numerous MarR homologs are converging to identify a shared ligand-binding "hot-spot". This review highlights recent research advances that point to shared features, yet at the same time highlights the remarkable flexibility with which members of this protein family implement responses to inducing signals. A more comprehensive understanding of protein function will pave the way towards the development of both antibacterial agents and biosensors that are based on MarR family proteins.

Assembly of Francisella novicida Cpf1 endonuclease in complex with guide RNA and target DNA

Bacteria and archaea use the CRISPR-Cas system as an adaptive response against infection by foreign nucleic acids. Owing to its remarkable flexibility, this mechanism has been harnessed and adopted as a powerful tool for genome editing. The CRISPR-Cas system includes two classes that are subdivided into six types and 19 subtypes according to conservation of the cas gene and loci organization. Recently, a new protein with endonuclease activity belonging to class 2 type V has been identified. This endonuclease, termed Cpf1, in complex with a single CRISPR RNA (crRNA) is able to recognize and cleave a target DNA preceded by a 5'-TTN-3' protospacer-adjacent motif (PAM) complementary to the RNA guide. To obtain structural insight into the inner workings of Cpf1, the crystallization of an active complex containing the full extent of the crRNA and a 31-nucleotide dsDNA target was attempted. The gene encoding Cpf1 from Francisella novicida was cloned, overexpressed and purified. The crRNA was transcribed and purified in vitro. Finally, the ternary FnCpf1-crRNA-DNA complex was assembled and purified by preparative electrophoresis before crystallization. Crystals belonging to space group C2221, with unit-cell parameters a = 85.2, b = 137.6, c = 320.5 Å, were obtained and subjected to preliminary diffraction experiments.

Structural Basis of MeCP2 Distribution on Non-CpG Methylated and Hydroxymethylated DNA

The Rett-syndrome-associated methyl-CpG-binding protein 2 (MeCP2) selectively binds methylated DNA to regulate transcription during the development of mature neurons. Like other members of the methyl-CpG-binding domain (MBD) family, MeCP2 functions through the recognition of symmetrical 5-methylcytosines in CpG (mCG) dinucleotides. Advances in base-level resolution epigenetic mapping techniques have revealed, however, that MeCP2 can bind asymmetrically methylated and hydroxymethylated CpA dinucleotides and that this alternative binding selectivity modifies gene expression in the developing mammalian brain. The structural determinants of binding to methylated CpA (mCA) and hydroxymethylated DNA have not been previously investigated. Here, we employ isothermal titration calorimetry and NMR spectroscopy to characterize MeCP2 binding to methylated and hydroxymethylated mCG and mCA DNA, examine the effects of Rett-syndrome-associated missense mutations, and make comparisons to the related and evolutionarily most ancient protein, MBD2. These analyses reveal that MeCP2 binds mCA with high affinity in a strand-specific and orientation-dependent manner. In contrast, MBD2 does not show high affinity or methyl-specific binding to mCA. The Rett-associated missense mutations (T158M, R106W, and P101S) destabilize the MeCP2 MBD and disrupt the recognition of mCG and mCA equally. Finally, hydroxymethylation of a high-affinity mCA site does not alter the binding properties, whereas hemi-hydroxylation of the equivalent cytosine in an mCG site decreases affinity and specificity. Based on these findings, we suggest that MeCP2 recognition of methylated/hydroxymethylated CpA dinucleotides functions as an epigenetic switch redistributing MeCP2 among mCG and mCA loci.

Cluster Analysis of p53 Binding Site Sequences Reveals Subsets with Different Functions

p53 is an important regulator of cell cycle arrest, senescence, apoptosis and metabolism, and is frequently mutated in tumors. It functions as a tetramer, where each component dimer binds to a decameric DNA region known as a response element. We identify p53 binding site subtypes and examine the functional and evolutionary properties of these subtypes. We start with over 1700 known binding sites and, with no prior labeling, identify two sets of response elements by unsupervised clustering. When combined, they give rise to three types of p53 binding sites. We find that probabilistic and alignment-based assessments of cross-species conservation show no strong evidence of differential conservation between types of binding sites. In contrast, functional analysis of the genes most proximal to the binding sites provides strong bioinformatic evidence of functional differentiation between the three types of binding sites. Our results are consistent with recent structural data identifying two conformations of the L1 loop in the DNA binding domain, suggesting that they reflect biologically meaningful groups imposed by the p53 protein structure.

Structure of the I-SceI nuclease complexed with its dsDNA target and three catalytic metal ions

Homing endonucleases are highly specific DNA-cleaving enzymes that recognize and cleave long stretches of DNA. The engineering of these enzymes provides instruments for genome modification in a wide range of fields, including gene targeting. The homing endonuclease I-SceI from the yeast Saccharomyces cerevisiae has been purified after overexpression in Escherichia coli and its crystal structure has been determined in complex with its target DNA. In order to evaluate the number of ions that are involved in the cleavage process, thus determining the catalytic mechanism, crystallization experiments were performed in the presence of Mn(2+), yielding crystals that were suitable for X-ray diffraction analysis. The crystals belonged to the orthorhombic space group P212121, with unit-cell parameters a = 80.11, b = 80.57, c = 130.87 Å, α = β = γ = 90°. The self-rotation function and the Matthews coefficient suggested the presence of two protein-DNA complexes in the asymmetric unit. The crystals diffracted to a resolution limit of 2.9 Å using synchrotron radiation. From the anomalous data, it was determined that three cations are involved in catalysis and it was confirmed that I-SceI follows a two-metal-ion DNA-strand cleavage mechanism.

Insights into the DNA-binding mechanism of a LytTR-type transcription regulator

A combination of surface plasmon resonance (SPR) spectroscopy and interaction map® (IM) analysis was used to characterize binding of the LytTR-type response regulator YpdB to promoter DNA. YpdB follows an ‘AB-BA’ mechanism involving sequential and cooperative DNA binding followed by rapid successive promoter clearance.

Most bacterial response regulators (RRs) make contact with DNA through a recognition α-helix in their DNA-binding domains. An emerging class of RRs interacts with DNA via a relatively novel type of binding domain, called the LytTR domain, which is mainly composed of β-strands. YpdB belongs to this latter class, is part of a nutrient-sensing network in Escherichia coli and triggers expression of its only target gene, yhjX, in response to extracellular pyruvate. Expression of yhjX mainly occurs in the late exponential growth phase, and in a pulsed manner. Although the DNA-binding sites for YpdB are well defined, exactly how YpdB initiates pulsed gene expression has remained elusive. To address this question, we measured the binding kinetics of wild-type YpdB and the phosphomimetic variant YpdB-D53E to the yhjX promoter region (P_yhjX) using surface plasmon resonance (SPR) spectroscopy combined with interaction map® (IM) analysis. Both YpdB and YpdB-D53E bound as monomers to the tandem-repeat sequences in the promoter, with YpdB-D53E displaying a higher maximal binding rate than YpdB. Furthermore, we identified a high-affinity (A-site) and a low-affinity binding site (B-site) within the yhjX promoter. Only YpdB-D53E utilizes an ‘AB-BA’ DNA-binding mechanism, involving sequential and cooperative promoter binding, and rapid, successive promoter clearance. We propose that response regulator phosphorylation, in combination with the cycle of cooperative DNA binding and rapid promoter clearance just described, can account for pulsed gene expression.

Functional architecture of the Reb1-Ter complex of Schizosaccharomyces pombe

Reb1 ofSchizosaccharomyces pomberepresents a family of multifunctional proteins that bind to specific terminator sites (Ter) and cause polar termination of transcription catalyzed by RNA polymerase I (pol I) and arrest of replication forks approaching the Ter sites from the opposite direction. However, it remains to be investigated whether the same mechanism causes arrest of both DNA transactions. Here, we present the structure of Reb1 as a complex with a Ter site at a resolution of 2.7 Å. Structure-guided molecular genetic analyses revealed that it has distinct and well-defined DNA binding and transcription termination (TTD) domains. The region of the protein involved in replication termination is distinct from the TTD. Mechanistically, the data support the conclusion that transcription termination is not caused by just high affinity Reb1-Ter protein-DNA interactions. Rather, protein-protein interactions between the TTD with the Rpa12 subunit of RNA pol I seem to be an integral part of the mechanism. This conclusion is further supported by the observation that double mutations in TTD that abolished its interaction with Rpa12 also greatly reduced transcription termination thereby revealing a conduit for functional communications between RNA pol I and the terminator protein.

Molecular stripping in the NF-κB/IκB/DNA genetic regulatory network

Genetic switches based on the [Formula: see text] system are master regulators of an array of cellular responses. Recent kinetic experiments have shown that [Formula: see text] can actively remove NF-κB bound to its genetic sites via a process called "molecular stripping." This allows the [Formula: see text] switch to function under kinetic control rather than the thermodynamic control contemplated in the traditional models of gene switches. Using molecular dynamics simulations of coarse-grained predictive energy landscape models for the constituent proteins by themselves and interacting with the DNA we explore the functional motions of the transcription factor [Formula: see text] and its various binary and ternary complexes with DNA and the inhibitor IκB. These studies show that the function of the [Formula: see text] genetic switch is realized via an allosteric mechanism. Molecular stripping occurs through the activation of a domain twist mode by the binding of [Formula: see text] that occurs through conformational selection. Free energy calculations for DNA binding show that the binding of [Formula: see text] not only results in a significant decrease of the affinity of the transcription factor for the DNA but also kinetically speeds DNA release. Projections of the free energy onto various reaction coordinates reveal the structural details of the stripping pathways.

Molecular determinants of the interactions between proteins and ssDNA

ssDNA binding proteins (SSBs) protect ssDNA from chemical and enzymatic assault that can derail DNA processing machinery. Complexes between SSBs and ssDNA are often highly stable, but predicting their structures is challenging, mostly because of the inherent flexibility of ssDNA and the geometric and energetic complexity of the interfaces that it forms. Here, we report a newly developed coarse-grained model to predict the structure of SSB-ssDNA complexes. The model is successfully applied to predict the binding modes of six SSBs with ssDNA strands of lengths of 6-65 nt. In addition to charge-charge interactions (which are often central to governing protein interactions with nucleic acids by means of electrostatic complementarity), an essential energetic term to predict SSB-ssDNA complexes is the interactions between aromatic residues and DNA bases. For some systems, flexibility is required from not only the ssDNA but also, the SSB to allow it to undergo conformational changes and the penetration of the ssDNA into its binding pocket. The association mechanisms can be quite varied, and in several cases, they involve the ssDNA sliding along the protein surface. The binding mechanism suggests that coarse-grained models are appropriate to study the motion of SSBs along ssDNA, which is expected to be central to the function carried out by the SSBs.

Spec-seq: determining protein-DNA-binding specificity by sequencing

The specificity of protein-DNA interactions can be determined directly by sequencing the bound and unbound fractions in a standard binding reaction. The procedure is easy and inexpensive, and the accuracy can be high for thousands of sequences assayed in parallel. From the measurements, simple models of specificity, such as position weight matrices, can be assessed for their accuracy and more complex models developed if useful. Those may provide more accurate predictions of in vivo binding sites and can help us to understand the details of recognition. As an example, we demonstrate new information gained about the binding of lac repressor. One can apply the same method to combinations of factors that bind simultaneously to a single DNA and determine both the specificity of the individual factors and the cooperativity between them.

High-resolution specificity from DNA sequencing highlights alternative modes of Lac repressor binding

Knowing the specificity of transcription factors is critical to understanding regulatory networks in cells. The lac repressor-operator system has been studied for many years, but not with high-throughput methods capable of determining specificity comprehensively. Details of its binding interaction and its selection of an asymmetric binding site have been controversial. We employed a new method to accurately determine relative binding affinities to thousands of sequences simultaneously, requiring only sequencing of bound and unbound fractions. An analysis of 2560 different DNA sequence variants, including both base changes and variations in operator length, provides a detailed view of lac repressor sequence specificity. We find that the protein can bind with nearly equal affinities to operators of three different lengths, but the sequence preference changes depending on the length, demonstrating alternative modes of interaction between the protein and DNA. The wild-type operator has an odd length, causing the two monomers to bind in alternative modes, making the asymmetric operator the preferred binding site. We tested two other members of the LacI/GalR protein family and find that neither can bind with high affinity to sites with alternative lengths or shows evidence of alternative binding modes. A further comparison with known and predicted motifs suggests that the lac repressor may be unique in this ability and that this may contribute to its selection.

Structure of a new DNA-binding domain which regulates pathogenesis in a wide variety of fungi

WOPR-domain proteins are found throughout the fungal kingdom where they function as master regulators of cell morphology and pathogenesis. Genetic and biochemical experiments previously demonstrated that these proteins bind to specific DNA sequences and thereby regulate transcription. However, their primary sequence showed no relationship to any known DNA-binding domain, and the basis for their ability to recognize DNA sequences remained unknown. Here, we describe the 2.6-Å crystal structure of a WOPR domain in complex with its preferred DNA sequence. The structure reveals that two highly conserved regions, separated by an unconserved linker, form an interdigitated β-sheet that is tilted into the major groove of DNA. Although the main interaction surface is in the major groove, the highest-affinity interactions occur in the minor groove, primarily through a deeply penetrating arginine residue. The structure reveals a new, unanticipated mechanism by which proteins can recognize specific sequences of DNA.

BuD, a helix-loop-helix DNA-binding domain for genome modification

Crystal structures of BurrH and the BurrH–DNA complex are reported.

DNA editing offers new possibilities in synthetic biology and biomedicine for modulation or modification of cellular functions to organisms. However, inaccuracy in this process may lead to genome damage. To address this important problem, a strategy allowing specific gene modification has been achieved through the addition, removal or exchange of DNA sequences using customized proteins and the endogenous DNA-repair machinery. Therefore, the engineering of specific protein–DNA interactions in protein scaffolds is key to providing ‘toolkits’ for precise genome modification or regulation of gene expression. In a search for putative DNA-binding domains, BurrH, a protein that recognizes a 19 bp DNA target, was identified. Here, its apo and DNA-bound crystal structures are reported, revealing a central region containing 19 repeats of a helix–loop–helix modular domain (BurrH domain; BuD), which identifies the DNA target by a single residue-to-nucleotide code, thus facilitating its redesign for gene targeting. New DNA-binding specificities have been engineered in this template, showing that BuD-derived nucleases (BuDNs) induce high levels of gene targeting in a locus of the human haemoglobin β (HBB) gene close to mutations responsible for sickle-cell anaemia. Hence, the unique combination of high efficiency and specificity of the BuD arrays can push forward diverse genome-modification approaches for cell or organism redesign, opening new avenues for gene editing.

Crystallization and preliminary X-ray diffraction analysis of the homing endonuclease I-CvuI from Chlorella vulgaris in complex with its target DNA

Homing endonucleases are highly specific DNA-cleaving enzymes that recognize long stretches of DNA. The engineering of these enzymes provides novel instruments for genome modification in a wide range of fields, including gene targeting, by inducing specific double-strand breaks. I-CvuI is a homing endonuclease from the green alga Chlorella vulgaris. This enzyme was purified after overexpression in Escherichia coli. Crystallization experiments of I-CvuI in complex with its DNA target in the presence of Mg(2+) yielded crystals suitable for X-ray diffraction analysis. The crystals belonged to the orthorhombic space group P212121, with unit-cell parameters a = 62.83, b = 83.56, c = 94.40 Å. The self-rotation function and the Matthews coefficient suggested the presence of one protein-DNA complex per asymmetric unit. The crystals diffracted to a resolution limit of 1.9 Å using synchrotron radiation.

Expression, purification, crystallization and preliminary X-ray diffraction analysis of the novel modular DNA-binding protein BurrH in its apo form and in complex with its target DNA

Different genome-editing strategies have fuelled the development of new DNA-targeting molecular tools allowing precise gene modifications. Here, the expression, purification, crystallization and preliminary X-ray diffraction of BurrH, a novel DNA-binding protein from Burkholderia rhizoxinica, are reported. Crystallization experiments of BurrH in its apo form and in complex with its target DNA yielded crystals suitable for X-ray diffraction analysis. The crystals of the apo form belonged to the primitive hexagonal space group P3(1) or its enantiomorph P3(2), with unit-cell parameters a = b = 73.28, c = 268.02 Å, α = β = 90, γ = 120°. The BurrH-DNA complex crystallized in the monoclinic space group P2(1), with unit-cell parameters a = 70.15, b = 95.83, c = 76.41 Å, α = γ = 90, β = 109.51°. The self-rotation function and the Matthews coefficient suggested the presence of two protein molecules per asymmetric unit in the apo crystals and one protein-DNA complex in the monoclinic crystals. The crystals diffracted to resolution limits of 2.21 and 2.65 Å, respectively, using synchrotron radiation.