1
SUPPORTING INFORMATION
2
S1 TEXT
3
Potent allosteric dengue virus NS5 polymerase inhibitors: mechanism of action and
4
resistance profiling
5 6
Siew Pheng Lim1,*, Christian Noble1, Cheah Chen Seh1, Tingjin Sherryl Soh1,2, Abbas El
7
Sahili2, Kar Yarn Grace Chan2, Julien Lescar2,3, Rishi Arora4, Timothy Benson4, Shahul
8
Nilar1, Ujjini Manjunatha1, Kah Fei Wan1, Hongping Dong1, Xie Xuping1,#, Pei-Yong Shi1,#,
9
Fumiaki Yokokawa1
10 11
1
12
138670
13
2
14
Singapore 637551
15
3
16
Maladies Infectieuses. Centre Hospitalier Universitaire Pitié-Salpêtrière, Faculté de Médecine
17
Pierre et Marie Curie, Paris, France
18
4
19
02139, USA
20
#
21
Structural Biology & Molecular Biophysics, University of Texas Medical Branch, Galveston,
22
TX 77555, USA.
Novartis Institute for Tropical Diseases, #05-01 Chromos, 10 Biopolis Road, Singapore
School of Biological Sciences, Nanyang Technological University. 60 Nanyang Drive,
UPMC UMRS CR7 - CNRS ERL 8255-INSERM U1135 Centre d’Immunologie et des
Novartis Institute for Biomedical Research, 250 Massachusetts Avenue, Cambridge, MA
Present address: Department of Biochemistry & Molecular Biology, Sealy Center for
1
1 2
*To whom correspondence may be addressed:
[email protected]; Tel: 65-
3
67222924.
4
2
1
Supporting Methods.
2
Compound testing with DENV dnI FAPA assays. Final dnI assay reaction mixtures
3
contained 100 nM DENV 5´UTR-L-3´UTR IVT RNA, 20 µM ATP, 20 µM GTP, 20 µM
4
UTP, 5 µM ATTO-CTP, and 100 nM of DENV FL NS5 protein in 15 µl in assay buffer
5
comprising 50 mM Tris/HCl, pH 7.5, 10 mM KCl, 1 mM MgCl2, 0.3 mM MnCl2 (DENV2
6
and 4) or 1 mM MnCl2 (DENV1 and 3), 0.001% Triton-X-100 and 10 µM cysteine [1]. To
7
perform dose-response inhibition measurements, compounds from 0-20 or -100 µM
8
concentrations are two-fold serially diluted into 384-well black opaque plates (Corning
9
Costar), after which 100 nM DENV FL NS5 protein was added into respective wells. The
10
plates were sealed, briefly shaken and centrifuged at 1500 rpm for 30 sec and incubated at RT
11
for 15 min. Thereafter, RNA and ATTO-CTP, ATP, GTP and UTP were added to the wells to
12
start the reactions. The plates were re-sealed, shaken and centrifuged as before, followed by
13
incubation RT for 120 min. Reactions were stopped by addition of 10 µl of 2.5X STOP
14
buffer (200 mM NaCl, 25 mM MgCl2, 1.5M DEA, pH 10; Promega) with 25 nM calf
15
intestinal phosphatase (CIP), re- incubation at RT for 60 min and read on a Tecan Safire II
16
microplate reader. For order-of-addition experiments, DENV4 FL NS5 was incubated for one
17
hour at RT with RNA, ATP, and GTP or RNA, ATP, GTP and ATTO-CTP, followed by
18
exposure to serially diluted compounds for 20 min at RT. The missing components (ATTO-
19
CTP and UTP or UTP alone) were added and the reactions continued for 120 min after which
20
STOP buffer was added as before. All datapoints were performed in duplicate wells. Each
21
compound was tested at least twice.
22 23
Compound testing in DENV1-4 infected cells by high content imaging. Approximately 7
24
× 103 A549 cells per well were seeded into 384-well plates in Ham’s F-12K medium 3
1
containing 2 % FBS, 1mM L-glutamine and 1 % penicillin-streptomycin and incubated over-
2
night at 37 °C in 5 % CO2. At 24 h post-seeding, cells were infected with DENV1-4 at
3
multiplicity of infection [MOI] of 0.3-1 (DENV1, strain MY97-10245, MOI = 0.5; DENV2,
4
strain MY97-10340, MOI = 0.3; DENV3 strain MY05-34640, MOI = 1; DENV4, strain
5
MY01-22713; MOI = 0.5) and treated immediately with 10-point, 3-fold serially diluted
6
compounds (at final 0.5% DMSO concentration; 4). At 48 h post-infection, cells were fixed
7
with 4 % paraformaldehye and stained with 4G2 antibody conjugated with Dylight 488 (for
8
DENV envelope protein ) and DRAQ5 (for nuclear DNA). Quantifications of DENV
9
envelope protein were determined using the Opera high content imaging system and analysis
10
software (Perkin Elmer, USA). Calculation of EC50 values were performed using Helios
11
software package (Novartis, Basel).
12 13
Compound testing in DENV2 and HCV sub-genomic replicon assays. A549, BHK-21 and
14
Huh7 DENV2 sub-genomic replicon cells [2] were seeded at a density of 3,000 cells per well
15
in a 384-well microplate. After incubation at 37 °C with 5 % CO2 overnight, the cells were
16
treated with serially diluted compounds in a 10-point dose response starting from 20 or 50
17
µM concentrations. At 48 hr of post-incubation, renilla luciferase activities were measured
18
with the EndurRen live-cell substrate (Promega, USA) according to the manufacturer’s
19
protocol. Following luciferase activity measurements, the CellTiter-Glo reagent (Promega,
20
USA) was added to determine the cytotoxic effects of the compounds. For the HCV replicon
21
assay, Huh-7.5 cells harboring the HCV replicon [3] were seeded at a density of 20,000 cells
22
per well in a 96-well microplate. At 48 hr after compound treatment, the cells were assayed
23
for firefly luciferase activity by using a Bright-Glo luciferase assay (Promega, USA). NITD-
24
008, a nucleoside inhibitor of DENV and HCV [4], was included as a control. Compounds
4
1
were tested up to 25 µM in HCV replicon cells due to limits in DMSO tolerability of these
2
cells, and up to 50 µM in other cell types, based on compound solubility.
3 4
Compound testing in cells transiently expressing DENV and WNV replicon. Renilla
5
luciferase sub-genomic replicon of DENV2 (strain New Guinea C) or WNV (New York
6
strain 3356) cloned into the pACYC plasmid, was in vitro transcribed using a T7 mMessage
7
mMachine kit (Ambion, Austin, TX) after linearization with ClaI, as described previously
8
[5]. A549 cells were electroporated with 10 μg of replicon RNA using a GenePulser Xcell
9
system (Bio-Rad, Hercules, CA) with 3 pulses, at 450 V, 25 µF, at 5-10 sec intervals.
10
Transfected cells were seeded into a 96-well plate at 1.5 × 104 cells per well, followed by
11
treatment with the compound or the dimethyl sulfoxide (DMSO) control and incubation at
12
37°C in 5 % CO2. At various time points post-transfection, cells were washed once with
13
phosphate-buffered saline (PBS) and lysed in 20 μl 1× lysis buffer (Promega, USA) and
14
assayed for luciferase signals with a Clarity luminescence microplate reader (BioTek, USA)
15
using the renilla luciferase assay system (Promega, USA). CellTiter-Glo reagent was added to
16
the lysates to determine the cytotoxicity effects of compounds. Duplicate wells were seeded
17
for each data point.
18 19
Compound testing in cells transiently expressing infectious DENV. Full length infectious
20
virus of DENV2 (strain New Guinea C) cloned into the pACYC plasmid was in vitro
21
transcribed using a T7 mMessage mMachine kit (Ambion, Austin, TX) after linearization
22
with ClaI, as described previously [5]. BHK-21 cells were electroporated with 10 μg of
23
replicon RNA using a GenePulser Xcell system (Bio-Rad, Hercules, CA) with 3 pulses, at
24
850 V, 25 µF, at 5-10 sec intervals.
25
105 cells per well), followed by treatment with the compound or the dimethyl sulfoxide
Transfected cells were seeded into 96-well plates (1 ×
5
1
(DMSO) control and incubated at 37 °C in 5 % CO2. At various time points post-transfection
2
(p.t.), cells were washed once with phosphate-buffered saline (PBS) and lysed in 20 μl 1×
3
lysis buffer (Promega, USA) and assayed for luciferase signals with a Clarity luminescence
4
microplate reader (BioTek, USA). Cell viability was measured by using CellTiter-Glo kit
5
(Promega, USA) according to the manufacturer's protocols. Absorbance was measured at 450
6
nm by using a microplate reader (Tecan). Duplicate wells were seeded for each data point.
7 8
Plaque assay. Virus stock was produced by harvesting the supernatant of DENV2 full-length
9
RNA-transfected BHK-21 cells at 24, 48, 72, 96, and 120 hours post transfection. Virus titer
10
and morphology were determined by standard plaque assay. Briefly, a series of 10-fold
11
dilutions was prepared by first diluting 50 µl virus stock with 450 µl RPMI 1640 media
12
containing 2 % FBS to obtain 10-1 dilution and then further diluted until a final 10-6 dilution
13
was achieved. Confluent BHK-21 cells (1 × 105 cells per well, plated 2 days in advance)
14
grown in 24-well plate was added with 200 µl of each dilution per well. Duplicates were
15
prepared for each time point and dilution factor. The infection was allowed to proceed at
16
30°C for 1 hour, followed by removing the media and adding 500 µl 0.8% methyl-cellulose
17
overlay (containing RPMI, 2 % FBS, 1 % P/S, 0.05 % NaHCO3, 25 mM Na Hepes, and
18
0.5 % DMSO) into each well. The plates were incubated for 5 days at 37 °C in 5 % CO2
19
before fixing in 3.7 % formaldehyde and staining with 1 % crystal violet. The viral titer was
20
calculated as plaque-forming unit (PFU) per ml.
21 22
Immunofluorescence assay (IFA). Cells harbouring DENV replicons or infectious virus
23
were seeded into an eight-well Lab-Tek chamber slide (Thermo Fisher Scientific) and
24
incubated at 37 °C in 5 % CO2. At various time-points, cells were washed twice in PBS, and
25
then fixed in cold methanol at -20 °C. Cells were blocked with PBS containing 1 % bovine 6
1
serum albumin, 1 % FBS and 0.05 % Tween-20, and washed thrice with PBS followed by
2
incubation with anti-rabbit DENV NS5 (GeneTex) or anti-mouse dsRNA (English &
3
Scientific Consulting Kft., Hungary), After 1 hr, cells were washed thrice with PBS, and
4
added with secondary antibodies, Alexa Fluor 488 goat anti-rabbit IgG or Alexa Fluor 568
5
donkey anti-mouse IgG (Invitrogen, USA). After washes with PBS, the cells were mounted in
6
mounting medium with 4′,6-diamidino-2-phenylindole (DAPI; Vector Laboratories, Inc.).
7
Fluorescence images were acquired with a Leica DM4000 B system.
8 9
Data visualization and analysis. Assay development data were generated in Microsoft
10
Excel. Data analysis and visualization was done with GraphPad Prism5 and Spotfire. Positive
11
(maximum RFU) and negative (minimum RFU) controls were averaged and subtracted from
12
sample data to calculate % inhibition by using equation 1 below:
13 14
% inhibition = (100 – ((sample signal – average Negative )/(average Positive – average
15
Negative )×100))
16
The four parameter logistic equation was used for curve fitting to calculate IC50 according to
17
equation 2 below:
18 19
Y = Bottom + (Top -Bottom)/(1+((10^( logIC50 - X)^Hill slope))
20
where Bottom is the minimum Y value, Top is the maximum Y value, and Hill slope is the
21
slope of the linear portion of the semi-log curve. IC50 was extrapolated from logIC50
22
according to the GraphPad algorithm.
23
7
1
Co-crystallization of compounds with DENV3 FL NS5. Co-crystallization of either
2
compounds 27 or 29 with DENV3 FL NS5 was performed by hanging drop method mixing
3
1µl of mother liquor (16-18 % PEG 3350, 0.2 M Mg Acetate, 1 mM TCEP) to 1 µL of
4
protein mixture (2-3 mg/ml of protein with final concentration of 1-5 mM of each compound
5
diluted in 10 % DMSO).
6
8
1
Supporting Figures.
2
Fig. 1A. Representative IC50 curves for N-pocket inhibitors tested in DENV4 FL NS5 de
3
novo initiation FAPA assay [1]. Briefly, compounds (10-point, 3-fold serially diluted
4
compounds from 0- 20 or 100 µM) were incubated at RT for 20 min with enzyme alone, in
5
384-well plate opaque plates, after which reactions were started with addition of ssRNA and
6
nucleotide substrate components, and allowed to proceed for 2 hr. Reactions were stopped by
7
addition of 10 µl of 2.5X STOP buffer with 25 nM CIP, re- incubated at RT for 60 min and
8
read on a Tecan Safire II microplate reader (excitationmax and emissionmax wavelengths 422
9
nm and 566 nm). Data was fitted to the four parameter logistic equation and IC50 curves
10
plotted using Graphpad® Prism software. Table contains average IC50 values and hill slopes
11
obtained. All data points were measured in duplicates. 3 dGTP 15000
Best-fit values Bottom Top LogIC50 HillSlope IC50
12000 8000 4000
= 1980 22108 -0.2898 -0.9099 0.5131
A v e ra g e R F U
A v e ra g e R F U
16000
10000
5000
0
PKF024-680 = 2262 12785 2.896 -1.205 787.7
Bottom Top LogIC50 HillSlope IC50
12000 9000
-2
-1
0
1
2
3000
-1
0
1
2
3
4
Bottom Top LogIC50 HillSlope IC50
4000
LNG748 1917 9325 1.394 -1.196 24.75
0 -2
-1
0
lo g [c p d ],u M
(E) 15
1
2
3
-2
-1
lo g [c p d ],u M
0
1
2
3
0
1
lo g [c p d ],u M
(G) 26
(F) 25
NVP LXT442 NX1
24000
10000
24000
8000
20000
20000
A v e ra g e R F U
A v e ra g e R F U
8000
0
-2
lo g [c p d ]/u M
6000 Bottom Top LogIC50 HillSlope IC50
4000
2000
LOW911 1016 9248 0.1249 -0.6021 1.333
A v e ra g e R F U
-3
LNI017 = 2262 15349 1.490 -1.606 30.93
Bottom Top LogIC50 HillSlope IC50
6000
0 -4
A v e ra g e R F U
15000
20000
A v e ra g e R F U
12000
18000
24000
16000
Best-fit values Bottom Top LogIC50 HillSlope IC50
12000 8000 4000
0
= 1980 20930 -0.5515 -1.009 0.2809
16000
8000 4000
0
-3
-2
-1
0
1
2
3
= 1980 22628 -0.9410 -0.8121 0.1146
0 -4
-3
-2
lo g [c p d ],u M
-1
0
1
-4
-3
lo g [c p d ]/u M
(I) 27 18000
Best-fit values Bottom Top LogIC50 HillSlope IC50
12000
-2
-1
0
1
lo g [C p d ] /u M
(K) 29i
(J) 29
18000
20000
Best-fit values Bottom Top LogIC50 HillSlope IC50
6000
= 1980 16922 -0.7745 -1.058 0.1681
0
15000
Best-fit values Bottom Top LogIC50 HillSlope IC50
10000
5000
= 1980 18355 -0.7359 -0.9193 0.1837
-3
-2
-1
lo g [C p d ] /u M
0
1
16000 12000 Best-f it values Bottom Top LogIC50 HillSlope IC50
8000 4000
= 1980 19844 -1.535 -0.9626 0.02915
-4
-3
-2
-1
lo g [C p d ] /u M
0
1
-4
-3
-2
12000 Best-f it values Bottom Top LogIC50 HillSlope IC50
6000
= 1980 17707 -1.205 -0.8516 0.06241
0
0
0
-4
R F U 4 2 2 /5 6 6 n m
12000
R F U 4 2 2 /5 6 6 n m
A v e ra g e R F U
A v e ra g e R F U
20000
-1
lo g [C p d ] /u M
0
1
-4
-3
-2
-1
lo g [C p d ] /u M
12 13
9
1
Fig. 1B. Representative EC50 and CC50 curves for N-pocket inhibitors tested in HuH-7
2
DENV-2 replicon cell-based assay [2]. Cells were seeded over-night in white opaque 384-
3
well plates followed by incubation with increasing compound concentrations (10-point, 2-
4
fold serially diluted compounds from 0-50 µM) for 48 hr, after which cellular renilla
5
luciferase (EC50) or ATP (CC50) levels, measured as relative light units (RLU), were
6
determined using ViviRen and Cell Titer-Glo (Promega) according to manufacturer’s
7
protocol. All data points were measured in duplicates.
8 (B) 26 800000
A v e ra g e R L U
600000 80000 400000 40000 200000
CTG
800000
120000
A v e ra g e R L U
120000
600000 80000 400000 40000
V iv iR e n
V iv iR e n 0 -2
-1
0
1
2
0
0
0 -3
200000
CTG
-3
3
-2
-1
0
1
2
lo g 1 0 [ c p d ], µ M
lo g 1 0 [ c p d ], µ M
800000
140000
800000
105000
600000
105000
600000
70000
400000
70000
400000
CTG
35000
200000
A v e ra g e R L U
A v e ra g e R L U
(D) 27 140000
35000
V iv iR e n
0
0
0 -3
-2
200000
CTG
V iv iR e n
-1
0
1
0 -3
2
-2
lo g 1 0 [ c p d ], µ M
-1
0
1
2
lo g 1 0 [ c p d ], µ M
(F) 29i
400000 40000 200000
CTG
A v e ra g e R L U
A v e ra g e R L U
600000 80000
800000
120000
800000
120000
600000 80000 400000 40000
V iv iR e n
0
0 -3
9
-2
200000
CTG
V iv iR e n -1
0
lo g 1 0 [ c p d ], µ M
1
2
0
0 -3
-2
-1
0
1
2
lo g 1 0 [ c p d ], µ M
10 10
1
Fig. 2. Enzyme inhibition kinetics of N-pocket compounds against DENV polymerase.
2
DENV4 FL NS5 dnI FAPA assays [1] were performed in increasing concentrations of N-
3
pocket inhibitor, 15 or 29 or 3’dGTP (control) with 0-500 nM RNA or 0-50 µM GTP to
4
determine mechanism of inhibition with respect to either suibstrates. Representative
5
Michaelis-Menton plots were derived from non-linear regression curve fitting using
6
Graphpad Prism software.
7 8 500
R F U p e r m in
400
3 `d G T P (2 µ M ) 3 `d G T P (1 µ M )
300
3 `d G T P ( 0 .5 µ M ) 200
3 `d G T P (0 .2 5 µ M ) 3 `d G T P (0 .1 2 5 µ M )
100
R F U p e r m in
400
300
3 `d G T P (2 µ M ) 3 `d G T P (1 µ M ) 3 `d G T P ( 0 .5 µ M )
200
3 `d G T P (0 .2 5 µ M ) 3 `d G T P (0 .1 2 5 µ M )
100
3 `d G T P (0 .0 6 2 5 µ M )
3 `d G T P (0 .0 6 2 5 µ M )
0
0 0
100
200
300
400
500
0
5
R N A (n M )
10
15
20
25
G T P (µ M )
250
200
1 5 (2 0 µ M ) 150
1 5 (1 0 µ M ) 1 5 (5 µ M )
100
1 5 ( 2 .5 µ M ) 1 5 (1 .2 5 µ M )
50
R F U p e r m in
R F U p e r m in
200
150
1 5 (2 0 u M ) 1 5 (1 0 µ M )
100
1 5 (5 µ M ) 1 5 ( 2 .5 µ M ) 1 5 ( 1 .2 5 µ M )
50
1 5 (0 .6 2 5 µ M )
1 5 (0 .6 2 5 µ M ) 0
0 0
100
200
300
400
500
0
5
R N A (n M )
10
15
20
25
G T P (µ M )
2 9 -R N A 600
250
2 9 (1 2 5 n M ) 2 9 (6 2 .5 n M ) 2 9 ( 3 1 .2 5 n M )
200
2 9 (1 5 .6 2 5 n M )
R F U p e r m in
R F U p e r m in
200
2 9 (2 5 0 n M )
400
2 9 (2 5 0 n M ) 150
2 9 (1 2 5 n M ) 2 9 (6 2 .5 n M )
100
2 9 ( 3 1 .2 5 n M ) 2 9 (1 5 .6 2 5 n M )
50
2 9 (7 .8 1 2 5 n M ) 0 0
9
2 9 (7 .8 1 2 5 n M ) 0
100
200
300
R N A (n M )
400
500
0
5
10
15
20
25
G T P (µ M )
10 11 11
1
Fig. 3. Crystal structures of DENV3 FL NS5 were obtained by co-crystallization according
2
to conditions from [6] as described in Supplementary Materials. (A) Overall view of the
3
DENV3 FL NS5 structure displayed as ribbon with the methyltransferase domain in red, the
4
linker region in orange, palm, thumb and fingers subdomains colored in olive, green and blue
5
respectively. Both compounds 23 (magenta sticks) and 29 (yellow sticks) are overlaid in the
6
polymerase domain. Magnified views of compound 29 (B) and 27 (C) with Fo-Fc contoured
7
at 4 σ where each compound was omitted from the phase calculation. (D) Superimposition of
8
the RdRp domain from FL NS5 (green ribbon) and RdRp (pink ribbon) with bound
9
compounds 29 and 23 in sticks showing the absence of conformational changes. (E)
10
Superimposition of compound 29 bound to FL NS5 structure (magenta sticks) and bound to
11
RdRp structure (pink sticks) and compound 23 (panel F) bound to FL NS5 (grey sticks) and
12
bound to RdRp (orange sticks). The compound conformations are closely superimposable.
13
PDB codes for the FL NS5 structures with compounds 27 and 29 are 5JJS and 5JJR,
14
respectively.
15
12
1 2 3 4
13
1
Fig. 4. Phylogenetic tree representing relatedness of N-pocket amino acid residues from
2
different members of the Flavivirus family, derived from Clustal Omega program [8].
3 4 5 6
Fig. 5. Effects of N-pocket compounds on DENV polymerase thermo-stability. Melting
7
temperature (Tm) was assessed by thermo-denaturation in presence of the SYPRO Orange
8
dye as described in Materials and Methods. (A-F) Representative melting curves of DENV4
9
RdRp domain (aa 266-900; A-C) and FL NS5 (D-F) in presence of 50 µM N-pocket
10
inhibitors or 5 % DMSO control. (G) Cellular thermal shifts assays [9] were performed with
11
BHK-21 DENV2-NGC replicon cell lysates and 40 mM of 27 or 29. Briefly, lysates were
12
incubated with compounds or 5 % DMSO for 1 hr at 4 oC, followed by heating at 30-70 oC.
13
Samples were spun and the supernatants loaded onto 12 % SDS-PAGE gels, followed by gel
14
electrophoresis and western blotting with anti-DENV2 NS5 antibody (GeneTex, USA). Table
15
shows the changes in protein melting temperatures in presence of compounds compared to
16
controls treated with DMSO. 14
1 (A)
(B)
0
0
10
20
30
40
50
60
70
80
90
100
o
T e m p e r a tu r e ( C )
-1 0 0
- d R F U /d T
-1 5 0
0 10
-5 0
20
30
40
50
60
70
80
90
100
o
T e m p e r a tu r e ( C )
-1 0 0 -1 5 0
150 100
50
0
0
-5 0
10
20
30
40
50
60
70
80
90
100
o
-1 0 0
T e m p e r a tu r e ( C )
-1 5 0
- d R F U /d T
50
10
20
30
40
50
60
70
T e m p e r a tu r e ( C )
-1 5 0
-2 5 0
-3 0 0
-3 0 0
70
80
90
-1 5 0
DENV4 FL NS5 + 29
80
90
100
10 -5 0
20
30
40
50
60
70
80
T e m p (o C )
-1 0 0 -1 5 0 -2 0 0 -2 5 0
3 (G) DMSO control 30
32.5
35
37.5 40
42.5
45
30
32.5
35
37.5
42.5
45
47.5
50
27-40 µM 40
47.5
50
29-40 µM 19 32.5 30
35
37.5
100
DENV4 FL NS5 + DM SO
o
-2 5 0
60
0
-5 0
-2 0 0
50
o
50
D EN V4 FL N S5 + 2 9
-1 0 0
-2 0 0
40
T e m p e r a tu r e ( C )
100
DENV4 FL NS5 + DM SO
100
D EN V4 FL N S5 + 1 5
30
(F)
150
DENV4 FL NS5 + DM SO
20
-3 0 0
(E)
(D)
10
-2 5 0
-3 0 0
-3 0 0
-5 0 -1 0 0
-2 0 0
-2 5 0
-2 5 0
D EN V4 R dR p + 2 9
50
-2 0 0
-2 0 0
- d R F U /d T
D EN V4 R dR p + 2 7
- d R F U /d T
- d R F U /d T
50
D EN V4 R dR p + D M SO
100
DENV4 RdRp + 15
50
-5 0
150
D EN V4 R dR p + D M SO
100
100
2
(C)
150
D EN V4 R dR p + D M SO
- d R F U /d T
150
40
42.5
45
47.5
50
4 5 6 7 8 9 10
15
90
100
1
Fig. 6. Activity profiles of DENV polymerase bearing resistant phenotype amino acid
2
changes. Recombinant DENV2 (A, C) and DENV4 (B, D) FL NS5 proteins bearing single or
3
double amino acid changes in the N-pocket were tested in de novo initiation (A, B) and
4
elongation (C, D) FAPA assays and compared against activities of WT DENV2 or DENV4
5
FL NS5 proteins. Reactions were conducted over 2 hr at RT and from average relative
6
fluorescence units (RFU) obtained from one experiment. All data points were measured in
7
triplicate. (B)
(A) 30000
50000 DENV4 W T
d e n o v o in itia tio n
L511V E 802D 20000
L 5 1 1 V /E 8 0 2 D
10000
Q 803N 30000
L 5 1 2 V /Q 8 0 3 N
20000
10000
0 0 .0
L512V
40000
A v e ra g e R F U
A v e ra g e R F U
d e n o v o in itia tio n
DENV2 W T
0 0 .5
1 .0
1 .5
2 .0
0 .0
0 .5
T im e ( h r )
15000
8500
L511V E 802D
9000
L 5 1 1 V /E 8 0 2 D
6000
3000
0 0 .0
6800
A v e ra g e R F U
E lo n g a tio n A c tiv ity
DENV2 W T 12000
A v e ra g e R F U
1 .5
2 .0
1 .5
2 .0
(C)
(C) E lo n g a tio n A c tiv ity
1 .0
T im e ( h r )
DENV4 W T L512V Q 803N
5100
L 5 1 2 V /Q 8 0 3 N
3400
1700
0 0 .5
1 .0
T im e ( h r )
1 .5
2 .0
0 .0
0 .5
1 .0
T im e ( h r )
8 9 10 11 12 13 16
1
Fig. 7. Effects of compounds on DENV polymerase thermo-stability. Melting temperature
2
(Tm) was assessed by thermo-denaturation in presence of the SYPRO Orange dye as
3
described in Materials and Methods. (A-D) Representative melting curves of in vitro
4
expressed recombinant DENV2 FL NS5 WT or mutant proteins in presence of 50 µM
5
compound or 5 % DMSO control. Table shows the melting temperatures of DENV-2 and -4
6
FL NS5 WT and mutant proteins as well as the changes in their melting temperatures in
7
presence of compounds compared to controls treated with 5 % DMSO.
8
(A)
(B) 150
150
L511V + DM SO
W T + DMSO 100
W T + 29
50
0 20
40
60
-5 0
80
0
T e m p e ra tu re ,
d R F U /d T
d R F U /d T
50
100
W T + 27
100
C
-1 0 0
L511V + 27 L 511V + 29
0 20
40
60
-5 0
0
T e m p e ra tu re ,
80
100
80
100
C
-1 0 0
-1 5 0
-1 5 0
-2 0 0
-2 0 0 -2 5 0
(C)
(D)
150
150
E802D + DM SO 100
9
E 80 2D + 27
L 5 1 1 V /E 8 0 2 D + D M S O L 5 1 1 V /E 8 0 2 D + 2 7 L 5 1 1 V /E 8 0 2 D + 2 9
E 80 2D + 29
50
0 20
40
60
-5 0
T e m p e ra tu re ,
0
80
d R F U /d T
d R F U /d T
50
100
100
C
-1 0 0
0 20
40
60
-5 0
T e m p e ra tu re ,
0
C
-1 0 0
-1 5 0
-1 5 0
-2 0 0
-2 0 0
-2 5 0
-2 5 0
Protein melting temperatures ,Tm (oC) [Tm change compared to DMSO control] DENV2 FL NS5
DENV4 FL NS5
WT
L511V
E802D
L511V/E802D
WT
L512V
Q802N
L512V/Q802N
+ DMSO
35.5
35.5
36
35
37.5
36.5
37
36
+ 27
37 [1.5]
35.5 [0]
36 [0]
35 [0]
38.5 [1]
36.5 [0]
37.5 [0.5]
35 [-1]
+ 29
40 [4.5]
36.5 [1.0]
38 [2.0]
35.5 [0.5]
40.5 [3]
38 [1.5]
38.5 [1.5]
35.5 [-0.5]
10 17
1
Fig. 8. Analysis of viral and NS5 protein expressions from DENV WT and mutant
2
replicons and virus. Immuno-fluorescence stainings for DENV dsRNA and NS5 protein
3
were performed on BHK-21 cells at days 1-4 (D1-4) after electroporation of DENV-2 WT
4
and mutant (A) replicon or (B) full length viral IVT RNA. Cells were fixed and the stated
5
time-points and probed with mouse monoclonal anti-dsRNA (red; Scicons, USA) and rabbit
6
polyclonal anti-NS5 (green; GeneTex, USA) primary antibodies, and goat anti-mouse-IgG-
7
Alexa Fluor568 and goat anti-rabbit-IgG-FITC secondary antibodies (Invitrogen, USA).
8
Nuclear DNA was stained with DAPI (blue; Thermofisher, USA).
9
(A)
WT
L512V
E802D
L512V/E802D
D1
D2
D4 10 11 12 13 14 15 16 18
1
(B)
WT
L511V
E802D
L512V/E802D
D1
D2
D3
D4 2
19
1
Fig. 9. WNV (New York strain 3356) replicon cDNA was electroporated in BHK-21 cells,
2
after which cells were seeded into 96-well plates and treated with compounds, 26i, 27, 29 and
3
29i (10-point, 3-fold serially diluted compounds from 0- 50 µM), for 2 days. EC50 values
4
from replicon cells were determined by measuring cellular renilla luciferase levels. All data
5
points were measured in duplicates.
6
(A) 100000
26i 27
RFU
80000 60000 40000 20000 0 -3
-2
-1
0
1
2
Log (µM)
(B) 100000
29 29i
RFU
80000 60000 40000 20000 0 -3
7
-2
-1
0
1
2
Log (µM)
8 9 20
1
Table 1. Data collection and refinement statistics of DENV3 FL NS5 co-crystals.
2
PDB code
5JJR
5JJS
18 % PEG 3350, 0.2 M
18 % PEG 3350, 0.2 M Mg
Mg Acetate, 1 mM TCEP
Acetate, 1 mM TCEP
Wavelength (Å)
1.0
1.0
Resolution (Å)
50-1.99 (2.11-1.99)
50-1.65 (1.75-1.65)
Space group
P21212
P21212
Unit cell (Å)
94.6/150.4/69.1
94.7/150.2/68.5
Measured reflections
44525 (64401)
784305 (127110)
Unique reflections
67993 (10422)
225431 (36133)
Multiplicity
6.5 (6.1)
3.5 (3.5)
Completeness (%)
99.0 (95.2)
99.2 (98.3)
Mean I/signa I (I)
10.04 (3.01)
13.13 (1.65)
Rmerge (%)
19.0 (58.6)
6.01 (74.0)
CC (1/2) (%)
99.2 (83.5)
99.9 (64.3)
R-work (%)
17.9
17.1
R-free (%)
22.5
19.8
Macromolecule
6858
6920
Ligands
34 (L3) / 26 (SAH)
33 (L9) / 26 (SAH)
Solvent
473
628
Crystallization conditions
Number of non-hydrogen atoms
Protein residues
845
RMS (bonds) (Å)
0.010
0.010
21
RMS (angles) (°)
1.09
1.01
Ramachandran favoured (%)
99.3
99.3
Ramachandran outliers (%)
0.7
0.7
Macromolecules
23.2
26.9
Ligands
23.8
33.5
Solvent
28.2
35.9
Average B-factor (Å2)
1 2
Values in parenthesis are those for the last (highest resolution) shell
3
MR: Molecular replacement
4
CC1/2= percentage of correlation between intensities from random half‐dataset [7]
5 6 7 8 9 10 11 12 13
22
1
Table 2. Analysis of DENV2 (strain NGC) WT and mutant replicons and virus
2
replication.
3 4
(A)
5
kit (Promega, USA) from BHK-21 cells at days 1-4 after electroporation with DENV2 WT
6
and mutant replicons IVT RNA. Levels of secreted DENV2 (pfu/ml) were determined by
7
plaque assay from supernatants of BHK-21 cells at days 1-4 after electroporation with
8
DENV2 WT and mutant infectious viral IVT RNA.
Absolute relative renilla luciferase light units (RLU) measured using renilla luciferase
Renilla luciferase levels in
DENV2 levels (pfu/ml)
DENV2 replicon (RLU)
in culture supernatant
PostL512V/
L511V/ electroporation
WT
L511V
WT
E802D
L512V
E802D E802D
E802D (hr) 2276562.9
3061250.9
684487.8
2323574.8
± 14766.5
± 42578.2
± 9525.1
± 66485
1701190.5
2028336.0
1084181
1718929.8
±105595.8
± 136472
± 44392.6
± 25087.1
458693.7
618389.6
166496.1
467361.7
± 278.6
± 2494
± 577.4
± 12258.4
60752.2
141721.1
19744.6
67222.3
± 1488.5
± 5521.3
± 1073
± 9150
24
48
72
96
52500
1200
14000
65000
132500
55000
425000
775000
825000
200000
4500000
5250000
7000000
3250000
9750000
3750000
9 10 11
23
1
(B)
2
secreted virus titers, respectively, compared to WT replicon or virus levels (expressed as fold
3
increases over WT levels).
Mutant DENV2 replicons or virus replication measured by renilla luciferase levels or
Fold increase in DENV2 levels Renilla luciferase levels, fold increase
(pfu/ml) in culture
compared to WT replicon levels
supernatant compared to WT DENV2 levels
PostL511V/
L511V/ electroporation
L511V
L511V
E802D
E802D E802D
E802D (hr) 24
26.72
6.0
19.87
0.02
0.27
1.24
48
3.38
1.80
2.83
0.42
3.21
5.85
72
7.14
1.92
5.30
0.24
5.45
6.36
96
14.47
2.02
6.20
0.46
1.39
0.54
4 5 6 7 8 9 10 11 12 13 14 15 24
1
References
2
1. Niyomrattanakit P, Wan KF, Chung KY, Abas SN, Seh CC, Dong H, et al. (2015)
3
Stabilization of dengue virus polymerase in de novo initiation assay provides advantages
4
for compound screening. Antiviral Res 119: 36-46.
5
2. Ng C, Gu F, Phong WY, Chen Y, Lim SP, Davidson A, et al. (2007) Construction and
6
Characterization of a Stable Subgenomic Dengue virus type 2 Replicon System for
7
antiviral compound and siRNA screening. Antiviral Res. 76(3):222-231.
8 9
3. Blight KJ, Kolykhalov AA, Rice CM. 2000. Efficient initiation of HCV RNA replication in cell culture. Science 290(5498):1972–1974.
10
4. Yin Z, Chen YL, Schul W, Wang QY, Gu F, Duraiswamy J, et al. (2009) An adenosine
11
nucleoside inhibitor of dengue virus. Proc Natl Acad Sci USA 106(48):20435-20439.
12
5. Wang QY, Dong H, Zou B, Karuna R, Wan KF, Zou J, et al. (2015) Discovery of Dengue
13 14
Virus NS4B Inhibitors. J Virol. 89(16):8233-8244. 6. Zhao Y, Soh S, Zheng J, Chan KWK, Phoo WW, Lee CC, et al. (2015). A crystal
15
structure of the dengue virus NS5 protein reveals a novel inter-domain interface essential
16
for protein flexibility and virus replication. PLOS Pathogens 11(3):e1004682.
17 18
7. Karplus PA, Diederichs K. (2012). Linking crystallographic model and data quality. Science 336, 1030–1033.
19
8. Sievers F, Wilm A, Dineen DG, Gibson TJ, Karplus K, Li W, et al. (2011). Fast, scalable
20
generation of high-quality protein multiple sequence alignments using Clustal Omega.
21
Molecular Systems Biology 7: 539.
22
9. Molina MD, Jafari R, Ignatushchenko M, Seki T, Larsson EA, Dan C, et al. (2013)
23
Monitoring drug target engagement in cells and tissues using the cellular thermal shift
24
assay. Science 341: 84-7. 25