Making the bend: DNA tertiary structure and protein-DNA interactions

FRET analysis of the unwrapping of nucleosomal DNA containing a sequence characteristic of the + 1 nucleosome

Sequence-dependent mechanical properties of DNA could play essential roles in nuclear processes by affecting histone-DNA interactions. Previously, we found that the DNA entry site of the first nucleosomes from the transcription start site (+ 1 nucleosome) in budding yeast enriches AA/TT steps, but not the exit site, and the biased presence of AA/TT in the entry site was associated with the transcription levels of yeast genes. Because AA/TT is a rigid dinucleotide step, we considered that AA/TT causes DNA unwrapping. However, our previous MNase-seq experiments with reconstituted nucleosomes left some doubt regarding this interpretation, owing to its high exonuclease activity. Furthermore, MNase cleavage did not provide direct evidence of its structural state. In this study, Förster resonance energy transfer (FRET) measurements were used to investigate salt-induced conformational changes in nucleosomal DNA containing AA/TT repeats at the entry site. We observed that the AA/TT region wrapped around the histone core was as likely as other DNA sequences at physiological salt concentrations. However, it unwrapped at a lower salt concentration, indicating weaker electrostatic interactions with the histone core. Ethidium-induced nucleosome disruption assay showed that the intercalator had greater access to DNA with AA/TT at the entry site. Taken together, these results suggest that AA/TT at the entry sites induces DNA unwrapping from the histone core on the promoter side, which may promote transcriptional activation in response to the approach of transcription-related proteins.

Bacteroides expand the functional versatility of a conserved transcription factor and transcribed DNA to program capsule diversity

The genomes of human gut bacteria in the genus Bacteroides include numerous operons for biosynthesis of diverse capsular polysaccharides (CPSs). The first two genes of each CPS operon encode a locus-specific paralog of transcription elongation factor NusG (called UpxY), which enhances transcript elongation, and a UpxZ protein that inhibits noncognate UpxYs. This process, together with promoter inversions, ensures that a single CPS operon is transcribed in most cells. Here, we use in-vivo nascent-RNA sequencing and promoter-less in-vitro transcription (PIVoT) to show that UpxY recognizes a paused RNA polymerase via sequences in both the exposed non-template DNA and the upstream duplex DNA. UpxY association is aided by 'pause-then-escape' nascent RNA hairpins. UpxZ binds non-cognate UpxYs to directly inhibit UpxY association. This UpxY-UpxZ hierarchical regulatory program allows Bacteroides to generate subpopulations of cells producing diverse CPSs for optimal fitness.

Structural dynamics of the Ca2+-regulated cutinase towards structure-based improvement of PET degradation activity

We previously revealed the structural basis of Ca2+ dependent regulation of a polyethylene terephthalate (PET)-degrading enzyme, Cut190, and proposed a unique reaction cycle in which the enzyme repeatedly binds and releases Ca2+. Here, we report crystal structures of Cut190 mutants with high thermal stability complexed with PET-like ligands that contain aromatic rings. The structural information has allowed us to perform further computational analyses using a PET-trimer bound model. Our multicanonical molecular dynamics simulations and subsequent analyses of the free energy landscapes revealed a novel intermediate form that occurs during the enzymatic reaction cycle. Furthermore, the computational analyses were used to investigate the effect of the point mutations F77L and F81L in the Ca2+-binding site, which showed that the former stabilizes the engaged and open forms to improve transition between the open and active forms, while the latter extremely increases the open form. Subsequent experiments showed that the F77L mutation increased the activity, while the F81L mutation decreased the activity. Our computational analysis has enabled us to explore the dynamics of Cut190 on a completely new level, providing key insights into how the balance between the various conformations influences the reaction cycle and ultimately how to improve the reaction cycle.

From straight to curved: A historical perspective of DNA shape

DNA is the macromolecule responsible for storing the genetic information of a cell and it has intrinsic properties such as deformability, stability and curvature. DNA Curvature plays an important role in gene transcription and, consequently, in the subsequent production of proteins, a fundamental process of cells. With recent advances in bioinformatics and theoretical biology, it became possible to analyze and understand the involvement of DNA Curvature as a discriminatory characteristic of gene-promoting regions. These regions act as sites where RNAp (ribonucleic acid-polymerase) binds to initiate transcription. This review aims to describe the formation of Curvature, as well as highlight its importance in predicting promoters. Furthermore, this article provides the potential of DNA Curvature as a distinguishing feature for promoter prediction tools, as well as outlining the calculation procedures that have been described by other researchers. This work may support further studies directed towards the enhancement of promoter prediction software.

Gene Editing with Artificial DNA Scissors

Artificial metallo-nucleases (AMNs) are small molecule DNA cleavage agents, also known as DNA molecular scissors, and represent an important class of chemotherapeutic with high clinical potential. This review provides a primary level of exploration on the concepts key to this area including an introduction to DNA structure, function, recognition, along with damage and repair mechanisms. Building on this foundation, we describe hybrid molecules where AMNs are covalently attached to directing groups that provide molecular scissors with enhanced or sequence specific DNA damaging capabilities. As this research field continues to evolve, understanding the applications of AMNs along with synthetic conjugation strategies can provide the basis for future innovations, particularly for designing new artificial gene editing systems.

Bacteroides expand the functional versatility of a universal transcription factor and transcribed DNA to program capsule diversity

Human gut Bacteroides species encode numerous (eight or more) tightly regulated capsular polysaccharides (CPS). Specialized paralogs of the universal transcription elongation factor NusG, called UpxY (Y), and an anti-Y UpxZ (Z) are encoded by the first two genes of each CPS operon. The Y-Z regulators combine with promoter inversions to limit CPS transcription to a single operon in most cells. Y enhances transcript elongation whereas Z inhibits noncognate Ys. How Y distinguishes among cognate CPS operons and how Z inhibits only noncognate Ys are unknown. Using in-vivo nascent-RNA sequencing and promoter-less in vitro transcription (PIVoT), we establish that Y recognizes a paused RNA polymerase via sequences in both the exposed non-template DNA and the upstream duplex DNA. Y association is aided by novel 'pause-then-escape' nascent RNA hairpins. Z binds non-cognate Ys to directly inhibit Y association. This Y-Z hierarchical regulatory program allows Bacteroides to create CPS subpopulations for optimal fitness.

Enforcing Local DNA Kinks by Sequence-Selective Trisintercalating Oligopeptides of a Tricationic Porphyrin: A Polarizable Molecular Dynamics Study

Bisacridinyl-bisarginyl porphyrin (BABAP) is a trisintercalating derivative of a tricationic porphyrin, formerly designed and synthesized in order to selectively target and photosensitize the ten-base pair palindromic sequence d(CGGGCGCCCG)2 . We resorted to the previously derived (Far et al., 2004) lowest energy-minimized (EM) structure of the BABAP complex with this sequence as a starting point. We performed polarizable molecular dynamics (MD) on this complex. It showed, over a 150 ns duration, the persistent binding of the Arg side-chain on each BABAP arm to the two G bases upstream from the central porphyrin intercalation site. We subsequently performed progressive shortenings of the connector chain linking the Arg-Gly backbone to the acridine, from n=6 methylenes to 4, followed by removal of the Gly backbone and further connector shortenings, from n=4 to n=1. These resulted into progressive deformations ('kinks') of the DNA backbone. In its most accented kinked structure, the DNA backbone was found to have a close overlap with that of DNA bound to Cre recombinase, with, at the level of one acridine intercalation site, negative roll and positive tilt values consistent with those experimentally found for this DNA at its own kinked dinucleotide sequence. Thus, in addition to their photosensitizing properties, some BABAP derivatives could induce sequence-selective, controlled DNA deformations, which are targets for cleavage by endonucleases or for repair enzymes.

Predictions of DNA mechanical properties at a genomic scale reveal potentially new functional roles of DNA flexibility

Mechanical properties of DNA have been implied to influence many of its biological functions. Recently, a new high-throughput method, called loop-seq, which allows measuring the intrinsic bendability of DNA fragments, has been developed. Using loop-seq data, we created a deep learning model to explore the biological significance of local DNA flexibility in a range of different species from different kingdoms. Consistently, we observed a characteristic and largely dinucleotide-composition-driven change of local flexibility near transcription start sites. In the presence of a TATA-box, a pronounced peak of high flexibility can be observed. Furthermore, depending on the transcription factor investigated, flanking-sequence-dependent DNA flexibility was identified as a potential factor influencing DNA binding. Compared to randomized genomic sequences, depending on species and taxa, actual genomic sequences were observed both with increased and lowered flexibility. Furthermore, in Arabidopsis thaliana, mutation rates, both de novo and fixed, were found to be associated with relatively rigid sequence regions. Our study presents a range of significant correlations between characteristic DNA mechanical properties and genomic features, the significance of which with regard to detailed molecular relevance awaits further theoretical and experimental exploration.

MucR from Sinorhizobium meliloti: New Insights into Its DNA Targets and Its Ability to Oligomerize

Proteins of the MucR/Ros family play a crucial role in bacterial infection or symbiosis with eukaryotic hosts. MucR from Sinorhizobium meliloti plays a regulatory role in establishing symbiosis with the host plant, both dependent and independent of Quorum Sensing. Here, we report the first characterization of MucR isolated from Sinorhizobium meliloti by mass spectrometry and demonstrate that this protein forms higher-order oligomers in its native condition of expression by SEC-MALS. We show that MucR purified from Sinorhizobium meliloti can bind DNA and recognize the region upstream of the ndvA gene in EMSA, revealing that this gene is a direct target of MucR. Although MucR DNA binding activity was already described, a detailed characterization of Sinorhizobium meliloti DNA targets has never been reported. We, thus, analyze sequences recognized by MucR in the rem gene promoter, showing that this protein recognizes AT-rich sequences and does not require a consensus sequence to bind DNA. Furthermore, we investigate the dependence of MucR DNA binding on the length of DNA targets. Taken together, our studies establish MucR from Sinorhizobium meliloti as a member of a new family of Histone-like Nucleoid Structuring (H-NS) proteins, thus explaining the multifaceted role of this protein in many species of alpha-proteobacteria.

The Cytomegalovirus M35 Protein Directly Binds to the Interferon-β Enhancer and Modulates Transcription of Ifnb1 and Other IRF3-Driven Genes

Induction of type I interferon (IFN) gene expression is among the first lines of cellular defense a virus encounters during primary infection. We previously identified the tegument protein M35 of murine cytomegalovirus (MCMV) as an essential antagonist of this antiviral system, showing that M35 interferes with type I IFN induction downstream of pattern-recognition receptor (PRR) activation. Here, we report structural and mechanistic details of M35's function. Determination of M35's crystal structure combined with reverse genetics revealed that homodimerization is a key feature for M35's immunomodulatory activity. In electrophoretic mobility shift assays (EMSAs), purified M35 protein specifically bound to the regulatory DNA element that governs transcription of the first type I IFN gene induced in nonimmune cells, Ifnb1. DNA-binding sites of M35 overlapped with the recognition elements of interferon regulatory factor 3 (IRF3), a key transcription factor activated by PRR signaling. Chromatin immunoprecipitation (ChIP) showed reduced binding of IRF3 to the host Ifnb1 promoter in the presence of M35. We furthermore defined the IRF3-dependent and the type I IFN signaling-responsive genes in murine fibroblasts by RNA sequencing of metabolically labeled transcripts (SLAM-seq) and assessed M35's global effect on gene expression. Stable expression of M35 broadly influenced the transcriptome in untreated cells and specifically downregulated basal expression of IRF3-dependent genes. During MCMV infection, M35 impaired expression of IRF3-responsive genes aside of Ifnb1. Our results suggest that M35-DNA binding directly antagonizes gene induction mediated by IRF3 and impairs the antiviral response more broadly than formerly recognized. IMPORTANCE Replication of the ubiquitous human cytomegalovirus (HCMV) in healthy individuals mostly goes unnoticed but can impair fetal development or cause life-threatening symptoms in immunosuppressed or -deficient patients. Like other herpesviruses, CMV extensively manipulates its hosts and establishes lifelong latent infections. Murine CMV (MCMV) presents an important model system as it allows the study of CMV infection in the host organism. We previously showed that during entry into host cells, MCMV virions release the evolutionary conserved protein M35 protein to immediately dampen the antiviral type I interferon (IFN) response induced by pathogen detection. Here, we show that M35 dimers bind to regulatory DNA elements and interfere with recruitment of interferon regulatory factor 3 (IRF3), a key cellular factor for antiviral gene expression. Thereby, M35 interferes with expression of type I IFNs and other IRF3-dependent genes, reflecting the importance for herpesviruses to avoid IRF3-mediated gene induction.

Extensive Bioinformatics Analyses Reveal a Phylogenetically Conserved Winged Helix (WH) Domain (Zτ) of Topoisomerase IIα, Elucidating Its Very High Affinity for Left-Handed Z-DNA and Suggesting Novel Putative Functions

The dynamic processes operating on genomic DNA, such as gene expression and cellular division, lead inexorably to topological challenges in the form of entanglements, catenanes, knots, "bubbles", R-loops, and other outcomes of supercoiling and helical disruption. The resolution of toxic topological stress is the function attributed to DNA topoisomerases. A prominent example is the negative supercoiling (nsc) trailing processive enzymes such as DNA and RNA polymerases. The multiple equilibrium states that nscDNA can adopt by redistribution of helical twist and writhe include the left-handed double-helical conformation known as Z-DNA. Thirty years ago, one of our labs isolated a protein from Drosophila cells and embryos with a 100-fold greater affinity for Z-DNA than for B-DNA, and identified it as topoisomerase II (gene Top2, orthologous to the human UniProt proteins TOP2A and TOP2B). GTP increased the affinity and selectivity for Z-DNA even further and also led to inhibition of the isomerase enzymatic activity. An allosteric mechanism was proposed, in which topoII acts as a Z-DNA-binding protein (ZBP) to stabilize given states of topological (sub)domains and associated multiprotein complexes. We have now explored this possibility by comprehensive bioinformatic analyses of the available protein sequences of topoII representing organisms covering the whole tree of life. Multiple alignment of these sequences revealed an extremely high level of evolutionary conservation, including a winged-helix protein segment, here denoted as Zτ, constituting the putative structural homolog of Zα, the canonical Z-DNA/Z-RNA binding domain previously identified in the interferon-inducible RNA Adenosine-to-Inosine-editing deaminase, ADAR1p150. In contrast to Zα, which is separate from the protein segment responsible for catalysis, Zτ encompasses the active site tyrosine of topoII; a GTP-binding site and a GxxG sequence motif are in close proximity. Quantitative Zτ-Zα similarity comparisons and molecular docking with interaction scoring further supported the "B-Z-topoII hypothesis" and has led to an expanded mechanism for topoII function incorporating the recognition of Z-DNA segments ("Z-flipons") as an inherent and essential element. We further propose that the two Zτ domains of the topoII homodimer exhibit a single-turnover "conformase" activity on given G(ate) B-DNA segments ("Z-flipins"), inducing their transition to the left-handed Z-conformation. Inasmuch as the topoII-Z-DNA complexes are isomerase inactive, we infer that they fulfill important structural roles in key processes such as mitosis. Topoisomerases are preeminent targets of anti-cancer drug discovery, and we anticipate that detailed elucidation of their structural-functional interactions with Z-DNA and GTP will facilitate the design of novel, more potent and selective anti-cancer chemotherapeutic agents.

A DNA origami-based device for investigating DNA bending proteins by transmission electron microscopy

The DNA origami technique offers precise positioning of nanoscale objects with high accuracy. This has facilitated the development of DNA origami-based functional nanomechanical devices that enable the investigation of DNA-protein interactions at the single particle level. Herein, we used the DNA origami technique to fabricate a nanoscale device for studying DNA bending proteins. For a proof of concept, we used TATA-box binding protein (TBP) to evaluate our approach. Upon binding to the TATA box, TBP causes a bend to DNA of ∼90°. Our device translates this bending into an angular change that is readily observable with a conventional transmission electron microscope (TEM). Furthermore, we investigated the roles of transcription factor II A (TF(II)A) and transcription factor II B (TF(II)B). Our results indicate that TF(II)A introduces additional bending, whereas TF(II)B does not significantly alter the TBP-DNA structure. Our approach can be readily adopted to a wide range of DNA-bending proteins and will aid the development of DNA-origami-based devices tailored for the investigation of DNA-protein interactions.

DeepBend: An interpretable model of DNA bendability

The bendability of genomic DNA impacts chromatin packaging and protein-DNA binding. However, we do not have a comprehensive understanding of the motifs influencing DNA bendability. Recent high-throughput technologies such as Loop-Seq offer an opportunity to address this gap but the lack of accurate and interpretable machine learning models still remains. Here we introduce DeepBend, a convolutional neural network model with convolutions designed to directly capture the motifs underlying DNA bendability and their periodic occurrences or relative arrangements that modulate bendability. DeepBend consistently performs on par with alternative models while giving an extra edge through mechanistic interpretations. Besides confirming the known motifs of DNA bendability, DeepBend also revealed several novel motifs and showed how the spatial patterns of motif occurrences influence bendability. DeepBend's genome-wide prediction of bendability further showed how bendability is linked to chromatin conformation and revealed the motifs controlling the bendability of topologically associated domains and their boundaries.

Pseudomonas-tailed lytic phages: genome mechanical analysis and putative correlation with virion morphogenesis yield

Aim: To unveil a putative correlation between phage genome flexibility and virion morphogenesis yield. Materials & methods: A deeper analysis of the mechanical properties of three Pseudomonas aeruginosa lytic phage genomes was undertaken, together with full genome cyclizability calculations. Results & conclusion: A putative correlation was established among phage genome flexibility, eclipse timeframe and virion particle morphogenesis yield, with a more flexible phage genome leading to a higher burst size and a more rigid phage genome leading to lower burst sizes. The results obtained are highly relevant to understand the influence of the phage genome plasticity on the virion morphogenesis yield inside the infected bacterial host cells and assumes particular relevance in the actual context of bacterial resistance to antibiotics.

Quantitative contribution of the spacer length in the supercoiling-sensitivity of bacterial promoters

DNA supercoiling acts as a global transcriptional regulator in bacteria, but the promoter sequence or structural determinants controlling its effect remain unclear. It was previously proposed to modulate the torsional angle between the -10 and -35 hexamers, and thereby regulate the formation of the closed-complex depending on the length of the 'spacer' between them. Here, we develop a thermodynamic model of this notion based on DNA elasticity, providing quantitative and parameter-free predictions of the relative activation of promoters containing a short versus long spacer when the DNA supercoiling level is varied. The model is tested through an analysis of in vitro and in vivo expression assays of mutant promoters with variable spacer lengths, confirming its accuracy for spacers ranging from 15 to 19 nucleotides, except those of 16 nucleotides where other regulatory mechanisms likely overcome the effect of this specific step. An analysis at the whole-genome scale in Escherichia coli then demonstrates a significant effect of the spacer length on the genomic expression after transient or inheritable superhelical variations, validating the model's predictions. Altogether, this study shows an example of mechanical constraints associated to promoter binding by RNA Polymerase underpinning a basal and global regulatory mechanism.

DNA Bending Force Facilitates Z-DNA Formation under Physiological Salt Conditions

Z-DNA, a noncanonical helical structure of double-stranded DNA (dsDNA), plays pivotal roles in various biological processes, including transcription regulation. Mechanical stresses on dsDNA, such as twisting and stretching, help to form Z-DNA. However, the effect of DNA bending, one of the most common dsDNA deformations, on Z-DNA formation is utterly unknown. Here, we show that DNA bending induces the formation of Z-DNA, that is, more Z-DNA is formed as the bending force becomes stronger. We regulated the bending force on dsDNA by using D-shaped DNA nanostructures. The B-Z transition was observed by single-molecule fluorescence resonance energy transfer. We found that as the bending force became stronger, Z-DNA was formed at lower Mg²⁺ concentrations. When dsDNA contained cytosine methylations, the B-Z transition occurred at 78 mM Mg²⁺ (midpoint) in the absence of the bending force. However, the B-Z transition occurred at a 28-fold lower Mg²⁺ concentration (2.8 mM) in the presence of the bending force. Monte Carlo simulation suggested that the B-Z transition stabilizes the bent form via the formation of the B-Z junction with base extrusion, which effectively releases the bending stress on DNA. Our results clearly show that the bending force facilitates the B-Z transition under physiological salt conditions.

A conformational study of the 10-23 DNAzyme via programmed DNA self-assembly

This communication measures the inter-helical angle of the 10-23 DNAzyme-substrate complex by atomic force microscopy (AFM). Herein, we have devised a strategy to assemble the DNAzyme-substrate complex into a periodic DNA 2D array, which allows reliable study of the conformation of the 10-23 DNAzyme by AFM imaging and fast Fourier transform (FFT). Specifically, the angle between the two flanking helical domains of the catalytic core has been determined via the repeating distance of the 2D array. We expect that the same strategy can generally be applicable for studying other nucleic acid structures.

Systematic mapping of DNAzymes research from 1995 to 2019

DNAzymes (catalytic DNA) have gained significant diagnostic and therapeutic applications with increasing research output over the years. Functional oligonucleotides are used as molecular recognition elements within biosensors for detection of analytes and viral infections such as SARS-CoV-2. DNAzymes are also applied for silencing and regulating cancer specific genes. However, there has not been any report on systematic analysis to track research status, reveal hotspots, and map knowledge in this field. Therefore, in the present study, research articles on DNAzymes from 1995 to 2019 were extracted from Web of Science (SCI-Expanded) after which, 1037 articles were imported into Rstudio (version 3.6.2) and analysed accordingly. The highest number of articles was published in 2019 (n = 138), while the least was in 1995 (n = 1). The articles were published across 216 journals by 2344 authors with 2337 multi-author and 7 single authors. The most prolific authors were Li Y (n = 47), Liu J (n = 46), Wang L (n = 33), Willner I (n = 33) and Zhang L (n = 33). The top three most productive countries were China (n = 2018), USA (n = 447) and Canada (n = 251). The most productive institutions were Hunan University, China (n = 141), University of Illinois, USA (n = 139) and Fuzhou University, China (n = 101). Despite the increasing interest in this field, international collaborations between institutions were very low which requires immediate attention to mitigate challenges such as limited funding, access to facilities, and existing knowledge gap.

Spectroscopic and in vitro investigations of Fe2+/α-Ketoglutarate-dependent enzymes involved in nucleic acid repair and modification

The activation of molecular oxygen for the highly selective functionalization and repair of DNA and RNA nucleobases is achieved by α-ketoglutarate (α-KG)/iron-dependent dioxygenases. Enzymes of special interest are the human homologs AlkBH of Escherichia coli EcAlkB and ten-eleven translocation (TET) enzymes. These enzymes are involved in demethylation or dealkylation of DNA and RNA, although additional physiological functions are continuously being revealed. Given their importance, studying enzyme-substrate interactions, turnover and kinetic parameters is pivotal for the understanding of the mode of action of these enzymes. Diverse analytical methods, including X-ray crystallography, UV/Vis absorption, electron paramagnetic resonance (EPR), circular dichroism (CD) and NMR spectroscopy have been employed to study the changes in the active site and the overall enzyme structure upon substrate, cofactor and inhibitor addition. Several methods are now available to assess activity of these enzymes. By discussing limitations and possibilities of these techniques for EcAlkB, AlkBH and TET we aim to give a comprehensive synopsis from a bioinorganic point of view, addressing researchers from different disciplines working in the highly interdisciplinary and rapidly evolving field of epigenetic processes and DNA/RNA repair and modification.

Characterization and in vitro testing of newly isolated lytic bacteriophages for the biocontrol of Pseudomonas aeruginosa

Aim: Two lytic phages were isolated using P. aeruginosa DSM19880 as host and fully characterized. Materials & methods: Phages were characterized physicochemically, biologically and genomically. Results & conclusion: Host range analysis revealed that the phages also infect some multidrug-resistant (MDR) P. aeruginosa clinical isolates. Increasing MOI from 1 to 1000 significantly increased phage efficiency and retarded bacteria regrowth, but phage ph0034 (reduction of 7.5 log CFU/ml) was more effective than phage ph0031 (reduction of 5.1 log CFU/ml) after 24 h. Both phages belong to Myoviridae family. Genome sequencing of phages ph0031 and ph0034 showed that they do not carry toxin, virulence, antibiotic resistance and integrase genes. The results obtained are highly relevant in the actual context of bacterial resistance to antibiotics.

Kinetics of DNA looping by Anabaena sensory rhodopsin transducer (ASRT) by using DNA cyclization assay

DNA cyclization assay together with single-molecule FRET was employed to monitor protein-mediated bending of a short dsDNA (~ 100 bp). This method provides a simple and easy way to monitor the structural change of DNA in real-time without necessitating prior knowledge of the molecular structures for the optimal dye-labeling. This assay was applied to study how Anabaena sensory rhodopsin transducer (ASRT) facilitates loop formation of DNA as a possible mechanism for gene regulation. The ASRT-induced DNA looping was maximized at 50 mM of Na⁺, while Mg²⁺ also played an essential role in the loop formation.

Bending of Canonical and G/T Mismatched DNAs

Mismatched base pairs alter the flexibility and intrinsic curvature of DNA. The role of such DNA features is not fully understood in the mismatch repair pathway. MutS/DNA complexes exhibit DNA bending, PHE intercalation, and changes of base-pair parameters near the mismatch. Recently, we have shown that base-pair opening in the absence of MutS can discriminate mismatches from canonical base pairs better than DNA bending. However, DNA bending in the absence of MutS was found to be rather challenging to describe correctly. Here, we present a computational study on the DNA bending of canonical and G/T mismatched DNAs. Five types of geometric parameters covering template-based bending toward the experimental DNA structure, global, and local geometry parameters were employed in biased molecular dynamics in the absence of MutS. None of these parameters showed higher discrimination than the base-pair opening. Only roll could induce a sharply localized bending of DNA as observed in the experimental MutS/DNA structure. Further, we demonstrated that the intercalation of benzene mimicking PHE decreases the energetic cost of DNA bending without any effect on mismatch discrimination.

Structures of heat shock factor trimers bound to DNA

Heat shock factor 1 (HSF1) and 2 (HSF2) play distinct but overlapping regulatory roles in maintaining cellular proteostasis or mediating cell differentiation and development. Upon activation, both HSFs trimerize and bind to heat shock elements (HSEs) present in the promoter region of target genes. Despite structural insights gained from recent studies, structures reflecting the physiological architecture of this transcriptional machinery remains to be determined. Here, we present co-crystal structures of human HSF1 and HSF2 trimers bound to DNA, which reveal a triangular arrangement of the three DNA-binding domains (DBDs) with protein-protein interactions largely mediated by the wing domain. Two structural properties, different flexibility of the wing domain and local DNA conformational changes induced by HSF binding, seem likely to contribute to the subtle differential specificity between HSF1 and HSF2. Besides, two more structures showing DBDs bound to "two-site" head-to-head HSEs were determined as additions to the published tail-to-tail dimer-binding structures.

Symphony of the DNA flexibility and sequence environment orchestrates p53 binding to its responsive elements

The p53 tumor suppressor protein maintains the genome fidelity and integrity by modulating several cellular activities. It regulates these events by interacting with a heterogeneous set of response elements (REs) of regulatory genes in the background of chromatin configuration. At the p53-RE interface, both the base readout and torsional-flexibility of DNA account for high-affinity binding. However, DNA structure is an entanglement of a multitude of physicochemical features, both local and global structure should be considered for dealing with DNA-protein interactions. The goal of current research work is to conceptualize and abstract basic principles of p53-RE binding affinity as a function of structural alterations in DNA such as bending, twisting, and stretching flexibility and shape. For this purpose, we have exploited high throughput in-vitro relative affinity information of responsive elements and genome binding events of p53 from HT-Selex and ChIP-Seq experiments respectively. Our results confirm the role of torsional flexibility in p53 binding, and further, we reveal that DNA axial bending, stretching stiffness, propeller twist, and wedge angles are intimately linked to p53 binding affinity when compared to homeodomain, bZIP, and bHLH proteins. Besides, a similar DNA structural environment is observed in the distal sequences encompassing the actual binding sites of p53 cistrome genes. Additionally, we revealed that p53 cistrome target genes have unique promoter architecture, and the DNA flexibility of genomic sequences around REs in cancer and normal cell types display major differences. Altogether, our work provides a keynote on DNA structural features of REs that shape up the in-vitro and in-vivo high-affinity binding of the p53 transcription factor.

The N-terminal Helix-Turn-Helix Motif of Transcription Factors MarA and Rob Drives DNA Recognition

DNA-binding proteins play an important role in gene regulation and cellular function. The transcription factors MarA and Rob are two homologous members of the AraC/XylS family that regulate multidrug resistance. They share a common DNA-binding domain, and Rob possesses an additional C-terminal domain that permits binding of low-molecular weight effectors. Both proteins possess two helix-turn-helix (HTH) motifs capable of binding DNA; however, while MarA interacts with its promoter through both HTH-motifs, prior studies indicate that Rob binding to DNA via a single HTH-motif is sufficient for tight binding. In the present work, we perform microsecond time scale all-atom simulations of the binding of both transcription factors to different DNA sequences to understand the determinants of DNA recognition and binding. Our simulations characterize sequence-dependent changes in dynamical behavior upon DNA binding, showcasing the role of Arg40 of the N-terminal HTH-motif in allowing for specific tight binding. Finally, our simulations demonstrate that an acidic C-terminal loop of Rob can control the DNA binding mode, facilitating interconversion between the distinct DNA binding modes observed in MarA and Rob. In doing so, we provide detailed molecular insight into DNA binding and recognition by these proteins, which in turn is an important step toward the efficient design of antivirulence agents that target these proteins.

The Physics-Biology continuum challenges darwinism: Evolution is directed by the homeostasis-dependent bidirectional relation between genome and phenotype

The physics-biology continuum relies on the fact that life emerged from prebiotic molecules. Here, I argue that life emerged from the coupling between nucleic acid and protein synthesis during which proteins (or proto-phenotypes) maintained the physicochemical parameter equilibria (or proto-homeostasis) in the proximity of their encoding nucleic acids (or proto-genomes). This protected the proto-genome physicochemical integrity (i.e., atomic composition) from environmental physicochemical constraints, and therefore increased the probability of reproducing the proto-genome without variation. From there, genomes evolved depending on the biological activities they generated in response to environmental fluctuations. Thus, a genome maintaining homeostasis (i.e., internal physicochemical parameter equilibria), despite and in response to environmental fluctuations, maintains its physicochemical integrity and has therefore a higher probability to be reproduced without variation. Consequently, descendants have a higher probability to share the same phenotype than their parents. Otherwise, the genome is modified during replication as a consequence of the imbalance of the internal physicochemical parameters it generates, until new mutation-deriving biological activities maintain homeostasis in offspring. In summary, evolution depends on feedforward and feedback loops between genome and phenotype, as the internal physicochemical conditions that a genome generates ─ through its derived phenotype in response to environmental fluctuations ─ in turn either guarantee its stability or direct its variation. Evolution may not be explained by the Darwinism-derived, unidirectional principle (random mutations-phenotypes-natural selection) but rather by the bidirectional relationship between genome and phenotype, in which the phenotype in interaction with the environment directs the evolution of the genome it derives from.

DNA-Targeted Metallodrugs: An Untapped Source of Artificial Gene Editing Technology

DNA binding metal complexes are synonymous with anticancer drug discovery. Given the array of structural and chemical reactivity properties available through careful design, metal complexes have been directed to bind nucleic acid structures through covalent or noncovalent binding modes. Several recognition modes - including crosslinking, intercalation, and oxidation - are central to the clinical success of broad-spectrum anticancer metallodrugs. However, recent progress in nucleic acid click chemistry coupled with advancement in our understanding of metal complex-nucleic acid interactions has opened up new avenues in genetic engineering and targeted therapies. Several of these applications are enabled by the hybridisation of oligonucleotide or polyamine probes to discrete metal complexes, which facilitate site-specific reactivity at the nucleic acid interface under the guidance of the probe. This Review focuses on recent advancements in hybrid design and, by way of an introduction to this topic, we provide a detailed overview of nucleic acid structures and metal complex-nucleic acid interactions. Our aim is to provide readers with an insight on the rational design of metal complexes with DNA recognition properties and an understanding of how the sequence-specific targeting of these interactions can be achieved for gene engineering applications.

Chromatin hierarchical branching visualized at the nanoscale by electron microscopy

Chromatin is spatially organized in a hierarchical manner by virtue of single nucleosomes condensing into higher order chromatin structures, conferring various mechanical properties and biochemical signals. These higher order chromatin structures regulate genomic function by organization of the heterochromatin and euchromatin landscape. Less is known about its transition state from higher order heterochromatin to the lower order nucleosome form, and there is no information on its physical properties. We have developed a facile method of electron microscopy visualization to reveal the interphase chromatin in eukaryotic cells and its organization into hierarchical branching structures. We note that chromatin hierarchical branching can be distinguished at four levels, clearly indicating the stepwise transition from heterochromatin to euchromatin. The protein-DNA density across the chromatin fibers decreases during the transition from compacted heterochromatin to dispersed euchromatin. Moreover, the thickness of the chromatin ranges between 10 to 270 nm, and the controversial 30 nm chromatin fiber exists as a prominent intermediate structure. This study provides important insights into higher order chromatin organization which plays a key role in diseases such as cancer.

Genome-wide analysis of DNA uptake across the outer membrane of naturally competent Haemophilus influenzae

The genomes of naturally competent Pasteurellaceae and Neisseriaceae have many short uptake sequences (USS), which allow them to distinguish self-DNA from foreign DNA. To fully characterize this preference we developed genome-wide maps of DNA uptake using both a sequence-based computational model and genomic DNA that had been sequenced after uptake by and recovery from competent Haemophilus influenzae cells. When DNA fragments were shorter than the average USS spacing of ∼1,000 bp, sharp peaks of uptake were centered at USS and separated by valleys with 1000-fold lower uptake. Long DNA fragments (1.5–17 kb) gave much less variation, with 90% of positions having uptake within 2-fold of the mean. All detectable uptake biases arose from sequences that fit the USS uptake motif. Simulated competition predicted that, in its respiratory tract environment, H. influenzae will efficiently take up its own DNA even when human DNA is present in 100-fold excess.

•

For short DNA fragments, an uptake sequence (USS) improves DNA uptake 1000-fold

•

Most longer H. influenzae fragments have USS, giving even uptake across the genome

•

Preferred USS are stiff, so strand melting may facilitate kinking for uptake

•

H. influenzae will take up its own DNA 100-fold better than human DNA

Orientation and dynamics of Cu2+ based DNA labels from force field parameterized MD elucidates the relationship between EPR distance constraints and DNA backbone distances

Pulsed electron paramagnetic resonance (EPR) based distance measurements using the recently developed Cu2+-DPA label present a promising strategy for measuring DNA backbone distance constraints. Herein we develop force field parameters for Cu2+-DPA in order to understand the features of this label at an atomic level. We perform molecular dynamics (MD) simulations using the force field parameters of Cu2+-DPA on four different DNA duplexes. The distance between the Cu2+ centers, extracted from the 2 μs MD trajectories, agrees well with the experimental distance for all the duplexes. Further analyses of the trajectory provide insight into the orientation of the Cu2+-DPA inside the duplex that leads to such agreement with experiments. The MD results also illustrate the ability of the Cu2+-DPA to report on the DNA backbone distance constraints. Furthermore, measurement of fluctuations of individual residues showed that the flexibility of Cu2+-DPA in a DNA depends on the position of the label in the duplex, and a 2 μs MD simulation is not sufficient to fully capture the experimental distribution in some cases. Finally, the MD trajectories were utilized to understand the key aspects of the double electron electron resonance (DEER) results. The lack of orientational selectivity effects of the Cu2+-DPA at Q-band frequency is rationalized in terms of fluctuations in the Cu2+ coordination environment and rotameric fluctuations of the label linker. Overall, a combination of EPR and MD simulations based on the Cu2+-DPA labelling strategy can contribute towards understanding changes in DNA backbone conformations during protein-DNA interactions.

Quantitative analysis of transcription factor binding and expression using calling cards reporter arrays

We report a tool, Calling Cards Reporter Arrays (CCRA), that measures transcription factor (TF) binding and the consequences on gene expression for hundreds of synthetic promoters in yeast. Using Cbf1p and MAX, we demonstrate that the CCRA method is able to detect small changes in binding free energy with a sensitivity comparable to in vitro methods, enabling the measurement of energy landscapes in vivo. We then demonstrate the quantitative analysis of cooperative interactions by measuring Cbf1p binding at synthetic promoters with multiple sites. We find that the cooperativity between Cbf1p dimers varies sinusoidally with a period of 10.65 bp and energetic cost of 1.37 K_BT for sites that are positioned ‘out of phase’. Finally, we characterize the binding and expression of a group of TFs, Tye7p, Gcr1p and Gcr2p, that act together as a ‘TF collective’, an important but poorly characterized model of TF cooperativity. We demonstrate that Tye7p often binds promoters without its recognition site because it is recruited by other collective members, whereas these other members require their recognition sites, suggesting a hierarchy where these factors recruit Tye7p but not vice versa. Our experiments establish CCRA as a useful tool for quantitative investigations into TF binding and function.

Cu2+-based distance measurements by pulsed EPR provide distance constraints for DNA backbone conformations in solution

Electron paramagnetic resonance (EPR) has become an important tool to probe conformational changes in nucleic acids. An array of EPR labels for nucleic acids are available, but they often come at the cost of long tethers, are dependent on the presence of a particular nucleotide or can be placed only at the termini. Site directed incorporation of Cu2+-chelated to a ligand, 2,2'dipicolylamine (DPA) is potentially an attractive strategy for site-specific, nucleotide independent Cu2+-labelling in DNA. To fully understand the potential of this label, we undertook a systematic and detailed analysis of the Cu2+-DPA motif using EPR and molecular dynamics (MD) simulations. We used continuous wave EPR experiments to characterize Cu2+ binding to DPA as well as optimize Cu2+ loading conditions. We performed double electron-electron resonance (DEER) experiments at two frequencies to elucidate orientational selectivity effects. Furthermore, comparison of DEER and MD simulated distance distributions reveal a remarkable agreement in the most probable distances. The results illustrate the efficacy of the Cu2+-DPA in reporting on DNA backbone conformations for sufficiently long base pair separations. This labelling strategy can serve as an important tool for probing conformational changes in DNA upon interaction with other macromolecules.

Dynamic regulation of Z-DNA in the mouse prefrontal cortex by the RNA-editing enzyme Adar1 is required for fear extinction

DNA forms conformational states beyond the right-handed double-helix; however, the functional relevance of these non-canonical structures in the brain remains unknown. We show that, in the prefrontal cortex of mice, the formation of one such structure, Z-DNA, is involved in the regulation of extinction memory. Z-DNA is formed during fear learning, and reduced during extinction learning, which is mediated, in part, by a direct interaction between Z-DNA and the RNA editing enzyme Adar1. Adar1 binds to Z-DNA during fear extinction learning which leads to a reduction in Z-DNA at sites where Adar1 is recruited. Knockdown of Adar1 leads to an inability to modify a previously acquired fear memory and blocks activity-dependent changes in DNA structure and RNA state; effects that are fully rescued by the introduction of full-length Adar1. These findings suggest a novel mechanism of learning-induced gene regulation dependent on both proteins which recognize DNA structure, and the state.

Specific subfamilies of transposable elements contribute to different domains of T lymphocyte enhancers

Transposable elements (TEs) compose nearly half of mammalian genomes and provide building blocks for cis-regulatory elements. Using high-throughput sequencing, we show that 84 TE subfamilies are overrepresented, and distributed in a lineage-specific fashion in core and boundary domains of CD8⁺ T cell enhancers. Endogenous retroviruses are most significantly enriched in core domains with accessible chromatin, and bear recognition motifs for immune-related transcription factors. In contrast, short interspersed elements (SINEs) are preferentially overrepresented in nucleosome-containing boundaries. A substantial proportion of these SINEs harbor a high density of the enhancer-specific histone mark H3K4me1 and carry sequences that match enhancer boundary nucleotide composition. Motifs with regulatory features are better preserved within enhancer-enriched TE copies compared to their subfamily equivalents located in gene deserts. TE-rich and TE-poor enhancers associate with both shared and unique gene groups and are enriched in overlapping functions related to lymphocyte and leukocyte biology. The majority of T cell enhancers are shared with other immune lineages and are accessible in common hematopoietic progenitors. A higher proportion of immune tissue-specific enhancers are TE-rich compared to enhancers specific to other tissues, correlating with higher TE occurrence in immune gene-associated genomic regions. Our results suggest that during evolution, TEs abundant in these regions and carrying motifs potentially beneficial for enhancer architecture and immune functions were particularly frequently incorporated by evolving enhancers. Their putative selection and regulatory cooption may have accelerated the evolution of immune regulatory networks.

Physicochemical Foundations of Life that Direct Evolution: Chance and Natural Selection are not Evolutionary Driving Forces

The current framework of evolutionary theory postulates that evolution relies on random mutations generating a diversity of phenotypes on which natural selection acts. This framework was established using a top-down approach as it originated from Darwinism, which is based on observations made of complex multicellular organisms and, then, modified to fit a DNA-centric view. In this article, it is argued that based on a bottom-up approach starting from the physicochemical properties of nucleic and amino acid polymers, we should reject the facts that (i) natural selection plays a dominant role in evolution and (ii) the probability of mutations is independent of the generated phenotype. It is shown that the adaptation of a phenotype to an environment does not correspond to organism fitness, but rather corresponds to maintaining the genome stability and integrity. In a stable environment, the phenotype maintains the stability of its originating genome and both (genome and phenotype) are reproduced identically. In an unstable environment (i.e., corresponding to variations in physicochemical parameters above a physiological range), the phenotype no longer maintains the stability of its originating genome, but instead influences its variations. Indeed, environment- and cellular-dependent physicochemical parameters define the probability of mutations in terms of frequency, nature, and location in a genome. Evolution is non-deterministic because it relies on probabilistic physicochemical rules, and evolution is driven by a bidirectional interplay between genome and phenotype in which the phenotype ensures the stability of its originating genome in a cellular and environmental physicochemical parameter-depending manner.

Sublimation thermodynamics of nucleobases derived from fast scanning calorimetry

Fast scanning calorimetry was utilized to measure the sublimation thermodynamics of nucleobases. The results were rationalized at the molecular level.

Signatures of Specific DNA Binding by the AT-Rich Interaction Domain of BAF250a

The AT-rich interaction domain (ARID) containing BAF250a is a subunit of the BAF-A class of SWI/SNF chromatin remodeling complexes. The ARID belongs to a family of conserved DNA binding domains found in several eukaryotic proteins; however, its exact contribution to BAF250a function and the mechanism of its DNA binding are not well understood. Here we have probed the interaction of the BAF250a ARID with three different double-stranded DNA (dsDNA) sequences to understand its DNA binding properties. A comprehensive biophysical and thermodynamic study using nuclear magnetic resonance (NMR) spectroscopy and isothermal titration calorimetry revealed the complex nature of BAF250a ARID-DNA interactions. The thermodynamic signatures of the BAF250a ARID with 12 A-T bp dsDNA (AT-12) are distinct from those of 12 G-C bp dsDNA (GC-12) or 12 bp Dickerson dodecamer DNA (DD-12) sequences. We observed that the binding of the BAF250a ARID with AT-12 DNA is enthalpically driven in a tested temperature range of 5-25 °C. BAF250a ARID/AT-12 DNA interaction exhibited a larger negative calorimetric specific heat change (ΔC_p) compared to that of BAF250a ARID/GC-12 DNA or BAF250a ARID/DD-12 DNA interactions. In the presence of salt (NaCl), ARID/AT-12 DNA binding was less perturbed than ARID/GC-12 DNA or ARID/DD-12 DNA binding. Overall, these results show that BAF250a ARID/AT-12 DNA interaction has signatures of "specific" binding. Furthermore, using NMR chemical shift perturbation experiments, we have identified DNA binding residues on the BAF250a ARID and generated a data-driven HADDOCK model of the ARID/DNA complex that was further supported by mutating key lysine residues that were found to be important for DNA binding.

Parallelized DNA tethered bead measurements to scrutinize DNA mechanical structure

Tethering beads to DNA offers a panel of single molecule techniques for the refined analysis of the conformational dynamics of DNA and the elucidation of the mechanisms of enzyme activity. Recent developments include the massive parallelization of these techniques achieved by the fabrication of dedicated nanoarrays by soft nanolithography. We focus here on two of these techniques: the Tethered Particle motion and Magnetic Tweezers allowing analysis of the behavior of individual DNA molecules in the absence of force and under the application of a force and/or a torque, respectively. We introduce the experimental protocols for the parallelization and discuss the benefits already gained, and to come, for these single molecule investigations.

Passenger Hotspot Mutations in Cancer

Current statistical models for assessing hotspot significance do not properly account for variation in site-specific mutability, thereby yielding many false-positives. We thus (i) detail a Log-normal-Poisson (LNP) background model that accounts for this variability in a manner consistent with models of mutagenesis; (ii) use it to show that passenger hotspots arise from all common mutational processes; and (iii) apply it to a ∼10,000-patient cohort to nominate driver hotspots with far fewer false-positives compared with conventional methods. Overall, we show that many cancer hotspot mutations recurring at the same genomic site across multiple tumors are actually passenger events, recurring at inherently mutable genomic sites under no positive selection.

DNA binding preferences of S. cerevisiae RNA polymerase I Core Factor reveal a preference for the GC-minor groove and a conserved binding mechanism

In Saccharomyces cerevisiae, Core Factor (CF) is a key evolutionarily conserved transcription initiation factor that helps recruit RNA polymerase I (Pol I) to the ribosomal DNA (rDNA) promoter. Upregulated Pol I transcription has been linked to many cancers, and targeting Pol I is an attractive and emerging anti-cancer strategy. Using yeast as a model system, we characterized how CF binds to the Pol I promoter by electrophoretic mobility shift assays (EMSA). Synthetic DNA competitors along with anti-tumor drugs and nucleic acid stains that act as DNA groove blockers were used to discover the binding preference of yeast CF. Our results show that CF employs a unique binding mechanism where it prefers the GC-rich minor groove within the rDNA promoter. In addition, we show that yeast CF is able to bind to the human rDNA promoter sequence that is divergent in DNA sequence and demonstrate CF sensitivity to the human specific Pol I inhibitor, CX-5461. Finally, we show that the human Core Promoter Element (CPE) can functionally replace the yeast Core Element (CE) in vivo when aligned by conserved DNA structural features rather than DNA sequence. Together, these findings suggest that the yeast CF and the human ortholog Selectivity Factor 1 (SL1) use an evolutionarily conserved, structure-based mechanism to target DNA. Their shared mechanism may offer a new avenue in using yeast to explore current and future Pol I anti-cancer compounds.

DNA Supercoiling: an Ancestral Regulator of Gene Expression in Pathogenic Bacteria?

DNA supercoiling acts as a global and ancestral regulator of bacterial gene expression. In this review, we advocate that it plays a pivotal role in host-pathogen interactions by transducing environmental signals to the bacterial chromosome and coordinating its transcriptional response. We present available evidence that DNA supercoiling is modulated by environmental stress conditions relevant to the infection process according to ancestral mechanisms, in zoopathogens as well as phytopathogens. We review the results of transcriptomics studies obtained in widely distant bacterial species, showing that such structural transitions of the chromosome are associated to a complex transcriptional response affecting a large fraction of the genome. Mechanisms and computational models of the transcriptional regulation by DNA supercoiling are then discussed, involving both basal interactions of RNA Polymerase with promoter DNA, and more specific interactions with regulatory proteins. A final part is specifically focused on the regulation of virulence genes within pathogenicity islands of several pathogenic bacterial species.

The C-terminal D/E-rich domain of MBD3 is a putative Z-DNA mimic that competes for Zα DNA-binding activity

The Z-DNA binding domain (Zα), derived from the human RNA editing enzyme ADAR1, can induce and stabilize the Z-DNA conformation. However, the biological function of Zα/Z-DNA remains elusive. Herein, we sought to identify proteins associated with Zα to gain insight into the functional network of Zα/Z-DNA. By pull-down, biophysical and biochemical analyses, we identified a novel Zα-interacting protein, MBD3, and revealed that Zα interacted with its C-terminal acidic region, an aspartate (D)/glutamate (E)-rich domain, with high affinity. The D/E-rich domain of MBD3 may act as a DNA mimic to compete with Z-DNA for binding to Zα. Dimerization of MBD3 via intermolecular interaction of the D/E-rich domain and its N-terminal DNA binding domain, a methyl-CpG-binding domain (MBD), attenuated the high affinity interaction of Zα and the D/E-rich domain. By monitoring the conformation transition of DNA, we found that Zα could compete with the MBD domain for binding to the Z-DNA forming sequence, but not vice versa. Furthermore, co-immunoprecipitation experiments confirmed the interaction of MBD3 and ADAR1 in vivo. Our findings suggest that the interplay of Zα and MBD3 may regulate the transition of the DNA conformation between B- and Z-DNA and thereby modulate chromatin accessibility, resulting in alterations in gene expression.

Determinants of Oligonucleotide Selectivity of APOBEC3B

APOBEC3B (A3B) is a prominent source of mutation in many cancers. To date, it has been difficult to capture the native protein-DNA interactions that confer A3B's substrate specificity by crystallography due to the highly dynamic nature of wild-type A3B active site. We use computational tools to restore a recent crystal structure of a DNA-bound A3B C-terminal domain mutant construct to its wild type sequence, and run molecular dynamics simulations to study its substrate recognition mechanisms. Analysis of these simulations reveal dynamics of the native A3Bctd-oligonucleotide interactions, including the experimentally inaccessible loop 1-oligonucleotide interactions. A second series of simulations in which the target cytosine nucleotide was computationally mutated from a deoxyribose to a ribose show a change in sugar ring pucker, leading to a rearrangement of the binding site and revealing a potential intermediate in the binding pathway. Finally, apo simulations of A3B, starting from the DNA-bound open state, experience a rapid and consistent closure of the binding site, reaching conformations incompatible with substrate binding. This study reveals a more realistic and dynamic view of the wild type A3B binding site and provides novel insights for structure-guided design efforts for A3B.

Identification of GCC-box and TCC-box motifs in the promoters of differentially expressed genes in rice (Oryza sativa L.): Experimental and computational approaches

The transcription factor selectively binds with the cis-regulatory elements of the promoter and regulates the differential expression of genes. In this study, we aimed to identify and validate the presence of GCC-box and TCC-box motifs in the promoters of upregulated differentially expressed genes (UR-DEGs) and downregulated differentially expressed genes (DR-DEGs) under anoxia using molecular beacon probe (MBP) based real-time PCR. The GCC-box motif was detected in UR-DEGs (DnaJ and 60S ribosomal protein L7 genes), whereas, the TCC-box was detected in DR-DEGs (DnaK and CPuORF11 genes). In addition, the mechanism of interaction of AP2/EREBP family transcription factor (LOC_Os03g22170) with GCC-box promoter motif present in DnaJ gene (LOC_Os06g09560) and 60S ribosomal protein L7 gene (LOC_Os08g42920); and TCC-box promoter motif of DnaK gene (LOC_Os02g48110) and CPuORF11 gene (LOC_Os02g01240) were explored using molecular dynamics (MD) simulations analysis including binding free energy calculations, principal component analyses, and free energy landscapes. The binding free energy analysis revealed that AP2/EREBP model residues such as Arg68, Arg72, Arg83, Lys87, and Arg90 were commonly involved in the formation of hydrogen bonds with GCC and TCC-box promoter motifs, suggesting that these residues are critical for strong interaction. The movement of the entire protein bound to DNA was restricted, confirming the stability of the complex. This study provides comprehensive binding information and a more detailed view of the dynamic interaction between proteins and DNA.

Long Time-Scale Atomistic Simulations of the Structure and Dynamics of Transcription Factor-DNA Recognition

Recent years have witnessed an explosion of interest in computational studies of DNA binding proteins, including both coarse-grained and atomistic simulations of transcription factor-DNA recognition, to understand how these transcription factors recognize their binding sites on the DNA with such exquisite specificity. The present study performs microsecond time scale all-atom simulations of the dimeric form of the lactose repressor (LacI), both in the absence of any DNA and in the presence of both specific and nonspecific complexes, considering three different DNA sequences. We examine, specifically, the conformational differences between specific and nonspecific protein-DNA interactions, as well as the behavior of the helix-turn-helix motif of LacI when interacting with the DNA. Our simulations suggest that stable LacI binding occurs primarily to bent A-form DNA, with a loss of LacI conformational entropy and optimization of correlated conformational equilibria across the protein. In addition, binding to the specific operator sequence involves a slightly larger number of stabilizing DNA-protein hydrogen bonds (in comparison to nonspecific complexes), which may account for the experimentally observed specificity for this operator. In doing so, our simulations provide a detailed atomistic description of potential structural drivers for LacI selectivity.

A Comprehensive Biophysical Analysis of the Effect of DNA Binding Drugs on Protamine-induced DNA Condensation

DNA condensation is a ubiquitous phenomenon in biology, yet the physical basis for it has remained elusive. Here, we have explored the mechanism of DNA condensation through the protamine-DNA interaction, and by examining on it the influence of DNA binding drugs. We observed that the DNA condensation is accompanied by B to Ψ-DNA transition as a result of DNA base pair distortions due to protamine binding, bringing about the formation of toroidal structure through coil-globule transition. The binding energetics suggested that electrostatic energy, bending energy and hydration energy must play crucial roles in DNA condensation. EtBr intercalation interferes with the protamine-DNA interaction, challenging the distortion of the DNA helix and separation of DNA base pairs by protamine. Thus, EtBr, by competing directly with protamine, resists the phenomenon of DNA condensation. On the contrary, netropsin impedes the DNA condensation by an allosteric mechanism, by resisting the probable DNA major groove bending by protamine. In summary, we demonstrate that drugs with distinct binding modes use different mechanism to interfere with DNA condensation.