Armelle Vigouroux

Thibault Meyer

Anaïs Naretto

Pierre Legrand

Magali Aumont-Nicaise

Aurélie Di Cicco

Sébastien Renoud

Jeanne Doré

Daniel Lévy

Ludovic Vial

Céline Lavire

Solange Moréra
email: solange.morera@i2bc.paris-saclay.fr

Characterization of the first tetrameric transcription factor of the GntR superfamily with allosteric regulation from the bacterial pathogen Agrobacterium fabrum

A species-specific region, denoted SpG8-1b allowing hydroxycinnamic acids (HCAs) degradation is important for the transition between the two lifestyles (rhizospheric versus pathogenic) of the plant pathogen Agrobacterium fabrum. Indeed, HCAs can be either use as trophic resources and/or as induced-virulence molecules. The SpG8-1b region is regulated by two transcriptional regulators, namely, HcaR (Atu1422) and Atu1419. In contrast to HcaR, Atu1419 remains so far uncharacterized. The highresolution crystal structures of two fortuitous citrate complexes, two DNA complexes and the apoform revealed that the tetrameric Atu1419 transcriptional regulator belongs to the VanR group of Pfam PF07729 subfamily of the large GntR superfamily. Until now, GntR regulators were dimeric. Here, we showed that Atu1419 represses three genes of the HCAs catabolic pathway. We characterized both the effector and DNA binding sites and identified key nucleotides in the target palindrome. From promoter activity measurement using defective gene mutants, structural analysis and gel-shift assays, we propose N5,N10-methylenetetrahydrofolate as the effector molecule, which is not a direct product/substrate of the HCA degradation pathway. The Zn 2+ ion present in the effector domain has both a structural and regulatory role. Overall, our work shed light on the allosteric mechanism of transcription employed by this GntR repressor.

INTRODUCTION

Agrobacterium fabrum has two lifestyles: it can interact with a large variety of plants as a rhizosphere inhabitant or as a pathogen when it harbors a tumor-inducing plasmid (which is a virulence plasmid), and transfers a portion of this to the plant cells upon infection, resulting in the crown-gall disease [START_REF] Nester | Agrobacterium: nature's genetic engineer[END_REF][START_REF] Dessaux | Niche Construction and Exploitation by Agrobacterium: How to Survive and Face Competition in Soil and Plant Habitats[END_REF][START_REF] Meyer | Ecological Conditions and Molecular Determinants Involved in Agrobacterium Lifestyle in Tumors[END_REF].

A. fabrum possesses a species-specific region, denoted SpG8-1b (a region present in strains of this species but absent from other Agrobacterium species), located in the circular chromosome and responsible for hydroxycinnamic acids (HCAs) degradation such as ferulic acid, caffeic acid and pcoumaric acid (4,5) (Figure 1). These latter compounds are common plant secondary metabolites being precursors of lignin incorporated into plant cell walls. They are abundantly released in soil during the decay of root cells and are significant environmental molecules for soil-and plant-interacting bacteria [START_REF] Bhattacharya | The roles of plant phenolics in defence and communication during Agrobacterium and Rhizobium infection[END_REF]. Although HCAs are generally a strong bacterial repellent, they appear to be chemoattractants in the case of rhizobia and agrobacteria for which they can be used as trophic resources and/or inducedvirulence molecules [START_REF] Guo | Is there any crosstalk between the chemotaxis and virulence induction signaling in Agrobacterium tumefaciens?[END_REF][START_REF] Parke | Chemotaxis to plant phenolic inducers of virulence genes is constitutively expressed in the absence of the Ti plasmid in Agrobacterium tumefaciens[END_REF][START_REF] Kape | Chemotaxis and nod Gene Activity of Bradyrhizobium japonicum in Response to Hydroxycinnamic Acids and Isoflavonoids[END_REF][START_REF] Kalogeraki | The phenolic vir gene inducer ferulic acid is O-demethylated by the VirH2 protein of an Agrobacterium tumefaciens Ti plasmid[END_REF][START_REF] Meyer | Regulation of hydroxycinnamic acid degradation drives Agrobacterium fabrum lifestyles[END_REF]. We have previously showed that HCAs degradation via the SpG8-1b region interferes with virulence genes expression suggesting that this metabolic pathway is important for the transition between the two lifestyles (rhizospheric versus pathogenic) of Agrobacterium [START_REF] Meyer | Regulation of hydroxycinnamic acid degradation drives Agrobacterium fabrum lifestyles[END_REF]. Such a transition requires a fine-tune regulation of gene expression to express the appropriate genes at the right time [START_REF] Duprey | Bacterial virulence and Fis: adapting regulatory networks to the host environment[END_REF][START_REF] Valentini | Lifestyle transitions and adaptive pathogenesis of Pseudomonas aeruginosa[END_REF][START_REF] Barton | Ecological and evolutionary dynamics of a model facultative pathogen: Agrobacterium and crown gall disease of plants[END_REF].

The SpG8-1b genomic region encodes six enzymes (Figure 1A). Five of them operate in this sequential order for the degradation of ferulic acid leading to different intermediates [START_REF] Campillo | Analysis of hydroxycinnamic acid degradation in Agrobacterium fabrum reveals a coenzyme A-dependent, beta-oxidative deacetylation pathway[END_REF]: Atu1416 (a feruloyl-CoA synthase), Atu1417 (an enoyl-CoA hydratase), Atu1415 (a phenylhydroxypropionyl-CoA dehydrogenase), Atu1421 (a 4-hydroxy-3-methoxyphenyl-β-ketopropionyl-CoA β-keto-thiolase) and Atu1420 (a tetrahydrofolate-dependent vanillate O-demethylase). Atu1420 shares 56% sequence identity with Sphingomonas paucimobilis and Sphingobium sp. SYK-6 homologues of known structures (PDB 5TL4 [START_REF] Kohler | Structure of aryl O-demethylase offers molecular insight into a catalytic tyrosine-dependent mechanism[END_REF] and PDB 5X1I [START_REF] Harada | The crystal structure of a new O-demethylase from Sphingobium sp. strain SYK-6[END_REF], respectively). This enzyme degrades vanillic acid into protocatechuic acid using the tetrahydrofolate (H4F) cofactor and producing the N5-methyltetrahydrofolate compound (MH4F) [START_REF] Campillo | Analysis of hydroxycinnamic acid degradation in Agrobacterium fabrum reveals a coenzyme A-dependent, beta-oxidative deacetylation pathway[END_REF]. Protocatechuic acid is then processed by pca genes to produce intermediates of the tricarboxylic acid (TCA) cycle [START_REF] Parke | Supraoperonic clustering of pca genes for catabolism of the phenolic compound protocatechuate in Agrobacterium tumefaciens[END_REF]. Based on sequence identity of 29% with the N5,N10-methylene tetrahydrofolate reductase MetF of Sphingobium sp. SYK-6, Atu1418 the sixth enzyme of the SpG8-1b genomic region, is likely a 5,10-methylene tetrahydrofolate reductase which transforms MH4F into N5,N10-methylenetetrahydrofolate (MEF), which would, in turn, allow the regeneration of H4F (4,18) (Figure 1A).

Although the SpG8-1b genomic region degrades both ferulic and p-coumaric acids, this latter induces the gene expression of the first part of the SpG8-1b region while ferulic acid induces that of the whole pathway [START_REF] Meyer | Regulation of hydroxycinnamic acid degradation drives Agrobacterium fabrum lifestyles[END_REF]. Therefore, two expression units are differently regulated suggesting the existence of two regulators. These regulators are Atu1419 (predicted to be a GntR member) and Atu1422 denoted HcaR (for hydroxycinnamic acid catabolic repressor belonging to the MarR family), respectively [START_REF] Lassalle | Genomic species are ecological species as revealed by comparative genomics in Agrobacterium tumefaciens[END_REF][START_REF] Meyer | Regulation of hydroxycinnamic acid degradation drives Agrobacterium fabrum lifestyles[END_REF] (Figure 1B). HcaR is the repressor of its own transcription and that of the first atu1416 and atu1417 genes of HCA degradation pathway [START_REF] Meyer | Regulation of hydroxycinnamic acid degradation drives Agrobacterium fabrum lifestyles[END_REF]. However, because atu1415 and atu1416 belong to the same transcription unit as well as atu1421 and hcaR do, HcaR also regulates atu1415 and atu1421 genes expression, totalizing the regulation of five genes (Figure 1B). In contrast to HcaR, the second regulatory protein Atu1419 likely involved in the second part of the pathway corresponding to the vanillate degradation remains uncharacterized so far.

Herein, we investigated the molecular role and structural aspects of Atu1419 combining in vitro and in vivo approaches. We first proved that Atu1419 was a transcription repressor of three genes of the second part of the HCA degradation pathway. None of HCA degradation intermediates could release Atu1419 from DNA binding. Nonetheless, structural analysis of five high-resolution crystal structures of Atu1419 in apoform, in complex with a fortuitous citrate molecule bound to the effector domain originated from the crystallization condition (2 structures), in complex with DNA and in complex with both DNA and citrate, helped us infer a possible effector molecule, which was confirmed by gel shift assays and microcalorimetry. The structures revealed Atu1419 to be a member of the VanR group of the FadR C-terminal domain (FCD; Pfam PF07729) subfamily of the large GntR superfamily of transcriptional factors (>93135 members in Pfam database) [START_REF] Hoskisson | Chapter 1: Variation in form and function the helix-turnhelix regulators of the GntR superfamily[END_REF][START_REF] Rigali | Subdivision of the helix-turn-helix GntR family of bacterial regulators in the FadR, HutC, MocR, and YtrA subfamilies[END_REF][START_REF] Suvorova | GntR Family of Bacterial Transcription Factors and Their DNA Binding Motifs: Structure, Positioning and Co-Evolution[END_REF][START_REF] Jain | Allosteric control of transcription in GntR family of transcription regulators: A structural overview[END_REF]. The FCD subfamily encompasses two groups of regulators namely FadR and VanR. So far, the dimeric FadR from Escherichia coli was the best characterized of the FadR group of FCD-GntR regulators shown to be regulated by the acyl-CoA effector and able to bind specific palindromic DNA through a winged Helix-Turn-Helix (wHTH) motif [START_REF] Van Aalten | The structural basis of acyl coenzyme A-dependent regulation of the transcription factor FadR[END_REF]. Like FadR, Atu1419 possesses a characteristic molecular architecture, composed of a conserved N-terminal DNA binding domain containing the wHTH motif and a C-terminal effector binding/oligomerization domain. Unlike FadR, this latter domain of Atu1419 displays six helices instead of seven observed in members of FadR group. Atu1419 is the first example of a transcriptional regulator of the whole GntR superfamily to be tetrameric. Our work brings new insights into mechanistic aspects of such repressor, which uses an induced-allosteric mechanism for DNA release upon effector binding.

MATERIAL AND METHODS

Bacterial strains and growth conditions

The bacteria and plasmids used for this study are listed in Supplementary Table S1. Escherichia coli were grown routinely, with shaking (150 rpm), at 37°C in LB medium. Growth media were supplemented with appropriate antibiotics (tetracycline, 10 µg/ml; gentamicin, 15 µg/ml; ampicillin, 100 µg/ml) when necessary. The A. fabrum strains were grown with shaking (160 rpm), at 28°C in YPG (Yeast Peptone Glucose)-rich medium or in AT minimal medium supplemented with 10 mM succinate and 10 mM ammonium sulfate. AT minimal medium was supplemented with 750 μM of ferulic acid or citrate and with the appropriate antibiotic (gentamicin, 10 µg/ml). Ferulic acid and MH4F were obtained from Sigma Aldrich (St. Louis, USA) and MEF from Merck Company (Switzerland).

Construction of the deletion mutant C58Δatu1419 and transcriptional fusions

The C58Δatu1419 and C58Δatu1420 strains were constructed according to a strategy as described [START_REF] Lassalle | Genomic species are ecological species as revealed by comparative genomics in Agrobacterium tumefaciens[END_REF][START_REF] Meyer | Regulation of hydroxycinnamic acid degradation drives Agrobacterium fabrum lifestyles[END_REF]. Vectors containing the recombinant region (amplified by PCR with specific primers listed in Supplementary Table S2), flanking downstream and upstream of the atu1419 or atu1420 genes, were introduced into A. fabrum C58 by electroporation. Single-crossover integration was selected by gentamycin resistance on YPG medium plates. Gentamycin-resistant colonies were spread on YPG plates containing 5% of sucrose to obtain plasmid excision and double-crossover events leading to nonpolar mutants. atu1419 and atu1420 deletions were confirmed by PCR analysis and DNA sequencing (GenoScreen, Lille, France). pOT1e [START_REF] Allaway | Use of differential fluorescence induction and optical trapping to isolate environmentally induced genes[END_REF] transcriptional fusions of the promoter regions of SpG8-1b genes, namely Patu1419 and Patu1420 were generated as described [START_REF] Meyer | Regulation of hydroxycinnamic acid degradation drives Agrobacterium fabrum lifestyles[END_REF] (Supplementary Table S2 for specific primers). Reporter constructions were introduced into A. fabrum C58 wild-type and derivatives by electroporation, and gentamycin-resistant colonies were selected.

Measurement of promoter activity

Genes expression was measured after 24 hours in A. fabrum with a pOT1e plasmid harboring an eGFP transcriptional fusion [START_REF] Allaway | Use of differential fluorescence induction and optical trapping to isolate environmentally induced genes[END_REF] as described [START_REF] Meyer | Regulation of hydroxycinnamic acid degradation drives Agrobacterium fabrum lifestyles[END_REF]. Results were normalized by dividing the fluorescence level by the optical density at 600 nm values. At least five technical replicates and two biological replicates were performed for each condition. Differences between conditions were determined with Tukey test (P-value= 0.05).

Oligonucleotides and DNA preparation

The synthetic palindromic oligodeoxyribonucleotide 5'-ATGTATACAT-3' was purchased from Sigma-Aldrich (Darmstadt, Germany). Oligonucleotide solution in sterile water at 2 mM was hybridized by heating to 90°C for 5 min and cooling in the crystallization room at 18°C overnight to produce a 10-mer DNA.

The atu1420 promoter regions were amplified with specific primer pairs listed in Supplementary Table S2. The atu1416-1417, atu1418-1419, virB and hcaR promoter regions were prepared as previously described [START_REF] Meyer | Regulation of hydroxycinnamic acid degradation drives Agrobacterium fabrum lifestyles[END_REF] and purified with a PCR Clean-up kit (Macherey-Nagel, Düren, Germany).

Cloning, expression and purification of Atu1419 and Atu1419-H3A mutant

Coding sequences for Atu1419 was amplified by PCR adding a C-terminal 6-Histidine tag and using atu1419F and atu1419R primers and was inserted into the NdeI/XhoI restriction sites of the pET-20b vector (Novagen, Merck Biosciences, France). The nucleotide sequence was confirmed by DNAsequence analysis (GATC Biotech, Mulhouse, France). E. coli BL21 competent cells transformed with pET-20b-Atu1419 were grown in LB media until OD600 of 0.8 and protein production was induced by 0.5 mM isopropyl β-D-thio-galactopyranoside (IPTG) for 3h at 37°C. Cells were centrifuged at 4000 g for 15 min at 4°C, resuspended in a buffer of 50 mM Tris-HCl pH 8, 300 mM NaCl and 20 mM imidazole and disrupted by sonication. After centrifugation at 25000 g for 30 min at 4°C, the filtrated supernatant was loaded onto a 5 ml His-Trap HP column (GE Healthcare, Chicago, Illinois, USA). After a washing step of 6% with 50 mM Tris-HCl pH 8, 300 mM NaCl and 300 mM imidazole (Buffer B), protein elution was performed with Buffer B. Protein fractions were loaded onto a gel filtration column (HiLoad 26/60 Superdex 200 prep grade, GE Healthcare) equilibrated with 50 mM Tris-HCl pH 8 and 150 mM NaCl.

The protein fractions were pooled, concentrated and stored at -80°C.

The synthetic gene (Genscript, Piscataway, New Jersey, USA) coding for Atu1419 mutant (H192A, H141A and H214A) namely Atu1419-H3A was inserted into pET-20b. The Atu1419-H3A mutant was expressed and purified as the wild-type protein described above.

Size exclusion chromatography (SEC) and SEC-MALS. For size exclusion chromatography (SEC)

and SEC coupled to multi-angle light scattering (SEC-MALS) analyses, apo Atu1419 was prepared at 25 M and Atu1419 in complex with the palindromic DNA was prepared with a ratio protein:DNA of 1:2 in a buffer containing 50 mM Tris pH 8 and 150 mM NaCl. The same buffer was used as the mobile phase for SEC using a Superdex 200 10/300 GL column on an AKTA FPLC system (GE Healthcare) and on a Shimadzu HPLC. Multiangle light scattering was detected with a MiniDAWN TREOS light scattering module and a refractometer Optilab T-rEX (Wyatt Technology).

Electrophoretic mobility shift assay (EMSA)

Atu1419 or Atu1419-H3A mutant were mixed with different promoter regions (Patu1416-1417, Patu1418-1419, PhcaR, Patu1420 and PvirB) and with variants of the atu1418-1419 and atu1420 regions. The intergenic regions of virB was used as a nonspecific control probe. Gel-shift assays were performed in 10 µl reaction mixture containing 30 nM of DNA probe without and with Atu1419 at different concentrations in 50 mM Tris-HCl pH 8 and 150 mM NaCl. 50 to 300 µM of MEF, MH4F, H4F or citrate were added for testing their influence on Atu1419-DNA complex formation. After incubation for 30 min at room temperature, the samples were separated by electrophoresis in TBE buffer (45 mM Tris-HCl pH 8, 45 mM boric acid and 1 mM EDTA) on non-denaturing 6% or 12% polyacrylamide gels at 150 V and 4°C for 2 h. Gels were then stained with either SYBR® Green EMSA nucleic acid gel stain (Invitrogen, Carlsbad, CA, USA) or ethidium bromide for 20 min. DNA was visualized under UV light (Fisher Bioblock Scientific, Illkirch, France or UVP BioDoc-it2 Imager, Analytic Jena, Germany).

Crystallization and structure determination of Atu1419

Crystallization conditions are summarized in Table 1. For all protein samples (140 µM tetrameric protein alone or in complex with 700 µM palindromic DNA), conditions were screened using QIAGEN kits (Valencia, CA, USA) with a Mosquito nanodrop robot (TTP Labtech, Melbourn, Great Britain) and were manually optimized at 20°C in hanging drop by mixing equal volumes of the protein or protein-DNA solution with precipitant solution. Crystals were transferred to a cryoprotectant solution (paraffin oil or mother liquor supplemented with 20% PEG 400) and flash-frozen in liquid nitrogen. Diffraction data were collected at 100 K on PROXIMA 1 and PROXIMA 2 beamlines at synchrotron SOLEIL (Saint-Aubin, France). Intensities were integrated using the XDS package (25) (Table 1).

The first structure of Atu1419 in P212121 space group was determined at 2.35 Å resolution by singlewavelength anomalous dispersion (SAD) method at the peak absorption energy of Zinc (Table 1). The presence of a zinc metal in crystals was found thanks to an X-ray fluorescence emission scan on the beamline. The steps of zinc ion substructure determination, phases calculation and density modification were performed using CRANK from CCP4 (Collaborative Computational Project, Number 4) and a partial model was built using BUCCANEER (CCP4). This model was then used to calculate the phases of a higher resolution dataset at 2 Å resolution leading to the complete polypeptide chain model. A tetramer is present in the asymmetric unit. Other structure determinations were performed by molecular replacement with PHASER (26) using the first refined structure of Atu1419 (monomer, dimer or tetramer). Because of the anisotropy of the diffraction of Atu1419-DNA complex in P6422 crystals and apo Atu1419, the DEBYE and STARANISO programs developed by Global phasing Ltd were applied to the data scaled with XDS using the STARANISO server (http://staraniso.globalphasing.or/). These programs perform an anisotropic cut-off of merge intensity data on the basis of an analysis of local I/(I), compute Bayesian estimates of structures amplitudes, taking into account their anisotropic fall-off, and apply an anisotropic correction to the data. The corrected anisotropic amplitudes were used for further refinement. Refinement of each structure was performed with BUSTER-2.10 (27) employing TLS groups and NCS restraints. Inspection of the density maps and manual rebuilding were performed using COOT [START_REF] Emsley | Coot: model-building tools for molecular graphics[END_REF]. Refinement details of each structure are shown in Table 1. Molecular graphics images were generated using PyMOL (http://www.pymol.org).

Circular dichroism experiments (CD)

Circular dichroïsm in the far-UV region was performed using a spectropolarimeter (Jasco J-810, Jasco, Lisses, France) equipped with a water-cooled Peltier unit (Jasco circular dichroïsm spectrometer model J810). Apo Atu1419 and Atu1419-H3A mutant were both concentrated at 25 M in 50 mM Tris pH 8 and 150 mM NaCl. Spectra were recorded in a cell width of 0.1 mm path length (121.QS, Hellma, Hellma Analytics, Müllheim, Baden Württemberg, Germany) from 190 to 240 nm at 20°C. Five consecutive scans from each sample were merged to produce an averaged spectrum; the spectra were corrected using buffer baselines measured under the same conditions. Data were recorded in mdeg and converted using the mean residues ellipticity method (deg.cm 2 mol -1 ). Secondary structure estimates were derived from the normalized spectra using the CDSSTR, SELCON3, CONTIN of the DICHROWEB server [START_REF] Sreerama | Estimation of protein secondary structure from circular dichroism spectra: comparison of CONTIN, SELCON, and CDSSTR methods with an expanded reference set[END_REF][START_REF] Whitmore | DICHROWEB, an online server for protein secondary structure analyses from circular dichroism spectroscopic data[END_REF].

Differential scanning calorimetry (auto PEAQ DSC)

Thermal stability of 10 μM apo Atu1419, 10 μM Atu1419 in the presence of 100 μM MEF and 20 μM Atu1419-H3A mutant was performed by DSC on an auto PEAQ DSC (Malvern, France) in a standard buffer. Each measurement was preceded by a baseline scan with the standard buffer. Scans were performed at 1 K.min -1 between 20 and 90°C. The heat capacity of the buffer was subtracted from that of the protein sample before analysis. Thermodynamic parameters were determined by fitting the data to the following equation:

∆𝑪 𝒑 (𝑻) = 𝑲 𝒅 (𝑻) ∆𝑯 𝒄𝒂𝒍 ∆𝑯 𝒗𝑯 [𝟏 + 𝑲 𝒅 (𝑻)] 𝟐 𝑹𝑻 𝟐
where Kd is the equilibrium constant for a two-state process, ΔHvh is the enthalpy calculated on the basis of a two-state process and ΔHcal is the measured enthalpy.

Isothermal titration microcalorimetry measurements (ITC)

Isothermal titration microcalorimetry experiments were performed with an ITC200 isothermal titration calorimeter from MicroCal (Malvern, Orsay, France). The experiments were carried out at 20°C. Protein concentration in the microcalorimeter cell (0.2 ml) was 100 µM. Nineteen injections of 2 µl MEF solution at 1.2 mM were performed at intervals of 180 s while stirring at 500 rpm. The experimental data were fitted to theoretical titration curves with software supplied by MicroCal (ORIGIN®).

Electron Microscopy

Atu1419 (150 nM) was incubated with equimolar Patu1418-1419 region (370 bp containing two palindromes separated by 190 bp) for 10 min at room temperature in a buffer containing 50 mM Tris pH 8 and 150 mM NaCl. Samples were deposited on a glow-discharged carbon coated grid and stained with 2% uranyl acetate. Images were recorded with a Lab6 Tecnai Spirit operating at 80kVolt and a Quemesa Olympus CCD camera. Nominal magnification was X 45 000 corresponding to 3.17 Å/pixel.

RESULTS AND DISCUSSION

The transcriptional regulator Atu1419 is the second repressor of the HCAs degradation pathway

To study the role of Atu1419 in transcriptional regulation of the HCAs degradation genes, we constructed reported fusions with pOT1e plasmids containing each promoter region of the SpG8-1b region (Figure 1) cloned upstream the egfp gene. This allowed reporting the transcription of each gene by measuring the eGFP fluorescence. The atu1419 and atu1420 plasmid reporter gene fusions were constructed in this study. Those corresponding to atu1416 and atu1418 were already available [START_REF] Meyer | Regulation of hydroxycinnamic acid degradation drives Agrobacterium fabrum lifestyles[END_REF].

The plasmid reporter fusions were introduced into the wild-type A. fabrum C58 strain to evaluate gene expression in response to ferulic acid. Compared to the fluorescence level measured in the absence of ferulic acid, the fluorescence level measured in the presence of ferulic acid was higher for the four reporter fusions (Table 2). Hence, atu1416, atu1418, atu1419 and atu1420 genes were induced in the presence of ferulic acid. The four plasmid reporter fusions were also introduced into the C58Δatu1419 mutant strain to compare genes expression with that of the wild-type C58 strain in the absence of ferulic acid. The fluorescence level measured after 24 hours for atu1418, atu1419 and atu1420 reporter fusions were respectively 3.67, 3.92 and 1.58 times higher in the C58Δatu1419 mutant strain than in the wild-type strain (Table 2). Thus, in the strain lacking Atu1419 regulatory protein (C58Δatu1419), the reporter fusions for atu1418, atu1419 and atu1420 genes were constitutively expressed. In contrast, the fold change between C58Δatu1419 and the wild type strains was 0.63 and 0.97 for the atu1416 gene fusion and the empty pOT1e vector, respectively indicating that atu1416 gene expression was not repressed by Atu1419 (Table 2). These results showed that Atu1419 was a transcriptional repressor that regulates atu1418, atu1420, and its own transcription.

DNA-binding targets

To characterize the target DNA sequence of Atu1419, we purified Atu1419 (theoretical molecular weight of ~26 542 Da for a monomer), which appeared as a tetramer in solution according to gel filtration/SEC-MALS chromatography (Supplementary Figure S1). Using the intergenic regions (between 300 and 400-base pair (bp)) of atu1416-1417, atu1418-1419, atu1420, hcaR and virB (negative control gene outside the SpG8-1b region) as DNA probes for binding assays (Figure 2A), we showed that Atu1419 bound to the atu1418-1419 and atu1420 promoter regions only, in agreement with the gene expression data described above.

In silico comparison of the promoter regions of atu1418-1419 and atu1420 using BPROM program [START_REF] Solovyev | Automatic Annotation of Microbial Genomes and Metagenomic Sequences[END_REF] revealed three 10-mer identical palindromic sequences (5'-ATGTATACAT-3'), two in atu1418-1419

promoter region and one in atu1420 promoter region (Figure 2B). A palindromic sequence overlaps the -10 box regulatory element of atu1418 and atu1420, whereas in the atu1418-1419 promoter sequence, an additional palindromic sequence is found between the -10 and the -35 regulatory elements of the atu1419 (Figure 2B). Using electrophoretic mobility shift assay (EMSA), we showed that Atu1419 was able to bind each palindrome of the atu1418-1419 promoter sequence without the presence of the other palindrome meaning that only one palindrome site was required for binding (Figure 2C). Nonetheless, the intensity of the shifted bands was greater for the DNA fragment containing both palindromes (Figure 2C). We also analyzed the quaternary structure of Atu1419 in complex with the 10-bp palindrome using gel filtration/SEC-MALS measurements (Supplementary Figure S1). Atu1419 was also tetrameric upon DNA binding.

Crystal structures of Atu1419 in complex with a fortuitous ligand

We solved the first structure of Atu1419 from crystals grown in sodium citrate buffer using SAD method at the peak absorption energy of Zinc at 2 Å resolution in the space group P212121 (Table 1). The zinc ion comes from E. coli protein expression because no metal was added during protein purification and crystallization. Four molecules (A, B, C and D) are present in the asymmetric unit and form a tetramer (a dimer of dimer) (Figure 3A), which is consistent with the observation that Atu1419 is a tetramer in solution (Supplementary Figure S1). Each monomer consists of an N-terminal DNA binding domain (residues 1-71) and a C-terminal all -helical effector binding domain (residues 76-244) (Figure 3B).

The DNA binding domain is composed of three helices (-3) with helix 2 and helix 3 forming the helix-turn-helix motif and two anti-parallel -strands (1-2) connected by a small loop designated as the wing motif. The secondary element with six helices (-9) and topology of the C-terminal domain place Atu1419 in the VanR group of the FCD subfamily of GntR transcriptional regulators.

Approximately 2776 Å 2 of accessible surface area is buried upon dimerization for dimers AB or CD corresponding to an average of 11.7% of the total surface area of each monomer and 41 amino acids per monomer. The interface between dimers AB and CD involves side chains of both domains: helix 3 and strand 1 (residues 49-67) of the DNA binding domain and the three helices 4 (residues 78-100), 7 (residues 150-163) and 9 (residues 208 and 212) of the effector binding domain (Figure 3B). Arg72

in the linker region participates to the dimer interface. Two salt bridges with Glu57-Arg63' (the prime refers to the second monomer in the dimer), and two polar interactions Glu57-Ser53' and Asn62-Arg155' are located in the N-terminal domain whereas fourteen H-bonds/salt-bridges including Asp81, Glu84, Arg95, Arg100, Ser150, Arg155, Glu156, Glu208 and Arg212 belong to the C-terminal domain. The contacts of dimer AD through eleven polar interactions/salt bridges bury 686 Å 2 per subunit and concern 23 amino acids, which are located in the loop between helices 5-6 (mainly residues 123-125), helices 8 (residues 183-194) and 9 (residues 203-216) (Figure 3B). The dimer BC displays a smaller interface of 557 Å 2 per subunit comprising 21 amino acids and the same structural elements as for dimer AD. The total surface contact area between two neighbouring subunits within the tetramer is characteristic of biological interactions [START_REF] Janin | Macromolecular recognition in the Protein Data Bank[END_REF].

Remarkably, the four subunits within the tetramer are not identical: molecule A shows large conformational changes as indicated by an average root mean square deviation (RMSD) of 1.28 Å for all C atoms compared with the subunits B/C/D, which are more similar (average RMSD of 0.7 Å for all C atoms). This is due to the presence of a citrate originated from the crystallization condition, which binds molecule A in the effector binding site in a different orientation (orientation A; Figures 3C and3D)

than that for citrates in subunits B/C/D (orientation B; Figures 3C and3E), which bind similarly (Figure 3C). Four loop regions (shown in red in Figure 3C), comprised of residues 64-67 of the wing between S2).

1

Crystal structures of Atu1419 in complex with the palindromic DNA

The structures of two palindromic DNA-Atu1419 complexes solved at 2.79 Å resolution (P6422 space group) and 2.05 Å resolution (C2221 space group) present a distinct neighboring crystal packing with an asymmetric unit containing a monomer bound to a single DNA strand and a dimer bound to the 10mer palindromic DNA, respectively (Table 1, Supplementary Figure S3). Nonetheless, a tetramer, where dimers AB and CD bind each to a DNA palindrome, was reconstituted by crystal symmetries (Figures 4A and4B) in line with the conservation of Atu1419 tetramer in solution upon DNA binding.

The three monomers of the asymmetric units (one in P6422 and two in C2221) are similar with an average RMSD value of 1 Å for all Cα atoms making almost identical interactions with DNA, although the positions of the DNA binding domains within the dimers between the two crystal structures do not completely overlap (Figure 4C). This observation explains the different crystal packing resulting of flexibility from both the DNA ligand and the DNA binding domains of Atu1419. Each monomer recognizes a half-site DNA with helices 2 and 3 of the HTH motif making nine polar interactions with both strands in the major groove of the DNA (Figures 4D and4E). Helix 3 via Ser44, Thr46, Arg45 and Arg49 side chains is responsible for six hydrogen bonds with both strands, half with oxygens of phosphate groups 5' of guanine (G3) at position 3 on one strand and of adenine (A5) at position 5 on the complementary strand. The other half is made by the guanidinium group of Arg45, which provides specific contacts with the N7 and O6 atoms of G3 on one strand and by OG1 of Thr46 with the N7 of A5 on the other strand (Figure 4C). The remaining H-bonds come from helix 2 via the Arg31 side chain and the main chain amino group of Glu34, and an additional interaction is present with the main chain NH of His9 in helix 1. They consist of phosphate contacts with A1 and T2 on one DNA strand and T4 on the other.

All Atu1419-DNA subunits contain a bound Zn 2+ ion at the same position as that observed in both Atu1419-citrate complexes. In contrast to the protein-DNA complex in C2221 space group, which crystallized in MES buffer, the P6422 protein-DNA complex reveals a bound citrate in the effector binding pocket. This citrate adopts the orientation B described in both Atu1419-citrate complexes, except for molecule A (orientation A) in the P212121 structure. Because both tetrameric Atu1419-DNA structures are similar, the presence of the bound citrate in orientation B in all subunits of one DNA complex has clearly no effect on Atu1419 for DNA binding and DNA release, in agreement with EMSA assay (Figure 4G) and promoter activity (Table S3).

We modified the specifically recognized guanine base G3 by an adenine base and its partner base C8 by a thymine on each strand within the unique palindrome of the atu1420 promoter region. Gel-shift assays showed that Atu1419 was no longer able to bind Patu1420s AT mutant (Figure 4F). Therefore, the guanine at position 3 has an essential role in the DNA palindrome binding by Atu1419 meaning that Arg45 has a key role in recognizing this guanine base.

Crystal structure of apo Atu1419 and structural comparison

We solved the structure of Atu1419 in the apoform at 2.7 Å resolution using Tris-HCl buffer instead of citrate buffer in the crystallization condition to avoid any bound citrate (Table 1). The asymmetric unit contains a tetramer with a bound Zn 2+ ion in each effector binding domain. Asn137 and three histidine (His141, His192 and His214) side chains form the metal binding site of Atu1419 and chelate the zinc ion with a classical tetrahedral coordination geometry and average distances around 2.2 Å (Figure 5A).

We mutated the three histidine (His141, His192 and His214) into alanine to make the triple point mutant Atu1419-H3A and alter the metal binding site. The zinc ion holds helices 6, 8 and 9 involved in the tetrameric interface far from the DNA binding domain, and likely interacts with the effector as shown by the structures solved with a citrate taking the place of the effector. Atu1419-H3A is correctly folded as checked by using circular dichroism (Supplementary Figure S4) and shares the same secondary structure content (about 58% α-helices and 1% β-sheets) as the wild-type protein in line with what is observed in all crystal structures of Atu1419. Similar to the wild-type protein, this mutant is tetrameric and capable of binding atu1418-1419 and atu1420 promoter regions.

The dimeric (dimers AB/CD) and tetrameric interfaces (dimers AD/BC) covering around 1400 Å 2 and 516 Å 2 per subunit, respectively, are comparable to those of Atu1419-citrate complexes and Atu1419-DNA complexes. Subunits A and D overlap well (RMSD of 0.53 Å), as well as subunits B and C (RMSD of 0.83 Å), whereas subunits A and B or C and D show a RMSD over 2 Å, mainly due to differences in the position of their HTH motif and the end of their C-terminal domain, evidencing large flexibility (Supplementary Figure S5A). This is confirmed by the structural comparison of the DNA binding domains of the apo dimer AB with those of the Atu1419-DNA complex (Figure 5B). While subunits A are rather similar, helix 3 of the apoform subunit B clashes into the DNA sugar-phosphate backbone.

Indeed, this helix is far from its optimized position for DNA binding, which is perpendicular to the helical axis of the DNA. A displacement up to 7 Å would be required to correctly places its wHTH motif. A similar structural analysis between the DNA binding domains of the Atu1419-citrate complexes and Atu1419-DNA complexes led to similar conclusions (Supplementary Figures S5B-D).

It is noteworthy that the conformation of the effector binding pocket in all subunits bound or not to DNA and containing a citrate in orientation B is similar to that of the empty pocket observed in the apo form and one DNA complex. This suggests that this citrate may mimic a non-productive-like binding in the effector binding site.

Comparison with other FCD members

Atu1419 belongs to the VanR subgroup of the FCD subfamily of the large GntR superfamily due to the number of helices in its effector binding domain. The FCD subfamily is divided into two groups; VanR and FadR. Members that are structurally characterized in FadR subgroup possess an additional helix after the linker region at the beginning of the effector binding domain [START_REF] Van Aalten | The structural basis of acyl coenzyme A-dependent regulation of the transcription factor FadR[END_REF][START_REF] Xu | The FadR.DNA complex. Transcriptional control of fatty acid metabolism in Escherichia coli[END_REF][START_REF] Shi | The 40residue insertion in Vibrio cholerae FadR facilitates binding of an additional fatty acyl-CoA ligand[END_REF][START_REF] Gao | Structural and functional characterization of the LldR from Corynebacterium glutamicum: a transcriptional repressor involved in L-lactate and sugar utilization[END_REF]. This major difference leads to a domain-swapped quaternary structure, in which the DNA binding domain of FadR subunit A crosses the dimerization interface to be in contact with the effector binding domain of subunit B [START_REF] Gao | Structural and functional characterization of the LldR from Corynebacterium glutamicum: a transcriptional repressor involved in L-lactate and sugar utilization[END_REF][START_REF] Van Aalten | Crystal structure of FadR, a fatty acid-responsive transcription factor with a novel acyl coenzyme A-binding fold[END_REF][START_REF] Lord | McbR/YncC: implications for the mechanism of ligand and DNA binding by a bacterial GntR transcriptional regulator involved in biofilm formation[END_REF]. This is not observed in the VanR subgroup where both domains from the same subunit are associated.

Atu1419 is now the seventh VanR member to be structurally characterized. S6B).

Like Atu1419, two other VanR members, which are TM0439 and PS5454 display a metal binding site with three histidine residues and an asparagine or an aspartate residue involved in the metal coordination [START_REF] Lord | McbR/YncC: implications for the mechanism of ligand and DNA binding by a bacterial GntR transcriptional regulator involved in biofilm formation[END_REF][START_REF] Zheng | Structure of Thermotoga maritima TM0439: implications for the mechanism of bacterial GntR transcription regulators with Zn2+-binding FCD domains[END_REF]. These histidine residues are conserved and structurally close (Supplementary Figure S6C). LldR (PDB 2DI3) from the FadR group binds a zinc ion chelated by these three conserved histidines and an aspartate [START_REF] Gao | Structural and functional characterization of the LldR from Corynebacterium glutamicum: a transcriptional repressor involved in L-lactate and sugar utilization[END_REF]. It was proposed that an additional subgroup within the FCD subfamily could be created for regulators from both FadR and VanR members, which are capable of binding to metal ions [START_REF] Gao | Structural and functional characterization of the LldR from Corynebacterium glutamicum: a transcriptional repressor involved in L-lactate and sugar utilization[END_REF][START_REF] Zheng | Structure of Thermotoga maritima TM0439: implications for the mechanism of bacterial GntR transcription regulators with Zn2+-binding FCD domains[END_REF]. Using thermal denaturation experiments by differential scanning calorimetry (DSC), Tm of 52.6°C and 34°C were measured with the wild-type Atu1419 and Atu1419-H3A mutant respectively (Supplementary Figure S7). The drastic difference of almost 20°C indicates that the metal assists in stabilizing the structure affording an explanation as to why Atu1419-H3A could not be concentrated over 50 M without precipitating, in contrast to the wild-type protein (140 M). All attempts to crystallize Atu1419-H3A were unsuccessful.

Atu1419 is the first VanR member to be structurally characterized in complex with DNA whereas several FadR-DNA complexes are available [START_REF] Van Aalten | The structural basis of acyl coenzyme A-dependent regulation of the transcription factor FadR[END_REF][START_REF] Xu | The FadR.DNA complex. Transcriptional control of fatty acid metabolism in Escherichia coli[END_REF][START_REF] Shi | The 40residue insertion in Vibrio cholerae FadR facilitates binding of an additional fatty acyl-CoA ligand[END_REF]. Each dimer of Atu1419 within the tetramer can bind an identical palindrome. Similarly to regulators with a wHTH motif such as MarR regulators [START_REF] Otani | The activity of CouR, a MarR family transcriptional regulator, is modulated through a novel molecular mechanism[END_REF][START_REF] Cogan | Structural basis of transcriptional regulation by CouR, a repressor of coumarate catabolism, in Rhodopseudomonas palustris[END_REF], FCD regulators bind to one half-site of the palindromic DNA, with the dimerization interface helping to establish the spacing between the two half-sites. The DNA binding domains of Atu1419, which are involved in the dimerization interface, allow Atu1419 to bind a short palindrome of 10 base pairs. As in E. coli FadR (PDB 1H9T [START_REF] Van Aalten | The structural basis of acyl coenzyme A-dependent regulation of the transcription factor FadR[END_REF] and 1HW2 [START_REF] Xu | The FadR.DNA complex. Transcriptional control of fatty acid metabolism in Escherichia coli[END_REF]) and in Vibrio cholerae FadR (PDB 4P9U(34)), the DNAbinding site of each Atu1419 monomer is recognized by conserved residues such as Glu34 in helix 2, Arg45 (specific DNA contacts), Thr46 and Arg49 (phosphate backbone contacts) in helix 3 that interact with the major groove and Gly66 in the tip of the wing (residues 64-68) between 1 and 2. Nonetheless, the wing is not crucial for DNA interaction of Atu1419 unlike FadR, which recognizes a pseudopalindrome of 17 base pairs. Indeed, His65, which specifically recognizes a DNA base in FadR is a glycine in Atu1419. DNA binding domains superposition of the Atu1419-DNA complex with those of E. coli and V. cholerae FadR-DNA complexes shows that although the positions of subunits A are rather similar, the second subunit of Atu1419 makes steric clashes with the DNA bound to FadR and that of FadR is too far from the DNA bound to Atu1419 (Supplementary Figure S8). Atu1419 displays only one specific interaction shared in FadR with the guanine present within the major groove of DNA via their conserved Arg45, a typical feature of the GntR family [START_REF] Jain | Allosteric control of transcription in GntR family of transcription regulators: A structural overview[END_REF]. The other specific DNA contacts of FadR (via Arg35, Thr44, Thr46 and His65) are not comparable with those of Atu1419. The recognized DNA sequence of Atu1419 which is 5'-ATGTATACAT-3' is in agreement with the predicted DNA signature for the GntR family: 5'-(N)yGT(N)xAC(N)y-3' where the number x and y vary [START_REF] Rigali | Subdivision of the helix-turn-helix GntR family of bacterial regulators in the FadR, HutC, MocR, and YtrA subfamilies[END_REF]. Later, this DNA binding signature was suggested to be modified as 5'-TNG(N)nCNA-3' on the basis of base-specific interactions [START_REF] Jain | Allosteric control of transcription in GntR family of transcription regulators: A structural overview[END_REF] or as 5'-TGGTNxACCA-3' for FadR subgroup [START_REF] Gao | Structural and functional characterization of the LldR from Corynebacterium glutamicum: a transcriptional repressor involved in L-lactate and sugar utilization[END_REF] but both sequences are not appropriate for Atu1419.

The DNA binding of Atu1419 requires conformational changes of the DNA binding domains. Indeed, the distance between the two DNA recognition helices (Thr46 of helix 3) narrows from 12 Å in the apoform to 9.3 Å in the Atu1419-DNA complex, for a productive interaction with DNA (Figure 5C). A similar observation was reported for Vibrio cholerae FadR (PDB 4P9U [START_REF] Shi | The 40residue insertion in Vibrio cholerae FadR facilitates binding of an additional fatty acyl-CoA ligand[END_REF]). The structures of FadR repressors in complex with their effector show significant conformational changes transmitted from the effector binding domain to the DNA binding domain, leading to a conformational state that is no longer favorable for interaction with DNA [START_REF] Jain | Allosteric control of transcription in GntR family of transcription regulators: A structural overview[END_REF][START_REF] Van Aalten | The structural basis of acyl coenzyme A-dependent regulation of the transcription factor FadR[END_REF][START_REF] Shi | The 40residue insertion in Vibrio cholerae FadR facilitates binding of an additional fatty acyl-CoA ligand[END_REF]. The fundamental process wherein the binding of a ligand or effector molecule alters the activity of the protein at a distant site is defined as an allosteric mechanism. In the case of Atu1419, the citrate in orientation A bound in the effector binding pocket induces large conformational changes from the C-terminal domain to the N-terminal domain. It affects four different protein regions: the wing of the DNA binding domain and three regions of the effector binding domain including the C-terminus protein which is in contact with the DNA binding domain (Figure 3C). The presence of this citrate clearly reveals the plasticity of the effector binding pocket and the large potential rearrangement of the repressor. It also shows that within the tetramer, one subunit can adopt a different conformation around the effector binding pocket and further away (Figure 3C) compared to the three others subunits, highlighting an allosteric mechanism. The orientation A of the citrate may mimic an active-like position of the physiological effector part bound to the metal ion.

Atu1419 effector is the N5,N10-methylenetetrahydrofolate (MEF)

The effector molecule bound to the members of the GntR superfamily are often related to catabolic substrates or intermediates of the pathway controlled by the transcription factor [START_REF] Hoskisson | Chapter 1: Variation in form and function the helix-turnhelix regulators of the GntR superfamily[END_REF][START_REF] Rigali | Subdivision of the helix-turn-helix GntR family of bacterial regulators in the FadR, HutC, MocR, and YtrA subfamilies[END_REF], and several molecules are produced during the ferulic acid degradation pathways [START_REF] Campillo | Analysis of hydroxycinnamic acid degradation in Agrobacterium fabrum reveals a coenzyme A-dependent, beta-oxidative deacetylation pathway[END_REF]. Because the citrate molecule was able to bind to the effector binding site of Atu1419, we searched for molecules resembling citrate among those described in the HCA pathway (Figure 1A), and we asked whether the substrate-cofactor H4F or the product-cofactor M4HF of the O-demethylase Atu1420 enzyme could be a potential effector.

Indeed, citrate compound with three carboxylate groups mimics two carboxylate groups located at one end of the H4F or MH4F molecules. Alternatively, MEF, the putative product of Atu1418 enzyme, which possesses the same two carboxylate groups could also be a potential effector of Atu1419.

To verify this, Patu1418, Patu1419 and Patu1420 plasmid reporter gene fusions were introduced into the C58Δatu1418 and C58Δatu1420 mutant strains. For each strain, the fluorescence level was monitored in the presence or absence of ferulic acid and compared to the wild-type A. fabrum C58 (Figures 6A-C). For the C58Δatu1418 strain, the overproduction of fluorescence observed in the presence of ferulic acid is significantly lower than that in the wild-type C58 and C58Δatu1420 strains.

These results indicate that the presence of atu1418 gene is important for the full induction of atu1418, atu1419 and atu1420 genes expression in the presence of ferulic acid. The accumulation of MH4F and the lack of MEF production in the defective C58Δatu1418 mutant suggests that MEF could be the effector of Atu1419. MEF is not stable in aqueous in vitro solutions as it enters a rapid reversible equilibrium into H4F [START_REF] Danenberg | Folates as adjuvants to anticancer agents: Chemical rationale and mechanism of action[END_REF][START_REF] Jagerstad | 5,10-Methylene-tetrahydrofolate dissociates into tetrahydrofolate and formaldehyde at physiological pH and acidic pH, typical conditions used during sample extraction and LC-MS/MS analysis of biological samples[END_REF]. Nonetheless, immediately using freshly-made MEF solution for gel-shift assays, we were able to show that MEF could relieve repression of Atu1419 (Figure 6D) whereas H4F and MH4F had no effect (Figure 6E). In the C58Δatu1418 mutant, induction of Patu1418, Patu1419 and Patu1420 was not totally abolished, probably due to functional redundancy of N5,N10-methylene tetrahydrofolate reductase activity in A. fabrum genome, which allows to keep a basal pool of folate and intermediates essential for the cell [START_REF] Feirer | A Pterin-Dependent Signaling Pathway Regulates a Dual-Function Diguanylate Cyclase-Phosphodiesterase Controlling Surface Attachment in Agrobacterium tumefaciens[END_REF].

The interaction between MEF and Atu1419 was confirmed by DSC with a Tm of 63.5°C for Atu1419 in complex with MEF compared to the Tm of 52.6°C for the apo Atu1419 (Figure 6F). This was also confirmed by isothermal titration microcalorimetry (ITC, Figure 6G). Because MEF instability did not allow a return to the baseline between two injections, we could not rigorously determine a dissociation constant KD value, which is around 20 M. Injection of MEF towards buffer solution was a control of its instability effect (Figure 6H). No interaction could be observed between MH4F and Atu1419 by DSC and ITC, in agreement with the results from EMSA assays.

Atu1419 tetramer is the biologically active form

We analyzed the quaternary structure of Atu1419 in complex with MEF using gel filtration measurement (Supplementary Figure S9). Atu1419 remains tetrameric when bound to MEF. Thus, Atu1419 is a tetramer in its apo form, upon DNA binding or effector binding suggesting that this quaternary structure is the biologically active form. Both crystal structures of Atu1419 in complex with DNA revealed that each dimer of Atu1419 within the tetramer can bind a palindromic site of 10 bp, separated by approximately 10 nM (Figure 3A). This separation distance is compatible with DNA loop formation allowing the simultaneous binding of two distant palindromes (190 bp apart) within the intergenic region of atu1418-1419. Thus, Atu1419 tetramer could repress gene expression via DNA looping, such as several negative regulators in prokaryotes [START_REF] Bylino | Evolution of Regulated Transcription[END_REF] with the example of the intensively studied tetrameric lac operon repressor [START_REF] Lewis | Crystal structure of the lactose operon repressor and its complexes with DNA and inducer[END_REF].

We searched for a second potential DNA binding site of Atu1419 in atu1420 region and found a degenerated palindrome (8 bp conserved over the 10 bp palindrome; Figure 7A). The degenerated palindrome is separated by 109 bp from the palindrome and is located downstream within the open reading frame of atu1420 gene (Figure 7A). We then performed EMSA assays with three oligonucleotides containing either the palindrome, the degenerated palindrome, or both and identified in the gel a retarded band for each oligonucleotide proving the formation of a stable non-covalent protein-DNA complex (Figure 7B). Thus, the repressor can recognize a palindrome and a degenerated palindrome providing it the opportunity to bind simultaneously to two DNA sites in each region (atu1418-1419 and atu1420), in agreement with its oligomeric state. Similarly, the structure of Cupriavidus necator CbnR transcriptional regulator (a tetramer composed of dimer of dimer) from LysR family (PDBs 1IXC and 1IZL (45)) is compatible with an interaction with two DNA binding sites on a bended DNA fragment, where the authors suggested that the quaternary structure of the CbnR tetramer should be changed to relax the bent DNA that binds to the CbnR tetramer [START_REF] Muraoka | Crystal structure of a full-length LysR-type transcriptional regulator, CbnR: unusual combination of two subunit forms and molecular bases for causing and changing DNA bend[END_REF]. In contrast to the quaternary organization of Lac repressor (a tetramer comprising a pair of dimers loosely associated with each other) which suggests that the structural change of each dimer seems to be able to occur independently, conformational changes of one subunit of Atu1419 and CbnR can be easily propagated to the other parts of the tetramer though inter-subunit interactions. Noteworthy, single molecule experiments revealed that a repressor securing a DNA loop is a stronger transcriptional roadblock for RNA polymerases than one bound to a single DNA site [START_REF] Voros | Proteins mediating DNA loops effectively block transcription[END_REF].

Finally, we investigated the binding between the tetrameric Atu1419 and the Patu1418-1419 region of 370 bp containing two palindromes (190 bp apart) using negative staining electron microscopy. We measured the size of the protein (~10nM) and the length of the DNA fragment (~125 nM), which are in good agreement with the expected values (10 nM and 129 nM, respectively Figure 8). We were able to observe three categories of binding: a tetramer bound to one DNA fragment, two tetramers bound to one DNA fragment and a tetramer forming a DNA loop (Figure 8).The loops were measured at ~60-65 nM compatible with the length of 190 bp, which separates two palindromes.

CONCLUSION

Hydroxycinnamic acids are involved in the initiation of the plasmid-encoded pathogenicity program of A. fabrum and can be degraded through the expression of the species-specific chromosomal gene cluster composed of eight genes, which are regulated by two repressors HcaR/Atu1422 and Atu1419.

We have previously shown that HcaR, a dimeric protein from the MarR family, is the transcriptional repressor of its own transcription, atu1421 gene and that of the first three genes (atu1415, atu1416 and atu1417) of HCA degradation pathway. Here, we characterized Atu1419 as a transcriptional repressor of its own repression in addition to two other genes (atu1418 and atu1420). We identified the short palindromic region composed of 10 base pairs and a degenerated palindrome bound to Atu1419 and the N5,N10-methylenetetrahydrofolate (MEF) as the effector molecule disrupting Atu1419-DNA interaction, allowing the expression of the second part of the HCA degradation pathway. Remarkably, Atu1419 repressor is not regulated by the direct catabolic substrates/intermediates of the pathway. An unexpected outcome was that Atu1419 is a tetrameric regulator whereas all GntR regulators so far studied have been shown to be dimeric [START_REF] Rigali | Subdivision of the helix-turn-helix GntR family of bacterial regulators in the FadR, HutC, MocR, and YtrA subfamilies[END_REF][START_REF] Suvorova | GntR Family of Bacterial Transcription Factors and Their DNA Binding Motifs: Structure, Positioning and Co-Evolution[END_REF][START_REF] Jain | Allosteric control of transcription in GntR family of transcription regulators: A structural overview[END_REF]. This quaternary structure allows the repressor to bind two distant DNA sites simultaneously making DNA loop repression (Figure 8) and providing a fine-tune mechanism in the transcriptional regulation of the appropriate genes at the right time. Our structural and biophysical study of Atu1419 revealed different structural rearrangements around the fortuitous citrate bound in the effector binding domain within the tetramer and a double role of the Zn 2+ ion in protein stability and effector binding. This suggests an induced-allosteric mechanism by MEF, involving conformational changes "in cascade" which enhance interactions between the DNA binding domain and the end of the effector binding domain to constrain the flexibility of the HTH motif to prevent it from DNA binding.

ACCESSION NUMBERS

The atomic coordinates and structure factors have been deposited at the Protein Data bank under PDB ID 6Z74 and 6ZA0 for Atu1419 in complex with a fortuitous citrate in P212121 and P21212, respectively, Citrate is shown in its Fo-Fc omit map contoured at 4σ in subunit A.

FIGURES LEGENDSFigure 1 .

 1 FIGURES LEGENDS Figure 1. (A) Ferulic acid degradation pathway: a coenzyme A is added to ferulic acid by Atu1416, a feruloyl-CoA synthase. Feruloyl-CoA is then converted by Atu1417, an enoyl-CoA hydratase, into 4hydroxy-3-methoxyphenyl-β-hydroxypropionyl (HMPHP)-CoA, which is in turn transformed into vanillic acid by Atu1415, a phenylhydroxypropionyl-CoA dehydrogenase, and then Atu1421, a 4-hydroxy-3methoxyphenyl-β-ketopropionyl-CoA (HMPKP)-CoA β-keto-thiolase. The O-demethylase Atu1420 degrades vanillic acid into protocatechuic acid using tetrahydrofolate (H4F) as a cofactor and produces N5-methyl-tetrahydrofolate (MH4F). Protocatechuic acid enters the cycle of Krebs. Atu1418 enzyme was proposed to be involved in the recycling of H4F from transformation of MH4F into MEF, a compound that can be enzymatically or spontaneously converted to H4F (18,50). Dashed arrows indicated putative reactions based on sequence similarities and references (4,18,44,50). (B) Organization of the SpG8-1b region for which genes expression is regulated by the two transcriptional factors HcaR and Atu1419 (their genes are framed in red). Promoters are shown with small arrows in black for HcaR dependenttranscription and in blue for putative Atu1419 dependent-transcription indicating the direction of gene transcription.

Figure 2 .

 2 Figure 2. Atu1419 regulation mechanism. (A) Gel shift analysis of 30 nM of each promoter region Patu1416-1417, Patu1418-1419, Patu1420, PhcaR (atu1422) incubated without and with 250 nM tetrameric Atu1419 (ratio protein:DNA of 8.3). PvirB was used as a promoter region control. (B) In silico analysis of the atu1418-1419 and atu1420 promoter regions performed with the BPROM program (31): -10 and -35 boxes are underlined, and the palindromic sequences are shown in red. The transcription initiation sites are shown with an arrow and the translation start sites are indicated in bold. Palindromic regions are separated by 190 bp in atu1418-1419 region. (C) Gel mobility shift assay analysis of 30 nM of three parts of the promoter region of atu1418-1419 containing each one palindrome (P1 or P2) or none (P3) incubated without and with 250 nM tetrameric Atu1419 (ratio protein:DNA of 8.3).

Figure 3 .

 3 Figure 3. (A) Wire representation of Atu1419 tetramer in complex with a fortuitous citrate molecule and a co-purified Zn 2+ ion, shown as sticks and ball, respectively, bound to the effector binding domain in the structure at 2 Å resolution (molecule A is in cyan whereas molecules B/C/D in gray). (B) Left: ribbon representation of dimer AB within the tetramer in A. The secondary structural elements are indicated in subunit A and the orange elements form the dimeric interface. The linker joining the N-terminal DNA binding domain and the C-terminal effector binding domain is shown in black. Right: ribbon representation of dimer AD within the tetramer in A. The secondary structural elements are indicated in subunit A and the orange elements form the other dimeric interface (C) Superposition of the subunits A in cyan and B in gray of the tetramer in A (the subunits B, C, D adopt the same fold). The red parts indicate major conformational changes between both subunits. A close-up view of the citrate and ion binding sites in the effector/metal binding domain shown in blue and orange for subunits A and B, respectively. (D) Interactions between the bound citrate in orientation A and in blue/Zn 2+ (blue ball) and subunit A in cyan. Hydrogen bonds are shown as dashed lines in black (distance below 3.2 Å) and metal contacts are shown as dashed lines in red. Residues involved in the interactions are labeled and

Figure 4 .

 4 Figure 4. Cartoon representation of the tetrameric Atu1419-DNA complex (A) in the P6422 structure with the monomer of the asymmetric unit shown in purple. (B) in the P212121 structure with the dimer of the asymmetric unit shown in magenta. (C) Superposition of the DNA binding domains of the dimers AB of the P6422 structure in purple and gray for subunits A and B, respectively and of the P212121 structure in magenta. The palindromic 10-mer DNA is in green and orange in the P6422 and P212121 structures, respectively. (D) View of the DNA binding domain of the dimer AB of the P6422 structure in purple and gray for subunits A and B, respectively bound to the palindromic 10-mer DNA shown in green. The secondary elements are indicated in subunit A. Arg45 in helix 3 is the key residue, which interacts with the essential guanine base at position 3. This residue is shown as sticks. (E) View showing a close-up of the Atu1419-DNA interface within a monomer and schematic diagram of Atu1419-DNA contacts. Nucleotide bases that interact with Atu1419 are shown in green. (F) Gel mobility shift assay analysis of Patu1420 incubated without and with Atu1419, and with Atu1419 and 100 µM citrate (ratio tetrameric protein:DNA of 8.3). (G) Gel mobility shift assay analysis of Patu1420s containing a unique palindrome or Patu1420s AT, in which the palindrome was mutated with the recognized guanine at position 3 replaced by an adenine and its cytosine partner replaced by a thymine, incubated without and with Atu1419 (ratio tetrameric protein:DNA of 4.6).

Figure 5 .

 5 Figure 5. (A) Cartoon representation of the C-terminal domain except for helix 5 shown in ribbon of Atu1419 with a co-purified Zn 2+ ion and a close-up view of the Zn 2+ binding site. The four amino acid residues bound to the metal ion (orange ball) are shown as sticks. (B) Superposition of dimers AB of apo Atu1419 (in orange) and Atu1419-DNA complex (in purple and gray). The DNA is in green. (C) Close-up view showing the superposition of their DNA binding domains colored as in B with a view at 90°.

Figure 6 .

 6 Figure 6. (A) (B) (C) Patu1418, Patu1419 or Patu1420 expression in A. fabrum C58 wild-type strain compared with that in C58atu1420 and C58atu1418 mutants in the presence of ferulic acid. Values correspond to normalized fluorescence intensity in the presence of ferulic acid corrected by subtracting values without ferulic acid. Different letters indicate statistical differences between conditions (one-way ANOVA and Tukey test; P-value=0.05). (D) Gel mobility shift assay analysis of Patu1420 DNA region incubated without and with Atu1419 (ratio tetrameric protein:DNA of 17) and N5,N10methylenetetrahydrofolate at different concentrations (0 to 300 M). (E) Gel mobility shift assay analysis of Patu1420 DNA region incubated without and with Atu1419, and with Atu1419 and H4F, MH4F or MEF at 100 M. (F) Differential scanning calorimetry thermograms of apo Atu1419 (black) and Atu1419 in the presence of MEF (magenta). The table below indicates the Tm. DSC experiments were performed

Figure 7 .

 7 Figure 7. (A) atu1420 promoter region: -10 and -35 boxes are underlined, the palindromic sequence and the degenerated palindromic sequence (109 bp apart) are in red. The transcription initiation site is shown with an arrow and the translation start site is indicated in bold. (B) Gel mobility shift assay analysis of four parts of the promoter region of atu1420 containing each one palindrome or one degenerated palindrome (P1 or P2), or both (P1+P2) or none (P3) incubated without and with 250 nM or 500 nM tetrameric Atu1419 (ratio protein:DNA of 8.3 or 16.6).

Figure 8 .

 8 Figure 8. Interactions between the tetrameric Atu1419 repressor and Patu1418-1419 region of 370 bp containing two palindromes separated by 190 bp visualized by electron microscopy and corresponding model. (A) Image of Atu1419 tetramer as a control, (B) Image of Patu1418-1419 region as a control, (C) Image of one tetramer bound to one palindrome only, (D) Image of two tetramers bound each to one palindrome, (E) Image of one tetramer and a DNA loop.

Table 2 .

 2 4 Å along b* and 0.894a* -0.447b* and 2.7 Å along c* † Two datasets collected from a crystal, which diffracted anisotropically to 3.2 Å along 0.81a* -0.58c*, 2.3 Å along b* and 2.2 Å along c* § One dataset collected from a crystal, which diffracted anisotropically to 2.02 Å along a*, 1.99 Å along b* and 1.95 Å along c* SpG8 1-b genes expression in the wild-type (WT) and C58Δatu1419 strains. Values refer to fold change at 24 hours obtained by comparison of genes expression in the WT C58 strain with and without 750 µM of ferulic acid and correspond to the mean of three biological replicates with at least six technical replicates. Fold change at 24 hours without ferulic acid was obtained by comparison of genes expression in the C58Δatu1419 and the WT strains and corresponds to the mean of at least two biological replicates and six technical replicates. Empty pOT1e corresponds to plasmid without any promoting region before the egfp gene in order to measure basal expression of the system.

  and 2 of the wHTH motif, residues 119-130 between helices 5 and 6, residues 169-184 between helices 7 and 8 and residues 225-244 corresponding to the end of the effector binding domain, can move between 3 and 10 Å. The loop region 169-184 can drastically rearrange upon ligand binding. In Thirteen polar interactions are observed between monomer A and the citrate molecule (orientation A), which is buried within the monomer, leaving only 26.6 Å 2 or 8.4% of the molecule surface exposed to solvent. In contrast, there are only ten polar interactions between monomers B/C/D and the citrate with a similar buried surface area.The orientation B of the citrate was also observed in another structure of Atu1419 solved at a higher resolution of 1.75 Å but in a different space group (P21212). Here, two similar monomers (RMSD of 0.98 Å for all C atoms) are in the asymmetric unit and form a dimer, which in turn form a tetramer by the crystal symmetry with globally unchanged dimer interfaces (1553 Å 2 and 701 Å 2 per subunit for dimers

	molecule A, the C atoms of Ser169, Lys175 and Arg183 are respectively, 0.92 Å, 6.84 Å and 7.34 Å
	away from those in molecules B/C/D, allowing the side chains of Lys175, Arg183 and the NH main
	chain of Ser169 to interact with the citrate molecule (Figures 3C and 3D). These latter protein-citrate
	interactions cannot exist in subunits B/C/D. Nonetheless, all bound citrates share interactions with
	Arg86, Tyr133, Asn137, His141, His214, Asn221 and a Zn 2+ ion that co-purified with Atu1419 (Figures
	3D and 3E).

AB and AD, respectively) compared with those from the P212121 structure. They resemble subunits B/C/D of the P212121 structure and bind similarly a citrate molecule. Nonetheless, their citrates are slightly shifted by 0.6 Å towards Arg86 losing the interaction with Asn221 compared with those in the P212121 structure (Supplementary Figure

Table 1 .

 1 Crystallographic data and refinement parametersValues for the highest resolution shell are in parentheses. CC1/2 = percentage of correlation between intensities from random hall-dataset a Calculated with MolProbity Numbers in italic account for statistical values after ellipsoidal mask application by Staraniso. # A dataset collected from a crystal, which diffracted anisotropically to 3.

			Zn 2+ SAD	Atu1419	Atu1419	Atu1419-	Atu1419-	Apo Atu1419 †
				Citrate	Citrate	DNA#	DNA §
	PDB code			6Z74	6ZA0	6ZAB	6ZA3	6ZA7
	Wavelength (Å)	1.282290	0.97541	0.978570	1	1	0.97934
	Crystallization	5% PEG	5% PEG 4K, 0.2 M	5% PEG	10%	25% PEG	12% PEG 4K, 0.2
	conditions		4K, 0.2 M	AS, 0.1 M Na-Citrate	4K, 0.2 M	Terbutanol,	400, 0.1 M	Ammonium
			AS, 0.1 M	pH 5.6	AS, 0.1 M	0.1 M	MES pH 6.7,	Sulfate, 0.1 M
			Na-Citrate		Na-Citrate	NaCitrate pH	0.2 M Na-	Tris/HCl pH 8.5
			pH 5.6		pH 5.6	5.6, 2 mM	Acetate
						MgCl2
	Za		4	4	2	1	2	4
			P212121	P212121	P21212	P6422	C2221	P21
	Space group	a= 77.4	a= 77.2	a= 72.2	a= 179.5	a= 62.4	a= 52.1
	Cell parameters	b= 114,1	b= 113.7	b= 145.6	b= 179.2	b= 112.4	b= 79.0
	(Å,°)		c= 130.5	c= 130.5	c= 41.2	c= 96.8	c= 179.6	c= 140.6
								β= 95.1
	Resolution (Å)	49.8-2.35	45.6-2	24.0-1.65	45.61-2.8	47.64-2.0	43.38-2.34
			(2.42-2.35)	(2.12-2)	(1.69-1.65)	(2.87-2.8)	(2.05-2.0)	(2.41-2.34)
	No. of observed	518098	467551	239555	444041	376669	504810
	reflections		(29677)	(72154)	(17376)	(30119)	(23170)	(31752)
	No. of unique	48018	78344	51917	23089	43067	47943
	reflections		(3262)	(12280)	(3674)	(1642)	(3095)	(3255)
	Completeness	99.4 (92.1)	99.6 (98.1)	97.9 (95.8) 99.7 (88.4)	99.9 (99.1)	99.4 (92.3)
	(%)
	Completeness				68.7 (11.9)	97.7 (79.3)	71.5 (30.6)
	Staraniso (%)
	Rsym (%)		8.8 (78.6)	12 (155.5)	8.1 (127.3) 11.6 (236.4) 8.6 (190.3)	10.9 (128.8)
	Rpim (%)		2.8 (26.6)	5.4 (78.1)	4.2 (64.4)	2.7 (57.7)	3.1 (79.3)	3.5 (41.5)
	I/(I)		18 (2.3)	10 (1.12)	10.7 (1.1)	18.8 (1.1)	13.7 (0.9)	13.0 (1.4)
	CC1/2		0.999	0.998 (0.582)	0.999	0.999 (0.469) 0.999 (0.421)	0.999 (0.696)
			(0.809)		(0.461)
	Rcryst (%)			19.1	16.9	18.8	19.2	19.4
	Rfree (%)			21.4	18.9	19.9	22.8	22.2
	rms	bond		0.01	0.001	0.009	0.01	0.009
	deviation (Å)
	rms	angle		0.99	0.95	1	1	1.02
	deviation (°)
	Average B (Å 2 )
	Protein A/B/C/D		42/53.1/60/57.5	27.9/29.2	86.4	50.7/50.1	66.2/74.3/73.2/80.8
	Zn 2+			35/41.4/39.7/45.2	18.3/18.6	106.9	43.3/46.2	46.2/63.2/69.3/66.6
	citrate			44.4/45.2/44.5/45.7	23.1/22	127.9
	DNA					65.3	60.2
	solvent			50.5	36.7	56.3	48.	50
	a Clashscore		4.56	1.6	3.2	2.11	2.88
	MolProbity score		1.51	0.93	2.07	1.24	1.38
	a Ramachandran
	plot (%)
	Favoured			99.36	98.73	94.94	98.73	98.2
	Outliers			0	0	0	0	0.21

Figure 1.

Figure 2.

Figure 3.

Figure 4.

Figure 5.

Figure 6.

FUNDING

AV, SM, TM, LV, CL was supported by CNRS (Mission pour l'interdisciplinarité, Agromics 2014-2016).

This work was submitted to fulfill the requirements for a doctorate of biology at ED341-E2M2 from Université de Lyon, granted from the French Ministère de l'Education Nationale, de l'Enseignement Supérieur et de la Recherche (T.M). This work has benefited from the I2BC crystallization and Proteinprotein interactions platforms supported by FRISBI ANR-10-INSB-05-01.

SUPPLEMENTARY DATA

Supplementary Data are available at NAR online: Supplementary Figures S1 to S9 and Tables S1 to S3.

ACKNOWLEDGEMENT

We acknowledge SOLEIL for provision of synchrotron radiation facilities (proposals ID 20130869, 20140774, 20150780 and 20170872) in using Proxima beamlines. We are grateful to Nadine Assrir (ICSN, Gif sur Yvette) and Pascal Retailleau (ICSN, Gif sur Yvette) for access to the CD equipment and help in data processing. We are grateful to Christophe Velours for performing SEC-MALS experiment on the Protein-protein interactions platform. We thank the Cell and Tissue Imaging (PICT IBiSA), Institut Curie, member of the French National Research Infrastructure France-BioImaging (ANR10-INBS-04).

We thank Andrew Saurin for a critical reading of the manuscript.

CONFLICT OF INTEREST

The authors declare that they have no competing interests.


Amélie Deglaire

Anne Blais

Juliane Calvez

Géraldine Lucchi

Karine Gourrat

Nadine Leconte

Romain Jeantet

Didier Dupont

Sophie Nicklaus

Anne-Marie Davila

The production of an infant formula with a minimally processed route impacts its nutritional, physiological and sensorial qualities

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

The production of an infant formula with a minimally processed route impacts its nutritional, physiological and sensorial qualities.

Amélie Deglaire 1 , Anne Blais 2 , Juliane Calvez 2 , Géraldine Lucchi 3 , Karine Gourrat 3 , Nadine Leconte 1 , Romain Jeantet 1 , Didier Dupont 1 , Sophie Nicklaus 3 & Anne-Marie Davila 2 . Many infants receive infant formulas (IFs), the most nutritionally adequate substitute to human milk. However, health effects still differ between these two foods. IF production requires multiple unit operations, particularly heat treatments leading among others to protein denaturation/aggregation. Recently, a minimal processing route (low heat load) allowed the semi-industrial production of a bacteriologically safe IF, with a high content of native proteins and stable over time. The objective of the present study was to evaluate its organoleptic, nutritional and physiological properties as compared with a high heat load IF.

Two complete powdered IFs were produced from microfiltered bovine milk and whey protein concentrate, with no heat treatment (T-) or successive ones (T+++), resulting in a protein denaturation extent of 4.4±0.4% and 60.0±0.1%, respectively. Their organoleptic properties were determined by aroma analysis (gas chromatography / mass spectrometry) and sensory evaluation in human adults (n =50). The digestive deconstruction of IFs was followed through in vitro dynamic digestion. The true protein digestibility and the impact on intestinal homeostasis and microbiota were determined in young rats (n = 6-8).

The organoleptic quality of the IFs significantly differed (triangular test, p<0.0001), with a higher aroma compound concentration in T-than in T+++. The IFs presented a different microstructure, leading to larger aggregates for T-during its in vitro gastric digestion. T-had a higher in vivo true protein digestibility than T+++ (96.6±5.4% vs. 91.9±5.4%, p-value=0.09) such as for plasma amino acid concentration 1-hour postprandially (p-value<0.001). Tincreased dietary intake and weight gain, decreased the caecal ammonium content, and modified the caecal microbiota species in rats.

This study demonstrated that the minimally processed IF (T-) had different organoleptic, nutritional and physiological properties, including an enhanced protein bioavailability. Such innovative IF may be a further step towards human milk biomimetics.

Quote : It is feasible to produce a native and safe infant formula with a minimally processed route. At the conference, I will present how this impacts its organoleptic, nutritional and physiological properties.