Yuri
Frosi‡
a,
Simon
Ng‡
a,
Yen-Chu
Lin‡
b,
Shimin
Jiang
a,
Siti Radhiah
Ramlan
a,
Dilraj
Lama
c,
Chandra S.
Verma
d,
Ignacio
Asial
*ef and
Christopher J.
Brown
*a
aDisease Intervention Technology Lab (DITL), IMCB (A*STAR), 8A Biomedical Grove, #06-04/05, Neuros/Immunos, 138648, Singapore. E-mail: cjbrown@imcb.a-star.edu.sg
bInsilico Medicine Taiwan Ltd., Suite 2013, No. 333, Sec.1, Keelung Rd., Xinyi Dist. 110, Taipei, Taiwan
cDepartment of Microbiology, Tumor and Cell Biology, Karolinska Institutet, Biomedicum Quarter 7B-C Solnavägen 9, 17165 Solna, Sweden
dBioinformatics Institute (A*STAR), 30 Biopolis Street, #07-01 Matrix, 138671, Singapore
eDotBio, 1 Research Link, 117604, Singapore. E-mail: ignacio.asial@dotbiopharma.com
fNanyang Technological University, School of Biological Sciences, Singapore
First published on 19th May 2022
Identifying new binding sites and poses that modify biological function are an important step towards drug discovery. We have identified a novel disulphide constrained peptide that interacts with the cap-binding site of eIF4E, an attractive therapeutic target that is commonly overexpressed in many cancers and plays a significant role in initiating a cancer specific protein synthesis program though binding the 5′cap (7′methyl-guanoisine) moiety found on mammalian mRNAs. The use of disulphide constrained peptides to explore intracellular biological targets is limited by their lack of cell permeability and the instability of the disulphide bond in the reducing environment of the cell, loss of which results in abrogation of binding. To overcome these challenges, the cap-binding site interaction motif was placed in a hypervariable loop on an VH domain, and then selections performed to select a molecule that could recapitulate the interaction of the peptide with the target of interest in a process termed Peptide Epitope Linker Evolution (PELE). A novel VH domain was identified that interacted with the eIF4E cap binding site with a nanomolar affinity and that could be intracellularly expressed in mammalian cells. Additionally, it was demonstrated to specifically modulate eIF4E function by decreasing cap-dependent translation and cyclin D1 expression, common effects of eIF4F complex disruption.
The translation of peptide compounds into suitable molecules for lead development and biological evaluation face several hurdles: (1) limited intracellular permeability and (2) proteolytic instability.1,5 The process to surmount this challenge is highly resource- and time-consuming. Therefore, to understand the biological and phenotypical consequences of the binding modes of these molecules further in a reasonable timeframe, a methodology is required to present these binding epitopes in a relevant manner to enable these studies. One attractive approach is to embed these binding epitopes within the context of a small and highly stable protein backbone (termed scaffold), where they will be doubly constrained at their N and C terminals. These modalities are usually referred to as peptide aptamers (PA).7 An advantageous property of PAs are that they are genetically encodable and can be expressed within mammalian cell lines and animal models to facilitate advanced studies.7,8 PAs are usually designed via 2 distinctive approaches either by insertion of a single sequence into a loop region (“loop on frame”) or via mutations of specific residues embedded in rigid secondary structural (“rigid motif”) elements within the selected scaffold protein.7,8
These approaches in addition to rational protein engineering techniques are both amenable to combinatorial display methods such as yeast and phage display. However, they have limitations: (1) specific or hypervariable loop sequences are inserted into scaffolds that are not optimized to stabilise their bound conformations, whilst (2) the use of “rigid motifs” leads to sampling limited areas of 3-dimensional space e.g., mutations that lie down one side of an α-helix. To circumvent these issues a 2-step discovery process (termed Peptide Epitope Linker Evolution, PELE) was utilized (Fig. 1A), whereby the surface of the target protein is mapped for novel peptides binders utilizing constrained and linear peptide phage libraries, where upon any isolated interaction motifs are inserted into a larger hypervariable loop located on a selected scaffold protein. Libraries are then constructed with the interaction motif (‘peptide epitope’) located at different positions within the hypervariable loop, and then the selection against the target protein re-performed to select for sequences (‘linker evolution’) that optimally present the interaction motif within the context of the scaffold. The advantage of this approach over generating naïve libraries is that the diversity of the libraries in both steps are now only focussed on finding the optimal solution for one variable each, ‘binding’ and ‘presentation’, respectively, rather than simultaneously in a single step. Thus, allowing for more effective searching of sequence space especially with techniques that attain lower thresholds of diversity (e.g., yeast ∼106–7vs. mRNA display ∼1012–13).
Fig. 1 (A) The schematic outlines the process termed Peptide Epitope Linker Evolution (PELE). Peptide phage display libraries are used to probe the surface of the target protein to discover new binding motifs and modalities (e.g., linear or constrained libraries). Upon identification, the novel motif or modality can be inserted into a larger hypervariable loop located on a selected scaffold protein. Several distinctive libraries can then be constructed with the interaction motif (‘peptide epitope’) located at different positions within the hypervariable loop, and then the selection against the target protein re-performed to select for sequences (‘linker evolution’) that optimally present the interaction motif within the context of the scaffold. (B) A brief outline of next generation sequencing (NGS) enhanced phage display. A selected phage library is panned against an immobilised target protein using three technical repeats and in parallel against corresponding negative control selections where the target protein is either removed or replaced with a different target protein. Bound phage is then eluted, amplified and sequenced using NGS protocols (NextSEQ, Illumina). The FASTQ file generated from the sequencing data was processed by in-house PYTHON scripts that identified the barcodes and constant flanking regions and extracted the reads of the correct length corresponding to the variable peptide library. The table presents the list of sequences identified from each selection with their associated abundance. The abundance is calculated by taking the copy number of each sequence and normalizing it by dividing the copy number by the total number of reads in each sequence. Sequences not observed in a specific replicate were assigned a copy number of zero. The enrichment ratio of each sequence in the target selection was calculated by determining the mean fraction from the target screen replicates and dividing it by the mean fraction from the selected control screen replicates. Since the denominator must not be a zero when taking the ratio, sequences found with a zero copy number are assigned an arbitrary copy number before performing the normalization. Significance of the ratio was assessed using one-tailed, unequal variance Welch test. A heat map (C and D) is then generated to identify the enriched peptides that have ratio and p-values corresponding to the parameters stated in the figure. Each individual colour coded block on the map represents the abundance of the unique sequence in each selection and the sequences are ordered by their ratio value. (C) Heatmap showing sequences enriched from the M13 disulphide constrained 7 mer library (C7C) against eIF4E, but not in the 2 control selections (Mdm2 and eIF4A). (D) Heatmap showing sequences enriched from the M13 linear 12 mer library against eIF4E, but not in the 3 control selections (Mdm2, K-RAS and DO-1). The sequence motif in (C) was generated from the enriched sequences using MEME,43 whilst in (D) it was generated from the sequences exhibiting the known eIF4E binding motif (YXXXL, X is any amino acid). |
The PELE process was used to identify novel cyclic peptides that interact with the eukarytic initiation factor 4E (eIF4E) at the capped-mRNA binding site and to evolve the identified interaction motif onto an autonomous, disulphide-free VH domain.9 eIF4E itself plays a critical role in cap-dependent translation and is an attractive target for anti-oncogenic therapy as it is frequently overexpressed and/or mis-regulated in many tumours, is associated with poor prognosis in patients, is implicated in chemoresistance10–12 and also orchestrates a cancer specific mRNA translational program e.g., c-myc, cyclin-D1.13,14 eIF4E binds the cap structure found at the 5’ end of mammalian mRNA as part of the eIF4F complex, whereupon the complex interacts with the initiation factor eIF3 via the eIF4G subunit, which enables the recruitment of the 43S pre-initiation complex (PIC). The PIC complex then transverses the 5′ UTR of the mRNA until it encounters the AUG start codon, where it initiates protein synthesis though formation of the 80s ribosomal complex.15 The formation of the eIF4F complex in normal cells is finely controlled by the 4E-BP proteins (4EBP1, 4EBP2 and 4EBP3), where they competitively inhibit the interaction of eIF4E with eIF4G. Both the 4E-BPs and eIF4G1 share a common eIF4E interaction motif (YXXXXLΦ, X = any amino acid and Φ = hydrophobic amino acid), which interacts at the same site on eIF4E.16 Disruption of eIF4G binding impairs the assembly of the eIF4F complex on the 5′ cap structure and prevents cap-dependent translation. The actions of 4EBPs are controlled by their phosphorylation status, which is under the specific control of mTOR.17 Aberrant eIF4E activity in cancer can be targeted using several approaches. One method is to disrupt eIF4F formation using inhibitors of the PI3K/AKT/mTORC pathway,18 which induce de-phosphorylation of the 4EBPs.15 Another strategy is to prevent eIF4E phosphorylation through the use of MNK1/2 kinase antagonists in the RAS/ERK signalling pathway.15 Other approaches that have been used to antagonise the activity of the eIF4F complex are ASOs (antisense oligonucleotides) to reduce eIF4E expression, inhibitors of the eIF4A helicase,15 and development of cap analogues and synthetic cap mimetics to prevent binding of capped mRNA.15,19,20 This includes the development of pro-nucleotide derivatives to overcome the cell permeability issues associated with cap-analogues.20–22 Additionally, small molecules and peptides can be used to target the interface between eIF4G and eIF4E.15,23,24
In this manuscript, we report the first peptide inhibitor of the eIF4E cap binding site and have used this to develop a novel peptide aptamer suitable for intracellular expression and therapeutic modelling studies. The structures of both the peptide and the peptide aptamer (PA) with eIF4E revealed a novel interaction with the cap-binding site of eIF4E that provides a unique template for future small molecule and peptidomimetic design strategies. Additionally, the VH domain-based PA was also used to construct a unique split luciferase with the potential to screen eIF4E cap antagonists directly in live cells and support future cap-analogue mimetic discovery.
In the absence of any sequence similarity with described eIF4E interacting proteins (eIF4G1 and the 4E-BP family) and the eIF4E interaction motif (enriched in the 12mer selection), competitive fluorescence anisotropy experiments were performed to delineate the binding sites of the EE-01 to EE-09 peptides on eIF4E utilising either a FAM labelled m7GTP (m7GTPFAM) or a FAM labelled eIF4G/4E-BP1 site interacting peptide (eIF4GFAM) (Fig. 2A–C). None of the disulphide constrained peptides were observed to displace eIF4GFAM from eIF4E, however several cyclic peptides did compete for binding at the cap-site with m7GTPFAM. EE-02 was determined to interact with eIF4E with a Kd of 406.2 ± 3.6 nM, 12-fold more strongly than the next best performing cyclic peptide (EE-09, Kd = 4,860 ± 633.4 nM). Alanine scanning mutagenesis experiments confirmed that the residues conserved in the disulphide constrained peptide motif (C2E3[M/L/T]4G5[F/Y]6[F/Y]7X8X9C, Fig. 1A) were all necessary for binding (Table 1). Mutation of Q8A in the EE-02 sequence had no effect on binding to eIF4E, whilst the replacement of D9A resulted in a small attenuation of the Kd (Table 1). Additionally, the dependence of EE-02 binding with eIF4E upon the disulphide constraint was demonstrated to be critical under reducing conditions.
Fig. 2 (A) A surface representation of eIF4E depicting the location of the m7GTP (capped mRNA) and eIF4G1/4E-BP1 binding sites. Locations of tryptophan residues whose intrinsic fluorescence is sensitive to binding by either m7GTP or peptides that interact with the eIF4G binding site are shown in green. (B) Competitive fluorescence anisotropy experiments with FAM labelled m7GTP assessing binding of the cyclic peptides to the cap-binding site. (C) Competitive fluorescence anisotropy experiments with FAM labelled eIF4G1 derived peptide assessing binding of the cyclic peptides to the eIF4G/4E-BP1 binding site. Apparent Kds (see Table 1) were determined by curve-fitting using Prism (Graphpad, Ltd). See materials and methods. (D) eIF4E intrinsic tryptophan fluorescence was assessed in response to titrations of m7GTP, PHAGESOL (Ac-KKRYSR*QLL*-NH2, * = 2-aminoisobutyric acid) and EE-02, respectively. |
Sequence | Competitive anisotropy (kd, nM) | ||
---|---|---|---|
m7GTPFAM | eIF4GFAM | ||
m 7 GDP | — | 200.7 ± 25.2 | No Binding |
m 7 GTP | — | 66.6 ± 7.2 | No Binding |
m 7 GTP BIOTIN | — | 34.7 ± 4.0 | No Binding |
EE-01 | NH2-AETGFFTGG-NH2 | 15030.0 ± 3078.0 | No Binding |
EE-02 | NH2-AEMGFFQDG-NH2 | 406.2 ± 3.6 | No Binding |
EE-03 | NH2-AELGYYNDG-NH2 | No binding | No Binding |
EE-04 | NH2-AETGFFLKG-NH2 | No binding | No Binding |
EE-05 | NH2-AELGFYRLG-NH2 | No binding | No Binding |
EE-06 | NH2-AETGFFLRG-NH2 | No binding | No Binding |
EE-07 | NH2-AETGYFSQG-NH2 | No binding | No Binding |
EE-08 | NH2-AETGFYKTG-NH2 | 16710 ± 3027.0 | No Binding |
EE-09 | NH2-AEMGYFGNG-NH2 | 4860 ± 633.4 | No Binding |
PHAGESOL | Ac-KKRYSR*QLL*-NH2 | No binding | 34.05 ± 6.42 |
EE-02 E3A | NH2-AAMGFFQDG-NH2 | No binding | ND |
EE-02 M4A | NH2-AEAGFFQDG-NH2 | No binding | ND |
EE-02 G5A | NH2-AEMAFFQDG-NH2 | No binding | ND |
EE-02 F6A | NH2-AEMGAFQDG-NH2 | No binding | ND |
EE-02 F7A | NH2-AEMGFAQDG-NH2 | No binding | ND |
EE-02 Q8A | NH2-AEMGFFADG-NH2 | 384.9 ± 62.1 | ND |
EE-02 D9A | NH2-AEMGFFQAG-NH2 | 1682 ± 537.4 | ND |
EE-02 N-Del | NH2-EMGFFQDG-NH2 | 66.5 ± 19.5 | ND |
EE-02 C-Del | NH2-AEMGFFQD–NH2 | 684.8 ± 54.8 | ND |
EE-02 NC-Del | NH2-EMGFFQD–NH2 | 525.9 ± 128.8 | ND |
Tryptophan quenching experiments were performed to compare the binding mechanisms of EE-02 with m7GTP and PHAGESOL (a phage modified eIF4E interacting peptide that binds at the eIF4G binding site)26,27 against eIF4E. It is well established that the binding of m7GTP to eIF4E results in significant tryptophan fluorescent quenching of W56 and W102 (Fig. 2A), both involved in recognising the m7G moiety of m7GTP (Fig. 2D), whilst peptides that interact with eIF4E at the eIF4G site also result in quenching of W73 (Fig. 2A and D).28 However, the reduction in total eIF4E fluorescence caused by peptide binding is significantly less compared to m7GTP binding at the cap-site (Fig. 2D). In contrast, the quenching of intrinsic tryptophan fluorescence by EE-02 produces a significantly different profile to both m7GTP and PHAGESOL (Fig. 2D), indicating that EE-02 binds via a mechanism substantially different to both.
Fig. 3 (A) The 2Fo–Fc electron density map (1.2σ) showing the EE-02 disulfide constrained peptide bound to eIF4E at the cap-binding site. EE-02 is highlighted in the blue mesh, structured waters in the red mesh and eIF4E interacting residues with green mesh. (B) Complex of eIF4E bound to m7GTP (PDB ID: 2V8W) indicating conformational differences with the EE-02:eIF4E complex structure. (C) EE-02 when bound to eIF4E forms a β-hairpin turn-like structure that is stabilized by intra hydrogen bonds between the backbone carbonyl of E3 and backbone amide of F6, and the backbone amide of E-03 and carbonyl of F6 (3.1 Å and 4.1 Å, respectively). The conformation of the cyclic peptide is further rigidified by hydrogen bonds between the C10 amide and the carbonyl of F7, and the backbone of D9 and the carbonyl of F7 (3.0 Å and 3.8 Å, respectively). The polypeptide backbone of EE-02 also forms a set of critical interactions with eIF4E (<3.2 Å) shown in yellow. (D–G) show the interactions made by the conserved residues of the cyclic peptide interaction motif in EE-02 (white) with eIF4E (green). (D) E3 electrostatically interacts with R112 and forms a water mediated hydrogen bond interaction with N155. (E) M4 forms a dipole interaction with the hydroxyl group of S92 and a variety of hydrophobic contacts with residues F48, W46, L60 and P100. (F) F6 forms hydrophobic contacts with the residues T203 A204, H200, W166 and W102 of eIF4E. (G) F7 forms stacking interactions with W56 and edge on face interactions with F48. Additionally, it forms a hydrophobic contact with P100. (H) The carbonyl group of M4 forms no direct interactions with eIF4E but forms a hydrogen bond with a structured water, which is part of a larger network of structured waters that facilitates the interaction of EE-02 with eIF4E. (I) Overlay of the EE-02:eIF4E complex with unbound eIF4E (PDB ID: 4BEA) and m7GTP bound eIF4E demonstrating the similarity of the EE-02 bound conformation to the apo structure. Ligands interacting at the cap binding site (EE-02 and m7GTP) are not shown for clarity. |
The EE-02 peptide forms a β-hairpin turn-like structure in the binding pocket that allows the side chains of the constrained peptide motif to efficiently interact with eIF4E (Fig. 3C). The glycine at position 5 due to its steric permissiveness enables optimal formation of a β-turn type structure, and in turn a stabilising intramolecular h-bond between the backbone carbonyl of E3 and backbone amide of F6. The E3 of the selected motif (E3MGFF7) forms direct electrostatic interactions with R112 of eIF4E (Fig. 3D), whilst M4 forms hydrophobic contacts with the back of the binding pocket and a specific hydrogen bond with S92 via its sulphur atom (Fig. 3E). Residue F6 forms a range of hydrophobic interactions (3.6 Å–4.2 Å) that include T203 A204, H200, W166 and W102 of eIF4E (Fig. 3F). In contrast, F7 forms stacking interactions with W56 and edge on face interactions with F48 (Fig. 3G). Additional main-chain interactions are also formed by EE-02 that contribute to the energetics of binding with eIF4E: the backbone carbonyl of C2 interacts with the R157 sidechain (Fig. 3C), the amide back bone of F7 form 2 water mediated interactions with the carboxylic group of E103 (Fig. 3C), and the backbone carbonyl of M4 interacts with a structure water network that involves h-bonds with N155 and R112 (Fig. 3H). Binding energy decomposition analysis from MD simulations of the eIF4E:EE-02 complex structure demonstrates that M4, F6 and F7 contribute a significant proportion of the binding energy (Fig. S1A, ESI†). The part of the cap-binding that recognizes the triphosphate tail (R157, K159 and K162) of m7GTP is not involved in binding EE-02, and undergoes negligible structural changes. Interestingly, the conformation of the eIF4E cap-binding site when bound to EE-02 is very similar to its unbound configuration,23,31,32 where the W102 side chain and the W56 loop also rotate and swing out of the cap-binding site (Fig. 3I). MD simulations of the unbound EE-02 were performed, indicating that its structure does not vary dramatically from the bound form and overall retains a similar fold to that observed in the crystal form (Fig. S2A and B, ESI†).
Fig. 4 (A) The CDR3 region (highlighted in red) of the VH-DIF scaffold (PDB ID: 7D8B) was selected for replacement by rationally designed loops. The engineered loops were designed to present the EE-02 motif in the correct conformation to interact with eIF4E using polyglycine linkers. However, the VH-DIF derived proteins, when tested, exhibited no binding to eIF4E. Peptide Epitope Loop Exchange (PELE) libraries were also constructed and inserted at the same site in the VH-DIF scaffold. Optimal linkers needed to present the EE-02 motif correctly for binding were selected by YSD (Yeast Surface Display). The YSD (yeast surface display) selection against eIF4E went through an initial round of selection performed with IMACs, followed by 2 rounds of in-solution selection using flow cytometry to enrich the population for high affinity eIF4E binders, where biotinylated eIF4E was detected using dye-labelled streptavidin. Insets show the enrichment in eIF4E cap-binders in the PELE library after rounds 2 and 3 of FACs selection. Negative control experiments were performed with the same library inputs that showed negligible non-specific binding within the enriched populations in the absence of eIF4E. (B) Samples from the final round input for YSD selection were co-incubated with either m7GTP, 4E-BP14ALA9 or VH-M49 in order to compete with the VH-DIF population enriched for eIF4E binding with biotinylated eIF4E (measured in A). A significant reduction in the enriched population interacting with eIF4E only occurred with m7GTP treatment indicating that the selected eIF4E binders were specific for the cap-binding site. (C) The table lists the 10 unique VH-DIF sequences identified from the 34 yeast clones sequenced in the final round of YSD selection, with their corresponding frequencies. A recognition motif was generated from the identified sequences using MEME,43 which in addition to showing the invariant cyclic peptide interaction motif, also identified that proline was preferentially enriched for at the position immediately preceding the motif. (D) Complex structure of eIF4E (green) with VH-DiFCAP-01 (cyan) highlighting the binding of the PELE selected motif presenting linkers to eIF4E. The 2Fo–Fc electron density map of the cap interacting loop structure is shown in blue (1.2σ). (E) Overlay of the cap binding motif of VH-DiFCAP-01 (E103MGFF107, shown in cyan and associated waters in red) with the equivalent residues in EE-02 (shown in white and associated waters in blue) highlighting the loss of the water mediated interactions between VH-DiFCAP-01 and E-103 and a small conformation change in E103, where the interaction with R112 and the structured water network are retained. However, it does result in an additional interaction with K162 not observed in the eIF4E:EE-02 complex. (F) The cap binding motif of VH-DiFCAP-01 (E103MGFF107, shown in magenta) forms a similar β-hairpin-like structure to that seen in the eIF4E:EE-02 complex. Additionally, the two intra backbone hydrogen bonds that formed to stabilize the bound structure of the EE-02 cyclic peptide (Fig. 3C) are also observed in the VH-DiFCAP-01 complex with eIF4E. (G) The hydrogen bond formed between C2 of EE-02 (cyan) with R157 is not observed in the eIF4E:VH-DiFCAP-01 complex, where it is replaced with a hydrogen bond between P102 (white) of the 100PLP102 linker and N50. (H) The PELE selected linker (T108NIPAMV114) forms 3 hydrogen bonds with eIF4E, where residues T108 and Q109 form 2 hydrogen bonds with the indole group of W102 of eIF4E (3.7 Å and 3.1 Å, respectively), and a hydrogen bond forms between the amide and carbonyl groups of A112 and eIF4E's A204. (I) Residues A112 and M113 of the linker region (T108NIPAMV114) form multiples hydrophobic contacts with residues T203, A204, T206 and G208 in eIF4E, which stabilize the α-helical secondary structure of the eIF4E region 201 to 205. (J) The conformation of E103MGFF107 (highlighted magenta) is stabilized by a hydrophobic cluster principally formed by I110 from the linker region (T108NIPAMV114) and a salt bridge between VH-DiFCAP-01 residues R51 and D36, which also interact with two buried structured waters. The water network in conjunction with R51 form hydrogen bonds with the polypeptide backbone of the PELE selected loop, helping to stabilize the conformation of the cap-site interaction motif for eIF4E binding. |
Initial designs where the cyclic eIF4E interaction motif (EMGFF) was introduced at different positions within a 15 mer loop in the CDR3 region of the VH-DIF scaffold were tested and none exhibited any binding to eIF4E (Fig. 4A and Fig. S3, ESI†). It was then hypothesized that the use of poly-glycine linkers for epitope presentation was too permissive and did not restrain the motif sufficiently in the correct structural conformation to be able to interact with eIF4E. To overcome this issue an alternative approach termed Peptide Epitope Linker Evolution (PELE) was adopted, wherein the linker regions were randomized, and a yeast surface display (YSD) library generated to select for linkers that optimally displayed the eIF4E interaction motif (Fig. 4A). To confirm that the PELE selection was successful, input samples of the enriched yeast library from the final selection round were used in competition experiments with either m7GTP, 4E-BP14ALA9 or VH-M49 (a VH domain that binds at the eIF4G interaction site of eIF4E) against biotinylated eIF4E, which confirmed that the binding of the selectants to eIF4E were specific to the cap-binding site (Fig. 4B). 34 yeast clones from the final round of YSD were sequenced, which yielded 10 unique VH-DIF peptide aptamers (termed VH-DiFCAP) (Fig. 4C). Sequence analysis revealed that eIF4E interactors were isolated from each of the PELE libraries used in the selection and that proline was preferentially selected for at the amino acid position preceding the interaction motif (Fig. 4C). The VH-DIFCAP peptide aptamers where then tested for bacterial expression, where upon those with good expression levels (VH-DIFCAP-01, VH-DIFCAP-02, VH-DIFCAP-05, VH-DIFCAP-09) were purified and screened for binding against eIF4E using the m7GTPFAM competition assay (Fig. S4, ESI† and Table 2). The VH-DiFCAP peptide aptamers that demonstrated binding in the competition assay including the constrained peptide EE-02, were then re-measured using ITC in direct binding titrations, which identified VH-DiFCAP-01 as the most potent eIF4E binder with a Kd of 35.3 ± 17.0 nM (Table 2). A Kd approximately equivalent to that determined for the constrained peptide EE-02.
Sequence | Competitive anisotropy (m7GTPFAM, kd, nM) | ITC | ||||
---|---|---|---|---|---|---|
K d (nM) | ΔH (kJ mol−1) | ΔS (kJ mol-1) | N | |||
EE-02 | NH2-AEMGFFQDG-NH2 | 406.2 ± 3.6 | 63.9 ± 18.8 | 25.5 ± 2.1 | −16.4 ± 3.1 | 1.3 ± 0.1 |
EE-02 M4A | NH2-AEAGFFQDG-NH2 | No binding | ND | ND | ND | ND |
VH-DiF CAP -05 | MHPSAICQDC–––– | No binding | 2152.5 ± 184.5 | −7.5 ± 1.4 | 83.4 ± 5.4 | 0.99 ± 0.1 |
VH-DiF CAP -09 | ––––MVPEPGLPSP | 3630 ± 300 | 421.9 ± 82.6 | −47.1 ± 2.6 | −35.6 ± 10.0 | 1.0 ± 0.1 |
VH-DiF CAP -01 | ––––PLPTNIPAMV | 501.0 ± 58.0 | 36.3 ± 17.0 | −31.0 ± 1.7 | 39.5 ± 10.0 | 1.0 ± 0.1 |
VH-DiF CAP -02 | ––PLYAPHVHHL-- | No binding | ND | ND | ND | ND |
VH-DiF CAP -01 MA | ––––PLPATNIPAMV | No binding | ND | ND | ND | ND |
The evolved linker regions of VH-DIFCAP-01, in addition to presenting the EE-02 interaction epitope optimally to interact with eIF4E, also forms multiple additional interactions with eIF4E (Fig. 4G–I). This contrasts sharply with the EE-02 cyclic peptide where only the Cys2 carbonyl group forms a hydrogen bond directly with R157 outside the residues critical for interacting with eIF4E. This hydrogen bond does not occur in the eIF4E:VH-DIFCAP-01:eIF4E complex structure where it is mimicked by a hydrogen bond with the carbonyl of the P102 in the P100LP102 linker region, which also forms another hydrogen bond to the side chain of N50 of eIF4e (Fig. 4G). The remainder of the N-terminal PLP linker forms no other interactions with eIF4E. However, the linker section (T108NIPAMV114) that occurs at the C-terminal of the interaction motif in VH-DIFCAP-01, forms most of the interactions made between eIF4E and the evolved linker regions. Residues T108 and Q109 form 2 hydrogen bonds with the indole group of W102 of eIF4E (3.7 Å and 3.1 Å, respectively). An additional hydrogen bond between the linker region and eIF4E is formed between the amide and carbonyl groups of A112 and eIF4E's A204, respectively (Fig. 4H). The remaining residues in the linker (110–114) apart from I110 and V114 make a multitude of hydrophobic contacts with residues T203, A204, T206 and G208 in eIF4E, resulting in stabilization of the α-helical secondary structure of this region of the protein (Fig. 4I). V114 and I110 in contrast are involved in interactions with the invariant part of the VH-DIF scaffold. VH-DIFCAP-01 also interacts weakly at a second site with eIF4E that is constituted from the CDR1 and CDR2 loops of the VH domain, where the residues S34 and S56 both form hydrogen bonds via their sidechains with the carbonyls of K52 and S53, respectively (Fig. S5, ESI†).
However, the most remarkable feature is the presence of two buried structure water molecules that allow the R51:D36 salt bridge to stabilize the fold of the CDR3 loop, when bound to eIF4E (Fig. 4J). These waters enable R51 and D36 to indirectly stabilize the polypeptide backbone of residue Ile110 and Q109 (Fig. 4J). Additionally, R51 also forms direct interactions with the carbonyl of F107 and the hydroxyl of the T108 sidechain, further rigidifying the presentation of the cap-site interaction motif (Fig. 4J). MD simulations also demonstrated that the unbound CDR3 loop of the VH-DiFCAP-01 domain undergoes a distinct structural change in comparison to the bound form (Fig. S2C and D, ESI†). The CDR3 loop undergoes a structural relaxation, whereby the β-hairpin structure associated with the ‘EMGFF’ motif is lost and instead there is a general movement of the CDR3 loop away from the body of the VH domain scaffold (Fig. S6, ESI†). Interestingly, this movement is underpinned by significant structural changes in the hydrophobic core of the CDR3 loop structure with the L110 sidechain rotating out and the M113 sidechain rotating in to replace it. In association with these changes, the two buried structured waters observed in the bound form also adopt new positions within the CDR3 loop, which help to stabilize the new conformation by forming 2 water mediated interactions between the amide backbones of Q109 and G98 with the D36 sidechain, respectively (Fig. S6, ESI†). In both simulations of the bound and unbound forms of VH-DiFCAP-01 (Fig. S7, ESI†), the two water positions remained predominately solvated indicating that water molecules found here exchanged slowly with the external solvent and formed enthalpically favorable interactions.
Fig. 5 (A) Anti-FLAG IP pull down of HEK293 cells transfected with either VH-DIFCAP-01, VH-DIFCAP-02, VH-DIFCAP-01 MA and VH-S4, a VH-domain that interacts with eIF4E at the eIF4G binding interface with high affinity.9 IP experiments were performed 24 hours post transfection. Whole cell lysate (WCL) was also blotted for the corresponding proteins and is shown on the left of the blot. (B) m7GTP pulldown of eIF4E containing complexes from HEK293 transfected with VH-DIFCAP01, VH-DIFCAP02, VH-DIFCAP01 MA (M104A) and VH-S4. Whole cell lysate (WCL) was also blotted for the corresponding proteins and is shown on the left of the blot. In the blot below an equivalent pull-down was performed but with the HEK293 cells transfected with increasing amounts of expression vector. (C) HEK293 cells were transfected with either empty vector, VH-DIFCAP01 or VH-DIFCAP01 MA (M104A), and eIF4E phosphorylation and cyclin D1 expressions levels assessed via western blot. Actin was used as a loading control, whilst anti-FLAG was used to assess expression of the transfected proteins. Protein levels were assessed 48 hours post transfection. (D) A bicistronic luciferase reporter, which measures the relative amount of cap-dependent translation (Renilla) to cap-independent translation (Firefly), was co-transfected with either empty vector (MOCK) or increasing amount of VH-DIFCAP-01, VH-DIFCAP-01 MA, VH-S4 plasmid vector into HEK293 cells (see material and method). Renilla and Firefly luciferase activity was measured 48h post transfection and plotted as a ratio-metric value. Statistical significance of differences in the ratio-metric values were determined by non-paired Student's t test. All statistical tests were two-sided and variance assumed to be equal. P-val < 0.05 was considered statistically significant. (E) Anti-His IP pulldown of purified VH-DIFCAP-01 exogenously added to HEK293 cell lysates either treated with CGP57380 or vehicle control. Input lysate is shown on left hand side of the western blot. |
The LgBIT-eIF4E and SmBIT-VH-DIFCAP-01 were identified as the transfection pair that reconstituted the highest luciferase signal without exhibiting background activity in the negative controls, thus confirming that the NanoBit reporter fragments were not spontaneously assembling under the experimental conditions (Fig. 6A and B). To further confirm the specificity of this assay for binding at the cap site, the NanoBit assay was re-performed with both the VH-DIFCAP-01M104A (termed VH-DIFCAP-01 “MA”) and VH-DIFCAP-01E103A (termed VH-DIFCAP-01 “EA) binding controls fused to smBIT, which resulted in the abrogation of the luciferase signal above background (Fig. 6B). Additionally, the ability of the assay to measure and differentiate between interactors that bound at the either the cap-binding or eIF4G binding sites of eIF4E was demonstrated by co-transfection of the NanoBit assay (termed NanoBIT eIF4ECAP) with either untagged VH-M49, a VH-domain designed to interact with eIF4E at the eIF4G binding site, or VH-DIFCAP-01 (Fig. 6C).
Fig. 6 (A) Inset shows how the interaction of proteins A and B fused to SmBiT and LgBiT (Promega, USA) enables reconstitution of the active NanoBit (Promega, USA) luciferase. Graph shows the reconstituted luminescence activity of the various combinations of either eIF4E or VH-DIFCAP-01 fused at either the N or C terminal of SmBiT and LgBiT, respectively, co-transfected into HEK293 cells. Individual N- and C-terminal LgBiT-linked eIF4E and VH-DIFCAP-01 constructs co-transfected with SmBiT-HALO served as negative controls. (B) To validate the specificity of the SmBIT-VH-DIFCAP-01 and LgBiT-eIF4E interaction pair, two VH-DIFCAP-01 point mutant controls were generated (E103A (EA) and M104A (MA), respectively) and co-transfected into HEK293 cells with LgBiT eIF4E, which resulted in loss of bioluminescence. Inset: Cell samples replicating the NanoBiT experimental conditions were assessed for their relative levels of LgBIT fused eIF4E to endogenous eIF4E (anti-eIF4E), and expression levels of the various SmBIT-VH-DIFCAP-01 constructs (anti-MYC). (C) (Right hand graph) The ability of the SmBIT-VH-DIFCAP-01:LgBiT-eIF4E (termed NanoBITCAP) interaction pair to discriminate between different classes of eIF4E binders was tested by co-expressing it with either VH-S4 (a VH domain that interacts specifically with the eIF4G interaction site) or VH-DIFCAP-01 not fused to SmBIT, where only VH-DIFCAP-01 caused a decrease in luminescence. (Left hand graph) The specificity of VH-DIFCAP-01 was further investigated by co-expressing either VH-DIFCAP-01 or VH-S4 with the NanoBit eIF4E:eIF4G604–646 system,41 which measures binding at the eIF4G interface and demonstrated that VH-DIFCAP-01 only interacts with the cap-binding interface. (D) HEK293 cells were transfected with the NanoBITCAP system and permeabilized with a sub-CMC (critical micelle concentration) concentration of digitionin. Cells were then treated with different titrations of small molecules that either specifically targeted the cap (m7GTP, m7GDP) or eIF4G (4EGI1) binding interfaces of eIF4E. (E) HEK293 cells were transfected with the NanoBit eIF4E:eIF4G604–646 system41 and again permeabilized with a sub-CMC (critical micelle concentration) concentration of digitionin. Cells were then treated with titration of the following compounds (m7GTP, m7GDP and 4EGI) to assess the specificity of the NanoBITCAP system. |
To validate that the NanoBIT eIF4ECAP system can measure small molecule mediated inhibition of the eIF4E cap-binding site, the system was also used to screen two known cap-analogue antagonists (m7GTP and m7GDP) and an established cell permeable inhibitor of the eIF4E:4G interface (4EGI1)24 as a negative control. Unfortunately, both cap-analogue molecules are cell impermeable. Therefore, to circumvent this issue, a sub-CMC (critical micelle concentration) treatment of digitonin was used to permeabilize and enable cellular entry of the cap-analogues into HEK293 cells transfected with the NanoBITCAP system. In the permeabilized cells, both m7GTP and m7GDP disrupted the interaction of LgBIT-eIF4E with SmBIT-VH-DIFCAP-01 with IC50s of 12.8 μM and 34.5 μM respectively, whilst 4EGI124 had a negligible effect on the NanoBIT signal, demonstrating that the assay system can measure cap-binding site antagonists (Fig. 6D). To highlight the specificity of the NanoBITCAP system further, it was also shown that only 4EGI124 and neither of the two cap-analogues were able to inhibit the luciferase signal in an alternative NanoBIT system (eIF4E:eIF4G604–646)41 that measures disruption of the eIF4E:4G interface (Fig. 6E). Both sets of described experiments were then repeated in non-permeabilized cell, where as expected the impermeable cap-analogues elicited no effects, and the cell permeable only 4EGI1 inhibited the NanoBIT signal in the NanoBit eIF4E:eIF4G604–646 system.41
An interesting feature of the presentation of the EE-02 binding motif in VH-DiFCAP-01 is the role played by two buried structured waters that facilitate the transition of the VH domain from its unbound form to the bound form with eIF4E. These waters in combination with the linkers isolated in the yeast selection allow the CDR3 loop to adopt the β-hairpin structure that the “EMGFF” epitope requires to bind to eIF4E and to accommodate the D38:R51 salt-bridge that the CDR3 loop sit upon. The CDR3 loop of VH-DiFCAP-01 folds back onto the former light-chain interaction surface, where the D38:R51 salt-bridge is located, in a manner highly like that seen in nanobodies (Fig. S9, ESI†). Also like many other VH and VHH structures and their interaction with their binding partner, the interaction with eIF4E is principally mediated through the CDR3 loop. However, this type of interaction does differ significantly from the reported VH domain interaction with VEGF, where both the CDR3 and the former light-chain interaction surface are involved in macromolecular recognition (Fig. S10, ESI†).42 The identified cap-binding site interaction pose can also be used as a template for the development of new therapeutic modalities. Furthermore, to support the discovery of drug-like cap-analogue molecules we have also used the VH-DiFCAP-01 peptide aptamer to develop a split-luciferase assay that can assess cellular target engagement of the cap binding site in live-cells. In addition, VH-DiFCAP-01 can also be used to further study the biological role of the eIF4F complex, especially in combination with the reported VH-S4 domain that disrupts eIF4F complex formation by binding eIF4E at its 4G interaction site. A suite of tools enabling discrimination of the effects of inhibition at precise sites on eIF4E in cellular studies. This contrasts with tools such as RNAi where mRNA levels are reduced and do not allow the study of specific functional domains within proteins, and in turn any insights may be therapeutically limited.
We also envision that the VH-DIFCAP-01 molecule and other such molecules, when combined with new and novel alternate delivery methods e.g., RNA, toxin, and cell-penetrating peptide mediated delivery methods, could also be used to target cap-dependent translation and constitute a potential therapeutic. These results demonstrate that the PELE process can enable the development of mini-proteins that interact at desired target interaction sites, that can model and assess the potential effects of drug inhibition and allow the construction of critical target engagement assays that can accelerate the identification of lead compounds for drug development.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d2cb00099g |
‡ Authors contributed equally. |
This journal is © The Royal Society of Chemistry 2022 |