Helix12-Stabilization Antagonist of Vitamin D Receptor - Bioconjugate

Jun 13, 2016 - (3) In order to recruit the coactivator, H12 of VDR adopts a stable active conformation by binding the VDR agonist. Various ..... Quali...
2 downloads 0 Views 2MB Size
Subscriber access provided by UNIV OF NEBRASKA - LINCOLN

Article

Helix12-stabilization antagonist of vitamin D receptor Akira Kato, Toshimasa Itoh, Yasuaki Anami, Daichi Egawa, and Keiko Yamamoto Bioconjugate Chem., Just Accepted Manuscript • DOI: 10.1021/acs.bioconjchem.6b00246 • Publication Date (Web): 13 Jun 2016 Downloaded from http://pubs.acs.org on June 14, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Bioconjugate Chemistry is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Bioconjugate Chemistry

Helix12-stabilization antagonist of vitamin D receptor

Akira Kato,1 Toshimasa Itoh,1 Yasuaki Anami,1 Daichi Egawa,1 Keiko Yamamoto1,*

1

Laboratory of Drug Design and Medicinal Chemistry, Showa Pharmaceutical University, 3-3165

Higashi-Tamagawagakuen, Machida, Tokyo 194-8543, Japan

*Corresponding author

Keiko Yamamoto: [email protected]; Phone, +81 42 721 1580; Fax, +81 42 721 1580

1

ACS Paragon Plus Environment

Bioconjugate Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Abstract To develop strong vitamin D receptor (VDR) antagonists and reveal their antagonistic mechanism, we designed and

synthesized vitamin D analogues with bulky side chains based on the “active antagonist” concept in which

antagonist prevents helix 12 (H12) folding. Of the synthesized analogues, compounds 3a and 3b showed strong antagonistic activity. Dynamic hydrogen/deuterium exchange coupled with mass spectrometry (HDX-MS) and

static X-ray crystal structure analyses indicated that compound 3a stabilizes H11-H12 but displaces H6-H7 so that

3a is a novel rather than “active” or “passive” type of antagonist. We classified 3a as a third type of antagonist and

called it “H11-H12 stabilization antagonist”. HDX-MS analysis indicated that antagonist 3b is an “active” antagonist. To date there are no reports relating to nuclear receptor antagonist that strongly stabilizes H12. In this

study, we found first VDR antagonist that stabilizes H12 and we showed that antagonistic mechanism is diverse

depending on each antagonist structure. Additionally, HDX-MS was proved to be very useful for investigations of

protein structure alterations resulting from ligand binding.

2

ACS Paragon Plus Environment

Page 2 of 35

Page 3 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Bioconjugate Chemistry

Introduction

1α,25-Dihydroxyvitamin D3 (1,25D3) is the active form of vitamin D3 and plays important roles in calcium metabolism, immunomodulation, and cellular differentiation and proliferation.1 1,25D3 exerts these actions by modulating gene transcription through binding to the vitamin D receptor (VDR), a member of the nuclear receptor (NR) superfamily.2 Helix 12 (H12) of the NR ligand-binding domain (LBD) plays a significant role in the

formation of the active form of NRs. The conformation of VDR changes from the inactive to active form by

binding 1,25D3 and then forms a heterodimer with retinoid X receptor (RXR). The VDR/RXR heterodimer binds to vitamin D response elements (VDREs) of the target gene and recruits the coactivator, resulting in transcription of the target gene.3 In order to recruit the coactivator, H12 of VDR adopts a stable active conformation by binding the

VDR agonist.

Various VDR ligands have been developed by many groups, and several vitamin D analogues are used in clinical applications, such as the treatment of metabolic bone diseases and skin diseases such as psoriasis.4,5 All of

the clinically used vitamin D analogues are VDR agonists. However, antagonists are needed to treat diseases of VDR hyperfunction, such as Paget’s disease of bone and osteopetrosis,6,7 but no VDR antagonists are currently

available for use in medicine.

In 2002, Greene et al. proposed that steroid hormone receptor antagonists could be categorized as either “active” or “passive”.8 Active antagonists have a bulky structure and inhibit the agonistic behavior of NRs in any

situation because a bulky substituent actively prevents H12 folding, which is a significant aspect of NR activation. 4-Hydroxytamoxifen,9 a selective estrogen receptor modulator, is a representative active antagonist. Indeed, X-ray 3

ACS Paragon Plus Environment

Bioconjugate Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

crystal structure analyses showed that 4-hydroxytamoxifen blocks H12 folding in estrogen receptor α (ERα).10 By contrast, a passive antagonist is comparatively small and fits into the ligand-binding pocket (LBP), permitting H12

folding but with an antagonistic effect. As the crystal structure of the complex of ERβ and its antagonist, 5,11-cis-diethyl-5,6,11,12-tetrahydrochrysene-2,8-diol (THC), adopts an agonist-binding form, THC is classified as a passive antagonist.8 Bulky side chain–containing VDR antagonists such as carboxylic ester (ZK series) compounds,11,12 adamantane compounds,13 and 1α,25-dihydroxyvitamin D3-26,23-lactams,14 are classified as active antagonists, and their antagonistic activity is believed to be caused by prevention of H12 folding by the bulky side chain.

However, this hypothesis has not been proven in structural biological studies. VDR analogues lacking a bulky side chain, such as TEI964715 and its derivatives15,16 and 22S-butyl vitamin D analogues17,18 (e.g., compound 1) (Figure

1), are classified as passive antagonists. The X-ray crystal structure showed that antagonist 1 fits into the LBP and permits H12 folding.17 In some situations, a passive antagonist may exhibit agonist activity due to the absence of

steric repulsion between H12 and the ligand.

Hydrogen/deuterium exchange coupled with mass spectrometry (HDX-MS) is often used to investigate protein-ligand interactions, ligand-binding modes, protein-protein interactions, and protein folding.19–22 HDX-MS

enables the analysis of protein dynamics in solution, unlike static X-ray crystal structure analysis. HDX-MS

utilizing deuterium labeling can be used to investigate the rate at which protein backbone amide hydrogens undergo exchange.19 The rate of hydrogen/deuterium exchange of backbone amide hydrogens reflects the conformational mobility, hydrogen bonding strength, and solvent accessibility in the protein structure.20 In fact, the in-solution 4

ACS Paragon Plus Environment

Page 4 of 35

Page 5 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Bioconjugate Chemistry

dynamics of several NRs, including VDR, have been investigated using HDX-MS.23–33

To develop strong antagonists and to reveal the mechanism of active antagonism of the VDR, we

designed and synthesized vitamin D analogues with bulky substituents at the side chain terminus and used them to

probe changes in the VDR-LBD structure upon binding to synthetic active antagonists. Here, we report the

synthesis and biological evaluation of compounds 2a,b-6a,b (Figure 1) and the results of static X-ray crystallographic and dynamic HDX-MS analyses of the VDR-LBD/antagonist complex.

5

ACS Paragon Plus Environment

Bioconjugate Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 35

OH

H

HO

OH

1α,25-(OH)2D3 (1,25D3)

ZK168281

22

TEI9647

25

2a : R=H 2b : R=Bu

1

3a : R=H 3b : R=Bu

R OH

H

HO

4a : R=H 4b : R=Bu

OH

5a : R=H 5b : R=Bu

Figure 1. Structures of 1,25-(OH)2D3 (1,25D3) and its analogues.

6

ACS Paragon Plus Environment

6a : R=H 6b : R=Bu

Page 7 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Bioconjugate Chemistry

Results

Design and Synthesis Vitamin D analogues 2a,b-6a,b were designed with a bulky substituent at the 25 position (Figure 1). These bulky substituents were expected to prevent H12 folding due to steric repulsion. In addition, a 22S-butyl group was introduced in 2b-6b because it was found to be an effective VDR antagonist.18

The synthesis of target compounds 2a,b-6a,b is shown in Scheme 1. Alcohol 7, which was synthesized in our previous study,18 was tosylated with TsCl to give 8. Treatment of tosylate 8 with KCN afforded cyanate 9,

which was then reduced to aldehyde 10. Aldehyde 10 was converted to compound 11a using a mild iodine oxidation method developed by our group previously.34 Compound 11b was prepared according to a previously reported procedure.18

Key compounds 11a,b were treated with n-BuLi to give dibutyl compounds 12a,b, which were then

deprotected with 10-camphorsulfonic acid (CSA) to afford target compounds 2a and 2b, respectively. Compounds

11a,b were treated with BnMgCl and then CSA to give compounds 3a and 3b, respectively. Diphenyl compounds

4a,b, dinaphthylmethyl compounds 5a,b, and dinaphthyl compounds 6a,b were obtained by similar synthetic procedures using appropriate Grignard reagents.

Scheme 1. Synthesis of compounds 2a,b-6a,b.

7

ACS Paragon Plus Environment

Bioconjugate Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

X

X

H

TBSO

H KCN

I2, KOH

DMSO

MeOH

OTBS

TBSO

OTBS

COOMe

H

TBSO

OTBS 11a (94%)

9 : X = CN (92%) 10 : X = CHO (63%) R1

7 : X = OH 8 : X = OTs (78%) R1

Bu OH Bu

CO2Me

H

TBSO

Page 8 of 35

H n - BuLi

CSA

THF

MeOH

OTBS

TBSO

11a : R1 = H 11b : R1 = Bu

2a : R1 = H (96%) 2b : R1 = Bu (88%)

OTBS

12a : R1 = H (89%) 12b : R1 = Bu (75%) R1 Bn OH Bn H

11a 11b

BnMgCl

CSA

THF

MeOH TBSO

3a : R1 = H (84%) 3b : R1 = Bu (92%)

OTBS 13a : R1 = H (64%) 13b : R1 = Bu (64%) R1 R2 OH R2 H

11a 11b

R2MgBr

CSA

THF

MeOH TBSO

OTBS

4a : R1 = H, R2 = Ph (89%) 4b : R1 = Bu, R2 = Ph (79%) 5a : R1 = H, R2 = 2-Np-CH 2- (75%) 5b : R1 = Bu, R2 = 2-Np-CH 2- (94%) 6a : R1 = H, R2 = 2-Np (92%) 6b : R1 = Bu, R2 = 2-Np (96%)

14a : R1 = H, R2 = Ph (59%) 14b : R1 = Bu, R2 = Ph (81%) 15a : R1 = H, R2 = 2-Np-CH 2- (50%) 15b : R1 = Bu, R2 = 2-Np-CH 2- (57%) 16a : R1 = H, R2 = 2-Np (61%) 8 16b : R1 = Bu, R2 = 2-Np (65%)

ACS Paragon Plus Environment

Page 9 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Bioconjugate Chemistry

Biological Activities VDR binding affinity was evaluated using a competitive binding assay according to a previously reported procedure.35 This assay was performed using [3H]-1,25D3 and recombinant human VDR-LBD expressed as the N-terminal GST-tagged protein using pGEX-VDR vector in Escherichia coli BL21. The results are summarized in

Table S1. All synthetic compounds showed significant binding to the VDR (IC50: 0.14-130 nM), indicating that they are VDR ligands. The 25-dibenzyl analogues 3a (IC50: 0.43 nM), 3b (IC50: 0.14 nM) and 25-dinaphthylmethyl analogues 5a (IC50: 7.4 nM), 5b (IC50: 8.0 nM) showed stronger affinity than the 25-diphenyl analogues 4a (IC50: 24 nM), 4b (IC50: 12 nM) and 25-naphthyl analogues 6a (IC50: 130 nM), 6b (IC50: 34 nM), respectively. These results indicate that compounds with flexible substituents at C-25 show stronger affinity than those with rigid

substituents. The 22S-butyl compounds 2b-4b and 6b (except 5b) showed stronger affinity than the corresponding

22-H compounds 2a-4a and 6a, respectively. That is, the 22S-butyl group increased VDR affinity.

The ability of synthetic compounds 2a,b-6a,b to induce transcription of vitamin D–responsive genes was examined using the mouse osteopontin luciferase reporter gene assay system in Cos7 and HEK293 cells, according to a previously reported procedure.35 The results are shown in Figure 2. Transcriptional activities in HEK293 cells were stronger than in Cos7 cells, which is the same as previous studies.17,36 Except for compounds 2a, 5b, and 6b, most compounds did not show transactivation activity in Cos7 cells but instead exhibited concentration-dependent

inhibition of transactivation induced by 1,25D3. In HEK293 cells, only compounds 3a,b did not exhibit transactivation activity and inhibited the transactivation induced by 1,25D3. From these results, we concluded that 9

ACS Paragon Plus Environment

Bioconjugate Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

compounds 3a and 3b are potent VDR antagonists.

10

ACS Paragon Plus Environment

Page 10 of 35

Page 11 of 35

Cos7

3b

4b

5b

6b

2b

3b

-log [M]

2a

3a

4a

5a

6a

1,25D3 10-8 M

4b

10 9 8 7 6 5

0

-log [M]

10 9 8 7 6 5

0.5

10 9 8 7 6 5

10 9 8 7 6 5

10 9 8 7 6 5

12 11 10 9 8 7 6 5

12 11 10 9 8 7 6 5

0.5

1

10 9 8 7 6 5

1

1.5

EtOH 1,25D3

relative luciferase activity

1.5

0

2b

d

EtOH 1,25D3

relative luciferase activity

c

12 11 10 9 8 7 6 5

6a

12 11 10 9 8 7 6 5

5a

12 11 10 9 8 7 6 5

12 11 10 9 8 7 6 5

4a

12 11 10 9 8 7 6 5

12 11 10 9 8 7 6 5

3a

0

-log [M]

12 11 10 9 8 7 6 5

12 11 10 9 8 7 6 5

2a

0.5

12 11 10 9 8 7 6 5

12 11 10 9 8 7 6 5

0

12 11 10 9 8 7 6 5

0.5

1

EtOH 1,25D3

1

12 11 10 9 8 7 6 5

relative luciferase activity

b

EtOH 1,25D3

relative luciferase activity

a

-log [M]

5b

6b

1,25D3 10-8 M

HEK293

12 11 10 9 8 7 6 5

3b

4b

5b

6b

2b

3b

-log [M]

2a

12 11 10 9 8 7 6 5

12 11 10 9 8 7 6 5

12 11 10 9 8 7 6 5

EtOH 1,25D3

0.5

3a

1

0.5

0

-log [M]

5a

EtOH 1,25D3

1

0

2b

h relative luciferase activity

g

12 11 10 9 8 7 6 5

6a

12 11 10 9 8 7 6 5

5a

12 11 10 9 8 7 6 5

12 11 10 9 8 7 6 5

4a

12 11 10 9 8 7 6 5

12 11 10 9 8 7 6 5

3a

0

-log [M]

12 11 10 9 8 7 6 5

12 11 10 9 8 7 6 5

2a

0.5

12 11 10 9 8 7 6 5

12 11 10 9 8 7 6 5

0

12 11 10 9 8 7 6 5

0.5

1

EtOH 1,25D3

1

12 11 10 9 8 7 6 5

relative luciferase activity

f

EtOH 1,25D3

relative luciferase activity

e

relative luciferase activity

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Bioconjugate Chemistry

1,25D3 10-9 M

4b 1,25D3 10-9 M

Figure 2. Transactivation of synthetic compounds. Transactivation of compounds 2a,b-6a,b was evaluated using Cos7 (a-d) and HEK293 (e-h) cells. Transcriptional activity was evaluated using a dual luciferase assay with a full-length human VDR expression plasmid (pCMX-hVDR), a reporter plasmid containing three copies of mouse osteopontin VDRE (SPPx3-TK-Luc), and an internal control plasmid containing sea pansy luciferase expression constructs (pRL-CMV), as described previously.35 In experiments in Cos7 cells, luciferase activity equal to 10−8 M 1,25D3 was defined as 1, and in experiments in HEK293 cells, activity equal to 10−9 M 1,25D3 was defined as 1. Inhibitory effects on transactivation were evaluated in the presence of 1,25D3 (c, d, g, h).

11

ACS Paragon Plus Environment

-log [M]

Bioconjugate Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

X-ray Crystal Structure Analysis To understand the mechanism of antagonism of compounds 3a and 3b, we attempted to analyze the X-ray

crystal structure of rat VDR-LBD complexed with 3a or 3b. In the absence of coactivator peptide, we could not obtain crystals. As all crystal structures of rat VDR-LBD/ligand complex reported so far were ternary complex with

DRIP205, we also added DRIP205. In the presence of coactivator peptide DRIP205, co-crystals of 3a were

obtained but not co-crystals of 3b. X-ray diffraction images of the 3a/VDR-LBD complex were obtained at 2.0-Å resolution. The data collection and refinement statistics are summarized in Table S2.

Crystals of the 3a/VDR-LBD complex were an asymmetric homodimer (Figure 3), whereas all crystals of ligand/VDR-LBD complexes reported to date are monomers. We termed the monomers “chain A” (red) and “chain

B” (yellow), respectively (Figure 3a). Chain A could be analyzed to Ser114 based on the clear 2Fo − Fc map at the N-terminus, whereas chain B was analyzed to Lys123 at the N-terminus (Figure 3a). In chain A, we found a new

1.5-pitch α-helix at the N-terminus and designated it “helix 0” (H0). H0 formed significant interactions with chain B and neighboring chain A (chain A’), as shown in Figure 3b, c. H0 appears to have been detected due to

stabilization caused by these interactions. N-terminus helix similar to H0 was observed in a complex of 1,25D3/VDR-LBD (R270L) (PDB code 3VT3),37 but the interactions differed from those of the 3a/VDR-LBD complex. This difference was considered reasonable because the crystal packing of the complexes differ.

12

ACS Paragon Plus Environment

Page 12 of 35

Page 13 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Bioconjugate Chemistry

Figure 3. X-ray crystal structure of the 3a/VDR-LBD complex. (a) Overall view of the 3a/VDR-LBD asymmetric homodimer. Chains A and B are shown as red and yellow ribbons, respectively. Only chain A contained a novel 1.5-pitch α-helix, “helix 0” (H0), at the N terminus. (b,c) Interactions between H0 of chain A and surrounding amino acid residues (pale blue shows neighboring chain A, chain A’). Hydrogen bonding and hydrophobic interactions are shown with red and green dashed lines, respectively.

As shown in Figure 3a, except for the N-terminus, almost no differences were observed between chain A

and chain B including the 3a conformation. Indeed, the Cα RMSD value between chain A and chain B was 0.24 Å. We therefore discuss chain A as being representative for this study. Figure 4a shows chain A (3a/VDR-LBD), and 13

ACS Paragon Plus Environment

Bioconjugate Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 4b shows the superimposition of chain A on the 1,25D3/VDR-LBD complex (PDB code 2ZLC).38 There are two significant differences between chain A and the 1,25D3/VDR-LBD complex. First, H0 was observed in chain A but not in the 1,25D3/VDR-LBD complex (Figure 4b). Second, because of steric repulsion between the benzyl group of 3a and loop 6-7, the main chain between Asp292 of H6 and Leu305 of H7 was markedly shifted from the canonical positon, and consequently, one pitch of the H7 N-terminus unraveled (Figure 4c). The shifted H6/loop

6-7/H7 and the unraveled H7 N-terminus were clearly observed in the electron density map (Figure 4d).

The A and CD rings of ligand 3a adopted a similar position to those of 1,25D3 (Figure 4e). The hydroxyl groups of 3a formed hydrogen bonds with five amino acid residues (Ser233, Arg270, Tyr143, Ser274, and His393),

just like 1,25D3, but interestingly, the 25-OH group of 3a formed a hydrogen bond with Gln396 of H10 instead of H301 of loop 6-7 (Figure 4e, f). Because of steric repulsion between 3a and loop 6-7, His301 shifted dramatically

from its original position, and Gln396 was rotated at Cβ to occupy empty space (Figure 4e). Figure 4g shows the electron density map of 3a and the benzyl group of 3a fits the electron density map well. The H12 of chain A adopted an almost identical position to that of the 1,25D3/VDR-LBD complex (Figure 4b). This result indicates that H12 of 3a/VDR-LBD assumes an active conformation similar to that of agonist/VDR-LBD.

14

ACS Paragon Plus Environment

Page 14 of 35

Page 15 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Bioconjugate Chemistry

15

ACS Paragon Plus Environment

Bioconjugate Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 4. X-ray crystal structure of chain A of 3a/VDR-LBD. (a) Overall view of chain A. (b) Superposition of 1,25D3/VDR-LBD (green) and chain A (red). (c) Conformational change in H6/loop 6-7/H7 of chain A. Compared with 1,25D3/VDR-LBD, the main chain was shifted from a canonical positon, and one pitch of H7 unraveled because the benzyl group of 3a pushes against loop 6-7. (d) The electron density map of H6/loop 6-7/H7 in chain A. (e) His301 was shifted and Gln396 was flipped because the benzyl group of 3a pushes against loop 6-7. The 25-OH group forms a hydrogen bond with Gln396 on H10 instead of with His301. Yellow dashed line shows the hydrogen bond between 1,25D3 and VDR-LBD; orange dashed line shows the hydrogen bond between 3a and VDR-LBD. (f) Hydrogen bonds between VDR-LBD and 3a. (g) The electron density map with refined 3a bound to VDR-LBD.

HDX-MS Analysis To investigate protein dynamics in solution, we performed HDX-MS analyses in the presence and

absence of synthetic antagonists. Rat (rather than human) VDR-LBD was used to enable direct comparison with the

crystal structure of rat VDR-LBD. HDX-MS experiments were performed in D2O buffer (10 mM Tris-HCl [pH 8.0], 200 mM NaCl, 2mM DTT) at six time points (10, 30, 60, 300, 900, and 3600 sec). Deuterated VDR-LBD was

digested with pepsin, and the resulting peptides were analyzed by HPLC-MS/MS. In these experiments, duplicates

were analyzed for each time point. The rate of hydrogen/deuterium exchange was calculated using the following

equation: deuterium level (%) = [(mt – m0%)/(m100% – m0%)] × 100, where mt, m0%, and m100% represent the centroid value of the partially deuterated, nondeuterated, and fully deuterated peptides, respectively.21 A total of 23 peptic

peptides covering 95% of the VDR-LBD sequence (residues 125-423) were analyzed.

Apo VDR-LBD First, we performed HDX-MS analysis of apo VDR-LBD (Figure 5 and Figure S1) to characterize the

protein’s stability. In addition, the resulting data enabled us to calculate the average difference in deuterium level 16

ACS Paragon Plus Environment

Page 16 of 35

Page 17 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Bioconjugate Chemistry

between VDR-LBD in the presence of ligand and apoVDR-LBD (differential HDX analysis). Results of the

HDX-MS analysis of apoVDR-LBD are shown in Figure 5. Peptic peptides are depicted according to color-coded

deuterium level for each exchange time point. The most unstable regions were H1 (residues 117-133), loop 2-3

(residues 151-215), H3 (residues 216-229), loop 11-12 (residues 408-415), and H12 (residues 416-423); >80% of

these regions were deuterated with a 10-sec exposure to D2O (Figure 5). The most stable region was a portion of the N-terminus of H10 (residues 375-381), which underwent only about 20% exchange to deuterium after 3600 sec (Figure 5). Similar results were observed in HDX-MS analyses of human VDR-LBD and some NR-LBDs.21,29

17

ACS Paragon Plus Environment

Bioconjugate Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Hinge

loop 1-2

H1

Page 18 of 35

loop 2-3

H2

106 116 120 130 140 150 160 215 GS HMGS P NS P L KD S L RP K L S E EQQH I I A I L L D A HHK T YD P T YAD F RD F RP P V RMDGS TGS V T L 10 sec 30 sec 60 sec 300 sec 900 sec 3600 sec 10 sec 30 sec 60 sec 300 sec 900 sec 3600 sec 10 sec 30 sec 60 sec 300 sec 900 sec 3600 sec

H3

H3’

H4

H5

S1

D (%) ≥ 90 < 90 < 80 < 70 < 60 < 50 < 40 < 30 < 20 < 10

216 220 230 240 250 260 270 280 D L S P L SM L P H L AD L V S YS I QK V I GF A KM I PGF RD L T SDDQ I V L L K S S A I E V I M L RS NQS F TMDD M

S2

S3

H6

loop 6-7

H7

loop 7-8

loop 8-9

H8

H9

281 290 300 310 320 330 340 347 SWDCGSQD YK YD V TD V S K AGH T L E L I E P L I K F QVG L K K L N L HE E E HV L L MA I C I V S PD R PGVQD A K L

H9

loop 9-10

H10

H11

loop 11-12

H12

348 360 370 380 390 400 410 420 423 V E A I QD R L S N T L Q T Y I RC RHP P PGS HQ L YA KM I QK L AD L RS L NE E HS KQ YR S L S F QP E N SMK L T P L V L E V F GN E I S

Figure 5. Amino acid sequences and deuterium levels of apoVDR-LBD peptic peptides. In this study, 23 VDR-LBD peptic peptides were analyzed, covering 95% of the sequence of VDR-LBD (residues 125-423). Horizontal blocks represent peptic peptides for which the deuterium levels (D) are shown in color code for the six different time points. Amino acid residues that form hydrogen bonds with 1,25D3 are identified by solid red circles.

VDR-LBD/Ligand Complex 18

ACS Paragon Plus Environment

Page 19 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Bioconjugate Chemistry

Next, we analyzed VDR-LBD by HDX-MS in the presence of ligand (1,25D3, 3a, and 3b) and calculated the deuterium levels (Figures S2 and S3). To investigate the ligand-induced conformational dynamics of VDR-LBD,

we then conducted differential HDX analyses of VDR-LBDs in the presence and absence of ligand. The average

differences in VDR-LBD deuterium levels are shown in Table 1. A positive value indicates an increase in the

deuterium level of a peptide region destabilized by ligand-binding. In contrast, a negative value represents a

decrease in deuterium exchange due to stabilization brought about by ligand-binding. The HDX profiles of

VDR-LBD in presence of 1,25D3, 3a, and 3b were then mapped onto the crystal structure of the 1,25D3/VDR-LBD complex (PDB code 2ZLC) (Figure 6).

19

ACS Paragon Plus Environment

Bioconjugate Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 35

Table 1. Change in HDX profile for each peptide in the presence of ligands. Region

z

Sequence

Structure

1,25D3

3a

3b

117-133 128-136 134-147 134-150 151-215 216-229 230-240 245-255 251-275 288-299 300-323 305-312 305-321 316-329 330-347 348-359 360-384 375-381 376-400 385-405 403-413 408-415 416-423

2 1 2 2 2 1 1 1 3 1 3 1 2 1 2 1 3 1 3 2 1 1 1

KDSLRPKLSEEQQHIIA QQHIIAILL ILLDAHHKTYDPTY ILLDAHHKTYDPTYADF RDFRPPVRMDGSTGSVTL DLSPLSMLPHLADL VSYSIQKVIGF PGFRDLTSDDQ TSDDQIVLLKSSAIEVIMLRSNQSF DYKYDVTDVSKA GHTLELIEPLIKFQVGLKKLNLHE LIEPLIKF LIEPLIKFQVGLKKLNL KKLNLHEEEHVLL MAICIVSPDRPGVQDAKL VEAIQDRLSNTL QTYIRCRHPPPGSHQLYAKMIQKLA LYAKMIQ YAKMIQKLADLRSLNEEHSKQYRSL DLRSLNEEHSKQYRSLSFQPE QPENSMKLTPL MKLTPLVL EVFGNEIS

Hinge/H1 H1 H1/Loop H1/Loop/H2 Loop Loop/H3 H3/H3’ H3’/H4 H4/H5/S1 S3/H6 Loop/H7/Loop H7 H7/Loop Loop/H8 H8/Loop/H9 H9 H9/Loop/H10 H10 H10/H11 H10/H11 Loop/H12 H12 H12

4 1 -4 -31 3 -52 -22 3 -33 -13 -22 -29 -42 -1 -1 -2 -22 6 -27 -20 -39 -25 -31

-2 1 -3 -29 -5 -56 -30 4 -27 5 13 8 -26 -12 -2 -1 -18 16 -25 -27 -34 -38 -57

-1 -2 4 -17 -4 -48 -24 2 -47 -2 7 1 -41 -14 4 0 -12 16 -23 7 6 0 -11

∆D (%) ≤ 20 ≤ 10 n.s. ≥ -10 ≥ -20 ≥ -30 ≥ -40 < -40

The values represent difference in percentage of deuterium level of apo VDR-LBD and VDR-LBD in the presence of ligand. Deuterium level is the average of six H/D exchange time experiments: 10, 30, 60, 300, 900 and 3600 sec. A positive value (red) indicates a region increased in deuterium exchange in that region of VDR-LBD in the presence of ligand. A negative value (blue) represents a decrease in exchange. z value is the charge state of the peptide ion.

20

ACS Paragon Plus Environment

Page 21 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Bioconjugate Chemistry

21

ACS Paragon Plus Environment

Bioconjugate Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 6. Differences in HDX profile for each ligand. (a) 1,25D3, (b) 3a, (c) 3b. The changes in HDX profile presented in Table 1 were mapped onto VDR-LBD (PDB: 2ZLC). Undetected regions are shown in black.

We would like to emphasize that, in the presence of 1,25D3, both sets of HDX data of ours and Grrifin's show good agreement, including the values of protection and deprotection, across the entire LBD.21 In the presence

of 1,25D3, deuterium exchange decreased in peptide regions constituting the VDR-LBP (H1-2, H3, H5, H6-7, and H10-12) (Table 1 and Figure 6a). These regions include amino acid residues that form hydrogen bonds with 1,25D3 [Tyr143 (H1), Ser233 (H3), Arg270 (H5), Ser274 (H5), His301 (loop 6-7), and His393 (H10)] and hydrophobic

residues (H12) that form hydrophobic interactions. These data indicate that the above-mentioned regions are

stabilized by 1,25D3 binding. The region encompassing residues 375-381 (H10), which is considerably resistant to exchange in apo VDR-LBD, nevertheless exhibited a slight increase in exchange upon 1,25D3 binding (Table 1). A similar observation was reported for human VDR-LBD.21 As this region is a part of the dimerization interface with

RXR, an increase in flexibility may be advantageous for heterodimerization with RXR.

Compound 3a caused a decrease in deuterium exchange similar to that associated with 1,25D3, except for H6-H7 (Table 1 and Figure 6b). H1-2, H3, H5, and H10-H12 were stabilized by hydrogen bonding and

hydrophobic interaction between 3a and VDR-LBD (Figure 4a, 4f). Interestingly, we found that the deuterium

levels in H6, loop 6-7, and H7 increased markedly in the presence of ligand 3a (Table 1), indicating that the

flexibility of these regions increases upon ligand 3a binding. This increased flexibility is consistent with the results

of analyses of the crystal structure of 3a/VDR-LBD, which showed that this region shifted substantially from the 22

ACS Paragon Plus Environment

Page 22 of 35

Page 23 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Bioconjugate Chemistry

canonical position because the benzyl group of 3a makes contact with loop 6-7 in the canonical position (Figure 4c).

Surprisingly, 3a decreased deuterium levels of H12 and its vicinity (-38% and -57%) and the decrease rate was

bigger than that by 1,25D3 binding, indicating that 3a strongly stabilizes H12 though it is VDR antagonist (Table 1 and Figure S2). Moreover, exchange in the region encompassing residues 375-381 (H10) was higher with 3a (16%) than 1,25D3 (6%). The HDX profile of compound 3b was similar to that of 3a except for H11, loop 11-12, and H12 (Table 1

and Figure 6c). In contrast to the results obtained with 3a, the regions encompassing residues 385-405 (H10-H11) and 403-413 (loop 11-12) exhibited a marked increase in deuterium level; that is, H11 and loop 11-12 became

unstable upon ligand 3b binding. In addition, the regions encompassing residues 408-415 (H12) and 416-423 (H12) exhibited an increase in deuterium level compared with 1,25D3 (Figures S2 and S3). Interestingly, the deuterium level of H12 was quite similar to that of apo VDR-LBD (Figure S3). These results indicate that ligand 3b does not contribute to the stabilization of H12 but instead destabilizes H11 and loop 11-12.

Docking Study Because we were unable to obtain crystals of the 3b/VDR-LBD complex, we conducted a docking study. As ligand 3b has a 22S-butyl group, VDR must form a butyl pocket.17,39,40 Therefore, we selected a VDR-LBD (PDB code 3WTQ), for which 22S-butyl-2-methylidene-19-nor-1α,25-dihydroxyvitamin D3 is a ligand.40 Crystals of this complex (3WTQ) exhibit the best resolution among crystals forming the butyl pocket. We performed a

normal docking study; however, 3b could not be docked in VDR-LBD. We hypothesized that 3b is too large to

dock in the LBP and therefore removed H12 (residues 412-420) and used this VDR-LBD (∆412-420) for the 23

ACS Paragon Plus Environment

Bioconjugate Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

docking study. The results of this analysis are shown in Figure S4a. Ligand 3b was then extracted from the data shown in Figure S4a and superimposed upon the original VDR-LBD, 3WTQ (Figure S4b), which revealed that the

22S-butyl group was accommodated in the butyl pocket (Figure S4c) while 25-benzyl group came into contact with

Val414 and Phe418 of H12 (Figure S4d). Repulsion between the benzyl group and H12 must therefore prevent the

proper folding of H12.

Recruitment of the Coactivator Peptide and RXRα α To elucidate the mechanism of antagonism of 3a and 3b, we examined ligand-dependent recruitment of a

coactivator peptide, SRC-1, and RXRα to the VDR using a mammalian two-hybrid assay in HEK293 cells according to a previously reported procedure.35 As shown in Figure S5, 3a exhibited moderate SRC-1 recruitment

and weak RXRα, whereas 3b did not recruit SRC-1 and exhibited only marginal recruitment of RXRα.

Discussion Synthetic compounds 2a,b-6a,b showed specific binding affinity for hVDR (Table S1). Bulky

compounds 5 and 6 exhibited less affinity for the VDR compared with 3 and 4, respectively, indicating that the

bulkiness of analogues is important for high VDR affinity. The 25-dibenzyl analogues 3a,b inhibited the activation

induced by 1,25D3 in Cos7 and HEK293 cells (Figure 2). The naphthyl analogues 5a,b and 6a,b, which had more bulkier substituents than 3a,b, showed agonistic activity, albeit weak, indicating that the appropriate bulkiness like

3a,b are critical for antagonistic activity. At the present time, we cannot explain why naphthyl analogues 5a,b and

6a,b show weak agonistic activity. 24

ACS Paragon Plus Environment

Page 24 of 35

Page 25 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Bioconjugate Chemistry

Extensive HDX-MS analyses of NRs have been conducted by Griffin et al. and other groups. Using the cutting-edge technology of HDX-MS, the behavior and mechanism of NRs such as RXRα,24,25 peroxisome proliferator-activated receptor γ (PPARγ),26-28 ERα,29,30 receptor-related orphan receptor γ (RORγ),31 liver receptor homolog-1 (LRH-1),32 and the progesterone receptor (PR)33 have been described. The VDR has also been studied,21-23 and those studies demonstrated that the stability of each region in apo hVDR-LBD differs and the

mobility of each region in hVDR-LBD or VDR/RXR heterodimer changes upon binding of agonists such as 1,25D3, ED-71, or alfacalcidol.21 Moreover, the stability of hVDR changes upon 1,25D3 binding, heterodimerization, VDRE binding, and coactivator recruitment.22

Using HDX-MS, we examined the effects of antagonists instead of agonists and found that three regions,

H6-H7, the N-terminus of H10, and H11/loop 11-12, become flexible upon ligand binding. 1,25D3 was shown to destabilize only the N-terminus of H10, similar to the results observed in NR-LBDs including hVDR-LBD.21,29

Antagonist 3a destabilized two regions, H6-H7 and the N-terminus of H10, and antagonist 3b destabilized three regions, loop 6-7/H7, H11/loop 11-12, and the N-terminus of H10 (Table 1 and Figure 6).

Antagonist 3a caused a dynamic conformational change in the main chain of H6-H7, as shown in the crystal structure (Figure 4c), and destabilized the same region, as shown by differential HDX analysis (Figure 6b).

Unexpectedly, H11-H12 of 3a/VDR-LBD adopted an almost canonical position (Figure 4b) and was stabilized more than that of 1,25D3/VDR-LBD (Table 1, Figure 6b). Thus, both crystal structure and HDX-MS analyses demonstrated that although 3a is a VDR antagonist, H11-H12 is stabilized. Surprisingly, H12 was strongly

stabilized by binding of antagonist 3a than 1,25D3. These structural analyses indicate 3a exhibits antagonistic 25

ACS Paragon Plus Environment

Bioconjugate Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

activity dependent on destabilization of the structure of H6-H7 but not H12. All VDR antagonists reported so far

are thought to destabilize H12 or be intact to H12, and no antagonist is known that stabilize H12. The same is

applicable for other NR antagonists. Therefore, compound 3a is a first antagonist against not only VDR but also

NRs that stabilizes H12. Thus, compound 3a stabilizes H11-H12 but crashes H6-H7 so that 3a is a novel rather than

“active” or “passive” type of antagonist. We classified 3a as a third type of antagonist and called it “H11-H12 stabilization antagonist”.

Antagonist 3b destabilized H11 and loop 11-12 and induced flexibility in H12 similar to apoVDR-LBD

(Table 1 and Figure 6c). Therefore, 3b is thought to inhibit VDR activity by preventing H12 folding. The docking study also demonstrated that it is reasonable that the 22-butyl group is accommodated in the butyl pocket despite

3b-induced prevention of H12 folding (Figure S4). This characteristic behavior (prevention of H12 folding)

explains why co-crystals of 3b were not obtained. As all VDR-LBD crystals reported to date have shown an H12-folding structure, H12 folding seems to be necessary for crystallization in terms of crystal packing. Antagonist

3b actively prevents H12 folding, thus preventing the formation of co-crystals of 3b. Binding of 3b also increased the deuterium levels in loop 6-7/H7. As reported in our previous studies, 22S-butyl compounds change the

side-chain conformation of residues in loop 6-7/H7 to induce butyl pocket formation but do not change the

conformation of the main chain. Therefore, 3b seems to destabilize loop 6-7/H7 without altering the conformation of the main chain. Taken together, these data indicate that the mechanism underlying destabilization of this region

upon 3b binding is distinct from that associated with 3a binding.

This study demonstrated that 3a is neither an “active” nor “passive” antagonist. Antagonist 3a is a third 26

ACS Paragon Plus Environment

Page 26 of 35

Page 27 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Bioconjugate Chemistry

type of antagonist, namely “H11-H12 stabilization antagonist”. In contrast, 3b is an active antagonist that prevents

H12 folding. It seems that 3a binds to the VDR-LBP, forming the activation function 2 (AF-2) surface upon H12

folding (Figures 3 and 4). Therefore, 3a binding might result in low-level coactivator recruitment (Figure S5). However, destabilization of regions H6-H7 and H10, which form part of the dimer interface, may inhibit

heterodimerization with RXR (Figure S5). Considering that the N-terminus of H10 becomes more unstable in

3a/VDR-LBD (16%) than 1,25D3/VDR-LBD (6%), H10 might adopt an unsuitable conformation for heterodimerization. In the case of antagonist 3b, 3b is thought to bind to the VDR-LBP but prevent H12 folding.

Therefore, the AF-2 surface is not formed and 3b inhibits coactivator recruitment (Figure S5). Moreover, 3b must

inhibit heterodimerization because, in addition to loop 6-7/H7, H10 becomes unstable, as in the case with 3a binding (Table 1, Figure S5).

The mechanism of the antagonistic activity of 3a differs from that of 3b. All VDR antagonists reported to date are side-chain analogues of 1,25D3. Those analogues must destabilize all or some of three regions via a similar or different mechanism than that described in this study. For example, binding of antagonist 1, which induces formation of the butyl pocket,17 must destabilize loop 6-7/H7 and the N-terminus of H10 by a mechanism

similar to that of 3b. However, unlike 3a and 3b, neither stabilization nor destabilization of H12 would occur

because there are no interactions between H12 and antagonist 1, as shown in the crystal structure of the 1/VDR-LBD complex.17 Thus, the flexibility of H12 in 1/VDR-LBD may be similar to that in the apo-form. These considerations are consistent with the observation that 1 does not recruit either RXR or the SRC-1 peptide.39 Based on these results, we conclude that the mechanism of VDR antagonistic behavior varies. 27

ACS Paragon Plus Environment

Bioconjugate Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

In conclusion, we designed and synthesized vitamin D analogues 2a,b-6a,b and evaluated their biological

activities. The 25-dibenzyl analogues 3a and 3b showed strong antagonistic activity. HDX-MS and X-ray crystal

structure analyses indicated that 3a is a “H11-H12 stabilization antagonist”. HDX-MS analysis also clearly showed

that 3b is an active antagonist. In this study we found first VDR antagonist 3a that strongly stabilizes H12. This finding must be applicable for other NRs and greatly affect drug design to treat NR related diseases. In addition, we

demonstrated that VDR antagonists have diverse mechanisms depending on each antagonist structure. Our results

also demonstrate that HDX-MS can provide information about dynamic changes in protein structure induced by

ligand binding that could prove useful in drug discovery research.

28

ACS Paragon Plus Environment

Page 28 of 35

Page 29 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Bioconjugate Chemistry

Experimental Section HDX-MS. A mixture of rVDR-LBD in buffer (10 mM Tris-HCl, pH 8.0, 200 mM NaCl, 2mM DTT) and ligand

(10 equiv.) was incubated for 1 h on ice. Each exchange reaction was initiated by incubating 8 µL of VDR/ligand

complex with 32 µL of D2O protein buffer for 10, 30, 60, 300, 900, or 3600 sec at 4°C. The exchange reaction was quenched by mixing with 60 µL of quench buffer (3 M urea, 1% TFA) on ice. Protein digestion was carried out with pepsin at an enzyme to substrate ratio of 1:1 for 5 min on ice. The digested samples were immediately frozen

in liquid N2. Peptide separation and mass spectrometry were carried out using an Agilent 6550 iFunnel Q-TOF LC/MS and MassHunter Workstation Qualitative Analysis Software (Agilent Technologies, CA, U.S.A.). Peptides

were separated using a JUPITER 4 µm Proteo 90 Å LC Column 50 × 4.6 mm (Phenomenex, CA, U.S.A.) with a linear gradient of 4–40% CH3CN/0.3% formic acid over 5 min at a flow rate of 1 mL/min. Peptide separation was performed within a thermal chamber held at 1°C to reduce D/H back exchange. The operating conditions for MS

analysis were as follows: positive ion mode, capillary voltage (VCap) 5000 V, gas temperature 220°C, drying gas

20 L/min, nebulizer gas 50 psig, fragmentor 400 V. Duplicates were performed for each on-exchange time point.

The intensity weighted average m/z value (centroid) of each deuterated peptide was calculated using MassHunter

Workstation Qualitative Analysis Software (Agilent Technologies) and Microsoft Excel. The deuterium level was

calculated according to the following equation: deuterium (%) = [(mt – m0%)/(m100% – m0%)] x 100, where mt, m0%, and m100% represent the centroid value of the partially deuterated, nondeuterated, and fully deuterated peptide, respectively.21

29

ACS Paragon Plus Environment

Bioconjugate Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Supporting Information The Supporting Information is available.

List of contents: Supplementary Figures and Tables, Supplementary Methods (Synthesis, Competitive Binding

Assay, Transfection and Transactivation Assay, Protein Expression and Purification, X-ray Crystallographic

Analysis, Docking Study, RXR and SRC-1 Recruitment), Supplementary References

Acknowledgement This work was financially supported by the Platform for Drug Discovery, Informatics, and Structural Life

Science from the Ministry of Education, Culture, Sports, Science and Technology (MEXT) and the Japan Agency

for Medical Research and Development (AMED), a Grant-in-Aid for Scientific Research (No. 26460155) from

MEXT, and the Takeda Science Foundation, Japan, to K.Y. Synchrotron-radiation experiments were performed at

the Photon Factory (Proposal No. 2013G656), and we are grateful for the assistance provided by the beamline

scientists at the Photon Factory.

AUTHOR CONTRIBUTIONS A.K. and K.Y. designed the research. K.Y. supervised the overall project. A.K. performed synthesis and biological

evaluation. A.K., T.I. and Y.A. performed x-ray crystal structure experiment and analyzed the data. A.K. and D.E.

performed the HDX-MS experiments and A.K. analyzed the data. A.K. and K.Y. wrote the paper. All authors

discussed the results and contributed to the final version of the manuscript.

30

ACS Paragon Plus Environment

Page 30 of 35

Page 31 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Bioconjugate Chemistry

References (1) Pike, J. W., Meyer, M. B., and Bishop, K. A. (2012) Regulation of target gene expression by the vitamin D receptor An update on mechanisms. Rev. Endocr. Metab. Disord. 13, 45–55. (2) Jurutka, P. W., Whitfield, G. K., Hsieh, J. C., Thompson, P. D., Haussler, C. A., and Haussler, M. R. (2001) Molecular nature of the vitamin D receptor and its role in regulation of gene expression. Rev. Endocr. Metab. Disord. 2, 203–216. (3) Orlov, I., Rochel, N., Moras, D., and Klaholz, B. P. (2011) Structure of the full human RXR/VDR nuclear receptor heterodimer complex with its DR3 target DNA. EMBO J. 31, 291–300. (4) Kubodera, N. (2010) Pharmaceutical studies on vitamin D derivatives and practical syntheses of six commercially available vitamin D derivatives that contribute to current clinical practice. Heterocycles 80, 83-98. (5) Nagpal, S., Na, S., and Rathnachalam, R. (2005) Noncalcemic actions of vitamin D receptor ligands. Endocr. Rev. 26, 662–687. (6) Menaa, C., Barsony, J., Reddy, S. V, Cornish, J., Cundy, T., and Roodman, G. D. (2000) 1,25-Dihydroxyvitamin D3 hypersensitivity of osteoclast precursors from patients with Paget’s disease. J. Bone Miner. Res. 15, 228–236. (7) Roodman, G., and Windle, J. (2005) Science in medicine-Paget disease of bone. J Clin Invest 115, 200–208. (8) Shiau, A. K., Barstad, D., Radek, J. T., Meyers, M. J., Nettles, K. W., Katzenellenbogen, B. S., Katzenellenbogen, J. A, Agard, D. A, and Greene, G. L. (2002) Structural characterization of a subtype-selective ligand reveals a novel mode of estrogen receptor antagonism. Nat. Struct. Biol. 9, 359–364. (9) Berry, M., Metzger, D., and Chambon, P. (1990) Role of the two activating domains of the oestrogen receptor in the cell-type and promoter-context dependent agonistic activity of the anti-oestrogen 4-hydroxytamoxifen. EMBO J. 9, 2811–2818. (10) Shiau, A. K., Barstad, D., Loria, P. M., Cheng, L., Kushner, P. J., Agard, D. A., and Greene, G. L. (1998) The structural basis of estrogen receptor/coactivator recognition and the antagonism of this interaction by tamoxifen. Cell 95, 927–937. (11) Bury, Y., Steinmeyer, A., and Carlberg, C. (2000) Structure activity relationship of carboxylic ester antagonists of the vitamin D3 receptor. Mol. Pharmacol. 58, 1067–1074. (12) Herdick, M., Steinmeyer, A., and Carlberg, C. (2000) Carboxylic ester antagonists of 1α,25-dihydroxyvitamin D3 show cell-specific actions. Chem. Biol. 7, 885–894. (13) Igarashi, M., Yoshimoto, N., Yamamoto, K., Shimizu, M., Ishizawa, M., Makishima, M., DeLuca, H. F., and Yamada, S. (2007) Identification of a highly potent vitamin D receptor antagonist: (25S)-26-adamantyl-25-hydroxy-2-methylene-22,23-didehydro-19,27-dinor-20-epi-vitamin D3 (ADMI3). Arch. Biochem. Biophys. 460, 240-253. (14) Kato, Y., Nakano, Y., Sano, H., Tanatani, A., Kobayashi, H., Shimazawa, R., Koshino, H., Hashimoto, Y., and Nagasawa, K. (2004) Synthesis of 1α,25-dihydroxyvitamin D3-26,23-lactams (DLAMs), a novel series of 1α,25-dihydroxyvitamin D3 antagonist. Bioorg. Med. Chem. Lett. 14, 2579–2583. (15) Miura, D., Manabe, K., Ozono, K., Saito, M., Gao, Q., Norman, A. W., and Ishizuka, S. (1999) Antagonistic action of novel 1α,25-dihydroxyvitamin D3-26, 23-lactone analogs on differentiation of human leukemia cells (HL-60) induced by 1α,25-dihydroxyvitamin D3. J. Biol. Chem. 274, 16392–16399.

31

ACS Paragon Plus Environment

Bioconjugate Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(16) Saito, N., Saito, H., Anzai, M., Yoshida, A., Fujishima, T., Takenouchi, K., Miura, D., Ishizuka, S., Takayama, H., and Kittaka, A. (2003) Dramatic enhancement of antagonistic activity on vitamin D receptor: A double functionalization of 1α-hydroxyvitamin D3 26,23-Lactones. Org. Lett. 5, 4859–4862. (17) Inaba, Y., Yoshimoto, N., Sakamaki, Y., Nakabayashi, M., Ikura, T., Tamamura, H., Ito, N., Shimizu, M., and Yamamoto, K. (2009) A new class of vitamin D analogues that induce structural rearrangement of the ligand-binding pocket of the receptor. J. Med. Chem. 52, 1438–1449. (18) Sakamaki, Y., Inaba, Y., Yoshimoto, N., and Yamamoto, K. (2010) Potent antagonist for the vitamin D receptor: Vitamin D analogues with simple side chain structure. J. Med. Chem. 53, 5813–5826. (19) Wales, T. E., and Engen, J. R. (2006) Hydrogen exchange mass spectrometry for the analysis of protein dynamics. Mass Spectrom. Rev. 25, 158–170. (20) Skinner, J. J., Lim, W. K., Bédard, S., Black, B. E., and Englander, S. W. (2012) Protein dynamics viewed by hydrogen exchange. Protein Sci. 21, 996–1005. (21) Zhang, J., Chalmers, M. J., Stayrook, K. R., Burris, L. L., Garcia-Ordonez, R. D., Pascal, B. D., Burris, T. P., Dodge, J. A., and Griffin, P. R. (2010) Hydrogen/deuterium exchange reveals distinct agonist/partial agonist receptor dynamics within vitamin D receptor/retinoid X receptor heterodimer. Structure 18, 1332–1341. (22) Zhang, J., Chalmers, M. J., Stayrook, K. R., Burris, L. L., Wang, Y., Busby, S. A, Pascal, B. D., Garcia-Ordonez, R. D., Bruning, J. B., Istrate, M. A, et al. (2011) DNA binding alters coactivator interaction surfaces of the intact VDR-RXR complex. Nat. Struct. Mol. Biol. 18, 556–563. (23) Carson, M. W., Zhang, J., Chalmers, M. J., Bocchinfuso, W. P., Holifield, K. D., Masquelin, T., Stites, R. E., Stayrook, K. R., Griffin, P. R., and Dodge, J. A. (2014) HDX reveals unique fragment ligands for the vitamin D receptor. Bioorganic Med. Chem. Lett. 24, 3459–3463. (24) Yan, X., Broderick, D., Leid, M. E., Schimerlik, M. I., and Deinzer, M. L. (2004) Dynamics and ligand-induced solvent accessibility changes in human retinoid X receptor homodimer determined by hydrogen deuterium exchange and mass spectrometry. Biochemistry 43, 909–917. (25) Boerma, L. J., Xia, G., Qui, C., Cox, B. D., Chalmers, M. J., Smith, C. D., Lobo-Ruppert, S., Griffin, P. R., Muccio, D. D., and Renfrow, M. B. (2014) Defining the communication between agonist and coactivator binding in the retinoid X receptor α ligand binding domain. J. Biol. Chem. 289, 814–826. (26) Hamuro, Y., Coales, S. J., Morrow, J. A, Molnar, K. S., Tuske, S. J., Southern, M. R., and Griffin, P. R. (2006) Hydrogen / deuterium-exchange (H/D-Ex) of PPARγ LBD in the presence of various modulators. Protein Sci. 15, 1883– 1892. (27) Bruning, J. B., Chalmers, M. J., Prasad, S., Busby, S. A., Kamenecka, T. M., He, Y., Nettles, K. W., and Griffin, P. R. (2007) Partial agonists activate PPARγ using a helix 12 independent mechanism. Structure 15, 1258–1271. (28) Hughes, T. S., Chalmers, M. J., Novick, S., Kuruvilla, D. S., Chang, M. R., Kamenecka, T. M., Rance, M., Johnson, B. A., Burris, T. P., Griffin, P. R., et al. (2012) Ligand and receptor dynamics contribute to the mechanism of graded PPARγ agonism. Structure 20, 139–150. (29) Dai, S. Y., Chalmers, M. J., Bruning, J., Bramlett, K. S., Osborne, H. E., Montrose-Rafizadeh, C., Barr, R. J., Wang, Y., Wang, M., Burris, T. P., et al. (2008) Prediction of the tissue-specificity of selective estrogen receptor modulators by using a single biochemical method. Proc. Natl. Acad. Sci. U. S. A. 105, 7171–7176.

32

ACS Paragon Plus Environment

Page 32 of 35

Page 33 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Bioconjugate Chemistry

(30) Dai, S. Y., Burris, T. P., Dodge, J. A., Montrose-Rafizadeh, C., Wang, Y., Pascal, B. D., Chalmers, M. J., and Griffin, P. R. (2009) Unique ligand binding patterns between estrogen receptor α and β revealed by hydrogen-deuterium exchange. Biochemistry 48, 9668–9676. (31) Kumar, N., Lyda, B., Chang, M. R., Lauer, J. L., Solt, L. A., Burris, T. P., Kamenecka, T. M., and Griffin, P. R. (2012) Identification of SR2211: A potent synthetic RORγ-selective modulator. ACS Chem. Biol. 7, 672–677. (32) Musille, P. M., Pathak, M. C., Lauer, J. L., Hudson, W. H., Griffin, P. R., and Ortlund, E. A. (2012) Antidiabetic phospholipid–nuclear receptor complex reveals the mechanism for phospholipid-driven gene regulation. Nat. Struct. Mol. Biol. 19, 532–537. (33) Kumar, R., Moure, C. M., Khan, S. H., Callaway, C., Grimm, S. L., Goswami, D., Griffin, P. R., and Edwards, D. P. (2013) Regulation of the structurally dynamic N-terminal domain of progesterone receptor by protein-induced folding. J. Biol. Chem. 288, 30285–30299. (34) Yamada, S., Morizono, D., and Yamamoto, K. (1992) Mild oxidation of aldehydes to the corresponding carboxylic acids and esters: alkaline iodine oxidation revisited. Tetrahedron Lett. 33, 4329–4332. (35) Inaba, Y., Yamamoto, K., Yoshimoto, N., Matsunawa, M., Uno, S., Yamada, S., and Makishima, M. (2007) Vitamin D3 derivatives with adamantane or lactone ring side chains are cell type-selective vitamin D receptor modulators. Mol. Pharmacol. 71, 1298–1311. (36) Anami, Y., Itoh, T., Egawa, D., Yoshimoto, N., and Yamamoto, K. (2014) A mixed population of antagonist and agonist binding conformers in a single crystal explains partial agonism against vitamin D receptor: Active vitamin D analogues with 22R-alkyl group. J. Med. Chem. 57, 4351–4367. (37) Nakabayashi, M., Tsukahara, Y., Iwasaki-Miyamoto, Y., Mihori-Shimazaki, M., Yamada, S., Inaba, S., Oda, M., Shimizu, M., Makishima, M., Tokiwa, H., et al. (2013) Crystal structures of hereditary vitamin D-resistant rickets-associated vitamin D receptor mutants R270L and W282R bound to 1,25-dihydroxyvitamin D3 and synthetic ligands. J. Med. Chem. 56, 6745–6760. (38) Shimizu, M., Miyamoto, Y., Takaku, H., Matsuo, M., Nakabayashi, M., Masuno, H., Udagawa, N., DeLuca, H. F., Ikura, T., and Ito, N. (2008) 2-Substituted-16-ene-22-thia-1α,25-dihydroxy-26,27-dimethyl-19-norvitamin D3 analogs: Synthesis, biological evaluation, and crystal structure. Bioorg. Med. Chem. 16, 6949–6964. (39) Inaba, Y., Nakabayashi, M., Itoh, T., Yoshimoto, N., Ikura, T., Ito, N., Shimizu, M., and Yamamoto, K. (2010) 22S-Butyl-1α,24R-dihydroxyvitamin D3: Recovery of vitamin D receptor agonistic activity. J. Steroid Biochem. Mol. Biol. 121, 146-150. (40) Yamamoto, K., Anami, Y., and Itoh, T. (2014) Development of vitamin D analogs modulating the pocket structure of vitamin D receptor. Curr. Top. Med. Chem. 14, 2378–2387.

33

ACS Paragon Plus Environment

Bioconjugate Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

34

ACS Paragon Plus Environment

Page 34 of 35

Page 35 of 35

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Bioconjugate Chemistry

39x21mm (300 x 300 DPI)

ACS Paragon Plus Environment