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.

An Evolutionarily Conserved Strategy for Ribosome Binding and Host Translation Inhibition by β-coronavirus Non-structural Protein 1

An important pathogenicity factor of SARS-CoV-2 and related coronaviruses is Non-structural protein 1 (Nsp1), which suppresses host gene expression and stunts antiviral signaling. SARS-CoV-2 Nsp1 binds the ribosome to inhibit translation through mRNA displacement and induces degradation of host mRNAs. Here we show that Nsp1-dependent host shutoff is conserved in diverse coronaviruses, but only Nsp1 from β-Coronaviruses (β-CoV) inhibits translation through ribosome binding. The C-terminal domain (CTD) of all β-CoV Nsp1s confers high-affinity ribosome binding despite low sequence conservation. Modeling of interactions of four Nsp1s with the ribosome identified the few absolutely conserved amino acids that, together with an overall conservation in surface charge, form the β-CoV Nsp1 ribosome-binding domain. Contrary to previous models, the Nsp1 ribosome-binding domain is an inefficient translation inhibitor. Instead, the Nsp1-CTD likely functions by recruiting Nsp1's N-terminal "effector" domain. Finally, we show that a cis-acting viral RNA element has co-evolved to fine-tune SARS-CoV-2 Nsp1 function, but does not provide similar protection against Nsp1 from related viruses. Together, our work provides new insight into the diversity and conservation of ribosome-dependent host-shutoff functions of Nsp1, knowledge that could aid future efforts in pharmacological targeting of Nsp1 from SARS-CoV-2 and related human-pathogenic β-CoVs. Our study also exemplifies how comparing highly divergent Nsp1 variants can help to dissect the different modalities of this multi-functional viral protein.

SASH1 interacts with TNKS2 and promotes human melanocyte stem cell maintenance

Both aging spots (hyperpigmentation) and hair graying (lack of pigmentation) are associated with aging, two seemingly opposite pigmentation phenotypes. It is not clear how they are mechanistically connected. This study investigated the underlying mechanism in a family with an inherited pigmentation disorder. Clinical examinations identified accelerated hair graying and skin dyspigmentation (intermixed hyper and hypopigmentation) in the family members carrying the SASH1 S519N variant. Cell assays indicated that SASH1 promoted stem-like characteristics in human melanocytes, and SASH1 S519N was defective in this function. Multiple assays showed that SASH1 binds to tankyrase 2 (TNKS2), which is required for SASH1's promotion of stem-like function. Further, the SASH1 S519N variant is in a bona fide Tankyrase-binding motif, and SASH1 S519N alters the binding kinetics and affinity. Results here indicate SASH1 as a novel protein regulating the appropriate balance between melanocyte stem cells (McSC) and mature melanocytes (MCs), with S519N variant causing defects. We propose that dysfunction of McSC maintenance connects multiple aging-associated pigmentation phenotypes in the general population.

Solution NMR Backbone Assignment of the C-Terminal Region of Human Dynein Light Intermediate Chain 2 (LIC2-C) Unveils Structural Resemblance with Its Homologue LIC1-C

Dynein, a homodimeric protein complex, plays a pivotal role in retrograde transportation along microtubules within cells. It consists of various subunits, among which the light intermediate chain (LIC) performs diverse functions, including cargo adaptor binding. In contrast to the vertebrate LIC homolog LIC1, LIC2 has received relatively limited characterization thus far, despite partially orthogonal functional roles. In this study, we present a near-to-complete backbone NMR chemical shift assignment of the C-terminal region of the light intermediate chain 2 of human dynein 1 (LIC2-C). We perform a comparative analysis of the secondary structure propensity of LIC2-C with the one previously reported for LIC1-C and show that the two transient helices in LIC1 that interact with motor adaptors are also present in LIC2.

The Structural Dynamics, Complexity of Interactions, and Functions in Cancer of Multi-SAM Containing Proteins

SAM domains are crucial mediators of diverse interactions, including those important for tumorigenesis or metastasis of cancers, and thus SAM domains can be attractive targets for developing cancer therapies. This review aims to explore the literature, especially on the recent findings of the structural dynamics, regulation, and functions of SAM domains in proteins containing more than one SAM (multi-SAM containing proteins, MSCPs). The topics here include how intrinsic disorder of some SAMs and an additional SAM domain in MSCPs increase the complexity of their interactions and oligomerization arrangements. Many similarities exist among these MSCPs, including their effects on cancer cell adhesion, migration, and metastasis. In addition, they are all involved in some types of receptor-mediated signaling and neurology-related functions or diseases, although the specific receptors and functions vary. This review also provides a simple outline of methods for studying protein domains, which may help non-structural biologists to reach out and build new collaborations to study their favorite protein domains/regions. Overall, this review aims to provide representative examples of various scenarios that may provide clues to better understand the roles of SAM domains and MSCPs in cancer in general.

Solution NMR backbone assignment of the SASH1 SLy proteins associated disordered region (SPIDER)

SASH1 is a scaffold protein with context-dependent biological functions in cell adhesion, tumor metastasis, lung development, and pigmentation. As a member of the SLy protein family, it contains the conserved SLY, SH3, and SAM domains. The 19 kDa SLY domain harbors over 70% of the SASH1 variants associated with pigmentation disorders. However, its solution structure or dynamics have not been investigated yet, and its exact position in the sequence is not clearly defined. Based on the bioinformatic and experimental evidence, we propose renaming this region to the SLy Proteins Associated Disordered Region (SPIDER) and defining the exact position to be amino acids 400-554 of SASH1. We have previously identified a variant in this region linked to a pigmentation disorder, S519N. Here, we used a novel deuteration technique, a suite of TROSY-based 3D NMR experiments, and a high-quality HNN to obtain near complete solution backbone assignment of SASH1's SPIDER. A comparison with the chemical shifts of non-variant (S519) SPIDER shows that the S519N substitution does not alter the free form solution structural propensities of SPIDER. This assignment is the first step to characterize the role of SPIDER in SASH1-mediated cellular functions and provides a model for the future study of sister SPIDER domains in the SLy protein family.

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.

Differential Structural Features of Two Mutant ADAR1p150 Zα Domains Associated with Aicardi-Goutières Syndrome

The Zα domain of ADARp150 is critical for proper Z-RNA substrate binding and is a key factor in the type-I interferon response pathway. Two point-mutations in this domain (N173S and P193A), which cause neurodegenerative disorders, are linked to decreased A-to-I editing in disease models. To understand this phenomenon at the molecular level, we biophysically and structurally characterized these two mutated domains, revealing that they bind Z-RNA with a decreased affinity. Less efficient binding to Z-RNA can be explained by structural changes in beta-wing, part of the Z-RNA-protein interface, and alteration of conformational dynamics of the proteins.

Proteolytic activation of human-specific Olduvai domains by the furin protease

Olduvai protein domains (formerly DUF1220) show the greatest human-specific increase in copy number of any coding region in the genome and are highly correlated with human brain evolution and cognitive disease. The majority of human copies are found within four NBPF genes organized in a variable number of a tandemly arranged three-domain blocks called Olduvai triplets. Here we show that these human-specific Olduvai domains are posttranslationally processed by the furin protease, with a cleavage site occurring once at each triplet. These findings suggest that all expanded human-specific NBPF genes encode proproteins consisting of many independent Olduvai triplet proteins which are activated by furin processing. The exceptional correlation of Olduvai copy number and brain size taken together with our new furin data, indicates the ultimate target of selection was a rapid increase in dosage of autonomously functioning Olduvai triplet proteins, and that these proteins are the primary active agent underlying Olduvai's role in human brain expansion.

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.

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.

Adoption of A-Z Junctions in RNAs by Binding of Zα Domains

While DNA and RNA helices often adopt the canonical B- or A-conformation, the fluid conformational landscape of nucleic acids allows for many higher energy states to be sampled. One such state is the Z-conformation of nucleic acids, which is unique in that it is left-handed and has a "zigzag" backbone. The Z-conformation is recognized and stabilized by Z-DNA/RNA binding domains called Zα domains. We recently demonstrated that a wide range of RNAs can adopt partial Z-conformations termed "A-Z junctions" upon binding to Zα and that the formation of such conformations may be dependent upon both sequence and context. In this chapter, we present general protocols for characterizing the binding of Zα domains to A-Z junction-forming RNAs for the purpose of determining the affinity and stoichiometry of interactions as well as the extent and location of Z-RNA formation.

Slow methyl axes motions in perdeuterated villin headpiece subdomain probed by cross-correlated NMR relaxation measurements

Protein methyl groups can participate in multiple motional modes on different time scales. Sub-nanosecond to nano-second time scale motions of methyl axes are particularly challenging to detect for small proteins in solutions. In this work we employ NMR relaxation interference between the methyl H-H/H-C dipole-dipole interactions [Sun&Tugarinov, J. Magn. Reason. 2012] to characterize methyl axes motions as a function of temperature in a small model protein villin headpiece subdomain (HP36), in which all non-exchangeable protons are deuterated with the exception of methyl groups of leucine and valine residues. The data points to the existence of slow motional modes of methyl axes on sub-nanosecond to nanosecond time scales. Further, at high temperatures for which the overall tumbling of the protein is on the order of 2 ns, we observe a coupling between the slow internal motion and the overall molecular tumbling, based on the anomalous order parameters and their temperature-dependent trends. The addition of 28%(w/w) glycerol-d8 increases the viscosity of the solvent and separates the timescales of internal and overall tumbling, thus permitting for another view of the necessity of the coupling assumption for these sites at high temperatures.

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.

Atomic resolution protein allostery from the multi-state structure of a PDZ domain

Recent methodological advances in solution NMR allow the determination of multi-state protein structures and provide insights into structurally and dynamically correlated protein sites at atomic resolution. This is demonstrated in the present work for the well-studied PDZ2 domain of protein human tyrosine phosphatase 1E for which protein allostery had been predicted. Two-state protein structures were calculated for both the free form and in complex with the RA-GEF2 peptide using the exact nuclear Overhauser effect (eNOE) method. In the apo protein, an allosteric conformational selection step comprising almost 60% of the domain was detected with an "open" ligand welcoming state and a "closed" state that obstructs the binding site by changing the distance between the β-sheet 2, α-helix 2, and sidechains of residues Lys38 and Lys72. The observed induced fit-type apo-holo structural rearrangements are in line with the previously published evolution-based analysis covering ~25% of the domain with only a partial overlap with the protein allostery of the open form. These presented structural studies highlight the presence of a dedicated highly optimized and complex dynamic interplay of the PDZ2 domain owed by the structure-dynamics landscape.

Promiscuity of response regulators for thioredoxin steers bacterial virulence

The exquisite specificity between a sensor kinase and its cognate response regulator ensures faithful partner selectivity within two-component pairs concurrently firing in a single bacterium, minimizing crosstalk with other members of this conserved family of paralogous proteins. We show that conserved hydrophobic and charged residues on the surface of thioredoxin serve as a docking station for structurally diverse response regulators. Using the OmpR protein, we identify residues in the flexible linker and the C-terminal β-hairpin that enable associations of this archetypical response regulator with thioredoxin, but are dispensable for interactions of this transcription factor to its cognate sensor kinase EnvZ, DNA or RNA polymerase. Here we show that the promiscuous interactions of response regulators with thioredoxin foster the flow of information through otherwise highly dedicated two-component signaling systems, thereby enabling both the transcription of Salmonella pathogenicity island-2 genes as well as growth of this intracellular bacterium in macrophages and mice.

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.

Ligand-specific conformational change drives interdomain allostery in Pin1

Pin1 is a two-domain cell regulator that isomerizes peptidyl-prolines. The catalytic domain (PPIase) and the other ligand-binding domain (WW) sample extended and compact conformations. Ligand binding changes the equilibrium of the interdomain conformations, but the conformational changes that lead to the altered domain sampling were unknown. Prior evidence has supported an interdomain allosteric mechanism. We recently introduced a magnetic resonance-based protocol that allowed us to determine the coupling of intra- and interdomain structural sampling in apo Pin1. Here, we describe ligand-specific conformational changes that occur upon binding of pCDC25c and FFpSPR. pCDC25c binding doubles the population of the extended states compared to the virtually identical populations of the apo and FFpSPR-bound forms. pCDC25c binding to the WW domain triggers conformational changes to propagate via the interdomain interface to the catalytic site, while FFpSPR binding displaces a helix in the PPIase that leads to repositioning of the PPIase catalytic loop.

On the use of residual dipolar couplings in multi-state structure calculation of two-domain proteins

Residual dipolar couplings (RDCs) are powerful nuclear magnetic resonance (NMR) probes for the structure calculation of biomacromolecules. Typically, an alignment tensor that defines the orientation of the entire molecule relative to the magnetic field is determined either before refinement of individual bond vectors or simultaneously with this refinement. For single-domain proteins this approach works well since all bond vectors can be described within the same coordinate frame, which is given by the alignment tensor. However, novel approaches are sought after for systems where no universal alignment tensor can be used. Here, we present an approach that can be applied to two-domain proteins that enables the calculation of multiple states within each domain as well as with respect to the relative positions of the two domains.

Solution NMR backbone assignments of disordered Olduvai protein domain CON1 employing Hα-detected experiments

Olduvai protein domains, encoded by the NBPF gene family, are responsible for the largest increase in copy number of any protein-coding region in the human genome. This has spawned various genetics studies which have linked these domains to human brain development and divergence from our primate ancestors, as well as currently relevant cognitive diseases such as schizophrenia and autism spectrum disorder (ASD). There are six separate Olduvai domains which together form the majority of the various protein products of the NBPF genes. The six domains include three conserved domains (CON1-3), and three human-lineage-specific domains (HLS1-3) which occur in triplet. Here, we present the solution nuclear magnetic resonance backbone assignments for the CON1 domain, which has been linked to the severity of ASD. The data confirm that CON1 is an intrinsically disordered protein (IDP). Additionally, we use innovative Hα-detected experiments which allow us to not only assign the Hα atoms and N atoms of proline residues, but also to assign residues where HN-experiments suffered from peak overlap or broadening.

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.

Solution NMR backbone assignments of the N-terminal Zα-linker-Zβ segment from Homo sapiens ADAR1p150

Adenosine-to-inosine (A-to-I) editing of a subset of RNAs in a eukaryotic cell is required in order to avoid triggering the innate immune system. Editing is carried out by ADAR1, which exists as short (p110) and long (p150) isoforms. ADAR1p150 is mostly cytoplasmic, possesses a Z-RNA binding domain (Zα), and is only expressed during the innate immune response. A structurally homologous domain to Zα, the Zβ domain, is separated by a long linker from Zα on the N-terminus of ADAR1 but its function remains unknown. Zβ does not bind to RNA in isolation, yet the binding kinetics of the segment encompassing Zα, Zβ and the 95-residue linker between the two domains (Zα-Zβ) are markedly different compared to Zα alone. Here we present the solution NMR backbone assignment of Zα-Zβ from H. Sapiens ADAR1. The predicted secondary structure of Zα-Zβ based on chemical shifts is in agreement with previously determined structures of Zα and Zβ in isolation, and indicates that the linker is intrinsically disordered. Comparison of the chemical shifts between the individual Zα and Zβ domains to the full Zα-Zβ construct suggests that Zβ may interact with the linker, the function of which is currently unknown.

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.

The Inherent Dynamics and Interaction Sites of the SARS-CoV-2 Nucleocapsid N-Terminal Region

The nucleocapsid protein is one of four structural proteins encoded by SARS-CoV-2 and plays a central role in packaging viral RNA and manipulating the host cell machinery, yet its dynamic behavior and promiscuity in nucleotide binding has made standard structural methods to address its atomic-resolution details difficult. To begin addressing the SARS-CoV-2 nucleocapsid protein interactions with both RNA and the host cell along with its dynamic behavior, we have specifically focused on the folded N-terminal domain (NTD) and its flanking regions using nuclear magnetic resonance solution studies. Studies performed here reveal a large repertoire of interactions, which includes a temperature-dependent self-association mediated by the disordered flanking regions that also serve as binding sites for host cell cyclophilin-A while nucleotide binding is largely mediated by the central NTD core. NMR studies that include relaxation experiments have revealed the complicated dynamic nature of this viral protein. Specifically, while much of the N-terminal core domain exhibits micro-millisecond motions, a central β-hairpin shows elevated inherent flexibility on the pico-nanosecond timescale and the serine/arginine-rich region of residues 176-209 undergoes multiple exchange phenomena. Collectively, these studies have begun to reveal the complexities of the nucleocapsid protein dynamics and its preferred interaction sites with its biological targets.

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.

Microbiota-derived butyrate is an endogenous HIF prolyl hydroxylase inhibitor

The gut microbiota is essential for human health. Microbial supply of short-chain fatty acids (SCFAs), particularly butyrate, is a well-established contributor to gut homeostasis and disease resistance. Reaching millimolar luminal concentrations, butyrate is sequestered and utilized in the colon as the favored energy source for intestinal epithelia. Given the steep oxygen gradient across the anoxic lumen and the highly oxygenated lamina propria, the colon provides a particularly interesting environment to study oxygen sensing. Previous studies have shown that the transcription factor hypoxia-inducible factor (HIF) is stabilized in healthy colonic epithelia. Here we show that butyrate directly inhibits HIF prolyl hydroxylases (PHDs) to stabilize HIF. We find that butyrate stabilizes HIF in vitro despite eliminating β-oxidation and resultant oxygen consumption. Using recombinant PHD protein in combination with nuclear magnetic resonance and enzymatic biochemical assays, we identify butyrate to bind and function as a unique, noncompetitive inhibitor of PHDs relative to other SCFAs. Butyrate inhibited PHD with a noncompetitive Ki of 5.3 ± 0.5 mM, a physiologically relevant concentration. We also confirm that microbiota-derived butyrate is necessary to stabilize HIF in mice colonic tissue through antibiotic-induced butyrate depletion and reconstitution experiments. Our results suggest that the co-evolution of mammals and mutualistic microbiota has selected for butyrate to impact a critical gene regulation pathway that can be extended beyond the mammalian gut. As PHDs are a major target for drug development in the stabilization of HIF, butyrate holds great potential as a well-tolerated endogenous inhibitor with far-reaching therapeutic impact.

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.

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 Motional Details Revealed by Complementary Structural Biology Techniques

Proteins depend on defined molecular plasticity for their functionality. How to comprehensively capture dynamics correctly is of ubiquitous biological importance. Approaches commonly used to probe protein dynamics include model-free elucidation of site-specific motion by NMR relaxation, molecular dynamics (MD)-based approaches, and capturing the substates within a dynamic ensemble by recent eNOE-based multiple-structure approaches. Even though MD is sometimes combined with ensemble-averaged NMR restraints, these approaches have largely been developed and used individually. Owing to the different underlying concepts and practical requirements, it has remained unclear how they compare, and how they cross-validate and complement each other. Here, we extract and compare the differential information contents of MD simulations, NMR relaxation measurements, and eNOE-based multi-state structures for the SH3 domain of chicken α-spectrin. The data show that a validated, consistent, and detailed picture is feasible both for timescales and actual conformational states sampled in the dynamic ensemble. This includes the biologically important side-chain plasticity, for which experimentally cross-validated assessment is a significant challenge.

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.

Deuteration of nonexchangeable protons on proteins affects their thermal stability, side-chain dynamics, and hydrophobicity

We have investigated the effect of deuteration of non-exchangeable protons on protein global thermal stability, hydrophobicity, and local flexibility using well-known thermostable model systems such as the villin headpiece subdomain (HP36) and the third immunoglobulin G-binding domain of protein G (GB3). Reversed-phase high-performance liquid chromatography (RP-HPLC) measurements as a function of temperature probe global thermal stability in the presence of acetonitrile, while differential scanning calorimetry determines thermal stability in solution. Both indicate small but measurable changes in the order of several degrees. RP-HPLC also permitted quantification of the effect of deuteration of just three core phenylalanine side chains of HP36. NMR dynamics investigation has focused on methyl axes motions using cross-correlated relaxation measurements. The analysis of order parameters provided a complex picture indicating that deuteration generally increases motional amplitudes of sub-nanosecond motion in GB3 but decreases those in HP36. Combined with earlier dynamics measurements at C_α -C_β sites and backbone sites of GB3, which probed slower time scales, the results point to the need to probe multiple atoms in the protein and variety of time scales to the discern the full complexity of the effects of deuteration on dynamics.

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.

Solution NMR backbone assignment reveals interaction-free tumbling of human lineage-specific Olduvai protein domains

Olduvai protein domains, encoded primarily by NBPF genes, have been linked to both human brain evolution and cognitive diseases such as autism and schizophrenia. There are six primary domains that comprise the Olduvai family: three conserved domains (CON1-3) and three human lineage-specific domains (HLS1-3), which typically occur as a triplet (HLS1, HLS2 and HLS3). Herein, we present the solution NMR assignment of the backbone chemical shifts of the separate HLS1, 2 and 3 domains of NBPF15. Our data suggest that there is no change in the structure of the separate domains when compared to the full-length triplet (HLS1-HLS2-HLS3). We also demonstrate that there is no direct interaction between the three domains.

Exact distance measurements for structure and dynamics in solid proteins by fast-magic-angle-spinning NMR

Fast-magic-angle-spinning solid-state NMR is a developing technique for determination of protein structure and dynamics. Proton-proton correlations usually lead to rough distance restraints, a serious hurdle towards high-resolution structures. Analogous to the "eNOE" concept in solution, an integrative approach for more accurate restraints enables improved structural accuracy with minimal analytical effort.

Backbone and side-chain chemical shift assignments of full-length, apo, human Pin1, a phosphoprotein regulator with interdomain allostery

Pin1 is a human peptidyl-prolyl cis-trans isomerase important for the regulation of phosphoproteins that are implicated in many diseases including cancer and Alzheimer's. Further biophysical study of Pin1 will elucidate the importance of the two-domain system to regulate its own activity. Here, we report near-complete backbone and side-chain ¹H, ¹³C and ¹⁵N NMR chemical shift assignments of full-length, apo Pin1 for the purpose of studying interdomain allostery and dynamics.

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.

A highly conserved mechanism links the microtubule minus end–directed motor dynein to structurally diverse cargo adaptors through its light intermediate chain; this interaction is crucial for dynein function in vivo.

The large size and complex organization of animal cells make the correct and efficient distribution of intracellular content a challenge. The solution is to use motor proteins, which harness energy from ATP hydrolysis to walk along actin filaments or microtubules, for directional transport of cargo. The multi-subunit motor cytoplasmic dynein 1 (dynein) is responsible for transport directed toward the minus ends of microtubules. An important question is how dynein is recruited to its diverse cargo, which includes organelles such as endosomes and mitochondria, proteins, and mRNA. In this study, we use nuclear magnetic resonance spectroscopy to show that the light intermediate chain (LIC) subunit of human dynein uses a short helix in its disordered C-terminal region to bind structurally distinct adaptor proteins that connect the motor to specific cargo. We then use genome editing in the animal model C. elegans to demonstrate the functional relevance of the C-terminal LIC helix for dynein-dependent cargo transport in neurons. Thus, dynein recruitment to cargo involves a highly conserved interaction between LIC and adaptor proteins.

Structure and dynamics conspire in the evolution of affinity between intrinsically disordered proteins

In every established species, protein-protein interactions have evolved such that they are fit for purpose. However, the molecular details of the evolution of new protein-protein interactions are poorly understood. We have used nuclear magnetic resonance spectroscopy to investigate the changes in structure and dynamics during the evolution of a protein-protein interaction involving the intrinsically disordered CREBBP (CREB-binding protein) interaction domain (CID) and nuclear coactivator binding domain (NCBD) from the transcriptional coregulators NCOA (nuclear receptor coactivator) and CREBBP/p300, respectively. The most ancient low-affinity "Cambrian-like" [540 to 600 million years (Ma) ago] CID/NCBD complex contained less secondary structure and was more dynamic than the complexes from an evolutionarily younger "Ordovician-Silurian" fish ancestor (ca. 440 Ma ago) and extant human. The most ancient Cambrian-like CID/NCBD complex lacked one helix and several interdomain interactions, resulting in a larger solvent-accessible surface area. Furthermore, the most ancient complex had a high degree of millisecond-to-microsecond dynamics distributed along the entire sequences of both CID and NCBD. These motions were reduced in the Ordovician-Silurian CID/NCBD complex and further redistributed in the extant human CID/NCBD complex. Isothermal calorimetry experiments show that complex formation is enthalpically favorable and that affinity is modulated by a largely unfavorable entropic contribution to binding. Our data demonstrate how changes in structure and motion conspire to shape affinity during the evolution of a protein-protein complex and provide direct evidence for the role of structural, dynamic, and frustrational plasticity in the evolution of interactions between intrinsically disordered proteins.

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 ¹³C'-¹⁵N and ¹³C_α-¹³C_β 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 ¹³C'-¹⁵N pairs and into the slow to intermediate regime for the ¹³C_α-¹³C_β 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 ¹³C_α-¹³C_β pairs.

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.

Extending the Applicability of Exact Nuclear Overhauser Enhancements to Large Proteins and RNA

Distance-dependent nuclear Overhauser enhancements (NOEs) are one of the most popular and important experimental restraints for calculating NMR structures. Despite this, they are mostly employed as semiquantitative upper distance bounds, and this discards the wealth of information that is encoded in the cross-relaxation rate constant. Information that is lost includes exact distances between protons and dynamics that occur on the sub-millisecond timescale. Our recently introduced exact measurement of the NOE (eNOE) requires little additional experimental effort relative to other NMR observables. So far, we have used eNOEs to calculate multistate ensembles of proteins up to approximately 150 residues. Here, we briefly revisit eNOE methodology and present two new directions for the use of eNOEs: applications to large proteins and RNA.

High-resolution small RNA structures from exact nuclear Overhauser enhancement measurements without additional restraints

RNA not only translates the genetic code into proteins, but also carries out important cellular functions. Understanding such functions requires knowledge of the structure and dynamics at atomic resolution. Almost half of the published RNA structures have been solved by nuclear magnetic resonance (NMR). However, as a result of severe resonance overlap and low proton density, high-resolution RNA structures are rarely obtained from nuclear Overhauser enhancement (NOE) data alone. Instead, additional semi-empirical restraints and labor-intensive techniques are required for structural averages, while there are only a few experimentally derived ensembles representing dynamics. Here we show that our exact NOE (eNOE) based structure determination protocol is able to define a 14-mer UUCG tetraloop structure at high resolution without other restraints. Additionally, we use eNOEs to calculate a two-state structure, which samples its conformational space. The protocol may open an avenue to obtain high-resolution structures of small RNA of unprecedented accuracy with moderate experimental efforts.

Efficient Stereospecific Hβ2/3 NMR Assignment Strategy for Mid-Size Proteins

We present a strategy for stereospecific NMR assignment of H^β2 and H^β3 protons in mid-size proteins (~150 residues). For such proteins, resonance overlap in standard experiments is severe, thereby preventing unambiguous assignment of a large fraction of β-methylenes. To alleviate this limitation, assignment experiments may be run in high static fields, where higher decoupling power is required. Three-bond H^α–H^β J-couplings (³J_Hα–Hβ) are critical for stereospecific assignments of β-methylene protons, and for determining rotameric χ₁ states. Therefore, we modified a pulse sequence designed to measure accurate ³J_Hα–Hβ couplings such that probe heating was reduced, while the decoupling performance was improved. To further increase the resolution, we applied non-uniform sampling (NUS) schemes in the indirect ¹H and ¹³C dimensions. The approach was applied to two medium-sized proteins, odorant binding protein 22 (OBP22; 14.4 kDa) and Pin1 (18.2 kDa), at 900 MHz polarizing fields. The coupling values obtained from NUS and linear sampling were extremely well correlated. However, NUS decreased the overlap of H^β2/3 protons, thus supplying a higher yield of extracted ³J_Hα-Hβ coupling values when compared with linear sampling. A similar effect could be achieved with linear prediction applied to the linearly sampled data prior to the Fourier transformation. Finally, we used ³J_Hα–Hβ couplings from Pin1 in combination with either conventional or exact nuclear Overhauser enhancement (eNOE) restraints to determine the stereospecific assignments of β-methylene protons. The use of eNOEs further increased the fraction of unambiguously assigned resonances when compared with procedures using conventional NOEs.

NOE-Derived Methyl Distances from a 360 kDa Proteasome Complex

Nuclear magnetic resonance spectroscopy is the prime tool to probe structure and dynamics of biomolecules at atomic resolution. A serious challenge for that method is the size limit imposed on molecules to be studied. Standard studies are typically restricted to ca. 30-40 kDa. More recent developments lead to spin relaxation measurements in methyl groups in single proteins or protein complexes as large as a mega-Dalton, which directly allow the extraction of angular information or experiments with paramagnetic samples. However, these probes are all of indirect nature in contrast to the most intuitive and easy-to-interpret structural/dynamics restraint, the internuclear distance, which can be measured by nuclear Overhauser enhancement (NOE). Herein, we demonstrate time-averaged distance measurements on the 360 kDa half proteasome from Thermoplasma acidophilium. The approach is based on exact quantification of the NOE (eNOE). Our findings open up an avenue for such measurements on very large molecules. These restraints will help in a detailed determination of conformational changes upon perturbation such as ligand binding.

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 ¹⁵N/¹H 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.