Genes encoding toxin–antitoxin (TA) systems are near ubiquitous in bacterial genomes and they play key roles in important aspects of bacterial physiology, including genomic stability, formation of persister cells under antibiotic stress, and resistance to phage infection. The CptIN locus from Eubacterium rectale is a member of the recently-discovered Type III class of TA systems, defined by a protein toxin suppressed by direct interaction with a structured RNA antitoxin. Here, we present the crystal structure of the CptIN protein–RNA complex to 2.2 Å resolution. The structure reveals a new heterotetrameric quaternary organization for the Type III TA class, and the RNA antitoxin bears a novel structural feature of an extended A-twist motif within the pseudoknot fold. The retention of a conserved ribonuclease active site as well as traits normally associated with TA systems, such as plasmid maintenance, implicates a wider functional role for Type III TA systems. We present evidence for the co-variation of the Type III component pair, highlighting a distinctive evolutionary process in which an enzyme and its substrate co-evolve.

Thymine DNA Glycosylase (TDG) performs essential functions in maintaining genetic integrity and epigenetic regulation. Initiating base excision repair, TDG removes thymine from mutagenic G·T mispairs caused by 5-methylcytosine (mC) deamination and other lesions including uracil (U) and 5-hydroxymethyluracil (hmU). In DNA demethylation, TDG excises 5-formylcytosine (fC) and 5-carboxylcytosine (caC), which are generated from mC by Tet (ten–eleven translocation) enzymes. Using improved crystallization conditions, we solved high-resolution (up to 1.45 Å) structures of TDG enzyme–product complexes generated from substrates including G·U, G·T, G·hmU, G·fC and G·caC. The structures reveal many new features, including key water-mediated enzyme–substrate interactions. Together with nuclear magnetic resonance experiments, the structures demonstrate that TDG releases the excised base from its tight product complex with abasic DNA, contrary to previous reports. Moreover, DNA-free TDG exhibits no significant binding to free nucleobases (U, T, hmU), indicating a Kd >> 10 mM. The structures reveal a solvent-filled channel to the active site, which might facilitate dissociation of the excised base and enable caC excision, which involves solvent-mediated acid catalysis. Dissociation of the excised base allows TDG to bind the beta rather than the alpha anomer of the abasic sugar, which might stabilize the enzyme–product complex.

The hexameric Minichromosome Maintenance (MCM) protein complex forms a ring that unwinds DNA at the replication fork in eukaryotes and archaea. Our recent crystal structure of an archaeal MCM N-terminal domain bound to single-stranded DNA (ssDNA) revealed ssDNA associating across tight subunit interfaces but not at the loose interfaces, indicating that DNA-binding is governed not only by the DNA-binding residues of the subunits (MCM ssDNA-binding motif, MSSB) but also by the relative orientation of the subunits. We now extend these findings by showing that DNA-binding by the MCM N-terminal domain of the archaeal organism Pyrococcus furiosus occurs specifically in the hexameric oligomeric form. We show that mutants defective for hexamerization are defective in binding ssDNA despite retaining all the residues observed to interact with ssDNA in the crystal structure. One mutation that exhibits severely defective hexamerization and ssDNA-binding is at a conserved phenylalanine that aligns with the mouse Mcm4(Chaos3) mutation associated with chromosomal instability, cancer, and decreased intersubunit association.

Fluorogenic RNAs that are based on the complex formed by 3,5-difluoro-4-hydroxybenzylidene imidazolinone (DFHBI) derivatives and the RNA aptamer named Spinach were used to engineer a new generation of in vitro and in vivo sensors for bioanalytics. With the resolved crystal structure of the RNA/small molecule complex, the engineering map becomes available, but comprehensive information regarding the thermodynamic profile of the molecule is missing. Here, we reconstructed the full thermodynamic binding and stability landscapes between DFHBI and a truncated sequence of first-generation Spinach. For this purpose, we established a systematic screening procedure for single- and double-point mutations on a microfluidic large-scale integrated chip platform for 87-nt long RNAs. The thermodynamic profile with single base resolution was used to engineer an improved fluorogenic spinach generation via a directed rather than evolutional approach.

Histone post-translational modifications are key contributors to chromatin structure and function, and participate in the maintenance of genome stability. Understanding the establishment and maintenance of these marks, along with their misregulation in pathologies is thus a major focus in the field. While we have learned a great deal about the enzymes regulating histone modifications on nucleosomal histones, much less is known about the mechanisms establishing modifications on soluble newly synthesized histones. This includes methylation of lysine 9 on histone H3 (H3K9), a mark that primes the formation of heterochromatin, a critical chromatin landmark for genome stability. Here, we report that H3K9 mono- and dimethylation is imposed during translation by the methyltransferase SetDB1. We discuss the importance of these results in the context of heterochromatin establishment and maintenance and new therapeutic opportunities in pathologies where heterochromatin is perturbed.

We have recently shown that the incorporation of modified nucleotides such as 5-N-carboxamide-deoxyuridines into random nucleic acid libraries improves success rates in SELEX experiments and facilitates the identification of ligands with slow off-rates. Here we report the impact of these modifications on the thermodynamic stability of both duplexes and intramolecular ‘single-stranded’ structures. Within duplexes, large, hydrophobic naphthyl groups were destabilizing relative to the all natural DNA duplex, while the hydrophilic groups exhibited somewhat improved duplex stability. All of the significant changes in stability were driven by opposing contributions from the enthalpic and entropic terms. In contrast, both benzyl and naphthyl modifications stabilized intramolecular single-stranded structures relative to their natural DNA analogs, consistent with the notion that intramolecular folding allows formation of novel, stabilizing hydrophobic interactions. Imino proton NMR data provided evidence that elements of the folded structure form at temperatures well below the Tm, with a melting transition that is distinctly less cooperative when compared to duplex DNA. Although there are no data to suggest that the unmodified DNA sequences fold into structures similar to their modified analogs, this still represents clear evidence that these modifications impart thermodynamic stability to the folded structure not achievable with unmodified DNA.

The in vivo potency of antisense oligonucleotides (ASO) has been significantly increased by reducing their length to 8–15 nucleotides and by the incorporation of high affinity RNA binders such as 2', 4'-bridged nucleic acids (also known as locked nucleic acid or LNA, and 2',4'-constrained ethyl [cET]). We now report the development of a novel ASO design in which such short ASO monomers to one or more targets are co-synthesized as homo- or heterodimers or multimers via phosphodiester linkers that are stable in plasma, but cleaved inside cells, releasing the active ASO monomers. Compared to current ASOs, these multimers and multi-targeting oligonucleotides (MTOs) provide increased plasma protein binding and biodistribution to liver, and increased in vivo efficacy against single or multiple targets with a single construct. In vivo, MTOs synthesized in both RNase H-activating and steric-blocking oligonucleotide designs provide 4–5-fold increased potency and 2-fold increased efficacy, suggesting broad therapeutic applications.

The mutational pattern for the TP53 tumour suppressor gene in lung tumours differs to other cancer types by having a higher frequency of G:C>T:A transversions. The aetiology of this differing mutation pattern is still unknown. Benzo[a]pyrene,diol epoxide (BPDE) is a potent cigarette smoke carcinogen that forms guanine adducts at TP53 CpG mutation hotspot sites including codons 157, 158, 245, 248 and 273. We performed molecular modelling of BPDE-adducted TP53 duplex sequences to determine the degree of local distortion caused by adducts which could influence the ability of nucleotide excision repair. We show that BPDE adducted codon 157 has greater structural distortion than other TP53 G:C>T:A hotspot sites and that sequence context more distal to adjacent bases must influence local distortion. Using TP53 trinucleotide mutation signatures for lung cancer in smokers and non-smokers we further show that codons 157 and 273 have the highest mutation probability in smokers. Combining this information with adduct structural data we predict that G:C>T:A mutations at codon 157 in lung tumours of smokers are predominantly caused by BPDE. Our results provide insight into how different DNA sequence contexts show variability in DNA distortion at mutagen adduct sites that could compromise DNA repair at well characterized cancer related mutation hotspots.

Development of an accurate protein–DNA recognition code that can predict DNA specificity from protein sequence is a central problem in biology. C2H2 zinc fingers constitute by far the largest family of DNA binding domains and their binding specificity has been studied intensively. However, despite decades of research, accurate prediction of DNA specificity remains elusive. A major obstacle is thought to be the inability of current methods to account for the influence of neighbouring domains. Here we show that this problem can be addressed using a structural approach: we build structural models for all C2H2-ZF–DNA complexes with known binding motifs and find six distinct binding modes. Each mode changes the orientation of specificity residues with respect to the DNA, thereby modulating base preference. Most importantly, the structural analysis shows that residues at the domain interface strongly and predictably influence the binding mode, and hence specificity. Accounting for predicted binding mode significantly improves prediction accuracy of predicted motifs. This new insight into the fundamental behaviour of C2H2-ZFs has implications for both improving the prediction of natural zinc finger-binding sites, and for prioritizing further experiments to complete the code. It also provides a new design feature for zinc finger engineering.

Here we describe our study of miRNA isoforms (isomiRs) in breast cancer (BRCA) and normal breast data sets from the Cancer Genome Atlas (TCGA) repository. We report that the full isomiR profiles, from both known and novel human-specific miRNA loci, are particularly rich in information and can distinguish tumor from normal tissue much better than the archetype miRNAs. IsomiR expression is also dependent on the patient's race, exemplified by miR-183-5p, several isomiRs of which are upregulated in triple negative BRCA in white but not black women. Additionally, we find that an isomiR's 5' endpoint and length, but not the genomic origin, are key determinants of the regulation of its expression. Overexpression of distinct miR-183-5p isomiRs in MDA-MB-231 cells followed by microarray analysis revealed that each isomiR has a distinct impact on the cellular transcriptome. Parallel integrative analysis of mRNA expression from BRCA data sets of the TCGA repository demonstrated that isomiRs can distinguish between the luminal A and luminal B subtypes and explain in more depth the molecular differences between them than the archetype molecules. In conclusion, our findings provide evidence that post-transcriptional studies of BRCA will benefit from transcending the one-locus-one-miRNA paradigm and taking into account all isoforms from each miRNA locus as well as the patient's race.

DNA-binding proteins (DBPs) rapidly search and specifically bind to their target sites on genomic DNA in order to trigger many cellular regulatory processes. It has been suggested that the facilitation of search dynamics is achieved by combining 3D diffusion with one-dimensional sliding and hopping dynamics of interacting proteins. Although, recent studies have advanced the knowledge of molecular determinants that affect one-dimensional search efficiency, the role of DNA molecule is poorly understood. In this study, by using coarse-grained simulations, we propose that dynamics of DNA molecule and its degree of confinement due to cellular crowding concertedly regulate its groove geometry and modulate the inter-communication with DBPs. Under weak confinement, DNA dynamics promotes many short, rotation-decoupled sliding events interspersed by hopping dynamics. While this results in faster 1D diffusion, associated probability of missing targets by jumping over them increases. In contrast, strong confinement favours rotation-coupled sliding to locate targets but lacks structural flexibility to achieve desired specificity. By testing under physiological crowding, our study provides a plausible mechanism on how DNA molecule may help in maintaining an optimal balance between fast hopping and rotation-coupled sliding dynamics, to locate target sites rapidly and form specific complexes precisely.

Recent studies have revealed significant roles of RNA structure in almost every step of RNA processing, including transcription, splicing, transport and translation. RNase footprint sequencing (RNase-seq) has emerged to dissect RNA structures at the genome scale. However, it remains challenging to analyze RNase-seq data because of the issues of signal sparsity, variability and correlations among various RNases. We present a probabilistic framework, joint Poisson-gamma mixture (JPGM), for integrative modeling of multiple RNase-seq profiles. Combining JPGM with hidden Markov model allows genome-wide inference of RNA structures. We apply the joint modeling approach for inferring base pairing states on simulated data sets and RNase-seq profiles of the double-strand specific RNase V1 and single-strand specific RNase S1 in yeast. We demonstrate that joint analysis of V1 and S1 profiles outputs interpretable RNA structure states, while approaches that analyze each profile separately do not. The joint modeling approach predicts the structure states of all nucleotides in 3196 transcripts of yeast without compromising accuracy, while the simple thresholding approach misses 43% of the nucleotides. Furthermore, the posterior probabilities outputted by our model are able to resolve the structural ambiguity of 300 000 nucleotides with overlapping V1 and S1 cleavage sites. Our model also generates RNA accessibilities, which are associated with three-dimensional conformations.

The transcriptional programs of ectothermic teleosts are directly influenced by water temperature. However, the cis- and trans-factors governing cold responses are not well characterized. We profiled transcriptional changes in eight zebrafish tissues exposed to mildly and severely cold temperatures using RNA-Seq. A total of 1943 differentially expressed genes (DEGs) were identified, from which 34 clusters representing distinct tissue and temperature response expression patterns were derived using the k-means fuzzy clustering algorithm. The promoter regions of the clustered DEGs that demonstrated strong co-regulation were analysed for enriched cis-regulatory elements with a motif discovery program, DREME. Seventeen motifs, ten known and seven novel, were identified, which covered 23% of the DEGs. Two motifs predicted to be the binding sites for the transcription factors Bcl6 and Jun, respectively, were chosen for experimental verification, and they demonstrated the expected cold-induced and cold-repressed patterns of gene regulation. Protein interaction modeling of the network components followed by experimental validation suggested that Jun physically interacts with Bcl6 and might be a hub factor that orchestrates the cold response in zebrafish. Thus, the methodology used and the regulatory networks uncovered in this study provide a foundation for exploring the mechanisms of cold adaptation in teleosts.

Mediator is a large multiprotein complex conserved in all eukaryotes, which has a crucial coregulator function in transcription by RNA polymerase II (Pol II). However, the molecular mechanisms of its action in vivo remain to be understood. Med17 is an essential and central component of the Mediator head module. In this work, we utilised our large collection of conditional temperature-sensitive med17 mutants to investigate Mediator's role in coordinating preinitiation complex (PIC) formation in vivo at the genome level after a transfer to a non-permissive temperature for 45 minutes. The effect of a yeast mutation proposed to be equivalent to the human Med17-L371P responsible for infantile cerebral atrophy was also analyzed. The ChIP-seq results demonstrate that med17 mutations differentially affected the global presence of several PIC components including Mediator, TBP, TFIIH modules and Pol II. Our data show that Mediator stabilizes TFIIK kinase and TFIIH core modules independently, suggesting that the recruitment or the stability of TFIIH modules is regulated independently on yeast genome. We demonstrate that Mediator selectively contributes to TBP recruitment or stabilization to chromatin. This study provides an extensive genome-wide view of Mediator's role in PIC formation, suggesting that Mediator coordinates multiple steps of a PIC assembly pathway.

Staphylococcus aureus pathogenesis is directed by regulatory proteins and RNAs. We report the case of an RNA attenuating virulence and host uptake, possibly to sustain commensalism. A S. aureus sRNA, SprC (srn_3610), reduced virulence and bacterial loads in a mouse infection model. S. aureus deleted for sprC became more virulent and increased bacterial dissemination in colonized animals. Conversely, inducing SprC expression lowered virulence and the bacterial load. Without sprC, S. aureus phagocytosis by monocytes and macrophages was higher, whereas bacteria were internalized at lower yields when SprC expression was stimulated. Without sprC, higher internalization led to a greater number of extracellular bacteria, facilitating colonization. SprC expression decreased after phagocytosis, concurring with the facilitated growth of bacteria lacking the sRNA in the presence of an oxidant. The major staphylococcal autolysin facilitates S. aureus uptake by human phagocytes. ATL proved to be negatively regulated by SprC. The SprC domains involved in pairing with atl mRNA were analyzed. The addition of ATL reduced phagocytosis of bacteria lacking sprC with no effects on wild-type bacterial uptake, implying that SprC influences phagocytosis, at least in part, by controlling ATL. Since the control of SprC on ATL was modest, other factors must contribute to atl regulation.

The gene is encoded in broad-host range and low-copy plasmids. It is genetically linked to antibiotic resistance genes of the major human pathogens of phylum Firmicutes. The homodimeric forms of (2) coordinate the plasmid copy number control, faithful partition (2 and 2) and better-than-random segregation (2) systems. The promoter (P) of the operon (P) transiently interacts with 2. Adding 2 facilitates the formation of stable 2·P complexes. Here we show that limiting 2 interacts with the N-terminal domain of the β’ subunit of the Bacillus subtilis RNA polymerase (RNAP-A) vegetative holoenzyme. In this way 2 recruits RNAP-A onto P DNA. Partial P occupancy by 2 increases the rate at which RNAP-A complex shifts from its closed (RPC) to open (RPO) form. This shift increases transcription activation. Adding 2 further increases the rate of P transcription initiation, perhaps by stabilizing the 2·P complex. In contrast, full operator occupancy by 2 facilitates RPC formation, but it blocks RPO isomerization and represses P utilization. The stimulation and inhibition of RPO formation is the mechanism whereby 2 mediates copy number fluctuation and stable plasmid segregation. By this mechanism, 2 also indirectly influences the acquisition of antibiotic resistance genes.