Sen1 architecture: RNA-DNA hybrid resolution, autoregulation, and insights into SETX inactivation in AOA2.

The senataxin (SETX, Sen1 in yeasts) RNA-DNA hybrid resolving helicase regulates multiple nuclear transactions, including DNA replication, transcription, and DNA repair, but the molecular basis for Sen1 activities is ill defined. Here, Sen1 cryoelectron microscopy (cryo-EM) reconstructions reveal an elongated inchworm-like architecture. Sen1 is composed of an amino terminal helical repeat Sen1 N-terminal (Sen1N) regulatory domain that is flexibly linked to its C-terminal SF1B helicase motor core (Sen1Hel) via an intrinsically disordered tether. In an autoinhibited state, the Sen1Sen1N domain regulates substrate engagement by promoting occlusion of the RNA substrate-binding cleft. The X-ray structure of an activated Sen1Hel engaging single-stranded RNA and ADP-SO4 shows that the enzyme encircles RNA and implicates a single-nucleotide power stroke in the Sen1 RNA translocation mechanism. Together, our data unveil dynamic protein-protein and protein-RNA interfaces underpinning helicase regulation and inactivation of human SETX activity by RNA-binding-deficient mutants in ataxia with oculomotor apraxia 2 neurodegenerative disease.

Coordinated DNA polymerization by Polγ and the region of LonP1 regulated proteolysis.

The replicative mitochondrial DNA polymerase, Polγ, and its protein regulation are essential for the integrity of the mitochondrial genome. The intricacies of Polγ regulation and its interactions with regulatory proteins, which are essential for fine-tuning polymerase function, remain poorly understood. Misregulation of the Polγ heterotrimer, consisting of (i) PolG, the polymerase catalytic subunit and (ii) PolG2, the accessory subunit, ultimately results in mitochondrial diseases. Here, we used single particle cryo-electron microscopy to resolve the structure of PolG in its apoprotein state and we captured Polγ at three intermediates within the catalytic cycle: DNA bound, engaged, and an active polymerization state. Chemical crosslinking mass spectrometry, and site-directed mutagenesis uncovered the region of LonP1 engagement of PolG, which promoted proteolysis and regulation of PolG protein levels. PolG2 clinical variants, which disrupted a stable Polγ complex, led to enhanced LonP1-mediated PolG degradation. Overall, this insight into Polγ aids in an understanding of mitochondrial DNA replication and characterizes how machinery of the replication fork may be targeted for proteolytic degradation when improperly functioning.

Structural insight and characterization of human Twinkle helicase in mitochondrial disease.

Twinkle is the mammalian helicase vital for replication and integrity of mitochondrial DNA. Over 90 Twinkle helicase disease variants have been linked to progressive external ophthalmoplegia and ataxia neuropathies among other mitochondrial diseases. Despite the biological and clinical importance, Twinkle represents the only remaining component of the human minimal mitochondrial replisome that has yet to be structurally characterized. Here, we present 3-dimensional structures of human Twinkle W315L. Employing cryo-electron microscopy (cryo-EM), we characterize the oligomeric assemblies of human full-length Twinkle W315L, define its multimeric interface, and map clinical variants associated with Twinkle in inherited mitochondrial disease. Cryo-EM, crosslinking-mass spectrometry, and molecular dynamics simulations provide insight into the dynamic movement and molecular consequences of the W315L clinical variant. Collectively, this ensemble of structures outlines a framework for studying Twinkle function in mitochondrial DNA replication and associated disease states.

Mechanism by which T7 bacteriophage protein Gp1.2 inhibits Escherichia coli dGTPase.

Levels of the cellular dNTPs, the direct precursors for DNA synthesis, are important for DNA replication fidelity, cell cycle control, and resistance against viruses. Escherichia coli encodes a dGTPase (2'-deoxyguanosine-5'-triphosphate [dGTP] triphosphohydrolase [dGTPase]; dgt gene, Dgt) that establishes the normal dGTP level required for accurate DNA replication but also plays a role in protecting E. coli against bacteriophage T7 infection by limiting the dGTP required for viral DNA replication. T7 counteracts Dgt using an inhibitor, the gene 1.2 product (Gp1.2). This interaction is a useful model system for studying the ongoing evolutionary virus/host "arms race." We determined the structure of Gp1.2 by NMR spectroscopy and solved high-resolution cryo-electron microscopy structures of the Dgt-Gp1.2 complex also including either dGTP substrate or GTP coinhibitor bound in the active site. These structures reveal the mechanism by which Gp1.2 inhibits Dgt and indicate that Gp1.2 preferentially binds the GTP-bound form of Dgt. Biochemical assays reveal that the two inhibitors use different modes of inhibition and bind to Dgt in combination to yield enhanced inhibition. We thus propose an in vivo inhibition model wherein the Dgt-Gp1.2 complex equilibrates with GTP to fully inactivate Dgt, limiting dGTP hydrolysis and preserving the dGTP pool for viral DNA replication.

Flipped over U: structural basis for dsRNA cleavage by the SARS-CoV-2 endoribonuclease.

Coronaviruses generate double-stranded (ds) RNA intermediates during viral replication that can activate host immune sensors. To evade activation of the host pattern recognition receptor MDA5, coronaviruses employ Nsp15, which is a uridine-specific endoribonuclease. Nsp15 is proposed to associate with the coronavirus replication-transcription complex within double-membrane vesicles to cleave these dsRNA intermediates. How Nsp15 recognizes and processes dsRNA is poorly understood because previous structural studies of Nsp15 have been limited to small single-stranded (ss) RNA substrates. Here we present cryo-EM structures of SARS-CoV-2 Nsp15 bound to a 52nt dsRNA. We observed that the Nsp15 hexamer forms a platform for engaging dsRNA across multiple protomers. The structures, along with site-directed mutagenesis and RNA cleavage assays revealed critical insight into dsRNA recognition and processing. To process dsRNA Nsp15 utilizes a base-flipping mechanism to properly orient the uridine within the active site for cleavage. Our findings show that Nsp15 is a distinctive endoribonuclease that can cleave both ss- and dsRNA effectively.

GADGETS: a genetic algorithm for detecting epistasis using nuclear families.

MOTIVATION: Epistasis may play an etiologic role in complex diseases, but research has been hindered because identification of interactions among sets of single nucleotide polymorphisms (SNPs) requires exploration of immense search spaces. Current approaches using nuclear families accommodate at most several hundred candidate SNPs. RESULTS: GADGETS detects epistatic SNP-sets by applying a genetic algorithm to case-parent or case-sibling data. To allow for multiple epistatic sets, island subpopulations of SNP-sets evolve separately under selection for evident joint relevance to disease risk. The software evaluates the identified SNP-sets via permutation testing and provides graphical visualization. GADGETS correctly identified epistatic SNP-sets in realistically simulated case-parent triads with 10 000 candidate SNPs, far more SNPs than competitors can handle, and it outperformed competitors in simulations with many fewer SNPs. Applying GADGETS to family-based oral-clefting data from dbGaP identified SNP-sets with possible epistatic effects on risk. AVAILABILITY AND IMPLEMENTATION: GADGETS is part of the epistasisGA package at https://github.com/mnodzenski/epistasisGA. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Method for the structural analysis of Twinkle mitochondrial DNA helicase by cryo-EM.

The mitochondrial replisome replicates the 16.6 kb mitochondria DNA (mtDNA). The proper functioning of this multicomponent protein complex is vital for the integrity of the mitochondrial genome. One of the critical protein components of the mitochondrial replisome is the Twinkle helicase, a member of the Superfamily 4 (SF4) helicases. Decades of research has uncovered common themes among SF4 helicases including self-assembly, ATP-dependent translocation, and formation of protein-protein complexes. Some of the molecular details of these processes are still unknown for the mitochondria SF4 helicase, Twinkle. Here, we describe a protocol for expression, purification, and single-particle cryo-electron microscopy of the Twinkle helicase clinical variant, W315L, which resulted in the first high-resolution structure of Twinkle helicase. The methods described here serve as an adaptable protocol to support future high-resolution studies of Twinkle helicase or other SF4 helicases.

Characterization of SARS2 Nsp15 nuclease activity reveals it's mad about U.

Nsp15 is a uridine specific endoribonuclease that coronaviruses employ to cleave viral RNA and evade host immune defense systems. Previous structures of Nsp15 from across Coronaviridae revealed that Nsp15 assembles into a homo-hexamer and has a conserved active site similar to RNase A. Beyond a preference for cleaving RNA 3' of uridines, it is unknown if Nsp15 has any additional substrate preferences. Here, we used cryo-EM to capture structures of Nsp15 bound to RNA in pre- and post-cleavage states. The structures along with molecular dynamics and biochemical assays revealed critical residues involved in substrate specificity, nuclease activity, and oligomerization. Moreover, we determined how the sequence of the RNA substrate dictates cleavage and found that outside of polyU tracts, Nsp15 has a strong preference for purines 3' of the cleaved uridine. This work advances our understanding of how Nsp15 recognizes and processes viral RNA, and will aid in the development of new anti-viral therapeutics.

Ubiquitin stimulated reversal of topoisomerase 2 DNA-protein crosslinks by TDP2.

Tyrosyl-DNA phosphodiesterase 2 (TDP2) reverses Topoisomerase 2 DNA-protein crosslinks (TOP2-DPCs) in a direct-reversal pathway licensed by ZATTZNF451 SUMO2 E3 ligase and SUMOylation of TOP2. TDP2 also binds ubiquitin (Ub), but how Ub regulates TDP2 functions is unknown. Here, we show that TDP2 co-purifies with K63 and K27 poly-Ubiquitinated cellular proteins independently of, and separately from SUMOylated TOP2 complexes. Poly-ubiquitin chains of ≥ Ub3 stimulate TDP2 catalytic activity in nuclear extracts and enhance TDP2 binding of DNA-protein crosslinks in vitro. X-ray crystal structures and small-angle X-ray scattering analysis of TDP2-Ub complexes reveal that the TDP2 UBA domain binds K63-Ub3 in a 1:1 stoichiometric complex that relieves a UBA-regulated autoinhibitory state of TDP2. Our data indicates that that poly-Ub regulates TDP2-catalyzed TOP2-DPC removal, and TDP2 single nucleotide polymorphisms can disrupt the TDP2-Ubiquitin interface.

Predicting tumor response to drugs based on gene-expression biomarkers of sensitivity learned from cancer cell lines.

BACKGROUND: Human cancer cell line profiling and drug sensitivity studies provide valuable information about the therapeutic potential of drugs and their possible mechanisms of action. The goal of those studies is to translate the findings from in vitro studies of cancer cell lines into in vivo therapeutic relevance and, eventually, patients' care. Tremendous progress has been made. RESULTS: In this work, we built predictive models for 453 drugs using data on gene expression and drug sensitivity (IC50) from cancer cell lines. We identified many known drug-gene interactions and uncovered several potentially novel drug-gene associations. Importantly, we further applied these predictive models to ~ 17,000 bulk RNA-seq samples from The Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression (GTEx) database to predict drug sensitivity for both normal and tumor tissues. We created a web site for users to visualize and download our predicted data ( https://manticore.niehs.nih.gov/cancerRxTissue ). Using trametinib as an example, we showed that our approach can faithfully recapitulate the known tumor specificity of the drug. CONCLUSIONS: We demonstrated that our approach can predict drugs that 1) are tumor-type specific; 2) elicit higher sensitivity from tumor compared to corresponding normal tissue; 3) elicit differential sensitivity across breast cancer subtypes. If validated, our prediction could have relevance for preclinical drug testing and in phase I clinical design.

Structural characterization of the virulence factor Sda1 nuclease from Streptococcus pyogenes.

Infection by Group A Streptococcus pyogenes (GAS) is a leading cause of severe invasive disease in humans, including streptococcal toxic shock syndrome and necrotizing fasciitis. GAS infections lead to nearly 163,000 annual deaths worldwide. Hypervirulent strains of S. pyogenes have evolved a plethora of virulence factors that aid in disease-by promoting bacterial adhesion to host cells, subsequent invasion of deeper tissues and blocking the immune system's attempts to eradicate the infection. Expression and secretion of the extracellular nuclease Sda1 is advantageous for promoting bacterial dissemination throughout the host organism, and evasion of the host's innate immune response. Here we present two crystal structures of Sda1, as well as biochemical studies to address key structural features and surface residues involved in DNA binding and catalysis. In the active site, Asn211 is observed to directly chelate a hydrated divalent metal ion and Arg124, on the putative substrate binding loop, likely stabilizes the transition state during phosphodiester bond cleavage. These structures provide a foundation for rational drug design of small molecule inhibitors to be used in prevention of invasive streptococcal disease.

A Structural Basis for Biguanide Activity.

Metformin is the most commonly prescribed treatment for type II diabetes and related disorders; however, molecular insights into its mode(s) of action have been limited by an absence of structural data. Structural considerations along with a growing body of literature demonstrating its effects on one-carbon metabolism suggest the possibility of folate mimicry and anti-folate activity. Motivated by the growing recognition that anti-diabetic biguanides may act directly upon the gut microbiome, we have determined structures of the complexes formed between the anti-diabetic biguanides (phenformin, buformin, and metformin) and Escherichia coli dihydrofolate reductase (ecDHFR) based on nuclear magnetic resonance, crystallographic, and molecular modeling studies. Interligand Overhauser effects indicate that metformin can form ternary complexes with p-aminobenzoyl-l-glutamate (pABG) as well as other ligands that occupy the region of the folate-binding site that interacts with pABG; however, DHFR inhibition is not cooperative. The biguanides competitively inhibit the activity of ecDHFR, with the phenformin inhibition constant being 100-fold lower than that of metformin. This inhibition may be significant at concentrations present in the gut of treated individuals, and inhibition of DHFR in intestinal mucosal cells may also occur if accumulation levels are sufficient. Perturbation of folate homeostasis can alter the pyridine nucleotide redox ratios that are important regulators of cellular metabolism.

A comprehensive genomic pan-cancer classification using The Cancer Genome Atlas gene expression data.

BACKGROUND: The Cancer Genome Atlas (TCGA) has generated comprehensive molecular profiles. We aim to identify a set of genes whose expression patterns can distinguish diverse tumor types. Those features may serve as biomarkers for tumor diagnosis and drug development. METHODS: Using RNA-seq expression data, we undertook a pan-cancer classification of 9,096 TCGA tumor samples representing 31 tumor types. We randomly assigned 75% of samples into training and 25% into testing, proportionally allocating samples from each tumor type. RESULTS: We could correctly classify more than 90% of the test set samples. Accuracies were high for all but three of the 31 tumor types, in particular, for READ (rectum adenocarcinoma) which was largely indistinguishable from COAD (colon adenocarcinoma). We also carried out pan-cancer classification, separately for males and females, on 23 sex non-specific tumor types (those unrelated to reproductive organs). Results from these gender-specific analyses largely recapitulated results when gender was ignored. Remarkably, more than 80% of the 100 most discriminative genes selected from each gender separately overlapped. Genes that were differentially expressed between genders included BNC1, FAT2, FOXA1, and HOXA11. FOXA1 has been shown to play a role for sexual dimorphism in liver cancer. The differentially discriminative genes we identified might be important for the gender differences in tumor incidence and survival. CONCLUSIONS: We were able to identify many sets of 20 genes that could correctly classify more than 90% of the samples from 31 different tumor types using TCGA RNA-seq data. This accuracy is remarkable given the number of the tumor types and the total number of samples involved. We achieved similar results when we analyzed 23 non-sex-specific tumor types separately for males and females. We regard the frequency with which a gene appeared in those sets as measuring its importance for tumor classification. One third of the 50 most frequently appearing genes were pseudogenes; the degree of enrichment may be indicative of their importance in tumor classification. Lastly, we identified a few genes that might play a role in sexual dimorphism in certain cancers.

Structure of Escherichia coli dGTP triphosphohydrolase: a hexameric enzyme with DNA effector molecules.

The Escherichia coli dgt gene encodes a dGTP triphosphohydrolase whose detailed role still remains to be determined. Deletion of dgt creates a mutator phenotype, indicating that the dGTPase has a fidelity role, possibly by affecting the cellular dNTP pool. In the present study, we have investigated the structure of the Dgt protein at 3.1-Šresolution. One of the obtained structures revealed a protein hexamer that contained two molecules of single-stranded DNA. The presence of DNA caused significant conformational changes in the enzyme, including in the catalytic site of the enzyme. Dgt preparations lacking DNA were able to bind single-stranded DNA with high affinity (Kd ∼ 50 nM). DNA binding positively affected the activity of the enzyme: dGTPase activity displayed sigmoidal (cooperative) behavior without DNA but hyperbolic (Michaelis-Menten) kinetics in its presence, consistent with a specific lowering of the apparent Km for dGTP. A mutant Dgt enzyme was also created containing residue changes in the DNA binding cleft. This mutant enzyme, whereas still active, was incapable of DNA binding and could no longer be stimulated by addition of DNA. We also created an E. coli strain containing the mutant dgt gene on the chromosome replacing the wild-type gene. The mutant also displayed a mutator phenotype. Our results provide insight into the allosteric regulation of the enzyme and support a physiologically important role of DNA binding.

Selective unfolding of one Ribonuclease H domain of HIV reverse transcriptase is linked to homodimer formation.

HIV-1 reverse transcriptase (RT), a critical enzyme of the HIV life cycle and an important drug target, undergoes complex and largely uncharacterized conformational rearrangements that underlie its asymmetric folding, dimerization and subunit-selective ribonuclease H domain (RH) proteolysis. In the present article we have used a combination of NMR spectroscopy, small angle X-ray scattering and X-ray crystallography to characterize the p51 and p66 monomers and the conformational maturation of the p66/p66' homodimer. The p66 monomer exists as a loosely structured molecule in which the fingers/palm/connection, thumb and RH substructures are connected by flexible (disordered) linking segments. The initially observed homodimer is asymmetric and includes two fully folded RH domains, while exhibiting other conformational features similar to that of the RT heterodimer. The RH' domain of the p66' subunit undergoes selective unfolding with time constant ∼6.5 h, consistent with destabilization due to residue transfer to the polymerase' domain on the p66' subunit. A simultaneous increase in the intensity of resonances near the random coil positions is characterized by a similar time constant. Consistent with the residue transfer hypothesis, a construct of the isolated RH domain lacking the two N-terminal residues is shown to exhibit reduced stability. These results demonstrate that RH' unfolding is coupled to homodimer formation.

Molecular mechanism of substrate specificity for heparan sulfate 2-O-sulfotransferase.

Heparan sulfate (HS) is an abundant polysaccharide in the animal kingdom with essential physiological functions. HS is composed of sulfated saccharides that are biosynthesized through a complex pathway involving multiple enzymes. In vivo regulation of this process remains unclear. HS 2-O-sulfotransferase (2OST) is a key enzyme in this pathway. Here, we report the crystal structure of the ternary complex of 2OST, 3'-phosphoadenosine 5'-phosphate, and a heptasaccharide substrate. Utilizing site-directed mutagenesis and specific oligosaccharide substrate sequences, we probed the molecular basis of specificity and 2OST position in the ordered HS biosynthesis pathway. These studies revealed that Arg-80, Lys-350, and Arg-190 of 2OST interact with the N-sulfo groups near the modification site, consistent with the dependence of 2OST on N-sulfation. In contrast, 6-O-sulfo groups on HS are likely excluded by steric and electrostatic repulsion within the active site supporting the hypothesis that 2-O-sulfation occurs prior to 6-O-sulfation. Our results provide the structural evidence for understanding the sequence of enzymatic events in this pathway.

Heparan sulfate selectively inhibits the collagenase activity of cathepsin K.

Cathepsin K (CtsK) is a cysteine protease with potent collagenase activity. CtsK is highly expressed by bone-resorbing osteoclasts and plays an essential role in resorption of bone matrix. Although CtsK is known to bind heparan sulfate (HS), the structural details of the interaction, and how HS regulates the biological functions of CtsK, remains largely unknown. In this report, we discovered that HS is a multifaceted regulator of the structure and function of CtsK. Structurally, HS forms a highly stable complex with CtsK and induces its dimerization. Co-crystal structures of CtsK with bound HS oligosaccharides reveal the location of the HS binding site and suggest how HS may support dimerization. Functionally, HS plays a dual role in regulating the enzymatic activity of CtsK. While it preserves the peptidase activity of CtsK by stabilizing its active conformation, it inhibits the collagenase activity of CtsK in a sulfation level-dependent manner. These opposing effects can be explained by our finding that the HS binding site is remote from the active site, which allows HS to specifically inhibit the collagenase activity without affecting the peptidase activity. At last, we show that structurally defined HS oligosaccharides effectively block osteoclast resorption of bone in vitro without inhibiting osteoclast differentiation, which suggests that HS-based oligosaccharide might be explored as a new class of selective CtsK inhibitor for many diseases involving exaggerated bone resorption.

Heparan sulfate selectively inhibits the collagenase activity of cathepsin K.

Cathepsin K (CtsK) is a cysteine protease with potent collagenase activity. CtsK is highly expressed by bone-resorbing osteoclasts and plays an essential role in bone remodeling. Although CtsK is known to bind heparan sulfate (HS), the structural details of the interaction, and how HS ultimately regulates the biological functions of CtsK, remains largely unknown. In this report, we determined that CtsK preferably binds to larger HS oligosaccharides, such as dodecasaccharides (12mer), and that the12mer can induce monomeric CtsK to form a stable dimer in solution. Interestingly, while HS has no effect on the peptidase activity of CtsK, it greatly inhibits the collagenase activity of CtsK in a manner dependent on sulfation level. By forming a complex with CtsK, HS was able to preserve the full peptidase activity of CtsK for prolonged periods, likely by stabilizing its active conformation. Crystal structures of Ctsk with a bound 12mer, alone and in the presence of the endogenous inhibitor cystatin-C reveal the location of HS binding is remote from the active site. Mutagenesis based on these complex structures identified 6 basic residues of Ctsk that play essential roles in mediating HS-binding. At last, we show that HS 12mers can effectively block osteoclast resorption of bone in vitro. Combined, we have shown that HS can function as a multifaceted regulator of CtsK and that HS-based oligosaccharide might be explored as a new class of selective CtsK inhibitor in many diseases that involve exaggerated bone resorption.

XDSGUI: a graphical user interface for XDS, SHELX and ARCIMBOLDO.

XDSGUI is a lightweight graphical user interface (GUI) for the XDS, SHELX and ARCIMBOLDO program packages that serves both novice and experienced users in obtaining optimal processing and phasing results for X-ray, neutron and electron diffraction data. The design of the program enables data processing and phasing without command line usage, and supports advanced command flows in a simple user-modifiable and user-extensible way. The GUI supplies graphical information based on the tabular log output of the programs, which is more intuitive, comprehensible and efficient than text output can be.