Shane M.
Devine
,
Christopher A.
MacRaild
,
Raymond S.
Norton
* and
Peter J.
Scammells
*
Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC 3052, Australia. E-mail: ray.norton@monash.edu; peter.scammells@monash.edu
First published on 4th November 2016
Malaria continues to frustrate humanity's attempts to eradicate this deadly disease. Although gains have been made over the last 15 years, drug resistance to malaria continues to be a major concern. The lack of new antimalarials with novel mechanisms of action continues to challenge the scientific community to find innovative targets to combat this persistent disease. One such target, apical membrane antigen 1 (AMA1), is an essential protein that helps the parasite invade host erythrocytes. Recently, a number of efforts have focused on the druggability of this target, aiming to block the interactions of AMA1 that mediate invasion of host cells. This review covers recent progress in drug discovery targeting this crucial protein–protein interaction in malaria.
000 deaths.1 Most concerning is the fact that 306
000 of these deaths were children under five years of age.
Artemisinin-based therapies have played a central role in recent successes against malaria,2 although resistance to these front-line treatments is now a serious threat, with detection of resistance in five countries in the Greater Mekong region of Asia, and the likelihood of further spread.1 This problem, coupled with the fact that the pipeline of novel antimalarial agents under development suffers from a lack of diversity, underscores the need for new targets to combat malaria.3 Despite recent progress towards a malaria vaccine,4,5 the most advanced of the current candidates possess only modest and short-lived efficacy, and there is no established timetable for widespread clinical roll-out.6 Therefore, there is a clear and present unmet need for novel targets and new antimalarials.
The AMA1–RON complex co-localises with a tight junction that forms between the apical tip of the parasite and its host as the parasite prepares to invade (Fig. 1A).10–12 This junction moves from the apex of the parasite to its posterior end as invasion proceeds, and so is known as the moving junction. By this mechanism, the parasite enters the host cell, forming the parasitophorous vacuole in which it will reside. Once the moving junction reaches the posterior end of the parasite, a membrane fusion event closes the parasitophorous vacuole, separating it from the host cell plasma membrane.
![]() | ||
| Fig. 1 (A) The role of AMA1 in host-cell invasion. Figure modified from Richard et al.21 (B) the conserved hydrophobic cleft (green) and DII loop (blue) on PfAMA1 (PDB: 1Z40).22 (C) The structure of AMA1 bound to RON2 residues 2021–2059 (purple), showing the displacement of the DII loop (top left; PDB: 3ZWZ).23 | ||
Although the molecular details of moving junction formation and the precise role of AMA1 in this process remain unclear, an extensive body of evidence suggests that AMA1 plays an essential role in invasion. Repeated attempts to generate AMA1 knock-outs in several Plasmodium species have failed to yield stable clonal lines, suggesting the gene is essential to parasite growth in blood-stage culture.13–15 Despite clear evidence of strong diversifying pressure acting on AMA1 from the host immune system,16,17 the AMA1–RON2 interaction remains highly conserved across the apicomplexan phylum.18 Furthermore, antibodies targeting AMA1, arising through natural infection or from experimental immunisations, block the AMA1–RON2 interaction, inhibit blood-stage culture of P. falciparum parasites, and are correlated with protection from disease in vivo.19,20
There has been some recent controversy about the essentiality of AMA1 and its role in the moving junction. On the basis of a series of transient AMA1 knock-down and knockout lines in Toxoplasma and P. bergei, Ménard and colleagues posited that AMA1 acts independently of the RON complex and is dispensable in any step of the invasion process.13,15,24 However, several factors complicate the analysis of these lines. In Toxoplasma, paralogues of both AMA1 and RON2 exist, which appear able to complement the loss of AMA1.25 Moreover, no stable clonal knockout lines were established in P. bergei, and the principal assay of host cell invasion used in this work may be unable to distinguish productive invasion from aborted invasion events in which the parasite remains associated with the host cell but is unable to undergo intracellular development.26 Irrespective of the dispensability of AMA1 at the genetic level, a wide range of antibody and peptide inhibitors of the AMA1–RON2 interaction inhibit host cell invasion and parasite growth in culture.27 Accordingly, the potential efficacy of AMA1 as a target for chemical inhibition is now accepted.13
Nonetheless, significant challenges must be overcome if anti-malarials targeting AMA1 are to be developed. The large and relatively flat interface between AMA1 and RON2 is typical of a protein–protein interaction (PPI), and presents a number of challenges shared by this type of target.28,29 Moreover, the formation of the moving junction and the AMA1–RON2 interaction itself is a transient event, occurring only briefly during each ∼48 hour asexual cycle. Together with the fact that AMA1 appears to be present in excess on the parasite surface,26 this may indicate that a saturating amount of inhibitor needs be present at the time of invasion in order to achieve a therapeutic outcome. AMA1 is, however, directly exposed to the host bloodstream prior to invasion, thus removing the difficulties associated with accessing intracellular targets and avoiding the potential for resistance by drug-efflux mechanisms.
![]() | ||
| Fig. 2 Small molecule AMA1 inhibitors (A) scaffolds identified via FBDD (1–4);43–46 (B) pyrrolopyrimidines (5–7) identified by Srinivasan et al.;47 (C) and diazepinone 8 identified by Pihan et al.48 | ||
Initially, a series of 4-aryl substituted 2-aminothiazoles was synthesised.43 However, when these were analysed by NMR and SPR, the structure–activity relationships (SAR) obtained were flat and inconsistent. We then reviewed our fragment screening campaigns at MIPS and found that in 14 screens against our fragment library at least one 2-aminothiazole hit all 14 targets by STD NMR.42 This led us to investigate all 2-aminothiazoles, of which there were 28 in our fragment library, for their potential promiscuity by SPR, NMR and computational methods. We also interrogated two high throughput-screening (HTS) libraries, one from academia and one from industry. Although the HTS libraries showed some signs of frequent hitter behaviour, the result was much more pronounced in the fragment set. The 28 fragments in our library were then examined against six unrelated protein targets and, although some patterns emerged, such as the association of the free 2-amino group with increased promiscuity, no clear mechanism of action was evident. Nonetheless, the SAR obtained was confusing and flat against several of the targets, as previously seen in the characterization of pan-assay interference compounds (PAINS) by HTS.49 We therefore concluded that this class of molecule should be excluded from fragment libraries, unless several orthogonal approaches are utilised to avoid these so called promiscuous 2-aminothiazoles (PrATs).45 The remaining scaffolds identified in the original fragment screen are currently under investigation and the outcome of these studies will be reported in due course.
![]() | ||
| Fig. 3 Peptide ligands R1 and RON2 make similar interactions with AMA1. R1 (orange) (PDB: 3SRJ)23 adopts a similar conformation to residues 2027–2047 of RON2 (blue) (PDB: 3ZWZ).23 Five residues in each peptide are identical (red), and make similar interactions with AMA1. Figure modified from Wang et al.35 | ||
A recent peptide engineering exercise identified two point mutations to RON2 that afford a 1–2 order of magnitude increase in binding affinity for AMA1. On the background of a 29-residue peptide based on residues 2027–2055 of RON2, these mutants give rise to KD values of 2.1 and 0.4 nM, respectively, for 3D7 and FVO AMA1.50 This is the most potent strain-transcending peptide identified to date. Nonetheless, in this peptide, as in the parental RON2 and R1, critical interactions are broadly distributed over the surface of AMA1 in ways that are unlikely to be effectively mimicked in a drug-like peptidomimetic inhibitor.
000 member library, in which a truncated RON2 peptide was conjugated to the AlphaScreen donor bead, and AMA1 to the AlphaScreen acceptor.47 This screen identified 20 initial hits, 14 of which were tested in blood-stage parasite culture, resulting in three hits that blocked merozoite invasion in vitro with IC50 values in the range 21–29 μM. Re-synthesis of pyrrolopyrimidine 5 (Fig. 2B) and chemical elaboration of this scaffold, produced two molecules (6 and 7) that showed enhanced inhibition, with reported IC50 values of 9.8 μM and 6 μM, respectively. The authors used immunoprecipitation assays to show that these compounds prevented the association of RON2 with AMA1 and electron microscopy of invading parasites to infer that the compounds blocked moving junction formation. Nonetheless, they were unable to quantify a direct interaction between these molecules and AMA1, impeded by poor compound solubility. Moreover, SPR data showed clear evidence of weak, super-stoichiometric binding. In an attempt to address these issues, we synthesized the three active compounds and examined them using our previously-described biophysical techniques based on NMR and SPR.43 Compounds 5–7 showed clear evidence of colloidal aggregation over the concentration range employed for in vitro testing, suggesting these assays may be confounded by the well-known non-specific effects of colloidal aggregates in a broad range of biological assays.52,53 We therefore attempted to improve some of the problematic physicochemical properties of this scaffold, including the high calculated partition co-efficients (clog
P),54 a measure of solubility, and the lipophilic ligand efficiency (LLEAT),55 an indicative measure of potency relative to molecular size and lipophilicity. To address these issues, we replaced the 7-cyclopentyl group with a methyl group to reduce the lipophilicity and explored substitutions on the 5-aryl ring, whilst maintaining the 4-amino group of the pyrrolo[2,3-d]pyrimidine scaffold. Although the Srinivasan compounds47 and our novel, more soluble, analogues56 showed modest activity as anti-malarials in a growth inhibition assay (GIA), their affinity for AMA1 (where measurable) was low and did not correlate with GIA activity. Likewise, a third group has been unable to quantify binding of 5 to AMA1, either by SPR or in a fluorescence anisotropy based R1 competition assay (described below).48 We therefore concluded that these compounds were acting in an AMA1-independent manner. As this scaffold is prevalent in the kinase literature, we believe another pathway, possibly the calcium-dependent protein kinase 1 (CDPK), may be the major target in the activity of these compounds.56
557
052 structures were selected from the initial 203
738
892 structures interrogated. Utilising the crystal structure of the PfAMA1–PfRON2sp1 (PDB: 3ZWZ)23 and the PLANTS program,58in silico docking experiments were conducted. The binding site was a 12 Å sphere centred on the oxygen of the side-chain of Tyr251, a key residue near the centre of the hydrophobic cleft.59 Ten poses were generated per structure and the best-scoring pose was selected for subsequent analysis. The binding poses of these structures were analysed for their interaction with three PfAMA1 residues, Val169, Phe183 and Tyr234. A final visualisation step reduced the structure-based analysis of 7161 down to eight compounds. However, it is unclear how this number was so significantly culled. These compounds included 5 diazepinones, an acetamide and two known drugs: Goserelin, for the treatment of breast and prostate cancer, and Indinavir, a retroviral treatment for HIV/AIDS. They also included the MMV Malaria Box in their computational workflow, resulting in 1 hit, MMV019881.
In order to test these compounds, a fluorescence anisotropy-binding assay was developed. This assay used a fluorescein-labelled peptide probe that was outcompeted by a putative small-molecule inhibitor binding to AMA1. A relatively weak probe is required as the nanomolar binder F*PfRON2sp1 was too potent to permit observation of competitive displacement. A truncated version of the potent binder R1 was utilised featuring a 5-carboxyfluoroscein (5-FAM) linked by an aminohexanoic acid (Ahx) to a 13-residue version with sequence: 5-FAM-Ahx-ALPLFSKFGSRMH, known as F*R1-13A1. The nine hits selected from the in silico screening, together with the Srinivasan compound 5,47 were tested for their ability to disrupt the interaction between PfAMA1 and F*R1-13A1. Diazepinone 8 and pyrrolopyrimidine 5 showed moderate competitive behaviour, with an IC50 determined for diazepinone 8 of 24 μM (Fig. 2C). Limited solubility of pyrrolopyrimidine 5 prevented the quantification of its potency in this assay. The compounds were also tested for their ability to bind AMA1 by SPR, with none showing saturable interaction with the target. Nonetheless, the fluorescence anisotropy assay developed in this study could find greater utility in the future to interrogate putative inhibitory effects at AMA1.
This computational approach followed work that used in silico strategies to identify putative peptidomimetics of AMA1-binding peptides identified by random phage display.60,61 These molecules were then docked into PfAMA1 (PDB: 1Z40)22 along the hydrophobic cleft.44 Then, using the ZINC database,62 small molecules that were similar to the five top scoring results were blindly docked into the surface of PfAMA1. These compounds were not predicted to interact with the polymorphic regions of PfAMA1 and their putative binding site in the hydrophobic cleft may suggest their potential to inhibit the AMA1/RON2 interaction. Nonetheless, experimental validation of this work is lacking, and would be required before the utility of these compounds became evident.
The PfAMA1 sequence is relatively polymorphic, with variation being driven by pressure from the host immune system. These polymorphisms are clustered on a single face of PfAMA1. As a result, the potent R1 peptide and inhibitory antibodies are strongly strain-specific; for example R1 binds 3D7 AMA1 with a KD of 100 nM, but to the divergent FVO AMA1 with a KD of ≥500 μM.64 Since our goal was to develop a strain-transcending inhibitor of AMA1, we sought to determine if the sequence divergence between these two strains results in structural change affecting binding. Accordingly the X-ray crystal structure of FVO AMA1 was obtained and compared with previous 3D7 structures.64 Overall the structures were highly similar and suggested that the strain-specificity of R1 and other inhibitors reflects the loss of a few key interactions rather than any significant structural change. Unlike two published structures of 3D7 PfAMA1 (PDB: 1Z40)22 and 1F9–3D7 PfAMA1 complex (PDB: 2Q8A),33 where the DII loop is partially ordered, the DII loop in FVO PfAMA1 is completely disordered in this structure.64 This difference appears to be a consequence of the distinct crystallisation conditions or crystal contacts, rather than any intrinsic difference in the flexibility of the DII loop. Indeed, the flexibility of the DII loop is associated with a conformational change that allows interaction between AMA1 and its endogenous partner RON2,23 and is therefore likely to be essential for AMA1 function.65
19F NMR spectroscopy has proven to be a valuable method for probing ligand-induced conformational changes and has found utility with a number of biologically relevant protein targets.66,67 Our group has utilised 19F NMR spectroscopy to detect ligand-induced movement of the DII loop of AMA1.44 We achieved this by mutating Phe367, near the centre of the DII loop to Trp (Fig. 4A) and replacing the introduced and native Trp residues with 5-fluorotryptophan (5-F-Trp) by expression in the presence of glyphosphate. The resulting 5-F-Trp F367W AMA1 was correctly folded and bound RON2L and R1 with affinities that compared well with those of wild-type AMA1. Additional Trp-to-Phe mutants were used to assign the four native Trp residues in AMA1, each of which lies outside the DII loop (Fig. 4A). The 19F NMR spectrum of the mutant F367W reveals five signals, corresponding to the four native Trp residues in AMA1, and the Trp introduced at residue 367 (Fig. 4B). The linewidth of the signal from Trp367 is more than twice that of the native AMA1 Trp signals, indicating the presence of conformational exchange centred on the DII loop, consistent with evidence from NMR46 and X-ray crystallography64 that the DII loop is quite flexible.
![]() | ||
| Fig. 4 19F NMR as a probe of ligand binding to AMA1. (A) The structure of AMA1 (domains I and II; PDB: 2Z8V)32 showing the RON2 binding site (green), the DII loop (blue), the site of the introduced 5-F-Trp367 (red) and the four native Trp residues (yellow); (B) 19F NMR spectrum of 5-F-Trp labelled F367W AMA1 (blue) contains signals from the four Trp residues in wt AMA1, plus an additional broad signal attributed to the introduced 5-F-Trp367. Deconvolution of the spectrum as the sum of five Lorentzian signals is shown (grey). Addition of the peptide ligands R1 (green) or RON2L (purple) causes the 5-F-Trp367 resonance to become sharp. A series of aminothiazole fragments at 0 (blue), 1 (purple) and 3 mM (red) bind to AMA1 and cause concentration-dependent sharpening of the 5-F-Trp367 resonance. Samples contained ∼100 μM 5-F-Trp F367W AMA1 in 20 mM sodium phosphate, pH 7.4, and spectra were recorded at 25 °C and a 19F frequency of 564 MHz without proton decoupling. Figure modified from Norton et al.67 | ||
In the presence of the peptide ligands RON2L and R1, a significant change in the 19F NMR of 5-F-Trp-AMA1[F367W] was observed, with the 5-F-Trp367 undergoing a 0.2 ppm change in chemical shift and a 2.5-fold decrease in linewidth, to become the sharpest peak in the spectrum (Fig. 4B). We interpret this as a significant increase in the flexibility of W367 relative to the rest of AMA1, consistent with the DII loop being displaced from the hydrophobic cleft and becoming highly flexible, as inferred from the crystallographic results. We then tested the effect of members of our fragment series on the 5-F-Trp-AMA1[F367W] spectrum.44 Of three series tested, one caused concentration-dependent sharpening of the 5-F-Trp 367 resonance, similar to that caused by the peptide ligands (Fig. 4D). This result suggests that members of this series of elaborated fragments displace the DII loop in much the same way as the peptide ligands do. In contrast, members of two other series caused no change to the 5-F-Trp-AMA1[F367W] spectrum, despite competing with R1 and therefore binding in the hydrophobic cleft. This suggests that these molecules bind at a distinct sub-site within the cleft, distal from the DII loop, and accordingly do not displace the loop. Thus, this approach has the potential to distinguish amongst the several potential binding sites along the length of the cleft, and to reveal the structural mode of action of AMA1–RON2 inhibitors.
To further our understanding of the effects of polymorphism in PfAMA1, the backbone NMR resonance assignments of two P. falciparum strains 3D7 and FVO were obtained by NMR spectroscopy.46 The spectra of FVO PfAMA1 were of higher quality than 3D7, reflecting differences in sample stability and optimal solution conditions for the two proteins. The spectral differences may also reflect differences in conformational dynamics, with 3D7 AMA1 experiencing more extensive line-broadening due to intermediate timescale conformational exchange. Approximately 84% of expected backbone resonances were assigned in FVO AMA1, compared with the less than 50% assignment coverage possible for 3D7. These assignments were used to map the binding site of several small molecules including thiazoles, pyrazoles and benzimidazoles. Of these, three benzimidazoles gave small chemical shift perturbations in 2D [1H–15N]-TROSY experiments and showed binding to a region near the hydrophobic cleft at the tetrapeptide region Lys177, Asp178, Gly179 and Gly180. This information will help guide future development of high-affinity ligands as PfAMA1 inhibitors. Our current research is utilising these findings to further elaborate the 2-phenylbenzimidazole fragment hits.
Footnote |
| † The authors declare no competing interests. |
| This journal is © The Royal Society of Chemistry 2017 |