Crystal structure of a human TATA box-binding protein/TATA element complex

ATOMDANCE: Kernel-based denoising and choreographic analysis for protein dynamic comparison

Comparative methods in molecular evolution and structural biology rely heavily upon the site-wise analysis of DNA sequence and protein structure, both static forms of information. However, it is widely accepted that protein function results from nanoscale nonrandom machine-like motions induced by evolutionarily conserved molecular interactions. Comparisons of molecular dynamics (MD) simulations conducted between homologous sites representative of different functional or mutational states can potentially identify local effects on binding interaction and protein evolution. In addition, comparisons of different (i.e., nonhomologous) sites within MD simulations could be employed to identify functional shifts in local time-coordinated dynamics indicative of logic gating within proteins. However, comparative MD analysis is challenged by the large fraction of protein motion caused by random thermal noise in the surrounding solvent. Therefore, properly denoised MD comparisons could reveal functional sites involving these machine-like dynamics with good accuracy. Here, we introduce ATOMDANCE, a user-interfaced suite of comparative machine learning-based denoising tools designed for identifying functional sites and the patterns of coordinated motion they can create within MD simulations. ATOMDANCE-maxDemon4.0 employs Gaussian kernel functions to compute site-wise maximum mean discrepancy between learned features of motion, thereby assessing denoised differences in the nonrandom motions between functional or evolutionary states (e.g., ligand bound versus unbound, wild-type versus mutant). ATOMDANCE-maxDemon4.0 also employs maximum mean discrepancy to analyze potential random amino acid replacements allowing for a site-wise test of neutral versus nonneutral evolution on the divergence of dynamic function in protein homologs. Finally, ATOMDANCE-Choreograph2.0 employs mixed-model analysis of variance and graph network to detect regions where time-synchronized shifts in dynamics occur. Here, we demonstrate ATOMDANCE's utility for identifying key sites involved in dynamic responses during functional binding interactions involving DNA, small-molecule drugs, and virus-host recognition, as well as understanding shifts in global and local site coordination occurring during allosteric activation of a pathogenic protease.

Re-pairing DNA: binding of a ruthenium phi complex to a double mismatch

We report a crystal structure at atomic resolution (0.9 Å) of a ruthenium complex bound to a consecutive DNA double mismatch, which results in a TA basepair with flipped out thymine, together with the formation of an adenine bulge. The structure shows a form of metalloinsertion interaction of the Λ-[Ru(phen)2phi]2+ (phi = 9,10-phenanthrenediimine) complex at the bulge site. The metal complex interacts with the DNA via the major groove, where specific interactions between the adenines of the DNA and the phen ligands of the complex are formed. One Δ-[Ru(phen)2phi]2+ complex interacts via the minor groove, which shows sandwiching of its phi ligand between the phi ligands of the other two ruthenium complexes, and no interaction of its phen ligands with DNA. To our knowledge, this binding model represents a new form of metalloinsertion in showing major rather than minor groove insertion.

Cryo-EM Structure of AAV2 Rep68 bound to integration site AAVS1: Insights into the mechanism of DNA melting

The Rep68 protein from Adeno-Associated Virus (AAV) is a multifunctional SF3 helicase that performs most of the DNA transactions required for the viral life cycle. During AAV DNA replication, Rep68 assembles at the origin and catalyzes the DNA melting and nicking reactions during the hairpin rolling replication process to complete the second-strand synthesis of the AAV genome. Here, we report the Cryo-EM structures of Rep68 bound to double-stranded DNA (dsDNA) containing the sequence of the AAVS1 integration site in different nucleotide-bound states. In the apo state, Rep68 forms a heptameric complex around DNA, with three Origin Binding Domains (OBDs) bound to the Rep Binding Site (RBS) sequence and three other OBDs forming transient dimers with them. The AAA+ domains form an open ring with no interactions between subunits and with DNA. We hypothesize the heptameric quaternary structure is necessary to load onto dsDNA. In the ATPγS-bound state, a subset of three subunits binds the nucleotide, undergoing a large conformational change, inducing the formation of intersubunit interactions interaction and interaction with three consecutive DNA phosphate groups. Moreover, the induced conformational change positions three phenylalanine residues to come in close contact with the DNA backbone, producing a distortion in the DNA. We propose that the phenylalanine residues can potentially act as a hydrophobic wedge in the DNA melting process.

Brewing COFFEE: A Sequence-Specific Coarse-Grained Energy Function for Simulations of DNA-Protein Complexes

DNA-protein interactions are pervasive in a number of biophysical processes ranging from transcription and gene expression to chromosome folding. To describe the structural and dynamic properties underlying these processes accurately, it is important to create transferable computational models. Toward this end, we introduce Coarse-grained Force Field for Energy Estimation, COFFEE, a robust framework for simulating DNA-protein complexes. To brew COFFEE, we integrated the energy function in the self-organized polymer model with side-chains for proteins and the three interaction site model for DNA in a modular fashion, without recalibrating any of the parameters in the original force-fields. A unique feature of COFFEE is that it describes sequence-specific DNA-protein interactions using a statistical potential (SP) derived from a data set of high-resolution crystal structures. The only parameter in COFFEE is the strength (λDNAPRO) of the DNA-protein contact potential. For an optimal choice of λDNAPRO, the crystallographic B-factors for DNA-protein complexes with varying sizes and topologies are quantitatively reproduced. Without any further readjustments to the force-field parameters, COFFEE predicts scattering profiles that are in quantitative agreement with small-angle X-ray scattering experiments, as well as chemical shifts that are consistent with NMR. We also show that COFFEE accurately describes the salt-induced unraveling of nucleosomes. Strikingly, our nucleosome simulations explain the destabilization effect of ARG to LYS mutations, which do not alter the balance of electrostatic interactions but affect chemical interactions in subtle ways. The range of applications attests to the transferability of COFFEE, and we anticipate that it would be a promising framework for simulating DNA-protein complexes at the molecular length-scale.

Mechanisms and Functions of the RNA Polymerase II General Transcription Machinery during the Transcription Cycle

Central to the development and survival of all organisms is the regulation of gene expression, which begins with the process of transcription catalyzed by RNA polymerases. During transcription of protein-coding genes, the general transcription factors (GTFs) work alongside RNA polymerase II (Pol II) to assemble the preinitiation complex at the transcription start site, open the promoter DNA, initiate synthesis of the nascent messenger RNA, transition to productive elongation, and ultimately terminate transcription. Through these different stages of transcription, Pol II is dynamically phosphorylated at the C-terminal tail of its largest subunit, serving as a control mechanism for Pol II elongation and a signaling/binding platform for co-transcriptional factors. The large number of core protein factors participating in the fundamental steps of transcription add dense layers of regulation that contribute to the complexity of temporal and spatial control of gene expression within any given cell type. The Pol II transcription system is highly conserved across different levels of eukaryotes; however, most of the information here will focus on the human Pol II system. This review walks through various stages of transcription, from preinitiation complex assembly to termination, highlighting the functions and mechanisms of the core machinery that participates in each stage.

High UV damage and low repair, but not cytosine deamination, stimulate mutation hotspots at ETS binding sites in melanoma

Noncoding mutation hotspots have been identified in melanoma and many of them occur at the binding sites of E26 transformation-specific (ETS) proteins; however, their formation mechanism and functional impacts are not fully understood. Here, we used UV (Ultraviolet) damage sequencing data and analyzed cyclobutane pyrimidine dimer (CPD) formation, DNA repair, and CPD deamination in human cells at single-nucleotide resolution. Our data show prominent CPD hotspots immediately after UV irradiation at ETS binding sites, particularly at sites with a conserved TTCCGG motif, which correlate with mutation hotspots identified in cutaneous melanoma. Additionally, CPDs are repaired slower at ETS binding sites than in flanking DNA. Cytosine deamination in CPDs to uracil is suggested as an important step for UV mutagenesis. However, we found that CPD deamination is significantly suppressed at ETS binding sites, particularly for the CPD hotspot on the 5' side of the ETS motif, arguing against a role for CPD deamination in promoting ETS-associated UV mutations. Finally, we analyzed a subset of frequently mutated promoters, including the ribosomal protein genes RPL13A and RPS20, and found that mutations in the ETS motif can significantly reduce the promoter activity. Thus, our data identify high UV damage and low repair, but not CPD deamination, as the main mechanism for ETS-associated mutations in melanoma and uncover important roles of often-overlooked mutation hotspots in perturbing gene transcription.

Molecular details of the CPSF73-CPSF100 C-terminal heterodimer and interaction with Symplekin

Eukaryotic pre-mRNA is processed by a large multiprotein complex to accurately cleave the 3' end, and to catalyse the addition of the poly(A) tail. Within this cleavage and polyadenylation specificity factor (CPSF) machinery, the CPSF73/CPSF3 endonuclease subunit directly contacts both CPSF100/CPSF2 and the scaffold protein Symplekin to form a subcomplex known as the core cleavage complex or mammalian cleavage factor. Here we have taken advantage of a stable CPSF73-CPSF100 minimal heterodimer from Encephalitozoon cuniculi to determine the solution structure formed by the first and second C-terminal domain (CTD1 and CTD2) of both proteins. We find a large number of contacts between both proteins in the complex, and notably in the region between CTD1 and CTD2. A similarity is also observed between CTD2 and the TATA-box binding protein (TBP) domains. Separately, we have determined the structure of the terminal CTD3 domain of CPSF73, which also belongs to the TBP domain family and is connected by a flexible linker to the rest of CPSF73. Biochemical assays demonstrate a key role for the CTD3 of CPSF73 in binding Symplekin, and structural models of the trimeric complex from other species allow for comparative analysis and support an overall conserved architecture.

Brewing COFFEE: A sequence-specific coarse-grained energy function for simulations of DNA-protein complexes

DNA-protein interactions are pervasive in a number of biophysical processes ranging from transcription, gene expression, to chromosome folding. To describe the structural and dynamic properties underlying these processes accurately, it is important to create transferable computational models. Toward this end, we introduce Coarse grained force field for energy estimation, COFFEE, a robust framework for simulating DNA-protein complexes. To brew COFFEE, we integrated the energy function in the Self-Organized Polymer model with Side Chains for proteins and the Three Interaction Site model for DNA in a modular fashion, without re-calibrating any of the parameters in the original force-fields. A unique feature of COFFEE is that it describes sequence-specific DNA-protein interactions using a statistical potential (SP) derived from a dataset of high-resolution crystal structures. The only parameter in COFFEE is the strength (λ D N A P R O ) of the DNA-protein contact potential. For an optimal choice of λ D N A P R O , the crystallographic B-factors for DNA-protein complexes, with varying sizes and topologies, are quantitatively reproduced. Without any further readjustments to the force-field parameters, COFFEE predicts the scattering profiles that are in quantitative agreement with SAXS experiments as well as chemical shifts that are consistent with NMR. We also show that COFFEE accurately describes the salt-induced unraveling of nucleosomes. Strikingly, our nucleosome simulations explain the destabilization effect of ARG to LYS mutations, which does not alter the balance of electrostatic interactions, but affects chemical interactions in subtle ways. The range of applications attests to the transferability of COFFEE, and we anticipate that it would be a promising framework for simulating DNA-protein complexes at the molecular length-scale.

Structural predictions of protein-DNA binding: MELD-DNA

Structural, regulatory and enzymatic proteins interact with DNA to maintain a healthy and functional genome. Yet, our structural understanding of how proteins interact with DNA is limited. We present MELD-DNA, a novel computational approach to predict the structures of protein-DNA complexes. The method combines molecular dynamics simulations with general knowledge or experimental information through Bayesian inference. The physical model is sensitive to sequence-dependent properties and conformational changes required for binding, while information accelerates sampling of bound conformations. MELD-DNA can: (i) sample multiple binding modes; (ii) identify the preferred binding mode from the ensembles; and (iii) provide qualitative binding preferences between DNA sequences. We first assess performance on a dataset of 15 protein-DNA complexes and compare it with state-of-the-art methodologies. Furthermore, for three selected complexes, we show sequence dependence effects of binding in MELD predictions. We expect that the results presented herein, together with the freely available software, will impact structural biology (by complementing DNA structural databases) and molecular recognition (by bringing new insights into aspects governing protein-DNA interactions).

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.

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.

Polyglutamine Ataxias: Our Current Molecular Understanding and What the Future Holds for Antisense Therapies

Polyglutamine (polyQ) ataxias are a heterogenous group of neurological disorders all caused by an expanded CAG trinucleotide repeat located in the coding region of each unique causative gene. To date, polyQ ataxias encompass six disorders: spinocerebellar ataxia types 1, 2, 3, 6, 7, and 17 and account for a larger group of disorders simply known as polyglutamine disorders, which also includes Huntington's disease. These diseases are typically characterised by progressive ataxia, speech and swallowing difficulties, lack of coordination and gait, and are unfortunately fatal in nature, with the exception of SCA6. All the polyQ spinocerebellar ataxias have a hallmark feature of neuronal aggregations and share many common pathogenic mechanisms, such as mitochondrial dysfunction, impaired proteasomal function, and autophagy impairment. Currently, therapeutic options are limited, with no available treatments that slow or halt disease progression. Here, we discuss the common molecular and clinical presentations of polyQ spinocerebellar ataxias. We will also discuss the promising antisense oligonucleotide therapeutics being developed as treatments for these devastating diseases. With recent advancements and therapeutic approvals of various antisense therapies, it is envisioned that some of the studies reviewed may progress into clinical trials and beyond.

The Long Road to Understanding RNAPII Transcription Initiation and Related Syndromes

In eukaryotes, transcription of protein-coding genes requires the assembly at core promoters of a large preinitiation machinery containing RNA polymerase II (RNAPII) and general transcription factors (GTFs). Transcription is potentiated by regulatory elements called enhancers, which are recognized by specific DNA-binding transcription factors that recruit cofactors and convey, following chromatin remodeling, the activating cues to the preinitiation complex. This review summarizes nearly five decades of work on transcription initiation by describing the sequential recruitment of diverse molecular players including the GTFs, the Mediator complex, and DNA repair factors that support RNAPII to enable RNA synthesis. The elucidation of the transcription initiation mechanism has greatly benefited from the study of altered transcription components associated with human diseases that could be considered transcription syndromes.

A Primer on the oxDNA Model of DNA: When to Use it, How to Simulate it and How to Interpret the Results

The oxDNA model of Deoxyribonucleic acid has been applied widely to systems in biology, biophysics and nanotechnology. It is currently available via two independent open source packages. Here we present a set of clearly documented exemplar simulations that simultaneously provide both an introduction to simulating the model, and a review of the model's fundamental properties. We outline how simulation results can be interpreted in terms of-and feed into our understanding of-less detailed models that operate at larger length scales, and provide guidance on whether simulating a system with oxDNA is worthwhile.

Tumor Susceptibility Gene 101 Regulates the Glucocorticoid Receptor through Disorder-Mediated Allostery

Tumor susceptibility gene 101 (TSG101) is involved in endosomal maturation and has been implicated in the transcriptional regulation of several steroid hormone receptors, although a detailed characterization of such regulation has yet to be conducted. Here we directly measure binding of TSG101 to one steroid hormone receptor, the glucocorticoid receptor (GR). Using biophysical and cellular assays, we show that the coiled-coil domain of TSG101 (1) binds and folds the disordered N-terminal domain of the GR, (2) upon binding improves the DNA binding of the GR in vitro, and (3) enhances the transcriptional activity of the GR in vivo. Our findings suggest that TSG101 is a bona fide transcriptional co-regulator of the GR and reveal how the underlying thermodynamics affect the function of the GR.

Sequence-Dependent Kink Formation in Short DNA Loops: Theory and Molecular Dynamics Simulations

Kink formation is essential in highly bent DNA complexed with gene regulatory proteins such as histones to release the bending stress stored within the DNA duplex. Local opening of the double-stranded DNA creates a sharp turn along the specific sequence, which leads to the global bending of the DNA strand. Despite the critical role of kink formation, it is still challenging to predict the position of kink formation for a given DNA sequence. In this study, we propose a theoretical model and perform molecular dynamics simulations to quantify the sequence-dependent kink probability of a strongly bent DNA. By incorporating the elastic bending energy and the sequence-specific thermodynamic parameters, we investigate the importance of the DNA sequence on kink formation. We find that the sequence with TA dinucleotide repeats flanked by GC steps increases the kink propensity by more than an order of magnitude under the same bending stress. The number of base pairs involved in the local opening is found to be coupled with the sequence-specific bubble formation free energy. Our study elucidates the molecular origin of the sequence heterogeneity on kink formation, which is fundamental to understanding protein-DNA recognition.

Kinetics of the Conformational Transformation between B- and A-Forms in the Drew-Dickerson Dodecamer

Some DNA sequences in crystals and in complexes with proteins can exist in the forms intermediate between the B- and A-DNA. Based on this, it was implied that the B-to-A transition for any DNA molecule should go through these intermediate forms also in kinetics. More precisely, the helix parameter Slide has to change first, and the molecule should take the E-form. After that, the Roll parameter changes. In the present work, we simulated the kinetics of the B-A transition in the Drew-Dickerson dodecamer, a known B-philic DNA oligomer. We used the "sugar" coarse-grained model that reproduces ribose flexibility, preserves sequence specificity, employs implicit water and explicit ions, and offers the possibility to vary friction. As the control parameter of the transition, we chose the volume available for a counterion and considered the change from a large to a small volume. In the described system, the B-to-A conformational transformation proved to correspond to a first-order phase transition. The molecule behaves like a small cluster in the region of such a transition, jumping between the A- and B-forms in a wide range of available volumes. The viscosity of the solvent does not affect the midpoint of the transition but only the overall mobility of the system. All helix parameters change synchronously on average, we have not observed the sequence "Slide first, Roll later" in kinetics, and the E-DNA is not a necessary step for the transition between the B- and A-forms in the studied system. So, the existence of the intermediate DNA forms requires specific conditions, shifting the common balance of interactions: certain nucleotide sequence in specific solution or/and the interaction with some protein.

Chromatin-modifying drugs and metabolites in cell fate control

For multicellular organisms, it is essential to produce a variety of specialized cells to perform a dazzling panoply of functions. Chromatin plays a vital role in determining cellular identities, and it dynamically regulates gene expression in response to changing nutrient metabolism and environmental conditions. Intermediates produced by cellular metabolic pathways are used as cofactors or substrates for chromatin modification. Drug analogues of metabolites that regulate chromatin-modifying enzyme reactions can also regulate cell fate by adjusting chromatin organization. In recent years, there have been many studies about how chromatin-modifying drug molecules or metabolites can interact with chromatin to regulate cell fate. In this review, we systematically discuss how DNA and histone-modifying molecules alter cell fate by regulating chromatin conformation and propose a mechanistic model that explains the process of cell fate transitions in a concise and qualitative manner.

The highly diverse TATA box-binding proteins among protists: A review

Transcription is the first step of gene expression regulation and is a fundamental mechanism for establishing the viability and development of a cell. The TATA box-binding protein (TBP) interaction with a TATA box in a promoter is one of the best studied mechanisms in transcription initiation. TBP is a transcription factor that is highly conserved from archaea to humans and is essential for the transcription initiated by each of the three RNA polymerases. In addition, the discovery of TBP-related factor 1 (TRF1) and other factors related to TBP shed light on the variability among transcription initiation complexes, thus demonstrating that the compositions of these complexes are, in fact, more complicated than originally believed. Despite these facts, the majority of studies on transcription have been performed on animal, plant and fungal cells, which serve as canonical models, and information regarding protist cells is relatively scarce. The aim of this work is to review the diversity of the TBPs that have been documented in protists and describe some of the specific features that differentiate them from their counterparts in higher eukaryotes.

Motif grammar: The basis of the language of gene expression

Collaboration of transcription factors (TFs) and their recognition motifs in DNA is the result of coevolution and forms the basis of gene regulation. However, the way how these short genomic sequences contribute to setting the level of gene products is not understood in sufficient detail. The biological problem to be solved by the cell is complex, because each gene requires a unique regulatory network in each cellular condition using the same genome. Thus far, only some components of these networks have been uncovered. In this review, we compiled the features and principles of the motif grammar, which dictates the characteristics and thus the likelihood of the interactions of the binding TFs and their coregulators. We present how sequence features provide specificity using, as examples, two major TF superfamilies, the bZIP proteins and nuclear receptors. We also discuss the phenomenon of “weak” (low affinity) binding sites, which appear to be components of several important genomic regulatory regions, but paradoxically are barely detectable by the currently used approaches. Assembling the complete set of regulatory regions composed of both weak and strong binding sites will allow one to get more comprehensive lists of factors playing roles in gene regulation, thus making possible the deeper understanding of regulatory networks.

A unified dinucleotide alphabet describing both RNA and DNA structures

By analyzing almost 120 000 dinucleotides in over 2000 nonredundant nucleic acid crystal structures, we define 96+1 diNucleotide Conformers, NtCs, which describe the geometry of RNA and DNA dinucleotides. NtC classes are grouped into 15 codes of the structural alphabet CANA (Conformational Alphabet of Nucleic Acids) to simplify symbolic annotation of the prominent structural features of NAs and their intuitive graphical display. The search for nontrivial patterns of NtCs resulted in the identification of several types of RNA loops, some of them observed for the first time. Over 30% of the nearly six million dinucleotides in the PDB cannot be assigned to any NtC class but we demonstrate that up to a half of them can be re-refined with the help of proper refinement targets. A statistical analysis of the preferences of NtCs and CANA codes for the 16 dinucleotide sequences showed that neither the NtC class AA00, which forms the scaffold of RNA structures, nor BB00, the DNA most populated class, are sequence neutral but their distributions are significantly biased. The reported automated assignment of the NtC classes and CANA codes available at dnatco.org provides a powerful tool for unbiased analysis of nucleic acid structures by structural and molecular biologists.

Atomic-resolution mapping of transcription factor-DNA interactions by femtosecond laser crosslinking and mass spectrometry

Transcription factors (TFs) regulate target genes by specific interactions with DNA sequences. Detecting and understanding these interactions at the molecular level is of fundamental importance in biological and clinical contexts. Crosslinking mass spectrometry is a powerful tool to assist the structure prediction of protein complexes but has been limited to the study of protein-protein and protein-RNA interactions. Here, we present a femtosecond laser-induced crosslinking mass spectrometry (fliX-MS) workflow, which allows the mapping of protein-DNA contacts at single nucleotide and up to single amino acid resolution. Applied to recombinant histone octamers, NF1, and TBP in complex with DNA, our method is highly specific for the mapping of DNA binding domains. Identified crosslinks are in close agreement with previous biochemical data on DNA binding and mostly fit known complex structures. Applying fliX-MS to cells identifies several bona fide crosslinks on DNA binding domains, paving the way for future large scale ex vivo experiments.

Revealing the Structural Contributions to Thermal Adaptation of the TATA-Box Binding Protein: Molecular Dynamics and QSPR Analyses

The TATA-box binding protein (TBP) is an important element of the transcription machinery in archaea and eukaryotic organisms. TBP is expressed in organisms adapted to different temperatures, indicating a robust structure, and experimental studies have shown that the mid-unfolding temperature (Tm) of TBP is directly correlated with the optimal growth temperature (OGT) of the organism. To understand which are the relevant structural requirements for its stability, we present the first structural and dynamic computational study of TBPs, combining molecular dynamics (MD) simulations and a quantitative structure-property relationship (QSPR) over a set of TBPs of organisms adapted to different temperatures. We found that the main structural properties of TBP used to adapt to high temperatures are an increase in the ease of desolvation of charged residues at the surface, an increase in the local resiliency, the presence of Leu clusters in the protein core, and an increase in the loss of hydrophobic packing in the N-terminal subdomain. In view of our results, we consider that TBP is a good model to study thermal adaptation, and our analysis opens the possibility of performing protein engineering on TBPs to study transcription at high or low temperatures.

Host Transcription Factors in Hepatitis B Virus RNA Synthesis

The hepatitis B virus (HBV) chronically infects over 250 million people worldwide and is one of the leading causes of liver cancer and hepatocellular carcinoma. HBV persistence is due in part to the highly stable HBV minichromosome or HBV covalently closed circular DNA (cccDNA) that resides in the nucleus. As HBV replication requires the help of host transcription factors to replicate, focusing on host protein–HBV genome interactions may reveal insights into new drug targets against cccDNA. The structural details on such complexes, however, remain poorly defined. In this review, the current literature regarding host transcription factors’ interactions with HBV cccDNA is discussed.

DROIDS 3.0-Detecting Genetic and Drug Class Variant Impact on Conserved Protein Binding Dynamics

The application of statistical methods to comparatively framed questions about the molecular dynamics (MD) of proteins can potentially enable investigations of biomolecular function beyond the current sequence and structural methods in bioinformatics. However, the chaotic behavior in single MD trajectories requires statistical inference that is derived from large ensembles of simulations representing the comparative functional states of a protein under investigation. Meaningful interpretation of such complex forms of big data poses serious challenges to users of MD. Here, we announce Detecting Relative Outlier Impacts from Molecular Dynamic Simulation (DROIDS) 3.0, a method and software package for comparative protein dynamics that includes maxDemon 1.0, a multimethod machine learning application that trains on large ensemble comparisons of concerted protein motions in opposing functional states generated by DROIDS and deploys learned classifications of these states onto newly generated MD simulations. Local canonical correlations in learning patterns generated from independent, yet identically prepared, MD validation runs are used to identify regions of functionally conserved protein dynamics. The subsequent impacts of genetic and/or drug class variants on conserved dynamics can also be analyzed by deploying the classifiers on variant MD simulations and quantifying how often these altered protein systems display opposing functional states. Here, we present several case studies of complex changes in functional protein dynamics caused by temperature, genetic mutation, and binding interactions with nucleic acids and small molecules. We demonstrate that our machine learning algorithm can properly identify regions of functionally conserved dynamics in ubiquitin and TATA-binding protein (TBP). We quantify the impact of genetic variation in TBP and drug class variation targeting the ATP-binding region of Hsp90 on conserved dynamics. We identify regions of conserved dynamics in Hsp90 that connect the ATP binding pocket to other functional regions. We also demonstrate that dynamic impacts of various Hsp90 inhibitors rank accordingly with how closely they mimic natural ATP binding.

In vitro reconstitution of yeast RNA polymerase II transcription initiation with high efficiency

Transcription initiation can be reconstituted from highly purified general transcription factors (GTFs), RNA polymerase II (pol II), and promoter DNA. However, earlier biochemical reconstitution systems had a serious technical limitation, namely very poor initiation efficiency. Due to the poor efficiency of the reaction and trace amounts of proteins involved in the pre-initiation complex (PIC) assembly, detection of transcription and PIC formation was only possible by the synthesis of a radiolabeled transcript and by immunoblotting for PIC components on templates. Here we describe a transcription system that is capable of initiating transcription with >90% efficiency of template usage using homogeneous, active yeast components including TFIIA, TFIIB, TBP, TFIIE, TFIIF, TFIIH, Sub1, and pol II. The abundant specifically assembled PICs on promoter DNA can be separated from free general transcription factors (GTFs) and pol II by density gradient sedimentation, irrespective of the length of promoter DNA. The system is robust, and can be modified to accommodate many other transcription factors, and the resulting complexes can be analyzed by SDS-PAGE followed by Coomassie Blue staining. This technical advance now paves the way to conduct definitive biochemical and structural studies of the complete process of pol II initiation from the PIC, through promoter escape, and finally to productive elongation.

The Structures of Eukaryotic Transcription Pre-initiation Complexes and Their Functional Implications

Transcription is a highly regulated process that supplies living cells with coding and non-coding RNA molecules. Failure to properly regulate transcription is associated with human pathologies, including cancers. RNA polymerase II is the enzyme complex that synthesizes messenger RNAs that are then translated into proteins. In spite of its complexity, RNA polymerase requires a plethora of general transcription factors to be recruited to the transcription start site as part of a large transcription pre-initiation complex, and to help it gain access to the transcribed strand of the DNA. This chapter reviews the structure and function of these eukaryotic transcription pre-initiation complexes, with a particular emphasis on two of its constituents, the multisubunit complexes TFIID and TFIIH. We also compare the overall architecture of the RNA polymerase II pre-initiation complex with those of RNA polymerases I and III, involved in transcription of ribosomal RNA and non-coding RNAs such as tRNAs and snRNAs, and discuss the general, conserved features that are applicable to all eukaryotic RNA polymerase systems.

Analyzing the Dynamics of Single TBP-DNA-NC2 Complexes Using Hidden Markov Models

Single-pair Förster resonance energy transfer (spFRET) has become an important tool for investigating conformational dynamics in biological systems. To extract dynamic information from the spFRET traces measured with total internal reflection fluorescence microscopy, we extended the hidden Markov model (HMM) approach. In our extended HMM analysis, we incorporated the photon-shot noise from camera-based systems into the HMM. Thus, the variance in Förster resonance energy transfer (FRET) efficiency of the various states, which is typically a fitted parameter, is explicitly included in the analysis estimated from the number of detected photons. It is also possible to include an additional broadening of the FRET state, which would then only reflect the inherent flexibility of the dynamic biological systems. This approach is useful when comparing the dynamics of individual molecules for which the total intensities vary significantly. We used spFRET with the extended HMM analysis to investigate the dynamics of TATA-box-binding protein (TBP) on promoter DNA in the presence of negative cofactor 2 (NC2). We compared the dynamics of two promoters as well as DNAs of different length and labeling location. For the adenovirus major late promoter, four FRET states were observed; three states correspond to different conformations of the DNA in the TBP-DNA-NC2 complex and a four-state model in which the complex has shifted along the DNA. The HMM analysis revealed that the states are connected via a linear, four-well model. For the H2B promoter, more complex dynamics were observed. By clustering the FRET states detected with the HMM analysis, we could compare the general dynamics observed for the two promoter sequences. We observed that the dynamics from a stretched DNA conformation to a bent conformation for the two promoters were similar, whereas the bent conformation of the TBP-DNA-NC2 complex for the H2B promoter is approximately three times more stable than for the adenovirus major late promoter.

Insights into an Extensively Fragmented Eukaryotic Genome: De Novo Genome Sequencing of the Multinuclear Ciliate Uroleptopsis citrina

Ciliated protists are a large group of single-celled eukaryotes with separate germline and somatic nuclei in each cell. The somatic genome is developed from the zygotic nucleus through a series of chromosomal rearrangements, including fragmentation, DNA elimination, de novo telomere addition, and DNA amplification. This unique feature makes them perfect models for research in genome biology and evolution. However, genomic research of ciliates has been limited to a few species, owing to problems with DNA contamination and obstacles in cultivation. Here, we introduce a method combining telomere-primer PCR amplification and high-throughput sequencing, which can reduce DNA contamination and obtain genomic data efficiently. Based on this method, we report a draft somatic genome of a multimacronuclear ciliate, Uroleptopsis citrina. 1) The telomeric sequence in U. citrina is confirmed to be C4A4C4A4C4 by directly blunt-end cloning. 2) Genomic analysis of the resulting chromosomes shows a "one-gene one-chromosome" pattern, with a small number of multiple-gene chromosomes. 3) Amino acid usage is analyzed, and reassignment of stop codons is confirmed. 4) Chromosomal analysis shows an obvious asymmetrical GC skew and high bias between A and T in the subtelomeric regions of the sense-strand, with the detection of an 11-bp high AT motif region in the 3' subtelomeric region. 5) The subtelomeric sequence also has an obvious 40 nt strand oscillation of nucleotide ratio. 6) In the 5' subtelomeric region of the coding strand, the distribution of potential TATA-box regions is illustrated, which accumulate between 30 and 50 nt. This work provides a valuable reference for genomic research and furthers our understanding of the dynamic nature of unicellular eukaryotic genomes.

Strong deformations of DNA: Effect on the persistence length

Extreme deformations of the DNA double helix attracted a lot of attention during the past decades. Particularly, the determination of the persistence length of DNA with extreme local disruptions, or kinks, has become a crucial problem in the studies of many important biological processes. In this paper we review an approach to calculate the persistence length of the double helix by taking into account the formation of kinks of arbitrary configuration. The reviewed approach improves the Kratky-Porod model to determine the type and nature of kinks that occur in the double helix, by measuring a reduction of the persistence length of the kinkable DNA.

Reconstruction of Atomistic Structures from Coarse-Grained Models for Protein-DNA Complexes

While coarse-grained (CG) simulations have widely been used to accelerate structure sampling of large biomolecular complexes, they are unavoidably less accurate and thus the reconstruction of all-atom (AA) structures and the subsequent refinement is desirable. In this study we developed an efficient method to reconstruct AA structures from sampled CG protein-DNA complex models, which attempts to model the protein-DNA interface accurately. First we developed a method to reconstruct atomic details of DNA structures from a three-site per nucleotide CG model, which uses a DNA fragment library. Next, for the protein-DNA interface, we referred to the side chain orientations in the known structure of the target interface when available. The other parts are modeled by existing tools. We confirmed the accuracy of the protocol in various aspects including the structure deviation in the self-reproduction, the base pair reproducibility, atomic contacts at the protein-DNA interface, and feasibility of the posterior AA simulations.

Probing the role of intercalating protein sidechains for kink formation in DNA

Protein binding can induce DNA kinks, which are for example important to enhance the specificity of the interaction and to facilitate the assembly of multi protein complexes. The respective proteins frequently exhibit amino acid sidechains that intercalate between the DNA base steps at the site of the kink. However, on a molecular level there is only little information available about the role of individual sidechains for kink formation. To unravel structural principles of protein-induced DNA kinking we have performed molecular dynamics (MD) simulations of five complexes that varied in their architecture, function, and identity of intercalated residues. Simulations were performed for the DNA complexes of wildtype proteins (Sac7d, Sox-4, CcpA, TFAM, TBP) and for mutants, in which the intercalating residues were individually or combined replaced by alanine. The work revealed that for systems with multiple intercalated residues, not all of them are necessarily required for kink formation. In some complexes (Sox-4, TBP), one of the residues proved to be essential for kink formation, whereas the second residue has only a very small effect on the magnitude of the kink. In other systems (e.g. Sac7d) each of the intercalated residues proved to be individually capable of conferring a strong kink suggesting a partially redundant role of the intercalating residues. Mutation of the key residues responsible for kinking either resulted in stable complexes with reduced kink angles or caused conformational instability as evidenced by a shift of the kink to an adjacent base step. Thus, MD simulations can help to identify the role of individual inserted residues for kinking, which is not readily apparent from an inspection of the static structures. This information might be helpful for understanding protein-DNA interactions in more detail and for designing proteins with altered DNA binding properties in the future.

Decomposing protein-DNA binding and recognition using simplified protein models

We analyze the role of different physicochemical factors in protein/DNA binding and recognition by comparing the results from all-atom molecular dynamics simulations with simulations using simplified protein models. These models enable us to separate the role of specific amino acid side chains, formal amino acid charges and hydrogen bonding from the effects of the low-dielectric volume occupied by the protein. Comparisons are made on the basis of the conformation of DNA after protein binding, the ionic distribution around the complex and the sequence specificity. The results for four transcription factors, binding in either the minor or major grooves of DNA, show that the protein volume and formal charges, with one exception, play a predominant role in binding. Adding hydrogen bonding and a very small number of key amino acid side chains at the all-atom level yields results in DNA conformations and sequence recognition close to those seen in the reference all-atom simulations.

Architecture of TAF11/TAF13/TBP complex suggests novel regulation properties of general transcription factor TFIID

General transcription factor TFIID is a key component of RNA polymerase II transcription initiation. Human TFIID is a megadalton-sized complex comprising TATA-binding protein (TBP) and 13 TBP-associated factors (TAFs). TBP binds to core promoter DNA, recognizing the TATA-box. We identified a ternary complex formed by TBP and the histone fold (HF) domain-containing TFIID subunits TAF11 and TAF13. We demonstrate that TAF11/TAF13 competes for TBP binding with TATA-box DNA, and also with the N-terminal domain of TAF1 previously implicated in TATA-box mimicry. In an integrative approach combining crystal coordinates, biochemical analyses and data from cross-linking mass-spectrometry (CLMS), we determine the architecture of the TAF11/TAF13/TBP complex, revealing TAF11/TAF13 interaction with the DNA binding surface of TBP. We identify a highly conserved C-terminal TBP-interaction domain (CTID) in TAF13, which is essential for supporting cell growth. Our results thus have implications for cellular TFIID assembly and suggest a novel regulatory state for TFIID function.

[KCl] Dependence of B-DNA Groove Bending Anisotropy

The energetics of B-DNA bending toward the major and minor grooves were quantified by free energy simulations at four different KCl concentrations. Increased [KCl] led to more flexible DNA, with persistence lengths that agreed well with experimental values. At all salt concentrations, major groove bending was preferred, although preferences for major and minor groove bending were similar for the A-tract containing sequence. Since the phosphate repulsions and DNA internal energy favored minor groove bending, the preference for major groove bending was thought to originate from differences in solvation. Water in the minor groove was tighter bound than water in the major groove, and harder to displace than major groove water, which favored the compression of the major groove upon bending. Higher [KCl] decreased the persistence length for both major and minor groove bending but did not greatly affect the free energy spacing between the minor and major groove bending curves. For sequences without A-tracts, salt affected major and minor bending to nearly the same degree, and did not change the preference for major groove bending. For the A-tract containing sequence, an increase in salt concentration decreased the already small energetic difference between major and minor groove bending. Since salts did not significantly affect the relative differences in bending energetics and hydration, it is likely that the increased bending flexibilities upon salt increase are simply due to screening.

Structural visualization of the p53/RNA polymerase II assembly

Singh et al. dissected the human p53/Pol II interaction via single-particle cryo-electron microscopy, structural docking, and biochemical analyses. These findings indicate that p53 may structurally regulate DNA-binding functions of Pol II via the clamp domain, thereby providing insights into p53-regulated Pol II transcription.

The master tumor suppressor p53 activates transcription in response to various cellular stresses in part by facilitating recruitment of the transcription machinery to DNA. Recent studies have documented a direct yet poorly characterized interaction between p53 and RNA polymerase II (Pol II). Therefore, we dissected the human p53/Pol II interaction via single-particle cryo-electron microscopy, structural docking, and biochemical analyses. This study reveals that p53 binds Pol II via the Rpb1 and Rpb2 subunits, bridging the DNA-binding cleft of Pol II proximal to the upstream DNA entry site. In addition, the key DNA-binding surface of p53, frequently disrupted in various cancers, remains exposed within the assembly. Furthermore, the p53/Pol II cocomplex displays a closed conformation as defined by the position of the Pol II clamp domain. Notably, the interaction of p53 and Pol II leads to increased Pol II elongation activity. These findings indicate that p53 may structurally regulate DNA-binding functions of Pol II via the clamp domain, thereby providing insights into p53-regulated Pol II transcription.

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.

A real-time study of the interaction of TBP with a TATA box-containing duplex identical to an ancestral or minor allele of human gene LEP or TPI

It is known that only a single-nucleotide substitution (SNP: a single nucleotide polymorphism) in the sequence of a TATA box can influence the affinity of the interaction of TBP with the TATA box and contribute to the pathogenesis of complex hereditary human diseases and sometimes may be a cause of monogenic diseases (for instance, β-thalassemia). In the present work, we studied the interaction of human TBP with a double-stranded oligodeoxyribonucleotide (ODN) 15 or 26 bp long identical to a TATA box of promoters of a real-life human gene, TPI or LEP, and labeled with fluorophores TAMRA and FAM. To analyze the interaction of TBP with a TATA box of an ancestral or minor allele (SNP in the TATA box) in real time, we used the stopped-flow method with detection of a Förster resonance energy transfer (FRET) signal. The nature of the resulting kinetic curves reflecting changes in the FRET signal (and therefore of DNA conformation during the interaction with TBP) pointed to a multistage mechanism of the formation of the TBP complex with the TATA-containing ODN. The results showed that with the increasing concentration and length of the ODN, heterogeneity of conformational changes (taking place during the first second of the interaction with TBP) in DNA also increases. In contrast to the initial nonspecific interaction, the subsequent phases strictly depend on TBP concentration: at the TBP:ODN ratio of 10:1, the velocity of change of the FRET signal increases approximately 100-fold.

Protein-DNA interfaces: a molecular dynamics analysis of time-dependent recognition processes for three transcription factors

We have studied the dynamics of three transcription factor-DNA complexes using all-atom, microsecond-scale MD simulations. In each case, the salt bridges and hydrogen bond interactions formed at the protein-DNA interface are found to be dynamic, with lifetimes typically in the range of tens to hundreds of picoseconds, although some interactions, notably those involving specific binding to DNA bases, can be a hundred times longer lived. Depending on the complex studied, this dynamics may or may not lead to the existence of distinct conformational substates. Using a sequence threading technique, it has been possible to determine whether DNA sequence recognition is sensitive or not to such conformational changes, and, in one case, to show that recognition appears to be locally dependent on protein-mediated cation distributions.

Anisotropy of B-DNA Groove Bending

DNA bending is critical for DNA packaging, recognition, and repair, and occurs toward either the major or the minor groove. The anisotropy of B-DNA groove bending was quantified for eight DNA sequences by free energy simulations employing a novel reaction coordinate. The simulations show that bending toward the major groove is preferred for non-A-tracts while the A-tract has a high tendency of bending toward the minor groove. Persistence lengths were generally larger for bending toward the minor groove, which is thought to originate from differences in groove hydration. While this difference in stiffness is one of the factors determining the overall preference of bending direction, the dominant contribution is shown to be a free energy offset between major and minor groove bending. The data suggests that, for the A-tract, this offset is largely determined by inherent structural properties, while differences in groove hydration play a large role for non-A-tracts. By quantifying the energetics of DNA groove bending and rationalizing the origins of the anisotropy, the calculations provide important new insights into a key biological process.

Uncovering ancient transcription systems with a novel evolutionary indicator

TBP and TFIIB are evolutionarily conserved transcription initiation factors in archaea and eukaryotes. Information about their ancestral genes would be expected to provide insight into the origin of the RNA polymerase II-type transcription apparatus. In obtaining such information, the nucleotide sequences of current genes of both archaea and eukaryotes should be included in the analysis. However, the present methods of evolutionary analysis require that a subset of the genes should be excluded as an outer group. To overcome this limitation, we propose an innovative concept for evolutionary analysis that does not require an outer group. This approach utilizes the similarity in intramolecular direct repeats present in TBP and TFIIB as an evolutionary measure revealing the degree of similarity between the present offspring genes and their ancestors. Information on the properties of the ancestors and the order of emergence of TBP and TFIIB was also revealed. These findings imply that, for evolutionarily early transcription systems billions of years ago, interaction of RNA polymerase II with transcription initiation factors and the regulation of its enzymatic activity was required prior to the accurate positioning of the enzyme. Our approach provides a new way to discuss mechanistic and system evolution in a quantitative manner.

Zooming in on Transcription Preinitiation

Class II gene transcription commences with the assembly of the Preinitiation Complex (PIC) from a plethora of proteins and protein assemblies in the nucleus, including the General Transcription Factors (GTFs), RNA polymerase II (RNA pol II), co-activators, co-repressors, and more. TFIID, a megadalton-sized multiprotein complex comprising 20 subunits, is among the first GTFs to bind the core promoter. TFIID assists in nucleating PIC formation, completed by binding of further factors in a highly regulated stepwise fashion. Recent results indicate that TFIID itself is built from distinct preformed submodules, which reside in the nucleus but also in the cytosol of cells. Here, we highlight recent insights in transcription factor assembly and the regulation of transcription preinitiation.

•

Architectural models of human and yeast PIC were proposed.

•

Mediator core–ITC complex structure reveals novel interactions.

•

TFIID submodule residing in the cytoplasm has been discovered.

•

Complex assembly emerges as key concept in transcription regulation.

Novel alleles of the VERNALIZATION1 genes in wheat are associated with modulation of DNA curvature and flexibility in the promoter region

BACKGROUND

In wheat, the vernalization requirement is mainly controlled by the VRN genes. Different species of hexaploid and tetraploid wheat are widely used as genetic source for new mutant variants and alleles for fundamental investigations and practical breeding programs. In this study, VRN-A1 and VRN-B1 were analysed for 178 accessions representing six tetraploid wheat species (Triticum dicoccoides, T. dicoccum, T. turgidum, T. polonicum, T. carthlicum, T. durum) and five hexaploid species (T. compactum, T. sphaerococcum, T. spelta, T. macha, T. vavilovii).

RESULTS

Novel allelic variants in the promoter region of VRN-A1 and VRN-B1 were identified based on the change in curvature and flexibility of the DNA molecules. The new variants of VRN-A1 (designated as Vrn-A1a.2, Vrn-A1b.2 - Vrn-A1b.6 and Vrn-A1i) were found to be widely distributed in hexaploid and tetraploid wheat, and in fact were predominant over the known VRN-A1 alleles. The greatest diversity of the new variants of VRN-B1 (designated as VRN-B1.f, VRN-B1.s and VRN-B1.m) was found in the tetraploid and some hexaploid wheat species. For the first time, minor differences within the sequence motif known as the VRN-box of VRN1 were correlated with wheat growth habit. Thus, vrn-A1b.3 and vrn-A1b.4 were revealed in winter wheat in contrast to Vrn-A1b.2, Vrn-A1b.5, Vrn-A1b.6 and Vrn-A1i. It was found that single nucleotide mutation in the VRN-box can influence the vernalization requirement and growth habit of wheat. Our data suggest that both the A-tract and C-rich segment within the VRN-box contribute to its functionality, and provide a new view of the hypothesised role of the VRN-box in regulating transcription of the VRN1 genes. Specifically, it is proposed that combination of mutations in this region can modulate vernalization sensitivity and flowering time of wheat.

CONCLUSIONS

New allelic variants of the VRN-A1 and VRN-B1 genes were identified in hexaploid and tetraploid wheat. Mutations in A-tract and C-rich segments within the VRN-box of VRN-A1 are associated with modulation of the vernalization requirement and flowering time. New allelic variants will be useful in fundamental investigations into the regulation of VRN1 expression, and provide a valuable genetic resource for practical breeding of wheat.

A reversible B–A transition of DNA duplexes induced by synthetic cationic copolymers

Reversible B–A transitions of DNA duplexes were induced by synthetic cationic and anionic polymers.