ABSTRACT
Meyerozyma guilliermondii is a rare opportunistic fungal pathogen that causes deadly invasive candidiasis in human. M. guilliermondii strain SO is a local yeast isolate that possesses huge industrial interests but also pathogenic towards zebrafish embryos. Enolases that bind to human extracellular matrix (ECM) proteins are among the fungal virulence factors. To understand its pathogenicity mechanism down to molecular level, especially in the rare M. guilliermondii, this study aimed to identify and characterize the potentially virulence-associated enolase in M. guilliermondii strain SO using bioinformatics approaches. Profile Hidden-Markov model was implemented to identify enolase-related sequences in the fungal proteome. Sequence analysis deciphered only one (MgEno4581) out of nine sequences exhibited potent virulence traits observed similarly in the pathogenic Candida albicans. MgEno4581 structure that was predicted via SWISS-MODEL using C. albicans enolase (CaEno1; PDB ID: 7vrd) as the homology modeling template portrayed a highly identical motif with CaEno1 that facilitates ECM proteins binding. Amino acid substitutions (D234K, K235A, Y238H, K239D, G243K, V248C and Y254F) in ECM-binding motif of Saccharomyces cerevisiae enolase (ScEno) compared to MgEno4581 and CaEno1 caused changes in motif's surface charges. Protein-protein docking indicated F253 in ScEno only interacted hydrophobically with human plasminogen (HPG). Hydrogen linkages were observed for both MgEno4581 and CaEno1, suggesting a stronger interaction with HPG in the hydrophilic host microenvironments. Thus, our in silico characterizations on MgEno4581 provided new perspectives on its potential roles in candidiasis (fungal-host interactions) caused by M. guilliermondii, especially M. guilliermondii strain SO on zebrafish embryos that mimic the immunocompromised individuals as previously evident.
