key: cord-0838962-3aerysms
authors: Andrade-Pavón, Dulce; Fernández-Muñoz, Vanessa; González-Ibarra, Wendy; Hernández-Rodríguez, César; Ibarra, J. Antonio; Villa-Tanaca, Lourdes
title: Point mutations in Candida glabrata 3-hydroxy-3-methylglutaryl-coenzyme A reductase (CgHMGR) decrease enzymatic activity and substrate/inhibitor affinity
date: 2021-10-21
journal: Sci Rep
DOI: 10.1038/s41598-021-00356-w
sha: 30f1c7850a14b8a4f3eeff6e742d078eb90e23ca
doc_id: 838962
cord_uid: 3aerysms

3-Hydroxy-3-methylglutaryl-coenzyme A reductase (HMGR) is a crucial enzyme in the ergosterol biosynthesis pathway. The aim of this study was to obtain, purify, characterize, and overexpress five point mutations in highly conserved regions of the catalytic domain of Candida glabrata HMGR (CgHMGR) to explore the function of key amino acid residues in enzymatic activity. Glutamic acid (Glu) was substituted by glutamine in the E680Q mutant (at the dimerization site), Glu by glutamine in E711Q (at the substrate binding site), aspartic acid by alanine in D805A, and methionine by arginine in M807R (the latter two at the cofactor binding site). A double mutation, E680Q-M807R, was included. Regarding recombinant and wild-type CgHMGR, in vitro enzymatic activity was significantly lower for the former, as was the in silico binding energy of simvastatin, alpha-asarone and the HMG-CoA substrate. E711Q displayed the lowest enzymatic activity and binding energy, suggesting the importance of Glu(711) (in the substrate binding site). The double mutant CgHMGR E680Q-M807R exhibited the second lowest enzymatic activity. Based on the values of the kinetic parameters K(M) and V(max), the mutated amino acids appear to participate in binding. The current findings provide insights into the role of residues in the catalytic site of CgHMGR.

The ergosterol biosynthesis pathway has been proposed as a target for identifying and designing new antifungals. In a rate-limiting step of the biosynthetic pathway of ergosterol in fungi and cholesterol in humans, the enzyme 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMGR) catalyzes the conversion of 3-hydroxy-3-methylglutaryl-coenzyme A (HMG-CoA) to mevalonate 1 . This pathway also generates other bioactive molecules, including a coenzyme (Q10) that participates in the transport of electrons in mitochondria 2 .

Phylogenetic analyses have revealed two types of HMGR: class I of eukaryotic organisms and some archaea and class II of prokaryotes and some archaea. Whereas the first class has three domains (catalytic, linker, and membrane-anchor), most members of the second one contain only the catalytic domain 1 .

The HMGR enzyme is the target of lipid-lowering drugs such as statins 3 . These drugs have been described as competitive inhibitors because they dispute the active site of HMGR with the substrate HMG-CoA 1 . Since statins can induce myopathy and hepatotoxicity 4 , new bioactive principles derived from natural sources (e.g., alpha-asarone) have been synthesized. In a murine model, alpha-asarone decreased cholesterol and triglyceride (S)-HMG-CoA + 2NADPH + 2H + → (R)-mevalonate + 2NADP + + CoASH

Conservation of the sites of interest in the catalytic domain of HMGRs from different organisms. A multiple alignment was carried out on the Clustal Omega program for the amino acid sequence of the soluble fraction of HMGR from plants, bacteria, insects, mammals (including Homo sapiens), and opportunistic pathogenic yeasts (e.g., C. glabrata) ( Supplementary Fig. S1 ). The access codes of the corresponding sequences can be consulted in the Materials and Methods section. The motif sequences in the sites of interest in the HMGR enzymes are highlighted in Fig. 1 : the dimerization site (ENVIG), the substrate binding site (EGCLVAS), and the cofactor binding site (DAMGMN). The amino acid sequences encoding the soluble fraction of the HMGR enzymes of different organisms are highly conserved, which emphasizes the importance of examining some of the amino acid residues in these sites in order to evaluate their role in the catalysis of the enzyme. Such an analysis has been performed in studies on P. mevalonii.

A phylogenetic tree was constructed by examining the amino acid sequences detected in the catalytic domain of the twenty-four organisms included in the alignment test (Fig. 2) . The phylogenetic analysis of CgHMGR demonstrated that it belongs to the clade of the yeast species of the Saccharomycotina subphylum and the Saccharomyces genus, not to the clade containing C. albicans. Molecular modeling of wild-type and mutant CgHMGR proteins. The homology modeling of the five mutants herein generated was carried out in the Modeller program, using the crystal structure of hHMGR (PDB: 1DQA) as the template. Of the fifteen models of each mutant provided, the one with the lowest DOPE score (the most thermodynamically stable) was selected. The selected models are portrayed in Fig. 3 , highlighting in red the amino acids that were changed in the transition from the wild-type (Fig. 3a) to the mutant proteins ( Fig. 3b-f ): for E680Q, Glu66 to Gln66; for E711Q, Glu97 to Gln97; for D805A, Asp191 to Ala191; and for M807R, Met193 to Arg193.

The models were validated with the PROCHECK program, which quantifies the loss of residues in the allowed regions of a Ramachandran plot ( Supplementary Fig. S2 ). In all the models, over 90% of the residues are found in the most favorable regions, with no residues in the disallowed regions (Supplementary Table 1 ). All the models constructed had Ramachandran statistics similar to those reported in the wild-type CgHMGR model. These results confirmed the high quality and reliability of the models generated.

Docking study of ligand binding to wild-type and mutant HMGR proteins. An in-silico analysis was carried out for the interaction of mutant and wild-type CgHMGR proteins with two of their inhibitors (simvastatin and alpha-asarone) and their substrate (HMG-CoA) (Fig. S5) . Each mutant showed a lower binding energy for each ligand than the wild-type CgHMGR (Fig. S5 , Table 2 ). The mutant protein E711Q displayed the lowest value (Fig. S5) . Glu97, Asp307, Lys321 and His403 are the residues most frequently involved in binding, being located in the catalytic domain of CgHMGR. However, not all of these amino acids participate in the binding of the ligands with each of the mutant peptides (Table 2) .

Regarding the binding of simvastatin, alpha-asarone, and HMG-CoA to the active site of the CgHMGR enzyme, the docking simulations gave insights into the binding energies as well as the residues, polar interactions, and hydrophobic interactions involved (Table 1) .

Generation and verification of the point mutants of HMGR genes. PCR reactions were carried out with specific oligonucleotides designed to obtain the following mutants: CgHMGRE680Q, CgHMGRE711Q, CgHMGRD805A, CgHMGRM807R and CgHMGRE680Q-M807R. The plasmid Rec-MBP-CgHMGR was used as the template DNA (where MBP refers to the maltose-binding protein). The amplified oligonucleotides had a size of 7.9 kbp ( Supplementary Fig. S3a ), which corresponds to the size of the plasmid pMAL-C2X (6.6 kbp) plus the gene encoding the soluble fraction of CgHMGR (1.3 kbp). The full-length gel original is included in Supplementary Information (Fig. S3c) . Additionally, the plasmids afforded by the transformation of each of the amplified products were subjected to double digestion ( Supplementary Fig. S3b ), finding the expected molecules of 6.6 and 1.3 kbp. The full-length blot is also in Supplementary Information (Fig. S3d ). www.nature.com/scientificreports/ Likewise, the mutations generated by sequencing were corroborated. The sequences obtained were aligned with the coding sequence of wild-type CgHMGR to confirm that the expected mutations had indeed been introduced. For CgHMGRE680Q and CgHMGRE711Q, the triplet changed from GAA to CAA, though in a distinct position in each case. The GAT triplet changed to GCT in CgHMGRD805A, and CAT to CCT in CgH-MGRM807R. For the double mutant CgHMGRE680Q-M807R, TTC changed to TTG and CAT to CCT (Fig.  S6 ). The sequencing revealed that the point mutations were introduced successfully.

Expression, detection, and purification of CgHMGR mutant proteins. The protein expression of the CgHMGRE680Q, CgHMGRE711Q, CgHMGRD805A, CgHMGRM807R and CgHMGRE680Q-M807R mutants was induced with IPTG. Samples taken before and after induction were visualized with SDS-PAGE. The gel in Fig. 4b shows the overexpression of the proteins of interest. The approximate molecular mass of 90 kDa corresponds to the fusion of CgHMGR + MBP (Fig. 4a) . The pMAL vector favored the expression of the recombinant proteins, which fused to the soluble MBP. The latter facilitated the isolation of the enzyme after its overexpression.

Hence, the mutants were overexpressed and the presence of MBP fused to the mutants was verified. The mutants found were the same as those previously identified by using anti-MBP antibodies (Fig. 4c) . Subsequently, the mutants were purified easily (in a few steps) by means of amylose affinity chromatography. In all cases, a single band appeared at the approximate molecular mass of 90 kDa, exhibiting the same behavior as wild-type CgHMGR.

It was demonstrated by a coupling study that the substrate (HMG-CoA) and probable inhibitors (simvastatin and alpha-asarone) of CgHMGR did not recognize a specific site in MBP. Based on the resulting low binding energy values, these compounds seem to have very limited affinity towards MBP. Moreover, HMG-CoA Reduced enzymatic activity of the mutants compared to wild-type CgHMGR. The enzymatic activity of the CgHMGRE680Q, CgHMGRE711Q, CgHMGRD805A, CgHMGRM807R and CgHMGRE680Q-M807R proteins was measured spectrophotometrically by reading the decrease in absorbance of NADPH at 340 nm during 10 min. For each mutant, this parameter was determined in triplicate and expressed as the percentage of activity of the wild-type protein (considered as 100%) ( Table 2 ). Compared to wild-type CgHMGR, all five mutants displayed a lower enzymatic activity. The smallest percentage of residual activity (23%) was observed for E711Q. The least affected enzyme was E680Q (mutated at the dimerization site), with 88.2% residual activity. The double mutant E680Q-M807R exhibited the second lowest activity (44%), while the residual activity of the corresponding single mutants, E680Q and M807R, was substantially higher (82% and 79%, respectively). The enzymatic activity of each mutant was significantly less than that of the wild-type protein, according to a paired Student's t-type analysis.

Kinetic parameters of the mutants. For the determination of kinetic parameters (K M and V max) , activity was assayed in a reaction mixture containing increasing amounts of HMG-CoA, ranging from 1 to 64 μM. Enzymatic activity was obtained from monitoring the conversion of HMGCoA-dependent oxidation of NADPH and taking the linear fraction for each of the substrate concentrations to obtain a Michaelis-Menten plot (Fig.  S7 ). The data were evaluated with the Lineweaver-Burk and Hanes-Woolf models (Fig. S7) . The Hanes-Woolf model was used because of explaining the results better than the other. Moreover, the Hanes-Woolf data analysis led to an outcome comparable to that described in a previous report on the wild-type protein 10 . The present data for the K M and V max was as follows: 8.31 µM and 0.28 µM min −1 for CgHMGRE680Q, 18 µM and 0.70 µM min −1 for CgHMGRE711Q, and 11.84 µM and 0.43 µM min −1 for CgHMGRD805A. 

HMGR is a highly conserved enzyme anchored to the endoplasmic reticulum membrane of eukaryotes and to the plasma membrane of prokaryotes. It is responsible for the synthesis of cholesterol in mammals, ergosterol in fungi, and isoprenoids in bacteria 1 . Intense research efforts have focused on hHMGR as the target of inhibitors such as statins and fibrates, which diminish the synthesis of cholesterol and reduce the risk of cardiac arrest in patients with hyperlipidemia. Consequently, hHMGR has been examined at many levels: transcription, translational regulation, synthesis, inhibition, crystallization, and tertiary structure. However, reports on point mutations in non-human HMGR are limited, to the best of our knowledge, to P. mevalonii and Syrian hamsters. The current investigation is the first to make point mutations in the HMGR enzyme of a yeast of medical interest.

The alignment presently carried out with the amino acid sequences of HMGRs from various species demonstrated that the motif sequences of the dimerization site (ENVIG), the substrate binding site (EGCLVAS), and the cofactor binding site (DAMGMN) are highly conserved, regardless of whether they are from fungal, mammalian, plant or bacterial proteins. Hence, these amino acid residues appear to have an essential function in the catalytic reaction of the enzyme and/or in its conformation and tertiary structure, which makes them candidates for the evaluation of point mutants in key sequences.

On the other hand, the phylogenetic analysis of the HMGR proteins from different yeasts revealed a classification in accordance with their taxonomy. The majority are of the Phylum Ascomycota, which includes the Candida and Saccharomyces genera. CgHMGR is grouped with the yeast proteins of the WGD clade, encompassing K. lactis and S. cerevisiae 17, 18 , and not with those of the CTG clade, containing C. albicans and other Candida species. Therefore, C. glabrata is a species phylogenetically closer to S. cerevisiae than to C. albicans. The HMGRs of Y. lipolytica, S. pombe and U. maydis (the latter being a fungus of the phylum Basidiomycota) belong to branches completely independent of Ascomycota. The HMGRs of other eukaryotes (e.g., plants, insects and mammals) www.nature.com/scientificreports/ are grouped into an independent clade and the HMGR of the bacterium P. mevalonii functions as an external group, further validating the previous results 1 . hHMGR, purified and crystallized in the year 2000, was used as a template for the molecular modeling of proteins 19 . Although the three-dimensional structure of the catalytic domain of hHMGR (426-888 aa) is a tetramer, it has been suggested that the protein dimer might be able to bind to the HMG-CoA substrate 19 . The molecular modeling of the five CgHMGR with point mutations generated proteins made up of two identical subunits (dimers). The structures of the mutated proteins provided an a priori approximation to the 3D structure of the corresponding peptides, a necessary step for their in-silico analysis 11 .

According to kinetic studies on hHMGR, statins compete for the HMG-CoA substrate without affecting binding to the NADPH cofactor 20 . Hence, the substitution of an amino acid residue in the conserved sequences could plausibly alter the affinity of the enzyme for its substrate HMG-CoA as well as for simvastatin and alphaasarone. To explore such a possibility, a molecular docking study was carried out. This is the first report, to our knowledge, of an evaluation of molecular coupling between mutant proteins of CgHMGR and reference inhibitors. Information is herein provided on the substantial effect of mutations to specific amino acids on the ligand-receptor interaction.

The amino acids Glu97, Asp307, Lys321 and His403 are part of the catalytic domain of CgHMGR and their plausible role in the ligand-HMGR interaction has been proposed 1,11 . According to previous docking studies, the amino acids affected by the mutations likely influence the recognition of the substrate/inhibitor and/or alter the native structure of the protein 21 .

Statins and alpha-asarone are competitive inhibitors of HMGRs, blocking access to the HMG-CoA substrate 1,3,10,11 . The CgHMGR mutant in the substrate binding motif, E711Q, displayed the lowest binding energy for the ligands (simvastatin, alpha-asarone, and the HMG-CoA substrate) compared to the other mutants herein examined. The limited enzymatic activity and binding energy of this mutant was expected, considering that the three ligands currently tested interact with the EGCLVAS motif (of the substrate binding site) 3,10,11 . Glutamic acid in the EGCLVAS region seems to have an important function in the catalysis of CgHMGR, as previously reported 12 . www.nature.com/scientificreports/ CgHMGR-M807R was mutated at the cofactor binding site, specifically at the methionine position 807 of the DAMGMN sequence 12 . The change was from methionine, an amino acid with a neutral and nonpolar charge, to arginine, one with a positive charge. Although both amino acids participate in protein methylation 22 , the positive charge of arginine could be detrimental to the binding of the cofactor to the enzyme. Additionally, arginine contains a guanidinium group, which when ionized has a lower charge density than other amino acids. While methionine is a weak nucleophile and cannot be protonated 22 , the positive charge of arginine might be able to alter the α-helix or β-folded structure of the protein and consequently its enzymatic activity. Finally, the larger Supplementary Information (Fig. S4a) . (c) Example of the verification of protein identity by Western blot based on anti-MBP antibodies: Lane 1, molecular weight marker; lane 2, uninduced cell extract; lane 3, induced cell extract. The full-length blot is also shown in Supplementary Information (Fig. S4b ). www.nature.com/scientificreports/ size of arginine than methionine may modify the structure of the protein. It is possible that these factors hindered the proper binding of the enzyme to its cofactor, which would contribute to the reduced enzymatic activity. The mutation made in the dimerization site caused glutamic acid to be replaced by glutamine in position 680. Although both amino acids are very similar in size, glutamic acid has a negative charge and glutamine a polar neutral charge. Glutamic acid has carboxylate (COO-) side chains that are potential proton acceptors, forming hydrogen bonds and thus the secondary structure of the protein (β-folded sheets or α-helix) 23 . This amino acid could possibly have a similar function as the Glu83 residue, which is highly conserved and participates in catalysis by transferring protons to Lys267. When glutamic acid carries out the decomposition of mevaldyl-CoA, its protonation and deprotonation confers distinct levels of energy to the reaction, and these are necessary for the enzymatic activity that gives rise to sterols 24 . Glutamine is a glutamic acid amide afforded by the replacement of the hydroxyl of glutamic acid with an amine group 24 . Since it has a polar neutral charge, however, it might have affected the interaction with the amino acids belonging to the ENVIG sequence. Therefore, the proper formation of the dimers may not have taken place, which would lead to the inadequate binding of the enzyme to the cofactor and consequently to limited enzymatic activity.

It is not surprising that the double mutant (CgHMGRE680Q-M807R) exhibited the second lowest enzymatic activity of the present mutants. The two mutations occurred at sites where amino acids are highly conserved: the sequence of the cofactor binding site (DAMGMN) and of the dimerization site (ENVIG).

There was a greater K M value found for each mutant than for the wild-type protein 10 . Hence, the data suggest that the affinity of HMG-CoA would be lower for the mutant versions than for the wild-type protein. Indeed, a marked decrease was detected in the enzymatic activity of the mutants herein examined compared to that of the wild-type protein.

C. glabrata has only one HMG1 gene, which is orthologous to the S. cerevisiae gene 7 . In the current contribution, the sensitivity of CgHMGR mutants to antifungals was not evaluated. According to a previous study, the mutation of this gene in Aspergillus fumigatus causes resistance to antifungal compounds of the triazole class 25, 26 . The mutation of the enzymes HMG1 and HMG2 of S. cerevisiae leads to a great sensitivity to compactin, an HMGR inhibitor 27 . There are no reports, to our knowledge, on a double mutant cdr1∆ hmg1∆. However, these mutants have been assessed separately, observing an increase in sensitivity to azoles caused by cdr1∆ in C. glabrata as well as resistance to the same drugs produced by hmg1∆ in Aspergillus fumigatus 25, 26, 28 .

The insights gained from the presently mutated proteins can be applied not only to the improvement of C. glabrata as a model for research on drug resistance, but also to the investigation of the HMGR of additional species, such as C. auris and C. haemulonii, which have recently emerged as pathogenic multi-drug resistant strains 29 . The emergence of Candida infections associated with SARS-Cov2 gives greater emphasis to the urgency of developing antifungal agents for alternative targets 30 . Table 2 . The enzymatic activity and binding energies of the point mutations made in the motifs of the dimerization site, substrate binding site, and cofactor binding site of CgHMGR. Point mutations were made in the catalytic domain of CgHMGR to express various recombinant proteins (Rec-CgHMGR) (as described in Material and Methods). The aa of the wild-type sequence is shown in green and that of the mutated protein in red. The site of the motif sequence is indicated. The binding energies of CgHMGR (wild-type and mutated) with its HMG-CoA substrate and probable inhibitors (simvastatin and alpha-asarone) were estimated by docking. The enzymatic activity of the recombinant enzymes is calculated as a percentage of the activity of the wild-type enzyme (considered as 100%). Data are expressed as the average of three independent replicates ± standard deviation (SD). ***P < 0.0005 **P < 0.0052, *P < 0.0153. Significant differences compared to the wild-type protein based on a paired Student's t-type analysis. 

Five point mutants were made at three highly conserved sites of the catalytic domain of HMGR: the dimerization, substrate binding and cofactor binding sites. After the mutants were expressed, they were detected and purified to verify their integrity. All the mutants exhibited significantly lower enzymatic activity and ligand affinity than wild-type CgHMGR. The in vitro and in silico testing of the mutants gave insights into the function of specific amino acid residues in the catalytic site of the enzyme, especially glutamic acid. The current findings should certainly facilitate the creation of new mutants of HMGR to further explore the catalysis, recognition, conformation, and dimerization mechanisms involved in the activity of the enzyme and the viability of C. glabrata and other opportunistic pathogenic fungi.

Multiple sequence alignment of HMGR enzymes from different organisms. The amino acid sequences were downloaded from the NCBI database (https:// www. ncbi. nlm. nih. gov/) for the HMGR enzymes of the following organisms: 

A phylogenetic tree was constructed with the sequences of the catalytic domain of the aforementioned twenty-four HMGRs in the MEGA6 program 34 , utilizing the maximum likelihood method and the WAG + G model. Additionally, a bootstrap replication of 100 was used to evaluate the reliability of the phylogenetic tree.

CgHMGR mutant proteins with Modeller 9.24 35 . The models with the lowest discrete optimized protein energy (DOPE) score were selected and edited in Discovery Studio 36 , then subjected to a validation process to assess their stereochemical quality with the PROCHECK 37 Microorganisms and plasmids. Transformation was achieved with the Escherichia coli DH10β strain and the overexpression of heterologous proteins with the E. coli BL21 strain. The plasmid pMAL-C2X-CgHMGR, which was previously obtained by our group and harbors the gene encoding the wild-type protein, served as the template DNA for generating the mutants 10 .

Culture media. E. coli harboring the plasmids was grown in liquid Luria Bertani (LB) medium (1.0% casein peptone, 0.5% yeast extract, and 0.5% NaCl) and solid LB medium with 1.5% bacteriological agar. These media were adjusted to pH 7.0 and supplemented with 100 μg/mL ampicillin. Culture media and antibiotics were acquired from Sigma Aldrich (St. Louis, MO, USA).

Mutagenic oligonucleotide design. The mutagenic primers were designed manually, taking as templates the sequences of the dimerization site, the substrate binding site, and the cofactor binding site within the catalytic domain of the HMGR enzyme. Consideration was given to 18 base pairs to the left and 18 base pairs to the right of the codon where the mutation was designed, resulting in a total of 39 nucleotides for each primer. Point mutations were made involving the following changes: glutamic acid to glutamine in the dimerization site (ENVIG) to provide E680Q, aspartic acid to glutamine in the substrate binding site (EGCLVAS) to afford E711Q, aspartic acid to alanine in the cofactor binding site (DAMGMN) to furnish (D805A), and methionine to arginine in the cofactor binding site to yield M807R. The double mutation of E680Q-M807R was formed by the respective single point mutations. Table 3 shows the original sequence (wild-type CgHMGR) and the mutagenic oligonucleotides. The codon of the amino acid to be mutated is illustrated in red, and the change made is underlined.

Generation of CgHMGR point mutants by PCR. The E680Q, E711Q and D805A point mutants were generated by PCR with a "Quick Change" commercial kit (Thermo Fisher Scientific, USA), as specified in the manufacturer's instructions with some modifications. This procedure is based on a PCR reaction with two complementary and antiparallel oligonucleotides and the appropriate DNA polymerase to achieve the desired mutation. Each oligonucleotide hybridizes with one of the two strands of the plasmid to be amplified. The plasmid obtained led to the sequence to be mutated (pMAL-c2X-CgHMGR), which was used as the template DNA. Mutants M807R and E680Q-M807R were produced by inverse PCR, with plasmid pMBP-CgHMGR as the template DNA 10, 11 . For the double mutant, the plasmid DNA of the rec-CgHMGR-E680Q mutant and the mutagenic oligonucleotide M807R served as the template DNA.

Once the reaction conditions were established and developed, the amplicon was digested with the Dpn-I enzyme (Thermo Fisher Scientific, USA). E. coli DH10β cells were then transformed with the digested PCR product. Subsequently, E. coli cells were plated on LB agar with ampicillin (100 µg/ml) and incubated at 37 °C for 24 h. The transforming colonies were selected for the extraction of plasmid DNA. Plasmids were purified with the ZymoPURE plasmid Miniprep Kit (Zymo Research, CA, USA) and confirmed by sequencing them according to the Sanger method, a procedure that demonstrated the presence of the desired mutations (Macrogen®, Seoul Korea).

Once the presence of the desired mutations was confirmed by sequencing the corresponding plasmids, the recombinant proteins were expressed. Briefly, E. coli BL21 cells were transformed with the mutated DNA by growing the cells containing the fusion plasmid in LB medium plus ampicillin (100 µg/mL) and 0.2% glucose until reaching an optical density of 0.6 at 600 nm. IPTG was added immediately (0.3 mM final concentration) to induce the expression of the protein, a process carried out at 37 °C for 4 h. Upon completion of this time, cells were harvested at 4000 × g for 20 min at 4 °C. The cells were mechanically broken by using glass beads and a vortex, and then the cell pack was dissolved in a column buffer (1 M Tris-HCl, pH 7.4, 11.7 g NaCl, 0.5 M EDTA, and 154 mg DTT) to obtain a cell-free extract. The cells were disrupted for 20 min, alternating every other 30 s on and off ice. The cell lysate was centrifuged at 18,000 g for 20 min and the supernatant (cell-free extract) was recovered. The protein was purified by affinity column chromatography with amylose resin (New England Biolabs, USA). The wild-type and mutant proteins were eluted with a buffer supplemented with 10 mM maltose. Fractions were collected to determine CgHMGR activity and protein concentration by the Lowry method. The expression and the successful purification of the mutant proteins was visualized by means of 10% SDS-PAGE stained with Coomassie blue. The mutated proteins were detected by Western blot with the anti-MBP monoclonal antibody (New England Biolabs, USA), according to the manufacturer's specifications.

The HMGR enzymatic activity assay. The enzymatic activity of the five mutants was quantified by the method of Bischoff and Rodwell 41 , beginning with the elaboration of a reaction mixture containing 10 mM HMG-CoA, 50 mM Tris-HCl regulator, injectable water, 9 µg of the protein of the respective mutant, and 10 mM NADPH. The cofactor (NADPH) is oxidized by the catalytic subunit of HMGR in the presence of the substrate HMG-CoA, causing a decrease in absorbance. Thus, absorbance was monitored at 340 nm for 10 min. Reagents were purchased from Sigma Aldrich (St Louis, MO, USA). One unit of specific activity was defined Table 3 . Oligonucleotides designed to obtain the following point mutants: E680Q, E711Q, D805A, M807R and E680Q-M807R. T m : melting temperature; %GC: Guanine and Cytosine Percentage. Substrate specificity and kinetic parameters. The steady-state kinetic constants for mutants were established by means of the standard activity assay under the same conditions described in the previous section. Activity was determined by incubating reaction mixtures with substrate concentrations ranging from 1 to 64 µM. The K M and V max constants were obtained from Hanes-Woolf plots and expressed as the mean of three different experiments 10, 42 .

Received: 12 June 2021; Accepted: 11 October 2021

The 3-hydroxy-3-methylglutaryl coenzyme-A (HMG-CoA) reductases

HMG-CoA reductase inhibitors and coenzyme Q10

Structural mechanism for statin inhibition of HMG-CoA reductase

Liver and statins: A critical appraisal of the evidence

Design, synthesis, and docking of highly hypolipidemic agents: Schizosaccharomyces pombe as a new model for evaluating alpha-asarone-based HMG-CoA reductase inhibitors

Synthesis and highly potent hypolipidemic activity of alpha-asarone-and fibrate-based 2-acyl and 2-alkyl phenols as HMG-CoA reductase inhibitors

The 3-hydroxy-3-methylglutaryl coenzyme-A reductases from fungi: A proposal as a therapeutic target and as a study model

Resistance of Candida to azoles and echinocandins worldwide

Statin resistance in Candida glabrata

Recombinant 3-Hydroxy 3-Methyl Glutaryl-CoA Reductase from Candida glabrata (Rec-CgHMGR) Obtained by heterologous expression, as a novel therapeutic target model for testing synthetic drugs

Inhibition of recombinant enzyme 3-hydroxy-3-methylglutaryl-CoA reductase from Candida glabrata by α-asarone-based synthetic compounds as antifungal agents

Identification of the principal catalytically important acidic residue of 3-hydroxy-3-methylglutaryl coenzyme A reductase

Identification of the catalytically important histidine of 3-hydroxy-3-methylglutarylcoenzyme A reductase

Role of cysteine residues in Pseudomonas mevalonii 3-hydroxy-3-methylglutaryl-CoA reductase: Site-directed mutagenesis and characterization of the mutant enzymes

His865 is the catalytically important histidyl residue of Syrian hamster 3-hydroxy-3-methylglutarylcoenzyme A reductase

Catalysis by Syrian hamster 3-hydroxy-3-methylglutaryl-coenzyme A reductase: Proposed roles of histidine 865, glutamate 558, and aspartate 766

Comparative genomics of emerging pathogens in the Candida glabrata clade

Genome structure and dynamics of the yeast pathogen Candida glabrata

Crystal structure of the catalytic portion of human HMG-CoA reductase: Insights into regulation of activity and catalysis

Competitive inhibition of 3-hydroxy-3-methylglutaryl-coenzyme A reductase by ML-236A and ML-236B fungal metabolites, having hipoclelesterolemic activity

Amino acids: Metabolism, functions, and nutrition

The role of methionine on metabolism, oxidative stress, and diseases

A protein functionalization platform based on selective reactions at methionine residues

Substrate-induced closure of the flap domain in the ternary complex structures provides insights into the mechanism of catalysis by 3-hydroxy-3-methylglutaril-CoA reductase

Mutations in hmg1, challenging the paradigm of clinical triazole resistance in Aspergillus fumigatus

Non-cyp51A Azole-Resistant Aspergillus fumigatus isolates with mutation in HMG-CoA Reductase

Saccharomyces cerevisiae contains two functional genes encoding 3-hydroxy-3-methylglutaryl-coenzyme A reductase

Unexpected effects of azole transporter inhibitors on antifungal susceptibility in Candida glabrata and other pathogenic Candida species

First report of a catheterrelated bloodstream infection by Candida haemulonii in a children´s hospital in Mexico City

Increased antimicrobial resistance during COVID-19 pandemic

Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega

ExPASy: The proteomics server for in-depth protein knowledge and analysis

Molecular evolutionary genetics analysis version 6.0

Comparative Protein structure modeling using MODELLER

PROCHECK: A program to check the stereochemical quality of protein structures

PubChem BioAssay: 2017 update

Open Babel: An open chemical toolbox

AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility

Biosynthesis and characterization of (S)-and (R)-3-hydroxy-3-methylglutaryl coenzyme A

3-hydroxy-3-methylglutaryl-CoA reductase from rat liver

The authors would like to thank Bruce Allan Larsen for proofreading the manuscript. 

The authors declare no competing interests.

Supplementary Information The online version contains supplementary material available at https:// doi. org/ 10. 1038/ s41598-021-00356-w.Correspondence and requests for materials should be addressed to L.V.-T.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.