CURVES+ web server for analyzing and visualizing the helical, backbone and groove parameters of nucleic acid structures

Structural dynamics of therapeutic nucleic acids with phosphorothioate backbone modifications

Antisense oligonucleotides (ASOs) offer ground-breaking possibilities for selective pharmacological intervention for any gene product-related disease. Therapeutic ASOs contain extensive chemical modifications that improve stability to enzymatic cleavage and modulate binding affinity relative to natural RNA/DNA. Molecular dynamics (MD) simulation can provide valuable insights into how such modifications affect ASO conformational sampling and target binding. However, force field parameters for chemically modified nucleic acids (NAs) are still underdeveloped. To bridge this gap, we developed parameters to allow simulations of ASOs with the widely applied phosphorothioate (PS) backbone modification, and validated these in extensive all-atom MD simulations of relevant PS-modified NA systems representing B-DNA, RNA, and DNA/RNA hybrid duplex structures. Compared to the corresponding natural NAs, single PS substitutions had marginal effects on the ordered DNA/RNA duplex, whereas substantial effects of phosphorothioation were observed in single-stranded RNA and B-DNA, corroborated by the experimentally derived structure data. We find that PS-modified NAs shift between high and low twist states, which could affect target recognition and protein interactions for phosphorothioated oligonucleotides. Furthermore, conformational sampling was markedly altered in the PS-modified ssRNA system compared to that of the natural oligonucleotide, indicating sequence-dependent effects on conformational preference that may in turn influence duplex formation.

Mechanism of Daunomycin Intercalation into DNA from Enhanced Sampling Simulations

Daunomycin is a widely used anticancer drug, yet the mechanism underlying how it binds to DNA remains contested. 469 all-atom trajectories of daunomycin binding to the DNA oligonucleotide d(GCG CAC GTG CGC) were collected using weighted ensemble (WE)-enhanced sampling. Mechanistic insights were revealed through analysis of the ensemble of trajectories. Initially, the binding process involves a ubiquitous hydrogen bond between the DNA backbone and the NH3+ group on daunomycin. During the binding process, most trajectories exhibited similar structural changes to DNA, including DNA base pair rise, bending, and minor groove width changes. Variability within the ensemble of binding trajectories illuminates differences in the orientation of daunomycin as it initially intercalates; around 10% of trajectories needed minimal rearrangement from intercalation to reaching the fully bound configuration, whereas most needed an additional 1-5 ns to rearrange. The results here emphasize the utility of generating an ensemble of trajectories to discern biomolecular binding mechanisms.

Caenorhabditis elegans Dicer acts with the RIG-I-like helicase DRH-1 and RDE-4 to cleave dsRNA

Invertebrates use the endoribonuclease Dicer to cleave viral dsRNA during antiviral defense, while vertebrates use RIG-I-like Receptors (RLRs), which bind viral dsRNA to trigger an interferon response. While some invertebrate Dicers act alone during antiviral defense, Caenorhabditis elegans Dicer acts in a complex with a dsRNA binding protein called RDE-4, and an RLR ortholog called DRH-1. We used biochemical and structural techniques to provide mechanistic insight into how these proteins function together. We found RDE-4 is important for ATP-independent and ATP-dependent cleavage reactions, while helicase domains of both DCR-1 and DRH-1 contribute to ATP-dependent cleavage. DRH-1 plays the dominant role in ATP hydrolysis, and like mammalian RLRs, has an N-terminal domain that functions in autoinhibition. A cryo-EM structure indicates DRH-1 interacts with DCR-1's helicase domain, suggesting this interaction relieves autoinhibition. Our study unravels the mechanistic basis of the collaboration between two helicases from typically distinct innate immune defense pathways.

Molecular choreography: Unveiling the dynamic landscape of type IIA DNA topoisomerases before T-segment passage through all-atom simulations

Type IIA DNA topoisomerases are molecular nanomachines responsible for controlling topological states of DNA molecules. Here, we explore the dynamic landscape of yeast topoisomerase IIA during key stages of its catalytic cycle, focusing in particular on the events preceding the passage of the T-segment. To this end, we generated six configurations of fully catalytic yeast topo IIA, strategically inserted a T-segment into the N-gate in relevant configurations and performed all-atom simulations. The essential motion of topo IIA protein dimer was characterized by rotational gyrating-like movement together with sliding motion within the DNA-gate. Both appear to be inherent properties of the enzyme and an inbuilt feature that allows passage of the T-segment through the cleaved G-segment. Coupled dynamics of the N-gate and DNA-gate residues may be particularly important for controlled and smooth passage of the T-segment and consequently the prevention of DNA double-strand breaks. QTK loop residue Lys367, which interacts with ATP and ADP molecules, is involved in regulating the size and stability of the N-gate. The unveiled features of the simulated configurations provide insights into the catalytic cycle of type IIA topoisomerases and elucidate the molecular choreography governing their ability to modulate the topological states of DNA topology.

Sequence Dependence in Nucleosome Dynamics

The basic packaging unit of eukaryotic chromatin is the nucleosome that contains 145-147 base pair duplex DNA wrapped around an octameric histone protein. While the DNA sequence plays a crucial role in controlling the positioning of the nucleosome, the molecular details behind the interplay between DNA sequence and nucleosome dynamics remain relatively unexplored. This study analyzes this interplay in detail by performing all-atom molecular dynamics simulations of nucleosomes, comparing the human α-satellite palindromic (ASP) and the strong positioning "Widom-601" DNA sequence at time scales of 12 μs. The simulations are performed at salt concentrations 10-20 times higher than physiological salt concentrations to screen the electrostatic interactions and promote unwrapping. These microsecond-long simulations give insight into the molecular-level sequence-dependent events that dictate the pathway of DNA unwrapping. We find that the "ASP" sequence forms a loop around SHL ± 5 for three sets of simulations. Coincident with loop formation is a cooperative increase in contacts with the neighboring N-terminal H2B tail and C-terminal H2A tail and the release of neighboring counterions. We find that the Widom-601 sequence exhibits a strong breathing motion of the nucleic acid ends. Coincident with the breathing motion is the collapse of the full N-terminal H3 tail and formation of an α-helix that interacts with the H3 histone core. We postulate that the dynamics of these histone tails and their modification with post-translational modifications (PTMs) may play a key role in governing this dynamics.

Understanding intercalative modulation of G-rich sequence folding: solution structure of a TINA-conjugated antiparallel DNA triplex

We present here the high-resolution structure of an antiparallel DNA triplex in which a monomer of para-twisted intercalating nucleic acid (para-TINA: (R)-1-O-[4-(1-pyrenylethynyl)phenylmethyl]glycerol) is covalently inserted as a bulge in the third strand of the triplex. TINA is a potent modulator of the hybridization properties of DNA sequences with extremely useful properties when conjugated in G-rich oligonucleotides. The insertion of para-TINA between two guanines of the triplex imparts a high thermal stabilization (ΔTM = 9ºC) to the structure and enhances the quality of NMR spectra by increasing the chemical shift dispersion of proton signals near the TINA location. The structural determination reveals that TINA intercalates between two consecutive triads, causing only local distortions in the structure. The two aromatic moieties of TINA are nearly coplanar, with the phenyl ring intercalating between the flanking guanine bases in the sequence, and the pyrene moiety situated between the Watson-Crick base pair of the two first strands. The precise position of TINA within the triplex structure reveals key TINA-DNA interactions, which explains the high stabilization observed and will aid in the design of new and more efficient binders to DNA.

Flexible structural arrangement and DNA-binding properties of protein p6 from Bacillus subtillis phage φ29

The genome-organizing protein p6 of Bacillus subtilis bacteriophage φ29 plays an essential role in viral development by activating the initiation of DNA replication and participating in the early-to-late transcriptional switch. These activities require the formation of a nucleoprotein complex in which the DNA adopts a right-handed superhelix wrapping around a multimeric p6 scaffold, restraining positive supercoiling and compacting the viral genome. Due to the absence of homologous structures, prior attempts to unveil p6's structural architecture failed. Here, we employed AlphaFold2 to engineer rational p6 constructs yielding crystals for three-dimensional structure determination. Our findings reveal a novel fold adopted by p6 that sheds light on its self-association mechanism and its interaction with DNA. By means of protein-DNA docking and molecular dynamic simulations, we have generated a comprehensive structural model for the nucleoprotein complex that consistently aligns with its established biochemical and thermodynamic parameters. Besides, through analytical ultracentrifugation, we have confirmed the hydrodynamic properties of the nucleocomplex, further validating in solution our proposed model. Importantly, the disclosed structure not only provides a highly accurate explanation for previously experimental data accumulated over decades, but also enhances our holistic understanding of the structural and functional attributes of protein p6 during φ29 infection.

Structure of the WYL-domain containing transcription activator, DriD, in complex with ssDNA effector and DNA target site

Transcription regulators play central roles in orchestrating responses to changing environmental conditions. Recently the Caulobacter crescentus transcription activator DriD, which belongs to the newly defined WYL-domain family, was shown to regulate DNA damage responses independent of the canonical SOS pathway. However, the molecular mechanisms by which DriD and other WYL-regulators sense environmental signals and recognize DNA are not well understood. We showed DriD DNA-binding is triggered by its interaction with ssDNA, which is produced during DNA damage. Here we describe the structure of the full-length C. crescentus DriD bound to both target DNA and effector ssDNA. DriD consists of an N-terminal winged-HTH (wHTH) domain, linker region, three-helix bundle, WYL-domain and C-terminal WCX-dimer domain. Strikingly, DriD binds DNA using a novel, asymmetric DNA-binding mechanism that results from different conformations adopted by the linker. Although the linker does not touch DNA, our data show that contacts it makes with the wHTH are key for specific DNA binding. The structure indicates how ssDNA-effector binding to the WYL-domain impacts wHTH DNA binding. In conclusion, we present the first structure of a WYL-activator bound to both effector and target DNA. The structure unveils a unique, asymmetric DNA binding mode that is likely conserved among WYL-activators.

C. elegans Dicer acts with the RIG-I-like helicase DRH-1 and RDE-4 to cleave dsRNA

Invertebrates use the endoribonuclease Dicer to cleave viral dsRNA during antiviral defense, while vertebrates use RIG-I-like Receptors (RLRs), which bind viral dsRNA to trigger an interferon response. While some invertebrate Dicers act alone during antiviral defense, C. elegans Dicer acts in a complex with a dsRNA binding protein called RDE-4, and an RLR ortholog called DRH-1. We used biochemical and structural techniques to provide mechanistic insight into how these proteins function together. We found RDE-4 is important for ATP-independent and ATP-dependent cleavage reactions, while helicase domains of both DCR-1 and DRH-1 contribute to ATP-dependent cleavage. DRH-1 plays the dominant role in ATP hydrolysis, and like mammalian RLRs, has an N-terminal domain that functions in autoinhibition. A cryo-EM structure indicates DRH-1 interacts with DCR-1's helicase domain, suggesting this interaction relieves autoinhibition. Our study unravels the mechanistic basis of the collaboration between two helicases from typically distinct innate immune defense pathways.

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.

Redox-Based Defect Detection in Packed DNA: Insights from Hybrid Quantum Mechanical/Molecular Mechanics Molecular Dynamics Simulations

The impact of an 8-oxoguanine (8oxoG) defect on the redox properties of DNA within the nucleosome core particle (NCP) was investigated employing hybrid quantum mechanical/molecular mechanics (QM/MM) molecular dynamics simulations of native and 8oxoG-containing NCP systems with an explicit representation of a biologically relevant environment. Two distinct NCP positions with varying solvent accessibility were considered for 8oxoG insertion. In both cases, it is found that the presence of 8oxoG drastically decreases the redox free energy of oxidation by roughly 1 eV, which is very similar to what was recently reported for free native and 8oxoG-containing DNA. In contrast, the effect of 8oxoG on the reorganization free energy is even smaller for packed DNA (decrease of 0.13 and 0.01 eV for defect-free and defect-containing systems, respectively) compared to the one for free DNA (0.25 eV), consistent with the increased rigidity of the NCP as compared to free DNA. Furthermore, the presence of an 8oxoG defect does not yield any significant changes in the packed DNA structure. Such a conclusion favors the idea that in the case of chromatin, defect-induced changes in DNA redox chemistry can also be exploited to detect damaged bases via DNA-mediated hole transfer.

One-Electron Oxidation Potentials and Hole Delocalization in Heterogeneous Single-Stranded DNA

The study of DNA processes is essential to understand not only its intrinsic biological functions but also its role in many innovative applications. The use of DNA as a nanowire or electrochemical biosensor leads to the need for a deep investigation of the charge transfer process along the strand as well as of the redox properties. In this contribution, the one-electron oxidation potential and the charge delocalization of the hole formed after oxidation are computationally investigated for different heterogeneous single-stranded DNA strands. We have established a two-step protocol: (i) molecular dynamics simulations in the frame of quantum mechanics/molecular mechanics (QM/MM) were performed to sample the conformational space; (ii) energetic properties were then obtained within a QM1/QM2/continuum approach in combination with the Marcus theory over an ensemble of selected geometries. The results reveal that the one-electron oxidation potential in the heterogeneous strands can be seen as a linear combination of that property within the homogeneous strands. In addition, the hole delocalization between different nucleobases is, in general, small, supporting the conclusion of a hopping mechanism for charge transport along the strands. However, charge delocalization becomes more important, and so does the tunneling mechanism contribution, when the reducing power of the nucleobases forming the strand is similar. Moreover, charge delocalization is slightly enhanced when there is a correlation between pairs of some of the interbase coordinates of the strand: twist/shift, twist/slide, shift/slide, and rise/tilt. However, the internal structure of the strand is not the predominant factor for hole delocalization but the specific sequence of nucleotides that compose the strand.

Probing the Nucleic Acid Flexibility to Disarm the Viral Counter-Defense Machinery: Design and Characterization of Potent p19 Inhibitors

Plant viruses are highly destructive and significant contributors to several global pandemics and epidemics in plants. A viral disease outbreak in plants can cause a scarcity of food supply and is a severe concern to humanity. The siRNA (small interfering RNA)-mediated RNA-induced silencing complex (RISC) formation is a primary defense mechanism in plants against viruses, where the RISC binds and degrades viral mRNAs. As a counter-defense, many viruses encode RNA-silencing suppressor proteins (e.g., the p19 protein from the Tombusviridae family) for viral proliferation in plants. The functional form of p19 (homodimer) binds to plant siRNA with high affinities, thereby interrupting the RISC formation and thus preventing the viral mRNA silencing in plants. By altering the RISC formation, the p19 protein helps the virus invasion in the plant and ultimately stunts host growth. In this study, we designed several modified siRNA-based molecules for p19 inhibition. The viral p19 protein is known to interact predominantly through H-bonds with 2'-OH and phosphates of the plant siRNA. We utilized this information and in silico-designed flexible substituents of siRNA, where we removed the C2'-C3' bond in each nucleotide unit. We performed all-atom explicit-solvent molecular dynamics simulations (400 ns, 3 replicates each) for control/modified siRNA─p19 complexes (8 in total) followed by energetic estimations. Strikingly, in a few modified complexes, the siRNA not only retained the double-helical structural integrity but also displayed remarkably enhanced p19 binding compared to the control siRNA; hence, we consider it important to perform biological and chemical in vitro and in vivo studies on proposed flexible nucleic acids as p19 inhibitors for crop protection.

Type II bacterial toxin-antitoxins: hypotheses, facts, and the newfound plethora of the PezAT system

Toxin-antitoxin (TA) systems are entities found in the prokaryotic genomes, with eight reported types. Type II, the best characterized, is comprised of two genes organized as an operon. Whereas toxins impair growth, the cognate antitoxin neutralizes its activity. TAs appeared to be involved in plasmid maintenance, persistence, virulence, and defence against bacteriophages. Most Type II toxins target the bacterial translational machinery. They seem to be antecessors of Higher Eukaryotes and Prokaryotes Nucleotide-binding (HEPN) RNases, minimal nucleotidyltransferase domains, or CRISPR-Cas systems. A total of four TAs encoded by Streptococcus pneumoniae, RelBE, YefMYoeB, Phd-Doc, and HicAB, belong to HEPN-RNases. The fifth is represented by PezAT/Epsilon-Zeta. PezT/Zeta toxins phosphorylate the peptidoglycan precursors, thereby blocking cell wall synthesis. We explore the body of knowledge (facts) and hypotheses procured for Type II TAs and analyse the data accumulated on the PezAT family. Bioinformatics analyses showed that homologues of PezT/Zeta toxin are abundantly distributed among 14 bacterial phyla mostly in Proteobacteria (48%), Firmicutes (27%), and Actinobacteria (18%), showing the widespread distribution of this TA. The pezAT locus was found to be mainly chromosomally encoded whereas its homologue, the tripartite omega-epsilon-zeta locus, was found mostly on plasmids. We found several orphan pezT/zeta toxins, unaccompanied by a cognate antitoxin.

Quantitative investigation of the effects of DNA modifications and protein mutations on MeCP2-MBD-DNA interactions

DNA methylation as an important epigenetic marker, has gained attention for the significance of three oxidative modifications (hydroxymethyl-C (hmC), formyl-C (fC), and carboxyl-C (caC)). Mutations occurring in the methyl-CpG-binding domain (MBD) of MeCP2 result in Rett. However, uncertainties persist regarding DNA modification and MBD mutation-induced interaction changes. Here, molecular dynamics simulations were used to investigate the underlying mechanisms behind changes due to different modifications of DNA and MBD mutations. Alanine scanning combined with the interaction entropy method was employed to accurately evaluate the binding free energy. The results show that MBD has the strongest binding ability for mCDNA, followed by caC, hmC, and fCDNA, with the weakest binding ability observed for CDNA. Further analysis revealed that mC modification induces DNA bending, causing residues R91 and R162 closer to the DNA. This proximity enhances van der Waals and electrostatic interactions. Conversely, the caC/hmC and fC modifications lead to two loop regions (near K112 and K130) closer to DNA, respectively. Furthermore, DNA modifications promote the formation of stable hydrogen bond networks, however mutations in the MBD significantly reduce the binding free energy. This study provides detailed insight into the effects of DNA modifications and MBD mutations on binding ability. It emphasizes the necessity for research and development of targeted Rett compounds that induce conformational compatibility between MBD and DNA, enhancing the stability and strength of their interactions.

Sequence-Dependent Properties of the RNA Duplex

Sequence-dependent properties of the DNA duplex have been accurately described using extensive molecular dynamics simulations. The RNA duplex meanwhile─which is typically represented as a sequence-averaged rigid rod─does not benefit from having equivalent molecular dynamics simulations. In this paper, we present a massive simulation effort using a set of ABC-optimized duplexes from which we derived tetramer-resolution properties of the RNA duplex and a simple mesoscopic model that can represent elastic properties of long RNA duplexes. Despite the extreme chemical similarity between DNA and RNA, the local and global elastic properties of the duplexes are very different. DNA duplexes show a complex and nonelastic pattern of flexibility, for instance, while RNA duplexes behave as an elastic system whose deformations can be represented by simple harmonic potentials. In RNA duplexes (RNA2), not only are intra- and interbase pair parameters (equilibrium and mechanical) different from those in the equivalent DNA duplex sequences (DNA2) but the correlations between movements also differ. Simple statements on the relative flexibility or stability of both polymers are meaningless and should be substituted by a more detailed description depending on the sequence and the type of deformation considered.

Thermodynamics and kinetics of DNA and RNA dinucleotide hybridization to gaps and overhangs

Hybridization of short nucleic acid segments (<4 nt) to single-strand templates occurs as a critical intermediate in processes such as nonenzymatic nucleic acid replication and toehold-mediated strand displacement. These templates often contain adjacent duplex segments that stabilize base pairing with single-strand gaps or overhangs, but the thermodynamics and kinetics of hybridization in such contexts are poorly understood because of the experimental challenges of probing weak binding and rapid structural dynamics. Here we develop an approach to directly measure the thermodynamics and kinetics of DNA and RNA dinucleotide dehybridization using steady-state and temperature-jump infrared spectroscopy. Our results suggest that dinucleotide binding is stabilized through coaxial stacking interactions with the adjacent duplex segments as well as from potential noncanonical base-pairing configurations and structural dynamics of gap and overhang templates revealed using molecular dynamics simulations. We measure timescales for dissociation ranging from 0.2-40 μs depending on the template and temperature. Dinucleotide hybridization and dehybridization involve a significant free energy barrier with characteristics resembling that of canonical oligonucleotides. Together, our work provides an initial step for predicting the stability and kinetics of hybridization between short nucleic acid segments and various templates.

Molecular dynamics of mismatch detection-How MutS uses indirect readout to find errors in DNA

The mismatch repair protein MutS safeguards genomic integrity by finding and initiating repair of basepairing errors in DNA. Single-molecule studies show MutS diffusing on DNA, presumably scanning for mispaired/unpaired bases, and crystal structures show a characteristic "mismatch-recognition" complex with DNA enclosed within MutS and kinked at the site of error. But how MutS goes from scanning thousands of Watson-Crick basepairs to recognizing rare mismatches remains unanswered, largely because atomic-resolution data on the search process are lacking. Here, 10 μs all-atom molecular dynamics simulations of Thermus aquaticus MutS bound to homoduplex DNA and T-bulge DNA illuminate the structural dynamics underlying the search mechanism. MutS-DNA interactions constitute a multistep mechanism to check DNA over two helical turns for its 1) shape, through contacts with the sugar-phosphate backbone, 2) conformational flexibility, through bending/unbending engineered by large-scale motions of the clamp domain, and 3) local deformability, through basepair destabilizing contacts. Thus, MutS can localize a potential target by indirect readout due to lower energetic costs of bending mismatched DNA and identify a site that distorts easily due to weaker base stacking and pairing as a mismatch. The MutS signature Phe-X-Glu motif can then lock in the mismatch-recognition complex to initiate repair.

Structural basis of sequence-specific cytosine deamination by double-stranded DNA deaminase toxin DddA

The interbacterial deaminase toxin DddA catalyzes cytosine-to-uracil conversion in double-stranded (ds) DNA and enables CRISPR-free mitochondrial base editing, but the molecular mechanisms underlying its unique substrate selectivity have remained elusive. Here, we report crystal structures of DddA bound to a dsDNA substrate containing the 5'-TC target motif. These structures show that DddA binds to the minor groove of a sharply bent dsDNA and engages the target cytosine extruded from the double helix. DddA Phe1375 intercalates in dsDNA and displaces the 5' (-1) thymine, which in turn replaces the target (0) cytosine and forms a noncanonical T-G base pair with the juxtaposed guanine. This tandem displacement mechanism allows DddA to locate a target cytosine without flipping it into the active site. Biochemical experiments demonstrate that DNA base mismatches enhance the DddA deaminase activity and relax its sequence selectivity. On the basis of the structural information, we further identified DddA mutants that exhibit attenuated activity or altered substrate preference. Our studies may help design new tools useful in genome editing or other applications.

Unveiling the interdomain dynamics of type II DNA topoisomerase through all-atom simulations: Implications for understanding its catalytic cycle

Type IIA DNA topoisomerases are complex molecular nanomachines that manage topological states of the DNA molecule in the cell and play a crucial role in cellular processes such as cell division and transcription. They are also established targets of cancer chemotherapy. Starting from the available crystal structure of a fully catalytic topoisomerase IIA homodimer from Saccharomyces cerevisiae, we constructed three states of this molecular motor primarily changing the configurations of the DNA segment bound in the DNA gate and performed μs-long all-atom molecular simulations. A comprehensive analysis revealed a sliding motion within the DNA gate and a teamwork between the N-gate and DNA gate that may be associated with the necessary molecular events that allow passage of the T-segment of DNA. The observed movement of the ATPase dimer relative to the DNA domain was reflected in different interaction patterns between the K-loops of the transducer domain and the B-A-B form of the bound DNA. Based on the obtained results, we mapped simulated configurations to the structures in the proposed catalytic cycle through which type IIA topoisomerases exert their function and discussed the possible transition events. The results extend our understanding of the mechanism of action of type IIA topoisomerases and provide an atomistic interpretation of some of the observed features of these molecular motors.

Structure, folding and flexibility of co-transcriptional RNA origami

RNA origami is a method for designing RNA nanostructures that can self-assemble through co-transcriptional folding with applications in nanomedicine and synthetic biology. However, to advance the method further, an improved understanding of RNA structural properties and folding principles is required. Here we use cryogenic electron microscopy to study RNA origami sheets and bundles at sub-nanometre resolution revealing structural parameters of kissing-loop and crossover motifs, which are used to improve designs. In RNA bundle designs, we discover a kinetic folding trap that forms during folding and is only released after 10 h. Exploration of the conformational landscape of several RNA designs reveal the flexibility of helices and structural motifs. Finally, sheets and bundles are combined to construct a multidomain satellite shape, which is characterized by individual-particle cryo-electron tomography to reveal the domain flexibility. Together, the study provides a structural basis for future improvements to the design cycle of genetically encoded RNA nanodevices.

Controlled sulfur-based engineering confers mouldability to phosphorothioate antisense oligonucleotides

Phosphorothioates (PS) have proven their effectiveness in the area of therapeutic oligonucleotides with applications spanning from cancer treatment to neurodegenerative disorders. Initially, PS substitution was introduced for the antisense oligonucleotides (PS ASOs) because it confers an increased nuclease resistance meanwhile ameliorates cellular uptake and in-vivo bioavailability. Thus, PS oligonucleotides have been elevated to a fundamental asset in the realm of gene silencing therapeutic methodologies. But, despite their wide use, little is known on the possibly different structural changes PS-substitutions may provoke in DNA·RNA hybrids. Additionally, scarce information and significant controversy exists on the role of phosphorothioate chirality in modulating PS properties. Here, through comprehensive computational investigations and experimental measurements, we shed light on the impact of PS chirality in DNA-based antisense oligonucleotides; how the different phosphorothioate diastereomers impact DNA topology, stability and flexibility to ultimately disclose pro-Sp S and pro-Rp S roles at the catalytic core of DNA Exonuclease and Human Ribonuclease H; two major obstacles in ASOs-based therapies. Altogether, our results provide full-atom and mechanistic insights on the structural aberrations PS-substitutions provoke and explain the origin of nuclease resistance PS-linkages confer to DNA·RNA hybrids; crucial information to improve current ASOs-based therapies.

Structural basis for water modulating RNA duplex formation in the CUG repeats of myotonic dystrophy type 1

Secondary structures formed by expanded CUG RNA are involved in the pathobiology of myotonic dystrophy type 1. Understanding the molecular basis of toxic RNA structures can provide insights into the mechanism of disease pathogenesis and accelerate the drug discovery process. Here, we report the crystal structure of CUG repeat RNA containing three U-U mismatches between C-G and G-C base pairs. The CUG RNA crystallizes as an A-form duplex, with the first and third U-U mismatches adopting a water-mediated asymmetric mirror isoform geometry. We found for the first time that a symmetric, water-bridged U-H₂O-U mismatch is well tolerated within the CUG RNA duplex, which was previously suspected but not observed. The new water-bridged U-U mismatch resulted in high base-pair opening and single-sided cross-strand stacking interactions, which in turn dominate the CUG RNA structure. Furthermore, we performed molecular dynamics (MD) simulations that complemented the structural findings and proposed that the first and third U-U mismatches are interchangeable conformations, while the central water-bridged U-U mismatch represents an intermediate state that modulates the RNA duplex conformation. Collectively, the new structural features provided in this work are important for understanding the recognition of U-U mismatches in CUG repeats by external ligands such as proteins or small molecules.

Synergistic binding of actinomycin D and echinomycin to DNA mismatch sites and their combined anti-tumour effects

Combination cancer chemotherapy is one of the most useful treatment methods to achieve a synergistic effect and reduce the toxicity of dosing with a single drug. Here, we use a combination of two well-established anticancer DNA intercalators, actinomycin D (ActD) and echinomycin (Echi), to screen their binding capabilities with DNA duplexes containing different mismatches embedded within Watson-Crick base-pairs. We have found that combining ActD and Echi preferentially stabilised thymine-related T:T mismatches. The enhanced stability of the DNA duplex-drug complexes is mainly due to the cooperative binding of the two drugs to the mismatch duplex, with many stacking interactions between the two different drug molecules. Since the repair of thymine-related mismatches is less efficient in mismatch repair (MMR)-deficient cancer cells, we have also demonstrated that the combination of ActD and Echi exhibits enhanced synergistic effects against MMR-deficient HCT116 cells and synergy is maintained in a MMR-related MLH1 gene knockdown in SW620 cells. We further accessed the clinical potential of the two-drug combination approach with a xenograft mouse model of a colorectal MMR-deficient cancer, which has resulted in a significant synergistic anti-tumour effect. The current study provides a novel approach for the development of combination chemotherapy for the treatment of cancers related to DNA-mismatches.

Computational Insights into the Dynamic Structural Features and Binding Characteristics of Recombinase UvsX Compared with RecA

RecA family recombinases are the core enzymes in the process of homologous recombination, and their normal operation ensures the stability of the genome and the healthy development of organisms. The UvsX protein from bacteriophage T4 is a member of the RecA family recombinases and plays a central role in T4 phage DNA repair and replication, which provides an important model for the biochemistry and genetics of DNA metabolism. UvsX shares a high degree of structural similarity and function with RecA, which is the most deeply studied member of the RecA family. However, the detailed molecular mechanism of UvsX has not been resolved. In this study, a comprehensive all-atom molecular dynamics simulation of the UvsX protein dimer complex was carried out in order to investigate the conformational and binding properties of UvsX in combination with ATP and DNA, and the simulation of RecA was synchronized with the property comparison learning for UvsX. This study confirmed the highly conserved molecular structure characteristics and catalytic centers of RecA and UvsX, and also discovered differences in regional conformation, volatility and the ability to bind DNA between the two proteins at different temperatures, which would be helpful for the subsequent understanding and application of related recombinases.

Direct monitoring of the thermodynamics and kinetics of DNA and RNA dinucleotide dehybridization from gaps and overhangs

Hybridization of short nucleic acid segments (<4 nucleotides) to single-strand templates occurs as a critical intermediate in processes such as non-enzymatic nucleic acid replication and toehold-mediated strand displacement. These templates often contain adjacent duplex segments that stabilize base pairing with single-strand gaps or overhangs, but the thermodynamics and kinetics of hybridization in such contexts are poorly understood due to experimental challenges of probing weak binding and rapid structural dynamics. Here we develop an approach to directly measure the thermodynamics and kinetics of DNA and RNA dinucleotide dehybridization using steady-state and temperature-jump infrared spectroscopy. Our results suggest that dinucleotide binding is stabilized through coaxial stacking interactions with the adjacent duplex segments as well as from potential non-canonical base pairing configurations and structural dynamics of gap and overhang templates revealed using molecular dynamics simulations. We measure timescales for dissociation ranging from 0.2 to 40 µs depending on the template and temperature. Dinucleotide hybridization and dehybridization involves a significant free energy barrier with characteristics resembling that of canonical oligonucleotides. Together, our work provides an initial step for predicting the stability and kinetics of hybridization between short nucleic acid segments and various templates.

Revealing intrinsic changes of DNA induced by spore photoproduct lesion through computer simulation

In bacterial endospores, a cross-linked thymine dimer, 5-thyminyl-5,6-dihydrothymine, commonly referred to as the spore photoproduct (SP), is found as the dominant DNA photo lesion under UV radiation. During spore germination, SP is faithfully repaired by the spore photoproduct lyase (SPL) for normal DNA replication to resume. Despite this general mechanism, the exact way in which SP modifies the duplex DNA structure so that the damaged site can be recognized by SPL to initiate the repair process is still unclear. A previous X-ray crystallographic study, which used a reverse transcriptase as a DNA host template, captured a protein-bound duplex oligonucleotide containing two SP lesions; the study showed shortened hydrogen bonds between the AT base pairs involved in the lesions and widened minor grooves near the damaged sites. However, it remains to be determined whether the results accurately reflect the conformation of SP-containing DNA (SP-DNA) in its fully hydrated pre-repair form. To uncover the intrinsic changes in DNA conformation caused by SP lesions, we performed molecular dynamics (MD) simulations of SP-DNA duplexes in aqueous solution, using the nucleic acid portion of the previously determined crystal structure as a template. After MD relaxation, our simulated SP-DNAs showed weakened hydrogen bonds at the damaged sites compared to those in the undamaged DNA. Our analyses of the MD trajectories revealed a range of local and global structural distortions of DNA induced by SP. Specifically, the SP region displays a greater tendency to adopt an A-like-DNA conformation, and curvature analysis revealed an increase in the global bending compared to the canonical B-DNA. Although these SP-induced DNA conformational changes are relatively minor, they may provide a sufficient structural basis for SP to be recognized by SPL during the lesion repair process.

An Attempt of Seeking Favorable Binding Free Energy Prediction Schemes Considering the Entropic Effect on Fis-DNA Binding

Protein-DNA binding mechanisms in a complex manner are essential for understanding many biological processes. Over the past decades, numerous experiments and calculations have analyzed the specificity of protein-DNA binding. However, the accuracy of binding free energy prediction for multi-base DNA systems still needs to be improved. Fis is a DNA-binding protein that regulates various transcription and recombination reactions. In the present work, we tested several methods of predict binding free energy based on this system to find a favorable prediction scheme and explore the binding mechanism of Fis protein and DNA. Two solvent models (explicit and implicit solvent models) were chosen for the dynamics process, and the predicted binding free energy was more accurate under the explicit solvent model. When different Poisson-Boltzmann/Generalized Born (PB/GB) models were tested for DNA force fields (BSC1 and OL15), it was found that the binding free energy predicted by the selected OL15 force field performed better and the correlation between predicted and experimental values was improved with the increasing interior dielectric constant (Dk). Finally, using Dk = 8, the GB^OBC1 model combined with interaction entropy (IE), which was calculated for entropic contribution (GB^OBC1_IE_8), was screened out for the binding free energy prediction and analysis of the Fis-DNA system, and the validity of the method was further verified by testing the Cren7-DNA system. By performing conformational analysis of the minor groove, it was found that mutation of the DNA central sequence A/T to C/G and deletion of the guanine 2-amino group would change the minor groove width and thus affect the formation of the major groove, altering the interaction and atomic contact between the protein and the major groove, thus changing the binding affinity of Fis and DNA. Hopefully, the series of tests in this work can shed some light on the related studies of protein and DNA systems.

Structures of NF-κB p52 homodimer-DNA complexes rationalize binding mechanisms and transcription activation

The mammalian NF-κB p52:p52 homodimer together with its cofactor Bcl3 activates transcription of κB sites with a central G/C base pair (bp), while it is inactive toward κB sites with a central A/T bp. To understand the molecular basis for this unique property of p52, we have determined the crystal structures of recombinant human p52 protein in complex with a P-selectin(PSel)-κB DNA (5'-GGGGTGACCCC-3') (central bp is underlined) and variants changing the central bp to A/T or swapping the flanking bp. The structures reveal a nearly two-fold widened minor groove in the central region of the DNA as compared to all other currently available NF-κB-DNA complex structures, which have a central A/T bp. Microsecond molecular dynamics (MD) simulations of free DNAs and p52 bound complexes reveal that free DNAs exhibit distinct preferred conformations, and p52:p52 homodimer induces the least amount of DNA conformational changes when bound to the more transcriptionally active natural G/C-centric PSel-κB, but adopts closed conformation when bound to the mutant A/T and swap DNAs due to their narrowed minor grooves. Our binding assays further demonstrate that the fast kinetics favored by entropy is correlated with higher transcriptional activity. Overall, our studies have revealed a novel conformation for κB DNA in complex with NF-κB and pinpoint the importance of binding kinetics, dictated by DNA conformational and dynamic states, in controlling transcriptional activation for NF-κB.

Electrical modulation properties of DNA drug molecules

DNA drug molecules are not only widely used in gene therapy, but also play an important role in controlling the electrical properties of molecular electronics. Covalent binding, groove binding and intercalation are all important forms of drug-DNA interaction. But its applications are limited due to a lack of understanding of the electron transport mechanisms after different drug-DNA interaction modes. Here, we used a combination of density functional theory calculations and nonequilibrium Green's function formulation with decoherence to study the effect of drug molecules on the charge transport property of DNA under three different binding modes. Conductance of DNA is found to decrease from 2.35E-5 G0 to 1.95E-6 G0 upon doxorubicin intercalation due to modifications of the density of states in the near-highest occupied molecular orbital region, δG = 1105.13%. Additionally, the conductance of DNA after cis-[Pt(NH3)2(py)Cl]+ covalent binding increases from 1.02E-6 G0 to 5.25E-5 G0, δG = 5047.06%. However, in the case of pentamidine groove binding, because there is no direct change in DNA molecular structure during drug binding, the conductance changes before and after drug binding is much smaller than in the two above cases, δG = 90.43%. Our theoretical calculations suggest that the conductance of DNA can be regulated by different drug molecules or switching the interaction modes between small molecules and DNA. This regulation opens new possibilities for their potential applications in controllable modulation of the electron transport property of DNA.

DNA-TCP complex structures reveal a unique recognition mechanism for TCP transcription factor families

Plant-specific TCP transcription factors are key regulators of diverse plant functions. TCP transcription factors have long been annotated as basic helix-loop-helix (bHLH) transcription factors according to remote sequence homology without experimental validation, and their consensus DNA-binding sequences and protein-DNA recognition mechanisms have remained elusive. Here, we report the crystal structures of the class I TCP domain from AtTCP15 and the class II TCP domain from AtTCP10 in complex with different double-stranded DNA (dsDNA). The complex structures reveal that the TCP domain is a distinct DNA-binding motif and the homodimeric TCP domains adopt a unique three-site recognition mode, binding to dsDNA mainly through a central pair of β-strands formed by the dimer interface and two basic flexible loops from each monomer. The consensus DNA-binding sequence for class I TCPs is a perfectly palindromic 11 bp (GTGGGNCCCAC), whereas that for class II TCPs is a near-palindromic 11 bp (GTGGTCCCCAC). The unique DNA binding mode allows the TCP domains to display broad specificity for a range of DNA sequences even shorter than 11 bp, adding further complexity to the regulatory network of plant TCP transcription factors.

Harmonizing Interstrand Electrostatic Repulsion by Conformational Rigidity in Counterion-Deprived Z-DNA: A Molecular Dynamics Study

Deoxyribonucleic acid (DNA) is a vital biomacromolecule. Although the right-handed B-DNA type helical structure is the most abundant and extensively studied form of DNA, several noncanonical forms, such as triplex, quadruplex, Z-DNA, A-DNA, and ss-DNA, have been probed from time to time to gain insights into the DNA's function. Z-DNA was recently found to be involved in cancer and several autoimmune diseases. In the present Article, we evaluated the conformational stability of locked-sugar-based Z-DNA via all-atom explicit-solvent molecular dynamics simulations and found that the modified DNA maintained the left-handed conformation even in the absence of counterions, wherein the structural rigidity dominates over the electrostatic repulsion between the complementary strands. The control Z-DNA without counterions, as expected, instantaneously resulted in unfolded states. The remarkable stability of the conformationally locked model system was thoroughly investigated via structural and energetic perspectives and was probably the result of the backbone widening in tandem with enhanced electrostatics between complementary strands. We believe that the design of the proposed modified Z-DNA construct could help understand the otherwise delicate Z-DNA conformation even in salt-deprived conditions. The design could also motivate the medicinal use of short segments of such modified nucleotides and could be utilized in more advanced modeling techniques, such as DNA origami which has gained popularity in recent years.

Molecular basis for DNA recognition by the maternal pioneer transcription factor FoxH1

Forkhead box H1 (FoxH1) is an essential maternal pioneer factor during embryonic development that binds to specific GG/GT-containing DNA target sequences. Here we have determined high-resolution structures of three FoxH1 proteins (from human, frog and fish species) and four DNAs to clarify the way in which FoxH1 binds to these sites. We found that the protein-DNA interactions extend to both the minor and major DNA grooves and are thus almost twice as extensive as those of other FOX family members. Moreover, we identified two specific amino acid changes in FoxH1 that allowed the recognition of GG/GT motifs. Consistent with the pioneer factor activity of FoxH1, we found that its affinity for nucleosomal DNA is even higher than for linear DNA fragments. The structures reported herein illustrate how FoxH1 binding to distinct DNA sites provides specificity and avoids cross-regulation by other FOX proteins that also operate during the maternal-zygotic transition and select canonical forkhead sites.

DNA Sequence-Dependent Binding of Linker Histone gH1 Regulates Nucleosome Conformations

Sequence-dependent binding between DNA and proteins in chromatin is an essential part of gene expression. Linker histone H1 is an important protein in the regulation of chromatin compartmentalization and compaction, and its binding with the nucleosome is sensitive to the DNA sequence. Although the interactions of H1 and DNA have been widely investigated, the mechanism of nucleosome conformation changes induced by the DNA-sequence-dependent binding with gH1 (globular H1.0) remains largely unclear at the atomic level. In the present molecular dynamics simulations, both linker and dyad DNAs were mutated to investigate the conformational changes of the nucleosome induced by the sequence-dependent binding of gH1 based on the on-dyad binding mode. Our results indicate that gH1 is insensitive to the DNA sequence of the dyad DNA but presents an apparent preference to linker DNA with an AT-rich sequence. Moreover, this specific binding induces the entry/exit region of a nucleosome to a tight conformation and regulates the accessibility of core histones. Considering that the entry/exit region of the nucleosome is a crucial binding site for many functional proteins related to gene expression, the conformational change at this region could represent an important gene regulation signal.

RNABPDB: Molecular Modeling of RNA Structure-From Base Pair Analysis in Crystals to Structure Prediction

The stable three-dimensional structure of RNA is known to play several important biochemical roles, from post-transcriptional gene regulation to enzymatic action. These structures contain double-helical regions, which often have different types of non-canonical base pairs in addition to Watson-Crick base pairs. Hence, it is important to study their structures from experimentally obtained or even predicted ones, to understand their role, or to develop a drug against the potential targets. Molecular Modeling of RNA double helices containing non-canonical base pairs is a difficult process, particularly due to the unavailability of structural features of non-Watson-Crick base pairs. Here we show a composite web-server with an associated database that allows one to generate the structure of RNA double helix containing non-canonical base pairs using consensus parameters obtained from the database. The database classification is followed by an evaluation of the central tendency of the structural parameters as well as a quantitative estimation of interaction strengths. These parameters are used to construct three-dimensional structures of double helices composed of Watson-Crick and/or non-canonical base pairs. Our benchmark study to regenerate double-helical fragments of many experimentally derived RNA structures indicate very high accuracy. This composite server is expected to be highly useful in understanding functions of various pre-miRNA by modeling structures of the molecules and estimating binding efficiency. The database can be accessed from http://hdrnas.saha.ac.in/rnabpdb .

Drug design and DNA structural research inspired by the Neidle laboratory: DNA minor groove binding and transcription factor inhibition by thiophene diamidines

The understanding of sequence-specific DNA minor groove interactions has recently made major steps forward and as a result, the goal of development of compounds that target the minor groove is an active research area. In an effort to develop biologically active minor groove agents, we are preparing and exploring the DNA interactions of diverse diamidine derivatives with a 5'-GAATTC-3' binding site using a powerful array of methods including, biosensor-SPR methods, and X-ray crystallography. The benzimidazole-thiophene module provides an excellent minor groove recognition component. A central thiophene in a benzimidazole-thiophene-phenyl aromatic system provides essentially optimum curvature for matching the shape of the minor groove. Comparison of that structure to one with the benzimidazole replaced with an indole shows that the two structures are very similar, but have some interesting and important differences in electrostatic potential maps, the DNA minor groove binding structure based on x-ray crystallographic analysis, and inhibition of the major groove binding PU.1 transcription factor complex. The binding K_D for both compounds is under 10 nM and both form amidine H-bonds to DNA bases. They both have bifurcated H-bonds from the benzimidazole or indole groups to bases at the center of the -AATT- binding site. Analysis of the comparative results provides an excellent understanding of how thiophene compounds recognize the minor groove and can act as transcription factor inhibitors.

Staggered intercalation of DNA duplexes with base-pair modulation by two distinct drug molecules induces asymmetric backbone twisting and structure polymorphism

The use of multiple drugs simultaneously targeting DNA is a promising strategy in cancer therapy for potentially overcoming single drug resistance. In support of this concept, we report that a combination of actinomycin D (ActD) and echinomycin (Echi), can interact in novel ways with native and mismatched DNA sequences, distinct from the structural effects produced by either drug alone. Changes in the former with GpC and CpG steps separated by a A:G or G:A mismatch or in a native DNA with canonical G:C and C:G base pairs, result in significant asymmetric backbone twists through staggered intercalation and base pair modulations. A wobble or Watson–Crick base pair at the two drug-binding interfaces can result in a single-stranded ‘chair-shaped’ DNA duplex with a straight helical axis. However, a novel sugar-edged hydrogen bonding geometry in the G:A mismatch leads to a ‘curved-shaped’ duplex. Two non-canonical G:C Hoogsteen base pairings produce a sharply kinked duplex in different forms and a four-way junction-like superstructure, respectively. Therefore, single base pair modulations on the two drug-binding interfaces could significantly affect global DNA structure. These structures thus provide a rationale for atypical DNA recognition via multiple DNA intercalators and a structural basis for the drugs’ potential synergetic use.

DNA Deformation Exerted by Regulatory DNA-Binding Motifs in Human Alkyladenine DNA Glycosylase Promotes Base Flipping

Human alkyladenine DNA glycosylase (AAG) is a key enzyme that corrects a broad range of alkylated and deaminated nucleobases to maintain genomic integrity. When encountering the lesions, AAG adopts a base-flipping strategy to extrude the target base from the DNA duplex to its active site, thereby cleaving the glycosidic bond. Despite its functional importance, the detailed mechanism of such base extrusion and how AAG distinguishes the lesions from an excess of normal bases both remain elusive. Here, through the Markov state model constructed on extensive all-atom molecular dynamics simulations, we find that the alkylated nucleobase (N3-methyladenine, 3MeA) everts through the DNA major groove. Two key AAG motifs, the intercalation and E131-N146 motifs, play active roles in bending/pressing the DNA backbone and widening the DNA minor groove during 3MeA eversion. In particular, the intercalated residue Y162 is involved in buckling the target site at the early stage of 3MeA eversion. Our traveling-salesman based automated path searching algorithm further revealed that a non-target normal adenine tends to be trapped in an exo site near the active site, which however barely exists for a target base 3MeA. Collectively, these results suggest that the Markov state model combined with traveling-salesman based automated path searching acts as a promising approach for studying complex conformational changes of biomolecules and dissecting the elaborate mechanism of target recognition by this unique enzyme.

Structural basis of transposon end recognition explains central features of Tn7 transposition systems

Tn7 is a bacterial transposon with relatives containing element-encoded CRISPR-Cas systems mediating RNA-guided transposon insertion. Here, we present the 2.7 Å cryoelectron microscopy structure of prototypic Tn7 transposase TnsB interacting with the transposon end DNA. When TnsB interacts across repeating binding sites, it adopts a beads-on-a-string architecture, where the DNA-binding and catalytic domains are arranged in a tiled and intertwined fashion. The DNA-binding domains form few base-specific contacts leading to a binding preference that requires multiple weakly conserved sites at the appropriate spacing to achieve DNA sequence specificity. TnsB binding imparts differences in the global structure of the protein-bound DNA ends dictated by the spacing or overlap of binding sites explaining functional differences in the left and right ends of the element. We propose a model of the strand-transfer complex in which the terminal TnsB molecule is rearranged so that its catalytic domain is in a position conducive to transposition.

Mechanism of RNA polymerase I selection by transcription factor UAF

Preribosomal RNA is selectively transcribed by RNA polymerase (Pol) I in eukaryotes. The yeast transcription factor upstream activating factor (UAF) represses Pol II transcription and mediates Pol I preinitiation complex (PIC) formation at the 35S ribosomal RNA gene. To visualize the molecular intermediates toward PIC formation, we determined the structure of UAF in complex with native promoter DNA and transcription factor TATA-box-binding protein (TBP). We found that UAF recognizes DNA using a hexameric histone-like scaffold with markedly different interactions compared with the nucleosome and the histone-fold-rich transcription factor IID (TFIID). In parallel, UAF positions TBP for Core Factor binding, which leads to Pol I recruitment, while sequestering it from DNA and Pol II/III-specific transcription factors. Our work thus reveals the structural basis of RNA Pol selection by a transcription factor.

Modulation of DNA conformation in electrolytic nanodroplets

The behavior of deoxyribonucleic acid (DNA) molecules in confinement is of profound importance in various bioengineering and medical applications. In the present study, all-atom molecular dynamics simulation is utilized to investigate the transition of the double-strand DNA (dsDNA) conformation in the electrolytic nanodroplet. Three typical conformations, i.e., C-shaped, folded S-shaped, and double C-shaped, are observed for different droplet sizes and ionic concentrations. To reveal the physics underlying this phenomenon, the characteristics of the dsDNA molecules, such as the overcharging intensity, the end-to-end distance, the radius of gyration, etc. are analyzed in detail based on the numerical results. It is found that the transition can be ascribed to the buckling of the polymer molecules under the compression due to the confinement of the nanodroplet, and it can be modulated by the ionic concentration in the electrolyte. Generally, nanoscale confinement dominates dsDNA behavior over the electrostatic effects in smaller nanodroplets, while the latter becomes more important for larger nanodroplets. This competition results in the persistence length increasing with the nanodroplet radii. Based on these discussions, a non-dimensional elasto-capillary number μ is proposed to classify the dsDNA conformations into three regions.

Deciphering the mechanism of p73 recognition of p53 response elements using the crystal structure of p73-DNA complexes and computational studies

P73 belongs to p53 family transcription factor activating more than 50% of cell fate p53 target genes involved in cell cycle, apoptosis, DNA damage response alongside neuronal system development and differentiation by binding to 20-bp response elements (REs) having sequence motif (PPPCA/T-T/AGYYY) where P-purines and Y-pyrimidines with each 10-bp separated by minimum 0 to 13-bp spacer. The promiscuous nature of recognizing both cell fate and development genes and the underlying RE selectivity mechanism by p73 is not well understood. Here, we report the molecular details of p73 recognizing the REs using the crystal structure of p73 DNA binding domain (DBD) in complex with 12 base pair DNA sequence 5'-cAGGCATGCCTg-3' and molecular dynamics simulations with six different p53 natural promoter sequences. Each 20-base pair natural promoter forms a different major/minor groove due to the presence of nucleotides A/T, A/C, G/G, T/T and G/T at positions 3, 8, 13, 18 uniquely recognized by p73 key residues Lys138 and Arg268. The loops L1 and L3 bearing these residues influence inter-and intra-dimer interfaces interactions and hence p73 forms a unique tetramer with each natural promoter sequence. Structural features of the DNA and the spacing between half-sites influence p73 tetramerization and its transactivation function.

Structure of the 5' untranslated region in SARS-CoV-2 genome and its specific recognition by innate immune system via the human oligoadenylate synthase 1

2'-5'-Oligoadenylate synthetase 1 (OAS1) is one of the key enzymes driving the innate immune system response to SARS-CoV-2 infection whose activity has been related to COVID-19 severity. OAS1 is a sensor of endogenous RNA that triggers the 2'-5'-oligoadenylate/RNase L pathway. Upon SARS-CoV-2 infection, OAS1 is responsible for the recognition of viral RNA and has been shown to possess a particularly high sensitivity for the 5'-untranslated (5'-UTR) RNA region, which is organized in a double-strand stem loop motif (SL1). Here we report the structure of the SL1/OAS1 complex also rationalizing the high affinity for OAS1.

Oxidation Enhances Binding of Extrahelical 5-Methyl-Cytosines by Thymine DNA Glycosylase

The DNA repair protein thymine DNA glycosylase (TDG) removes mispaired or damaged bases, such as oxidized methyl-cytosine, from DNA by cleavage of the glycosidic bond between the sugar and the target base flipped into the enzyme's active site. The enzyme is active against formyl-cytosine and carboxyl-cytosine, whereas the lower oxidized hydroxymethyl-cytosine and methyl-cytosine itself are not processed by the enzyme. Molecular dynamics simulations with thermodynamic integration of TDG complexed to DNA carrying one of four different (oxidized) methyl-cytosine bases in extrahelcial conformation, methyl-cytosine (mC), hydroxymethyl-cytosine (hmC), formyl-cytosine (fC), or carboxyl-cytosine (caC), show a more favorable binding affinity of the higher oxidized forms, fC and caC, than the nonsubstrate bases hmC and mC. Despite rather comparable, reaction-competent conformations of the flipped bases in the active site of the enzyme, more and stronger interactions with active site residues account for the preferred binding of the higher oxidized bases. Binding of the negatively charged caC and the neutral fC are strengthened by interactions with positively charged His151. Our calculated proton affinities find this protonation state of His151 the preferred one in the presence of caC and conceivable in the presence of fC as well as increasing the binding affinity toward the two bases. Discrimination of the substrate bases is further achieved by the backbone of Tyr152 that forms a strong hydrogen bond to the carboxyl and formyl oxygen atoms of caC and fC, respectively, a contact that is completely lacking in mC and much weaker in hmC. Overall, our computational results indicate that the enzyme discriminates the different oxidation forms of methyl-cytosine already at the formation of the extrahelical complexes.

Role of helical structure and dynamics in oligoadenylate synthetase 1 (OAS1) mismatch tolerance and activation by short dsRNAs

The 2'-5'-oligoadenylate synthetases (OAS) are innate immune sensors of cytosolic double-stranded RNA (dsRNA) that play a critical role in limiting viral infection. How these proteins are able to avoid aberrant activation by cellular RNAs is not fully understood, but adenosine-to-inosine (A-to-I) editing has been proposed to limit accumulation of endogenous RNAs that might otherwise cause stimulation of the OAS/RNase L pathway. Here, we aim to uncover whether and how such sequence modifications can restrict the ability of short, defined dsRNAs to activate the single-domain form of OAS, OAS1. Unexpectedly, we find that all tested inosine-containing dsRNAs have an increased capacity to activate OAS1, whether in a destabilizing (I•U) or standard Watson-Crick-like base pairing (I-C) context. Additional variants with strongly destabilizing A•C mismatches or stabilizing G-C pairs also exhibit increased capacity to activate OAS1, eliminating helical stability as a factor in the relative ability of the dsRNAs to activate OAS1. Using thermal difference spectra and molecular dynamics simulations, we identify both increased helical dynamics and specific local changes in helical structure as important factors in the capacity of short dsRNAs to activate OAS1. These helical features may facilitate more ready adoption of the distorted OAS1-bound conformation or stabilize important structures to predispose the dsRNA for optimal binding and activation of OAS1. These studies thus reveal the molecular basis for the greater capacity of some short dsRNAs to activate OAS1 in a sequence-independent manner.

Molecular basis of Arginine and Lysine DNA sequence-dependent thermo-stability modulation

We have used a variety of theoretical and experimental techniques to study the role of four basic amino acids-Arginine, Lysine, Ornithine and L-2,4-Diaminobutyric acid-on the structure, flexibility and sequence-dependent stability of DNA. We found that the presence of organic ions stabilizes the duplexes and significantly reduces the difference in stability between AT- and GC-rich duplexes with respect to the control conditions. This suggests that these amino acids, ingredients of the primordial soup during abiogenesis, could have helped to equalize the stability of AT- and GC-rich DNA oligomers, facilitating a general non-catalysed self-replication of DNA. Experiments and simulations demonstrate that organic ions have an effect that goes beyond the general electrostatic screening, involving specific interactions along the grooves of the double helix. We conclude that organic ions, largely ignored in the DNA world, should be reconsidered as crucial structural elements far from mimics of small inorganic cations.