Interactome of Glyceraldehyde-3-Phosphate Dehydrogenase Points to the Existence of Metabolons in Paracoccidioides lutzii

Paracoccidioides is a dimorphic fungus, the causative agent of paracoccidioidomycosis. The disease is endemic within Latin America and prevalent in Brazil. The treatment is based on azoles, sulfonamides and amphotericin B. The seeking for new treatment approaches is a real necessity for neglected infections. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) is an essential glycolytic enzyme, well known for its multitude of functions within cells, therefore categorized as a moonlight protein. To our knowledge, this is the first approach performed on the Paracoccidioides genus regarding the description of PPIs having GAPDH as a target. Here, we show an overview of experimental GAPDH interactome in different phases of Paracoccidioides lutzii and an in silico analysis of 18 proteins partners. GAPDH interacted with 207 proteins in P. lutzii. Several proteins bound to GAPDH in mycelium, transition and yeast phases are common to important pathways such as glycolysis and TCA. We performed a co-immunoprecipitation assay to validate the complex formed by GAPDH with triose phosphate isomerase, enolase, isocitrate lyase and 2-methylcitrate synthase. We found GAPDH participating in complexes with proteins of specific pathways, indicating the existence of a glycolytic and a TCA metabolon in P. lutzii. GAPDH interacted with several proteins that undergoes regulation by nitrosylation. In addition, we modeled the GAPDH 3-D structure, performed molecular dynamics and molecular docking in order to identify the interacting interface between GAPDH and the interacting proteins. Despite the large number of interacting proteins, GAPDH has only four main regions of contact with interacting proteins, reflecting its ancestrality and conservation over evolution.

Introduction

Paracoccidioides is a dimorphic fungi and the genus includes species that cause paracoccidioidomycosis (PCM). The disease is endemic in Latin America and prevalent in Brazil (Brummer et al., 1993). The genus Paracoccidioides comprises five distinct species Paracoccidioides brasiliensis, Paracoccidioides lutzii, Paracoccidioides americana, Paracoccidioides venezuelensis, and Paracoccidioides restrepiensis (Muñoz et al., 2016; Turissini et al., 2017). The mycelium phase grows in the environment or in culture at 22°C. The fungus undergoes differentiation into the yeast form at 37°C in culture or within the host lungs and then disseminates into several organs (Brummer et al., 1993). The morphogenetic shift from mycelia to yeast deserves special attention once it denotes a virulence strategy of the pathogen to establish the infection (Rooney and Klein, 2002). The PCM treatment is arduous (Travassos et al., 2008) and the most common drugs used are based on azoles, sulfonamides and amphotericin B and they present a considerable degree of toxicity. Moreover, the period of treatment can be extended to years in some cases (Bocca et al., 2013).

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) is an essential glycolytic enzyme, also well known for its multitude of functions within cells. Therefore, it is categorized as a moonlight protein. GAPDH is widely explored virtually in all organisms and it has been shown to play vital roles in eukaryotic metabolism (Seidler, 2013), production of exporting vesicles (Tisdale, 2001), regulating proteins (Kornberg et al., 2010), related to the cell cycle and apoptosis (Tarze et al., 2007) and protecting cells against reactive oxygen species (ROS) (Hara et al., 2005). GAPDH was also found to be forming glycolytic multiprotein complexes in Saccharomyces cerevisiae (Brandina et al., 2006), with the heat shock protein (HSP) family in hepatitis B virus (Liu et al., 2009) and in plants (Bustos and Iglesias, 2003), besides being exported and contributing to Streptococcus pyogenes virulence (Jin et al., 2011). GAPDH has been intensely investigated in Paracoccidioides. It is differentially expressed in several approaches such as copper and zinc deprivation in P. lutzii (Parente et al., 2013; de Oliveira et al., 2014), up regulated in P. brasiliensis yeast cells (Rezende et al., 2011), acting as adhesin (Barbosa et al., 2006), present in the secretome of mycelia and yeast cells (Weber et al., 2012) and interacting with malate synthase (de Oliveira et al., 2013). As a moonlight macromolecule, GAPDH interacts with a large variety of proteins, forming multiprotein complexes that play unique roles within cells (Brandina et al., 2006; Ferreira et al., 2013). Interactions among proteins have been described by innumerous techniques such as immunoprecipitation (Ferreira et al., 2013), yeast two-hybrid and pull down (de Oliveira et al., 2013), tandem affinity purification (Xu et al., 2010), affinity chromatography (Wong et al., 2004) and Blue Native PAGE (BN-PAGE) (Brandina et al., 2006).

The study of protein-protein interactions (PPIs) adds knowledge related to protein function and helps to characterize protein complexes and the pathway they are involved, giving insight into new roles of proteins or adding more information to biological known processes. GAPDH establishes interactions with other proteins (Brandina et al., 2006; Liu et al., 2009; Kornberg et al., 2010; Araiza-Olivera et al., 2013; Ferreira et al., 2013) and due to the importance of GAPDH over cell dynamics, here we show an overview of PPIs of GAPDH in P. lutzii in order to explore its functions and participation in different metabolic processes.

Here, we show an overview of experimental GAPDH interactome in different phases of P. lutzii and an in silico validation of 18 proteins that bound to GAPDH through the experimental approaches. GAPDH interacts with 207 proteins in P. lutzii. We modeled the GAPDH 3-D structure, performed molecular dynamics and molecular docking in order to identify the interacting interface between GAPDH and the interacting proteins. Moreover, we found that molecular dynamics significantly improved GAPDH 3-D model and that GAPDH has four main regions of contact with interacting proteins. Several proteins bound to GAPDH in mycelium, transition and yeast phases are common to important pathways such as glycolysis and (tricarboxylic citric acid) TCA. The interaction between GAPDH and triose phosphate isomerase, enolase, 2-methylcitrate synthase and isocitrate lyase was validated by co-immunoprecipitation assay. GAPDH of P. lutzii is a moonlight protein and its interactome shows that GAPDH participate in a large variety of biological processes in the pathogen.

Materials and Methods

P. lutzii Mycelium and Yeast Culture Maintenance

Paracoccidioides lutzii (ATCC-MYA-826) was used in all experiments in the present study. This isolate is the most representative in the region where this work was performed and responsible for most of the PCM cases in central-west of Brazil. It was cultivated on Fava-Netto solid medium (1.0% w/v peptone, 0.5% w/v yeast extract, 0.3% w/v proteose peptone, 0.5% w/v beef extract, 0.5% w/v NaCl, 4% w/v glucose, and 1.4% w/v agar, pH 7.2) (Fava-Netto and Raphael, 1961) at 22°C for the growth of P. lutzii mycelium phase or 37°C for the growth of the yeast phase. The media was renewed every 7 days.

P. lutzii Mycelium to Yeast Transition

Paracoccidioides lutzii mycelium to yeast transition was performed in Fava-Netto liquid medium (Fava-Netto and Raphael, 1961) by a change in the culture temperature from 22°C to 37°C (Da Silva et al., 1994). P. Lutzii mycelium cells were previously grown in liquid medium for 48 h at 22°C before the temperature shift to 37°C. The samples were collected at 72 h of growth in order to perform morphological monitoring through a Neubauer chamber.

P. lutzii Protein Extraction Procedure

Mycelium cells were collected after 14 days of growth, mycelium-to-yeast transition cells after 3 days post temperature shift and yeast cells after 3 days of growth. Each sample phase was submitted to total protein extraction. Cells were collected by centrifugation at 10,000 × g for 10 min at 4°C and washed in sterile phosphate buffered saline (PBS; 1.4 mM KH₂PO₄, 8 mM Na₂HPO₄, 140 mM NaCl, and 2.5 mM KCl at pH 7.2). The cells were then resuspended in a buffer with 20 mM Tris–HCl at pH 8.8 and 2 mM CaCl₂ and added a protease inhibitor mix (GE Healthcare, Uppsala County, Uppsala, Sweden) in order to prevent protein degradation. The samples were conditioned in proper tubes with glass beads and applied to a bead beater apparatus (BioSpec, Bartlesville, OK, United States) for 5 cycles of 30 s each, in order to disrupt cells and expose proteins to the solution. Samples were kept on ice to avoid protein denaturation. Then, the samples were centrifuged at 10,000 × g for 15 min at 4°C. The supernatant was collected and the protein concentration was determined using Bradford reagent (Sigma-Aldrich, St. Louis, MO, United States). Bovine serum albumin (BSA) was used as standard.

P. lutzii GAPDH Heterologous Expression and Recombinant Protein Purification

Glyceraldehyde-3-phosphate dehydrogenase recombinant protein was obtained as described by Barbosa et al. (2006). Briefly, we used the TOPO-pET-100-GAPDH construct in order to produce recombinant GAPDH. The recombinant protein synthesis started by the addition of isopropyl-β-D-thiogalactoside (IPTG) (Sigma-Aldrich) to a final concentration of 1 mM. BL21 cells expressing recombinant P. lutzii GAPDH were incubated with 1 μ/mL lysozyme and underwent subsequent sonication on ice. For the recombinant GAPDH with his-tag at the N-terminal end purification, we used the nickel-nitrilotriacetic acid resin (Ni-NTA) (Invitrogen, Carlsbad, CA, United States). For 1 mL of lysate, we used 250 μL of resin. The lysate was incubated with the resin on ice under gentle agitation for 2 h. Then the resin was washed 5 times with native wash buffer (50 mM Na₂HPO₄, 20 mM imidazol, pH 8.0) in order to eliminate contaminants. The recombinant GAPDH bound to the resin was eluted with native elution buffer (50 mM Na₂HPO₄, 250 mM imidazol, pH 6.0). The purity and size of the protein were assessed by a 12% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) followed by Coomassie staining.

P. lutzii GAPDH Pull Down Assay

Nickel-nitrilotriacetic acid resin pull down assay is an affinity chromatography procedure that uses a bait protein, P. lutzii GAPDH in our case, immobilized and incubated with a protein source containing putative protein preys. The assay was performed in native conditions and the results included direct and indirect associations between bait and prey proteins. Firstly, recombinant GAPDH was immobilized onto Ni-NTA resin following the purification assay without the elution step. Then, 300 μg of P. lutzii cell lysate containing the protein extract was incubated for 3 h on ice and under gentle agitation. Next, the column containing bait and prey proteins was washed with native wash buffer 5 times to reduce unspecific interactions or contaminants and then, the complex bait-preys was eluted with native elution buffer. The eluted sample underwent tryptic digestion and the digested peptides were separated further via NanoUPLC-MS^E and analyzed using a nanoACQUITY system (Waters Corporation, Milford, Manchester, United Kingdom) in order to identify the proteins that possibly interacted with GAPDH.

A control sample was prepared similarly but 300 μg of P. lutzii protein extract was incubated with Ni-NTA without the immobilized GAPDH and the proteins identified in both experiments were excluded from the results.

Blue Native PAGE Approach

Blue Native-PAGE was performed according to the protocol established by Schägger et al. (1994) with some modifications. The samples loaded in the gels were P. lutzii mycelium, mycelium-to-yeast and yeast water-soluble proteins. The samples (300 μg of proteins) were dissolved in 10% (w/v) glycerol and 50 mM Bistris/HCl at pH 7.0. A 5–18% (w/v) polyacrylamide gradient gel was casted and a gel buffer (150 mM Bistris/HCl, 1.5 M aminocaproic acid at pH 7.0), a cathode buffer (50 mM tricine, 15 mM Bistris/HCl, 0.02% Coomassie blue G-250, at pH 7.0) and an anode buffer (50 mM Bistris/HCl at pH 7.0) were used to conduct native electrophoresis. The electrophoresis was performed on a vertical Hoefer SE600 ruby apparatus (GE Healthcare) at 12°C, with a starting voltage of 100 V until the loaded samples were inside the stacking gel, then a constant current limited to 300 V and 15 mA were applied.

Co-immunoprecipitation

The assay was performed according to Araiza-Olivera et al. (2013) with some modifications. Briefly, target antibodies were incubated with protein A sepharose 4B (Invitrogen, Waltham, MA, United States) for 3 h. Then, the solution antibodies and resin was incubated with 3 mg/mL of total protein extract for 3 h. Washing steps were performed with wash solution, 10 m mM tris at pH 7.4, 1% Triton and 1 mM EDTA (ethylenediaminetetraacetic acid). Elution was performed with 0.2 M glycine at pH 2.6 and equal amount of 0.2 M tris at pH 8.0. Eluted samples underwent SDS (sodium dodecyl sulfate)-PAGE followed by Western blot.

STRING Database Analysis

The STRING (Search Tool for the Retrieval of Interacting Genes/Proteins) database focuses on the prediction and unification of PPIs, direct and indirect associations, in order to cover even organisms where experiments are unavailable. Here, the accession numbers of identified proteins were analyzed for their interactions using STRING database version 9.1¹ (Szklarczyk et al., 2011). In order to perform the analysis we included active interaction sources such as text mining, experiments, databases, and co-expression with the highest confidence score (0.900).

In-Gel Enzymatic Digestion

The in-gel digestion was performed in accordance with (Rezende et al., 2011). The Coomassie blue-stained spots were manually excised from the gels, washed twice with MilliQ H₂O, dehydrated in 100 μL of acetonitrile (ACN) and dried in a speed vacuum (Eppendorf, Hamburg, Germany). Following, the reduction with 30 μL of 10 mM dithiothreitol (DTT) solution (GE Healthcare, Piscataway, NJ, United States), incubated in a dry bath at 56°C for 1 h and alkylated with 30 μL of 55 mM iodoacetamide (GE Healthcare, Piscataway, NJ, United States) for 45 min at room temperature in the dark. The supernatant was removed, and the gels were washed with 100 μL of a 25 mM ammonium bicarbonate solution by vigorous mixing for 10 min, and dehydrated twice in 100 μL of a 25 mM ammonium bicarbonate and 50% (v/v) ACN solution, vortexed for 5 min and centrifuged. The gel pieces were dried, and a 0.01 μg/μL trypsin solution (Promega, Madison, WI, United States) was added, followed by rehydration on ice at 4°C for 10 min. The supernatant was removed, and 25 μL of 25 mM ammonium bicarbonate solution was added to the gel pieces, followed by incubation at 37°C for 16 h. Then, the supernatant was placed into a clean tube. Fifty microliters of 50% (v/v) ACN, 5% (v/v) trifluoroacetic acid (TFA) solution was added to the gel pieces. The samples were vortexed for 10 min and the supernatant was then combined with the solution retrieved after digestion. The samples were dried in a speed vacuum, and the peptides were solubilized in 10 μL MilliQ H₂O. The concentration and purification was performed using a pipette tip (ZipTips^®C18 Pipette Tips, Milipore, Bedford, MA, United States).

Digestion of Protein Extracts

Enzymatic digestion of proteins was processed according to Murad et al. (2011) and Murad and Rech (2012) with some modifications. Briefly, 50 μg of total protein was added to 10 μL of 50 mM ammonium bicarbonate pH 8.5, in a microcentrifuge tube. Then 25 μL of RapiGEST^TM SF Surfactant (0.2% v/v) (Waters, Milford, MA, United States) was added, and the sample was vortexed and incubated in a dry bath at 80°C for 15 min. Following the reduction of disulfide bonds with 2.5 μL of 100 mM dithiothreitol (DTT) (GE Healthcare, Piscataway, NJ, United States) at 60°C for 30 min and alkylation of cysteine with 2.5 μL of 300 mM iodoacetamide (GE Healthcare, Piscataway, NJ, United States) for 30 mim at room temperature in the dark. The proteins were digested with 10 μL of trypsin 0.05 μg/μL (Promega, Madison, WI, United States) at 37°C in dry bath for 16 h. To precipitate the RapiGEST, 10 μL of a 5% trifluoroacetic acid solution (Sigma-Aldrich, St. Louis, MO, United States) was added, the samples were vortexed, and incubated for 90 min at 37°C in a dry bath, and centrifuged at 18,000 × g at 4°C for 30 min. The supernatant was dried in a speed vaccum apparatus (Eppendorf, Hamburg, Germany). All obtained peptides were suspended in 50 μL of a solution containing 20 mM of ammonium formate and 100 fmol/μL of phosphorylase B (Waters Corporation, Manchester, United Kingdom) (MassPREP^TM protein).

NanoLC-MS^E Acquisition

Nanoscale LC separation of tryptic peptides was performed using a ACQUITY UPLC^®M-Class system (Waters Corporation, United States) equipped with a XBridge^®Peptide 5 μm BEH130 C18 300 μm × 50 mm pre-column; Trap, 2D Symmetry^®5 μm BEH100 C18, 180 μm × 20 mm column and Peptide CSH^TM BEH130 C18 1.7 μm, 100 μm × 100 mm analytical reversed-phase column (Waters Corporation, United States). The peptides obtained from the digestion of proteins in the pull down assay were separated using a gradient of 11.4, 14.7, 17.4, 20.7, and 50% of ACN and 0.1% (v/v) formic acid, with a flow rate of 2000 μL/min. For samples obtained from the gel pieces, 4 μL full loop injections were initially transferred with an aqueous 0.1% formic acid solution to the pre-column. After desalting and pre-concentration, the peptides were eluted from the pre-column to the analytical column and separated with a gradient of 3–45% mobile phase B (0.1% formic acid in ACN) over 40 min at a flow rate of 2000 μL/min, followed by a 10 min rinse with 99.9% of mobile phase B.

The lock mass was used for calibration of the apparatus, using a constant flow rate of 0.5 μL/min at a concentration of 200 fmol of GFP ([Glu1]-Fibrinopeptide B human (Sigma-Aldrich, St. Louis, MO, United States). Mass spectrometry analysis was performed on a Synapt G1 MS^TM (Waters, United States) equipped with a nanoelectrospray source and two mass analyzers: a quadrupole and a time-of-flight (TOF) operating in TOF V-mode. Data were obtained using the instrument in the MS^E mode, which switches the low energy (6 V) and elevated energy (40 V) acquisition modes every 0.4 s. Samples were analyzed from three replicates.

Data Processing and Protein Identification

The mass spectrometer data obtained from the LC-MS^E analysis were processed using the ProteinLynx Global Server version 3.0.2 (Waters, Manchester, United Kingdom). The data were subjected to automatic background subtraction, deisotoping and charge state deconvolution. The processed spectra were searched against Paracoccidioides Pb01 protein sequences (Broad Institute)^2,3 . The mass error tolerance for peptide identification was under 50 ppm. The protein identification criteria also included the detection of at least two fragment ions per peptide, five fragments per protein and the determination of at least one peptide per protein; the identification of the protein was allowed with a maximum 4% false positive discovery rate in at least three technical replicate injections. The searches were performed with fixed modification of carbamidomethyl-C, and variable modifications were phosphorylation of serine, threonine and tyrosine. One missed cleavage site was allowed. The minimum repeat rate number for each protein in all replicates was two. Protein tables generated by ProteinLynx Global Server were merged and the dynamic range of the experiment was calculated using the software program MassPivot v1.0.1. The peptide and protein tables were compared using the Spotfire v8.0 software, suitable graphics were generated for all data as previously described (Murad et al., 2011). Microsoft Excel (Microsoft, United States) was used for table manipulations. UniProt⁴ and Pedant on MIPS⁵ databases were used for functional classification. NCBI database was employed for annotation of uncharacterized proteins⁶. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE (Vizcaíno et al., 2016) partner repository with the dataset identifiers PXD008253, PXD008308, PXD008309, and PXD008317.

Preparation of 3D Structures, Molecular Dynamics and Molecular Docking

The I-TASSER server⁷ modeled GAPDH and interacting proteins. The I-TASSER modeling is based on templates of homologous proteins with experimental structure available on the PDB database (Yang et al., 2015). All proteins under study had their three-dimensional structure assessed for quality using the MolProbity server⁸ (Cheng et al., 2007).

In order to obtain a more stable structure and similar to the native model, the target protein (GAPDH) was submitted to molecular dynamics (MD). The MD was performed using the software GROMACS 4.5.5, AMBER force field (ff99SB-ILDM), with the presence of explicit solvent (water TIP3P) in order to get a more stable protein structure (Pronk et al., 2013). In the first stage of MD, the overall free-energy was minimized in order to remove conflicting contacts and the process was completed only when the tolerance of 1000 kJ/mol was not exceeded. Seeking the balance of the thermodynamic variables of the system, a 100 ps simulation of NVT (amount of substance in moles, volume and temperature) and NPT (amount of substance in moles, pressure and temperature) simulation were applied. For the NVT phase, where there is variation in the pressure of the system, the temperature was set to 300 K by the thermostat and the velocities were calculated through Maxwell’s equations. For the NPT simulation, the volume is allowed to vary and the pressure was maintained constant by the Parrinello-Rahman barostat (Iannuzzi et al., 2003; Hutter, 2011).

Next, the protein was submitted to a simulation of 100 ns, 300 K, 1 atm and time interval of 2 fs, without restriction of the protein conformation (Pronk et al., 2013). MD analysis of GAPDH trajectories was performed by root mean square deviation (RMSD) in relation to the initial structure by the gromos algorithm (Daura et al., 2004). The g_cluster program (GROMACS package) was used to determine the conformations that were most frequent during the simulation. We used a cut-off of 0.3 nm in order to distinguish the conformational sets based on the RMSD profile. The cluster analysis and the RMSD allow the analysis of the protein profile throughout the simulation and the most representative conformational mode was selected to undergo molecular docking (Pronk et al., 2013).

We used the GRAMM-X protein-protein anchor server⁹ to determine the best protein complex conformations between GAPDH and interacting proteins (Tovchigrechko and Vakser, 2006). After the assembly of the complexes, the next step was to identify the amino acids involved in the interaction. We used the KFC2 server¹⁰ to recognize all residues from the interaction interface from GAPDH and interacting protein (Zhu and Mitchell, 2011).

We used the CoCoMAPS¹¹ server after the most frequent GAPDH residues involved in the interactions were identified. CoCoMAPS enables the analysis and visualization of the interface of interaction in protein complexes by making use of intermolecular contact maps in order to identify hot spot residues (Vangone et al., 2011). After identifying the interacting residues we used Pymol molecular visualization program to highlight the GAPDH regions most frequently involved in the interactions.

Identification of Nitrosylation Target Signal

We used the software GPS-SNO 1.0 (Xue et al., 2010) for the prediction of nitrosylation sites. A nitrosylation signal peptide comprises a short sequence of amino acids with a cysteine residue, normally placed in the middle of the sequence. As peptides with similar amino acid sequences may present similar structures, a substitution matrix is used in order to compute similarity between two peptide sequences. The results undergo clusterization steps and a score is generated along with a cut-off value and the determination in which cluster the signal peptide was identified.

Results and Discussion

Overview of GAPDH Interacting Proteins

The PPI assays performed here identified 207 GAPDH interacting proteins. Among those interactions, 30 (14.5%) occurred in mycelium, 69 (33.3%) in mycelium-to-yeast transition and 108 (52.2%) in P. lutzii yeast cells. Some of the GAPDH interactions were identified in more than one morphological phase of the fungus (Figure 1). The difference in GAPDH binding patterns among those phases results from two main factors. Firstly, certain proteins are differentially expressed according to the life cycle of the pathogen (Rezende et al., 2011) and secondly, some of the interactions may be lost during experimental procedures (Hayes et al., 2016).

FIGURE 1

Figure 1. Venn diagram of GAPDH interacting proteins. We found 207 proteins interacting with GAPDH through two different approaches and performed independently for the three stages of P. lutzii life cycle. The numbers on the diagram represent the number of proteins interacting with GAPDH in mycelium (M), mycelium-to-yeast transition (T) and yeast (Y) P. lutzii phases.

According to the pull down assay, 18 proteins from P. lutzii mycelium phase bound to the recombinant GAPDH (Supplementary Table S1). Their biological classification includes metabolism (28%), energy (28%), protein synthesis (11%), protein fate (5%), cell rescue, defense and virulence (5%), biogenesis of cellular components (5%) and hypothetical proteins (18%) (Figure 2). The multiprotein complex identified by BN-PAGE in mycelium phase, included proteins related to metabolism (14%), energy (29%), protein synthesis (7%) and protein fate (50%) (Supplementary Table S2). Interestingly, enolase (ENO) and isocitrate lyase (ICL) interaction with GAPDH was found in both PPI approaches. We performed a co-immunoprecipitation assay to validate the complex formed by GAPDH and ENO and ICL (Figure 3).

FIGURE 2

Figure 2. Schematic model of the metabolic processes involving GAPDH through PPIs in P. lutzii mycelium phase. All proteins shown in the figures interacted with GAPDH as a result of the chromatographic and BN-PAGE approaches. GAPDH (glyceraldehyde 3-phosphate dehydrogenase), SNO (nitrosylation), CS (citrate synthase), ACO (aconitase), ICL (isocitrate lyase), MCS (2-methylcitrate synthase), MCD (2-methylcitrate dehydratase), HSP60 (heat shock protein 60), HSP70 (heat shock protein 70), NQO4 (NAD(P)H:quinone oxidoreductase type IV) and ENO (enolase). Asterisks indicate proteins identified as GAPDH partners by both the chromatographic and BN-PAGE approaches.

FIGURE 3

Figure 3. Validation of GAPDH interacting proteins by co-immunoprecipitation assay. The figure shows the confirmation of GAPDH being complexed with proteins that we found as a result of chromatographic and co-immunoprecipitation assays TPI (triose phosphate isomerase), ENO (enolase), MCS (2-methylcitrate synthase), ICL (isocitrate lyase) and C-(negative control).

The pull down assay, regarding mycelium-to-yeast transition, identified 54 proteins interacting with GAPDH (Supplementary Table S3). They were related to metabolism (19%), energy (7%), cell cycle and DNA processing (11%), transcription (7%), protein synthesis (32%), protein fate (9%) and hypothetical proteins (15%). For the BN-PAGE assay (Supplementary Table S4), GAPDH interacting proteins were related to metabolism (33%), energy (17%), cell cycle and DNA processing (22%), protein fate (17%) and hypothetical proteins (11%) (Figure 4). NAD(P)H:quinone oxidoreductase (NQO4) and an isoform of HSP70 interacted with GAPDH in both assays (Supplementary Tables S3, S4). We found six proteins bound to GAPDH that were up-regulated during mycelium-to-yeast transition in comparison to mycelium and yeast phases in P. brasiliensis (Rezende et al., 2011) and those proteins were malate dehydrogenase (MDH), ATP synthase β, two variants of HSP70, HSP90, alcohol dehydrogenase (ADH) and adenosylhomocysteinase (Figure 2). These findings show that GAPDH is likely to play important roles during Paracoccidioides spp. mycelium-to-yeast transition. Due to its probable role in substrate channeling, the interaction between GAPDH and triose phosphate isomerase (TPI) (Sweetlove and Fernie, 2018) were validated by a co-immunoprecipitation assay (Figure 3).

FIGURE 4

Figure 4. Schematic model of the metabolic processes involving GAPDH through PPIs in P. lutzii mycelium-to-yeast transition. All proteins shown in the figures interacted with GAPDH. GAPDH (glyceraldehyde-3-phosphate dehydrogenase), SNO (nitrosylation), MDH (malate dehydrogenase), TPI (triosephosphate isomerase), NQO4 (NAD(P)H:quinone oxidoreductase type IV), tubulin α-1 chain (TUB), HSP70 (heat shock protein 70), histone H2A (H2A) and ADH (alcohol dehydrogenase). The figure also shows the proteins 14-3-3, Rad24 and ATP synthase β. Asterisks indicate proteins identified as GAPDH partners by both the chromatographic and BN-PAGE approaches.

The pull down approach identified 88 target proteins interacting with GAPDH in the yeast phase. We found proteins related to metabolism (23%), energy (18%), cell cycle and DNA processing (8%), protein synthesis (39%), protein fate (11%) and hypothetical proteins (1%) (Supplementary Table S5). For the BN-PAGE approach (Supplementary Table S6), we found 26 GAPDH interacting proteins related to metabolism (19%), energy (23%), cell cycle and DNA processing (4%), protein synthesis (11%), protein fate (12%), cellular transport (4%), signal transduction (4%) and hypothetical proteins (23%) (Figure 5). GAPDH interaction with ICL, 2-methylcitrate synthase (MCS), 2-methylcitrate dehydratase (MCD), acetyl-CoA acetyltransferase (ACA) and HSP90 were found in both PPI approaches performed (Supplementary Tables S5, S6).

FIGURE 5

Figure 5. Schematic model of the metabolic processes involving GAPDH through PPI in P. lutzii yeast phase. All proteins shown in the figures interacted with GAPDH. GAPDH (glyceraldehyde-3-phosphate dehydrogenase), SNO (nitrosylation), ACAD (acyl-CoA dehydrogenase), ECHD (enoyl-CoA hydratase), KAT (3-ketoacyl-CoA thiolase), CS (citrate synthase), IDH (isocitrate dehydrogenase), ICL (isocitrate lyase), KGD (dihydrolipoyl dehydrogenase), SDHA (succinate dehydrogenase flavoprotein subunit), MDH (malate dehydrogenase), MCS (2-methylcitrate synthase), MCD (2-methylcitrate dehydratase), TPI (triosephosphate isomerase), PGK (phosphoglycerate kinase), PK (pyruvate kinase), PDH (pyruvate dehydrogenase E1 component subunit β), PDC (pyruvate decarboxylase), ADH (alcohol dehydrogenase), ALDH (aldehyde dehydrogenase), ACA (acetyl-CoA acetyltransferase), AHTL (o-acetylhomoserine thiol-lyase), CTL (cystathionine γ-lyase), AAT (aspartate aminotransferase), IVCD (isovaleryl-CoA dehydrogenase), ASS (argininosuccinate synthase), CIMS (cobalamin-independent methionine synthase), OND (oxidoreductase 2-nitropropane dioxygenas), LTD (L-threonine 3-dehydrogenase), CCP (cytochrome c peroxidase) and NQO4 (NAD(P)H:quinone oxidoreductase type IV). The figure also shows the proteins HSP70 (heat shock protein 70) and ATP synthase β. Asterisks indicate proteins identified as GAPDH partners by both the chromatographic and BN-PAGE approaches.

Among the proteins that interacted with GAPDH, 15 were up regulated in the yeast morphological phase of P. brasiliensis compared to the mycelium and transition phases (Supplementary Table S5) (Rezende et al., 2011), showing that GAPDH has important roles regarding specific activities played by P. lutzii yeast cells in a variety of biological processes (Figure 5). In addition, interaction between GAPDH and MCS was validated by co-immunoprecipitation (Figure 3).

GAPDH Integrates Metabolons in Substrate Channeling

Glyceraldehyde-3-phosphate dehydrogenase is a moonlight protein with a multitude of functions (Barbosa et al., 2006; Rezende et al., 2011; Weber et al., 2012; de Oliveira et al., 2013, 2014; Parente et al., 2013) due to interactions with other proteins. Here, GAPDH interacts with several proteins related to energy metabolism via glycolysis, TCA, respiration, glyoxylate and methylcitrate pathways. It has been shown that mycelium P. lutzii cells use aerobic metabolism to obtain energy while the yeast cells use anaerobic metabolism (Felipe et al., 2005). We found GAPDH interacting with proteins from aerobic metabolism such as ENO, citrate synthase (CS), aconitase (ACO) and (NQO4) with proteins from anaerobic metabolism, such as ICL. Other studies demonstrate that among several moonlight function of GAPDH is to be a sensor of cellular stress (Colell et al., 2009; Tristan et al., 2011) and control of the aerobic and anaerobic metabolic shift (Schuppe-Koistinen et al., 1994).

Glyceraldehyde-3-phosphate dehydrogenase interacts with four glycolytic enzymes, ENO (Supplementary Tables S1, S2), TPI (Supplementary Tables S4, S5), pyruvate kinase (PK) and phosphoglycerate kinase (PGK) (Supplementary Table S5). The interactions between GAPDH and ENO and TPI were identified in the affinity chromatography assay (Supplementary Table S1), BN-PAGE (Supplementary Table S2), co-immunoprecipitation (Figure 2) and in the STRING PPI database. GAPDH interacting with PK was identified in the affinity chromatography assay and in the STRING database. Several PPI studies corroborate our findings regarding GAPDH interacting with glycolytic enzymes ENO (Kumar et al., 2004; Ferreira et al., 2013), TPI (Stephan et al., 1986; Liaud et al., 2000; Kumar et al., 2004; Graham et al., 2007; Araiza-Olivera et al., 2013; Menard et al., 2014), PGK (Prabhakar et al., 1999; Tomokuni et al., 2010; Araiza-Olivera et al., 2013; Wan et al., 2015) and PK (Kumar et al., 2004; Araiza-Olivera et al., 2013; Ferreira et al., 2013; Wan et al., 2015; Das et al., 2016). Coupling catalytic function of enzymes, such as GAPDH and other glycolytic enzymes, affords optimization of regulation and combined functionality of the reactions. They associate in order to accomplish their metabolic functions (Auger et al., 2015) and to channel pyruvate toward TCA cycle, establishing a tangible connection between glycolysis and TCA.

Protein-protein interactions studies in S. cerevisiae through BN-PAGE and immunoprecipitation approaches resulted in 13 proteins associated in a macromolecular complex (Brandina et al., 2006). They found three proteins from the glycolytic pathway, including ENO and GAPDH. PPI studies on Escherichia coli found GAPDH interacting with ENO and other glycolytic enzymes (Kumar et al., 2004; Ferreira et al., 2013). In S. cerevisiae, glycolytic enzymes associate among themselves and with actin filaments in order to stabilize and protect enzymes in a glycolytic metabolon (Araiza-Olivera et al., 2013). This is the first time a metabolon is identified in P. lutzii and we found GAPDH interacting with actin cytoskeleton protein (VIP1) in the pull down assay (Supplementary Table S3), which strengthens the idea of the existence of a glycolytic metabolon in P. lutzii. The interaction between GAPDH and TPI is a coherent example of substrate channeling because the product of the TPI reaction is the substrate of GAPDH. Metabolons increase the concentration of enzymatic intermediates locally and decrease the amount of enzyme needed to maintain the intermediate flux in a pathway (Sweetlove and Fernie, 2018). Hence, the final product is directed in a scaffold manner to specific subcellular spots and the metabolon prevents loss of those intermediates.

Glyceraldehyde-3-phosphate dehydrogenase interacts with several enzymes of TCA. We also found GAPDH bound to two enzymes from the pyruvate dehydrogenase complex (pyruvate dehydrogenase E1 component subunit β and dihydrolipoyl dehydrogenase; Supplementary Table S5). These enzymatic complex converts pyruvate into acetyl-CoA and the multiprotein complex formed with GAPDH could drive acetyl-CoA into its destiny by substrate channeling. GAPDH and TCA enzymes have been found to form protein complexes with glycolytic enzymes and mitochondrial membrane carriers in S. cerevisiae (Brandina et al., 2006).

We found GAPDH bound to MDH (Figure 4) and the presence of MDH among up regulated proteins in P. brasiliensis mycelium-to-yeast transition is noteworthy (Rezende et al., 2011). Another study detected, through a pull down assay, potential interactions between Porphyromonas gingivalis proteins and Streptococcus oralis recombinant GAPDH (Maeda et al., 2013). MDH was among the proteins investigated. Moreover, they showed that P. gingivalis MDH is a surface associated protein and the interaction between the recombinant GAPDH and MDH inhibited P. gingivalis and S. oralis biofilm formation. In addition, GAPDH is on the cell surface of several organisms, including P. gingivalis where it is related to the regulation of host invasion (Tribble et al., 2006) and P. brasiliensis where it is involved in fungal adhesion to extracellular proteins (Barbosa et al., 2006; Weber et al., 2012). These results indicate that an interaction between GAPDH and MDH may play a role in the host-pathogen relationship.

Glyceraldehyde-3-phosphate dehydrogenase interacted with isocitrate dehydrogenase (IDH) (Supplementary Table S5). IDH catalyzes the conversion of isocitrate into α-ketoglutarate and is regulated by acetylation. A deacetylase family protein named sirtuin (SIRT) stimulates its activity (Someya et al., 2010). Interestingly, SIRT is negatively regulated by GAPDH through nitrosylation (Kornberg et al., 2010). GAPDH and IDH interaction has been shown by other study and they found GAPDH and IDH in a multiprotein complex that used thioredoxin as protein bait (Kumar et al., 2004). Another protein found in this thioredoxin multiprotein complex was succinate dehydrogenase flavoprotein subunit (SDHA) (Supplementary Table S5), which corroborates our findings regarding the pull down assay.

The proteins ICL, MCS (Supplementary Tables S1, S5, S6) and MCD (Supplementary Tables S2, S5, S6) interactions with GAPDH were identified in both PPI approaches. The former takes part in the glyoxylate cycle and the other two enzymes belong to the methylcitrate pathway. ICL and MCD were up regulated in P. lutzii yeast cells, which utilize anaerobic metabolism driving two-carbon sources to feed glyoxylate cycle and three-carbon sources, such as propionate, to feed methylcitrate cycle. GAPDH was found in trypanosome peroxisomes (Opperdoes and Borst, 1977) and here we found GAPDH interacting with peroxin-19. Peroxins takes part in the recognition of proteins with peroxisomal or glyoxysomal targeting signals and consequently mark them to be transported into peroxisomes. This explains that P. lutzii GAPDH may be found in peroxisomes allowing its interaction with proteins from the glyoxylate cycle. Moreover, it has been shown through affinity chromatography and in silico approaches that glyoxylate cycle enzymes interact with GAPDH in P. lutzii (de Oliveira et al., 2013).

Regarding cellular respiration, we found GAPDH interacting with NQO4 (Supplementary Tables S1, S3–S5), ATP synthase β (Supplementary Tables S3, S5) and cytochrome c peroxidase (Supplementary Table S5), showing a probable regulation of ATP intracellular levels. Ferreira et al. (2013) showed that GAPDH may have high affinity for ATP-binding domains. They found GAPDH interacting with ATP-binding proteins and ATP synthase subunits α and β. The latter was also found in studies regarding multiprotein complexes related to the protein thioredoxin on E. coli (Kumar et al., 2004). The relation between GAPDH and cellular energy levels has been investigated in S. pyogenes where mutants incapable of export GAPDH to the cellular surface, increased the intracellular ATP concentration significantly compared to wild strains. Moreover, they showed that up regulation of GAPDH also led to up regulation of ATPases genes (Jin et al., 2011). A study performed on E. coli showed through immunoprecipitation and pull down assays that GAPDH interacts with ATP synthase β (Ferreira et al., 2013). They concluded that the interaction between those proteins indicates an additional process for the regulation of ATP dependent processes. In S. cerevisiae, NQO4 and GAPDH were in the same multiprotein complex (Gavin et al., 2002) and another group found GADPH interacting with cytochrome c oxidase in bovine heart tissue (Ramzan et al., 2013). These interactions altogether show that GADPH influences intracellular ATP levels.

GAPDH Participates in Nitrosylation Processes

Glyceraldehyde-3-phosphate dehydrogenase takes part in protein regulation via gasotransmitters (Mustafa et al., 2009). Hydrogen sulfide (H₂S), nitric oxide (NO) and carbon monoxide (CO) mediate cellular processes regulating several protein functions (Nguyen et al., 2011). Some proteins interact via thiol and disulfide bonds, normally at cysteine residues. A study on gasotransmitters showed that around 25% of cellular GAPDH is sulfhydrated and exerts functions apart from the glycolytic pathway (Hara et al., 2005). This is one of the reasons why GAPDH has several partners, being categorized as a moonlight protein. GAPDH acts as a nitrosylase (Hara et al., 2005) and also undergo nitrosylation (Kohr et al., 2014), increasing considerably the number of interacting partners. Interestingly, many of the proteins that were identified in the present work as interacting with GAPDH has a target signal for nitrosylation (Gould et al., 2013) (Supplementary Table S7). To support the hypothesis that GAPDH is a nitrosylase of mitochondria proteins, it has been shown that nitrosylated proteins that are expressed when GAPDH is overexpressed were down regulated after GAPDH knockdown in heart cells of mice (Kohr et al., 2014).

The nitrosylation features a covalent addition of nitric oxide to a cysteine thiol in a target protein and exerts regulatory function from bacteria to mammalian cells (Anand and Stamler, 2012). Co-immunoprecipitation experiments in mice cardiac cells showed that GAPDH regulates TCA enzymes (Kohr et al., 2014), including ACO that we found to be one of the GAPDH interacting protein in P. lutzii mycelium cells (Figure 2). Another study performing the same methodological experiment on S. cerevisiae also found GAPDH and TCA enzymes in the same complex (CS and MDH) (Kohr et al., 2011; Gould et al., 2013), showing that interactions among those enzymes must be playing important roles in regulatory pathways, metabolite channeling and maintaining a link between glycolysis and TCA (Brandina et al., 2006). In addition, the glyoxysome enzyme, ICL, undergoes nitrosylation and it has been shown that its activity is reversibly reduced after nitrosylation (Morisse et al., 2014).

The presence of GAPDH in mitochondria is related to its potential function as a nitrosylase of several mitochondrial protein targets (Kohr et al., 2014). We found GAPDH interacting with proteins from β-oxidation, among them acyl-CoA dehydrogenase (ACAD) (Doulias et al., 2013), 3-ketoacyl-CoA thiolase (KAT) (Supplementary Table S5) and enoyl-CoA hydratase (ECHD) (Supplementary Tables S4, S6) (Piantadosi, 2012). Nitrosylation also targets proteins related to amino acid metabolism, such as the gasotransmitter protein argininosuccinate synthase (ASS). We found this protein interacting with GAPDH in the pull down approach (Supplementary Table S5) and a large scale PPI study also found ASS bound to GAPDH (Wan et al., 2015). The biosynthesis pathway of the amino acid arginine is regulated by the nitrosylation of ASS, which inhibits the activity of this enzyme. The NO from the nitrosylation process acts as an important signaling molecule, and the regulation takes place at cysteine residues present at the surface of the interacting proteins (Hao et al., 2004).

Cell cycle proteins interacted with GAPDH in both approaches employed (Supplementary Tables S3, S4). The protein 14-3-3 has a multitude of functions, including regulation of signaling pathways and cell cycle. GAPDH and the protein 14-3-3 interacted in the BN-PAGE approach (Supplementary Table S4). The interaction between GAPDH and 14-3-3 proteins have been shown in several organisms such as plants (Bustos and Iglesias, 2003), human (Pozuelo Rubio et al., 2004) and hydra (Pauly et al., 2007). In the latter, this interaction increased sensitivity to inhibition by ATP, thus being regulated, which maintain safe levels of energy within cytoplasm (Pauly et al., 2007). In P. brasiliensis, 14-3-3 protein has been reported to be an adhesin (Assato et al., 2015), to be related to apoptosis (Silva Jde et al., 2015) and to be involved in fungus pathogenesis (Marcos et al., 2016). All those activities are also played by GAPDH (Barbosa et al., 2006; Tarze et al., 2007) in Paracoccidioides spp. and other microorganisms such as yeast (Carmona-Gutierrez et al., 2010).

Glyceraldehyde-3-phosphate dehydrogenase is present in the nucleus of cells (Sawa et al., 1997) in several organisms and one of the many functions played by this protein is DNA repair (Ferreira et al., 2013). We found GAPDH interacting with DNA damage checkpoint protein rad24 (Supplementary Table S4), which is related to DNA repair prior to mitosis. Some studies have classified rad24 as belonging to the 14-3-3 protein family (Ozoe et al., 2002). GAPDH, rad24 and the protein 14-3-3. We also found GAPDH interacting with histone H2A in the pull down approach (Supplementary Tables S3, S5). Histone H2A is also related to DNA repair when phosphorylated (Jakob et al., 2011). Several studies reported the relationship between GADPH and histone H2A expression, which depends on a redox status that could be guaranteed by GAPDH in multiprotein complexes (Shi and Shi, 2004; Dai et al., 2008; He et al., 2013).

Here, GAPDH interacted with ADH and aldehyde dehydrogenase (AHDH) (Supplementary Tables S4, S6). ADH acts in the interconversion of ethanol to acetaldehyde, producing NADH, while GAPDH glycolytic reaction promotes the reduction of NAD⁺ to NADH. Thus, both proteins within a multiprotein complex are able to provide or deprive the reducing equivalent NADH in response to metabolic changes. ADH is immunogenic in several microorganisms, including Candida albicans (Bertram et al., 1996) and P. brasiliensis (Dantas et al., 2009). GAPDH and ADH have been found in the same protein complex in large-scale PPI approaches (Gavin et al., 2002). A study performed on S. cerevisiae found GAPDH and ADH interacting in two complexes along with other proteins that we have also identified in our approaches, such as proteins from the cell division control (Supplementary Table S3) and the HSP (Supplementary Tables S3, S4) family (Gavin et al., 2002). A different large-scale study on Drosophila melanogaster identified 556 protein complexes and GAPDH and ADH were among the interacting proteins (Guruharsha et al., 2011).

The interaction between GAPDH and ADH or GAPDH and ALDH is found in other organisms as well (Gavin et al., 2002; Wong et al., 2004; Guruharsha et al., 2011; Ferreira et al., 2013). ADH is overexpressed in P. brasiliensis yeast phase undergoing anaerobic conditions (Rezende et al., 2011), which leads to a metabolic shift and ethanol production via pyruvate (Felipe et al., 2005). Ethanol is used to damage tissue and facilitates host invasion besides producing energy (Grahl et al., 2011; Luong et al., 2015). ALDH catalyzes the oxidation of aldehydes and participates in detoxification of such compounds (Marchitti et al., 2008). Through the activity of ALDH, pathogens evade osmotic, oxidative and other cellular stressors besides host immune response (Singh et al., 2004). The interaction between GAPDH and ALDH may optimize host invasion and protection against stress conditions (Abdul et al., 2018).

We found several HSPs complexing with GAPDH in P. lutzii (Supplementary Tables S2–S6). HSPs act as molecular chaperones and control protein biogenesis and correct folding (Craig et al., 1994). The protein HSP60 is mostly mitochondrial but is present in the cytoplasm in lower levels. It is related to the correct folding of proteins, stabilization of unfolded proteins and assembly of multimeric complexes (Soltys and Gupta, 1999). HSP70 and its isoforms bind to incompletely folded proteins, such as polypeptides on ribosomes and proteins that undergo transportation from cytosol into mitochondria, such as GAPDH (Murphy et al., 2014). HSP90 acts on modulation of protein activity (Craig et al., 1994).

The BN-PAGE approach was performed to study the HSP machinery involved in hepatitis B infection and it was found a protein complex with GAPDH, HSP60, HSP70, and HSP90 (Liu et al., 2009). The authors stated that proteins presented in the complexes they analyzed influenced the viral infection. HSPs play a key role maintaining the protein complex functional and HSP 70 and 90 could be candidates for antiviral therapy once they are related to the host-pathogen interaction (Liu et al., 2009).

Molecular Dynamics Improved GAPDH 3-D Model

After predicting three-dimensional structure of GAPDH through the I-TASSER server, the stereochemistry of the three-dimensional model were assessed by the Ramachandran plo (Supplementary Figures S1A,B). We performed GAPDH molecular dynamics in 100 ns. Supplementary Figure S1A shows that before the MD 77.1% residues were within favored regions and 90.8% of the residues were within allowed regions. The figure also shows that 9.22% of the amino acids in the initial conformation were at unfavorable regions. Supplementary Figure S1B shows that after MD, 89.0% residues were within favored regions and 97.3% of the residues were within allowed regions. Moreover, only 2.7% of the residues were within unfavorable regions and it is mostly likely that those regions are not present in the binding and catalytic site of GAPDH and may refer to unconserved regions. The number of residues with atoms in unfavorable angular positions decreased sharply. Before the MD, we had 31 residues in a forbidden position and after the MD, this number lowered to 9. Ramachandran plots consistently showed that GAPDH structure is refined after MD (Supplementary Figures S1A,B).

Protein stability is an important trait for protein interactions. Hydrogen bonds, hydrophobic and electrostatic interactions are key to determine the stability of a protein and a protein complex. Rational design of peptides and drugs in order to inhibit PPIs tries to improve stabilizing interactions between protein and ligand. We used RMSD to plot the analyses of MD trajectories to establish the equilibration period and quality of the biomolecular simulations besides clustering similar conformations. For RMSD at a value of 0.4 nm, the variation of the simulations becomes more stable (Supplementary Figure S1C). According to the trajectories of MD simulations, we identified at least seven conformational clusters. Cluster analysis showed that the structure of the first three clusters had a more stable conformation. Moreover, the first cluster remained stable for a longer period than the other clusters (Supplementary Figure S1D). The Molprobity score of the total structure improved from 2.78 to 1.62, and the clash-score from 6.44 to 0. The sharp reduction of those scores indicate the number of structurally superimposed atoms in the main GAPDH model. Overall, we observed a significant improvement of the GAPDH structure model after the MD (Supplementary Figure S1).

GAPDH Has Four Main Regions of Contact With Interacting Proteins

We proposed 36 GAPDH complexes (Supplementary Figure S2). For each GAPDH partner, we selected two conformations with the lowest free energy scores. Supplementary Figure S2 shows GAPDH in gray, the interaction interface in dark-blue and the GAPDH partner proteins in different colors to highlight the two most stable conformations of the interaction. We chose the GAPDH interacting partners for the in silico validation according to the main biological function categories they participate or according to their role as regulatory proteins. There are 11 proteins related to the energy metabolism (Supplementary Figures S2A–K), three to amino acid metabolism (Supplementary Figures S2L–N), two related to virulence (Supplementary Figures S2O,P), one to fatty acid metabolism (Supplementary Figure S2Q) and one related to protein fate and complex stabilization (Supplementary Figure S2R).

Molecular docking provided the prediction of ligand orientation and possible conformations of the interacting proteins. There are at least two steps to carry out the molecular docking process. Firstly, a conformational search of the ligand proteins within a well-defined grid box around the interacting surface of GAPDH to represent different possible conformations. Then, scores of each conformation are checked and the first two scores with the lowest free energy were selected. KFC-2 server analyzes the residues within the interaction interface between GAPDH and an interacting protein to identify the contact residues. We found that four regions define the main PPIs among all the 36 analyzed complexes (Figure 6A) more frequently. The areas represented in red are the regions that contain GAPDH residues that interact with a protein partner more frequently. The identified interactions may occur in vivo provided GAPDH is found in the same subcellular location of its interacting partners and most of them present nitrosylation signal. The prediction of protein interacting sites has been performed successfully by several similar approaches (Aumentado-Armstrong et al., 2015; Laine and Carbone, 2015; Maheshwari and Brylinski, 2015).

FIGURE 6

Figure 6. Frequent contact regions of GAPDH protein-protein interaction. (A) GAPDH (gray) and its four more frequent contact regions (red) of interaction with partner proteins. (B) The graphic shows the frequency of interaction for 18 GAPDH protein partners. Each letter and color represents a different protein that interacted with GAPDH in silico and the green bars represent the overall of the interactions in each residue.

Figure 6B represents an overview of the contact preference regions for all of the 18 interacting proteins selected for the present in silico study. The graphic shows that the GAPDH interacting region near residue 202 has the most frequent hits of interaction followed by the region 239 and 313. A large-scale study analyzed the contact preferences in 621 protein-protein complexes. They found out that those residues in the binding PPI sites are governed by hydrophobic contacts that help to maintain a stable complex conformational structure (Glaser et al., 2001). The contact preference residues on GAPDH PPI interface establish hydrophobic interactions more frequently but there are also non-polar and polar interactions involved.

Conclusion

Here, we showed an overview of experimental GAPDH interactome in different phases of P. lutzii. Several proteins bound to GAPDH from mycelium (Figure 2), transition (Figure 4) and yeast (Figure 5) phases are common to important pathways such as glycolysis and TCA. We performed affinity chromatography and BN-PAGE in order to identify GAPDH partners. We found 207 proteins interacting with GAPDH, which reflects its broad spectrum of functionality and its characterization as a moonlight protein. We performed a co-immunoprecipitation assay to validate the complex formed by GAPDH with TPI, ENO, ICL and MCS. The former has a probable role in substrate channeling, the others were found in both pull down and BN-PAGE approaches. We found GAPDH participating in complexes with proteins of specific pathways, indicating the existence of a glycolytic and a TCA metabolon in P. lutzii. GAPDH interacted with several proteins that undergoes regulation by nitrosylation, indicating that this proteins in P. lutzii is involved in several metabolic process such as sensitivity to ATP intracellular levels, stabilization of multiprotein complex conformation with HSPs and regulation of mitochondrial and nuclear proteins. Moreover, in silico analysis showed that GAPDH has four main regions of interaction with other proteins, reflecting its conservation across evolution. Identifying GAPDH interacting proteins may shed light on the main regions of this moonlight protein that interacts with other protein partners. Finding non-conserved regions to Homo sapiens GAPDH can drive the search for new therapeutic approaches against PCM. Our next step is to perform the screening of new inhibiting compounds and design small molecules (peptides) to modulate pathogen GAPDH.

Data Availability

The datasets generated for this study can be found in the ProteomeXchange Consortium via the PRIDE, PXD008253, PXD008308, PXD008309, and PXD008317.

Author Contributions

KS contributed to the chomatography, BN-PAGE, proteomics, in silico analysis, and the writing of the manuscript. RL contributed to the in silico analysis and the writing of the manuscript. LB and PL contributed to the proteomics and writing of the manuscript. TC and SC contributed to the BN-PAGE and writing of the manuscript. RS contributed to the in silico analysis and writing of the manuscript. CS contributed to the writing and editing of the manuscript. MP contributed to the experimental procedures, and the writing and revision of the manuscript.

Funding

This work was performed at Universidade Federal de Goiás and supported by MCTI/CNPq (Ministério da Ciência e Tecnologia/Conselho Nacional de Desenvolvimento Científico e Tecnológico), FNDCT (Fundo Nacional de Desenvolvimento Científico e Tecnológico), FAPEG (Fundação de Amparo à Pesquisa do Estado de Goiás), CAPES (Coordenação de Aperfeiçoamento de Pessoal de Nível Superior), Finance Code 001, FINEP (Financiadora de Estudos e Projetos), INCT-IF (Instituto Nacional de Ciência e Tecnologia para Inovação Farmacêutica), and PRONEX (Programa de Apoio a Núcleos de Excelência). Additionally, KS and RL were supported by fellowship from CAPES.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.01537/full#supplementary-material

FIGURE S1 | Molecular Dynamics result evaluation. Ramachandran plots showing the analysis of residues in allowed, favored and forbidden regions (A) before and (B) after GAPDH molecular dynamics. The light blue regions are favorable and the dark blue regions are allowed. Regions outside blue lines are forbidden regions and may characterize unconserved regions of GAPDH. (C) RMSD trajectories over 100 ns of GAPDH MD simulation. (D) The cluster plot shows seven sets of conformations along the GAPDH MD trajectories and the first cluster shows the most stable conformation.

FIGURE S2 | Schematic representation of GAPDH interaction with 18 protein partners. Each image shows two conformations of GAPDH (gray), its protein partner (different colors) and the interacting interface (dark-blue). (A) malate dehydrogenase; (B) succinate dehydrogenase; (C) triosephosphate isomerase; (D) phosphoglycerate kinase; (E) cytochrome c peroxidase; (F) 2-methylcitrate synthase; (G) 2-methylcitrate dehydratase; (H) aconitase; (I) pyruvate kinase; (J) isocitrate lyase; (K) enolase; (L) L-threonine 3-dehydrogenase; (M) ribose- phosphate pyrophosphokinase; (N) fumarylacetoacetase; (O) alcohol dehydrogenase; (P) aldehyde dehydrogenase; (Q) enoyl-CoA hydratase; (R) HSP70.

TABLE S1 | Potential GAPDH target proteins identified in P. lutzii mycelium through a pull down assay.

TABLE S2 | Potential GAPDH target proteins identified in P. lutzii mycelium phase through BN-PAGE.

TABLE S3 | Potential GAPDH target proteins identified in P. lutzii mycelium to yeast transition through a pull down assay.

TABLE S4 | Potential GAPDH target proteins identified in P. lutzii mycelium to yeast transition through BN-PAGE.

TABLE S5 | Potential GAPDH target proteins identified in P. lutzii yeast phase through a pull down assay.

TABLE S6 | Potential GAPDH target proteins identified in P. lutzii yeast phase through BN-PAGE.

TABLE S7 | Prediction of potential GAPDH target proteins with S-nitrosylation sites in P. lutzii.

Footnotes

1. ^http://string.embl.de
2. ^http://www.broadinstitute.org/annotation/genome/paracoccidioides_brasiliensis/Multiome.html
3. ^http://www.uniprot.org/uniprot/?query=macrophage&sort=score
4. ^http://www.uniprot.org
5. ^http://mips.helmholtz-muenchen.de/funcatDB/
6. ^https://www.ncbi.nlm.nih.gov/
7. ^http://zhanglab.ccmb.med.umich.edu/I-TASSER/
8. ^http://molprobity.biochem.duke.edu/
9. ^http://vakser.compbio.ku.edu/resources/gramm/grammx/
10. ^https://mitchell-lab.biochem.wisc.edu/KFC_Server/index.php
11. ^https://www.molnac.unisa.it/BioTools/cocomaps/

References

Abdul, W., Aliyu, S. R., Lin, L., Sekete, M., Chen, X., Otieno, F. J., et al. (2018). Family-four aldehyde dehydrogenases play an indispensable role in the pathogenesis of Magnaporthe oryzae. Front. Plant Sci. 9:980. doi: 10.3389/fpls.2018.00980

PubMed Abstract | CrossRef Full Text | Google Scholar

Anand, P., and Stamler, J. S. (2012). Enzymatic mechanisms regulating protein S-nitrosylation: implications in health and disease. J. Mol. Med. 90, 233–244. doi: 10.1007/s00109-012-0878-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Araiza-Olivera, D., Chiquete-Felix, N., Rosas-Lemus, M., Sampedro, J. G., Peña, A., Mujica, A., et al. (2013). A glycolytic metabolon in Saccharomyces cerevisiae is stabilized by F-actin. FEBS J. 280, 3887–3905. doi: 10.1111/febs.12387

PubMed Abstract | CrossRef Full Text | Google Scholar

Assato, P. A., da Silva Jde, F., de Oliveira, H. C., Marcos, C. M., Rossi, D., et al. (2015). Functional analysis of Paracoccidioides brasiliensis 14-3-3 adhesin expressed in Saccharomyces cerevisiae. BMC Microbiol. 15:256. doi: 10.1186/s12866-015-0586-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Auger, C., Appanna, N. D., Alhasawi, A., and Appanna, V. D. (2015). Deciphering metabolic networks by blue native polyacrylamide gel electrophoresis: a functional proteomic exploration. EuPA Open Proteomics 7, 64–72. doi: 10.1016/j.euprot.2015.05.003

CrossRef Full Text | Google Scholar

Aumentado-Armstrong, T. T., Istrate, B., and Murgita, R. A. (2015). Algorithmic approaches to protein-protein interaction site prediction. Algorithms Mol. Biol. 10:7. doi: 10.1186/s13015-015-0033-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Barbosa, M. S., Báo, S. N., Andreotti, P. F., de Faria, F. P., Felipe, M. S. S., dos Santos Feitosa, L., et al. (2006). Glyceraldehyde-3-phosphate dehydrogenase of Paracoccidioides brasiliensis is a cell surface protein involved in fungal adhesion to extracellular matrix proteins and interaction with cells. Infect. Immun. 74, 382–389. doi: 10.1128/IAI.74.1.382-389.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

Bertram, G., Swoboda, R. K., Gooday, G. W., Gow, N. A., and Brown, A. J. (1996). Structure and regulation of the Candida albicans ADH1 gene encoding an immunogenic alcohol dehydrogenase. Yeast 12, 115–127. doi: 10.1002/(sici)1097-0061(199602)12:2<115::aid-yea889>3.0.co;2-e

PubMed Abstract | CrossRef Full Text | Google Scholar

Bocca, A. L., Amaral, A. C., Teixeira, M. M., Sato, P. K., Sato, P., Shikanai-Yasuda, M. A., et al. (2013). Paracoccidioidomycosis: eco-epidemiology, taxonomy and clinical and therapeutic issues. Future Microbiol. 8, 1177–1191. doi: 10.2217/fmb.13.68

PubMed Abstract | CrossRef Full Text | Google Scholar

Brandina, I., Graham, J., Lemaitre-Guillier, C., Entelis, N., Krasheninnikov, I., Sweetlove, L., et al. (2006). Enolase takes part in a macromolecular complex associated to mitochondria in yeast. Biochim. Biophys. Acta 1757, 1217–1228. doi: 10.1016/j.bbabio.2006.07.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Brummer, E., Castaneda, E., and Restrepo, A. (1993). Paracoccidioidomycosis: an update. Clin. Microbiol. Rev. 6, 89–117. doi: 10.1128/cmr.6.2.89

PubMed Abstract | CrossRef Full Text | Google Scholar

Bustos, D. M., and Iglesias, A. A. (2003). Phosphorylated non-phosphorylating glyceraldehyde-3-phosphate dehydrogenase from heterotrophic cells of wheat interacts with 14-3-3 proteins. Plant Physiol. 133, 2081–2088. doi: 10.1104/pp.103.030981

PubMed Abstract | CrossRef Full Text | Google Scholar

Carmona-Gutierrez, D., Eisenberg, T., Büttner, S., Meisinger, C., Kroemer, G., and Madeo, F. (2010). Apoptosis in yeast: triggers, pathways, subroutines. Cell Death Differ. 17, 763–773. doi: 10.1038/cdd.2009.219

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, A. C., Coleman, R. G., Smyth, K. T., Cao, Q., Soulard, P., Caffrey, D. R., et al. (2007). Structure-based maximal affinity model predicts small-molecule druggability. Nat. Biotechnol. 25, 71–75. doi: 10.1038/nbt1273

PubMed Abstract | CrossRef Full Text | Google Scholar

Colell, A., Green, D. R., and Ricci, J.-E. (2009). Novel roles for GAPDH in cell death and carcinogenesis. Cell Death Differ. 16, 1573–1581. doi: 10.1038/cdd.2009.137

PubMed Abstract | CrossRef Full Text | Google Scholar

Craig, E. A., Weissman, J. S., and Horwich, A. L. (1994). Heat shock proteins and molecular chaperones: mediators of protein conformation and turnover in the cell. Cell 78, 365–372. doi: 10.1016/0092-8674(94)90416-2

CrossRef Full Text | Google Scholar

Da Silva, S. P., Felipe, M. S. S., Pereira, M., Azevedo, M. O., and Soares, C. M. A. (1994). Phase transition and stage-specific protein synthesis in the dimorphic fungus Paracoccidioides brasiliensis. Exp. Mycol. 18, 294–299. doi: 10.1016/S0147-5975(06)80002-2

CrossRef Full Text | Google Scholar

Dai, R.-P., Yu, F.-X., Goh, S.-R., Chng, H.-W., Tan, Y.-L., Fu, J.-L., et al. (2008). Histone 2B (H2B) expression is confined to a proper NAD+/NADH redox status. J. Biol. Chem. 283, 26894–26901. doi: 10.1074/jbc.M804307200

PubMed Abstract | CrossRef Full Text | Google Scholar

Dantas, S. F. I. M., Vieira de Rezende, T. C., Bailão, A. M., Taborda, C. P., da Silva Santos, R., Pacheco de Castro, K., et al. (2009). Identification and characterization of antigenic proteins potentially expressed during the infectious process of Paracoccidioides brasiliensis. Microbes Infect. 11, 895–903. doi: 10.1016/j.micinf.2009.05.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Das, M. R., Bag, A. K., Saha, S., Ghosh, A., Dey, S. K., Das, P., et al. (2016). Molecular association of glucose-6-phosphate isomerase and pyruvate kinase M2 with glyceraldehyde-3-phosphate dehydrogenase in cancer cells. BMC Cancer 16:152. doi: 10.1186/s12885-016-2172-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Daura, X., Gademann, K., Jaun, B., Seebach, D., van Gunsteren, W. F., and Mark, A. E. (2004). Peptide folding: when simulation meets experiment. Angewandte Chem. Int. Edit. 38, 236–240. doi: 10.1039/c0cp01773f

PubMed Abstract | CrossRef Full Text | Google Scholar

de Oliveira, H. C., da Silva Jde, F., Matsumoto, M. T., Marcos, C. M., Peres da Silva, R., Moraes da Silva, R. A., et al. (2014). Alterations of protein expression in conditions of copper-deprivation for Paracoccidioides lutzii in the presence of extracellular matrix components. BMC Microbiol. 14:302. doi: 10.1186/s12866-014-0302-7

PubMed Abstract | CrossRef Full Text | Google Scholar

de Oliveira, K. M., da Silva Neto, B. R., Parente, J. A., da Silva, R. A., Quintino, G. O., Voltan, A. R., et al. (2013). Intermolecular interactions of the malate synthase of Paracoccidioides spp. BMC Microbiol. 13:107. doi: 10.1186/1471-2180-13-107

PubMed Abstract | CrossRef Full Text | Google Scholar

Doulias, P.-T., Tenopoulou, M., Greene, J. L., Raju, K., and Ischiropoulos, H. (2013). Nitric oxide regulates mitochondrial fatty acid metabolism through reversible protein S-nitrosylation. Sci. Signal. 6:rs1. doi: 10.1126/scisignal.2003252

PubMed Abstract | CrossRef Full Text | Google Scholar

Fava-Netto, C., and Raphael, A. A. (1961). Reação intradérmica com polissacarídeo do Paracaccidioides brasiliensis, na blastomicose sul-americana. Rev. Inst. Med. Trop. 4, 161–165.

Google Scholar

Felipe, M. S. S., Andrade, R. V., Arraes, F. B. M., Nicola, A. M., Maranhão, A. Q., Torres, F. A. G., et al. (2005). Transcriptional profiles of the human pathogenic fungus Paracoccidioides brasiliensis in mycelium and yeast cells. J. Biol. Chem. 280, 24706–24714. doi: 10.1074/jbc.M500625200

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferreira, E., Giménez, R., Aguilera, L., Guzmán, K., Aguilar, J., Badia, J., et al. (2013). Protein interaction studies point to new functions for Escherichia coli glyceraldehyde-3-phosphate dehydrogenase. Res. Microbiol. 164, 145–154. doi: 10.1016/j.resmic.2012.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Gavin, A.-C., Bösche, M., Krause, R., Grandi, P., Marzioch, M., Bauer, A., et al. (2002). Functional organization of the yeast proteome by systematic analysis of protein complexes. Nature 415, 141–147. doi: 10.1038/415141a

PubMed Abstract | CrossRef Full Text | Google Scholar

Glaser, F., Steinberg, D. M., Vakser, I. A., and Ben-Tal, N. (2001). Residue frequencies and pairing preferences at protein-protein interfaces. Proteins 43, 89–102. doi: 10.1002/1097-0134(20010501)43:2<89::aid-prot1021>3.0.co;2-h

CrossRef Full Text | Google Scholar

Gould, N., Doulias, P.-T., Tenopoulou, M., Raju, K., and Ischiropoulos, H. (2013). Regulation of protein function and signaling by reversible cysteine S-nitrosylation. J. Biol. Chem. 288, 26473–26479. doi: 10.1074/jbc.R113.460261

PubMed Abstract | CrossRef Full Text | Google Scholar

Graham, J. W. A., Williams, T. C. R., Morgan, M., Fernie, A. R., Ratcliffe, R. G., and Sweetlove, L. J. (2007). Glycolytic enzymes associate dynamically with mitochondria in response to respiratory demand and support substrate channeling. Plant Cell 19, 3723–3738. doi: 10.1105/tpc.107.053371

PubMed Abstract | CrossRef Full Text | Google Scholar

Grahl, N., Puttikamonkul, S., Macdonald, J. M., Gamcsik, M. P., Ngo, L. Y., Hohl, T. M., et al. (2011). In vivo hypoxia and a fungal alcohol dehydrogenase influence the pathogenesis of invasive pulmonary aspergillosis. PLoS Pathog. 7:e1002145. doi: 10.1371/journal.ppat.1002145

PubMed Abstract | CrossRef Full Text | Google Scholar

Guruharsha, K. G., Rual, J.-F., Zhai, B., Mintseris, J., Vaidya, P., Vaidya, N., et al. (2011). A protein complex network of Drosophila melanogaster. Cell 147, 690–703. doi: 10.1016/j.cell.2011.08.047

PubMed Abstract | CrossRef Full Text | Google Scholar

Hao, G., Xie, L., and Gross, S. S. (2004). Argininosuccinate synthetase is reversibly inactivated by S-nitrosylation in vitro and in vivo. J. Biol. Chem. 279, 36192–36200. doi: 10.1074/jbc.M404866200

PubMed Abstract | CrossRef Full Text | Google Scholar

Hara, M. R., Agrawal, N., Kim, S. F., Cascio, M. B., Fujimuro, M., Ozeki, Y., et al. (2005). S-nitrosylated GAPDH initiates apoptotic cell death by nuclear translocation following Siah1 binding. Nat. Cell Biol. 7, 665–674. doi: 10.1038/ncb1268

PubMed Abstract | CrossRef Full Text | Google Scholar

Hayes, S., Malacrida, B., Kiely, M., and Kiely, P. A. (2016). Studying protein-protein interactions: progress, pitfalls and solutions. Biochem. Soc. Trans. 44, 994–1004. doi: 10.1042/BST20160092

PubMed Abstract | CrossRef Full Text | Google Scholar

He, H., Lee, M.-C., Zheng, L.-L., Zheng, L., and Luo, Y. (2013). Integration of the metabolic/redox state, histone gene switching, DNA replication and S-phase progression by moonlighting metabolic enzymes. Biosci. Rep. 33:e00018. doi: 10.1042/BSR20120059

PubMed Abstract | CrossRef Full Text | Google Scholar

Hutter, J. (2011). Car–Parrinello molecular dynamics. Wiley Interdiscipl. Rev. Comput. Mol. Sci. 2, 604–612. doi: 10.1002/wcms.90

CrossRef Full Text | Google Scholar

Iannuzzi, M., Laio, A., and Parrinello, M. (2003). Efficient exploration of reactive potential energy surfaces using Car-Parrinello molecular dynamics. Phys. Rev. Lett. 90:238302. doi: 10.1103/PhysRevLett.90.238302

PubMed Abstract | CrossRef Full Text | Google Scholar

Jakob, B., Splinter, J., Conrad, S., Voss, K.-O., Zink, D., Durante, M., et al. (2011). DNA double-strand breaks in heterochromatin elicit fast repair protein recruitment, histone H2AX phosphorylation and relocation to euchromatin. Nucleic Acids Res. 39, 6489–6499. doi: 10.1093/nar/gkr230

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, H., Agarwal, S., Agarwal, S., and Pancholi, V. (2011). Surface export of GAPDH/SDH, a glycolytic enzyme, is essential for Streptococcus pyogenes virulence. mBio 2:e00068-11. doi: 10.1128/mBio.00068-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Kohr, M. J., Aponte, A. M., Sun, J., Wang, G., Murphy, E., Gucek, M., et al. (2011). Characterization of potential S-nitrosylation sites in the myocardium. Am. J. Physiol. Heart Circ. Physiol. 300, H1327–H1335. doi: 10.1152/ajpheart.00997.2010

PubMed Abstract | CrossRef Full Text | Google Scholar

Kohr, M. J., Murphy, E., and Steenbergen, C. (2014). Glyceraldehyde-3-phosphate dehydrogenase acts as a mitochondrial trans-S-nitrosylase in the heart. PLoS One 9:e111448. doi: 10.1371/journal.pone.0111448

PubMed Abstract | CrossRef Full Text | Google Scholar

Kornberg, M. D., Sen, N., Hara, M. R., Juluri, K. R., Nguyen, J. V. K., Snowman, A. M., et al. (2010). GAPDH mediates nitrosylation of nuclear proteins. Nat. Cell Biol. 12, 1094–1100. doi: 10.1038/ncb2114

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, J. K., Tabor, S., and Richardson, C. C. (2004). Proteomic analysis of thioredoxin-targeted proteins in Escherichia coli. Proc. Natl. Acad. Sci. U.S.A. 101, 3759–3764. doi: 10.1073/pnas.0308701101

PubMed Abstract | CrossRef Full Text | Google Scholar

Laine, E., and Carbone, A. (2015). Local geometry and evolutionary conservation of protein surfaces reveal the multiple recognition patches in protein-protein interactions. PLoS Comput. Biol. 11:e1004580. doi: 10.1371/journal.pcbi.1004580

PubMed Abstract | CrossRef Full Text | Google Scholar

Liaud, M. F., Lichtlé, C., Apt, K., Martin, W., and Cerff, R. (2000). Compartment-specific isoforms of TPI and GAPDH are imported into diatom mitochondria as a fusion protein: evidence in favor of a mitochondrial origin of the eukaryotic glycolytic pathway. Mol. Biol. Evol. 17, 213–223. doi: 10.1093/oxfordjournals.molbev.a026301

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, K., Qian, L., Wang, J., Li, W., Deng, X., Chen, X., et al. (2009). Two-dimensional blue native/SDS-PAGE analysis reveals heat shock protein chaperone machinery involved in hepatitis B virus production in HepG2.2.15 cells. Mol. Cell Proteomics 8, 495–505. doi: 10.1074/mcp.M800250-MCP200

PubMed Abstract | CrossRef Full Text | Google Scholar

Luong, T. T., Kim, E.-H., Bak, J. P., Nguyen, C. T., Choi, S., Briles, D. E., et al. (2015). Ethanol-induced alcohol dehydrogenase E (AdhE) potentiates pneumolysin in Streptococcus pneumoniae. Infect. Immun. 83, 108–119. doi: 10.1128/IAI.02434-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Maeda, K., Nagata, H., Kuboniwa, M., Ojima, M., Osaki, T., Minamino, N., et al. (2013). Identification and characterization of Porphyromonas gingivalis client proteins that bind to Streptococcus oralis glyceraldehyde-3-phosphate dehydrogenase. Infect. Immun. 81, 753–763. doi: 10.1128/IAI.00875-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Maheshwari, S., and Brylinski, M. (2015). Predicting protein interface residues using easily accessible on-line resources. Brief. Bioinformatics 16, 1025–1034. doi: 10.1093/bib/bbv009

PubMed Abstract | CrossRef Full Text | Google Scholar

Marchitti, S. A., Brocker, C., Stagos, D., and Vasiliou, V. (2008). Non-P450 aldehyde oxidizing enzymes: the aldehyde dehydrogenase superfamily. Expert Opin. Drug Metab. Toxicol. 4, 697–720. doi: 10.1517/17425255.4.6.697

PubMed Abstract | CrossRef Full Text | Google Scholar

Marcos, C. M., Silva Jde, F., Oliveira, H. C., Assato, P. A., Singulani Jde, L., Lopez, A. M., et al. (2016). Decreased expression of 14-3-3 in Paracoccidioides brasiliensis confirms its involvement in fungal pathogenesis. Virulence 7, 72–84. doi: 10.1080/21505594.2015.1122166

PubMed Abstract | CrossRef Full Text | Google Scholar

Menard, L., Maughan, D., and Vigoreaux, J. (2014). The structural and functional coordination of glycolytic enzymes in muscle: evidence of a metabolon? Biology 3, 623–644. doi: 10.3390/biology3030623

PubMed Abstract | CrossRef Full Text | Google Scholar

Morisse, S., Zaffagnini, M., Gao, X.-H., Lemaire, S. D., and Marchand, C. H. (2014). Insight into protein S-nitrosylation in Chlamydomonas reinhardtii. Antioxid. Redox Signal. 21, 1271–1284. doi: 10.1089/ars.2013.5632

PubMed Abstract | CrossRef Full Text | Google Scholar

Muñoz, J. F., Farrer, R. A., Desjardins, C. A., Gallo, J. E., Sykes, S., Sakthikumar, S., et al. (2016). Genome diversity, recombination, and virulence across the major lineages of paracoccidioides. mSphere 1, e213–e216. doi: 10.1128/mSphere.00213-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Murad, A. M., and Rech, E. L. (2012). NanoUPLC-MSE proteomic data assessment of soybean seeds using the Uniprot database. BMC Biotechnol 12:82. doi: 10.1186/1472-6750-12-82

PubMed Abstract | CrossRef Full Text | Google Scholar

Murad, A. M., Souza, G. H., Garcia, J. S., and Rech, E. L. (2011). Detection and expression analysis of recombinant proteins in plant-derived complex mixtures using nanoUPLC-MS(E). J. Sep. Sci. 34, 2618–2630. doi: 10.1002/jssc.201100238

PubMed Abstract | CrossRef Full Text | Google Scholar

Murphy, E., Kohr, M., Menazza, S., Nguyen, T., Evangelista, A., Sun, J., et al. (2014). Signaling by S-nitrosylation in the heart. J. Mol. Cell. Cardiol. 73, 18–25. doi: 10.1016/j.yjmcc.2014.01.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Mustafa, A. K., Gadalla, M. M., Sen, N., Kim, S., Mu, W., Gazi, S. K., et al. (2009). H2S signals through protein S-sulfhydration. Sci. Signal. 2:ra72. doi: 10.1126/scisignal.2000464

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, T. T., Stevens, M. V., Kohr, M., Steenbergen, C., Sack, M. N., and Murphy, E. (2011). Cysteine 203 of cyclophilin D is critical for cyclophilin D activation of the mitochondrial permeability transition pore. J. Biol. Chem. 286, 40184–40192. doi: 10.1074/jbc.M111.243469

PubMed Abstract | CrossRef Full Text | Google Scholar

Opperdoes, F. R., and Borst, P. (1977). Localization of nine glycolytic enzymes in a microbody-like organelle in Trypanosoma brucei: the glycosome. FEBS Lett. 80, 360–364. doi: 10.1016/0014-5793(77)80476-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Ozoe, F., Kurokawa, R., Kobayashi, Y., Jeong, H. T., Tanaka, K., Sen, K., et al. (2002). The 14-3-3 proteins Rad24 and Rad25 negatively regulate Byr2 by affecting its localization in Schizosaccharomyces pombe. Mol. Cell. Biol. 22, 7105–7119. doi: 10.1128/mcb.22.20.7105-7119.2002

PubMed Abstract | CrossRef Full Text | Google Scholar

Parente, A. F. A., de Rezende, T. C. V., de Castro, K. P., Bailão, A. M., Parente, J. A., Borges, C. L., et al. (2013). A proteomic view of the response of Paracoccidioides yeast cells to zinc deprivation. Fungal Biol. 117, 399–410. doi: 10.1016/j.funbio.2013.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Pauly, B., Lasi, M., MacKintosh, C., Morrice, N., Imhof, A., Regula, J., et al. (2007). Proteomic screen in the simple metazoan Hydra identifies 14-3-3 binding proteins implicated in cellular metabolism, cytoskeletal organisation and Ca2+ signalling. BMC Cell Biol. 8:31. doi: 10.1186/1471-2121-8-31

PubMed Abstract | CrossRef Full Text | Google Scholar

Piantadosi, C. A. (2012). Regulation of mitochondrial processes by protein S-nitrosylation. Biochim. Biophys. Acta 1820, 712–721. doi: 10.1016/j.bbagen.2011.03.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Pozuelo Rubio, M., Geraghty, K. M., Wong, B. H. C., Wood, N. T., Campbell, D. G., Morrice, N., et al. (2004). 14-3-3-affinity purification of over 200 human phosphoproteins reveals new links to regulation of cellular metabolism, proliferation and trafficking. Biochem. J. 379, 395–408. doi: 10.1042/BJ20031797

PubMed Abstract | CrossRef Full Text | Google Scholar

Prabhakar, P., Malhotra, O. P., and Kayastha, A. M. (1999). Phosphoglycerate kinase–glyceraldehyde-3-phosphate dehydrogenase interaction: reaction rate studies. Indian J. Biochem. Biophys. 36, 88–100.

Google Scholar

Pronk, S., Páll, S., Schulz, R., Larsson, P., Bjelkmar, P., Apostolov, R., et al. (2013). GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit. Bioinformatics 29, 845–854. doi: 10.1093/bioinformatics/btt055

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramzan, R., Weber, P., Linne, U., and Vogt, S. (2013). GAPDH: the missing link between glycolysis and mitochondrial oxidative phosphorylation? Biochem. Soc. Trans. 41, 1294–1297. doi: 10.1042/BST20130067

PubMed Abstract | CrossRef Full Text | Google Scholar

Rezende, T. C. V., Borges, C. L., Magalhães, A. D., de Sousa, M. V., Ricart, C. A. O., Bailão, A. M., et al. (2011). A quantitative view of the morphological phases of Paracoccidioides brasiliensis using proteomics. J. Proteomics 75, 572–587. doi: 10.1016/j.jprot.2011.08.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Rooney, P. J., and Klein, B. S. (2002). Linking fungal morphogenesis with virulence. Cell. Microbiol. 4, 127–137. doi: 10.1046/j.1462-5822.2002.00179.x

CrossRef Full Text | Google Scholar

Sawa, A., Khan, A. A., Hester, L. D., and Snyder, S. H. (1997). Glyceraldehyde-3-phosphate dehydrogenase: nuclear translocation participates in neuronal and nonneuronal cell death. Proc. Natl. Acad. Sci. U.S.A. 94, 11669–11674. doi: 10.1073/pnas.94.21.11669

PubMed Abstract | CrossRef Full Text | Google Scholar

Schägger, H., Cramer, W. A., and von Jagow, G. (1994). Analysis of molecular masses and oligomeric states of protein complexes by blue native electrophoresis and isolation of membrane protein complexes by two-dimensional native electrophoresis. Anal. Biochem. 217, 220–230. doi: 10.1006/abio.1994.1112

PubMed Abstract | CrossRef Full Text | Google Scholar

Schuppe-Koistinen, I., Moldéus, P., Bergman, T., and Cotgreave, I. A. (1994). S-thiolation of human endothelial cell glyceraldehyde-3-phosphate dehydrogenase after hydrogen peroxide treatment. Eur. J. Biochem. 221, 1033–1037. doi: 10.1111/j.1432-1033.1994.tb18821.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Seidler, N. W. (2013). GAPDH and intermediary metabolism. Adv. Exp. Med. Biol. 985, 37–59. doi: 10.1007/978-94-007-4716-6_2

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, Y., and Shi, Y. (2004). Metabolic enzymes and coenzymes in transcription–a direct link between metabolism and transcription? Trends Genet. 20, 445–452. doi: 10.1016/j.tig.2004.07.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Silva Jde, F., Vicentim, J., Oliveira, H. C., Marcos, C. M., Assato, P. A., Andreotti, P. F., et al. (2015). Influence of the Paracoccidioides brasiliensis 14-3-3 and gp43 proteins on the induction of apoptosis in A549 epithelial cells. Mem. Inst. Oswaldo Cruz 110, 476–484. doi: 10.1590/0074-02760150057

PubMed Abstract | CrossRef Full Text | Google Scholar

Singh, P., Salih, M., Leddy, J. J., and Tuana, B. S. (2004). The muscle-specific calmodulin-dependent protein kinase assembles with the glycolytic enzyme complex at the sarcoplasmic reticulum and modulates the activity of glyceraldehyde-3-phosphate dehydrogenase in a Ca2+/calmodulin-dependent manner. J. Biol. Chem. 279, 35176–35182. doi: 10.1074/jbc.M402282200

PubMed Abstract | CrossRef Full Text | Google Scholar

Soltys, B. J., and Gupta, R. S. (1999). Mitochondrial-matrix proteins at unexpected locations: are they exported? Trends Biochem. Sci. 24, 174–177. doi: 10.1016/s0968-0004(99)01390-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Someya, S., Yu, W., Hallows, W. C., Xu, J., Vann, J. M., Leeuwenburgh, C., et al. (2010). Sirt3 mediates reduction of oxidative damage and prevention of age-related hearing loss under caloric restriction. Cell 143, 802–812. doi: 10.1016/j.cell.2010.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Stephan, P., Clarke, F., and Morton, D. (1986). The indirect binding of triose-phosphate isomerase to myofibrils to form a glycolytic enzyme mini-complex. Biochim. Biophys. Acta 873, 127–135. doi: 10.1016/0167-4838(86)90198-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Sweetlove, L. J., and Fernie, A. R. (2018). The role of dynamic enzyme assemblies and substrate channelling in metabolic regulation. Nat. Commun. 9:2136. doi: 10.1038/s41467-018-04543-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Szklarczyk, D., Franceschini, A., Kuhn, M., Simonovic, M., Roth, A., Minguez, P., et al. (2011). The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored. Nucleic Acids Res. 39, D561–D568. doi: 10.1093/nar/gkq973

PubMed Abstract | CrossRef Full Text | Google Scholar

Tarze, A., Deniaud, A., Le Bras, M., Maillier, E., Molle, D., Larochette, N., et al. (2007). GAPDH, a novel regulator of the pro-apoptotic mitochondrial membrane permeabilization. Oncogene 26, 2606–2620. doi: 10.1038/sj.onc.1210074

PubMed Abstract | CrossRef Full Text | Google Scholar

Tisdale, E. J. (2001). Glyceraldehyde-3-phosphate dehydrogenase is required for vesicular transport in the early secretory pathway. J. Biol. Chem. 276, 2480–2486. doi: 10.1074/jbc.M007567200

PubMed Abstract | CrossRef Full Text | Google Scholar

Tomokuni, Y., Goryo, K., Katsura, A., Torii, S., Yasumoto, K., Kemnitz, K., et al. (2010). Loose interaction between glyceraldehyde-3-phosphate dehydrogenase and phosphoglycerate kinase revealed by fluorescence resonance energy transfer-fluorescence lifetime imaging microscopy in living cells. FEBS J. 277, 1310–1318. doi: 10.1111/j.1742-4658.2010.07561.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Tovchigrechko, A., and Vakser, I. A. (2006). GRAMM-X public web server for protein-protein docking. Nucleic Acids Res. 34, W310–W314. doi: 10.1093/nar/gkl206

PubMed Abstract | CrossRef Full Text | Google Scholar

Travassos, L. R., Taborda, C. P., and Colombo, A. L. (2008). Treatment options for paracoccidioidomycosis and new strategies investigated. Expert Rev. Anti Infect. Ther. 6, 251–262. doi: 10.1586/14787210.6.2.251

PubMed Abstract | CrossRef Full Text | Google Scholar

Tribble, G. D., Mao, S., James, C. E., and Lamont, R. J. (2006). A Porphyromonas gingivalis haloacid dehalogenase family phosphatase interacts with human phosphoproteins and is important for invasion. Proc. Natl. Acad. Sci. U.S.A. 103, 11027–11032. doi: 10.1073/pnas.0509813103

PubMed Abstract | CrossRef Full Text | Google Scholar

Tristan, C., Shahani, N., Sedlak, T. W., and Sawa, A. (2011). The diverse functions of GAPDH: views from different subcellular compartments. Cell Signal. 23, 317–323. doi: 10.1016/j.cellsig.2010.08.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Turissini, D. A., Gomez, O. M., Teixeira, M. M., McEwen, J. G., and Matute, D. R. (2017). Species boundaries in the human pathogen Paracoccidioides. Fungal Genet. Biol. 106, 9–25. doi: 10.1016/j.fgb.2017.05.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Vangone, A., Spinelli, R., Scarano, V., Cavallo, L., and Oliva, R. (2011). COCOMAPS: a web application to analyze and visualize contacts at the interface of biomolecular complexes. Bioinformatics 27, 2915–2916. doi: 10.1093/bioinformatics/btr484

PubMed Abstract | CrossRef Full Text | Google Scholar

Vizcaíno, J. A., Csordas, A., del-Toro, N., Dianes, J. A., Griss, J., Lavidas, I., et al. (2016). 2016 update of the PRIDE database and its related tools. Nucleic Acids Res. 44, D447–D456. doi: 10.1093/nar/gkv1145

PubMed Abstract | CrossRef Full Text | Google Scholar

Wan, C., Borgeson, B., Phanse, S., Tu, F., Drew, K., Clark, G., et al. (2015). Panorama of ancient metazoan macromolecular complexes. Nature 525, 339–344. doi: 10.1038/nature14877

PubMed Abstract | CrossRef Full Text | Google Scholar

Weber, S. S., Parente, A. F. A., Borges, C. L., Parente, J. A., Bailão, A. M., and de Almeida Soares, C. M. (2012). Analysis of the secretomes of Paracoccidioides mycelia and yeast cells. PLoS One 7:e52470. doi: 10.1371/journal.pone.0052470

PubMed Abstract | CrossRef Full Text | Google Scholar

Wong, J. H., Cai, N., Balmer, Y., Tanaka, C. K., Vensel, W. H., Hurkman, W. J., et al. (2004). Thioredoxin targets of developing wheat seeds identified by complementary proteomic approaches. Phytochemistry 65, 1629–1640. doi: 10.1016/j.phytochem.2004.05.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, X., Song, Y., Li, Y., Chang, J., Zhang, H., and An, L. (2010). The tandem affinity purification method: an efficient system for protein complex purification and protein interaction identification. Protein Expr. Purif. 72, 149–156. doi: 10.1016/j.pep.2010.04.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Xue, Y., Liu, Z., Gao, X., Jin, C., Wen, L., Yao, X., et al. (2010). GPS-SNO: computational prediction of protein S-Nitrosylation sites with a modified GPS algorithm. PLoS One 5:e11290. doi: 10.1371/journal.pone.0011290

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J., Yan, R., Roy, A., Xu, D., Poisson, J., and Zhang, Y. (2015). The I-TASSER suite: protein structure and function prediction. Nat. Methods 12, 7–8. doi: 10.1038/nmeth.3213

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, X., and Mitchell, J. C. (2011). KFC2: a knowledge-based hot spot prediction method based on interface solvation, atomic density, and plasticity features. Proteins 79, 2671–2683. doi: 10.1002/prot.23094

PubMed Abstract | CrossRef Full Text | Google Scholar