Twisting right to left: A…A mismatch in a CAG trinucleotide repeat overexpansion provokes left-handed Z-DNA conformation

Z-DNA at the crossroads: untangling its role in genome dynamics

DNA can fold into noncanonical left-handed Z-DNA conformation beyond the right-handed B-DNA. While its crystal structure was discovered nearly four decades ago, it was predominantly considered a structural curiosity. Recent evidence suggests that Z-DNA formation occurs in nuclear and mitochondrial DNA (mtDNA), with significant biological implications. However, our understanding of its roles remains in its infancy, primarily due to a lack of study tools. In this review we summarize the structure and function of Z-DNA within the genome while addressing the difficulties associated with identifying and investigating its role(s). We then critically evaluate several intracellular factors that can modulate and regulate Z-DNA. Additionally, we discuss the recent technological and methodological advances that may overcome the challenges and enhance our understanding of Z-DNA.

RNA Editing by ADAR Adenosine Deaminases in the Cell Models of CAG Repeat Expansion Diseases: Significant Effect of Differentiation from Stem Cells into Brain Organoids in the Absence of Substantial Influence of CAG Repeats on the Level of Editing

Expansion of CAG repeats in certain genes is a known cause of several neurodegenerative diseases, but exact mechanism behind this is not yet fully understood. It is believed that the double-stranded RNA regions formed by CAG repeats could be harmful to the cell. This study aimed to test the hypothesis that these RNA regions might potentially interfere with ADAR RNA editing enzymes, leading to the reduced A-to-I editing of RNA and activation of the interferon response. We studied induced pluripotent stem cells (iPSCs) derived from the patients with Huntington's disease or ataxia type 17, as well as midbrain organoids developed from these cells. A targeted panel for next-generation sequencing was used to assess editing in the specific RNA regions. Differentiation of iPSCs into brain organoids led to increase in the ADAR2 gene expression and decrease in the expression of protein inhibitors of RNA editing. As a result, there was increase in the editing of specific ADAR2 substrates, which allowed identification of differential substrates of ADAR isoforms. However, comparison of the pathology and control groups did not show differences in the editing levels among the iPSCs. Additionally, brain organoids with 42-46 CAG repeats did not exhibit global changes. On the other hand, brain organoids with the highest number of CAG repeats in the huntingtin gene (76) showed significant decrease in the level of RNA editing of specific transcripts, potentially involving ADAR1. Notably, editing of the long non-coding RNA PWAR5 was nearly absent in this sample. It could be stated in conclusion that in most cultures with repeat expansion, the hypothesized effect on RNA editing was not confirmed.

Non-canonical DNA structures: Diversity and disease association

A complete understanding of DNA double-helical structure discovered by James Watson and Francis Crick in 1953, unveil the importance and significance of DNA. For the last seven decades, this has been a leading light in the course of the development of modern biology and biomedical science. Apart from the predominant B-form, experimental shreds of evidence have revealed the existence of a sequence-dependent structural diversity, unusual non-canonical structures like hairpin, cruciform, Z-DNA, multistranded structures such as DNA triplex, G-quadruplex, i-motif forms, etc. The diversity in the DNA structure depends on various factors such as base sequence, ions, superhelical stress, and ligands. In response to these various factors, the polymorphism of DNA regulates various genes via different processes like replication, transcription, translation, and recombination. However, altered levels of gene expression are associated with many human genetic diseases including neurological disorders and cancer. These non-B-DNA structures are expected to play a key role in determining genetic stability, DNA damage and repair etc. The present review is a modest attempt to summarize the available literature, illustrating the occurrence of non-canonical structures at the molecular level in response to the environment and interaction with ligands and proteins. This would provide an insight to understand the biological functions of these unusual DNA structures and their recognition as potential therapeutic targets for diverse genetic diseases.

Biophysical interaction between lanthanum chloride and (CG)n or (GC)n repeats: A reversible B-to-Z DNA transition

The transition from right-handed to left-handed DNA is not only acts as the controlling factor for switching gene expression but also has equal importance in designing nanomechanical devices. The (CG)_n and (GC)_n repeat sequences are well known model molecules to study B-Z transition in the presence of higher concentration of monovalent cations. In this communication, we report a cyclic transition in (CG)₆ DNA using millimolar concentration of trivalent lanthanide salt LaCl₃. The controlled and reversible transition was seen in (CG)₁₂, and (GC)₁₂ DNA employing CD spectroscopy. While LaCl₃ failed to induce B-Z transition in shorter oligonucleotides such as (CG)₃ and (GC)₃, a smooth B-Z transition was recorded for (CG)₆, (CG)₁₂ and (GC)₁₂ sequences. Interestingly, the phenomenon was reversible (Z-B transition) with addition of EDTA. Particularly, two rounds of cyclic transition (B-Z-B-Z-B) have been noticed in (CG)₆ DNA in presence of LaCl₃ and EDTA which strongly suggest that B-Z transition is reversible in short repeat sequences. Thermal melting and annealing behaviour of B-DNA are reversible while the thermal melting of LaCl₃-induced Z-DNA is irreversible which suggest a stronger binding of LaCl₃ to the phosphate backbone of Z-DNA. This was further supported by isothermal titration calorimetric study. Molecular dynamics (MD) simulation indicates that the mode of binding of La³⁺ (of LaCl₃) with d(CG)₈.d(CG)₈ is through the minor groove, wherein, 3 out of 11 La³⁺ bridge the anionic oxygens of the complementary strands. Such a tight coordination of La³⁺ with the anionic oxygens at the minor groove surface may be the reason for the experimentally observed irreversibility of LaCl₃-induced Z-DNA seen in longer DNA fragments. Thus, these results indicate LaCl₃ can easily be adopted as an inducer of left-handed DNA in other short oligonucleotides sequences to facilitate the understanding of the molecular mechanism of B-Z transition.

Construction of ssDNA-Attached LR-Chimera Involving Z-DNA for ZBP1 Binding Analysis

The binding of proteins to Z-DNA is hard to analyze, especially for short non-modified DNA, because it is easily transferred to B-DNA. Here, by the hybridization of a larger circular single-stranded DNA (ssDNA) with a smaller one, an LR-chimera (involving a left-handed part and a right-handed one) with an ssDNA loop is produced. The circular ssDNAs are prepared by the hybridization of two ssDNA fragments to form two nicks, followed by nick sealing with T4 DNA ligase. No splint (a scaffold DNA for circularizing ssDNA) is required, and no polymeric byproducts are produced. The ssDNA loop on the LR-chimera can be used to attach it with other molecules by hybridization with another ssDNA. The gel shift binding assay with Z-DNA specific binding antibody (Z22) or Z-DNA binding protein 1 (ZBP1) shows that stable Z-DNA can form under physiological ionic conditions even when the extra ssDNA part is present. Concretely, a 5'-terminal biotin-modified DNA oligonucleotide complementary to the ssDNA loop on the LR-chimera is used to attach it on the surface of a biosensor inlaid with streptavidin molecules, and the binding constant of ZBP1 with Z-DNA is analyzed by BLI (bio-layer interferometry). This approach is convenient for quantitatively analyzing the binding dynamics of Z-DNA with other molecules.

Nonalternating purine pyrimidine sequences can form stable left-handed DNA duplex by strong topological constraint

In vivo, left-handed DNA duplex (usually refers to Z-DNA) is mainly formed in the region of DNA with alternating purine pyrimidine (APP) sequence and plays significant biological roles. It is well known that d(CG)_n sequence can form Z-DNA most easily under negative supercoil conditions, but its essence has not been well clarified. The study on sequence dependence of Z-DNA stability is very difficult without modification or inducers. Here, by the strong topological constraint caused by hybridization of two complementary short circular ssDNAs, left-handed duplex part was generated for various sequences, and their characteristics were investigated by using gel-shift after binding to specific proteins, CD and T_m analysis, and restriction enzyme cleavage. Under the strong topological constraint, non-APP sequences can also form left-handed DNA duplex as stable as that of APP sequences. As compared with non-APP sequences, the thermal stability difference for APP sequences between Z-form and B-form is smaller, which may be the reason that Z-DNA forms preferentially for APP ones. This result can help us to understand why nature selected APP sequences to regulate gene expression by transient Z-DNA formation, as well as why polymer with chirality can usually form both duplexes with left- or right-handed helix.

Conformational distortions induced by periodically recurring A…A in d(CAG).d(CAG) provide stereochemical rationale for the trapping of MSH2.MSH3 in polyQ disorders

CAG repeat instability causes a number of neurodegenerative disorders. The unusual hairpin stem structure formed by the CAG repeats in DNA traps the human mismatch repair MSH2.MSH3 (Mutsβ) complex. To understand the mechanism behind the abnormal binding of Mutsβ with the imperfect hairpin stem structure formed by CAG repeats, molecular dynamics simulations have been carried out for Mutsβ-d(CAG)₂(CAG)(CAG)₂.d(CTG)₂(CAG)(CTG)₂ (1 A…A mismatch) and Mutsβ-d(CAG)₅.d(CAG)₅ (5 mismatches, wherein, A…A occurs periodically) complexes. The interaction of MSH3 residue Tyr₂₄₅ at the minor groove side of A…A, an essential interaction responsible for the recognition by Mutsβ, are retained in both the cases. Nevertheless, the periodic unwinding caused by the nonisostericity of A…A with the flanking canonical base pairs in d(CAG)₅.d(CAG)₅ distorts the regular B-form geometry. Such an unwinding exposes one of the A…A mismatches (that interacts with Tyr₂₄₅) at the major groove side and also facilitates the on and off hydrogen bonding interaction with Lys₅₄₆ sidechain (MSH2-domain-IV). In contrast, kinking of the DNA towards the major groove in Mutsβ-d(CAG)₂(CAG)(CAG)₂.d(CTG)₂(CAG)(CTG)₂ doesn’t facilitate such an exposure of the bases at the major groove. Further, the unwinding of the helix in d(CAG)₅.d(CAG)₅ enhances the tighter binding between MSH2-domain-I and d(CAG)₅.d(CAG)₅ at the major groove side as well as between MSH3-domain-I and MSH3-domain-IV. Markedly, such enhanced interactions are absent in Mutsβ-d(CAG)₂(CAG)(CAG)₂.d(CTG)₂(CAG)(CTG)₂ that has a single A…A mismatch. Thus, the above-mentioned enhancement in intra- and inter- molecular interactions in Mutsβ-d(CAG)₅.d(CAG)₅ provide the stereochemical rationale for the trapping of Mutsβ in CAG repeat expansion disorders.

Secondary structural choice of DNA and RNA associated with CGG/CCG trinucleotide repeat expansion rationalizes the RNA misprocessing in FXTAS

CGG tandem repeat expansion in the 5'-untranslated region of the fragile X mental retardation-1 (FMR1) gene leads to unusual nucleic acid conformations, hence causing genetic instabilities. We show that the number of G…G (in CGG repeat) or C…C (in CCG repeat) mismatches (other than A…T, T…A, C…G and G…C canonical base pairs) dictates the secondary structural choice of the sense and antisense strands of the FMR1 gene and their corresponding transcripts in fragile X-associated tremor/ataxia syndrome (FXTAS). The circular dichroism (CD) spectra and electrophoretic mobility shift assay (EMSA) reveal that CGG DNA (sense strand of the FMR1 gene) and its transcript favor a quadruplex structure. CD, EMSA and molecular dynamics (MD) simulations also show that more than four C…C mismatches cannot be accommodated in the RNA duplex consisting of the CCG repeat (antisense transcript); instead, it favors an i-motif conformational intermediate. Such a preference for unusual secondary structures provides a convincing justification for the RNA foci formation due to the sequestration of RNA-binding proteins to the bidirectional transcripts and the repeat-associated non-AUG translation that are observed in FXTAS. The results presented here also suggest that small molecule modulators that can destabilize FMR1 CGG DNA and RNA quadruplex structures could be promising candidates for treating FXTAS.

Spontaneous and frequent conformational dynamics induced by A…A mismatch in d(CAA)·d(TAG) duplex

Base pair mismatches in DNA can erroneously be incorporated during replication, recombination, etc. Here, the influence of A…A mismatch in the context of 5′CAA·5′TAG sequence is explored using molecular dynamics (MD) simulation, umbrella sampling MD, circular dichroism (CD), microscale thermophoresis (MST) and NMR techniques. MD simulations reveal that the A…A mismatch experiences several transient events such as base flipping, base extrusion, etc. facilitating B–Z junction formation. A…A mismatch may assume such conformational transitions to circumvent the effect of nonisostericity with the flanking canonical base pairs so as to get accommodated in the DNA. CD and 1D proton NMR experiments further reveal that the extent of B–Z junction increases when the number of A…A mismatch in d(CAA)·d(T(A/T)G) increases (1–5). CD titration studies of d(CAA)·d(TAG)_n=5 with the hZα_ADAR1 show the passive binding between the two, wherein, the binding of protein commences with B–Z junction recognition. Umbrella sampling simulation indicates that the mismatch samples anti…+ syn/+ syn…anti, anti…anti & + syn…+ syn glycosyl conformations. The concomitant spontaneous transitions are: a variety of hydrogen bonding patterns, stacking and minor or major groove extrahelical movements (with and without the engagement of hydrogen bonds) involving the mismatch adenines. These transitions frequently happen in anti…anti conformational region compared with the other three regions as revealed from the lifetime of these states. Further, 2D-NOESY experiments indicate that the number of cross-peaks diminishes with the increasing number of A…A mismatches implicating its dynamic nature. The spontaneous extrahelical movement seen in A…A mismatch may be a key pre-trapping event in the mismatch repair due to the accessibility of the base(s) to the sophisticated mismatch repair machinery.

Sequence dependent influence of an A…A mismatch in a DNA duplex: An insight into the recognition by hZαADAR1 protein

Base pair mismatches can erroneously be incorporated in the DNA. An adenine pairing with another adenine is one of the eight possible mismatches. The atomistic insights about the structure and dynamics of an A…A mismatch in a DNA (unbound form) is not yet accessible to any experimental technique. Earlier molecular dynamics (MD) simulations have shown that A…A mismatch in the midst of 5'CAG/3'GAC, 5'GAC/3'CAG and 5'CAA/3'GAT (underline represents the mismatch) are highly dynamic in nature. By employing MD simulation, the influence of an A…A mismatch in the midst of 5'GAA/3'CAT, 5'GAG/3'CAC, 5'AAC/3'TAG, 5'AAG/3'TAC, 5'TAA/3'AAT, 5'TAT/3'AAA and 5'AAT/3'TAA sequences have been investigated here. The results indicate that irrespective of the flanking sequences, the mismatch samples a variety of transient conformations, including a B-Z junction. Further, circular dichroism studies have been carried out to explore the ability of these sequences to bind with hZα_ADAR1 which specifically recognizes B-Z junction/Z-DNA. The results indicate that hZα_ADAR1 could not lead to a complete B to Z transition in the above sequences. Notably, a complete transition to Z-form has been reported earlier for 5'GAC/3'CAG upon titrating with hZα_ADAR1. Intriguingly, 5'AAC/3'TAG, 5'AAG/3'TAC and 5'GAA/3'CAT exhibit a B-Z junction formation rather than a complete transition to Z-form, similar to the situation of 5'CAA/3'GAT. These indicate that although A…A mismatch could induce a local B-Z junction transiently, hZα_ADAR1 requires the presence of a G…C/C…G base pair adjacent to the A…A mismatch for the binding. Additionally, the extent of B-Z junction has enhanced upon binding with hZα_ADAR1 in the presence of the A…A mismatch (specifically when CG, CA, AC, GA and AG steps occur), but not in the presence of the canonical base pairs. These confirm the inclination of A…A mismatch towards the B-Z junction.

Somatic mutations in neurodegeneration: An update

Mosaicism, the presence of genomic differences between cells due to post-zygotic somatic mutations, is widespread in the human body, including within the brain. A role for this in neurodegenerative diseases has long been hypothesised, and technical developments are now allowing the question to be addressed in detail. The rapidly accumulating evidence is discussed in this review, with a focus on recent developments. Somatic mutations of numerous types may occur, including single nucleotide variants (SNVs), copy number variants (CNVs), and retrotransposon insertions. They could act as initiators or risk factors, especially if they arise in development, although they could also result from the disease process, potentially contributing to progression. In common sporadic neurodegenerative disorders, relevant mutations have been reported in synucleinopathies, comprising somatic gains of SNCA in Parkinson's disease and multiple system atrophy, and in Alzheimer's disease, where a novel recombination mechanism leading to somatic variants of APP, as well as an excess of somatic SNVs affecting tau phosphorylation, have been reported. In Mendelian repeat expansion disorders, mosaicism due to somatic instability, first detected 25 years ago, has come to the forefront. Brain somatic SNVs occur in DNA repair disorders, and there is evidence for a role of several ALS genes in DNA repair. While numerous challenges, and need for further validation, remain, this new, or perhaps rediscovered, area of research has the potential to transform our understanding of neurodegeneration.

Characterization and application of fluidic properties of trinucleotide repeat sequences by wax-on-plastic microfluidics

Trinucleotide repeat (TNR) sequences introduce sequence-directed flexibility in the genomic makeup of all living species leading to unique non-canonical structure formation. In humans, the expansions of TNR sequences are responsible for almost 24 neurodegenerative and neuromuscular diseases because their unique structures disrupt cell functions. The biophysical studies of these sequences affect their electrophoretic mobility and spectroscopic signatures. Here, we demonstrate a novel strategy to characterize and discriminate the TNR sequences by monitoring their capillary flow in the absence of an external driving force using wax-on-plastic microchannels. The wax-on-plastic microfluidic system translates the sequence-directed flexibility of TNR into differential flow dynamics. Several variables were used to characterize sequences including concentration, single- vs. double-stranded samples, type of repeat sequence, length of the repeat sequence, presence of mismatches in duplex, and presence of metal ion. All these variables were found to influence the flow velocities of TNR sequences as these factors directly affect the structural flexibility of TNR at the molecular level. An overall trend was observed as the higher flexibility in the TNR structure leads to lower capillary flow. After testing samples derived from relevant cells harboring expanded TNR sequences, it is concluded that this approach may transform into a reagent-free and pump-free biosensing platform to detect microsatellite expansion diseases.

MSH3 modifies somatic instability and disease severity in Huntington's and myotonic dystrophy type 1

The mismatch repair gene MSH3 has been implicated as a genetic modifier of the CAG·CTG repeat expansion disorders Huntington's disease and myotonic dystrophy type 1. A recent Huntington's disease genome-wide association study found rs557874766, an imputed single nucleotide polymorphism located within a polymorphic 9 bp tandem repeat in MSH3/DHFR, as the variant most significantly associated with progression in Huntington's disease. Using Illumina sequencing in Huntington's disease and myotonic dystrophy type 1 subjects, we show that rs557874766 is an alignment artefact, the minor allele for which corresponds to a three-repeat allele in MSH3 exon 1 that is associated with a reduced rate of somatic CAG·CTG expansion (P = 0.004) and delayed disease onset (P = 0.003) in both Huntington's disease and myotonic dystrophy type 1, and slower progression (P = 3.86 × 10-7) in Huntington's disease. RNA-Seq of whole blood in the Huntington's disease subjects found that repeat variants are associated with MSH3 and DHFR expression. A transcriptome-wide association study in the Huntington's disease cohort found increased MSH3 and DHFR expression are associated with disease progression. These results suggest that variation in the MSH3 exon 1 repeat region influences somatic expansion and disease phenotype in Huntington's disease and myotonic dystrophy type 1, and suggests a common DNA repair mechanism operates in both repeat expansion diseases.

Z-DNA in the genome: from structure to disease

The scope of studies investigating the architecture of genomic DNA has progressed steadily since the elucidation of the structure of B-DNA. In recent years, several non-canonical DNA structures including Z-DNA, G-quadruplexes, H-DNA, cruciform DNA, and i-motifs have been reported to form in genomic DNA and are closely related to the evolution and development of disease. The ability of these structures to form in genomic DNA indicates that they might have important cellular roles and are therefore retained during evolution. Understanding the impact of the formation of these secondary structures on cellular processes can enable identification of new targets for therapeutics. In this review, we report the state of understanding of Z-DNA structure and formation and their implication in disease. Finally, we state our perspective on the potential of Z-DNA as a therapeutic target.

The central role of DNA damage and repair in CAG repeat diseases

Diseases such as Huntington's disease and certain spinocerebellar ataxias are caused by the expansion of genomic cytosine-adenine-guanine (CAG) trinucleotide repeats beyond a specific threshold. These diseases are all characterised by neurological symptoms and central neurodegeneration, but our understanding of how expanded repeats drive neuronal loss is incomplete. Recent human genetic evidence implicates DNA repair pathways, especially mismatch repair, in modifying the onset and progression of CAG repeat diseases. Repair pathways might operate directly on repeat sequences by licensing or inhibiting repeat expansion in neurons. Alternatively, or in addition, because many of the genes containing pathogenic CAG repeats encode proteins that themselves have roles in the DNA damage response, it is possible that repeat expansions impair specific DNA repair pathways. DNA damage could then accrue in neurons, leading to further expansion at repeat loci, thus setting up a vicious cycle of pathology. In this review, we consider DNA damage and repair pathways in postmitotic neurons in the context of disease-causing CAG repeats. Investigating and understanding these pathways, which are clearly relevant in promoting and ameliorating disease in humans, is a research priority, as they are known to modify disease and therefore constitute prevalidated drug targets.

A B-Z junction induced by an A … A mismatch in GAC repeats in the gene for cartilage oligomeric matrix protein promotes binding with the hZαADAR1 protein

GAC repeat expansion from five to seven in the exonic region of the gene for cartilage oligomeric matrix protein (COMP) leads to pseudoachondroplasia, a skeletal abnormality. However, the molecular mechanism by which GAC expansions in the COMP gene lead to skeletal dysplasias is poorly understood. Here we used molecular dynamics simulations, which indicate that an A … A mismatch in a d(GAC)6·d(GAC)6 duplex induces negative supercoiling, leading to a local B-to-Z DNA transition. This transition facilitates the binding of d(GAC)7·d(GAC)7 with the Zα-binding domain of human adenosine deaminase acting on RNA 1 (ADAR1, hZαADAR1), as confirmed by CD, NMR, and microscale thermophoresis studies. The CD results indicated that hZαADAR1 recognizes the zigzag backbone of d(GAC)7·d(GAC)7 at the B-Z junction and subsequently converts it into Z-DNA via the so-called passive mechanism. Molecular dynamics simulations carried out for the modeled hZαADAR1-d(GAC)6d(GAC)6 complex confirmed the retention of previously reported important interactions between the two molecules. These findings suggest that hZαADAR1 binding with the GAC hairpin stem in COMP can lead to a non-genetic, RNA editing-mediated substitution in COMP that may then play a crucial role in the development of pseudoachondroplasia.

Structure and Dynamics of DNA and RNA Double Helices of CAG and GAC Trinucleotide Repeats

CAG trinucleotide repeats are known to cause 10 late-onset progressive neurodegenerative disorders as the repeats expand beyond a threshold, whereas GAC repeats are associated with skeletal dysplasias and expand from the normal five to a maximum of seven repeats. The TR secondary structure is believed to play a role in CAG expansions. We have carried out free energy and molecular dynamics studies to determine the preferred conformations of the A-A noncanonical pairs in (CAG)_n and (GAC)_n trinucleotide repeats (n = 1, 4) and the consequent changes in the overall structure of the RNA and DNA duplexes. We find that the global free energy minimum corresponds to A-A pairs stacked inside the core of the helix with anti-anti conformations in RNA and (high-anti)-(high-anti) conformations in DNA. The next minimum corresponds to anti-syn conformations, whereas syn-syn conformations are higher in energy. Transition rates of the A-A conformations are higher for RNA than DNA. Mechanisms for these various transitions are identified. Additional structural and dynamical aspects of the helical conformations are explored, with a focus on contrasting CAG and GAC duplexes. The neutralizing ion distribution around the noncanonical pairs is described.

3D-NuS: A Web Server for Automated Modeling and Visualization of Non-Canonical 3-Dimensional Nucleic Acid Structures

The inherent conformational flexibility of nucleic acids facilitates the formation of a range of conformations such as duplex, triplex, quadruplex, etc. that play crucial roles in biological processes. Elucidation of the influence of non-canonical base pair mismatches on DNA/RNA structures at different sequence contexts to understand the mismatch repair, misregulation of alternative splicing mechanisms and the sequence-dependent effect of RNA-DNA hybrid in relevance to antisense strategy demand their three-dimensional structural information. Furthermore, structural insights about nucleic acid triplexes, which are generally not tractable to structure determination by X-ray crystallography or NMR techniques, are essential to establish their biological function(s). A web server, namely 3D-NuS (http://iith.ac.in/3dnus/), has been developed to generate energy-minimized models of 80 different types of triplexes, 64 types of G-quadruplexes, left-handed Z-DNA/RNA duplexes, and RNA-DNA hybrid duplex along with inter- and intramolecular DNA or RNA duplexes comprising a variety of mismatches and their chimeric forms for any user-defined sequence and length. It also generates an ensemble of conformations corresponding to the modeled structure. These structures may serve as good starting models for docking proteins and small molecules with nucleic acids, NMR structure determination, cryo-electron microscope modeling, DNA/RNA nanotechnology applications and molecular dynamics simulation studies.

Selective Preference of Parallel DNA Triplexes Is Due to the Disruption of Hoogsteen Hydrogen Bonds Caused by the Severe Nonisostericity between the G*GC and T*AT Triplets

Implications of DNA, RNA and RNA.DNA hybrid triplexes in diverse biological functions, diseases and therapeutic applications call for a thorough understanding of their structure-function relationships. Despite exhaustive studies mechanistic rationale for the discriminatory preference of parallel DNA triplexes with G*GC & T*AT triplets still remains elusive. Here, we show that the highest nonisostericity between the G*GC & T*AT triplets imposes extensive stereochemical rearrangements contributing to context dependent triplex destabilisation through selective disruption of Hoogsteen scheme of hydrogen bonds. MD simulations of nineteen DNA triplexes with an assortment of sequence milieu reveal for the first time fresh insights into the nature and extent of destabilization from a single (non-overlapping), double (overlapping) and multiple pairs of nonisosteric base triplets (NIBTs). It is found that a solitary pair of NIBTs, feasible either at a G*GC/T*AT or T*AT/G*GC triplex junction, does not impinge significantly on triplex stability. But two overlapping pairs of NIBTs resulting from either a T*AT or a G*GC interruption disrupt Hoogsteen pair to a noncanonical mismatch destabilizing the triplex by ~10 to 14 kcal/mol, implying that their frequent incidence in multiples, especially, in short sequences could even hinder triplex formation. The results provide (i) an unambiguous and generalised mechanistic rationale for the discriminatory trait of parallel triplexes, including those studied experimentally (ii) clarity for the prevalence of antiparallel triplexes and (iii) comprehensive perspectives on the sequence dependent influence of nonisosteric base triplets useful in the rational design of TFO's against potential triplex target sites.