Mononuclear binding and catalytic activity of europium(III) and gadolinium(III) at the active site of the model metalloenzyme phosphotriesterase

The crystal structures of europium(III)- and gadolinium(III)-bound phosphotriesterase are presented along with esterase activity data for the lanthanide-bound enzymes.

The 4f valence shell of Ln(III) ions is shielded by the outer 5s and 5p shells, resulting in a 4f shell that lies inert, deep in the cation (Bu ¨nzli, 2006;Mikami et al., 2002;Evans, 2000).This buried 4f orbital offers poor shielding of the outer orbitals from nuclear charge, resulting in the relatively small size of the lanthanides, referred to as the 'lanthanide contraction' (Cotruvo, 2019).Another consequence of the buried 4f orbitals is a lack of orbital overlap with coordinating ligands with respect to d-block transition-metal coordination, resulting in organolanthanide interactions that are predominantly ionic in nature (Cotruvo, 2019;Mikami et al., 2002).These ionic interactions result in little steric requirement for coordination, in stark contrast to d-block transition metals; thus, Ln(III) ions have a preference for coordination numbers (CNs) of 8-12 (Jahn et al., 2018;Deng et al., 2018;Pol et al., 2014;Bu ¨nzli, 2006;Franklin, 2001;Shannon, 1976).A further consequence of the buried and diffuse 4f orbitals, and the predominantly ionic interactions, is that Ln(III) ions have little ligand-field splitting effect; thus, Ln(III) complexes have fast ligandexchange kinetics (Bu ¨nzli, 2006;Franklin, 2001).

Lanthanides as catalysts
Ln(III) chemistry, specifically the high Lewis acidity, high charge, fast ligand-exchange kinetics, high CN, flexible coordination geometry and redox inertness, make Ln(III) ions ideal for catalysis (Cotruvo, 2019;Lim & Franklin, 2004).For example, the NPAC-La(III) complex is capable of hydrolysing the phosphodiester bonds of DNA and RNA (Baykal & Akkaya, 1998).This complex uses neutral ligands so as to not reduce the Lewis acidity of La(III) (Baykal & Akkaya, 1998).The Lewis acidity and high CN of the metal are required for catalysis to occur (Franklin, 2001;Komiyama et al., 1999).
With their unique chemistry, it may be that Ln(III) ions could be used to augment enzyme chemistry, permitting the turnover of synthetic non-natural substrates, whilst retaining the high regioselectivity and stereoselectivity of enzymes.Indeed, haem-containing metalloenzymes have been engineered via the introduction of non-natural metals in place of iron, which has resulted in non-natural catalysis of synthetic substrates (Dydio et al., 2017;Key et al., 2016Key et al., , 2017)).

Lanthanides in biology
Previous proposals have suggested that lanthanoenzymes should be more common than they are, and that Ln(III) is better suited for catalysis than Ca(II), for instance (Lim & Franklin, 2004).Despite this, there are few known examples of lanthanoenzymes and lanthanoproteins.
The first natural lanthanoenzyme to be described was XoxF-MDH in 2011, a methyl dehydrogenase (MDH) containing a pyrroloquinoline quinone (PQQ) cofactor along with a Ln(III) cofactor, which catalyses the two-electron oxidation of methanol (Good et al., 2020;Vu et al., 2016;Nakagawa et al., 2012;Fitriyanto et al., 2011;Hibi et al., 2011).There is also a Ca(II)-binding MDH analogue, MxaF-MDH, which catalyses the same oxidation of methanol (Chistoserdova & Lidstrom, 1997;Anthony & Zatman, 1964a,b).The two enzymes are very similar, differing only in an additional aspartate in the Ln(III) coordination sphere of XoxF-MDH, increasing the CN from seven to nine, which satisfies the high CN preference and stabilizes the 3+ oxidation state (Good et al., 2020).As a result of the extra ligand, the Ln(III) ion has a higher binding affinity for XoxF-MDH than Ca(II) does for MxaF-MDH (Nakagawa et al., 2012).The high charge density and Lewis acidity of Ln(III) with respect to Ca(II) results in XoxF-MDH having a higher catalytic efficiency than MxaF-MDH (Nakagawa et al., 2012;Fitriyanto et al., 2011;Hibi et al., 2011).This may suggest that Ln(III) ions are indeed preferable to Ca(II) in enzymes, as was initially suggested in 2004, and indicates that the development of lanthanoenzymes would be a fruitful line of enquiry that could produce catalytically efficient enzymes (Lim & Franklin, 2004).Since this discovery, further Ln(III)bound alcohol oxidases have been described such as ExaF, an ethanol dehydrogenase (Good et al., 2016).
In 2018, an Ln(III)-binding protein was discovered called lanmodulin (LanM), which binds Ln(III) ions using EF-hand motifs (Cook et al., 2019;Cotruvo et al., 2018).The function of this protein has yet to be determined, but clearly indicates the potential for a greater contribution of Ln(III) ions to protein fold and function in the bacterial kingdom (Featherston & Cotruvo, 2021;Cotruvo et al., 2018).Lanmodulin provides lessons in designing Ln(III) chelators, as the LanM EF-hand has an additional aspartate with respect to calmodulin EF-hands, thus satisfying the high CN preference of Ln(III) ions and the high charge, but also indicating the similarity in binding domains between the metals (Cook et al., 2019;Cotruvo, 2019;Cotruvo et al., 2018).Notably, this difference provides lanmodulin with higher binding affinities for Ln(III) than calmodulin has for Ca(II); these are in the micromolar range (VanScyoc et al., 2002;Linse et al., 1991).
To date there has been one example of an engineered lanthanoenzyme, the chimeric 33-mer metallopeptide 'P3W' formed by the fusion of an EF-hand motif of calmodulin and the helix-turn-helix motif of the engrailed homeodomain (Sirish & Franklin, 2002;Kim et al., 2001;Welch et al., 2001).The design of this chimera was achieved by superposition of the structures of calmodulin and the engrailed homeodomain, which have very similar structures.The chimera permitted DNA binding and hydrolysis of supercoiled DNA, due to the proximal positioning of the Ln(III) ion to the phosphodiester bond.Overall, these enzymes provide proof-of-principle that lanthanoenzymes can be engineered and that Ln(III) ions maintain sufficient Lewis acidity for catalysis when coordinated by negatively charged ligands.
Most work coordinating Ln(III) ions to proteins has been for the purpose of X-ray crystallography or NMR, and has been executed using synthetic organic polydentate ligands that are attached via non-natural amino acids or cysteine residues (Herath et al., 2021;Loh et al., 2013;Allen & Imperiali, 2010;Silvaggi et al., 2007).Ln(III) ions scatter X-rays extremely well and thus can aid in solving the phase problem when solving crystal structures (Silvaggi et al., 2007;Harker, 1956).Ln(III) ions are also useful for NMR studies due to their paramagnetism, which produces pseudocontact shifts in the protein of up to 40 A ˚from the ion; thus, tagging a protein with an Ln(III) ion provides long-range structural information (Pilla et al., 2017;Yagi et al., 2013;Schmitz et al., 2012;Allen & Imperiali, 2010).However, for crystallography the Ln(III) ion must be well ordered in relation to the protein, and for NMR Ln(III) mobility results in decreased information; therefore, a high-affinity Ln(III)-binding tag is required (Allen & Imperiali, 2010).Hence, the Ln(III)-binding tags used in protein NMR experiments have the requisite features to accommodate Ln(III) ions: a high CN and hard Lewis base ligands (Herath et al., 2021;Loh et al., 2013).
Development of lanthanoenzymes could result in biocatalysts that are capable of catalysing non-natural reactions using non-natural substrates, resulting in potentially useful products being produced at low temperatures under gentle conditions.In this work, we describe a lanthanide-bound model enzyme.Phosphotriesterase (PTE), a binuclear Zn(II) metalloenzyme from Pseudomonas diminuta, has previously undergone directed evolution to change its activity from a phosphotriesterase to an esterase of 2-naphthyl hexanoate (2NH); the round 18 variant (PTE-R18) in this engineering pathway had an additional Zn(II)-binding site within the active site, indicating a loss of specificity for two divalent metals (Tokuriki et al., 2012).PTE-R18 was generated from PTE by the mutations H254R, D233E, F306I, I274S, T172I, S269T, M138I, T199I, L272M, A80V, S111R, A204G, L130V, L271F, A49V, K77E, L140M and I313F (Tokuriki et al., 2012).Given that PTE is naturally promiscuous for divalent metal coordination, combined with this additional loss of specificity for binding two cations, it was hypothesized that PTE-R18 could have the potential to bind unusual metals in its active site.Ln(III) ions were an interesting candidate, as their larger size compared with Zn(II) ions could result in mononuclear occupancy of the unusual 'third' binding site in PTE-R18.Here, we describe the coordination and affinity of these mononuclear Ln(III)-PTE-R18 complexes using X-ray crystallography and isothermal titration calorimetry, and we show that esterase activity is maintained by the enzyme.

Production and purification of PTE-R18
PTE-R18 was overexpressed from pETMCSI-PTE-R18 (see the supporting information).Production and purification were performed as described previously (Campbell et al., 2016).

Formation and verification of apo PTE-R18
Zn(II) was removed from PTE-R18 by two rounds of dialysis against chelation buffer (5 mM HEPES, 100 mM NaCl, 5 mM 1,10-phenanthroline pH 8.0) in a 1:100 ratio through SnakeSkin Dialysis Tubing 10K MWCO 35 mm (Thermo Fisher Scientific, USA) for 24 h at 4 � C. The chelator 1,10-phenanthroline was subsequently removed by four exchanges against dialysis buffer (5 mM HEPES, 100 mM NaCl pH 8.0) at a 1:100 ratio through SnakeSkin Dialysis Tubing 10K MWCO 35 mm for 8 h at 4 � C. PTE-R18 retention was determined by SDS-PAGE.The formation of apo PTE-R18 was confirmed via incubation of 100 nM enzyme with 2.5 mM p-nitrophenyl butyrate in assay buffer [20 mM Tris, 100 mM NaCl, 2.5%(v/v) methanol pH 8.5] with the positive control containing 100 mM ZnCl 2 , and analysis of the absorbance at 405 nm over 20 min at room temperature.

Crystallization
All proteins were buffer-exchanged into 20 mM HEPES, 50 mM NaCl pH 8.0 via centrifugation at 3000g and 4 � C using an Amicon Ultra-15 Centrifugal Filter 10 kDa MWCO prior to crystallization.All crystals were grown at 4 � C using 100 mM sodium cacodylate, 14%(v/v) MPD pH 6.5 as mother liquor (ML) by hanging-drop diffusion.Zn-PTE-R18 (22 mg ml À 1 ) was crystallized at a 2:1 ML:protein ratio.These crystals were manually pulverized for use as a Zn seed stock.

Fluorescence scans
To assess the presence of lanthanides in the crystal samples, the fluorescence emission spectrum was measured using Se Kedge detection with input energy 12 857.8eV and a scan time of 30 s using BL_3ID1 on the MX2 beamline at the Australian Synchrotron (Araga ˜o et al., 2018).

Structure determination
Diffraction data were collected on beamline MX2 at the Australian Synchrotron (Araga ˜o et al., 2018).Data for the apo PTE-R18 structures were collected at 13 and 9.5 keV, and data for the Ln-PTE-R18 structures were collected at 9.5 keV to allow the collection of anomalous data resulting from the lanthanides.All data were collected with 70% attenuation of the beam, with detector distances of 130 mm for the apo and Eu(III)-bound structures and 150 mm for the Gd(III)-bound structure.

Isothermal titration calorimetry
To determine the binding stoichiometries of Eu(III) and Gd(III) to apo PTE-R18, isothermal titration calorimetry was performed using a Nano ITC low-volume calorimeter (TA Instruments).All experiments were performed at 20 � C with stirring at 200 rev min À 1 , with 23 � 2 ml serial injections of 1.5 mM EuCl 3 and 2 mM GdCl 3 into 238 and 360 mM apo PTE-R18 (monomer concentration), respectively, with 360 s intervals.
Data were integrated using NITPIC version 1.2, followed by global analysis using SEDPHAT version 12.1 (Scheuermann & Brautigam, 2015;Zhao et al., 2015).Data were analysed using the 'A + B $ AB Hetero-Association' model and the Marquardt-Levenberg equation to fit a function to the data.The data were plotted and visualized using GUSSI.

Ln(III)-PTE-R18 activity assays
To determine the activities of Eu(III)-and Gd(III)-bound PTE-R18, hydrolysis of 2NH was quantified.Assays consisted of 500 mM EuCl 3 or GdCl 3 , 50 mM apo PTE-R18, 10 mM 2NH substrate and 5%(v/v) methanol and were performed in 20 mM Tris, 150 mM NaCl pH 7.0 or 8.0.Reactions were incubated at room temperature with mixing for 42 h.Control reactions were performed without apo PTE-R18.Reactions were treated with a 4.5� volume of acetonitrile with 2.2 mM 1,3,5-trimethoxybenzene, followed by filtration through a 0.22 mm filter and concentration under nitrogen gas.The formation of 2-naphthol (2N) was determined by 1 H NMR at 400 MHz using CDCl 3 .

Formation of apo PTE-R18
The strong chelating agent 1,10-phenanthroline was used to remove the co-purified and tightly bound Zn(II) ions to permit lanthanide binding in the active site of PTE-R18.The removal of Zn(II) was confirmed using a colorimetric assay and X-ray diffraction experiments on apo PTE-R18 crystals.
3.1.1.Hydrolysis assay.To test for the removal of Zn(II), hydrolytic enzyme activity against para-nitrophenyl butyrate (pNPB) was assessed.Previous work has demonstrated the ability of PTE-R18 to rapidly hydrolyse para-nitrophenyl acetate, yielding a yellow product (Tokuriki et al., 2012).The product of the hydrolysis reaction, para-nitrophenol (paranitrophenolate at pH 8.5), was measured by absorbance at 405 nm (Fig. 1).The dialysed enzyme showed a complete lack of para-nitrophenol production, with an absorbance at 405 nm similar to those of the buffer controls, whereas the Zn(II)supplemented apo PTE-R18 sample demonstrated a rapid increase in absorbance at 405 nm, indicating a loss of activity of the dialysed protein and thus successful Zn(II) removal.The reconstituted enzyme achieved complete hydrolysis of 2.5 mM pNPB after �30 s, resulting in an estimated k cat of 809 � 165 s À 1 , which is comparable to the reported k cat for hydrolysis of the ester 2NH of 880 � 30 s À 1 , demonstrating that the dialysed protein can be reconstituted with metal ions to give native holoenzyme activity (Kaltenbach et al., 2015).Previous work has also shown that 1,10-phenanthroline is capable of chelating Zn(II) from PTE-R18 and that PTE-R18 is stable in the apo form (Tokuriki et al., 2012).
3.1.2.Apo PTE-R18 crystallography.X-ray diffraction data were collected from apo PTE-R18 crystals at 13 and 9.5 keV to maximize the anomalous scattering from Zn(II) and Ln(III) ions, respectively (Table 1).This was performed to confirm the absence of the metal ions from the active site of apo PTE-R18.Both the 13 keV (1.50A ˚resolution; PDB entry 8uqw) and 9.5 keV (1.52A ˚resolution; PDB entry 8uqx) structures show a lack of anomalous signal at the expected Zn(II) positions in the active site, indicating the successful removal of Zn(II) and the lack of any anomalous scatterer which could be confused for an Ln(III) signal (Figs.2a and 2c).However, in the 9.5 keV structure an anomalous signal visible to 6� was present in the active site; this signal was not present in the 13 keV structure above noise.The anomalous signal was suggested to be a result of the Cl K edge (2.8 keV) and thus was modelled as a chloride in both structures given the proximity to Lys169 and the abundance of NaCl in the sample from the purification buffers (Merritt, 2012).
In both PDB entries 8uqw and 8uqx the active-site residue Lys169 was modelled in two separate conformations (Figs.2a  and 2c); normally Lys169 undergoes spontaneous carboxylation upon metal binding and directly coordinates to the activesite Zn(II) ions (Tokuriki et al., 2012).However, in the apo structure Lys169 is decarboxylated and thus is capable of adopting two conformations.
In both structures residues 260-275 of chain A were unable to be modelled due to a lack of electron density, whereas in chain B these residues were modelled into clear density.This loop (loop 7) inherently undergoes conformational fluctuations during catalysis (Jackson et al., 2009;Caldwell et al., 1991).Previous directed evolution of wild-type PTE into PTE-R18 resulted in stabilization of loop 7, and loop 7 was able to be built into both chains A and B of Zn-PTE-R18 crystal structures (Campbell et al., 2016;Tokuriki et al., 2012).This suggests that the removal of Zn(II) resulted in increased mobility of loop 7, which was observed in chain A of the crystal structure, but loop 7 of chain B was rigidified by stabilizing crystal contacts and was able to be modelled.
Previous structures of wild-type PTE demonstrated that the apo PTE crystal structure was substantially different from the crystal structure of the holoenzyme (Benning et al., 1994(Benning et al., , 1995)).In contrast, our apo PTE-R18 structure is very similar to the holoenzyme crystal structures, with r.m.s.d.values of 0.184 and 0.191 A ˚for PDB entries 8uqw and 8uqx, respectively, against PDB entry 4e3t (Tokuriki et al., 2012).This suggests that the enzyme structure was stabilized during the course of laboratory directed evolution performed previously, Para-nitrophenyl butyrate hydrolysis assay.Para-nitrophenol absorbance at 405 nm, measured over 5 min, with buffer consisting of 20 mM Tris, 100 mM NaCl, 2.5%(v/v) methanol pH 8.5 (magenta) and with buffer with 100 mM ZnCl 2 (black), with 100 nM dialysed protein (cyan) and with 100 nM dialysed protein and 100 mM ZnCl 2 (green) (N = 3).
resulting in an apoenzyme crystal structure that is similar to that of the holoenzyme.

Ln(III)-PTE-R18 crystallography.
A range of Ln(III) ions were chosen for co-crystallization to provide a variation in radius and Lewis acidity, with La(III), Pr(III), Sm(III), Eu(III) and Gd(III) being co-crystallized with apo PTE-R18.Diffraction data were collected at 9.5 keV from the resulting crystals to maximize the anomalous scattering from any bound Ln(III) ions.The structures of crystals grown with La(III), Pr(III) and Sm(III) were solved to 1.50, 2.16 and 1.84 A resolution, respectively, and showed a lack of electron density and anomalous signal within the active site, indicating that the desired Ln(III) was not bound within the active site of PTE-R18.Both the Eu(III)-bound (1.78A ˚resolution; PDB entry 8uqy) and Gd(III)-bound (1.61A ˚resolution; PDB entry 8uqz) structures showed a strong single spherical anomalous signal within the active site, indicating mononuclear binding of Ln(III) to PTE-R18 (Figs. 3a and 3c).The greater ionic radius of Eu(III) and Gd(III) compared with Zn(II) would explain the observation of only one ion being accommodated within the active site.These anomalous signals were present to over 12� and 40� for Eu(III) and Gd(III), respectively, supporting their placement in the structure.However, these anomalous signals were only present in chain A of the PTE-R18 dimer in both structures.Metals were modelled into this density and the occupancies were refined to 16% and 24% for Eu(III) and Gd(III), respectively.These low occupancy values and the presence of the metal in chain A only could be the result of weak binding to apo PTE-R18 or due to Ln(III) coordination by abundant solvent molecules, such as the high-concentration cacodylate or MPD, which would sequester Ln(III) from PTE-R18.Furthermore, the presence of the Ln(III) ions in only chain A of the dimer could be the result of altered conformations or dynamics between the chains; in chain A of the Gd-PTE-R18 structure residues 260-274 (260-275 in Eu-PTE-R18) were unable to be modelled due to a lack of electron density, whereas in chain B these residues were modelled; this was also the case for the apo PTE-R18 structures.
We expected both chains A and B to have bound Ln(III).However, it may be that the binding of Ln(III) to PTE-R18 does not stabilize loop 7 as Zn(II) does, which would explain the missing electron density for loop 7 of chain A in the Ln(III)-bound structures.The lack of Ln(III) in chain B could be explained as follows: when loop 7 of chain B adopts the rigid conformation that is required for crystallization, the monomer forms a conformation that is incapable of coordinating Ln(III), potentially ejecting any large Ln(III) ion from the active site and resulting in chain B being metal-free in the crystal.Alternatively, the binding of an Ln(III) ion to chain A of PTE-R18 could result in slight conformational changes to chain B that make subsequent Ln(III) binding unfavourable, resulting in the binding of Ln(III) ions to only one chain of the dimer.
Nonetheless, the fluorescence scans displayed X-ray emission lines at the expected values for Eu(III) and Gd(III) (Figs. 4b and 4d), further indicating that Eu(III) and Gd(III) are the coordinating metals (Bearden, 1967;Kortright & Thompson, 2009).
The Eu(III)-and Gd(III)-bound structures show identical coordination within the PTE-R18 active site, being coordinated by His55, His57, Asp301, H 2 O (252 in PDB entry 8uqy and 158 in PDB entry 8uqz) and a Cl À bridging from Lys169 (Figs. 3a and 3c).The positioning of the Ln(III) ion is nearidentical to the �-Zn(II) ion in previous PTE-R18 structures (PDB entries 4e3t, 4gy0 and 4gy1) with very similar coordination, varying by the Ln(III) ions having two coordinate bonds from Asp301 rather than one and by the Ln(III) ions having a bridging Cl À rather than a carboxylated Lys169 (Tokuriki et al., 2012).This was rather unexpected, as structures of Zn-PTE-R18 show that Lys169 is carboxylated, providing two hard coordinating O atoms; however, this is not the case in the Eu(III)-and Gd(III)-bound structures, despite these metals favouring hard Lewis bases and high coordination numbers.Also, the bridging Cl À would be less entropically favourable when compared with direct amino-acid coordination, but Cl À was present in the apo structures, suggesting that the entropic cost to binding would not be substantial with regard to the Cl À ion.
Interestingly, both Ln(III) ions are coordinated to histidine residues, which has not previously been observed in protein structures, and have a CN of six, with five coordinate bonds from the protein.This low CN, combined with coordination from slightly softer histidines, could be the reason for the low metal occupancy in the structures.Typically, lanthanide-bound protein structures have CNs of 9-10, although some have been reported to be as low as seven, and these coordination groups are hard Lewis bases, often oxygen.Indeed, a de novo lanthanide-binding protein has previously been designed by the Zeymer group with four coordinating Glu residues and one water, and an engineered chimeric nuclease from the Franklin group was coordinated using one Glu, three Asp and one Thr residues (Caldwell et al., 2020;Sirish & Franklin, 2002;Kim et al., 2001;Welch et al., 2001).Furthermore, the natural lanthanide-binding proteins XoxF and LanM have been described to coordinate using two Asp, one Glu and one Asn residues and three groups of the cofactor PQQ, and one Glu, four Asp and one Thr residues, respectively (Cook et al., 2019;Cotruvo, 2019;Cotruvo et al., 2018;Pol et al., 2014).
These structures illustrate the ability of PTE-R18 to bind Eu(III) and Gd(III) in the active site of one chain of the dimeric enzyme and show that histidines can be ligands for the coordination of metal ions in lanthanoproteins.

Isothermal titration calorimetry.
To further substantiate mononuclear binding, and to provide bindingaffinity data, isothermal titration calorimetry (ITC) was performed using apo PTE-R18 with EuCl 3 and GdCl 3 .Both the Eu(III) (N ' 1.1) and the Gd(III) (N ' 1.1) titrations suggested mononuclear binding of the metal to each monomer of the apo PTE-R18 dimer (Figs.4a and 4b).This suggests that Eu(III) and Gd(III) do bind to both chains A and B of PTE-R18 and that the lack of metal in chain B is likely to be the result of crystallization.The isotherms also indicated K d values of 11 mM (95% confidence interval 7.2-17.2mM) and 20 mM (95% confidence interval 14.2-26.8mM) for Eu(III) and Gd(III), respectively.Based upon these dissociation constants and the use of 2.77 mM Ln(III) during the crystallization experiments (100 times the K d ), we could expect occupancy values of 90%, rather than 16% and 24%.However, these dissociation constants are based upon solution chemistry and differences would be expected as a result of the crystal contacts.
These K d values are relatively high when compared with other lanthanide-binding proteins, which could be a result of the low CN and the relatively weak ligands.For example, Zeymer's de novo lanthanide-binding protein, with numerous hard ligands, was described as having subfemtomolar affinity (Caldwell et al., 2020).Moreover, the equally numerous and hard ligands of LanM resulted in picomolar affinity values for lanthanides (Cook et al., 2019;Cotruvo, 2019;Cotruvo et al.,

2018
).Therefore, it could be concluded that the weak binding of the metals was due to the low coordination number and the softness of the histidine ligands.Indeed, EF4 of LanM has an Asn ligand in one position where EF1, EF2 and EF3 have Asp, and this reduces the affinity to micromolar values (Cotruvo et al., 2018).
This, however, does not rule out Eu-PTE-R18 and Gd-PTE-R18 as potential catalysts, as the Franklin chimeric nucleases were determined to have affinity values for Eu(III) of 3-20 mM whilst being catalytically active (Sirish & Franklin, 2002;Kim et al., 2001).Also, the high-affinity Zeymer protein and LanM completely saturate the primary coordination sphere of the metal, leaving no space for a substrate to coordinate; thus, it may not be reasonable to expect such high affinities from a substrate-accessible lanthanoenzyme.Furthermore, these high-affinity proteins use negatively charged ligands, which would reduce the Lewis acidity of the metal, thus reducing a catalytically favourable characteristic of Ln(III).

Esterase activity of Ln(III)-PTE-R18
The coordination of a catalytically relevant metal ion into the active site of PTE-R18 led us to hypothesize that there should be catalytic activity.Given that apo PTE-R18 showed no esterase activity against the more reactive pNPB, we tested for esterase activity against 2NH (Table 2), which PTE-R18 has been documented to hydrolyse (Tokuriki et al., 2012).
This kinetic analysis indicated that aqueous Eu(III) and Gd(III) had no detectable esterase activity against 2NH under these conditions, whereas Eu(III)-and Gd(III)-bound PTE-R18 displayed esterase activity against 2NH, as determined by the detection of 2N and hexanoic acid.Both Eu(III)-PTE-R18 and Gd(III)-PTE-R18 showed increased conversion at pH 8.0 compared with pH 7.0, which is likely to be a result of Eu(III) and Gd(III) having pK a values of 7.66 and 7.87, respectively, resulting in increased deprotonation of a metal-ion-coordinated water at pH 8.0 to form a stronger nucleophile for attack at the carbonyl of 2NH, consistent with the proposed mechanism of PTE-R18 (Brown & Ekberg, 2016;Tokuriki et al., 2012).Thus, despite coordinating a single large metal cation, PTE-R18 was capable of binding 2NH within the active site and the Ln(III) holoenzyme was capable of catalysing ester-bond hydrolysis.

Conclusion
To conclude, this work demonstrated the ability of PTE-R18 to coordinate a single Eu(III) or Gd(III) ion within the active site and described the first observation of Ln(III) ions being coordinated by histidine residues in proteins.The ability of this enzyme to coordinate catalytically important lanthanide metals into the accessible active site suggests potential for lanthanoenzyme-based catalysis, as was demonstrated by the esterase activity against 2NH.Furthermore, this work provides the opportunity for mutagenesis experiments to increase the CN and Lewis base hardness of coordinating groups, which could result in an altered affinity of PTE-R18 for Eu(III) and Gd(III).

Figure 2
Figure 2 Crystallographic study of apo PTE-R18.(a) Structure of the apo PTE-R18 chain A active site from data collected at 13 keV (PDB entry 8uqw).(b) Fluorescence emission of the apo PTE-R18 crystal used to collect diffraction data at 13 keV.(c) Structure of the apo PTE-R18 chain A active site from data collected at 9.5 keV (PDB entry 8uqx) with the anomalous map at 3� (cyan) and 5� (magenta).(d) Fluorescence emission of the apo PTE-R18 crystal used to collect diffraction data at 9.5 keV.

Figure 4 ITC
Figure 4 ITC analysis of Ln(III) binding to PTE.ITC data for metal binding to apo PTE-R18 processed using the A + B $ AB binding model (N = 1).(a) Isotherm for EuCl 3 .(b) Isotherm for GdCl 3 .

Table 1
Crystallographic statistics.Data-collection, processing and refinement statistics for the crystal structures deposited as PDB entries 8uqw, 8uqx, 8uqy and 8uqz.Values in parentheses are for the highest resolution shell.