# Involvement of arginine 878 together with Ca2+ in aminopeptidase A substrate specificity for N-terminal acidic amino-acid residues

Pierre Couvineau, Hugo de Almeida, Bernard Maigret, Catherine Llorens-Cortes, Xavier Iturrioz

► To cite this version:

Pierre Couvineau, Hugo de Almeida, Bernard Maigret, Catherine Llorens-Cortes, Xavier Iturrioz. Involvement of arginine 878 together with Ca2+ in aminopeptidase A substrate specificity for N-terminal acidic amino-acid residues. PLoS ONE, 2017, 12 (9), pp.e0184237. 10.1371/journal.pone.0184237. hal-01580832

HAL Id: hal-01580832

https://inria.hal.science/hal-01580832

Submitted on 10 Jan 2024

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_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Involvement of arginine 878 together with Ca2+ in aminopeptidase A substrate specificity for N-terminal acidic amino-acid residues

Pierre Couvineau, Hugo de Almeida, Bernard Maigret, Catherine Llorens-Cortes, Xavier Iturrioz

► To cite this version:

Pierre Couvineau, Hugo de Almeida, Bernard Maigret, Catherine Llorens-Cortes, Xavier Iturrioz. Involvement of arginine 878 together with Ca2+ in aminopeptidase A substrate specificity for N-terminal acidic amino-acid residues. PLoS ONE, 2017, 12 (9), pp.e0184237. 10.1371/journal.pone.0184237. hal-01580832

HAL Id: hal-01580832

https://inria.hal.science/hal-01580832v1

Submitted on 10 Jan 2024

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_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Involvement of arginine 878 together with Ca²⁺ in mouse aminopeptidase A substrate specificity for N-terminal acidic amino-acid residues

Pierre Couvineau¹,²,³, Hugo de Almeida¹,²,³, Bernard Maigret¹,²,³,⁴, Catherine Llorens-Cortes¹,²,³*, Xavier Iturrioz¹,²,³*

¹ Laboratory of central neuropeptides in the regulation of body fluid homeostasis and cardiovascular functions, INSERM U1050, Paris, France, ² Center for Interdisciplinary Research in Biology (CIRB), Collège de France, Paris, France, ³ CNRS, UMR 7241, Paris, France, ⁴ Laboratoire Lorrain de Recherche en Informatique et ses applications (LORIA), Vandoeuvre Les Nancy, France

‡ These authors are co-last authors on this work.

* c.llorens-cortes@college-de-france.fr (CLC); xavier.iturrioz@college-de-france.fr (XI)

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

## Abstract

Aminopeptidase A (APA) is a membrane-bound zinc metalloprotease cleaving, in the brain, the N-terminal aspartyl residue of angiotensin II to generate angiotensin III, which exerts a tonic stimulatory effect on the control of blood pressure in hypertensive animals. Using a refined APA structure derived from the human APA crystal structure, we docked the specific and selective APA inhibitor, EC33 in the presence of Ca²⁺. We report the presence in the S1 subsite of Arg-887 (Arg-878 in mouse APA), the guanidinium moiety of which established an interaction with the electronegative sulfonate group of EC33. Mutagenic replacement of Arg-878 with an alanine or a lysine residue decreased the affinity of the recombinant enzymes for the acidic substrate, α-L-glutamyl-β-naphthylamide, with a slight decrease in substrate hydrolysis velocity either with or without Ca²⁺. In the absence of Ca²⁺, the mutations modified the substrate specificity of APA for the acidic substrate, the mutated enzymes hydrolyzing more efficiently basic and neutral substrates, although the addition of Ca²⁺ partially restored the acidic substrate specificity. The analysis of the 3D models of the Arg-878 mutated APAs revealed a change in the volume of the S1 subsite, which may impair the binding and/or the optimal positioning of the substrate in the active site as well as its hydrolysis. These findings demonstrate the key role of Arg-878 together with Ca²⁺ in APA substrate specificity for N-terminal acidic amino acid residues by ensuring the optimal positioning of acidic substrates during catalysis.

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

## OPEN ACCESS

**Citation:** Couvineau P, de Almeida H, Maigret B, Llorens-Cortes C, Iturrioz X (2017) Involvement of arginine 878 together with Ca²⁺ in mouse aminopeptidase A substrate specificity for N-terminal acidic amino-acid residues. PLoS ONE 12 (9): e0184237. https://doi.org/10.1371/journal.pone.0184237

**Editor:** Eugene A. Permyakov, Russian Academy of Medical Sciences, RUSSIAN FEDERATION

**Received:** May 22, 2017

**Accepted:** August 21, 2017

**Published:** September 6, 2017

**Copyright:** © 2017 Couvineau et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This work was supported by the Institut National de la Santé et de la Recherche Médicale, the Centre National de la Recherche Scientifique, the Collège de France. P. C. and H. D. A. were supported by fellowships from the Agence Nationale pour la Recherche RPIB and LabCom. P.C. has been also supported by a fellowship from

## Introduction

Aminopeptidase A (APA; EC 3.4.11.7) is a 160 kDa homodimeric type II Zn²⁺ membrane-bound aminopeptidase [1, 2]. APA cleaves the N-terminal glutamyl or aspartyl residue from peptide substrates such as angiotensin II (AngII) and cholecystokinin-8 and is activated by Ca²⁺
RESEARCH ARTICLE

# Involvement of arginine 878 together with Ca²⁺ in mouse aminopeptidase A substrate specificity for N-terminal acidic amino-acid residues

Pierre Couvineau¹,²,³, Hugo de Almeida¹,²,³, Bernard Maigret¹,²,³,⁴, Catherine Llorens-Cortes¹,²,³*, Xavier Iturrioz¹,²,³*

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

¹ Laboratory of central neuropeptides in the regulation of body fluid homeostasis and cardiovascular functions, INSERM U1050, Paris, France, ² Center for Interdisciplinary Research in Biology (CIRB), Collège de France, Paris, France, ³ CNRS, UMR 7241, Paris, France, ⁴ Laboratoire Lorrain de Recherche en Informatique et ses applications (LORIA), Vandoeuvre Les Nancy, France

‡ These authors are co-last authors on this work.

* c.llorens-cortes@college-de-france.fr (CLC); xavier.iturrioz@college-de-france.fr (XI)

## Abstract

Aminopeptidase A (APA) is a membrane-bound zinc metalloprotease cleaving, in the brain, the N-terminal aspartyl residue of angiotensin II to generate angiotensin III, which exerts a tonic stimulatory effect on the control of blood pressure in hypertensive animals. Using a refined APA structure derived from the human APA crystal structure, we docked the specific and selective APA inhibitor, EC33 in the presence of Ca²⁺. We report the presence in the S1 subsite of Arg-887 (Arg-878 in mouse APA), the guanidinium moiety of which established an interaction with the electronegative sulfonate group of EC33. Mutagenic replacement of Arg-878 with an alanine or a lysine residue decreased the affinity of the recombinant enzymes for the acidic substrate, α-L-glutamyl-β-naphthylamide, with a slight decrease in substrate hydrolysis velocity either with or without Ca²⁺. In the absence of Ca²⁺, the mutations modified the substrate specificity of APA for the acidic substrate, the mutated enzymes hydrolyzing more efficiently basic and neutral substrates, although the addition of Ca²⁺ partially restored the acidic substrate specificity. The analysis of the 3D models of the Arg-878 mutated APAs revealed a change in the volume of the S1 subsite, which may impair the binding and/or the optimal positioning of the substrate in the active site as well as its hydrolysis. These findings demonstrate the key role of Arg-878 together with Ca²⁺ in APA substrate specificity for N-terminal acidic amino acid residues by ensuring the optimal positioning of acidic substrates during catalysis.

## Introduction

Aminopeptidase A (APA; EC 3.4.11.7) is a 160 kDa homodimeric type II Zn²⁺ membrane-bound aminopeptidase [1, 2]. APA cleaves the N-terminal glutamyl or aspartyl residue from peptide substrates such as angiotensin II (AngII) and cholecystokinin-8 and is activated by Ca²⁺

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

### OPEN ACCESS

**Citation:** Couvineau P, de Almeida H, Maigret B, Llorens-Cortes C, Iturrioz X (2017) Involvement of arginine 878 together with Ca²⁺ in mouse aminopeptidase A substrate specificity for N-terminal acidic amino-acid residues. PLoS ONE 12 (9): e0184237. https://doi.org/10.1371/journal.pone.0184237

**Editor:** Eugene A. Permyakov, Russian Academy of Medical Sciences, RUSSIAN FEDERATION

**Received:** May 22, 2017

**Accepted:** August 21, 2017

**Published:** September 6, 2017

**Copyright:** © 2017 Couvineau et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This work was supported by the Institut National de la Santé et de la Recherche Médicale, the Centre National de la Recherche Scientifique, the Collège de France. P. C. and H. D. A. were supported by fellowships from the Agence Nationale pour la Recherche RPIB and LabCom. P.C. has been also supported by a fellowship from
Labex MemoLife (Paris Sciences et Lettres). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared that no competing interests exist.

[3, 4]. Ca²⁺ not only enhances the hydrolysis by APA of N-terminal acidic amino acid residues from substrates, but also decreases the hydrolysis of N-terminal neutral or basic residues [5]. APA is expressed in various tissues, including the brush border of intestinal and renal epithelial cells and the vascular endothelium [6]. This enzyme is also expressed in several brain nuclei involved in the control of body fluid homeostasis and cardiovascular function, together with other components of the brain renin-angiotensin system [7]. Studies with the specific and selective APA inhibitor EC33 [(S)-3-amino-4-mercapto-butyl sulfonic acid] [8] have shown that, in the brain, APA generates angiotensin III (AngIII) [9] which exerts a tonic stimulatory action on the control of blood pressure in hypertensive animals and increases arginine vasopressin release [9–11]. Thus, the inhibition of brain APA with specific and selective inhibitors normalizes blood pressure in conscious hypertensive rats [10–13]. These studies suggest that brain APA constitutes a promising target for the treatment of hypertension [14–16]. This justifies the development of potent and selective APA inhibitors crossing the blood-brain barrier after oral administration for use as centrally-acting antihypertensive agents. With this aim, we have pursued our studies on the APA active site organization by molecular modeling and site-directed mutagenesis. First, using the crystal structure of leukotriene-A4 hydrolase (LTA4H; EC 3.3.2.6) [17] as a template and the functional information collected from previous site-directed mutagenesis studies on APA [18–24], we built a three-dimensional (3D) model of the mouse APA (mAPA) ectodomain from residue 79 to 559 [16, 24], including the catalytic domain. We subsequently docked into the APA active site, specific and selective APA inhibitors, such as Glu-PO₃H₂ (4-amino-4-phosphonobutyric acid) a transition state analog [25].

In the 3D model of APA, the Zn²⁺ atom is coordinated by the two histidine residues (His385 and His389) of the HEXXH motif, a water molecule, glutamate 408 (Glu408) of the WLNEG motif and either by one of the oxygen atoms of the phosphate of GluPO₃H₂ [24]. The water molecule bound to the Zn²⁺atom also interacts with Glu-352 and Glu-386 side chains [19, 21]. Ca²⁺ was then introduced in the APA 3D model. The Ca²⁺ atom was localized at the bottom of the S1 subsite, as suggested previously by Danielsen *et al*. [1]. The Ca²⁺ atom interacts with the acidic side chains of Asp-213 and Asp-218, the carbonyl group of Glu-215 and three water molecules, one of them being engaged in a hydrogen bond with the negatively charged carboxylate side chains of the inhibitors [26]. Site-directed mutagenesis of Asp-213 and Asp-218 (Asp-221 in human APA (hAPA)) showed that their negative charges are required for the presence of the Ca²⁺ atom that ensures acidic APA substrate specificity [26, 27]. More recently, the X-ray structure of human APA (residues 76 to 956) was resolved [28] and its comparison with the APA 3D model showed a perfect overlap of the catalytic regions.

In the present study, we have therefore taken advantage of the X-ray-derived structure of human APA to build a new system comprising the complete refined APA 3D structure complexed with EC33, an APA inhibitor exhibiting an inhibitory potency of 10⁻⁷ M on APA and 100-fold less active on APN [8]. We visualized, in the S1 subsite, Arg-878 (Arg-887 in hAPA), which interacts with the electronegative sulfonate moiety of EC33 as observed by Yang *et al* in the crystal structure of human APA bound to a Glu, exhibiting an inhibitory potency of 7 x 10⁻³ M on APA [29]. The first data obtained by these authors on the role of Arg-887 in the absence of Ca²⁺ suggested an involvement in hAPA substrate specificity [28]. However, the Ca²⁺ increases APA preference for acidic substrates and this Ca²⁺-modulated APA substrate specificity is deemed physiologically relevant since the concentrations of Ca²⁺ that modulate APA activity, are in the same range as those found in brain fluid (i.e. 1–2 mM) [30]. Taking into account that Ca²⁺ plays a major role in APA substrate specificity, the aim of our work was to deepen the role of Arg-878 of mAPA together with Ca²⁺ in substrate/inhibitor binding and substrate specificity of APA for N-terminal acidic amino acid residues. For this purpose, we replaced Arg-878 with an alanine or a lysine residue and checked that the mutated enzymes
Labex MemoLife (Paris Sciences et Lettres). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared that no competing interests exist.

[3, 4]. Ca²⁺ not only enhances the hydrolysis by APA of N-terminal acidic amino acid residues from substrates, but also decreases the hydrolysis of N-terminal neutral or basic residues [5]. APA is expressed in various tissues, including the brush border of intestinal and renal epithelial cells and the vascular endothelium [6]. This enzyme is also expressed in several brain nuclei involved in the control of body fluid homeostasis and cardiovascular function, together with other components of the brain renin-angiotensin system [7]. Studies with the specific and selective APA inhibitor EC33 [(S)-3-amino-4-mercapto-butyl sulfonic acid] [8] have shown that, in the brain, APA generates angiotensin III (AngIII) [9] which exerts a tonic stimulatory action on the control of blood pressure in hypertensive animals and increases arginine vasopressin release [9–11]. Thus, the inhibition of brain APA with specific and selective inhibitors normalizes blood pressure in conscious hypertensive rats [10–13]. These studies suggest that brain APA constitutes a promising target for the treatment of hypertension [14–16]. This justifies the development of potent and selective APA inhibitors crossing the blood-brain barrier after oral administration for use as centrally-acting antihypertensive agents. With this aim, we have pursued our studies on the APA active site organization by molecular modeling and site-directed mutagenesis. First, using the crystal structure of leukotriene-A4 hydrolase (LTA4H; EC 3.3.2.6) [17] as a template and the functional information collected from previous site-directed mutagenesis studies on APA [18–24], we built a three-dimensional (3D) model of the mouse APA (mAPA) ectodomain from residue 79 to 559 [16, 24], including the catalytic domain. We subsequently docked into the APA active site, specific and selective APA inhibitors, such as Glu-PO₃H₂ (4-amino-4-phosphonobutyric acid) a transition state analog [25].

In the 3D model of APA, the Zn²⁺ atom is coordinated by the two histidine residues (His385 and His389) of the HEXXH motif, a water molecule, glutamate 408 (Glu408) of the WLNEG motif and either by one of the oxygen atoms of the phosphate of GluPO₃H₂ [24]. The water molecule bound to the Zn²⁺atom also interacts with Glu-352 and Glu-386 side chains [19, 21]. Ca²⁺ was then introduced in the APA 3D model. The Ca²⁺ atom was localized at the bottom of the S1 subsite, as suggested previously by Danielsen *et al*. [1]. The Ca²⁺ atom interacts with the acidic side chains of Asp-213 and Asp-218, the carbonyl group of Glu-215 and three water molecules, one of them being engaged in a hydrogen bond with the negatively charged carboxylate side chains of the inhibitors [26]. Site-directed mutagenesis of Asp-213 and Asp-218 (Asp-221 in human APA (hAPA)) showed that their negative charges are required for the presence of the Ca²⁺ atom that ensures acidic APA substrate specificity [26, 27]. More recently, the X-ray structure of human APA (residues 76 to 956) was resolved [28] and its comparison with the APA 3D model showed a perfect overlap of the catalytic regions.

In the present study, we have therefore taken advantage of the X-ray-derived structure of human APA to build a new system comprising the complete refined APA 3D structure complexed with EC33, an APA inhibitor exhibiting an inhibitory potency of 10⁻⁷ M on APA and 100-fold less active on APN [8]. We visualized, in the S1 subsite, Arg-878 (Arg-887 in hAPA), which interacts with the electronegative sulfonate moiety of EC33 as observed by Yang *et al* in the crystal structure of human APA bound to a Glu, exhibiting an inhibitory potency of 7 x 10⁻³ M on APA [29]. The first data obtained by these authors on the role of Arg-887 in the absence of Ca²⁺ suggested an involvement in hAPA substrate specificity [28]. However, the Ca²⁺ increases APA preference for acidic substrates and this Ca²⁺-modulated APA substrate specificity is deemed physiologically relevant since the concentrations of Ca²⁺ that modulate APA activity, are in the same range as those found in brain fluid (i.e. 1–2 mM) [30]. Taking into account that Ca²⁺ plays a major role in APA substrate specificity, the aim of our work was to deepen the role of Arg-878 of mAPA together with Ca²⁺ in substrate/inhibitor binding and substrate specificity of APA for N-terminal acidic amino acid residues. For this purpose, we replaced Arg-878 with an alanine or a lysine residue and checked that the mutated enzymes
displayed similar processing and subcellular distribution to wild-type mAPA. We then bio-
chemically and kinetically characterized the purified recombinant wild-type and mutated
enzymes with various substrates and determined their sensitivity to Ca²⁺ and various inhibi-
tors exhibiting different side chains targeting the S1 subsite of APA.

**Materials and methods**

Materials

Restriction endonucleases Dpn1 was obtained from New England Biolabs Inc (Evry, France) and was used according to the manufacturer's instructions. The PfuUltra High-fidelity DNA Polymerase was purchased from Agilent (les Ulis, France). The liposomal transfection reagent, Lipofectamine 2000, the pcDNA 3.1-His vector and the monoclonal anti-Xpress antibody were purchased from Life Technologies (Cergy-Pontoise, France). CHO-K1 cells were obtained from American Type Culture Collection (Manassas, VA, USA). The horseradish peroxidase-conjugated sheep anti-mouse antibody was purchased from Sigma-Aldrich (Saint Quentin Fal-lavier, France). The cOmplete, EDTA-free Protease Inhibitor Cocktail was purchased from Roche (Mannheim, Germany). Immobilized cobalt affinity columns (Talon) were obtained from Clontech (Heidelberg, Germany). The synthetic substrates, GluNA, AspNA, AlaNA and LysNA were purchased from Bachem (Bunderdorf, Switzerland).

Molecular docking and molecular dynamics simulations

The Crystallographic structure of human aminopeptidase A complexed with glutamate and calcium (PDB ID: 4KXD) was used to perform docking and molecular dynamics simulations. The 3D-structure was first treated to remove the Glu ligand, as well as the N-Acetyl-D-Glucosamine chains and non-conserved water molecules. Zn²⁺ and Ca²⁺ atoms were kept, and a loop was added for residues missing from crystallographic data (residues 608 to 611) using the software MODELLER [31]. Reference inhibitors developed by our group were docked into our 3D-model using GOLD v5.1 [32], with the ChemPLP scoring function selected and several penalties added based on the analysis of available crystallographic structures. The list of docked molecules and the applied penalties can be found in S3 and S4 Tables, respectively. In order to analyze if the predicted interactions would be conserved in a physiological and dynamical molecular system, compounds EC33 and GluPO₃H₂ were chosen to be parametrized, and poses obtained from the automated docking process above were submitted to molecular dynamics (MD) simulations using NAMD. Briefly, small ligands were parametrized with Discovery Studio Visualizer 4.5 (BIOVIA) using CHARMM 36 forcefields [33, 34], with atom types being assigned manually. Ligand-protein complexes were generated with VMD autopsf plug-in using the BIOVIA provided topology files. The MUTATOR plugin was used to build systems with mutated hAPAs (R887A and R887K). The obtained systems were then included in a water-box using the VMD solvate plug-in, followed by charge neutralization through the ionize plug-in. The final systems comprising ~80,000 atoms were then subjected to three stages of conjugate gradient energy minimization of 6400 steps each. The first stage had the ligand-protein complex fixed, and allowed only water to move freely. The second stage allowed ligand motions, and the last stage consisted of the energy minimization of the entire system. Those systems were then subjected to a short MD simulation of 500 ps in order to equilibrate pressure and the temperatures before proceeding to production runs. Each MD production run simulation was carried out in the NPT ensemble, with periodic boundary conditions, Particle Mesh Ewald for treating electrostatic interactions and 1 fs per time step for 50,000,000 steps, generating a trajectory of 50 ns. Trajectory analyses were performed with VMD, pair interaction calculations with NAMD and ligand binding site volume with MDPOCKET.
displayed similar processing and subcellular distribution to wild-type mAPA. We then bio-
chemically and kinetically characterized the purified recombinant wild-type and mutated
enzymes with various substrates and determined their sensitivity to Ca²⁺ and various inhibi-
tors exhibiting different side chains targeting the S1 subsite of APA.

**Materials and methods**

Materials

Restriction endonucleases Dpn1 was obtained from New England Biolabs Inc (Evry, France) and was used according to the manufacturer's instructions. The PfuUltra High-fidelity DNA Polymerase was purchased from Agilent (les Ulis, France). The liposomal transfection reagent, Lipofectamine 2000, the pcDNA 3.1-His vector and the monoclonal anti-Xpress antibody were purchased from Life Technologies (Cergy-Pontoise, France). CHO-K1 cells were obtained from American Type Culture Collection (Manassas, VA, USA). The horseradish peroxidase-conjugated sheep anti-mouse antibody was purchased from Sigma-Aldrich (Saint Quentin Fal-lavier, France). The cOmplete, EDTA-free Protease Inhibitor Cocktail was purchased from Roche (Mannheim, Germany). Immobilized cobalt affinity columns (Talon) were obtained from Clontech (Heidelberg, Germany). The synthetic substrates, GluNA, AspNA, AlaNA and LysNA were purchased from Bachem (Bunderdorf, Switzerland).

Molecular docking and molecular dynamics simulations

The Crystallographic structure of human aminopeptidase A complexed with glutamate and calcium (PDB ID: 4KXD) was used to perform docking and molecular dynamics simulations. The 3D-structure was first treated to remove the Glu ligand, as well as the N-Acetyl-D-Glucosamine chains and non-conserved water molecules. Zn²⁺ and Ca²⁺ atoms were kept, and a loop was added for residues missing from crystallographic data (residues 608 to 611) using the software MODELLER [31]. Reference inhibitors developed by our group were docked into our 3D-model using GOLD v5.1 [32], with the ChemPLP scoring function selected and several penalties added based on the analysis of available crystallographic structures. The list of docked molecules and the applied penalties can be found in S3 and S4 Tables, respectively. In order to analyze if the predicted interactions would be conserved in a physiological and dynamical molecular system, compounds EC33 and GluPO₃H₂ were chosen to be parametrized, and poses obtained from the automated docking process above were submitted to molecular dynamics (MD) simulations using NAMD. Briefly, small ligands were parametrized with Discovery Studio Visualizer 4.5 (BIOVIA) using CHARMM 36 forcefields [33, 34], with atom types being assigned manually. Ligand-protein complexes were generated with VMD autopsf plug-in using the BIOVIA provided topology files. The MUTATOR plugin was used to build systems with mutated hAPAs (R887A and R887K). The obtained systems were then included in a water-box using the VMD solvate plug-in, followed by charge neutralization through the ionize plug-in. The final systems comprising ~80,000 atoms were then subjected to three stages of conjugate gradient energy minimization of 6400 steps each. The first stage had the ligand-protein complex fixed, and allowed only water to move freely. The second stage allowed ligand motions, and the last stage consisted of the energy minimization of the entire system. Those systems were then subjected to a short MD simulation of 500 ps in order to equilibrate pressure and the temperatures before proceeding to production runs. Each MD production run simulation was carried out in the NPT ensemble, with periodic boundary conditions, Particle Mesh Ewald for treating electrostatic interactions and 1 fs per time step for 50,000,000 steps, generating a trajectory of 50 ns. Trajectory analyses were performed with VMD, pair interaction calculations with NAMD and ligand binding site volume with MDPOCKET.
## Cloning and site-directed mutagenesis

The cDNA encoding mouse APA [22] was used as a template for the generation of mutants by PCR-based site-directed mutagenesis. The arginine residue at position 878 was replaced with either an alanine or a lysine residue. The mutated primers A (forward) and B (reverse) were:

A (5'-CAATCAATGACGCATACCTTGGCCGGATCG-3') and B (5'-GGCCAAGGTATGCGTCATTGATTGTAAATC-3') for R878A

A (5'-CAATCAATGACAAATACCTTGGCCGGATCG-3') and B (5'-GGCCAAGGTATTTGTCATTGATTGTAAATC-3') for R878K

The underlined bases encode the new amino acid residue and replace the naturally occurring codon at position 2634. Nucleotide numbering corresponds to the mAPA sequence deposited in the GenBank/EBI database (accession number M29961). PCR was performed with PfuUltra High-fidelity DNA Polymerase to make a full-length plasmid PCR product. The methylated plasmid used as the template was then specifically digested with DpnI and Top 10 E. coli were transformed to obtain the mutated cDNA. Sequences of all constructs were checked by DNA Sequencing (GATC Biotech, Mulhouse, France).

## Cell culture and establishment of stable CHO-K1 cell lines producing wild-type and mutated His-mAPAs

Chinese Hamster Ovarian (CHO—American Type Culture Collection) cells were cultured in Ham's F-12 medium supplemented with 10% fetal calf serum, 0.5 mM glutamine, 100 units/ml penicillin and 100 µg/ml streptomycin (all from Life Technologies) at 37°C in a humidified atmosphere of 5% CO₂, 95% air. The cells were transfected with either 2.5 µg of plasmid containing the Xpress-polyHis-tagged wild-type or the Xpress-polyHis-tagged mutated mAPA cDNAs using Lipofectamine 2000 following manufacturer instructions. Stable cell lines producing the wild-type and mutated His-APAs were established as previously described [23].

## Immunofluorescence analysis of stably transfected CHO cells

Stable transfected CHO cells producing wild-type and mutated His-mAPAs were seeded (100,000 cells) onto 14-mm-diameter coverslips and grown for 24h. Then, the cells were fixed and permeabilized in 100% ice-cold methanol for 5 min on ice. After blocking with 0.5% BSA in 0.1 M PBS, coverslips were incubated with the anti-Xpress antibody (1:1000) in 0.5% BSA in PBS for 2h at 20°C, washed thrice with cold PBS and then incubated with a 1:500 dilution of AlexaFluor 488-conjugated polyclonal anti-rabbit antibody in 0.5% BSA in PBS for 2 h at 20°C. The coverslips were then mounted in Mowiol. The cells were examined with a Leica TCS SP V (Leica Microsystem, Heidelberg, Germany) confocal laser scanning microscope equipped with an argon/krypton laser and configured with a Leica DM IRBE inverted microscope. AlexaFluor was detected after 30% excitation at 488 nm. The images (1024 x1024 pixels) were obtained with a 63x magnification oil-immersion objective. Each image corresponded to a cross-section of the cells.

## Purification and western blot analysis of wild-type and mutated recombinant His-mAPAs

Stable transfected CHO cells (3x10⁹ cells) were harvested and a crude membrane preparation was obtained, as previously described [23]. Crude membrane preparation was solubilized overnight in 50 ml of 20 mM Tris-HCL buffer (pH 8) containing 100 mM NaCl and 1% nonidet P-40 (loading buffer, 1 ml/10⁶ cells ratio). The solubilized membrane preparation was centrifuged at 100,000 x g for 90 min and the supernatant was subjected to metal affinity
## Cloning and site-directed mutagenesis

The cDNA encoding mouse APA [22] was used as a template for the generation of mutants by PCR-based site-directed mutagenesis. The arginine residue at position 878 was replaced with either an alanine or a lysine residue. The mutated primers A (forward) and B (reverse) were:

A (5'-CAATCAATGACGCATACCTTGGCCGGATCG-3') and B (5'-GGCCAAGGTATGCGTCATTGATTGTAAATC-3') for R878A

A (5'-CAATCAATGACAAATACCTTGGCCGGATCG-3') and B (5'-GGCCAAGGTATTTGTCATTGATTGTAAATC-3') for R878K

The underlined bases encode the new amino acid residue and replace the naturally occurring codon at position 2634. Nucleotide numbering corresponds to the mAPA sequence deposited in the GenBank/EBI database (accession number M29961). PCR was performed with PfuUltra High-fidelity DNA Polymerase to make a full-length plasmid PCR product. The methylated plasmid used as the template was then specifically digested with DpnI and Top 10 E. coli were transformed to obtain the mutated cDNA. Sequences of all constructs were checked by DNA Sequencing (GATC Biotech, Mulhouse, France).

## Cell culture and establishment of stable CHO-K1 cell lines producing wild-type and mutated His-mAPAs

Chinese Hamster Ovarian (CHO—American Type Culture Collection) cells were cultured in Ham's F-12 medium supplemented with 10% fetal calf serum, 0.5 mM glutamine, 100 units/ml penicillin and 100 µg/ml streptomycin (all from Life Technologies) at 37°C in a humidified atmosphere of 5% CO₂, 95% air. The cells were transfected with either 2.5 µg of plasmid containing the Xpress-polyHis-tagged wild-type or the Xpress-polyHis-tagged mutated mAPA cDNAs using Lipofectamine 2000 following manufacturer instructions. Stable cell lines producing the wild-type and mutated His-APAs were established as previously described [23].

## Immunofluorescence analysis of stably transfected CHO cells

Stable transfected CHO cells producing wild-type and mutated His-mAPAs were seeded (100,000 cells) onto 14-mm-diameter coverslips and grown for 24h. Then, the cells were fixed and permeabilized in 100% ice-cold methanol for 5 min on ice. After blocking with 0.5% BSA in 0.1 M PBS, coverslips were incubated with the anti-Xpress antibody (1:1000) in 0.5% BSA in PBS for 2h at 20°C, washed thrice with cold PBS and then incubated with a 1:500 dilution of AlexaFluor 488-conjugated polyclonal anti-rabbit antibody in 0.5% BSA in PBS for 2 h at 20°C. The coverslips were then mounted in Mowiol. The cells were examined with a Leica TCS SP V (Leica Microsystem, Heidelberg, Germany) confocal laser scanning microscope equipped with an argon/krypton laser and configured with a Leica DM IRBE inverted microscope. AlexaFluor was detected after 30% excitation at 488 nm. The images (1024 x1024 pixels) were obtained with a 63x magnification oil-immersion objective. Each image corresponded to a cross-section of the cells.

## Purification and western blot analysis of wild-type and mutated recombinant His-mAPAs

Stable transfected CHO cells (3x10⁹ cells) were harvested and a crude membrane preparation was obtained, as previously described [23]. Crude membrane preparation was solubilized overnight in 50 ml of 20 mM Tris-HCL buffer (pH 8) containing 100 mM NaCl and 1% nonidet P-40 (loading buffer, 1 ml/10⁶ cells ratio). The solubilized membrane preparation was centrifuged at 100,000 x g for 90 min and the supernatant was subjected to metal affinity
chromatography with a metal chelate resin column (Talon Co²⁺, Clontech). The wild-type or mutated His-mAPAs were eluted with loading buffer containing 250 mM imidazole. The eluted fractions were loaded into Pierce Protein Concentrators, 150K MWCO (ThermoFisher) to concentrate the fraction and to replace the elution buffer with 20 mM Tris-HCl buffer (pH 7.4) containing 100 mM NaCl and 0.1% nonidet P-40. Protein concentrations were determined with the BCA assay kit (Calbiochem). The purity of the final preparation was assessed by NuPAGE (4–12% polyacrylamide gels, Invitrogen) and SilverQuest Silver staining (all from Life Technologies). Purified wild-type and mutated His-mAPAs were resolved by 4–12% SDS-PAGE and the proteins were transferred on a nitrocellulose membrane. His-tagged recombinant proteins were detected with the anti-Xpress antibody (1:5000 dilution). The resulting immune complex was detected with an anti-mouse antibody coupled to HRP (GE Healthcare, 1:10,000) and revealed with the ECL Plus Western blotting detection system (GE Healthcare) and image was acquired with an ImageQuant LAS 4000 (GE Healthcare).

## Enzyme assay

Enzyme assays were carried out for wild-type and mutated His-mAPAs in 50 mM Tris-HCl buffer, pH 7.4, with or without 4 mM CaCl₂, by monitoring the rate of hydrolysis of GluNA, as previously described [19]. The sensitivity to Ca²⁺ of purified wild-type and mutated His-mAPA was determined by incubating 0.5 mM GluNA with 0 to 4 mM CaCl₂ in a final volume of 100 µl of 50 mM Tris-HCl buffer, pH 7.4. All the assays were performed in black 96-well plates (solid black 96-well flat-bottomed non-treated plates with no lid; Corning Costar). All of the enzymatic studies were performed such that substrate hydrolysis rates were maintained below 10% (initial rate conditions). The kinetic curves were monitored by following the increase in fluorescence at 460 nm ($l_{ex}$ = 330 nm), induced by the release of the β-naphthylamide fluorogenic part of the substrate by APA. Fluorescence signals were monitored by counting photons with a spectrophotometer (Fusion universal microplate analyser; Packard) equipped with a temperature control device and a plate shaker. The kinetic parameters ($K_m$ and $k_{cat}$) were determined from Michaelis-Menten plots, using GraphPad Prism 4 software (GraphPad) with a final concentration of substrate (GluNA) from 0.005 to 5 mM with or without 4 mM Ca²⁺. For evaluation of the substrate specificity of APA with or without 4 mM Ca²⁺, the rate of hydrolysis of two additional substrates, α-L-aspartate-β-naphthylamide (AspNA), α-L-alanine-β-naphthylamide (AlaNA) and basic α-L-lysine-β-naphthylamide (LysNA), was determined under the same experimental conditions as described for GluNA. The sensitivity of wild-type and mutated His-mAPAs to inhibition by [(3S)-3-amino-4-sulfanyl-butane-1-sulfonic acid] (EC33), glutamate-thiol (GluSH), glutamate-phosphonic-acid (GluPO₃H₂), methionine-thiol (MetSH) and lysine-thiol (LysSH) was determined by establishing dose-dependent inhibition curves for a final GluNA concentration of 0.5 mM in the presence or absence of 4 mM CaCl₂. These compounds being linear competitive inhibitors, their $K_i$ were calculated from the formula $K_i = IC_{50}/(1 + [GluNA]/K_m)$.

## Statistical analysis

Statistical comparisons were carried out with Student's unpaired t-test. Differences were considered significant if *p* was less than 0.05.

## Results and discussion

APA belongs to the monozinc aminopeptidase family, which requires Zn²⁺ for catalytic activity. In addition, the substrate specificity of APA, unlike that of other monozinc aminopeptidases, is enhanced by Ca²⁺. Indeed, APA preferentially hydrolyzes peptides with N-terminal
chromatography with a metal chelate resin column (Talon Co²⁺, Clontech). The wild-type or mutated His-mAPAs were eluted with loading buffer containing 250 mM imidazole. The eluted fractions were loaded into Pierce Protein Concentrators, 150K MWCO (ThermoFisher) to concentrate the fraction and to replace the elution buffer with 20 mM Tris-HCl buffer (pH 7.4) containing 100 mM NaCl and 0.1% nonidet P-40. Protein concentrations were determined with the BCA assay kit (Calbiochem). The purity of the final preparation was assessed by NuPAGE (4–12% polyacrylamide gels, Invitrogen) and SilverQuest Silver staining (all from Life Technologies). Purified wild-type and mutated His-mAPAs were resolved by 4–12% SDS-PAGE and the proteins were transferred on a nitrocellulose membrane. His-tagged recombinant proteins were detected with the anti-Xpress antibody (1:5000 dilution). The resulting immune complex was detected with an anti-mouse antibody coupled to HRP (GE Healthcare, 1:10,000) and revealed with the ECL Plus Western blotting detection system (GE Healthcare) and image was acquired with an ImageQuant LAS 4000 (GE Healthcare).

## Enzyme assay

Enzyme assays were carried out for wild-type and mutated His-mAPAs in 50 mM Tris-HCl buffer, pH 7.4, with or without 4 mM CaCl₂, by monitoring the rate of hydrolysis of GluNA, as previously described [19]. The sensitivity to Ca²⁺ of purified wild-type and mutated His-mAPA was determined by incubating 0.5 mM GluNA with 0 to 4 mM CaCl₂ in a final volume of 100 µl of 50 mM Tris-HCl buffer, pH 7.4. All the assays were performed in black 96-well plates (solid black 96-well flat-bottomed non-treated plates with no lid; Corning Costar). All of the enzymatic studies were performed such that substrate hydrolysis rates were maintained below 10% (initial rate conditions). The kinetic curves were monitored by following the increase in fluorescence at 460 nm ($l_{ex}$ = 330 nm), induced by the release of the β-naphthylamide fluorogenic part of the substrate by APA. Fluorescence signals were monitored by counting photons with a spectrophotometer (Fusion universal microplate analyser; Packard) equipped with a temperature control device and a plate shaker. The kinetic parameters ($K_m$ and $k_{cat}$) were determined from Michaelis-Menten plots, using GraphPad Prism 4 software (GraphPad) with a final concentration of substrate (GluNA) from 0.005 to 5 mM with or without 4 mM Ca²⁺. For evaluation of the substrate specificity of APA with or without 4 mM Ca²⁺, the rate of hydrolysis of two additional substrates, α-L-aspartate-β-naphthylamide (AspNA), α-L-alanine-β-naphthylamide (AlaNA) and basic α-L-lysine-β-naphthylamide (LysNA), was determined under the same experimental conditions as described for GluNA. The sensitivity of wild-type and mutated His-mAPAs to inhibition by [(3S)-3-amino-4-sulfanyl-butane-1-sulfonic acid] (EC33), glutamate-thiol (GluSH), glutamate-phosphonic-acid (GluPO₃H₂), methionine-thiol (MetSH) and lysine-thiol (LysSH) was determined by establishing dose-dependent inhibition curves for a final GluNA concentration of 0.5 mM in the presence or absence of 4 mM CaCl₂. These compounds being linear competitive inhibitors, their $K_i$ were calculated from the formula $K_i = IC_{50}/(1 + [GluNA]/K_m)$.

## Statistical analysis

Statistical comparisons were carried out with Student's unpaired t-test. Differences were considered significant if *p* was less than 0.05.

## Results and discussion

APA belongs to the monozinc aminopeptidase family, which requires Zn²⁺ for catalytic activity. In addition, the substrate specificity of APA, unlike that of other monozinc aminopeptidases, is enhanced by Ca²⁺. Indeed, APA preferentially hydrolyzes peptides with N-terminal
acidic amino acids and this is due in part to the presence of Ca²⁺ in the S1 subsite of the APA active site [22, 35]. In the present work, we used the recent crystal structure of human APA to build a 3D model of APA complexed with EC33, a specific and selective APA inhibitor targeting the S1 subsite to explore if other residues, in addition to Ca²⁺, participate to APA acidic substrate specificity.

## Molecular modeling of human APA and docking of the APA inhibitor EC33

Using the recent crystal structure of hAPA (PDB ID 4KXD, [28]), with a bound Glu (Fig 1A and 1C), we built a 3D model of hAPA complexed with EC33 (Fig 1B and 1D). EC33 interacted in the S1 subsite of the APA active site with the Zn²⁺ atom via its thiol group and Glu-215, Glu-352 and Glu-408 (Glu-223, Glu-360 and Glu-416 in hAPA) via their free alpha amino

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

E

887

<table><tr><td>Human APA</td><td>YLVNRYTLNNR</td><td>NLGRIVTI.AEPFNOTE</td><td>895</td></tr><tr><td>Mouse APA</td><td>YLVSRTINDR</td><td>YLGRIVTI.AEPFNOTE</td><td>886</td></tr><tr><td>Human APN</td><td>KLFNDYGGGSF</td><td>SFSNLIQAVTRRFSTE</td><td>911</td></tr><tr><td>Mouse APN</td><td>KLFENYGGGSF</td><td>SFANLIQGVTRRFSSSE</td><td>910</td></tr><tr><td>Human APB</td><td></td><td></td><td></td></tr><tr><td>Mouse APB</td><td></td><td></td><td></td></tr><tr><td>Human LTA4h</td><td></td><td></td><td></td></tr><tr><td>Mouse LTA4h</td><td></td><td></td><td></td></tr><tr><td>Human TRHDE</td><td>I L N T R Y G E A L F</td><td>F M N S K L I S G V T E F L N T E</td><td>963</td></tr><tr><td>Mouse TRHDE</td><td>I L N T R Y G E A L F</td><td>F M N S K L I S G V T E F L N T E</td><td>964</td></tr><tr><td>Human IRAP2</td><td>KL V Q K F P L GSV</td><td>TI Q N I V A G S T Y L F S T K</td><td>972</td></tr><tr><td>Mouse IRAP2</td><td>KL V H K F H L GSV</td><td>TI Q S I V A G S T H L F S T K</td><td>972</td></tr></table>

**Fig 1.** Human Aminopeptidase A with bound glutamate in the crystallographic structure and the proposed binding mode of the inhibitor EC33 after molecular docking in the 3D-model of human APA.

(A and C), S1 subsite visualized in the crystallographic structure of APA co-crystallized with the glutamate (green) (PDB ID 4KXD). (B and D), S1 subsite visualized in the 3D-model of human APA after molecular docking of EC33 (orange). Protein residues either blocking the visualization or interacting with ligands are represented with a transparent surface. Predicted hydrogen or metallic bonds are drawn as dashed lines. Interacting protein residues are labeled (C and D) accordingly to the human APA numeration. (E), Alignment of the human and mouse APA amino acid sequences with the sequences of other monozinc aminopeptidases from the M1 family. Residues in bold are the Arg-887 of human APA and the Arg-878 of mouse APA. APA, Aminopeptidase A; APN, Aminopeptidase N; APB, Aminopeptidase B; LTA4H, Leukotriene A4 hydrolase; TRHDE, Thyrotropin releasing degrading enzyme; IRAP, Insulin regulated aminopeptidase.

https://doi.org/10.1371/journal.pone.0184237.g001
acidic amino acids and this is due in part to the presence of Ca²⁺ in the S1 subsite of the APA active site [22, 35]. In the present work, we used the recent crystal structure of human APA to build a 3D model of APA complexed with EC33, a specific and selective APA inhibitor targeting the S1 subsite to explore if other residues, in addition to Ca²⁺, participate to APA acidic substrate specificity.

## Molecular modeling of human APA and docking of the APA inhibitor EC33

Using the recent crystal structure of hAPA (PDB ID 4KXD, [28]), with a bound Glu (Fig 1A and 1C), we built a 3D model of hAPA complexed with EC33 (Fig 1B and 1D). EC33 interacted in the S1 subsite of the APA active site with the Zn²⁺ atom via its thiol group and Glu-215, Glu-352 and Glu-408 (Glu-223, Glu-360 and Glu-416 in hAPA) via their free alpha amino

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

E

887

<table><tr><td>Human APA</td><td>YLVNRYTLNNR</td><td>NLGRIVTI.AEPFNOTE</td><td>895</td></tr><tr><td>Mouse APA</td><td>YLVSRTINDR</td><td>YLGRIVTI.AEPFNOTE</td><td>886</td></tr><tr><td>Human APN</td><td>KLFNDYGGGSF</td><td>SFSNLIQAVTRRFSTE</td><td>911</td></tr><tr><td>Mouse APN</td><td>KLFENYGGGSF</td><td>SFANLIQGVTRRFSSSE</td><td>910</td></tr><tr><td>Human APB</td><td></td><td></td><td></td></tr><tr><td>Mouse APB</td><td></td><td></td><td></td></tr><tr><td>Human LTA4h</td><td></td><td></td><td></td></tr><tr><td>Mouse LTA4h</td><td></td><td></td><td></td></tr><tr><td>Human TRHDE</td><td>I L N T R Y G E A L F</td><td>F M N S K L I S G V T E F L N T E</td><td>963</td></tr><tr><td>Mouse TRHDE</td><td>I L N T R Y G E A L F</td><td>F M N S K L I S G V T E F L N T E</td><td>964</td></tr><tr><td>Human IRAP2</td><td>KL V Q K F P L GSV</td><td>T I Q N I V A G S T Y L F S T K</td><td>972</td></tr><tr><td>Mouse IRAP2</td><td>KL V H K F H L GSV</td><td>T I Q S I V A G S T H L F S T K</td><td>972</td></tr></table>

**Fig 1.** Human Aminopeptidase A with bound glutamate in the crystallographic structure and the proposed binding mode of the inhibitor EC33 after molecular docking in the 3D-model of human APA. (A and C), S1 subsite visualized in the crystallographic structure of APA co-crystallized with the glutamate (green) (PDB ID 4KXD). (B and D), S1 subsite visualized in the 3D-model of human APA after molecular docking of EC33 (orange). Protein residues either blocking the visualization or interacting with ligands are represented with a transparent surface. Predicted hydrogen or metallic bonds are drawn as dashed lines. Interacting protein residues are labeled (C and D) accordingly to the human APA numeration. (E), Alignment of the human and mouse APA amino acid sequences with the sequences of other monozinc aminopeptidases from the M1 family. Residues in bold are the Arg-887 of human APA and the Arg-878 of mouse APA. APA, Aminopeptidase A; APN, Aminopeptidase N; APB, Aminopeptidase B; LTA4H, Leukotriene A4 hydrolase; TRHDE, Thyrotropin releasing degrading enzyme; IRAP, Insulin regulated aminopeptidase.

https://doi.org/10.1371/journal.pone.0184237.g001
group. [24]. In addition, the sulfonate group of EC33 was embedded at the bottom of the S1 subsite formed by Asp 221 and Asp 226, which bind the Ca²⁺ atom. All these interactions were previously visualized in the 3D model of the mouse APA ectodomain from residues 79 to 559 that we generated by homology using the X-ray crystal structure of leukotriene A4 hydrolase (LTA4H; EC 3.3.2.6) as a template and the functional data collected from site-directed mutagenesis studies on APA [17].

The crystal structure of APA with a bound Glu showed the presence of Arg-878 (Arg-887 in hAPA) at the bottom of the S1 subsite, interacting with the electronegative carboxylate of Glu. We also visualized the presence of this residue in the S1 subsite of the 3D model of human APA/EC33 complex, the guanidinium moiety of which established an interaction with the sulfonate groups of EC33. Moreover, multiple sequence alignments of monozinc aminopeptidases (Fig 1E) demonstrated the presence of this Arg residue exclusively in APA.

Yang *et al* substituted this residue by an Ala and showed a decrease in GluNA hydrolysis and an increase in ArgNA hydrolysis by the mutated APA versus wild-type APA, all these experiments being performed in the absence of Ca²⁺ [28], suggesting that Arg-878 could participate to APA substrate specificity. However, when conducting these experiments in the absence of Ca²⁺ might influence the results since one of the key features of APA is its stimulation by Ca²⁺ [1]. The enzymatic activity of APA with substrates bearing acidic amino acids at their N-terminal part increases 5 to 7-fold when the Ca²⁺ concentration rises from 0 to 4 mM [2, 22]. Ca²⁺ not only enhances the hydrolysis by APA of N-terminal acidic residues from substrates, but also reduces the hydrolysis of neutral or basic residues [5]. Furthermore, we have shown, by molecular modeling and site-directed mutagenesis studies, that Asp-218 and Asp-213 bind the Ca²⁺ atom into the S1 subsite [26] and that Ca²⁺ interacts with the side chain carboxylate of the P1 residue, ensuring the substrate specificity of APA for N-terminal acidic amino acid residues. The high degree of sequence identity (77%) between hAPA and mAPA (S1 Fig) suggested that both proteins have closely related structures. We therefore investigated the role of Arg-878 of mAPA together with Ca²⁺ by molecular modeling and site-directed mutagenesis.

## Site-directed mutagenesis, expression and purification of recombinant wild-type and mutated His-mAPAs

For this purpose, we replaced Arg-878 with an alanine (R878A) or a lysine (R878K) residue by site-directed mutagenesis. Wild-type and mutated His-mAPAs were then stably expressed in CHO cells. We first verified that the mutations did not affect the production and the processing of the recombinant proteins by investigating the subcellular distribution of wild-type and mutated His-mAPAs in stably transfected CHO cells by immunofluorescence staining. Confo-cal microscopy analysis of wild-type and mutated His-mAPAs showed that all recombinant His-mAPAs were located at the plasma membrane consistent with the correct processing of these enzymes (Fig 2A). Wild-type and mutated His-mAPAs stably expressed in CHO cells were then purified by metal affinity chromatography and purified proteins were subjected to SDS-PAGE/silver staining analysis (Fig 2B) and Western blot analysis (Fig 2C). We showed that all recombinant His-mAPAs displayed a major mature form of 180-kDa (Fig 2B and 2C) that corresponded to the mature monomeric glycosylated APA sorted from the Golgi appara-tus. SDS-PAGE and silver staining analysis showed that the enzymes reached a purity higher than 90% with an overall yield of 72% and a purification factor of 118- fold.

## Enzymatic activity of purified wild-type and mutated His-mAPAs

We then determined the Ca²⁺ activation profile of purified recombinant wild-type and mutated His-mAPAs, using an acidic substrate, α-L-glutamyl-β-naphthylamide (GluNA) as
group. [24]. In addition, the sulfonate group of EC33 was embedded at the bottom of the S1 subsite formed by Asp 221 and Asp 226, which bind the Ca²⁺ atom. All these interactions were previously visualized in the 3D model of the mouse APA ectodomain from residues 79 to 559 that we generated by homology using the X-ray crystal structure of leukotriene A4 hydrolase (LTA4H; EC 3.3.2.6) as a template and the functional data collected from site-directed mutagenesis studies on APA [17].

The crystal structure of APA with a bound Glu showed the presence of Arg-878 (Arg-887 in hAPA) at the bottom of the S1 subsite, interacting with the electronegative carboxylate of Glu. We also visualized the presence of this residue in the S1 subsite of the 3D model of human APA/EC33 complex, the guanidinium moiety of which established an interaction with the sulfonate groups of EC33. Moreover, multiple sequence alignments of monozinc aminopeptidases (Fig 1E) demonstrated the presence of this Arg residue exclusively in APA.

Yang *et al* substituted this residue by an Ala and showed a decrease in GluNA hydrolysis and an increase in ArgNA hydrolysis by the mutated APA versus wild-type APA, all these experiments being performed in the absence of Ca²⁺ [28], suggesting that Arg-878 could participate to APA substrate specificity. However, when conducting these experiments in the absence of Ca²⁺ might influence the results since one of the key features of APA is its stimulation by Ca²⁺ [1]. The enzymatic activity of APA with substrates bearing acidic amino acids at their N-terminal part increases 5 to 7-fold when the Ca²⁺ concentration rises from 0 to 4 mM [2, 22]. Ca²⁺ not only enhances the hydrolysis by APA of N-terminal acidic residues from substrates, but also reduces the hydrolysis of neutral or basic residues [5]. Furthermore, we have shown, by molecular modeling and site-directed mutagenesis studies, that Asp-218 and Asp-213 bind the Ca²⁺ atom into the S1 subsite [26] and that Ca²⁺ interacts with the side chain carboxylate of the P1 residue, ensuring the substrate specificity of APA for N-terminal acidic amino acid residues. The high degree of sequence identity (77%) between hAPA and mAPA (S1 Fig) suggested that both proteins have closely related structures. We therefore investigated the role of Arg-878 of mAPA together with Ca²⁺ by molecular modeling and site-directed mutagenesis.

## Site-directed mutagenesis, expression and purification of recombinant wild-type and mutated His-mAPAs

For this purpose, we replaced Arg-878 with an alanine (R878A) or a lysine (R878K) residue by site-directed mutagenesis. Wild-type and mutated His-mAPAs were then stably expressed in CHO cells. We first verified that the mutations did not affect the production and the processing of the recombinant proteins by investigating the subcellular distribution of wild-type and mutated His-mAPAs in stably transfected CHO cells by immunofluorescence staining. Confo-cal microscopy analysis of wild-type and mutated His-mAPAs showed that all recombinant His-mAPAs were located at the plasma membrane consistent with the correct processing of these enzymes (Fig 2A). Wild-type and mutated His-mAPAs stably expressed in CHO cells were then purified by metal affinity chromatography and purified proteins were subjected to SDS-PAGE/silver staining analysis (Fig 2B) and Western blot analysis (Fig 2C). We showed that all recombinant His-mAPAs displayed a major mature form of 180-kDa (Fig 2B and 2C) that corresponded to the mature monomeric glycosylated APA sorted from the Golgi appara-tus. SDS-PAGE and silver staining analysis showed that the enzymes reached a purity higher than 90% with an overall yield of 72% and a purification factor of 118- fold.

## Enzymatic activity of purified wild-type and mutated His-mAPAs

We then determined the Ca²⁺ activation profile of purified recombinant wild-type and mutated His-mAPAs, using an acidic substrate, α-L-glutamyl-β-naphthylamide (GluNA) as
<figure><img src="image_5.png" /><figcaption>**Fig 2.** Expression and purification of wild-type and mutated mouse recombinant His-APAs. (A), Transfected CHO cells stably expressing the wild-type or mutated Xpress-His-mAPA construct were fixed and immunolabeled with a mouse monoclonal anti-Xpress primary antibody, which was detected with an Alexa 488-conjugated anti-mouse secondary antibody. Immunofluorescence was monitored by confocal microscopy. The *bar* indicates 10 µm. (B and C), Purified recombinant wild-type and mutated Xpress-His-mAPA were analyzed by SDS-PAGE 4–12% and proteins were silver stained (B) or transferred into nitrocellulose membrane and subjected to Western blot analysis using mouse monoclonal anti-X-press antibody (C). The arrows indicate the monomeric form of APA.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0184237.g002

previously described [22]. We observed that wild-type mAPA was strongly activated by Ca²⁺ (EC₅₀ = 22 ± 5 µM), with enzymatic activity increasing by 7-fold as the Ca²⁺ concentration raised from 0 to 4 mM, as previously reported [22]. The replacement of Arg-878 by a lysine residue, conserving the positive charge, did not change the ability of this enzyme to be activated by Ca²⁺ when compared to that of the wild-type. By contrast, the replacement of Arg-878 by an alanine, a neutral residue devoid of charge, very slightly reduced activation by Ca²⁺ (Fig 3A). This finding is consistent with the localization of Arg-878 (Arg-887 in hAPA) into the APA active site 3D model, but showed that Arg-878 does not contribute to Ca²⁺ binding in the S1 subsite. However, the maximal enzymatic activities of R878A and R878K mutated mAPAs measured in the presence of a supramaximal concentration of Ca²⁺ were significantly decreased (by a factor of 2.4 and 1.7, respectively, (Fig 3B) when compared to wild-type mAPA. Since these mutated enzymes exhibited a similar pattern of membrane expression when compared to wild-type mAPA, these defect suggest that Arg-878 could contribute to substrate recognition and hydrolysis.

## Kinetic parameters for purified recombinant wild-type and mutated His-mAPAs

This led us to determine the kinetic parameters (*K*<sub>m</sub> and *k*<sub>cat</sub>) for wild-type and mutated His-APAs in the absence or presence of 4 mM Ca²⁺ using GluNA as a substrate. The results are
<figure><img src="image_5.png" /><figcaption>**Fig 2.** Expression and purification of wild-type and mutated mouse recombinant His-APAs. (A), Transfected CHO cells stably expressing the wild-type or mutated Xpress-His-mAPA construct were fixed and immunolabeled with a mouse monoclonal anti-Xpress primary antibody, which was detected with an Alexa 488-conjugated anti-mouse secondary antibody. Immunofluorescence was monitored by confocal microscopy. The *bar* indicates 10 µm. (B and C), Purified recombinant wild-type and mutated Xpress-His-mAPA were analyzed by SDS-PAGE 4–12% and proteins were silver stained (B) or transferred into nitrocellulose membrane and subjected to Western blot analysis using mouse monoclonal anti-X-press antibody (C). The arrows indicate the monomeric form of APA.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0184237.g002

previously described [22]. We observed that wild-type mAPA was strongly activated by Ca²⁺ (EC₅₀ = 22 ± 5 µM), with enzymatic activity increasing by 7-fold as the Ca²⁺ concentration raised from 0 to 4 mM, as previously reported [22]. The replacement of Arg-878 by a lysine residue, conserving the positive charge, did not change the ability of this enzyme to be activated by Ca²⁺ when compared to that of the wild-type. By contrast, the replacement of Arg-878 by an alanine, a neutral residue devoid of charge, very slightly reduced activation by Ca²⁺ (Fig 3A). This finding is consistent with the localization of Arg-878 (Arg-887 in hAPA) into the APA active site 3D model, but showed that Arg-878 does not contribute to Ca²⁺ binding in the S1 subsite. However, the maximal enzymatic activities of R878A and R878K mutated mAPAs measured in the presence of a supramaximal concentration of Ca²⁺ were significantly decreased (by a factor of 2.4 and 1.7, respectively, (Fig 3B) when compared to wild-type mAPA. Since these mutated enzymes exhibited a similar pattern of membrane expression when compared to wild-type mAPA, these defect suggest that Arg-878 could contribute to substrate recognition and hydrolysis.

## Kinetic parameters for purified recombinant wild-type and mutated His-mAPAs

This led us to determine the kinetic parameters (*K*<sub>m</sub> and *k*<sub>cat</sub>) for wild-type and mutated His-APAs in the absence or presence of 4 mM Ca²⁺ using GluNA as a substrate. The results are
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /><figcaption>**Fig 3.** Effects of the divalent cation Ca²⁺ on recombinant wild-type and mutated mouse APAs enzymatic activities. Purified mAPAs (2 µg) were incubated at 37°C for 30 min with 0.5 mM GluNA with various concentrations of CaCl₂ in a final volume of 100 µl of 50 mM Tris-HCl buffer (pH7.4). (A), APA enzymatic activities expressed as a percentage of maximal hydrolysis velocity obtained with 4 mM Ca²⁺ for each purification vs log of Ca²⁺ concentration. (B), Effect of increasing Ca²⁺ concentration on APA enzymatic activities.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0184237.g003

summarized in Table 1. In the absence of Ca²⁺, the cleavage efficiency (kcat/Km) of the mutated enzymes for the acidic substrate, GluNA was clearly reduced when compared to that of the wild-type, but less markedly (by a factor 25) for R878A than for R878K (by a factor 281). This was due to a net decrease in the affinity of the mutated enzymes for GluNA together with a slight decrease in hydrolysis velocity. The addition of Ca²⁺ enhanced wild-type APA activity by increasing the affinity of the enzyme for GluNA as well as substrate hydrolysis velocity, increasing cleavage efficiency (kcat/Km) by 25-fold. In the same conditions, a lower cleavage

<table><caption>**Table 1.** Kinetic parameters for wild type and mutated mAPAs in the absence or presence of Ca²⁺ using HGluβNA as a substrate.</caption><thead><tr><th rowspan="2">Recombinant APAs</th><th colspan="3">No Ca<sup>2+</sup></th><th colspan="3">4 mM Ca<sup>2+</sup></th></tr><tr><th>K<sub>m</sub><br>(µM)</th><th>k<sub>cat</sub><br>(s-1)</th><th>k<sub>cat</sub>/K<sub>m</sub><br>(s-1/mM)</th><th>K<sub>m</sub><br>(µM)</th><th>k<sub>cat</sub><br>(s-1)</th><th>k<sub>cat</sub>/K<sub>m</sub><br>(s-1/mM)</th></tr></thead><tbody><tr><td>Wild type</td><td>165 ± 13.5</td><td>13.9 ± 0.1</td><td>84.3</td><td>39.1 ± 3.9</td><td>81.7 ± 4.0</td><td>2090</td></tr><tr><td>R878A</td><td>1500 ± 79.6***</td><td>4.95 ± 0.1***</td><td>3.3</td><td>221 ± 4.0***</td><td>22.4 ± 0.1***</td><td>102</td></tr><tr><td>R878K</td><td>12500 ± 1450***</td><td>3.12 ± 0.4***</td><td>0.3</td><td>269 ± 5.1***</td><td>135 ± 1.3***</td><td>502</td></tr></tbody></table>

K_m and k_cat values are the mean ± S.E.M from six to eight separate experiments performed in duplicate

*** p<0.001, significant when compared to the corresponding wild type value

https://doi.org/10.1371/journal.pone.0184237.t001
A

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

B

<figure><img src="image_7.png" /><figcaption>**Fig 3.** Effects of the divalent cation Ca²⁺ on recombinant wild-type and mutated mouse APAs enzymatic activities. Purified mAPAs (2 µg) were incubated at 37°C for 30 min with 0.5 mM GluNA with various concentrations of CaCl₂ in a final volume of 100 µl of 50 mM Tris-HCl buffer (pH7.4). (A), APA enzymatic activities expressed as a percentage of maximal hydrolysis velocity obtained with 4 mM Ca²⁺ for each purification vs log of Ca²⁺ concentration. (B), Effect of increasing Ca²⁺ concentration on APA enzymatic activities.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0184237.g003

summarized in Table 1. In the absence of Ca²⁺, the cleavage efficiency ($k_{cat}$/K_m) of the mutated enzymes for the acidic substrate, GluNA was clearly reduced when compared to that of the wild-type, but less markedly (by a factor 25) for R878A than for R878K (by a factor 281). This was due to a net decrease in the affinity of the mutated enzymes for GluNA together with a slight decrease in hydrolysis velocity. The addition of Ca²⁺ enhanced wild-type APA activity by increasing the affinity of the enzyme for GluNA as well as substrate hydrolysis velocity, increasing cleavage efficiency ($k_{cat}$/K_m) by 25-fold. In the same conditions, a lower cleavage

<table><caption>**Table 1.** Kinetic parameters for wild type and mutated mAPAs in the absence or presence of Ca²⁺ using HGluβNA as a substrate.</caption><thead><tr><th rowspan="2">Recombinant APAs</th><th colspan="3">No Ca<sup>2+</sup></th><th colspan="3">4 mM Ca<sup>2+</sup></th></tr><tr><th>K<sub>m</sub><br>(µM)</th><th>k<sub>cat</sub><br>(s-1)</th><th>k<sub>cat</sub>/K<sub>m</sub><br>(s-1/mM)</th><th>K<sub>m</sub><br>(µM)</th><th>k<sub>cat</sub><br>(s-1)</th><th>k<sub>cat</sub>/K<sub>m</sub><br>(s-1/mM)</th></tr></thead><tbody><tr><td>Wild type</td><td>165 ± 13.5</td><td>13.9 ± 0.1</td><td>84.3</td><td>39.1 ± 3.9</td><td>81.7 ± 4.0</td><td>2090</td></tr><tr><td>R878A</td><td>1500 ± 79.6***</td><td>4.95 ± 0.1***</td><td>3.3</td><td>221 ± 4.0***</td><td>22.4 ± 0.1***</td><td>102</td></tr><tr><td>R878K</td><td>12500 ± 1450***</td><td>3.12 ± 0.4***</td><td>0.3</td><td>269 ± 5.1***</td><td>135 ± 1.3***</td><td>502</td></tr></tbody></table>

$K_m$ and $k_{cat}$ values are the mean ± S.E.M from six to eight separate experiments performed in duplicate

*** p<0.001, significant when compared to the corresponding wild type value

https://doi.org/10.1371/journal.pone.0184237.t001
<table><caption>Table 2. Ki values (nM) for several inhibitors with wild-type and mutated recombinant mAPAs targeting the S1 subsite in absence of calcium.</caption><thead><tr><th></th><th>WT</th><th>R878A</th><th>R878K</th></tr></thead><tbody><tr><th>GluSH</th><td>1570 ± 270</td><td>1720 ± 190<sup>n.s</sup></td><td>> 10000</td></tr><tr><th>MetSH</th><td>750 ± 46</td><td>475 ± 140<sup>***</sup></td><td>1250 ± 57<sup>***</sup></td></tr><tr><th>LysSH</th><td>317 ± 30</td><td>118 ± 13<sup>***</sup></td><td>378 ± 20<sup>n.s</sup></td></tr></tbody></table>

Ki values are the mean ± S.E.M from seven separate experiments with duplicate determinations

*** p<0.001 and n.s (not significant) when compared with the corresponding wild type values

https://doi.org/10.1371/journal.pone.0184237.t002

efficiency of the mutated enzymes for GluNA was observed when compared to wild-type but less markedly than in the absence of Ca²⁺, by a factor 21 for R878A and a factor 4 for R878K. This was due to a decrease in the affinity of both mutated enzymes for GluNA (by a factor 6) but also to a decrease in hydrolysis velocity for R878A (by a factor 4) whereas the *k*<sub>cat</sub> for R878K was unchanged. These findings suggested that Arg-878 of mAPA participates in the establishment of the Michaelis complex as well as in catalysis either in the presence or in the absence of Ca²⁺.

## Inhibitory potency of various classes of compounds against GluNA on purified recombinant wild-type and mutated His-mAPAs

To further investigate the functional role of Arg-878, we evaluated the binding mode of various inhibitors mimicking the Michaelis complex such as the β-amino thiol inhibitors, EC33, GluSH, LysSH and MetSH by measuring their inhibitory potencies (*K*<sub>i</sub>) against GluNA hydrolysis by wild-type and mutated His-mAPAs in the absence or in the presence of 4 mM Ca²⁺. These compounds interact with the amine binding site, within the S1 subsite via their acidic, basic, or neutral side chains, and with the Zn²⁺ atom via their chelating thiol group. We also used the inhibitor, GluPO₃H₂ which behaves as an analog of the transition state, mimicking the interaction of the substrate with the enzyme at the transition state by its phosphonate group that mimics the tetrahedral intermediate formed during catalysis. The structures of the inhibitors and the results are summarized in Table 2 and Table 3.

In the absence of Ca²⁺, we showed that R878A was slightly better inhibited (maximally inhibited by a factor of 2) by basic (LysSH) or neutral (MetSH) inhibitors than wild-type APA whereas R878K was less inhibited by acidic (GluSH) or neutral inhibitors than wild-type. The presence of 4 mM Ca²⁺ enhanced GluSH inhibitory potency on wild-type mAPA, R878A and R878K mutants by 5.0, 1.7 and 6.6-fold, respectively, when compared to the corresponding *K*<sub>i</sub> values obtained in the absence of Ca²⁺. In contrast, the presence of Ca²⁺ decreased the inhibitory potency of MetSH and LysSH on wild-type and mutated APAs by a factor of 12 and 42 for wild-type, 7 and 64 for R878A and 4 and 37 for R878K when compared to the *K*<sub>i</sub> values obtained in the absence of Ca²⁺. This showed that the presence of Ca²⁺ partially restored

<table><caption>Table 3. Ki values (nM) of several inhibitors with wild-type and mutated mAPAs targeting the S1 subsite in presence of calcium.</caption><thead><tr><th></th><th>WT</th><th>R878A</th><th>R878K</th></tr></thead><tbody><tr><th>GluSH</th><td>314 ± 33</td><td>983 ± 160<sup>***</sup></td><td>2100 ± 440<sup>***</sup></td></tr><tr><th>EC33</th><td>268 ± 22</td><td>2000 ± 200<sup>***</sup></td><td>1520 ± 126<sup>***</sup></td></tr><tr><th>GluPO<sub>3</sub>H<sub>2</sub></th><td>74 ± 7.7</td><td>6500 ± 240<sup>***</sup></td><td>619 ± 21<sup>***</sup></td></tr><tr><th>MetSH</th><td>9180 ± 770</td><td>3220 ± 200<sup>***</sup></td><td>5140 ± 560<sup>***</sup></td></tr><tr><th>LysSH</th><td>13400 ± 1700</td><td>7510 ± 660<sup>***</sup></td><td>13900 ± 390<sup>n.s</sup></td></tr></tbody></table>

Ki values are the mean ± S.E.M from seven separate experiments with duplicate determinations

*** p<0.001 and n.s (not significant) when compared with the corresponding wild-type values

https://doi.org/10.1371/journal.pone.0184237.t003
<table><caption>Table 2. Ki values (nM) for several inhibitors with wild-type and mutated recombinant mAPAs targeting the S1 subsite in absence of calcium.</caption><thead><tr><th></th><th>WT</th><th>R878A</th><th>R878K</th></tr></thead><tbody><tr><th>GluSH</th><td>1570 ± 270</td><td>1720 ± 190<sup>n.s</sup></td><td>> 10000</td></tr><tr><th>MetSH</th><td>750 ± 46</td><td>475 ± 140<sup>***</sup></td><td>1250 ± 57<sup>***</sup></td></tr><tr><th>LysSH</th><td>317 ± 30</td><td>118 ± 13<sup>***</sup></td><td>378 ± 20<sup>n.s</sup></td></tr></tbody></table>

Ki values are the mean ± S.E.M from seven separate experiments with duplicate determinations

*** p<0.001 and n.s (not significant) when compared with the corresponding wild type values

https://doi.org/10.1371/journal.pone.0184237.t002

efficiency of the mutated enzymes for GluNA was observed when compared to wild-type but less markedly than in the absence of Ca²⁺, by a factor 21 for R878A and a factor 4 for R878K. This was due to a decrease in the affinity of both mutated enzymes for GluNA (by a factor 6) but also to a decrease in hydrolysis velocity for R878A (by a factor 4) whereas the *k*<sub>cat</sub> for R878K was unchanged. These findings suggested that Arg-878 of mAPA participates in the establishment of the Michaelis complex as well as in catalysis either in the presence or in the absence of Ca²⁺.

## Inhibitory potency of various classes of compounds against GluNA on purified recombinant wild-type and mutated His-mAPAs

To further investigate the functional role of Arg-878, we evaluated the binding mode of various inhibitors mimicking the Michaelis complex such as the β-amino thiol inhibitors, EC33, GluSH, LysSH and MetSH by measuring their inhibitory potencies (*K*<sub>i</sub>) against GluNA hydrolysis by wild-type and mutated His-mAPAs in the absence or in the presence of 4 mM Ca²⁺. These compounds interact with the amine binding site, within the S1 subsite via their acidic, basic, or neutral side chains, and with the Zn²⁺ atom via their chelating thiol group. We also used the inhibitor, GluPO₃H₂ which behaves as an analog of the transition state, mimicking the interaction of the substrate with the enzyme at the transition state by its phosphonate group that mimics the tetrahedral intermediate formed during catalysis. The structures of the inhibitors and the results are summarized in Table 2 and Table 3.

In the absence of Ca²⁺, we showed that R878A was slightly better inhibited (maximally inhibited by a factor of 2) by basic (LysSH) or neutral (MetSH) inhibitors than wild-type APA whereas R878K was less inhibited by acidic (GluSH) or neutral inhibitors than wild-type. The presence of 4 mM Ca²⁺ enhanced GluSH inhibitory potency on wild-type mAPA, R878A and R878K mutants by 5.0, 1.7 and 6.6-fold, respectively, when compared to the corresponding *K*<sub>i</sub> values obtained in the absence of Ca²⁺. In contrast, the presence of Ca²⁺ decreased the inhibitory potency of MetSH and LysSH on wild-type and mutated APAs by a factor of 12 and 42 for wild-type, 7 and 64 for R878A and 4 and 37 for R878K when compared to the *K*<sub>i</sub> values obtained in the absence of Ca²⁺. This showed that the presence of Ca²⁺ partially restored

<table><caption>Table 3. Ki values (nM) of several inhibitors with wild-type and mutated mAPAs targeting the S1 subsite in presence of calcium.</caption><thead><tr><th></th><th>WT</th><th>R878A</th><th>R878K</th></tr></thead><tbody><tr><th>GluSH</th><td>314 ± 33</td><td>983 ± 160<sup>***</sup></td><td>2100 ± 440<sup>***</sup></td></tr><tr><th>EC33</th><td>268 ± 22</td><td>2000 ± 200<sup>***</sup></td><td>1520 ± 126<sup>***</sup></td></tr><tr><th>GluPO<sub>3</sub>H<sub>2</sub></th><td>74 ± 7.7</td><td>6500 ± 240<sup>***</sup></td><td>619 ± 21<sup>***</sup></td></tr><tr><th>MetSH</th><td>9180 ± 770</td><td>3220 ± 200<sup>***</sup></td><td>5140 ± 560<sup>***</sup></td></tr><tr><th>LysSH</th><td>13400 ± 1700</td><td>7510 ± 660<sup>***</sup></td><td>13900 ± 390<sup>n.s</sup></td></tr></tbody></table>

Ki values are the mean ± S.E.M from seven separate experiments with duplicate determinations

*** p<0.001 and n.s (not significant) when compared with the corresponding wild-type values

https://doi.org/10.1371/journal.pone.0184237.t003
inhibitor specificity of the mutated mAPAs for inhibitors with a N-terminal acidic amino acid residue (GluSH>MetSH>LysSH), demonstrating that the Ca²⁺ atom is a major component of APA acidic inhibitor specificity. However, the affinity of EC33 and GluSH for the mutated enzymes (R878A and R878K) remains weaker (factor 3 to 9) than that for the wild-type suggesting that the interaction between the acidic side-chain of the inhibitors with APA required the positive charge of Ca²⁺ in the S1 subsite but not only. Indeed, our data showed that the conservative charge of the lysine substitution is not enough to maintain an optimal positioning of the acidic inhibitors into the active site versus wild-type mAPA. This can be explained by a different orientation of the inhibitor due to the longer side chain of the lysine residue when compared to that of the arginine. Since our data showed that Arg-878 plays a role not only in the establishment of the Michaelis complex but also in catalysis, we investigated the effects of the mutations of Arg-878 on GluPO₃H₂ binding. GluPO₃H₂, in the presence of Ca²⁺, inhibited R878A 87-times less efficiently than wild-type, whereas it inhibited R878K 8-times less efficiently than wild-type. This different behavior of GluPO₃H₂ regarding the mutated enzymes is in agreement with data obtained on the hydrolysis of GluNA. Altogether, our data suggested that Arg-878 played a key role in substrate specificity of mAPA towards substrates with a N-terminal acidic amino acid residue by contributing with the Ca²⁺ in the S1 subsite to interactions with the electronegative charge of the P1 substrate side chain.

## Molecular dynamics simulation analysis

To better understand the effects of mutations on inhibitor binding and substrate accessibility, we performed MD simulations with wild-type and mutated hAPAs in the presence of Ca²⁺ and various inhibitors. Pair interaction analysis revealed that the docked inhibitors (EC33 and GluPO₃H₂) were kept stable at the active site of wild-type hAPA during the entire simulation (Fig 4A) and the average interaction energies between the inhibitor and the protein versus those obtained between the inhibitor and the solvent were determined (Fig 4A). Wild-type displayed an average interaction energy of -270.1 ± 17.5 kcal/mol for EC33, and -299.5 ± 17.6 kcal/mol for GluPO₃H₂. Following virtual mutation, those energies became less spontaneous (less negative), but still favorable for inhibitor binding, agreeing with our experimental data. The substitution of Arg-887 by an alanine (R878A mutant in mAPA) decreased the binding of the inhibitors with an average interaction energy of -278.1 ± 8.9 kcal/mol for EC33 and of -254.5 ± 13.3 kcal/mol for GluPO₃H₂. Regarding the more conservative R887K mutation, interaction energies were close to those of the wild-type model (-225.6 ± 10.9 kcal/mol for EC33 and -298.9 ± 11.6 kcal/mol for GluPO₃H₂). Indeed, the substitution of Arg-887 by an alanine leads to a strong decrease in the energy of interaction of GluPO₃H₂ with the hAPA active site that can be due to a major modification of the electropositive environment of the S1 subsite. In addition, the higher interaction between GluPO₃H₂ and the solvent observed with the R887A mutant led to a decrease in the affinity of the inhibitor towards the R887A mutated enzyme. On the other hand, the fact that the inhibitory potency of GluPO₃H₂ towards the R878K mutant is less impacted than that towards R878A mutant, could be explained by the fact that the ε amine moiety of the lysine residue maintains the electropositive environment of the S1 subsite. This hypothesis is further supported by the fact that isolated pair interaction analysis performed between single residues and the inhibitors revealed that the Lys-887 substitute is able to interact spontaneously with both EC33 and GluPO₃H₂ (Fig 4B). For the EC33 inhibitor, the calculated energies were -59.64 ± 5.05 kcal/mol for the Arg-887 of the wild-type hAPA, -6.58 ± 1.35 kcal/mol for the Ala-887 of the R887A mutant and -27.19 ± 6.35 kcal/mol for the Lys-887 of the R887K mutant. For the GluPO₃H₂ inhibitor, the calculated values were -59.26 ± 9.32 kcal/mol for the Arg-887 of the wild-type hAPA, -11.15 ± 1.03 kcal/mol for
inhibitor specificity of the mutated mAPAs for inhibitors with a N-terminal acidic amino acid residue (GluSH>MetSH>LysSH), demonstrating that the Ca²⁺ atom is a major component of APA acidic inhibitor specificity. However, the affinity of EC33 and GluSH for the mutated enzymes (R878A and R878K) remains weaker (factor 3 to 9) than that for the wild-type suggesting that the interaction between the acidic side-chain of the inhibitors with APA required the positive charge of Ca²⁺ in the S1 subsite but not only. Indeed, our data showed that the conservative charge of the lysine substitution is not enough to maintain an optimal positioning of the acidic inhibitors into the active site versus wild-type mAPA. This can be explained by a different orientation of the inhibitor due to the longer side chain of the lysine residue when compared to that of the arginine. Since our data showed that Arg-878 plays a role not only in the establishment of the Michaelis complex but also in catalysis, we investigated the effects of the mutations of Arg-878 on GluPO₃H₂ binding. GluPO₃H₂, in the presence of Ca²⁺, inhibited R878A 87-times less efficiently than wild-type, whereas it inhibited R878K 8-times less efficiently than wild-type. This different behavior of GluPO₃H₂ regarding the mutated enzymes is in agreement with data obtained on the hydrolysis of GluNA. Altogether, our data suggested that Arg-878 played a key role in substrate specificity of mAPA towards substrates with a N-terminal acidic amino acid residue by contributing with the Ca²⁺ in the S1 subsite to interactions with the electronegative charge of the P1 substrate side chain.

## Molecular dynamics simulation analysis

To better understand the effects of mutations on inhibitor binding and substrate accessibility, we performed MD simulations with wild-type and mutated hAPAs in the presence of Ca²⁺ and various inhibitors. Pair interaction analysis revealed that the docked inhibitors (EC33 and GluPO₃H₂) were kept stable at the active site of wild-type hAPA during the entire simulation (Fig 4A) and the average interaction energies between the inhibitor and the protein versus those obtained between the inhibitor and the solvent were determined (Fig 4A). Wild-type displayed an average interaction energy of -270.1 ± 17.5 kcal/mol for EC33, and -299.5 ± 17.6 kcal/mol for GluPO₃H₂. Following virtual mutation, those energies became less spontaneous (less negative), but still favorable for inhibitor binding, agreeing with our experimental data. The substitution of Arg-887 by an alanine (R878A mutant in mAPA) decreased the binding of the inhibitors with an average interaction energy of -278.1 ± 8.9 kcal/mol for EC33 and of -254.5 ± 13.3 kcal/mol for GluPO₃H₂. Regarding the more conservative R887K mutation, interaction energies were close to those of the wild-type model (-225.6 ± 10.9 kcal/mol for EC33 and -298.9 ± 11.6 kcal/mol for GluPO₃H₂). Indeed, the substitution of Arg-887 by an alanine leads to a strong decrease in the energy of interaction of GluPO₃H₂ with the hAPA active site that can be due to a major modification of the electropositive environment of the S1 subsite. In addition, the higher interaction between GluPO₃H₂ and the solvent observed with the R887A mutant led to a decrease in the affinity of the inhibitor towards the R887A mutated enzyme. On the other hand, the fact that the inhibitory potency of GluPO₃H₂ towards the R878K mutant is less impacted than that towards R878A mutant, could be explained by the fact that the ε amine moiety of the lysine residue maintains the electropositive environment of the S1 subsite. This hypothesis is further supported by the fact that isolated pair interaction analysis performed between single residues and the inhibitors revealed that the Lys-887 substitute is able to interact spontaneously with both EC33 and GluPO₃H₂ (Fig 4B). For the EC33 inhibitor, the calculated energies were -59.64 ± 5.05 kcal/mol for the Arg-887 of the wild-type hAPA, -6.58 ± 1.35 kcal/mol for the Ala-887 of the R887A mutant and -27.19 ± 6.35 kcal/mol for the Lys-887 of the R887K mutant. For the GluPO₃H₂ inhibitor, the calculated values were -59.26 ± 9.32 kcal/mol for the Arg-887 of the wild-type hAPA, -11.15 ± 1.03 kcal/mol for
<figure><img src="image_8.png" /><figcaption>**Fig 4.** Analysis of the molecular dynamics simulations trajectories reveals loss of inhibitor affinity for the mutated hAPAs and a larger active site pocket for the R887K mutant. (A), The pair interaction energies analysis between inhibitor and protein (blue) or inhibitor and solvent (red) indicate less spontaneous interactions between inhibitors and the mutant R887A when compared to wild-type hAPA. (B) Pair interaction energies analysis between residue 887 in the wild-type or mutated hAPAs and the inhibitors shows that the Arg-887 from the wild-type APA and the Lys-887 of the R887K mutant spontaneously interact with EC33 and GluPO₃H₂. (C), Pocket volume was extracted from the trajectories after removing the docked inhibitors, showing that the R887K mutant exhibits a larger binding pocket. Each point in the scatter plots correspond to the values calculated from a different frame in the MD trajectories, vertical bars illustrate the minimum and maximum values for each calculation, horizontal bars correspond to the means.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0184237.g004

the Ala-887 of the R887A mutant, and -20,22 ± 3.58 kcal/mol for the Lys-887 of the R887K mutant. Furthermore, the distance between the ε nitrogen atom from the Lys-887 sidechain of the R887K mutant and the sulfonate moiety of the EC33 inhibitor fluctuated from 3.06 to 6.39 Å during the entire simulation, suggesting that this substitution was not engaged in a stable conformation with the EC33 inhibitor (S2 Fig), which might explain why this mutant was still impacted regarding inhibitor binding. Together those results imply that, despite still being
<figure><img src="image_8.png" /><figcaption>**Fig 4.** Analysis of the molecular dynamics simulations trajectories reveals loss of inhibitor affinity for the mutated hAPAs and a larger active site pocket for the R887K mutant. (A), The pair interaction energies analysis between inhibitor and protein (blue) or inhibitor and solvent (red) indicate less spontaneous interactions between inhibitors and the mutant R887A when compared to wild-type hAPA. (B) Pair interaction energies analysis between residue 887 in the wild-type or mutated hAPAs and the inhibitors shows that the Arg-887 from the wild-type APA and the Lys-887 of the R887K mutant spontaneously interact with EC33 and GluPO₃H₂. (C), Pocket volume was extracted from the trajectories after removing the docked inhibitors, showing that the R887K mutant exhibits a larger binding pocket. Each point in the scatter plots correspond to the values calculated from a different frame in the MD trajectories, vertical bars illustrate the minimum and maximum values for each calculation, horizontal bars correspond to the means.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0184237.g004

the Ala-887 of the R887A mutant, and -20,22 ± 3.58 kcal/mol for the Lys-887 of the R887K mutant. Furthermore, the distance between the ε nitrogen atom from the Lys-887 sidechain of the R887K mutant and the sulfonate moiety of the EC33 inhibitor fluctuated from 3.06 to 6.39 Å during the entire simulation, suggesting that this substitution was not engaged in a stable conformation with the EC33 inhibitor (S2 Fig), which might explain why this mutant was still impacted regarding inhibitor binding. Together those results imply that, despite still being
able to interact with the EC33 inhibitor, Lys-878 is not as efficient as Arg-887 from the wild-
type APA. We then analyzed the effect of the mutations on the volume of the S1 subsite. For
this purpose, we extracted the binding site volume information from our MD simulations tra-
jectories (Fig 4C, S1, S2 and S3 Movies). Although the lysine being a bulkier amino acid than
alanine, the binding site volume of this mutant was significantly higher than that of R887A or
wild-type hAPA, regardless of the inhibitor (EC33 or GluPO3H2). The binding pocket volume
averages for the APA + EC33 simulation were, 1194 ± 196 Å³, 821 ± 130 Å³ and 1828 ± 253 Å³
for wild-type, R887A and R887K, respectively. For the hAPA + GluPO3H2 system, those values
were 831 ± 137 Å³, 833 ± 149 Å³ and 1634 ± 206 Å³ for wild-type, R887A and R887K, respec-
tively. The change in the volume of the S1 subsite observed with the R878K mutant, may
impair the binding and/or the optimal positioning of the substrate as well as its hydrolysis,
resulting in a lower cleavage efficiency (kcat/Km) of GluNA and a lower inhibitory potency of
GluPO3H2 with the R878K mutated enzyme when compared to the wild-type enzyme. This
strengthened the conclusion that Arg-878 is located in the mAPA S1 subsite and therefore
may drive APA substrate specificity for N-terminal acidic amino acid residues. These changes
were more marked in the absence than in the presence of Ca²⁺, supporting that Arg-878
together with Ca²⁺ plays a crucial role in APA acidic substrate specificity.

Hydrolysis of various substrates by wild-type and mutated His-mAPAs

To strengthen the hypothesis that Arg-878 is involved in substrate specificity for N-terminal acidic amino acid residues, we evaluated the efficiency of hydrolysis of the wild-type and mutated His-mAPAs with various substrates differing by their side-chains. We therefore determined the kinetic parameters of wild-type and the mutated enzymes for the substrates, GluNA, AspNA, AlaNA and LysNA, in the absence or presence of Ca²⁺ and we compared their substrate hydrolysis profiles (Fig 5, S1 and S2 Tables). In the absence of Ca²⁺, the substrate hydrolysis profile for wild-type mAPA was GluNA = AspNA > LysNA > AlaNA. This profile was changed by Arg-878 substitutions. The substrate profile for R878A was LysNA > AlaNA > GluNA = AspNA and for R878K AlaNA = LysNA >> GluNA > AspNA. The analysis of the kinetics parameters of wild-type and mutated mAPAs showed that the modifications in the substrate hydrolysis profile observed for the mutated mAPAs were due to an increased substrate hydrolysis velocity and to a better affinity of the mutated enzymes for AlaNA and LysNA when compared to that for GluNA or AspNA.

The addition of Ca²⁺ amplifies the preference of wild-type mAPA for acidic substrates. Thus, GluNA and AspNA were much more efficiently hydrolyzed than LysNA, which was hydrolyzed more efficiently than AlaNA (GluNA > AspNA >> LysNA > AlaNA). This profile was partially maintained in R878K with the conservative mutation. In contrast, the R878A mutant in the presence of Ca²⁺, similarly hydrolyzed GluNA, AspNA, AlaNA and LysNA. Altogether, these data confirmed that Arg-878 participates with Ca²⁺ to APA acidic substrate specificity. Additionally, in the presence of Ca²⁺, the hydrolysis efficiency of GluNA by R878K and R878A remained 4- and 21-fold weaker than that of the wild-type, respectively. Since in the cerebrospinal fluid, the Ca²⁺ concentration fluctuates from 4.3 mM to 5.9 mM [36], which is 40 times higher than the EC₅₀ of Ca²⁺ for APA. This suggests the constant presence of Ca²⁺ in the S1 subsite of endogenous APA, which consequently recognizes only APA substrates or inhibitors with an N-terminal acidic side chain, highlighting the key role of Arg-878 with Ca²⁺ in APA substrate specificity.

In conclusion, we showed that the guanidium moiety of Arg-878 in the mAPA S1
subsite interacts with the N-terminal acidic amino acid residue of the substrate or inhibitor,
participating with Ca²⁺ to the substrate specificity of APA for N-terminal acidic amino acid resi-
dues. We also showed that Arg-878 contributes to ensure an optimal positioning of the substrate
able to interact with the EC33 inhibitor, Lys-878 is not as efficient as Arg-887 from the wild-
type APA. We then analyzed the effect of the mutations on the volume of the S1 subsite. For
this purpose, we extracted the binding site volume information from our MD simulations tra-
jectories (Fig 4C, S1, S2 and S3 Movies). Although the lysine being a bulkier amino acid than
alanine, the binding site volume of this mutant was significantly higher than that of R887A or
wild-type hAPA, regardless of the inhibitor (EC33 or GluPO3H2). The binding pocket volume
averages for the APA + EC33 simulation were, 1194 ± 196 Å³, 821 ± 130 Å³ and 1828 ± 253 Å³
for wild-type, R887A and R887K, respectively. For the hAPA + GluPO3H2 system, those values
were 831 ± 137 Å³, 833 ± 149 Å³ and 1634 ± 206 Å³ for wild-type, R887A and R887K, respec-
tively. The change in the volume of the S1 subsite observed with the R878K mutant, may
impair the binding and/or the optimal positioning of the substrate as well as its hydrolysis,
resulting in a lower cleavage efficiency (kcat/Km) of GluNA and a lower inhibitory potency of
GluPO3H2 with the R878K mutated enzyme when compared to the wild-type enzyme. This
strengthened the conclusion that Arg-878 is located in the mAPA S1 subsite and therefore
may drive APA substrate specificity for N-terminal acidic amino acid residues. These changes
were more marked in the absence than in the presence of Ca²⁺, supporting that Arg-878
together with Ca²⁺ plays a crucial role in APA acidic substrate specificity.

Hydrolysis of various substrates by wild-type and mutated His-mAPAs

To strengthen the hypothesis that Arg-878 is involved in substrate specificity for N-terminal acidic amino acid residues, we evaluated the efficiency of hydrolysis of the wild-type and mutated His-mAPAs with various substrates differing by their side-chains. We therefore determined the kinetic parameters of wild-type and the mutated enzymes for the substrates, GluNA, AspNA, AlaNA and LysNA, in the absence or presence of Ca²⁺ and we compared their substrate hydrolysis profiles (Fig 5, S1 and S2 Tables). In the absence of Ca²⁺, the substrate hydrolysis profile for wild-type mAPA was GluNA = AspNA > LysNA > AlaNA. This profile was changed by Arg-878 substitutions. The substrate profile for R878A was LysNA > AlaNA > GluNA = AspNA and for R878K AlaNA = LysNA >> GluNA > AspNA. The analysis of the kinetics parameters of wild-type and mutated mAPAs showed that the modifications in the substrate hydrolysis profile observed for the mutated mAPAs were due to an increased substrate hydrolysis velocity and to a better affinity of the mutated enzymes for AlaNA and LysNA when compared to that for GluNA or AspNA.

The addition of Ca²⁺ amplifies the preference of wild-type mAPA for acidic substrates. Thus, GluNA and AspNA were much more efficiently hydrolyzed than LysNA, which was hydrolyzed more efficiently than AlaNA (GluNA > AspNA >> LysNA > AlaNA). This profile was partially maintained in R878K with the conservative mutation. In contrast, the R878A mutant in the presence of Ca²⁺, similarly hydrolyzed GluNA, AspNA, AlaNA and LysNA. Altogether, these data confirmed that Arg-878 participates with Ca²⁺ to APA acidic substrate specificity. Additionally, in the presence of Ca²⁺, the hydrolysis efficiency of GluNA by R878K and R878A remained 4- and 21-fold weaker than that of the wild-type, respectively. Since in the cerebrospinal fluid, the Ca²⁺ concentration fluctuates from 4.3 mM to 5.9 mM [36], which is 40 times higher than the EC₅₀ of Ca²⁺ for APA. This suggests the constant presence of Ca²⁺ in the S1 subsite of endogenous APA, which consequently recognizes only APA substrates or inhibitors with an N-terminal acidic side chain, highlighting the key role of Arg-878 with Ca²⁺ in APA substrate specificity.

In conclusion, we showed that the guanidium moiety of Arg-878 in the mAPA S1
subsite interacts with the N-terminal acidic amino acid residue of the substrate or inhibitor,
participating with Ca²⁺ to the substrate specificity of APA for N-terminal acidic amino acid resi-
dues. We also showed that Arg-878 contributes to ensure an optimal positioning of the substrate
<figure><img src="image_9.png" /><figcaption>**Fig 5. Hydrolysis of various substrates, GluNA, AspNA, AlaNA and LysNA, by recombinant mAPAs.**
Recombinant wild-type and mutated mAPAs hydrolysis efficiencies were evaluated using GluNA, AspNA, AlaNA and LysNA as a substrate in the absence (A) or presence (B) of 4 mM Ca²⁺. The data are the mean ± S.E of five to seven experiments performed in duplicate. *** p<0.001; ** p<0.01; * p<0.05 and n.s (not significant) when compared with the value obtained for GluNA.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0184237.g005

during catalysis, thereby optimizing the hydrolysis of the substrate scissile peptide bond. In addition, our work showed that, whereas Arg-878 is far from the consensus sequence HEXXH and other active site contributing residues in term of primary sequence, this residue is very close from the active site with regard to the 3D structure. Indeed, the X-ray structure of APA shows that the C-terminal domain that contains the Arg-878 residue is in close contact to the region of the active site and complement by this way the structure of the S1 subsite. These new insights into the organization of the APA S1 subsite showed that this subsite is at the interface of the N- and C-terminal domains of APA and that the C-terminal domain has not only a chaperone activity as previously demonstrated but also a functional role in the S1 subsite, in substrate recognition [22, 37].

**Supporting information**

**S1 Table. Kinetic parameters for wild type and mutated mAPAs in absence of calcium, using different synthetic substrates.** *K*<sub>m</sub> and *k*<sub>cat</sub> values are the mean ± S.E.M from six to eight separate experiments performed in duplicate.
(DOCX)
<figure><img src="image_9.png" /><figcaption>**Fig 5. Hydrolysis of various substrates, GluNA, AspNA, AlaNA and LysNA, by recombinant mAPAs.**
Recombinant wild-type and mutated mAPAs hydrolysis efficiencies were evaluated using GluNA, AspNA, AlaNA and LysNA as a substrate in the absence (A) or presence (B) of 4 mM Ca²⁺. The data are the mean ± S.E of five to seven experiments performed in duplicate. *** p<0.001; ** p<0.01; * p<0.05 and n.s (not significant) when compared with the value obtained for GluNA.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0184237.g005

during catalysis, thereby optimizing the hydrolysis of the substrate scissile peptide bond. In addition, our work showed that, whereas Arg-878 is far from the consensus sequence HEXXH and other active site contributing residues in term of primary sequence, this residue is very close from the active site with regard to the 3D structure. Indeed, the X-ray structure of APA shows that the C-terminal domain that contains the Arg-878 residue is in close contact to the region of the active site and complement by this way the structure of the S1 subsite. These new insights into the organization of the APA S1 subsite showed that this subsite is at the interface of the N- and C-terminal domains of APA and that the C-terminal domain has not only a chaperone activity as previously demonstrated but also a functional role in the S1 subsite, in substrate recognition [22, 37].

**Supporting information**

**S1 Table. Kinetic parameters for wild type and mutated mAPAs in absence of calcium, using different synthetic substrates.** *K*<sub>m</sub> and *k*<sub>cat</sub> values are the mean ± S.E.M from six to eight separate experiments performed in duplicate.
(DOCX)
**S2 Table.** Kinetic parameters for wild type and mutated mAPAs in presence of calcium, using different synthetic substrates. *K*<sub>m</sub> and *k*<sub>cat</sub> values are the mean ± S.E.M from six to eight separate experiments performed in duplicate.
(DOCX)

**S3 Table.** List of docked molecules used to identify key residues of the APA S1 subsite.
(DOCX)

**S4 Table.** X-Ray based score penalties applied during molecular docking.
(DOCX)

**S1 Fig.** Sequence alignment of hAPA and mAPA. Conserved residues are shaded in black, semi-conserved residues are shaded in grey. The secondary structure information was extracted from the 4KXD pdb structure, where red, blue and yellow indicate α-helix, 3₁₀-helix, and extended structure (β-sheets), respectively. Active site residues are highlighted by an orange box. Residues responsible for Zn²⁺ binding are marked by orange diamonds, the catalytic glutamate is marked by a green triangle, and the glutamate responsible for N-terminus recognition is marked by a blue triangle. The mouse Arg-878 and the homologous residue in the human sequence Arg-887 are shaded in red and marked by an arrow.
(TIF)

**S2 Fig.** Evolution of the distance between the epsilon amide from Lys-887 of the R887K mutant and the sulfonate moiety of the EC33 inhibitor. Snapshot ‘A’ corresponds to the one where the lowest distance was observed. Snapshot ‘B’ corresponds to a conformation where the distances were the same as the average of all frames. Snapshot ‘C’ corresponds to the largest distance observed. In A, B and C, the inhibitor carbon atoms are colored in orange, Zn²⁺ is colored in ochre, Ca²⁺ ion is colored in purple. The Lys-887 residue is depicted as translucent surface representation over a licorice representation in order to display bond distances with the sulfonate moiety of the EC33 inhibitor.
(TIF)

**S1 Movie.** Analysis of the binding site volume through the MD simulation trajectory (50 ns) of the wild-type humanAPA.
(MPG)

**S2 Movie.** Analysis of the binding site volume through the MD simulation trajectory (50 ns) of the R887A humanAPA.
(MPG)

**S3 Movie.** Analysis of the binding site volume through the MD simulation trajectory (50 ns) of the R887K humanAPA.
(MPG)

**Acknowledgments**

We thank Jérémie Teillon of the CIRB for microscopy core facility. We acknowledge GENCI, as part of the molecular dynamics simulations were performed thanks to an allocation of computer time on the CINES supercomputer OCCIGEN.

**Author Contributions**

**Conceptualization:** Hugo de Almeida.
**S2 Table.** Kinetic parameters for wild type and mutated mAPAs in presence of calcium, using different synthetic substrates. *K*<sub>m</sub> and *k*<sub>cat</sub> values are the mean ± S.E.M from six to eight separate experiments performed in duplicate.

(DOCX)

**S3 Table.** List of docked molecules used to identify key residues of the APA S1 subsite.

(DOCX)

**S4 Table.** X-Ray based score penalties applied during molecular docking.

(DOCX)

**S1 Fig.** Sequence alignment of hAPA and mAPA. Conserved residues are shaded in black, semi-conserved residues are shaded in grey. The secondary structure information was extracted from the 4KXD pdb structure, where red, blue and yellow indicate α-helix, 3₁₀-helix, and extended structure (β-sheets), respectively. Active site residues are highlighted by an orange box. Residues responsible for Zn²⁺ binding are marked by orange diamonds, the catalytic glutamate is marked by a green triangle, and the glutamate responsible for N-terminus recognition is marked by a blue triangle. The mouse Arg-878 and the homologous residue in the human sequence Arg-887 are shaded in red and marked by an arrow.

(TIF)

**S2 Fig.** Evolution of the distance between the epsilon amide from Lys-887 of the R887K mutant and the sulfonate moiety of the EC33 inhibitor. Snapshot ‘A’ corresponds to the one where the lowest distance was observed. Snapshot ‘B’ corresponds to a conformation where the distances were the same as the average of all frames. Snapshot ‘C’ corresponds to the largest distance observed. In A, B and C, the inhibitor carbon atoms are colored in orange, Zn²⁺ is colored in ochre, Ca²⁺ ion is colored in purple. The Lys-887 residue is depicted as translucent surface representation over a licorice representation in order to display bond distances with the sulfonate moiety of the EC33 inhibitor.

(TIF)

**S1 Movie.** Analysis of the binding site volume through the MD simulation trajectory (50 ns) of the wild-type humanAPA.

(MPG)

**S2 Movie.** Analysis of the binding site volume through the MD simulation trajectory (50 ns) of the R887A humanAPA.

(MPG)

**S3 Movie.** Analysis of the binding site volume through the MD simulation trajectory (50 ns) of the R887K humanAPA.

(MPG)

**Acknowledgments**

We thank Jérémie Teillon of the CIRB for microscopy core facility. We acknowledge GENCI, as part of the molecular dynamics simulations were performed thanks to an allocation of computer time on the CINES supercomputer OCCIGEN.

**Author Contributions**

**Conceptualization:** Hugo de Almeida.
**Data curation:** Pierre Couvineau, Hugo de Almeida, Xavier Iturrioz.

**Formal analysis:** Pierre Couvineau, Hugo de Almeida, Xavier Iturrioz.

**Funding acquisition:** Catherine Llorens-Cortes.

**Investigation:** Pierre Couvineau, Hugo de Almeida.

**Methodology:** Hugo de Almeida.

**Supervision:** Bernard Maigret, Catherine Llorens-Cortes, Xavier Iturrioz.

**Writing – original draft:** Pierre Couvineau, Hugo de Almeida, Catherine Llorens-Cortes, Xavier Iturrioz.

**Writing – review & editing:** Pierre Couvineau, Hugo de Almeida, Catherine Llorens-Cortes, Xavier Iturrioz.

## References

1. Danielsen EM, Noren O, Sjostrom H, Ingram J, Kenny AJ. Proteins of the kidney microvillar membrane. Aspartate aminopeptidase: purification by immunoadsorbent chromatography and properties of the detergent- and proteinase-solubilized forms. The Biochemical journal. 1980; 189(3):591–603. PMID: 7011318; PubMed Central PMCID: PMCPMC1162040.

2. Glenner GG, Mc MP, Folk JE. A mammalian peptidase specific for the hydrolysis of N-terminal alpha-L-glutamyl and aspartyl residues. Nature. 1962; 194:867. PMID: 13899213.

3. Healy DP, Wilk S. Localization of immunoreactive glutamyl aminopeptidase in rat brain. II. Distribution and correlation with angiotensin II. Brain research. 1993; 606(2):295–303. PMID: 8490722.

4. Nagatsu I, Nagatsu T, Yamamoto T, Glenner GG, Mehl JW. Purification of aminopeptidase A in human serum and degradation of angiotensin II by the purified enzyme. Biochimica et biophysica acta. 1970; 198(2):255–70. PMID: 4313532.

5. Feracci H, Benajiba A, Gorvel JP, Doumeng C, Maroux S. Enzymatic and immunological properties of the protease form of aminopeptidase N and A from pig and rabbit intestinal brush border. Biochimica et biophysica acta. 1981; 658(1):148–57. PMID: 6163468.

6. Lojda Z, Gossrau R. Study on aminopeptidase A. Histochemistry. 1980; 67(3):267–90. PMID: 6772609.

7. Zini S, Masdehors P, Lenke Z, Fournie-Zaluski MC, Roques BP, Corvol P, et al. Aminopeptidase A: distribution in rat brain nuclei and increased activity in spontaneously hypertensive rats. Neuroscience. 1997; 78(4):1187–93. PMID: 9174084.

8. Chauvel EN, Llorens-Cortes C, Coric P, Wilk S, Roques BP, Fournie-Zaluski MC. Differential inhibition of aminopeptidase A and aminopeptidase N by new beta-amino thiols. Journal of medicinal chemistry. 1994; 37(18):2950–7. PMID: 7915326.

9. Zini S, Fournie-Zaluski MC, Chauvel E, Roques BP, Corvol P, Llorens-Cortes C. Identification of metabolic pathways of brain angiotensin II and III using specific aminopeptidase inhibitors: predominant role of angiotensin III in the control of vasopressin release. Proceedings of the National Academy of Sciences of the United States of America. 1996; 93(21):11968–73. PMID: 8876246; PubMed Central PMCID: PMC38167.

10. Reaux A, Fournie-Zaluski MC, David C, Zini S, Roques BP, Corvol P, et al. Aminopeptidase A inhibitors as potential central antihypertensive agents. Proceedings of the National Academy of Sciences of the United States of America. 1999; 96(23):13415–20. PMID: 10557335; PubMed Central PMCID: PMC23962.

11. Fournie-Zaluski MC, Fassot C, Valentin B, Djordijijevic D, Reaux-Le Goazigo A, Corvol P, et al. Brain renin-angiotensin system blockade by systemically active aminopeptidase A inhibitors: a potential treatment of salt-dependent hypertension. Proceedings of the National Academy of Sciences of the United States of America. 2004; 101(20):7775–80. https://doi.org/10.1073/pnas.0402312101 PMID: 15136730; PubMed Central PMCID: PMC419682.

12. Bodineau L, Frugiere A, Marc Y, Inguimbert N, Fassot C, Balavoine F, et al. Orally active aminopeptidase A inhibitors reduce blood pressure: a new strategy for treating hypertension. Hypertension. 2008; 51(5):1318–25. https://doi.org/10.1161/HYPERTENSIONAHA.107.098772 PMID: 18362226.

13. Marc Y, Gao J, Balavoine F, Michaud A, Roques BP, Llorens-Cortes C. Central antihypertensive effects of orally active aminopeptidase A inhibitors in spontaneously hypertensive rats. Hypertension. 2012; 60 (2):411–8. https://doi.org/10.1161/HYPERTENSIONAHA.112.190942 PMID: 22710644.
**Data curation:** Pierre Couvineau, Hugo de Almeida, Xavier Iturrioz.

**Formal analysis:** Pierre Couvineau, Hugo de Almeida, Xavier Iturrioz.

**Funding acquisition:** Catherine Llorens-Cortes.

**Investigation:** Pierre Couvineau, Hugo de Almeida.

**Methodology:** Hugo de Almeida.

**Supervision:** Bernard Maigret, Catherine Llorens-Cortes, Xavier Iturrioz.

**Writing – original draft:** Pierre Couvineau, Hugo de Almeida, Catherine Llorens-Cortes, Xavier Iturrioz.

**Writing – review & editing:** Pierre Couvineau, Hugo de Almeida, Catherine Llorens-Cortes, Xavier Iturrioz.

## References

1. Danielsen EM, Noren O, Sjostrom H, Ingram J, Kenny AJ. Proteins of the kidney microvillar membrane. Aspartate aminopeptidase: purification by immunoadsorbent chromatography and properties of the detergent- and proteinase-solubilized forms. The Biochemical journal. 1980; 189(3):591–603. PMID: 7011318; PubMed Central PMCID: PMCPMC1162040.

2. Glenner GG, Mc MP, Folk JE. A mammalian peptidase specific for the hydrolysis of N-terminal alpha-L-glutamyl and aspartyl residues. Nature. 1962; 194:867. PMID: 13899213.

3. Healy DP, Wilk S. Localization of immunoreactive glutamyl aminopeptidase in rat brain. II. Distribution and correlation with angiotensin II. Brain research. 1993; 606(2):295–303. PMID: 8490722.

4. Nagatsu I, Nagatsu T, Yamamoto T, Glenner GG, Mehl JW. Purification of aminopeptidase A in human serum and degradation of angiotensin II by the purified enzyme. Biochimica et biophysica acta. 1970; 198(2):255–70. PMID: 4313532.

5. Feracci H, Benajiba A, Gorvel JP, Doumeng C, Maroux S. Enzymatic and immunological properties of the protease form of aminopeptidase N and A from pig and rabbit intestinal brush border. Biochimica et biophysica acta. 1981; 658(1):148–57. PMID: 6163468.

6. Lojda Z, Gossrau R. Study on aminopeptidase A. Histochemistry. 1980; 67(3):267–90. PMID: 6772609.

7. Zini S, Masdehors P, Lenke Z, Fournie-Zaluski MC, Roques BP, Corvol P, et al. Aminopeptidase A: distribution in rat brain nuclei and increased activity in spontaneously hypertensive rats. Neuroscience. 1997; 78(4):1187–93. PMID: 9174084.

8. Chauvel EN, Llorens-Cortes C, Coric P, Wilk S, Roques BP, Fournie-Zaluski MC. Differential inhibition of aminopeptidase A and aminopeptidase N by new beta-amino thiols. Journal of medicinal chemistry. 1994; 37(18):2950–7. PMID: 7915326.

9. Zini S, Fournie-Zaluski MC, Chauvel E, Roques BP, Corvol P, Llorens-Cortes C. Identification of metabolic pathways of brain angiotensin II and III using specific aminopeptidase inhibitors: predominant role of angiotensin III in the control of vasopressin release. Proceedings of the National Academy of Sciences of the United States of America. 1996; 93(21):11968–73. PMID: 8876246; PubMed Central PMCID: PMC38167.

10. Reaux A, Fournie-Zaluski MC, David C, Zini S, Roques BP, Corvol P, et al. Aminopeptidase A inhibitors as potential central antihypertensive agents. Proceedings of the National Academy of Sciences of the United States of America. 1999; 96(23):13415–20. PMID: 10557335; PubMed Central PMCID: PMC23962.

11. Fournie-Zaluski MC, Fassot C, Valentin B, Djordijijevic D, Reaux-Le Goazigo A, Corvol P, et al. Brain renin-angiotensin system blockade by systemically active aminopeptidase A inhibitors: a potential treatment of salt-dependent hypertension. Proceedings of the National Academy of Sciences of the United States of America. 2004; 101(20):7775–80. https://doi.org/10.1073/pnas.0402312101 PMID: 15136730; PubMed Central PMCID: PMC419682.

12. Bodineau L, Frugiere A, Marc Y, Inguimbert N, Fassot C, Balavoine F, et al. Orally active aminopeptidase A inhibitors reduce blood pressure: a new strategy for treating hypertension. Hypertension. 2008; 51(5):1318–25. https://doi.org/10.1161/HYPERTENSIONAHA.107.098772 PMID: 18362226.

13. Marc Y, Gao J, Balavoine F, Michaud A, Roques BP, Llorens-Cortes C. Central antihypertensive effects of orally active aminopeptidase A inhibitors in spontaneously hypertensive rats. Hypertension. 2012; 60 (2):411–8. https://doi.org/10.1161/HYPERTENSIONAHA.112.190942 PMID: 22710644.
14. Xavier Iturrioz AR, Le Goazigo, Catherine Llorens-Cortes. Aminopeptidase Inhibitors as Anti-Hypertensive Drugs. Proteases in Biology and Disease 2004; 2(Aminopeptidases in Biology and Disease):229-50.

15. Ferreira AJ, Raizada MK. Aminopeptidase A: could it be a novel target for neurogenic hypertension? Hypertension. 2008; 51(5):1273-4. https://doi.org/10.1161/HYPERTENSIONAHA.107.109561 PMID: 18362227.

16. Gao J, Marc Y, Iturrioz X, Leroux V, Balavoine F, Llorens-Cortes C. A new strategy for treating hypertension by blocking the activity of the brain renin-angiotensin system with aminopeptidase A inhibitors. Clin Sci (Lond). 2014; 127(3):135-48. https://doi.org/10.1042/CS20130396 PMID: 24697296.

17. Thunnissen MM, Nordlund P, Haeggstrom JZ. Crystal structure of human leukotriene A(4) hydrolase, a bifunctional enzyme in inflammation. Nature structural biology. 2001; 8(2):131-5. https://doi.org/10.1038/84117 PMID: 11175901.

18. Wang J, Cooper MD. Histidine residue in the zinc-binding motif of aminopeptidase A is critical for enzymatic activity. Proceedings of the National Academy of Sciences of the United States of America. 1993; 90(4):1222-6. PMID: 8433982 PubMed Central PMCID: PMC45844.

19. Vazeux G, Wang J, Corvol P, Llorens-Cortes C. Identification of glutamate residues essential for catalytic activity and zinc coordination in aminopeptidase A. The Journal of biological chemistry. 1996; 271(15):9069-74. PMID: 8621556.

20. Vazeux G, Iturrioz X, Corvol P, Llorens-Cortes C. A tyrosine residue essential for catalytic activity in aminopeptidase A. The Biochemical journal. 1997; 327 (Pt 3):883-9. PMID: 9581570 PubMed Central PMCID: PMC1218871.

21. Vazeux G, Iturrioz X, Corvol P, Llorens-Cortes C. A glutamate residue contributes to the exopeptidase specificity in aminopeptidase A. The Biochemical journal. 1998; 334 (Pt 2):407-13. PMID: 9716499 PubMed Central PMCID: PMC1219703.

22. Iturrioz X, Vazeux G, Celerier J, Corvol P, Llorens-Cortes C. Histidine 450 plays a critical role in catalysis and, with Ca2+, contributes to the substrate specificity of aminopeptidase A. Biochemistry. 2000; 39(11):3061-8. PMID: 10715127.

23. Iturrioz X, Rozenfeld R, Michaud A, Corvol P, Llorens-Cortes C. Study of asparagine 353 in aminopeptidase A: characterization of a novel motif (GXMEN) implicated in exopeptidase specificity of monozinc aminopeptidases. Biochemistry. 2001; 40(48):14440-8. PMID: 11724556.

24. Rozenfeld R, Iturrioz X, Maigret B, Llorens-Cortes C. Contribution of molecular modeling and site-directed mutagenesis to the identification of two structural residues, Arg-220 and Asp-227, in aminopeptidase A. The Journal of biological chemistry. 2002; 277(32):29242-52. https://doi.org/10.1074/jbc.M204406200 PMID: 12042323.

25. Lejczak B, De Choszczak MP, Kafarski P. Inhibition of aminopeptidases by phosphonic acid and phosphinic acid analogues of aspartic and glutamic acids. Journal of enzyme inhibition. 1993; 7(2):97-103. PMID: 7509873.

26. Claperon C, Rozenfeld R, Iturrioz X, Inguimbert N, Okada M, Roques B, et al. Asp218 participates with Asp213 to bind a Ca2+ atom into the S1 subsite of aminopeptidase A: a key element for substrate specificity. The Biochemical journal. 2008; 416(1):37-46. https://doi.org/10.1042/BJ20080471 PMID: 18598240.

27. Goto Y, Hattori A, Mizutani S, Tsujimoto M. Asparatic acid 221 is critical in the calcium-induced modulation of the enzymatic activity of human aminopeptidase A. The Journal of biological chemistry. 2007; 282(51):37074-81. https://doi.org/10.1074/jbc.M707251200 PMID: 17932029.

28. Yang Y, Liu C, Lin YL, Li F. Structural insights into central hypertension regulation by human aminopeptidase A. The Journal of biological chemistry. 2013; 288(35):25638-45. https://doi.org/10.1074/jbc.M113.494955 PMID: 23888046 PubMed Central PMCID: PMC3757224.

29. Lalu K, Lampelo S, Vanha-Perttula T. Characterization of three aminopeptidases purified from maternal serum. Biochimica et biophysica acta. 1986; 873(2):190-7. PMID: 2875736.

30. Jones HC, Keep RF. Brain fluid calcium concentration and response to acute hypercalcaemia during development in the rat. J Physiol. 1988; 402:579-93. PMID: 3236250 PubMed Central PMCID: PMCPMC1191910.

31. Webb B, Sali A. Protein structure modeling with MODELLER. Methods Mol Biol. 2014; 1137:1-15. https://doi.org/10.1007/978-1-4939-0366-5_1 PMID: 24573470.

32. Jones G, Willett P, Glen RC, Leach AR, Taylor R. Development and validation of a genetic algorithm for flexible docking. Journal of molecular biology. 1997; 267(3):727-48. https://doi.org/10.1006/jmbi.1996.0897 PMID: 9126849.

33. Best RB, Zhu X, Shim J, Lopes PE, Mittal J, Feig M, et al. Optimization of the additive CHARMM all-atom protein force field targeting improved sampling of the backbone phi, psi and side-chain chi(1) and
14. Xavier Iturrioz AR, Le Goazigo, Catherine Llorens-Cortes. Aminopeptidase Inhibitors as Anti-Hypertensive Drugs. Proteases in Biology and Disease 2004; 2(Aminopeptidases in Biology and Disease):229-50.

15. Ferreira AJ, Raizada MK. Aminopeptidase A: could it be a novel target for neurogenic hypertension? Hypertension. 2008; 51(5):1273-4. https://doi.org/10.1161/HYPERTENSIONAHA.107.109561 PMID: 18362227.

16. Gao J, Marc Y, Iturrioz X, Leroux V, Balavoine F, Llorens-Cortes C. A new strategy for treating hypertension by blocking the activity of the brain renin-angiotensin system with aminopeptidase A inhibitors. Clin Sci (Lond). 2014; 127(3):135-48. https://doi.org/10.1042/CS20130396 PMID: 24697296.

17. Thunnissen MM, Nordlund P, Haeggstrom JZ. Crystal structure of human leukotriene A(4) hydrolase, a bifunctional enzyme in inflammation. Nature structural biology. 2001; 8(2):131-5. https://doi.org/10.1038/84117 PMID: 11175901.

18. Wang J, Cooper MD. Histidine residue in the zinc-binding motif of aminopeptidase A is critical for enzymatic activity. Proceedings of the National Academy of Sciences of the United States of America. 1993; 90(4):1222-6. PMID: 8433982 PubMed Central PMCID: PMC45844.

19. Vazeux G, Wang J, Corvol P, Llorens-Cortes C. Identification of glutamate residues essential for catalytic activity and zinc coordination in aminopeptidase A. The Journal of biological chemistry. 1996; 271(15):9069-74. PMID: 8621556.

20. Vazeux G, Iturrioz X, Corvol P, Llorens-Cortes C. A tyrosine residue essential for catalytic activity in aminopeptidase A. The Biochemical journal. 1997; 327 (Pt 3):883-9. PMID: 9581570 PubMed Central PMCID: PMC1218871.

21. Vazeux G, Iturrioz X, Corvol P, Llorens-Cortes C. A glutamate residue contributes to the exopeptidase specificity in aminopeptidase A. The Biochemical journal. 1998; 334 (Pt 2):407-13. PMID: 9716499 PubMed Central PMCID: PMC1219703.

22. Iturrioz X, Vazeux G, Celerier J, Corvol P, Llorens-Cortes C. Histidine 450 plays a critical role in catalysis and, with Ca2+, contributes to the substrate specificity of aminopeptidase A. Biochemistry. 2000; 39(11):3061-8. PMID: 10715127.

23. Iturrioz X, Rozenfeld R, Michaud A, Corvol P, Llorens-Cortes C. Study of asparagine 353 in aminopeptidase A: characterization of a novel motif (GXMEN) implicated in exopeptidase specificity of monozinc aminopeptidases. Biochemistry. 2001; 40(48):14440-8. PMID: 11724556.

24. Rozenfeld R, Iturrioz X, Maigret B, Llorens-Cortes C. Contribution of molecular modeling and site-directed mutagenesis to the identification of two structural residues, Arg-220 and Asp-227, in aminopeptidase A. The Journal of biological chemistry. 2002; 277(32):29242-52. https://doi.org/10.1074/jbc.M204406200 PMID: 12042323.

25. Lejczak B, De Choszczak MP, Kafarski P. Inhibition of aminopeptidases by phosphonic acid and phosphinic acid analogues of aspartic and glutamic acids. Journal of enzyme inhibition. 1993; 7(2):97-103. PMID: 7509873.

26. Claperon C, Rozenfeld R, Iturrioz X, Inguimbert N, Okada M, Roques B, et al. Asp218 participates with Asp213 to bind a Ca2+ atom into the S1 subsite of aminopeptidase A: a key element for substrate specificity. The Biochemical journal. 2008; 416(1):37-46. https://doi.org/10.1042/BJ20080471 PMID: 18598240.

27. Goto Y, Hattori A, Mizutani S, Tsujimoto M. Asparatic acid 221 is critical in the calcium-induced modulation of the enzymatic activity of human aminopeptidase A. The Journal of biological chemistry. 2007; 282(51):37074-81. https://doi.org/10.1074/jbc.M707251200 PMID: 17932029.

28. Yang Y, Liu C, Lin YL, Li F. Structural insights into central hypertension regulation by human aminopeptidase A. The Journal of biological chemistry. 2013; 288(35):25638-45. https://doi.org/10.1074/jbc.M113.494955 PMID: 23888046 PubMed Central PMCID: PMC3757224.

29. Lalu K, Lampelo S, Vanha-Perttula T. Characterization of three aminopeptidases purified from maternal serum. Biochimica et biophysica acta. 1986; 873(2):190-7. PMID: 2875736.

30. Jones HC, Keep RF. Brain fluid calcium concentration and response to acute hypercalcaemia during development in the rat. J Physiol. 1988; 402:579-93. PMID: 3236250 PubMed Central PMCID: PMCPMC1191910.

31. Webb B, Sali A. Protein structure modeling with MODELLER. Methods Mol Biol. 2014; 1137:1-15. https://doi.org/10.1007/978-1-4939-0366-5_1 PMID: 24573470.

32. Jones G, Willett P, Glen RC, Leach AR, Taylor R. Development and validation of a genetic algorithm for flexible docking. Journal of molecular biology. 1997; 267(3):727-48. https://doi.org/10.1006/jmbi.1996.0897 PMID: 9126849.

33. Best RB, Zhu X, Shim J, Lopes PE, Mittal J, Feig M, et al. Optimization of the additive CHARMM all-atom protein force field targeting improved sampling of the backbone phi, psi and side-chain chi(1) and
chi(2) dihedral angles. J Chem Theory Comput. 2012; 8(9):3257-73. https://doi.org/10.1021/ct300400x
PMID: 23341755; PubMed Central PMCID: PMCPMC3549273.

34. MacKerell AD, Bashford D, Bellott M, Dunbrack RL, Evanseck JD, Field MJ, et al. All-atom empirical potential for molecular modeling and dynamics studies of proteins. The journal of physical chemistry B. 1998; 102(18):3586-616. https://doi.org/10.1021/jp973084f PMID: 24889800.

35. Claperon C, Banegas-Font I, Iturrioz X, Rozenfeld R, Maigret B, Llorens-Cortes C. Identification of three- onine 348 as a residue involved in aminopeptidase A substrate specificity. The Journal of biological chemistry. 2009; 284(16):10618-26. https://doi.org/10.1074/jbc.M806783200 PMID: 19228697; PubMed Central PMCID: PMC2667749.

36. Robertson WG, Marshall RW. Ionized calcium in body fluids. Crit Rev Clin Lab Sci. 1981; 15(2):85-125. https://doi.org/10.3109/10408368109105869 PMID: 7026165.

37. Rozenfeld R, Muller L, El Messari S, Llorens-Cortes C. The C-terminal domain of aminopeptidase A is an intramolecular chaperone required for the correct folding, cell surface expression, and activity of this monozinc aminopeptidase. The Journal of biological chemistry. 2004; 279(41):43285-95. https://doi.org/10.1074/jbc.M404369200 PMID: 15263000.
chi(2) dihedral angles. J Chem Theory Comput. 2012; 8(9):3257–73. https://doi.org/10.1021/ct300400x
PMID: 23341755; PubMed Central PMCID: PMCPMC3549273.

34. MacKerell AD, Bashford D, Bellott M, Dunbrack RL, Evanseck JD, Field MJ, et al. All-atom empirical potential for molecular modeling and dynamics studies of proteins. The journal of physical chemistry B. 1998; 102(18):3586–616. https://doi.org/10.1021/jp973084f PMID: 24889800.

35. Claperon C, Banegas-Font I, Iturrioz X, Rozenfeld R, Maigret B, Llorens-Cortes C. Identification of three- onine 348 as a residue involved in aminopeptidase A substrate specificity. The Journal of biological chemistry. 2009; 284(16):10618–26. https://doi.org/10.1074/jbc.M806783200 PMID: 19228697; PubMed Central PMCID: PMC2667749.

36. Robertson WG, Marshall RW. Ionized calcium in body fluids. Crit Rev Clin Lab Sci. 1981; 15(2):85–125. https://doi.org/10.3109/10408368109105869 PMID: 7026165.

37. Rozenfeld R, Muller L, El Messari S, Llorens-Cortes C. The C-terminal domain of aminopeptidase A is an intramolecular chaperone required for the correct folding, cell surface expression, and activity of this monozinc aminopeptidase. The Journal of biological chemistry. 2004; 279(41):43285–95. https://doi.org/10.1074/jbc.M404369200 PMID: 15263000.