# Constrained G4 structures unveil topology specificity of known and new G4 binding proteins

A. Pipier, A. Devaux, T. Lavergne, A. Adrait, Y. Couté, S. Britton, P. Calsou,

J. Riou, E. Defrancq, D. Gomez

► To cite this version:

A. Pipier, A. Devaux, T. Lavergne, A. Adrait, Y. Couté, et al. Constrained G4 structures unveil topology specificity of known and new G4 binding proteins. Scientific Reports, 2021, 11 (1), 10.1038/s41598-021-92806-8. [https://doi.org/10.1038/s41598-021-92806-8](https://doi.org/10.1038/s41598-021-92806-8)

HAL Id: hal-03369749

https://hal.science/hal-03369749

Submitted on 11 Oct 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à 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.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

# Constrained G4 structures unveil topology specificity of known and new G4 binding proteins

A. Pipier<fnref n="1" />,², A. Devaux³, T. Lavergne D³, A. Adrait D⁴, Y. Couté D⁴, S. Britton D<fnref n="1" />,², P. Calsou D<fnref n="1" />,², J. F. Riou D⁵, E. Defrancq D³ & D. Gomez D<fnref n="1" />,²✉

G-quadruplexes (G4) are non-canonical secondary structures consisting in stacked tetrads of hydrogen-bonded guanines bases. An essential feature of G4 is their intrinsic polymorphic nature, which is characterized by the equilibrium between several conformations (also called topologies) and the presence of different types of loops with variable lengths. In cells, G4 functions rely on protein or enzymatic factors that recognize and promote or resolve these structures. In order to characterize new G4-dependent mechanisms, extensive researches aimed at identifying new G4 binding proteins. Using G-rich single-stranded oligonucleotides that adopt non-controlled G4 conformations, a large number of G4-binding proteins have been identified in vitro, but their specificity towards G4 topology remained unknown. Constrained G4 structures are biomolecular objects based on the use of a rigid cyclic peptide scaffold as a template for directing the intramolecular assembly of the anchored oligonucleotides into a single and stabilized G4 topology. Here, using various constrained RNA or DNA G4 as baits in human cell extracts, we establish the topology preference of several well-known G4-interacting factors. Moreover, we identify new G4-interacting proteins such as the NELF complex involved in the RNA-Pol II pausing mechanism, and we show that it impacts the clastogenic effect of the G4-ligand pyridostatin.

In the last twenty years, G-quadruplex structures (G4) emerged as cis-acting factors impacting almost all DNA and RNA transactions. G4 in telomeric sequences were first shown to play essential roles in telomeres capping and telomeres replication by telomerase<fnref n="1" />. Now, G4 are associated with the firing of DNA replication origins²,³, transcription initiation and termination, mRNA processing, mRNA transport⁴⁻⁶, translation⁷ and mitochondrial maintenance⁸.

G4 are noncanonical secondary structures formed by stacked tetrads of Hoogsteen hydrogen-bonded guanines bases, which are stabilized through the coordination of physiologically relevant cations (Na⁺, K⁺). G4 can result from the intramolecular folding of a unique G-rich sequence or from the intermolecular assembly of different G-rich containing strands⁹. An essential feature of G4 is their intrinsic polymorphic nature: numerous in vitro studies have revealed their ability to adopt different conformations, also called topologies¹⁰. Indeed, depending on the length and the composition of the sequence, as well as the environmental conditions (including the nature and concentration of metal cations, and local molecular crowding), a G4-forming sequence can adopt different topologies, in which the strands are in parallel, antiparallel or hybrid conformations, with the co-existence of different types of loops (lateral, diagonal or propeller) of variable lengths⁹⁻¹¹. In particular, this polymorphism is exacerbated for the human telomeric sequence and leads to intricate structural mixtures¹².

In cells, the impact of G4 on cellular metabolism is mainly associated with protein or enzymatic factors that bind, stabilize or resolve these structures. The folding of G-rich sequences into a G4, on DNA and RNA molecules, is associated with the formation of DSBs, transcription and translation repression and the alteration of the RNA processing¹³⁻¹⁵. To handle these major threats, cells use a battery of DNA and RNA helicases to control G4 formation¹⁶,¹⁷. Notably, most of helicases resolving DNA G4 are associated, when mutated, with genetic disorders, progeria and cancer progression (WRN, BLM, FANCI, RTEL), underlying the major impact of G4 structures on cell fitness¹⁸,¹⁹. In addition to helicases, the formation of G4 structures in cells is counteracted by proteins that bind single-stranded nucleic acids¹⁹,²⁰ through their OB-fold, RRM or RGG interacting

<footnote n="1">¹Institut de Pharmacologie et Biologie Structurale, IPBS, Université de Toulouse, CNRS, UPS, Toulouse, France. ²Equipe Labellisée Ligue Contre Le Cancer 2018, Toulouse, France. ³Département de Chimie Moléculaire, UMR CNRS 5250, Université Grenoble Alpes, 38058 Grenoble, France. ⁴CEA, INSERM, IRIG, BGE, Université Grenoble Alpes, 38000 Grenoble, France. ⁵Structure et Instabilité des Génomes, Muséum National d'Histoire Naturelle, CNRS, INSERM, CP 26, 75005 Paris, France. ✉email: dennis.gomez@ipbs.fr</footnote>
motifs<sup>21-23</sup>. Interestingly, RGG motif containing proteins also promote G4 stabilization<sup>24</sup> and control mRNAs localization through interaction with G4<sup>25</sup>. A major impact of G4 structures in cells is related to transcription<sup>4</sup>. Found enriched on promoters and transcriptional start sites (TSS)<sup>26,27</sup>, G4 structures have been shown to act predominantly as transcriptional repressors<sup>4,13,15</sup>, although some G4 have been also described as involved in transcription activation<sup>28,29</sup>. Furthermore, the presence of G4 motifs in the TSS proximal regions is associated with RNA-Pol II pausing sites and R-loops formation, two different factors promoting RNA-Pol II arrests and transcription-dependent DNA breaks<sup>30-34</sup>.

Given the increasing roles of G4 structures in cellular metabolism, extensive researches have been conducted in the last years in order to identify new G4-dependent mechanisms. Notably classical pull-down approaches identified hundreds of proteins associated to G-rich oligonucleotides forming G4 structures<sup>35-40</sup>. In solution, G-rich single-stranded molecules are in equilibrium between unfolded and folded states, and thus numerous identified G4 binding proteins are also able to recognize unfolded G-rich sequences<sup>20</sup>. In addition, G4 derived from single-stranded oligonucleotides can adopt different topologies<sup>9,11,41</sup>, that precludes to establish the specific contribution of each G4 topology to protein binding.

In this context, we have developed an approach to constrain the G4 into a single well-defined topology. The strategy is based on the use of a rigid cyclic peptide scaffold as a template for directing the intramolecular assembly of the anchored RNA or DNA oligonucleotides<sup>42-45</sup>. Moreover, such locked G4 display a thermal stability significantly higher than unconstrained G4 that strongly reduces the possibility to form unfolded single-stranded sequences. These constrained systems represent original tools, that we have used here for the identification and characterization of proteins interacting with a well-defined RNA or DNA G4 topology. In this study we identified through affinity purifications coupled to mass spectrometry (MS)-based quantitative proteomics a set of human proteins associated to locked G4 structures. Notably, this approach allowed us to identify NELF proteins as a new G4-interacting complex, leading us to investigate the impact of RNA-Pol II pausing mechanism into the response to G4 stabilization by G4 ligands.

## Results

**Identification and characterization of G4 associated proteins.** To identify human proteins interacting with G4 structures we performed classical pull-down assays followed by MS-based quantitative proteomic analysis. Various constrained G4 topologies based on the telomeric sequence (excepted for 5') were used (Supplementary Figure S1): systems **1, 6** and **7** depict a parallel topology and systems **2–5** have an antiparallel topology (Supplementary Figure S2)<sup>42–45</sup>. In our approach, biotin-functionalized G4-constrained molecules **1–7** and the biotin-functionalized duplex-DNA control **8**, were incubated individually with a semi-total human protein extract prepared from HeLa cells<sup>46</sup>, before being trapped using streptavidin magnetic beads to isolate interacting proteins (Fig. 1B). In a first-round assay, and in order to validate our approach, western-blotting analyses were performed to test the interaction and the binding specificity of some well-established and depicted G4-binding proteins to constrained G4 constructions. From these analyses we observed that eIF4G, WRN, Nucleolin, Mre11, DHX36, hnRNP A1 and CNBP, all well-known G4-interacting proteins<sup>20</sup>, were enriched using constrained G4 structures compared to the duplex control **8**. Conversely, the KU heterodimer, one of the most abundant human duplex-DNA-binding proteins<sup>47</sup>, was found enriched using duplex control **8** but was barely detectable on constrained G4 structures (Supplementary Figure S3A-B). Comparative analysis of the G4-interacting proteins enrichment on the six different constrained G4 structures **1a-3** and **5-7** shows a differential binding for these human proteins. Indeed, while eIF4G and Mre11 proteins are particularly enriched on constrained G4 **2** and **3** (i.e. antiparallel topology with two lateral loops), DHX36 and hnRNP A1 proteins are abundant on constrained G4 **1a, 6** and **7** (i.e. parallel topology without loops) (Supplementary Fig. 3A-B). Altogether, these first assays confirm that our pull-down strategy using constrained G4 structures allowed to both identify G4 binding proteins and to discriminate important aspects of their structural interaction with G4. Thus, our approach represents a powerful tool to find new proteins recognizing particular G4 conformations and prompted us to proceed to an extensive MS-based quantitative proteomic analysis of human proteins interacting with two particular constrained G4 structures, systems **1a** (i.e. parallel without loops) and **2** (i.e. antiparallel with two lateral loops), compared to duplex control **8** (Fig. 1A-B).

Based on three independent experiments, MS-based proteomic analyses identified a total of 1021 proteins interacting with constrained structures (see Materials and Methods for identification conditions) (Supplementary Table 1). This total corresponds to the sum of proteins interacting with **1a-2** and control **8** (duplex) (Fig. 1C and Supplementary Table 1). Using a label-free quantification and statistical filtering to compare the abundances of the proteins eluted from different constructions (see Materials and Methods for filter conditions), we identified 425 proteins enriched on constrained G4 structures **1a-2** compared to duplex control **8** (fold change ≥2 and *p*-value < 0.05, allowing to reach a false discovery rate (FDR) inferior to 5%) (Fig. 1C and Supplementary Table 2). These proteins belong to six significant enriched KEGG pathways cluster (*p* < 0.05): Spliceosome, RNA transport, RNA degradation, mRNA surveillance, DNA replication and Homologous Recombination (Fig. 2A). To go further, enriched GO Biological Processes and Molecular Functions terms were also determined. This analysis revealed that proteins enriched on constrained-G4 structures are mainly associated with DNA and RNA transactions (Fig. 2B-C), in agreement with current knowledge on genomic localization and biological function of G4. In line with these data, 214 out of 425 proteins enriched on G4 structures have been described as nucleic acid interacting factors, as indicated by terms from GO Molecular functions data analyses (Fig. 1C and Supplementary Table 3). Furthermore, the KEGG pathways clusters from these 214 nucleic acid binding proteins correspond to almost the same biological processes defined by the complete set of G4-interacting proteins (Supplementary Figure S4). An additional analysis of nucleic acid binding proteins enriched on constrained G4 structures defines five functional groups covering (i) ATP dependent DNA/RNA helicases activities, (ii) hnRNP
<figure><img src="image_3.png" /><figcaption>B</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>C</figcaption></figure>

<figure><img src="image_5.png" /></figure>

**Figure 1.** (A) Schematic representation of constrained DNA structures used in the pull-down assay. (B) Global strategy to identify constrained G4 interacting proteins from human cells. Biotin-functionalized G4-constrained molecules (1a and 2) and the biotin-functionalized duplex-DNA control 8 were individually mixed with a semi-total human protein extract from HeLa cells, then trapped by streptavidin magnetic beads to isolate interacting proteins. Protein identification was obtained from MS-based quantitative proteomic analysis and further characterized by western blotting (arrow), or directly by western blotting (dashed arrow). (C) Diagram showing the differential enrichment of human proteins on constrained G4 structures relative to control duplex DNA. G4 enriched proteins refer to proteins found enriched on 1a and/or 2 G4 constructions relative to the duplex control 8. 214 out of 425 proteins found enriched on constrained G4 have been shown to interact with nucleic acids. Differentially interacting proteins were sorted out using a fold change ≥2 and *p*-value <0.05, allowing to reach a false discovery rate (FDR) inferior to 5% according to the Benjamini–Hochberg procedure.

proteins, (iii) proteins involved in the polyadenylation process, (iv) a large group of proteins involved in splicing and (v) proteins related to small nuclear ribonucleoprotein complexes (Table 1).

Next, using UniprotKB, gene Ontology (GO), G4IPDB⁴⁸ databases and G4 search terms in PubMed, we determined that eighteen proteins found enriched on constrained G4 structures **1a-2** have been already implicated in G4 biology (Table 2).

We compared our results with two recent studies that identified proteins associated with RNA G4 structures³⁵,³⁶ and found that 98 out of the 425 proteins identified in our study have been already shown to interact with G4 structures (Fig. 3A). Fourteen factors are in common in the three studies, including DHX36 and DDX3 proteins, two major G4 resolvases (Fig. 3A). Finally, in order to further investigate the association of constrained G4 interacting proteins identified here with G4 biological functions, we compared our list of 425 proteins with the list of 758 G4 sensitizers genes, the deficiency of which leads to increased sensitivity to G4 ligands, established by Zyner et al.⁴⁹. From this analysis, we determined that at least 62 out of the 425 proteins enriched on constrained G4 structures **1a-2** were reported as G4 sensitiser proteins (Fig. 3B).

**Different sets of proteins are enriched on particular G4 conformations.** Statistical analysis of the abundances of proteins enriched on G4 structures **1a-2** shows differential interactions for some of them with the two constrained G4 conformations. Indeed, among the proteins with log2 (fold change **1a/2**) ≥ 1 and *p*-values < 0.05, we established three groups of constrained G4 interacting proteins (Fig. 4A). The first and largest group is composed of 204 proteins found significantly enriched on construction **2** (i.e. antiparallel topology with two lateral loops) compared to construction **1a** (i.e. parallel topology without loop). The second group comprises 190 proteins without significant enrichment on one particular conformation (−2 < fold change **1a/2** < 2). Finally, we found a third group with only 31 proteins significantly enriched on construction **1a** relatively to structure **2** (Fig. 4A and Supplementary Tables 4.1, 4.2 and 4.3). KEGG analysis of the first group and Common proteins showed that both groups define similar significant enriched pathways than those covered by total constrained G4 interacting proteins (Spliceosome, mRNA surveillance, RNA degradation, DNA replication and repair) (Supplementary Figure S5). Surprisingly, the third group of proteins enriched on constrained G4 structure **1a** defines
<figure><img src="image_6.png" /><figcaption>B</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>C</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 2.** Most significant pathways and processes covered by constrained G4 interacting proteins. (A) Enriched KEGG pathways. Gene Ontology terms, (B) GO-Biological processes and (C) GO-Molecular Functions for the 425 proteins found enriched on constrained G4 structures. A right-sided (Enrichment) test based on the hyper-geometric distribution was performed on the corresponding Entrez gene IDs for each gene list and the Bonferroni adjustment (*p* < 0.05).</figcaption></figure>

a unique and highly enriched KEGG pathway cluster (*p* < 0.0005) consisting of factors involved in aminoacyl-tRNA biosynthesis⁵⁰. Further analysis based on the STRING protein-protein interactions database⁵¹ unveiled that all of the components of the multi-tRNA synthetase complex (MSC) are enriched on constrained G4 **1a** (Fig. 4B-C). Indeed, in addition to eight cytoplasmic aminoacyl tRNA synthetase enzymes (methionyl MARS, glutaminyl DARS, lysyl KARS, arginyl RARS, isoleucyl IARS, leucyl LARS, aspartyl DARS and glutamyl-, prolyl EPRS-tRNA synthetase) composing this complex, our pull down assay also isolated the three non-enzyme components (AIMP1, AIMP2 and AIMP3, also known as EEF1E, proteins) of the MSC complex (Fig. 4B-C). Furthermore, western blotting analysis showed that AIMP1 interacts with constrained and unconstrained G4 structures in vitro (Fig. 5).

**Impact of orientation and nucleotide composition of connecting loops on the differential association of proteins on G4 structures.** Since most of the proteins characterized in this study are associated with construction **2** (i.e. antiparallel with two lateral loops), we explored the impact of the terminal single-stranded extension, loop sequence and loops orientation on the binding of proteins to constrained G4 molecules. First, western blotting analysis of the binding of nine selected factors, found enriched on constrained G4 structures relatively to duplex control **8**, indicates that the relative orientation or the nucleotide loop composition has a not major impact on the binding of these proteins to constrained G4 structures (Fig. 5). Indeed, we found that protein signals obtained with systems **2** (i.e. with both loops oriented in the same 5'-3' sense), **3** (i.e. with loops oriented in the 5'-3' and 3'-5' senses, respectively), and **4** in which the ATT sequence of the external loops was replaced by a TCT, are not significantly different (Fig. 5). Next, in order to explore the impact of the extension length, we constructed a new system **1b** with a three-nucleotide extension consisting of the 5' TTA 3' sequence. Western-blotting analysis shows that the addition of a supplementary nucleotide on the terminal single-stranded extension of construction **1b** significantly improves the binding of some factors (WRN, DHX36, Mre11, NELF-E and AIMP1). However, intensity signals observed on constrained G4 structures with connecting loops (i.e. systems **2-4**) remain considerably stronger relatively to the signals obtained with both **1a** and **1b** constructions, except for AIMP1 (Fig. 5). These results indicate that longer single-stranded extensions improve
<table><caption>**Table 1.** Manually curated functional groups from nucleic acid binding proteins enriched on constrained G4 structures.</caption><thead><tr><th>Mapped IDs</th><th>Gene Name</th><th>Mapped IDs</th><th>Gene Name</th></tr></thead><tbody><tr><td>DDX1</td><td>ATP-dependent RNA helicase DDX1</td><td>SF1</td><td>Splicing factor 1</td></tr><tr><td>DDX20</td><td>Probable ATP-dependent RNA helicase DDX20</td><td>SF3A1</td><td>Splicing factor 3A subunit 1</td></tr><tr><td>DDX23</td><td>Probable ATP-dependent RNA helicase DDX23</td><td>SF3A2</td><td>Splicing factor 3A subunit 2</td></tr><tr><td>DDX3X</td><td>ATP-dependent RNA helicase DDX3X</td><td>SF3A3</td><td>Splicing factor 3A subunit 3</td></tr><tr><td>DDX41</td><td>Probable ATP-dependent RNA helicase DDX41</td><td>SF3B1</td><td>Splicing factor 3B subunit 1</td></tr><tr><td>DDX42</td><td>ATP-dependent RNA helicase DDX42</td><td>SF3B2</td><td>Splicing factor 3B subunit 2</td></tr><tr><td>DDX52</td><td>Probable ATP-dependent RNA helicase DDX52</td><td>SF3B3</td><td>Splicing factor 3B subunit 3</td></tr><tr><td>DDX6</td><td>Probable ATP-dependent RNA helicase DDX6</td><td>SF3B6</td><td>Splicing factor 3B subunit 6</td></tr><tr><td>DHX16</td><td>Pre-mRNA-splicing factor ATP-dependent RNA helicase DHX16</td><td>SRSF1</td><td>Serine/arginine-rich splicing factor 1</td></tr><tr><td>DHX29</td><td>ATP-dependent RNA helicase DHX29</td><td>SRSF10</td><td>Serine/arginine-rich splicing factor 10</td></tr><tr><td>DHX30</td><td>Putative ATP-dependent RNA helicase DHX30</td><td>SRSF11</td><td>Serine/arginine-rich splicing factor 11</td></tr><tr><td>DHX36</td><td>ATP-dependent RNA helicase DHX36</td><td>SRSF2</td><td>Serine/arginine-rich splicing factor 2</td></tr><tr><td>DHX38</td><td>Pre-mRNA-splicing factor ATP-dependent RNA helicase DHX38</td><td>SRSF6</td><td>Serine/arginine-rich splicing factor 6</td></tr><tr><td>DHX40</td><td>Probable ATP-dependent RNA helicase DHX40</td><td>SRSF7</td><td>Serine/arginine-rich splicing factor 7</td></tr><tr><td>hnRNP A1</td><td>Heterogeneous nuclear ribonucleoprotein A1</td><td>SRSF9</td><td>Serine/arginine-rich splicing factor 9</td></tr><tr><td>hnRNP A2B1</td><td>Heterogeneous nuclear ribonucleoproteins A2-B1</td><td>SNRNP200</td><td>U5 small nuclear ribonucleoprotein 200 kDa helicase</td></tr><tr><td>hnRNP A3</td><td>Heterogeneous nuclear ribonucleoprotein A3</td><td>SNRNP40</td><td>U5 small nuclear ribonucleoprotein 40 kDa protein</td></tr><tr><td>hnRNP F</td><td>Heterogeneous nuclear ribonucleoprotein F</td><td>SNRNP70</td><td>U1 small nuclear ribonucleoprotein 70 kDa</td></tr><tr><td>hnRNP H1</td><td>Heterogeneous nuclear ribonucleoprotein H1</td><td>SNRPA</td><td>U1 small nuclear ribonucleoprotein A</td></tr><tr><td>hnRNP H3</td><td>Heterogeneous nuclear ribonucleoprotein H3</td><td>SNRPA1</td><td>U2 small nuclear ribonucleoprotein A</td></tr><tr><td>hnRNP L</td><td>Heterogeneous nuclear ribonucleoprotein L</td><td>SNRPC</td><td>U1 small nuclear ribonucleoprotein C</td></tr><tr><td>hnRNP R</td><td>Heterogeneous nuclear ribonucleoprotein R</td><td>SNRPD1</td><td>Small nuclear ribonucleoprotein Sm D1</td></tr><tr><td></td><td></td><td>SNRPD2,SNRPD1</td><td>Small nuclear ribonucleoprotein Sm D2</td></tr><tr><td>CPSF1</td><td>Cleavage and polyadenylation specificity factor subunit 1</td><td>SNRPD3</td><td>Small nuclear ribonucleoprotein Sm D3</td></tr><tr><td>CPSF2</td><td>Cleavage and polyadenylation specificity factor subunit 2</td><td>SNRPE</td><td>Small nuclear ribonucleoprotein E</td></tr><tr><td>CPSF3</td><td>Cleavage and polyadenylation specificity factor subunit 3</td><td>SNRPN</td><td>Small nuclear ribonucleoprotein-associated protein N;SNRPN</td></tr><tr><td>CPSF4</td><td>Cleavage and polyadenylation specificity factor subunit 4</td><td></td><td></td></tr><tr><td>CRNKL1</td><td>Crooked neck-like protein 1</td><td></td><td></td></tr><tr><td>CSTF1</td><td>Cleavage stimulation factor subunit 1</td><td></td><td></td></tr><tr><td>CSTF3</td><td>Cleavage stimulation factor subunit 3</td><td></td><td></td></tr></tbody></table>

the binding of proteins to type 1 constructions. Unconstrained G4 structures (free) formed by the telomeric and the c-myc promoter sequences were also used to further characterize their interaction (Fig. 5, right part). Most of the factors enriched on constrained G4 structures also present a significant interaction with non-constrained G4-forming sequences, relative to the scramble sequence, indicating a selective binding of these proteins to G4 structures. However, proteins showing a significant interaction with 1a and 1b systems (hnRNP A1, AIMP1, WRN and at lesser extent DHX36) also display a strong interaction with the scramble non-G4 forming oligonucleotide, confirming their ability to bind both G4 structures and single-stranded unfolded G-rich DNA.

**Constrained G4 structures identify the NELF complex as a new G4 interacting factor.** MS-based quantitative proteomic analyses identified NELF-A, NELF-B, NELF-E and NELF-C/D, all members of the NELF complex³², as enriched on constrained G4 structures relatively to the duplex control construction 8 (Fig. 6A, Supplementary Table 2). Selective interaction of NELF complex proteins with G4 structures was further confirmed through western blotting analysis performed using both constrained and unconstrained G4 structures (Fig. 5). In order to further characterize the interaction of the NELF complex with G4 structures, we performed the reverse experiment in which an immunoprecipitated NELF complex was used to investigate its interaction with constrained G4 structures. For that, the NELF complex was immunoprecipitated from a HeLa cell line overexpressing an ectopic Flag-tagged form of the NELFE protein. After elution by competition with a Flag peptide (Supplementary Figure S6), the NELF complex was used in pull-down experiments using
<table><caption>**Table 2.** Constrained G4 interacting factors found related to G4 on UniprotKB, gene Ontology (GO), PubMed abstract and G4IPDB⁴⁸ data bases.</caption><thead><tr><td>Mapped IDs</td><td>Gene Name</td><td>PUBMED ID</td></tr></thead><tbody><tr><td>ADAR</td><td>Double-stranded RNA-specific adenosine deaminase</td><td>24813121 23381195</td></tr><tr><td>CNBP</td><td>Cellular nucleic acid-binding protein</td><td>23774591 28329689 24594223 26332732 31219592</td></tr><tr><td>DDX1</td><td>ATP-dependent RNA helicase</td><td>29731414</td></tr><tr><td>DDX42</td><td>ATP-dependent RNA helicase</td><td>31287417</td></tr><tr><td>DHX36</td><td>ATP-dependent RNA helicase</td><td>29269411 28069994 25653156 25611385 24151078 22238380 21149580 18842585 16150737</td></tr><tr><td>DNMT1</td><td>DNA (cytosine-5)-methyltransferase 1</td><td>30275516</td></tr><tr><td>EWSR1</td><td>RNA-binding protein</td><td>21244633 21561087 22214309</td></tr><tr><td>FUS</td><td>RNA-binding protein</td><td>18776329 19749353 23521792 24251952 28575444 29434328 29800261</td></tr><tr><td>hnRNP A1</td><td>Heterogeneous nuclear ribonucleoprotein A1</td><td>9188487 19282454 20213319 24371143 24831962 26930004 28510424 29361764 30247678 31311954</td></tr><tr><td>hnRNP A2B1</td><td>Heterogeneous nuclear ribonucleoproteins A2/B1</td><td>15302914 17716999</td></tr><tr><td>hnRNP A3</td><td>Heterogeneous nuclear ribonucleoprotein A3</td><td>27623008 23381195</td></tr><tr><td>hnRNP F</td><td>Heterogeneous nuclear ribonucleoprotein F</td><td>29269483</td></tr><tr><td>hnRNP H1</td><td>Heterogeneous nuclear ribonucleoprotein H</td><td>26930004 27623008</td></tr><tr><td>MID1</td><td>E3 ubiquitin-protein ligase Midline-1</td><td>21930711</td></tr><tr><td>Mre11A (*)</td><td>Double-strand break repair protein MRE11</td><td>16116037</td></tr><tr><td>RIF1</td><td>Telomere-associated protein</td><td>26436827 29348174 29357064 30510058 31197198</td></tr><tr><td>SF3B3</td><td>Splicing factor 3B subunit 3</td><td>23381195</td></tr><tr><td>SRSF1</td><td>Serine/arginine-rich splicing factor 1</td><td>24771345</td></tr></tbody></table>

<figure><img src="image_9.png" /><figcaption>**Figure 3.** Constrained-G4 interacting factors are associated with RNA-G4 binding activities and with the sensitisation to small molecules that stabilize G4 structures. (A) Venn diagram showing the overlap of our study (orange) with the RNA-G4 interacting proteins identified in Herdy³⁵ (green) and Herviou³⁶ (blue). 98 out of 425 proteins identified in our study were known to interact with RNA-G4 structures, with 14 indicated factors common to three studies. (B) Schematic representation of constrained G4 interacting proteins identified in our study that are associated to an increased sensitivity to G4 ligands, established by Zyner et al.⁴⁹. From this analysis we determined that 62 out of 425 proteins were reported as G4 ligands sensitisers.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>B</figcaption></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /><figcaption>C</figcaption></figure>

<table><thead><tr><th rowspan="2">Protein names</th><th rowspan="2">Gene names</th><th rowspan="2">Razor + unique peptides</th><th colspan="2">2 vs 8</th><th colspan="2">1a vs 8</th><th colspan="2">2 vs 1a</th></tr><tr><th>Significant</th><th>p-value</th><th>log2(Fold Change)</th><th>Significant</th><th>p-value</th><th>log2(Fold Change)</th><th>Significant</th><th>p-value</th><th>log2(Fold Change)</th></tr></thead><tbody><tr><td>Aminoacyl tRNA synthase complex-interacting multifunctional protein 1; Endothelial monocyte-activating polypeptide 2</td><td>AIMP1</td><td>4</td><td></td><td>5.99E-01</td><td>-0.26</td><td>+</td><td>1.48E-04</td><td>2.90</td><td>+</td><td>7.74E-05</td><td>-3.15</td></tr><tr><td>Aminoacyl tRNA synthase complex-interacting multifunctional protein 2</td><td>AIMP2</td><td>2</td><td>+</td><td>9.78E-03</td><td>1.44</td><td>+</td><td>4.34E-05</td><td>3.22</td><td>+</td><td>3.03E-03</td><td>-1.77</td></tr><tr><td>Isoleucine-tRNA ligase, cytoplasmic</td><td>IARS</td><td>11</td><td>+</td><td>9.30E-04</td><td>1.94</td><td>+</td><td>3.92E-06</td><td>3.92</td><td>+</td><td>7.94E-04</td><td>-1.98</td></tr><tr><td>Leucine-tRNA ligase, cytoplasmic</td><td>LARS</td><td>11</td><td>+</td><td>7.45E-03</td><td>1.97</td><td>+</td><td>8.74E-05</td><td>3.82</td><td>+</td><td>1.05E-02</td><td>-1.85</td></tr><tr><td>Methionine-tRNA ligase, cytoplasmic</td><td>MARS</td><td>5</td><td></td><td>1.00E+00</td><td>1.14</td><td>+</td><td>2.49E-06</td><td>4.22</td><td>+</td><td>3.36E-05</td><td>-3.09</td></tr><tr><td>Bifunctional glutamate/proline-tRNA ligase; Glutamate-tRNA ligase; Proline-tRNA ligase</td><td>EPRS</td><td>12</td><td>+</td><td>3.49E-02</td><td>1.23</td><td>+</td><td>5.31E-04</td><td>2.58</td><td>+</td><td>2.33E-02</td><td>-1.35</td></tr><tr><td>Arginine-tRNA ligase, cytoplasmic</td><td>RARS</td><td>17</td><td></td><td>1.48E-01</td><td>0.71</td><td>+</td><td>2.25E-05</td><td>3.54</td><td>+</td><td>1.33E-04</td><td>-2.82</td></tr><tr><td>Eukaryotic translation elongation factor 1 epsilon-1</td><td>EEF1E1; EEF1E1-BLOC1S5</td><td>2</td><td></td><td>1.55E-01</td><td>0.45</td><td>+</td><td>3.63E-05</td><td>2.14</td><td>+</td><td>2.24E-04</td><td>-1.69</td></tr><tr><td>Lysine-tRNA ligase</td><td>KARS</td><td>11</td><td></td><td>4.94E-02</td><td>1.08</td><td>+</td><td>2.18E-04</td><td>2.79</td><td>+</td><td>5.52E-03</td><td>-1.71</td></tr><tr><td>Glutamine-tRNA ligase</td><td>QARS</td><td>11</td><td></td><td>2.11E-01</td><td>0.51</td><td>+</td><td>2.33E-05</td><td>2.98</td><td>+</td><td>1.04E-04</td><td>-2.46</td></tr><tr><td>Aspartate-tRNA ligase, cytoplasmic</td><td>DARS</td><td>13</td><td>+</td><td>1.78E-02</td><td>1.14</td><td>+</td><td>2.86E-05</td><td>3.01</td><td>+</td><td>9.92E-04</td><td>-1.87</td></tr></tbody></table>

**Figure 4.** Differential interaction of human proteins with constrained G4 structures adopting different topologies. (A) Diagram showing the differential enrichment of human proteins on **1a** (green) and **2** (purple) constructions. Differential enrichment of proteins on structures **1a** or **2** was determined through statistical analysis using the fold change ≥ 2 and *p*-value < 0.05, allowing to reach a false discovery rate (FDR) inferior to 5%. (B) KEGG pathway covered by the 31 proteins found enriched on constrained G4 structure **1a** and functional interaction network analysis using STRING⁵¹ for the 31 proteins found enriched on construction **1a**. A right-sided (Enrichment) test based on the hyper-geometric distribution was performed on the corresponding Entrez gene IDs for each gene list and the Bonferroni adjustment (*p* < 0.05). (C) MS-based quantitative proteomic analysis of the interaction of MSC-complex proteins with constrained-G4 structures (extracted from Supplementary Table 1). Differentially interacting proteins were sorted out using a fold change ≥ 2 and *p*-value < 0.05, allowing to reach a false discovery rate (FDR) inferior to 5% according to the Benjamini-Hochberg procedure.

constrained constructions and western blotting analyses were performed to quantify the relative signal of the bound NELFE protein. As shown in the Fig. 6B, the NELFE protein is significantly enriched on constrained G4 structures relatively to the cyclopeptide (CP-T23) and duplex control 8, indicating a selective binding of the NELFE protein to G4 structures. In addition, hybridization with an anti-NELFA antibody indicates that at least two proteins from the immunoprecipitated NELF complex are found enriched on constrained G4 structures relatively to the controls (CP-T23 and system 8) (data not shown).

**The NELF complex facilitates DNA double-strand breaks induction by Pyridostatin.** In human cells the NELF complex plays an essential role in the RNA-Pol II pausing mechanism⁵³,⁵⁴. RNA-Pol II pausing is a highly controlled mechanism regulating gene expression in eukaryotic cells⁵³. During the pausing state, RNA-Pol II remains tightly associated with nascent RNA molecule in a promoter proximal region. Bioinformatic and pangenomic studies clearly establish that the formation of G4 structures, which are significantly enriched in promoter regions in human cells, correlates with the formation of R-loops and RNA-Pol II pausing sites, two processes associated with RNA-Pol II arrests and the induction of transcription dependent double-stranded DNA breaks (DSBs)³⁰⁻³²,⁵⁵⁻⁶⁰. Pyridostatin (PDS), one of the most selective and potent G4 stabilizers described so far, provokes a rapid induction of transcription dependent DSBs in human cells (⁶¹,⁶² and unpublished data). In order to investigate the impact of the NELF complex on PDS-induced DSBs, we quantified γH2AX signals (DSBs marker) through immunofluorescence studies in NELF proficient and deficient (siRNA-depleted NELF-E) HeLa cells. As shown in Fig. 6C, depletion of the NELFE protein, which also leads to the reduction of other NELF proteins in human cells⁵²,⁶³, provokes a significant reduction of PDS-induced DSBs signals relatively to
<figure><img src="image_13.png" /><figcaption>**Figure 5.** Impact of the orientation and nucleotide composition of connecting loops on the differential enrichment of proteins on G4 structures. Western-blotting analysis and quantification of the interaction of proteins found enriched on constrained G4 structures with modified molecules (1a, 1b, 2, 3, 4) and with unconstrained G4 (c-myc and 21 T). Arrows indicate the 5'-3' strand orientation of single-stranded extensions or connecting loops present on different systems. The modification of the nucleotide composition of connecting loops in the system 4 is indicated by the sequence TCT. Construct 8 and scramble sequence were used as control for pull-down performed with constrained of free-G4 structures, respectively.</figcaption></figure>
A

<table><thead><tr><th rowspan="2">Protein names</th><th rowspan="2">Gene names</th><th rowspan="2">Razor + unique peptides</th><th colspan="3">2 vs 8</th><th colspan="3">1a vs 8</th><th colspan="3">2 vs 1a</th></tr><tr><th>Significant</th><th>p-value</th><th>log2(Fold Change)</th><th>Significant</th><th>p-value</th><th>log2(Fold Change)</th><th>Significant</th><th>p-value</th><th>log2(Fold Change)</th></tr></thead><tbody><tr><td>Negative elongation factor A</td><td>NELFA</td><td>5</td><td>+</td><td>2.37E-04</td><td>3.09</td><td>+</td><td>6.71E-04</td><td>2.67</td><td>4.55E-01</td><td>0.42</td><td></td></tr><tr><td>Negative elongation factor B</td><td>NELFB</td><td>9</td><td>+</td><td>2.04E-05</td><td>2.07</td><td>+</td><td>6.82E-05</td><td>1.78</td><td>2.93E-01</td><td>0.29</td><td></td></tr><tr><td>Negative elongation factor C/D</td><td>NELFCD; TH<sub>1L</sub></td><td>8</td><td>+</td><td>2.34E-06</td><td>3.51</td><td>+</td><td>5.25E-06</td><td>3.19</td><td>3.74E-01</td><td>0.32</td><td></td></tr><tr><td>Negative elongation factor E</td><td>NELFE</td><td>7</td><td>+</td><td>2.22E-07</td><td>3.15</td><td>+</td><td>2.31E-07</td><td>3.13</td><td>9.48E-01</td><td>0.02</td><td></td></tr></tbody></table>

B

<figure><img src="image_14.png" /></figure>

C

<figure><img src="image_15.png" /><figcaption>**Figure 6.** NELF complex interact with G4 structures and modulates the cellular response to G4 ligands. (A) MS-based quantitative proteomic analysis of the interaction of the NELF-complex proteins with constrained-G4 structures (extracted from Supplementary Table 1). Differentially interacting proteins were sorted out using a fold change $\ge$ 2 and *p*-value < 0.05, allowing to reach a false discovery rate (FDR) inferior to 5% according to the Benjamini-Hochberg procedure. (B) Quantification of the interaction of immunoprecipitated Flag NELF-E protein with constrained G4 structures (1a, 2) relative to cyclopeptide (CP-T23) and duplex control (8) constructions. Error bars represent SD from the means, $n \ge$ 3 independent experiments. *p* values were calculated using unpaired t-tests (without corrections for multiple comparisons). ns: *p* > 0.05; \*: *p* < 0.05; \*: \*: *p* < 0.01; \*: \*: *p* < 0.001; \*: \*: *p* < 0.0001. ns non-significant difference. (C) Quantification and representative images $\gamma$H2AX foci fluorescence signal (red) detected HeLa cells transfected with control (Ctrl), or two NELF-E siRNAs (independent sequences) and treated with PDS (20 $\mu$M) for 4 h. Error bars represent SD from the means, $n \ge$ 3 independent experiments. *p* values were calculated using an unpaired multiple Student's *t* test. ns: *p* > 0.05; \*: *p* < 0.05; \*: \*: *p* < 0.01; \*: \*: *p* < 0.001; \*: \*: *p* < 0.0001. signals. Western-blotting analysis of NELF-E depletion in HeLa cells following siRNA treatment is shown.</figcaption></figure>

the control cells. This result indicates that the NELF complex and RNA Pol II pausing favours DSB induction following G4 stabilization.

Discussion

The relationship between G4 structures/motifs and the progressive discovery of proteins that modulate the dynamic of G4 formation, such as helicases or other proteins involved in DNA and RNA transactions has expanded our knowledge to understand the ubiquitous function of G4 structures on cellular metabolism and cell fate<sup>5,6,13,15,17-20,64</sup>. In the present study, we identified through a MS-based quantitative proteomic analysis 425 human proteins significantly enriched on constrained G4 structures relative to a duplex DNA control, in which well-described G4 interacting factors were present. However, we found that some of well-known G4-interacting proteins, such as BLM and WRN, are not enriched on constrained G4 relatively to duplex control, thereby confirming the previously reported non-selective interaction of these proteins with G4s relative to duplex or hairpin DNAs<sup>65</sup>. Interestingly, numerous proteins interacting with constrained G4-DNA structures identified in this
work are involved in RNA transactions: splicing, RNA degradation, transport and RNA surveillance pathways. Moreover, about one hundred proteins enriched from 425 total proteins identified here on constrained G4-DNA structures have been already identified in similar works using RNA G4 structures to trap G4-RNAs binding factors³⁵⁻³⁷. It is noteworthy that most RNA Binding Proteins (RBP) interact with RNAs through OB-fold, RGG, RRM and GARD motifs that have been also involved in the recognition of G4 structures²¹,²²,⁶⁶. Most G4 interacting proteins unfold G4 in order to counteract their impact on DNA and RNA transactions. Interestingly, 62 out of the 425 proteins identified through our approach have also been characterized as playing an important role in the resistance to cytotoxicity induced by G4-ligands⁴⁹, and for many of them our work provides an evidence for a direct or indirect interaction with G4.

An important finding from our results concerns the differential enrichment of G4 interacting proteins on constrained G4 molecules mimicking parallel (systems 1, 6 and 7) versus anti-parallel conformations (systems 2–5). MS-based quantitative proteomic analysis identified 204 proteins significantly enriched on the G4 construction 2 relative to the G4 construction 1, and only 31 proteins with an inverse selectivity. Surprisingly, relative strand orientation and nucleotide composition of terminal connecting loops had no major impact on the binding of proteins to loops containing constructions (2–4), suggesting either a minor role of loops for the interaction with proteins recognizing antiparallel constrained G4, or a role related to their folding into a particular conformation independent on nucleotide composition or strand orientation. In agreement with the second hypothesis, proteins found enriched on system 2–4 relatively to system 1a-b, bind unconstrained G4 structures formed by the c-myc and telomeric sequences. In solution, both sequences have been shown to adopt different topologies and exhibit loops with different nucleotide compositions and strand orientations⁶⁷,⁶⁸. Altogether, these results suggest that while the folding of G-rich sequences into a G4 structure is required for the interaction with identify proteins, the presence of connecting loops is an important determinant of this interaction. In agreement with our data, although selective binding of proteins and small molecules to G4 structures are driven mainly through π-π interactions with tetrad faces, additional contacts with lateral grooves and connecting loops have been shown to stabilize both molecules and protein-G4 complexes²²,⁶⁹. Finally, although the addition of a supplementary nucleotide to construction 1a significantly increases the binding of system 2 enriched proteins on the system 1, we hypothesize that a part of the differential binding of human proteins to loops-containing systems 2–4 relative to loop-free constructions 1, 6 and 7 would be dependent on the particular conformations adopted by connecting loops.

No significant difference was found for 190 out of the 425 proteins interacting with systems 1a and 2, a result that would indicate that these proteins interact with common structural motifs such as tetrad faces or lateral grooves. Although for systems 2–4 the interaction of proteins with tetrad faces could be hindered by lateral loops, future studies with constrained structures containing propeller loops will be needed to fully characterize this impact and the role of lateral grooves.

An intriguing finding from our study is the characterization of the multi-tRNA synthetase complex (MSC) as a G4-interacting complex. The interaction of the MSC complex with G4 structures was confirmed by a recent publication that identified the MSC proteins in a molecular screen identifying RNA G4-interacting factors³⁶. The MSC complex consists of eight Aminoacyl-tRNA Synthetases (ARSs) and three non-enzymatic ARS-interacting multi-functional proteins (AIMP1/p43, AIMP2/p38, and AIMP3/p18) that play an essential role in the protein synthesis by catalysing the activation of amino acids and linking them to their cognate transfer RNAs (tRNAs)⁵⁰. AIMP1/p43 protein is a multifunction protein involved in various physiological and pathological processes. AIMP1 is the precursor of EMAP II, which was released after AIMP1 cleavage⁷⁰. Surprisingly, the C-terminal segment of EMAP II are 50% identical to residues 205–364 of p42 protein from S. cerevisiae that was firstly identified as G4p1 a protein that shows a high and specific affinity for G4 nucleic acids⁷¹.

Finally, the interaction of NELF complex with G4 structures is a major finding of our study. To the best of our knowledge, this is the first time that a physical interaction between the NELF complex and G4 is reported. The NELF complex comprises the NELFA, NELFB, NELFC/D and NELFE proteins. Whether the interaction of NELF with G4 structures is direct or indirect needs to be elucidated. Nevertheless, the NELFE protein contains a RRM-domain (RNA Recognition Motif)⁵² that is present in different RNA binding proteins and interacts with G-rich tracts. These RRM motifs have also been involved in the interaction of several proteins with G4 structures, such as hnRNP A1, and in agreement with our findings, a NELF sub-complex formed by NELFA and NELFE proteins has been shown to selectively bind GC rich sequences⁷². Numerous bioinformatic and genomics analysis have established a strong correlation between G4 and RNA Pol II pausing occurring at the promoter proximal regions³¹,⁷³. Transcriptional pausing consists in the arrest of RNA-Pol II mainly in a promoter proximal position and acts as a transcriptional checkpoint regulating gene expression⁵³. It is induced by the association of two complexes, NELF and DSIF, with RNA-Pol II at the early stages of transcriptional elongation⁵⁴. Although G4 motifs are found at pausing sites³¹, it is still unclear whether G4s act as signals for RNA-Pol II pausing. However, during transcription the impact of G4 on RNA Pol II progression could reflect their capacity to act as physical barriers, to promote the formation of other secondary structures such as R-loops or finally to drive the binding of protein complexes regulating RNA Pol II progression. In human cells, G4 motifs and RNA Pol II pausing have been also associated with transcription-dependent DNA breaks³¹,³²,⁵⁵, and the stabilisation of G4 structures by small molecules has been shown to provoke the formation of double-stranded DNA breaks that are, in part, dependent on RNA-Pol II transcription⁶¹. In this study, we show that the NELF complex regulates the formation of double-stranded DNA breaks induced by pyridostatin, a potent and selective G4 ligand. Altogether, our results suggest that G4 formation could facilitate the binding of the NELF complex to the chromatin to promote RNA Pol II pausing, and thus to act as a mediator of the response to G4 stabilization in human cells.

In this study, we show that biotin-functionalized constrained nucleic acids structures are powerful tools to identify proteins interacting with non-canonical secondary structures such as G4. We validated our approach by the identification of well-known G4 interacting factors and the functional characterization of new protein
complexes related to G4 metabolism. Especially, the identification of NELF complex as a G4 interacting factor establishes a physical link between G4 structures and RNA Pol II pausing mechanism.

Future directions of our approach will concern the construction of constrained G4 structures mimicking intramolecular parallel G4 structures in order to refine the impact of propeller loops on protein binding. Finally, other constrained nucleic acid constructions may represent powerful tools to identify proteins interacting with other non-canonical secondary structures such as the i-motif or R-loops.

## Material and methods

**Constrained G4 and HP synthesis.** The different constrained systems **1–8** were prepared according to the previously reported protocols (42–45). Systems **1a,b**, and **6** with a parallel orientation of the G4-forming strands, were synthetized through a single oxime ligation reaction onto the cyclopeptide of the appropriate DNA or RNA sequences, respectively (45, 43). For system **7**, two successive ligations (i.e. oxime and CuAAC reactions) were used to anchor on the cyclopeptide scaffold, DNA and RNA sequences, respectively (43). Systems **2–5** with an antiparallel orientation of the G4-forming strands, were prepared by anchoring on the cyclopeptide scaffold through the aforementioned ligations the appropriate bi-functionalized oligonucleotides from telomeric sequence (44) and HIV sequence (42). Controls **8** and **CP-T23** were synthetized according to the oxime method by anchoring the appropriate sequences on the cyclopeptide scaffold (45). Biotinylated c-myc (5'GGA-GGG-TGG-GGA-GGG-TGG-GGA-A-TEG-biot), 21T (5' TTA-GGG-TTA-GGG-TTA-GGG-TTA-GGG-TT-TEG-biot) and scramble (5' AAG-TGT-GTG-TGT-GTG-TGT-GTG-TGA-AG-TEG-biot) sequences were purchased from Eurogentec.

The CD spectra of the systems **1–7** which show typical CD signatures of parallel (positive peak at 263 nm and a negative peak around 240 nm for **1a, 1b, 6** and **7**) an antiparallel (two characteristic positive peaks at 242 nm and 294 nm and a negative peak at 262 nm for **2, 3, 4** and **5**) G-quadruplex are depicted in Supplementary Figure S2.

**Cell culture.** HeLa cells were grown in humidified atmosphere with 5% CO2 at 37 °C, in Dulbecco's Modified Eagle Medium (Gibco). Culture medium was supplemented with 10% foetal bovine serum (Eurobio), 100 U/mL penicillin (Gibco) and 100 µg/mL streptomycin (Gibco).

**Plasmid constructions cell transfection and transduction.** The pCMVsVg (envelope plasmid) were kindly provided by E. Gilson. pLPC-puro-N-Flag was a gift from Titia de Lange (Addgene plasmid #12521; http://n2t.net/addgene:12521; RRID: Addgene_12521). The pLPC-puro-N-Flag-NELFE plasmid was obtained by amplification of the human NELFE cDNA (obtained from gene synthesis, GeneArt, LifeTechnologies) using oligonucleotides: Fwd 5'-GACGATGACGATAAAGGATCCTTGGTGATACCCCCCGGACT-3' and Rev- 5'CCCTCTAGATGCATGCTCGAGCTAGAAGCCATCCACAAGGTTTTCC-3' and inserted in the pLPC-puro-N-Flag plasmid between XhoI and BamH1 restriction sites. Retroviral production was performed by transient transfection of HEK-GP2 293 cells (Clontech) with 0.8 µg pCMVsVg and 1.2 µg of pLPC-puro-N-Flag-NELFE plasmids using JetPrime reagent (polyplus). Retroviral particles were recovered from culture supernatant of HEK-GP2 293 cells 48 h and 72 h post transfection. For transduction, fifty thousand Hela cells were plated on six-well plates 24 h prior to transduction. Cell population was selected by puromycin resistance (1 µg/ml). The expression of Flag-NELF-E protein was verified by western blotting with Flag and NELF-E antibodies.

**Pull-down from total protein extract with constrained G4s.** First, proteins were bound to constrained G4s. Briefly, 1 mg of NHEJ protein extract (prepared as previously described³⁶) was incubated with 10 µM of one of constrained G4s (Figure A) in 100 µL of binding buffer (20 mM Hepes (pH 7.5), 50 mM KCl, 0.01% NP40 and 0.5 mM EDTA) for 1.5 h at 4 °C under intermittent shaking (10 s at 1400 rpm every 2 min). During protein-G4 binding step, 1 mL of streptavidin-coupled magnetic beads (Promega, Z5481) per condition were washed three times during 10 min at 4 °C (under intermittent shaking) with binding buffer. Then, 100 µL of protein-G4 binding solution was put onto washed streptavidin-coupled magnetic beads for 30 min at 4 °C under intermittent shaking. After this step, supernatant was stored at –80 °C as "unbound fraction" and beads were washed during 10 min three times at 4 °C (under intermittent shaking) with 0.1% NP40, 150 mM NaCl PBS. Finally, beads were incubated with 100 µL of 0.01% bromophenol blue, 15% glycerol, 2% SDS, 60 mM Tris-HCl (pH 8) for 10 min at 95 °C and supernatant were collected and stored at -80 °C before being used in mass spectrometry and western blotting assays.

**MS-based quantitative proteomic analysis.** Eluted proteins were stacked in a single band in the top of a SDS-PAGE gel (4–12% NuPAGE, Life Technologies) and stained with Coomassie blue R-250 before in-gel digestion using modified trypsin (Promega, sequencing grade) as previously described⁷⁴. Resulting peptides were analyzed by online nanoliquid chromatography coupled to tandem MS (UltiMate 3000 and LTQ-Orbitrap Velos Pro, Thermo Scientific). Peptides were sampled on a 300 µm × 5 mm PepMap C18 precolumn and separated on a 75 µm × 250 mm C18 column (PepMap, Thermo Scientific) using a 120-min gradient. MS and MS/MS data were acquired using Xcalibur (Thermo Scientific).

Peptides and proteins were identified and quantified using MaxQuant (version 1.6.2.10,⁷⁵) using the Uniprot database (*Homo sapiens* reference proteome, October 22nd 2018 version) and the frequently observed contaminant database embedded in MaxQuant. Trypsin was chosen as the enzyme and 2 missed cleavages were allowed. Peptide modifications allowed during the search were: carbamidomethylation (C, fixed), acetyl (Protein N-ter, variable) and oxidation (M, variable). Minimum peptide length was set to 7 amino acids. Minimum number of peptides and razor + unique peptides were set to 1. Maximum false discovery rates—calculated by employing a reverse database strategy—were set to 0.01 at peptide and protein levels. The mass spectrometry proteomics
data have been deposited to the ProteomeXchange Consortium via the PRIDE<sup>76</sup> partner repository with the dataset identifier PXD 021003.

Statistical analysis were performed using ProStaR<sup>77</sup>. Proteins identified in the reverse and contaminant databases, proteins only identified by site, proteins identified with only 1 peptide and proteins exhibiting less than 3 intensity values in one condition were discarded from the list. After log2 transformation, intensity values were normalized by median centering before missing value imputation (slsa algorithm for partially observed values in the condition and DetQuantile algorithm set to first percentile for totally absent values in the condition); statistical testing was conducted using limma test. Differentially interacting proteins were sorted out using a log2 (fold change) cut-off of 1 and a *p*-value cut-off allowing to reach an FDR inferior to 5% according to the Benjamini–Hochberg procedure.

**Immunoprecipitation of Flag-NELF-E and pull-down with constrained G4s.** Four days after seeding, around 30.10<sup>6</sup> Flag-NELF-E expressing HeLa cells were harvested by scrapping in cold PBS. After centrifugation at 2000 g for 5 min at 4 °C, pelleted cells were lysed for 40 min at 4 °C with 1 mL of non-denaturing lysis buffer (50 mM Hepes (pH 7.5), 150 mM NaCl, 0.01% NP40) complemented with 1 mM DTT, 1X Halt Protease Inhibitor Cocktail (ThermoScientific), 23 U/mL benzonase. Cells were then centrifugated at 12,000 rpm at 4 °C for 10 min and protein concentrations in total extract were determined by measuring absorbance at 280 nm (Nanodrop). From total extract, 1 mg of proteins was used to performed incubation overnight at 4 °C with 5 µg of anti-Flag mouse antibody (Sigma) in 250 µL of complemented non-denaturing lysis buffer. Protein A/G-coupled magnetic beads (Pierce, 88802) were washed with complemented non-denaturing lysis buffer and 25 µL of beads were used to immunoprecipitated Flag antibody for 1 h at 4 °C. Supernatant was collected as “flag unbound fraction” and protein concentration was determined as described. Immunoprecipitated proteins were eluted overnight at 4 °C with 50 µg of Flag peptide (Sigma-Aldrich) in 100 µL complemented non-denaturing lysis buffer. Supernatant was collected as “flag eluted fraction” and protein concentration was determined as described. For protein-G4 binding step, 7 µg of flag eluted fraction were incubated with 1 µM of constrained G4s in 100 µL of binding buffer overnight at 4 °C, and pull-down was performed as described above. Finally, beads were incubated with 30 µL of 0.01% bromophenol blue, 15% glycerol, 2% SDS, 60 mM Tris (pH 8) for 10 min at 95 °C and supernatant were collected and stored at -80 °C before being used for western blotting assays.

**Western blotting.** For total pulled-down proteins and for Flag-eluted pulled-down proteins, respectively 45 µL and 10 µL of proteins were loaded for each condition and separated on gradient 4–12% polyacrylamide TGX Stain-Free pre-cast gels (Biorad) and transferred onto nitrocellulose membrane (Biorad). Before blocking (0.1% Tween-20, non-fat dry milk 5% and PBS), UV exposition of membrane was used to confirm homogeneous loading and to quantify transfer signal. The membrane was successively probed with primary antibodies and appropriate goat secondary antibodies coupled to horseradish peroxidase (described in table below). Chemidoc imager (Biorad) was used to perform UV and Clarity ECL (Biorad) detection. Digital data were processed and quantified using ImageLab (Biorad) or ImageJ softwares. Quantifications of antibody signal are relative to transfer signal.

<table><thead><tr><td>Target</td><td>Dilution</td><td>Species</td><td>Class</td><td>Reference</td><td>Manufacturer</td></tr></thead><tbody><tr><td>KU70</td><td>0.2 µg/mL</td><td>Mouse</td><td>Monoclonal</td><td>MA5-13110</td><td>Invitrogen</td></tr><tr><td>eIF4G</td><td>1/1000</td><td>Rabbit</td><td>Polyclonal</td><td>2498</td><td>Cell Signaling Technology</td></tr><tr><td>WRN</td><td>1/2000</td><td>Mouse</td><td>Monoclonal</td><td>W0393</td><td>Sigma-Aldrich</td></tr><tr><td>Mre11</td><td>1/1000</td><td>Mouse</td><td>Monoclonal</td><td>611366</td><td>BD bioscience</td></tr><tr><td>hnRNP A1</td><td>1/1000</td><td>Rabbit</td><td>Polyclonal</td><td>GTX106208</td><td>Gene Tex</td></tr><tr><td>DHX36</td><td>1/1000</td><td>Rabbit</td><td>Polyclonal</td><td>HPA035399</td><td>Sigma-Aldrich</td></tr><tr><td>CNBP</td><td>8/1000</td><td></td><td></td><td></td><td>Kindly provided by N. Calcaterra's group</td></tr><tr><td>NCL23</td><td>1/1000</td><td>Rabbit</td><td>Polyclonal</td><td>ab50279</td><td>Abcam</td></tr><tr><td>NELFE</td><td>1 µg/mL</td><td>Rabbit</td><td>Monoclonal</td><td>A301-913A</td><td>Bethyl laboratories</td></tr><tr><td>NELFB (COBRA)</td><td>1 µg/mL</td><td>Rabbit</td><td>Monoclonal</td><td>A301-911A</td><td>Bethyl laboratories</td></tr><tr><td>Anti-rabbit</td><td>1:10,000</td><td>Goat</td><td>Polyclonal</td><td>111-035-003</td><td>Jackson Immunoresearch</td></tr><tr><td>Anti-mouse</td><td>1:10,000</td><td>Goat</td><td>Polyclonal</td><td>115-035-003</td><td>Jackson Immunoresearch</td></tr></tbody></table>

**RNA interferences.** HeLa cells were seeded at 250.000 cells per well in a 6-wells plate. siRNA (Table) were transfected twice at 50 nM final concentration per well with Lipofectamine RNAiMax Reagent (Invitrogen) according to manufacturer's instructions. Cells were treated and proteic extracts are realized 72 h after the first transfection.

<table><thead><tr><td>Target</td><td>Name</td><td>Sequence</td><td>Manufacturer</td></tr></thead><tbody><tr><td>Luciferase</td><td>siLUC</td><td>5'-CUUACGCUGAGUACUUCGATT-3'</td><td>Eurofins</td></tr><tr><td>NELFE</td><td>siNELFE</td><td>5'-AAGAUGGAGUCAGCAGAUCAG-3'</td><td>Eurofins</td></tr><tr><td>NELFE</td><td>siNELFE_2</td><td>5'-GACCUUCUGGAGAAGAGCUTT-3'</td><td>Eurofins</td></tr></tbody></table>
**Immunofluorescence.** HeLa cells were seeded in 24-wells plate at 100.000 cells/well on glass coverslips (VWR, #631-0150). Twenty-four hours later, HeLa cells were treated with 20 µM pyridostatin (Sigma-Aldrich; CAS number 1085412-37-8) for 4 h, and then washed with PBS and fixed with 2% paraformaldehyde in PBS at room temperature for 10 min, washed with PBS and permeabilized for 15 min at room temperature with 10 mM Tris-HCl pH 7.5, 120 mM KCl, 20 mM NaCl, 0.1% Triton-X 100. Then, cells were washed with PBS and incubated for about 1 h at 37 °C in blocking buffer (20 mM Tris-HCl pH 7.5, 150 mM NaCl, 2% BSA, 0.2% fish gelatin, 0.1% Triton-X 100) prior to incubation overnight at 4 °C with γH2AX (Phospho S139) antibody (Abcam, #81299) diluted at 0.7 µg/mL in blocking buffer. Cells were then washed with 0.1% Tween20-PBS and incubated with secondary goat anti-rabbit antibody coupled to AlexaFluor 488 (Thermo Fisher Scientific) diluted at 2 µg/mL in blocking buffer for 1 h at room temperature. At last, cells were washed with 0.1% Tween20-PBS and stained with 0.1 µg/mL DAPI for 20 min at room temperature, and coverslips were mounted with Vectashield mounting medium (Vector Laboratories). Nuclear γH2AX foci staining overlapping with DAPI staining were quantified with ImageJ software. Quantifications of nuclear γH2AX foci induced by pyridostatin are represented normalized to non-treated (NT) conditions.

**KEGG pathway, gene ontology.** The ClueGO v2.3.3<sup>78,79</sup> plugin for Cytoscape<sup>80</sup> (v3.8) was used to determine networks of enriched KEGG pathways and Gene Ontology terms (Biological Process and Molecular Function). A right-sided (Enrichment) test based on the hyper-geometric distribution was performed on the corresponding Entrez gene IDs for each gene list and the Bonferroni adjustment (*p*<0.05) was applied to correct for multiple hypothesis testing. The Kappa-statistics score threshold was set to 0.4 and GO term fusion was used to diminish redundancy of terms shared by similar proteins. Other parameters include: GO level intervals (3–8 genes) and Group Merge (50%).

**Statistical analysis.** All results provide from at least three independent experiments. Statistical analyses were performed with GraphPad Prism Software (version 8). For γH2AX quantifications analyses, multiple unpaired t-tests (without corrections for multiple comparisons) were performed between pairs of conditions. On all figures, significant differences between specified pairs of conditions are shown by asterisks (*: *p*-value < 0.05; **: *p*-value < 0.01; ***: *p*-value < 0.0005; ****: *p*-value < 0.0001). ns means non-significant difference.

Received: 6 April 2021; Accepted: 11 June 2021

Published online: 29 June 2021

**References**

1. Riou, J. F., Gomez, D., Lemarteleur, T. & Trentesaux, C. G-quadruplex DNA: myth or reality? *Bull. Cancer* **90**, 305–313 (2003).

2. Prorok, P. et al. Involvement of G-quadruplex regions in mammalian replication origin activity. *Nat. Commun.* **10**, 3274 (2019).

3. Valton, A. L. et al. G4 motifs affect origin positioning and efficiency in two vertebrate replicators. *EMBO J.* **33**, 732–746 (2014).

4. Kim, N. The interplay between G-quadruplex and transcription. *Curr. Med. Chem.* **26**, 2898–2917 (2019).

5. Cammas, A. & Millevoi, S. RNA G-quadruplexes: emerging mechanisms in disease. *Nucleic Acids Res.* **45**, 1584–1595 (2017).

6. Kharel, P., Balaratnam, S., Beals, N. & Basu, S. The role of RNA G-quadruplexes in human diseases and therapeutic strategies. *Wiley Interdiscip. Rev. RNA* **11**, e1568 (2020).

7. Song, J., Perreault, J. P., Topisirovic, I. & Richard, S. RNA G-quadruplexes and their potential regulatory roles in translation. *Translation (Austin)* **4**, e1244031 (2016).

8. Falabella, M., Fernandez, R. J., Johnson, F. B. & Kaufman, B. A. Potential roles for G-Quadruplexes in Mitochondria. *Curr. Med. Chem.* **26**, 2918–2932 (2019).

9. Neidle, S. & Balasubramanian, S. *Quadruplex Nucleic Acids. RSC Biomolecular Sciences* (Cambridge University Press, 2006).

10. Ma, Y., Iida, K. & Nagasawa, K. Topologies of G-quadruplex: Biological functions and regulation by ligands. *Biochem. Biophys. Res. Commun.* **531**(1), 3–17 (2020).

11. Lightfoot, H. L., Hagen, T., Tatum, N. J. & Hall, J. The diverse structural landscape of quadruplexes. *FEBS Lett.* **593**, 2083–2102 (2019).

12. Dai, J., Carver, M. & Yang, D. Polymorphism of human telomeric quadruplex structures. *Biochimie* **90**, 1172–1183 (2008).

13. Hansel-Hertsch, R., Di Antonio, M. & Balasubramanian, S. DNA G-quadruplexes in the human genome: detection, functions and therapeutic potential. *Nat. Rev. Mol. Cell. Biol.* **18**, 279–284 (2017).

14. Murat, P. & Balasubramanian, S. Existence and consequences of G-quadruplex structures in DNA. *Curr. Opin. Genet. Dev.* **25**, 22–29 (2014).

15. Varshney, D., Spiegel, J., Zyner, K., Tannahill, D. & Balasubramanian, S. The regulation and functions of DNA and RNA G-quadruplexes. *Nat. Rev. Mol. Cell Biol.* **21**, 459–474 (2020).

16. Sauer, M. & Paeschke, K. G-quadruplex unwinding helicases and their function in vivo. *Biochem. Soc. Trans.* **45**, 1173–1182 (2017).

17. Mendoza, O., Bourdoncle, A., Boule, J. B., Brosh, R. M. Jr. & Mergny, J. L. G-quadruplexes and helicases. *Nucleic Acids Res.* **44**, 1989–2006 (2016).

18. Maizels, N. G4-associated human diseases. *EMBO Rep.* **16**, 910–922 (2015).

19. Wu, Y. & Brosh, R. M. Jr. G-quadruplex nucleic acids and human disease. *FEBS J.* **277**, 3470–3488 (2010).

20. Brazda, V., Haronikova, L., Liao, J. C. & Fojta, M. DNA and RNA quadruplex-binding proteins. *Int. J. Mol. Sci.* **15**, 17493–17517 (2014).

21. Flynn, R. L. & Zou, L. Oligonucleotide/oligosaccharide-binding fold proteins: A growing family of genome guardians. *Crit. Rev. Biochem. Mol. Biol.* **45**, 266–275 (2010).

22. McRae, E. K. S., Booy, E. P., Padilla-Meier, G. P. & McKenna, S. A. On characterizing the interactions between proteins and guanine quadruplex structures of nucleic acids. *J. Nucleic Acids Res.* **2017**, 9675348 (2017).

23. Takahama, K. et al. G-Quadruplex DNA- and RNA-specific-binding proteins engineered from the RGG domain of TLS/FUS. *ACS Chem. Biol.* **10**, 2564–2569 (2015).

24. Gonzalez, V., Guo, K., Hurley, L. & Sun, D. Identification and characterization of nucleolin as a c-myc G-quadruplex-binding protein. *J. Biol. Chem.* **284**, 23622–23635 (2009).
25. Ishiguro, A., Kimura, N., Watanabe, Y., Watanabe, S. & Ishihama, A. TDP-43 binds and transports G-quadruplex-containing mRNAs into neurites for local translation. *Genes Cells* **21**, 466–481 (2016).

26. Chambers, V. S. et al. High-throughput sequencing of DNA G-quadruplex structures in the human genome. *Nat. Biotechnol.* **33**, 877–881 (2015).

27. Hansel-Hertsch, R., Spiegel, J., Marsico, G., Tannahill, D. & Balasubramanian, S. Genome-wide mapping of endogenous G-quadruplex DNA structures by chromatin immunoprecipitation and high-throughput sequencing. *Nat. Protoc.* **13**, 551–564 (2018).

28. David, A. P. et al. CNBP controls transcription by unfolding DNA G-quadruplex structures. *Nucleic Acids Res* **47**, 7901–7913 (2019).

29. Fleming, A. M., Ding, Y. & Burrows, C. J. Oxidative DNA damage is epigenetic by regulating gene transcription via base excision repair. *Proc. Natl. Acad. Sci. U.S.A.* **114**, 2604–2609 (2017).

30. De Magis, A. et al. DNA damage and genome instability by G-quadruplex ligands are mediated by R loops in human cancer cells. *Proc. Natl. Acad. Sci. U.S.A.* **116**, 816–825 (2019).

31. Eddy, J. et al. G4 motifs correlate with promoter-proximal transcriptional pausing in human genes. *Nucleic Acids Res.* **39**, 4975–4983 (2011).

32. Puget, N., Miller, K. M. & Legube, G. Non-canonical DNA/RNA structures during transcription-coupled double-strand break repair: Roadblocks or bona fide repair intermediates?. *DNA Repair (Amst.)* **81**, 102661 (2019).

33. Kotsantis, P. et al. RTEL1 regulates G4/R-loops to avert replication-transcription collisions. *Cell Rep.* **33**, 108546 (2020).

34. Miglietta, G., Russo, M. & Capranico, G. G-quadruplex-R-loop interactions and the mechanism of anticancer G-quadruplex binders. *Nucleic Acids Res.* **48**, 11942–11957 (2020).

35. Herdy, B. et al. Analysis of NRAS RNA G-quadruplex binding proteins reveals DDX3X as a novel interactor of cellular G-quadruplex containing transcripts. *Nucleic Acids Res.* **46**, 11592–11604 (2018).

36. Herviou, P. et al. hnRNP H/F drive RNA G-quadruplex-mediated translation linked to genomic instability and therapy resistance in glioblastoma. *Nat. Commun.* **11**, 2661 (2020).

37. Serikawa, T. et al. Comprehensive identification of proteins binding to RNA G-quadruplex motifs in the 5' UTR of tumor-associated mRNAs. *Biochimie* **144**, 169–184 (2018).

38. Vlasenok, M. et al. Data set on G4 DNA interactions with human proteins. *Data Brief* **18**, 348–359 (2018).

39. von Hacht, A. et al. Identification and characterization of RNA guanine-quadruplex binding proteins. *Nucleic Acids Res.* **42**, 6630–6644 (2014).

40. Williams, P., Li, L., Dong, X. & Wang, Y. Identification of SLIRP as a G quadruplex-binding protein. *J. Am. Chem. Soc.* **139**, 12426–12429 (2017).

41. Burge, S., Parkinson, G. N., Hazel, P., Todd, A. K. & Neidle, S. Quadruplex DNA: sequence, topology and structure. *Nucleic Acids Res.* **34**, 5402–5415 (2006).

42. Bonnat, L. et al. Template-mediated stabilization of a DNA G-quadruplex formed in the HIV-1 promoter and comparative binding studies. *Chemistry* **23**, 5602–5613 (2017).

43. Bonnat, L. et al. Templated formation of discrete RNA and DNA:RNA hybrid G-quadruplexes and their interactions with targeting ligands. *Chemistry* **22**, 3139–3147 (2016).

44. Bonnet, R., Lavergne, T., Gennaro, B., Spinelli, N. & Defrancq, E. Construction of anti-parallel G-quadruplexes through sequential templated click. *Chem. Commun. (Camb.)* **51**, 4850–4853 (2015).

45. Murat, P. et al. A novel conformationally constrained parallel g quadruplex. *ChemBioChem* **9**, 2588–2591 (2008).

46. Buck, D. et al. Cernunnos, a novel nonhomologous end-joining factor, is mutated in human immunodeficiency with microcephaly. *Cell* **124**, 287–299 (2006).

47. Walker, J. R., Corpina, R. A. & Goldberg, J. Structure of the Ku heterodimer bound to DNA and its implications for double-strand break repair. *Nature* **412**, 607–614 (2001).

48. Mishra, S. K., Tawani, A., Mishra, A. & Kumar, A. G4IPDB: A database for G-quadruplex structure forming nucleic acid interacting proteins. *Sci Rep* **6**, 38144 (2016).

49. Zyner, K. G. et al. Genetic interactions of G-quadruplexes in humans. *Elife* **8**, e46793 (2019).

50. Mirande, M. The Aminoacyl-tRNA synthetase complex. *Subcell Biochem.* **83**, 505–522 (2017).

51. Szklarczyk, D. et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. *Nucleic Acids Res.* **47**, D607–D613 (2019).

52. Narita, T. et al. Human transcription elongation factor NELF: identification of novel subunits and reconstitution of the functionally active complex. *Mol. Cell Biol.* **23**, 1863–1873 (2003).

53. Core, L. & Adelman, K. Promoter-proximal pausing of RNA polymerase II: a nexus of gene regulation. *Genes Dev.* **33**, 960–982 (2019).

54. Yamaguchi, Y., Shibata, H. & Handa, H. Transcription elongation factors DSIF and NELF: promoter-proximal pausing and beyond. *Biochim. Biophys. Acta* **1829**, 98–104 (2013).

55. Chen, L. et al. R-ChIP using inactive RNase H reveals dynamic coupling of R-loops with transcriptional pausing at gene promoters. *Mol. Cell* **68**, 745–757 (2017).

56. Kim, J. J. et al. Systematic bromodomain protein screens identify homologous recombination and R-loop suppression pathways involved in genome integrity. *Genes Dev.* **33**, 1751–1774 (2019).

57. Aguilera, A. & Garcia-Muse, T. R loops: from transcription byproducts to threats to genome stability. *Mol. Cell* **46**, 115–124 (2012).

58. Aguilera, A. & Gomez-Gonzalez, B. DNA-RNA hybrids: the risks of DNA breakage during transcription. *Nat. Struct. Mol. Biol.* **24**, 439–443 (2017).

59. Singh, S. et al. Pausing sites of RNA polymerase II on actively transcribed genes are enriched with DNA double-stranded breaks. *J. Biol. Chem.* **295**(12), 3990–4000 (2020).

60. Sollier, J. & Cimprich, K. A. Breaking bad: R-loops and genome integrity. *Trends Cell Biol.* **25**, 514–522 (2015).

61. Rodriguez, R. et al. Small-molecule-induced DNA damage identifies alternative DNA structures in human genes. *Nat. Chem. Biol.* **8**, 301–310 (2012).

62. Rodriguez, R. et al. A novel small molecule that alters shelterin integrity and triggers a DNA-damage response at telomeres. *J. Am. Chem. Soc.* **130**, 15758–15759 (2008).

63. Sun, J. et al. Deregulation of cofactor of BRCA1 expression in breast cancer cells. *J. Cell Biochem.* **103**, 1798–1807 (2008).

64. Zell, J., Rota Sperti, F., Britton, S. & Monchaud, D. DNA folds threaten genetic stability and can be leveraged for chemotherapy. *RSC Chem. Biol.* **2**, 47–76 (2020).

65. Kamath-Loeb, A., Loeb, L. A. & Fry, M. The Werner syndrome protein is distinguished from the Bloom syndrome protein by its capacity to tightly bind diverse DNA structures. *PLoS ONE* **7**, e30189 (2012).

66. Hentze, M. W., Castello, A., Schwarzl, T. & Preiss, T. A brave new world of RNA-binding proteins. *Nat. Rev. Mol. Cell Biol.* **19**, 327–341 (2018).

67. Ambrus, A. et al. Human telomeric sequence forms a hybrid-type intramolecular G-quadruplex structure with mixed parallel/antiparallel strands in potassium solution. *Nucleic Acids Res.* **34**, 2723–2735 (2006).

68. Ambrus, A., Chen, D., Dai, J., Jones, R. A. & Yang, D. Solution structure of the biologically relevant G-quadruplex element in the human c-MYC promoter. Implications for G-quadruplex stabilization. *Biochemistry* **44**, 2048–2058 (2005).
69. Haider, S. M., Neidle, S. & Parkinson, G. N. A structural analysis of G-quadruplex/ligand interactions. *Biochimie* **93**, 1239–1251 (2011).

70. Zhou, Z., Sun, B., Huang, S., Yu, D. & Zhang, X. Roles of aminoacyl-tRNA synthetase-interacting multi-functional proteins in physiology and cancer. *Cell Death Dis.* **11**, 579 (2020).

71. Frantz, J. D. & Gilbert, W. A novel yeast gene product, G4p1, with a specific affinity for quadruplex nucleic acids. *J. Biol. Chem.* **270**, 20692–20697 (1995).

72. Vos, S. M. et al. Architecture and RNA binding of the human negative elongation factor. *Elife* **5**, e14981 (2016).

73. Szlachta, K. et al. Alternative DNA secondary structure formation affects RNA polymerase II promoter-proximal pausing in human. *Genome Biol.* **19**, 89 (2018).

74. Salvetti, A. et al. Nuclear functions of nucleolin through global proteomics and interactomic approaches. *J. Proteome Res.* **15**, 1659–1669 (2016).

75. Tyanova, S., Temu, T. & Cox, J. The MaxQuant computational platform for mass spectrometry-based shotgun proteomics. *Nat. Protoc.* **11**, 2301–2319 (2016).

76. Perez-Riverol, Y. et al. The PRIDE database and related tools and resources in 2019: improving support for quantification data. *Nucleic Acids Res.* **47**, D442–D450 (2019).

77. Wieczorek, S. et al. DAPAR & ProStaR: software to perform statistical analyses in quantitative discovery proteomics. *Bioinformatics* **33**, 135–136 (2017).

78. Bindea, G., Galon, J. & Mlecnik, B. CluePedia Cytoscape plugin: pathway insights using integrated experimental and in silico data. *Bioinformatics* **29**, 661–663 (2013).

79. Bindea, G. et al. ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. *Bioinformatics* **25**, 1091–1093 (2009).

80. Shannon, P. et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. *Genome Res.* **13**, 2498–2504 (2003).

## Acknowledgements

This work was funded by grants from ANR (ANR-16-CE11-0006-01), and La Ligue Nationale Contre le Cancer (Equipe Labellisée 2018). Patrick Calsou is a researcher from INSERM. Proteomic experiments were partly supported by the ProFi grant (ANR-10-INBS-08-01). The Nanobio-ICMG platform (UAR 2607) is acknowledged for providing facilities for the synthesis and purification of oligonucleotides (R. Lartia) as well as for mass spectrometry analyses (A. Durand, L. Fort and R. Guéret).

## Author contributions

D.G., J.F.R. and E.D. proposed the idea A.P., A.D., A.A., Y.C. and T.L. performed experiments. D.G., J.F.R., E.D., T.L., Y.C., S.B. and P.C., wrote the initial draft. All authors confirmed the final manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-021-92806-8.

**Correspondence** and requests for materials should be addressed to D.G.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_16.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or

format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021