Novel Z-DNA binding domains in giant viruses.

Z-nucleic acid structures play vital roles in cellular processes and have implications in innate immunity due to their recognition by Zα domains containing proteins (Z-DNA/Z-RNA binding proteins, ZBPs). Although Zα domains have been identified in six proteins, including viral E3L, ORF112, and I73R, as well as, cellular ADAR1, ZBP1, and PKZ, their prevalence across living organisms remains largely unexplored. In this study, we introduce a computational approach to predict Zα domains, leading to the revelation of previously unidentified Zα domain-containing proteins in eukaryotic organisms, including non-metazoan species. Our findings encompass the discovery of new ZBPs in previously unexplored giant viruses, members of the Nucleocytoviricota phylum. Through experimental validation, we confirm the Zα functionality of select proteins, establishing their capability to induce the B-to-Z conversion. Additionally, we identify Zα-like domains within bacterial proteins. While these domains share certain features with Zα domains, they lack the ability to bind to Z-nucleic acids or facilitate the B-to-Z DNA conversion. Our findings significantly expand the ZBP family across a wide spectrum of organisms and raise intriguing questions about the evolutionary origins of Zα-containing proteins. Moreover, our study offers fresh perspectives on the functional significance of Zα domains in virus sensing and innate immunity and opens avenues for exploring hitherto undiscovered functions of ZBPs.

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.

PBRM1 bromodomains associate with RNA to facilitate chromatin association.

PBRM1 is a subunit of the PBAF chromatin remodeling complex, which is mutated in 40-50% of clear cell renal cell carcinoma patients. It is thought to largely function as a chromatin binding subunit of the PBAF complex, but the molecular mechanism underlying this activity is not fully known. PBRM1 contains six tandem bromodomains which are known to cooperate in binding of nucleosomes acetylated at histone H3 lysine 14 (H3K14ac). Here, we demonstrate that the second and fourth bromodomains from PBRM1 also bind nucleic acids, selectively associating with double stranded RNA elements. Disruption of the RNA binding pocket is found to compromise PBRM1 chromatin binding and inhibit PBRM1-mediated cellular growth effects.

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.

Targeting tumor-derived NLRP3 reduces melanoma progression by limiting MDSCs expansion.

Interleukin-1ß (IL-1ß)-mediated inflammation suppresses antitumor immunity, leading to the generation of a tumor-permissive environment, tumor growth, and progression. Here, we demonstrate that nucleotide-binding domain, leucine-rich containing family, pyrin domain-containing-3 (NLRP3) inflammasome activation in melanoma is linked to IL-1ß production, inflammation, and immunosuppression. Analysis of cancer genome datasets (TCGA and GTEx) revealed greater NLRP3 and IL-1ß expression in cutaneous melanoma samples (n = 469) compared to normal skin (n = 324), with a highly significant correlation between NLRP3 and IL-1ß (P < 0.0001). We show the formation of the NLRP3 inflammasome in biopsies of metastatic melanoma using fluorescent resonance energy transfer analysis for NLRP3 and apoptosis-associated speck-like protein containing a CARD. In vivo, tumor-associated NLRP3/IL-1 signaling induced expansion of myeloid-derived suppressor cells (MDSCs), leading to reduced natural killer and CD8+ T cell activity concomitant with an increased presence of regulatory T (Treg) cells in the primary tumors. Either genetic or pharmacological inhibition of tumor-derived NLRP3 by dapansutrile (OLT1177) was sufficient to reduce MDSCs expansion and to enhance antitumor immunity, resulting in reduced tumor growth. Additionally, we observed that the combination of NLRP3 inhibition and anti-PD-1 treatment significantly increased the antitumor efficacy of the monotherapy by limiting MDSC-mediated T cell suppression and tumor progression. These data show that NLRP3 activation in melanoma cells is a protumor mechanism, which induces MDSCs expansion and immune evasion. We conclude that inhibition of NLRP3 can augment the efficacy of anti-PD-1 therapy.

Advances in the exact nuclear Overhauser effect 2018-2022.

The introduction of the exact nuclear Overhauser enhancement (eNOE) methodology to solution-state nuclear magnetic resonance (NMR) spectroscopy results in tighter distance restraints from NOEs than in convention analysis. These improved restraints allow for higher resolution in structure calculation and even the disentanglement of different conformations of macromolecules. While initial work primarily focused on technical development of the eNOE, structural studies aimed at the elucidation of spatial sampling in proteins and nucleic acids were published in parallel prior to 2018. The period of 2018-2022 saw a continued series of technical innovation, but also major applications addressing biological questions. Here, we review both aspects, covering topics from the implementation of non-uniform sampling of NOESY buildups, novel pulse sequences, adaption of the eNOE to solid-state NMR, advances in eNOE data analysis, and innovations in structural ensemble calculation, to applications to protein, RNA, and DNA structure elucidation.

Structure and Formation of Z-DNA and Z-RNA.

Despite structural differences between the right-handed conformations of A-RNA and B-DNA, both nucleic acids adopt very similar, left-handed Z-conformations. In contrast to their structural similarities and sequence preferences, RNA and DNA exhibit differences in their ability to adopt the Z-conformation regarding their hydration shells, the chemical modifications that promote the Z-conformation, and the structure of junctions connecting them to right-handed segments. In this review, we highlight the structural and chemical properties of both Z-DNA and Z-RNA and delve into the potential factors that contribute to both their similarities and differences. While Z-DNA has been extensively studied, there is a gap of knowledge when it comes to Z-RNA. Where such information is lacking, we try and extend the principles of Z-DNA stability and formation to Z-RNA, considering the inherent differences of the nucleic acids.

SAM1 domain of SASH1 harbors distinctive structural heterogeneity.

The sterile alpha motif (SAM) domains are among the most versatile protein domains in biology, and the variety of the oligomerization states contribute to their diverse roles in many diseases. A better understanding of the structure and dynamics of various SAM domains will provide a scientific basis for drug development targeting them. Here, we used SEC-MALS, HPLC, NMR, and other biophysical techniques to characterize the structural features and dynamics of the SAM1 domain in SASH1. SASH1 is a scaffold protein belonging to the same family as SASH3. Unlike the dimerization seen in SASH3's SAM domain, our SEC-MALS and SE-HPLC showed that SAM1 exists primarily as a less compact monomer with a minor oligomer. NMR assignment, relaxation, and exchange experiments revealed the presence of both a disordered monomer and a more structured oligomer with multiple timescale exchange regimes in solution. Mutagenesis and SE-HPLC showed that D663A/T664K substitutions in SAM1 increased its oligomerization. In sum, this study is the first to characterize a disordered structure for a SAM domain, provides additional evidence and framework for the diversity of SAM domains, and identifies a region in SAM1 as a potential starting point to further characterize the structural mechanism of oligomerization of the domain.

Microbial-derived indoles inhibit neutrophil myeloperoxidase to diminish bystander tissue damage.

During episodes of acute inflammation, polymorphonuclear leukocytes (PMNs) are actively recruited to sites of inflammation or injury where they provide anti-microbial and wound-healing functions. One enzyme crucial for fulfilling these functions is myeloperoxidase (MPO), which generates hypochlorous acid from Cl- and hydrogen peroxide. The potential exists, however, that uncontrolled the extracellular generation of hypochlorous acid by MPO can cause bystander tissue damage and inhibit the healing response. Previous work suggests that the microbiota-derived tryptophan metabolites 1H-indole and related molecules ("indoles") are protective during intestinal inflammation, although their precise mechanism of action is unclear. In the present work, we serendipitously discovered that indoles are potent and selective inhibitors of MPO. Using both primary human PMNs and recombinant human MPO in a cell-free system, we revealed that indoles inhibit MPO at physiologic concentrations. Particularly, indoles block the chlorinating activity of MPO, a reliable marker for MPO-associated tissue damage, as measured by coulometric-coupled HPLC. Further, we observed direct interaction between indoles and MPO using the established biochemical techniques microscale thermophoresis and STD-NMR. Utilizing a murine colitis model, we demonstrate that indoles inhibit bystander tissue damage, reflected in decreased colon 3-chlorotyrosine and pro-inflammatory chemokine expression in vivo. Taken together, these results identify microbiota-derived indoles that acts as endogenous immunomodulatory compounds through their actions on MPO, suggesting a symbiotic association between the gut microbiota and host innate immune system. Such findings offer exciting new targets for future pharmacological intervention.

A transient helix in the disordered region of dynein light intermediate chain links the motor to structurally diverse adaptors for cargo transport.

All animal cells use the motor cytoplasmic dynein 1 (dynein) to transport diverse cargo toward microtubule minus ends and to organize and position microtubule arrays such as the mitotic spindle. Cargo-specific adaptors engage with dynein to recruit and activate the motor, but the molecular mechanisms remain incompletely understood. Here, we use structural and dynamic nuclear magnetic resonance (NMR) analysis to demonstrate that the C-terminal region of human dynein light intermediate chain 1 (LIC1) is intrinsically disordered and contains two short conserved segments with helical propensity. NMR titration experiments reveal that the first helical segment (helix 1) constitutes the main interaction site for the adaptors Spindly (SPDL1), bicaudal D homolog 2 (BICD2), and Hook homolog 3 (HOOK3). In vitro binding assays show that helix 1, but not helix 2, is essential in both LIC1 and LIC2 for binding to SPDL1, BICD2, HOOK3, RAB-interacting lysosomal protein (RILP), RAB11 family-interacting protein 3 (RAB11FIP3), ninein (NIN), and trafficking kinesin-binding protein 1 (TRAK1). Helix 1 is sufficient to bind RILP, whereas other adaptors require additional segments preceding helix 1 for efficient binding. Point mutations in the C-terminal helix 1 of Caenorhabditis elegans LIC, introduced by genome editing, severely affect development, locomotion, and life span of the animal and disrupt the distribution and transport kinetics of membrane cargo in axons of mechanosensory neurons, identical to what is observed when the entire LIC C-terminal region is deleted. Deletion of the C-terminal helix 2 delays dynein-dependent spindle positioning in the one-cell embryo but overall does not significantly perturb dynein function. We conclude that helix 1 in the intrinsically disordered region of LIC provides a conserved link between dynein and structurally diverse cargo adaptor families that is critical for dynein function in vivo.

Integrating NMR and simulations reveals motions in the UUCG tetraloop.

We provide an atomic-level description of the structure and dynamics of the UUCG RNA stem-loop by combining molecular dynamics simulations with experimental data. The integration of simulations with exact nuclear Overhauser enhancements data allowed us to characterize two distinct states of this molecule. The most stable conformation corresponds to the consensus three-dimensional structure. The second state is characterized by the absence of the peculiar non-Watson-Crick interactions in the loop region. By using machine learning techniques we identify a set of experimental measurements that are most sensitive to the presence of non-native states. We find that although our MD ensemble, as well as the consensus UUCG tetraloop structures, are in good agreement with experiments, there are remaining discrepancies. Together, our results show that (i) the MD simulation overstabilize a non-native loop conformation, (ii) eNOE data support its presence with a population of ≈10% and (iii) the structural interpretation of experimental data for dynamic RNAs is highly complex, even for a simple model system such as the UUCG tetraloop.

Reconstruction of Coupled Intra- and Interdomain Protein Motion from Nuclear and Electron Magnetic Resonance.

Proteins composed of multiple domains allow for structural heterogeneity and interdomain dynamics that may be vital for function. Intradomain structures and dynamics can influence interdomain conformations and vice versa. However, no established structure determination method is currently available that can probe the coupling of these motions. The protein Pin1 contains separate regulatory and catalytic domains that sample "extended" and "compact" states, and ligand binding changes this equilibrium. Ligand binding and interdomain distance have been shown to impact the activity of Pin1, suggesting interdomain allostery. In order to characterize the conformational equilibrium of Pin1, we describe a novel method to model the coupling between intra- and interdomain dynamics at atomic resolution using multistate ensembles. The method uses time-averaged nuclear magnetic resonance (NMR) restraints and double electron-electron resonance (DEER) data that resolve distance distributions. While the intradomain calculation is primarily driven by exact nuclear Overhauser enhancements (eNOEs), J couplings, and residual dipolar couplings (RDCs), the relative domain distribution is driven by paramagnetic relaxation enhancement (PREs), RDCs, interdomain NOEs, and DEER. Our data support a 70:30 population of the compact and extended states in apo Pin1. A multistate ensemble describes these conformations simultaneously, with distinct conformational differences located in the interdomain interface stabilizing the compact or extended states. We also describe correlated conformations between the catalytic site and interdomain interface that may explain allostery driven by interdomain contact.

Special Issue: A, B and Z: The Structure, Function and Genetics of Z-DNA and Z-RNA.

It is now difficult to believe that a biological function for the left-handed Z-DNA and Z-RNA conformations was once controversial. The papers in this Special Issue, "Z-DNA and Z-RNA: from Physical Structure to Biological Function", are based on presentations at the ABZ2021 meeting that was held virtually on 19 May 2021 and provide evidence for several biological functions of these structures. The first of its kind, this international conference gathered over 200 scientists from many disciplines to specifically address progress in research involving Z-DNA and Z-RNA. These high-energy left-handed conformers of B-DNA and A-RNA are associated with biological functions and disease outcomes, as evidenced from both mouse and human genetic studies. These alternative structures, referred to as "flipons", form under physiological conditions, regulate type I interferon responses and induce necroptosis during viral infection. They can also stimulate genetic instability, resulting in adaptive evolution and diseases such as cancer. The meeting featured cutting-edge science that was, for the most part, unpublished. We plan for the ABZ meeting to reconvene in 2022.

Reducing the measurement time of exact NOEs by non-uniform sampling.

We have previously reported on the measurement of exact NOEs (eNOEs), which yield a wealth of additional information in comparison to conventional NOEs. We have used these eNOEs in a variety of applications, including calculating high-resolution structures of proteins and RNA molecules. The collection of eNOEs is challenging, however, due to the need to measure a NOESY buildup series consisting of typically four NOESY spectra with varying mixing times in a single measurement session. While the 2D version can be completed in a few days, a fully sampled 3D-NOESY buildup series can take 10 days or more to acquire. This can be both expensive as well as problematic in the case of samples that are not stable over such a long period of time. One potential method to significantly decrease the required measurement time of eNOEs is to use non-uniform sampling (NUS) to decrease the number of points measured in the indirect dimensions. The effect of NUS on the extremely tight distance restraints extracted from eNOEs may be very pronounced. Therefore, we investigated the fidelity of eNOEs measured from three test cases at decreasing NUS densities: the 18.4 kDa protein human Pin1, the 4.1 kDa WW domain of Pin1 (both in 3D), and a 4.6 kDa 14mer RNA UUCG tetraloop (2D). Our results show that NUS imparted negligible error on the eNOE distances derived from good quality data down to 10% sampling for all three cases, but there is a noticeable decrease in the eNOE yield that is dependent upon the underlying sparsity, and thus complexity, of the sample. For Pin1, this transition occurred at roughly 40% while for the WW domain and the UUCG tetraloop it occurred at lower NUS densities of 20% and 10%, respectively. We rationalized these numbers through reconstruction simulations under various conditions. The extent of this loss depends upon the number of scans taken as well as the number of peaks to be reconstructed. Based on these findings, we have created guidelines for choosing an optimal NUS density depending on the number of peaks needed to be reconstructed in the densest region of a 2D or 3D NOESY spectrum.

Protein Allostery at Atomic Resolution.

Protein allostery is a phenomenon involving the long range coupling between two distal sites in a protein. In order to elucidate allostery at atomic resoluion on the ligand-binding WW domain of the enzyme Pin1, multistate structures were calculated from exact nuclear Overhauser effect (eNOE). In its free form, the protein undergoes a microsecond exchange between two states, one of which is predisposed to interact with its parent catalytic domain. In presence of the positive allosteric ligand, the equilibrium between the two states is shifted towards domain-domain interaction, suggesting a population shift model. In contrast, the allostery-suppressing ligand decouples the side-chain arrangement at the inter-domain interface thereby reducing the inter-domain interaction. As such, this mechanism is an example of dynamic allostery. The presented distinct modes of action highlight the power of the interplay between dynamics and function in the biological activity of proteins.

Distance-independent Cross-correlated Relaxation and Isotropic Chemical Shift Modulation in Protein Dynamics Studies.

Cross-correlated relaxation (CCR) in multiple-quantum coherences differs from other relaxation phenomena in its theoretical ability to be mediated across an infinite distance. The two interfering relaxation mechanisms may be dipolar interactions, chemical shift anisotropies, chemical shift modulations or quadrupolar interactions. These properties make multiple-quantum CCR an attractive probe for structure and dynamics of biomacromolecules not accessible from other measurements. Here, we review the use of multiple-quantum CCR measurements in dynamics studies of proteins. We compile a list of all experiments proposed for CCR rate measurements, provide an overview of the theory with a focus on protein dynamics, and present applications to various protein systems.

Activity and Affinity of Pin1 Variants.

Pin1 is a peptidyl-prolyl isomerase responsible for isomerizing phosphorylated S/T-P motifs. Pin1 has two domains that each have a distinct ligand binding site, but only its PPIase domain has catalytic activity. Vast evidence supports interdomain allostery of Pin1, with binding of a ligand to its regulatory WW domain impacting activity in the PPIase domain. Many diverse studies have made mutations in Pin1 in order to elucidate interactions that are responsible for ligand binding, isomerase activity, and interdomain allostery. Here, we summarize these mutations and their impact on Pin1's structure and function.

The Sign of Nuclear Magnetic Resonance Chemical Shift Difference as a Determinant of the Origin of Binding Selectivity: Elucidation of the Position Dependence of Phosphorylation in Ligands Binding to Scribble PDZ1.

The use of nuclear magnetic resonance chemical shift perturbation to monitor changes taking place around the binding site of a ligand-protein interaction is a routine and widely applied methodology in the field of protein biochemistry. Shifts are often acquired by titrating various concentrations of ligand to a fixed concentration of the receptor and may serve the purpose, among others, of determining affinity constants, locating binding surfaces, or differentiating between binding mechanisms. Shifts are quantified by the so-called combined chemical shift difference. Although the directionality of shift changes is often used for detailed analysis of specific cases, the approach has not been adapted in standard chemical shift monitoring. This is surprising as it would not require additional effort. Here, we demonstrate the importance of the sign of the chemical shift difference induced by ligand-protein interaction. We analyze the sign of the 15N/1H shift changes of the PDZ1 domain of Scribble upon interaction with two pairs of phosphorylated and unphosphorylated peptides. We find that detailed differences in the molecular basis of this PDZ-ligand interaction can be obtained from our analysis to which the classical method of combined chemical shift perturbation analysis is insensitive. In addition, we find a correlation between affinity and millisecond motions. Application of the methodology to Cyclophilin a, a cis-trans isomerase, reveals molecular details of peptide recognition. We consider our directionality vector chemical shift analysis as a method of choice when distinguishing the molecular origin of binding specificities of a class of similar ligands, which is often done in drug discovery.

Correlated motions of C'-N and Cα-Cß pairs in protonated and per-deuterated GB3.

We investigated correlated µs-ms time scale motions of neighboring 13C'-15N and 13Cα-13Cß nuclei in both protonated and perdeuterated samples of GB3. The techniques employed, NMR relaxation due to cross-correlated chemical shift modulations, specifically target concerted changes in the isotropic chemical shifts of the two nuclei associated with spatial fluctuations. Field-dependence of the relaxation rates permits identification of the parameters defining the chemical exchange rate constant under the assumption of a two-site exchange. The time scale of motions falls into the intermediate to fast regime (with respect to the chemical shift time scale, 100-400 s-1 range) for the 13C'-15N pairs and into the slow to intermediate regime for the 13Cα-13Cß pairs (about 150 s-1). Comparison of the results obtained for protonated and deuterated GB3 suggests that deuteration has a tendency to reduce these slow scale correlated motions, especially for the 13Cα-13Cß pairs.