Theoretical Insights into N-Glycoside Bond Cleavage of 5-Carboxycytosine by Thymine DNA Glycosylase: A QM/MM Study.

Thymine DNA glycosylase (TDG)-mediated excision of 5-formylcytosine and 5-carboxylcytosine (5-caC) is a critical step in active DNA demethylation. Herein, we employed a combined quantum mechanics/molecular mechanics approach to investigate the reaction mechanism of TDG-catalyzed N-glycosidic bond cleavage of 5-caC. The calculated results show that TDG-catalyzed 5-caC excision follows a concerted (SN2) mechanism in which glycosidic bond dissociation is coupled with nucleophile attack. Protonation of the 5-caC anion contributes to the cleavage of the N-glycoside bond, in which the N3-protonated zwitterion and imino tautomers are more favorable than carboxyl-protonated amino tautomers. This is consistent with the experimental data. Furthermore, our results reveal that the configuration rearrangement process of the protonated 5-caC would lower the stability of the N-glycoside bond and substantially reduce the barrier height for the subsequent C1'-N1 bond cleavage. This should be attributed to the smaller electrostatic repulsion between the leaving base and the negative phosphate group as a result of the structural rearrangement.

A DFTB study on the electronic response of encapsulated DNA nucleobases onto chiral CNTs as a sequencer.

Sequencing the DNA nucleobases is essential in the diagnosis and treatment of many diseases related to human genes. In this article, the encapsulation of DNA nucleobases with some of the important synthesized chiral (7, 6), (8, 6), and (10, 8) carbon nanotubes were investigated. The structures were modeled by applying density functional theory based on tight binding method (DFTB) by considering semi-empirical basis sets. Encapsulating DNA nucleobases on the inside of CNTs caused changes in the electronic properties of the selected chiral CNTs. The results confirmed that van der Waals (vdW) interactions, π-orbitals interactions, non-bonded electron pairs, and the presence of high electronegative atoms are the key factors for these changes. The result of electronic parameters showed that among the CNTs, CNT (8, 6) is a suitable choice in sequencing guanine (G) and cytosine (C) DNA nucleobases. However, they are not able to sequence adenine (A) and thymine (T). According to the band gap energy engineering approach and absorption energy, the presence of G and C DNA nucleobases decreased the band gap energy of CNTs. Hence selected CNTs suggested as biosensor substrates for sequencing G and C DNA nucleobases.

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.

Computational analysis of antiviral drugs using topological descriptors.

Many health challenges are attributed to viral infections, which represent significant concerns in public health. Among these infections, diseases such as herpes simplex virus (HSV), cytomegalovirus (CMV), and varicella-zoster virus (VZV) infections have garnered attention due to their prevalence and impact on human health. There are specific antiviral medications available for the treatment of these viral infections. Drugs like Cidofovir, Valacyclovir, and Acyclovir are commonly prescribed. These antiviral drugs are known for their efficacy against herpesviruses and related viral infections, leveraging their ability to inhibit viral DNA polymerase. A molecular descriptor is a numerical value that correlates with specific physicochemical properties of a molecular graph. This article explores the calculation of distance-based topological descriptors, including the Trinajstic, Mostar, Szeged, and PI descriptors for the aforementioned antiviral drugs. These descriptors provide insights into these drugs' structural and physicochemical characteristics, aiding in understanding their mechanism of action and the development of new therapeutic agents.

The bisulfite reaction with cytosine and genomic DNA structure.

The bisulfite reaction with native DNA has been extensively employed in the detection of non-B DNA structures that can form spontaneously in DNA. These sequences are dynamic in that they can adopt both normal Watson-Crick paired B-DNA or unusual structures like the Triplex, G-Quadruplex, i-motif and Cruciform or Hairpin. Considerable evidence now suggests that these dynamic sequences play roles in both epigenetics and mutagenesis. The bisulfite reaction with native DNA offers a key approach to their detection. In this application whole cells, isolated nuclei or isolated DNA are treated with bisulfite under non-denaturing conditions in order to detect bisulfite accessible regions DNA that are associated with these structures. Here I review the stereochemistry of the bisulfite reaction, the electronic structure of its DNA cytosine substrates and its application in the detection of unusual structures in native DNA.

Excited States in Single-Stranded and i-Motif DNA with Silver Ions.

We have studied the excited states and structural properties for the complexes of cytosine (dC)10 chains with silver ions (Ag+) in a wide range of the Ag+ to DNA ratio (r) and pH conditions using circular dichroism, steady-state absorption, and fluorescence spectroscopy along with the ultrafast fluorescence upconversion technique. We also calculated vertical electronic transition energies and determined the nature of the corresponding excited states in some models of the cytosine-Ag+ complexes. We show that (dC)10 chains in the presence of silver ions form a duplex stabilized by C-Ag+-C bonds. It is also shown that the i-motif structure formed by (dC)10 chains is destabilized in the presence of Ag+ ions. The excited-state properties in the studied complexes depend on the amount of binding ions and the binding sites, which is supported by the calculations. In particular, new low-lying excited states appear when the second Ag+ ion interacts with the O atom of cytosine in the C-Ag+-C pairs. A similar picture is observed in the case when one Ag+ ion interacts with one cytosine via the N7 atom.

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.

Plasmon-emitter coupling in cytosine-rich hairpin DNA-templated silver nanoclusters: Thermal reversibility, white light emission, and dynamics inside live cells.

Bio-templated luminescent noble metal nanoclusters (NCs) have attracted great attention for their intriguing physicochemical properties. Continuous efforts are being made to prepare NCs with high fluorescence quantum yield (QY), good biocompatibility, and tunable emission properties for their widespread practical applications as new-generation environment-friendly photoluminescent materials in materials chemistry and biological systems. Herein, we explored the unique photophysical properties of silver nanoclusters (AgNCs) templated by cytosine-rich customized hairpin DNA. Our results indicate that a 36-nucleotide containing hairpin DNA with 20 cytosine (C20) in the loop can encapsulate photostable red-emitting AgNCs with an absolute QY of ∼24%. The luminescent properties in these DNA-templated AgNCs were found to be linked to the coupling between the surface plasmon and the emitter. These AgNCs exhibited excellent thermal sensitivity and were employed to produce high-quality white light emission with an impressive color rendering index of 90 in the presence of dansyl chloride. In addition, the as-prepared luminescent AgNCs possessing excellent biocompatibility can effectively mark the nuclear region of HeLa cells and can be employed as a luminescent probe to monitor the cellular dynamics at a single molecular resolution.

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.

5-Formylcytosine mediated DNA-peptide cross-link induces predominantly semi-targeted mutations in both Escherichia coli and human cells.

Histone proteins can become trapped on DNA in the presence of 5-formylcytosine (5fC) to form toxic DNA-protein conjugates. Their repair may involve proteolytic digestion resulting in DNA-peptide cross-links (DpCs). Here, we have investigated replication of a model DpC comprised of an 11-mer peptide (NH2-GGGKGLGK∗GGA) containing an oxy-lysine residue (K∗) conjugated to 5fC in DNA. Both CXG and CXT (where X = 5fC-DpC) sequence contexts were examined. Replication of both constructs gave low viability (<10%) in Escherichia coli, whereas TLS efficiency was high (72%) in HEK 293T cells. In E. coli, the DpC was bypassed largely error-free, inducing only 2 to 3% mutations, which increased to 4 to 5% with SOS. For both sequences, semi-targeted mutations were dominant, and for CXG, the predominant mutations were G→T and G→C at the 3'-base to the 5fC-DpC. In HEK 293T cells, 7 to 9% mutations occurred, and the dominant mutations were the semi-targeted G → T for CXG and T → G for CXT. These mutations were reduced drastically in cells deficient in hPol η, hPol ι or hPol ζ, suggesting a role of these TLS polymerases in mutagenic TLS. Steady-state kinetics studies using hPol η confirmed that this polymerase induces G → T and T → G transversions at the base immediately 3' to the DpC. This study reveals a unique replication pattern of 5fC-conjugated DpCs, which are bypassed largely error-free in both E. coli and human cells and induce mostly semi-targeted mutations at the 3' position to the lesion.

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

Gas-phase and aqueous vertical ionization potentials, vIPgas and vIPaq respectively and measurements of the molecular electrostatic and local ionization maps calculated at the DFT/B3LYP-D3/ 6-311+G** level of theory and the C-PCM reaction field model for single- and double-stranded CpG and 5MeCpG pairs show that the vIPaq for single- and double-stranded pairs of C-G and 5MeC-G are practically the same, in the range of 5.79 to 5.81 eV. The aqueous adiabatic ionization potentials for single-stranded CpG and 5MeCpG are 5.52 eV and 5.51 eV respectively and they reflect the nuclear reorganization that takes place after the abstraction of the electron. The aqueous adiabatic ionization energy values that correspond to the CpG+. radical cation and the hydrated electron, e-,, being at infinite distance, adIPaq+Vo, are 3.92 eV and 3.91 eV respectively with (Vo=-1.6 eV) Analysis of data suggest that the HOMO-LUMO energy gap in the hard/soft-acid/base (HSAB) concept cannot be used a priori to determine the effect of cytosine methylation on the guanine enhanced oxidative damage in DNA.

Prediction of DNA i-motifs via machine learning.

i-Motifs (iMs), are secondary structures formed in cytosine-rich DNA sequences and are involved in multiple functions in the genome. Although putative iM forming sequences are widely distributed in the human genome, the folding status and strength of putative iMs vary dramatically. Much previous research on iM has focused on assessing the iM folding properties using biophysical experiments. However, there are no dedicated computational tools for predicting the folding status and strength of iM structures. Here, we introduce a machine learning pipeline, iM-Seeker, to predict both folding status and structural stability of DNA iMs. The programme iM-Seeker incorporates a Balanced Random Forest classifier trained on genome-wide iMab antibody-based CUT&Tag sequencing data to predict the folding status and an Extreme Gradient Boosting regressor to estimate the folding strength according to both literature biophysical data and our in-house biophysical experiments. iM-Seeker predicts DNA iM folding status with a classification accuracy of 81% and estimates the folding strength with coefficient of determination (R2) of 0.642 on the test set. Model interpretation confirms that the nucleotide composition of the C-rich sequence significantly affects iM stability, with a positive correlation with sequences containing cytosine and thymine and a negative correlation with guanine and adenine.

Evolved Readers of 5-Carboxylcytosine CpG Dyads Reveal a High Versatility of the Methyl-CpG-Binding Domain for Recognition of Noncanonical Epigenetic Marks.

Mammalian genomes are regulated by epigenetic cytosine (C) modifications in palindromic CpG dyads. Including canonical cytosine 5-methylation (mC), a total of four different 5-modifications can theoretically co-exist in the two strands of a CpG, giving rise to a complex array of combinatorial marks with unique regulatory potentials. While tailored readers for individual marks could serve as versatile tools to study their functions, it has been unclear whether a natural protein scaffold would allow selective recognition of marks that vastly differ from canonical, symmetrically methylated CpGs. We conduct directed evolution experiments to generate readers of 5-carboxylcytosine (caC) dyads based on the methyl-CpG-binding domain (MBD), the widely conserved natural reader of mC. Despite the stark steric and chemical differences to mC, we discover highly selective, low nanomolar binders of symmetric and asymmetric caC-dyads. Together with mutational and modelling studies, our findings reveal a striking evolutionary flexibility of the MBD scaffold, allowing it to completely abandon its conserved mC recognition mode in favour of noncanonical dyad recognition, highlighting its potential for epigenetic reader design.

Selective Formation of Pd-DNA Hybrids Using Tailored Palladium-Mediated Base Pairs: Towards Heteroleptic Pd-DNA Systems.

The formation of highly organized metal-DNA structures has significant implications in bioinorganic chemistry, molecular biology and material science due to their unique properties and potential applications. In this study, we report on the conversion of single-stranded polydeoxycytidine (dC15 ) into a Pd-DNA supramolecular structure using the [Pd(Aqa)] complex (Aqa=8-amino-4-hydroxyquinoline-2-carboxylic acid) through a self-assembly process. The resulting Pd-DNA assembly closely resembles a natural double helix, with continuous [Pd(Aqa)(C)] (C=cytosine) units serving as palladium-mediated base pairs, forming interbase hydrogen bonds and intrastrand stacking interactions. Notably, the design of the [Pd(Aqa)] complex favours the interaction with cytosine, distinguishing it from our previously reported [Pd(Cheld)] complex (Cheld=chelidamic acid). This finding opens possibilities for creating heteroleptic Pd-DNA hybrids where different complexes specifically bind to nucleobases. We confirmed the Pd-DNA supramolecular structural assembly and selective binding of the complexes using NMR spectroscopy, circular dichroism, mass spectrometry, isothermal titration calorimetry, and DFT calculations.

Intermolecular hydrogen bonding in DNA base pairs interacting with different numbers of bare and hydrated Li+: NBO, QTAIM, and computational spectroscopic studies.

In this study, the effect of different numbers of Li+ interacting with various sites of DNA base pairs (adenine-thymine (AT) and cytosine-guanine (GC)) on the base pair structures, the strength of hydrogen bonding between the bases, and spectroscopic properties (IR and absorption spectra) of the base pairs was investigated. Two quantum computational analyses, the natural bonding orbitals (NBO) and the quantum theory of atoms in molecules (QTAIM), were used to follow the change in the strength of hydrogen bonds between the bases in each pair. The type of base pair's site interacting with Li+ showed different effects on the change in the strength of the hydrogen bonds between the bases. The IR and absorption spectra of the lithiated base pairs were calculated and compared with those of bare base pairs. This comparison provided the changes in the spectra as a fingerprint for the structural identification of different lithiated base pairs. Also, the determination of the change in the strength of hydrogen bonds in the lithiated base pairs compared to their bare base pairs. In the other part of this study, the effect of the hydration of the attached Li+ in the structure of lithiated base pairs on the strength of their hydrogen bonds and spectra was investigated.

Site specifically probing the unfolding process of human telomere i-motif DNA using vibrationally enhanced alkynyl stretch.

The microscopic unfolding process of a cytosine-rich DNA forming i-motif by hemi-protonated base pairs is related to gene regulation. However, the detailed thermal unfolding mechanism and the protonation/deprotonation status of site-specific cytosine in DNA in a physiological environment are still obscure. To address this issue, a vibration-enhanced CC probe tagged on 5'E terminal cytosine of human telomere i-motif DNA was examined using linear and nonlinear infrared (IR) spectroscopies and quantum-chemistry calculations. The CC probe extended into the major groove of the i-motif was found using nonlinear IR results only to introduce a minor steric effect on both steady-state structure and local structure dynamics; however, its IR absorption profile effectively reports the cleavage of the hemi-protonated base pair of C1-C13 upon the unfolding with C1 remaining protonated. The temperature mid-point (Tm) of the local transition reported using the CC tag was slightly lower than the Tm of global transition, and the enthalpy of the former exceeds 60% of the global transition. It is shown that the base-pair unraveling is noncooperative, with outer base pairs breaking first and being likely the rate limiting step. Our results offered an in-depth understanding of the macroscopic unfolding characteristics of the i-motif DNA and provided a nonlinear IR approach to monitoring the local structural transition and dynamics of DNA and its complexes.

Triplet Excimer Formation in a DNA Duplex with Silver Ion-Mediated Base Pairs.

The dynamics of excited electronic states in self-assembled structures formed between silver(I) ions and cytosine-containing DNA strands or monomeric cytosine derivatives were investigated by time-resolved infrared (TRIR) spectroscopy and quantum mechanical calculations. The steady-state and time-resolved spectra depend sensitively on the underlying structures, which change with pH and the nucleobase and silver ion concentrations. At pH ∼ 4 and low dC20 strand concentration, an intramolecularly folded i-motif is observed, in which protons, and not silver ions, mediate C-C base pairing. However, at the higher strand concentrations used in the TRIR measurements, dC20 strands associate pairwise to yield duplex structures containing C-Ag+-C base pairs with a high degree of propeller twisting. UV excitation of the silver ion-mediated duplex produces a long-lived excited state, which we assign to a triplet excimer state localized on a pair of stacked cytosines. The computational results indicate that the propeller-twisted motifs induced by metal-ion binding are responsible for the enhanced intersystem crossing that populates the triplet state and not a generic heavy atom effect. Although triplet excimer states have been discussed frequently as intermediates in the formation of cyclobutane pyrimidine dimers, we find neither computational nor experimental evidence for cytosine-cytosine photoproduct formation in the systems studied. These findings provide a rare demonstration of a long-lived triplet excited state that is formed in a significant yield in a DNA duplex, demonstrating that supramolecular structural changes induced by metal ion binding profoundly affect DNA photophysics.

Electron interaction with DNA constituents in aqueous phase.

Electron driven chemistry of biomolecules in aqueous phase presents the realistic picture to study molecular processes. In this study we have investigated the interactions of electrons with the DNA constituents in their aqueous phase in order to obtain the quantities useful for DNA damage assessment. We have computed the inelastic mean free path (IMFP), mass stopping power (MSP) and absorbed dose (D) for the DNA constituents (Adenine, Cytosine, Guanine, Thymine and Uracil) in the aqueous medium from ionisation threshold to 5000 eV. We have modified complex optical potential formalism to include band gap of the systems to calculate inelastic cross sections which are used to estimate these entities. This is the maiden attempt to report these important quantities for the aqueous DNA constituents. We have compared our results with available data in gas and other phase and have observed explicable accord for IMFP and MSP. Since these are the first results of absorbed dose (D) for these compounds, we have explored present results vis-a-vis dose absorption in water.

Modular cytosine base editing promotes epigenomic and genomic modifications.

Prokaryotic and eukaryotic adaptive immunity differ considerably. Yet, their fundamental mechanisms of gene editing via Cas9 and activation-induced deaminase (AID), respectively, can be conveniently complimentary. Cas9 is an RNA targeted dual nuclease expressed in several bacterial species. AID is a cytosine deaminase expressed in germinal centre B cells to mediate genomic antibody diversification. AID can also mediate epigenomic reprogramming via active DNA demethylation. It is known that sequence motifs, nucleic acid structures, and associated co-factors affect AID activity. But despite repeated attempts, deciphering AID's intrinsic catalytic activities and harnessing its targeted recruitment to DNA is still intractable. Even recent cytosine base editors are unable to fully recapitulate AID's genomic and epigenomic editing properties. Here, we describe the first instance of a modular AID-based editor that recapitulates the full spectrum of genomic and epigenomic editing activity. Our 'Swiss army knife' toolbox will help better understand AID biology per se as well as improve targeted genomic and epigenomic editing.