Reversing protonation of weakly basic drugs greatly enhances intracellular diffusion and decreases lysosomal sequestration

  1. Debabrata Dey
  2. Shir Marciano
  3. Anna Poryval
  4. Ondřej Groborz
  5. Lucie Wohlrabova
  6. Tomás Slanina
  7. Gideon Schreiber  Is a corresponding author
  1. Department of Biomolecular Sciences, Weizmann Institute of Science, Israel
  2. Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Czech Republic
25 figures, 3 tables and 2 additional files

Figures

Chemical structures and pKa values (prediction from ChemAxon software) of the small molecules and drugs used in this study.
Schematic representation of Line-FRAP and diffusion coefficients of proteins in the HeLa cells.

(A) Line-FRAP performed in a fluorescent drug labeled HeLa cell. The white line represents the scanning line; The red circular area denotes the bleaching area. (B) A single scanning line as a function of time (in the vertical direction) including the bleach is shown. The white vertical rectangular region of interest (ROI), which was used for analysis. A close-up view of the Line-FRAP profile, is also shown where the horizontal rectangular area (white colored marks ROI) denotes the fluorescence intensity as a function of length across the scanning line. (C) Average of 30 FRAP recovery curves as a function of time (N=30; Correlation coefficient of exponential fit: R=0.98) and (D) Averaged bleach size profile with gaussian fits (N=30; R=0.99) are shown. (E) Comparison of diffusion coefficients (Dconfocal in blue circles) and percentage of recoveries (in red squares) of bacterial proteins and BSA as measured in HeLa cell cytoplasm are shown. Error bars represent SE calculated from fitting the FRAP progression curves, which are averaged over at least 30 independent measurements.

Figure 3 with 3 supplements
Fluorescein, CCF2 and CF514 diffusion in PBS and inside HeLa cells.

Comparative merged FRAP profiles (N=30; R=0.99 for each of the fits) in PBS and inside HeLa cells with exponential fits are shown for (A) Fluorescein, (E) CCF2 and (H) CF514. Comparative Dconfocal values for (C) Fluorescein (F) CCF2 and (I) CF514 are also shown. (B) Dconfocal values as calculated from individual FRAP curves for Fluorescein (see also Figure 3—figure supplement 2). (C) Dconfocal values as calculated from the merged FRAP curves in panel A is shown (D) HeLa cells after Micro-injection of Fluorescein; (G) CCF2; and (J) for CF514. Error bars represent SE calculated from fitting the FRAP progression curves of 30 merged independent measurements. Statistical significance calculations are detailed in the Materials and methods section.

Figure 3—figure supplement 1
FRAP measurements for Fluorescein following different incubation times.

(A) Micrographs of HeLa cells instantaneously after micro-injection vs (B) FRAP after 24 hr incubation. (C) Comparative merged FRAP profiles with fits (N=30; R=0.99 for each of the fits) and (D) their Dconfocal values. Error bars represent SE calculated from fitting the FRAP progression curves averaged over at least 30 independent measurements.

Figure 3—figure supplement 2
Dconfocal values calculated from individual FRAP and merged FRAP curves and diffusion rates and percent recovery as a function of Fluorescein concentration inside HeLa Cells.

FRAP curves and Dconfocal values for Fluorescein (A–D) inside HeLa cell cytoplasm are shown. Three representative individual FRAP profiles with exponential fits with R=0.80–0.90 are shown in (A). Individual Dconfocal values and their geometric means from individual FRAP curves are shown for Fluorescein in (B). (C) Shows the merged FRAP curves over 30 independent measurements. This provides an improved fit (R=0.99). Panel (D) show Dconfocal values calculated from the merged FRAP profiles. Error bars represent SE calculated from fitting the FRAP progression curves. The diffusion rates and percent recoveries in (E–H) were calculated from individual FRAP progression curves. The drug concentration is linearly related to the fluorescence intensity. (E) show FRAP rates as a function of the fluorescence signal for Fluorescein in the cytoplasm and (F) in the nucleus. (G–H) are as (E–F) but the y-axis denotes percent recovery after FRAP.

Figure 3—figure supplement 3
UV spectrum of CF514 Dye on pH variations.

Comparative UV spectrum of CF514 labelling dye in 20 mM sodium phosphate buffer at different pH values are shown in (A–B), from which a pKa value of 3 was calculated.

Figure 4 with 2 supplements
Diffusion of GSK3 inhibitor in PBS and inside HeLa cells.

Comparative (A) averaged FRAP recovery profiles with exponential fits for GSK3 inhibitor (SB216763; N=30; R=0.99 for each of the fits). Dots show prebleach fluorescence, Black line for DMEM, blue for PBS buffer +50 mg/mL of BSA, green for HeLa cells and orange for HeLa extract. (B) Dconfocal values calculated from individual FRAP curves. (C) Dconfocal values calculated from merger of 30 independent FRAP curves. (D–E) Colocalization of GSK3 inhibitor in lysosomal compartments. HeLa cells were treated with Lyso Tracker 633 for 30 min. GSK3 inhibitor treatment for (D) 45 mins and (E) 2 hr respectively at 10 µM concentrations are also shown. Error bars represent SE calculated from fitting the FRAP progression curves, which are averaged over at least 30 independent measurements. Statistical significance calculations are detailed in the Materials and methods section.

Figure 4—figure supplement 1
Diffusion rates at different doses of GSK3 inhibitor in HeLa cytoplasm.

Effect of drug dosage (2–10 µM) on (A) averaged FRAP profiles with fits (N=30; R=0.99 for each of the fits) and (B) Dconfocal values. Error bars represent SE calculated from fitting the FRAP progression curves, which are averaged over at least 30 independent measurements.

Figure 4—figure supplement 2
Dconfocal values calculated from individual FRAP and merged FRAP curves and diffusion rates and percent recovery as a function of GSK3 inhibitor concentration inside HeLa cells.

FRAP curves and Dconfocal values for GSK3 inhibitor (A–D) inside HeLa cell cytoplasm are shown. Three representative individual FRAP profiles with exponential fits with R=0.90–0.95 are shown in (A). Individual Dconfocal values and their geometric means from individual FRAP curves are shown for GSK3 inhib. in (B). (C) Shows the merged FRAP curves over 30 independent measurements. This provides an improved fit (R=0.99). Panel (D) show Dconfocal values calculated from the merged FRAP profiles. Error bars represent SE calculated from fitting the FRAP progression curves. The diffusion rates and percent of recoveries in (E–F) were calculated from individual FRAP progression curves. The drug concentration is linearly related to the fluorescence intensity. (E) show FRAP rates as a function of the fluorescence signal for GSK3 inhibitor in the cytoplasm. (F) is as (E) but the y-axis denotes percent recovery after FRAP.

Figure 5 with 3 supplements
Diffusion of Quinacrine in PBS and in HeLa cells.

Comparative (A) averaged FRAP recovery profiles with exponential fits for Quinacrine dihydrochloride (N=30; R=0.99 for each of the fits), (B) Dconfocal values calculated from individual FRAP curves. (C) Dconfocal values calculated from merger of 30 independent FRAP curves measured in PBS buffer, 50 mg/mL of BSA, HeLa cell extract and in live HeLa cells. (D–E) Colocalization of Quinacrine DHC in lysosomal compartments are shown. HeLa cells were treated with Lyso Tracker 633 for 30 min. Quinacrine DHC treatment was done for (D) 45 min and (E) 2 hr respectively at 10 µM concentration. Error bars represent SE calculated from fitting the FRAP progression curves, which are averaged over at least 30 independent measurements. Statistical significance calculations are detailed in the Materials and methods section.

Figure 5—figure supplement 1
Dconfocal values calculated from individual FRAP and merged FRAP curves and diffusion rates and percent recovery as a function of Quinacrine concentration inside HeLa cells.

FRAP curves and Dconfocal values for Quinacrine (A–D) inside HeLa cell cytoplasm are shown. Three representative individual FRAP profiles with exponential fits with R=0.90–0.95 are shown in (A). Individual Dconfocal values and their geometric means from individual FRAP curves are shown for Quinacrine in (B). (C) Shows the merged FRAP curves over 30 independent measurements. This provides an improved fit (R=0.99). Panel (D) shows Dconfocal values calculated from the merged FRAP profiles. Error bars represent SE calculated from fitting the FRAP progression curves. The diffusion rates and percent of recoveries in (E–F) were calculated from individual FRAP progression curves. The drug concentration is linearly related to the fluorescence intensity. (E) shows FRAP rates as a function of the fluorescence signal for Quinacrine in the cytoplasm. (F) is as (E) but the y-axis denotes percent recovery after FRAP.

Figure 5—figure supplement 2
Quinacrine diffusion in HeLa after micro-injection or incubation.

Comparison of FRAP immediately following micro-injection (instant FRAP) and FRAP following 2–24 hr incubation. HeLa cell micrographs following (A) micro-injection and (B) 2 hr of incubation. Comparative (C) averaged FRAP profiles with exponential curve fits (N=30; R=0.99 for each of the fits) and (D) calculated Dconfocal values. Error bars represent SE calculated from fitting the FRAP progression curves, which are averaged over at least 30 independent measurements.

Figure 5—figure supplement 3
Dose and time dependency of Quinacrine diffusion in HeLa cells.

HeLa cell micrographs after 2 hr of incubation with (A) 2 µM and (B) 6 µM of Quinacrine. (C) and (D) are micrographs of HeLa cells after 24 hr of incubation with 2 µM and 6 µM Quinacrine, respectively. (E) Effect of drug dose and time dependency on comparative FRAP profiles, with fits (N=30; R=0.99 for each of the fits) and (F) Calculated Dconfocal values from the fits in (E). Error bars represent SE calculated from fitting the FRAP progression curves, which are averaged over at least 30 independent measurements.

Figure 6 with 1 supplement
Mitoxantrone and Primaquine diffusion inside HeLa cells.

Averaged XY- FRAP recovery profiles with exponential fits (N=20; R=0.98 for each of the fits) for (A) Mitoxantrone and (E) Primaquine in HeLa cells. Fluorescent and Transmission channel images of HeLa cells after 45 min incubations with (B–C) Mitoxantrone and (F–G) Primaquine are shown. Time lapse micrographs of a portion of HeLa cell using a classical rectangular XY-FRAP protocol for (D) Mitoxatrone and, (H) Primaquine are shown.

Figure 6—figure supplement 1
Colocalization of Primaquine and Amidoqiuine with Lyso-tracker.

Colocalization of Primaquine (A–B) and Amidoquine (C–D) with Lyso Tracker 633 in HeLa cells. Time-dependent measurements show higher aggregation of both drugs in the lysosomes over time.

Figure 7 with 2 supplements
Effect of sodium azide on the diffusion of GSK3 inhibitor and Quinacrine inside HeLa cells.

(A) GSK3 inhibitor and (E) Quinacrine-dihydrochloride-treated HeLa cells with or without sodium azide. Comparison of merged FRAP profiles with exponential fits (N=30; R=0.99 for each of the fits) and calculated Dconfocal values from individual fits and combined averaged fit for (B–D) GSK3 inhibitor and (F–H) Quinacrine dihydrochloride are shown. Dconfocal estimations from individual FRAP curves are shown in (C,G). Error bars represent SE calculated from fitting the FRAP progression curves, which are averaged over at least 30 independent measurements. Statistical significance calculations are detailed in the Materials and methods section.

Figure 7—figure supplement 1
Inhibition of lysosomal accumulation of GSK3 inhibitor by Bafilomycin A1.

(A–B) Colocalization of GSK3 inhibitor with Lyso Tracker 633 in HeLa cells. (A) Untreated and (B) treated with 100 nM Bafilomycin A1 and corresponding colocalizations are shown. (C) Comparative averaged FRAP profiles with fits (N=30; R=0.99 for each of the fits) and (D) DconfocalDconfocal diffusion values with and without the Bafilomycin A1 treatment. Error bars represent SE calculated from fitting the FRAP progression curves, which are averaged over at least 30 independent measurements.

Figure 7—figure supplement 2
Inhibition of lysosomal accumulation of Quinacrine by Bafilomycin A1.

(A–B) Colocalization of Quinacrine with Lyso Tracker 633 in HeLa cells. HeLa cells (A) untreated and (B) treated with 100 nM Bafilomycin A1 and corresponding colocalizations. (C) Comparative averaged FRAP profiles with fits (N=30; R=0.99 for each of the fits) and (D) Dconfocal values with and without the Bafilomycin A1 treatment. Error bars represent SE calculated from fitting the FRAP progression curves, which are averaged over at least 30 independent measurements.

Figure 8 with 2 supplements
FRAP of Fluorescein analogues inside HeLa cells.

(A) Structures of Fluorescein analogues with pKa values. (B) FRAP recoveries with exponential fits for Fluorescein, CCF2 and 5-amino Fluorescein, with/without Na-Azide treatments in HeLa cells (N=30; R=0.99 for each of the fits). (C) Line FRAP profiles with time lapses. (D) Micrograph images of treated HeLa cells after small molecule incubations.

Figure 8—figure supplement 1
Comparative FRAP of Primaquine analogues in HeLa cell cytoplasm.

(A) Structures of Primaquine analogues. Comparative (B–D) XY FRAP micrographs with time lapses in HeLa cell. (E–H) Micrographs of treated HeLa cells either with micro-injection or after drug incubations with Prim or Prim-Ac. Comparative (I) XY-FRAP recovery profiles with fits (N=20; R=0.98 for each of the fits) and (J) percentage of recoveries for Primaquine analogues in HeLa cytoplasm (N=20).

Figure 8—figure supplement 2
Comparative FRAP of BODIPY analogues in HeLa cell cytoplasm.

(A) Structures of BODIPY analogues. (B–C) Comparative XY FRAP micrographs with time lapses in HeLa cell cytoplasm. (D–G) Micrographs of treated HeLa cells after incubations (with BDP/BDP-Ac) and, Comparative (H) averaged XY-FRAP recovery profiles with fits (N=30; R=0.99 for each of the fits), (I) averaged half-life and percentage of recoveries for BODIPY analogues in HeLa cytoplasm (H) (N=30).

Figure 9 with 1 supplement
Dconfocal of small molecules in PBS buffer and inside HeLa cells.

Comparison of Dconfocal values of small molecule drugs in (A) PBS buffer and in (B) HeLa cells. Error bars represent SE calculated from fitting the FRAP progression curves, which are averaged over at least 30 independent measurements. Dependency of Dconfocal values on (C) pKa and (D) logP values in HeLa cells. # in (A) denotes DMEM media (instead of PBS), * in (B) denotes measurements done in the nucleus.

Figure 9—figure supplement 1
Dependency of Dconfocal rates on MW, number of aromatic rings.

(A) Dependency of Dconfocal on molecular weight (MW). (B) Dependency of Dconfocal on the number of aromatic rings in the drugs.

Appendix 1—scheme 1
Chemical structures of Primaquine and N-acetylated Primaquine.
Appendix 1—scheme 2
Chemical synthesis scheme of BODIPY analogues.
Appendix 1—figure 1
Time dependent colocalization of GSK3 inhibitor and Nile Red in HeLa cell cytoplasm.

Panels (A–C) show Nile red accumulating in lipid droplets, while GSK3 inhibitor does not colocalize.

Appendix 1—figure 2
Time dependent colocalization of GSK3 inhibitor and ER marker in HeLa cell cytoplasm.

Panels (A–D) show BFP/mCherry fused antibody marker accumulating in ER, but GSK3 inhibitor does not colocalize. HeLa cells were first transiently transfected with BFP/mCherry fused antibody markers, then treated with GSK3 inhibitor.

Appendix 1—figure 3
Time dependent colocalization study of Mitoxantrone and SYTO blue marker in HeLa cells.

Panels (A–B) show SYTO blue marker accumulating in nucleic acids, but Mitoxantrone is not colocalized. HeLa cells were first treated with SYTO blue, then treated with Mitoxantrone.

Appendix 1—figure 4
Super resolution images of GSK3 inhibitor treated HeLa cells.

Micrographs were taken using a Zeiss LSM900 microscope.

Appendix 1—figure 5
Colocalization within Lysosomes.

(A–B) Comparison of percentage of colocalization among GSK3 inhibitor vs Mitoxantrone vs Quinacrine in acidic lysosomes are shown. Each dot represents 4–5 stained cells in a single frame acquisition. Colocalization study for GSK3 inhibitor and Quinacrine was done with commercial Lyso tracker-633, and for Mitoxantrone using Quinacrine itself used as a Lyso tracker dye.

Appendix 1—figure 6
Evaluating channel leakage.

Control experiments showing that channel leakage did not occur during colocalization studies. (A) Lyso-Tracker only, (B) Quinacrine only and (C) Mitoxantrone only.

Appendix 1—figure 7
Characterization of Primaquine (Prim) and N-acetylated Primaquine (PrimAc).

(A) Absorption (c ≈ 1 ∙ 10–3) and (B) emission (c ≈ 3 ∙ 10–4) spectra. Both experiments were measured in the mixture of 10% DMSO in PBS or in Glycine/NaOH buffer. (C) Measurement of hydrolytic stability by HPLC. (D) HPLC-MS measurement of PrimAc are shown. According to HPLC-MS data the peak at 7.125 corresponds to 302 m/z (ESI +) so as the main peak at 7.296. The molar mass of PrimAc corresponds to the measured value.

Appendix 1—figure 8
N-acetylated Primaquine.

(A) 1H NMR spectrum is shown. The presence of impurities is caused by the degradation in the used solvent. (B) 13C NMR spectrum is shown.

Appendix 1—figure 9
N-acetylated Primaquine.

HRMS spectrum (+EI) is shown.

Appendix 1—figure 10
Structure, photo-physical properties & Absorption/Emission spectra of BODIPY analogues.

(A) Structure and photo-physical properties of BODIPY analogues. (B) Absorption (full line) and emission spectra (dashed line) of BDP-NH2 in DMSO (black) and PBS (red) are shown. (C) Absorption (full line) and emission spectra (dashed line) of BDP-NHAc in DMSO (black) and PBS (red) are shown.

Appendix 1—figure 11
Characterization of BODIPY-NH2.

(A) 1H NMR spectrum and (B) 13C NMR spectrum are shown.

Appendix 1—figure 12
BODIPY-NH2.

HRMS spectrum (+EI) is shown.

Appendix 1—figure 13
Characterization of BODIPY-NHAc.

(A) 1H NMR spectrum and (B) 13C NMR spectrum are shown.

Appendix 1—figure 14
BODIPY-NHAc.

HRMS spectrum (+EI) is shown.

Tables

Table 1
LogP, pKa and lysosomal sequestration of drugs used in this study.

(Predicted by ChemAxon software, and obtained from drug bank data base Zhitomirsky and Assaraf, 2015).

Therapeutic indicationDrug/small moleculesCalculated LogPCalculated pKa(strongest basic)Lysosomal sequestration
Anti-cancermitoxantrone1.199.08+
Antimalarialquinacrine5.1510.33+
primaquine-biphosphate1.6410.4+
amodiaquine4.838.66+
Other small molecules5-Amino fluorescein3.029.34+
fluorescein-Sodium Salt2.652.85/3.7-
GSK3 inhibitor (SB216763)4.57+
BODIPY-NH2_9.00_
BODIPY-NHAc___
CCF24.55.1_
CF514_~3.0 (from UV spectrum)_
Values of pKa were estimated using ChemAxon. Lysosomal sequestration was taken from the in-cell micrographs in this work.
Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Gene (E.coli)nfuA, nadK, fabG, sodA, tpiA, gpmA, can, upp, speB, ispD, acuL, thiD, pyrF, nadE, crpExtracted by qPCR from BL21 E. coli K12Marciano et al., 2022
Cell line (Human)HeLa cell lineEpithelial cellRRID: CVCL_0030
Peptide, recombinant proteinnfuA, nadK, fabG, sodA, tpiA, gpmA, can, upp, speB, ispD, acuL, thiD, pyrF, nadE, crpGenes were cloned into pET-His- bdSumo and purified as describedMarciano et al., 2022
Chemical compound, drugDulbecco’s phosphate-buffered saline (PBS 1 X)Biological IndustriesCat# 02-023-1A
Chemical compound, drugHEPES bufferFisher BioReagentsLOT# 170358
Chemical compound, drugDMEM(1 X) GibcoLife Technologies LimitedREF# 41965–039
Chemical compound, drugFetal Bovine Serum (Gibco)Life Technologies LimitedREF# 12657–029
Chemical compound, drugTrypsin-EDTA Solution ABiological IndustriesREF# 03-050-1B
Chemical compound, drugPenicillin/ StreptomycinBiological IndustriesREF# 03-031-1B
Chemical compound, drugSodium Pyruvate SolutionBiological IndustriesREF# 03-042-1B
Chemical compound, drug35 mm glass-bottomed dishesMatTek CorporationP35G-0–14 C
Chemical compound, drugPierce Dye Removal ColumnsThermo-Fisher scientificCat# 22858
Chemical compound, drugGebaflex tubes
3.5 KDa (GeBa)
TIVAN BIOTECHCat# MIDI3-100
Chemical compound, drugCapillary Glass TubingWarner InstrumentsModel No. G120TF-4
203-776-0664
Chemical compound, drugwell Black F-Bottom PlateGreiner-bio-oneREF# 655076
Chemical compound, drugDisposable cuvettesFisher ScientificZEN0040
Chemical compound, drugBSA (Albumin Bovine,
fraction V)
MP Biomedicals, LLCCAS# 9048-46-8
Chemical compound, drugHEWL (Lysozyme from chicken egg white)MerckCAS# 12650-88-3
Chemical compound, drugMyoglobin (from
equine heart)
MerckCAS# 100684-32-0
Chemical compound, drugDoxorubicinAdooQ BioscienceCat# A14403
Chemical compound, drugFluorescein disodium saltchemcruzCat# sc-206026
Chemical compound, drugGSK3 inhibitor SB216763AbcamCat# ab120202
Chemical compound, drugQuinacrine dihydrochlorideAbcamCat# ab120749
Chemical compound, drugPrimaquine biphosphateMerckCAS#
63-45-6
Chemical compound, drugCCF2-FAThermo-Fisher ScientificCat# K1039a
Chemical compound, drugCF514 Labelling dyeBiotiumCF514 Dye Alternative green fluorescent dye
Chemical compound, drugAmodiaquineAdooQ BioscienceCat# A17660
Chemical compound, drugSodium AzideMerckCAS#
26628-22-8
Chemical compound, drugMitoxantroneAbcamCat# ab141041
Chemical compound, drugBafilomycin A1AbcamCat# ab120497
Chemical compound, drug5-amino FluoresceinAbcamCat#ab145305
Appendix 1—table 1
Oligomeric state, Isoelectric point and Net charge at pH 7.4.
Protein entryIsoelectric point: pINet charge at pH 7.4:ZDconfocalOligomeric stateUniprot(kDa)Monomer
NfuA4.534–15.76724.8±1.1M/D20.93
NadK6.387–4.80829.8±1.4M/D/T32.57
FabG6.250–1.6630.33±1.4D/T/H25.56
E-fts5.228–9.52722.5±1.0M/D30.29
SodA6.511–2.48125.5±1.2D22.97
tpiA5.688–7.85127.9±1.7
GpmA5.913–4.26221.7±1.0D28.43
Can6.230–5.90818.7±1.4T25.1
Upp5.335–6.82923.0±1.1M/D22.53
SpeB5.160–14.78830.4±1.4Tri/H33.56
IspD6.257–5.44524.1±1.2D/T/H25.61
BaeR5.554–6.03221.9±1.1M/D27.66
Acul5.669–5.98628.0±1.3M/D/T/H34.73
PyrF5.874–3.926.0±1.20M/D26.35
ThiD5.793–6.84228.7±1.3M/D28.64
NadE5.603–4.99323.2±1.1D27.16
CRP7.979+1.38830.1±1.3D/T23.64
DES3im-4rp45.041–5.65929.4±1.4
DES7im-3rk64.599–10.65428.8±1.3
6E5C_denovobeta306.336–0.71430.0±1.4
BSA4.5–4.8–1529.6±1.4M/D66.5
  1. Isoelectric points and net-charges were predicted using Prot Pi. Oligomeric states were taken from native mass spectrometry measurements.

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Debabrata Dey
  2. Shir Marciano
  3. Anna Poryval
  4. Ondřej Groborz
  5. Lucie Wohlrabova
  6. Tomás Slanina
  7. Gideon Schreiber
(2024)
Reversing protonation of weakly basic drugs greatly enhances intracellular diffusion and decreases lysosomal sequestration
eLife 13:RP97255.
https://doi.org/10.7554/eLife.97255.3