- Research article
- Open Access
Synthetic organotelluride compounds induce the reversal of Pdr5p mediated fluconazole resistance in Saccharomyces cerevisiae
© Reis de Sá et al.; licensee BioMed Central Ltd 2014
Received: 28 April 2014
Accepted: 18 July 2014
Published: 26 July 2014
Resistance to fluconazole, a commonly used azole antifungal, is a challenge for the treatment of fungal infections. Resistance can be mediated by overexpression of ABC transporters, which promote drug efflux that requires ATP hydrolysis. The Pdr5p ABC transporter of Saccharomyces cerevisiae is a well-known model used to study this mechanism of antifungal resistance. The present study investigated the effects of 13 synthetic compounds on Pdr5p.
Among the tested compounds, four contained a tellurium-butane group and shared structural similarities that were absent in the other tested compounds: a lateral hydrocarbon chain and an amide group. These four compounds were capable of inhibiting Pdr5p ATPase activity by more than 90%, they demonstrated IC50 values less than 2 μM and had an uncompetitive pattern of Pdr5p ATPase activity inhibition. These organotellurides did not demonstrate cytotoxicity against human erythrocytes or S. cerevisiae mutant strains (a strain that overexpress Pdr5p and a null mutant strain) even in concentrations above 100 μM. When tested at 100 μM, they could reverse the fluconazole resistance expressed by both the S. cerevisiae mutant strain that overexpress Pdr5p and a clinical isolate of Candida albicans.
We have identified four organotellurides that are promising candidates for the reversal of drug resistance mediated by drug efflux pumps. These molecules will act as scaffolds for the development of more efficient and effective efflux pump inhibitors that can be used in combination therapy with available antifungals.
The last decades have seen an increase in the immunocompromised population for several reasons including as a result of treatment of malignant diseases, HIV infection, as well as advances in organ transplantation procedures. In this scenario opportunistic infections, especially those caused by fungi, have become a serious public health problem –. Candidiasis is the most common fungal infection among patients with a condition that leads to immunosuppression ,.
Azoles, especially fluconazole, have been commonly used to treat fungal infections . However, overexpression of membrane efflux pumps by fungal cells is an important mechanism that causes azole resistance . Some of these efflux pumps belong to the Pleiotropic Drug Resistance (PDR) sub-family of ATP-Binding Cassette (ABC) transporters, and they lead to active transport of drugs using energy derived from ATP hydrolysis .
Saccharomyces cerevisiae can express several ABC transporters, and of these, Pdr5p has been the best studied . This efflux pump causes the extrusion of several drugs that are used to treat fungal infections. Also, it exhibits a profile of substrates and inhibitors that is similar to those of other ABC transporters that are expressed by pathogenic fungi . These features make Pdrp5 a good experimental model for the study of antifungal resistance mediated by ABC transporters.
One strategy for overcoming drug resistance mediated by efflux pumps is the use of compounds that can function as chemosensitizers. These compounds potentiate the efficacy of existing azoles, such as fluconazole, by inhibiting these ABC transporters . Thus, the development of novel azole chemosensitizers that increase the potency of these drugs against both sensitive and resistant fungi may allow the use of previously ineffective antifungal to treat fungal infections . Some studies have already reported compounds that are capable of reversing the resistance phenotype, such as D-Octapeptides , enniatin , isonitrile  and gallic acid derivatives .
Recently, interest in organic compounds containing tellurium (Te) or selenium (Se) has increased and several studies have been published demonstrating biological properties for both elements. Despite the relative toxicity conferred by organic compounds containing tellurium , some studies have shown that these molecules may have immunomodulatory and anti-inflammatory properties , antioxidant abilities , and anti-proliferative actions against certain tissues . Selenium is a nutritionally essential trace element for mammals. Studies have shown that some organic compounds derived from this chalcogenide exhibit antinociceptive, hepatoprotective, neuroprotective, anti-inflammatory and anti-carcinogenic properties . Furthermore, some organochalcogenides containing Te or Se are capable of inhibiting the ATPase activity of the Na+/K+ ATPase that is present in rat brains  and can inhibit the ATPase activity of P-Glycoprotein and vinblastine efflux mediated by this neoplasic cell multidrug transporter . Finally, Te and Se containing compounds can inhibit the plasma membrane H+-ATPase from S. cerevisiae.
Although several biological properties have already been described in the literature for chalcogenides and their derivatives, molecules containing selenium or tellurium with the capacity to reverse efflux pump-mediated azole resistance have not yet been reported. We were interested in studying the effects of organic compounds containing tellurium or selenium on Pdr5p, which is a well-known experimental model for the study of fungal resistance mediated by efflux pumps. In this study, we evaluated 13 synthetic compounds; some of which contained tellurium (Te) or selenium (Se), and others that were devoid of both chalcogenides.
Reagents were purchased from Sigma-USA (ATP-Sodium) or Tecoland-USA (FK 506-tacrolimus) unless otherwise stated. All reagents purchased were of highest available standard.
Synthetic compounds used in this study
Strains and culture conditions
In this study, two mutant strains of Saccharomyces cerevisiae were used. The first strain AD124567 (Pdr5p+) overexpresses Pdr5p, while the genes encoding the Pdr3p regulator and the other five ABC transporters (Yor1p, Snq2p, Pdr10p, Pdr11p and Ycf1p) have been deleted. The second one AD1234567 (Pdr5p-) contains deletions of the same six genes, as well as the gene that encodes the Pdr5p transporter . The yeast strains were grown in YPD medium (2% glucose, 1% yeast extract, 2% peptone) at 30°C with agitation and were harvested in the exponential phase of growth. One fluconazole resistant strain of Candida albicans, isolated from urine sample, was also used (approved by Instituto de Estudos em Saúde Coletiva – IESC/UFRJ – Protocol N° 030/2001). In this case, the yeast were cultivated in Sabouraud medium (4% glucose and 1% peptone), at 37°C under agitation (150 rpm).
Preparation of plasma membranes
Yeast plasma membrane isolates from the S. cerevisiae mutant strain Pdrp5+ and from the null mutant Pdr5p- were obtained as previously described by Rangel et al. . The plasma membrane preparations were stored in liquid nitrogen and thawed immediately prior to use in the Pdr5p ATPase activity assays.
ATPase activity assay
The effect of the compounds on the ATPase activity of Pdr5p was quantified by incubating Pdr5p-containing membranes (0.013 mg/mL final concentration) in a 96-well plate at 37°C for 60 min in a reaction medium containing 100 mM Tris–HCl (pH 7.5), 4 mM MgCl2, 75 mM KNO3, 7.5 mM NaN3, 0.3 mM ammonium molybdate and 3 mM ATP in the presence of the synthetic compounds. After incubation, the reaction was stopped by the addition of 1% SDS, as described previously by Dulley . The amount of released inorganic phosphate (Pi) was measured as previously described by Fiske & Subbarrow . Preparations containing plasma membranes obtained from the null mutant strain AD1234567 (Pdr5p- membranes) were used as controls. The difference between the ATPase activity of the Pdr5p + and Pdr5p- membranes represents the ATPase activity that is mediated by Pdr5p.
Effect of compounds on the growth of S. cerevisiae strains
This assay was conducted according to Niimi et al. . The effect of the compounds on the growth of both mutant strains of S.cerevisiae used in this work was determined by microdilution assays using 96-well microplates. The cells were inoculated into YPD medium at a concentration of 1 × 104 cells per well and incubated at 30°C for 48 h with agitation (150 rpm) in the presence of different concentrations of the compounds. Controls were performed using DMSO at a final concentration of 1% to verify the toxicity of the solvent used to solubilize the compounds. Cell growth was determined using a microplate reader at 600 nm (Fluostar Optima, BMG Labtech, Offenburg, Germany).
Lytic effect of compounds on human erythrocytes
This assay was conduct as described by Niimi et al. . Human erythrocytes were previously washed three times and resuspended in phosphate-buffered saline (PBS-pH 7.2). Red blood cells (final density 0.5%) were then incubate in the presence of different concentrations of the synthetic compounds for 60 min at 37°C. After incubation, the cells were pelleted by centrifugation at 3,000 g for 5 min and aliquots of 100 μL of the supernatant were transferred to the wells of a microplate. The absorbance of the hemoglobin released from the erythrocytes was measured at 540 nm. A control of 100% hemolysis was performed incubating the cells in the presence of PBS containing 1% Triton X-100.
Evaluation of fluconazole resistance reversion by the synthetic compounds
The “spot test” was used as a measure of growth as previously described by Rangel et al. . For S. cerevisiae strain Pdr5p+, 5 μL samples of fivefold serially diluted yeast cultures (initially suspended to an OD of 0.1) were spotted on YPD agar in 6 well sterile polystyrene plates. They were incubated in the presence of synthetic compounds (100 μM) only or associated with fluconazole at 120 μg/mL. Controls were performed using YPD alone and YPD supplemented with: 120 μg/mL fluconazole, 120 μg/mL fluconazole + 0.5% DMSO, 120 μg/mL fluconazole + 10 μM FK506. Plates were incubated at 30°C for 48 h.
In the case of C. albicans, the same methodology was used, but with some adaptations: 5 μL of a five-fold serial dilution from a yeast suspension containing 6 × 105 cells/mL was spotted on Sabouraud agar supplemented with the compounds at 100 μM alone or combined with fluconazole at 64 μg/mL. The incubation of the six well plates was carried at 37°C for 48 h.
Checkerboard assay with compounds and fluconazole using Candida strain from clinical isolate
Candida albicans cells, in exponential growth phase (2.5 × 103 cells/mL) were incubated in presence of different combinations of fluconazole and compound at 37°C for 48 hours in RPMI 1640 (Sigma) using 96-well plates under stirring. Cell growth was determined using a plate reader (Fluostar Optima, BMG Labtech, Germany) at a wavelength of 600 nm. The MIC value was referred to concentration capable of causing 80% growth inhibition (MIC 80). Possible synergism between fluconazole and tested compounds was determined based on the fractional inhibition concentration index (FICI). Synergic, indifferent and antagonistic interactions were defined by a FICI of <0.5, 0.5-4.0 or 4.0 respectively .
All experiments were performed in triplicate. Data were presented as mean ± standard error. A probability level of 5% (p < 0.05) in Student’s t -test was considered significant.
Results and discussion
Pdr5p is an ABC transporter and as such the inhibition of its ATPase activity could significantly affect the efflux of fluconazole and contribute to the reversal of resistance against this antifungal. Thus, a screening assay was performed to identify synthetic compounds that could promote inhibition of ATP hydrolysis catalyzed by Pdr5p (at 100 μM final concentration).
The IC 50 values of the compounds against the ATPase activity of Pdr5p
1.14 ± 0.21
1.45 ± 0.49
1.74 ± 0.91
1.48 ± 0.32
Until now, there have been no reports in the literature of organic synthetic compounds containing tellurium that inhibit Pdr5p ATPase activity. However, many other molecules, of synthetic or natural origin, also exhibit this ability. Silva et al.  demonstrated that oroidin, a derivative of a compound from a sponge, is able to inhibit the catalytic activity of this multidrug transporter with an IC50 of 20 μM. Rangel et al. , while studying gallic acid derivatives, observed that decyl gallate has an IC50 value of 13.5 μM. Both compounds competitively inhibit the enzyme activity of Pdr5p. Competitive inhibition is a more common characteristic than the uncompetitive inhibition shown by the four organotellurides.
As mentioned by Cannon et al. , inhibition of plasma membrane H+-ATPase activity could contribute to the reversal of ABC transporter-mediated azole resistance, by depleting the intracellular ATP concentration. To investigate this, the effects of the four organotellurides (1, 2, 3 and 5) on the plasma membrane H+-ATPase of S. cerevisiae were evaluated. The organotellurides leaded a powerful inhibition of the H+-ATPase activity (more than 90%) and exhibited IC50 values of approximately 2.7 μM (data not shown). Chan and colleagues  previously demonstrated that Ebselen, a well-known organoselenium compound, was also able to inhibit the activity of S. cerevisiae plasma membrane H+-ATPase in a dose dependent manner. Ebselen was also shown to be toxic for S. cerevisiae at a concentration of 10 μM, unlike the organotellurides investigated in this study.
Effect of the compounds on the growth of Pdr5p+ and Pdr5p- mutant S. cerevisiae strains
Evaluation of cytotoxicity against human erythrocytes
Fluconazole resistance reversion by the synthetic organotellurides
Checkerboard assay* using Candida albicans strain
MIC combined (μg/mL)
MIC combined (μg/mL)
Compounds 1, 2, 3 and 5 are synthetic compounds that have some similarities. Firstly, all they contain a butyl tellurium residue, secondly, they have a lateral hydrocarbon chain and finally, they all possess an amide group. All they were able to reverse the fluconazole resistance mediated by Pdr5p from S. cerevisiae. A likely mechanism for this reversal is the direct inhibition of the ATPase activity of Pdr5p and the indirect inhibition of the plasma membrane H+-ATPase.
Furthermore, the compounds could also overcome fluconazole resistance expressed by a clinical isolate of Candida albicans. The reversal of fluconazole resistance was obtained using 100 μM of the compounds. This concentration did not demonstrate toxicity against human erythrocytes or fungal cells. In conclusion, these compounds could be promising candidates for the reversal of resistance mediated by drug efflux pumps, act synergistically with fluconazole and could serve as prototypes for the synthesis of other molecules that could be capable of inhibiting efflux pumps with greater efficiency.
Availability of supporting data
The data sets supporting the results of this article are included within the article.
The authors thank FAPERJ (E-26/111.338/2013), FAPESP (2005/59572-7, 2008/55401-1, 2010/17228-6, 2011/03244-2, 2011/11613-8 and 2012/17093-9), CNPq (470360/2012-7) and CAPES for financial support and scholarships. The authors are grateful for the financial and structural support offered by the University of São Paulo through the NAP-CatSinQ (Research Core in Catalysis and Chemical Synthesis). The authors thank also to our lab assistant, Mrs. Geralda Rodrigues Almeida for her great support and Dr. Louise Kemp for your critical reading of this manuscript.
- Brown GD, Meintjes G, Kolls JK, Gray C, Horsnell W: AIDS-related mycoses: the way forward. Trends Microbiol. 2014, 22 (3): 107-109. 10.1016/j.tim.2013.12.008.PubMed CentralView ArticlePubMedGoogle Scholar
- Calton EA, Le Doaré K, Appleby G, Chisholm JC, Sharland M, Ladhani SN, CABIN Participants: Invasive bacterial and fungal infections in paediatric patients with cancer: incidence, risk factors, aetiology and outcomes in a UK regional cohort 2009–2011.Pediatr Blood Cancer 2014, doi:10.1002/pbc.,Google Scholar
- Kauffman CA, Freifeld AG, Andes DR, Baddley JW, Herwaldt L, Walker RC, Alexander BD, Anaissie EJ, Benedict K, Ito JI, Knapp KM, Lyon GM, Marr KA, Morrison VA, Park BJ, Patterson TF, Schuster MG, Chiller TM, Pappas PG: Endemic fungal infections in solid organ and hematopoietic cell transplant recipients enrolled in the Transplant-Associated Infection Surveillance Network (TRANSNET). Transpl Infect Dis. 2014, 0: 1-12.Google Scholar
- Yapar N: Epidemiology and risk factors for invasive candidiasis. Ther Clin Risk Manag. 2014, 10: 95-105. 10.2147/TCRM.S40160.PubMed CentralView ArticlePubMedGoogle Scholar
- Wille MP, Guimarães T, Furtado GHC, Colombo AL: Historical trends in the epidemiology of candidaemia: analysis of an 11-year period in a tertiary care hospital in Brazil. Mem Inst Oswaldo Cruz. 2013, 108 (3): 288-292. 10.1590/S0074-02762013000300005.PubMed CentralView ArticleGoogle Scholar
- Odds FC, Brown AJ, Gow NA: Antifungal agents: mechanisms of action. Trends Microbiol. 2003, 11: 272-279. 10.1016/S0966-842X(03)00117-3.View ArticlePubMedGoogle Scholar
- Martinez L, Falson P: Multidrug resistance ATP-binding cassette membrane transporters as targets for improving oropharyngeal candidiasis treatment. Adv Cell Mol Otolaryngol. 2014, 2: 1-8. 10.3402/acmo.v2.23955.View ArticleGoogle Scholar
- Prasad R, Goffeau A: Yeast ATP-binding cassette transporters conferring multidrug resistance. Annu Rev Microbiol. 2012, 66: 39-63. 10.1146/annurev-micro-092611-150111.View ArticlePubMedGoogle Scholar
- Morschhäuser J: Regulation of multidrug resistance in pathogenic fungi. Fungal Genet Biol. 2010, 47: 94-106. 10.1016/j.fgb.2009.08.002.View ArticlePubMedGoogle Scholar
- Wolfger H, Mamnun Y, Kuchler MK: Fungal ABC proteins: pleiotropic drug resistance, stress response and cellular detoxification. Res Microbiol. 2001, 152: 375-389. 10.1016/S0923-2508(01)01209-8.View ArticlePubMedGoogle Scholar
- Cannon RD, Lamping E, Holmes AR, Niimi K, Baret PV, Keniya V, Tanabe K, Niimi M, Goffeau A, Monk BC: Efflux-mediated antifungal drug resistance. Clin Microbiol Rev. 2009, 22: 291-321. 10.1128/CMR.00051-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Niimi K, Harding DRK, Parshot R, King L, Decottignies A, Niimi M, Lin S, Cannon RD, Goffeau A, Monk BC: Chemosensitization of fluconazole resistance in Saccharomyces cerevisiae and pathogenic fungi by a D-octapeptide derivative. Antimicrob Agents Chemother. 2004, 48: 1256-1271. 10.1128/AAC.48.4.1256-1271.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Hiraga K, Yamamoto S, Fukuda HN, Oda K: Enniatin has a new function as an inhibitor of Pdr5p, one of the ABC transporters in Saccharomyces cerevisiae. Biochem Biophys Res Commun. 2005, 328: 1119-11125. 10.1016/j.bbrc.2005.01.075.View ArticlePubMedGoogle Scholar
- Yamamoto S, Hiraga K, Abiko A, Hamanaka N, Oda K: A new function of isonitrile as an inhibitor of the Pdr5p multidrug ABC transporter in Saccharomyces cerevisiae. Biochem Biophys Res Commun. 2005, 330: 622-628. 10.1016/j.bbrc.2005.03.009.View ArticlePubMedGoogle Scholar
- Rangel LP, Fritzen M, Yunes RA, Leal PC, Creczynski-Pasa TB, Ferreira-Pereira A, Fritzen M, Yunes RA, Leal PC, Creczynski-Pasa TB, Ferreira-Pereira A: Inibitory effects of gallic acid ester derivates on Saccharomyces cerevisiae multidrug resistance protein Pdr5. FEMS Yeast Res. 2010, 10: 244-251. 10.1111/j.1567-1364.2009.00603.x.View ArticleGoogle Scholar
- Schiar VP, Dos-Santos DB, Paixão MW, Nogueira CW, Rocha JBT, Zeni G: Human erythrocite hemolysis induced by selenium and tellurium compounds increased by GSH or glucose: a possible envolvement of ractive oxygen species. Chem Biol Interact. 2009, 177: 28-33. 10.1016/j.cbi.2008.10.007.View ArticlePubMedGoogle Scholar
- Sredni-Kenigsbunch D, Shohat M, Shohat B, Ben-Amitai D, Chan CC, David M: The novel tellurium immunomodulator AS 101 inhibits interleukin-10 production and p 38 MAPK expression in atopic dermatidis. J Dermato Sci. 2008, 50: 232-235. 10.1016/j.jdermsci.2007.12.007.View ArticleGoogle Scholar
- Ren X, Xue Y, Liu JK, Zheng J, Luo G, Guo C, Mu Y, Shen J: A novel cyclodextrin-derived tellurium compound with glutathione peroxidase activity. Chembiochem. 2002, 3: 356-363. 10.1002/1439-7633(20020402)3:4<356::AID-CBIC356>3.0.CO;2-O.View ArticlePubMedGoogle Scholar
- Kalechman Y, Gaffer U, Weinstein T, Chagnac A, Freidkin I, Tobar A, Albec M, Sredni B: Inhibition of interleukin-10 by the immunomodulator AS 101 reduces mesangial cell proliferation in experimental mesangioproliferative gromerulonephrits: association with dephosphorilation of STAT 3. J Biol Chem. 2004, 279: 24784-24732. 10.1074/jbc.M312006200.View ArticleGoogle Scholar
- Nogueira CW, Zeni G, Rocha JBT: Organoselenium and organotellurium compounds: toxicology and pharmacology. Chem Rev. 2004, 104: 6255-6286. 10.1021/cr0406559.View ArticlePubMedGoogle Scholar
- Borges VC, Rocha JBT, Nogueira CW: Effect of diphenyl diselenide, diphenyl ditelluride and ebselen on cerebral Na+, K+-ATPase activity in rats. Toxicology. 2005, 25: 191-197. 10.1016/j.tox.2005.07.002.View ArticleGoogle Scholar
- Sawada GA, Raub TJ, Higgins JW, Brenan NK, Moore TM, Tombline G, Detty MR: Chalcogenopyrylium dyes as inhibitors/modulators of P-glycoprotein in multidrug-resistant cells. Bioorg Med Chem. 2008, 16: 9745-9756. 10.1016/j.bmc.2008.09.065.View ArticlePubMedGoogle Scholar
- Chan G, Hardej D, Santoro M, Lau-Cam C, Billack B: Evaluation of the antimicrobial activity of ebselen: role of the yeast plasma membrane H+-ATPase. J Biochem Molecular Toxicology. 2007, 21: 252-264. 10.1002/jbt.20189.View ArticleGoogle Scholar
- Keppler AF, Gariani RA, Lopes DG, Comasseto JV: Lithium butylchalcogenolate induced Michael-aldol tandem sequence: Easy and rapid access to highly functionalized organochalcogenides and unsaturated compounds. Tetrahedron Lett. 2009, 50: 2181-2184. 10.1016/j.tetlet.2009.02.158.View ArticleGoogle Scholar
- Vargas F, Toledo FT, Comasseto JV:N-Functionalized organolithium compounds via tellurium/lithium exchange reaction. J Braz Chem Soc. 2010, 21: 2072-2078. 10.1590/S0103-50532010001100007.View ArticleGoogle Scholar
- Sousa BA, Keppler AF, Gariani RA, Comasseto JV, Dos Santos AA: Metallic chalcogenolates mediated modular Michael-aldol cascade reaction: an easy route to multi-functionalized chalcogenides and Morita-baylis-Hillman adducts. Tetrahedron. 2012, 68: 10406-10413. 10.1016/j.tet.2012.09.027.View ArticleGoogle Scholar
- Sousa BA, Dos Santos AA: A facile, versatile, and mild Morita-baylis-Hillman-type reaction for the modular one-pot synthesis of highly functionalized MBH adducts. Eur J Org Chem. 2012, 2012 (18): 3431-3436. 10.1002/ejoc.201200371.View ArticleGoogle Scholar
- Decottignies A, Grant AM, Nichols JW, de Wet H, McIntosh DB, Goffeau A: ATPase and multidrug transport activities of the overexpressed yeast ABC protein Yor1p. J Biol Chem. 1998, 273: 12612-12622. 10.1074/jbc.273.20.12612.View ArticlePubMedGoogle Scholar
- Dulley J: Determination of inorganic phosphate in the presence of detergents or protein. Anal Biochem. 1965, 67: 91-96. 10.1016/0003-2697(75)90275-4.View ArticleGoogle Scholar
- Fiske CH, Subbarow YJ: The colorimetric determination of phosphorus. J Biol Chem. 1925, 66: 375-400.Google Scholar
- Mukherjee PK, Sheehan DJ, Hitchcock CA, Ghannoum MA: Combination treatment of invasive fungal infections. Clin Microbiol Rev. 2005, 18 (1): 163-194. 10.1128/CMR.18.1.163-194.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Silva FR, Tessis AC, Ferreira PF, Rangel LP, Garcia-Gomes AS, Pereira FR, Berlinck RGS, Muricy G, Ferreira-Pereira A: Oroidin inhibits the activity of the multidrug resistance target Pdr5p from yeast plasma membranes. J Nat Prod. 2011, 74: 279-282. 10.1021/np1006247.View ArticlePubMedGoogle Scholar
- Egner R, Bauer BE, Kuchler K: The transmembrane domain 10 of the yeast Pdr5p ABC antifungal efflux pump determines both substrate specificity and inhibitor susceptibility. Mol Microbiol. 2000, 35 (5): 1255-1263. 10.1046/j.1365-2958.2000.01798.x.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.