Fluorescent proteins such as eGFP lead to catalytic oxidative stress in cells Academic research paper on "Biological sciences"

BIOLOGY

Author's Accepted Manuscript

Fluorescent proteins such as eGFP lead to catalytic oxidative stress in cells

Douglas Ganini, Fabian Leinisch, Ashutosh Kumar, JinJie Jiang, Erik Tokar, Christine C. Malone, Robert M. Petrovich, Ronald P. Mason

www.elsevier.com/locate/redox

PII: S2213-2317(17)30153-2

DOI: http ://dx. doi.org/ 10.1016/j. redox.2017.03.002

Reference: REDOX597

To appear in: Redox Biology

Received date: 28 February 2017 Accepted date: 3 March 2017

Cite this article as: Douglas Ganini, Fabian Leinisch, Ashutosh Kumar, JinJii Jiang, Erik Tokar, Christine C. Malone, Robert M. Petrovich and Ronald P Mason, Fluorescent proteins such as eGFP lead to catalytic oxidative stress i cells, Redox Biology, http://dx.doi.org/10.10167j.redox.2017.03.002

This is a PDF file of an unedited manuscript that has been accepted fo publication. As a service to our customers we are providing this early version o the manuscript. The manuscript will undergo copyediting, typesetting, an< review of the resulting galley proof before it is published in its final citable form Please note that during the production process errors may be discovered whic could affect the content, and all legal disclaimers that apply to the journal pertain

Fluorescent proteins such as eGFP lead to catalytic

oxidative stress in cells'

Douglas Ganinia,i , Fabian Leinischa2,i, Ashutosh Kumara, JinJie Jianga, Erik Tokarb, Christine C. Malonec, Robert M. Petrovichc and Ronald P. Masona*

a. Free Radical Biology, Immunity, Inflammation & Disease Laboratory, National Institute of Environmental Health Sciences, National Institutes of Health, Research Triangle Park, NC 27709, USA.

b. Stem Cell Toxicology Group, National Toxicology Program Laboratory, National Institute of Environmental Health Sciences, National Institutes of Health, Research Triangle Park, NC 27709, USA.

c. Protein Expression Core Facility, Genome Integrity and Structural Biology Laboratory, National Institute of Environmental Health Sciences, National Institutes of Health, Research Triangle Park, NC 27709, USA.

Corresponding author: Ronald P. Mason, Free Radical Metabolites Group, Immunity, Inflammation & Disease Laboratory, National Institute of Environmental Health Sciences, National Institutes of Health, Research Triangle Park, NC 27709, USA. Telephone number: (919) 541-3910. mason4@niehs.nih.gov ABSTRACT

Fluorescent proteins are an important tool that has become omnipresent in life sciences research. They are frequently used for localization of proteins and monitoring of cells1, 2. Green fluorescent protein (GFP) was the first and has been the most used fluorescent protein. Enhanced GFP (eGFP) was optimized from wild-type GFP for increased fluorescence yield and improved expression in mammalian systems3. Many GFP-like fluorescent proteins have been discovered, optimized or created, such as the red fluorescent protein TagRFP4. Fluorescent proteins are expressed colorless and immature and, for eGFP, the conversion to the fluorescent form, mature, is known to produce one equivalent of hydrogen peroxide (H2O2) per molecule of chromophore5, 6. Even though it has been proposed that this process is non-catalytic and

* Footnote: This research was supported by the Intramural Research Program of the NIEHS, National Institute of Environmental Health Sciences/NIH.

1 Authors equally contributed to this work.

2 Current address: University of Copenhagen, Department of Biomedical Sciences, Cellular and

Metabolic Research Section, Blegdamsvej 3, 2200 Kobenhavn N, 4.5, Copenhagen, Denmark.

generates nontoxic levels of H2O2 6, this study investigates the role of fluorescent proteins in generating free radicals and inducing oxidative stress in biological systems. Immature eGFP and TagRFP catalytically generate the free radical superoxide anion (O2'-) and H2O2 catalytically in the presence of NADH. Generation of the free radical O2*- and H2O2 by eGFP in the presence of NADH affects the gene expression of cells. Many biological pathways are altered, such as a decrease in HIFIa stabilization and activity. The biological pathways altered by eGFP are known to be implicated in the pathophysiology of many diseases associated with oxidative stress; therefore, it is critical that such experiments using fluorescent proteins are validated with alternative methodologies and the results are carefully interpreted. Since cells inevitably experience oxidative stress when fluorescent proteins are expressed, the use of this tool for cell labeling and in vivo cell tracing also requires validation using alternative methodologies.

Keywords: GFP, H2O2 (hydrogen peroxide), O2*~ (superoxide free radical anion), free radicals, oxidative stress, redox biology.

INTRODUCTION

Fluorescent proteins are commonly used in biology for protein labeling and cell tracing. Green fluorescent protein (GFP) was the first and has been the most used fluorescent protein. Drs. Shimomura, Chalfie and Tsien were awarded a Nobel Prize in chemistry "for the discovery and development of the green fluorescent protein, GFP," in 20087. Enhanced GFP (eGFP) was developed for increasing the fluorescence yield and improving its expression in mammalian systems3. More recently, many optimized GFP-like fluorescent proteins have been discovered, developed or created, such as the red fluorescent protein TagRFP4.

The fluorescent chromophore of fluorescent proteins is formed through an intramolecular reaction between the side chains of certain amino acids localized inside the barrel structure of the protein5, 8. Even though it is known that chromophore formation requires cyclization and airmediated oxidation, proposed mechanisms for this complex reaction remain to be proved. It has been calculated that one molecule of hydrogen peroxide (H2O2) is generated independently of NAD(P)H during the maturation of each chromophore of eGFP. This non-catalytic H2O2 generation is within the cellular baseline level and, consequently, should not be toxic to cells or organisms6. Maturation of the fluorescent chromophore of TagRFP has been proposed to generate two molecules of H2O2 independently of NAD(P)H9, which is believed to be low and

nontoxic to cells. However, many reports show that eGFP expression is cytotoxic10-14 and, even though it is believed that it does not induce lethality in animals, other GFP-like fluorescent proteins and eGFP are shown to induce abnormalities in skeletal muscle15 and heart16. Since chemical fluorescent probes such as DCFH are known to catalyze the formation of superoxide anion free radical (O2*-) and H2O2 in biological systems by a variety of mechanisms17, we have searched for the formation of these reactive oxygen species by fluorescent proteins.

MATERIAL AND METHODS Chemicals and buffer

All chemicals were purchased from Sigma-Aldrich (St. Louis, MO) unless stated otherwise. Chelex-100 resin was purchased from Biorad (Hercules, CA). DMPO was purchased from Dojindo (Rockville, MD). Amplex Red was purchased from Invitrogen (Life Technologies, Grand Island, NY). To suppress trace metal contamination, phosphate buffer was Chelex-100-treated and contained 25 pM diethylene triamine pentaacetic acid (DTPA).

eGFP and TagRFP preparation in E. coli

Gateway-technology (Invitrogen, Thermo Fisher Scientific, Carlsbad, CA) was used to generate the protein-expression vectors. The pDEST vectors with eGFP and TagRFP genes were transformed by heat shock into Rosetta 2 DE3 pLacI competent cells. Ampicillin- and chloramphenicol-resistant colonies were grown in LB with the selecting agents. After cultures of OD 0.05 reached OD 0.6-0.8 (2h, 220 rpm at 37°C), protein expression was induced by the addition of 0.1 mM isopropyl ß-D-1-thiogalactopyranoside (IPTG); and then cultures were incubated for 4h at RT (220 rpm). eGFP and TagRFP were purified using an ÄKTA FPLC system (GE Healthcare Life Science, Pittsburgh, PA). Concentrations of the fluorescent proteins were determined based on protein content using the BCA Protein Assay kit from Pierce® (Rockford, IL) following the manufacturer's instructions. The molecular mass of eGFP was adopted as 27333.8 g/mol and that of TagRFP was adopted as 26732.5 g/mol. No difference was detected when eGFP or TagRFP concentrations were calculated using protein content or absorbance at 280 nm (e for eGFP = 21,890 M-1.cm-1 and £ for TagRFP = 27,640 M-1.cm-1). eGFP was expressed and purified 7 independent times.

Quantification of H2O2 in vitro

Samples containing eGFP and NADH were kept under continuous agitation (700 rpm) at room temperature, and H2O2 was quantified using the FOX1 assay18 adapted to a microplate format. Briefly, 5 pL of samples or authentic H2O2 standards (1-100 pM, £ at 240 nm = 43.6 M-1.cm-1) were added to 200 pL of FOX assay reagent (100 pM xylenol orange, 250 pM ammonium ironN sulfate, 100 mM sorbitol in 25 mM H2SO4) in wells of 96-well microplates. After 30 min, microplates were read for absorbance at 560 nm. The concentration of H2O2 was calculated after subtracting the background absorbance of wells containing 200 pL FOX1 assay reagent added of the 5 pL buffer blank.

eGFP preparation in HEK 293 cells

We prepared a pDEST26-eGFP plasmid using Gateway® technology from Invitrogen (Thermo Fisher Scientific, Carlsbad, CA). This plasmid makes an N-terminal 6xHis tagged eGFP under the control of the cytomegalovirus (CMV) promoter. FreestyleTM293-F cells (Thermo Fisher Scientific, Carlsbad, CA) were grown in FreestyleTM293 expression media as a suspension culture using Thomson Optimum Growth™ flasks (Oceanside, CA). For eGFP expression, one liter of 1 x106 cells/mL was transfected with 1 mg of pDEST26-eGFP in 2 mL of 1.6 mg/mL PEI MAX (Polysciences, Inc., Warrington, PA). Cultures were allowed to express eGFP for 48 hours. Cells were harvested by centrifugation, washed with PBS, flash frozen and stored at -80°C. For protein extraction, M-Per extraction reagent was added to frozen cell pellets (Thermo Fisher Scientific, Carlsbad, CA) containing 0.2 M NaCl, 5 mM imidazole and EDTA-free protease inhibitors (Roche, Indianapolis, IN). Lysates were prepared with 15 strokes of a Dounce homogenizer and sonication on ice, and insoluble material was removed by centrifugation at 28,000 g, for 30 min at 4°C. Cleared lysate was applied to a NiNTA-agarose column (QIAGEN, Valencia, CA) pre-equilibrated in wash buffer (50 mM Hepes, 100 mM NaCl and 10 mM imidazole, pH 7.5). An ÄKTA FPLC system was used (GE Healthcare Life Science, Pittsburgh, PA). The column was then washed with 10 column volumes of wash buffer containing 20 mM imidazole. Proteins bound to the column were eluted with a step gradient of imidazole in wash buffer (50, 100, and 200 mM). Fractions were analyzed by SDS-PAGE. Selected fractions were pooled and dialyzed into anion exchange buffer (50 mM Hepes and 10 mM NaCl, pH 7.5) before being applied to a HiTrap Q column (GE Healthcare Life Science, Pittsburgh, PA). eGFP was eluted with a linear gradient from 50 to 1000 mM NaCl over 20 column volumes. Pure fractions of eGFP were buffer exchanged into PBS, pH 7.4 and concentrated using 10 kDa cutoff centrifugal concentrators (Millipore, Billerica, MA). The purity of eGFP was examined by looking at the homogeneity of peaks in the FPLC-chromatograms and by SDS-PAGE.

ESR spectroscopy

ESR spectra were recorded with a Bruker ElexSys E-500 spectrometer (Bruker, Billerica, MA). The instrument operates at 9.80 GHz (X-Band). We used 20 mW microwave power, a modulation frequency of 100 kHz and a modulation amplitude of 1.0 G. The time constant and conversion time were 163.84 ms and the receiver gain was 1x104. For anaerobic measurements, the solutions were bubbled with nitrogen for at least 20 minutes. All experiments were repeated twice (n = 2) on three different days.

Expression of eGFP and MBP in E. coli

The gene for eGFP was recombined from pDONR 221 into a pDEST 527 expression plasmid (destination vector kindly supplied by Dr. Dominic Esposito, SAIC, MD). The pDEST 527 vector was transformed by heat shock into BL21-AI™ competent cells. The mbp gene was recombined from pDONR 221 into a pDEST 566-expressing plasmid and transformed by heat shock into Origami B DE3™ competent cells (Novagen, Darmstadt, Germany). Bacteria were routinely grown in Luria-Bertani broth with ampicillin (100 pg/mL). Cultures were prepared with Luria-Bertani broth at OD 0.05. After the cultures reached OD 0.6-0.8 (2h, 220 rpm at 37°C), IPTG (0.1 mM) was added and protein expression was induced by the addition of 0.002 and 0.02% arabinose (w/v). For the expression of MBP, IPTG alone (0.1 mM) was used as the inducer for high expression in bacteria. Cultures were then incubated for 2h at 37°C (220 rpm). Three independent cultures were prepared, and cell pellets were collected and frozen at -80°C for RNA isolation or Amplex Red assay for the determination of H2O2 production.

Bacterial RNA isolation and real-time quantitative PCR

Total RNA was isolated from three biological replicate cultures of bacterial cells with the use of the TRIzol Max Bacterial RNA Isolation Kit from Ambion (Life Technologies, Grand Island, NY) and purified with the use of RNeasy Mini Kit Columns from QIAGEN (Valencia, CA) according to the manufacturer's instructions. Purified RNA was reverse transcribed to complementary DNA with the use of random hexamers and Moloney Murine Leukemia Virus Reverse Transcriptase from Applied Biosystems (Foster City, CA). Gene-specific primers (Supplementary Data Table S2) were designed using the NCBI/Primer BLAST Designing Tool

(https://www.ncbi.nlm.nih.gov/tools/primer-blast/). The ABsolute SYBR Green ROX Mix from Abgene (Rockford, IL) was used for amplifications. Conditions for amplification were as follows: 15 minutes at 950C, followed by 40 cycles of 950C for 1 minute and 600C for 1 minute. Cycle time (Ct) values for the selected genes were normalized to values for 16S rRNA (n = 3 in pseudo triplicates). Ct values for all controls were set at 100%.

Mammalian cell culture

The HeLa cell line stably expressing eGFP (HeLa/eGFP) was purchased from Cell Biolabs, Inc. (San Diego, CA) and grown in parallel with wild-type HeLa cells. Those cell lines were certified to be free of mycoplasma. HeLa.tet and HeLa.tet.eGFP were purchased from GenTarget Inc. (San Diego, CA) and were also certified to be free of mycoplasma. Cells were cultured in high

glucose Dulbecco's Modified Eagle's Medium, supplemented with 10% heat-inactivated fetal bovine serum (v/v), 50 U/mL penicillin, 50 pg/mL streptomycin, 2 mM L-glutamine and 0.1 mM MEM nonessential amino acids (NEAA). Cells were maintained in a humidified atmosphere of 5% CO2 at 37°C and cultured in T-75 Falcon tissue culture flasks at an initial density of 1x106 cells/flask. The culture medium was replaced every 2 days and passaged when confluence was reached. Three independent cell cultures were prepared. Cultures of HeLa.tet and HeLa.tet.eGFP at 70% confluency were exposed to 1 ug/mL doxycycline for 48 h before collection. Confluent cells were collected for RNA isolation by trypsinization, washed with PBS and stored at -80°C.

Amplex Red assay for determination of H2O2 production by HeLa cells

Mammalian cells were trypsinized and 1x105 cells were seeded to a black 96-well plate, incubated overnight and washed with PBS. Wells with cells (n = 8 for each cell line) were filled with 100 pL of 50 pM Amplex Red, 0.2 U/mL HRP in PBS and 0.5% DMSO (v/v) from the Amplex Red stock solution. Catalase was added to control wells (1 pl/well, 1 kU/mL). Exposure to light was avoided by preparing the assay in a dark room19. A microplate plate reader was prewarmed and kept at 37°C, and fluorescence intensities (filter: Aex = 530 nm, Aem = 590 nm) were measured from the top of each well every 10 minutes for 120 minutes. E. coli was also used for this assay, but Luria Brentani broth medium was used at 1.0 of OD at 600 nm/mL of bacterial cells. Experiments were repeated three independent times with mammalian cells and two independent times with bacterial cells.

RNA isolation from HeLa cells, microarray and GSEA analyses

Total RNA was isolated from three independent triplicates using the QIAshredder and RNeasy Midi® Kit (QIAGEN, Valencia, CA) with in-column DNAse treatment (QIAGEN, Valencia, CA). Gene expression analysis was conducted using Agilent Whole Human Genome 4x44 Multiplex Arrays (014850, Agilent Technologies, CA). Total RNA samples were labeled with Cy3 according to the manufacturer's protocol. Data were obtained using the Agilent Feature Extraction Software (v. 12). Data was further processed using R (version 3.1.3) and RStudio (version 0.98.1103). Raw microarray data was normalized and the background subtracted using the limma package (version 3.22.6). Annotation was made using the hgug4112a database. Heat maps were prepared in R using the gplot package (version 2.16.0).

Two different 2-contrast GSEA analyses20 were conducted: 1. HeLa/eGFP versus HeLa and 2. the differences in the log2-transformed intensities between HeLa.tet.eGFP incubated with and

without doxycycline versus the differences of HeLa.tet incubated with and without doxycycline. Both analyses were performed with the full table using the hallmark database, h.all.v5.0.symbols, 1,000 permutations, and defaults for other parameters. The microarray raw data for (a) HeLa cells stably expressing eGFP can be accessed at

https://figshare.com/s/2523da57680cffaf1a0c; (b) HeLa cells expressing eGFP in the presence of doxyxycline can be accessed at https://figshare.com/sZ5dbf0a84b2662ea718f4; and (c) the tables (list of all probe IDs and their intensities) used in the GSEA analyses can be accessed at https://figshare.com/s/2c99caae4c27d410ce6b.

Statistical analyses

Experiments were routinely prepared in three independent replicates (n = 3) and repeated on two different days. Mean and error bars as standard deviation are shown in graphs. ANOVA statistics and the post-hoc two-sided t-test were used for comparisons between groups. Comparisons with p-value lower than 0.05 were considered statistically significant.

RESULTS

Generation of O2'- and H2O2 by eGFP and TagRFP in the presence of NAD(P)H

Initially we used purified recombinant eGFP expressed in E. coli. As shown in Fig. 1a, we detected for the first time the generation of the free radical O2'- in samples of eGFP and NADH using ESR spin trapping with 5,5-dimethyl-1-pyrroline W-oxide (DMPO). The signal was a superimposition of the adducts DMPO/'OH and DMPO/'OOH (Fig. 1b; analysis using spectral simulation21), was present for over 30 min and was totally dependent on eGFP and NADH. DMPO/'OOH is known to decay to DMPO/'OH by a number of pathways. The signal was sensitive to superoxide dismutase (SOD), but not to catalase, as expected. We demonstrated that molecular oxygen is required for O2'- formation (Fig. 1c) since the consumption of NADH did not occur in the absence of oxygen (Fig. 1d). No change in the light absorption for the chromophore of mature eGFP was detectable in samples either aerobically or anaerobically prepared (Fig. 1d). The reaction restarted as soon as the sample was oxygenated.

Generation of O2'- by eGFP in the presence of NADH was also determined independently by the cytochrome c reduction assay (Fig. 2a, O2'--mediated reduction of cytochrome c is followed optically). As expected, cytochrome c reduction was inhibited by SOD but was insensitive to catalase (Fig. 2a, red bars). Cytochrome c did not interfere with the rate of NADH consumption by eGFP (Fig. 2a, black bars). Formation of O2'- was calculated over a range of NADH concentrations and was formed at approximately one third of the rate of NADH consumed (Fig. 2b). By using the same eGFP preparation, we also showed that the consumption of NADH was nearly equimolar to the formation of H2O2 (412 nM/min of NADH consumption and 388 nM/min of H2O2 production, Fig. 2c). The catalytical activity of eGFP in this reaction was further tested in samples containing 25 pM eGFP and 1 mM NADH (Supplemental Data Fig. S1a). The concentration of H2O2 was much greater than stoichiometric at the end of this experiment (478 pM H2O2 or 19 equivalents to total eGFP). In conclusion, eGFP catalytically consumes NADH to produce O2'- and H2O2 with an overall stoichiometry of one NADH oxidized to NAD+ and one O2 reduced to H2O2. It is noteworthy that in all experiments, NADH could be substituted for NADPH with no changes in the chemical reaction (data not shown).

Interestingly, different preparations of our extensively purified eGFP at the same total protein concentration showed different rates of NADH consumption (Supplemental Data Fig. S1b). These preparations had similar quantities of mature eGFP per total protein content (Supplemental Data Fig. S1c) and showed NADH consumption dependent on both NADH and

eGFP (Supplemental Data Fig. S1d). By calculating the amount of mature eGFP chromophore per protein, we found that approximately 75% of the total eGFP is in its chromophore-mature state. Apparently, the mature eGFP is not the generator of O2'-.

We next used extensively purified eGFP expressed in the mammalian cell line HEK 293 (Fig. 2d). This protein was found approximately 85% in its chromophore-mature form. eGFP expressed in HEK 293 also generated O2*- in the presence NADH, as shown by the inhibition of cytochrome c reduction when SOD was present (first and second columns, Fig. 2d). Formation of O2*- was dependent on NADH (third column, Fig. 2d).

We also tested the formation of O2*- and consumption of NADH by a different fluorescent protein, the red fluorescent protein TagRFP. Highly purified TagRFP produced in E. coli consumed NADH and generated O2*- (Supplemental Data Fig. S1e). Purified TagRFP was approximately 90% in its chromophore-mature form (£ at 555 nm = 100,000 M-1cm-1 4). The yield of O2- generation by TagRFP lies between the higher yields found for eGFP produced in E. coli and the lower yields for the eGFP produced in HEK 293.

Cells expressing eGFP show oxidative stress and changes in many biological processes

In order to investigate how eGFP can alter gene expression associated with oxidative stress in living systems, we used first the bacterium E. coli to monitor the expression of genes under the control of the transcription factors SoxR and OxyR. SoxR is activated by direct reaction with superoxide anion free radical22, whereas OxyR induces gene expression when key cysteines are oxidized as a result of an oxidizing intracellular redox status23. Experiments were performed with E. coli expressing eGFP under low and high levels. All genes expressed upon activation of SoxR (Fig. 3a) and OxyR (Fig. 3b) were highly upregulated in these cells, proportionally to eGFP expression. A bacterium expressing maltose binding protein (MBP) was used as a control for the possible oxidative stress involved in overexpressing a protein. Notably, the high level of expression of MBP (Fig. 3c, first two columns) did not lead to changes in expression for any of the genes activated by SoxR or OxyR (Fig. 3c, left two panels). Even though the differences in gene expression among the bacteria expressing eGFP were striking, we did not detect differences in extracellular H2O2 production (Supplementary Data Fig. S2). This negative result can be explained by the fact that E. coli is very resistant to H2O2. Even at 1 mM H2O2, E. coli does not lose viability, presumably because this bacterium is equipped with very efficient systems for H2O2 decomposition24.

Next, we studied the effect of eGFP expression in mammalian cells using the HeLa cell lines. We discovered that HeLa cells stably expressing eGFP (HeLa/eGFP, Fig. 4a) show

significantly higher rates of H2O2 efflux compared to regular HeLa cells (7.3 and 5.15 nM H2O2/min, respectively), indicating a higher generation of intracellular H2O2 (Fig. 4b). The H2O2-decomposing enzyme catalase inhibits the Amplex Red oxidation to levels that do not differ between cell lines (1.8 and 2.3 nM H2O2/min, respectively; black traces on Fig. 4b).

Production of H2O2 was further investigated using a HeLa cell line engineered to express eGFP under the control of a tet-mediated promoter (HeLa.tet.eGFP). In this cell model, eGFP is expressed only when doxycycline is added to the cell culture medium (Supplemental Data Fig. S3a). Possible non-specific effects of tet and/or doxycycline were taken into account by analyzing the response of its isogenic parental cell line (HeLa.tet) exposed to doxycycline. We compared HeLa.tet.eGFP cells with induced expression of eGFP to the isogenic parental cell line (HeLa.tet) exposed to doxycycline and found no detectable difference in the low rate of H2O2 generation in both cultures. The level of eGFP in the Hela.tet.eGFP exposed to doxycycline was approximately 13 times lower than the level in the stable eGFP-expressing cell line, HeLa/eGFP. However, when we compared the gene expression and GSEA (Gene Set Enrichment Analyses20) of the HeLa/eGFP versus the control HeLa cells, and the specific effects of eGFP expression in HeLa.tet.eGFP versus Hela.tet, we saw that both eGFP-expressing cell models showed significant gene set enrichment for nine biological processes (gene sets with p-values lower than 0.001, Fig. 4c). Three gene sets associated with pathogen response and inflammation are upregulated and are related to interferon-a response, interferon-Y response, and TNF-a signalling via NFkp (Fig. 4c and Supplemental Data Fig. 3b-d), which are classically known to be activated when cells are producing intracellular H2O2 and suffering from oxidative stress25. Another six gene sets are downregulated in HeLa cells that express eGFP (Fig. 4c and Supplemental Data Figure 3e-j), including hypoxia. In agreement with our analyses of global gene expression, we found three times lower levels of HIF-1a, the master regulator for hypoxia, in HeLa/eGFP compared to HeLa (Fig. 4d and e). It has been reported that under normoxia, the intracellular level of HIF-1a is decreased in response to the toxic effects of higher intracellular H2O2 and O2'- formation26, 27, such as observed in HeLa/eGFP. This cell line secreted nearly 500 nM of H2O2 into the medium in the course of one hour (Fig. 4b).

CONCLUSIONS

Our in vitro experiments show that eGFP in the presence of NAD(P)H generates O2*-and H2O2. Collectively, the absence of changes in the visible spectrum of the eGFP-chromophore, the lack of correlation between the concentration of mature chromophore and rates for the reaction with NADH in different batches of eGFP lead us to conclude that an inevitable small portion of non-fluorescent or immature chromophore present in eGFP is responsible for the catalytical consumption of NADH with O2*- and H2O2 formations. This is consistent with the finding that the mimic probe for the mature eGFP chromophore, 3,5'-difluoro-4-hydroxybenzyldene imidazolinone (DFHBI, Lucerna, Inc., New York, NY), fails to consume NAD(P)H (data not shown).

Even though we do not have direct experimental evidence to support a definitive mechanism for this reaction, we speculate that the last immature intermediate in the proposed chromophore maturation process28 reacts with NAD(P)H. This reaction would result in the regeneration of an earlier immature intermediate, which ultimately leads to the catalytical consumption of NAD(P)H and generation of H2O2, as shown in Scheme 1 for the chromophore of eGFP. It is very likely that the maturation of the last immature intermediate is slow since tertiary hydroperoxides are known to be highly stable, e.g., tert-butyl hydroperoxide.

More importantly, this chemistry likely occurs inside of bacteria and mammalian cells expressing eGFP. Our experiments were carried out with 50, 25 and 10 pM eGFP or 10 pM TagRFP, which is comparable to the concentration found in different cellular systems29, 30. NAD(P)H is present at mM concentrations in biological systems; therefore, the chemistry shown here is likely to occur in most biological experiments, which, in turn, leads to oxidative stress and cellular adaptations in gene expression.

The level of extracellular H2O2 detected in HeLa cells stably expressing eGFP appears well in the range of concentration in which H2O2 acts as a global messenger. As a global messenger, H2O2 affects many biological processes beyond oxidative stress response, such as cell metabolism, cellular fate, immunity and homeostasis31, 32. The mechanism shown here for the generation of O2*- and H2O2 by fluorescent proteins explains the cytotoxicity and tissue abnormalities seen in cells and animals expressing fluorescent proteins10-12, 15, 16. Given the extensive use of fluorescent proteins, especially eGFP, in studies of different human diseases (Supplemental Data Table S1) and the known role of O2*- and H2O2 in the pathophysiology of those conditions, we call for a careful evaluation of experimental design, validation and interpretation of results in investigations using fluorescent proteins. AUTHORS' CONTRIBUTIONS

D.G. and F.L. designed experiments, performed experiments and wrote the paper. A.K. performed experiments with bacterial cells and wrote the paper. J.J. supervised the ESR experiments. E.T. supervised the experiments with bacterial cells. C.C.M. and R.M.P. expressed and purified eGFP and RFP. R.P.M. provided guidance and wrote the paper.

COMPETING FINANCIAL INTERESTS

The authors declare no competing financial interests.

MATERIALS & CORRESPONDENCE

Correspondence and materials requests should be addressed to Dr. Ronald P. Mason (mason4@niehs.nih.gov) and Dr. Douglas Ganini da Silva (ganinidasilvad@niehs.nih.gov).

ACKNOWLEDGEMENTS

We acknowledge Dr. Ann Motten and Ms. Mary Mason for editing the manuscript. We thank Dr. Tom Eling and Dr. Thomas Van'T Erve for their suggestions regarding the manuscript. We are also in debt to Ms. Jean Corbett for her invaluable technical assistance. We declare no conflicts of interest. This research was supported by the Intramural Research Program of the NIEHS, National Institute of Environmental Health Sciences/NIH.

REFERENCES

1. Shaner, N.C., Steinbach, P.A. & Tsien, R.Y. A guide to choosing fluorescent proteins. Nature methods 2, 905-909 (2005).

2. Zhang, J., Campbell, R.E., Ting, A.Y. & Tsien, R.Y. Creating new fluorescent probes for cell biology. Nat Rev Mol Cell Biol 3, 906-918 (2002).

3. Zhang, G., Gurtu, V. & Kain, S.R. An enhanced green fluorescent protein allows sensitive detection of gene transfer in mammalian cells. Biochemical and biophysical research communications 227, 707-711 (1996).

4. Merzlyak, E.M. et al. Bright monomeric red fluorescent protein with an extended fluorescence lifetime. Nature methods 4, 555-557 (2007).

5. Zhang, L., Patel, H.N., Lappe, J.W. & Wachter, R.M. Reaction progress of chromophore biogenesis in green fluorescent protein. Journal of the American Chemical Society 128, 4766-4772 (2006).

6. Lu, C., Albano, C.R., Bentley, W.E. & Rao, G. Quantitative and kinetic study of oxidative stress regulons using green fluorescent protein. Biotechnol Bioeng 89, 574-587 (2005).

7. Tsien, R.Y. Constructing and Exploiting the Fluorescent Protein Paintbox (Nobel Lecture). Angewandte Chemie-International Edition 48, 5612-5626 (2009).

8. Ormo, M. et al. Crystal structure of the Aequorea victoria green fluorescent protein. Science 273, 1392-1395 (1996).

9. Subach, O.M. et al. Structural characterization of acylimine-containing blue and red chromophores in mTagBFP and TagRFP fluorescent proteins. Chemistry & biology 17, 333-341 (2010).

10. Taghizadeh, R.R. & Sherley, J.L. CFP and YFP, but not GFP, provide stable fluorescent marking of rat hepatic adult stem cells. Journal of biomedicine & biotechnology 2008, 453590 (2008).

11. Ansari, A.M. et al. Cellular GFP Toxicity and Immunogenicity: Potential Confounders in in Vivo Cell Tracking Experiments. Stem cell reviews 12, 553-559 (2016).

12. Koike, M., Yutoku, Y. & Koike, A. Ku80 attentuates cytotoxicity induced by green fluorescent protein transduction independently of non-homologous end joining. FEBS open bio 3, 46-50 (2013).

13. Liu, H.S., Jan, M.S., Chou, C.K., Chen, P.H. & Ke, N.J. Is green fluorescent protein toxic to the living cells? Biochemical and biophysical research communications 260, 712-717 (1999).

14. Goto, H. et al. Transduction of green fluorescent protein increased oxidative stress and enhanced sensitivity to cytotoxic drugs in neuroblastoma cell lines. Molecular cancer therapeutics 2, 911-917 (2003).

15. Wallace, L.M., Moreo, A., Clark, K.R. & Harper, S.Q. Dose-dependent Toxicity of Humanized Renilla reniformis GFP (hrGFP) Limits Its Utility as a Reporter Gene in Mouse Muscle. Molecular therapy. Nucleic acids 2, e86 (2013).

16. Huang, W.Y., Aramburu, J., Douglas, P.S. & Izumo, S. Transgenic expression of green fluorescence protein can cause dilated cardiomyopathy. Nat Med 6, 482-483 (2000).

17. Wardman, P. Fluorescent and luminescent probes for measurement of oxidative and nitrosative species in cells and tissues: progress, pitfalls, and prospects. Free radical biology & medicine 43, 995-1022 (2007).

18. Wolff, S.P. Ferrous ion oxidation in presence of ferric ion indicator xylenol orange for measurement of hydroperoxides. Methods in enzymology 233, 182-189 (1994).

19. Summers, F.A., Zhao, B., Ganini, D. & Mason, R.P. Photooxidation of Amplex Red to resorufin: implications of exposing the Amplex Red assay to light. Methods in enzymology 526, 1-17 (2013).

20. Subramanian, A. et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proceedings of the National Academy of Sciences of the United States of America 102, 15545-15550 (2005).

21. Duling, D.R. Simulation of multiple isotropic spin-trap EPR spectra. Journal of Magnetic Resonance Series B 104, 105-110 (1994).

22. Koo, M.S. et al. A reducing system of the superoxide sensor SoxR in Escherichia coli. The EMBO journal 22, 2614-2622 (2003).

23. Michan, C., Manchado, M., Dorado, G. & Pueyo, C. In vivo transcription of the Escherichia coli oxyR regulon as a function of growth phase and in response to oxidative stress. Journal of bacteriology 181, 2759-2764 (1999).

24. Dukan, S. & Touati, D. Hypochlorous acid stress in Escherichia coli: resistance, DNA damage, and comparison with hydrogen peroxide stress. Journal of bacteriology 178, 6145-6150 (1996).

25. Dichi, I., Bregano, J.W., Colado Simao, A.a.N. & Cecchini, R. Role of oxidative stress in chronic diseases.

26. Bunn, H.F., Gu, J., Huang, L.E., Park, J.W. & Zhu, H. Erythropoietin: a model system for studying oxygen-dependent gene regulation. The Journal of experimental biology 201, 1197-1201 (1998).

27. Qutub, A.A. & Popel, A.S. Reactive oxygen species regulate hypoxia-inducible factor 1alpha differentially in cancer and ischemia. Molecular and cellular biology 28, 51065119 (2008).

28. Barondeau, D.P., Kassmann, C.J., Tainer, J.A. & Getzoff, E.D. The case of the missing ring: radical cleavage of a carbon-carbon bond and implications for GFP chromophore biosynthesis. Journal of the American Chemical Society 129, 3118-3126 (2007).

29. Niswender, K.D., Blackman, S.M., Rohde, L., Magnuson, M.A. & Piston, D.W. Quantitative Imaging of Green Fluorescent Protein in Cultured-Cells - Comparison of Microscopic Techniques, Use in Fusion Proteins and Detection Limits. Journal of Microscopy-Oxford 180, 109-116 (1995).

30. Tsien, R.Y. The green fluorescent protein. Annu Rev Biochem 67, 509-544 (1998).

31. Stone, J.R. & Yang, S.P. Hydrogen peroxide: A signaling messenger. Antioxidants & Redox Signaling 8, 243-270 (2006).

32. Murrell, G.A., Francis, M.J. & Bromley, L. Modulation of fibroblast proliferation by oxygen free radicals. Biochem J 265, 659-665 (1990).

33. Wood, E.J. Data for biochemical research (third edition) by R M C Dawson, D C Elliott, W H Elliott and K M Jones, pp 580. Oxford Science Publications, OUP, Oxford, 1986. ISBN 0-19-855358-7. Biochemical Education 15, 122 (1987).

34. O'Brien, P.J. Superoxide production. Methods in enzymology 105, 370-378 (1984).

Immature, colorless intermediates of GFP chromophore

R'' OH

Mature GFP chromophore

NAD(P)+ NAD(P)H

Scheme 1. GFP chromophore maturation and proposed mechanism for the catalytical consumption of NAD(P)H and generation of O2'-/H2O2 by GFP immature intermediates.

The GFP chromophore maturation process is shown in black (adapted from Barondeau et al.28). Shown in red, the proposed reaction of the last immature intermediate of GFP with NAD(P)H leads to the catalytical generation of O2*-/H2O2 and consumption of NAD(P)H. The eGFP structure is shown (PDB# 2y0g).

30 min

(a) - eGFP

(a) + SOD

catalase

3440 3460 3480 3500 3520 3540 Magnetic field (G)

Experimental

Simulation

35% DMPO/ 'OH

65% DMPO/ OOH

3440 3460 3480 3500 3520 3540 Magnetic field (G)

N-purged

3460 3480 3500 3520 Magnetic field (G)

1.95-, 1. 1.851.80 1.75 1.70-

N2-purged Re-oxygenated

120 ¡5 <t z

10 20 30 40 50 60 Time (min)

Fig. 1. Reaction of NADH with eGFP leads to the formation of superoxide free radical anion. (a) Spin trapping of superoxide free radical anion (O2'-) generated by eGFP in the presence of NADH. S amples containing 50 pM GFP, 500 pM NADH and 100 mM DMPO showed a multiline-signal, trace (1). The effect of incubation time and eGFP or NADH was also studied. SOD was added to samples at 500 U/mL and catalase at 1 kU/mL. (b) Spectral simulation of trace (1) from panel A. The signal (1) is a superposition of (a) 65% DMPO/•OOH and (b) 35 % DMPO/^OH. (c) ESR spin trapping for O2*- was prepared with samples containing 50 pM eGFP, 500 pM NADH and 100 mM DMPO using air-equilibrated solutions (aerobic) or anaerobic solutions (nitrogen-purged solutions). (d) A sample containing 25 pM eGFP was prepared without oxygen and tested for the ability to consume NADH at 250 pM (sample purged with nitrogen, •). This sample was then aerated with atmospheric air, and NADH consumption was followed (re-oxygenated sample, ▲). NADH was measured by its absorbance at 340 nm (£34onm = 6,220 M-1cm-1 33). Black filled symbols were used to show the absorbance of the fluorescence-active, mature eGFP (• for the nitrogen-purged sample, and ■ for the sample after re-oxygenation, £489nm = 55,000 M-1cm-1).

Figure 2

NADH consumption (pM/min)

000000 o 1 2 3 4 5

O*2 formation (pM/min)

NADH consumption (|M/min) 000000000 bbbbo1111

bl24CT>Wbl24CT>

O*" formation (|M/min)

H2O2 (pM)

>. 6 8 C

Fig. 2. Generation of superoxide free radical anion (O2'~) and H2O2 by eGFP in the presence of NADH. (a) Using eGFP expressed by E. coli, rates of NADH consumption and O2*-formation were determined in samples containing 15 pM eGFP, 150 pM NADH and 20 pM cytochrome c. The effects of catalase (1 kU/mL) and SOD (250 U/mL) were also studied. O2*~ formation was determined by the cytochrome c reduction measured at 550 nm feonm = 21,000 M-1cm-1 34). (b) Rates of NADH consumption and O2*~ formation were studied in samples containing 15 pM GFP, 20 pM cytochrome c and different concentrations of NADH. (c) H2O2 production was determined in samples containing 25 pM GFP and 250 pM NADH. Every 5-6 min, an aliquot was withdrawn and H2O2 measured. Samples to which catalase (1 kU/mL) had been added were assayed in parallel. (d) Using eGFP expressed by the mammalian cell line HEK 293, the rate of O2*~ formation was determined in samples containing 15 pM eGFP, 150 pM NADH and 20 pM cytochrome c, such as described in (a). The effect of SOD (250 U/mL) was also studied. * is used to show statistically significant differences between samples and controls.

Gene expression (fold change)

c . i . i . i . i . i . i . i . i

eGFP " *

katG " *

oxyR " *

gor - *

trxC - h*

grxA " 1*

ahpC "

oxyS " * *

Gene expression (fold change)

O INJ O)

acnA fumC " zwf -soxR sodA soxS "

wwwwwwwwwwW^ $$

^ w^^htwjo

--h ^ o

Gene expression (fold change)

mbp katG-^ oxyR-^1 gor-g trxC-^ grxA-^ ahpC-^ oxyS-^ acnA-^ fumC-^¡

zwf soxR-^1 sodA-^ soxS-^

Fig. 3. Oxidative stress response in E. coli follows the expression of eGFP. Oxidative stress response was determined by the quantification of mRNAs by real-time PCR of genes activated by the SoxR/S and OxyR regulons in E. coli. (a) eGFP expression (fist three bars on the left) and quantification of genes that are targets for the activation of the SoxR/S regulon are shown. (b) Quantification of gene expression for different genes activated by OxyR. (c) Quantification of mRNAs for MBP (maltose binding protein) and targets of SoxR/S and OxyR were determined in an E. coli strain engineered for the controlled expression of high levels of MBP. * is used to show statistically significant differences between samples and controls, and $ shows statistically significant differences between samples of high and low expression.

Figure 4

HeLa HeLa/eGFP

Anti-eGFP Anti-GAPDH

O 20 40 60 80 100 120 Time (min)

HeLa/eGFP vs HeLa

(HeLa.tet.eGFP+DOX - HeLa.tet.eGFP) vs (HeLa.tet+DOX - HeLa.tet)

HALLMARK INTERFERON GAMMA RESPONSE

172 ■ .2292497

HALLMARK INTERFERON GAMMA RESPONSE

172 ■.1225555

HALLMARK TN FA SIGNALING VIA NFKB

180 1.2171624

HALLMARK TNFA SIGNALING VIA NFKB

180 1.1431295

HALLMARK_PROTEIN_SECRETION KL^IARK INFLAMMATORY RESPONSE HALLMARK TGF BETA SIGNALING

83 1.21278 179 1.1497971 51 1.1435704

HALLMARK INTERFERON ALPHA RESPONSE

82 1.1414292

HALLMARK INTERFERON ALPHA RESPONSE

82 1.1135279

HALLMARK_COMPLEMENT HALLMARK KRAS SIGNALING UP

180 0884131 179 -1.0829203

HALLMARK EPITHELIAL MESENCHYMAL TRANSITION 189 -1

HALLMARK EPITHELIAL MESENCHYMAL TRANSITION

Fig. S1D

HALLMARK HEDGEHOG SIGNALING

HALLMARK HEDGEHOG SIGNALING

34 -1.1624191 Fig. S1E

HALLMARK ANDROGEN RESPONSE

HALLMARK ANDROGEN RESPONSE

1 3241451 Fig. S1 F

HALLMARK ANGIOGENESIS

HALLMARK_MYC_TARGETS_V1 HALLMARK_ESTROGEN_RESPONSE_LATE HALLMARK E2F TARGETS

167 183 159

V/S////S,

1.2143934 1.2313852 1.2363398

^488971 _Fic^1^

HALLMARK UV RESPONSE DN

HALLMARK UV RESPONSE DN

HALLMARK_ESTROGEN_RESPONSE_EARLY

HALLMARK_APOPTOSIS

HALLMARK NOTCH SIGNALING

170 -1 155 -1 28 -1

2759069 2941011 2985458

HALLMARK HYPOXIA

HALLMARK HYPOXIA

180 -1.2566458 Fig. S1H

HALLMARK UV RESPONSE UP

3130034 HALLMARK UV RESPONSE UP

148 -1.3328159

Fig Sil

HALLMARK HALLMARK HALLMARK HALLMARK HALLMARK HALLMARK

G2M_CHECKPOINT

GLYCOLYSIS

DNA_REPAIR

MITOTIC_SPINDLE

ADIPOGENESIS

MTORC1 SIGNALING

167 -1 175 -1 135 -1 163 -1 171 -1 173 ■

3177447 3468486 3512396 3532993 3537997 4249561

HALLMARK PANCREAS BETA CELLS

HeLa + HeLa/eGFP -

HIF1a GAPDH

HeLa HeLa/eGFP

Fig. 4. eGFP in cells induces gene expression associated with oxidative stress and other biological processes. Images of HeLa cells and stable eGFP-expressing HeLa cells (HeLa/eGFP) under a fluorescent microscope with a green filter can be seen in (a). eGFP was also detected by Western blotting, with GAPDH staining as the loading control. In (b) extracellular H2O2 was determined in HeLa and HeLa/eGFP. Catalase was added to cultures at 500 U/mL. (c) Microarray-based experiments were designed to show the effect of GFP expression in the gene expression of two independent cell models. The two GSEA (Gene Set Enrichment Analyses20) were performed using the hallmark dataset (h.all.v5.0.symbols) and 1,000 permutations. The contrasts were HeLa/eGFP versus HeLa, and Hela.tet.eGFP + doxycycline subtracted from HeLa.tet.eGFP versus HeLa.tet + doxycycline subtracted from HeLa.tet. The enriched gene sets with p-values lower than 0.001 were ranked, and matches between the two independent analyses are shown highlighted in yellow. The top genes contributing to the enrichment for the commonly enriched gene sets are depicted as heat maps in Supplementary Data Fig. S3. (d) Intracellular levels of HIFIa were determined by Western blotting using whole cell lysates of HeLa cells and HeLa cells stably expressing eGFP. In (e), GAPDH densitometry was used as loading control to normalize the densitometry of HIFIa. * is used to show statistically significant differences between HeLa and HeLa/eGFP.

Highlights

- Immature eGFP and TagRFP catalytically generate O2* - and H2O2 in the presence of NADH.

- Cells expressing eGFP show gene expression associated with oxidative stress.

- Experiments using fluorescent proteins should be alternatively validated.

Cells expressing

NAD(P)H

NAD(P)

Gene expression

adaptation to tH202 and |02'~

Immature eGFP