Next Article in Journal
Discovering the RNA-Binding Proteome of Plant Leaves with an Improved RNA Interactome Capture Method
Previous Article in Journal
Selenium and Selenoproteins in Adipose Tissue Physiology and Obesity
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biomolecules
Volume 10
Issue 4
10.3390/biom10040659
biomolecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

P1′ Residue-Oriented Virtual Screening for Potent and Selective Phosphinic (Dehydro) Dipeptide Inhibitors of Metallo-Aminopeptidases

by
Michał Talma
and
Artur Mucha
*
Department of Bioorganic Chemistry, Wrocław University of Science and Technology, Wybrzeże Wyspiańskiego 27, 50-370 Wrocław, Poland
*
Author to whom correspondence should be addressed.
Biomolecules 2020, 10(4), 659; https://doi.org/10.3390/biom10040659
Submission received: 20 March 2020 / Revised: 17 April 2020 / Accepted: 20 April 2020 / Published: 24 April 2020
Browse Figures
Versions Notes

Abstract

:
Designing side chain substituents complementary to enzyme binding pockets is of great importance in the construction of potent and selective phosphinic dipeptide inhibitors of metallo-aminopeptidases. Proper structure selection makes inhibitor construction more economic, as the development process typically consists of multiple iterative preparation/bioassay steps. On the basis of these principles, using noncomplex computation and modeling methodologies, we comprehensively screened 900 commercial precursors of the P1′ residues of phosphinic dipeptide and dehydrodipeptide analogs to identify the most promising ligands of 52 metallo-dependent aminopeptidases with known crystal structures. The results revealed several nonproteinogenic residues with an improved energy of binding compared with the best known inhibitors. The data are discussed taking into account the selectivity and stereochemical implications of the enzymes. Using this approach, we were able to identify nontrivial structural elements substituting the recognized phosphinic peptidomimetic scaffold of metallo-aminopeptidase inhibitors.

1. Introduction

Phosphinic dipeptides are generally acknowledged as transition-state analog inhibitors of metallo-aminopeptidases [1,2,3], a large and homogenous class of exopeptidases that specifically cleave the N-terminal amino acid residues of polypeptides and proteins and contain metal ion(s) in their active sites. Being ubiquitously distributed in subcellular organelles, cytoplasm, and cellular membranes throughout all of the biological kingdoms, these enzymes perform critical roles in physiological processes and diseases. In humans, they are responsible for processing peptides and proteins involved in, for example, angiogenesis, antigen presentation, cancer and metastasis, blood pressure regulation, and bacterial and protozoal infections [3,4,5]. The inhibitory potential of phosphinic dipeptides originates from two fundamental structural features that determine the competitive mechanism of binding, and govern the potency and specificity of the ligands. First, the presence of the central phosphinic functional group ensures steric and electronic resemblance to the scissile amide bond of hydrolyzed substrates in the tetrahedral transition state. This also involves the ability to coordinate with catalytic metal ion(s). Second, the optimized structures of the P1 and P1′ side chains guarantee complementarity to the corresponding S1 and S1′ binding pockets (encryptions according to the nomenclature proposed by Schechter and Berger to define the N-terminal, non-primed and C-terminal, primed residues of the scissile amide bond in peptide substrates, and the corresponding binding pockets [6]).
Mammalian, protozoan, and bacterial enzymes of the M01 and M17 families have probably been the most comprehensively studied metallo-aminopeptidases with phosphinic dipeptides [2,5]. Alanine aminopeptidase (aminopeptidase N, APN, M01.001, EC 3.4.11.2) and leucine aminopeptidase (LAP, M17.001, EC 3.4.11.1) can be considered as the prototypical representatives. Importantly, the greatest results of inhibition were achieved for compounds comprising side chains originating from nonproteinogenic amino acids. Such fragments of optimal ligands are not trivial and can hardly be predicted by directly translating the results of studies on peptidyl substrates. These residues can be conveniently designed with the aid of molecular modeling methods; however, this process may demand a few iterations of design, synthesis, and evaluation. For example, it took us ten years to inhibit (Ki) mammalian APN and LAP by phosphinic dipeptides at a low nanomolar/subnanomolar level [7,8].
The synthesis of each individually modified metallo-aminopeptidase inhibitor requires the development of multistep synthetic pathways with orthogonal protective strategies, leading to fundamental building blocks. In the most typical procedure, P1-specific N-protected aminoalkyl-H-phosphinic acid and P1′-specific α-substituted acrylate are such components (Scheme 1) [1]. Aldehydes, whose structures encode the P1 and P1′ residues, are primary substrates for the preparation of these synthons. Thus, the proper and reliable selection of the fundamental starting materials is of invaluable significance as it results in either fruitful or vain efforts.
In the current work, we address the design issue to a virtual screening approach. In particular, the P1′ residues of active and specific phosphinic dipeptide inhibitors of metallo-aminopeptidases were targeted. A total of 900 commercially available aldehydes were selected as potential precursors of these P1′ residues, while the P1 fragment (mimicking homophenylalanine) remained conserved. Docking analysis was performed with 52 crystal structures of 27 enzymes. In addition, P1′-dehydrodipeptide analogs were included in the study, as the aldehydes may also be the substrates of the P1′-unsaturated counterparts. One of the methods for their preparation involves substitution of P1-specific N-protected aminoalkyl-H-phosphinic acids with Morita–Baylis–Hillman acetates (Scheme 1) [9]. The latter are obtained by acetylation of allyl alcohols, which, in turn, are synthesized from aldehydes and alkyl acrylates. Similar to the saturated compounds, phosphinic dehydrodipeptides were reported as inhibitors of metallo-aminopeptidases, although to a lesser extent [10].
As the main outcome of this study, we defined a set of novel nonproteinogenic P1′ substituents that are potential structural fragments of active pseudodipeptide inhibitors of metallo-aminopeptidases. The details of interactions with the target proteins were analyzed at the molecular level by docking of the individual stereoisomers in the enzyme active sites. Stereochemical considerations allowed implementation of favorable configurations of the α-carbon atoms and geometry of the double bonds in the optimized structures.

2. Materials and Methods

2.1. Target Selection and Preparation

The metallo-aminopeptidases from the MEROPS database were selected for this study (ebi.ac.uk/merops/) [11]. The fundamental criterion was an available crystal structure in the Research Collaboratory for Structural Bioinformatics Protein Data Bank (RCSB-PDB, rcsb.org) [12]. For the enzymes with more than one deposited structure for a particular species, the structure with the highest resolutions was chosen. Ten of 62 crystal structures downloaded from the PDB lacking either metal ion(s) in the active centers or putative peptidases were excluded from the study. The remaining 52 enzymes are listed in Table 1. Using options implemented in Protein Preparation Wizard the structures were prepared for calculation: all of the water molecules and the ligands were removed, and the proteins were protonated at the experimental pH [13].

2.2. Ligand Preparation

The first step of structural optimization of phosphinic dipeptide and dehydrodipeptide inhibitors (Figure 1) involved selection of aldehydes as the building blocks coding the R. The set of commercially available carbonyl compounds was taken from the Sigma-Aldrich/Merck online catalog. The full list of the aldehydes in this study is included in the Supplementary Materials. LigPrep [66] was used to create the stereoisomers of the final compounds and to pronate them at the experimental pH. Epik and OPLS3e force fields were used for the geometry optimization, and the metal binding states were added to the calculations [13]. The procedure was also applied to known inhibitors, which were considered as the reference compounds.

2.3. Docking Studies

All of the stereomeric compounds were docked (Figure 2) in the grid created for the selected enzymes (900 × 4 × 52 experiments). Receptor Grid Generation used the centroid position of metal or metals for site creation. The site of the enclosing box was set at 20 Ȧ (36 Ȧ for PDB: 2EK8), and all of the rotatable groups of amino acids in the range were selected. Other settings were set to default, and no constraints or excluded volumes were added. The ligands were docked by Glide-HTVS (High Throughput Virtual Screening algorithm) with flexible ligands, nitrogen inversions and ring conformation sampling. Bias sampling of torsions was performed for the amides only. The Epik state penalties were added to the docking score. Five ligand poses were included in scoring optimization, but the pose with the highest score was written out.
The top 100 inhibitors with the highest Glide-HTVS docking score were selected for the next step: redocking with the Glide-SP (standard precision) algorithm. The settings were kept the same as in the previous step. The top 10 inhibitors interacting with metal ion(s) by the phosphinic group were conservatively used in the next step. The inhibitors that did not interact with metal(s), but received relatively high scores were also included (maximum of 10). Thus, the number of inhibitors ranged from 10 (only phosphinic) to 20 (all top 10 metal-interacting inhibitors with relatively low scores compared with the top 10 nonmetal-interacting inhibitors). The penultimate stage was induced fit docking (IFD) [67]. The box center was set on the center of the metal or metals with size of 20 Ȧ to keep the size similar to the grid from the first step (36 Ȧ for PDB: 2EK8). The VSGB (variable-dielectric generalized born) model, which incorporates residue-dependent effects, was used. The solvent was water. The side chains were optimized within 5.0 Ȧ of ligand poses, and Glide redocking was carried out with the XP (extra precision) algorithm. The top pose for each ligand was saved. The last steps, rePrime refinement [67] and MM-GBSA (molecular mechanics-generalized born surface area) calculations, were performed to calculate the Gibbs free energies with protein flexibility, and the distance from the ligand was also set to 5.0 Ȧ.

2.4. ADMET Studies

ADMET was performed online with SwissADME [68] online tools. The following properties were selected for publication: molecular weight (MW), lipophilicity score (logP), water solubility (WS), number of possible H-bond donors (Hdon) and acceptors (Hacc), Lipinski’s properties for drug likeness, and gastrointestinal absorption (GI).

3. Results and Discussion

As the main outcome of the screening, the top results in three categories are presented and discussed. The optimized structures emerged based on the lowest calculated Gibbs free energy of ligand–protein complexes, taking into account compound stereochemistry. The first category covered potentially the most active phosphinic (dehydro)dipeptide inhibitors. In the cases of enzymes for which the most favorable ligand did not interact with the metal ion(s), the best metal-complexing inhibitors were also selected. The structures that were specifically bound were also proposed as the third category. The specificity reflected with the particular P1′ substituent presents exclusively among the top 10 inhibitors of a single aminopeptidase.
The binding energies obtained for virtual structures were compared with the values calculated for potent known inhibitors that were reported in the literature, and the most frequent inhibitor was bestatin. Only selected data are included in the body of the main manuscript. The M1 family, the most numerous and intensively studied family, is presented in detail. Then, the results for arbitrarily chosen enzymes important for medicinal chemistry applications are also given. The full list of preferred inhibitors for each enzyme is given in the Supplementary Materials (Table S1).

3.1. M1 Family

M1 is a family of metallo-aminopeptidases containing the ‘HEXXH-(X18)-E’ zinc-binding motif and ‘GXMEN’ exopeptidase motif in their sequence and one zinc cation in their active site [3,69]. The family contains enzymes widely occurring in living organisms, including alanyl, cystinyl, and lysyl aminopeptidases; leukotriene A4 hydrolase; and endoplasmic reticulum aminopeptidases I and II [70]. These enzymes play different physiological roles and are also targets for treating various pathological states. Nine M1 aminopeptidases have been identified in humans, six of which have been associated with disorders [69]. For example, their altered activity has been connected with regulation of blood pressure controlling systems, resulting in hypertension and strokes [71]. M1 aminopeptidases are also involved in cancer development [69]. In particular, APN has been shown to be overexpressed in leukemia blasts in acute myeloid leukemia. APN is active in tumor angiogenesis and metastasis [72]. LTA4H has been associated with acute lung damage, while dysregulation of ERAP1, ERAP2, and IRAP can reflect the risk of autoimmune diseases. M1 aminopeptidases are receptors for viruses, for example, coronaviruses [73]. They are also connected with several infections caused by bacteria, for example, meningitis caused by Neisseria meningitides [74], or protozoans, for example, malaria transmitted by Plasmodium falciparum. The parasite aminopeptidase of this species is involved in the digestion of short hemoglobin-derived peptides and has been rated as a highly prospective therapeutic target [75,76].
The S1′ binding pockets of M1 aminopeptidases are available for substrates with a wide range of sizes (Table 2). Principally, these cavities are formed from functional (performing catalytic roles) and hydrophobic amino acid residues. Among the functional amino acids, each enzyme comprises tyrosine, which stabilizes the gem-diol transition state, and histidine, which participates in metal ion complexation. The carbonyl oxygen atom of the P1′ residue of the cleaved substrates interacts with the positively charged side chains of arginine and lysine and noncharged chains of asparagine and threonine or the amino group of amino acids located on the border of the S1–S1′ cavities. Acidic side chains of glutamic and aspartic acids are responsible for controlling the flexibility of the active site. The position of these amino acids may change inside or outside the cavity.
In accordance with the overall character of S1′ pockets, the best docking scores for M1 aminopeptidases were achieved for bulky aromatic and heteroaromatic P1′ substituents (Table 3). Although π-conjugated systems are naturally hydrophobic, certain modifications with heteroatoms/functional groups tighten binding with the enzyme active site by making specific contacts. In some cases of sterically extended residues, the fit appeared so close that the S1′ pocket seemed virtually “stoppered” in docking experiments. Rotational freedom in phosphinic dipeptides positively influenced the energy of binding. These analogs were more frequently represented among the top compounds compared with their P1′-unsaturated counterparts. Nevertheless, a couple of enzymes appeared to favor dehydropeptides among the dozen listed in Table 3. Regarding stereochemistry, the R,S diastereoisomers of phosphinic dipeptides prevailed, and this result is consistent with the relative configuration l,l of natural substrates. For the dehydro analogs, the E geometry was somewhat favored.
Considering the selected examples, human APN showed the minimal energy of binding with phosphinic dipeptide derived from vanillin (R = 4-hydroxy-3-methoxyphenyl) (Figure 3), yet the compound could not be considered specific. Phosphinic inhibitors with P1′ aromatic rings functionalized with hydroxyl groups have been recognized as good inhibitors of metallo-aminopeptidases [7], as OH in the para position is well positioned to form hydrogen contacts with the carboxylate of Glu418. Flexibility of the side chain directs the aromatic ring towards Val385 and His388. Stacking of the ring with imidazole is facilitated by oxygen-containing groups, which change its character by the electron-withdrawal effect. The distance of 4.04 Å to the aromatic ring of Tyr477 causes additional stabilization through hydrophobic interactions. The C-terminal carboxyl group is pointed to Gly352. The compound has a nonnatural configuration on the α carbon of the P1 fragment. This is not significant, as both arrangements gave rise to very similar inhibition potential. It is also typical to observe the 2-phenylethyl fragment of hPhe-based inhibitors; because of its high flexibility, the absolute configuration does not determine the activity [84]. The most specific compound, whose structure is based on the 5-amino-thiophen-2-yl residue, reproduces an interaction with Glu418 (involving NH2), but the heteroaromatic ring is less conveniently positioned compared with His388 (Figure S1). The P1′ fragment is brought closer to the surface of the cavity, strengthening the influence of Thr384 and Val385. To our satisfaction, the binding energy calculated for a phosphinic dipeptide that was rationally designed, synthesized, and described in our previous studies (Ki = 0.69 nM) outscored the current structures [8]. However, starting materials with the specific P1 fragment (4-aminomethylphenyl) are not commercially available, and compound development demands a complex and multistep preparation. This can be considered an exception. In general, the structure proposed here surpasses known inhibitors.
As the overall architecture of the S1′ binding sites of porcine, bacterial, and protozoal APNs is rather conserved, the advantageous P1′ residues of these peptidases appeared similar to those of the human counterpart. They all contained two oxygen atoms at the remote positions of the phenyl ring. The oxygen atoms come from either boronic acid functionality (porcine APN) or 3,4-hydroxy/methoxy groups (E. coli and P. falciparum APNs). Interestingly, the bacterial enzymes showed a preference towards the constrained structure of phosphinic dehydrodipeptides. The energy of binding for this compound was equal to a very potent known pseudotripeptidyl inhibitor [79].
Groups that may be used in the development of molecular probes, including 3-(furan-2-carbonyl)quinolin-2-yl, 6-ethyl-4-oxo-4H-chromen-3-yl, phenanthren-4-yl, and quinolin-4-yl, are well scored for particular enzymes. Thus, they have the potential to serve as substituents with two roles as inhibitors and affinity-based probes. One of the most exciting results from the binding energy calculations was achieved for a 3-(furan-2-carbonyl)quinolin-2-yl-based inhibitor of cystinyl aminopeptidase. The obtained energy value was double that calculated for a very potent ligand. The residue tightly fits to the S1′ pocket and fully fills the cavity (Figure 4). The quinoline ring is situated between two aromatic rings of Tyr549 and S1 Tyr961, which highly stabilize the structure. The furan ring is located near His464, Ile461, and Glu494 at the opposed wall of the pocket.
The quinoline system installed on the dehydropeptide scaffold was shown to be well accepted by aminopeptidase from N. meningitidis. The ring is centrally located between aromatic residues of Tyr377 and His293, and these interactions cover the opposing walls of the S1′ pocket (Figure S2). The specific biphen-2-yl inhibitor of a comparable size is shifted more to the bottom to form numerous hydrophobic contacts with Phe271, Tyr372, and Tyr377 and long-range contacts with Val290 (Figure S3).
The aminopeptidase ERAP2 also prefers spatial aromatic P1′ substituents; however, these bulky residues alter the canonical binding mode. The lowest energy inhibitor is supposed to interact with the metal ion by the C-terminal carboxyl group, not by the phosphinic moiety (Figure S4). The two aromatic rings of 2-benzyloxyphenyl surround Tyr455 owing to the flexibility of the methoxy connection. This limits the opportunity for typical complexation. Similarly, the most specific phosphinic dipeptide bearing a diphenylmethyl residue does not employ the central phosphinate to coordinate with the zinc atom (Figure S5). Strong stacking interactions between aromatic systems shift the overall ligand–protein architecture. One of the P1′ phenyl rings is inserted between tyrosines 455 and 892 to form a sandwich structure, while the other one occupies the bottom-central position of the S1′ pocket between Tyr892 and Trp363. Such a bulky group as phenanthrene can also be accepted in the P1′ position. This aromatic system substituting the dehydrodipeptide scaffold of the Z configuration forms a constrained system that perfectly fills the corresponding cavity (Figure 5). Despite a visibly higher energy compared with the abovementioned cases, these compounds preserve phosphinate–zinc complexation abilities.
Large S1′ cavities could also be effectively explored by alkylaryl substituents. tert-Butyl groups were found to be particularly advantageous for aminopeptidase from S. cerevisiae (3,5-di-tert-butylphenyl) and cold-active aminopeptidase (2-tert-butylthiophenyl). In the second case, the sulfur atom participates in the interaction with Lys586, while the alkyl wraps around Tyr410 (Figure S6). The structure of the phosphinic dipeptide based on this P1′ substituent can be promising because of the low energy calculated for the ligand–protein complex, which is accompanied by a high selectivity.
For particular aminopeptidases with relatively small S1′ binding pockets, namely aminopeptidase A and human leukotriene A4 hydrolase, ambiguous results were obtained. Aminopeptidase A prefers substrates with acidic amino acid residues, whereas the best computed structures comprise a positively charged nitrogen atom. The docking results revealed the possibility of nontypical binding. Phosphinic dipeptide with a 3-isobutylaziridin-2-yl P1′ fragment may interact with the zinc ion via a phosphinate oxygen atom, but both side chains are presumably located on the surface of the S1′ pocket (Figure S7). The most specific derivative (comprising 1-methyl-1H-imidazol) binds in a reverse direction. Although the phosphinate moiety maintains the central position, the P1′ residue occupies the S1 cavity, and vice versa. The calculated binding energies are not remarkable. For leukotriene A4 hydrolase, none of the docked inhibitors exhibit canonical metal complexation. In the case of the best inhibitors, this ability could be assigned to the boronyl group of the P1′ substituent.
In conclusion, for M1 metallo-aminopeptidases, energetically favored results were obtained for compounds with large side chains filling the active center. Interactions with the amino acids on the borders of consecutive pockets are also of great importance, in particular, for catalytic Tyr and His residues involved in ligand binding.

3.2. Human Aminopeptidases other than M1

Metallo-aminopeptidase activity in humans involves extensive catalytic machinery working with different types of peptides. It is indispensable to produce vast amounts of enzymes that can cope with this task under different conditions. The M1 metallopeptidase family is 1 of 29 formally classified families and one of only three families in which the crystal structures of human enzymes were determined (Table 1). The other two families are M18 and M24. The functions of human M18 aspartyl aminopeptidase are still poorly understood. It participates in blood pressure regulation and maintenance of proper hydroelectrolytic balance, and it has also been associated with angiogenesis in breast cancer [85]. The active center of M18 aspartyl aminopeptidase contains two metal ions and is open to the central tunnel formed of 12 subunits. This gives the protein the opportunity to bind long peptides in a number of consecutive pockets. The S1 pocket is arranged to bind asparagine or glutamine, while the S1′ pocket prefers small alkyl chain-type amino acid residues (Table 4). Human M24 aminopeptidases can be divided into two subgroups, namely methionine metallo-aminopeptidases and Xaa-Pro aminopeptidases. The two subgroups are differentiated based on the presence or absence of a common fold called the “pitta bread”. Similar to other aminopeptidases, their substrate specificity is broad [86]. The catalytic centers contain two cobalt or manganese ions. The enzymes differ in amino acid preferences and active site architecture (Table 4). For example, Xaa-Pro aminopeptidase 1 has a very small pocket specific for proline, while Xaa-Pro aminopeptidase 3, despite playing a similar function, has evolved to accept larger substrates.
The best phosphinic dipeptide inhibitors found for human M18 and M24 aminopeptidases appeared unexpectedly extended considering the limited size of the binding pockets (Table 5). The detailed molecular modeling analysis revealed that none of these compounds were able to bind in a substrate-like manner. For example, nominal P1′ substituents found for aspartyl aminopeptidase (3,4-dibenzyloxyphenyl), methionyl aminopeptidase 2 (4-borono-2,3-difluorophenyl), and aminopeptidase P3 (3-phenylaziridin-2-yl) had a tendency to occupy the S1 pockets, not the S1′ pockets, showing a distorted binding mode. This distortion involved interactions of the carboxylate group with the metal ions, and the C-terminal position was held by the phosphinate. The S1 cleft of aspartyl aminopeptidase was able to adopt unexpectedly bulky polyaromatic residues, which were surrounded by the alkyl residues of Lys358, Leu405, Val407, Pro412, Cys413, and His359 (π–π stacking) (Figure S8). Tyr381 formed a hydrogen bond with an ether oxygen atom. For aminopeptidase P3, 3-phenylaziridin interacted with the catalytic residues Pro301 and Asp331 (NH of aziridine) and His413 and Arg438 (Figure S9). For methionyl aminopeptidase 2, the ligand did not typically point towards Met384, but docked in the shallow part of the S1 cleft, close to catalytic amino acids (Figure S10). The phenyl ring and fluorine atoms interacted with His231 and His382, while the boronyl group formed numerous hydrogen bonds with Ser224, Lys249, and Gln457. In all three discussed cases, the P1 phenylethyl fragment was nonspecifically arranged.
Aromatic residues functionalized with the boronyl group appeared the most universal in the construction of the best inhibitors of M18 and M24 aminopeptidases in different categories. As mentioned, the 4-borono-2,3-difluorophenyl substituent of phosphinic dipeptide was the most energetically beneficial for methionyl aminopeptidase 2, whereas the 4-boronylphenyl group of the dedydrodipeptide scaffold was the most beneficial for methionyl aminopeptidase 1, and the 4-borono-2-fluorophenyl-modified phosphinic dipeptide was the most beneficial for aminopeptidase P1. Borono-functionalized fragments were also present among the most specific or the best metal-interacting compounds. However, boronic acid was the most frequently engaged in metal complexation (together with the carboxylate, or alone), which caused a non-classical inclusion of this portion out of the S1′ clefts. Only individual borono-inhibitors exhibited the assumed binding modes, for example, 3-borono-5-isopropoxyphenyl phosphinate with aminopeptidase P3 (Figure 6). The P1′ acidic system of the ligand formed hydrogen bonds with His424, His431, and Arg438. The isopropyl group interacted with Trp163, while the ether oxygen atom was in hydrogen bond distance to the His421 amino group.
Alterations were also observed for modes that involved metal–phosphinate coordination. 5-Ethoxymethylfuran-based phosphinic dipeptide inhibitor of aspartyl aminopeptidase and 4-(6-(hydroxymethyl)pyridinyl)phenyl-based methionyl aminopeptidase 1 showed, for instance, a reversed orientation of binding, with the P1 and P1′ positions flipped. For the first-mentioned case, the flexible N-terminal homophenylalanine fragment easily fit into the S1′ subsite, maintaining the ability of the amino group to interact with acidic residues of the zinc complexation system, in particular Glu301 and Glu302 (Figure 7). The phosphinic group coordinates with the metal cations, while the C-terminal carboxylate occupies the position of the acidic side chain of substrates preferred in the S1 pocket. The 5-ethoxymethylfuran group is extended to the S2 part and forms contacts with His349, Met439, and His440.

3.3. Aminopeptidases from Pathogenic Bacteria

Various pathological conditions can be caused by invading (micro)organisms, usually bacteria, protozoans, or viruses, that colonize humans and other species. These pathogens disrupt the body homeostasis through many mechanisms; from introducing toxins and competing for substrates to digesting peptides playing a significant role in the host organism. Microbial peptidases, including metallo-aminopeptidases, were adapted for the last-mentioned purpose [92]. They have been classified in a wide spectrum of families, but here, we focus on a few notorious examples in humans. Disabling their activity is extremely difficult because of the similarities in the active center structure and substrate specificity of microbial peptidases and human metallo-aminopeptidases. On the other hand, finding specific inhibitors offers specialized tools for highly targeted antimicrobial intervention.
The crystal structures of enzymes of microorganisms pathogenic to humans were determined in six families: M1, M17, M18, M24, M29, and M42. The M1 family (Table 3) includes, for example, the enzyme from Plasmodium falciparum, the most fatal of the four protozoan species that cause malaria in humans. This species also produces an M17 aminopeptidase [93], an aspartyl aminopeptidase classified in the M18 family [36], and the M24 aminopeptidase P [51]. Metallo-aminopeptidases are etiological factors of microbial infections, for example, meningitis caused by Neisseria meningitides, or yeast species. Accordingly, Saccharomyces cerevisiae can cause alterations in the production of B4 leukotriene and is a target in the treatment of bowel diseases, candidiasis, or psoriasis [94]. The full description of the active centers and substrate specificities of some microbial metallo-aminopeptidases is not always possible because of the lack of exhaustive data (Table 6).
M17 mixed type A/I aminopeptidase from Helicobacter pylori, in addition to urease (Nobel Prize in Medicine for Marshall and Warren in 2005 [95]), is among the potential targets for organophosphorus inhibitors. This bizinc-dependent enzyme has a wide active center arranged for peptides with long side chains [51], but the S1′ preferences have not yet been specified. Nevertheless, good docking results were obtained for compounds with boronyl groups (Table 7). The 3-boronyl-6-ethoxyphenyl substituent proved the best dehydrodipeptide analog inhibitor, with a binding mode similar to that of bestatin (Figure 8). The P1 group of hPhe is arranged in S1 and interacts with Met278 and Tyr433 of another enzyme subunit. The P1′ substituent also fills the corresponding cavity very well. Tyr274 interacts with the boronyl group by a hydrogen bond, while the aromatic ring and ethoxy group form hydrophobic interactions. Lys270 plays a role in stabilizing the transition state of peptide bond hydrolysis; accordingly, it binds with its amino group to the central phosphinate moiety.
The M17 enzymes are also characteristic of bacteria of the Pseudomonas genus, which are responsible for endogenous infections of individuals with reduced immunity, including pneumonia, among other microorganisms that express metallo-aminopeptidases and are involved in respiratory diseases, for example, Pseudomonas aeruginosa (M18), Mycobacterium tuberculosis (M24), or Streptococcus pneumoniae (M42). The best pseudodipeptide that interacts with metal ions of glutamyl M42 aminopeptidase via phosphinate functionality is based on the 4-(difluoromethoxy)-3-hydroxyphenyl P1′ fragment. Its orientation is reversed in the active center, similar to the aspartyl-specific enzyme (Figure 9). The C-terminal carboxylate is exposed to Arg257 and Ser238. The P1′ aromatic ring fits perfectly in the S1 pocket by multiple hydrophobic interactions with Pro239, Ile317, and catalytic His318. The docking is strengthened by additional halogen bonding of carbonyl Gly90.
The situation for methionine aminopeptidase 1 from Staphylococcus aureus is complex as three cobalt ions in the active center are responsible for catalytic functions. This bacterium is responsible for purulent infections of the skin, subcutaneous and soft tissues, systemic infections, and poisoning associated with the production of specific toxins. The ΔG binding energies of predicted phosphinic dipeptide inhibitors are significantly lower than those of known ligands. For metal–phosphinate-interacting compounds, the large catalytic center allows both acidic groups to point towards the three cobalt regions (Figure 10). The spatial hydrophobic S1 pocket shows a preference towards the large, heteroaromatic system. Both aromatic rings of the ligand interact with abundant histidine residues.

4. Conclusions

The study demonstrated a noncomplex, practical, and inexpensive method of determining favorable P1′ substituents in the structure of phosphinic dipeptide inhibitors of metallo-aminopeptidases. The approach involved virtual screening of the inhibitors with commonly available molecular modeling tools. The optimized residues were not generated computationally, but were taken from the structures of commercially available synthetic precursors. This supposition should ensure the feasibility of the subsequent synthesis. The obtained dataset was analyzed considering stereochemical implications and selectivity of interactions in the enzyme active sites. The most potent, the most potent metal-interacting, and the most selective diastereomeric inhibitors that were calculated (by the lowest free binding energy) were thus identified. The structures of the optimized substituents appeared not evident and, quite surprisingly, not similar to those of proteinogenic residues. Two types of these structures were predominantly present. The first type included spatially extended P1′ fragments of aromatic character, principally hydrophobic, such as aryl, heteroaryl, and arylakil. These residues were characterized by extremely tight insertion into the large S1′ pockets. Substituents of a moderate size, but modified with two oxygen atom-containing group(s), for example, boronyl and dihydroxyl, represented the second most abundant group. Oxygen-comprising groups tended to form a favorable net of hydrogen bonds with the S1′ residues of certain metallo-aminopeptidases. However, the presence of alternatives to phosphinate metal binding groups, such as boronyl and C-terminal carboxyl groups, gave rise to alterations in the typical substrate-like binding mode. Such changes were the most frequent for aminopeptidases with small substrate-oriented S1′ subsites. Reorientation of the overall direction of binding was also observed in some cases. Enriched by such binding mode details, the method can be recommended as general to design specified residues of inhibitors of peptidomimetic structures.
Phosphinic dipeptides do not show toxicity [5] and meet the basic rules that would make them orally active drugs in humans. Selected ADMET properties predicted for the inhibitors of M1 aminopeptidases validate their druglikeness (Table 8, the full table is included in Supplementary Information, Table S2). Lipinski’s rule of five is fulfilled, as the compounds are relatively small (MW only in rare cases exceeds 500) and moderately lipophilic (logP from −1.6 to 3.3). The pseudopeptidic phosphinic backbone is a source of four hydrogen bond acceptors and one donor. These basic numbers may be increased for the heteroatoms comprised by the substituents. Water solubility and gastrointestinal absorption are variable. In summary, in addition to the biological activity, phosphinic dipeptide analogs possess the physicochemical properties needed for drug candidates, opening opportunities to use phosphinic dipeptide analogs in the treatment of a spectrum of diseases.
The results of our screening can be considered as qualitative. Globally, there is no significant correlation between experimental affinities of the described inhibitors and the corresponding values of the computed free Gibbs energy. Too many variables influence the reported data based on experimental procedures (e.g., different protocols of inhibition measurements and calculations) and the data based on calculations (structural, stereochemical, and conformational implications, for both enzymes and ligands). Nevertheless, regressions that involved limited datasets showed clearer tendencies. For the first relationship, we selected the universal inhibitor (bestatin), which was assayed with several aminopeptidases, and prepared the plot of ΔG versus pKi (Table 9). Although the correlation coefficient for the linear regression is not remarkable (R2 = 0.65), the trend is visible. Much more pronounced correlation (R2 = 0.91) was envisaged for the specified type of inhibitors (phosphinic (dehydro)dipeptide analogs) and a single aminopeptidase (human aminopeptidase N, M01.001, Figure 11). Although these compounds were tested as a diastereomeric mixture [8,10], while the binding energies were calculated for the best single diastereoisomers, the obtained fit to linear regression is remarkable. Thus, the applied procedure can be recommended for inhibitors of a particular structure.

Supplementary Materials

The following are available online at https://www.mdpi.com/2218-273X/10/4/659/s1, List of commercially available aldehydes considered in the study; Table S1. The full list of enzymes with preferred P1′ substituents that were the most favored, had the best interaction with the catalytic metal ion via the phosphinate group, and were the most specific; Table S2. Selected ADMET properties for the best inhibitors; Figure S1. Modeled binding mode of a phosphinic dipeptide analog containing a 5-amino-thiophen-2-yl P1′ fragment in the active center of human aminopeptidase N (PDB: 4FYT); Figure S2. Modeled binding mode of a phosphinic dipeptide analog containing a quinolin-4-yl P1′ fragment in the active center of Neisseria meningitides aminopeptidase (PDB: 2GTQ); Figure S3. Modeled binding mode of a phosphinic dipeptide analog containing quinolin-4-yl (carbon atoms are colored light green) and biphenyl-2-yl (carbon atoms are colored gray) P1′ fragments in the active center of Neisseria meningitides aminopeptidase (PDB: 2GTQ); Figure S4. Modeled binding mode of a phosphinic dipeptide analog containing a 2-benzyloxyphenyl P1′ fragment in the active center of ERAP2 aminopeptidase (PDB: 5AB0); Figure S5. Modeled binding mode of a phosphinic dehydrodipeptide analog containing a diphenylmethyl P1′ fragment in the active center of ERAP2 aminopeptidase (PDB: 5AB0); Figure S6. Modeled binding mode of a phosphinic dipeptide analog containing a 2-tert-butylthiophenyl P1′ fragment in the active center of cold-active aminopeptidase (PDB: 5AB0); Figure S7. Modeled binding mode of a phosphinic dipeptide analog containing a 3-isobutylaziridin-2-yl P1′ fragment in the active center of aminopeptidase A (PDB: 4KX7); Figure S8. Modeled binding mode of a phosphinic dipeptide analog containing a 3,4-dibenzyloxyphenyl P1′ fragment in the active center of aspartyl aminopeptidase (PDB: 4DYO); Figure S9. Modeled binding mode of a phosphinic dipeptide analog containing a 3-phenylaziridin-2-yl P1′ fragment in the active center of aminopeptidase P3 (PDB: 5 × 49); Figure S10. Modeled binding mode of a phosphinic dipeptide analog containing a 4-borono-2,3-difluorophenyl P1′ fragment in the active center of methionyl aminopeptidase 2 (PDB: 1B6A).

Author Contributions

M.T. and A.M., conceptualization. M.T., conducted theoretical studies and methodology. M.T. and A.M., analyzed the data. M.T. and A.M., wrote the original manuscript. Both authors reviewed and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

The work was financed by a statutory activity subsidy from the Polish Ministry of Science and Higher Education for the Faculty of Chemistry of Wrocław University of Science and Technology.

Acknowledgments

The authors acknowledge the Schrödinger, Inc., American multinational biopharmaceutical and technology company for providing a trial license. The Discovery Studio Visualizer was used for preparation of the figures.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Mucha, A. Synthesis and modifications of phosphinic dipeptide analogues. Molecules 2012, 17, 13530–13568. [Google Scholar] [CrossRef] [PubMed]
  2. Talma, M.; Maślanka, M.; Mucha, A. Recent developments in the synthesis and applications of phosphinic peptide analogs. Bioorg. Med. Chem. Lett. 2019, 29, 1031–1042. [Google Scholar] [CrossRef] [PubMed]
  3. Mucha, A.; Drag, M.; Dalton, J.P.; Kafarski, P. Metallo-aminopeptidase inhibitors. Biochimie 2010, 92, 1509–1529. [Google Scholar] [CrossRef] [PubMed]
  4. Hooper, N.M.; Lendeckel, U. Aminopeptidases in Biology and Disease; Red.; Springer: Boston, MA, USA, 2004; ISBN 978-1-4613-4698-2. [Google Scholar]
  5. Skinner-Adams, T.S.; Stack, C.M.; Trenholme, K.R.; Brown, C.L.; Grembecka, J.; Lowther, J.; Mucha, A.; Drag, M.; Kafarski, P.; McGowan, S.; et al. Plasmodium falciparum neutral aminopeptidases: New targets for anti-malarials. Trends Biochem. Sci. 2010, 35, 53–61. [Google Scholar] [CrossRef]
  6. Schechter, I.; Berger, A. On the size of the active site in proteases. I. Papain. Biochem. Biophys. Res. Commun. 1967, 27, 157–162. [Google Scholar] [CrossRef]
  7. Grembecka, J.; Mucha, A.; Cierpicki, T.; Kafarski, P. The most potent organophosphorus inhibitors of leucine aminopeptidase. Structure-based design, chemistry, and activity. J. Med. Chem. 2003, 46, 2641–2655. [Google Scholar] [CrossRef]
  8. Vassiliou, S.; Węglarz-Tomczak, E.; Berlicki, Ł.; Pawełczak, M.; Nocek, B.; Mulligan, R.; Joachimiak, A.; Mucha, A. Structure-guided, single-point modifications in the phosphinic dipeptide structure yield highly potent and selective inhibitors of neutral aminopeptidases. J. Med. Chem. 2014, 57, 8140–8151. [Google Scholar] [CrossRef]
  9. Talma, M.; Mucha, A. P-C bond formation in reactions of Morita-Baylis-Hillman adducts with phosphorus nucleophiles. Arkivoc 2016, 2017, 324–344. [Google Scholar]
  10. Talma, M. Phosphinic dehydrodipeptides: Diversification of the P1′ residue with the Morita–Baylis–Hillman acetates and inhibition of alanyl aminopeptidases. Int. J. Pept. Res. Ther. 2020. accepted. [Google Scholar] [CrossRef]
  11. Rawlings, N.D.; Waller, M.; Barrett, A.J.; Bateman, A. MEROPS: The database of proteolytic enzymes, their substrates and inhibitors. Nucleic Acids Res. 2014, 42, D503–D509. [Google Scholar] [CrossRef] [Green Version]
  12. Burley, S.K.; Berman, H.M.; Bhikadiya, C.; Bi, C.; Chen, L.; Di Costanzo, L.; Christie, C.; Dalenberg, K.; Duarte, J.M.; Dutta, S.; et al. RCSB Protein Data Bank: Biological macromolecular structures enabling research and education in fundamental biology, biomedicine, biotechnology and energy. Nucleic Acids Res. 2019, 47, D464–D474. [Google Scholar] [CrossRef] [Green Version]
  13. Schrödinger, L.L.C. Schrödinger Release 2018-4: Schrödinger Suite 2018-4 Protein Preparation Wizard; Epik: New York, NY, USA, 2018. [Google Scholar]
  14. Wong, A.H.; Zhou, D.; Rini, J.M. The X-ray crystal structure of human aminopeptidase N reveals a novel dimer and the basis for peptide processing. J. Biol. Chem. 2012, 287, 36804–36813. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Chen, L.; Lin, Y.L.; Peng, G.; Li, F. Structural basis for multifunctional roles of mammalian aminopeptidase N. Proc. Natl. Acad. Sci. USA 2012, 109, 17966–17971. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Yang, Y.; Liu, C.; Lin, Y.L.; Li, F. Structural insights into central hypertension regulation by human aminopeptidase A. J. Biol. Chem. 2013, 288, 25638–25645. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Stsiapanava, A.; Olsson, U.; Wan, M.; Kleinschmidt, T.; Rutishauser, D.; Zubarev, R.A.; Samuelsson, B.; Rinaldo-Matthis, A.; Haeggstrom, J.Z. Binding of Pro-Gly-Pro at the active site of leukotriene A4 hydrolase/aminopeptidase and development of an epoxide hydrolase selective inhibitor. Proc. Natl. Acad. Sci. USA 2014, 111, 4227–4232. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Addlagatta, A.; Gay, L.; Matthews, B.W. Structural basis for the unusual specificity of Escherichia coli aminopeptidase N. Biochemistry 2008, 47, 5303–5311. [Google Scholar] [CrossRef]
  19. Mpakali, A.; Saridakis, E.; Harlos, K.; Zhao, Y.; Kokkala, P.; Georgiadis, D.; Giastas, P.; Papakyriakou, A.; Stratikos, E. Ligand-induced conformational change of insulin-regulated aminopeptidase: Insights on catalytic mechanism and active site plasticity. J. Med. Chem. 2017, 60, 2963–2972. [Google Scholar] [CrossRef]
  20. Kyrieleis, O.J.P.; Goettig, P.; Kiefersauer, R.; Huber, R.; Brandstetter, H. Crystal structures of the tricorn interacting factor F3 from Thermoplasma acidophilum, a zinc aminopeptidase in three different conformations. J. Mol. Biol. 2005, 349, 787–800. [Google Scholar] [CrossRef]
  21. Mpakali, A.; Giastas, P.; Mathioudakis, N.; Mavridis, I.M.; Saridakis, E.; Stratikos, E. structural basis for antigenic peptide recognition and processing by endoplasmic reticulum (Er) aminopeptidase 2. J. Biol. Chem. 2015, 290, 26021–26032. [Google Scholar] [CrossRef] [Green Version]
  22. McGowan, S.; Porter, C.J.; Lowther, J.; Stack, C.M.; Golding, S.J.; Skinner-Adams, T.S.; Trenholme, K.R.; Teuscher, F.; Donnelly, S.M.; Grembecka, J.; et al. Structural basis for the inhibition of the essential Plasmodium falciparum M1 neutral aminopeptidase. Proc. Natl. Acad. Sci. USA 2009, 106, 2537–2542. [Google Scholar] [CrossRef] [Green Version]
  23. Bauvois, C.; Jacquamet, L.; Huston, A.L.; Borel, F.; Feller, G.; Ferrer, J.L. Crystal structure of the cold-active aminopeptidase from Colwellia psychrerythraea, a close structural homologue of the human bifunctional leukotriene A4 hydrolase. J. Biol. Chem. 2008, 283, 23315–23325. [Google Scholar] [CrossRef] [Green Version]
  24. Helgstrand, C.; Hasan, M.; Usyal, H.; Haeggstrom, J.Z.; Thunnissen, M.M.G.M. A leukotriene A(4) hydrolase-related aminopeptidase from yeast undergoes induced fit upon inhibitor binding. J. Mol. Biol. 2011, 406, 120–134. [Google Scholar] [CrossRef]
  25. Nocek, B.; Mulligan, R.; Bargassa, M.; Collart, F.; Joachimiak, A. Crystal structure of aminopeptidase N from human pathogen Neisseria meningitidis. Proteins 2007, 70, 273–279. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Strater, N.; Lipscomb, W.N. Two-metal ion mechanism of bovine lens leucine aminopeptidase: Active site solvent structure and binding mode of L-leucinal, a gem-diolate transition state analogue, by X-ray crystallography. Biochemistry 1995, 34, 14792–14800. [Google Scholar] [CrossRef] [PubMed]
  27. Duprez, K.; Scranton, M.A.; Walling, L.L.; Fan, L. Structure of tomato wound-induced leucine aminopeptidase sheds light on substrate specificity. Acta Crystallogr. Sect. D 2014, 70, 1649–1658. [Google Scholar] [CrossRef] [PubMed]
  28. Straeter, N.; Sherratt, D.J.; Colloms, S.D. X-Ray structure of aminopeptidase a from Escherichia coli and a model for the nucleoprotein complex in xer site-specific recombination. Embo J. 1999, 18, 4513–4522. [Google Scholar] [CrossRef] [Green Version]
  29. Modak, J.K.; Rut, W.; Wijeyewickrema, L.C.; Pike, R.N.; Drag, M.; Roujeinikova, A. Structural basis for substrate specificity of Helicobacter pylori M17 aminopeptidase. Biochimie 2016, 121, 60–71. [Google Scholar] [CrossRef]
  30. Zhan, C.; Patskovsky, Y.; Wengerter, B.C.; Ramagopal, U.; Milstein, S.; Vidal, M.; Almo, S.C. Crystal structure and function of Caenorhabditis elegans leucine aminopeptidase. To be published. Deposited in the RCSB PDB as 2HC9 in 2006.
  31. Kale, A.; Pijning, T.; Sonke, T.; Dijkstra, B.W.; Thunnissen, A.M. Crystal structure of the leucine aminopeptidase from Pseudomonas putida reveals the molecular basis for its enantioselectivity and broad substrate specificity. J. Mol. Biol. 2010, 398, 703–714. [Google Scholar] [CrossRef] [Green Version]
  32. Natarajan, S.; Huynh, K.-H.; Kang, L.W. Crystal structure of leucyl aminopeptidase (pepA) from Xoo0834, Xanthomonas oryzae pv. oryzae KACC10331. To be published. Deposited in the RCSB PDB as 3JRU in 2009.
  33. Chen, Y.; Farquhar, E.R.; Chance, M.R.; Palczewski, K.; Kiser, P.D. Insights into substrate specificity and metal activation of mammalian tetrahedral aspartyl aminopeptidase. J. Biol. Chem. 2012, 287, 13356–13370. [Google Scholar] [CrossRef] [Green Version]
  34. Chaikuad, A.; Pilka, E.S.; Riso, A.D.; Delft, F.V.; Kavanagh, K.L.; Venien-Bryan, C.; Oppermann, U.; Yue, W.W. Structure of human aspartyl aminopeptidase complexed with substrate analogue: Insight into catalytic mechanism, substrate specificity and M18 peptidase family. BMC Struct. Biol. 2012, 12, 14–23. [Google Scholar] [CrossRef] [Green Version]
  35. Nguyen, D.D.; Pandian, R.; Kim, D.; Ha, S.C.; Yoon, H.J.; Kim, K.S.; Yun, K.H.; Kim, J.H.; Kim, K.K. Structural and kinetic bases for the metal preference of the M18 aminopeptidase from Pseudomonas aeruginosa. Biochem. Biophys. Res. Commun. 2014, 447, 101–107. [Google Scholar] [CrossRef] [PubMed]
  36. Sivaraman, K.K.; Oellig, C.A.; Huynh, K.; Atkinson, S.C.; Poreba, M.; Perugini, M.A.; Trenholme, K.R.; Gardiner, D.L.; Salvesen, G.; Drag, M.; et al. X-ray crystal structure and specificity of the Plasmodium falciparum malaria aminopeptidase PfM18AAP. J. Mol. Biol. 2012, 422, 495–507. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Min, T.; Shapiro, L. Crystal structure of aminopeptidase I from Clostridium acetobutylicum. To be published. Deposited in the RCSB PDB as 2GLJ in 2006.
  38. Min, T.; Shapiro, L. Crystal structure of Aminopeptidase (M18 family) from Thermotoga maritima. To be published. Deposited in the RCSB PDB as 2GLF in 2006.
  39. Evdokimov, A.G.; Pokross, M.; Walter, R.L.; Mekel, M.; Barnett, B.L.; Amburgey, J.; Seibel, W.L.; Soper, S.J.; Djung, J.F.; Fairweather, N.; et al. Serendipitous discovery of novel bacterial methionine aminopeptidase inhibitors. Proteins 2007, 66, 538–546. [Google Scholar] [CrossRef] [PubMed]
  40. Addlagatta, A.; Quillin, M.L.; Omotoso, O.; Liu, J.O.; Matthews, B.W. Identification of an SH3-binding motif in a new class of methionine aminopeptidases from Mycobacterium tuberculosis suggests a mode of interaction with the ribosome. Biochemistry 2005, 44, 7166–7174. [Google Scholar] [CrossRef]
  41. Helgren, T.R.; Chen, C.; Wangtrakuldee, P.; Edwards, T.E.; Staker, B.L.; Abendroth, J.; Sankaran, B.; Housley, N.A.; Myler, P.J.; Audia, J.P.; et al. Rickettsia prowazekii methionine aminopeptidase as a promising target for the development of antibacterial agents. Bioorg. Med. Chem. 2017, 25, 813–824. [Google Scholar] [CrossRef] [Green Version]
  42. Liu, S.; Widom, J.; Kemp, C.W.; Crews, C.M.; Clardy, J. Structure of human methionine aminopeptidase-2 complexed with fumagillin. Science 1998, 282, 1324–1327. [Google Scholar] [CrossRef]
  43. Addlagatta, A.; Matthews, B.W. Structure of the angiogenesis inhibitor ovalicin bound to its noncognate target, human type 1 methionine aminopeptidase. Protein Sci. 2006, 15, 1842–1848. [Google Scholar] [CrossRef]
  44. Tahirov, T.H.; Oki, H.; Tsukihara, T.; Ogasahara, K.; Yutani, K.; Ogata, K.; Izu, Y.; Tsunasawa, S.; Kato, I. Crystal structure of methionine aminopeptidase from hyperthermophile, Pyrococcus furiosus. J. Mol. Biol. 1998, 284, 101–124. [Google Scholar] [CrossRef]
  45. Douangamath, A.; Dale, G.E.; D’Arcy, A.; Almstetter, M.; Eckl, R.; Frutos-Hoener, A.; Henkel, B.; Illgen, K.; Nerdinger, S.; Schulz, H.; et al. Crystal structures of Staphylococcus aureus methionine aminopeptidase complexed with keto heterocycle and aminoketone inhibitors reveal the formation of a tetrahedral intermediate. J. Med. Chem. 2004, 47, 1325–1328. [Google Scholar] [CrossRef]
  46. Alvarado, J.J.; Nemkal, A.; Sauder, J.M.; Russell, M.; Akiyoshi, D.E.; Shi, W.; Almo, S.C.; Weiss, L.M. Structure of a microsporidian methionine aminopeptidase type 2 complexed with fumagillin and TNP-470. Mol. Biochem. Parasitol. 2009, 168, 158–167. [Google Scholar] [CrossRef] [Green Version]
  47. Graham, S.C.; Bond, C.S.; Freeman, H.C.; Guss, J.M. Structural and functional implications of metal ion selection in aminopeptidase P, a metalloprotease with a dinuclear metal center. Biochemistry 2005, 44, 13820–13836. [Google Scholar] [CrossRef] [PubMed]
  48. Li, X.; Lou, Z.; Li, X.; Zhou, W.; Ma, M.; Cao, Y.; Geng, Y.; Bartlam, M.; Zhang, X.C.; Rao, Z. Structure of human cytosolic X-prolyl aminopeptidase: A double Mn(II)-dependent dimeric enzyme with a novel three-domain subunit. J. Biol. Chem. 2008, 283, 22858–22866. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Singh, R.; Jamdar, S.N.; Goyal, V.D.; Kumar, A.; Ghosh, B.; Makde, R.D. Structure of the human aminopeptidase XPNPEP3 and comparison of its in vitro activity with Icp55 orthologs: Insights into diverse cellular processes. J. Biol. Chem. 2017, 292, 10035–10047. [Google Scholar] [CrossRef] [Green Version]
  50. Mizutani, H.; Kunishima, N. Crystal structure of X-Pro aminopeptidase from Thermotoga maritima MSB8. To be published. Deposited in the RCSB PDB as 2ZSG in 2008.
  51. Drinkwater, N.; Sivaraman, K.K.; Bamert, R.S.; Rut, W.; Mohamed, K.; Vinh, N.B.; Scammells, P.J.; Drag, M.; McGowan, S. Structure and substrate fingerprint of aminopeptidase P from Plasmodium falciparum. Biochem. J. 2016, 473, 3189–3204. [Google Scholar] [CrossRef]
  52. Kumar, A.; Are, V.; Ghosh, B.; Jamdar, S.; Makde, R. R372A mutant of Xaa-Pro dipeptidase from Xanthomonas campestris at 1.85 Angstrom resolution. To be published. Deposited in the RCSB PDB as 5CDE in 2015.
  53. Gilboa, R.; Greenblatt, H.M.; Perach, M.; Spungin-Bialik, A.; Lessel, U.; Wohlfahrt, G.; Schomburg, D.; Blumberg, S.; Shoham, G. Interactions of Streptomyces griseus aminopeptidase with a methionine product analogue: A structural study at 1.53 A resolution. Acta Crystallogr. Sect. D 2000, 56, 551–558. [Google Scholar] [CrossRef] [PubMed]
  54. Joint Center for Structural Genomics. Crystal structure of putative peptidase (NP_348812.1) from Clostridium acetobutylicum at 2.37 A resolution. To be published. Deposited in the RCSB PDB as 3K9T in 2009.
  55. Joint Center for Structural Genomics. Crystal structure of a hypothetical Zn-dependent exopeptidase (BDI_3547) from Parabacteroides distasonis ATCC 8503 at 1.06 A resolution. To be published. Deposited in the RCSB PDB as 3TC8 in 2011.
  56. Desmarais, W.; Bienvenue, D.L.; Bzymek, K.P.; Petsko, G.A.; Ringe, D.; Holz, R.C. The high-resolution structures of the neutral and the low pH crystals of aminopeptidase from Aeromonas proteolytica. J. Biol. Inorg. Chem. 2006, 11, 398–408. [Google Scholar] [CrossRef] [PubMed]
  57. Akioka, M.; Nakano, H.; Tsujimoto, Y.; Matsui, H.; Nakatsu, T.; Kato, H.; Watanabe, K. 2EK8: aminopeptidase from aneurinibacillus sp. strain AM-1. To be published. Deposited in the RCSB PDB as 2EK8 in 2007.
  58. Odintsov, S.G.; Sabala, I.; Bourenkov, G.; Rybin, V.; Bochtler, M. Substrate access to the active sites in aminopeptidase T, a representative of a new metallopeptidase clan. J. Mol. Biol. 2005, 354, 403–412. [Google Scholar] [CrossRef]
  59. Odintsov, S.G.; Sabala, I.; Bourenkov, G.; Rybin, V.; Bochtler, M. Staphylococcus aureus aminopeptidases is a founding member of a new peptidase clan. J. Biol. Chem. 2005, 280, 27792–27799. [Google Scholar] [CrossRef] [Green Version]
  60. Badger, J.; Sauder, J.M.; Adams, J.M.; Antonysamy, S.; Bain, K.; Bergseid, M.G.; Buchanan, S.G.; Buchanan, M.D.; Batiyenko, Y.; Christopher, J.A.; et al. Structural analysis of a set of proteins resulting from a bacterial genomics project. Proteins 2005, 60, 787–796. [Google Scholar] [CrossRef]
  61. Kim, D.; San, B.H.; Moh, S.H.; Park, H.J.; Kim, D.Y.; Lee, S.; Kim, K.K. Structural basis for the substrate specificity of PepA from Streptococcus pneumoniae, a dodecameric tetrahedral protease. Biochem. Biophys. Res. Commun. 2010, 391, 431–436. [Google Scholar] [CrossRef]
  62. Schoehn, G.; Vellieux, F.M.D.; Dura, M.A.; Receveur-Brechot, V.; Fabry, C.M.S.; Ruigrok, R.W.H.; Ebel, C.; Roussel, A.; Franzetti, B. An archaeal peptidase assembles into two different quaternary structures: A tetrahedron and a giant octahedron. J. Biol. Chem. 2006, 281, 36327–36337. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Borissenko, L.; Groll, M. Crystal structure of TET protease reveals complementary protein degradation pathways in prokaryotes. J. Mol. Biol. 2005, 346, 1207–1219. [Google Scholar] [CrossRef] [PubMed]
  64. Dura, M.A.; Rosenbaum, E.; Larabi, A.; Gabel, F.; Vellieux, F.M.; Franzetti, B. The structural and biochemical characterizations of a novel tet peptidase complex from Pyrococcus horikoshii reveal an integrated peptide degradation system in hyperthermophilic archaea. Mol. Microbiol. 2009, 72, 26–40. [Google Scholar] [CrossRef] [PubMed]
  65. Remaut, H.; Bompard-Gilles, C.; Goffin, C.; Frere, J.M.; Van Beeumen, J. Structure of the Bacillus subtilis D-aminopeptidase Dppa reveals a novel self-compartmentalizing protease. Nat. Struct. Mol. Biol. 2001, 8, 674–678. [Google Scholar] [CrossRef]
  66. Schrödinger LLC. Schrödinger Release 2018-4: LigPrep; Schrödinger LLC: New York, NY, USA, 2018. [Google Scholar]
  67. Schrödinger Release 2018-4: Induced Fit; Docking Protocol; Induced Fit Docking protocol; Glide, Schrödinger, LLC: New York, NY, USA, 2016; Prime, Schrödinger, LLC: New York, NY, USA, 2018.
  68. Daina, A.; Michielin, O.; Zoete, V. SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci. Rep. 2017, 7, 42717. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Drinkwater, N.; Lee, J.; Yang, W.; Malcolm, T.R.; McGowan, S. M1 aminopeptidases as drug targets: Broad applications or therapeutic niche? FEBS J. 2017, 284, 1473–1488. [Google Scholar] [CrossRef]
  70. Peer, W.A. The role of multifunctional M1 metallopeptidases in cell cycle progression. Ann. Bot. 2011, 107, 1171–1181. [Google Scholar] [CrossRef] [Green Version]
  71. Tsujimoto, M.; Goto, Y.; Maruyama, M.; Hattori, A. Biochemical and enzymatic properties of the M1 family of aminopeptidases involved in the regulation of blood pressure. Heart Fail. Rev. 2008, 13, 285–291. [Google Scholar] [CrossRef]
  72. Mina-Osorio, P. The moonlighting enzyme CD13: Old and new functions to target. Trends Mol. Med. 2008, 14, 361–371. [Google Scholar] [CrossRef]
  73. Yeager, C.L.; Ashmun, R.A.; Williams, R.K.; Cardellichio, C.B.; Shapiro, L.H.; Look, A.T.; Holmes, K.V. Human aminopeptidase N is a receptor for human coronavirus 229E. Nature 1992, 357, 420–422. [Google Scholar] [CrossRef] [Green Version]
  74. Takahashi, H.; Hirose, K.; Watanabe, H. Necessity of meningococcal-glutamyl aminopeptidase for Neisseria meningitidis growth in rat cerebrospinal fluid (CSF) and CSF-like medium. J. Bacteriol. 2004, 186, 244–247. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Gonzalez-Bacerio, J.; Fando, R.; Monte-Martinez, A.; Charli, J.-L.; Chávez, M. Plasmodium falciparum M1-aminopeptidase: A promising target for the development of antimalarials. Curr. Drug Targets 2014, 15, 1144–1165. [Google Scholar] [CrossRef] [PubMed]
  76. Gardner, M.J.; Hall, N.; Fung, E.; White, O.; Berriman, M.; Hyman, R.W.; Carlton, J.M.; Pain, A.; Nelson, K.E.; Bowman, S.; et al. Genome sequence of the human malaria parasite Plasmodium falciparum. Nature 2002, 419, 498–511. [Google Scholar] [CrossRef] [PubMed]
  77. Su, L.; Jia, Y.; Zhang, L.; Xu, Y.; Fang, H.; Xu, W. Design, synthesis and biological evaluation of novel amino acid ureido derivatives as aminopeptidase N/CD13 inhibitors. Bioorg. Med. Chem. 2012, 20, 3807–3815. [Google Scholar] [CrossRef]
  78. Orning, L.; Krivi, G.; Fitzpatrick, F.A. Leukotriene A4 hydrolase. Inhibition by bestatin and intrinsic aminopeptidase activity establish its functional resemblance to metallohydrolase enzymes. J. Biol. Chem. 1991, 266, 1375–1378. [Google Scholar]
  79. Fournié-Zaluski, M.C.; Poras, H.; Roques, B.P.; Nakajima, Y.; Ito, K.; Yoshimoto, T. Structure of aminopeptidase N from Escherichia coli complexed with the transition-state analogue aminophosphinic inhibitor PL250. Acta Crystallogr. Sect. D Biol. Crystallogr. 2009, 65, 814–822. [Google Scholar] [CrossRef] [Green Version]
  80. Nikolaou, A.; Van Den Eynde, I.; Tourwé, D.; Vauquelin, G.; Tóth, G.; Mallareddy, J.R.; Poglitsch, M.; Van Ginderachter, J.A.; Vanderheyden, P.M.L. [3H]IVDE77, a novel radioligand with high affinity and selectivity for the insulin-regulated aminopeptidase. Eur. J. Pharmacol. 2013, 702, 93–102. [Google Scholar] [CrossRef]
  81. Papakyriakou, A.; Zervoudi, E.; Tsoukalidou, S.; Mauvais, F.-X.; Sfyroera, G.; Mastellos, D.C.; van Endert, P.; Theodorakis, E.A.; Vourloumis, D.; Stratikos, E. 3,4-Diaminobenzoic acid derivatives as inhibitors of the oxytocinase subfamily of M1 aminopeptidases with immune-regulating properties. J. Med. Chem. 2015, 58, 1524–1543. [Google Scholar] [CrossRef] [Green Version]
  82. Deprez-Poulain, R.; Flipo, M.; Piveteau, C.; Leroux, F.; Dassonneville, S.; Florent, I.; Maes, L.; Cos, P.; Deprez, B. Structure–activity relationships and blood distribution of antiplasmodial aminopeptidase-1 inhibitors. J. Med. Chem. 2012, 55, 10909–10917. [Google Scholar] [CrossRef]
  83. Wȩglarz-Tomczak, E.; Porȩba, M.; Byzia, A.; Berlicki, Ł.; Nocek, B.; Mulligan, R.; Joachimiak, A.; Dra̧g, M.; Mucha, A. An integrated approach to the ligand binding specificity of Neisseria meningitidis M1 alanine aminopeptidase by fluorogenic substrate profiling, inhibitory studies and molecular modeling. Biochimie 2013, 95, 419–428. [Google Scholar] [CrossRef] [Green Version]
  84. Mucha, A.; Lämmerhofer, M.; Lindner, W.; Pawełczak, M.; Kafarski, P. Individual stereoisomers of phosphinic dipeptide inhibitor of leucine aminopeptidase. Bioorg. Med. Chem. Lett. 2008, 18, 1550–1554. [Google Scholar] [CrossRef]
  85. Geng, N.; Zhang, W.; Li, Y.; Li, F. Aspartyl aminopeptidase suppresses proliferation, invasion, and stemness of breast cancer cells via targeting CD44. Anat. Rec. 2019, 302, 2178–2185. [Google Scholar] [CrossRef]
  86. Folkman, J. Angiogenesis in cancer, vascular, rheumatoid and other disease. Nat. Med. 1995, 1, 27–30. [Google Scholar] [CrossRef] [PubMed]
  87. Kelly, J.A.; Neidle, E.L.; Neidle, A. An aminopeptidase from mouse brain cytosol that cleaves N-terminal acidic amino acid residues. J. Neurochem. 1983, 40, 1727–1734. [Google Scholar] [CrossRef] [PubMed]
  88. Farrell, P.J.; Zopf, C.J.; Huang, H.J.; Balakrishna, D.; Holub, C.; Bilakovics, J.; Fanjul, A.; Matuszkiewicz, J.; Plonowski, A.; Rolzin, P.; et al. Using target engagement biomarkers to predict clinical efficacy of MetAP2 inhibitors. J. Pharmacol. Exp. Ther. 2019, 371, 299–308. [Google Scholar] [CrossRef]
  89. Haldar, M.K.; Scott, M.D.; Sule, N.; Srivastava, D.K.; Mallik, S. Synthesis of barbiturate-based methionine aminopeptidase-1 inhibitors. Bioorg. Med. Chem. Lett. 2008, 18, 2373–2376. [Google Scholar] [CrossRef] [Green Version]
  90. Sulpizio, A.C.; Pullen, M.A.; Edwards, R.M.; Louttit, J.B.; West, R.; Brooks, D.P. Mechanism of vasopeptidase inhibitor-induced plasma extravasation: Comparison of omapatrilat and the novel neutral endopeptidase 24.11/angiotensin-converting enzyme inhibitor GW796406. J. Pharmacol. Exp. Ther. 2005, 315, 1306–1313. [Google Scholar] [CrossRef] [PubMed]
  91. Prechel, M.M.; Orawski, A.T.; Maggiora, L.L.; Simmons, W.H. Effect of a new aminopeptidase P inhibitor, apstatin, on bradykinin degradation in the rat lung. J. Pharmacol. Exp. Ther. 1995, 275, 1136–1142. [Google Scholar]
  92. Gonzales, T.; Robert-Baudouy, J. Bacterial aminopeptidases: Properties and functions. FEMS Microbiol. Rev. 1996, 18, 319–344. [Google Scholar] [CrossRef]
  93. Stack, C.M.; Lowther, J.; Cunningham, E.; Donnelly, S.; Gardiner, D.L.; Trenholme, K.R.; Skinner-Adams, T.S.; Teuscher, F.; Grembecka, J.; Mucha, A.; et al. Characterization of the Plasmodium falciparum M17 leucyl aminopeptidase: A protease involved in amino acid regulation with potential for antimalarial drug development. J. Biol. Chem. 2007, 282, 2069–2080. [Google Scholar] [CrossRef] [Green Version]
  94. Eppinga, H.; Thio, H.B.; Schreurs, M.W.J.; Blakaj, B.; Tahitu, R.I.; Konstantinov, S.R.; Peppelenbosch, M.P.; Fuhler, G.M. Depletion of Saccharomyces cerevisiae in psoriasis patients, restored by Dimethylfumarate therapy (DMF). PLoS ONE 2017, 12, e0176955. [Google Scholar] [CrossRef] [PubMed]
  95. Marshall, B.J.; Armstrong, J.A.; McGechie, D.B.; Glancy, R.J. Attempt to fulfil Koch’s postulates for pyloric campylobacter. Med. J. Aust. 1985, 142, 436–439. [Google Scholar] [CrossRef] [PubMed]
  96. Dong, L.; Cheng, N.; Wang, M.W.; Zhang, J.; Shu, C.; Zhu, D.X. The leucyl aminopeptidase from Helicobacter pylori is an allosteric enzyme. Microbiology 2005, 151, 2017–2023. [Google Scholar] [CrossRef] [Green Version]
  97. Dutoit, R.; Brandt, N.; Legrain, C.; Bauvois, C. Functional characterization of two M42 aminopeptidases erroneously annotated as cellulases. PLoS ONE 2012, 7, e50639. [Google Scholar] [CrossRef] [Green Version]
  98. Leyhausen, G.; Schuster, D.K.; Vaith, P.; Zahn, R.K.; Umezawa, H.; Falke, D.; Müller, W.E.G. Identification and properties of the cell membrane bound leucine aminopeptidase interacting with the potential immunostimulant and chemotherapeutic agent bestatin. Biochem. Pharmacol. 1983, 32, 1051–1057. [Google Scholar] [CrossRef]
  99. Stamper, C.C.; Bienvenue, D.L.; Bennett, B.; Ringe, D.; Petsko, G.A.; Holz, R.C. Spectroscopic and X-ray crystallographic characterization of bestatin bound to the aminopeptidase from Aeromonas (Vibrio) proteolytica. Biochemistry 2004, 43, 9620–9628. [Google Scholar] [CrossRef]
Biomolecules 10 00659 sch001 550
Scheme 1. Retrosynthetic analysis of phosphinic dipeptide and dehydrodipeptide analogs. The P1 phenylethyl residue (homophenylalanine) is conserved as canonical. The pathways involve P1′-specific electrophiles that can be obtained from aldehydes.
Scheme 1. Retrosynthetic analysis of phosphinic dipeptide and dehydrodipeptide analogs. The P1 phenylethyl residue (homophenylalanine) is conserved as canonical. The pathways involve P1′-specific electrophiles that can be obtained from aldehydes.
Biomolecules 10 00659 sch001
Biomolecules 10 00659 g001 550
Figure 1. A general procedure for phosphinic ligand preparation. Following the R substituent selection, four possible stereoisomers (either R/S, R/S or R/S, E/Z combinations) were considered in the calculations.
Figure 1. A general procedure for phosphinic ligand preparation. Following the R substituent selection, four possible stereoisomers (either R/S, R/S or R/S, E/Z combinations) were considered in the calculations.
Biomolecules 10 00659 g001
Biomolecules 10 00659 g002 550
Figure 2. Generalized docking scheme. HVTS, high throughput virtual screening; SP, standard precision; IFD, induced fit docking; XP, extra precision.
Figure 2. Generalized docking scheme. HVTS, high throughput virtual screening; SP, standard precision; IFD, induced fit docking; XP, extra precision.
Biomolecules 10 00659 g002
Biomolecules 10 00659 g003 550
Figure 3. Modeled binding mode of the phosphinic dipeptide analog containing a 4-hydroxy-3-methoxyphenyl P1′ fragment in the active center of human aminopeptidase N (Protein Data Bank (PDB): 4FYT). Amino acid residues of the inhibitor and enzyme are shown as sticks, while zinc ions are shown as gray spheres. Hydrogen bonds and nonpolar interactions are shown as thin dashed lines.
Figure 3. Modeled binding mode of the phosphinic dipeptide analog containing a 4-hydroxy-3-methoxyphenyl P1′ fragment in the active center of human aminopeptidase N (Protein Data Bank (PDB): 4FYT). Amino acid residues of the inhibitor and enzyme are shown as sticks, while zinc ions are shown as gray spheres. Hydrogen bonds and nonpolar interactions are shown as thin dashed lines.
Biomolecules 10 00659 g003
Biomolecules 10 00659 g004 550
Figure 4. Modeled binding mode of a phosphinic dipeptide analog containing a 3-(furan-2-carbonyl)quinolin-2-yl P1′ fragment in the active center of cystinyl aminopeptidase (PDB: 5MJ6).
Figure 4. Modeled binding mode of a phosphinic dipeptide analog containing a 3-(furan-2-carbonyl)quinolin-2-yl P1′ fragment in the active center of cystinyl aminopeptidase (PDB: 5MJ6).
Biomolecules 10 00659 g004
Biomolecules 10 00659 g005 550
Figure 5. Modeled binding mode of a phosphinic dehydrodipeptide analog containing a phenanthren-4-yl P1′ fragment in the active center of ERAP2 aminopeptidase (PDB: 5AB0), with specified interactions as dashed lines: hydrogen bonds (green), π-alkyl (light pink), π–π stacking (dark pink).
Figure 5. Modeled binding mode of a phosphinic dehydrodipeptide analog containing a phenanthren-4-yl P1′ fragment in the active center of ERAP2 aminopeptidase (PDB: 5AB0), with specified interactions as dashed lines: hydrogen bonds (green), π-alkyl (light pink), π–π stacking (dark pink).
Biomolecules 10 00659 g005
Biomolecules 10 00659 g006 550
Figure 6. Modeled binding mode of a phosphinic dipeptide analog containing a 3-borono-5-isopropoxyphenyl P1′ fragment in the active center of aminopeptidase P3 (PDB: 5 × 49). Amino acid residues of the inhibitor and enzyme are shown as sticks, while cobalt ions are shown as purple spheres. Hydrogen bonds and nonpolar interactions are shown as thin dashed lines.
Figure 6. Modeled binding mode of a phosphinic dipeptide analog containing a 3-borono-5-isopropoxyphenyl P1′ fragment in the active center of aminopeptidase P3 (PDB: 5 × 49). Amino acid residues of the inhibitor and enzyme are shown as sticks, while cobalt ions are shown as purple spheres. Hydrogen bonds and nonpolar interactions are shown as thin dashed lines.
Biomolecules 10 00659 g006
Biomolecules 10 00659 g007 550
Figure 7. Modeled binding mode of a phosphinic dipeptide analog in the active center of aspartyl aminopeptidase (PDB: 4DYO) with the reverse orientation and a 5-(ethoxymethyl)furan-2-yl fragment (nominally P1′) extended to the S2 pocket.
Figure 7. Modeled binding mode of a phosphinic dipeptide analog in the active center of aspartyl aminopeptidase (PDB: 4DYO) with the reverse orientation and a 5-(ethoxymethyl)furan-2-yl fragment (nominally P1′) extended to the S2 pocket.
Biomolecules 10 00659 g007
Biomolecules 10 00659 g008 550
Figure 8. Modeled binding mode of a phosphinic dehydrodipeptide analog containing a 3-borono-6-ethoxyphenyl P1′ fragment in the active center of Helicobacter pylori aminopeptidase A/I (PDB: 4ZLA). The arrangement of the designed inhibitor was compared with that of bestatin, whose carbon atoms are colored in purple. Similar coloration was applied to amino acids interacting with the selected ligand. Carbon atoms of Tyr433 of another protein subunit are colored in dark blue. Amino acid residues of the inhibitor and enzyme are shown as sticks, while zinc ions are shown as gray spheres. Hydrogen bonds and nonpolar interactions are shown as thin dashed lines.
Figure 8. Modeled binding mode of a phosphinic dehydrodipeptide analog containing a 3-borono-6-ethoxyphenyl P1′ fragment in the active center of Helicobacter pylori aminopeptidase A/I (PDB: 4ZLA). The arrangement of the designed inhibitor was compared with that of bestatin, whose carbon atoms are colored in purple. Similar coloration was applied to amino acids interacting with the selected ligand. Carbon atoms of Tyr433 of another protein subunit are colored in dark blue. Amino acid residues of the inhibitor and enzyme are shown as sticks, while zinc ions are shown as gray spheres. Hydrogen bonds and nonpolar interactions are shown as thin dashed lines.
Biomolecules 10 00659 g008
Biomolecules 10 00659 g009 550
Figure 9. Modeled binding mode of a phosphinic dipeptide analog containing a 4-(difluoromethoxy)-3-hydroxyphenyl P1′ fragment in the active center of Streptococcus pneumoniae glutamyl aminopeptidase (PDB: 3KL9). Amino acid residues of the inhibitor and enzyme are shown as sticks, while zinc ions are shown as gray spheres. Hydrogen bonds and nonpolar interactions are shown as thin dashed lines.
Figure 9. Modeled binding mode of a phosphinic dipeptide analog containing a 4-(difluoromethoxy)-3-hydroxyphenyl P1′ fragment in the active center of Streptococcus pneumoniae glutamyl aminopeptidase (PDB: 3KL9). Amino acid residues of the inhibitor and enzyme are shown as sticks, while zinc ions are shown as gray spheres. Hydrogen bonds and nonpolar interactions are shown as thin dashed lines.
Biomolecules 10 00659 g009
Biomolecules 10 00659 g010 550
Figure 10. Modeled binding mode of a phosphinic dipeptide analog containing a 2-methylindol-7-yl P1′ fragment in the active center of Streptococcus pneumoniae glutamyl aminopeptidase (PDB: 3KL9). Amino acid residues of the inhibitor and enzyme are shown as sticks, while cobalt ions are shown as gray spheres. Hydrogen bonds and nonpolar interactions are shown as thin dashed lines. The surface is colored according to the hydrophobicity of amino acid side chains.
Figure 10. Modeled binding mode of a phosphinic dipeptide analog containing a 2-methylindol-7-yl P1′ fragment in the active center of Streptococcus pneumoniae glutamyl aminopeptidase (PDB: 3KL9). Amino acid residues of the inhibitor and enzyme are shown as sticks, while cobalt ions are shown as gray spheres. Hydrogen bonds and nonpolar interactions are shown as thin dashed lines. The surface is colored according to the hydrophobicity of amino acid side chains.
Biomolecules 10 00659 g010
Biomolecules 10 00659 g011 550
Figure 11. Correlation plot of inhibitory activity of phosphinic dipeptide and dehydrodipeptide analogs (pKi) reported for human aminopeptidase N [8,10] with the corresponding values of the calculated free Gibbs energy (R2 = 0.91).
Figure 11. Correlation plot of inhibitory activity of phosphinic dipeptide and dehydrodipeptide analogs (pKi) reported for human aminopeptidase N [8,10] with the corresponding values of the calculated free Gibbs energy (R2 = 0.91).
Biomolecules 10 00659 g011
Table
Table 1. Metallo-aminopeptidases included in the MEROPS database and selected for this study; the originating organisms and the crystal structures are referenced. PDB, Protein Data Bank.
Table 1. Metallo-aminopeptidases included in the MEROPS database and selected for this study; the originating organisms and the crystal structures are referenced. PDB, Protein Data Bank.
MEROPS IDNameOrganismPDB CodeMetal Ion(s)
M01.001Aminopeptidase NHomo sapiens4FYT [14]Zn
M01.001Aminopeptidase NSus scrofa4FKE [15]Zn
M01.003Aminopeptidase AHomo sapiens4KX7 [16]Zn
M01.004Leukotriene A4 hydrolaseHomo sapiens4MKT [17]Zn
M01.005Aminopeptidase NEscherichia coli3B34 [18]Zn
M01.011Cystinyl aminopeptidaseHomo sapiens5MJ6 [19]Zn
M01.021Tricorn interacting factor F3Thermoplasma acidophilum1Z5H [20]Zn
M01.024ERAP2 aminopeptidaseHomo sapiens5AB0 [21]Zn
M01.029M1 aminopeptidase EyPlasmodium falciparum3EBH [22]Zn
M01.031Cold-active aminopeptidaseColwellia psychrerythraea3CIA [23]Zn
M01.034Leukotriene A4 hydrolaseSaccharomyces cerevisiae2XQ0 [24]Zn
M01.unsAminopeptidase NNeisseria meningitidis2GTQ [25]Zn
M17.001Leucyl aminopeptidaseBos taurus1LAM [26]Zn, Zn
M17.002Leucyl aminopeptidaseSolanum lycopersicum4KSI [27]Mg, Mg
M17.003Aminopeptidase AEscherichia coli1GYT [28]Zn, Zn
M17.016Aminopeptidase A/IHelicobacter pylori4ZLA [29]Zn, Zn
M17.A05Leucyl aminopeptidaseCaenorhabditis elegans2HC9 [30]Zn, Zn
M17.unsLeucyl aminopeptidasePseudomonas putida3H8G [31]Mn, Zn
M17.unsLeucyl aminopeptidaseXanthomonas oryzae3JRU [32]Zn, Zn
M18.002Aspartyl aminopeptidaseBos taurus3VAT [33]Mg, Zn
M18.002Aspartyl aminopeptidaseHomo sapiens4DYO [34]Zn, Zn
M18.002Aspartyl aminopeptidasePseudomonas aeruginosa4NJR [35]Zn, Zn
M18.003Aspartyl aminopeptidasePlasmodium falciparum4EME [36]Zn, Zn
M18.004Aminopeptidase IClostridium acetobutylicum2GLJ [37]Mn, Mn
M18.004Aminopeptidase IThermotoga maritima2GLF [38]Mn, Mn
M24.001Methionyl aminopeptidase 1Escherichia coli2GGC [39]Co, Co
M24.001Methionyl aminopeptidase 1Mycobacterium tuberculosis1YJ3 [40]Co, Co
M24.001Methionyl aminopeptidase 1Rickettsia prowazekii3MX6 [41]Mn, Mn
M24.002Methionyl aminopeptidase 2Homo sapiens1B6A [42]Co, Co
M24.017Methionyl aminopeptidase 1Homo sapiens2GZ5 [43]Co, Co
M24.035Methionyl aminopeptidasePyrococcus furiosus1XGS [44]Co, Co
M24.036Methionyl aminopeptidase 1Staphylococcus aureus1QXY [45]Co, Co, Co
M24A.unMethionyl aminopeptidase 2Encephalitozoon cuniculi3FM3 [46]Fe, Fe
M24.004Aminopeptidase PEscherichia coli1WL9 [47]Mn, Mn
M24.009Aminopeptidase P1Homo sapiens3CTZ [48]Mn, Mn
M24.026Aminopeptidase P3Homo sapiens5 × 49 [49]Mn, Mn
M24.031Aminopeptidase PThermotoga maritima2ZSG [50]Zn
M24.038Aminopeptidase PPlasmodium falciparum5JR6 [51]Mn, Mn
M24.A11Aminopeptidase PXanthomonas campestris5CDE [52]Zn, Zn
M28.003Aminopeptidase SStreptomyces griseus1CP7 [53]Zn, Zn
M28.unsNon-peptidase homologuesClostridium acetobutylicum3K9T [54]Zn
M28.unsNon-peptidase homologuesParabacteroides distasonis3TC8 [55]Zn
M28.002Aminopeptidase Ap1Vibrio proteolyticus1RTQ [56]Zn, Zn
M28.020AM-1 aminopeptidaseAneurinibacillus sp. AM-12EK8 [57]Zn, Zn
M29.001Aminopeptidase TThermus thermophilus2AYI [58]Zn, Zn
M29.005Aminopeptidase SStaphylococcus aureus1ZJC [59]Co, Co
M42.001Glutamyl aminopeptidaseBacillus subtilis1VHE [60]Zn
M42.001Glutamyl aminopeptidaseStreptococcus pneumoniae3KL9 [61]Zn, Zn
M42.003PhTET1 aminopeptidasePyrococcus horikoshii2CF4 [62]Co, Co
M42.004PhTET2 aminopeptidasePyrococcus horikoshii1Y0Y [63]Zn, Zn
M42.009PhTET3 aminopeptidasePyrococcus horikoshii2VPU [64]Zn, Zn
M55.001D-aminopeptidase DppABacillus subtilis1HI9 [65]Zn, Zn
Table
Table 2. Selected M1 aminopeptidases with the P1′-preferred amino acid in the S1′ binding pocket with residues forming the cavity.
Table 2. Selected M1 aminopeptidases with the P1′-preferred amino acid in the S1′ binding pocket with residues forming the cavity.
Merops IDEnzyme NameP1′Amino Acids of S1′
M01.001Aminopeptidase NF, G, WGly352, Ala353, Arg381, Thr384, Val385, His388, Glu389, Ser415, Glu418, Tyr477
M01.003Aminopeptidase AE, M, RGly357, Ala358, Asn371, Gln385, Thr389, Val390, His393, Glu394, Ser420, Glu423, Phe424, Asp444, Tyr479, Ser480
M01.004Leukotriene A4 hydrolaseG, HGly268, Gly269, Asp291, Val292, His295, Glu296, Val322, Glu325, Tyr383, Arg563, Lys565
M01.005Aminopeptidase N
(bacterial)		G, NGly261, Ala262, Tyr275, Arg293, Val294, His297, Asp327, Tyr381, Glu382
M01.011Cystinyl aminopeptidaseR, T, VGly428, Ala429, Thr442, Leu457, Lys460, Ile461, His464, Glu465, Thr491, Glu494, Tyr549
M01.024ERAP2 aminopeptidaseT, Y, VGly334, Ala335, Ser348, Trp363, Arg366, Val376, His370, Glu371, Lys397, Glu400, Tyr455
Table
Table 3. Selected members of the M1 family with preferred P1′ substituents (R, compare Figure 1) that were the most favored, had the best interaction with the catalytic metal ion via the phosphinate group, and were the most specific (compare explanation in the text) are shown. The stereochemistry of the ligands and the calculated Gibbs free energy of binding to the active center are included. Known inhibitors are given as the references, with inhibition potency (Ki or IC50) and calculated ΔG values.
Table 3. Selected members of the M1 family with preferred P1′ substituents (R, compare Figure 1) that were the most favored, had the best interaction with the catalytic metal ion via the phosphinate group, and were the most specific (compare explanation in the text) are shown. The stereochemistry of the ligands and the calculated Gibbs free energy of binding to the active center are included. Known inhibitors are given as the references, with inhibition potency (Ki or IC50) and calculated ΔG values.
EnzymeInhibitorP1′ Substituent/Known InhibitorsConfigurationΔG [kJ/mol]
M01.001
Homo sapiens
4FYT		The best found4-Hydroxy-3-methoxyphenylS,S−53.82
Metal-interacting
The most specific5-Amino-thiophen-2-ylR,S−34.19
Ki = 0.69 nMR = 4-(Aminomethyl)phenyl [8] −64.11
M01.001
Sus scrofa
4FKE		The best found4-Boronophenyl MIDA esterR,R−39.09
Metal-interacting
The most specific
IC50 = 8.10 µMBestatin [77] −45.38
M01.003
Homo sapiens
4KX7		The best found3-Isobutylaziridin-2-yl R,S−38.80
Metal-interacting
The most specific1-Methyl-1H-imidazol-5-ylR,S−26.34
Ki = 75 µMBestatin [16] −7.54
M01.004
Homo sapiens
4MKT		The best found2-Borono-5-methoxyphenylR,S−57.39
Metal-interactingnone
The most specific2-Butyl-5-chloro-1H-imidazol-4-ylR,S−40.94
Ki = 200 nMBestatin [78] −53.73
M01.005
Escherichia coli
3B34		The best found3,4-DihydroxyphenylR,E−49.91
Metal-interacting
The most specificnone
Ki = 1.5 nMAla-P[(O)(OH)CH2]-Phe-Phe [79] −50.47
M01.011
Homo sapiens
5MJ6		The best found3-(Furan-2-carbonyl)quinolin-2-ylR,S−117.58
Metal-interacting
The most specificnone
Ki = 1.7 nM IVDE771 [80] −53.58
M01.021
Thermoplasma acidophilum
1Z5H		The best found4-(6-(Hydroxymethyl)pyridin-2-yl)phenylR,S−37.43
Metal-interacting6-Ethyl-4-oxo-4H-chromen-3-ylR,S−37.15
The most specific−37.15
UnknownBestatin −6.42
M01.024
Homo sapiens
5AB0		The best found2-BenzyloxyphenylR,S−76.49
Metal-interactingPhenanthren-4-ylR,Z−36.15
The most specificDiphenylmethylR,Z−52.41
IC50 = 237 nMDABA-Trp derivative2 [81] −54.74
M01.029
Plasmodium falciparum
3EBH		The best found3-Hydroxy-4-methoxyphenylR,S−63.40
Metal-interacting
The most specific4-(2-Pyridyl)phenylR,E−63.35
IC50 = 6 nMBDM144713 [82] −37.38
M01.031
Colwellia psychrerythraea
3CIA		The best found2-(tert-Butylthio)phenylR,R−76.28
Metal-interacting
The most specific
Unknown valueBestatin −66.93
M01.034
Saccharomyces cerevisiae
2XQ0		The best found3,5-Di-tert-butylphenylS,R−53.62
Metal-interacting
The most specific3-phenyl-1H-pyrazol-4-ylR,R−44.73
Unknown valueBestatin −47.06
M01.un
Neisseria meningitides
2GTQ		The best foundQuinolin-4-ylR,E−54.98
Metal-interacting
The most specificBiphenyl-2-ylR,E−45.78
Ki = 54 nMR = 4-Hydroxyphenyl [83] −26.61
1 IVDE77, peptidomimetic analogue β2hVal-Tyr-Ile-Aia-Nva-Phe; β2HVal, β2-homovaline; Aia, 4-amino-1,2,4,5-tetrahydro-indolo [2,3-c]azepin-3-on, a constrained Trp analogue. 2 (S)-2-(3-Amino-4-((S)-2-aminohexanamido)benzamido)-3-(1H-indol-3-yl)propanoic acid. 3 BDM14471, (2Z)-2-Benzylidene-N-(4-fluorobenzyl)-N′-hydroxymalonamide.
Table
Table 4. Homo sapiens M18 and M24 aminopeptidases with the P1′-preferred amino acid in the S1′ binding pocket with residues forming the cavity.
Table 4. Homo sapiens M18 and M24 aminopeptidases with the P1′-preferred amino acid in the S1′ binding pocket with residues forming the cavity.
Merops IDEnzyme NameP1′Amino Acids of S1′
M18.002Aspartyl aminopeptidaseA, H, LGly123, Glu301, Glu302, Val303, Gly304, Ser305, Gly414, Thr415, His440
M24.002Methionyl aminopeptidase 2A, S, VAla230, Asn327, Leu328, Asn329, Gly330, His331, His339, Thr343, Glu364, Phe366, Leu447
M24.017Methionyl aminopeptidase 1A, S, VCys211, Tyr300, Gly302, His303, His310, Val315, His317, Glu336, Met338, Gln365
M24.009Xaa-Pro aminopeptidase 1PHis485, His489, Glu523, Tyr527, Arg535, Glu537
M24.026Xaa-Pro aminopeptidase 3A, S, PAsp110, Trp163, Leu313, His420, His421, Val422, Gly423, His424, His431, Glu451, Pro452, Arg473
Table
Table 5. Selected human metallo-aminopeptidases with P1′ substituents that were the most favored, had the best interaction with the catalytic metal ion via the phosphinate group, and were the most specific are shown. The stereochemistry of the ligands and the calculated Gibbs free energy of binding in the active center are included. Known inhibitors are given as references.
Table 5. Selected human metallo-aminopeptidases with P1′ substituents that were the most favored, had the best interaction with the catalytic metal ion via the phosphinate group, and were the most specific are shown. The stereochemistry of the ligands and the calculated Gibbs free energy of binding in the active center are included. Known inhibitors are given as references.
EnzymeInhibitorP1′ Substituent/Known InhibitorsConfigurationΔG [kJ/mol]
M18.002
Homo sapiens
4DYO		The best found3,4-DibenzyloxyphenylR,S−88.05
The most specificNone
Metal-interacting5-(Ethoxymethyl)furan-2-ylS,S−46.49
C = 50 μg/mL – 50%Phenantroline [87] −16.74
M24.002
Homo sapiens
1B6A		The best found4-Borono-2,3-difluorophenylR,S−84.15
The most specific3-(3-Chlorophenyl)propylS,S−55.00
Metal-interactingPhenylpropylS,S−59.33
Ki = 26 nMBeloranib [88] −28.39
M24.017
Homo sapiens
2GZ5		The best found4-BoronophenylR,E−85.18
The most specific4-Borono-3-fluorophenylR,R−51.19
Metal-interacting4-(6-(Hydroxymethyl)pyridinyl)phenylS,R−77.96
Ki = 10 nMBarbiturate-based inhibitor1 [89] −53.46
M24.009
Homo sapiens
3CTZ		The best found4-Borono-2-fluorophenylR,R−51.91
The most specific2-BoronophenylS,S−41.61
Metal-interactingNone
Ki = 260 nMGW7964062 [90] −42.11
M24.026
Homo sapiens
5 × 49		The best found3-Phenylaziridin-2-ylR,S−102.22
The most specific2-Borono-6-fluorophenylR,R−68.93
Metal-interacting3-Borono-5-isopropoxyphenylR,R−83.35
Ki = 640 nMApstatin [91] −103.24
1 5-[(p-Dimethylaminostyryl)methylidene]barbituric acid. 2 GW796406, N-[(2S)-2-(mercaptomethyl)-3-methylbutanoyl]-4-(1H-pyrazol-1-yl)-l-phenylalanine.
Table
Table 6. Selected bacterial aminopeptidases with the P1′-preferred amino acid in the S1′ binding pocket with residues forming the cavity.
Table 6. Selected bacterial aminopeptidases with the P1′-preferred amino acid in the S1′ binding pocket with residues forming the cavity.
Merops IDEnzyme NameP1′Amino Acids of S1′
M17.016Aminopeptidase A/IND1Lys258, Lys270, Asn338, Asp340, Ala341, Arg344. Leu368, Thr369, Gly370, Ala371, Ile428, Gly435, Thr438, Ala439
M24.036Aminopeptidase SND1Ala75, Asp104, Leu165, Thr166, His168, His175, Ala179, His180, Glu202, Phe204, Phe221
M42.001Glutamyl aminopeptidaseND1Gly89, Gly90, Glu213, Glu214, Gly216, Leu217, Gly292, Thr293, His318
1. The enzyme specificities with kcat or KM were not determined (ND).
Table
Table 7. Selected bacterial metallo-aminopeptidases with P1′ substituents that were the most favored, had the best interaction with the catalytic metal ion via the phosphinate group, and were the most specific are shown. The stereochemistry of the ligands and the calculated Gibbs free energy of binding in the active center are included. Known inhibitors are given as references.
Table 7. Selected bacterial metallo-aminopeptidases with P1′ substituents that were the most favored, had the best interaction with the catalytic metal ion via the phosphinate group, and were the most specific are shown. The stereochemistry of the ligands and the calculated Gibbs free energy of binding in the active center are included. Known inhibitors are given as references.
EnzymeInhibitorP1′ Substituent/Known Inhibitors ConfigurationΔG [kJ/mol]
M17.016
Helicobacter pylori
4ZLA		The best found3-Borono-6-ethoxyphenylS,Z−81.94
Metal-interacting
The most specificNone
IC50 = 49 nMBestatin [96] −26.60
M24.036
Staphylococcus aureus
1QXY		The best foundBenzo [1,3]dioxol-5-ylR,R−93.88
Metal-interacting
The most specific 2-Methylindol-7-ylR,S−59.33
IC50 = 7 µM3-Amino-1-(cyclopropylamino)
heptan-2-one [45]		 −34.55
M42.001
Streptococcus pneumoniae
3KL9		The best found2-Chloro-5-phenylpyridin-3-ylR,Z−67.27
Metal-interacting
The most specific4-(Difluoromethoxy)-3-hydroxyphenylR,S−56.17
Ki = 292 nMBestatin1 [97] −41.97
1 Inhibition measured for the Thermotoga maritima tetrahedral-shaped dodecameric aminopeptidase of the M42 family.
Table
Table 8. Selected ADMET properties for M1 predicted inhibitors: MW, molecular weight; H-acc, H-bond acceptor number; H-don, H-bond donor number; LogP, logarithm of partition coefficient; GI, gastrointestinal absorption; Lipinski, fulfilling the rule of five; otherwise, the exceptions are given in parentheses.
Table 8. Selected ADMET properties for M1 predicted inhibitors: MW, molecular weight; H-acc, H-bond acceptor number; H-don, H-bond donor number; LogP, logarithm of partition coefficient; GI, gastrointestinal absorption; Lipinski, fulfilling the rule of five; otherwise, the exceptions are given in parentheses.
P1′ GroupMWH-accH-donLogPWater SolubilityGILipinski
4-Hydroxy-3-methoxyphenyl406.39620.5Moderately solubleLowYes
5-Amino-thiophen-2-yl381.41420.16SolubleLowYes
4-Boronophenyl MIDA ester515.2881−1.61Moderately solubleLow1 (MW)
3-Isobutylaziridin-2-yl368.4142−1.13SolubleHighYes
1-Methyl-1H-imidazol-5-yl365.3642−0.64SolubleHighYes
2-Borono-5-methoxyphenyl434.2173−1.07SolubleLowYes
2-Butyl-5-chloro-1H-imidazol-4-yl440.88521.19Poorly solubleLowYes
3,4-Dihydroxyphenyl390.35630.12SolubleLowYes
3-(Furan-2-carbonyl)quinolin-2-yl505.48711.62Poorly solubleLow1 (MW)
4-(6-(Hydroxymethyl)pyridin-2-yl)phenyl467.47621.02Poorly solubleLowYes
6-Ethyl-4-oxo-4H-chromen-3-yl456.45611.33Poorly solubleLowYes
2-Benzyloxyphenyl466.49512.1Poorly solubleHighYes
Diphenylmethyl448.47412.34Poorly solubleHighYes
Phenanthren-4-yl458.47412.67Poorly solubleHighYes
3-Hydroxy-4-methoxyphenyl406.39620.38Moderately solubleLowYes
4-(2-Pyridyl)benzaldehyde435.43511.38Poorly solubleHighYes
2-(tert-Butylthio)phenyl448.54412.09Poorly solubleLowYes
3,5-Di-tert-butylphenyl472.58413.28Poorly solubleHighYes
3-Phenyl-1H-pyrazol-4-yl426.43520.77Poorly solubleLowYes
Quinolin-4-yl409.39511.04Moderately solubleHighYes
Biphenyl-2-yl434.44412.27Poorly solubleHighYes
Table
Table 9. Correlation plot of inhibitory activity of bestatin (pKi) reported for different metallo-aminopeptidases with the corresponding values of the calculated free Gibbs energy (R2 = 0.65).
Table 9. Correlation plot of inhibitory activity of bestatin (pKi) reported for different metallo-aminopeptidases with the corresponding values of the calculated free Gibbs energy (R2 = 0.65).
EnzymeKi [M]ΔG [kJ/mol] Biomolecules 10 00659 i001
M01.0013.03 · 10−6 [98]−36.74
M01.0014.05 · 10−6 [77]−45.38
M01.0037.50 · 10−6 [78]−7.54
M01.0042.00 · 10−7 [79]−53.73
M01.0294.81 · 10−7 [22]−55.83
M17.0019.00 · 10−9 [98]−59.77
M28.0021.80 · 10−8 [99]−71.01
M42.0012.92 · 10−7 [97]−41.97

Share and Cite

MDPI and ACS Style

Talma, M.; Mucha, A. P1′ Residue-Oriented Virtual Screening for Potent and Selective Phosphinic (Dehydro) Dipeptide Inhibitors of Metallo-Aminopeptidases. Biomolecules 2020, 10, 659. https://doi.org/10.3390/biom10040659

AMA Style

Talma M, Mucha A. P1′ Residue-Oriented Virtual Screening for Potent and Selective Phosphinic (Dehydro) Dipeptide Inhibitors of Metallo-Aminopeptidases. Biomolecules. 2020; 10(4):659. https://doi.org/10.3390/biom10040659

Chicago/Turabian Style

Talma, Michał, and Artur Mucha. 2020. "P1′ Residue-Oriented Virtual Screening for Potent and Selective Phosphinic (Dehydro) Dipeptide Inhibitors of Metallo-Aminopeptidases" Biomolecules 10, no. 4: 659. https://doi.org/10.3390/biom10040659

APA Style

Talma, M., & Mucha, A. (2020). P1′ Residue-Oriented Virtual Screening for Potent and Selective Phosphinic (Dehydro) Dipeptide Inhibitors of Metallo-Aminopeptidases. Biomolecules, 10(4), 659. https://doi.org/10.3390/biom10040659

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop