Abstract
This study was designed to understand the molecular mechanisms of azole resistance in Candida tropicalis using genetic and bioinformatics approaches. Thirty-two azole-resistant and 10 azole-susceptible (S) clinical isolates of C. tropicalis were subjected to mutation analysis of the azole target genes including ERG11. Inducible expression analysis of 17 other genes potentially associated with azole resistance was also evaluated. Homology modeling and molecular docking analysis were performed to study the effect of amino acid alterations in mediating azole resistance. Of the 32 resistant isolates, 12 (37.5%) showed A395T and C461T mutations in the ERG11 gene. The mean overexpression of CDR1, CDR3, TAC1, ERG1, ERG2, ERG3, ERG11, UPC2, and MKC1 in resistant isolates without mutation (R-WTM) was significantly higher (p<0.05) than those with mutation (R-WM) and the sensitive isolates (3.2-11 vs. 0.2-2.5 and 0.3-2.2 folds, respectively). Although the R-WTM and R-WM had higher (p<0.05) CDR2 and MRR1 expression compared to S isolates, noticeable variation was not seen among the other genes. Protein homology modelling and molecular docking revealed that the mutations in the ERG11 gene were responsible for structural alteration and low binding efficiency between ERG11p and ligands. Isolates with ERG11 mutations also presented A220C in ERG1 and together T503C, G751A mutations in UPC2. Nonsynonymous mutations in the ERG11 gene and coordinated overexpression of various genes including different transporters, ergosterol biosynthesis pathway, transcription factors, and stress-responsive genes are associated with azole resistance in clinical isolates of C. tropicalis.