Chemogenomic profiling to understand the antifungal action of a bioactive aurone compound
Authors:
Fatmah M. Alqahtani aff001; Brock A. Arivett aff001; Zachary E. Taylor aff002; Scott T. Handy aff002; Anthony L. Farone aff001; Mary B. Farone aff001
Authors place of work:
Department of Biology, Middle Tennessee State University, Murfreesboro, Tennessee, United States of America
aff001; Department of Chemistry, Middle Tennessee State University, Murfreesboro, Tennessee, United States of America
aff002
Published in the journal:
PLoS ONE 14(12)
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pone.0226068
Summary
Every year, more than 250,000 invasive candidiasis infections are reported with 50,000 deaths worldwide. The limited number of antifungal agents necessitates the need for alternative antifungals with potential novel targets. The 2-benzylidenebenzofuran-3-(2H)-ones have become an attractive scaffold for antifungal drug design. This study aimed to determine the antifungal activity of a synthetic aurone compound and characterize its mode of action. Using the broth microdilution method, aurone SH1009 exhibited inhibition against C. albicans, including resistant isolates, as well as C. glabrata, and C. tropicalis with IC50 values of 4–29 μM. Cytotoxicity assays using human THP-1, HepG2, and A549 human cell lines showed selective toxicity toward fungal cells. The mode of action for SH1009 was characterized using chemical-genetic interaction via haploinsufficiency (HIP) and homozygous (HOP) profiling of a uniquely barcoded Saccharomyces cerevisiae mutant collection. Approximately 5300 mutants were competitively treated with SH1009 followed by DNA extraction, amplification of unique barcodes, and quantification of each mutant using multiplexed next-generation sequencing. Barcode post-sequencing analysis revealed 238 sensitive and resistant mutants that significantly (FDR P values ≤ 0.05) responded to aurone SH1009. The enrichment analysis of KEGG pathways and gene ontology demonstrated the cell cycle pathway as the most significantly enriched pathway along with DNA replication, cell division, actin cytoskeleton organization, and endocytosis. Phenotypic studies of these significantly enriched responses were validated in C. albicans. Flow cytometric analysis of SH1009-treated C. albicans revealed a significant accumulation of cells in G1 phase, indicating cell cycle arrest. Fluorescence microscopy detected abnormally interrupted actin dynamics, resulting in enlarged, unbudded cells. RT-qPCR confirmed the effects of SH1009 in differentially expressed cell cycle, actin polymerization, and signal transduction genes. These findings indicate the target of SH1009 as a cell cycle-dependent organization of the actin cytoskeleton, suggesting a novel mode of action of the aurone compound as an antifungal inhibitor.
Keywords:
Saccharomyces cerevisiae – Cell cycle and cell division – Antimicrobial resistance – Candida albicans – Antifungals – Actins – Cytoskeleton – Endocytosis
Introduction
Life-threatening fungal infections have been increasing due to the difficulties with diagnosis and treatment that accelerate mortality rates associated with fungal infections, which now exceed deaths caused by malaria [1]. Candida albicans is the most frequently isolated opportunistic fungal pathogen and is implicated in superficial mucosal infections, or candidiasis of the oral cavity and genitalia of humans, particularly in immunocompromised patients [2]. In healthy individuals, Candida spp. are a commensals of the mucosal surfaces of genitalia, oral cavity, and gastrointestinal tract. However, with the introduction of antibacterial antibiotics as medical therapy in the 1940s, a gradual increase in the number of invasive candidiasis cases has been reported due to antibiotic-associated loss of the bacterial biota and subsequent colonization of Candida spp. on epithelial surfaces, a requirement for pathogenesis [3]. Several risk factors contribute to the pathogenesis of invasive candidiasis, including organ transplantation, prolonged hospitalization in an intensive care unit, catheterization, and intensive utilization of antibiotics and immunosuppressive agents. These factors could lead Candida spp. to colonize mucosal surfaces, resulting in superficial infections. The fungus can also advance to candidemia, or invasion of the bloodstream, and from there disseminate to different organs. Certain virulence factors are attributed to the pathogenicity of Candida spp., including adherence to epithelial surfaces, dimorphic growth, biofilm formation, and production of tissue-damaging enzymes [4, 5].
For treating candidiasis, there are five groups of antifungal agents as defined by their mode of action and for which mechanisms of resistance have been described. Group I: polyenes (amphotericin B) bind to ergosterol in the cell membrane and form pores in it, while Group II: echinocandins (caspofungin) inhibit β(1,3)-glucan synthase in the cell wall. Group III: azoles (fluconazole) inhibit lanosterol 14 α-demethylase in the ergosterol biosynthesis pathway. Group IV: synthetic pyrimidines (5-fluorocytosine) inhibit DNA synthesis and disturb protein synthesis, and Group V: allylamines (terbinafine) inhibit squalene epoxidase in the ergosterol biosynthesis pathway [6]. Resistance mechanisms have been described to these antifungals as cellular determinants that lead to drug extrusion by active efflux, altered drug targets, or drug target overexpression. However, novel drug resistance mechanisms have been recently reported as robust responses that enhance antifungal tolerance by pathways such as regulation of the oxidative or thermal stress responses [7].
Even with treatment by commercially-available antifungal agents, the mortality rate from disseminated candidiasis has surged to ~40–60%, representing a 20-fold increase compared to only two decades ago [3]. Every year, more than 250,000 invasive candidiasis infections are reported with 50,000 deaths worldwide [4]. Furthermore, in the USA alone, the cost of combating candidiasis was estimated to be $2–4 billion annually in the year 2000 [8]. Candidiasis has recently been reported as the third-to-fourth most frequent healthcare-acquired infection globally [9]. Although the majority of candidiasis cases in humans are attributed to C. albicans, other Candida species have not only emerged as causative agents of candidiasis but have also developed resistance to antifungal drugs. These species most often include C. glabrata, C. tropicalis, C. parapsilosis, and C. krusei [3]. The expanding immunosuppressed population, the limited number of fundamental antifungal agents along with their resistances and toxicity issues, and the emergence of non-albicans pathogenic strains all necessitate the need to seek alternative antifungal agents with potential novel targets.
To achieve this goal of seeking alternative antifungals, the exploitation of natural products, particularly those derived from plants, appears to be a promising source for antifungal compound development [10]. Because plants have their own fungal pathogens, these interactions between plants and fungi have resulted in the origination of diverse chemical entities within the plants intended to enhance not only their protection from fungal pathogens, but their survival and competitiveness as well [11]. One of the most promising classes of natural products are the secondary metabolites, aurones, which are ubiquitously present in plants. Aurones, or 2-benzylidenebenzofuran-3-(2H)-ones, are structural isomers of flavonoids that naturally occur as yellow-color pigments in plants [12–14]. In addition to their roles in pigmentation, aurones possess a variety of protective roles in the plant, including insect antifeedant [13], antiparasitic [14], and antifungal activities [15, 16]. Moreover, activities of aurones as anticancer [16, 17], antiparasitic [14], antileishmanial [18, 19], and antifungal agents have been reported [20]. Since the bioactive properties and therapeutic prospective of natural and synthetic aurones are promising, these bioactive compounds can be considered as an attractive scaffold for antifungal drug design and development.
We have previously reported the synthesis and anti-Candida activities of non-natural aurone derivatives containing different functional groups, including aurone SH1009 (Fig 1), which exhibited significant inhibition of C. albicans when compared to the other derivatives [21]. The reported disruption of biofilm formation by aurone SH1009 emphasizes the importance of understanding the SH1009 mode of action. In the present study, the antifungal activity of aurone SH1009 was determined against different standard and clinical Candida spp., including resistant isolates, applying a modified synthetic strategy based on an acid-mediated condensation between the appropriate benzofuranone and aldehyde. The mode of antifungal action of aurone SH1009 was characterized using chemogenomic approach in Saccharomyces cerevisiae mutant collections and validated in C. albicans SC5314.
Results
Aurone SH1009 is selectively inhibitory for Candida spp.
To assess the antifungal activity of aurone SH1009 against Candida strains (listed in S1 Table), antifungal susceptibility testing was performed using the Clinical Laboratory Standards Institute (CLSI) broth microdilution protocol [22]. SH1009 exhibited promising antifungal activity against all Candida spp. tested (Table 1). C. albicans strains M2, M3, M5, and M7 are all fluconazole resistant and, as such, represent a serious risk for patients since fluconazole is the principle antifungal drug for treating candidiasis [2]. In this study, fluconazole had no inhibitory activity on these isolates even at high concentrations (MIC > 64 μg/mL), while the IC50 of aurone SH1009 is considerably lower [11, 12, 4, and 8 μM) for C. albicans M2, M3, M5, and M7 strains, respectively. Similarly, another aggressive, multi-drug resistant isolate, C. albicans ATCC 64124, exhibited resistance to a high concentrations of amphotericin B (16 μg/mL), caspofungin (4 μg/mL), and fluconazole (64 μg/mL), although 5-fluorocytosine was inhibitory (MIC < 1 μg/mL). This strain was also sensitive to SH1009 (IC50 = 21 μM). Although C. albicans is recovered from at least 50% of candidiasis cases, C. glabrata has recently emerged as responsible for approximately 25% of cases in the last two decades, and tends to infect elderly populations [23]. Exposure of C. glabrata to SH1009 resulted in inhibition with a low concentration of aurone (IC50 < 3.125 μM), suggesting it could be a successful treatment, especially since the development of C. glabrata resistance to caspofungin and lower susceptibility to fluconazole have been reported [24]. C. tropicalis, a causative agent in 10–20% of candidiasis cases in the USA and 35–40% of cases in tropical regions, was also susceptible to SH1009 [3]. Because aurone SH1009 is not affected by the resistance mechanisms of these isolates with known resistances to antifungals, the compound presents new possibilities for further exploration as a potential antifungal agent.
Although SH1009 was inhibitory, the aurone was not fungicidal even at the highest concentration of aurone tested (200 μM), as indicated by colony growth on solid media. Therefore, a more sensitive cell viability assay, compatible with yeast cells, was required to more precisely measure viable cell number after treatment. Along with flow cytometry, the live/dead Fungalight Yeast Viability Kit was used to quantitate viability after SH1009 treatment. After 48 h of treatment with a high SH1009 concentration (200 μM), C. albicans SC5314 strain (Fig 2A) revealed two divided subpopulations after gating for the stained cells in the dot plot graphs. When compared to isopropanol-treated cells (Fig 2B), a significant fraction (Fig 2D) of SH1009-treated cells clustered in the upper-left quadrant, indicating the bioactivity of SH1009 in eliminating C. albicans SC5314 growth. In contrast, the viable cells in the SH1009-treated sample (Fig 2A) were reduced by approximately half compared to the untreated-cells sample (Fig 2C).
Growth of SH1009-treated C. albicans SC5314 was monitored at 30-min intervals over a 46 h incubation using dilutions of the aurone. A representative set of growth curves of C. albicans SC5314 is given in Fig 3A. The curves indicate that there is no regrowth of C. albicans SC5314 in the presence of SH1009 above 12.5 μM. Additionally, inhibition of C. albicans SC5314 by SH1009 is highly dose dependent. The potency of any bioactive-compound is usually measured by the concentration of compound that inhibits 50% (IC50) of a fixed-pathogenic inoculum in a dose–response assay in vitro. However, dose–response curve slope is another clinically independent criterion that can be used as an indicator of the expected therapeutic effectiveness [25]. Typically, a dose–response curve with a particularly steep slope (Hill coefficient >1) indicates that a small increase in concentration of a drug above the IC50 causes extraordinarily high-level inhibition. SH1009 exhibited a steeper dose-response IC50 curve in inhibition of C. albicans SC5314 as is indicated by the higher slope factor (~2.241) (Fig 3B). The steep slope of the IC50 curve reflects the growth curve results for treated cells, in that there is significant loss of viability with aurone concentrations above the IC50 concentration of 16.28 μM. This predicts a possible therapeutic potency for SH1009 with a selective 50% inhibitory concentration required to inhibit C. albicans SC5314 cells (0.5–2.5×103 CFU/mL).
A time-kill assay was employed to investigate the killing kinetics of aurone SH1009 against C. albicans SC5314 along with fluconazole and amphotericin B treatments (Fig 3C). As expected, amphotericin B affected the growth curve at 8 μg/mL after 2 h, regardless of the initial inoculum of C. albicans SC5314 cells (0.5–0.25×105 CFU/mL), demonstrating fungicidal activity. Treatment with SH1009 at a concentration representing an approximate five-fold increase of the IC50 calculated from the two-fold dilution series (500 μM) reduced the growth of the yeast similarly to fluconazole (16 μg/mL) with no significant differences until 12h of incubation, indicating fungistatic activity. Given the initial inoculum, the SH1009 treatment yielded an approximate > 3 log10 decrease in CFU/mL after each time point of treatment compared with untreated culture, exhibiting a significant reduction (P value ≤0.0001) in the colony count. This reduction in the CFU/mL resulted in no appreciably significant increase in growth rate after SH1009 treatment over 24h, whereas the untreated culture recorded > 4.4 log10 significant increase in the CFU/mL between 0 and 24h. These results indicate a fungistatic activity which might be less desirable for immunosuppressed patients to remedy opportunistic infections. Consequently, we measured the cytotoxicity of SH1009 for three different human cell lines.
Toxicity assays to determine the cytotoxic concentration of aurone that reduced cell viability by 50% (CC50) were carried out based on the reduction of resazurin to resorufin by metabolically active cells as a sensitive method for detecting viable cells [26]. The CC50 values for THP-1, HepG2, A549 cells after 24h treatment with two-fold serial dilutions of SH1009 (3.125–200 μM) were 140, 168, and >200 μM, respectively (Table 2). Selectivity index (SI) values were calculated to correlate the antifungal activity of aurone SH1009 against C. albicans SC5314 cells (IC50 16.28 μM) with the CC50 concentration for the human cell lines. The CC50/IC50 ratio was an ~8.6->12-fold difference in the concentrations that resulted in 50% loss of cell viability, suggesting a selectivity of SH1009 for the pathogenic yeast cells over human cells. Treatment of C. albicans SC5314 cells with SH1009 at 25 μM resulted in a significant reduction (P ≤ 0.001) in cell viability when compared to the untreated control group, while treatment of human cell lines with the same concentration resulted in no appreciably significant reduction with 70% cell viability for THP-1 cells and 100% viability for both HepG2 and A549 cells. With increasing SH1009 concentrations, the cell viability of C. albicans SC5314 was reduced significantly (P ≤ 0.01), compared to the cell viability of human cells (Fig 3D). These results indicate that aurone SH1009 has a selective toxicity for C. albicans cells. To assist with assessing the therapeutic potential of aurone SH1009 as an antifungal, we next sought to determine the bioactivity of the compound by defining its mode of toxicity.
Chemogenomic profiling of aurone SH1009-treated yeast cells identifies roles for genes involved in the cell cycle, cell division, and the actin cytoskeleton
Chemogenomic profiling was used to characterize the mode of action for aurone SH1009. Haploinsufficiency profiling (HIP) and homozygous profiling (HOP) allow for paralleled assessment of the sensitivity and resistance of the pooled genome-wide set of S. cerevisiae deletion mutants. First, two pools of S. cerevisiae heterozygous (HIP) and homozygous (HOP) deletion mutant collections were treated with aurone SH1009 at the concentration (~500 uM) that inhibited the growth of the wild type S. cerevisiae-S288C (1.25× 109 cells /mL) by 20% for 48h. After purifying genomic DNA from the mutants, the synthetic UPTAG DNA barcodes (20 bp) were amplified using uniquely indexed primers to distinguish each sample. Because each mutant is uniquely identified with DNA barcode, multiplexed-next generation sequencing as a highly robust technique was employed to quantitate the abundance of each mutant.
Haploinsufficiency profiling was performed on ~1056 heterozygous mutants that are essential for growth and express only 50% of gene dosage because one functional copy of that particular gene in the diploid organism has been deleted. Whereby, the identification of the direct target of a certain bioactive compound can be identified in the presence of that compound as the mutant that has a large fitness defect compared to the other mutants that do not encode the drug target. Conversely, homozygous profiling was performed on ~4244 of homozygous mutants that are non-essential for growth and express 0% of gene dosage because both copies of the particular gene are deleted in the diploid organism. With the HOP assay, it is possible to suppress drug sensitivity due to the complete loss-of-function alleles, allowing identification of pathways that confer the drug sensitivity or identification of the direct target of drugs following the principle that deletion of drug target will render the cells insensitive to the compound [27].
The post-sequencing data analysis of the aurone-treated deletion mutants revealed 3923 mutants, which included 3,133 mutants from the homozygous deletion pool (non-essential genes) and 790 mutants from the heterozygous deletion pool (essential genes), representing ~ 75% of the mutant population with usable read counts and 0.90 correlation between samples, indicating high sample quality and agreement. The chemical-genetic interaction of the positive control methyl methanesulfonate (MMS), a well-characterized antifungal agent that damages DNA [28], demonstrated a highly significant enrichment (P value ≤ 0.001) for cellular response to DNA damage stimulus, confirming a successful HIP-HOP assay procedure and accurate post-sequencing data analysis. The chemical-genetic interaction identified 238 gene deletion mutants that were significantly responsive to aurone SH1009 (FDR P values ≤ 0.05) for both HIP-HOP profiles (Fig 4A and 4B). The sensitive and resistant genes with that P value (approximately ≥ 1.5 fold change) for HIP and HOP independently were used for gene ontology (GO) enrichment analysis using ClueGO [29] and Yeast Gene Ontology Slim Term Mapper at Saccharomyces Genome Database (SGD) [30].
For essential genes (HIP profile), GO enrichment (P ≤ 0.0001) for mutants that have deleted genes encoding cell cycle proteins were found (e.g. CDC28Δ, CDC34Δ, CDC25Δ, CDC13Δ, ORC4Δ, MCM2Δ, MCM7Δ, CBK1Δ, and RFC1Δ) (Fig 4C). Moreover, two other biological responses involved in the cell cycle showed a significant GO enrichment (P ≤ 0.01). These included genes for meiosis-yeast (ORC4Δ, MCM2Δ, MCM7Δ, APC4Δ, APC1Δ, CDC28Δ, and CDC5Δ) and DNA replication (POL3Δ, MCM2Δ, MCM7Δ, RFC3Δ and RFC1Δ). For non-essential genes (HOP profile), among both sensitive and resistant mutants, GO enrichment (P ≤ 0.01) detected sensitive mutants (ARK1Δ, NIP100Δ, GIC1Δ, SLK19Δ, KIP3Δ, and ACF4Δ) and resistant mutants (BOI2Δ, EDEΔ, END3Δ, and GIC2Δ) with deleted genes that are associated with actin cytoskeleton organization (Fig 3D). Also, mutants that are associated with cell division were significantly enriched from the HOP profile (LDB19Δ, SLK19Δ, LTE1Δ, ZIP2Δ, BUB1Δ, ELM1Δ, TOF1Δ, TOF2Δ, SWI4Δ, CSM3Δ, CLB3Δ, and CLN2Δ). Lastly, there was an enrichment of deletion mutants for the biological response of endocytosis, which also relates to the actin cytoskeleton, that included mutants that were among the top 10 sensitive in the HOP profile (FEN2Δ, LTE1Δ, LDB19Δ, and ARK1Δ) (S2 and S3 Tables).
For several phenotypes described above, the cell cycle pathway appeared to be the primary target pathway that was significantly enriched by SH1009 treatment for the essential genes profile and non-essential genes profile, separately. To obtain a broader view of the changes in growth patterns of heterozygous and homozygous mutants after SH1009 treatment, the significantly responsive mutants from both HIP-HOP profiles were combined and used for the enrichment analysis. For all 238 deleted protein-coding genes with (FDR P values ≤ 0.05 and fold change ≥ 1.5) from both profiles, the clustered analysis of KEGG pathways and gene ontology terms that significantly enriched (FDR P values ≤ 0.05) were plotted for the number of genes that are associated with each term in a histogram (Fig 5). Notably, both the profiles for the 80 essential genes (HIP-profile) and 158 non-essential genes (HOP-profile) are largely clustered for most the significantly enriched category (S4 Table). These significantly responsive genes were mapped to six pathways in the KEGG pathway database, with the cell cycle pathway as the most significantly enriched pathway (P ≤ 0.0001) along with other pathways that are completely overlapped with cell cycle pathway (meiosis and DNA replication). For gene ontology categories, within the biological process category 18 terms were enriched in differentially sensitive or resistant mutants, including cell division, cytoskeleton organization, and regulation of endocytosis. Also, the nucleotide-binding, aminoacyl-tRNA ligase activity, and DNA-binding terms were significantly enriched in the molecular function category. Nucleus, cytoskeleton, cellular bud, as well as the site of polarized growth terms were significantly enriched in the cellular component category.
In order to gain insight into the chemical-genetic interaction of genes targeted by aurone SH1009 and visualize the connections between targeted biological responses and pathways identified by both the HIP and HOP profiles, ClueGO and CluePedia along with Cytoscape software were used to extract and map non-redundant biological responses for a large set of functionally clustered genes using GO terms and KEGG pathway, simultaneously [29, 31]. Fig 6A depicts the functional annotation network of clustered essential and non-essential genes. This overlapped network revealed relatively connected biological categories that started from the highest significantly enriched pathway (cell cycle) to the less significant pathways (meiosis, DNA replication, endocytosis, and RNA biogenesis). When focusing on the cell cycle pathway only (Fig 6B), the non-essential genes (HOP profile) were strongly clustered with the essential genes (HIP profile) for significantly enriched biological processes associated with actin cytoskeleton organization and endocytosis. Eight clustered genes for actin cytoskeleton and endocytosis were found among the top 20 sensitive mutants of both the HIP and HOP profiles, implicating cell-cycle-dependent organization of the actin cytoskeleton and endocytosis as targets of aurone SH1009. Additionally, the most significantly enriched molecular function for the SH1009-responsive mutants was the nucleotide-binding protein. The annotation network revealed 51 nucleotide-binding proteins that significantly responded to aurone SH1009; 31 of these genes encode proteins that act as ATP-binding proteins, while the other 10 encode GTP-binding proteins. Taken together, the functional enrichment analysis of the deletion pool of essential and non-essential genes illustrates that aurone SH1009 treatment might target nucleotide binding proteins, leading to a series of cellular defects that belong to cell cycle pathway, actin cytoskeleton organization, and endocytosis.
Before assessing biological responses to the aurone, growth curves were used to confirm the sensitivity or resistance responses of the individual mutants to SH1009 [32]. Twelve mutants that had more than or equal to a two-fold change and were also involved in the suggested biological responses from both HIP and HOP profiles were tested individually using dilutions of SH1009 (3.125–200 μM). The differences in the log-phase of the growth between the wild type S. cerevisiae-S288C and the mutants treated with SH1009 was compared. Fig 7 and S1 Fig depict the growth curves of individual S. cerevisiae deletion mutants for confirmation of the growth patterns of the HIP and HOP profiles.
The individual, heterozygous sensitive strains TSC11Δ, YPP1Δ, CDC42Δ, CDC25Δ, CDC28Δ, and RSP5Δ from chemical-genetic profiling all showed reduced growth with SH1009 treatment when compared to the log-phase growth of S. cerevisiae-S288C wild type, which supports results from chemical-genetic interaction analysis (Fig 7A). For the homozygous mutant strains, the same cellular effects were observed with decreased growth for sensitive mutants (LTE1Δ, FEN2Δ, and ARK1Δ) and increased growth for resistant mutants (CLN2Δ, LIN1Δ, and ZIP2Δ) in the presence of SH1009 when compared to the wild type S. cerevisiae-S288C (Fig 7B). Additionally, the growth of each mutant in the IC50 concentration of SH1009 (16 μM) was compared independently to the growth of S. cerevisiae-S288C to assess significant differences between the growth of each mutant and the wild type. All mutants demonstrated significant differences in growth from the wild type (P ≤ 0.05–0.0001), confirming their expected growth patterns as sensitive or resistant mutants to the aurone SH1009 and supporting the enrichment analysis results since these mutants harbor deletions of genes that are associated with cell cycle progression (CDC25Δ, CDC28Δ, CLN2Δ, LIN1Δ, and ZIP2Δ), actin cytoskeleton organization (TSC11Δ, CDC42Δ, RSP5Δ, LTE1Δ, and ARK1Δ), and endocytosis (TSC11Δ, YPP1Δ, FEN2Δ, and RSP5Δ) (Fig 7C). Once we confirmed the growth patterns of individual S. cerevisiae mutants to aurone SH1009, C. albicans was treated with the aurone to detect phenotypic changes associated with the suggested biological responses in the pathogen.
Aurone SH1009 blocks cell cycle progression in Candida albicans
Cell cycle was the most significantly enriched pathway in the chemical-genetic interaction analysis (Figs 5 and 6). The other enriched pathways, meiosis and DNA replication, also contribute to the cell cycle pathway. Accordingly, if aurone SH1009 targets cell cycle gene-encoded proteins, the distribution of cell cycle phases in the fungal cell population during exponential phase should be altered compared to the normal distribution, indicating cell cycle arrest. For these experiments, changes in DNA content throughout different cell cycle phases was assessed by flow cytometry, which allows quantitative single cell detection. By labeling cellular DNA with propidium iodide (PI), the cells can be quantitatively discriminated in different phases of the cell cycle based on the fluorescence intensity, which is proportional to specific cell cycle phase [33]. Since G1 phase cells have a single copy of the genome, thus having the lowest amounts of DNA, whereas S phase cells are actively involved in DNA replication and will have increased amounts of DNA, and G2/M cells have two times the nuclear DNA of G1 phase cells, flow cytometric analyses of propidium iodide-stained nuclei can effectively differentiate G0/G1, S, and G2/M populations.
Before treating C. albicans SC5314 cells with SH1009 or cytochalasin D (CytoD) as a positive control for cell cycle arrest, a sample of early exponential culture was harvested and processed as described in the Materials and Methods to ensure cells were actively growing. The histogram in Fig 8A. indicates a rapidly dividing culture by having a significant fraction of cells (~70%) in the S phase. The histogram depicts an optimal flow-cytometric distribution for an actively growing yeast culture in rich media during early exponential phase, and was compatible with previous observations [33, 34]. After an additional 3 h of incubation without treatment, the untreated C. albicans SC5314 cells were still dividing with approximately similar cell cycle-distributed phases before three additional hours of incubation (Fig 8B). However, contrary to the cell cycle progression in untreated samples, the cell cycle distribution for C. albicans SC5314 cells that were incubated for 3 h with the SH1009 IC50 concentration were distinctly perturbed (Fig 8C). There was a significant decrease in the fraction of cells in the S phase (46.88% compared to 70% in untreated cells, P ≤ 0.001) and a significant increase in the proportion of cells in the G1 phase (46% compared to 16.53%, P ≤ 0.0001), indicating accumulating cells in G1 phase (Fig 8E). As expected, 3 h of treatment with CytoD, an anticancer drug that inhibits the assembly and disassembly of actin subunits, led to a delay in the progression of G1 phase [35]. The DNA histogram revealed fewer cells in S phase and more in the G1 phase (52.68%, 39.17.0% and 2.13% in G0/G1, S and G2/M phases, respectively) (Fig 8D).
The accumulation of SH1009-treated cells in G1 phase implies that aurone SH1009 arrests cell cycle progression, supporting the chemical-genetic interaction results in which the most significantly responsive mutants to SH1009 essentially possess deletions for cell cycle-encoding proteins. The chemical-genetic analysis showed that homozygous deletion strain CLN2Δ, which lacks the G1 cyclin gene, was resistant to SH1009 (Fig 7B), and heterozygous deletion strain CDC28Δ, cyclin-dependent kinase (CDK), was sensitive to SH1009 (Fig 7A). In S. cerevisiae, cyclin Cln2 activates CdcC28 in late G1 phase, resulting in regulation of actin cytoskeleton polarization, which is crucial for bud emergence and G1 to S phase transition during cell cycle progression [36]. Therefore, we next sought to determine whether aurone SH1009 affected the disruption of the actin cytoskeleton.
Aurone SH1009 perturbs actin cytoskeleton dynamics in C. albicans
The chemical-genetic interaction analysis indicated that cell-cycle-dependent organization of the actin cytoskeleton and endocytosis are targets of aurone SH1009 due to highly significant growth perturbation of 18 heterozygous and homozygous deletion mutants annotated as involved in actin cytoskeleton organization and endocytosis (Fig 6B). Fluorescent staining of actin and confocal laser scanning microscopy were used to visualize the effects of SH1009 on polarization of the actin cytoskeleton. SH1009-treated and untreated C. albicans SC5314 cells along with the S. cerevisiae CDC42Δ heterozygous deletion strain with a deletion for an actin-regulatory gene were examined microscopically for phenotypes indicative of disruption of actin distribution. In late G1 phase of the cell cycle of budding yeast, the actin patches should normally be assembled in the tip of the small-growing bud with clear actin cables that are polarized from the bud tip towards the mother cell. Actin cables serve as tracks to transport mRNAs, proteins, mitochondria, and ribosomes from the mother cell to the daughter cell. Actin patches are clustered in the secretion and endocytosis sites for critical roles in plasma membrane invagination during endocytosis [37] and cell wall remodeling through orienting the polarized secretion of cell wall constituents and enzymes towards the actin patches, promoting bud emergence [38, 39].
The images in Fig 9 depict typical actin polarization for untreated C. albicans SC5314 cells during exponential phase in which ~76.47% of cells were observed with red-fluorescent actin patches concentrated in the buds along with polarized-apical cell growth. However, after 3 h of SH1009 treatment, C. albicans cells appeared abnormally round and enlarged compared to untreated cells, indicating depolarized growth (Fig 9, DIC panel). In addition, only ~34.17% of SH1009-treated cells retained actin patches in the buds while the remaining had actin patches that were scattered randomly in both the mother cell and buds, signifying actin depolarization. In addition, SH1009-treated cells displayed a distinctively distorted actin assembly resulting in a considerable number of large aggregates of the actin (Fig 9, Actin panel). This phenotype was previously associated with endocytic mutants such as, but not limited to, ARK1Δ [40], END3Δ [41], and RPS5Δ [42]. In growth curve studies with SH1009-individually treated deletion mutants, these strains were either sensitive (ARK1Δ and RPS5Δ) or resistant (END3Δ) to SH1009 (Fig 7A and 7B) and were also identified in our chemical-genetic interaction analysis (Fig 6B). A previous study that investigated the composition of the actin clumps using immune-electron microscopy detected an accumulation of endocytic vesicles surrounded with actin filaments and a mixture of actin patches and endocytic proteins, suggesting a failure to mature the endocytic vesicles properly as a consequence of the inability to disassemble the actin-associated endocytic complexes [43].
The actin distribution in the S. cerevisiae CDC42Δ deletion mutant was also investigated due to its sensitivity to SH1009 (Fig 7A). In the majority of S. cerevisiae CDC42Δ mutant cells, there was an accumulation of large, round, unbudded cells with distributed actin patches (Fig 9, DIC panel), which supports previous reports [36]. The depolarized growth as well as the abnormal distribution of the actin patches in SH1009-treated C. albicans cells support inhibition of cell cycle arrest by preventing G1 phase progression (Fig 8D). In yeast cells, it has been evidenced in several studies that failure in rearrangement of actin patches at the bud site in G1 phase will prevent the emergence of the bud and arrest the cell cycle at G1 phase [44, 45]. These findings are compatible with previous observations that have established a definitive link between intact actin cytoskeleton organization and cell cycle progression [35].
Aurone SH1009 alters C. albicans expression of genes involved in the cell cycle, actin polarization, and endocytosis
To study changes in gene expression in response to aurone SH1009, quantitative reverse transcription-PCR was employed to confirm mRNA abundance of a set of C. albicans SC5314 genes homologous to the S. cerevisiae-S288C genes in the heterozygous and homozygous deletion mutants identified as significantly enriched in pathways from the chemical-genetic interaction analysis (Fig 5). The cell cycle-associated genes were the most responsive to aurone SH1009. Successful cell cycle progression requires expression of checkpoint genes that guarantee sequential execution of certain cellular processes. For instance, before chromosome segregation, chromosomal DNA must first be replicated. These processes are ordered by activation and inactivation of (CDK) Cdc28, an ATP-binding protein, which complexes with activating subunits called cyclins [44]. HGC1, is a C. albicans gene homologous to S. cerevisiae G1 cyclin gene CLN2, encoding a G1 cyclin that binds with Cdc28 to form a cyclin-CDK complex with a central role in G1/S transition during the cell cycle [46, 47]. SH1009-treated C. albicans SC5314 results in a ~2.3 and 4.9-fold respective downregulation of CDC28 and HGC1 gene expression after 1.5 h of treatment (Fig 10), supporting a role for SH1009 in arresting the C. albicans cell cycle in G1 phase (Fig 8C).
Binding of Cdc28 to Hgc1 is an important regulator step for controlling cell cycle progression. This complex regulates polarity by phosphorylating and blocking Rga2 (a negative regulator of Cdc42), thus activating Cdc42 (GTP-binding protein), leading to sustained actin polarization and hyphal growth [48]. Failure in the relocation of Cdc42 by the Cdc28-Hgc1 complex to either the bud site or hyphal tip in late G1 phase results in a haphazard distribution of actin cytoskeleton and subsequently accumulation of unbudded, enlarged cells [36, 45]. The S. cerevisiae CDC42Δ deletion mutant was hypersensitive to aurone SH1009 treatment (>70% inhibition) at a concentration of 16 μM (Fig 7A). In SH1009-treated C. albicans SC5314, CDC42 expression was downregulated two-fold (Fig 10), which would confirm the downstream effect of SH1009 on expression of CDC28 and HGC1.
Cdc28 also controls the cell cycle by regulation of DNA replication. The MCM2 gene encodes an ATP-binding protein that is a part of the pre-replicative complex known as mini-chromosome maintenance, consisting of proteins Mcm2-7, which acts as a helicase to unwind DNA and initiate DNA replication [49]. In S. cerevisiae, after assembly of the MCM2-7 complex in G1 phase, MCM2-7 requires phosphorylation at the end of G1 by a Cdc28-kinase complex in order to recruit DNA polymerase and initiate DNA replication [50]. SH1009 treatment resulted in a significant five-fold decrease in MCM2 by 1.5 h after treatment (Fig 10), which could represent an additional downstream effect of SH1009 on Cdc28. Mmc2 is also an ATP-binding protein and, along with 12 other proteins enriched by chemical-genetic interaction analysis, contains the widely-distributed P-loop motif, supporting nucleotide-binding proteins as potential targets for SH1009 [51].
In addition to CDC28, other genes that encode for nucleotide-binding proteins with roles in signal transduction pathways were identified in deletion mutants with differential sensitivities to aurone SH1009. The C. albicans CDC25 and RAS1 genes were downregulated by an approximate three-fold and two-fold, respectively (Fig 10). CDC25 encodes a GTP-binding protein that acts as a guanine nucleotide exchange factor and is the upstream activator of the Ras/cAMP signaling pathway, responsible for activating the GTP-binding protein, Ras1, which then activates the cAMP synthesis required for cell cycle progression [52, 53]. Additionally, Ras1 activation of cAMP synthesis leads to maintenance of the hyphal growth, a virulence factor in systemic candidiasis, by regulating the G1 cyclin Hgc1 which binds to Cdc28 to localize Cdc42 to the hyphal tip [46]. As evidenced by a previous study, SH1009 was the only aurone compound that disrupted the biofilm formation [21], which would require inhibiting hyphal growth as a key component of biofilm formation.
In contrast to decreased expression of RAS1 transcripts, RAS2 transcripts were upregulated two-fold after 3 h of SH1009 treatment (Fig 10). RAS2 encodes a GTPase in C. albicans, and a previous study has shown cAMP levels in a RAS1Δ mutant declined 20-fold yet increased by ~10% in a RAS2Δ mutant, indicating an antagonistic activity of C. albicans Ras2 on the cAMP levels through an unknown mechanism. In S. cerevisiae, both Ras1 and Ras2 are required for activating cAMP synthesis. However, unlike in S. cerevisiae in which Ras1 and Ras2 have been well studied and share sequence homology, C. albicans RAS2 has poor sequence homology with S. cerevisiae RAS sequences [54]. Although the response of RAS2 to SH1009 is unclear due to the lack of studies on C. albicans Ras2 activity or its interaction with other members of the Ras/cAMP signaling pathway, the upregulation of transcripts may be a response to downregulation of upstream molecules affected by SH1009.
Another transcriptional upregulation in SH1009-treated C. albicans was increased expression of TSC11 after 1.5 h of aurone treatment (Fig 10). TSC11 encodes a GTP-binding protein characterized as regulating actin cytoskeletal dynamics during polarized growth and endocytosis [55, 56]. The increased expression of TSC11 could be a detoxification mechanism for C. albicans to counter the negative effects of SH1009 on actin cytoskeletal dynamics and endocytosis. Previous studies reported phenotypes of TSC11Δ mutants (also known as AVO3 or RICTOR in mammals) as having abnormal actin polarization [57], reduced endocytosis rate [58], and arrested cell cycle progression [59]. Another gene involved in endocytosis exhibiting altered expression in response to SH1009 was ERG3 which was downregulated three-fold (Fig 10). ERG3 is involved in the ergosterol biosynthesis pathway, and gene deletion has been documented as attenuating the endocytosis rate [60]. Reduced expression of ERG3 could also be attributed to the upstream effects of SH1009 on the RAS/cAMP signaling pathway that positively regulates the expression of the ERG gene family [61].
Gene expression of RPS4a, which encodes a 40S ribosomal protein, decreased two-fold with SH1009 treatment (Fig 10). A genome-wide study of haploinsufficient C. alibcans deletion mutants that were fractionated for abnormal cell size revealed genes associated with ribosome biogenesis and cell cycle, suggesting a correlation between ribosome biogenesis rate and size-dependent cell cycle progression [62]. In S. cerevisiae, repressing the synthesis of this conserved gene, RPS4a, results in arresting G1 phase or a significantly prolonged G1 phase, along with the phenotype of increased cell size [63], which could explain the enlarged cell phenotype observed microscopically in SH1009-treated C. albicans cells (Fig 9).
Discussion
In this study, we investigated in vitro antifungal activity of the bioactive compound, aurone SH1009, against widely used CLSI reference resistance tester strains and additional strains of Candida spp. using the CLSI broth microdilution method [22]. Aurone SH1009 was not only inhibitory for susceptible strains, but was also inhibitory for resistant-clinical isolates. The clinical isolates, C. albicans M4: Gu2 and C. albicans M6: F1 were isolated initially as fluconazole-susceptible strains from two HIV-infected patients who were suffering from recurrent oropharyngeal candidiasis. After a two-year period of fluconazole treatment, the last isolates in each series were fluconazole-resistant counterparts known as C. albicans M5: Gu5 and C. albicans M7: F5 [64]. C. albicans M5 and M7 have acquired gain of function mutations, leading to homozygous alleles for the transcriptional factors MRR1 and TAC1 which have been documented for overexpressing the CDR1/2 (Candida Drug Resistance) and MDR1 (Multi-Drug Resistance) genes for azole-efflux pumps, respectively [65]. Additionally, C. albicans M2: ScTAC1R34A and C. albicans M3: ScMRR1R34A are fluconazole-resistant strains that have been mutated in C. albicans SC5314 background to encode MRR1P683S and TAC1G980E homozygous activating alleles, respectively. As a consequence of these mutations, the expression of efflux pumps are induced, causing a reduction in intracellular accumulation of azoles and ultimately high fluconazole resistance [66, 67].
Loss-of-heterozygosity events as genetic alterations are commonly associated with fluconazole resistance, reflecting the capacity of C. albicans to generate adaptive-homozygous mutations, which is attributed to the extraordinary plasticity of the C. albicans genome [68, 69]. Aurone SH1009 exhibited considerably low IC50 values in the inhibition of these fluconazole-resistant isolates, indicating that the highly expressed efflux pumps could be modulated by aurone SH1009. The efflux pumps CDR1/2 are mainly ATP-binding cassette (ABC) transporter proteins that harness the energy from ATP hydrolysis in order to extrude drug out of the yeast cell [7]. In fact, as consistent with previous studies on cancer cells, aurones selectively inhibited the pumping action of these transporters by binding with the C-terminal nucleotide-binding domain of P-glycoprotein which belongs to the ATP-binding cassette (ABC) superfamily [70, 71]. The modulatory effect of aurones on drug-efflux pumps allows the intracellular accumulation of the compound and its associated toxicity, which reverses the resistance mechanism and renders the resistant cells as susceptible again, making these bioactive compounds attractive candidates as reversal agents to control drug resistance [72].
Multidrug resistant isolate C. albicans ATCC 64124 is known as Darlington’s strain, the name of the patient suffering from chronic mucocutaneous candidiasis [73]. Darlington strain, 64124, is defined as resistant to all antifungal agents except 5-fluorocytosine at high concentration [74]. The resistance mechanisms toward multiple classes of antifungal agents of Darlington’s strain have not been studied, but evidence of the evolution of multidrug resistance has been published recently, demonstrating the capacity of clinical C. albicans isolates to adapt during long-term antifungal treatment [75]. Multidrug resistance is recognized as a multistep process resulting from a gradual accumulation of mutations during long-term antifungal treatment. Briefly, the resistance began with fluconazole treatment that causes an activation mutation in TAC1, followed by caspofungin exposure that leads to a mutation in FKS1. Then, amphotericin B treatment causes a loss of function mutation in ERG2. Consequently, the final strain was a multidrug resistant strain that acquired all the three mutations and exhibited resistance to three classes of the standard antifungal drugs [7].
Unlike C. albicans that acquires the resistance, C. glabrata is resistant to azole drugs due to ~67 gain of function mutations in the transcriptional factor CgPDR1 gene, resulting in overexpression of CgCDR1/2-encoded efflux pumps [76]. This could account for the dramatic increase in C. glabrata infections since the introduction of azole drugs in the 1980s, identifying C. glabrata as the second most common non-albicans cause of candidiasis. Drug resistance of C. glabrata has not only been found for azoles, but also for all known antifungal drugs, posing a serious challenge to antifungal therapy for the organism [77]. Our results indicated that aurone SH1009 effectively inhibited Candida spp. growth and exhibited acceptable pharmacodynamics properties (time- and concentration-dependent killing) that were similar to fluconazole, as a valuable estimation for a prospective antimicrobial agent [78].
After establishing SH1009 antifungal activity, we proceeded to characterize the antifungal mechanism. Studying the mode of action of a new drug is a critical prerequisite to improve therapeutic action and avoid unwanted side effects, but it has been a challenge in drug discovery due to the high consumption of cost and time. However, since the first yeast deletion collection was created in 2002 [79], chemogenomic profiling has been proven as a powerful system for characterizing the mode of action and molecular targets of several drugs [80], as well as thousands of additional bioactive molecules [81, 82]. A considerable number of studies have been recently reviewed the broad spectrum of biological activities of aurone compounds as anticancer agents identified using phenotype-based approaches, exploiting the existing knowledge of a given cellular process inhibited by aurones, including inhibition of CDKs, DNA scissoring, histone deacetylase, topoisomerase, ATP-binding cassette transporter, and tubulin polymerization [83]. In contrast, here, a genetic-based approach was applied to comprehensively reveal all the biological effects of aurone treatment as a potential antimicrobial. The chemogenomic approach that was used in this study is a reverse-genetics technique that allows for comprehensive identification of gene products that functionally interact with bioactive molecules by exposing previously constructed libraries of S. cerevisiae mutants to the condition of interest, which in this study was exposure to the novel aurone SH1009. These yeast mutants harbor deletions of the majority (96%) of the yeast genome with only one mutant per deletion-gene. Due to nearly complete coverage of the yeast genome, this approach guarantees unbiased results because it completely maps all targets (biological responses) simultaneously under the term of interconnected pathways of eukaryotic cell [84].
In this study, chemogenomic profiling of S. cerevisiae mutants with aurone SH1009 identified roles for genes involved in cell cycle, DNA replication, cell division, actin cytoskeleton, and endocytosis. Enrichment analysis for biological responses in heterozygous and homozygous mutants was first performed separately to confirm significant responses in each pool. The responses were then combined to map the functional annotation network of KEGG pathways, biological responses, and molecular functions simultaneously using GO enrichment analysis to find that the top responsive heterozygous and homozygous mutants clustered for actin cytoskeleton and endocytosis processes. GO enrichment analysis was also used to map the cell cycle pathway as the most significantly enriched pathway, also encompassing the other enriched biological processes and revealing nucleotide-binding proteins as the molecular function of the most of significantly enriched genes.
Nucleotide-binding proteins play a central role in a variety of pivotal cellular processes, including cell signaling, proliferation, cytoskeletal assembly, protein synthesis, and apoptosis [85]. The low number of GTP-binding proteins that were identified in this study could be attributed to a relatively low abundance of GTP-binding proteins in nature compared to ATP-binding proteins [86]. Nevertheless, because purine nucleoside triphosphates share a similar structure, some ATP-binding proteins may utilize GTP as the phosphate donor and vice versa [87]. In this respect, aurone SH1009 might possess an affinity for the same binding site in different but related nucleotide-binding proteins, hence, affecting a broader spectrum of cellular processes. The functional enrichment of genes encoding nucleotide-binding proteins during analysis of the deletion pool of essential and non-essential genes suggests that aurone SH1009 potentially interacts with nucleotide-binding proteins, leading to a series of cellular defects impacting the cell cycle, actin cytoskeleton, and endocytosis. The benzofuranone ring of SH1009 may dictate the biological activity due to mimicry of adenine, subsequently, inhibiting the activity of ATP-dependent proteins [88].
Our chemical genetic interaction and enrichment analysis also suggest that treatment with aurone SH1009 targets nucleotide-binding proteins, producing alterations in cell cycle, actin cytoskeleton organization, and endocytosis in S. cerevisiae. These phenotypic changes in response to SH1009 were confirmed in C. albicans SC5314 using cell sorting analysis by flow cytometry, which delineated increased populations of cells in G1 phase after 3 h of aurone SH1009 treatment. In S. cerevisiae, heterozygous mutants with gene deletions for proteins affecting the cell division cycle (CDC5Δ, CDC34Δ, CDC42Δ, CDC25Δ, and CDC13Δ) that sensitively responded to aurone SH1009 are all controlled by the master regulator of the cell cycle, (CDK) Cdc28, a catalytic kinase that couples with activating subunits such as Clb3 and Cln2 [89], the genes for which both responded to SH1009 as homozygous-resistant mutants. Also, the accumulation of SH1009-treated C. albicans cells in G1 phase implies that the cells are in cell cycle arrest and do not yet engage in DNA synthesis, also supporting the chemical-genetic interaction results in which DNA replication mutants POL3Δ, MCM2Δ, MCM7Δ, RFC1Δ, RFC3Δ, and ORC4Δ from the HIP profile were sensitive to SH1009. In C. albicans, gene expression during different cell cycle stages has been previously characterized using microarray technology [90]. This study reported expression of several genes that are required for G1/S transition and some of these genes are identified as sensitive or resistant mutants from our chemical-genetic profiles, including POL3 for DNA polymerase subunit, RFC1 and RFC2 for DNA replication factor elements, TOF1 gene for a DNA replication checkpoint, CSM3 gene for accurate chromosome segregation, SWI4 transcription factor for activation of G1 transition (Fig 6).
Our results are also in agreement with a prior study reporting that aurones arrest the cell cycle via targeting (CDK) protein Cdc28, making aurones a promising candidate for cancer therapy [12]. Among different aurone derivatives synthesized as potential anticancer agents, aurone compounds that possess a methoxy group at position 4 on ring B, like aurone SH1009, were found to produce an enhanced antiproliferative activity by significantly arresting the cell cycle at G2/M phase [83]. In our study, fluorescent imaging was performed to detect actin cytoskeleton depolarization as a consequence of aurone SH1009 perturbation of rearrangement of actin patches in the bud site at the end of G1 phase, leading to cell cycle arrest at G1 phase. These phenotypic findings support the chemogenomic responses of S. cerevisiae genes CDC28Δ, CLN2Δ, and CDC42Δ to SH1009, because it has been well characterized that the G1 cyclin Cln2 binds to Cdc28 to localize Cdc42 to promote actin polarization and G1/S phase transition [36].
Chemogenomic profiling does have limitations; for instance, cellular mechanisms affected by SH1009 were identified using growth as a sole endpoint measurement, potentially resulting in selection for false positive sensitivities due to slow growth of some deletion mutants or, in the worst scenario, loss of indispensable mutants. Another limitation is that responses in S. cerevisiae differ from responses in C. albicans because of the ability of C. albicans to produce true hyphae and differences in sensitivity and resistance to antifungals. To confirm specific gene responses in C. albicans, relative levels of mRNA of genes identified as differentially responsive were quantified in C. albicans SC5314. A sizable subset of genes that function in cell cycle progression was differentially regulated in response to SH1009, leading to a potentially hierarchal pathway beginning with Cdc25 and Ras1 that regulate G1 cyclin Hgc1, which binds to Cdc28 in order to localize Cdc42, leading to actin polarization and cell cycle progression. Cdc25 and Ras1 regulate the RAS/cAMP signaling pathway which controls Erg3 maintenance of normal endocytosis rates. The transcriptional changes in the genes for these molecules support the biological responses observed in phenotypic studies.
Two genes, CBK1 and RSP5, identified by chemogenomic profiling in S. cerevisiae, did not demonstrate significant differential expression to SH1009 treatment in C. albicans in spite of their respective documented roles in polarized growth [91] and organized actin patches with normal endocytosis [42]. CBK1Δ mutant was among the top 20 sensitive heterozygous deletion mutants (S2 Table), and the RSP5Δ mutant was confirmed using growth curves for its sensitivity to SH1009 (Fig 7A). Failure to detect changes in gene expression could be related to the timing of expression or discordance between regulation and gene expression. In S. cerevisiae, the TSC11Δ mutant was hypersensitive to SH1009, however, this gene was transcriptionally upregulated in C albicans, which could indicate a difference in survival responses between the two yeasts. This highlights that although this study has validated chemical-genetic interaction using S. cerevisiae mutant collections to characterize the mode of action for a novel bioactive aurone compound against C. albicans, chemical-genetic interaction studies could also be used to elucidate potential differences in drug responses between organisms.
Conclusion
In summary, we investigated aurone SH1009 for antifungal activity, human cell line cytotoxicity, and pharmacodynamics properties. Aurone SH1009 exhibited promising selectively fungistatic-inhibitory activity against Candida spp., particularly resistant isolates, highlighting the possibility of a broad-spectrum property for aurone SH1009 against other pathogenic fungi. We also provided the first, comprehensive genome-wide study of haploinsufficiency and homozygous screening in S. cerevisiae for a bioactive aurone compound. Chemical-genetic interaction analysis predicted that SH1009 targets cell cycle-dependent organization of the actin cytoskeleton and endocytosis, suggesting a novel mode of action for this aurone compound. Phenotypic studies for these significantly enriched biological responses were completed in C. albicans demonstrating G1 phase-arrested cells with abnormal, depolarized actin cytoskeleton. Differential expression of genes, identified by chemical-genetic interaction, in response of SH1009 in C. albicans coincided with enrichment of cells in G1 phase and actin depolarization producing a model that offers a significant improvement in understanding the mechanism for toxicity. The current study provides experimental frameworks for future mechanistic studies that could be used to investigate the bioactivity of a large number of aurone compounds using chemical-genetic interaction and phenotypic-based methods.
Materials and methods
Materials and reagents
Candida strains (listed in S1 Table) were provided by the Dr. P. David Rogers lab of the University of Tennessee Health Science Center, Memphis, TN. Human cell lines; THP-1 (ATCC, TIB-202), HepG2 (ATCC, HB-8065), and A549 (ATCC, CCL-185) were purchased from the American Type Culture Collection (Manassas, VA, USA). Dulbecco′s Modified Eagle′s culture medium (DMEM), 1X trypsin-EDTA solution, fetal bovine serum (FBS), 100X penicillin/streptomycin solution, Amphotericin B, caspofungin, fluconazole, itraconazole, 5-fluorocytosine, YPD agar and broth, cytochalasin D, 3-(N-morpholino) propanesulfonic acid (MOPS) buffer, and methyl methanesulfonate (MMS) were all purchased form Sigma-Aldrich (St. Louis, MO, USA). RPMI-1640 medium was purchased from Corning Incorporated (Coring, NY, USA). Geneticin selective antibiotic (G-148 sulfate), phosphate buffered saline (PBS), rhodamine phalloidin, propidium iodide, PrestoBlue, and RNaseA enzyme were purchased from Life Technologies Corporation (Carlsbad, CA, USA). 4% Paraformaldehyde was purchased from Alfa Aesar (Ward Hill, MA, USA) and zymolase 20T was purchased from MP Biomedicals, LLC (Solon, OH, USA). The yeast deletion collections (~1,056 heterozygous mutants and ~4,320 homozygous mutants) were purchased from GE Healthcare Life Sciences (Pittsburg, PA, USA) and ThermoFisher Scientific (Waltham, MA, USA), respectively.
Antifungal susceptibility testing
Preparation of stock solutions
Aurone SH1009 was synthesized as described in supporting information (S1 File). The powder of SH1009 was dissolved in dimethyl sulfoxide (DMSO) to a high concentration of 20 mM. Two-fold serial dilutions were prepared at concentrations from 200 to 3.125 μM using RPMI–1640 medium that was buffered previously to pH 7.0 with MOPS and sterilized by filtration. Using 96-wells microtiter plates, 100 μL of each SH1009 concentration was added to respective wells with four replicates for each concentration. Amphotericin B and caspofungin were used as positive controls to ensure 100% growth inhibition of the yeast at concentrations of 16 μg/mL and 8 μg/mL, respectively. In addition, fluconazole, itraconasole, and 5-fluorocytosine were prepared following manufacturer instructions, serially diluted two-fold according to CLSI protocol concentrations [22], and used as a reference for some of the strains to confirm their resistance profiles.
Preparation of inocula
Strains for this study (S1 Table) were cultured on YPD agar and incubated at 35°C for 24 h. The inoculum suspensions for each strain was prepared according to the CLSI broth microdilution protocol for antifungal susceptibility testing of yeasts [22]. The suspension of 5–6 colonies was vortexed in approximately 4 mL of sterile saline solution (0.85% NaCl) and adjusted spectrophotometrically to optical densities at a 530 nm wavelength (OD530) that ranged from 0.12 to 0.15. The inocula were then diluted 1:1000 in RPMI 1640 medium resulting in 1×103 to 5×103 CFU/mL working concentration. Volumes of 100 μL of each inoculum for each strain were added to the wells of its respective plate. For each isolate, there were drug-free wells and media-control wells with and without 1% DMSO to detect any contamination in the media and for use as an optical blank for optical density and fluorescence measurements.
Inhibition assay
After 24 h of incubation at 35°C, 20 μL of PrestoBlue reagent were added to each microtiter well to a final concentration of 10% after which plates were incubated at 35°C for an additional 60–70 min [21]. The minimum inhibitory concentrations (MICs) were defined as the concentration of aurone 1009 that reduced the growth by 90%. The MIC value was determined quantitatively by measuring the fluorescence that results from reducing blue-nonfluorescent resazurin, to red-fluorescent resorufin as a result of metabolic activity of the active cells at 560 nm excitation and 590 nm emission with a SpectraMax M5e spectrophotometer (Molecular Devices, LLC, USA). The percentages of yeast growth were calculated by comparing the fluorescent readings of the drug-containing wells with that of the drug-free wells to calculate MIC90. The assay for each strain was performed in duplicate. The IC50 values were calculated using GraphPad Prism (GraphPad Software, USA).
Cell viability assay
Using 24-well plates, C. albicans SC5314 cells were cultured and treated with 200 μM SH1009 in 2 mL of RPMI 1640 medium as described above. After incubation at 35°C for 48 h, 1 mL of SH1009-treated cells, untreated cells, and previously prepared isopropanol-killed cells were washed with PBS. The samples were then diluted to approximately 1 x 106 cells/mL in PBS and stained according the live/dead Fungalight Yeast Viability Kit protocol (ThermoFisher, Waltham, MA, USA). Stained cells were investigated for their viability using flow cytometry (Guava Millipore, Burlington, MA, USA) with Guava InCyte Software. To distinguish between live and dead cells in a dot plot, gating the yeast population was determined based on red and green fluorescence levels of isopropanol-killed cells.
Growth rate assay
In 100-well Bioscreen honeycomb plates, C. albicans M1:SC5314 cells were cultured and treated with aurone SH1009 (3.125–200 μM) in RPMI 1640 medium as described above. Plates were loaded into the Bioscreen C instrument with Bioscreen software (Growth Curves USA, Piscataway, NJ, USA) at a temperature of 35°C with continuous shaking and 30 min interval measurements at an OD530 for 40 h. Growth curves were used to compare growth of untreated-C. albicans cells with SH1009-treated cells at different aurone concentrations. Validation of cellular sensitivity and resistance responses of S. cerevisiae mutants to SH1009 were also performed with the Bioscreen C instrument as the same manner for C. albicans except the media was YPD broth and incubation was at 30°C.
Time-kill assay
A previously described and evaluated antifungal time-kill method was utilized to evaluate the fungicidal activity of aurone SH1009 [92]. An initial inoculum of C. albicans SC5314 ranging from 0.5–1.0 ×105 CFU/mL was treated with an aurone concentration approximately five-fold higher (500 μM) than the dilution producing the IC50 of SH1009 (16.25 μM) in the two-fold dilution series. Fluconazole (16 μg/mL) was used as fungistatic control, and Amphotericin B (8 μg/mL) was used as fungicidal control. After incubation at 35°C, 10 μL of each treatment was spread onto YPD agar after 0, 6, 12, 24, and 30 h of treatment and plates were incubated at 35°C for 24 h to determine viable cell numbers. The fungicidal activity was determined as ≥ 3 –log10 which is equivalent to 99.9% reduction in CFU/mL from the working concentration 0.25–0.5×105 CFU/mL.
Cytotoxicity assay
The A549 human lung carcinoma epithelial cell line and human monocytic THP-1 cell line were grown in RPMI-1640 culture medium, while HepG2 human liver carcinoma epithelial cells were grown in DMEM culture medium. Both media were supplemented with 10% FBS and 1% penicillin-streptomycin antibiotics. After maintaining the cell growth at 37°C with 5% CO2 in a humidified incubator until reaching 90% confluency, the A549 and HepG2 cells were trypsinized with 1X trypsin-EDTA and resuspended in fresh medium. The cells were seeded into 96-well microtiter plates at a density of 10,000 viable cells/well and grown overnight, while the suspension THP-1 cells were seeded directly into 96-well microtiter plates at the same density before treatment. The final concentrations of aurone SH1009 were prepared in two-fold serial dilutions (3.125 μM– 200 μM) as described above for the antifungal susceptibility assay. The media containing A549 or HepG2 cells were then replaced after overnight incubation by 200μL of fresh culture media containing the final concentrations of SH1009. The cells were then incubated for additional 24h at 37°C with 5% CO2 in a humidified incubator. To evaluate cell viability, each well was treated with 20 μL of PrestoBlue for 3–6h. Metabolically active cells converted the blue non-fluorescent dye resazurin to the pink fluorescent dye resorufin, which can be measured by plate reader as described above in antifungal inhibition assay. Triton X-100 (1%, v/v) was used as a positive control to give a complete loss of cell viability. Percentages of cell viability were calculated as follows: [(negative control value–treated value) × 100]/negative control value. The assay for each cell line was performed in triplicate. The CC50 values were calculated using GraphPad Prism (GraphPad Software, USA).
Chemogenomic profiling in Saccharomyces cerevisiae
Combination of individual mutants into a single pool
The yeast deletion collections were obtained as individual mutants in 96-well plates that had been stored at −80°C. Both the HIP deletion pool and HOP deletion pool were created separately as previously described [84, 93]. The 96-well plates of mutant stocks were thawed completely, after which a 96-well transfer pin was used to transfer a small volume of mutants to a Nunc Omni Tray containing YPD agar with geneticin antibiotic. Between transfers, the 96-well transfer pin was sterilized with ethanol and flamed three times. After growing the cells for 48 h at 30°C, the missing and slow growing mutants were recorded and two times the cell mass of these mutants were added separately. Working in a sterile hood, each tray was flooded with ~10 mL of YPD broth and all formed colonies were gently scraped by sterile cell spreader. The resuspended colonies were transferred into a sterile 1000 mL flask with a sterile stir bar. The suspension was mixed for 5 min on a stir plate to obtain an homogenized pool. The concentration of the freshly prepared pool was adjusted to 125–250 cells/mutant/μL by centrifugation at 500 × g. Once the concentration was adjusted, sterile glycerol was added to 15% (vol/vol) and 200 μL aliquots of the pool were stored in PCR strip tubes at −80°C.
Pooled competition with aurone SH1009
Before exposing the pooled deletion mutants to SH1009, the inhibitory concentration for roughly 20% of the S. cerevisiae S288C wild type parent strain of these mutants was determined (500 μM) as previously described [93]. Using 96-well plates, 8 wells were filled with 198 μL of SH1009 diluted in YPD broth at concentration 500 μM with no more than 1% of DMSO. For the positive control, 8 wells were filled with 198 μL of MMS diluted in YPD broth at concentration 0.01 μg/mL. For the negative control, 8 wells were filled with 198 μL of YPD broth with 1% DMSO only. Two aliquots that were prepared from the previous step of the HOP deletion pool (~4,320 mutants), representing nonessential genes, and HIP deletion pool (~1,056 mutants), representing essential genes, were thawed completely. A 2 μL volume at a concertation of 125–250 cells/mutant/μL of the non-essential deletion pool were added to every 12 wells containing SH1009, MMS, and 1% DMSO. In the same manner, 2 μL of the essential deletion pool at a concertation of 125–250 cells/mutant/μL were added to every remaining 12 wells of SH1009, MMS, and 1% DMSO. After a 48h incubation at 30°C, cells from each well were harvested independently by pipetting up and down and centrifuging at ~20,000 x g for 3 min. The supernatant was removed and the pellet was processed for genomic DNA extraction.
Construction of the DNA library
Pellets of 24 samples were resuspended individually in 125 μL of Zymolyase solution and incubated for 1 h at 37°C. The DNA was extracted from all 24 samples according to Maxwell 16 LEV Plant DNA Kit manual (Promega Corporation, Madison, WI, USA). To amplify the UPTAG unique 20 bp DNA barcodes as previously described [93], 24 PCR reactions were prepared independently, such that there was one PCR reaction for each sample in a total volume of 25 μL as follows: 21.5 μL of Taq mix, 0.5 μL of reverse common primer at 0.5 μM, 0.5 μL of indexed primer at 0.5 μM (for each sample, a distinct indexed primer was used, S5 Table), and 2.5 μL of genomic DNA at ~ 100 ng. PCR conditions were as follows: 5 min at 95°C for an initial denaturation, followed by 30 cycles of 1 min at 95°C, 30 s at 55 °C, 45 s at 68°C, then, 10 min at 68°C for a final extension. After PCR reactions, 25 μL of all PCR products were pooled together from individual PCR tubes into one tube library. This library was purified by separation on a 2% TAE agarose gel for 50 min at 120V. The desired band (267 bp) containing the amplified UPTAG DNA barcodes was cut and purified from the gel using a QIAGEN Quick Gel Extraction Kit (QIAGEN, Germantown, MD, USA). The library was diluted to 1:5,000, 1:10,000, and 1:20,000 and quantified with the KAPA Library Quantification Kit (KAPA Biosystems, Wilmington, MA, USA) and Bio-Rad CFX96 real-time PCR system (Hercules, CA, USA). After quantifying the correct concentration of the library, the library was prepared as a DNA template at the final concentration 15 nM with a 5% Spike-In of PhiX control according to the MiSeq System Denature and Dilute Library Guide (Illumina, San Diego, CA, USA). The Illumina MiSeq sequencer was used to run the DNA template for 1 × 50 cycles to yield a cluster density of 700–900 k/mm2.
Post-sequencing data analysis
The Illumina sequencer generated a Fastq file that was converted to Fasta file using a converter tool due to the ease of manipulating the Fasta file [78]. In order to process and analyze sequence reads, Perl scripts were created (https://github.com/fma3b/Barcode_Seq_Analysis) [94]. Two raw-count files for HIP-HOP profiles were imported into excel sheet (Microsoft Corporation, US) to normalize the absolute counts and calculate fitness scores, Z-scores, P-values, and FDR values following calculations previously reported [84]. The raw sequence reads were deposited in Sequence Read Archive (SRA) under project number PRJNA491750. Enrichment analysis of KEGG pathway and gene ontology (GO) analysis was conducted using hypergeometric testing through ClueGo software to find the significantly enriched KEGG/GO terms using GO categories in the Saccharomyces cerevisiae-S288C as a background. To visualize the interactive annotation network between significant genes, ClueGO and CluePedia apps [29, 31] along with Cytoscape (Cytoscape Consortium, USA) were used.
Flow cytometry
C. albicans SC5314 cells were grown at 30°C until reaching exponential phase, diluted to ~1.2–3×106 CFU/mL, and treated with SH1009 at the IC50 concentration followed by incubation for an additional 3 h. After harvesting by centrifugation, supernatants were discarded and pellets were washed with PBS and then fixed with cold 70% ethanol at -20°C for 2 h. Fixed cells were washed in PBS and resuspended in 500 uL of PBS containing 20 μg/mL RNase and incubated at 37°C for 2 h. 200 μL of PBS containing 20 μg/mL propidium iodide (PI) were added to the treated cells. Using the Millipore Guava flow cytometer, 5000 events were counted, and the fluorescent intensity of PI measured. After acquiring the data using Guava PCA-96 software, the data was gated to exclude debris or aggregates. Experiments were performed in triplicate using cytochalasin D, which has been reported to arrest the cell cycle [35], as a positive control.
Confocal microscopy
C. albicans SC5314 cells were grown at 35°C until reaching exponential phase, then diluted to ~ 1.2–3×106 CFU/mL. Aurone SH1009 was added at the IC50 concentration and cells were incubated for 3 h. Cells were fixed by addition of 4% paraformaldehyde and subsequently incubated for 2 h at room temperature. After pelleting and washing the cells with PBS, cells were incubated with 1% Triton-X100 for 1 h at room temperature. Rhodamine Phalloidin (RP) was added to the cells followed by incubation in the dark at 4°C for 1 h. After two PBS washes, cells were imaged using confocal microscopy (Zeiss, Thornwood, NJ, USA) at 60x magnification. Cells was assessed for actin distribution by considering that the actin is depolarized if more than five patches were observed in the mother cell [95]. Approximately 100 cells were counted per experiment in triplicate experiments. The S. cerevisiae CDC42Δ mutant was used as a positive control.
RT-qPCR
After growing C. albicans SC5314 cells at 30°C in YPD broth until reaching exponential phase, the culture was treated with aurone SH1009 at the IC50 concentration followed by incubation for additional 1.5, 3, and 6 h. After harvesting the cells by centrifugation, the RNA was extracted according to the instructions of the Maxwell 16 LEV Plant RNA Kit. Total RNA from treated and untreated samples were normalized to 1 μg. cDNA was constructed by following the manufacturer protocol of SuperScript IV First-Strand Synthesis System kit (ThermoFisher, Waltham, MA, USA) using 10 ng of RNA. The RT-qPCR was preformed using 2× iQ SYBR green supermix (Bio-Rad, Hercules, CA, USA) under the recommended cycle conditions. All reactions were performed in triplicate using listed primer pairs (S6 Table). Transcript levels were normalized to the expression level of the housekeeping gene GAPDH and compared to the untreated sample using ΔΔCT method [96].
Supporting information
S1 Table [docx]
Strains used in this study.
S2 Table [xlsx]
Haploinsufficiency data analysis.
S3 Table [xlsx]
Homozygous data analysis.
S4 Table [docx]
Enrichment analysis of HIP-HOP profiles.
S5 Table [docx]
Primers used in chemogenomic analysis.
S6 Table [docx]
Primers used in RT-qPCR.
S1 Fig [tif]
Growth rate of . mutants.
S1 File [docx]
Synthesis of SH1009.
Zdroje
1. Brown GD, Denning DW, Gow NAR, Levitz SM, Netea MG, White TC. Hidden killers: human fungal infections. Sci Transl Med. 2012;4(165):165rv13, 10 pp. doi: 10.1126/scitranslmed.3004404 23253612
2. Sardi JCO, Scorzoni L, Bernardi T, Fusco-Almeida AM, Giannini MJSM. Candida species: current epidemiology, pathogenicity, biofilm formation, natural antifungal products and new therapeutic options. J Med Microbiol. 2013;62(1):10–24.
3. Rajendra P. Candida albicans: Cellular and Molecular Biology. 2nd ed: Springer Nature; 2017. doi: 10.14715/cmb/2017.63.7.12 28838343
4. Kullberg BJ, Arendrup MC. Invasive candidiasis. N Engl J Med. 2015;373(15):1445–56. doi: 10.1056/NEJMra1315399 26444731
5. Liken HB, Kaufman DA. Candida. In: Cantey JB, editor. Neonatal Infections: Pathophysiology, Diagnosis, and Management. Cham: Springer International Publishing; 2018. p. 33–49.
6. Bondaryk M, KurzAtkowski W, Staniszewska M. Antifungal agents commonly used in the superficial and mucosal candidiasis treatment: mode of action and resistance development. Postepy Dermatol Alergol. 2013;30(5):293–301. doi: 10.5114/pdia.2013.38358 24353489
7. Prasad R, Nair R, Banerjee A. Emerging Mechanisms of Drug Resistance in Candida albicans. In: Sá-Correia I, editor. Yeasts in Biotechnology and Human Health: Physiological Genomic Approaches. Cham: Springer International Publishing; 2019. p. 135–53.
8. Wilson Leslie S, Reyes Carolina M, Stolpman M, Speckman J, Allen K, Beney J. The direct cost and incidence of systemic fungal infections. Value Health. 2002;5(1):26–34. doi: 10.1046/j.1524-4733.2002.51108.x 11873380
9. Wisplinghoff H, Ebbers J, Geurtz L, Stefanik D, Major Y, Edmond MB, et al. Nosocomial bloodstream infections due to Candida spp. in the USA: species distribution, clinical features and antifungal susceptibilities. Int J Antimicrob Agents. 2014;43(1):78–81. doi: 10.1016/j.ijantimicag.2013.09.005 24182454
10. Aqil F, Zahin M, Ahmad I, Owais M, Sajjad M, Khan A, et al. Antifungal Activity of Medicinal Plant Extractsand Phytocompounds: A Review: Springer-Verlag Berlin Heidelberg; 2010.
11. Mishra BB, Tiwari VK. Natural products: An evolving role in future drug discovery. European Journal of Medicinal Chemistry. 2011;46(10):4769–807. doi: 10.1016/j.ejmech.2011.07.057 21889825
12. Boumendjel A. Aurones: A subclass of flavones with promising biological potential. Current Medicinal Chemistry. 2003;10(23):2621–30. doi: 10.2174/0929867033456468 14529476
13. Morimoto M, Fukumoto H, Nozoe T, Hagiwara A, Komai K. Synthesis and insect antifeedant activity of aurones against Spodoptera litura larvae. Journal of Agricultural and Food Chemistry. 2007;55(3):700–5. doi: 10.1021/jf062562t 17263463
14. Kayser O, Kiderlen AF, Folkens U, Kokodziej H. In vitro leishmanicidal activity of aurones. Planta Medica. 1999;65(4):316–9. doi: 10.1055/s-1999-13993 10364835
15. Bandgar BP, Patil SA, Korbad BL, Biradar SC, Nile SN, Khobragade CN. Synthesis and biological evaluation of a novel series of 2,2-bisaminomethylated aurone analogs as anti-inflammatory and antimicrobial agents. Eur J Med Chem. 2010;45(7):3223–7. doi: 10.1016/j.ejmech.2010.03.045 20430485
16. Lawrence NJ, Rennison D, McGown AT, Hadfield JA. The total synthesis of an aurone isolated from Uvaria hamiltonii: Aurones and flavones as anticancer agents. Bioorganic & Medicinal Chemistry Letters. 2003;13(21):3759–63.
17. Huang W, Liu MZ, Li Y, Tan Y, Yang GF. Design, syntheses, and antitumor activity of novel chromone and aurone derivatives. Bioorganic & Medicinal Chemistry. 2007;15(15):5191–7.
18. Ferreira EO, Salvador MJ, Pral EMF, Alfieri SC, Ito IY, Dias DA. A new heptasubstituted (E)-aurone glucoside and other aromatic compounds of Gomphrena agrestis with biological activity. Zeitschrift Fur Naturforschung C-a Journal of Biosciences. 2004;59(7–8):499–505.
19. Pare PW, Dmitrieva N, Mabry TJ. PHYTOALEXIN AURONE INDUCED IN CEPHALOCEREUS-SENILIS LIQUID SUSPENSION-CULTURE. Phytochemistry. 1991;30(4):1133–5.
20. Lopez-Lazaro M, Martin-Cordero C, Toro MV, Ayuso MJ. Flavonoids as DNA topoisomerase I poisons. J Enzyme Inhib Med Chem. 2002;17(1):25–9. doi: 10.1080/14756360290011744 12365457
21. Sutton CL, Taylor ZE, Farone MB, Handy ST. Antifungal activity of substituted aurones. Bioorg Med Chem Lett. 2017;27(4):901–3. doi: 10.1016/j.bmcl.2017.01.012 28094180
22. CLSI. Reference method for broth dilution antifungal susceptibility testing of yeasts. CLSI standard M27 4th ed ed: Wayne, PA: Clinical and Laboratory Standards Institute; 2017.
23. Pfaller MA, Diekema DJ. Epidemiology of invasive candidiasis: a persistent public health problem. Clin Microbiol Rev. 2007;20(1):133–63. doi: 10.1128/CMR.00029-06 17223626
24. Silva S, Negri M, Henriques M, Oliveira R, Williams DW, Azeredo J. Candida glabrata, Candida parapsilosis and Candida tropicalis: biology, epidemiology, pathogenicity and antifungal resistance. FEMS Microbiol Rev. 2012;36(2):288–305. doi: 10.1111/j.1574-6976.2011.00278.x 21569057
25. Webb NE, Montefiori DC, Lee B. Dose-response curve slope helps predict therapeutic potency and breadth of HIV broadly neutralizing antibodies. Nat Commun. 2015;6:8443. doi: 10.1038/ncomms9443 26416571
26. Boncler M, Różalski M, Krajewska U, Podsędek A, Watala C. Comparison of PrestoBlue and MTT assays of cellular viability in the assessment of anti-proliferative effects of plant extracts on human endothelial cells. J Pharmacol Toxicol Methods. 2014;69(1):9–16. doi: 10.1016/j.vascn.2013.09.003 24103906
27. Hoon S, St. Onge RP, Giaever G, Nislow C. Yeast chemical genomics and drug discovery: an update. Trends Pharmacol Sci. 2008;29(10):499–504. doi: 10.1016/j.tips.2008.07.006 18755517
28. Chang M, Bellaoui M, Boone C, Brown GW. A genome-wide screen for methyl methanesulfonate-sensitive mutants reveals genes required for S phase progression in the presence of DNA damage. Proceedings of the National Academy of Sciences. 2002;99(26):16934–9.
29. Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, et al. ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics. 2009;25(8):1091–3. doi: 10.1093/bioinformatics/btp101 19237447
30. Harris MA, Clark J, Ireland A, Lomax J, Ashburner M, Foulger R, et al. The Gene Ontology (GO) database and informatics resource. Nucleic acids research. 2004;32(Database issue):D258–D61. doi: 10.1093/nar/gkh036 14681407
31. Bindea G, Galon J, Mlecnik B. CluePedia Cytoscape plugin: pathway insights using integrated experimental and in silico data. Bioinformatics. 2013;29(5):661–3. doi: 10.1093/bioinformatics/btt019 23325622
32. Ericson E, Hoon S, St. Onge RP, Giaever G, Nislow C. Exploring gene function and drug action using chemogenomic dosage assays. Methods Enzymol. 2010;470(Guide to Yeast Genetics (2nd Edition)):233–55.
33. Rosebrock AP. Analysis of the Budding Yeast Cell Cycle by Flow Cytometry. Cold Spring Harb Protoc. 2017;2017(1).
34. Slater ML, Sharrow SO, Gart JJ. Cell cycle of Saccharomycescerevisiae in populations growing at different rates. Proc Natl Acad Sci U S A. 1977;74(9):3850–4. doi: 10.1073/pnas.74.9.3850 333447
35. Heng YW, Koh CG. Actin cytoskeleton dynamics and the cell division cycle. Int J Biochem Cell Biol. 2010;42(10):1622–33. doi: 10.1016/j.biocel.2010.04.007 20412868
36. Park HO, Bi E. Central roles of small GTPases in the development of cell polarity in yeast and beyond. Microbiol Mol Biol Rev. 2007;71(1):48–96. doi: 10.1128/MMBR.00028-06 17347519
37. Mishra M, Huang J, Balasubramanian MK. The yeast actin cytoskeleton. FEMS Microbiology Reviews. 2014;38(2):213–27. doi: 10.1111/1574-6976.12064 24467403
38. Irazoqui JE, Lew DJ. Polarity establishment in yeast. Journal of Cell Science. 2004;117(11):2169–71.
39. Kopecká M, Yamaguchi M, Kawamoto S. Effects of the F-actin inhibitor latrunculin A on the budding yeast Saccharomyces cerevisiae. Microbiology. 2015;161(7):1348–55. doi: 10.1099/mic.0.000091 25858300
40. Ayscough KR. Coupling actin dynamics to the endocytic process in Saccharomyces cerevisiae. Protoplasma. 2005;226(1–2):81–8. doi: 10.1007/s00709-005-0107-5 16231104
41. Gourlay CW, Ayscough KR. Identification of an upstream regulatory pathway controlling actin-mediated apoptosis in yeast. J Cell Sci. 2005;118(Pt 10):2119–32. doi: 10.1242/jcs.02337 15855235
42. Kamińska J, Gajewska B, Hopper AK, Zoładek T. Rsp5p, a new link between the actin cytoskeleton and endocytosis in the yeast Saccharomyces cerevisiae. Mol Cell Biol. 2002;22(20):6946–8. doi: 10.1128/MCB.22.20.6946-6958.2002 12242276
43. Sekiya-Kawasaki M, Groen AC, Cope MJ, Kaksonen M, Watson HA, Zhang C, et al. Dynamic phosphoregulation of the cortical actin cytoskeleton and endocytic machinery revealed by real-time chemical genetic analysis. J Cell Biol. 2003;162(5):765–72. doi: 10.1083/jcb.200305077 12952930
44. Howell AS, Lew DJ. Morphogenesis and the cell cycle. Genetics. 2012;190(1):51–77. doi: 10.1534/genetics.111.128314 22219508
45. Ushinsky SC, Harcus D, Ash J, Dignard D, Marcil A, Morchhauser J, et al. CDC42 is required for polarized growth in human pathogen Candida albicans. Eukaryot Cell. 2002;1(1):95–104. doi: 10.1128/EC.1.1.95-104.2002 12455975
46. Zheng X, Wang Y. Hgc1, a novel hypha-specific G1 cyclin-related protein regulates Candida albicans hyphal morphogenesis. EMBO J. 2004;23(8):1845–56. doi: 10.1038/sj.emboj.7600195 15071502
47. Wang Y. Hgc1-Cdc28-how much does a single protein kinase do in the regulation of hyphal development in Candida albicans? J Microbiol. 2016;54(3):170–7. doi: 10.1007/s12275-016-5550-9 26920877
48. Zheng XD, Lee RT, Wang YM, Lin QS, Wang Y. Phosphorylation of Rga2, a Cdc42 GAP, by CDK/Hgc1 is crucial for Candida albicans hyphal growth. EMBO J. 2007;26(16):3760–9. doi: 10.1038/sj.emboj.7601814 17673907
49. Tye BK. MCM proteins in DNA replication. Annu Rev Biochem. 1999;68:649–86. doi: 10.1146/annurev.biochem.68.1.649 10872463
50. Enserink JM, Kolodner RD. An overview of Cdk1-controlled targets and processes. Cell Div. 2010;5:11. doi: 10.1186/1747-1028-5-11 20465793
51. Saraste M R. Sibbald P, Wittinghofer A. The P-loop—a common motif in ATP- and GTP-binding proteins. Trends in Biochemical Sciences. 1990;15(11):430–4. doi: 10.1016/0968-0004(90)90281-f 2126155
52. Toda T, Uno I, Ishikawa T, Powers S, Kataoka T, Broek D, et al. In yeast, RAS proteins are controlling elements of adenylate cyclase. Cell. 1985;40(1):27–36. doi: 10.1016/0092-8674(85)90305-8 2981630
53. Weeks G, Spiegelman GB. Roles played by Ras subfamily proteins in the cell and developmental biology of microorganisms. Cell Signal. 2003;15(10):901–9. doi: 10.1016/s0898-6568(03)00073-1 12873703
54. Zhu Y, Fang HM, Wang YM, Zeng GS, Zheng XD, Wang Y. Ras1 and Ras2 play antagonistic roles in regulating cellular cAMP level, stationary-phase entry and stress response in Candida albicans. Mol Microbiol. 2009;74(4):862–75. doi: 10.1111/j.1365-2958.2009.06898.x 19788542
55. deHart AK, Schnell JD, Allen DA, Tsai JY, Hicke L. Receptor internalization in yeast requires the Tor2-Rho1 signaling pathway. Mol Biol Cell. 2003;14(11):4676–84. doi: 10.1091/mbc.E03-05-0323 14593073
56. Schmidt A, Kunz J, Hall MN. TOR2 is required for organization of the actin cytoskeleton in yeast. Proc Natl Acad Sci U S A. 1996;93(24):13780–5. doi: 10.1073/pnas.93.24.13780 8943012
57. Ho H-L, Shiau Y-S, Chen M-Y. Saccharomyces cerevisiaeTSC11/AVO3 participates in regulating cell integrity and functionally interacts with components of the Tor2 complex. Curr Genet. 2005;47(5):273–88. doi: 10.1007/s00294-005-0570-8 15809876
58. Grahammer F, Helmstaedter M, Osenberg D, Kuhne L, Kretz O, Wanner N, et al. mTOR Regulates Endocytosis and Nutrient Transport in Proximal Tubular Cells. J Am Soc Nephrol. 2017;28(1):230–41. doi: 10.1681/ASN.2015111224 27297946
59. Hietakangas V, Cohen SM. TOR complex 2 is needed for cell cycle progression and anchorage-independent growth of MCF7 and PC3 tumor cells. BMC Cancer. 2008;8:No pp given.
60. Heese-Peck A, Pichler H, Zanolari B, Watanabe R, Daum G, Riezman H. Multiple functions of sterols in yeast endocytosis. Mol Biol Cell. 2002;13(8):2664–80. doi: 10.1091/mbc.E02-04-0186 12181337
61. Bahn YS, Molenda M, Staab JF, Lyman CA, Gordon LJ, Sundstrom P. Genome-wide transcriptional profiling of the cyclic AMP-dependent signaling pathway during morphogenic transitions of Candida albicans. Eukaryot Cell. 2007;6(12):2376–90. doi: 10.1128/EC.00318-07 17951520
62. Chaillot J, Cook MA, Corbeil J, Sellam A. Genome-Wide Screen for Haploinsufficient Cell Size Genes in the Opportunistic Yeast. G3 (Bethesda). 2017;7(2):355–60.
63. Thapa M, Bommakanti A, Shamsuzzaman M, Gregory B, Samsel L, Zengel JM, et al. Repressed synthesis of ribosomal proteins generates protein-specific cell cycle and morphological phenotypes. Mol Biol Cell. 2013;24(23):3620–33. doi: 10.1091/mbc.E13-02-0097 24109599
64. Franz R, Kelly SL, Lamb DC, Kelly DE, Ruhnke M, Morschhäuser J. Multiple molecular mechanisms contribute to a stepwise development of fluconazole resistance in clinical Candida albicans strains. Antimicrob Agents Chemother. 1998;42(12):3065–72. 9835492
65. Coste A, Ferrari S, Sanglard D. Antifungal drug resistance mechanisms in fungal pathogens from the perspective of transcriptional gene regulation. FEMS Yeast Research. 2009;9(7):1029–50. doi: 10.1111/j.1567-1364.2009.00578.x 19799636
66. Sasse C, Schillig R, Dierolf F, Weyler M, Schneider S, Mogavero S, et al. The transcription factor Ndt80 does not contribute to Mrr1-, Tac1-, and Upc2-mediated fluconazole resistance in Candida albicans. PloS one. 2011;6(9):e25623-e.
67. Schubert S, Barker KS, Znaidi S, Schneider S, Dierolf F, Dunkel N, et al. Regulation of efflux pump expression and drug resistance by the transcription factors Mrr1, Upc2, and Cap1 in Candida albicans. Antimicrobial agents and chemotherapy. 2011;55(5):2212–23. doi: 10.1128/AAC.01343-10 21402859
68. Pais P, Galocha M, Viana R, Cavalheiro M, Pereira D, Teixeira MC. Microevolution of the pathogenic yeasts Candida albicans and Candida glabrata during antifungal therapy and host infection. Microbial cell (Graz, Austria). 2019;6(3):142–59.
69. Ford CB, Funt JM, Abbey D, Issi L, Guiducci C, Martinez DA, et al. The evolution of drug resistance in clinical isolates of Candida albicans. eLife. 2015;4:e00662-e.
70. Boumendjel A, Beney C, Deka N, Mariotte AM, Lawson MA, Trompier D, et al. 4-Hydroxy-6-methoxyaurones with high-affinity binding to cytosolic domain of P-glycoprotein. Chem Pharm Bull (Tokyo). 2002;50(6):854–6.
71. Sim HM, Loh KY, Yeo WK, Lee CY, Go ML. Aurones as Modulators of ABCG2 and ABCB1: Synthesis and Structure–Activity Relationships. ChemMedChem. 2011;6(4):713–24. doi: 10.1002/cmdc.201000520 21302361
72. Pradines B. P-Glycoprotein-Like Transporters in Leishmania: A Search for Reversal Agents: Consequences, Molecular Mechanisms and Possible Treatments. 2018. p. 319–40.
73. Ryley JF, Wilson RG, Barrett-Bee KJ. Azole resistance in Candida albicans. Sabouraudia. 1984;22(1):53–63. 6322363
74. Basso V, Garcia A, Q. Tran D, Schaal J, Tran P, Ngole D, et al. Fungicidal Potency and Mechanisms of -Defensins against Multidrug-Resistant Candida Species2018. AAC.00111-18 p.
75. Jensen RH, Astvad KMT, Silva LV, Sanglard D, Jørgensen R, Nielsen KF, et al. Stepwise emergence of azole, echinocandin and amphotericin B multidrug resistance in vivo in Candida albicans orchestrated by multiple genetic alterations. The Journal of antimicrobial chemotherapy. 2015;70(9):2551–5. doi: 10.1093/jac/dkv140 26017038
76. Ferrari S, Ischer F, Calabrese D, Posteraro B, Sanguinetti M, Fadda G, et al. Gain of function mutations in CgPDR1 of Candida glabrata not only mediate antifungal resistance but also enhance virulence. PLoS pathogens. 2009;5(1):e1000268-e.
77. Pais P, Galocha M, Teixeira MC. Genome-Wide Response to Drugs and Stress in the Pathogenic Yeast Candida glabrata. In: Sá-Correia I, editor. Yeasts in Biotechnology and Human Health: Physiological Genomic Approaches. Cham: Springer International Publishing; 2019. p. 155–93.
78. Pfaller MA, Sheehan DJ, Rex JH. Determination of fungicidal activities against yeasts and molds: lessons learned from bactericidal testing and the need for standardization. Clin Microbiol Rev. 2004;17(2):268–80. doi: 10.1128/CMR.17.2.268-280.2004 15084501
79. Giaever G, Chu AM, Ni L, Connelly C, Riles L, Veronneau S, et al. Functional profiling of the Saccharomyces cerevisiae genome. Nature (London, U K). 2002;418(6896):387–91.
80. Parsons AB, Brost RL, Ding H, Li Z, Zhang C, Sheikh B, et al. Integration of chemical-genetic and genetic interaction data links bioactive compounds to cellular target pathways. Nat Biotechnol. 2004;22(1):62–9. doi: 10.1038/nbt919 14661025
81. Hoepfner D, Helliwell SB, Sadlish H, Schuierer S, Filipuzzi I, Brachat S, et al. High-resolution chemical dissection of a model eukaryote reveals targets, pathways and gene functions. Microbiol Res. 2014;169(2–3):107–20. doi: 10.1016/j.micres.2013.11.004 24360837
82. Lee AY, St. Onge RP, Proctor MJ, Wallace IM, Nile AH, Spagnuolo PA, et al. Mapping the Cellular Response to Small Molecules Using Chemogenomic Fitness Signatures. Science (Washington, DC, U S). 2014;344(6180):208–11.
83. Alsayari A, Muhsinah AB, Hassan MZ, Ahsan MJ, Alshehri JA, Begum N. Aurone: A biologically attractive scaffold as anticancer agent. Eur J Med Chem. 2019;166:417–31. doi: 10.1016/j.ejmech.2019.01.078 30739824
84. Pierce SE, Davis RW, Nislow C, Giaever G. Genome-wide analysis of barcoded Saccharomyces cerevisiae gene-deletion mutants in pooled cultures. Nat Protoc. 2007;2(11):2958–74. doi: 10.1038/nprot.2007.427 18007632
85. Xiao Y, Wang Y. Global discovery of protein kinases and other nucleotide-binding proteins by mass spectrometry. Mass Spectrom Rev. 2016;35(5):601–19. doi: 10.1002/mas.21447 25376990
86. Xiao Y, Guo L, Jiang X, Wang Y. Proteome-wide discovery and characterizations of nucleotide-binding proteins with affinity-labeled chemical probes. Anal Chem. 2013;85(6):3198–206. doi: 10.1021/ac303383c 23413923
87. Jakobi R, Traugh JA. Analysis of the ATP/GTP binding site of casein kinase II by site-directed mutagenesis. Physiol Chem Phys Med NMR. 1995;27(4):293–301. 8768785
88. Sicheri F, Moarefi I, Kuriyan J. Crystal structure of the Src family tyrosine kinase Hck. Nature. 1997;385(6617):602–9. doi: 10.1038/385602a0 9024658
89. Suryadinata R, Sadowski M, Sarcevic B. Control of cell cycle progression by phosphorylation of cyclin-dependent kinase (CDK) substrates. Bioscience Reports. 2010;30(4):243–55. doi: 10.1042/BSR20090171 20337599
90. Cote P, Hogues H, Whiteway M. Transcriptional analysis of the Candida albicans cell cycle. Mol Biol Cell. 2009;20(14):3363–73. doi: 10.1091/mbc.E09-03-0210 19477921
91. Song Y, Cheon SA, Lee KE, Lee SY, Lee BK, Oh DB, et al. Role of the RAM network in cell polarity and hyphal morphogenesis in Candida albicans. Mol Biol Cell. 2008;19(12):5456–77. doi: 10.1091/mbc.E08-03-0272 18843050
92. Klepser ME, Ernst EJ, Lewis RE, Ernst ME, Pfaller MA. Influence of test conditions on antifungal time-kill curve results: proposal for standardized methods. Antimicrob Agents Chemother. 1998;42(5):1207–12. 9593151
93. Piotrowski JS, Simpkins SW, Li SC, Deshpande R, McIlwain SJ, Ong IM, et al. Chemical Genomic Profiling via Barcode Sequencing to Predict Compound Mode of Action. Methods Mol Biol (N Y, NY, U S). 2015;1263(Chemical Biology):299–318.
94. Tisdall J. Beginning Perl for Bioinformatics. First Edition ed: O’REILLY; 2001.
95. Helliwell SB, Howald I, Barbet N, Hall MN. TOR2 is part of two related signaling pathways coordinating cell growth in Saccharomyces cerevisiae. Genetics. 1998;148(1):99–112. 9475724
96. Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008;3(6):1101–8. doi: 10.1038/nprot.2008.73 18546601
Článek vyšel v časopise
PLOS One
2019 Číslo 12
- Řekněte mi, co jíte, a já vám řeknu, jaká je úroveň vaší kognice
- Původní český test ALBA pro odhalení mírných kognitivních deficitů získal certifikaci MZ ČR
- Duševní onemocnění se u teenagerů přenáší v rámci školní třídy
- Od kostiček k lidské kůži. Jak si postavit 3D biotiskárnu?
- Prof. Libor Grubhoffer: Nová varianta mpox je příležitostí pokročit v poznání původce varioly
Nejčtenější v tomto čísle
- Methylsulfonylmethane increases osteogenesis and regulates the mineralization of the matrix by transglutaminase 2 in SHED cells
- Oregano powder reduces Streptococcus and increases SCFA concentration in a mixed bacterial culture assay
- The characteristic of patulous eustachian tube patients diagnosed by the JOS diagnostic criteria
- Parametric CAD modeling for open source scientific hardware: Comparing OpenSCAD and FreeCAD Python scripts