The Origin of Left-Handed Poly[d(G-C)].

The discovery of a reversible transition in the helical sense of a double-helical DNA was initiated by the first synthesis in 1967 of the alternating sequence poly[d(G-C)]. In 1968, exposure to high salt concentration led to a cooperative isomerization of the double helix manifested by an inversion in the CD spectrum in the 240-310 nm range and in an altered absorption spectrum. The tentative interpretation, reported in 1970 and then in detailed form in a 1972 publication by Pohl and Jovin, was that the conventional right-handed B-DNA structure (R) of poly[d(G-C)] transforms at high salt concentration into a novel, alternative left-handed (L) conformation. The historical course of this development and its aftermath, culminating in the first crystal structure of left-handed Z-DNA in 1979, is described in detail. The research conducted by Pohl and Jovin after 1979 is summarized, ending with an assessment of "unfinished business": condensed Z*-DNA; topoisomerase IIα (TOP2A) as an allosteric ZBP (Z-DNA-binding protein); B-Z transitions of phosphorothioate-modified DNAs; and parallel-stranded poly[d(G-A)], a double helix with high stability under physiological conditions and potentially also left-handed.

Characterization of Z-DNA Using Circular Dichroism.

The B-DNA to Z-DNA transition is a remarkable conformational change in DNA, which was originally observed in poly-GC DNA in the presence of high salt concentration. This eventually prompted the observation of the crystal structure of Z-DNA, a left-handed double-helical DNA, at atomic resolution. Despite advances in Z-DNA research, the application of circular dichroism (CD) spectroscopy as the fundamental technique to characterize this unique DNA conformation has remained constant. In this chapter, we describe a CD spectroscopic method for characterizing the B-DNA to Z-DNA transition of a CG-repeat double-stranded DNA fragment formed from a protein or chemical inducer.

BZ Junctions and Its Application as Probe (2AP) to Detect Z-DNA Formation and Its Effector.

The left-handed Z-DNA is surrounded by right-handed canonical B-DNA, and thus the junction between B- and Z-DNA has been occurred during temporal Z-DNA formation in the genome. The base extrusion structure of the BZ junction may help detect Z-DNA formation in DNAs. Here we describe the BZ junction structural detection by using 2-aminopurine (2AP) fluorescent probe. BZ junction formation can be measured in solution by this method.

Thermogenomic Analysis of Left-Handed Z-DNA Propensities in Genomes.

The initial discovery of left-handed Z-DNA was met with great excitement as a dramatic alternative to the right-handed double-helical conformation of canonical B-DNA. In this chapter, we describe the workings of the program ZHUNT as a computational approach to mapping Z-DNA in genomic sequences using a rigorous thermodynamic model for the transition between the two conformations (the B-Z transition). The discussion starts with a brief summary of the structural properties that differentiate Z- from B-DNA, focusing on those properties that are particularly relevant to the B-Z transition and the junction that splices a left- to right-handed DNA duplex. We then derive the statistical mechanics (SM) analysis of the zipper model that describes the cooperative B-Z transition and show that this analysis very accurately simulates this behavior of naturally occurring sequences that are induced to undergo the B-Z transition through negative supercoiling. A description of the ZHUNT algorithm and its validation are presented, followed by how the program had been applied for genomic and phylogenomic analyses in the past and how a user can access the online version of the program. Finally, we present a new version of ZHUNT (called mZHUNT) that has been parameterized to analyze sequences that contain 5-methylcytosine bases and compare the results of the ZHUNT and mZHUNT analyses on native and methylated yeast chromosome 1.

Synthesis of Aminotroponyl-/Difluoroboronyl Aminotroponyl Deoxyuridine Phosphoramidites.

This report describes the chemical synthesis of aminotroponyl-conjugated deoxyuridine analog (at-dU) and its difluoroboron complex (dfbat-dU) and their phosphoramidites by using the versatile phosphorylating reagent 2-Cyanoethyl N,N-diisopropylchlorophosphoramidite. Tropolone is a non-benzenoid aromatic bioactive natural fluorescent molecule, possessing intramolecular charge transfer and metal chelating properties with transition metal ions such as Cu2+/ Zn2+/ Ni2+ . Its synthetic derivatives, 2-aminotropones also exhibit unique bioactivities and are considered potential therapeutic drug candidate. Recently, the fluorescence properties of aminotropone has improved by complexing with difluoroboron residue that generates aminotroponyl-BODIPY analog. These could be employed for the synthesis of at-dU/dfbat-dU containing DNA oligonucleotides for designing the 11 B/19 F-NMR/fluorescence-based DNA probes. © 2022 Wiley Periodicals LLC. Basic Protocol 1: Synthesis of N-propargyl-2-aminotropone (2) and difluoroboronyl N-propargyl-2-aminotropone (3) molecules. Basic Protocol 2: Synthesis of N-propargyl-2-aminotroponyl deoxyuridinyl (at-dU) phosphoramidites (7). Basic Protocol 3: Synthesis of difluoroboronyl N-propargyl-2-aminotroponyl deoxyuridinyl (dfbat-dU) phosphoramidites (10).

G-Quadruplex DNA and Other Non-Canonical B-Form DNA Motifs Influence Productive and Latent HIV-1 Integration and Reactivation Potential.

The integration of the HIV-1 genome into the host genome is an essential step in the life cycle of the virus and it plays a critical role in the expression, long-term persistence, and reactivation of HIV expression. To better understand the local genomic environment surrounding HIV-1 proviruses, we assessed the influence of non-canonical B-form DNA (non-B DNA) on the HIV-1 integration site selection. We showed that productively and latently infected cells exhibit different integration site biases towards non-B DNA motifs. We identified a correlation between the integration sites of the latent proviruses and non-B DNA features known to potently influence gene expression (e.g., cruciform, guanine-quadruplex (G4), triplex, and Z-DNA). The reactivation potential of latent proviruses with latency reversal agents also correlated with their proximity to specific non-B DNA motifs. The perturbation of G4 structures in vitro using G4 structure-destabilizing or -stabilizing ligands resulted in a significant reduction in integration within 100 base pairs of G4 motifs. The stabilization of G4 structures increased the integration within 300-500 base pairs from G4 motifs, increased integration near transcription start sites, and increased the proportion of latently infected cells. Moreover, we showed that host lens epithelium-derived growth factor (LEDGF)/p75 and cleavage and polyadenylation specificity factor 6 (CPSF6) influenced the distribution of integration sites near several non-B DNA motifs, especially G4 DNA. Our findings identify non-B DNA motifs as important factors that influence productive and latent HIV-1 integration and the reactivation potential of latent proviruses.

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.

Fast Polarizable Water Model for Atomistic Simulations.

Simulating water accurately has been a major challenge in atomistic simulations for decades. Inclusion of electronic polarizability effects holds considerable promise, yet existing approaches suffer from significant computational overheads compared to the widely used nonpolarizable water models. We have developed a globally optimal polarizable water model, OPC3-pol, that explicitly accounts for electronic polarizability with minimal impact on the computational efficiency. OPC3-pol reproduces five key bulk water properties at room temperature with an average relative error of 0.6%. In atomistic simulations, OPC3-pol's computational efficiency is in between that of 3- and 4-point nonpolarizable models; the model supports increased (4 fs) integration time step. OPC3-pol is tested in simulations of globular protein ubiquitin and a B-DNA dodecamer with several AMBER force fields, ff99SB, ff14SB, ff19SB, and OL15, demonstrating structure stability close to reference on multi-microsecond time scale. Simulation of an intrinsically disordered amyloid ß-peptide yields an ensemble with the radius of gyration of a random coil. The proposed water model can be trivially adopted by any package that supports standard nonpolarizable force fields and water models; its intended use is in long classical atomistic simulations where water polarization effects are expected to be important.

Fragment-Based Design of Small Molecules to Study DNA Minor Groove Recognition.

DNA-protein interactions are ubiquitous in cellular processes. Impeding unwanted nucleic acid interactions and protein recognition have therapeutic implications. Therefore, new chemical scaffolds and studies related to their structural basis of nucleic acid recognition are essential for developing high-affinity DNA binders. In this study, we have employed a fragment-based strategy to design and synthesize benzimidazole-guanidinium hybrid compounds and study the individual fragment's role in imparting selectivity and specificity in DNA recognition. The fragments were extensively studied using thermal denaturation, circular dichroism, UV-vis absorption spectroscopy, and molecular docking techniques. The results indicate an interdependent role of the benzimidazole core, polar ends, and the DNA composition in imparting sequence-selective binding of the benzimidazole-guanidinium hybrid compounds in the DNA minor groove. Circular dichroism and molecular docking studies indicated minor groove binding analogous to classical minor groove binders such as DAPI and Hoechst 33258. Thermal denaturation studies indicated that the best binder (compound 8) gave similar thermal stabilization to B-DNA as given by DAPI.

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.

Gene regulation on extrachromosomal DNA.

Oncogene amplification on extrachromosomal DNA (ecDNA) is prevalent in human cancer and is associated with poor outcomes. Clonal, megabase-sized circular ecDNAs in cancer are distinct from nonclonal, small sub-kilobase-sized DNAs that may arise during normal tissue homeostasis. ecDNAs enable profound changes in gene regulation beyond copy-number gains. An emerging principle of ecDNA regulation is the formation of ecDNA hubs: micrometer-sized nuclear structures of numerous copies of ecDNAs tethered by proteins in spatial proximity. ecDNA hubs enable cooperative and intermolecular sharing of DNA regulatory elements for potent and combinatorial gene activation. The 3D context of ecDNA shapes its gene expression potential, selection for clonal heterogeneity among ecDNAs, distribution through cell division, and reintegration into chromosomes. Technologies for studying gene regulation and structure of ecDNA are starting to answer long-held questions on the distinct rules that govern cancer genes beyond chromosomes.

Visualizing "Alternative Isoinformational Engineered" DNA in A- and B-Forms at High Resolution.

A fundamental property of DNA built from four informational nucleotide units (GCAT) is its ability to adopt different helical forms within the context of the Watson-Crick pair. Well-characterized examples include A-, B-, and Z-DNA. For this study, we created an isoinformational biomimetic polymer, built (like standard DNA) from four informational "letters", but with the building blocks being artificial. This ALternative Isoinformational ENgineered (ALIEN) DNA was hypothesized to support two nucleobase pairs, the P:Z pair matching 2-amino-imidazo-[1,2a]-1,3,5-triazin-[8H]-4-one with 6-amino-3-5-nitro-1H-pyridin-2-one and the B:S pair matching 6-amino-4-hydroxy-5-1H-purin-2-one with 3-methyl-6-amino-pyrimidin-2-one. We report two structures of ALIEN DNA duplexes at 1.2 Å resolution and a third at 1.65 Å. All of these are built from a single self-complementary sequence (5'-CTSZZPBSBSZPPBAG) that includes 12 consecutive ALIEN nucleotides. We characterized the helical, nucleobase pair, and dinucleotide step parameters of ALIEN DNA in these structures. In addition to showing that ALIEN pairs retain basic Watson-Crick pairing geometry, two of the ALIEN DNA structures are characterized as A-form DNA and one as B-form DNA. We identified parameters that map differences effecting the transition between the two helical forms; these same parameters distinguish helical forms of isoinformational natural DNA. Collectively, our analyses suggest that ALIEN DNA retains essential structural features of natural DNA, not only its information density and Watson-Crick pairing but also its ability to adopt two canonical forms.

Rapid Molecular Diagnostics in the Field and Laboratory to Detect Plant Pathogen DNA in Potential Insect Vectors.

A variety of sensitive and specific molecular diagnostic assays has been described for detecting nucleic acids in biological samples that may harbor pathogens of interest. These methods include very rapid, isothermal nucleic acid amplification methods that can be deployed outside of the laboratory environment, such as loop-mediated isothermal DNA amplification (LAMP) and recombinase-polymerase amplification (RPA). However, all molecular diagnostic assays must be preceded by nucleic acid extraction from the biological samples of interest, which provides suitable template molecules for the assays. To exploit the features of the amplification assays and be utilized outside of the lab, these methods must be rapid and avoid the need for typical laboratory chemicals and equipment. We describe a protocol for the extraction of DNA from field-collected insects that can be implemented at the point of collection and used to detect the presence of DNA sequences from potential plant pathogens that may be vectored by the insects. This protocol provides template DNA that is suitable for PCR, LAMP, and RPA. The FTA PlantSaver card-based DNA extraction product was also confirmed to amplify the mitochondrial cytochrome oxidase 1 (CO1) universal barcode that could later be sequenced to identify any insect. Lastly, we provide an example using field-collected insects, Neokolla (Graphocephala) heiroglyphica, and demonstrate the detection of the plant pathogen Xylella fastidiosa in carrier insects using PCR, RPA, and LAMP.

Presenting B-DNA as macromolecular crowding agent to improve efficacy of cytochrome c under various stresses.

Existence of numerous biomolecules results in biological fluids to be extremely crowded. Thus, Macromolecular crowding is an essential phenomenon to sustain active conformation of proteins in biological systems. Herein, double helical deoxyribonucleic acid (B-DNA) is presented for the first time as a biomacromolecular crowding system for sustainable packaging of cytochrome c (Cyt C). The peroxidase activity of Cyt C was investigated in the presence of various concentrations of B-DNA (from salmon milt). At an optimized concentration of 0.125 mg/mL B-DNA, an 11-fold higher catalytic activity was found than in native Cyt C with improved stability. Molecular docking and spectroscopic analyses revealed that electrostatic and H-bonding are the main interactions between DNA and Cyt C that affect the structural stability and activity of the protein. Moreover, the catalytic activity and stability of the protein were further investigated in the presence of severe process conditions by UV-visible, circular dichroism, and Fourier-transform infrared spectroscopies. Molecularly crowded Cyt C showed significantly higher activity and stability under severe environments such as high temperature (110 °C), oxidative stress, high pH (pH 10) and biological (trypsin) and chemical denaturants (urea) compared to bare Cyt C. The observed results support the suitability of DNA-based macromolecular crowding media as a viable and effective stabilizer of proteins against multiple stresses.

Understanding the role of non-Watson-Crick base pairs in DNA-protein recognition: Structural and energetic aspects using crystallographic database analysis and quantum chemical calculation.

Specific recognition of DNA base sequences by proteins is vital for life-cycles of all organisms. In a large number of crystal structures of protein-DNA complexes, DNA conformation significantly deviates from the canonical B-DNA structure. A key question is whether such alternate conformations exist prior to protein binding and one is selected for complexation or the structure observed is induced by protein binding. Non-canonical base pairs, such as Hoogsteen base pairs, are often observed in crystal structures of protein-DNA complexes. We decided to explore whether the occurrence of such non-canonical base pairs in protein-DNA complexes is induced by the protein or is selected from pre-existing conformations. Detailed quantum chemical calculations with dispersion-corrected density functional theory (DFT-D) indicated that most of the non-canonical base pairs with DNA bases are stable even in the absence of the interacting amino acids. However, the G:G Hoogsteen base pair, which also appears in the telomere structure, appears to be unstable in the absence of other stabilizing agents, such as positively charged amino acids. Thus, the stability of many of the non-canonical base pair containing duplexes may be close to the canonical B-DNA structure and hence energetically accessible in the ground state; suggesting that the selection from pre-existing conformations may be an important mechanism for observed non-canonical base pairs in protein-DNA complexes.

Extrachromosomal DNA formation enables tumor immune escape potentially through regulating antigen presentation gene expression.

Extrachromosomal DNA (ecDNA) is a type of circular and tumor specific genetic element. EcDNA has been reported to display open chromatin structure, facilitate oncogene amplification and genetic material unequal segregation, and is associated with poor cancer patients' prognosis. The ability of immune evasion is a typical feature for cancer progression, however the tumor intrinsic factors that determine immune evasion remain poorly understood. Here we show that the presence of ecDNA is associated with markers of tumor immune evasion, and obtaining ecDNA could be one of the mechanisms employed by tumor cells to escape immune surveillance. Tumors with ecDNA usually have comparable TMB and neoantigen load, however they have lower immune cell infiltration and lower cytotoxic T cell activity. The microenvironment of tumors with ecDNA shows increased immune-depleted, decreased immune-enriched fibrotic types. Both MHC class I and class II antigen presentation genes' expression are decreased in tumors with ecDNA, and this could be the underlying mechanism for ecDNA associated immune evasion. This study provides evidence that ecDNA formation is an immune escape mechanism for cancer cells.

Decrypting the interaction pattern of Piperlongumine with calf thymus DNA and dodecamer d(CGCGAATTCGCG)2 B-DNA: Biophysical and molecular docking analysis.

The mechanisms of interaction of DNA with pharmacological molecules are critical to understanding their therapeutic actions on physiological systems. Piperlongumine is widely studied for its anticancer potential. Multi-spectrometry, calorimetry and in silico studies were employed to study the interaction of piperlongumine and calf thymus DNA. UV-Vis spectroscopy illustrated a hyperchromic pattern in spectra of the calf thymus DNA-piperlongumine complex, while fluorescent quenching was observed in emission spectral studies. Competitive displacement assay demonstrated higher displacement and binding constant for DNA-rhodamine B complex by piperlongumine than DNA-methylene blue complex. Differential scanning calorimetry presented non-significant changes in melting temperature and molecular docking presented the precise interaction site of piperlongumine with calf thymus DNA at minor groove. Further, piperlongumine treatment did not result in pBluescript KS plasmid DNA cleavage as revealed from the DNA topology assay. All these experiments confirmed the binding of piperlongumine with DNA through minor groove binding mode.

B-DNA Structure and Stability: The Role of Nucleotide Composition and Order.

Invited for this month's cover are the groups of Célia Fonseca Guerra at the Vrije Universiteit Amsterdam and Leiden University, Giampaolo Barone from the Università degli Studi di Palermo, and F. Matthias Bickelhaupt at Vrije Universiteit Amsterdam and Radboud University Nijmegen. The cover picture shows the four primary interaction components (hydrogen bonding, cross-terms, base stacking, and solvation) that determine the stability of B-DNA duplexes. Quantum chemical analyses identify an interplay between the stabilizing hydrogen bonds between nucleotides that drive the formation of the DNA double-strand, and the destabilizing loss of stacking interactions within individual strands combined with partial desolvation. The sequence-dependence in the duplex stability originates mainly from the cross-terms, which can be attractive or repulsive. Read the full text of their Research Article at 10.1002/open.202100231.

B-DNA Structure and Stability: The Role of Nucleotide Composition and Order.

We have quantum chemically analyzed the influence of nucleotide composition and sequence (that is, order) on the stability of double-stranded B-DNA triplets in aqueous solution. To this end, we have investigated the structure and bonding of all 32 possible DNA duplexes with Watson-Crick base pairing, using dispersion-corrected DFT at the BLYP-D3(BJ)/TZ2P level and COSMO for simulating aqueous solvation. We find enhanced stabilities for duplexes possessing a higher GC base pair content. Our activation strain analyses unexpectedly identify the loss of stacking interactions within individual strands as a destabilizing factor in the duplex formation, in addition to the better-known effects of partial desolvation. Furthermore, we show that the sequence-dependent differences in the interaction energy for duplexes of the same overall base pair composition result from the so-called "diagonal interactions" or "cross terms". Whether cross terms are stabilizing or destabilizing depends on the nature of the electrostatic interaction between polar functional groups in the pertinent nucleobases.

Assessing the DNA structural integrity via selective annihilation of Watson-Crick hydrogen bonds: Insights from molecular dynamics simulations.

Understanding the role of base pairing and stacking displayed by polynucleotide chains interwind together resulting in a double-helical B-DNA type structure is crucial to gaining access to the sophisticated structural arrangement of DNA. Several computational and experimental studies hinted towards the dominance of base pairing over stacking for duplex stability. To find out how significant the individual Watson-Crick hydrogen bonds are in maintaining the double-helical integrity of the DNA, in the present article, we selectively switched off the hydrogen bonds (one specific bond or their combinations in all the base pairs at a time) via manipulating the force fields for A-T and G-C base pairs. We studied 12 systems in total via all-atom explicit-solvent molecular dynamics simulations (200 ns each). The MD output structures were compared with the control system by means of structural, dynamic, and energetic properties to monitor the overall consequences of removing H-bond(s) on the B-DNA characteristics of the model systems. Our findings suggest that all the individual hydrogen bonds involved in base pairing are vital for maintaining the DNA structural integrity as any possible alteration in Watson-Crick hydrogen bond(s) leads to the disintegration/collapse of DNA strands resulting in unfolded states.