Ionization Patterns and Chemical Reactivity of Cytosine-Guanine Watson-Crick Pairs.

The front cover artwork is provided by Prof. Papadantonakis' group. The image shows a Watson-Crick Guanine-Cytosine pair, and the difference between vertical and adiabatic ionization potentials. Read the full text of the Research Article at 10.1002/cphc.202300946.

The effect of pseudoknot base pairing on cotranscriptional structural switching of the fluoride riboswitch.

A central question in biology is how RNA sequence changes influence dynamic conformational changes during cotranscriptional folding. Here we investigated this question through the study of transcriptional fluoride riboswitches, non-coding RNAs that sense the fluoride anion through the coordinated folding and rearrangement of a pseudoknotted aptamer domain and a downstream intrinsic terminator expression platform. Using a combination of Escherichia coli RNA polymerase in vitro transcription and cellular gene expression assays, we characterized the function of mesophilic and thermophilic fluoride riboswitch variants. We showed that only variants containing the mesophilic pseudoknot function at 37°C. We next systematically varied the pseudoknot sequence and found that a single wobble base pair is critical for function. Characterizing thermophilic variants at 65°C through Thermus aquaticus RNA polymerase in vitro transcription showed the importance of this wobble pair for function even at elevated temperatures. Finally, we performed all-atom molecular dynamics simulations which supported the experimental findings, visualized the RNA structure switching process, and provided insight into the important role of magnesium ions. Together these studies provide deeper insights into the role of riboswitch sequence in influencing folding and function that will be important for understanding of RNA-based gene regulation and for synthetic biology applications.

High frequency of transition to transversion ratio in the stem region of RNA secondary structure of untranslated region of SARS-CoV-2.

Introduction: The propensity of nucleotide bases to form pairs, causes folding and the formation of secondary structure in the RNA. Therefore, purine (R): pyrimidine (Y) base-pairing is vital to maintain uniform lateral dimension in RNA secondary structure. Transversions or base substitutions between R and Y bases, are more detrimental to the stability of RNA secondary structure, than transitions derived from substitutions between A and G or C and T. The study of transversion and transition base substitutions is important to understand evolutionary mechanisms of RNA secondary structure in the 5'  and 3'  untranslated (UTR) regions of SARS-CoV-2. In this work, we carried out comparative analysis of transition and transversion base substitutions in the stem and loop regions of RNA secondary structure of SARS-CoV-2. Methods: We have considered the experimentally determined and well documented stem and loop regions of 5' and 3' UTR regions of SARS-CoV-2 for base substitution analysis. The secondary structure comprising of stem and loop regions were visualized using the RNAfold web server. The GISAID repository was used to extract base sequence alignment of the UTR regions. Python scripts were developed for comparative analysis of transversion and transition frequencies in the stem and the loop regions. Results: The results of base substitution analysis revealed a higher transition (ti) to transversion (tv) ratio (ti/tv) in the stem region of UTR of RNA secondary structure of SARS-CoV-2 reported during the early stage of the pandemic. The higher ti/tv ratio in the stem region suggested the influence of secondary structure in selecting the pattern of base substitutions. This differential pattern of ti/tv values between stem and loop regions was not observed among the Delta and Omicron variants that dominated the later stage of the pandemic. It is noteworthy that the ti/tv values in the stem and loop regions were similar among the later dominant Delta and Omicron variant strains which is to be investigated to understand the rapid evolution and global adaptation of SARS-CoV-2. Conclusion: Our findings implicate the lower frequency of transversions than the transitions in the stem regions of UTRs of SARS-CoV-2. The RNA secondary structures are associated with replication, translation, and packaging, further investigations are needed to understand these base substitutions across different variants of SARS-CoV-2.

From Base Pairs to City Squares: Comprehensive Precision Oncology for the Future.

SUMMARY: An increased understanding of the role of the social determinants of health in cancer prevention, cancer care, and outcomes can lead to their integration into genetics and genomics as well as informing interventions and clinical trials, creating a comprehensive precision oncology framework.

Advancing Genetic Alphabet Expansion: Synthesis of 7-(2-Thienyl)-Imidazo[4,5-b]pyridine (Ds) and 4-(4-Pentyne-1,2-diol)-1-Propynyl-2-Nitropyrrole (Diol-Px) for Use in Replicable Unnatural Base Pairs for PCR Applications.

Expanding the genetic alphabet enhances DNA recombinant technologies by introducing unnatural base pairs (UBPs) beyond the standard A-T and G-C pairs, leading to biomaterials with novel and increased functionalities. Recent developments include UBPs that effectively function as a third base pair in replication, transcription, and/or translation processes. One such UBP, Ds-Px, demonstrates extremely high specificity in replication. Chemically synthesized DNA fragments containing Ds bases are amplified by PCR with the 5'-triphosphates of Ds and Px deoxyribonucleosides (dDsTP and dPxTP). The Ds-Px pair system has applications in enhanced DNA data storage, generation of high-affinity DNA aptamers, and incorporation of functional elements into RNA through transcription. This protocol describes the synthesis of the amidite derivative of Ds (dDs amidite), the triphosphate dDsTP, and the diol-modified dPxTP (Diol-dPxTP) for PCR amplifications involving the Ds-Px pair. © 2024 Wiley Periodicals LLC. Basic Protocol 1: Synthesis of Ds deoxyribonucleoside (dDs) Basic Protocol 2: Synthesis of dDs amidite Basic Protocol 3: Synthesis of dDs triphosphate (dDsTP) Basic Protocol 4: Synthesis of Pn deoxyribonucleoside (4-iodo-dPn) Basic Protocol 5: Synthesis of acetyl-protected diol-modified Px deoxyribonucleoside (Diol-dPx) Basic Protocol 6: Synthesis of Diol-dPx triphosphate (Diol-dPxTP) Basic Protocol 7: Purification of triphosphates Support Protocol 1: Synthesis of Hoffer's chlorosugar Support Protocol 2: Preparation of 0.5 M pyrophosphate in DMF Support Protocol 3: Preparation of 2 M TEAB buffer.

Insights into elastic properties of coarse-grained DNA models: q-stiffness of cgDNA vs cgDNA.

Coarse-grained models have emerged as valuable tools to simulate long DNA molecules while maintaining computational efficiency. These models aim at preserving interactions among coarse-grained variables in a manner that mirrors the underlying atomistic description. We explore here a method for testing coarse-grained vs all-atom models using stiffness matrices in Fourier space (q-stiffnesses), which are particularly suited to probe DNA elasticity at different length scales. We focus on a class of coarse-grained rigid base DNA models known as cgDNA and its most recent version, cgDNA+. Our analysis shows that while cgDNA+ closely follows the q-stiffnesses of the all-atom model, the original cgDNA shows some deviations for twist and bending variables, which are rather strong in the q → 0 (long length scale) limit. The consequence is that while both cgDNA and cgDNA+ give a suitable description of local elastic behavior, the former misses some effects that manifest themselves at longer length scales. In particular, cgDNA performs poorly on twist stiffness, with a value much lower than expected for long DNA molecules. Conversely, the all-atom and cgDNA+ twist are strongly length scale dependent: DNA is torsionally soft at a few base pair distances but becomes more rigid at distances of a few dozen base pairs. Our analysis shows that the bending persistence length in all-atom and cgDNA+ is somewhat overestimated.

Base pair dynamics, electrostatics, and thermodynamics at the LTR-III quadruplex:duplex junction.

G-quadruplexes (GQs) play key regulatory roles within the human genome and have also been identified to play similar roles in other eukaryotes, bacteria, archaea, and viruses. Human immunodeficiency virus 1, the etiological agent of acquired immunodeficiency syndrome, can form two GQs in its long terminal repeat (LTR) promoter region, each of which act to regulate viral gene expression in opposing manners. The major LTR GQ, called LTR-III, is a distinct hybrid GQ containing a 12-nucleotide duplex loop attached to the quadruplex motif. The resulting quadruplex:duplex junction (QDJ) has been hypothesized to serve as a selective drug targeting site. To better understand the dynamics of this QDJ, we performed conventional and enhanced-sampling molecular dynamics simulations using the Drude-2017 force field. We observed unbiased and reversible formation of additional base pairs in the QDJ, between Ade4:Thy14 and Gua3:Thy14. Both base pairs were electrostatically favored, but geometric constraints within the junction may drive the formation of, and preference for, the Ade4:Thy14 base pair. Finally, we demonstrated that the base pairs are separated only by small energy barriers that may enable transitions between both base-paired states. Together, these simulations provide new insights into the dynamics, electrostatics, and thermodynamics of the LTR-III QDJ.

Chemical Repair of Radical Damage to the GC Base Pair by DNA-Bound Bisbenzimidazoles.

The migration of an electron-loss center (hole) in calf thymus DNA to bisbenzimidazole ligands bound in the minor groove is followed by pulse radiolysis combined with time-resolved spectrophotometry. The initially observed absorption spectrum upon oxidation of DNA by the selenite radical is consistent with spin on cytosine (C), as the GC• pair neutral radical, followed by the spectra of oxidized ligands. The rate of oxidation of bound ligands increased with an increase in the ratio (r) ligands per base pair from 0.005 to 0.04. Both the rate of ligand oxidation and the estimated range of hole transfer (up to 30 DNA base pairs) decrease with the decrease in one-electron reduction potential between the GC• pair neutral radical of ca. 1.54 V and that of the ligand radicals (E0', 0.90-0.99 V). Linear plots of log of the rate of hole transfer versus r give a common intercept at r = 0 and a free energy change of 12.2 ± 0.3 kcal mol-1, ascribed to the GC• pair neutral radical undergoing a structural change, which is in competition to the observed hole transfer along DNA. The rate of hole transfer to the ligands at distance, R, from the GC• pair radical, k2, is described by the relationship k2 = k0 exp(constant/R), where k0 includes the rate constant for surmounting a small barrier.

Comprehensive Comparison and Critical Assessment of RNA-Specific Force Fields.

Molecular dynamics simulations play a pivotal role in elucidating the dynamic behaviors of RNA structures, offering a valuable complement to traditional methods such as nuclear magnetic resonance or X-ray. Despite this, the current precision of RNA force fields lags behind that of protein force fields. In this work, we systematically compared the performance of four RNA force fields (ff99bsc0χOL3, AMBERDES, ff99OL3_CMAP1, AMBERMaxEnt) across diverse RNA structures. Our findings highlight significant challenges in maintaining stability, particularly with regard to cross-strand and cross-loop hydrogen bonds. Furthermore, we observed the limitations in accurately describing the conformations of nonhelical structural motif, terminal nucleotides, and also base pairing and base stacking interactions by the tested RNA force fields. The identified deficiencies in existing RNA force fields provide valuable insights for subsequent force field development. Concurrently, these findings offer recommendations for selecting appropriate force fields in RNA simulations.

Expanding the Horizon of the Xeno Nucleic Acid Space: Threose Nucleic Acids with Increased Information Storage.

Xeno nucleic acids (XNAs) constitute a class of synthetic nucleic acid analogues characterized by distinct, non-natural modifications within the tripartite structure of the nucleic acid polymers. While most of the described XNAs contain a modification in only one structural element of the nucleic acid scaffold, this work explores the XNA chemical space to create more divergent variants with modifications in multiple parts of the nucleosidic scaffold. Combining the enhanced nuclease resistance of α-l-threofuranosyl nucleic acid (TNA) and the almost natural-like replication efficiency and fidelity of the unnatural hydrophobic base pair (UBP) TPT3:NaM, novel modified nucleoside triphosphates with a dual modification pattern were synthesized. We investigated the enzymatic incorporation of these nucleotide building blocks by XNA-compatible polymerases and confirmed the successful enzymatic synthesis of TPT3-modified TNA, while the preparation of NaM-modified TNA presented greater challenges. This study marks the first enzymatic synthesis of TNA with an expanded genetic alphabet (exTNA), opening promising opportunities in nucleic acid therapeutics, particularly for the selection and evolution of nuclease-resistant, high-affinity aptamers with increased chemical diversity.

Investigation of the dynamical structures of double-chain deoxyribonucleic acid model in biological sciences.

The present research investigates the double-chain deoxyribonucleic acid model, which is important for the transfer and retention of genetic material in biological domains. This model is composed of two lengthy uniformly elastic filaments, that stand in for a pair of polynucleotide chains of the deoxyribonucleic acid molecule joined by hydrogen bonds among the bottom combination, demonstrating the hydrogen bonds formed within the chain's base pairs. The modified extended Fan sub equation method effectively used to explain the exact travelling wave solutions for the double-chain deoxyribonucleic acid model. Compared to the earlier, now in use methods, the previously described modified extended Fan sub equation method provide more innovative, comprehensive solutions and are relatively straightforward to implement. This method transforms a non-linear partial differential equation into an ODE by using a travelling wave transformation. Additionally, the study yields both single and mixed non-degenerate Jacobi elliptic function type solutions. The complexiton, kink wave, dark or anti-bell, V, anti-Z and singular wave shapes soliton solutions are a few of the creative solutions that have been constructed utilizing modified extended Fan sub equation method that can offer details on the transversal and longitudinal moves inside the DNA helix by freely chosen parameters. Solitons propagate at a consistent rate and retain their original shape. They are widely used in nonlinear models and can be found everywhere in nature. To help in understanding the physical significance of the double-chain deoxyribonucleic acid model, several solutions are shown with graphics in the form of contour, 2D and 3D graphs using computer software Mathematica 13.2. All of the requisite constraint factors that are required for the completed solutions to exist appear to be met. Therefore, our method of strengthening symbolic computations offers a powerful and effective mathematical tool for resolving various moderate nonlinear wave problems. The findings demonstrate the system's potentially very rich precise wave forms with biological significance. The fundamentals of double-chain deoxyribonucleic acid model diffusion and processing are demonstrated by this work, which marks a substantial development in our knowledge of double-chain deoxyribonucleic acid model movements.

A single base pair substitution in zebrafish distinguishes between innate and acute startle behavior regulation.

Behavioral thresholds define the lowest stimulus intensities sufficient to elicit a behavioral response. Establishment of baseline behavioral thresholds during development is critical for proper responses throughout the animal's life. Despite the relevance of such innate thresholds, the molecular mechanisms critical to establishing behavioral thresholds during development are not well understood. The acoustic startle response is a conserved behavior whose threshold is established during development yet is subsequently acutely regulated. We have previously identified a zebrafish mutant line (escapist) that displays a decreased baseline or innate acoustic startle threshold. Here, we identify a single base pair substitution on Chromosome 25 located within the coding sequence of the synaptotagmin 7a (syt7a) gene that is tightly linked to the escapist acoustic hypersensitivity phenotype. By generating animals in which we deleted the syt7a open reading frame, and subsequent complementation testing with the escapist line, we demonstrate that loss of syt7a function is not the cause of the escapist behavioral phenotype. Nonetheless, escapist mutants provide a powerful tool to decipher the overlap between acute and developmental regulation of behavioral thresholds. Extensive behavioral analyses reveal that in escapist mutants the establishment of the innate acoustic startle threshold is impaired, while regulation of its acute threshold remains intact. Moreover, our behavioral analyses reveal a deficit in baseline responses to visual stimuli, but not in the acute regulation of responses to visual stimuli. Together, this work eliminates loss of syt7a as causative for the escapist phenotype and suggests that mechanisms that regulate the establishment of behavioral thresholds in escapist larvae can operate independently from those regulating acute threshold regulation.

Utility of intercalator displacement assays for screening of ligands for i-motif DNA structures.

Cytosine rich sequences can form intercalated, i-motif DNA structures stabilized by hemi-protonated cytosine:cytosine base pairing. These sequences are often located in regulatory regions of genes such as promoters. Ligands targeting i-motif structures may provide potential leads for treatments for genetic disease. The focus on ligands interacting with i-motif DNA has been increasing in recent years. Here, we describe the fluorescent intercalator displacement (FID) assay using thiazole orange binding i-motif DNA and assess the binding affinity of a ligand to the i-motif DNA by displacing thiazole orange. This provides a time and cost-effective high throughput screening of ligands against secondary DNA structures for hit identification.

Potentiometric titrations to study ligand interactions with DNA i-motifs.

i-Motifs are non-canonical secondary structures of DNA formed by mutual intercalation of hemi-protonated cytosine-cytosine base pairs, most typically in slightly acidic conditions (pH<7.0). These structures are well-studied in vitro and have recently been suggested to exist in cells. Despite nearly a decade of active research, the quest for small-molecule ligands that could selectively bind to and stabilize i-motifs continues, and no reference, bona fide i-motif ligand is currently available. This is, at least in part, due to the lack of robust methods to assess the interaction of ligands with i-motifs, since many techniques well-established for studies of other secondary structures (such as CD-, UV-, and FRET-melting) may generate artifacts when applied to i-motifs. Here, we describe an implementation of automated, potentiometric (pH) titrations as a robust isothermal method to assess the impact of ligands or cosolutes on thermodynamic stability of i-motifs. This approach is validated through the use of a cosolute previously known to stabilize i-motifs (PEG2000) and three small-molecule ligands that are able to stabilize, destabilize, or have no effect on the stability of i-motifs, respectively.

A beginner's handbook to identify and characterize i-motif DNA.

Genomic DNA exhibits an innate ability to manifest diverse sequence-dependent secondary structures, serving crucial functions in gene regulation and cellular equilibrium. While extensive research has confirmed the formation of G-quadruplex structures by guanine-rich sequences in vitro and in cells, recent investigations have turned the quadruplex community's attention to the cytosine (C)-rich complementary strands that can adopt unique tetra-stranded conformation, termed as intercalated motif or i-motif. I-motifs are stabilized by hemi-protonated C:CH+ base pairs under acidic conditions. Initially, the in vivo occurrence of i-motifs was underestimated because their formation is favored at non-physiological pH. However, groundbreaking research utilizing the structure-specific iMab antibody and high-throughput sequencing have recently detected their conserved dispersion throughout the genome, challenging previous assumptions. Given the evolving nature of this research field, it becomes imperative to conduct independent in vitro experiments aimed at identifying potential i-motif formation in C-rich sequences and consolidating the findings to address the properties of i-motifs. This chapter serves as an introductory guide for the swift identification of novel i-motifs, where we present an experimental framework for investigating and characterizing i-motif sequences in vitro. In this chapter, we selected a synthetic oligonucleotide (C7T3) sequence and outlined appropriate methodologies for annealing the i-motif structure into suitable buffers. Then, we validated its formation by CD (Circular Dichroism) and NMR (Nuclear Magnetic Resonance) spectroscopy. Finally, we provided a thorough account of the step-by-step procedures to investigate the effect of i-motif formation on the stalling or retardation of DNA replication using high resolution primer extension assays.

Accurate prediction of RNA secondary structure including pseudoknots through solving minimum-cost flow with learned potentials.

Pseudoknots are key structure motifs of RNA and pseudoknotted RNAs play important roles in a variety of biological processes. Here, we present KnotFold, an accurate approach to the prediction of RNA secondary structure including pseudoknots. The key elements of KnotFold include a learned potential function and a minimum-cost flow algorithm to find the secondary structure with the lowest potential. KnotFold learns the potential from the RNAs with known structures using an attention-based neural network, thus avoiding the inaccuracy of hand-crafted energy functions. The specially designed minimum-cost flow algorithm used by KnotFold considers all possible combinations of base pairs and selects from them the optimal combination. The algorithm breaks the restriction of nested base pairs required by the widely used dynamic programming algorithms, thus enabling the identification of pseudoknots. Using 1,009 pseudoknotted RNAs as representatives, we demonstrate the successful application of KnotFold in predicting RNA secondary structures including pseudoknots with accuracy higher than the state-of-the-art approaches. We anticipate that KnotFold, with its superior accuracy, will greatly facilitate the understanding of RNA structures and functionalities.

Investigation of the stability of D5SIC-DNAM-incorporated DNA duplex in Taq polymerase binary system: a systematic classical MD approach.

DNA polymerases are fundamental enzymes that play a crucial role in processing DNA with high fidelity and accuracy ensuring the faithful transmission of genetic information. The recognition of unnatural base pairs (UBPs) by polymerases, enabling their replication, represents a significant and groundbreaking discovery with profound implications for genetic expansion. Romesberg et al. examined the impact of DNA containing 2,6-dimethyl-2H-isoquiniline-1-thione: D5SIC (DS) and 2-methoxy-3-methylnaphthalene: DNAM (DN) UBPs bound to T. aquaticus DNA polymerase (Taq) through crystal structure analysis. Here, we have used polarizable and nonpolarizable classical molecular dynamics (MD) simulations to investigate the structural aspects and stability of Taq in complex with a DNA duplex including a DS-DN pair in the terminal 3' and 5' positions. Our results suggest that the flexibility of UBP-incorporated DNA in the terminal position is arrested by the polymerase, thus preventing fraying and mispairing. Our investigation also reveals that the UBP remains in an intercalated conformation inside the active site, exhibiting two distinct orientations in agreement with experimental findings. Our analysis pinpoints particular residues responsible for favorable interactions with the UBP, with some relying on van der Waals interactions while other on Coulombic forces.

The origin of different bending stiffness between double-stranded RNA and DNA revealed by magnetic tweezers and simulations.

The subtle differences in the chemical structures of double-stranded (ds) RNA and DNA lead to significant variations in their biological roles and medical implications, largely due to their distinct biophysical properties, such as bending stiffness. Although it is well known that A-form dsRNA is stiffer than B-form dsDNA under physiological salt conditions, the underlying cause of this difference remains unclear. In this study, we employ high-precision magnetic-tweezer experiments along with molecular dynamics simulations and reveal that the relative bending stiffness between dsRNA and dsDNA is primarily determined by the structure- and salt-concentration-dependent ion distribution around their helical structures. At near-physiological salt conditions, dsRNA shows a sparser ion distribution surrounding its phosphate groups compared to dsDNA, causing its greater stiffness. However, at very high monovalent salt concentrations, phosphate groups in both dsRNA and dsDNA become fully neutralized by excess ions, resulting in a similar intrinsic bending persistence length of approximately 39 nm. This similarity in intrinsic bending stiffness of dsRNA and dsDNA is coupled to the analogous fluctuations in their total groove widths and further coupled to the similar fluctuation of base-pair inclination, despite their distinct A-form and B-form helical structures.

A 76-base pair duplication within the enhancer region of the HMX1 gene causes sheep microtia.

Sheep congenital microtia is characterized by underdeveloped ears and provides an ideal basis for studying human microtia. This study identified the causal mutation and regulatory mechanisms underlying this disorder. Whole-genome association analysis was conducted using 23 ear tissue samples from sheep with microtia and 28 samples from normal-eared sheep. A significant correlation was found between microtia and a 76-base pair duplication in the enhancer region of the HMX1 gene. Further analysis of offspring phenotypes confirmed an autosomal dominant inheritance pattern. Genotypic analysis showed that individuals that are homozygous for this duplication were earless, heterozygous individuals exhibited shortened ears, and wild-type individuals had normal ears. Moreover, luciferase assays confirmed that this duplication increased HMX1 gene expression, and duplication knock-in mice also exhibited shorter and narrower external ears compared to wild-type mice. Transcriptomic analysis further demonstrated that this duplication enhanced HMX1 gene expression in animal models. This study characterized the causal regulatory mutation underlying sheep microtia.