The structural basis of DNA flexibility

Cooperative Gsx2-DNA binding requires DNA bending and a novel Gsx2 homeodomain interface

The conserved Gsx homeodomain (HD) transcription factors specify neural cell fates in animals from flies to mammals. Like many HD proteins, Gsx factors bind A/T-rich DNA sequences prompting the following question: How do HD factors that bind similar DNA sequences in vitro regulate specific target genes in vivo? Prior studies revealed that Gsx factors bind DNA both as a monomer on individual A/T-rich sites and as a cooperative homodimer to two sites spaced precisely 7 bp apart. However, the mechanistic basis for Gsx-DNA binding and cooperativity is poorly understood. Here, we used biochemical, biophysical, structural and modeling approaches to (i) show that Gsx factors are monomers in solution and require DNA for cooperative complex formation, (ii) define the affinity and thermodynamic binding parameters of Gsx2/DNA interactions, (iii) solve a high-resolution monomer/DNA structure that reveals that Gsx2 induces a 20° bend in DNA, (iv) identify a Gsx2 protein-protein interface required for cooperative DNA binding and (v) determine that flexible spacer DNA sequences enhance Gsx2 cooperativity on dimer sites. Altogether, our results provide a mechanistic basis for understanding the protein and DNA structural determinants that underlie cooperative DNA binding by Gsx factors.

Dynamics in Cre-loxP site-specific recombination

Cre recombinase is a phage-derived enzyme that has found utility for precise manipulation of DNA sequences. Cre recognizes and recombines pairs of loxP sequences characterized by an inverted repeat and asymmetric spacer. Cre cleaves and religates its DNA targets such that error-prone repair pathways are not required to generate intact DNA products. Major obstacles to broader applications are lack of knowledge of how Cre recognizes its targets, and how its activity is controlled. The picture emerging from high resolution methods is that the dynamic properties of both the enzyme and its DNA target are important determinants of its activity in both sequence recognition and DNA cleavage. Improved understanding of the role of dynamics in the key steps along the pathway of Cre-loxP recombination should significantly advance our ability to both redirect Cre to new sequences and to control its DNA cleavage activity in the test tube and in cells.

Molecular Characterization and Genome Mechanical Features of Two Newly Isolated Polyvalent Bacteriophages Infecting Pseudomonas syringae pv. garcae

Coffee plants have been targeted by a devastating bacterial disease, a condition known as bacterial blight, caused by the phytopathogen Pseudomonas syringae pv. garcae (Psg). Conventional treatments of coffee plantations affected by the disease involve frequent spraying with copper- and kasugamycin-derived compounds, but they are both highly toxic to the environment and stimulate the appearance of bacterial resistance. Herein, we report the molecular characterization and mechanical features of the genome of two newly isolated (putative polyvalent) lytic phages for Psg. The isolated phages belong to class Caudoviricetes and present a myovirus-like morphotype belonging to the genuses Tequatrovirus (PsgM02F) and Phapecoctavirus (PsgM04F) of the subfamilies Straboviridae (PsgM02F) and Stephanstirmvirinae (PsgM04F), according to recent bacterial viruses' taxonomy, based on their complete genome sequences. The 165,282 bp (PsgM02F) and 151,205 bp (PsgM04F) genomes do not feature any lysogenic-related (integrase) genes and, hence, can safely be assumed to follow a lytic lifestyle. While phage PsgM02F produced a morphogenesis yield of 124 virions per host cell, phage PsgM04F produced only 12 virions per host cell, indicating that they replicate well in Psg with a 50 min latency period. Genome mechanical analyses established a relationship between genome bendability and virion morphogenesis yield within infected host cells.

Cooperative Gsx2-DNA Binding Requires DNA Bending and a Novel Gsx2 Homeodomain Interface

The conserved Gsx homeodomain (HD) transcription factors specify neural cell fates in animals from flies to mammals. Like many HD proteins, Gsx factors bind A/T-rich DNA sequences prompting the question - how do HD factors that bind similar DNA sequences in vitro regulate specific target genes in vivo? Prior studies revealed that Gsx factors bind DNA both as a monomer on individual A/T-rich sites and as a cooperative homodimer to two sites spaced precisely seven base pairs apart. However, the mechanistic basis for Gsx DNA binding and cooperativity are poorly understood. Here, we used biochemical, biophysical, structural, and modeling approaches to (1) show that Gsx factors are monomers in solution and require DNA for cooperative complex formation; (2) define the affinity and thermodynamic binding parameters of Gsx2/DNA interactions; (3) solve a high-resolution monomer/DNA structure that reveals Gsx2 induces a 20° bend in DNA; (4) identify a Gsx2 protein-protein interface required for cooperative DNA binding; and (5) determine that flexible spacer DNA sequences enhance Gsx2 cooperativity on dimer sites. Altogether, our results provide a mechanistic basis for understanding the protein and DNA structural determinants that underlie cooperative DNA binding by Gsx factors, thereby providing a deeper understanding of HD specificity.

DNA conformational dynamics in the context-dependent non-CG CHH methylation by plant methyltransferase DRM2

DNA methylation provides an important epigenetic mechanism that critically regulates gene expression, genome imprinting, and retrotransposon silencing. In plants, DNA methylation is prevalent not only in a CG dinucleotide context but also in non-CG contexts, namely CHG and CHH (H = C, T, or A) methylation. It has been established that plant non-CG DNA methylation is highly context dependent, with the +1- and +2-flanking sequences enriched with A/T nucleotides. How DNA sequence, conformation, and dynamics influence non-CG methylation remains elusive. Here, we report structural and biochemical characterizations of the intrinsic substrate preference of DOMAINS REARRANGED METHYLTRANSFERASE 2 (DRM2), a plant DNA methyltransferase responsible for establishing all cytosine methylation and maintaining CHH methylation. Among nine CHH motifs, the DRM2 methyltransferase (MTase) domain shows marked substrate preference toward CWW (W = A or T) motifs, correlating well with their relative abundance in planta. Furthermore, we report the crystal structure of DRM2 MTase in complex with a DNA duplex containing a flexible TpA base step at the +1/+2-flanking sites of the target nucleotide. Comparative structural analysis of the DRM2-DNA complexes provides a mechanism by which flanking nucleotide composition impacts DRM2-mediated DNA methylation. Furthermore, the flexibility of the TpA step gives rise to two alternative DNA conformations, resulting in different interactions with DRM2 and consequently temperature-dependent shift of the substrate preference of DRM2. Together, this study provides insights into how the interplay between the conformational dynamics of DNA and temperature as an environmental factor contributes to the context-dependent CHH methylation by DRM2.

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.

Flexibility of flanking DNA is a key determinant of transcription factor affinity for the core motif

Selective gene regulation is mediated by recognition of specific DNA sequences by transcription factors (TFs). The extremely challenging task of searching out specific cognate DNA binding sites among several million putative sites within the eukaryotic genome is achieved by complex molecular recognition mechanisms. Elements of this recognition code include the core binding sequence, the flanking sequence context, and the shape and conformational flexibility of the composite binding site. To unravel the extent to which DNA flexibility modulates TF binding, in this study, we employed experimentally guided molecular dynamics simulations of ternary complex of closely related Hox heterodimers Exd-Ubx and Exd-Scr with DNA. Results demonstrate that flexibility signatures embedded in the flanking sequences impact TF binding at the cognate binding site. A DNA sequence has intrinsic shape and flexibility features. While shape features are localized, our analyses reveal that flexibility features of the flanking sequences percolate several basepairs and allosterically modulate TF binding at the core. We also show that lack of flexibility in the motif context can render the cognate site resistant to protein-induced shape changes and subsequently lower TF binding affinity. Overall, this study suggests that flexibility-guided DNA shape, and not merely the static shape, is a key unexplored component of the complex DNA-TF recognition code.

Evolutionary Invariant of the Structure of DNA Double Helix in RNAP II Core Promoters

Eukaryotic and archaeal RNA polymerase II (POL II) machinery is highly conserved, regardless of the extreme changes in promoter sequences in different organisms. The goal of our work is to find the cause of this conservatism. The representative sets of aligned promoter sequences of fifteen organisms belonging to different evolutional stages were studied. Their textual profiles, as well as profiles of the indexes that characterize the secondary structure and the mechanical and physicochemical properties, were analyzed. The evolutionarily stable, extremely heterogeneous special secondary structure of POL II core promoters was revealed, which includes two singular regions—hexanucleotide “INR” around TSS and octanucleotide “TATA element” of about −28 bp upstream. Such structures may have developed at some stage of evolution. It turned out to be so well matched for the pre-initiation complex formation and the subsequent initiation of transcription for POL II machinery that in the course of evolution there were selected only those nucleotide sequences that were able to reproduce these structural properties. The individual features of specific sequences representing the singular region of the promoter of each gene can affect the kinetics of DNA-protein complex formation and facilitate strand separation in double-stranded DNA at the TSS position.

Structural determinants of DNA recognition by the NO sensor NsrR and related Rrf2-type [FeS]-transcription factors

Several transcription factors of the Rrf2 family use an iron-sulfur cluster to regulate DNA binding through effectors such as nitric oxide (NO), cellular redox status and iron levels. [4Fe-4S]-NsrR from Streptomyces coelicolor (ScNsrR) modulates expression of three different genes via reaction and complex formation with variable amounts of NO, which results in detoxification of this gas. Here, we report the crystal structure of ScNsrR complexed with an hmpA1 gene operator fragment and compare it with those previously reported for [2Fe-2S]-RsrR/rsrR and apo-IscR/hyA complexes. Important structural differences reside in the variation of the DNA minor and major groove widths. In addition, different DNA curvatures and different interactions with the protein sensors are observed. We also report studies of NsrR binding to four hmpA1 variants, which indicate that flexibility in the central region is not a key binding determinant. Our study explores the promotor binding specificities of three closely related transcriptional regulators.

The crystal structure of the iron-sulfur protein NsrR from Streptomyces coelicolor bound to a gene operator fragment is reported and compared with other structures, giving insight into the structural determinants of DNA recognition by the NO sensor.

Superhelical DNA liquid crystals from dendrimer-induced DNA compaction

Electrostatic compaction of double stranded DNA induced by a positively charged poly(amidoamine) (PAMAM) dendrimer of generation four (G4) was found to produce two unique types of DNA mesophases, in which the DNA bent into superhelices packed in a tetragonal or hexagonal lattice. The structure formed at a lower dendrimer charge density was three-dimensionally (3D) ordered, as characterized by the P4₁2₁2 space group with a 4₁ screw axis in a tetragonal arrangement, showing that the weakly bent DNA superhelices with a pitch length of ca. 5.0 nm possessed both identical handedness and phase conservation. The 3D ordered structure transformed into a 2D mesophase at a higher dendrimer charge density, wherein the strongly bent superhelices with a pitch length of ca. 4.0 nm organized in a hexagonal lattice without lateral coherence of helical trajectory. The counterion valency of the protonic acid that is used to charge the dendrimer was found to influence the phase diagram. Under a given dendrimer charge density, the complex with a multivalent acid-protonated dendrimer tended to form structures with less curved DNA, attesting that the driving force of charge matching was reduced by increasing the counterion valency of the dendrimer.

A molecular view of DNA flexibility

DNA dynamics can only be understood by taking into account its complex mechanical behavior at different length scales. At the micrometer level, the mechanical properties of single DNA molecules have been well-characterized by polymer models and are commonly quantified by a persistence length of 50 nm (~150 bp). However, at the base pair level (~3.4 Å), the dynamics of DNA involves complex molecular mechanisms that are still being deciphered. Here, we review recent single-molecule experiments and molecular dynamics simulations that are providing novel insights into DNA mechanics from such a molecular perspective. We first discuss recent findings on sequence-dependent DNA mechanical properties, including sequences that resist mechanical stress and sequences that can accommodate strong deformations. We then comment on the intricate effects of cytosine methylation and DNA mismatches on DNA mechanics. Finally, we review recently reported differences in the mechanical properties of DNA and double-stranded RNA, the other double-helical carrier of genetic information. A thorough examination of the recent single-molecule literature permits establishing a set of general 'rules' that reasonably explain the mechanics of nucleic acids at the base pair level. These simple rules offer an improved description of certain biological systems and might serve as valuable guidelines for future design of DNA and RNA nanostructures.

Xenogeneic Silencing and Bacterial Genome Evolution: Mechanisms for DNA Recognition Imply Multifaceted Roles of Xenogeneic Silencers

Horizontal gene transfer (HGT) is a major driving force for bacterial evolution. To avoid the deleterious effects due to the unregulated expression of newly acquired foreign genes, bacteria have evolved specific proteins named xenogeneic silencers to recognize foreign DNA sequences and suppress their transcription. As there is considerable diversity in genomic base compositions among bacteria, how xenogeneic silencers distinguish self- from nonself DNA in different bacteria remains poorly understood. This review summarizes the progress in studying the DNA binding preferences and the underlying molecular mechanisms of known xenogeneic silencer families, represented by H-NS of Escherichia coli, Lsr2 of Mycobacterium, MvaT of Pseudomonas, and Rok of Bacillus. Comparative analyses of the published data indicate that the differences in DNA recognition mechanisms enable these xenogeneic silencers to have clear characteristics in DNA sequence preferences, which are further correlated with different host genomic features. These correlations provide insights into the mechanisms of how these xenogeneic silencers selectively target foreign DNA in different genomic backgrounds. Furthermore, it is revealed that the genomic AT contents of bacterial species with the same xenogeneic silencer family proteins are distributed in a limited range and are generally lower than those species without any known xenogeneic silencers in the same phylum/class/genus, indicating that xenogeneic silencers have multifaceted roles on bacterial genome evolution. In addition to regulating horizontal gene transfer, xenogeneic silencers also act as a selective force against the GC to AT mutational bias found in bacterial genomes and help the host genomic AT contents maintained at relatively low levels.

A perspective on the molecular simulation of DNA from structural and functional aspects

As genetic material, DNA not only carries genetic information by sequence, but also affects biological functions ranging from base modification to replication, transcription and gene regulation through its structural and dynamic properties and variations. The motion and structural properties of DNA involved in related biological processes are also multi-scale, ranging from single base flipping to local DNA deformation, TF binding, G-quadruplex and i-motif formation, TAD establishment, compartmentalization and even chromosome territory formation, just to name a few. The sequence-dependent physical properties of DNA play vital role in all these events, and thus it is interesting to examine how simple sequence information affects DNA and the formation of the chromatin structure in these different hierarchical orders. Accordingly, molecular simulations can provide atomistic details of interactions and conformational dynamics involved in different biological processes of DNA, including those inaccessible by current experimental methods. In this perspective, which is mainly based on our recent studies, we provide a brief overview of the atomistic simulations on how the hierarchical structure and dynamics of DNA can be influenced by its sequences, base modifications, environmental factors and protein binding in the context of the protein-DNA interactions, gene regulation and structural organization of chromatin. We try to connect the DNA sequence, the hierarchical structures of DNA and gene regulation.

GapR binds DNA through dynamic opening of its tetrameric interface

GapR is a nucleoid-associated protein that is an essential regulator of chromosome replication in the cell cycle model Caulobacter crescentus. Here, we demonstrate that free GapR is a homotetramer, but not a dimer as previously reported (Guo et al., Cell 175: 583–597, 2018). We have determined the crystal structure of GapR in complex with a 10-bp A-tract DNA, which has an open tetrameric conformation, different from the closed clamp conformation in the previously reported crystal structure of GapR/DNA complex. The free GapR adopts multiple conformations in dynamic exchange equilibrium, with the major conformation resembling the closed tetrameric conformation, while the open tetrameric conformation is a representative of minor conformers. As it is impossible for the circular genomic DNA to get into the central DNA binding tunnel of the major conformation, we propose that GapR initially binds DNA through the open conformation, and then undergoes structural rearrangement to form the closed conformation which fully encircles the DNA. GapR prefers to bind DNA with 10-bp consecutive A/T base pairs nonselectively (K_d ∼12 nM), while it can also bind GC-rich DNA sequence with a reasonable affinity of about 120 nM. Besides, our results suggest that GapR binding results in widening the minor groove of AT-rich DNA, instead of overtwisting DNA.

Best available technologies and treatment trains to address current challenges in urban wastewater reuse for irrigation of crops in EU countries

Conventional urban wastewater treatment plants (UWTPs) are poorly effective in the removal of most contaminants of emerging concern (CECs), including antibiotics, antibiotic resistant bacteria and antibiotic resistance genes (ARB&ARGs). These contaminants result in some concern for the environment and human health, in particular if UWTPs effluents are reused for crop irrigation. Recently, stakeholders' interest further increased in Europe, because the European Commission is currently developing a regulation on water reuse. Likely, conventional UWTPs will require additional advanced treatment steps to meet water quality limits yet to be officially established for wastewater reuse. Even though it seems that CECs will not be included in the proposed regulation, the aim of this paper is to provide a technical contribution to this discussion as well as to support stakeholders by recommending possible advanced treatment options, in particular with regard to the removal of CECs and ARB&ARGs. Taking into account the current knowledge and the precautionary principle, any new or revised water-related Directive should address such contaminants. Hence, this review paper gathers the efforts of a group of international experts, members of the NEREUS COST Action ES1403, who for three years have been constructively discussing the efficiency of the best available technologies (BATs) for urban wastewater treatment to abate CECs and ARB&ARGs. In particular, ozonation, activated carbon adsorption, chemical disinfectants, UV radiation, advanced oxidation processes (AOPs) and membrane filtration are discussed with regard to their capability to effectively remove CECs and ARB&ARGs, as well as their advantages and drawbacks. Moreover, a comparison among the above-mentioned processes is performed for CECs relevant for crop uptake. Finally, possible treatment trains including the above-discussed BATs are discussed, issuing end-use specific recommendations which will be useful to UWTPs managers to select the most suitable options to be implemented at their own facilities to successfully address wastewater reuse challenges.

DNA Structural Changes Induced by Intermolecular Triple Helix Formation

DNase I footprints of intermolecular DNA triplexes are often accompanied by enhanced cleavage at the 3'-end of the target site at the triplex-duplex junction. We have systematically studied the sequence dependence of this effect by examining oligonucleotide binding to sites flanked by each base in turn. For complexes with a terminal T.AT triplet, the greatest enhancement is seen with ApC, followed by ApG and ApT, with the weakest enhancement at ApA. Similar DNase I enhancements were observed for a triplex with a terminal C+.GC triplet, though with little difference between the different GpN sites. Enhanced reactivity to diethylpyrocarbonate was observed at As that flank the triplex-duplex junction at AAA or AAC but not AAG or AAT. Fluorescence melting experiments demonstrated that the flanking base affected the stability with a 4 °C difference in T m between a flanking C and G. Sequences that produced the strongest enhancement correlated with those having the lower thermal stability. These results are interpreted in terms of oligonucleotide-induced changes in DNA structure and/or flexibility.

DNA Modification Readers and Writers and Their Interplay

Genomic DNA is modified in a postreplicative manner and several modifications, the enzymes responsible for their deposition as well as proteins that read these modifications, have been described. Here, we focus on the impact of DNA modifications on the DNA helix and review the writers and readers of cytosine modifications and how they interplay to shape genome composition, stability, and function.

Nanomechanics of Diaminopurine-Substituted DNA

2,6-diaminopurine (DAP) is a nucleobase analog of adenine. When incorporated into double-stranded DNA (dsDNA), it forms three hydrogen bonds with thymine. Rare in nature, DAP substitution alters the physical characteristics of a DNA molecule without sacrificing sequence specificity. Here, we show that in addition to stabilizing double-strand hybridization, DAP substitution also changes the mechanical and conformational properties of dsDNA. Thermal melting experiments reveal that DAP substitution raises melting temperatures without diminishing sequence-dependent effects. Using a combination of atomic force microscopy (AFM), magnetic tweezer (MT) nanomechanical assays, and circular dichroism spectroscopy, we demonstrate that DAP substitution increases the flexural rigidity of dsDNA yet also facilitates conformational shifts, which manifest as changes in molecule length. DAP substitution increases both the static and dynamic persistence length of DNA (measured by AFM and MT, respectively). In the static case (AFM), in which tension is not applied to the molecule, the contour length of DAP-DNA appears shorter than wild-type (WT)-DNA; under tension (MT), they have similar dynamic contour lengths. At tensions above 60 pN, WT-DNA undergoes characteristic overstretching because of strand separation (tension-induced melting) and spontaneous adoption of a conformation termed S-DNA. Cyclic overstretching and relaxation of WT-DNA at near-zero loading rates typically yields hysteresis, indicative of tension-induced melting; conversely, cyclic stretching of DAP-DNA showed little or no hysteresis, consistent with the adoption of the S-form, similar to what has been reported for GC-rich sequences. However, DAP-DNA overstretching is distinct from GC-rich overstretching in that it happens at a significantly lower tension. In physiological salt conditions, evenly mixed AT/GC DNA typically overstretches around 60 pN. GC-rich sequences overstretch at similar if not slightly higher tensions. Here, we show that DAP-DNA overstretches at 52 pN. In summary, DAP substitution decreases the overall stability of the B-form double helix, biasing toward non-B-form DNA helix conformations at zero tension and facilitating the B-to-S transition at high tension.

Flexibility and structure of flanking DNA impact transcription factor affinity for its core motif

Spatial and temporal expression of genes is essential for maintaining phenotype integrity. Transcription factors (TFs) modulate expression patterns by binding to specific DNA sequences in the genome. Along with the core binding motif, the flanking sequence context can play a role in DNA-TF recognition. Here, we employ high-throughput in vitro and in silico analyses to understand the influence of sequences flanking the cognate sites in binding of three most prevalent eukaryotic TF families (zinc finger, homeodomain and bZIP). In vitro binding preferences of each TF toward the entire DNA sequence space were correlated with a wide range of DNA structural parameters, including DNA flexibility. Results demonstrate that conformational plasticity of flanking regions modulates binding affinity of certain TF families. DNA duplex stability and minor groove width also play an important role in DNA-TF recognition but differ in how exactly they influence the binding in each specific case. Our analyses further reveal that the structural features of preferred flanking sequences are not universal, as similar DNA-binding folds can employ distinct DNA recognition modes.

How bacterial xenogeneic silencer rok distinguishes foreign from self DNA in its resident genome

Bacterial xenogeneic silencers play important roles in bacterial evolution by recognizing and inhibiting expression from foreign genes acquired through horizontal gene transfer, thereby buffering against potential fitness consequences of their misregulated expression. Here, the detailed DNA binding properties of Rok, a xenogeneic silencer in Bacillus subtilis, was studied using protein binding microarray, and the solution structure of its C-terminal DNA binding domain was determined in complex with DNA. The C-terminal domain of Rok adopts a typical winged helix fold, with a novel DNA recognition mechanism different from other winged helix proteins or xenogeneic silencers. Rok binds the DNA minor groove by forming hydrogen bonds to bases through N154, T156 at the N-terminal of α3 helix and R174 of wing W1, assisted by four lysine residues interacting electrostatically with DNA backbone phosphate groups. These structural features endow Rok with preference towards DNA sequences harboring AACTA, TACTA, and flexible multiple TpA steps, while rigid A-tracts are disfavored. Correspondingly, the Bacillus genomes containing Rok are rich in A-tracts and show a dramatic underrepresentation of AACTA and TACTA, which are significantly enriched in Rok binding regions. These observations suggest that the xenogeneic silencing protein and its resident genome may have evolved cooperatively.

The Interaction of the Metallo-Glycopeptide Anti-Tumour Drug Bleomycin with DNA

The cancer chemotherapeutic drug, bleomycin, is clinically used to treat several neoplasms including testicular and ovarian cancers. Bleomycin is a metallo-glycopeptide antibiotic that requires a transition metal ion, usually Fe(II), for activity. In this review, the properties of bleomycin are examined, especially the interaction of bleomycin with DNA. A Fe(II)-bleomycin complex is capable of DNA cleavage and this process is thought to be the major determinant for the cytotoxicity of bleomycin. The DNA sequence specificity of bleomycin cleavage is found to at 5′-GT* and 5′-GC* dinucleotides (where * indicates the cleaved nucleotide). Using next-generation DNA sequencing, over 200 million double-strand breaks were analysed, and an expanded bleomycin sequence specificity was found to be 5′-RTGT*AY (where R is G or A and Y is T or C) in cellular DNA and 5′-TGT*AT in purified DNA. The different environment of cellular DNA compared to purified DNA was proposed to be responsible for the difference. A number of bleomycin analogues have been examined and their interaction with DNA is also discussed. In particular, the production of bleomycin analogues via genetic manipulation of the modular non-ribosomal peptide synthetases and polyketide synthases in the bleomycin gene cluster is reviewed. The prospects for the synthesis of bleomycin analogues with increased effectiveness as cancer chemotherapeutic agents is also explored.

Sequencing of semiflexible polymers of varying bending rigidity using patterned pores

We study the translocation of a semiflexible polymer through extended pores with patterned stickiness, using Langevin dynamics simulations. We find that the consequence of pore patterning on the translocation time dynamics is dramatic and depends strongly on the interplay of polymer stiffness and pore-polymer interactions. For heterogeneous polymers with periodically varying stiffness along their lengths, we find that variation of the block size of the sequences and the orientation results in large variations in the translocation time distributions. We show how this fact may be utilized to develop an effective sequencing strategy. This strategy involving multiple pores with patterned surface energetics can predict heteropolymer sequences having different bending rigidity to a high degree of accuracy.

The sequence specificity of UV-induced DNA damage in a systematically altered DNA sequence

The sequence specificity of UV-induced DNA damage was investigated in a specifically designed DNA plasmid using two procedures: end-labelling and linear amplification. Absorption of UV photons by DNA leads to dimerisation of pyrimidine bases and produces two major photoproducts, cyclobutane pyrimidine dimers (CPDs) and pyrimidine(6-4)pyrimidone photoproducts (6-4PPs). A previous study had determined that two hexanucleotide sequences, 5'-GCTC*AC and 5'-TATT*AA, were high intensity UV-induced DNA damage sites. The UV clone plasmid was constructed by systematically altering each nucleotide of these two hexanucleotide sequences. One of the main goals of this study was to determine the influence of single nucleotide alterations on the intensity of UV-induced DNA damage. The sequence 5'-GCTC*AC was designed to examine the sequence specificity of 6-4PPs and the highest intensity 6-4PP damage sites were found at 5'-GTTC*CC nucleotides. The sequence 5'-TATT*AA was devised to investigate the sequence specificity of CPDs and the highest intensity CPD damage sites were found at 5'-TTTT*CG nucleotides. It was proposed that the tetranucleotide DNA sequence, 5'-YTC*Y (where Y is T or C), was the consensus sequence for the highest intensity UV-induced 6-4PP adduct sites; while it was 5'-YTT*C for the highest intensity UV-induced CPD damage sites. These consensus tetranucleotides are composed entirely of consecutive pyrimidines and must have a DNA conformation that is highly productive for the absorption of UV photons.

Multivalent polyrotaxane vectors as adaptive cargo complexes for gene therapy

The philosophy to design and construct polyrotaxane carriers, as efficient gene delivery systems.

The Sequence Preference of Gamma-Radiation-Induced Damage in End-Labeled DNA after Heat Treatment

In this work, we examined the DNA sequence preference of gamma-radiation-induced DNA damage in purified DNA sequences after heat treatment. DNA was fluorescently end-labeled and gamma-radiation-induced DNA cleavage was examined using capillary electrophoresis with laser-induced fluorescence detection. Our findings provide evidence that gamma-radiation-induced DNA damage to end-labeled DNA is nonrandom and has a sequence preference. The degree of cleavage was quantified at each nucleotide, and we observed that preferential cleavage occurred at C nucleotides with lesser cleavage at G nucleotides, while being very low at T nucleotides. The differences in percentage cleavage at individual nucleotides ranged up to sixfold. The DNA sequences surrounding the most intense radiation-induced DNA cleavage sites were examined and a consensus sequence 5'-AGGC*C (where C* is the cleavage site) was found. The highest intensity gamma-radiation-induced DNA cleavage sites were found at the dinucleotides, 5'-GG*, 5'-GC*, 5'-C*C and 5'-G*G and at the trinucleotides, 5'-GG*C, 5'-TC*A, 5'-GG*G and 5'-GC*C. These findings have implications for our understanding of ionizing radiation-induced DNA damage.

Sequence-Dependent Persistence Length of Long DNA

Using a high-throughput genome-mapping approach, we obtained circa 50 million measurements of the extension of internal human DNA segments in a 41  nm×41  nm nanochannel. The underlying DNA sequences, obtained by mapping to the reference human genome, are 2.5-393 kilobase pairs long and contain percent GC contents between 32.5% and 60%. Using Odijk's theory for a channel-confined wormlike chain, these data reveal that the DNA persistence length increases by almost 20% as the percent GC content increases. The increased persistence length is rationalized by a model, containing no adjustable parameters, that treats the DNA as a statistical terpolymer with a sequence-dependent intrinsic persistence length and a sequence-independent electrostatic persistence length.

The nuclear higher-order structure defined by the set of topological relationships between DNA and the nuclear matrix is species-specific in hepatocytes

During the interphase the nuclear DNA of metazoan cells is organized in supercoiled loops anchored to constituents of a nuclear substructure or compartment known as the nuclear matrix. The stable interactions between DNA and the nuclear matrix (NM) correspond to a set of topological relationships that define a nuclear higher-order structure (NHOS). Current evidence suggests that the NHOS is cell-type-specific. Biophysical evidence and theoretical models suggest that thermodynamic and structural constraints drive the actualization of DNA-NM interactions. However, if the topological relationships between DNA and the NM were the subject of any biological constraint with functional significance then they must be adaptive and thus be positively selected by natural selection and they should be reasonably conserved, at least within closely related species. We carried out a coarse-grained, comparative evaluation of the DNA-NM topological relationships in primary hepatocytes from two closely related mammals: rat and mouse, by determining the relative position to the NM of a limited set of target sequences corresponding to highly-conserved genomic regions that also represent a sample of distinct chromosome territories within the interphase nucleus. Our results indicate that the pattern of topological relationships between DNA and the NM is not conserved between the hepatocytes of the two closely related species, suggesting that the NHOS, like the karyotype, is species-specific.

Structural features of DNA that determine RNA polymerase II core promoter

Background

The general structure and action of all eukaryotic and archaeal RNA polymerases machinery have an astonishing similarity despite the diversity of core promoter sequences in different species. The goal of our work is to find common characteristics of DNA region that define it as a promoter for the RNA polymerase II (Pol II).

Results

The profiles of a large number of physical and structural characteristics, averaged over representative sets of the Pol II minimal core promoters of the evolutionary divergent species from animals, plants and unicellular fungi were analysed. In addition to the characteristics defined at the base-pair steps, we, for the first time, use profiles of the ultrasonic cleavage and DNase I cleavage indexes, informative for internal properties of each complementary strand.

Conclusions

DNA of the core promoters of metazoans and Schizosaccharomyces pombe has similar structural organization. Its mechanical and 3D structural characteristics have singular properties at the positions of TATA-box. The minor groove is broadened and conformational motion is decreased in that region. Special characteristics of conformational behavior are revealed in metazoans at the region, which connects the end of TATA-box and the transcription start site (TSS). The intensities of conformational motions in the complementary strands are periodically changed in opposite phases. They are noticeable, best of all, in mammals. Such conformational features are lacking in the core promoters of S. pombe. The profiles of Saccharomyces cerevisiae core promoters significantly differ: their singular region is shifted down thus pointing to the uniqueness of their structural organization. Obtained results may be useful in genetic engineering for artificial modulation of the promoter strength.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3292-z) contains supplementary material, which is available to authorized users.

Why Cortical Neurons Cannot Divide, and Why Do They Usually Die in the Attempt?

Cortical neurons are prime examples of terminally differentiated, postmitotic cells. However, under experimental or pathological conditions, they can re-enter the cell cycle and replicate DNA but are unable to divide, dying by apoptosis or becoming either polyploid or aneuploid. Any cellular state that depends on the action of genes and their products can be reverted or bypassed by spontaneous or induced mutations, yet there are currently no reports of dividing cortical neurons. Thus, it seems unlikely that the remarkably stable postmitotic condition of cortical neurons depends on specific gene functions. This Review summarizes evidence that the postmitotic state of cortical neurons depends on the high stability of its underlying nuclear structure that results from an entropy-driven process aimed at dissipating the intrinsic structural stress present in chromosomal DNA in such a way that the structural stability of the neuronal nucleus becomes an insurmountable energy barrier for karyokinesis and mitosis. From this perspective, the integral properties of the nuclear higher order structure in neurons provide an explanation not only for why cortical neurons cannot divide but also for why they usually die if they happen to replicate their DNA. © 2016 Wiley Periodicals, Inc.

Reducible Dimeric Conjugates of Small Internally Segment Interfering RNA for Efficient Gene Silencing

The condensation of nucleic acids into compact nanoparticles with cationic carriers is a powerful tool for translocating exogenous nucleic acids into cells. To date, most efforts have been focused on the development of novel gene carriers for safe and efficient gene delivery. However, small interfering RNA (siRNA) is generally not strongly associated with cationic carriers due to its stiff structure and low spatial charge density. To overcome this limitation, this work introduces a well-defined dimeric conjugate of small internally segment interfering RNA (sisiRNA) linked via a disulfide bond for enhanced cellular uptake and gene silencing. Dimeric sisiRNA is synthesized through oxidizing two monomeric sisiRNA molecules, each of which consists of a sense strand carrying a nick and an antisense strand modified with a thiol group at the 3'-end. The nick in the sense strand enables the dimeric sisiRNA to be more effectively condensed into nanosized complexes due to the increased structural flexibility, which results in a higher gene silencing efficiency compared with the dimeric siRNA containing the intact sense strands. The results indicate that the discontinuity of the sense strands is a simple method of adding more flexibility to various siRNA-based nanostructures for enhanced gene silencing.

A historical account of Hoogsteen base-pairs in duplex DNA

In 1957, a unique pattern of hydrogen bonding between N3 and O4 on uracil and N7 and N6 on adenine was proposed to explain how poly(rU) strands can associate with poly(rA)-poly(rU) duplexes to form triplexes. Two years later, Karst Hoogsteen visualized such a noncanonical A-T base-pair through X-ray analysis of co-crystals containing 9-methyladenine and 1-methylthymine. Subsequent X-ray analyses of guanine and cytosine derivatives yielded the expected Watson-Crick base-pairing, but those of adenine and thymine (or uridine) did not yield Watson-Crick base-pairs, instead favoring "Hoogsteen" base-pairing. More than two decades ensued without experimental "proof" for A-T Watson-Crick base-pairs, while Hoogsteen base-pairs continued to surface in AT-rich sequences, closing base-pairs of apical loops, in structures of DNA bound to antibiotics and proteins, damaged and chemically modified DNA, and in polymerases that replicate DNA via Hoogsteen pairing. Recently, NMR studies have shown that base-pairs in duplex DNA exist as a dynamic equilibrium between Watson-Crick and Hoogsteen forms. There is now little doubt that Hoogsteen base-pairs exist in significant abundance in genomic DNA, where they can expand the structural and functional versatility of duplex DNA beyond that which can be achieved based only on Watson-Crick base-pairing. Here, we provide a historical account of the discovery and characterization of Hoogsteen base-pairs, hoping that this will inform future studies exploring the occurrence and functional importance of these alternative base-pairs.

Organization of DNA in a bacterial nucleoid

BACKGROUND

It is unclear how DNA is packaged in a bacterial cell in the absence of nucleosomes. To investigate the initial level of DNA condensation in bacterial nucleoid we used in vivo DNA digestion coupled with high-throughput sequencing of the digestion-resistant fragments. To this end, we transformed E. coli cells with a plasmid expressing micrococcal nuclease. The nuclease expression was under the control of AraC repressor, which enabled us to perform an inducible digestion of bacterial nucleoid inside a living cell.

RESULTS

Analysis of the genomic localization of the digestion-resistant fragments revealed their non-random distribution. The patterns observed in the distribution of the sequenced fragments indicate the presence of short DNA segments protected from the enzyme digestion, possibly because of interaction with DNA-binding proteins. The average length of such digestion-resistant segments is about 50 bp and the characteristic repeat in their distribution is about 90 bp. The gene starts are depleted of the digestion-resistant fragments, suggesting that these genomic regions are more exposed than genomic sequences on average. Sequence analysis of the digestion-resistant segments showed that while the GC-content of such sequences is close to the genome-wide value, they are depleted of A-tracts as compared to the bulk genomic DNA or to the randomized sequence of the same nucleotide composition.

CONCLUSIONS

Our results suggest that DNA is packaged in the bacterial nucleoid in a non-random way that facilitates interaction of the DNA binding factors with regulatory regions of the genome.

Elastic Energy Partitioning in DNA Deformation and Binding to Proteins

We study the elasticity of DNA based on local principal axes of bending identified from over 0.9-μs all-atom molecular dynamics simulations of DNA oligos. The calculated order parameters describe motion of DNA as an elastic rod. In 10 possible dinucleotide steps, bending about the two principal axes is anisotropic yet linearly elastic. Twist about the centroid axis is largely decoupled from bending, but DNA tends to overtwist for unbending beyond the typical range of thermal motion, which is consistent with experimentally observed twist-stretch coupling. The calculated elastic stiffness of dinucleotide steps yield sequence-dependent persistence lengths consistent with previous single-molecule experiments, which is further analyzed by performing coarse-grained simulations of DNA. Flexibility maps of oligos constructed from simulation also match with those from the precalculated stiffness of dinucleotide steps. These support the premise that base pair interaction at the dinucleotide-level is mainly responsible for the elasticity of DNA. Furthermore, we analyze 1381 crystal structures of protein-DNA complexes. In most structures, DNAs are mildly deformed and twist takes the highest portion of the total elastic energy. By contrast, in structures with the elastic energy per dinucleotide step greater than about 4.16 kBT (kBT: thermal energy), the major bending becomes dominant. The extensional energy of dinucleotide steps takes at most 35% of the total elastic energy except for structures containing highly deformed DNAs where linear elasticity breaks down. Such partitioning between different deformational modes provides quantitative insights into the conformational dynamics of DNA as well as its interaction with other molecules and surfaces.

The interphase mammalian chromosome as a structural system based on tensegrity

Each mammalian chromosome is constituted by a DNA fiber of macroscopic length that needs to be fitted in a microscopic nucleus. The DNA fiber is subjected at physiological temperature to random thermal bending and looping that must be constrained so as achieve structural stability thus avoiding spontaneous rupturing of the fiber. Standard textbooks assume that chromatin proteins are primarily responsible for the packaging of DNA and so of its protection against spontaneous breakage. Yet the dynamic nature of the interactions between chromatin proteins and DNA is unlikely to provide the necessary long-term structural stability for the chromosomal DNA. On the other hand, longstanding evidence indicates that stable interactions between DNA and constituents of a nuclear compartment commonly known as the nuclear matrix organize the chromosomal DNA as a series of topologically constrained, supercoiled loops during interphase. This results in a primary level of DNA condensation and packaging within the nucleus, as well as in protection against spontaneous DNA breakage, independently of chromatin proteins which nevertheless increase and dynamically modulate the degree of DNA packaging and its role in the regulation of DNA function. Thus current evidence, presented hereunder, supports a model for the organization of the interphase chromosome as resilient system that satisfies the principles of structural tensegrity.

Modulations of DNA Contacts by Linker Histones and Post-translational Modifications Determine the Mobility and Modifiability of Nucleosomal H3 Tails

Post-translational histone modifications and linker histone incorporation regulate chromatin structure and genome activity. How these systems interface on a molecular level is unclear. Using biochemistry and NMR spectroscopy, we deduced mechanistic insights into the modification behavior of N-terminal histone H3 tails in different nucleosomal contexts. We find that linker histones generally inhibit modifications of different H3 sites and reduce H3 tail dynamics in nucleosomes. These effects are caused by modulations of electrostatic interactions of H3 tails with linker DNA and largely depend on the C-terminal domains of linker histones. In agreement, linker histone occupancy and H3 tail modifications segregate on a genome-wide level. Charge-modulating modifications such as phosphorylation and acetylation weaken transient H3 tail-linker DNA interactions, increase H3 tail dynamics, and, concomitantly, enhance general modifiability. We propose that alterations of H3 tail-linker DNA interactions by linker histones and charge-modulating modifications execute basal control mechanisms of chromatin function.

Dynamic DNA Nanotubes: Reversible Switching between Single and Double-Stranded Forms, and Effect of Base Deletions

DNA nanotubes hold great potential as drug delivery vehicles and as programmable templates for the organization of materials and biomolecules. Existing methods for their construction produce assemblies that are entirely double-stranded and rigid, and thus have limited intrinsic dynamic character, or they rely on chemically modified and ligated DNA structures. Here, we report a simple and efficient synthesis of DNA nanotubes from 11 short unmodified strands, and the study of their dynamic behavior by atomic force microscopy and in situ single molecule fluorescence microscopy. This method allows the programmable introduction of DNA structural changes within the repeat units of the tubes. We generate and study fully double-stranded nanotubes, and convert them to nanotubes with one, two and three single-stranded sides, using strand displacement strategies. The nanotubes can be reversibly switched between these forms without compromising their stability and micron-scale lengths. We then site-specifically introduce DNA strands that shorten two sides of the nanotubes, while keeping the length of the third side. The nanotubes undergo bending with increased length mismatch between their sides, until the distortion is significant enough to shorten them, as measured by AFM and single-molecule fluorescence photobleaching experiments. The method presented here produces dynamic and robust nanotubes that can potentially behave as actuators, and allows their site-specific addressability while using a minimal number of component strands.

Understanding the structural and dynamic consequences of DNA epigenetic modifications: computational insights into cytosine methylation and hydroxymethylation

We report a series of molecular dynamics (MD) simulations of up to a microsecond combined simulation time designed to probe epigenetically modified DNA sequences. More specifically, by monitoring the effects of methylation and hydroxymethylation of cytosine in different DNA sequences, we show, for the first time, that DNA epigenetic modifications change the molecule's dynamical landscape, increasing the propensity of DNA toward different values of twist and/or roll/tilt angles (in relation to the unmodified DNA) at the modification sites. Moreover, both the extent and position of different modifications have significant effects on the amount of structural variation observed. We propose that these conformational differences, which are dependent on the sequence environment, can provide specificity for protein binding.

Association of new deletion/duplication region at chromosome 1p21 with intellectual disability, severe speech deficit and autism spectrum disorder-like behavior: an all-in approach to solving the DPYD enigma

We describe an as yet unreported neocentric small supernumerary marker chromosome (sSMC) derived from chromosome 1p21.3p21.2. It was present in 80% of the lymphocytes in a male patient with intellectual disability, severe speech deficit, mild dysmorphic features, and hyperactivity with elements of autism spectrum disorder (ASD).

Several important neurodevelopmental genes are affected by the 3.56 Mb copy number gain of 1p21.3p21.2, which may be considered reciprocal in gene content to the recently recognized 1p21.3 microdeletion syndrome. Both 1p21.3 deletions and the presented duplication display overlapping symptoms, fitting the same disorder category. Contribution of coding and non-coding genes to the phenotype is discussed in the light of cellular and intercellular homeostasis disequilibrium. In line with this the presented 1p21.3p21.2 copy number gain correlated to 1p21.3 microdeletion syndrome verifies the hypothesis of a cumulative effect of the number of deregulated genes - homeostasis disequilibrium leading to overlapping phenotypes between microdeletion and microduplication syndromes.

Although miR-137 appears to be the major player in the 1p21.3p21.2 region, deregulation of the DPYD (dihydropyrimidine dehydrogenase) gene may potentially affect neighboring genes underlying the overlapping symptoms present in both the copy number loss and copy number gain of 1p21. Namely, the all-in approach revealed that DPYD is a complex gene whose expression is epigenetically regulated by long non-coding RNAs (lncRNAs) within the locus. Furthermore, the long interspersed nuclear element-1 (LINE-1) L1MC1 transposon inserted in DPYD intronic transcript 1 (DPYD-IT1) lncRNA with its parasites, TcMAR-Tigger5b and pair of Alu repeats appears to be the “weakest link” within the DPYD gene liable to break. Identification of the precise mechanism through which DPYD is epigenetically regulated, and underlying reasons why exactly the break (FRA1E) happens, will consequently pave the way toward preventing severe toxicity to the antineoplastic drug 5-fluorouracil (5-FU) and development of the causative therapy for the dihydropyrimidine dehydrogenase deficiency.

Atomic force microscopy study of DNA flexibility on short length scales: smooth bending versus kinking

The apparently anomalous flexibility of DNA on short length scales has attracted a lot of attention in recent years. We use atomic force microscopy (AFM) in solution to directly study the DNA bending statistics for small lengths down to one helical turn. The accuracy of experimental estimates could be improved due to a large data volume and a refined algorithm for image processing and measuring bend angles. It is found that, at length scales beyond two helical turns (7 nm), DNA is well described by the harmonic worm-like chain (WLC) model with the bending persistence length of 56 nm. Below this threshold, the AFM data are also described by the WLC model assuming that the accuracy of measured bend angles is limited by the physical width of the double helix. We conclude that the double helical DNA behaves as a uniform elastic rod even at very short length scales. Strong bends due to kinks, melting bubbles and other deviations from the WLC model are statistically negligible.

A novel indirect sequence readout component in the E. coli cyclic AMP receptor protein operator

The cyclic AMP receptor protein (CRP) from Escherichia coli has been extensively studied for several decades. In particular, a detailed characterization of CRP interaction with DNA has been obtained. The CRP dimer recognizes a consensus sequence AANTGTGANNNNNNTCACANTT through direct amino acid nucleobase interactions in the major groove of the two operator half-sites. Crystal structure analyses have revealed that the interaction results in two strong kinks at the TG/CA steps closest to the 6-base-pair spacer (N6). This spacer exhibits high sequence variability among the more than 100 natural binding sites in the E. coli genome, but the exact role of the N6 region in CRP interaction has not previously been systematic examined. Here we employ an in vitro selection system based on a randomized N6 spacer region to demonstrate that CRP binding to the lacP1 site may be enhanced up to 14-fold or abolished by varying the N6 spacer sequences. Furthermore, on the basis of sequence analysis and uranyl (UO2(2+)) probing data, we propose that the underlying mechanism relies on N6 deformability.

Non-specificity and synergy at the binding site of the carboplatin-induced DNA adduct via molecular dynamics simulations of the MutSα-DNA recognition complex

MutSα is the most abundant mismatch-binding factor of human DNA mismatch repair (MMR) proteins. MMR maintains genetic stability by recognizing and repairing DNA defects. Failure to accomplish their function may lead to cancer. In addition, MutSα recognizes at least some types of DNA damage making it a target for anticancer agents. Here, complementing scarce experimental data, we report unique hydrogen-bonding motifs associated with the recognition of the carboplatin induced DNA damage by MutSα. These data predict that carboplatin and cisplatin induced damaging DNA adducts are recognized by MutSα in a similar manner. Our simulations also indicate that loss of base pairing at the damage site results in (1) non-specific binding and (2) changes in the atomic flexibility at the lesion site and beyond. To further quantify alterations at MutSα-DNA interface in response to damage recognition, non-bonding interactions and salt bridges were investigated. These data indicate (1) possible different packing and (2) disruption of the salt bridges at the MutSα-DNA interface in the damaged complex. These findings (1) underscore the general observation of disruptions at the MutSα-DNA interface and (2) highlight the nature of the anticancer effect of the carboplatin agent. The analysis was carried out from atomistic simulations.