Protein grafting of p53TAD onto a leucine zipper scaffold generates a potent HDM dual inhibitor

Abstract

Reactivation of the p53 pathway by a potential therapeutic antagonist, which inhibits HDM2 and HDMX, is an attractive strategy for drug development in oncology. Developing blockers towards conserved hydrophobic pockets of both HDMs has mainly focused on small synthetic compounds; however, this approach has proved challenging. Here we describe an approach to generate a potent HDM dual inhibitor, p53LZ2, by rational protein grafting of the p53 transactivation domain onto a homodimeric leucine zipper. p53LZ2 shows tight binding affinity to both HDMs compared with wild-type p53 in vitro. X-ray crystallographic, comparative modelling and small-angle X-ray scattering studies of p53LZ2-HDM complexes show butterfly-shaped structures. A cell-permeable TAT-p53LZ2 effectively inhibits the cancer cell growth in wild-type but not mutant p53 by arresting cell cycle and inducing apoptosis in vitro. Thus, p53LZ2, designed by rational grafting, shows a potential therapeutic approach against cancer.

ARTICLE

Received 4 Oct 2013 | Accepted 7 Apr 2014 | Published 7 May 2014

Jung-Hoon Lee1, Eunji Kang1, Jungmin Lee1, Jungmin Kim1, Kyoung Hu Lee1, Jieun Han1, Hye Yoon Kang1, Soshin Ahn1, Youngmi Oh1, Dongkyu Shin1, Kyeyeon Hur1, Su Young Chae1, Paul H. Song1, Yong-In Kim1,2, Jae Chan Park2 & Jae Il Lee1

Reactivation of the p53 pathway by a potential therapeutic antagonist, which inhibits HDM2 and HDMX, is an attractive strategy for drug development in oncology. Developing blockers towards conserved hydrophobic pockets of both HDMs has mainly focused on small synthetic compounds; however, this approach has proved challenging. Here we describe an approach to generate a potent HDM dual inhibitor, p53LZ2, by rational protein grafting of the p53 transactivation domain onto a homodimeric leucine zipper. p53LZ2 shows tight binding afnity to both HDMs compared with wild-type p53 in vitro. X-ray crystallographic, comparative modelling and small-angle X-ray scattering studies of p53LZ2-HDM complexes show buttery-shaped structures. A cell-permeable TAT-p53LZ2 effectively inhibits the cancer cell growth in wild-type but not mutant p53 by arresting cell cycle and inducing apoptosis in vitro. Thus, p53LZ2, designed by rational grafting, shows a potential therapeutic approach against cancer.

DOI: 10.1038/ncomms4814

Protein grafting of p53TAD onto a leucine zipper scaffold generates a potent HDM dual inhibitor

1 Bio Therapeutics Laboratory, Samsung Advanced Institute of Technology (SAIT), Samsung Electronics Co. Ltd., Suwon 443-803, South Korea. 2 Bio Research Center, Samsung Advanced Institute of Technology (SAIT), Samsung Electronics Co. Ltd., Suwon 443-803, South Korea. Correspondence and requests for materials should be addressed to J.C.P. (email: mailto:[email protected]

Web End [email protected] ) or to J.I.L. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 5:3814 | DOI: 10.1038/ncomms4814 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4814

The cellular level of p53 in unstressed normal cells is tightly controlled by an auto-regulatory feedback loop that regulates the activity of p53 protein and the expression

level of a specic gene encoding a potent p53-negative regulator HDM2 (ref. 1). The transcription factor p53 binds to the E3 ubiquitin ligase HDM2 gene and the elevated intracellular level of HDM2 protein in turn decreases the endogenous p53 by a direct interaction between the N-terminal domains of both proteins resulting in proteasome-mediated degradation of polyubiquitinated p53 (ref. 2). While several previously identied E3 ubiquitin ligases out of HDM2 (for example, Pirh2, COP1 and ARF-BP1)35 target p53 for ubiquitin-mediated degradation, there is no direct evidence that any of these ligases can replace HDM2 to regulate p53 stability in cancer cells6. On the other hand, a close relative of HDM2, HDMX (also called HDM4), is known to control p53 status in cancer cells through a direct interaction that neutralizes the transcriptional activation potency of the transactivation domain of endogenous p53 (hereafter p53TAD) in a similar way as observed in the p53-HDM2 interaction. HDMX per se does not act as the E3 ubiquitin ligase of p53, but rather enhances HDM2 ligase activity mediated by the C-terminal RINGRING interaction between HDM2 and HDMX, albeit the function of HDM2-HDMX heterodimerization has not been completely claried yet710. Interestingly, HDM2 is not only a bona de E3 ubiquitin ligase for p53 but also ubiquitinates HDMX11. Thus, HDM2 serves as a key regulator that orchestrates HDMX stability as well as p53 activity in vivo.

In terms of a therapeutic point of view, tumour suppressor p53 is a potent anti-proliferative, anti-apoptotic guardian molecule that induces cell cycle arrest, apoptosis, autophagy and senescence in cancer cells12. Almost 11 million patients suffering from various cancer types worldwide are known to bear wild-type p53 in tumours13. Thus, reactivation of the p53 pathway by displacing the negative regulator HDM2 and/or its homologue HDMX from endogenous wild-type p53 in cancer cells has long been considered an important therapeutic strategy in oncology. Accordingly, numerous innovative efforts from academic laboratories as well as pharmaceutical industries have been

focused on developing synthetic compounds, peptidomimetics or miniature protein antagonists of the p53HDM2 (and HDMX) interaction1423. Among them, the rst reported HDM2-specic inhibitor Nutlin-3, which mimics three essential hydrophobic residues of wild-type p53 (F19, W23 and L26; human p53 numbering) required for HDM2 recognition, demonstrated that this small compound disrupting the p53HDM2 interaction recovered p53 function and resensitized cancer cells to apoptosis in vitro and in vivo18. In fact, the p53TAD is largely unstructured under physiological conditions, while transforming into an a-helical conformation when interacting with target proteins in the presence of binding partners such as RPA and HDMs24. Recent studies of a stapled p53-based peptide, which contains a chemical hydrocarbon crosslink within a synthetic peptide (stapling), demonstrated that the p53 peptidomimetic restored the intrinsic a-helix-forming tendency of p53TAD and conferred the intracellular protease resistance2527. Such seminal structural insights into the stabilized a-helix (SAH) of p53 inspired us to challenge this attractive target in oncology and design a new molecular entity of a p53TAD mimetic.

Here we describe a new approach in generating a potent HDM2 and HDMX dual antagonist by rational protein grafting of p53TAD onto the GCN4 leucine zipper domain (hereafter GCN4LZ), a well-known homodimeric a-helical scaffold having a completely helical structure extensively stabilized by the intra-chain salt bridges as well as the inter-chain hydrophobic interactions28,29. By transferring the three key hydrophobic residues of p53TAD, F14, W18 and L21 (new numbering in this study, see Fig. 1a) onto the ideal positions of each helix in the GCN4LZ homodimer and also adding two more substitutions, G10 and S15 (new numbering), which further ne-tune the binding afnity of a designed p53TAD mimetic against HDM2 and HDMX, we generate a potent dual inhibitor with high binding afnity against both HDMs in comparison with wild-type p53 in vitro. We also report the X-ray crystallographic, comparative protein structure modelling and small-angle X-ray scattering (SAXS) studies of coiled-coil architecture of the GCN4LZ-driven p53TAD mimetic (hereafter p53LZ) and its complex with the N-terminal domain of either HDM2 or HDMX.

a

Heptad position

1st Repeat 2nd Repeat 3rd Repeat 4th Repeat

Residue No. GCN4LZ p53LZ1 p53LZ2 p53LZ3

1

R R R R

a 2 M M M M

b 3 K K K K

c 4 Q Q Q Q

d 5 L L L L

e 6 E E E E

f 7 D D D D

g 8 K L K K

a 9 V V V V

b 10 E E G S

c 11 E E E E

d 12 L L L L

e 13 L L L L

T Q

f 14 S F F F

g 15 K K S A

a 16 N N N E

b 17 Y Y Y Y

Y L L

c 18 H W W W

W W W

d 19 L L L L

e 20 E E E L

f 21 N L L L

L L L

g 22 E E E S

S X P

a 23 V V V V

b 24 A A A T

Q E

c 25 R R R R

d 26 L L L L

e 27 K K K K

f 28 K K K K

g 29 L L L L

a 30 V V V V

pMI SAH-p53-8 p53TAD

Q

SS E

S T T

Q

F F F

A X S

E N D

A R K

L L L

Q

N N

b

R H

E Q

L E K K

V

V N

V

K

E

L E

A

Y E K

Q

E H

R W

E Q

L E K K

V V

N

K

E

L E

A

Y E K

Q

E W

R W

E Q

L E S K

V

V N

K

E

L E

A

Y G

K

Q

E W

R W

E Q

L S A K

V

V E V

K

L L E

T

Y S K

Q

E W

L

L L L

M

V N

V V

L

L L L

M

V N

V V

L

L L L

M

V N

V V

L

L L L

M

V E V V

K N

S D f

b e

K

LF D

K

LF D

K

LF D

c g

d

a d

L g c

e

g

b

c

g

c

V e c g

d

a

b

e b

f

a

d a

V e b

f

d

a

a

d

M

L L L

M

L

L L L

M

L L L

M

L

L L L

f

D

N

K

S

f

D

L K

F

f

D

L K

F

f

D

L K

F

d

d

b e

a

g

c

b e

a

L g c

f

b

c

K E Y

A

e

g

E

L E K

K

K E

L

K E Y

A

E

L E K

K

K E

L

K G

Y A

E

L E K

K

S E

L

K S Y

T

E

L L K

K

A S

L

R

R

R

R

GCN4LZ p53LZ1 p53LZ2 p53LZ3

Figure 1 | p53LZ design concept. (a) Sequence alignment of GCN4LZ and designed p53LZs. Two different peptides showing high afnities against HDM2 and HDMX (pMI and SAH-p53-8) together with naturally occurring WT p53 (p53TAD) aligned to GCN4LZ and designed p53LZs. Three key hydrophobic residues, F, W and L, are highlighted in blue. Additional substitutions of amino acids in p53LZ1, p53LZ2 and p53LZ3 are shown in green. Each heptad repeat is marked at the top of GCN4LZ sequence. (b) Helical wheel representations of GCN4LZ and grafted p53LZs.

2 NATURE COMMUNICATIONS | 5:3814 | DOI: 10.1038/ncomms4814 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4814 ARTICLE

Further biochemical and cell-based studies exhibit an impressive on-target biological efcacy of a cell penetrating p53LZ (TAT-p53LZ2) in various cancer cell lines expressing wild-type but not mutant p53. Reactivation of the p53 pathway by our designed p53LZ2 shows a potential therapeutic possibility against cancer and provides a framework to advance the development of peptide-based biologics in drug discovery technology.

ResultsRational protein grafting of p53TAD onto the GCN4LZ. The GCN4LZ encompassing 30 residues shows the characteristic seven-residue (heptad) repeat (af)n, where hydrophobic or charged residues are positioned at the (a, d) and (e, g) slots, respectively (Fig. 1a). The GCN4LZ core, which consists of (a, d) and (e, g) positions, is not easily accessible by its binding partners, whereas the (c, h, f) side-chains protruding from the GCN4LZ surface can be considered as desired areas for p53TAD grafting because of the large solvent accessible area of homodimeric GCN4LZ (PDB ID 2ZTA, B5,101 2). Therefore, these latter residues were primarily tested for the grafting of p53TAD onto the GCN4LZ (Fig. 1b). The side-chain orientations of the three key hydrophobic residues, F19, W23 and L26 (human p53 numbering), of p53 directs towards the conserved hydrophobic pockets on HDM2 and HDMX, showing complementarity to the HDM pockets30,31. We chose (c, f) positions of the 2nd and 3rd heptad repeats (a9g22) to graft F, W and L residues onto the GCN4LZ because disturbing the 1st and 4th (trigger site of the GCN4LZ) repeats decreases the overall protein stability more than other repeats (Fig. 1b). These grafted p53TAD mimetics, which are named p53LZ1, p53LZ2 and p53LZ3, contain 3, 5 and 9 mutation points, respectively. p53LZ1 has the simplest amino-acid substitution on the GCN4LZ, whereas p53LZ3, bearing more a complicated substitution, is designed to be more similar to the pMI sequence (TSFAEYWALLS), one of the previously identied high-afnity peptides against the HDM2 and HDMX17. Two additional changes onto the p53LZ2, E10G and K15S, in addition

to the three key hydrophobic residues, were carefully identied to strengthen the proteinprotein interaction between the designed p53LZ and HDM proteins based on the structural analysis by comparative structure modelling. We also calculated the free energy of binding of the designed p53LZ to the HDM2 N-terminal domain by using the available X-ray structures (PDB ID 1YCR and 2ZTA) and GROMACS32 software to simulate the possible binding mode of each complex. Molecular dynamic simulations coupled with MM/PBSA (Molecular Mechanics/PoissonBoltzmann Surface Area) approach were carried out for 60 ns. The calculated binding energy (DGcomp)

of the p53LZ2HDM2 complex exhibited the lowest value of

41.5 kJ mol 1, whereas the HDM2 N-terminal domain in complex with either p53LZ1 or p53LZ3 revealed a relatively higher binding energy ( 9.0 or 6.9 kJ mol 1) compared with the

p53LZ2-HDM2 complex, suggesting that the E10G and K15S substitutions in the p53LZ2 sequence may contribute to the higher afnity binding (Supplementary Table 1). Two mutations, E20L and E22S, in p53LZ3, which originally interacts with K94 and K51 of HDM2, respectively, likely impair the p53LZ3-HDM2 interaction by disrupting electrostatic interactions between both proteins, revealing signicant differences of the calculated coulombic (DGcoul 710.8 kJ mol 1) and polar solvation

(DGps 262.6 kJ mol 1) terms, relative to those of p53LZ1

and p53LZ2 (Supplementary Fig. 1 and Supplementary Table 1).

Characterization of HDM2 and HDMX interaction with p53LZ. To assess the binding of HDM2 and HDMX to p53LZ in vitro, we expressed and puried recombinant GST-p53LZ proteins (p53LZ1, p53LZ2 and p53LZ3) and N-terminal domains of His6-HDM2 and His6-HDMX corresponding to residues 1125 and 1132, respectively (hereafter HDM2NTD and HDMXNTD) and then performed an in vitro GST pull-down assay (Fig. 2a). GST-p53LZ1 and GST-p53LZ2 effectively pulled down both HDM2NTD and HDMXNTD, whereas GST and GST-p53LZ3 did not show a relevant interaction with either

a b

800

600

400

200

37 25

20

15 10

GST-p53LZ3 GST-p53LZ2 + His6-HDM2NTD GST-p53LZ2 + His6-HDMXNTD

GST-p53LZ2

Normalized A 280(mAU)

+ +

+ +

+ +

GST-p53LZ1

His6-HDM2NTD His6-HDMXNTD

+

+

+

+ His6-HDMXNTD

+ + + +

His6-HDM2NTD

+

+

GST

GST-p53LZ2

His6-HDMXNTD

His6-HDM2NTD

37

25 20

15 10

GST-p53LZ2 GST-p53LZ1 GST

His6-HDMXNTD His6-HDM2NTD

GST-p53LZ3 GST-p53LZ2

GST-p53LZ1 GST

His6-HDMXNTD His6-HDM2NTD

Input Precipitate

37

25 20

15

10

Lane 1 2 3 4 5 6 7 8

GST-p53LZ3

0 0 2 4 6

Retention volume (ml)

8 10 12 14 16 18

c

His6-HDM2NTD + p53LZ2 Mw = 39.8 kDa

His6-HDMXNTD + p53LZ2 Mw = 40.6 kDa

Relative light scattering

1.2

0.8

0.4

80,000

40.2 kDa (0.12%) 1.2

0.8

0.4

80,000

Relative light scattering

41.6 kDa (0.07%)

60,000

60,000

40,000

Mw (Da)

40,000

Mw (Da)

20,000

20,000

0

8.7 8.9 9.1 9.3 Retention volume (ml)

9.5 9.7

0

0

8.7 8.9 9.1 9.3 Retention volume (ml)

9.5 9.7

0

Figure 2 | p53LZ1 and p53LZ2 interact with both HDM2 and HDMX in vitro. (a) In vitro GST pull-down assay of designed p53LZs against both HDM2 and HDMX N-terminal domains (HDM2NTD and HDMXNTD). GST or GST-tagged proteins were incubated with HDM2NTD or HDMXNTD and then the protein mixtures were precipitated with glutathione resin after washing to remove unbound proteins. (b) Gel ltration analysis on GST-p53LZ2 in the absence and presence of HDM2NTD or HDMXNTD. The retention volumes and SDSPAGE results are indicated. (c) Multi-angle light scattering (MALS) measurement of p53LZ2 bound to HDM2NTD or HDMXNTD showing the relative light scattering signal as a functionof elution volume. The measured molecular mass of each peak is shown in blue.

NATURE COMMUNICATIONS | 5:3814 | DOI: 10.1038/ncomms4814 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4814

HDM2NTD or HDMXNTD, demonstrating a direct physical interaction of HDM2 and HDMX with the designed p53LZ1 and p53LZ2 in vitro. Apparently, GST-p53LZ2, compared with GST-p53LZ1 binds preferentially to both HDM2NTD and HDMXNTD, suggesting that the two additional substitutions in the p53LZ2, E10G and K15S contribute to a tight interaction with HDM2 and HDMX (Fig. 3). Further binding studies with analytical size-exclusion chromatographic (SEC) analysis demonstrated a direct interaction between p53LZ2 and both HDM N-terminal domains, revealing co-migration of p53LZ2 and its binding partner in gel ltration (Fig. 2b). To precisely determine the stoichiometry between p53LZ2 and both HDM2NTD and HDMXNTD, we used multi-angle light scattering (MALS) analysis in conjunction with refractive index, which measures molecular mass independent of its molecular

shape. The result indicates the existence of a homogeneous species for p53LZ2HDM2NTD or p53LZ2HDMXNTD heterotetrameric complex in solution with a corresponding molecular mass of 40.2 or 41.6 kDa, respectively, conrming a 1:2 stoichiometry between p53LZ2 homodimer and both HDM N-terminal domains (Fig. 2c).

In vitro binding of wild-type p53 and its mimetics to HDMs.

To further investigate the designed p53LZs in interaction with either HDM2 or HDMX, we assessed the binding of MBP-fused HDM proteins, full-length (FL) HDM2, HDM2NTD and HDMXNTD to GST, GST wild-type (WT) p53TAD, GST-p53LZ1 and GST-p53LZ2, puried from E. coli in ELISA. In addition, a well-known high-afnity variant of p53TAD (pMI)17

a d

Time (min)

50 60

FL HDM2 binding to WT p53TAD and its mimetics

1.0

0.5

0.00

0 10 20 30 40

Absorbance at 450 nm

GST-p53LZ2

cal s1

kcal mol1 of injectant

0.20

0.40

0.60

0.800.0

2.0

4.0

6.0

8.0

10.0

12.0

0.0 0.5

Molar ratio

Time (min) 0

0.0

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6

GST-p53LZ1

GST-WT p53

GST

GST-pMI

GST or GST proteins (M)

b

2.0

HDM2NTD binding to WT p53TAD and its mimetics

p53LZ2/HDM2NTD

N = 0.65

KD = 17.75 nMH = 11.76 kcal mol1

TS = 1.48 kcal mol1

Absorbance at 450 nm

1.5

1.0

0.5

GST-p53LZ2

1.0 1.5

10 20 30 40 50 60

0.0

0.00

0.50

1.00

1.50

2.0

4.0

6.0

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6

GST-p53LZ1

GST-WT p53

GST

GST-pMI

GST or GST protein (M)

cal s1

kcal mol1 of injectant

c

2.0

HDMXNTD binding to WT p53TAD and its mimetics

Absorbance at 450 nm

1.5

1.0

0.5

0.0

10.0

GST-p53LZ2

p53LZ2/HDMXNTD

N = 0.65KD = 80 nMH = 9.66 kcal mol1

TS = 1.16 kcal mol1

8.0

0.0

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6

GST-p53LZ1

GST-WT p53

GST

GST-pMI

GST or GST protein (M)

0.0 0.5 1.0 Molar ratio

1.5 2.0

Figure 3 | p53LZ2 exhibits high-afnity binding to both HDM2 and HDMX, relative to WT p53 in vitro. GST or GST proteins were used in ELISA assays of (a) MBP-FL HDM2, (b) MBP-HDM2NTD and (c) MBP-HDMXNTD binding. GST-p53LZ2 harbouring ve grafted residues, G10, F14, S15, W18 and L21, showed signicantly higher binding afnities to HDM2 and HDMX in comparison with WT p53TAD (residues 172) and GST-p53LZ1. All data are presented as means.d. of three independent experiments. (d) ITC titration and tting curves for p53LZ2 binding to HDM2NTD or HDMXNTD. The numbers show the tting of the data, with N (stoichiometry), KD (dissociation constant), DH (enthalpy change) and TDS (entropy change).

4 NATURE COMMUNICATIONS | 5:3814 | DOI: 10.1038/ncomms4814 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4814 ARTICLE

previously identied by phage display was also produced as GST-tagged proteins to compare its relative binding afnity to HDM2 and HDMX. The ELISA plate coated with MBP-FL HDM2, MBP-HDM2NTD or MBP-HDMXNTD was incubated with increased amounts of GST or GST-fused proteins. A horseradish peroxidase-conjugated anti-GST antibody was used to detect the bound protein. p53LZ2, harbouring a total of ve grafted residues, bound preferentially to both HDM2 and HDMX compared with others, suggesting that the two more grafting (G10 and S15) in addition to the essential hydrophobic residues (F14, W18 and L21) onto the GCN4LZ is critical for tight interaction with both HDM proteins (Fig. 3ac). Furthermore, the quantitative ELISA studies revealed a relatively similar binding of GST-p53LZ2 and GST-pMI to both HDM2 and HDMX. To quantify the target-binding strength of grafted p53LZ2 in vitro, we measured the direct binding afnity of p53LZ2 to HDM2NTD and HDMXNTD using isothermal titration calorimetry (ITC) (Fig. 3d). p53LZ2 strongly bound to HDM2NTD (KD 18 nM), exceeding a previously reported value

of WT p53TAD for HDM2NTD (KD 61.3 nM, measured in

PBS containing 0.1 mM TCEP and 0.01% NaN3)17, whereas the p53LZ2 homodimer exhibited a moderately higher afnity for HDMXNTD (KD 80 nM) in comparison to a previously

determined ITC data (KD 202 nM, measured in PBS containing

0.1 mM TCEP and 0.01% NaN3)17. The enthalpyentropy compensation in both cases was revealed by favourable enthalpy (DH) and unfavourable entropy ( TDS) changes,

which contribute to a high-afnity binding of p53LZ2 to both HDM N-terminal domains. In addition, the binding stoichiometry (N) determined by the ITC analysis further supports that homodimeric p53LZ2 preferentially interacts with both HDM proteins in 1:2 complexes (Fig. 3d).

Structural analysis of HDM2 and HDMX recognition by p53LZ2. Next, we performed three distinct experiments to investigate the structural mechanism of HDM2 and HDMX recognition by p53LZ2. First, we conducted circular dichroism (CD) experiments to evaluate the structural impact of protein grafting onto the GCN4LZ. The recombinant p53LZ2 and its parent molecule GCN4LZ displayed highly a-helical conformations (490% helicity) in physiological pH and almost identical Tm values (B63 C), suggesting that the protein grafting strategy of p53TAD onto this homodimeric scaffold to facilitate target binding did not alter the overall structural integrity of GCN4LZ (Supplementary Fig. 2). Second, we determined a high-resolution(1.2 ) crystal structure of p53LZ2 harbouring substitution in ve residues on the GCN4LZ (Table 1). The initial electron density map corresponding to homodimeric parallel coiled-coil architecture of the grafted p53LZ2 was readily generated in high quality by molecular replacement method using the GCN4LZ as search model (Fig. 4a and Supplementary Fig. 3). The rened high-resolution structure of p53LZ2 is a twisted elliptical cylinder protruding the key HDM binding residues F14, W18 and L21 from the main-chains of both helices. Furthermore, two parallel a-helices are associated with tightly packed conformations stabilized by the inter-chain salt bridges, R1E6, S15N16, N16N16 and E22K27 as well as the extensive inter-chain hydrophobic interactions, L5L5, V9V9, L12L12, L12L13, L19L19, V23V23, L26L26, L29V30 and V30V30. Interestingly, K15S mutation in the p53LZ2 likely compensated the loss of K15E20 interaction between the two helices of GCN4LZ (Fig. 1b) by creating a new inter-chain S15N16 interaction. To understand the HDM2 and HDMX recognition by p53LZ2 homodimer, we generated a complex structure for p53LZ2 bound to HDM2NTD or HDMXNTD using HADDOCK33. The

sophisticated structural superposition of WT p53TAD-HDM2 or HDMX and the HADDOCK structure were also conducted to evaluate proteinprotein docking. The key HDM-binding residues, F14, W18 and L21, of p53LZ2 were situated along one face of each a-helix and well positioned into the three hydrophobic pockets of both HDM2NTD and HDMXNTD, as similarly observed in the wild-type complexes (Fig. 4b,c). In contrast, E10G and K15S substitutions of p53LZ2 contributed towards reducing possible steric clash between adjacent side-chains of GCN4LZ and HDM proteins (for example, GCN4LZ E10HDM2 Q69, Q72 and GCN4LZ K15HDM2 M62) and optimizing spatial arrangement in p53LZ2-HDM interaction for enhanced binding. The p53LZ2HDM2 interaction is further reinforced by numerous electrostatic interactions between both proteins for example, p53LZ2 K3, E6, D7, E20, K27 and HDM2 E69, K70, R65, K94, E52, respectively (Supplementary Fig. 4), supporting the notion that a high-afnity binding of p53LZ2 to HDMs is likely attributed to large binding areas through the N-terminal and C-terminal segments of p53LZ2 in part. As a third experiment, given the available crystal structures and possible protein complex models, we co-puried homodimeric p53LZ2 in complex with either HDM2NTD or HDMXNTD fromE. coli and further investigated the resulting protein complexes using SAXS analysis3436. All SAXS measurements of the p53LZ2HDM complexes, in PBS buffer, were done in triplicate exposures and analysed to calculate the nal composite scattering curves (Fig. 4d, left). By the Guinier analysis of the SAXS scattering curves, we conrmed the absence of protein aggregation and further processed the data (Supplementary Fig. 5). The maximum diameter (Dmax) from the

pair distance distribution function p(r) and the radius-of-gyration Rg of p53LZ2HDM2NTD complex are B60 and B22 , respectively (Fig. 4d, right). The ab initio reconstitution of p53LZ2 in complex with either HDM2NTD or HDMXNTD using GASBOR37 represents the buttery-shaped SAXS envelope displaying pseudo twofold symmetry (Fig. 4e, Supplementary Figs 6 and 7). These reproducible envelopes show a good

Table 1 | Data collection and renement statistics.

p53LZ2 (native)

Data collection

Space group P21Cell dimensionsa, b, c () 30.48, 29.39, 32.70 a, b, g () 90.00, 98.25, 90.00

Resolution () 301.2 (1.231.2)* Rsym (%) 2.2 (5.5)

I/sI 59.73 (22.52) Completeness (%) 97.5 (88.5)

Redundancy 6.7

Renement

Resolution () 301.2No. of reections (total/test) 17,557/879 Rwork/Rfree (%) 16.2/19.5

No. of atoms

Protein 548 Water 126 B factors

Protein 7.5 Water 22.2 R.m.s.d.

Bond lengths () 0.011 Bond angles (o) 1.157

*Highest-resolution shell is shown in parenthesis.

NATURE COMMUNICATIONS | 5:3814 | DOI: 10.1038/ncomms4814 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4814

d

Tetramer

a

1.2

100 10

1 0 0.1

F14 W18

L21

Experimental

Tetramer

10

1

Intensity ()

0.8

1

Dimer

p(r)

0.6

0.1

F14 W18 L21

0.4

0.01

Experimental

Dimer ([afii9851]2 = 3.8)

Tetramer ([afii9851]2 = 1.2)

b

0.2

N C

HDM2NTD 0.001

0

0 0.1 0.2 0.3 0

20 40 60 S (1) r ()

L21

W18

e

SAXS envelope (GASBOR)p53LZ2HDM2NTD complex HDM2NTD HDM2NTD

p53LZ2

F14

HDM2NTD

p53LZ2 N

N

90

c

L21

L21

60

F14

W18

W18

F14

WT p53TAD-HDM2NTD p53LZ2-HDM2NTD

50

Figure 4 | Structural studies of p53LZ2 homodimer and its complex with HDM2NTD. (a) Overall structure of p53LZ2 homodimer at 1.2 (left). Three key hydrophobic residues, F14, W18 and L21, are indicated. Twofold axis of the coiled-coil is shown by dotted line. A 2Fo Fc map contoured at

1.0 s is superimposed with p53LZ2 (right). (b) HADDOCK-based quaternary structure of p53LZ2-HDM2NTD complex. HDM2NTD molecules (marine) bound to p53LZ2 homodimer (green and yellow) are represented as surface or cartoon. Essential hydrophobic residues, F14, W18 and L21, for recognizing

HDM2NTD are shown as stick. (c) Structural comparison in HDM2 binding between WT p53TAD and p53LZ2. Close-up view of the hydrophobic pockets of HDM2 (marine) for the amino acid triad, F, W and L of WT p53TAD (magenta) and p53LZ2 (yellow) are shown as surface. (d) Experimental SAXS scattering curve (black, s.e.) plotted versus the theoretical CRYSOL calculations for the dimer (red dots) and tetramer (blue dots), and normalized distribution function of p53LZ2HDM2NTD quaternary complex are shown on the left and right, respectively. An expanded plot of the low-angledata (0oSo0.1) is shown in the top right insert (left). (e) HADDOCK-based structure, represented as a cartoon, was tted into the SAXS envelope (blue)

with dimensions.

correlation between experimental and calculated scattering patterns (sqrt(w) 1.347 (p53LZ2HDM2NTD) and 1.554

(p53LZ2HDMXNTD)). Finally, superposition of the HADDOCK structures of p53LZ2HDM complexes onto the corresponding SAXS envelopes reveals almost perfect spatial arrangement of the two HDM2NTD or HDMXNTD wings relative to the central p53LZ2 body (w2 1.2 and w2free 2.8

for p53LZ2HDM2NTD complex, w2 1.5 and w2free 3.3

for p53LZ2HDMXNTD complex) further supporting the 1:2 p53LZ2HDM complexes in solution (Fig. 4e and Supplementary Fig. 6).

p53LZ targets endogenous HDM2 and HDMX. The protein transduction strategy using the positively charged cell-penetrating peptide transactivator of transcription (TAT) was applied to increase the cell permeability of designed p53LZ2, thereby overcoming the membrane barrier by endocytosis and/or direct penetration into cells38. To investigate the ability of p53LZ2 to target endogenous HDM2 and HDMX, we rst produced the cell-permeable TAT-p53LZ2 to high purity from E. coli and then conducted immunoprecipitation studies in MCF7 (human breast carcinoma) cells. After 6 h incubation of cells treated with 10 mM of each LZ molecule (His-tagged TAT-p53LZ2 or His-tagged TAT-GCN4LZ), the cell lysates were prepared and subjected to

Ni-NTA pull-down assay, followed by HDM2 and HDMX western blot analysis. The intracellular HDM2 and HDMX specically co-precipitated with TAT-p53LZ2 but not with its non-functional equivalent TAT-GCN4LZ (Fig. 5a). To further evaluate the cell-penetrating TAT-p53LZ2-mediated inhibition of p53 binding by both HDM proteins in cells, we treated MCF7 with 10 mM of each LZ (TAT-p53LZ2 or TAT-GCN4LZ) for 6 h and subsequently determined the amount of endogenous HDM2 or HDMX associated with p53 in cell lysates by immunoprecipitation and western blot analysis. Our data clearly conrmed that the TAT-p53LZ2 transferred across the cell membrane induced a substantial depletion of both HDM proteins bound to the intracellular p53 in cell lysates compared with those of TAT-GCN4LZ (Fig. 5b). Taken together, these data demonstrate that TAT-p53LZ2 targets both HDM2 and HDMX in situ, showing the expandability of our research scope to the cellular level.

Reactivation of the p53 pathway in cancer cells by p53LZ2. To examine the effect of p53LZ2 on growth and viability in cancer cells, we treated various exponentially growing cancer cells expressing the wild-type p53 (osteosarcoma SJSA-1, colorectal carcinoma HCT 116, breast carcinoma MCF7 and lung carcinoma A549) or mutant p53 (epidermoid carcinoma A-431 and

6 NATURE COMMUNICATIONS | 5:3814 | DOI: 10.1038/ncomms4814 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4814 ARTICLE

a b

TAT-GCN4LZ

TAT-p53LZ2

kDa

His-TAT-p53LZ2 His-TAT-GCN4LZ

HDM2 HDMX MG132

+ + +

p53

HDM2

p53

HDMX

HDM2

HDMX

p53

+ +

+ +

+ +

+ +

- 50

- 50

- 50

IP: HDM2

IP: HDMX

Lysates

kDa

kDa

- 80

- 80

80

60

50

40

30

80

60

50

40

30

- 60

- 60

Ni-NTA pull down Ni-NTA pull down

Input Input

Figure 5 | TAT-p53LZ2 inhibits direct interactions between p53 and HDM2 or HDMX in situ. (a) MCF7 cells were treated with His-tagged TAT-p53LZ2 or His-tagged TAT-GCN4LZ for 6 h. A total of 10 mM of MG132 was treated for 4 h prior to cell lysis. Cell lysates were subjected to

Ni-NTA pull-down followed by immunoblotting with anti-HDM2 or anti-HDMX antibody. (b) Cell lysates were prepared as in a, and then subjected to immunoprecipitation with anti-HDM2 and anti-HDMX antibody followed by immunoblotting with indicated antibodies.

a

b

120

Mutant p53

Cell viability (% of control)

100

SJSA-1 HCT 116 MCF7 A-431 MDA-MB-231

80

Nutlin-3

TAT

GCN4LZ

Control

TAT-

p53LZ2

Nutlin-3

TAT

GCN4LZ

Control

TAT-

p53LZ2

TAT

GCN4LZ

Nutlin-3

Control

TAT-

p53LZ2

Nutlin-3

TAT

GCN4LZ

Control

TAT-

p53LZ2

Nutlin-3

TAT

GCN4LZ

60

MCF7 SJSA-1 A549HCT 116 A-431 MDA-MB-231

40

5 10 20 10 20 M 0 5 10 20 10 20 M 0 5 10 20 10 10 M 0 5 10 20 10 20 M 0 5 10 20 10 20 M

kDa

p53

0

20

p21

HDM2

GAPDH

Control

TAT-

p53LZ2

0

HDMX

TAT-p53LZ2 (M)

c d

SJSA-1 HCT 116 A-431 MDA-MB-231

60

1,250 1,000

750 500 250

60

40

20

TAT-p53LZ2

Doxorubicin

50

mRNA (fold change)

HCT 116

MCF7

p53 p21 HDM2 MIC-1

50

20 80

80

40

SJSA-1 MCF7 A549 JEG-3

SJSA-1

HCT 116

40

+ + + + + + + +

kDa

30

p53

p-p53S15

p21

HDM2

GAPDH

0.1 1 10

50

50

20 80

40

20

10

0 5 10 20 20 10 (M)

5 10 20 20 10 5 10 20 20 10 5 10 20 20 10

TAT-GCN4LZ

Nutlin-3

TAT-p53LZ2

TAT-GCN4LZ

Nutlin-3

0

TAT-p53LZ2

TAT-GCN4LZ

Nutlin-3

TAT-p53LZ2

TAT-GCN4LZ

Nutlin-3

Figure 6 | TAT-p53LZ2 reactivates p53 signalling pathway in cancer cells expressing wild-type p53. (a) Cell viability of four wild-type p53, SJSA-1, HCT 116, MCF7 and A549, and two mutant p53, A-431 and MDA-MB-231, cancer cell lines was measured by the CellTiter-Glo (CTG) assay after 5-d incubation with TAT-p53LZ2. Data represent the means.d. of three independent experiments. (b) TAT-p53LZ2 stabilizes endogenous p53 and elevates protein levels of p53 targets, p21CIP1/WAF1 and HDM2. Three wild-type p53, SJSA-1, HCT 116 and MCF7, and two mutant p53, A-431 and MDA-MB-231,

cancer cells were incubated with TAT-p53LZ2, TAT-GCN4LZ or Nutlin-3 and then analysed by western blotting. (c) Exogenously introduced TAT-p53LZ2 does not induce genotoxic response in cancer cells. Cells were incubated with 10 mM of TAT-p53LZ2 or 1 mM of doxorubicin for 24 h and then the phosphorylated p53 at S15 were determined by western blotting. (d) TAT-p53LZ2 but not TAT-GCN4LZ induces p53 target genes in cancer cells bearing

WT p53 in a dose-dependent manner. mRNA levels of p53, p21CIP1/WAF1, HDM2 and MIC-1 in cancer cells were determined by the TaqMan RT-qPCR

assays. Data represent the means.d. of three independent experiments.

breast carcinoma MDA-MB-231) with the highly pure recombinant TAT-p53LZ2 produced from E. coli. A dose-dependent anti-proliferative and cytotoxic efcacy of TAT-p53LZ2 accessed by CTG (CellTiter-Glo) assay demonstrated that the cell-permeable TAT-p53LZ2 signicantly diminished the cancer cell viability in the wild-type but not mutant p53 cell lines (IC50,

SJSA-1 3.12 mM, HCT 116 6.13 mM, MCF7 2.20 mM and

A549 4.94 mM) (Fig. 6a). In addition, we carried out western

blot analysis on cancer cell lines incubated with TAT-p53LZ2, TAT-GCN4LZ or Nutlin-3 for 24 h, demonstrating that TAT-p53LZ2 effectively rescued the endogenous p53 and increased the protein expression level of its downstream transcription targets, p21CIP1/WAF1 and HDM2, in a dose-dependent manner (Fig. 6b). On the other hand, p53 phosphorylation, which

NATURE COMMUNICATIONS | 5:3814 | DOI: 10.1038/ncomms4814 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4814

represents the genotoxic stress and p53 activation in both carcinoma and non-carcinoma by a synthetic compound (for example, doxorubicin), was not detectable in the TAT-p53LZ2-treated cancer cell lines, suggesting the non-genotoxic mechanism of p53 reactivation by delivering exogenously introduced TAT-p53LZ2 into the cancer cells (Fig. 6c). Further studies of the real-time quantitative PCR (RTqPCR) of the wild-type (SJSA-1 and HCT 116) and mutant p53 (A-431 and MDA-MB-231) cancer cell lines revealed that the accumulation of endogenous p53 in cancer cells expressing wild-type p53 by exogenous TAT-p53LZ2 was due to decreased proteasome-mediated p53 degradation rather than increased level of p53 gene expression, while the dose-dependent mRNA induction of three p53 transcriptional targets, p21CIP1/WAF1, HDM2 and MIC-1 (macrophage inhibitory cytokine-1), was strongly correlated with elevated p53 levels that restored p53 pathway in cancer cells expressing wild-type p53 (Fig. 6d).

Functional analysis of HDM inhibition by p53LZ2. The consequence of elevated expression levels of the cyclin-dependent kinase (Cdk) inhibitor p21CIP1/WAF1 mediated by TAT-p53LZ2 is to arrest cell cycle of cancer cells in G1 and G2 phases. Treatment of 10 mM of TAT-GCN4LZ, TAT-p53LZ2 or Nutlin-3 for 24 h on various cancer cell lines (wild-type p53: SJSA-1, MCF7, A549 and HCT 116, mutant p53: A-431) demonstrated that EdU (5-ethynyl-20-deoxyuridine)-labelled cancer cells harbouring wild-type p53 exhibited statistically remarkable cell cycle arrest in G1 and G2M phases and depletion of S phase population in

comparison with other cell line expressing mutant p53 (Fig. 7a). Especially, we observed most effective accumulation in G1 and

G2M phases of HCT 116 upon TAT-p53LZ2 treatment (35.8% more in G1 and 18.7% more in G2M), relative to the other WT p53 cells. Further studies by the quantitative ow cytometry (FACS) analysis also indicated that TAT-p53LZ2 but not TAT-GCN4LZ induces G1 and G2M cell cycle arrest with a decreased fraction of cells in S phase, conrming that the cell-permeable TAT-p53LZ2 reactivates key machineries of p53 pathway (Supplementary Fig. 8).

Next, we investigated the apoptotic activity of TAT-p53LZ2 in cancer cells with wild-type or mutant p53. All cancer cells were treated with 10 mM of Nutlin-3, TAT-GCN4LZ or TAT-p53LZ2 for 48 h, and then analysed by 7-AAD (7-aminoactinomycin D) and Annexin V staining in FACS. After 48 h of incubation with TAT-p53LZ2, 31.4, 22.1 and 18.9% of the cell populations in SJSA-1, MCF7 and HCT 116, respectively, were induced to apoptotic cell death, whereas any cell lines expressing mutant p53 (A-431 and MDA-MB-231) were indistinguishable from the cells treated with TAT-GCN4LZ (Supplementary Fig. 9).

Finally, TAT-p53LZ2 was evaluated for its potential anti-tumour efcacy in the human xenograft model. We chose to perform the study with human colorectal carcinoma HCT 116 (p53 / ). Previous in vitro and/or in vivo studies of HCT 116

with either Nutlin-3 or pMI peptide treatment exhibited a relatively low sensitivity17,39, encouraging us to tackle this tumour model to assess in vivo potency of TAT-p53LZ2 homodimer, harbouring B2.5-fold higher afnity against

HDMX in vitro compared with WT p53TAD as well as a

50

40

30

a b

TAT-GCN4LZ

Nutlin-3

TAT-p53LZ2

Annexin V positive cells (% increase from control)

20

10

0

S: 14.0%

S: 0.4% S: 0.6%

TAT-GCN4LZ TAT-p53LZ2

Nutlin-3

SJSA-1

MCF7

A549 HCT 116 A-431

EdU staining

G1: 81.8% G2/M: 4.1% G1: 85.5% G2/M: 14.1% G2/M: 12.7%

G1: 86.7%

S: 43.1% S: 10.3% S: 10.0%

S: 0.4%

G1: 57.4%

G1: 40.7% G2/M: 32.3%

G2/M: 16.2%

G1: 53.0% G2/M: 14.3%

G1: 14.2% G2/M: 24.1%

G1: 38.0% G2/M: 17.6% G1: 38.1% G2/M: 17.6% G1: 34.1% G2/M: 15.9%

DNA content

G1: 33.0%

G2/M: 57.0%

S: 32.7%

S: 61.4%

S: 0.8%

G1: 72.1% G1: 74.9%

G2/M: 27.1% G2/M: 24.7%

SJSA-1 MCF7 HCT 116 A-431 MDA-MB-231

(Mutant p53)

S: 5.8% S: 7.2%

G1: 54.6% G1: 50.0%

G2/M: 39.8% G2/M: 42.8%

S: 44.4% S: 44.3% S: 49.9%

Figure 7 | TAT-p53LZ2 induces cell cycle arrest and apoptosis in cancer cells. (a) For cell cycle analysis, ve cancer cell lines were treated with TAT-p53LZ2, TAT-GCN4LZ or Nutlin-3 for 24 h. Cells were incubated with 10 mM of EdU followed by Alexa uor 488 labelling and nuclear DNA staining with Hoechst 33342. Cell cycle distribution was calculated from the image analysis and expressed as the percentage of the total population rounded to the full percentage value. (b) Cell lines expressing wild-type or mutant p53 were cultured with 10 mM of TAT-p53LZ2, TAT-GCN4LZ or Nutlin-3 and the percentage of apoptotic cells (s.d.) was determined by the Annexin V assay after 48 h.

8 NATURE COMMUNICATIONS | 5:3814 | DOI: 10.1038/ncomms4814 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4814 ARTICLE

superb in vitro binding to HDM2. An orally administered Nutlin-3 (200 mg per kg twice a day for 3 weeks) strongly suppressed tumour growth (Po0.005 versus control; two-sided

Students t-test), while its intravenous administration was largely ineffective compared with vehicle control (P 0.37 versus

control; two-sided Students t-test). In contrast, TAT-p53LZ2 treated intravenously (0.6 mM once a day for 3 weeks) showed notably inhibited tumour growth of HCT 116 xenograft, relative to the vehicle control (Fig. 8). Neither a signicant weight loss in the mice nor any abnormalities upon autopsy was observed at the end of the treatment.

DiscussionThe signicant role of the E3 ubiquitin ligase HDM2 and subsequently identied HDMX in regulating tumour suppressor p53 was rst revealed by the nding that both HDMs can directly interact with p53 and more signicantly HDM2 in turn mediates mono-ubiquitination of p53 in vitro and in vivo1,2,11. A major oncogenic aspect of HDM2 and HDMX is the ability to abolish the transcriptional activity of tumour suppressor p53 (ref. 40). Paradoxically, HDM2 itself is the product of the p53 inducible gene, in which it is tightly associated through the auto-regulatory feedback loop aimed at effectively downregulating p53 activity in the absence of stimulus such as hypoxia or DNA damage. Furthermore, the emerging picture in which HDMX plays an independent but also cooperative role with HDM2 in regulating p53 stability makes the regulatory mechanism governing the cellular level of p53 seem more complicated9. In addition, recent studies suggest that a dynamic control of both p53 and HDM2 by HAUSP-mediated de-ubiquitination also contribute towards regulating the p53-HDM2 pathway41. Therefore, in terms of therapeutic needs in the oncology eld, the inactivation of p53 pathway occurred in most human cancers, promotes an awareness of a clear value proposition that restoring wild-type p53 in cancer cells by targeting both HDM2 and HDMX with a potential therapeutic antagonist would be a highly promising strategy for treating cancers.

The goal of our studies presented here is to strive towards molecular minimization of p53 while preserving functional activity. Myriad structural and functional versatilities of cellular proteins in biology allow us the way to imitate a particular form and/or function of natural proteins. In our case, the concept of how to properly design structural and functional mimicry performing mostly analogous roles compared with p53, but showing advanced structural and functional integrities was necessary. Accordingly, we came up with biomimetic approach, which is a bio-inspired methodology of creating and constructing new articial biology systems. Our rational grafting approach has successfully elicited the molecular mimicry of the N-terminal domain of p53 (p53TAD) in the grafted p53LZ2, enforcing the disruption of p53-HDM2 and/or HDMX interaction in vitro and in vivo. The wild-type p53TAD exists with B11% a-helices in solution with the majority as random coil25. Several pioneering concepts suggested that conformation-stabilized polypeptides that capture their bioactive forms have benets in vitro and in vivo4245. The p53LZ2 homodimer is stabilized by the LZ motif, through which the close packing of hydrophobic leucine and valine residues between each helices is formed, increasing molecular a-helicity to over 90%. The designed p53LZ2 occupies key accessible positions on HDM2 for the inhibition of protein protein interaction (PPI). A comparison of in vitro binding characteristics of p53LZ2 shows improvement of its binding to both HDMs compared with WT p53TAD (quantitatively B3.4-and B2.5-fold higher afnities to HDM2NTD and HDMXNTD, respectively). Considering the relatively larger binding area of grafted p53LZ2 homodimer in comparison with small synthetic compounds, the binding potency of p53LZ2 may be enhanced further. Moreover, the newly created electrostatic interactions between both proteins, p53LZ2 K3, E6, E20, E22 and HDM2NTD E69, K70, K94, K51, respectively, demonstrated by our structural studies, further reinforce the proteinprotein interaction compared with that of WT p53TAD. Given these ndings based on our in vitro biochemical and biophysical studies in this early stage of discovery research, we then focused on how to augment the intracellular efcacy of grafted p53LZ2 to eventually overcome a tremendous hurdle of HDM2 and HDMX overexpression in cancer.

Numerous efforts have shown that TAT-mediated peptide or protein delivery into cells is capable of modulating various functional activities in vitro and in vivo, including an autophagy-inducing peptide, TAT-beclin-1, and a novel class of neuroprotectant, TAT-NR2B9c (refs 46,47). TAT-p53LZ2, composed of the HIV-1 TAT protein transduction domain (PTD) attached via (GS)2 linker to 30 amino acids of each helix of parallel homodimeric p53LZ, as an antagonist to both HDM2 and HDMX, successfully targets endogenous HDM2 and HDMX in cancer cells, preventing p53-HDM interaction, bypassing proteasome-mediated degradation of p53, and thus, leading to the intracellular p53 stabilization. A subsequent activation of the Cdk inhibitor p21CIP1/WAF1 by p53-signalling cascade in response to exogenous treatment with TAT-p53LZ2 relies on the accumulation of endogenous p53 in cancer cells expressing wild-type but not mutant p53. Unfortunately, the cell-permeable TAT-p53LZ2 is likely more destructive in interfering with p53-HDM2 interaction than with p53-HDMX in our cultured cancer cells, even though TAT-p53LZ2 itself is capable of enhanced binding to HDMX compared with WT p53TAD in vitro. We could not observe statistically signicant difference in efcacy between HDM2-overexpressing (for example, SJSA-1) and HDMX-overexpressing (for example, JEG3) p53 / cancer

cells upon either HDM2-specic Nutlin-3 or TAT-p53LZ2 treatment. Only the anti-proliferative assay of MCF7 upon TAT-p53LZ2 treatment exhibited that the cell viability was

1,000

800

Mean tumour volume (mm3 )

400

200

0 0 4 8 12 Days post-treatment

Vehicle

TAT-p53LZ2 (IV)

Nutlin-3 (IV)

Nutlin-3 (PO)

600

*

**

16 19

Figure 8 | TAT-p53LZ2 inhibits tumour growth in a HCT 116 xenograft model. Nude mice (n 8 per group) were implanted subcutaneously (s.c.)

with HCT 116 and groups of mice (N 8) received vehicle, TAT-p53LZ2

(0.6 mM) or Nutlin-3 (10 mg kg 1) intravenously once a day for 3 weeks or Nutlin-3 (200 mg kg 1) orally twice a day for 3 weeks. Tumour size was assessed by direct caliper measurement. Data are represented as mean tumour volume s.e.m. (*Po0.05, **Po0.005 versus vehicle in two-sided Students t-test).

NATURE COMMUNICATIONS | 5:3814 | DOI: 10.1038/ncomms4814 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4814

reduced by 410% compared with that of Nutlin-3 (IC50 3.74 mM for Nutlin-3 versus 2.20 mM for TAT-p53LZ2)

(Supplementary Fig. 10). One of possible reasons for this is the elevated level of HDM2 in response to increased levels of wild-type p53 in cells by exogenously introduced TAT-p53LZ2, suggesting that a buildup of the intracellular HDM2 level may compel the equilibrium between both HDMs to move towards p53LZ2-mediated HDM2 homo-dimerization rather than HDMXHDMX and HDM2HDMX dimerization, thus decreasing the effective concentration of TAT-p53LZ2 available for HDMX targeting.

Our nding that p53 reactivation by antagonizing p53-HDM interaction with exogenously treated TAT-p53LZ2 is capable of inducing not only a cell cycle arrest but also an apoptotic outcome in cancer cells is clearly demonstrated. First, TAT-p53LZ2 and Nutlin-3 but not TAT-GCN4LZ exclusively induced both responses upon pharmacological activation of p53. Second, in contrast to mutant p53, defective in p53-mediated transcriptional activation, cancer cells bearing wild-type p53 were susceptible to both marked changes to regulate cell fate only upon p53 activation by TAT-p53LZ2. Finally, intravenous administration of TAT-p53LZ2 to subcutaneous xenograft HCT 116 tumour model in vivo moderated substantial tumour expansion in mice, compared with the vehicle control, suggesting a possible therapeutic use of TAT-p53LZ2. Future ne-tuned studies will further investigate how strongly the cell-permeable p53LZ2 exhibits its in vivo efcacy.

Previously identied high afnity peptide inhibitors (for example, pMI and pDI) disrupting WT p53-HDM interactions are known to possess limited therapeutic scopes by the fact that most of them are less stable and poorly internalized into cells17,48. It is likely that the largely unstructured nature of those single-chain peptide inhibitors leads to a relative instability in vivo. In addition, these naked peptide antagonists in the absence of any particular cell-penetrating moiety that facilitates cellular uptake of fusion partner across the cell membrane is mostly inefcient in cell internalization. Furthermore, stapled p53 peptides, which have to be modied through in vitro chemical synthesis, still have been under development and evaluation processes for therapeutic use25,27,49. Our alternative approach to develop a new molecular entity of HDM dual antagonist, which is readily producible from bacterial expression system, demonstrates that TAT-mediated internalization of highly stabilized a-helical scaffold of p53LZ2 exhibits on-target efcacy in cells. In conclusion, modulation of p53 response with our designed p53LZ2, generated by a rational protein grafting strategy and protein engineering, shows a therapeutic possibility against various cancers and provides a framework to further advance the development of peptide-based drug discovery platform.

Methods

Protein preparation. All constructs were generated by a standard PCR method (Supplementary Table 2). All p53LZ constructs and GCN4LZ were subcloned into NdeI-XhoI of pET21b, BamHI-XhoI of pGEX-4T3 or NdeI-XhoI sites of pDB-HisGST. Full-length (FL) WT p53TAD (residues 172) was subcloned into BamHI-XhoI of pGEX-4T3. FL HDM2, HDM2NTD (residues 1125) and HDMXNTD (residues 1132) were subcloned into NdeI-XhoI sites of pDB-HisMBP, pET21b or pET28a. Each plasmid was transformed into E. coli BL21(DE3)CodonPlus-RIPL (Agilent Technologies). Protein expression was induced by addition of 0.5 mM of IPTG at OD600B0.8 and then cells were further grown at 18 C for overnight. For co-purication of p53LZ2-HDM complexes, GST-p53LZ2 mixed with either HDM2NTD or HDMXNTD at 1:2 molar ratio was subjected to GST afnity chromatography followed by anion-exchange chromatography (Mono Q) and size-exclusion chromatography (superdex-75). GST-, MBP- or His-tagged individual proteins were puried by an afnity chromatography, followed by additional purication steps through Mono-Q and gel ltration in the absence or presence of protease cleavage (thrombin or TEV). For all cell-based assays and in vivo studies, the protein storage buffer of TAT-p53LZ2 and TAT-GCN4LZ was exchanged with RPMI-1640 (Gibco).

In vitro GST pull-down assay. GST, GST-p53LZ1, GST-p53LZ2 or GST-p53LZ3 (3 nmol) was incubated with either His-HDM2NTD or His-HDMXNTD (1:2 molar ratio) for 1 h at 4 C. A total of 20 ml of glutathione resin was used to pull down GST-fused proteins mixed with HDM. The binding reaction was performed at 4 C with gentle agitation. The GST-bound resin was washed three times with 500 ml of PBS (pH 7.4) buffer. Proteins were eluted in SDS sample buffer and separated by 12% SDSPAGE.

Multi-angle light scattering analysis. The molecular mass of p53LZ2 bound to HDM2NTD or HDMXNTD was determined by the quantitative analysis using MALS. Each sample (5 mg ml 1, PBS pH 7.4) was injected into a WTC-030S5 column (Wyatt Technology) coupled to 18-angle light scattering detector (DAWN

HELEOS II) and a refractive index detector (Optilab T-rEX) (Wyatt Technology). All data were collected at 25 C (0.5 ml min 1 of ow-rate) and analysed using

ASTRA 6 software.

Enzyme-linked immunosorbent assay. ELISA microtiter plates (Nunc) were coated with 100 ml of MBP-FL HDM2, MBP-HDM2NTD or MBP-HDMXNTD (500 nM) supplemented with 50 mM carbonate/bicarbonate (pH 9.6) for 16 h at4 C. After washing the wells with PBS containing 0.05% Tween-20 (PBST buffer), ELISA plates were incubated in blocking buffer (PBS, 0.1% Tween-20 and 3% skim milk) for 2 h at room temperature and then increasing amounts of binding protein (GST or GST-tagged proteins) were added to the wells for 1 h at room temperature. The plates were washed three times with PBST buffer and subsequently incubated with horseradish peroxidase-conjugated anti-GST antibody in blocking buffer (1:2,000) for 1 h. After three times washing, the plates were developed with TMB peroxidase substrate and the absorbance was measured at 450 nm.

Isothermal titration calorimetry. HDM2NTD, HDMXNTD and p53LZ2 homodimer were extensively dialyzed against PBS (pH 7.4) before titration. The calorimetric titrations were carried out at 18 C on MicroCal Auto-ITC200 system with 20 consecutive injections of p53LZ2 (0.3 mM, 2 ml each), spaced 200 s apart, into the cell containing either HDM2NTD (0.03 mM) or HDMXNTD (0.025 mM). The raw data were processed with ORIGIN software. The following equations were used to calculate; KD 1Ka 1 (KD dissociation constant, Ka association con

stant), RTlnKa DG DH TDS (DG free energy change, DH enthalpy

change and DS entropy change)

Circular dichroism spectroscopy. CD spectra and thermal melt analyses of p53LZ2 and GCN4LZ were performed with 20 mM of each protein in PBS (pH 7.4) using a Jasco J-715 spectrometer. The wavelength dependence of [y] was monitored from 190 to 260 nm at 50 nm min 1 speed. All spectra, accumulated three times, were corrected for solvent contributions. The melting proles were obtained at 220 nm with a heating rate of 2 C per min.

Molecular dynamics simulations. MD simulations were performed using GROMACS 4.5.5 coupled with GROMOS-96 53A6 forceeld50. The initial structures in this study were generated by the least-square tting between modied GCN4LZ (PDB ID 2ZTA) and WT p53TAD bound to HDM2 (PDB ID 1YCR) using LSQKAB51. Each system was solvated with SPC water and neutralized at the salt concentration of 150 mM by adding appropriate sodium and chloride ions. The initial structures were relaxed through energy minimization followed by two-step equilibrium simulations to relax the solvent and ions while restraints were applied to protein complexes. These calculations were conducted in two phases: 100 ps NVT (N moles, V volume, T temperature) followed by 100 ps NPT

(N moles, P pressure, T temperature) ensemble dynamics (T 310 K). MD

simulations in production phase were performed in the NPT ensemble for 60 ns with a time step of 2 fs. The long-range electrostatic interactions were treated with the Particle Mesh Ewald (PME) technique while non-bond interactions were truncated with an 1.2-nm cutoff. The MM/PBSA approach was applied to binding energy calculations using GMXAPBS tool52, which is a Bash/Perl script to facilitate the interface between GROMACS and APBS53. The average of last 15 ns MD simulations (600 frames) was considered to calculate the binding energy and the entropic term was neglected.

Structural modelling. The initial crystallization trials using a high-throughput nanolitre dispensing robot, the Mosquito (TTP labtech), were performed with hanging drop vapour diffusion method at 18 C. The optimized p53LZ2 crystals were grown against crystallization buffers containing 0.1 M TrisHCl (pH 9.0) and 11% ethanol. The drops were set up using 1 ml of the protein solution (3 mg ml 1 in 20 mM Tris HCl (pH 7.4) and 100 mM NaCl) and 1 ml of the reservoir solution. X-ray diffraction data were collected at 170 C with a single crystal ash-frozen

in crystallization buffer supplemented with nal 30% glycerol. A native crystal data set was collected at the PAL (Pohang Accelerator Laboratory, South Korea) 5C SB II beamline. Integration, scaling and merging of the diffraction data were performed using the XDS54. X-ray structure was determined by molecular replacement method with PHASER55 using the structure of GCN4LZ (PDB ID

10 NATURE COMMUNICATIONS | 5:3814 | DOI: 10.1038/ncomms4814 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4814 ARTICLE

2ZTA) as a search model. Model rebuilding was carried out using COOT56, and then modied model was rened by a rigid body protocol followed by energy minimization and simulated-annealing methods in PHENIX57. The nal structure was assessed with PROCHECK58. Comparative protein structure modelling was performed with HADDOCK using crystal structures of p53LZ2 homodimer, HDM2NTD (PDB ID 1YCR) and HDMXNTD (PDB ID 3U15). Active residues were based on solvent accessibility, as well as the structural superposition between p53LZ2 homodimer and WT p53TAD-HDM2NTD (or HDMXNTD) complex.

Small-angle X-ray scattering. SAXS data collection for p53LZ2-HDM2NTD and p53LZ2-HDMXNTD complexes in PBS pH 7.4 was conducted at the PAL 4C SAXS II beamline (0.675 ). All SAXS measurements were done in triplicate exposures at 18 C and the scattering from the buffer alone were subtracted from the solution scattering data of each protein sample. Exposure times were typically 10 min for individual measurements. The controlled sample-to-detector distance (14 m) was applied to collect the complete data set of both complexes. All data were processed and analysed using the ATSAS software suite including PRIMUS, GNOM and CRYSOL5961. Ab initio SAXS envelopes for both p53LZ2-HDM complexes were calculated using GASBOR37.

Immunoprecipitation and pull-down. For immunoprecipitation, MCF7 cells were solubilized in lysis buffer containing protease inhibitor and phosphatase inhibitor (Roche). Cell lysates were incubated with indicated antibody at 4 C overnight. The antigenantibody complex was added to protein A/G resins (Thermo) and incubated for 2 h at 4 C. The resins were collected by centrifugation, and then washed three times with lysis buffer and twice with buffer A supplemented with 20 mM TrisHCl pH 7.4, 150 mM NaCl and 0.1 M EDTA. For Ni-NTA pull-down, the cell lysates were incubated with Ni-NTA beads for 2 h at 4 C and then washed three times with buffer B containing 20 mM TrisHCl pH 7.4, 150 mM NaCl,1 mM EDTA, 0.2% Tween-20 and 20 mM imidazole. Bound proteins were eluted by 0.2% SDS and separated by SDSPAGE, followed by immunoblotting with the appropriate antibodies.

Cell assays and western blots. All human cell lines were purchased from the American Type Culture Collection (ATCC). Cells were cultured in the indicated media supplemented with 10% fetal bovine serum (Gibco) and 1% Zell Shield (Minerva biolabs) at 37 C in 5% CO2. For cell viability assay, cells (1 103

2 103 cells) were seeded in 96-well plates and treated with TAT-p53LZ2 or

TAT-GCN4LZ at the indicated concentrations for 5 days. Cell viability assays were performed in triplicate and assessed with CellTiter-Glo (Promega) by using EnVision Multilabel Reader (PerkinElmer). For western blotting, cells were grown in six-well plates (4 105 cells per well) and lysed in 0.10.2 ml of lysis buffer

(cOmplete lysis kit, Roche). The following antibodies were used; p53 (Santa Cruz sc-126, 1:200 dilution), HDMX (Santa Cruz sc-374147, 1:100 dilution), GAPDH (Santa Cruz sc-25778, 1:200 dilution), p21 (Cell Signaling no. 2947, 1:500 dilution), p-p53 (Cell Signaling no. 9284, 1:1,000 dilution), HDM2 (Abcam ab-16895, 1:50 dilution), anti-mouse IgG-HRP (Cell Signaling no. 7076, 1:2,000 dilution) and anti-rabbit IgG-HRP (Cell Signaling no. 7074, 1:2,000 dilution). Uncropped versions of representative western blots used in the current study are shown in Supplementary Figs 11 and 12.

Quantitative RTPCR. Total RNA was extracted using RNeasy Mini Kit (Qiagen) from cells treated with vehicle, TAT-p53LZ2, TAT-GCN4LZ or Nutlin-3 (Sigma) at 37 C for 24 h. cDNA was prepared from total RNA using the high-capacity RNA-to-cDNA kit (Applied Biosystems). RT-qPCR was performed by using TaqMan Gene Expression Master Mix (Applied Biosystems) and gene-specic primer/probe set (Applied Biosystems, TP53 (Hs00153349_m1), HDM2 (Hs00242813_m1), p21 (CDKN1A) (Hs00355782_m1), MIC-1 (GDF15) (Hs00171132_m1) and GAPDH (Hs99999905_m1)) on a LightCycler 480 (Roche). Data were analysed with the LightCycler 480 software (Roche) determiningthe threshold cycle (Ct). Relative quantities of mRNA were calculated using comparative cycle threshold method after normalization to GAPDH.

Cell cycle arrest and apoptotic assays. For cell cycle analysis, SJSA-1, MCF7, A549, HCT 116 and A-431 cells were plated in a m-bottom 96-well plate (ibidi GmbH). After 24 h, all cells were incubated with 10 mM of TAT-p53LZ2, TAT-GCN4LZ or Nutlin-3 for 24 h. For EdU uptake, cells were incubated with 10 mM of EdU for 1.5 h. EdU labelling was performed using Click-reaction with

Alexa uor 488 in EdU staining kit (C10350, Invitrogen). Nuclear DNA was stained with Hoechst 33342. Cell imaging was conducted using Incell Analyzer 6000 (GE healthcare). Captured images were analysed using MTA (multi-target analysis) module of Incell communicator (image analysis software, GE healthcare) to obtain parameters and cell cycle population. Nuclear integrated intensity of Hoechst 33342 and EdU as parameters were used to plot cell population. For apoptotic cell measurement by Annexin V-PE and 7-AAD staining, SJSA-1, MCF7, HCT 116, A-431 and MDA-MB-231 cells were incubated with 10 mM of

TAT-p53LZ2, TAT-GCNLZ or Nutlin-3 (Sigma). After 48 h, cells were harvested and resuspended in 0.2 ml of binding buffer and then stained with Annexin

V-PE/7-AAD (BD bioscience) for 15 min at room temperature, followed by FACS analysis.

Growth inhibition of human HCT 116 xenograft model. All animals were dosed and monitored according to guidelines from the Institutional Animal Care and Use Committee (IACUC) at ChemOn Inc. Female BALB/c nude mice (68 weeks old) were purchased from Japan SLC Inc (Shizuoka, Japan). Mice were injected (right ank) subcutaneously with HCT 116 cells (3 105 cells per mouse; suspended in

100 ml of PBS mixed with 100 ml of BD Matrigel). Once tumours reached a mean volume of B150 mm3, all tumour-bearing animals were randomized into four study groups according to tumour sizes. Mice were treated via tail vein with vehicle control (RPMI-1640) or TAT-p53LZ2. Nutin-3 was treated via either intravenously (dissolved in 5% DMSO and 5% dextrose) or orally (dissolved in 2% HPMC and0.2% Tween-80). Animals were monitored and tumours were measured every 23 days. The length (L) and width (W) of the tumour mass were measured, and the tumour volume (TV) was calculated as: TV (L W2) 2 1.

References

1. Picksley, S. M. & Lane, D. P. The p53-mdm2 autoregulatory feedback loop: a paradigm for the regulation of growth control by p53? Bioessays 15, 689690 (1993).

2. Honda, R., Tanaka, H. & Yasuda, H. Oncoprotein MDM2 is a ubiquitin ligase E3 for tumor suppressor p53. FEBS Lett. 420, 2527 (1997).

3. Leng, R. P. et al. Pirh2, a p53-induced ubiquitin-protein ligase, promotes p53 degradation. Cell 112, 779791 (2003).

4. Dornan, D. et al. The ubiquitin ligase COP1 is a critical negative regulator of p53. Nature 429, 8692 (2004).

5. Chen, D. et al. ARF-BP1/Mule is a critical mediator of the ARF tumor suppressor. Cell 121, 10711083 (2005).

6. Vassilev, L. T. MDM2 inhibitors for cancer therapy. Trends Mol. Med. 13, 2331 (2006).

7. Poyurovsky, M. V. et al. The Mdm2 RING domain C-terminus is required for supramolecular assembly and ubiquitin ligase activity. EMBO J. 26, 90101 (2007).

8. Linke, K. et al. Structure of the MDM2/MDMX RING domain heterodimer reveals dimerization is required for their ubiquitylation in trans. Cell Death Differ. 15, 841848 (2008).

9. Wade, M., Wang, Y. V. & Wahl, G. M. The p53 orchestra: Mdm2 and Mdmx set the tone. Trends Cell Biol. 20, 299309 (2010).

10. Kawai, H., Lopez-Pajares, V., Kim, M. M., Wiederschain, D. & Yuan, Z. M. RING domain-mediated interaction is a requirement for MDM2s E3 ligase activity. Cancer Res. 67, 60266030 (2007).

11. Brooks, C. L. & Gu, W. p53 ubiquitination: Mdm2 and beyond. Mol. Cell 21, 307315 (2006).

12. Okada, H. & Mak, T. W. Pathways of apoptotic and non-apoptotic death in tumour cells. Nat. Rev. Cancer 4, 592603 (2004).

13. Gomes, N. P. & Espinosa, J. M. Disparate chromatin landscapes and kinetics of inactivation impact differential regulation of p53 target genes. Cell Cycle 9, 34283437 (2010).

14. Grasberger, B. L. et al. Discovery and cocrystal structure of benzodiazepinedione HDM2 antagonists that activate p53 in cells. J. Med. Chem. 24, 909912 (2005).

15. Kritzer, J. A., Lear, J. D., Hodsdon, M. E. & Schepartz, A. Helical beta-peptide inhibitors of the p53-hDM2 interaction. J. Am. Chem. Soc. 126, 94689469 (2004).

16. Garca-Echeverra, C., Chne, P., Blommers, M. J. & Furet, P. Discovery of potent antagonists of the interaction between human double minute 2 and tumor suppressor p53. J. Med. Chem. 43, 32053208 (2000).

17. Pazgier, M. et al. Structural basis for high-afnity peptide inhibition of p53 interactions with MDM2 and MDMX. Proc. Natl Acad. Sci. USA 106, 46654670 (2009).

18. Vassilev, L. T. et al. In vivo activation of the p53 pathway by small-molecule antagonists of MDM2. Science 303, 844848 (2004).

19. Bttger, A. et al. Design of a synthetic Mdm2-binding mini protein that activates the p53 response in vivo. Curr. Biol. 7, 860869 (1997).

20. Kritzer, J. A. et al. Miniature protein inhibitors of the p53-hDM2 interaction. Chem. Bio. Chem. 7, 2931 (2006).

21. Li, C. et al. Turning a scorpion toxin into an antitumor miniprotein. J. Am. Chem. Soc. 130, 1354613548 (2008).

22. Li, C., Pazgier, M., Liu, M., Lu, W.-Y. & Lu, W. Apamin as a template for structural-based rational design of potent peptide activators of p53. Angew. Chem. Int. Ed. 48, 87128715 (2009).

23. Ji, Y. et al. In vivo activation of the p53 tumor suppressor pathway by an engineered cycotide. J. Am. Chem. Soc. 135, 1162311633 (2013).

24. Lee, H. et al. Local structural elements in the mostly unstructured transcriptional activation domain of human p53. J. Biol. Chem. 275, 2942629432 (2000).

NATURE COMMUNICATIONS | 5:3814 | DOI: 10.1038/ncomms4814 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 11

& 2014 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms4814

25. Bernal, F., Tyler, A. F., Korsmeyer, S. J., Walensky, L. D. & Verdine, G. L. Reactivation of the p53 tumor suppressor pathway by a stapled p53 peptide.J. Am. Chem. Soc. 129, 24562457 (2007).26. Baek, S. et al. Structure of the stapled p53 peptide bound to Mdm2. J. Am. Chem. Soc. 134, 103106 (2011).

27. Bernal, F. et al. A stapled p53 helix overcomes HDMX-mediated suppression of p53. Cancer Cell 18, 411422 (2010).

28. Gonzalez, Jr L., Woolfson, D. N. & Alber, T. Buried polar residues and structural specicity in the GCN4 leucine zipper. Nat. Struct. Biol. 3, 10111018 (1996).

29. OShea, E. K., Klemm, J. D., Kim, P. S. & Alber, T. X-ray structure of the GCN4 leucine zipper, a two-stranded, parallel coiled coil. Science 254, 539544 (1991).

30. Kussie, P. H. et al. Structure of the MDM2 oncoprotein bound to the p53 tumor suppressor transactivation domain. Science 274, 948953 (1996).

31. Popowicz, G. M., Czarna, A. & Holak, T. A. Structure of the human Mdmx protein bound to the p53 tumor suppressor transactivation domain. Cell Cycle 7, 24412443 (2008).

32. Pronk, S. et al. GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit. Bioinformatics 29, 845854 (2013).

33. de Vries, S. J., van Dijk, M. & Bonvin, A. M. J. J. The HADDOCK web server for data-driven biomolecular docking. Nat. Protoc. 5, 883897 (2010).

34. Hura, G. L. et al. Robust, high-throughput solution structural analyses by small angle X-ray scattering (SAXS). Nat. Methods 6, 606612 (2009).

35. Rambo, R. P. & Tainer, J. A. Super-resolution in solution X-ray scattering and its applications to structural systems biology. Annu. Rev. Biophys. 42, 415441 (2013).

36. Rambo, R. P. & Tainer, J. A. Bridging the solution divide: comprehensive structural analyses of dynamic RNA, DNA, and protein assemblies by small-angle X-ray scattering. Curr. Opin. Struct. Biol. 20, 128137 (2010).

37. Svergun, D. I., Petoukhov, M. V. & Koch, M. H. Determination of domain structure of proteins from X-ray solution scattering. Biophys. J. 80, 29462953 (2001).

38. Ho, A., Schwarze, S. R., Mermelstein, S. J., Waksman, G. & Dowdy, S. F. Synthetic protein transduction domains: enhanced transduction potential in vitro and in vivo. Cancer Res. 61, 474477 (2001).

39. Kitagawa, M., Aonuma, M., Lee, S. H., Fukutake, S. & McCormick, F. E2F-1 transcriptional activity is a critical determinant of Mdm2 antagonist-induced apoptosis in human tumor cell lines. Oncogene 27, 53035314 (2008).

40. Kubbutat, M. H., Jones, S. N. & Vousden, K. H. Regulation of p53 stability by Mdm2. Nature 387, 299303 (1997).

41. Li, M., Brooks, C. L., Kon, N. & Gu, W. A dynamic role of HAUSP in the p53-Mdm2 pathway. Mol. Cell 13, 879886 (2004).

42. Craik, D. J. Circling the enemy: cyclic proteins in plant defence. Trends Plant Sci. 14, 328335 (2009).

43. Craik, D. J. Chemistry. Seamless proteins tie up their loose ends. Science 311, 15631564 (2006).

44. Schafmeister, C., Po, J. & Verdine, G. L. An all-hydrocarbon cross-linking system for enhancing the helicity and metabolic stability of peptide. J. Am. Chem. Soc. 122, 58915892 (2000).

45. Walensky, L. D. et al. Activation of apoptosis in vivo by a hydrocarbon-stapled BH3 helix. Science 305, 14661470 (2004).

46. Aarts, M. et al. Treatment of ischemic brain damage by perturbing NMDA receptor-PSD-95 protein interactions. Science 298, 846850 (2002).

47. Shoji-Kawata, S. et al. Identication of a candidate therapeutic autophagy-inducing peptide. Nature 494, 201206 (2013).

48. Hu, B., Gilkes, D. M. & Chen, J. Efcient p53 activation and apoptosis by simultaneous disruption of binding to MDM2 and MDMX. Cancer Res. 67, 88108817 (2007).

49. Chang, Y. S. et al. Stapled a-helical peptide drug development: a potent dual inhibitor of MDM2 and MDMX for p53-dependent cancer therapy. Proc. Natl Acad. Sci. USA 110, 34453454 (2013).

50. Oostenbrink, C., Villa, A., Mark, A. E. & van Gunsteren, W. F. A biomolecular force eld based on the free enthalpy of hydration and solvation: the GROMOS force-eld parameter sets 53A5 and 53A6. J. Comput. Chem. 25, 16561676 (2004).

51. Kabsch, W. A solution for the best rotation to relate sets of vectors. Acta Crystallogr. A 32, 922923 (1976).

52. Spiliotopoulos, D., Spitaleri, A. & Musco, G. Exploring PHD ngers and H3K4me0 interactions with molecular dynamics simulations and binding free energy calculations: AIRE-PHD1, a comparative study. PLoS ONE 7, e46902 (2012).

53. Baker, N. A., Sept, D., Joseph, S., Holst, M. J. & McCammon, J. A. Electrostatics of nanosystems: application to microtubules and the ribosome. Proc. Natl Acad. Sci. USA 98, 1003710041 (2001).

54. Kabsch, W. Xds. Acta Crystallogr. D Biol. Crystallogr. 66, 125132 (2010).55. McCoy, A. et al. Phaser crystallographic software. J. Appl. Crystallogr. 40,

658674 (2007).

56. Emsley, P. & Cowtan, K. Coot: model-building tools for molecular graphics. Acta Crystallogr. D Biol. Crystallogr. 60, 21262132 (2004).

57. Adams, P. et al. PHENIX: building new software for automated crystallographic structure determination. Acta Crystallogr. D 58, 19481954 (2002).

58. Laskowski, R. A., MacArthur, M. W., Moss, D. S. & Thornton, J. M. PROCHECKa program to check the stereochemical quality of protein structures. J. App. Cryst. 26, 283291 (1993).

59. Konarev, P. V., Volkov, V. V., Sokolova, A. V., Koch, M. H. J. & Svergun, D. I. PRIMUS: a Windows PC-based system for small-angle scattering data analysis.J. Appl. Cryst. 36, 12771282 (2003).60. Svergun, D. I. Determination of the regularization parameter in indirect-transform methods using perceptual criteria. J. Appl. Cryst. 25, 495503 (1992).

61. Svergun, D. I., Barberato, C. & Koch, M. H. J. CRYSOLa program to evaluate X-ray solution scattering of biological macromolecules from atomic coordinates. J. Appl. Cryst 28, 768773 (1995).

Acknowledgements

We thank Kyoung-Sik Jin at PAL and Jung-Hye Ha at KBSI for excellent technical assistance. We thank Yunje Cho at POSTECH for sharing beam time and Dimitrios Spiliotopoulos for providing us with GMXAPBS tool. This work was supported by funds from Samsung Electronics Co. Ltd.

Author contributions

J.-H.L. conceived the study and designed most experiments. J.-H.L. performed protein purication, in vitro binding (ITC) and structural studies. E.K., J.K. and J.H. carried out cloning, cell culture, protein purication, in vitro binding assay (GST pull-down and ELISA), CD analysis, immunoprecipitation and western blot experiments. Y.O. conducted and supervised immunoprecipitation and western blot studies. J.L., K.H.L., H.Y.K. and S.A. performed cell viability assay, western blot, cell cycle analysis, and apoptosis experiment. S.A. designed and organized human xenograft experiment. D.S. analysed molecular dynamics simulation. K.H. provided bioinformatics support. S.Y.C., P.H.S. and Y.-I.K. discussed the study. J.C.P. provided comments. J.I.L. participated in project design. All authors analysed the data and contributed the manuscript preparation. J.-H.L. wrote the manuscript, which was reviewed and edited by the other coauthors.

Additional information

Accession codes: The atomic coordinate of p53LZ2 has been deposited in the Protein Data Bank under the accession code 4OWI.

Supplementary Information accompanies this paper at http://www.nature.com/nature

Web End =http://www.nature.com/ http://www.nature.com/nature

Web End =nature communications

Competing nancial interests: All authors are employees of Samsung Electronics Co. Ltd.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions

Web End =reprintsandpermissions /

How to cite this article: Lee, J.-H. et al. Protein grafting of p53TAD onto a leucine zipper scaffold generates a potent HDM dual inhibitor. Nat. Commun. 5:3814doi: 10.1038/ncomms4814 (2014).

12 NATURE COMMUNICATIONS | 5:3814 | DOI: 10.1038/ncomms4814 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2014 Macmillan Publishers Limited. All rights reserved.