Aikaterini C Tsika

Efstathios Melekis

Sofia-Antigoni Tsatsouli

Nicolas Papageorgiou

Maria Mate

Bruno Canard

Bruno Coutard
email: bruno.coutard@univ-amu.fr

Detlef Bentrop

Georgios A Spyroulias
email: g.a.spyroulias@upatras.gr@url

Maria J Maté

Prof B Coutard

Prof G A Spyroulias

Deciphering the Nucleotide and RNA Binding Selectivity of the Mayaro Virus Macro Domain

Keywords: Mayaro virus, Alphavirus, macro domain, ADP-ribose, RNA, NMR MAYV, Mayaro virus, CHIKV, Chikungunya virus, VEEV, Venezuelan Equine Encephalitis virus, SINV, Sindbis virus, MD, macro domain, ss(+) RNA, positive-sense single stranded RNA, PTM, Post Translational Modification, ADP-r, Adenosine diphosphate ribose, ADP, Adenosine diphosphate, ATP, Adenosine triphosphate, AMP, Adenosine monophosphate, NAD, Nicotinamide adenine dinucleotide, NMR, Nuclear Magnetic Resonance, HSQC, Heteronuclear Single Quantum Coherence, NOESY, Nuclear Overhauser Effect Spectroscopy, CSP, Chemical Shift Perturbation, ITC Isothermal Titration Calorimetry, EMSA, Electrophoretic Mobility Shift Assay

à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

Introduction

Alphaviruses belong to the Togaviridae family, whose members form enveloped particles characterized by icosahedral symmetry enclosing their genome. They cause mainly mosquito-borne diseases in vertebrates such as birds and rodents, which act as reservoir hosts. The genus Alphavirus is divided in two subcategories depending on the viruses' primary areas of endemicity and clinical manifestations 1 2 . Old World Alphaviruses were initially identified in Asia and Africa, where they have been considered as causal agent of arthropathies 3 4 . By contrast, New World Alphaviruses are usually implicated in neurological diseases with complications such as dysarthria, motor disorders and abnormal reflexes, and were detected solely in the Americas [START_REF] Forrester | Genome-scale phylogeny of the alphavirus genus suggests a marine origin[END_REF] . Mayaro virus (MAYV) is an exception to this classification. Despite its endemic distribution in South America, the induced disease has a clinical picture similar to the Old World Alphaviruses [START_REF] Tesh | Mayaro virus disease: an emerging mosquitoborne zoonosis in tropical South America[END_REF] . In particular, symptoms include mildto-severe Dengue-like febrile syndrome, characterized by fever, headache, rash, myalgia, arthralgia and arthritis, which can be debilitating and can persist for months 4 7 8 . These clinical manifestations are reported also in the Chikungunya virus (CHIKV) disease.

The resemblance of MAYV and CHIKV diseases, the worldwide outbreaks of CHIKV, and the ability of MAYV to be reproduced into mosquitoes inhabiting urban environments underline the importance of this "neglected virus" 6 9 10 11 . It is worth mentioning that MAYV is mainly transmitted by mosquitos of Haemagogus and Aedes genera, the latter being also responsible for the distribution of Zika, Dengue and Chikungunya viruses. The recent Alphavirus outbreaks in concert with the absence of vaccines and targeted antiviral therapy emphasize the necessity to understand their replication mechanisms in detail 3 12 . The global distribution and the constantly reported outbreaks of these viral infections point out the significance of preventing the potentially chronic, disabling and long-lasting complications caused by these viruses, in order to reduce the social and economic impacts of these zoonoses.

The replication of Alphavirus occurs in spherules 13 14 containing four non-structural proteins (nsP1 to nsP4) [START_REF] Hellström | Partially uncleaved alphavirus replicase forms spherule structures in the presence and absence of RNA template[END_REF] resulting from processing of the precursors P123 and P1234. NsP1 is carrying the enzymatic activities for the mRNA cap formation, nsP2 has a N-terminal domain with multiple enzymatic activities (e.g. ATPase, RNA helicase, RNA triphosphatase) 16 17 18 . Its C-terminal part is carrying a protease activity for polyprotein processing with a papain-like domain and is followed by a classical 2'O methyltransferase fold 19 20 . NsP4 is carrying the RNA dependent RNA polymerase activity [START_REF] Chen | Chikungunya virus nsP4 RNA-dependent RNA polymerase core domain displays detergent-sensitive primer extension and terminal adenylyltransferase activities[END_REF] in charge of the replication/transcription. Many functions have been attributed to nsP3 in recent years [START_REF] Götte | The Enigmatic Alphavirus Non-Structural Protein 3 (nsP3) Revealing Its Secrets at Last[END_REF] . Its C-terminal region, named HVR (hypervariable region), is a hub for interactions with several partners while its N-terminal region contains a globular macro domain (MD).

MDs are conserved and widely distributed versatile structural motifs of 130-190 amino acids. They are found both in prokaryotes and eukaryotes as well as in viruses [START_REF] Feijs | Macrodomain-containing proteins: regulating new intracellular functions of mono (ADP-ribosyl) ation[END_REF] [START_REF] Han | The macro domain protein family: structure, functions, and their potential therapeutic implications[END_REF] 25 26 27 .

Only few genera of positive-sense single stranded RNA ((+) ss RNA) viruses contain MD folds, e.g. Alphavirus, Coronovirus (CoV), Rubella virus and Hepatitis E virus. The Alphaviruses MDs are located at the N-terminal region of nsP3 and adopt a β/α/β mixed sandwich topology.

Initially, MD's biological functions were linked to their ability to bind ADP-ribose (ADP-r), its derivative (Poly-) ADP-r and RNA. In addition, Alphaviruses MDs were found to be adenosine di-phosphoribose 1''-phosphate phosphatases, but the role of this activity in viral replication remained elusive [START_REF] Malet | The crystal structures of Chikungunya and Venezuelan equine encephalitis virus nsP3 macro domains define a conserved adenosine binding pocket[END_REF] . Subsequently, several studies proved its significance to viral replication and virulence, while recent studies proposed that the function of MDs is associated with their activity as ADP-ribosylhydrolases, suggesting their implication in counteracting responses to infection 28 29 30 31 32 33 34 35 36 . In general, viral MDs are postulated to be critically involved in replication complex formation, cell specific viral replication and ADPribosylation dependent interactions with host proteins [START_REF] Abraham | ADP-ribosyl-binding and hydrolase activities of the alphavirus nsP3 macrodomain are critical for initiation of virus replication[END_REF] . It remains to be investigated if and how distinct properties of different Alphaviruses MDs such as sequence, structure, dynamics and ligand binding can be correlated to different functions and different viral replication profiles. ADP-ribosylation is a fundamental nucleotide-based post-translational modification (PTM) that is vital to cellular processes including cell signaling, DNA repair, RNA regulation and probably many others in various organisms. The addition as well as the removal of an ADP-r moiety to or from cellular targets are strongly regulated processes leading to functional changes through mono-ADP-ribosylation (MARylation) and/or poly-ADP-ribosylation (PARylation). The induced alterations of enzymatic activities and the ability of ADPribosylated proteins to serve as "anchors" for biomolecular interactions are connected with their assigned biological roles. These include DNA damage response, transcription signaling, cellular stress response, protein degradation and immune response. Most of the aforementioned functions require the cross-talk of proteins participating in the ADPribosylation and in other PTMs, such as ubiquitination 38 39 40 39 41 42 . Recently, it was reported that ADP-ribosylation is induced upon Alphavirus infection in an interferonindependent manner [START_REF] Abraham | ADP-ribosyl-binding and hydrolase activities of the alphavirus nsP3 macrodomain are critical for initiation of virus replication[END_REF] .

The enzymes and protein domains involved in ADP-ribosylation are characterized as writers, erasers and readers 42; 43; 44; 45 based on their ability to transfer, remove and recognize ADP-r or ADP-ribosylated substrates. Among proteins binding ADP-r, MDs are classified as readers or/and erasers [START_REF] Luscher | ADP-Ribosylation, a multifaceted posttranslational modification involved in the control of cell physiology in health and disease[END_REF] .

Herein, we address the structure-function relationship of the MAYV MD through the determination of its high-resolution NMR structure and dynamics. Furthermore, interaction studies with ADP-r, other adenosine derivatives and RNA oligonucleotides were performed through heteronuclear NMR spectroscopy, ITC and EMSA analysis. The characterization on both a structural and a biochemical level elucidates physicochemical features of the MAYV MD, reveals that non-conserved loops are crucial for the structural plasticity and substrate selectivity of the domain, and provides novel mechanistic insights into the binding of various adenosine-based ligands.

Results & Discussion

NMR solution structure of MAYV MD

The solution structure of MAYV MD was determined from more than 1600 NMR-derived geometrical restraints (PDB ID: 5IQ5; BMRB accession no. 30043) (Figure 1a and Table S2) [START_REF] Melekis | NMR study of non-structural proteins-part I: 1H, 13C, 15N backbone and side-chain resonance assignment of macro domain from Mayaro virus (MAYV)[END_REF] ID: 3GQE Z score: 24, RMSD: 1.7 Å), all of them exhibit significant topological similarities with MAYV MD. Structural differences between these proteins are mainly localized in four loop regions (β2-α1, α1-β3, β5-α3, β6-α4) (Figure 1c), which are believed to be responsible for the catalytic properties [START_REF] Zapata-Pérez | Structural and functional analysis of Oceanobacillus iheyensis macrodomain reveals a network of waters involved in substrate binding and catalysis[END_REF] and the different affinities of MDs for ADP-r [START_REF] Mcpherson | ADP-ribosylhydrolase activity of Chikungunya virus macrodomain is critical for virus replication and virulence[END_REF] . A striking example is the loop β2-α1 with only two conserved glycine residues, which is considered crucial for the catalysis [START_REF] Zapata-Pérez | Structural and functional analysis of Oceanobacillus iheyensis macrodomain reveals a network of waters involved in substrate binding and catalysis[END_REF] The other three loops are important for ligand binding and exhibit ~30% sequence identity (Figure 1c). Loop α1-β3 contains only a single conserved residue at position 45 (Phe) and exhibits different conformations in the four MDs. The only conserved residues of the loop β5-α3 are at its C-terminal end (Lys118, Asp119). The MDs of MAYV and CHIKV differ in three positions of this loop: 113 (Ile vs. Val) 114 (Phe vs. Tyr), and 116 (Ala vs. Gly). Finally, the β6-α4 loop is the one with the highest number of conserved residues (four out of six residues) but plays only a minor role in ligand binding and catalysis.

Analysis of the electrostatic surface potential of MAYV MD identified a solvent accessible cleft surrounded by positively charged surface areas. Helices α1 and α4 and the loops connecting the secondary structure elements β1-β2-α1-β3 and β6-α4 contribute to the formation of this cleft that could serve as a binding site for negatively charged molecules such as RNA or (Poly-) ADP-r. The opposite side of MAYV MD is negatively charged (Figure 1d).

Dynamic properties of MAYV MD

The dynamic properties of the polypeptide on the ps/ns timescale were derived from [START_REF] Hellström | Partially uncleaved alphavirus replicase forms spherule structures in the presence and absence of RNA template[END_REF] N R1 and R2 relaxation rates and the heteronuclear { 1 H N }- [START_REF] Hellström | Partially uncleaved alphavirus replicase forms spherule structures in the presence and absence of RNA template[END_REF] N NOEs (Figure S1) (average values of 1.53 ± 0.19 s -1 , 14.14 ± 2.52 s -1 and 0.84 ± 0.04, respectively, for residues 2-159). The correlation time for isotropic tumbling in solution as calculated from the R2/R1 ratio is 8.8 ns which roughly corresponds to MW ~17.6 kDa (theoretical MW 18.1 kDa) and clearly indicates that the protein is a monomer in solution. Calculation of the S 2 order parameter (Figure 1e) showed that the N-and C-terminal regions of MAYV MD are as rigid as most of the protein. Increased mobility on the ps/ns timescale was only found for three segments which comprise the loops β2-α1 (Gly32-Val33), β4-α2 (Ser77-Ala79), and β5-α3 (Ser115-Gly117). Note that signals of the neighboring residues Gly30-Asp31, Cys34-Ala36, Ala38, Gly112-Phe114, and Asp119-Arg120 are broadened beyond detection in the 1 H, [START_REF] Hellström | Partially uncleaved alphavirus replicase forms spherule structures in the presence and absence of RNA template[END_REF] N -HSQC spectrum of the apo form of MAYV MD. Thus, their 1 H, [START_REF] Hellström | Partially uncleaved alphavirus replicase forms spherule structures in the presence and absence of RNA template[END_REF] N correlation peaks could not be assigned, strongly suggesting that the loops upstream to helices α1 and α3 undergo conformational exchange processes that may be relevant for the selectivity of substrate binding [START_REF] Ziarek | Binding site identification and structure determination of protein-ligand complexes by NMR[END_REF] . These data support the plasticity of two functionally important loops (Figure 1c).

According to our NMR data, the MAYV MD loop β2-α1 seems to be more dynamic (lower S 2 values and negative values { 1 HN}- [START_REF] Hellström | Partially uncleaved alphavirus replicase forms spherule structures in the presence and absence of RNA template[END_REF] N NOEs while at the same region CHIKV exhibits only positive values; data not shown) than the corresponding segment of CHIKV MD (Figure S2).

In addition, the fact that the HN resonances of the neighboring residues mentioned above (Cys34-Ala36 and Ala38) remain unidentified while the HNs of the corresponding segment in CHIKV apo MD could be assigned 53 (Figure S2c), shows that an exchange process on an intermediate NMR time scale is present in the β2-α1 loop of MAYV MD, but not in CHIKV MD.

An H/D exchange experiment in which a [START_REF] Hellström | Partially uncleaved alphavirus replicase forms spherule structures in the presence and absence of RNA template[END_REF] N labeled sample was lyophilized and dissolved in 100% D2O provided additional information on the mobility of MAYV MD on the time scale of minutes to hours. All loops, helix α1 and the outer part of helix α4 as well as the first residues of strands β1 and β3 show complete exchange of the backbone amide protons in the first 11 minutes of the experiment (Figure 1f) [START_REF] Kleckner | An introduction to NMR-based approaches for measuring protein dynamics[END_REF] . After an exchange time of 1125 minutes the HN resonances of residues that belong to the solvent-inaccessible core of the protein which is formed by specific parts of i) strands β2 (Ala18-Val20), β4 (Ile65-Ile66, Ala68), β5 (Val104-Ile106), and β6 (Thr139-Tyr142) and ii) helices α2 (Leu85-Ile96), α3 (Leu128-Met132), and α4 (Leu156-Gln157), are still visible (Figure 1f). Hence, there is a stable hydrogen bond network in this rigid core region of MAYV MD, in contrast to the mobility of the majority of residues on the min/h time scale (Figure 1f).

MAYV MD nucleotide selectivity and binding

A series of adenosine-based molecules was used to explore the ability of MAYV MD to bind adenosine derivatives as well as to elucidate the impact of individual chemical moieties on the MAYV MD -ligand interaction and the role of certain residues in ligand recognition. Titrations of MAYV MD with ADP-r, ADP, ATP, AMP and NAD+ were monitored by acquisition of 1 H, [START_REF] Hellström | Partially uncleaved alphavirus replicase forms spherule structures in the presence and absence of RNA template[END_REF] N-HSQC spectra and showed chemical exchange between free and ligand-bound MD ranging from fast to slow on the NMR time scale. NMR titration of the MAYV MD with ADP-r suggests a strong interaction of both molecules and complex formation at 1:1 molar ratio (Figure 2a). Furthermore, the analysis of the backbone amide chemical shift perturbations (CSPs) identifies the residues that are involved in this process (Figure 2c andd). To this end a complete backbone and side chain resonance assignment of the MAYV MD-ADP-r complex was obtained (BMRB accession no. 27612).

Specifically, residues Asp10, Ala23, Asn24, Gln28, Gly32, Val33, Lys56, Val69, Ser110 and Ala116 exhibited the largest CSPs while undergoing chemical exchange between the free and the bound form in the slow exchange regime (~0.1-1 s, kex <<|Δν| and usually Kd ≈ μΜ or lower [START_REF] Ziarek | Binding site identification and structure determination of protein-ligand complexes by NMR[END_REF] ). Additionally, adjacent amino acids like Ala9, Ala12, Thr13, Gly27, Val29, Asn72, Phe73, Leu109, Ser115 and G117 were also affected exhibiting a slow-exchange process during the interaction. The backbone HN resonances of Ile11, His25, Val37, Phe45, Ala68, and Lys146 broadened beyond detection while those of Ala22, Gly30, Asp31, Cys34, Arg35, Ala36, Ala38, Gly112, Cys143, and Arg144, which were unassigned in the spectra of the apo MAYV MD, could be assigned in the 1 H, [START_REF] Hellström | Partially uncleaved alphavirus replicase forms spherule structures in the presence and absence of RNA template[END_REF] N HSQC spectrum of the MAYV MD-ADP-r complex.

These data suggest that the majority of the residues, which participate in ligand recognition and binding are located near the positively charged ligand-binding cleft (Figure 1d and2d).

It is worth mentioning that upon binding of ADP-r the backbone and side chain resonances of residues in the loops that comprise the positively charged cavity are observable, while they were not detectable in the free state of the protein. All these findings prove the high affinity interaction of MAYV MD with ADP-r and indicate ligand binding via loops β1-β2, β2-α1, β4-α2, and β5-α3, which should comprise or be located nearby the anticipated binding cleft (Figure 2d) [START_REF] Ziarek | Binding site identification and structure determination of protein-ligand complexes by NMR[END_REF] . Most residues in these loops show relatively low order parameters (Figure 1e) or conformational averaging in the apo state.

Then NMR data showed that the complex formation with ADP required a 1:3 (protein:ligand) molar ratio (no further chemical shift changes observed at higher ligand concentrations up to a 1:4 molar ratio) and that the exchange occurred in the intermediate time regime (kex ~|Δν|), strongly suggesting that ADP interacts more weakly than ADP-r.

Similar to ADP-r the addition of ADP to MAYV MD induced the largest chemical shift perturbations in the loops connecting the β1 and the β2 strands (Ala9-Thr13), the β2 strand and the α1 helix (Ala23-His25), the β4 strand and the α2 helix (Asn72, Phe73), and the β5 strand and the α3 helix (Leu108 -Thr111) (Figure 2e). During the titration with ADP the HN resonances of Ala9, Asp10, Ala12, Thr13 (Figure 2e ATP was subsequently tested as a putative ligand for MAYV MD to explore the role of the negatively charged phosphate groups. The complex formation is achieved at a molar ratio of 1:4 (protein:ligand) and led to smaller CSP values than the intermediate or slow chemical exchange during the titrations with ADP and ADP-r, respectively. However, the residues exhibiting CSPs are practically the same as in the case of ADP. Similarly, the 1 H, [START_REF] Hellström | Partially uncleaved alphavirus replicase forms spherule structures in the presence and absence of RNA template[END_REF] N HSQC cross peaks of Ile11, Ala23, Asn24, Gly32, Val33, Val69, Leu108, Ser110, and Ala116 remained broadened beyond detection even after the last addition of ATP (Figure 2f). On the other hand, Asp10, a key residue for the ligand binding, is in the fast exchange regime upon addition of ATP to MAYV MD. The same behavior was observed for the other key residues Thr13 (Figure 2f) and Leu109. Therefore, the interaction can be characterized best as fast exchange on the NMR time scale (kex >>|Δν|). Such an analysis of the binding kinetics can also provide an estimate of the binding affinity 52 52 55 , which in this case implies a weak binding with a Kd of the order of ~10 -3 M. Therefore, it may be concluded that ATP binds less tightly than ADP. Titration of MAYV MD with AMP, the simplest adenosine derivative used in this study, was characterized by relatively small CSPs and fast exchange, while the affected regions are the same as with the other ligands. (Figure 2g). As with ATP, Asp10 in the β1-β2 loop is in fast exchange and shows a prominent chemical shift change. Note that Ala23 in the β2-α1 loop whose HN resonance was broadened beyond detection during the ATP titration was observed up to the last stage of the titration with AMP. Moreover, the CSP of Val37 in the first turn of helix α1 could be followed at all titration steps for the first time. Ala116 in the β5-α3 loop showed slow exchange with ADP-r, intermediate exchange in the titrations with ADP and ATP, but fast exchange in the case of AMP. These observations exemplify the adaptability of MAYV MD to adenosine-based ligands. Titration of MAYV MD with NAD+ up to a molar ratio of 1:4 (protein:NAD+) had only a minor influence on the 1 H, [START_REF] Hellström | Partially uncleaved alphavirus replicase forms spherule structures in the presence and absence of RNA template[END_REF] N resonances of the protein indicating that NAD+ exhibits the weakest interaction among the five nucleotides tested. The intensities of the cross peaks of Ala23, Asn24, Gly32, and Ser110 are decreased or broadened beyond detection at the final step of titration. These observations indicate that despite the weak protein-NAD+ interaction these residues are affected more by the addition of NAD+ than by AMP. The segments influenced by NAD+ are the β1-β2 loop, the β5-α3 linker, the C-terminal end of the β6-α4 loop, and the first α4 residues (Figure 2h). Although NAD+ interacts only weakly with the protein, it is the only ligand that induces CSPs above threshold (lower than the threshold of any of the other four nucleotides) for the backbone HN resonances of Gln43 and Ala44 at the C-terminus of helix α1, suggesting a different interaction interface for the relatively bulky NAD+ (Figure S3).

The screening of adenosine-based molecules for their potential binding to MAYV MD revealed that ADP-r exhibits the strongest interaction and forms a 1:1 complex. In addition, the role of individual chemical moieties of the ADP-r molecule can be evaluated through a comparative analysis of the ligand-induced chemical shift changes (Figure 2 andS3). The ligands used (ADP-r, ADP, ATP, AMP, and NAD+) differ in each case in the chemical group connected to adenosine, thus allowing to elucidate the contribution of each moiety to the chemical shift changes. Absence of the terminal ribose moiety has a significant effect on the interaction, since ADP-r is the only MAYV MD ligand leading to a slow exchange process on the NMR time scale (kex << |Δν|). In general, ADP and AMP interact with the same residues of MAYV MD, but interestingly, the residues that interact with the β-phosphate group of ADP can be identified clearly. According to the literature [START_REF] Zapata-Pérez | Structural and functional analysis of Oceanobacillus iheyensis macrodomain reveals a network of waters involved in substrate binding and catalysis[END_REF] there are two structural water molecules which are H-bonded to the β-phosphate. In the crystal structures of the macro domains from Oceanobacillus iheyensis (PDB ID: 5L9K) and CHIKV (PDB ID: 3GPO) in complexes with ADP-r these water molecules interact with the conserved residues Ala29, Gly91 and Ala23, Gly70, respectively. Ala23 and Gly70 in MAYV MD exhibit much smaller CSPs (Figure 2g ande) and undergo a fast exchange process upon the addition of AMP compared to ADP, indicating an active role in the binding of the second phosphate group.

AMP and especially NAD+ exhibit weak interaction. The addition of a substituted pyridine ring to the distal ribose of ADP-r seems to impede the access of NAD+ to the binding cavity of the smaller adenosine derivatives.

Identification of MAYV MD residues critical for ADP-r binding through NMR and ITC

To further characterize the MAYV MD-ADP-r interaction, MD mutants were designed considering the most affected regions in the CSP analysis (Figure 2c) as well as available crystal structures of MD -ADP-r complexes 25 27 28 50 56 57 from various organisms. These structures reveal conserved patterns of amino acids that mediate the interaction of ADP-r with the protein (Figure 2b). According to these models, the NH2 group of adenine forms a hydrogen bond with the carboxyl group of an aspartate (Figure 2b, indicated with "*") at the beginning of the β1-β2 loop. The hydrophobic residue (Ile/Val) following this aspartate is also hydrogen-bonded to adenine via its backbone amide group. In several homologous structures (e.g. PDB ID: 2BFQ), an aromatic residue (Tyr/Phe) (indicated with "$" in Figure 2b) is found to stabilize further the adenine through π-π stacking to its purine ring. The phosphate groups are connected through their oxygen atoms to the backbone amides of residues in two conserved sequence motifs around almost invariable Gly residues. The first motif is a Gly followed by Val/Leu just before helix α1 (Figure 2b, indicated with "+") and the second one is Ser/Gly -X/Thr -Gly -Ile/Val/Asn -Tyr/Phe/Leu/Ala in the β5-α3 loop (Figure 2b, indicated with "#"). Finally, the distal ribose is stabilized by a hydrogen bond to an Asn residue in the β2-α1 loop (Figure 2b, indicated with "˄") [START_REF] Haikarainen | Proximal ADP-ribose hydrolysis in trypanosomatids is catalyzed by a macrodomain[END_REF] 57 58 59 . These residues, namely Asp10/Ile11, Gly32/Val33, Gly112, and Asn24, are also present in MAYV MD (Figure 2b) and interact with ADP-r as clearly shown by our CSP data (Figure 2c and Figure S3a). Surprisingly, the weakly conserved residues Asn72 and Phe73, in the β4-α2 loop, are also strongly affected by the binding of ADP-r indicating a rearrangement of this loop.

The impact of residues Asp10, Asn24, and Phe114 on the binding of ADP-r to MAYV MD was examined by means of three mutants: D10A, N24A/F114A, and the triple mutant D10A/N24A/F114A. Titrations of all mutants with ADP-r were followed by the acquisition of 1 H, [START_REF] Hellström | Partially uncleaved alphavirus replicase forms spherule structures in the presence and absence of RNA template[END_REF] N HSQC spectra (Figure 3a). The comparison of the number of shifted backbone amide signals qualitatively showed that all mutants bind ADP-r more weakly than WT (with the triple mutant losing completely the ability to bind ADP-r), strongly supporting the implication of the three mutated residues in the interaction.

Then, the ADP-r binding properties of three MAYV MD variants were quantified using ITC (Figure 3b). For MAYV MD WT the results confirm a 1:1 stoichiometry. The obtained dissociation constant Kd is 10.9 ± 0.7 µΜ, while the NMR titration data show a slow exchange process between free and ADP-r-bound MD. Thus, the slow exchange on the NMR time scale coincides with a relatively strong binding event (Figure 2c). The binding is enthalpically favored (∆Η = -3993 ± 68.55 cal/mol, ∆S = 9.09 cal/mol/K) indicating the formation of a complex through hydrogen bonds and hydrophobic interactions. The Kd for the N24A/F114A mutant is 52.1 ± 3.8 µΜ (∆Η = -2853 ±123.4 cal/mol, ∆S = 9.86 cal/mol/K).

In the case of the D10A mutant and ADP-r, the titration curve could not be fitted by a binding isotherm and no Kd value was obtained. This fact emphasizes the key role of the conserved residue Asp10 for ADP-r recognition and binding (Figure 3b).

Interaction of MAYV MD with RNA oligonucleotides

To gain insights into the RNA binding properties of MAYV MD, NMR titration studies and CSP analyses (Figure 4a, b, c, d ande) were performed to probe the interaction of 15 N-labeled MAYV MD with the following RNA substrates: 1) an artificial poly-A tail mimicking natural mRNA (N9A15, a 24mer with the sequence 5'CUCCGCCCCA153') and 2) 12mers of adenine (A12), cytosine (C12), guanine (G12), and uracil (U12) to investigate the specificity of RNA base recognition. Each titration was conducted up to a 1:1 molar ratio of protein:RNA oligonucleotide.

The CSP analyses were carried out for four out of five tested RNA oligonucleotides. There were no clear results for the NMR titration with the G12 oligonucleotide due to the intensity loss of all HN resonances at the final step of the titration and precipitation of the sample.

However, it should be noted that the first affected signals, at a MAYV MD:G12 ratio of 1:0.25, correspond to amino acids located at the periphery of the ADP-r binding pocket (His25-Glu28, Lys40-Trp41, Glu43, Val52, Asn97, Ser110, Gly117, Arg138, Lys146, Gln157).

Denaturing PAGE analysis (data not shown) revealed that the precipitate consisted of both protein and RNA. Similar observations have been reported for SUD-M and SUD-MC domains of SARS-CoV which have also an MD fold [START_REF] Johnson | SARS coronavirus unique domain: three-domain molecular architecture in solution and RNA binding[END_REF] .

The NMR data showed that the addition of the four RNA oligonucleotides affects the same regions of the polypeptide. In particular the largest CSPs were observed for amino acids in or at the periphery of the positively charged ligand-binding cavity. It is worth mentioning that the presence of the RNA ligands affects residues in this binding cleft, but not in other surface areas of the MD molecule, suggesting specificity of the interaction which is most likely driven by electrostatic attraction between positively charged surface residues of MD and negative charges of the RNA phosphate groups (Figure 4a, b, c, d, andS4).

The interactions between MAYV MD and A12 and N9A15 influence more amino acids (14 residues with CSPs above threshold and 10 HNs broadened beyond detection for A12; 17 residues with CSPs above threshold and 10 HNs broadened for N9A15) compared to U12 ( 9CSPs above threshold and 2 HNs broadened) and C12 (10 CSPs above threshold and 8 HNs broadened. In agreement with literature data [START_REF] Malet | The crystal structures of Chikungunya and Venezuelan equine encephalitis virus nsP3 macro domains define a conserved adenosine binding pocket[END_REF] , this finding strongly suggests that MD exhibits higher affinity to the adenine oligonucleotides (N9A15 and A12) than to U12 and C12.

Note that the latter oligonucleotides also influence MAYV MD residues in the two loops that seem to surround the negatively charged phosphate groups in the case of ADP-r binding (β2-α1 and β5-α3) 25 27 28 50 56 57 .

The MD binding surface of A12 is mainly defined by connecting loops and parts of helices α1 and α4. Specifically, the interaction surface is formed by the β1-β2 loop (Lys7, Asp10, Ile11, Ala12, Thr13, Ile15), residue Ala23 in the β2-α1 loop, the amino acids Val37 and Trp41 in the α1-helix, Ala48 in the α1-β3 loop, the residues Leu109, Ser110, Ile111, Ala116, Lys118 in the β5-α3 linker and the amino acids just before and at the beginning of helix α4 (Lys146, Trp148, Glu149). The same areas are affected by N9A15.

A comparison of the interaction data of A12 with the crystal structure of CHIKV MD (PDB ID: 3GPQ) [START_REF] Malet | The crystal structures of Chikungunya and Venezuelan equine encephalitis virus nsP3 macro domains define a conserved adenosine binding pocket[END_REF] in complex with an A3 RNA oligonucleotide reveals that in both cases the same protein regions participate in the interaction (Figure 4e). The higher number of amino acids that are influenced in the MAYV MD is due to the bigger size of the RNA oligonucleotide used here.

Complementary to the NMR titrations electrophoretic mobility shift assays (EMSA) were conducted to verify the above-mentioned results. Three out of five single-strand RNA oligonucleotides, namely N9A15 (5'CUCCGCCCCA153'), A12 and G12 were tested. Since MAYV MD does not enter the polyacrylamide gel under the tested condition (theoretical pI: 9.1), the EMSA results were evaluated in terms of the complex detection. The experiments showed that there is a concentration-dependent complex formation of MAYV MD with N9A15 (Figure 4f). Moreover, binding was observed for A12 and G12 (Figure S6).

To characterize further the selectivity of MAYV MD for nucleotides we monitored titrations of the protein with GDP, CDP and UDP through 1 H, [START_REF] Hellström | Partially uncleaved alphavirus replicase forms spherule structures in the presence and absence of RNA template[END_REF] N-HSQC spectra (CSP plots in Figure S5). While the addition of CDP and UDP barely causes shifts to the 1 H, [START_REF] Hellström | Partially uncleaved alphavirus replicase forms spherule structures in the presence and absence of RNA template[END_REF] N resonances of MAYV MD suggesting that there is no interaction, GDP is apparently bound weakly. The interaction with GDP occurs in the fast exchange regime and there are only five residues with significant chemical shift changes: Ile11, Gly32, Leu109, Tyr142, and Lys146 (Figure S5b). Compared to ADP binding (Figure S5a) the CSPs are much smaller, no resonances are broadened beyond detection, there are no meaningful CSPs for residues Ala23-Asn24 and Val69-Asn72, and the key residue Asp10 is only minimally involved in the interaction.

According to the available X-ray models 25 50 , in the presence of ADP the carboxyl group of Asp10 of MAYV MD is most likely hydrogen-bonded to the NH2 group of adenine. However, if the orientation of the two nucleotides in the MAYV MD binding cleft were the same this carboxyl group would repel the carbonyl oxygen on the C6 atom of guanine, which may explain the weaker interaction of GDP observed here. On the other hand, Ile11 usually forms a hydrogen bond to the base moiety of ADP and is still found to interact with GDP.

The data presented here indicate that both nucleotides and RNA oligonucleotides interact with MAYV MD. Our NMR data clearly identify the MD surface region involved in these biomolecular interactions. Adenine-rich RNA and ADP-r make use of the same binding cleft to dock onto the protein, but in the case of RNA binding, the α4 helix is additionally involved in the interactions. The observed interactions of MAYV MD with C12 and U12 RNA oligonucleotides are driven putatively by electrostatic attraction. This hypothesis is based on the absence of interaction between CDP or UDP and MAYV MD (Figure S5c andd). The oligonucleotides C12 and U12 seem to bind to the positively charged cleft of the protein (Figure 4b andc) due to their increased negative charge compared to the diphosphate nucleotides.

Concluding Remarks

Overall, this work reports the first solution structure of MAYV MD and identifies the protein residues involved in selective binding of ADP-r and adenine-containing ssRNA oligonucleotides to MAYV MD. This selectivity characterizes a number of MDs, and data reported herein shows that MAYV MD shares structural and functional features with other Alphavirus MDs. The MAYV MD polypeptide folds into a well-ordered structure in solution, characterized by the β/α/β "sandwich" topology and peptide segments with functionally relevant intrinsic flexibility. These segments surround a well-defined binding pocket, which provides high selectivity and variable affinity towards adenosine derivatives and adeninerich RNA oligonucleotides, with a preference for ADP-r. The critical role [START_REF] Malet | The crystal structures of Chikungunya and Venezuelan equine encephalitis virus nsP3 macro domains define a conserved adenosine binding pocket[END_REF] of Asp10 for ligand binding is confirmed by a dramatic reduction of binding affinity upon its mutation to Ala. Moreover, this study demonstrates that non-conserved residues in loop regions, like β2-α1, α1-β3, and β5-α3 induce different local conformations and backbone dynamics in MAYV MD compared to other homologous MDs, thus customizing its ligand binding properties. This finding may represent a general principle in the structure-function relationship of viral MDs and precludes the prediction of ligand interactions and function on the sole basis of sequence homology 30 51 . Note that ADP-r binding has recently been correlated in different cell lines with immediate or slow viral replication [START_REF] Luscher | ADP-Ribosylation, a multifaceted posttranslational modification involved in the control of cell physiology in health and disease[END_REF] .

According to the literature, there are striking examples of highly homologous MDs with distinct function, such as those from (a) CHIKV and Semliki Forest Virus (SFV) and (b) group 1 and group 3 CoV. In the former case the two MDs have ~77% sequence identity, but SFV does not bind ADP-r [START_REF] Neuvonen | Differential activities of cellular and viral macro domain proteins in binding of ADP-ribose metabolites[END_REF] . In the latter case a group 3 CoV MD exhibits great structural similarity to a group 1 CoV MD, but does not bind ADP-r whereas the group 1 CoV MD binds ADP-r with a Kd in the μΜ-range 61 62 . A similar finding has also been reported for MDs of the Alphavirus family. Some of us have observed different binding properties for the two MDs from VEEV and CHIKV, when they attempted to form the a MD -ADP-r complex [START_REF] Malet | The crystal structures of Chikungunya and Venezuelan equine encephalitis virus nsP3 macro domains define a conserved adenosine binding pocket[END_REF] . VEEV MD which has 59% sequence identity with CHIKV MD, was observed that it resists in ADP-r binding due to an unfavorable reorientation of the N24 and F114 side chains, when VEEV crystals soaked with ADPr. This finding suggests that even minor conformational differences in the loops β2-α1 and β5-α3 (Figure 1c), can fine-tune the substrate binding capacity of closely related MDs. The differences in structure and mobility of the MAYV MD loops β2-α1, α1-β3, and β5-α3 with respect to CHIKV (Figure 1c, S1 and S2) provides new insights into the structural basis of the functional role of these MD polypeptide segments.

Additionally, the structure elucidation of the MAYV MD along with the NMR driven interaction studies between the MD and ssRNA defines the positively charged surface illustrated in Figure 1d, which represents another binding site, distinct from the ADP-r binding crevice, and its role in binding negatively charged biopolymer(s) (Figure 4). These results are in line with the different roles of MD during the early and late stages of viral replication, in which a multi-functional domain is likely required [START_REF] Abraham | ADP-ribosyl-binding and hydrolase activities of the alphavirus nsP3 macrodomain are critical for initiation of virus replication[END_REF] . The observations presented here enlighten the structural and physicochemical determinants of MD-ligand interactions and provide valuable functional insights into the MD of this neglected virus reinforcing the efforts to design novel structure-based inhibitors of this Alphavirus protein domain.

(TCI) (acquisition parameters in Table S1). All chemical shifts were referenced to internal 2,2-dimethyl-2-silapentane-5-sulfonic acid (DSS) [START_REF] Melekis | NMR study of non-structural proteins-part I: 1H, 13C, 15N backbone and side-chain resonance assignment of macro domain from Mayaro virus (MAYV)[END_REF] .

The NMR data were processed with the TOPSPIN 3.3 software and analyzed with XEASY [START_REF] Bartels | The program XEASY for computer-supported NMR spectral analysis of biological macromolecules[END_REF] and CARA [START_REF] Keller | The Computer Aided Resonance Assignment Tutorial[END_REF] . Distance constraints for structure determination were obtained from [START_REF] Hellström | Partially uncleaved alphavirus replicase forms spherule structures in the presence and absence of RNA template[END_REF] N-and [START_REF] Kääriäinen | Functions of alphavirus nonstructural proteins in RNA replication[END_REF] C-edited 3D NOESY-HSQC experiments and 1 H-1 H 2D NOESY spectra, all recorded on the 700 MHz NMR spectrometer. 84 torsion angles constraints used in the structural calculations were derived from TALOS+ [START_REF] Shen | TALOS+: a hybrid method for predicting protein backbone torsion angles from NMR chemical shifts[END_REF] . Structure calculations were then performed through iterative cycles of DYANA [START_REF] Güntert | Torsion angle dynamics for NMR structure calculation with the new program Dyana[END_REF] , followed by restrained energy minimization with AMBER [START_REF] Pearlman | AMBER 5.0[END_REF] applied to each member of the final family of 20 DYANA models and the average model.

The NMR solution structure is characterized by a target function tf=0.73±7.87*10 -2 Å 2 and RMSD values of 1.07±0.14 Å and 1.84±0.14 Å for backbone and heavy atoms, respectively (Table S2). The coordinates of the energy-refined models (21 models, including the energy minimized average model) along with the structural constraints were deposited in the RCSB Protein Data Bank (PDB ID: 5IQ5). A small number of backbone and side-chain dihedral angles exhibit non-favorable values are due to minor gaps in the assignment and missing experimental.

All figures showing structures were generated using PyMol 66 .

N-Relaxation studies

Relaxation experiments ( 15 N T1, T2 and { 1 H}-15 N NOE) were conducted using the aforementioned spectrometers at 298 K. The delays for the [START_REF] Hellström | Partially uncleaved alphavirus replicase forms spherule structures in the presence and absence of RNA template[END_REF] N T1 were 7, 18, 40, 85, 150, 230, 350, 500, 680, and 900 ms, while delays of 16.99, 33.98, 50.98, 67.97, 84.96, 101.95, 135.94, 186.91, and 220.9 ms were used in the [START_REF] Hellström | Partially uncleaved alphavirus replicase forms spherule structures in the presence and absence of RNA template[END_REF]

NMR binding studies of ADP-r, nucleotides and RNA

The interactions of MAYV MD with adenosine-based ligands, GDP, CDP, UDP and ss RNAs were evaluated by comparing the 2D 1 H, [START_REF] Hellström | Partially uncleaved alphavirus replicase forms spherule structures in the presence and absence of RNA template[END_REF] N-HSQC spectra in the presence and absence of each tested molecule. The titrations were incremental and a 1 H, [START_REF] Hellström | Partially uncleaved alphavirus replicase forms spherule structures in the presence and absence of RNA template[END_REF] N-HSQC spectrum was recorded after each addition. The titrations were stopped when no further spectral changes were observed. Combined amide chemical shift perturbations were calculated using the equation: Δδ = [(ΔδHN) 2 +(ΔδN/5) 2 ] 1/2 67 . In order to apply a consensus approach for the identification of shifts that are large enough to be considered as indicators of the binding site in more than10 interaction studies between MAYV MD and small molecules or RNA oligos, we always applied the widely used method to calculate a threshold value 68 69 70 71 . This threshold is based on the calculation of the average CSP value and the standard deviation σ and sets an unbiased criterion to detect the residues that are significantly involved in an interaction.

Isothermal Titration Calorimetry (ITC)

ITC measurements were performed at 293 K using a MicroCal ITC200 instrument (Malvern).

Experiments were conducted in a buffer containing 25 mM HEPES, pH 7.5, 150 mM NaCl.

The protein concentration in the cell was between 80 µM (wt protein) and 120 µM (mutant protein) whereas the ADP-r concentration in the syringe was 600 to 800 µM. Heats of dilution were measured by injecting the ligand into the protein solution. Titration curves were fitted with the MicroCal Origin software, assuming a one-site binding model, and enthalpy (ΔH), entropy changes (∆S), dissociation equilibrium constants (Kd), and stoichiometry were extracted.

Electrophoretic mobility shift assays (EMSA)

Purified MAYV MD was incubated with ssRNA oligomers in reaction buffer containing 20 mM Tris base, pH 7, 100 mM KCl, 3 mM MgCl2, 0.5 mΜ EDTA, 1 mM DTT in DEPC H2O. The protein-nucleic acid mixtures were incubated at 37°C for 15 min, followed by 30 min incubation at 25°C. Then 6% glycerol was added, and the mixtures remained on ice for at least 15 min. The reactions were analyzed on 6% non-denaturing PAGE in the presence of 3 mM Tris-borate buffer, pH 8, at 80 V for 1 hour at 4 °C.

Nucleic acid was detected by GelRed (Biotium) staining and visualized using a UV light source equipped with a digital camera. GelRed was washed out and protein was subsequently detected by Coomasie staining (Bio-Rad).

The concentration of each ssRNA transcript was 15 µM. MAYV MD was used in increasing concentrations of 30, 60, 120, 150, 180, 210, and 300 µM.

MAYV MD does not enter the polyacrylamide gel under the conditions used here due to its basic pI (calculated pI = 9.1). Therefore, it could not be detected by PAGE analysis. However, its complexes with RNA carry enough negative charge to enter the gel and to be observed as discrete bands.

Accession numbers

The backbone assignments of the Mayaro virus macro domain have been deposited to the BMRB (Biological Magnetic Resonance Data Bank) under accession numbers 30043 (apoprotein) & 27612 (ADP-ribose complex), while the 3D structural models have been deposited in the Protein Data Bank (ID 5IQ5).

47 48 .

 48 The MAYV MD, spanning residues 1335 to 1493 of P1234, is characterized by a sixstranded β-sheet (β1 4-6, β2 18-21, β3 56-60, β4 63-67, β5 103-107, β6 138-142) surrounded at both sides by four major α-helices (α1 34-43, α2 78-98, α3 121-134, α4 149-157) (Figure 1b). The secondary structure elements are arranged in a ββαββαβαβα topology with β3 being the only antiparallel β-strand. This type of folding has also been encountered in other MDs exhibiting the typical three-layered β/α/β sandwich topology. As expected, its structural homologs according to DALI database 49 are the MDs of Chikungunya virus (CHIKV -PDB ID: 3GPG, Z score: 27, RMSD: 1.3 Å), Sindbis virus (SINV-PDB ID: 4GUA, Z score: 27, RMSD: 1.4 Å), and Venezuelan Equine Encephalitis virus (VEEV -PDB

  and adopts different orientations in the four structures. Specifically, MAYV, CHIKV and SINV MDs have at position 31 Asp or Glu instead of Gly in the VEEV sequence, while position 35 is usually occupied by charged amino acids, such as Arg (MAYV and CHIKV) or Lys (SINV), in contrast to a Gly residue in the VEEV sequence. Additionally, MAYV MD has Val and Ala at positions 29 and 38, respectively, whereas these positions are occupied by Pro and Tyr in the other three MDs.

  inset), Ala23-Gly27, Val37, Gly70, Leu108, Leu109, Thr111, Ser115, and Gly117 are shifted and broadened, exhibiting intermediate-to-fast exchange, different from the slow exchange process that they undergo during interaction with ADP-r. The HN resonances of the residues Ile11, Gly32, Val33, Val37, Val69, Ser110, and Ala116 could not be identified at the final step of the titration due to broadening of the corresponding peaks, indicating chemical exchange processes on an intermediate NMR time scale (Figure2e).

  N T2 experiments. The model free approach in the TENSOR CEA/CNRS60 program was used in order to obtain the S 2 values. The S 2 values of MAYV MD are presented in Figure 1e, while comparison of the S 2 values of the MDs from CHIKV and MAYV is shown in Figure S2. A H/D exchange experiment was performed with MAYV MD by dissolving a lyophilized sample of fully protonated uniformly 15 N labeled protein in 100% D2O and recording a series of 1 H, 15 N-HSQC spectra at 298 K on the 600 MHz NMR spectrometer. Each HSQC spectrum lasted 14 min 40 s and the midpoint of acquisition of the first spectrum was 11 min after dissolution of the sample in D2O.

Figure 1 .

 1 Figure 1. a. Characteristic NOE restraints derived from a 15 N-edited 3D NOESY-HSQC experiment. b. Ribbon representation of the energy minimized average NMR model of MAYV macro domain. c. Backbone superposition of macro domains from MAYV (grey-PDB 5IQ5), CHIKV (orange-PDB 3GPG), SINV (blue-PDB 4GUA), and VEEV (beige-PDB 3GQE) showing also the degree of sequence conservation in four loops. d. Representation of the electrostatic surface potential of MAYV macro domain (red: -4, blue: +4; calculated with the PyMol APBS extension). e. S 2 order parameters of the backbone of MAYV MD. Residues with values close to one have a decreased mobility on the ps/ns time scale. f. Mapping of residues with slowly exchanging backbone amide protons in the core region of MAYV macro domain (dark red). The corresponding 1 H, 15 N HSQC spectra at time points 11 min and 1125 min of the H/D exchange experiment are shown below.

Figure 2 a.

 2 Figure 2 a. Superposition of the 1 H, 15 N HSQC spectra of MAYV MD in the absence (black) and presence of ADP-r at a 1:1 molar ratio (red). Arrows indicate the changes of the peak positions of characteristic amino acids. b. Sequence alignment of PDB deposited MD structures in complex with ADP-r; symbols below indicate the key residues that interact with the ADP-r molecule. c. e. f. g. h. Plots of combined amide chemical shift changes, Δδ, induced by binding of adenine-based molecules to MAYV MD versus amino acid sequence at the given molar ratios of protein:ligand. In each plot the dotted line indicates the applied threshold of chemical shift perturbation calculated as the average value plus one standard deviation. An asterisk indicates an unassigned residue of apo-MD, p indicates a proline residue and a dot indicates a residue whose HN resonance could not be detected at each titration step. The Δδ values of Asn24 and Val33 in panel c (binding of ADP-r) are 2.2 and 3.8, respectively. Inset at right of each CSP plot: behavior of the amide peak of Thr13 during the titration indicating different binding regimes for the tested ligands (color code for the molar ratio protein:ligand: 1:0 black, 1:0.25 blue, 1:0.5 orange, 1:0.75 purple, 1:1 red, 1:3 green). d. Surface representation of MAYV MD depicting in dark red the amino acids exhibiting CSP values above the threshold in the presence of ADP-r.

Figure 3 .

 3 Figure3. a. Superposition of the 1 H,[START_REF] Hellström | Partially uncleaved alphavirus replicase forms spherule structures in the presence and absence of RNA template[END_REF] N HSQC spectra of MAYV MD mutants (D10A, N24A/F114A, D10A/N24A/F114A) in the absence (black) and presence of ADP-r at a 1:3 molar ratio (cyan). b. ITC data demonstrating the binding of ADP-r to MAYV MD (WT and mutants D10A and N24A/F114A). The corresponding Kd values were 10.9 µM for WT and 52.1 μM for N24A/F114A. For the D10A mutant a Kd value could not be calculated.

Figure 4 .

 4 Figure 4. a. b. c. d. Plots of combined amide chemical shift changes, Δδ, induced by the binding of RNA oligos to MAYV MD versus amino acid sequence at a 1:1 molar ratio of protein:RNA oligo (pink-A12, light yellow-U12, dark cyan-C12, grey-N9A15). In each plot the dotted line indicates the applied threshold of chemical shift perturbation calculated as average value plus one standard deviation. The surface representations above each plot show the amino acids exhibiting CSP values above the threshold upon binding of the respective RNA oligo to MAYV MD (color code: pink-A12, light yellow-U12, dark cyan-C12, grey-N9A15) and those with broadened amide cross peaks (in dark grey). An asterisk indicates an unassigned residue of apo-MD, p indicates a proline residue, a dot indicates a residue whose HN resonance could not be detected at each titration step. e. Left: The amino acids of MAYVMD that are critical for binding of the A12 oligo were mapped on its solution structure using only the amino acids with CSPs above the threshold (shown as pink sticks) and those, whose signals broadened beyond detection during the titration (in grey). Right: Representation of the residues (green sticks) that bind the A3 molecule in the CHIKV MD-A3 complex (PDB 3GPQ). f. Gel shift assay (EMSA) probing the interaction of MAYV MD with the 24mer 5'CUCCGCCCCA153' (=N9A15). The protein concentrations used are indicated above the gel.

Acknowledgments

We thank Prof. C. Stathopoulos and Dr. M. Apostolidi (Dept of Medicine, U of Patras) for useful discussions and their support and help with EMSA, and Dr. B. Fakler (U of Freiburg) for continued support. EU FP7 REGPOT CT-2011-285950 -"SEE-DRUG" project is acknowledged for financial support (DB, GAS) as well as for the purchase of UPAT's 700 MHz NMR equipment. The research work was supported by the Hellenic Foundation for Research and Innovation (HFRI) and the General Secretariat for Research and Technology (GSRT), under the HFRI PhD Fellowship grant (GA. no. 2430) and European program H2020 under the EVAg Research Infrastructure (grant agreement No. 653316)

Methods

Expression and purification of MAYV MD and its variants

The construction of the recombinant vector containing the gene of MAYV MD, the purification procedure of the protein, as well as the preparation of the labeled samples are described elsewhere [START_REF] Melekis | NMR study of non-structural proteins-part I: 1H, 13C, 15N backbone and side-chain resonance assignment of macro domain from Mayaro virus (MAYV)[END_REF] . The TRVL 4675 strain used in this study has a Thr at position 5 of the polypeptide sequence (GeneBank MK070492.1).

In order to obtain high purity for the binding assays a size exclusion chromatography was added, using a Superose 12 10/300 GL column (GE Healthcare) on an Äkta Purifier system (GE Healthcare) with 10 mM HEPES, pH 7.0, 20 mM NaCl, as elution buffer.

The fractions containing MAYV MD were pooled and concentrated to a final volume of 500 μL. The solution was then supplemented with 50 μL of D2O and DTT was added to a final concentration of 1-2 mM. The protein concentration in the NMR samples that were used for the titrations varied from 0.1 mM to 0.6 mM.

The MAYV MD mutants (D10A, N24A/F114A; D10A/N24A/F114A) were obtained by PCR mutagenesis reaction, using the plasmid containing the MAYV MD native gene as template and primers carrying the appropriate mutated codons. The expression and purification of the mutants were done according to the protocol for the native form of the protein.

Structure Determination

All NMR experiments were recorded at 298 K on a Bruker Avance 600 MHz NMR spectrometer, equipped with a cryogenically cooled pulsed-field gradient triple-resonance probe (TXI), and on a Bruker Avance III High-Definition four-channel 700 MHz NMR spectrometer equipped with a cryogenically cooled 5 mm 1 H/ 13 C/ 15 N/D Z-gradient probe