{'fraction_non_alphanumeric': 0.03713042930088886, 'fraction_numerical': 0.029912374708441027, 'mean_word_length': 4.7227633477633475, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 1, 'https://': 1, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
RNA G-quadruplex forming regions from SARS-2, SARS-1 and MERS coronoviruses

COVID-19 (Corona Virus Disease 2019), SARS (Severe Acute Respiratory Syndrome) and MERS (Middle East Respiratory Syndrome) are infectious diseases each caused by coronavirus outbreaks. Small molecules and other therapeutics are rapidly being developed to treat these diseases, but the threat of new variants and outbreaks argue for the identification of additional viral targets. Here we identify regions in each of the three coronavirus genomes that are able to form G-quadruplex (G4) structures. G4s are structures formed by DNA or RNA with a core of two or more stacked planes of guanosine tetrads. In recent years, numerous DNA and RNA G4s have emerged as promising pharmacological targets for the treatment of cancer and viral infection. We use a combination of bioinformatics and biophysical approaches to identify conserved RNA G4 regions from the ORF1A and S sequences of SARS-CoV, SARS-CoV-2 and MERS-CoV. Although a general depletion of G4-forming regions is observed in coronaviridae, the preservation of these selected G4 sequences support a significance in viral replication. Targeting these RNA structures may represent a new antiviral strategy against these viruses distinct from current approaches that target viral proteins.


Introduction SARS, MERS and SARS-2 coronaviruses
the COrona VIrus Disease 2019  is caused by infection of the beta-coronavirus SARS-CoV-2.
Coronaviruses are spherical viruses of 0.1 μm diameter, which contain an encapsidated RNA genome of positive polarity ( Figures  1A,B). The viral genome (around 30 Kb) encodes three types of viral proteins: structural, non-structural and accessory. In the case of the structural proteins, the RNA is coated by the nucleocapsid protein (N). The membrane (M) protein is responsible for binding both to the nucleocapsid and envelope (E) viral proteins to induce the necessary curvature for viral assembly (Schoeman and Fielding, 2020). The S protein is responsible of the specific binding to the angiotensin-converting enzyme 2 (ACE2) cellular receptor (Letko et al., 2020), (Walls et al., 2020) abundantly present in lung epithelia (Ziegler et al., 2020). The non-structural proteins also called viral replication and transcription complex (RTC) are responsible of the replication of the viral RNA and are encoded by two open-reading frames (ORF1a and ORF1b). The accessory proteins (3, 4, 5, 6, 7, 8 and 9) are not all conserved in the three viruses and are not essential for the viral cycle but play a role in virus host interaction enabling virus adaptation to host (Liu et al., 2014).

Current pharmacological treatments against COVID-19
Symptoms of the COVID-19 virus vary from asymptomatic to very severe and serious symptoms including fever, cough, chest tightness, headaches and other critical symptoms that might lead to death from the disease (Cui et al., 2019) (Ganesh et al., 2021) (WHO reported that 14.9 million excess deaths are associated with the . New variants of this virus are emerging through the countries making this virus a difficult threat even after more than 2 years of the date of its emergence (Zhou and Wang, 2021). Although many vaccines have been developed in record times no pharmaceutical therapeutic medications have been proven effective to treat the COVID-19 virus (Zhu F. C et al., 2020).
Several strategies are currently used to treat COVID-19 and other coronavirus diseases. These treatments are mainly based on repurposed drugs that block the viral cycle at different stages. For instance, Imatinib (Sisk et al., 2018), an Abelson kinase inhibitor is used to inhibit viral endocytosis. The viral proteases papainlike protease (PL pro ) and chymotrypsin-like protease (3CL pro ) (Thiel et al., 2003), responsible for the proteolytic processing of the viral proteins, have been targeted using a combination of Lopinavir and Ritonavir initially used for treating HIV, but have shown poor results against COVID-19 (Dong et al., 2020) Molnupiravir (Kabinger et al., 2021) and Remdesivir (de Wit et al., 2020),  are synthetic nucleoside analogues that inhibit viral RNA polymerases are also used against COVID-19. Among the very recent treatments, Nirmatrelvir is used as a specific 3C-like protease inhibitor (Owen et al., 2021). Finally, several Anti-SARS-CoV-2 monoclonal antibodies that target the spike protein have been developed. They have shown potent effects against SARS-CoV-2 infections, but these effects are strongly dependent on the circulating variant (Taylor et al., 2021). Despite these new approaches, COVID-19 is still progressing. Finding new strategies against this virus is still high priority. In this context, targeting RNA G4 structures may become an interesting avenue.
Structure and biological functions of G-quadruplex nucleic acids G-quadruplex nucleic acid structures (or G4s) are noncanonical nucleic acid secondary structures that can be formed by guanine rich DNA or RNA sequences (Jana et al., 2021) in which guanines are organized in several guanines arrays ( Figure 1C). To be able to form a G-quadruplex, 2 or more square planar arrangements of 4 guanines (G-quartets) connected through a Hoogsteen-type hydrogen bond needs to be stacked on top of each other ( Figures 1D,E). The abundant presence of K + or Na + cations in the cellular environment stabilizes these structures due to specific electrostatic interactions with the oxygen of the carbonyl groups of the four guanines. Depending on the direction of each guanine column, G4s can adopt different conformations: The parallel conformation when all the G-columns point in the same direction, the anti-parallel one, when 2 adjacent G-columns point in opposite direction and the hybrid one when 3 G-columns point in the same direction and one in opposite direction ( Figure 1F).
G4s have been gaining more interest in recent years. Many studies proved their impact on various biological processes such as gene expression mechanisms (Lago et al., 2021). The presence of G4s in the human genome and the role they play in disease modulation have been inspected and reviewed significantly in recent years mainly focused on their biological impact in telomeres (Bryan, 2020) and oncogene promoters (Brázda et al., 2021). In addition, it has been shown in vitro as well as in cultured cells that G4s play an important role in translational regulation when present in RNA coding regions that stimulate ribosomal frame-shifts (Cammas and Millevoi, 2016).
G-quadruplexes have already been identified in other viral genomes  and many studies have emphasized their role as pharmacological targets (Abiri et al., 2021), (Métifiot et al., 2014). Over a thousand G4 ligands have been portrayed in the G4 ligand database (Li et al., 2013) where at least one was tested in clinical trials for anticancer strategies (Onel et al., 2014) These small ligands that are able to bind the G4s present in viral genomes can be a potential new class of therapeutic agent to combat viral infections. Several G4 ligand families such as bisquinoliniums, naphthalene diimides and acridiniums had an inhibiting effect over DNA viruses such as herpesvirus HSV (Callegaro et al., 2017), EBV (Murat et al., 2014) and Hepatitis B virus (Biswas et al., 2017) as well as the RNA viruses HIV-1 and Hepatitis C virus (Jaubert et al., 2018). In contrast to classical antiviral approaches aimed at targeting viral proteins, the advantage of this strategy is that the viral RNA, the actual source of the disease, is directly targeted.
In this study, we identified evolutionary conserved G4 forming sequences in the genomes and negative RNA strands of SARS-CoV-2, SARS-CoV and MERS-CoV by using a combination of bioinformatics, Circular Dichroism, UVmelting, and NMR. Targeting these viral RNA structures may represent a promising new antiviral strategy.

Results
In silico detection of putative G4 forming sequences One important feature of RNA viruses is that they evolve extremely rapidly with up to 30% divergence at the sequence level between closely related subspecies (Firth, 2014). Sequence conservation relates to important residues in the encoded proteins, but also to conserved RNA secondary structures (Kiening et al., 2019) that have key functions during the virus life cycle. In HIV-1, such RNA structures include the TAR, DIS or RRE structural elements. In this study, we set out to identify possible G4-containing structural elements in the three coronavirus genomes based on bioinformatic prediction, sequence conservation and biophysical characterization.
The first step was to conduct an in silico screening of coronavirus genomes to detect the presence of putative G4 forming sequences (putative G4FS). We used the G4-Hunter scoring algorithm to search for putative G4FS over alignments of 2500 SARS-CoV-2, 308 SARS-CoV and 459 MERS-CoV viral genomes retrieved from the NCBI database. To avoid missing any possible G4 forming sequences (G4FS), we used high sensitivity settings for the G4-hunter algorithm with a threshold score of 1 and a sliding window of 20 nucleotides. As a result, we detected 16 putative G4FS in SARS-CoV-2, 21 in SARS-CoV and 40 in MERS-CoV ( Figure 2 and Supplementary Tables S1-S3). All identified sequences are highly conserved in the alignments as indicated by the logo representation of each putative G4FS (Supplementary Figures S1-S3). Due to emphasizing high sensitivity there is a corresponding drop in accuracy of the detection with an increased rate of false positives. Therefore, we manually inspected the 77 putative G4FS from the three coronaviruses (Supplementary Tables S1-S3) and discarded the sequences that presented single runs of 4-5 guanines. Indeed, despite a G4H-score>1, these sequences are unlikely to form mono-molecular G4 structures but will tend to form tetra-molecular G4s. We finally selected 24 putative G4FS from the 3 viruses ( Figure 2; Table 1) to undergo in vitro biophysical validation. In order to evaluate the actual propensity of these sequences to form G4s, we Frontiers in Chemistry frontiersin.org 03 synthesized the putative G4FS and analysed them by using the complementary techniques of Circular Dichroism (CD), Thermal Differential Spectra (TDS), UV-melting and 1 H-1D NMR. The results are summarized in Table 1  Analysis of putative G4 forming sequences by circular dichroism spectroscopy Circular Dichroism (CD) spectroscopy is a suitable tool to study the folding of nucleic acid structures. In the case of DNA structures, this technique successfully assigns specific CD spectra signatures to a given DNA topology including B-DNA, G4s or I-motifs (Villar-Guerra et al., 2018). In the case of RNA structures, CD spectra unfortunately show less discrimination. For instance, a stem-loop (or hairpin) RNA structure will give a similar spectral signature as for a G4 with a minimum at 245 nm and a maximum at 264 nm. However, since the CD signal is very sensitive to base stacking, it will still give valuable information on whether the structure is folded or unfolded. The CD spectra are summarized in the first column in Figure   In order to visualize all the genes, their sizes are not scaled to their actual lengths. The positions of the putative G4FSs are indicated by red triangles for sequences located on the (+) RNA and blue triangles indicate those located on the (−) RNA strand. Green triangles indicate putative G4FS detected by other softwares (Table 1). Putative G4FSs are detected on sequence alignments (SARS-2 = 2500 sequences, SARS = 308 sequences and MERS = 459 sequences) using G4Hunter. G4Hunter is a scoring algorithm that assigns each nucleotide with a score over a chosen window (20 nucleotides in this study). A positive score is given for the Gs and a negative for the Cs, whereas the other bases have a score of 0. The G4Hunter score is then calculated for each 20 nucleotides window. The average of the scores for each window of the alignment is then assigned to the first nucleotide of the sequence. Regions with scores greater than or equal to 1 (in absolute value) can form a G4, whereas those with scores below 1 (in absolute value) are unlikely to form a G4.
Frontiers in Chemistry frontiersin.org 04 tested display similar CD spectra with a single positive peak at 264 nm that indicates significant base stacking and suggests the formation of RNA secondary structures. Some sequences (S2-2 and M-5) also display a weaker positive peak at 290 nm which is not typically observed in RNA structures. Thermal Differential Spectra (TDS) are obtained by a simple arithmetic subtraction of the absorbance spectrum of the sequence at a high temperature (unfolded) state from that of the low temperature (folded) state (Mergny, 2005). This technique enables the identification of various nucleic acid structures such as duplexes, triplexes, i-motifs or G4s. For example, the observation of a negative peak at 295 nm would  Figures S4-S27).
Frontiers in Chemistry frontiersin.org 06 reveal the presence of a G4 forming sequence, whereas this peak is absent for hairpin duplexes. The second column of Figure 3 shows examples of G4-forming sequences from the three viruses. These spectra present a negative minimum at 297 nm and two or three positive maxima between 246 and 272 nm. These features are consistent with the formation of G4 structures. Conversely, in the case of S2-6 (Supplementary Figure S9), S2-7 (Supplementary Figure S10) S1-1 (Supplementary Figure S11), S1-7 (Supplementary Figure S17) Analysis of putative G4 forming sequences by UV-melting UV-melting experiments record the cooperative transition between the folded and unfolded forms of a given oligonucleotide. At 295 nm, a UV-melting-curve of a G4 structure usually results in a hypo-chromic (decrease of absorbance) melting process presenting a cooperative sigmoidal transition (Mergny and Lacroix, 2009). For instance, S2-1, S1-2, S1-4, M-3 and M-5 oligonucleotides present a reversible transition because the heating and cooling processes are superimposed. This reversibility indicates the folding and unfolding kinetics of the structure are relatively fast which strongly suggest the formation of a mono-molecular (formed by one strand) G4 structure (Mergny, 2005). The UV-melting experiments also allows us to estimate the stability of the G4 structures by measuring the melting temperature (T m ): 40°C for S2-1, 41°C for S1-2, 42°C for S1-4, 38°C for M-3 and 63°C for M-5. However, in most cases, the melting process is not reversible and presents a hysteresis (the heating and cooling profiles have two different half transitions: T 1/2 ), indicating relatively slow folding and unfolding kinetics. Important hysteresis (for instance 4°C for S2-2 and 25°C for S1-3) results from slower folding kinetics due to the formation of multimeric G4s (formed by two or more strands). This is observed for S1-3 (Supplementary Figure S13), S1-5 and S1-6 (Supplementary Figures S15, S16) Figures  S5-S7). In the case of the non-G4-forming sequences S2-5, S2-6, S2-7 (Supplementary Figures S8-S10), S1-1 (Supplementary Figure S11), S1-7 (Supplementary Figure S17 Figures  S8-S11, S13, S16-S18, S21) with T m ranging from 25°C to 65°C. In particular, S1-3, S1-5 and M1 show distinct behavior. At 260 nm, the melting profiles are reversible, whereas at 295 nm, a significant hysteresis is present. We interpret this behavior as an equilibrium between a multi-molecular G4 and a hairpin-like structure.

Analysis of putative G4 forming sequences by 1 H-1D NMR
The final characterization to confirm the formation of G4 structures used 1 H-1D NMR experiments. The imino proton H1 from guanine ( Figure 4A), or H3 in uracil, can normally only be observed by NMR when they are involved in hydrogen-bonds or are otherwise protected from exchanging with water molecules. Furthermore, their chemical shift is indicative of the type of baseparing involved . For standard base pairs, a GC guanine H1 and a UA uracil H3 have peaks at [12-13] ppm and [13][14] ppm, respectively. In contrast, the four imino H1 protons from guanines involved in an RNA G-tetrad corresponds to four peaks in the region from [10-12] ppm. However, it should also be noted that the guanine H1 and uracil H3 protons in GU or GG base-pairs will also present two peaks in the same region as the G-tetrad.
We present the 1 H-1D NMR spectra for a selection of G4FS in Figure 4. In the case of M-3, the presence of a set of well resolved peaks in the [10-12] ppm region confirms the formation of G4s, with a limited conformational polymorphism that is suggested by a second set of minor peaks. The presence of 12 major peaks also suggests the formation of a three-tetrad G-quadruplex. M-5, S2-4, S2-2 and S2-3 (Figures 4B-D) also exhibit G4 signatures with quite well resolved peaks in the [10-12] ppm region but the signal broadening indicate a higher polymorphism. S2-1, M-2, M-7, M-8, S1-2, S1-4 and S1-5 present a very broad NMR signal in the [10-12] ppm suggesting the formation of highly polymorphic G4 structures. Finally, very weak peaks were detected in the case of M10, indicating an aggregation or that only a small fraction of the sequences forms a G4. Peaks of higher intensity only appeared when the spectrum was recorded at 4°C.

Monomeric G4s are mainly located in the ORF1A and S genes
The results of the biophysical analysis of the 24 putative G4FS found in the three coronaviruses are summarized in Table 1 and Figure 5. By combining the results from the separate techniques, we found that 16 sequences were most likely to form G4 structures ( Figure 5A), with 8 forming other secondary Frontiers in Chemistry frontiersin.org 07 structures like hairpins. The 16 G4s are all highly conserved as suggested by the logo representations of the sequence alignments ( Figure 5B). Among the G4s, 11 sequences form multi-molecular G4s and 5 form mono-molecular G4s. Interestingly, monomolecular G4s are mainly present in ORF1A of the three viruses and S genes of SARS-1 and MERS.

Discussion
In this study, we wanted to contribute to the worldwide research effort to fight coronaviruses and in particular SARS-CoV-2. Since RNA G4s are becoming a promising antiviral target (Abiri et al., 2021), (Ruggiero and Richter, 2018), we wanted to investigate if G4s are formed in SARS-CoV-2, SARS-CoV-1 and MERS-CoV. In previous investigations, this effort has mainly been carried out through bioinformatics analyses (Zhao et al., 2021) , (Belmonte-Reche et al., 2020), (Bartas et al., 2020b) that could lead to the detection of false positives. In our approach, we believed that it was necessary to confront the predictions made in silico with an in vitro validation by biophysical methods. Indeed, we found that among the 24 sequences that we predicted and tested, only 5 were actually forming stable mono-molecular G4s. The biophysical validation of the predictions is indeed important; however, it is also essential to use separate biophysical techniques to obtain reliable conclusions. In fact, Circular Dichroism alone cannot distinguish a G4 RNA structure from a hairpin RNA because they both present a maximum around 260 nm and a minimum around 240 nm. The circular dichroism analysis in combination with the TDS and UV-melting techniques allows to better discriminate the two structures. For instance, the sequences S2-6 (Bartas et al., 2020b) and S2-7 (Belmonte-Reche et al., 2020) (Zhao et al., 2021) have been presented in other investigations as forming G4s. In our analysis they do not form G4s but hairpin structures instead. Moreover, in our approach we have stressed the reversibility of the thermal transition obtained by UV-melting analysis. This allowed us to get an idea of the kinetics of structure formation and therefore to have an insight on its molecularity. Indeed, the presence of an hysteresis in the melting profile argue for the presence of a multi-molecular form. In most previous studies on SARS-CoV-2 G4s, only thermal denaturation ramps (heating) were performed/shown and not renaturations (cooling), thus missing the information on the molecularity. Here, we found that the majority of the G4 structures detected in SARS-CoV-2 (11/ 16) were not mono-molecular, but multi-molecular which is a crucial information in the case of G4 RNA structures.
Notably the RNA G4s we validated are mostly present in the ORF1A gene of the positive strand and in the negative RNA strand of the S gene. As a result, the three coronavirus genomes show a significantly reduced frequency of G4 motifs. This trend has also been observed in several bioinformatics studies. For instance, G4-hunter algorithm has previously been applied to detect putative G4FS in several Nidoviral genomes (Bartas et al., 2020b). The highest frequency was obtained in Arteriviridae and the lowest in Coronaviridae. SARS-CoV-2 genome showed the most pronounced scarcity of G4s with a putative G4FS density in the lower part of the Coronaviridae family (Belmonte-Reche et al., 2021). Conversely, SARS-CoV-2 genome is rich in inverted repeats that can adopt RNA hairpins or cruciform structures (Goswami et al., 2021). These structural motifs are important for proper genome organization and are often specifically recognized by regulatory proteins. Hence, G4s have been effectively eliminated during evolution in order to favor hairpins or cruciform structures. The existence of these very rare but conserved G4 structures, despite this negative constraint, highlights their importance.
Hence, targeting these G4s and stabilizing them with small molecules could alter viral fitness. Indeed, two different studies have recently demonstrated SARS-CoV-2 inhibition by two well-known G4-ligands: PDS and TMPyP4. PDS compound, a G4-specific stabilizer, has been shown to attenuate SARS-CoV-2 infectivity in both pseudovirus cell systems and mouse models . More recently, TMPyP4 G4-ligands, already known for its antiviral activities against HIV-1, also exhibited antiviral activity against SARS-CoV-2. Remarkably, in Syrian hamster and transgenic mouse models of SARS-CoV-2 infection, administration of TMPyP4 resulted in a significant reduction in viral loads and lung injury (Qin et al., 2022). In the same study, the activity of TMPyP4 proved to be significantly stronger than that of remdesivir. These findings highlight the high potential of an antiviral strategy based on targeting G4 RNAs present in the SARS-CoV-2 genome. Inhibition mechanism is still unknown but it may be mediated by various mechanisms, as shown in Figure 6. Targeting ORF1a RNA G4s during translation Besides defining the amino acid sequence of proteins, an mRNA can also regulate gene expression through the formation of RNA secondary structures (Ding et al., 2014). These structures can have a direct impact on protein translation. Several studies have already proposed novel mechanisms of G4-mediated gene regulation at the translational level. For instance, RNA G4s can act as obstacles to the progression of ribosomes on their mRNA templates (Endoh and Sugimoto, 2016). Furthermore, depending on its stability, an RNA G4 located in an open reading frame can also inhibit protein expression (Endoh et al., 2013) by promoting termination and/ or ribosome stalling. In the case of Epstein Barr virus, the stabilization of RNA G4s by G4 ligands (pyridostatin derivatives) slows down the translation of viral mRNA into maintenance proteins (Murat et al., 2014).
Ribosomal frame-shifting is a well-known re-encoding process that occurs during translation elongation and is stimulated by secondary mRNA structures such as pseudo-knot stem-loop or G4s (Yu et al., 2014), (Howard et al., 2004). For instance, Ribosomal -1-frameshift (-1-FS) is a translational mechanism whereby ribosomes are forced to move back one nucleotide, leading to the decoding of a second open reading frame (ORF). It has been shown that G4 RNAs are able to stimulate this displacement in vitro and in cultured cells (Yu et al., 2014). The efficiency of the displacement is directly correlated with the stability of the G4 structure. Indeed, the stabilization of the G4 by PhenDC3 or PDS (Zhao et al., 2021) G4 ligand, induced an increase of the displacement rates. It has also been proposed that ribosomal frame-shift is stimulated by RNA dimerisation through the formation of kissing complexes in SARS-CoV (Ishimaru et al., 2013). Viral RNA dimerization is also possible via the formation of dimeric G4 structures as has already been proposed in HIV-1 (Shen et al., 2011). Such dimeric G4 intermediates could also be formed in SARS-CoV-2. One of the most important steps of the coronaviruses replication cycle is the primary translation of ORF1AB ( Figure 6 step C). It occurs in the early steps of the infection and results in the production of the RTC: the viral replication and transcription complex. A properly produced RTC is essential for the correct functioning of all downstream stages of the virus cycle. In our opinion, altering ribosomal progression by G4 ligands or changing the ORF could lead to deleterious consequences for the virus.

Targeting RNA G4s during RNA transcription
G4s not only act as roadblocks for ribosomes, but also perturb the progression of RNA polymerases and reverse transcriptases on their RNA template. For instance, in the case of HSV-1, in vitro assays revealed the ability of G4 ligands to stabilise HSV-1 G4s and block polymerase progression (Callegaro et al., 2017). Likewise, the stalling of viral RNA-dependent RNA polymerase, induced by the G4 ligands TMPyP4 and a PDS derivative (PDP), resulted in a decrease of HCV core protein expression (Wang et al., 2016). In the case of coronaviruses, such stalling could occur in the step D, E or F of the viral cycle ( Figure 6).
Several RNA binding proteins like hnRNP A1 or nucleolin have been shown to interact with G4s and some of them are involved in the replication of other viruses such as HCV (Bian et al., 2019), (Lista et al., 2017). Among the 16 non-structural proteins (Nsp1-16) encoded by ORF1a and ORF1b, Nsp3 protein is involved in the assembly of the viral RTC (Lei et al., 2018). NSP3 contains a domain called SUD, that is conserved in SARS-1 and SARS-2 and is responsible of the high pathogenicity of these two viruses. Interestingly, the SUD domain is able to interact with DNA and RNA G4s . Therefore, binding of NSP3 to the RNA G4s present in the viral genome may occur during the early steps of the viral cycle. Disrupting this interaction using G4 ligands can affect the correct progression of the replication cycle.
We have identified 5 RNA G4s in SARS-1, SARS-2 and MERS that are mostly present in ORF1A gene. The presence of these G4 structures, despite an overall negative evolutionary constraint, highlights their importance. Targeting these RNA structures may represent a promising new antiviral avenue with limited risks of developing resistance as suggested by the high conservation of the viral G4s.

Materials and methods
Bioinformatics detection of putative G4 forming sequences Several programs were used to detect the presence of guanine rich sequences that are likely to form G-quadruplexes. In the first step, Seaview software (Galtier et al., 1996) was used for the sequence visualization and alignment. Alignments have been made on all SARS-CoV-2 (2500 sequences), SARS-CoV (308 sequences) and MERS-CoV complete genomes (459 sequences) available on the National Center for Biotechnology Information [https://www.ncbi.nlm.nih.gov/ ]. From the alignments, the G4hunter software (Bedrat et al., 2016) was used to scan the alignment and detect sequences that can form G4 structures. G4H is a scoring algorithm that assigns each nucleotide with a score over a chosen window (20 in this study). A positive score is given for the Guanines and a negative for the Cytosines, whereas the other bases have a score of 0. The G4H score is then calculated for each 20 nucleotides window. Then an average of all the G4H scores of the alignment is calculated. To select the putative G4 forming sequences we chose to use a threshold of 1 (in absolute value) which is a good compromise between sensitivity and accuracy. The positive scores will predict G4FS that are present in the (+) strand while the negative scores will predict G4FS in the (-) strand. From the identified sequences, the WebLogo software (http://weblogo.berkeley. edu) (Crooks et al., 2004), (Schneider and Stephens, 1990) was used to visualize the conservation. In addition, the possibility of RNA structures other than G4s were predicted by using RNAfold (Gruber et al., 2008)

Oligonucleotides Synthesis
Oligonucleotides Synthesis: Oligonucleotide were synthesized in the lab according to the standard beta-phosphoramidite methodology using an H8 automated synthesizer with. The DMT protection group was left on the 5' terminus. Poly-pak purification procedure was performed to de-protect the synthetic oligonucleotides using poly-pak 2 cartridges. Oligonucleotide strand concentrations were determined by measuring absorbance at 260 nm with the extinction coefficient determined by an algorithm based on the nearest neighbour method using the IDT oligo calc online software.

Circular dichroism spectroscopy
CD spectra were measured in 10 mm path-length quartz cells using a (Jasco J-1500 instrument). Each final spectrum was the average of two scans between 220 and 335 nm at a rate of 50 nm/ min, with 0.5 nm data pitch, 2 nm bandwidth, and 2 s response time. Spectra were measured in 10 mm path-length quartz cells Frontiers in Chemistry frontiersin.org 11 with a 4 µM strand concentration in a buffer composed of 10 mM lithium cacodylate (pH 7.0) and 120 mM KCl. Spectra were recorded at 20°C, 15°C or 4°C.
Ultraviolet melting measurements and thermal difference spectra UV melting assays were measured on a SAFAS UVmc2 double beam spectrophotometer equipped with a 10cell holder and Peltier controller. Profiles were measured in 10 mm path-length quartz cells at the wavelengths of 260, 273 and 295 nm, with a 4 µM strand concentration in a buffer composed of 10 mM lithium cacodylate (pH 7.0) and 120 mM KCl. The heating and cooling curves were recorded between 90°C and 10°C with a temperature gradient of 0.4°C/min. Thermal difference spectra TDS were calculated by subtracting the UV spectrum at 10°C from the UV spectrum at 90°C.

NMR experiments
A 700 MHz Bruker spectrometer was used to perform the NMR experiments. Sample strand concentrations was around 100 µM dissolved in 20 mM potassium phosphate pH 7 and 120 mM KCl in 90% H2O/10% in 5 mm tubes. The 1D-1H-NMR spectra were recorded using a double pulse field gradient perfect spin echo (zgesgppe) pulse sequence to suppress the water signal. Spectra were recorded with a spectral width of 19 ppm, an acquisition time of 1.2s and a relaxation delay of 2s. Experiments were performed at 15°C and 4°C.

Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.

Author contributions
SA designed the study. SA, AK, PB, and CM co-wrote the paper. AK performed the bioinformatics analysis and the biophysical analysis. JM performed the NMR analysis. BV synthesized the oligonucleotides.

Funding
This work was supported by the Agence Nationale de Recherche contre le SIDA (ANRS), the Centre National de la Recherche Scientifique (CNRS), the Institut National de la Santé et de la Recherche Médicale (INSERM) and the Université de Bordeaux. JM benefits from a post-doctora fellowship from the ANRS (ECTZ103899 and ECTZ190071).