Next Article in Journal
The Effects of Tamoxifen on Tolerogenic Cells in Cancer
Previous Article in Journal
The First Comprehensive Biodiversity Study of Culturable Fungal Communities Inhabiting Cryoconite Holes in the Werenskiold Glacier on Spitsbergen (Svalbard Archipelago, Arctic)
Previous Article in Special Issue
An Ensemble Deep Learning Model with a Gene Attention Mechanism for Estimating the Prognosis of Low-Grade Glioma
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Putative Protein Discovery from Microalgal Genomes as a Synthetic Biology Protein Library for Heavy Metal Bio-Removal

1
Department of Biology, Faculty of Science, Chiang Mai University, Chiang Mai 50200, Thailand
2
Department of Entomology and Plant Pathology, Faculty of Agriculture, Chiang Mai University, Chiang Mai 50200, Thailand
3
Research Center in Bioresources for Agriculture, Industry and Medicine, Chiang Mai University, Chiang Mai 50200, Thailand
4
Research Center of Microbial Diversity and Sustainable Utilization, Faculty of Science, Chiang Mai University, Chiang Mai 50200, Thailand
*
Author to whom correspondence should be addressed.
Biology 2022, 11(8), 1226; https://doi.org/10.3390/biology11081226
Submission received: 13 July 2022 / Revised: 6 August 2022 / Accepted: 12 August 2022 / Published: 17 August 2022
(This article belongs to the Special Issue Synthetic Biology Aided by Artificial Intelligence)

Abstract

:

Simple Summary

Nowadays, heavy metal polluted wastewater is one of the global challenges that leads to an insufficient supply of clean water. Taking advantage of what nature has to offer, several organisms, including microalgae, can natively bioremediate these heavy metals. However, the effectiveness of such processes does not meet expectations, especially with the increasing amount of pollution in today’s world. Therefore, with the goal of creating effective strains, synthetic biology via bioengineering is widely used as a strategy to enhance the heavy metal bio-removing capability, either by directly engineering the native ability of organisms or by transferring the ability to a more suitable host. In order to do so, a list of genes or proteins involved in the processes is crucial for stepwise engineering. Yet, a large amount of information remains to be discovered. In this work, a comprehensive library of putative proteins that are involved in heavy metal bio-removal from microalgae was constructed. Moreover, with the development of machine learning, the 3D structures of these proteins are also predicted, using machine learning-based methods, to aid the use of synthetic biology further.

Abstract

Synthetic biology is a principle that aims to create new biological systems with particular functions or to redesign the existing ones through bioengineering. Therefore, this principle is often utilized as a tool to put the knowledge learned to practical use in actual fields. However, there is still a great deal of information remaining to be found, and this limits the possible utilization of synthetic biology, particularly on the topic that is the focus of the present work—heavy metal bio-removal. In this work, we aim to construct a comprehensive library of putative proteins that might support heavy metal bio-removal. Hypothetical proteins were discovered from Chlorella and Scenedesmus genomes and extensively annotated. The protein structures of these putative proteins were also modeled through Alphafold2. Although a portion of this workflow has previously been demonstrated to annotate hypothetical proteins from whole genome sequences, the adaptation of such steps is yet to be done for library construction purposes. We also demonstrated further downstream steps that allow a more accurate function prediction of the hypothetical proteins by subjecting the models generated to structure-based annotation. In conclusion, a total of 72 newly discovered putative proteins were annotated with ready-to-use predicted structures available for further investigation.

1. Introduction

Heavy metal contaminated wastewater has been a major global concern that directly affects the human population [1,2]. A considerable amount of heavy metals is released into the environment by several industries. The steel industry, for example, is a well-known source of heavy metal contamination [3]. Not only do industrial processes cause pollution, but the products from many industries, such as batteries, are also a major source of contamination [4,5]. Microalgae are promising bio-removers of contaminated heavy metals from wastewater effluents [6]. In our latest review, we summarized three mechanisms used by microalgae to bio-remove heavy metals from wastewater: biosorption, bioaccumulation and biotransformation [7]. In brief, biosorption is a process by which microalgae absorb heavy metal ions onto their cell surface; the anionic composition plays a role in attracting positive heavy metal ions (Figure 1A). This allows the cells to remove heavy metals from aqueous phases by simply attaching them to the cell surface. Interestingly, this process does not need the cells to be in an active form, meaning that dead cells can also function similarly to the living cells as long as the heavy metals can attach to the cell surface. Bioaccumulation is another process by which heavy metals are taken up by microalgal cells and later sequestered in vacuoles or specialized compartments. Lastly, biotransformation allows the conversion of toxic heavy metals to lesser or non-toxic forms using intracellular enzymes. Living cells are required for the latter two mechanisms to work continuously. Moreover, other native cellular functions, such as metal transportation and cellular responses to stress, could also impact the cell ability to bio-remove heavy metals. Heavy metal transporters also play a role in these processes as they determine the efficiency of heavy metal uptake for bioaccumulation and biotransformation. To put it simply, the more heavy metals are taken up, the greater is the possibility that they can be accumulated and transformed. Moreover, cellular responses to heavy metal-induced stresses have been reported to improve cells’ ability to remove heavy metals [8,9]. This is because the cell itself is able to tolerate higher levels of heavy metals while continuing to function at a relatively higher concentration of heavy metals.
Several microalgae have been reported to bio-remove heavy metals through the aforementioned mechanisms [6]. Chlorella and Scenedesmus are two groups of microalgae exhibiting their potential as heavy metal bio-removers. Chlorella has been well-known for its ability as a heavy metal biosorbent and, in one study, could effectively remove 220 mg of Cu+2 per one gram of cells [10]. Similarly, Scenedesmus is another group of microalgae that has been practically used for wastewater treatment in several parts of the world [11], though its application may have not been as extensive; a bio-removal capability of 574 mg/g of cells was reported from this group of microalgae [12,13]. Moreover, other microalgae such as Tetraselmis, Porphyridium, and Pseudochlorococcum have been reported to effectively bio-remove heavy metals [6,14,15], but since they are relatively underexplored, only one genome sequence of Tetraselmis and Porphyridium was reported for each genus, and none was reported for Pseudochlorococcum in the NCBI database (www.ncbi.nlm.nih.gov, accessed on 20 June 2022). As the aim of this study is to construct a comprehensive library of proteins based on reported genome sequences, Chlorella and Scenedesmus were selected, given the relatively high number of their available genome sequences.
Even though these processes in microalgae may sound as though they hold great promise, there are still limitations in the bio-removing capabilities of microalgae [16]. Thus, there is potential for improvement, particularly in the development of a supreme microalgal chassis for heavy metal contaminated wastewater treatment. Synthetic biology is currently one of the most rapidly growing biological disciplines, with the aim of using integrated knowledge of biology and bioengineering to develop or redesign biological systems with specific functionalities. To provide a related example, an engineered Chlorella sp. DT was constructed by expressing a heterologous mercuric reductase from Bacillus megaterium MB1 and showed 2-fold increase in mercury bio-removal compared with the wild type [17]. Such techniques allow us to design novel biological systems as solutions to tackle unsolved challenges.
We propose that microalgal native abilities to bio-remove heavy metals could be enhanced via synthetic biology-based approaches with the proposed strategies previously listed [7]. Generally, the study of synthetic biology follows a guideline called “Design-Build-Test-Learn” or DBTL cycle (Figure 1B), which allows systematic implementation of the new system. Interestingly, the work on synthetic biology to enhance heavy metal bio-removal is not well-established compared with other microorganisms (i.e., Escherichia coli and cyanobacteria), and indeed this is true of microalgae in general. This may be related to the low availability of genetic tools and techniques for microalgae as synthetic biology builds on genetic manipulation. More importantly, synthetic biology often applies bottom-up approaches to facilitate the design, in which genetic parts are handpicked and put together to construct novel systems. In this case, the availability of genetic parts is crucial for stepwise implementation. In addition to regulatory genetics (e.g., promoter, ribosome binding site, etc.), a library of genes responsible for significant properties of the systems is also important. To illustrate this, Figure 1C depicts the simplified workflow of a bottom-up approach for synthetic biology and the importance of genetic part libraries.
As the genetic part library—a collection of genetic parts and related genes—is essential for structured design, genes involved in bio-removal should be identified to improve microalgal abilities or even to employ such abilities in other chassis. Previously, a unicellular microalga, Chlamydomonas, has been genetically engineered to enhance its abilities as a heavy metal bio-remover. As a result, the engineered strains showed significantly improved bio-removing capability [8,9,18,19,20]. Certainly, this demonstrates the potential of employing synthetic biology to enhance the cell’s native ability for heavy metal bio-removal. However, on the basis of extremely limited data, systematic investigation is yet to be achieved, which emphasizes the importance of this work. Thus, we aim to construct a library of potential proteins involved in heavy metal bio-removal for further use via synthetic biology.
With a rapid growth of bioinformatics, a number of tools have been developed to facilitate the analysis of biological data. Basic Local Alignment Search Tool or BLAST from National Center for Biotechnology Information (NCBI) is one of the pioneer tools with a great number of users up to today. Many web interfaces designed for function prediction of the proteins from amino acid sequences and/or protein structures were also developed to allow stepwise prediction of the proteins. Pfam, InterProScan, SUPERFAMILY 2.0, and CATH are among the commonly use online tools [21,22]. Not only were the tools for functional annotation developed, well-known web servers such as SignalP were also developed to predict whether the proteins contain signal peptides or not [23]. Function prediction of hypothetical proteins from whole genome sequences has been demonstrated previously. One early example is the use of the aforementioned tools to annotate the function of all hypothetical proteins from Haemophilus influenzae Rd KW20 genome aiming to use the information retrieved for medical applications [22]. Later, in another work, hypothetical proteins from a Gram-negative bacterium, Litorilituus sediminis, were computationally predicted to be tumor-suppressors [21]. Apart from the medical point of view, hypothetical proteins that play a role in cold adaptation in Pseudomonas sp. were characterized using the similar computational workflow. However, all of these works follow the same general pattern in that they aim to annotate and characterize hypothetical proteins from a single organism, and/or only proteins involved in specific functions were targeted. This is in contrast with the present study as we aim to construct a library of putative proteins from whole groups of microalgae and to annotate as many proteins that contribute to heavy metal bio-removal as possible using our proposed approach. Moreover, the development of Alphafold [24], an accurate computational tool for protein folding prediction, caused attention to shift from traditional methods to machine learning to study protein structure. Therefore, to make the information ready to use, in this work, Alphafold2 was also used to predict the 3-dimensional (3D) structures of all putative proteins, and the structures were deposited in a protein structure database for future applications. In summary, a comprehensive library of putative proteins from microalgal genomes was constructed to allow stepwise engineering approaches to enhance heavy metal bio-removal in microalgae and potentially in other organisms.

2. Materials and Methods

A simplified workflow and computational online tools used in this study are outlined in Figure 2. The URLs of databases and tools are listed in Table S1.
To validate the annotation, three positive and negative controls each were annotated along with the putative proteins (Phase II, Figure 2) [25]. These controls were selected based on their experimental validation. The positive controls were validated proteins from microalgae that are involved in heavy metal bio-removal: natural resistance-associated macrophage protein (accession no. XP_001691702.1), ascorbate peroxidase (accession no. AAY26385.1), and glutathione peroxidase (accession no. AFI55004.1), while the negative controls are general characterized proteins that function in cellular processes: elongation factor EF-2 (accession no. NP_001321033.1), transcriptional regulator (accession no. WP_097343503.1), and cytochrome C peroxidase (accession no. WP_016160016.1).

2.1. Target Protein Identification and Sequence Retrieval

The terms given in Table S2 were used as inputs to search for proteins in UniprotKB database. The filter “reviewed” is also applied to obtain only the proteins that had been characterized or computationally curated. All search hits were manually confirmed that all proteins selected from this stage possess the expected function, as the search hits may also result in different proteins with the same gene/protein abbreviation. All selected sequences were exported as template sequences for the next step in a FASTA format.

2.2. Sequence Similarity Search

NCBI BLASTp function was used for sequence similarity search, using the obtained sequences as templates against the genomes of two microalgal groups namely Chlorella and Scenedesmus (Table 1). Non-redundant protein sequences (nr) database and blastp (protein-protein BLAST) algorithm were selected. Once the search results from BLASTp showed hits for each group of proteins, only hypothetical protein sequences with an E-value of less than 1 × 10−10 were selected. When searching a database, the E-value is a parameter that indicates how many hits are likely to occur by chance: the lower the E-value, the more significant the match is. It is recommended that the E-value between 1 × 10−10–1 × 10−50 should, at least, allow a domain match. The E-value between 1 × 10−50–1 × 10−100 indicates almost identical sequences, and if the E-value is less than 1 × 10−100, it suggests identical sequences as described in the Qiagen handbook [26]. However, it is important to note that E-value cutoffs should be considered on a case-by-case basis as they are dependent on the length of the sequence and the size of the databases.

2.3. Function Prediction

To primarily predict the functions of the retrieved protein sequences, four different webservers, namely, Pfam, InterProScan, SUPERFAMILY 2.0, CATH, were selected (Table S1). All webservers were previously used to annotate hypothetical proteins from whole genome sequences with a comparative assessment of each tool [29].

2.4. Signal Peptide Prediction

Signal peptide prediction is used to confirm the primary annotation as the target proteins in this work localize specifically, either as transporters, which localize on the cell membrane, or enzymes that function intracellularly. SignalP is one of the most used webservers for signal peptide identification. This server has recently released its latest version SignalP 6.0 early this year and claimed its ability to detect all types of signal peptides [30].

2.5. Structure Modeling

Alphafold2 is used for structure prediction in this study. Alphafold is a recent innovative tool to accurately model the protein structures based on the provided protein sequences [24]. With an increased interest in Alphafold, a recent work aiming to make the structure modeling available to all has recently been published [31] and the online tool is free to use under the name ColabFold. In this study, a related version of Alphafold via Google Colab notebook was used (Table S1). The analysis was performed with default parameters and Amber-Relax applied. The genetic database, mmseqs2, was selected. No custom MSA was uploaded to any runs. The filter option was left unchanged since the raw hypothetical protein was run without any trims. The Alphafold was run using pLDDT metric with 512:1024 max msa. Five models were run for each template with the use of ptm to fine-tune the model parameters. Each model was fed back to the neutral network three times. All generated protein structures were deposited in ModelArchive with a list of unique DOI as shown in Table S4.

2.6. Structure-Based Functional Annotation

All protein structures were investigated further to validate the primary annotation in the previous step. ProFunc was used for this purpose. This tool was designed to annotate the proteins based on sequence and structure inputs [32]. In our case, as mentioned, the structures generated from Alphafold2 were used as the template. Option “reverse templates” was examined and the E-value of at least 1 × 10−6 was considered, according to the server.

3. Results

3.1. Target Protein Identification

Microalgae perform three main mechanisms that allow them to bio-remove heavy metals from wastewater effluents. Previously, we have compiled a list of genes with reported use for genetic engineering to enhance microalgal bio-removing capacity [7]. Thus, in this work, we utilize the list as a primary source for protein targets. As part of our effort to investigate all relevant aspects toward heavy metal bio-removal, we have broadened the scope of the search by including more proteins from the literature. A total of 26 protein targets are shown in Table S2, listed separately based on their function in each mechanism. As for bioaccumulation, phytochelatins play an important role in binding with heavy metals, and enzymes involved in the synthesis of phytochelatins namely glutamate cysteine ligase, phytochelatin synthase, and glutathione synthetase [33] were therefore targeted. Although another type of heavy metal binding proteins, metallothioniens, also play a role in the same mechanism, they are directly transcribed and translated from nucleotide sequences as cysteine-rich short peptides [33], which are difficult to search against the genomes. Therefore, this type of heavy metal binding protein is not included in this study. Reductases are a major group of enzymes that detoxify heavy metals by converting them into their less or non-toxic derivatives [34,35]. In our list, three reductases specific to mercury, chromium, and arsenic were used to search against the microalgal genomes, as these enzymes have been reported to alleviate the toxicity of heavy metals in microalgae [18,36]. As previously described, enhanced metal transportation is another strategy that allows better performance of bio-removal, especially to facilitate the bioaccumulation and biotransformation in which enzymes function intracellularly. Several types of transporters have been reported to be responsible for microalgal metal transportation [37]. In this work, we selected a few representatives of heavy metal transporters as our templates. Similarly, five different enzymes were selected as representatives of cellular stress responses. All of these enzymes were reported to take part in cellular responses to heavy metal toxicity in microalgae [38]. As mentioned, biosorption is a mechanism by which heavy metal ions are attracted onto the cell surface of microalgae and, in turn, removed upon removal of the microalgae. However, this mechanism is not a direct effect from proteins synthesized by the cells, but rather indirectly from the composition of the cell surface. To enhance biosorption efficiency, the strategy may include cell surface modification to make the surface components more ionic and attract more heavy metal ions. Such modifications could be achieved by, for example, addition of CXXEE motif onto the cell surface components [39] or inducing biofilm formation [40]. Therefore, this mechanism is not included in our list. The proteins and keywords used are listed in Table S2. However, it should be noted that when using UniprotKB database or other keyword-based tools, the results are not exclusive for some keywords as different meanings can be inferred. To state the obvious, the keyword “inorganic phosphate transporter”, which is a family name of the protein transporters involved in import-export systems of phosphate [41], resulted in a total of 712 hits from UniprotKB search. However, manual selection only showed 52 hits of the inorganic phosphate transporters, the rest of the hit results showed other related proteins such as vacuolar transporter chaperones, glucose-6-phosphate exchangers, and other ATP-binding cassettes.
According to the UniprotKB database, we collectively report the number of hits for each group of proteins (bioaccumulation, biotransformation, heavy metal transporters, and cellular stress responses) as shown in Figure 3.

3.2. Sequence Similarity Search

The sequences of all 27 groups of protein targets were obtained and used as templates for BLASTp against two microalgal genomes, Chlorella (taxid 3071) and Scenedesmus (taxid 3087). Hypothetical proteins with an E-value lower than 1 × 10−10 were obtained. The number of matches is shown in Table 2.

3.3. Protein Function Prediction

The NCBI accession numbers for all protein matches from BLASTp are listed in Table 3. Function and signal peptide prediction for each hypothetical protein were investigated. As we aim to only annotate the uncharacterized proteins, only sequences described as hypothetical proteins were selected and run through function prediction servers as described above. The results for all hypothetical proteins are listed in Table S3. Each bioinformatic server has its own strength and weakness. Thus, using multiple servers would result in a more accurate prediction. A recent publication has compared several webservers for protein prediction and the results revealed that Pfam and InterPro showed the highest scores among the compared webservers [29]. Pfam is a widely used protein family database and tool that is still active with a recent update in 2021 [42]. InterProScan is another well-known functional classification tool that is based on several databases [43]. SUPERFAMILY, another webserver used in this study, was also considered the second-best server according to the aforementioned study. This tool was designed to predict both superfamily and family of the protein sequences [44]. In the same comparative report [29], SBase was also used and demonstrated to have high overall scores for accuracy, sensitivity, specificity, and ROC analysis, yet the server was last updated in 2006; therefore, SBase was not included in the present study. Lastly, although CATH did not score much in the aforementioned comparison especially regarding the specificity, according to our predictions, CATH produced relatively specific results (Table S3). To provide an explanation, when CATH was used to predict the function of a template (accession no. XP_005845237.1) expected to be a manganese transporter, CATH predicted the protein as “ABC transporter G family member 22”, whilst Pfam predicted as “ABC transporter”. Moreover, CATH is the only webtool that allows a specific annotation of copper-transporting ATPase (Table S3). Pfam, on the other hand, predicted the function of ATPases separately for different protein domains. For example, copper-transporting ATPase shows 4 predicted functions: (1) cation transporting ATPase, C-terminus, (2) E1-E2 ATPase, (3) haloacid dehalogenase-like hydrolase (CL0137), and (4) cation transporter/ATPase, N-terminus, which are basic structures for ATPase activity [45]. Similarly, InterProScan also resulted in a list of molecular functions that could be inferred as ATPase rather than stating that the proteins are ATPase (Table S3).
Signal peptide prediction was used to double-validate the predicted function of the hypothetical proteins. SignalP is a well-reputed webserver for signal peptide prediction. A report on the comparison of signal peptide prediction suggests that the use of SignalP 4.1 is most consistent compared with the other versions [46]. However, the most recent version, SignalP 6.0, was released after that [30]. Interestingly, the results showed that most of the proteins did not contain signal peptides, even though some of them were expected to be transporters (Table 3).

3.4. Homology Modeling

All hypothetical protein sequences primarily characterized to have the same function and signal peptide as their templates were modeled using Alphafold2. Amber-relax was applied to generate more accurate models. Examples of protein structures generated are presented in Figure 4. It should be noted that a hypothetical protein from Scenedesmus sp. PABB004 (accession no. KAF8061310.1) is 1553 amino acid long; therefore, this protein was not modeled through Alphafold2 as the recommended longest sequence was 1400 amino acids. When using Alphafold, the confidence measure of the models can be evaluated using a pLDDT score (0–100), by which the guidance is as follows; regions with a pLDDT score of more than 90 are considered highly accurate. Regions with a pLDDT score between 70 and 90 are considered generally good and regions with a pLDDT score between 50 and 70 are considered with low confidence. In our case, out of 72 models, 31 models showed the average pLDDT scores of more than 90, 31 models showed the average scores between 70 and 90, and 10 models showed the average scores between 50 and 70 (Table S4). All generated models were deposited in ModelArchive with the accession ID listed in Table S4. It is also important to note that the commonly used Protein Data Bank (PDB) currently only accepts experimental model depositions [47].

3.5. Structure-Based Annotation

Using the created models as a resource, ProFunc—a structure-based functional annotation—was used to annotate the proteins with structure inputs. ProFunc revealed that the majority of the annotated proteins possess the same functions as primarily annotated in Table 3. However, 7 structures out of 72 protein structures were annotated with different functions as shown in Table S5. To elaborate, two structures with manganese transporter domain (MntA) putative functions were annotated as ATP-bound human transporter found in retina and human sterol transporter, respectively. Two zinc-regulated/zinc transporters were annotated as NH3 transporter from Nitrosomonas europaea and multidrug transporter from Lactobacillus lactis. One ascorbate peroxidase was annotated as cytochrome C peroxidase and two glutathione S-transferases were annotated as apo-dehydroascorbate reductase and glutathione-bound dehydroascorbate reductase, respectively.

4. Discussion

The rise in the number of whole genome sequences prompts a rapid development of computational tools for the analysis of this available data. A webtool was previously developed to help predict the functions of genes from microalgal genomes and is currently active via http://pathways.mcdb.ucla.edu/algal/index.html (accessed on 23 May 2022); however, the database’s scope is limited only to two microalgae, Chlamydomonas reinhardtii and Chlorella NC64A. Thus, other means of methods are still useful, particularly if the genes or proteins from other microalgal strains are in question.
In our work, it is obvious that BLASTp resulted in greater numbers of hits against Chlorella than Scenedesmus. The most straightforward explanation for this is because the number of reported whole genome sequences of Chlorella is significantly higher (21 genomes; www.ncbi.nlm.nih.gov, accessed on 23 May 2022), in comparison with that of Scenedesmus (6 genomes; accessed on 23 May 2022). Moreover, Chlorella is considered a frequently used microalga for wastewater treatment, despite the increased attention toward Scenedesmus. Noticeably, no enzymes implicated in the bioaccumulation mechanism were found from Scenedesmus genomes. This could be due to the explanation stated above, or it may suggest that the Scenedesmus does not actually use this mechanism. However, it should also be emphasized that glutathione peroxidase, which is a common enzyme found in most organisms to protect the cell from oxidative damage [49], was also not found in Scenedesmus. This may suggest that the use of BLASTp to retrieve proteins from microalgal genomes could be further improved.
Available webtools for functional annotation are useful, especially to primarily screen for protein sequences with particular functions. Interestingly, protein sequences retrieved when using reductase enzymes as templates were not explicitly identified as reductases. Although two proteins (accession no. KAG7668560.1 and XP_005845177.1) were identified by CATH as thioredoxin reductase, which is often seen coupled with arsenate reductases [50,51], the other servers identified otherwise. Moreover, it is undeniably challenging to predict the functions of transporters as they consist of several domains to form a functional protein [37,52] and they are often broadly specific toward different substrates [53]. This was illustrated clearly from BLASTp search where all cadmium-, zinc- and lead-transporting ATPase showed the same hit results with copper-transporting ATPase (Table S3), suggesting that their structures are similar, and they are broadly specific to several heavy metals. All heavy metal transporting ATPases are classified as P-type ATPases and contain (1) E1-E2 ATPase (2) hydrolase (3) cation ATPase domains [54]. In this regard, Pfam and InterProScan broadly annotate the molecular functions of different regions on the hypothetical proteins as ATP hydrolysis activity, ATP binding, transporter activity, nucleotide binding (Table S5). Though this fits the definition of P-ATPase, it suggests that the specificity of these tools could be further improved. We suggest that the use of these servers is still applicable, but further interpretation is required. Noticeably, the results from SUPERFAMILY 2.0 and CATH showed several calcium-transporting ATPase hits for all heavy metal transporting ATPase, which is because heavy metal transporting ATPase could also function as a Ca2+ pump [55].
Alphafold is a breakthrough computational method for protein structure prediction with atomic level accuracy. Using coding-based platform, this makes it challenging for researchers outside of the computational fields to apply Alphafold to their studies. Recently, a few reports, including from the Alphafold creators themselves, have developed relatively more user-friendly Alphafold on Google Colab platform [31]. With such widely accessible platforms, the number of models generated from Alphafold is anticipated to sharply increase in the near future. In our work, we observed that a majority of the predicted protein structures showed long amino acid chains that do not form secondary structures with the rest of the protein region (Figure 4), which could be because the actual start points of the protein sequences retrieved from BLASTp were not accurately identified. This instantly brought us back to the signal peptide prediction from SignalP server, as this tool only allows the prediction of the first 70 amino acids from the N-terminus, which means that if the hypothetical protein sequences retrieved from NCBI contain more than 70 amino acids upstream of the actual start point, the prediction using SignalP would not be accurate.
Though Alphafold has its own protein database (https://alphafold.ebi.ac.uk/ (accessed on 20 June 2022)), the submission is not made available. Therefore, in this study, we submitted our Alphafold models to ModelArchive, which is the only database that allows submission of modeling-based protein structures.
Structure-based annotation is considered a more accurate method for predicting protein function than homology-based prediction, as the latter often considers only the sequence similarity of the input proteins and their homologs, while the former also considers other factors, including protein binding pocket [56]. In this work, ProFunc is used as a method to validate the protein prediction from primary screening. The results were not surprising, as most of the proteins were annotated to have the same function as from homology-based prediction. This assures that the homology-based annotation could be used to some extent. Interestingly, when looking at seven proteins with different annotation results, it was noticed that the transporters were also annotated as transporters but with specified substrates/ligands. It is important to point out that since homology-based annotation only allows a broad annotation of transporters as unspecified transporters, structure-based annotation allows the substrates/ligands of the transporters to also be identified. Unfortunately, in these cases, the substrates identified were not heavy metals as expected. In the case of enzymes, however, predicted ascorbate peroxidase from homology-based prediction was predicted as cytochrome C peroxidase in this step. This could be linked to their highly similar protein structures, even though they have different substrate-binding sites [57]. Therefore, this observation demonstrates that structure-based annotation allows a more accurate prediction of the protein function, especially when the proteins have similar structures. Two structures of glutathione S-transferases were annotated as apo-dehydroascorbate reductase and glutathione-bound dehydroascorbate reductase. The dehydroascorbate reductase (DHRA) is a member of glutathione S-transferase superfamily [58] suggesting that ProFunc could identify subgroups of the proteins in some cases. Moreover, the challenge encountered when using Pfam and InterProscan in the primary annotation was not presented when using ProFunc. To give an example, instead of predicting a heavy metal transporting ATPase (accession no. XP_005851032.1) separately for each region as stated above, ProFunc predicted the protein to be similar to zinc-transporting PIB-type ATPase. Overall, the use of structure-based annotation is a useful approach that accurately annotates the protein sequences, especially in our case, when investigating the proteins with similar structures. Furthermore, it allows the annotation of transporters as a whole protein rather than separated domains. However, it should also be taken into account that out of all 72 protein structures used as templates for structure-based annotation via ProFunc, only 31 of them were considered highly accurate predictions according to the Alphafold2 confidence measures.
To validate the proposed approach, natural resistance-associated macrophage protein, ascorbate peroxidase, and glutathione peroxidase with experimentally confirmed activities were used as positive controls. The results showed that after extensive annotation following our workflow, the protein sequences were annotated as expected. Experimentally validated elongation factor EF-2, transcriptional regulator, and cytochrome C peroxidase were selected as negative controls. These proteins function in general cellular processes, and in the case of cytochrome C peroxidase, it has a very similar structure with ascorbate peroxidase and is often shown as matches when ascorbate peroxidase is blasted. Therefore, this choice of negative control would allow the assessment of our approach even when influenced by these factors. The results showed that these proteins were annotated to possess their validated functions, which demonstrates that the workflow presented here could distinguish between the two very similar-structured proteins.
Ultimately, the applications of our protein library could range from selection of the proteins from our library for a single engineering design to combinatorial optimization of all genetic parts. At present, the construction of synthetic biology-based systems mostly relies on previously reported genes or proteins, which limits the possibility of the design and the optimization. For example, engineering of a microalga to enhance heavy metal bioremediation was recently reported and the work utilized gene sequences from Arabidopsis that were experimentally characterized prior to the study [18]. This limits the selection of the genes to only the reported ones and raises the question whether the reported genes are the best homologs that there are. Therefore, selection of the genes from our library could provide more choices and allow a stepwise design [59]. Moreover, the predicted protein structures from our work could be used for further computational analysis to predict the activity of the proteins and compare homologs prior to the actual engineering step. Molecular docking is a recent field of protein study that investigates the compatibility of proteins and their ligands or substrates, which can be used to predict the activity of the proteins. To provide a related example, a recent work performed molecular docking to confirm the binding between chromium ion and cell surface proteins in biosorption mechanism [60].

5. Conclusions

In this work, we aim to construct a comprehensive library of putative proteins from Chlorella and Scenedemus genomes to facilitate synthetic biology-based engineering for heavy metal bio-removal. We selected 27 different groups of protein targets based on their reported capability as proteins involved in heavy metal bio-removing mechanisms. We found a total of 72 putative proteins. Among these, 65 were exclusively annotated to possess the same functions as their templates. The protein structures of all annotated proteins were also generated and deposited in the protein structure archive for any further use. This study, therefore, provides a putative protein library that could be used as a database for synthetic biologists to handpick the proteins for engineering purposes with readily available structures for additional investigations. Nonetheless, it should be kept in mind that some limitations remain when using our approach for putative protein discovery. First, as encountered in our case, the relatively low number of reported microalgal genomes limits the number of discoverable proteins. Second, the current platform of Alphafold2 only allows up to 1400 amino acid long proteins as templates for modeling. This restriction is another factor to be considered if large proteins are in question. If these limitations are mitigated, a larger number of putative proteins are anticipated to be discovered and annotated using our approach.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biology11081226/s1, Table S1: List of bioinformatic webtools used in this study; Table S2: A list of target proteins in this study and their search results from UniprotKB; Table S3: Function prediction using different webservers; Table S4: Protein structure archive ID; Table S5: Structure-based annotation of the putative proteins.

Author Contributions

Conceptualization, P.S.; investigation, T.U., N.M., N.C. and P.S.; validation, T.U. and P.S.; writing—original draft preparation, P.S.; writing—review and editing, T.U., N.M., N.C., W.P.-A., J.P., C.P. and P.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research work was funded by The Murata Science Foundation and partially supported by Chiang Mai University.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data associated with this research is included in this article and its supplementary information. Any further information is available upon reasonable request.

Acknowledgments

This research work was supported by The Murata Science Foundation and Chiang Mai University. The protein structures were visualized by UCSF ChimeraX, developed by the Resource for Biocomputing, Visualization, and Informatics at the University of California, San Francisco.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jaishankar, M.; Tseten, T.; Anbalagan, N.; Mathew, B.B.; Beeregowda, K.N. Toxicity, Mechanism and Health Effects of Some Heavy Metals. Interdiscip. Toxicol. 2014, 7, 60–72. [Google Scholar] [CrossRef] [PubMed]
  2. Briffa, J.; Sinagra, E.; Blundell, R. Heavy Metal Pollution in the Environment and Their Toxicological Effects on Humans. Heliyon 2020, 6, e04691. [Google Scholar] [CrossRef] [PubMed]
  3. Namuhani, N.; Kimumwe, C. Soil Contamination with Heavy Metals around Jinja Steel Rolling Mills in Jinja Municipality, Uganda. J. Health Pollut. 2015, 5, 61–67. [Google Scholar] [CrossRef] [PubMed]
  4. Islam, M.T.; Huda, N.; Baumber, A.; Hossain, R.; Sahajwalla, V. Waste Battery Disposal and Recycling Behavior: A Study on the Australian Perspective. Environ. Sci. Pollut. Res. 2022. [Google Scholar] [CrossRef]
  5. Krekeler, M.P.S.; Barrett, H.A.; Davis, R.; Burnette, C.; Doran, T.; Ferraro, A.; Meyer, A. An Investigation of Mass and Brand Diversity in a Spent Battery Recycling Collection with an Emphasis on Spent Alkaline Batteries: Implications for Waste Management and Future Policy Concerns. J. Power Sources 2012, 203, 222–226. [Google Scholar] [CrossRef]
  6. Suresh Kumar, K.; Dahms, H.U.; Won, E.J.; Lee, J.S.; Shin, K.H. Microalgae—A Promising Tool for Heavy Metal Remediation. Ecotoxicol. Environ. Saf. 2015, 113, 329–352. [Google Scholar] [CrossRef]
  7. Sattayawat, P.; Yunus, I.S.; Noirungsee, N.; Mukjang, N.; Pathom-Aree, W.; Pekkoh, J.; Pumas, C. Synthetic Biology-Based Approaches for Microalgal Bio-Removal of Heavy Metals from Wastewater Effluents. Front. Environ. Sci. 2021, 9, 1–12. [Google Scholar] [CrossRef]
  8. Wei, Y.Y.; Zheng, Q.; Liu, Z.P.; Yang, Z.M. Regulation of Tolerance of Chlamydomonas Reinhardtii to Heavy Metal Toxicity by Heme Oxygenase-1 and Carbon Monoxide. Plant Cell Physiol. 2011, 52, 1665–1675. [Google Scholar] [CrossRef]
  9. Siripornadulsil, S.; Traina, S.; Verma, D.P.S.; Sayre, R.T. Molecular Mechanisms of Proline-Mediated Tolerance to Toxic Heavy Metals in Transgenic Microalgae. Plant Cell 2002, 14, 2837–2847. [Google Scholar] [CrossRef]
  10. Doshi, H.; Ray, A.; Kothari, I.L.; Gami, B. Spectroscopic and Scanning Electron Microscopy Studies of Bioaccumulation of Pollutants by Algae. Curr. Microbiol. 2006, 53, 148–157. [Google Scholar] [CrossRef]
  11. Peña-Castro, J.M.; Martínez-Jerónimo, F.; Esparza-García, F.; Cañizares-Villanueva, R.O. Heavy Metals Removal by the Microalga Scenedesmus incrassatulus in Continuous Cultures. Bioresour. Technol. 2004, 94, 219–222. [Google Scholar] [CrossRef] [PubMed]
  12. Terry, P.A.; Stone, W. Biosorption of Cadmium and Copper Contaminated Water by Scenedesmus abundans. Chemosphere 2002, 47, 249–255. [Google Scholar] [CrossRef]
  13. Monteiro, C.M.; Castro, P.M.L.; Malcata, X.X. Use of the Microalga Scenedesmus obliquus to Remove Cadmium Cations from Aqueous Solutions. World J. Microbiol. Biotechnol. 2009, 25, 1573–1578. [Google Scholar] [CrossRef]
  14. da Costa, A.C.A.; de França, F.P. The Behaviour of the Microalgae Tetraselmis chuii in Cadmium-Contaminated Solutions. Aquac. Int. 1998, 6, 57–66. [Google Scholar] [CrossRef]
  15. Shanab, S.; Essa, A.; Shalaby, E. Bioremoval Capacity of Three Heavy Metals by Some Microalgae Species (Egyptian Isolates). Plant Signal. Behav. 2012, 7, 392–399. [Google Scholar] [CrossRef]
  16. Yan, C.; Qu, Z.; Wang, J.; Cao, L.; Han, Q. Microalgal Bioremediation of Heavy Metal Pollution in Water: Recent Advances, Challenges, and Prospects. Chemosphere 2022, 286, 131870. [Google Scholar] [CrossRef]
  17. Huang, C.C.; Chen, M.W.; Hsieh, J.L.; Lin, W.H.; Chen, P.C.; Chien, L.F. Expression of Mercuric Reductase from Bacillus megaterium MB1 in Eukaryotic Microalga Chlorella Sp. DT: An Approach for Mercury Phytoremediation. Appl. Microbiol. Biotechnol. 2006, 72, 197–205. [Google Scholar] [CrossRef]
  18. Ibuot, A.; Webster, R.E.; Williams, L.E.; Pittman, J.K. Increased Metal Tolerance and Bioaccumulation of Zinc and Cadmium in Chlamydomonas reinhardtii Expressing a AtHMA4 C-Terminal Domain Protein. Biotechnol. Bioeng. 2020, 117, 2996–3005. [Google Scholar] [CrossRef]
  19. Ibuot, A.; Dean, A.P.; McIntosh, O.A.; Pittman, J.K. Metal Bioremediation by CrMTP4 Over-Expressing Chlamydomonas reinhardtii in Comparison to Natural Wastewater-Tolerant Microalgae Strains. Algal Res. 2017, 24, 89–96. [Google Scholar] [CrossRef]
  20. Han, S.H.; Hu, Z.L.; Lei, A.P. Expression and Function Analysis of the Metallothionein-like (MT-like) Gene from Festuca rubra in Chlamydomonas reinhardtii Chloroplast. Sci. China C Life Sci. 2008, 51, 1076–1081. [Google Scholar] [CrossRef]
  21. Kader, M.A.; Ahammed, A.; Khan, M.S.; al Ashik, S.A.; Islam, M.S.; Hossain, M.U. Hypothetical Protein Predicted to Be Tumor Suppressor: A Protein Functional Analysis. Genom. Inform. 2022, 20, e6. [Google Scholar] [CrossRef] [PubMed]
  22. Shahbaaz, M.; Hassan, M.I.; Ahmad, F. Functional Annotation of Conserved Hypothetical Proteins from Haemophilus influenzae Rd KW20. PLoS ONE 2013, 8, e84263. [Google Scholar] [CrossRef] [PubMed]
  23. Nielsen, H.; Krogh, A. Prediction of Signal Peptides and Signal Anchors by a Hidden Markov Model SignalP View Project. In Proceedings of the Sixth International Conference on Intelligent Systems for Molecular Biology, Montreal, QC, Canada, 28 June–1 July 1998; AAAI Press: Palo Alto, CA, USA; Volume 6, pp. 122–130. [Google Scholar]
  24. Jumper, J.; Evans, R.; Pritzel, A.; Green, T.; Figurnov, M.; Ronneberger, O.; Tunyasuvunakool, K.; Bates, R.; Žídek, A.; Potapenko, A.; et al. Highly Accurate Protein Structure Prediction with AlphaFold. Nature 2021, 596, 583–589. [Google Scholar] [CrossRef] [PubMed]
  25. Kootery, K.P.; Sarojini, S. Structural and Functional Characterization of a Hypothetical Protein in the RD7 Region in Clinical Isolates of Mycobacterium tuberculosis—An in Silico Approach to Candidate Vaccines. J. Genet. Eng. Biotechnol. 2022, 20, 55. [Google Scholar] [CrossRef]
  26. QIAGEN. CLC Genomics Workbench User Manual. Available online: https://resources.qiagenbioinformatics.com/manuals/clcgenomicsworkbench/801/User_Manual.pdf (accessed on 2 August 2022).
  27. Mubashar, M.; Naveed, M.; Mustafa, A.; Ashraf, S.; Baig, K.S.; Alamri, S.; Siddiqui, M.H.; Zabochnicka-światek, M.; Szota, M.; Kalaji, H.M. Experimental Investigation of Chlorella Vulgaris and Enterobacter Sp. Mn17 for Decolorization and Removal of Heavy Metals from Textile Wastewater. Water 2020, 12, 3034. [Google Scholar] [CrossRef]
  28. Shivaji, S.; Dronamaraju, S.V.L. Scenedesnus Rotundus Isolated from the Petroleum Effluent Employs Alternate Mechanisms of Tolerance to Elevated Levels of Cadmium and Zinc. Sci. Rep. 2019, 9, 8485. [Google Scholar] [CrossRef]
  29. Da Rosa Tavares, B.A.; Paes, J.A.; Zaha, A.; Ferreira, H.B. Reannotation of Mycoplasma Hyopneumoniae Hypothetical Proteins Revealed Novel Potential Virulence Factors. Microb. Pathog. 2022, 162, 105344. [Google Scholar] [CrossRef]
  30. Teufel, F.; Almagro Armenteros, J.J.; Johansen, A.R.; Gíslason, M.H.; Pihl, S.I.; Tsirigos, K.D.; Winther, O.; Brunak, S.; von Heijne, G.; Nielsen, H. SignalP 6.0 Predicts All Five Types of Signal Peptides Using Protein Language Models. Nat. Biotechnol. 2022, 40, 1023–1025. [Google Scholar] [CrossRef]
  31. Mirdita, M.; Schütze, K.; Moriwaki, Y.; Heo, L.; Ovchinnikov, S.; Steinegger, M. ColabFold: Making Protein Folding Accessible to All. Nat. Methods 2022, 19, 679–682. [Google Scholar] [CrossRef]
  32. Laskowski, R.A.; Watson, J.D.; Thornton, J.M. ProFunc: A Server for Predicting Protein Function from 3D Structure. Nucleic Acids Res. 2005, 33, W89–W93. [Google Scholar] [CrossRef]
  33. Balzano, S.; Sardo, A.; Blasio, M.; Chahine, T.B.; Dell’Anno, F.; Sansone, C.; Brunet, C. Microalgal Metallothioneins and Phytochelatins and Their Potential Use in Bioremediation. Front. Microbiol. 2020, 11, 517. [Google Scholar] [CrossRef] [PubMed]
  34. Wang, H.; He, Z.; Luo, L.; Zhao, X.; Lu, Z.; Luo, T.; Li, M.; Zhang, Y. An Aldo-Keto Reductase, Bbakr1, Is Involved in Stress Response and Detoxification of Heavy Metal Chromium but Not Required for Virulence in the Insect Fungal Pathogen Beauveria bassiana. Fungal Genet. Biol. 2018, 111, 7–15. [Google Scholar] [CrossRef] [PubMed]
  35. Singh, S.; Kumar, V. Mercury Detoxification by Absorption, Mercuric Ion Reductase, and Exopolysaccharides: A Comprehensive Study. Environ. Sci. Pollut. Res. 2020, 27, 27181–27201. [Google Scholar] [CrossRef] [PubMed]
  36. Danouche, M.; el Ghachtouli, N.; el Arroussi, H. Phycoremediation Mechanisms of Heavy Metals Using Living Green Microalgae: Physicochemical and Molecular Approaches for Enhancing Selectivity and Removal Capacity. Heliyon 2021, 7, e07609. [Google Scholar] [CrossRef]
  37. Blaby-Haas, C.E.; Merchant, S.S. The Ins and Outs of Algal Metal Transport. Biochim. Biophys. Acta-Mol. Cell Res. 2012, 1823, 1531–1552. [Google Scholar] [CrossRef]
  38. Danouche, M.; el Ghatchouli, N.; Arroussi, H. Overview of the Management of Heavy Metals Toxicity by Microalgae. J. Appl. Phycol. 2022, 34, 475–488. [Google Scholar] [CrossRef]
  39. Kotrba, P.; Ruml, T. Surface Display of Metal Fixation Motifs of Bacterial Pi-Type ATPases Specifically Promotes Biosorption of Pb2+ by Saccharomyces cerevisiae. Appl. Environ. Microbiol. 2010, 76, 2615–2622. [Google Scholar] [CrossRef]
  40. Ma, X.; Chen, Y.; Liu, F.; Zhang, S.; Wei, Q. Enhanced Tolerance and Resistance Characteristics of Scenedesmus obliquus FACHB-12 with K3 Carrier in Cadmium Polluted Water. Algal Res. 2021, 55, 102267. [Google Scholar] [CrossRef]
  41. Bøttger, P.; Pedersen, L. Evolutionary and Experimental Analyses of Inorganic Phosphate Transporter PiT Family Reveals Two Related Signature Sequences Harboring Highly Conserved Aspartic Acids Critical for Sodium-Dependent Phosphate Transport Function of Human PiT2. FEBS J. 2005, 272, 3060–3074. [Google Scholar] [CrossRef]
  42. Mistry, J.; Chuguransky, S.; Williams, L.; Qureshi, M.; Salazar, G.A.; Sonnhammer, E.L.L.; Tosatto, S.C.E.; Paladin, L.; Raj, S.; Richardson, L.J.; et al. Pfam: The Protein Families Database in 2021. Nucleic Acids Res. 2021, 49, D412–D419. [Google Scholar] [CrossRef]
  43. Blum, M.; Chang, H.Y.; Chuguransky, S.; Grego, T.; Kandasaamy, S.; Mitchell, A.; Nuka, G.; Paysan-Lafosse, T.; Qureshi, M.; Raj, S.; et al. The InterPro Protein Families and Domains Database: 20 Years on. Nucleic Acids Res. 2021, 49, D344–D354. [Google Scholar] [CrossRef] [PubMed]
  44. Pandurangan, A.P.; Stahlhacke, J.; Oates, M.E.; Smithers, B.; Gough, J. The SUPERFAMILY 2.0 Database: A Significant Proteome Update and a New Webserver. Nucleic Acids Res. 2019, 47, D490–D494. [Google Scholar] [CrossRef] [PubMed]
  45. Wang, K.; Sitsel, O.; Meloni, G.; Autzen, H.E.; Andersson, M.; Klymchuk, T.; Nielsen, A.M.; Rees, D.C.; Nissen, P.; Gourdon, P. Structure and Mechanism of Zn2+-Transporting P-Type ATPases. Nature 2014, 514, 518–522. [Google Scholar] [CrossRef] [PubMed]
  46. Garcion, C.; Béven, L.; Foissac, X. Comparison of Current Methods for Signal Peptide Prediction in Phytoplasmas. Front. Microbiol. 2021, 12, 661524. [Google Scholar] [CrossRef] [PubMed]
  47. Schwede, T. Protein Modeling: What Happened to the “Protein Structure Gap”? Structure 2013, 21, 1531–1540. [Google Scholar] [CrossRef] [PubMed]
  48. Pettersen, E.F.; Goddard, T.D.; Huang, C.C.; Meng, E.C.; Couch, G.S.; Croll, T.I.; Morris, J.H.; Ferrin, T.E. UCSF ChimeraX: Structure Visualization for Researchers, Educators, and Developers. Protein Sci. 2021, 30, 70–82. [Google Scholar] [CrossRef] [PubMed]
  49. Tamaki, S.; Mochida, K.; Suzuki, K. Diverse Biosynthetic Pathways and Protective Functions against Environmental Stress of Antioxidants in Microalgae. Plants 2021, 10, 1250. [Google Scholar] [CrossRef]
  50. Shi, L.D.; Chen, Y.S.; Du, J.J.; Hu, Y.Q.; Shapleigh, J.P.; Zhao, H.P. Metagenomic Evidence for a Methylocystis Species Capable of Bioremediation of Diverse Heavy Metals. Front. Microbiol. 2019, 10, 3297. [Google Scholar] [CrossRef]
  51. Li, Y.; Hu, Y.; Zhang, X.; Xu, H.; Lescop, E.; Xia, B.; Jin, C. Conformational Fluctuations Coupled to the Thiol-Disulfide Transfer between Thioredoxin and Arsenate Reductase in Bacillus subtilis. J. Biol. Chem. 2007, 282, 11078–11083. [Google Scholar] [CrossRef]
  52. Pedersen, B.P.; Buch-Pedersen, M.J.; Preben Morth, J.; Palmgren, M.G.; Nissen, P. Crystal Structure of the Plasma Membrane Proton Pump. Nature 2007, 450, 1111–1114. [Google Scholar] [CrossRef]
  53. Kühlbrandt, W. Biology, Structure and Mechanism of P-Type ATPases. Nat. Rev. Mol. Cell Biol. 2004, 5, 282–295. [Google Scholar] [CrossRef]
  54. Migocka, M. Copper-Transporting ATPases: The Evolutionarily Conserved Machineries for Balancing Copper in Living Systems. IUBMB Life 2015, 67, 737–745. [Google Scholar] [CrossRef] [PubMed]
  55. Moreno, I.; Norambuena, L.; Maturana, D.; Toro, M.; Vergara, C.; Orellana, A.; Zurita-Silva, A.; Ordenes, V.R. AtHMA1 Is a Thapsigargin-Sensitive Ca2+/Heavy Metal Pump. J. Biol. Chem. 2008, 283, 9633–9641. [Google Scholar] [CrossRef] [PubMed]
  56. Sael, L.; Chitale, M.; Kihara, D. Structure- and Sequence-Based Function Prediction for Non-Homologous Proteins. J. Struct. Funct. Genom. 2012, 13, 111–123. [Google Scholar] [CrossRef] [PubMed]
  57. Meharenna, Y.T.; Oertel, P.; Bhaskar, B.; Poulos, T.L. Engineering Ascorbate Peroxidase Activity into Cytochrome c Peroxidase. Biochemistry 2008, 47, 10324–10332. [Google Scholar] [CrossRef] [PubMed]
  58. Islam, S.; Rahman, I.A.; Islam, T.; Ghosh, A. Genome-Wide Identification and Expression Analysis of Glutathione S-Transferase Gene Family in Tomato: Gaining an Insight to Their Physiological and Stress-Specific Roles. PLoS ONE 2017, 12, e0187504. [Google Scholar] [CrossRef]
  59. Naseri, G.; Koffas, M.A.G. Application of Combinatorial Optimization Strategies in Synthetic Biology. Nat. Commun. 2020, 11, 2446. [Google Scholar] [CrossRef]
  60. Li, N.; Qin, L.; Jin, M.; Zhang, L.; Geng, W.; Xiao, X. Extracellular Adsorption, Intracellular Accumulation and Tolerance Mechanisms of Cyclotella Sp. to Cr(VI) Stress. Chemosphere 2021, 270, 128662. [Google Scholar] [CrossRef]
Figure 1. Microalgal heavy metal bio-removing mechanisms and the use of synthetic biology for system enhancement. (A) Mechanisms and cellular functions used by microalgae to bio-remove heavy metals (HMs) from wastewater, adapted from Sattayawat et al., 2021 [7]; (i) biosorption, (ii) bioaccumulation (iii) biotransformation and (iv) cellular uptake of heavy metals via transporters (v) cellular responses to heavy metal-induced stresses. (B) Synthetic biology Design-Build-Test-Learn (DBTL) cycle. (C) Genetic part library construction to facilitate synthetic biology-based approaches for heavy metal bio-removal enhancement. The figure was created using https://biorender.com (accessed on 5 July 2022).
Figure 1. Microalgal heavy metal bio-removing mechanisms and the use of synthetic biology for system enhancement. (A) Mechanisms and cellular functions used by microalgae to bio-remove heavy metals (HMs) from wastewater, adapted from Sattayawat et al., 2021 [7]; (i) biosorption, (ii) bioaccumulation (iii) biotransformation and (iv) cellular uptake of heavy metals via transporters (v) cellular responses to heavy metal-induced stresses. (B) Synthetic biology Design-Build-Test-Learn (DBTL) cycle. (C) Genetic part library construction to facilitate synthetic biology-based approaches for heavy metal bio-removal enhancement. The figure was created using https://biorender.com (accessed on 5 July 2022).
Biology 11 01226 g001
Figure 2. A schematic computational workflow to construct a genetic part library for heavy metal bio-removal. Phase I focuses on protein sequence retrieval and sequence similarity search for protein homologs against available microalgal genomes, whereas Phase II focuses on functional annotation and protein structure modeling to validate their functions.
Figure 2. A schematic computational workflow to construct a genetic part library for heavy metal bio-removal. Phase I focuses on protein sequence retrieval and sequence similarity search for protein homologs against available microalgal genomes, whereas Phase II focuses on functional annotation and protein structure modeling to validate their functions.
Biology 11 01226 g002
Figure 3. Search results for each group of proteins presented separately by their fundamental roles in heavy metal bio-removing mechanisms. The green bars represent a total number of hits when each keyword was used to search against UniprotKB database, and the orange bars represent the number of hits after manual selection of the proteins.
Figure 3. Search results for each group of proteins presented separately by their fundamental roles in heavy metal bio-removing mechanisms. The green bars represent a total number of hits when each keyword was used to search against UniprotKB database, and the orange bars represent the number of hits after manual selection of the proteins.
Biology 11 01226 g003
Figure 4. Protein structures predicted using Alphafold2. (A) Glutamate-cysteine ligase (accession no. KAG7671258.1) and (B) Copper-transporting ATPase (accession no. XP_005845243.1). The structures are visualized by UCSF ChimeraX [48].
Figure 4. Protein structures predicted using Alphafold2. (A) Glutamate-cysteine ligase (accession no. KAG7671258.1) and (B) Copper-transporting ATPase (accession no. XP_005845243.1). The structures are visualized by UCSF ChimeraX [48].
Biology 11 01226 g004
Table 1. Microalgae with reported HM bio-removal capability.
Table 1. Microalgae with reported HM bio-removal capability.
MicroalgaReported MechanismReference
Chlorella
(taxid 3071)
Biosorption, bioaccumulation, biotransformation[27]
Scenedesmus
(taxid 3087)
Bioaccumulation, biotransformation, cellular stress response[28]
Table 2. Number of hypothetical proteins obtained from BLASTp analysis.
Table 2. Number of hypothetical proteins obtained from BLASTp analysis.
ProteinMicroalgal GenomeHypothetical ProteinHypothetical Protein after
Primary Function Prediction
Bioaccumulation
γ-Glutamylcysteine synthetase or Glutamate—cysteine ligaseChlorella42
Scenedesmus10
Phytochelatin synthaseChlorella22
Scenedesmus00
Glutathione synthetaseChlorella22
Scenedesmus00
Inorganic phosphate transporterChlorella107
Scenedesmus00
Biotransformation
Mercuric reductaseChlorella100
Scenedesmus50
Arsenate reductaseChlorella60
Scenedesmus10
Chromate reductaseChlorella00
Scenedesmus00
Metal transportation
Copper-transporting ATPaseChlorella269
Scenedesmus82
Manganese-transporting ATPaseChlorella90
Scenedesmus10
Cadmium-transporting ATPaseChlorella220
Scenedesmus40
Zinc-transporting ATPaseChlorella250
Scenedesmus60
Lead-transporting ATPaseChlorella110
Scenedesmus10
Manganese transporter domain (MntA) Chlorella431
Scenedesmus101
Manganese transporter domain (MntB)Chlorella271
Scenedesmus61
Manganese transporter domain
(MntC)
Chlorella00
Scenedesmus00
Manganese transporter domain
(MntD)
Chlorella
Scenedesmus
0
0
0
0
Mercuric transport protein (cytoplasmic/transmembrane)Scenedesmus00
Scenedesmus00
Mercuric transport protein
(periplasmic)
Chlorella00
Scenedesmus00
Zinc-regulated transporter (ZRT)Chlorella51
Scenedesmus10
Zinc transporter (ZIP)Chlorella131
Scenedesmus41
Natural resistance-associated macrophage proteinChlorella22
Scenedesmus00
Copper transporterChlorella00
Scenedesmus00
Cellular tolerance
Ascorbate peroxidaseChlorella98
Scenedesmus00
Superoxide dismutaseChlorella109
Scenedesmus44
Heme oxygenaseChlorella54
Scenedesmus22
Glutathione peroxidaseChlorella1510
Scenedesmus00
Glutathione S-transferaseChlorella4020
Scenedesmus112
Table 3. List of selected hypothetical protein matches when the template sequences were used to BLASTp against the selected microalgal genomes.
Table 3. List of selected hypothetical protein matches when the template sequences were used to BLASTp against the selected microalgal genomes.
GenomeAccession No.Signal Peptide *Putative Function
Bioaccumulation
Chlorella desiccata (nom. nud.)KAG7671258.1NGlutamate-cysteine ligase
Chlorella variabilisXP_005844806.1NGlutamate-cysteine ligase
Chlorella desiccata (nom. nud.)KAG7668718.1NPhytochelatin synthase
Chlorella variabilisXP_005845668.1NPhytochelatin synthase
Chlorella desiccata (nom. nud.)KAG7673317.1NGlutathione synthetase
Chlorella variabilisXP_005847003.1NGlutathione synthetase
Metal transportation
Chlorella variabilisXP_005845243.1YHeavy metal transporting ATPase
Chlorella variabilisXP_005851032.1YHeavy metal transporting ATPase
Scenedesmus sp. NREL 46B-D10KAF6264708.1YManganese transporter domain (MntA) **
Chlorella desiccata (nom. nud.)KAG7670010.1YManganese transporter domain (MntB)
Chlorella variabilisXP_005845281.1YManganese transporter domain (MntB) **
Chlorella variabilisXP_005844148.1YZinc-regulated transporter (ZRT)/Zinc transporter (ZIP) **
Chlorella variabilisXP_005846850.1YZinc-regulated transporter (ZRT)/Zinc transporter (ZIP) **
Chlorella desiccata (nom. nud.)KAG7667456.1YZinc-regulated transporter (ZRT)/Zinc transporter (ZIP)
Chlorella desiccata (nom. nud.)KAG7675010.1NNatural resistance-associated macrophage protein
Chlorella variabilisXP_005847346.1NNatural resistance-associated macrophage protein
Cellular tolerance
Chlorella variabilisXP_005842918.1NAscorbate peroxidase
Chlorella variabilisXP_005847371.1NAscorbate peroxidase
Chlorella desiccata (nom. nud.)KAG7671272.1NAscorbate peroxidase **
Chlorella desiccata (nom. nud.)KAG7672626.1NAscorbate peroxidase
Chlorella variabilisXP_005842951.1NAscorbate peroxidase
Chlorella desiccata (nom. nud.)KAG7671850.1NAscorbate peroxidase
Chlorella variabilisXP_005851196.1NAscorbate peroxidase
Chlorella desiccata (nom. nud.)KAG7671273.1NAscorbate peroxidase
Chlorella variabilisXP_005852313.1NSuperoxide dismutase
Chlorella variabilisXP_005852314.1NSuperoxide dismutase
Chlorella variabilisXP_005850331.1NSuperoxide dismutase
Chlorella variabilisXP_005850533.1NSuperoxide dismutase
Chlorella variabilisXP_005850825.1NSuperoxide dismutase
Chlorella variabilisXP_005851580.1NSuperoxide dismutase
Chlorella desiccata (nom. nud.)KAG7672127.1NSuperoxide dismutase
Chlorella desiccata (nom. nud.)KAG7672915.1NSuperoxide dismutase
Chlorella desiccata (nom. nud.)KAG7673432.1NSuperoxide dismutase
Scenedesmus sp. NREL 46B-D3KAF6253844.1NSuperoxide dismutase
Scenedesmus sp. PABB004KAF8054759.1NSuperoxide dismutase
Scenedesmus sp. PABB004KAF8070899.1NSuperoxide dismutase
Scenedesmus sp. PABB004KAF8072345.1NSuperoxide dismutase
Chlorella variabilisXP_005851913.1NHeme oxygenase 1
Chlorella desiccata (nom. nud.)KAG7671693.1NHeme oxygenase 1
Chlorella variabilisXP_005845884.1NHeme oxygenase 1
Chlorella variabilisXP_005842792.1NHeme oxygenase 1
Scenedesmus sp. NREL 46B-D3KAF6256065.1NHeme oxygenase 1
Scenedesmus sp. PABB004KAF8061310.1NHeme oxygenase 1
Chlorella variabilisXP_005852198.1NGlutathione peroxidase
Chlorella desiccata (nom. nud.)KAG7666639.1NGlutathione peroxidase
Chlorella variabilisXP_005847444.1NGlutathione peroxidase
Chlorella variabilisXP_005848232.1NGlutathione peroxidase
Chlorella variabilisXP_005851691.1NGlutathione peroxidase
Chlorella desiccata (nom. nud.)KAG7666823.1NGlutathione peroxidase
Chlorella variabilisXP_005850288.1NGlutathione peroxidase
Chlorella desiccata (nom. nud.)KAG7675006.1NGlutathione peroxidase
Chlorella variabilisXP_005844151.1NGlutathione peroxidase
Chlorella desiccata (nom. nud.)KAG7671063.1NGlutathione peroxidase
Chlorella desiccata (nom. nud.)KAG7667083.1NGlutathione S-transferase **
Chlorella desiccata (nom. nud.)KAG7667402.1NGlutathione S-transferase
Chlorella desiccata (nom. nud.)KAG7667544.1NGlutathione S-transferase
Chlorella desiccata (nom. nud.)KAG7667774.1NGlutathione S-transferase
Chlorella desiccata (nom. nud.)KAG7667817.1NGlutathione S-transferase
Chlorella desiccata (nom. nud.)KAG7669352.1NGlutathione S-transferase
Chlorella desiccata (nom. nud.)KAG7669598.1NGlutathione S-transferase **
Chlorella desiccata (nom. nud.)KAG7670514.1NGlutathione S-transferase
Chlorella desiccata (nom. nud.)KAG7675170.1NGlutathione S-transferase
Chlorella variabilisXP_005843180.1NGlutathione S-transferase
Chlorella variabilisXP_005845006.1NGlutathione S-transferase
Chlorella variabilisXP_005845127.1NGlutathione S-transferase
Chlorella variabilisXP_005845396.1NGlutathione S-transferase
Chlorella variabilisXP_005847002.1NGlutathione S-transferase
Chlorella variabilisXP_005848700.1NGlutathione S-transferase
Chlorella variabilisXP_005849485.1NGlutathione S-transferase
Chlorella variabilisXP_005849684.1NGlutathione S-transferase
Chlorella variabilisXP_005850654.1NGlutathione S-transferase
Chlorella variabilisXP_005852104.1NGlutathione S-transferase
Scenedesmus sp. NREL 46B-D3KAF6265595.1NGlutathione S-transferase
* Note that Y indicates Yes and N indicates No. The cutoff value is 0.1 for Signal Peptide (Sec/SPI). ** Note that these proteins were annotated differently when using ProFunc structure-based annotation.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Uttarotai, T.; Mukjang, N.; Chaisoung, N.; Pathom-Aree, W.; Pekkoh, J.; Pumas, C.; Sattayawat, P. Putative Protein Discovery from Microalgal Genomes as a Synthetic Biology Protein Library for Heavy Metal Bio-Removal. Biology 2022, 11, 1226. https://doi.org/10.3390/biology11081226

AMA Style

Uttarotai T, Mukjang N, Chaisoung N, Pathom-Aree W, Pekkoh J, Pumas C, Sattayawat P. Putative Protein Discovery from Microalgal Genomes as a Synthetic Biology Protein Library for Heavy Metal Bio-Removal. Biology. 2022; 11(8):1226. https://doi.org/10.3390/biology11081226

Chicago/Turabian Style

Uttarotai, Toungporn, Nilita Mukjang, Natcha Chaisoung, Wasu Pathom-Aree, Jeeraporn Pekkoh, Chayakorn Pumas, and Pachara Sattayawat. 2022. "Putative Protein Discovery from Microalgal Genomes as a Synthetic Biology Protein Library for Heavy Metal Bio-Removal" Biology 11, no. 8: 1226. https://doi.org/10.3390/biology11081226

APA Style

Uttarotai, T., Mukjang, N., Chaisoung, N., Pathom-Aree, W., Pekkoh, J., Pumas, C., & Sattayawat, P. (2022). Putative Protein Discovery from Microalgal Genomes as a Synthetic Biology Protein Library for Heavy Metal Bio-Removal. Biology, 11(8), 1226. https://doi.org/10.3390/biology11081226

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop