Marianne Duquenne

Hélène Prost

Joachim Schöpfel

Bioeconomy and Grey Literature -Empirical Evidence

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.

This work is part of a larger, interdisciplinary research project funded by the French National Research Agency (RECABIO ANR-16-IDEX-0004 ULNE).

Acknowledgment

The purpose of this poster is to provide empirical evidence on the role of grey literature in bioeconomy, as one of the emergent research fields at the crossroads of several disciplines and with high societal and industrial impact.

Abstract Methodology

Based on scientometric studies and on results from a survey in the field with research laboratories, corporate companies and other organizations.

QUALITY, VISIBILITY AND RETRIEVAL OF GREY LITERATURE ARE (STILL) ISSUES.

Conclusions

Assignment of persistent identifiers and description with rich metadata

Improving FINDABILITY

Recommendations

Raising AWARENESS

Creating TRUST

Results

• Grey literature on bioeconomy is more or less invisible in bibliographic databases, except for conference papers ; • Grey literature on bioeconomy, especially reports, preprints and PhD theses, is more visible in open repositories ; • Identification and retrieval of grey literature on bioeconomy with the usual discovery tools is not easy ; • Most respondents were rather reluctant to use grey literature (reports, preprints, PhD theses). The perceived advantage is novelty and diversity ; but they are concerned about poor (uncertain) quality.

Mass of information «

We have to look elsewhere for information, that is to say all the grey literature, reports and activity reports, databases.

L E A R N M O R E A B O U T T H I S R E S E A R C H

The purpose of our research is the analysis of the production, the usage and the communication of scientific and technical information in the field of bioeconomy and in the context of open science, in the French region Hauts-de-France.

Research project Bioeconomy ?

Bioeconomy is an "economy based on the sustainable production and conversion of renewable biomass into a range of bio-based products, chemicals and energy".

Figure 2 .Figure 1 .Figure 3 .Figure 4 .

 2134 Figure 2. Bioeconomy stakeholders interactive map

(1

  Marianne Duquenne 1 , Hélène Prost 2 , Joachim Schöpfel 1 Univ. Lille, EA 4073 -GERiiCO -Groupement d'Etudes et de Recherche Interdisciplinaire en Information et Communication, F-59000 Lille, France. 2 CNRS, EA 4073 -GERiiCO -Groupement d'Etudes et de Recherche Interdisciplinaire en Information et Communication, F-59000 Lille, France.


Yoan Monneau

Cyrielle Arrault

Coraline Duroux

Marie Martin

Fabien Chirot

Luke Mac Aleese

Marion Girod

Clothilde Comby-Zerbino

Agnès Hagège

Olivier Walker

Maggy Hologne

Structural and dynamical insights into SilE silver binding from combined analytical probes

ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

Introduction

Silver-based antiseptics have a wide application in health care due to their high availability, broad spectrum, and safety for humans. [1][2][3] However, silver-resistant Gram-negative bacteria were reported and the resistance mainly relies on the presence of a plasmid-encoded specific silver resistance sil operon 4 which is homologous to the copper resistance Cus operon. Both operons include the tripartite efflux pump CBA that acts as a proton antiport, the RS system that senses the metal present in the periplasm and enhances the operon gene expression accordingly, and a metal chaperon F, that routes metal ions to the efflux pump. The SilE protein does not have any counterpart in the Cus operon but has a homolog in E. coli, PcoE, that displays 48% identity in sequence with SilE. A previous work has reported that the binding of PcoE with Cu + and Ag + induces its dimerization and partial structural folding while it has been shown that PcoE can bind up to 6 equivalents of Ag + before precipitation. 5 Gene deletion studies show a complete loss of bacterial silver resistance capacity when the SilE gene is deleted along with the Cus operon. 6 SilE has its own promoter and is early expressed while bacteria are exposed to silver. 4 Considering that new silver-based drugs are under development to fight multi-resistant bacteria 1,7,8 , the proper design of such new drugs relies on our understanding of the molecular mechanism underlying silver resistance. Thus, providing insights into silver-bound SilE structure and dynamics will improve our knowledge of the role of this essential sil operon component concerning the efflux pump machinery. In a previous work of Asiani et al, SilE was described as an intrinsically disordered protein (IDP) that acts as a "molecular sponge" dedicated to the early sequestration of silver ions. 9 Nevertheless, these previous studies concerning SilE and PcoE lack a clear identification of the different silver binding sites at the atomic level and no experimental findings regarding a possible tertiary structure have been reported. Recently, we have improved our knowledge in the field by deriving NMR (Nuclear Magnetic Resonance) structures of SilE mimicking peptides and shown that sequences encompassing two successive motifs HxxM or MxxH fold into helices in the presence of silver ions. 10 In the same order of idea, we have determined low micromolar dissociation constants from NMR titration analysis. From the vast landscape of biophysical techniques dedicated to the study of IDPs, NMR spectroscopy is particularly well suited especially in the case of complexes where disorder is at least partly retained. [START_REF] Prestel | Intrinsically Disordered Proteins[END_REF] Structure and dynamics can be explored at the atomic level even for low-affinity complexes or low-populated states. Although NMR is powerful on its own, the efforts toward our detailed understanding of the conversion between fully IDP to partially disordered complexes introduce significant technical challenges making an integrative structural biology approach necessary. For instance, the stoichiometry for multisite ligand binding is often difficult to assess by NMR. Moreover spectral quality can significantly decrease due to dramatic line broadening when interconversion rates are at the time regime of chemical shifts. [START_REF] Ferreon | [END_REF] In the present work, we have advanced our knowledge regarding the structural organization of SilE compared to the previous studies. Therefore, we propose an original combination of analytical techniques to characterize the structure and dynamics of silver binding to SilE, including Nuclear Magnetic Resonance (NMR) and Circular Dichroism (CD) for local structure and dynamics characterization; Small-Angle X-Ray Scattering (SAXS) and Ion-Mobility-Mass-Spectrometry (IM-MS) for global shape analysis; Liquid Chromatography High-Resolution Mass-Spectrometry (LC-HRMS) and Capillary-Electrophoresis Inductively-Coupled-Plasma Mass-Spectrometry (CE-ICP-MS) to assess the stoichiometry of the metal-protein complexes. In this regard, we will show how these complementary techniques allow the characterization of the SilE binding sites at an atomic level. An assessment of the stoichiometry of the complexes formed but also the structural modifications that SilE undergoes when complexed with Ag + will be provided. In a final step, we will use a combination of AlphaFold structural prediction, SAXS, and NMR spin relaxation measurements to propose for the first time an insightful structure of the silverbound SilE protein. Overall, this ensemble of results will allow a deeper overview of the role of SilE in the general context of the silver efflux pump machinery.

Results and discussion

The SilE sequence (Fig. S1) used in this work was defined according to the wild-type sequence of the Salmonella Typhimurium SilE (UniProt accession number Q9Z4N3). The first twenty residues of the signal peptide that routes SilE to the periplasm have been removed. Therefore, the SilE construct has been successfully expressed and isotopically labeled with E. coli cells. As indicated in blue in the sequence, SilE is composed of 6 HxxM, 2 MxxH, and one HxM motifs. As previously described by Chabert et al, 10,13 silver ions can bind His and Met residues and association binding constants have been calculated for each motif.

Silver-bound SilE complexes exist at different stoichiometries

As indicated above, the SilE binding sites are widely distributed over its amino acid sequence from M56 to G143, with a minor site around the position M40. Since the global stoichiometry of the silver-bound SilE complex needs to be determined, mass spectrometry is particularly well suited for such analysis. Solutions of SilE containing 0 to 10 equivalents of Ag + relative to SilE, were analyzed in positive electrospray ionization highresolution mass spectrometry (ESI-HRMS). The same stock solution of silver nitrate (oxidation state I) has been used throughout all titration studies presented hereafter and the accuracy of the Ag + solution concentration was assessed before any assays using a calibration curve and silver ion selective electrode. SilE protein was observed with a charge state distribution from 6 + to 21 + protonated adducts (Fig. S2A, SilE in the free state) for HR-MS and IM-MS experiments. With an increase in the concentration of Ag + solution, protonated [SilE:nAg + ] complex ions were also detected. As an example, Fig. S2B shows the ESI-MS spectrum of SilE with 5 equivalents of Ag + solution, where SilE protonated ions, as well as peaks corresponding to [SilE + nAg + (z-n)H] z+ complexes with n = 1 to 5, were detected (see insert in Fig. S2B). Some hybrid adducts with Ag + and 1 to 3 Na + ions were also observed in the MS spectra. To simplify and better observe the different stoichiometries of [SilE:nAg + ] complexes, all the recorded ESI-MS spectra were deconvoluted. Deconvoluted spectra, averaged over the end of the elution band in order to avoid chemical artefacts, are presented in Fig. 1A for the different molar equivalents of Ag + in solution. The number of Ag + bound to one molecule of SilE increases with the concentration of Ag + present in the solution. A high abundance of 2 to 6 Ag + per SilE molecule was observed, while complexes with 7 or 8 Ag + were only observed at low abundance (even with a higher concentration of Ag + ). However, higher concentrations of Ag + did not allow for the observation of more Ag + complexes, suggesting the saturation of the SilE complex sites. More quantitatively, the relative abundance of the different complexes [SilE-nAg + ] were determined by integration of the deconvoluted spectra over the full elution band and plotted in Fig. 1B as a function of the molar equivalents Ag + in solution. The free SilE protein is the most intense species detected for 1 to 6 equivalents Ag + in the solution (Fig. 1A). However, starting from 4 equivalents of Ag + , the distribution shifts significantly towards the higher stoichiometries. Remarkably, only little 1:2 / 1:3 and 1:4 complexes seem to be formed before 3, 4, and 5 equivalents of Ag + are respectively added to the solution (while 1:5 and 1:6 already appear at 5 and 6 equivalents). After 4 molar equivalents of Ag + in the solution, the abundances of the 1:1 complex starts decreasing while the 1:2 and 1:3 soon reach a plateau and the 1:4 complex progressively takes over. Up to 10 equivalents of Ag + , the abundance of the 1:4 complex remains maximal, pointing towards 4 strong complexations sites.

Please do not adjust margins

Please do not adjust margins

To further confirm the observed stoichiometry in ESI-HRMS, the samples were analyzed using CE-ICP-MS. From 1 to 4 equivalents of Ag + in solution, the quantity of Ag + bound to SilE linearly increases (Fig. 1C,D). The nearly quantitative formation of the complex indicates a high affinity of SilE for Ag + . However, for higher concentrations of Ag + , no increase of the silver bound to SilE was detected but only free Ag + (Fig. 1D). This result indicates that either we have reached the saturation of SilE binding sites or additional labile sites exist. In the latter case, the potential sites could not be detected by using this CE-ICP-MS strategy. Indeed, the CE separation relies on the electric field applied along the capillary and the analyte must join the anode to be analyzed while silver ions are positive. We hypothesized that labile sites are prompted to dissociate while the migration front is quickly depleted from silver. To challenge the potential presence of such labile sites, the same experiments have been performed but with increasing electrolyte Ag + concentrations (1 -3 µM). Results are displayed in Fig. 1E,F. Increasing the electrolyte Ag + concentration led to an increase of the SilE-bound Ag + peak areas until a Please do not adjust margins Please do not adjust margins concentration close to 1 µM. For the higher Ag + electrolyte concentrations, the [SilE:nAg + ] peak area remains constant, which reflects the saturation of the protein. Quantification of the signal at the plateau was performed by external calibration and revealed the complexation of 7-8 Ag + per protein in these non-dissociating conditions. These results agree with those of our ESI-HRMS studies and further corroborate the maximum number of Ag + that can specifically bind to SilE. As a consequence, it can be assumed that SilE contains 4 kinetically inert complexation sites, which can be seen by directly injecting the complex, and 3-4 more labile sites for Ag + , which require an equilibrium state all along the separation to be detected.

SilE is mainly disordered in the free state

Using 3D standard NMR experiments, the backbone assignment of SilE has been performed (HN, Ha, Ca, Cb, and N). The 2D 1 H-15 N HSQC NMR spectrum including amino acid numbering is displayed in Fig. 2A.

ARTICLE

Please do not adjust margins Please do not adjust margins Resonances are dispersed in a narrow spectral range from 8 to 8.7 ppm, which is the hallmark of an intrinsically disordered protein. 14 The observed chemical shifts result from the average of a conformational ensemble the protein may adopt if the exchange time scale is much faster than the chemical shift time scale. Consequently, the propensity to form different secondary structures can be assessed from the analysis of the backbone chemical shifts. 15 For the amino acid regions 76-88, 112-117, and 132-139, we observe a propensity for transient a-helix structure between 10% and 27% while the rest of the SilE sequence appeared in a fully extended conformation (Fig. 2C). To support these findings, we have set up NMR spin relaxation experiments that are closely related to the global and local motion of the protein.

For the region that exhibits a propensity for a-helix, the mean value of R2 reaches 5.1 s -1 while it is 3.8 s -1 for the rest of the sequence (excluding N-and C-ter). Additionally, the average hetNOEs show a value of 0.30 and 0.1 for the same regions respectively. The latter one is a typical value for IDP. 16 This effect is less pronounced in R1 experiments due to its higher sensitivity for local motion compared to global motion. The HetNOEs are mainly dependent on the local N-H bond rigidity or flexibility, while the R2 are also affected by the global motion. The N-and C-termini present highly negative hetNOE values indicating high flexibility. This corroborates that more rigid conformations are sampled by the regions 76-88, 112-117, and 132-139 compared to the purely intrinsically disordered behavior of the rest of the protein (Fig. 2D).

SilE has several silver binding sites

To investigate the silver binding areas, we have performed NMR titration experiments, recording a series of 1 H-15 N HSQC spectra on 15 N-SilE by gradually increasing the Ag + solution concentration (Fig. S3). Along the titration experiments, we can make different observations that correspond to different exchange regimes. First and foremost, a set of 40 peaks experienced a slow exchange rate. Second, 33 other peaks exhibit a very strong decrease of their signal intensity so that they completely disappeared from the different spectra in saturating conditions. At this level, we can rule out a paramagnetic effect due to Ag + binding. Indeed, Ag + has an electron configuration of 5s 0 4d 10 so that it does not show any unpaired electron and paramagnetic properties. 17 Consequently, signal disappearance may be ascribed to intermediate exchange regime or oligomerization. 18 This statement is also corroborated by the very strong increase of the transverse relaxation rate in this region (see below) and is likely due to the significant contribution of Rex to the global R2. Therefore, we can rule out any oligomerization effect in the present case. The expected number of peaks was never recovered even with the addition of a large excess of Ag + (up to 50 eq). As can be seen on Fig S4, any attempt to recover these peaks was also unsuccessful either by modifying magnetic field (600 and 900 MHz), temperature (from 283 to 323 K), pH (5.2 to 7.8), ionic strength (20 to 300 mM NaF) or pressure (High pressure NMR experiments up to 2250 bars). Therefore, we have assigned chemical shifts of SilE in silver saturating conditions (HN, Ha, Ca, Cb, and N). All visible 1 H- 15 N correlations correspond to residues 19-58, 76, and 93-143 (Fig. 2B). Unfortunately, residues 59 to 92 are then out of reach for studying binding, structures, and dynamics by NMR. It is noteworthy that both spectra in the free and silver-bound states are in stark contrast compared to the ones previously published where neither assignment nor structural information have been derived. 9 The NMR chemical shift perturbations (CSP) upon the addition of 6 equivalents of Ag + were calculated for each amino acid (Fig. 3). From S104 to A142, free and bound states are visible due to the slow exchange regime. The proportion of the bound state can be estimated from the NMR signal ratio between both states (Fig. S5). For 1 equivalent of Ag + , 35 % of the SilE region encompassing residues S104 to A142 is in its bound conformation and reaches 91 % for 6 equivalents of Ag + . No significant changes occur for 6 to 9 equivalents of Ag + . Large variations of CSP are observed for residues E110 to N122 and residues H129 to S141. In these two regions, several His (H111, H118, H129, H136) and Met (M121, M132, M139) residues are present and their associated perturbations are in good agreement with the perturbations seen on the SilEmimicking peptides. 10 Additionally, we have detected measurable shifts in the region H80 to M90 and M59 to M72 before signal disappearance where 3 HxxM and 1 MxxH motifs are present. The latter one ( 59 MDQH 63 ) harbors visible signals Please do not adjust margins Please do not adjust margins until 4 equivalents of Ag + while the others disappeared after the addition of 2 equivalents of Ag + (Fig. S5B). It is also noteworthy that the 59 MDQH 63 motif is involved in a less stable helical formation. 9 The motif H38 to M40 also binds silver but displays significant CSPs from 5 to 9 equivalents of Ag + and their corresponding signals are in a fast exchange regime. Overall, we have identified two major areas of SilE that bind silver ions with different dynamics, one in a slow exchange regime, and the second one in an intermediate exchange regime. The third one is the weakest area of binding and has fast exchange regime dynamics. All the identified binding motifs of the entire SilE protein are in perfect agreement with the ones found in previous peptide-mimicking SilE 10 and suggest that (i) the individual SilE-mimicking peptides may be used individually to describe the SilE/Ag + interaction and (ii) each helical binding motif can be treated as independent segments in the SilE protein.

Although a general analysis of CSP is widely used to describe the protein/ligand binding interface, the discrimination between direct binding and a local conformational rearrangement upon ligand binding may be problematic. This is particularly true for IDP proteins which have no defined structure and possibly fold to stabilize their ligands. 19 To further confirm this binding site mapping onto the SilE sequence, we have performed a differential digestion analysis monitored by mass spectroscopy. The SilE protein in the free state or in presence of 6 equivalents of Ag + was digested and the resulting peptides were analyzed by LC-MS/MS. Among the digested peptides, 3 peptides containing the 59 MDQH 63 , 69 HETM 72 , 108 MNEH 111 , 118 HEFM 121 , and 129 HQAM 132 motifs were detected as protonated adducts from the digested solution of SilE (Table S1A). When the digestion was performed on SilE protein in interaction with Ag + , only the peptides containing the 59 MDQH 63 , 69 HETM 72 , and 108 MNEH 111 motifs were detected (Table S1B). This result could indicate that the strong binding between silver ions and the two motifs 118 HEFM 121 and 129 HQAM 132 may disturb the digestion process due to a particular structural organization (or would not allow to observe the protonated form of the peptide). On the contrary, 59 MDQH 63 , 69 HETM 72 , and 108 MNEH 111 motifs are observed in protonated peptides in both conditions and do not perturb the digestion process when Ag + is present.

Silver-bound SilE displays secondary structure folds

A series of CD spectra were recorded for the free SilE and in the presence of 1 to 8 equivalents of silver ions relative to SilE (Fig. 4A). The ellipticity significantly decreases at 222 nm and increases at 190 nm upon the addition of up to 4 equivalents of Ag + in solution. These data demonstrate the formation of ahelices upon the addition of silver. We can also point out the presence of a clear isodichroic point at 206 nm indicating a twostep process between the free-form and the silver-bound form for SilE. From 5 to 8 equivalents of Ag + in solution, the ellipticity only slightly changes toward more a-helix folds. This structuration is further confirmed by NMR since new peaks appeared out of the 8 to 8.7 ppm area of the 1 H-15 N HSQC spectra upon the addition of 2 to 9 equivalents of Ag + in solution, which is a strong indicator of a protein structuration. The secondary structure propensities 15 calculated from backbone chemical shifts assignment for each residue of the silver-bound SilE in presence of 6 equivalents of Ag + in solution further confirmed the a-helix folding upon silver addition (Fig. 4B). The folded area encompasses more residues compared to the free form (segments 110-123 + 128-142 vs 112-117 + 132-139 without silver) and the a-helix propensity increases from 18 % in the free state up to 65 % in the bound state. To gain insight into the structuration of the SilE part that undergoes strong NMR signal attenuation, we designed a truncated SilE 57-95 . A series of CD experiments were performed and confirmed the formation of a-helix upon silver binding for this shorter construct (Fig. S6) as observed for the full-length SilE. We can also point out the presence of a clear isodichroic point at 206 nm indicating a two-step process between the apoform and the silver-bound form for SilE 57-95 . To get further insight into the structural changes and their diversity as a function of the stoichiometry of the complexes, ion mobility coupled to mass spectrometry (IM-MS) spectra were acquired using solutions with 0 to 12 equivalents of silver ions relative to the SilE protein. Two factors lead to the sorting of the [SilE:nAg + ] complexes in IM-MS (with n the number of bound Ag + ): the mass of the complexes and their arrival times. The latter can be linked to the overall shape of the complexes. We detect species from n = 0 to 6, which is consistent with HRMS data, considering that low-abundant complexes with n = 7 and 8 have not been detected due to the lower sensitivity of the IM-MS instrument compared to the HRMS instrument. The species with n = 1 was not detected with sufficient intensity, whatever the protein:silver ratio in the solution. Arrival time distributions (ATDs) were extracted for the different complexes detected. A selection of ATDs for different stoichiometries is displayed in Fig. 4C. All ATDs are undoubtedly bimodal indicating the coexistence of two ion populations with different mobilities. We interpreted those two populations as two conformational families since no evidence for isobaric SilE oligomers with higher charge states was found in the mass spectra. As generally observed for proteins, each peak is broader than expected from ion diffusion, 20 which indicates that each family is a collection of different conformers with a close structure. 21 The overall shape of the ATD would be compatible with relatively slow interconversion between the two families during the course of the separation due to the plateau between the two peaks. 22 However, no sign of such process was observed when selecting either of the population in a tandem-IMS scheme. 23 Consequently, the two observed peaks can be attributed to the existence of two distinct conformational families, either pre-existing in solution, or more likely resulting from structural reorganization upon desolvation. In the following, we focus on the relative intensity of the main two populations. The ATDs in Fig. 4C, clearly indicate that the relative intensity of the two forms depends on the number of Ag bound in the complex. The peak at lower arrival times corresponds to species with a relatively more compact structure Please do not adjust margins Please do not adjust margins leading to a DT CCSHe value of 1385 Å 2 and the second peak corresponds to an extended form characterized by a DT CCSHe value of 1693 Å 2 . Both forms are nevertheless detected for the bare SilE, with a predominance of the extended form. As the number of bound Ag + increases, the relative population of the compact form increases to become predominant, as visible from Fig. 4D. The data present the average ratio from ATDs for the same ionic species originating from solutions with different relative stoichiometry in Ag + . Importantly, the relative intensity of the two forms does not seem to be affected by the relative concentration of Ag + in the solution, but only depends on the stoichiometry of the complex. Consequently, the observed evolution is not driven by a change in the conformational preferences of SilE due to a change in its environment, but rather reflects a structural change induced by Ag + binding. As already mentioned, the gas phase structures of the complexes probably differ from their original solution structures. However, the observed ATDs can be understood as an indirect mapping of the conformational space sampled by the complexes in solution. 24,25 Following this line, each gas phase conformer corresponds to the endpoint of a distinct portion of the solution-phase conformational space, after conformational relaxation. The observed changes can then be attributed to a modification in the solution-phase conformational space of SilE upon silver binding. In this context, we tentatively interpret the relative increase of the population of the most compact conformer for complexes with an increased silver content as global compaction of the complexes in solution.

ARTICLE

Please do not adjust margins

Please do not adjust margins

This interpretation is supported by small angle scattering data (SAXS) of free SilE and SilE + 6 equivalents of Ag + in solution (Fig. S7). 26 Through a Guinier plot analysis, we derived a radius of gyration Rg = 37.7 ± 0.7 Å for free SilE and Rg = 32.9 ± 0.2 Å for SilE + 6 equivalents Ag + . Therefore, we can conclude the compactness of the structure in presence of silver as observed in IM-MS experiments. At this level, it is important to clarify the meaning of compactness according to IM-MS and SAXS techniques. Compactness is due to the formation of a-helices and is different from a globular structure which refers to a tertiary structure. Interestingly, HRMS enables the observation of an inverse relationship between the number of silver ions observed on SilE and the charge state under which it is detected. This is independent of the concentration of Ag + in solution, and is illustrated in Fig S2B in the case of 5 equivalents of Ag + : at high charge states (18+, 19+) the dominant forms are SilE:Ag + complexes 1:0 and 1.1, i.e. relatively depleted in silver. On the contrary, at low charge states (9+, 10+, 11+) the dominant complexes are 1:3 and 1:4. This indicates that SilE bound to many silver ions is less prone to bearing many charges. A general trend in MS is that disordered species can unfold and extend to accommodate more charges during electrospray ionization and tend to be observed with larger and higher charge state distributions. Along this line, the HRMS observation suggests that silver-bound SilE is generally less disordered than silver-free SilE, consistent with NMR and CD data. This is also consistent with the IM-MS and SAXS conclusion on the preference for more compact structures for silver binding SilE. Overall is a strong indication that silver binding yields important structural adjustments in SilE. Complementary experimental observations finally allow drawing a global multiscale picture of the changes in SilE structure upon silver binding. At the secondary structure level, NMR and CD data suggest that a-helices are formed as the relative silver concentration increases. Moreover, SAXS, IM-MS, and HRMS data show that this local structuration is accompanied by global compaction. Finally, IM-MS and HRMS data provide evidence that this compaction is correlated with the number of silver cations bound to the complex, rather than in the relative silver concentration.

Silver-bound SilE complexes adopt different conformations

To understand how SilE may accommodate up to eight silver ions, we have built a putative model of SilE based on ColabFold 27 , running the AlphaFold prediction protocol. 28 Indeed, AlphaFold has demonstrated its capability to accurately predict protein structures by using deep learning methods to analyze coevolutionary information. Based on the SilE sequence, we have modeled five structures that adopt a compact or an elongated organization. All the predicted structures display four helical segments that perfectly match the structure of SilE-mimicking peptides (Fig. S8). 10 However, we can rule out the AlphaFold structures that exhibit the weakest pLDDT for the helical regions (Fig. S9A,C). The rank1 structure also suggests that AlphaFold is quite confident about the relative position of a1 with respect to a2 or a3 with respect to a4, due to its low PAE (Fig. S9B). The relative positioning of the different helices along the SilE sequence is also confirmed by our spin relaxation measurements of the silver-bound SilE (Fig. 5). From residues T21 to T52 and Q94 to L105, R1, R2, and hetNOE parameters exhibit similar or close values compared to

Please do not adjust margins

Please do not adjust margins the ones seen in the case of SilE in the free state. This result indicates that these parts remain as flexible as they were in the free form. For the predicted helical regions a3 and a4, our data show a slight decrease in the average R1 from the free to the silver-bound state (1.67 to 1.53 s -1 respectively). This value is concomitant with a strong increase of the average R2 and hetNOE in the silver-bound state compared to the free state (from 5.0 to 10.7 s -1 and from 0.32 to 0.63 respectively). This effect could be explained by a drastic change in the microdynamic parameters with a strong increase of the order parameter S 2 along with a decrease in the local motion tloc (Fig. S10). The modification of these two parameters undoubtedly reflects rigidification and formation of helices a3 and a4. Despite the absence of NMR signals in the a1/a2 regions, it is noticeable that there is a tendency for a significant increase of the R2 and hetNOE parameters, further denoting the presence of a helical formation and a strong conformational exchange. This fact is further supported by the CD spectra recorded for the SilE 57-95 construct that concludes to the formation of an a-helix for this region (Fig. S6). The local increase of the hetNOE values that exactly follow the secondary structure rather than a general increase led also us to conclude that silver-bound SilE has no defined tertiary structure. Therefore, the prediction of different helical regions added to our relaxation data is in good agreement with IM-MS and HRMS that conclude in higher compactness of SilE when bound to silver ions. At this level, it is important to highlight that the compactness is due to the silver-induced helical folding of SilE compared to the disordered state of the free SilE. Whether SilE forms a globular or elongated structure still presents an unresolved conundrum, but a realistic lowresolution tertiary structure may be inferred from SAXS data. Indeed, Alphafold predicts the correct secondary structure of SilE with a perfect match of the four helices but is not suitable to establish the highly flexible tertiary structure of SilE in the presence of Ag + . 29 Therefore, we have used MultiFoXS 30 for computing N-states (N = 1 to 5) models of the [SilE:nAg + ] complex combined with a divide and conquer approach. As a starting point, we have used the compact rank 1 [SilE:nAg + ] complex structure predicted by Alphafold that provides the lowest pLDDT and we have progressively increased the number of possible flexible regions until reaching full flexibility of the SilE intervening linkers (Table S2 for the detailed flexible regions under study). Please do not adjust margins

Please do not adjust margins

This strategy allows the sampling of an extensive conformational space accessible to SilE. In the case of the rank1 globular structure, the only allowed flexible part has been defined to be the N-terminus. As can be seen in Fig. S11A, MultiFoXS could not find a suitable structural ensemble capable of reproducing the SAXS data and shows a c score of 37.4 (Table S2). Additionally, the associated radius of gyration equals 24.6 ± 2.0 Å in complete disagreement with the one experimentally determined from SAXS data (32.9 ± 0.2 Å). This important result indicates that the [SilE:nAg + ] complex does not adopt a globular structure in the liquid state and that the rank1 structure, predicted by AlphaFold, contradicts our experimental data. To keep the positioning of a1,2 and a3,4, we have added the possibility for SilE to remain flexible in a region spanning the linker between a2 and a3. The c score decreases from a onestate model to a three-states model (from 7.3 to 6.3 respectively) while a further increase in the number of conformations to 4 and 5 states did not improve the c score. It is also noteworthy that the c score significantly decreases compared to the one obtained for the globular structure and therefore supports the fact that the [SilE:nAg + ] complex structure adopts a more elongated conformation in the liquid state (Fig. S11B). Thus, we have added another degree of flexibility by allowing the linker between a1,2 or a3,4 region to move freely. While the c score significantly decreases, we did not observe any differences when flexibility affects the linker between a1 and a2 or a3 and a4 with a c score of 1.9 for this ensemble of structures in the case of a 3-states model. In both cases, the different ratios that represent the conformational sampling are rather similar (Fig. S11 C,D). The derived average Rg for the structure that presents a flexible a1,2 intervening linker is 28.9 ± 1.4 Å compared to an average Rg of 31.9 ± 4.0 Å for the structure exhibiting flexibility between the a3,4 linker. Finally, the c score decreases to 0.9 for a 3-states model when the different linkers of the [SilE:nAg + ] complex structure are allowed to move freely (Fig. 6). For the 3-states model, the average Rg of 33.6 ± 3.7 Å (Table S2) is in good agreement with the experimental Rg derived from SAXS data (32.9 ± 0.2 Å). This result allows accurate mapping of the conformational landscape of [SilE:nAg + ] complex and supports the fact that it could be described by different conformations co-existing in solution, from elongated to intermediate states, as concluded by IM-MS data analysis. Overall, the combined use of AlphaFold structure prediction with the NMR-derived structures of silverbound SilE peptides and SAXS data support the fact that (i) the [SilE:nAg + ] complex secondary structure exhibits four helical regions and (ii) different structural organizations of SilE co-exist in solution where Met and His residues may be favorably oriented and accessible to bind silver ions.

Conclusions

Several metals remain essential to life in extant organisms and play a diversity of roles in many different physiological processes. Sometimes, the toxicity of a given element may vary significantly because of speciation, so different chemical species containing the same metal may have a very different impact on living organisms. In the case of silver, Gram-negative bacteria have developed a resistant system based on a specific efflux pump dedicated to the removal of silver ions. Within this system, SilE has been recognized as an essential component for the proper function of the efflux pump by sequestering silver ions. From the vast panorama of the known metal binding proteins, SilE displays specific features in the sense that it is disordered in its free state and is prone to bind several silver ions while folding upon binding. In the present article, we have advanced our knowledge regarding the unknown SilE structure, dynamics, and interaction with silver ions. By combining experimental methods and AlphaFold prediction, we have determined up to eight binding sites that localize in four well-defined helical segments. Beyond these identified binding sites, four of them exhibit strong interaction with silver ions. In particular, we have demonstrated that SilE does not adopt a globular structure upon silver binding but rather samples a large conformational space from elongated to more compact structures. This significant result deviates from the hypothesis established previously 9 and leads us to ask this dangling question: "How SilE may accommodate up to eight silver ions and what is the advantage of such a structural organization?". To answer this crucial question, we must resort to the state of the art and recall that SilE is synthesized only during growth in the presence of silver 4 and is mandatory to avoid the disruption of the efflux pump machinery. At this level, we may hypothesize that SilE acts as a regulator that retains silver ions through His and Met residues and avoids metal saturation in the periplasm. Therefore, to be efficient, this process needs high solvent accessibility of those residues. To lend credence to our hypothesis, we have calculated the solvent accessible surface area (SASA) for (i) the globular single structure and (ii) the ensemble of the three structures found in the case of fully flexible SilE (Fig. 6). As seen on the Table S3, the total percentage of SASA contrasts between the two cases. For the globular case, the average SASA on the interacting His and Met residues is 37 % while it reaches 70 % for the ensemble of elongated structures. Consequently, this structural organization facilitates silver binding through much higher accessibility of the involved His and Met residues while the globular structure strongly prevents close contact with silver ions. The most elongated structure is an amazing way to bind silver and allow the bacteria to remain active even at high silver concentrations. To the best of our knowledge, this is the first example of an intrinsically disordered protein that folds into helices to bind several metal ions at the same time. For a deeper understanding of the cytotoxic mechanism of silver and how its antimicrobial properties are used in medicine, we refer the reader to an extensive review 3 . Nevertheless, some Gram-negative bacteria have developed resistance to silver and impair the proper functioning of its antimicrobial properties. Thus, we may hypothesize that this resistance process is Please do not adjust margins Please do not adjust margins triggered by a high silver concentration in the cell. According to our aforementioned results, SilE can sequestrate up to eight silver ions. The SilCBA pump is also located in the periplasm and oversees the silver efflux out of the cell. Since we have identified two different binding affinities with silver ions (4 strong and 3-4 labile sites), it allows a fine tune of the silver release to the periplasm before extrusion by the efflux pump. This last hypothesis needs to be confirmed by further studies of the different interacting partners in the whole system.

Figure 1 .

 1 Figure 1. Stoichiometry of the Ag-SilE complex. A) Deconvoluted spectra, over the end of the elution band (in order to avoid chemical artefacts), of the 5 µM solutions of SilE and in presence of 1 to 10 molar equivalents of silver ions. The red lines represent the masses of the different stoichiometries [SilE:nAg + ] complexes. B) Relative abundances or branching ratios of the different [SilE:nAg + ] complexes as a function of the molar equivalents Ag + in solution. (C) Electropherograms obtained in CE-ICP-MS for 5 µM SilE and increasing Ag + concentrations in the sample for (1) [SilE:nAg] complex and (2) sorbed Ag+ released by HNO3. Electropherograms are shifted in 107 Ag signal intensity for sake of clarity, background signal is around 50 counts for all. (D) Peak areas measured from C) for Ag species as a function of the Ag + equivalents added in the sample for [SilE:nAg] complex (black) and sorbed Ag + released by HNO3 (blue). (E) Electropherograms obtained in CE-ICP-MS for 5 µM SilE and increasing Ag + concentrations in the electrolyte determined using CE-ICP/MS analyses. Electropherograms are shifted in 107 Ag signal intensity for sake of clarity. (F) Peak areas measured from E) for [SilE:nAg] complex as a function of the Ag + concentration in the electrolyte.

Figure 2 .

 2 Figure 2. Backbone NMR assignments and dynamics of SilE. 2D 1 H-15 N HSQC NMR spectra of SilE at 140 µM A) in the free state and (B) in the presence of 6 equivalents of silver ions. C) Secondary structure propensity calculated from NMR backbone assignments of free SilE. Positive values indicate a-helical structure propensity. D) NMR relaxation data 15 N-R1, and -R2 and hetNOEs for the free SilE.

Figure 3 .

 3 Figure 3. Silver binding sites. CSP observed at the end of the titration for 6 Ag + equivalent on a 140 µM 15 N-SilE (20 mM MES buffer pH 6.8 and 20 mM NaF). Residues 59-92 (except 76) are missing due to the large signal broadening in 1 H-15 N HSQC spectra.

Figure 4 .

 4 Figure 4. Structuration of SilE upon silver addition. (A) CD curves recorded at 25°C for the free SilE and in presence of 1 to 8 equivalents of silver ions ([SilE] = 20 µM in 20 mM MES buffer pH 6.8 implemented with 20 mM NaF). The decrease of the ellipticity at 222 nm upon silver ions addition demonstrate the progressive formation of a-helices. (B) Secondary structure propensity calculated from NMR backbone assignments of SilE + 6 equivalents of Ag + . (C) Extracted arrival time distributions for selected [SilE:nAg] 6+ complexes with different stoichiometries emphasizing the relative increase of the compact (low arrival times) species as n increases. (D) Evolution of the relative intensity of the compact form as a function of the number n of Ag bound to SilE 6+ . Error bars correspond to 95% confidence intervals from repeated measurements with different relative Ag + concentrations in the sample. The full line represents a fit using an exponential decay function.

Figure 5 .

 5 Figure 5. NMR spin relaxation data of the [SilE:nAg] complex. On top, the [SilE:nAg] complex secondary structure deduced from AlphaFold is sketched. NMR spin relaxation data 15 N-R1, and -R2 and hetNOE for SilE in presence of 6 equivalents of Ag + in solution are displayed in red. For the sake of comparison, the NMR relaxation data of the free SilE have been added in black.

Figure 6 .

 6 Figure 6. SAXS analysis of [SilE:nAg] complex. On top, the [SilE:nAg] complex secondary structure deduced from Alphafold is sketched. The different flexible segments used in the MultiFoXS calculation are coloured in blue. SAXS data (blue circle), back-calculated SAXS data (red line), and residuals are represented in the left panel. On the right panel, the representation of the top-scoring N-ensemble of structures and their respective weight that best fit the SAXS data is presented close to each structure and are determined using the MultiFoXS server. Residuals per ith SAXS data point have been calculated according to the following expression: (I(qi) exp -I(qi) calc )/si. I(qi) exp is the measured scattering intensity, I(qi) calc is the intensity calculated using MultiFoXS for scattering vector qi and si is the experimental error.

This journal is © The Royal Society of Chemistry 20xx J. Name., 2013, 00, 1-3 | 11

Acknowledgements

This work was supported by the CNRS through the MITI interdisciplinary program Metallomix and the French ANR (grant # ANR-17-CE29-0013-01). Financial support from the IR INFRANALYTICS FR2054 for conducting the research is gratefully acknowledged. The 900 MHz and the high-pressure NMR experiments have been performed at IMEC-IS-UCCS Lille and ICSN Gif sur Yvette respectively.

Author Contributions

MH coordinated, managed the overall project, and wrote the manuscript. MH analyzed NMR spin relaxation. YM engineered the SilE plasmid, optimized the production and purification protocols, performed NMR assignments and titration analysis. CA and MM prepared the SilE samples and performed CD analysis. AH and CD performed and analyzed the CE-ICP-MS experiments. FC and CCZ performed and analyzed IM-MS experiments. MG performed the HRMS experiments and in collaboration with LMA analyzed the MS data. OW performed MultiFoXs analysis by combining AlphaFold simulation and SAXS data interpretation. YM, FC, LMA, MG, AH, CCZ and OW assisted in writing the manuscript, reviewed the final version, and approved the content and submission.

Conflicts of interest

There are no conflicts to declare.