key: cord-0746230-u4w5q0ar
authors: Hammond, Robert G.; Tan, Xuan; Johnson, Margaret A.
title: SARS‐unique fold in the Rousettus bat coronavirus HKU9
date: 2017-06-15
journal: Protein Science
DOI: 10.1002/pro.3208
sha: 7bd017deb1253a77c74ca78c836acc2d3ae2d179
doc_id: 746230
cord_uid: u4w5q0ar

The coronavirus nonstructural protein 3 (nsp3) is a multifunctional protein that comprises multiple structural domains. This protein assists viral polyprotein cleavage, host immune interference, and may play other roles in genome replication or transcription. Here, we report the solution NMR structure of a protein from the “SARS‐unique region” of the bat coronavirus HKU9. The protein contains a frataxin fold or double‐wing motif, which is an α + β fold that is associated with protein/protein interactions, DNA binding, and metal ion binding. High structural similarity to the human severe acute respiratory syndrome (SARS) coronavirus nsp3 is present. A possible functional site that is conserved among some betacoronaviruses has been identified using bioinformatics and biochemical analyses. This structure provides strong experimental support for the recent proposal advanced by us and others that the “SARS‐unique” region is not unique to the human SARS virus, but is conserved among several different phylogenetic groups of coronaviruses and provides essential functions.

Coronaviruses are single-stranded, positive-sense, enveloped RNA viruses that infect both humans and animals. Coronavirus infections have a range of severity and include upper and lower respiratory symptoms, with a low frequency of acute lung injury and acute respiratory distress syndrome. 1 Acute gastrointestinal, hepatic, and neurological symptoms have also been observed. 2 Since 2002, the human coronaviruses (CoVs) have emerged as significant public health threats. The severe acute respiratory syndrome (SARS) virus is the etiological agent of the 2003-2005 pandemic that affected more than 30 countries. 3 In 2012, the Middle East respiratory syndrome (MERS) virus emerged in the Middle East, followed by the spread of the virus to other countries (e.g., the UK, South Korea). As of 2016, there had been 1728 confirmed cases of MERS affecting persons in 27 countries. 4 Prior to these outbreaks, CoVs were known to be responsible for mild upper and lower respiratory infections. For example, human CoV 229E and OC43 cause a minority of respiratory tract infections. 2 Based on phylogenetic and serological analyses, the International Committee for Taxonomy of Viruses has placed the CoVs in four genera, namely the Alphacoronaviruses, Betacoronaviruses, Gammacoronaviruses and Deltacoronaviruses. 5 Under this classification, the betacoronavirus genus has been divided into groups a to d, whereby the SARS-like CoVs are found in group B and MERS-like CoVs in group C. The group D so far has been detected only in bats. 6 Bats are reservoir hosts of multiple zoonotic viruses, including CoVs. Surveillance studies and phylogenetic analyses have shown that high genetic diversity exists among the SARS-like viruses present in bats, allowing for the possibility of recombination and the evolution of new variants. 7 A bat virus with 96% nucleotide sequence identity to the human SARS-CoV was shown to be capable of using the human ACE2 enzyme as a receptor. This demonstrates the same mode of cell entry as the human SARS-CoV. 8 The bat SL-CoV-WIV1 could grow on human epithelial cells and Vero E6 cells, and was neutralized by human SARS convalescent sera. This virus is a possible direct progenitor of the human SARS-CoV. 8, 9 Several group c betacoronaviruses, such as the HKU4, HKU5, and PREDICT/PDF-2180, have been identified in bats from distinct locations around the world. Some genome regions in these bat viruses are highly conserved with respect to the human MERS virus; for example, PREDICT/PDF-2180 shares 97% sequence identity with the MERS virus in ORF1B. 10 It is hypothesized that RNA recombination either in the bat or in an intermediate animal host gave rise to the MERS-CoV. 10 The HKU4 virus, which is derived from the lesser bamboo bat (Tylonycteris pachypus), shares 92.4% RNA polymerase, 67.4% spike protein, and 72.3% nucleocapsid amino acid identity with the MERS CoV and is able to use the same receptor for attachment and entry (the cell surface protein DPP4). 11, 12 The group D betacoronavirus Hong Kong University 9 (HKU9) is also widely distributed, and has been detected in diverse species including Rousettus leschenaulti, Hipposidereos commersoni, Eidolon helvum, and Rousettus aegyptiacus from Asia to Africa. [13] [14] [15] [16] Whether bat CoVs undergo adaptation to intermediate hosts, or are transmitted directly to humans, it is clear that they pose a threat to human health. Hence, it is imperative to understand bat CoV biochemical and biological functions. At present, only one high-resolution structure of a BatCoV HKU9 protein domain is known, the spike protein external receptor-binding domain (RBD). 17 This structure revealed critical new information such as the external subdomain adopting a helical fold versus the beta-sheet topology observed in other beta-CoV receptor domains. As a result, the HKU9 RBD does not bind to the other betaCoV receptors, ACE2 and CD26, underlining the importance of carrying out structural studies on bat proteins. Hence, we have initiated a program to explore bat protein structure-function relationships, with the goal of determining conserved versus divergent functions.

The CoV virion is composed of four structural proteins, which are believed to assist genome packaging, cell entry and virus spread. 2 In contrast, the replicase gene directs the expression of two large nonstructural polyproteins, pp1a and pp1ab, that become mature nonstructural proteins (nsps) after cleavage by viral proteases. These proteins assemble into a replicase-transcriptase complex (RTC) that is responsible for RNA genome replication, processing and transcription of sub-genomic RNAs. Interference with the innate immune system, and other interactions with functions of the host cell also localize to the nsps. Several of these functions are essential for viral replication, growth and virulence. [18] [19] [20] [21] [22] [23] [24] [25] The nonstructural protein 3 (nsp3) is a multifunctional protein consisting of sixteen functional domains and 1,922 amino acid residues. 18, 21, [26] [27] [28] [29] [30] [31] [32] This protein is the largest component of the RTC. Nsp3 is one of the most divergent regions of the CoV genome. 33 The domain structure of nsp3 is variable among CoVs, 32 with one or two papain-like cysteine proteases, transmembrane regions, RNA-binding proteins, and one or more macrodomains. 27, 34, 35 Key functions of the nsp3 include protein/protein interactions involved in replicase assembly and function; 36 polyprotein processing by the papain-like cysteine protease domain; 37 and deubiquitinase activity involved in innate immune system interference. 38 There are one or more macrodomains in the protein, for which roles in countering the host cell innate immunity have been demonstrated 21, 39 and roles in viral RNA synthesis have been proposed. 40 A "SARSunique region" with a three-domain structure was identified in the nsp3 of SARS. 35 The macrodomains in the SARS-unique region were shown to be Gquadruplex binding proteins, and to interact with the RCHY ubiquitin ligase to target p53 for degradation. 35, 41, 42 The smaller C-terminal domain in this region adopts a frataxin-like fold and has been shown to bind purine-rich RNA sequences. 35 In the human SARS-CoV, the functions of this region were essential for viral replication. 43 However, based on discoveries since 2002 and the emergence of other viruses, it has been hypothesized that the "SARS-unique region" is in fact conserved in other viruses, in particular in the group B, C, and D betacoronaviruses.

We are investigating the "SARS-unique region" of bat CoVs. Here, we report the solution structure of the small C-terminal domain of this region, which we term HKU9 C. We describe for the first time the structural and functional analysis of a nonstructural protein domain from the betacoronavirus lineage D. We also discuss the conserved elements of the nsp3 C domain compared to other proteins in the frataxin fold family; including a possible functional site that is conserved relative to the human SARS-CoV.

NMR experiments were performed with uniformly 15 N, 13 C-labeled HKU9 C expressed and purified from E. coli. The construct used contains the entire predicted C domain spanning the residues 573-646 of the nonstructural protein 3 (nsp3), with an additional N-terminal segment Ser-His-Met derived from fusion tag cleavage. These residues correspond to the residues 1345-1418 of the replicase polyprotein 1ab of BatCoV HKU9 (Uniprot ID: P0C6W5). The numbering differs because the viral polyprotein is cleaved by the viral protease PLpro to yield the mature viral nsp3. 36, [44] [45] [46] We use the numbering of the mature nsp3 herein. Multidimensional NMR experiments were performed to assign 96% of the observable resonances of the peptide backbone and amino acid sidechains. All backbone 15 N and 1 H N resonances were assigned (Fig. 1) . The structure determination was carried out based on 3D 15 Table I displays the statistics of the structure calculation, indicating a high-quality structure determination. A dense network of long-range NOEs was observed and the sequential and medium-range NOE pattern was consistent with the secondary structures in the protein (Supporting Information Fig. S1 ). The ensemble of 20 conformers representing the solution structure of the HKU9 C domain (RMSD 0.34 Å ) is well-defined with the exception of the N-terminal expression tag residues Ser 23 and His 22, and the C-terminal residue Lys 646.

A fold consisting of six b-strands arranged in an antiparallel b-sheet, together with two a-helices at the N-and C-termini that pack on one side of the sheet is observed (Fig. 2 ). The fold is described as a double-wing motif or frataxin-like fold 48 and is classified as similar to the N-terminal domain of CyaY, a bacterial regulatory protein. 49 The helices rest in the same plane antiparallel to each other and contribute to one side of the hydrophobic core [ Fig.  2 (A)]. The two helices, a1 and a2, are comprised of residues 574-585 and 636 2 644, respectively. The first beta strands b1 (591 2 592) and b2 (596 2 599) follow an extended loop after a1 and lead to the first b hairpin. The remaining beta strands b3-b6 span the residues 602 2 609, 613 2 616, 622 2 626, and 629 2 632 forming a curved b-sheet. The topology of the frataxin fold is shown in Figure 2 (C).

The hydrophobic core is primarily defined by residues from the a-helices and b-strands 

Structural alignment of HKU9 C to other proteins using the programs TM-Align 50 and Dali 51 revealed structural similarity to betacoronavirus (b-CoV) C domains, frataxins, and hypothetical proteins (Table  IIA) . The most structurally similar proteins originate from other b-CoV C domains, namely those of the human SARS-CoV and murine hepatitis virus (MHV) C. 30 The HKU9 C fold is similar to these viral domains, with a similar topology and overall backbone RMSD values of 1.7 Å and 2.2 Å , respectively. These viral domains have conserved residues and a highly similar fold despite their low sequence identity. Similarity to the frataxins is also evident, Predicted binding regions of the HKU9 C protein are described: residues from the a-helices (a), residues from the solvent-accessible b-sheet (b), and the conserved polar face containing the Arg-Asp-Trp and Lys-Arg-Gly motif (CF). d The confidence score (C-score) is used to evaluate the reliability of the prediction. The binding site score (BS-score) evaluates how significant is the match between the predicted binding site and the template binding site. Alignment score values range from 0 to 1, with higher values having greater significance.

with RMSD values of approximately 3 Å and 1-10% sequence identity. These proteins also show slightly different topologies, with longer loops and secondary structure insertions between several secondary structure elements [ Fig. 3(C) ]. Functional predictions of HKU9 C were based on an analysis of b-CoV C domain structure-function relationships, together with COACH meta-server results. 56 COACH creates a complementary profile and binding site prediction from TM-SITE and S-SITE and utilizes multiple structure-based programs (COFACTOR, FINDSITE, and Concavity) 56 to derive ligand binding predictions. We used this consensus server approach to predict functional characteristics of HKU9-C (Table IIB) .

Based on similarities to the human SARS-CoV C domain, a possible function for HKU9 C is nucleic acid binding. 35 To investigate this possibility, we conducted electrophoretic mobility shift assays (EMSA) with a panel of RNA and DNA oligonucleotides including purine-rich, pyrimidine-rich and G-quadruplex sequences. However, no oligonucleotide binding was detected. A second possibility is that HKU9 C functions in concert with the neighboring macrodomains, which are binding proteins and enzymes acting on ADP-ribose and related metabolites. [57] [58] [59] Structural similarity and binding site similarity to adenylate-binding proteins is also present. Chemical shift perturbation analysis was employed by titrating to 20 times the protein concentration of ADP and ADP-ribose, which are known ligands for macrodomains. 57 Functional predictions based on binding site analysis suggested other possible ligands. To investigate, chemical shift perturbation experiments were repeated with cyanocobalamin (vitamin B12), zinc (II) ions, EDTA, and peptides. Again, no changes in the spectrum were observed, suggesting other likely functions for HKU9 C.

The structure determination of HKU9 C revealed unexpected structural similarity with the corresponding SARS-unique domain in the human SARS-CoV. These two sequences share only 18% sequence similarity. An area of strong conservation is present around the residues Arg 588-Trp 590 in the loop joining a1 to b1, where the residues are conserved [indicated by stars, Fig. 3(C) ] and the protein Fig. 3(A,B) ]. Additional similarity is present around the residues Lys 609-Gly 611. In particular, the residues Arg 588-Asp 589 -Trp 590 in HKU9 and the residues Arg 670-Asp 671-Trp 672 in SARS adopt nearly identical side chain orientations (Fig. 4) . This suggests a possible conserved function between the two viruses. We describe this surface as the conserved 55 The alignment is based on structural alignments obtained with TM-Align. 50 PDB codes are included after each protein name. The residue numbers for HKU9-C are indicated. Alpha helix regions are displayed in red (cylinders) and beta strands are shown in blue (arrows). Gaps are shown as dashes (-) and insertions where additional secondary structures are present are indicated by forward slash marks (//). Residues indicated by stars (*) discussed in the text are involved in potential functional sites. The corresponding Dali scores for the pairwise alignment of each protein with HKU9 C and the percent amino acid identity between each protein and HKU9 C domain are listed. Dali scores of 2.0 and higher indicate significant sequence identity 51 face (CF) of the protein. This surface is defined by the loop connecting a1 and b1 and the beta turn between b3 and b4, near the C-terminus of the protein.

In contrast, the corresponding region of the MHV C domain has acidic and hydrophobic character [ Fig. 3(A,B) ]. This is a consequence of the substitution of the sequence Arg 588-Asp 589-Trp 590 by Thr-Asp-Trp and Lys 609-Arg 610-Gly 611 by Glu-Cys-Pro. Since the three proteins share a low level of overall sequence identity (15-18%), this difference would not have been apparent without a structural comparison.

This structure represents clear evidence that the SARS-unique domain is also conserved in bat CoVs. The overall structural similarity between the HKU9 C domain and the SARS C domain, from betacoronavirus lineage B, was assessed by the program DALI. 51 The resulting RMSD value was 1.66 Å with a DALI score of 8.2, indicating a strong match. The RMSD value for the MHV C domain from the betacoronavirus lineage A was 2.16 Å , with a DALI score of 8.7. Since the b-CoV HKU9 belongs to lineage D, this analysis reinforces the hypothesis advanced by us and others that the unique region of SARS nsp3 is actually conserved across multiple b-CoVs. 30, 32, 43, 60, 61 A structure-based sequence alignment of the C domains and related proteins is shown in Figure 3 (C). Residues such as Phe 578, Val 581, and Trp 590 are conserved when compared to the sequence and structure of the SARS and MHV orthologues. In contrast, Trp 590 is replaced by other aromatic or hydrophobic residues in the frataxins. Based on their low solvent accessibility, we conclude that these residues are likely to be important in stabilizing the fold, rather than for intermolecular interactions. However, other conserved residues that contribute to the surface potential such as the side chains of Arg 588, Asp 589, Lys 609, Arg 610, and Gly 611, described above, that are oriented to the same face of the protein, are likely to be responsible for a shared function between the CoV groups 2b and 2d (Fig. 3) . Correspondingly, these residues are not conserved throughout the protein family.

The sequence alignment of Figure 3 (C) reveals the conserved topology in the frataxin fold family. It also reveals differences between the viral domains and more distantly related proteins. The viral proteins retain a similar sequence length, have conserved residues in both helices and the b-sheet, and align with high DALI scores of 8.2 and above, indicating a strong match. These features are not conserved in the distantly related frataxin-like folds. For example, the adenylate-binding AcsD domain (PDB ID: 2W04) 55 has an extended loop with an alpha turn insertion between b1 and b2 and another long loop between b3 and b4. DALI scores for the alignment of the HKU9 C domain to the human frataxin (PDB: 3T3X) and to the bacterial frataxin (PDB: 4HS5) are 3.8 and 3.5, respectively. These scores are also significant (>2.0) and indicate a conserved fold, but with some structural variability. 51 This is underscored by the presence of structural insertions relative to the viral proteins.

We employed bioinformatics analysis with the COACH 56 meta-server to predict possible functions for the bat CoV HKU9. Several possible functions emerged from this analysis. One possible function is as a nucleic acid-binding protein, predicted by the COACH and COFACTOR 62 servers (Table II) with low confidence score values of 0.02 and 0.01. This is also highlighted by the sequence and structure alignment of the SARS and HKU9 C domains. Several residues involved in the binding of SARS C to RNA are conserved in HKU9-C. 35 The RNA-binding residues from the SARS-CoV protein, such as His 695 (b52b6 loop), Gly 707 (b62b7 loop) and Val 709 (b7strand), align to Phe 617 (b42b5 loop), Gly 627 (b52b6 loop), and Val 630 (b6) in HKU9 C [ Fig.  3(C) ]. Additionally, a distantly related viral frataxin, the C terminal domain of the T4 activator MotA, (PDB ID: 1KAF), 53 binds an E. coli DNA promoter sequence. The MotA CTD double-wing b-sheet utilizes asparagine residues to bind DNA. These residues are not conserved; for instance, one (Asn 187) aligns to Gly 627 in HKU9 C [ Fig. 3(C) ]. Consistent with this lack of conservation, no nucleic acid binding was observed for HKU9 C. However, it is possible that this function is present but requires the presence of neighboring nsp3 domains.

A second possible function for the HKU9 C domain is protein/protein interaction. In the SARS-CoV, the SUD region interacts with host cell proteins to enhance p53 degradation. 41 The frataxins also have protein binding partners, where the interaction is mediated by side chains that are exposed on the planar face of the b-sheet. It is notable that in the viral proteins, the b-sheet face is smaller and less planar than that of the frataxins. The latter proteins have the b12b2 hairpin in the same plane as the b-sheet [ Fig. 3(A) ]. However, in the b-CoV domains, the b12b2 hairpin wraps over the b-sheet, obscuring the side chains in b22b6 from the protein surface. In addition, the b-sheet side chains that are important to frataxin binding and catalysis such as Trp 155 and Arg 165 (human) or Arg 53 and Trp 61 (Psychromonas ingrahamii) are not conserved in the b-CoV proteins. 52, 54 Protein or peptide binding is another function that is predicted by the bioinformatic analysis of HKU9 C (Table II) . A potential protein-binding site is predicted to be present on the conserved face (CF) of the protein. The site shows structural and chemical similarity to that of the AAA 1 delivery protein 63 and the Nsl1 protein. 64 The surface identified by this prediction includes the residues Arg 588-Asp 589-Trp 590 and Lys 609 -Arg 610 -Gly 611 that we have identified as a conserved functional site.

Analysis of the bioinformatics results displays a theme with respect to HKU9 C surface regions. The conserved face of the fold [ Fig. 3(B) ] is the only region of the protein that was predicted to have protein-protein interactions, while the other surface regions predicted metal ion and small molecule recognition [Table II(B)]. The b-sheet sidechains are not solvent exposed, but side chains from the b22b3 and b42b5 loops and from the b6 strand could potentially bind small molecules. Interestingly, the metal ion ligands such as Ca 21 and Zn 21 were predicted to bind to the a-helices. To date, we were not able to experimentally confirm any metal ion binding activity or nucleic acid binding activity for HKU9 C. The prediction of a possible protein/protein interaction function is intriguing and awaits further experiment.

We hypothesize that the conserved face of the HKU9 C domain is a likely interface for HKU9 C binding partners. Based on our FFAS (Fold and Function Assignment System) analysis, 65 and on the experimental results reported here, we predict structural conservation between the nsp3 proteins of the human SARS-CoV and bat HKU9. We used this structural alignment to predict the linker regions that would join the HKU9 C domain to the neighboring domains in nsp3. At the N-terminus of HKU9 C, a short, three-residue linker is predicted to join the domain to the neighboring M domain; while at the C-terminus, on the "conserved face" of the protein, a seven-residue linker joins the C domain to the papain-like protease of the virus. 36, 44, 66 A longer linker would provide flexibility to accommodate binding partners and interactions. This would coincide with our hypothesis that the conserved face of the protein near the C-terminus may harbor a potential functional site for reactivity or binding to other biomolecules.

The frataxin or double-wing fold of the bat HKU9 nsp3 C domain reported here has high structural similarity to the human SARS-CoV C domain. Although there is low sequence similarity to the other CoV nsp3 proteins, some residues are structurally conserved. The conservation of specific surface polar residues relative to the human SARS virus may indicate a conserved function among certain betacoronaviruses.

The DNA sequence encoding the central region of nsp3 (37 2 1037) was obtained as a codon-optimized synthetic gene from Genscript (Piscataway, NJ). The residues 573 2 646 of nsp3, corresponding to the recombinant HKU9 C domain, were cloned into the vector pET-15b-TEV 67 vector from the Northeast Structural Genomics Consortium (DNASU). The construct was expressed in E. coli strain BL21 (DE3) with a 6xHis tag. The Ser-His-Met sequence at the N-terminus of the proteins remained after tag cleavage with the tobacco etch virus protease. The sample was prepared in both LB medium for natural isotopic abundance and in minimal medium for uniform 15 N-and 13 C-labeling. These samples were used for functional analysis and structure determination, respectively. Sample conditions such as buffer, pH, and salt concentration were optimized based on peak intensity and linewidth in the [ 15 N, 1 H]-HSQC spectrum, leading to the selection of 20 mM sodium phosphate (pH 6.0), 150 mM NaCl, and 5 mM DTT. The protein was monomeric as assessed by sizeexclusion chromatography on a GE Healthcare HiLoad 26/600 Superdex TM 200 pg column. 

The C domain structure was determined based on multidimensional NMR experiments using uniformly 15 N-or [ 15 N, 13 C]-labeled protein solution with 97% H 2 O/3% D 2 O (v/v). All experiments were conducted on Bruker Avance III HD spectrometers (600 and 850 MHz) equipped with Bruker 5 mm TCI cryoprobes and on a Bruker Avance II 700 MHz spectrometer equipped with a CP TCI H-C/N-D cryoprobe. The sequence-specific backbone assignments were based on 3D HNCACB, CBCA(CO)NH, HNCA, HNHA, and HNCO experiments. Aliphatic and aromatic side chain assignments were determined using the ASCAN 68 69 All assignments were verified manually using the CARA and CCPNmr Analysis programs. 70,71 1 H chemical shifts were calibrated from internal 3-(trimethylsilyl)propane-1-sulfonic acid (DSS) and 15 N and 13 C shifts were referenced indirectly. 72 The ATNOS and CANDID 73, 74 algorithms in the J-UNIO suite were used to pick the NOESY spectra and to calculate the structure of the C domain. A globular fold obtained after the first cycle remained consistent throughout the calculation with a steady decrease in RMSD of the ensemble. A set of 132 tight dihedral angle restraints were obtained from the program Talos1. 75 A set of loose /, w, and v 1 restraints produced by the HABAS algorithm in CYANA 2.0 76 based on intraresidual and sequential NOEs provided an additional 347 dihedral angle restraints for the structure calculation. 73, 74 The set of unambiguous NOE assignments obtained in the final cycle of calculation included 1828 restraints or 24 restraints/ residue ( Table I) . The 20 structures with the lowest CYANA target function values in cycle 7 were further refined by explicit solvent minimization using the AMBER03 force field in explicit solvent (TIP3P-BOX) with a 10 Å box geometry using the webbased AMBER interface AMPS-NMR in the WeNMR portal. 77 Structure validation of the final ensemble employed the Protein Data Bank validation suite, MolMol 2K.2, and ProCheck 3.5.4 from the PSVS suite 1.5. [78] [79] [80] [81] Validation also employed input and output of the Unio calculation and agreement of the structure with NMR observables (Supporting Information Table S2 ).

The atomic coordinates of the ensemble of conformers of Figure 2 (B) have been deposited in the Protein Data Bank with accession number 5UTV. The sequence-specific resonance assignments have been deposited in the BioMagResBank with accession number 30247. 

Binding assays were conducted by incubating purified HKU9 C protein with a set of DNA and RNA oligonucleotides available in our laboratory for studying protein-nucleic acid interactions. Proteinoligonucleotide mixtures were incubated at 258C for 1 h in EMSA buffer: 20 mM sodium phosphate (pH 6.0), 75 mM NaCl, and 3% glycerol. G-quadruplex oligomers were annealed by heating to 958C for 5 min and slowly cooling to 188C overnight in buffer: 20 mM sodium phosphate (pH 6.0), 75 mM NaCl. The mixtures were resolved by native electrophoresis on 10% TBE gels (Invitrogen) for 1 h at 48C. Gels were stained with SYBR Gold stain (Invitrogen) and visualized by the Safe Imager 2.0 Blue-Light Transilluminator (Invitrogen).

Molecular pathology of emerging coronavirus infections

Coronaviruses: an overview of their replication and pathogenesis

Severe acute respiratory syndrome

SARS and MERS: Recent insights into emerging coronaviruses

Comparative analysis of twelve genomes of three novel group 2c and group 2d coronaviruses reveals unique group and subgroup features

Bats are natural reservoirs of SARS-like coronaviruses

Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor

Isolation and characterization of a novel bat coronavirus closely related to the direct progenitor of severe acute respiratory syndrome coronavirus

Further evidence for bats as the evolutionary source of Middle East respiratory syndrome coronavirus

Coronaviruses in bats from Mexico

Genetic characterization of Betacoronavirus lineage C viruses in bats reveals marked sequence divergence in the spike protein of pipistrellus bat coronavirus HKU5 in Japanese pipistrelle: implications for the origin of the novel Middle East respiratory syndrome coronavirus

Coexistence of different genotypes in the same bat and serological characterization of Rousettus bat coronavirus HKU9 belonging to a novel betacoronavirus subgroup

Genomic characterization of seven distinct bat coronaviruses in Kenya

Detection of novel SARS-like and other coronaviruses in bats from Kenya

Metagenomic analysis of viruses from bat fecal samples reveals many novel viruses in insectivorous bats in China

Putative receptor binding domain of bat-derived coronavirus HKU9 spike protein: evolution of betacoronavirus receptor binding motifs

Characterization of a critical interaction between the coronavirus nucleocapsid protein and nonstructural protein 3 of the viral replicase-transcriptase complex

An interaction between the nucleocapsid protein and a component of the replicase-transcriptase complex is crucial for the infectivity of coronavirus genomic RNA

The conserved coronavirus macrodomain promotes virulence and suppresses the innate immune response during severe acute respiratory syndrome coronavirus infection

The nsp3 macrodomain promotes virulence in mice with coronavirusinduced encephalitis

The nsp1, nsp13, and M proteins contribute to the hepatotropism of murine coronavirus JHM.WU

Coronavirus nonstructural protein 1 is a major pathogenicity factor: implications for the rational design of coronavirus vaccines

Mutagenesis of S-adenosyl-l-methionine-binding residues in coronavirus nsp14 N7-methyltransferase demonstrates differing requirements for genome translation and resistance to innate immunity

Murine coronavirus ubiquitin-like domain is important for papain-like protease stability and viral pathogenesis

Severe acute respiratory syndrome coronavirus nonstructural proteins 3, 4, and 6 induce double-membrane vesicles

Topology and membrane anchoring of the coronavirus replication complex: not all hydrophobic domains of nsp3 and nsp6 are membrane spanning

Structural basis of severe acute respiratory syndrome coronavirus ADP-ribose-1 00 -phosphate dephosphorylation by a conserved domain of nsP3

Crystal structure of the Middle East respiratory syndrome coronavirus (MERS-CoV) papain-like protease bound to ubiquitin facilitates targeted disruption of deubiquitinating activity to demonstrate its role in innate immune suppression

X-ray structural and functional studies of the three tandemly linked domains of nonstructural protein 3 (nsp3) from murine hepatitis virus reveal conserved functions

Dissection of amino-terminal functional domains of murine coronavirus nonstructural protein 3

Proteomics analysis unravels the functional repertoire of coronavirus nonstructural protein 3

The Nidoviruses: Toward Control of SARS and Other Nidovirus Diseases

Deubiquitination, a new function of the severe acute respiratory syndrome coronavirus papain-like protease?

SARS coronavirus unique domain: three-domain molecular architecture in solution and RNA binding

The molecular biology of coronaviruses

Replication of murine hepatitis virus is regulated by papain-like proteinase 1 processing of nonstructural proteins 1, 2, and 3

Nidovirus papain-like proteases: multifunctional enzymes with protease, deubiquitinating and deISGylating activities

The ADP-ribose-1 00 -monophosphatase domains of severe acute respiratory syndrome coronavirus and human coronavirus 229E mediate resistance to antiviral interferon responses

The crystal structures of Chikungunya and Venezuelan equine encephalitis virus nsP3 macro domains define a conserved adenosine binding pocket

down-regulates SARS coronavirus replication and is targeted by the SARS-unique domain and PLpro via E3 ubiquitin ligase RCHY1

Nuclear magnetic resonance structure shows that the severe acute respiratory syndrome coronavirus-unique domain contains a macrodomain fold

A G-quadruplex-binding macrodomain within the "SARS-unique domain" is essential for the activity of the SARS-coronavirus replication-transcription complex

Identification and characterization of severe acute respiratory syndrome coronavirus replicase proteins

Proteolytic processing, deubiquitinase and interferon antagonist activities of Middle East respiratory syndrome coronavirus papain-like protease

Assessing activity and inhibition of Middle East respiratory syndrome coronavirus papain-like and 3C-like proteases using luciferase-based biosensors

The J-UNIO protocol for automated protein structure determination by NMR in solution

The structure and function of frataxin

Crystal structure of Escherichia coli CyaY protein reveals a previously unidentified fold for the evolutionarily conserved frataxin family

TM-align: a protein structure alignment algorithm based on the TM-score

Dali server: conservation mapping in 3D

Structure-function analysis of Friedreich's ataxia mutants reveals determinants of frataxin binding and activation of the Fe-S assembly complex

The MotA transcription factor from bacteriophage T4 contains a novel DNA-binding domain: the 'double wing' motif

Understanding the role of dynamics in the iron sulfur cluster molecular machine

AcsD catalyzes enantioselective citrate desymmetrization in siderophore biosynthesis

Protein-ligand binding site recognition using complementary binding-specific substructure comparison and sequence profile alignment

The macro domain is an ADP-ribose binding module

Differential activities of cellular and viral macro domain proteins in binding of ADP-ribose metabolites

Identification of macrodomain proteins as novel Oacetyl-ADP-ribose deacetylases

Computational and experimental studies of mono-and poly-ADP-ribosylation of peptides

Characterization and biochemical analysis of noncanonical coronavirus macrodomains

COFACTOR: an accurate comparative algorithm for structure-based protein function annotation

Structure of a delivery protein for an AAA1 protease in complex with a peptide degradation tag

Modular assembly of RWD domains on the Mis12 complex underlies outer kinetochore organization

FFAS server: novel features and applications

Severe acute respiratory syndrome coronavirus papain-like protease: structure of a viral deubiquitinating enzyme

Robotic cloning and protein production platform of the northeast structural genomics consortium

Automated amino acid side-chain NMR assignment of proteins using (13)C-and (15)N-resolved 3D [(1)H, (1)H]-NOESY

Heteronuclear multidimensional NMR experiments for the structure determination of proteins in solution employing pulsed field gradients

Optimizing the process of nuclear magnetic resonance spectrum analysis and computer aided resonance assignment. [Optimizing the process of nuclear magnetic resonance spectrum analysis and computer aided resonance assignment

The CCPN data model for NMR spectroscopy: development of a software pipeline

) 1H, 13C and 15N chemical shift referencing in biomolecular NMR

Protein NMR structure determination with automated NOE assignment using the new software CANDID and the torsion angle dynamics algorithm DYANA

Protein NMR structure determination with automated NOEidentification in the NOESY spectra using the new software ATNOS

TALOS1: a hybrid method for predicting protein backbone torsion angles from NMR chemical shifts

Automated stereospecific proton NMR assignments and their impact on the precision of protein structure determinations in solution

A Grid-enabled web portal for NMR structure refinement with AMBER

Evaluating protein structures determined by structural genomics consortia

MOLMOL: a program for display and analysis of macromolecular structures

PROCHECK: a program to check the stereochemical quality of protein structures

MolProbity: all-atom structure validation for macromolecular crystallography

We thank Chong Tian, Pamlea N. Brady and members of the Johnson laboratory for technical assistance.

The authors declare that they have no conflicts of interest with the contents of this article.