Structural and Biochemical Characterization of an
Active Arylamine N-Acetyltransferase Possessing a
Non-canonical Cys-His-Glu Catalytic Triad

Xavier Kubiak, Inès Li de La Sierra-Gallay, Alain F. Chaffotte, Benjamin Pluvinage, Patrick Weber, Ahmed Haouz, Jean-Marie Dupret, Fernando Rodrigues-Lima

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

To cite this version:

Xavier Kubiak, Inès Li de La Sierra-Gallay, Alain F. Chaffotte, Benjamin Pluvinage, Patrick Weber, et al.. Structural and Biochemical Characterization of an Active Arylamine N-Acetyltransferase Possessing a Non-canonical Cys-His-Glu Catalytic Triad. Journal of Biological Chemistry, 2013, 288 (31), pp.22493 - 22505. 10.1074/jbc.m113.468595. pasteur-03525397

HAL Id: pasteur-03525397

https://pasteur.hal.science/pasteur-03525397

Submitted on 13 Jan 2022

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

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

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

Distributed under a Creative Commons Attribution 4.0 International License
# Structural and Biochemical Characterization of an Active Arylamine N-Acetyltransferase Possessing a Non-canonical Cys-His-Glu Catalytic Triad*§

Received for publication, March 11, 2013, and in revised form, June 13, 2013 Published, JBC Papers in Press, June 16, 2013, DOI 10.1074/jbc.M113.468595

Xavier Kubiak<sup>†‡</sup>, Inès Li de la Sierra-Gallay<sup>§</sup>, Alain F. Chaffotte<sup>¶</sup>, Benjamin Pluvinage<sup>‡</sup>,<sup>‡</sup> Patrick Weber<sup>||</sup>, Ahmed Haouz<sup>||</sup>, Jean-Marie Dupret<sup>‡</sup>, and Fernando Rodrigues-Lima<sup>‡</sup><sup>§</sup>

From the <sup>‡</sup>Université Paris Diderot, Sorbonne Paris Cité, Unité de Biologie Fonctionnelle et Adaptative, CNRS EAC4413, 75013 Paris, France, <sup>§</sup>Université Paris-Sud, Institut de Biochimie et Biophysique Moléculaire et Cellulaire, CNRS UMR8619, 91405 Orsay, France, <sup>¶</sup>Institut Pasteur, Unité de Résonance Magnétique Nucléaire des Biomolécules, 75015 Paris, France, and <sup>||</sup>Institut Pasteur, Plateforme 6, CNRS-UMR3528, 25 rue du Dr. Roux, F-75724 Paris, France

**Background:** Catalytic activity of prokaryotic and eukaryotic arylamine N-acetyltransferases (NATs) relies on a strictly conserved catalytic Cys-His-Asp triad.

**Results:** Structural and biochemical studies identified a functional NAT with a Cys-His-Glu catalytic triad.

**Conclusion:** The catalytic triad of these acetyltransferases is more plastic than previously believed.

**Significance:** (BACCR)NAT3 represents the first known functional acetyltransferase with a non-canonical Cys-His-Glu catalytic triad.

Arylamine N-acetyltransferases (NATs), a class of xenobiotic-metabolizing enzymes, catalyze the acetylation of aromatic amine compounds through a strictly conserved Cys-His-Asp catalytic triad. Each residue is essential for catalysis in both prokaryotic and eukaryotic NATs. Indeed, in (HUMAN)NAT2 variants, mutation of the Asp residue to Asn, Gln, or Glu dramatically impairs enzyme activity. However, a putative atypical NAT harboring a catalytic triad Glu residue was recently identified in *Bacillus cereus* ((BACCR)NAT3) but has not yet been characterized. We report here the crystal structure and functional characterization of this atypical NAT. The overall fold of (BACCR)NAT3 and the geometry of its Cys-His-Glu catalytic triad are similar to those present in functional NATs. Importantly, the enzyme was found to be active and to acetylate prototypic arylamine NAT substrates. In contrast to (HUMAN) NAT2, the presence of a Glu or Asp in the triad of (BACCR)NAT3 did not significantly affect enzyme structure or function. Computational analysis identified differences in residue packing and steric constraints in the active site of (BACCR)NAT3 that allow it to accommodate a Cys-His-Glu triad. These findings overturn

the conventional view, demonstrating that the catalytic triad of this family of acetyltransferases is plastic. Moreover, they highlight the need for further study of the evolutionary history of NATs and the functional significance of the predominant Cys-His-Asp triad in both prokaryotic and eukaryotic forms.

Arylamine N-acetyltransferases (NATs<fnref n="4" />; EC. 2.3.1.5) are phase II xenobiotic-metabolizing enzymes. These xenobiotic-metabolizing enzymes play an important role in the biotransformation and biological fates of aromatic amine compounds like drugs, carcinogens, and toxins (1, 2). NATs are structurally related to transglutaminases and cysteine proteases and possess a strictly conserved Cys-His-Asp catalytic triad (3–5). In functional NATs, this catalytic triad allows the N-acetylation of arylamines and/or O-acetylation of arylhydroxylamines through a Ping Pong Bi Bi mechanism (6, 7).

Site-directed mutagenesis experiments have shown that the NAT catalytic reaction unequivocally requires each of the Cys, His, and Asp catalytic residues (8–10). Accordingly, all known functional prokaryotic and eukaryotic NATs harbor a Cys-His-Asp catalytic triad (11, 12). Natural variations in the catalytic triad residues are poorly documented among the NAT family. So far, one unique natural modification of a catalytic triad residue has been reported and characterized. A single nucleotide polymorphism (364G→A) in a (HUMAN)NAT2 isozyme (NAT2*12D) leads to an inactive NAT2 variant in which the catalytic triad Asp residue is replaced by Asn (10). Surprisingly, mutation of the catalytic Asp to the chemically analogous Glu in (HUMAN)NAT2 produced a variant with very low enzyme

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

<footnote>* This work was supported in part by the Université Paris Diderot-Paris 7, Délégation Générale de l'Armement (DGA), the Institut Pasteur Paris, the CNRS and French Infrastructure for Integrated Structural Biology (FRISBI) Grant ANR-10-INSB-05-01.</footnote>

<footnote>§ This article contains supplemental Fig. S1.</footnote>

<footnote>The atomic coordinates and structure factors (code 4DMO) have been deposited in the Protein Data Bank (http://wwpdb.org/).</footnote>

<footnote n="1">¹ Supported by a fellowship from the Université Paris Diderot-Paris 7. Present address: Molecular Neuropharmacology and Genetics Laboratory, Dept. of Neuroscience and Pharmacology, University of Copenhagen, DK-2200 Copenhagen, Denmark.</footnote>

<footnote n="2">² Supported by a fellowship from the DGA. Present address: Biochemistry and Microbiology, University of Victoria, Victoria, British Columbia BC V8P 5C2, Canada.</footnote>

<footnote n="3">³ To whom correspondence should be addressed: Université Paris Diderot, Sorbonne Paris Cité, Unité de Biologie Fonctionnelle et Adaptative, CNRS EAC4413, 4 Rue Marie Andrée Lagroua Weill Halle, 75013 Paris, France. Tel.: 33-1-57-27-83-32; E-mail: fernando.rodrigues-lima@univ-paris-diderot.fr.</footnote>

<footnote n="4">⁴ The abbreviations used are: NAT, arylamine N-acetyltransferase (unless specified, all NAT enzymes referred to in the text are named according to the current official NAT nomenclature); (BACCR)NAT3, *B. cereus* NAT3; AcCoA, acetyl coenzyme A; PNPA, *para*-nitrophenyl acetate; IAA, iodoacetamide; GdmCI, guanidine hydrochloride; CAPS, 3-(cyclohexylamino)propanesulfonic acid; r.m.s.d., root mean square deviation.</footnote>
A Functional NAT with a Non-canonical Catalytic Triad

activity (2% of the wild-type activity). This Asp-to-Glu variant is
also expressed at very low levels, suggesting that the substitu-
tion alters the structure of the enzyme (10).

Analyses of functional prokaryotic and eukaryotic NATs and
those with triad variants prompted the hypothesis that a strict Cys-
His-Asp catalytic triad is required for catalysis by this family of
acetyltransferases (10, 11). Interestingly, however, among the
three *nat* genes present in *Bacillus cereus*, one encodes for a
putative NAT isoenzyme, arylamine *N*-acetyltransferase 3
((BACCR)NAT3), with an atypical Cys-His-Glu catalytic triad (9).
The catalytic activity and structure of this atypical NAT have not
been studied presumably because it was assumed to be inactive.

Here, x-ray crystallography, site-directed mutagenesis, bio-
physical, and computational approaches were combined to
characterize this atypical NAT protein. We show that, despite
the presence of a Glu residue in its catalytic triad, this atypical
NAT is functional and able to acetylate prototypic NAT sub-
strates. In contrast to (HUMAN)NAT2, (BACCR)NAT3 can
support either a Glu or Asp residue in its catalytic triad with no
significant structural or catalytic impact. Our data show that,
contrary to current knowledge, the plasticity of the active site of
this family of acetyltransferases also extends to their catalytic
triad residues. More broadly, it raises the question of the evo-
lutionary history of NAT enzymes and in particular of the func-
tional significance of the predominance of the Cys-His-Asp
triad in both prokaryotic and eukaryotic isozymes.

EXPERIMENTAL PROCEDURES

*Chemical Compounds, Strains, and Culture Conditions*—
4-Aminosalicylic acid, 2-aminofluorene, isoniazid, acetyl
coenzyme A (AcCoA), *para*-nitrophenyl acetate (PNPA), iodo-
acetamide (IAA), isopropyl *β*-D-1-thiogalactopyranoside, imida-
zole, lysozyme, kanamycin, and guanidine hydrochloride
(GdmCl) were purchased from Sigma. Substrates were pre-
pared in Milli-Q water or dimethyl sulfoxide (DMSO; Eurome-
dex) at 100 *mM* and stored at −20 °C until use. Protease
inhibitor mixture (100×) was prepared with 100 *mM* 4-(2-ami-
noethyl)benzenesulfonyl fluoride hydrochloride and 2 *mM* leu-
peptin (Euromedex). Other compounds were provided by
Sigma unless specified. *Escherichia coli* BL21 strains were cul-
tured in Luria broth (LB) Miller medium at 37 °C. *B. cereus*
strain ATCC 14579 was kindly provided by Pr. I. Martin-Vers-
traete (Institut Pasteur, Paris, France).

*Amino Acid Sequence Analysis of B. cereus NAT3*—The *nat3*
open reading frame (ORF) coding for the (BACCR)NAT3 was
identified using BLASTp. Multiple sequence alignments with
other NATs were achieved with Muscle at the European Bioin-
formatics Institute. All parameters were set to defaults. Align-
ment outputs were generated with ESPript V2.2.

*Cloning of B. cereus nat3 ORF and Construction of an E123D*
*Mutant*—Cloning of the ORF coding for (BACCR)NAT3
wild type (WT) was reported in Kubiak *et al.* (13). The
(BACCR)NAT3 E123D mutant was generated by site-directed
mutagenesis. Two independent PCRs were conducted on genomic
DNA, one using primers (BACCR)NAT3_for (5’-CCGGATCC-
ATGACCGACTTTCAAAAACA-3’) and E123D_rev (5’-
CGAATCCTACGTCAATTAAAT-3’) and the other using
primers E123_for (5’-ATTTAATTGACGTAGGATTCG-3’)

and (BACCR)NAT3_rev (5'-GGGAGCTCTTATAATGTAAT-
TCCGTT-3'). A third PCR using primers (BACAN)NAT3_for
and (BACAN)NAT3_rev was used to obtain a single amplicon
containing the E123D mutation. PCR products were double
digested using BamHI and XhoI (New England Biolabs) and
cloned into a pET28(a) plasmid for further expression of His₆-
tagged fusion recombinant proteins. Cloning was verified by
DNA sequencing.

*Expression and Purification of Recombinant (BACCR)NAT3* *WT and E123D Enzymes*—*E. coli* BL21 strains containing the (BACCR)NAT3 WT or E123D constructions were grown at 37 °C until absorbance at 600 nm reached 0.7. Addition of 500 µM isopropyl *β*-D-1-thiogalactopyranoside for 5 h at 25 °C started the induction of NAT protein overexpression. Bacteria were then resuspended in phosphate-buffered saline (PBS), pH 7.4, 10 mg of lysozyme, 0.1% Triton X-100, and 1× protease inhibitor mixture prior to sonication at 4 °C (3 min; 8-s pulse on, 30-s pulse off; Branson Digital Sonicator model 450). Bacterial debris were pelleted (30 min, 12,000 × g, 4 °C), and supernatants were incubated with Chromatrix nickel-nitrilotriacetic acid beads (Jena Bioscience) for a minimum of 2 h at 4 °C. Extensive washes were achieved on columns successively with PBS with 0, 20, and 30 mM imidazole before elution of His₆-tagged NAT enzymes with PBS with 300 mM imidazole. Finally, proteins were reduced with 10 mM dithiothreitol (DTT) and dialyzed overnight against 25 mM Tris-HCl, pH 7.5. Proteins were digested with thrombin and further purified by ion exchange chromatography (Mono Q column, ÄKTA system, GE Healthcare) as described previously (13). Purification yield was estimated by SDS-PAGE, and protein concentration was determined by absorbance at 280 nm.

*Crystal Structure Determination*—Crystallization and data collection of the (BACCR)NAT3 wild-type protein were reported previously (13). Briefly, crystals were obtained with the hanging drop-vapor diffusion method with a 1:1 mixture of 30 mg·ml⁻¹ (BACCR)NAT3 protein with 1.6 *m* sodium citrate, pH 6.5, 0.28 *m* NDSB-221 as additive (Hampton Research) at 18 °C. The best crystals were flash cooled in liquid nitrogen with a 1:1 volume ratio of Paratone and paraffin oil as cryoprotectant before diffraction data collection at the Swiss Light Source synchrotron (Beamline PX3, Paul Scherrer Institute, Switzerland). The crystallographic parameters and data statistics were reported in Kubiak *et al.* (13). Briefly, data collection was achieved at 0.984 Å with a 1° oscillation per image on a charge-coupled device mar225 detector. XDS was used for data integration (14). The (BACCR)NAT3 structure was resolved by molecular replacement using the program Phaser (15) implemented in the CCP4 program suite (16). A model of (BACCR)NAT3 was built using CHAINSAW (17) with (BACAN)NAT1 structure (Protein Data Bank code 3LNB) as the search model. The protein structure was obtained by iterative manual rebuild using Coot (18) and refined by Twin Lattice Symmetry and restrained refinement using REFMAC5 (19). After addition of water molecules, a final refinement step was achieved using BUSTER 2.0 software (20). The stereochemical quality of the model was followed during the refinement process using the program MolProbity (21). Water molecules were added with a 3.0σ cut-off and manual verification before final validation was per-
<caption>**TABLE 1**
Refinement statistics and model parameters of (BACCR)NAT3 WT</caption>

SLS, Swiss Light Source.

<table><thead><tr><th></th><th></th></tr></thead><tbody><tr><td><strong>Data collection</strong></td><td></td></tr><tr><td>Beam line</td><td>X06DA, SLS</td></tr><tr><td>Wavelength (Å)</td><td>0.984</td></tr><tr><td>Space group</td><td>C121</td></tr><tr><td>Cell dimensions</td><td></td></tr><tr><td>a, b, c (Å)</td><td>90.43, 44.52, 132.97</td></tr><tr><td>α, β, γ (°)</td><td>90.0, 103.8, 90.0</td></tr><tr><td>Resolution (Å)</td><td>20.0-2.14 (2.27-2.14)<sup>a</sup></td></tr><tr><td>R<sub>meas</sub> (%)<sup>b</sup></td><td>10.9 (103.1)</td></tr><tr><td>I/σ(I)</td><td>14.1 (2.0)</td></tr><tr><td>Completeness (%)</td><td>99.2 (96.0)</td></tr><tr><td>Multiplicity</td><td>7.4 (7.3)</td></tr><tr><td>No. of reflections</td><td>208,589 (29,418)</td></tr><tr><td>No. of unique reflections</td><td>28,411 (4,348)</td></tr><tr><td>Wilson B-factor (Å<sup>2</sup>)</td><td>37.8</td></tr><tr><td><strong>Refinement</strong></td><td></td></tr><tr><td>Resolution (Å)</td><td>20.0-2.14</td></tr><tr><td>R/R<sub>free</sub> (%)</td><td>20.4/25.5</td></tr><tr><td>No. of non-hydrogen atoms</td><td></td></tr><tr><td>Protein</td><td>4,064</td></tr><tr><td>Solvent</td><td>128</td></tr><tr><td>Ligands</td><td>42</td></tr><tr><td>B-factor (Å<sup>2</sup>)</td><td></td></tr><tr><td>Protein</td><td>27.8</td></tr><tr><td>Water molecules</td><td>24.1</td></tr><tr><td>Ligands</td><td>60.4</td></tr><tr><td>r.m.s.d. bond lengths (Å)</td><td>0.012</td></tr><tr><td>r.m.s.d. bond angle (°)</td><td>1.2</td></tr><tr><td>Ramachandran plot (%)</td><td></td></tr><tr><td>Favored</td><td>97.0</td></tr><tr><td>Allowed</td><td>3.0</td></tr><tr><td>Disallowed</td><td>0.0</td></tr></tbody></table>

a Numbers in parentheses are for the highest resolution shell.

b $R_{meas} = \Sigma_h \sqrt{(n/n-1)} \Sigma_i |I_{hi} - I_{ph}| / \Sigma_h \Sigma_i |I_{hi}|$ where $I_{hi}$ is the ith observation of the reflection h, $I_{ph}$ is the mean intensity of reflection h, and n is the multiplicity of reflection h.

formed with PROCHECK (22). The final refinement statistics and model parameters are shown in Table 1. All structural representations were generated with PyMOL (23). Atomic coordinates and structure factors of the (BACCR)NAT3 protein have been deposited in the Protein Data Bank under accession code 4DMO.

**Determination of Apparent Kinetic Parameters of (BACCR)NAT3 and the E123D Variant Enzymes**—Steady-state kinetic experiments were conducted using the 5,5'-dithiobis(nitrobenzoic acid) assay (24) to determine the apparent kinetics of (BACCR)NAT3 and E123D variant enzymes parameters $K_m^{app}$ and $k_{cat}$ toward 4-aminosalicylic acid, 2-aminofluorene, and isoniazid. Initial velocities ($V_i$) were determined in the presence of arylamine substrate at various concentrations, 400 µM AcCoA, and 0.5 µg·ml⁻¹ enzyme in 25 mM Tris-HCl, pH 7.5. Alternatively, various concentrations of acetyl donors were used with 500 µM 4-aminosalicylic acid to determine the apparent kinetic parameters toward AcCoA. The PNPA assay described by Cleland and Hengge (25) was used to obtain kinetic parameters toward the acetyl donor PNPA. Control experiments were conducted in the absence of arylamine substrate. $V_i$ values were plotted against substrate concentration and fitted to the Michaelis-Menten equation using Kaleidagraph 3.5 (Synergy Software). Results of triplicate experiments are shown including standard deviations.

**Circular Dichroism (CD) Spectra and Chemical Unfolding Followed by CD**—Recombinant enzyme samples were extensively dialyzed against 10 mM sodium phosphate, pH 7.5 buffer before recording CD spectra in an Aviv215 spectropolarimeter. Far-UV CD spectra (180–260 nm) were acquired with 505 and

463 µg·ml⁻¹ final concentrations of (BACCR)NAT3 and E123D variant, respectively, through a cylindrical cell with a 0.02-cm path length. The ellipticity signal was recorded each 0.5 nm (1-nm bandwidth) with an integration time of 1 s per step. The average spectra of four successive scans was normalized to the protein concentration, and ellipticity was converted to differential molar extinction coefficient $\Delta\epsilon$ per residue. Secondary structure quantitative composition was deduced from the decomposition of the spectra using the CONTIN analysis included in the CDPro package with the 29-protein basis of Johnson and co-worker (26). Near-UV CD spectra (250–350 nm) were obtained from an average of five successive scans between 250 and 350 nm through a 1-cm path length rectangular cell with 2.016 and 1.852 mg·ml⁻¹ final concentration of (BACCR)NAT3 and E123D variant, respectively. Spectra were normalized to the protein concentration as $\Delta\epsilon$ per chain.

Native and 5 mM GdmCl totally unfolded proteins both at 0.1 mg·ml⁻¹ were prepared in 20 mM Tris-HCl, pH 7.8. Various amounts of native and unfolded proteins were mixed in a 1-ml final volume to obtain GdmCl concentrations ranging from 0 to 5 mM with a 0.1 mM steps. The samples were held for 24 h at 25 °C before recording ellipticity at 222 nm using a 0.5-cm path length rectangular cell (10-s integration time, 2-nm bandwidth). The ellipticity at 222 nm was expressed as a function of GdmCl concentration and then fitted to the following two-state transition model equation.

$$ \theta = \frac{\theta_N + aC + (\theta_U + bC)\exp[-(\Delta G_{H2O} + mC)/RT]}{1 + \exp[-(\Delta G_{H2O} + mC)/RT]} \quad (\text{Eq. 1}) $$

where C is the concentration of GdmCl (M), $\Delta G_{H2O}$ is the unfolding free energy in the absence of denaturant, and m is the dependence of $\Delta G$ on denaturant concentration. a and b are the slopes of the ellipticity base line (supposedly linear) of the native and unfolded states, respectively. Finally, the ellipticity signal was normalized to the total amplitude difference between the fully unfolded protein and native protein (27).

**pH-dependent Inactivation by Iodoacetamide**—Recombinant enzyme was mixed with either 50 mM MES, pH 6.0–6.5; 50 mM Tris-HCl, pH 7–9; or 50 mM CAPS, pH 9.5–10. Addition of 1 mM IAA to a solution containing 15 µM enzyme started the inactivation reaction at 25 °C. The percentage of residual activity was measured at different incubation times by adding 5-µl aliquots to a solution containing 25 mM Tris-HCl, pH 7.5, 2 mM PNPA, and 500 µM 4-aminosalicylic acid as described above. The percentage of residual activity fitted a single exponential curve allowing us to determine the first-order inactivation constant $k_{inact}$, which was corrected by the $k_{inact}$ determined in the absence of IAA. ln($k_{inact}$) was plotted as a function of pH. The catalytic cysteine $pK_a$ value was determined using the Henderson-Hasselbalch equation (28).

$$ k_{\text{inact}} = k_{\text{min}} + \frac{(k_{\text{max}} - k_{\text{min}})}{1 + 10^{pK_a - \text{pH}}} \quad (\text{Eq. 2}) $$

where $k_{max}$ is the highest and $k_{min}$ is the lowest measured $k_{inact}$ value, respectively. Standard deviations were calculated from three independent experiments.

<figure><img src="image_4.png" /></figure>
A Functional NAT with a Non-canonical Catalytic Triad

*Phylogenetic Analysis of NAT Sequences*—Amino acid sequences were retrieved with BLASTp on the non-redundant sequence databases using the (BACCR)NAT3 sequence as a query. The 250 protein sequences with the highest scores were aligned. Of these, one strain representative of each species with different NAT isoform sequences was selected. Sequences with lower scores but presenting a putative Cys-His-Glu catalytic triad were added to this selection. All sequences were aligned using Muscle with the (BACCR)NAT1, (BACCR)NAT2, and (BACCR)NAT3 sequences and sequences of previously characterized NATs. The phylogenetic tree was constructed using a distance method implemented in MEGA (29) through a neighbor-joining analysis computed using the Dayhoff model. A total of 500 bootstrap replications were performed to determine statistical support for clades.

RESULTS AND DISCUSSION

*Analysis of B. cereus Strain ATCC 14579 NAT Sequences Reveals an Isoform with a Cys-His-Glu Catalytic Triad—B. cereus* strain ATCC14579 possesses three NAT homologue sequences (BACCR)NAT1, (BACCR)NAT2, and (BACCR)NAT3, which share 94, 83, and 49% sequence identity with the characterized *Bacillus anthracis* NAT1, NAT2, and NAT3, respectively, indicating that they are orthologous (5, 9, 30). We aligned the sequences of these three NATs with characterized eukaryotic and prokaryotic NATs (Fig. 1); all three *B. cereus* NAT sequences possess the characteristic NAT functional motifs, suggesting that these proteins belong to the NAT family of acetyltransferases (11, 31). Importantly, as reported previously, one of the *B. cereus* NAT sequences, (BACCR)NAT3, contains Glu in place of the canonical catalytic triad Asp residue (5, 9).

*Expression of (BACCR)NAT3 and an E123D (BACCR)NAT3 Variant*—The only characterized NAT enzyme harboring a putative Cys-His-Glu catalytic triad is (BACCAN)NAT3, a short NAT homologue devoid of NAT activity (5, 30). Previous evidence suggests that NAT sequences harboring a Cys-His-Glu catalytic triad are devoid of NAT activity (9, 10, 30). To characterize the (BACCR)NAT3, we cloned and expressed the protein in *E. coli* and purified it to homogeneity. In parallel, a variant of (BACCR)NAT3 in which the catalytic Glu was replaced by an Asp residue (E123D variant) was also cloned, expressed, and purified. Both proteins were soluble and had identical purification yields (5 mg/liter of culture). SDS-PAGE revealed a unique band of ~34 kDa consistent with the predicted (ProtParam) molecular mass of 33,465 and 33,450 Da for the wild-type and mutant proteins, respectively (Fig. 2A). Altogether, these data suggest that, in contrast to (HUMAN)NAT2 isozyme, the presence of a Glu or Asp residue in the catalytic triad of (BACCR)NAT3 does not impact protein stability and does not lead to aggregation during expression (10).

Further analysis of (BACCR)NAT3 and its E123D variant by circular dichroism revealed that their secondary and tertiary structures were very similar to one another (Fig. 2, B and C). The unfolding profiles of both proteins showed an unusual transient dramatic ellipticity decrease centered on 0.6–0.7 M GdmCl (Fig. 2D). At higher denaturant concentration, a regain in ellipticity was observed, and at GdmCl concentrations above 1.1 M, the signal change appears as a classical unfolding transition. Both proteins showed very similar unfolding profiles (including the large transient 0.1–1 M GdmCl ellipticity loss). Fitting the ellipticity change between 1.2 and 5 M denaturant to the classical two-state model led to very similar stability characteristics with a midpoint at 1.83 M ($\Delta G_{H2O}$ = 6377.0; *m* = 3486.3) and 1.72 M ($\Delta G_{H2O}$ = 8043.5; *m* = 4683.0) GdmCl for the wild-type and E123D isozymes, respectively. The large ellipticity loss preceding the standard unfolding transition coincided with a very significant transient increase of light scattering intensity reflected by the dynode high voltage during CD measurement. The marginal loss in ellipticity at small denaturant concentration may denote the transient accumulation of ordered or unordered soluble aggregates (fibrils?) induced by a helix → coil or helix → β-structure conformational transition. This atypical phenomenon could be related to the “denaturation trough” initially reported by Goldberg and co-workers (32) and then by other authors (33–35). Although the (HUMAN)NAT2 D122N protein has been suggested to retain its global folding over molecular dynamics experiments, mutagenesis studies showed that the (HUMAN)NAT2 D122E variant was unstable and prematurely degraded (10, 36). In contrast, our data demonstrate that the presence of either an Asp or a Glu residue in the catalytic triad appears not to affect the folding properties of (BACCR)NAT3.

*A Cys-His-Glu Catalytic Triad Is Consistent with NAT Activity in (BACCR)NAT3 Wild Type: a Kinetic Comparison with the E123D Mutant*—As both the (BACCR)NAT3 and (BACCR)NAT3 E123D proteins appeared properly folded, we investigated their potential to N-acetylate three model aromatic amine substrates (4-aminosalicylic acid, 2-aminofluorene, and isoniazid). (BACCR)NAT3 displayed significant N-acetylation activity toward these three NAT substrates (Table 2). Restoring the Asp residue at position 123 slightly decreased the turnover of the enzyme ($k_{cat}$) by 1.3–1.7 times and the affinity constant ($K_{m}^{app}$) by 2.1–6.8 times. (BACCR)NAT3 wild-type catalytic efficiency ($k_{cat}/K_{m}^{app}$) was 1.6–4.4 times lower than that of the E123D mutant, but both enzymes were in the range of mean catalytic efficiency values observed for other bacterial NATs (30, 37). Similarly, the $K_{m}^{app}$ value of (BACCR)NAT3 for the acetyl donor PNPA was 2 times higher than that of the mutant, and the catalytic efficiency was 2.5 times higher. However, only small differences were observed between the two enzymes in

FIGURE 1. Alignment of *B. cereus* NAT sequences with sequences of functionally characterized NATs. All three NAT sequence homologues from *B. cereus* ((BACCR)NAT1, (BACCR)NAT2, and (BACCR)NAT3) possess the two conserved NAT functional motifs (black bars). The catalytic Cys and His residues (★) are fully conserved throughout all *B. cereus* NAT homologues. However, (BACCR)NAT3 and (BACAN)NAT3 harbor a Glu residue instead of the canonical Asp residue at catalytic position 123 (▲). The alignment suggests a three-domain (α-helical bundle (1–85), β-barrel (86–192), and C-terminal α/β lid (193–263)) organization of the *B. cereus* NAT homologues by comparison with the secondary structure of (SALTY)NAT1 indicated on top of the alignment. The region corresponding to a mammalian-like insertion loop is underlined in gray. The alignment was performed by the Muscle algorithm, and the output was achieved with ESPript 2.0. Numbers indicate amino acid position in (SALTY)NAT1.

<figure><img src="image_5.png" /></figure>
A Functional NAT with a Non-canonical Catalytic Triad

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

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

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

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

FIGURE 2. Expression yields and structural similarities between (BACCR)NAT3 wild-type and E123D enzymes. A, separation profile of purified recombinant (BACCR)NAT3 wild-type and E123D proteins (1 µg) on a 12.5% SDS-PAGE gel stained with Coomassie Blue. Lane M, molecular mass markers. B, far-UV CD spectra. Decomposition of the spectra (in box) reveals identical fraction of secondary structure content (expressed in percent) between the wild-type (black line) and the E123D mutant enzymes (gray line). H, helix; S, strand; Unrd, unordered. C, near-UV CD spectra (black line, wild type; gray line, E123D mutant). D, guanidinium chloride-induced unfolding of (BACCR)NAT3 wild-type and E123D enzymes followed at 222 nm. The solid lines correspond to the simulated transition curved built from the fit of the data to a two-state transition between 1.2 and 5 M GdmCl (black boxes and line, wild type; open circles and gray line, E123D mutant). The experimental conditions are detailed under "Experimental Procedures."

<caption>TABLE 2</caption>

Kinetic parameters of (BACCR)NAT3 WT and E123D enzymes for prototypic arylamines and acetyl donors

4-AS, 4-aminosalicylic acid; 2-AF, 2-aminofluorene; INH, isoniazid.

<table><thead><tr><th>Compound</th><th></th><th>K<sub>m</sub><sup>app</sup>, mM</th><th>k<sub>cat</sub>, s<sup>-1</sup></th><th>k<sub>cat</sub>/K<sub>m</sub><sup>app</sup>, x10<sup>6</sup> M<sup>-1</sup>s<sup>-1</sup></th></tr></thead><tbody><tr><td rowspan="3">Acetyl donor</td><td rowspan="3">Arylamine</td><td></td><td>(BACCR)NAT3 wt</td><td></td></tr><tr><td>4-AS</td><td>1.34 ± 0.3</td><td>56.6 ± 6.7</td><td>42.4 ± 7.7</td></tr><tr><td>2-AF</td><td>5.8 ± 1.1</td><td>67.4 ± 9.2</td><td>11.5 ± 1.8</td></tr><tr><td rowspan="2">(400 µM AcCoA)</td><td>1NH</td><td>76.9 ± 9.8</td><td>289.3 ± 30.5</td><td>3.76 ± 0.44</td></tr><tr><td></td><td>0.050 ± 0.006</td><td>19.3 ± 0.5</td><td>390 ± 28</td></tr><tr><td>AcCoA</td><td rowspan="2">(500 µM 4-AS)</td><td>11.09 ± 0.88</td><td>7.6 ± 0.5</td><td>0.69 ± 0.04</td></tr><tr><td>PNPA</td><td></td><td></td><td></td></tr><tr><td rowspan="3">(400 µM AcCoA)</td><td>4-AS</td><td>0.64 ± 0.04</td><td>43.9 ± 1.2</td><td>68.4 ± 3.0</td></tr><tr><td>2-AF</td><td>0.86 ± 0.15</td><td>44.1 ± 3.3</td><td>51.3 ± 6.4</td></tr><tr><td>1NH</td><td>14.55 ± 0.12</td><td>167.1 ± 0.7</td><td>11.5 ± 0.1</td></tr><tr><td>AcCoA</td><td rowspan="2">(500 µM 4-AS)</td><td>0.052 ± 0.004</td><td>14.3 ± 0.3</td><td>275 ± 13</td></tr><tr><td>PNPA</td><td>4.74 ± 0.99</td><td>8.0 ± 1.2</td><td>0.17 ± 0.12</td></tr></tbody></table>

kinetic parameters for the physiological acetyl donor AcCoA (Table 2).

Consistent with previous reports for other NATs, PNPA was a less efficient acetyl donor than AcCoA for both NATs we studied (37). AcCoA has been shown to possess numerous stabilizing interactions all over the active site cleft in contrast to

aromatic amines that interact deeply in the active site to reach the acetyl-cysteine intermediate for acetylation (38–41). Thus, the differences in kinetic parameters observed for aromatic amine substrates and the aromatic amine structural analogue PNPA suggest that the shape of the active site proximal to the catalytic triad might be affected; a similar mechanism has been proposed for the (HUMAN)NAT2 D122N mutant (36). Glu and Asp residues share very similar chemical properties and have similar side chain pK<sub>a</sub> values of 4.15 and 3.71, respectively. The third catalytic residue of the triad has been suggested to act as a base assisting the deprotonation of the two other catalytic residues; a change in the catalytic mechanism would therefore be expected to modify the reactivity of the catalytic cysteine (42, 43). Hence, we conducted a pH-dependent inactivation experiment by IAA to determine the reactivity/pK<sub>a</sub> of the catalytic Cys<sup>69</sup> in both (BACCR)NAT3 and the E123D variant. The pK<sub>a</sub> values of the catalytic cysteine of the wild-type and E123D mutant enzymes were found to be similar (8.26 ± 0.11 and

<figure><img src="image_10.png" /></figure>
A Functional NAT with a Non-canonical Catalytic Triad

7.80 ± 0.14, respectively) (Fig. 3). This suggests that the presence of a catalytic Asp or Glu residue has no significant impact on the reactivity of the catalytic cysteine of (BACCR)NAT3. Moreover, the IAA/pH inactivation profile is more similar to that of *Mycobacterium tuberculosis* NAT (whose catalytic cysteine p*K*<sub>a</sub> lies above 10) than to that of hamster NAT2 for which the catalytic cysteine p*K*<sub>a</sub> is 5.23. Therefore, the catalytic mechanism of (BACCR)NAT3 is unlikely to rely on a thiolate/imida-

zolium ion pair as shown for hamster NAT2 (42) but rather on a general base catalysis mechanism in which the catalytic Cys⁶⁹ is deprotonated in the presence of substrate at physiological pH as proposed for *M. tuberculosis* NAT (Fig. 4) (43).

Mutation of individual catalytic residues led to a dramatic loss of activity in NATs (10, 38). For the first time, we demonstrate the existence of a non-canonical but functional Cys-His-Glu catalytic triad in an NAT. Our data suggest that this non-canonical triad has no significant impact on activity or substrate-specific catalytic mechanisms when compared with other characterized prokaryotic NATs. This is an important finding for documenting the catalytic diversity of NAT-related enzymes. Indeed, only two catalytic triad variants of the papain-like cysteine peptidases cathepsins L have been reported and are associated with normal activity levels but greatly altered substrate specificity (44). The catalytic triads of transglutaminases and deubiquitinating enzymes are strictly conserved (45). Nevertheless, two esterases (serine proteases) have been reported to have normal activity levels and substrate specificity despite an Asp-to-Glu mutation in the Cys-His-Asp catalytic triad (46). Alternatively, normal activity levels have been reported in particular inteins with a modified catalytic residue but with an accompanying change in the catalytic mechanism (47, 48).

*Crystal Structure of (BACCR)NAT3 Wild-type Enzyme: Overall Structure and Catalytic Triad Geometry*—Although mutagenesis and computational studies have been conducted on catalytic triad variants of NAT enzymes, no structural determination studies have been reported. Thus, we have solved the

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

**FIGURE 3.** pH-dependent inactivation of (BACCR)NAT3 by IAA. The first-order inactivation constant (k_inact) was determined between pH 6 and 9.5 by measuring the ln of residual activity after 1 mM IAA inactivation of the enzyme. The non-linear fit of data to the Henderson-Hasselbalch equation (Equation 2) indicates pKₐ values of 8.26 and 7.80 for the Cys⁶⁹ catalytic residue of the wild type (black curve) and E123D mutant (red curve), respectively. Standard error bars are shown for three independent experiments.

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

**FIGURE 4.** Putative catalytic mechanism for the Glu- and Asp-containing (BACCR)NAT3 enzymes. The pKₐ values of the catalytic cysteine residue suggest that it exists mostly as a thiol at physiological pH (≤8), which would be inconsistent with the nucleophilic attack on the AcCoA molecule (first reaction step). Thus, as suggested by Sikora *et al*. (43) for the *M. tuberculosis* NAT enzyme, we propose that the Glu/Asp¹²³ and His¹⁰⁸ dyad act as general base that facilitates the formation of a thiolate. In the absence of AcCoA in the active site, the thiol proton is shared between Cys⁶⁰ and His¹⁰⁸ (dashed line). The AcCoA binding in the active site could displace this equilibrium and help in the deprotonation of the thiol, leading to the nucleophilic attack of the thiolate on the AcCoA. The following catalysis steps are likely to involve the canonical tetrahedral intermediate states (in brackets) that are known to occur in the Ping Pong Bi Bi mechanism of NAT enzymes. The active site is represented by a gray dashed curve, the His¹⁰⁸ and Glu/Asp¹²³ residues are buried, and the Cys⁶⁹ is solvent-accessible.

<figure><img src="image_13.png" /></figure>
A Functional NAT with a Non-canonical Catalytic Triad

<figure><img src="image_14.png" /><figcaption>FIGURE 5. **Three-dimensional structure of the (BACCR)NAT3 wild-type enzyme.** A, schematic and surface representation of the two (BACCR)NAT3 molecules in the asymmetric unit, solved at 2.14 Å. Each enzyme chain has a typical three-domain fold including an α-helical bundle (red), a β-barrel (yellow), and a C-terminal α/β lid (blue). The structure confirms that residues Cys⁶⁹, His¹⁰⁸, and Glu¹²³ form a catalytic triad (green ball and stick). No electron density was attributed to the region spanning from Lys¹⁶⁷ and Asp¹⁷⁸ in both molecules, indicating the presence of a flexible loop corresponding to the mammalian-like insertion found in human NAT enzymes. B, 2F~o~ − F~c~ density map at 1.5σ for (BACCR)NAT3 catalytic residues Cys⁶⁹, His¹⁰⁸, and Glu¹²³. Hydrogen bonds are shown as black dashes, and distances (molecule A in the asymmetric unit) are expressed in Å. Distances in molecule B are shown in parentheses. C, superimposition of (BACCR)NAT3 (red) with (SALTY)NAT1 (yellow), (PSEAE)NAT1 (blue), (MYCMR)NAT1 (orange), and (RHILO)NAT1 (green) catalytic triads with r.m.s.d. values of 0.869, 0.672, 0.825, and 0.812 Å, respectively. It reveals a conserved geometry for the non-canonical Cys-His-Glu catalytic triad in (BACCR)NAT3 compared with other canonical NAT catalytic triads. The (BACCR)NAT3 backbone is shown in orange.</figcaption></figure>

crystal structure of (BACCR)NAT3 at 2.14 Å. The Matthews
coefficient value of 2.16 Å³·Da⁻¹ and solvent content of 43.03%
are consistent with the presence of two molecules in the
asymmetric unit after molecular replacement. According to
PROCHECK, 97% of the residues are in favored regions of the
Ramachandran plot, and 3% are in additional allowed regions,
and the *R*-factor and *R*<sub>free</sub> values are 20.4 and 25.5%, respec-
tively (Table 1). Three molecules of NDSB-221 additive have
been identified in the asymmetric unit with no biological signif-
icance. No electron density was observed for the two residues
Gly⁻¹ and Ser�0 that remained after the His₆ tag thrombin cleav-
age. Clear electron density was observed for all residues except
for side chains of surfaced-distributed residues 6, 15–18, 23, 54,
88, 230, 241, and 249. The enzyme has a typical three-domain
fold consisting of an α-helical bundle (1–85), a β-barrel (86–
192), and a C-terminal α/β lid (193–263) (Fig. 5A) (11). Inter-
estingly, no electron density was observed for residues Asn¹⁶⁸–

Ala¹⁷⁷, which align with the “mammalian insertion” present in
human NATs and inferred in (BACAN)NAT1 (Fig. 1) (40, 41).
This suggests the presence of a similar highly flexible loop in
(BACCR)NAT3 and reinforces the idea that some *Bacillus*
NAT isoforms may have undergone a particular evolutionary
history (e.g. horizontal gene transfer from a eukaryotic host,
giving rise to eukaryotic features in these bacterial NATs). The
crystal structure of (BACCR)NAT3 shows unambiguous and
strong electron densities for residues Cys⁶⁹, His¹⁰⁸, and more
importantly the non-canonical Glu¹²³, shaping a catalytic triad
(Fig. 5B). Furthermore, this catalytic triad is highly superimpos-
able with the canonical Cys-His-Asp triad of other NATs as
shown by r.m.s.d. values between 0.627 and 0.869 Å. These
values are ~4-fold lower compared with r.m.s.d. values
between the catalytic triads of canonical NATs likely due to the
presence of a larger Glu residue in (BACCR)NAT3 rather than
to a different triad conformation as all catalytic triads are sim-
A Functional NAT with a Non-canonical Catalytic Triad

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

FIGURE 6. Comparison of structural properties between (BACCR)NAT3 and (HUMAN)NAT2. A, structural alignment of (BACCR)NAT3 (orange) and (HUMAN)NAT2 (cyan) catalytic triads reveals a very similar geometry. Electrostatic interactions, hydrogen bonds, and distances are represented as indicated in Fig. 5 legend. B, amino acid network interactions of the third catalytic residues in (BACCR)NAT3 (upper panel) and (HUMAN)NAT2 (lower panel). In contrast to Asp<sup>122</sup> in (HUMAN)NAT2, the Glu<sup>123</sup> main chain interacts with two water molecules (red spheres) that are further stabilized through hydrogen bonds with residues Leu<sup>121</sup>, Leu<sup>132</sup>, and Ala<sup>133</sup>. C, structural alignment representing the main differences between (BACCR)NAT3 (light orange) and (HUMAN)NAT2 (light blue) conformations. The two structures were aligned using PyMOL (The PyMOL Molecular Graphics System, Version 0.99, Schrödinger, LLC.); additional residues found in one enzyme but not the other are colored in red and blue for amino acid extensions in (BACCR)NAT3 and (HUMAN)NAT2, respectively. The structural alignment was further supported by structural alignment of the two sequences using DaliLite at the European Bioinformatics Institute.

ilarly oriented (Fig. 5C). The side chains of Cys⁶⁹, His¹⁰⁸, and Glu¹²³ are similarly oriented compared with the Asp and Cys residues of canonical NATs, and the Cys⁶⁹-His¹⁰⁸ and His¹⁰⁸-Glu¹²³ in (BACCR)NAT3 are also comparable with those found in other NATs. Altogether, these observations depict clear structural bases for the catalytic activity measured for (BACCR)NAT3 and demonstrate that the geometry of the non-canonical Cys-His-Glu catalytic triad fulfills the structural and chemical requirements of NAT catalysis in the (BACCR)NAT3 enzyme.

*Protein Packing and Catalytic Residue Environment Likely Explain the Accommodation of a Glu in the Catalytic Triad of (BACCR)NAT3*—As previously reported by Zang *et al.* (10), the loss of activity in (HUMAN)NAT2 enzyme likely results from disruption of the protein structure integrity. The structures of (BACCR)NAT3 and (HUMAN)NAT2 are highly superimposable with r.m.s.d. values of 1.483 and 0.539 Å for the whole structure and the catalytic triad, respectively. Both enzymes have a similar catalytic triad geometry, although His¹⁰⁸ appears slightly remote in (BACCR)NAT3 compared with its position in (HUMAN)NAT2 (Fig. 6A). In (BACCR)NAT3, the catalytic Glu¹²³ residue is stabilized by a complex network of hydrogen bonds established with residues Asn⁷³, His¹⁰⁸, Gly¹²⁵, and Tyr¹⁸⁴ and unexpectedly with two water molecules that are further stabilized through three hydrogen bonds to residues Leu¹²¹, Leu¹³², and Ala¹³³ (Fig. 6B, *upper panel*). In (HUMAN)NAT2, the catalytic Asp¹²² interacts with the same residues (except for Ser¹²⁵), but the main chain is stabilized by a single hydrogen

bond with Gln¹³⁰ rather than with water molecules (Fig. 6B, *lower panel*). The topology of (BACCR)NAT3 is similar to that of (HUMAN)NAT2. Despite 27 additional residues in the (HUMAN)NAT2 enzyme compared with (BACCR)NAT3, the solvent-accessible surface of both proteins is very similar with values of 13,552 and 12,485 Å², respectively (Fig. 6C). In addition, (BACCR)NAT3 has nine unresolved/absent surfaced-exposed residues in its inferred mammalian-like loop. Fifteen of the additional residues in (HUMAN)NAT2 compose a C-terminal extension that is stabilized by a 2.8-Å hydrogen bond between Thr²⁸⁹ and Ser¹²⁷ (Fig. 7A) (40). In contrast, the C-terminal amino acid of (BACCR)NAT3 is 19 Å far from the corresponding Ser¹²⁸ position. In both enzymes, this position is located on the turn of the loop carrying the third catalytic residue. Ten of the additional residues are located in the “mammalian-like insertion,” which is stabilized by several interactions with residues on domain III and hence is much more rigid than the unresolved flexible loop in (BACCR)NAT3 (Fig. 6B). This suggests a higher packing of the (HUMAN)NAT2 enzyme that could prevent structural rearrangements of the enzyme and thus preclude the accommodation of the larger Glu residue in the catalytic triad. Additionally, several differences are found in the vicinity (<5 Å) of the third catalytic residue. In particular, the large residues Gln¹³⁰ and Met¹³¹ in (HUMAN)NAT2 are replaced by the smaller Leu¹³⁰ and Pro¹³¹ residues, respectively, in (BACCR)NAT3 (Fig. 7A). These amino acid differences lead to a larger space available around the Glu residue in (BACCR)NAT3 (167 Å³) compared with the Asp residue in

<figure><img src="image_16.png" /></figure>
A Functional NAT with a Non-canonical Catalytic Triad

A (BACCR)NAT3

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

(HUMAN)NAT2

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

B

(BACCR)NAT3

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

Wild-type

(HUMAN)NAT2

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

Wild-type

$$E123D \Delta\Delta G = 2.69 \text{ kcal/mol}$$

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

3.2 Å

4.7 Å

$$D122E \Delta\Delta G = 9.05 \text{ kcal/mol}$$

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

3.1 Å

2.5 Å

**FIGURE 7.** A model proposed for the differential accommodation of Glu as a catalytic residue in (BACCR)NAT3 and (HUMAN)NAT2. A, the cavity formed by residues surrounding (<5 Å) the Glu<sup>123</sup> or Asp<sup>122</sup> catalytic residue shows a smaller size with values of 167 and 101 Å<sup>3</sup> (CastP calculation (53)) for (BACCR)NAT3 (green pocket) and (HUMAN)NAT2 (red pocket), respectively. This difference is likely due to the presence of the Pro<sup>131</sup> residue in (BACCR)NAT3 instead of Met<sup>131</sup> residue in (HUMAN)NAT2. Both enzymes are represented in the same orientation, and a vertical 180° rotation view of the cavity in (HUMAN)NAT2 is shown in the right panel. The surface of Asn<sup>72</sup> has been subtracted in (HUMAN)NAT2 for increased visibility of the cavity, but it was included for the cavity volume calculation. The hydrogen bond in (HUMAN)NAT2 between the C-terminal residue Thr<sup>289</sup> and Ser<sup>127</sup> is represented in dashes (distances in Å); the corresponding Ser<sup>128</sup> does not interact with the C-terminal region in (BACCR)NAT3. B, FoldX energy minimization models of the E123D and D122E mutants for (BACCR)NAT3 and (HUMAN)NAT2, respectively. The difference in free energy variation (ΔΔG) generated upon mutation was calculated by FoldX. The wild-type and mutant catalytic triads are represented in sticks, and the catalytic His are represented as Van der Waals spheres (1-Å radius; white). The Asp or Glu residue surface is also represented in Van der Waals spheres (1-Å radius) for (BACCR)NAT3 (orange) and (HUMAN)NAT2 (cyan) models. These models show that a steric clash occurs between His and Glu upon D122E mutation in (HUMAN)NAT2.

(HUMAN)NAT2 (101 Å<sup>3</sup>). It is likely that the interaction between the Asp<sup>122</sup> main chain and Gln<sup>130</sup> in (HUMAN)NAT2 allows fewer structural rearrangement of the two residues compared with the conformations available for the Glu<sup>123</sup> main chain and the water molecules in (BACCR)NAT3. Thus, considering the similar chemistry of the Asp and Glu residues and the results of the kinetic study of (BACCR)NAT3 enzymes, we propose that the accommodation of a Glu as a catalytic residue in NAT enzymes relies on a combination of low steric constraints, favorable protein packing, and amino acid environment.

Energy calculations and energy minimization models of the (BACCR)NAT3 E123D and (HUMAN)NAT2 D122E mutants using FoldX 3.0 and Popmusic 2.1 support this model (49, 50). The ΔΔG energy values obtained using the FoldX and Popmusic algorithms are +2.69 and 0.31 kcal/mol upon E123D mutation in (BACCR)NAT3 and +9.05 and +2.51 kcal/mol upon D122E mutation in (HUMAN)NAT2, respectively. This clearly indicates that the presence of a glutamate is energetically favorable in (BACCR)NAT3 in contrast to its destabilizing effect in (HUMAN)NAT2. Energy minimization models suggest that such an increase in ΔG is concomitant with a steric clash between the His<sup>107</sup> and Glu<sup>122</sup> catalytic residues in (HUMAN)NAT2 (Fig. 7B). Although the model seems robust, the absence of protein backbone rearrangements upon the mutation leads to a His<sup>108</sup>-Asp<sup>123</sup> distance of 4.7 Å, which is incompatible with the catalytic activity observed experimentally in the (BACCR)NAT3 E123D mutant enzyme (Fig. 7B).

The geometry of the His<sup>108</sup>, Phe<sup>125</sup>, and Tyr<sup>184</sup> residues suggests the existence of a strong network of π-π interactions (Fig. 6B). In hamster NAT2, residue Tyr<sup>190</sup> has been suggested to be important for maintaining the catalytic triad conformation by interacting with the catalytic Asp<sup>122</sup> and His<sup>107</sup> residues through hydrogen bonding and π-π interactions, respectively (51). The presence of a Ser instead of a Phe residue at position 125 in (HUMAN)NAT2 is likely to weaken this network. One

can only speculate that local structural rearrangements upon insertion of a larger Glu residue may trigger the disruption of this network. Coupled with a lower degree of freedom of the Asp<sup>122</sup>-containing loop (Ser<sup>127</sup>-Thr<sup>287</sup> interaction) and a greater compaction, this could nevertheless decrease the ability of the protein to "find" an energetically favorable conformation, leading to the disruption of the protein structure integrity and its premature degradation. Molecular dynamics simulations have only been performed on the (HUMAN)NAT2 D122N mutant, but these suggest that the enzyme possesses fewer hydrogen bonds, which could decrease its overall stability (36).

**Phylogenetic Analysis of Glu Residue Distribution in the Catalytic Triads of NAT Enzymes**—To date, only the (BACCR)NAT3 and (BACAN)NAT3 orthologous enzymes and the (HUMAN)NAT2 D122N variant have been shown to naturally harbor a mutation in a catalytic residue. A closer investigation of the diversity of the NAT catalytic triad was achieved by BLASTing the non-redundant protein with (BACCR)NAT3 as query. The 250 highest scoring sequences retrieved all correspond to NAT homologues belonging to *Bacillus* species (supplemental Fig. S1). These sequences cluster into three main clades because of sequence identity with *B. cereus* NAT1, NAT2, and NAT3 enzymes (Fig. 8). Interestingly, every strain possessing several NAT sequences systematically harbors one with a Glu residue instead of an Asp at the catalytic position 123. More strikingly, seven NAT homologue sequences are present in the low score range and also present a putative Cys-His-Glu catalytic triad, including three sequences from *Dictyostelium* species and four from other bacteria (Fig. 8). Glenn *et al.* (5) showed that NAT sequences from *B. anthracis* and *B. cereus nat* genes were closely related to NATs from slime molds, suggesting the possible involvement of horizontal *nat* gene transfer events. Thus, the presence of a non-canonical NAT enzyme in the majority of *Bacillus* strains might stem from early horizontal gene transfer events from a eukaryotic host. Still, other unrelated bacteria harbor such atypical NAT
A Functional NAT with a Non-canonical Catalytic Triad

<figure><img src="image_23.png" /><figcaption>FIGURE 8. Distribution of the non-canonical Cys-His-Glu catalytic triad throughout species. NAT homologue sequences from *Bacillus* species cluster differently from canonical bacterial NATs and are closely related to lower eukaryotes as reported previously (5). The *Bacillus* sequences cluster in three main clades corresponding to (BACCR)NAT1, (BACCR)NAT2, and (BACCR)NAT3 orthologous sequences. All NAT homologues present in the (BACCR)NAT3-containing clade have a Glu residue instead of the Asp at catalytic position 123 (in *green*). Three sequences from lower eukaryotes and four bacterial sequences also possess a putative Glu as the third catalytic residue (in *red*). Characterized NAT sequences are indicated according to the NAT nomenclature; other sequences are indicated as *genus* species strain followed by the UniProt identifier in brackets.</figcaption></figure>
# A Functional NAT with a Non-canonical Catalytic Triad

homologues. Considering the environmental adaptation role proposed for prokaryotic NATs, such findings suggest that other non-canonical NAT enzymes may be functional and that such enzymes might have particular properties, e.g. substrate specificity, conferring specific features to these organisms. The characterization of the different *B. cereus* NAT isoforms and other non-canonical enzymes will help in answering this question.

**Concluding Remarks**—Although the activity of the prokaryotic and eukaryotic NAT enzymes was known to rely on the presence of a strictly conserved Cys-His-Asp catalytic triad, we report here the characterization of a fully functional NAT enzyme possessing a non-canonical Cys-His-Glu catalytic triad. The crystal structure of this NAT enzyme from *B. cereus* ((BACCR)NAT3) suggests that the accommodation of a Glu residue might be achieved only in certain NAT enzymes that possess favorable packing and steric constraints around the catalytic position rather than depending on chemical properties of the third catalytic residue, although it remains unclear which regions or particular residues are involved. The identification of both prokaryotic and eukaryotic NAT sequences adds clues in this direction. However, it raises the question of the evolutionary history of the catalytic triad of NAT enzymes and the significance of the predominance of Cys-His-Asp catalytic triads over the non-canonical Cys-His-Glu catalytic triad, which can be catalytically active as shown in this study. A similar question has been proposed for the conservation of an Asp residue instead of the equivalent Glu residue in the catalytic triad of acetylcholine esterase (52). Finally, our work reinforces recent studies demonstrating greater structural and functional diversity than expected in the NAT enzyme family of acetyltransferases. The systematic functional and structural characterization of NATs harboring atypical features, in particular the homologues carrying the non-canonical Cys-His-Glu catalytic triad, will be necessary to better understand the basis of this emerging diversity.

**Acknowledgments**—We thank Dr. Karine Moncoq and Ximing Xu for helpful discussions on crystallogenesis assays. We acknowledge the core platform "Biophysique des macromolécules et de leurs intégrations" (Institut Pasteur) for provision of CD facilities. We also acknowledge Synchrotron Swiss Light Source (Villigen, Switzerland) for the provision of synchrotron radiation facilities, and we thank Meitian Wang for assistance with the use of the Beamline X06DA.

## REFERENCES

1. Hein, D. W. (2002) Molecular genetics and function of NAT1 and NAT2: role in aromatic amine metabolism and carcinogenesis. *Mutat. Res.* **506–507**, 65–77

2. Sim, E., Fakis, G., Laurieri, N., and Boukouvala, S. (2012) Arylamine N-acetyltransferases—from drug metabolism and pharmacogenetics to identification of novel targets for pharmacological intervention. *Adv. Pharmacol.* **63**, 169–205

3. Sinclair, J. C., Sandy, J., Delgoda, R., Sim, E., and Noble, M. E. (2000) Structure of arylamine N-acetyltransferase reveals a catalytic triad. *Nat. Struct. Biol.* **7**, 560–564

4. Rodrigues-Lima, F., Deloménie, C., Goodfellow, G. H., Grant, D. M., and Dupret, J. M. (2001) Homology modelling and structural analysis of human arylamine N-acetyltransferase NAT1: evidence for the conservation

of a cysteine protease catalytic domain and an active-site loop. *Biochem. J.* **356**, 327–334

5. Glenn, A. E., Karagianni, E. P., Ullndreaj, A., and Boukouvala, S. (2010) Comparative genomic and phylogenetic investigation of the xenobiotic metabolizing arylamine N-acetyltransferase enzyme family. *FEBS Lett.* **584**, 3158–3164

6. Riddle, B., and Jencks, W. P. (1971) Acetyl-coenzyme A:arylamine N-acetyltransferase. Role of the acetyl-enzyme intermediate and the effects of substituents on the rate. *J. Biol. Chem.* **246**, 3250–3258

7. Weber, W. W., and Hein, D. W. (1985) N-Acetylation pharmacogenetics. *Pharmacol. Rev.* **37**, 25–79

8. Dupret, J. M., and Grant, D. M. (1992) Site-directed mutagenesis of recombinant human arylamine N-acetyltransferase expressed in *Escherichia coli*. Evidence for direct involvement of Cys⁶⁸ in the catalytic mechanism of polymorphic human NAT2. *J. Biol. Chem.* **267**, 7381–7385

9. Sandy, J., Mushtaq, A., Holton, S. J., Schartau, P., Noble, M. E., and Sim, E. (2005) Investigation of the catalytic triad of arylamine N-acetyltransferases: essential residues required for acetyl transfer to arylamines. *Biochem. J.* **390**, 115–123

10. Zang, Y., Zhao, S., Doll, M. A., Christopher States, J., and Hein, D. W. (2007) Functional characterization of the A411T (L137F) and G364A (D122N) genetic polymorphisms in human N-acetyltransferase 2. *Pharmacogenet. Genomics* **17**, 37–45

11. Sim, E., Lack, N., Wang, C.-J., Long, H., Westwood, I., Fullam, E., and Kawamura, A. (2008) Arylamine N-acetyltransferases: structural and functional implications of polymorphisms. *Toxicology* **254**, 170–183

12. Kubiak, X., Dairou, J., Dupret, J.-M., and Rodrigues-Lima, F. (2013) Crystal structure of arylamine N-acetyltransferases: insights into the mechanisms of action and substrate selectivity. *Expert Opin. Drug Metab. Toxicol.* **9**, 349–362

13. Kubiak, X., Pluvinage, B., Li de la Sierra-Gallay, I., Weber, P., Haouz, A., Dupret, J.-M., and Rodrigues-Lima, F. (2012) Purification, crystallization and preliminary x-ray characterization of Bacillus cereus arylamine N-acetyltransferase 3 [(BACCR)NAT3]. *Acta Crystallogr. Sect. F Struct. Biol. Cryst. Commun.* **68**, 196–198

14. Kabsch, W. (2010) XDS. *Acta Crystallogr. D Biol. Crystallogr.* **66**, 125–132

15. McCoy, A. J., Grosse-Kunstleve, R. W., Adams, P. D., Winn, M. D., Stononi, L. C., and Read, R. J. (2007) Phaser crystallographic software. *J. Appl. Crystallogr.* **40**, 658–674

16. Winn, M. D., Ballard, C. C., Cowtan, K. D., Dodson, E. J., Emsley, P., Evans, P. R., Keegan, R. M., Krissinel, E. B., Leslie, A. G., McCoy, A., McNicholas, S. J., Murshudov, G. N., Pannu, N. S., Potterton, E. A., Powell, H. R., Read, R. J., Vagin, A., and Wilson, K. S. (2011) Overview of the CCP 4 suite and current developments. *Acta Crystallogr. D Biol. Crystallogr.* **67**, 235–242

17. Stein, N. (2008) CHAINSAW: a program for mutating pdb files used as templates in molecular replacement. *J. Appl. Crystallogr.* **41**, 641–643

18. Emsley, P., Lohkamp, B., Scott, W. G., and Cowtan, K. (2010) Features and development of Coot. *Acta Crystallogr. D Biol. Crystallogr.* **66**, 486–501

19. Murshudov, G. N., Skubåk, P., Lebedev, A. A., Pannu, N. S., Steiner, R. A., Nicholls, R. A., Winn, M. D., Long, F., and Vagin, A. A. (2011) REFMAC5 for the refinement of macromolecular crystal structures. *Acta Crystallogr. D Biol. Crystallogr.* **67**, 355–367

20. Bricogne, G., Blanc, E., Brnadl, M., Flensburg, C., Keller, P., Paciorek, W., Roversi, P., Sharff, A., Smart, O. S., Vonrhein, C., and Womack, T. O. (2011) BLISTER, Version 2.10.0, Global Phasing Ltd., Cambridge, UK

21. Chen, V. B., Arendall, W. B., 3rd, Headd, J. J., Keedy, D. A., Immormino, R. M., Kapral, G. J., Murray, L. W., Richardson, J. S., and Richardson, D. C. (2010) MolProbity: all-atom structure validation for macromolecular crystallography. *Acta Crystallogr. D Biol. Crystallogr.* **66**, 12–21

22. Laskowski, R. A., MacArthur, M. W., Moss, D. S., and Thornton, J. M. (1993) PROCHECK: a program to check the stereochemical quality of protein structures. *J. Appl. Crystallogr.* **26**, 283–291

23. DeLano, W. L. (2010) *The PyMOL Molecular Graphics System*, Version 0.99rc6, Schrödinger, LLC, New York

24. Brooke, E. W., Davies, S. G., Mulvaney, A. W., Pompeo, F., Sim, E., and Vickers, R. J. (2003) An approach to identifying novel substrates of bacterial arylamine N-acetyltransferases. *Bioorg. Med. Chem.* **11**, 1227–1234

25. Cleland, W. W., and Hengge, A. C. (1995) Mechanisms of phosphoryl and
acyl transfer. *FASEB J.* **9**, 1585–1594

26. King, S. M., and Johnson, W. C. (1999) Assigning secondary structure from protein coordinate data. *Proteins* **35**, 313–320

27. Le, H.-T., Chaffotte, A. F., Demey-Thomas, E., Vinh, J., Friguet, B., and Mary, J. (2009) Impact of hydrogen peroxide on the activity, structure, and conformational stability of the oxidized protein repair enzyme methionine sulfoxide reductase A. *J. Mol. Biol.* **393**, 58–66

28. Krekel, F., Samland, A. K., Macheroux, P., Amrhein, N., and Evans, J. N. (2000) Determination of the pKa value of C115 in MurA (UDP-N-acetylglucosamine enolpyruvyltransferase) from *Enterobacter cloacae*. *Biochemistry* **39**, 12671–12677

29. Tamura, K., Peterson, D., Peterson, N., Stecher, G., Nei, M., and Kumar, S. (2011) MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. *Mol. Biol. Evol.* **28**, 2731–2739

30. Pluvinage, B., Dairou, J., Possot, O. M., Martins, M., Fouet, A., Dupret, J.-M., and Rodrigues-Lima, F. (2007) Cloning and molecular characterization of three arylamine N-acetyltransferase genes from *Bacillus anthracis*: identification of unusual enzymatic properties and their contribution to sulfamethoxazole resistance. *Biochemistry* **46**, 7069–7078

31. Rodrigues-Lima, F., and Dupret, J.-M. (2002) In silico sequence analysis of arylamine N-acetyltransferases: evidence for an absence of lateral gene transfer from bacteria to vertebrates and first description of paralogs in bacteria. *Biochem. Biophys. Res. Commun.* **293**, 783–792

32. London, J., Skrzynia, C., and Goldberg, M. E. (1974) Renaturation of *Escherichia coli* tryptophanase after exposure to 8 M urea. Evidence for the existence of nucleation centers. *Eur. J. Biochem.* **47**, 409–415

33. Mitraki, A., Betton, J.-M., Desmadril, M., and Yon, J. M. (1987) Quasi-irreversibility in the unfolding-refolding transition of phosphoglycerate kinase induced by guanidine hydrochloride. *Eur. J. Biochem.* **163**, 29–34

34. Hammarström, P., Persson, M., Freskgård, P.-O., Mårtensson, L.-G., Andersson, D., Jonsson, B.-H., and Carlsson, U. (1993) Structural mapping of an aggregation nucleation site in a molten globule intermediate. *J. Biol. Chem.* **274**, 32897–32903

35. Xu, Q., Xie, Z., Ding, J., Lin, S.-X., and Xu, G. (2004) Monoclonal antibodies assisting refolding of firefly luciferase. *Protein Sci.* **13**, 1851–1858

36. Rajasekaran, M., Abirami, S., and Chen, C. (2011) Effects of single nucleotide polymorphisms on human N-acetyltransferase 2 structure and dynamics by molecular dynamics simulation. *PLoS One* **6**, e25801

37. Westwood, I. M., and Sim, E. (2007) Kinetic characterisation of arylamine N-acetyltransferase from *Pseudomonas aeruginosa*. *BMC Biochem.* **8**, 3

38. Sandy, J., Holton, S., Fullam, E., Sim, E., and Noble, M. (2005) Binding of the anti-tubercular drug isoniazid to the arylamine N-acetyltransferase protein from *Mycobacterium smegmatis*. *Protein Sci.* **14**, 775–782

39. Fullam, E., Westwood, I. M., Anderton, M. C., Lowe, E. D., Sim, E., and Noble, M. E. (2008) Divergence of cofactor recognition across evolution: coenzyme A binding in a prokaryotic arylamine N-acetyltransferase. *J. Mol. Biol.* **375**, 178–191

40. Wu, H., Dombrovsky, L., Tempel, W., Martin, F., Loppnau, P., Goodfel-

low, G. H., Grant, D. M., and Plotnikov, A. N. (2007) Structural basis of substrate-binding specificity of human arylamine N-acetyltransferases. *J. Biol. Chem.* **282**, 30189–30197

41. Pluvinage, B., Li de la Sierra-Gallay, I., Kubiak, X., Xu, X., Dairou, J., Dupret, J.-M., and Rodrigues-Lima, F. (2011) The *Bacillus anthracis* arylamine N-acetyltransferase ((BACAN)NAT1) that inactivates sulfamethoxazole, reveals unusual structural features compared with the other NAT isoenzymes. *FEBS Lett.* **585**, 3947–3952

42. Wang, H., Vath, G. M., Gleason, K. J., Hanna, P. E., and Wagner, C. R. (2004) Probing the mechanism of hamster arylamine N-acetyltransferase 2 acetylation by active site modification, site-directed mutagenesis, and pre-steady state and steady state kinetic studies. *Biochemistry* **43**, 8234–8246

43. Sikora, A. L., Frankel, B. A., and Blanchard, J. S. (2008) Kinetic and chemical mechanism of arylamine N-acetyltransferase from *Mycobacterium tuberculosis*. *Biochemistry* **47**, 10781–10789

44. Pillay, D., Boulangé, A. F., and Coetzer, T. H. (2010) Expression, purification and characterisation of two variant cysteine peptidases from *Trypanosoma congolense* with active site substitutions. *Protein Expr. Purif.* **74**, 264–271

45. Mehta, D. K., and Eckert, R. L. (2005) *Transglutaminases: Family of Enzymes with Diverse Functions*, Karger Publishers, Basel, Switzerland

46. Chu, X., He, H., Guo, C., and Sun, B. (2008) Identification of two novel esterases from a marine metagenomic library derived from South China Sea. *Appl. Microbiol. Biotechnol.* **80**, 615–625

47. Southworth, M. W., Benner, J., and Perler, F. B. (2000) An alternative protein splicing mechanism for inteins lacking an N-terminal nucleophile. *EMBO J.* **19**, 5019–5026

48. Saleh, L., Southworth, M. W., Considine, N., O'Neill, C., Benner, J., Bollinger, J. M., Jr., and Perler, F. B. (2011) Branched intermediate formation is the slowest step in the protein splicing reaction of the Ala1 KlbA intein from *Methanococcus jannaschii*. *Biochemistry* **50**, 10576–10589

49. Guerois, R., Nielsen, J. E., and Serrano, L. (2002) Predicting changes in the stability of proteins and protein complexes: a study of more than 1000 mutations. *J. Mol. Biol.* **320**, 369–387

50. Gonnelli, G., Rooman, M., and Dehouck, Y. (2012) Structure-based mutant stability predictions on proteins of unknown structure. *J. Biotechnol.* **161**, 287–293

51. Zhou, X., Zhang, N., Liu, L., Walters, K. J., Hanna, P. E., and Wagner, C. R. (2009) Probing the catalytic potential of the hamster arylamine N-acetyltransferase 2 catalytic triad by site-directed mutagenesis of the proximal conserved residue, Tyr190. *FEBS J.* **276**, 6928–6941

52. Soreq, H., and Seidman, S. (2001) Acetylcholinesterase—new roles for an old actor. *Nat. Rev. Neurosci.* **2**, 294–302

53. Dundas, J., Ouyang, Z., Tseng, J., Binkowski, A., Turpaz, Y., and Liang, J. (2006) CASTp: computed atlas of surface topography of proteins with structural and topographical mapping of functionally annotated residues. *Nucleic Acids Res.* **34**, W116–W118

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