Introduction
Natural products and compounds directly derived from them continue to play a significant role in modern drug discovery and development (Newman and Cragg, 2007). The most effective compounds targeting chitinases, the enzymes that catalyze the hydrolysis of β(1,4)-linked N-acetylglucosamine (chitin), are natural products (reviewed in Andersen et al., 2005). Chitin is an important structural component of the cell wall of fungi (Munro and Gow, 2001), the microfilarial sheath in parasitic nematodes (Harris et al., 2000) as well as the cuticles of insects (Merzendorfer, 2006). Although chitin is absent from mammalian physiology, fungi, nematodes, and insects must all metabolize chitin at important stages of their life cycle; hence, small-molecule chitinase inhibitors have the potential to act as biocontrol agents in both medicine and agriculture. Several natural product chitinase inhibitors have indeed been shown to exhibit interesting biological activities. For example, the pseudotrisaccharide allosamidin, the most extensively studied and generally the most potent chitinase inhibitor, has been shown to inhibit cell separation in fungi (Kuranda and Robbins, 1991; Sakuda et al., 1990), block malaria parasite penetration into the mosquito midgut (Tsai et al., 2001), and inhibit insect larvae development (Sakuda, 1996). Additionally, this compound was shown to decrease the number of inflammatory cells in a mouse model of asthma by targeting the acidic mammalian chitinase (AMCase) (Zhu et al., 2004). Although total syntheses of allosamidin have been reported (Berecibar et al., 1999), the length and complexity of the synthetic routes involved severely limit both the availability of this compound and the scope for structure-based design of novel allosamidin-derived inhibitors.
Recently, there has been considerable interest in the ability of peptides to mimic carbohydrate-protein interactions and bind to enzymes and lectins (e.g., Johnson and Pinto, 2002; Vyas et al., 2003; Dharmasena et al., 2007). Synthetic compounds that can mimic such peptide ligands, i.e., peptidomimetics, could provide potential drug candidates with improved stability toward proteolytic breakdown but retaining high potency. Two classes of peptide-based chitinase inhibitors are known: the cyclic proline-containing dipeptides isolated from broth of a marine bacterium (Izumida et al., 1996a) and the cyclic pentapeptides argadin and argifin secreted from fungal strains (Omura et al., 2000; Arai et al., 2000). The synthesis of the natural product cyclo-(L-Arg-D-Pro), or CI-4, is straightforward (Izumida et al., 1996b; Houston et al., 2002a) and in principle further analogs are readily available by substitution of either amino acid residue. However, although several cyclic dipeptide derivatives have been synthesized, their chitinase inhibition is moderate (high micromolar range) (Houston et al., 2002a, 2004). The total syntheses of both argadin and argifin have been achieved (Dixon et al., 2005, 2006) using a combination of solid phase and solution chemistry. For argifin (1, Figure 1A), IC50 values of 27 nM (Dixon et al., 2005), 33 μM (Houston et al., 2002b), and 4.5 μM (Rao et al., 2005c) have been determined for chitinase B1 from Aspergillus fumigatus (AfChiB1), chitinase B from Serratia marcescens, and human chitotriosidase (hCHT), respectively. However, although argifin is a potent chitinase inhibitor, its molecular weight and number of hydrogen-bond donors/acceptors well exceed the Lipinski criteria for drug-likeliness (Lipinski, 2000), while its peptide-based framework is also likely to be prone to proteolytic degradation.
Here, we use the available AfChiB1-argifin complex crystal structure to design linearized fragments of this natural product cyclopentapeptide. Using a combination of X-ray crystallographic analysis of the peptide-protein complexes and enzymology, we show that these argifin fragments maintain their binding mode to the enzyme and continue to inhibit potently. Surprisingly, the (tiny) dimethylguanylurea fragment appears to be responsible for forming most of the key interactions and represents an unusually efficient binder and an attractive starting point for structure-based optimization.
Results and Discussion
Design and Synthesis of Linear Peptide Inhibitors
The natural product argifin 1 (Figure 1A) is a potent inhibitor of several family 18 chitinases, acting through extensive (peptide) mimicry of the natural carbohydrate substrate (Figure 2; Houston et al., 2002b; Rao et al., 2005b). The structure of the compound in complex with chitinase B from A. fumigatus (Rao et al., 2005b; Figure 2) was used as a basis to design a range of peptide fragments that all contain the dimethylguanylurea group previously shown to deeply penetrate the active site pocket (Houston et al., 2002b; Rao et al., 2005b; Figure 2). Hence, the linear tetra-, tri-, and dipeptides 2–4 were generated (Figure 1), along with an isolated Arg derivative 5 with an N-methylcarbamoyl-derivatized side chain, and the dimethylguanylurea 6 itself (Figure 1).
For the preparation of 2–4, Arg-containing di-, tri-, and tetrapeptide precursors were assembled by manual Fmoc SPPS on 2-chlorotrityl chloride polystyrene resin. After cleavage from the resin and side chain deprotection, the Arg(MC) modification was introduced by acylation with MeNHCO2Su/DBU/DMF at 40°C, as described for the synthesis of argifin. 2–4 were isolated by preparative HPLC and were fully characterized by NMR and high resolution mass spectrometry. 5 was prepared from Fmoc-L-Arg(Pmc)-OH by standard solution chemistry, but using methyl isocyanate to install the Arg(MC) modification, and was isolated by preparative HPLC and characterized as previously. 6 was prepared from methylguanidine hydrochloride by neutralization with ion exchange resin (Dowex hydroxide form), followed by treatment with methyl isocyanate and isolation as for 5. Full details of the syntheses are available in the Supplemental Data available with this article online.
Investigation of Inhibition and Binding Mode of the Argifin Peptide Fragments
Argifin has been shown to be a competitive inhibitor of family 18 chitinases, binding with affinities from in the nanomolar to micromolar range (Arai et al., 2000; Houston et al., 2002b; Rao et al., 2005b). To compare the linear peptides synthesized in this study, inhibition was tested against a fungal chitinase (AfChiB1), the human chitotriosidase (hCHT), and total chitinase activity from a lung homogenate (LH) from a mouse model of chronic asthma. This lung homogenate has previously been shown to mainly contain the acidic mammalian chitinase (AMCase) that has recently been proposed as an antiasthma drug target (McMillan and Lloyd, 2004; Zhu et al., 2004; Schüttelkopf et al., 2006; Table 1). In parallel, AfChiB1 crystals were soaked in solutions containing compounds 2–6 and the resulting crystal structures investigated. The five structures were solved to a maximum resolution of 1.90–2.2 Å (Table 2), similar to the published (2.0 Å) AfChiB1-argifin complex, and thus the structures possess a comparable amount of solvent detail. The data yielded electron density maps showing clear positions for the ligand, confirming the structure and stereochemistry of the synthesized linear peptides (Figure 2).
Argifin itself binds to family 18 chitinases through subsites (i.e., where the individual sugars of the chitin polymer are known to bind) −1, +1, and +2 of the active site, following standard glycoside hydrolase subsite nomenclature (Davies et al., 1997), with hydrolysis taking place between the −1 and +1 sugar. In the previously published AfChiB1-argifin crystal structure (Figure 2), all inhibitor peptide bonds are in the trans configuration except for the cis L-Arg to N-MePhe bond (Rao et al., 2005b). The N-methyl carbamoyl nitrogen hydrogen bonds to the side chain of Asp175, whereas its oxygen forms a hydrogen bond to the hydroxyl group of Tyr245. The arginine guanidinium group forms two hydrogen bonds with Glu177 (the catalytic acid) and one with the side chain of Asp246. Additionally, the first isoAsp side chain generates a hydrogen bond with the Trp137 side chain, and the argifin arginine carbonyl forms two hydrogen bonds with the Arg301 side chain. One intramolecular hydrogen bond is observed between the first isoAsp backbone N and the D-Ala carbonyl O, associated with a β-turn centered on the L-Arg-L-MePhe dipeptide moiety (Figure 2). Numerous water-mediated hydrogen bonds are observed involving Arg301, Glu322, and Asn323. The side chains of Tyr48, Met243, and Trp384 generate hydrophobic interactions and surround the N-methyl carbamoyl group, whereas the side chains of Trp137 and Phe251 form extensive stacking interactions with the argifin L-Phe. Almost all of these residues are highly conserved in family 18 chitinases, explaining the potency of this inhibitor reported for several enzymes.
Argifin Linearization and Truncation Generates an Active Tetrapeptide
Linearization of the cyclopentapeptide argifin 1 by removal of one of the Asp residues results in tetrapeptide 2 and surprisingly leads to retention of significant inhibitory activity (IC50 = 4.3 μM against AfChiB, Table 1). In agreement with this, the crystal structure (Figure 2) reveals that it binds the chitinase with the Arg(MC) side chain in an identical position to that observed in the AfChiB1-argifin complex, with a similar conformation of the peptide backbone (shifts up to 0.6 Å for equivalent argifin backbone atoms). Hence, the inhibitor fills the −1, +1, and +2 subsites, and all direct hydrogen bonds and intramolecular hydrogen bonds observed in the argifin complex are conserved (Figure 2). The only significant differences are the loss of water-mediated hydrogen bonds to Glu322/Asn323 and a cis amide bond in the N-terminal N-acetyl group, allowing an additional internal hydrogen bond between the C-terminal carboxylate and the backbone nitrogen of D-Ala (Figure 2). The total surface area buried by argifin and 2 are identical (141 Å2). Thus, it seems likely that removal of the second isoAsp and linearization of the natural product lowers the affinity of binding to the chitinases through entropic penalties incurred upon complexation of the more flexible, acyclic structure of 2 and/or loss of several water-mediated hydrogen bonds.
Further Truncation of the Tetrapeptide Generates an Active Tripeptide
It is apparent from the AfChiB-argifin/tetrapeptide complexes that the D-Ala side chain does not make significant interactions with the enzyme active site (Figure 2). Indeed, the tripeptide 3, lacking the D-Ala residue, inhibits chitinase activity of AfChiB1, hCHT, and LH with IC50 values similar to those of the tetrapeptide 2 (Figure 1B; Table 1). Again, crystallographic analysis of a complex with AfChiB (Figure 2) shows that the Arg(MC) side chain binds in its usual position and that the tripeptide backbone adopts a conformation almost identical to that observed in the AfChiB1-argifin complex (shifts up to 0.5 Å for equivalent argifin backbone atoms). The differences from the tetrapeptide are the displacement of the isoAsp C terminus, the loss of the extra internal hydrogen bond, and an extra water-mediated hydrogen bond with Ser250.
An Argifin Dipeptide Fragment Showing Micromolar Inhibition
To further narrow down the pharmacophore of the argifin natural product, the tripeptide-protein interactions of 3 were evaluated for the individual amino acids to design the most potent dipeptide. Both the derivatized argifin and the N-methylated phenylalanine make tight interactions with the active site (Figure 2) and truncation of these amino acids was expected to significantly compromise binding. The only interaction the remaining isoAsp side chain makes with the protein is to accept a hydrogen bond from the Trp137 indole, although this interaction is fully conserved in the argifin, tetrapeptide and tripeptide complexes, and mutation of the tryptophan is known to completely abolish argifin inhibition (Rao et al., 2005b; Figure 2). Nevertheless, it was decided to synthesize dipeptide 4, which inhibits the chitinases with IC50 values less than an order of magnitude higher than those measured for the tri/tetrapeptides 2 and 3 (Figure 1B; Table 1). The AfChiB1-dipeptide complex (Figure 2) reveals that the terminal Arg(MC) moiety generates the same five hydrogen bonds with Asp175, Glu177, Tyr245, and Asp246 as seen in the AfChiB1 complexes with argifin. However, the loss of an internal hydrogen bond and the hydrogen bond with the side chain of Trp137 appears to induce significant conformational changes in the backbone (shifts up to 3.4 Å for equivalent argifin backbone atoms), with the L-Arg-N-methyl-L-Phe peptide bond adopting a trans configuration (Figure 2). Consequently, the Arg301 to L-Arg carbonyl hydrogen bond is replaced with a hydrogen bond from Arg301 to the MePhe carboxylate, and the MePhe side chain is displaced to only partially occupy the +2 subsite through a reduced interaction with Phe251 (Figure 2). The shift of the inhibitor backbone also causes Trp137 to settle in a dual conformation, with an additional conformation pointing toward the indole group of Trp384, which itself is also displaced. However, these conformational changes also result in the generation of three new water-mediated hydrogen bonds that may partially compensate for the loss of direct hydrogen bonds. In addition, the considerably smaller dipeptide generates a protein contact surface of equivalent size to that observed in the AfChiB1-tripeptide complex (Table 2).
The N-Methylcarbamoylated Arginine Is the Pharmacophoric Argifin Amino Acid
Despite the linearization and truncation of argifin down to a dipeptide, significant activity is retained and the pharmacophore is thus likely to be an even smaller fragment. We next investigated the activity of the acetylated single amino acid, “monopeptide” 5. Strikingly, this monopeptide still inhibits AfChiB1 with an IC50 value of 81 μM (Figure 2; Table 1), although this is 3 orders of magnitude higher than inhibition by the intact cyclopentapeptide 1 and 40-fold higher than the linear tetrapeptide 2. The crystal structure of the AfChiB-monopeptide complex shows that despite the removal of the MePhe residue and concomitant loss of the hydrophobic interactions with Trp137/Phe 251, the modified arginine side chain still establishes the same five hydrogen bonds as observed in the argifin complex (Figure 2). Similar to the dipeptide complex, Trp137 is flipped with a 120° rotation around the Cα-Cβ bond, allowing stacking interactions with the L-Arg backbone and the N-terminal acetyl group.
Dimethylguanylurea Is an Active Chitinase Inhibitor Fragment
Because the terminal atoms of the modified arginine residue of 1–5 appear to make the majority of the interactions with the protein, we decided to investigate inhibition of the dimethylguanylurea 6 itself. Strikingly, this small fragment, consisting of only 9 heavy atoms and with a molecular weight of 130 Da, still shows detectable inhibition of AfChiB1 (IC50 = 500 μM) and other chitinases (Figure 1B; Table 1). In the dimethyl-guanylurea complex crystal structure with AfChiB1, the fragment is well ordered, occupies the same position as observed in the argifin and linear peptide complexes, and is able to make five hydrogen bonds (Figure 2). The reduction in potency compared to the monopeptide can be explained by the loss of the L-Arg carbonyl to Arg301 hydrogen bond as well as a considerably smaller buried surface area (Table 2) due to loss of interactions with Trp137 (Figure 2).
Concluding Remarks
In order to establish the minimum elements of the argifin structure that are required for activity and provide the basis for the design of new inhibitors derived from this lead, five acyclic fragments of the natural product have been synthesized, inspired by the previously published structures of argifin-chitinase complexes (Houston et al., 2002b; Rao et al., 2005b). The mode of binding of these fragments with AfChiB has been elucidated and allowed differences in inhibitor potency observed upon incremental structural simplification to be explained. Although such analyses of potency changes upon systematic deconstruction of optimized synthetic inhibitors have previously been reported (Hajduk, 2006), to the best of our knowledge this is the first example of a combined kinetic and crystallographic investigation of the binding of a systematically dissected natural product.
The argifin peptide fragments can be ranked tetrapeptide > tripeptide > dipeptide > monopeptide > dimethylguanylurea with respect to their potencies (Table 1). This rank can be explained by the progressive loss of protein-inhibitor contacts observed in the crystal structures of the complexes in conjunction with predicted entropic penalties linked to linearization of a conformationally restricted cyclic peptide. However, all peptide fragments, including the minimal dimethylguanylurea, maintain a core minimum of interactions, namely the hydrogen bonds to Asp175, Glu177 (the catalytic acid), Tyr245, and Asp246 and stacking with Trp384. With the exception of Asp246, these residues are conserved in family 18 chitinases throughout all kingdoms of life. The ranking of the argifin fragments is also in good agreement with a previously published mutagenesis study to identify key AfChiB-argifin interactions, which showed large effects on argifin affinity if amino acids either hydrogen bonding to the modified argifin side chain or showing stacking interactions were mutated (Rao et al., 2005b). Other known family 18 chitinase inhibitors, such as allosamidin, diketopiperazines, and purine derivatives, also derive their potencies from interacting with these active site residues (Terwisscha van Scheltinga et al., 1994; Houston et al., 2002a; Rao et al., 2005a).
The surprising result of this study is that the pharmacophore of the large cyclopentapeptide argifin 1 is a 9 atom fragment, the dimethylguanylurea 6. In absolute terms, this is a poor inhibitor, with an IC50 of 500 μM compared to 27 nM for the natural product (Table 1). However, viewing the same data in terms of ligand efficiency (Kuntz et al., 1999; Hopkins et al., 2004; Abad-Zapatero and Metz, 2005), here defined as binding efficiency index (BEI) (Abad-Zapatero and Metz, 2005), dimethylguanylurea is the most efficient argifin fragment with a BEI of 27.7, a significantly more efficient binder than argifin itself (BEI = 11.7 with a molecular weight of 676 Da). As evidenced by the structures of argifin-chitinase complexes (Houston et al., 2002b; Rao et al., 2005b; Figure 2), there are relatively large complementary buried surfaces on the inhibitor and protein, suggesting that while not very efficient, argifin may be a very selective inhibitor. Conversely, it is likely that the small dimethylguanylurea fragment with its limited contacts will also fit a number of other active sites.
In a given series, inhibitor potency commonly correlates with molecular weight, and hence ligand efficiency is a more useful measure than potency when selecting a lead compound for further optimization into a drug or monitoring progress in a drug development program. Recently, a retrospective analysis of several drug discovery campaigns suggested that selecting relatively small, efficient binders as a starting framework increases the chances of obtaining a final optimized ligand with high potency that does not violate the Lipinski rule-of-five (Hajduk, 2006; Lipinski, 2000) and fits the restrictions for acceptable absorption and permeability properties. The dimethylguanylurea fits all Lipinski rules, with 4 hydrogen bond donors, 1 hydrogen bond acceptor, a molecular mass of 130 Da, and clogP = −0.48. In addition, removal of all peptide bonds eliminates argifin's intrinsic susceptibility to proteolytic breakdown and is thus likely to give more favorable pharmacokinetic properties. Using the observation that, in an ideal optimization program, ligand efficiency remains constant throughout lead optimization (Hopkins et al., 2004; Hajduk, 2006), dimethylguanylurea serves as an attractive fragment lead with a speculative potential to be optimized to a compound with a Ki of 60 nM at a molecule weight of 260 Da or a Ki of 16 pM at a weight of 390 Da.