3DNA: a versatile, integrated software system for the analysis, rebuilding and visualization of three-dimensional nucleic-acid structures

Localization and discrimination of GG mismatch in duplex DNA by synthetic ligand-enhanced protein nanopore analysis

Mismatched base pairs in DNA are the basis of single-nucleotide polymorphism, one of the major issues in genetic diseases. However, the changes of physical and chemical properties of DNA caused by single-site mismatches are often influenced by the sequence and the structural flexibility of the whole duplex, resulting in a challenge of direct detection of the types and location of mismatches sensitively. In this work, we proposed a synthetic ligand-enhanced protein nanopore analysis of GG mismatch on DNA fragment, inspired by in silico investigation of the specific binding of naphthyridine dimer (ND) on GG mismatch. We demonstrated that both the importing and unzipping processes of the ligand-bound DNA duplex can be efficiently slowed down in α-hemolysin nanopore. This ligand-binding induced slow-down effect of DNA in nanopore is also sensitive to the relative location of the mismatch on DNA duplex. Especially, the GG mismatch close to the end of a DNA fragment, which is hard to be detected by either routine nanopore analysis or tedious nanopore sequencing, can be well differentiated by our ND-enhanced nanopore experiment. These findings provide a promising strategy to localize and discriminate base mismatches in duplex form directly at the single-molecule level.

RNAtango: Analysing and comparing RNA 3D structures via torsional angles

RNA molecules, essential for viruses and living organisms, derive their pivotal functions from intricate 3D structures. To understand these structures, one can analyze torsion and pseudo-torsion angles, which describe rotations around bonds, whether real or virtual, thus capturing the RNA conformational flexibility. Such an analysis has been made possible by RNAtango, a web server introduced in this paper, that provides a trigonometric perspective on RNA 3D structures, giving insights into the variability of examined models and their alignment with reference targets. RNAtango offers comprehensive tools for calculating torsion and pseudo-torsion angles, generating angle statistics, comparing RNA structures based on backbone torsions, and assessing local and global structural similarities using trigonometric functions and angle measures. The system operates in three scenarios: single model analysis, model-versus-target comparison, and model-versus-model comparison, with results output in text and graphical formats. Compatible with all modern web browsers, RNAtango is accessible freely along with the source code. It supports researchers in accurately assessing structural similarities, which contributes to the precision and efficiency of RNA modeling.

Nucleoprotein phase-separation affinities revealed via atomistic simulations of short peptide and RNA fragments

Liquid-liquid phase separation of proteins and nucleic acids into condensate phases is a versatile mechanism for ensuring compartmentalization of cellular biochemistry. RNA molecules play critical roles in these condensates, particularly in transcriptional regulation and stress responses, exhibiting a wide range of thermodynamic and dynamic behaviors. However, deciphering the molecular grammar that governs the stability and dynamics of protein-RNA condensates remains challenging due to the multicomponent and heterogeneous nature of these biomolecular mixtures. In this study, we employ atomistic simulations of twenty distinct mixtures containing minimal RNA and peptide fragments to dissect the phase-separating affinities of all twenty amino acids in the presence of RNA. Our findings elucidate chemically specific interactions, hydration profiles, and ionic effects that synergistically promote or suppress protein-RNA phase separation. We map a ternary phase diagram of interactions, identifying four distinct groups of residues that promote, maintain, suppress, or disrupt protein-RNA clusters.

Relative genotoxicity of polycyclic aromatic hydrocarbons inferred from free energy perturbation approaches

Utilizing molecular dynamics and free energy perturbation, we examine the relative binding affinity of several covalent polycyclic aromatic hydrocarbon - DNA (PAH-DNA) adducts at the central adenine of NRAS codon-61, a mutational hotspot implicated in cancer risk. Several PAHs classified by the International Agency for Research on Cancer as probable, possible, or unclassifiable as to carcinogenicity are found to have greater binding affinity than the known carcinogen, benzo[a]pyrene (B[a]P). van der Waals interactions between the intercalated PAH and neighboring nucleobases, and minimal disruption of the DNA duplex drive increases in binding affinity. PAH-DNA adducts may be repaired by global genomic nucleotide excision repair (GG-NER), hence we also compute relative free energies of complexation of PAH-DNA adducts with RAD4-RAD23 (the yeast ortholog of human XPC-RAD23) which constitutes the recognition step in GG-NER. PAH-DNA adducts exhibiting the greatest DNA binding affinity also exhibit the least RAD4-RAD23 complexation affinity and are thus predicted to resist the GG-NER machinery, contributing to their genotoxic potential. In particular, the fjord region PAHs dibenzo[a,l]pyrene, benzo[g]chrysene, and benzo[c]phenanthrene are found to have greater binding affinity while having weaker RAD4-RAD23 complexation affinity than their respective bay region analogs B[a]P, chrysene, and phenanthrene. We also find that the bay region PAHs dibenzo[a,j]anthracene, dibenzo[a,c]anthracene, and dibenzo[a,h]anthracene exhibit greater binding affinity and weaker RAD4-RAD23 complexation affinity than B[a]P. Thus, the study of PAH genotoxicity likely needs to be substantially broadened, with implications for public policy and the health sciences. This approach can be broadly applied to assess factors contributing to the genotoxicity of other unclassified compounds.

From straight to curved: A historical perspective of DNA shape

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

Utilizing biological experimental data and molecular dynamics for the classification of mutational hotspots through machine learning

Motivation

Benzo[a]pyrene, a notorious DNA-damaging carcinogen, belongs to the family of polycyclic aromatic hydrocarbons commonly found in tobacco smoke. Surprisingly, nucleotide excision repair (NER) machinery exhibits inefficiency in recognizing specific bulky DNA adducts including Benzo[a]pyrene Diol-Epoxide (BPDE), a Benzo[a]pyrene metabolite. While sequence context is emerging as the leading factor linking the inadequate NER response to BPDE adducts, the precise structural attributes governing these disparities remain inadequately understood. We therefore combined the domains of molecular dynamics and machine learning to conduct a comprehensive assessment of helical distortion caused by BPDE-Guanine adducts in multiple gene contexts. Specifically, we implemented a dual approach involving a random forest classification-based analysis and subsequent feature selection to identify precise topological features that may distinguish adduct sites of variable repair capacity. Our models were trained using helical data extracted from duplexes representing both BPDE hotspot and nonhotspot sites within the TP53 gene, then applied to sites within TP53, cII, and lacZ genes.

Results

We show our optimized model consistently achieved exceptional performance, with accuracy, precision, and f1 scores exceeding 91%. Our feature selection approach uncovered that discernible variance in regional base pair rotation played a pivotal role in informing the decisions of our model. Notably, these disparities were highly conserved among TP53 and lacZ duplexes and appeared to be influenced by the regional GC content. As such, our findings suggest that there are indeed conserved topological features distinguishing hotspots and nonhotpot sites, highlighting regional GC content as a potential biomarker for mutation.

Availability and implementation

Code for comparing machine learning classifiers and evaluating their performance is available at https://github.com/jdavies24/ML-Classifier-Comparison, and code for analysing DNA structure with Curves+ and Canal using Random Forest is available at https://github.com/jdavies24/ML-classification-of-DNA-trajectories.

Quantum mechanical electronic and geometric parameters for DNA k-mers as features for machine learning

We are witnessing a steep increase in model development initiatives in genomics that employ high-end machine learning methodologies. Of particular interest are models that predict certain genomic characteristics based solely on DNA sequence. These models, however, treat the DNA as a mere collection of four, A, T, G and C, letters, dismissing the past advancements in science that can enable the use of more intricate information from nucleic acid sequences. Here, we provide a comprehensive database of quantum mechanical (QM) and geometric features for all the permutations of 7-meric DNA in their representative B, A and Z conformations. The database is generated by employing the applicable high-cost and time-consuming QM methodologies. This can thus make it seamless to associate a wealth of novel molecular features to any DNA sequence, by scanning it with a matching k-meric window and pulling the pre-computed values from our database for further use in modelling. We demonstrate the usefulness of our deposited features through their exclusive use in developing a model for A->C mutation rates.

NMR structures and magnetic force spectroscopy studies of small molecules binding to models of an RNA CAG repeat expansion

RNA repeat expansions fold into stable structures and cause microsatellite diseases such as Huntington's disease (HD), myotonic dystrophy type 1 (DM1), and spinocerebellar ataxias (SCAs). The trinucleotide expansion of r(CAG), or r(CAG)exp, causes both HD and SCA3, and the RNA's toxicity has been traced to its translation into polyglutamine (polyQ; HD) as well as aberrant pre-mRNA alternative splicing (SCA3 and HD). Previously, a small molecule, 1, was discovered that binds to r(CAG)exp and rescues aberrant pre-mRNA splicing in patient-derived fibroblasts by freeing proteins bound to the repeats. Here, we report the structures of single r(CAG) repeat motif (5'CAG/3'GAC where the underlined adenosines form a 1×1 nucleotide internal loop) in complex with 1 and two other small molecules via nuclear magnetic resonance (NMR) spectroscopy combined with simulated annealing. Compound 2 was designed based on the structure of 1 bound to the RNA while 3 was selected as a diverse chemical scaffold. The three complexes, although adopting different 3D binding pockets upon ligand binding, are stabilized by a combination of stacking interactions with the internal loop's closing GC base pairs, hydrogen bonds, and van der Waals interactions. Molecular dynamics (MD) simulations performed with NMR-derived restraints show that the RNA is stretched and bent upon ligand binding with significant changes in propeller-twist and opening. Compound 3 has a distinct mode of binding by insertion into the helix, displacing one of the loop nucleotides into the major groove and affording a rod-like shape binding pocket. In contrast, 1 and 2 are groove binders. A series of single molecule magnetic force spectroscopy studies provide a mechanistic explanation for how bioactive compounds might rescue disease-associated cellular phenotypes.

CRISPR-Cas12a bends DNA to destabilize base pairs during target interrogation

RNA-guided endonucleases are involved in processes ranging from adaptive immunity to site-specific transposition and have revolutionized genome editing. CRISPR-Cas9, -Cas12 and related proteins use guide RNAs to recognize ~20-nucleotide target sites within genomic DNA by mechanisms that are not yet fully understood. We used structural and biochemical methods to assess early steps in DNA recognition by Cas12a protein-guide RNA complexes. We show here that Cas12a initiates DNA target recognition by bending DNA to induce transient nucleotide flipping that exposes nucleobases for DNA-RNA hybridization. Cryo-EM structural analysis of a trapped Cas12a-RNA-DNA surveillance complex and fluorescence-based conformational probing show that Cas12a-induced DNA helix destabilization enables target discovery and engagement. This mechanism of initial DNA interrogation resembles that of CRISPR-Cas9 despite distinct evolutionary origins and different RNA-DNA hybridization directionality of these enzyme families. Our findings support a model in which RNA-mediated DNA engineering begins with local helix distortion by transient CRISPR-Cas protein binding.

Formation of left-handed helices by C2'-fluorinated nucleic acids under physiological salt conditions

Recent findings in cell biology have rekindled interest in Z-DNA, the left-handed helical form of DNA. We report here that two minimally modified nucleosides, 2'F-araC and 2'F-riboG, induce the formation of the Z-form under low ionic strength. We show that oligomers entirely made of these two nucleosides exclusively produce left-handed duplexes that bind to the Zα domain of ADAR1. The effect of the two nucleotides is so dramatic that Z-form duplexes are the only species observed in 10 mM sodium phosphate buffer and neutral pH, and no B-form is observed at any temperature. Hence, in contrast to other studies reporting formation of Z/B-form equilibria by a preference for purine glycosidic angles in syn, our NMR and computational work revealed that sequential 2'F…H2N and intramolecular 3'H…N3' interactions stabilize the left-handed helix. The equilibrium between B- and Z- forms is slow in the 19F NMR time scale (≥ms), and each conformation exhibited unprecedented chemical shift differences in the 19F signals. This observation led to a reliable estimation of the relative population of B and Z species and enabled us to monitor B-Z transitions under different conditions. The unique features of 2'F-modified DNA should thus be a valuable addition to existing techniques for specific detection of new Z-binding proteins and ligands.

High-resolution RNA tertiary structures in Zika virus stem-loop A for the development of inhibitory small molecules

Flaviviruses such as Zika (ZIKV) and dengue virus (DENV) are positive-sense RNA viruses belonging to Flaviviridae The flavivirus genome contains a 5' end stem-loop promoter sequence known as stem-loop A (SLA) that is recognized by the flavivirus polymerase NS5 during viral RNA synthesis and 5' guanosine cap methylation. The crystal structures of ZIKV and DENV SLAs show a well-defined fold, consisting of a bottom stem, side loop, and top stem-loop, providing unique interaction sites for small molecule inhibitors to disrupt the promoter function. To facilitate the identification of small molecule binding sites in flavivirus SLA, we determined high-resolution structures of the bottom and top stems of ZIKV SLA, which contain a single U- or G-bulge, respectively. Both bulge nucleotides exhibit multiple orientations, from folded back on the adjacent nucleotide to flipped out of the helix, and are stabilized by stacking or base triple interactions. These structures suggest that even a single unpaired nucleotide can provide flexibility to RNA structures, and its conformation is mainly determined by the stabilizing chemical environment. To facilitate discovery of small molecule inhibitors that interfere with the functions of ZIKV SLA, we screened and identified compounds that bind to the bottom and top stems of ZIKV SLA.

Quantitative and systematic NMR measurements of sequence-dependent A-T Hoogsteen dynamics uncovers unique conformational specificity in the DNA double helix

The propensities to form lowly-populated short-lived conformations of DNA could vary with sequence, providing an important source of sequence-specificity in biochemical reactions. However, comprehensively measuring how these dynamics vary with sequence is challenging. Using 1H CEST and 13C R 1 ρ NMR, we measured Watson-Crick to Hoogsteen dynamics for an A-T base pair in thirteen trinucleotide sequence contexts. The Hoogsteen population and exchange rate varied 4-fold and 16-fold, respectively, and were dependent on both the 3'- and 5'-neighbors but only weakly dependent on monovalent ion concentration (25 versus 100 mM NaCl) and pH (6.8 versus 8.0). Flexible TA and CA dinucleotide steps exhibited the highest Hoogsteen populations, and their kinetics rates strongly depended on the 3'-neighbor. In contrast, the stiffer AA and GA steps had the lowest Hoogsteen population, and their kinetics were weakly dependent on the 3'-neighbor. The Hoogsteen lifetime was especially short when G-C neighbors flanked the A-T base pair. The Hoogsteen dynamics had a distinct sequence-dependence compared to duplex stability and minor groove width. Thus, our results uncover a unique source of sequence-specificity hidden within the DNA double helix in the form of A-T Hoogsteen dynamics and establish the utility of 1H CEST to quantitively measure sequence-dependent DNA dynamics.

Systematic assessment of the flexibility of uracil damaged DNA

Uracil is a common DNA lesion which is recognized and removed by uracil DNA-glycosylase (UDG) as a part of the base excision repair pathway. Excision proceeds by base flipping, and UDG efficiency is thought to depend on the ease of deformability of the bases neighboring the lesion. We used molecular dynamics simulations to assess the flexibility of a large library of dsDNA strands, containing all tetranucleotide motifs with U:A, U:G, T:A or C:G base pairs. Our study demonstrates that uracil damaged DNA largely follows trends in flexibility of undamaged DNA. Measured bending persistence lengths, groove widths, step parameters and base flipping propensities demonstrate that uracil increases the flexibility of DNA, and that U:G base paired strands are more flexible than U:A strands. Certain sequence contexts are more deformable than others, with a key role for the 3' base next to uracil. Flexibilities are large when this base is an A or G, and repressed for a C or T. A 5' T adjacent to the uracil strongly promotes flexibility, but other 5' bases are less influential. DNA bending is correlated to step deformations and base flipping, and bending aids flipping. Our study implies that the link between substrate flexibility and UDG efficiency is widely valid, helps explain why UDG prefers to bind U:G base paired strands, and suggests that the DNA bending angle of the UDG-substrate complex is optimal for base flipping.Communicated by Ramaswamy H. Sarma.

Site-Specific Investigation of DNA Holliday Junction Dynamics and Structure with 6-Methylisoxanthopterin, a Fluorescent Guanine Analog

DNA Holliday Junction (HJ) formation and resolution is requisite for maintaining genomic stability in processes such as replication fork reversal and double-strand break repair. If HJs are not resolved, chromosome disjunction and aneuploidy result, hallmarks of tumor cells. To understand the structural features that lead to processing of these four-stranded joint molecule structures, we seek to identify structural and dynamic features unique to the central junction core. We incorporate the fluorescent guanine analog 6-methylisoxanthopterin (6-MI) at ten different locations throughout a model HJ structure to obtain site-specific information regarding the structure and dynamics of bases relative to those in a comparable sequence context in duplex DNA. These comparisons were accomplished through measuring fluorescence lifetime, relative brightness, fluorescence anisotropy, and thermodynamic stability, along with fluorescence quenching assays. These time-resolved and steady-state fluorescence measurements demonstrate that the structural distortions imposed by strand crossing result in increased solvent exposure, less stacking of bases and greater extrahelical nature of bases within the junction core. The 6-MI base analogs in the junction reflect these structural changes through an increase in intensity relative to those in the duplex. Molecular dynamics simulations performed using a model HJ indicate the primary sources of deformation are in the shift and twist parameters of the bases at the central junction step. These results suggest that junction-binding proteins may use the unique structure and dynamics of the bases at the core for recognition.

Rugged fitness landscapes minimize promiscuity in the evolution of transcriptional repressors

How a protein's function influences the shape of its fitness landscape, smooth or rugged, is a fundamental question in evolutionary biochemistry. Smooth landscapes arise when incremental mutational steps lead to a progressive change in function, as commonly seen in enzymes and binding proteins. On the other hand, rugged landscapes are poorly understood because of the inherent unpredictability of how sequence changes affect function. Here, we experimentally characterize the entire sequence phylogeny, comprising 1,158 extant and ancestral sequences, of the DNA-binding domain (DBD) of the LacI/GalR transcriptional repressor family. Our analysis revealed an extremely rugged landscape with rapid switching of specificity, even between adjacent nodes. Further, the ruggedness arises due to the necessity of the repressor to simultaneously evolve specificity for asymmetric operators and disfavors potentially adverse regulatory crosstalk. Our study provides fundamental insight into evolutionary, molecular, and biophysical rules of genetic regulation through the lens of fitness landscapes.

Structural insights into the morpholino nucleic acid/RNA duplex using the new XNA builder Ducque in a molecular modeling pipeline

The field of synthetic nucleic acids with novel backbone structures [xenobiotic nucleic acids (XNAs)] has flourished due to the increased importance of XNA antisense oligonucleotides and aptamers in medicine, as well as the development of XNA processing enzymes and new XNA genetic materials. Molecular modeling on XNA structures can accelerate rational design in the field of XNAs as it contributes in understanding and predicting how changes in the sugar-phosphate backbone impact on the complementation properties of the nucleic acids. To support the development of novel XNA polymers, we present a first-in-class open-source program (Ducque) to build duplexes of nucleic acid analogs with customizable chemistry. A detailed procedure is described to extend the Ducque library with new user-defined XNA fragments using quantum mechanics (QM) and to generate QM-based force field parameters for molecular dynamics simulations within standard packages such as AMBER. The tool was used within a molecular modeling workflow to accurately reproduce a selection of experimental structures for nucleic acid duplexes with ribose-based as well as non-ribose-based nucleosides. Additionally, it was challenged to build duplexes of morpholino nucleic acids bound to complementary RNA sequences.

Site-specific incorporation of a fluorescent nucleobase analog enhances i-motif stability and allows monitoring of i-motif folding inside cells

The i-motif is an intriguing non-canonical DNA structure, whose role in the cell is still controversial. Development of methods to study i-motif formation under physiological conditions in living cells is necessary to study its potential biological functions. The cytosine analog 1,3-diaza-2-oxophenoxazine (tCO) is a fluorescent nucleobase able to form either hemiprotonated base pairs with cytosine residues, or neutral base pairs with guanines. We show here that when tCO is incorporated in the proximity of a G:C:G:C minor groove tetrad, it induces a strong thermal and pH stabilization, resulting in i-motifs with Tm of 39ºC at neutral pH. The structural determination by NMR methods reveals that the enhanced stability is due to a large stacking interaction between the guanines of the tetrad with the tCO nucleobase, which forms a tCO:C+ in the folded structure at unusually-high pHs, leading to an increased quenching in its fluorescence at neutral conditions. This quenching is much lower when tCO is base-paired to guanines and totally disappears when the oligonucleotide is unfolded. By taking profit of this property, we have been able to monitor i-motif folding in cells.

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.

Structural and dynamical aspect of DNA motif sequence specific binding of AP-1 transcription factor

Activator protein-1 (AP-1) comprises one of the largest and most evolutionary conserved families of ubiquitous eukaryotic transcription factors that act as a pioneer factor. Diversity in DNA binding interaction of AP-1 through a conserved basic-zipper (bZIP) domain directs in-depth understanding of how AP-1 achieves its DNA binding selectivity and consequently gene regulation specificity. Here, we address the structural and dynamical aspects of the DNA target recognition process of AP-1 using microsecond-long atomistic simulations based on the structure of the human AP-1 FosB/JunD bZIP-DNA complex. Our results show the unique role of DNA shape features in selective base specific interactions, characteristic ion population, and solvation properties of DNA grooves to form the motif sequence specific AP-1-DNA complex. The TpG step at the two terminals of the AP-1 site plays an important role in the structural adjustment of DNA by modifying the helical twist in the AP-1 bound state. We addressed the role of intrinsic motion of the bZIP domain in terms of opening and closing gripper motions of DNA binding helices, in target site recognition and binding of AP-1 factors. Our observations suggest that binding to the cognate motif in DNA is mainly accompanied with the precise adjustment of closing gripper motion of DNA binding helices of the bZIP domain.

Structure-based prediction and characterization of photo-crosslinking in native protein-RNA complexes

UV-crosslinking of protein and RNA in direct contacts has been widely used to study protein-RNA complexes while our understanding of the photo-crosslinking mechanisms remains poor. This knowledge gap is due to the challenge of precisely mapping the crosslink sites in protein and RNA simultaneously in their native sequence and structural contexts. Here we systematically analyze protein-RNA interactions and photo-crosslinking by bridging crosslinked nucleotides and amino acids mapped using different assays with protein-RNA complex structures. We developed a computational method PxR3D-map which reliably predicts crosslink sites using structural information characterizing protein-RNA interaction interfaces. Analysis of the informative features revealed that photo-crosslinking is facilitated by base stacking with not only aromatic residues, but also dipeptide bonds that involve glycine, and distinct mechanisms are utilized by different RNA-binding domains. Our work suggests protein-RNA photo-crosslinking is highly selective in the cellular environment, which can guide data interpretation and further technology development for UV-crosslinking-based assays.

Systematic Comparison of Atomistic Force Fields for the Mechanical Properties of Double-Stranded DNA

The response of double-stranded DNA to external mechanical stress plays a central role in its interactions with the protein machinery in the cell. Modern atomistic force fields have been shown to provide highly accurate predictions for the fine structural features of the duplex. In contrast, and despite their pivotal function, less attention has been devoted to the accuracy of the prediction of the elastic parameters. Several reports have addressed the flexibility of double-stranded DNA via all-atom molecular dynamics, yet the collected information is insufficient to have a clear understanding of the relative performance of the various force fields. In this work, we fill this gap by performing a systematic study in which several systems, characterized by different sequence contexts, are simulated with the most popular force fields within the AMBER family, bcs1 and OL15, as well as with CHARMM36. Analysis of our results, together with their comparison with previous work focused on bsc0, allows us to unveil the differences in the predicted rigidity between the newest force fields and suggests a roadmap to test their performance against experiments. In the case of the stretch modulus, we reconcile these differences, showing that a single mapping between sequence-dependent conformation and elasticity via the crookedness parameter captures simultaneously the results of all force fields, supporting the key role of crookedness in the mechanical response of double-stranded DNA.

Enhanced sensitivity in electrochemical detection of ochratoxin A within food samples using ferrocene- and aptamer-tethered gold nanoparticles on disposable electrodes

Ensuring food security is crucial for public health, and the presence of mycotoxins, produced by fungi in improperly stored processed or unprocessed food, poses a significant threat. This research introduces a novel approach - a disposable aptasensing platform designed for the detection of ochratoxin A (OTA). The platform employs gold-nanostructured screen-printed carbon electrodes functionalized with a ferrocene derivative, serving as an integrated faradaic transducing system, and an anti-OTA aptamer as a bioreceptor site. Detection relies on the ferrocene electrochemical signal changes induced by the aptamer folding in the presence of the target molecule. Remarkably sensitive, the platform detects OTA within the range of 0.5 to 70 ng mL-1 and a detection limit of 11 pg mL-1. This limit is approximately 200 times below the levels stipulated by the European Commission for agricultural commodities. Notably, the sensing device exhibits efficacy in detecting OTA in complex media, such as roasted coffee beans and wine, without the need for sample pretreatment, yielding accurate recoveries. Furthermore, while label-free electrochemical aptasensors have proliferated, this study addresses a gap in understanding the binding mechanisms of some aptasensors. To enhance the experimental findings, a theoretical study was conducted to underscore the specificity of the anti-OTA aptamer as a donor for OTA detection. The molecular docking technique was employed to unveil the key binding region of the aptamer, providing valuable insights into the aptasensor specificity.

The structure and physical properties of a packaged bacteriophage particle

A string of nucleotides confined within a protein capsid contains all the instructions necessary to make a functional virus particle, a virion. Although the structure of the protein capsid is known for many virus species1,2, the three-dimensional organization of viral genomes has mostly eluded experimental probes3,4. Here we report all-atom structural models of an HK97 virion5, including its entire 39,732 base pair genome, obtained through multiresolution simulations. Mimicking the action of a packaging motor6, the genome was gradually loaded into the capsid. The structure of the packaged capsid was then refined through simulations of increasing resolution, which produced a 26 million atom model of the complete virion, including water and ions confined within the capsid. DNA packaging occurs through a loop extrusion mechanism7 that produces globally different configurations of the packaged genome and gives each viral particle individual traits. Multiple microsecond-long all-atom simulations characterized the effect of the packaged genome on capsid structure, internal pressure, electrostatics and diffusion of water, ions and DNA, and revealed the structural imprints of the capsid onto the genome. Our approach can be generalized to obtain complete all-atom structural models of other virus species, thereby potentially revealing new drug targets at the genome-capsid interface.

AMBERff at Scale: Multimillion-Atom Simulations with AMBER Force Fields in NAMD

All-atom molecular dynamics (MD) simulations are an essential structural biology technique with increasing application to multimillion-atom systems, including viruses and cellular machinery. Classical MD simulations rely on parameter sets, such as the AMBER family of force fields (AMBERff), to accurately describe molecular motion. Here, we present an implementation of AMBERff for use in NAMD that overcomes previous limitations to enable high-performance, massively parallel simulations encompassing up to two billion atoms. Single-point potential energy comparisons and case studies on model systems demonstrate that the implementation produces results that are as accurate as running AMBERff in its native engine.

Incorporating Sequence-Dependent DNA Shape and Dynamics into Transcriptome Data Analysis

Differentially expressed genes in a cellular context may be co-regulated by the same transcription factor. However, in the absence of a concurrent transcription factor binding data, such interactions are difficult to detect, especially at the single cell expression level. Motif enrichments in such genes can be used to gain insight into differential expressions caused by the shared upstream TFs. However, it is now established that many genes are co-regulated by the same TF due to a shared DNA shape or sequence-dependent conformational dynamics instead of sequence motif. In this work, we demonstrate how, starting from a gene expression data, such DNA shape and dynamics signatures can be potentially detected using publicly available tools, including DynaSeq, developed in our group for predicting the sequence-dependent components of these DNA shape features.

Insight into the on/off switch that regulates expression of the MSMEG-3762/63 efflux pump in Mycobacterium smegmatis

Drug resistance is one of the most difficult challenges facing tuberculosis (TB) control. Drug efflux is among the mechanisms leading to drug resistance. In our previous studies, we partially characterized the ABC-type MSMEG-3762/63 efflux pump in Mycobacterium smegmatis, which shares high percentage of identity with the Mycobacterium tuberculosis Rv1687/86c pump. MSMEG-3762/63 was shown to have extrusion activity for rifampicin and ciprofloxacin, used in first and second-line anti-TB treatments. Moreover, we described the functional role of the TetR-like MSMEG-3765 protein as a repressor of the MSMEG_3762/63/65 operon and orthologous Rv1687/86/85c in M. tuberculosis. Here we show that the operon is upregulated in the macrophage environment, supporting a previous observation of induction triggered by acid-nitrosative stress. Expression of the efflux pump was also induced by sub-inhibitory concentrations of rifampicin or ciprofloxacin. Both these drugs also prevented the binding of the MSMEG-3765 TetR repressor protein to its operator in the MSMEG_3762/63/65 operon. The hypothesis that these two drugs might be responsible for the induction of the efflux pump operon was assessed by bioinformatics analyses. Docking studies using a structural model of the regulator MSMEG-3765 showed that both antibiotics abolished the ability of this transcriptional repressor to recognize the efflux pump operon by interacting with the homodimer at different binding sites within the same binding pocket. Reduced binding of the repressor leads to induction of the efflux pump in M. smegmatis, and reduced efficacy of these two anti-mycobacterial drugs.

DNA damage by radiation as a function of electron energy and interaction at the atomic level with Monte Carlo simulation

In radiotherapy, X-ray or heavy ion beams target tumors to cause damage to their cell DNA. This damage is mainly induced by secondary low energy electrons. In this paper, we report the DNA molecular breaks at the atomic level as a function of electron energy and types of electron interactions using of Monte Carlo simulation. The number of DNA single and double strand breaks are compared to those from experimental results based on electron energies. In recent years, DNA atomistic models were introduced but still the simulations consider energy deposition in volumes of DNA or water equivalent material. We simulated a model of atomistic B-DNA in vacuum, forming 1122 base pairs of 30 nm in length. Each atom has been represented by a sphere whose radius equals the radius of van der Waals. We repeatedly simulated 10 million electrons for each energy from 4 eV to 500 eV and counted each interaction type with its position x, y, z in the volume of DNA. Based on the number and types of interactions at the atomic level, the number of DNA single and double strand breaks were calculated. We found that the dissociative electron attachment has the dominant effect on DNA strand breaks at energies below 10 eV compared to excitation and ionization. In addition, it is straightforward with our simulation to discriminate the strand and base breaks as a function of radiation interaction type and energy. In conclusion, the knowledge of DNA damage at the atomic level helps design direct internal therapeutic agents of cancer treatment.

Hydroxyl Radical vs. One-Electron Oxidation Reactivities in an Alternating GC Double-Stranded Oligonucleotide: A New Type Electron Hole Stabilization

We examined the reaction of hydroxyl radicals (HO•) and sulfate radical anions (SO4•-), which is generated by ionizing radiation in aqueous solutions under anoxic conditions, with an alternating GC doubled-stranded oligodeoxynucleotide (ds-ODN), i.e., the palindromic 5'-d(GCGCGC)-3'. In particular, the optical spectra of the intermediate species and associated kinetic data in the range of ns to ms were obtained via pulse radiolysis. Computational studies by means of density functional theory (DFT) for structural and time-dependent DFT for spectroscopic features were performed on 5'-d(GCGC)-3'. Comprehensively, our results suggest the addition of HO• to the G:C pair moiety, affording the [8-HO-G:C]• detectable adduct. The previous reported spectra of one-electron oxidation of a variety of ds-ODN were assigned to [G(-H+):C]• after deprotonation. Regarding 5'-d(GCGCGC)-3' ds-ODN, the spectrum at 800 ns has a completely different spectral shape and kinetic behavior. By means of calculations, we assigned the species to [G:C/C:G]•+, in which the electron hole is predicted to be delocalized on the two stacked base pairs. This transient species was further hydrated to afford the [8-HO-G:C]• detectable adduct. These remarkable findings suggest that the double-stranded alternating GC sequences allow for a new type of electron hole stabilization via delocalization over the whole sequence or part of it.

Asymmetric Periodic Boundary Conditions for All-Atom Molecular Dynamics and Coarse-Grained Simulations of Nucleic Acids

Periodic boundary conditions are commonly applied in molecular dynamics simulations in the microcanonical (NVE), canonical (NVT), and isothermal-isobaric (NpT) ensembles. In their simplest application, a biological system of interest is placed in the middle of a solvation box, which is chosen 'sufficiently large' to minimize any numerical artifacts associated with the periodic boundary conditions. This practical approach brings limitations to the size of biological systems that can be simulated. Here, we study simulations of effectively infinitely long nucleic acids, which are solvated in the directions perpendicular to the polymer chain, while periodic boundary conditions are also applied along the polymer chain. We study the effects of these asymmetric periodic boundary conditions (APBC) on the simulated results, including the mechanical properties of biopolymers and the properties of the surrounding solvent. To get some further insights into the advantages of using the APBC, a coarse-grained worm-like chain model is first studied, illustrating how the persistence length can be extracted from the local properties of the polymer chain, which are less affected by the APBC than some global averages. This is followed by all-atom molecular dynamics simulations of DNA in ionic solutions, where we use the APBC to investigate sequence-dependent properties of DNA molecules and properties of the surrounding solvent.

Modular Protein-DNA Cocrystals as Precise, Programmable Assembly Scaffolds

High-precision nanomaterials to entrap DNA-binding molecules are sought after for applications such as controlled drug delivery and scaffold-assisted structural biology. Here, we engineered protein-DNA cocrystals to serve as scaffolds for DNA-binding molecules. The designed cocrystals, isoreticular cocrystals, contain DNA-binding protein and cognate DNA blocks where the DNA-DNA junctions stack end-to-end. Furthermore, the crystal symmetry allows topology preserving (isoreticular) expansion of the DNA stack without breaking protein-protein contacts, hence providing larger solvent channels for guest diffusion. Experimentally, the resulting designed isoreticular cocrystal adopted an interpenetrating I222 lattice, a phenomenon previously observed in metal-organic frameworks (MOFs). The interpenetrating lattice crystallized dependably in the same space group despite myriad modifications at the DNA-DNA junctions. Assembly was modular with respect to the DNA inserted for expansion, providing an interchangeable DNA sequence for guest-specified scaffolding. Also, the DNA-DNA junctions were tunable, accommodating varied sticky base overhang lengths and terminal phosphorylation. As a proof of concept, we used the interpenetrating scaffold crystals to separately entrap three distinct guest molecules during crystallization. Isoreticular cocrystal design offers a route to a programmable scaffold for DNA-binding molecules, and the design principles may be applied to existing cocrystals to develop scaffolding materials.

Conformational dynamics and mechanical properties of biomimetic RNA, DNA, and RNA-DNA hybrid nanotubes: an atomistic molecular dynamics study

With the nanotechnology boom, artificially designed nucleic acid nanotubes have aroused interest due to their practical applications in nanorobotics, vaccine design, membrane channels, drug delivery, and force sensing. In this paper, computational study was performed to investigate the structural dynamics and mechanical properties of RNA nanotubes (RNTs), DNA nanotubes (DNTs), and RNA-DNA hybrid nanotubes (RDHNTs). So far, the structural and mechanical properties of RDHNTs have not been examined in experiments or theoretical calculations, and there is limited knowledge regarding these properties for RNTs. Here, the simulations were carried out using the equilibrium molecular dynamics (MD) and steered molecular dynamics (SMD) approaches. Using in-house scripting, we modeled hexagonal nanotubes composed of six double-stranded molecules connected by four-way Holliday junctions. Classical MD analyses were performed on the collected trajectory data to investigate structural properties. Analyses of the microscopic structural parameters of RDHNT indicated a structural transition from the A-form to a conformation between the A- and B-forms, which may be attributable to the increased rigidity of RNA scaffolds compared to DNA staples. Comprehensive research on the elastic mechanical properties was also conducted based on spontaneous thermal fluctuations of nanotubes and employing the equipartition theorem. The Young's modulus of RDHNT (E = 165 MPa) and RNT (E = 144 MPa) was found to be almost the same and nearly half of that found for DNT (E = 325 MPa). Furthermore, the results showed that RNT was more resistant to bending, torsional, and volumetric deformations than DNT and RDHNT. We also used non-equilibrium SMD simulations to acquire comprehensive knowledge of the mechanical response of nanotubes to tensile stress.

Molecular Mechanism of Nucleosome Recognition by the Pioneer Transcription Factor Sox

Pioneer transcription factors (PTFs) have the remarkable ability to directly bind to chromatin to stimulate vital cellular processes. In this work, we dissect the universal binding mode of Sox PTF by combining extensive molecular simulations and physiochemistry approaches, along with DNA footprinting techniques. As a result, we show that when Sox consensus DNA is located at the solvent-facing DNA strand, Sox binds to the compact nucleosome without imposing any significant conformational changes. We also reveal that the base-specific Sox:DNA interactions (base reading) and Sox-induced DNA changes (shape reading) are concurrently required for sequence-specific nucleosomal DNA recognition. Among three different nucleosome positions located on the positive DNA arm, a sequence-specific reading mechanism is solely satisfied at the superhelical location 2 (SHL2). While SHL2 acts transparently for solvent-facing Sox binding, among the other two positions, SHL4 permits only shape reading. The final position, SHL0 (dyad), on the other hand, allows no reading mechanism. These findings demonstrate that Sox-based nucleosome recognition is essentially guided by intrinsic nucleosome properties, permitting varying degrees of DNA recognition.

Uracil-DNA glycosylase efficiency is modulated by substrate rigidity

Uracil DNA-glycosylase (UNG) is a DNA repair enzyme that removes the highly mutagenic uracil lesion from DNA using a base flipping mechanism. Although this enzyme has evolved to remove uracil from diverse sequence contexts, UNG excision efficiency depends on DNA sequence. To provide the molecular basis for rationalizing UNG substrate preferences, we used time-resolved fluorescence spectroscopy, NMR imino proton exchange measurements, and molecular dynamics simulations to measure UNG specificity constants (kcat/KM) and DNA flexibilities for DNA substrates containing central AUT, TUA, AUA, and TUT motifs. Our study shows that UNG efficiency is dictated by the intrinsic deformability around the lesion, establishes a direct relationship between substrate flexibility modes and UNG efficiency, and shows that bases immediately adjacent to the uracil are allosterically coupled and have the greatest impact on substrate flexibility and UNG activity. The finding that substrate flexibility controls UNG efficiency is likely significant for other repair enzymes and has major implications for the understanding of mutation hotspot genesis, molecular evolution, and base editing.

Molecular dynamics simulation-based trinucleotide and tetranucleotide level structural and energy characterization of the functional units of genomic DNA

Genomes of most organisms on earth are written in a universal language of life, made up of four units - adenine (A), thymine (T), guanine (G), and cytosine (C), and understanding the way they are put together has been a great challenge to date. Multiple efforts have been made to annotate this wonderfully engineered string of DNA using different methods but they lack a universal character. In this article, we have investigated the structural and energetic profiles of both prokaryotes and eukaryotes by considering two essential genomic sites, viz., the transcription start sites (TSS) and exon-intron boundaries. We have characterized these sites by mapping the structural and energy features of DNA obtained from molecular dynamics simulations, which considers all possible trinucleotide and tetranucleotide steps. For DNA, these physicochemical properties show distinct signatures at the TSS and intron-exon boundaries. Our results firmly convey the idea that DNA uses the same dialect for prokaryotes and eukaryotes and that it is worth going beyond sequence-level analyses to physicochemical space to determine the functional destiny of DNA sequences.

Crystal structures of 'ALternative Isoinformational ENgineered' DNA in B-form

The first structural model of duplex DNA reported in 1953 by Watson & Crick presented the double helix in B-form, the form that genomic DNA exists in much of the time. Thus, artificial DNA seeking to mimic the properties of natural DNA should also be able to adopt B-form. Using a host-guest system in which Moloney murine leukemia virus reverse transcriptase serves as the host and DNA as the guests, we determined high-resolution crystal structures of three complexes including 5'-CTTBPPBBSSZZSAAG, 5'-CTTSSPBZPSZBBAAG and 5'-CTTZZPBSBSZPPAAG with 10 consecutive unnatural nucleobase pairs in B-form within self-complementary 16 bp duplex oligonucleotides. We refer to this ALternative Isoinformational ENgineered (ALIEN) genetic system containing two nucleobase pairs (P:Z, pairing 2-amino-imidazo-[1,2-a]-1,3,5-triazin-(8H)-4-one with 6-amino-5-nitro-(1H)-pyridin-2-one, and B:S, 6-amino-4-hydroxy-5-(1H)-purin-2-one with 3-methyl-6-amino-pyrimidin-2-one) as ALIEN DNA. We characterized both position- and sequence-specific helical, nucleobase pair and dinucleotide step parameters of P:Z and B:S pairs in the context of B-form DNA. We conclude that ALIEN DNA exhibits structural features that vary with sequence. Further, Z can participate in alternative stacking modes within a similar sequence context as captured in two different structures. This finding suggests that ALIEN DNA may have a larger repertoire of B-form structures than natural DNA. This article is part of the theme issue 'Reactivity and mechanism in chemical and synthetic biology'.

i-Motif folding intermediates with zero-nucleotide loops are trapped by 2'-fluoroarabinocytidine via F···H and O···H hydrogen bonds

G-quadruplex and i-motif nucleic acid structures are believed to fold through kinetic partitioning mechanisms. Such mechanisms explain the structural heterogeneity of G-quadruplex metastable intermediates which have been extensively reported. On the other hand, i-motif folding is regarded as predictable, and research on alternative i-motif folds is limited. While TC₅ normally folds into a stable tetrameric i-motif in solution, we report that 2'-deoxy-2'-fluoroarabinocytidine (araF-C) substitutions can prompt TC₅ to form an off-pathway and kinetically-trapped dimeric i-motif, thereby expanding the scope of i-motif folding landscapes. This i-motif is formed by two strands, associated head-to-head, and featuring zero-nucleotide loops which have not been previously observed. Through spectroscopic and computational analyses, we also establish that the dimeric i-motif is stabilized by fluorine and non-fluorine hydrogen bonds, thereby explaining the superlative stability of araF-C modified i-motifs. Comparative experimental findings suggest that the strength of these interactions depends on the flexible sugar pucker adopted by the araF-C residue. Overall, the findings reported here provide a new role for i-motifs in nanotechnology and also pose the question of whether unprecedented i-motif folds may exist in vivo.

pH-Dependent Capping Interactions Induce Large-Scale Structural Transitions in i-Motifs

We study here a DNA oligonucleotide having the ability to form two different i-motif structures whose relative stability depends on pH and temperature. The major species at neutral pH is stabilized by two C:C+ base pairs capped by two minor groove G:C:G:C tetrads. The high pH and thermal stability of this structure are mainly due to the favorable effect of the minor groove tetrads on their adjacent positively charged C:C+ base pairs. At pH 5, we observe a more elongated i-motif structure consisting of four C:C+ base pairs capped by two G:T:G:T tetrads. Molecular dynamics calculations show that the conformational transition between the two structures is driven by the protonation state of key cytosines. In spite of large conformational differences, the transition between the acidic and neutral structures can occur without unfolding of the i-motif. These results represent the first case of a conformational switch between two different i-motif structures and illustrate the dramatic pH-dependent plasticity of this fascinating DNA motif.

3D RNA-scaffolded wireframe origami

Hybrid RNA:DNA origami, in which a long RNA scaffold strand folds into a target nanostructure via thermal annealing with complementary DNA oligos, has only been explored to a limited extent despite its unique potential for biomedical delivery of mRNA, tertiary structure characterization of long RNAs, and fabrication of artificial ribozymes. Here, we investigate design principles of three-dimensional wireframe RNA-scaffolded origami rendered as polyhedra composed of dual-duplex edges. We computationally design, fabricate, and characterize tetrahedra folded from an EGFP-encoding messenger RNA and de Bruijn sequences, an octahedron folded with M13 transcript RNA, and an octahedron and pentagonal bipyramids folded with 23S ribosomal RNA, demonstrating the ability to make diverse polyhedral shapes with distinct structural and functional RNA scaffolds. We characterize secondary and tertiary structures using dimethyl sulfate mutational profiling and cryo-electron microscopy, revealing insight into both global and local, base-level structures of origami. Our top-down sequence design strategy enables the use of long RNAs as functional scaffolds for complex wireframe origami.

Structural and Energetic Features of Base-Base Stacking Contacts in RNA

Nucleobase π-π stacking is one of the crucial organizing interactions within three-dimensional (3D) RNA architectures. Characterizing the structural variability of these contacts in RNA crystal structures will help delineate their subtleties and their role in determining function. This analysis of different stacking geometries found in RNA X-ray crystal structures is the largest such survey to date; coupled with quantum-mechanical calculations on typical representatives of each possible stacking arrangement, we determined the distribution of stacking interaction energies. A total of 1,735,481 stacking contacts, spanning 359 of the 384 theoretically possible distinct stacking geometries, were identified. Our analysis reveals preferential occurrences of specific consecutive stacking arrangements in certain regions of RNA architectures. Quantum chemical calculations suggest that 88 of the 359 contacts possess intrinsically stable stacking geometries, whereas the remaining stacks require the RNA backbone or surrounding macromolecular environment to force their formation and maintain their stability. Our systematic analysis of π-π stacks in RNA highlights trends in the occurrence and localization of these noncovalent interactions and may help better understand the structural intricacies of functional RNA-based molecular architectures.

Site-specific investigation of DNA Holliday Junction dynamics and structure with 6-Methylisoxanthopterin, a fluorescent guanine analog

DNA Holliday Junction (HJ) formation and resolution is requisite for maintaining genomic stability in processes such as replication fork reversal and double-strand break repair. If HJs are not resolved, chromosome disjunction and aneuploidy result, hallmarks of tumor cells. To understand the structural features that lead to processing of these four-stranded joint molecule structures, we seek to identify structural and dynamic features unique to the central junction core. We incorporated the fluorescent guanine analog 6-methylisoxanthopterin (6-MI) at ten different locations throughout a model HJ structure to obtain site-specific information regarding the structure and dynamics of bases relative to those in a comparable sequence context in duplex DNA. These comparisons were accomplished through measuring fluorescence lifetime, relative brightness, fluorescence anisotropy, and quenching assays. These time-resolved and steady-state fluorescence measurements demonstrate that the structural distortions imposed by strand crossing result in increased solvent exposure, less stacking of bases and greater extrahelical nature of bases within the junction core. The 6-MI base analogs in the junction reflect these structural changes through an increase in intensity relative to those in the duplex. Molecular dynamics simulations performed using a model HJ indicate that the primary sources of deformation are in the shift and twist parameters of the bases at the central junction step. These results suggest that junction-binding proteins may use the unique structure and dynamics of the bases at the core for recognition.

StackCirRNAPred: computational classification of long circRNA from other lncRNA based on stacking strategy

BACKGROUND

CircRNAs are essential for the regulation of post-transcriptional gene expression, including as miRNA sponges, and play an important role in disease development. Some computational tools have been proposed recently to predict circRNA, since only one classifier is used, there is still much that can be done to improve the performance.

RESULTS

StackCirRNAPred was proposed, the computational classification of long circRNA from other lncRNA based on stacking strategy. In order to cope with the potential problem that a single feature might not be able to distinguish circRNA well from other lncRNA, we first extracted features from different sources, including nucleic acid composition, sequence spatial features and physicochemical properties, Alu and tandem repeats. We innovatively apply the stacking strategy to integrate the more advantageous classifiers of RF, LightGBM, XGBoost. This allows the model to incorporate these features more flexibly. StackCirRNAPred was found to be significantly better than other tools, with precision, accuracy, F1, recall and MCC of 0.843, 0.833, 0.831, 0.819 and 0.666 respectively. We tested it directly on the mouse dataset. StackCirRNAPred was still significantly better than other methods, with precision, accuracy, F1, recall and MCC of 0.837, 0.839, 0.839, 0.841, 0.677.

CONCLUSIONS

We proposed StackCirRNAPred based on stacking strategy to distinguish long circRNAs from other lncRNAs. With the test results demonstrating the validity and robustness of StackCirRNAPred, we hope StackCirRNAPred will complement existing circRNA prediction methods and is helpful in down-stream research.

Duplex DNA Retains the Conformational Features of Single Strands: Perspectives from MD Simulations and Quantum Chemical Computations

Molecular dynamics simulations and geometry optimizations carried out at the quantum level as well as by quantum mechanical/molecular mechanics methods predict that short, single-stranded DNA oligonucleotides adopt conformations very similar to those observed in crystallographic double-stranded B-DNA, with rise coordinates close to ≈3.3 Å. In agreement with the experimental evidence, the computational results show that DNA single strands rich in adjacent purine nucleobases assume more regular arrangements than poly-thymine. The preliminary results suggest that single-stranded poly-cytosine DNA should also retain a substantial helical order in solution. A comparison of the structures of single and double helices confirms that the B-DNA motif is a favorable arrangement also for single strands. Indeed, the optimal geometry of the complementary single helices is changed to a very small extent in the formation of the duplex.

Sequence-Specific Structural Features and Solvation Properties of Transcription Factor Binding DNA Motifs: Insights from Molecular Dynamics Simulation

Sequence-specific recognition of transcription factor (TF) binding motifs in the target site of DNA over the vast amount of non-target DNA is of primary importance for the transcriptional regulation of gene expression by the TFs. Binding of TFs to the target site of DNA relies not only on the direct contact formation but also on the structural and conformational features of DNA. Recognition of DNA structural features or shape readout by proteins is an important factor in the context of TF-DNA interaction. Based on the atomistic molecular simulation, here we report the sequence-dependent unique structural features, solvation, and ion-binding properties of biologically relevant AT- and GC-rich human TF binding motifs in DNA. Counterion and water distribution around the motif is found to be sensitive to the motif sequence, which is accompanied with the DNA shape features. The motif sequence affects the electrostatic potential along the grooves, and cytosine methylation alters the DNA shape features. Characteristic solvation properties of TF binding motif DNA fragments infer that an ionic environment and hydration influences are essential to describe TF-DNA interactions.

Structures of a mobile intron retroelement poised to attack its structured DNA target

Group II introns are ribozymes that catalyze their self-excision and function as retroelements that invade DNA. As retrotransposons, group II introns form ribonucleoprotein (RNP) complexes that roam the genome, integrating by reversal of forward splicing. Here we show that retrotransposition is achieved by a tertiary complex between a structurally elaborate ribozyme, its protein mobility factor, and a structured DNA substrate. We solved cryo-electron microscopy structures of an intact group IIC intron-maturase retroelement that was poised for integration into a DNA stem-loop motif. By visualizing the RNP before and after DNA targeting, we show that it is primed for attack and fits perfectly with its DNA target. This study reveals design principles of a prototypical retroelement and reinforces the hypothesis that group II introns are ancient elements of genetic diversification.

Gold nanoparticle design for RNA compaction

RNA-based therapeutics hold a great promise in treating a variety of diseases. However, double-stranded RNAs (dsRNAs) are inherently unstable, highly charged, and stiff macromolecules that require a delivery vehicle. Cationic ligand functionalized gold nanoparticles (AuNPs) are able to compact nucleic acids and assist in RNA delivery. Here, we use large-scale all-atom molecular dynamics simulations to show that correlations between ligand length, metal core size, and ligand excess free volume control the ability of nanoparticles to bend dsRNA far below its persistence length. The analysis of ammonium binding sites showed that longer ligands that bind deep within the major groove did not cause bending. By limiting ligand length and, thus, excess free volume, we have designed nanoparticles with controlled internal binding to RNA's major groove. NPs that are able to induce RNA bending cause a periodic variation in RNA's major groove width. Density functional theory studies on smaller models support large-scale simulations. Our results are expected to have significant implications in packaging of nucleic acids for their applications in nanotechnology and gene delivery.

HDOCK update for modeling protein-RNA/DNA complex structures

Protein-nucleic acid interactions are involved in various cellular processes. Therefore, determining the structures of protein-nucleic acid complexes can provide insights into the mechanisms of the interactions and thus guide the rational drug design to modulate these interactions. Due to the high cost and technical difficulties of solving complex structures experimentally, computational modeling such as molecular docking has been playing an important role in the study of molecular interactions. In order to make it easier for researchers to obtain biomolecular complex structures through molecular docking, we developed the HDOCK server for protein-protein and protein-RNA/DNA docking (accessed at http://hdock.phys.hust.edu.cn/). Since its first release in 2017, HDOCK has been widely used in the scientific community. As nucleic acids may include single-stranded (ss) RNA/DNA and double-stranded (ds) RNA/DNA, we now present an updated version of HDOCK, which offers new options for structural modeling of ssRNA, ssDNA, dsRNA, and dsDNA. We hope this update will better help the scientific community solve important biological problems, thereby advancing the field. In this article, we describe the general protocol of HDOCK with emphasis on the new functions on RNA/DNA modeling. Several application examples are also given to illustrate the usage of the new functions.

Helicene Structure between DNA and Cyanuric Acid: The Role of Noncovalent Interactions

The self-assembly between DNA and small organic molecules can expand the structural space of and introduce novel functionalities to DNA nanomaterials. In particular, it was demonstrated that poly(adenosine) DNA self-assembles with cyanuric acid (CA) to form a triplex helical structure. Previous molecular dynamics simulations showed that the DNA-CA assemblies adopt a novel noncovalent helicene structure that has a continuous helical hydrogen bond network. This article explores why the assemblies adopt the helicene geometry instead of an alternative planar hexameric rosette geometry. Analysis of the hydrogen bonding and stacking interaction energies indicates that constraining the system to the hexameric rosette geometry strains the hydrogen bonds without significantly improving the interaction energy. Molecular dynamics simulations for the assemblies between adenosine nucleosides and CA confirm that the formation of helicene structure is primarily driven by base-pair interactions and not because of the DNA backbone.

Predicting accurate ab initio DNA electron densities with equivariant neural networks

One of the fundamental limitations of accurately modeling biomolecules like DNA is the inability to perform quantum chemistry calculations on large molecular structures. We present a machine learning model based on an equivariant Euclidean neural network framework to obtain accurate ab initio electron densities for arbitrary DNA structures that are much too large for conventional quantum methods. The model is trained on representative B-DNA basepair steps that capture both base pairing and base stacking interactions. The model produces accurate electron densities for arbitrary B-DNA structures with typical errors of less than 1%. Crucially, the error does not increase with system size, which suggests that the model can extrapolate to large DNA structures with negligible loss of accuracy. The model also generalizes reasonably to other DNA structural motifs such as the A- and Z-DNA forms, despite being trained on only B-DNA configurations. The model is used to calculate electron densities of several large-scale DNA structures, and we show that the computational scaling for this model is essentially linear. We also show that this machine learning electron density model can be used to calculate accurate electrostatic potentials for DNA. These electrostatic potentials produce more accurate results compared with classical force fields and do not show the usual deficiencies at short range.

pyDockDNA: A new web server for energy-based protein-DNA docking and scoring

Proteins and nucleic acids are essential biological macromolecules for cell life. Indeed, interactions between proteins and DNA regulate many biological processes such as protein synthesis, signal transduction, DNA storage, or DNA replication and repair. Despite their importance, less than 4% of total structures deposited in the Protein Data Bank (PDB) correspond to protein-DNA complexes, and very few computational methods are available to model their structure. We present here the pyDockDNA web server, which can successfully model a protein-DNA complex with a reasonable predictive success rate (as benchmarked on a standard dataset of protein-DNA complex structures, where DNA is in B-DNA conformation). The server implements the pyDockDNA program, as a module of pyDock suite, thus including third-party programs, modules, and previously developed tools, as well as new modules and parameters to handle the DNA properly. The user is asked to enter Protein Data Bank files for protein and DNA input structures (or suitable models) and select the chains to be docked. The server calculations are mainly divided into three steps: sampling by FTDOCK, scoring with new energy-based parameters and the possibility of applying external restraints. The user can select different options for these steps. The final output screen shows a 3D representation of the top 10 models and a table sorting the model according to the scoring function selected previously. All these output files can be downloaded, including the top 100 models predicted by pyDockDNA. The server can be freely accessed for academic use (https://model3dbio.csic.es/pydockdna).

Three's a crowd - stabilisation, structure, and applications of DNA triplexes

DNA is a strikingly flexible molecule and can form a variety of secondary structures, including the triple helix, which is the subject of this review. The DNA triplex may be formed naturally, during homologous recombination, or can be formed by the introduction of a synthetic triplex forming oligonucleotide (TFO) to a DNA duplex. As the TFO will bind to the duplex with sequence specificity, there is significant interest in developing TFOs with potential therapeutic applications, including using TFOs as a delivery mechanism for compounds able to modify or damage DNA. However, to combine triplexes with functionalised compounds, a full understanding of triplex structure and chemical modification strategies, which may increase triplex stability or in vivo degradation, is essential - these areas will be discussed in this review. Ruthenium polypyridyl complexes, which are able to photooxidise DNA and act as luminescent DNA probes, may serve as a suitable photophysical payload for a TFO system and the developments in this area in the context of DNA triplexes will also be reviewed.