{'fraction_non_alphanumeric': 0.03696863177336233, 'fraction_numerical': 0.028597813458845957, 'mean_word_length': 4.773993808049536, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
First description of betalains biosynthesis in an aquatic organism: characterization of 4,5‐DOPA‐extradiol‐dioxygenase activity in the cyanobacteria Anabaena cylindrica

This paper describes molecular and functional characterization of a novel dioxygenase enzyme from the aquatic cyanobacterium Anabaena cylindrica. The enzyme showed a strong inhibition by excess of the precursor L‐DOPA. However, its heterologous expression has allowed detecting the formation of the main compounds in the biosynthetic pathway of betalains. The discovery of this novel dioxygenase in the phylum cyanobacteria shows that betalains' formation is more diverse than expected.


Introduction
Betalains are a class of hydrophilic, nitrogen-containing pigments present in plants of the order Caryophyllales which lack the most common pigments anthocyanins (Gand ıa-Herrero and Garc ıa-Carmona, 2013). These two families of pigments are mutually exclusive, and they have never located together in the same plant (Castellanos-Santiago and Yahia, 2008;Brockington et al., 2011). Betalains have high antioxidant and free radical scavenging activities in vitro (Gand ıa-Herrero et al., 2010) and they show a health-promoting effect in vivo, expanding the lifespan of Caenorhabditis elegans (Guerrero-Rubio et al., 2019a). Several studies have probed its importance as chemopreventive molecules in different lines of cancer cells (Sreekanth et al., 2007) and a reduction of tumours in mice after their intake (Lechner et al., 2010). In addition, betalain-rich extracts have a positive activity in humans (Pietrzkowski et al., 2010;Rahimi et al., 2019) which is probably related to the antioxidant and anti-inflammatory activities of betalain pigments (Vidal et al., 2014).
Our last studies revealed that betalamic acid is not only present in plants of order Caryophyllales and fungi Amanita and Hygrocybe but it is also present in bacteria as Gluconacetobacter diazotrophicus, the first betalainproducing bacterium described (Contreras-Llano et al., 2019). G. diazotrophicus contains the smallest DODA enzyme (17.8 kDa) able to produce betalamic acid and muscaflavin. Besides, the high stability of this enzyme has led its heterologous expression in a recombinant system, the fully characterization and time evolution of intermediates 2,3-seco-DOPA and 4,5-seco-DOPA and the development of biofactories for betalains' production (Guerrero-Rubio et al., 2019b).
In the present study, we expand the range of betalainproducing bacteria describing the cloning, expression, purification, molecular and functional characterization of a DOPA-extradiol-dioxygenase enzyme from Anabaena cylindrica, the first cyanobacterium described producing betalamic acid.

Results and discussion
Anabaena cylindrica cultures produce betalamic acid The discovery of Gluconacetobacter diazotrophicus as the first betalain-forming bacterium (Contreras-Llano et al., 2019) has opened a new field in the search for novel biological systems able to produce betalains. Previous phylogenetic analysis showed sequence homology for the dioxygenase enzyme of Gluconacetobacter in several proteobacteria, such as Bradyrhizobium, Komagataeibacter, Mesorhizobium and Inquilinus. An extended phylogenetic analysis performed in this work has shown the presence of homologous proteins not only in proteobacteria genus but also in the phylum cyanobacteria (Fig. S1). Within the sequences found in this analysis, several genera belonging to the order Nostocales, such as Anabaena, Nostoc, Trichormus or Aphanizomenon, showed sequence similarity with the dioxygenase enzyme of G. diazotrophicus. Among them, A. cylindrica, an autotrophic filamentous bacterium with nitrogen-fixing capabilities, was identified and chosen as a representative species of the phylum cyanobacteria to study the possible betalamic acid-producing abilities in this group of water organisms. Anabaena cylindrica is a worldwide-distributed species which is well known for its use in nitrogen-fixing studies since many decades (Allen and Arnon, 1955). However, there is no evidence on the existence of the betalains pathway in this species and a betalamic acid-forming activity has never been described. Thanks to the interest in this aquatic microorganism its genome has been completely sequenced (Shih et al., 2013). A. cylindrica cultures enriched with increasing concentrations of L-DOPA up to its solubility limit at 7.6 mM were employed. Preliminary assays of A. Fig. 1. Scheme of activity of DOPA-extradiol-dioxygenase (DODA) enzymes. DODA enzymes catalyse the cleavage of L-DOPA between C 2 and C 3 , yielding 2,3-seco-DOPA, or between C 4 and C 5 , yielding 4,5-seco-DOPA. Both compounds are intermediates which, by spontaneous cyclization, give rise to muscaflavin and betalamic acid respectively. The condensation of betalamic acid with L-DOPA yields dopaxanthin.
ª 2020 The Authors. Microbial Biotechnology published by John Wiley & Sons Ltd and Society for Applied Microbiology., Microbial Biotechnology, 13, 1948-1959 cylindrica cultures supplemented with L-DOPA as a precursor in the formation of betalamic acid did not produce any yellow pigmentation, but they produced the darkening of the medium. It was hypothesized that this was due to the production of DOPA-chrome, a product derived from L-DOPA in the production of melanins by oxidase enzymes. Previous studies have reported that ascorbic acid or an analogous reducing agent is necessary in the biosynthesis pathway to recover the pool of L-DOPA, thus avoiding the accumulation of DOPAchrome (Gand ıa-Herrero and Garc ıa-Carmona, 2013). The addition of sodium ascorbate 15 mM to the media did not produce the expected results, and the conversion of L-DOPA to DOPA-chrome was again evident from the darkening of culture media. Then, it was considered that this might be due to the presence of an ascorbate peroxidase particularly active in A. cylindrica, demonstrated to be involved in the scavenging of hydrogen peroxide (Miyake et al., 1991). This enzyme is able to capture sodium ascorbate molecules, by reducing them and limiting their bioavailability to revert the conversion of L-DOPA to DOPA-chrome. Taking these considerations into account, EDTA 100 mM was also added to the reaction media in order to avoid the reduction of sodium ascorbate. 7-day cultures of A. cylindrica were supplemented with L-DOPA at concentrations 0.23, 0.76, 1.90, 3.80 and 7.60 mM, 15 mM sodium ascorbate and 100 mM EDTA. Non-supplemented cultures were employed as controls. After two days, supernatants were analysed by HPLC-ESI/TOF/ MS and pure pigment dopaxanthin was employed as a standard to check out the activity of a hypothetical DODA enzyme. This betaxanthin is the result of the condensation of betalamic acid with L-DOPA, and its presence reveals the formation of betalamic acid from L-DOPA through the presence of a 4,5-DOPA-extradioldioxygenase enzyme. Control pigment analysis yielded a retention time of 12.1 min with an exact mass value of 391.1136 m/z (Fig. 2). The same peak was obtained in cultures of A. cylindrica supplemented with 3.8 and 7.6 mM of L-DOPA but it was not found at lower concentrations of the added substrate neither in its absence. Therefore, the existence of a betalamic acid-forming DODA activity in A. cylindrica has been demonstrated. This metabolism does not seem to be active or relevant in the development of A. cylindrica but it is able to produce betalamic acid when the precursor amino acid is present in the medium. The betalamic acid produced condenses with L-DOPA to yield its derived betaxanthin, dopaxanthin. A regulatory system for the 4,5-DOPA-extradiol-dioxygenase activity or the biosynthesis of betalains in microorganisms has not been described. However, in A. cylindrica it is here demonstrated that ascorbate peroxidase plays a role in limiting the access of the main enzyme in the pathway to the initial substrate L-DOPA. Betalamic acid and the final pigment dopaxanthin can only be produced by the cyanobacteria when the activity of the peroxidase enzyme is lowered. This may point out to a yet unknown regulatory system that control expression and activity of betalamic acid-forming dioxygenases in microorganisms. In plants, expression of 4,5-DOPA-extradiol dioxygenase is highly controlled and it only appears in specific organs and developmental stages (Sasaki et al., 2009;Bean et al., 2018).
A. cylindrica 4,5-DODA sequence, expression and purification The search for non-plant L-DOPA dioxygenases able to produce betalamic acid -and thus to start the biosynthetic pathway of betalains in prokaryotes -has led us to search for sequences related to that described in G. diazotrophicus. Phylogenetic analysis shows a sequence susceptible to produce betalains in the genome of A. cylindrica. The genome of this autotrophic cyanobacterium is deposited in the NCBI database with the access number NZ_AP018166.1 (project PRJDB5665). Its chromosome presents a size of 7.06 Mb where the sequence for the protein WP_015213489.1, corresponding to the 4,5-DODA sequence, is located with the locus tag CA722_RS27340. The sequence for the protein WP_015213489.1 is accessible from the databases and was here used as a template to obtain a synthetic sequence optimized for expression in E. coli. This optimization introduced mutations which did not change the amino acid sequence coded. The modified sequence was inserted into the multiple cloning site of the expression vector pRSETA, obtaining the recombinant plasmid pRSETA-AcDODA.
This plasmid was used to transform thermocompetent E. coli BL21 (DE3) cells. For an increased expression of activity in transformed E. coli cultures, optimization conditions of IPTG concentration and temperature of induction was performed using IPTG at final concentrations of 0.1, 0.5 and 1 mM and temperatures of 20 and 25°C. In all cases, 20 h after induction L-DOPA 7.6 mM and ascorbate 15 mM were added to the cultures and dopaxanthin formation was evaluated by HPLC. This measure allowed to establish the optimal conditions for the activity, and thus protein expression. The highest dopaxanthin content was achieved after induction with 0.1 mM IPTG at 20°C, as shown in Fig. S2. These conditions were then employed for protein purification.
The purified enzyme was run on an SDS-PAGE electrophoresis gel, revealing that the recombinant A. cylindrica DODA was purified to homogeneity after affinity separation. A higher amount of pure protein was obtained when the protein solution was eluted at pH 7.5 (Fig. 3A). The results of the purification process are shown in Table 1, where it can be seen that a purification yield of 75% and a purification fold of 1.6 were obtained, in accordance with the electrophoretic estimation of DODA expression.

Molecular characterization
The band corresponding to the recombinant protein obtained by SDS-PAGE yields an estimated molecular mass of 17 kDa (Fig. 3), consistent with the molecular weight calculated according to the protein sequence (17.8 kDa). The purified protein was analysed by HPLC-ESI-MS-TOF in order to determine its accurate molecular mass. Mass spectra showed a single peak with a molecular mass of 17.806 kDa, which agrees with the above-mentioned values obtained by SDS-PAGE and for the theorical molecular weight. To fully characterize the protein, its peptide mass fingerprint was determined by MALDI-TOF analysis after trypsin digestion. The main peptides identified corresponded to the The purified protein was subjected to gel filtration chromatography on a Superdex 200 column under native conditions in order to determine whether this DODA enzyme is a monomer or if it forms oligomers. In this sense, different concentrations of purified protein were employed to investigate the possibility of an equilibrium among the monomer and the possible oligomers. A single peak eluting at the same volume was obtained in all cases, with a molecular mass estimated of 32.9 kDa (Fig. 3). This value doubles the expected mass These results strongly support the sequence homology found by a search performed on the databases for sequences with known structures. The A. cylindrica DODA has 71.3% identity (83.0% similarity) with a putative dioxygenase (PDB entry 2PEB) enzyme characterized from Nostoc punctiforme PCC 73102 (additional information can be found in Appendix S1). X-ray crystallography at 1.46 A resolution has demonstrated that this protein is a dimer composed of two monomers of 122 amino acids each. The similarity between these two enzymes agrees with the molecular characterization performed for AcDODA in the present study. In addition, A. cylindrica DODA residues can be assigned to specific secondary motifs also found in the enzyme from the only prokaryote known to produce betalains, G. diazotrophicus, by structurally assisted sequence comparison (Armougom et al., 2006;Robert and Gouet, 2014). Figure 3C shows the structural similitude between both enzymes with respect to the residues involved in the formation of b-strands and a-helixes. Taking into account the molecular similarities between the A. cylindrica DODA and G. diazotrophicus enzymes and with the protein crystalized from N. punctiforme, a three-dimensional model of the novel DODA from A. cylindrica was performed (Biasini et al., 2014;Bienert et al., 2017). The resulting structure can be found in Fig. 3D where a dimer modelled by the bioinformatic techniques is shown. The predicted dimeric protein finds experimental support in Fig. 3B where the multimeric nature of the novel enzyme is unambiguously determined.

Kinetic characterization
The activity of the recombinant DODA enzyme from A. cylindrica was characterized spectrophotometrically by its addition to a reaction medium with the amino acid L-DOPA as substrate. The addition of the enzyme yielded a yellow coloration with a k max of 414 nm. Spectral changes were not observed in the absence of the enzyme and, therefore, they were considered proof of the activity of AcDODA. The described activity agrees with the absorbance properties reported for betalamic acid and muscaflavin extracted from plants and fungi (Trezzini and Zr€ yb, 1991;Mueller et al., 1997; and for those formed in vitro by the enzyme from the fungus A. muscaria (Girod and Zryd, 1991), from the plants B. vulgaris (Gand ıa-Herrero and Garc ıa-Carmona, 2012) and M. jalapa (Sasaki et al., 2009) and from the bacterium G. diazotrophicus (Contreras-Llano et al., 2019).
The highest activity for the novel DODA enzyme was determined at pH 7.0 (Fig. 4C). Dependence of the a. Activity was determined using 100 µl of protein solution under the assay conditions. b. Crude extract was obtained from a cellular paste harvested from a 0.5 l culture. activity rate on substrate concentration was analysed at the optimum pH determined. In these assays, AcDODA enzyme showed a reduction of its activity at high L-DOPA concentrations, and thus a strong inhibition by excess of L-DOPA as clearly shown in Fig. 4A. This phenomenon is known as inhibition by excess of substrate. In this kinetic model, the plot of the double-reciprocal (inverse activity rate as a function of the inverse of L-DOPA concentration) provides a curve like that shown in Fig. 4B (Segel, 1975). The equation for the obtained curve showed the kinetic parameters estimated as K m = 53 M and V max = 105.8 mM min À1 , and a substrate inhibition constant of 31.9 lM. The kinetic mechanism and its corresponding equation are shown in Fig. S3. The values determined are far from those obtained for previously characterized DODA enzymes, which did not show substrate inhibition for L-DOPA. Previous works reported that values determined for K m were 6.9 mM for B. vulgaris 4,5-DODA (Gand ıa-Herrero and Garc ıa-Carmona, 2012), 7.9 mM for E. coli YgiD (Gand ıa-Herrero and Garc ıa-Carmona, 2014), 3.9 mM for A. Muscaria dioxygenase (Girod and Zryd, 1991) and 1.4 mM for G. diazotrophicus dioxygenase (Contreras-Llano et al., 2019). Thus, the novel enzyme here described presents a singular behaviour with respect to the activity in the presence of L-DOPA as substrate. Figure 3 shows how AcDODA and GdDODA have a close structural similarity with respect to the residues involved in the formation of b-strands and a-helixes. However, the residues that differ between these enzymes could be involved in the kinetic properties of AcDODA. Some of these residues could produce a three-dimensional conformation that affects the union of L-DOPA to the active site. Kinetically the meaning of the inhibition by excess of substrate indicates the possibility of the blockage of the catalytic active site by the access of more than one molecule of substrate. Although the biological significance is unclear, this may function as a relevant regulatory mechanism that has not been described in plants and fungi (Reed et al., 2010).

Functional characterization
Escherichia coli (pRSETA-AcDODA) expressing the novel DODA from A. cylindrica was used to study the evolution of the intermediates of the biosynthetic pathway of betalains and for the biotechnological obtention of their products. The induced expression of the enzyme was performed in LB media under optimal conditions of temperature and IPTG concentration as mentioned above. Then, the reaction was followed for 100 h by the analysis of culture supernatants by HPLC. As Fig. 5 shows, the formation of intermediates 2,3-and 4,5-seco-DOPA occurs in the first hours of reaction and these reached the maximum concentrations after 7 h. They were detected at k = 360 nm with retention times of 6.7 and 7.1 min respectively. Betalamic acid was detected at k = 405 nm with a retention time of 13.8 min and a minor peak was detected with a k max of 403 nm and a retention time of 15.5 min. This peak corresponds to muscaflavin. HPLC analysis also showed a peak with a k max of 470 nm and a retention time of 13.2 min, corresponding to dopaxanthin, which is accumulated due to the condensation of betalamic acid and L-DOPA.
DODA enzyme from A. cylindrica was compared, in terms of sequence homology and phylogeny, with the dioxygenase enzymes previously characterized to produce betalains and with the sequences of potential dioxygenase enzymes related to the G. diazotrophicus sequence. The conserved block of G. diazotrophicus among residues His91 and Asp122 was used for sequence evaluation of the A. cylindrica enzyme (Chang et al., 2015;Li et al., 2015). This block includes one of the three strictly conserved histidines present in the plant enzymes, which is present in GdDODA and in the novel AcDODA as His 101. It also presents acid (Asp100) and hydrophobic (Leu102) residues flanking out the His101, as occurs in the enzymes described from plants (Gand ıa-Herrero and Garc ıa-Carmona, 2014). Phylogenetic analysis (Fig. 6) of the region conserved showed that there is a compact clade formed by plant sequences and the homologue YgiD protein separated from A. muscaria dioxygenase and another compact clade formed by proteobacteria enzymes. However, A. cylindrica appears as an individual clade apart, showing that although its enzyme produces betalamic acid and muscaflavin, AcDODA is not closely related to these two branches.
DOPA-extradiol-dioxygenase, the key enzyme involved in the betalains' biosynthetic pathway of Caryophyllales plants, is also present in the cyanobacterium Anabaena cylindrica. The protein sequence has been optimized for its recombinant expression in E. coli, purified to homogeneity and characterized for the first time as a protein of 17.8 kDa which forms dimers under native conditions. Reactions performed in one-pot experiments show how this enzyme is able to produce betalamic acid and muscaflavin from the amino acid L-DOPA due to a dual activity 2,3-and 4,5-DOPA-extradiol-dioxygenase. This work provides evidence of the presence of DODA enzymes in the phylum cyanobacteria and expands the presence of betalain-forming capabilities in prokaryotes along with the previously described enzymes from E. coli (Gand ıa-Herrero and Garc ıa-Carmona, 2014) and G. diazotrophicus (Contreras-Llano et al., 2019). The sequence analysis of the known enzymes has shown little phylogenetic relationship of the A. cylindrica DODA with those previously described. This seems to indicate that dioxygenase enzymes able to cleave the ring of L-DOPA to produce betalamic acid may be more diverse and frequent than expected. This opens the possibility of finding novel betalamic acid-forming dioxygenases in organisms of different nature with no apparent relationship among them.

Chemicals, bacterial strains, plasmids and enzymes
Distilled water was purified using a Milli-Q system (Millipore, Bedford, MA, USA). HPLC-grade acetonitrile was obtained from Fisher Scientific UK (Leicestershire, UK).

Anabaena cylindrica culture
Anabaena cylindrica was acquired from Spanish Bank of Algae (BEA -Banco Español de Algas, Gran Canaria, Spain). Active growing microorganisms were inoculated in 50 ml cultures of BG-11 medium (Videau and Cozy, 2019) that contains, per litre, 1.5 g NaNO 3 , 0.036 g CaCl 2 -2H 2 O, 0.012 g FeNH 4 -citrate, 0.001 g Na 2 -EDTA, 0.04 g K 2 HPO 4 , 0.075 g MgSO 4 -7H 2 O, 0.02 g Na 2 CO 3 and 1 ml of trace metals solution. Components of trace metal solution are, per litre, 2.86 g H 3 Bo 3 , 1.81 g MnCl 2 -4H 2 O, 0.222 g ZnSO 4 -7H 2 O, 0.391 g Na 2 MoO 4 -2H 2 O, 0.079 g CuSO 4 -5H 2 O and 0.0494 g Co(NO 3 ) 2 . These cultures were maintained seven days at 20°C with orbital shaken at 150 rpm and light/dark cycles of 14:8 h. Afterwards, cultures media were supplemented with varying concentrations of L-DOPA (from 0 to 7.6 mM), sodium ascorbate 15 mM and 100 mM EDTA and further incubated at 20°C. Media were collected after two days and analysed Fig. 6. Phylogenetic analysis of characterized betalamic acid-forming DODAs and related proteobacteria sequences. The analysis was performed with the conserved block of G. diazotrophicus which includes the His 101. This residue is conserved in all characterized 4,5-DOPAdioxygenases known to produce betalamic acid (green dot) and several proteobacteria. The conserved residues were also detected in the novel Anabaena cylindrica enzyme here described to produce dopaxanthin (pink dot), and it conforms an independent clade. Multiple sequence alignment was performed using Clustal Omega, and the phylogenetic tree was inferred by using the phylogeny tool from the ClustalW2 package and plotted using ITOL.

Anabaena cylindrica DODA sequence and cloning
The sequence of the protein WP_015213489 from Anabaena cylindrica, deposited at the National Center for Biotechnology Information (NCBI, Bethesda, MD, USA), was used as a template to synthetically obtain the 4,5-DODA sequence from A. cylindrica. The sequence was enhanced for E. coli expression (Geneart, Regensburg, Germany), and the synthetic DODA gene yielded a 357 bp product that was expressed into the recombinant plasmid pRSETA, which encodes for an additional N-terminal sequence containing a 6 His-tag. The new plasmid pRSETA-AcDODA was transformed into E. coli BL21 (Invitrogen) thermocompetent cells and plated onto LB agar plates containing ampicillin (Amp) 50 lg ml À1 . The resulting colonies were then analysed and confirmed by PCR. PCR amplification was performed using Taq DNA polymerase and the following primers: AcDODA-F (5T GGGGACATATGAAAAGCAAAAC) and AcDODA-R (5Ǵ CTGCAGATCTCGAGTTACAGAC). This amplification yielded a 384 bp product, which coincides with the entire DODA synthetic gene plus the additional 6 His-tag sequence. The plasmid pRSETA-AcDODA was subsequently used in further experiments.

Expression and purification
The AcDODA protein derived from plasmid pRSETA-AcDODA was expressed in E. coli BL21 (Invitrogen), grown at 37°C in LB medium containing Amp 50 lg ml À1 up to an A 600 of 0.8-1.2. Induction was performed for 20 h with different concentrations of isopropyl-1-thio-b-D -galactopyranoside (IPTG) at different temperatures. Protein production for the following steps was performed in 0.5 l cultures of LB medium under the optimized conditions. In the optimal conditions, cells were harvested by centrifugation and resuspended in sodium phosphate buffer 50 mM, pH 8.0, with 0.3 M sodium chloride. Cell lysis was performed by sonication in a Cole-Parmer 4710 series ultrasonic homogenizer (Chicago, IL, USA). Recombinant A. cylindrica DODA was purified from E. coli cells by His-select nickel affinity gel (Sigma) according to the manufacturer's instructions. After non-tagged proteins were washed, the His-tagged protein was eluted with a buffer supplemented with 250 mM imidazole. Immediately after the elution, the protein solution was desalted using PD10 columns (General Electric Healthcare, Milwaukee, USA) and eluted into phosphate buffer 20 mM pH 7.5. Purified protein was quantified using the Bradford protein assay (Bio-Rad, Hercules, CA, USA;Bradford, 1976), and bovine serum albumin was used as standard to obtain a calibration curve. Samples were analysed by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE), by application to 15% polyacrylamide gels and stained using a standard Coomassie Blue method.

Gel filtration
Samples of pure recombinant protein were applied to a Superdex 200 10/300 GL column equilibrated with sodium phosphate buffer 50 mM, pH 7.5, with 150 mM NaCl. The protein was eluted with the same buffer at a flow rate of 0.5 ml min À1 . Elutions were performed in an € Akta purifier apparatus (General Electric Healthcare) and monitored at 280 nm. Column calibration was performed with the following protein markers (Sigma): cytochrome C (12.4 kDa), carbonic anhydrase (29 kDa), albumin (66 kDa), alcohol dehydrogenase (150 kDa) and b-amylase (200 kDa).

MALDI-TOF-MS protein analysis
Matrix solution for peptide analyses was a-cyano-4-hydroxycinnamic acid (20 mg ml À1 ) in ACN/water/TFA (70:30:0.1). The peptide sample was dissolved in 0.1% TFA and mixed with the matrix solution. One microlitre of this mixture was applied to the atmospheric pressure matrix-assisted laser desorption ionization (AP-MALDI) target plate and allowed to dry. Experiments were carried out with an Agilent Time of Flight (TOF) Mass Spectrometer (Agilent Technologies, Santa Clara, CA, USA), equipped with an AP-MALDI Ion Source with an N 2 laser (337 nm). Samples were measured in reflectron mode to identify molecular formulas based on precise mass measurements in positive mode. External calibration of the spectrometer was performed with standard peptides from the ProteoMass TM Peptide MALDI-MS Calibration Kit (Sigma). Data were recorded and processed with Agilent MassHunter Workstation Software. Peptide mass fingerprint determination was carried out using Agilent Spectrum Mill Software. Determination of protein absolute molecular mass was carried out using an HPLC-ESI-MS-TOF system. This system comprises an HPLC Agilent VL 1100 apparatus equipped with an autosampler l-well plate and a capillary pump, connected to an Agilent 6100 Mass Spectrometer Time of Flight (MS-TOF), and an electrospray ionization interface (ESI) was used. The column employed was a Zorbax Poroshell 300SB-C18, 1 9 75 mm with 5 lm (Agilent Technologies, Santa Clara, CA, USA). The column was operated at 60°C, and the samples were injected with a flux of 0.2 ml min À1 . The protein was eluted with a linear gradient using water/ACN/formic acid (95:4.9:0.1) as solvent A, and water/ACN/formic acid (10:89.9:0.1) as solvent B.
A linear gradient was performed for 30 min from 0% to 90% of B. Protein separation was monitored at 210 and 280 nm using a multiple wavelength detector. The mass spectrometer was operated in positive mode in the range of 100-2200 m/z, using a capillary voltage of 3.5 kV. Nebulizer gas pressure was 30 psi, and drying gas flux was 8 l min À1 at a temperature of 350°C. External spectrometer calibration was carried out using the Pro-teoMass TM Peptide MALDI-MS Calibration Kit (Sigma-Aldrich, St. Louis, MO, USA). Two different peptides were used as controls (cytochrome C and carbonic anhydrase, Sigma-Aldrich). All data were recorded and processed through the software Agilent MassHunter Workstation Qualitative Analysis Software (Agilent Technologies), and the intact molecular weight of the protein was obtained using the deconvolution algorithm from this software.

Trypsin digestion
The protein sample was prepared in 100 ll of buffer NH 4 HCO 3 50 mM, pH 8.0, with 0.02% ProteaseMAX TM Surfactant (Promega, Madison, WI, USA). Then, the sample was reduced with DTT 10 mM for 20 min at 56°C and alkylated with iodoacetamide 50 mM at room temperature in the dark for 20 min. One microgram of proteomics grade trypsin (Promega) was added and the sample was incubated for 4 h at 37°C. After that, the sample was centrifuged at 15 000 g for 1 min to collect the condensate and the digestion was stopped by adding 0.5% TFA. Peptides were cleaned up with C18 Zip-Tips (Millipore) and evaporated using an Eppendorf vacuum concentrator model 5301.

HPLC analysis
A Shimadzu LC-20AD apparatus (Kyoto, Japan) equipped with an SPD-M20A photodiode array detector was used for analytical HPLC separations performed with a 250 9 4.6 mm Kromasil 100 C-18 column packed with 5 lm particles (Teknokroma, Barcelona, Spain) (Gand ıa-Herrero et al., 2005). A linear gradient was performed using water with 0.05% TFA as solvent A, and acetonitrile with 0.05% TFA as solvent B. The linear gradient was performed for 24 min from 0% B to 35% B with the flow was 1 ml min À1 , and the injection volume was 50 ll.

Absorbance spectroscopy
Enzyme activity was determined using a continuous spectrophotometric method by measuring the absorbance at k = 414 nm due to the presence of betalamic acid and muscaflavin (Girod and Zryd, 1991; Gand ıa-Herrero and Garc ıa-Carmona, 2012). Unless otherwise stated, the reaction medium contained sodium phosphate buffer 50 mM, pH 7, L-DOPA 7.6 mM, and sodium ascorbate 100 mM at a final volume of 300 ll. Measurements were performed at 25°C in 96-well plates in a Synergy HT plate reader (Bio-Tek Instruments, Winooski, USA). Betalamic acid solutions of known concentration were employed to calibrate the plate reader detector signal. The molar extinction coefficient of betalamic acid at 424 nm, Ɛ = 24 000 M -1 cm -1 , was employed for the quantification of betalamic acid (Trezzini and Zr€ yb, 1991). Kinetic data analysis was performed by using non-linear regression fitting with the SigmaPlot Scientific Graphing for Windows version 10.0 (2006, Systat Software, San Jose, CA, USA).
E. coli (pRSETA-AcDODA) culture conditions Escherichia coli BL21 cells harbouring pRSETA-AcDODA were grown in 50 ml LB supplemented with Amp 50 lg ml À1 until reached an O.D. 600 0.8-1.0. Then, heterologous expression of the DODA enzyme was performed by addition of 0.1 mM IPTG. The culture was kept at 20°C for 20 h and then centrifuged 10 min at 7 500 rpm. Supernatant was discarded, and cells were diluted in 50 ml of distilled water containing L-DOPA 7.6 mM and sodium ascorbate 15 mM. Culture was maintained at 20°C with orbital shaken for 200 h, and samples were taken at different hours in order to follow the formation of derived products. Samples were centrifuged 1 min at 7 500 rpm, and the supernatant was immediately analysed by HPLC.

Phylogenetic analysis of novel betalamic acid-forming DODAs
The DODA enzyme from G. diazotrophicus was used as a template to search similar enzymes. The obtention of a neighbour-joining phylogenetic tree (Altschul et al., 2005;Li et al., 2015) found the DODA enzyme from A. cylindrica by sequence similarity. The relationship of the novel enzymes with betalamic acid-forming dioxygenases previously reported was performed from the conserved block of G. diazotrophicus which includes the His 101. In all cases, multiple sequence alignment was performed using Clustal Omega from the EMBL-EBI (European Molecular Biology Laboratory -European Bioinformatics Institute) (Madeira et al., 2019) and the phylogenetic tree was inferred by using the phylogeny tool from the ClustalW2 package of the same platform and plotted using iTOL (Interactive Tree of Life; Letunic and Bork, 2019). assistance in ESI-MS/MS, MALDI-TOF and TOF/Q-TOF mass spectrometry experiments. This work was supported by 'Ministerio de Ciencia e Innovaci on', project AGL2017-86526-P (MCI/AEI/FEDER, UE) and by 'Programa de Ayudas a Grupos de Excelencia de la Regi on de Murcia, Fundaci on S eneca, Agencia de Ciencia y Tecnolog ıa de la Regi on de Murcia' (Project 19893/ GERM/15). M.A.G.-R. holds a contract financed by MEC-FSE (Spain). Fig. S1. Phylogenetic analysis of dioxygenase enzymes related to Gluconacetobacter diazotrophicus (orange). These enzymes correspond to the phyllum proteobacteria (yellow) and cyanobacteria (blue), including Anabaena cylindrica (navy blue). Fig. S2. Determination of the optimal conditions for the heterologous expression of AcDODA in E. coli cells. A: Comparison of the obtention of the products from the activity of AcDODA and derived compounds in E. coli cultures supplemented with different concentrations of IPTG and maintained under different temperatures. The reactions were followed by HPLC in order to evaluate the presence of the enzyme. B: Macroscopic images of E. coli (pRSETA-AcDODA) cultures, after 24 hours of reaction, supplemented with 1.0 (left), 0.5 (centre) or 0.1 (right) mM IPTG and maintained at 20°C (top) or 25°C (bottom).