NMR Dynamics Study Reveals the Zα Domain of Human ADAR1 Associates with and Dissociates from Z-RNA More Slowly than Z-DNA

Preferential formation of Z-RNA over intercalated motifs in long noncoding RNA

Secondary structure is a principal determinant of lncRNA function, predominantly regarding scaffold formation and interfaces with target molecules. Noncanonical secondary structures that form in nucleic acids have known roles in regulating gene expression and include G-quadruplexes (G4s), intercalated motifs (iMs), and R-loops (RLs). In this paper, we used the computational tools G4-iM Grinder and QmRLFS-finder to predict the formation of each of these structures throughout the lncRNA transcriptome in comparison to protein-coding transcripts. The importance of the predicted structures in lncRNAs in biological contexts was assessed by combining our results with publicly available lncRNA tissue expression data followed by pathway analysis. The formation of predicted G4 (pG4) and iM (piM) structures in select lncRNA sequences was confirmed in vitro using biophysical experiments under near-physiological conditions. We find that the majority of the tested pG4s form highly stable G4 structures, and identify many previously unreported G4s in biologically important lncRNAs. In contrast, none of the piM sequences are able to form iM structures, consistent with the idea that RNA is unable to form stable iMs. Unexpectedly, these C-rich sequences instead form Z-RNA structures, which have not been previously observed in regions containing cytosine repeats and represent an interesting and underexplored target for protein-RNA interactions. Our results highlight the prevalence and potential structure-associated functions of noncanonical secondary structures in lncRNAs, and show G4 and Z-RNA structure formation in many lncRNA sequences for the first time, furthering the understanding of the structure-function relationship in lncRNAs.

Zα Domain of ADAR1 Binds to an A-Form-like Nucleic Acid Duplex with Low Micromolar Affinity

The left-handed Z-conformation of nucleic acids can be adopted by both DNA and RNA when bound by Zα domains found within a variety of viral and innate immune response proteins. While Z-form adoption is preferred by certain sequences, such as the commonly studied (CpG)n repeats, Zα has been reported to bind to a wide range of sequence contexts. Studying how Zα interacts with B-/A-form helices prior to their conversion to the Z-conformation is challenging as binding coincides with Z-form adoption. Here, we studied the binding of Zα fromHomo sapiens ADAR1 to a locked "A-type" version of the (CpG)3 construct (LNA (CpG)3) where the sugar pucker is locked into the C3'-endo/C2'-exo conformation, which prevents the duplex from adopting the alternating C2'/C3'-endo sugar puckers found in the Z-conformation. Using NMR and other biophysical techniques, we find that ZαADAR1 binds to the LNA (CpG)3 using a similar interface as for Z-form binding, with a dissociation constant (KD) of ∼4 μM. In contrast to Z-DNA/Z-RNA, where two ZαADAR1 bind to every 6 bp stretch, our data suggests that ZαADAR1 binds to multiple LNA molecules, indicating a completely different binding mode. Because ZαADAR1 binds relatively tightly to a non-Z-form model, its binding to B/A-form helices may need to be considered when experiments are carried out which attempt to identify the Z-form targets of Zα domains. The use of LNA constructs may be beneficial in experiments where negative controls for Z-form adoption are needed.

Z-Form Adoption of Nucleic Acid is a Multi-Step Process Which Proceeds through a Melted Intermediate

The left-handed Z-conformation of nucleic acids can be adopted by both DNA and RNA when bound by Zα domains found within a variety of innate immune response proteins. Zα domains stabilize this higher-energy conformation by making specific interactions with the unique geometry of Z-DNA/Z-RNA. However, the mechanism by which a right-handed helix contorts to become left-handed in the presence of proteins, including the intermediate steps involved, is poorly understood. Through a combination of nuclear magnetic resonance (NMR) and other biophysical measurements, we have determined that in the absence of Zα, under low salt conditions at room temperature, d(CpG) and r(CpG) constructs show no observable evidence of transient Z-conformations greater than 0.5% on either the intermediate or slow NMR time scales. At higher temperatures, we observed a transient unfolded intermediate. The ease of melting a nucleic acid duplex correlates with Z-form adoption rates in the presence of Zα. The largest contributing factor to the activation energies of Z-form adoption as calculated by Arrhenius plots is the ease of flipping the sugar pucker, as required for Z-DNA and Z-RNA. Together, these data validate the previously proposed "zipper model" for Z-form adoption in the presence of Zα. Overall, Z-conformations are more likely to be adopted by double-stranded DNA and RNA regions flanked by less stable regions and by RNAs experiencing torsional/mechanical stress.

Z-RNA biology: a central role in the innate immune response?

Z-RNA is a higher-energy, left-handed conformation of RNA, whose function has remained elusive. A growing body of work alludes to regulatory roles for Z-RNA in the immune response. Here, we review how Z-RNA features present in cellular RNAs-especially containing retroelements-could be recognized by a family of winged helix proteins, with an impact on host defense. We also discuss how mutations to specific Z-contacting amino acids disrupt their ability to stabilize Z-RNA, resulting in functional losses. We end by highlighting knowledge gaps in the field, which, if addressed, would significantly advance this active area of research.

NMR Titration Studies in Z-DNA Dynamics

Chemical shift perturbation (CSP) is a simple NMR technique for studying the DNA binding of proteins. Titration of the unlabeled DNA into the ¹⁵N-labeled protein is monitored by acquiring a two-dimensional (2D) heteronuclear single-quantum correlation (HSQC) spectrum at each step of the titration. CSP can also provide information on the DNA-binding dynamics of proteins, as well as protein-induced conformational changes in DNA. Here, we describe the titration of DNA for the ¹⁵N-labeled Z-DNA-binding protein, monitored via 2D HSQC spectra. NMR titration data can be analyzed with the active B-Z transition model to provide the protein-induced B-Z transition dynamics of DNA.

Thermodynamic analysis of Zα domain-nucleic acid interactions

DNA/RNA molecules adopting the left-handed conformation (Z-form) have been attributed with immunogenic properties. However, their biological role and importance has been a topic of debate for many years. The discovery of Z-DNA/RNA binding domains (Zα domains) in varied proteins that are involved in the innate immune response, such as the interferon inducible form of the RNA editing enzyme ADAR1 (p150), Z-DNA binding protein 1 (ZBP1), the fish kinase PKZ and the poxvirus inhibitor of interferon response E3L, indicates important roles of Z-DNA/RNA in immunity and self/non-self-discrimination. Such Zα domain-containing proteins recognise left-handed Z-DNA/RNA in a conformation-specific manner. Recent studies have implicated these domains in virus recognition. Given these important emerging roles for the Zα domains, it is pivotal to understand the mechanism of recognition of the Z-DNA/Z-RNA by these domains. To this end, we assessed the binding thermodynamics of Zα domain from ORF112 and ADAR1 on T(CG)3 and T(CG)6 oligonucleotides which have high propensity to adopt the Z-conformation. Our study highlights important differences in the mode of oligonucleotide binding by the two Zα domains originating from different proteins. Site-directed mutagenesis was employed together with isothermal titration calorimetry to tease apart finer details of the binding thermodynamics. Our work advances the understanding on binding thermodynamics of Zα domains to their cognate nucleic acid substrates and paves the ground for future efforts to gain a complete appreciation of this process.

Conformational exchange of the Zα domain of human RNA editing enzyme ADAR1 studied by NMR spectroscopy

Z-DNA binding proteins (ZBPs) play important roles in RNA editing, innate immune responses, and viral infections. Numerous studies have implicated a role for conformational motions during ZBPs binding upon DNA, but the quantitative intrinsic conformational exchanges of ZBP have not been elucidated. To understand the correlation between the biological function and dynamic feature of the Zα domains of human ADAR1 (hZα_ADAR1), we have performed the ¹⁵N backbone amide Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion experiments on the free hZα_ADAR1 at two different magnetic fields at 35 °C. The robust inter-dependence of parameters in the global fitting process using multi-magnetic field CPMG profiles allows us characterizing the dynamic properties of conformational changes in hZα_ADAR1. This study found that free hZα_ADAR1 exhibited the conformational exchange with a k_ex of 5784 s^-1 between the states "A" (89% population) and "B" (11% population). The different hydrophobic interactions among helices α1, α2, and α3 between these two states might correlate with efficient Z-DNA binding achieved by the hydrogen bonding interactions between its side-chains and the phosphate backbone of Z-DNA.

Regulatory role of Non-canonical DNA Polymorphisms in human genome and their relevance in Cancer

DNA has the ability to form polymorphic structures like canonical duplex DNA and non-canonical triplex DNA, Cruciform, Z-DNA, G-quadruplex (G4), i-motifs, and hairpin structures. The alteration in the form of DNA polymorphism in the response to environmental changes influences the gene expression. Non-canonical structures are engaged in various biological functions, including chromatin epigenetic and gene expression regulation via transcription and translation, as well as DNA repair and recombination. The presence of non-canonical structures in the regulatory region of the gene alters the gene expression and affects the cellular machinery. Formation of non-canonical structure in the regulatory site of cancer-related genes either inhibits or dysregulate the gene function and promote tumour formation. In the current article, we review the influence of non-canonical structure on the regulatory mechanisms in human genome. Moreover, we have also discussed the relevance of non-canonical structures in cancer and provided information on the drugs used for their treatment by targeting these structures.

Systematic Approach to Find the Global Minimum of Relaxation Dispersion Data for Protein-Induced B-Z Transition of DNA

Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion spectroscopy is commonly used for quantifying conformational changes of protein in μs-to-ms timescale transitions. To elucidate the dynamics and mechanism of protein binding, parameters implementing CPMG relaxation dispersion results must be appropriately determined. Building an analytical model for multi-state transitions is particularly complex. In this study, we developed a new global search algorithm that incorporates a random search approach combined with a field-dependent global parameterization method. The robust inter-dependence of the parameters carrying out the global search for individual residues (GSIR) or the global search for total residues (GSTR) provides information on the global minimum of the conformational transition process of the Zα domain of human ADAR1 (hZαADAR1)-DNA complex. The global search results indicated that a α-helical segment of hZαADAR1 provided the main contribution to the three-state conformational changes of a hZαADAR1-DNA complex with a slow B-Z exchange process. The two global exchange rate constants, kex and kZB, were found to be 844 and 9.8 s-1, respectively, in agreement with two regimes of residue-dependent chemical shift differences-the "dominant oscillatory regime" and "semi-oscillatory regime". We anticipate that our global search approach will lead to the development of quantification methods for conformational changes not only in Z-DNA binding protein (ZBP) binding interactions but also in various protein binding processes.

Recognition of non-CpG repeats in Alu and ribosomal RNAs by the Z-RNA binding domain of ADAR1 induces A-Z junctions

Adenosine-to-inosine (A-to-I) editing of eukaryotic cellular RNAs is essential for protection against auto-immune disorders. Editing is carried out by ADAR1, whose innate immune response-specific cytoplasmic isoform possesses a Z-DNA binding domain (Zα) of unknown function. Zα also binds to CpG repeats in RNA, which are a hallmark of Z-RNA formation. Unexpectedly, Zα has been predicted — and in some cases even shown — to bind to specific regions within mRNA and rRNA devoid of such repeats. Here, we use NMR, circular dichroism, and other biophysical approaches to demonstrate and characterize the binding of Zα to mRNA and rRNA fragments. Our results reveal a broad range of RNA sequences that bind to Zα and adopt Z-RNA conformations. Binding is accompanied by destabilization of neighboring A-form regions which is similar in character to what has been observed for B-Z-DNA junctions. The binding of Zα to non-CpG sequences is specific, cooperative and occurs with an affinity in the low micromolar range. This work allows us to propose a model for how Zα could influence the RNA binding specificity of ADAR1.

ADAR1 is an interferon-induced enzyme that catalyzes editing of adenine to inosine across the transcriptome as part of the immune response. Here the authors establish how ADAR1 recognizes non-CpG RNA sequences to facilitate the formation of A-Z junctions.

The Role of the Z-DNA Binding Domain in Innate Immunity and Stress Granules

Both DNA and RNA can maintain left-handed double helical Z-conformation under physiological condition, but only when stabilized by Z-DNA binding domain (ZDBD). After initial discovery in RNA editing enzyme ADAR1, ZDBD has also been described in pathogen-sensing proteins ZBP1 and PKZ in host, as well as virulence proteins E3L and ORF112 in viruses. The host-virus antagonism immediately highlights the importance of ZDBD in antiviral innate immunity. Furthermore, Z-RNA binding has been shown to be responsible for the localization of these ZDBD-containing proteins to cytoplasmic stress granules that play central role in coordinating cellular response to stresses. This review sought to consolidate current understanding of Z-RNA sensing in innate immunity and implore possible roles of Z-RNA binding within cytoplasmic stress granules.

Protein-induced B-Z transition of DNA duplex containing a 2'-OMe guanosine

Structural transformation of the canonical right-handed helix, B-DNA, to the non-canonical left-handed helix, Z-DNA, can be induced by the Zα domain of the human RNA editing enzyme ADAR1 (hZα_ADAR1). To characterize the site-specific preferences of binding and structural changes in DNA containing the 2'-O-methyl guanosine derivative (_mG), titration of the imino proton spectra and chemical shift perturbations were performed on hZα_ADAR1 upon binding to Z-DNA. The structural transition between B-Z conformation as the changing ratio between DNA and protein showed a binding affinity of the modified DNA onto the Z-DNA binding protein similar to wild-type DNA or RNA. The chemical shift perturbation results showed that the overall structure and environment of the modified DNA revealed DNA-like properties rather than RNA-like characteristics. Moreover, we found evidence for two distinct regimes, "Z-DNA Sensing" and "Modification Sensing", based on the site-specific chemical shift perturbation between the DNA (or RNA) binding complex and the modified DNA-hZα_ADAR1 complex. Thus, we propose that modification of the sugar backbone of DNA with 2'-O-methyl guanosine promotes the changes in the surrounding α3 helical structural segment as well as the non-perturbed feature of the β-hairpin region.

Dynamic regulation of Z-DNA in the mouse prefrontal cortex by the RNA-editing enzyme Adar1 is required for fear extinction

DNA forms conformational states beyond the right-handed double-helix; however, the functional relevance of these non-canonical structures in the brain remains unknown. We show that, in the prefrontal cortex of mice, the formation of one such structure, Z-DNA, is involved in the regulation of extinction memory. Z-DNA is formed during fear learning, and reduced during extinction learning, which is mediated, in part, by a direct interaction between Z-DNA and the RNA editing enzyme Adar1. Adar1 binds to Z-DNA during fear extinction learning which leads to a reduction in Z-DNA at sites where Adar1 is recruited. Knockdown of Adar1 leads to an inability to modify a previously acquired fear memory and blocks activity-dependent changes in DNA structure and RNA state; effects that are fully rescued by the introduction of full-length Adar1. These findings suggest a novel mechanism of learning-induced gene regulation dependent on both proteins which recognize DNA structure, and the state.

Dynamics Studies of DNA with Non-canonical Structure Using NMR Spectroscopy

The non-canonical structures of nucleic acids are essential for their diverse functions during various biological processes. These non-canonical structures can undergo conformational exchange among multiple structural states. Data on their dynamics can illustrate conformational transitions that play important roles in folding, stability, and biological function. Here, we discuss several examples of the non-canonical structures of DNA focusing on their dynamic characterization by NMR spectroscopy: (1) G-quadruplex structures and their complexes with target proteins; (2) i-motif structures and their complexes with proteins; (3) triplex structures; (4) left-handed Z-DNAs and their complexes with various Z-DNA binding proteins. This review provides insight into how the dynamic features of non-canonical DNA structures contribute to essential biological processes.