Aminopeptidases (APs) constitute a group of exopeptidases with closely related activities, able to remove amino-acids from unblocked N-termini of peptide and protein substrates. Among them, metallo-aminopeptidases (metallo-APs) are widely distributed in organisms with representative members in animal cells, plant cells, fungi, parasites and bacteria. Metallo-APs have been broadly sorted into two important subgroups on the basis of the number of zinc ions in their active site that are involved in catalysis as well as on their protein signatures and folding. While some of these enzymes contain only one zinc, others possess a binuclear metal centre, usually referred as a co-catalytic unit. Members of the M1 family belong to the subclan MA(E) of metallopeptidases, also called Glu-zincins and possess a single essential catalytic zinc ion and a common fold, related to their common ancestry that is related to thermolysin. Present mainly as monomers or homo-dimers (in mammals, mainly), metallo-APs occupy various cell compartments and although some are secreted, the vast majority operates as membrane bound ectoenzymes or cytosolic catalysts. Their evolutionary history indicates that they are related to a complex and ancestral gene family that subsequently yielded series of more divergent sequences. For this reason, their biochemical characterization and respective biological importance has, for a long time, been hindered by considerable confusion, due in part to their often large and sometimes overlapping substrate specificities as well as their very close biochemical properties. For the time being, we have a more precise picture of their highly specialized roles depending on their localization with many related paralogous members in mammals (12 in H. sapiens, 11 in rodents), whereas only a handful are known to be encoded in other species (3 or less; only one in P. falciparum).
Human M1 paralogous members have been named Leucyl-cystinyl aminopeptidase (HsIRAP), endoplasmic reticulum aminopeptidase 1 (HsERAP1), endoplasmic reticulum aminopeptidase 2 (HsERAP2), aminopeptidase Q (HsAPQ), puromycin-sensitive aminopeptidase (HsPSA), thyrotropin-releasing hormone-degrading ectoenzyme (HsTRHDE), aminopeptidase A (HsAPA), aminopeptidase N (HsAPN), aminopeptidase O (HsAPO), leukotriene A4 hydrolase (HsLTA4H), arginyl aminopeptidase B (HsAPB) and arginyl aminopeptidase B like 1 (HsRNPL1) and their orthologues in E. coli (EcPepN) and P. falciparum (PfAM1). Their so far recognized biological implications are summarized in Table 1. Broadly, human M1 members play important biological roles in the renin-angiotensin system, the immune system as well as during the inflammatory process, while the bacterial and malarial members are involved in nutrient acquisition.
As a consequence from the definition of the M1 family, all members (the MEROPS database has so far identified 14718 amino acid sequences belonging to this M1 family) share critical characteristics. Sequences alignment of a selection of M1 aminopeptidases highlights the conserved active site residues for substrate recognition and peptide bond hydrolysis that is, the zinc binding motif HExxHx18E, the substrate binding motif, including the GxMEN domain and an extra E or Q distal residue, responsible for the recognition and binding of the key free amino terminus of peptide substrates, as well as the catalytic E residue (part of HExxH motif) and a distal Y residue essential for transition state stabilization in the closed conformations of M1 aminopeptidases (Figure S1).
Although overall amino acid sequence similarity may fall down to a mean value of 20%, a highly conserved 3D structure, reminding that of a sea-horse with 4 domains, has been found for all these M1 APs. Domain I folds as a twisted β-barrel of 200 aa, domain II corresponds to the ancestral thermolysin-like fold, domain III is a beta sandwich built with 2 β-sheets (absent in HsLTA4H) and domain IV is a prominent bowl-like structure composed of α-helices (330 aa) which contains the active site channel entrance (about 40 Å long and 12 Å large).
Dynamic conformational ensembles, involving inter-domains movements, are crucial for catalysis and various conformational states are particularly marked for mammalian enzymes. They include active site opening and closure, particularly important for HsERAP1, HsIRAP and mammalian APN, suggesting that the peptide-binding channel could accommodate variously exposed N-termini of peptides and/or proteins. The influences of protein dimerization on active site accessibility and specificity as well as on various moonlighting functions have still to be more extensively studied for all M1 family of APs.
Substrate specificities have also been responsible for the confusion in biological roles attributed to these enzymes; as illustrated in Table 2, these substrate specificities are rarely discriminant with overlapping specificities being the rule and specific specificities the rare exceptions.
This common substrate specificity, together with a very conserved structure and catalytic mechanism, makes the development of highly selective inhibitors of a given metallo-aminopeptidase obviously challenging. Thereby the most notable inhibitors, depicted in Figure 1, consist in tetrahedral intermediate mimics, as for the widely used bestatin and phosphinic/phosphonic acid derivatives, or zinc-chelating group inserted in a peptide-like scaffold, as for hydroxamic acids.
In Table 3, we report inhibition values for several members of the M1 aminopeptidase family and for three selected bi-metallic aminopeptidases, since selectivity towards these other metallo-aminopeptidases remains of major concern and is indeed currently poorly documented.
Inhibitors targeting mammalian APN have been the most extensively developed, since mammalian APN was one of the first metallo-aminopeptidase to be biochemically and enzymatically characterized (although many of its first assigned biological functions were later demonstrated to involve other related M1 paralogous enzymes). Among them, Bestatin 1 has been very often misused as a “specific” inhibitor of this peculiar APN enzyme, since it turns out to be a highly potent inhibitor of bimetallic enzyme families and a quite broad “unspecific” inhibitor of monometallic M1 aminopeptidases. But overall, it appears obvious that the development of M1 aminopeptidase inhibitors remains mainly in its infancy.
Our research group has previously reported the development of potent and selective inhibitors of APN, a human M1 aminopeptidase involved in angiogenesis and tumour metastasis. This work allowed us to identify aminobenzosuberone moiety as a promising scaffold that demonstrates an exceptional selectivity towards mono-metallic aminopeptidases. In the present work, we have closely considered the evaluation of these inhibitors on a number of potential interesting targets within the M1 family of aminopeptidases and established their inhibition profiles. We have reported herein an improvement of the aminobenzosuberone synthetic route and the biological evaluation of this family of non-peptidic compounds against 6 human and 2 bacterial/parasitic M1 aminopeptidases, plus members of the M17 and M28 families. We have also verified the observed selectivity on a panel of proteinases, which includes representative members of aspartic, cysteine, metallo and serine endopeptidases.
The aminobenzosuberone’s synthetic route, based on our previously described procedure, was optimized to improve reproducibility and avoid the appearance of side-products (Scheme 1). Indeed, silyl enol ether 15, obtained in four steps from commercially available product 14, was submitted to a Rubottom oxidation to yield the key intermediate silyloxy-ketone 16. This reaction was classically performed by action of m-CPBA but regular and unpredictable failure with this reagent led us to consider alternative oxidizing agents. The best results were obtained with in situ generated dimethyl dioxirane (DMDO), which gave only the expected product 16 in a nearly quantitative and reproducible yield (up to 95%). The second optimized step was the conversion of the hydroxy-oxime 17 to its corresponding amino-alcohol 18. Our previously described method for this reduction consisted of a hydrogenation over Raney nickel. Although the procedure worked efficiently, it was very dependent on the quality and supply of Raney nickel. Moreover, in the case of brominated compounds it required a careful monitoring to avoid any aromatic debromination. Thus, several conditions were investigated to obtain compound 18 as Mg/ammonium formate, SmI2, BH3·THF, LiAlH4 or NaBH4/NiCl2·6H2O. Promising results were obtained with this latter but a small amount of debrominated product was still observed. Therefore, addition of NiCl2·6H2O was successfully replaced by addition of CoCl2·6H2O and the expected amino compound was solely obtained with no effect on aromatic bromine, even without strict monitoring of the reaction progress and in acceptable yields (60–80%).
2.2. In Vitro Inhibition of A Representative Panel of M1 Aminopeptidases
We have confirmed and extended the outstanding selectivity of aminobenzosuberone derivatives for M1 aminopeptidases over most other classes of proteinases (aspartic, cysteine, metallo and serine endopeptidase) with IC50 value >> 30 µM in all testes cases (Table S1). As expected, the bimetallic aminopeptidases from the M17 (PfAM17) and M28 (AAP) families were not or only very weakly inhibited by the tested compounds (Table 4).
Among M1 aminopeptidases, the most potent inhibitory values remained on mammalian APN, with a subnanomolar Ki value for compound 21i (Ki = 60 pM). Note that HsLTA4H was not inhibited by this low-molecular weight scaffold, except for 21c (Ki = 19 µM). The SAR determined for mammalian APN was very similar for its M1 orthologues from microorganisms, EcPepN and P. falciparum PfAM1, with an improved affinity when the aminobenzosuberone scaffold was substituted in C-1 or C-4 position with a hydrophobic functional group (for instance 21i exhibited Ki values of 50 nM and 30 nM against EcPepN and PfAM1, respectively). The influence of bromine in C-1 position remained similar for the enzymes of both microorganisms with a 10 fold increase in potency when compared to the unsubstituted derivative 21a. The main differences relied on the lack of synergistic effect (100-fold more potent for APN) when both positions were occupied by these functions (Ph and/or Br). Compounds 21c and 21i showed similar potencies on PfAM1 (Ki values of 50 nM and 30 nM, respectively) and 21i was only 3 times more potent than 21c on EcPepN.
As regard to human enzymes, HsPSA is well-inhibited with Ki values of 55 and 21 nM for 21c and 21i, respectively. Within the oxytocinase subfamily, compounds were less active with inhibitory potencies in the micromolar range, except for 21h on HsIRAP (Ki = 34 nM). Generally, the inhibitors were less potent on HsERAP2 and HsERAP1 even if a similar trend emerged in the SAR studies. Actually, substitution at C-4 position seemed to be more efficient than substitution in the C-1 position, 21c was usually more potent than 21f and 1,4-disubstituted derivatives 21h and 21i displayed the most interesting binding affinity with a precise substitution pattern to achieve an excellent potency.
3.1. Structure-Based Analysis of the Structure-Activity Relationships
In the present work, we focused on subtle amino acid sequence differences in the active site that may explain the various inhibition values experimentally determined for the selected M1 aminopeptidases. In addition, conformational dynamics have been considered for some members of this family among the structures captured during various crystallization processes. Data remain scarce due to the small number of solved 3D structures and no general statements could be outlined, nevertheless our SAR studies enabled to precise the structural environment required for potent inhibition by the aminobenzosuberone derivatives of a given enzyme within the M1 family. Our discussion will essentially focus on “closed” conformation of different aminopeptidases.
The determined structures from various species indicate that the active site of a M1 aminopeptidase is buried inside the protein and its zinc binding site may be located in a spacious cavity with wide openings, or located in a small and gated compartment as it is the case for bacterial/parasite orthologues: this was in agreement with the active site’s cavity volume of different M1 aminopeptidases analysed by KVFinder (Table S2 and Figure S2). Local or intradomain/interdomain movements can cause enlargement or contraction of the active site as encountered for EcPepN or mammalian APN.
However, the architecture and amino acid composition of the active site are highly conserved within the whole family, which permits to predict an identical binding mode of the aminobenzosuberone ligand in the various M1 aminopeptidases similar to the previous model reported for mammalian APN or EcPepN-aminobenzosuberone co-crystallized structures (Figure 2).
3.2. S1 Subsite
The S1 subsite is always described as a cylinder formed by hydrophobic residues and capped by two polar amino acids. In Table 5, key residues of this constitutive recognition site are summarized. Interestingly, a highly variable residue (highlighted in “red”), both in steric size and polarity is located at the entrance of this pocket and may contribute to the diversity of S1 specificities (Table S3 and Figure S3).
The aforementioned variable residue may be flexible upon ligand binding or might be sterically hindered, restricting thus the S1 subsite entrance (Figure 3 and for more details, Figure S3). For instance, in EcPepN it corresponds to Met260 (PDB 2DQM or 5MFS) with a temperature factor value of 28–37 Å2, in agreement with a local significant flexibility of its side chain, which acts as a cushion to accommodate various side chains in P1 position of a substrate/ligand. This local movement of Met residue can cause the contraction of the active site cavity. KVFinder analysis revealed that the cavity volume of EcPepN-bestatin 1 complex (PDB 2DQM) was higher (638.5 Å3), compared to 531.25 Å3 for the EcPepN-aminobenzuberone 21c complex (PDB 5MFS) (Table S2). In this latter complex, this Met260 residue was found at two closed positions and moved slightly inward the S1 subsite, interrupting and shortening its typical cylindrical shape. The distance between the sulphur atom of Met260 and the centroid of the aromatic ring of the aminobenzosuberone scaffold is 5.3–5.7Å, which is consistent with sulphur-carbon van der Waals contacts and the sulphur lone pairs might interact with the edge of the aromatic core, providing additional stability (Figure S4). On the other hand, Val459 in PfAM1 (PDB 3EBH) and Pro333 in HsERAP2 (PDB 4JBS) might be considered as straight and bulky, shortening the S1 pocket. In HsLTA4H, Tyr267 (PDB 1HS6) in this position induces the formation of a narrow S1 tunnel. Moreover, domain III is missing in this latter enzyme and the presence of Tyr267, Lys565 and particularly Tyr378, which extends into the active site cavity pointing inward toward the Zn ion, reduced substantially the active site cavity volume (Figure S5). As a result of this spatial arrangement, steric clashes with the aromatic core of aminobenzosuberone might explain the very weak inhibitory activities of this scaffold on HsLTA4H.
The aminobenzosuberone derivatives are markedly more potent on M1 aminopeptidases exhibiting a shallow S1 cavity with small to moderate size (EcPepN, HsAPN, PfAM1 and a borderline case for HsIRAP). At the top of the elongated S1 subsite cylinder are located two or three cap residues (presented in bold type in Table 5 and depicted in stick in Figure 3a and Figure S3) from the C-terminal domain of the proteins. These residues are known to influence S1 subsite specificity. Indeed, depending on the nature of the residues and through the conformational flexibility of their side chains, they can mitigate any potential unfavourable steric clashes and favour van der Waals, H-bond or electrostatic interactions with the substrates/inhibitors. Moreover, the closing of the S1 subsite induces the expulsion of water molecules from the active site into the bulk (entropically favourable) and may strengthen both van der Waals and electrostatic interactions. The analysis of superimposed structures showed that our aminobenzosuberone inhibitors were not expected to penetrate into the S1 cavity as deeply as other inhibitors (like bestatin) (Figure S4a and Figure S6), indicating they should therefore interact more efficiently with shallow S1 subsite.
However, depending on the nature of the variable amino acid residue, the aminobenzosuberone core might slightly move upward to minimize any potential steric hindrance with this special ‘gatekeeper’ residue. Substituents on position-1 most likely induced steric conflicts with residues at the entrance of S1 cavities. More suitable substitutions in position-9 should offer new scopes for novel series of more potent and selective aminobenzosuberone derivatives (Figure S6).
3.3. S1′ Subsite
In contrast to S1, the S1′ subsite is an open cavity of 9 to 10 residues for all enzymes. For instance, in EcPepN, Gly261, Ala262, Tyr275, Arg293, Val294, His297, Glu298, Val324, Asp327, Tyr381 compose this subsite. The careful analysis of the various 3D structures showed that only 2 to 3 residues (highlighted in blue and red in Table 6) seem to be involved in S1′ plasticity, thereby controlling the cavity width and depth (Figure 4a and Figure S7). Two amino acids, which are Tyr275 (located on the sidewall opposite to the catalytic zinc ion) and Arg293 (partially capping the S1′ cavity) for EcPepN, vary in size and polarity within this enzyme family, as summarized in Table 6 and Figure 4a.
For example, PfAM1 Arg489 is found more buried in the cavity than HsAPN Arg381, reducing thus slightly the width of S1′ subsite, whereas the EcPepN Tyr275 and HsAPN Arg381 are positioned at a similar location, modifying thereby the local polarity. Also, the EcPepN Arg293 side chain is located further inside the cavity reducing the cavity depth, whereas the less bulky and polar PfAM1 Thr492 and HsAPN Thr384 occupy a similar position. Finally, local flexibility was observed within captured structures. Indeed, the side chain of residue HsAPN Arg381, HsERAP2 Arg366 and EcPepN Arg293 can adopt different positions upon ligand binding. A similar observation might also be noticed for EcPepN Asp327, this carboxylate side chain is either orientated in- or outside S1′ cavity, with the inside orientation potentially involved in further interactions as exemplified by N. meningitides alanyl aminopeptidase inhibition study. In its close vicinity, a Lys residue is found in HsERAP1 and HsERAP2, Lys380 and Lys397, respectively, with their side chains extending in the cavity (residue in black bold in Figure S6). Developing interactions with this particular residue via para- or meta-substituted phenyl ring at position 4 of the aminobenzosuberone core should improve the selectivity and affinity for both members of this oxytocinase family.
All this information clarifies the inhibition observed for HsAPN and its microbial orthologues, EcPepN and PfAM1. Despite PfAM1 S1′ cavity volume being slightly larger than HsAPN and EcPepN (Table S4), the S1′ subsite cavity width of PfAM1 is shorter, leading to a lack of binding affinity for the bulky ferrocenyl substituent of 21d which cannot fit appropriately the PfAM1 S1′ subsite (Ki > 50 µM). With regard to the SAR experimentally determined, it highlights the importance of the substitution in position-4 of the aminobenzosuberone core, which, according to our model, almost fully occupies this subsite in the evaluated enzymes. Most aminopeptidases can accommodate bromine or a phenyl group with an improved affinity of the inhibitor, when compared to the unsubstituted derivative. Nevertheless, a phenyl group is mostly preferred probably due to a possible π-stacking with EcPepN His297, leading to the low Ki value of 21c when compared to 21a or 21b.
3.4. Intra/Inter-Domain Movements
An important feature to be considered for M1 aminopeptidases molecular recognition is their inter- or intra-domain movements. Various conformations have been captured during crystallization processes and open, intermediate and closed conformations have been determined for various given enzymes, suggesting that these dynamics are important for ligand binding and peptide hydrolysis. A large inter-domains motion has been markedly observed for HsERAP1, associated with local changes for the highly conserved Tyr residue (Tyr438 for HsERAP1) in the active site (corresponding to the position Tyr381 of EcPepN in Figure S1). In an “open” form of the protein, this residue undergoes a 6 Å shift away from the zinc ion, whereas in a closed form, this Tyr residue is proposed to both assist in the activation of the peptide bond and in the stabilization of the tetrahedral intermediate during the proteolytic catalytic cycle. It has been shown, by site directed mutagenesis, that mutation of this Tyr residue into Phe completely abolished the enzyme activity.
As a result of this dynamic, the open form of the enzyme is likely catalytically inactive and a subsequent allosteric activation upon substrate binding might trigger a catalytically active closed form. This dynamic transition between those different states might help the substrates/products to bind to or to be released from the active site. In solution, a number of conformers might exist, from open to closed states, along with different position of the Tyr residue. We have previously suggested that this Tyr residue probably activates the carbonyl function of the aminobenzosuberone towards a nucleophilic attack of the catalytic water to form a hydrated form of this ketone function and generates in situ a potent chelator of the zinc ion and a perfect mimic of the transition state. So, our low-molecular weight reversible inhibitors would be likely to bind to a closed form of HsERAP1 for optimum interactions. Interestingly, HsIRAP was reported to undergo a similar conformational change with, this time, complete different orientation of the GAMEN motif upon inhibitor binding, leading to a closed conformation of HsIRAP (PDB 5MJ6) and highlighting some significant structural plasticity of the GAMEN loop. Since those enzymes in solution should exist in a conformational equilibrium between different forms (even so small angle X-ray scattering experiments have suggested that the HsERAP1 average structure in solution adopts an open conformation), these processes could account for the moderate inhibition activity values observed with the HsERAP sub-family (IC50 = 1.6 to 150 µM) and, to a lesser extent, HsIRAP (IC50 = 0.034 to 2.6 µM), indicating that the aminobenzosuberone scaffolds are not able to induce domain movements to a more closed form of the enzyme.
Intra-domain mobility is also to be considered. For example, besides a more compact folding, a particularity of the human and porcine APN structures determined up to now is the plasticity of the S1 subsite through conformational changes. It looks like a small/shallow hydrophobic pocket, compared to the deeper S1 subsites observed in other crystallized M1 aminopeptidases. Indeed, HsAPN Phe896 which serves to cap the S1 pocket belongs to a flexible loop (as indicated by an average temperature factor of 25–35 Å2) of 8 residues (Y891GGGSFSF898) in domain IV and might act as a gatekeeper restricting access and closing the S1 subsite. Upon ligand binding, this S1 pocket might undergo a major conformational change via movement of the flexible 891–898 loop, revealing thus a deeper cylindrical hydrophobic S1 subsite. The cavity volumes of the active site of HsAPN varied from 572.12 Å3 to 811.87 Å3 (Figure S8 and Table S2). Our series of compounds are very potent on this human enzyme, probably through interactions with the “Phe-In” conformation of this loop via the bromine substituent in position-1 of the aminobenzosuberone core and the aromatic ring of Phe896 (8i, Ki = 60 pM). These interactions might induce or stabilize the closure of domain IV, especially that of mammalian APN and hence a closed APN form, which is an important feature to block porcine coronavirus cell entry via its porcine APN receptor.
These different conformational dynamic processes seem to be an important aspect of catalytic regulation for some M1 enzymes, although it has yet to be more documented with a larger number of structures. These inter- or intra-domain movements are probably only one remarkable feature of aminopeptidase M1 dynamics, according to their biological functions and their cellular localization.
3.5. ADME-Tox Properties
To gain insight into the physicochemical and ADME-Tox properties of these aminobenzosuberone derivatives, several parameters were experimentally determined or computed (Table S5). The experimental logD7.4 values of the whole aminobenzosuberone series complied with Lipinski’s “rule of five” and the calculated TPSA value is comprised between 44.71 and 68.10 Å2 for the ketone and hydrate forms, respectively. These compounds are predicted to penetrate the blood-brain barrier, although some possible P-glycoprotein (P-gp) efflux issues should be taken into consideration (predicted probability to be P-gp substrate is 0.5). They could be considered as valuable chemical tools to target M1 aminopeptidases implicated in CNS diseases, mostly involved in the brain renin-angiotensin system (HsIRAP, HsAPN, HsAPA) and in the pain management (HsAPN). However, the aminobenzosuberone scaffold requires careful optimization and increasing selectivity toward targeted M1 aminopeptidase versus HsPSA, which plays major neuroprotective roles. Analysis of the calculated probabilities of drug-drug interaction with different CYP450 isozymes shows that this series of molecules should inhibit CYP1A2 and to a lesser extent CYP2D6. In addition, the probability of toxicity is high and is mostly associated with cardiotoxicity with the inhibition of hERG channel and also hepatotoxicity. The uncertainties in the predicted outcomes should be taken into account and these potential toxicity liabilities should be investigated further. Thus, considering ADMET properties and inhibitory activities, the most promising compound for subsequent studies is the 1,4-dibromo derivative 21h, which is also the most synthetically feasible molecule. However, according to the predictive method, decreasing the lipholicity (clogP ou logD7.4) by introducing heteroatoms into the phenyl ring of compound 21c might also be a direction of modification since it should improve the ADME properties and have beneficial effects on toxicity.
4.1. General Information
Reactants were purchased from usual provider. Usual solvents were freshly distilled, dry MeOH distilled over Mg/MgI2, dry DME over Na and benzophenone, dry Et2O was distilled and stored over Na, DCM was distilled over P2O5 and stored over anhydrous K2CO3.
Flash chromatography: silica gel (Merck 60, 230–400 mesh). TLC: Al-roll silica gel (Merck 60, F254). Mp: Kofler hot bench, corrected. 1H- and 13C-NMR (400 MHz and 100.6 MHz resp.) spectra: Bruker Avance 400; δ in ppm and J in Hertz.
All these racemic substituted 7-amino-5,7,8,9-tetrahydrobenzocyclohepten-6-one hydrochloride salts were already fully described in literature and were resynthesized to determine their selectivity profiles. Their inhibitory activities against pAPN, hLTA4H and AAP were already reported.
4.2. General Procedure for Rubottom Oxidation
To an ice-cold mixture of water and acetone (70:35 mL) were added NaHCO3 (6.18 g, 20 eq.) and Oxone® (11.3 g, 5 eq.). The suspension was stirred at 0 °C for 30 min and then a solution of silyl enol ether 15 (1.3 g, 1 eq.) in DCM (70 mL) was dropwise added. The mixture was warmed to r.t. and stirred for 3 h (TLC monitoring). Layers were separated and aqueous layer was extracted with DCM. Combined organic layers were washed with brine, dried on MgSO4, filtered and concentrated to give silyl-oxy ketone 16, which was used without further purification.
4.3. General Procedure for Oxime Reduction
To a solution of hydroxy-oxime 17 (220 mg, 1 eq.) in methanol (13 mL) was added CoCl2·6H2O (388 mg, 2 eq.). The mixture was cooled to −30 °C and NaBH4 (462 mg, 15 eq.) was carefully added. The reaction was slowly warmed to r.t. and stirred for 2 h (TLC monitoring). The mixture was diluted with water and extracted by AcOEt. Organic layers were washed with brine, dried on MgSO4, filtered and concentrated to give amino-alcohol 18, which was N-protected without further purification.
4.4.1. PfAM1 and PfAM17
Genes encoding residues 192–1085 of native Plasmodium falciparum alanyl aminopeptidase PfAM1 (PlasmoDB PF3D7_1311800) and residues 84–605 of native Plasmodium falciparum leucyl aminopeptidase PfAM17 (PlasmoDB PF3D7_1446200) were synthetized by Genecust (Luxembourg), with the help of algorithms developed for optimizing, DNA sequences to improve proteins expression in Escherichia coli. iSynthetic genes were cloned into the T7 expression pET45b(+) vector (Novagen), which appended a N-terminal hexahistidine tag (KnpI and SalI sites for PfA-M1–BamHI and SalI sites for PfA-M17).
Escherichia coli Rosetta 2 (DE3) bacteria (Novagen) were transformed with these recombinant plasmids, after their validation by Sanger sequencing (Beckman Coulter Genomics). Bacterial cultures were grown in auto-induced LB medium (Merck) supplemented with carbenicillin (50 µg/mL) and chloramphenicol (34 µg/mL), during 24 h at 25 °C, prior to bacterial extract preparations with BugBusterTM (Novagen, Darmstadt, Germany). The clarified lysates were loaded onto Ni2-charged HisTrap column (GE Healthcare) equilibrated in 20 mM imidazole phosphate buffer and washed in the same buffer. Bound recombinant proteins were then eluted in 80 mM imidazole for PfAM1 and 200 mM imidazole for PfAM17, in phosphate buffer. Eluted fractions were finally purified by size exclusion chromatography on a Superdex 200 10 300 (equilibrated with Tris HCl 50 mM, NaCl 200 mM, ZnCl2 10 µM, pH 7.4 for PfAM1 and with Hepes 50 mM, NaCl 300 mM, ZnCl2 10 µM, 5% glycerol, pH 8.5 for PfAM17) using an Äkta purifier chromatography system (GE Healthcare Life Science, Little Chalfont, England).
4.4.2. HsERAP1, HsERAP2 and HsIRAP
The expression and purification of recombinant human endoplasmic reticulum aminopeptidase 1 (HsERAP1), endoplasmic reticulum aminopeptidase 2 (HsERAP2) and insulin-regulated aminopeptidase (HsIRAP) have been described before. Briefly, recombinant baculovirus containing the gene for each enzyme was produced in sf9 cells according to the manufacturer’s instructions (Bac-to-Bac baculovirus expression system, Invitrogen). Recombinant proteins were expressed in Hi5 cells after infection with the appropriate recombinant baculovirus and purified by NiNTA affinity chromatography as previously described (PMID: 21314638). Proteins were aliquoted and stored at −80 °C in a buffer containing 10 mM Hepes pH 7.0, 100 mM NaCl, 10% glycerol, until needed.
4.4.3. EcPepN and HsPSA
Both the enzymes were purified as reported earlier.
4.5.1. PfAM1 and PfAM17
Tests were carried out on HP/Agilent UV-Visible diode array spectrophotometer 8453 (HP/Agilent, Santa Clara, CA, USA), at 30 °C in Tris HCl 50 mM pH 7.4 for rPfAM1 and pH 8 for rPfAM17, with a final DMSO concentration of 1%. Enzyme activities of rPfAM1 (10 nM) and rPfAM17 (60 nM) were determined by continuously measuring the release of para-nitroaniline at 405 nm with alanine para-nitroanilide as substrate for rPfAM1 (Km = 1.5 mM) and leucine para-nitroanilide for rPfAM17 (Km = 0.5 mM). After 10 min of incubation, CI50 measurements were conducted during 30 min for rPfAM1 and 3 h for rPfAM17. Initial velocities were measured with increasing concentrations of inhibitors and Ki values were determined by Dixon plots in triplicate.
5 µg of rLTA4H (Sigma) were used to measure the activity of this protein through the release of alanine para-nitroanilide at 405 nm (Km = 2 mM). Experiments were performed on HP/Agilent UV-Visible diode array spectrophotometer 8453 at 30 °C, during 30 min, in Tris HCl 10 mM, KCl 0.1 mM, pH 7.5 to obtain IC50 values.
4.5.3. EcPepN and HsPSA
Enzyme activity assays were carried out in 10 mM Tris-HCl buffer at pH 7.5 with 150 mM NaCl, enzyme, 100 µM substrate (Leu-pNA) in a 100 μL reaction at 37 °C. The absorbance was read at 405 nm. Concentrations ranging from 1 nM to 100 μM of inhibitors (dissolved in DMSO) were incubated with each of the enzymes (75 µM of EcPepN, 120 µM of HsPSA) for 30 min at 37 °C in a 96 well plate. After addition of the substrate, immediate increase in absorption was measured. All inhibition assays were performed in triplicates and initial velocities were measured. Percentage inhibition and log values of inhibitor concentrations were plotted using SigmaPlot to determine IC50 values. Ki values were determined by the Dixon method. The assay consists of reaction buffer, enzyme and varied concentrations of Leu-pNA (25, 50, 100, 200 and 400 µM). Based on the IC50 values, at each substrate concentration, inhibitor concentration was varied in the range of 50 nM to 500 nM or 0.5 µM to 25 µM. Data were plotted as 1/rate versus inhibitor concentration for each substrate concentration and a linear fit was calculated by non-linear regression using SigmaPlot12.5. All the reactions were performed in triplicates, SD values were reported.
4.5.4. HsERAP1, HsERAP2 and HsIRAP
The enzymatic activity of HsERAP1, HsERAP2 and HsIRAP was determined by following the time-dependent increase in fluorescence at 460 nm (excitation: 380 nm) of the fluorigenic substrates L-Leucine-7-amido-4-methyl coumarin (L-AMC; Sigma, Darmstatd, Germany) for HsERAP1 and HsIRAP and L-arginyl-7-amido-4-methyl coumarin (R-AMC; Sigma) for HsERAP2. Measurements were performed on a 96-well plate (150 μL total volume in each well) on a TECAN infinite M200 microplate fluorescence reader. 30 nM of HsERAP1, 6 nM of HsERAP2 or 6 nM of HsIRAP was added in each well, along with 50 μM of substrate and varied concentrations of compound (within the range of 1 nM to 200 μM). The reaction was followed for 10 min at room temperature and the remaining enzymatic activity was calculated by measuring the slope of the time course. For calculation of the in vitro IC50 values, experimental data were fit to the following equation using the GraphPad Prism software package: Y = Bottom + (Top − Bottom)/(1 + 10^((LogIC50 − X) * HillSlope))(1)
where Y is the enzymatic activity and X the inhibitor concentration.
4.6. Measurement of Octanol/Water Partition Coefficient (logD7.4)
The lipophilicity (logD7.4 values) of the amino-benzosuberones was determined by reversed-phase HPLC on a C-18 column (Eclipse XDB-C18, 5 µm, 4.6 × 150 mm, from Agilent, France) according to the method previously described by Minick and Pomper. Measurement of the chromatographic capacity factors (k′) for each compound was done at various concentrations in the range of 90–60% methanol (containing 0.25% octanol) and an aqueous phase consisting of 0.15% n-decylamine in 0.02 M MOPS (3-(N-morpholino)propanesulfonic acid) buffer pH 7.4 (prepared in 1-octanol-satured water). The flow was 0.9 mL/min. The aminobenzosuberone derivatives were dissolved at 0.4 µM in methanol. The standards were also dissolved in methanol. Column void volume was estimated from the retention time of uracil, which was included as a non-retained internal reference with each injection. The capacity factors (k′) are extrapolated to 100% of the aqueous component given the value of k′w. The logD7.4 is then obtained by the formula:logD7.4 = 0.13418 + 0.98452k′w(2)
Positive control with β-oestradiol, experimental logD7.4 3.043 ± 0.093 (n = 9), literature value 3.261.
4.7. In Silico Prediction of ADMET Properties
The study of ADMET properties was carried out on the website http://admet.scbdd.com.
M1 family aminopeptidases have broad and overlapping substrate specificity; hence small-molecule inhibitors may not always be specific, especially regarding their selectivity toward bimetallic enzymes. The aminobenzosuberone scaffold demonstrated exclusive selectivity for the monometallic M1 aminopeptidase family with particular potent inhibitory activities against mammalian APN and its microbial orthologues EcPepN and PfAM1 (Ki values in the nanomolar range). It retains interesting potency on HsPSA and also on HsIRAP (Ki values in the nano- to sub-micromolar range) but it is seemingly less active on both oxytocinase enzymes HsERAP1/2 (IC50 values in the micromolar range) and totally inactive against HsLTA4H.
This promising scaffold represents an attractive new class of potent antimalarial compounds targeting PfAM1, a crucial enzyme involved in at least the last committed step of haemoglobin degradation by Plasmodium falciparum. In addition, growing evidences highlight the crucial role played by mammalian aminopeptidases in a great variety of cancer types, especially HsAPN and HsERAP1/2 via their proteolytic activity or their ability to modulate protein-protein interactions. The aminobenzosuberone core is an attractive starting point to design triple inhibitor of all these three enzymes, acting by interfering with endothelial cell morphogenesis and cell motility and by modulating antigen processing to trigger cancer immunotherapy.
The next challenge is to rationally design selective inhibitors for individual M1 aminopeptidase, to study their biological roles and precise their functions, or to avoid any potential adverse effects following in vivo treatment with aminobenzosuberone derivatives by targeting other members of this diverse family. To achieve this goal, the design process efforts should take into account the plasticity of the active site and the conformational dynamics of these M1 aminopeptidases. Interesting hints suggest deeper interactions into the S1 subsite through a substitution on position-9 of our scaffold should offer new opportunities to improve both activity and selectivity. Another approach for achieving selectivity is to look for the cellular/subcellular localization and/or tissue distribution of the targeted aminopeptidase, to design targeted-prodrug to improve site-specific drug delivery. These different strategies are currently under investigation and will be reported in due course.