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

Concentration and time-dependent amyloidogenic characteristics of intrinsically disordered N-terminal region of Saccharomyces cerevisiae Stm1

Saccharomyces cerevisiae Stm1 protein is a ribosomal association factor, which plays an important role in preserving ribosomes in a nutrition-deprived environment. It is also shown to take part in apoptosis-like cell death. Stm1 N-terminal region (Stm1_N1-113) is shown to recognize purine motif DNA triplex and G-quadruplex. Circular dichroism (CD) spectra of Stm1_N1-113 (enriched in positively-charged Lysine and Arginine; negatively-charged Aspartate; polar-uncharged Threonine, Asparagine, Proline and Serine; hydrophobic Alanine, Valine, and Glycine) collected after 0 and 24 h indicate that the protein assumes beta-sheet conformation at the higher concentrations in contrast to intrinsically disordered conformation seen for its monomeric form found in the crystal structure. Thioflavin-T kinetics experiments indicate that the lag phase is influenced by the salt concentration. Atomic force microscopy (AFM) images collected for a variety of Stm1_N1-113 concentrations (in the range of 1-400 μM) in the presence of 150 mM NaCl at 0, 24, and 48 h indicate a threshold concentration requirement to observe the time-dependent amyloid formation. This is prominent seen at the physiological salt concentration of 150 mM NaCl with the fibrillation observed for 400 μM concentration at 48 h, whereas oligomerization or proto-fibrillation is seen for the other concentrations. Such concentration-dependent fibrillation of Stm1_N1-113 explains that amyloid fibrils formed during the overexpression of Stm1_N1-113 may act as a molecular device to trigger apoptosis-like cell death.

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.

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.

Z-DNA and Z-RNA in human disease

Left-handed Z-DNA/Z-RNA is bound with high affinity by the Zα domain protein family that includes ADAR (a double-stranded RNA editing enzyme), ZBP1 and viral orthologs regulating innate immunity. Loss-of-function mutations in ADAR p150 allow persistent activation of the interferon system by Alu dsRNAs and are causal for Aicardi-Goutières Syndrome. Heterodimers of ADAR and DICER1 regulate the switch from RNA- to protein-centric immunity. Loss of DICER1 function produces age-related macular degeneration, a different type of Alu-mediated disease. The overlap of Z-forming sites with those for the signal recognition particle likely limits invasion of primate genomes by Alu retrotransposons.