Reagents
MDCK cells were obtained from the American Type Culture Collection and cultured in DME (Invitrogen) containing 10% fetal bovine serum and 1% penicillin (Invitrogen). The immortalized HBECs (
Ramirez et al., 2004) were obtained originally from J. Minna (University of Texas Southwestern Medical Center, Dallas, TX) and were cultured in keratinocyte serum-free medium (Invitrogen). The following antibodies were used for immunoblot analysis: influenza virus proteins (B65141G; Biodesign/Meridian Life Science), NP (influenza virus nucleoprotein; ab20343; Abcam), NS1 (
Solórzano et al., 2005), M1 (
Bourmakina and García-Sastre, 2005), PB1 (SC-17601; Santa Cruz Biotechnology, Inc.), β-actin (Sigma-Aldrich), NXF1/TAP (tip-associated protein; Proteintech), Sp1 (Millipore), tubulin (Sigma-Aldrich), HIF1-α polyclonal antibody (obtained from J. Brugarolas, University of Texas Southwestern Medical Center, Dallas, TX; Bethyl Laboratories, Inc.), Nup98 (
Radu et al., 1995), and Rae1 (
Pritchard et al., 1999). Antibodies against VSV G, N, and P proteins were a gift from G. Barber (University of Miami School of Medicine, Miami, FL;
Balachandran and Barber, 2004). We generated the VSV M protein rabbit polyclonal antibody from full-length recombinant M protein. The plasmids pCMV-luciferase, expressing luciferase, and pCAGGS-NS1, expressing NS1, were used in the luciferase reporter gene assays. pEGFPN3-M-GFP was a gift from E. Izaurralde (Max Planck Institute for Developmental Biology, Tübingen, Germany). The NXF1 siRNA pool and control siRNA were obtained from Thermo Fisher Scientific.
Luciferase reporter gene assay
293T cells were transfected with a 10:1 ratio of plasmid pCMV-luciferase encoding luciferase and/or pCAGGS-NS1 encoding NS1 using Lipofectamine 2000 (Invitrogen). Cells were transfected with the luciferase plasmid alone and empty vector as a positive control. Luciferase activity was measured in a luminometer (Fluoroskan Ascent; Thermo Fisher Scientific;
Mata et al., 2011).
Synthesis and structure of quinoline carboxylic acids
Reactions were performed under air or dry argon in oven-dried glassware. All fine chemicals were obtained from Sigma-Aldrich.
1H and
13C spectra were recorded on a spectrometer (500; Varian) in DMSO-d
6 with tetramethylsilyl as the internal standard, chemical shifts (δ) are reported in parts per million downfield of tetramethylsilyl, and coupling constants (
J) are expressed in hertz. All the compounds were synthesized after a modified literature protocol (
Mueller and Stobaugh, 1950).
Synthesis of 6-chloro-2-(4-ethoxyphenyl)quinoline-4-carboxylic acid (1)
Into a flame-dried round-bottomed flask fitted with a reflux condenser was weighed 501 mg 4-chloroisatin (2.76 mmol, 1.0 equivalents). 5.5 ml ethanol was added to the flask. An ethanolic solution of 474 mg 4-ethoxyacetophenone (2.90 mmol, 1.05 equivalents) was added into the flask followed by the gradual addition of 447 mg KOH (dissolved in 500 µl H2O; 8.0 mmol, 2.9 equivalents). The flask was refluxed under an atmosphere of argon for 48 h. It was then cooled to room temperature, and the ethanol was removed using a rotary evaporator (Buchi). Into the residue was added 1 M aqueous hydrochloride acid until the pH became 2.0. The precipitate was collected by filtration, washed with brine, and dried under a vacuum. The crude product was purified via recrystallization from acetic acid obtained as light yellow powder with a yield of 452.3 mg (50%). 1H nuclear magnetic resonance (NMR; 500 MHz; DMSO-d6) was calculated as 1.38 (doublet [d], J = 8.1 Hz, 3H), 4.14 (quadruplet [q], J = 8.3 Hz, 2H), 7.10 (d, J = 8.0 Hz, 2H), 7.84 (d, J = 9.1 Hz, 1H), 8.12 (d, J = 8.3 Hz, 1H), 8.25 (d, J = 8.3 Hz, 2H), 8.48 (singlet [s], 1H), and 8.76 (s, 1H); 13C NMR (125 MHz; DMSO-d6) was calculated as 15.3, 64.0, 115.5, 120.7, 124.8, 125.1, 129.4, 130.6, 131.3, 132.2, 132.4, 137.2, 147.6, 156.6, 161.0, and 168.0. High resolution mass spectrometry (MS; HRMS) electrospray ionization (ESI−) calculated for C18H13ClNO3 (mass of protonated molecule [MH]) was 326.0589; found HRMS was 326.0593.
Synthesis of 6-chloro-2-(4-fluorophenyl)quinoline-4-carboxylic acid (1-4)
Into a flame-dried round-bottomed flask fitted with a reflux condenser was weighed 501 mg 4-chloroisatin (2.76 mmol, 1.0 equivalents). 5.5 ml ethanol was added to the flask. An ethanolic solution of 350 µl 4-fluoroacetophenone (2.90 mmol, 1.05 equivalents) was added into the flask followed by the gradual addition of 447 mg KOH (dissolved in 500 µl H2O; 8.0 mmol, 2.9 equivalents). The flask was refluxed under an atmosphere of argon for 48 h. It was then cooled to room temperature, and the ethanol was removed using a rotary evaporator. Into the residue was added 1 M aqueous hydrochloride acid until the pH became 2.0. The precipitate was collected by filtration, washed with brine, and dried under a vacuum. The crude product was purified via recrystallization from acetic acid obtained as light brown powder with a yield of 349 mg (42%). 1H NMR (500 MHz; DMSO-d6) was calculated as 7.40 (triplet [t], J = 8.8 Hz, 2H), 7.85 (dd, J = 9.1 and 2.3 Hz, 1H), 8.13 (d, J = 8.9 Hz, 1H), 8.34 (dd, J = 8.5 and 5.7 Hz, 2H), 8.50 (s, 1H), and 8.77 (s, 1H); 13C NMR (125 MHz; DMSO-d6) was calculated as 116.4 (d, J = 16.5 Hz), 120.7, 124.9, 125.2, 130.1 (d, J = 7.08 Hz), 131.1, 132.3, 132.8, 134.7, 138.0, 147.4, 155.8, 164.2 (d, J = 241 Hz), and 167.9. HRMS (ESI−) calculated for C16H8ClFNO2 (MH) was 300.0233; found HRMS was 300.0244.
Synthesis of 6-chloro-2-(4-(trifluoromethoxy)phenyl)quinoline-4-carboxylic acid (1-9)
Into a flame-dried round-bottomed flask fitted with a reflux condenser was weighed 501 mg 4-chloroisatin (2.76 mmol, 1.0 equivalents). 5.5 ml ethanol was added to the flask. An ethanolic solution of 463 µl 4-trifluoromethoxyacetophenone (2.90 mmol, 1.05 equivalents) was added into the flask followed by the gradual addition of 447 mg KOH (dissolved in 500 µl H2O; 8.0 mmol, 2.9 equivalents). The flask was refluxed under an atmosphere of argon for 48 h. It was then cooled to room temperature, and the ethanol was removed using a rotary evaporator. Into the residue was added 1 M aqueous hydrochloride acid until the pH became 2.0. The precipitate was collected by filtration, washed with brine, and dried under a vacuum. The crude product was purified via recrystallization from acetic acid obtained as light pink powder with a yield of 416.1 mg (41%). 1H NMR (500 MHz; DMSO-d6) was calculated as 7.50 (d, J = 8.3 Hz, 2H), 7.71 (d, J = 9.1 Hz, 1H), 8.03 (d, J = 9.1 Hz, 1H), 8.29 (s, 1H), 8.35 (d, J = 8.4 Hz, 2H), and 8.98 (s, 1H); 13C NMR (125 MHz; DMSO-d6) was calculated as 23.0, 118.6, 121.8, 126.4, 127.1, 129.8, 130.4, 131.1, 131.8, 138.4, 147.4, 148.6, 150.0, 155.5, and 173.8. HRMS (ESI−) calculated for C17H8ClF3NO3 (MH) was 366.0150; found HRMS was 366.0154.
Synthesis of 6-chloro-2-p-tolylquinoline-4-carboxylic acid (1-10)
Into a flame-dried round-bottomed flask fitted with a reflux condenser was weighed 501 mg 4-chloroisatin (2.76 mmol, 1.0 equivalents). 5.5 ml ethanol was added to the flask. An ethanolic solution of 388 mg 4-methylacetophenone (2.90 mmol, 1.05 equivalents) was added into the flask followed by the gradual addition of 447 mg KOH (dissolved in 500 µl H2O; 8.0 mmol, 2.9 equivalents). The flask was refluxed under an atmosphere of argon for 48 h. It was then cooled to room temperature, and the ethanol was removed using a rotary evaporator. Into the residue was added 1 M aqueous hydrochloride acid until the pH became 2.0. The precipitate was collected by filtration, washed with brine, and dried under a vacuum. The crude product was purified via recrystallization from acetic acid obtained as yellow powder with a yield of 434.5 mg (53%). 1H NMR (500 MHz; DMSO-d6) was calculated as 2.41 (s, 3H), 7.39 (d, J = 8.0 Hz, 2H), 7.86 (d, J = 9.0 Hz, 1H), 8.16 (d, J = 9.0 Hz, 1H), 8.20 (d, J = 8.5 and 5.7 Hz, 2H), 8.53 (s, 1H), and 8.77 (s, 1H); 13C NMR (125 MHz; DMSO-d6) was calculated as 21.6, 121.0, 124.8, 125.0, 127.8, 130.3, 131.2, 132.4, 132.8, 135.4, 136.5, 140.7, 147.6, 156.9, and 167.7. HRMS (ESI−) calculated for C17H11ClNO2 (MH) was 296.0484; found HRMS was 296.0498.
Synthesis of 6-fluoro-2-(4-phenoxyphenyl)quinoline-4-carboxylic acid (1-14)
Into a flame-dried round-bottomed flask fitted with a reflux condenser was weighed 500 mg 4-fluoroisatin (3.03 mmol, 1.0 equivalents). 5.5 ml ethanol was added to the flask. An ethanolic solution of 675 mg 4-phenoxyacetophenone (3.18 mmol, 1.05 equivalents) was added into the flask followed by the gradual addition of 494 mg KOH (dissolved in 500 µl H2O; 8.8 mmol, 2.9 equivalents). The flask was refluxed under an atmosphere of argon for 48 h. It was then cooled to room temperature, and the ethanol was removed using a rotary evaporator. Into the residue was added 1 M aqueous hydrochloride acid until the pH became 2.0. The precipitate was collected by filtration, washed with brine, and dried under a vacuum. The crude product was purified via recrystallization from acetic acid obtained in light brown powder with a yield of 468.2 mg (43%). 1H NMR (500 MHz; DMSO-d6) was calculated as 7.13 (multiplet [m], 4H), 7.21 (triplet [t], J = 7.4 Hz, 1H), 7.44 (t, J = 7.7 Hz, 2H), 7.72 (dt, J = 8.5 and 2.8 Hz, 1H), 8.16 (dd, J = 9.3 and 5.9 Hz, 1H), 8.27 (d, J = 8.7 Hz, 2H), and 8.47 (m, 2H); 13C NMR (125 MHz; DMSO-d6) was calculated as 110.1 (d, J = 28.2 Hz), 119.0, 120.0, 120.5, 120.6 (d, J = 26.0 Hz), 125.1, 125.2 (d, J = 12.5 Hz), 129.6, 130.9, 132.9 (d, J = 12.1 Hz), 133.4, 146.4, 155.3, 156.6, 159.2, 160.9 (d, J = 244.3 Hz), and 168.3. HRMS (ESI−) calculated for C22H13FNO3 (MH) was 358.0885; found HRMS was 358.0891.
DHODH purification
The gene encoding an N-terminal truncation of human DHODH was amplified from a cDNA library derived from human pituitary gland (Takara Bio Inc.;
Baldwin et al., 2002) and cloned into pET-28b (EMD) protein overexpression vector to generate a construct with a His
6 C-terminal tag. Proteins were expressed in Terrific Broth medium containing 50 µg/ml kanamycin (pET28b). Cells were grown to 0.8 OD
600 at 37°C, 0.2 mM IPTG was added to induce protein expression, and cells were grown overnight at 16°C. Cells were pelleted by centrifugation (4,000
g) and resuspended in lysis buffer (100 mM Hepes, pH 8.0, 150 mM NaCl, 10% glycerol, and 0.05% Thesit detergent [Fluka]), containing protease inhibitor cocktail for the His tag protein (Sigma-Aldrich). Cells were lysed by three passes through a high pressure homogenizer (EmulsiFlex-C5; Avestin, Inc.), the lysate was clarified by centrifugation (20,000
g), and the resulting supernatant was applied to a high performance purification column (HisTrap; GE Healthcare) precharged with Ni
+2. The column was sequentially washed with lysis buffer and lysis buffer containing 20 mM imidazole. DHODH was eluted from the column using a linear gradient from 20 to 400 mM imidazole. Fractions containing DHODH were pooled, concentrated with a concentrator (Amicon Ultra; Millipore), and then purified by gel filtration column chromatography on a Superdex 200 column (HiLoad 16/60; GE Healthcare) equilibrated with crystallization buffer (10 mM Hepes, pH 7.8, 100 mM NaCl, 1 mM
N,
N-dimethyldodecylamine
N-oxide [Fluka], 5% glycerol, and 10 mM DTT).
Enzyme assay and IC50 determination protocol
A dye-based assay that couples the oxidation of CoQ to the reduction of 2,6-dichloroindophenol (DCIP) was followed at 600 nm (molar extinction coefficient [e] = 18.8 mM−1cm−1). Enzyme stocks were diluted into an assay buffer (100 mM Hepes, pH 8.0, 150 mM NaCl, 10% glycerol, and 0.1% Triton X-100) containing 0.1 mM BSA to make a 100× working stock solution. Assays were initiated by adding enzyme into assay buffer containing substrates and inhibitors (final volume of 0. 5 ml) using the following conditions: DHODH (total enzyme concentration [ET] = 5–10 nM), substrates (0.2 mM l-dihydroorotate and 0.02 mM CoQD), 0.12 mM DCIP, and an inhibitor at 20°C. To determine half maximal inhibitory concentrations (IC50’s), the inhibitor concentration was varied in a threefold dilution series (0.01–100 µM). Data were collected on a spectrophotometer (DU 650; Beckman Coulter) using a six-cell auto sampler and an interval time of 12 s. The assay was followed for 6 min, and rate data (antibodies per minute) were extracted from the initial rates portion of the curve. Rate data (v) were converted to micromoles per minute using the molar extinction coefficient for DCIP. A77 1726 was a gift from P. Rathod (University of Washington, Seattle, WA).
Calculations. The percentage of inhibition relative to the no inhibitor control was determined (vi/vo × 100), and data were fitted to the log[I] versus response (three parameters) equation in Prism (GraphPad Software) to determine the IC50. Reported error represents the 95% confidence interval of the fit.
Model of 1-14 bound to human DHODH generated by DOCK
The structure of
1-14 was drawn with Sketcher through the CCP4 suite (
Collaborative Computational Project, Number 4, 1994) and then verified with COOT (Crystallographic Object-Oriented Toolkit;
Emsley et al., 2010). The x-ray structure coordinates (
2B0M;
Hurt et al., 2006) of human DHODH bound to the brequinar derivative ABBA was used for the docking analysis. First, ABBA was removed from the coordinates, and the receptor file of 2B0M(-ABBA) and the ligand file of inhibitor were prepared using program Chimera (
Pettersen et al., 2004). Then, the surface file of 2B0M(-ABBA) was generated with Chimera, and the sphere file around 2B0M(-ABBA) was generated by the program Sphgen under DOCK (University of California, San Francisco;
Kuntz et al., 1982). Next, a subset of spheres (clusters) of 2B0M(-ABBA) representing the binding site were selected and ranked using the program Showsphere under DOCK. The first and largest cluster was selected, and a box around the active site was generated by Showbox. Files of the rapid score evaluation for contact and energy scoring were generated using Grid. The Flexible Ligand option was used to run DOCK generating the structure of
1-14 docked to human DHODH as displayed in . The structure was displayed using PyMol.
Measurements of cellular UMP, UDP, and UTP levels
Metabolites were extracted from HBECs using cold (−20°C) 50% methanol. Extracts were dried down, resuspended in 10 mM tributylamine/15 mM acetic acid, pH 5.0, and filtered through a 0.2-µm spin column. The liquid chromatography MS/MS instrumentation used for quantitation of glutamine was an HPLC (Prominence LC20/SIL-20AC; Shimadzu) coupled to a triple quadruple mass spectrometer (3200 QTRAP; AB SCIEX). Metabolites were separated chromatographically on a C18-based column with polar embedded groups (150 × 2.0 mm at 4 µm; Synergi Fusion; Phenomenex) using a tributylammonium acetate/methanol gradient and subjected to MS/MS analysis—flow rate was 0.5 ml/min using the following method: buffer A (10 mM tributylamine adjusted with 15 mM acetic acid to pH 5.0) and buffer B (100% methanol). At the following times, concentrations were as follows: t = 0 min, 0% buffer B; t = 2 min, 0% buffer B; t = 13 min, 100% buffer B; t = 15 min, 100% buffer B; t = 16 min, 0% buffer B, t = 20 min, 0% buffer B; and t = 20.1 min, stop. The best multiple reaction monitoring (MRM) transitions in negative mode for UMP, UDP, and UTP were identified through quantitative optimization and indicated in . UMP, UDP, and UTP standards were injected to confirm MRM transitions and elution time. Peaks were quantitated and normalized against the total ion count using Analyst 1.4 software (Applied Biosystems).
Influenza virus replication
MDCK or HBEC cells were infected with A/WSN/1933 at an MOI of 0.001 or an MOI of 0.01, respectively, for 48 h in the absence or presence of the depicted compounds. Compounds were added 1 h after infection. Virus titers were measured by plaque assays. Assays were performed in triplicate from independent infections.
VSV replication assay
MDCK or HBEC cells were infected with VSV-GFP (Indiana strain) at an MOI of 0.001 in the absence or presence of 1 or its analogues at concentrations depicted in the figures. At 48 h after infection, supernatants were used for titration on BHK cells.
Oligo d(T) in situ hybridization
Poly(A)
+ RNA distribution was analyzed in HeLa cells or HBECs. To perform oligo (dT) in situ hybridization simultaneously with immunofluorescence (
Chakraborty et al., 2006), cells were fixed with paraformaldehyde, permeabilized with Triton X-100, and labeled with primary antibodies. The primary antibodies were diluted in PBS containing 0.2% Triton X-100, 1 mM DTT, and 200 U/ml RNasin (Promega). Cells were again fixed with paraformaldehyde and washed with PBS. Oligo (dT) in situ hybridization was then performed at 42°C overnight with a biotinylated oligo (dT) probe. Samples were washed with 2 and 0.5× SSC at 42°C. Cells were then fixed with formaldehyde, washed in PBS, and incubated for 30 min at room temperature with CY3-streptavidin (Sigma-Aldrich) and FITC-donkey anti–rabbit antibodies (Jackson ImmunoResearch Laboratories, Inc.), to detect NS1 protein, or FITC-donkey anti–goat antibodies (Jackson ImmunoResearch Laboratories, Inc.), to detect influenza proteins. Cells were washed in PBS, stained with DAPI, and mounted on glass slides using mounting medium (Dako). Images were acquired at room temperature with AxioVision 4.8.2.0 software on a camera (AxioCam MRm model 0445–554) in an inverted microscope (Axiovert 200M; Carl Zeiss). The Plan Neofluar 100×/1.30 NA objective oil lens was used. Images were processed using Photoshop and Illustrator (Adobe).
Generation of WSN:delNS1 virus
The NS1 segment carrying deletion of entire NS1 open reading frame (deletion of 58–516 nucleotides) was generated using fusion PCR and cloned into pDZ vector using the SapI site. WSN:del NS1 virus was generated using a reverse genetics system. In brief, 0.5 µg of each of eight pDZ plasmids representing the eight segments of the influenza A virus genome was transfected into 293T cells using Lipofectamine 2000. 24 h after transfection, the 293T cells were co-cultured with an NS1 complementing MDCK cells stably expressing NS1-GFP. The cell culture supernatants were passaged into fresh MDCK cells every 2 d and monitored for cytopathic effect. The successful rescue of virus was confirmed by performing a hemagglutination assay on the supernatants with chicken red blood cells. After three rounds of plaque purification and sequence confirmation, viral stocks were grown and titered in MDCK NS1-GFP cells.
NXF1 knockdown, nuclear/cytoplasmic fractionation, and quantitative PCR
293T cells were cotransfected with control siRNA and empty vector, control siRNA plasmids encoding NS1, NXF1 siRNA and empty vector, or NXF1 siRNA and plasmid encoding NS1 for 36 h. Then, cells were untreated or treated with 1 µM
1-14 for 24 h. RNA was purified from total cell extracts, nuclear fractions, and cytoplasmic fractions and subjected to real-time RT-PCR (
Wang et al., 2006).
Cytotoxicity measurements
Compound 1 was dissolved in sterile DMSO (Sigma-Aldrich) at a stock concentration of 10 mM. Compound 1 was then diluted in MEM (Invitrogen) to final concentrations depicted in . MDCK cells were infected with A/WSN/1033 at an MOI of 0.001. 1 h after infection, compound 1 was added. After 48 h of compound addition, cells were lysed, and ATP levels were measured by luminescence using a cell viability kit (CellTiter-Glo; Promega) following the manufacturer’s instructions.
Minigenome assay
293T cells were seeded in 24-well plates 24 h before transfection. Cells were transfected using Lipofectamine 2000 according to the manufacturer’s instructions with 100 ng pCAGGS-PB1, -PB2, and -PA, 200 ng pCAGGS-NP (of A/WSN/1933 or A/Viet Nam/1203/2004, respectively), 200 ng pLuc-Firefly-spNP, and 100 ng pCMV-renilla-luciferase. Compounds were added immediately after transfection. Activity of the reconstituted WSN polymerase complex was determined by luciferase reporter activity 24 h after transfection by dual luciferase assay (Promega) according to the manufacturer’s instruction.
We have constructed a hybrid
Homo sapiens interaction and reaction network by combining protein–protein interactions with directional signal transduction and metabolic reactions. Interaction information from IntAct (
Kerrien et al., 2007), NetworKin (
Linding et al., 2007), the Human Protein Reference Database (
Keshava Prasad et al., 2009), and from Palsson’s group have been used yielding a network of ~21,000 nodes (genes, proteins, and small chemicals) as well as ~132,000 interactions (gene–protein and protein–protein) and reactions (chemical, protein phosphorylation, etc.). The dataset was filtered to include only direct and physical interactions between human proteins, and all loops and duplicate edges were removed. Although, duplicate edges from different data sources and different property (e.g., an interaction identified as generic protein–protein interaction in one dataset and predicted as phosphorylation of a protein by a kinase in another dataset) were kept to emphasize the importance/validity of such interactions.
Microarray analysis
DHODH inhibitor did not alter bulk gene expression (Table S9). HBECs were incubated with 5 µM compound 1 for 3, 6, or 12 h or with compound 1-14 for 6 or 12 h. RNA was isolated using the RNeasy Mini kit (Invitrogen) following the manufacturer’s instructions. cDNA was synthesized, labeled, and hybridized to a HumanRef-8 BeadChip 22K (Illumina). After baseline correction and normalization of the expression data from the Illumina bead array, we further filtered the expression profile by omitting entries with a p-value of 0.05 or greater. Fold changes between test and reference sets have been calculated. For genes with multiple oligonucleotide probes, fold changes have been calculated before the calculation of means. Only test and reference pairs for each probe were used when both p-values were at or below the 0.05 cutoff value. The results are presented in the attached table Excel (Microsoft) file. Microarray data were deposited in GEO Datasets under accession no. GSE35117.
DHODH inhibitor regulates nucleocytoplasmic distribution of subsets of mRNAs in the presence of virulence factors (Table S10 and Table S11). Five million 293T cells were nontransfected or transfected with 6 µg pCAGGS-NS1 or pEGFPN3-M-GFP for 16 h. Then, cells were untreated or treated with 5 µM
1 for 24 h. RNA from total cell extracts or from nuclear or cytoplasmic fractions was obtained (
Wang et al., 2006). cDNA was synthesized, labeled, and hybridized to HumanRef-8 BeadChip 22K. After baseline correction and normalization of the expression data from the Illumina bead array, we further filtered the expression profile by restricting entries to possess expression values within 10% of duplicate measurements. For genes with multiple oligonucleotide probes, ratios have been calculated before the calculation of means. Only test and reference pairs for each probe were used when duplicate measurements for both expression values were within 10%. We then focused our analysis on mRNAs that did not have their levels altered by NS1 in the total cell extract but showed increased nuclear/cytoplasmic ratios in the presence of NS1 and reduced nuclear/cytoplasmic ratios upon compound treatment, indicating release of mRNA export block. We then subjected the postprocessed gene expression data to Gene Set Enrichment Analysis (GSEA;
Subramanian et al., 2005) using fold changes as ranks within the Prerank algorithm of the GSEA software with default parameters and 1,000 permutations. The curated C2/CP canonical pathways set v3.0 from the Molecular Signature Database at the Broad Institute has been used as reference gene sets for GSEA, consisting of sets of genes known to function in 880 pathways.
In parallel, we subjected highly responding genes for response network analysis. For this purpose, we have developed a computational method to identify response networks in large biological networks based on quantitative data (
Cabusora et al., 2005;
Mawuenyega et al., 2005). By superimposing gene expression data upon a large network, this method and the corresponding computer program NetworkExpress identifies
k-shortest paths (
Eppstein, 1998;
Hershberger et al., 2007) between a set of seed genes, scoring the subnetwork spanned by the set of
k-shortest paths that are shorter than a predefined maximum weighted length (
l) and finding the enriched response networks. The initial seed genes have been selected from a list of genes with the highest (positive) fold release of mRNAs into the cytoplasm by compound treatment in cells transfected with NS1 as compared with cells untreated with compound and transfected with NS1. The fold release of mRNAs into the cytoplasm was calculated between nuclear/cytoplasmic ratios of mRNAs in both aforementioned conditions. Response networks have been calculated by restricting the human network to genes with fold release > 1 (no change in nuclear/cytoplasmic ratio upon compound treatment). All other genes and corresponding interactions have been deleted from the network, yielding 3,400 nodes and 9,400 interactions. We scored this network with fold release > 1 values and calculated a corresponding response network by using 100 genes with the highest fold release ratio.
X-ray structure of 1-14
Crystals grew as pale yellow prisms by slow evaporation from acetic acid. The data crystal was cut from a larger crystal and had approximate dimensions of 0.33 × 0.23 × 0.08 mm. The data were collected on a diffractometer (SCXmini; Rigaku) with a charge-coupled device (Mercury) using a graphite monochromator with MoK-α radiation (λ = 0.71073 Å). A total of 1,800 frames of data were collected using ω scans with a scan range of 0.5° and a counting time of 19 s per frame. The data were collected at 223 K using a low temperature device (Tech50; Rigaku). Details of crystal data, data collection, and structure refinement are listed in
Tables S1, S2, S3, S4, S5, S6, and S7. Data reduction was performed using CrystalClear version 1.40 (Rigaku). The structure was solved by direct methods using SIR97 and refined by full-matrix least-squares on F
2 with anisotropic displacement parameters for the non-H atoms using SHELXL-97. The hydrogen atoms on carbon were calculated in ideal positions with isotropic displacement parameters set to 1.2× equivalent isotopic value (Ueq) of the attached atom (1.5× Ueq for methyl hydrogen atoms). The acidic hydrogen atoms on O2 and O2a were observed in a ΔF amp and refined with isotropic displacement parameters. The function, Σw(|F
o|
2 − |F
c|
2)
2, was minimized, in which w = 1/[(σ(F
o))
2 + (0.0611 × P)
2 + (0.5767 × P)] and P = (|F
o|
2 + 2|F
c|
2)/3. R
w(F
2) refined to 0.167, with R(F) equal to 0.0626 and a goodness of fit (S) = 1.02. Definitions used for calculating R(F), R
w(F
2), and the goodness of fit are given in this paragraph. The data were corrected for secondary extinction. The correction takes the form: F
corr = kF
c/[1 + (1.5(6) × 10
−6) × F
c2 λ
3/(sin2θ)]
0.25, in which k is the overall scale factor. Neutral atom scattering factors and values used to calculate the linear absorption coefficient are from the International Tables for x-ray Crystallography. All figures were generated using SHELXTL/PC. Tables of positional and thermal parameters, bond lengths and angles, torsion angles, and figures are found in the supplemental material. R
w(F
2) = {Σw(|F
o|
2 − |F
c|
2)
2/Σw(|F
o|)
4}
1/2, in which w is the weight given each reflection. R(F) = {Σ(|F
o| − |F
c|)/Σ|F
o|} was used for reflections with F
o > 4(σ(F
o)). S = [Σw(|F
o|
2 − |F
c|
2)
2/(n − p)]
1/2, in which n is the number of reflections and p is the number of refined parameters. F
o is the observed structure factor, and F
c is the calculated structure factor.
Online supplemental material
Fig. S1 shows the effects of
1 on cytotoxicity, viral protein expression, viral replication, and virus polymerase activity. Fig. S2 demonstrates the effects of
1 on influenza virus–mediated mRNA export block. Fig. S3 depicts the x-ray structure of
1-14. Table S1 shows crystallographic data for
1-14. Table S2 shows fractional coordinates and equivalent isotropic thermal parameters for the nonhydrogen atoms of
1-14. Table S3 shows bond lengths and angles for the nonhydrogen atoms of
1-14. Table S4 shows anisotropic thermal parameters for the nonhydrogen atoms of
1-14. Table S5 shows fractional coordinates and isotropic thermal parameters for the hydrogen atoms of
1-14. Table S6 demonstrates torsion angles for the nonhydrogen atoms of
1-14. Table S7 shows hydrogen bond lengths and angles for
1-14. Table S8 shows structure activity relationship experiments of quinoline carboxylic acids with respect to virus replication. Table S9 shows that the DHODH inhibitor did not alter bulk gene expression. Table S10 shows that DHODH inhibitor regulates nucleocytoplasmic distribution of subsets of mRNAs in the presence of NS1 protein. Table S11 shows that DHODH inhibitor regulates nucleocytoplasmic distribution of subsets of mRNAs in the presence of M protein. Online supplemental material is available at
http://www.jcb.org/cgi/content/full/jcb.201107058/DC1.