Difference between revisions of "Part:BBa K4768000"

Line 33: Line 33:
 
</div>
 
</div>
  
<p>The CALIPSO part (BBa_K4768000) consists of the transcriptional repression factor, DhdR, which has been isolated from the bacterium Achromobacter denitrificans. The codon sequence has been optimized for expression in E.coli. Additionally, the presence of a T7 promoter and terminator enables its inducible expression by IPTG in E.coli BL21 (DE3). Finally, a Histidine tag is included in the sequence to facilitate the purification of DhdR.</p>
+
<p>The CALIPSO part (BBa_K4768000) consists of the transcriptional repression factor, DhdR, which has been isolated from the bacterium Achromobacter denitrificans. The codon sequence has been optimized for expression in <i>E.coli</i>. Additionally, the presence of a T7 promoter and terminator enables its inducible expression by IPTG in <i>E.coli</i> BL21 (DE3). Finally, a Histidine tag is included in the sequence to facilitate the purification of DhdR.</p>
<p> This transcriptional repressor was employed in our biosensing system, which was inspired by the work of Ping Xu et al. <a href=”https://www.nature.com/articles/s41467-021-27357-7”>[1]</a> and the work conducted by iGEM Duke 2021. Our objective was to use the affinity between 2-Hydroxyglutarate, an oncometabolite, and DhdR to trigger the production of our drug-activating enzyme when liposomes are situated in a tumoral environment. To achieve this, we positioned our gene of interest under the control of the operator site of DhdR, referred to as dhdO.</p>
+
<p> This transcriptional repressor was employed in the biosensing system, which was inspired by the work of Ping Xu <i>et al</i>. <a href=”https://www.nature.com/articles/s41467-021-27357-7”>[1]</a> and the work conducted by iGEM Duke 2021. The objective was to use the affinity between 2-Hydroxyglutarate, an oncometabolite, and DhdR to trigger the production of the drug-activating enzyme when liposomes are situated in a tumoral environment. To achieve this, we positioned the gene of interest under the control of the operator site of DhdR, referred to as dhdO.</p>
  
 
<div align="center">
 
<div align="center">
Line 41: Line 41:
 
             style="width:60%;"
 
             style="width:60%;"
 
             src="https://static.igem.wiki/teams/4768/wiki/module-2/unnamed-2.png">  
 
             src="https://static.igem.wiki/teams/4768/wiki/module-2/unnamed-2.png">  
             <figcaption class="normal"><span class="titre-image"><i><b>Figure 2: Operating principle of our biosensor.</b> In presence of 2-HG, repression is removed and the gene is expressed leading to the production of the thymidine phosphorylase which converts the Tegafur prodrug into an active anticancer drug, 5-fluorouracile.
+
             <figcaption class="normal"><span class="titre-image"><i><b>Figure 2: Operating principle of the biosensor.</b> In presence of 2-HG, repression is removed and the gene is expressed leading to the production of the thymidine phosphorylase which converts the Tegafur prodrug into an active anticancer drug, 5-fluorouracile.
 
</i></span></figcaption>
 
</i></span></figcaption>
 
         </figure>
 
         </figure>
Line 97: Line 97:
 
<h3>2) Functionality tests in PURE system </h3>
 
<h3>2) Functionality tests in PURE system </h3>
  
<p>The aim of our experiments was to establish that the binding of the repressor DhdR to its operator site, <i>dhdO</i>, effectively inhibits transcription of a gene of interest regulated by <i>dhdO</i>. Then, we wanted to show that the presence of 2-HG leads to the de-repression of that gene in PURE system.</p>
+
<p>The aim of the experiments was to establish that the binding of the repressor DhdR to its operator site, <i>dhdO</i>, effectively inhibits transcription of a gene of interest regulated by <i>dhdO</i>. Then, we wanted to show that the presence of 2-HG leads to the de-repression of that gene in PURE system.</p>
  
<p>Inhibition was tested on the <i>sfgfp</i> reporter gene by fluorescence measurements. To determine the minimal concentration of DhdR required to obtain strong repression, sfGFP was synthesized in the presence of different concentrations of DhdR. Our biochemical network model predicted a range of DhdR concentrations expected to lead to different sfGFP levels, which we experimentally tested.</p>
+
<p>Inhibition was tested on the <i>sfgfp</i> reporter gene by fluorescence measurements. To determine the minimal concentration of DhdR required to obtain strong repression, sfGFP was synthesized in the presence of different concentrations of DhdR. The biochemical network model predicted a range of DhdR concentrations expected to lead to different sfGFP levels, which we experimentally tested.</p>
  
 
<div align="center">
 
<div align="center">
Line 193: Line 193:
 
<p>The impact of DhdR on cell-free expression was modeled through a set of differential equations inspired from iGEM Teams Duke 2021 and Delft 2021, and implemented on  <a href="https://gitlab.igem.org/2023/toulouse-insa-ups/-/tree/main/CALIPSO_model/Pre-calibration?ref_type=heads">COPAS v(4.38)I.</a></p>
 
<p>The impact of DhdR on cell-free expression was modeled through a set of differential equations inspired from iGEM Teams Duke 2021 and Delft 2021, and implemented on  <a href="https://gitlab.igem.org/2023/toulouse-insa-ups/-/tree/main/CALIPSO_model/Pre-calibration?ref_type=heads">COPAS v(4.38)I.</a></p>
  
<p>The responsiveness of our anticancer liposome to a tumoral environment is a key design feature for limiting toxic side-effects. We thus used the model to evaluate whether 5-FU was produced at higher concentration in the presence of 2-HG (synthesized by cancer cells) compared to its absence. In absence of 2-HG, we expect the expression of TYPH to be repressed and, therefore, only a small amount of 5-FU would be produced. Consistently, simulations shown in Fig. 3 predicted that 1.5 μM of DhdR repressed gene expression, resulting in a ~2-fold lower concentration of 5-FU outside the liposome. Adding 100 μM of 2-HG, corresponding to typical concentrations in surrounding tumors <a href=”https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7533892/”>[2]</a>, derepressed the system and boosted 5-FU production. The model therefore met the behavior expected for our liposomes and confirmed that our targeted drug production strategy could help tackle adverse effects of current chemotherapies.</p>
+
<p>The responsiveness of the anticancer liposome to a tumoral environment is a key design feature for limiting toxic side-effects. We thus used the model to evaluate whether 5-FU was produced at higher concentration in the presence of 2-HG (synthesized by cancer cells) compared to its absence. In absence of 2-HG, we expect the expression of TYPH to be repressed and, therefore, only a small amount of 5-FU would be produced. Consistently, simulations shown in Fig. 3 predicted that 1.5 μM of DhdR repressed gene expression, resulting in a ~2-fold lower concentration of 5-FU outside the liposome. Adding 100 μM of 2-HG, corresponding to typical concentrations in surrounding tumors <a href=”https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7533892/”>[2]</a>, derepressed the system and boosted 5-FU production. The model therefore met the behavior expected for liposomes and confirmed that the targeted drug production strategy could help tackle adverse effects of current chemotherapies.</p>
  
 
<div align="center">
 
<div align="center">
Line 207: Line 207:
 
<h3> 2) Guiding experiments and validating the model </h3>
 
<h3> 2) Guiding experiments and validating the model </h3>
  
<p>We used the model to determine the optimal concentration of DhdR to be encapsulated in our liposome. To that aim, we built a reduced model containing only the biosensing and gene expression modules, along with a gene encoding the sfGFP protein that was used in our experiments as a reporter. To verify model predictions, we carried out the same experiments as done in silico using the same concentration range of DhdR and 2-HG. </p>
+
<p>We used the model to determine the optimal concentration of DhdR to be encapsulated in the liposome. To that aim, we built a reduced model containing only the biosensing and gene expression modules, along with a gene encoding the sfGFP protein that was used in the experiments as a reporter. To verify model predictions, we carried out the same experiments as done in silico using the same concentration range of DhdR and 2-HG. </p>
  
<p>A gradual decrease in the fluorescence intensity, i.e., in the expression of <i>sfgfp</i>, was experimentally verified when increasing the concentration of DhdR. Importantly, the measurements are consistent with the model predictions qualitatively, confirming the parameter values and structure of the model.</p>
+
<p>A gradual decrease in the fluorescence intensity, i.e., in the expression of <i>sfgfp</i>, was experimentally verified when increasing the concentration of DhdR. Importantly, measurements are consistent with the model predictions qualitatively, confirming the parameter values and structure of the model.</p>
  
 
<p>Experimental results confirm that increasing 2-HG concentration progressively lifts the repression caused by DhdR. However, the model poorly predicted the gradual increase of sfGFP levels measured experimentally, as it was more sensitive to lower concentrations of 2HG.</p>
 
<p>Experimental results confirm that increasing 2-HG concentration progressively lifts the repression caused by DhdR. However, the model poorly predicted the gradual increase of sfGFP levels measured experimentally, as it was more sensitive to lower concentrations of 2HG.</p>
  
<p>We thus decided to use our measurements to refine the model. We estimated the values of different parameters of the biosensing module (k<sub>D2</sub>, k<sub>-D2</sub>, k<sub>D0</sub>, k<sub>-D0</sub>, k<sub>DH</sub> and k<sub>-DH</sub>) by fitting our experimental data. In comparison to the predictions of the <i>ab initio</i> model constructed from the literature, the simulations with the refined parameter values were in quantitative agreement with our measurements. Importantly, the refined model with the calibrated biosensor module was also able to accurately simulate protein expression levels in the presence of 2-HG, corroborating its improved predictive power (Fig. 11).</p>
+
<p>We thus decided to use measurements to refine the model. We estimated the values of different parameters of the biosensing module (k<sub>D2</sub>, k<sub>-D2</sub>, k<sub>D0</sub>, k<sub>-D0</sub>, k<sub>DH</sub> and k<sub>-DH</sub>) by fitting the experimental data. In comparison to the predictions of the <i>ab initio</i> model constructed from the literature, the simulations with the refined parameter values were in quantitative agreement with measurements. Importantly, the refined model with the calibrated biosensor module was also able to accurately simulate protein expression levels in the presence of 2-HG, corroborating its improved predictive power (Fig. 11).</p>
  
 
<div align="center">
 
<div align="center">
Line 227: Line 227:
  
 
<h2>Conclusion and Perspectives</h2>
 
<h2>Conclusion and Perspectives</h2>
<p>These experiments provide evidence that our 2-HG biosensor that relies on the part BBa_K4768000 is functional in bulk reactions at concentrations of oncometabolite that are physiologically relevant. Our biochemical network model was used to predict DhdR concentrations, and then optimized according to experimental results. Moreover,we established a protocol for encapsulating DhdR, the expression product of the part BBa_K4768000, inside liposomes. Although DhdR-based repression was not clearly demonstrated in vesicles, we gave recommendations and provided image analysis tools for future investigations.</p>
+
<p>These experiments provide evidence that the 2-HG biosensor that relies on the part BBa_K4768000 is functional in bulk reactions at concentrations of oncometabolite that are physiologically relevant. The biochemical network model was used to predict DhdR concentrations, and then optimized according to experimental results. Moreover,we established a protocol for encapsulating DhdR, the expression product of the part BBa_K4768000, inside liposomes. Although DhdR-based repression was not clearly demonstrated in vesicles, we gave recommendations and provided image analysis tools for future investigations.</p>
  
 
<h2>References</h2>
 
<h2>References</h2>

Revision as of 10:44, 9 October 2023


DhdR repressor

Part for expression and purification of our transcriptional repressor DhdR for biosensing.


Sequence and Features


Assembly Compatibility:
  • 10
    INCOMPATIBLE WITH RFC[10]
    Illegal XbaI site found at 47
  • 12
    COMPATIBLE WITH RFC[12]
  • 21
    INCOMPATIBLE WITH RFC[21]
    Illegal BglII site found at 493
    Illegal BglII site found at 622
  • 23
    INCOMPATIBLE WITH RFC[23]
    Illegal XbaI site found at 47
  • 25
    INCOMPATIBLE WITH RFC[25]
    Illegal XbaI site found at 47
  • 1000
    INCOMPATIBLE WITH RFC[1000]
    Illegal BsaI.rc site found at 756


Introduction

Figure 1:DhdR part

The CALIPSO part (BBa_K4768000) consists of the transcriptional repression factor, DhdR, which has been isolated from the bacterium Achromobacter denitrificans. The codon sequence has been optimized for expression in E.coli. Additionally, the presence of a T7 promoter and terminator enables its inducible expression by IPTG in E.coli BL21 (DE3). Finally, a Histidine tag is included in the sequence to facilitate the purification of DhdR.

This transcriptional repressor was employed in the biosensing system, which was inspired by the work of Ping Xu et al. [1] and the work conducted by iGEM Duke 2021. The objective was to use the affinity between 2-Hydroxyglutarate, an oncometabolite, and DhdR to trigger the production of the drug-activating enzyme when liposomes are situated in a tumoral environment. To achieve this, we positioned the gene of interest under the control of the operator site of DhdR, referred to as dhdO.

Figure 2: Operating principle of the biosensor. In presence of 2-HG, repression is removed and the gene is expressed leading to the production of the thymidine phosphorylase which converts the Tegafur prodrug into an active anticancer drug, 5-fluorouracile.

Construction

The CALIPSO part BBa_K4768000 comprises the transcriptional repression factor DhdR fused with a histidine tag at its N-terminus. The synthesis of this gBlock was performed and provided by IDT.

The gBlock was then cloned into the pET_21a(+) plasmid and transformed into Stellar competent cells.

Primers used to clone this part in the pET21: (from 5' to 3'):
  • DhdR-pET21-F: AGCAGCCGGATCTCATCATGACGTCTGACGCGC
  • DhdR-pET21-R: GAAGGAGATATACATATGGGCCATCATCATCATCATC

Figure 3 shows the enzymatic restriction pattern of the resulting clones. Clone 4 was digested using EcoRV and NdeI. Two bands were expected at 1.3 kb and 4.8 kb, as experimentally measured (lane 5).

Figure 3: Enzymatic digestion of the plasmid part extracted from clone 4. Patterns from simple digestion by EcoRV or by NdeI, and double digestion by both enzymes are shown.

Clone 4 was sequence verified.

Characterization

1) Production and purification of DhdR

The pET21a(+) vector including the dhdR insert was transformed into E.coli strain BL21 (DE3). This strain was provided by Cédric Montanier (researcher at TBI). When DO reaches 0.5-0.6, expression of the recombinant protein was induced overnight at 16°C using IPTG. The His-tagged protein was then purified on TALON® Metal Affinity Resin. Pure fractions were assessed by SDS-PAGE.Results are shown in Figure 4.

Figure 4: SDS-PAGE analysis of the different steps of the DhdR purification: protein ladder (SE250 Mighty Small II Mini Vertical Protein Electrophoresis Unit), pellet DhdR, cell-free extract (CFE), flowthrough (FT), wash (W), elution with 50 mM imidazole, 100 mM imidazole, 250 mM imidazole and 500 mM imidazole (respectively E150, E1100, E1250 and E1500). The band corresponding to the protein was clearly identified by Coomassie Blue staining.

The expected size for DhdR is 28.27 kDa. Clear bands were observed for the four elution fractions E150, E1100, E1250 and E1 500, but not in the negative control sample, as expected. These data show that DhdR was efficiently purified and can be used for subsequent assays.

The E1250 fraction was dialysed leading to a concentration of 14.4 µM (> 95% pure protein). Fractions E1100 and E1500 were pooled and dialysed, resulting in a concentration of 7.63 µM (> 95% pure protein).

2) Functionality tests in PURE system

The aim of the experiments was to establish that the binding of the repressor DhdR to its operator site, dhdO, effectively inhibits transcription of a gene of interest regulated by dhdO. Then, we wanted to show that the presence of 2-HG leads to the de-repression of that gene in PURE system.

Inhibition was tested on the sfgfp reporter gene by fluorescence measurements. To determine the minimal concentration of DhdR required to obtain strong repression, sfGFP was synthesized in the presence of different concentrations of DhdR. The biochemical network model predicted a range of DhdR concentrations expected to lead to different sfGFP levels, which we experimentally tested.

Figure 5: Effect of different concentrations of DhdR on the expression of a fluorescent reporter gene. Experiments 1 and 2 were performed with the same batch of PCR product from clone 8, while a new batch of PCR product from the same clone was used in experiment 3. PUREfrex2.0 was used in all conditions. The intensity value of sfGFP without DhdR was used for normalization. Excitation and emission wavelengths were 488 nm and 510 nm, respectively.

As expected, the higher the concentration of DhdR, the stronger the repression in all three experiments (Figure 5). With the new batch of linear DNA, repression was consistently stronger. We deduced from these results that the optimal concentration of DhdR to efficiently repress expression of a gene under transcriptional control of a dhdO operator sequence was 1.5 µM, validating the predictions of the biochemical network model.

Induction of gene expression that was repressed by 1.5 µM of DhdR was then assayed using physiological concentrations of 2-HG found around tumor cells, i.e., between 10 and 100 µM. A higher concentration was also tested, corresponding to full saturation of the DhdR repressor. The results demonstrate that 2-HG de-represses transcription of DhdR-bound DNA in a concentration-dependent manner (Figure 6). Up to 48% of sfGFP signal was recovered at a saturating concentration of 2-HG. The reason why protein production is not fully restored remains to be investigated.

Figure 6: Effect of different concentrations of 2-HG on DhdR repression. PUREfrex2.1 was used. The intensity value of sfGFP without DhdR and with 10 µM 2-HG was used for normalization. Each concentration was corrected taking into account the inactivation effect of 2-HG on PURE system.

3) In-liposome expression of the sfgfp gene

  • DhdR repression in liposomes
  • Two experiments were conducted in which the PURE system solution was encapsulated in liposomes along with the sfgfp gene, either with 1.5 µM of DhdR or without DhdR. Liposomes were imaged by optical microscopy.

    Figure 7a displays a population of liposomes localized by the membrane dye Topfluor594. A zoom-in image of liposomes showed the fluorescent rim characteristic of membrane-labeled vesicles (Fig. 7b). The line intensity profile generated with ImageJ confirmed that the intensity was higher at the membrane and lower inside the liposome (Fig. 7c).

    Figures 7: a) Liposomes were localized with a fluorescent membrane dye. b) Zoom-in image of the liposome area depicted with a red arrow in a. A yellow line crossing the liposome has been appended. c) Fluorescence intensity profile along the line appended in b. The two peaks correspond to the two regions of membrane crossing.

    Figure 8a displays a population of liposomes expressing the sfGFP gene. In the liposome shown in Fig. 8b, one can clearly see the distribution of GFP fluorescence inside the lumen of the liposome. A quantitative analysis is represented in Fig. 8c. Analysis of the two samples with or without DhdR did not reveal notable differences neither in the occurrence of liposomes exhibiting GFP nor in the intensity level of GFP inside individual liposomes.

    Figures 8: Fluorescence microscopy image of liposomes in the GFP channel. Expressed sfGFP signal was localized in the liposome lumen. a) Large field-of-view. b) Blow-up of the liposome depicted with the red arrow in a. c) Intensity profile along the yellow line appended in b.
  • Liposomes are capable of expressing GFP in the presence of living cancer cells
  • Two experiments were conducted in which the PURE system solution was encapsulated in liposomes along with the sfgfp gene, Tegafur and 1.5 µM of DhdR. Liposomes were also coated with anti-HER2-nb and folate ligands. Two conditions were tested for exposing liposomes to Caco-2 cancer cells. In the first protocol, liposomes were incubated in a thermocycler for gene expression prior to their functionalization with anti-HER2-nb and injection on top of cancer cells (sample 1). In a second protocol, liposomes pre-coated with anti-HER2-nb were injected in the growth medium on top of Caco-2 cells, where they have been incubated for in situ gene expression (sample 2). The latter protocol more closely mimics the in vivo conditions for drug delivery. Fluorescence microscopy was used to image living cells, liposomes and sfGFP expression.

    In sample 1, in the field of view displayed in Figures 9, two liposomes were localized using the fluorescent membrane dye (Fig. 9b) but only one exhibits sfGFP signal (Fig. 9c).

    Figures 9: Fluorescence microscopy images of anti-HER2-nb- and folate-decorated liposomes on top of Caco-2 cells (sample 1). a) Imaging of Caco-2 cells in the Brightfield channel. b) Imaging of liposomes localized with a fluorescent membrane dye. c) Imaging of liposomes in the GFP channel.

    Similar results were obtained with sample 2, as shown in Figures 10.

    Figures 10: Fluorescence microscopy images of anti-HER2-nb- and folate-decorated liposomes on top of Caco-2 cells (sample 2).a) Imaging of Caco-2 cells in the Brightfield channel. b) Imaging of a liposome localized with a fluorescent membrane dye. c) Imaging of the same liposome in the GFP channel.

    No difference was observed between liposomes incubated at 37°C and those incubated directly on cancer cells. In both cases, we obtained some liposomes able to produce sfGFP. Follow-up experiments will be necessary to ascertain that gene expression was enabled by 2-HG and not by insufficient repression by DhdR. For instance, optimizing the relative and absolute amounts of DNA and DhdR in liposomes will allow for a better discrimination between repressing and non-repressing conditions.

Molecular Modeling

1) Model-driven feasibility assessment of the 2-HG biosensor

The impact of DhdR on cell-free expression was modeled through a set of differential equations inspired from iGEM Teams Duke 2021 and Delft 2021, and implemented on COPAS v(4.38)I.

The responsiveness of the anticancer liposome to a tumoral environment is a key design feature for limiting toxic side-effects. We thus used the model to evaluate whether 5-FU was produced at higher concentration in the presence of 2-HG (synthesized by cancer cells) compared to its absence. In absence of 2-HG, we expect the expression of TYPH to be repressed and, therefore, only a small amount of 5-FU would be produced. Consistently, simulations shown in Fig. 3 predicted that 1.5 μM of DhdR repressed gene expression, resulting in a ~2-fold lower concentration of 5-FU outside the liposome. Adding 100 μM of 2-HG, corresponding to typical concentrations in surrounding tumors [2], derepressed the system and boosted 5-FU production. The model therefore met the behavior expected for liposomes and confirmed that the targeted drug production strategy could help tackle adverse effects of current chemotherapies.

Figure 11: Concentration of 5-FU after 6 hours in presence or absence of DhdR and 2-HG (concentrations are as indicated), with an initial Tegafur concentration of 1000 μM.

2) Guiding experiments and validating the model

We used the model to determine the optimal concentration of DhdR to be encapsulated in the liposome. To that aim, we built a reduced model containing only the biosensing and gene expression modules, along with a gene encoding the sfGFP protein that was used in the experiments as a reporter. To verify model predictions, we carried out the same experiments as done in silico using the same concentration range of DhdR and 2-HG.

A gradual decrease in the fluorescence intensity, i.e., in the expression of sfgfp, was experimentally verified when increasing the concentration of DhdR. Importantly, measurements are consistent with the model predictions qualitatively, confirming the parameter values and structure of the model.

Experimental results confirm that increasing 2-HG concentration progressively lifts the repression caused by DhdR. However, the model poorly predicted the gradual increase of sfGFP levels measured experimentally, as it was more sensitive to lower concentrations of 2HG.

We thus decided to use measurements to refine the model. We estimated the values of different parameters of the biosensing module (kD2, k-D2, kD0, k-D0, kDH and k-DH) by fitting the experimental data. In comparison to the predictions of the ab initio model constructed from the literature, the simulations with the refined parameter values were in quantitative agreement with measurements. Importantly, the refined model with the calibrated biosensor module was also able to accurately simulate protein expression levels in the presence of 2-HG, corroborating its improved predictive power (Fig. 11).

Figures 12: Comparison of the predictive capabilities of the kinetic model before (in pink) and after calibration (in green) using experimental data as inputs (in red). (A) Concentrations of DhdR were as indicated and 2-HG was not present. (B) Concentration of DhdR was fixed to 1.5 µM and different concentrations of 2-HG were used as indicated.

Conclusion and Perspectives

These experiments provide evidence that the 2-HG biosensor that relies on the part BBa_K4768000 is functional in bulk reactions at concentrations of oncometabolite that are physiologically relevant. The biochemical network model was used to predict DhdR concentrations, and then optimized according to experimental results. Moreover,we established a protocol for encapsulating DhdR, the expression product of the part BBa_K4768000, inside liposomes. Although DhdR-based repression was not clearly demonstrated in vesicles, we gave recommendations and provided image analysis tools for future investigations.

References

  1. [1]Xiao, D., Zhang, W., Guo, W., Liu, Y., Hu, C., Guo, S., Kang, Z., Xu, X., Ma, C., Gao, C., & Xu, P. 2021. A D-2-hydroxyglutarate biosensor based on specific transcriptional regulator DhdR. Nature Communications 12, 7108.
  2. [2]Jezek, P. 2020. 2-Hydroxyglutarate in Cancer Cells. Antioxid Redox Signal, 33(13),903-926.