Insight Into the Molecular Mechanisms Underpinning the Mycoremediation of Multiple Metals by Proteomic Technique

We investigated the fungus Aspergillus fumigatus PD-18 responses when subjected to the multimetal combination (Total Cr, Cd2+, Cu2+, Ni2+, Pb2+, and Zn2+) in synthetic composite media. To understand how multimetal stress impacts fungal cells at the molecular level, the cellular response of A. fumigatus PD-18 to 30 mg/L multimetal stress (5 mg/L of each heavy metal) was determined by proteomics. The comparative fungal proteomics displayed the remarkable inherent intracellular and extracellular mechanism of metal resistance and tolerance potential of A. fumigatus PD-18. This study reported 2,238 proteins of which 434 proteins were exclusively expressed in multimetal extracts. The most predominant functional class expressed was for cellular processing and signaling. The type of proteins and the number of proteins that were upregulated due to various stress tolerance mechanisms were post-translational modification, protein turnover, and chaperones (42); translation, ribosomal structure, and biogenesis (60); and intracellular trafficking, secretion, and vesicular transport (18). In addition, free radical scavenging antioxidant proteins, such as superoxide dismutase, were upregulated upto 3.45-fold and transporter systems, such as protein transport (SEC31), upto 3.31-fold to combat the oxidative stress caused by the multiple metals. Also, protein–protein interaction network analysis revealed that cytochrome c oxidase and 60S ribosomal protein played key roles to detoxify the multimetal. To the best of our knowledge, this study of A. fumigatus PD-18 provides valuable insights toward the growing research in comprehending the metal microbe interactions in the presence of multimetal. This will facilitate in development of novel molecular markers for contaminant bioremediation.


INTRODUCTION
Heavy metals, such as copper (Cu), chromium (Cr), cadmium (Cd), zinc (Zn), lead (Pb), and nickel (Ni), occur in the water bodies, such as river water, and drains, in developing countries that are above the permissible mandates as prescribed by the Food and Agriculture Organization (Bhattacharya et al., 2015). These heavy metals are released from numerous smalland medium-scale enterprises (pesticide, textile, electroplating, fertilizer, batteries, etc.) due to lack of or improper sewage treatment plant systems (Bhattacharya et al., 2015;Bhardwaj et al., 2017;Gola et al., 2020). These hazardous contaminants are discharged into various water bodies via irrigation during agricultural activities. This results in bioaccumulation of these contaminants that enter the food chain and cause detrimental health ailments, such as cancer in human beings (Kou et al., 2018;Li et al., 2019;Abed et al., 2020;Chauhan et al., 2021;Irawati et al., 2021). Thus, these heavy metals are the environmental priority contaminants threatening the environment and therefore must be remediated before discharge into the environment.
The conventional physico-chemical methods, such as chemical precipitation, ion exchange, adsorption, membrane filtration, coagulation-flocculation, and flotation, are usually utilized to remediate these harmful contaminants. These techniques have the disadvantages of being expensive, having low selectivity, production of additional sludge, and further treatment is required for better results (Dey et al., 2021;Zamora-Ledezma et al., 2021).
On the other hand, living and actively growing microbial cells can be a lucrative option for bioremediation (Malik, 2004). Of the bioremediation techniques, mycoremediation has shown to be a promising technology having the potential to ameliorate these hazardous chemicals Sabuda et al., 2020). The superiority of fungi over singlecelled microbes, such as bacteria, to remediate these recalcitrant heavy metals is well-documented (Deshmukh et al., 2016). In addition, the fungi are omnipresent, multifarious, and have a wider arsenal to acclimatize to environmental limitations, such as immoderations of temperature, extremes of pH, higher metal concentrations, and low nutrient accessibility, due to their morphological diversity (Anand et al., 2006). Besides, fungal mycelia have enhanced enzymatic and mechanical contact with the pollutant due to a greater cell to the surface ratio (Sagar and Singh, 2011).
Further, in microbial cells, such as fungus, heavy metals are key components in the number of catalytic and structural proteins that are integral to biochemical processes. The outcome of these processes differs depending on the type of metal involved and its concentration inside the cell (Gadd, 1994). Particularly, the species of Aspergillus have a high metal uptake capacity for metals, such as Cu (Dusengemungu et al., 2020). Also, filamentous fungi develop signature metabolic pathways which are species-specific to survival in the harsh environment of heavy metals and other contaminants that are utilized as nutrients and energy sources (Kuhn and Käufer, 2003). Accordingly, proteins expressed in cells under such diverse conditions and at different times are dissimilar (Boopathy, 2000;Keller, 2015;Wisecaver and Rokas, 2015).
Also, there is evidence that fungal resistance toward one element does not necessarily infer resistance to another element even though the elements possess similar valency charges (Høiland, 1995). For example, the toxicity of Pb toward the microbe is less compared to other toxic metals, such as Cd, As, and Hg (Jiang et al., 2020). Further, in fungi, the metallothioneins formation is predominantly induced by the heavy metal Cu (Jaeckel et al., 2005). Thus, understanding the altered heavy metal uptake in the fungus in presence of diverse environmental conditions needs to be evaluated.
To gain such mechanistic insight, high throughput techniques, such as proteomic analysis by LC-MS/MS, can be used to identify and characterize proteins involved in the multimetal resistance mechanism in a filamentous fungus (Pandey and Mann, 2000;Kraut et al., 2009;Zhang et al., 2010). Moreover, proteomics facilitates the development of new and important protein biomarkers that specify and monitor metal contamination in the environments as protein type and also their copy numbers are estimated by the translational regulation (Ohno et al., 2014).
Several researchers have studied the change in proteomes under various metal conditions. Selamoglu et al. (2020) have assessed the change in the differential proteome expression of the prokaryotic organism Rhodobacter capsulatus in presence of 5 µM Cu by nano-LC-MS/MS. About 75 proteins were significantly regulated. Most of the proteins present were responsible for maintaining Cu homeostasis. Further, Lotlikar et al. (2020) enumerated the variable expression of proteins in a eukaryotic fungus Penicillium chrysogenum in the presence of 100 mg/L and 500 mg/L Cu. Several key proteins related to genetic information, carbohydrate metabolism, glycan biosynthesis and metabolism, amino acid metabolism, and energy metabolism were expressed. Cherrad et al. (2012) highlighted the overaccumulation of proteins of the oxidoreductase family when exposed to Cd, Cu, and Ni but not when exposed to Zn. Thus, the secretion of proteins in fungus is extremely dynamic and its production depends on different environmental triggers.
However, despite these growing proteomic studies dealing with the bioremediation of single heavy metal, there is a dearth of information on the modulated proteins triggered by the cumulative toxicity of hexametals, namely, Cd, Cr, Cu, Ni, Pb, and Zn, specifically in a filamentous fungus. Interestingly, Aspergillus fumigatus has high adaptability when subjected to an altered environment (Bakti et al., 2018). Also, as established in our previous works (Dey et al., 2016Bhattacharya et al., 2020), A. fumigatus PD-18 is a filamentous fungus capable of removal of 30 mg/L multimetal exceptionally well. Furthermore, there are structural and functional similarities between the numerous genes of lower eukaryotes, such as fungi and mammals (Bae and Chen, 2004).
Thus, this fungus A. fumigatus PD-18 would be a good eukaryotic model for helping us understand how cells adopt various cellular strategies and lay a foundation study to decipher the enzymes produced in the presence of simultaneous effects of multimetal cocktail on interaction with a fungus for scaleup process.

Chemicals and Reagents
The stock solutions (10 g/L) of different individual metals were made by dissolving their respective salts, viz. K 2 Cr 2 O 7 , Cd(NO 3 ) 2 , Ni(NO 3 ) 2 , Cu(NO 3 ) 2 , Zn(NO 3 ) 2 , and Pb(CH 3 COO) 2 in double-distilled water and were diluted to the concentrations that were required for the experiments. For preparing the reagents and calibration standards, deionized ultrapure water (RIONS Ultra 370 series) was used. Rest all other chemicals utilized were of analytical grade and were obtained from Merck, Sigma, and Qualigens.

Microorganism and Culture Media Composition
The fungal strain used was A. fumigatus PD-18 which was isolated from the polluted banks of the river Yamuna, New Delhi, India, and characterized with accession number KX365202 after depositing the sequence to the Genbank (NCBI) (Dey et al., 2016).

Methodology
The methodology adopted in this study was to estimate the differential expression of proteins in A. fumigatus PD18 brought about by 30 mg/L multimetal (MM) viz. Five milligram per liter of each of the individual Cd, Total Cr, Cu, Ni, Pb, and Zn amended in the composite media in addition to 1% of glucose. The concentration of 5 mg/L of each heavy metal was chosen to take into consideration the permissible mandates for heavy metal occurrence in water that can be utilized for irrigation according to Food and Agriculture Organization (FAO). The prescribed limit for each heavy metal is as follows: Cd, 0.01; Total Cr, 0.1; Cu, 0.2; Ni, 0.2; Pb, 5.0; Zn, 2.0. Also, the typical concentration of heavy metals occurring in a mixture in the Yamuna river was considered (Bhattacharya et al., 2015). The biotic control included composite media added with only 1% of glucose.
Studies were conducted in a series of Erlenmeyer flasks (250 mL) comprising 100 mL of composite growth media. One milliliter of spore suspension (having a concentration of 10 7 spores eluted with sterile distilled water containing 0.01% Tween 80) was inoculated in the flasks and incubated at 30 • C and 150 rpm agitation for 72 h to ensure complete uptake of metal ions. The flasks were withdrawn after the late log phase incubation period. Analysis was done in three technical replicates for biotic control and multimetal treated samples separately.

Cell Lysis and Protein Extraction
The fungal pellets were separated from the media by inversion and followed by centrifugation at 4,000 g for 10 min at 4 • C. The mycelia were further snap-freezed in liquid nitrogen and then lyophilized and stored as a dry powder at -20 • C before the subsequent steps. The total protein from the dried powdered fungal cells was extracted in a buffer of composition (8M Urea and 2M Thiourea) as per the modified protocol (Rughöft et al., 2020). The supernatant was collected by centrifugation at 14,000 rpm, 4 • C for 10 min, and precipitated overnight in five volumes of ice-cold acetone and the pellet was stored at -20 • C for further use. The concentration of protein to be measured was estimated according to the Bradford assay of protein quantification (Bradford, 1976). In brief, a 2-mg/mL concentration stock solution of Bovine Serum Albumin (BSA) was prepared in water. BSA standard curve was prepared using the following concentration range (0-2.00 µg/µL) followed by the addition of Coomassie Brilliant Blue dye reagent in 95% ethanol and 100 mL 85% (w/v) phosphoric acid and incubated for 5 min. The absorbance was measured at 595 nm.

Protein Separation by 1D-SDS PAGE and in Gel Tryptic Digestion of the Proteins
The separation of the proteins was carried out by onedimensional SDS-PAGE and the separated protein bands were visualized by Coomassie brilliant blue staining. The individual gel bands were cut into small pieces and washed with 10 mM ammonium bicarbonate followed by reduction with 10 mM 1,4-Dithiothreitol and 100 mM 2-Iodoacetamide. The proteins lysates were subjected to tryptic digestion and incubated overnight and kept at 37 • C. The peptides were extracted from the gel pieces in 5 mM ammonium bicarbonate and further purified by C-18 ZipTip columns according to the protocol (Haange et al., 2012). Thereafter, the peptide lysates were dried in a vacuum centrifuge and stored at -20 • C until further analysis.

Mass Spectrometric Analysis
The dried peptides lysates were reconstituted in 0.1% (v/v) formic acid and loaded and separated by reversed-phase chromatography before being measured on mass spectrometer UPLC-coupled (Waters) LTQ Orbitrap Velos MS/MS (Thermo Fisher Scientific).

Data Analysis (Identification and Characterization of the Proteins)
The raw MS-MS ion spectra data from the instrument were processed using the Proteome Discoverer software (v1.0 build 43, Thermo Fisher Scientific) (Orsburn, 2021). The fixed and variable modifications (Carbamidomethylation at cysteines was given as fixed and oxidation of methionines as variable modifications) were taken into consideration. Only ranked one peptide hits with ≥1 high confidence and less than 1% false discovery rate (FDR) were considered as identified and taken up for further analysis. The intensity data were further processed by logarithmic transformation and normalization for the determination of changes in the abundance of proteins (at least 1.5-fold upregulated or downregulated). Only the values of p < 0.05 were considered as statistically significant. Triplicate gels were run for biotic control and MM samples. For functional annotation of the proteins, the generated protein lists were analyzed by prophane. 1 The creation of the VENN diagrams was done by the online tool VENNY. 2 The heatmaps and volcano plots were constructed using the heatmap and volcano programs from R tool version 3.2.0. The WoLF PSORT tool 3 was used for the prediction of the subcellular localization of the proteins. The proteinprotein interactions between the differentially expressed proteins were studied by the STRING (Search Tool for the Retrieval of Interacting Genes/Proteins) database. 4

RESULTS
Fungi usually acclimate by making dynamic changes in the cell structure and composition when exposed to heavy metals. In our previous study on the toxicity of individual metals, that is, Cd, Cr, Cu, Ni, Pb, and Zn (500 mg/L) in A. fumigatus PD-18, we determined the highest tolerance index for each metal in solid media. We calculated the cube root growth (k) constant of A. fumigatus PD-18 when subjected to 30 mg/L multimetal in liquid composite media. This fungal strain had exceptional multimetal removal ability after ∼72 h. The growth profile of the fungus was altered when multimetal was added to the composite media and substantial variations were observed when compared with individual metals (Dey et al., 2016).
Further, the morphological changes in this fungus in response to the 30 mg/L multimetal were determined by scanning electron microscope, the localization of the heavy metals inside the fungal cell by transmission electron microscopy, and chelation of the heavy metals with the functional groups occurring in the fungus by Fourier-transform infrared spectroscopy .

Expression of Proteins in Biotic Control and Multimetal Extracts of Aspergillus fumigatus PD-18
From the six sample preparations (conditions, including one biological control and one multimetal exposure and their three technical replicates), 434 proteins were uniquely present in multimetal extracts and 400 proteins in biotic control of A. fumigatus PD-18.
From Figure 2A, the proteins that were expressed in higher abundance in A. fumigatus when treated with 30 mg/L multimetal belonged to the main metabolic functional classes of RNA processing & modification, post-translational modification & protein turnover & chaperones, nuclear structure, intracellular trafficking & secretion & vesicular transport, translation & ribosomal structure & biogenesis, and nucleotide transport & metabolism. Further, hypothetical proteins with general function prediction were also expressed which could play role in tackling the multimetal stress. The afflicted metabolic functions classes where production of proteins decreased belonged to amino acid transport & metabolism, carbohydrate transport & metabolism, lipid transport & metabolism, and secondary metabolites biosynthesis & transport and catabolism. This proteomic study displayed proteins in A. fumigatus PD-18 that showed a significant fold increase and decrease of differentially expressed proteins during growth in 30 mg/L multimetal supplemented composite media as depicted in the volcano plot in Figure 2B. Among the differentially expressed proteins identified, 399 proteins were differentially upregulated while 266 proteins were differentially downregulated. The highly upregulated protein identified hydroxymethylglutaryl-CoA synthase upregulated upto 5.2-fold times. Other important proteins with general function prediction were upregulation of serine/threonine-protein phosphatase upto 3.42-fold times and G-protein beta subunit SfaD upto 2.98-fold. The highly upregulated and downregulated proteins by fold changes and their locations inside the fungal cell are depicted in Table 1. Furthermore, it was found that the majority of the upregulated proteins were located in the cytoplasm (18), mitochondria (12) followed by extracellular (6), nuclear (5), plasma membrane (2), peroxisome (2), nuclear and cytoplasmic (1), and endoplasmic reticulum lumen (1).

Enriched Protein Network and Pathways in Aspergillus fumigatus PD-18
The interaction network as constructed by the STRING database of the highly upregulated proteins in A. fumigatus is depicted in Figure 3. The strength of the association is depicted by the thickness of the line. The number of nodes is 263 that represents the proteins, the number of edges is 467 that represents associations, the average node degree is 3.55, and the average local clustering coefficient is 0.412. The STRING analysis identified three clusters of protein interactions. The first major cluster of proteins (green) was involved in ribosome and preinitiation factors (60S and 40S ribosomal proteins). The predicted partners functionally associated with 60S and 40S ribosomal proteins include translation elongation factor, fibrillarin of class translation, ribosomal structure, and biogenesis. The second cluster   The third cluster proteins (blue) are proteasome regulatory particles, Hsp70 chaperone, protein transport protein (sec31), and V-type proton ATPase. The significant junction proteins in protein-protein interaction identified were cytochrome c, oxidase, and septin suggest that energy production and conversion and cell cycle control, cell division, and chromosome partitioning that play a pivotal role in the resistance toward multimetal exposure. Figure 4 summarizes the salient mechanism of the multimetal detoxification process by A. fumigatus PD-18.

DISCUSSION
In this study, we identified the chief classes of proteins that are crucial for the resistance and tolerance of A. fumigatus PD-18 for multimetals, including intracellular and extracellular mechanisms of metal uptake. These two processes are further elaborated.

Intracellular and Extracellular Mechanism of Metal Uptake by Fungus
The mechanism of the fungal resistance and tolerance could be attributed to its occurrence at the metal-contaminated site (Baker, 1987;Gadd and White, 1993). The fungal detoxification process involves strategies, such as intracellular bioaccumulation, extracellular precipitation, biotransformation, biomineralization, and biosorption, that involve several signaling pathways . In addition, the reduction of heavy metals from the cell comprises of activation of different metabolic processes in the cell that are stimulated by heavy metals (Goyal et al., 2003).
Fungi largely counter heavy metals in two ways. The first way involves averting the metal uptake and its passage inside the fungal cell. This happens chiefly by a decrease in metal uptake or increased efflux of metals, metal biosorption to the impermeable cell walls by metal binding polysaccharides, peptides, extracellular formation of complexes, and the release of organic acids that chelates heavy metals outside the cell. Thus, the extracellular mechanism operates in the cell to counter heavy metals by circumventing the entry of heavy metals. Particularly, the secondary metabolites, such as citric acid, oxalic acids, succinate, and fumarate, which are low molecular weight compounds (<900 daltons) are secreted by fungal species namely A. niger specifically in response to heavy metals exposure and they bind the heavy metals extracellularly (Kolen, 2013). Also, the secondary metabolite oxalic acid is produced as an intermediate compound in the biochemical tricarboxylic acid cycle (TCA) (Munir et al., 2001). Here, we found the upregulation of the key enzymes that are involved in the TCA cycle viz. phosphoglycerate kinase (2.91-fold), glyceraldehyde-3-phosphate dehydrogenase (2.38-fold), enolase (1.8-fold), and pyruvate kinase (1.86-fold) which mediates the enhanced secretion of oxalic acids.
In the second way, the fungus subsists the high concentration of metals inside the cell by tolerance after the process of detoxification via metal chelation by synthesizing ligands, such as metallothioneins and phytochelatins, that bind heavy metals intracellularly or by compartmentalization of heavy metals within the cell organelles of vacuoles by polyphosphates. The three main classes of intracellular peptides binding metal ions are phytochelatins (PCs), metallothioneins (MTs), and glutathione (GSH). MTs are low molecular weight cysteine-rich metalbinding proteins that have high affinity toward both the essential metal ions, such as Cu and Zn, non-essential metal ions, such as Cd, Hg, and Ag, and also have large metal-binding capacities (Reddy et al., 2014). Further, MTs chelate heavy metals by forming thiolate bonds with the heavy metals. Glutathione S transferases (GSTs) are enzymes that metabolize heavy metals and other contaminants by catalyzing the binding of glutathione to non-polar compounds comprising of electrophilic nitrogen, carbon, and sulfur atom (Morel et al., 2009). Usually, the metals Cd, Cu, Pb, and Zn are removed via glutathione (GSH)-mediated sequestration. However, in this study, there was no evidence of the production of proteins glutathione, metallothioneins, and phytochelatins despite the presence of these heavy metals Cd, Cu, Pb, and Zn in the multimetal mixture. The reason for this phenomenon could be attributed to the dynamics of individual heavy metals when present in a mixture. As the expression of GST is related to the type of heavy metal, its concentration, and the extent of treatment time of the heavy metal (Shen et al., 2015).
Further, it is reported that heavy metals induce oxidative damage to the cell membranes of fungi by the generation of reactive oxygen species (ROS). These ROS are detoxified by the production of antioxidants that are components from the thioredoxin system, such as peroxiredoxins, NADPH dehydrogenases, catalase, superoxide dismutase, and peroxidase, that enables the fungus to confront the reactive-oxygen species that accumulate in the cell on exposure to the metals (Zhang et al., 2015). Thus, in principle, intracellular mechanisms decrease the metal load in the cytosol (Sandau et al., 1996;Gadd, 2000Gadd, , 2007Gadd, , 2010. In this study, we found 3.45-fold upregulation of the antioxidant protein Cu-Zn superoxide dismutase. Other important functional groups were detected that expressed amino acid metabolism, lipid metabolism, energy metabolism, and also the proteins involved in signal transduction, transcription, translation, or DNA repair. In general, the upregulated proteins are stimulated to display the fungal resistance against the contaminants' stress, while the downregulated proteins are suppressed by the action of the pollutants' toxicity. The highly upregulated protein hydroxymethylglutaryl-CoA synthase with 5.19-fold upregulation is responsible for the production of secondary metabolite carotenoid from the precursor molecule of acetyl-CoA when stimulated by the heavy metal stress (Bhosale, 2004). The protein serine/threonine phosphatase with 3.42-fold upregulation is responsible for maintaining the conformation of cell organelles and proteasomes (Dias et al., 2019).
The other proteins, such as G-protein beta subunit SfaD, with 2.98-fold upregulation depicted the role of G-proteincoupled receptors (GPCRs) in heavy metal bioremediation. These are the largest transmembrane receptors that aid in communicating the extracellular signals, such as stresses of heavy metals into the intracellular sites. G-protein-coupled receptors (GPCRs) regulate the important effector molecules, such as adenylate cyclase and phospholipase C, and regulate the function of kinase and ion channel by producing secondary messengers, such as cAMP, thereby inducing signaling cascades (El-Defrawy and Hesham, 2020).
From Figures 2A, discussion of the highly regulated selected KOG classes are as follows:

Post-translational Modification, Protein Turnover, Chaperones
Protein homeostasis is crucial for the cell proliferation and viability of all organisms. Further, cellular signaling is greatly affected by the protein homeostasis under different physiological conditions and environmental stresses, such as heavy metals, and therefore they can be suitable biomarkers. Molecular chaperones aid in the delivery of metal ions to the cell organelles and metalloproteins (Lotlikar and Damare, 2018). The two vital regulators of molecular chaperones in the proteostasis network are heat shock transcription factor Hsf1 and heat shock protein Hsp90. Under a stressful environment, heat shock protein enables the folding of newly synthesized proteins and helps in the degradation of damaged/misfolded proteins with the help of the ubiquitin-proteasome system (Hossain et al., 2020). Heat shock proteins bind to the denatured proteins, compelling them to refold into their native conformation and regain their original structure (Feng et al., 2018). In this study, we observed 42 proteins upregulated that belonged to posttranslational modification, protein turnover, and chaperones. There was upto 4-fold upregulation of proteasome regulatory particle subunit Rpt3 (KOG0727), upto 2-fold upregulation of Hsp70 chaperone (HscA) (KOG0101), and upto 2.5-fold upregulation of protein geranylgeranyltransferase (KOG1439). The enzyme geranylgeranyltransferase I (GGTase I) aids in the catalysis of the post-translational transfer of lipophilic diterpenoid geranylgeranyl molecule to the cysteine residue of proteins with the termination at CaaX motif (Rho1p and Cdc42p). This alteration helps in the membrane localization of the protein and thereby rendering it biologically active. Rho1p is a regulatory subunit of 1,3-β-D-glucan synthesis and contributes to the cell wall synthesis in fungi which is vital for cell viability under stressful condition of excess metals (Singh et al., 2005).

Translation, Ribosomal Structure, and Biogenesis
Different proteins related to protein translation under multimetal stress were overexpressed. Here, we found 60 proteins of translation, ribosomal structure and biogenesis upregulated. There was upto 3.0-fold increase in glutamyl-tRNA synthetase (KOG1147), upto 2.5-fold increase in eukaryotic translation initiation factor 3 subunit B (KOG2314), upto 2.0-fold increase in 60S ribosomal protein L23 (KOG1751), upto 2.3-fold increase in mitochondrial translation initiation factor IF 2 (KOG1144), and 2.0-fold increase 40S ribosomal protein S10b (KOG3344). Similar elements of protein synthesis, such as translation initiation factor 5A, elongation factor 2, 40S and 60S ribosomal proteins, ATPdependent RNA helicase, and aspartyl-tRNA synthetase, were overexpressed in Phanerochaete chrysosporium under Cu stress as a result of the need for production of new proteins or renewal of the damaged proteins (Okay et al., 2020).

Intracellular Trafficking and Secretion and Vesicular Transport
Proteins such as ion transporters and other solutes are crucial for processes such as detoxification, cell nutrition, cell signaling, cellular homeostasis, and resistance toward metal stress. These polytopic transmembrane proteins are translated altogether and folded in the endoplasmic reticulum (ER) of the eukaryotic cells that are later ultimately arranged to their respective membrane location through vesicular secretion. During any physiological or stressful environment, transporters undertake several regulated turnovers. Thus, in the process, transporters briefly interact dynamically with multiple proteins (Dimou et al., 2021). In this study, we found 18 proteins of intracellular trafficking and secretion and vesicular transport upregulated. The levels of the proteins were expressed in higher amounts (SEC31) (KOG0307) by 3.3-fold, endosomal cargo receptor (P24) (KOG1692) by 2.9fold, and mitochondrial inner membrane translocase (KOG2580) by 2.4-fold in A. fumigatus under the effect of multimetal stress. The vesicle (Ves) are tissues composed of a lipid bilayer whose size varies ∼nanometers to micrometers. The Ves structures fuse with the plasma membrane of the cell and eject the trapped materials either inside or outside of the cytoplasm. There are three types of intracellular Ves viz. protein complex I (COPI)-coated Ves, protein complex II (COPII)-coated Ves, and BAR-domain protein Ves. These Ves proteins aid in physiological processes, such as the exchange of proteins and RNA intercellularly (Jiang et al., 2020).

Energy Production and Conversion
Heavy metals, such as Cd, Cu, Ni, and Zn, function as cofactors in bacteria and fungi. However, excess amounts of these metals are toxic to these cells and also produce reactive oxygen species (Liu et al., 2017). The need for metabolic energy in the fungal increases during abiotic stress, such as exposure to excess heavy metals. Thus, ATPases are responsible for the biochemical and physiological processes by the production of energy. Heavy metal ATPases (HMAs) or P-type ATPases can be categorized into three groups namely, A, B, and C. Further, P-type ATPases are utilized by numerous organisms to facilitate the transport of cations viz. Na + , K + , and Ca 2+ . To eliminate these excess metals, fungal HMA Saccharomyces cerevisiae CCC2 (Group A) localizes metals to metal-containing proteins, for example, in the case of copper metal, copper-containing protein FET3 in trans-golgi compartment transports metals to the cell membrane via efflux pumps, such as cadmium efflux pumps, encoded by fungal HMA Saccharomyces cerevisiae PCA1 (Group B and Group C) (Adle et al., 2007). Saitoh et al. (2009) studied the CCC2-type HMA gene that targets copper-containing proteins from the fungus Cochliobolus heterostrophus by cloning. There was upto 0.78-fold upregulation in the production of V-type ATPase and upto 1fold upregulation in the production of mitochondrial ATPase subunit ATP4. This gene has other multifarious roles, such as in the formation of dark brown colored melanin pigment located in fungal cell walls, that also sequester metals (Chang et al., 2019).

Amino Acid Transport and Metabolism
The nitrogen cycle is essential for nitrogen assimilation and transformation and also for stress tolerance. Heavy metals impact the enzymes that play important role in nitrogen metabolism (Khouja et al., 2014). There was upto 4-fold upregulation in the production of aminotransferase (KOG1549) in response to multimetal. A similar response was observed by Okay et al. (2020), where the production of enzyme aspartate aminotransferase was enhanced in the fungus Phanerochaete chrysosporium to tackle the Cu stress. This enzyme has a possible role in the renewal of the mitochondrial NAD/NADH imbalance.

Analysis of Biological Pathways and Protein-Protein Interactions
In addition to these mechanisms, there are contributions from other regulatory systems, such as cross-talks in various pathways, interconnection amongst these different pathways, and regulation of different genes. These regulatory networks of the microbial proteins are intricate and play crucial roles in resistance to metal contaminants by modifying the series of specific functional proteins/non-proteins and altering the different metabolic enzymes at the cellular level. The metabolic processes related to the detoxification of contaminants are usually regulated by the complete set of proteins and their networks instead of a single enzyme (Zhao and Poh, 2008;Feng et al., 2018).
The STRING analysis displayed protein-protein interaction networks that are related to the resistance and tolerance mechanism A. fumigatus PD-18 for multimetals. Protein interactions in the first cluster (green) are involved in ribosome and preinitiation factors (60S and 40S ribosomal proteins). Ribosomal proteins form the protein part of ribosomes and participate in protein synthesis in cells in conjunction with rRNA. Thus, the increased expression of the large subunit of ribosome renders resistance against abiotic stresses, such as heavy metals, radiation, cold, and salt (Liu et al., 2014).
The important hub proteins expressed were cytochromec oxidase that is mitochondrial proteins and catalyst in the electron transport chain and is responsible for the transport of heavy metals particularly copper (Dias et al., 2019). The second cluster (blue) had hub protein ATP synthase. The important interconnecting protein in the network is the septin protein of KOG class cell cycle control, cell division, and chromosome partitioning, and is involved in vesicular trafficking and countering the apoptotic cell death initiated by the toxicity of heavy metals. This energy-intensive process is mediated by ATP synthase to maintain the cellular structure and function under the lethal environment of heavy metals. Excess to or above the permissible limit of heavy metals exposure can substantially delimit the growth of organisms. Here, more growth is initiated in the fungal cell as a result of the enhanced activity of ATP synthase (Yıldırım et al., 2011). This was evident in our study by the upregulation of tubulin by 1.1-fold and actin proteins by 2.2-fold which are responsible for cellular division and growth (Marks et al., 1986). This enhanced growth in the fungus also corroborated with our previous study where there was an increased dry weight of the fungal biomass under 30 mg/L multimetal (Dey et al., 2016). These results showed that the proteins in this network played important functions in cell functioning under heavy metals stress.
The mechanism of hexametal uptake by A. fumigatus PD-18 has been summarized in Figure 4.

CONCLUSION
In this study, we found the fungus A. fumigatus PD-18 developed stress coping strategies by secreting a suite of proteins that were either unique or upregulated/overexpressed when compared to the control. The proteomics study revealed that maximum proteins that were upregulated belonged to KOG class translation, ribosomal structure, and biogenesis. This study also highlighted the enhanced expression of antioxidants superoxide dismutase, molecular chaperone heat shock proteins, and involvement of proton transporter, such as ATPase. These proteins are involved in the tolerance and detoxification of multimetals by the fungus. A. fumigatus PD-18. Therefore, this investigation on the response of cellular proteomes to multimetal stress enabled us to better understand the cellular mechanism regarding the cumulative effect of the inorganic heavy metals stress on microbes. Further, it will be conducive to screening the key genes coding for enzymes that have higher resistance to these inorganic pollutants with enhanced capability to transform pollutants.

DATA AVAILABILITY STATEMENT
The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE [1] partner repository with the dataset identifier PXD031741.

AUTHOR CONTRIBUTIONS
PD, AM, and NJ contributed to conception and design of the study. PD organized the database and wrote the first draft of the manuscript. PD and S-BH performed the statistical analysis. PD, AM, NJ, DS, and MB wrote sections of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.

FUNDING
PD thankfully acknowledges the fellowship from the Ministry of Human Resource and Development (Government of India) and German Academic Exchange Service (DAAD). MB and NJ are grateful for the funding of the UFZ for the ProMetheus platform for proteomics and metabolomics.