Effect of Conformational Variability on the Drug Resistance of Candida auris ERG11p and FKS1

Abstract

Candida auris infection has been recognized as an urgent threat to antifungal drug resistance, and the Eagle effect of C. auris FKS1 (1,3-β-d-glucan synthase) wild-type isolates has also been noted. The Eagle effect, namely, where higher concentrations of antifungals reduce fungicidal activity relative to lower concentrations, is a confounding factor of apparent antifungal resistance, but the detailed mechanism remains unclear. Here, we present the conformational variability of mutation sites for ERG11p (lanosterol 14α-demethylase) and FKS1 from deep neural network-based prediction along with the reported X-ray crystallographic and cryo-electron microscopy (cryo-EM) structures of antifungals. The sequence variability maps provide valuable insights into the inconsistent correlation between azole resistance and the mysterious Eagle effect with the dispersion of minimal inhibitory concentration (MIC) for echinocandin resistance. The conformational variability prediction supports the hypothesis that mutations K143R of clade I, VF125AL of clade III, and Y132F of clade IV for C. auris ERG11p make the corresponding site variable and that an increased population of invisible variable conformations potentially contributes to triazole resistance. In contrast, the predicted rigid conformation by the S639F mutation of hot spot region 1 (HS1) for FKS1 suggests that caspofungin (CAS) is involved in an uncompetitive inhibition, and a decreased population of the CAS-bound state of FKS1 with Rho1 leads to drug resistance. The predicted variable HS1 region for FKS1 WT isolates and the rigid one for FKS1 S639F mutants support the in vivo drug response and the in vitro MIC dispersion. A plausible mechanism of the Eagle effect is hereby proposed, namely, that a high concentration of CAS with a high membrane affinity reduces the population of the CAS-bound state of FKS1 with Rho1, as well as accompanying events such as aggregation or association depending on the conformational variability of HS1.