key: cord-0850786-l2fsoedg
authors: Ni, Xiaomin; Schröder, Martin; Olieric, Vincent; Sharpe, May E.; Hernandez-Olmos, Victor; Proschak, Ewgenij; Merk, Daniel; Knapp, Stefan; Chaikuad, Apirat
title: Structural Insights into Plasticity and Discovery of Remdesivir Metabolite GS-441524 Binding in SARS-CoV-2 Macrodomain
date: 2021-03-16
journal: ACS Med Chem Lett
DOI: 10.1021/acsmedchemlett.0c00684
sha: b2840e000a2e9fc33c3a6e452a65c4b34034af8b
doc_id: 850786
cord_uid: l2fsoedg

[Image: see text] The nsP3 macrodomain is a conserved protein interaction module that plays essential regulatory roles in the host immune response by recognizing and removing posttranslational ADP-ribosylation sites during SARS-CoV-2 infection. Thus targeting this protein domain may offer a therapeutic strategy to combat current and future virus pandemics. To assist inhibitor development efforts, we report here a comprehensive set of macrodomain crystal structures complexed with diverse naturally occurring nucleotides, small molecules, and nucleotide analogues including GS-441524 and its phosphorylated analogue, active metabolites of remdesivir. The presented data strengthen our understanding of the SARS-CoV-2 macrodomain structural plasticity and provide chemical starting points for future inhibitor development.

T he recently emerging coronavirus SARS-CoV-2 has become a major global health concern. SARS-Cov-2 is an enveloped (+) single-stranded RNA ((+)ssRNA) coronavirus, encoding several functionally relevant enzymes and protein interaction domains that are necessary for the viral life cycles and deregulation of the innate immune response of the host. 1,2 Nonstructural protein 3 (nsP3) is one of the most complex coronavirus proteins comprising several domains and activities. 1,3−5 The coding sequence includes a macrodomain (MD) module, which is structurally conserved in coronaviruses and generally binds adenosine diphosphate ribose (ADPribose). 6 The viral MDs have been shown to play a pivotal role in recognizing and removing posttranslational ADP-ribosylation from specific sites. 6−9 This posttranslational modification, which may occur as poly-ADP-ribosylation or mono-ADPribosylation (MARylation), is catalyzed by poly-ADP-ribose polymerase/diphtheria toxin-like ADP-ribosyltransferase (PARP/ARTD) enzymes, and it is associated with stress signaling and the innate immune response against pathogens. 10, 11 Several MARylating PARP enzymes are activated in response to pathogen-associated molecular patterns and interferons, 9, 11 therefore suggesting that MARylation is a key modification of the innate immune response against viruses.

Several studies have provided strong evidence that viral MDs, including those of several coronaviruses, such as SARS-CoV and MERS-CoV, and different alphaviruses, such as Chikungunya virus and O'nyong'nyong virus, function as hydrolases capable of removing ADP-ribose from proteins. 6,9,12−14 Similarly, such de-mono-ADP-ribosylating activ-ity (de-MARylation) has recently been demonstrated for the first nsP3 MD, also known as the Mac1, of SARS-CoV-2, 15, 16 whereas the other two MDs likely exert distinct activities. 6 The MD and its de-MARylation function is important for efficient viral replication, as viruses harboring inactive MDs exhibit reduced replication ability and pathogenicity 10, 17 and are sensitive to interferon pretreatment. 18, 19 In addition, the de-MARylation activity of the viral MD has been convincingly linked, directly or indirectly, to the function of viral proteins involved in viral genome replication or in antagonizing the innate immune response. 9, 10 The central role of the SARS-CoV-2 MD in viral replication and in mediating inflammatory responses in the host suggests that targeting this protein module may offer a therapeutic strategy for combating severe respiratory disease. Recent studies have elucidated the crystal structure of this protein module, 15,16,20−22 demonstrating that the MD of SARS-CoV-2 exhibits a macro-H2A topology highly similar to homologues of other viruses including SARS-CoV. 9, 16, 23, 24 The pocket is lined with a set of highly conserved amino acids and catalytic residues, and thus its interaction with ADP-ribose resembles binding modes reported for other viral MD homologues. 16, 20 Nevertheless, some subtle structural differences have been observed, suggesting that SARS-CoV-2 MD may possess some degree of plasticity, 20 opening a possibility for the binding of diverse small molecules. 15 Thus to strengthen our understanding of the nature and druggability of the SARS-CoV-2 MD, we report a collective set of crystal structures, including an apo state as well as highresolution structures of complexes with diverse naturally occurring nucleosides and nucleotides. A comparison of these structures highlighted the ligand-induced domain plasticity of the ADP-ribose binding site. In addition, we describe the unprecedented discovery of the binding of GS-441524 and its phosphorylated analogue, active metabolites of remdesivir.

To understand the plasticity of the MD, we first attempted to crystallize an apo state and obtained two orthorhombic crystals that diffracted to 2.16 to 2.20 Å resolution. The presence of five and three molecules in their asymmetric units therefore suggested the potential flexibility of the protein. Indeed, a structural comparison of all chains revealed that whereas the core β-sheet and helices superimposed well, some differences in the conformations of residues lining the ADPribose binding pocket were evident ( Figure 1A ,B). By comparison, the adenosine monophosphate (AMP) binding site was more rigid than the ribose-1-phosphate binding site, in which all three loops connecting β3-α2, β5-α4, and β6-α5 that lined the pocket exhibited various degrees of flexibility. Similar to a previous report, 20 we observed diverse conformations of Phe132 and Ile131 located in the phosphate binding site ( Figure 1B ). The flexibility of the phenylalanine and isoleucine was rather interesting because it potentially determines the availability of the pocket to diverse ligands. This role was supported by conformational changes observed in the complexes with MES and HEPES molecules that were used as buffers in crystallization reagents and occupied the ribose-1phosphate binding site ( Figure 1C ,D). Both buffer molecules bound within the phosphate binding site similarly to a previous observation 20 and caused conformational changes to the loop β6-α5 residuesin essence, an out-swing of Phe132 and a rotation of the Ile131 side chain. In comparison, a slight difference was observed for the binding of MES that also induced the Ala129-Gly130 main-chain flip, enabling two hydrogen bonds that were absent in the HEPES complex ( Figure 1C ,D). Overall, the structural changes seen in the apo state and both complexes with these weakly interacting ligands suggested by their high B-factors therefore revealed the potential intrinsic plasticity of the MD, especially within its ribose-phosphate binding pocket, which may form a determinant factor for the binding of ligands.

To further understand the conformational changes necessary for the binding of ADP-ribose, we performed soaking of the ligand into the apo crystals and compared the ligand-bound form (2.50 Å resolution) with the apo state. Interestingly, a number of significant side-chain alterations were observed, primarily at two regions: loop β6-α5 as well as loop β3-α2, which constructed the ribose-phosphate binding site ( Figure  2A ). Substantial alterations were noted for the former structural motif, which included an ∼56°outswing of Phe132, vacating the phosphate binding pocket, in addition to a rotation of the Ile131 side chain to pack on top of the terminal ribose and a main-chain carbonyl flip of A129 to enable hydrogen-bond interactions between the Gly130 amine and the phosphate moiety of the ADP-ribose. At the opposing site of the binding pocket, small alterations to the glycine-rich loop β3-α2 were also evident, involving the rearrangement of Gly46−47 main chains and bringing the loop closer to the ligand. Nevertheless, two configurations of this glycine-rich loop were observed for the ADP-ribose complex: one had a Gly47 carbonyl atom flipped "in" (two of five monomers) and the other conformation had an "out"-oriented carbonyl group (in the other three monomers; Figure 2B ,C). These structural changes provided further evidence of the high flexibility of this glycine-rich loop, which contains several catalytic residues that are important for enzymatic function. 16 Nonetheless, the two conformations of the loop did not alter the contacts between the terminal ribose and the protein, which were also highly conserved compared with similar complexes that have been previously reported. 16, 20 The plasticity of the ADP-ribose binding pocket of the MD prompted us to question whether it might be able to accommodate other nucleosides and nucleotides. Thus we attempted to cocrystallize the protein with selected naturally occurring ligands, including AMP, GDP-glucose, ribose-1phosphate, β-NAD (nicotinamide adenine dinucleotide), and β-NADP (nicotinamide adenine dinucleotide phosphate), all of which (except for ribose-1-phosphate) led to the successful determination of complexed structures at 1.55 to 2.05 Å resolution. An examination of the electron density suggested, however, the modification of the bound ligands ( Figure 3 ). In the case of AMP, only the adenosine moiety, the product of autohydrolysis lacking the phosphate group, was presented in the structure. Interestingly, the adenosine ligand exerted two binding modes: One resembled that of the adenosine moiety in the ADP-ribose complex, whereas a slight tilt of the adenine ring in the second binding mode positioned the ribose group outside in the solvent-exposed region ( Figure 3A,B) . Although the hydrogen bond between the adenine amine group and Asp22 was still maintained, different water-mediated contacts between the ribose moiety and the protein were evident.

For GDP-glucose, the electron density allowed modeling of only guanosine monophosphate (GMP), a product of hydrolysis. Nonetheless, the binding of GMP in the ribose-1phosphate binding pocket of the MD was intriguing ( Figure  3C ). The guanine ring surprisingly occupied the cavity reserved for the terminal ribose moiety of ADP-ribose, not the AMP-binding site, likely due to a poor fit between the guanine 6-hydroxyl and 2-amine groups and Asp22 and the backbone amide of Val155 and Phe156 of loop β7-α6. The bound GMP was stabilized mainly through interactions of the guanine ring, which was sandwiched between Ile131 and the glycine-rich loop β3-α2, forming a network of hydrogen bonds to Lys44, Val49, Ala50, and Ala38. No contribution to binding was observed from the ribose and phosphate groups that were positioned outside the binding pocket in the solvent region.

For β-NAD and β-NADP, the observed electron density suggested that the entities of the bound ligands were likely nicotinamide-cleaved products, resulting in the complex of ADP-ribose and ADP-ribose-2′-phosphate (ADPRP), respectively. Although hydrolase activity catalyzing NAD hydrolysis 25 has been documented for some MDs, we assumed that these cleavage products were likely a consequence of autohydrolysis of the compounds because no significant increase in ADPribose or ADPRP traces was observed in a high-performance liquid chromatography (HPLC) assay upon an incubation of β-NAD or β-NADP with SARS-CoV-2 MD (data not shown). A structural comparison demonstrated that both the NADcleaved product ADP-ribose and the NADP-cleaved product ADPRP assumed a binding mode that was highly identical to that of ADP-ribose (Figures 2 and 3D ). Nevertheless, a difference was noted for ADPRP in that the 2′-phosphate group, which was located in the solvent region, could engage additional magnesium-mediated contacts with an inward side chain of Asp157 ( Figure 3D ).

The flexibility of the ADP-ribose binding pocket with its ability to accept various ligands prompted us to speculate the potential binding of antiviral nucleoside analogues. We therefore selected a set of diverse clinically used antivirals, including abacavir, entecavir, acyclovir, gemcitabine, remdesivir, and remdesivir metabolite GS-441524, and tested the binding of these drugs using cocrystallization. Among the set, an examination of electron density maps revealed that the remdesivir metabolite GS-441524 was the only ligand that showed binding in the crystal structures (2.15 Å resolution; Figure 4A ,B). The interaction of this compound was rather unexpected because this metabolite and its active triphosphorylated form have been designed to target the RNAdependent RNA polymerase (RdRp) of several viruses, including Ebola virus, MERS-CoV, SARS-CoV, and potentially SARS-CoV-2. 26−29 In addition, a recent study has further demonstrated that the direct use of GS-441524 can potently inhibit the replication of SAR-CoV-2 in a mouse model. 30 On the basis of our crystal structure, a structural comparison demonstrated that the binding mode of GS-441524 in the MD highly resembled that of the adenosine moiety of ADP-ribose. The amine group of the pyrrolotriazine ring, which was sandwiched between Val49 and Phe159, formed a hydrogen bond to Asp22, whereas the hydroxyl groups of the ribose moiety engaged two contacts with the backbone carbonyl atoms of Leu126 and Val155 ( Figure 4C ). The 1′-nitrile group was accommodated in the space adjacent to loop β7-α6. With a distance of ∼3.3 Å, direct contacts to the backbone amine atoms of Phe156 and Asp157 might be weak, yet charge compatibility with this binding pocket was provided by the negative electrostatic potential of the nitrile group.

In comparison with the adenosine substrate, GS-441524 showed improved shape complementarity with the binding pocket ( Figures 3A and 4C) . To assess this, we performed isothermal titration calorimetry (ITC) to determine the affinities of this ligand, adenosine, AMP, and ADP-ribose.

Unfortunately, adenosine and AMP did not yield interpretable ITC binding isotherms. Nonetheless, we observed that the K D values of GS-441524 and ADP-ribose were remarkably comparable (10.8 and 7.3 μM, respectively). This was a surprising result considering that the two phosphate moieties and the terminal ribose are missing in GS-441524. Thus we synthesized the 5′-monophosphate derivative of this compound (GS-441524 monophosphate) and characterized the binding. Interestingly, an installation of 5′-phosphate led to a slightly improved affinity with a similar binding strength compared with ADP-ribose (K D of 8.6 μM), suggesting a contribution of the phosphate group for binding. Nevertheless, the different thermodynamics of the three ligands were of particular note. The presence of the phosphate group in GS-441524 monophosphate and ADP-ribose as well as the additional terminal ribose-1′-phosphate moiety in the latter resulted in large negative binding enthalpy changes, which were counteracted by highly unfavorable entropy changes (TΔS) that were likely attributable to the constrained ligand and potentially the pocket in the bound state.

In conclusion, targeting the MD offers an attractive target for the development of antiviral agents against SARS-CoV-2 and other viruses. To assist the inhibitor discovery efforts, we report a collective set of crystal structures, extending our structural knowledge of this protein. A structural comparison demonstrated that the ADP-ribose binding site, essentially the pocket for ribose-1-phosphate, of the MD possessed high structural plasticity. This was demonstrated by its adaptability to diverse naturally occurring nucleosides and nucleotides, including ADP-ribose-phosphate (ADPRP). This flexibility offers opportunities for the rational targeting of this protein module by small-molecule inhibitors. Indeed, a recent study has reported a number of potential small-molecule binders discovered through crystallographic fragment screening. 31 In line with this, our unprecedented discovery of the binding of GS-441524, an active metabolite of remdesivir, supported this hypothesis that the ADP-ribose pocket of the MD indeed represents a druggable site. The low micromolar affinity of GS-441524 and its small molecular weight offer a versatile starting point for ligand design.

Protein Production. The DNA encoding the SARS-CoV-2 MD was commercially synthesized and subcloned into a pET-28a(+) vector (Genscript, Supplementary Table S1 ). The expression of the recombinant protein harboring N-terminal His 6 tags was carried out in E. coli Rosetta, which was cultured in TB media. The culture was grown at 37°C until an OD 600 of 1.6 was reached prior to cooling to 18°C. At an OD 600 of 2.6 to 2.8, the protein expression was induced by adding 0.5 μM IPTG, and the expression was continued overnight. Cells were harvested by centrifugation and resuspended in buffer containing 50 mM Tris, pH 8.0, 500 mM NaCl, 10 mM imidazole pH 7, 5% glycerol, and 1 mM TCEP. Lysis was performed by sonication, and the supernatant was clarified by centrifugation. The recombinant protein was initially purified using Ni 2+ -affinity chromatography, and subsequent cleavage of the histidine tag was performed by TEV protease treatment. The cleaved protein was passed through a Ni 2+ -NTA column and was further purified by size exclusion chromatography using a Superdex S75 column. The pure protein was stored at −80°C in 25 mM Tris pH 8.0, 150 mM NaCl, and 5% glycerol buffer.

Crystallization and Structure Determination. The MD was buffer-exchanged into 25 mM Tris pH 8.0 and 150 mM NaCl and concentrated to ∼8.5−10 mg/mL. Crystallization was performed using the sitting drop vapor diffusion method at 20°C, and the conditions are listed in Supplementary Table S2 . Microseeding was employed in all crystallization experiments using the initial crystals of the MD, which were stored in 30% broad-molecular-weight PEG smears, 32 0.1 M MgCl 2 , 0.1 M tris pH 7.0. For cocrystallization, the protein was mixed with a 10-fold molar excess of the ligands. In the case of ADP-ribose, soaking was performed overnight using the ligand at 10 mM concentration. Viable crystals were cryoprotected with mother liquor supplemented with 20% ethylene glycol prior to flash cooling in liquid nitrogen.

Diffraction data were collected at the Swiss Light Source and were processed and scaled with XDS 33 and aimless, 34 respectively. The initial structure solution was obtained with the molecular replacement method using Phaser 35 and the coordinate of the MD (PDB ID: 6WEN). Manual model rebuilding alternated with refinement was performed in COOT 36 and REFMAC5. 37 The geometry of the final models was verified using MOLPROBITY. 38 Omitted electron density maps for the bound ligands are shown in Supplementary Figure S1 . The summary of the data collection and refinement statistics are in Supplementary Table S2 .

Synthesis of GS-441524 Monophosphate. The chemical synthetic procedure is described in the Supporting Information.

Isothermal Calorimetry. The ITC experiment was performed using a NanoITC instrument (TA Instrument) at 25°C in the buffer containing 20 mM Tris pH 8.0, 100 mM NaCl, 0.5 mM TCEP, and 5% glycerol. For ADP-ribose and GS-441524 monophosphate, the ligand at 500 μM in the syringe was titrated into the reaction cell containing the protein at 35−106 μM, whereas for GS-441524, the protein at 250 μM was titrated into the reaction cell containing the compound at 22 μM. Data analyses were performed with Nano-Analyze software (TA Instrument) using an independent binding model from which the affinities and thermodynamic parameters were calculated.

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acsmedchemlett.0c00684.

Supplementary figure and table for the electron density map for the bound ligands, protein sequence, crystallographic data collection, and refinement statistics and chemical synthesis (PDF)

A new coronavirus associated with human respiratory disease in China

The species Severe acute respiratory syndrome-related coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2

Genomic characterization of the 2019 novel humanpathogenic coronavirus isolated from a patient with atypical pneumonia after visiting Wuhan

The Chinese SARS Molecular Epidemiology Consortium. Molecular evolution of the SARS coronavirus during the course of the SARS epidemic in China

Nsp3 of coronaviruses: Structures and functions of a large multi-domain protein

The macro domain is an ADP-ribose binding module

Macrodomains: Structure, Function, Evolution, and Catalytic Activities

The Viral Macrodomain Counters Host Antiviral ADP-Ribosylation

The impact of PARPs and ADP-ribosylation on inflammation and host-pathogen interactions

ADP-Ribosylation, a Multifaceted Posttranslational Modification Involved in the Control of Cell Physiology in Health and Disease

Viral Macro Domains Reverse Protein ADP-Ribosylation

The conserved macrodomains of the non-structural proteins of Chikungunya virus and other pathogenic positive strand RNA viruses function as mono-ADP-ribosylhydrolases

Ahel, I. Reversible ADP-ribosylation of RNA

Ahel, I. Viral macrodomains: a structural and evolutionary assessment of the pharmacological potential

Biophysical, and Biochemical Elucidation of the SARS-CoV-2 Nonstructural Protein 3 Macro Domain

ADP-ribose-1″-monophosphatase: a conserved coronavirus enzyme that is dispensable for viral replication in tissue culture

Viral Macrodomains: Unique Mediators of Viral Replication and Pathogenesis

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

Crystal structures of SARS-CoV

The SARS-CoV-2 Conserved Macrodomain Is a Mono-ADP-Ribosylhydrolase

Molecular Basis for ADP-Ribose Binding to the Mac1 Domain of SARS-CoV-2 nsp3

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

Evolutionary and structural analysis of SARS-CoV-2 specific evasion of host immunity

The ARH and Macrodomain Families of alpha-ADP-riboseacceptor Hydrolases Catalyze alpha-NAD(+) Hydrolysis

A Review of Its Discovery and Development Leading to Emergency Use Authorization for Treatment of COVID-19

Remdesivir for the Treatment of Covid-19 -Final Report

Comparative therapeutic efficacy of remdesivir and combination lopinavir, ritonavir, and interferon beta against MERS-CoV

Remdesivir Metabolite GS-441524 Effectively Inhibits SARS-CoV-2 Infection in Mouse Models

Defined PEG smears as an alternative approach to enhance the search for crystallization conditions and crystal-quality improvement in reduced screens

How good are my data and what is the resolution?

Phaser crystallographic software

Features and development of Coot

Overview of refinement procedures within REFMAC5: utilizing data from different sources

MolProbity: More and better reference data for improved all-atom structure validation