Supervised Chromatin Loop Detection Using Peakachu Version 2.

Peakachu is a supervised-learning-based approach that identifies chromatin loops from chromatin contact data. Here, we present Peakachu version 2, an updated version that significantly improves extensibility, usability, and computational efficiency compared to its predecessor. It features pretrained models tailored for a wide range of experimental platforms, such as Hi-C, Micro-C, ChIA-PET, HiChIP, HiCAR, and TrAC-loop. This chapter offers a step-by-step tutorial guiding users through the process of training Peakachu models from scratch and utilizing pretrained models to predict chromatin loops across various platforms.

Mechanistic analysis of Riboswitch Ligand interactions provides insights into pharmacological control over gene expression.

Riboswitches are structured RNA elements that regulate gene expression upon binding to small molecule ligands. Understanding the mechanisms by which small molecules impact riboswitch activity is key to developing potent, selective ligands for these and other RNA targets. We report the structure-informed design of chemically diverse synthetic ligands for PreQ1 riboswitches. Multiple X-ray co-crystal structures of synthetic ligands with the Thermoanaerobacter tengcongensis (Tte)-PreQ1 riboswitch confirm a common binding site with the cognate ligand, despite considerable chemical differences among the ligands. Structure probing assays demonstrate that one ligand causes conformational changes similar to PreQ1 in six structurally and mechanistically diverse PreQ1 riboswitch aptamers. Single-molecule force spectroscopy is used to demonstrate differential modes of riboswitch stabilization by the ligands. Binding of the natural ligand brings about the formation of a persistent, folded pseudoknot structure, whereas a synthetic ligand decreases the rate of unfolding through a kinetic mechanism. Single round transcription termination assays show the biochemical activity of the ligands, while a GFP reporter system reveals compound activity in regulating gene expression in live cells without toxicity. Taken together, this study reveals that diverse small molecules can impact gene expression in live cells by altering conformational changes in RNA structures through distinct mechanisms.

Genetic polymorphism in untranslated regions of PRKCZ influences mRNA structure, stability and binding sites.

BACKGROUND: Variations in untranslated regions (UTR) alter regulatory pathways impacting phenotype, disease onset, and course of disease. Protein kinase C Zeta (PRKCZ), a serine-threonine kinase, is implicated in cardiovascular, neurological and oncological disorders. Due to limited research on PRKCZ, this study aimed to investigate the impact of UTR genetic variants' on binding sites for transcription factors and miRNA. RNA secondary structure, eQTLs, and variation tolerance analysis were also part of the study. METHODS: The data related to PRKCZ gene variants was downloaded from the Ensembl genome browser, COSMIC and gnomAD. The RegulomeDB database was used to assess the functional impact of 5' UTR and 3'UTR variants. The analysis of the transcription binding sites (TFBS) was done through the Alibaba tool, and the Kyoto Encyclopaedia of Genes and Genomes (KEGG) was employed to identify pathways associated with PRKCZ. To predict the effect of variants on microRNA binding sites, PolymiRTS was utilized for 3' UTR variants, and the SNPinfo tool was used for 5' UTR variants. RESULTS: The results obtained indicated that a total of 24 variants present in the 3' UTR and 25 variants present in the 5' UTR were most detrimental. TFBS analysis revealed that 5' UTR variants added YY1, repressor, and Oct1, whereas 3' UTR variants added AP-2alpha, AhR, Da, GR, and USF binding sites. The study predicted TFs that influenced PRKCZ expression. RNA secondary structure analysis showed that eight 5' UTR and six 3' UTR altered the RNA structure by either removal or addition of the stem-loop. The microRNA binding site analysis highlighted that seven 3' UTR and one 5' UTR variant altered the conserved site and also created new binding sites. eQTLs analysis showed that one variant was associated with PRKCZ expression in the lung and thyroid. The variation tolerance analysis revealed that PRKCZ was an intolerant gene. CONCLUSION: This study laid the groundwork for future studies aimed at targeting PRKCZ as a therapeutic target.

DNA domino circuits based on a hairpin exonuclease assistance signal transmission architecture for temporal logic operations.

DNA circuits are important fundamental tools for performing temporal logic operations with complex structures, but they lack sequence orthogonality. Here, we developed a simple and orthogonal hairpin exonuclease assistance signal (H-EAST) architecture to construct DNA domino circuits with time-delay characteristics and temporal logic operations, which has potential applications in biomolecular computing.

RNA dynamics from experimental and computational approaches.

Conformational dynamics is crucial for the biological function of RNA molecules and for their potential as therapeutic targets. This meeting report outlines key "take-home" messages that emerged from the presentations and discussions during the CECAM workshop "RNA dynamics from experimental and computational approaches" in Paris, June 26-28, 2023.

Beyond the "spine of hydration": Chiral SFG spectroscopy detects DNA first hydration shell and base pair structures.

Experimental methods capable of selectively probing water at the DNA minor groove, major groove, and phosphate backbone are crucial for understanding how hydration influences DNA structure and function. Chiral-selective sum frequency generation spectroscopy (chiral SFG) is unique among vibrational spectroscopies because it can selectively probe water molecules that form chiral hydration structures around biomolecules. However, interpreting chiral SFG spectra is challenging since both water and the biomolecule can produce chiral SFG signals. Here, we combine experiment and computation to establish a theoretical framework for the rigorous interpretation of chiral SFG spectra of DNA. We demonstrate that chiral SFG detects the N-H stretch of DNA base pairs and the O-H stretch of water, exclusively probing water molecules in the DNA first hydration shell. Our analysis reveals that DNA transfers chirality to water molecules only within the first hydration shell, so they can be probed by chiral SFG spectroscopy. Beyond the first hydration shell, the electric field-induced water structure is symmetric and, therefore, precludes chiral SFG response. Furthermore, we find that chiral SFG can differentiate chiral subpopulations of first hydration shell water molecules at the minor groove, major groove, and phosphate backbone. Our findings challenge the scientific perspective dominant for more than 40 years that the minor groove "spine of hydration" is the only chiral water structure surrounding the DNA double helix. By identifying the molecular origins of the DNA chiral SFG spectrum, we lay a robust experimental and theoretical foundation for applying chiral SFG to explore the chemical and biological physics of DNA hydration.

Sensing the structural and conformational properties of single-stranded nucleic acids using electrometry and molecular simulations.

Inferring the 3D structure and conformation of disordered biomolecules, e.g., single stranded nucleic acids (ssNAs), remains challenging due to their conformational heterogeneity in solution. Here, we use escape-time electrometry (ETe) to measure with sub elementary-charge precision the effective electrical charge in solution of short to medium chain length ssNAs in the range of 5-60 bases. We compare measurements of molecular effective charge with theoretically calculated values for simulated molecular conformations obtained from Molecular Dynamics simulations using a variety of forcefield descriptions. We demonstrate that the measured effective charge captures subtle differences in molecular structure in various nucleic acid homopolymers of identical length, and also that the experimental measurements can find agreement with computed values derived from coarse-grained molecular structure descriptions such as oxDNA, as well next generation ssNA force fields. We further show that comparing the measured effective charge with calculations for a rigid, charged rod-the simplest model of a nucleic acid-yields estimates of molecular structural dimensions such as linear charge spacings that capture molecular structural trends observed using high resolution structural analysis methods such as X-ray scattering. By sensitively probing the effective charge of a molecule, electrometry provides a powerful dimension supporting inferences of molecular structural and conformational properties, as well as the validation of biomolecular structural models. The overall approach holds promise for a high throughput, microscopy-based biomolecular analytical approach offering rapid screening and inference of molecular 3D conformation, and operating at the single molecule level in solution.

Solution-based biophysical characterization of conformation change in structure-switching aptamers.

Structure-switching aptamers have become ubiquitous in several applications, notably in analytical devices such as biosensors, due to their ease of supporting strong signaling. Aside from their ability to bind specifically with their respective target, this class of aptamers also undergoes a conformational rearrangement upon target recognition. While several well-studied and early-developed aptamers (e.g., cocaine, ATP, and thrombin) have been found to have this structure-switching property, the vast majority do not. As a result, it is common to try to engineer aptamers into switches. This proves challenging in part because of the difficulty in obtaining structural and functional information about aptamers. In response, we review various readily available biophysical characterization tools that are capable of assessing structure switching of aptamers. In doing so, we delve into the fundamentals of these different techniques and detail how they have been utilized in characterizing structure-switching aptamers. While each of these biophysical techniques alone has utility, their real power to demonstrate the occurrence of structural change with ligand binding is when multiple techniques are used. We hope that through a deeper understanding of these techniques, researchers will be better able to acquire biophysical information about their aptamer-ligand systems and accelerate the translation of aptamers into biosensors.

AlphaFold 3 - Aided Design of DNA Motifs To Assemble into Triangles.

Self-assembly of biomolecules provides a powerful tool for a wide range of applications in nanomedicine, biosensing and imaging, vaccines, computation, nanophotonics, etc. The key is to rationally design building blocks and the intermolecule interactions. Along this line, structural DNA nanotechnology has rapidly developed by limiting DNA secondary structures to primarily well-established, B-form DNA duplexes, which can be readily and reliably predicted. As the field evolves, more sophisticated structural elements must be introduced. While increasing the structural complexity, they bring challenges to predicting DNA nanostructures. In the past, a brutal and tedious error-and-trial approach has often been used to solve this problem. Here, we report a case study of applying AlphaFold 3 to model the structural elements to facilitate DNA nanostructure design. This protocol is expected to be generally applicable and greatly facilitates the further development of structural DNA nanotechnology.

Peptide nucleic acids can form hairpins and bind RNA-binding proteins.

RNA-binding proteins (RBPs) are a major class of proteins that interact with RNAs to change their fate or function. RBPs and the ribonucleoprotein complexes they constitute are involved in many essential cellular processes. In many cases, the molecular details of RBP:RNA interactions differ between viruses, prokaryotes and eukaryotes, making prokaryotic and viral RBPs good potential drug targets. However, targeting RBPs with small molecules has so far been met with limited success as RNA-binding sites tend to be extended, shallow and dynamic with a mixture of charged, polar and hydrophobic interactions. Here, we show that peptide nucleic acids (PNAs) with nucleic acid-like binding properties and a highly stable peptide-like backbone can be used to target some RBPs. We have designed PNAs to mimic the short RNA stem-loop sequence required for the initiation of prokaryotic signal recognition particle (SRP) assembly, a target for antibiotics development. Using a range of biophysical and biochemical assays, the designed PNAs were demonstrated to fold into a hairpin structure, bind the targeted protein and compete with the native RNA hairpin to inhibit SRP formation. To show the applicability of PNAs against other RBPs, a PNA was also shown to bind Nsp9 from SARS-CoV-2, a protein that exhibits non-sequence-specific RNA binding but preferentially binds hairpin structures. Taken together, our results support that PNAs can be a promising class of compounds for targeting RNA-binding activities in RBPs.

Isolation and genetic characterization of novel bovine parechoviruses from Japanese black cattle.

Novel bovine parechoviruses (Bo ParVs) were isolated from the feces of Japanese black cattle. Phylogenetic analysis revealed that the novel Bo ParVs formed an independent cluster, exhibiting 72.2-75.6% nucleotide sequence identity to previous Bo ParVs, suggesting that they represent a new genotype. Bo ParVs, including the novel Bo ParVs, shared sequence similarity with each other in the 3' untranslated region (3'UTR) and exhibited low sequence similarity (<38.9% identity) to other parechoviruses. However, a secondary structure prediction of the 3'UTR revealed that the Bo ParVs shared conserved motifs in domain 2 with parechovirus B and E, suggesting some evolutionary constrains in this region.

Molecular Dynamics Simulations of the miR-155 Duplex: Impact of Ionic Strength on Structure and Na+ and Cl- Ion Distribution.

MiR-155 is a multifunctional microRNA involved in many biological processes. Since miR-155 is overexpressed in several pathologies, its detection deserves high interest in clinical diagnostics. Biosensing approaches often exploit the hybridization of miR-155 with its complementary strand. Molecular Dynamics (MD) simulations were applied to investigate the complex formed by miR-155 and its complementary strand in aqueous solution with Na+ and Cl- ions at ionic strengths in the 100-400 mM range, conditions commonly used in biosensing experiments. We found that the main structural properties of the duplex are preserved at all the investigated ionic strengths. The radial distribution functions of both Na+ and Cl- ions around the duplex show deviation from those of bulk with peaks whose relative intensity depends on the ionic strength. The number of ions monitored as a function of the distance from the duplex reveals a behavior reminiscent of the counterion condensation near the duplex surface. The occurrence of such a phenomenon could affect the Debye length with possible effects on the sensitivity in biosensing experiments.

Nucleosomal DNA unwinding pathway through canonical and non-canonical histone disassembly.

The nucleosome including H2A.B, a mammalian-specific H2A variant, plays pivotal roles in spermatogenesis, embryogenesis, and oncogenesis, indicating unique involvement in transcriptional regulation distinct from canonical H2A nucleosomes. Despite its significance, the exact regulatory mechanism remains elusive. This study utilized solid-state nanopores to investigate DNA unwinding dynamics, applying local force between DNA and histones. Comparative analysis of canonical H2A and H2A.B nucleosomes demonstrated that the H2A.B variant required a lower voltage for complete DNA unwinding. Furthermore, synchronization analysis and molecular dynamics simulations indicate that the H2A.B variant rapidly unwinds DNA, causing the H2A-H2B dimer to dissociate from DNA immediately upon disassembly of the histone octamer. In contrast, canonical H2A nucleosomes unwind DNA at a slower rate, suggesting that the H2A-H2B dimer undergoes a state of stacking at the pore. These findings suggest that nucleosomal DNA in the H2A.B nucleosomes undergoes a DNA unwinding process involving histone octamer disassembly distinct from that of canonical H2A nucleosomes, enabling smoother unwinding. The integrated approach of MD simulations and nanopore measurements is expected to evolve into a versatile tool for studying molecular interactions, not only within nucleosomes but also through the forced dissociation of DNA-protein complexes.

3dRNA/DNA: 3D Structure Prediction from RNA to DNA.

There is an increasing need for determining 3D structures of DNAs, e.g., for increasing the efficiency of DNA aptamer selection. Recently, we have proposed a computational method of 3D structure prediction of DNAs, called 3dDNA, which has been integrated into our original web server 3dRNA, now renamed 3dRNA/DNA (http://biophy.hust.edu.cn/new/3dRNA). Currently, 3dDNA can only output the predicted DNA 3D structures for users but cannot rank them as an energy function for assessing DNA 3D structures is still lacking. Here, we first provide a brief introduction to 3dDNA and then introduce a new energy function, 3dDNAscore, for the assessment of DNA 3D structures. 3dDNAscore is an all-atom knowledge-based potential by integrating 86 atomic types from nucleic acids. Benchmarks demonstrate that 3dDNAscore can effectively identify near-native structures from the decoys generated by 3dDNA, thus enhancing the completeness of 3dDNA.

RiboVision2: A Web Server for Advanced Visualization of Ribosomal RNAs.

RiboVision2 is a web server designed to visualize phylogenetic, structural, and evolutionary properties of ribosomal RNAs simultaneously at the levels of primary, secondary, and three-dimensional structure and in the context of full ribosomal complexes. RiboVision2 instantly computes and displays a broad variety of data; it has no login requirements, is open-source, free for all users, and available at https://ribovision2.chemistry.gatech.edu.

Controlled mechanochemical coupling of anti-junctions in DNA origami arrays.

Allostery is a hallmark of cellular function and important in every biological system. Still, we are only starting to mimic it in the laboratory. Here, we introduce an approach to study aspects of allostery in artificial systems. We use a DNA origami domino array structure which-upon binding of trigger DNA strands-undergoes a stepwise allosteric conformational change. Using two FRET probes placed at specific positions in the DNA origami, we zoom in into single steps of this reaction cascade. Most of the steps are strongly coupled temporally and occur simultaneously. Introduction of activation energy barriers between different intermediate states alters this coupling and induces a time delay. We then apply these approaches to release a cargo DNA strand at a predefined step in the reaction cascade to demonstrate the applicability of this concept in tunable cascades of mechanochemical coupling with both spatial and temporal control.

Cryo-EM: A window into the dynamic world of RNA molecules.

RNAs are critical for complex cellular functions, characterized by their structural versatility and ability to undergo conformational transitions in response to cellular cues. The elusive structures of RNAs are being unraveled with unprecedented clarity, thanks to the technological advancements in structural biology, including nuclear magnetic resonance (NMR), X-ray crystallography, cryo-electron microscopy (cryo-EM) etc. This review focuses on examining the revolutionary impact of cryo-EM on our comprehension of RNA structural dynamics, underscoring the technique's contributions to structural biology and envisioning the future trajectory of this rapidly evolving field.

Spectral-based detection of chromatin loops in multiplexed super-resolution FISH data.

Involved in mitotic condensation, interaction of transcriptional regulatory elements and isolation of structural domains, loop formation has become a paradigm in the deciphering of chromatin architecture and its functional role. Despite the emergence of increasingly powerful genome visualization techniques, the high variability in cell populations and the randomness of conformations still make loop detection a challenge. We introduce an approach for determining the presence and frequency of loops in a collection of experimental conformations obtained by multiplexed super-resolution imaging. Based on a spectral approach, in conjunction with neural networks, this method offers a powerful tool to detect loops in large experimental data sets, both at the population and single-cell levels. The method's performance is confirmed on experimental FISH data where Hi-C and other loop detection results are available. The method is then applied to recently published experimental data, where it provides a detailed and statistically quantified description of the global architecture of the chromosomal region under study.

Unveiling the solution structure of a DNA duplex with continuous silver-modified Watson-Crick base pairs.

The challenge of transforming organized DNA structures into their metallized counterparts persists in the scientific field. In this context, utilizing DNA molecules modified with 7-deazapurine, provides a transformative solution. In this study, we present the solution structure of a DNA duplex that can be transformed into its metallized equivalent while retaining the natural base pairing arrangement through the creation of silver-modified Watson-Crick base pairs. Unlike previously documented X-ray structures, our research demonstrates the feasibility of preserving the intrinsic DNA self-assembly while incorporating AgI into the double helix, illustrating that the binding of silver does not disrupt the canonical base-pairing organization. Moreover, in our case, the uninterrupted AgI chain deviates from forming conventional straight linear chains; instead, it adheres to a helical arrangement dictated by the underlying DNA structure. This research challenges conventional assumptions and opens the door to precisely design structures based on the organization of highly stable Ag-DNA assemblies.